Dipole patterns in tropical precipitation were pervasive across landmasses throughout Marine Isotope Stage 5

Most of Earth’s rain falls in the tropics, often in highly seasonal monsoon rains, which are thought to be coupled to the inter-hemispheric migrations of the Inter-Tropical Convergence Zone in response to the seasonal cycle of insolation. Yet characterization of tropical rainfall behaviour in the geologic past is poor. Here we combine new and existing hydroclimate records from six large-scale tropical regions with fully independent model-based rainfall reconstructions across the last interval of sustained warmth and ensuing climate cooling between 130 to 70 thousand years ago (Marine Isotope Stage 5). Our data-model approach reveals large-scale heterogeneous rainfall patterns in response to changes in climate. We note pervasive dipole-like tropical precipitation patterns, as well as different loci of precipitation throughout Marine Isotope Stage 5 than recorded in the Holocene. These rainfall patterns cannot be solely attributed to meridional shifts in the Inter-Tropical Convergence Zone. Tropical rainfall patterns varied throughout Marine Isotope Stage 5, beyond the influence of meridional shifts in the Inter-Tropical Convergence Zone, suggest analyses of tropical hydroclimate proxy records and model-based precipitation reconstructions.

T he hydrological cycle is expected to intensify as Earth warms, owing to increased moisture availability (following the Clausius-Clapeyron equation) 1 . Of crucial importance is an accurate projection of precipitation change in tropical regions (30°N-30°S), where dense human populations live with intimate socioeconomic links with climate. The majority of Earth's rain falls in the tropics and in some regions it is highly seasonal (e.g., the monsoons), being coupled to the interhemispheric migrations of the Inter-Tropical Convergence Zone (ITCZ) in response to the seasonal cycle of insolation. Therefore, the tropical precipitation forms an integral part of the global climate system. Yet, despite both its regional and global significance, the response of rainfall to future climate has been proven challenging to predict 2 and there remains a poor characterization of rainfall in the geologic past. Reconstructions of past rainfall are limited by an unequal spatial and temporal coverage of proxy records, augmented by chronological uncertainties. Moreover, the principal challenge in gaining an accurate quantification of past rainfall amount resides in the capacity of geochemical climate proxies to faithfully record a single, specific climate variable of interest 3 (e.g., local rainfall amount). Paleoclimate records inferring past rainfall variability have nevertheless typically converged on attributing changes in precipitation to fluctuations in the meridional position of the ITCZ, with the ITCZ position across the Late Pleistocene primarily being regulated by insolation, Northern Hemisphere ice-sheet size and/or alterations of the oceanic overturning circulation [4][5][6] . A notable example is the identified antiphase trend between speleothem oxygen isotope (δ 18 O) records from China with those in South America during periods of Northern Hemisphere cooling and warming 7 . Attributing north-south ITCZ shifts have seemed a sufficient explanation for describing tropical paleoclimate variability given that the ITCZ is sensitive to energy fluxes, and thus ice-sheet build-up in the northern hemisphere should be associated with a parallel decrease in tropical precipitation in the same hemisphere 8 . More recently, however, the view of large-scale uniform shifts of the ITCZ position for determining past rainfall responses has been cast into doubt by modelling studies, suggesting variable regional 9,10 and minimal observed meridional ITCZ movements despite applying significant climate forcing 11 . Regional contractions and expansions of the tropical precipitation belt have been inferred for the late Holocene, based on some proxy studies, as better depicting the nature of past rainfall changes rather than a simple uniform meridional shift of the ITCZ in response to past climate forcings 12,13 . This suggested regional heterogeneity in tropical precipitation highlights the rapidly increasing spatial coverage of proxy records. This has resultant implications on understanding of the "Global Monsoon", a concept that links the regional monsoon systems (seasonal rainfall affecting large parts of Asia, Africa, Australia and the Americas) in an attempt to convey their collective sensitivity to the annual cycle of solar radiation and atmospheric overturning circulation 14 . Thus, application of this concept in past monsoon rainfall reconstructions may be inadequate due to its inability to encapsulate the potential diversity in the response of each regional monsoon system to climate forcing 15 .
There have been meaningful and successful attempts within the paleoclimate community to compile global and large-scale regional records of sea surface temperatures (SSTs) across various time periods in the past [16][17][18] , although equivalent studies of past rainfall have largely remained limited 19,20 . To depict the evolving characteristics of geographic variations of rainfall to past changes in climate, a greater geographic coverage of proxy records and more integrated approaches are needed 19 . Here we investigate transient tropical past hydroclimate variability in context to the "Global Monsoon" during the warm interval of Marine Isotope Stage (MIS) 5 (130-70 ka) by combining newly generated proxy records of the Indian Summer Monsoon (ISM) from the core ISM rainfall region of the Bay of Bengal (Methods and Fig. 1) with existing paleoclimate records that are inferred to capture hydroclimate (Methods). We employ the broad term hydroclimate to encompass the various nuances in paleoclimate approaches to reconstruction of "wet" or "dry" past environments. To evaluate the robustness of these signals, Fig. 2 compares the proxy records with fully independent reconstructions of precipitation that are derived from PaleoPGEM 21 , a spatiotemporal emulator of the intermediate-complexity atmosphere-ocean GCM PLASIM-GENIE 22 . Although investigations to capture a specific time snapshot using high-resolution equilibrium simulations 19 are useful, our approach, which explores past hydroclimate using a time series of hydroclimate emulations at intermediate resolution with uncertainty quantification, permits a more holistic assessment of the behaviour of the tropical rain belt in response to evolving climate forcing.

Results and discussion
We have employed a simple approach of allocating "wetter" and "drier" conditions relative to the mean signal (across the study interval, 70-130 ka), to present a diverse suite of tropical multiproxy records in a unified manner for comparison ( Fig. 2) (Methods). These tropical multi-proxy records, inferred to capture past hydroclimate, include various caveats in their ability to accurately record past hydroclimate associated with measurement (instrument/user) and calibration errors in addition to the interpretation of the reconstructed signal, of which the latter is by far the largest source of uncertainty. Furthermore, these records are predominantly confined to capturing a terrestrial rainfall signal owing to the nature of sample locations (Fig. 1). However, the collations of the available proxy-based hydroclimate records show a relative consistency among themselves and with the model-based (yellow line, Fig. 2) average precipitation record across MIS 5 (Fig. 2). This evidence for regional wetting and  (Fig. 2), red circles represent proxy records used only in the anomaly map (Fig. 3) and stars represent proxy records generated in this study (red, U1446 and orange, U1448) (see Methods and Supplementary Data 1 for detail and references). Precipitation data available from https://psl.noaa.gov/data/gridded/data.cmap.html. Map produced using Panoply (source: https://www.giss.nasa.gov/tools/panoply/). drying across the warm sub-stages (MIS 5e, c and a) and cooler sub-stages (MIS 5d, b) lends support to the ability of the collated proxy records to capture a representative signal of hydroclimate across MIS 5, in spite of the range of complications associated with each respective proxy record. This is further supported by good coherence found between the collation of our proxy records with a recently published record of the Dole effect (ΔDE*), proposed as a measure of the low-latitude hydrological cycle 23 ( Supplementary Fig. 1). However, deviations between proxy records nevertheless exist at both the orbital and millennial scales (Fig. 2), highlighting the influence of local climate processes, differing sensitivities and time-transgressive responses to boundary conditions 24 , and complication with the proxy signal itself (Methods). Regions 3 (Fig. 2e) and 6 ( Fig. 2f) exhibit clear divergences among the compiled proxy records within each respective region manifested by contemporaneous wetting and drying episodes across MIS 5 (Fig. 2e, f). Region 3 encompasses the two Asian Monsoon subsystems, the ISM and East Asian Summer Monsoon (EASM), in addition to the maritime continental regions of the Indo-Pacific Warm Pool (IPWP), which also overlaps in Region 6 and further includes the Australian-Indonesian Monsoon. Discrepancies in temporal coupling between the ISM and EASM have been documented at both orbital and millennial scales [24][25][26][27] , and are attributed to the differing sensitivities of these two subsystems to boundary conditions and climatological differences 28 . The observed discrepancies may also be attributed to the interpretation of the monsoon rainfall signal reconstructed from speleothems in the regions, which remains contentious (Methods). However, the spatiotemporal emulation of Region 3 exhibits a lack of correlation at the local scale (mean correlation of 0.19 between individual grid cells and the regional mean in PaleoPGEM; Supplementary  Fig. 2), highlighting that the observed diversity in response to climate forcing in this region captured by the proxy records is The emulation-derived precipitation is shown in yellow with the lighter yellow shading representing the 95% confidence limit. The black line represents the average of all proxy data in each respective region. The z-scores represent the deviation from the mean of each proxy record, with blue representing "wetter" (more positive z-score) and red "drier" (more negative z-score) conditions. Each plot contains an inset map for each region with red squares for the location of published proxy records and location of the two new records generated in this study are shown as coloured stars (see Methods and Supplementary Data 1). Grey shaded bars represent the warm MIS 5 sub-stages (MIS 5e, c and a). associated with significant local sensitivities to climate forcing. Separation of the proxy records into their respective ISM and EASM-dominated realms ( Supplementary Fig. 3) highlights that there is a distinct heterogeneous response of each monsoonal subsystem to climate forcing across MIS 5. Not only is Northern Hemisphere Solar Insolation (NHSI) thought to be a key driver in the monsoonal circulation over the Asian landmass 29,30 but the effects of sea level 31 , ice volume and greenhouse gas concentrations 32 change in Atlantic Meridional Overturning Circulation strength 33 and bi-hemispheric forcing [24][25][26]34 are all implicated to exert an influence in dictating past monsoonal strength in the Asian Monsoon regions across glacial-interglacial timescales, although proxy records from the IPWP infer climate to be insensitive to high-latitude forcing in response to a dominance of local insolation forcing [35][36][37] or sea-level changes 31 . Thus, the observed deviations among the individual proxy records in Regions 3 and 6 at the millennial scale highlight the complex interplay of regional and local sensitivities in dictating hydroclimate within these realms requiring proxy-based paleoclimate interpretations to recognize these complexities.
In comparison, Regions 1, 2, 4 and 5 ( Fig. 2a-d) exhibit a more coherent relationship between the compiled proxy records within each region at the millennial scale with periods of extreme wet/ dry conditions akin to the warm/cold sub-stages of MIS 5 being identified. This is further reinforced by a good correlation at the local scale as emulated by PaleoPGEM ( Supplementary Fig. 2). Modelling studies have shown that rainfall within the Atlantic Ocean and its margins correlate well with the zonal mean ITCZ 9,10 . Proxy records from these regions also invoke an intimate link of hydroclimate variations within Central and South America 5,7 , northern and southern Africa [38][39][40] to be linked with meridional ITCZ shifts as a result of changes in the interhemispheric Atlantic SST gradient coupled with precession-paced changes in summer insolation 4 ( Supplementary Fig. 4). However, temporal discrepancies between precession, inter-hemispheric heat transport and emulated precipitation do exist across MIS 5 ( Supplementary Fig. 4), highlighting the additional importance of local processes in determining hydroclimate variability. The existence of intervals containing limited anti-phasing between respective regional northern and southern hemisphere records ( Supplementary Fig. 5) supports limited hemispheric antiphasing found in studies from the Holocene 10-12 and in recent times 41 , therefore, underscoring that hydroclimate variability cannot be explained by ITCZ shifts alone.
Meridional and zonal precipitation variability. Figure 3 shows the qualitative values of "wetter" and "drier" assigned to proxy data (blue and red data points, respectively) (Methods), to simplify and further disentangle past patterns in hydroclimate variability that is limited by the transient approach (Fig. 3). Results are shown for each MIS 5 sub-stage with respect to the average Holocene (Methods) and are superimposed on our emulation of spatial rainfall. The results reveal pronounced tropical regional heterogeneities in response to climate forcing across MIS 5 (Fig. 3). Overall, there is good proxy-model data convergence, particularly during both the warmest time-slice MIS 5e (Fig. 3a) and coolest time-slice MIS 5b (Fig. 3e). Discrepancies between the proxy and model anomalies with respect to the Holocene may be the result of the intermediate-resolution nature of the model 21 (Methods) or due to the various factors associated with the proxies (Methods).
In general, the Southern Hemisphere tropics (0-30°S) experience an intensified rainfall regime during the cooler substages of MIS 5d and b, especially over land masses (Fig. 3d, e), and a reduction in rainfall during the warmer sub-stages ( Fig. 3a-c) with respect to the Holocene. This intensification of the Southern Hemisphere tropical hydroclimate is consistent with a shift of the ITCZ southwards (northwards) during cool (warm) Northern (Southern) Hemisphere episodes 4 . Although this hemispheric antiphase pattern in hydroclimate variability across MIS 5 may be associated with ITCZ shifts, it is insufficient for explaining our reconstructed spatial pattern of precipitation, which is coherent from both our proxy collations and emulation (Fig. 3). The modelled precipitation anomalies and proxy data reveal intensified rainfall regimes both north and south of the equator in South America and Africa during MIS 5e (Fig. 3a), whereas distinct zonal wet-dry patterns exist throughout MIS 5 30   5c-Holocene, c MIS 5a-Holocene, d MIS 5d-Holocene, e MIS 5b-Holocene. Proxy data are superimposed; blue square = "wetter" than Holocene average, red square = "drier" than Holocene average and black square = no change (Methods and Supplementary Data 1). Emulated data are integrated from native model resolution (64 × 32 grid cells) onto 32 × 16 resolution, to smooth out grid-scale variability. Maps were produced using Panoply (source: https://www.giss.nasa.gov/tools/panoply/). (Fig. 3). Across MIS 5, the western Amazon retained a consistently "wet" climate (Fig. 3). In contrast, the eastern Amazon, under the influence of the South American Monsoon, experienced drier conditions during the warmer MIS 5 sub-stages (MIS 5e, c and a) (Fig. 3a-c) and wetter, stronger monsoonal conditions during MIS 5d and b (Fig. 3d, e).
Regional precipitation gradients similarly characterize the African continent across MIS 5. Intensified rainfall in northeasternequatorial Africa during MIS 5e contrasts with dry SE African conditions (Fig. 3a). This contrasting inter-hemispheric climatic regime during MIS 5e is thought to have provided conditions favourable for promoting human dispersal out of southern Africa northwards 42,43 . Evidence for expansive ancient lake systems 44,45 and cave deposits 46 in the presently arid regions of northern Africa, Arabia and the Levant (Fig. 3a) point to MIS 5e being associated with an intensified and more expanded rainfall regime in comparison to the Holocene. The penetration of monsoonal rains into northern Africa is typically impeded, associated with the high surface albedo of the Sahara Desert 47 . However, the wet conditions of MIS 5e, primarily influenced by increased NHSI during precession min 48 , would have facilitated a northward intrusion of monsoonal rain via land-surface feedbacks associated with the subsequent vegetation of the Sahara and expanded lake systems 49 (Fig. 3a). By MIS 5b (Fig. 3e), these wet conditions had largely dissipated in northeastern Africa and the Arabian Peninsula, highlighting the unique climatic window associated with the earlier half of MIS 5. The African continent was further characterized by a dipole-like precipitation pattern south of the equator, which was prevalent across all MIS 5 time slices, highlighting the competing influence of the Atlantic ITCZ to the west vs. the influence of Indian Ocean dynamics to the east 50,51 .
In contrast to the intensified hydroclimate regime in the Levant and Arabian Peninsula during MIS 5e, both subsystems of the Asian Monsoon, the ISM and EASM, experience reduced monsoonal precipitation compared to the Holocene (Fig. 3a). This may be indicative of a lack of homogenous meridional shift of the zonal mean ITCZ or influenced by biases within the proxy signal, which remains complex within this region (Methods). The overall reduced monsoonal strength associated with MIS 5e in comparison to the Holocene over Asia (Fig. 3a) highlights discrepancies in application of the "Global Monsoon" term due to the potentially greater regional influence of land-surface feedbacks in modulating and sustaining the response of rainfall over the northeastern African tropics 52 compared to that over the Asian domain. In addition, the differential rainfall response to warming during the Holocene and MIS 5e may highlight the competing influence of thermodynamic and dynamic responses, and the respective regional importance under warm climates controlling the loci of rainfall as is observed with the greater amount of rainfall falling over the equatorial Indian Ocean compared to the Holocene (Fig. 3a).
Concluding remarks. The integrated approach applied in this study, exploiting both multi-proxy data and model output, indicates that the spatial pattern of tropical hydroclimate in response to past climate changes cannot be explained solely by commonly ascribed mechanisms related to a meridional framework, specifically ITCZ shifts, within paleoclimate studies. Dipole-like precipitation patterns are pervasive throughout MIS 5, existing on all land masses straddling the tropics (Fig. 3), highlighting the competing influences of different ocean basins 53 in influencing the local internal climate dynamic response to forcing. This is further reinforced by the differing response of the regional monsoon systems, such as extreme "dry" conditions across the Indian Monsoon region, to climate warming and cooling across MIS 5 in comparison to the Holocene, indicating the importance of competing internal and external boundary conditions in preconditioning monsoonal strength and loci of monsoon rainfall. Consequentially, the "Global Monsoon" term is undermined. Importantly, these regional nuances in the response of tropical hydroclimate to climate forcing has had important implications for human migration "out-of-Africa" 42,43 and Amazonian biodiversity across MIS 5 54,55 . The underappreciated complexity in the response of tropical rainfall to past climate forcing warrants more concentrated efforts in compiling proxy records coupled with model outputs across large spatial scales. Encouragingly, the results of this study highlight that an approximation of past hydroclimate can be discerned only when such multi-location and multi-proxy approaches are taken. A need for more concentrated proxy-model integrated studies of past warm intervals is needed considering that, regionally, hydroclimate does not respond uniformly with respect to climate forcing and thereby has important implications for future predictions of tropical rainfall patterns in response to climate change. Details of the procedure for cleaning prior to trace-element analysis and subsequent derivation of temperature (SST) and seawater δ 18 O (δ 18 O sw-IVC ) have been described in ref. 24 ; the same approaches are applied here. Samples from Sites U1446 and U1448 were subjected to a full trace-element cleaning method 56 with slight modifications, namely an extended clay-removal step was performed to ensure removal of fine clays, which might bias Mg/Ca results 57 . Samples were analysed at the Open University using an Agilent Technologies Triple-Quad Inductively Coupled Plasma-Mass Spectrometer 8800. Samples from Site U1446 and Site U1448 were run at 10 and 25 p.p.m. Ca, respectively. Reproducibility of a synthetic monitor standard of Mg/Ca ratio of 3.33 mmol/mol was 0.06 mmol/mol during the analytical window. Following ref. 58 , elements were normalized to Ca and ratios were calculated using calibration derived from matrix-matched standards with varying trace-element concentrations. Contaminant ratios (Al/Ca and Fe/Ca) were monitored to assess any contamination by clays and organic material. Samples that returned values of Al/Ca > 300 µmol/mol were discarded; in total, 20 G. ruber ss samples with high Al/Ca values, including two samples with anomalously low Mg/Ca ratios from Site U1446, were rejected. Linear interpolation between data points was applied to provide a first-order estimation of missing Mg/Ca (SST) values. It is inferred that this is sufficient for deriving the local δ 18 O sw-IVC due to U1446 G. ruber ss δ 18 O c being overwhelmed by the effects of the local hydrology (precipitation and runoff) as exhibited by the significant magnitudes of change throughout MIS 5 ranging up to 2.5‰ (Supplementary Fig. 6). Furthermore, annual temperature changes within the modern-day Bay of Bengal are limited by a couple of degrees 59 and, therefore, the estimated Mg/Ca values should be within error. SST was estimated using the widely applied multi-species equation from the North Atlantic 60 with a correction for a 10% reduction in Mg/Ca following ref. 61 , associated with the reductive method; Mg/Ca = 0.342(±0.02)exp(0.09(±0.03) × T). To isolate the local δ 18 O sw , both a temperature and an ice volume correction was applied to the calcite δ 18 O; T C ¼ 14:9ð ± 0:1Þ À 4:8ð ± 0:08Þ ðδ 18 O c À ðδ 18 O sw À 0:27ÞÞ 62 . A sea-level correction was applied by exploiting the benthic δ 18 O of both U1446 and U1448, which were scaled to the Waelbroeck sea-level curve 63 .

Methods
Compilation of low-latitude hydroclimate proxy records. A selection of existing paleoclimate records that cover a latitudinal band of~30°N to 30°S, extending wholly or partially across MIS 5 (~130-74 ka) and have been inferred to capture past changes in hydroclimate, were compiled (Supplementary Data 1).
Paleoclimate records, which had a resolution of 2 ka or better, with publicly archived data, were gathered to perform the time-series comparisons (Supplementary Data 1). Application of this criteria resulted in a remaining database of 55 hydroclimate records, including the 2 new records generated in this study (Supplementary Fig. 1). All records were normalized and linear interpolation was performed to resample datasets onto a common age scale at a 1 ka resolution (i.e., the same resolution as the model time series). Records were subsequently assigned "Regions", to compare with the regional time series of the emulation of precipitation. Largely, these followed the latitude-longitude of the proxy site location, which denoted the corresponding region (Fig. 1). However, records within the IPWP 35,36 and northeast South American margin 64 were assigned to Regions 6 and 4 rather than Regions 3 and 2, respectively.
Further, we apply a qualitative approach 19,31,65 to compare regional precipitation variability with model derived anomaly maps (MIS 5 time-slice average-Holocene (12-0 ka average signal)), where values 1 (wetter than Holocene, blue), 0 (no change, black) and −1 (drier than Holocene, red) were assigned to paleoclimate records based on the original interpretation and/or assessment with the average Holocene (12-0 ka) values where available (Supplementary Data 1). This strategy was followed due to the complex nature of proxies applied to reconstruct past hydroclimate (Methods, see "Reconstructing past hydroclimate") and consequent limitations in producing and comparing any meaningful quantitative estimates (e.g., rainfall (mm)) unlike studies that synthesize past SST changes (°C) 16 . We believe the time-slice strategy applied (i.e., average signal across the interstadials (MIS 5a, c and e) and stadial (MIS 5b, d) periods compared to the average Holocene) permits a representative qualification of the proxy signal with chronological and proxy errors suppressed compared to studies that employ a single time-point approach (e.g., see ref. 19 ).
Age model. The age model and tuning strategy of Site U1446 has been initially described in Nilsson-Kerr et al. 24 and has subsequently undergone further refinement for this study. U1448 benthic foraminifera (Cibicidoides wuellerstorfi and Uvigerina spp.) oxygen isotope (δ 18 O) (S. C. Clemens, unpublished) was similarly transferred onto the Antarctic Ice Core Chronology (AICC2012) 66,67 ( Supplementary Fig. 7).
To reduce potential offsets induced by differing tuning strategies, a common age model was applied, where possible, to all marine records that were extracted for time-series analysis. U1446 benthic δ 18 O was used as a reference target for benthic δ 18 O records to transfer onto AICC2012 chronology. Benthic δ 18 O of sediment cores GeoB7920-2 68 , GeoB9528-3 69 (Supplementary Fig. 8). Linear interpolation was subsequently applied between tie-points to construct the new-age models.
Precipitation emulation. Emulated precipitation fields were derived using PaleoPGEM 21 . The approach applies Gaussian process emulation of the singular value decomposition of ensembles of runs from the intermediate-complexity atmosphere-ocean GCM PLASIM-GENIE 22 with varied boundary-condition forcing (CO 2 , orbit and ice volume). The underlying simulations were performed at spectral T21 and atmospheric resolution (5.625°) with 10 vertical layers, and a matching ocean grid with 16 logarithmically spaced depth levels, using the optimized parameter set 21 . Spatial fields of annual average precipitation were then emulated at 1000-year intervals, driven by the time series of scalar boundarycondition forcing, and assuming the climate is in quasi-equilibrium. Twenty Gaussian process emulators were built and applied, and the mean is presented in Figs. 2 and 3. In Fig. 3, emulated data are integrated from the native model resolution (64 × 32 grid cells; Supplementary Fig. 11) onto 32 × 16 resolution, to smooth out the grid-scale structure. The spatial distribution and magnitude of PaleoPGEM climate and corresponding uncertainties have been shown to be comparable to model inter-comparisons (Section 9 of ref. 21 ) in state-of-the-art climate models as measured by spatial fields of mean and variance across ensembles of the Paleoclimate Modelling Intercomparison Project (PMIP) simulations of the Last Glacial Maximum and Mid Holocene 21 , noting that PaleoPGEM uncertainties derive from the emulation step of the reconstructions, and thus underestimate the full simulation uncertainty. Supplementary Fig. 12b compares the PaleoPGEM precipitation anomaly between 127 ka and Preindustrial with the results from the lig12k simulation for CMIP6-PMIP4 ( Supplementary Fig. 12a). Considering the intermediate resolution of PaleoPGEM, there are discernible similarities in the large-scale precipitation anomalies reconstructed with the multimodel ensemble average (Supplementary Fig 12a). Furthermore, the PLASIM-GENIE simulation of present day precipitation shares similar large-scale spatial patterns with that of measured precipitation ( Fig. 1 and Supplementary Fig. 13). Our compilation of the available proxy data in time-series format for low-latitude hydroclimate (30°N-30°S) (Supplementary Fig. 1d) and for both northern hemisphere low latitudes (0-30°N) (Supplementary Fig. 1f) and southern hemisphere low latitudes (0-30°S) ( Supplementary Fig. 1h), as well as the time-varying PaleoPGEM emulations, correspond well with the CCSM3 model output for terrestrial rainfall ( Supplementary Fig. 1c, e, g) and to the record of the ΔDE* ( Supplementary Fig. 1a,i), inferred as a measure of the large-scale low-latitude hydroclimate signal 23 .
Reconstructing past hydroclimate. Sea-level changes altering the transport of riverine material to the ocean can bias applications of grain size, clay mineral content and elemental concentrations of marine sediment cores, inferred to track riverine runoff and continental weathering, as a result of changes in continental precipitation intensity 91,92 . Furthermore, normalization of X-Ray Fluorescencederived element concentrations to calcium (e.g., Fe/Ca and Ti/Ca), widespread within the paleoclimate studies (and hence, several examples can be found in Supplementary Data 1), may be biased by the effect of biogenic dilution 93 , whereas proxies inferred to trace past vegetation changes (e.g., pollen assemblage and plant leaf waxes) in response to changes in rainfall regime are similarly influenced by changes in coastal environments 94,95 . Inferences of past vegetation as representing hydroclimate may further be obscured by an increasing amount of anthropogenic activities occurring during the Holocene (e.g., anthropogenic fire and burning). Isotopic tracers (δ 18 O, δD) of past rainfall hosted in carbonates (e.g., speleothems and microfossils (foraminifera)) or plant leaf waxes are commonly applied to capture changes in past rainfall intensity (e.g., see refs. 69,96,97 ). However, the isotopic composition of precipitation received at a site is fundamentally determined by several processes including (but not limited to) evaporation (evapotranspiration), moisture transport and the associated amalgamation of differing moisture sources 3,98-101 . The isotopic signal is further determined by a mixture of local processes inherent to the proxy archive. The isotopic composition of leaf wax can be complicated by different plant types 102,103 , whereas speleothem δ 18 O can be influenced by cave temperature, soil and karstic processes 104,105 . Residuals of foraminifera δ 18 O (seawater δ 18 O) (having corrected for the effects of temperature and sea level, which carry their own respective uncertainties [106][107][108][109] ) commonly used as a means of capturing changes in surface freshening in marginal locations involve large uncertainties if a conversion to salinity is applied 110 and any freshening by rainfall alone is expected to lack any discernible signal. Modern empirical studies 104,111,112 , application of isotope-enabled models 99,101,113 and refined calibration approaches 114,115 have been employed, to disentangle the varying complications to permit confidence in paleoclimate reconstructions, but with limited success.
In addition to the difficulties associated with the acquisition of an authentic record of past rainfall variability, each proxy record is further hampered by the time-averaging of sediments and associated chronological uncertainties. The timeslice approach taken in this study associated with investigation of past hydroclimate patterns across MIS 5 compared to the Holocene (Fig. 3) alleviates some chronological biases due to averaging of several thousands of years (e.g., MIS 5e, 130-115 ka). Thereby, ensuring a representative signal from that time slice has been acquired unlike single time-point approaches. However, the time-series approach to regional hydroclimate reconstruction (Fig. 1) will suffer discrepancies associated with chronological errors, although these have been minimized by ensuring marine records are all placed on a consistent chronology. Chronological errors exert a large influence on paleoclimate proxy records and, thus, our ability to constrain timing/phasing will always be limited to an extent. This is further impacted by the proxy amalgamating a system-wide response to climate forcing (e.g., river runoff and sea-level changes).
Considering the aforementioned factors, it is promising that the average combined proxy signal (i.e., the combined average of all the proxy record z-scores from each respective region) approximates past hydroclimate as that estimated by the emulation for the different regions (Fig. 2). This underlines the importance of multi-proxy, basin-wide strategies, to circumvent locally inflicted biases, both climatically and carrier induced, when inferring past hydroclimate.
Holocene vs. MIS 5e. MIS 5e has largely been regarded as representing one of the warmer interglacials of the Late Pleistocene experiencing similar CO 2 levels to PI 116 ( Supplementary Fig. 14d) with sea level 6-9 m higher than present 117 . Primarily, it was in the high latitudes of the Northern Hemisphere where anomalous warming with reconstructed temperatures was witnessed, suggesting >5°C warmer than PI 118 , whereas the low latitudes remained similar to the present-day SSTs 16 . The climatic conditions of MIS 5e are mainly a result of preconditioning by a combination of high eccentricity ( Supplementary Fig. 14b,i) enhancing the effects of precession min (Supplementary Fig. 14b,ii) 119 . This orbital configuration resulted in enhanced seasonality (stronger insolation received during Northern Hemisphere summers) compared to the Holocene but the overall mean annual insolation was lower 120 . Obliquity maximum was reached prior to MIS 5e at~131 ka (Supplementary Fig. 14 a). Meanwhile, precession peaked at~127 ka prior to its decline bỹ 117 ka (Supplementary Fig. 14b,ii). Across MIS 5e, the orbital configuration varied (Supplementary Fig. 14); therefore, latitudinal temperature gradients and seasonality were not stationary throughout this interval, and subsequently would have been accompanied by variations in the tropical hydroclimate. The averaged time-slice approach taken, i.e., MIS 5e 130-116 ka and Holocene 12-0 ka, encompasses both the increase into and attainment of peak interglacial conditions, as well as the descent from peak interglacial conditions and, thus, ideally represent the average tropical hydroclimate response under a generally warm climatic state.

Data availability
Proxy data generated in this study (U1446 and U1448) and regional time-series precipitation emulation data can be found archived at https://doi.org/10.1594/ PANGAEA.920668. All proxy data used in the analysis are publicly available and detailed information can be found in Supplementary Data 1.