Global storm-resolving models (GSRMs) that can explicitly resolve some of deep convection are now being integrated for climate timescales. GSRMs are able to simulate more realistic precipitation distributions relative to traditional Coupled Model Intercomparison Project 6 (CMIP6) models. In this study, we present results from two-year-long integrations of a GSRM developed at Geophysical Fluid Dynamics Laboratory, eXperimental System for High-resolution prediction on Earth-to-Local Domains (X-SHiELD), for the response of precipitation to sea surface temperature warming and an isolated increase in CO2 and compare it to CMIP6 models. At leading order, X-SHiELD's response is within the range of the CMIP6 models. However, a close examination of the precipitation distribution response reveals that X-SHiELD has a different response at lower percentiles and the response of the extreme events are at the lower end of the range of CMIP6 models. A regional decomposition reveals that the difference is most pronounced for midlatitude land, where X-SHiELD shows a lower increase at intermediate percentiles and drying at lower percentiles.
The climate simulation frontier of a global storm-resolving model (GSRM; or k-scale model because of its kilometer-scale horizontal resolution) is deployed for climate change simulations. The climate sensitivity, effective radiative forcing, and relative humidity changes are assessed in multiyear atmospheric GSRM simulations with perturbed sea-surface temperatures and/or carbon dioxide concentrations. Our comparisons to conventional climate model results can build confidence in the existing climate models or highlight important areas for additional research. This GSRM’s climate sensitivity is within the range of conventional climate models, although on the lower end as the result of neutral, rather than amplifying, shortwave feedbacks. Its radiative forcing from carbon dioxide is higher than conventional climate models, and this arises from a bias in climatological clouds and an explicitly simulated high-cloud adjustment. Last, the pattern and magnitude of relative humidity changes, simulated with greater fidelity via explicitly resolving convection, are notably similar to conventional climate models.
Changes in tropical deep convection with global warming are a leading source of uncertainty for future climate projections. A comparison of the responses of active sensor measurements of cloud ice to interannual variability and next-generation global storm-resolving model (also known as k-scale models) simulations to global warming shows similar changes for events with the highest column-integrated ice. The changes reveal that the ice loading decreases outside the most active convection but increases at a rate of several percent per Kelvin surface warming in the most active convection. Disentangling thermodynamic and vertical velocity changes shows that the ice signal is strongly modulated by structural changes of the vertical wind field towards an intensification of strong convective updrafts with warming, suggesting that changes in ice loading are strongly influenced by changes in convective velocities, as well as a path toward extracting information about convective velocities from observations.
Atmospheric macroturbulence transports energy down the equator-to-pole gradient. This is represented by diffusion in energy balance models (EBMs), and EBMs have proven valuable for understanding and quantifying the pattern of surface temperature change. They typically assume climate-state-independent diffusivity, chosen to well represent the current climate, and find that this is sufficient to emulate warming response in general circulation models (GCMs). Meanwhile, model diagnoses of GCM simulations have shown that the diffusivity changes with climate. There is also ongoing development for diffusivity theories based on atmospheric dynamics. Here, we examine the role that changes in diffusivity play in the large-scale equator-to-pole contrast in surface warming in EBMs, building on previous analytic EBM theories for polar-amplified warming. New analytic theories for two formulations of climate-state-dependent diffusivity capture the results of numerical EBM solutions. For reasonable choices of parameter values, the success of the new analytic theories reveals why the change of diffusivity is limited in response to radiative forcing and does not eliminate polar-amplified warming.
An idealized aquaplanet moist global atmospheric model with realistic radiative transfer but no clouds and no convective parameterization is found to possess multiple climate equilibria. When forced symmetrically about the equator, in some cases the Inter Tropical Convergence Zone (ITCZ) migrates to an off-equatorial equilibrium position. Mechanism denial experiments prescribing relative humidity imply that radiation-circulation coupling is essential to this instability. The cross-equatorial asymmetry occurs only when the underlying slab ocean is sufficiently deep and the atmosphere's spectral dynamical core is sufficiently coarse (∼T170 or less with our control parameters). At higher resolutions, initializing with an asymmetric state indicates metastability with very slow (thousands of days) return to hemispheric symmetry. There is some sensitivity to the model timestep, which affects the time required to transition to the asymmetric state, with little effect on the equilibrium climate. The instability is enhanced when the planetary boundary layer scheme favors deeper layers or by a prescribed meridional heat transport away from the equator within the slab. The instability is not present when the model is run with a convective parameterization scheme commonly utilized in idealized moist models. We argue that the instability occurs when the asymmetric heating associated with a spontaneous ITCZ shift drives a circulation that rises poleward of the perturbed ITCZ. These results serve as a warning of the potential for instability and non-uniqueness of climate that may complicate studies with idealized models of the tropical response to perturbations in forcing.
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 use a kilometer-scale global storm resolving model (GSRM) 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. Comparisons show that the increased SST enhances the frequency of intense convection globally with large spatial and seasonal variations. Changes in the spatial pattern of intense convection are associated with changes in planetary circulation. Increases in the intense convection frequency do not necessarily reflect increases in convective available potential energy. The GSRM results are also compared with previously published traditional climate model projections.
Tropical cyclone (TC) potential intensity (PI) theory has a well-known form, consistent with a Carnot cycle interpretation of TC energetics, which relates PI to mean environmental conditions: the difference between surface and TC outflow temperatures and the air–sea enthalpy disequilibrium. PI has also been defined as a difference in convective available potential energy (CAPE) between two parcels, and quantitative assessments of future changes make use of a numerical algorithm based on this definition. Here, an analysis shows the conditions under which these Carnot and CAPE-based PI definitions are equivalent. There are multiple conditions, not previously enumerated, which in particular reveal a role for irreversible entropy production from surface evaporation. This mathematical analysis is verified by numerical calculations of PI’s sensitivity to large changes in surface-air relative humidity. To gain physical insight into the connection between the CAPE and Carnot formulations of PI, we use a recently developed analytic theory for CAPE to derive, starting from the CAPE-based definition, a new approximate formula for PI that nearly recovers the previous Carnot PI formula. The derivation shows that the difference in undilute buoyancies of saturated and environmental parcels that determines CAPE PI can in fact be expressed as a difference in the parcels’ surface moist static energy, providing a physical link between the Carnot and CAPE formulations of PI. This combination of analysis and physical interpretation builds confidence in previous numerical CAPE-based PI calculations that use climate model projections of the future tropical environment.
Purpose of Review:
Tropical cyclones (TCs) are strongly influenced by the large-scale environment of the tropics and will, therefore, be modified by climate changes. Numerical simulations designed to understand the sensitivities of TCs to environmental changes have typically followed one of two approaches: single-storm domain sizes with convection-permitting resolution and uniform thermal boundary conditions or comprehensive global high-resolution (about 50 km in the horizontal) atmospheric general circulation model (GCM) simulations. The approaches reviewed here rest between these two and are an important component of hierarchical modelling of the atmosphere: aquaplanet TC simulations.
Recent Findings:
Idealized model configurations have revealed controls on equilibrium TC size in large-domain simulations of rotating radiative-convective equilibrium. Simulations that include differential rotation (spherical geometry) but retain uniform thermal forcing have revealed a new mechanism of TC propagation change via storm-scale dynamics and show a poleward shift in genesis in response to warming. Simulations with Earth-like meridional thermal forcing gradients have isolated competing influences on TC genesis via shifts in the atmospheric general circulation and the temperature dependence of TC genesis in the absence of mean circulation changes.
Summary:
Aquaplanet simulations of TCs with variants that include or inhibit certain processes have recently emerged as a research methodology that has advanced the understanding of the climatic controls on TC activity. Looking forward, idealized boundary condition model configurations can be used as a bridge between GCM resolution and convection-permitting resolution models and as a tool for identifying additional mechanisms through which climate changes influence TC activity.
Tropical cyclone (TC)-permitting general circulation model simulations are performed with spherical geometry and uniform thermal forcing, including uniform sea surface temperature (SST) and insolation. The dependence of the TC number and TC intensity on SST is examined in a series of simulations with varied SST. The results are compared to corresponding simulations with doubly periodic f-plane geometry, rotating radiative convective equilibrium. The turbulent equilibria in simulations with spherical geometry have an inhomogenous distribution of TCs with the density of TCs increasing from low-to-high latitudes. The preferred region of TC genesis is the subtropics, but genesis shifts poleward and becomes less frequent with increasing SST. Both rotating radiative convective equilibrium and spherical geometry simulations have decreasing TC number and increasing TC intensity as SST is increased.
The sensitivity of global tropical cyclone (TC) activity to changes in a zonally-symmetric sea surface temperature (SST) distribution and the associated large-scale atmospheric circulation are investigated. High-resolution (~50-km horizontal grid spacing) atmospheric general circulation model simulations with maximum SST away from the equator are presented. Simulations with both fixed SST and slab ocean lower boundary conditions are compared.
The simulated TCs that form on the poleward flank of the Intertropical Convergence Zone (ITCZ) are tracked and changes in the frequency and intensity of those storms are analyzed between the different experiments. The total accumulated cyclone energy (ACE) increases as the location of the maximum SST shifts further away from the equator. The location of the ITCZ also shifts in conjunction with changes to the SST profile, and this plays an important role in mediating the frequency and intensity of the TCs that form within this modeling framework.
Coupled climate model simulations of volcanic eruptions and abrupt changes in CO2 concentration are compared in multiple realizations of Geophysical Fluid Dynamics Laboratory’s (GFDL) CM2.1. The change in global-mean surface temperature (GMST) is analyzed to determine whether a fast component of the climate sensitivity of relevance to the transient climate response (TCR, defined with the 1% yr−1 CO2-increase scenario) can be estimated from shorter-timescale climate changes. The fast component of the climate sensitivity estimated from the response of the climate model to volcanic forcing is similar to that of the simulations forced by abrupt CO2 changes, but is 5–15% smaller than the TCR. In addition, the partition between the top-of-atmosphere radiative restoring and ocean heat uptake is similar across radiative forcing agents. The possible asymmetry between warming and cooling climate perturbations, which may affect the utility of volcanic eruptions for estimating the TCR, is assessed by comparing simulations of abrupt CO2 doubling to abrupt CO2 halving. There is slightly less (~5%) GMST change in 0.5 × CO2 simulations than in 2 × CO2 simulations on the short (~10-yr) timescales relevant to the fast component of the volcanic signal. However, inferring the TCR from volcanic eruptions is more sensitive to uncertainties from internal climate variability and the estimation procedure.
The response of the GMST to volcanic eruptions is similar in GFDL’s CM2.1 and CM3, even though the latter has a higher TCR associated with a multidecadal timescale in its response. This is consistent with the expectation that the fast component of the climate sensitivity inferred from volcanic eruptions is a lower bound for the TCR.
Merlis, Timothy M., Tapio Schneider, S Bordoni, and I Eisenman, February 2013: Hadley Circulation Response to Orbital Precession. Part I: Aquaplanets. Journal of Climate, 26(3), DOI:10.1175/JCLI-D-11-00716.1. Abstract
The response of the monsoonal and annual-mean Hadley circulation to orbital precession is examined in an idealized atmospheric general circulation model with an aquaplanet slab-ocean lower boundary. Contrary to expectations, the simulated monsoonal Hadley circulation is weaker when perihelion occurs at the summer solstice than when aphelion occurs at the summer solstice. The angular momentum balance and energy balance are examined to understand the mechanisms that produce this result. That the summer with stronger insolation has a weaker circulation is the result of an increase in the atmosphere’s energetic stratification, the gross moist stability, which increases more than the amount required to balance the change in atmospheric energy flux divergence necessitated by the change in top-of-atmosphere net radiation. The solstice-season changes result in annual-mean Hadley circulation changes (e.g., changes in circulation strength).
Merlis, Timothy M., Tapio Schneider, S Bordoni, and I Eisenman, February 2013: Hadley Circulation Response to Orbital Precession. Part II: Subtropical Continent. Journal of Climate, 26(3), DOI:10.1175/JCLI-D-12-00149.1. Abstract
The response of the monsoonal and annual-mean Hadley circulation to orbital precession is examined in an idealized atmospheric general circulation model with a simplified representation of land surface processes in subtropical latitudes. When perihelion occurs in the summer of a hemisphere with a subtropical continent, changes in the top-of-atmosphere energy balance, together with a poleward shift of the monsoonal circulation boundary, lead to a strengthening of the monsoonal circulation. Spatial variations in surface heat capacity determine whether radiative perturbations are balanced by energy storage or by atmospheric energy fluxes. Although orbital precession does not affect annual-mean insolation, the annual-mean Hadley circulation does respond to orbital precession because its sensitivity to radiative changes varies over the course of the year: the monsoonal circulation in summer is near the angular momentum-conserving limit and responds directly to radiative changes; whereas in winter, the circulation is affected by the momentum fluxes of extratropical eddies and is less sensitive to radiative changes.
Orbital precession changes the seasonal distribution of insolation at a given latitude but not the annual mean. Hence, the correlation of paleoclimate proxies of annual-mean precipitation with orbital precession implies a nonlinear rectification in the precipitation response to seasonal solar forcing. It has previously been suggested that the relevant nonlinearity is that of the Clausius–Clapeyron relationship. Here it is argued that a different nonlinearity related to moisture advection by the atmospheric circulation is more important. When perihelion changes from one hemisphere’s summer solstice to the other’s in an idealized aquaplanet atmospheric general circulation model, annual-mean precipitation increases in the hemisphere with the brighter, warmer summer and decreases in the other hemisphere, in qualitative agreement with paleoclimate proxies that indicate such hemispherically antisymmetric climate variations. The rectification mechanism that gives rise to the precipitation changes is identified by decomposing the perturbation water vapor budget into “thermodynamic” and “dynamic” components. Thermodynamic changes (caused by changes in humidity with unchanged winds) dominate the hemispherically antisymmetric annual-mean precipitation response to precession in the absence of land–sea contrasts. The nonlinearity that enables the thermodynamic changes to affect annual-mean precipitation is a nonlinearity of moisture advection that arises because precession-induced seasonal humidity changes correlate with the seasonal cycle in low-level convergence. This interpretation is confirmed using simulations in which the Clausius–Clapeyron relationship is explicitly linearized. The thermodynamic mechanism also operates in simulations with an idealized representation of land, although in these simulations the dynamic component of the precipitation changes is also important, adding to the thermodynamic precipitation changes in some latitudes and offsetting it in others.
The response of hurricane frequency to climate changes in an aquaplanet configuration of a 50-km resolution atmospheric general circulation model is examined. The lower boundary condition is an energetically consistent slab ocean with a prescribed cross-equatorial ocean heat flux, which breaks the hemispheric symmetry and moves the Intertropical Convergence Zone (ITCZ) off the equator. In this idealized configuration, hurricane frequency increases in response to radiatively forced warming. The ITCZ shifts poleward when the model is warmed with fixed cross-equatorial ocean heat flux, and it is argued that the increase in hurricane frequency results from this poleward shift. Varying the imposed cross-equatorial ocean heat flux amplitude with fixed radiative forcing can isolate the effect of ITCZ shifts. If an increase in radiative forcing is accompanied by a reduction in the ocean heat flux amplitude such that the position of the ITCZ is unchanged, the simulated hurricane frequency decreases under warmed conditions.
