|Home | About | Journals | Submit | Contact Us | Français|
Millennial-scale cooling events termed Heinrich Stadials punctuated Northern Hemisphere climate during the last glacial period. Latitudinal shifts of the intertropical convergence zone (ITCZ) are thought to have rapidly propagated these abrupt climatic signals southward, influencing the evolution of Southern Hemisphere climates and contributing to major reorganisation of the global ocean-atmosphere system. Here, we use neodymium isotopes from a marine sediment core to reconstruct the hydroclimatic evolution of subtropical Australia between 90 to 20 thousand years ago. We find a strong correlation between our sediment provenance proxy data and records for western Pacific tropical precipitations and Australian palaeolakes, which indicates that Northern Hemisphere cooling phases were accompanied by pronounced excursions of the ITCZ and associated rainfall as far south as about 32°S. Comparatively, however, each of these humid periods lasted substantially longer than the mean duration of Heinrich Stadials, overlapping with subsequent warming phases of the southern high-latitudes recorded in Antarctic ice cores. In addition to ITCZ-driven hydroclimate forcing, we infer that changes in Southern Ocean climate also played an important role in regulating late glacial atmospheric patterns of the Southern Hemisphere subtropical regions.
During the last glacial, episodes of massive iceberg discharges in the Northern Hemisphere were accompanied by substantial cooling and weakening of the Atlantic Meridional Oceanic Circulation (AMOC)1, leading ultimately to out-of-phase warming periods in Antarctica2,3. These abrupt climate change events, known as Heinrich Stadials (HS), were also rapidly transferred southward via latitudinal migrations of the ITCZ, as directly inferred from proxy records indicating reduced precipitation in northern tropics and subtropical regions4,5, and increased monsoon rainfall in the Southern Hemisphere6,7. This interhemispheric atmospheric seesaw is thought to have caused global displacement of air masses, with presumably major impact on higher-latitude climates in the Southern Hemisphere8,9 and global ocean circulation patterns2,10,11. To date however, it is unclear whether the evolution of Southern Hemisphere climates during the last glacial period was mainly driven by northern (via ITCZ latitudinal variations) and/or southern high-latitude climate forcing.
To provide new insights on this issue, we have analysed a well-dated12 marine sediment core (MD03-2607; 36°57.64′S, 137°24.39′E; 865m water depth) recovered from near the mouth of the River Murray, the end point of the Murray-Darling Basin (MDB), Australia’s largest river system (Fig. 1). The MDB lies in the subtropical climate zone between about 25°S and 37°S, displaying a marked latitudinal gradient of contrasting geological and climatic settings. In the northern part of the basin (~25–32°S), an area with weak dominance of ITCZ-driven summer monsoon rainfall, the Darling River sub-basin drains large outcrops of Mesozoic and Cenozoic terrains. In contrast, the southern River Murray sub-basin (~32–37°S) is dominated by Palaeozoic rocks, and is more strongly influenced by winter precipitation associated with Southern Hemisphere westerly winds (SHWW). As a consequence, the use of geochemical proxies allowing discrimination between Darling versus Murray sources in sedimentary records of ancient MDB discharge can provide continuous and integrated information on past hydroclimate variability in subtropical Australia, and its potential link to ITCZ and SHWW strength and/or latitudinal shifts.
We reconstructed the composition of past sediment discharge using neodymium isotopic ratios (143Nd/144Nd, or εNd in epsilon notation), focusing on the fine-grained (<2μm) clay-rich detrital fraction of the sediment. In addition to being exported from river basins with presumably minimum transfer time, clay-size fractions are also less prone than coarser sedimentary particles to mineralogical sorting and hence better suited for provenance studies, especially near upwind arid source areas. The three main potential contributors of fine-grained sediments to the studied site, i.e. the Darling River, the River Murray and southern Australian dust-source regions, are characterized by distinctive Nd isotopic signatures13,14 (εNd=−2.4±2.4, 1SD; −9.5±0.9 and −12.0±2.9, respectively; Fig. 2; Table S1), which makes Nd isotopes particularly well-suited for provenance studies in this area (Fig. 1). Neodymium isotopes remain largely unaffected by erosion and transport processes15, so that the average εNd composition of each potential provenance area is expected to have remained relatively unchanged over the last glacial period.
The εNd detrital record of core MD03-2607 indicates highly variable contributions from MDB sources and southern Australian dusts, starting, from about 86 to 68 thousand years before present (kyr BP), with a pronounced shift towards more radiogenic Nd isotopic compositions (from about εNd=−10 to −6; Fig. 2). This trend reflects enhanced contribution from Darling river-borne material and reduced inputs from southern Australia dusts and/or River Murray particles. During that period, sea-level dropped by about 50meters (Fig. 2), resulting in a progressive ‘migration’ of the proximity of the studied area to the Murray mouth, and producing a more efficient sediment transfer from the MDB to site MD03-2607 (ref. 16). An inverse trend from εNd ~−5 to −13 also characterizes the last deglaciation period at the nearby site MD03-2611 (Fig. 2), coinciding with an almost 90-m steep sea-level rise between ~17 and 9kyr BP and also interpreted in terms of overwhelming contribution of aeolian dust particles relative to riverine detrital inputs from the MDB (ref. 16). Overall, the observed εNd trend between about 90 and 20ka correlates well with palaeoprecipitation records for the tropical western Pacific region5,17,18,19 (Fig. 3a,b,c). Rainfall patterns in the tropics are strongly controlled by the position and intensity of the ITCZ, which varies on seasonal timescales from north to south of the Equator (Fig. 1). During the last ice age, large southward ITCZ migrations have accompanied both the onset of MIS-4 and the abrupt millennial-scale HS events20. In the northwest Pacific region, these ITCZ shifts are illustrated by abrupt drops in precipitation recorded both by Borneo speleothems5 (Fig. 1; Fig. 3a) and, farther south (~3°S), at the northern Papua New Guinea (PNG) margin (site MD05-2920; Fig. 1), as inferred from substantial decreases of sedimentary Ti/K ratios during HS-3-6-7b events17,18 (Fig. 3b). At site MD05-2920, Ti/K ratios have been interpreted as a proxy for the presence of coarse-grained river-borne terrigenous material, presumably increasing during past wet periods when intensifying rainfall caused more intense erosion and export of coarser-grained riverine particles to the ocean, and vice versa18. For the same time intervals, foraminifera trace element proxy data suggest episodes of enhanced river discharge in southern PNG (~9°S; Fig. 1). The use of Nd/Ca ratios in planktonic foraminifera at this latter southern PNG site (MD05-2925; Fig. 3c), taken as an indicator for the intensity of REE-rich dissolved riverine inputs from Papua New Guinean rivers, shows that the sudden drops in precipitation recorded at Borneo and northern PNG were accompanied by more intense rainfall in the southerly tropical regions19. This pattern agrees well with results obtained from dust flux reconstructions in the central equatorial Pacific21, which provide evidence for a major latitudinal shift of the ITCZ by at least 4 degrees in this area during an earlier HS event. Within the uncertainties of dating (see Supplementary Information), our εNd provenance profile for palaeo-MDB sediments matches remarkably well with the southern PNG river-discharge record. In core MD03-2607, the periods of HS events and glacial maxima (Marine Isotope Stage 4 and MIS-2) correspond to more radiogenic εNd signatures (Fig. 3d), indicating enhanced sediment input from the Darling River watershed. While a reduction in SHWW in the River Murray sub-basin could possibly account for the observed trend, the good agreement between MDB and southern PNG proxy records seems more consistent with southward ITCZ migrations leading to intensifying rainfall in the northern Darling River sub-basin relative to the southern MDB area.
Our proxy data are also consistent with palaeoclimatic records for Australian lakes22,23,24 (Fig. 3e). The palaeoshoreline reconstruction of Lake Eyre, the largest central Australia inland basin draining summer-monsoon-fed watersheds (Fig. 1), indicate a trend towards generally dryer conditions since the last interglacial (MIS-5) period, punctuated with alternating periods of high and low (or dry) lake levels22. Lake Frome, a separate depocentre of the same basin, also displays some high lake shorelines, though not always synchronous with Lake Eyre levels23. Despite the paucity of dates and the poor constraints on the timing of upper and lower boundaries of wet/dry lake phases (inherent to low-precision luminescence-derived chronologies), the presumed periods of high lake stands are consistent with inferred episodes of southward ITCZ migrations, as far as about 32°S south (i.e. the boundary between the two Darling and Murray sub-basins), and associated shifts towards more radiogenic εNd signatures recorded at site MD03-2607.
In contrast to the abrupt seesaw-like tropical precipitation records of Borneo during HS events, the εNd profile for MDB sediment provenance suggests instead a more gradual evolution of rainfall in subtropical Australia. In particular, the εNd shifts towards more radiogenic values associated with every major southward excursions of the ITCZ appear to have started earlier, and persisted later, than the corresponding HS events as recorded in Borneo (Fig. 4). There are clear inherent uncertainties when establishing the chronostratigraphy of marine sediment cores beyond the limit of radiocarbon (~50kyr BP), which prevent any detailed comparison with more precise U-Th-dated speleothem records. In addition, possible temporal biases in our detrital εNd record may exist due to temporary sediment storage within the drainage basin and/or post-depositional bioturbation effect. Nevertheless, bearing in mind these potential limitations, the relatively good agreement observed between our εNd profile and PNG river-discharge/Australian lake records still suggests that wet episodes in late glacial Australia lasted substantially longer (between ~3 to >5ka each, as visually estimated from both MD03-2607 and palaeolake data22,23,24) than the average duration of the southward ITCZ migrations associated with HS events (<2ka; as inferred from Borneo5 and Chinese6 cave monsoon records; Fig. 4). All the above suggests therefore that rainfall patterns in subtropical Australia during the last glacial period were modulated by additional mechanisms, rather than being simply driven by ITCZ-driven Northern Hemisphere forcing.
In fact, our εNd proxy record for hydroclimate variability in Australia also exhibits strong similarities with the late glacial evolution of Antarctic temperatures25 (Fig. 4d), hence suggesting a possible link to southern high-latitude forcing. Except for MIS-4, the most radiogenic εNd values observed between 35 and 65kyr BP (i.e. indicative of higher sediment contributions from the Darling River and tributaries relative to the River Murray sub-basin) appear to coincide with the warm phases recorded in Antarctic ice-cores that followed the onset of HS events (Fig. 4d). Antarctic climate variability over the last ice age is thought to have been closely connected to Northern Hemisphere forcing too, with gradual warming periods being generally attributed to enhanced heat storage in the Southern Hemisphere, in response to AMOC collapse during HS events and other abrupt cold stadials3. During these periods, Antarctic warming may have also been promoted through atmospheric teleconnections, via ITCZ-driven southward shifts of the SHWW belt2,10. In turn, the evolution of Antarctic climate during the last ice age most probably played a major role in controlling the position of the Subtropical Front in the Southern Ocean26, with presumably direct impact on the position/intensity of SHWW and associated rainfall. This complex interhemispheric interplay is illustrated by a δ18O stalagmite record from New Zealand’s South Island (~42°S)9, which reveals that wet periods in the Southern Hemisphere mid-latitudes prevailed at times of southward ITCZ migrations, but also when cooler conditions in Antarctica led to northward shifts of the STF and strengthening of the SHWW (Fig. 4c). Previous studies have already argued for a connection between high-latitude climate changes in the Southern Hemisphere and the hydrologic evolution of monsoon-dominated regions during the Late Quaternary period7,27,28,29. Using high-resolution elemental profiles in a series of sediment cores from the Timor Sea, Kuhn et al. (ref. 27) proposed that the onset of Antarctic warming that followed the HS1 climate event led to a southward shift of the ITCZ over northwestern Australia. A detailed comparison of the climate variability recorded in the Hulu Cave and polar ice-cores also led Rohling et al. (ref. 28) to suggest that millennial-scale monsoon variability was dominated by Southern Hemisphere climate change during glacial times, when monsoon was weak overall. Taken together, the above-mentioned studies provide support that the late glacial hydroclimate variability of subtropical Australia, as inferred from our εNd profile for core MD03-2607, was probably driven by combined influences of both northern- (via ITCZ) and southern (via SHWW) high-latitude forcings. For instance, it is possible that the poleward shifts of the Subtropical Front during warm periods in Antarctica, which led to temporary southward atmospheric displacement of the SHWW belt, also acted as a factor delaying the contemporaneous northward retreat of the ITCZ. Alternatively, enhanced greenhouse forcing associated with the episodic rises of atmospheric CO2 that accompanied the Antarctic warming phases could have resulted in enhanced heat low-pressure cells over Australia, thereby accentuating the southward ‘pull’ of the ITCZ across the SH mid-latitudes27. While the exact mechanism remains elusive, these interactions between northern and southern high-latitude climates most likely drove important fluctuations of interhemispheric temperature/pressure gradients28, which may have resulted in protracted rainfall episodes in the Southern Hemisphere subtropical regions during the last glacial period.
The detailed description of chemical preparation and analytical methods can be found elsewhere15. Briefly, about 3g of bulk sediment were treated successively with 10% (v/v) acetic acid (AA), mixed 15% (v/v) AA and 0.05M hydroxylamine hydrochloride (HH), and 5% hydrogen peroxide (H2O2) solutions, in order to remove any carbonate, Fe-Mn oxyhydroxide and organic components, respectively. The sequential leaching procedure used in this study was developed to minimize partial dissolution of the silicate component30. Recent work has shown however that easily alterable silicate minerals such as fresh volcanic components could undergo partial dissolution with the use of mixed AA-HH solutions31. In the case of our study, one cannot exclude that a minor proportion of the volcanogenic component hosted by Darling River-borne sediments dissolved during our second leaching step. Considering a simple mixing model between various proportions of Darling- versus Murray-borne sediments, and assuming arbitrarily that 10% of the Darling sediment endmember dissolved during the AA+HH leaching phase, the obtained Nd isotope composition of the resulting mixed sediment would be biased towards Murray endmember by less than 0.5 εNd units. Importantly, this bias, while being proportional to the proportion of Darling River sediment, would not affect the observed downcore trend for Nd isotopes, and hence the conclusions of our study.
Clay-size fractions (about <2μm) were then separated from the residual detritus by centrifugation. For Nd isotopic analyses, about 100mg of clay-size sediment powder were digested by alkaline fusion. Neodymium isotopic measurements were performed at the Pôle Spectrométrie Océan (Brest, France) using a Thermo Scientific Neptune multi-collector ICPMS, after Nd purification by conventional ion chromatography. Mass bias corrections on Nd were made with the exponential law, using 146Nd/144Nd=0.7219, and 143Nd/144Nd corrected values were normalized to a JNdi-1 value of 143Nd/144Nd=0.512115. Repeated analyses of a JNdi-1 standard solution during the two measurement sessions of this study gave 143Nd/144Nd of 0.512106±0.000005 (2SD, n=22) and 0.512121±0.000009 (2SD, n=33), hence corresponding to an external reproducibility of ~±0.10ε and±0.17ε (2SD), respectively (Table S1). Epsilon Nd values (εNd) were calculated using 143Nd/144Nd=0.512630 (ref. 32).
How to cite this article: Bayon, G. et al. Extensive wet episodes in Late Glacial Australia resulting from high-latitude forcings. Sci. Rep. 7, 44054; doi: 10.1038/srep44054 (2017).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
We thank the crews of R/V Marion Dufresne, especially Yvon Balut, and all participants of the AUSCAN cruise (MD131; PI: P.D.D.) for their assistance at sea. This work was funded via an IEF Marie Curie fellowship to G.B. (Grant No. FP7-PEOPLE-2012-IEF 327778).
The authors declare no competing financial interests.
Author Contributions G.B. and P.D.D. conceived the project and wrote the article. P.D.D. led the AUSCAN cruise, to which J.W.M. participated. P.D.D. collected the marine sediment and dust samples. M.D.N. analysed dust samples. G.B. and Y.G. analysed marine sediment samples. All authors contributed to discussions, interpretation of the results and manuscript writing.