Journal topic
Clim. Past, 15, 449–462, 2019
https://doi.org/10.5194/cp-15-449-2019
Clim. Past, 15, 449–462, 2019
https://doi.org/10.5194/cp-15-449-2019

Research article 19 Mar 2019

Research article | 19 Mar 2019

# The response of tropical precipitation to Earth's precession: the role of energy fluxes and vertical stability

The response of tropical precipitation to Earth's precession: the role of energy fluxes and vertical stability
Chetankumar Jalihal1,2, Joyce Helena Catharina Bosmans3,a,b, Jayaraman Srinivasan2, and Arindam Chakraborty1,2 Chetankumar Jalihal et al.
• 1Centre for Atmospheric and Oceanic Sciences, Indian Institute of Science, Bangalore, India
• 2Divecha Centre for Climate Change, Indian Institute of Science, Bangalore, India
• 3Department of Environmental Science, Radboud University, Nijmegen, the Netherlands
• aformerly at: Faculty of Geosciences, Utrecht University, Utrecht, the Netherlands
• bformerly at: Royal Netherlands Meteorological Institute (KNMI), De Bilt, the Netherlands

Correspondence: Chetankumar Jalihal (jalihal@iisc.ac.in)

Abstract

The changes in Earth's precession have an impact on the tropical precipitation. This has been attributed to the changes in seasonal solar radiation at the top of the atmosphere. The primary mechanism that has been proposed is the change in thermal gradient between the two hemispheres. This may be adequate to understand the zonal mean changes, but cannot explain the variations between land and oceans. We have used a simple model of the intertropical convergence zone (ITCZ) to unravel how precipitation changes with precession. Our model attributes the changes in precipitation to the changes in energy fluxes and vertical stability. We include the horizontal advection terms in this model, which were neglected in the earlier studies. The final response of the land and oceans is a result of complex feedbacks triggered by the initial changes in the insolation. We find that the changes in precipitation over the land are mainly driven by changes in insolation, but over the oceans, precipitation changes on account of changes in surface fluxes and vertical stability. Hence insolation can be a trigger for changes in precipitation on orbital timescales, but surface energy and vertical stability play an important role too. The African monsoon intensifies during a precession minimum (higher summer insolation). This intensification is mainly due to the changes in vertical stability. The precipitation over the Bay of Bengal decreases for minimum precession. This is on account of a remote response to the enhanced convective heating to the west of the Bay of Bengal. This weakens the surface winds and thus leads to a decrease in the surface latent heat fluxes and hence the precipitation.

1 Introduction

The most dominant temporal mode in insolation and tropical precipitation is the 23 000-year precession cycle of the Earth . Both proxy (; ; ) and model (; ; ; ; ; ) based studies suggest that the intensities of monsoons have varied in proportion to insolation on orbital timescales. When changes in precession increase the insolation in the Northern Hemisphere, the zonal mean precipitation band shifts northward on account of the increase in thermal gradient between the two hemispheres (). This mechanism cannot explain the longitudinal changes in precipitation . The simulation of climate models shows that precipitation over land and oceans responds differently to precessional forcing (; ; ). This has been observed in the idealized as well as realistic precession experiments with climate models .

It is attributed to the land–sea contrast theory in the previous studies . The land warms more than the surrounding ocean due to its lower thermal inertia. Hence a low pressure develops over land and a monsoon circulation is established. The increase in insolation leads to deeper thermal lows over land, which enhance the onshore flow of moisture-laden winds. This leads to stronger ascent over land and an increase in precipitation. This thermal contrast, however, disappears after the onset of monsoon due to the cooling of land by precipitation and cloud cover. In fact, in good monsoon years, the land surface temperature is lower (Gadgil2018).

Some studies have used the changes in energy balance to understand the response of precipitation to precession . suggested that the net energy in the atmosphere over land and adjacent oceans changes due to precession. The atmosphere then acts to redistribute the excess energy, thereby setting up a land–ocean difference in precipitation. showed that the precipitation changes due to precession are related to the changes in the total column energy, which drives changes in vertical velocity. on the other hand has argued that the stability over oceans changes, whereas land regions respond by transporting the excess moist static energy. These are, however, generalizations for the entire tropics. The role of local processes and feedbacks might be important in driving regional changes in precipitation. Thus, individual regions need to be studied separately to understand the cause of the changes. For example, suggested that higher summer insolation leads to a migration of the near-surface moist static energy from the Bay of Bengal to India, before the onset of monsoon. They argued that hence the precipitation centroid shifts to India.

Moisture and MSE equations can be used separately to understand the dynamics of monsoon under different climate scenarios . In this work, we follow and demonstrate the advantage of combining the two equations. The resulting simple model attributes precipitation to energy fluxes and vertical stability of the atmosphere. This model, however, can only be used for regions where moisture and temperature gradients are weak. In this paper, we propose a modified version of the simple model which takes into account the horizontal gradients as well. We have used time-slice experiments in a high-resolution general circulation model (GCM), EC-Earth . This GCM was run in two orbital configurations which correspond to the extremes in precession (Fig. 1). The advantage of doing this is that the amplitude of the response is large, while the spatial pattern is similar to a simulation of realistic precession such as the Mid-Holocene (MH).

Figure 1The schematic diagram showing the orbital configuration of minimum precession (Pmin) and maximum precession (Pmax). In Pmin, summer solstice (SS) occurs at perihelion, while in Pmax, winter solstice (WS) coincides with the perihelion. AE and VE are the autumn and vernal equinoxes, respectively.

The paper is organized as follows. The next section describes the model and the experimental setup, and outlines the derivation of a simple model for the intertropical convergence zone (ITCZ). We have used this simple ITCZ model to understand the factors leading to the shift in precipitation between land and oceans, at the regional scale. The results are described in Sect. 3. It is followed by a discussion about the precipitation response to MH and obliquity forcing with the help of the simple ITCZ model.

2 Experimental design and analysis method

## 2.1 Climate model description

EC-Earth is a fully coupled ocean–atmosphere GCM . We have used model version 2.2. The Integrated Forecasting System (IFS) was the atmospheric component. The spectral resolution was T159 (roughly $\mathrm{1.125}{}^{\circ }×\mathrm{1.125}{}^{\circ }$) with 62 vertical levels. The convective scheme was used along with the H-TESSEL land surface scheme, including surface runoff. The Nucleus for European Modeling of the Ocean (NEMO, version 2) was the ocean component. The horizontal resolution was 1 with 42 vertical levels . NEMO includes sea-ice model LIM2. The OASIS3 coupler couples the ocean, sea ice, land, and atmosphere. EC-Earth performs well for the present day when compared to CMIP3 models in terms of climatology as well as inter-annual, spatial, and temporal variability .

## 2.2 Experimental designs

The two precession extremes, precession minima Pmin and precession maxima Pmax, correspond to summer solstice at perihelion and winter solstice at perihelion, respectively (Fig. 1). Table 1 shows the orbital configurations used. This leads to a stronger seasonal cycle in the Northern Hemisphere (NH) and a weaker seasonal cycle in the Southern Hemisphere (SH) in Pmin (Fig. 2). On the other hand, the seasonal cycle is weaker in the NH and stronger in the SH in Pmax.

Figure 2The difference in the incoming solar radiation at the top of atmosphere between Pmin and Pmax as a function of latitude and month.

Table 1The orbital configuration used for the extremes in precession, precession minima Pmin, precession maxima Pmax, and the pre-industrial. “e” represents eccentricity, δ is the tilt, and ω is the longitude of perihelion. The values of these have been chosen based on the extremes in the precession parameter $e*\mathrm{sin}\left(\mathit{\pi }+\mathit{\omega }\right)$ that have occurred in the last 1 Myr (Berger1978). Pre-industrial values are shown for comparison.

The model is run separately for each of the orbital configurations. The length of each simulation is 100 years, with the first 50 years being considered spinup. We have used the climatology of the last 50 years, for all our analysis. The orbital parameters remain constant throughout the simulation. All other boundary conditions (e.g., the solar constant, greenhouse gas concentrations, orography, ice sheets, vegetation) were kept constant at the pre-industrial levels. Vernal equinox has been fixed at 21 March, and the present-day calendar is used. Since the length of the season and the dates of equinoxes change along the precession cycle, the autumn equinoxes do not coincide. This is known as the “calendar effect”. It introduces some errors due to the phasing of insolation. We do not make any corrections in order to be consistent with previous studies. Further details about the experiments are provided in .

## 2.3 Diagnostic methodology

The Hadley cell is a thermally direct overturning circulation in the tropics. It takes energy away from the tropics and transports it towards the poles. The Hadley cell has a rising branch in the deep tropics and a descending branch in the extra-tropics. This leads to moisture convergence near the rising branch. The ITCZ coincides with the rising branch of the Hadley cell and is responsible for the zone of heaviest precipitation in the tropics. The characteristics of the ITCZ can be described by using the conservation equations for moist static energy (MSE) and moisture. Using this approach, proposed a simple model for ITCZ in terms of net energy input into the atmosphere and vertical stability. This is a diagnostic model that has been used to explain variations in rainfall due to global warming and the impacts of aerosols . In this section, we have discussed this simple model in detail. Equations (1) and (2) correspond to the conservation of MSE and moisture in a vertical column of the atmosphere. The first term in both the equations is horizontal divergence, with the second term being the vertical divergence of MSE and moisture fluxes, respectively.

The quantities on the right-hand side are the sum of all sources and sinks. Further details on the derivation of Eq. (1) can be found in . The time derivatives have been dropped in these equations because the climate is assumed to be in a steady state. The angle brackets (〈〉) indicate vertical integral.

$\begin{array}{}\text{(1)}& & 〈\mathrm{\nabla }\cdot m\mathbit{U}〉+〈\frac{\partial m\mathit{\omega }}{\partial p}〉={Q}_{\mathrm{div}},\text{(2)}& & 〈\mathrm{\nabla }\cdot q\mathbit{U}〉+〈\frac{\partial q\mathit{\omega }}{\partial p}〉=E-P,\text{(3)}& & 〈A〉=-\underset{{P}_{b}}{\overset{{P}_{t}}{\int }}A\frac{dp}{g},\end{array}$

where P is the precipitation rate (mm day−1). E is the evaporation rate (mm day−1). Qdiv is the total column energy, i.e., the sum of all the energy fluxes into the atmosphere (Eq. 4) (in mm day−1; taking the latent heat of vaporization as 2.26×106 J kg−1, we get 1 mm day${}^{-\mathrm{1}}=\mathrm{26.16}$ W m−2). Over land since the storage term is small, the sum of all the energy fluxes at the surface is small. Hence, Qdiv is mainly governed by the fluxes at the top of atmosphere (TOA). However, over oceans the contribution of surface fluxes is large. q is specific humidity (kg kg−1). m is moist static energy (J kg−1), which is the sum of internal energy, potential energy, and moist energy (${C}_{\mathrm{p}}T+gZ+{L}_{\mathrm{v}}q$). Pb is the pressure at the bottom of the atmospheric column (Pa). Pt is the pressure at the top of the atmospheric column (Pa). g is the acceleration due to gravity (m s−2). The full equation for Qdiv is

$\begin{array}{ll}{Q}_{\mathrm{div}}& =\underset{\text{bottom fluxes}}{\underbrace{\mathrm{LHF}+\mathrm{SHF}+\mathrm{Net}\mathrm{_}\mathrm{Sfc}\mathrm{_}\mathrm{Rad}}}\\ \text{(4)}& & +\underset{\text{TOA Fluxes}}{\underbrace{\mathrm{Net}\mathrm{_}\mathrm{TOA}\mathrm{_}\mathrm{LW}+\mathrm{Net}\mathrm{_}\mathrm{TOA}\mathrm{_}\mathrm{SW}}}.\end{array}$

where LHF is surface latent heat flux (mm day−1). SHF is surface sensible heat flux (mm day−1). Net_Sfc_Rad is net surface radiation (long wave + short wave) (mm day−1). Net_TOA_LW is net top of atmosphere longwave radiation (mm day−1). Net_TOA_SW is net top of atmosphere shortwave radiation (mm day−1). Clubbing all the radiation fluxes together into one quantity “Qrad”, we get

$\begin{array}{}\text{(5)}& {Q}_{\mathrm{div}}=\mathrm{LHF}+\mathrm{SHF}+{Q}_{\mathrm{rad}}.\end{array}$

Assuming ω=0 at the top as well as the surface leaves us with the horizontal terms only. The governing equations can be combined and simplified as

$\begin{array}{}\text{(6)}& & P-E=\frac{{Q}_{\mathrm{div}}}{\mathrm{GMS}},\text{(7)}& & \mathrm{GMS}=\frac{{m}_{\mathrm{1}}-{m}_{\mathrm{2}}}{{L}_{\mathrm{v}}\left({q}_{\mathrm{2}}-{q}_{\mathrm{1}}\right)},\text{(8)}& & {m}_{\mathrm{1}}=\frac{{\int }_{{P}_{\mathrm{t}}}^{{P}_{\mathrm{m}}}m\mathrm{\nabla }\cdot \mathbit{U}\mathrm{d}p/g}{{\int }_{{P}_{t}}^{{P}_{\mathrm{m}}}\mathrm{\nabla }\cdot \mathbit{U}\mathrm{d}p/g},\text{(9)}& & {m}_{\mathrm{2}}=\frac{{\int }_{{P}_{\mathrm{m}}}^{{P}_{b}}m\mathrm{\nabla }\cdot \mathbit{U}\mathrm{d}p/g}{{\int }_{{P}_{\mathrm{m}}}^{{P}_{b}}\mathrm{\nabla }\cdot \mathbit{U}\mathrm{d}p/g},\end{array}$

where GMS is the gross moist stability, as obtained by taking the ratio of Eqs. (2.11) and (2.12) from . m1 and m2 are, respectively, the total MSE in the upper (mid troposphere to top) and lower troposphere (surface to mid troposphere), normalized by the divergence of that layer. Thus, GMS is mainly a function of vertical profiles of MSE, and it provides a measure of vertical stratification of the atmosphere. Pm is pressure at the mid-troposphere level. Similarly, q1 and q2 represent the total moisture in the upper and lower troposphere, normalized by divergence. The mass convergence in the lower troposphere is the same as the mass divergence in the upper troposphere. The horizontal variations of temperature and moisture are assumed to be weak within the tropics. This implies that the horizontal advection of temperature and moisture is small. This simple model attributes the changes in PE to either the changes in total column energy or the vertical stability.

Figure 3The dependence of PE on (a) Qdiv and (b) GMS, for three regions: central India (15–25 N; 73–83 E), the Bay of Bengal (10–20 N; 85–95 E), and Africa (5–15 N; 20 W–0 E). The months Jun–Jul–Aug are taken separately. The blue and orange symbols correspond to Pmin and Pmax, respectively.

Figure 3a shows a scatter of PE as a function of Qdiv for the 3 summer months June, July, and August taken separately. The scatter is made for central India, the Bay of Bengal, and northern Africa for each of the precession extremes. We chose these three regions to highlight that neglecting the role of horizontal advection may not always be appropriate. The plot shows that the two are nearly linear, as indicated by the simple model (Eq. 6). The slight deviations from linearity are due to variations in GMS. As we go from Pmax to Pmin (low to high insolation in NH summer), both Qdiv and PE increase over central India and northern Africa (land regions). However, both these quantities decrease over the Bay of Bengal (oceanic region). The net energy input into the atmosphere and thus Qdiv is positive for all these regions during the summer.

We have shown in Fig. 3b a scatter of PE vs. GMS for the same regions. The figure shows that there is no definite relation between the two. Equation (6) suggests that all values for GMS should be positive since PE and Qdiv are both positive. There are, however, some points in the scatter where GMS is negative. This indicates that the assumption about the horizontal advection being small is not always valid. Hence, we modify the definition of GMS to include the horizontal advection terms.

By taking the ratio of Eqs. (1) and (2), after multiplying Eq. (2) by Lv (the latent heat of vaporization for water), we get

$\begin{array}{}\text{(10)}& & P-E=\frac{{Q}_{\mathrm{div}}}{\mathrm{TGMS}},\text{(11)}& & \mathrm{TGMS}=\frac{〈\mathrm{\nabla }\cdot m\mathbit{U}+\frac{\partial m\mathit{\omega }}{\partial p}〉}{-{L}_{\mathrm{v}}〈\mathrm{\nabla }\cdot q\mathbit{U}+\frac{\partial q\mathit{\omega }}{\partial p}〉},\end{array}$

where TGMS stands for total GMS (the term “total” indicates inclusion of all advection terms). TGMS is based on only one assumption, that the time derivatives of m and q are negligible. This is a good assumption for seasonal mean conditions. TGMS is particularly useful for smaller regions, where horizontal advection can be large. TGMS represents how efficiently an atmospheric column can diverge MSE per unit moisture converged into the column. TGMS is an extension of the concept of GMS, with horizontal advection terms included. This suggests that, along with the energy fluxes and vertical stratification of a column, the lateral transport of MSE and moisture determine the precipitation. A value of TGMS similar in magnitude to GMS indicates that the horizontal transport of MSE is negligible. A change in TGMS between two climates would suggest that the transport of MSE has changed. We have used the equivalence in Eqs. (1) and (2) to estimate TGMS. Since our goal is not to estimate the changes in PE but to diagnose the cause of these changes, there is no need to make an independent estimate of TGMS.

To quantify the relative contribution of Qdiv and TGMS to the changes in PE, we do the following. Writing Eq. (10) for Pmax,

$\begin{array}{}\text{(12)}& P-E=\frac{Q}{G},\end{array}$

where P, E, Q, and G are precipitation, evaporation, Qdiv, and TGMS, respectively. Considering Pmax as the reference case and Pmin as the perturbed case, we can write the following for Pmin:

$\begin{array}{}\text{(13)}& \left(P+\mathrm{\Delta }P\right)-\left(E+\mathrm{\Delta }E\right)=\frac{Q+\mathrm{\Delta }Q}{G+\mathrm{\Delta }G},\end{array}$

where Δ represents the perturbation from Pmax. Now dividing by PE, we get

$\begin{array}{}\text{(14)}& \mathrm{1}+\frac{\mathrm{\Delta }\left(P-E\right)}{\left(P-E\right)}=\frac{\mathrm{1}+\mathrm{\Delta }Q/Q}{\mathrm{1}+\mathrm{\Delta }G/G}.\end{array}$

This equation can further be modified as

$\begin{array}{}\text{(15)}& \underset{\text{Change in}\phantom{\rule{0.125em}{0ex}}P-E}{\underbrace{\mathrm{\Delta }\left(P-E\right)}}=\underset{\text{Contribution from}\phantom{\rule{0.125em}{0ex}}{Q}_{\mathrm{div}}}{\underbrace{\frac{\frac{\mathrm{\Delta }Q}{Q}}{\mathrm{1}+\frac{\mathrm{\Delta }G}{G}}\left(P-E\right)}}+\underset{\text{Contribution from TGMS}}{\underbrace{\frac{-\frac{\mathrm{\Delta }G}{G}}{\mathrm{1}+\frac{\mathrm{\Delta }G}{G}}\left(P-E\right)}}.\end{array}$
3 Results

In this section, we have explained the changes in PE between Pmin and Pmax in terms of Qdiv and TGMS. We start by giving an overview of the entire tropics and then we look at the South Asian monsoon in detail.

Figure 4The difference in precipitation (PminPmax) for all the tropical land and ocean taken separately (30 S–30 N).

## 3.1 Response of tropical precipitation to precession

### 3.1.1 Spatial patterns of the response

Figure 4 shows the difference in precipitation between Pmin and Pmax averaged over the tropical land and oceans separately. Precipitation change over the tropical land is out of phase with the changes in precipitation over the oceans. The amplitude of the change is higher over land than over oceans. Furthermore, over land, the change is of a higher magnitude during the boreal summer than the austral summer. This implies that the Northern Hemisphere monsoons are more sensitive to precession than the Southern Hemisphere monsoons. The vernal equinoxes during Pmin and Pmax occur on 21 March. Therefore, the difference in insolation between the two cases is very small during March. Hence the changes in land and ocean precipitation have a zero crossing during this month. Since we are interested in regions where there is moisture convergence, our analysis will focus on PE instead of precipitation.

Figure 5The difference in PE (a, c) and Qdiv (b, d). (a) and (b) are for the JJA mean and (c) and (d) are for the DJF mean.

In Fig. 5, the spatial patterns of the changes in PE and Qdiv are shown averaged over JJA (panels (a) and (b)) and DJF (panels (c) and (d)). First, we discuss the changes in precipitation during JJA. Most of the land regions in the Northern Hemisphere show an increase in PE. The African monsoon has strengthened substantially in Pmin, with an increase of about 10 mm day−1. PE has, in general, decreased over the oceans. However, there are many regions over the oceans (e.g., the Arabian Sea) where PE has increased. Hence, the amplitude of the changes in PE is small when averaged over all the tropical oceans (Fig. 4). The changes in Qdiv have a pattern similar to that of PE, with positive values over most of the land regions, and both positive and negative values over the oceanic regions. This is due to the direct relation between PE and Qdiv as suggested by the simple model. There are, however, some exceptions like the Arabian Sea and Africa. Qdiv has decreased over the Arabian Sea, but PE has increased. The regions of the largest increase in PE and Qdiv are not co-located over Africa. These are on account of the changes in TGMS.

During DJF, Pmin has less insolation (Fig. 2), and correspondingly a decrease in PE and Qdiv is seen over the land regions (Fig. 5c and d). Over oceans, there are regions of both positive and negative changes in PE during DJF as well. The magnitude of changes in Qdiv is on a similar order during JJA and DJF. However, the changes in PE are larger during JJA compared to DJF.

Figure 6The contribution of Qdiv and TGMS to the changes in PE. (a) is for the JJA mean and regions in the Northern Hemisphere, while (b) is for regions in the Southern Hemisphere and averaged over DJF. The blue bar is the change in PE, while pink and red bars represent the contribution from Qdiv and TGMS. The abbreviations used in (a) are N Land: northern tropics (land only), N Ocean: northern tropics (ocean only), CI: central India, BoB: the Bay of Bengal, SE Asia (Lnd): Southeast Asia (land only), SE Asia (Ocn): Southeast Asia (ocean only), N. Af: northern Africa, and AS: Arabian Sea, and in (b), S Land: southern tropics (land only), S Ocean: southern tropics (ocean only), S. At: South Atlantic, S. Af: South Africa, and N. Aus: northern Australia. The coordinates of these regions are provided in Table 2.

### 3.1.2 Dominant factors determining the response of tropics

In this section, we look at the various terms in Eq. (15) for different regions of the tropics (Fig. 6). Panels (a) and (b) are for the Northern Hemisphere and Southern Hemisphere, respectively. The analysis was done for the summer months of the respective hemispheres (JJA for the Northern Hemisphere and DJF for the Southern Hemisphere). The blue bar represents the changes in PE, whereas the light red and dark red bars are contributions from Qdiv and TGMS. Qdiv explains most of the changes in PE when all the land regions in the northern tropics are taken together. This need not be true in smaller regions. For example, TGMS contributes most to the changes in PE over Africa. Because PE has a different sign over various oceanic regions, the change in PE, averaged over all the tropical oceans, is small. The contributions from Qdiv and TGMS are in opposite directions, thus almost canceling each other out. The contribution from TGMS is, however, slightly higher. The Arabian Sea shows an increase in PE, due to a change in TGMS. The decrease in PE over the Bay of Bengal is, however, mainly due to changes in Qdiv, with the changes in TGMS being small.

In the southern tropics the dominant contribution is from changes in Qdiv over land and changes in TGMS over oceans. In the cases of South Africa and Brazil changes in TGMS and Qdiv make an equal contribution. TGMS drives most of the changes in PE over northern Australia and the South Atlantic. Figure 6 highlights the fact that the mechanisms for the changes in precipitation are region specific. Hence, each region must be studied separately to understand the physical mechanism that caused the changes in PE. Both the Indian land mass and the Bay of Bengal are part of the Indian monsoon system, yet they demonstrate a different response to the precessional forcing. Hence, we discuss this asymmetric response in detail in the following subsection. Such an asymmetry also exists within the East Asian monsoon, which has been discussed in a separate subsection.

Figure 7The seasonal cycle of near-surface equivalent potential temperature (θe) for India and the Bay of Bengal for the (a) Pmin configuration and (b) Pmax.

## 3.2 Explaining the response of the Indian monsoon to precession

suggested that different response of the Indian land mass and the Bay of Bengal is due to migration of near-surface equivalent potential temperature from the Bay of Bengal over to India. This is because the rate of increase in insolation is higher in the high insolation (similar to Pmin) experiment. This causes the equivalent potential temperature θe to rise rapidly over India. It is known that the location of ITCZ coincides with that of the surface energy maxima . Hence ITCZ migrates over India quickly and remains there. However, EC-Earth simulates higher near-surface equivalent potential temperature θe over the Bay of Bengal, in both Pmin and Pmax (Fig. 7). In this section, we propose an alternate mechanism for the different response of the Indian land mass and the Bay of Bengal to the changes in precession.

We have shown earlier that in Pmin there is an increase in Qdiv over the Indian land mass and a decrease over the Bay of Bengal with respect to Pmax (Figs. 3 and 6). Here we examine the factors that caused the changes in Qdiv. Splitting Qdiv into its component fluxes (Eq. 5) will help us to determine which flux contributed the most. Figure 8 is a spatial map of the differences in PE, Qdiv and its component fluxes Qrad, LHF, and SHF. Qdiv has a good spatial coherence with PE, over most of the regions except the Arabian Sea. As was discussed earlier, this is due to the changes in TGMS, which is able to counter the effect of reduced Qdiv. PE has decreased along the southern parts of the Western Ghats but has increased in the northern parts of the Western Ghats. Qrad bears a resemblance to PE. This suggests that radiative feedbacks from clouds are present. Changes in Qrad are not large enough to counter the decrease in latent heat flux (LHF) over the Arabian Sea and the Bay of Bengal. Thus, the decrease in LHF over these regions reduces Qdiv there. In fact, Qdiv and LHF have similar spatial patterns over the oceanic regions. The changes in sensible heat flux (SHF) are small in most places.

Figure 8The JJA mean difference (PminPmax) in (a) PE, (b) Qdiv (sum of energy fluxes at the top and bottom of the atmosphere), (c) Qrad (sum of all radiative fluxes at the top and bottom of the atmosphere), (d) latent heat flux, and (e) sensible heat flux. The boxes shown in (a) and (b) are the regions chosen for this study: central India (15–25 N; 73–83 E), the Bay of Bengal (10–20 N; 85–95 E), and Southeast Asia (0–25 N; 100–125 E). (f) shows the decomposition of Qdiv into radiative, latent, and sensible heat fluxes for the two regions.

We take two regions: one over central India and the other over the Bay of Bengal, to identify the flux which contributes most to the changes in Qdiv. These regions are outlined with black boxes in Fig. 8a. The changes in the three components of Qdiv over these two regions are depicted in the bar chart (Fig. 8f). It shows the dominance of the radiative terms over India, and LHF over the Bay of Bengal, respectively.

LHF is a function of surface wind speed, sea surface temperature (SST), and near-surface relative humidity. LHF increases with an increase in SST and wind speed. SST has increased over the Bay of Bengal and southern Arabian Sea by about 2 C (Supplement Fig. S1). Thus, it cannot explain the decrease in LHF. Hence, we look at the changes in wind speed (Fig. 9). Fig. 9a and b show the mean winds at 850 hPa. The shading indicates wind speed and the streamlines show the direction of flow. The axis of the low-level jet (LLJ) has shifted to the north, and this has led to a decrease in winds over the Bay of Bengal. Due to LLJ, deep oceanic water upwells along the coast of Somalia. This leads to cooler SSTs over the western parts of the Arabian Sea. Since in Pmin, LLJ has shifted further north, the region of upwelling also shifts north, thus leading to cooler SSTs in the western coast of the northern Arabian Sea (Fig. S1). Hence, LHF over the Arabian Sea decreases due to weaker winds in the southern parts and colder SST in the northern parts.

Figure 9The JJA mean wind speed (850 hPa) in shading for (a) Pmin and (b) Pmax, with streamlines of the wind vector field superimposed. The difference of the winds between Pmin and Pmax is shown in (c); 40 E longitude has a convergence at the Equator and cyclonic circulation over the Middle East. An anti-cyclonic circulation exists in the Southern Hemisphere over Madagascar. This is similar to the response of the atmosphere to equatorial plus off-equatorial heating (Gill1980).

The shift in LLJ leads to lesser moisture flux along the southern part of the Western Ghats. Hence, PE decreases there. At the same time, the LLJ brings more moisture into the northern parts of the Western Ghats, leading to increase in PE. The shift of the LLJ can be seen more clearly in Fig. 9c, where the difference in winds between Pmin and Pmax is shown. Along the Equator, there exists an anomalous low-level easterly over the Indian Ocean. This meets an anomalous westerly from over equatorial Africa, at around 40 E. This indicates low-level convergence. Furthermore, on the same meridian, there exists a cyclonic circulation to the north (over the Middle East) and an anti-cyclonic circulation to the south (over Madagascar). This resembles the response of the winds to the heating of an atmospheric column as shown by Gill (1980).

Gill (1980) proposed a simple shallow water model on an equatorial β plane to elucidate the role of latent heating on surface winds. In order to represent convective heating due to latent heat release, he introduced mass divergence in the atmospheric column. When this model was forced with “heating” over a region at the Equator and another region to the north of the Equator, it produced a Kelvin wave and a mixed Rossby-gravity wave. The Kelvin wave leads to an anomalous low-level easterly and an anomalous low-level westerly along the Equator. The easterly is to the east of the heat source and the westerly to the west of the heat source. These anomalous winds thus lead to low-level convergence at the Equator, near the region of the heat source. The mixed Rossby-gravity wave has a cyclonic circulation to the north of the Equator and an anti-cyclonic circulation to the south of the Equator. The wind response of EC-Earth hence, suggests that the wind patterns over the Indian subcontinent, are driven by atmospheric heating near the Equator and off-Equator. Examining Fig. 5a shows that the heat sources correspond to convective heating of the column due to increased precipitation over the West Equatorial Indian Ocean (WEIO) and over the Middle East (particularly the Red Sea). There are, however, some important differences between EC-Earth and the Gill model. EC-Earth is a full GCM with non-zero mean background winds, whereas Gill model is linearized with respect to zero mean background winds. Thus, the EC-Earth’s response includes non-linear terms as well.

Figure 10The difference between Pmin and Pmax in (a), Qdiv and (b), PE along with streamlines of change in the wind, for the month of May. Note that the large increase in Qdiv over Africa leads to an early onset of African monsoon. Thus, influencing the winds over the equatorial Indian Ocean.

To summarize, the decrease in Qdiv over the Bay of Bengal is due to lower wind speeds. The winds decrease because of convective heating over west equatorial Indian Ocean and the Red Sea. The convection over the Red Sea is an extension of the African monsoon. Hence, we examine the factors which lead to an increase in precipitation over these regions. The prevailing conditions in the pre-monsoon month of May, leads to enhanced convection over these regions, later in the summer. Figure 10a and b, show the difference in Qdiv and PE for the month of May. Figure 10b shows changes in PE in shading and the streamlines represent the changes in the wind direction. Qdiv is higher over Africa, and this causes early onset of the African monsoon (Fig. 10b) and changes the low-level winds along the eastern coast of Africa. The SST along the eastern coast depends on the coastal upwelling. The changes in winds thus reduce upwelling and increase SST. This enhances convection over the West Equatorial Indian Ocean, further leading to low-level convergence. This positive feedback is responsible for the convective heating that persists through the summer months. As the season advances from May onwards, the African monsoon propagates northward. The region of convection over the eastern side of Africa moves over to the Red Sea. This becomes the off-equatorial heat source.

## 3.3 Factors determining the response of Southeast Asian monsoon to precession

showed that the Southeast Asian monsoon and the Northeast Asian monsoon are out of phase owing to the El Niño-like SST pattern in Pmin. Here we are addressing the differences in the precipitation changes over Southeast Asia (land) and the adjacent ocean. The domain for Southeast Asia is shown in Fig. 8b. Based on the analysis using Eq. (15), we find that the increase (decrease) in PE over the land (ocean) grids is mainly due to the increase (decrease) in Qdiv (Fig. 6). Even though Qdiv is dominant, the contribution of TGMS is higher over Southeast Asia (oceanic regions) when compared to the Bay of Bengal. Once again decomposing Qdiv into its component fluxes suggests a similar mechanism that leads to the Indian Bay of Bengal redistribution of precipitation (Fig. 8f). The increase in insolation leads to an increase in Qdiv over the Southeast Asian land, whereas a decrease in LHF over the oceanic regions leads to a decrease in Qdiv. The convective heating over WEIO and the Red Sea leads to reduced winds, and hence decreased LHF in the northwestern Pacific.

4 Discussion

In this section, we have discussed the similarities between the sets of idealized experiments (Pmin, Pmax) vs. (Mid-Holocene (MH), Pre-Industrial (PI)). The MH and PI experiments were conducted with the same model EC-Earth, the details of which are available in . The difference in solar forcing between MH and PI is similar to that between Pmin and Pmax, albeit with a smaller amplitude (Fig. S2). Moreover, MH has an obliquity 0.66 higher than PI, and hence it contributes little to the total forcing (Fig. S2b). Previous research with models has shown that the climate response to precession is independent of obliquity . The climate of MH is therefore mainly driven by precession. The peak in the insolation difference between MH and PI is delayed by a month with respect to the insolation difference between Pmin and Pmax. Hence the largest precipitation changes in MH occur about a month later than in Pmin (Figs. 4 and S3). Therefore, we consider Jul–Aug–Sep averages for MH. The land-ocean shift in precipitation in MH is qualitatively explained by changes in Qdiv (Fig. S4). Particularly, the displacement of precipitation from the Bay of Bengal to India is due to the same mechanism that drives these changes in Pmin (Figs. S5, S6, and S7). The SE Asian monsoon also exhibits a land-ocean shift in rainfall. This is due to radiative heating over land as well as the ocean. This suggests that the cloud radiative feedbacks are stronger for the SE Asian monsoon. The changes in LHF are, however, due to the same reason as in Pmin. We also performed the analysis for a set of obliquity experiments Tmax and Tmin, corresponding to the maximum and minimum tilt, with eccentricity set to zero (). The tropical precipitation shows a land-ocean shift in precipitation, but the amplitude of change is small compared to the precession experiments (Fig. S9). The mechanisms leading to this shift are different for obliquity and precession (Figs. S10 and S11).

Models with different levels of complexities: QTCM , Quasi-geostrophic model EC-Bilt , GCM with slab ocean and finally the fully coupled model EC-Earth have all shown a shift in precipitation between land and ocean, when subjected to the precessional forcing. However, there are no proxies for precipitation over oceans to verify this. Since the climate over islands is influenced by the surrounding oceans, proxies obtained from islands can be thought of as a representation of climate over the surrounding ocean. A speleothem chronology from the Baratang cave in the Andaman Islands in this regard, represents precipitation over the Bay of Bengal. This chronology goes back to 4000 years before present and shows a long-term decreasing trend in precipitation as we move back in time. The time period corresponding to 4 ka being closer to MH has higher summer insolation and proxies over Indian continent register an increase in precipitation . This suggests that the GCMs and observations indicate the response of Indian land mass is different from the response in the Bay of Bengal.

5 Summary and conclusions

Using a simple model for ITCZ, we have interpreted the response of a high resolution fully coupled model EC-Earth to precession. The changes in precipitation can be attributed to either the changes in total energy fluxes going into the column (Qdiv) or the changes in vertical stability of the atmosphere (TGMS). We have included the horizontal advection terms in the calculation of TGMS, which were originally assumed to be small . This allows us to use the simple ITCZ model for relatively smaller domains, where horizontal advection terms can be large. TGMS represents the total transport of the MSE. In places where the horizontal transport is weak, TGMS is the same as GMS. Changes in precession provide an initial forcing. The final response of the precipitation is due to this initial forcing and the consequent feedbacks. These feedbacks are in the form of changes in surface energy fluxes and changes in stability of the atmosphere. In agreement with , we find that precipitation changes between precession extremes over the whole tropics are, due to changes in Qdiv over land and due to TGMS over the ocean. This generalization is, however, not valid for smaller regions. Within the domain of the South Asian monsoon, insolation drives changes in Qdiv over the land, whereas latent heat fluxes contribute most over the oceans. Particularly, the decrease in LHF over the Bay of Bengal and the northwestern Pacific is associated with the weakening of the low-level westerlies over these regions. These changes in westerlies are driven by convective heating of the atmospheric column over the western equatorial Indian Ocean and the Middle East. There are, however, regions where the changes in TGMS is the main cause of the changes in precipitation (e.g., Africa and the Arabian Sea). We have demonstrated that the simple ITCZ model can be used to explain the precipitation response for any orbital configuration (e.g., MH, maximum and minimum obliquity experiments).

Data availability
Data availability.

Explanation for data not being publicly available: the data are being used for publication, and hence cannot be made public as yet. The data are available upon request.

Supplement
Supplement.

Author contributions
Author contributions.

CJ, JS and AC analysed and interpreted the GCM output. JHCB designed and ran the experiments. CJ wrote the manuscript with input from all authors. All authors reviewed the manuscript.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

We thank A. Nikumbh for useful comments. The authors acknowledge support from the Centre for Excellence in the Divecha Centre for Climate Change (DCCC). This work was partially funded by DST India.

Review statement
Review statement.

This paper was edited by Qiuzhen Yin and reviewed by two anonymous referees.

References

Balsamo, G., Beljaars, A., Scipal, K., Viterbo, P., van den Hurk, B., Hirschi, M., and Betts, A. K.: A revised hydrology for the ECMWF model: Verification from field site to terrestrial water storage and impact in the Integrated Forecast System, J. Hydrometeorol., 10, 623–643, 2009. a

Battisti, D., Ding, Q., and Roe, G.: Coherent pan-Asian climatic and isotopic response to orbital forcing of tropical insolation, J. Geophys. Res.-Atmos., 119, 11997–12020, 2014. a, b, c, d

Bechtold, P., Köhler, M., Jung, T., Doblas-Reyes, F., Leutbecher, M., Rodwell, M. J., Vitart, F., and Balsamo, G.: Advances in simulating atmospheric variability with the ECMWF model: From synoptic to decadal time-scales, Q. J. Roy. Meteor. Soc., 134, 1337–1351, 2008. a

Berger, A.: Long-term variations of daily insolation and Quaternary climatic changes, J. Atmos. Sci., 35, 2362–2367, 1978. a, b

Boos, W. R. and Kuang, Z.: Mechanisms of poleward propagating, intraseasonal convective anomalies in cloud system–resolving models, J. Atmos. Sci., 67, 3673–3691, 2010. a

Bordoni, S. and Schneider, T.: Monsoons as eddy-mediated regime transitions of the tropical overturning circulation, Nat. Geosci., 1, 515–519, 2008. a

Bosmans, J. H. C., Drijfhout, S. S., Tuenter, E., Lourens, L. J., Hilgen, F. J., and Weber, S. L.: Monsoonal response to mid-holocene orbital forcing in a high resolution GCM, Clim. Past, 8, 723–740, https://doi.org/10.5194/cp-8-723-2012, 2012. a, b

Bosmans, J., Drijfhout, S., Tuenter, E., Hilgen, F., and Lourens, L.: Response of the North African summer monsoon to precession and obliquity forcings in the EC-Earth GCM, Clim. Dynam., 44, 279–297, 2015. a, b, c, d

Bosmans, J., Erb, M., Dolan, A., Drijfhout, S., Tuenter, E., Hilgen, F., Edge, D., Pope, J., and Lourens, L.: Response of the Asian summer monsoons to idealized precession and obliquity forcing in a set of GCMs, Quaternary Sci. Rev., 188, 121–135, 2018. a

Braconnot, P., Marzin, C., Grégoire, L., Mosquet, E., and Marti, O.: Monsoon response to changes in Earth's orbital parameters: comparisons between simulations of the Eemian and of the Holocene, Clim. Past, 4, 281–294, https://doi.org/10.5194/cp-4-281-2008, 2008. a, b, c

Caley, T., Roche, D. M., and Renssen, H.: Orbital Asian summer monsoon dynamics revealed using an isotope-enabled global climate model, Nat. Commun., 5, 5371, https://doi.org/10.1038/ncomms6371, 2014. a

Chamales, K. A.: The effects of orbital precession on tropical precipitation, MS thesis, Dept. of Meteorology and Physical Oceanography, University of Miami, 63 pp., 2014. a, b, c, d

Chou, C. and Neelin, J. D.: Mechanisms of global warming impacts on regional tropical precipitation, J. Climate, 17, 2688–2701, 2004. a

Chou, C., Neelin, J. D., Lohmann, U., and Feichter, J.: Local and remote impacts of aerosol climate forcing on tropical precipitation, J. Climate, 18, 4621–4636, 2005. a

Chou, C., Neelin, J. D., Tu, J.-Y., and Chen, C.-T.: Regional tropical precipitation change mechanisms in ECHAM4/OPYC3 under global warming, J. Climate, 19, 4207–4223, 2006. a

Clement, A. C., Hall, A., and Broccoli, A.: The importance of precessional signals in the tropical climate, Clim. Dynam., 22, 327–341, 2004. a

Cruz Jr., F. W., Burns, S. J., Karmann, I., Sharp, W. D., Vuille, M., Cardoso, A. O., Ferrari, J. A., Dias, P. L. S., and Viana Jr., O.: Insolation-driven changes in atmospheric circulation over the past 116,000 years in subtropical Brazil, Nature, 434, 63–66, 2005. a

Donohoe, A., Marshall, J., Ferreira, D., and Mcgee, D.: The relationship between ITCZ location and cross-equatorial atmospheric heat transport: From the seasonal cycle to the Last Glacial Maximum, J. Climate, 26, 3597–3618, 2013. a

Gadgil, S.: The monsoon system: Land–sea breeze or the ITCZ?, J. Earth Syst. Sci., 127, 1–29, https://doi.org/10.1007/s12040-017-0916-x, 2018. a

Gill, A.: Some simple solutions for heat-induced tropical circulation, Q. J. Roy. Meteor. Soc., 106, 447–462, 1980. a, b, c

Hazeleger, W., Severijns, C., Semmler, T., Ştefǎnescu, S., Yang, S., Wang, X., Wyser, K., Dutra, E., Baldasano, J. M., Bintanja, R., Bougeault, P., Caballero, R., Ekman, A. M. L., Christensen, J. H., van den Hurk, B., Jimenez, P., Jones, C., Kållberg, P., Koenigk, T., McGrath, R., Miranda, P., van Noije, T., Palmer, T., Parodi, J. A., Schmith, T., Selten, F., Storelvmo, T., Sterl, A., Tapamo, H., Vancoppenolle, M., Viterbo, P., and Willén, U.: EC-Earth: a seamless earth-system prediction approach in action, B. Am. Meteorol. Soc., 91, 1357–1363, 2010. a, b

Hazeleger, W., Wang, X., Severijns, C., Ştefănescu, S., Bintanja, R., Sterl, A., Wyser, K., Semmler, T., Yang, S., van den Hurk, B., van Noije, T., van der Linden, E., and van der Wiel, K.: EC-Earth V2. 2: description and validation of a new seamless earth system prediction model, Clim. Dynam., 39, 2611–2629, 2012. a, b

Hsu, Y.-H., Chou, C., and Wei, K.-Y.: Land–ocean asymmetry of tropical precipitation changes in the mid-Holocene, J. Climate, 23, 4133–4151, 2010. a, b, c

Kang, S. M., Held, I. M., Frierson, D. M. W., and Zhao, M.: The response of the ITCZ to extratropical thermal forcing: Idealized slab-ocean experiments with a GCM, J. Climate, 21, 3521–3532, 2008. a

Kathayat, G., Cheng, H., Sinha, A., Yi, L., Li, X., Zhang, H., Li, H., Ning, Y., and Edwards, R. L.: The Indian monsoon variability and civilization changes in the Indian subcontinent, Science Advances, 3, e1701296, https://doi.org/10.1126/sciadv.1701296, 2017. a

Kutzbach, J., Liu, X., Liu, Z., and Chen, G.: Simulation of the evolutionary response of global summer monsoons to orbital forcing over the past 280,000 years, Clim. Dynam., 30, 567–579, 2008. a

Kutzbach, J. E.: Monsoon climate of the early Holocene: climate experiment with the earth's orbital parameters for 9000 years ago, Science, 214, 59–61, 1981. a, b

Laskar, A. H., Yadava, M., Ramesh, R., Polyak, V., and Asmerom, Y.: A 4 kyr stalagmite oxygen isotopic record of the past Indian Summer Monsoon in the Andaman Islands, Geochem. Geophy. Geosy., 14, 3555–3566, 2013. a

Madec, G.: NEMO, the ocean engine, Note du Pole de modelisation, Institut Pierre-Simon Laplace (IPSL), France, No. 27 ISSN No. 1288–1619, 2008. a

Merlis, T. M., Schneider, T., Bordoni, S., and Eisenman, I.: Hadley circulation response to orbital precession. Part II: Subtropical continent, J. Climate, 26, 754–771, 2013. a

Mohtadi, M., Prange, M., and Steinke, S.: Palaeoclimatic insights into forcing and response of monsoon rainfall, Nature, 533, 191–199, 2016. a

Neelin, J. D. and Held, I. M.: Modeling tropical convergence based on the moist static energy budget, Mon. Weather Rev., 115, 3–12, 1987. a, b, c, d, e

Patnaik, R., Gupta, A. K., Naidu, P. D., Yadav, R., Bhattacharyya, A., and Kumar, M.: Indian monsoon variability at different time scales: marine and terrestrial proxy records, Proceedings of the Indian National Science Academy, 78, 535–547, 2012. a

Pokras, E. M. and Mix, A. C.: Earth's precession cycle and Quaternary climatic change in tropical Africa, Nature, 326, 486–487, 1987. a

Privé, N. C. and Plumb, R. A.: Monsoon dynamics with interactive forcing. Part I: Axisymmetric studies, J. Atmos. Sci., 64, 1417–1430, 2007a. a

Privé, N. C. and Plumb, R. A.: Monsoon dynamics with interactive forcing. Part II: Impact of eddies and asymmetric geometries, J. Atmos. Sci., 64, 1431–1442, 2007b. a

Ramesh, R.: High resolution Holocene monsoon records from different proxies: An assessment of their consistency, Curr. Sci. India, 81, 1432–1436, 2001. a

Schneider, T., Bischoff, T., and Haug, G. H.: Migrations and dynamics of the intertropical convergence zone, Nature, 513, 45–53, 2014. a

Shi, Z.: Response of Asian summer monsoon duration to orbital forcing under glacial and interglacial conditions: Implication for precipitation variability in geological records, Quaternary Sci. Rev., 139, 30–42, 2016. a

Shi, Z., Liu, X., and Cheng, X.: Anti-phased response of northern and southern East Asian summer precipitation to ENSO modulation of orbital forcing, Quaternary Sci. Rev., 40, 30–38, 2012. a

Sterl, A., Bintanja, R., Brodeau, L., Gleeson, E., Koenigk, T., Schmith, T., Semmler, T., Severijns, C., Wyser, K., and Yang, S.: A look at the ocean in the EC-Earth climate model, Clim. Dynam., 39, 2631–2657, 2012. a

Sun, Y., Zhou, T., Ramstein, G., Contoux, C., and Zhang, Z.: Drivers and mechanisms for enhanced summer monsoon precipitation over East Asia during the mid-Pliocene in the IPSL-CM5A, Clim. Dynam., 46, 1437–1457, 2016. a

Sun, Y., Ramstein, G., Li, L. Z., Contoux, C., Tan, N., and Zhou, T.: Quantifying East Asian summer monsoon dynamics in the ECP4. 5 scenario with reference to the mid-Piacenzian warm period, Geophys. Res. Lett., 45, 12–523, 2018.  a

Tuenter, E., Weber, S., Hilgen, F., and Lourens, L.: The response of the African summer monsoon to remote and local forcing due to precession and obliquity, Global Planet. Change, 36, 219–235, 2003. a, b, c

Tuenter, E., Weber, S., Hilgen, F., Lourens, L., and Ganopolski, A.: Simulation of climate phase lags in response to precession and obliquity forcing and the role of vegetation, Clim. Dynam., 24, 279–295, 2005. a

Valcke, S. and Morel, T.: OASIS and PALM, the CERFACS couplers, Tech. rep., CERFACS, 2006. a

Wang, X., Auler, A. S., Edwards, R., Cheng, H., Ito, E., Wang, Y., Kong, X., and Solheid, M.: Millennial-scale precipitation changes in southern Brazil over the past 90,000 years, Geophys. Res. Lett., 34, L23701, https://doi.org/10.1029/2007GL031149, 2007. a

Wang, Y., Cheng, H., Edwards, R. L., Kong, X., Shao, X., Chen, S., Wu, J., Jiang, X., Wang, X., and An, Z.: Millennial-and orbital-scale changes in the East Asian monsoon over the past 224,000 years, Nature, 451, 1090–1093, 2008. a

Weber, S. and Tuenter, E.: The impact of varying ice sheets and greenhouse gases on the intensity and timing of boreal summer monsoons, Quaternary Sci. Rev., 30, 469–479, 2011. a

Zhang, E., Wang, Y., Sun, W., and Shen, J.: Holocene Asian monsoon evolution revealed by a pollen record from an alpine lake on the southeastern margin of the Qinghai-Tibetan Plateau, China, Clim. Past, 12, 415–427, https://doi.org/10.5194/cp-12-415-2016, 2016. a

Zhao, Y. and Harrison, S.: Mid-Holocene monsoons: a multi-model analysis of the inter-hemispheric differences in the responses to orbital forcing and ocean feedbacks, Clim. Dynam., 39, 1457–1487, 2012. a, b