An official website of the United States government
Abstract We design a new strategy to load‐balance high‐intensity sub‐grid atmospheric physics calculations restricted to a small fraction of a global climate simulation's domain. We show why the current parallel load balancing infrastructure of Community Earth System Model (CESM) and Energy Exascale Earth Model (E3SM) cannot efficiently handle this scenario at large core counts. As an example, we study an unusual configuration of the E3SM Multiscale Modeling Framework (MMF) that embeds a binary mixture of two separate cloud‐resolving model grid structures that is attractive for low cloud feedback studies. Less than a third of the planet uses high‐resolution (MMF‐HR; sub‐km horizontal grid spacing) relative to standard low‐resolution (MMF‐LR) cloud superparameterization elsewhere. To enable MMF runs with Multi‐Domain cloud resolving models (CRMs), our load balancing theory predicts the most efficient computational scale as a function of the high‐intensity work's relative overhead and its fractional coverage. The scheme successfully maximizes model throughput and minimizes model cost relative to precursor infrastructure, effectively by devoting the vast majority of the processor pool to operate on the few high‐intensity (and rate‐limiting) high‐resolution (HR) grid columns. Two examples prove the concept, showing that minor artifacts can be introduced near the HR/low‐resolution CRM grid transition boundary on idealized aquaplanets, but are minimal in operationally relevant real‐geography settings. As intended, within the high (low) resolution area, our Multi‐Domain CRM simulations exhibit cloud fraction and shortwave reflection convergent to standard baseline tests that use globally homogenous MMF‐LR and MMF‐HR. We suggest this approach can open up a range of creative multi‐resolution climate experiments without requiring unduly large allocations of computational resources.
more »
« less
On the Dependence of Simulated Convection on Domain Size in CRMs
Abstract We present a heuristic model to explain the suppression of deep convection in convection‐resolving models (CRMs) with a small number of grid columns, such as those used in super‐parameterized or multi‐scale modeling framework (MMF) general circulation models (GCM) of the atmosphere. Domains with few grid columns require greater instability to sustain convection because they force a large convective fraction, driving strong compensating subsidence warming. Updraft dilution, which is stronger for reduced horizontal grid spacing, enhances this effect. Thus, suppression of deep convection in CRMs with few grid columns can be reduced by increasing grid spacing. Radiative‐convective equilibrium simulations using standalone CRM simulations with the System for Atmospheric Modeling (SAM) and using GCM‐coupled CRM simulations with the Energy Exascale Earth System Model (E3SM)‐MMF confirm the heuristic model results.
more »
« less
- Award ID(s):
- 1743753
- PAR ID:
- 10588249
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Advances in Modeling Earth Systems
- Volume:
- 17
- Issue:
- 3
- ISSN:
- 1942-2466
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The Walker circulation connects the regions with deep atmospheric convection in the western tropical Pacific to the shallow‐convection, tropospheric subsidence, and stratocumulus cloud decks of the eastern Pacific. The purpose of this study is to better understand the multi‐scale interactions between the Walker circulation, cloud systems, and interactive radiation. To do this we simulate a mock‐Walker Circulation with a full‐physics general circulation model using idealized boundary conditions. Our experiments use a doubly‐periodic domain with grid‐spacing of 1, 2, 25, and 100 km. We thus span the range from General Circulation Models (GCMs) to Cloud‐system Resolving Models (CRMs). Our model is derived from the Geophysical Fluid Dynamics Laboratory atmospheric GCM (AM4.0). We find substantial differences in the mock‐Walker circulation simulated by our GCM‐like and CRM‐like experiments. The CRM‐like experiments have more upper level clouds, stronger overturning circulations, and less precipitation. The GCM‐like experiments have a low‐level cloud fraction that is up to 20% larger. These differences leads to opposite atmospheric responses to changes in the longwave cloud radiative effect (LWCRE). Active LWCRE leads to increased precipitation for our GCMs, but decreased precipitation for our CRMs. The LWCRE leads to a narrower rising branch of the circulation and substantially increases the fraction of precipitation from the large‐scale cloud parameterization. This work demonstrates that a mock‐Walker circulation is a useful generalization of radiative convective equilibrium that includes a large‐scale circulation.more » « less
-
Abstract A severe derecho impacted the Midwestern United States on 10 August 2020, causing over $12 billion (U.S. dollars) in damage, and producing peak winds estimated at 63 m s−1, with the worst impacts in Iowa. The event was not forecast well by operational forecasters, nor even by operational and quasi-operational convection-allowing models. In the present study, nine simulations are performed using the Limited Area Model version of the Finite-Volume-Cubed-Sphere model (FV3-LAM) with three horizontal grid spacings and two physics suites. In addition, when a prototype of the Rapid Refresh Forecast System (RRFS) physics is used, sensitivity tests are performed to examine the impact of using the Grell–Freitas (GF) convective scheme. Several unusual results are obtained. With both the RRFS (not using GF) and Global Forecast System (GFS) physics suites, simulations using relatively coarse 13- and 25-km horizontal grid spacing do a much better job of showing an organized convective system in Iowa during the daylight hours of 10 August than the 3-km grid spacing runs. In addition, the RRFS run with 25-km grid spacing becomes much worse when the GF convective scheme is used. The 3-km RRFS run that does not use the GF scheme develops spurious nocturnal convection the night before the derecho, removing instability and preventing the derecho from being simulated at all. When GF is used, the spurious storms are removed and an excellent forecast is obtained with an intense bowing echo, exceptionally strong cold pool, and roughly 50 m s−1surface wind gusts.more » « less
-
Abstract We examine the influence of convective organization on extreme tropical precipitation events using model simulation data from the Radiative‐Convective Equilibrium Model Intercomparison Project (RCEMIP). At a given SST, simulations with convective organization have more intense precipitation extremes than those without it at all scales, including instantaneous precipitation at the grid resolution (3 km). Across large‐domain simulations with convective organization, models with explicit convection exhibit better agreement in the response of extreme precipitation rates to warming than those with parameterized convection. Among models with explicit convection, deviations from the Clausius‐Clapeyron scaling of precipitation extremes with warming are correlated with changes in organization, especially on large spatiotemporal scales. Though the RCEMIP ensemble is nearly evenly split between CRMs which become more and less organized with warming, most of the models which show increased organization with warming also allow super‐CC scaling of precipitation extremes. We also apply an established precipitation extremes scaling to understand changes in the extreme condensation events leading to extreme precipitation. Increased organization leads to greater increases in precipitation extremes by enhancing both the dynamic and implied efficiency contributions. We link these contributions to environmental variables modified by the presence of organization and suggest that increases in moisture in the aggregated region may be responsible for enhancing both convective updraft area fraction and precipitation efficiency. By leveraging a controlled intercomparison of models with both explicit and parameterized convection, this work provides strong evidence for the amplification of tropical precipitation extremes and their response to warming by convective organization.more » « less
-
Abstract. The radiative–convective equilibrium (RCE) model intercomparison project (RCEMIP) leveraged the simplicity of RCE to focus attention on moist convective processes and their interactions with radiation and circulation across a wide range of model types including cloud-resolving models (CRMs), general circulation models (GCMs), single-column models, global cloud-resolving models, and large-eddy simulations. While several robust results emerged across the spectrum of models that participated in the first phase of RCEMIP (RCEMIP-I), two points that stand out are (1) the strikingly large diversity in simulated climate states and (2) the strong imprint of convective self-aggregation on the climate state. However, the lack of consensus in the structure of self-aggregation and its response to warming is a barrier to understanding. Gaining a deeper understanding of convective aggregation and tropical climate will require reducing the degrees of freedom with which convection can vary. Therefore, we propose phase II of RCEMIP (RCEMIP-II) that utilizes a prescribed sinusoidal sea surface temperature (SST) pattern to provide a constraint on the structure of convection and move one critical step up the model hierarchy. This so-called “mock-Walker” configuration generates features that resemble observed tropical circulations. The specification of the mock-Walker protocol for RCEMIP-II is described, along with example results from one CRM and one GCM. RCEMIP-II will consist of five required simulations: three simulations with the same three mean SSTs as in RCEMIP-I but with an SST gradient and two additional simulations at one of the mean SSTs with different values of the SST gradients. We also test the sensitivity to the imposed SST gradient and the domain size. Under weak SST gradients, unforced self-aggregation emerges across the entire domain, similar to what was found in RCEMIP. As the SST gradient increases, the convective region narrows and is more confined to the warmest SSTs. At warmer mean SSTs and stronger SST gradients, low-frequency variability in the convective aggregation emerges, suggesting that simulations of at least 200 d may be needed to achieve robust equilibrium statistics in this configuration. Simulations with different domain sizes generally have similar mean statistics and convective structures, depending on the value of the SST gradient. The prescribed SST boundary condition is the only difference in the set-up between RCEMIP-II and RCEMIP-I, which enables comparison between the two; however, we also welcome participation in RCEMIP-II from models that did not participate in RCEMIP-I.more » « less
