Search NASA⌕ Search

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.

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

At least 271 records · Page 15

Measuring rapid ocean tidal earth orientation variations with very long baseline interferometry

Ocean tidal effects on universal time and polar motion (UTPM) are investigated at four nearly diurnal (K(sub 1), P(sub 1), O(sub 1), and Q(sub 1)) and four nearly semidiurnal (K(sub 2), S(sub 2), M(sub 2), and N(sub 2)) frequencies by analyzing very long baseline interferometry (VLBI) data extending from 1978 to 1992. We discuss limitations of comparisons between experiment and theory for the retograde nearly diurnal polar motion components due to their degeneracy with prograde components of the nutation model. Estimating amplitudes of contributions to the modeled VLBI observables at these eight frequencies produces a statistically highly significant improvement of 7 mm to the residuals of a fit to the observed delays. Use of such an improved UTPM model also reduces the 14-30 mm scatter of baseline lengths about a time-linear model of tectonic motion by 3-14 mm, also withhigh significance levels. A total of 28 UTPM ocean tidal amplitudes can be unambiguously estimated from the data, with resulting UTI and PM magnitudes as large as 21 micro secs and 270 microarc seconds and formal uncertainties of the order of 0.3 micro secs and 5 microarc secs for UTI and PM, respectively. Empirically determined UTPM amplitudes and phases are com1pared to values calculated theoretically by Gross from Seiler's global ocean tide model. The discrepancy between theory and experiment is larger by a factor of 3 for UTI amplitudes (9 micro secs) than for prograde PM amplitudes (42 microarc secs). The 14-year VLBI data span strongly attenuates the influence of mismodeled effects on estimated UTPM amplitudes and phases that are not coherent with the eight frequencies of interest. Magnitudes of coherent and quasi-coherent systematic errors are quantified by means of internal consistency tests. We conclude that coherent systematic effects are many times larger than the formal uncertainties and can be as large as 4 micro secs for UTI and 60 microarc secs for polar motion. On the basis of such ealistic error estimates, 22 of the 31 fitted UTPM ocean tidal amplitudes differ from zero by more than 2 sigma.

Sovers, O. J.↗

Observations of Ocean Primary Productivity Using MODIS

Measuring the magnitude and variability of oceanic net primary productivity (NPP) represents a key advancement toward our understanding of the dynamics of marine ecosystems and the role of the ocean in the global carbon cycle. MODIS observations make two new contributions in addition to continuing the bio-optical time series begun with Orbview-2's SeaWiFS sensor. First, MODIS provides weekly estimates of global ocean net primary productivity on weekly and annual time periods, and annual empirical estimates of carbon export production. Second, MODIS provides additional insight into the spatial and temporal variations in photosynthetic efficiency through the direct measurements of solar-stimulated chlorophyll fluorescence. The two different weekly productivity indexes (first developed by Behrenfeld & Falkowski and by Yoder, Ryan and Howard) are used to derive daily productivity as a function of chlorophyll biomass, incident daily surface irradiance, temperature, euphotic depth, and mixed layer depth. Comparisons between these two estimates using both SeaWiFS and MODIS data show significant model differences in spatial distribution after allowance for the different integration depths. Both estimates are strongly dependence on the accuracy of the chlorophyll determination. In addition, an empirical approach is taken on annual scales to estimate global NPP and export production. Estimates of solar stimulated fluorescence efficiency from chlorophyll have been shown to be inversely related to photosynthetic efficiency by Abbott and co-workers. MODIS provides the first global estimates of oceanic chlorophyll fluorescence, providing an important proof of concept. MODIS observations are revealing spatial patterns of fluorescence efficiency which show expected variations with phytoplankton photo-physiological parameters as measured during in-situ surveys. This has opened the way for research into utilizing this information to improve our understanding of oceanic NPP variability. Deriving the ocean bio-optical properties places severe demands on instrument performance (especially band to band precision) and atmospheric correction. Improvements in MODIS instrument characterization and calibration over the first 16 mission months have greatly improved the accuracy of the chlorophyll input fields and FLH, and therefore the estimates of NPP and fluorescence efficiency. Annual estimates now show the oceanic NPP accounts for 40-50% of the global total NPP, with significant interannual variations related to large scale ocean processes. Spatial variations in ocean NPP, and exported production, have significant effects on exchange of CO2 between the ocean and atmosphere. Further work is underway to improve both the primary productivity model functions, and to refine our understanding of the relationships between fluorescence efficiency and NPP estimates. We expect that the MODIS instruments will prove extremely useful in assessing the time dependencies of oceanic carbon uptake and effects of iron enrichment, within the global carbon cycle.

Esaias, Wayne E.↗

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↗

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↗

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.↗

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.↗

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.↗

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↗

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↗

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↗

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.↗

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.↗

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↗

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.↗

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.↗

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.↗

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.↗

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↗