Southern hemisphere monsoonal system during superinterglacial stages: MIS5e, MIS11c and MIS31

The current study investigates changes in Austral Summer Monsoon based on numerical experiments conducted with the coupled ICTP-CGCM model. The interannual variability and intensity of the monsoonal system have been analyzed from vorticity indices and air-sea interaction in Africa, Australia and South America. We focus on interglacial stages MIS5e (127 ka), MIS11c (409 ka) and MIS31 (1072 ka). Results show orbitally-driven decreased summer precipitation and slightly shifted monsoon onset and demise with respect to present day conditions. Sensitivity experiments indicate that monsoons are forced not only by the dominant effect of insolation, but also by remote teleconnections, such as the equatorial Atlantic and Pacific ocean basins. During those interglacial stages, cooling occurs in the Southern Hemisphere whereas Northern Hemisphere substantially warms. This induces meridional displacement of oceanic subtropical high pressure systems and the equatorial convergence zone. Regionally, these mechanisms contribute to droughts over the Amazon and northeastern Brazil, northern Australia and southern Africa. Monsoonal rainfall shows different responses to precessional forcing, as well as the relationship between the monsoon and Niño 3.4 differs among the interglacial stages. Results also indicate a weaker influence of the equatorial Pacific Ocean on the Austral summer monsoon for the MIS31 interglacial stage as compared to current climate conditions across Africa and Australia. On the other hand, South America monsoon is strongly influenced by Niño 3.4 and tropical Atlantic.


Introduction
Previous studies have demonstrated that the strength of global hydrological cycle increases in a warmer climate (Held and Soden 2006), but additional scrutiny is needed to verify hydroclimatic changes under precessional interhemispheric seesaw temperature anomalies, with Boreal warming and Austral cooling, as typical of interglacial stages (Justino et al. 2017;Melles et al. 2012;Braconnot et al. 2008;Lisiecki et al. 2005). Despite the large effort reported in the literature to investigate past climatic features, such as precipitation and temperature during distant interglacial stages, evaluation of the spatio-temporal patterns of Southern Hemisphere (SH) monsoon related to changes in external factors is still needed. The SH monsoon characteristics are more noticeable during the austral summer, in which the inversion of the zonal wind component from east to west are highly correlated to significant increases in precipitation rates (Yim et al. 2014;Geen et al. 2020). The seasonal zonal wind shifts in the lower troposphere is a response to the land-ocean thermal contrast which causes rainy summers and dry winters. Investigation of the SH monsoon system under very different climate conditions are very useful to foresee adaptive measures to cope with future changes in the anthropocene. Past changes in the radiative forcing and decreased net energy input have also been responsible 1 3 for the weakening of SH monsoons in the mid-Holocene (D'Agostino et al. 2020).
Simulations of past climate is an excellent tool to provide information on the impact of external forcing in the SH tropical monsoons. Interglacial climates such as occurred during the Marine Isotope Stages (MIS) MIS5e, MIS11c, and MIS31, at approximately 127 ka, 409 ka and 1072 ka (ka = thousand years ago), are characterized by Northern Hemispheric increases in oceanic and terrestrial temperatures. The MIS5e, also known as the Eemian, experienced sudden changes in atmosphere and ocean dynamics in response to strong insolation in the Northern Hemisphere (NH), but weaker in the SH (Yin and Berger 2012;Berger and Loutre 1991;Siccha et al. 2015; Bard et al. 1990). Additional evidences indicate that the northward migration of the Intertropical Convergence Zone (ITCZ) during the Eemian, favor more intense NH summer monsoon as compared with the present climate (Montoya et al. 2000).
The MIS11c ( ∼ 409 ka) also experienced global increase in temperatures in response to enhanced Boreal insolation. However, during this interglacial, the Asian summer monsoon (ASM), as revealed by data from the Yongxing cave, China, does not differ remarkably from monsoon characteristics in the Late Holocene (Zhao et al. 2019). It should be noted that the intensity of the MIS11c warming, is regionally dependent (Yin and Berger 2012;Oppo et al. 1998;Rousseau et al. 1992), rising questions on how MIS11c warming or slight cooling may have induced changes in local precipitation. The MIS31 ( ∼ 1072 ka) is considered the warmest marine isotope stage interval that contributed to substantial melting of NH polar glaciers compared to the present day (Lisiecki et al. 2005;Melles et al. 2012). Modeling studies have demonstrated that orbitally-driven enhanced seasonality affected the global atmospheric and oceanic circulation, which includes the El Niño-Southern Oscillation (ENSO) magnitude and periodicity (Clement et al. 2001;Timmermann et al. 2007;Justino et al. 2017Justino et al. , 2019. Insofar as monsoonal patterns are concerned, it is found that the Indian monsoon was enhanced but no correlation with the ENSO is identified in the MIS31 climate, in contrast with presentday conditions (Trenberth 2011;Yim et al. 2014). In addition, the MIS31 climate shows weak correlation between the equatorial Pacific SST and precipitation indices for Asia, Australia, and Western North Pacific monsoons during this interglacial period (Justino et al. 2017).
Thus, past interglacial temperature contrasts between land and oceanic regions, and the out-of-phase SST pattern induced by the precessional cycle, between both hemispheres might have triggered important changes in seasonal precipitation. Intercomparison among SH monsoon is evaluated based on a series of coupled climate model simulations conducted for different interglacial stages, including all forcing factors such as changes in orbital parameters and greenhouse gases characteristics. The study describes in Sect. 2 the coupled model specifications and experimental modelling setup. Section 3 focuses on results under different perspectives: (1) the control climate is compared to ERA5, (2) changes in the radiative components and temperature are analyzed; and (3) the impact of those changes on atmospheric dynamics responsible for anomalous magnitude, onset and demise of monsoonal precipitation during the interglacials, are shown. Finally, Sect. 4 presents the conclusions in the light of previous investigations and paleo-proxies perspectives.

ICTP-CGCM and modeling experimental design
The present study is based on simulations conducted with the Global Coupled Atmosphere-Ocean Circulation Model developed at the International Centre for Theoretical Physics (ICTP-CGCM). The GCM model consists of the global atmospheric climate model "SPEEDY" (version 41) coupled to the Nucleus for European Modeling of the Ocean (NEMO) model (Madec et al. 1998;Madec 2008), and uses the OASIS3 as coupler (Valcke 2013). In computational terms, SPEEDY proofs to be effective in reproducing the main characteristics of the climate system of tropical and extratropical latitudes (Molteni 2003;Kucharski et al. 2006;Justino et al. 2021). The atmospheric component runs at eight vertical levels in T30 horizontal resolution (Kucharski et al. 2016). The oceanic model NEMO runs at 31 vertical levels with thicknesses ranging from 10 to 5000 m of the ocean floor, with 16 levels in the first 200 m. The current version uses a tri-polar ORCA2 configuration with a horizontal grid resolution of 2 • and a tropical refinement of 0.5 • (Madec 2012).
To investigate the differences between the current climate and the interglacial stages in a systematical way, a series of sensitivity experiments are performed. The first experiment is the control climate (CTRL). It employs an atmospheric CO 2 concentration of 380 ppm and present-day albedo, vegetation, orbital forcing, and topographic boundary conditions as described in Justino et al. (2017Justino et al. ( , 2019. Three experiments have been carried out for 500 model years to characterize the interglacial intervals. Moreover, the orbital parameters have been set originally according to Berger (1978) and also adopted by Coletti et al. (2015). The atmospheric greenhouse gas concentration is extracted from Bereiter et al. (2015); Lüthi et al. (2008). The study centered the interglacial time at three fixed periods, 127 ka (MIS5e), 409 ka (MIS11c) and 1080 ka (MIS31), however, the orbital forcing is characteristics of the vernal equinox as in the control climate (Table 1). In fact, the position of the vernal Equinox will define when the fall and spring season occur and the length of the season. According to Timm et al. (2008) the season varies in length by up to ±6 days. According to their findings, studies on monsoon activity should be compared with fixed angular seasons. For process studies and climate sensitivity investigation the use of fixed angular seasons is imperative. For all experiments the first 400 years of the simulation are not included in the analyses of results. Discussion herein is based on the last 100 years of the 500 years model run.
Several studies have discussed procedure to verify the climatic equilibrium in modeling experiments (Foley 1995;Peltier and Solheim 2004;Otto-Bliesner et al. 2021). In the present case the focus is on atmospheric processes on seasonal timescales, thus, it is reasonable to assume that surface oceanic conditions play the lead role as boundary conditions, with respect to the impact of deep ocean temperatures that vary in much longer timescales (Valdes et al. 2021). It has to be noted that setting up just a fixed time slide do not cover all aspects of the interglacial stages. Notwithstanding, our selected dates are associated with the maxima temperature values throughout those intervals according to Coletti et al. (2015). Verification of the CTRL climate in reproducing present day conditions is based on comparison with the ERA5 reanalysis from 1979 to 2020 (Hersbach 2016). Comparison is conducted for precipitation, zonal wind component (u), sea surface temperatures (SST) and Sea Level Pressure (SLP).

CTRL radiative components and temperature
There are no substantial changes of the atmospheric composition among the interglacials (Table 1), but reduced concentrations favour slightly cooling conditions with respect to CTRL, due to lower atmospheric CO 2 concentration and resulted weaker greenhouse capacity of the atmosphere. These features contrast with modifications of the Earth's orbit during the interglacial epochs, in particular the precessional cycle. Figure 1 shows zonally averaged radiative components at the top of atmosphere (TOA) based on ERA5 and the CTRL simulation.
Overall, the coupled model is able to represent the energy balance at the TOA and surface. The ICTP-CGCM underestimates ERA5 values by about 5% at the TOA (Fig. 1a). Good match between ERA5 and the model is found in the zonal distribution, but the presence of the ITCZ in the equatorial region leads to more pronounced differences across 10 • S and 10 • N. Atmospheric models running under lower resolution show drawbacks in reproducing meso-scale convective system in the tropics leading to modified radiative balance (Tomita et al. 2005;Satoh et al. 2019). In the polar regions and extratropics as expected (Fig. 1a), more energy is emitted to space via longwave (LW) radiation than is provided by the shortwave (SW) component. In general grounds, the modeled radiative budget (Rn) is well represented, that in fact, is responsable for the large-scale atmospheric and oceanic circulations. It should be emphasized that estimates of the radiative balance has been for many years a challenging task associated with improvements of coupled models (Wild 2020).
Differences of short wave radiation at the top of the atmosphere during the MIS5e, MIS11c and MIS31, with respect to the CTRL climate ( Fig. 1b-d), show the large influence of modified orbital parameters from June to November in the SH tropics, with an increase of 60 Wm −2 in the MIS5e, and 20 Wm −2 in MIS11c ( Fig. 1b-c). The most significant increase is, however, observed in NH where the shortwave radiation is larger than the CTRL values in summer. The MIS31 is among the interglacial periods the most affected by orbital changes from March to October. From November to March, the SH experiences reduction in the amount of incoming energy, which is the dominant period of the Austral monsoonal system (Fig. 1b-d). This may indicate that the attribution of warmer planetary or hemispheric climate in past interglacials, with respect to today's climate, in particular based on paleoreconstructions analyses, should be taken with care due to enhanced seasonal variability and regional characteristics. Figure 2 shows that modified orbital parameters and subsequent changes of incident solar radiation, result in anomalous near surface air temperatures. Differences between the interglacials and current climate across the continents are characterized by cooling from November to February (NDJF, Fig. 2a-c-e). This gradually turns to warmer conditions in Boreal spring with largest anomalies in summer, which are reproduced by the annual anomalous pattern  Fig. 2b-d-f). The SH summer monsoons, which are dominant in NDJF months, can be severely affected by surface thermal changes not only causing a reduction in the amount of rainfall, but also by inducing distinct onset and demise periods. Surface condition are determinant for SH monsoons because soil moisture and evapotranspiration play an important role at early stages of the monsoon (D'Agostino et al. 2020; Collini et al. 2008). The substantial warming delivered over the continents as shown by annual time averaged temperatures ( Fig. 2b-f), induces a northward migration of the ITCZ, as well as a weakening of the southern trade winds (Justino et al. 2019). The modified orbital forcing leads to reduced radiative balance in the Southern Hemisphere (SH) and associated warming across the Northern Hemisphere (NH). This pattern induces an inter-hemispheric thermal gradient that shift northward the ITCZ with respect to its regular position. This is also related to the anomalous pattern of SLP. Indeed, during the interglacials stages the SH (NH) experiences higher (lower) SLP with respect to the CTRL climate. This intensified middle-latitudes high pressure centers, modify the strength of the trade winds, prevailing southwesterly and weakening northeasterly, driving the northward migration of the ITCZ. It is not speculative to argue that these anomalous features will be associated with modifications in the thermal ocean-land characteristics, impacting the monsoonal system.

Annual distribution of precipitation and SH monsoonal characteristics
To identify the SH monsoons spatial domain, the amplitude of the first harmonic of precipitation has been computed ( Fig. 3a, b). The first harmonic characterizes regions where the annual cycle is dominant explaining more than 80% of the variance. This fits nicely with monsoonal features. Harmonic analyses also known as Fourier transformation, are computed on a series of sine and cosine functions which may be used to verify differences in the precipitation magnitude between the rainy and dry seasons, for instance. Figure 3b shows the first harmonic amplitude of precipitation, based on a 100-year CTRL simulation. The ICTP-CGCM is capable to reasonable locate the larger amplitude of precipitation over the South African Monsoon System (SAF), the Australian Monsoon System (AUS) and the South American Monsoon System (SAM) regions (Fig. 3a). The CTRL experiment (Fig. 3b), however, overestimates the 1st harmonic amplitudes by 5 mm/day in particular across  (Fig. 3a, b). Drawbacks in simulating tropical dynamics in coupled models, arise due to limitations in reproducing ocean-atmosphere interaction due to their lower resolution, as compared to processes in the extratropics. This occurs becausec onvection, turbulent heat fluxes and diabatic processes are more difficult to be properly simulated (Flato et al. 2014;Neelin et al. 1994;Marengo et al. 2012).
Furthermore, it should be mentioned that the CTRL climate discussed here shows differences from ERA5, in part because time intervals used are not the same. The CTRL is based on a 100 years simulation forced with initial conditions correspondent to 1950. Certainly, these different intervals drive slightly distinct climates, also due to the ICTP-CGCM interannual and interdecadal variability, which may differ with respect to ERA5 conditions computed from 1979 to 2020.
The CTRL climate delivers higher precipitation as reproduced by areal averages over SAMS and AUSMS and lower values are found for SAFMS (Fig. 3). Among the monsoons, higher daily precipitation amount is delivered for the AUSMS ( ≈ 10 mm/day), whereas lower values are observed in the SAFMS ( ≈ 6 mm/day). This is highlighted from November to March. Regionally, 70% of annual rainfall occurs during the NDJF months. ERA5 shows along the SAMS, SAFMS and AUSMS domains, magnitudes by about 7, 9 and 8 mm/day (Fig. 3c-e).
Correlation analyses based on monsoon indices as represented by the zonal circulation and the areal averaged rainfall show reasonable correspondence between ERA5 (b-d-f) is the same as (a-c-e) but for mean annual differences. Dotted regions are statistically significant at the 95% confidence level and the CTRL run (Table 2). Values up to 0.86 and 0.74 are found for AUSMS and SAFMS, but much lower correlation is delivered for SAMS (0.37). This indicates that insofar as the SAMS is concerned, the vorticity index in the CTRL run does not match the precipitation temporal characteristics related to the monsoon. The presence of the Andes, the recurrent cold fronts and the Amazonian contribution to the SAMS through evapotranspiration, are complex features that may weak the correlation between the vorticity index and SAMS precipitation in the ICTP-CGCM. It may be argue that nearby oceanic processes exert an important role to define the SAM. These processes involve SST changes in the southwestern Atlantic Ocean, which impact in two-ways interaction the strength of the South Atlantic Convergence Zone (SACZ), by modifying the maritime water vapor transport as well as the heat flux onto continental regions (Pezzi et al. 2022;Jorgetti et al. 2014). The SST and oceanic conditions impact on the monsoonal system are discussed in more detail in the following sections.

Monsoonal precipitation
Changes in insolation during the interglacials lead to drastically modified monsoonal precipitation amounts (Fig. 4).
With respect to CTRL, results show that increased precipitation occurs in the dry season (June, July and August) in particular for SAFMS and SAMS (Fig. 4a, c). These features are enhanced in the MIS11c and MIS5e, and may represent a slightly shift in the rainy period over the two monsoons domains. In AUSMS, larger changes are depicted from October to February, with a drop in rainfall (Fig. 4e), therefrom an increase of the rainfall is found until April. The MIS31 delivers remarkable changes of precipitation, with drier conditions in the SAMS, with precipitation as low as 60% in comparison to CTRL (Fig. 4f). Reduction is also found for the other SH systems, the weakening of the SAFMS precipitation occurs during the summer months, whereas for the AUSMS drier conditions are found from June to December (Fig. 4d, e). Changes in surface conditions such as shown by evapotranspiration (ET, Fig. 4g-i), in line with those depicted by precipitation in the SAFMS, indicate the importance of surface processes and precipitation, in particular soil moisture. Analyses for the other two SH monsoon do not demonstrate a close relationship between ET and precipitation. However, reduced precipitation over the SH monsoon regions, in turn, also leads to less water available for evapotranspiration. In the AUSMS, enhanced monsoon rainfall is accompanied by reduced ET, which lead to assume that the AUSMS during these interglacial, experiences large contribution from water vapor transport originating from outside the monsoonal domain.
The SAMS pattern shows a time lag relationship, in which the reduction of ET leads by two months negative anomalies of precipitation in the MIS11c and MIS5e. The . d-f show differences of precipitation ( Prec ) between the MIS5e, MIS11c, and MIS31 with respect to CTRL conditions. g-i are the same as d-f but for evapotranspiration (mm/ day, Evapo.) delay in the rainy season can be triggered mainly by the low rates of evapotranspiration in the Amazon and Cerrado by the end of the dry season. The seasonal rainfall reduction related to evapotranspiration, is also forced by the weakening mechanism of the SAMS. In the case of MIS31, changes in ET influence in lower degree modification of precipitation, with respect to the other interglacials (Fig. 4g-i).
Changes in surface conditions such as shown by evapotranspiration (ET, Fig. 4g-i) are in line with those depicted by precipitation in the SAFMS, indicating the importance of interaction between surface processes and precipitation, in particular soil moisture. Analyses for the other two SH monsoon do not demonstrate a close relationship between ET and precipitation. However, reduced precipitation leads to lower soil moisture and reduced evapotranspiration. In the AUSMS, enhanced monsoonal rainfall is accompanied by reduced ET, which may indicate that the AUSMS during these interglacial, experiences large contribution from water vapor transport originating from outside regions of the monsoonal domain. The pattern over the SAMS shows a time lag relationship, in which lowest ET leads by two months negative anomalies of precipitation in the MIS11c and MIS5e. The delay in the rainy season can be triggered mainly by low rates of evapotranspiration in the Amazon and Cerrado, at the end of the dry season. In the case of MIS31, changes in ET are not closely related to precipitation, with respect to the other interglacials (Fig. 4g-i).
The onset and demise of the rainy season related to monsoonal characteristics is also evaluated by the Hovmöller diagram (Fig. 5), which shows the inversion of the zonal wind; eastward winds in the dry season and westward winds for the rainy season. In general, changes in atmospheric circulation via the inversion of the zonal wind are well simulated by the ICTP-CGCM (Fig. 5a-c). Similar pattern is depicted for the SAFMS and AUSMS by indicating an east-west wind shift during the summer months (Fig. 5a, b). However, different pattern emerges for the SAMS where wind transition is not clearly seen. In the SAMS wind features are longitudinal and confined to the eastbound of the monsoon between 60 • -30 • W (Fig. 5c). Despite reasonably reproducing the monsoon annual cycle, the ICTP-CGCM does not relate the wind features to precipitation in the SAMS domain. This is associated with the topographic barrier imposed by the Andes, which channels throughout the year the easterly Atlantic flow, hampering the seasonal shift, as discussed by Nogués-Paegle et al. (2002).
The MIS5e and MIS11c deliver a semi-annual component that is most evident in the SAFMS and AUSMS but is absent in the MIS31, which exhibits a well defined annual cycle (Fig. 5j, k). Indeed, stronger annual cycle in the MIS31 has been found by Justino et al. (2019). Weaker interannual ENSO may act in line to enhance the orbital forcing (Carré et al. 2021). The SAMS is different because it experiences large changes related to the implementation of the orbital forcing (Fig. 5f, i, j). For this monsoon, the MIS5e Fig. 6 Correlation coefficient between regional summer monsoons precipitation indices and SST for SAF (a-d), AUS (e-h) and SAM (i-l). m-p display the correlation between the ENSO and SST anomalies. Also shown is the 95% significant areas (dotted grids) and 11c differ substantially from the CRTL with increased zonal wind shear, and subsequently increased convergence over 45 • W, from December to March (Fig. 5j, i). Moreover, it is found drastically modified circulation from August to November with respect to CTRL. Strong westerly flow is evident in the MIS5e and MIS11c climates from 60-30 • W, which results from intensified sub-tropical Atlantic high.
In order to evaluate in more detail mechanisms responsible for interglacial monsoonal changes, it is shown in Fig. 6, the correlation pattern between the SAFMS, AUSMS and SAMS vorticity indices and SST changes. It highlights that similarities between the SAFMS and AUSMS under CTRL conditions, do not propagate through interglacials stages (Fig. 6a-h). The MIS5e exhibits a pattern, in particular across the Pacific Ocean, that is between CTRL and the MIS11c and MIS31. For instances, both monsoon seem to be well correlated with Pacific and Indian Ocean SST in the CTRL climate, however, in MIS11c and MIS31, the SAFMS is highly related to Indian Ocean only (Fig. 6c, d). In the AUSMS case, changes over the warm pool Pacific region are dominant insofar as SST are concerned because weak correlation are found over the eastern Equatorial Pacific region (Fig. 6g, h).
As demonstrated by previous discussion, the SAMS ( Fig. 6i-l) differs substantially from other monsoons. The Atlantic Ocean exerts an important role which is dependent upon the interglacial investigated. Indeed, under CTRL conditions higher positive correlations are found between precipitation and warmer Atlantic SST, under La Niña-like background (Fig. 6i). Turning to MIS5e, the presence of the Tropical Atlantic variability (TAV, (Cabos et al. 2019) is enhanced with respect to CTRL, with a well defined dipole (Fig. 6j). The SST correlation pattern may indicate a southward position of the ITCZ, due to northern Atlantic cooling, but it has also to be noted that other forcing such as the meridional pressure gradient infleunces the ITCZ position. Positive correlation also in the southwestern Atlantic is other feature that indicates intensification of the South Atlantic Convergence Zone (SACZ). The relationship between SST and SAMS precipitation for the MIS11c and MIS31 demonstrated that the positive correlation is weaker in the southern Atlantic in the MIS11c, whereas an overall anticorrelation is found in the MIS31 (Fig. 6l). In general, the MIS31 exhibits contrasting values with respect to the CTRL (Fig. 6i, l), insofar as the SAM is concerned.
The eastern tropical Pacific exhibits a strong relationship with continental rainfall, particularly due to ENSO characteristics. ENSO events have been defined based on Empirical Orthogonal Function (EOF) of monthly SST anomalies over the tropical Pacific. El Niño and La Niña are the first principal component (PC) time series over the equatorial Pacific domain (20 N-20 S). Table 3 indicates that during El Niño events, monsoonal precipitation of the three analyzed sectors is reduced with larger changes for AUSM in both ERA5 and in CTRL. These significant correlations, also reveal that ENSO can drive interannual variability over the summer monsoons for the MIS5e and MIS11c, with exception of the AUSSM sector. On the other hand, the AUSSM is negatively correlated with ENSO in MIS31 at a 95% confidence level. Although, the spatio-temporal variability patterns of rainfall are different for each interglacial, it depends on the intensity of individual ENSO events. The correlations between SH continental monsoon rainfall and TAV shows that under CTRL conditions most correlations are not significant. Turning to the interglacial intervals, the analyses demonstrate weaker correlations with respect to ENSO, but significant over the SAMS in the MIS5e. Interesting is that correlations between TAV and SAMS (− 0.33) are higher than correlations based on ENSO (− 0.29) ( Table 3).

Spatial patterns of Precipitation and SLP
Figures 7 and 8 show differences in precipitation and SLP (interglacial minus control) overlaid with significant correlation (dotted), between the regional monsoon vorticity index; and precipitation and SLP. Analyses of SLP can be further used to determine convergence and divergent flows which may increase or suppress precipitation.
During the DJFM there is a reduction in rainfall in southern/southeast Africa and much of South America for MIS5e and MIS11c, with values between − 2 and − 5 mm/ day (7a-f, with respect to CTRL. Precipitation increases are found in the Indonesian Archipelago and southern Australia by up to 4 mm/day (Fig. 7b-e). The dominant effect of the orbital forcing in all 3 interglacials is highlighted as reduced precipitation across most of continental regions. In the MIS31, reduced rainfall in Africa and South America show values between − 2 and − 6 mm/day (Fig. 7g-i). Rainfall increases in northern Amazon and semiarid Northeast Brazil (3 mm/day), but there is no evidence that this is linked to monsoonal changes. It should be noted that precipitation in Australia are less affected by the implementation of the orbital forcing. Most anomalies range between ± 2 mm/day, and are not statistically significant in MIS31.
Analyses of differences of SLP may clarify causes of precipitation changes. The SLP anomalies between the interglacials experiments (MIS5e, MIS11c e MIS31) and CTRL, exhibit statistically significant increased pressure in line with precipitation reduction over SAFMS, AUSMS and SAMS regions (Fig. 8). Differences of MIS11c and CTRL show slight reduction in SLP across southern Australia and south-central Brazil (Fig. 8e, f). Those low pressure anomalies induce precipitation (Fig. 7e) in Australia by enhancing convection and upward motion. It has to be noticed that this drop in SLP during the MIS11c Fig. 7 Austral monsoon precipitation anomalies between MIS5e minus CTRL (a-c); MIS11c minus CTRL (d-f) and MIS31 minus CTRL (g-i). Dots indicate significant correlation between the regional monsoonal vorticity index and precipitation in southern Brazil/South America, theoretically would increase the frequency of subtropical frontal systems, passing onto the SAMS region, but the high pressure anomaly associated with the subtropical Atlantic high, leads to a blocking situation hampering the subtropical system to move northwards. Thus, leading to dryness in most of Brazil during the interglacial. In the MIS31, most striking feature is related to cooling in western Pacific by up to − 3 • C, in particular over the AUSMS domain. This results in positive SLP anomalies characteristics of permanent El Niño conditions (Figs. 2f and 8h), reduced AUSMS monsoonal precipitation (Fig. 7h).

Intercomparison between modeled precipitation and proxy data
Paleoclimate simulations are able to simulate Earth's past climate, but comparison with reconstruction are still necessary to reveal model-proxies differences, in regions where models may struggle to represent particular atmospheric/ oceanic features. Thus, indicating climate mechanisms that should be addressed in more detail (Braconnot et al. 2012). Justino et al. (2017) provides an evaluation of modeled MIS31 surface temperature and proxies, in which demonstrated that the ICTP-CGCM is able to reproduce reconstructions, differing only by ± 1 • C, but larger values up to 3 • C are found in the extratropics. These large differences have been attributed to reduced NH sea-ice cover, and therefore higher SSTs. The three interglacials are warmer (colder) in the NH (SH) than the CTRL climate (Table 4). It is also worth mentioning that the ICTP-CGCM model simulates lower temperatures than the reconstructions primarily in the Austral extratropical region for the MIS5e and MIS11c. These two interglacials show warmer temperatures in the NH (2-6 • C). The paleoproxies also exhibits caveats because in some cases it may reproduce long term changes that are dictated by a particular seasonal strength, such as demonstrated by Liu et al. (2014). This argumentation is not used to justify differences between model and reconstructions, but usefully serves to indicate that reconstructions represent local conditions that can be smoother in a grid box representation.
Hydroclimatic reconstructions from MIS5e and MIS11c are shown in Table 5 and Fig. 9, across the 5 • N-33 • S domain. This allows to identify the precipitation and evaporation relationship, and therefore atmospheric humidity in both time slices, with respect to CTRL climate. By evaluating proxies for individual latitudinal bands shows that between 5 • N-5 • S, drier conditions prevailed during NDJF months in the MIS5e, with exception of African sites in reconstruction and model results (red circles in Fig. 9). These conditions are present not only during the monsoonal season bit it is distributed throughout the year (Fig. 9a-c).
From 6 • S to 15 • S most reconstructions are marked by increased humidity (blue circles in Fig. 9), differing from simulated values, that during the monsoonal season show dryer conditions, in particular in central Africa and Australia. Agreement is found over Indonesian Archipelago, east Africa and the Andes. By comparing reconstructions and Fig. 9 Map showing precipitation differences between MIS5e and CRTL simulations overlaid by proxy records (wet-blue dots and dry-red dots) during the MIS5e, with respect the present day; for NDJF (a), MAMJ (b) and JASO (c) model results, good agreement between these data is found for the MAMJ and JASO months, indicating that reconstructions primary reproduce enhanced precipitation in Australia, Andes and eastern Africa, from March to October during the MIS5e (Fig. 9b, c). Southward to 16 • S most proxies indicate wetter conditions in Australia and southern Africa. In contrast, the simulated climate is characterized across land surfaces in southern Africa by drier conditions (Fig. 9a,  c). The MIS5e climate shows that the orbital forcing leads to decreased austral summer precipitation, due to anomalous high pressure center which may explain the simulated drought in Southern Africa (Fig. 8a, c).
The ICTP-CGCM simulation matches nicely the reconstructions over Australia and western South America indicating more humid conditions during the MIS5e with respect to CTRL (Fig. 9). Some conditions persist along the months, such as the weakening of the ITCZ, that is supported by both modelled precipitation and reconstructions. In fact, based on reconstruction this depends on the ocean basin. In the Atlantic, it is very clear the most proxies are characterized by lower precipitation during those interglacials. The paterrn is different across the Pacific, where precipitation increases in the equatorial belt along the South America coast. Little information is available for MIS11c regarding moisture. Proxies from Africa and South America allow us to infer that during MIS11c conditions were wetter than during the current climate. For MIS31 there is scarce information which hampers to derive conclusions from comparison between model outputs in this study. But Justino et al. (2019) argue that the link between the ENSO and the Australian monsoon is weakened with respect to the CTRL characteristics, and indicates that changes in AUSM during the MIS31 are more closely connected to hemispherical features than to the equatorial Pacific.

Concluding remarks
Limitation about our conclusions arises because results are based on single time-slice simulation which may reproduce the interglacial climates characteristics of the summer peak. It is well known that this approach does not capture the temporal variability of those interglacials in response to slow changes of the orbital forcing. For the time being, these caveats related to our technological capacity that does not support in a suitable time millennial scales climate implementations, the experiments are, however, very useful to investigate seasonal peak conditions allowing to understand how the climate system is affected by strengthened radiative forcing, leading to NH warming and/or SH cooling. This study indicates that the Austral summer monsoons were severely affected during the MIS5e, MIS11c and the MIS31 interglacials due to changes in orbital insolation, which further modifies the SSTs and global teleconnection patterns. Indeed, changes in the precessional cycle triggered significant cooling across Southern Oceans. In the SAMS region, for example, MIS5e and MIS11c show a decrease in the transport of moisture from the Amazon to southeastern South America. The contribution of continental dry air mass favored the decrease of summer monsoon rainfall, and the lower SST during these interglacials, leads to a northward displacement of the SASM, as well as of the SACZ.
Changes of precipitation in summer for the MIS5e and MIS11c stages are dominated by enhanced rainfall in the SH winter and early spring for SAFM, and winter/autumm for AUSM. Reduced precipitation is related to colder conditions during the SH summer, in particular over the continents. Correlation analyses indicate that large-scale processes, such as ENSO are weaker when compared to control run, in exception of increased significant correlation values over the SAMS. Although ICTP-CGCM performs well in reproducing the current climate, the lack of paleoreconstruction data makes it difficult to analyze interglacial periods prior to MIS5e.
Comparison between palaeoreconstructions and model simulations also revealed a good agreement for the analyzed interglacials. It is worth noting that the palaeoproxies deliver the main response to a particular forcing, but does not indicate the preferential season of occurrence. The use of paleo-proxies is always subject of uncertainties because it reproduces local conditions, where climate models deliver values based on areal averages. In addition, results do not show transient climate as may be reproduced by proxies. However, the current study provided an attempt to compare the modeled interglacial climates with paleoreconstruction, but discrepancies can unfortunately be found, because those intervals, which occurred under very particular conditions and high seasonality, may have been dominated by vegetation patterns different than today, modifying the global evaporation rates and the hydrological cycle.
Climate simulations aiming at reproducing the monsoonal system have indicated that enhanced Northern Hemisphere warming can modify the temporal and spatial pattern of Southern Hemisphere precipitation. In a palaeoclimate perspective, this may be similar to conditions during Dansgaard-Oeschger events by means of the two hemispheres seesaw Stocker and Johnsen (2003). To some extent, recent global warming (GW) also shows warmer Boreal climate as delivered by interglacials stages. Therefore, future GW may also result in modified SH monsoonal precipitation with potential to impact the environment and vulnerable population. The understanding of past events could bring important insight on the expected behavior of the climate system and monsoonal precipitation in the future.