Orbital-scale Variability Research Articles

North Pacific warmth synchronous with the Miocene Climatic Optimum

Peak Neogene warmth and minimal polar ice volumes occurred during the Miocene Climatic Optimum (MCO, ca. 16.95−13.95 Ma) followed by cooling and ice sheet expansion during the Middle Miocene Climate Transition (MMCT, ca. 13.95−12.8 Ma). Previous records of northern high-latitude sea surface temperatures (SSTs) during these global climatic transitions are limited to Atlantic sites, and none resolve orbital-scale variability. Here, we present an orbital-resolution alkenone SST proxy record from the subpolar North Pacific that establishes a local maximum of SSTs during the MCO as much as 16 °C warmer than modern with rapid warming initiating the MCO, cooling synchronous with Antarctic ice sheet expansion during the MMCT, and high variability on orbital time scales. Persistently cooler North Pacific SST anomalies than in the Atlantic at equivalent latitudes throughout the Miocene suggest enhanced Atlantic northward heat transport under a globally warm climate. We conclude that a global forcing mechanism, likely elevated greenhouse gas concentrations, is the most parsimonious explanation for synchronous global high-latitude warmth during the Miocene.

CO2-driven and orbitally driven oxygen isotope variability in the Early Eocene

Abstract. Paleoclimate reconstructions of the Early Eocene provide important data constraints on the climate and hydrologic cycle under extreme warm conditions. Available terrestrial water isotope records have been primarily interpreted to signal an enhanced hydrologic cycle in the Early Eocene associated with large-scale warming induced by high atmospheric CO2. However, orbital-scale variations in these isotope records have been difficult to quantify and largely overlooked, even though orbitally driven changes in solar irradiance can impact temperature and the hydrologic cycle. In this study, we fill this gap using water isotope–climate simulations to investigate the orbital sensitivity of Earth's hydrologic cycle under different CO2 background states. We analyze the relative difference between climatic changes resulting from CO2 and orbital changes and find that the seasonal climate responses to orbital changes are larger than CO2-driven changes in several regions. Using terrestrial δ18O and δ2H records from the Paleocene–Eocene Thermal Maximum (PETM), we compare our modeled isotopic seasonal range to fossil evidence and find approximate agreement between empirical and simulated isotopic compositions. The limitations surrounding the equilibrated snapshot simulations of this transient event and empirical data include timing and time interval discrepancies between model and data, the preservation state of the proxy, analytical uncertainty, the relationship between δ18O or δ2H and environmental context, and vegetation uncertainties within the simulations. In spite of the limitations, this study illustrates the utility of fully coupled, isotope-enabled climate models when comparing climatic changes and interpreting proxy records in times of extreme warmth.

Resolving Cenozoic climate pattern debate in East Asia: Insights from orbital-scale oscillations

Paleoclimate reconstruction and aeolian deposits in China indicate a major environmental reorganization – the formation of a nonzonal monsoon-dominated climate pattern – in East Asia in the late Oligocene and the early Miocene. In contrast, a zonal planetary-wind-dominated climate pattern has prevailed during the Paleocene and Eocene. However, this classical paradigm describing the climate transition has been questioned by recent studies. Here, our study emphasizes the importance of considering orbital-scale climate oscillations to better understand and reassess the paradigm. By examining the halite and mudstone sequence in the Jianghan Basin and conducting numerical experiments using the Norwegian Earth System model, we propose the concept of a “green Yangtze” phenomenon, which utilizes orbital-scale climate oscillations to reconcile the wet and dry environmental indicators in Eastern China during the Paleogene. It is essential to differentiate the Paleogene climate, characterized by intense orbital-scale wet and arid oscillations, from the Neogene monsoon-dominated climate. The classic paradigm remains valuable in understanding the long-term mean environment, albeit without capturing the orbital-scale details. Future research should focus on uncovering orbital-scale variations to enhance our understanding of climate evolution and dynamics in East Asia during the Cenozoic.

Millennial and orbital-scale variability in a 54 000-year record of total air content from the South Pole ice core

Abstract. The total air content (TAC) of polar ice cores has long been considered a potential proxy for past ice sheet elevation. Recent work, however, has shown that a variety of other factors also influence this parameter. In this paper we present a high-resolution TAC record from the South Pole ice core (SPC14) covering the last 54 000 years and discuss the implications of the data for interpreting TAC from ice cores. The SPC14 TAC record shows multiple features of interest, including (1) long-term orbital-scale variability, (2) millennial-scale variability in the Holocene and last glacial period, and (3) a period of stability from 35 to 25 ka. The longer, orbital-scale variations in TAC are highly correlated with integrated summer insolation (ISI), corroborating the potential of TAC to provide an independent dating tool via orbital tuning. Large millennial-scale variability in TAC during the last glacial period is positively correlated with past accumulation rate reconstructions as well as δ15N-N2, a firn thickness proxy. These TAC variations are too large to be controlled by direct effects of temperature and too rapid to be tied to elevation changes. We propose that grain size metamorphism near the firn surface explains these changes. We note, however, that at sites with different climate histories than the South Pole, TAC variations may be dominated by other processes. Our observations of millennial-scale variations in TAC show a different relationship with accumulation rate than observed at sites in Greenland.

A transient coupled general circulation model (CGCM) simulation of the past 3 million years

Abstract. Driven primarily by variations in the earth's axis wobble, tilt, and orbit eccentricity, our planet experienced massive glacial/interglacial reorganizations of climate and atmospheric CO2 concentrations during the Pleistocene (2.58 million years ago (Ma)–11.7 thousand years ago (ka)). Even after decades of research, the underlying climate response mechanisms to these astronomical forcings have not been fully understood. To further quantify the sensitivity of the earth system to orbital-scale forcings, we conducted an unprecedented quasi-continuous coupled general climate model simulation with the Community Earth System Model version 1.2 (CESM1.2, ∼3.75∘ horizontal resolution), which covers the climatic history of the past 3 million years (3 Myr). In addition to the astronomical insolation changes, CESM1.2 is forced by estimates of CO2 and ice-sheet topography which were obtained from a simulation previously conducted with the CLIMBER-2 earth system model of intermediate complexity. Our 3 Ma simulation consists of 42 transient interglacial/glacial simulation chunks, which were partly run in parallel to save computing time. The chunks were subsequently merged, accounting for spin-up and overlap effects to yield a quasi-continuous trajectory. The computer model data were compared against a plethora of paleo-proxy data and large-scale climate reconstructions. For the period from the Mid-Pleistocene Transition (MPT, ∼1 Ma) to the late Pleistocene we find good agreement between simulated and reconstructed temperatures in terms of phase and amplitude (−5.7 ∘C temperature difference between Last Glacial Maximum and Holocene). For the earlier part (3–1 Ma), differences in orbital-scale variability occur between model simulation and the reconstructions, indicating potential biases in the applied CO2 forcing. Our model-proxy data comparison also extends to the westerlies, which show unexpectedly large variance on precessional timescales, and hydroclimate variables in major monsoon regions. Eccentricity-modulated precessional variability is also responsible for the simulated changes in the amplitude and flavors of the El Niño–Southern Oscillation. We further identify two major modes of planetary energy transport, which played a crucial role in Pleistocene climate variability: the first obliquity and CO2-driven mode is linked to changes in the Equator-to-pole temperature gradient; the second mode regulates the interhemispheric heat imbalance in unison with the eccentricity-modulated precession cycle. During the MPT, a pronounced qualitative shift occurs in the second mode of planetary energy transport: the post-MPT eccentricity-paced variability synchronizes with the CO2 forced signal. This synchronized feature is coherent with changes in global atmospheric and ocean circulations, which might contribute to an intensification of glacial cycle feedbacks and amplitudes. Comparison of this paleo-simulation with greenhouse warming simulations reveals that for an RCP8.5 greenhouse gas emission scenario, the projected global mean surface temperature changes over the next 7 decades would be comparable to the late Pleistocene glacial-interglacial range; but the anthropogenic warming rate will exceed any previous ones by a factor of ∼100.

Palynological implications for the paleoclimate and paleoceanographic reconstruction of the East Sea since the early Pleistocene at IODP site U1430

Palynomorphs, including pollen and dinoflagellate cysts from marine sediments, can provide a history of paleoclimate and paleoceanography. We present new palynological profiles based on a core from the Integrated Ocean Drilling Program Site U1430, which was drilled in the Ulleung Plateau East Sea. This core covers the past 2 Ma and spans the Mid-Pleistocene Transition (MPT, 1.2–0.8 Ma) and Mid-Brunhes Event (MBE, ∼0.43 Ma). Terrestrial and marine palynomorph records indicate orbital-scale variations in paleoclimate and paleoceanography, reflecting changes in the interconnections between vegetation, climate, and the ocean. Before the MPT, the high pollen taxa ratio and abundant temperate pollen suggest a relatively warm temperature associated with astronomical forcing-induced insolation changes. At ca. 1.3 Ma, boreal coniferous pollen increased, suggesting intensification of the East Asian winter monsoon (EAWM). The long-term trend of boreal coniferous pollen concentration was similar to that of benthic δ18O. This orbital scale boreal conifer pollen assemblage variation indicates Quaternary EAWM evolution controlled by changes in global ice volume and climatic cooling. The marine palynological data reflect the ocean circulation history of the study area, where sea level fluctuations and tectonic activity played important roles around the marginal sea. The abundance of warm-water dinoflagellate species after ca. 1.9 Ma, including Tuberculodinium vancampoae, indicates Tsushima Warm Current intrusion related to the expansion of the Korea (Tsushima) Strait. Nematosphaeropsis labyrinthus and other cold-water species increased markedly during the MPT, suggesting a low sea surface temperature and weakened vertical mixing of the East Sea due to the cooled climate and sea ice during the MPT. The marked sea level fluctuation associated with the glacial–interglacial cycling after the MBE resulted in an assemblage composition shift from cold to warm water species. Tsuga and Betula thrived during mild winters, primarily influenced by winter insolation. These genera were particularly prosperous in the early Pleistocene, coinciding with a strong intrusion of the Tsushima Warm Current (TWC). Our findings should be considered in developing a comprehensive understanding of the paleoclimate and paleoceanographic history of northern East Asia across the major climate transitions of the Pleistocene.

Milankovitch-paced South Asian monsoons during Marine Isotope Stage 5

The orbital-scale variations and associated driving mechanisms of the South Asian summer monsoon (SASM) and South Asian winter monsoon (SAWM) remain poorly understood. In this study, we present the first record of a 4.3-m loess-paleosol sequence from the southern Tibetan Plateau, which provides insight into the changes of these two monsoons over most of the last interglacial-glacial cycle. The age framework of this sequence was established using twenty-one samples with absolute luminescence techniques. The SASM was represented by frequency-dependent susceptibility in percentage and redness, while the SAWM was reflected by the recently proposed sorting coefficient of grain size. The results suggest an anti-phase relationship between the SASM and SAWM. During marine isotope stage (MIS) 5, variations of the SASM and SAWM showed a strong correlation with the Northern Hemisphere summer (July) insolation (NHSI) at 65°N, in terms of both phase and amplitude. Stronger SASM (SAWM) corresponded to higher (lower) July NHSI at 65°N, supporting the Milankovitch theory. The variation amplitude of the SASM record in this study significantly differs from the that of stalagmite δ18O from the nearby Tianmen Cave, which presented no distinct difference between MIS 5e, MIS 5c, and MIS 5a. During MIS 4 and MIS 3, the larger global ice volume and lower atmospheric CO2 could have enhanced the SAWM but weakened the SASM, which might have restrained the response amplitude of the South Asian monsoons to the July NHSI at 65°N.

Decoupled orbital-scale variability of late Pleistocene-Holocene monsoonal circulation and rainfall in East Asia

Decoupled orbital-scale variability of late Pleistocene-Holocene monsoonal circulation and rainfall in East Asia

Nature and origin of variations in pelagic carbonate production in the tropical ocean since the mid-Miocene (ODP Site 927)

Abstract. Marine plankton is an important component of the global carbon cycle. Whereas the production and seafloor export of organic carbon produced by the plankton, the biological pump, has received much attention, the long-term variability in plankton calcification, controlling the carbonate counter pump, remains less well understood. However, it has been shown that changes in pelagic calcification (biological compensation) could affect the ocean's buffering capacity and thus regulate global carbon budget on geological timescales. Here we use Neogene pelagic sediments deposited on the Ceara Rise in the tropical Atlantic to characterize the variability in pelagic carbonate production with a focus on warm climates. A re-evaluation of published records of carbonate accumulation at the Ceara Rise reveals a systematic increase in sedimentation rates since the late Miocene, but the carbonate accumulation rate does not show a clear trend. Instead, we observe substantial orbital timescale variability in carbonate accumulation, combined with a trend towards less carbonate on average at sites located below 4 km, likely due to the effect of carbonate dissolution. To evaluate long-term changes against possible orbital-scale variability, we generated new high-resolution records of carbonate accumulation rate at Ocean Drilling Program (ODP) Site 927 across two Quaternary interglacials (MIS 5 and MIS 9), the Pliocene warm period (MIS KM5) and the Miocene Climatic Optimum (MCO). We observe that the highest carbonate accumulation rates occurred during the Pliocene but that each of the studied intervals was characterized by large-magnitude orbital variability. Prominent variations in carbonate accumulation prior to the Quaternary preservation cycles appear to follow Earth obliquity and eccentricity. These results imply that pelagic carbonate accumulation in the tropical ocean, buffered from large temperature changes, varied on orbital timescales. The magnitude of the orbital-scale variability was similar or even higher than the long-term mean differences among the studied intervals. Since preservation can be excluded as a driver of these changes prior to the Quaternary, the observed variations must reflect changes in the export flux of pelagic biogenic carbonate. We conclude that the overall carbonate production by pelagic calcifiers responded to local changes in light, temperature, and nutrients delivered by upwelling, which followed long orbital cycles, as well as to long-term shifts in climate and/or ocean chemistry. The inferred changes on both timescales were sufficiently large such that when extrapolated on a global scale, they could have played a role in the regulation of the carbon cycle and global climate evolution during the transition from the Miocene warm climates into the Quaternary icehouse.

Orbital forcing of tropical climate dynamics in the Early Cambrian

According to modern atmospheric circulation models, the Intertropical Convergence Zone (ITCZ), as the Earth's meteorological equator, plays an essential role in the low-latitude hydrologic cycles. The limited availability of high-resolution tropical climate archives, especially from the Early Paleozoic Era, severely limits our understanding of ITCZ migration dynamics in deep time. Here we present high-resolution climate-proxy records (i.e., magnetic susceptibility (MS) and Zirconium/Aluminum (Zr/Al)) from tropical marine sediments of the ~526-million-year-old Qiongzhusi Formation in South China to investigate the link between orbitally forced insolation changes, ITCZ migration dynamics, and low-latitude climate processes. These orbital-scale variations in MS and Zr/Al series are interpreted as alternations between wet and dry cycles, controlled by monsoon intensity under the orbitally forced ITCZ-related paleo-Hadley Cell dynamics. Our results show that combined precession and obliquity orbital cycles had an impact on the Early Cambrian ITCZ migration. Specifically, the precession and obliquity forcing shift the mean position of the ITCZ latitudinally by changing the interhemispheric pressure contrasts, thus affecting the low latitude hydroclimate cycle. We report semi-precession cycles of 8.3–7.9 kyr, which were probably associated with the twice-annual passage of the ITCZ across the intertropical zone, consistent with the paleogeographical location of South China near the equator during the Early Cambrian. Observed ~1.1 – ~1.5 Myr eccentricity amplitude modulation (AM) cycles and ~ 1.0 – ~1.2 Myr obliquity AM cycles may provide geological evidence for the chaotic motion between Earth and Mars in the Early Cambrian. • Combined precession and obliquity cycles regulated the Early Cambrian ITCZ migration. • We report the semi-precession cycles of 7.9–8.3 kyr in the Early Cambrian. • Our study provides geological evidence for the chaotic Solar System in the Early Cambrian.

Antiphase response of the Indonesian–Australian monsoon to millennial-scale events of the last glacial period

Antiphase behaviour of monsoon systems in alternate hemispheres is well established at yearly and orbital scales in response to alternating sensible heating of continental landmasses. At intermediate timescales without a sensible heating mechanism both in-phase and antiphase behaviours of northern and southern hemisphere monsoon systems are recorded at different places and timescales. At present, there is no continuous, high resolution, precisely dated record of millennial-scale variability of the Indonesian–Australian monsoon during the last glacial period with which to test theories of paleomonsoon behaviour. Here, we present an extension of the Liang Luar, Flores, speleothem δ18O record of past changes in southern hemisphere summer monsoon intensity back to 55.7 kyr BP. Negative δ18O excursions (stronger monsoon) occur during Heinrich events whereas positive excursions (weaker monsoon) occur during Dansgaard-Oeschger interstadials—a first order antiphase relationship with northern hemisphere summer monsoon records. An association of negative δ18O excursions with speleothem growth phases in Liang Luar suggests that these stronger monsoons are related to higher rainfall amounts. However, the response to millennial-scale variability is inconsistent, including a particularly weak response to Heinrich event 3. We suggest that additional drivers such as underlying orbital-scale variability and drip hydrology influence the δ18O response.

Indian monsoon variability in the Mahanadi Basin over the last two glacial cycles and its implications on the Indonesian throughflow

The orbital-scale variability of the Indian summer monsoon (ISM) has been influenced by multiple factors, such as atmospheric CO2 concentration, global ice volume, and insolation. Proxies for weathering activity and paleo-productivity provide potential insights into the driving forces of its variability. We documented multi-proxy data at IODP Site U1445, located in the Mahanadi Basin of the northwestern Bay of Bengal, to find out ISM variability over the last 200 ka. The proxy records, such as Nd/Sr isotopes of detrital particles, clay mineral compositions of the fine-grained sediments, biogenic opal and CaCO3, organic carbon contents, and carbon isotopes of organic matter, represent sediment sources, weathering patterns, and paleo-productivity related to the ISM variability. Detrital Nd/Sr isotope data and clay mineral compositions suggest that the sediments at Site U1445 originated mainly from the Ganges, Brahmaputra, and Meghna rivers without dramatic provenance change between the glacial and interglacial periods. The weathering activity inferred from clay mineral compositions and the paleo-productivity shift reconstructed by biogenic opal and CaCO3 contents suggest that the land-sea interactions were closely linked to the ISM precipitation between the glacial and interglacial periods. High precipitation by the strong ISM resulted in intense chemical weathering and dominant biogenic opal deposition during the interglacial periods. In contrast, low precipitation by the weak ISM led to reduced chemical weathering and predominant CaCO3 deposition during the glacial periods. Further, the ISM variability driving the land-sea interactions in the Mahanadi Basin was modulated by the Indonesian throughflow (ITF). Our study emphasizes the role of low-latitude forcing of climatic changes in the strong relationship between the ISM and ITF over orbital periods, providing a base for future investigations.

Model-based orbital-scale precipitation δ18O variations and distinct mechanisms in Asian monsoon and arid regions

The past Asian precipitation δ18O (δ18Op) records from stalagmites and other deposits have shown significant orbital-scale variations, but their climatic implications and regional differences are still not fully understood. This study, as the first attempt of a 300-kyr transient stable isotope-enabled simulation, investigated the characteristics and mechanisms of the orbital-scale δ18Op variations in three representative regions of Asia: arid Central Asia (CA), monsoonal South Asia (SA) and monsoonal East Asia (EA). The modelling results showed that the variations in the CA, SA and EA annual δ18Op exhibited significant but asynchronous 23-kyr precession cycles. Further analyses revealed that although the precession-induced insolation variation was the ultimate cause of the δ18Op variation in all three regions, the dominant mechanisms and the involved physical processes were distinct among them. For the CA region, the rainy-season (November-March) temperature effect and water vapour transport by the westerly circulation were identified as the key precession-scale processes linking the October-February boreal mid-latitude insolation to the rainy-season or annual δ18Op. In the SA region, the rainy-season (June-September) precipitation amount effect and upstream depletion of the monsoonal water vapour δ18O served as the main mechanisms linking the rainy-season or annual δ18Op to the April-July insolation variation at the precession scale. For the EA region, however, the precession-scale annual δ18Op was mainly controlled by the late-monsoon (August-September) and pre-monsoon (April-May) water vapour transport patterns, which were driven by the July-August insolation and the global ice volume, respectively. These results suggest that the climatic implications of the orbital-scale Asia δ18Op variations are sensitive to their geographic locations as determined by the combined effects of insolation and regional circulation patterns associated with the respective rainy seasons. This study provides new insights into understanding the regional differences and formation mechanisms of the Asian orbital-scale δ18Op variations.

Orbital-scale hydroclimate variations in the southern Tibetan Plateau over the past 414,000 years

Speleothem oxygen isotope (δ18O) constitutes one of key proxies to reconstruct the Indian and East Asian summer monsoon (ISM and EASM) hydroclimate history on a wide range of timescales. Here, we use three new well-dated speleothem δ18O records from Tianmen cave, in the south-central Tibetan Plateau (TP), together with previous records from the same cave, to reconstruct a composite Tianmen record spanning over the past 414,000 years. The Tianmen δ18O record follows Northern Hemisphere summer insolation (NHSI), broadly coherent with speleothem δ18O records from both ISM and EASM regimes. In comparison with previous model simulations and modern hydroclimate analyses, the Tianmen δ18O was interpreted primarily as a proxy indicating the ISM intensity. The Tianmen δ18O record comprises both high and low NHSI time periods, which reveals an exceptional large amplitude (∼14‰) of δ18O variations at the precession band, possibly resulted from a combination effect of precipitation amount, moisture source, transport pathway, high elevation and rainout. Besides, a broad consistency on orbital-scale variability among Tianmen δ18O, δ13C and trace element records is found, suggesting that the ISM intensity regulates the effective rainfall at the cave site as well. Notably, Tianmen speleothem growth occurred merely during warm-wet interglacial periods, even during low NHSI times near the precession maximum, and in contrast, hiatuses persistently occurred throughout cold-dry glacial periods. This pattern is similar to that found in travertine formations in the region, both suggesting that the interglacial temperature might be a prerequisite for speleothem/travertine formations in the region, while the strong ISM precipitation in turn provides a sufficient condition. Our new Tianmen record, combined with a large set of speleothem records from both Asian Westerlies and Asian monsoon domains, gives further evidence to reinforce a previous hypothesis: a consistent temporal pattern in the δ18O of precipitation on orbital and millennial scales over the vast Asian monsoon and Westerlies regimes. As such, in order to reconcile the precipitation δ18O variations between speleothem and TP ice cores, we argue that the Guliya ice-core chronology may have a basal age of ∼70 ka BP, which is supported by recent dating results from the Guliya ice cap.

Secular and orbital-scale variability of equatorial Indian Ocean summer monsoon winds during the late Miocene

Abstract. In the modern northern Indian Ocean, biological productivity is intimately linked to near-surface oceanographic dynamics forced by the South Asian, or Indian, monsoon. In the late Pleistocene, this strong seasonal signal is transferred to the sedimentary record in the form of strong variance in the precession band (19–23 kyr), because precession dominates low-latitude insolation variations and drives seasonal contrast in oceanographic conditions. In addition, internal climate system feedbacks (e.g. ice-sheet albedo, carbon cycle, topography) play a key role in monsoon variability. Little is known about orbital-scale monsoon variability in the pre-Pleistocene, when atmospheric CO2 levels and global temperatures were higher. In addition, many questions remain open regarding the timing of the initiation and intensification of the South Asian monsoon during the Miocene, an interval of significant global climate change that culminated in bipolar glaciation. Here, we present new high-resolution (<1 kyr) records of export productivity and sediment accumulation from International Ocean Discovery Program Site U1443 in the southernmost part of the Bay of Bengal spanning the late Miocene (9 to 5 million years ago). Underpinned by a new orbitally tuned benthic isotope stratigraphy, we use X-ray fluorescence-derived biogenic barium variations to discern productivity trends and rhythms. Results show strong eccentricity-modulated precession-band productivity variations throughout the late Miocene, interpreted to reflect insolation forcing of summer monsoon wind strength in the equatorial Indian Ocean. On long timescales, our data support the interpretation that South Asian monsoon winds were already established by 9 Ma in the equatorial sector of the Indian Ocean, with no apparent intensification over the latest Miocene.

Fractionation of O<sub>2</sub>∕N<sub>2</sub> and Ar∕N<sub>2</sub> in the Antarctic ice sheet during bubble formation and bubble–clathrate hydrate transition from precise gas measurements of the Dome Fuji ice core

Abstract. The variations of δO2/N2 and δAr/N2 in the Dome Fuji ice core were measured from 112 m (bubbly ice) to 2001 m (clathrate hydrate ice). Our method, combined with the low storage temperature of the samples (−50 ∘C), successfully excludes post-coring gas-loss fractionation signals from our data. From the bubbly ice to the middle of the bubble–clathrate transition zone (BCTZ) (112–800 m) and below the BCTZ (>1200 m), the δO2/N2 and δAr/N2 data exhibit orbital-scale variations similar to local summer insolation. The data in the lower BCTZ (800–1200 m) have large scatter, which may be caused by millimeter-scale inhomogeneity of air composition combined with finite sample lengths. The insolation signal originally recorded at the bubble close-off remains through the BCTZ, and the insolation signal may be reconstructed by analyzing long ice samples (more than 50 cm for the Dome Fuji core). In the clathrate hydrate zone, the scatter around the orbital-scale variability decreases with depth, indicating diffusive smoothing of δO2/N2 and δAr/N2. A simple gas diffusion model was used to reproduce the smoothing and thus constrain their permeation coefficients. The relationship between δAr/N2 and δO2/N2 is markedly different for the datasets representing bubble close-off (slope ∼ 0.5), bubble–clathrate hydrate transformation (∼1), and post-coring gas loss (∼0.2), suggesting that the contributions of the mass-independent and mass-dependent fractionation processes are different for those cases. The method and data presented here may be useful for improving the orbital dating of deep ice cores over the multiple glacial cycles and further studying non-insolation-driven signals (e.g., atmospheric composition) of these gases.

Benthic foraminiferal response to orbital-scale variability in primary productivity off the Portuguese margin over the last 1.3 Myr

Detailed analyses of benthic foraminiferal fauna and planktonic foraminiferal stable carbon isotopes at IODP Site U1391 were performed to reconstruct orbital-scale variability in primary productivity off the Portuguese margin over the last 1.3 Myr. Q-mode factor analysis of benthic foraminiferal census data from the >250 μm fraction reveals four assemblages, which is mainly controlled by the organic matter flux to the seafloor deriving from surface water productivity. The Sphaeroidina bulloides assemblage accounts for nearly half of the total variance, consisting of S. bulloides, the agglutinated taxa and Lenticulina as major components and Cibicidoides robertsonianus, Planulina ariminensis and Uvigerina peregrina parva as subsidiary species. Cibicidoides wuellerstorfi and Bulimina mexicana are important accessory species of the Cibicidoides pachyderma-dominated assemblage. The miliolids assemblage is dominated by miliolids and B. mexicana with Uvigerina mediterranea as the important associated species. The Praeglobobulimina ovata assemblage is represented by P. ovata along with Globobulimina affinis, Melonis pompiloides, Melonis barleeanum and U. peregrina parva as important accessory species. High productivity often occurs during the interglacial stages and is characterized by the S. bulloides assemblage with subsidiary P. ovata assemblage, while low productivity often occurs during the glacial stages and is marked by the miliolids assemblage and the C. pachyderma assemblage. In addition to the distinct glacial/interglacial variations, there is a long-term decreasing trend in primary productivity recorded by benthic foraminifera, which is supported by a broadly increasing trend in the sediment redness (a*). Continuous wavelet analysis of high productivity-related Factor 1 loading (the S. bulloides assemblage) suggests that variations in primary productivity have a middle Pleistocene transition (MPT)-related shift in cyclicity from the ~41-kyr dominant cycle prior to the MPT to the ~100-kyr dominant cycle after the MPT. Cross-spectral analyses between the Factor 1 loading record and δ18O of Globigerina bulloides reveal the control of global ice volume changes on primary productivity. The high coherency between the Factor 1 loading record and δ13C of G. bulloides on the glacial/interglacial timescale suggests that primary productivity off the Portuguese margin may be largely affected by nutrient availability which was related to the Iberian surface hydrographic changes associated with changes in global ice volume.

Impacts of the Middle Miocene Climatic Transition on deep-water carbonate preservation and oxygenation in the equatorial Indian Ocean

We reconstructed the orbital-scale variability of deep-water mass properties in the equatorial Indian Ocean during the late early to late Miocene (17.9 to 8.2 Ma) using benthic foraminiferal oxygen (δ18O) and carbon (δ13C) isotope records, in combination with X-ray fluorescence scanner elemental records at International Ocean Discovery Program (IODP) Site U1443. These records suggest that hyperthermal-like features (negative excursions in δ18O and δ13C coupled to suboxic and acidic conditions in the deep ocean) developed at eccentricity maxima during the Miocene Climatic Optimum (MCO; 16.9 to 14.7 Ma) on an almost ice-free Earth. Elevated Mn and U concentrations indicate that bottom and pore waters in the Indian Ocean became progressively suboxic during the later phase of the MCO. A fundamental change in deep-water circulation occurred after the expansion of the East Antarctic Ice Sheet during the Middle Miocene Climatic Transition (MMCT; 14.7 to ~13.8 Ma). Stepwise increases in δ18O at ~13.8 and ~13.1 Ma concurred with dampening of deep-water δ18O and δ13C variability, as Antarctic ice cover expanded and became more stable. A marked improvement in deep-water oxygenation after ~13.5 Ma coincided with the end of the last δ13C maximum of the Monterey Excursion. Carbonate dissolution in the equatorial Indian Ocean intensified and prolonged episodes of reduced carbonate deposition during the Carbonate Crash interval (~12.6 to 8.6 Ma) were not restricted to eccentricity maxima as during the MCO. Rising Mn and U concentrations after ~9 Ma indicate decreased oxygenation of bottom and pore waters, associated with the onset of the Biogenic Bloom in the Indian Ocean.

Dominant Role of Sea Level on the Sedimentary Environmental Evolution in the Bohai and Yellow Seas Over the Last 1 Million Years

Transgression and regression deposits from the shallow continental margin provide information on orbital-scale variations in sea level, climate change, and local tectonics. In this study, we conduct a high-resolution chronological and sedimentological analysis of a 125-m core (NHH01) drilled in the southern Yellow Sea. We developed a high-resolution age model at the orbital timescale over the last ∼1 Myr by the astronomical tuning of the sediment grain size. Sedimentological analysis and environmental proxies reveal that the study area is characterized by cyclic alternations of neritic and littoral/fluvial deposits controlled by glacial–interglacial sea-level changes. These results confirm the earlier assumption that sea-level fluctuations play a dominant role in the sedimentary architecture of the southern Yellow Sea. Moreover, only low-frequency sea-level fluctuations (∼100 kyr) were preserved in the NHH01 sequence; however, additional high-frequency (∼40 kyr) sea-level variations were also present in the sediments of the shallower Bohai Sea. Despite the large spatial difference, this finding implies that the sedimentary environment in the eastern marginal seas of China had been influenced by the sea level as a whole over the last 1 Myr. The comparison of the sedimentary environment with other cores from the Bohai Sea, Yellow Sea, and coastal region, as well as records from Indonesia and Japan, reveals that a consistent initiation of large-scale marine transgression occurred at ∼0.8–1 Ma. This phenomenon was likely caused by the long duration of glacio-eustatic high sea-level stands following the Middle Pleistocene transition, which is characterized by a shift in glacial cycles from 40 to 100 kyr.

Orbital-scale Asian summer monsoon variations: Paradox and exploration

The Asian summer monsoon (ASM) is a vast climate system, whose variability is critical to the livelihoods of billions of people across the Asian continent. During the past half-century, much progress has been made in understanding variations on a wide range of timescales, yet several significant issues remain unresolved. Of note are two long-standing problems concerning orbital-scale variations of the ASM. (1) Chinese loess magnetic susceptibility records show a persistent glacial-interglacial dominated ~100 kyr (thousand years) periodicity, while the cave oxygen-isotope (δ18O) records reveal periodicity in an almost pure precession band (~20 kyr periodicity)—the “Chinese 100 kyr problem”. (2) ASM records from the Arabian Sea and other oceans surrounding the Asian continent show a significant lag of 8–10 kyr to Northern Hemisphere summer insolation (NHSI), whereas the Asian cave δ18O records follow NHSI without a significant lag—a discrepancy termed the “sea-land precession-phase paradox”. How can we reconcile these differences? Recent and more refined model simulations now provide spatial patterns of rainfall and wind across the precession cycle, revealing distinct regional divergences in the ASM domain, which can well explain a large portion of the disparities between the loess, marine, and cave proxy records. Overall, we also find that the loess, marine, and cave records are indeed complementary rather than incompatible, with each record preferentially describing a certain aspect of ASM dynamics. Our study provides new insight into the understanding of different hydroclimatic proxies and largely reconciles the “Chinese 100 kyr problem“ and “sea-land precession-phase paradox”.