Abstract
During El Niño events, increased precipitation occurs over the equatorial central eastern Pacific, corresponding to enhanced convective heating that modulates global climate by exciting atmospheric teleconnections. These precipitation anomalies are projected to shift and extend eastward in response to global warming. We show that this predicted change is caused by narrowing of the meridional span of the underlying El Niño–related sea surface temperature (SST) anomalies that leads to intensification of the meridional gradient of the SST anomalies, strengthening boundary-layer moisture convergence over the equatorial eastern Pacific, and enhancing local positive precipitation anomalies. The eastward shift and extension of these anomalies also intensify and extend eastward negative precipitation anomalies over the tropical western North Pacific, by strengthening equatorward advection of low mean moist enthalpy. Changes in El Niño–induced tropical precipitation anomalies suggest that, under global warming, El Niño events decay faster after their peak phase, thus shortening their duration.
INTRODUCTION
The El Niño–Southern Oscillation (ENSO), the strongest interannual variability signal in the climate system, is generated by air-sea interactions in the tropical Pacific (1–5). ENSO-driven convective heating anomalies over the tropical Pacific influence global climate system and thus extreme weather events, agriculture, and ecosystem around the world through atmospheric teleconnections (4, 6–9).
Global climate model (GCM) projections suggest that ENSO intensity could be highly sensitive to planetary scale future warming (10, 11). Uncertainties though, due to incomplete knowledge about forcing and feedback processes, are large (10–18). However, despite these uncertainties, projected changes in ENSO-driven precipitation anomalies are surprisingly robust, with intensification of both positive and negative anomalies over the equatorial central eastern Pacific (CEP) and tropical western Pacific, respectively (19–24). From this perspective, ENSO intensity will increase under global warming, both within the Pacific (19, 20) and beyond (25, 26). The ENSO-driven precipitation changes have been attributed to the warming of background mean sea surface temperature (SST) and associated changes in the atmospheric circulations and moisture (19, 20, 27), as well as changes in ENSO-related circulation anomalies due to pattern changes in background SST and ENSO-related SST anomalies (20–23).
In addition to the intensity changes, ENSO-driven positive precipitation anomalies over the equatorial Pacific were projected to be shifted and extended eastward (19–22) and thus lead to eastward shifts of the dominant variability mode of the Walker circulation (28) and the atmospheric teleconnection pattern over the North Pacific and North America forced by ENSO (29–33). However, the causation of this shift and extension has remained unexplained. In this study, we demonstrate that the change in ENSO-driven precipitation anomalies results from the narrowing of the meridional span of ENSO-related positive SST anomalies based on a 57-member earth system parameter-perturbed ensemble conducted with the state-of-the-art GCM HadCM3C (hereafter ESPPE; Materials and Methods).
RESULTS
During El Niño mature winters [December(0)-January-February(1), when El Niño events typically peak], the tropical Pacific is dominated by a dipole of positive and negative precipitation anomalies over the equatorial CEP and tropical western Pacific (11, 34), respectively, corresponding to the eastward shift of the Walker circulation caused by changes in the zonal gradient of SST in the equatorial Pacific. The historical simulations from the ESPPE reproduce these features (Fig. 1, A and D) well, with overall simulation skill notably higher than that of the multimodel ensemble from the phase 5 of the Coupled Model Intercomparison Project (CMIP5) (fig. S1), partly due to its high skill in background mean states (figure not shown) (35–37). Meanwhile, member spreads in simulation skill from the ESPPE and CMIP5 are comparable (fig. S1). The encouraging skill of the ESPPE suggests that it is suitable for examining projections of future ENSO-driven precipitation.
(A and B) Ensemble-mean ENSO-driven precipitation anomalies for the historical and RCP 8.5 simulations, respectively. (C) Differences between (A) and (B). A Student’s t test is conducted for the change in precipitation and SST anomalies. Values reaching 5% significance level are dotted. Yellow (green) line denotes the area where change in precipitation anomalies is less than −0.5 mm day−1 (greater than 0.5 mm day−1). (D to F) As in (A) to (C), but for ENSO-related SST anomalies.
For the climate projections under the representative concentration pathway (RCP) 8.5 scenario (hereafter RCP8.5 simulations) from the ESPPE, the ENSO-driven precipitation anomalies over both the equatorial CEP and tropical western North Pacific (WNP) are greatly intensified (Fig. 1B). The center of the positive anomalous precipitation over the equatorial CEP is displaced eastward by approximately 10° with a concomitant extension of its tail (Fig. 1B). As a result, the difference in the precipitation anomalies between the RCP8.5 and historical simulations (ΔPr′) shows a zonal dipole pattern, with a positive pole centered close to 150°W and along the equator and a negative pole centered at 155°E, 6°N (Fig. 1C). Here, the precipitation anomalies are regressed onto the original (unnormalized) Niño3.4 index time series for both the historical and RCP8.5 simulations, respectively (Materials and Methods). Thus, the obtained ENSO-driven precipitation anomalies correspond to 1°C Niño3.4 index for both simulations. A change in precipitation anomalies driven by changes in ENSO SST amplitude is excluded in consideration of the large uncertainties and controversies in the ENSO amplitude change (10–18, 38).
In contrast to the zonal dipole pattern of the ΔPr′, a change is seen in the structure of the ENSO-related SST anomalies. For ENSO events of similar Niño3.4 indices, warm SST anomalies in the equatorial CEP in the RCP8.5 simulations significantly shrink toward equator relative to that in the historical simulations (Fig. 1, D and E). Correspondingly, differences in the ENSO-related SST anomalies between the two simulations (ΔSST′) show a sandwich structure in the CEP, with a very narrow positive band along the equator, bounded by negative values on both sides (Fig. 1F).
To investigate the mechanisms responsible for the dipole pattern of the ΔPr′, we diagnose the moisture budgets for the positive and negative poles, respectively (Materials and Methods). The change in the precipitation anomalies (ΔPr′) is mainly caused by the change in vertically integrated vertical advection of climatological moisture by anomalous vertical motions (
The positive ΔPr′ over the equaotiral CEP is dominated by the
To understand the change in the three-dimensional anomalous vertical motions (Δω′), we simplify the problem by separating its horizontal distribution from the vertical structure. In the tropics, vertical motions can be approximately described as two vertical modes, deep and shallow modes (Ωd and Ωs) (Materials and Methods) (39, 40). Vertical profiles of both modes show bow structures, with the former being top-heavy and nearly having no boundary-layer convergence between 850 hPa and surface, while the latter being bottom-heavy and having a strong boundary-layer convergence (right panels of Fig. 2, B and C, and fig. S4, B and C). For both the historical and RCP8.5 simulations, the two modes account for about 77% of the total variances of vertical profiles of horizontal convergence in the tropics. On the basis of the two modes, the
(A) Vertically integrated vertical advection of climatological specific humidity by the differences in anomalous vertical motions between the RCP8.5 and the historical simulations, averaged over the ESPPE members. (B and C) Contributions of changes in deep- and shallow-mode amplitudes to that shown in (A), respectively. Red lines attached to the right of (B) and (C) are the normalized vertical profiles of the deep and shallow modes, respectively. Thick yellow and green lines denote the area where change in precipitation anomalies is less than −0.5 mm day−1 (greater than 0.5 mm day−1), as shown in Fig. 1C.
Shallow-mode vertical motions are tightly linked to boundary-layer convergence patterns, the formation of which can be understood through an analytical mixed layer model (Supplementary Materials and Methods) (41). Surface wind anomalies and associated convergence anomalies are driven by both free-tropospheric and boundary-layer processes, which can be distinguished by decomposing surface pressure anomalies into free-tropospheric (above mean height of 850 hPa) and boundary-layer components (Supplementary Materials and Methods, contours in fig. S5) (41). It is found that the change in the latter (hereafter PBL) dominates the change in the surface convergence anomalies over the equatorial CEP (fig. S5).
The PBL gradient can be simplified as a function of the SST, assuming that the boundary layer’s air temperature varies linearly from the surface to the top of the boundary layer (mean height of 850 hPa), and the 850 hPa temperature anomalies are uniform (horizontal gradients are zero). A theoretical model similar to the classical Lindzen-Nigam model (LN-like model) (42) can then be constructed for understanding the simulations from the ESPPE, in which only the simulated SST is specified as an input (Supplementary Materials and Methods). The meridional sandwich-like structure of the boundary-layer convergence can be reproduced by the LN-like model (Fig. 3, A and B). The structure corresponds to an underlying hill-shaped meridional structure of zonal-averaged SST anomalies over the CEP, which reaches maximum on the equator and has steepest gradients between 5°N and 5°S (Fig. 3C). Furthermore, the magnitudes of divergence of ΔSST′ gradient (
(A) Differences in surface convergence (shading, 10−6 s−1) anomalies between the RCP8.5 and historical simulations, averaged over the ESPPE members. (B) As in (A), but from the LN-like model driven by the simulated changes in SST anomalies only. (C) Meridional distribution of the differences in SST anomalies between two simulations zonally averaged over 155° to 110°W. Thick line shows the ensemble mean of the ESPPE. The light blue shading is the spread corresponding to 5th to 95th percentiles. (D) Scatter diagram for changes in divergence of meridional gradient of SST anomalies averaged in 4°N to 4°S, 155° to 110°W (ordinate axis) versus change in local surface convergence anomalies (abscissa axis). Their correlation coefficient is noted on the upper right corner.
The meridional width of ENSO is associated with the strength of poleward advection of oceanic temperature anomalies from the equator by mean meridional overturning circulation of the subtropical cell (STC) (43–48). The simulated STC is weakened by 30% under global warming, as a response to the robust weakening of the mean easterly trade wind over the equatorial CEP (Fig. 4) (48, 49). The weakened mean STC causes the narrowing of ENSO meridional width and further forces stronger boundary-layer convergence anomalies over the equatorial eastern Pacific, even if ENSO amplitude does not vary.
(A and B) Longitude-depth plots for ensemble-mean climatological meridional ocean current (10−2 m s−1) in D(0)JF(1) averaged over 160°E to 90°W for the historical and RCP 8.5 simulations, respectively. (C) Differences between (A) and (B). (D to F) As in (A) to (C), but for horizontal distribution of climatological 10-m wind (m s−1).
In addition to the change in the shallow-mode amplitude, the increase in the deep-mode amplitude also has a contribution to the intensification of ENSO-driven positive precipitation anomalies over the equatorial CEP (Fig. 2B). In the tropics, deep-mode vertical motion is constrained by moist static energy (MSE) budget in the atmospheric column (50). We diagnose the change in the MSE budget between the two simulations (Materials and Methods). The results indicate that the increase in the amplitude of deep-mode upward motion anomalies in the positive pole is dominated by the increase in ENSO-related net MSE flux anomalies from top of atmosphere and surface (
We now turn to focus on the formation of the negative pole of ΔPr′ over the tropical WNP (yellow contour in Fig. 1C). The negative ΔPr′ is caused by the strengthening of ENSO-related descending anomalies over the tropical WNP under global warming, while the increase in the background mean moisture almost has no contribution (fig. S3). Further decomposition of the vertical motion anomalies into deep and shallow modes indicates that the amplitude changes in the two modes have approximately equal contributions (Fig. 2 and fig. S2B). The strengthening of the shallow-mode amplitude is not caused by the change in underlying SST anomalies gradient as the case for the positive pole (Fig. 3, A and B) but mainly by free-tropospheric processes including horizontal temperature gradient on top of the boundary layer and downward momentum mixing (fig. S5). It suggests that the change in the boundary-layer convergence anomalies and associated shallow-mode amplitude is driven by intensified local anomalous descending motion (suppressed deep convections) associated with the change in the deep-mode amplitude.
The deep-mode amplitude is constrained by the MSE budget, as noted above. The suppressed convection over the tropical WNP during El Niño mature winter is primarily driven by negative horizontal advection of mean moist enthalpy by anomalous northerly anomalies (
The physical processes responsible for the change in the ENSO-driven tropical precipitation anomalies under global warming are summarized as follows: The weakened climatological STC associated with the weakened mean easterly trade wind leads to a meridional narrowing of ENSO-related warm SST anomalies by reducing poleward advection of anomalous warm water from the equator. This narrowing leads to stronger meridional temperature gradient anomalies, which enhance the ENSO-driven boundary-layer convergence anomalies over the equatorial CEP. This process plays a key role in the eastward shift and extension of ENSO-driven positive precipitation anomalies over the CEP. In addition, the changes in the mean states, including moisture content, vertical motions, and vertical stratification of MSE, and positive convection-cloud-longwave radiation feedback have contributions to the intensification of the ENSO-driven positive precipitation anomalies. The change in the positive precipitation anomalies over the CEP stimulates a cyclonic change in low-level wind anomalies to the west. The northerly component to the western flank of the cyclone intensifies ENSO-driven negative precipitation anomalies over the tropical WNP by enhancing advection of low climatological moist enthalpy from the subtropics (schematic diagram in Fig. 5).
(A) Historical simulation; (B) RCP8.5 simulation.
We also checked the results of the 31 CMIP5 models, which support the mechanism from the ESPPE, namely, that the change in the ENSO-driven precipitation anomalies over the CEP is modulated by the underlying meridional span of the ENSO-related SST anomalies. In ensemble mean of the 31 CMIP5 models, the change in ENSO-driven positive precipitation anomalies over the equatorial CEP between RCP8.5 and historical runs (ΔPr′) corresponds to the eastward shift and extension of the positive precipitation anomalies, similar to that in the ESPPE (fig. S9C and Fig. 1C). For all these models, the changes in the underlying surface convergence anomalies are highly correlated with the change in divergence of meridional gradient of the underlying SST anomalies, with their correlation coefficient reaching −0.85 (fig. S10).
The changes in the ENSO-driven precipitation anomalies have great climate impacts globally. First, the eastward shift and extension of the positive precipitation anomalies lead to the eastward shift of the Pacific–North American teleconnection pattern, thus the change in ENSO impacts on the North American winter temperature (Fig. 6A) (29–33). Second, they lead to the eastward shift of ENSO-driven anomalous Walker circulation (Fig. 6B), which is associated with the eastward shift of the dominant variability mode of the Walker circulation (28). Third, the negative ΔPr′ over the tropical WNP strengthens WNP anomalous anticyclone (WNPAC) during El Niño mature winter (Fig. 6C). The strengthened WNPAC leads to southeastern China becoming wetter in El Niño winter under global warming by transporting more moisture northward (Fig. 6). Meanwhile, the easterly anomalies to the southern flank of the WNPAC are intensified (Fig. 6C). The equatorial easterly wind anomalies stimulate eastward-propagating oceanic upwelling Kelvin waves, which contribute to the decay of El Niño (52–56). The intensified equatorial easterly anomalies tend to accelerate the decay of El Niño, which shorten its duration (Fig. 6, C and D).
DISCUSSION
(A) Ensemble mean of ENSO-related D(0)JF(1)-mean 500-hPa geopotential height anomalies (m) for the historical (contour) and RCP8.5 simulations (shading). Black dots and red crosses denote the centers of geopotential height anomalies in historical and RCP8.5 simulations, respectively. (B) As in (A), but for 200-hPa velocity potential anomalies (106 m2 s−1). (C) Simulated changes in ENSO-driven precipitation anomalies (shading; mm day−1) and 925-hPa wind anomalies (vector; m s−1) under global warming. Values reaching 5% significance level are dotted. (D) Temporal evolutions of Niño3.4 indices for El Niño events in the historical (black) and RCP8.5 (red) simulations. Thick lines denote ensemble means. Shadings denote 5th to 95th percentiles. Thin vertical lines denote the timing when Niño3.4 index transforms from positive to negative.
In summary, this study attributes the change in the ENSO-induced precipitation anomalies over the equatorial CEP under the global warming to the change in the meridional span of the ENSO-related SST anomalies based on the ESPPE. The ESPPE shows that mean STC would weaken under the global warming, corresponding to the narrowing of the ENSO-related SST anomalies, which is consistent with theoretical model results (47). In addition, we cannot rule out the possibility that the narrowing of the ENSO-related SST anomalies, in turn, is also contributed by the change in the overlying precipitation anomalies. These open questions deserve further study in the future.
MATERIALS AND METHODS
Model and experiments
This study is based on the ESPPE, a large ensemble of climate simulations produced by the HadCM3C model developed by the Met Office Hadley Centre. The HadCM3C is a state-of-the-art coupled GCM, with the regular resolutions of 2.5° × 3.75° and 1.25° × 1.25° in the atmospheric and oceanic components, respectively (57, 58). The ensemble has 57 members from different model versions, generated by simultaneously perturbing prescribed combinations of multiple parameters in the atmosphere, ocean, sulfur cycle, and terrestrial ecosystem components of HadCM3C. In the atmospheric component, 30 parameters were perturbed including, for example, the relative humidity threshold required for cloud formation, cloud droplet to rain conversion rates, and surface roughness lengths. In the ocean component, perturbed parameters included those controlling the horizontal mixing of heat and momentum and the penetration of shortwave fluxes into the ocean. Sulfur cycle parameter perturbations were designed to model uncertainties in the conversion of SO2 to sulfate aerosols and their removal. Ecosystem component parameters included perturbations to the rate of soil respiration and factors controlling the sensitivity of plant stomata to atmospheric CO2. Details about the spinning up processes and perturbed parameter combinations can be found in (59). Standard CMIP5 historical and RCP8.5 simulations were used in this study. We use the simulations from the model years 1951–2000 and 2047–2096 following historical and RCP8.5 emission pathways, respectively.
Analysis processes
Model results were analyzed through the following steps. First, model variables in December(0)-January-February(1) [D(0)JF(1)] were detrended. Second, ENSO-related anomalies were obtained by regressing detrended variables onto D(0)JF(1)-mean Niño3.4 index (area-averaged SST anomalies in 5°N to 5°S, 170° to 120°W) for the historical and RCP8.5 simulations, respectively. Third, the difference in an ENSO-related anomaly between the RCP8.5 and historical simulations represented its change under global warming (denoted by Δ). Analyses were conducted for each ensemble member, and then ensemble means were calculated. Significances of ensemble-mean differences between the two simulations were tested by Student’s t test.
Moisture and MSE equation diagnosis
Neglecting the small time-scale tendency term (∂t〈q〉′), the linearized anomalous moisture equation can be written as
Similarly, by neglecting the time tendency of anomalous MSE, the linearized anomalous MSE equation can be written as
Changes in the moisture (MSE) budget under global warming were calculated as the difference in Eq. 1 (Eq. 2) between RCP8.5 and historical simulations. For both the moisture and MSE equations, changes in advection terms under global warming can be decomposed into two linear components, which are associated with changes in circulations and state variables, respectively. For example
Deep and shallow modes of vertical motion profiles
Following (40), we obtained the deep and shallow modes of tropical vertical motions and associated horizontal convergences by following three steps. First, by applying a principal components analysis (PCA) on horizontal convergence profiles over tropical oceans (20°S to 20°N) in the historical simulations, we extracted the first two dominant modes (fig. S4A). Second, the two orthogonal modes were linearly combined to construct two new orthogonal basis modes under the condition that only one basis has projection on the near-surface (below 925 hPa) convergences (fig. S4B). The new mode having no (strong) surface convergence is referred to as a deep (shallow) mode. Third, the amplitudes of the deep and shallow modes (σd and σs) were calculated for the historical and RCP8.5 simulations, respectively. The vertical motion profiles associated with the deep and shallow modes (Ωd and Ωs) were subsequently computed by integrating the horizontal convergence downward from 100 hPa, at which vertical motion was assumed to be zero. Vertical velocity on the mean surface pressure was set to zero, and the integration residual was redistributed between the surface pressure and 100 hPa to satisfy mass continuity (fig. S4C). Hence, anomalous three-dimensional vertical motion can be approximately described by the two modes and written as
Considering that the vertical structures of the deep and shallow modes vary less under global warming (fig. S4C), from Eq. 4, the change in vertical advection of a state variable by anomalous vertical motion caused by the change in the anomalous vertical motion can be further decomposed into two subcomponents associated with amplitude changes in deep and shallow modes, respectively. For example, the first term on the right-hand side of Eq. 3 can be decomposed as
The deep-mode vertical motion was mainly constrained by the MSE budget (50). Then, the equation for the change in MSE budget can be written as
Selection of strong El Niño and La Niña events
Empirical orthogonal function analyses were applied to the monthly tropical Pacific SST anomalies (20°S to 20°N, 120°E to 80°W) for each ESPPE member. Strong El Niño events were selected on the basis of the criteria that (i) D(0)JF(1)-mean PC time series are greater than 1 SD (61) and (ii) the time series of Niño3.4 index reach the maximum between October(0) and March(1) in two consecutive years. Considering that La Niña events tend to maintain in two consecutive years, much longer than the El Nino events, the La Niña events were selected on the basis of one more criterion; that is, the PC1 index in preceding winter [D(−1)JF(0)] is larger than −0.5 SD.
SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/2/eaax4177/DC1
Supplementary Materials and Methods
Fig. S1. Comparisons of the simulation skill for the ENSO-related precipitation and SST anomalies between ESPPE and CMIP5 models by Taylor diagrams.
Fig. S2. Budget analysis of the moisture (eq. S1).
Fig. S3. Linear decomposition of change in vertically integrated vertical advection of climatological humidity by anomalous vertical motion under global warming.
Fig. S4. Derivation of deep and shallow modes of vertical motions and horizontal convergences over the tropical oceans (20°S to 20°N).
Fig. S5. Changes in ENSO-driven boundary-layer wind (vector, m s−1) and convergence anomalies (shading, 10−6 s−1) under global warming estimated by the MLM.
Fig. S6. Budget analysis of the MSE (eq. S6).
Fig. S7. Physical processes responsible for the change in the vertically integrated vertical advection of mean MSE by anomalous vertical motions associated with the change in deep-mode amplitude (
Fig. S8. Decomposition of the change in vertically integrated horizontal advection of mean moist enthalpy by ENSO-related wind anomalies (
Fig. S9. Simulated changes in the ENSO-driven D(0)JF(1)-mean precipitation and SST anomalies under global warming in 31 CMIP5 models.
Fig. S10. Scatter diagram for changes in divergence of meridional gradient of SST anomalies averaged in 4°N to 4°S, 170° to 120°W (ordinate axis) versus change in local surface convergence anomalies (abscissa axis) derived from the 31 CMIP5 models.
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.
REFERENCES AND NOTES
- Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).