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 379 records · Page 21

ICESat Observations of Arctic Sea Ice: A First Look

Analysis of near-coincident ICESat and RADARSAT imagery shows that the retrieved elevations from the laser altimeter are sensitive to new openings (containing thin ice or open water) in the sea ice cover as well as to surface relief of old and first-year ice. The precision of the elevation estimates, measured over relatively flat sea ice, is approx. 2 cm. Using the thickness of thin-ice in recent openings to estimate sea level references, we obtain the sea-ice freeboard along the altimeter tracks. This step is necessitated by the large uncertainties in the sea surface topography compared to that required for accurate determination of freeboard. Unknown snow depth introduces the largest uncertainty in the conversion of freeboard to ice thickness. Surface roughness is also derived, for the first time, from the variability of successive elevation estimates along the altimeter track. Overall, these ICESat measurements provide an unprecedented view of the Arctic Ocean ice cover at length scales at and above the spatial dimension of the altimeter footprint of approx. 70 m.

Arctic Ocean↗

Comparative Analysis of Upper Ocean Heat Content Variability from Ensemble Operational Ocean Analyses

Upper ocean heat content (HC) is one of the key indicators of climate variability on many time-scales extending from seasonal to interannual to long-term climate trends. For example, HC in the tropical Pacific provides information on thermocline anomalies that is critical for the longlead forecast skill of ENSO. Since HC variability is also associated with SST variability, a better understanding and monitoring of HC variability can help us understand and forecast SST variability associated with ENSO and other modes such as Indian Ocean Dipole (IOD), Pacific Decadal Oscillation (PDO), Tropical Atlantic Variability (TAV) and Atlantic Multidecadal Oscillation (AMO). An accurate ocean initialization of HC anomalies in coupled climate models could also contribute to skill in decadal climate prediction. Errors, and/or uncertainties, in the estimation of HC variability can be affected by many factors including uncertainties in surface forcings, ocean model biases, and deficiencies in data assimilation schemes. Changes in observing systems can also leave an imprint on the estimated variability. The availability of multiple operational ocean analyses (ORA) that are routinely produced by operational and research centers around the world provides an opportunity to assess uncertainties in HC analyses, to help identify gaps in observing systems as they impact the quality of ORAs and therefore climate model forecasts. A comparison of ORAs also gives an opportunity to identify deficiencies in data assimilation schemes, and can be used as a basis for development of real-time multi-model ensemble HC monitoring products. The OceanObs09 Conference called for an intercomparison of ORAs and use of ORAs for global ocean monitoring. As a follow up, we intercompared HC variations from ten ORAs -- two objective analyses based on in-situ data only and eight model analyses based on ocean data assimilation systems. The mean, annual cycle, interannual variability and longterm trend of HC have been analyzed

Xue, Yan↗

A preliminary evaluation of ocean topography from the TOPEX/POSEIDON mission

We have analyzed 50 ten-day cycles of TOPEX/POSEIDON (T/P) altimeter data to evaluate the ocean dynamic topography and its temporal variations. We have employed data from both the U.S. and French altimeters along with the NASA precision orbits in this analysis. Errors in the diurnal and semidiurnal components of the Cartwright-Ray tide model have been significantly reduced using a correction developed from a combination of JGM-2 and OSU91A was employed, as well as a geoid model based solely on OSU91A. The long wavelengths of the comparisons to historical data, although geoid error still corrupts the dynamic topography for wavelengths shorter than 2500 km. The root mean square (RMS) variability is similar to previous results from Geosat, with bakground 'noise' approaching 3 cm RMS. The computed annual and semiannual variations are also similar to previous Geosat results, although the hemispheric distribution of the annual heating cycle is much better presented in the T/P results. They also compare reasonably well with the Levitus hydrographic compilation in the northern hemisphere, although the T/P variations generally have larger amplitudes. Ten-day average maps of variations in sea level compare well with simulations measurements at ocean tide gauges, with RMS differences of less than 4 cm and correlations greater than 0.6 for most of the island gauges. Time-longitude plots of these sea level variations at different latitudes in the Pacific clearly show the presence of equatorial Kelvin waves and Rossby waves, with the wave speeds agreeing well with theoretical and observed values. Measurement of variations in global sea level over cycles 2-51 have an RMS variability of 6.3 mm and a rate of change of -3.5 +/- 8 mm/yr, the uncertainty primarily due to insufficient averaging of the interannual and periodic sea level variations. These results show that the accuracy of the T/P measurements of sea level has dramatically improved over previous missions, with estimated time variable errors of 4 cm or less. Although geographically correlated orbit errors have also been reduced to the few centimeter level, further improvement in determinations of the mean dynamic topography will be difficult to obtain until a more accurate model of the marine geoid is available.

Nerem, R. S.↗

Variability in Global Top-of-Atmosphere Shortwave Radiation Between 2000 and 2005

Measurements from various instruments and analysis techniques are used to directly compare changes in Earth-atmosphere shortwave (SW) top-of-atmosphere (TOA) radiation between 2000 and 2005. Included in the comparison are estimates of TOA reflectance variability from published ground-based Earthshine observations and from new satellite-based CERES, MODIS and ISCCP results. The ground-based Earthshine data show an order-of-magnitude more variability in annual mean SW TOA flux than either CERES or ISCCP, while ISCCP and CERES SW TOA flux variability is consistent to 40%. Most of the variability in CERES TOA flux is shown to be dominated by variations global cloud fraction, as observed using coincident CERES and MODIS data. Idealized Earthshine simulations of TOA SW radiation variability for a lunar-based observer show far less variability than the ground-based Earthshine observations, but are still a factor of 4-5 times more variable than global CERES SW TOA flux results. Furthermore, while CERES global albedos exhibit a well-defined seasonal cycle each year, the seasonal cycle in the lunar Earthshine reflectance simulations is highly variable and out-of-phase from one year to the next. Radiative transfer model (RTM) approaches that use imager cloud and aerosol retrievals reproduce most of the change in SW TOA radiation observed in broadband CERES data. However, assumptions used to represent the spectral properties of the atmosphere, clouds, aerosols and surface in the RTM calculations can introduce significant uncertainties in annual mean changes in regional and global SW TOA flux.

Loebe, Norman G.↗

The Flux-Anomaly-Forced Model Intercomparison Project (FAFMIP) Contribution to CMIP6: Investigation of Sea-Level and Ocean Climate Change in Response to CO2 Forcing

The Flux-Anomaly-Forced Model Intercomparison Project (FAFMIP) aims to investigate the spread in simulations of sea-level and ocean climate change in response to CO2 forcing by atmosphere-ocean general circulation models (AOGCMs). It is particularly motivated by the uncertainties in projections of ocean heat uptake, global-mean sealevel rise due to thermal expansion and the geographical patterns of sea-level change due to ocean density and circulation change. FAFMIP has three tier-1 experiments, in which prescribed surface flux perturbations of momentum, heat and freshwater respectively are applied to the ocean in separate AOGCM simulations. All other conditions are as in the pre-industrial control. The prescribed fields are typical of pattern and magnitude of changes in these fluxes projected by AOGCMs for doubled CO2 concentration. Five groups have tested the experimental design with existing AOGCMs. Their results show diversity in the pattern and magnitude of changes, with some common qualitative features. Heat and water flux perturbation cause the dipole in sea-level change in the North Atlantic, while momentum and heat flux perturbation cause the gradient across the Antarctic Circumpolar Current. The Atlantic meridional overturning circulation (AMOC) declines in response to the heat flux perturbation, and there is a strong positive feedback on this effect due to the consequent cooling of sea-surface temperature in the North Atlantic, which enhances the local heat input to the ocean. The momentum and water flux perturbations do not substantially affect the AMOC. Heat is taken up largely as a passive tracer in the Southern Ocean, which is the region of greatest heat input, while the weakening of the AMOC causes redistribution of heat towards lower latitudes. Future analysis of these and other phenomena with the wider range of CMIP6 FAFMIP AOGCMs will benefit from new diagnostics of temperature and salinity tendencies, which will enable investigation of the model spread in behaviour in terms of physical processes as formulated in the models.

thermosteric sea-level change↗

Detecting a Terrestrial Biosphere Sink for Carbon Dioxide: Interannual Ecosystem Modeling for the Mid-1980s

There is considerable uncertainty as to whether interannual variability in climate and terrestrial ecosystem production is sufficient to explain observed variation in atmospheric carbon content over the past 20-30 years. In this paper, we investigated the response of net CO2 exchange in terrestrial ecosystems to interannual climate variability (1983 to 1988) using global satellite observations as drivers for the NASA-CASA (Carnegie-Ames-Stanford Approach) simulation model. This computer model of net ecosystem production (NEP) is calibrated for interannual simulations driven by monthly satellite vegetation index data (NDVI) from the NOAA Advanced Very High Resolution Radiometer (AVHRR) at 1 degree spatial resolution. Major results from NASA-CASA simulations suggest that from 1985 to 1988, the northern middle-latitude zone (between 30 and 60 degrees N) was the principal region driving progressive annual increases in global net primary production (NPP; i.e., the terrestrial biosphere sink for carbon). The average annual increase in NPP over this predominantly northern forest zone was on the order of +0.4 Pg (10 (exp 15) g) C per year. This increase resulted mainly from notable expansion of the growing season for plant carbon fixation toward the zonal latitude extremes, a pattern uniquely demonstrated in our regional visualization results. A net biosphere source flux of CO2 in 1983-1984, coinciding with an El Nino event, was followed by a major recovery of global NEP in 1985 which lasted through 1987 as a net carbon sink of between 0.4 and 2.6 Avg C per year. Analysis of model controls on NPP and soil heterotrophic CO2 fluxes (Rh) suggests that regional warming in northern forests can enhance ecosystem production significantly. In seasonally dry tropical zones, periodic drought and temperature drying effects may carry over with at least a two-year lag time to adversely impact ecosystem production. These yearly patterns in our model-predicted NEP are consistent in magnitude with the estimated exchange of CO2 by the terrestrial biosphere with the atmosphere, as determined by previous isotopic (delta (sup 13 C) convolution analysis. Ecosystem simulation results can help further target locations where net carbon sink fluxes have occurred in the past or may be verified in subsequent field studies.

Potter, Christopher S.↗

The Geoscience Laser Altimeter System (GLAS) for the ICESAT Mission

Accurate measurements of surface heights and atmospheric backscatter have been demonstrated with the SLA, MOLA and LITE space lidar. Recent MOLA measurements of the Mars surface have 40 cm resolution and have reduced the global uncertainty in Mars topography from a few km to approx. 10 m. GLAS is a next generation lidar being developed as part of NASA's Icesat Mission for Earth orbit . The GLAS design combines a 10 cm precision surface lidar with a sensitive dual wavelength cloud and aerosol lidar. GLAS will precisely measure the heights of the Earth's polar ice sheets, determine the height profiles of the Earth's land topography, and profile the vertical backscatter of clouds and aerosols on a global scale. GLAS will fly on a small dedicated spacecraft in a polar orbit at 598 km altitude with an inclination of 94 degrees. GLAS is scheduled to launch in summer 2001 and to operate continuously for a minimum of 3 years with a goal of 5 years. The primary mission for GLAS is to measure the seasonal and annual changes in the heights of the Greenland and Antarctic ice sheets. GLAS will measure the vertical distance to the ice sheet from orbit with 1064 nm pulses from a Nd:Yag laser at 40 Hz. Each 5 nsec wide laser pulse is used for a single range measurement. When over land GLAS will profile the heights of the topography and vegetation. The GLAS receiver uses a I m diameter telescope and a Si APD detector. The detector signal is sampled by an all digital receiver which records each surface echo waveform with I nsec resolution and a stored echo record lengths of either 200, 400, or 600 samples. Analysis of the echo waveforms within the instrument permits discrimination between cloud and surface echoes. Ground based echo analysis permits precise ranging, determining the roughness or slopes of the surface as well as the vertical distributions of vegetation illuminated by the laser, Errors in knowledge of the laser beam pointing angle can bias height measurements of sloped surfaces. For surfaces with 2 deg. slopes, knowledge of pointing angle of the beam centroid to about 8 urad is required to achieve 10 cm height accuracy. GLAS uses a stellar reference system (SRS) to determine the pointing angle of each laser firing relative to inertial space. The SRS uses a high precision star camera oriented toward local zenith whose measurements are combined with a gyroscope to determine the inertial orientation of the SRS optical bench. The far field pattern of each laser pulse is measured with a laser reference system (LRS). Optically measuring each laser far field pattern relative to the star camera and gyroscope permits the angular offsets of each laser pulse to be determined. GLAS will also determine the vertical distributions of clouds and aerosols by measuring atmospheric backscatter profiles at both 1064 and 532 nm. The 1064 nm measurements use an analog detector and profile the height and vertical structure of thicker clouds. Measurements at 532 nm use new highly sensitive photon counting detectors, and measure the height distributions of very thin clouds and aerosol layers. With averaging these can be used to determine the height of the planetary boundary layer. The instrument design and expected performance will be discussed.

Abshire, James B.↗

Evaluation of long-term Northern Hemisphere snow water equivalent products

Nine gridded Northern Hemisphere snow water equivalent (SWE) products were evaluated as part of the European Space Agency (ESA) Satellite Snow Product Intercomparison and Evaluation Exercise (SnowPEx). Three categories of datasets were assessed: (1) those utilizing some form of reanalysis (the NASA Global Land Data Assimilation System version 2 – GLDAS-2; the European Centre for Medium-Range Weather Forecasts (ECMWF) interim land surface reanalysis – ERA-Interim/Land and ERA5; the NASA Modern-Era Retrospective Analysis for Research and Applications version 1 (MERRA) and version 2 (MERRA-2); the Crocus snow model driven by ERA-Interim meteorology – Crocus); (2) passive microwave remote sensing combined with daily surface snow depth observations (ESA GlobSnow v2.0); and (3) stand-alone passive microwave retrievals (NASA AMSR-E SWE versions 1.0 and 2.0) which do not utilize surface snow observations. Evaluation included validation against independent snow course measurements from Russia, Finland, and Canada and product intercomparison through the calculation of spatial and temporal correlations in SWE anomalies. The stand-alone passive microwave SWE products (AMSR-E v1.0 and v2.0 SWE) exhibit low spatial and temporal correlations to other products and RMSE nearly double the best performing product. Constraining passive microwave retrievals with surface observations (GlobSnow) provides performance comparable to the reanalysis-based products; RMSE over Finland and Russia for all but the AMSR-E products is ∼50 mm or less, with the exception of ERA-Interim/Land over Russia. Using a seven-dataset ensemble that excluded the stand-alone passive microwave products reduced the RMSE by 10 mm (20 %) and increased the correlation from 0.67 to 0.78 compared to any individual product. The overall performance of the best multiproduct combinations is still at the margins of acceptable uncertainty for scientific and operational requirements; only through combined and integrated improvements in remote sensing, modeling, and observations will real progress in SWE product development be achieved.

Northern Hemisphere snow water equivalent (SWE) pr↗

Assessing the Variability of Radiation, Water and Energy in the Deep Tropics Over the Last 3 Decades

Global surface energy closure and its variability depends heavily upon the surface radiative energy budget in the deep tropics. A key process observed is the change in convection regimes from disorganized convection to organized deep convection. During these periods an imbalance is observed between energy and water fluxes in terms of energy closure in some atmospheric reanalysis systems, implying that processes are not simulated well. A recent paper by Hsaio et al (2023, pre-published version) found that the transition process has wind shear and longwave cloud radiative feedback signatures. In this presentation, we assess the surface radiative budget in the deep tropics by comparing multiple data products (GEWEX SRB Rel4IP, CERES Ed4.1, ISSCP FH, etc.) and describing the variability across the deep tropics for the period from 1988 to near present. We assess this variability of the radiative flux anomalies (including the net TOA, surface and the atmospheric divergence fluxes) against water vapor divergence, cloud properties (include ISCCP “Weather States, Tselioudis et al. 2021) and larger scale wind shear. We include further analysis contrasting 6 key tropical oceanic regions (Indian Ocean east and west, tropical western, central and eastern Pacific Ocean, and Tropical Atlantic Ocean) and the sensitivity of the fluxes as a function of fluctuations in cloud types in response to various larger scale atmospheric processes (e.g. El Niño, Indian Ocean Dipole). These are contrasted against the modeled radiative fluxes from ERA-5 and MERRA-2. Ocean buoy measurements of radiative fluxes are utilized to help assess data radiative flux uncertainty over the nearly 40 year period. In general, the variability for overlap periods of these various data products agrees well, but there are significant differences in the net fluxes that vary according to the rendering of surface, atmospheric, aerosol and cloud properties. We conclude with recommendations for continuing work.

Radiation budget↗

Assimilation of Sea Ice Thickness Derived from Cryosat-2 Along-Track Freeboard Measurements into the Met Office's Forecast Ocean Assimilation Model (FOAM)

The feasibility of assimilating sea ice thickness (SIT) observations derived from CryoSat-2 along-track measurements of sea ice freeboard is successfully demonstrated using a 3D-Var assimilation scheme, NEMOVAR, within the Met Office's global, coupled ocean–sea-ice model, Forecast Ocean Assimilation Model (FOAM). The CryoSat-2 Arctic freeboard measurements are produced by the Centre for Polar Observation and Modelling (CPOM) and are converted to SIT within FOAM using modelled snow depth. This is the first time along-track observations of SIT have been used in this way, with other centres assimilating gridded and temporally averaged observations. The assimilation leads to improvements in the SIT analysis and forecast fields generated by FOAM, particularly in the Canadian Arctic. Arctic-wide observation-minus-background assimilation statistics for 2015–2017 show improvements of 0.75 m mean difference and 0.41 m root-mean-square difference (RMSD) in the freeze-up period and 0.46 m mean difference and 0.33 m RMSD in the ice break-up period. Validation of the SIT analysis against independent springtime in situ SIT observations from NASA Operation IceBridge (OIB) shows improvement in the SIT analysis of 0.61 m mean difference (0.42 m RMSD) compared to a control without SIT assimilation. Similar improvements are seen in the FOAM 5 d SIT forecast. Validation of the SIT assimilation with independent Beaufort Gyre Exploration Project (BGEP) sea ice draft observations does not show an improvement, since the assimilated CryoSat-2 observations compare similarly to the model without assimilation in this region. Comparison with airborne electromagnetic induction (Air-EM) combined measurements of SIT and snow depth shows poorer results for the assimilation compared to the control, despite covering similar locations to the OIB and BGEP datasets. This may be evidence of sampling uncertainty in the matchups with the Air-EM validation dataset, owing to the limited number of observations available over the time period of interest. This may also be evidence of noise in the SIT analysis or uncertainties in the modelled snow depth, in the assimilated SIT observations, or in the data used for validation. The SIT analysis could be improved by upgrading the observation uncertainties used in the assimilation. Despite the lack of CryoSat-2 SIT observations available for assimilation over the summer due to the detrimental effect of melt ponds on retrievals, it is shown that the model is able to retain improvements to the SIT field throughout the summer months due to prior, wintertime SIT assimilation. This also results in regional improvements to the July modelled sea ice concentration (SIC) of 5 % RMSD in the European sector, due to slower melt of the thicker sea ice.

Emma K. Fiedler↗

Reduction of Aircraft Noise Uncertainty for a Notional Supersonic Business Jet

NASA supported a study by the International Civil Aviation Organization on the environmental impact of adding supersonic aircraft to the existing global fleet by designing a notional 55-tonne supersonic business jet. The aircraft is referred to as the Supersonic Technology Concept Aeroplane (STCA) and has been used in a multitude of studies over the past few years. One of the many studies on this aircraft was the noise produced by the STCA during landing and take-off (LTO) operations. The LTO noise for the STCA was assessed using contemporary noise prediction tools that have been primarily developed for and utilized by subsonic aircraft, resulting in a high level of uncertainty for the STCA noise predictions. NASA has recently been conducting research and developing tools to reduce the uncertainty for airport noise predictions of supersonic aircraft. The initial focus of the uncertainty reduction was on the jet and inlet-radiated fan noise produced by supersonic engine systems since these sources are dominant during take-off and landing operations for supersonic jets. The results of these efforts are discussed and applied to the STCA model to update the airport noise predictions and associated uncertainty metrics. It is shown that by using new noise source prediction models, the overall system-level cumulative noise uncertainty of the study vehicle is reduced from a standard deviation of 7.8 EPNdB to 2.0 EPNdB.

uncertainty↗

Composition Reanalyses at NASA’s Global Modeling and Assimilation Office

Scientific Assessments of Ozone Depletion 2018 and 2022 report the recovery of Antarctic ozone, consistent with diminishing concentrations of halogens in the stratosphere. Both the size and the depth of the springtime ozone holes exhibit downward trends, particularly in September. Superimposed on the recovery trend are significant year-to-year fluctuations arising from dynamical variability potentially influenced by climate forcing and volcanic eruptions. Since 2018 the Antarctic stratosphere has experienced a wild parade of austral spring seasons with a sudden stratospheric warming and exceptionally high ozone in 2019 followed by large and long-lasting ozone holes in the following years, prompting questions about the speed and detectability of Antarctic ozone recovery and its potential connections to human-induced climate change. From the Australian New Years’ wildfires of 2020 to volcanic aerosols to the unprecedented injection of water vapor by the Hunga volcanic eruption, recent scientific literature on the subject investigates many factors that may or may not have contributed to this remarkable interannual variability. This presentation summarizes dynamical influences on springtime Antarctic stratospheric ozone in recent years and discusses these results in the context of long-term trends. In particular, this work seeks to isolate dynamical and chemical contributions to the Antarctic springs since 2019 using a consistent methodological framework. Our analysis uses constituent and meteorological fields from the MERRA-2 Stratospheric Composition Reanalysis of Aura MLS (M2-SCREAM) in combination with data from the Atmospheric Chemistry Experiment – Fourier Transform Spectrometer and from ozone sondes. Using the high spatial and temporal resolution of M2-SCREAM, we analyze the daily evolution of southern polar ozone, HCl, N2O, and water vapor in dynamical flow-following coordinates relative to the location, size, and shape of the polar vortex. For vertically integrated stratospheric ozone we emphasize the importance of the wave-influenced geometry of the polar vortex and its vertical alignment. The ultimate purpose of this project is to contribute to the quantification of the speed of the Antarctic ozone recovery and associated uncertainties.

Krzysztof Wargan↗

A USCLIVAR Project to Assess and Compare the Responses of Global Climate Models to Drought-Related SST Forcing Patterns: Overview and Results

The USCLI VAR working group on drought recently initiated a series of global climate model simulations forced with idealized SST anomaly patterns, designed to address a number of uncertainties regarding the impact of SST forcing and the role of land-atmosphere feedbacks on regional drought. Specific questions that the runs are designed to address include: What are the mechanisms that maintain drought across the seasonal cycle and from one year to the next? What is the role of the leading patterns of SST variability, and what are the physical mechanisms linking the remote SST forcing to regional drought, including the role of land-atmosphere coupling? The runs were carried out with five different atmospheric general circulation models (AGCM5), and one coupled atmosphere-ocean model in which the model was continuously nudged to the imposed SST forcing. This paper provides an overview of the experiments and some initial results focusing on the responses to the leading patterns of annual mean SST variability consisting of a Pacific El Nino/Southern Oscillation (ENSO)-like pattern, a pattern that resembles the Atlantic Multi-decadal Oscillation (AMO), and a global trend pattern. One of the key findings is that all the AGCMs produce broadly similar (though different in detail) precipitation responses to the Pacific forcing pattern, with a cold Pacific leading to reduced precipitation and a warm Pacific leading to enhanced precipitation over most of the United States. While the response to the Atlantic pattern is less robust, there is general agreement among the models that the largest precipitation response over the U.S. tends to occur when the two oceans have anomalies of opposite sign. That is, a cold Pacific and warm Atlantic tend to produce the largest precipitation reductions, whereas a warm Pacific and cold Atlantic tend to produce the greatest precipitation enhancements. Further analysis of the response over the U.S. to the Pacific forcing highlights a number of noteworthy and to some extent unexpected results. These include a seasonal dependence of the precipitation response that is characterized by signal-to-noise ratios that peak in spring, and surface temperature signal-to-noise ratios that are both lower and show less agreement among the models than those found for the precipitation response. Another interesting result concerns what appears to be a substantially different character in the surface temperature response over the U.S. to the Pacific forcing by the only model examined here that was developed for use in numerical weather prediction. The response to the positive SST trend forcing pattern is an overall surface warming over the world's land areas with substantial regional variations that are in part reproduced in runs forced with a globally uniform SST trend forcing. The precipitation response to the trend forcing is weak in all the models.

Schubert, Siegfried↗

Chemical Ionization Mass Spectrometry-Based Measurements of HO2 and RO2 During TRACE-P

The Transport and Chemical Evolution over the Pacific (TRACE-P) mission extends NASA's Global Tropospheric Experiment (GTE) series of campaigns. TRACE-P was an aircraft-based campaign that was part of a larger ground-based and aircraft-based program (APARE) under the auspices of the International Global Atmospheric Chemistry (IGAC) program. TRACE-P was designed to (1) determine the chemical composition of Asian outflow over the western Pacific, and to (2) determine the chemical evolution of the Asian outflow. These objectives were addressed through a variety of observations and numerical modeling exercises. In particular, the goals included sampling strategies that would improve understanding of the budgets of odd hydrogen species (OH and HO2), budgets of NOx (NO, NO2, and their reservoirs), and impacts of oxidants produced in the outflow on air quality in the United States. The NASA DC-8 and P-3B aircraft were deployed in the March and April, 2001 out of primary bases of operation in Hong Kong and Yokota Air Base in Japan. These two aircraft have complementary capabilities which allow high altitude and long range impacts, as well as low altitude, local impacts to be assessed. In order to quantify the composition and evolution of Asian outflow, it is important to quantify as many species as possible including photochemically active species (e.g. NO2, CH2O, O3, acetone, etc.), sources species (VOCs, CO, NOx, SO2, aerosols), reactive intermediates including free radicals (OH, HO2, RO2, and their reservoirs), and end products (nitric acid, sulfuric acid, secondary aerosols, etc.). The more complete the measurement suite, the more tightly constrained the numerical modeling can be (within the uncertainties of the measurements). The numerical models range in sophistication from simple steady state box models (as employed in this study) to multi-dimensional chemical transport models. Data were collected on approximately 20 flights of the DC-8 and 21 flights of the P-3B. Observations from both aircraft were used in the present analysis, but primarily focused on the P-3B flights since that was the platform on which the peroxy radical instrumentation was based.

Cantrell, Christopher A.↗

The Geoscience Laser Altimeter System (GLAS) for the ICESAT Mission

The Laser In space Technology Experiment, Shuttle Laser Altimeter and the Mars Observer Laser Altimeter have demonstrated accurate measurements of atmospheric backscatter and Surface heights from space. The recent MOLA measurements of the Mars surface have 40 cm vertical resolution and have reduced the global uncertainty in Mars topography from a few km to about 5 m. The Geoscience Laser Altimeter System (GLAS) is a next generation lidar for Earth orbit being developed as part of NASA's Icesat Mission. The GLAS design combines a 10 cm precision surface lidar with a sensitive dual wavelength cloud and aerosol lidar. GLAS will precisely measure the heights of the Earth's polar ice sheets, establish a grid of accurate height profiles of the Earth's land topography, and profile the vertical backscatter of clouds and aerosols on a global scale. GLAS is being developed to fly on a small dedicated spacecraft in a polar orbit with a 590 630 km altitude at inclination of 94 degrees. GLAS is scheduled to launch in the summer 2001 and to operate continuously for a minimum of 3 years with a goal of 5 years. The primary mission for GLAS is to measure the seasonal and annual changes in the heights of the Greenland and Antarctic ice sheets. GLAS will continuously measure the vertical distance from orbit to the Earth's surface with 1064 nm pulses from a ND:YAG laser at a 40 Hz rate. Each 5 nsec wide laser pulse is used to produce a single range measurement, and the laser spots have 66 m diameter and about 170 m center-center spacings. When over land GLAS will profile the heights of the topography and vegetation. The GLAS receiver uses a 1 m diameter telescope and a Si APD detector. The detector signal is sampled by an all digital receiver which records each surface echo waveform with I nsec resolution and a stored echo record lengths of either 200, 400, or 600 samples. Analysis of the echo waveforms within the instrument permits discrimination between cloud and surface echoes. Ground based echo analysis permits precise ranging, determining the roughness or slopes of the surface as well as the vertical distributions of vegetation illuminated by the laser. Accurate knowledge of the laser beam's pointing angle is needed to prevent height biases when over sloped surfaces. For surfaces with 2 deg. slopes, knowledge of pointing angle of the beam's centroid to about 8 urad is needed to achieve 10 cm height accuracy. GLAS uses a stellar reference system (SRS) to determine the pointing angle of each laser firing relative to inertial space. The SRS uses a high precision star camera oriented toward local zenith and a gyroscope to determine the inertial orientation of the SRS optical bench. The far field pattern of each laser is measured pulse relative to the star camera with a laser reference system (LRS). Optically measuring each laser far field pattern relative to the orientation of the star camera and gyroscope permits the precise pointing angle of each laser pulse to be determined. GLAS will also determine the vertical distributions of clouds and aerosols by measuring the vertical profile of laser energy backscattered by the atmosphere at both 1064 and 532 nm. The 1064 nm measurements use the Si APD detector and profile the height and vertical structure of thicker clouds. The measurements at 532 nm use new highly sensitive photon counting, detectors, and measure the height distributions of very thin Clouds and aerosol layers. With averaging these can be used to determine the height of the planetary boundary layer. The instrument design and expected performance will be discussed.

Abshire, James B.↗

Understanding the Uncertainty and Variability of Earth’s Radiation Budget Components Over the Past 3 Decades from the NASA’s GEWEX Surface Radiation Budget Release 4 Integrated Product

The NASA/GEWEX Surface Radiation Budget (SRB) project produced a 30+ year time series of 3-hourly Top-of-Atmosphere (TOA) and surface radiative flux components at a 1x1 degree resolution. This SRB Release 4 Integrated Productare derived using imager data products from the cross-calibrated from ISCCP (International Satellite Cloud Climatology Project) “HSeries” radiance, cloud and surface property data products. The ISCCP HXS product replaces ISCCP DX with a ninefold increase in pixel count (10km instead of 30km). ISCCP also provides an atmospheric temperature and moisture dataset known as nnHIRS which we use and discuss radiative flux sensitivities to in this presentation. Total solar irradiance is now variable consistent with SORCE measurements. Sea surface (SeaFlux v2.0) and land surface products (LandFlux) are utilized to achieve an radiative flux estimatesintegrated with global data sets from these GEWEX projects. The variable aerosol composition are specified using a detailed aerosol history from the Max Planck Institute Aerosol Climatology (MAC) also used for background aerosol in the ISCCP HX. Numerous other improvements were made to the SW and LW algorithms including improved cloud optical property treatments and specification of surface spectral albedo and emissivity for various surface conditions. Here we analyze the new SWand LWSRB datasets, comparing them to the previous Release 3, BSRN, GEBA and PMEL surface measurements, and ERBE and CERES satellite datasets. For surface flux validation besides ensemble comparisons, we show the variability of SRB vs surface measurements from BSRN beginning in 1992 and GEBA SW from 1983. We assess the variability of the surface flux anomalies, particularly the net fluxes and the cloud radiative effects, over the entire time period including the period containing CERES overlap. For the TOA fluxes, the anomalies of the fluxes in the early period are compared to the MEaSURES ERBE and the latter period to CERES EBAF/SYN1Deg datasets for a benchmark. Long-term changes in the surface radiation budget components and cloud radiative effects are shown and discussed relative to CERES and surface measurements. This includes a change point analysis that helps associate various changes in the ISCCP GEO/LEO configurations and other inputs to variations in the radiative budget components.

Earth Radiation Budget↗

NASA Earth Systems Digital Twins (ESDT)

"Similarly to artificial intelligence, which is now revolutionizing many aspects of our daily lives, Earth system digital twin technologies have the potential to revolutionize the way Earth Science research will be conducted in the future, and how results and knowledge from this research will provide information to support decision making and yield impactful societal benefits. An Earth System Digital Twin or ESDT is a dynamic and interactive information system that first provides a digital replica of the past and current states of the Earth or Earth system as accurately and timely as possible; second, allows for computing forecasts of future states under nominal assumptions and based on the current replica; and third, offers the capability to investigate many hypothetical scenarios under varying impact assumptions. In other words, an ESDT provides the integrated What-Now, What-Next, and What-If pictures of the Earth or Earth system, by continuously ingesting newly observed data and by leveraging multiple interconnected models, machine learning as well advanced computing and visualization capabilities. Digital twins have been developed in engineering since 2002, but the interest in digital twins for the Earth domain is more recent and stems from the convergence of several developments: - The huge amount of diverse data that has now been collected continuously for more than 50 years, and that is becoming more and more difficult to access, understand, and utilize. - At the same time, because of climate change and its impacts the information produced by all of this data is becoming of interest to many new non-traditional users for analyzing and predicting various phenomena. - Because of advances in computational and visualization capabilities and the parallel unprecedented development of machine learning (ML), extracting relevant information from these large amounts of data and running complex models faster has become possible. As a result, it is becoming necessary and possible to build intuitive and interactive frameworks that will enable users with various skill levels and/or organizational hierarchy levels to easily access large amounts of targeted information along with the relevant tools and models (Earth system and human activity models), to support them in analyzing and visualizing this information, to help them understand interactions among models, to visualize the potential outcomes of various impacts, and to support decision or policy making. The full power of digital twins is that, through an integrated representation and standardized tools and software technologies, the same digital replica can address the needs of multiple users at various resolutions (spatial and temporal) and for various applications (science, economic, policy, etc.) – “from farmer to scientist”. With all these interests at stake, the challenges of building optimal digital twins are many and complex. The first challenge is to determine if a Digital Twin should be global or local, and multi-domain or thematic. For example, some domains such as Climate or Weather will require a global Digital Twin or Digital Twin capabilities while science areas such as Biodiversity might be more local. We can also envision that multiple thematic ESDTs, e.g., Air Quality, Wildfires, Hydrology could be federated or provide input to other ESDTs, either on a regional level or to a more global ESDT. Overall, we can imagine a future “web” of Digital Twins co-existing in a hierarchy or in a network, and capable of being connected or federated depending on the needs. This last point brings up the very important challenge of interoperability, including standards and protocols that will need to be built into these systems from the beginning. Each individual digital twin would have full flexibility in internal construction but would need standards-based interfaces (input and output) or hooks to make it compatible with others. Another challenge when building digital twins will be to decide how to organize each digital replica. Based on the applications targeted by the DT under implementation, various amounts and types of raw data, Analysis Ready Data (ARD) and information will need to be incorporated. Depending on the required latencies and needs of the users, various solutions can be considered, including Data Cubes, Data Lakes, pointers, or computing information on demand. We envision that each ESDT will choose a solution adapted to its specific objectives. Another important challenge is the type(s) of visualization that will be used, as well as the level of interactivity and refresh rate that will be required. Again, this will depend on the objectives of the ESDT, but also on the various users’ needs. In most cases, several types of visualizations and human interfaces will need to be offered depending on the projected users of that system. In parallel to the challenges highlighted above, there are also many tools and technologies that will need to be developed or improved for all types of digital twins. Among those are improved machine learning technologies, for example providing explainability, but also ML techniques for causality and providing a better integration of physics models. Additionally, reliable uncertainty quantification methods will be needed for all ESDT components, from validating data fusion and assimilation to assessing the accuracy of ML models and weighing the values of decisions supported by those systems. This presentation introduces the ESDT concept, presents several ESDT use cases, and a proposed ESDT architecture framework, as well as various technologies being developed by the Advanced Information Systems Technology (AIST) Program."

Earth Science Remote Sensing; Information Systems↗

Forecasting Ocean Chlorophyll in the Equatorial Pacific

Using a global ocean biogeochemical model combined with a forecast of physical oceanic and atmospheric variables from the NASA Global Modeling and Assimilation Office, we assess the skill of a chlorophyll concentrations forecast in the Equatorial Pacific for the period 2012-2015 with a focus on the forecast of the onset of the 2015 El Niño event. Using a series of retrospective 9-month hindcasts, we assess the uncertainties of the forecasted chlorophyll by comparing the monthly total chlorophyll concentration from the forecast with the corresponding monthly ocean chlorophyll data from the Suomi-National Polar-orbiting Partnership Visible Infrared Imaging Radiometer Suite (S-NPP VIIRS) satellite. The forecast was able to reproduce the phasing of the variability in chlorophyll concentration in the Equatorial Pacific, including the beginning of the 2015-2016 El Niño. The anomaly correlation coefficient (ACC) was significant (p less than 0.05) for forecast at 1-month (R=0.33), 8-month (R=0.42) and 9-month (R=0.41) lead times. The root mean square error (RMSE) increased from 0.0399 microgram ch1 L(exp -1) for the 1-month lead forecast to a maximum of 0.0472 microgram ch1 L(exp -1) for the 9-month lead forecast indicating that the forecast of the amplitude of chlorophyll concentration variability was getting worse. Forecasts with a 3-month lead time were on average the closest to the S-NPP VIIRS data (23% or 0.033 microgram ch1 L(exp -1)) while the forecast with a 9-month lead time were the furthest (31% or 0.042 microgram ch1 L(exp -1)). These results indicate the potential for forecasting chlorophyll concentration in this region but also highlights various deficiencies and suggestions for improvements to the current biogeochemical forecasting system. This system provides an initial basis for future applications including the effects of El Niño events on fisheries and other ocean resources given improvements identified in the analysis of these results.

Anomaly Correlation Coefficient (ACC)↗