Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earth System Models”

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 829 records · Page 46

Investigating coastal backwater effects and flooding in the coastal zone using a global river transport model on an unstructured mesh

Abstract. Coastal backwater effects are caused by the downstream water level increase as a result of elevated sea level, high river discharge and their compounding influence. Such effects have crucial impacts on floods in densely populated regions but have not been well represented in large-scale river models used in Earth system models (ESMs), partly due to model mesh deficiency and oversimplifications of river hydrodynamics. Using two mid-Atlantic river basins as a testbed, we perform the first attempt to simulate the backwater effects comprehensively over a coastal region using the MOSART river transport model under an ESM framework, i.e., Energy Exascale Earth System Model (E3SM) configured on a regionally refined unstructured mesh, with a focus on understanding the backwater drivers and their long-term variations. By including sea level variations at the river downstream boundary, the model performance in capturing backwaters is greatly improved. We also propose a new flood event selection scheme to facilitate the decomposition of backwater drivers into different components. Our results show that while storm surge is a key driver, the influence of extreme discharge cannot be neglected, particularly when the river drains to a narrow river-like estuary. Compound flooding, while not necessarily increasing the flood peaks, exacerbates the flood risk by extending the duration of multiple coastal and fluvial processes. Furthermore, our simulations and analysis highlight the increasing strength of backwater effects due to sea level rise and more frequent storm surge during 1990–2019. Thus, backwaters need to be properly represented in ESMs to improve the predictive understanding of coastal flooding.

54 ENVIRONMENTAL SCIENCES↗

The Polar System Analysis Package (Ver. 1 Specifications)

The Polar System Analysis Package (POLeSTAr) is a modern data fusion tool being designed to benchmark,evaluate and compare the polar physics and biogeochemistry of Earth System Models (ESMs). It is being written in Python to make use of libraries including Xarray, Matplotlib, SciPy, NumPy and PyTorch within algorithms developed specifically for polar system modeling. Initial POLeSTAr algorithms will be adapted from the MATLAB Ridgepack package that was designed to analyze the Model for Prediction Across Scales (MPAS) sea ice component of the Energy Exascale Earth System Model (E3SM) and sea ice model CICE within the Community Earth System Model (CESM).

97 MATHEMATICS AND COMPUTING↗

Complex functionality with minimal computation: Promise and pitfalls of reduced‐tracer ocean biogeochemistry models

Abstract Earth System Models increasingly include ocean biogeochemistry models in order to predict changes in ocean carbon storage, hypoxia, and biological productivity under climate change. However, state‐of‐the‐art ocean biogeochemical models include many advected tracers, that significantly increase the computational resources required, forcing a trade‐off with spatial resolution. Here, we compare a state‐of‐the art model with 30 prognostic tracers (TOPAZ) with two reduced‐tracer models, one with 6 tracers (BLING), and the other with 3 tracers (miniBLING). The reduced‐tracer models employ parameterized, implicit biological functions, which nonetheless capture many of the most important processes resolved by TOPAZ. All three are embedded in the same coupled climate model. Despite the large difference in tracer number, the absence of tracers for living organic matter is shown to have a minimal impact on the transport of nutrient elements, and the three models produce similar mean annual preindustrial distributions of macronutrients, oxygen, and carbon. Significant differences do exist among the models, in particular the seasonal cycle of biomass and export production, but it does not appear that these are necessary consequences of the reduced tracer number. With increasing CO 2 , changes in dissolved oxygen and anthropogenic carbon uptake are very similar across the different models. Thus, while the reduced‐tracer models do not explicitly resolve the diversity and internal dynamics of marine ecosystems, we demonstrate that such models are applicable to a broad suite of major biogeochemical concerns, including anthropogenic change. These results are very promising for the further development and application of reduced‐tracer biogeochemical models that incorporate “sub‐ecosystem‐scale” parameterizations.

58 GEOSCIENCES↗

Intensified Soil Moisture Extremes Decrease Soil Organic Carbon Decomposition: A Mechanistic Modeling Analysis

Earth system models have predicted that there will be more frequent and severe precipitation and drought events in terrestrial ecosystems. Microbially mediated decomposition of soil organic carbon (SOC) tends to increase as soils wet and decrease as soils dry. However, the long-term SOC change under intensified moisture extremes remains poorly known as it depends on the frequency and intensity of soil drying and wetting. In this study, we explored long-term SOC dynamics under scenarios of alternating drying-wetting cycles using the Microbial-ENzyme Decomposition model, a mechanistic microbial model. The model was parameterized with 11 years of observations from a temperate deciduous broadleaf forest site, showing satisfactory model performance in both model calibration (R 2 = 0.67) and validation (R 2 = 0.69) against heterotrophic respiration. We then used the model to simulate the long-term SOC dynamics under five scenarios of alternating drying-wetting cycles with different frequencies and severities over a period of 100 years. Results showed that the changes in active microbial biomass C and the corresponding turnover rates of SOC pools were more sensitive to soil drying than soil wetting. As a result, the cumulative soil carbon emission from microbial respiration decreased by 433.7 g C m -2 after the 100-year simulation in the highest frequency and intensity moisture scenario, but was not significantly affected by the lowest frequency and intensity scenario. This study emphasizes the nonlinear response of SOC decomposition to soil moisture changes, which causes decreased decomposition by microbes under drying that is, not compensated by increased decomposition under wetting conditions.

58 GEOSCIENCES↗

Disentangling the Effects of Vapor Pressure Deficit and Soil Water Availability on Canopy Conductance in a Seasonal Tropical Forest During the 2015 El Niño Drought

In land models used in Earth system models, vegetation water stress is heuristically related to soil water matric potential, and root water uptake in soil layers is simply proportional to root density to meet the transpiration demand. The efficiency of roots utilizing deeper soil water under water stress cannot be well represented, but it is of critical importance for transpiration regulation and carbon flux. To improve modeling of tropical forest response to drought, a plant hydrodynamic scheme (HYDRO), adapted from a trait-based model in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), has been implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) to represent the hydraulic dynamics and root water uptake strategy. The adapted model conforms to the big-leaf photosynthesis model in ELM. Vegetation water stress in the model is a function of prognostic leaf water potential and is used to regulate stomatal conductance. Point simulations were performed at two tropical forest sites in the Amazon and Panama. Model sensitivity to the numerical solver, plant hydraulic traits, soil properties, and atmospheric forcing has been evaluated. Our results showed that, even though water stress factors calculated at our study sites using HYDRO are comparable to that simulated by the heuristic water stress model, soil moisture prediction is improved with soil water redistributions through roots. Simulated results using HYDRO is not as sensitive to plant hydraulic traits measured for different tree species as to the soil properties and atmospheric forcing. Simplified numerical solver may cause inaccuracy as to where roots uptake water. This study suggests that plant hydrodynamic representation is critical for improved simulations of soil moisture under water stresses.

54 ENVIRONMENTAL SCIENCES↗

Assessing Two Approaches for Enhancing the Range of Simulated Scales in the E3SMv1 and the Impact on the Character of Hourly US Precipitation

Abstract Improving the representation of precipitation in Earth system models is essential for understanding and projecting water cycle changes across scales. Progress has been hampered by persistent deficiencies in representing precipitation frequency, intensity, and timing in current models. Here, we analyze simulated US precipitation in the low‐resolution (LR) configuration of the Energy Exascale Earth System Model (E3SMv1) and assess the effect of two approaches to enhance the range of explicitly resolved scales: high‐resolution (HR) and multiscale modeling framework (MMF), which incur similar computational expense. Both E3SMv1‐MMF and E3SMv1‐HR capture more intense and less frequent precipitation on hourly and daily timescales relative to E3SMv1‐LR. E3SMv1‐HR improves the intensity over the Eastern and Northwestern US during winter, while E3SMv1‐MMF improves the intensity over the Eastern US and summer diurnal timing over the Central US. These results indicate that both methods may be needed to improve simulations of different storm types, seasons, and regions.

58 GEOSCIENCES↗

A modular biogeochemical modeling suite for next-generation ocean models (Final Report)

The Marine Biogeochemistry Library (MARBL) is a self-contained implementation of ocean biogeochemistry suitable for incorporation in Earth system models. The default configuration of MARBL represents multiple nutrient co-limitation (N, P, Si, and Fe), includes three explicit phytoplankton functional groups (diatoms, diazotrophs, and pico/nano phytoplankton), and one implicit group (calcifiers); there is one zooplankton group. MARBL includes fully prognostic carbonate chemistry and simulates sinking particulate organic matter implicitly, subject to ballasting by mineral dust, biogenic mineral carbonate and opal. MARBL represents subgrid-scale variations in light, variable C:P stoichiometry, semi-labile and refractory dissolved organic material pools, prognostic oceanic emission of ammonia, and an explicit ligand tracer that complexes iron. MARBL is fully runtime-configurable, and is flexible enough to accommodate a variable number of plankton functional types and ecosystem configurations. Critically, MARBL is implemented to be independent of any particular ocean general circulation model (OGCM); it is a standalone library that interfaces with an OGCM's via an explicit "driver" layer. MARBL can thus be coupled to multiple OGCM's, enabling efficient maintenance, development, and cross-model inter-comparison. MARBL is the ocean biogeochemistry component for the Community Earth System Model and the DOE’s Energy Exascale Earth System Model.

54 ENVIRONMENTAL SCIENCES↗

Urban Landscape and Climate Change

A workshop on how earth system models incorporate the influences of urban environments on climate, with an emphasis on how large observation studies can be used to aid the modeling community, was held on August 27 - 28, 2013 at Argonne National Laboratory. The goal of the workshop was to identify urban processes important in influencing the predictive capability of regional climate models and global earth system models, the impacts of climate change that are important for future urban centers, data available or needed to develop a representative model of urban ecosystems, and how cities can use that information to improve the lives of their residents.

54 ENVIRONMENTAL SCIENCES↗

Climate Base State Influences on South Asian Monsoon Processes Derived From Analyses of E3SMv2 and CESM2

Abstract The effects of differences in climate base state are related to processes associated with the present‐day South Asian monsoon simulations in the Energy Exascale Earth System Model version 2 (E3SMv2) and the Community Earth System Model version 2 (CESM2). Though tropical Pacific and Indian Ocean base state sea surface temperatures (SSTs) are over 1°C cooler in E3SMv2 compared to CESM2, and there is an overall reduction of Indian sector precipitation, the pattern of South Asian monsoon precipitation is similar in the two models. Monsoon‐ENSO teleconnections, dynamically linked by the large‐scale east‐west atmospheric circulation, are reduced in E3SMv2 compared to CESM2. In E3SMv2, this is related to cooler tropical SSTs and ENSO amplitude that is less than half that in CESM2. Comparison to a tropical Pacific pacemaker experiment shows, to a first order, that the base state SSTs and ENSO amplitude contribute roughly equally to lower amplitude monsoon‐ENSO teleconnections in E3SMv2.

54 ENVIRONMENTAL SCIENCES↗

Rainbows and Climate Change: A tutorial on climate model diagnostics and parameterization

Earth System Models (ESMs) must represent processes below the grid scale of a model using representations (parameterizations) of physical and chemical processes. As a tutorial exercise to understand diagnostics and parameterization, this work presents a representation of rainbows for an ESM: the Community Earth System Model version 2 (CESM2). Using the 'state' of the model, basic physical laws, and some assumptions, we generate a representation of this unique optical phenomena as a diagnostic output. Rainbow occurrence and it's possible changes are related to cloud occurrence and rain formation which are critical uncertainties for climate change prediction. The work highlights issues which are typical of many diagnostics parameterizations such as assumptions, uncertain parameters and the difficulty of evaluation against uncertain observations. Results agree qualitatively with limited available global 'observations' of Rainbows. Rainbows are seen in expected locations in the sub-tropics over the ocean where broken clouds and frequent precipitation occurs. The diurnal peak is in the morning over ocean and in the evening over land. The representation of rainbows is found to be quantitatively sensitive to the assumed amount of cloudiness and the amount of stratiform rain. Rainbows are projected to have decreased, mostly in the Northern Hemisphere, due to aerosol pollution effects increasing cloud coverage since 1850. In the future, continued climate change is projected to decrease cloud cover, associated with a positive cloud feedback. As a result the rainbow diagnostic projects that rainbows will increase in the future, with the largest changes at mid-latitudes. The diagnostic may be useful for assessing cloud parameterizations, and is an exercise in how to build and test parameterizations of atmospheric phenomena.

54 ENVIRONMENTAL SCIENCES↗

Evaluating the Simulation of CONUS Precipitation by Storm Type in E3SM

Abstract Conventional low‐resolution (LR) climate models, including the Energy Exascale Earth System Model (E3SMv1), have well‐known biases in simulating the frequency, intensity, and timing of precipitation. Approaches to next‐generation E3SM, whether the high‐resolution (HR) or multiscale modeling framework (MMF) configuration, improve the simulation of the intensity and frequency of precipitation, but regional and seasonal deficiencies still exist. Here we apply a methodology to assess the contribution of tropical cyclones (TCs), extratropical cyclones (ETCs), and mesoscale convective systems (MCSs) to simulated precipitation in E3SMv1‐HR and E3SMv1‐MMF relative to E3SMv1‐LR. Across the United States, E3SMv1‐MMF provides the best simulation in terms of precipitation accumulation, frequency and intensity from MCSs and TCs compared to E3SMv1‐LR and E3SMv1‐HR. All E3SMv1 configurations overestimate precipitation amounts from and the frequency of ETCs over CONUS, with conventional E3SMv1‐LR providing the best simulation compared to observations despite limitations in precipitation intensity within these events.

54 ENVIRONMENTAL SCIENCES↗

Graphics Processing Unit (GPU) Acceleration of the Goddard Earth Observing System Atmospheric Model

The Goddard Earth Observing System 5 (GEOS-5) is the atmospheric model used by the Global Modeling and Assimilation Office (GMAO) for a variety of applications, from long-term climate prediction at relatively coarse resolution, to data assimilation and numerical weather prediction, to very high-resolution cloud-resolving simulations. GEOS-5 is being ported to a graphics processing unit (GPU) cluster at the NASA Center for Climate Simulation (NCCS). By utilizing GPU co-processor technology, we expect to increase the throughput of GEOS-5 by at least an order of magnitude, and accelerate the process of scientific exploration across all scales of global modeling, including: The large-scale, high-end application of non-hydrostatic, global, cloud-resolving modeling at 10- to I-kilometer (km) global resolutions Intermediate-resolution seasonal climate and weather prediction at 50- to 25-km on small clusters of GPUs Long-range, coarse-resolution climate modeling, enabled on a small box of GPUs for the individual researcher After being ported to the GPU cluster, the primary physics components and the dynamical core of GEOS-5 have demonstrated a potential speedup of 15-40 times over conventional processor cores. Performance improvements of this magnitude reduce the required scalability of 1-km, global, cloud-resolving models from an unfathomable 6 million cores to an attainable 200,000 GPU-enabled cores.

Putnam, Williama↗

sympl (v. 0.4.0) and climt (v. 0.15.3) – towards a flexible framework for building model hierarchies in Python

sympl (System for Modelling Planets) andclimt (Climate Modelling and Diagnostics Toolkit) are anattempt to rethink climate modelling frameworks from the ground up. The aimis to use expressive data structures available in the scientific Pythonecosystem along with best practices in software design to allow scientists toeasily and reliably combine model components to represent the climate systemat a desired level of complexity and to enable users to fully understandwhat the model is doing. sympl is a framework which formulates the model in terms of astate that gets evolved forward in time or modified within a specifictime by well-defined components. sympl's design facilitates buildingmodels that are self-documenting, are highly interoperable, and providefine-grained control over model components and behaviour. symplcomponents contain all relevant information about the input they expect andoutput that they provide. Components are designed to be easily interchanged,even when they rely on different units or array configurations.sympl provides basic functions and objects which could be used inany type of Earth system model. climt is an Earth system modelling toolkit that contains scientificcomponents built using sympl base objects. These include both purePython components and wrapped Fortran libraries. climt providesfunctionality requiring model-specific assumptions, such as stateinitialization and grid configuration. climt's programming interfacedesigned to be easy to use and thus appealing to a wide audience. Model building, configuration and execution are performed through a Pythonscript (or Jupyter Notebook), enabling researchers to build an end-to-endPython-based pipeline along with popular Python data analysis andvisualization tools.

97 MATHEMATICS AND COMPUTING↗

Forcing data (CESM2/CMIP6) for projection of drought impacts (2015-2100) at the K34 site in Manaus, Brazil

Historical and projected output data variables extracted and derived from Community Earth System Model 2 (CESM2) runs from the Coupled Model Intercomparison Project Phase 6 (CMIP6) archive. CESM2 is a fully coupled Earth system model used in simulations of Earth's past, present, and future climates (Danabasoglu et al., 2020). Variables in this dataset are in six-hourly resolution, and include air temperature (both in K and ℃), specific humidity (kg kg-1), air pressure (Pa), relative humidity (%), and vapor pressure deficit (kPa). CESM2 was the only CMIP6 model that provided VPD at the high temporal resolution required for this analysis. Data are included in .csv files, and the text file CESM2-CMIP6_forcing_K34-Manaus_headers.txt provides descriptions of data file headers.

54 ENVIRONMENTAL SCIENCES↗

E3SM‐Arctic: Regionally Refined Coupled Model for Advanced Understanding of Arctic Systems Interactions

Earth system models are essential tools for climate projections, but coarse resolutions limit regional accuracy, especially in the Arctic. Regionally refined meshes (RRMs) enhance resolution in key areas while maintaining computational efficiency. This paper provides an overview of the United States (U.S.) Department of Energy's (DOE's) Energy Exascale Earth System Model version 2.1 with an Arctic RRM, hereafter referred to as E3SMv2.1-Arctic, for the atmosphere (25 km), land (25 km), and ocean/ice (10 km) components. We evaluate the atmospheric component and its interactions with land, ocean, and cryosphere by comparing the RRM (E3SM2.1-Arctic) historical simulations (1950–2014) with the uniform low-resolution (LR) counterpart, reanalysis products, and observational data sets. The RRM generally reduces biases in the LR model, improving simulations of Arctic large-scale mean fields, such as precipitation, atmospheric circulation, clouds, atmospheric river frequency, and sea ice thickness. However, it introduces a seasonally dependent surface air temperature bias, reducing the LR cold bias in summer but enhancing the LR warm bias in winter, which contributes to the underestimated winter sea ice area and volume. Radiative feedback analysis shows similar climate feedback strengths in both model configurations, with the RRM exhibiting a more positive surface albedo feedback and contributing to a stronger surface warming than LR. These findings underscore the importance of high-resolution modeling for advancing our understanding of Arctic climate changes and their broader global impacts, although some persistent biases appear to be independent of model resolution at 10–100 km scales.

Energy Exascale Earth System Model (E3SM)↗

An evaluation of the ocean and sea ice climate of E3SM using MPAS and interannual CORE-II forcing

The Energy Exascale Earth System Model (E3SM) is a new coupled Earth system model sponsored by the US Department of Energy. Here we present E3SM global simulations using active ocean and sea ice that are driven by the CORE–II inter–annual atmospheric forcing data set. The E3SM ocean and sea–ice components are MPAS–Ocean and MPAS–Seaice, which use the Model for Prediction Across Scales (MPAS) framework and run on unstructured horizontal meshes. For this study, grid cells vary from 30 to 60 km for the low resolution mesh and 6 to 18 km at high resolution. The vertical grid is a structured z–star coordinate and uses 60 and 80 layers for low and high resolution, respectively. The lower resolution simulation was run for five CORE cycles (310 years) with little drift in sea surface temperature or heat content. The meridional heat transport is within observational range, while the meridional overturning circulation at 26.5° N is low compared to observations. The largest temperature biases occur in the Labrador Sea and western boundary currents, and the mixed layer is deeper than observations at northern high latitudes in the winter months. In the Antarctic, maximum mixed layer depths (MLD) compare well with observations, but the spatial MLD pattern is shifted relative to observations. Sea–ice extent, volume and concentration agree well with observations. At high resolution, the sea surface height compares well with satellite observations in mean and variability.

54 ENVIRONMENTAL SCIENCES↗

The International Land Model Benchmarking (ILAMB) System: Design, Theory, and Implementation

Abstract The increasing complexity of Earth system models has inspired efforts to quantitatively assess model fidelity through rigorous comparison with best available measurements and observational data products. Earth system models exhibit a high degree of spread in predictions of land biogeochemistry, biogeophysics, and hydrology, which are sensitive to forcing from other model components. Based on insights from prior land model evaluation studies and community workshops, the authors developed an open source model benchmarking software package that generates graphical diagnostics and scores model performance in support of the International Land Model Benchmarking (ILAMB) project. Employing a suite of in situ, remote sensing, and reanalysis data sets, the ILAMB package performs comprehensive model assessment across a wide range of land variables and generates a hierarchical set of web pages containing statistical analyses and figures designed to provide the user insights into strengths and weaknesses of multiple models or model versions. Described here is the benchmarking philosophy and mathematical methodology embodied in the most recent implementation of the ILAMB package. Comparison methods unique to a few specific data sets are presented, and guidelines for configuring an ILAMB analysis and interpreting resulting model performance scores are discussed. ILAMB is being adopted by modeling teams and centers during model development and for model intercomparison projects, and community engagement is sought for extending evaluation metrics and adding new observational data sets to the benchmarking framework.

54 ENVIRONMENTAL SCIENCES↗

OCEANFILMS (Organic Compounds from Ecosystems to Aerosols: Natural Films and Interfaces via Langmuir Molecular Surfactants) sea spray organic aerosol emissions – implementation in a global climate model and impacts on clouds

Sea spray aerosol is one of the major sources of atmospheric particulate matter globally. It has increasingly been recognized that organic matter derived from ocean biological precursors contributes significantly to the composition of submicron sea spray and may modify sea spray aerosol impacts on clouds and climate. This paper describes the implementation of the OCEANFILMS (Organic Compounds from Ecosystems to Aerosols: Natural Films and Interfaces via Langmuir Molecular Surfactants) parameterization for sea spray organic aerosol emissions in a global Earth system model, the Energy Exascale Earth System Model (E3SM). OCEANFILMS is a physically based model that links sea spray chemistry with ocean biogeochemistry using a Langmuir partitioning approach. We describe the implementation details of OCEANFILMS within E3SM, compare simulated aerosol fields with observations, and investigate impacts on simulated clouds and climate. Four sensitivity cases are tested, in which organic emissions either strictly add to or strictly replace sea salt emissions (in mass and number) and are either fully internally or fully externally mixed with sea salt. The simulation with internally mixed, added organics agrees reasonably well with observed seasonal cycles of organic matter in marine aerosol and has been selected as the default configuration of the E3SM. In this configuration, sea spray organic aerosol contributes an additional source of cloud condensation nuclei, adding up to 30 cm -3 to Southern Ocean boundary-layer cloud condensation nuclei concentrations (supersaturation = 0.1 %). The addition of this new aerosol source strengthens shortwave radiative cooling by clouds by -0.36 W m -2 in the global annual mean and contributes more than -3.5 W m -2 to summertime zonal mean cloud forcing in the Southern Ocean, with maximum zonal mean impacts of about -4 W m -2 around 50–60°S. This is consistent with a previous top-down, satellite-based empirical estimate of the radiative forcing by sea spray organic aerosol over the Southern Ocean. Through its mechanistic approach, OCEANFILMS offers a path towards improved understanding of the feedbacks between ocean biology, sea spray organic matter, and climate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗