Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earth system modeling”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Ensemble-Based Characterization of Historical Storm Types in ERA5 over CONUS

We present Storm-Type Labeled Precipitation, a gridded dataset that classifies ERA5 precipitation over CONUS (20–50° N, 125–66° W) by storm type at 6-hourly resolution. Each grid cell/time step is assigned to one of five classes—mesoscale convective system (MCS), extratropical cyclone (ETC), hurricane (HUR), atmospheric river (AR), or other convective, with an “unidentified” code reserved for cases with no detected type. Two products are provided: (i) a TempestExtremes (TE)–only version, and (ii) a hybrid TRACK–TE version that uses TRACK for ETCs and HURs and TE for MCSs and ARs. Outputs are integer masks aligned to ERA5 precipitation, supporting storm-type attribution of rainfall, event compositing, trend analysis, and model evaluation. Future work will extend detection to high-resolution downscaled Earth system model projections to support evaluations of projected changes in hydrologic hazards that may threaten critical water and energy infrastructures.

Rastogi, Deeksha [ORNL] (ORCID:0000000204624027)↗
NASA NTRSDate not supplied

Challenges in Observing, Modeling, and Forecasting the June 2023 Smoke Event Over the Northeast United States

In the first week of May 2023, boreal Canada began an early start to the biomass burning season when fires erupted across Alberta and western Saskatchewan. Smoke entered the troposphere and spread across the United States as biomass burning emissions in Canada quadrupled the prior maximum since the onset of the MODIS satellite record. Adding to the already existing smoke, wildfires ignited in Quebec on June 2, 2023, producing a heavy plume of aerosol over the northeastern part of the continent. Meteorological factors, including an area of low pressure situated, nearly stationary, over Maine, transported the smoke into a densely populated corridor containing the cities of Washington, DC, Baltimore, Philadelphia, and New York City. Air quality alerts became widespread as near-surface levels of fine particulate matter (PM2.5) exceeded harmful levels due the low altitude of the smoke plume. Forecasting the smoke transport and impacts of such events requires a complex combination of observed initial conditions for meteorology and aerosols, knowledge of the emissions from wildfires, and a state-of-the-art Earth System model that couples these components together. Using a case-study perspective with the Goddard Earth Observing System (GEOS), challenges associated with near real time forecasting of the smoke plume are presented. Analyzed and forecasted aerosol optical depth will be evaluated using available satellite and AERONET observations, while surface aerosol will be assessed relative to observations of PM2.5. Lacking real time estimates of biomass burning emissions for inputs to our forecast model, particular attention is given to the use of day-old biomass burning emissions throughout the forecast period covering the first ten days in June 2023, and the resulting reduction in forecast skill will be quantified.

Allison Collow↗
NASA NTRSDate not supplied

Challenges in Observing, Modeling, and Forecasting the June 2023 Smoke Event over the Northeast United States

In the first week of May 2023, boreal Canada began an early start to the biomass burning season when fires erupted across Alberta and western Saskatchewan. Smoke entered the troposphere and spread across the United States as biomass burning emissions in Canada quadrupled the prior maximum since the onset of the MODIS satellite record. Adding to the already existing smoke, wildfires ignited in Quebec on June 2, 2023, producing a heavy plume of aerosol over the northeastern part of the continent. Meteorological factors, including an area of low pressure situated, nearly stationary, over Maine, transported the smoke into a densely populated corridor containing the cities of Washington, DC, Baltimore, Philadelphia, and New York City. Air quality alerts became widespread as near-surface levels of fine particulate matter (PM2.5) exceeded harmful levels due the low altitude of the smoke plume. Forecasting the smoke transport and impacts of such events requires a complex combination of observed initial conditions for meteorology and aerosols, knowledge of the emissions from wildfires, and a state-of-the-art Earth System model that couples these components together. Using a case-study perspective with the Goddard Earth Observing System (GEOS ), challenges associated with near real time forecasting of the smoke plume are presented. Analyzed and forecasted aerosol optical depth will be evaluated using available satellite and AERONET observations, while surface aerosol will be assessed relative to observations of PM2.5. Lacking real time estimates of biomass burning emissions for inputs to our forecast model, particular attention is given to the use of day-old biomass burning emissions throughout the forecast period covering the first ten days in June 2023, and the resulting reduction in forecast skill will be quantified.

Allison Collow↗
DOE OSTI2025

Evaluation of CMIP6 Streamflow in the Arctic

Earth system models (ESMs) serve as the primary basis for projecting future streamflow changes, but they are biased in terms of the their ability to reproduce historical observations of streamflow. Rigorous downscaling and bias-correction procedures, which are time-consuming and introduce uncertainties, are commonly used to address these biases in streamflow or runoff. However, the applicable limits of streamflow projections remain largely untested through direct comparisons with observations. In this work, we compare historical streamflow time series observations with a suite of Coupled Model Intercomparison Project phase 6 (CMIP6) ESMs to assess the space and time scales at which each ESM represents historic streamflow observations. We focus our analysis on streamflow across the Arctic given the higher degree of warming experienced over this region, underscoring its importance for understanding the future changes. A series of metrics—volume, seasonality, extreme capturing streamflow variability, drought, and flooding is tested to assess ESM simulation skill for total events and overall distributions. Our results indicate that although improvements are necessary in ESM runoff and streamflow projections, the use of temporal averaging and targeting the best-fit model for comparison to observations when a suite of ESMs is considered can provide robust streamflow projections. Moreover, this work provides a basis for contextualizing the limits of ESM applicability for studies aiming to assess streamflow in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Improving thermodynamic nudging in the E3SM Atmosphere Model version 2 (EAMv2): strategy and hindcast skills on weather systems

Nudging techniques are commonly employed to constrain atmospheric simulations toward observed states, facilitating model evaluation and sensitivity studies. However, if applied improperly – particularly to thermodynamic variables such as temperature and humidity – nudging can distort physical processes and introduce spurious biases, undermining the credibility of the simulations. This study presents an improved nudging implementation that applies vertically modulated tendencies to reduce adverse impacts on model physics. The framework is tested in version 2 of the Energy Exascale Earth System Model (EAMv2) using a suite of hindcast simulations nudged toward ERA5 reanalysis. We systematically evaluate the individual and combined effects of nudging wind, temperature, and humidity fields on the model's ability to represent large-scale atmospheric states and high-impact weather systems. Results show that the revised strategy – particularly when nudging temperature and humidity at selected levels – enhances hindcast skill by improving agreement with ERA5 without degrading the hydrological cycle or precipitation processes. Additional improvements in surface temperature, outgoing longwave radiation, and precipitation biases are achieved through targeted nudging of land-surface variables. The proposed approach strengthens the representation of large-scale conditions relevant to tropical cyclones, atmospheric rivers, and extratropical cyclones in the low-resolution EAMv2. These findings demonstrate that carefully designed thermodynamic nudging, especially of temperature and humidity, improves the realism of constrained simulations and broadens the utility of nudged EAMv2 for atmospheric modeling, machine learning, and high-impact weather research.

Atmospheric river↗
NASA NTRSDate not supplied

Earth System Digital Twins (ESDT) Technology for NASA Earth Science

For NASA's Advanced Information Systems Technology (AIST) Program, an Earth System Digital Twin (ESDT) is defined as an interactive and integrated multidomain, multiscale, digital replica of the state and temporal evolution of Earth systems. It dynamically integrates: relevant Earth system models and simulations; other relevant models (e.g., related to the world's infrastructure); continuous and timely (including near real time and direct readout) observations (e.g., space, air, ground, over/underwater, Internet of Things (IoT), socioeconomic); long-time records; as well as analytics and artificial intelligence tools. Effective ESDTs enable users to run hypothetical scenarios to improve the understanding, prediction of and mitigation/response to Earth system processes, natural phenomena and human activities as well as their many interactions. An ESDT is a type of integrated information system that, for example, enables continuous assessment of impact from naturally occurring and/or human activities on physical and natural environments. AIST ESDT strategic goals are to: 1. Develop information system frameworks to provide continuous and accurate representations of systems as they change over time; 2. Mirror various Earth Science systems and utilize the combination of Data Analytics, Artificial Intelligence, Digital Thread, and state-of-the-art models to help predict the Earth’s response to various phenomena; 3. Provide the tools to conduct "what if" investigations that can result in actionable predictions. The AIST ESDT thrust is developing capabilities toward the development of future digital twins of the Earth or of subcomponents of the Earth. This will enable the development of an overarching framework that will integrate New Observing Strategies (NOS) to enable new observation measurements, i.e., multi-source, coordinated, dynamic and responsive to needs and requests defined by Analytic Collaborative Frameworks (ACF) that enable agile science investigations fusing and analyzing very large amounts of diverse data. NOS and ACF capabilities along with open access to various science, infrastructure and human data, interconnected modeling, data assimilation, simulations, surrogate modeling, high-performance computing and advanced visualization, will define a powerful framework that could be utilized for local, regional or global and/or thematic digital twins. This presentation will describe a general overview of the AIST ESDT vision including prior work done in the areas of NOS and ACF as well as current and upcoming ESDT projects.

Jacqueline Le Moigne↗
DOE OSTI2021

Evaluation and modification of ELM seasonal deciduous phenology against observations in a southern boreal peatland forest

Phenological transitions determine the timing of changes in land surface properties and the seasonality of exchanges of biosphere-atmosphere energy, water, and carbon. Accurate mechanistic modeling of phenological processes is therefore critical to understand and correctly predict terrestrial ecosystem feedbacks with changing atmospheric and climate conditions. However, the phenological components in the land model of the US Department of Energy's (DOE) Energy Exascale Earth System Model (ELM of E3SM) were previously unable to accurately capture the observed phenological responses to environmental conditions in a well-studied boreal peatland forest. In this research, we introduced new seasonal-deciduous phenology schemes into version 1.0 of ELM and evaluated their performance against the PhenoCam observations at the Spruce and Peatland Responses Under Changing Environments (SPRUCE) experiment in northern Minnesota from 2015 to 2018. We found that phenology simulated by the revised ELM (i.e., earlier spring onsets and stronger warming responses of spring onsets and autumn senescence) was closer to observations than simulations from the original algorithms for both the deciduous conifer (Larix laricina) and mixed shrub layers. Moreover, the revised ELM generally produced higher carbon and water fluxes (e.g., photosynthesis and evapotranspiration) during the growing season and stronger flux responses to warming than the default ELM. A parameter sensitivity analysis further indicated the significant contribution of phenology parameters to uncertainty in key carbon and water cycle variables, underscoring the importance of precise phenology parameterization. Furthermore, this phenological modeling effort demonstrates the potential to enhance the E3SM representation of land-climate interactions at broader spatiotemporal scales, especially under anticipated elevated CO 2 and warming conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Equation‐Free Surrogate Modeling of Geophysical Flows at the Intersection of Machine Learning and Data Assimilation

Abstract There is a growing interest in developing data‐driven reduced‐order models for atmospheric and oceanic flows that are trained on data obtained either from high‐resolution simulations or satellite observations. The data‐driven models are non‐intrusive in nature and offer significant computational savings compared to large‐scale numerical models. These low‐dimensional models can be utilized to reduce the computational burden of generating forecasts and estimating model uncertainty without losing the key information needed for data assimilation (DA) to produce accurate state estimates. This paper aims at exploring an equation‐free surrogate modeling approach at the intersection of machine learning and DA in Earth system modeling. With this objective, we introduce an end‐to‐end non‐intrusive reduced‐order modeling (NIROM) framework equipped with contributions in modal decomposition, time series prediction, optimal sensor placement, and sequential DA. Specifically, we use proper orthogonal decomposition (POD) to identify the dominant structures of the flow, and a long short‐term memory network to model the dynamics of the POD modes. The NIROM is integrated within the deterministic ensemble Kalman filter (DEnKF) to incorporate sparse and noisy observations at optimal sensor locations obtained through QR pivoting. The feasibility and the benefit of the proposed framework are demonstrated for the NOAA Optimum Interpolation Sea Surface Temperature (SST) V2 data set. Our results indicate that the NIROM is stable for long‐term forecasting and can model dynamics of SST with a reasonable level of accuracy. Furthermore, the prediction accuracy of the NIROM gets improved by almost one order of magnitude by the DEnKF algorithm.

Pawar, Suraj↗
DOE OSTI2021

Multi-Scale Land-Atmosphere Interactions: Modeling Convective Processes from Plants to Planet

Research accomplishments include: 1)Two case studies of the effects of heterogeneous soil moisture and surface energy budgets on organization and propagation of convective precipitation during MC3E. 2) Development and evaluation of a new approach for simulating the effect of heterogeneous soil moisture at ARM-SGP using an innovative modeling approach. 3) Investigation of the effects of spatial coupling scale using multidecade global simulations in CESM with the multiscale modeling framework. 4) Exploration of changes to future precipitation intensity resulting from two different climate change scenarios using the multiscale model. 5) Provision of the new cloud-scale coupled multiscale Earth System Model to the larger community through the CESM process.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Science & Technology Review: Climate Change Comes into Focus

At Lawrence Livermore National Laboratory, we focus on science and technology research to ensure our nation’s security. We also apply that expertise to solve other important national problems in energy, bioscience, and the environment. Science & Technology Review is published eight times a year to communicate, to a broad audience, the Laboratory’s scientific and technological accomplishments in fulfilling its primary missions. The publication’s goal is to help readers understand these accomplishments and appreciate their value to the individual citizen, the nation, and the world. The cover image shows a simulation run on the Department of Energy’s (DOE’s) state-of-the-science Energy Exascale Earth Systems Model (E3SM), which leverages DOE’s advanced supercomputers for climate and earth science research. As the article beginning on p. 4 describes, the newly released version 2 will be used to simulate aspects of Earth’s variability at weather-scale resolution and investigate decadal changes in climate that will critically impact the United States in coming years.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Development of inter-grid-cell lateral unsaturated and saturated flow model in the E3SM Land Model (v2.0)

Abstract. The lateral transport of water in the subsurface is important in modulating terrestrial water energy distribution. Although a few land surface models have recently included lateral saturated flow within and across grid cells, it is not a default configuration in the Climate Model Intercomparison Project version 6 experiments. In this work, we developed the lateral subsurface flow model within both unsaturated and saturated zones in the Energy Exascale Earth System Model (E3SM) Land Model version 2 (ELMv2.0). The new model, called ELMlat, was benchmarked against PFLOTRAN, a 3D subsurface flow and transport model, for three idealized hillslopes that included a convergent hillslope, divergent hillslope, and tilted V-shaped hillslope with variably saturated initial conditions. ELMlat showed comparable performance against PFLOTRAN in terms of capturing the dynamics of soil moisture and groundwater table for the three benchmark hillslope problems. Specifically, the mean absolute errors (MAEs) of the soil moisture in the top 10 layers between ELMlat and PFLOTRAN were within 1 %±3 %, and the MAEs of water table depth were within ±0.2 m. Next, ELMlat was applied to the Little Washita experimental watershed to assess its prediction of groundwater table, soil moisture, and soil temperature. The spatial pattern of simulated groundwater table depth agreed well with the global groundwater table benchmark dataset generated from a global model calibrated with long-term observations. The effects of lateral groundwater flow on the energy flux partitioning were more prominent in lowland areas with shallower groundwater tables, where the difference in simulated annual surface soil temperature could reach 0.3–0.4 ∘C between ELMv2.0 and ELMlat. Incorporating lateral subsurface flow in ELM improves the representation of the subsurface hydrology, which will provide a good basis for future large-scale applications.

58 GEOSCIENCES↗
DOE OSTI2023

Effective climate sensitivity distributions from a 1D model of global ocean and land temperature trends, 1970–2021

Abstract Current theoretically based Earth system models (ESMs) produce Effective Climate Sensitivities (EffCS) that range over a factor of three, with 80% of those models producing stronger global warming trends for 1970–2021 than do observations. To make a more observationally based estimate of EffCS, a 1D time-dependent forcing-feedback model of temperature departures from energy equilibrium is used to match measured ranges of global-average surface and sub-surface land and ocean temperature trends during 1970–2021. In response to two different radiative forcing scenarios, a full range of three model free parameters are evaluated to produce fits to a range of observed surface temperature trends (± 2σ) from four different land datasets and three ocean datasets, as well as deep-ocean temperature trends and borehole-based trend retrievals over land. Land-derived EffCS are larger than over the ocean, and EffCS is lower using the newer Shared Socioeconomic Pathways (SSP245, 1.86 °C global EffCS, ± 34% range 1.48–2.15 °C) than the older Representative Concentration Pathway forcing (RCP6, 2.49 °C global average EffCS, ± 34% range 2.04–2.87 °C). The strongest dependence of the EffCS results is on the assumed radiative forcing dataset, underscoring the role of radiative forcing uncertainty in determining the sensitivity of the climate system to increasing greenhouse gas concentrations from observations alone. The results are consistent with previous observation-based studies that concluded EffCS during the observational period is on the low end of the range produced by current ESMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Comprehensive Analysis of Uncertainties in Warm-Rain Parameterizations in Climate Models Based on In Situ Measurements

Abstract Because of the coarse grid size of Earth system models (ESMs), representing warm-rain processes in ESMs is a challenging task involving multiple sources of uncertainty. Previous studies evaluated warm-rain parameterizations mainly according to their performance in emulating collision–coalescence rates for local droplet populations over a short period of a few seconds. The representativeness of these local process rates comes into question when applied in ESMs for grid sizes on the order of 100 km and time steps on the order of 20–30 min. We evaluate several widely used warm-rain parameterizations in ESM application scenarios. In the comparison of local and instantaneous autoconversion rates, the two parameterization schemes based on numerical fitting to stochastic collection equation (SCE) results perform best. However, because of Jessen’s inequality, their performance deteriorates when grid-mean, instead of locally resolved, cloud properties are used in their simulations. In contrast, the effect of Jessen’s inequality partly cancels the overestimation problem of two semianalytical schemes, leading to an improvement in the ESM-like comparison. In the assessment of uncertainty due to the large time step of ESMs, it is found that the rainwater tendency simulated by the SCE is roughly linear for time steps smaller than 10 min, but the nonlinearity effect becomes significant for larger time steps, leading to errors up to a factor of 4 for a time step of 20 min. After considering all uncertainties, the grid-mean and time-averaged rainwater tendency based on the parameterization schemes is mostly within a factor of 4 of the local benchmark results simulated by SCE.

Meteorology & Atmospheric Sciences↗
DOE OSTI2026

Gridded Vegetation and Land Unit Datasets over North America (1km x 1km) for E3SM Land Model, Version 2

Gridded 1 km land-surface-property dataset for Energy Exascale Earth System Model Land Model (ELM) simulations over North America. Thirty-metre NALCMS land cover is aggregated to the Daymet grid; vegetation classes are cross-walked to ELM plant functional types using mean-temperature-of-the-coldest-month (MTCO) rules. NetCDF includes PFT and land-unit fractions, counts, and MTCO.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Model Data Archive Associated with Manuscript "Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon"

This data package supports the publication “Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon” by Li et al. (2026). The package contains processed model inputs, configuration files, restart files, simulation outputs, scripts, and visualization products used to evaluate post-fire dissolved organic carbon (DOC) dynamics in the Naches River Watershed, Washington, USA, following the 2021 Schneider Springs Fire. The modeling workflow couples ELM-BGC, the biogeochemistry-enabled Energy Exascale Earth System Model Land Model; ATS, the Advanced Terrestrial Simulator for integrated surface-subsurface hydrology; and PFLOTRAN, a reactive transport model for multicomponent aqueous geochemistry. Together, these models simulate how wildfire-induced changes in vegetation, litter, coarse woody debris, and soil organic matter influence DOC production, transport, and reaction from burned hillslopes to stream networks. The archive includes preprocessed meteorological, geospatial, hydrologic, and biogeochemical forcing data; ELM-BGC-derived DOC source terms; ATS mesh files; PFLOTRAN reactive-transport inputs; model configuration files; spin-up and transient restart files; watershed-scale diagnostic outputs; stream concentration time series; and figures or visualization files used to inspect and reproduce key results. File types include Hierarchical Data Format 5 (HDF5) files for gridded forcing and model-coupling data, model input and configuration files for ELM-BGC, ATS, and PFLOTRAN, restart and simulation-output files generated by the modeling workflow, tabular or time-series diagnostic outputs, scripts for post-processing and figure generation, and image or visualization products associated with the manuscript. Use of the package depends on the intended task. Re-running the simulations requires the relevant modeling software, including ELM-BGC, ATS, and PFLOTRAN as ATS's geochemical engine. Inspecting outputs and reproducing figures requires Python with scientific plotting libraries such as Matplotlib, and three-dimensional model outputs may be viewed with ParaView. Geographic information system files or maps may be inspected with ArcGIS Pro or comparable GIS software. The data package is intended to enable traceability, reuse, and partial reproduction of the coupled land-to-watershed hydro-biogeochemical modeling workflow used to test how wildfire disturbance affects terrestrial carbon pools and downstream DOC dynamics.

ATS↗
DOE OSTI2025

How initial conditions-, structural-, and parameter-based model uncertainty interact and influence predictions in permafrost ecosystems: Modeling Archive

This dataset contains model output and input data, as well as source code examples for the Terrestrial Ecosystem Model with the Dynamic Vegetation Model and Dynamic Organic Soil (DVM-DOS-TEM) for the field sites Imnavait creek and the Bonanza creek Long Term Ecological Research Network (LTER). The data covers simulations from the last glacial maximum (LGM) until 2100 for a selection of paleo scenarios, setting the mean temperature of the LGM up to 10°C lower than pre-industrial conditions. The model structure was modulated to represent various model versions, and this dataset contains the relevant changes in the source code. The raw output data, the processed statistical data, the setup and processing scripts as well as parameter value distribution files from a parameter sensitivity analysis are included as well. Model outputs include active layer depth, organic soil carbon, soil layer depths, gross primary productivity (GPP) with and without nitrogen limitation, net primary productivity (NPP), soil liquid water content, heterotrophic, maintenance, and growth respiration, soil temperature, and vegetation carbon (*.nc files). The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research.Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A hierarchy of global ocean models coupled to CESM1

We develop a hierarchy of simplified ocean models for coupled ocean, atmosphere, and sea ice climate simulations using the Community Earth System Model version 1 (CESM1). The hierarchy has four members: a slab ocean model, a mixed-layer model (MLM) with entrainment and detrainment, an Ekman MLM, and an ocean general circulation model (OGCM). Flux corrections of heat and salt are applied to the simplified models ensuring that all hierarchy members have the same climatology. We diagnose the needed flux corrections from auxiliary simulations in which we restore the temperature and salinity to the daily climatology obtained from a target CESM1 simulation. The resulting three-dimensional corrections contain the interannual variability fluxes that maintain the correct vertical gradients of temperature and salinity in the tropics. We find that the inclusion of mixed-layer entrainment and Ekman flow produces sea surface temperature and surface air temperature fields whose means and variances are progressively more similar to those produced by the target CESM1 simulation. We illustrate the application of the hierarchy to the problem of understanding the response of the climate system to the loss of Arctic sea ice. We find that the shifts in the positions of the mid-latitude westerly jet and of the Inter-tropical Convergence Zone (ITCZ) in response to sea-ice loss depend critically on upper ocean processes. Specifically, heat uptake associated with the mixed-layer entrainment influences the shift in the westerly jet and ITCZ. Moreover, the shift of ITCZ is sensitive to the form of Ekman flow parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

MPAS-Seaice (v1.0.0): sea-ice dynamics on unstructured Voronoi meshes

Abstract. We present MPAS-Seaice, a sea-ice model which uses the Model for Prediction Across Scales (MPAS) framework and spherical centroidal Voronoi tessellation (SCVT) unstructured meshes. As well as SCVT meshes, MPAS-Seaice can run on the traditional quadrilateral grids used by sea-ice models such as CICE. The MPAS-Seaice velocity solver uses the elastic–viscous–plastic (EVP) rheology and the variational discretization of the internal stress divergence operator used by CICE, but adapted for the polygonal cells of MPAS meshes, or alternatively an integral (“finite-volume”) formulation of the stress divergence operator. An incremental remapping advection scheme is used for mass and tracer transport. We validate these formulations with idealized test cases, both planar and on the sphere. The variational scheme displays lower errors than the finite-volume formulation for the strain rate operator but higher errors for the stress divergence operator. The variational stress divergence operator displays increased errors around the pentagonal cells of a quasi-uniform mesh, which is ameliorated with an alternate formulation for the operator. MPAS-Seaice shares the sophisticated column physics and biogeochemistry of CICE and when used with quadrilateral meshes can reproduce the results of CICE. We have used global simulations with realistic forcing to validate MPAS-Seaice against similar simulations with CICE and against observations. We find very similar results compared to CICE, with differences explained by minor differences in implementation such as with interpolation between the primary and dual meshes at coastlines. We have assessed the computational performance of the model, which, because it is unstructured, runs with 70 % of the throughput of CICE for a comparison quadrilateral simulation. The SCVT meshes used by MPAS-Seaice allow removal of equatorial model cells and flexibility in domain decomposition, improving model performance. MPAS-Seaice is the current sea-ice component of the Energy Exascale Earth System Model (E3SM).

58 GEOSCIENCES↗