Years-long global storm-resolving simulations of the present and warmed climates are conducted with 3.25km horizontal grid spacing. We focus on wintertime precipitation in the coastal western United States, a notably water-sensitive region with complex topography. The model generates more realistic orographic precipitation compared with a coarser-resolution model having the same dynamical core. In response to uniform sea surface warming, the increase in extreme precipitation rates together with the better resolved orographic precipitation lead to persistence of snowpack in parts of the coastal western United States. In contrast, snowpack in the coarser-resolution model is largely eliminated in the warmed climate. Whether snowpack persists influences the regional surface energy budget due to the snow-albedo feedback. The results highlight the importance of resolving orographic precipitation and the large-scale circulation response to warming in a consistent framework.

Observations reveal a clear difference in the intensity of deep convection over tropical land and ocean. This observed land-ocean contrast provides a natural benchmark for evaluating the fidelity of global storm-resolving models (GSRMs; global models with horizontal resolution on the order of kilometers), and GSRMs provide a potentially valuable tool for probing unresolved scientific questions about the origin of the observed land-ocean contrast. However, land-ocean differences in convective intensity have received relatively little attention in GSRM research. Here, we show that the strength of the land-ocean contrast simulated by GSRMs is strongly sensitive to details of GSRM implementations, and not clearly governed by any of several hypothesized drivers of the observed land-ocean contrast. We first examine DYAMOND Summer GSRM simulations, and show that only a subset produce a clear land-ocean contrast in the frequency of strong updrafts. We then show that the use of a sub-grid shallow convection scheme can determine whether or not the GSRM X-SHiELD produces a clear land-ocean contrast. Finally, we show that three hypothesized drivers of the observed land-ocean contrast all fail to explain why a land-ocean contrast is present in X-SHiELD simulations with sub-grid shallow convection disabled. These results provide encouraging evidence that GSRMs can mimic the observed land-ocean convective intensity contrast. However, they also show that their ability to do so can be sensitive to uncertain sub-grid parameterizations, and suggest that existing theory may not fully capture drivers of the land-ocean contrast simulated by some GSRMs.

Global storm-resolving model (GSRM) simulations (kilometer-scale horizontal resolution) of the atmosphere can capture the interaction between the scales of deep cumulus convection and the large-scale dynamics and thermodynamic properties of the atmosphere. Here, we assess the vertical structure of tropical temperature change in the GSRM X-SHiELD, developed by the Geophysical Fluid Dynamics Laboratory, perturbed by a uniform sea surface temperature (SST) warming and/or increased CO2 concentration. The simulated response to SST warming shows weakly amplified warming from the surface through the mid-troposphere before increasing to a factor of about 2.5 near the tropopause. This combination of muted warming in the mid-troposphere and amplified warming aloft is within the range of CMIP6 models at individual pressure levels but, taken together, is distinctive behavior. The response to CO2 increase with unchanged SST is an approximately vertically uniform warming, comparable to CMIP6 models, and is linearly additive with the SST-induced warming in X-SHiELD.

We investigate the representation of individual supercells and intriguing tornado-like vortices in a simplified, locally-refined global atmosphere model. The model, featuring grid stretching, can locally enhance the model resolution and economically reach the cloud-resolving scale. Given an unstable sheared environment, the model can simulate supercells realistically, with a near-ground vortex and funnel cloud at the center of a rotating updraft reminiscent of a tornado. An analysis of the vorticity budget suggests that the updraft core of the supercell tilts environmental horizontal vorticity into the tornado-like vortex. The updraft also acts to amplify the vortex through vertical stretching. Results suggest that the simulated vortex is dynamically similar to observed tornadoes and modeling studies at much higher horizontal resolution.

Recent advances have allowed for integration of global storm resolving models (GSRMs) to a timescale of several years. These short simulations are sufficient for studying characteristics and statistics of short- and small-scale phenomena; however, it is questionable what we can learn from these integrations about the large-scale climate response to perturbations. To address this question, we use the response of X-SHiELD (a GSRM) to uniform SST warming and CO$_2$ increase in a two-year integration and compare it to similar CMIP6 experiments. Specifically, we assess the statistical meaning of having two years in one model outside the spread of another model or model ensemble. This is of particular interest because X-SHiELD shows a distinct response of the global mean precipitation to uniform warming, and the northern hemisphere jet shift response to isolated CO$_2$ increase. We use the CMIP6 models to estimate the probability of two years in one model being more than one standard deviation away from another model (ensemble) mean, knowing the mean of two models. For example, if two years in one model are more than one standard deviation away from the other model’s mean, we find that the chances for these models’ means to be within one standard deviation are $\sim 25\%$. We find that for some large-scale metrics, there is an important base-state dependence that, when taken into account, can qualitatively change the interpretation of the results. We note that a year-to-year comparison is physically meaningful due to the use of prescribed sea-surface-temperature simulations.

We present the System for High-resolution prediction on Earth-to-Local Domains (SHiELD), an atmosphere model coupling the nonhydrostatic FV3 Dynamical Core to a physics suite originally taken from the Global Forecast System. SHiELD is designed to demonstrate new capabilities within its components, explore new model applications, and to answer scientific questions through these new functionalities. A variety of configurations are presented, including short-to-medium-range and subseasonal-to-seasonal (S2S) prediction, global-to-regional convective-scale hurricane and contiguous US precipitation forecasts, and global cloud-resolving modeling. Advances within SHiELD can be seamlessly transitioned into other Unified Forecast System (UFS) or FV3-based models, including operational implementations of the UFS. Continued development of SHiELD has shown improvement upon existing models. The flagship 13-km SHiELD demonstrates steadily improved large-scale prediction skill and precipitation prediction skill. SHiELD and the coarser-resolution S-SHiELD demonstrate a superior diurnal cycle compared to existing climate models; the latter also demonstrates 28 days of useful prediction skill for the Madden-Julian Oscillation. The global-to-regional nested configurations T-SHiELD (tropical Atlantic) and C-SHiELD (contiguous United States) shows significant improvement in hurricane structure from a new tracer advection scheme and promise for medium-range prediction of convective storms, respectively.

Intense convection (updrafts exceeding 10 m∙s-1) plays an essential role in severe weather and Earth’s energy balance. Despite its importance, how the global pattern of intense convection changes in response to warmed climates remains unclear, as simulations from traditional climate models are too coarse to simulate intense convection. Here we take advantage of a kilometer-scale global storm resolving model and conduct year-long simulations of a control run, forced by analyzed sea surface temperature (SST), and one with a 4-K increase in SST for comparison. Comparisons show that the increased SST enhances the frequency of intense convection globally with large spatial and seasonal variations. Increases in the intense convection frequency do not necessarily reflect increases in convective available potential energy (CAPE). Results are also compared with traditional climate model projections. Changes in the spatial pattern of intense convection are associated with changes in planetary circulation.

We present the characteristics of rotating convective updrafts in the 2021 version of GFDL’s Experimental System for High-resolution prediction on Earth-to-Local Domains (X-SHiELD), a kilometer scale global storm resolving model (GSRM). Rotation is quantified using 2–5 km Updraft Helicity (UH) in a year-long integration forced by analyzed SSTs. Updrafts with UH magnitudes above 50 m2/s2; are common over the mid-latitude continents, especially in the warm seasons where they are associated with severe weather but are also common over most tropical ocean basins. In nearly all areas cyclonically rotating convection dominates, with larger UH values increasingly preferring cyclonic rotation. The ratio of cyclonic to anticyclonic updrafts is largest in the subtropical and mid-latitude oceans and is slightly lower over mid-latitude continents. The ratio of cyclonic to anticyclonic updrafts can be substantively explained by the mean storm-relative helicity (SRH) in convective regions, indicating the importance for environmental controls on the sense of storm rotation, although internal storm dynamics also plays a role in the generation of anticyclonic updrafts.