skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • Advancing precipitation prediction using a new-generation storm-resolving model framework – SIMA-MPAS (V1.0): a case study over the western United States
Title: Advancing precipitation prediction using a new-generation storm-resolving model framework – SIMA-MPAS (V1.0): a case study over the western United States
Abstract. Global climate models (GCMs) have advanced in many ways ascomputing power has allowed more complexity and finer resolutions. As GCMsreach storm-resolving scales, they need to be able to produce realisticprecipitation intensity, duration, and frequency at fine scales withconsideration of scale-aware parameterization. This study uses astate-of-the-art storm-resolving GCM with a nonhydrostatic dynamical core – theModel for Prediction Across Scales (MPAS), incorporated in the atmosphericcomponent (Community Atmosphere Model, CAM) of the open-source CommunityEarth System Model (CESM), within the System for Integrated Modeling of theAtmosphere (SIMA) framework (referred to as SIMA-MPAS). At uniform coarse (here, at 120 km) gridresolution, the SIMA-MPAS configuration is comparable to the standardhydrostatic CESM (with a finite-volume (FV) dynamical core) with reasonableenergy and mass conservation on climatological timescales. With thecomparable energy and mass balance performance between CAM-FV (workhorse dynamical core) and SIMA-MPAS (newly developed dynamical core), it gives confidence inSIMA-MPAS's applications at a finer resolution. To evaluate this, we focuson how the SIMA-MPAS model performs when reaching a storm-resolving scale at3 km. To do this efficiently, we compose a case study using a SIMA-MPASvariable-resolution configuration with a refined mesh of 3 km covering thewestern USA and 60 km over the rest of the globe. We evaluated the modelperformance using satellite and station-based gridded observations withcomparison to a traditional regional climate model (WRF, the WeatherResearch and Forecasting model). Our results show realistic representationsof precipitation over the refined complex terrains temporally and spatially.Along with much improved near-surface temperature, realistic topography, andland–air interactions, we also demonstrate significantly enhanced snowpackdistributions. This work illustrates that the global SIMA-MPAS atstorm-resolving resolution can produce much more realistic regional climatevariability, fine-scale features, and extremes to advance both climate andweather studies. This next-generation storm-resolving model could ultimatelybridge large-scale forcing constraints and better inform climate impactsand weather predictions across scales.  more » « less
Award ID(s):
2005137
PAR ID:
10434395
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Geoscientific Model Development
Volume:
15
Issue:
21
ISSN:
1991-9603
Page Range / eLocation ID:
8135 to 8151
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, Journal of Geophysical Research: Atmospheres)
    Abstract Global Storm Resolving Models (GSRMs) provide a way to understand weather and climate events across scales for better‐informed climate impacts. In this work, we apply the recently developed and validated CAM (Community Atmosphere Model)—MPAS (Model for Prediction Across Scales) modeling framework, based on the open‐source Community Earth System Model (CESM2), to examine the tropical convection features at the storm resolving scale over the Maritime Continent region at 3 km horizontal spacing. We target two global numerical experiments during the winter season of 2018 for comparison with observation in the region. We focus on the investigation of the representations of the convective systems, precipitation statistics, and tropical cyclone behaviors. We found that regional‐refined experiments show more accurate precipitation distributions, diurnal cycles, and better agreement with observations for tropical cyclone features in terms of intensity and strength statistics. We expect the exploration of this work will further advance the development and use of the storm‐resolving model in precipitation predictions across scales. 
    more » « less
  2. ; ; ; ; ; (, Journal of Advances in Modeling Earth Systems)
    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
  3. ; ; ; ; ; ; ; ; ; ; et al (, Journal of Advances in Modeling Earth Systems)
    Abstract Global modeling of aerosol‐particle number and size is important for understanding aerosol effects on Earth's climate and air quality. Fine‐resolution global models are desirable for representing nonlinear aerosol‐microphysical processes, their nonlinear interactions with dynamics and chemistry, and spatial heterogeneity. However, aerosol‐microphysical simulations are computationally demanding, which can limit the achievable global horizontal resolution. Here, we present the first coupling of the TwO‐Moment Aerosol Sectional (TOMAS) microphysics scheme with the High‐Performance configuration of the GEOS‐Chem model of atmospheric composition (GCHP), a coupling termed GCHP‐TOMAS. GCHP's architecture allows massively parallel GCHP‐TOMAS simulations including on the cloud, using hundreds of computing cores, faster runtimes, more memory, and finer global horizontal resolution (e.g., 25 km × 25 km, 7.8 × 105model columns) versus the previous single‐node capability of GEOS‐Chem‐TOMAS (tens of cores, 200 km × 250 km, 1.3 × 104model columns). GCHP‐TOMAS runtimes have near‐ideal scalability with computing‐core number. Simulated global‐mean number concentrations increase (dominated by free‐tropospheric over‐ocean sub‐10‐nm‐diameter particles) toward finer GCHP‐TOMAS horizontal resolution. Increasing the horizontal resolution from 200 km × 200–50 km × 50 km increases the global monthly mean free‐tropospheric total particle number by 18.5%, and over‐ocean sub‐10‐nm‐diameter particles by 39.8% at 4‐km altitude. With a cascade of contributing factors, free‐tropospheric particle‐precursor concentrations increase (32.6% at 4‐km altitude) with resolution, promoting new‐particle formation and growth that outweigh coagulation changes. These nonlinear effects have the potential to revise current understanding of processes controlling global aerosol number and aerosol impacts on Earth's climate and air quality. 
    more » « less
  4. ; ; ; ; ; ; (, The Cryosphere)
    Abstract. Earth system models are essential tools for understandingthe impacts of a warming world, particularly on the contribution of polarice sheets to sea level change. However, current models lack full couplingof the ice sheets to the ocean and are typically run at a coarse resolution(1∘ grid spacing or coarser). Coarse spatial resolution isparticularly a problem over Antarctica, where sub-grid-scale orography iswell-known to influence precipitation fields, and glacier models requirehigh-resolution atmospheric inputs. This resolution limitation has beenpartially addressed by regional climate models (RCMs), which must be forcedat their lateral and ocean surface boundaries by (usually coarser) globalatmospheric datasets, However, RCMs fail to capture the two-way couplingbetween the regional domain and the global climate system. Conversely,running high-spatial-resolution models globally is computationallyexpensive and can produce vast amounts of data. Alternatively, variable-resolution grids can retain the benefits of highresolution over a specified domain without the computational costs ofrunning at a high resolution globally. Here we evaluate a historicalsimulation of the Community Earth System Model version 2 (CESM2)implementing the spectral element (SE) numerical dynamical core (VR-CESM2)with an enhanced-horizontal-resolution (0.25∘) grid over theAntarctic Ice Sheet and the surrounding Southern Ocean; the rest of theglobal domain is on the standard 1∘ grid. We compare it to1∘ model runs of CESM2 using both the SE dynamical core and thestandard finite-volume (FV) dynamical core, both with identical physics andforcing, including prescribed sea surface temperatures (SSTs) and sea ice concentrations fromobservations. Our evaluation reveals both improvements and degradations inVR-CESM2 performance relative to the 1∘ CESM2. Surface massbalance estimates are slightly higher but within 1 standard deviation ofthe ensemble mean, except for over the Antarctic Peninsula, which isimpacted by better-resolved surface topography. Temperature and windestimates are improved over both the near surface and aloft, although theoverall correction of a cold bias (within the 1∘ CESM2 runs) hasresulted in temperatures which are too high over the interior of the icesheet. The major degradations include the enhancement of surface melt aswell as excessive cloud liquid water over the ocean, with resultant impactson the surface radiation budget. Despite these changes, VR-CESM2 is avaluable tool for the analysis of precipitation and surface mass balanceand thus constraining estimates of sea level rise associated with theAntarctic Ice Sheet. 
    more » « less
  5. ; ; ; ; ; ; (, Geophysical Research Letters)
    Abstract Ocean warming is a key factor impacting future changes in climate. Here we investigate vertical structure changes in globally averaged ocean heat content (OHC) in high‐ (HR) and low‐resolution (LR) future climate simulations with the Community Earth System Model (CESM). Compared with observation‐based estimates, the simulated OHC anomalies in the upper 700 and 2,000 m during 1960–2020 are more realistic in CESM‐HR than ‐LR. Under RCP8.5 scenario, the net surface heat into the ocean is very similar in CESM‐HR and ‐LR. However, CESM‐HR has a larger increase in OHC in the upper 250 m compared to CESM‐LR, but a smaller increase below 250 m. This difference can be traced to differences in eddy‐induced vertical heat transport between CESM‐HR and ‐LR in the historical period. Moreover, our results suggest that with the same heat input, upper‐ocean warming is likely to be underestimated by most non‐eddy‐resolving climate models. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email