# The effects of anthropogenic and volcanic aerosols and greenhouse gases on twentieth century Sahel precipitation

Rebecca Jean Herman, Alessandra Giannini, Michela Biasutti, Yochanan Kushnir

► To cite this version:

Rebecca Jean Herman, Alessandra Giannini, Michela Biasutti, Yochanan Kushnir. The effects of anthropogenic and volcanic aerosols and greenhouse gases on twentieth century Sahel precipitation. Scientific Reports, 2020, 10, 10.1038/s41598-020-68356-w. insu-03726978

HAL Id: insu-03726978

https://insu.hal.science/insu-03726978

Submitted on 28 Jul 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# The effects of anthropogenic and volcanic aerosols and greenhouse gases on twentieth century Sahel precipitation

Rebecca Jean Herman* ID<fnref n="1" />,²,³, Alessandra Giannini ID²,³, Michela Biasutti ID⁴ & Yochanan Kushnir ID⁴

There is little scientific consensus on the importance of external climate forcings—including anthropogenic aerosols, volcanic aerosols, and greenhouse gases (GHG)—relative to each other and to internal variability in dictating past and future Sahel rainfall. We address this query by relating a 3-tiered multi-model mean (MMM) over the Climate Model Intercomparison Project phase 5 “twentieth century” and pre-Industrial control simulations to observations. The comparison of single-forcing and historical simulations highlights the importance of anthropogenic and volcanic aerosols over GHG in generating forced Sahel rainfall variability in models. However, the forced MMM only accounts for a small fraction of observed variance. A residual consistency test shows that simulated internal variability cannot explain the residual observed multidecadal variability, and points to model deficiency in simulating multidecadal variability in the forced response, internal variability, or both.

The Sahel—the boundary between the North African Savanna and the Sahara Desert—experienced dramatic, long-term rainfall variability in the twentieth century which was unparalleled in the rest of the world. This variability was marked by a striking decline in rainfall between about 1960 and the early 1980s, including devastating droughts and famine in the early 1970s and 80s, which left 100,000 people dead and 750,000 dependent on food aid<fnref n="1" />. Scientific work immediately began to explore potential relationships between Sahel rainfall and a wide variety of local²,³ and global⁴,⁵ climatic factors. Giannini, et al.⁶ confirmed the importance of global over local processes by showing that an atmospheric model forced with observed global sea surface temperature (SST) alone could reproduce the profile of the first principal component of Sahel twentieth century rainfall variability, if not the amplitude, at a correlation of approximately 0.7. Studies since then have continued to focus on various global processes, reinforcing the connections between the Sahel and the temperature of ocean basins across the globe, and establishing links to internal variability—such as the El Niño-Southern Oscillation (ENSO)⁷⁻⁹ and the Atlantic Multidecadal Oscillation (AMO)⁹,¹⁰—and external forcing—such as greenhouse gases (GHG)¹¹⁻¹⁶ and volcanic and anthropogenic aerosols¹³,¹⁷,¹⁸. However, the relative importance of internal variability and different sources of external forcing remain unclear.

There is a developing consensus in the literature that anthropogenic aerosols have contributed to the Sahel drought, though there is disagreement over the prominence of this contribution and the physical mechanism that governs it¹⁰,¹¹,¹³,¹⁶⁻²⁶. The magnitude of the contribution is somewhat contentious because of disagreement about the strength of the indirect aerosol effects²⁷⁻²⁹, which may influence SSTs and global precipitation much more than the direct radiative effect³⁰⁻³², and which may cause non-linear interactions affecting both the spatial pattern (i.e. Polson et al.¹⁸ on the Asian monsoon) and even the mean³³ of the precipitation and temperature responses to other sources of forcing. The role of greenhouse gases (GHG) is even more widely debated—not just in the twentieth century¹²⁻¹⁷,²²,²⁶,³⁰, but even in the future when GHG forcing dominates¹¹,¹²,¹⁵,¹⁹. Some argue that there are also non-linear interactions between different effects of increasing GHG³⁴,³⁵ or between GHG and

<footnote n="1">¹Department of Earth and Environmental Sciences, Columbia University, New York, NY, USA. ²International Research Institute for Climate and Society, The Earth Institute at Columbia University, New York, NY, USA. ³Laboratoire de Météorologie Dynamique/IPSL, École Normale Supérieure, PSL Research University, Sorbonne Université, École Polytechnique, IP Paris, CNRS, Paris, France. ⁴Lamont-Doherty Earth Observatory of Columbia University, Palisades, NY, USA. *email: rebecca.herman@columbia.edu</footnote>
other external forcings¹⁶ and internal processes³⁶. Finally, many studies claim that SST and Sahel rainfall variation are primarily of internal origin³⁷⁻³⁹.

Many of the above studies on the Sahel focus on one or two types of forcing or on one model, or are limited to CMIP3⁴⁰, in which most models did not include indirect aerosol effects. Some, such as Giannini and Kaplan¹⁶, use a storyline approach—focusing on proposing physically-consistent pathways in order to avoid underestimating regional impacts⁴¹. Others¹⁸,²³ use fingerprinting, extracting distinct spatial and/or temporal patterns associated with different forcings and scaling the model response to match observations in order to correct sensitivity biases and avoid compensating errors in the models⁴².

We attempt to enrich the debate about the influence of external forcing and internal variability on Sahel rainfall over the twentieth century by performing an attribution study using the Coupled Model Intercomparison Project phase 5 (CMIP5)⁴³, which is the first large ensemble of coupled models to include aerosol indirect effects and run “single-forcing” model simulations, in which one external source of radiative forcing—such as greenhouse gases (GHG), anthropogenic aerosols (AA), or natural forcing (which includes volcanic aerosols and solar and orbital variations, NAT)—varies historically while the other external forcings are held at constant pre-Industrial values. We compare the evolution of spatially- and seasonally-averaged July–September (JAS) observed Sahel rainfall to that of the twentieth century single-forcing and “historical” simulations, in which all external forcings vary simultaneously (ALL). We determine the forced responses via a weighted, tiered multi-model mean (MMM) of CMIP5 simulations (see “Methods” for details), and then calculate correlation coefficients and root mean squared errors (RMSE, expressed as fraction of observed variance) to estimate the contributions of different external forcings to observed precipitation variability, using bootstrapping methods to estimate uncertainty in those statistics. To estimate noise in the MMM and significance, we use the long preindustrial control (piC) simulations, in which all external radiative forcings are held at constant preindustrial values. We employ spectral analysis of individual twentieth century and piC simulations to estimate the contribution of internal variability to observed precipitation variability at multidecadal time scales.

## Results

**Multi-model mean performance.** For each forcing experiment, we compute the MMM as follows: (1) an average across individual runs gives the ensemble mean (EM) for each CMIP5 model, (2) a weighted average across EMs gives the institution mean (IM) for each participating research institution, (3) and a weighted average across IMs gives the multi-model mean (MMM). The weights are designed to counteract attenuation of noise in ensemble and institution means that include more runs and EMs, respectively, so that they will not be underrepresented in the MMM relative to their noisier and more variable counterparts. For a formula and its derivation, see “The multi-model mean” under “Methods”.

In Fig. 1a,b, we compare Sahel twentieth century precipitation anomalies for the ALL MMM (blue line) to individual ALL runs (blue-grey lines, background) and IMs (cyan lines), and to observations from the Global Precipitation Climatology Center⁴⁴ (GPCC, black line) and the Climatic Research Unit⁴⁵ (CRU, red dotted line). Despite disagreement in the first three years, the spatial averages of the two observational records look similar enough that uncertainty in area-averaged Sahelian precipitation is considered small, and only GPCC is used throughout the rest of the paper. Despite the spread of the IMs, the standardized anomalies (Fig. 1a) reveal a striking similarity between observations and the MMM, which captures much of the observations’ multi-decadal variation by reproducing the drought of the 70s and 80s and its recovery, and even many episodes of dramatic interannual rainfall changes, most notably near 1984, the driest year in observations. Assuming the averaging was successful in preferentially filtering out internal variability present in individual model simulations, the MMM represents a consensus, forced Sahelian rainfall profile which is recognizable in the observations (Fig. 1a). However, the actual rainfall anomalies (Fig. 1b) reveal substantial attenuation of variance in the ALL MMM compared to individual simulations and to the observations.

The remaining panels of Fig. 1 display the correlations (Fig. 1c) and the RMSE (Fig. 1d) of individual IMs (cyan histogram) and of the MMM (blue dot) with observations. The blue curves show probability density functions (PDF) from bootstrapping over the IMs (see “Methods”), and represent how those statistics might change with a slightly different set of models. The correlation measures the similarity in the shape of one time series with respect to the other but is independent of relative amplitude, whereas the RMSE estimates the difference in amplitude of the simulated and observed yearly rainfall time series.

The MMM performs as well as or better than most individual IMs in both metrics, consistent with previous studies which compared other versions of multi-model means to individual models⁴⁶. Though some research institutions may appear to outperform the MMM in correlation and RMSE with twentieth century observations (notably, GISS outperforms the MMM in both), as we are comparing only one variable (precipitation) to one realization of observations in which forced and internal variability are indistinguishable, it is unclear whether these models truly capture the underlying mechanisms better than the ensemble. The RMSE values for the MMM and the IMs are near 100% of observed variance, partially reflecting the severe attenuation seen in Fig. 1b.

**Model response to different forcing experiments.** Figure 2 displays the MMMs for the three different single-forcing experiments: AA for anthropogenic aerosols (pink, Fig. 2b), NAT for natural forcing (brown, Fig. 2c), and GHG for greenhouse gases (green, Fig. 2d); and compares them to observations (black). Figure 2a again displays the ALL MMM (blue). Note that the observations correspond to the black ordinates on the left, while forced and piC model outputs (colors, including yellow) correspond to the colored ordinates on the right, which have a scale a quarter the range to facilitate comparison. The blue, pink, brown, and green shaded areas are the 95% range of bootstrapped forced MMMs. They represent agreement in the forced signal between the institutions, even though, due to small sample size, they do not fully capture the magnitude of noise in the
<figure><img src="image_3.png" /><figcaption>**Figure 1.** MMM Performance: Standardized (a) and actual (b) departures from climatology of twentieth century Sahel precipitation in individual ALL runs (blue–grey solid lines), ALL institution means (IMs, cyan), the ALL MMM (blue), and observations from GPCC (black) and CRU (red dotted line). Histogram (cyan) of correlations (c) and RMSE (d) between GPCC observations and the IMs, actual correlation (e) and RMSE (f) of the MMM with GPCC observations (blue dot), and the bootstrapping PDFs (blue curve) of the correlation (g) and RMSE (h) between the ALL MMM and GPCC observations.</figcaption></figure>

MMM caused by coincident simulated internal variability (see "Uncertainty and significance: bootstrapping and randomized bootstrapping" under "Methods"). The yellow shaded areas are also a 95% confidence interval, but they are obtained using randomly-chosen continuous subsequences of the piC runs in place of the historical simulations, where the piC simulations are taken from the same set of research institutions which provided simulations for that historical forcing experiment. The yellow shading captures the magnitude of noise in the MMM; however, it is worth note that it may be a slight overestimate since the MMM is less effective at filtering noise when there are fewer runs per model, and there is almost always only one piC run per model.

The variance of the forced MMMs over time (solid lines) and of the bootstrapped forced MMMs and randomized bootstrapped piC MMMs in a given year (shaded areas) vary from panel to panel inversely (though not proportionally) with the square root of the number of research institutions which simulated each forcing experiment (N), and are all roughly a quarter of observed variance, consistent with many precipitation fingerprinting studies, which often scale simulated precipitation up by a factor of $3-5^{18,23,42,47}$. Aside from a few exceptions, the yearly magnitudes of the forced MMMs are not significant, as they do not surpass the yellow zone consistent with noise in the MMM; this limits the detail with which we can examine the MMM directly. However, NAT (Fig. 2c) and ALL (Fig. 2a) are both significantly dry in 1982 (the year of the El Chichón eruption, near the driest observed year in 1984), and AA (Fig. 2b) and ALL both display multi-decadal variability in the second half of the century (including a partial recovery) which is characteristic of the observations and uncharacteristic of NAT and GHG (Fig. 2d).

Figure 3 displays the mean padded power spectra (PS, lines) and 95% confidence intervals (shaded areas) of the bootstrapped forced MMMs (colors other than yellow), and compares them to that of the randomized bootstrapped piC MMMs (yellow). We calculate the piC MMM using the reduced set of models that contributed the AA experiment. With only 8 contributing research institutions, the AA MMMs filter out less noise from modelled internal variability—and thus have more power at all frequencies—than the MMMs associated with the other experiments. Thus, using this reduced set of models provides a conservative estimate of the spectral noise in all four forcing experiments. Figure 3 shows that the multi-decadal variability in AA (pink) and ALL
<figure><img src="image_4.png" /><figcaption>a. ALL, N=20</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>b. AA, N=8</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>c. NAT, N=15</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>d. GHG, N=15</figcaption></figure>

**Figure 2.** Forced MMMs: Forced MMM Sahel precipitation anomalies (colored lines; right, colored ordinates) and their yearly 95% confidence intervals from bootstrapping (colored shaded areas; right, colored ordinates) over observed Sahel precipitation anomalies (black lines; left, black ordinates) and the 95% confidence interval of the piC runs from randomized bootstrapping (yellow shaded areas; right, colored ordinates). N are the number of research institutions which performed each forcing experiment. Panel (c) additionally identifies the dates of large volcanic eruptions which had different effects on the aerosol optical depth in the northern and the southern hemispheres, as well as the sign of that difference¹⁷.

(blue) is distinct from noise (yellow). It also confirms that the high-frequency variability in GHG is consistent with noise. Episodic volcanic forcing should not give rise, per se, to spectral peaks, though the observed pattern of large eruptions at the beginning and at the end of the century (see Fig. 2c) may induce some spectral power at multidecadal timescales. Since we do not detect any meaningful spectral peak in the NAT PS (brown) associated with solar variability at 11 years, we interpret the NAT MMM to be mostly the result of volcanic aerosols.

Figure 4 again displays the values (dots) and PDFs (curves) of correlation (Fig. 4a) and RMSE (Fig. 4b) between observations and the bootstrapped ALL MMMs from Fig. 1c, d (blue) and compares them to the values (dots) and PDFs (curves) for individual forcing experiments (solid curves distinguished by color) and the piC PDFs associated with the ALL experiment (dotted yellow curves). The piC PDFs corresponding to the three individual forcing experiments (which make use of only the models contributing to that experiment) are sufficiently similar to the ALL piC PDF that they are not plotted separately, with the exception of the AA piC RMSE PDF (pink dotted curve in Fig. 4b), which is wider and centered at a higher RMSE than those of the other experiments, reflecting the high variance in the yearly values seen in the yellow shaded area in Fig. 2b. Despite this difference, the *p* = 0.05 significance levels are still sufficiently similar for all four experiments for both correlation and RMSE that they are represented by a single vertical grey dashed line at the *p* = 0.05 significance level of the ALL experiment. As the NAT and GHG MMMs contain mostly high-frequency variability—which is difficult to distinguish from noise remaining in the MMM (see Fig. 3)—their PDFs are wider than the PDFs for the AA and ALL MMMs, which exhibit low-frequency variability uncharacteristic of noise in the MMM.

While the GHG MMM is not significantly better than noise at matching observations in both correlation (r = 0.06) and RMSE (100% of observed variance), ALL (r = 0.36, RMSE = 0.96), AA (r = 0.26, RMSE = 0.97), and NAT (r = 0.23, RMSE = 0.98) all achieve significance at *p* = 0.05. The discrepancy between the ALL MMM and NAT and AA individually under both metrics suggests that both anthropogenic and volcanic aerosols contribute
<figure><img src="image_8.png" /><figcaption>**Figure 3.** Forced MMM power spectra: mean (lines) and 95% confidence intervals (shaded areas) of padded power spectra (PS) of bootstrapped forced MMMs (ALL—blue, NAT—brown, AA—pink, GHG—green) and randomized bootstrapped AA piC MMMs (yellow).</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Performance of forced MMMs: probability density function (PDF) of correlations (a) and RMSE (b) of bootstrapped forced MMM twentieth century Sahel precipitation (colored curves: blue = ALL, pink = AA, brown = NAT, green = GHG) and of randomized bootstrapped piC MMM Sahel precipitation corresponding to the ALL experiment (dotted yellow curves) and the AA experiment (dotted pink curve, b) with observed twentieth century Sahel precipitation. Actual forced MMM values are represented with colored dots on the PDFs. One-sided 95% confidence level represented with grey vertical dashed lines.</figcaption></figure>

substantially to the performance of the ALL MMM. Because the metrics for both AA and NAT fall within the other's bootstrapping confidence interval, according to this analysis, AA and NAT contribute roughly equally to the performance of the observed ALL MMM.

The ALL MMM has limited predictive power as it is nearly constant, and, according to the RMSE, it leaves 96% of the variance unexplained. This unexplained variance could be due either to model deficiency or internal variability, since the MMM is designed to filter out internal variability—which will have similar characteristics but different phase across individual simulations—in favor of forced variability. Since observations include both internal and forced variability, no MMM would be able to match observations exactly. In this light, the ALL MMM correlation with observations of 0.36 is substantial. For comparison, we may liken this to simulations forced with observed SST, which reflect as best as possible observed internal climate variability as well as forced variability. As reported in Giannini et al.⁶, the correlation of the unsmoothed observations with the unsmoothed mean over version 1 of the atmospheric general circulation model developed at NASA's Goddard Space Flight Center in the framework of the Seasonal-to-Interannual Prediction Project (NSIPP1) from 1930–2000 is 0.60; the correlation of the ALL MMM with observations over the same period is not far behind at 0.47, suggesting that a large fraction of the variability that SST-forced climate models can capture is externally forced.
<figure><img src="image_10.png" /><figcaption>**Figure 5.** Residual consistency: power spectra (PS) of observed twentieth century Sahel rainfall (solid black, **a** and **c**) and the residual after removing the ALL MMM (black dotted-dashed, **b** and **d**). (**a**) and (**b**): Mean PS by model of individual ALL (**a**) and piC (**b**) runs, colored by average JAS rainfall bias of the ALL runs compared to twentieth century observations, where observed rainfall is grey, wet models are turquoise, and dry models are brown. piC PS (**b**) are additionally averaged over multiple subsections of the runs. (**c**) Tiered mean (blue dashed line) and 66% and 95% range (blue shading) of mean PS by model of individual ALL runs which were first rescaled to match twentieth century observed JAS rainfall. Also displayed are the tiered means over PS of individual forced AA, NAT, and GHG runs (colored dashed lines). The black dashed line shows the sum of the tiered mean piC PS (from panel **d**) and the ALL MMM PS (i.e. Fig. 3). (**d**): Tiered mean (orange dashed line) and 66% and 95% range (yellow shading) of mean PS by model of individual piC runs which were first rescaled so their corresponding ALL runs match twentieth century observed yearly rainfall, as in (**c**).</figcaption></figure>

**Residual consistency.** To test the role of internal variability in the CMIP5 fully coupled models, we cannot use the MMM, because internal variability will have differing phase across different simulations. Instead, we examine power at different frequencies in individual coupled runs. Figure 5a compares the padded power spectrum (PS) of twentieth century observed Sahel precipitation (solid black) to the padded PS of the ALL simulations, first estimated for individual runs, then averaged across ensemble members for each model. They are colored by the difference in the modelled and observed rainfall climatology from 1901 to 2003, where brown is used for models which are drier than observations, grey is used for models whose climatologies are near the observed climatology, and turquoise is used for models which are wetter than observations. As the individual ALL runs are single realizations, compounding forced and internal variability like observations, they are directly comparable to observations.

While there are three models (MIROC-ESM p1, MIROC-ESM-CHEM p1, and GFDL-ESM2G p1) which nearly reach the high power of the observations at a period of 100 years, these models are biased wet, and also exhibit over-estimates of high-frequency variability. Figure 5b displays the PS of the estimate of observed internal variations implied by the MMM, calculated as the residual of observations with respect to the ALL MMM (black dashed-dotted line), and compares it to its modeled counterpart, estimated as the mean PS by model of the individual piC runs, colored by the same rainfall biases used in Fig. 5a. Since there is often only one piC simulation per model, in order to reduce uncertainty in the PS, the long piC runs are divided into continuous, non-overlapping sections, and PS are taken separately for each section and then averaged together. We again see
that wet models overestimate high-frequency variability, and no model matches the low-frequency power of the residual, pointing to inconsistency between model simulations, their MMM, and observations. If the models underestimate forced variability, or if the MMM underestimates the magnitude of the modelled forced variability, this will cause the estimate of observed internal variability to be too large; so while this comparison allows us to make a statement about consistency, it does not determine whether it is simulated internal variability or our estimate of forced variability that is incorrect. However, it is clear that modelled internal variability does not contribute substantial power at low frequencies.

The PS for both the forced and piC runs are clearly stratified by modelled precipitation climatology. To investigate whether any of the models capture the observed distribution of power across different frequencies, in Fig. 5c and d we rescale the simulations by model before taking the PS and the mean by model so that the climatology of each model's ALL simulations matches observed rainfall climatology. This mostly destroys the stratification in the previous panels (see Figure S1). The distribution of model-mean scaled ALL and piC PS are represented by blue and yellow shaded areas in Fig. 5c and d, respectively. The blue and orange dashed lines in Fig. 5c and d mark the centers of these distributions with 3-tiered, unweighted means over the PS of the ALL and piC runs, respectively. The other colored dashed lines in Fig. 5c mark the tiered means over the PS of all runs in each of the three individual forcing experiments (pink = AA, brown = NAT, green = GHG) for comparison.

The black dashed line in Fig. 5c shows the sum of the tiered mean piC PS (orange dashed line from Fig. 5d) and the PS of the ALL MMM (i.e. the blue line in Fig. 3). If the MMM accurately represented the simulated forced power when scaled to the observed climatology, we would expect this sum to match the tiered mean ALL PS (blue dashed line). Instead, it falls short at low frequency, suggesting that the ratio of the variance of the ALL MMM to observed climatology underestimates the ratios of simulated forced variance to climatological Sahelian precipitation in CMIP5 models. This may be because the ensemble is biased dry, or because differing responses to forcing between models cause the consensus forced response to have lower variance than exhibited in individual models. In addition to any implications for the RMSE calculations displayed earlier, this means that the residual spectrum in Fig. 5d is an overestimate of internal variability in observations as implied by the CMIP5 ensemble.

However, it is still clear that even scaled piC simulations do not exhibit any increase in power at low frequency (Fig. 5d). Even though the inclusion of external forcing introduces low-frequency variance (Fig. 5c), the CMIP5 models are unable to capture the scale of the increase in power at low frequency in the observed PS, which exceeds the 95th percentile of rescaled ALL PS at periods longer than 65 years. Of the different forced experiments, ALL and AA are the only ones that exhibit substantial multi-decadal variability. Thus, while the variance of the ALL MMM is somewhat underestimated, the vast majority of the discrepancy in low-frequency power between simulations and observations is not due to attenuation in the MMM, but rather to model deficiency, whether in capturing the full magnitude of the forced response to AA, or in detailing the true character and magnitude of the other forced responses, low-frequency internal variability, and their interactions.

## Discussion

The analysis in this study shows that the consensus response of Sahelian precipitation to twentieth century external forcing in CMIP5 simulations, as defined by the 3-tiered multi-model mean (MMM), correlates significantly with observations. It further shows that both anthropogenic aerosols (AA) and volcanic aerosols (NAT) contribute significantly and substantially to making the CMIP5 MMM similar to observations, with AA mostly responsible for the multidecadal forced variability. Given that the performance of the ALL MMM can apparently be explained with AA and NAT alone, we conclude that GHG do not contribute to the consensus forced response of Sahel seasonal precipitation in CMIP5 models. This does not mean that GHG do not influence Sahelian precipitation in any way, or that GHG will not play a significant role in the future as the magnitude of the forcing increases. While some individual models have indicated a role for GHG in the recovery since the mid 1990s¹⁵, it is possible that the models as an ensemble do not yet capture the effects of GHG on Sahelian rainfall because the magnitude of the forcing is still too small over the historical period. Alternatively, competition between the mechanisms linking GHG forcing to Sahelian rainfall may have masked the effects of GHG by cancelling within individual simulations²⁶ or between models³⁴ in the MMM. Finally, it has been suggested that the response to GHG is inherently non-linear (e.g. different circulation responses to different magnitudes of warming in Neupane and Cook³⁶), or interacts non-linearly with other forcings (e.g. the interaction of an “upped ante” and changing moisture supply, as suggested by Giannini and Kaplan¹⁶). These non-linearities are difficult to test without the ability to compare the ALL MMM to “all but GHG” simulations, which are not widely available in CMIP5.

While the root mean squared error (RMSE) of the ALL MMM with observations is also significantly different from noise, it is 96% of the observed rainfall variance, meaning that modelled forced variability can hardly account for observed variability since the ALL MMM is hardly better than a constant prediction. Our residual consistency test showed that while the MMM is somewhat over-attenuated relative to the forced response in individual models, the discrepancy between total observed variability and total modelled variability is an order of magnitude larger, and modelled internal variability cannot account for the difference between the simulated forced response and observations.

Since modelled internal variability does not show substantial low-frequency variability while the AA MMM does, it is tempting to attribute the full magnitude of observed multi-decadal variability to AA, as many previous studies have done by focusing only on standardized trends¹¹, correlations¹⁶, or detectability in a fingerprinting framework¹⁸,²³. However, such a claim would rely heavily on assumed grid-point linearity of the climate response to different forcings, which is disputed for tropical rainfall (i.e. Giannini and Kaplan¹⁶ on GHG and anthropogenic aerosols, Polson, et al.¹⁸ on the indirect aerosol effect and spatial trend patterns in the Asian Monsoon, Lohmann and Feichter³³ on feedbacks involving the indirect aerosol effects, Neupane and Cook³⁶ on GHG-induced circulation changes over Africa, and Meehl, et al.⁴⁸ on non-linear feedbacks between solar
forcing and GHG), as well as on the accuracy of simulated forced and internal variability. In fact, it is not possible to say without further investigation into the physical pathways influencing Sahelian precipitation whether the model deficiency is in the modelled response to forcing or in modelled internal variability. Given the strong link between Sahelian rainfall and North Atlantic SST<sup>13,14,16,38,49</sup>, it is perhaps not a coincidence that models lack strong low-frequency variability both in Sahel rainfall and in internally-generated Atlantic Multidecadal Variability in SST (AMV)<sup>50</sup>. The community is currently still debating whether the observed AMV is forced by AA<sup>20,22,30</sup> or is an internal phenomenon which is linked to ocean circulation variability<sup>37,51-54</sup> and is dramatically underestimated in most models<sup>50</sup>.

Future work that focuses on characterizing and quantifying the mechanisms of influence on Sahelian precipitation in simulations and observations and using the next generation of climate models<sup>55</sup> might shed new light on whether the model/observation discrepancy documented here is due to an underestimate in the strength of the precipitation response to AA or a failure of CMIP5-class climate models to capture low-frequency internal variability.

## Methods

**Data.** Our index of Sahel rainfall variability is land-averaged precipitation anomalies for the monsoon season (July–September; JAS) over the region 12°–18°N, 20°W–40°E. For precipitation observations we use the Global Precipitation Climatology Center (GPCC) dataset<sup>44</sup>, which is quite similar to the Climate Research Unit (CRU)<sup>45</sup> dataset in average precipitation over the Sahel. The two are compared in Fig. 1, and GPCC is used for the rest of the paper. Model simulations come from the Coupled Model Intercomparison Project phase 5 (CMIP5)<sup>43</sup>, which includes simulations by over 50 models from 20 different research institutions. Not all models contribute simulations to all four historical experiments; we use all available runs (between 1 and 10 for a given model) from all models (distinct name or physics number) and research institutions that have complete data from 1901 (where the observed rainfall record begins) to 2003 (where some models stop their historical simulations). There are 14 models from 8 institutions that contributed model simulations to the AA experiment, 21 models from 15 institutions that contributed to the GHG experiment, 22 models from 15 institutions that contributed to the NAT experiment, and 51 models from 20 institutions that contributed to the ALL experiment (Table S1). Here, if the physics number is changed, it is treated as a different model under the same institution.

**The multi-model mean.** The MMM is defined as a 3-tiered, weighted average: (1) across individual simulations (runs) to get an ensemble mean (EM) for each model, (2) across EMs to get an institution mean (IM) for each research institution, and (3) across IMs to get the multi-model mean (MMM) for that experiment. While any averaging helps to filter internal variability from the MMM, the first tier focuses on reducing internal variability present in the individual runs, the second tier focuses on reducing variability between models from uncertainty in parameter values, and the third tier focuses on reducing variability between institutions from uncertainty in parameterization. A simple mean across all model simulations is very similar to the tiered mean (not shown), but tiers are used to prevent over-representation of particular parameterizations and parameter choices in the MMM and in the uncertainty and significance calculations (which will be described below under “Uncertainty and significance: bootstrapping and randomized bootstrapping”).

If a random variable (such as the internal variability component of yearly JAS Sahel precipitation) has a variance of $\sigma^2$, then the mean over $n$ realizations of that variable will have a variance of $\sigma^2/n$. The forced variability component may experience some attenuation as well due to differences in the simulated response to forcing between models. Given that the forced signal ought to be similar across simulations, we expect attenuation of internal variability to overwhelm attenuation in forced variability. Thus, means over models with more runs or over institutions with more models will have a higher signal (forced variability) to noise (internal variability) ratio than their counterparts. However, they will also have less total variability, causing them to (counterproductively) contribute less to the MMM than their more noisy counterparts. We counteract this by using weights which are inversely proportional to the expected attenuation of noise in the MMM tiers.

For a weighted mean $\sum_i w_i X_i$ between independent random variables $X_i$ with mean $\mu_i$, variance $a_i\sigma^2$, and weight $w_i$, where $\sum_i w_i = 1$, we find that:

$$ \sigma^2 \sum_i w_i X_i = E \left[ \left( \sum_i w_i X_i \right)^2 \right] - E \left[ \sum_i w_i X_i \right]^2 = \sum_i w_i^2 (E[X_i^2] - \mu_i^2) = \sigma^2 \sum_i w_i^2 a_i $$

Thus, to counteract the attenuation from a previous tier, captured in $a_i$, we define the weights as $w_i = a_i^{-5} / \sum_i a_i^{-5} \propto a_i^{-5}$. Specifically, let $f, i, m, r, N_f, N_{fim}$, and $N_{fim}$ be such that each forcing experiment $f$ is simulated by $N_f$ institutions, with $N_{fim}$ models from each institution $i$, and $N_{fim}$ runs from each model $m$, and assume that the JAS Sahel precipitation in a given year for each run $r$ has a variance of $\sigma^2$. In the first tier, where $a_{fimr} = 1$ and $w_{fimr} = \frac{1}{N_{fim}}$ (an unweighted mean), we find that the variances of the EMs are $\sigma^2_{EM_{fim}} = \sigma^2 \sum_r \frac{1}{N_{fim}^2} = \frac{\sigma^2}{N_{fim}}$, giving $a_{fim} = \frac{1}{N_{fim}}$ for the second tier. To combat this attenuation, in the second tier we define weights $w_{fim} = \frac{\sqrt{N_{fim}}}{\sum_m \sqrt{N_{fim}}} = \frac{\sqrt{N_{fim}}}{M_{fi}} \propto \sqrt{N_{fim}}$, where $M_{fi} = \sum_m \sqrt{N_{fim}}$ is the normalization constant for those weights. Using these weights, the variances of the IMs are $\sigma^2_{IM_{fi}} = \sigma^2 \sum_m \frac{N_{fim}}{M_{fi}^2} = \frac{N_f}{M_{fi}^2} \sigma^2$, giving $a_{fi} = \frac{N_f}{M_{fi}^2}$ for the third tier. Then in the third tier, $w_{fi} \propto \frac{M_{fi}}{\sqrt{N_{fi}}}$.
**Approach.** MMMs are compared to observations using correlations, which capture similarity in frequency and phase, and root mean squared errors (RMSE), which capture differences in magnitude and are expressed as a fraction of observed variance. When comparing the observations to themselves, the correlation would be 1 and the RMSE would be 0; when comparing the observations to a constant prediction, the correlation would be 0 and the RMSE would be 1 (or 100% of observed variance).

**Uncertainty and significance: bootstrapping and randomized bootstrapping.** Estimates of sampling uncertainty over all possible model parameterizations are obtained by bootstrapping (resampling with replacement) available forced IMs before calculating the MMM and corresponding correlations and RMSE, yielding probability density functions (PDF) around the MMM correlation and RMSE. This PDF can also be interpreted as a measure of agreement between CMIP5 models.

To test the null hypothesis—that all results from the forced experiments are consistent with noise in the MMM derived from modelled internal variability alone—we measure uncertainty in the MMM by repeating the bootstrapping procedure once for each of the four forced experiments, using the long, constant-forcing preindustrial control (piC) runs from the set of models contributing historical simulations to that experiment, choosing random, continuous, 103-year subsets before each bootstrap (referred to as “randomized bootstrapping”).

In addition to uncertainty derived from model parameterization, the MMM still contains noise from lingering coincident internal variability, and because bootstrapping underestimates variance when sample size is small, this procedure does not capture the full magnitude of that uncertainty (when randomizing is not used while bootstrapping the piC MMMs, for comparison, the piC confidence interval contains high-frequency variability similar to that seen around the forced MMMs in Fig. 2, not pictured). However, the length of the piC runs allows us to effectively increase the sample size of 103-year runs in the randomized bootstrapping method enough to give an accurate estimate of noise uncertainty in the MMM: this is evident from the nearly-uniform confidence intervals of the piC MMMs (yellow, Fig. 2), which contain no time-varying signal by construction.

**Residual consistency.** We evaluate consistency between modelled and observed internal and externally-forced variability by examining and comparing the power spectra (PS) of individual ALL and piC simulations. For increased sampling of frequencies we zero-pad the time series before taking the PS, and for clarity and decreased uncertainty, we average across PS from the same model before presenting the PS. This is less effective for the piC simulations, which usually contain only one (long) simulation per model. To help reduce uncertainty in the piC PS, we divide them into consecutive, non-overlapping segments of 103 years, calculate the PS of the segments separately, and average them together. To calculate the rescaled PS, we scale the individual ALL and piC runs from a given model by (mean observed twentieth century precipitation)/(mean precipitation from the ALL runs for that model) before taking the PS. We then average the PS by model and present the 66% and 95% range of the PS. We also present 3-tiered, unweighted means over all simulations of the rescaled ALL, AA, NAT, and GHG PS. The mean PS are calculated using an unweighted mean because different realizations of internal variability in the simulations do not cause attenuation of the spectral peaks characterizing that variability.

**Data availability**

Global Precipitation Climatology Center (GPCC)⁴⁴ observational data is freely available online (see https://www.esrl.noaa.gov/psd/data/gridded/data.gpcc.html) and CMIP5⁴³ model data is freely available through the Earth System Grid (see https://esgf-node.llnl.gov/projects/esgf-llnl/).

Received: 19 November 2019; Accepted: 22 June 2020

Published online: 22 July 2020

**References**

1. Bird, G. & Medina, S. *Africa Environment Outlook: Past, Present and Future Perspectives* (UNEP, Paris, 2002).

2. Charney, J. G. Dynamics of deserts and drought in the Sahel. *Q. J. R. Meteorol. Soc.* **101**, 193–202. https://doi.org/10.1002/qj.49710142802 (1975).

3. Taylor, C. M., Lambin, E. F., Stephenne, N., Harding, R. J. & Essery, R. L. H. The influence of land use change on climate in the Sahel. *J. Clim.* **15**, 3615–3629. https://doi.org/10.1175/1520-0442(2002)015<3615:TIOLUC>2.0.CO;2 (2002).

4. Folland, C. K., Palmer, T. N. & Parker, D. E. Sahel rainfall and worldwide sea temperatures, 1901–85. *Nature* **320**, 602–607. https://doi.org/10.1038/320602a0 (1986).

5. Rowell, D. P., Folland, C. K., Maskell, K. & Ward, M. N. Variability of summer rainfall over tropical North Africa (1906–92): observations and modelling. *Q. J. R. Meteorol. Soc.* **121**, 669–704. https://doi.org/10.1002/qj.49712152311 (1995).

6. Giannini, A., Saravanan, R. & Chang, P. Oceanic forcing of Sahel rainfall on interannual to interdecadal time scales. *Science* **302**, 1027–1030 (2003).

7. Pomposi, C., Giannini, A., Kushnir, Y. & Lee, D. E. Understanding Pacific Ocean influence on interannual precipitation variability in the Sahel. *Geophys. Res. Lett.* **43**, 9234–9242 (2016).

8. Parhi, P., Giannini, A., Gentine, P. & Lall, U. Resolving contrasting regional rainfall responses to El Niño over tropical Africa. *J. Clim.* **29**, 1461–1476 (2016).

9. Okonkwo, C. *et al*. Combined effect of El Niño southern oscillation and Atlantic multidecadal oscillation on Lake Chad level variability. *Cogent Geosci.* **1**, 1117829 (2015).

10. Pomposi, C., Kushnir, Y. & Giannini, A. Moisture budget analysis of SST-driven decadal Sahel precipitation variability in the twentieth century. *Clim. Dyn.* **44**, 3303–3321 (2015).

11. Held, I. M., Delworth, T. L., Lu, J., Findell, K. U. & Knutson, T. Simulation of Sahel drought in the 20th and 21st centuries. *Proc. Nat. Acad. Sci.* **102**, 17891–17896 (2005).

12. Haarsma, R. J., Selten, F. M., Weber, S. L. & Kliphuis, M. Sahel rainfall variability and response to greenhouse warming. *Geophys. Res. Lett.* **32**, L17702 (2005).
13. Ackerley, D. et al. Sensitivity of twentieth-century Sahel rainfall to sulfate aerosol and CO2 forcing. *J. Clim.* **24**, 4999–5014. https://doi.org/10.1175/JCLI-D-11-00019.1 (2011).

14. Giannini, A. et al. A unifying view of climate change in the Sahel linking intra-seasonal, interannual and longer time scales. *Environ. Res. Lett.* **8**, 024010 (2013).

15. Dong, B. & Sutton, R. Dominant role of greenhouse-gas forcing in the recovery of Sahel rainfall. *Nat. Clim. Change* **5**, 757. https://doi.org/10.1038/nclimate2664 (2015).

16. Giannini, A. & Kaplan, A. The role of aerosols and greenhouse gases in Sahel drought and recovery. *Clim. Change* **152**, 449–466 (2019).

17. Haywood, J. M., Jones, A., Bellouin, N. & Stephenson, D. Asymmetric forcing from stratospheric aerosols impacts Sahelian rainfall. *Nat. Clim. Change* **3**, 660. https://doi.org/10.1038/nclimate1857 (2013).

18. Polson, D., Bollasina, M., Hegerl, G. & Wilcox, L. Decreased monsoon precipitation in the Northern Hemisphere due to anthropogenic aerosols. *Geophys. Res. Lett.* **41**, 6023–6029 (2014).

19. Biasutti, M. & Giannini, A. Robust Sahel drying in response to late 20th century forcings. *Geophys. Res. Lett.* **33**, 706 (2006).

20. Rotstayn, L. D. & Lohmann, U. Tropical rainfall trends and the indirect aerosol effect. *J. Clim.* **15**, 2103–2116. https://doi.org/10.1175/1520-0442(2002)015<2103:TRTATI>2.0.CO;2 (2002).

21. Hwang, Y.-T., Frierson, D. M. W. & Kang, S. M. Anthropogenic sulfate aerosol and the southward shift of tropical precipitation in the late 20th century. *Geophys. Res. Lett.* **40**, 2845–2850. https://doi.org/10.1002/grl.50502 (2013).

22. Chang, C.-Y., Chiang, J., Wehner, M., Friedman, A. & Ruedy, R. Sulfate aerosol control of tropical Atlantic climate over the twentieth century. *J. Clim.* **24**, 2540–2555 (2011).

23. Undorf, S. et al. Detectable impact of local and remote anthropogenic aerosols on the 20th century changes of West African and South Asian monsoon precipitation. *J. Geophys. Res. Atmos.* **123**, 4871–4889 (2018).

24. Robock, A. & Liu, Y. The volcanic signal in Goddard Institute for Space Studies three-dimensional model simulations. *J. Clim.* **7**, 44–55 (1994).

25. Iles, C. E. & Hegerl, G. C. The global precipitation response to volcanic eruptions in the CMIP5 models. *Environ. Res. Lett.* **9**, 104012 (2014).

26. Kawase, H. et al. Physical mechanism of long-term drying trend over tropical North Africa. *Geophys. Res. Lett.* **37**, L21706 (2010).

27. McCoy, D. T. et al. The global aerosol-cloud first indirect effect estimated using MODIS, MERRA, and AeroCom. *J. Geophys. Res. Atmos.* **122**, 1779–1796. https://doi.org/10.1002/2016JD026141 (2017).

28. Penner, J. E. et al. Model intercomparison of indirect aerosol effects. *Atmos. Chem. Phys.* **6**, 3391–3405. https://doi.org/10.5194/acp-6-3391-2006 (2006).

29. Stevens, B. & Feingold, G. Untangling aerosol effects on clouds and precipitation in a buffered system. *Nature* **461**, 607. https://doi.org/10.1038/nature08281 (2009).

30. Booth, B. B., Dunstone, N. J., Halloran, P. R., Andrews, T. & Bellouin, N. Aerosols implicated as a prime driver of twentieth-century North Atlantic climate variability. *Nature* **484**, 228 (2012).

31. Wang, Y., Jiang, J. H. & Su, H. Atmospheric responses to the redistribution of anthropogenic aerosols. *J. Geophys. Res. Atmos.* **120**, 9625–9641 (2015).

32. Lin, L. et al. Changes in extreme rainfall over India and China attributed to regional aerosol-cloud interaction during the late 20th century rapid industrialization. *Geophys. Res. Lett.* **45**, 7857–7865 (2018).

33. Lohmann, U. & Feichter, J. Global indirect aerosol effects: a review. *Atmos. Chem. Phys.* **5**, 715–737. https://doi.org/10.5194/acp-5-715-2005 (2005).

34. Biasutti, M. Forced Sahel rainfall trends in the CMIP5 archive. *J. Geophys. Res. Atmos.* **118**, 1613–1623 (2013).

35. Biasutti, M., Held, I. M., Sobel, A. H. & Giannini, A. SST forcings and Sahel rainfall variability in simulations of the twentieth and twenty-first centuries. *J. Clim.* **21**, 3471–3486 (2008).

36. Neupane, N. & Cook, K. H. A nonlinear response of Sahel rainfall to Atlantic warming. *J. Clim.* **26**, 7080–7096. https://doi.org/10.1175/JCLI-D-12-00475.1 (2013).

37. Sutton, R. T. & Hodson, D. L. Atlantic Ocean forcing of North American and European summer climate. *Science* **309**, 115–118 (2005).

38. Hoerling, M., Hurrell, J., Eischeid, J. & Phillips, A. Detection and attribution of twentieth-century northern and southern African rainfall change. *J. Clim.* **19**, 3989–4008 (2006).

39. Ting, M., Kushnir, Y., Seager, R. & Li, C. Forced and internal twentieth-century SST trends in the North Atlantic. *J. Clim.* **22**, 1469–1481 (2009).

40. Meehl, G. et al. The global coupled model dataset: a new era in climate change research. *Bull. Am. Meteorol. Soc.* **88**, 1383–1394 (2007).

41. Shepherd, T. G. Storyline approach to the construction of regional climate change information. *Proc. R. Soc. A* **475**, 20190013 (2019).

42. Hegerl, G. & Zwiers, F. Use of models in detection and attribution of climate change. *Wiley Interdiscip. Rev. Clim. Change* **2**, 570–591 (2011).

43. Taylor, K. E., Stouffer, R. J. & Meehl, G. A. An overview of CMIP5 and the experiment design. *Bull. Am. Meteorol. Soc.* **93**, 485–498 (2012).

44. Becker, A. et al. A description of the global land-surface precipitation data products of the Global Precipitation Climatology Centre with sample applications including centennial (trend) analysis from 1901-present. *Earth Syst. Sci. Data* **5**, 71–99. https://doi.org/10.5194/essd-5-71-2013 (2013).

45. Harris, I., Jones, P. D., Osborn, T. J. & Lister, D. H. Updated high-resolution grids of monthly climatic observations—the CRU TS3.10 Dataset. *Int. J. Climatol.* **34**, 623–642 (2014).

46. Gillett, N. et al. Detecting anthropogenic influence with a multi-model ensemble. *Geophys. Res. Lett.* **29**, 31–34 (2002).

47. Zhang, X. et al. Detection of human influence on twentieth-century precipitation trends. *Nature* **448**, 461 (2007).

48. Meehl, G. A., Washington, W. M., Wigley, T., Arblaster, J. M. & Dai, A. Solar and greenhouse gas forcing and climate response in the twentieth century. *J. Clim.* **16**, 426–444 (2003).

49. Martin, E. R., Thorncroft, C. & Booth, B. B. The multidecadal Atlantic SST—Sahel rainfall teleconnection in CMIP5 simulations. *J. Clim.* **27**, 784–806 (2014).

50. Yan, X., Zhang, R. & Knutson, T. R. Underestimated AMOC variability and implications for AMV and predictability in CMIP models. *Geophys. Res. Lett.* **45**, 4319–4328. https://doi.org/10.1029/2018GL077378 (2018).

51. Zhang, R. et al. Have aerosols caused the observed Atlantic multidecadal variability?. *J. Atmos. Sci.* **70**, 1135–1144 (2013).

52. Zhang, R. et al. Comment on “The atlantic multidecadal oscillation without a role for ocean circulation”. *Science* **352**, 1527–1527 (2016).

53. Zhang, R. On the persistence and coherence of subpolar sea surface temperature and salinity anomalies associated with the Atlantic multidecadal variability. *Geophys. Res. Lett.* **44**, 7865–7875 (2017).

54. Yan, X., Zhang, R. & Knutson, T. R. A multivariate AMV index and associated discrepancies between observed and CMIP5 externally forced AMV. *Geophys. Res. Lett.* **46**, 4421–4431. https://doi.org/10.1029/2019GL082787 (2019).

55. Eyring, V. et al. Overview of the coupled model intercomparison project phase 6 (CMIP6) experimental design and organization. *Geosci. Model Dev. (Online)* **9**, 1937–1958 (2016).
## Acknowledgements

We acknowledge the World Climate Research Programme's Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modeling groups for producing and making available their model output (listed in Table S1 of this paper). For CMIP the U.S. Department of Energy's Program for Climate Model Diagnosis and Intercomparison provides coordinating support and led development of software infrastructure in partnership with the Global Organization for Earth System Science Portals. We thank Haibo Liu for preparing the CMIP5 data for use and for helping with data access throughout this project. This research was supported by the U.S. National Science Foundation Grant AGS-1612904. A.G. benefited from French state support managed by the Agence Nationale de la Recherche under Investissements d'Avenir Program contract ANR-17-MPGA-0015, and M.B. and Y.K. received funding from DOE Grant DESC0014423. The Lamont publication number is 8420.

## Author contributions

A.G. conceived of the original question and the basic approach, and supervised the project along with M.B. and Y. Kushnir. R.J.H. developed the approach, designed the methods, and built the code. All authors contributed to writing the paper, with R.J.H. taking the lead.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** is available for this paper at https://doi.org/10.1038/s41598-020-68356-w.

Correspondence and requests for materials should be addressed to R.J.H.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020