SEARCH · Engineering Papers

Results for “Earth System Model”

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 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Using a surrogate-assisted Bayesian framework to calibrate the runoff-generation scheme in the Energy Exascale Earth System Model (E3SM) v1

Abstract. Runoff is a critical component of the terrestrial water cycle, and Earth system models (ESMs) are essential tools to study its spatiotemporal variability. Runoff schemes in ESMs typically include many parameters so that model calibration is necessary to improve the accuracy of simulated runoff. However, runoff calibration at a global scale is challenging because of the high computational cost and the lack of reliable observational datasets. In this study, we calibrated 11 runoff relevant parameters in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) using a surrogate-assisted Bayesian framework. First, the polynomial chaos expansion machinery with Bayesian compressed sensing is used to construct computationally inexpensive surrogate models for ELM-simulated runoff at 0.5∘ × 0.5∘ for 1991–2010. The error metric between the ELM simulations and the benchmark data is selected to construct the surrogates, which facilitates efficient calibration and avoids the more conventional, but challenging, construction of high-dimensional surrogates for the ELM simulated runoff. Second, the Sobol' index sensitivity analysis is performed using the surrogate models to identify the most sensitive parameters, and our results show that, in most regions, ELM-simulated runoff is strongly sensitive to 3 of the 11 uncertain parameters. Third, a Bayesian method is used to infer the optimal values of the most sensitive parameters using an observation-based global runoff dataset as the benchmark. Our results show that model performance is significantly improved with the inferred parameter values. Although the parametric uncertainty of simulated runoff is reduced after the parameter inference, it remains comparable to the multimodel ensemble uncertainty represented by the global hydrological models in ISMIP2a. Additionally, the annual global runoff trend during the simulation period is not well constrained by the inferred parameter values, suggesting the importance of including parametric uncertainty in future runoff projections.

58 GEOSCIENCES↗

DOE OSTI2024

The atmospheric effect of aerosols on future tropical cyclone frequency and precipitation in the Energy Exascale Earth System Model

This study uses experiments from the Energy Exascale Earth System Model (E3SM) to compare the influence on tropical cyclone (TC) activity of: (i) the atmospheric effect of aerosols under specified sea-surface temperatures (SSTs); and (ii) the net effect of greenhouse gases (GhGs) (including changes in SSTs). The experiments were performed using the CMIP6 Shared Socioeconomic Pathway SSP5-8.5 emissions scenario with GhG-induced SST warming specified and atmospheric aerosol effects simulated but without explicit ocean coupling. Insignificant changes in global TC frequency are found in response to the atmospheric effect of future aerosols and GhGs, as significant regional responses in TC frequency counteract each other. Future GhGs contribute to more frequent TCs in the North Atlantic, and reductions over the Northwestern Pacific and Southern Indian Ocean. The atmospheric effect of future aerosols drives more frequent TCs over the Northwestern Pacific and reductions over the Northeast Pacific and North Atlantic. Along with increases in TC intensity, global TC precipitation (TCP) is projected to increase by 52.8% (14.1%/K) due to the combined effect of future aerosols and GhGs. Although both forcings contribute to TCP increases (14.7–19.3% from reduced aerosols alone and 28.1–33.3% from increased GhGs alone), they lead to different responses in the spatial structure of TCP. TCP increases preferentially in the inner-core due to increased GhGs, whereas TCP decreases in the inner-core and increases in the outer-bands in response to the atmospheric effects of decreased aerosols. These changes are distinct from those caused by aerosol-induced SST changes, which have been considered in other studies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Is land use producing robust signals in future projections from Earth system models, all else being equal?

Abstract We use six Earth system models (ESMs) run under SSP3-7.0, a scenario characterized by a relatively large land use change (LUC) over the 21st century, and under a variant of the same scenario where a significantly different pattern of LUC, taken from SSP1-2.6, was used, all else being equal. Our goal is to identify changes in climate extremes between the two scenarios that are statistically significant and robust across the ESMs. The motivation for this study is to test a long-held assumption of the shared socio-economic pathway-representative concentration pathway (SSP-RCP) scenario framework: that the signal from LUC can be safely disregarded when pairing different SSPs to the compatible RCPs, where compatibility only considers global radiative forcing, predominantly determined by well-mixed greenhouse gasses emissions. We analyze extremes of daily minimum and maximum temperatures and precipitation, after fitting non-stationary generalized extreme value distributions in a way that borrows strength along the length of the simulation (2015–2100) and across initial condition ensembles. We consider changes in the 20 year return levels (RL20s) of these metrics by 2100, and focus on eight locations where LUC is large within each scenario, and strongly differs between scenarios, averaging the RL20s over a neighborhood characterized by the same LUC to enhance the signal to noise. We find that precipitation extremes do not show significant differences attributable to LUC differences. For temperature extremes (cold and hot) results are mixed, with some location-index combination showing significant results for some of the ESMs but not all, and not many coherent changes appearing for indices across regions, or regions across indices. These ESMs are representative of what is typically adopted as the source of climate information for impact studies, when the SSP-RCP framework is put to use. Overall, our analysis suggests that the hypothesis to pair SSPs to RCPs in a flexible fashion is overall defensible. However, the appearance of some coherence in a few locations and for some indices invites further investigation.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRSDate not supplied

Ocean and Earth System Modelling

Petascale supercomputing infrastructure + modelling and analysis capabilities + interdisciplinary upper-ocean expertise Multiscale ocean turbulence simulation Physical-biogeochemical interactions Analysis of large data sets from remote sensing and Earth system model ensembles

Ocean↗

DOE OSTI2021

Projecting climate change in South America using variable‐resolution Community Earth System Model: An application to Chile

Abstract We introduce variable‐resolution enabled Community Earth System Model (VR‐CESM) results simulating historical and future climate conditions at 28 km over South America and 14 km over the Andes. Three 30‐year simulations are performed: a historic (1985–2014), a near future (2030–2059), and an end‐century (2070–2099) simulation under the RCP8.5 scenario. Historic results compare favourably to several temperature and precipitation reanalysis products, though local biases are present, particularly during austral summer. Future simulations highlight broad warming patterns (+3–6°C by end‐century) and heterogeneous precipitation responses across South America that qualitatively agree with prior modelling efforts. Our results reveal that the interaction between temperature and precipitation changes produce shifts in several Köppen–Geiger climates. Notable changes include the near‐elimination of the Andean Tundra or Alpine climates, a 15% decrease in Tropical Rainforests and a Tropical Savannah expansion of 20%. To provide a regionally focused analysis of projected climate change and to illustrate the benefits of variable resolution modelling, we analyse changes in the magnitude and trend in seasonal and daily temperature and precipitation in Chile. We also examined several metrics [e.g., snow water equivalent (SWE), temperatures on wet days, and days below 0°C] to evaluate potential impacts of climate change on the Chilean cryosphere between the end‐of‐century and historic periods, finding wide‐ranging indications of cryospheric decline. These changes are interpreted through reductions in the timing (1–2.5 months earlier peak SWE) and magnitude (200–1,000 mm SWE decreases) of water stored as snow in the Andes, a 10–30% decrease in number of cool season wet days with temperatures below 1°C, and 50–200 fewer days (annually) with minimum temperatures below 0°C. Our aim in producing a high‐resolution dataset of climate projections from VR‐CESM is to support analyses of climate change throughout South America but especially in vulnerable montane regions and to provide additional results for comparison with previous, ongoing, and upcoming modelling efforts.

variable-resolution global climate models↗

DOE OSTI2026

Importance of Spatially Continuous Urban Surface Properties in Urban‐Resolving Earth System Modeling

Accurate representation of urban properties and processes at higher resolutions in global modeling systems is essential for advancing our ability to capture the complexities of urban systems and informing effective resilience strategies. However, the prescription of coarse global-scale urban properties in most state-of-the-art Earth system models (ESMs) is limiting their potential for capturing urban signals as they advance toward kilometer-scale simulation capabilities. To bridge this gap in inadequate urban property representation and to advance urban-resolving Earth system modeling, this work integrates the newly-developed global 1 km-resolution facet-level urban surface property data set, U-Surf, into the land component of Community Earth System Model (CESM)—Community Terrestrial System Model (CTSM). The land-only CTSM simulations are validated against satellite measurements, ground-based urban weather stations, flux tower observations, and reanalysis data. Results demonstrate that the enhanced urban properties allow improved simulations of urban meteorology and surface energy fluxes compared to the default coarse-resolution categorical urban canopy parameters. Spatial scaling analysis reveals regime-dependent information loss during resolution aggregation, as well as substantial scale-dependent variations in urban surface energy flux representation. Furthermore, these findings have critical implications for coupled Earth system modeling when including the effect of land-atmosphere interaction. This work establishes a foundation for future urban-resolving kilometer-scale ESM development, which will enable systematic intra- and inter-city comparisons that inform urban adaptation strategies across diverse global urban environments.

Cheng, Yifan [University of Illinois Urbana-Champa↗

DOE OSTI2023

Simulating Southern Ocean Aerosol and Ice Nucleating Particles in the Community Earth System Model Version 2

Southern Ocean (SO) low-level mixed phase clouds have been a long-standing challenge for Earth system models to accurately represent. While improvements to the Community Earth System Model version 2 (CESM2) resulted in increased supercooled liquid in SO clouds and improved model radiative biases, simulated SO clouds in CESM2 now contain too little ice. Previous observational studies have indicated that marine particles are major contributor to SO low-level cloud heterogeneous ice nucleation, a process that initiates a number of cloud processes that govern cloud radiative properties. In this study, we utilize detailed aerosol and ice nucleating particle (INP) measurements from two recent measurement campaigns to assess simulated aerosol abundance, number size distributions, and composition and INP parameterizations for use in CESM2. Our results indicate that CESM2 has a positive bias in simulated surface-level total aerosol surface area at latitudes north of 58°S. Measured INP populations were dominated by marine INPs and we present evidence of refractory INPs present over the SO assumed here to be mineral dust INPs. Results highlight a critical need to assess simulated mineral dust number and size distributions in CESM2 in order to adequately represent SO INP populations and their response to long-term changes in atmospheric transport patterns and land use change. Furthermore, we also discuss important cautions and limitations in applying a commonly used mineral dust INP parameterization to remote regions like the pristine SO.

58 GEOSCIENCES↗

DOE OSTI2022

Aircraft Observations of Turbulence in Cloudy and Cloud-Free Boundary Layers Over the Western North Atlantic Ocean From ACTIVATE and Implications for the Earth System Model Evaluation and Development

This study examines boundary layer turbulence derived from high temporal resolution meteorological measurements from 40 research flights over the western North Atlantic Ocean during the 2020 deployments of ACTIVATE. Frequency distributions of various turbulent quantities reveal stronger turbulence during the winter deployment than in summer and for cloud‐topped than in cloud‐free boundary layers during the summer deployment. Maximum turbulence kinetic energy (TKE) is most often within cloud from observations in winter and summer, whereas it is mostly below cloud in both seasons by a global model turbulence parameterization. Further, bivariate frequency distributions are consistent with the bivariate Gaussian probability distribution functions assumed for the closure of higher‐order turbulence/shallow convection parameterizations used by some global models. Turbulence simulated by the Community Atmosphere Model version 6 and the Energy Exascale Earth System Model Atmosphere Model version 2 using such parameterizations is not as strong as observed, with more TKE going into vertical wind perturbations rather than into zonal wind perturbations as observed, suggesting that the treatment of turbulence in Earth system models still needs to be further improved.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRSDate not supplied

Climate Change Impact on Surface Ozone Based on CMIP6 Earth System Models

It is presented an analysis of the effect of climate change on surface ozone (O3) discussing the related penalties and benefits around the globe from the global modeling perspective based on simulations with five CMIP6 (Coupled Model Intercomparison Project Phase 6) Earth System Models. All models conducted simulation experiments considering future climate (ssp370SST) and present-day climate (ssp370pdSST) under the same future emissions scenario (SSP3-7.0). Over regions remote from pollution sources, there is a robust decline in mean surface ozone concentration varying spatially from -0.2 to -2 ppbv oC-1, with strongest decline over tropical oceanic regions, which is mainly linked to the dominating role of enhanced ozone chemical loss with higher water vapour abudances under a warmer climate. However, ozone increases over regions close to anthropogenic pollution sources or close to enhanced natural Biogenic Volatile Organic Compounds (BVOC) emission sources with a rate ranging regionally from 0.2 to 2 ppbv oC-1, implying a regional surface ozone penalty due to global warming. The individual models show this robustly for south-eastern China and India as well as for regions of Africa but there are inter-model differences in areas within Europe and the United States (US) as well as in South America. The future climate change enhances the efficiency of precursor emissions to generate surface ozone in polluted regions and thus the magnitude of this effect depends on the regional emission changes considered in this study within the SSP3_7.0 scenario. The comparison of the climate change impact effect on surface ozone versus the combined effect of climate and emission changes indicates the dominant role of precursor emission changes in projecting surface ozone concentrations under future climate change scenarios.

climate change↗

DOE OSTI2022

The Seasonal-to-Multiyear Large Ensemble (SMYLE) prediction system using the Community Earth System Model version 2

Abstract. The potential for multiyear prediction of impactful Earth system change remains relatively underexplored compared to shorter (subseasonal to seasonal) and longer (decadal) timescales. In this study, we introduce a new initialized prediction system using the Community Earth System Model version 2 (CESM2) that is specifically designed to probe potential and actual prediction skill at lead times ranging from 1 month out to 2 years. The Seasonal-to-Multiyear Large Ensemble (SMYLE) consists of a collection of 2-year-long hindcast simulations, with four initializations per year from 1970 to 2019 and an ensemble size of 20. A full suite of output is available for exploring near-term predictability of all Earth system components represented in CESM2. We show that SMYLE skill for El Niño–Southern Oscillation is competitive with other prominent seasonal prediction systems, with correlations exceeding 0.5 beyond a lead time of 12 months. A broad overview of prediction skill reveals varying degrees of potential for useful multiyear predictions of seasonal anomalies in the atmosphere, ocean, land, and sea ice. The SMYLE dataset, experimental design, model, initial conditions, and associated analysis tools are all publicly available, providing a foundation for research on multiyear prediction of environmental change by the wider community.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Enhanced Simulation of Atmospheric Blocking in a High‐Resolution Earth System Model: Projected Changes and Implications for Extreme Weather Events

Atmospheric blocking is closely linked to the occurrence of extreme weather events. However, low-resolution Earth system models often underestimate the frequency of blocking, undermining confidence in future projections. Here, in this study, we use the high-resolution Community Earth System Model (CESM-HR; 25 km atm and 10 km ocean) to show that CESM-HR reduces biases in atmospheric blocking for both winter and summer, particularly for events lasting longer than 10 days. This improvement is partly due to reduced sea surface temperature biases at higher resolution. Additionally, applying a bias correction to the 500 hPa geopotential height further enhances blocking frequency simulations, highlighting the crucial role of the mean state. Under the Representative Concentration Pathway 8.5 scenario, CESM-HR projects a decrease in wintertime blocking over regions such as the Euro-Atlantic and Chukchi-Alaska, consistent with previous studies. In contrast, summer blocking is expected to become more frequent and persistent, driven by weakened zonal winds. The blocking center shifts from historical locations over Scandinavia and eastern Russia to central Eurasia, significantly increasing blocking over the Ural region. Summer blocking frequency over the Scandinavia-Ural region may eventually surpass historical winter blocking over the Euro-Atlantic. This increase in summer blocking could exacerbate summer heatwaves in a warming climate, making severe heatwaves, like those observed recently, more common in the future.

Atmospheric blocking↗

DOE OSTI2024

Advancing earth system model calibration: a diffusion-based method

Understanding of climate impact on ecosystems globally requires site-specific model calibration. Here we introduce a novel diffusion-based uncertainty quantification (DBUQ) method for efficient model calibration. DBUQ is a score-based diffusion model that leverages Monte Carlo simulation to estimate the score function and evaluates a simple neural network to quickly generate samples for approximating parameter posterior distributions. DBUQ is stable, efficient, and can effectively calibrate the model given diverse observations, thereby enabling rapid and site-specific model calibration on a global scale. This capability significantly advances Earth system modeling and our understanding of climate impacts on Earth systems. We demonstrate DBUQ's capability in E3SM land model calibration at the Missouri Ozark AmeriFlux forest site. Both synthetic and real-data applications indicate that DBUQ produces accurate parameter posterior distributions similar to those generated by Markov Chain Monte Carlo sampling but with 30X less computing time. This efficiency marks a significant stride in model calibration, paving the way for more effective and timely climate impact analyses.

Liu, Yanfang↗

DOE OSTI2023

Earth System Model Aerosol–Cloud Diagnostics (ESMAC Diags) package, version 2: assessing aerosols, clouds, and aerosol–cloud interactions via field campaign and long-term observations

Abstract. Poor representations of aerosols, clouds, and aerosol–cloud interactions (ACIs) in Earth system models (ESMs) have long been the largest uncertainties in predicting global climate change. Huge efforts have been made to improve the representation of these processes in ESMs, and the key to these efforts is the evaluation of ESM simulations with observations. Most well-established ESM diagnostics packages focus on the climatological features; however, they lack process-level understanding and representations of aerosols, clouds, and ACIs. In this study, we developed the Earth System Model Aerosol–Cloud Diagnostics (ESMAC Diags) package to facilitate the routine evaluation of aerosols, clouds, and ACIs simulated the Energy Exascale Earth System Model (E3SM) from the US Department of Energy (DOE). This paper documents its version 2 functionality (ESMAC Diags v2), which has substantial updates compared with version 1 (Tang et al., 2022a). The simulated aerosol and cloud properties have been extensively compared with in situ and remote-sensing measurements from aircraft, ship, surface, and satellite platforms in ESMAC Diags v2. It currently includes six field campaigns and two permanent sites covering four geographical regions: the eastern North Atlantic, the central US, the northeastern Pacific, and the Southern Ocean. These regions produce frequent liquid- or mixed-phase clouds, with extensive measurements available from the DOE Atmospheric Radiation Measurement user facility and other agencies. ESMAC Diags v2 generates various types of single-variable and multivariable diagnostics, including percentiles, histograms, joint histograms, and heatmaps, to evaluate the model representation of aerosols, clouds, and ACIs. Select examples highlighting the capabilities of ESMAC Diags are shown using E3SM version 2 (E3SMv2). In general, E3SMv2 can reasonably reproduce many observed aerosol and cloud properties, with biases in some variables such as aerosol particle and cloud droplet sizes and number concentrations. The coupling of aerosol and cloud number concentrations may be too strong in E3SMv2, possibly indicating a bias in processes that control aerosol activation. Furthermore, the liquid water path response to a perturbed cloud droplet number concentration behaves differently in E3SMv2 and observations, which warrants further study to improve the cloud microphysics parameterizations in E3SMv2.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Earth System Model Aerosol–Cloud Diagnostics (ESMAC Diags) package, version 1: assessing E3SM aerosol predictions using aircraft, ship, and surface measurements

Abstract. An Earth system model (ESM) aerosol–cloud diagnostics package is developed to facilitate the routine evaluation of aerosols, clouds, and aerosol–cloud interactions simulated by the Energy Exascale Earth System Model (E3SM) from the US Department of Energy (DOE). The first version focuses on comparing simulated aerosol properties with aircraft, ship, and surface measurements, which are mostly measured in situ. The diagnostics currently cover six field campaigns in four geographical regions: eastern North Atlantic (ENA), central US (CUS), northeastern Pacific (NEP), and Southern Ocean (SO). These regions produce frequent liquid- or mixed-phase clouds, with extensive measurements available from the Atmospheric Radiation Measurement (ARM) program and other agencies. Various types of diagnostics and metrics are performed for aerosol number, size distribution, chemical composition, cloud condensation nuclei (CCN) concentration, and various meteorological quantities to assess how well E3SM represents observed aerosol properties across spatial scales. Overall, E3SM qualitatively reproduces the observed aerosol number concentration, size distribution, and chemical composition reasonably well, but it overestimates Aitken-mode aerosols and underestimates accumulation-mode aerosols over the CUS and ENA regions, suggesting that processes related to particle growth or coagulation might be too weak in the model. The current version of E3SM struggles to reproduce the new particle formation events frequently observed over both the CUS and ENA regions, indicating missing processes in current parameterizations. The diagnostics package is coded and organized in a way that can be extended to other field campaign datasets and adapted to higher-resolution model simulations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Energy Exascale Earth System Model v2.0

First release of version 2 of the Energy Exascale Earth System Model. The atmosphere component remains EAM. Major changes since version 1 include: all column-physics parameterizations are computed on a separate grid that has approximately half the number of points of the dynamics grid, a new nonhydrostatic dynamical core (running in hydrostatic mode) with semi-Lagrangian tracer transport, CLUBB updated from v1 to v2, a new convective trigger (dCAPE/ULL) based on the dynamic Convective Available Potential Energy (CAPE) (dCAPE) and the Unrestricted Launch Level (ULL) concepts is used in ZM. minimum cloud droplet number changed, gravity wave drag energy conservation fixed and new tunings used, dust emission size distribution changed to emit more coarse dust particles The land component is still ELM. Major changes since version 1 include: using SNICAR-AD for radiation in snow to match the sea-ice model and fixing bugs in snow compaction and water state calculation. The ocean component remains MPAS-ocean. Major change since version 1 include: Redi isopycnal mixing has been updated, tested, and tuned in combination with the Gent-McWilliams parameterization, a sign error was fixed in the 3rd-order flux routines, the mesh used in low-resolution coupled cases was modified and the time steps adjusted, new regionally refined meshes were created, one focused on North America and another on the Southern Ocean, including ice shelf cavities. The sea-ice component remains MPAS-seaice. Major changes since version 1 include: A new heat- and freshwater-conserving coupling of frazil ice, turning off of SSH filtering, addition of SNICAR-AD and snow grain aging. The land-ice component remains MPAS-Albany-landIce (MALI) and is a static ice sheet. There is more out-of-the-box support for cryosphere configurations including the new regionally refined configuration around Antarctica with ice shelf cavities. The river model is MOSART. The half-degree river mesh was redone so it no longer treats Black and Caspian seas as ocean. The coupler remains cpl7/MCT. Major changes since version 1 include handling of ice shelf melt fluxes (heat / freshwater exchange with the ocean), and data icebergs. All components allow regional refinement of their meshes and two separate example of refinement, one in and around North America and one around Antarctica and the Southern Ocean, are provided.

E3SM Project, DOE↗

DOE OSTI2023

Energy Exascale Earth System Model v2.1.0

First release of version 2.1 of the Energy Exascale Earth System Model. [ATM] The atmosphere component remains EAM. There are no major changes in the default configuration since 2.0. New features include: A semi-lagrangian tracer transport for theta-l dycore, a new algorithm for finding the tropopause, new RRM mesh configurations. Add and update SSP370 and SSP585 cases. Restore the FIDEAL case. [LAND] The land component is ELM. There are no major changes in the default configuration since 2.0. Several option features have been added including: implementation of topography-based subgrid structure (topounits) and accompanying parameterizations and atmospheric forcing downscaling methods; a new plant hydraulics scheme; two-way land-river hydrological coupling through the infiltration of floodplain water; an implementation of perennial crops; updates to the SNICAR-AD snow radiative transfer model; and implementation of soil erosion and sediment yield in ELM-Erosion. Each of these new changes is modular in design and can be turned on or off as the user specifies; they are currently being tested in different “BGC” configurations. [OCEAN] The ocean component remains MPAS-Ocean. Major change since version 2.0 include the addition of the Fox-Kemper et al. 2011 parameterization for submesocale eddies, a correction for barotropic thickness consistency that reduces divergence noise, and the addition of an ocean carbon conservation analysis member. [SEAICE] The sea-ice component remains MPAS-Seaice. Major changes since version 2.0 include: A correction to how shortwave parameters are interpolated in the snicar-ad 5-band radiation scheme, the addition of a sea ice carbon conservation analysis member, updates to the default sea ice biogeochemistry namelist parameters to be consistent with version 2.0 improvements to nitrogen cycling and a correction in the ice-ocean dissolved organic nitrogen coupling. [LAND ICE] The land-ice component remains MPAS-Albany-landIce (MALI). Major changes since 2.0 include an update to the MALI version and the Greenland mesh.[RIVER] The river model is MOSART. There are no major changes in the default configuration since 2.0. A major new optional feature is two-way river-ocean hydrological coupling between MOSART and MPAS-O. This change can be turned on or off as the user specifies, and is being tested in different configurations. [COUPLER] The coupler remains cpl7/MCT. Major changes since version 2.0 include: Carbon budget calculated when heat/water budgets active. Fix a bug in land-atm fluxes for tri-grid configurations. [OTHER] a small bug in the zenith angle calculation was fixed in the data models.

ECP↗

DOE OSTI2022

Assessing the Sensitivity of the Tropical Cyclone Boundary Layer to the Parameterization of Momentum Flux in the Community Earth System Model

Recent studies have demonstrated that high-resolution (~25 km) Earth System Models (ESMs) have the potential to skillfully predict tropical cyclone (TC) occurrence and intensity. However, biases in ESM TCs still exist, largely due to the need to parameterize processes such as boundary layer (PBL) turbulence. Building on past studies, we hypothesize that the depiction of the TC PBL in ESMs is sensitive to the configuration of the PBL parameterization scheme, and that the targeted perturbation of tunable parameters can reduce biases. The Morris one-at-a-time (MOAT) method is implemented to assess the sensitivity of the TC PBL to tunable parameters in the PBL scheme in an idealized configuration of the Community Atmosphere Model, version 6 (CAM6). The MOAT method objectively identifies several parameters in an experimental version of the Cloud Layers Unified by Binormals (CLUBB) scheme that appreciably influence the structure of the TC PBL. We then perturb the parameters identified by the MOAT method within a suite of CAM6 ensemble simulations and find a reduction in model biases compared to observations and a high-resolution, cloud-resolving model. Importantly, we demonstrate that the high-sensitivity parameters are tied to PBL processes that reduce turbulent mixing and effective eddy diffusivity, and that in CAM6 these parameters alter the TC PBL in a manner consistent with past modeling studies. In this way, we provide an initial identification of process-based input parameters that, when altered, have the potential to improve TC predictions by ESMs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Metrics for Intercomparison of Remapping Algorithms (MIRA) protocol applied to Earth system models

Abstract. Strongly coupled nonlinear phenomena such as those described by Earth system models (ESMs) are composed of multiple component models with independent mesh topologies and scalable numerical solvers. A common operation in ESMs is to remap or interpolate component solution fields defined on their computational mesh to another mesh with a different combinatorial structure and decomposition, e.g., from the atmosphere to the ocean, during the temporal integration of the coupled system. Several remapping schemes are currently in use or available for ESMs. However, a unified approach to compare the properties of these different schemes has not been attempted previously. We present a rigorous methodology for the evaluation and intercomparison of remapping methods through an independently implemented suite of metrics that measure the ability of a method to adhere to constraints such as grid independence, monotonicity, global conservation, and local extrema or feature preservation. A comprehensive set of numerical evaluations is conducted based on a progression of scalar fields from idealized and smooth to more general climate data with strong discontinuities and strict bounds. We examine four remapping algorithms with distinct design approaches, namely ESMF Regrid (Hill et al., 2004), TempestRemap (Ullrich and Taylor, 2015), generalized moving least squares (GMLS) (Trask and Kuberry, 2020) with post-processing filters, and WLS-ENOR (Li et al., 2020). By repeated iterative application of the high-order remapping methods to the test fields, we verify the accuracy of each scheme in terms of their observed convergence order for smooth data and determine the bounded error propagation using challenging, realistic field data on both uniform and regionally refined mesh cases. In addition to retaining high-order accuracy under idealized conditions, the methods also demonstrate robust remapping performance when dealing with non-smooth data. There is a failure to maintain monotonicity in the traditional L2-minimization approaches used in ESMF and TempestRemap, in contrast to stable recovery through nonlinear filters used in both meshless GMLS and hybrid mesh-based WLS-ENOR schemes. Local feature preservation analysis indicates that high-order methods perform better than low-order dissipative schemes for all test cases. The behavior of these remappers remains consistent when applied on regionally refined meshes, indicating mesh-invariant implementations. The MIRA intercomparison protocol proposed in this paper and the detailed comparison of the four algorithms demonstrate that the new schemes, namely GMLS and WLS-ENOR, are competitive compared to standard conservative minimization methods requiring computation of mesh intersections. The work presented in this paper provides a foundation that can be extended to include complex field definitions, realistic mesh topologies, and spectral element discretizations, thereby allowing for a more complete analysis of production-ready remapping packages.

58 GEOSCIENCES↗