SEARCH · Engineering Papers

Results for “ensemble simulations”

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

SourceResearch connectionFrom yearThrough year

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

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Development and Evaluation of Ensemble Learning-based Environmental Methane Detection and Intensity Prediction Models

The environmental impacts of global warming driven by methane (CH 4 ) emissions have catalyzed significant research initiatives in developing novel technologies that enable proactive and rapid detection of CH 4 . Several data-driven machine learning (ML) models were tested to determine how well they identified fugitive CH 4 and its related intensity in the affected areas. Various meteorological characteristics, including wind speed, temperature, pressure, relative humidity, water vapor, and heat flux, were included in the simulation. We used the ensemble learning method to determine the best-performing weighted ensemble ML models built upon several weaker lower-layer ML models to (i) detect the presence of CH 4 as a classification problem and (ii) predict the intensity of CH 4 as a regression problem. The classification model performance for CH 4 detection was evaluated using accuracy, F1 score, Matthew’s Correlation Coefficient (MCC), and the area under the receiver operating characteristic curve (AUC ROC), with the top-performing model being 97.2%, 0.972, 0.945 and 0.995, respectively. The R 2 score was used to evaluate the regression model performance for CH 4 intensity prediction, with the R 2 score of the best-performing model being 0.858. The ML models developed in this study for fugitive CH 4 detection and intensity prediction can be used with fixed environmental sensors deployed on the ground or with sensors mounted on unmanned aerial vehicles (UAVs) for mobile detection.

Majumder, Reek↗

NASA NTRS2014

Soil Moisture Data Assimilation in the NASA Land Information System for Local Modeling Applications and Improved Situational Awareness

As part of the NASA Soil Moisture Active Passive (SMAP) Early Adopter (EA) program, the NASA Shortterm Prediction Research and Transition (SPoRT) Center has implemented a data assimilation (DA) routine into the NASA Land Information System (LIS) for soil moisture retrievals from the European Space Agency's Soil Moisture Ocean Salinity (SMOS) satellite. The SMAP EA program promotes application‐driven research to provide a fundamental understanding of how SMAP data products will be used to improve decision‐making at operational agencies. SPoRT has partnered with select NOAA/NWS Weather Forecast Offices (WFOs) that use output from a real‐time regional configuration of LIS, without soil moisture DA, to initialize local numerical weather prediction (NWP) models and enhance situational awareness. Improvements to local NWP with the current LIS have been demonstrated; however, a better representation of the land surface through assimilation of SMOS (and eventually SMAP) retrievals is expected to lead to further model improvement, particularly during warm‐season months. SPoRT will collaborate with select WFOs to assess the impact of soil moisture DA on operational forecast situations. Assimilation of the legacy SMOS instrument data provides an opportunity to develop expertise in preparation for using SMAP data products shortly after the scheduled launch on 5 November 2014. SMOS contains a passive L‐band radiometer that is used to retrieve surface soil moisture at 35‐km resolution with an accuracy of 0.04 cu cm cm (exp -3). SMAP will feature a comparable passive L‐band instrument in conjunction with a 3‐km resolution active radar component of slightly degraded accuracy. A combined radar‐radiometer product will offer unprecedented global coverage of soil moisture at high spatial resolution (9 km) for hydrometeorological applications, balancing the resolution and accuracy of the active and passive instruments, respectively. The LIS software framework manages land surface model (LSM) simulations and includes an Ensemble Kalman Filter for conducting land surface DA. SPoRT has added a module to read, quality‐control and bias‐correct swaths of Level II SMOS soil moisture retrievals prior to assimilation within LIS. The impact of SMOS DA is being tested using the Noah LSM. Experiments are being conducted to examine the impacts of SMOS soil moisture DA on the resulting LISNoah fields and subsequent NWP simulations using the Weather Research and Forecasting (WRF) model initialized with LIS‐Noah output. LIS‐Noah soil moisture will be validated against in situ observations from Texas A&M's North American Soil Moisture Database to reveal the impact and possible improvement in soil moisture trends through DA. WRF model NWP case studies will test the impacts of DA on the simulated near‐surface and boundary‐layer environments, and precipitation during both quiescent and disturbed weather scenarios. Emphasis will be placed on cases with large analysis increments, especially due to contributions from regional irrigation patterns that are not represented by precipitation input in the baseline LIS‐Noah run. This poster presentation will describe the soil moisture DA methodology and highlight LIS‐Noah and WRF simulation results with and without assimilation.

Case, Jonathan L.↗

DOE OSTI2020

Robust late twenty-first century shift in the regional monsoons in RegCM-CORDEX simulations

Abstract We use an unprecedented ensemble of regional climate model (RCM) projections over seven regional CORDEX domains to provide, for the first time, an RCM-based global view of monsoon changes at various levels of increased greenhouse gas (GHG) forcing. All regional simulations are conducted using RegCM4 at a 25 km horizontal grid spacing using lateral and lower boundary forcing from three General Circulation Models (GCMs), which are part of the fifth phase of the Coupled Model Inter-comparison Project (CMIP5). Each simulation covers the period from 1970 through 2100 under two Representative Concentration Pathways (RCP2.6 and RCP8.5). Regional climate simulations exhibit high fidelity in capturing key characteristics of precipitation and atmospheric dynamics across monsoon regions in the historical period. In the future period, regional monsoons exhibit a spatially robust delay in the monsoon onset, an increase in seasonality, and a reduction in the rainy season length at higher levels of radiative forcing. All regions with substantial delays in the monsoon onset exhibit a decrease in pre-monsoon precipitation, indicating a strong connection between pre-monsoon drying and a shift in the monsoon onset. The weakening of latent heat driven atmospheric warming during the pre-monsoon period delays the overturning of atmospheric subsidence in the monsoon regions, which defers their transitioning into deep convective states. Monsoon changes under the RCP2.6 scenario are mostly within the baseline variability.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Nonane and Hexanol Adsorption in the Lamellar Phase of a Nonionic Surfactant: Molecular Simulations and Comparison to Ideal Adsorbed Solution Theory

Adsorption of n-nonane/1-hexanol (C9/C6OH) mixtures into the lamellar phase formed by a 50/50 w/w triethylene glycol mono-n-decyl ether (C10E3)/water system was studied using configurational-bias Monte Carlo simulations in the osmotic Gibbs ensemble. The interactions were described by the Shinoda–Devane–Klein coarse-grained force field. Prior simulations probing single-component adsorption indicated that C9 molecules preferentially load near the center of the bilayer increasing the bilayer thickness, whereas C6OH molecules are more likely to be found near the interface of the polar and non-polar moieties swelling the bilayer in the lateral dimension. Here, we extend this work to binary C9/C6OH adsorption to probe whether the difference in the spatial preferences may lead to a synergistic effect and enhanced loadings for the mixture. Comparing loading trends and the thermodynamics of binary adsorption to unary adsorption reveals that C9–C9 interactions lead to the largest enhancement, whereas C9–C6OH and C6OH–C6OH interactions are less favorable for this bilayer system. As a result, ideal adsorbed solution theory yields satisfactory predictions of the binary loading.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

A Deep Learning-Driven Sampling Technique to Explore the Phase Space of an RNA Stem-Loop

The folding and unfolding of RNA stem-loops are critical biological processes; however, their computational studies are often hampered by the ruggedness of their folding landscape, necessitating long simulation times at the atomistic scale. Here, we adapted DeepDriveMD (DDMD), an advanced deep learning-driven sampling technique originally developed for protein folding, to address the challenges of RNA stem-loop folding. Although tempering- and order parameter-based techniques are commonly used for similar rare-event problems, the computational costs or the need for a priori knowledge about the system often present a challenge in their effective use. DDMD overcomes these challenges by adaptively learning from an ensemble of running MD simulations using generic contact maps as the raw input. DeepDriveMD enables on-the-fly learning of a low-dimensional latent representation and guides the simulation toward the undersampled regions while optimizing the resources to explore the relevant parts of the phase space. We showed that DDMD estimates the free energy landscape of the RNA stem-loop reasonably well at room temperature. Our simulation framework runs at a constant temperature without external biasing potential, hence preserving the information on transition rates, with a computational cost much lower than that of the simulations performed with external biasing potentials. Here, we also introduced a reweighting strategy for obtaining unbiased free energy surfaces and presented a qualitative analysis of the latent space. This analysis showed that the latent space captures the relevant slow degrees of freedom for the RNA folding problem of interest. Finally, throughout the manuscript, we outlined how different parameters are selected and optimized to adapt DDMD for this system. We believe this compendium of decision-making processes will help new users adapt this technique for the rare-event sampling problems of their interest.

Gupta, Ayush↗

NASA NTRS2003

Simulated Radar Characteristics of LBA Convective Systems: Easterly and Westerly Regimes

The 3D Goddard Cumulus Ensemble (GCE) model was used to simulate convection that occurred during the TRMM LBA field experiment in Brazil. Convection in this region can be categorized into two different regimes. Low-level easterly flow results in moderate to high CAPE and a drier environment. Convection is more intense like that seen over continents. Low-level westerly flow results in low CAPE and a moist environment. Convection is weaker and more widespread characteristic of oceanic or monsoon-like systems. The GCE model has been used to study both regimes n order to provide cloud datasets that are representative of both environments in support of TRMM rainfall and heating algorithm development. Two different cases are analyzed: Jan 26, 1999, an eastely regime case, and Feb 23, 1999, a westerly regime case. The Jan 26 case is an organized squall line, while the Feb 23 case is less organized with only transient lines. Radar signatures, including CFADs, from the two simulated cases are compared to each other and with observations. The microphysical processes simulated in the model are also compared between the two cases.

Lang, Stephen E.↗

DOE OSTI2021

Unraveling the implications of finite specimen size on the interpretation of dynamic experiments for polycrystalline aluminum through direct numerical simulations

Normal and Pressure-shear plate impact (NPI and PSPI) experiments are popular experimental techniques for studying the mean-field macroscopic behavior of polycrystalline metals under high-rate dynamic loading. However, since both configurations rely upon geometry for subjecting the specimen to high strain rates, these experiments often involve a limited specimen size. Moreover, because of the inherent heterogeneities present within polycrystalline metals, it is difficult to ascertain if the size of the specimen and/or regions where measurements are made are sufficiently large for making representative inferences about the mean-field macroscopic properties from single-point velocity measurements. In the present study, we quantify the expected measurement variability on observable point measurements in NPI and PSPI experiments by carrying out direct numerical simulations (DNS) of statistically representative polycrystalline microstructures subjected to dynamic compression and compression-shear loading. In particular, we consider the role of specific material heterogeneities (e.g. the grain-to-grain difference in size, crystallographic orientation) on dispersion in the normal and transverse particle velocity records and on local fluctuations in key state variables (e.g. velocity, accumulated plastic strain) by incorporating these effects directly into a representative synthetic microstructure geometry and crystalline description of pure polycrystalline aluminum. The form of the present study is a large parametric investigation, consisting of ten ensembles of one hundred simulations. Each of the thousand simulations reflects a randomly realized synthetic microstructure in one of five cases of decreasing average grain size for the two loading configurations. Our analysis of the DNS results demonstrates that for both of these experimental configurations, the grain size directly correlates with the coefficient of variation (CV) in simulated point measurements, showing a convergent decrease in CV to zero (i.e. particle velocity record approaches the mean-field value) with decreasing grain size. Remarkably, the magnitude of variations in the particle velocity record is shown to be largest where the deviatoric stresses are most significant. In the case of NPI, this occurs at the elastic and plastic wavefront, whereas, in the case of PSPI, the magnitude of fluctuations are approximately constant throughout the experimental window time. The reasoning for the scatter in particle velocity due to the heterogeneous microstructure is demonstrated to be dependent on the mechanisms for accommodating deformation and on the interaction of reflection waves generated at sites of heterogeneities occurring at the scale of grains. Lastly, we develop a power-law description for the magnitude of scattering versus characteristic length, which provides a statistical framework for assessing the required number of grains per characteristic specimen dimension for minimizing scatter within these two experimental configurations (NPI, PSPI).

36 MATERIALS SCIENCE↗

NASA NTRSDate not supplied

Improving the Representation of Land Surface Processes using the Data Assimilation Research Testbed (DART)

The land surface is a critical part of the earth system as processes related to water, carbon, energy and nitrogen cycling have important implications for climate forcing, air quality, water availability and seasonal atmospheric forecasting. Despite advances in land surface modeling, land surface model performance is often limited because of errors related to initial and boundary conditions, model structure, and parameters. Data assimilation (DA) techniques combined with an expanding network of earth system observations present an opportunity to reduce these errors and improve simulations. Here we apply an Ensemble Kalman Filter DA system as part of the Data Assimilation Research Testbed to a variety of land surface simulations. First, we describe the use of remotely sensed biomass observations to provide improved simulations of plant phenology, carbon and water cycling for regions highly sensitive to climate change (Western US, China, and Arctic). We discuss approaches to account for systemic biases between models and observations, including the use of spatially-varying adaptive ensemble inflation as an alternative approach to re-scaling soil moisture observations. Finally, we discuss a strategy to incorporate complementary observations (snow water equivalent, solar-induced fluorescence) to better constrain the representation of carbon and water cycling across complex terrain.

DART↗

DOE OSTI2023

The 1996 Mid-Atlantic Winter Flood: Exploring Climate Risk through a Storyline Approach

Abstract This article explores the application of thermodynamic perturbations to a historical midlatitude, wintertime, rain-on-snow flood event to evaluate how similar events may evolve under different climate forcings. In particular, we generate a hindcast of the 1996 Mid-Atlantic flood using an ensemble of 14-km variable-resolution simulations completed with the U.S. Department of Energy’s global Energy Exascale Earth System Model (E3SM). We show that the event is skillfully reproduced over the Susquehanna River Basin (SRB) by E3SM when benchmarked against in situ observational data and high-resolution reanalyses. In addition, we perform five counterfactual experiments to simulate the flood under preindustrial conditions and four different levels of warming as projected by the Community Earth System Model Large Ensemble. We find a nonlinear response in simulated surface runoff and streamflow as a function of atmospheric warming. This is attributed to changing contributions of liquid water input from a shallower initial snowpack (decreased snowmelt), increased surface temperatures and rainfall rates, and increased soil water storage. Flooding associated with this event peaks from around +1 to +2 K of global average surface warming and decreases with additional warming beyond this. There are noticeable timing shifts in peak runoff and streamflow associated with changes in the flashiness of the event. This work highlights the utility of using storyline approaches for communicating climate risk and demonstrates the potential nonlinearities associated with hydrologic extremes in areas that experience ephemeral snowpack, such as the SRB.

Meteorology & Atmospheric Sciences↗

DOE OSTI2022

Simulation of ENSO Teleconnections to Precipitation Extremes over the United States in the High-Resolution Version of E3SM

Abstract We evaluate the simulated teleconnection of El Niño–Southern Oscillation (ENSO) to winter season precipitation extremes over the United States in a long (98 years) 1950 control high-resolution version (HR; 25-km nominal atmosphere model horizontal resolution) of the U.S. Department of Energy’s (DOE) Energy Exascale Earth System Model version 1 (E3SMv1). The model bias and spatial pattern of ENSO teleconnections to mean and extreme precipitation in HR overall are similar to the low-resolution model’s (LR; 110 km) historical simulation (four-member ensemble, 1925–59). However, over the southeastern United States (SE-U.S.), HR produces stronger El Niño–associated extremes, reducing LR’s model bias. Both LR and HR produce weaker than observed increase in storm track activity during El Niño events there, but HR improves the ENSO-associated variability of moisture transport over SE-U.S. During El Niño, stronger vertical velocities in HR produce stronger large-scale precipitation, causing larger latent heating of the troposphere that pulls in more moisture from the Gulf of Mexico into the SE-U.S. This positive feedback also contributes to the stronger mean and extreme precipitation response in HR. Over the Pacific Northwest, LR’s bias of stronger than observed La Niña associated extremes is amplified in HR. Both models simulate stronger than observed moisture transport from the Pacific Ocean into the region during La Niña years. The amplified HR bias there is due to stronger orographically driven vertical updrafts that create stronger large-scale precipitation, despite weaker La Niña–induced storm track activity. Significance Statement New high-resolution Earth system models (ESMs) solve mathematical equations of fluid flow at much smaller spatial scales than prevalent ESMs, and thus are prohibitively expensive to compute. However, they can be useful for simulating accurate details of regional climate extremes that are driven by naturally occurring climate oscillations like El Niño–Southern Oscillation (ENSO). Here, we evaluate the simulation of ENSO-driven precipitation extremes over the United States in the high-resolution version of the U.S. Department of Energy’s new Energy Exascale Earth System Model version 1. We find that the high-resolution model improves upon its low-resolution counterpart over the southeastern United States by producing a better transport of moisture into the region from the Gulf of Mexico during El Niño. Over the U.S. Pacific Northwest, the high-resolution model simulates the atmospheric flow in more detail over the complex mountainous terrain. However, it also brings in more moisture from the Pacific Ocean just like the low-resolution model. This causes it to produce precipitation extremes during La Niña years there that are stronger than that observed in the real world.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2016

The Effect of Future Ambient Air Pollution on Human Premature Mortality to 2100 Using Output from the ACCMIP Model Ensemble

Ambient air pollution from ground-level ozone and fine particulate matter (PM(sub 2.5)) is associated with premature mortality. Future concentrations of these air pollutants will be driven by natural and anthropogenic emissions and by climate change. Using anthropogenic and biomass burning emissions projected in the four Representative Concentration Pathway scenarios (RCPs), the ACCMIP ensemble of chemistry climate models simulated future concentrations of ozone and PM(sub 2.5) at selected decades between 2000 and 2100. We use output from the ACCMIP ensemble, together with projections of future population and baseline mortality rates, to quantify the human premature mortality impacts of future ambient air pollution. Future air-pollution-related premature mortality in 2030, 2050 and 2100 is estimated for each scenario and for each model using a health impact function based on changes in concentrations of ozone and PM(sub 2.5) relative to 2000 and projected future population and baseline mortality rates. Additionally, the global mortality burden of ozone and PM(sub 2.5) in 2000 and each future period is estimated relative to 1850 concentrations, using present-day and future population and baseline mortality rates. The change in future ozone concentrations relative to 2000 is associated with excess global premature mortality in some scenarios/periods, particularly in RCP8.5 in 2100 (316 thousand deaths per year), likely driven by the large increase in methane emissions and by the net effect of climate change projected in this scenario, but it leads to considerable avoided premature mortality for the three other RCPs. However, the global mortality burden of ozone markedly increases from 382000 (121000 to 728000) deaths per year in 2000 to between 1.09 and 2.36 million deaths per year in 2100, across RCPs, mostly due to the effect of increases in population and baseline mortality rates. PM(sub 2.5) concentrations decrease relative to 2000 in all scenarios, due to projected reductions in emissions, and are associated with avoided premature mortality, particularly in 2100: between 2.39 and 1.31 million deaths per year for the four RCPs. The global mortality burden of PM(sub 2.5) is estimated to decrease from 1.70 (1.30 to 2.10) million deaths per year in 2000 to between 0.95 and 1.55 million deaths per year in 2100 for the four RCPs due to the combined effect of decreases in PM(sub 2.5) concentrations and changes in population and baseline mortality rates. Trends in future air-pollution-related mortality vary regionally across scenarios, reflecting assumptions for economic growth and air pollution control specific to each RCP and region. Mortality estimates differ among chemistry climate models due to differences in simulated pollutant concentrations, which is the greatest contributor to overall mortality uncertainty for most cases assessed here, supporting the use of model ensembles to characterize uncertainty. Increases in exposed population and baseline mortality rates of respiratory diseases magnify the impact on premature mortality of changes in future air pollutant concentrations and explain why the future global mortality burden of air pollution can exceed the current burden, even where air pollutant concentrations decrease.

Silva, Raquel A.↗

NASA NTRS2022

Model Evaluation of Short-Lived Climate Forcers for the Arctic Monitoring and Assessment Programme: A Multi-Species, Multi-Model Study

While carbon dioxide is the main cause for global warming, modeling short-lived climate forcers (SLCFs) such as methane, ozone, and particles in the Arctic allows us to simulate near-term climate and health impacts for a sensitive, pristine region that is warming at 3 times the global rate. Atmospheric modeling is critical for understanding the long-range transport of pollutants to the Arctic, as well as the abundance and distribution of SLCFs throughout the Arctic atmosphere. Modeling is also used as a tool to determine SLCF impacts on climate and health in the present and in future emissions scenarios. In this study, we evaluate 18 state-of-the-art atmospheric and Earth system models by assessing their representation of Arctic and Northern Hemisphere atmospheric SLCF distributions, considering a wide range of different chemical species (methane, tropospheric ozone and its precursors, black carbon, sulfate, organic aerosol, and particulate matter) and multiple observational datasets. Model simulations over 4 years (2008-2009 and 2014-2015) conducted for the 2022 Arctic Monitoring and Assessment Programme (AMAP) SLCF assessment report are thoroughly evaluated against satellite, ground, ship, and aircraft-based observations. The annual means, seasonal cycles, and 3-D distributions of SLCFs were evaluated using several metrics, such as absolute and percent model biases and correlation coefficients. The results show a large range in model performance, with no one particular model or model type performing well for all regions and all SLCF species. The multi-model mean (mmm) was able to represent the general features of SLCFs in the Arctic and had the best overall performance. For the SLCFs with the greatest radiative impact (CH4, 03, BC, and SO(sup 2-)(sub 4)), the mmm was within ±25 % of the measurements across the Northern Hemisphere. Therefore, we recommend a multi-model ensemble be used for simulating climate and health impacts of SLCFs. Of the SLCFs in our study, model biases were smallest for C"4 and greatest for OA. For most SLCFs, model biases skewed from positive to negative with increasing latitude. Our analysis suggests that vertical mixing, long-range transport, deposition, and wildfires remain highly uncertain processes. These processes need better representation within atmospheric models to improve their simulation of SLCFs in the Arctic environment. As model development proceeds in these areas, we highly recommend that the vertical and 3-D distribution of SLCFs be evaluated, as that information is critical to improving the uncertain processes in models.

Arctic Monitoring and Assessment Programme↗

DOE OSTI2021

The Climate Response to Emissions Reductions Due to COVID-19: Initial Results From CovidMIP

Many nations responded to the COVID-19 pandemic by restricting travel and other activities during 2020, resulting in temporarily reduced emissions of CO 2 , other greenhouse gases and ozone and aerosol precursors. We perform a coordinated Intercomparison, CovidMIP, of Earth System model simulations to assess the impact on climate of these emissions reductions. Eleven models performed multiple initial-condition ensembles to produce over 280 simulations spanning both initial condition and model structural uncertainty. We find model consensus on reduced aerosol amounts (particularly over East Asia) and associated increases in surface shortwave radiation levels. However, any impact on near-surface temperature or rainfall during 2020-2024 is extremely small and is not detectable in this initial analysis. Regional analyses on a finer scale, and closer attention to extremes (especially linked to changes in atmospheric composition and air quality) are required to test the impact of COVID-19-related emission reductions on near-term climate.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2012

The Influence of Internal Model Variability in GEOS-5 on Interhemispheric CO2 Exchange

An ensemble of eight atmospheric CO2 simulations was completed employing the National Aeronautics and Space Administration (NASA) Goddard Earth Observation System, Version 5 (GEOS-5) for the years 2000-2001, each with initial meteorological conditions corresponding to different days in January 2000 to examine internal model variability. Globally, the model runs show similar concentrations of CO2 for the two years, but in regions of high CO2 concentrations due to fossil fuel emissions, large differences among different model simulations appear. The phasing and amplitude of the CO2 cycle at Northern Hemisphere locations in all of the ensemble members is similar to that of surface observations. In several southern hemisphere locations, however, some of the GEOS-5 model CO2 cycles are out of phase by as much as four months, and large variations occur between the ensemble members. This result indicates that there is large sensitivity to transport in these regions. The differences vary by latitude-the most extreme differences in the Tropics and the least at the South Pole. Examples of these differences among the ensemble members with regard to CO2 uptake and respiration of the terrestrial biosphere and CO2 emissions due to fossil fuel emissions are shown at Cape Grim, Tasmania. Integration-based flow analysis of the atmospheric circulation in the model runs shows widely varying paths of flow into the Tasmania region among the models including sources from North America, South America, South Africa, South Asia and Indonesia. These results suggest that interhemispheric transport can be strongly influenced by internal model variability.

Allen, Melissa↗

NASA NTRS2021

The Climate Response to Emissions Reductions due to COVID-19: Initial Results from CovidMIP

Many nations responded to the COVID-19 pandemic by restricting travel and other activities during 2020, resulting in temporarily reduced emissions of CO2, other greenhouse gases and ozone and aerosol precursors. We present the initial results from a coordinated Intercomparison, CovidMIP, of Earth system model simulations which assess the impact on climate of these emissions reductions. Twelve models performed multiple initial-condition ensembles to produce over 300 simulations spanning both initial condition and model structural uncertainty. We find model consensus on reduced aerosol amounts (particularly over southern and eastern Asia) and associated increases in surface shortwave radiation levels. However, any impact on near-surface temperature or rainfall during 2020-2024 is extremely small and is not detectable in this initial analysis. Regional analyses on a finer scale, and closer attention to extremes (especially linked to changes in atmospheric composition and air quality) are required to test the impact of COVID- 19-related emission reductions on near-term climate.

Chris D. Jones↗

DOE OSTI2020

A generalized approximate control variate framework for multifidelity uncertainty quantification

We describe and analyze a variance reduction approach for Monte Carlo (MC) sampling that accelerates the estimation of statistics of computationally expensive simulation models using an ensemble of models with lower cost. These lower cost models — which are typically lower fidelity with unknown statistics — are used to reduce the variance in statistical estimators relative to a MC estimator with equivalent cost. We derive the conditions under which our proposed approximate control variate framework recovers existing multifidelity variance reduction schemes as special cases. We demonstrate that existing recursive/nested strategies are suboptimal because they use the additional low-fidelity models only to efficiently estimate the unknown mean of the first low-fidelity model. As a result, they cannot achieve variance reduction beyond that of a control variate estimator that uses a single low-fidelity model with known mean. However, there often exists about an order-of-magnitude gap between the maximum achievable variance reduction using all low-fidelity models and that achieved by a single low-fidelity model with known mean. We show that our proposed approach can exploit this gap to achieve greater variance reduction by using non-recursive sampling schemes. The proposed strategy reduces the total cost of accurately estimating statistics, especially in cases where only low-fidelity simulation models are accessible for additional evaluations. Several analytic examples and an example with a hyperbolic PDE describing elastic wave propagation in heterogeneous media are used to illustrate the main features of the methodology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Reducing uncertainty of polar to midlatitude linkages using DOE’s E3SM in a coordinated model-experiment setting

This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1997

An Application of the "Virtual Spacecraft" Concept in Evaluation of the Mars Pathfinder Lander Low Gain Antenna

The virtual spacecraft concept is embodied in a set of subsystems, either in the form of hardware or computational models, which together represent all, or a portion of, a spacecraft. For example, the telecommunications transponder may be a hardware prototype while the propulsion system may exist only as a simulation. As the various subsystems are realized in hardware, the spacecraft becomes progressively less virtual. This concept is enabled by JPL's Mission System Testbed which is a set of networked workstations running a message passing operating system called "TRAMEL" which stands for Task Remote Asynchronous Message Exchange Layer. Each simulation on the workstations, which may in fact be hardware controlled by the workstation, "publishes" its operating parameters on TRAMEL and other simulations requiring those parameters as input may "subscribe" to them. In this manner, the whole simulation operates as a single virtual system. This paper describes a simulation designed to evaluate a communications link between the earth and the Mars Pathfinder Lander module as it descends under a parachute through the Martian atmosphere toward the planet's surface. This link includes a transmitter and a low gain antenna on the spacecraft and a receiving antenna and receiver on the earth as well as a simulation of the dynamics of the spacecraft. The transmitter, the ground station antenna, the receiver and the dynamics are all simulated computationally while the spacecraft antenna is implemented in hardware on a very simple spacecraft mockup. The dynamics simulation is a record of one output of the ensemble of outputs of a Monte Carlo simulation of the descent. Additionally, the antenna/spacecraft mock-up system was simulated using APATCH, a shooting and bouncing ray code developed by Demaco, Inc. The antenna simulation, the antenna hardware, and the link simulation are all physically located in different facilities at JPL separated by several hundred meters and are linked via the local area network (LAN).

Pogorzelski, R. J.↗