SEARCH · Engineering Papers
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.
The grant supported the development of new tools to accelerate the spin up of the marine biogeochemical component of earth system models. The grant also supported the development of new diagnostic tools based on the application of these tools.
Antarctic coastal polynyas produce dense shelf water, a primary source of Antarctic Bottom Water that contributes to the global overturning circulation. This paper investigates Antarctic dense water formation in the high-resolution version of the Energy Exascale Earth System Model (E3SM-HR). The model is able to reproduce the main Antarctic coastal polynyas, although the polynyas are smaller in area compared to observations. E3SM-HR also simulates several occurrences of open-ocean polynyas (OOPs) in the Weddell Sea at a higher rate than what the last 50 years of the satellite sea ice observational record suggests, but similarly to other high-resolution Earth system model simulations. Furthermore, the densest water masses in the model are formed within the OOPs rather than on the continental shelf as is typically observed. Biases related to the lack of dense water formation on the continental shelf are associated with overly strong atmospheric polar easterlies, which lead to a strong Antarctic Slope Front and too little exchange between on- and off-continental shelf water masses. Strong polar easterlies also produce excessive southward Ekman transport, causing a build-up of sea ice over the continental shelf and enhanced ice melting in the summer season. This, in turn, produces water masses on the continental shelf that are overly fresh and less dense relative to observations. Our results indicate that high resolution alone is insufficient for models to properly reproduce Antarctic dense water; the large-scale polar atmospheric circulation around Antarctica must also be accurately simulated.
Understanding the interaction between humans and the Earth system is a computationally daunting task, with many possible approaches depending on resources available and questions of interest. For example, state-of-the-art impact models require decade-long time series of relatively high frequency, spatially resolved and often multiple variables representing climatic impact-drivers (Ruane et al., 2022). Most commonly these are derived from the outputs of detailed, computationally expensive Earth System Models (ESMs) run according to a standard, limited set of future scenarios, the latest being the SSP-RCPs run under CMIP6/ScenarioMIP (Eyring et al., 2016; O’Neill et al., 2016). At the time of writing, O’Neill et al. (2016) has been cited more than 1750 times and Eyring et al. (2016) more than 5000 times, highlighting the broad, general applications of this data. Often, however, impact modeling seeks to explore new scenarios that were not part of the ScenarioMIP protocol, and/or needs a larger set of initial condition ensemble members than are typically available to quantify the effects of ESM internal variability. In addition, the recognition that the human and Earth systems are fundamentally intertwined, and may feature potentially significant feedback loops, is making integrated, simultaneous modeling of the coupled human-Earth system increasingly necessary, if computationally challenging with most existing tools (Thornton et al., 2017). For impact modelers, climate model emulators can be the answer to meet both the needs of: 1) creating realizations for novel scenarios and 2) achieving a simplified, computationally tractable representation of ESM behavior in a coupled human-Earth system modeling framework. We proposed a new, comprehensive approach to such emulation of gridded, multivariate ESM outputs for novel scenarios without the computational cost of a full ESM, STITCHES (Tebaldi et al., 2022). The approach outlined in Tebaldi et al. (2022) should be extensible to future CMIP eras, although the STITCHES software at present is strictly focused on CMIP6/ScenarioMIP data hosted on Pangeo (https://gallery.pangeo.io/repos/pangeo-gallery/cmip6/). The corresponding STITCHES Python package uses existing archives of ESMs’ scenario experiments from CMIP6/ScenarioMIP to construct gridded, multivariate realizations of new scenarios provided by reduced complexity climate models (Hartin et al., 2015; Meinshausen et al., 2011; Smith et al., 2018), or to enrich existing initial condition ensembles. Its output provides the same characteristics as the emulated ESM output: multivariate (spanning potentially all variables that the ESM has saved), spatially resolved (down to the native grid of the ESM), and preserving the same high frequency as the original data. A new realization of multiple variables can be generated on the order of minutes with STITCHES, rather than the hours or sometimes days that ESMs require.
Abstract The update to the deep convection parameterization in the energy exascale Earth system model version 2 (E3SMv2) makes the simulated quasi‐biennial oscillation (QBO) have a shorter period than without the update when other tunable parameters the same since the update makes convection more intense but less frequent while leaving the time‐mean convective heating tendency almost unchanged. Amplitudes of momentum fluxes of parameterized gravity waves (GWs) are intensified since they are determined by the square of convective heating tendency. In addition, stronger planetary waves are simulated in E3SMv2 with the convection scheme update, partially due to the tropospheric precipitation change. Furthermore, there is evidence that planetary waves are intensified in the stratosphere due to the dissipation of enhanced parameterized GWs. These factors are found to be responsible for modulating the QBO simulated in E3SMv2.
Abstract Machine learning (ML) is a revolutionary technology with demonstrable applications across multiple disciplines. Within the Earth science community, ML has been most visible for weather forecasting, producing forecasts that rival modern physics‐based models. Given the importance of deepening our understanding and improving predictions of the Earth system on all time scales, efforts are now underway to develop Earth‐system models (ESMs) capable of representing all components of the coupled Earth system (or their aggregated behavior) and their response to external changes over long timescales. Building trust in ESMs is a much more difficult problem than for weather forecast models, not least because the model must represent the alternate (e.g., future or paleoclimatic) coupled states of the system for which there are no direct observations. Given that the physical principles that enable predictions about the response of the Earth system are often not explicitly coded in these ML‐based models, demonstrating the credibility of ML‐based ESMs thus requires us to build evidence of their consistency with the physical system. To this end, this paper puts forward five recommendations to enhance comprehensive, standardized, and independent evaluation of ML‐based ESMs to strengthen their credibility and promote their wider use.
Surface soil moisture (mrsos) and vertically integrated soil moisture (mrsol) over the top 10 cm should, by definition, be physically consistent in Earth System Models (ESMs). However, an evaluation of nine CMIP6 models reveals substantial inconsistencies: in some models, mrsos and integrated mrsol agree globally; in others, they align only in specific regions; and in a few, they diverge across all grid cells. These discrepancies arise from a combination of factors, including metadata errors, inconsistent variable definitions, or diagnostic sequencing within the model. We demonstrate how such issues can lead to significant biases, even when both variables are present and seemingly well-defined. As model complexity increases and multi-model comparisons become more common, assumptions about variable equivalence may lead to flawed conclusions. This study highlights the need for routine consistency checks, improved metadata standards, and community-wide practices that ensure reliability of derived variables across ESM outputs, particularly in preparation for CMIP7.
Bias correction is a crucial step in using Earth system model outputs for assessments, as it adjusts systematic errors by comparing the model to observations. However, standard methods – ranging from mean-based linear scaling to distribution-based quantile mapping typically treat bias correction as a single-scale process, overlooking the fact that biases can manifest differently across daily, seasonal, and annual timescales. In this study, we propose a novel, timescale-aware bias-correction approach built on Empirical Mode Decomposition. By decomposing the meteorological signal into multiple oscillatory components and aggregating them to represent distinct timescales, we apply targeted corrections to each component, thereby preserving both short- and long-term structure in the data. Experimental illustrations show that the timescale-aware EMDBC framework matches the performance of conventional quantile-delta mapping (QDM) at the native daily scale and achieves progressively larger bias reductions at bi-weekly, seasonal, and annual scales. As a result, the proposed approach offers a more robust path to accurate and reliable Earth system projections, strengthening their utility for resilience and adaptation planning.
Lateral groundwater flow (LGF) is an important hydrologic process in controlling water table dynamics. Due to the relatively coarse spatial resolutions of land surface models, the representation of this process is often overlooked or overly simplified. In this study, we developed a hillslope-based lateral groundwater flow model. Specifically, we first developed a hillslope definition model based on an existing watershed delineation model to represent the subgrid spatial variability in topography. Building upon this hillslope definition, we then developed a physical-based lateral groundwater flow using Darcy’s equation. This model explicitly considers the relationships between the groundwater table along the hillslope and the river water table levels. We coupled this intra-grid model to the land component (E3SM Land Model: ELM) and river component (MOdel for Scale Adaptive River Transport: MOSART) of the Energy Exascale Earth System Model (E3SM). We tested both the hillslope definition model and the lateral groundwater flow model and performed sensitivity experiments using different configurations. Simulations for a single grid cell at 0.5°×0.5° within the Amazon basin show that the definition of hillslope is the key to modeling lateral flow processes and the runoff partition between surface and subsurface can be dramatically changed using the hillslope approach. Although our method provides a pathway to improve the lateral flow process, future improvements are needed to better capture the subgrid structure to account for the spatial variability in hillslopes within the simulated grid of land surface models.
The second-generation Sup3rCC dataset provides high-resolution meteorological data generated through the downscaling of multiple earth system models (ESMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6). This downscaling is performed through application of a generative machine learning approach called Super-Resolution for Renewable Resource Data (sup3r). This dataset builds on the first-generation Sup3rCC data by applying improved bias correction methods and adding downscaled precipitation to the output variables. In this presentation, we explore the output characteristics of the dataset and various validation analyses. We also present and discuss plans for the integration of this data into power system planning models using a decision-making under deep uncertainty (DMDU) methodology.
Abstract The water cycle is an important component of the earth system and it plays a key role in many facets of society, including energy production, agriculture, and human health and safety. In this study, the Energy Exascale Earth System Model version 1 (E3SMv1) is run with low‐resolution (roughly 110 km) and high‐resolution (roughly 25 km) configurations—as established by the High Resolution Model Intercomparison Project protocol—to evaluate the atmospheric and terrestrial water budgets over the conterminous United States (CONUS) at the large watershed scale. The warm season water cycle slows down in the HR experiment relative to the LR, with decreasing fluxes of precipitation, evapotranspiration, atmospheric moisture convergence, and runoff. The reductions in these terms exacerbate biases for some watersheds, while reducing them in others. For example, precipitation biases are exacerbated at HR over the Eastern and Central CONUS watersheds, while precipitation biases are reduced at HR over the Western CONUS watersheds. The most pronounced changes with resolution to the water cycle come from reductions in precipitation and evapotranspiration. The reduction in evapotranspiration reduces the biases across nearly all of the CONUS. Additional exploratory metrics show improvements to water cycle extremes (both in precipitation and streamflow), fractional contributions of different storm types to total precipitation, and mountain snowpack.
Numerical water tracers are implemented in the Energy Exascale Earth System Model version 2. Simulations performed with the water‐tag‐enabled model for both pre‐industrial and future greenhouse gas concentrations reveal a marked increase in the role of mid‐latitude and southern subtropical regions as exporters of atmospheric moisture—to the extratropical upper troposphere and the tropical free troposphere. For the latter, the northward shift of the Intertropical Convergence Zone increases cross‐hemispheric transport of subtropical water vapor to the Northern Hemisphere. In the polar regions, most of the lower tropospheric moistening instead arises from increases in local evaporation. These findings illustrate the utility of the water tags, underscore critical changes in global hydrologic cycle, and provide insight into atmospheric dynamics under future climate scenarios. For applications when a global grid is desired, we additionally propose a novel statistical reconstruction, based on copula modeling, of the joint distribution of origin of water vapor, which reduces the number of tracers from order $\mathcal{O}\left({n}^{2}\right)$to order $\mathcal{O}(n)$, substantially ameliorating the considerable computational cost of water tracers. This statistical reconstruction is particularly beneficial to the interpretation of the relationship between latitude and longitude of origin of moisture over the tropical oceans and in the lower troposphere over land.
This paper describes the atmospheric component of the US Department of Energy's Energy Exascale Earth System Model (E3SM) version 3. Significant updates have been made to the atmospheric physics compared to earlier versions. Specifically, interactive gas chemistry has been implemented, along with improved representations of aerosols and dust emissions. A new stratiform cloud microphysics scheme more physically treats ice processes and aerosol‐cloud interactions. The deep convection parameterization has been largely improved with sophisticated microphysics for convective clouds, making model convection sensitive to large‐scale dynamics, and incorporating the dynamical and physical effects of organized mesoscale convection. Improvements in aerosol wet removal processes and parameter re‐tuning of key aerosol and cloud processes have improved model aerosol radiative forcing. The model's vertical resolution has increased from 72 to 80 layers with the extra eight layers added in the lower stratosphere to better simulate the Quasi‐Biennial Oscillation. These improvements have enhanced E3SM's capability to couple aerosol, chemistry, and biogeochemistry and reduced some long‐standing biases in simulating tropical variability. Compared to its predecessors, the model shows a much stronger signal for the Madden‐Julian Oscillation, Kelvin waves, mixed Rossby‐gravity waves, and eastward inertia‐gravity waves. Aerosol radiative forcing has been considerably reduced and is now better aligned with community best estimates, leading to significantly improved skill in simulating historical temperature records. Its simulated mean‐state climate is largely comparable to E3SMv2, but with some notable degradation in shortwave cloud radiative effect, precipitable water, and surface wind stress, which will be addressed in future updates.
Increasing concentrations of CO 2 in the atmosphere influence climate both through CO 2 ’s role as a greenhouse gas and through its impact on plants. Plants respond to atmospheric CO 2 concentrations in several ways that can alter surface energy and water fluxes and thus surface climate, including changes in stomatal conductance, water use, and canopy leaf area. These plant physiological responses are already embedded in most Earth system models, and a robust literature demonstrates that they can affect global-scale temperature. However, the physiological contribution to transient warming has yet to be assessed systematically in Earth system models. Here this gap is addressed using carbon cycle simulations from phases 5 and 6 of the Coupled Model Intercomparison Project (CMIP) to isolate the radiative and physiological contributions to the transient climate response (TCR), which is defined as the change in globally averaged near-surface air temperature during the 20-yr window centered on the time of CO 2 doubling relative to preindustrial CO 2 concentrations. In CMIP6 models, the physiological effect contributes 0.12°C ( σ : 0.09°C; range: 0.02°–0.29°C) of warming to the TCR, corresponding to 6.1% of the full TCR ( σ : 3.8%; range: 1.4%–13.9%). Moreover, variation in the physiological contribution to the TCR across models contributes disproportionately more to the intermodel spread of TCR estimates than it does to the mean. The largest contribution of plant physiology to CO 2 -forced warming—and the intermodel spread in warming—occurs over land, especially in forested regions.
We posit that AI methods can be leveraged to significantly enhance the predictive skill of forward Earth system modeling (ESM) activities. A hybrid framework incorporating traditional ESM modeling, inference methods, and AI techniques could make better use of both measured and computed information as well as computational resources by targeting inference tasks at program priorities, such as the predictability of precipitation extremes. This runtime pathway to closing the simulation/analysis-data/model improvement loop will streamline the traditional offline pathway to model improvement, which is based on domain science expertise, while suggesting guidance for further observations and measurements.
The focus of this project is to improve the fidelity of the radiation parameterization used in earth system models to deliver increased accuracy with a cheaper computational approach.
This annual report describes the activities of PCMDI – An Earth System Model Evaluation Project for the period November 1, 2019 – June 30, 2020. PCMDI is a project at Lawrence Livermore National Laboratory that is funded by the Climate and Environmental Sciences Division (CESD) of the Department of Energy. PCMDI is funded by two efforts within CESD: The Regional and Global Model Analysis program area (Dr. Renu Joseph, program manager) and the Atmospheric Systems Research program (Dr. Shaima Nasiri and Jeffrey Stehr, program managers). As instructed by Drs. Joseph, Nasiri, and Stehr, we provide an annual report addressing sections 4 (Progress and Plans) and 5 (Staffing and Budget Summary) of the SFA Program Management and Performance Reporting Criteria.
This subcontract was used to evaluate the convection-permitting version of the Energy Exascale Earth System Model (E3SM), generally referred to as the simple Cloud-resolving E3SM Atmosphere Model (SCREAM). Three projects were addressed and briefly summarized..
Abstract This paper provides a comprehensive derivation of the total energy equations for the atmospheric components of Earth System Models (ESMs). The assumptions and approximations made in this derivation are motivated and discussed. In particular, it is emphasized that closing the energy budget is conceptually challenging and hard to achieve in practice without resorting to ad hoc fixers. As a concrete example, the energy budget terms are diagnosed in a realistic climate simulation using a global atmosphere model. The largest total energy errors in this example are spurious dynamical core energy dissipation, thermodynamic inconsistencies (e.g., coupling parameterizations with the host model) and missing processes/terms associated with falling precipitation and evaporation (e.g., enthalpy flux between components). The latter two errors are not, in general, reduced by increasing horizontal resolution. They are due to incomplete thermodynamic and dynamic formulations. Future research directions are proposed to reconcile and improve thermodynamics formulations and conservation principles.