Search NASA⌕ Search

SEARCH · Search NASA

Results for “global uncertainty analysis”

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 325 records · Page 18

Tuning and Robustness Analysis for the Orion Absolute Navigation System

The Orion Multi-Purpose Crew Vehicle (MPCV) is currently under development as NASA's next-generation spacecraft for exploration missions beyond Low Earth Orbit. The MPCV is set to perform an orbital test flight, termed Exploration Flight Test 1 (EFT-1), some time in late 2014. The navigation system for the Orion spacecraft is being designed in a Multi-Organizational Design Environment (MODE) team including contractor and NASA personnel. The system uses an Extended Kalman Filter to process measurements and determine the state. The design of the navigation system has undergone several iterations and modifications since its inception, and continues as a work-in-progress. This paper seeks to show the efforts made to-date in tuning the filter for the EFT-1 mission and instilling appropriate robustness into the system to meet the requirements of manned space ight. Filter performance is affected by many factors: data rates, sensor measurement errors, tuning, and others. This paper focuses mainly on the error characterization and tuning portion. Traditional efforts at tuning a navigation filter have centered around the observation/measurement noise and Gaussian process noise of the Extended Kalman Filter. While the Orion MODE team must certainly address those factors, the team is also looking at residual edit thresholds and measurement underweighting as tuning tools. Tuning analysis is presented with open loop Monte-Carlo simulation results showing statistical errors bounded by the 3-sigma filter uncertainty covariance. The Orion filter design uses 24 Exponentially Correlated Random Variable (ECRV) parameters to estimate the accel/gyro misalignment and nonorthogonality. By design, the time constant and noise terms of these ECRV parameters were set to manufacturer specifications and not used as tuning parameters. They are included in the filter as a more analytically correct method of modeling uncertainties than ad-hoc tuning of the process noise. Tuning is explored for the powered-flight ascent phase, where measurements are scarce and unmodelled vehicle accelerations dominate. On orbit, there are important trade-off cases between process and measurement noise. On entry, there are considerations about trading performance accuracy for robustness. Process Noise is divided into powered flight and coasting ight and can be adjusted for each phase and mode of the Orion EFT-1 mission. Measurement noise is used for the integrated velocity measurements during pad alignment. It is also used for Global Positioning System (GPS) pseudorange and delta- range measurements during the rest of the flight. The robustness effort has been focused on maintaining filter convergence and performance in the presence of unmodeled error sources. These include unmodeled forces on the vehicle and uncorrected errors on the sensor measurements. Orion uses a single-frequency, non-keyed GPS receiver, so the effects due to signal distortion in Earth's ionosphere and troposphere are present in the raw measurements. Results are presented showing the efforts to compensate for these errors as well as characterize the residual effect for measurement noise tuning. Another robustness tool in use is tuning the residual edit thresholds. The trade-off between noise tuning and edit thresholds is explored in the context of robustness to errors in dynamics models and sensor measurements. Measurement underweighting is also presented as a method of additional robustness when processing highly accurate measurements in the presence of large filter uncertainties.

Holt, Greg N.↗

South Asian Monsoon Response to Regional Aerosol Emission Reductions: Insights From RAMIP

Anthropogenic aerosol and associated precursor gas emissions are already declining in many regions and are likely to decline significantly by 2050, with major implications for regional climate. Unlike greenhouse gases, aerosol impacts are spatially heterogeneous and can influence climate both near emission sources and through remote teleconnections. This is particularly important for the South Asian monsoon system, where both local and remote aerosol changes can significantly affect precipitation patterns. Using simulations from the Regional Aerosol Model Intercomparison Project, we examine how local and remote aerosol emission reductions influence South Asian climate across both pre-monsoon and monsoon seasons, including weakening of the elevated heat pump (EHP) mechanism during pre-monsoon months. Our analysis employs 10-member ensembles from 10 CMIP6-era models to compare three experiments with global, South Asian, and East Asian aerosol reductions relative to a high-emission baseline (SSP3-7.0). This experimental design allows us to isolate and quantify the distinct impacts of regional emission changes. Results reveal that global aerosol reductions produce a larger all-India precipitation increase (+0.28 mm day −1 ) than South Asian reductions alone (+0.19 mm day −1 ), with the strongest regional responses over the northern Bay of Bengal, the Western Ghats and Indo-Gangetic Plains. East Asian reductions show uncertain but modest precipitation reductions of 0.2-0.6 mm day −1 over parts of west-central and eastern Indian regions. Pre-monsoon carbonaceous aerosol reductions systematically weaken the EHP, cooling the mid-troposphere over the Himalayan foothills by up to 0.4 K and producing localised anomalous descent that opposes the climatological pre-monsoon ascent. Substantial inter-model diversity exists in the spatial patterns and magnitudes of these responses, highlighting key uncertainties in aerosol–monsoon interactions. The South Asian precipitation response is driven by enhanced land–sea thermal contrast, while the pre-monsoon EHP weakening is attributable exclusively to local South Asian carbonaceous aerosol forcing rather than remote teleconnections. These findings have direct implications for air quality and climate adaptation planning across South Asia as regional aerosol emissions diverge under different development pathways.

cmip6↗

A View of the Lunar Interior Through Lunar Laser Range Analysis

Laser ranges between observatories on the Earth and retroreflectors on the Moon started in 1969 and continue to the present. Recent range accuracies are 2 cm while earliest ranges are an C, order of magnitude less certain. Four retroreflectors are ranged: three located at the Apollo 11, 14, and 15 sites and one on the Lunakhod 2 rover. Accurate analysis of the range data determines a number of lunar science parameters. The lunar interior variables include a fluid core parameter. The Lunar Laser Ranging effort is reviewed elsewhere. Many parameters are detected through their influence on rotation. Also detected are solid-body tides and accurate selenocentric reflector locations. Determined through the rotation are moment-of-inertia differences, gravitational harmonics, potential Love number, and dissipation effects due to tides and molten core. The rotation of the Moon is not at its minimum energy state; some recently active process has caused free librations. The moment differences contributed to the recent improvement of the Moon's moment of inertia from the Lunar Prospector gravity field. The Love numbers provide bulk elastic properties. Future possibilities for measurement include oblateness of the core-mantle boundary and core moment. A study of dissipation signatures in the rotation determines tidal Q vs. frequency and concludes that the Moon has a molten core. At 1 month the tidal Q is 37 and at 1 yr it is 60. The core radius is < or = 352 km for Fe and < or = 374 km for the Fe-FeS eutectic. The core detection exceeds 3x its uncertainty. The spin of the core is not aligned with the spin of the mantle and torque arises from the velocity difference at the boundary. Yoder's turbulent boundary layer theory is used to compute the radii. The present heat generation from tides and core interaction is minor compared to radiogenic heating. The heating for ancient times is more interesting. Peale and Cassen investigated lunar tidal heating while the lunar orbit expanded due to tides on Earth. Their calculations predate the measurement of Q and should be multiplied by 3.45 to match the lunar-laser-determined Love number and monthly Q. Tidal-heating computations depend on how fast the lunar orbit evolved and whether the tidal dissipation is localized. Neither is known, but plausible assumptions lead to early central region temperature increases of several hundred degrees. Most of the energy is deposited early in the Moon's history. The turbulent boundary layer theory allows a prediction of energy dissipated at the core-mantle boundary during orbit evolution. Under the assumption that the properties of the early core are the same as at present, the energy dissipated by core-mantle interaction is about the same as for tidal dissipation, but it is deposited in a smaller volume. This source of energy is capable of promoting convection in an early fluid core and driving a dynamo. This is a transient phase with a duration depending on the rate of orbit evolution. Plausible assumptions lead to a duration of a few hundred million years. Thus, the remnant magnetization of many lunar rocks is compatible with a brief global magnetic field powered by dynamical energy dissipation. Analysis of the lunar laser ranges is providing information on lunar geophysics. Future data will improve accuracies of present solution parameters, and several more interior effects should be detectable.

Williams, J. G.↗

Towards an Optimal Estimation Retrieval of Cirrus Cloud Optical and Microphysical Properties Using Hyperspectral Shortwave Instruments and A Fast Radiative Transfer Algorithm

Cirrus cloud retrieval products (here, cloud optical depth, effective particle size, and cloud top height) are important inputs into numerical weather and climate models. Uncertainties in such retrieval products, as a matter of course, propagate downstream, impacting model calculations. Improvements in high-quality global cirrus cloud optical and microphysical data products from satellite observations are needed to understand and reduce retrieval uncertainties. Hyperspectral shortwave instruments produce high-resolution and information-dense spectra, thus offering the opportunity to reduce uncertainties in retrieval products. We are in the process of developing a retrieval that uses the very fast Principal Component Radiative Transfer Model in the solar spectral region (PCRTM-Solar) in the forward model calculations. This retrieval will use measured reflectances from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) and the forthcoming Climate Absolute Radiance and Refractivity Observatory Pathfinder (CLARREO-Pathfinder) instruments. In this manuscript we present progress towards a reference retrieval employing a widely used, verified, accurate, yet computationally slower radiative transfer modeling technique. The reference retrieval, while too slow for using the complete hyperspectral measurement, will allow us to study the behavior of retrieval products and help us verify results from our in-development fast retrieval. Both retrievals will use the optimal estimation retrieval framework. In this manuscript we present results from an uncertainty analysis considering three uncertainty sources for a cirrus cloud retrieval in the form of error covariance matrices: reflectance uncertainty due to water vapor, the reflectance uncertainty due to ice crystal scattering assumptions, and the instrument measurement uncertainty. Results show that the uncertainty due to habit selection is the largest, while that due to water vapor is at most 0.6% relative to channel reflectance. As a first step, the retrieval is being designed for single layer ice clouds over open ocean water.

Cirrus cloud↗

Daily and Hourly Variability in Global Fire Emissions and Consequences for Atmospheric Model Predictions of Carbon Monoxide

Attribution of the causes of atmospheric trace gas and aerosol variability often requires the use of high resolution time series of anthropogenic and natural emissions inventories. Here we developed an approach for representing synoptic- and diurnal-scale temporal variability in fire emissions for the Global Fire Emissions Database version 3 (GFED3). We distributed monthly GFED3 emissions during 2003-2009 on a daily time step using Moderate Resolution Imaging Spectroradiometer (MODIS)-derived measurements of active fires from Terra and Aqua satellites. In parallel, mean diurnal cycles were constructed from Geostationary Operational Environmental Satellite (GOES) active fire observations. We found that patterns of daily variability in fires varied considerably across different biomes, with short but intense periods of daily emissions in boreal ecosystems and lower intensity (but more continuous) periods of bunting in savannas. On diurnal timescales, our analysis of the GOES active fires indicated that fires in savannas, grasslands, and croplands occurred earlier in the day as compared to fires in nearby forests. Comparison with Total Carbon Column Observing Network (TCCON) and Measurements of Pollution in the Troposphere (MOPITT) column CO observations provided evidence that including daily variability in emissions moderately improved atmospheric model simulations, particularly during the fire season and near regions with high levels of biomass burning. The high temporal resolution estimates of fire emissions developed here may ultimately reduce uncertainties related to fire contributions to atmospheric trace gases and aerosols. Important future directions include reconciling top-down and bottom up estimates of fire radiative power and integrating burned area and active fire time series from multiple satellite sensors to improve daily emissions estimates.

Mu, M.↗

SWIPE: Spectral Water Inversion Processor and Emulator

Degradation of Earth’s inland water resources due to anthropogenic perturbations and climate anomalies at both local and global scales continues to place human health at substantial risk. There is now a growing necessity to develop pragmatic approaches that allow timely and effective extrapolation of local processes, to spatially resolved global products, and to promote operational and sustainable resource policy management. This presentation will be discussing the progress made developing SWIPE: Spectral Water Inversion Processor and Emulator. SWIPE is a platform for advanced modeling of coastal and inland aquatic habitats. The goal is create a comprehensive and cohesive system to leverage recent advancements in computation and machine learning to develop a synthetic training ground for sensitivity studies and algorithm development. The four principal facets of SWIPE include: 1. Advanced two-layer coated sphere bio-optical modeling and GPU radiative transfer modeling, 2. Big Data involving massive synthetic spectral libraries of optical properties of various global aquatic particles, surface reflectance, and top-of-atmosphere reflectance, all at hyperspectral resolution leveraging high-end computing systems at NASA Ames Research Center, 3. Deep Learning for algorithm development for water quality inversion of concentrations of common biogeophysical variables as well as optics, full uncertainty characterization by water type, and forward emulation, and lastly, 4. Image Processing for application of developed retrieval algorithms for both hyperspectral and multispectral sensors with experimental corrections for global adjacency, noise, sunglint, and benthic reflectance. This presentation will demonstrate the Equivalent Algal Populations (EAP) two-layer coated sphere scattering model which has been used develop spectral libraries of hyperspectral inherent optical properties of roughly 80 species of phytoplankton, covering 15 different classes and nine taxonomic functional types. The EAP model was also used to derive spectral properties of 10 different non-algal particle functional types. Examples of how the SMART-G (Speed-up Monte-carlo Advanced Radiative Transfer using GPU) radiative transfer code is used to model optically complex aquatic signals will be presented and discussed in the context of creating a massive synthetic database which can leverage the full power of next generation machine learning techniques and high end computing for water quality inversion. We will discuss our active investigation in things like appropriate model architectures, dimensionality reduction techniques such as PCA and autoencoders, uncertainty quantification and abstaining, and which variables actually benefit most from hyperspectral information versus multispectral resolution. We are also curious about questions relating to cost/benefit analysis in terms of computation resources, neural network complexity, and data volumes. Answers to these questions will hopefully elaborate on cost efficiency for potential future sensor design considerations.

SWIPE↗

Expected Performance of Ozone Climate Data Records from Ozone Mapping and Profiler Suite Limb Profiler

The Ozone Mapping and Profiler Suite Limb Profiler (OMPS/LP) was launched on board of the Soumi NPP space platform in late October 2011. It provides ozone-profiling capability with high-vertical resolution from 60 Ian to cloud top. In this study, an end-to-end Observing System Simulation Experiment (OSSE) of OMPS/LP ozone is discussed. The OSSE was developed at NASA's Global Modeling and Assimilation Office (GMAO) using the Goddard Earth Observing System (GEOS-5) data assimilation system. The "truth" for this OSSE is built by assimilating MLS profiles and OMI ozone columns, which is known to produce realistic three-dimensional ozone fields in the stratosphere and upper troposphere. OMPS/LP radiances were computed at tangent points computed by an appropriate orbital model. The OMPS/LP forward RT model, Instrument Models (IMs) and EDR retrieval model were introduced and pseudo-observations derived. The resultant synthetic OMPS/LP observations were evaluated against the "truth" and subsequently these observations were assimilated into GEOS-5. Comparison of this assimilated dataset with the "truth" enables comparisons of the likely uncertainties in 3-D analyses of OMPS/LP data. This study demonstrated the assimilation capabilities of OMPS/LP ozone in GEOS-5, with the monthly, zonal mean (O-A) smaller than 0.02ppmv at all levels, the nns(O-A) close to O.lppmv from 100hPa to 0.2hPa; and the mean(O-B) around the 0.02ppmv for all levels. The monthly zonal mean analysis generally agrees to within 2% of the truth, with larger differences of 2-4% (0.1-0.2ppmv) around 10hPa close to North Pole and in the tropical tropopause region, where the difference is above 20% due to the very low ozone concentrations. These OSSEs demonstrated that, within a single data assimilation system and the assumption that assimilated MLS observations provide a true rendition of the stratosphere, the OMPS/LP ozone data are likely to produce accurate analyses through much of the stratosphere. Additional studies are needed to evaluate the importance of the extremely high vertical resolution of the OMPS/LP ozone data

Xu, P. Q.↗

Version 2 Ozone Monitoring Instrument SO2 Product (OMSO2 V2): New Anthropogenic SO2 Vertical Column Density Dataset

The Ozone Monitoring Instrument (OMI) has been providing global observations of SO2 pollution since 2004. Here we introduce the new anthropogenic SO2 vertical column density (VCD) dataset in the version 2 OMI SO2 product (OMSO2 V2). As with the previous version (OMSO2 V1.3), the new dataset is generated with an algorithm based on principal component analysis of OMI radiances, but features several updates. The most important among those is the use of expanded lookup tables and model a priori profiles to estimate SO2 Jacobians for individual OMI pixels, in order to better characterize pixel-to-pixel variations in SO2 sensitivity, including over snow and ice. Additionally, new data screening and spectral fitting schemes have been implemented to improve the quality of the spectral fit. As compared with the planetary boundary layer SO2 dataset in OMSO2 V1.3, the new dataset has substantially better data quality, especially over areas that are relatively clean or affected by the south Atlantic anomaly. The updated retrievals over snow/ice yield more realistic seasonal changes in SO2 at high latitudes and offer enhanced sensitivity to sources during wintertime. An error analysis has been conducted to assess uncertainties in SO2 VCDs from both the spectral fit and Jacobian calculations. The uncertainties from spectral fitting are reflected in SO2 slant column densities (SCDs) and largely depend on the signal-to-noise ratio of the measured radiances, as implied by the generally smaller SCD uncertainties over clouds or for smaller solar zenith angles. The SCD uncertainties for individual pixels are estimated to be~0.15-0.3 DU (Dobson Units) between ~40°S and ~40°N and to be~0.2-0.5 DU at higher latitudes. The uncertainties from the Jacobians are approximately ~50-100% over polluted areas, and primarily attributed to errors in SO2 a priori profiles and cloud pressures, as well as the lack of explicit treatment for aerosols. Finally, the daily mean and median SCDs over the presumably SO2-free equatorial East Pacific have increased by only~0.0035 DU and ~0.003 DU respectively over the entire 15-year OMI record; while the standard deviation of SCDs has grown by only~0.02 DU or ~10%. Such remarkable long-term stability makes the new dataset particularly suitable for detecting regional changes in SO2 pollution.

OMI, SO2, Remote Sensing↗

TPSAS-NF1676L-10863-DND

One important aspect of successful climate monitoring is satellite calibration. There is at least a 25-year of operational geostationary (GEO) record available for climate studies, if properly calibrated. The Global Space-based Inter-Calibration System (GSICS) project is another effort to provide the climate community consistent calibration across multiple platforms. Also, the Clouds and the Earth's Radiant Energy System (CERES) project, provides consistent GEO retrieved cloud properties and derived broadband fluxes across multiple platforms. This year GSICS is focusing on geostationary visible channel calibration. The GSICS strategy is to use multiple approaches, such as cross-calibrating the GEO radiances against MODIS as a reference, deep convective clouds as bright stable targets, deserts, sun-glint, stars, moon and other model based calibration methods. Each of these methods is independent of each other and can assess the stability of the GEO visible instrument with differing uncertainties. This study focuses on the inter-calibration of the GEO radiances against MODIS as a reference by using collocated, and coincident ray-matched radiances. The advantage of the ray-matching technique is that it provides a full dynamic range of radiances to check the linearity of the GEO visible instrument. This presentation will focus on the technique, some error analysis, comparisons with other methods, and an approach to inter-calibrate historical GEOs before MODIS.

D R Doelling↗

Estimates of Ground Temperature and Atmospheric Moisture from CERES Observations

A method is developed to retrieve surface ground temperature (Tg) and atmospheric moisture using clear sky fluxes (CSF) from CERES-TRMM observations. In general, the clear sky outgoing long-wave radiation (CLR) is sensitive to upper level moisture (q(sub h)) over wet regions and Tg over dry regions The clear sky window flux from 800 to 1200 /cm (RadWn) is sensitive to low level moisture (q(sub j)) and Tg. Combining these two measurements (CLR and RadWn), Tg and q(sub h) can be estimated over land, while q(sub h) and q(sub t) can be estimated over the oceans. The approach capitalizes on the availability of satellite estimates of CLR and RadWn and other auxiliary satellite data. The basic methodology employs off-line forward radiative transfer calculations to generate synthetic CSF data from two different global 4-dimensional data assimilation products. Simple linear regression is used to relate discrepancies in CSF to discrepancies in Tg, q(sub h) and q(sub t). The slopes of the regression lines define sensitivity parameters that can be exploited to help interpret mismatches between satellite observations and model-based estimates of CSF. For illustration, we analyze the discrepancies in the CSF between an early implementation of the Goddard Earth Observing System Data Assimilation System (GEOS-DAS) and a recent operational version of the European Center for Medium-Range Weather Prediction data assimilation system. In particular, our analysis of synthetic total and window region SCF differences (computed from two different assimilated data sets) shows that simple linear regression employing (Delta)Tg and broad layer (Delta)q(sub l) from 500 hPa to surface and (Delta)q(sub h) from 200 to 500 hPa provides a good approximation to the full radiative transfer calculations, typically explaining more than 90% of the 6-hourly variance in the flux differences. These simple regression relations can be inverted to "retrieve" the errors in the geophysical parameters. Uncertainties (normalized by standard deviation) in the monthly mean retrieved parameters range from 7% for (Delta)T to about 20% for (Delta)q(sub t). Our initial application of the methodology employed an early CERES-TRMM data set (CLR and Radwn) to assess the quality of the GEOS2 data. The results showed that over the tropical and subtropical oceans GEOS2 is, in general, too wet in the upper troposphere (mean bias of 0.99 mm) and too dry in the lower troposphere (mean bias of -4.7 mm). We note that these errors, as well as a cold bias in the Tg, have largely been corrected in the current version of GEOS-2 with the introduction of a land surface model, a moist turbulence scheme and the assimilation of SSTM/I total precipitable water.

Wu, Man Li C.↗

The NASA AfriSAR Campaign: Airborne SAR and Lidar Measurements of Tropical Forest Structure and Biomass in Support of Current and Future Space Missions

In 2015 and 2016, the AfriSAR campaign was carried out as a collaborative effort among international space and National Park agencies (ESA, NASA, ONERA, DLR, ANPN and AGEOS) in support of the upcoming ESA BIOMASS, NASA-ISRO Synthetic Aperture Radar (NISAR) and NASA Global Ecosystem Dynamics Initiative (GEDI) missions. The NASA contribution to the campaign was conducted in 2016 with the NASA LVIS (Land Vegetation and Ice Sensor) Lidar, the NASA L-band UAVSAR (Uninhabited Aerial Vehicle Synthetic Aperture Radar). A central motivation for the AfriSAR deployment was the common AGBD estimation requirement for the three future spaceborne missions, the lack of sufficient airborne and ground calibration data covering the full range of ABGD in tropical forest systems, and the intercomparison and fusion of the technologies. During the campaign, over 7000 km2 of waveform Lidar data from LVIS and 30000 km2 of UAVSAR data were collected over 10 key sites and transects. In addition, field measurements of forest structure and biomass were collected in sixteen 1 hectare sized plots. The campaign produced gridded Lidar canopy structure products, gridded aboveground biomass and associated uncertainties, Lidar based vegetation canopy cover profile products, Polarimetric Interferometric SAR and Tomographic SAR products and field measurements. Our results showcase the types of data products and scientific results expected from the spaceborne Lidar and SAR missions; we also expect that the AfriSAR campaign data will facilitate further analysis and use of waveform Lidar and multiple baseline polarimetric SAR datasets for carbon cycle, biodiversity, water resources and more applications by the greater scientific community.

AfriSAR↗

Evaluation of IGS Orbits with Satellite Laser Ranging

The accuracy with which orbits for the Global Positioning System (GPS) spacecraft, can be computed directly affects the accuracy of the resulting site coordinates and polar motion. Several groups routinely analyze GPS ground tracking data to compute precise orbits and terrestrial reference frame solutions. In this paper, we infer the accuracy of the orbits of two of the GPS satellites by comparing to independent laser ranges of subcentimeter accuracy obtained by a small but reasonably well distributed network of tracking sites. We find that all seven International GPS Service for Geodynamics (IGS) analysis centers achieve range residual root mean square (rms) errors at or below the 100 mm level. The best orbit solutions, from JPL, CODE, and the IGS combined product, yield a residual rms of about 50 mm. These residuals are consistent with three dimensional orbit errors of less than 150 mm. Estimating yaw rates for the spacecraft during shadow events, and using these estimates to compute the laser residual, significantly improves the fit. A small mean residual value of -15 to -30 mm seems to exist for most centers and laser sites which is not fully explained at present, but may be due to uncertainties in the corrections to the laser data, such as the reflector to spacecraft center of mass vector or small reference frame differences between the SLR sites and the GPS orbits.

Watkins, M. M.↗

Estimates of Ground Temperature and Atmospheric Moisture from CERES Observations

A method is developed to retrieve surface ground temperature (T(sub g)) and atmospheric moisture using clear sky fluxes (CSF) from CERES-TRMM observations. In general, the clear sky outgoing longwave radiation (CLR) is sensitive to upper level moisture (q(sub l)) over wet regions and (T(sub g)) over dry regions The clear sky window flux from 800 to 1200/cm (RadWn) is sensitive to low level moisture (q(sub t)) and T(sub g). Combining these two measurements (CLR and RadWn), Tg and q(sub h) can be estimated over land, while q(sub h) and q(sub l) can be estimated over the oceans. The approach capitalizes on the availability of satellite estimates of CLR and RadWn and other auxiliary satellite data. The basic methodology employs off-line forward radiative transfer calculations to generate synthetic CSF data from two different global 4-dimensional data assimilation products. Simple linear regression is used to relate discrepancies in CSF to discrepancies in T(sub g), q(sub h) and q(sub l). The slopes of the regression lines define sensitivity parameters that can be exploited to help interpret mismatches between satellite observations and model-based estimates of CSF. For illustration, we analyze the discrepancies in the CSF between an early implementation of the Goddard Earth Observing System Data Assimilation System (GEOS-DAS) and a recent operational version of the European Center for Medium-Range Weather Prediction data assimilation system. In particular, our analysis of synthetic total and window region SCF differences (computed from two different assimilated data sets) shows that simple linear regression employing Delta(T(sub g)) and broad layer Delta(q(sub l) from .500 hPa to surface and Delta(q(sub h)) from 200 to .300 hPa provides a good approximation to the full radiative transfer calculations. typically explaining more than 90% of the 6-hourly variance in the flux differences. These simple regression relations can be inverted to "retrieve" the errors in the geophysical parameters. Uncertainties (normalized by standard deviation) in the monthly mean retrieved parameters range from 7% for Delta(T(sub g)) to about 20% for Delta(q(sub l)). Our initial application of the methodology employed an early CERES-TRMM data set (CLR and Radwn) to assess the quality of the GEOS2 data. The results showed that over the tropical and subtropical oceans GEOS2 is, in general, too wet in the upper troposphere (mean bias of 0.99 mm) and too dry in the lower troposphere (mean bias of -4.7 min). We note that these errors, as well as a cold bias in the T(sub g). have largely been corrected in the current version of GEOS-2 with the introduction of a land surface model, a moist turbulence scheme and the assimilation of SSM/I total precipitable water.

Wu, Man Li C.↗

Simulating PACE Global Ocean Radiances

The NASA PACE mission is a hyper-spectral radiometer planned for launch in the next decade. It is intended to provide new information on ocean biogeochemical constituents by parsing the details of high resolution spectral absorption and scattering. It is the first of its kind for global applications and as such, poses challenges for design and operation. To support pre-launch mission development and assess on-orbit capabilities, the NASA Global Modeling and Assimilation Office has developed a dynamic simulation of global water-leaving radiances, using an ocean model containing multiple ocean phytoplankton groups, particulate detritus, particulate inorganic carbon (PIC), and chromophoric dissolved organic carbon (CDOC) along with optical absorption and scattering processes at 1 nm spectral resolution. The purpose here is to assess the skill of the dynamic model and derived global radiances. Global bias, uncertainty, and correlation are derived using available modern satellite radiances at moderate spectral resolution. Total chlorophyll, PIC, and the absorption coefficient of CDOC (aCDOC), are simultaneously assimilated to improve the fidelity of the optical constituent fields. A 5-year simulation showed statistically significant (P < 0.05) comparisons of chlorophyll (r = 0.869), PIC (r = 0.868), and a CDOC (r =0.890) with satellite data. Additionally, diatoms (r = 0.890), cyanobacteria (r = 0.732), and coccolithophores (r = 0.716) were significantly correlated with in situ data. Global assimilated distributions of optical constituents were coupled with a radiative transfer model (Ocean-Atmosphere Spectral Irradiance Model, OASIM) to estimate normalized water-leaving radiances at 1 nm for the spectral range 250-800 nm. These unassimilated radiances were within 0.074 mW/sq cm/micron/sr of MODIS-Aqua radiances at 412, 443, 488, 531, 547, and 667 nm. This difference represented a bias of 10.4% (model low). A mean correlation of 0.706 (P < 0.05) was found with global distributions of MODIS radiances. These results suggest skill in the global assimilated model and resulting radiances. The reported error characterization suggests that the global dynamical simulation can support some aspects of mission design and analysis. For example, the high spectral resolution of the simulation supports investigations of band selection. The global nature of the radiance representations supports investigations of satellite observing scenarios. Global radiances at bands not available in current and past missions support investigations of mission capability. PACE, ocean color, water-leaving radiances, biogeochemical model, radiative transfer model

PIC↗

Using MERRA Gridded Innovations for Quantifying Uncertainties in Analysis Fields and Diagnosing Observing System Inhomogeneities

MERRA is a NASA reanalysis for the satellite era using a major new version of the Goddard Earth Observing System Data Assimilation System Version 5 (GEOS-5). The project focuses on historical analyses of the hydrological cycle on a broad range of weather and climate time scales and places the NASA EOS suite of observations in a climate context. The characterization of uncertainty in reanalysis fields is a commonly requested feature by users of such data. While intercomparison with reference data sets is common practice for ascertaining the realism of the datasets, such studies typically are restricted to long term climatological statistics and seldom provide state dependent measures of the uncertainties involved. In principle, variational data assimilation algorithms have the ability of producing error estimates for the analysis variables (typically surface pressure, winds, temperature, moisture and ozone) consistent with the assumed background and observation error statistics. However, these "perceived error estimates" are expensive to obtain and are limited by the somewhat simplistic errors assumed in the algorithm. The observation minus forecast residuals (innovations) by-product of any assimilation system constitutes a powerful tool for estimating the systematic and random errors in the analysis fields. Unfortunately, such data is usually not readily available with reanalysis products, often requiring the tedious decoding of large datasets and not so-user friendly file formats. With MERRA we have introduced a gridded version of the observations/innovations used in the assimilation process, using the same grid and data formats as the regular datasets. Such dataset empowers the user with the ability of conveniently performing observing system related analysis and error estimates. The scope of this dataset will be briefly described. We will present a systematic analysis of MERRA innovation time series for the conventional observing system, including maximum-likelihood estimates of background and observation errors, as well as global bias estimates. Starting with the joint PDF of innovations and analysis increments at observation locations we propose a technique for diagnosing bias among the observing systems, and document how these contextual biases have evolved during the satellite era covered by MERRA.

da Silva, Arlindo↗

Observation of Deep Convective Cloud-Top Height and Vertical Temperature Structure of Hurricane Using Hyperspectral Infrared Sounder and its Single-Field-View Retrieval Products

Hurricanes, severe tropical cyclones (TC), or typhoons are significant natural disasters that often result in substantial loss of life and property damage. Numerous studies have indicated that changes in TC intensity are closely linked to deep convective clouds (DCC), with stronger TCs typically exhibiting higher cloud top heights (CTH) compared to weaker TCs. The CTH can help determine if a tropical depression is at the onset of rapid intensification based on case studies. Therefore, accurate determination of TC CTH will be greatly helpful for monitoring TC development and studying TC dynamics. One traditional and most common method to derive CHT from satellite observations is using the thermal brightness temperature in atmospheric channels to match the sounding temperature profile. However, it was found that thermally derived CTH has a lower bias of approximately 1 km, and this bias tends to worsen for the tallest clouds. A new method using the hyperspectral infrared sounder will be presented. From the measurements of Cross-track Infrared Sounder (CrIS) on S-NPP and J-1, along with radiative transfer simulations, we identified the inverted-V spectral feature in the ozone (O3) band (near 9.6 μm) corresponding to high clouds. The depth of the inverted-V can be used to estimate the CTH. Since the depth is computed using the peak absorption O3 channel and the nearby most transparent O3 channel in this O3 band, the uncertainties associated with cloud emissivity and scattering by cloud particles in the traditional method can be ignored. From several hurricane case studies, we found that the CHT derived using this method can accurately capture the structure of the cloud tops in the eyewall, spiral rainbands, and surrounding regions. For example, Hurricane Dorian on September 2, 2019, showed a nicely outward-sloping and circular shape eye cloud in the early morning, but the circular shape of the eyewall cloud became distorted in the afternoon. For various hurricanes we examined, the distribution of CHT for the eyewall clouds differed significantly. To better study the thermodynamic structure of hurricane clouds, this research will analyze the vertical temperature profiles from a new single Field of View (SFOV) Sounder Atmospheric Products (SiFSAP), derived using CrIS and the Advanced Technology Microwave Sounder (ATMS) onboard SNPP and JPSS-1. SiFSAP has a spatial resolution of 15 km at nadir, which surpasses most global weather and climate models and other current operational sounding products. The combined use of ATMS and CrIS allows for retrievals near hurricane eyewalls and spiral rainbands. Wind fields from NASA’s Modern-Era Retrospective Analysis for Research and Applications Version-2 (MERRA-2) and ERA5 will be used to characterize transport, and comparisons between the model temperature and water vapor profiles with the corresponding SiFSAP products will also be provided.

SiFSAP↗

Quantifying the Observability of CO2 Flux Uncertainty in Atmospheric CO2 Records Using Products from Nasa's Carbon Monitoring Flux Pilot Project

NASAs Carbon Monitoring System (CMS) Flux Pilot Project (FPP) was designed to better understand contemporary carbon fluxes by bringing together state-of-the art models with remote sensing datasets. Here we report on simulations using NASAs Goddard Earth Observing System Model, version 5 (GEOS-5) which was used to evaluate the consistency of two different sets of observationally constrained land and ocean fluxes with atmospheric CO2 records. Despite the strong data constraint, the average difference in annual terrestrial biosphere flux between the two land (NASA Ames CASA and CASA-GFED) models is 1.7 Pg C for 2009-2010. Ocean models (NOBM and ECCO2-Darwin) differ by 35 in their global estimates of carbon flux with particularly strong disagreement in high latitudes. Based upon combinations of terrestrial and ocean fluxes, GEOS-5 reasonably simulated the seasonal cycle observed at northern hemisphere surface sites and by the Greenhouse gases Observing SATellite (GOSAT) while the model struggled to simulate the seasonal cycle at southern hemisphere surface locations. Though GEOS-5 was able to reasonably reproduce the patterns of XCO2 observed by GOSAT, it struggled to reproduce these aspects of AIRS observations. Despite large differences between land and ocean flux estimates, resulting differences in atmospheric mixing ratio were small, typically less than 5 ppmv at the surface and 3 ppmv in the XCO2 column. A statistical analysis based on the variability of observations shows that flux differences of these magnitudes are difficult to distinguish from natural variability, regardless of measurement platform.

FPP↗

The Relation Between Atmospheric Humidity and Temperature Trends for Stratospheric Water

We analyze the relation between atmospheric temperature and water vapor-a fundamental component of the global climate system-for stratospheric water vapor (SWV). We compare measurements of SWV (and methane where available) over the period 1980-2011 from NOAA balloon-borne frostpoint hygrometer (NOAA-FPH), SAGE II, Halogen Occultation Experiment (HALOE), Microwave Limb Sounder (MLS)/Aura, and Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS) to model predictions based on troposphere-to-stratosphere transport from ERA-Interim, and temperatures from ERA-Interim, Modern Era Retrospective-Analysis (MERRA), Climate Forecast System Reanalysis (CFSR), Radiosonde Atmospheric Temperature Products for Assessing Climate (RATPAC), HadAT2, and RICHv1.5. All model predictions are dry biased. The interannual anomalies of the model predictions show periods of fairly regular oscillations, alternating with more quiescent periods and a few large-amplitude oscillations. They all agree well (correlation coefficients 0.9 and larger) with observations for higherfrequency variations (periods up to 2-3 years). Differences between SWV observations, and temperature data, respectively, render analysis of the model minus observation residual difficult. However, we find fairly well-defined periods of drifts in the residuals. For the 1980s, model predictions differ most, and only the calculation with ERA-Interim temperatures is roughly within observational uncertainties. All model predictions show a drying relative to HALOE in the 1990s, followed by a moistening in the early 2000s. Drifts to NOAA-FPH are similar (but stronger), whereas no drift is present against SAGE II. As a result, the model calculations have a less pronounced drop in SWV in 2000 than HALOE. From the mid-2000s onward, models and observations agree reasonably, and some differences can be traced to problems in the temperature data. These results indicate that both SWV and temperature data may still suffer from artifacts that need to be resolved in order to answer the question whether the large-scale flow and temperature field is sufficient to explain water entering the stratosphere.

Fueglistaler, S.↗