1. Introduction
The intertropical convergence zone (ITCZ) is a planetary-scale band of low-level mass convergence and intense rainfall (Waliser and Gautier 1993). The ITCZ migrates north and south across the equator in response to seasonal changes in insolation. Given its dominant influence on tropical hydroclimate and importance for societies at low latitudes, the physical mechanisms controlling the ITCZ have received considerable attention (e.g., Philander et al. 1996; Broccoli et al. 2006; Kang et al. 2008; Frierson et al. 2013; Donohoe et al. 2013; Bischoff and Schneider 2014; Byrne and Schneider 2016a). The bulk of ITCZ research has focused on the mechanisms controlling its mean position or latitude, and these mechanisms are now well understood at least in the zonal-mean context [see reviews by Chiang and Friedman (2012), Schneider et al. (2014), Donohoe and Voigt (2017), and Kang et al. (2018)].
By comparison, the processes determining the width of the ITCZ are understudied despite its importance for regional and global climate. Recent work has shown that ITCZ width influences global precipitation and top-of-the-atmosphere energy balance (Su et al. 2017; Byrne and Schneider 2018) and is dynamically linked to the poleward extent of the Hadley cells and the latitudes of the storm tracks (Watt-Meyer and Frierson 2019). Furthermore, changes in ITCZ width in response to radiative forcing are strongly anticorrelated with changes in ITCZ strength (Byrne et al. 2018), suggesting that a better understanding of the processes controlling ITCZ width could help constrain ITCZ strength and tropical rainfall intensity.
We currently lack a closed theory for ITCZ width and its response to climate perturbations. In recent years, however, progress has been made in understanding how different processes affect ITCZ width. A theory derived from atmospheric mass and moist static energy budgets quantifies how the ITCZ narrows or widens in response to changes in gross moist stability, cloud radiative effects, energy/moisture advection, and moisture transport by transient eddies (Byrne and Schneider 2016a); other studies have also identified that these processes can affect ITCZ width (Bretherton and Sobel 2002; Sobel 2003; Chou and Neelin 2004; Sobel and Neelin 2006; Chou et al. 2009; Harrop and Hartmann 2016; Popp and Silvers 2017; Dixit et al. 2018; Emanuel 2019). This energetic theory has been used to investigate the processes contributing to projected narrowing of the ITCZ in coupled climate models (Byrne and Schneider 2016b) and could potentially be applied to understand observed ITCZ narrowing over recent decades (Wodzicki and Rapp 2016).
Despite these recent advances in our mechanistic understanding of ITCZ width, there are opportunities for further progress by considering the dynamics of the atmospheric boundary layer. The boundary layer momentum budget has been successfully invoked to explain the spatial distribution of tropical low-level convergence and precipitation (Lindzen and Nigam 1987; Waliser and Somerville 1994; Tomas et al. 1999; Stevens et al. 2002; Sobel 2007; Schneider and Bordoni 2008; Gonzalez et al. 2016) but has not been applied specifically to the problem of ITCZ width. There has been a long-standing debate in the atmospheric dynamics community regarding whether tropical convergence is driven by boundary layer dynamics or by thermodynamic processes higher in the atmosphere such as latent heat release [see section 2 of Sobel and Neelin (2006) for a summary of these competing arguments]. Several influential studies argue that low-level convergence—and hence ITCZ width—is primarily controlled by surface pressure gradients and boundary layer dynamics (Lindzen and Nigam 1987; Sobel and Neelin 2006; Back and Bretherton 2009). It is therefore important to develop an understanding of ITCZ width from a dynamics perspective to complement the recent studies based on energetics and thermodynamics (Chou et al. 2009; Byrne and Schneider 2016a; Harrop and Hartmann 2016; Popp and Silvers 2017; Dixit et al. 2018).
Here we use the boundary layer zonal momentum budget together with idealized simulations to investigate the dynamics of ITCZ width. In particular, we examine the momentum budget near the ITCZ edge, assess the dominant terms, and quantify the processes controlling low-level convergence and vertical velocity in that region. We also derive a range of scalings for ITCZ width and assess the importance of Ekman versus non-Ekman processes. Finally, we link ITCZ width to sea surface temperature (SST) and boundary layer moist entropy using simple extensions to the theories of Lindzen and Nigam (1987) and Emanuel (1995).
2. Simulations
a. Idealized GCM
We perform simulations using the idealized, gray-radiation GCM of O’Gorman and Schneider (2008b) based on the model developed by Frierson et al. (2006) and Frierson (2007). Insolation is set to an annual-mean profile and the longwave optical thickness is specified as a function of latitude and pressure. All the idealized simulations discussed here are run in an aquaplanet configuration with a mixed layer depth of 1 m and no horizontal energy transports (zero q fluxes). Surface fluxes of momentum, moisture and sensible heat are calculated using bulk aerodynamic formulas, with a gustiness factor of 1 m s−1 added to the surface velocities. Monin–Obukhov similarity theory is used to calculate the drag coefficients and the roughness lengths are 5 × 10−3 m for momentum and 1 × 10−5 m for moisture and sensible heat. The boundary layer height is calculated dynamically and a k-profile scheme is used to parameterize boundary layer convection. The aquaplanet configuration has no interhemispheric asymmetries in statistical steady state and consequently the ITCZ is centered on the equator. The model has a horizontal spectral resolution of T127, 20 vertical σ levels, and an integration time step of 150 s. Simulations are spun up for 700 days with averages taken over the subsequent 1000 days. See Byrne and Schneider (2016a) for an extended description of the model setup. As the model is symmetric about the equator, the ITCZ-edge latitudes discussed later have been averaged between the Northern and Southern Hemispheres.
In the set of simulations analyzed here, the longwave optical thickness is varied so as to mimic changes in greenhouse-gas concentrations. We analyze a suite of 11 simulations: The longwave optical thickness for each simulation is based on a reference profile that has been rescaled by different factors α.1
b. CMIP5 models
We also examine ITCZ width in three fixed-SST aquaplanet simulations from phase 5 of the Coupled Model Intercomparison Project (CMIP5) (Taylor et al. 2012): aquaControl, aqua4K, and aqua4xCO2. Simulations from the CNRM-CM5, FGOALS-g2, HadGEM2-A, IPSL-CM5A-LR, IPSL-CM5B-LR, MIROC5, MPI-ESM-LR, MPI-ESM-MR, and MRI-CGCM3 models are analyzed.
3. Boundary layer dynamics near the ITCZ edge
a. Definition of ITCZ width
We define ITCZ width as the distance in degrees of latitude between the northern and southern ITCZ edges. For the majority of the analyses in this paper the edges of the ITCZ are defined as the latitudes closest to the equator (north and south) at which the time-mean vertical velocity at the boundary layer top is zero. For analysis purposes we specify the boundary layer top as the σ = 0.8 level in the idealized GCM simulations, though note that the model determines the boundary layer height dynamically. This is similar to the definition of Byrne and Schneider (2016a) but in their analysis the midtropospheric mass streamfunction is used to define the ITCZ edges. We choose to focus on the boundary layer ITCZ width because it is more directly connected to the boundary layer momentum budget, which is the tool used here to investigate the ITCZ. Although we will mostly focus on the boundary layer definition of ITCZ width, below we compare how different measures of ITCZ width scale with global-mean surface temperature.
b. ITCZ width versus global-mean SST in idealized GCM simulations
The average ITCZ width (defined using boundary layer vertical velocity) across the set of simulations over which longwave optical thickness is varied is 9.6° latitude; the width varies by 4° as global-mean SST is increased from 272 to 306 K (Fig. 1a). As discussed by Byrne and Schneider (2016a), vertical velocity–based definitions of ITCZ width are considerably wider than definitions based on brightness temperature or heavy precipitation (e.g., Dias and Pauluis 2011; Wodzicki and Rapp 2016). Differences between definitions of ITCZ width based on vertical velocity versus precipitation are likely to be explained in part by convergence at the top of the boundary layer sometimes driving shallow, nonprecipitating convection rather than deep, precipitating convection—see section 4.4.1 of Sobel (2007) for a discussion of this issue.
Different measures of ITCZ width scale differently with global-mean SST (Fig. 1a). Boundary layer ITCZ width does not change monotonically with global warming, rather the ITCZ widens with warming in cooler climate and narrows with warming in hotter climates. This nonmonotonic dependence of ITCZ width on temperature in the idealized GCM is reminiscent of the behavior of other quantities including mean available potential energy (O’Gorman and Schneider 2008a), strength of the Hadley cell (O’Gorman and Schneider 2008b), tropical meridional wind (Byrne and Schneider 2016a), and the land–ocean warming contrast (Byrne and O’Gorman 2013a). In contrast, ITCZs defined in terms of the midtropospheric streamfunction and precipitation minus evaporation (P − E) generally widen with warming across the simulations (Fig. 1a). The differing magnitudes and scaling behaviors of mass-based versus moisture-based definitions of ITCZ width have been discussed previously and linked to moisture transport by transient eddies (Byrne and Schneider 2016a). We will now use boundary layer dynamics to analyze variations in ITCZ width, defined using the boundary layer vertical velocity, over a wide range of climates.
c. Zonal momentum budget
In the reference idealized GCM simulation (α = 1.0), the Coriolis and frictional terms dominate the zonal momentum budget in the vicinity of the ITCZ edge (Figs. 1b and 2a). The vertical and horizontal momentum advection terms have similar magnitudes but are substantially smaller than the Coriolis and frictional terms; the metric term is negligible (Fig. 2a) and will not be discussed further. It is common to assume Ekman balance in the tropical atmospheric boundary layer, that is, a balance between Coriolis and frictional forces (e.g., Lindzen and Nigam 1987; Emanuel 1995; Held 2001). This force balance holds approximately at the ITCZ edge across the wide range of climates simulated by the idealized GCM (Fig. 1b), suggesting that Ekman balance is a natural starting point for understanding the dynamics of ITCZ width. A recent study found that Ekman balance breaks down near the ITCZ, at least in a dry equatorial β-plane model of the tropical boundary layer, because the horizontal advection term becomes large (Gonzalez et al. 2016). In our idealized GCM simulations we note that the vertical and horizontal advection terms become relatively more important at the ITCZ edge in warmer climates as the magnitudes of the Coriolis and frictional terms decrease (Fig. 1b). However, the zeroth-order balance is between Coriolis and frictional forces across all simulations.
In the boundary layer momentum budget the frictional force provides an eastward acceleration to the zonal wind at the ITCZ edge because the tropical surface winds are easterly (Fig. 2b); there is a transfer of zonal momentum from the surface to the overlying atmosphere. This eastward acceleration is largely balanced by westward acceleration associated with equatorward advection of air with relatively low planetary angular momentum by the surface branch of the Hadley cell (the Coriolis term). For all except the coldest two simulations (Fig. 1b), both the vertical and horizontal advection terms provide a westward acceleration to the boundary layer zonal wind at the ITCZ edge by advecting low–angular momentum air. The horizontal advection term accelerates the easterly zonal wind at the ITCZ edge provided the maximum in tropical surface easterlies is poleward of the ITCZ edge, as is the case in the reference simulation (Fig. 2b). Equivalently, if the latitude of zero surface relative vorticity (where ∂usfc/∂ϕ = 0) lies poleward of the ITCZ edge, the horizontal advection term will accelerate the easterly zonal flow there. For the vertical advection term, if the latitude at which the flow is nondivergent (i.e., where ∂ω/∂p = 0) did not vary with height, the vertical advection term in (2) would be identically zero at the ITCZ edge. However, the latitude of nondivergence generally tilts poleward with height in the idealized GCM simulations, which is consistent with a midtropospheric definition of the ITCZ being wider than a definition based on vertical velocity at the top of the boundary layer (see Fig. 3 and discussion in section 4). Consequently, the vertical advection term in the zonal momentum budget is nonzero at the ITCZ edge and accelerates the easterly zonal wind for the majority of simulations (Fig. 1b).
4. Scalings for ITCZ width based on Ekman balance
We have demonstrated that the boundary layer at the ITCZ edge is approximately in Ekman balance in the idealized GCM (Fig. 1b). We now use Ekman balance to derive simple dynamical scalings for ITCZ width and apply these scalings to the idealized GCM and CMIP5 simulations.
a. Relationship between ITCZ width, wind stress, and vorticity
Testing this simple prediction in the idealized GCM, we find that the latitude of zero relative vorticity does not scale monotonically with the ITCZ edge defined in terms of boundary layer vertical velocity (black circles in Fig. 3). This implies, perhaps not surprisingly, that frictional Ekman convergence alone—where we have also neglected the component due to meridional gradients in the reciprocal of the Coriolis parameter—cannot capture the behavior of ITCZ width over the full range of idealized GCM simulations. Interestingly, there is a strong correlation between the latitude of zero relative vorticity and the ITCZ edge defined in terms of the midtropospheric streamfunction (red circles in Fig. 3) though the latitude of zero relative vorticity tends to overestimate the ITCZ width. The variation of the ITCZ edge with height and the contrasting scaling relationships between different ITCZ-width definitions (Figs. 1a and 3) are due to the momentum advection terms becoming increasingly important for driving ascent higher in the atmosphere versus at the top of the boundary layer (not shown). In summary, surface relative vorticity is closely tied to ITCZ width in the midtroposphere but the connection to ITCZ width at the top of the boundary layer is more complex.
Related to this simple estimate for the ITCZ-edge latitude, Tomas and Webster (1997) argue that the barrier between boundary layer convergence and divergence in the tropics (which is the ITCZ edge by our definition) coincides with the zero contour of the vertical component of absolute vorticity, where absolute vorticity in the zonal mean is η = f − (1/a)(∂usfc/∂ϕ). If the product of the Coriolis parameter and absolute vorticity is negative (fη < 0) this implies that the flow is inertially unstable; Tomas and Webster (1997) find evidence of such unstable regions close to the equator and argue that they are associated with cross-equatorial advection of absolute vorticity by large pressure gradients at the equator. In our idealized GCM simulations, however, there are no interhemispheric asymmetries by construction and so there is no time-mean pressure gradient at the equator and the zero contour of absolute vorticity lies at the equator. Consequently, the Tomas and Webster (1997) theory cannot explain the off-equatorial ITCZ edges in our simulations though could be relevant for solstitial conditions or other hemispherically asymmetric climate states.
A more complete description of Ekman balance that includes the GCM-simulated wind stress τx,sfc, but still neglects the meridional gradient of 1/f, suggests that the ITCZ edge should scale with the latitude where the easterly wind stress maximizes (i.e., where ∂τx,sfc/∂ϕ = 0). Estimates of the ITCZ edge based on where wind stress maximizes are very similar to those based on the latitude of zero relative vorticity (Fig. 3). This agreement indicates that, at least for the idealized GCM simulations, the simplifying assumptions of linear drag and a constant drag coefficient are reasonable when investigating ITCZ width.
Figure 3 demonstrates that the surface wind stress alone cannot fully capture the dynamics of boundary layer ITCZ width. Before considering how other terms in the momentum budget affect ITCZ width (specifically momentum advection), we first explore in more depth the implications of perfect Ekman balance for ITCZ width.
b. An Ekman scaling for ITCZ width
c. Impact of wind stress perturbations on ITCZ width: Investigations with a toy model
The diagnostic scaling (4) shows how, assuming Ekman balance, the ITCZ-edge latitude relates to the magnitude and gradient of the zonal wind stress. But how do perturbations in wind stress, for example due to changes in wind speed or surface roughness, impact ITCZ width? To develop intuition for how ITCZ width is controlled by the wind stress, consider an idealized, reference zonal wind stress profile: τref(ϕ) = υmax × (fΔp/g) sin(9ϕ) (Fig. 7a). The reference wind stress profile is specified to have a sin(9ϕ) shape so that τref = 0 at 20°N, which is roughly consistent with the observed wind stress. The wind stress is defined to be in Ekman balance with a reference boundary layer meridional wind, υ(ϕ) = υmax sin(9ϕ) (Fig. 7b). We specify υmax = −3 m s−1 and Δp = 200 hPa. For these reference profiles, the meridional wind maximizes at 10°N; the ITCZ-edge latitude is shifted marginally equatorward of this maximum (9.9°N) because of the cosϕ weighting when calculating convergence and vertical velocity (Fig. 7c). The maximum wind stress is on the poleward side of the ITCZ edge at 12.9°N (Fig. 7a) as is also the case in the idealized GCM simulations (Fig. 3). The displacement between the latitudes of maximum wind stress and ITCZ edge in this toy model is caused by meridional gradients in the reciprocal of the Coriolis parameter.
We now perturb the reference stress profile in three ways and assess how these perturbations affect the ITCZ width:
Add or subtract a constant stress: Adding a constant westerly wind stress perturbation of δτ = 0.03 N m−2 at all latitudes shifts the ITCZ edge poleward relative to its reference position by 0.8°. Adding a constant easterly wind stress of δτ = −0.03 N m−2 shifts the ITCZ edge equatorward by 1.3° (Figs. 7a–c).
Multiply or divide stress by a constant factor: Multiplying or dividing the wind stress by a constant factor μ = 2 does not shift the ITCZ edge (Figs. 7d–f), despite the imposed changes in stress being substantially larger than in the previous case where we added a constant stress.
Shift latitude of maximum stress: Manipulating the reference stress profile so as to shift the latitude of maximum stress equatorward and poleward leads to large changes in the ITCZ-edge latitude (Figs. 7g–i). Specifically, shifting the stress maximum equatorward by 2° causes an equatorward ITCZ shift of 3.4°; shifting the stress maximum poleward by 2° shifts the ITCZ edge poleward by 2.6°.
The expression (6) reveals that the ITCZ-edge shift is directly proportional to the perturbation stress δτ, and inversely proportional to the curvature of the reference τx,sfc profile. For the wind stress profiles considered in this toy model (Figs. 7a,d,g), the curvature term is positive in the vicinity of the ITCZ edge (∂2τx,sfc/∂ϕ2 > 0). Consequently, adding a westerly stress (e.g., through weakened surface easterlies) widens the ITCZ whereas adding an easterly stress narrows the ITCZ (Fig. 7c). The asymmetry in the magnitudes of ITCZ-width responses to adding a westerly versus an easterly stress perturbation (Figs. 7a–c) is not captured by the first-order Taylor expansion used to derive (6), but could be investigated with a higher-order expansion.
Perhaps intuitively, shifting the latitude of maximum wind stress toward the equator narrows the ITCZ whereas shifting it poleward widens the ITCZ (Fig. 7g). Simple analytical expressions for changes in ITCZ width due to meridional shifts in the wind stress profile are not as straightforward as for the previous two cases. Nevertheless, physical insights can be gained by examining Fig. 7g. Shifting the wind stress profile equatorward increases the stress magnitude at the reference ITCZ-edge latitude (marked by the black dot) and decreases its meridional gradient; both effects contribute to ITCZ narrowing. Similar arguments apply to shifting the wind stress profile poleward. Interestingly, the ITCZ widening is smaller in magnitude than the narrowing (2.6° versus 3.4°) despite equal prescribed shifts in the latitude of maximum wind stress. This widening versus narrowing asymmetry is caused by curvature in the wind stress profile near the ITCZ edge, further highlighting how the meridional shape of the zonal wind stress strongly influences ITCZ width.
5. Boundary layer convergence: Ekman versus non-Ekman processes
Ekman dynamics cannot fully capture the behavior of ITCZ width across all simulations (Fig. 5), implying that additional processes need to be accounted for. By defining ITCZ width as a function of vertical velocity, the influences of Ekman versus non-Ekman processes can be quantified by assessing their respective contributions to ωBL near the ITCZ edge. To perform this decomposition, we first convert the zonal momentum budget, (2), into an equation for boundary layer mass convergence and vertical velocity.
Although the Coriolis and frictional terms approximately balance one another in the zonal momentum budget (Fig. 1b), Ekman balance is insufficient for accurately estimating vertical velocity close to and within the ITCZ (Fig. 8). It is therefore clear why the latitude of the ITCZ edge cannot be fully understood through Ekman dynamics alone.
The vertical and horizontal advection terms have a relatively weak influence on the zonal momentum budget near the ITCZ edge (Fig. 2a), but it is the meridional gradients of these advection terms that determine convergence [see (8)]. The meridional gradients of the advection and frictional terms are comparable at the ITCZ edge, implying that the advection components of the vertical velocity are nonnegligible and contribute substantially to low-level convergence and ascent (Fig. 8). Interestingly, the vertical and horizontal advection terms in the zonal momentum budget at the ITCZ edge are strongly correlated with one another across the idealized GCM simulations (r = 0.99; Fig. 9). The importance of vertical advection of momentum between the boundary layer and free troposphere for determining tropical surface winds has been discussed previously (Stevens et al. 2002; Back and Bretherton 2009) so it is not surprising that this process affects ITCZ width. The horizontal advection component of vertical velocity is comparable to the vertical advection component across the idealized GCM simulations (Fig. 10). This horizontal advection term is often neglected in studies of tropical dynamics (e.g., Lindzen and Nigam 1987) but its role in shaping tropical convergence has been highlighted previously (Holton 1975; Tomas et al. 1999; Gonzalez et al. 2016). Horizontal momentum advection has also been found to be important for determining the location and spatial pattern of solstitial convergence zones (Schneider and Bordoni 2008), and its influence increases the farther the ITCZ migrates off the equator (Faulk et al. 2017). There is Ekman-driven descent at the ITCZ edge in the reference simulation (Fig. 8) and across all the idealized GCM simulations (Fig. 10). In the majority of simulations the vertical and horizontal advection terms drive convergence and ascent at the ITCZ edge (Figs. 8 and 10) and thus tend to widen the ITCZ relative to a hypothetical boundary layer in perfect Ekman balance. Figure 8 shows that vertical and horizontal momentum advection are crucial for the spatial distribution of convergence and rainfall in the tropics, at least in the idealized GCM—we will return to the influence of these advection processes on ITCZ width in the next section.
Vertical velocity is driven by components associated with meridional gradients in (i) zonal wind stress and momentum advection and (ii) the reciprocal of the Coriolis parameter [see (8)]. It is interesting to compare these components: Both are large near the ITCZ (cf. the solid and dotted black lines in Fig. 8), indicating that ITCZ width is controlled by a delicate balance between the two terms. Meridional structure in the reciprocal of the Coriolis parameter is of zeroth-order importance for ITCZ width; the strong gradient in this function close to the equator drives descent and therefore acts to narrow the ITCZ (Fig. 8). The influence of the Coriolis parameter on boundary layer convergence decays rapidly (~1/tanϕ) moving poleward away from the ITCZ.
6. Scaling for ITCZ width based on the full momentum budget
The magnitude and pattern of zonal wind stress, via its influence on frictional Ekman convergence, is a strong constraint on ITCZ width (Figs. 4–7). However, in the previous section we demonstrated that in the idealized GCM ITCZ width is controlled by convergence driven by wind stress and horizontal/vertical advection of momentum. Below we introduce an additional dynamical scaling to quantitatively investigate these competing physical controls on ITCZ width.
7. Relationships between ITCZ width and surface temperature, moist entropy
Up to this point we have studied ITCZ width in terms of dynamical quantities, specifically wind stress and momentum advection. It is also useful (and somewhat traditional in the field of tropical atmospheric dynamics; Sobel 2007) to construct theories connecting the circulation to surface thermodynamic quantities such as SST. Below we propose two such theories for ITCZ width based on boundary layer Ekman balance.
a. Lindzen–Nigam theory
Lindzen and Nigam (1987) argued that tropical surface winds and convergence are closely controlled by the SST distribution via its influence on boundary layer pressure gradients. To make the problem more tractable and link surface winds directly to SST, Lindzen and Nigam (1987) neglected horizontal and vertical momentum advection and assumed that pressure gradients vanish above the boundary layer. Despite the validity of these assumptions being challenged (e.g., Battisti et al. 1999), observations suggest that surface convergence in the tropics is largely driven by pressure gradients within (rather than above) the boundary layer (Back and Bretherton 2009). This is a strong indication that SST gradients are dynamically tied to ITCZ width.
In the idealized GCM simulations, there is reasonable agreement between the simulated ITCZ-edge latitudes and the latitudes where ∂2SST/∂ϕ2 = 0 (Fig. 11a), with a correlation coefficient of r = 0.85. The ITCZ-edge estimates from this adapted Lindzen–Nigam theory do not lie on the one-to-one line, and there are several reasons for this. Most importantly we have assumed Ekman balance, which, as discussed, cannot fully capture ITCZ width across the simulations. We have also assumed zero horizontal pressure gradients above the boundary layer and ignored vertical momentum transport, though these assumptions are not generally valid (Back and Bretherton 2009). Finally, we have neglected meridional gradients in the drag coefficient and Coriolis parameter. Nevertheless, the prediction from Lindzen–Nigam theory that the ITCZ-edge latitude should scale with the latitude where the Laplacian of SST is zero is found to approximately hold in idealized GCM simulations and provides a simple framework for understanding ITCZ-width changes in response to changing SST patterns.
b. Emanuel (1995) theory
A complementary but distinct relationship between ITCZ width and surface thermodynamics can be derived following Emanuel (1995). Emanuel considered the transition between an atmosphere in radiative–convective equilibrium (RCE) and an atmosphere with a thermally direct overturning circulation. Starting with the assumption that a RCE state is not sustainable if the product of the Coriolis parameter and absolute vorticity at the tropopause is negative (fη < 0),4 Emanuel found that a sufficient condition for the RCE state to breakdown and a circulation to emerge is for the meridional gradient in subcloud moist entropy to exceed a critical threshold. Beyond this critical threshold, thermal wind becomes inconsistent with Hide’s theorem (Hide 1969) and a circulation spins up so as to reduce the meridional moist entropy gradient. Assuming active convection and moist adiabatic lapse rates, Emanuel derived an expression for this critical moist entropy gradient and further argued that tropical surface winds can be estimated by the degree to which the actual moist entropy distribution departs from the critical gradient. Shaw and Voigt (2016) applied this theory to understand shifts in the Hadley circulation and storm tracks under climate change, but it has not been applied to ITCZ width to our knowledge.
8. Summary
The dynamical processes controlling ITCZ width have been investigated using idealized simulations together with the boundary layer zonal momentum budget and various scalings. A key result of this study is that Ekman processes place a strong constraint on ITCZ width: The boundary layer is approximately in Ekman balance (Fig. 2a) and the Ekman component of vertical velocity is large at the ITCZ edge (Fig. 8). Scalings for ITCZ width based on Ekman dynamics are useful for understanding the behavior of ITCZ width across simulations (Figs. 3–7). In particular, ITCZ width defined in terms of the midtropospheric streamfunction scales with the latitude where the surface relative vorticity is zero and the latitude where the meridional gradient in zonal wind stress is zero (Fig. 3). However, the boundary layer definition of ITCZ width does not scale with either the latitude of zero relative vorticity or with definitions of ITCZ width based on the midtropospheric streamfunction and precipitation minus evaporation. Understanding why different measures of ITCZ width respond differently to increasing temperature is a topic for future work.
Although the ITCZ edge scales with the latitude where the wind stress maximizes, the ITCZ edge robustly lies equatorward of the stress maximum due to meridional variations in the Coriolis parameter (Fig. 3). An analytical Ekman scaling for the ITCZ-edge latitude, (4), suggests that ITCZ width is directly proportional to the magnitude of the wind stress and inversely proportional to its meridional gradient. The link between wind stress magnitude and ITCZ width is also found in the idealized GCM simulations (Fig. 4). The strong control of Ekman dynamics on ITCZ width is examined in more detail using a toy model of the atmospheric boundary layer (Fig. 7). The ITCZ-width responses to simple wind stress perturbations reveal nonintuitive behavior: Adding a constant westerly (easterly) wind stress widens (narrows) the ITCZ (Figs. 7a–c). The size of the change in ITCZ width depends not only on the magnitude of the stress perturbation but also on the curvature of the reference stress profile. Multiplying a reference wind stress profile by a constant factor does not shift the ITCZ edge (Figs. 7d–f) but changes in the latitude at which the wind stress maximizes lead to large shifts in the ITCZ edge (Figs. 7g–i). Given that the angular momentum imparted by the surface to the atmosphere in the tropics is transported poleward and ultimately returned to the solid Earth by the midlatitude surface westerlies (e.g., Held 2000), it is clear that the dynamics of ITCZ width are connected to atmospheric circulations at higher latitudes (Watt-Meyer and Frierson 2019) and the problem is therefore inherently nonlocal.
Although Ekman processes are the primary influence on ITCZ width, a scaling that includes horizontal and vertical momentum advection [see (9)] is needed to quantitatively capture ITCZ width across the full range of idealized GCM simulations (Fig. 5). The magnitudes of the advection terms in the zonal momentum budget are small (Fig. 2) but it is the meridional gradients of these terms that control vertical velocity and ITCZ width, and advection contributes substantially to vertical velocity at the ITCZ edge (Fig. 8). Therefore a full, prognostic theory for ITCZ width in the idealized GCM simulations requires not only a prediction for wind stress changes but also for changes in horizontal and vertical momentum advection (though helpfully these advection terms are strongly correlated; Fig. 9). In the more comprehensive CMIP5 simulations, ITCZ width is only weakly controlled by momentum advection (Fig. 6) raising the possibility of a more dominant influence of Ekman processes on ITCZ width than suggested by the idealized GCM simulations. Additional research is required to assess the relative influences of Ekman processes versus momentum advection on ITCZ width in more realistic GCMs and in the observed atmosphere.
Because of the fast time scales of atmospheric processes relative to the slow time scales on which SST evolves, there has been a long-standing interest in developing theories for tropical circulation and precipitation given the SST distribution (e.g., Sobel 2007). Extending the Lindzen and Nigam (1987) theory, which itself is based on the assumption of Ekman balance, we find that ITCZ width scales with the latitude at which the Laplacian of SST is zero (Fig. 11). This powerful result connects ITCZ width to a well-observed and fundamental physical feature of the climate system (SST) and provides a simple dynamical framework with which to understand the responses of ITCZ width to SST variability, for example, between phases of El Niño–Southern Oscillation (Dias and Pauluis 2011) and under global warming (e.g., Byrne and Schneider 2016b). An alternative model of the tropical circulation (Emanuel 1995) suggests that ITCZ width scales with the latitude where the Laplacian of the departure of boundary layer moist entropy from a critical distribution is zero. Our extensions to the Lindzen and Nigam (1987) and Emanuel (1995) theories have some skill in capturing the behavior of ITCZ width across the idealized GCM simulations.
The analyses and scalings presented here establish a dynamical framework with which to understand the processes determining ITCZ width. This dynamical framework complements recent efforts to understand ITCZ width based on the atmospheric energy budget (Chou et al. 2009; Byrne and Schneider 2016a; Harrop and Hartmann 2016; Dixit et al. 2018) and has the potential to be more informative and predictive, given that surface convergence and the ITCZ are thought to be driven predominantly by boundary layer pressure gradients and dynamics rather than by thermodynamic effects higher in the atmosphere (Lindzen and Nigam 1987; Sobel and Neelin 2006; Back and Bretherton 2009). Ultimately, however, the tropical circulation obeys both momentum balance and energy balance. The challenge for future theoretical work on ITCZ width will be to combine dynamical and energetic perspectives into a unified theory, perhaps through the use of a quasi-equilibrium framework following Sobel and Neelin (2006). In conclusion, our investigations into the dynamics of ITCZ width have advanced fundamental understanding of the atmospheric circulation and are likely to be useful for examining ITCZ biases in more complex climate models, differences in ITCZ width over land versus ocean regions, and the past and future evolution of tropical rainfall.
Acknowledgments
This project has received funding from the EU’s Horizon 2020 Research and Innovation Programme under the Marie Skłodowska-Curie Grant Agreement 794063. We also acknowledge support from the Imperial College London Research Fellowship Scheme and useful discussions with Arnaud Czaja, David Battisti, Kerry Emanuel, and Marty Singh. We appreciate constructive comments from Spencer Hill and two anonymous reviewers that improved the manuscript. This work used JASMIN, the U.K. collaborative data analysis facility.
REFERENCES
Back, L. E., and C. S. Bretherton, 2009: On the relationship between SST gradients, boundary layer winds, and convergence over the tropical oceans. J. Climate, 22, 4182–4196, https://doi.org/10.1175/2009JCLI2392.1.
Battisti, D. S., E. S. Sarachik, and A. C. Hirst, 1999: A consistent model for the large-scale steady surface atmospheric circulation in the tropics. J. Climate, 12, 2956–2964, https://doi.org/10.1175/1520-0442(1999)012<2956:ACMFTL>2.0.CO;2.
Bischoff, T., and T. Schneider, 2014: Energetic constraints on the position of the intertropical convergence zone. J. Climate, 27, 4937–4951, https://doi.org/10.1175/JCLI-D-13-00650.1.
Bretherton, C. S., and A. H. Sobel, 2002: A simple model of a convectively coupled Walker circulation using the weak temperature gradient approximation. J. Climate, 15, 2907–2920, https://doi.org/10.1175/1520-0442(2002)015<2907:ASMOAC>2.0.CO;2.
Broccoli, A. J., K. A. Dahl, and R. J. Stouffer, 2006: Response of the ITCZ to Northern Hemisphere cooling. Geophys. Res. Lett., 33, L01702, https://doi.org/10.1029/2005GL024546.
Byrne, M. P., and P. A. O’Gorman, 2013a: Land–ocean warming contrast over a wide range of climates: Convective quasi-equilibrium theory and idealized simulations. J. Climate, 26, 4000–4016, https://doi.org/10.1175/JCLI-D-12-00262.1.
Byrne, M. P., and P. A. O’Gorman, 2013b: Link between land–ocean warming contrast and surface relative humidities in simulations with coupled climate models. Geophys. Res. Lett., 40, 5223–5227, https://doi.org/10.1002/grl.50971.
Byrne, M. P., and T. Schneider, 2016a: Energetic constraints on the width of the intertropical convergence zone. J. Climate, 29, 4709–4721, https://doi.org/10.1175/JCLI-D-15-0767.1.
Byrne, M. P., and T. Schneider, 2016b: Narrowing of the ITCZ in a warming climate: Physical mechanisms. Geophys. Res. Lett., 43, 11 350–11 357, https://doi.org/10.1002/2016GL070396.
Byrne, M. P., and T. Schneider, 2018: Atmospheric dynamics feedback: Concept, simulations, and climate implications. J. Climate, 31, 3249–3264, https://doi.org/10.1175/JCLI-D-17-0470.1.
Byrne, M. P., A. G. Pendergrass, A. D. Rapp, and K. R. Wodzicki, 2018: Response of the intertropical convergence zone to climate change: Location, width, and strength. Curr. Climate Change Rep., 4, 355–370, https://doi.org/10.1007/s40641-018-0110-5.
Chiang, J. C. H., and A. R. Friedman, 2012: Extratropical cooling, interhemispheric thermal gradients, and tropical climate change. Annu. Rev. Earth Planet. Sci., 40, 383–412, https://doi.org/10.1146/annurev-earth-042711-105545.
Chou, C., and J. D. Neelin, 2004: Mechanisms of global warming impacts on regional tropical precipitation. J. Climate, 17, 2688–2701, https://doi.org/10.1175/1520-0442(2004)017<2688:MOGWIO>2.0.CO;2.
Chou, C., J. D. Neelin, C.-A. Chen, and J.-Y. Tu, 2009: Evaluating the “rich-get-richer” mechanism in tropical precipitation change under global warming. J. Climate, 22, 1982–2005, https://doi.org/10.1175/2008JCLI2471.1.
Dias, J., and O. Pauluis, 2011: Modulations of the phase speed of convectively coupled Kelvin waves by the ITCZ. J. Atmos. Sci., 68, 1446–1459, https://doi.org/10.1175/2011JAS3630.1.
Dixit, V., O. Geoffroy, and S. C. Sherwood, 2018: Control of ITCZ width by low-level radiative heating from upper-level clouds in aquaplanet simulations. Geophys. Res. Lett., 45, 5788–5797, https://doi.org/10.1029/2018GL078292.
Donohoe, A., and A. Voigt, 2017: Why future shifts in tropical precipitation will likely be small: The location of the tropical rain belt and the hemispheric contrast of energy input to the atmosphere. Climate Extremes: Patterns and Mechanisms, Vol. 226, Amer. Geophys. Union, 115–137, https://doi.org/10.1002/9781119068020.ch8.
Donohoe, A., J. Marshall, D. Ferreira, and D. McGee, 2013: The relationship between ITCZ location and cross-equatorial atmospheric heat transport: From the seasonal cycle to the Last Glacial Maximum. J. Climate, 26, 3597–3618, https://doi.org/10.1175/JCLI-D-12-00467.1.
Emanuel, K. A., 1995: On thermally direct circulations in moist atmospheres. J. Atmos. Sci., 52, 1529–1536, https://doi.org/10.1175/1520-0469(1995)052<1529:OTDCIM>2.0.CO;2.
Emanuel, K. A., 2019: Inferences from simple models of slow, convectively coupled processes. J. Atmos. Sci., 76, 195–208, https://doi.org/10.1175/JAS-D-18-0090.1.
Faulk, S., J. Mitchell, and S. Bordoni, 2017: Effects of rotation rate and seasonal forcing on the ITCZ extent in planetary atmospheres. J. Atmos. Sci., 74, 665–678, https://doi.org/10.1175/JAS-D-16-0014.1.
Frierson, D. M. W., 2007: The dynamics of idealized convection schemes and their effect on the zonally averaged tropical circulation. J. Atmos. Sci., 64, 1959–1976, https://doi.org/10.1175/JAS3935.1.
Frierson, D. M. W., I. M. Held, and P. Zurita-Gotor, 2006: A gray-radiation aquaplanet moist GCM. Part I: Static stability and eddy scale. J. Atmos. Sci., 63, 2548–2566, https://doi.org/10.1175/JAS3753.1.
Frierson, D. M. W., and Coauthors, 2013: Contribution of ocean overturning circulation to tropical rainfall peak in the Northern Hemisphere. Nat. Geosci., 6, 940–944, https://doi.org/10.1038/ngeo1987.
Gonzalez, A. O., C. J. Slocum, R. K. Taft, and W. H. Schubert, 2016: Dynamics of the ITCZ boundary layer. J. Atmos. Sci., 73, 1577–1592, https://doi.org/10.1175/JAS-D-15-0298.1.
Harrop, B. E., and D. L. Hartmann, 2016: The role of cloud radiative heating in determining the location of the ITCZ in aquaplanet simulations. J. Climate, 29, 2741–2763, https://doi.org/10.1175/JCLI-D-15-0521.1.
Held, I. M., 2000: The general circulation of the atmosphere. Woods Hole Oceanographic Institution Fellows Project Rep. WHOI-2001-03, 54 pp.
Held, I. M., 2001: The partitioning of the poleward energy transport between the tropical ocean and atmosphere. J. Atmos. Sci., 58, 943–948, https://doi.org/10.1175/1520-0469(2001)058<0943:TPOTPE>2.0.CO;2.
Held, I. M., and A. Y. Hou, 1980: Nonlinear axially symmetric circulations in a nearly inviscid atmosphere. J. Atmos. Sci., 37, 515–533, https://doi.org/10.1175/1520-0469(1980)037<0515:NASCIA>2.0.CO;2.
Hide, R., 1969: Dynamics of the atmospheres of the major planets with an appendix on the viscous boundary layer at the rigid bounding surface of an electrically-conducting rotating fluid in the presence of a magnetic field. J. Atmos. Sci., 26, 841–853, https://doi.org/10.1175/1520-0469(1969)026<0841:DOTAOT>2.0.CO;2.
Holton, J. R., 1975: On the influence of boundary layer friction on mixed Rossby-gravity waves. Tellus, 27, 107–115, https://doi.org/10.3402/tellusa.v27i2.9893.
Kang, S. M., I. M. Held, D. M. W. Frierson, and M. Zhao, 2008: The response of the ITCZ to extratropical thermal forcing: Idealized slab-ocean experiments with a GCM. J. Climate, 21, 3521–3532, https://doi.org/10.1175/2007JCLI2146.1.
Kang, S. M., Y. Shin, and S.-P. Xie, 2018: Extratropical forcing and tropical rainfall distribution: Energetics framework and ocean Ekman advection. npj Climate Atmos. Sci., 1, 20172, https://doi.org/10.1038/s41612-017-0004-6.
Lindzen, R. S., and S. Nigam, 1987: On the role of sea surface temperature gradients in forcing low-level winds and convergence in the tropics. J. Atmos. Sci., 44, 2418–2436, https://doi.org/10.1175/1520-0469(1987)044<2418:OTROSS>2.0.CO;2.
O’Gorman, P. A., and T. Schneider, 2008a: Energy of midlatitude transient eddies in idealized simulations of changed climates. J. Climate, 21, 3815–3832, https://doi.org/10.1175/2007JCLI2065.1.
O’Gorman, P. A., and T. Schneider, 2008b: The hydrological cycle over a wide range of climates simulated with an idealized GCM. J. Climate, 21, 5797–5806, https://doi.org/10.1175/2008JCLI2099.1.
Peixóto, J. P., and A. H. Oort, 1984: Physics of climate. Rev. Mod. Phys., 56, 365, https://doi.org/10.1103/RevModPhys.56.365.
Philander, S. G. H., D. Gu, G. Lambert, T. Li, D. Halpern, N. C. Lau, and R. C. Pacanowski, 1996: Why the ITCZ is mostly north of the equator. J. Climate, 9, 2958–2972, https://doi.org/10.1175/1520-0442(1996)009<2958:WTIIMN>2.0.CO;2.
Popp, M., and L. G. Silvers, 2017: Double and single ITCZs with and without clouds. J. Climate, 30, 9147–9166, https://doi.org/10.1175/JCLI-D-17-0062.1.
Schneider, T., and S. Bordoni, 2008: Eddy-mediated regime transitions in the seasonal cycle of a Hadley circulation and implications for monsoon dynamics. J. Atmos. Sci., 65, 915–934, https://doi.org/10.1175/2007JAS2415.1.
Schneider, T., T. Bischoff, and G. H. Haug, 2014: Migrations and dynamics of the intertropical convergence zone. Nature, 513, 45–53, https://doi.org/10.1038/nature13636.
Shaw, T. A., and A. Voigt, 2016: What can moist thermodynamics tell us about circulation shifts in response to uniform warming? Geophys. Res. Lett., 43, 4566–4575, https://doi.org/10.1002/2016GL068712.
Sobel, A. H., 2003: On the coexistence of an evaporation minimum and precipitation maximum in the warm pool. J. Climate, 16, 1003–1009, https://doi.org/10.1175/1520-0442(2003)016<1003:OTCOAE>2.0.CO;2.
Sobel, A. H., 2007: Simple models of ensemble-averaged precipitation and surface wind, given the sea surface temperature. The Global Circulation of the Atmosphere, T. Schneider and A. H. Sobel, Eds., Princeton University Press, 219–251.
Sobel, A. H., and J. D. Neelin, 2006: The boundary layer contribution to intertropical convergence zones in the quasi-equilibrium tropical circulation model framework. Theor. Comput. Fluid Dyn., 20, 323–350, https://doi.org/10.1007/s00162-006-0033-y.
Stevens, B., J. Duan, J. C. McWilliams, M. Münnich, and J. D. Neelin, 2002: Entrainment, Rayleigh friction, and boundary layer winds over the tropical Pacific. J. Climate, 15, 30–44, https://doi.org/10.1175/1520-0442(2002)015<0030:ERFABL>2.0.CO;2.
Su, H., and Coauthors, 2017: Tightening of tropical ascent and high clouds key to precipitation change in a warmer climate. Nat. Commun., 8, 15771, https://doi.org/10.1038/ncomms15771.
Taylor, K. E., R. J. Stouffer, and G. A. Meehl, 2012: An overview of CMIP5 and the experiment design. Bull. Amer. Meteor. Soc., 93, 485–498, https://doi.org/10.1175/BAMS-D-11-00094.1.
Tomas, R. A., and P. J. Webster, 1997: The role of inertial instability in determining the location and strength of near-equatorial convection. Quart. J. Roy. Meteor. Soc., 123, 1445–1482, https://doi.org/10.1002/qj.49712354202.
Tomas, R. A., J. R. Holton, and P. J. Webster, 1999: The influence of cross-equatorial pressure gradients on the location of near-equatorial convection. Quart. J. Roy. Meteor. Soc., 125, 1107–1127, https://doi.org/10.1002/qj.1999.49712555603.
Waliser, D. E., and C. Gautier, 1993: A satellite-derived climatology of the ITCZ. J. Climate, 6, 2162–2174, https://doi.org/10.1175/1520-0442(1993)006<2162:ASDCOT>2.0.CO;2.
Waliser, D. E., and R. C. J. Somerville, 1994: Preferred latitudes of the intertropical convergence zone. J. Atmos. Sci., 51, 1619–1639, https://doi.org/10.1175/1520-0469(1994)051<1619:PLOTIC>2.0.CO;2.
Watt-Meyer, O., and D. M. W. Frierson, 2019: ITCZ width controls on Hadley cell extent and eddy-driven jet position, and their response to warming. J. Climate, 32, 1151–1166, https://doi.org/10.1175/JCLI-D-18-0434.1.
Wei, H.-H., and S. Bordoni, 2018: Energetic constraints on the ITCZ position in idealized simulations with a seasonal cycle. J. Adv. Model. Earth Syst., 10, 1708–1725, https://doi.org/10.1029/2018MS001313.
Wodzicki, K. R., and A. D. Rapp, 2016: Long-term characterization of the Pacific ITCZ using TRMM, GPCP, and ERA-Interim. J. Geophys. Res. Atmos., 121, 3153–3170, https://doi.org/10.1002/2015JD024458.
For the idealized GCM simulations we rescale the longwave optical thickness of a reference simulation by the following factors: α = 0.4, 0.5, 0.7, 0.8, 1.0, 1.2, 1.5, 1.8, 2.0, 2.5, 3.0.
The steady-state continuity equation written in pressure coordinates is (1/a cosϕ)∂(υ cosϕ)/∂ϕ + ∂ω/∂p = 0.
Note that Lindzen and Nigam (1987) presented their theoretical model as being applicable primarily to the zonally asymmetric tropical flow, though they did consider the zonally symmetric case in their section 5. Here, in contrast, our objective is to derive a theory linking the zonal-mean ITCZ width to the SST distribution.
We recall from section 4 that Tomas and Webster (1997), using a theory based on the dynamically equilibrated surface absolute vorticity field [and therefore distinct from Emanuel (1995)], argued that the ITCZ is bounded by the contour along which fη = 0.