Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ensemble Impacts”

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 361 records · Page 20

The Influence of Summertime Convection Over Southeast Asia on Water Vapor in the Tropical Stratosphere

The relative contributions of Southeast Asian convective source regions during boreal summer to water vapor in the tropical stratosphere are examined using Lagrangian trajectories. Convective sources are identified using global observations of infrared brightness temperature at high space and time resolution, and water vapor transport is simulated using advection-condensation. Trajectory simulations are driven by three different reanalysis data sets, GMAO MERRA, ERA-Interim, and NCEP/NCAR, to establish points of consistency and evaluate the sensitivity of the results to differences in the underlying meteorological fields. All ensembles indicate that Southeast Asia is a prominent boreal summer source of tropospheric air to the tropical stratosphere. Three convective source domains are identified within Southeast Asia: the Bay of Bengal and South Asian subcontinent (MON), the South China and Philippine Seas (SCS), and the Tibetan Plateau and South Slope of the Himalayas (TIB). Water vapor transport into the stratosphere from these three domains exhibits systematic differences that are related to differences in the bulk characteristics of transport. We find air emanating from SCS to be driest, from MON slightly moister, and from TIB moistest. Analysis of pathways shows that air detrained from convection over TIB is most likely to bypass the region of minimum absolute saturation mixing ratio over the equatorial western Pacific; however, the impact of this bypass mechanism on mean water vapor in the tropical stratosphere at 68 hPa is small 0.1 ppmv). This result contrasts with previously published hypotheses, and it highlights the challenge of properly quantifying fluxes of atmospheric humidity.

Wright, J. S.↗

Subseasonal-to-Seasonal Hindcast Skill Assessment of Ridging Events Related to Drought Over the Western United States

Persistent atmospheric ridging events centered near the western United States are associated with widespread precipitation deficits and meteorological drought. Due to the relatively low skill of- dynamical models in forecasting precipitation on subseasonal‐to‐seasonal (S2S) time scales across this region, forecasts of ridging are explored in this study as a potential bridge for early warning drought prediction. To assess skill, we evaluate deterministic and probabilistic S2S hindcasts (out to 6 week lead time) of ridging events and geopotential height anomalies over the western United States in five ensemble hindcast systems. Prediction skill for ridging events is shown to be highly variable across models. For some models, longer‐time‐averaged patterns of geopotential height anomalies across the first 6 weeks can be skillfully simulated, when evaluated probabilistically against climatology. The most skillful models show modest skill in forecasting above normal ridging occurrences at lead times of Weeks 3–4 and Weeks 5–6,with some sensitivity to the specific ridge location and method for determining skill. Using the European Centre for Medium‐Range Weather Forecasts (ECMWF) model as a case study, longer lead time forecast busts are shown to often occur under extended periods of highly active La Nina‐like tropical convection. Model errors at simulating these tropical convection features may have a disproportionately large impact and degrade downstream forecasts over the western United States. Despite the documented model shortcomings, our results highlight an opportunity to target skillful ridging forecasts from dynamical models for improving early warning drought forecasting at S2S lead times.

S2S↗

Soil Moisture and Snow Cover: Active or Passive Elements of Climate

A key question is the extent to which surface effects such as soil moisture and snow cover are simply passive elements or whether they can affect the evolution of climate on seasonal and longer time scales. We have constructed ensembles of predictability studies using the NCAR CCM3 in which we compared the relative roles of initial surface and atmospheric conditions over the central and western U.S. in determining the subsequent evolution of soil moisture and of snow cover. Results from simulations with realistic soil moisture anomalies indicate that internal climate variability may be the strongest factor, with some indication that the initial atmospheric state is also important. Model runs with exaggerated soil moisture reductions (near-desert conditions) showed a much larger effect, with warmer surface temperatures, reduced precipitation, and lower surface pressures; the latter indicating a response of the atmospheric circulation. These results suggest the possibility of a threshold effect in soil moisture, whereby an anomaly must be of a sufficient size before it can have a significant impact on the atmospheric circulation and climate. Results from simulations with realistic snow cover anomalies indicate that the time of year can be crucial. When introduced in late winter, these anomalies strongly affected the subsequent evolution of snow cover. When introduced in early winter, however, little or no effect is seen on the subsequent snow cover. Runs with greatly exaggerated initial snow cover indicate that the high reflectivity of snow is the most important process by which snow cover can impact climate, through lower surface temperatures and increased surface pressures. The results to date were obtained for model runs with present-day conditions. We are currently analyzing runs made with projected forcings for the 21st century to see if these results are modified in any way under likely scenarios of future climate change. An intriguing new statistical technique involving 'clustering' is developed to assist in this analysis.

Oglesby, Robert J.↗

Version 3 of the SMAP Level 4 Soil Moisture Product

The NASA Soil Moisture Active Passive (SMAP) Level 4 Soil Moisture (L4_SM) product provides 3-hourly, 9-km resolution, global estimates of surface (0-5 cm) and root zone (0-100 cm) soil moisture as well as related land surface states and fluxes from 31 March 2015 to present with a latency of 2.5 days. The ensemble-based L4_SM algorithm is a variant of the Goddard Earth Observing System version 5 (GEOS-5) land data assimilation system and ingests SMAP L-band (1.4 GHz) Level 1 brightness temperature observations into the Catchment land surface model. The soil moisture analysis is non-local (spatially distributed), performs downscaling from the 36-km resolution of the observations to that of the model, and respects the relative uncertainties of the modeled and observed brightness temperatures. Prior to assimilation, a climatological rescaling is applied to the assimilated brightness temperatures using a 6 year record of SMOS observations. A new feature in Version 3 of the L4_SM data product is the use of 2 years of SMAP observations for rescaling where SMOS observations are not available because of radio frequency interference, which expands the impact of SMAP observations on the L4_SM estimates into large regions of northern Africa and Asia. This presentation investigates the performance and data assimilation diagnostics of the Version 3 L4_SM data product. The L4_SM soil moisture estimates meet the 0.04 m3m3 (unbiased) RMSE requirement. We further demonstrate that there is little bias in the soil moisture analysis. Finally, we illustrate where the assimilation system overestimates or underestimates the actual errors in the system.

GEOS-↗

Strong Regional Influence of Climatic Forcing Datasets on Global Crop Model Ensembles

We present results from the Agricultural Model Intercomparison and Improvement Project (AgMIP) Global Gridded Crop Model Intercomparison (GGCMI) Phase I, which aligned 14 global gridded crop models (GGCMs) and 11 climatic forcing datasets (CFDs) in order to understand how the selection of climate data affects simulated historical crop productivity of maize, wheat, rice and soybean. Results show that CFDs demonstrate mean biases and differences in the probability of extreme events, with larger uncertainty around extreme precipitation and in regions where observational data for climate and crop systems are scarce. Countries where simulations correlate highly with reported FAO national production anomalies tend to have high correlations across most CFDs, whose influence we isolate using multi-GGCM ensembles for each CFD. Correlations compare favorably with the climate signal detected in other studies, although production in many countries is not primarily climate-limited (particularly for rice). Bias-adjusted CFDs most often were among the highest model-observation correlations, although all CFDs produced the highest correlation in at least one top-producing country. Analysis of larger multi-CFD-multi-GGCM ensembles (up to 91 members) shows benefits over the use of smaller subset of models in some regions and farming systems, although bigger is not always better. Our analysis suggests that global assessments should prioritize ensembles based on multiple crop models over multiple CFDs as long as a top-performing CFD is utilized for the focus region.

Agricultural Model Intercomparison and Improvement↗

Understanding and Utilizing PBL Height Data from Multiple Observing Systems in the GEOS System

The accuracy of PBL height simulation is a key issue in many applications including forecasting near surface meteorology and air quality, however, it is a very challenging problem due to the lack of not only comprehensive, global Planetary Boundary Layer (PBL) observations but also a strategy and infrastructure to utilize PBL height data from a variety of sensors. Following the designation of PBL as an incubation class observable in the 2017 Decadal Survey, the PBL Incubation Study Team Report [14] made clear that “a future global PBL observing system requires modeling and data assimilation as essential components.” There is an urgent need for global modeling development in order to utilize Program of Record (POR) observations, assess their impacts, and identify gaps to be filled by future PBL missions. Our overall objective is to develop PBL data assimilation capabilities in the NASA Global Earth Observing System (GEOS), focusing on PBL height from multiple observing systems, to support the assessment and use of future PBL observations. The NASA GEOS system is composed of the GEOS global atmospheric general circulation model (AGCM) and the atmospheric data assimilation system (ADAS). The PBL parameterizations include the “Lock” K-profile scheme driven by surface and cloud-top buoyancy fluxes ([4]), and the “Louis” local scheme for stable conditions based on the Richardson number ([5]). Above the mixed layer defined by the Lock surface plume, shallow cumulus convection is represented by the mass flux scheme of [9]. Additional parameterizations are summarized in [1]. The ADAS employs the hybrid 4D Ensemble- Variational (EnVar) configuration ([15]), with the ensemble providing flow-dependent background error covariance information. The resultant analysis increments are fed back to the forecast model through the 4D incremental analysis update (IAU) approach ([11]). In this study, PBL height data are being or have been generated from radiosondes, GNSS RO, satellite (CATS, CALIPSO and ICESat-2) and ground-based (MPLNET) lidars, and wind profiler. Investigations have been conducted to specify quality marks for PBL height retrievals for the data assimilation purpose. These PBL height data have different strengths and weaknesses ([2], [3], [6], [7], [8], [10]), and the satellite PBL height data provide better global coverage and complement in-situ PBL height data. Radiosondes offer high accuracy and in situ measurement of temperature and humidity profiles, but with poor spatio-temporal sampling. The in-situ observing systems like MPLNET and wind profiler provide long history of PBL height records at each station. The GNSS RO based PBL height is retrieved based on the sharp gradients in refractivity profile that represent the fine vertical structure of temperature and moisture changes above the PBL. However, not all RO refractivity profiles reach the surface depending on location and regime, and RO refractivity retrievals can be negatively biased below 2km. The PBL height data from satellite lidars provide high resolution along track PBL height retrievals, but over land they are affected by previous day convective PBL aerosol and strongly associated with mixing layer and retrievals cannot be made below thick, attenuating clouds. A successful assimilation of PBL height data requires a thorough understanding of the observing method and the retrieval algorithm for each observing system in order to use the PBL height data from multiple observing systems properly. Due to the sensitivity of PBL height data to the observing method and choice of algorithm, it is important to use a model definition appropriate for each observation type to compute differences between PBL height data and model PBL height (OmFs). The GEOS model currently includes two PBL height definitions suitable for direct comparison with observed PBL height, and additional definitions are being added in this study. Evaluation of different model PBL height definitions is underway. Meanwhile, efforts have been made in the GEOS data assimilation system to develop PBL height data assimilation capability. PBL height data can be assimilated using two different approaches. The traditional approach is to construct an observation operator and its tangent linear and adjoint, which link control variables to PBL height data from each observing system. This observation operator can be very complicated, e.g., the lidar-based PBL height observation operator includes the backscatter lidar forward observation operator, the algorithm to derive PBL height from attenuated total backscatter, interpolation, and calculations handling the mismatch between observed and model scales. The other approach is to augment PBL height to the control variable vector, and it is adopted in this study. The latter approach was also used in previous studies, e.g., the assimilation of PBL height data from radiosonde and aircraft in the Real Time Mesoscale Analysis (RTMA) system for a dispersion modelling study ([13]); the PBL height assimilation study using lidar PBL height data at Greensburg, Kansas for a field campaign ([12]). The PBL height assimilation from multiple observing systems in this study allows us to take advantage of the diverse PBL height data that provide much better global coverage collectively under different meteorological conditions and with different temporal and spatial scales. As all the PBL heights are tightly coupled with the PBL thermodynamic variables, the strong correlations, which are provided by the 4D ensemble forecast, enable PBL height data from various sources to interact and combine coherently and provide additional information for PBL temperature and moisture fields. The results of comparisons among PBL height data from different sources and the evaluation of the model PBL height definitions with the PBL height data will be presented, and the PBL height data synergy strategies and preliminary results will also be discussed at the conference.

Y. Zhu↗

Acoustic Testing of Flight Hardware Using Loudspeakers: How Much do We Know About This Method of Testing?

Loudspeakers have been used for acoustic qualification of spacecrafts, reflectors, solar panels, and other acoustically responsive structures for more than a decade. Even though a lot of hardware has been acoustic tested using this method, the nature of the acoustic field generated by controlling an ensemble of speakers with and without the hardware in the test volume has not been thoroughly investigated. Limited measurements from some of the recent speaker tests used to qualify flight hardware have indicated significant spatial variation of the acoustic field within the test volume. Also structural responses have been reported to differ when similar tests were performed using reverberant chambers. Unlike the reverberant chamber acoustic test, for which the acoustic field in most chambers is known to be diffuse except below several tens of Hz where acoustic standing waves and large spatial variations exist, the characteristics of the acoustic field within the speaker test volume has not been quantified. It has only been recently that a detailed acoustic field characterization of speaker testing has been made at Jet Propulsion Laboratory (JPL) with involvement of various organizations. To address the impact of non-uniform acoustic field on structures, a series of acoustic tests were performed using a flat panel and a 3-ft cylinder exposed to the field controlled by speakers and repeated in a reverberant chamber. The analysis of the data from this exercise reveals that there are significant differences both in the acoustic field and in the structural responses. In this paper the differences between the two methods are reviewed in some detail and the over- or under-testing of articles that could pose un-anticipated structural and flight qualification issues are discussed. A framework for discussing the validity of the speaker acoustic testing method with the current control system and a path forward for improving it will be provided.

Acoustics↗

Satellite-Based Evaluation of AeroCom Model Bias in Biomass Burning Regions

Global models are widely used to simulate biomass burning aerosol (BBA). Exhaustive evaluations on model representation of aerosol distributions and properties are fundamental to assess health and climate impacts of BBA. Here we conducted a comprehensive comparison of Aerosol Comparisons between Observations and Models (AeroCom) project model simulations with satellite observations. A total of 59 runs by 18 models from three AeroCom Phase-III experiments (i.e., biomass burning emissions, CTRL16, and CTRL19) and 14 satellite products of aerosols were used in the study. Aerosol optical depth (AOD) at 550 nm was investigated during the fire season over three key fire regions reflecting different fire dynamics (i.e., deforestation-dominated Amazon, Southern Hemisphere Africa where savannas are the key source of emissions, and boreal forest burning in boreal North America). The 14 satellite products were first evaluated against AErosol RObotic NETwork (AERONET) observations, with large uncertainties found. But these uncertainties had small impacts on the model evaluation that was dominated by modeling bias. Through a comparison with Polarization and Directionality of the Earth’s Reflectances measurements with the Generalized Retrieval of Aerosol and Surface Properties algorithm (POLDER-GRASP), we found that the modeled AOD values were biased by −93 % to 152 %, with most models showing significant underestimations even for the state-of-the-art aerosol modeling techniques (i.e., CTRL19). By scaling up BBA emissions, the negative biases in modeled AOD were significantly mitigated, although it yielded only negligible improvements in the correlation between models and observations, and the spatial and temporal variations in AOD biases did not change much. For models in CTRL16 and CTRL19, the large diversity in modeled AOD was in almost equal measures caused by diversity in emissions, lifetime, and the mass extinction coefficient (MEC). We found that in the AeroCom ensemble, BBA lifetime correlated significantly with particle deposition (as expected) and in turn correlated strongly with precipitation. Additional analysis based on Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) aerosol profiles suggested that the altitude of the aerosol layer in the current models was generally too low, which also contributed to the bias in modeled lifetime. Modeled MECs exhibited significant correlations with the Ångström exponent (AE, an indicator of particle size). Comparisons with the POLDER-GRASP-observed AE suggested that the models tended to overestimate the AE (underestimated particle size), indicating a possible underestimation of MECs in models. The hygroscopic growth in most models generally agreed with observations and might not explain the overall underestimation of modeled AOD. Our results imply that current global models contain biases in important aerosol processes for BBA (e.g., emissions, removal, and optical properties) that remain to be addressed in future research.

biomass burning aerosols↗

Dust Model Sensitivity to Dust Source Mask, Sandblasting Efficiency, Air Density, and Land Use: Implications for Model Improvement

This study compares dust storm simulations using two commonly adopted methods for representing four important dust emission parameters. Compared to a dynamic dust source mask based on land use and vegetation cover, a static mask based solely on land use overestimates dust concentration and optical depth by a factor of 2, besides generating spurious emissions. The results reinforce that seasonal variations in vegetation cover can significantly affect dust emissions. For sandblasting efficiency, a clay-dependent semiempirical expression produces 10 times more dust than a physics-based expression. Simulations using model-predicted versus a fixed constant for air density differ by only 8%. However, this difference could range between 12 and 22% for annual simulations over global dust source regions. Simulations with updated versus old land use data, using the same dust source mask, differ twofold, indicating the significant impact of land use change on regional dust emission in central Arizona. The differences in the pairs of these simulations are generally larger than the uncertainty due to meteorology. The simulations align better with observation when using the dynamic dust source mask, the physics-based sandblasting efficiency, and the up-to-date land use data. Given the high sensitivity of dust to surface conditions, the results discussed have implications for improving the dust cycle in weather and climate models and for interpreting model intercomparisons.

Dust-emission source↗

The Influence of Soil Moisture, Coastline Curvature, and Land-Breeze Circulations on Sea-Breeze Initiated Precipitation

Idealized numerical simulations are performed with a coupled atmosphere/land-surface model to identify the roles of initial soil moisture, coastline curvature, and land breeze circulations on sea breeze initiated precipitation. Data collected on 27 July 1991 during the Convection and Precipitation Electrification Experiment (CAPE) in central Florida are used. The 3D Goddard Cumulus Ensemble (GCE) cloud resolving model is coupled with the Goddard Parameterization for Land-Atmosphere-Cloud Exchange (PLACE) land surface model, thus providing a tool to simulate more realistically land-surface/atmosphere interaction and convective initiation. Eight simulations are conducted with either straight or curved coast-lines, initially homogeneous soil moisture or initially variable soil moisture, and initially homogeneous horizontal winds or initially variable horizontal winds (land breezes). All model simulations capture the diurnal evolution and general distribution of sea-breeze initiated precipitation over central Florida. The distribution of initial soil moisture influences the timing, intensity and location of subsequent precipitation. Soil moisture acts as a moisture source for the atmosphere, increases the connectively available potential energy, and thus preferentially focuses heavy precipitation over existing wet soil. Strong soil moisture-induced mesoscale circulations are not evident in these simulations. Coastline curvature has a major impact on the timing and location of precipitation. Earlier low-level convergence occurs inland of convex coastlines, and subsequent precipitation occurs earlier in simulations with curved coastlines. The presence of initial land breezes alone has little impact on subsequent precipitation. however, simulations with both coastline curvature and initial land breezes produce significantly larger peak rain rates due to nonlinear interactions.

Baker, David R.↗

Uncertainty Models for the Hybrid Parametric Variation Method of Uncertainty Quantification; Analysis

There is some level of uncertainty in every finite element model (FEM), which flows to a level of uncertainty in predicted results. The purpose of uncertainty quantification (UQ) is to provide statistical bounds on prediction accuracy based on model uncertainty. This is distinct from model updating, which attempts to modify models to improve their accuracy. UQ does not improve the accuracy of models, but accepts that the models are inherently inaccurate and attempts to quantify the impact of that inaccuracy on predicted results. Previously, an alternate method for UQ, called the Hybrid Parametric Variation (HPV) method, was applied to Space Launch System (SLS) Hurty/Craig-Bampton (HCB) components to predict system-level statistics for launch vehicle attitude control transfer functions and core stage section loads due to buffet. The HPV method combines a parametric variation of the HCB fixed-interface (FI) modal frequencies with a nonparametric variation (NPV) method that randomly varies the HCB mass and stiffness matrices as Wishart random matrix distributions using random matrix theory (RMT). Alternatively, the most common method for modeling uncertainty in the structural dynamics community is a parametric approach, which varies physical parameters in the model. However, there are several disadvantages associated with the parametric method. Determining a reduced set of parameters that have a significant impact on the system response can be time consuming, and the selected parameter probability distributions are rarely reliably known. Therefore, in practice, the parameters are surrogates for the actual errors, and the link to parameter uncertainty is unknown. Another major drawback is that the uncertainty that can be represented is limited to the form of the nominal FEM. It is the experience of the authors that based on numerous aerospace programs, almost all FEM errors are in form rather than parameter values. This hypothesis is supported by the observation of the authors that it is almost never possible to ‘tune’ a FEM to match modal test results by only modifying model parameters. Model-form uncertainty cannot be directly represented by FEM input parameters nor included in a parametric approach. However, model-form uncertainty can be modeled using RMT, where a probability distribution is developed for the matrix ensemble of interest. The major advantage of the NPV method is that it covers errors in model form. The HPV method anchors uncertainty at the HCB component level to component modal test results by matching the HCB and test modes based on mode descriptions or other methods, and then applying differing levels of frequency variation. The specific variations depend on the confidence to which a component FEM has been validated through modal testing. The NPV method is layered on the frequency variation to match modal test self-orthogonality and cross-orthogonality (XO) results. Once the component uncertainty models are identified, they are assembled, and the uncertainty is propagated to the system level using a Monte Carlo (MC) analysis approach that generates statistics for system-level predictions This provides a UQ method that can be traced to test data, which can be updated as additional data and improved correlated models become available. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The second example is an application to a component that is design specific to the SLS. Based on this work and other assessments, the HPV method provides another tool to the toolset used for complex system UQ analysis. From experience gathered to date using the HPV method, additional design specific applications must be investigated to provide further confidence in the validity of the HPV method of UQ analysis.

Uncertainty quantification↗

Understanding Climate Policy Data Needs. NASA Carbon Monitoring System Briefing: Characterizing Flux Uncertainty, Washington D.C., 11 January 2012

Climate policy in the United States is currently guided by public-private partnerships and actions at the local and state levels. This mitigation strategy is made up of programs that focus on energy efficiency, renewable energy, agricultural practices and implementation of technologies to reduce greenhouse gases. How will policy makers know if these strategies are working, particularly at the scales at which they are being implemented? The NASA Carbon Monitoring System (CMS) will provide information on carbon dioxide fluxes derived from observations of earth's land, ocean and atmosphere used in state of the art models describing their interactions. This new modeling system could be used to assess the impact of specific policy interventions on CO2 reductions, enabling an iterative, results-oriented policy process. In January of 2012, the CMS team held a meeting with carbon policy and decision makers in Washington DC to describe the developing modeling system to policy makers. The NASA CMS will develop pilot studies to provide information across a range of spatial scales, consider carbon storage in biomass, and improve measures of the atmospheric distribution of carbon dioxide. The pilot involves multiple institutions (four NASA centers as well as several universities) and over 20 scientists in its work. This pilot study will generate CO2 flux maps for two years using observational constraints in NASA's state-of -the-art models. Bottom-up surface flux estimates will be computed using data-constrained land and ocean models; comparison of the different techniques will provide some knowledge of uncertainty in these estimates. Ensembles of atmospheric carbon distributions will be computed using an atmospheric general circulation model (GEOS-5), with perturbations to the surface fluxes and to transport. Top-down flux estimates will be computed from observed atmospheric CO2 distributions (ACOS/GOSAT retrievals) alongside the forward-model fields, in conjunction with an inverse approach based on the CO2 model of GEOS ]Chem. The forward model ensembles will be used to build understanding of relationships among surface flux perturbations, transport uncertainty and atmospheric carbon concentration. This will help construct uncertainty estimates and information on the true spatial resolution of the top-down flux calculations. The relationship between the top-down and bottom-up flux distributions will be documented. Because the goal of NASA CMS is to be policy relevant, the scientists involved in the flux modeling pilot need to understand and be focused on the needs of the climate policy and decision making community. If policy makers are to use CMS products, they must be aware of the modeling effort and begin to design policies that can be evaluated with information. Improving estimates of carbon sequestered in forests, for example, will require information on the spatial variability of forest biomass that is far more explicit than is presently possible using only ground observations. Carbon mitigation policies being implemented by cities around the United States could be designed with the CMS data in mind, enabling sequential evaluation and subsequent improvements in incentives, structures and programs. The success of climate mitigation programs being implemented in the United States today will hang on the depth of the relationship between scientists and their policy and decision making counterparts. Ensuring that there is two-way communication between data providers and users is important for the success both of the policies and the scientific products meant to support them..

Brown, Molly E.↗

An Overview of the GEOS-5 Aerosol Reanalysis

GEOS-5 is the latest version of the NASA Global Modeling and Assimilation Office (GMAO) earth system model. GEOS-5 contains components for atmospheric circulation and composition (including data assimilation), ocean circulation and biogeochemistry, and land surface processes. In addition to traditional meteorological parameters, GEOS-5 includes modules representing the atmospheric composition, most notably aerosols and tropospheric/stratospheric chemical constituents, taking explicit account of the impact of these constituents on the radiative processes of the atmosphere. MERRA is a NASA meteorological reanalysis for the satellite era (1979-present) using GEOS-5. This project focuses on historical analyses of the hydrological cycle on a broad range of weather and climate time scales. As a first step towards an integrated Earth System Analysis (IESA), the GMAO is extending MERRA with reanalyses for other components of the earth system: land, ocean, bio-geochemistry and atmospheric constituents. In this talk we will present results from the MERRA-driven aerosol reanalysis covering the Aqua period (2003-present). The assimilation of Aerosol Optical Depth (AOD) in GEOS-5 involves very careful cloud screening and homogenization of the observing system by means of a Neural Net scheme that translates MODIS radiances into AERONET calibrated AOD. These measurements are further quality controlled using an adaptive buddy check scheme, and assimilated using the Local Displacement Ensemble (LDE) methodology. For this reanalysis, GEOS-5 runs at a nominal 50km horizontal resolution with 72 vertical layers (top at approx. 8Skm). GEOS-5 is driven by daily biomass burning emissions derived from MODIS fire radiative power retrievals. We will present a summary of our efforts to validate such dataset. The GEOS-5 assimilated aerosol fields are first validated by comparison to independent in-situ measurements (AERONET and PM2.5 surface concentrations). In order to asses aerosol absorption on a global scale, we perform a detailed radiative transfer calculation to simulate the UV aerosol index, comparing our results to OMI measurements. By simulating aerosol attenuated backscatter, we use CALIPSO measurements to evaluate the vertical structure of our aerosol estimates, in particular in regions where we have larger discrepancies with OMI. Finally, the consistency of our AOD estimates with estimates from MISR, MODIS/Deep Blue, OMI and PARASOL will be briefly discussed.

da Silva, Arlindo↗

Volcanic Eruption Signatures in the Isotope Enabled Last Millennium Ensemble

Explosive volcanic eruptions are one of the largest natural climate perturbations, but few observational constraints exist on either the climate responses to eruptions or the properties (size, hemispheric aerosol distribution, etc.) of the eruptions themselves. Paleoclimate records are thus important sources of information on past eruptions, often through the measurement of oxygen isotopic ratios (𝛿18O) in natural archives. However, since many processes affect 𝛿18O, the dynamical interpretation of these records can be quite complex. Here we present results from new, isotope-enabled members of the Community Earth System Model Last Millennium Ensemble, documenting eruption-induced 𝛿18O variations throughout the climate system. Eruptions create significant perturbations in the 𝛿18O of precipitation and soil moisture in central/eastern North America, via excitation of the Atlantic Multidecadal Oscillation. Monsoon Asia and Australia also exhibit strong precipitation and soil 𝛿18O anomalies; in these cases, 𝛿18O may reflect changes to El Niño-Southern Oscillation phase following eruptions. Salinity and seawater 𝛿18O patterns demonstrate the importance of both local hydrologic shifts and the phasing of the El Niño-Southern Oscillation response, both along the equator and in the subtropics. In all cases, the responses are highly sensitive to eruption latitude, which points to the utility of isotopic records in constraining aerosol distribution patterns associated with past eruptions. This is most effective using precipitation 𝛿18O; all Southern eruptions and the majority (66%) of Northern eruptions can be correctly identified. This work thus serves as a starting point for new, quantitative uses of isotopic records for understanding volcanic impacts on climate.

Stevenson, S.↗

CMEs and SEPs During November–December 2020: A Challenge for Real-Time Space Weather Forecasting

Predictions of coronal mass ejections (CMEs) and solar energetic particles (SEPs) are a central issue in space weather forecasting. In recent years, interest in space weather predictions has expanded to include impacts at other planets beyond Earth as well as spacecraft scattered throughout the heliosphere. In this sense, the scope of space weather science now encompasses the whole heliospheric system, and multipoint measurements of solar transients can provide useful insights and validations for prediction models. In this work, we aim to analyze the whole inner heliospheric context between two eruptive flares that took place in late 2020, that is, the M4.4 flare of 29 November and the C7.4 flare of 7 December. This period is especially interesting because the STEREO-A spacecraft was located ∼60° east of the Sun–Earth line, giving us the opportunity to test the capabilities of “predictions at 360°” using remote-sensing observations from the Lagrange L1 and L5 points as input. We simulate the CMEs that were ejected during our period of interest and the SEPs accelerated by their shocks using the WSA–Enlil–SEPMOD modeling chain and four sets of input parameters, forming a “mini-ensemble.” We validate our results using in situ observations at six locations, including Earth and Mars. We find that, despite some limitations arising from the models' architecture and assumptions, CMEs and shock-accelerated SEPs can be reasonably studied and forecast in real time at least out to several tens of degrees away from the eruption site using the prediction tools employed here.

Coronal mass ejections↗

The Effects of Lightning NO(x) Production during the July 21 EULINOX Storm studied with a 3-D Cloud-scale Chemical Transport Model

The July 21,1998 thunderstonn observed during the European Lightning Nitrogen Oxides Project (EULINOX) project was simulated using the three-dimensional Goddard Cumulus Ensemble (GCE) model. The simulation successfully reproduced a number of observed storm features including the splitting of the original cell into a southern cell which developed supercell characteristics, and a northern cell which became multicellular. Output from the GCE simulation was used to drive an offline cloud-scale chemical transport model which calculates tracer transport and includes a parameterization of lightning NO(x) production which uses observed flash rates as input. Estimates of lightning NO(x) production were deduced by assuming various values of production per intracloud and production per cloud-to-ground flash and comparing the results with in-cloud aircraft observations. The assumption that both types of flashes produce 360 moles of NO per flash on average compared most favorably with column mass and probability distribution functions calculated from observations. This assumed production per flash corresponds to a global annual lightning NOx source of 7 Tg N per yr. Chemical reactions were included in the model to evaluate the impact of lightning NO(x), on ozone. During the storm, the inclusion of lightning NOx in the model results in a small loss of ozone (on average less than 4 ppbv) at all model levels. Simulations of the chemical environment in the 24 hours following the storm show on average a small increase in the net production of ozone at most levels resulting from lightning NO(x), maximizing at approximately 5 ppbv per day at 5.5 km. Between 8 and 10.5 km, lightning NO(x) causes decreased net ozone production.

Ott, Lesley E.↗

On the Role of Trend and Variability of Hydroxyl Radical (OH) in the Global Methane Budget

Decadal trends and interannual variations in the hydroxyl radical (OH), while poorly constrained at present, are critical for understanding the observed evolution of atmospheric methane (CH4). Through analyzing the OH fields simulated by the model ensemble of the Chemistry-Climate Model Initiative (CCMI), we find (1) the negative OH anomalies during the El Niño years mainly corresponding to the enhanced carbon monoxide (CO) emissions from biomass burning and (2) a positive OH trend during 1980-2010 dominated by the elevated primary production and the reduced loss of OH due to decreasing CO after 2000. Both two-box model inversions and variational 4D inversions suggest that ignoring the negative anomaly of OH during the El Niño years leads to a large overestimation of the increase in global CH4 emissions by up to10±3Tg yr-1to match the observed CH4 increase over these years. Not accounting for the increasing OH trends given by the CCMI models leads to an underestimation of the CH4 emission increase by 23±9Tg yr-1from 1986 to 2010. The variational inversion estimated CH4 emissions show that the tropical regions contribute most to the uncertainties related to OH. This study highlights the significant impact of climate and chemical feedbacks related to OH on the top-down estimates of the global CH4 budget.

Hydroxyl radical (OH)↗

Airborne Observations of the Spatial and Temporal Variability of Tropospheric Carbon Dioxide during the INTEX-B Campaign

The Intercontinental Chemical Transport Experiment-North America (INTEX-NA) is an international field campaign envisioned to investigate the transport and transformation of gases and aerosols on transcontinental/intercontinental scales and assess their impact on air quality and climate. Phase B (INTEX-B) of the mission was conducted during a 10- week period from March 1 to May 15, 2006 and focused initially on pollution outflow from the Mexico City Metropolitan Area, later addressing the transport of pollution from Asia to North America during springtime meteorological conditions. During the deployment, fast-response (1-s resolution) CO2 measurements were recorded aboard the NASA DC-8 providing valuable regional-scale information on carbon sources and sinks over sparsely sampled areas of North America and adjacent ocean basins. When coupled with the enormously sophisticated chemistry payload on the DC-8, these measurements collectively afford extremely powerful multi-tracer constraints for carbon source/sink attribution. Preliminary examination of the two data sets from the INTEX-B campaign, acquired one month apart, reveals not only the influence of the CO2 seasonal cycle, but also the preponderance of human population and industrial activity in the northern hemisphere. In this presentation, a synergy of the ensemble of airborne and surface observations, bottomup emission inventories, as well as transport history are invoked in a GIS framework to elucidate the source/sink processes reflected in the observations. The airborne CO2 data, along with simultaneous surface measurements (e.g. NOAA ESRL), are examined to establish the vertical distribution and variability of CO2 as a function of location. The role of localized sources, long-range transport, the biosphere, stratospheric exchange, and dynamical processes on the CO2 spatial variability observed throughout the tropospheric column will be discussed.

Vay, Stephanie A.↗