# Comparison of the green-to-desert Sahara transitions between the Holocene and the last interglacial

Huan Li, Hans Renssen, Didier M. Roche

► To cite this version:

Huan Li, Hans Renssen, Didier M. Roche. Comparison of the green-to-desert Sahara transitions between the Holocene and the last interglacial. Climate of the Past, 2022, 18, pp.2303-2319. 10.5194/cp-18-2303-2022. insu-03993101

HAL Id: insu-03993101

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

Submitted on 16 Feb 2023

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>

# Comparison of the green-to-desert Sahara transitions between the Holocene and the last interglacial

**Huan Li**<fnref n="1" />,<fnref n="2" />, **Hans Renssen**<fnref n="2" />,<fnref n="3" />, and **Didier M. Roche**<fnref n="2" />,<fnref n="4" />

**Correspondence:** Hans Renssen (hans.renssen@usn.no)

Received: 8 March 2022 – Discussion started: 9 March 2022

Revised: 16 September 2022 – Accepted: 17 September 2022 – Published: 18 October 2022

**Abstract.** The desertification and vegetation feedbacks of the “green Sahara” during the last interglacial (LIG) and the Holocene have been investigated by many studies. Yet the abruptness of climate and vegetation changes and their interactions are still under discussion. In this study, we apply an Earth system model of intermediate complexity (iLOVECLIM) in combination with two dynamical vegetation models (VECODE and LPJ-GUESS) to simulate climate–vegetation changes during the Holocene and the LIG to compare the patterns of North African vegetation evolutions and mechanisms of their feedbacks during these two interglacials. Our results confirmed the existence of the green Sahara during the early LIG, which is as an analogue to the green Sahara during the Holocene. During both interglacials, an overall consistent transition from vegetated Sahara to desert is shown in our results, but the amplitudes of these transitions vary. These simulated Sahara vegetation transitions are nearly linearly related to the summer insolation declines at 20° N, resulting in faster declines of vegetation cover during the LIG than in the Holocene. The decline of vegetation cover peaks at 25% kyr⁻<fnref n="1" /> at around 122 ka, while during the Holocene the steepest vegetation cover decline is 10% kyr⁻<fnref n="1" /> at around 6 ka. Our results suggest net positive vegetation feedbacks to climate during the two interglacials. During the early LIG and Holocene, vegetation strengthens precipitation by a factor of 2 to 3 through the vegetation–albedo feedback when the vegetation cover is greater than

60%. Vegetation cover decreases with declines of the incoming moisture transport by the atmosphere due to the reduced summer insolation at 20° N, weakening the summer monsoon during both interglacials. This desertification is accelerated when the positive vegetation–albedo feedback cannot offset the reduction of precipitation due to a weaker summer monsoon. The impacts of this positive vegetation feedback on precipitation decrease with decreased vegetation cover, during which the impacts of negative vegetation–evaporation feedbacks increase, accelerating the loss of soil moisture and vegetation cover. Overall, the net positive vegetation feedback is strong during the early phases of both interglacials, but the vegetation transition is more abrupt during the LIG than during the Holocene due to the more rapid changes in summer insolation during the LIG. The main difference between the two interglacials is the rate of precipitation change, which is relatively gradual during the Holocene, leading to a more gradual vegetation transition in comparison to the LIG.

## 1 Introduction

From about 14.8 to 5.5 ka (thousands of years before present), relatively high levels of summer insolation in the Northern Hemisphere (NH) led to an enhanced African monsoon, resulting in a humid and green Sahara (Kutzbach and Street-Perrot, 1985; Jolly et al., 1998; Prentice and Jolly,

<footnote n="1">¹School of Geographic Science, Nantong University, Tongjingdadao 999, Nantong 226007, China</footnote>

<footnote n="2">²Vrije Universiteit Amsterdam, Faculty of Science, Cluster Earth and Climate, de Boelelaan 1085, Amsterdam, the Netherlands</footnote>

<footnote n="3">³Department of Natural Sciences and Environmental Health, University of South-Eastern Norway, Gullbringvegen 36, 3800 Bø i Telemark, Norway</footnote>

<footnote n="4">⁴Laboratoire des Sciences du Climat et de l'Environnement, LSCE/IPSL, CEA-INSU-UVSQ-CNRS, Université Paris-Saclay, Gif-sur-Yvette, France</footnote>
2000; Tuenter et al., 2007; Rachmayani et al., 2016). This period is therefore referred to as the African Humid Period (AHP, Ritchie et al., 1985), featured by shrub- and grass-covered land surface in northern Africa where there is desert today (Hoelzmann et al., 2004; Lézine, 2017). Forced by the decline in summer insolation in the NH, the termination of the AHP occurred at about 5.5 ka, featured by a rapid reduction in precipitation and a transition from the vegetated "green" Sahara to desert (Jolly et al., 1998; deMenocal et al., 2000; Prentice and Jolly, 2000; Levis et al., 2004; Liu et al., 2006, 2007; Renssen et al., 2006; Kröpelin et al., 2008; Notaro et al., 2008; Wang et al., 2008; Lezine et al., 2011; Claussen et al., 2013; Francus et al., 2013; Rachmayani et al., 2015). Several studies indicate considerable spatial variations in the nature of this transition (Shanahan et al., 2015; Tierney et al., 2017; Dallmeyer et al., 2020). Moreover, the transition rate varies among data and models, and thus the character of the transition remains controversial (e.g., Pausata et al., 2016; Tierney et al., 2017). DeMenocal et al. (2000) inferred from marine sedimentary data that dramatic desertification in North Africa occurred abruptly within centuries, while several other studies have suggested slower desertification in response to the gradual decline in summer insolation (Wang et al., 2005; Renssen et al., 2006; Kröpelin et al., 2008; Francus et al., 2013; Shanahan et al., 2015). One reason for these contrasting rates of desertification could be related to variations in the strength of vegetation feedbacks that amplify the regional climate response to gradually declining orbital forcing. This is supported by a recent modeling study that showed that the abruptness of the green-desert transition is highly dependent on the model setup (Hopcroft and Valdes, 2021).

Understanding the termination of the AHP thus sheds light on the nonlinear response of the African monsoon to orbital forcing and the role of vegetation feedbacks. Several studies (e.g., Tuenter et al., 2007; Rachmayani et al., 2015, 2016) have investigated to what extent the AHP can be attributed to the long-term changes in the orbital forcing and the feedback processes, respectively. Tuenter et al. (2007) and Rachmayani et al. (2016) demonstrated that the long-term varying precession is the pacemaker of the African monsoon and thereby of precipitation over northern Africa. However, a factor analysis by Rachmayani et al. (2015) found an amplifying effect of vegetation feedback on mid-Holocene precipitation over West Africa. Moreover, the importance of surface albedo processes due to vegetation changes was supported by several studies (e.g., Levis et al., 2004; Liu et al., 2006, 2007; Notaro et al., 2008; Wang et al., 2008; Fischer and Jungclaus, 2010; Vamborg et al., 2011; Rachmayani et al., 2015). Still, the mechanism of the termination of the AHP remains controversial, and only a few studies have focused on the abruptness of this termination. We consider the desertification-vegetation transition to be abrupt when the rate of the transition is faster than the gradual decline in orbital forcing and discuss the role of orbital forcing and vegetation

feedback in the termination of the AHP through comparisons with two sets of conceptual models (Liu et al., 2006; Claussen et al., 1999, 2017). Liu et al. (2006) proposed two possibilities: either the abrupt vegetation transition in North Africa was caused by a strong vegetation feedback combined with abrupt precipitation reduction (so-called unstable collapse or UC) or by a weak vegetation feedback combined with a strong low-frequency variability of gradual decline in precipitation (so-called stable collapse or SC). In the model experiments of Liu et al. (2006), the vegetation transition in northern Africa was caused by SC. In another modeling study (Claussen et al., 1999), support for the first possibility (UC) was found, with vegetation collapse resulting from strong positive vegetation feedback combined with abrupt precipitation changes. In order to understand the rate of vegetation change during the termination of the AHP, the magnitude of vegetation feedbacks is thus the key factor under the gradual change in summer insolation in NH during the Holocene.

The magnitudes of climate-vegetation feedbacks have been investigated by several studies. For example, a positive vegetation-albedo feedback was found to be an important factor amplifying changes in North African precipitation during the Holocene (Charney, 1975; Claussen and Gayler, 1997; Texier et al., 1997; Knorr and Schnitzler, 2006). Compared to desert, vegetated surface is characterized by a lower albedo (e.g., 0.20 vs. 0.33 in our iLOVECLIM model), which enhances the absorption of shortwave radiative energy and therefore leads to a warmer surface. Model studies suggest that vegetated surfaces in North Africa are warmer than desert, inducing larger land-sea temperature contrasts and a stronger African summer monsoon, which promotes vegetation development because of the enhanced precipitation. Therefore, the key factors adjusting the magnitude of this vegetation-albedo feedback are vegetation cover and differences in surface albedo between vegetated and bare soil surfaces. However, a vegetated surface affects climate not only through a lower surface albedo compared to bare soil, but also through its impacts on evaporation and transpiration. Compared to desert, on the one hand, vegetated surfaces promote evapotranspiration, leading to evaporative cooling of the surface, resulting in lower temperature and a reduced tendency for convective precipitation. On the other hand, an increased latent heat flux (evapotranspiration) leads to more humidity in the air, which can lead to enhanced precipitation through local recirculation. Whether the net vegetation feedback through changes in evaporation and transpiration is positive or negative to precipitation is related to several factors, e.g., soil moisture availability, surface albedo, the percentage of organic matter, and the amount of evaporation from different types of land surfaces. For example, a dark wet soil (Levis et al., 2004) and the inclusion of canopy transpiration (Rachmayani et al., 2015) positively impact precipitation, while the feedback negatively affects precipitation when the evaporation from bare soil is stronger than transpiration from vegetated surfaces (Liu et al., 2007; Wang et al., 2008; Notaro
et al., 2008). The net effect on the surface temperature of vegetation depends on the balance between the impact of the relatively low albedo, promoting shortwave radiation absorption and warming, and the effect of an enhanced latent heat flux, leading to evaporative cooling.

The last interglacial (LIG) is an interesting test period with insolation changes between 129 and 120 ka that were stronger than in the Holocene (Fig. 1; Berger, 1978) but with similar greenhouse gas (GHG) concentrations and continental configurations (Otto-Bliesner et al., 2017). In addition, proxy-based evidence indicates that the Sahara was also green and humid in the early LIG and experienced a transition from green to desert during the LIG. This evidence for a humid period during the LIG is based on marine sediment records (Rohling et al., 2002; Osborne et al., 2008), speleothem records (Fleitmann et al., 2011), and terrestrial sediment data (Ehrmann et al., 2017; Scussolini et al., 2019). The impacts of vegetation feedback on climate during the LIG were also studied by several modeling studies (e.g., Gröger et al., 2007; Schurgers et al., 2007). For example, Schurgers et al. (2007) suggested as much as 20% local albedo changes in northwestern Africa during the early LIG and more than 10% increases in global latent heat loss compared to the present day. Gröger et al. (2007) simulated a more than doubling of precipitation in the Sahara (North Africa) due to vegetation feedbacks during the LIG. In a previous study, we performed a series of model experiments assessing the vegetation impacts on the LIG climate (Li et al., 2020). Our global results generally compared favorably to proxy-based reconstructions, including clear warming over North Africa (Turney and Jones, 2010; McKay et al., 2011; Li et al., 2020). In addition, we found consistent magnitudes of vegetation feedback with previous modeling studies, including about twice as large amounts of precipitation (~60 cm yr⁻¹) and higher surface temperatures by about 2.5° in North Africa in response to dynamical vegetation evolution during the early LIG compared to simulations with fixed pre-industrial vegetation.

The controversial abruptness of vegetation transitions and vegetation feedback during the AHP could be related to the gradual changes in external orbital forcing during the Holocene, which may not be strong enough to counterbalance the effects of internal vegetation-climate interactions. During the LIG, conditions were different from the Holocene, and therefore a comparison of North African vegetation transition between the LIG and the Holocene could help us disentangle the relative contribution of the external and internal processes to North African desertification.

In this paper, we first compare the LIG and Holocene transient simulations from an intermediate-complexity model (iLOVECLIM), including different vegetation components with different complexities (VECODE and LPJ-GUESS). The LIG results cover 127 to 116 ka and are the same as reported in our previous study (Li et al., 2020). The Holocene simulations cover 8.5 to 1 ka. The simulations of the two

interglacials with the same climate model offer the opportunity to analyze the impacts of the external forcings without model dependence. In addition, the use of two vegetation models increases our insight into the model dependency of the simulated vegetation states. Thus, by performing transient climate model experiments on two different interglacials and by applying two different vegetation models of different complexity, we are able to provide a thorough evaluation of the abruptness and the most important factors contributing to the green-desert transitions. In our simulations, VECODE was synchronously coupled, while LPJ-GUESS was asynchronously coupled to iLOVECLIM. We then assess the North African vegetation feedback by switching on and off interactive dynamic vegetation in each coupled version. Finally, we compare the climate and vegetation simulations with the same standardized spatial scale between the two interglacial periods. We analyze the vegetation transitions during both interglacials with the aim of addressing the following questions.

1. Are the evolutions of North African (NA) vegetation during the two interglacials simulated by iLOVECLIM with different vegetation components consistent?

2. Are the evolutions of NA vegetation different during the Holocene and LIG, and what are the dominating factors?

3. How do the NA vegetation geophysical feedbacks affect the climate during the Holocene and LIG?

# 2 Methods

## 2.1 Model descriptions

We applied the climate model, iLOVECLIM version 1.0, which is an updated version of LOVECLIM 1.2 (Goosse et al., 2010; Roche et al., 2014). This intermediate-complexity model includes the main climatic system components and has been successfully applied to analyze climate-vegetation interactions during the Holocene AHP (Renssen et al., 2003, 2006) and the LIG (e.g., Li et al., 2020). In the present study, we make use of three dynamical components in iLOVECLIM representing the atmosphere (ECBilt), the oceans (CLIO), and land-based vegetation (either VECODE or LPJ-GUESS). ECBilt is a three-level quasi-geostrophic atmospheric model at T21 resolution, corresponding to ~625 km at the Equator (Opsteegh et al., 1998). ECBilt calculates the climatic forcing for the applied vegetation component (either VECODE or LPJ-GUESS). Cloud cover is prescribed according to present-day climatology (ISCCP D2 dataset, Rossow et al., 1996). In ECBilt, the hydrological cycle is closed by applying a bucket for soil moisture (Goosse et al., 2010). Every bucket is connected and drains to a nearby oceanic grid cell to calculate river runoff. The surface albedo is a function of the fraction of the grid box covered by
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Evolution of 20° N July insolation, North African precipitation, surface temperature, and values of evaporation minus precipitation during the LIG (a, b, c, d) and Holocene (e, f, g, h). Our target area is delimited by 10° W to 35° E and 15 to 30° N. The blue lines and red lines are simulated results from transient experiments with dynamical vegetation and fixed pre-industrial vegetation, respectively. The green dots with error bars are simulated results from experiments with asynchronously coupled “dynamical” vegetation.</figcaption></figure>

the different surface types, i.e., ocean, sea ice, trees, desert, and grass (Goosse et al., 2010). CLIO consists of a three-dimensional, free surface ocean general circulation model coupled to a dynamic-thermodynamic sea ice model (Goosse et al., 2010). This oceanic component has a horizontal resolution of 3° latitude by 3° longitude and 20 unevenly spaced vertical layers in the ocean.

As in Li et al. (2020), we apply three model configurations here: (1) ECBilt-CLIO fully coupled to VECODE (ECBilt-CLIO-VECODE), (2) ECBilt-CLIO asynchronously coupled to LPJ-GUESS (ECBilt-CLIO_LPJ-GUESS), and (3) ECBilt-CLIO with prescribed pre-industrial vegetation (ECBilt-CLIO-LUH2). The difference among these three model configurations is their vegetation component and thereby the surface albedo and soil moisture.

In ECBilt-CLIO-VECODE, with climate forcing from ECBilt, VECODE calculates cover for two plant functional

types (PFTs: grasses and trees) and bare soil as a dummy type at an annual time step (Brovkin et al., 1997). The climate forcing for VECODE includes annual mean temperature, precipitation, and GDD0 (growing degree days above 0°). We note that using annual mean climate input implies that the vegetation in VECODE is not directly affected by seasonal variations in precipitation and temperature, although some seasonal information is provided through the GDD0. This simplification could potentially lead to an underestimation of the impact of orbital forcing on vegetation in North Africa.

Similar to ECBilt-CLIO-VECODE, LPJ-GUESS also calculates PFT cover in ECBilt-CLIO_LPJ-GUESS based on climate input from ECBilt, but LPJ-GUESS requires more specific climate input and calculates vegetation dynamics in more detail than VECODE. We apply version 3.1 of LPJ-GUESS (Smith et al., 2001, 2014). The climate in-
put required by LPJ-GUESS includes monthly mean surface temperature, precipitation, and cloud cover. The model employs biophysical and physiological process parameterizations identical to the equilibrium model BIOME3 (Haxeltine and Prentice, 1996). Vegetation dynamics in LPJ-GUESS result from growth and competition for light, space, and soil resources among 11 PFTs (Smith et al., 2001). These 11 PFTs are boreal needle-leaved evergreen trees, boreal needle-leaved evergreen shade-intolerant trees, boreal needle-leaved summer-green trees, temperate broadleaved summer-green trees, boreal-temperate broadleaved summer-green trees, temperate broadleaved evergreen trees, tropical broadleaved evergreen shade-intolerant trees, temperate broadleaved rain-green trees, C₃ grass, and C₄ grass. Physiological processes are simulated with a daily time step. The net primary production (NPP) is accrued at the end of a simulation year (Smith et al., 2001), resulting in changes in height, diameter, and biomass growth of each PFT. In ECBilt-CLIO_LPJ-GUESS, LPJ-GUESS is asynchronously coupled to iLOVECLIM and calculates the cover of the 11 PFTs. This asynchronous coupling procedure is identical to the study analyzing climate-vegetation interactions during the LIG (Li et al., 2020). Here, monthly climate inputs from the fully coupled iLOVECLIM model (ECBilt-CLIO-VECODE) are fed to LPJ-GUESS in an initial step, which simulates vegetation distributions. Then, the resulting vegetation distributions are given back to ECBilt as a fixed vegetation component in iLOVECLIM (ECBilt-CLIO_LPJ-GUESS) during the next round of climate simulation. A new climatology from this integration is simulated and subsequently used offline as climate forcing by LPJ-GUESS to produce a new global vegetation distribution that is subsequently used as a boundary condition by iLOVECLIM, and so on. A frequency of vegetation simulation for 10 years is chosen, and four iterations are performed for the asynchronous couplings with this frequency. We treated vegetation results from the last iteration as vegetation responses to climate and analyze vegetation feedbacks to climate. As explained in Li et al. (2019a), the 11 PFTs of LPJ-GUESS were converted to the three vegetation types (trees, grasses, and bare soil) used in iLOVECLIM.

In the configuration with prescribed pre-industrial vegetation (ECBilt-CLIO-LUH2), pre-industrial vegetation from the CMIP LUH2 (2012) dataset is upscaled to the same resolution as ECBilt. The LUH2 (Land-Use Harmonization, Phase 2) project prepared a harmonized set of land use scenarios that smoothly connects the historical reconstructions of land use with future projections in the format required for Earth system models. We applied the land use data from the CMIP LUH2 dataset at 850 CE (Fig. S1 in the Supplement). The land surface albedo is based on this prescribed vegetation and is therefore fixed at a constant value in ECBilt during the simulations.

In the three applied configurations, the land surface albedo and soil hydrology are computed in a simple land surface and

bucket model (LBM) embedded in ECBilt, in which land surface albedo and the maximum water volume of the bucket are a function of vegetation cover. The land surface albedo depends on the percentage of PFT cover in the LBM, which is negatively related to tree cover. In contrast, the maximum water of the bucket is positively related to total vegetation cover. In our experiments, the fractional surface albedo of trees, grassland, and desert are seasonally fixed at a constant value in the LBM, which is at 0.13, 0.20, and 0.33, respectively. The land surface albedo is thus calculated in the LBM based on dynamical vegetation cover in either VECODE or LPJ-GUESS and is then given back to ECBilt in both the fully and asynchronously coupled versions of iLOVECLIM. In ECBilt-CLIO-LUH2, however, the land surface albedo is fixed to prescribed pre-industrial level.

In both ECBilt-CLIO-VECODE and ECBilt-CLIO_LPJ-GUESS, soil hydrology depends on dynamic vegetation modeled in either VECODE or LPJ-GUESS, while soil hydrology depends on prescribed pre-industrial vegetation in ECBilt-CLIO-LUH2. The water exchanges calculated in the full hydrology model in LPJ-GUESS only impact vegetation within LPJ-GUESS, but we keep the soil hydrology calculated in the LBM in the asynchronous coupling of ECBilt-CLIO to LPJ-GUESS. In other words, the hydrology model in LPJ-GUESS affects the vegetation dynamics but the calculation of vegetation feedbacks to soil hydrology is outside LPJ-GUESS, which is in the LBM embedded in ECBilt in the asynchronously coupled configuration (ECBilt-CLIO_LPJ-GUESS).

## 2.2 Experimental design

The LIG transient experiments are from Li et al. (2020). Our Holocene experiments were set up in a similar way as the LIG experiments described in detail by Li et al. (2020). We performed transient simulations (Table 1) with two versions of our model configuration (ECBilt-CLIO-VECODE and ECBilt-CLIO-LUH2) for the Holocene (8.5–1 ka). The Holocene transient experiments start from the equilibrium simulation of 8.5 ka by Li et al. (2019b). The name of each experiment reflects both the period of interest and vegetation component used. For example, HOLVEC and HOLFIX refer to transient experiments during the Holocene with the vegetation component of VECODE and fixed pre-industrial vegetation, respectively. In addition, we performed a series of time-slice simulations (HOL_LPJ) with the asynchronously coupled version (ECBilt-CLIO_LPJ-GUESS) for the Holocene (8–1 ka) at 1 kyr intervals. These simulations were run for 1000 model years per time slice, the last 30 years of which we used for the analysis. The initial climatic forcing for each time-slice simulation is a 100-year monthly mean climatology (surface temperature, precipitation, and cloud cover) from fully coupled transient simulations (HOLVEC), and then the asynchronous coupling starts with vegetation from LPJ-GUESS. Although these
time-slice simulations do not provide a continuous view of
the transitions, the results at 1 kyr intervals cover the full evo-
lution from green to desert in North Africa during both the
LIG and the Holocene.

The transient forcings for experiments during the LIG and the Holocene include variations in orbital parameters and greenhouse gas (GHG) concentrations. The orbital forcings for the Holocene are identical to the PMIP4 protocol (Otto-Bliesner et al., 2017). Changes in insolation received by the Earth are calculated in iLOVECLIM during transient simulations according to Berger (1978). From 8.5 ka to about 1 ka, the NH received more insolation during boreal summer and less insolation during boreal fall and winter compared to pre-industrial values. Both the Holocene and the LIG are characterized by substantial changes in the spatial and temporal distribution of insolation. The magnitude of the climatic precession changes during the LIG is stronger than the Holocene. The mean levels of eccentricity are also higher during the LIG, while the obliquity is higher during the Holocene. These differences result in higher summer insolation at 20° N during the early periods of the LIG than during the Holocene and stronger changes in the annual, seasonal, and latitudinal insolation during the LIG than the Holocene (Fig. 1a, e).

The GHG concentrations during both interglacials were similar to those of the pre-industrial period. The greenhouse gas (GHG) concentrations for the Holocene in our simulations are taken from the PMIP3 protocol (https://pmip3. lsce.ipsl.fr, last access: 3 March 2022) to keep the GHG forcing consistent with the 8.5 ka equilibrium simulation (8.5 ka_VEC) from which the initial conditions were derived. In contrast, the GHG for the LIG is identical to the PMIP4 protocol (Otto-Bliesner et al., 2017). Compared to PMIP3, the GHG concentration has been updated in PMIP4 using new data and a revised chronology that provides a consistent history of the evolution of these gases across the Holocene (Otto-Bliesner et al., 2017). The difference in GHG forcing between PMIP4 and PMIP3 was estimated to be −0.8 W m⁻² by Otto-Bliesner et al. (2017), inducing a mean global cooling by 0.24 ± 0.04°. A linear interpolation was applied to GHG data to get annual values for transient simulations, and 100-year mean values of these interpolated values were used for asynchronously coupled simulations, e.g., mean values from 7.95 to 8.05 ka for the simulation HOL_LPJ in 8 ka.

In our experiments we kept the other boundary conditions
identical to the pre-industrial setup. These include the config-
urations of the ice sheets, the land–sea distribution, and the
solar constant.

3 Results and discussion

3.1 North African climate changes

3.1.1 The LIG

All the LIG experiments simulate relatively large amounts
of precipitation in North Africa during the early LIG
(Fig. 1b), corresponding to the relatively high summer in-
solation at 20° N (Fig. 1a). Background information is pro-
vided in Figs. S2 to S6. During the early LIG (about 127–
125 ka), LIGVEC and LIG_LPJ simulate about twice as
large amounts of precipitation (58 cm yr⁻¹) as LIG_FIX
(28 cm yr⁻¹). The two simulations with dynamical vege-
tation (LIGVEC and LIG_LPJ) thus indicate compara-
ble stronger declines in precipitation than the simulation
with fixed pre-industrial vegetation (LIG_FIX). From 125
to 117 ka, precipitation decreases by 48 cm yr⁻¹ from about
58 cm yr⁻¹ to less than 10 cm yr⁻¹ in North Africa in
LIGVEC/LIG_LPJ, while this decline in LIG_FIX is only
about 15 cm yr⁻¹.

A relatively cooler early LIG (29.0°) in North Africa com-
pared to pre-industrial (30.2°) is simulated in LIG_FIX, and
an upward temperature trend during the LIG is seen in this
simulation (Fig. 1c), which is consistent with multi-model
simulations (Lunt et al., 2013; Bakker et al., 2013, 2014) and
PMIP4 experiments for the early LIG during which North
Africa was cooler by about 1.0° (Williams et al., 2020). Even
without the albedo effect, the precipitation in northern Africa
was still significantly higher in the early part of the inter-
glacial due to the enhanced summer monsoon, forced by el-
evated insolation values. This high precipitation resulted in
relatively humid soils and enhanced evaporation, leading to
evaporative cooling in the first part of the LIG relative to
the end of the LIG_FIX experiment. This created the posi-
tive temperature trend. Different from this positive temper-
ature trend in LIG_FIX, the surface temperature in North
Africa decreases during the LIG in LIGVEC and LIG_LPJ.
A warmer early LIG (about 31.5°) than pre-industrial is seen
in experiments LIGVEC and LIG_LPJ, in particular in the
northwestern part of NA (Li et al., 2020), which confirms
the positive impacts of dynamical vegetation on surface tem-
perature. After 121 ka, the temperature differences between
LIGVEC/LIG_LPJ and LIG_FIX are small, indicating the
end of desertification in North Africa. During the desertifi-
cation phase, between about 123 and 121 ka, the changing
rates of surface temperature reach a maximum in all three
LIG simulations.

3.1.2 The Holocene

Precipitation declines in all simulations (Fig. 1f) during the Holocene, corresponding to the reduction in July insolation at 20° N (Fig. 1e). Similar precipitation trends in North Africa during the Holocene were suggested by previous studies (Renssen et al., 2003; Schurgers et al., 2006). Com-
<table><caption>**Table 1.** Overview of the performed experiments.</caption><thead><tr><th></th><th></th></tr></thead><tbody><tr><td colspan="2">Pre-industrial experiments</td></tr><tr><td>PIVEC</td><td>Pre-industrial run with dynamical vegetation in VECODE (ECBilt-CLIO-VECODE)</td></tr><tr><td>PIFIX</td><td>Pre-industrial run with fixed vegetation from CMIP LUH2 (ECBilt-CLIO-LUH2)</td></tr><tr><td>PLPJ</td><td>Pre-industrial run with asynchronous dynamical vegetation in LPJ-GUESS (ECBilt-CLIO-LPJ-GUESS)</td></tr><tr><td colspan="2">LIG experiments</td></tr><tr><td>127KVEC</td><td>127 ka ctrl run with dynamical vegetation in VECODE (ECBilt-CLIO-VECODE)</td></tr><tr><td>LIGVEC</td><td>The LIG (127–116 ka) run with transient insolation forcing and dynamical vegetation in VECODE (ECBilt-CLIO-VECODE)</td></tr><tr><td>LIGFIX</td><td>The LIG (127–116 ka) run with transient insolation forcing and fixed vegetation from CMIP LUH2 (ECBilt-CLIO-LUH2)</td></tr><tr><td>LIGLPJ</td><td>The LIG (125–117 ka) run with 1 kyr interval insolation forcing and asynchronous dynamical vegetation in LPJ-GUESS (ECBilt-CLIO-LPJ-GUESS)</td></tr><tr><td colspan="2">Holocene experiments</td></tr><tr><td>8.5 kaVEC</td><td>8.5 ka ctrl run with dynamical vegetation in VECODE (ECBilt-CLIO-VECODE)</td></tr><tr><td>HOLVEC</td><td>The Holocene (8.5–1 ka) run with transient insolation forcing and dynamical vegetation in VECODE (ECBilt-CLIO-VECODE)</td></tr><tr><td>HOLFIX</td><td>The Holocene (8.5–1 ka) run with transient insolation forcing and fixed vegetation from CMIP LUH2 (ECBilt-CLIO-LUH2)</td></tr><tr><td>HOLLPJ</td><td>The Holocene (8–1 ka) run with 1 kyr interval insolation forcing and asynchronous dynamical vegetation in LPJ-GUESS (ECBilt-CLIO-LPJ-GUESS)</td></tr></tbody></table>

parable to precipitation features during the LIG, precipitation also saw larger magnitudes of decline in experiments with dynamical vegetation (HOLVEC) than in the experiment with fixed pre-industrial vegetation (HOLFIX). Precipitation gradually decreases from 20 to 16 cm yr⁻¹ during the period (8.5–2 ka) in HOLFIX, corresponding to the gradual reduction in summer insolation at 20° N (Fig. 1e). From 8.5 ka to about 7 ka, the amount of precipitation in North Africa remains at a relatively high level in HOLVEC and HOLLPJ (45 cm yr⁻¹) in comparison to the value (22 cm yr⁻¹) in HOLFIX. Compared to this early period, the amount of precipitation decreases from 40 to 30 cm yr⁻¹ with increased variability during the period between 7 and 5 ka in HOLVEC, and this precipitation trend is more gradual in HOLFIX than in HOLVEC. Subsequently, the annual precipitation decreases gradually to a relatively low level (around 20 cm yr⁻¹) at 2 ka, in particular in experiments HOLVEC and HOLFIX. However, the reduction in precipitation in HOLLPJ is much more modest (from about 40 to 35 cm yr⁻¹) over the 8.5 to 2 ka period, which could be related to the simulated vegetation distribution and spatial divergence in North Africa (Fig. 3).

During the period from 8.5 to 7 ka, a warmer North Africa (31.2°) compared to pre-industrial (30.2°) is suggested by HOLVEC and HOLLPJ (Fig. 1g), but a lower surface temperature is shown in HOLFIX (29.0°). Compared to the temperature evolution in HOLFIX, a subsequent cooling trend from 7 to 2 ka is found in HOLVEC, indicating positive impacts of vegetation feedbacks on temperature consistent with earlier studies (Renssen et al., 2003; Bakker et al., 2014). In contrast, no clear temperature change is simulated

in HOLLPJ, similar to HOLFIX. This difference in surface temperature between HOLVEC and HOLLPJ reflects different magnitudes of vegetation feedbacks under different dynamical vegetation conditions.

## 3.2 North African vegetation evolution and feedbacks

### 3.2.1 The LIG

Vegetated North Africa experienced desertification during the LIG in experiments LIGVEC and LIGLPJ. Compared to LIGFIX, the magnitudes of vegetation feedbacks to climate in both experiments vary with simulated vegetation distributions through vegetation-induced changes in albedo and surface evaporation. During the early LIG (127–125 ka), the Sahara region is covered with vegetation (> 70%) in both LIGVEC and LIGLPJ (Fig. 2a). Vegetation cover mainly declines in the northwestern part of the region from 123 ka onwards (Fig. 3). The decreases in vegetation cover accelerate from 123 ka in these two experiments (Fig. 2a), with a steeper trend in LIGVEC. The sharp decline in vegetation in LIGVEC is simultaneous with a strong reduction in precipitation (Fig. 1b), but also with stronger cooling (Fig. 1c), showing that feedbacks between vegetation and climate are behind the enhanced desertification. The rate of this decline reaches a peak at 122 ka in LIGVEC, while it is later in LIGLPJ (at about 121 ka). This faster vegetation transition in LIGVEC than in LIGLPJ is related to the different complexities of vegetation components, as we discussed in detail in earlier studies (Li et al., 2019a, 2020). We found that the complexity of VECODE and LPJ-GUESS affects vegetation simulations mainly through diversity when
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Total vegetation cover during the LIG (a) and the Holocene (b). The blue lines are simulated results from transient experiments with dynamical vegetation (LIGVEC and HOLVEC). The light green dots with error bars are simulated results from experiments with asynchronously coupled “dynamical” vegetation (LIG LPJ and HOL LPJ). The dark green dots with error bars are vegetation results from LPJ-GUESS forced offline with climatic forcing from LIGVEC and HOLVEC. The blue, red, and light green horizontal lines indicate pre-industrial vegetation cover from PIVEC, PIFIX, and PILPJ, respectively.</figcaption></figure>

the atmospheric CO₂ level is around the pre-industrial level (280 ppmv), while the difference in complexity affects vegetation simulations mainly through ecophysiological processes when the atmospheric CO₂ level is largely different from this pre-industrial value. This finding that higher vegetation diversity leads to a more gradual vegetation change is consistent with previous studies on vegetation transitions during the Holocene (Claussen et al., 2013; Hély et al., 2014; Li et al., 2019a).

Corresponding to the vegetation transitions in the LIG experiments, the impacts of vegetation-climate feedbacks are revealed by comparisons between LIGVEC/LIG LPJ and LIG_FIX. Enhanced vegetation cover (Fig. 2a) accompanied by relatively high levels of precipitation (Fig. 1a, b) is simulated in LIGVEC and LIG LPJ, indicating the effects of vegetation feedbacks on precipitation. These effects are positive with maximum magnitudes during the early LIG due to high vegetation cover (Fig. 2a). This positive feedback of vegetation to precipitation has been associated with changes in albedo (Charney, 1975; Eltahir, 1996; Eltahir and Gong, 1996; Braconnot et al., 1999; Schurgers et al., 2007; Levis et al., 2004; Vamborg et al., 2011). During the LIG, the effect of vegetation feedback on precipitation decreases with the decreasing vegetation cover due to the increasing surface albedo, and this effect on precipitation becomes close to 0 cm yr⁻¹ at 121 kyr after desertification.

### 3.2.2 The Holocene

Vegetation cover (Fig. 2b) decreases gradually from 60% to 5% during the period between 8.5 and 2 ka in HOLVEC, while more fluctuations occur in HOLLPJ. A minor decline in vegetation cover around 6 ka is shown in HOLLPJ, and this decline continues with large variability after 6 ka. Such gradual vegetation declines have also been found in terrestrial records (e.g., Kröpelin et al., 2008; Amaral et al., 2013) and simulations (e.g., Renssen et al., 2006; Schurg-

ers et al., 2006), but they are inconsistent with the abrupt decline suggested by dust records (deMenocal et al., 2000) and some simulation studies (Liu et al., 2007; Notaro et al., 2008). One important reason for this discrepancy is the spatial heterogeneity during the vegetation transition. For example, Renssen et al. (2003) found a gradual vegetation decline in most of the Sahara, yet Liu et al. (2007) suggested abrupt vegetation decrease in the eastern part of the Sahara. The other reason for the gradual vegetation decline and large variability in HOLLPJ could be the asynchronous coupling of LPJ-GUESS to iLOVECLIM. In such coupling, the simulated vegetation is averaged every 10 years and is then passed back to ECBilt as static land surface parameters. Thus, the vegetation transition could be easily smoothed. Likewise, the discrepancy between the vegetation decline simulated by LPJ-GUESS offline (directly forced by climate from HOLVEC) and HOLLPJ also increases from 6 ka, confirming the increasing impacts of the asynchronous coupling on vegetation during the Holocene when the summer insolation and precipitation decline gradually (Fig. 1). The large error bars for the HOLLPJ snapshot experiments for 5 and 4 ka (Fig. 2b) indicate that there was a strong increase in the interannual variability in the simulation of PFT cover. This is likely related to one or more PFTs being very close to their climatic limit such that relatively small variations in precipitation cause large shifts between the area covered by the involved PFT. In this case, bare ground was also involved.

To investigate the spatial heterogeneity of Holocene vegetation transitions in the Sahara, we calculated regional-scale (the western and eastern parts of North Africa) features for climate and vegetation changes in HOLVEC and HOLLPJ. Similar to the evolution of vegetation in the whole Sahara, the vegetation cover declines gradually during the Holocene in the western part of the Sahara (Fig. 4a), but it experiences a slight drop around 6 ka in the eastern part, especially in HOLVEC (Fig. 4b). Consistently, an earlier end of the AHP in the eastern than the western part of the Sa-
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Total vegetation cover anomalies during the LIG (left two columns) and Holocene (right two columns) in the top panel compared to simulated pre-industrial vegetation (PI, bottom panels) from VECODE (left) and LPJ-GUESS (right).</figcaption></figure>

hara is seen in HOLVEC (Fig. 4), indicating considerable variations in the AHP termination in both space and time. This spatial and temporal complexity implies that the termination of the AHP could have been locally abrupt. Regionally different responses to the Holocene insolation decrease are also suggested by Shanahan et al. (2015) and Dallmeyer et

al. (2020, 2021). Moreover, the rate of decline at the end of the AHP indicates complexity. The gradual decline in vegetation cover over western North Africa is consistent with simulations by Renssen et al. (2006) in response to the gradual declines in precipitation and summer insolation at 20° N. In contrast, although the slight drop around 6 ka over eastern
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Evolution of western (a) and eastern (b) North African precipitation, surface temperature, evaporation, soil moisture, and total vegetation cover during the Holocene. The longitudinal boundaries were between 10° W–10° E (western part) and 10–35° E (eastern part), both with latitudinal limits at 15 and 30° N.</figcaption></figure>

North Africa is also seen in simulations by Liu et al. (2007), the magnitude of our vegetation decline is much weaker than in their study that shows a decline of about 30% from 6 to 5 ka. The different magnitudes could be related to the differences in model complexity and also to the selection of a specific area in eastern North Africa. This spatial heterogeneity suggests that the abruptness of vegetation decline is a regional-scale feature (Brovkin and Claussen, 2008).

The magnitudes of vegetation feedback to climate during the Holocene depend on differences in vegetation cover. The difference between HOLVEC and HOLFIX, interpreted here as the impact of dynamical vegetation on precipitation, declines gradually from 27 cm yr⁻¹ at 8.5 ka to 5 cm yr⁻¹ at 2 ka in response to the reduced vegetation cover (Figs. 11, 3).

A similar impact of 25 cm yr⁻¹ on precipitation is shown for HOL_LPJ (Fig. 1f) at 8 ka due to the similar level of vegetation cover as HOLVEC. However, the trends of this vegetation feedback indicate large fluctuations in HOL_LPJ (Figs. 1f, 2b) due to its complex vegetation cover changes after 6 ka. The increased precipitation (roughly a factor of 2) at around 8 ka in HOL_LPJ and HOLVEC compared to the fixed pre-industrial vegetation counterpart (HOLFIX) is consistent with strong vegetation feedbacks suggested by previous simulations (Claussen and Gayler, 1997; Renssen et al., 2003; Rachmayani et al., 2015) and mid-Holocene reconstructions (Bartlein et al., 2011). The strong vegetation feedback to precipitation shows that the pure insolation forcing (as in HOLFIX) is not sufficient to reproduce the reconstructed rainfall, implying that vegetation feedbacks substantially amplify the rainfall response in North Africa in the early Holocene.

## 3.3 Comparison between the LIG and Holocene

The climate and vegetation changes in North Africa during the LIG and the Holocene share many common aspects (Figs. 1 and 2): (1) the amount of precipitation decreases in response to the reduction in summer insolation at 20° N, (2) the land surface in this region experienced a transition from vegetated cover to desert cover, and (3) the positive impacts of vegetation feedback on precipitation reach a maximum during the early period of both interglacials. In addition to these common features, the main differences between the two interglacials are the stronger amplitudes of climate and vegetation changes during the LIG compared to the Holocene.

Based on the results of our experiments, during the early period of both the LIG and Holocene, land surface was dominantly covered by vegetation with only a small proportion of desert in North Africa. This is in agreement with the “green Sahara” suggested by reconstructions (e.g., Prennice and Jolly, 2000) and simulations (e.g., Renssen et al., 2003; Schurgers et al., 2006, 2007). The summer insolation at 20° N around 125 ka is higher by 15 W m⁻² than around 8.5 ka (490 W m⁻²), resulting in a higher amount of precipitation in the LIG over North Africa by about 12 cm yr⁻¹ and slightly higher vegetation cover by 14% around 125 ka than the early Holocene. Compared to the pre-industrial, the higher summer insolation during the early LIG and Holocene induces a stronger summer monsoon, resulting in more moisture transported from the Atlantic to the continent and therefore supporting larger vegetation cover (Fig. 2). In addition to the amount of precipitation induced by the high insolation, vegetation feedbacks also positively affect the amount of precipitation in North Africa. The large vegetation cover amplified regional climate through changes in surface albedo and evaporation, acting as biogeophysical feedbacks to climate. The amounts of precipitation in experiments with dynamical vegetation are about twice the values in experiments
with fixed pre-industrial vegetation during both interglacials when their total vegetation cover is above 60%. This fact indicates that vegetation cover plays a central role in the magnitude of vegetation feedbacks, which is consistent with the highlighted role of vegetation in climate by studies based on fully coupled Earth system models (Braconnot et al., 2012; Shanahan et al., 2015; Dallmeyer et al., 2020).

During the transition period from vegetation to desert, vegetation declined faster during the LIG (about 12% kyr⁻¹) than in the Holocene (about 4% kyr⁻¹), corresponding to a faster reduction in summer insolation at 20° N (compare Fig. 2 with Fig. 1a and e). Indeed, in the LIG between 120 and 118 ka, the summer insolation decreases to a value well below that of the Holocene. This faster transition during the LIG can also be seen in the development of precipitation, evaporation, and soil moisture in all our simulations (Fig. 1). In the Sahara, the amount of precipitation is the most important limiting factor for vegetation growth. In our simulations, the vegetated Sahara is sustained when the amount of precipitation is larger than about 40 cm yr⁻¹ (Fig. 5), which is consistent with earlier studies suggesting a requirement for precipitation of more than 37 cm yr⁻¹ to sustain vegetation in the region between 18–23° N and 11–36° E (Liu et al., 2006; Hopcroft et al., 2017). During the vegetation transition in both interglacials, the reduction in vegetation cover is associated with an increasing surface albedo, leading to decreased amounts of precipitation (Fig. 5), showing the active role of surface albedo in the vegetation feedback to precipitation.

In the transition periods during both interglacials, the amount of precipitation decreases with the reduction in summer insolation at 20° N, and the slope of the relationship between precipitation and summer insolation at 20° N increases when this insolation is greater than about 470 W m⁻² (Fig. 6). The nearly linear relationship between precipitation and summer insolation at 20° N indicates that insolation is the driver of precipitation change. Moreover, the enhanced slopes in the experiments with dynamical vegetation (green and blue circles in Fig. 6) suggest that the positive biogeochemical vegetation feedbacks adjust the amount of precipitation in North Africa in response to the high vegetation cover (Fig. 5a, d). In contrast, several studies found negative vegetation feedbacks to precipitation during the Holocene through the mechanism whereby evaporation from the bare ground appears to be stronger than evapotranspiration from grassland under wet conditions, during which the differences in surface albedo between wet soil and vegetation are small (Liu et al., 2007; Notaro et al., 2008; Wang et al., 2008). However, this negative feedback is not seen in our Holocene experiments, which could be related to the fixed surface albedo of bare soil that is not a function of water content in both our dynamical vegetation experiments, leading to an overestimation of the positive vegetation–albedo feedback. In addition, no seasonal variations in albedo were applied, which likely also affected the strength of this feedback.

In the LIG the simulated precipitation was lower than in the Holocene, even with similar insolation values at 20° N (Fig. 6). This is likely related to differences in the thermal contrast between the Atlantic Ocean and the continent, as well as to differences in the strength of vegetation feedbacks. The thermal contrast determines the strength of the monsoonal precipitation, and since the Atlantic Ocean was slightly warmer in the LIG than in the Holocene, this could partly explain stronger precipitation in the Holocene with the same summer insolation. In addition, comparing 122 ka with 8 ka (Fig. 3), the vegetation cover is more extensive in the early Holocene, which enforces the precipitation through the vegetation–climate feedbacks.

Although the positive vegetation feedbacks enhance the amount of precipitation in North Africa during both interglacials, the mechanism of vegetation transitions during the two interglacials could be different, in particular during the transition periods. During the transition in the LIG (Figs. 1b, 2a), the declines in vegetation and precipitation are characterized by large variability, but these reductions during the Holocene are relatively gradual (Figs. 1f, 2b). Based on a series of sensitivity experiments by Liu et al. (2006), vegetation collapse in North Africa could occur either when the vegetation feedback to precipitation is relatively strong with rapid declines in precipitation or when the feedback is weak with gradually decreasing precipitation. In their study, the vegetation feedback to precipitation is strong when the amount of vegetation-induced precipitation is greater than the lowest requirement of precipitation for a vegetated Sahara. In our case, the vegetation feedbacks to precipitation are greater than the lowest requirement of precipitation for a vegetated Sahara during the early periods of both interglacials. Therefore, the abrupt vegetation transition during the LIG is likely to be due to the combination of the strong vegetation feedback and rapidly decreasing precipitation, while the gradual vegetation transition during the Holocene is a result of the combination of strong vegetation feedback with gradually decreasing precipitation. Moreover, during both interglacials, larger variability in precipitation and vegetation is seen in experiments with dynamic vegetation in VECODE than in LPJ-GUESS. One reason for this difference is the higher complexity of LPJ-GUESS, which simulates more gradual vegetation responses when more plant types interact with climate, consistent with other studies (Claussen et al., 2013; Hély et al., 2014; Li et al., 2019a).

## 3.4 Uncertainties

One source of uncertainty is related to the spatial vegetation variability in North Africa, in particular during the Holocene when the reduction in summer insolation at 20° N is rather gradual compared to the LIG. In order to find the signal of abrupt vegetation decline in Sahara during the Holocene, we investigated the evolutions of vegetation and climate over the western and eastern part of North Africa. We find vegetation
<figure><img src="image_7.png" /><figcaption>**Figure 5.** The relationship between (a) precipitation and vegetation cover; (b) vegetation cover and surface albedo; (c) surface albedo and precipitation with dynamical vegetation from VECODE during two interglacials; (d) precipitation and vegetation cover; (e) vegetation cover and surface albedo; (f) surface albedo and precipitation with dynamical vegetation from LPJ-GUESS during two interglacials.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 6.** The relationship between insolation at 20° N and precipitation with dynamic vegetation from LPJ-GUESS (a) and VECODE (b) during two interglacials, with blue dots representing the LIG and green dots representing the Holocene. For reference, the results of HOL_FIX, LIG_FIX, and PI are also included.</figcaption></figure>

declines consistent with previous studies and regional-scale features of vegetation changes. However, the coarse resolution of our climate model limits our scope of spatial variations in vegetation distribution, implying, for instance, that regional abrupt vegetation changes are not well represented. This tends to smooth the simulated results and affects the abruptness of transitions, especially when the external forcing (insolation) is around the threshold of vegetation transition. In addition, Bathiany et al. (2012) pointed out the importance of spatial interactions; i.e., abrupt changes in one region can be induced by critical transitions in the neighboring region. These abrupt changes can occur in one region triggered by strong vegetation–atmosphere interaction in a neighboring region, which is called “inducing tipping of

ecosystems” by Bathiany et al. (2013a, b). However, this inducing tipping of ecosystems is not seen in our simulations. We therefore recommend more studies on regional-scale vegetation–climate interactions by fully coupled AOVGCMs (atmosphere–ocean–vegetation general circulation models), as well as comparisons with results from EMICs (Earth system models of intermediate complexity) to explore the regional-scale features and to improve our understanding of climate and vegetation sensitivity to these small-scale features in North Africa. In addition, these model results should be thoroughly evaluated using proxy-based reconstructions of the green–desert transition during the Holocene and LIG.

A second major source of uncertainty is related to the land surface parameters (surface albedo and vegetation cover)
provided to ECBilt, resulting in vegetation feedbacks to climate. In our study, the value of albedo for each PFT is prescribed, implying that this value in each grid cell is decided by fractions of PFTs and is passed back to ECBilt through full or asynchronous coupling. The prescribed albedo value for bare soil (0.33) could exaggerate the simulated positive vegetation–albedo feedback during the interglacials, especially during the transition periods when the soil moisture could be higher than typically found in deserts. This overestimation is obvious, in particular during the Holocene when the summer insolation at 20°N declines rather gradually. Therefore, more studies on albedo values of bare soil with different water retentions and the parameterizations of these mechanisms in climate models could improve the understanding of vegetation–climate interactions.

In addition, there are two vegetation models in this study. VECODE, fully coupled to ECBilt, suggests more abrupt vegetation and climate changes than LPJ-GUESS, which is asynchronously coupled to ECBilt during both interglacials. The main reason for this difference in our results could be related to their different complexity. The faster vegetation transition in VECODE confirms the hypothesis that the low diversity implies a high likelihood of abrupt transitions (Scherer-Lorenzen, 2005; Claussen et al., 2013; Hély et al., 2014; Li et al., 2019a). Moreover, the simulated vegetation in the asynchronously coupled LPJ-GUESS is averaged every 10 years and is then passed back to ECBilt as static land surface parameters. Hence, the transition from green Sahara to desert is easily smoothed, and the feedback is thereby underestimated. In order to narrow down the uncertainty in terms of the asynchronous coupling during simulations, we prefer the full dynamical coupling of LPJ-GUESS to ECBilt in the future, which could enhance the robustness in climate–vegetation interaction simulations.

## 4 Conclusions

In this study, we compare the climate and vegetation changes over North Africa during the Holocene and the LIG using the climate model iLOVECLIM in which the vegetation component is either the fully coupled VECODE model or the asynchronously coupled LPJ-GUESS model. We address three scientific questions. (1) Are the evolutions of NA vegetation during the two interglacials simulated by iLOVECLIM with different vegetation components consistent? (2) Are the evolutions of NA vegetation different during the Holocene and LIG, and what are the dominating factors? (3) How do the NA vegetation geophysical feedbacks affect the climate during the Holocene and LIG?

During both interglacials, we find overall agreement with a transition from a vegetated Sahara (with more than about 60% vegetation cover) to desert and positive vegetation feedback to climate. The simulations with a relatively complex vegetation component (LPJ-GUESS) suggest a gradual veg-

estation transition compared to simulations with a simple vegetation component (VECODE). This more gradual vegetation transition in simulations with LPJ-GUESS is a result of the smoothness of the desertification trend due to the increasing vegetation diversity and features of different PFTs in the more complex vegetation component. Except for the impacts of the applied vegetation model, the simulated Sahara vegetation transitions are nearly linearly related to the decreasing summer insolation at 20°N. Due to the faster decrease in 20°N summer insolation during the LIG than in the Holocene, the vegetation cover experienced a stronger reduction in the LIG. The decline of vegetation cover peaks at 25% kyr⁻¹ at around 122 ka, while the steepest vegetation cover change is 10% kyr⁻¹ at around 6 ka.

In both interglacials, the abruptness of vegetation changes is nearly linearly related to the changes in summer insolation at 20°N, and the effects of vegetation feedback to climate positively depend on vegetation cover. Vegetation feeds back to precipitation by a factor of 2 to 3 when the vegetation cover is greater than 60% during the early LIG and Holocene through vegetation–albedo feedbacks. Due to the different magnitudes of declines in summer insolation at 20°N, the abruptness of vegetation transitions during the LIG and the Holocene is different. During the LIG, the abrupt transition is a result of strong vegetation feedback combined with rapidly decreasing precipitation, while the gradual transition during the Holocene is related to the strong vegetation feedback and gradually decreasing precipitation.

**Code and data availability.** The data are available from the authors upon request.

**Supplement.** The supplement related to this article is available online at: https://doi.org/10.5194/cp-18-2303-2022-supplement.

**Author contributions.** The study was designed by all authors. HL performed the simulations and wrote the paper with contributions of DMR and HR. All authors contributed to the analysis and interpretation of the model results.

**Competing interests.** The contact author has declared that none of the authors has any competing interests.

**Disclaimer.** Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Acknowledgements.** This work was supported by the China Scholarship Council, the National Natural Science Foundation of China (4210010473), the Natural Science Foundation of Jiangsu
Province, China (BK20210848), the Natural Science Foundation of Nantong (JC2021161), and the Jiangsu Provincial Double-Innovation Doctor Program (JSSCBS20211094). DMR is supported by CNRS/INSU and by the Vrije Universiteit Amsterdam. HR and DMR are supported through the European Union's Horizon 2020 research and innovation program within the TERRANOVA project (grant no. 813904).

**Financial support.** This research has been supported by the China Scholarship Council and the National Natural Science Foundation of China (grant no. 4210010473), the Natural Science Foundation of Jiangsu Province (grant no. BK20210848), the Natural Science Foundation of Nantong (grant no. JC2021161), the Jiangsu Provincial Double-Innovation Doctor Program (grant no. JSSCBS20211094), and the European Union's Horizon 2020 research and innovation program within the TERRANOVA project (grant no. 813904).

**Review statement.** This paper was edited by Martin Claussen and reviewed by Qiong Zhang and one anonymous referee.

## References

Amaral, P. G. C., Vincens, A., Guiot, J., Buchet, G., Deschamps, P., Doumnang, J.-C., and Sylvestre, F.: Palynological evidence for gradual vegetation and climate changes during the African Humid Period termination at 13°N from a Mega-Lake Chad sedimentary sequence, Clim. Past, 9, 223–241, https://doi.org/10.5194/cp-9-223-2013, 2013.

Bakker, P., Stone, E. J., Charbit, S., Gröger, M., Krebs-Kanzow, U., Ritz, S. P., Varma, V., Khon, V., Lunt, D. J., Mikolajewicz, U., Prange, M., Renssen, H., Schneider, B., and Schulz, M.: Last interglacial temperature evolution – a model inter-comparison, Clim. Past, 9, 605–619, https://doi.org/10.5194/cp-9-605-2013, 2013.

Bakker, P., Masson-Delmotte, V., Martrat, B., Charbit, S., Renssen, H., Gröger, M., Krebs-Kanzow, U., Lohmann, G., Lunt, D. J., Pfeiffer, M., Phipps, S. J., Prange, M., Ritz, S. P., Schulz, M., Stenni, B., Stone, E. J., and Varma, V.: Temperature trends during the Present and Last Interglacial periods – a multi-model-data comparison, Quaternary Sci. Rev., 99, 224–243, 2014.

Bartlein, P. J., Harrison, S. P., Brewer, S., Connor, S., Davis, B. A. S., Gajewski, K., Guiot, J., Harrison-Prentice, T. I., Henderson, A., Peyron, O., Prentice, I. C., Scholze, M., Seppä, H., Shuman, B., Sugita, S., Thompson, R. S., Viau, A. E., Williams, J., and Wu, H.: Pollen-based continental climate reconstructions at 6 and 21 ka: a global synthesis, Clim. Dynam., 37, 775–802, 2011.

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, 2012.

Bathiany, S., Claussen, M., and Fraedrich, K.: Detecting hotspots of atmosphere-vegetation interaction via slowing down – Part 1: A stochastic approach, Earth Syst. Dynam., 4, 63–78, 2013a.

Bathiany, S., Claussen, M., and Fraedrich, K.: Detecting hotspots of atmosphere-vegetation interaction via slowing down – Part 2:

Application to a global climate model, Earth Syst. Dynam., 4, 79–93, https://doi.org/10.5194/esd-4-79-2013, 2013b.

Berger, A.: Long-term variations of caloric insolation resulting from the earth's orbital elements, Quaternary Res., 9, 139–167, 1978.

Braconnot, P., Joussaume, S., Marti, O., and de Noblet, N.: Syner-gistic feedbacks from ocean and vegetation on the African Mon-soon response to Mid-Holocene insolation, Geophys. Res. Lett., 26, 2481–2484, 1999.

Braconnot, P., Harrison, S. P., Kageyama, M., bartlein, P. J., Masson-Delmotte, V., Abe-Ouchi, A., Otto-Bliesner, B., and Zhao, Y.: Evaluation of climate models using palaeoclimatic data, Nat. Clim. Change 2, 417–424, 2012.

Brovkin, V. and Claussen, M.: Comment on “Climate-Driven Ecosystem Succession in the Sahara: The Past 6000 Years”, Sci-ence, 322, 1326–1326, 2008.

Brovkin, V., Ganopolski, A., and Svirezhev, Y.: A continuous climate-vegetation classification for use in climate-biosphere studies, Ecol. Model., 101, 251–261, 1997.

Charney, J. G.: Dynamics of deserts and drought in Sahel, Q. J. Roy. Meteor. Soc., 101, 193–202, 1975.

Claussen, M. and Gayler, V.: The greening of the Sahara during the mid-Holocene: results of an interactive atmosphere-biome model, Global Ecol. Biogeogr., 6, 369–377, 1997.

Claussen, M., Kubatzki, K., Brovkin, V., Ganopolski, A., Hoelz-mann, P., and Pachur, H.-J.: Simulation of an abrupt change in Saharan vegetation in the Mid-Holocene, Geophys. Res. Lett., 26, 2037–2040, 1999.

Claussen, M., Bathiany, S., Brovkin, V., and Kleinen, T.: Simulated climate-vegetation interaction in semi-arid regions affected by plant diversity, Nat. Geosci., 6, 954–958, 2013.

Claussen, M., Dallmeyer, A., and Bader, J.: Theory and mod-eling of the African humid period and the green Sa-hara, Oxford Research Encyclopedia of Climate Science, https://doi.org/10.1093/acrefore/9780190228620.013.532, 2017.

Dallmeyer, A., Claussen, M., Lorenz, S. J., and Shanahan, T.: The end of the African humid period as seen by a transient comprehensive Earth system model simulation of the last 8000 years, Clim. Past, 16, 117–140, https://doi.org/10.5194/cp-16-117-2020, 2020.

Dallmeyer, A., Claussen, M., Lorenz, S. J., Sigl, M., Toohey, M., and Herzschuh, U.: Holocene vegetation transitions and their climatic drivers in MPI-ESM1.2, Clim. Past, 17, 2481–2513, https://doi.org/10.5194/cp-17-2481-2021, 2021.

deMenocal, P., Ortiz, J., Guilderson, T., Adkins, J., Sarnthein, M., Baker, L., and Yarusinsky, M.: Abrupt onset and termination of the African Humid Period: Rapid climate responses to gradual insolation forcing, Quaternary Sci. Rev., 19, 347–361, 2000.

Ehrmann, W., Schmiedl, G., Beuscher, S., and Krüger, S.: Intensity of African Humid Periods Estimated from Saharan Dust Fluxes, PLoS ONE, 12, e0170989, https://doi.org/10.1371/journal.pone.0170989, 2017.

Eltahir, E. A. B.: Role of vegetation in sustaining large-scale at-mospheric circulations in the tropics., J. Geophys. Res.-Atmos., 101, 4255–4268, 1996.

Eltahir, E. A. B. and Gong, C.: Dynamics of wet and dry years in West Africa, J. Climate, 9, 1030–1042, 1996.

Fischer, N. and Jungclaus, J. H.: Effects of orbital forcing on at-mosphere and ocean heat transports in Holocene and Eemian climate simulations with a comprehensive Earth system model,
Clim. Past, 6, 155–168, https://doi.org/10.5194/cp-6-155-2010, 2010.

Fleitmann, D., Burns, S. J., Pekala, M., Mangini, A., Al-Subbary, A., Al-Aowah, M., Kramers, J., and Matter, A.: Holocene and Pleistocene pluvial periods in Yemen, southern Arabia, Quaternary Sci. Rev., 30, 783–787, 2011.

Francus, P., von Suchodoletz, H., Dietze, M., Donner, R. V., Bouchard, F., Roy, A. -J., Fagot, M., Verschuren, D., Kröpelin, S., and Ariztegui, D.: Varved sediments of Lake Yoa (Ounianga Kebir, Chad) reveal progressive drying of the Sahara during the last 6100 years, Sedimentology, 60, 911–934, 2013.

Goosse, H., Brovkin, V., Fichefet, T., Haarsma, R., Huybrechts, P., Jongma, J., Mouchet, A., Selten, F., Barriat, P.-Y., Campin, J.-M., Deleersnijder, E., Driesschaert, E., Goelzer, H., Janssens, I., Loutre, M.-F., Morales Maqueda, M. A., Opsteegh, T., Mathieu, P.-P., Munhoven, G., Pettersson, E. J., Renssen, H., Roche, D. M., Schaeffer, M., Tartinville, B., Timmermann, A., and Weber, S. L.: Description of the Earth system model of intermediate complexity LOVECLIM version 1.2, Geosci. Model Dev., 3, 603–633, https://doi.org/10.5194/gmd-3-603-2010, 2010.

Gröger, M., Maier-Reimer, E., Mikolajewicz, U., Schurgers, G., Vizcaíno, M., and Winguth, A.: Vegetation-climate feedbacks in transient simulations over the last interglacial (128000–113000 yr BP), in: The climate of past interglacials, edited by: Sirocko, F., Claussen, M., Sanchez-Goni, M. F., and Litt, T., Develop. Quaternary Sci., 7, 563–572, 2007.

Haxeltine, A. and Prentice, I.C.: BIOME3: An equilibrium terrestrial biosphere model based on ecophysiological constraints, resource availability, and competition among plant functional types, Glob. Biogeochem. Cy., 10, 693–709, https://doi.org/10.1029/96GB02344, 1996.

Hély, C., Lézine, A.-M., and contributors, A.: Holocene changes in African vegetation: tradeoff between climate and water availability, Clim. Past, 10, 681–686, https://doi.org/10.5194/cp-10-681-2014, 2014.

Hoelzmann, P., Gasse, F., Dupont, L. M., Salzmann, U., Staubwasser, M., Leuschner, D. C., and Sirocko, F.: Paleoenvironmental changes in the arid and subarid belt (Sahara–Sahel–Arabian Peninsula) from 150 kyr to present, in: Past Climate Variability through Europe and Africa, edited by: Battarbee, R. W., Gasse, F., and Stickley, C. E., Springer, Dordrecht, 219–256, 2004.

Hopcroft, P. O. and Valdes, P. J.: Paleoclimate-conditioning reveals a North Africa land–atmosphere tipping point, P. Natl. Acad. Sci. USA, 118, e2108783118, https://doi.org/10.1073/pnas.2108783118, 2021.

Hopcroft, P. O., Valdes, P. J., Harper, A. B., and Beerling, D.J.: Multi vegetation model evaluation of the Green Sahara climate regime, Geophys. Res. Lett. 44, 6804–6813, 2017.

Jolly, D., Prentice, I. C., Bonnefille, R., Ballouche, A., Bengo, M., Brenac, P., Buchet, G., Burney, D., Cazet, J. P., Cheddadi, R., Edorh, T., Elenga, H., Elmoutaki, S., Guiot, J., Laarif, F., Lamb, H., Lezine, A. M., Maley, J., Mbenza, M., Peyron, O., Reille, M., Reynaud-Farrera, I., Riollet, G., Ritchie, J. C., Roche, E., Scott, L., Semmanda, I., Straka, H., Umer, M., Van Campo, E., Vilmumbalo, S., Vincens, A., and Waller, M.: Biome reconstruction from pollen and plant macrofossil data for Africa and the Arabian peninsula at 0 and 6000 years, J. Biogeogr., 25, 1007–1027, 1998.

Knorr, W. and Schnitzler, K.-G.: Enhanced albedo feedback in North Africa from possible combined vegetation and soil-formation processes, Clim. Dynam., 26, 55–63, 2005

Kröpelin, S., Verschuren, D., Lézine, A. M., Eggermont, H., Cocquyt, C., Francus, P., Cazet, J.-P., Fagot, M., Rumes, B., Russell, J. M., Darius, F., Conley, D. J., Schuster, M., von Suchodoletz, H., and Engstrom, D. R.: Climate-driven ecosystem succession in the Sahara: The past 6000 years, Science, 320, 765–768, 2008.

Kutzbach, J. E. and Street-Perrott, F. A.: Milankovitch forcing of fluctuations in the level of tropical lakes from 18 to 0 kyr BP, Nature, 317, 130–134, 1985.

Levis, S., Bonan, G. B., and Bonfils, C.: Soil feedback drives the mid-Holocene North African monsoon northward in fully coupled CCSM2 simulations with a dynamic vegetation model, Clim. Dynam., 23, 791–802, 2004

Lézine, A. M.: Vegetation at the Time of the African Humid Period. Regional and Local Climates, Oxford University Press, Oxford Research Encyclopedia of Climate Science, https://doi.org/10.1093/acrefore/9780190228620.013.530, 2017.

Lézine, A. M., Hély, C., Grenier, C., Braconnot, P., and Krinner, G.: Sahara and Sahel vulnerability to climate changes, lessons from Holocene hydrological data, Quaternary Sci. Rev., 30, 3001–3012, 2011.

Li, H., Renssen, H., and Roche, D. M.: Global vegetation distribution driving factors in two Dynamic Global Vegetation Models of contrasting complexities, Global Planet. Change, 180, 51–65, 2019a.

Li, H., Renssen, H., Roche, D. M., and Miller, P. A.: Modelling the vegetation response to the 8.2 ka BP cooling event in Europe and Northern Africa, J. Quaternary Sci., 34, 650–661, https://doi.org/10.1002/jqs.3157, 2019b.

Li, H., Renssen, H., and Roche, D. M.: Modelling climate-vegetation interactions during the Last Interglacial: the impact of biogeophysical feedbacks in North Africa, Quaternary Sci. Rev., 249, 106609, https://doi.org/10.1016/j.quascirev.2020.106609, 2020.

Liu, Z., Wang, Y., Gallimore, R., Notaro, M. and Prentice, I. C.: On the cause of abrupt vegetation collapse in North Africa during the Holocene: Climate variability vs. vegetation feedback, Geophys. Res. Lett., 33, L22709, https://doi.org/10.1029/2006GL028062, 2006.

Liu, Z., Wang, Y., Gallimore, R., Gasse, F., Johnson, T., deMenocal, P., Adkins, J., Notaro, M., Prentice, I. C., Kutzbach, J., Jacob, R., Behling, P., Wang, L., and Ong, E.: Simulating the transient evolution and abrupt change of Northern Africa atmosphere-ocean-terrestrial ecosystem in the Holocene, Quaternary Sci. Rev., 26, 1818–1837, 2007.

LUH2: Land-Use Harmonisation 2., World Climate research Program Coupled Intercomparison Project (CMIP6), https://luh.umd.edu (last access: 3 March 2022), 2012.

Lunt, D. J., Abe-Ouchi, A., Bakker, P., Berger, A., Braconnot, P., Charbit, S., Fischer, N., Herold, N., Jungclaus, J. H., Khon, V. C., Krebs-Kanzow, U., Langebroek, P. M., Lohmann, G., Nisan-cioglu, K. H., Otto-Bliesner, B. L., Park, W., Pfeiffer, M., Phipps, S. J., Prange, M., Rachmayani, R., Renssen, H., Rosenbloom, N., Schneider, B., Stone, E. J., Takahashi, K., Wei, W., Yin, Q., and Zhang, Z. S.: A multi-model assessment of last interglacial temperatures, Clim. Past, 9, 699–717, https://doi.org/10.5194/cp-9-699-2013, 2013.
McKay, N. P., Overpeck, J. T., and Otto-Bliesner, B. L.: The role of ocean thermal expansion in Last Interglacial sea level rise, Geophys. Res. Lett., 38, L14605, https://doi.org/10.1029/2011GL048280, 2011.

Notaro, M., Wang, Y., Liu, Z., Gallimore, R., and Levis, S.: Combined statistical and dynamical assessment of simulated vegetation-rainfall interactions in North Africa during the mid-Holocene1, Glob. Change Biol., 14, 347–368, 2008.

Opsteegh, J. D., Haarsma, R. J., Selten, F. M., and Kattenberg, A.: ECBILT: a dynamic alternative to mixed boundary conditions in ocean models, Tellus A, 50, 348–367, 1998.

Osborne, A. H., Vance, D., Rohling, E. J., Barton, N., Rogerson, M., and Fello, N.: A humid corridor across the Sahara for the migration of early modern humans out of Africa 120,000 years ago, P. Natl. Acad. Sci. USA, 105, 16444–16447, 2008.

Otto-Bliesner, B. L., Braconnot, P., Harrison, S. P., Lunt, D. J., Abe-Ouchi, A., Albani, S., Bartlein, P. J., Capron, E., Carlson, A. E., Dutton, A., Fischer, H., Goelzer, H., Govin, A., Haywood, A., Joos, F., LeGrande, A. N., Lipscomb, W. H., Lohmann, G., Mahowald, N., Nehrbass-Ahles, C., Pausata, F. S. R., Peterschmitt, J.-Y., Phipps, S. J., Renssen, H., and Zhang, Q.: The PMIP4 contribution to CMIP6 – Part 2: Two interglacials, scientific objective and experimental design for Holocene and Last Interglacial simulations, Geosci. Model Dev., 10, 3979–4003, https://doi.org/10.5194/gmd-10-3979-2017, 2017.

Pausata, F. S. R., Messori, G., and Zhang, Q.: Impacts of dust reduction on the northward expansion of the African monsoon during the Green Sahara period, Earth Planet. Sci. Lett., 434, 298–307, 2016.

Prentice, I. C. and Jolly, D.: Biome 6000 participants: Mid-Holocene and glacial-maximum vegetation geography of the northern continents and Africa, J. Biogeogr., 27, 507–519, 2000.

Rachmayani, R., Prange, M., and Schulz, M.: North African vegetation-precipitation feedback in early and mid-Holocene climate simulations with CCSM3-DGVM, Clim. Past, 11, 175–185, https://doi.org/10.5194/cp-11-175-2015, 2015.

Rachmayani, R., Prange, M., and Schulz, M.: Intra-interglacial climate variability: model simulations of Marine Isotope Stages 1, 5, 11, 13, and 15, Clim. Past, 12, 677–695, https://doi.org/10.5194/cp-12-677-2016, 2016.

Renssen, H., Brovkin, V., Fichefet, T., and Goosse, H.: Holocene climate instability during the termination of the African Humid Period, Geophys. Res. Lett., 30, 1184, https://doi.org/10.1029/2002GL016636, 2003.

Renssen, H., Brovkin, V., Fichefet, T., and Goosse, H.: Simulation of the Holocene climate evolution in Northern Africa: The termination of the African Humid Period, Quaternary Int., 150, 95–102, 2006.

Ritchie, J. C., Eyles, C. H., and Haynes, C. V.: Sediment and pollen evidence for an early to mid-Holocene humid period in the eastern Sahara, Nature, 314, 352–355, 1985.

Roche, D. M., Dumas, C., Bügelmayer, M., Charbit, S., and Ritz, C.: Adding a dynamical cryosphere to iLOVECLIM (version 1.0): coupling with the GRISLI ice-sheet model, Geosci. Model Dev., 7, 1377–1394, https://doi.org/10.5194/gmd-7-1377-2014, 2014.

Rohling, E. J., Cane, R., Cooke, S., Sprovieri, M., Bouloubassi, I., Emeis, K. C., Schiebel, R., Kroon, D., Jorissen, F. J., Lorre, A., and Kemp, A. E. S.: African monsoon variability during the pre-

vious interglacial maximum, Earth Planet. Sci. Lett., 202, 61–75, 2002.

Rossow, W. B., Walker, A. W., Beuschel, D. E., and Roiter, M. D.: International satellite cloud climatology project (ISCCP) documentation of new cloud datasets, WMO/TD-No 737, World Meteorological Organisation, https://isccp.giss.nasa.gov/docs/documents.html (last access: 3 March 2022), 76 pp., 1996.

Scherer-Lorenzen, M.: Biodiversity and ecosystem functioning: basic principles, in: Biodiversity: Structure and Function, edited by: Barthlott, W., Linsenmair, K. E., and Porembski, S., Encyclopedia of Life Support Systems (EOLSS), http://www.eolss.net/ (last access: 3 March 2022), 2005.

Schurgers, G., Mikolajewicz, U., Gröger, M., Maier-Reimer, E., Vizcaíno, M., and Winguth, A.: Dynamics of the terrestrial biosphere, climate and atmospheric CO2 concentration during interglacials: a comparison between Eemian and Holocene, Clim. Past, 2, 205–220, https://doi.org/10.5194/cp-2-205-2006, 2006.

Schurgers, G., Mikolajewicz, U., Gröger, M., Maier-Reimer, E., Vizcaíno, M., and Winguth, A.: The effect of land surface changes on Eemian climate. Clim. Dynam., 29, 357–373, 2007.

Scussolini, P., Bakker, P., Guo, C., Stepanek, C., Zhang, Q., Braconnot, P., Cao, J., Guarino, M.-V., Coumou, D., Prange, M., Ward, P. J., Renssen, H., Kageyama, M., Otto-Bliesner, B., and Aerts, J. C. J. H.: Agreement between reconstructed and modeled boREAL precipitation of the Last Interglacial, Sci. Adv., 5, eaax7047, https://doi.org/10.1126/sciadv.aax7047, 2019.

Shanahan, T. M., McKay, N. P., Hughen, K. A., Overpeck, J. T., Otto-Bliesner, B., Heil, C. W., King, J., Scholz, C. A., and Peck, J.: The time-transgressive termination of the African Humid Period, Nat. Geosci., 8, 140–144, https://doi.org/10.1038/NGEO2329, 2015.

Smith, B., Prentice, I. C., and Sykes, M. T.: Representation of vegetation dynamics in the modelling of terrestrial ecosystems: comparing two contrasting approaches within European climate space, Global Ecol. Biogeogr., 10, 621–637, 2001.

Smith, B., Wårlind, D., Arneth, A., Hickler, T., Leadley, P., Siltberg, J., and Zaehle, S.: Implications of incorporating N cycling and N limitations on primary production in an individual-based dynamic vegetation model, Biogeosciences, 11, 2027–2054, https://doi.org/10.5194/bg-11-2027-2014, 2014.

Texier, D., de Noblet, N., Harrison, S. P., Haxeltine, A., Jolly, D., Joussaume, S., Laarif, F., Prentice, I. C., and Tarasov, P.: Quantifying the role of biosphere-atmosphere feedbacks in climate change: coupled model simulations for 6000 years BP and comparison with palaeodata for northern Eurasia and northern Africa, Clim. Dynam., 13, 865–882, 1997.

Tierney, J., Pausata, F. S. R., and DeMenocal, P. B.: Rainfall regimes of the Green Sahara, Science Advances, 3, e1601503, https://doi.org/10.1126/sciadv.1601503, 2017.

Tuenter, E., Weber, S. L., Hilgen, F. J., and Lourens, L. J.: Simulating sub-Milankovitch climate variations associated with vegetation dynamics, Clim. Past, 3, 169–180, https://doi.org/10.5194/cp-3-169-2007, 2007.

Turney, C. S. M. and Jones, R. T.: Does the Agulhas Current amplify global temperatures during super-interglacials? J. Quaternary Sci., 25, 839–843, 2010.

Vamborg, F. S. E., Brovkin, V., and Claussen, M.: The effect of a dynamic background albedo scheme on Sahel/Sahara pre-
cipitation during the mid-Holocene, Clim. Past, 7, 117–131,
https://doi.org/10.5194/cp-7-117-2011, 2011.

Wang, Y., Mysak, L. A., Wang, Z., and Brovkin, V.: The green-
ing of the McGill Paleoclimate Model. Part II: Simulation of
Holocene millennial-scale natural climate changes, Clim. Dy-
nam., 24, 481–496, 2005.

Wang, Y., Notaro, M., Liu, Z., Gallimore, R., Levis, S., and Kutzbach, J. E.: Detecting vegetation-precipitation feedbacks in mid-Holocene North Africa from two climate models, Clim. Past, 4, 59–67, https://doi.org/10.5194/cp-4-59-2008, 2008.

Williams, C. J. R., Guarino, M.-V., Capron, E., Malmierca-Vallet, I., Singarayer, J. S., Sime, L. C., Lunt, D. J., and Valdes, P. J.: CMIP6/PMIP4 simulations of the mid-Holocene and Last Interglacial using HadGEM3: comparison to the pre-industrial era, previous model versions and proxy data, Clim. Past, 16, 1429–1450, https://doi.org/10.5194/cp-16-1429-2020, 2020.