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 1,351 records · Page 75

Evolution of Ocean Heat Content Related to ENSO

As the strongest interannual perturbation to the climate system, El Niño–Southern Oscillation (ENSO) dominates the year-to-year variability of the ocean energy budget. Here we combine ocean observations, reanalyses, and surface flux data with Earth system model simulations to obtain estimates of the different terms affecting the redistribution of energy in the Earth system during ENSO events, including exchanges between ocean and atmosphere and among different ocean basins, and lateral and vertical rearrangements. This comprehensive inventory allows better understanding of the regional and global evolution of ocean heat related to ENSO and provides observational metrics to benchmark performance of climate models. Results confirm that there is a strong negative ocean heat content tendency (OHCT) in the tropical Pacific Ocean during El Niño, mainly through enhanced air–sea heat fluxes Q into the atmosphere driven by high sea surface temperatures. In addition to this diabatic component, there is an adiabatic redistribution of heat both laterally and vertically (0–100 and 100–300 m) in the tropical Pacific and Indian oceans that dominates the local OHCT. Heat is also transported and discharged from 20°S–5°N into off-equatorial regions within 5°–20°N during and after El Niño. OHCT and Q changes outside the tropical Pacific Ocean indicate the ENSO-driven atmospheric teleconnections and changes of ocean heat transport (i.e., Indonesian Throughflow). The tropical Atlantic and Indian Oceans warm during El Niño, partly offsetting the tropical Pacific cooling for the tropical oceans as a whole. While there are distinct regional OHCT changes, many compensate each other, resulting in a weak but robust net global ocean cooling during and after El Niño.

54 ENVIRONMENTAL SCIENCES↗

Less Atmospheric Radiative Heating By Dust Due to the Synergy of Coarser Size and Aspherical Shape

Mineral dust aerosols cool and warm the atmosphere by scattering and absorbing solar (shortwave: SW) and thermal (longwave: LW) radiation. However, significant uncertainties remain in dust radiative effects, largely due to differences in the dust size distribution and spectral optical properties simulated in Earth system models. Dust models typically underestimate the coarse dust load (more than 2.5 µm in diameter) and assume a spherical shape, which leads to an overestimate of the fine dust load (less than 2.5 µm) after the dust emissions in the models are scaled to match observed dust aerosol optical depth at 550 nm (DAOD550). Here, we improve the simulated dust properties with data sets that leverage measurements of size-resolved dust concentration, asphericity factor, and refractive index in a coupled global chemical transport model with a radiative transfer module. After the adjustment of size-resolved dust concentration and spectral optical properties, the global and annual average of DAOD550 from the simulation increases from 0.023 to 0.029 and falls within the range of a semi-observationally based estimate (0.030 ± 0.005). The reduction of fine dust load after the adjustment leads to a reduction of the SW cooling at the top of the atmosphere (TOA). To improve agreement against a semi-observationally based estimate of the radiative effect efficiency at TOA, we find that a less absorptive SW dust refractive index is required for coarser aspherical dust. Thus, only a minor difference is estimated for the net global dust radiative effect at TOA (−0.08 vs. −0.00 W m−2 on a global scale). Conversely, our sensitivity simulations reveal that the surface warming is substantially enhanced near the strong dust source regions (less cooling to −0.23 from −0.60 W m−2 on a global scale). Thus, less atmospheric radiative heating is estimated near the major source regions (less heating to 0.15 from 0.59 W m−2 on a global scale), because of enhanced LW warming at the surface by the synergy of coarser size and aspherical shape.

simulated dust properties↗

Changes in Sea Salt Emissions Enhance ENSO Variability

Two 150-year pre-industrial simulations with and without interactive sea salt emissions from the Community Earth System Model (CESM) are performed to quantify the interactions between sea salt emissions and El Niño–Southern Oscillation (ENSO). Variations in sea salt emissions over the tropical Pacific Ocean are affected by changing wind speed associated with ENSO variability. ENSO-induced interannual variations in sea salt emissions result in decreasing (increasing) aerosol optical depth (AOD) by 0.03 over the equatorial central-eastern (western) Pacific Ocean during El Niño events compared to those during La Niña events. These changes in AOD further increase (decrease) radiative fluxes into the atmosphere by +0.2 W m -2 (-0.4 W m -2 ) over the tropical eastern (western) Pacific. Thereby, sea surface temperature increases (decreases) by 0.2–0.4 K over the tropical eastern (western) Pacific Ocean during El Niño compared to La Niña events and enhances ENSO variability by 10%. The increase in ENSO amplitude is a result of systematic heating (cooling) during the warm (cold) phase, of ENSO in the eastern Pacific. Interannual variations in sea salt emissions then produce the anomalous ascent (subsidence) over the equatorial eastern (western) Pacific between El Niño and La Niña events, which is a result of heating anomalies. Due to variations in sea salt emissions, the convective precipitation is enhanced by 0.6–1.2 mm day -1 over the tropical central-eastern Pacific Ocean and weakened by 0.9–1.5 mm day -1 over the Maritime Continent during El Niño compared to La Niña events, enhancing the precipitation variability over the tropical Pacific.

58 GEOSCIENCES↗

Significant Rural Contributions to Severe Coastal Urban Flooding Through Flow Connectivity

Coastal cities are increasingly threatened by hurricane-induced compound flooding. However, the relative contributions of rainfall-surge interactions and other flood drivers, such as runoff from surrounding rural and local impervious areas, in shaping urban flooding remain unclear. Here, we use a multi-scale Earth system modeling framework to disentangle the contributions of diverse flooding drivers. The results show that urban flooding is primarily affected by topographic factors (low elevation, flat slopes, and localized flow constrictions), rather than increased runoff from their impervious surfaces. Importantly, runoff from rural areas contributes more to severely flooded urban areas with inundation depth exceeding 1.7 m than local runoff. Although rainfall–surge interactions and sea-level rise exacerbate inundation along the coastline, over 90% of the compounding effects are absorbed by wetlands. By resolving flooding processes across scales and components, this study emphasizes the role of runoff from surrounding rural areas in shaping zones of severe urban flooding.

58 GEOSCIENCES↗

Size-resolved process understanding of stratospheric sulfate aerosol following the Pinatubo eruption

Stratospheric sulfate aerosol produced by volcanic eruptions plays important roles in atmospheric chemistry and the global radiative balance of the atmosphere. The simulation of stratospheric sulfate concentrations and optical properties is highly dependent on the chemistry scheme and microphysical treatment. In this work, we implemented a sophisticated gas-phase chemistry scheme (full chemistry, FC) and a 5-mode version of the Modal Aerosol Module with Prognostic Stratospheric Aerosol (MAM5-PSA) for the interactive treatment of stratospheric sulfate aerosol in the Department of Energy's Energy Exascale Earth System Model version 2 (E3SMv2) model to better simulate the chemistry-aerosol feedback following the Pinatubo eruption, and to compare it against a simulation using simplified chemistry (SC) and the default 4-mode version of the Modal Aerosol Module (MAM4). MAM5-PSA experiments were found to better capture the stratospheric sulfate burden from the eruption of the volcano to the end of 1992 as compared to the High-resolution Infrared Sounder (HIRS) observations, and the formation of sulfate in MAM5-PSA with FC (with an additional OH replenishment reaction) was significantly faster than in MAM4 with FC. Analyses of microphysical processes indicate that more sulfate aerosol mass was generated in total in FC experiments than in SC experiments. MAM5-PSA performs better than MAM4 in simulation of aerosol optical depth (AOD); AOD anomalies from the MAM5-PSA experiment have better agreement with observations. The simulated largest changes in global mean net radiative flux at the top of the atmosphere following the eruption were about −3 W m −2 in MAM5-PSA experiments and roughly −1.5 W m −2 in MAM4 experiments.

AEROSOL↗

The Freundlich isotherm equation best represents phosphate sorption across soil orders and land use types in tropical soils of Puerto Rico

Biomass production in the lowland wet tropical forest is greater than in any other biome, and it is typically limited by soil phosphorus (P) availability. However, the mechanisms involved in the P cycle remain poorly represented in Earth System Models (ESMs). Soil P sorption processes are key in the P cycle and for understanding the extent of P limitation for plant productivity. Currently, a few ESMs include isotherm equations to model these processes. Although the Langmuir equation is widely cited, other isotherm equations may better describe sorption in tropical soils. Here, we use a diverse range of soil samples from Puerto Rico to test the validity of the Langmuir, Freundlich, and Temkin equation. We found that across four soil orders (Inceptisols, Mollisols, Oxisols, Ultisols), and forested and cultivated land use types, the Freundlich equation best represented soil P sorption. Furthermore, the Langmuir and the Temkin equations poorly represent soil P adsorption, especially at low P concentrations. Specifically, the Langmuir equation underestimated soil P adsorption by 40% and the Temkin equation overestimated adsorption by 76%. We also found, as expected, that soil clay content and pH were the most important parameters explaining the variability of the Freundlich (K f ) constant. Greater clay content and lower pH, common in highly weathered Ultisols and Oxisols which are abundant in the tropics, led to greater K f values. Overall, our results suggest that a diversity of soils can prompt underestimation of P sorption when using the Langmuir isotherm, which leads to an overestimation of available P that can have repercussions on ESM predictions of the P cycle and tropical forest productivity.

Earth system models↗

NCAR Contribution to A U.S. National Multi-Model Ensemble (NMME) ISI Prediction System

NCAR brought the latest version of the Community Earth System Model (version 1, CESM1) into the mix of models in the NMME effort. This new version uses our newest atmospheric model CAM5 and produces a coupled climate and ENSO that are generally as good or better than those of the Community Climate System Model version 4 (CCSM4). Compared to CCSM4, the new coupled model has a superior climate response with respect to low clouds in both the subtropical stratus regimes and the Arctic. However, CESM1 has been run to date using a prognostic aerosol model that more than doubles its computational cost. We are currently evaluating a version of the new model using prescribed aerosols and expect it will be ready for integrations in summer 2012. Because of this NCAR has not been able to complete the hindcast integrations using the NCAR loosely-coupled ensemble Kalman filter assimilation method nor has it contributed to the current (Stage I) NMME operational utilization. The expectation is that this model will be included in the NMME in late 2012 or early 2013. The initialization method will utilize the Ensemble Kalman Filter Assimilation methods developed at NCAR using the Data Assimilation Research Testbed (DART) in conjunction with Jeff Anderson’s team in CISL. This methodology has been used in our decadal prediction contributions to CMIP5. During the course of this project, NCAR has setup and performed all the needed hindcast and forecast simulations and provide the requested fields to our collaborators. In addition, NCAR researchers have participated fully in research themes (i) and (ii). Specifically, i) we have begun to evaluate and optimize our system in hindcast mode, focusing on the optimal number of ensemble members, methodologies to recalibrate individual dynamical models, and accessing our forecasts across multiple time scales, i.e., beyond two weeks, and ii) we have begun investigation of the role of different ocean initial conditions in seasonal forecasts. The completion of the calibration hindcasts for Seasonal to Interannual (SI) predictions and the maintenance of the data archive associated with the NCAR portion of this effort has been the responsibility of the Project Scientist I (Alicia Karspeck) that was partially supported on this project.

54 ENVIRONMENTAL SCIENCES↗

Computational Reproducibility of Scientific Workflows at Extreme Scales

We propose an approach for improved reproducibility that includes capturing and relating provenance characteristics and performance metrics, in a hybrid queriable system, the ProvEn server. The system capabilities are illustrated on two use cases: scientific reproducibility of results in the ACME climate simulations and performance reproducibility in molecular dynamics workflows on HPC computing platforms.

provenance, Reproducibility, Climate & Earth Syste↗

A multivariate approach to ensure statistical reproducibility of climate model simulations

Effective utilization of novel hybrid architectures of pre-exascale and exascale machines requires transformations to global climate modeling systems that may not reproduce the original model solution bit-for-bit. Round-off level differences grow rapidly in these non-linear and chaotic systems. This makes it difficult to isolate bugs/errors from innocuous growth expected from round-off level differences. Here, we apply two modern multivariate two sample equality of distribution tests to evaluate statistical reproducibility of global climate model simulations using standard monthly output of short (~ 1-year) simulation ensembles of a control model and a modified model of US Department of Energy's Energy Exascale Earth System Model (E3SM). Both the tests are able to identify changes induced by modifications to some model tuning parameters. We also conduct formal power analyses of the tests by applying them on designed suites of short simulation ensembles each with an increasingly different climate from the control ensemble. The results are compared against those from another such test. These power analyses provide a framework to quantify the degree of differences that can be detected confidently by the tests for a given ensemble size (sample size). This will allow model developers using the tests to make an informed decision when accepting/rejecting an unintentional non-bit-for-bit change to the model solution.

Mahajan, Salil↗

Improving the Parameterization of Cloud and Rain Microphysics in E3SM using Novel Observationally-Constrained Bayesian Approach (Final Technical Report)

In this project, we sought to develop new cloud and rain microphysics frameworks within the Energy Exascale Earth System Model (E3SM). This work encompassed two primary avenues of research: 1) Further development of a Bayesian-based scheme called BOSS (Bayesian Observationally-constrained Statistical-physical Scheme) to represent cloud and rain microphysics, testing it in realistic high-resolution cloud models, and implementing it in E3SM; 2) Development of a methodology utilizing machine learning to enable computationally tractable use of tractable use of Markov chain Monte Carlo sampling for Bayesian parameter estimation in Earth system and cloud models. In this project, we adapted the BOSS microphysics scheme, originally formulated for rain-only, to include all liquid-phase microphysical processes for cloud and rain, in particular the processes that mediate between these two categories, for example the conversion from cloud to rain through collision and coalescence of drops. We constrained the scheme via comparison and testing against a detailed model that explicitly represents the evolution of cloud and rain particles, called a bin microphysics scheme.

54 ENVIRONMENTAL SCIENCES↗

A Framework for Deep Learning Emulation of Numerical Models With a Case Study in Satellite Remote Sensing

Numerical models based on physics represent the state of the art in Earth system modeling and comprise our best tools for generating insights and predictions. Despite rapid growth in computational power, the perceived need for higher model resolutions overwhelms the latest generation computers, reducing the ability of modelers to generate simulations for understanding parameter sensitivities and characterizing variability and uncertainty. Thus, surrogate models are often developed to capture the essential attributes of the full-blown numerical models. Recent successes of machine learning methods, especially deep learning (DL), across many disciplines offer the possibility that complex nonlinear connectionist representations may be able to capture the underlying complex structures and nonlinear processes in Earth systems. A difficult test for DL-based emulation, which refers to function approximation of numerical models, is to understand whether they can be comparable to traditional forms of surrogate models in terms of computational efficiency while simultaneously reproducing model results in a credible manner. A DL emulation that passes this test may be expected to perform even better than simple models with respect to capturing complex processes and spatiotemporal dependencies. Here, we examine, with a case study in satellite-based remote sensing, the hypothesis that DL approaches can credibly represent the simulations from a surrogate model with comparable computational efficiency. Our results are encouraging in that the DL emulation reproduces the results with acceptable accuracy and often even faster performance. We discuss the broader implications of our results in light of the pace of improvements in high-performance implementations of DL and the growing desire for higher resolution simulations in the Earth sciences.

Bayesian Deep Learning↗

The Race To Understand A Changing Planet

The Earth's climate is changing rapidly. In some respects, the rate of change is outpacing the predictions of only a few years ago. The challenge to Earth Science is to put forward credible projections of possible future climates so that the public and policy makers can make science-based decisions about energy development strategies. Models, observations and experiments all play strong roles in improving knowledge and increasing confidence in our predictions. The models have progressed from simple, coarse-resolution descriptions of atmospheric dynamics and physics only twenty years ago, to full-up Earth System models (ESMs) that include complete descriptions of the oceans and cryosphere. It has been convincingly argued that such complexity - the construction of realistic "toy" Earth's - is necessary to address the complex processes involved in climate change, including not only the physical atmosphere, oceans and cryosphere, but also the carbon cycle - both its natural and anthropogenic components - and the biosphere. Observations, particularly satellite observations, have more or less kept pace with the demands of the modelers, being able to observe progressively more and different facets of the Earth system, but the global satellite fleet is in need of an overhaul very soon. Lastly, field experiments and process studies confront the models with facts and allow us to develop more sophisticated and accurate satellite data algorithms. The challenges facing our relatively small Earth and planetary science communities are considerable and the stakes are significant. The stakeholders, now numbering 7 billion but soon to be 10 billion, will be relying on our results and capabilitie's to guide them into the future.

Sellers, Piers J.↗

Composition and physical properties of the Asian Tropopause Aerosol Layer and the North American Tropospheric Aerosol Layer: Composition of ATAL and NATAL

Recent studies revealed layers of enhanced aerosol scattering in the upper troposphere and lower stratosphere over Asia (Asian Tropopause Aerosol Layer (ATAL)) and North America (North American Tropospheric Aerosol Layer (NATAL)). We use a sectional aerosol model (Community Aerosol and Radiation Model for Atmospheres (CARMA)) coupled with the Community Earth System Model version 1 (CESM1) to explore the composition and optical properties of these aerosol layers. The observed aerosol extinction enhancement is reproduced by CESM1/CARMA. Both model and observations indicate a strong gradient of the sulfur-to-carbon ratio from Europe to the Asia on constant pressure surfaces. We found that the ATAL is mostly composed of sulfates, surface-emitted organics, and secondary organics; the NATAL is mostly composed of sulfates and secondary organics. In conclusion, the model also suggests that emission increases in Asia between 2000 and 2010 led to an increase of aerosol optical depth of the ATAL by 0.002 on average which is consistent with observations.

54 ENVIRONMENTAL SCIENCES↗

Future projections of wind patterns in California with the variable-resolution CESM: a clustering analysis approach

Wind energy production is expected to be affected by shifts in wind patterns that will accompany climate change. However, many questions remain on the magnitude and character of this impact, especially on regional scales. In this study, clustering is used to group and analyze large-scale wind patterns in California using model simulations from the variable-resolution Community Earth System Model (VR-CESM). Specifically, simulations have been produced that cover historical (1980–2000), mid-century (2030–2050), and end-of-century (2080–2100) time periods. Once clustered, observed changes to wind patterns can be analyzed in terms of both the change in frequency of those clusters and changes to winds within-clusters. Statistically significant capacity factors changes have been found at all five wind plant sites. Further, decomposition of the capacity factor changes into frequency changes and within-cluster changes enables a better understanding of their drivers. A further examination of the synoptic-scale fields associated with each cluster then provides a better understanding of how changes to large-scale meteorological fields are important for driving changes in localized wind speeds.

17 WIND ENERGY↗

Region and cloud regime dependence of parametric sensitivity in E3SM atmosphere model

Abstract The Department of Energy (DOE)’s Energy Exascale Earth System Model (E3SM), including its atmosphere model (EAM), has many relatively new features. In a previous study we conducted a systematic parametric sensitivity analysis for EAM based on short, perturbed parameter ensemble (PPE) simulations, mainly focusing on global mean climate features and metrics. While parameter values in global climate models are generally invariant in space and time, model response to parameters perturbation may vary by regions and climate regimes, which motivates the need to better understand the EAM model behaviors and physics at regional scale and process level. In this study, using the same set of PPE simulations and a similar sensitivity analysis framework, we identify parameters that cause largest sensitivities over different regions and compare model responses in fast atmospheric processes to the parameters across different cloud regimes for several important cloud-related fidelity metrics. We find that cloud forcing has opposite response to some parameters over mid-latitude vs. tropical land. We also analyze how the parametric sensitivity varies as stratocumulus transitions to shallow convection and to deep convection over ocean. Low cloud forcing and shortwave cloud forcing in the subtropical eastern Pacific are most sensitive to the parameters controlling the width of the probability density function (PDF) of the subgrid vertical velocity ( w’ ) ( gamma ) and the damping of the w’ skewness ( c8 ) near the coast but become more sensitive to the parameter affecting the damping of the w’ variance ( c1 ) further offshore. Detailed interpretation of the spatial dependence of parametric sensitivity is provided. We also investigate how the parametric sensitivity evolves with prediction duration. This study improves our process-level understanding of cloud physics and parameterization and provides insights for developing more advanced regime-aware parameterization schemes in global climate model.

54 ENVIRONMENTAL SCIENCES↗

Advancing quantitative models of soil microbiology, ecology, and biochemistry

Rapid advances in soil biology are increasingly reflected in mathematical models. From molecular to global scales, models contribute to fundamental understanding and prediction of critical soil processes, such as carbon sequestration. Recent models explicitly incorporate microbial control over soil carbon and nutrient cycling, an approach that offers greater biological realism but raises new challenges in model parameterization and stability. Although they have proliferated at community and ecosystem scales, microbial-explicit models have yet to be incorporated into most Earth system models. Applying soil models across scales requires integration with data to support parameterization and validation using approaches such as Bayesian data assimilation. In conclusion, to catalyze future progress, we recommend breaking down barriers between modeling and empirical disciplines while broadening access to the computational infrastructure that supports soil biological modeling.

Xie, Wally↗

Evaluation of an Improved Convective Triggering Function: Observational Evidence and SCM Tests

This study provides a strong observational support for a recently developed convective triggering function that uses the large-scale dynamic convective available potential energy (dCAPE) as a constraint combined with an unrestricted air parcel launch level (ULL) to relax the unrealistic strong coupling of convection to surface heating and capture nocturnal elevated convection. Both case study and statistical analysis are conducted using the observations collected from the Department of Energy (DOE)'s Atmospheric Radiation Measurement (ARM) program at its Southern Great Plains (SGP) and Manaus (MAO) sites. They show that dCAPE has a much stronger correlation with precipitation than CAPE and ULL is essential to detect elevated convection above boundary layer under both midlatitude and tropical conditions and for both afternoon and nighttime deep convection regimes. Sensitivity tests with the single-column model (SCM) of DOE's Energy Exascale and Earth System Model (E3SM) indicate that the role of dCAPE in suppressing daytime convection is more effective for tropical convection than midlatitude convection. Even though the dCAPE can suppress the overestimated convection, ULL plays a much bigger role in improving the diurnal cycle of precipitation than dCAPE. It not only helps capture nocturnal elevated convection, but also significantly removes the spurious morning precipitation seen in the default model, due to the release of unstable energy at night. However, the use of ULL has led to an overestimation of light-to-moderate precipitation (1-10 mm/day) due to more convection being triggered above boundary layer.

54 ENVIRONMENTAL SCIENCES↗

Vertically Resolved Analysis of the Madden-Julian Oscillation Highlights the Role of Convective Transport of Moist Static Energy

We simulate the Madden-Julian oscillation (MJO) over an aquaplanet with uniform surface temperature using the multiscale modeling framework (MMF) configuration of the Energy Exascale Earth System Model (E3SM-MMF). The model produces MJO-like features that have a similar spatial structure and propagation behavior to the observed MJO. To explore the processes involved in the propagation and maintenance of these MJO-like features, we perform a vertically resolved moist static energy (MSE) analysis for the MJO (Yao et al., 2022, https://doi.org/10.1175/jas-d-20-0254.1). Unlike the column-integrated MSE analysis, our method emphasizes the local production of MSE variance and quantifies how individual physical processes amplify and propagate the MJO's characteristic vertical structure. We find that radiation, convection, and boundary layer (BL) processes all contribute to maintaining the MJO, balanced by the large-scale MSE transport. Furthermore, large-scale dynamics, convection, and BL processes all contribute to the propagation of the MJO, while radiation slows the propagation. Additionally, we perform mechanism-denial experiments to examine the role of radiation and associated feedbacks in simulating the MJO. We find that the MJO can still self-emerge and maintain its characteristic structures without radiative feedbacks. This study highlights the role of convective MSE transport in the MJO dynamics, which was overlooked in the column-integrated MSE analysis.

54 ENVIRONMENTAL SCIENCES↗