The role of aerosols and greenhouse gases in Sahel drought and recovery.

We exploit the multi-model ensemble produced by phase 5 of the Coupled Model Intercomparison Project (CMIP5) to synthesize current understanding of external forcing of Sahel rainfall change, past and future, through the lens of oceanic influence. The CMIP5 multi-model mean simulates the twentieth century evolution of Sahel rainfall, including the mid-century decline toward the driest years in the early 1980s and the partial recovery since. We exploit a physical argument linking anthropogenic emissions to the change in the temperature of the sub-tropical North Atlantic Ocean relative to the global tropical oceans to demonstrate indirect attribution of late twentieth century Sahel drought to the unique combination of aerosols and greenhouse gases that characterized the post-World War II period. The subsequent reduction in aerosol emissions around the North Atlantic that resulted from environmental legislation to curb acid rain, occurring as global tropical warming continued unabated, is consistent with the current partial recovery and with projections of future wetting. Singular Value Decomposition (SVD) applied to the above-mentioned sea surface temperature (SST) indices provides a succinct description of oceanic influence on Sahel rainfall and reveals the near-orthogonality in the influence of emissions between twentieth and twenty-first centuries: the independent effects of aerosols and greenhouse gases project on the difference of SST indices and explain past variation, while the dominance of greenhouse gases projects on their sum and explains future projection. This result challenges the assumption that because anthropogenic warming had a hand in past Sahel drought, continued warming will result in further drying. In fact, the twenty-first century dominance of greenhouse gases, unchallenged by aerosols, results in projections consistent with warming-induced strengthening of the monsoon, a response that has gained in coherence in CMIP5 compared to prior multi-model exercises.

Introduction

The Sahel witnessed an outstanding climatic shift with the abrupt onset of drought in the late 1960s. The magnitude and spatial extent of the downward trend in precipitation over the twentieth century (Greene et al. 2009) and the persistence of years of deficient rainfall through the 1970s and 1980s (Lamb 1982; Nicholson 1983; Dai et al. 2004; Ali and Lebel 2009) led scientists to hypothesize that Sahelian drought may be a sign of anthropogenic change. Some pointed to the local pressure of rapid population growth on vegetation cover and the consequent impact of an increase in albedo on the atmospheric energy budget (Charney 1975). Others linked drought to the nascent preoccupation with the impact of global anthropogenic emissions, of greenhouse gases (GHGs), and other pollutants, on the general circulation of the atmosphere (Bryson 1973).

We reconcile the just-mentioned arguments about the role of anthropogenic emissions in Sahel drought, and more recent arguments attributing global precipitation changes to aerosols, natural and anthropogenic (Robock and Liu 1994; Gillett et al. 2004; Liepert et al. 2004; Lambert et al. 2004; Zhang et al. 2007; Haywood et al. 2013) with demonstration that the dominant cause of Sahel drought lies in changes in the surface temperatures of the global oceans (Folland et al. 1986; Giannini et al. 2003). We do so by linking the unique post-World War II combination of aerosol and greenhouse gas emissions to the patterns of sea surface temperature (SST) known to affect Sahel rainfall.

On interannual to millennial time scales, Sahel rainfall variations have been ascribed to interhemispheric differences in SST, whether global, restricted to the Atlantic, or to the tropical Atlantic Ocean (Lough 1986; Shanahan et al. 2009; Park et al. 2015; Lindzen and Nigam 1987; Kang et al. 2008; Schneider et al. 2014). Recent work attempting to explain the disagreement in the sign of projected twenty-first century rainfall change (Biasutti and Giannini 2006) with differences in projected SST change added the explicit consideration of oceanic warming, resulting in a linear regression model that describes Sahel rainfall as a function of two regional SST indices: the difference of tropical North and tropical South Atlantic SST averages and the average of tropical Indo-Pacific sector SSTs (Biasutti et al. 2008). This regression model was further distilled into a single predictor, the North Atlantic Relative Index (NARI; Giannini et al. 2013), computed as the difference between SST averages over the sub-tropical North Atlantic (10°N–40°N, 75°W–15°W) and the global tropical oceans (20°S–20°N). Here, we denote these as NA and GT, respectively, so that NARI = NA − GT. In this difference, the tropical SST average (GT) captures global vertical stability (Chou and Neelin 2004): the stabilization effect of warming SSTs is communicated vertically through deep convection and is spread laterally by upper-atmosphere wave dynamics, given the weak temperature gradient constraint (Chiang and Sobel 2002). The sub-tropical North Atlantic SST average (NA) counters stabilization through moisture supply (Giannini et al. 2008; Seth et al. 2013). To the extent that Sahel rainfall variations on different time scales are the result of the underlying oceanic forcing’s preferred time scales of variation, our physical interpretation of NARI follows the paradigm that convection is in quasi-equilibrium with its environment (Emanuel et al. 1994). This assumption eliminates the need to consider time scales of variation separately, as long as these are longer than the time scales of convective adjustment (including the adjustment in the coupled ocean-atmosphere system described in Chiang and Sobel (2002)). In addition, it should be noted that a “relative index” effectively detrends variations (Vecchi et al. 2008), in the sense that insofar as both SST time series express GHG-induced warming, taking their difference largely removes it. In observations, NARI explains ~ 50% of the variance in twentieth century Sahel rainfall variability, all time scales included (Giannini et al. 2013).

In the “CMIP5 simulation of Sahel rainfall” section, we show that the CMIP5 multi-model mean reproduces the twentieth century evolution of Sahel rainfall, rendering discussion of its attribution to external forcing possible. In the “Climate change in the Sahel through the lens of oceanic influence” section, we propose an explanation that reconciles studies of variability, which typically seek to relate oceanic forcing and regional rainfall response, with studies of change, which typically seek to attribute regional response to external forcing, natural or anthropogenic, by developing an argument for indirect attribution of Sahel rainfall, through the influence of anthropogenic emissions on sea surface temperatures. In the “Transformation of predictors” sub-section, we exploit Singular Value Decomposition (SVD) to describe the linearly independent linear combinations of the two SST predictors defined above, i.e., the averages of sub-tropical North Atlantic and global tropical SSTs. In the “Linear regressions of Sahel rainfall” sub-section, we use the resulting singular vectors in bivariate linear regressions that describe their influence on Sahel rainfall. We present explicit formulas for these SVD-based predictors in the bivariate case; although, as noted by Preisendorfer (1988), such formulas had been available since a publication by Pearson (1901), their presentation in terms of non-dimensional parameters, introduced here, helps to visualize the underlying relationships and facilitates the interpretation of results. In particular, we show that when the ratio of standard deviations of the original predictors is close to one, the bivariate SVD rotation results in the sum and difference of these original predictor time series. In our case, this means that NARI, the difference of NA and GT, which we defined above based on a physical argument, is nearly proportional to one of the SVD predictors and is exhaustively complemented by the other, which is nearly proportional to the sum of the same SST indices, and represents global warming. The derivation of all necessary formulas is relegated to the Appendix. In the “Conclusions: Past is not prologue” section, we conclude with a synthesis that finds coherence in attributing past Sahel drought to anthropogenic emissions, while lending credence to projections of a wetter future.

CMIP5 simulation of Sahel rainfall

The multi-model ensemble produced by the Coupled Model Intercomparison Project in support of the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC), referred to as CMIP5 (Taylor et al. 2012), reproduces the observed twentieth century evolution of Sahel rainfall and projects wetter end of twenty-first century conditions more coherently than its CMIP3 predecessor (e.g., compare Biasutti and Giannini (2006) with Biasutti (2013)).

We restrict our analysis to 29 models participating in CMIP5, because only these have all the thermodynamic and dynamical variables necessary to evaluate their moisture budget (the subject of a parallel study). The models used are named in the tables in Supplementary Material. We focus on ensemble means, because we are interested in attribution of observed twentieth century Sahel rainfall variability, i.e., in describing its externally forced component, whether the external forcing is natural (i.e., due to variations in incoming top of the atmosphere insolation and spikes in aerosol concentration from volcanic eruptions) or anthropogenic (most notably, emissions of aerosols and GHGs from fossil fuel burning). Single-model ensemble sizes range between 1 and 10 in the historical simulations of the twentieth century and in the RCP8.5 scenario simulations of the twenty-first century. (In the case of the pre-Industrial control, typically, one simulation is run per model, albeit of varying length.) For each model and type of simulation, when more than one is available, we average realizations in the model’s ensemble mean. We then compute the multi-model mean as the average of the single models’ ensemble mean, a procedure which further suppresses the manifestation of internal variability.

The top panels in Fig. 1 display standardized Sahel rainfall (the precipitation average for land points in 10°N−20°N, 20°W−40°E, in each model’s original resolution, over July–September (Nicholson 1983; Giannini et al. 2003)): twentieth century (1900–1999) time series are depicted on the left-hand side and twenty-first century (2006–2099) time series on the right-hand side. To facilitate comparison, we standardize all the time series: each model’s ensemble mean, in the thin green lines, is standardized by its own mean and standard deviation. The multi-model mean, in the thick yellow line, is first computed as the average of the original (non-standardized) single-model ensemble mean precipitation time series and then standardized. The reason for proceeding in this order is that we wish to characterize the externally forced component of Sahel rainfall. Since the models’ ability to capture the relevant physical processes varies, if we computed the multi-model mean on the standardized single-model ensemble means, we would give all models, “good” and “bad,” equal weight. Conversely, if the multi-model mean is taken on rainfall anomalies, the internal variability that dominates models that do not respond to external forcing is weighted out, because by definition it does not have a preferred temporal sequencing. This allows the temporally coherent response to external forcing to emerge. For a comparison to non-standardized time series, i.e., Sahel rainfall anomalies in units of mm/day, the reader is referred to Fig. 2 in Biasutti (2013): there it can be seen that (i) the amplitude of variations in the multi-model mean is smaller than that in observations, as expected, given that the single real-world realization expresses the combination of internal variability and externally forced change. However, (ii) the observed evolution falls within the envelope defined by individual ensemble member simulations. The thick red line in the left panel of our Fig. 1 represents observed Sahel precipitation (Becker et al. 2013). Its correlation with the multi-model mean is 0.37 and is corroborated by a visible qualitative match: the first half of the century is characterized by variations around the mean, and is followed by a precipitous decline that lasts until the mid-1980s, and a partial recovery since. Standardization, because it normalizes the amplitude of variation, makes it possible to clearly discern the minimum in the multi-model mean, which occurs in 1982 and coincides with the first of three driest consecutive years in observations, 1982–1984, and with the eruption of El Chichón (Robock and Liu 1994). In contrast, the twenty-first century multi-model mean (top, right panel of Fig. 1) shows an upward trend, albeit with periods characterized by interannual-to-decadal variations that mask the trend, most notably at the end of the century.

Fig. 1

Time series of standardized Sahel rainfall (top) and time series of the dominant and trailing SVD modes of the predictors’ space p1 (middle) and p2 (bottom) for the single-model ensemble means of 29 CMIP5 models (thin green, turquoise, and orange lines) and for the multi-model ensemble mean (thick yellow, brown, and blue lines), from the historical simulations (1900–1999) in the left panels, and from the RCP8.5 simulations (2006–2099) in the right panels. The red lines represent observations

Fig. 2

Scatterplots of sub-tropical North Atlantic and global tropical ocean temperatures in observations (1901–1999) colored by standardized values of Sahel rainfall, with depictions of the angle of rotation ϕ, standard deviations of the original (σ1, σ2) and SVD-based predictors (λ11/2, λ21/2) and the concentration ellipse (dotted): see text for details

Climate change in the Sahel through the lens of oceanic influence

We build an argument for indirect attribution based on the influence of external forcing—external to the coupled ocean-atmosphere system—on the SST patterns that affect Sahel rainfall. In the “Transformation of predictors” sub-section, we concisely describe oceanic influence using SVD applied to our set of two original predictors for Sahel rainfall, sub-tropical North Atlantic and global tropical SSTs. In the “Linear regressions of Sahel rainfall” sub-section, we use the resulting singular vectors as predictors in bivariate regression. Results of the analyses of the multi-model means and observations are reported in Tables 1 and 2. Results of the analyses of the single-model ensemble means are collected in tables in Supplementary Material. The Appendices provide a brief review of the nomenclature and definitions for SVD in the general m-dimensional case (Singular Value Decomposition of an m-dimensional predictors’ space) and for linear regression (Use of time-dimensional singular vectors as predictors), as well as the derivation of the explicit form for the parameters of SVD in the bivariate case (Singular Value Decomposition of a two-dimensional predictors’ space; derivation of Eqs. (1)–(5)).

Table 1

Predictors’ space characteristics for the multi-model means from CMIP5 simulations and for twentieth century observations

Century	CMIP5 simulation	σ1, °CGT	σ1, °CNA	ϱ	k	ϕ, arc°	λ 1 °C2	λ 2 °C2	λ1/(λ1 + λ2)%	λ2/(λ1 + λ2)%
19th	PI control	0.03234	0.04195	0.369	0.771	62.75	0.0020	0.0008	71.90	28.10
20th	Historical	0.13150	0.12475	0.955	1.054	43.42	0.0321	0.0007	97.74	2.26
21st	RCP8.5	0.86711	0.83118	0.999	1.043	43.79	1.4422	0.0005	99.96	0.04
1901–1999	Observations	0.19565	0.22748	0.619	0.860	51.86	0.0734	0.0166	81.52	18.48

Columns from the left to right show the simulated century and the type of CMIP simulation (or observational reference period), the standard deviations of the original predictors, σ1 and σ2, in °C (for GT and NA, respectively), their correlation coefficient, ρ, the ratio of their standard deviations, k = σ1/σ2, angle ϕ from Eq. (3), in arc degree, variances λ1 and λ2, in (°C)2, respectively, explained by the leading and trailing modes, corresponding to time series p1 and p2, and the same in percent of the total variance λ1 + λ2

Table 2

Regression coefficients and skill for Sahel rainfall in CMIP5 multi-model means and in twentieth century observations

Century	CMIP5simulation	a	b	Correlation coefficient of y, ŷ
19th	PI-control	0.054	0.312	0.317
20th	Historical	−0.320	0.477	0.574
21st	RCP8.5	0.783	0.171	0.801
1901–1999	Observations	−0.182	0.607	0.633

Columns from the left to right show the simulated century and the type of CMIP simulation (or observational reference period), regression coefficients a and b of Sahel rainfall on the time series p1 and p2, respectively, and the correlation coefficient between y, Sahel rainfall simulated in CMIP5, or observed, and its predicted values from regression model (6): ŷ = ap1 + bp2

Transformation of predictors

We relate differences between our predictors’ covariance structures to differences in the external forcing applied to the pre-Industrial control and twentieth and twenty-first century simulations. The latter two types of simulations are run with time-varying external forcing, natural and anthropogenic: variations in each realization of a single model’s ensemble are a combination of internal variability and externally forced change, and their ensemble mean is taken to filter out internal variability. The former are run with constant external forcing, including CO2 concentrations held fixed at pre-Industrial levels (280 ppm), the intent being to provide a truthful estimation of each model’s internal variability. We repeat analyses on each model’s first 100 years of the pre-Industrial control simulation (typically, there is only one such simulation per model, of variable length) and on each model’s ensemble mean (with ensemble sizes varying between one and ten) for the twentieth and twenty-first century simulations.

When there are only two predictors—in our case, x1 = GT and x2 = NA—standardized and mutually uncorrelated SVD-based predictors p1, p2 are obtained from the original predictors x1, x2 by orthogonal rotation in the x1x2 plane followed by rescaling:andwhere the rotation angle ϕ and squared rescaling coefficient λ1, λ2 are given by:and

In Eqs. (3)–(4), k = σ1/σ2 is the ratio of standard deviations σ1 and σ2 of x1 and x2, respectively, and ρ is their correlation coefficient. Incidentally, k and ρ are the two non-dimensional parameters, characterizing the variance-covariance matrix of x1 and x2 up to a constant factor; to represent the same in terms of p1 and p2, one such parameter is ϕ, while the relative value of either of λ1 or λ2 where λ1 > λ2, can be another:

Derivation of Eqs. (1)–(5) is given in “Singular Value Decomposition of a two-dimensional predictors’ space; derivation of Eqs. (1)–(5)”.

Figure 2 illustrates these relationships in twentieth century observations (see the last row of Table 1 for the values of relevant parameters). The observed values of the original predictors are the coordinates of the color triangles on the (x1, x2) plane, with the corresponding standardized values of Sahel rainfall shown by color. The thick segments shown on the x1, x2 axes indicate sample standard deviations, σ1 and σ2, of these original predictors. Their values are approximately 0.20 and 0.23 °C, respectively, and their ratio is k = σ1/σ2 = 0.86. The correlation coefficient of the original predictors is ρ = 0.62. SVD-based predictors p1, p2, whose corresponding axes are shown by dash lines, are obtained by orthogonal rotation of the original predictors x1, x2 with rotation angle about ϕ = 52°, indicated by the arc connecting the x1 and p1 axes. The thick segments on the p1 and p2 axes show standard deviations (in our notation, λ11/2 = 0.27 °C and λ21/2 = 0.13 °C) of the data projections on these axes. The variance of projections is maximized for p1, among all possible directions on the x1, x2 plane, therefore λ1 > λ2, necessarily. Also by construction, the in-sample correlation coefficient between p1 and p2 is zero; since there are only two predictors, p2, being the “last” of them, minimizes the variance among all directions on the x1, x2 plane. Therefore, p2 only explains λ2/(λ1 + λ2) = 18.5% of variance in the sample (while p1 explains the other 81.5%). These features of the data and predictors’ pairs are illustrated by the shape of “concentration ellipses” (von Storch and Zwiers 1999; see their example 2.8.12, p. 43). The 95% concentration ellipse is shown by the dotted line in Fig. 2; its major semiaxes are 2.45λ11/2 and 2.45λ21/2. It should be noted that λ1 and λ2 do not necessarily determine the relative importance of the predictors p1 and p2 in modeling other variables: as color changes among triangles in Fig. 2 demonstrate, p2 is much more important than p1 in describing Sahel rainfall changes over 1901–1999. This is confirmed by the calculation of the regression coefficients in the last row of Table 2.

In models, since the ratio, k, of standard deviations of global tropical and sub-tropical North Atlantic SSTs is also close to 1 in the twentieth and twenty-first centuries, and the correlation coefficient of the SST indices, ρ, is positive (Table 1), the angle ϕ in Eq. (3) is close to 45°. That the diagonal is the direction along which is expressed the greatest covariance between x1 = GT and x2 = NA can also be clearly seen in Fig. 3. The remaining variance, orthogonal to it by construction, is larger in the twentieth (blue dots in Fig. 3) than in the twenty-first century (red squares). Since sinϕ = cosϕ for ϕ = 45o, p1 and p2, respectively, in Eqs. (1) and (2), become nearly proportional to the sum and difference of the original predictors, and specifically, p2 approximates NARI. The time series of p1 are depicted in the middle panels of Fig. 1: indeed, p represents a “global warming” mode, evident in both sub-tropical North Atlantic and global tropical ocean temperatures. p2 (~NARI) explains more variance in the twentieth than in the twenty-first century. Its standardized form is depicted in the bottom panels of Fig. 1. In the middle and bottom left-hand side panels, thick red lines represent the principal components in observations (Kaplan et al. 1998).

Fig. 3

Scatterplots of sub-tropical North Atlantic and global tropical ocean temperatures in 100 years of pre-Industrial control (green circles), in the twentieth century/historical (blue dots), and in the twenty-first century/RCP8.5 (red squares) simulations. The first 29 panels correspond to single-model ensemble means, and the last panel (in the lower right corner) corresponds to the multi-model mean

Our original predictors, GT and NA, are highly correlated in the twentieth and twenty-first century simulations, less so in the pre-Industrial control simulations (Fig. 3 and Table 1). As a consequence, the percent of variance explained by the leading mode, λ1 in Eq. (4), increases in the simulations with time-dependent external forcing, as compared to those with constant forcing. The correlation further increases from twentieth to twenty-first century and consequently so does the relative variance of the leading mode, λ1/(λ1 + λ2) = 1 − λ2/(λ1 + λ2), see Eq. (5). This behavior represents one prominent way in which external forcing expresses itself: it forces a positive correlation between the two predictors, through their common, GHG-induced warming. In addition, the relative cooling of the North Atlantic associated with aerosols contributes to the co-variability of indices with a second, preferred direction along which variance is expressed, which by construction is orthogonal to the first. This direction expresses more variance, visible in the scatter plots in Fig. 3, in the twentieth (blue dots) than in the twenty-first century (red squares).

Linear regressions of Sahel rainfall

Our SVD analysis of the original predictors confines any trend entirely to the leading mode time series, p1, i.e., the “sum” or “warming” mode, while the residual mode, p2, i.e., the “difference” or NARI, is completely devoid of any significant trend. We model standardized Sahel rainfall y, our predictand, by linear regression on p1 and p2 (see Section “Use of time-dimensional singular vectors as predictors” in the Appendix):

Regression coefficients a and b and the correlation coefficients between simulated and predicted Sahel rainfall for the multi-model mean and for observations are reported in Table 2. Since the SVD-based predictor time series are standardized and mutually uncorrelated, and predictand time series are standardized, the regression coefficients, a and b, are precisely the correlation coefficients of the predictand and corresponding predictors. In the twentieth century, correlation of the multi-model mean simulation of Sahel rainfall with its SST-based regression is 0.574 and correlation in observations is 0.633. Both terms in the regression contribute, but p2 contributes more (since b > a): late twentieth century Sahel drought is explained by the absence of North Atlantic warming, relative to warming of the global tropical oceans. The (positive) sign of b, the regression coefficient multiplying p2 (NARI), is consistent across centuries. However, its magnitude and explanatory power vis à vis Sahel rainfall are drastically diminished in the twenty-first century. Over the twenty-first century, the correlation between the simulated Sahel rainfall and its in-sample prediction is 0.801. The multi-model ensemble achieves this high a correlation essentially through a positive correlation between Sahel rainfall and warming, captured by p1. This happens because the coherence among models in projections of twenty-first century precipitation change has increased considerably from CMIP3 to CMIP5 (Biasutti 2013): a majority of models in CMIP5 project a future increase in Sahel rainfal to go along with warming (Schewe and Levermann 2017), in analogy to paleoclimate reconstructions of the “Green Sahara”, when summer insolation was stronger than now (de Menocal et al. 2000). Among the models with a twenty-first century correlation coefficient between p1 and Sahel rainfall larger than 0.45 (thus explaining ~ 20% of the variance with the single predictor p1), seven of nine project future wetting (Table S3 in Supplementary Materials). In sum, regression models between centuries are nearly orthogonal: p2, analogous to the difference of predictors or NARI, explains a larger fraction of twentieth century variation, while p1, analogous to the sum of predictors or warming, alone explains the externally forced twenty-first century change.

Conclusions: Past is not prologue

Our analysis of the multi-model mean by design emphasizes the externally forced response and leads us to propose the following indirect attribution argument for Sahelian climate change, past and future. In the second half of the twentieth century, as global dimming (Stanhill and Cohen 2001; Liepert 2002) opposed global warming in the northern hemisphere, a unique combination of anthropogenic emissions contributed to the late twentieth century drying of the Sahel through their effect on sea surface temperatures: aerosols by cooling the North Atlantic and greenhouse gases by warming the tropical oceans, especially the Indian Ocean. Warming of the global tropical oceans “upped the ante” for deep convection (Chou and Neelin 2004; Held et al. 2005), while the absence of warming in the North Atlantic reduced the moisture supply to the monsoon and thus its  potential to meet the “upped ante” and to trigger vertical instability (Giannini et al. 2008). Our argument is distinct from others previously proposed, which attributed late twentieth century Sahel drought solely to aerosols, whether through cooling of the North Atlantic or of the entire northern hemisphere (Rotstayn and Lohmann 2002; Kawase et al. 2010; Ackerley et al. 2011; Booth et al. 2012; Hwang et al. 2013; Park et al. 2015; Wang et al. 2016), in three ways. First of all, we argue for an indirect effect of anthropogenic emissions, i.e., mediated by sea surface temperatures. Linking emissions to Sahel rainfall through SSTs has the added advantage that it allows the synthesis in a single physical explanation of natural variability and anthropogenic change. Secondly, our argument is not about the role of interhemispheric gradients in SST shifting the latitudinal location of the Intertropical Convergence Zone and by extension rainfall in the Sahel. Rather, it is about quasi-equilibrium in convection as the world warms. The observed late twentieth century drying of the Sahel was more profound—longer lasting and of greater amplitude—than, e.g., the early twentieth century drought around 1910, despite the North Atlantic being cooler during the early period. This observation justifies the search for additional factors in drought, which we identify in greenhouse gas-induced tropical warming, which in the twentieth century occurred simultaneously with North Atlantic cooling. Therefore, thirdly, we argue for the combined drying effect of the two contrasting anthropogenic influences in twentieth century: drought was not caused by aerosols alone or by the cooling of the North Atlantic alone. Neither could it have been caused by greenhouse gases alone, as is made evident in simulations which include only the latter forcing (Biasutti 2013; Dong and Sutton 2015; Gaetani et al. 2017; Hill et al. 2017). Rather, drought resulted from the combination of aerosols and greenhouse gases. One influence cooled and reduced the moisture supply, while the other, warming, raised the threshold for convection—a double jeopardy. In the twenty-first century, external influence is dominated by GHG-induced warming, of both the North Atlantic and global tropical oceans. NARI, their difference, tends to zero, leaving all response to external forcing to be explained by their sum.

SVD rotation of the original predictors to their weighted sum and difference and the relative sizes of their regression coefficients (Table 2) expose the near-orthogonality of the twentieth and twenty-first centuries in the response of SSTs and rainfall to external forcing. The oceans’ translation of external forcing of predominantly anthropogenic nature leads to contrasting but not inconsistent outcomes in the past and the future. Late twentieth century drying and twenty-first century projections of wetting are consistent with different balances in the influences of aerosol and greenhouse gas emissions on regional energy and moisture budgets, underlined by one coherent physical explanation rooted in the dynamics of quasi-equilibrium in convection. The fact that emissions need to be considered holistically, not additively, explains how it is possible that while anthropogenic warming had a role in past drought, continued warming in the twenty-first century gives rise to a strengthening of the monsoon. As aerosol emissions have abated, especially around the North Atlantic, it is the combination of warmings of both North Atlantic and global tropical oceans that explains the strengthening of the monsoon: the upped ante can now be met. This situation is consistent with understanding of monsoon response to warming on paleoclimate time scales: when the ocean warms enough to contribute to an “upped ante” with increased moisture supply but not enough to shift the favored locus of convection from land to ocean (Chou et al. 2001), it reinforces the monsoon, as is the case for West Africa (Braconnot et al. 2012). With the usual caveats related to uncertainty in multi-model mean projections, that further warming may result in wetting of the Sahel is a conclusion worth reinforcing, in equal parts because it has gained in coherence between CMIP3 and CMIP5 (Biasutti 2013), and because the recovery of the rains is being experienced on the ground, in a fashion that is remarkably consistent with expectation from anthropogenic warming (West et al. 2008; Lodoun et al. 2013).

Electronic supplementary material

(PDF 98.1 kb)