We present the System for High‐resolution prediction on Earth‐to‐Local Domains (SHiELD), an atmosphere model developed by the Geophysical Fluid Dynamics Laboratory (GFDL) 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 prediction, global‐to‐regional convective‐scale hurricane and contiguous U.S. precipitation forecasts, and global cloud‐resolving modeling. Advances within SHiELD can be seamlessly transitioned into other Unified Forecast System or FV3‐based models, including operational implementations of the Unified Forecast System. 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) show significant improvement in hurricane structure from a new tracer advection scheme and promise for medium‐range prediction of convective storms.

We describe the third version of the Geophysical Fluid Dynamics Laboratory cloud microphysics scheme (GFDL MP v3) implemented in the System for High‐resolution prediction on Earth‐to‐Local Domains (SHiELD). Compared to the GFDL MP v2, the GFDL MP v3 is entirely reorganized, optimized, and modularized into functions. The particle size distribution (PSD) of all hydrometeor categories is redefined to better mimic observations, and the cloud droplet number concentration (CDNC) is calculated from the Modern‐Era Retrospective analysis for Research and Applications, version 2 (MERRA2) aerosol data. In addition, the GFDL MP has been redesigned so all processes use the redefined PSD to ensure overall consistency and easily permit the introduction of new PSDs and microphysical processes. A year's worth of global 13‐km, 10‐day weather forecasts were performed with the new GFDL MP. Compared to the GFDL MP v2, the GFDL MP v3 significantly improves SHiELD's predictions of geopotential height, air temperature, and specific humidity in the Troposphere, as well as the high, middle and total cloud fractions and the liquid water path. The predictions are improved even further by the use of reanalysis aerosol data to calculate CDNC, and also by using the more realistic PSD available in GFDL MP v3. However, the upgrade of the GFDL MP shows little impact on the precipitation prediction. Degradations caused by the new scheme are discussed and provide a guide for future GFDL MP development.

High‐resolution atmospheric models are powerful tools for hurricane track and intensity predictions. Although using high resolution contributes to better representation of hurricane structure and intensity, its value in the prediction of steering flow and storm tracks is uncertain. Here we present experiments suggesting that biases in the predicted North Atlantic hurricane tracks in a high‐resolution (approximately 3 km grid‐spacing) model originates from the model's explicit simulation of deep convection. Differing behavior of explicit convection leads to changes in the synoptic‐scale pattern and thereby to the steering flow. Our results suggest that optimizing small‐scale convection activity, for example, through the model's horizontal advection scheme, can lead to significantly improved hurricane track prediction (∼10% reduction of mean track error) at lead times beyond 72 hr. This work calls attention to the behavior of explicit convection in high‐resolution models, and its often overlooked role in affecting larger‐scale circulations and hurricane track prediction.

In this study, the authors numerically investigate the response of an axisymmetric tropical cyclone (TC) vortex to the vertical fluxes of momentum, heat, and moisture induced by roll vortices (rolls) in the boundary layer. To represent the vertical fluxes induced by rolls, a two-dimensional high-resolution Single-Grid Roll-Resolving Model (SRM) is embedded at multiple horizontal grid points in the mesoscale COAMPS for Tropical Cyclones (COAMPS-TC) model domain. Idealized experiments are conducted with the SRM embedded within 3 times the radius of maximum wind of an axisymmetric TC. The results indicate that the rolls induce changes in the boundary layer wind distribution and cause a moderate (approximately 15%) increase in the TC intensification rate by increasing the boundary layer convergence in the eyewall region and induce more active eyewall convection. The numerical experiments also suggest that the roll-induced tangential momentum flux is most important in contributing to the TC intensification process, and the rolls generated at different radii (within the range considered in this study) all have positive contributions. The results are not qualitatively impacted by the initial TC vortex or the setup of the vertical diffusivity in COAMPS-TC.

We demonstrate that two‐way nesting significantly improves the structure of simulated hurricane in an Atmospheric General Circulation Model (AGCM). Two sets of 30‐day hindcast experiments are conducted, one with the global‐uniform‐resolution (approximately 25‐km nominal horizontal resolution) and the other with a regionally refined two‐way nest (approximately 8‐km over the tropical North Atlantic). The increase in the horizontal resolution on the nested grid improves the representation of storm intensity and intensification rate. When normalized by the radius of maximum wind (RMW), composite hurricane structures are generally similar in both simulations and compare well to observations. However, the hurricanes in the globally‐uniform configuration have much larger RMWs than observed, while those in the two‐way‐nested configuration have more realistic RMWs. We also find that the representation of the RMW has a critical impact on the simulation of inertial stability and boundary‐layer convergence in the inner‐core region. The more realistic inner‐core size (indicated by RMW) and structure are possible reasons for the improved intensification rates in the two‐way‐nested configuration.