Search NASA⌕ Search

SEARCH · Search NASA

Results for “circulation 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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

MITgcm-AD v2: Open source tangent linear and adjoint modeling framework for the oceans and atmosphere enabled by the Automatic Differentiation tool Tapenade

The Massachusetts Institute of Technology General Circulation Model (MITgcm) is widely used by the climate science community to simulate planetary atmosphere and ocean circulations. A defining feature of the MITgcm is that it has been developed to be compatible with an algorithmic differentiation (AD) tool, TAF, enabling the generation of tangent-linear and adjoint models. These provide gradient information which enables dynamics-based sensitivity and attribution studies, state and parameter estimation, and rigorous uncertainty quantification. Importantly, gradient information is essential for computing comprehensive sensitivities and performing efficient large-scale data assimilation, ensuring that observations collected from satellites and in-situ measuring instruments can be effectively used to optimize a large uncertain control space. As a result, the MITgcm forms the dynamical core of a key data assimilation product employed by the physical oceanography research community: Estimating the Circulation and Climate of the Ocean (ECCO) state estimate. Although MITgcm and ECCO are used extensively within the research community, the AD tool TAF is proprietary and hence inaccessible to a large proportion of these users. The new version 2 (MITgcm-AD v2) framework introduced here is based on the source-to-source AD tool Tapenade, which has recently been open-sourced. Another feature of Tapenade is that it stores required variables by default (instead of recomputing them) which simplifies the implementation of efficient, AD-compatible code. The framework has been integrated with the MITgcm model’s main branch and is now freely available.

Adjoints↗
DOE OSTI2020

Evaluation of the Quasi‐Biennial Oscillation in global climate models for the SPARC QBO‐initiative

Abstract Quasi‐biennial oscillations (QBOs) in thirteen atmospheric general circulation models forced with both observed and annually repeating sea surface temperatures (SSTs) are evaluated. In most models the QBO period is close to, but shorter than, the observed period of 28 months. Amplitudes are within ±20 % of the observed QBO amplitude at 10 hPa, but typically about half of that observed at lower altitudes (50 and 70 hPa). For almost all models, the oscillation's amplitude profile shows an overall upward shift compared to reanalysis and its meridional extent is too narrow. Asymmetry in the duration of eastward and westward phases is reasonably well captured, though not all models replicate the observed slowing of the descending westward shear. Westward phases are generally too weak, and most models have an eastward time mean wind bias throughout the depth of the QBO. The intercycle period variability is realistic and in some models is enhanced in the experiment with observed SSTs compared to the experiment with repeated annual cycle SSTs. Mean periods are also sensitive to this difference between SSTs, but only when parametrized non‐orographic gravity wave (NOGW) sources are coupled to tropospheric parameters and not prescribed with a fixed value. Overall, however, modelled QBOs are very similar whether or not the prescribed SSTs vary interannually. A portrait of the overall ensemble performance is provided by a normalized grading of QBO metrics. To simulate a QBO, all but one model used parametrized NOGWs, which provided the majority of the total wave forcing at altitudes above 70 hPa in most models. Hence the representation of NOGWs either explicitly or through parametrization is still a major uncertainty underlying QBO simulation in these present‐day experiments.

Bushell, A. C.↗
DOE OSTI2022

Convective Momentum Transport and Its Impact on the Madden‐Julian Oscillation in E3SM‐MMF

Abstract Convective momentum transport (CMT) is the process of vertical redistribution of horizontal momentum by small‐scale turbulent flows from moist convection. Traditional general circulation models (GCMs) and their multiscale modeling framework (MMF) versions poorly represent CMT due to insufficient information of subgrid‐scale flows at each GCM grid. Here the explicit scalar momentum transport (ESMT) scheme for representing CMT is implemented in the Energy Exascale Earth System Model‐Multiscale Modeling Framework (E3SM‐MMF) with embedded 2‐D cloud‐resolving models (CRMs), and verified against E3SM‐MMF simulations with 3‐D CRMs and observations. The results show that representing CMT by ESMT helps reduce climatological mean precipitation model bias over the western Pacific and the ITCZ regions, which is attributed to the weakened mean easterlies over the Pacific. Also, CMT from simulations with 2‐D and 3‐D CRMs impose a similar impact on Kelvin waves by reducing their variability and slowing down their phase speed, but opposite impacts on the Madden‐Julian Oscillation (MJO) variability. The ESMT scheme readily captures the climatological mean spatial patterns of the zonal and meridional components of CMT and their variability across multiple time scales, but shows some differences in estimating its magnitude. CMT mainly affects the MJO by decelerating its winds in the free troposphere, but accelerates its near‐surface winds. This study serves as a prototype for implementing CMT scheme in the MMF simulations, highlighting its crucial role in reducing model bias in mean state and spatiotemporal variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A Moist Physics Parameterization Based on Deep Learning

Abstract Current moist physics parameterization schemes in general circulation models (GCMs) are the main source of biases in simulated precipitation and atmospheric circulation. Recent advances in machine learning make it possible to explore data‐driven approaches to developing parameterization for moist physics processes such as convection and clouds. This study aims to develop a new moist physics parameterization scheme based on deep learning. We use a residual convolutional neural network (ResNet) for this purpose. It is trained with 1‐year simulation from a superparameterized GCM, SPCAM. An independent year of SPCAM simulation is used for evaluation. In the design of the neural network, referred to as ResCu, the moist static energy conservation during moist processes is considered. In addition, the past history of the atmospheric states, convection, and clouds is also considered. The predicted variables from the neural network are GCM grid‐scale heating and drying rates by convection and clouds, and cloud liquid and ice water contents. Precipitation is derived from predicted moisture tendency. In the independent data test, ResCu can accurately reproduce the SPCAM simulation in both time mean and temporal variance. Comparison with other neural networks demonstrates the superior performance of ResNet architecture. ResCu is further tested in a single‐column model for both continental midlatitude warm season convection and tropical monsoonal convection. In both cases, it simulates the timing and intensity of convective events well. In the prognostic test of tropical convection case, the simulated temperature and moisture biases with ResCu are smaller than those using conventional convection and cloud parameterizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Hierarchy of Global Ocean Models Coupled to CESM1

Abstract 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

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 OSTI2021

Characterization of oceanic post-cold frontal clouds and their model representation. Final report

Low-level clouds are ubiquitous over the oceans, and their impact on the Earth’ radiative budget is important. However, general circulation models (GCMs) that are used to represent our climate and its evolution still experience problems to represent their correct amount and radiative impact. In the extratropical latitudes (~30-60 N/S), cloud properties are strongly modulated by the occurrence of anti-cyclones (high pressure systems) and cyclones (low pressure systems). While deep clouds dominate in cyclones near the center and along the cold and warm fronts, low-level clouds populate the colder regions of the cyclones, and more specifically the region in the wake of the cold fronts, the post-cold frontal region (PCF). The GCMs mentioned above were found to underestimate the amount of these PCF clouds, and this caused errors in their representation of the amount of solar radiation at the surface of the southern hemisphere oceans, causing an excess in absorption and errors in long-term predictions. To better understand which processes are ill-represented in the GCMs to cause such issues, we used ground-based observations from the Eastern North Atlantic (ENA) ARM site and a high resolution regional model (called WRF) to examine the properties of PCF clouds. With the WRF model, we simulated the passage of a cold front at the ENA site using theoretically distinct representations of convection and the physics of the boundary layer. While these physical representations had little impact on the timing, structure and circulation of the cold front, we found a high sensitivity to how the large scale information was fed into the model. Errors rapidly developed if the incoming flow information was injected too far from the site, but 1000 km was found to be an optimal distance between the outer boundary and the ENA site for a realistic cold front passage. With this model, we then conducted a series of experiment where convection and boundary layer schemes were changed, and revealed that 1) the clouds were more sensitive to the convection than boundary layer schemes and 2) the interaction between the two has a large impact. Furthermore, changing the representation of the physics in the model had greater impacts on PCF clouds than changing the initial conditions. As large-scale models increase their spatial resolution, convection will eventually be resolved but it remains that their representation of the physics of the boundary layer will matter for the correct representation of shallow (low-level) clouds. Observations were used to find which large-scale information had a clear relationship with PCF cloud and precipitation properties, in view of helping with GCM parameterizations and evaluation. Both lower troposphere stability and surface forcing were found to be important for PCF clouds that experience more vigorous dynamics than their tropical or subtropical counterparts. For cloud base and top heights, a simple measure of the change in potential temperature between the surface and the 800 hPa was found to be a good predictor. This relationship was found to hold as well when using data from campaigns recently conducted in the southern ocean south of Tasmania. A recent GCM was found to succeed in representing similar relations. For cloud depth, water content, and different characteristics of precipitation both the temperature contrast between the surface and near-air and the near-surface wind speed had a strong relationship. To summarize, through observational constraints and regional model experiments, the project has provided guidance to evaluate GCMs and tools to improve their representation of low-level clouds in the extratropics. It has also highlighted differences in behavior between these clouds and those typically found in the subtropics, thus opening new avenues of research for shallow clouds in general.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Response of the Quasi–Biennial Oscillation to a warming climate in global climate models

We compare the response of the Quasi–Biennial Oscillation (QBO) to a warming climate in eleven atmosphere general circulation models that performed time–slice simulations for present–day, doubled, and quadrupled CO 2 climates. No consistency was found among the models for the QBO period response, with the period decreasing by 8 months in some models and lengthening by up to 13 months in others in the doubled CO 2 simulations. In the quadrupled CO 2 simulations, a reduction in QBO period of 14 months was found in some models, whereas in several others the tropical oscillation no longer resembled the present–day QBO, although it could still be identified in the deseasonalized zonal mean zonal wind timeseries. In contrast, all the models projected a decrease in the QBO amplitude in a warmer climate with the largest relative decrease near 60 hPa. In simulations with doubled and quadrupled CO 2 , the multi–model mean QBO amplitudes decreased by 36 and 51%, respectively. Across the models the differences in the QBO period response were most strongly related to how the gravity wave momentum flux entering the stratosphere and tropical vertical residual velocity responded to the increases in CO 2 amounts. Likewise it was found that the robust decrease in QBO amplitudes was correlated across the models to changes in vertical residual velocity, parametrized gravity wave momentum fluxes, and to some degree the resolved upward wave flux. Here, we argue that uncertainty in the representation of the parameterized gravity waves is the most likely cause of the spread among the eleven models in the QBO's response to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Nitrate loading projection is sensitive to freeze-thaw cycle representation

Climate change can have substantial impacts on nitrogen runoff, which is a major cause of eutrophication, harmful algal blooms, and hypoxia in freshwaters and coastal regions. We examined responses of nitrate loading to climate change in the Upper Mississippi River Basin (UMRB) with an enhanced Soil and Water Assessment Tool with physically based Freeze-Thaw cycle representation (SWAT-FT), as compared with the original SWAT model that employs an empirical equation. Driven by future climate projections from five General Circulation Models (GCMs) from 1960 to 2099 under the Representative Concentrations Pathways (RCP) 8.5 scenario, we analyzed changes in riverine nitrate loadings, as well as terrestrial surface and subsurface contributions of the UMRB in the 21st century relative to the baseline period of 1960–1999. By the end of the 21st century, the original SWAT model predicted about a 50% increase in riverine nitrate loadings which is nearly twice as much as that estimated by SWAT-FT (ca. 25%). Such a large difference in projected nitrate changes can potentially mislead mitigation strategies that aim to reduce nitrogen runoff from the UMRB. Further analysis shows that the difference between the original SWAT model and SWAT-FT led to substantial discrepancies in the spatial distribution of surface and subsurface nitrate loadings in the UMRB. In general, SWAT-FT predicted more nitrate leaching for northwestern parts of the UMRB which are more sensitive to freeze-thaw cycle, mainly because SWAT-FT simulated less frequent frozen soils. This study highlights the importance of using physically based freeze-thaw cycle representation in water quality modeling. Design of future nitrogen runoff reduction strategies should include careful assessment of effects that land management has on the freeze-thaw cycles to provide reliable projection of water quality under climate change.

Wang, Qianfeng↗
DOE OSTI2022

A Decentralized Approach for Modeling Organized Convection Based on Thermal Populations on Microgrids

Abstract In this study, a spectral model for convective transport is coupled to a thermal population model on a two‐dimensional horizontal “microgrid,” covering the typical gridbox size of general circulation models. The goal is to explore new ways of representing impacts of spatial organization in cumulus cloud fields. The thermals are considered the smallest building block of convection, with thermal life cycle and movement represented through binomial functions. Thermals interact through two simple rules, reflecting pulsating growth and environmental deformation. Long‐lived thermal clusters thus form on the microgrid, exhibiting scale growth and spacing that represent simple forms of spatial organization and memory. Size distributions of cluster number are diagnosed from the microgrid through an online clustering algorithm, and provided as input to a spectral multiplume eddy‐diffusivity mass flux scheme. This yields a decentralized transport system, in that the thermal clusters acting as independent but interacting nodes that carry information about spatial structure. The main objectives of this study are (a) to seek proof of concept of this approach, and (b) to gain insight into impacts of spatial organization on convective transport. Single‐column model experiments demonstrate satisfactory skill in reproducing two observed cases of continental shallow convection. Metrics expressing self‐organization and spatial organization match well with large‐eddy simulation results. We find that in this coupled system, spatial organization impacts convective transport primarily through the scale break in the size distribution of cluster number. The rooting of saturated plumes in the subcloud mixed layer plays a key role in this process.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Dynamically Downscaled Hourly Future Weather Data with 12-km Resolution Covering Most of North America

This is an hourly future weather dataset for energy modeling applications. The dataset is primarily based on the output of a regional climate model (RCM), i.e., the Weather Research and Forecasting (WRF) model version 3.3.1. The WRF simulations are driven by the output of a general circulation model (GCM), i.e., the Community Climate System Model version 4 (CCSM4). This dataset is in the EPW format, which can be read or translated by more than 25 building energy modeling programs (e.g., EnergyPlus, ESP-r, and IESVE), energy system modeling programs (e.g., System Advisor Model (SAM)), indoor air quality analysis programs (e.g., CONTAM), and hygrothermal analysis programs (e.g., WUFI). It contains 13 weather variables, which are the Dry-Bulb Temperature, Dew Point Temperature, Relative Humidity, Atmospheric Pressure, Horizontal Infrared Radiation Intensity from Sky, Global Horizontal Irradiation, Direct Normal Irradiation, Diffuse Horizontal Irradiation, Wind Speed, Wind Direction, Sky Cover, Albedo, and Liquid Precipitation Depth. The weather data is created for two emissions scenarios: RCP4.5 and RCP8.5 and spans two 10-year time slices in the future: 2045 - 2054 and 2085 - 2094. It offers a spatial resolution of 12 km by 12 km with extensive coverage across most of North America. Due to the enormous size of the entire dataset, in the first stage of its distribution, we provide 20 years of future weather data for the centroid of each Public Use Microdata Area (PUMA), excluding Hawaii. PUMAs are non-overlapping, statistical geographic areas that partition each state or equivalent entity into geographic areas containing no fewer than 100,000 people each. The 2,378 PUMAs as a whole cover the entirety of the U.S. The weather data can be utilized alongside the large-scale energy analysis tools, ResStock and ComStock, developed by National Renewable Energy Laboratory, whose smallest resolution is at the PUMA scale. The data for RCP4.5 is still being processed and will be published soon.

Array↗
DOE OSTI2020

Reduced Complexity Model Intercomparison Project Phase 1: introduction and evaluation of global-mean temperature response

Reduced-complexity climate models (RCMs) are critical in the policy and decision making space, and are directly used within multiple Intergovernmental Panel on Climate Change (IPCC) reports to complement the results of more comprehensive Earth system models. To date, evaluation of RCMs has been limited to a few independent studies. Here we introduce a systematic evaluation of RCMs in the form of the Reduced Complexity Model Intercomparison Project (RCMIP). We expect RCMIP will extend over multiple phases, with Phase 1 being the first. In Phase 1, we focus on the RCMs' global-mean temperature responses, comparing them to observations, exploring the extent to which they emulate more complex models and considering how the relationship between temperature and cumulative emissions of CO2 varies across the RCMs. Our work uses experiments which mirror those found in the Coupled Model Intercomparison Project (CMIP), which focuses on complex Earth system and atmosphere–ocean general circulation models. Using both scenario-based and idealised experiments, we examine RCMs' global-mean temperature response under a range of forcings. We find that the RCMs can all reproduce the approximately 1°C of warming since pre-industrial times, with varying representations of natural variability, volcanic eruptions and aerosols. We also find that RCMs can emulate the global-mean temperature response of CMIP models to within a root-mean-square error of 0.2°C over a range of experiments. Furthermore, we find that, for the Representative Concentration Pathway (RCP) and Shared Socioeconomic Pathway (SSP)-based scenario pairs that share the same IPCC Fifth Assessment Report (AR5)-consistent stratospheric-adjusted radiative forcing, the RCMs indicate higher effective radiative forcings for the SSP-based scenarios and correspondingly higher temperatures when run with the same climate settings. In our idealised setup of RCMs with a climate sensitivity of 3°C, the difference for the ssp585–rcp85 pair by 2100 is around 0.23°C(±0.12°C) due to a difference in effective radiative forcings between the two scenarios. Phase 1 demonstrates the utility of RCMIP's open-source infrastructure, paving the way for further phases of RCMIP to build on the research presented here and deepen our understanding of RCMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Greenhouse-gas forced changes in the Atlantic meridional overturning circulation and related worldwide sea-level change

Abstract The effect of anthropogenic climate change in the ocean is challenging to project because atmosphere-ocean general circulation models (AOGCMs) respond differently to forcing. This study focuses on changes in the Atlantic Meridional Overturning Circulation (AMOC), ocean heat content ( $$\Delta$$ Δ OHC), and the spatial pattern of ocean dynamic sea level ( $$\Delta \zeta$$ Δ ζ ). We analyse experiments following the FAFMIP protocol, in which AOGCMs are forced at the ocean surface with standardised heat, freshwater and momentum flux perturbations, typical of those produced by doubling $$\hbox {CO}_{{2}}$$ CO 2 . Using two new heat-flux-forced experiments, we find that the AMOC weakening is mainly caused by and linearly related to the North Atlantic heat flux perturbation, and further weakened by a positive coupled heat flux feedback. The quantitative relationships are model-dependent, but few models show significant AMOC change due to freshwater or momentum forcing, or to heat flux forcing outside the North Atlantic. AMOC decline causes warming at the South Atlantic-Southern Ocean interface. It does not strongly affect the global-mean vertical distribution of $$\Delta$$ Δ OHC, which is dominated by the Southern Ocean. AMOC decline strongly affects $$\Delta \zeta$$ Δ ζ in the North Atlantic, with smaller effects in the Southern Ocean and North Pacific. The ensemble-mean $$\Delta \zeta$$ Δ ζ and $$\Delta$$ Δ OHC patterns are mostly attributable to the heat added by the flux perturbation, with smaller effects from ocean heat and salinity redistribution. The ensemble spread, on the other hand, is largely due to redistribution, with pronounced disagreement among the AOGCMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

pyTCR: A tropical cyclone rainfall model for python

pyTCR is a climatology software package developed in the Python programming language. It integrates the capabilities of several legacy physical models and increases computational efficiency to allow rapid estimation of tropical cyclone (TC) rainfall consistent with the large-scale environment. Specifically, pyTCR implements a horizontally distributed and vertically integrated model [Zhu et al., 2013] for simulating rainfall driven by TCs. Along storm tracks, rainfall is estimated by computing the cross-boundary-layer, upward water vapor transport caused by different mechanisms including frictional convergence, vortex stretching, large-scale baroclinic effect (i.e., wind shear), topographic forcing, and radiative cooling [Lu et al., 2018]. The package provides essential functionalities for modeling and interpreting spatio-temporal TC rainfall data. pyTCR requires a limited number of model input parameters, making it a convenient and useful tool for analyzing rainfall mechanisms driven by TCs. To sample rare (most intense) rainfall events that are often of great societal interest, pyTCR adapts and leverages outputs from a statistical-dynamical TC downscaling model [Lin et al., 2023] capable of rapidly generating a large number of synthetic TCs given a certain climate. As a result, pyTCR significantly reduces computational effort and improves the efficiency in capturing extreme TC rainfall events at the tail of the distributions from limited datasets. Furthermore, the TC downscaling model is forced entirely by large-scale environmental conditions from reanalysis data or coupled General Circulation Models (GCMs), simplifying the projection of TC-induced rainfall and wind speed under future climate using pyTCR. Finally, pyTCR can be coupled with hydrological and wind models to assess risks associated with independent and compound events (e.g., storm surges and freshwater flooding).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Towards multiscale modeling of ocean surface turbulent mixing using coupled MPAS-Ocean v6.3 and PALM v5.0

A multiscale modeling approach for studying the ocean surface turbulent mixing is explored by coupling an ocean general circulation model (GCM) MPAS-Ocean with the Parallelized Large Eddy Simulation Model (PALM). The coupling approach is similar to the superparameterization approach that has been used to represent the effects of deep convection in atmospheric GCMs. However, the focus of this multiscale modeling approach is on the small-scale turbulent mixing and their interactions with the larger-scale processes in the ocean, so that a more flexible coupling strategy is used. To reduce the computational cost, a customized version of PALM is ported on the general-purpose graphics processing unit (GPU) with OpenACC, achieving 10–16 times overall speedup as compared to running on a single CPU. Even with the GPU-acceleration technique, a superparameterization-like approach to represent the ocean surface turbulent mixing in GCMs using embedded high fidelity and three-dimensional large eddy simulations (LESs) over the global ocean is still computationally intensive and infeasible for long simulations. However, running PALM regionally on selected MPAS-Ocean grid cells is shown to be a promising approach moving forward. The flexible coupling between MPAS-Ocean and PALM allows further exploration of the interactions between the ocean surface turbulent mixing and larger-scale processes, as well as future development and improvement of ocean surface turbulent mixing parameterizations for GCMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Combining regional mesh refinement with vertically enhanced physics to target marine stratocumulus biases as demonstrated in the Energy Exascale Earth System Model version 1

Abstract. In this paper we develop a novel framework aimed to significantly reduce biases related to marine stratocumulus clouds in general circulation models (GCMs) while circumventing excessive computational cost requirements. Our strategy is to increase the horizontal resolution using a regionally refined mesh (RRM) over our region of interest in addition to using the Framework for Improvement by Vertical Enhancement (FIVE) to increase the vertical resolution only for specific physical processes that are important for stratocumulus. We apply the RRM off the coast of Peru in the southeastern Pacific, a region that climatologically contains the most marine stratocumulus in the subtropics. We find that our new modeling framework is able to replicate the results of our high-resolution benchmark simulation with much fidelity, while reducing the computational cost by several orders of magnitude. In addition, this framework is able to greatly reduce the long-standing biases associated with marine stratocumulus in GCMs when compared to the standard-resolution control simulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A protocol and analysis of year-long simulations of global storm-resolving models and beyond

We propose a protocol to evaluate and analyze year-long simulations of global storm-resolving models (GSRMs). The proposed protocol complements an earlier 40-day simulation protocol under the DYAMOND (DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains) project to allow the analysis of the seasonal cycle and associated climatic relevant phenomena. This intercomparison aims to reveal how GSRMs, which can simulate mesoscale convective systems (MCSs) in the global domain, reproduce atmospheric large-scale structures related to convection beyond month-long simulations. The intercomparison for one-year simulations is conducted by either atmosphere-only models or atmosphere–ocean coupled models with atmospheric horizontal mesh sizes less than 5 km. We recommend the continuous four seasons from March 2020 to February 2021 as a target period for the intercomparison but with options for many groups to join more flexibly. The output variables are collected at 0.25° resolution, and archives of a small set of native grid variables are encouraged to analyze tropical cyclones and MCSs. Through the proposed global storm-resolving simulation, we will evaluate the climatological distributions of the atmospheric large-scale circulations, such as the Intertropical Convergence Zone (ITCZ), monsoon, midlatitude jets, their time evolution, and the upscale impacts on them. We present sample analyses from a one-year simulation using the 3.5 km mesh Nonhydrostatic Icosahedral Atmospheric Model (NICAM), revealing the realistic zonal contrast of tropical precipitation, no double ITCZ structure, the reasonable midlatitude jet position and intensity but a weak bias of storm track activities, and a warm bias over the Eurasia during boreal winter. We also clarify the cross-scale interaction, such as the effects of cold pools on mean precipitation over the Maritime Continent through the precipitation diurnal cycle and the effects of resolved gravity waves on midlatitude mean flows. The proposed one-year simulation protocol is referred to as the “Sendai Protocol.” This protocol is not unique or definite for evaluating GSRMs; we prospect a hierarchical set of experiments from short-term to multi-year simulations as GSRM intercomparisons.

54 ENVIRONMENTAL SCIENCES↗