Journal topic
Earth Syst. Dynam., 10, 631–650, 2019
https://doi.org/10.5194/esd-10-631-2019
Earth Syst. Dynam., 10, 631–650, 2019
https://doi.org/10.5194/esd-10-631-2019

Research article 18 Oct 2019

Research article | 18 Oct 2019

# Tipping the ENSO into a permanent El Niño can trigger state transitions in global terrestrial ecosystems

Tipping the ENSO into a permanent El Niño can trigger state transitions in global terrestrial ecosystems
Mateo Duque-Villegas, Juan Fernando Salazar, and Angela Maria Rendón Mateo Duque-Villegas et al.
• GIGA, Escuela Ambiental, Facultad de Ingeniería, Universidad de Antioquia, Medellín, Colombia

Correspondence: Juan Fernando Salazar (juan.salazar@udea.edu.co)

Abstract

Some large-scale components of the Earth's climate system have been identified as policy-relevant “tipping elements”, meaning that anthropogenic forcing and perturbations may push them across a tipping point threshold, with potential global scale impact on ecosystems and concomitant environmental and social phenomena. A pronounced change in the amplitude and/or frequency of the El Niño–Southern Oscillation (ENSO) is among such tipping elements. Here, we use the Planet Simulator (PlaSim), an Earth system model of intermediate complexity, to investigate the potential impact on global climate and terrestrial ecosystems of shifting the current dynamics of the ENSO into a permanent El Niño. When forced with sea surface temperature (SST) derived from observations, the PlaSim model yields a realistic representation of large-scale climatological patterns, including realistic estimates of the global energy and water balances, and gross primary productivity (GPP). In a permanent El Niño state, we found significant differences in the global distribution of water and energy fluxes, and associated impacts on GPP, indicating that vegetation production decreases in the tropics, whereas it increases in temperate regions. We identify regions in which these El Niño-induced changes are consistent with potential state transitions in global terrestrial ecosystems, including potential greening of western North America, dieback of the Amazon rainforest, and further aridification of south-eastern Africa and Australia.

1 Introduction

Earth's current climatology can be considered a steady state of a non-equilibrium system (Kleidon2010) in which constant influx of solar radiation and internal feedback mechanisms stall equilibrium, so that properties of the climate system can remain relatively close to mean values. However, it is becoming increasingly clear that the Earth's climate system as a whole, as well as some of its large-scale components or subsystems, have tipping points (i.e. critical thresholds in forcing and/or a feature of the system) at which the system shifts from one state (e.g. current climate) to another (e.g. future climate) . Perturbations related to human-induced climate change could push several Earth's subsystems past a tipping point , potentially changing the modes of natural climatic variability .

identified several tipping elements in the Earth's climate system. A tipping element is a large-scale component of the Earth system that may pass over a tipping point with the potential to alter global climate. A pronounced change in the amplitude and/or frequency of the El Niño–Southern Oscillation (ENSO), including its shift into a permanent El Niño state, is among the policy-relevant tipping elements identified by . Although there is no definitive evidence at present for specific changes in ENSO behaviour in response to climate change , global climate modelling studies do indicate the possibility of ENSO intensification during the 21st century due to anthropogenic forcing . Greenhouse warming in climate projections can lead to more frequent El Niño events, because it diminishes temperature gradients in the eastern equatorial Pacific Ocean, which facilitate the changes in convection zones that are observed during this phenomenon .

Through atmospheric teleconnections, ENSO exerts strong effects on manifold environmental and social systems at global to regional scales , including major implications for the functioning of terrestrial ecosystems . For instance, in Amazonia, temperature and precipitation changes associated with ENSO have been found to alter the net carbon flux to the atmosphere and seasonal inundation of floodplains . Moreover, El Niño-induced droughts in Amazonia increase the likelihood of forest fires that modify landscape features and have top-down effects that damage not only vegetation but also faunal populations . Similarly, in Australia's drylands, ENSO is a modulator of biomass removal through fires and grazing . Such land-use impacts are particularly important in places where the population largely depends on agriculture, such as south-eastern Africa, where ENSO also affects rainfall patterns . Therefore, and despite the inherent uncertainty of climate change projections, the potential impact of shifting the state of the ENSO is of global concern . Advancing the quantitative understanding of the range of expected physical climate changes (let alone tipping events) and their potential impact on different sectors and/or regions is crucial for decision-making related to global change .

Here, we used the Planet Simulator (PlaSim; Fraedrich et al.2005a, b) to investigate the potential impact on global climate of shifting the dynamics of the ENSO into a permanent El Niño state and examine the potential consequences of this tipping event on global terrestrial ecosystems. Our aim was to see what Earth's climate could be if the climatology of sea surface temperature (SST) were that observed during a warm phase of ENSO. Thus, we simulated a control scenario and compared it to current climate using data sets derived from observations, in order to understand the model's limitations. Then, we set up a permanent El Niño experiment and compared it to said control simulation, focusing on the impacts on water deficit and terrestrial ecosystems. We also briefly discuss how this permanent El Niño experiment relates to data and simulations from the warm Pliocene period, when it has been suggested that sustained El Niño conditions occurred . The exploratory nature of our modelling approach is more about identifying potential impacts than about estimating their likelihood.

2 Model description

PlaSim is a climate model of intermediate complexity. This kind of model is well suited for our goals because of its ability to parsimoniously simulate feedback mechanisms between many components of the climate system over long-term simulation ensembles covering a global domain . The dynamical core of PlaSim is based on PUMA , an atmospheric general circulation model that solves the moist primitive equations representing the conservation of momentum, mass, and energy in the Earth's atmosphere. It uses a vertical σ coordinate and the spectral transform method with triangular truncation over a sphere . Unresolved subgrid-scale processes such as radiation, moist and dry convection, large-scale precipitation, surface fluxes, and vertical and horizontal diffusion are parameterized .

The atmospheric model is coupled to a five-layer land surface model that includes parameterizations for soil temperature, snow cover, and a bucket-type runoff transport scheme . A dynamic vegetation model (SimBA) included in PlaSim allows to represent vegetation–atmosphere feedback mechanisms through climate-induced changes in large-scale surface parameters such as albedo, roughness, and water-holding capacity (Kleidon2006). The ocean can be represented through prescribed climatological values of SST or using dynamical models . For more details about the model description, we refer the reader to and .

PlaSim has been successfully used to study the evolution of the Earth's climate system under different physical scenarios and forcings. One of the first applications of PlaSim was to estimate the maximum possible influence of vegetation on the global climate through simulation of two extreme scenarios: “green planet” and “desert world” . More recent applications have studied the planetary climate system's response to, for instance, vegetation–atmosphere feedback mechanisms , large-scale changes in orography , and variations in astronomical forcings . The model has also been used for paleoclimate studies and snowball Earth experiments . Some of these studies have also explored the existence of alternative equilibrium states of the planetary climate system, as well as the mechanisms behind potential state transitions .

3 Experimental setup

For all simulations, PlaSim was implemented with a spectral horizontal resolution of T21 (about 5.6 on a Gaussian grid) and 10 non-equally-spaced σ vertical levels. Though relatively coarse, this resolution allows a realistic representation of large-scale atmospheric circulation and baroclinic structures that are important in the energy transport poleward in a slowly rotating and phase-locked planet like Earth . All simulations were integrated for 500 years at 40 min time steps, with model output being stored as daily averages. Such long-term simulations allow the model to approach an equilibrium state , which we characterized using the last 30 years of each simulation. The dynamic vegetation model SimBA was used to allow variations in terrestrial ecosystems. The solar constant value was set to 1365 W m−2. For atmospheric CO2, we used values of 360 and 415 ppm; however, differences among simulations were modest enough that here we present only the results obtained with the latter, which is closer to current estimates of this variable . All other model physical constants, planetary values and tuning parameters were kept at their default setting .

The ocean was represented through prescribed SST values that were used to define two scenarios: control (CTL) and permanent El Niño (PEN). Differences between scenarios are entirely determined by differences in the prescribed SST annual cycle (Figs. 1 and S1 in the Supplement). Due to the inherent chaotic nature of atmospheric dynamics, each scenario was studied using an equally weighted ensemble of three simulations with different initial conditions, which were varied through a random noise disturbance in the surface pressure field (e.g. Kunz et al.2009). Though we found little sensitivity to these variations in initial conditions (Fig. S2), we nevertheless used the ensemble means in our analyses. Average results from different scenarios were compared using a Student's t test at a significance level of α=0.05 and 30 degrees of freedom (corresponding to a 30-year simulation period).

Figure 1Annual mean of prescribed SST climatology in the equatorial Pacific Ocean region for CTL (a) and PEN (b) simulations, and their absolute differences (c).

## 3.1 Control scenario and model evaluation

The CTL scenario is based on present-day climatology and therefore was first used to assess the model's ability to realistically represent the Earth's climate system and then used as reference for comparison with the PEN scenario. The CTL scenario is characterized by a climatological annual cycle of SST values obtained from the second phase of the Atmospheric Model Intercomparison Project (AMIP II) for the period 1979–2010 . Regridding to model resolution and adjusting the least possible temperature to be 271.38 K was done following recommendations by .

Observational gridded data sets used to compare with the model results for the CTL scenario are 1987–2016 monthly values of precipitation from the Global Precipitation Climatology Project (GPCP; Adler et al.2003); 1988–2017 monthly values of near-surface air temperature from the Hadley Centre – Climatic Research Unit (HadCRUT4) temperature data set ; 1979–2004 monthly values of mean sea level pressure from the Hadley Centre Sea Level Pressure (HadSLP2) data set ; 2000–2015 monthly values of gross primary productivity (GPP) from the Moderate Resolution Imaging Spectroradiometer (MODIS) product MOD17A2 ; and 1987–2016 monthly values of free air temperature, as well as zonal and meridional winds at different pressure levels (1000–50 hPa) from the European Centre for Medium-Range Weather Forecasts (ECMWF) Interim reanalysis (ERA-Interim; Dee et al.2011). The global energy and water balances produced by PlaSim were compared with estimates based on ERA-Interim and those reported by . More details about all data sets used in this work are shown in Tables S1 and S2 as part of the Supplement.

## 3.2 Permanent El Niño

In this scenario, the prescribed climatological annual cycle of SST was replaced with the annual cycle of El Niño that occurred between June 2015 and May 2016, using monthly values from the AMIP II data set. This is the strongest El Niño recorded so far, surpassing strong events from last century: 1982–1983 and 1997–1998 . Figure 1c shows the annual mean differences between the PEN and CTL scenarios for the equatorial Pacific region. It can be seen that this El Niño event depicted well the first main mode of spatiotemporal variability of SST anomalies (the mature phase of ENSO) and did not exhibit the El Niño Modoki second main mode . The PEN climatology has a global annual warm bias of +0.4C, with both boreal winter (DJF) and autumn (SON) seasons having the largest positive deviations (Fig. S1). The warm bias is greater specifically in the equatorial Pacific Ocean, where it reaches an annual mean value of +0.8C.

suggest the possibility that permanent El Niño-like conditions have occurred during the warm early Pliocene period. Using environmental reconstructions from isotopes and bioindicators, they found that the zonal west-to-east gradient of SST stayed very close to the one observed during modern El Niño events. Although more recent studies suggest that their reconstruction may underestimate the Pacific warm pool temperature and its variability , it is still interesting to compare PEN scenario with data from the warm Pliocene period. Thus, we compared the SST forcing of PEN with the paleoenvironmental reconstruction data set PRISM3 , which covers said period and has been used in several Pliocene modelling experiments . PEN forcing is on average about 1 C warmer in the tropics, though differences in this region are not all significant (Fig. S3). PEN is cooler than PRISM3 elsewhere. The west-to-east zonal gradient was computed for CTL and PEN scenarios, as well as for PRISM3, at the sites shown in Fig. S3. For PEN, the gradient is 1.5 C, whereas for PRISM3 it is 1.8 C, and for CTL it is 3.1 C. Since this gradient is a good indicator of strength of the Walker circulation , it is expected that circulation will be weaker in PEN than in the mid-Pliocene and much weaker than in the CTL experiment.

## 3.3 Climatological water deficit estimation

In this study, we focus on the consequences that such a PEN climate scenario could have for terrestrial ecosystems. We used the maximum climatological water deficit (MCWD) as defined by , which is an indicator of accumulated water stress during dry seasons and has been linked to ecosystem degradation via hydraulic failure in plants , fires , and deforestation . We followed an approach similar to , who plotted MCWD versus mean annual precipitation (MAP) to study the possibility of land cover transitions under climate change projections in the Amazon rainforest. found that most savannah biomes in Amazonia (i.e. Cerrado) had rainfall rates below 1500 mm year−1 and accumulated water deficits less than −300 mm. Unlike , our aim was global and therefore, in order to be able to calculate MCWD, we could not use a single fixed evapotranspiration rate (ET) value for the full extent of our domain. Instead, we used annual mean ET for every grid point from the CTL scenario as a proxy of expected ET with current land cover. These reference values (Fig. S4) agree with the approximation used by of 100 mm per month average ET for Amazonia. Likewise, we assumed soil saturation at the wettest time of the year for each grid point and set the cumulative water deficit to zero for this month (Fig. S5).

Additionally, we judged useful to combine, in a single global map, statistically significant differences between scenarios in MAP and MCWD, in order to support our discussion about the impact that a PEN scenario could have on terrestrial ecosystems. For this map, we defined 10 % of the CTL value as a threshold to identify places where differences between scenarios were the largest. In this way, what we mean by “much greater/less” is that the difference exceeds the said threshold. Otherwise, even though both MAP and MCWD changed significantly for a particular grid point, this change was not as pronounced and only “greater/less” are used. Different colours represent the possible combinations of changes in both variables.

4 Results

## 4.1 Control climate

A comparison between PlaSim results for the CTL scenario and estimates from ERA-Interim and Trenberth et al. (2011, henceforth TFM11) shows that simulations produce realistic estimates of the global energy and water budgets (Fig. 2). At the top of the atmosphere (TOA), the magnitude of the net energy imbalance (i.e. the difference Δ between incoming and outgoing radiation) is lower for PlaSim ($\mathrm{\Delta }=+\mathrm{0.4}$ W m−2) than it is for ERA-Interim ($\mathrm{\Delta }=-\mathrm{1.4}$ W m−2) and TFM11 ($\mathrm{\Delta }=+\mathrm{0.9}$ W m−2). Previous studies have reported that PlaSim's energy imbalance is closer to zero than it is for many state-of-the-art climate models . Differences between PlaSim and observations indicate that the model overestimates the outgoing solar radiation at TOA, as a result of more solar radiation being reflected both at the surface and atmosphere, and less solar radiation being absorbed by the atmosphere. However, the model's estimate of solar radiation absorbed at the surface falls in between those of ERA-Interim and TFM11. The surface energy balance in the model shows realistic estimates of thermal, latent heat (LH), and sensible heat (SH) fluxes, with a slightly overestimated Bowen ratio (SH  LH) of 0.23 as compared to 0.21 in both ERA-Interim and TFM11. These biases in the surface fluxes are likely to be related to prescribed values of aerodynamic roughness length , which has also been related to ERA-Interim biases . The model's representation of the greenhouse effect through the estimate of thermal energy emitted by the atmosphere falls also within estimates from ERA-Interim and TFM11.

Figure 2Global average annual mean energy (a) and water (b) budgets for the CTL simulation compared with ERA-Interim and . Units are shown in the title of each panel.

As for the water budget, PlaSim presents a sensible distribution of water among ocean, land, and atmosphere, preserving key hydrological cycle characteristics: greater evaporation than precipitation over the oceans and greater precipitation than evapotranspiration over continental areas (Fig. 2b). As compared to TFM11, atmospheric transport from ocean to land, as well as continental runoff, are underestimated both in PlaSim and ERA-Interim, but unlike the reanalysis, the water balance for the model shows conservation. This is partly because the model values are not affected by any assimilation scheme, as is the reanalysis . Overall, PlaSim estimates of global water fluxes fall within the uncertainty ranges proposed by TFM11 for ERA-Interim and other reanalyses.

There are statistically significant differences between PlaSim and observational data in the global annual mean distribution of near-surface air temperature (Fig. 3c), precipitation rate (Fig. 3f), GPP (Fig. 3i) and mean sea level pressure (Fig. 3l). However, these biases are within typical ranges obtained for more comprehensive state-of-the-art climate models and cannot be entirely attributed to model deficiencies because of uncertainties in the observational data sets themselves (e.g. Prein and Gobiet2017). Importantly, large-scale spatial patterns of near-surface air temperature, precipitation rate, GPP, and mean sea level pressure (Fig. 3), as well as their seasonality (Figs. S6–S9), are realistically represented by PlaSim.

Figure 3Annual mean results for the CTL simulation compared with observational data: near-surface temperature in CTL (a) and HadCRUT4 (b) and model bias (c); total precipitation in CTL (d) and GPCP (e) and model bias (f); GPP in CTL (g) and MODIS (h) and model bias (i); and mean sea level pressure in CTL (j) and HadSLP2 (k) and model bias (l). In the bias panels, white depicts grid points with statistically non-significant difference (α=0.05). Grid lines are spaced every 30 from the Equator and 90 from Greenwich, in latitude and longitude, respectively.

Differences between PlaSim results and HadCRUT4 (Fig. 3a–c) indicate that the model exhibits a global average cold bias of about −0.4C. This global bias combines an average warm bias of +0.9C over continental areas and an average cold bias of −1.4C over the oceans. The global average annual mean temperature in the CTL simulation (13.9 C) is similar to that from observations (14.3 C). Seasonal variations of temperature are sensibly simulated too (Fig. S6). Some of the largest temperature biases occur on land near polar regions, especially in Greenland and Antarctica, and regions with complex terrain such as the Himalayas mountain range. Biases over complex terrain are related to known limitations even in state-of-the-art climate models, in part because the orography representation has to be smoothed to keep dynamical stability . Warm biases in western Africa and Australia and slightly cold biases in the Sahara desert and Amazonia are also observed. Strongest biases take place in boreal summer (JJA) and autumn (SON) seasons, when the global average biases on land are +1.7 and +1.6C, respectively. Over the oceans, temperatures from the model are colder than observations year-round, with a maximum of −1.6C in boreal winter (DJF).

Annual mean daily precipitation from the CTL simulation is depicted in Fig. 3d. The large-scale patterns exhibited in GPCP (Fig. 3e) are present: the equatorial maximum in the Intertropical Convergence Zone (ITCZ), middle latitude peaks of synoptic-scale weather systems, and low precipitation rates in the polar regions. There are high precipitation rates over the western equatorial Pacific and the Western Hemisphere warm pools. The PlaSim global average annual mean precipitation rate reaches 2.8 mm d−1, whereas for observational data, it is 2.7 mm d−1. Concerning seasonal variation (Fig. S7), tropical precipitation shows latitudinal displacement and has maximum peaks in south-eastern Asia in the boreal winter (DJF) and summer (JJA) seasons. In northern South America, the maximum of precipitation crosses the Equator twice a year, unlike over the eastern Pacific where the ITCZ stays north of the Equator for all seasons. PlaSim also displays the double ITCZ band that has been described for other global climate models (Lin2007), especially during the boreal spring season (MAM).

Figure 3g shows the dynamic vegetation results of CTL for GPP on continental areas. Global average annual mean for the CTL simulation reaches 0.7 kg C m−2 year−1, whereas for MODIS data it is 0.8 kg C m−2 year−1. Although large deviations with respect to MODIS can be seen in Fig. 3i, these are similar to those reported in other modelling studies and could be related to the climate data sets used as input in the simulations . Differences between this model and other data sets (whether derived from observations or not) occur because parameterization schemes to obtain GPP indirectly from climate variables are different. Nevertheless, the model shows a realistic estimate of the GPP distribution with high productivity near the Equator and low productivity elsewhere. Also the location of the subtropical deserts is consistent with observations (Fig. 3h). Furthermore, the model is able to represent seasonal changes in primary productivity (Fig. S8). There is a year-round high productivity around the Equator, whereas for the Northern Hemisphere and Southern Hemisphere it is at the maximum in their respective summer season and it is at the minimum in their respective winter season.

Figure 4Zonal average annual means results for the CTL simulation compared with ERA-Interim at different pressure levels for air temperature in CTL (a), ERA-Interim (b), and model bias (c); zonal wind in CTL (d), ERA-Interim (e), and model bias (f); and meridional mass stream function in CTL (g), ERA-Interim (h), and model bias (i). In the bias panels, white depicts grid points with statistically non-significant difference (α=0.05).

When compared with observations, PlaSim also displays a very similar annual mean pattern of high- and low-pressure cells in the subtropics (Fig. 3j and k). The largest differences are seen near Antarctica, where the model overestimates sea level pressure (Fig. 3l). Seasonally, the model is also able to represent the appearance of semi-permanent pressure systems throughout the year (Fig. S9). For instance, in the Northern Hemisphere, the Aleutian and Icelandic lows and the Pacific, Azores, and Siberian highs are seen in winter, whereas only the Pacific and Bermuda highs are seen in the summer season.

Regarding general circulation, the temperature gradient in the atmosphere for the CTL simulation has a very similar structure as the one from ERA-Interim reanalysis (Fig. 4a and b). Differences are relatively small below 500 hPa and away from the poles (Fig. 4c). Strongest biases are found around 200 hPa. In the upper troposphere, the model is much warmer in the tropics and much cooler near the North Pole. This behaviour is common to most general circulation models, which have difficulties parameterizing moist thermodynamic processes . Convective heating is overestimated in the Kuo-type convection scheme implemented in PlaSim (Arakawa2004).

The distribution of temperature is closely related to the zonal wind velocity (Fig. 4d and e). The location of the subtropical jet streams in the CTL simulation is in agreement with ERA-Interim, with the zonal wind speeds of the jet streams being higher in PlaSim. As with temperature, main differences are found in the upper troposphere and the lower stratosphere above 200 hPa (Fig. 4f). This is a known problem with some general circulation models that display too-fast zonal jet streams due to a coarse representation of the planet's orography and simple representation of the low-level drag . The trade winds in PlaSim prevail around the Equator but are slower than in the reanalysis below 300 hPa and faster above. As compared to ERA-Interim, zonal wind speed is underestimated eastward in PlaSim at 60 S at all pressure levels.

PlaSim adequately simulates the thermally induced Hadley cell and its displacement across the Equator for the annual mean (Fig. 4g and h), with clockwise rotation when in the Southern Hemisphere and anti-clockwise in the Northern Hemisphere. As described for observations, when in the Southern Hemisphere, the cell is stronger. Similarly, the indirect Ferrell cell can be seen moving in the opposite direction as the former (different colour in Fig. 4g and h). Differences in the stream function occur for the Hadley cell at both hemispheres, which is slower and shorter in PlaSim (Fig. 4i). Because of this, the indirect Ferrell cell is also smaller and is almost absent when in the Northern Hemisphere. Polar cells are present but not visible in any of these figures because of the chosen contour levels.

Altogether, the PlaSim CTL simulation realistically represents the broad-scale patterns seen in observational data of present-day climate. In spite of PlaSim's simplicity and coarse resolution of simulations, the model produces a sensible picture of the current climatological state of the Earth's climate system. Our results generally agree with previous studies showing that PlaSim is able to simulate many dominant features of the present-day climate . Global averages of annual mean and seasonal mean values for CTL and observational data sets are also summarized in Table S3.

## 4.2 Permanent El Niño climate

To be able to judge this scenario on account of its different ENSO state with respect to CTL, we checked that model bias was similar in PlaSim during El Niño conditions and CTL. Consequently, we compared the PEN scenario to the observational data sets during very strong El Niño events, except for MODIS data which only include moderate ones (Fig. S10). Model bias for PEN variables is very similar to that shown in Fig. 3, and also there is good agreement in the broad-scale patterns seen in PEN and observations during El Niño events. Figure S11 shows more clearly that model bias is very similar in both experiments; hence, we can attribute most of the differences between CTL and PEN to their ENSO condition.

Warmer-than-usual equatorial Pacific Ocean temperatures, as set in the PEN scenario, lead to statistically significant changes in the global energy and water budgets (Fig. 5). For the energy balance, a warmer ocean surface results in higher long-wave radiation being emitted, as seen in the heat fluxes and the thermal upward radiation. It also diminishes surface snow and sea ice formation, which is followed by a low planetary effective albedo, and therefore less solar incoming radiation is reflected. Like the ocean, the atmosphere also becomes warmer and its thermal emission increases. Ultimately, this leads to less incoming solar radiation being reflected and more long-wave radiation being emitted back to space. These two processes do not balance each other in the model, and therefore the net imbalance at TOA becomes slightly more negative in PEN than in the CTL simulation (Fig. 5a). As expected, the PEN scenario shows much warmer surface and atmosphere than CTL.

Figure 5Global average annual mean energy (a) and water (b) budgets for the PEN simulation compared with the CTL simulation. Units are shown in the title of each panel. En dashes (–) show statistically non-significant differences (α=0.05).

According to PEN simulation outputs, water is redistributed throughout the world (Fig. 5b). A warmer atmosphere stores more water, but the net transport of water from the ocean to land is reduced. This may be explained by changes in wind patterns over the Pacific Ocean, linked to the Walker circulation, which keep most water vapour above the ocean surface, as is typical during El Niño conditions . Consistently, evaporation and precipitation increase over the ocean, while they decrease on land. The values in Fig. 5 are also presented in tabular form in Tables S4 and S5, including the results of the significance test.

There are statistically significant differences in the near-surface air temperatures for the PEN scenario as compared to the CTL simulation (Fig. 6a). PEN has an average global annual mean increase of approximately 0.5 C. This temperature change is greater on land with an average 0.8 C increase, than over the oceans, where the change reaches 0.5 C. In the tropics, variation occurs mainly in the equatorial Pacific region (as expected) but also in Central America, northern South America, Amazonia, the Sahel region in Africa, and northern Australia. Moreover, air over the North Atlantic Ocean is slightly colder, while the Arctic Ocean near Canada and Siberia is warmer. Air over some ocean areas in the Southern Hemisphere, close to Antarctica, becomes colder in PEN than in CTL. No significant differences in temperature were seen in most of North America. Seasonally, boreal winter (DJF) shows the greatest warming on land, especially for Australia (Fig. S12). South America is the warmest during boreal summer (JJA).

Figure 6Differences in annual means results for the PEN simulation compared with CTL simulation (PEN minus CTL) for near-surface temperature (a), daily precipitation (b), GPP (c), and mean sea level pressure (d). In all panels, white depicts grid points with statistically non-significant difference (α=0.05). Grid lines are spaced every 30 from the Equator and 90 from Greenwich, in latitude and longitude, respectively.

Regarding precipitation, PEN leads to an average rise in global mean precipitation of approximately 0.1 mm d−1 and a decrease in continental areas of about the same value (Fig. 6b). Precipitation over the tropical Atlantic Ocean is also slightly reduced in the PEN scenario. On land, the most conspicuous differences occur near the Equator, in Central America and in Amazonia, where precipitation is reduced. Conversely, in the Northern Hemisphere above midlatitudes, land areas experience slightly higher precipitation rates. Similar to temperature, boreal winter (DJF) exhibits the largest deviations (Fig. S12), with above- (below)-CTL rainfall in equatorial eastern (western) Pacific Ocean.

Figure 7Meridional divergent circulation in the equatorial Pacific for the CTL simulation (a), climatological ERA-Interim (b), and CTL bias in vertical velocity values (c), as well as for the PEN simulation (d), ERA-Interim mean of 1997–1998 and 2015–2016 (e), and PEN bias in vertical velocity values (f). Vectors are plotted using the divergent meridional wind component and the negative vertical velocity in pressure units (ω), averaged over latitudes 5 S to 5 N. Filled contours show the magnitude of the vertical velocity. In all bias panels, white is for statistically non-significant differences (α=0.05).

Higher temperatures and precipitation rates in temperate zones result in higher productivity rates (Fig. 6c). The opposite happens in regions where temperatures increase, while precipitation rates decrease, as in Central America and northern South America and Amazonia. In Africa, near the Middle East, there is a relatively small increase in productivity related to increased precipitation. On the contrary, below the Sahel region and south-eastern Africa, GPP diminishes where precipitation declined with respect to the CTL simulation. The same happens in Australia. Globally, a decrease in productivity of about 0.1 kg C m−2 year−1 occurs in the PEN scenario. PEN results also show decreased productivity in South America year-round (Fig. S13).

PEN simulation shows the Southern Oscillation hemispheric zonal dipole in surface pressure (Fig. 6d). In the tropics, there is an overall increase in surface pressure over the oceans (except in the eastern Pacific), which is related to drier air and reduced transport of moisture into continental areas. As expected, lower (higher) surface pressure conditions remain year-round in eastern (western) equatorial Pacific Ocean (Fig. S13). The North Atlantic region has lower surface pressure throughout the year but most markedly in boreal winter (DJF).

The atmosphere becomes warmer in PEN than it is in the CTL simulation at most pressure levels and latitudes (Fig. S14a). There is only a slight cold bias in the Northern Hemisphere above 200 hPa. The hotspot in the tropics in the upper troposphere is much warmer in the PEN simulation. This result has been previously reported and seems to be related to changes in the radiative forcing of the troposphere layer . Wind speed increases eastward (or lessens westward) in the model atmosphere (Fig. S14b). Notably, in the midlatitudes in the middle troposphere, the zonal velocity rises eastward, moving the subtropical jet streams towards the Equator as commonly seen during El Niño conditions. As for the Hadley and Ferrell cells (Fig. S14c), their magnitude increases slightly in both hemispheres, as has been previously stated .

Figure 8Annual mean results for the PEN simulation compared with the CTL simulation for MAP in CTL (a) and PEN (b) and their absolute difference (c), and MCWD in CTL (d) and PEN (e) and their absolute difference (f). In the difference panels, white depicts grid points with statistically non-significant difference (α=0.05). Grid lines are spaced every 30 from the Equator and 90 from Greenwich, in latitude and longitude, respectively.

The Walker circulation in the equatorial Pacific is shown in Fig. 7 for CTL and PEN simulations. These are compared to ERA-Interim climatology and ERA-Interim during very strong El Niño conditions, respectively. We observe very similar broad-scale structures in simulations and the reanalysis. In CTL, air rises in the west Pacific near longitude 150 E, flows eastward and sinks around longitude 120 W. This is the expected behaviour during normal conditions and is also displayed in ERA-Interim climatology. It seems like PlaSim underestimates slightly both upward and downward motions as compared to the reanalysis. In the PEN simulation, convection is displaced towards the central Pacific, while subsidence in the east is decreased. This is very similar to what happens in ERA-Interim for very strong El Niño events (mean of 1997–1998 and 2015–2016), although in the reanalysis subsidence in the eastern Pacific ceases almost completely. Overall, PlaSim is able to adequately show the weakening of the Walker circulation during El Niño-like conditions.

## 4.3 Climatological water deficit

Results for both MAP and MCWD are shown in Fig. 8. In terms of MAP, the PEN simulation clearly shows rainfall decline in the tropics, whereas it increases in northern temperate latitudes. A similar situation is observed for MCWD. As described in Sect. 3.3, we combined the differences in MAP and MCWD between scenarios in a single map that is shown in Fig. 9. This map resembles others that display the global effects of El Niño (e.g. Holmgren et al.2001). According to these results, tropical ecosystems seem to be the most vulnerable to a PEN scenario, whereas northern temperate ecosystems could enhance their primary productivity due to greater water availability. We identify the following regions as the most vulnerable: Central America, northern and central South America, south-eastern Africa, and Australia. On the other hand, some parts of western North America and a wide band in Eurasia have greater precipitation, and cumulative water deficit has diminished. Figure 10 shows how conditions are different in CTL and PEN scenarios, in terms of both MAP and MCWD, in each of these world regions. In this figure, for every grid point in a region of interest in Fig. 9, markers are placed at their PEN value and a solid line shows the distance to what they are in the CTL scenario. Most grid points in the vulnerable regions (Fig. 10b, c, e, and f) shift towards much less precipitation conditions with greater water deficits (see the shift to the left). Using a similar approach as for Amazonia, but in a global context near the tropics, these results indicate that a PEN scenario can shift terrestrial ecosystems towards more arid conditions. In the case of ecosystems in North America and Eurasia (Fig. 10a and d), a PEN scenario drives them towards wetter conditions (see the shift to the right).

Figure 9Map of grid points where the PEN scenario had the largest differences in terms of both MAP and MCWD, with respect to the CTL simulation. Grid lines are spaced every 30 from the Equator and 90 from Greenwich, in latitude and longitude, respectively.

## 4.4 Warm Pliocene period

The possibility that sustained El Niño conditions might have occurred during the Pliocene epoch compels us to compare the PEN experiment with previous simulations of the warm Pliocene period. simulated this period (using a CO2 concentration of 400 ppm) and found global surface temperature to rise. Similarly, near-surface temperature increases in PEN scenario with respect to CTL. In both cases, global warming is greater on land than on the oceans. Regarding average precipitation rate, it increases slightly both in PEN and in the Pliocene simulation with respect to their control scenarios. Another interesting correspondence is the fact that in and in PEN, annual mean precipitation decreases in Central America and South America (tropical regions), whereas it increases in the Northern Hemisphere. Particularly for Amazonia, both studies show that average temperature rises and precipitation decreases. report that in the mid-Pliocene warm period there might have been strengthening of the meridional circulation, leading to aridification of Africa and Australia, which is similar to what we observe in PEN results. However, they also report possible aridification of North America, and this is the opposite of what is seen in PEN (Fig. 9). In spite of the similarities between PEN and the mid-Pliocene warm period simulations, a biome reconstruction of the Late Pliocene period seems to disagree with our results, specially for Africa and Australia, since it shows more vegetated regions than we would expect with the water deficit changes observed in PEN.

5 Discussion

Global terrestrial ecosystems can have alternative stable equilibrium sates and can shift between them when external perturbations and/or internal dynamics force them past critical thresholds . Alternative equilibria have already been identified for some biomes in the world. In tropical South America, in the eastern Amazon River basin, where there is currently a rainforest, simulation results suggest the presence of an alternative savannah (Cerrado) state . In both the Sahara and Sahel regions in northern Africa, geological records and modelling results suggest that environmental conditions could have supported a “green Sahara” and a “wet Sahel”, but climate shifts and atmosphere–vegetation feedback mechanisms could have stabilized their current contrasting states . Similarly for the wet tropics in Australia, bi-stability between rainforest and savannah-like sclerophyll land cover has also been studied . Forest, savannah, and desert seem to be three possible alternative stable states for terrestrial ecosystems .

Figure 10Changes in MAP and MCWD in the PEN scenario compared to CTL for different world regions. For every grid point shown in the zoomed map axes, we use markers (triangles) to show the value of both variables in the PEN scenario, and a straight line shows the distance to the value they have in the CTL scenario.

studied the possibility of a climate-change-induced state transition in the Amazon rainforest. Estimating water availability and accumulated water stress along a dry season, they identified transition areas between rainforest and savannah alternative states for this large-scale ecosystem. Here, we extended this approach to other global ecosystems in order to evaluate whether a PEN scenario could trigger state transitions not only in Amazonia but also in other global ecosystems. The most noticeable effects suggested by the PlaSim model include a potential greening of western North America, dieback of the Amazon rainforest, drying of south-eastern Africa, and further aridification of Australia.

Increased precipitation in parts of North America is a known effect of strong El Niño events and can be seen in PEN results, as well as reductions in water deficit. In Fig. 10a, we observe that significant changes in these variables (with respect to CTL) were greater in grid points around the Great Basin region, which is an endorheic watershed that includes arid and semi-arid biomes, and is located close to the Mojave and Sonoran deserts. Several stable states for vegetation in said ecosystems include shrublands, short-grass steppes, and desert grasslands (Laycock1991). Recently, it has been found that El Niño is linked to increased primary production in North America , as well as to growth and persistence of vegetation in these desert areas . Thus, in a PEN scenario, greater water availability could help avoid further aridification in this region and even trigger greening of some parts of it.

The opposite situation occurs in Central America and northern South America as seen in Fig. 10b. In this case, there is less annual precipitation and the water deficit increases. Such results agree with observations during El Niño events which show negative anomalies in several hydrological variables like rainfall and river flows for tropical South America and Central America . In spite of many factors being at play, including the meridional displacement of the ITCZ, complex topography, and two different neighbouring oceans, ENSO has been found to be the greatest cause of interannual variability in this region , and therefore sustained El Niño conditions are very likely to have a noticeable impact as depicted in Fig. 10b. Since this region is considered a biodiversity hotspot in the world , ecosystem degradation could be related to a demise in its number of biological species due to the drastic changes in water availability observed in PEN. Moreover, hydropower generation is a common practice in this part of the world and a change of over 10 % in annual precipitation could hinder energy production.

Our results show that PEN can be a driver of forest dieback in the Amazon (Fig. 10c). El Niño reduces annual precipitation and increases the water deficit across the dry season for this region . This rainfall reduction could perhaps threaten the existence of parts of the rainforest because resilience to dry conditions is limited and a long-lasting drought can increase tree mortality . Most forest species operate within a narrow margin of hydraulic safety and therefore are vulnerable to drought-induced embolisms which reduce productivity and increase the risk of forest decline . Additionally, El Niño-induced fires (not simulated by PlaSim) could further exacerbate forest degradation or transition from evergreen to seasonal tropical forests . Consequences of this are manifold, but one important aspect to consider is the global impact the Amazon rainforest has as a carbon sink . A PEN-induced dieback in the Amazon could tip the rainforest from a sink of carbon to a possible source . This could in turn amplify the magnitude of extreme hydrological and weather events that occur during El Niño, such as droughts, in a positive-feedback-like manner , as well as disrupt mechanisms of river flow regulation and continental precipitation distribution . Near-synchronous effects of this sort between different tipping elements could give rise to a “tipping cascade” that could contribute to a planetary-scale state shift . Despite potential forest collapse in most of the Amazon, our results also agree with that some areas in western Amazonia (adjacent to the Andes) are less vulnerable and could act as a refuge for biodiversity. Consequences of a PEN scenario in the Amazon rainforest would have global impact, since it is the largest remaining intact forest in the world, which provide globally significant environmental benefits .

In Eurasia (Fig. 10d), water availability rises in the PEN scenario across central Europe, the Middle East around the Caspian Sea, and central Asia. This is a known feature of El Niño which increases precipitation rates in the region . Increased primary production is another behaviour observed during El Niño years in Eurasia that is present in PEN. Particularly in central Asia, a PEN scenario shows wetter conditions near Karakum and Kyzylkum deserts, where it has been found recently that shrublands are extremely sensitive to climatic variations and that drought is the key driver of vegetation loss . Therefore, we believe that under sustained El Niño conditions, desertification processes in central Asia could be inhibited or slowed down. Compared to North America, changes are not as pronounced in precipitation and water deficit; hence, we believe greening mechanisms are less likely for Eurasia. The fact that both North America and Eurasia exhibit similar results in the PEN scenario can be related to warmer winter and spring seasons and a cooler summer season that leads to greater primary productivity in both places, which is observed in temperate forests during El Niño years .

The influence of ENSO on rainfall patterns in south-eastern Africa was identified long ago . This happens directly via atmospheric teleconnections and also through warming of the Indian and Atlantic oceans . El Niño in this region is linked to below average rainfall and droughts . This is present in our results south of the Great Horn of Africa, where the PEN scenario creates a greater water deficit and much less precipitation (Fig. 10e). Our results agree with , who discussed the impact on Africa of increased frequency of El Niño episodes using climate projections and state-of-the-art general circulation models. They also found south-eastern Africa to be particularly vulnerable. Earlier modelling studies found that such savannah ecosystems in Africa are closely related to climate variability, which keeps the system in between desert-like and forest-like equilibrium states . More recent work also associates the savannah equilibrium state to rainfall variability . Hence, tipping the ENSO permanently could trigger land cover transitions from forest to savannah and/or savannah to desert-like biomes. This is problematic not only in terms of ecosystem degradation but also in terms of putting at risk food security and hydropower generation (Lyon2014). In the case of west Africa, the impacts of ENSO have been contested because of the many other climatic drivers that come together in this area . However, ENSO events have been linked to climatic variability in both the Sahel and Guinea coastal regions . These two regions are connected via the west African monsoon in boreal summer, which brings moisture inland from the Atlantic Ocean . During El Niño conditions, a coherent drying pattern for the Sahel has been identified . PlaSim simulation results also support the role of the warm ENSO phase in modulating rainfall in west Africa and show it leads to a drier state in southern Sahel and Guinea coast regions (Fig. 8c). This southward expansion of the Sahel dry conditions into Guinea Coast, agrees with the notion of potential enlargement of the world's deserts due to climate change . Though this region only has a single grid point highlighted in Fig. 9 (due to small changes in water deficit), we believe a PEN scenario could also push its arid and semi-arid ecosystems into vegetation–atmosphere positive feedback mechanisms of desertification.

El Niño events have had different effects on rainfall in Australia, depending on how far east the Walker circulation moves . In this case, we show that a PEN based on the strong El Niño of 2015–2016 is able to alter climate with potentially severe ecological consequences (Fig. 10f). Likewise, in spite of observed weakening of teleconnections between Australia rainfall and ENSO , our results indicate that a PEN scenario could cause significant widespread aridification of Australia. Coupling this process with land cover changes could give rise to synergistic effects that could potentially tip Australia's climate towards a much drought-prone state . Furthermore, extreme weather events (such as droughts) and changes in local hydrology have been identified as key drivers of terrestrial and marine Australian ecosystems towards tipping points . Hence, following our results, a PEN could originate a “tipping cascade” in this continent.

Thus, the response of tropical and temperate ecosystems to PEN seems to be asymmetric: while tropical ecosystems are degraded (GPP is generally reduced), temperate ecosystems are enhanced (GPP increases). This, combined with land cover changes (human or PEN induced), could turn tropical biomes into global net sources of carbon that, according to current estimates , could not be counterbalanced by extratropical forests.

6 Conclusions

This study has provided additional evidence that PlaSim can realistically represent the current state of the Earth's climate system. Despite its coarse resolution and simplified parameterizations, the model is able to display the large-scale patterns of atmospheric dynamics seen in observational data sets. In describing some of its limitations, we found that it displays common biases seen in other state-of-the-art climate models. Therefore, PlaSim provides an adequate framework to study the dynamics of the Earth system (as a whole), including complex climate feedback mechanisms, using long simulations.

Our key finding is that a permanent tipping of the ENSO into a very strong warm phase (El Niño) can trigger state transitions in global terrestrial ecosystems, with asymmetric effects between the tropics and extratropical regions. We identified some regions that seem to be more susceptible to shifting into more arid biomes: Central America, the Amazon rainforest, south-eastern Africa, and Australia, whereas parts of North America and Eurasia could become more productive. These results also suggest the possibility that a permanent El Niño can trigger tipping cascades. It is important to mention that despite there being other sources of variability for these ecosystems (for instance, the Atlantic Ocean temperatures variation), here, we focused only on the consequences that an observed very strong El Niño-based SST climatology could have on continental areas, without linking any simulated anomalies to a specific region in space. This means that despite the Pacific Ocean showing the largest anomalies (as in El Niño), we cannot attribute all of our results to changes only in this region. Even though a modelling study of this sort is not enough to prove that such state shifts will occur in a climate change scenario, it does show a consistent picture and should raise concern about the conservation of global ecosystems and related water resources management practices in a world with more frequent and/or intense El Niño events .

Data availability
Data availability.

All observational data sets and reanalyses used in this study are publicly available online. Model results are readily available upon request.

Supplement
Supplement.

Author contributions
Author contributions.

MDV performed the simulations. All authors participated in the design of the experiments, the discussion of the results, and the writing of the paper.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

We thank Daniel Ruiz Carrascal and Luis F. Salazar for insightful comments on an early version of the manuscript. We also thank Stephen Sitch and Lina Mercado for their valuable remarks. Finally, we thank Rui A. P. Perdigão and two anonymous reviewers for their helpful observations.

Financial support
Financial support.

This research has been supported by the Universidad de Antioquia (Estudiante Instructor program) and the UT-GRA (Programa de investigación en la gestión de riesgo asociado con cambio climático y ambiental en cuencas hidrográficas) (grant no. 543-2011 Colciencias).

Review statement
Review statement.

This paper was edited by Rui A. P. Perdigão and reviewed by two anonymous referees.

References

Adler, R. F., Huffman, G. J., Chang, A., Ferraro, R., Xie, P.-P., Janowiak, J., Rudolf, B., Schneider, U., Curtis, S., Bolvin, D., Gruber, A., Susskind, J., Arkin, P., and Nelkin, E.: The Version-2 Global Precipitation Climatology Project (GPCP) Monthly Precipitation Analysis (1979–Present), J. Hydrometeorol., 4, 1147–1167, https://doi.org/10.1175/1525-7541(2003)004<1147:TVGPCP>2.0.CO;2, 2003. a

Alencar, A., Nepstad, D., and Diaz, M. C. V.: Forest understory fire in the Brazilian Amazon in ENSO and non-ENSO years: area burned and committed carbon emissions, Earth Interact., 10, 1–17, https://doi.org/10.1175/EI150.1, 2006. a

Allan, R. and Ansell, T.: A New Globally Complete Monthly Historical Gridded Mean Sea Level Pressure Dataset (HadSLP2): 1850–2004, J. Climate, 19, 5816–5842, https://doi.org/10.1175/JCLI3937.1, 2006. a

Andersen, T., Carstensen, J., Hernandez-Garcia, E., and Duarte, C. M.: Ecological Thresholds And Regime Shifts: Approaches To Identification, Trends Ecol. Evol., 24, 49–57, https://doi.org/10.1016/j.tree.2008.07.014, 2009. a

Aragao, L. E. O., Malhi, Y., Roman-Cuesta, R. M., Saatchi, S., Anderson, L. O., and Shimabukuro, Y. E.: Spatial Patterns And Fire Response Of Recent Amazonian Droughts, Geophys. Res. Lett., 34, L07701, https://doi.org/10.1029/2006GL028946, 2007. a

Arakawa, A.: The cumulus parameterization problem: past, present, and future, J. Climate, 17, 2493–2525, 2004. a

Ardö, J.: Comparison between remote sensing and a dynamic vegetation model for estimating terrestrial primary production of Africa, Carbon Balance Manage., 10, 1–15, 2015. a

Ashcroft, L., Gergis, J., and Karoly, D. J.: Long-term Stationarity Of El Niño–Southern Oscillation Teleconnections In Southeastern Australia, Clim. Dynam., 46, 2991–3006, https://doi.org/10.1007/s00382-015-2746-3, 2016. a

Ashok, K., Behera, S. K., Rao, S. A., Weng, H., and Yamagata, T.: El Niño Modoki and its possible teleconnection, J. Geophys. Res.-Oceans, 112, 1–27, 2007. a

Barlow, J. and Peres, C. A.: Ecological Responses To El Niño-Induced Surface Fires In Central Brazilian Amazonia: Management Implications For Flammable Tropical Forests, Philos. T. Roy. Soc. B, 359, 367–380, https://doi.org/10.1098/rstb.2003.1423, 2004. a

Barnosky, A. D., Hadly, E. A., Bascompte, J., Berlow, E. L., Brown, J. H., Fortelius, M., Getz, W. M., Harte, J., Hastings, A., Marquet, P. A., Martinez, N. D., Mooers, A., Roopnarine, P., Vermeij, G., Williams, J. W., Gillespie, R., Kitzes, J., Marshall, C., Matzke, N., Mindell, D. P., Revilla, E., and Smith, A. B.: Approaching A State Shift In Earth's Biosphere, Nature, 486, 52–58, https://doi.org/10.1038/nature11018, 2012. a, b

Bathiany, S., Claussen, M., and Fraedrich, K.: Implications Of Climate Variability For The Detection Of Multiple Equilibria And For Rapid Transitions In The Atmosphere-Vegetation System, Clim. Dynam., 38, 1775–1790, https://doi.org/10.1007/s00382-011-1037-x, 2012. a, b

Betts, R. A., Malhi, Y., and Roberts, J. T.: The Future Of The Amazon: New Perspectives From Climate, Ecosystem And Social Sciences, Philos. T. Roy. Soc. B, 363, 1729–1735, https://doi.org/10.1098/rstb.2008.0011, 2008. a

Boschi, R., Lucarini, V., and Pascale, S.: Bistability Of The Climate Around The Habitable Zone: A Thermodynamic Investigation, Icarus, 226, 1724–1742, https://doi.org/10.1016/j.icarus.2013.03.017, 2013. a, b

Brando, P. M., Balch, J. K., Nepstad, D. C., Morton, D. C., Putz, F. E., Coe, M. T., Silvério, D., Macedo, M. N., Davidson, E. A., Nóbrega, C. C., Alencar, A., and Soares-Filho, B. S.: Abrupt Increases In Amazonian Tree Mortality Due To Drought-Fire Interactions, P. Natl. Acad. Sci. USA, 111, 6347–6352, https://doi.org/10.1073/pnas.1305499111, 2014. a, b

Brovkin, V., Claussen, M., Petoukhov, V., and Ganopolski, A.: On The Stability Of The Atmosphere-Vegetation System In The Sahara/Sahel Region, J. Geophys. Res.-Atmos., 103, 31613–31624, https://doi.org/10.1029/1998JD200006, 1998. a

Buermann, W., Anderson, B., Tucker, C. J., Dickinson, R. E., Lucht, W., Potter, C. S., and Myneni, R. B.: Interannual covariability in Northern Hemisphere air temperatures and greenness associated with El Niño–Southern Oscillation and the Arctic Oscillation, J. Geophys. Res.-Atmos., 108, 4396, https://doi.org/10.1029/2002JD002630, 2003. a

Cai, W., Borlace, S., Lengaigne, M., Van Rensch, P., Collins, M., Vecchi, G., Timmermann, A., Santoso, A., McPhaden, M. J., Wu, L., England, M. H., Wang, G., Guilyardi, E., and Jin, F.-F.: Increasing Frequency Of Extreme El Niño Events Due To Greenhouse Warming, Nat. Clim. Change, 4, 111–116, https://doi.org/10.1038/nclimate2100, 2014. a, b, c

Cai, W., Santoso, A., Wang, G., Yeh, S.-W., An, S.-I., Cobb, K. M., Collins, M., Guilyardi, E., Jin, F.-F., Kug, J.-S., Lengaigne, M., McPhaden, M. J., Takahashi, K., Timmermann, A., Vecchi, G., Watanabe, M., and Wu, L.: ENSO And Greenhouse Warming, Nat. Clim. Change, 5, 849–859, https://doi.org/10.1038/nclimate2743, 2015. a, b, c

Choat, B., Jansen, S., Brodribb, T. J., Cochard, H., Delzon, S., Bhaskar, R., Bucci, S. J., Feild, T. S., Gleason, S. M., Hacke, U. G., Jacobsen, A. L., Lens, F., Maherali, H., Martínez-Vilalta, J., Mayr, S., Mencuccini, M., Mitchell, P. J., Nardini, A., Pittermann, J., Pratt, R. B., Sperry, J. S., Westoby, M., Wright, I. J., and Zanne, A. E.: Global Convergence In The Vulnerability Of Forests To Drought, Nature, 491, 752–755, https://doi.org/10.1038/nature11688, 2012. a

Claussen, M., Mysak, L., Weaver, A., Crucifix, M., Fichefet, T., Loutre, M.-F., Weber, S., Alcamo, J., Alexeev, V., Berger, A., Calov, R., Ganopolski, A., Goosse, H., Lohmann, G., Lunkeit, F., Mokhov, I., Petoukhov, V., Stone, P., and Wang, Z.: Earth System Models Of Intermediate Complexity: Closing The Gap In The Spectrum Of Climate System Models, Clim. Dynam., 18, 579–586, https://doi.org/10.1007/s00382-001-0200-1, 2002. a

Collins, M., An, S.-I., Cai, W., Ganachaud, A., Guilyardi, E., Jin, F.-F., Jochum, M., Lengaigne, M., Power, S., Timmermann, A., Vecchi, G., and Wittenberg, A.: The impact of global warming on the tropical Pacific Ocean and El Niño, Nat. Geosci., 3, 391–397, https://doi.org/10.1038/ngeo868, 2010. a

Costa, M. H. and Pires, G. F.: Effects Of Amazon And Central Brazil Deforestation Scenarios On The Duration Of The Dry Season In The Arc Of Deforestation, Int. J. Climatol., 30, 1970–1979, https://doi.org/10.1002/joc.2048, 2010. a

Cox, P. M., Betts, R., Collins, M., Harris, P. P., Huntingford, C., and Jones, C.: Amazonian Forest Dieback Under Climate-Carbon Cycle Projections For The 21st Century, Theor. Appl. Climatol., 78, 137–156, https://doi.org/10.1007/s00704-004-0049-4, 2004. a, b

Davidson, E. A., de Araújo, A. C., Artaxo, P., Balch, J. K., Brown, I. F., Bustamante, M. M., Coe, M. T., DeFries, R. S., Keller, M., Longo, M., Munger, J. W., Schroeder, W., Soares-Filho, B. S., Souza, C. M., and Wofsy, S. C.: The Amazon Basin In Transition, Nature, 481, 321–328, https://doi.org/10.1038/nature10717, 2012. a

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim Reanalysis: Configuration And Performance Of The Data Assimilation System, Q. J. Roy. Meteorol. Soc., 137, 553–597, https://doi.org/10.1002/qj.828, 2011. a

Dekker, S., De Boer, H., Brovkin, V., Fraedrich, K., Wassen, M., and Rietkerk, M.: Biogeophysical Feedbacks Trigger Shifts In The Modelled Vegetation–Atmosphere System At Multiple Scales, Biogeosciences, 7, 1237–1245, https://doi.org/10.5194/bg-7-1237-2010, 2010. a, b

Diaz, H. F. and Markgraf, V.: El Niño and the Southern Oscillation: Multiscale variability and global and regional impacts, Cambridge University Press, Cambridge, 2000. a, b

Dlugokencky, E. and Tans, P.: Globally averaged marine surface monthly mean data, available at: https://www.esrl.noaa.gov/gmd/ccgg/trends/, last access: 8 August 2019. a

Dowsett, H., Robinson, M., and Foley, K.: Pliocene three-dimensional global ocean temperature reconstruction, Clim. Past, 5, 769–783, https://doi.org/10.5194/cp-5-769-2009, 2009. a

Ellison, D., Morris, C. E., Locatelli, B., Sheil, D., Cohen, J., Murdiyarso, D., Gutierrez, V., Van Noordwijk, M., Creed, I. F., Pokorny, J., Gaveau, D., Spracklen, D. V., Bargués Tobella, A., Ilstedt, U., Teuling, A. J., Gebrehiwot, S. G., Sands, D. C., Muys, B., Verbist, B., Springgay, E., Sugandi, Y., and Sullivan, C. A.: Trees, forests and water: Cool insights for a hot world, Global Environ. Change, 43, 51–61, https://doi.org/10.1016/j.gloenvcha.2017.01.002, 2017. a

Fedorov, A., Brierley, C., Lawrence, K. T., Liu, Z., Dekens, P., and Ravelo, A.: Patterns and mechanisms of early Pliocene warmth, Nature, 496, 43–49, https://doi.org/10.1038/nature12003, 2013. a

Flato, G., Marotzke, J., Abiodun, B., Braconnot, P., Chou, S. C., Collins, W. J., Cox, P., Driouech, F., Emori, S., Eyring, V., Forest, C., Gleckler, P., Guilyardi, E., Jakob, C., Kattsov, V., Reason, C., and Rummukaine, M.: Evaluation Of Climate Models, in: Climate Change 2013: the Physical Science Basis, Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, 741–866, https://doi.org/10.1017/CBO9781107415324.020, 2013. a

Foley, J. A., Botta, A., Coe, M. T., and Costa, M. H.: El Niño–Southern Oscillation and the climate, ecosystems and rivers of Amazonia, Global Biogeochem. Cy., 16, 1132, https://doi.org/10.1029/2002GB001872, 2002. a

Foley, J. A., Coe, M. T., Scheffer, M., and Wang, G.: Regime Shifts In The Sahara And Sahel: Interactions Between Ecological And Climatic Systems In Northern Africa, Ecosystems, 6, 524–532, https://doi.org/10.1007/s10021-002-0227-0, 2003. a

Fraedrich, K.: A Suite Of User-Friendly Global Climate Models: Hysteresis Experiments, Eur. Phys. J. Plus, 127, 1–9, https://doi.org/10.1140/epjp/i2012-12053-7, 2012. a

Fraedrich, K., Jansen, H., Kirk, E., Luksch, U., and Lunkeit, F.: The Planet Simulator: Towards A User Friendly Model, Meteorol. Z., 14, 299–304, https://doi.org/10.1127/0941-2948/2005/0043, 2005a. a, b

Fraedrich, K., Jansen, H., Kirk, E., and Lunkeit, F.: The Planet Simulator: Green Planet And Desert World, Meteorol. Z., 14, 305–314, https://doi.org/10.1127/0941-2948/2005/0044, 2005b. a, b, c

Fraedrich, K., Kirk, E., Luksch, U., and Lunkeit, F.: The Portable University Model Of The Atmosphere (PUMA): Storm Track Dynamics And Low-Frequency Variability, Meteorol. Z., 14, 735–745, https://doi.org/10.1127/0941-2948/2005/0074, 2005c. a

Garreaud, R. D., Molina, A., and Farias, M.: Andean Uplift, Ocean Cooling And Atacama Hyperaridity: A Climate Modeling Perspective, Earth Planet. Sc. Lett., 292, 39–50, https://doi.org/10.1016/j.epsl.2010.01.017, 2010. a, b

Giannini, A., Biasutti, M., Held, I. M., and Sobel, A. H.: A Global Perspective On African Climate, Climatic Change, 90, 359–383, https://doi.org/10.1007/s10584-008-9396-y, 2008. a

Gizaw, M. S. and Gan, T. Y.: Impact Of Climate Change And El Niño Episodes On Droughts In Sub-Saharan Africa, Clim. Dynam., 49, 665–682, https://doi.org/10.1007/s00382-016-3366-2, 2017. a

Guilyardi, E.: El Niño-mean state-seasonal cycle interactions in a multi-model ensemble, Clim. Dynam., 26, 329–348, 2006. a

Haywood, A. M. and Valdes, P. J.: Modelling Pliocene warmth: contribution of atmosphere, oceans and cryosphere, Earth Planet. Sc. Lett., 218, 363–377, https://doi.org/10.1016/S0012-821X(03)00685-X, 2004. a, b

Haywood, A. M., Dowsett, H. J., and Dolan, A. M.: Integrating geological archives and climate models for the mid-Pliocene warm period, Nat. Commun., 7, 10646, https://doi.org/10.1038/ncomms10646, 2016. a

Henrot, A.-J., François, L., Brewer, S., and Munhoven, G.: Impacts Of Land Surface Properties And Atmospheric CO2 On The Last Glacial Maximum Climate: A Factor Separation Analysis, Clim. Past, 5, 183–202, https://doi.org/10.5194/cp-5-183-2009, 2009. a

Hirota, M., Holmgren, M., Van Nes, E. H., and Scheffer, M.: Global resilience of tropical forest and savanna to critical transitions, Science, 334, 232–235, https://doi.org/10.1126/science.1210657, 2011. a

Holmgren, M., Scheffer, M., Ezcurra, E., Gutiérrez, J. R., and Mohren, G. M.: El Niño Effects On The Dynamics Of Terrestrial Ecosystems, Trends Ecol. Evol., 16, 89–94, https://doi.org/10.1016/S0169-5347(00)02052-8, 2001. a, b

Holmgren, M., Stapp, P., Dickman, C. R., Gracia, C., Graham, S., Gutiérrez, J. R., Hice, C., Jaksic, F., Kelt, D. A., Letnic, M., Lima, M., López, B., Meserve, P. L., Milstead, W. B., Polis, G. A., Previtali, M. A., Richter, M., Sabaté, S., and Squeo, F. A.: A synthesis of ENSO effects on drylands in Australia, North America and South America, Adv. Geosci., 6, 69–72, https://doi.org/10.5194/adgeo-6-69-2006, 2006. a

Hoskins, B. and Simmons, A.: A multi-layer spectral model and the semi-implicit method, Q. J. Roy. Meteorol. Soc., 101, 637–655, https://doi.org/10.1002/qj.49710142918, 1975. a

Hu, L., Andrews, A. E., Thoning, K. W., Sweeney, C., Miller, J. B., Michalak, A. M., Dlugokencky, E., Tans, P. P., Shiga, Y. P., Mountain, M., Nehrkorn, T., Montzka, S. A., McKain, K., Kofler, J., Trudeau, M., Michel, S. E., Biraud, S. C., Fischer, M. L., Worthy, D. E. J., Vaughn, B. H., White, J. W. C., Yadav, V., Basu, S., and van der Velde, I. R.: Enhanced North American carbon uptake associated with El Niño, Sci. Adv., 5, eaaw0076, https://doi.org/10.1126/sciadv.aaw0076, 2019. a

Hurrell, J. W., Hack, J. J., Shea, D., Caron, J. M., and Rosinski, J.: A New Sea Surface Temperature And Sea Ice Boundary Dataset For The Community Atmosphere Model, J. Climate, 21, 5145–5153, https://doi.org/10.1175/2008JCLI2292.1, 2008. a

Huxman, T. E. and Smith, S. D.: Photosynthesis in an invasive grass and native forb at elevated CO2 during an El Niño year in the Mojave Desert, Oecologia, 128, 193–201, https://doi.org/10.1007/s004420100658, 2001. a

Indeje, M., Semazzi, F. H., and Ogallo, L. J.: ENSO signals in East African rainfall seasons, Int. J. Climatol., 20, 19–46, https://doi.org/10.1002/(SICI)1097-0088(200001)20:1<19::AID-JOC449>3.0.CO;2-0, 2000. a

Ineson, S. and Scaife, A.: The role of the stratosphere in the European climate response to El Niño, Nat. Geoscience, 2, 32–36, https://doi.org/10.1038/ngeo381, 2009. a

Jiang, L., Jiapaer, G., Bao, A., Guo, H., and Ndayisaba, F.: Vegetation dynamics and responses to climate change and human activities in Central Asia, Sci. Total Environ., 599, 967–980, https://doi.org/10.1016/j.scitotenv.2017.05.012, 2017. a

Jiménez-Muñoz, J. C., Mattar, C., Barichivich, J., Santamaría-Artigas, A., Takahashi, K., Malhi, Y., Sobrino, J. A., and van der Schrier, G.: Record-Breaking Warming And Extreme Drought In The Amazon Rainforest During The Course Of El Niño 2015–2016, Scient. Rep., 6, 33130, https://doi.org/10.1038/srep33130, 2016. a

Joly, M. and Voldoire, A.: Influence Of ENSO On The West African Monsoon: Temporal Aspects And Atmospheric Processes, J. Climate, 22, 3193–3210, https://doi.org/10.1175/2008JCLI2450.1, 2009. a

Karmalkar, A. V., Bradley, R. S., and Diaz, H. F.: Climate change in Central America and Mexico: regional climate model validation and climate change projections, Clim. Dynam., 37, 605–629, https://doi.org/10.1007/s00382-011-1099-9, 2011. a

Kleidon, A.: The Climate Sensitivity To Human Appropriation Of Vegetation Productivity And Its Thermodynamic Characterization, Global Planet. Change, 54, 109–127, https://doi.org/10.1016/j.gloplacha.2006.01.016, 2006. a

Kleidon, A.: Non-Equilibrium Thermodynamics, Maximum Entropy Production And Earth-System Evolution, Philos. T. Roy. Soc. Lond. A, 368, 181–196, https://doi.org/10.1098/rsta.2009.0188, 2010. a

Kunz, T., Fraedrich, K., and Lunkeit, F.: Response Of Idealized Baroclinic Wave Life Cycles To Stratospheric Flow Conditions, J. Atmos. Sci., 66, 2288–2302, https://doi.org/10.1175/2009JAS2827.1, 2009. a

Latif, M., Semenov, V. A., and Park, W.: Super El Niños In Response To Global Warming In A Climate Model, Climatic Change, 132, 489–500, https://doi.org/10.1007/s10584-015-1439-6, 2015. a, b

Laurance, W. F., Dell, B., Turton, S. M., Lawes, M. J., Hutley, L. B., McCallum, H., Dale, P., Bird, M., Hardy, G., Prideaux, G., Gawne, B., McMahon, C. R., Yu, R., Hero, J.-M., Schwarzkopf, L., Krockenberger, A., Douglas, M., Silvester, E., Mahony, M., Vella, K., Saikia, U., Wahren, C.-H., Xu, Z., Smith, B., and Cocklin, C.: The 10 Australian Ecosystems Most Vulnerable To Tipping Points, Biol. Conserv., 144, 1472–1480, https://doi.org/10.1016/j.biocon.2011.01.016, 2011. a

Laycock, W. A.: Stable states and thresholds of range condition on North American rangelands: a viewpoint, Rangeland Ecol. Manage./J. Range Manage. Arch., 44, 427–433, 1991. a

Leitold, V., Morton, D. C., Longo, M., dos Santos, M. N., Keller, M., and Scaranello, M.: El Niño Drought Increased Canopy Turnover In Amazon Forests, New Phytol., 219, 959–971, https://doi.org/10.1111/nph.15110, 2018. a

Lenton, T. M.: Early Warning Of Climate Tipping Points, Nat. Clim. Change, 1, 201–209, https://doi.org/10.1038/nclimate1143, 2011. a

Lenton, T. M. and Ciscar, J. C.: Integrating Tipping Points Into Climate Impact Assessments, Climatic Change, 117, 585–597, https://doi.org/10.1007/s10584-012-0572-8, 2013. a

Lenton, T. M. and Williams, H. T.: On The Origin Of Planetary-Scale Tipping Points, Trends Ecol. Evol., 28, 380–382, https://doi.org/10.1016/j.tree.2013.06.001, 2013. a

Lenton, T. M., Held, H., Kriegler, E., Hall, J. W., Lucht, W., Rahmstorf, S., and Schellnhuber, H. J.: Tipping Elements In The Earth's Climate System, P. Natl. Acad. Sci. USA, 105, 1786–1793, https://doi.org/10.1073/pnas.0705414105, 2008. a, b, c, d

Lin, J.-L.: The double-ITCZ problem in IPCC AR4 coupled GCMs: Ocean-atmosphere feedback analysis, J. Climate, 20, 4497–4525, 2007. a

Lin, P., Paynter, D., Ming, Y., and Ramaswamy, V.: Changes of the tropical tropopause layer under global warming, J. Climate, 30, 1245–1258, 2017. a

Linsenmeier, M., Pascale, S., and Lucarini, V.: Climate Of Earth-Like Planets With High Obliquity And Eccentric Orbits: Implications For Habitability Conditions, Planet. Space Sci., 105, 43–59, https://doi.org/10.1016/j.pss.2014.11.003, 2015. a, b, c

Lucarini, V., Fraedrich, K., and Lunkeit, F.: Thermodynamic Analysis Of Snowball Earth Hysteresis Experiment: Efficiency, Entropy Production And Irreversibility, Q. J. Roy. Meteorol. Soc., 136, 2–11, https://doi.org/10.1002/qj.543, 2010. a

Lucarini, V., Pascale, S., Boschi, R., Kirk, E., and Iro, N.: Habitability And Multistability In Earth-Like Planets, Astron. Nachr., 334, 576–588, https://doi.org/10.1002/asna.201311903, 2013. a, b

Lucarini, V., Ragone, F., and Lunkeit, F.: Predicting Climate Change Using Response Theory: Global Averages And Spatial Patterns, J. Stat. Phys., 166, 1036–1064, https://doi.org/10.1007/s10955-016-1506-z, 2017. a, b

Lunkeit, F., Borth, H., Böttinger, M., Fraedrich, K., Jansen, H., Kirk, E., Kleidon, A., Luksch, U., Paiewonsky, P., Schubert, S., Sielmann, S., and Wan, H.: Planet Simulator Reference Manual Version 16, Meteorologisches Institut, Universität Hamburg, Hamburg, 2011. a, b

Lyon, B.: Seasonal Drought In The Greater Horn Of Africa And Its Recent Increase During The March–May Long Rains, J. Climate, 27, 7953–7975, https://doi.org/10.1175/JCLI-D-13-00459.1, 2014. a

Malhi, Y., Baldocchi, D., and Jarvis, P.: The Carbon Balance Of Tropical, Temperate And Boreal Forests, Plant Cell Environ., 22, 715–740, https://doi.org/10.1046/j.1365-3040.1999.00453.x, 1999. a

Malhi, Y., Aragão, L. E., Galbraith, D., Huntingford, C., Fisher, R., Zelazowski, P., Sitch, S., McSweeney, C., and Meir, P.: Exploring The Likelihood And Mechanism Of A Climate-Change-Induced Dieback Of The Amazon Rainforest, P. Natl. Acad. Sci. USA, 106, 20610–20615, https://doi.org/10.1073/pnas.0804619106, 2009. a, b, c, d, e, f, g

McAlpine, C., Syktus, J., Ryan, J., Deo, R., McKeon, G., McGowan, H., and Phinn, S.: A Continent Under Stress: Interactions, Feedbacks And Risks Associated With Impact Of Modified Land Cover On Australia's Climate, Global Change Biol., 15, 2206–2223, https://doi.org/10.1111/j.1365-2486.2009.01939.x, 2009. a

McPhaden, M. J., Zebiak, S. E., and Glantz, M. H.: ENSO As An Integrating Concept In Earth Science, Science, 314, 1740–1745, https://doi.org/10.1126/science.1132588, 2006. a, b

Mercado-Bettin, D., Salazar, J., and Villegas, J.: Long-term water balance partitioning explained by physical and ecological characteristics in world river basins, Ecohydrology, 12, e2072, https://doi.org/10.1002/eco.2072, 2019. a

Molina, R. D., Salazar, J. F., Martínez, J. A., Villegas, J. C., and Arias, P. A.: Forest-induced exponential growth of precipitation along climatological wind streamlines over the Amazon, J. Geophys. Res.-Atmos., 124, 2589–2599, https://doi.org/10.1029/2018JD029534, 2019. a

Morice, C. P., Kennedy, J. J., Rayner, N. A., and Jones, P. D.: Quantifying Uncertainties In Global And Regional Temperature Change Using An Ensemble Of Observational Estimates: The HadCRUT4 Data Set, J. Geophys. Res.-Atmos., 117, 1–22, https://doi.org/10.1029/2011jd017187, 2012. a

Myers, N., Mittermeier, R. A., Mittermeier, C. G., Da Fonseca, G. A., and Kent, J.: Biodiversity hotspots for conservation priorities, Nature, 403, 853–858, https://doi.org/10.1038/35002501, 2000. a

Nicholson, S. E.: The West African Sahel: A Review Of Recent Studies On The Rainfall Regime And Its Interannual Variability, ISRN Meteorol., 2013, 453521, https://doi.org/10.1155/2013/453521, 2013. a

Nicholson, S. E. and Kim, J.: The Relationship Of The El Niño–Southern Oscillation To African Rainfall, Int. J. Climatol., 17, 117–135, https://doi.org/10.1002/(SICI)1097-0088(199702)17:2<117::AID-JOC84>3.0.CO;2-O, 1997. a

Nicholson, S. E., Funk, C., and Fink, A. H.: Rainfall Over The African Continent From The 19th Through The 21st Century, Global Planet. Change, 165, 114–127, https://doi.org/10.1016/j.gloplacha.2017.12.014, 2017. a

Oort, A. H. and Yienger, J. J.: Observed Interannual Variability In The Hadley Circulation And Its Connection To ENSO, J. Climate, 9, 2751–2767, https://doi.org/10.1175/1520-0442(1996)009<2751:OIVITH>2.0.CO;2, 1996. a

Orszag, S. A.: Transform Method For The Calculation Of Vector-Coupled Sums: Application To The Spectral Form Of The Vorticity Equation, J. Atmos. Sci., 27, 890–895, https://doi.org/10.1175/1520-0469(1970)027<0890:TMFTCO>2.0.CO;2, 1970. a

Oyama, M. D. and Nobre, C. A.: A New Climate-Vegetation Equilibrium State for Tropical South America, Geophys. Res. Lett., 30, 2199, https://doi.org/10.1029/2003gl018600, 2003. a

Pan, Y., Birdsey, R. A., Fang, J., Houghton, R., Kauppi, P. E., Kurz, W. A., Phillips, O. L., Shvidenko, A., Lewis, S. L., Canadell, J. G., Ciais, P., Jackson, R. B., Pacala, S. W., McGuire, A. D., Piao, S., Rautiainen, A., Sitch, S., and Hayes, D.: A Large And Persistent Carbon Sink In The World's Forests, Science, 333, 988–993, https://doi.org/10.1126/science.1201609, 2011. a, b

Phillips, O. L., Aragão, L. E., Lewis, S. L., Fisher, J. B., Lloyd, J., López-González, G., Malhi, Y., Monteagudo, A., Peacock, J., Quesada, C. A., van der Heijden, G., Almeida, S., Amaral, I., Arroyo, L., Aymard, G., Baker, T. R., Bánki, O., Blanc, L., Bonal, D., Brando, P., Chave, J., Alves de Oliveira, Á. C., Cardozo, N. D., Czimczik, C. I., Feldpausch, T. R., Freitas, M. A., Gloor, E., Higuchi, N., Jiménez, E., Lloyd, G., Meir, P., Mendoza, C., Morel, A., Neill, D. A., Nepstad, D., Patiño, S., Peñuela, M. C., Prieto, A., Ramírez, F., Schwarz, M., Silva, J., Silveira, M., Thomas, A. S., ter Steege, H., Stropp, J., Vásquez, R., Zelazowski, P., Alvarez Dávila, E., Andelman, S., Andrade, A., Chao, K.-J., Erwin, T., Di Fiore, A., Honorio, E. C., Keeling, H., Killeen, T. J., Laurance, W. F., Peña Cruz, A., Pitman, N. C. A., Núñez Vargas, P., Ramírez-Angulo, H., Rudas, A., Salamão, R., Silva, N., Terborgh, J., and Torres-Lezama, A.: Drought Sensitivity Of The Amazon Rainforest, Science, 323, 1344–1347, https://doi.org/10.1126/science.1164033, 2009. a

Pithan, F., Shepherd, T. G., Zappa, G., and Sandu, I.: Climate Model Biases In Jet Streams, Blocking And Storm Tracks Resulting From Missing Orographic Drag, Geophys. Res. Lett., 43, 7231–7240, https://doi.org/10.1002/2016gl069551, 2016. a

Pomposi, C., Giannini, A., Kushnir, Y., and Lee, D. E.: Understanding Pacific Ocean Influence On Interannual Precipitation Variability In The Sahel, Geophys. Res. Lett., 43, 9234–9242, https://doi.org/10.1002/2016GL069980, 2016. a

Poveda, G., Mesa, O. J., and Waylen, P. R.: Nonlinear forecasting of river flows in Colombia based upon ENSO and its associated economic value for hydropower generation, in: Climate and water, Springer, 351–371, https://doi.org/10.1007/978-94-015-1250-3_15, 2003. a

Poveda, G., Waylen, P. R., and Pulwarty, R. S.: Annual and inter-annual variability of the present climate in northern South America and southern Mesoamerica, Palaeogeogr. Palaeocl. Palaeoecol., 234, 3–27, https://doi.org/10.1016/j.palaeo.2005.10.031, 2006. a, b

Prein, A. F. and Gobiet, A.: Impacts Of Uncertainties In European Gridded Precipitation Observations On Regional Climate Analysis, Int. J. Climatol., 37, 305–327, https://doi.org/10.1002/joc.4706, 2017. a

Ratnam, J., Behera, S., Masumoto, Y., and Yamagata, T.: Remote Effects Of El Niño And Modoki Events On The Austral Summer Precipitation Of Southern Africa, J. Climate, 27, 3802–3815, https://doi.org/10.1175/JCLI-D-13-00431.1, 2014. a

Ropelewski, C. F. and Halpert, M. S.: North American precipitation and temperature patterns associated with the El Niño/Southern Oscillation (ENSO), Mon. Weather Rev., 114, 2352–2362, https://doi.org/10.1175/1520-0493(1986)114<2352:NAPATP>2.0.CO;2, 1986. a

Rouault, M. and Richard, Y.: Intensity And Spatial Extent Of Droughts In Southern Africa, Geophys. Res. Lett., 32, L15702, https://doi.org/10.1029/2005GL022436, 2005. a

Salazar, J. F., Villegas, J. C., Rendón, A. M., Rodríguez, E., Hoyos, I., Mercado-Bettín, D., and Poveda, G.: Scaling properties reveal regulation of river flows in the Amazon through a forest reservoir, Hydrol. Earth Syst. Sci., 22, 1735–1748, https://doi.org/10.5194/hess-22-1735-2018, 2018. a, b

Salzmann, U., Williams, M., Haywood, A. M., Johnson, A. L., Kender, S., and Zalasiewicz, J.: Climate and environment of a Pliocene warm world, Palaeogeogr. Palaeocl. Palaeoecol., 309, 1–8, https://doi.org/10.1016/j.palaeo.2011.05.044, 2011. a

Santer, B. D., Solomon, S., Pallotta, G., Mears, C., Po-Chedley, S., Fu, Q., Wentz, F., Zou, C.-Z., Painter, J., Cvijanovic, I., and Bonfils, C.: Comparing tropospheric warming in climate models and satellite data, J. Climate, 30, 373–392, 2017. a

Sausen, R., Schubert, S., and Dümenil, L.: A model of river runoff for use in coupled atmosphere-ocean models, J. Hydrol., 155, 337–352, https://doi.org/10.1016/0022-1694(94)90177-5, 1994. a

Scheffer, M., Carpenter, S., Foley, J. A., Folke, C., and Walker, B.: Catastrophic Shifts In Ecosystems, Nature, 413, 591–596, https://doi.org/10.1038/35098000, 2001. a

Scheffer, M., Bascompte, J., Brock, W. A., Brovkin, V., Carpenter, S. R., Dakos, V., Held, H., Van Nes, E. H., Rietkerk, M., and Sugihara, G.: Early-Warning Signals For Critical Transitions, Nature, 461, 53–59, https://doi.org/10.1038/nature08227, 2009. a

Schmittner, A., Silva, T. A., Fraedrich, K., Kirk, E., and Lunkeit, F.: Effects Of Mountains And Ice Sheets On Global Ocean Circulation, J. Climate, 24, 2814–2829, https://doi.org/10.1175/2010JCLI3982.1, 2011. a

Seager, R. and Henderson, N.: Diagnostic computation of moisture budgets in the ERA-Interim reanalysis with reference to analysis of CMIP-archived atmospheric model data, J. Climate, 26, 7876–7901, https://doi.org/10.1175/JCLI-D-13-00018.1, 2013. a

Semtner Jr., A. J.: A model for the thermodynamic growth of sea ice in numerical investigations of climate, J. Phys. Oceanog., 6, 379–389, https://doi.org/10.1175/1520-0485(1976)006<0379:AMFTTG>2.0.CO;2, 1976. a

Sheen, K., Smith, D., Dunstone, N., Eade, R., Rowell, D., and Vellinga, M.: Skilful Prediction Of Sahel Summer Rainfall On Inter-Annual And Multi-Year Timescales, Nat. Commun., 8, 14966, https://doi.org/10.1038/ncomms14966, 2017. a

Spiegl, T., Paeth, H., and Frimmel, H.: Evaluating Key Parameters For The Initiation Of A Neoproterozoic Snowball Earth With A Single Earth System Model Of Intermediate Complexity, Earth Planet. Sc. Lett., 415, 100–110, https://doi.org/10.1016/j.epsl.2015.01.035, 2015. a

Staver, A. C., Archibald, S., and Levin, S. A.: The Global Extent And Determinants Of Savanna And Forest As Alternative Biome States, Science, 334, 230–232, https://doi.org/10.1126/science.1210465, 2011. a, b

Stevens, B., Giorgetta, M., Esch, M., Mauritsen, T., Crueger, T., Rast, S., Salzmann, M., Schmidt, H., Bader, J., Block, K., Brokopf, R., Fast, I., Kinne, S., Kornblueh, L., Lohmann, U., Pincus, R., Reichler, T., and Roeckner, E.: Atmospheric Component Of The MPI-M Earth System Model: ECHAM6, J. Adv. Model. Earth Syst., 5, 146–172, https://doi.org/10.1002/jame.20015, 2013. a

Taylor, K. E., Williamson, D., and Zwiers, F.: The sea surface temperature and sea-ice concentration boundary conditions for AMIP II simulations, Tech. Rep. 60, Program for Climate Model Diagnosis and Intercomparison, Lawrence Livermore National Laboratory, University of California, California, 2000. a

Timmermann, A., Oberhuber, J., Bacher, A., Esch, M., Latif, M., and Roeckner, E.: Increased El Niño Frequency In A Climate Model Forced By Future Greenhouse Warming, Nature, 398, 694–697, https://doi.org/10.1038/19505, 1999. a

Tippett, M. K. and Giannini, A.: Potentially Predictable Components Of African Summer Rainfall In An SST-forced GCM Simulation, J. Climate, 19, 3133–3144, https://doi.org/10.1175/JCLI3779.1, 2006. a

Trenberth, K. E., Fasullo, J. T., and Mackaro, J.: Atmospheric Moisture Transports From Ocean To Land And Global Energy Flows In Reanalyses, J. Climate, 24, 4907–4924, https://doi.org/10.1175/2011jcli4171.1, 2011. a, b, c

Wang, G. and Hendon, H. H.: Sensitivity Of Australian Rainfall To Inter-El Niño Variations, J. Climate, 20, 4211–4226, https://doi.org/10.1175/JCLI4228.1, 2007. a

Wara, M. W., Ravelo, A. C., and Delaney, M. L.: Permanent El Niño-like conditions during the Pliocene warm period, Science, 309, 758–761, https://doi.org/10.1126/science.1112596, 2005. a, b, c

Warman, L. and Moles, A. T.: Alternative stable states in Australia's Wet Tropics: a theoretical framework for the field data and a field-case for the theory, Landsc. Ecol., 24, 1–13, https://doi.org/10.1007/s10980-008-9285-9, 2009. a

Watson, J. E., Evans, T., Venter, O., Williams, B., Tulloch, A., Stewart, C., Thompson, I., Ray, J. C., Murray, K., Salazar, A., McAlpine, C., Potapov, P., Walston, J., Robinson, J. G., Painter, M., Wilkie, D., Filardi, C., Laurance, W. F., Houghton, R. A., Maxwell, S., Grantham, H., Samper, C., Wang, S., Laestadius, L., Runting, R. K., Silva-Chávez, G. A., Ervin, J., and Lindenmayer, D.: The Exceptional Value Of Intact Forest Ecosystems, Nat. Ecol. Evol., 2, 599–610, https://doi.org/10.1038/s41559-018-0490-x, 2018. a

Weng, W., Luedeke, M. K., Zemp, D. C., Lakes, T., and Kropp, J. P.: Aerial and surface rivers: downwind impacts on water availability from land use changes in Amazonia, Hydrol. Earth Syst. Sci., 22, 911–927, https://doi.org/10.5194/hess-22-911-2018, 2018.  a

Wu, Z., Boke-Olén, N., Fensholt, R., Ardö, J., Eklundh, L., and Lehsten, V.: Effect of climate dataset selection on simulations of terrestrial GPP: Highest uncertainty for tropical regions, PloS One, 13, e0199383, https://doi.org/10.1371/journal.pone.0199383, 2018. a

Yang, K., Koike, T., Fujii, H., Tamagawa, K., and Hirose, N.: Improvement of Surface Flux Parametrizations With A Turbulence-Related Length, Q. J. Roy. Meteorol. Soc., 128, 2073–2087, https://doi.org/10.1256/003590002320603548, 2002. a

Yeh, S.-W., Cai, W., Min, S.-K., McPhaden, M. J., Dommenget, D., Dewitte, B., Collins, M., Ashok, K., An, S.-I., Yim, B.-Y., and Kug, J.-S.: ENSO Atmospheric Teleconnections And Their Response To Greenhouse Gas Forcing, Rev. Geophys., 56, 185–206, https://doi.org/10.1002/2017RG000568, 2018. a

Zebiak, S. E., Orlove, B., Muñoz, Á. G., Vaughan, C., Hansen, J., Troy, T., Thomson, M. C., Lustig, A., and Garvin, S.: Investigating El Niño–Southern Oscillation And Society Relationships, Wiley Interdisciplin. Rev.: Clim. Change, 6, 17–34, https://doi.org/10.7916/D84X577H, 2015. a

Zemp, D. C., Schleussner, C.-F., Barbosa, H. M., Hirota, M., Montade, V., Sampaio, G., Staal, A., Wang-Erlandsson, L., and Rammig, A.: Self-amplified Amazon Forest Loss Due To Vegetation-Atmosphere Feedbacks, Nat. Commun., 8, 14681, https://doi.org/10.1038/ncomms14681, 2017. a

Zeng, N. and Neelin, J. D.: The Role Of Vegetation-Climate Interaction And Interannual Variability In Shaping The African Savanna, J. Climate, 13, 2665–2670, https://doi.org/10.1175/1520-0442(2000)013<2665:TROVCI>2.0.CO;2, 2000. a

Zeng, N. and Yoon, J.: Expansion Of The World's Deserts Due To Vegetation-Albedo Feedback Under Global Warming, Geophys. Res. Lett., 36, L17401, https://doi.org/10.1029/2009GL039699, 2009. a

Zhang, Y. G., Pagani, M., and Liu, Z.: A 12-million-year temperature history of the tropical Pacific Ocean, Science, 344, 84–87, https://doi.org/10.1126/science.1246172, 2014. a

Zhao, M., Heinsch, F. A., Nemani, R. R., and Running, S. W.: ImprovementS Of The MODIS Terrestrial Gross And Net Primary Production Global Data Set, Remote Sens. Environ., 95, 164–176, https://doi.org/10.1016/j.rse.2004.12.011, 2005. a

Zhou, C. and Wang, K.: Evaluation Of Surface Fluxes In ERA-Interim Using Flux Tower Data, J. Climate, 29, 1573–1582, https://doi.org/10.1175/JCLI-D-15-0523.1, 2016. a