Pliocene-Pleistocene biostratigraphic datums of ODP Site 184-1148

Based on the stable isotopic analysis of planktonic and benthic foraminifers from Ocean Drilling Program Core 1148 of the northern South China Sea (SCS), Pliocene-Pleistocene isotope stratigraphy and events have been reconstructed. The benthic foraminiferal delta18O record shows that the Pacific intermediate water had a greater influence upon the SCS or the Pacific deep water above ~2600 m was warmer before ~3.2Ma than at present. After that, the benthic delta18O conspicuously increased during the ~3.2-2.5 Ma period, in correspondence to the formation of the Northern Hemisphere ice sheet, whereas the planktonic delta18O signal suggests a stepwise overall decrease of sea surface temperature during the ~2.2-0.9 Ma period. Compared to the equatorial Pacific records, the decrease in planktonic (Globigerinoides ruber) delta13C during the ~3.2-2.2 Ma period is particularly striking, suggesting that fertility of surface water increased noticeably. According to the modern delta13C distribution of G. ruber in the northern SCS, it is inferred that the East Asian winter monsoon strengthened during this interval. Afterwards, there were several conspicuous decreases of G. ruber delta13C at ~1.7, 1.3, 0.9, 0.45 and 0.15 Ma BP, that is, about every 0.4 Ma, suggesting that the East Asian winter monsoon became episodically stronger. This is confirmed by changes in relative abundance of planktonic foraminifer species Neogloboquadrina dutertrei, a typical East Asian winter monsoon proxy. The deepwater delta13C of the SCS is close to that of the Pacific, but lighter than that of the Atlantic, implying that the pattern of deep water originating mainly from the Atlantic and through the Pacific entering the SCS existed at least since the early Pliocene. After 1.4 Ma, the benthic delta13C signal decreased conspicuously but with a periodicity of ~100 ka, suggesting that the deep-water ventilation of the SCS was reduced, probably corresponding to a decrease of the North Atlantic Deep Water and/or further isolation of the SCS deep basin from the Pacific during glaciations.

Age model of ODP Site 184-1146 (Table 1)

Clay mineral assemblages at ODP Site 1146 in the northern South China Sea are used to investigate sediment source and transport processes and to evaluate the evolution of the East Asian monsoon over the past 2 Myr. Clay minerals consist mainly of illite (22-43%) and smectite (12-48%), with associated chlorite (10-30%), kaolinite (2-18%), and random mixed-layer clays (5-22%). Hydrodynamic and mineralogical studies indicate that illite and chlorite sources include Taiwan and the Yangtze River, that smectite and mixed-layer clays originate predominantly from Luzon and Indonesia, and that kaolinite is primarily derived from the Pearl River. Mineral assemblages indicate strong glacial-interglacial cyclicity, with high illite, chlorite, and kaolinite content during glacials and high smectite and mixed-layer clay content during interglacials. During interglacials, summer enhanced monsoon (southwesterly) currents transport more smectite and mixed-layer clays to Site 1146 whereas during glacials, enhanced winter monsoon (northerly) currents transport more illite and chlorite from Taiwan and the Yangtze River. The ratio (smectite+mixed layers)/(illite+chlorite) was adopted as a proxy for East Asian monsoon variability. Higher ratios indicate strengthened summer-monsoon winds and weakened winter-monsoon winds during interglacials. In contrast, lower ratios indicate a strongly intensified winter monsoon and weakened summer monsoon during glacials. Spectral analysis indicates the mineral ratio was dominantly forced by monsoon variability prior to the development of large-scale glaciation at 1.2 Myr and by both monsoon variability and the effects of changing sea level in the interval 1.2 Myr to present.

(Table 1) Element concentrations in detrital sediments of ODP Site 184-1148

The chemical index of alteration (CIA) and elemental ratios that are sensitive to chemical weathering, such as Ca/Ti, Na/Ti, Al/Ti, Al/Na, Al/K, and La/Sm, were analyzed for detrital sediments at Ocean Drilling Program Site 1148 from the northern South China Sea to reveal information of chemical weathering in the source regions during the early Miocene. High CIA values of ~80, coupled with high Al/Ti and Al/Na and low Na/Ti and Ca/Ti, are observed for the sediments at ~23 Ma, indicating a high chemical weathering intensity in the north source region, i.e., south China. This was followed by gradual decreases in Al/Ti, Al/Na, La/Sm, and Al/K ratios, as well as the CIA values, and increases in Ca/Ti and Na/Ti ratios. These records together with other paleoclimate proxies, such as black carbon d13C and benthic foraminifer d18O, give reliable information on the climate changes in south China. Our results show that the climate in south China was warm and humid in the early Miocene (~23 Ma) according to the chemical weathering records. The humidity in south China decreased from the early Miocene to Present with several fluctuations centering at approximately 15.7 Ma, 8.4 Ma, and 2.5 Ma, coincident with the global cooling since the middle Miocene. These climate changes implied that the summer east Asian monsoon has dramatically affected south China in the early Miocene, whereas the influence of the summer monsoon on this region has decreased continuously since that time, probably because of the intensification of the winter monsoon. Such an evolution for the east Asian monsoon is different from that for the Indian monsoon.

(Table 1) Age determination of sediment core SO139-74KL

A multiproxy record has been acquired from a piston core (SO139-74KL) taken offshore southern Sumatra, an area which is situated in the southwestern sector of the tropical Indo-Pacific Warm Pool. The high-resolution data sets (X-ray fluorescence, total organic carbon, and C37 alkenones) were used to track changes in paleoproductivity, freshwater budget, and sea surface temperature (SST) of the tropical climate system at orbital time scales over the past 300 ka. Our paleoclimatic data show that enhanced marine paleoproductivity was directly related to strengthening of coastal upwelling during periods of increased boreal summer insolation and associated SE monsoon strength with a precessional cyclicity. Changes in freshwater supply were primarily forced by precession-controlled changes in boreal NW winter monsoon rainfall enclosing an additional sea level component. SST variations of 2°-5°C occurred at eccentricity and precessional cyclicity. We suggest that the sea surface temperature variability off southern Sumatra is predominantly related to three major causes: (1) variations in upwelling intensity; (2) an elevated freshwater input into the southern Makassar Strait leading to reduced supply of warmer surface waters from the western Pacific and increased subsurface water transport via the Indonesian Throughflow into the Indian Ocean; and (3) long-term changes in the intensity or frequency of low-latitude climate phenomena, such as El Niño-Southern Oscillation.

Water temperature reconstructions based on paired Uk'37 and Tex86 records, and paired surface and subsurface Mg/Ca records

This dataset contains the collection of available published paired Uk'37 and Tex86 records spanning multi-millennial to multi-million year time scales, as well as a collection of Mg/Ca-derived temperatures measured in parallel on surface and subsurface dwelling foraminifera, both used in the analyses of Ho and Laepple, Nature Geoscience 2016. As the signal-to-noise ratios of proxy-derived Holocene temperatures are relatively low, we selected records that contain at least the last deglaciation (oldest sample >18kyr BP).

(Appendix B) Benthic foraminifera in Indian Ocean surface sediments

Multivariate analysis was performed on percentages of 46 species of unstained deep-sea benthic foraminifera from 131 core-top to near-core-top samples (322-5013 m) from across the Indian Ocean. Faunal data are combined with GEOSECS geochemical data to investigate any relationship between benthic foraminifera (assemblages and species) and deep-sea properties. In general, benthic foraminifera show a good correlation to surface productivity, organic carbon flux to the sea floor, deep-sea oxygenation and, to a lesser extent, to bottom temperature, without correlation with the water depths. The foraminiferal census data combined with geochemical data has enabled the division of the Indian Ocean into two faunal provinces. Province A occupies the northwestern Indian Ocean (Arabian Sea region) where surface primary production has a major maximum during the summer monsoon season and a secondary maximum during winter monsoon season that leads to high organic flux to the seafloor, making the deep-sea one of the most oxygen-deficient regions in the world ocean, with a pronounced oxygen minimum zone (OMZ). This province is dominated by benthic foraminifera characteristic of low oxygen and high organic food flux including Uvigerina peregrina, Robulus nicobarensis, Bolivinita pseudopunctata, Bolivinita sp., Bulimina aculeata, Bulimina alazanensis, Ehrenbergina carinata and Cassidulina carinata. Province B covers southern, southeastern and eastern parts of the Indian Ocean and is dominated by Nuttallides umbonifera, Epistominella exigua, Globocassidulina subglobosa, Uvigerina proboscidea, Cibicides wuellerstorfi, Cassidulina laevigata, Pullenia bulloides, Pullenia osloensis, Pyrgo murrhina, Oridorsalis umbonatus, Gyroidinoides (= Gyroidina) soldanii and Gyroidinoides cf. gemma suggesting well-oxygenated, cold deep water with low (oligotrophic) and pulsed food supply.

Planktonic foraminiferal faunas from surface-sediment samples from the South China Sea and the western Pacific

We present 30 new planktonic foraminiferal census data of surface sediment samples from the South China Sea, recovered between 630 and 2883 m water depth. These new data, together with the 131 earlier published data sets from the western Pacific, are used for calibrating the SIMMAX-28 transfer function to estimate past sea-surface temperatures. This regional SIMMAX method offers a slightly better understanding of the marginal sea conditions of the South China Sea than the linear transfer function FP-12E, which is based only on open-ocean data. However, both methods are biased toward the tropical temperature regime because of the very limited data from temperate to subpolar regions. The SIMMAX formula was applied to sediment core 17940 from the northeastern South China Sea, with sedimentation rates of 20-80 cm/ka. Results revealed nearly unchanged summer temperatures around 28°C for the last 30 ky, while winter temperatures varied between 19.5°C in the last glacial maximum and 26°C during the Holocene. During Termination 1A, the winter estimates show a Younger Dryas cooling by 3°C subsequent to a temperature optimum of 24°C during the Bölling=Alleröd. Estimates of winter temperature differences between 0 and 100 m water depth document the seasonal variations in the thickness of the mixed layer and provide a new proxy for estimating past changes in the strength of the winter monsoon.

Magnetic susceptibility and particle size for the Shilou red clay section on the eastern Chinese Loess Plateau

The Chinese Loess Plateau red clay sequences display a continuous alternation of sedimentary cycles that represent recurrent climatic fluctuations from 2.58 Ma to the Miocene. Deciphering such a record can provide us with vital information on global and Asian climatic variations. Lack of fossils and failure of absolute dating methods made magnetostratigraphy a leading method to build age models for the red clay sequences. Here we test the magnetostratigraphic age model against cyclostratigraphy. For this purpose we investigate the climate cyclicity recorded in magnetic susceptibility and sedimentary grain size in a red clay section previously dated 11Myr old with magnetostratigraphy alone. Magnetostratigraphy dating based on only visual correlation could potentially lead to erroneous age model. In this study the correlation is executed through the iteration procedure until it is supported by cyclostratigraphy; i.e., Milankovitch cycles are resolved in the best possible manner. Our new age model provides an age of 5.2Ma for the Shilou profile. Based on the new age model, wavelet analysis reveals the well-preserved 400 kyr and possible 100 kyr eccentricity cycles on the eastern Chinese Loess Plateau. Further, paleomonsoon evolution during 2.58-5.2Ma is reconstructed and divided into three intervals (2.58-3.6Ma, 3.6-4.5Ma, and 4.5-5.2Ma). The upper part, the youngest stage, is characterized by a relatively intensified summer monsoon, the middle stage reflects an intensification of the winter monsoon and aridification in Asia, and the earliest stage indicates that summer and winter monsoon cycles may have rapidly altered. The use of cyclostratigraphy along withmagnetostratigraphy gives us an effectivemethod of dating red clay sequences, and our results imply that many presently published age models for the red clay deposits should be perhaps re-evaluated.

(Table DR1) Biofacies, percentages of important benthic foraminifera and isotope data at ODP Site 758

The Indian monsoon system, as recorded by ocean-floor biota (benthic foraminifera) at Ocean Drilling Program Site 758 in the eastern equatorial Indian Ocean, has varied dramatically over the past 5.5 m.y., long after the onset of the monsoons at 10-8 Ma. Benthic foraminifera that thrive with high productivity year-round were common before the formation of Northern Hemisphere continental ice sheets ca. 3.1-2.5 Ma, indicating that the summer (southwest) monsoon had high intensity and long seasonal duration. Ca. 2.8 Ma benthic faunas became dominated by taxa that flourish with a seasonally strongly fluctuating food supply, indicating that the northeast (winter) monsoon, during which primary productivity is relatively low, increased in duration and strength to form a system similar to that of today. The change occurred coeval with the initiation of the Northern Hemisphere glaciation, documenting a close link between the development of the Indian monsoon and Northern Hemisphere glaciation.

Sea surface temperature and primary productivity reconstruction of sediment core SO130-275KL from the northeastern Arabian Sea during the late Holocene

Variability in the oceanic environment of the Arabian Sea region is strongly influenced by the seasonal monsoon cycle of alternating wind directions. Prominent and well studied is the summer monsoon, but much less is known about late Holocene changes in winter monsoon strength with winds from the northeast that drive convective mixing and high surface ocean productivity in the northeastern Arabian Sea. To establish a high-resolution record of winter monsoon variability for the late Holocene, we analyzed alkenone-derived sea surface temperature (SST) variations and proxies of primary productivity (organic carbon and d15N) in a well-laminated sediment core from the Pakistan continental margin. Weak winter monsoon intensities off Pakistan are indicated from 400 B.C. to 250 A.D. by reduced productivity and relatively high SST. At about 250 A.D., the intensity of the winter monsoon increased off Pakistan as indicated by a trend to lower SST. We infer that monsoon conditions were relatively unstable from ~500 to 1300 A.D., because primary production and SST were highly variable. Declining SST and elevated biological production from 1400 to 1900 A.D. suggest invigorated convective winter mixing by strengthening winter monsoon circulation, most likely a regional expression of colder climate conditions during the Little Ice Age on the Northern Hemisphere. The comparison of winter monsoon intensity with records of summer monsoon intensity suggests that an inverse relationship between summer and winter monsoon strength exists in the Asian monsoon system during the late Holocene, effected by shifts in the Intertropical Convergence Zone.

Age determination of sediment core GIK16523-1

High-resolution studies of a planktonic foraminifer core record from the South China Sea (SCS) (31KL: 18°45.4'N, 115°52.4'E, water depth 3360 m) reveal changes driven by ice-volume forcings in the climate of the East Asian monsoon in the western Pacific marginal sea during the late Quaternary. The analyses of planktonic foraminifer faunal abundance data from the core indicate significant variations in the relative abundances of the dominant taxa over the past 100,000 years since the isotope stage 5. The transfer function estimates of faunal sea surface temperatures (SST) correlate well with a long-term (104-105 years) trend of global glaciation. About 65,000 years ago, there was an observable change in the mode of SST variability as many low-latitude records have shown. These findings suggest that the SCS surface circulation and the East Asian winter monsoon systems are mainly driven by variations in global glaciation levels. The association of surface ocean cooling in the SCS with global climatic events suggests that fluctuations in the strength of the East Asian winter monsoon may be linked to shifts in the latitudinal position of the westerly winds and the Siberian high-pressure system.

Age model and palaeoclimate records over the past 22,000 years of sediment core GeoB10053-7

The Australian-Indonesian monsoon is an important component of the climate system in the tropical Indo-Pacific region. However, its past variability, relation with northern and southern high-latitude climate and connection to the other Asian monsoon systems are poorly understood. Here we present high-resolution records of monsoon-controlled austral winter upwelling during the past 22,000 years, based on planktic foraminiferal oxygen isotopes and faunal composition in a sedimentary archive collected offshore southern Java. We show that glacial-interglacial variations in the Australian-Indonesian winter monsoon were in phase with the Indian summer monsoon system, consistent with their modern linkage through cross-equatorial surface winds. Likewise, millennial-scale variability of upwelling shares similar sign and timing with upwelling variability in the Arabian Sea. On the basis of element composition and grain-size distribution as precipitation-sensitive proxies in the same archive, we infer that (austral) summer monsoon rainfall was highest during the Bølling-Allerød period and the past 2,500 years. Our results indicate drier conditions during Heinrich Stadial 1 due to a southward shift of summer rainfall and a relatively weak Hadley cell south of the Equator. We suggest that the Australian-Indonesian summer and winter monsoon variability were closely linked to summer insolation and abrupt climate changes in the northern hemisphere.

Benthic foraminifera assemblages of the Gulf of Aden

Benthic and selected planktic foraminifera and stable isotope records were determined in a piston core from the Gulf of Aden, NW Arabian Sea that spans the last 530 ka. The benthic foraminifera were grouped into four principal assemblages using Q-mode Principal Component Analyses. Comparison of each of these assemblages with the fauna of the nearby regions enabled us to identify their specific environmental requirements as a function of variability in food supply and strength of the oxygen minimum zone and by that to use them as indicators of surface water productivity. The benthic foraminiferal productivity indicators coupled with the record of Globigerina bulloides, a planktic foraminifer known to be sensitive to productivity changes in the region, all indicate higher productivity during glacial intervals and productivity similar to present or even reduced during interglacial stages. This trend is opposite to the productivity pattern related to the SW summer monsoon of the Arabian Sea and indicates the role of the NE winter monsoon on the productivity of the Gulf of Aden. A period of exceptionally enhanced productivity is recognized in the Gulf of Aden region between ~60 and 13 kyr indicating the intensification of the NE winter monsoon to its maximal activity. Contemporaneous indication of increased productivity in other parts of the Arabian Sea, unexplained so far by the SW summer monsoon variability, might be related to the intensification of the NE winter monsoon. Another prominent event of high productivity, second in its extent to the last glacial productivity event is recognized between 430 and 460 kyr. These two events seem to correspond to periods of similar orbital positioning of rather low precession (and eccentricity) amplitude for a relatively long period. Glacial boundary conditions seem to control to a large extent the NE winter monsoon variability as also indicated by the dominance of the 100 ka cycle in the investigated time series. Secondary in their importance are the 23 and 41 ka cycles which seem also to contribute to the NE monsoonal variability. Following the identification of productivity events related to the NE winter monsoon in the Gulf of Aden, it is possible now to extend this observation to other parts of the Arabian Sea and consider the contribution of this monsoonal system to the productivity fluctuations preserved in the sedimentary records.

Climate variability of the last 1000 years in the NW Pacific: high resolution, multi-biomarker records from Lake Toyoni

The East Asian Monsoon (EAM) is an active component of the global climate system and has a profound social and economic impact in East Asia and its surrounding countries. Its impact on regional hydrological processes may influence society through industrial water supplies, food productivity and energy use. In order to predict future rates of climate change, reliable and accurate reconstructions of regional temperature and rainfall are required from all over the world to test climate models and better predict future climate variability. Hokkaido is a region which has limited palaeo-climate data and is sensitive to climate change. Instrumental data show that the climate in Hokkaido is influenced by the East Asian Monsoon (EAM), however, instrumental data is limited to the past ~150 years. Therefore down-core climate reconstructions, prior to instrumental records, are required to provide a better understanding of the long-term behaviour of the climate drivers (e.g. the EAM, Westerlies, and teleconnections) in this region. The present study develops multi-proxy reconstructions to determine past climatic and hydrologic variability in Japan over the past 1000 years and aid in understanding the effects of the EAM and the Westerlies independently and interactively. A 250-cm long sediment core from Lake Toyoni, Hokkaido was retrieved to investigate terrestrial and aquatic input, lake temperature and hydrological changes over the past 1000-years within Lake Toyoni and its catchment using X-Ray Fluorescence (XRF) data, alkenone palaeothermometry, the molecular and hydrogen isotopic composition of higher plant waxes (δD(HPW)). Here, we conducted the first survey for alkenone biomarkers in eight lakes in the Hokkaido, Japan. We detected the occurrence of alkenones within the sediments of Lake Toyoni. We present the first lacustrine alkenone record from Japan, including genetic analysis of the alkenone producer. C37 alkenone concentrations in surface sediments are 18µg C37 g−1 of dry sediment and the dominant alkenone is C37:4. 18S rDNA analysis revealed the presence of a single alkenone producer in Lake Toyoni and thus a single calibration is used for reconstructing lake temperature based on alkenone unsaturation patterns. Temperature reconstructions over the past 1000 years suggest that lake water temperatures varies between 8 and 19°C which is in line with water temperature changes observed in the modern Lake Toyoni. The alkenone-based temperature reconstruction provides evidence for the variability of the EAM over the past 1000 years. The δD(HPW) suggest that the large fluctuations (∼40‰) represent changes in temperature and source precipitation in this region, which is ultimately controlled by the EAM system and therefore a proxy for the EAM system. In order to complement the biomarker reconstructions, the XRF data strengthen the lake temperature and hydrological reconstructions by providing information on past productivity, which is controlled by the East Asian Summer monsoon (EASM) and wind input into Lake Toyoni, which is controlled by the East Asian Winter Monsoon (EAWM) and the Westerlies. By combining the data generated from XRF, alkenone palaeothermometry and the δD(HPW) reconstructions, we provide valuable information on the EAM and the Westerlies, including; the timing of intensification and weakening, the teleconnections influencing them and the relationship between them. During the Medieval Warm Period (MWP), we find that the EASM dominated and the EAWM was suppressed, whereas, during the Little Ice Age (LIA), the influence of the EAWM dominated with time periods of increased EASM and Westerlies intensification. The El Niño Southern Oscillation (ENSO) significantly influenced the EAM; a strong EASM occurred during El Niño conditions and a strong EAWM occurred during La Niña. The North Atlantic Oscillation, on the other hand, was a key driver of the Westerlies intensification; strengthening of the Westerlies during a positive NAO phase and weakening of the Westerlies during a negative NAO phase. A key finding from this study is that our data support an anti-phase relationship between the EASM and the EAWM (e.g. the intensification of the EASM and weakening of the EAWM and vice versa) and that the EAWM and the Westerlies vary independently from each other, rather than coincide as previously suggested in other studies.