Search for: All records

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.

 

Subkilometer processes are critical to the physics of aerosol‐cloud interaction (ACI) but have been dependent on parameterizations in global model simulations. We thus report the strength of ACI in the Ultra‐Parameterized Community Atmosphere Model (UPCAM), a multiscale climate model that uses coarse exterior resolution to embed explicit cloud‐resolving models with enough resolution (250 m horizontal, 20 m vertical) to quasi‐resolve subkilometer eddies. To investigate the impact on ACIs, UPCAM's simulations are compared to a coarser multiscale model with 4 km horizontal resolution. UPCAM produces cloud droplet number concentrations (N_d) and cloud liquid water path (LWP) values that are higher than the coarser model but equally plausible compared to observations. Our analysis focuses on the Northern Hemisphere (20–50°N) oceans, where historical aerosol increases have been largest. We find similarities in the overall radiative forcing from ACIs in the two models, but this belies fundamental underlying differences. The radiative forcing from increases in LWP is weaker in UPCAM, whereas the forcing from increases inN_dis larger. Surprisingly, the weaker LWP increase is not due to a weaker increase in LWP in raining clouds, but a combination of weaker increase in LWP in nonraining clouds and a smaller fraction of raining clouds in UPCAM. The implication is that as global modeling moves toward finer than storm‐resolving grids, nuanced model validation of ACI statistics conditioned on the existence of precipitation and good observational constraints on the baseline probability of precipitation will become key for tighter constraints and better conceptual understanding.

 

Abstract Tropical cyclogenesis (TCG) is a multiscale process that involves interactions between large-scale circulation and small-scale convection. A near-global aquaplanet cloud-resolving model (NGAqua) with 4-km horizontal grid spacing that produces tropical cyclones (TCs) is used to investigate TCG and its predictability. This study analyzes an ensemble of three 20-day NGAqua simulations, with initial white-noise perturbations of low-level humidity. TCs develop spontaneously from the northern edge of the intertropical convergence zone (ITCZ), where large-scale flows and tropical convection provide necessary conditions for barotropic instability. Zonal bands of positive low-level absolute vorticity organize into cyclonic vortices, some of which develop into TCs. A new algorithm is developed to track the cyclonic vortices. A vortex-following framework analysis of the low-level vorticity budget shows that vertical stretching of absolute vorticity due to convective heating contributes positively to the vorticity spinup of the TCs. A case study and composite analyses suggest that sufficient humidity is key for convective development. TCG in these three NGAqua simulations undergoes the same series of interactions. The locations of cyclonic vortices are broadly predetermined by planetary-scale circulation and humidity patterns associated with ITCZ breakdown, which are predictable up to 10 days. Whether and when the cyclonic vortices become TCs depend on the somewhat more random feedback between convection and vorticity.