SEARCH · Search NASA

Results for “ocean model comparisons”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS2009

Combining Satellite Microwave Radiometer and Radar Observations to Estimate Atmospheric Latent Heating Profiles

In this study, satellite passive microwave sensor observations from the TRMM Microwave Imager (TMI) are utilized to make estimates of latent + eddy sensible heating rates (Q1-QR) in regions of precipitation. The TMI heating algorithm (TRAIN) is calibrated, or "trained" using relatively accurate estimates of heating based upon spaceborne Precipitation Radar (PR) observations collocated with the TMI observations over a one-month period. The heating estimation technique is based upon a previously described Bayesian methodology, but with improvements in supporting cloud-resolving model simulations, an adjustment of precipitation echo tops to compensate for model biases, and a separate scaling of convective and stratiform heating components that leads to an approximate balance between estimated vertically-integrated condensation and surface precipitation. Estimates of Q1-QR from TMI compare favorably with the PR training estimates and show only modest sensitivity to the cloud-resolving model simulations of heating used to construct the training data. Moreover, the net condensation in the corresponding annual mean satellite latent heating profile is within a few percent of the annual mean surface precipitation rate over the tropical and subtropical oceans where the algorithm is applied. Comparisons of Q1 produced by combining TMI Q1-QR with independently derived estimates of QR show reasonable agreement with rawinsonde-based analyses of Q1 from two field campaigns, although the satellite estimates exhibit heating profile structure with sharper and more intense heating peaks than the rawinsonde estimates. 2

Grecu, Mircea↗

NASA NTRS2019

Atmospheric Correction for Coastal Waters Based on Multi-Angle Polarimetric Observations

Atmospheric correction is an essential procedure of ocean color remote sensing, which obtains water-leaving signals from satellite or airborne sensors. Due to the small percentage of the water leaving signals in the total measurement, atmospheric correction requires precise evaluation of the radiometric contributions from the aerosol and ocean surface. It is often challenging over coastal waters when absorbing aerosols are present and when water leaving signals in the near infrared spectral region are non-negligible. In this work, we report on a joint retrieval algorithm that determines aerosol and ocean optical properties using the polarimetric measurements based on a coupled atmosphere and ocean radiative transfer model. Two bio-optical models are designed for ocean water properties. One model parameterizes the open ocean optical properties in terms of the concentration of chlorophyll a. The other is a generalized bio-optical model for coastal waters that describes the absorption and scattering by phytoplankton, colored dissolved organic matter and non-algal particles using multiple parameters. We applied the algorithm to the airborne Research Scanning Polarimeter (RSP) measurements acquired over both open and coastal ocean waters. Through the comparisons with in situ ocean color measurements, the flexibility and accuracy of the retrieval algorithm in retrieving aerosol and water leaving radiance properties are demonstrated under various aerosol and ocean water conditions.

Gao, Meng↗

NASA NTRS1976

Ocean properties

The author has identified the following significant results. Results of testing the CP program indicate that the best results can be obtained in the near infrared water bands. The absorption due to water vapor and carbon dioxide in the thermal infrared band appeared to be less reliable in comparison to spacecraft-acquired data and band models. Comparisons of laboratory carbon dioxide transmission in the thermal infrared band show good agreement except in regions where lines are known to be missing. The comparison of ozone transmission at a wavelength of 9.6 micrometers to laboratory data showed unexceptedly large differences.

Korb, C. L.↗

NASA NTRS1983

An accurate radiative transfer model for use in the direct physical inversion of HIRS2 and MSU temperature sounding data

The direct computation of atmospheric transmittance and clear column radiances for the channels of HIRS2 and MSU as a function of atmospheric and surface conditions is described in detail. A comparison is made between the observations and the calculated radiances derived from colocated oceanic radiosondes. It is found that under clear conditions, calculated brightness temperatures for the HIRS2 have a standard deviation of the order of 0.7 C compared with observations, whereas MSU channels have a standard deviation of approximately 1 C. In some channels, small biases are found that can be removed by an empirical 'tuning' with coefficients that can be successfully transferred from one season to another. Less satisfactory agreement is obtained from a comparison of calculations with 'reconstructed' clear radiances, which are used in analyzing sounding data under partially cloudy conditions.

Susskind, J.↗

DOE OSTI2020

Control of Vibration Systems with Mechanical Motion Rectifier and their Applications to Vehicle Suspension and Ocean Energy Harvester

Vibration control is a large branch in control research, because all moving systems may induce desired or undesired vibration. Due to the limitation of passive system’s adaptability and changing excitation input, vibration control brings the solution to change system dynamic with desired behavior to fulfill control targets. According to preference, vibration control can be separated into two categories: vibration reduction and vibration amplification. Lots of research papers only examine one aspect in vibration control. The thesis investigates the control development for both control targets with two different control applications: vehicle suspension and ocean wave energy converter. It develops control methods for both systems with simplified modeling setup, then followed by the application of a novel mechanical motion rectifier (MMR) gearbox that uses mechanical one-way clutches in both systems. The flow is from the control for common system to the control design for a specifically designed system. In the thesis, active (model predictive control: MPC), semi-active (Skyhook, skyhook-power driven damper: SHPDD, hybrid model predictive control: HMPC), and passive control (Latching Control) methods are developed for different applications or control performance comparison on single system. The thesis also studies about new type of system with switching mechanism, in which other papers do not talk too much and possible control research direction to deal with such complicated system in vibration control. The state-space modeling for both systems are provided in the thesis with detailed model of the MMR gearbox. From the simulation, it can be shown that in the vehicle suspension application, the controlled MMR gearbox can be effective in improving vehicle ride comfort by 29.2% compared to that of the traditional hydraulic suspension. In the ocean wave energy converter, the controlled MMR WEC with simple latching control can improve the power generation by 57% compared to the passive MMR WEC. Besides, the passive MMR WEC also shows its advantage on the passive direct drive WEC in power generation improvement. From the control development flow for the MMR system, the limitation of the MMR gearbox is also identified, which introduces the future work in developing active-MMR gearbox by using an electromagnetic clutch. Some possible control development directions on the active-MMR is also mentioned at the end of the thesis to provide reference for future works.

16 TIDAL AND WAVE POWER↗

NASA NTRS1988

Microscale ocean dynamics

The detailed dynamics of the micro-scale ocean surface phenomena are studied along with the relationships among the surface signatures with the underlying dynamical processes. The approach to advance the understanding in this area is as follows: (1) Conduct rigorous theoretical studies of the ocean surface wave dynamics and statistical properties; (2) Conduct process-oriented laboratory experiments to verify the theoretical results, and to provide guidance for further studies; and (3) Prepare testable hypotheses for field verifications and comparisons during the ONR/NASA sponsored Surface Wave Dynamics Experiment (SWADE). An analytic model was established for wave breaking probability to study the influence of wave breaking on the spectrum shape in both deep and finite-depth waters. The spectrum was processed along with the structures of the water surface under the influence of wind and existing waves.

Huang, N.↗

NASA NTRS2011

Spectral Kernel Approach to Study Radiative Response of Climate Variables and Interannual Variability of Reflected Solar Spectrum

The radiative kernel approach provides a simple way to separate the radiative response to different climate parameters and to decompose the feedback into radiative and climate response components. Using CERES/MODIS/Geostationary data, we calculated and analyzed the solar spectral reflectance kernels for various climate parameters on zonal, regional, and global spatial scales. The kernel linearity is tested. Errors in the kernel due to nonlinearity can vary strongly depending on climate parameter, wavelength, surface, and solar elevation; they are large in some absorption bands for some parameters but are negligible in most conditions. The spectral kernels are used to calculate the radiative responses to different climate parameter changes in different latitudes. The results show that the radiative response in high latitudes is sensitive to the coverage of snow and sea ice. The radiative response in low latitudes is contributed mainly by cloud property changes, especially cloud fraction and optical depth. The large cloud height effect is confined to absorption bands, while the cloud particle size effect is found mainly in the near infrared. The kernel approach, which is based on calculations using CERES retrievals, is then tested by direct comparison with spectral measurements from Scanning Imaging Absorption Spectrometer for Atmospheric Cartography (SCIAMACHY) (a different instrument on a different spacecraft). The monthly mean interannual variability of spectral reflectance based on the kernel technique is consistent with satellite observations over the ocean, but not over land, where both model and data have large uncertainty. RMS errors in kernel ]derived monthly global mean reflectance over the ocean compared to observations are about 0.001, and the sampling error is likely a major component.

Jin, Zhonghai↗

NASA NTRSDate not supplied

Distribution and Sources of Tropospheric Aerosols Over the Western North Atlantic During ACTIVATE (February-March 2020)

The Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE) is a five-year (2019-2023) NASA Earth-Venture Suborbital-3 (EVS-3) mission to robustly characterize aerosol-cloud-meteorology interactions during February-June, with a focus on marine boundary layer clouds. This characterization requires understanding of aerosol composition, distribution, transport pathways, and sources. Here we use the GEOS-Chem chemical transport model driven by the MERRA-2 reanalysis to simulate tropospheric carbon monoxide (CO) and aerosols over the western North Atlantic Ocean (WNAO) during the first field deployment of ACTIVATE (February-March 2020). CO is used as a tracer to diagnose transport pathways for continental pollution outflow to the WNAO. The aerosol simulation represents all major aerosols including sulfate-nitrate-ammonium (SNA), mineral dust, sea salt, black and organic carbon aerosols. Model-simulated mixing ratios of SNA and organic carbon are compared with aircraft observations. We show that intensive aerosol mass measurements in the lower troposphere over the North American outflow region during ACTIVATE provide strong constraints on model aerosol wet scavenging. Comparisons of model aerosol extinctions with the airborne High Spectral Resolution Lidar-2 (HSRL-2) measurements indicate that the model generally captures continental outflow of aerosols and enhanced extinctions associated with sea salt. Simulated aerosol optical depths (AODs) and extinction vertical profiles are also compared with satellite retrievals from MODIS/Aqua and CALIOP/CALIPSO, respectively. CALIOP sometimes observes enhanced aerosol extinction at altitudes between ~1.5km and ~2.0km south of 32°N. The model reproduces this enhancement and attributes it to coarse-mode sea salt associated with high relative humidity. The relative contributions of continental and oceanic sources to the aerosol loading, AOD, and their distributions over the study domain will also be discussed.

aerosols↗

NASA NTRS2004

Lunar and Planetary Science XXXV: Moon and Mercury

The session" Moon and Mercury" included the following reports:Helium Production of Prompt Neutrinos on the Moon; Vapor Deposition and Solar Wind Implantation on Lunar Soil-Grain Surfaces as Comparable Processes; A New Lunar Geologic Mapping Program; Physical Backgrounds to Measure Instantaneous Spin Components of Terrestrial Planets from Earth with Arcsecond Accuracy; Preliminary Findings of a Study of the Lunar Global Megaregolith; Maps Characterizing the Lunar Regolith Maturity; Probable Model of Anomalies in the Polar Regions of Mercury; Parameters of the Maximum of Positive Polarization of the Moon; Database Structure Development for Space Surveying Results by Moon -Zond Program; CM2-type Micrometeoritic Lunar Winds During the Late Heavy Bombardment; A Comparison of Textural and Chemical Features of Spinel Within Lunar Mare Basalts; The Reiner Gamma Formation as Characterized by Earth-based Photometry at Large Phase Angles; The Significance of the Geometries of Linear Graben for the Widths of Shallow Dike Intrusions on the Moon; Lunar Prospector Data, Surface Roughness and IR Thermal Emission of the Moon; The Influence of a Magma Ocean on the Lunar Global Stress Field Due to Tidal Interaction Between the Earth and Moon; Variations of the Mercurian Photometric Relief; A Model of Positive Polarization of Regolith; Ground Truth and Lunar Global Thorium Map Calibration: Are We There Yet?;and Space Weathering of Apollo 16 Sample 62255: Lunar Rocks as Witness Plates for Deciphering Regolith Formation Processes.

Source record↗

NASA NTRS2012

Evaluation of Cloud Microphysics Simulated using a Meso-Scale Model Coupled with a Spectral Bin Microphysical Scheme through Comparison with Observation Data by Ship-Borne Doppler and Space-Borne W-Band Radars

Equivalent radar reflectivity factors (Ze) measured by W-band radars are directly compared with the corresponding values calculated from a three-dimensional non-hydrostatic meso-scale model coupled with a spectral-bin-microphysical (SBM) scheme for cloud. Three case studies are the objects of this research: one targets a part of ship-borne observation using 95 GHz Doppler radar over the Pacific Ocean near Japan in May 2001; other two are aimed at two short segments of space-borne observation by the cloud profiling radar on CloudSat in November 2006. The numerical weather prediction (NWP) simulations reproduce general features of vertical structures of Ze and Doppler velocity. A main problem in the reproducibility is an overestimation of Ze in ice cloud layers. A frequency analysis shows a strong correlation between ice water contents (IWC) and Ze in the simulation; this characteristic is similar to those shown in prior on-site studies. From comparing with the empirical correlations by the prior studies, the simulated Ze is overestimated than the corresponding values in the studies at the same IWC. Whereas the comparison of Doppler velocities suggests that large-size snowflakes are necessary for producing large velocities under the freezing level and hence rules out the possibility that an overestimation of snow size causes the overestimation of Ze. Based on the results of several sensitivity tests, we conclude that the source of the overestimation is a bias in the microphysical calculation of Ze or an overestimation of IWC. To identify the source of the problems needs further validation research with other follow-up observations.

Iguchi, T.↗

DOE OSTI2021

Benchmark estimate of the effect of anthropogenic emissions on the ocean surface

Abstract Investigations into the role of anthropogenic emissions in the occurrence of extreme weather often use a method that compares simulations of atmospheric climate models run under a factual scenario of historical boundary conditions observed during the period of the event against simulations run under a counterfactual scenario of what those boundary conditions might naturally have been over that same period in the absence of anthropogenic emissions. A particular requirement for this experiment design is an accurate estimation of ocean surface boundary conditions for use by the counterfactual natural simulations. Here we use output from the CMIP5 multi‐climate‐model archive to develop a robust estimate of sea surface temperatures and sea ice conditions for use in counterfactual natural simulations, intended as a benchmark estimate to facilitate comparison across climate models and across studies. This development includes tests to ensure that the final estimate is stable from year‐to‐year and stable against other perturbations to the methodology, as well as consideration of the strengths and weaknesses in comparison to other available attributable warming estimates. While this estimate is tailored specifically for the International CLIVAR C20C+ Detection and Attribution Project, it can be used by related projects as well.

Stone, Dáithí A.↗

NASA NTRS1982

The Seasat-A satellite scatterometer - The geophysical evaluation of remotely sensed wind vectors over the ocean

A description is given of the algorithm used to convert Seasat-A satellite microwave scatterometer measurements of ocean normalized radar cross section to the neutral stability vector wind at 19.5 m height, as well as to compare these winds with high-quality surface observations. The wind vector algorithm used an empirical normalized radar cross section model function to describe the ocean normalized radar cross section's dependence on the 19.5-m neutral stability wind vector. In addition, two model functions were evaluated by means of an independent set of in situ surface wind observations from the Joint Air Sea Interaction experiment (JASIN). Better results were produced by these comparisons than the stipulated Seasat wind speed and direction accuracy specifications of + or - 2 m/sec and + or - 20 deg, respectively, over the 0-16 m/sec range of winds observed during JASIN.

Jones, W. L.↗

NASA NTRS1992

Analysis of a general circulation model. II - Distribution of kinetic energy in the South Atlantic and Kuroshio/Oyashio systems

It was found (Garzoli et al., 1992) that the general circulation model of Semtner and Chervin (1992) provides accurate descriptions of the Brazil-Malvinas and the Kuroshio/Oyashio confluence systems, except for the fact that the model prediction shows less variability than that present in observations. This paper investigates the problem of model variability by analyzing the mean and the eddy kinetic energy from the model and comparing the values with the Geosat altimeter observations for the South Atlantic Ocean and for the Kuroshio system. It is found that, while the model shows transient eddy activity in the areas that overlap the Geosat observations, the energy level of the model transient motions is considerably smaller following an arch along the bottom topography. The same was found from the comparisons made with values obtained from FGGE and surface drifters. It is suggested that the model is poorly resolving instabilities in the confluence front, and is not resolving other transients appearing in regions of marked topography.

Garraffo, Zulema↗

NASA NTRS1998

The Development of a Degree 360 Expansion of the Dynamic Ocean Topography of the POCM_4B Global Circulation Model

This paper documents the development of a degree 360 expansion of the dynamic ocean topography (DOT) of the POCM_4B ocean circulation model. The principles and software used that led to the final model are described. A key principle was the development of interpolated DOT values into land areas to avoid discontinuities at or near the land/ocean interface. The power spectrum of the POCM_4B is also presented with comparisons made between orthonormal (ON) and spherical harmonic magnitudes to degree 24. A merged file of ON and SH computed degree variances is proposed for applications where the DOT power spectrum from low to high (360) degrees is needed.

Rapp, Richard H.↗

NASA NTRS2007

Uncertainty Analysis of Historical Hurricane Data

An analysis of variance (ANOVA) study was conducted for historical hurricane data dating back to 1851 that was obtained from the U. S. Department of Commerce National Oceanic and Atmospheric Administration (NOAA). The data set was chosen because it is a large, publicly available collection of information, exhibiting great variability which has made the forecasting of future states, from current and previous states, difficult. The availability of substantial, high-fidelity validation data, however, made for an excellent uncertainty assessment study. Several factors (independent variables) were identified from the data set, which could potentially influence the track and intensity of the storms. The values of these factors, along with the values of responses of interest (dependent variables) were extracted from the data base, and provided to a commercial software package for processing via the ANOVA technique. The primary goal of the study was to document the ANOVA modeling uncertainty and predictive errors in making predictions about hurricane location and intensity 24 to 120 hours beyond known conditions, as reported by the data set. A secondary goal was to expose the ANOVA technique to a broader community within NASA. The independent factors considered to have an influence on the hurricane track included the current and starting longitudes and latitudes (measured in degrees), and current and starting maximum sustained wind speeds (measured in knots), and the storm starting date, its current duration from its first appearance, and the current year fraction of each reading, all measured in years. The year fraction and starting date were included in order to attempt to account for long duration cyclic behaviors, such as seasonal weather patterns, and years in which the sea or atmosphere were unusually warm or cold. The effect of short duration weather patterns and ocean conditions could not be examined with the current data set. The responses analyzed were the storm latitude, longitude and intensity, as recorded in the data set, 24 or 120 hours beyond the current state. Several ANOVA modeling schemes were examined. Two forms of validation were used: 1) comparison with official hurricane prediction performance metrics and 2) cases studies conducted on hurricanes from the 2005 season, which were not included within the model construction and ANOVA assessment. In general, the ANOVA technique did not perform as well as the established official prediction performance metrics published by NOAA; still, the technique did remarkably well in this demonstration with a difficult data set and could probably be made to perform better with more knowledge of hurricane development and dynamics applied to the problem. The technique provides a repeatable prediction process that eliminates the need for judgment in the forecast.

Green, Lawrence L.↗

NASA NTRS1998

Wind Forcing of the North Sea Pole Tide

The Chandler wobble of the earth's rotation has a period near 14 months and sets up the 0.5 cm amplitude pole tide in the deep oceans, However, the pole tide is anomalously large in the North Sea, where the amplitude increases sharply up to 3 cm along the continental coast. It is shown here that the sea levels are well correlated with winds at the pole tide frequency. The Princeton Ocean Model is used to investigate the response of the North Sea to wind forcing. The barotropic numerical ocean model depicts realistic coastlines and bathymetry at 5 ft x 5 ft resolution, with 97 x 73 grid points. The monthly mean wind fields for the 40-year period (1958-1997) from the National Centers for Atmospheric Prediction (NCEP) reanalysis were used to force the model. The winds were converted to stress with a neutral drag coefficient that varied linearly with windspeed (instead of using the NCEP windstress). A 5-day simulation was made for each month until the resulting flow regime came into equilibrium, and model water levels at various station locations were saved for comparison with tidal-gauge observed sea levels from the Permanent Service for Mean Sea Level (PSMSL). The comparison is made for 10 North Sea stations with high quality tide gauge data. We find: (1) good agreement in annual and semi-annual phases and in the trend of amplitude w.r.t. latitude; (2) more importantly for this study, the model-predicted and observed non-seasonal sea level variations show a very significant temporal correlation as well as spectral coherence. However, a large amplitude difference exists between the two sets -- the overall amplitude variability of the observed is generally a factor of 2-3 larger than the model prediction (this same phenomenon has been reported in ocean circulation studies, although the cause is not yet clear.) Our results indicate that the wind forcing is the main cause of the observed large pole tide in the North Sea.

OConnor, W.↗

NASA NTRS2001

Measurements of Aerosol Vertical Profiles and Optical Properties during INDOEX 1999 Using Micro-Pulse Lidars

Micro-pulse lidar systems (MPL) were used to measure aerosol properties during the Indian Ocean Experiment (INDOEX) 1999 field phase. Measurements were made from two platforms: the NOAA ship RN Ronald H. Brown, and the Kaashidhoo Climate Observatory (KCO) in the Maldives. Sunphotometers were used to provide aerosol optical depths (AOD) needed to calibrate the MPL. This study focuses on the height distribution and optical properties (at 523 nm) of aerosols observed during the campaign. The height of the highest aerosols (top height) was calculated and found to be below 4 km for most of the cruise. The marine boundary layer (MBL) top was calculated and found to be less than 1 km. MPL results were combined with air mass trajectories, radiosonde profiles of temperature and humidity, and aerosol concentration and optical measurements. Humidity varied from approximately 80% near the surface to 50% near the top height during the entire cruise. The average value and standard deviation of aerosol optical parameters were determined for characteristic air mass regimes. Marine aerosols in the absence of any continental influence were found to have an AOD of 0.05 +/- 0.03, an extinction-to-backscatter ratio (S-ratio) of 33 +/- 6 sr, and peak extinction values around 0.05/km (near the MBL top). The marine results are shown to be in agreement with previously measured and expected values. Polluted marine areas over the Indian Ocean, influenced by continental aerosols, had AOD values in excess of 0.2, S-ratios well above 40 sr, and peak extinction values approximately 0.20/km (near the MBL top). The polluted marine results are shown to be similar to previously published values for continental aerosols. Comparisons between MPL derived extinction near the ship (75 m) and extinction calculated at ship-level using scattering measured by a nephelometer and absorption using a PSAP were conducted. The comparisons indicated that the MPL algorithm (using a constant S-ratio throughout the lower troposphere) calculates extinction near the surface in agreement with the ship-level measurements only when the MBL aerosols are well mixed with aerosols above. Finally, a review of the MPL extinction profiles showed that the model of aerosol vertical extinction developed during an earlier INDOEX field campaign (at the Maldives) did not correctly describe the true vertical distribution over the greater Indian Ocean region. Using the average extinction profile and AOD obtained during marine conditions, a new model of aerosol vertical extinction was determined for marine atmospheres over the Indian Ocean. A new model of aerosol vertical extinction for polluted marine atmospheres was also developed using the average extinction profile and AOD obtained during marine conditions influenced by continental aerosols.

Welton, Ellsworth J.↗

NASA NTRS2021

Evaluation of Aqua MODIS Thermal Emissive Bands Stability through Radiative Transfer Modeling

Moderate Resolution Imaging Spectroradiometer (MODIS) on Aqua has been inoperation providing continuous global observations for science research and applications since2002. The long-term stability of thermal emissive bands (TEBs) of Aqua MODIS was monitoredthrough inter-comparisons with measurements by hyperspectral or multi-spectral infraredsensors or through vicarious monitoring over cold targets such as Dome-C and deep convectiveclouds. The radiative transfer modeling (RTM)-based simulation models are developed to per-form the long-term monitoring of the stability of Aqua MODIS TEBs. Long-term EuropeanCentre for Medium-Range Weather Forecasts global atmospheric reanalysis data are used asinputs to the RTM simulation. By confining ocean in low to middle latitudes as the area ofinterest, the long-term stabilities of Aqua MODIS TEBs are monitored through observationminus background (O-B) brightness temperature (BT) bias between MODIS measurements andsimulations from the RTM. In general, the RTM-based O-B analysis shows that the long-termstability of Aqua MODIS TEBs are maintained well. It is shown that Aqua MODIS surfacebands are all radiometrically stable with yearly BT bias drifts <0.005 K∕yearfor B20, B22,and B23 and∼0.01 K∕yearfor B31 and B32. The carbon dioxide (CO2) absorption channelsof Aqua MODIS, e.g., B33 to B36, are stable with a BT bias drift <0.005 K∕year. It is also foundthat after accounting for the long-term growth of greenhouse gasN2O, the O-B bias trends forB24 and B25 are quite stable with yearly BT bias drifts around 0.0002 and0.0055 K∕year,respectively. The stabilities of the moisture-sensitive channels B27 to B29 and the ozone channelB30 of Aqua MODIS are also evaluated and the remnant small variations in the O-B bias trend-ing of these channels are further discussed. The RTM-based inter-comparison of MODIS TEBmeasurements with global atmospheric and climate re-analysis data provides validation of thelong-term radiometric stability of Aqua MODIS TEBs. The analysis in this paper also demon-strates that the comparison of well-calibrated MODIS TEB measurements with RTM simulationhelps quantify the long-term impacts of global greenhouse gas concentration growth on the BTdecreases in the MODIS greenhouse gas-sensitive channels.

MODIS↗