Search

The Indonesian Throughflow (ITF) controls the oceanic flux of heat and salt between the Pacific and Indian Oceans and therewith plays an important role in modulating the meridional overturning circulation and low latitude hydrological cycle. Here, we report new sea surface temperature and aridity records from the west coast of Australia (IODP Site U1460), which allow us to assess the sensitivity of the eastern Indian Ocean to the major reorganization of Earth’s climate that occurred during the Mid-Pleistocene Transition. Our records indicate glacial coolings at 1.55 and 0.65 million years ago that are best explained by a weakening of the ITF as a consequence of global sea level and tectonic changes. These coincide with the development of pronounced gradients in the carbon isotope composition of the different ocean basins and with substantial changes in regional aridity, suggesting that the restrictions of the ITF influenced both the evolution of global ocean circulation and the development of the modern hydrological cycle in Western Australia., Scientific Reports, 9 (1), ISSN:2045-2322

The International Ocean Discovery Programme (IODP) and its predecessors generated a treasure trove of Cenozoic climate and carbon cycle dynamics. Yet, it remains unclear how climate and carbon cycle interacted under changing geologic boundary conditions. Here, we present the carbon isotope (δ13C) megasplice, documenting deep-ocean δ13C evolution since 35 million years ago (Ma). We juxtapose the δ13C megasplice with its δ18O counterpart and determine their phase-difference on ~100-kyr eccentricity timescales. This analysis reveals that 2.4-Myr eccentricity cycles modulate the δ13C-δ18O phase relationship throughout the Oligo-Miocene (34-6 Ma), potentially through changes in continental weathering. At 6 Ma, a striking switch from in-phase to anti-phase behaviour occurs, signalling a reorganization of the climate-carbon cycle system. We hypothesize that this transition is consistent with Arctic cooling: Prior to 6 Ma, low-latitude continental carbon reservoirs expanded during astronomically-forced cool spells. After 6 Ma, however, continental carbon reservoirs contract rather than expand during cold periods due to competing effects between Arctic biomes (ice, tundra, taiga). We conclude that, on geologic timescales, System Earth experienced state-dependent modes of climate–carbon cycle interaction., Climate and carbon cycle interactions during major Earth system changes through the Cenozoic remain unclear. Here, the authors present a combined δ13C-δ18O megasplice for the last 35 Ma which allows them to identify three marked intervals of distinct climate–carbon cycle interactions.

Ocean gateways facilitate circulation between ocean basins, thereby impacting global climate. The Indonesian Gateway transports water from the Pacific to the Indian Ocean via the Indonesian Throughflow (ITF) and drives the strength and intensity of the modern Leeuwin Current, which carries warm equatorial waters along the western coast of Australia to higher latitudes. Therefore, ITF dynamics are a vital component of global thermohaline circulation. Plio‐Pleistocene changes in ITF behavior and Leeuwin Current intensity remain poorly constrained due to a lack of sedimentary records from regions under its influence. Here, organic geochemical proxies are used to reconstruct sea surface temperatures on the northwest Australian shelf at IODP Site U1463, downstream of the ITF outlet and under the influence of the Leeuwin Current. Our records, based on TEX86 and the long‐chain diol index, provide insight into past ITF variability (3.5–1.5 Ma) and confirm that sea surface temperature exerted a control on Australian continental hydroclimate. A significant TEX86 cooling of ~5°C occurs within the mid‐Pliocene Warm Period (3.3–3.1 Ma) suggesting that this interval was characterized by SST fluctuations at Site U1463. A major feature of both the TEX86 and long‐chain diol index records is a strong cooling from ~1.7 to 1.5 Ma. We suggest that this event reflects a reduction in Leeuwin Current intensity due to a major step in ongoing ITF constriction, accompanied by a switch from South to North Pacific source waters entering the ITF inlet. Our new data suggest that an additional ITF constriction event may have occurred in the Pleistocene.

Pre-Quaternary paleoclimate studies in Australia mainly focus on terrestrial records from the southeastern part of the continent. IODP Expedition 356 drilled on the northwestern Australian shelf, yielding Miocene-Pleistocene paleoclimate records in an area where climate archives are scarce. Postexpedition research revealed a dry-to-humid transition across the latest Miocene and early Pliocene (start of the "Humid Interval"). However, the complex tectonic history of the area makes these interpretations challenging. In this study, we investigate late Miocene to early Pliocene sediment cores from two sites that are only 100 km apart but situated in two adjacent basins (Northern Carnarvon and Roebuck Basins). Combining lithofacies study, time series analysis of potassium content (K wt%), and calcareous nannofossil abundance counts (N/g), this work disentangles the complex interplay between basin evolution and climate change between 6.1 and 4 Ma. Overall, the investigated proxies show high correlation between both sites, except during 6.1-5.7 Ma. During this interval, Site U1463 records a gradual increase in K wt%, correlated with basin deepening, whereas Site U1464 records an abrupt rise in K wt% at similar to 6 Ma. We explain this diachronicity by differential basin subsidence. The tectonic interplay with our paleorecords makes it difficult to pinpoint the exact onset of the "Humid Interval," but we conclude that K wt% and coccolith abundances at Site U1464 indicate that a fluvial deposition system was already established since at least 6 Ma. This age is consistent with data supporting a southward movement of the Intertropical Convergence Zone rain belt at similar to 7 Ma.

International audience; The geologic time scale for the Cenozoic Era has been notably improved over the last decades by virtue of integrated stratigraphy, combining high-resolution astrochronologies, biostratigraphy and magnetostratigraphy with high-precision radioisotopic dates. However, the middle Eocene remains a weak link. The so-called “Eocene time scale gap” reflects the scarcity of suitable study sections with clear astronomically-forced variations in carbonate content, primarily because large parts of the oceans were starved of carbonate during the Eocene greenhouse. International Ocean Discovery Program (IODP) Expedition 369 cored a carbonate-rich sedimentary sequence of Eocene age in the Mentelle Basin (Site U1514, offshore southwest Australia). The sequence consists of nannofossil chalk and exhibits rhythmic clay content variability. Here, we show that IODP Site U1514 allows for the extraction of an astronomical signal and the construction of an Eocene astrochronology, using 3-cm resolution X-Ray fluorescence (XRF) core scans. The XRF-derived ratio between calcium and iron content (Ca/Fe) tracks the lithologic variability and serves as the basis for our U1514 astrochronology. We present a 16 million-year-long (40-56 Ma) nearly continuous history of Eocene sedimentation with variations paced by eccentricity and obliquity. We supplement the high-resolution XRF data with low-resolution bulk carbon and oxygen isotopes, recording the long-term cooling trend from the Paleocene-Eocene Thermal Maximum (PETM – ca. 56 Ma) into the middle Eocene (ca. 40 Ma). Our early Eocene astrochronology corroborates existing chronologies based on deep-sea sites and Italian land sections. For the middle Eocene, the sedimentological record at U1514 provides a single-site geochemical backbone and thus offers a further step towards a fully integrated Cenozoic geologic time scale at orbital resolution.

Much of our understanding of Earth's past climate comes from the measurement of oxygen and carbon isotope variations in deep-sea benthic foraminifera. Yet, long intervals in existing records lack the temporal resolution and age control needed to thoroughly categorize climate states of the Cenozoic era and to study their dynamics. Here, we present a new, highly resolved, astronomically dated, continuous composite of benthic foraminifer isotope records developed in our laboratories. Four climate states-Hothouse, Warmhouse, Coolhouse, Icehouse-are identified on the basis of their distinctive response to astronomical forcing depending on greenhouse gas concentrations and polar ice sheet volume. Statistical analysis of the nonlinear behavior encoded in our record reveals the key role that polar ice volume plays in the predictability of Cenozoic climate dynamics.

An unusually short glaciation interrupted the warm Pliocene around 3.3 Ma (Marine Isotope Stage (MIS) M2). Different hypotheses exist to explain why this glaciation event was so pronounced, and why the global climate system returned to warm Pliocene conditions relatively quickly afterwards. One of these proposed mechanisms is a reduced equator-to-pole heat transfer, in response to a tectonically reduced Indonesian Throughflow (ITF). The ITF is a critical part of the global thermohaline ocean circulation, transporting heat from the Indo-Pacific Warm Pool to the Indian Ocean. When ITF connectivity is reduced, the water and heat supply for the Leeuwin Current, flowing poleward along Australia's west coast, is also diminished. To assess the possible relationship between mid-Pliocene glaciations and latitudinal heat transport through the Indonesian Throughflow, we constructed a multi-proxy orbital-scale record for the 3.7–2.8 Ma interval from International Ocean Discovery Program (IODP) Site U1463, off northwest Australia. The comparison of the Site U1463 record with paleoclimate records from nearby Site 763 and West Pacific Warm Pool Site 806 allows for a detailed regional reconstruction of Pliocene paleoceanography and thus for testing the proposed hypothesis. An astronomically-paced decrease in potassium content characterizes the late Pliocene interval of U1463. This record documents the increasing aridity of northwest Australia, periodically alleviated by reinforced summer monsoon precipitation under summer insolation maxima. The δ 18 O record of the planktonic foraminifer Globigerinoides sacculifer correlates exceptionally well with the sea surface temperature (SST) record from Site 806 in the West Pacific Warm Pool, even during MIS M2. Hence, Site U1463 preserves an uninterrupted ITF signal even during Pliocene glaciations. However, the U1463 δ 18 O G.sacculifer record exhibits a 0.5‰ offset with the nearby Site 763A record around MIS M2. This implies that Site 763A, about 500 km west of U1463, more closely tracks Indian Ocean SST records across MIS M2. The U1463 data reveal that heat-transport through the Indonesian Throughflow did not shut down completely during MIS M2, but rather its intensity decreased prior to and during MIS M2, causing Site 763A to temporarily reflect an Indian Ocean, rather than an ITF signal. We conclude that ITF variability significantly influenced latitudinal heat transport by means of the Leeuwin Current and hence contributed to the relative intensity of MIS M2. We propose the ITF valve between the Pacific and Indian Ocean as a positive feedback mechanism, in which an initial sea level lowering reduces ITF heat transport, in turn amplifying global cooling by advancing the thermal isolation of Antarctica.

A high-resolution potassium content record of sediments from Integrated Ocean Drilling Program (IODP) Site U1422 drilled in the northern Sea of Japan was employed to yield an astronomical timescale for this ∼205 m thick sedimentary archive. The K content was quantified using the natural gamma radiation (NGR) data routinely measured on DV JOIDES Resolution. The U1422 K (wt.%) series shows an increasing trend with time, which parallels the growth of North Hemisphere ice sheets since ∼4 Ma, as revealed by the global benthic foraminifer's oxygen isotope stack (LR04). We propose that K content variations reflect changes in the relative contributions of Asian eolian dust and volcanic weathering products, in response to changes in global ice volume. Using the shipboard age model, constrained by palaeomagnetism and nanofossils datums, we tuned the U1422 K content variations to the LR04 stack. Our tuned age model for site U1422 goes back to 3.9 Ma. The studied record indicates that Asian dust generation is enhanced during periods of global cooling. The synchronous variations between the U1422 K record and the LR04 stack also suggest that global cooling played a dominant role in promoting the central Asian aridification since ∼3.9 Ma, while the uplift of Tibetan Plateau plays a secondary role, considering many disagreements still exists between the timing and amplitude of the Qinghai-Tibetan uplift and the evolution of central Asian aridity since the early Pliocene.

North Atlantic Deep Water (NADW) currently redistributes heat and salt between Earth's ocean basins, and plays a vital role in the ocean-atmosphere CO2 exchange. Despite its crucial role in today's climate system, vigorous debate remains as to when deep-water formation in the North Atlantic started. Here, we present datasets from carbonate-rich middle Eocene sediments from the Newfoundland Ridge, revealing a unique archive of paleoceanographic change from the progressively cooling climate of the middle Eocene. Well-defined lithologic alternations between calcareous ooze and clay-rich intervals occur at the ∼41-kyr beat of axial obliquity. Hence, we identify obliquity as the driver of middle Eocene (43.5–46 Ma) Northern Component Water (NCW, the predecessor of modern NADW) variability. High-resolution benthic foraminiferal δ 18 O and δ 13 C suggest that obliquity minima correspond to cold, nutrient-depleted, western North Atlantic deep waters. We thus link stronger NCW formation with obliquity minima. In contrast, during obliquity maxima, Deep Western Boundary Currents were weaker and warmer, while abyssal nutrients were more abundant. These aspects reflect a more sluggish NCW formation. This obliquity-paced paleoceanographic regime is in excellent agreement with results from an Earth system model, in which obliquity minima configurations enhance NCW formation.

The Middle Pleistocene Transition (MPT) represents a major change in Earth's climate state, exemplified by the switch from obliquity-dominated to ∼100-kyr glacial/interglacial cycles. To date, the causes of this significant change in Earth's climatic response to orbital forcing are not fully understood. Nonetheless, this transition represents an intrinsic shift in Earth's response to orbital forcing, without fundamental changes in the astronomical rhythms. This study presents new high-resolution records of International Ocean Discovery Program (IODP) Site U1460 (eastern Indian Ocean, 27°S), including shallow marine productivity and organic matter flux reconstructions. The proxy series covers the interval between 1.1 and 0.6 Ma and provides insights into Pleistocene Leeuwin Current dynamics along the West Australian shelf. The large >45 m global sea level drop during the marine isotope stage (MIS) 22–24 is marked in our data, suggesting that the MPT led to large-scale changes in Indian Ocean circulation patterns and surface water conditions. We consider shelf exposure (and thus the “Sahul-Indian Ocean Bjerknes mechanism”) as a possible key process to increase the upwelling of nutrient-rich sub-Antarctic Mode waters through the Leeuwin Undercurrent along the Australian shelf. We conclude that the shoaling of nutrient-rich lower-thermocline waters enhanced mid-latitude productivity patterns in the eastern Indian Ocean across the 900-ka event.

During International Ocean Discovery Program (IODP) expeditions, shipboard-generated data provide the first insights into the cored sequences. The natural gamma radiation (NGR) of the recovered material, for example, is routinely measured on the ocean drilling research vessel DV JOIDES Resolution. At present, only total NGR counts are readily available as shipboard data, although full NGR spectra (counts as a function of gamma-ray energy level) are produced and archived. These spectra contain unexploited information, as one can estimate the sedimentary contents of potassium (K), thorium (Th), and uranium (U) from the characteristic gamma-ray energies of isotopes in the ^(40)K, ^(232)Th, and ^(238)U radioactive decay series. Dunlea et al. [2013] quantified K, Th and U contents in sediment from the South Pacific Gyre by integrating counts over specific energy levels of the NGR spectrum. However, the algorithm used in their study is unavailable to the wider scientific community due to commercial proprietary reasons. Here, we present a new MATLAB algorithm for the quantification of NGR spectra that is transparent and accessible to future NGR users. We demonstrate the algorithm's performance by comparing its results to shore-based inductively coupled plasma-mass spectrometry (ICP-MS), inductively coupled plasma-emission spectrometry (ICP-ES), and quantitative wavelength-dispersive X-ray fluorescence (XRF) analyses. Samples for these comparisons come from eleven sites (U1341, U1343, U1366-U1369, U1414, U1428-U1430, U1463) cored in two oceans during five expeditions. In short, our algorithm rapidly produces detailed high-quality information on sediment properties during IODP expeditions at no extra cost.

In this paper we formulate answers to three important questions related to Givetian carbonate records and their use for reconstructing million-year-scale past palaeoenvironmental changes. First, we provide detailed illustrations of the fascinating diversity that shaped a significant shallow reefal platform during early to late Givetian time in the Rhenohercynian Ocean; secondly we improve the sedimentological model of the extensive Givetian carbonate platform in the Dinant Basin; and thirdly we evaluate the application of magnetic susceptibility as a tool for long-term trend correlations and palaeoenvironmental reconstructions. These goals are reached by making a sedimentological, geophysical and geochemical study of the La Thure section. Through the early–late Givetian interval we discerned 18 microfacies ranging from a homoclinal ramp to a discontinuously rimmed shelf and then a drowning shelf. The comparison of these sedimentological results with those published for the south of the Dinant Syncline allowed us to provide an up to date model of the vertical and lateral environmental development of one of the largest Givetian carbonate platforms in Europe. This comparison also increased the knowledge on the distribution of facies belts in the Dinant Basin and allowed us to highlight the Taghanic Event. Palaeoredox proxies reveal a substantial change in the oxygenation level, from oxygen-depleted to more oxic conditions, between middle and late Givetian time. We demonstrated the relationship between variation in magnetic susceptibility values and proxies for siliciclastic input (such as Si, Al). The La Thure section is considered a key section for the understanding of internal shelf settings bordering Laurussia9s southeastern margin.

A detailed paleoenvironment reconstruction from the Mozambique Channel, western Indian Ocean, based on the calcareous nannoplankton assemblages was conducted for the interval between 2.85 and 1.85 Myr. This study covers the period during which the successive extinction of the last five species of discoasters occurred. New productivity data obtained from the abundances of the Discoaster species (Discoaster brouweri, D. triradiatus, D. pentaradiatus, D. surculus, and D. tamalis) and other indicative calcareous nannoplankton taxa showed abundance variations, which were at paced with the 100, 41, and 23 kyr astronomical periodicities. A shift in the productivity and water-column stratification proxies occurred at ~2.4 Ma, after the onset of the Northern Hemisphere glaciation. Here we propose that the variability recorded at International Ocean Discovery Program Site U1476 reflects the interplay between forcing associated with warm tropical Pacific and cold southern ocean influences. The former is shown by consistent occurrence of warm water taxa (Calcidiscus leptoporus, Oolithotus spp., Rhabdosphaera clavigera, Syracosphaera spp., Umbellosphaera spp.), typical of Indonesian Throughflow surface waters. On the other hand, the occurrence of Coccolithus pelagicus indicates the influence of cold, nutrient-rich sub-Antarctic surface waters. A more mixed water column initiated at ~2.4 Ma, and a consequent productivity increase led to the gradual reduction of the Discoaster species, until their extinction at 1.91 Ma. This period was characterized by the low values of the Florisphaera profunda index and high abundances of upper photic zone flora, indicative of nutrient-rich surface water conditions. High productivity at the location during this period could have also been amplified by localized upwelling events driven by the Mozambique Channel eddies.

Around the early–middle Eocene boundary, the first occurrence of contourite drift sediments and widespread deep ocean erosion indicate changes in the North Atlantic ocean circulation. Interestingly, these changes coincide with the first steps of Cenozoic cooling from the Paleogene greenhouse climate towards the modern icehouse. The cause for this ocean circulation reorganization is poorly understood since modern water mass tracers may have worked fundamentally different in the past and the paleoceanographic proxy record is limited in both time and space. As a result, it is challenging to reliably reconstruct the climatic and tectonic boundary conditions e.g. atmospheric greenhouse gas concentration and the depth and geometry of developing and closing passages between ocean basins. In this study, we attempt to identify thresholds in tectonic gateway passages and atmospheric CO2 concentration, using the fully coupled Earth System Model COSMOS. Indeed, the simulation of Earth's past climates can unravel the physical processes driving deep-water formation in a greenhouse world. Specifically, we use COSMOS to evaluate the impact of changes in the North Atlantic gateways at the early–middle Eocene boundary on the North Atlantic Deep Western Boundary Currents under low obliquity configuration. We find that Northern Component Waters start to form when the Greenland Scotland Ridge reaches a threshold depth of deeper than 200 m, while the Arctic Ocean is still shut off from the North Atlantic. In this scenario, the relatively deep Greenland Scotland Ridge allows for sufficient inflow of warm, salty Atlantic surface waters into the Nordic Seas to initiate convection during winter cooling. Opening the seaway towards the Arctic leads to a cessation of Northern Component Water formation as it allows for inflow of brackish surface waters into the northern Nordic Seas, hindering Northern Component Water formation.

The Frasnian–Famennian boundary records one of the most catastrophic mass extinctions of the Phanerozoic Eon. Several possible causes for this extinction have been suggested, including extra-terrestrial impacts and large-scale volcanism. However, linking the extinction with these potential causes is hindered by the lack of precise dating of either the extinction or volcanic/impact events. In this study, a bentonite layer in uppermost-Frasnian sediments from Steinbruch Schmidt (Germany) is re-analysed using CA-ID-TIMS U-Pb zircon geochronology in order to constrain the date of the Frasnian–Famennian extinction. A new age of 372.36 ± 0.053 Ma is determined for this bentonite, confirming a date no older than 372.4 Ma for the Frasnian–Famennian boundary, which can be further constrained to 371.93–371.78 Ma using a pre-existing Late Devonian age model. This age is consistent with previous dates, but is significantly more precise. When compared with published ages of the Siljan impact crater and basalts produced by large-scale volcanism, there is no apparent correlation between the extinction and either phenomenon, not clearly supporting them as a direct cause for the Frasnian–Famennian event. This result highlights an urgent need for further Late Devonian geochronological and chemostratigraphic work to better understand the cause(s) of this extinction.

Cyclostratigraphic analyses rely on techniques that trace astronomical components in paleoclimate signals. These techniques have demonstrated their value but rely on certain assumptions on the presence and quality of the astronomical imprint. Here we explore a new conceptual approach to time series analysis. Specifically, we evaluate the potential of spectral moments to characterize the full spectral characteristics of a record and thus not only the frequency ranges of interpreted astronomical components. Mathematically speaking, moments are unique quantities describing a specific set of points. In the case of spectral moments, we apply the concept of moments on the distribution of spectral power in a signal's periodogram. We present four case studies that illustrate the advantages and disadvantages of the spectral moment approach in gaining insight in the (astronomical) features of a particular data record. We discuss the effects of outliers in a series, variable sedimentation rate, and changing climate dynamics on the spectral moments of a power spectrum. At the same time, we carry out a sedimentation rate reconstruction based on the spectral moment approach and compare that reconstruction to results obtained through classic approaches. Compared to classic approaches, the spectral moments approach is robust and requires less strict assumptions to obtain similar first‐order information. Yet such assumptions are still necessary to achieve more precise sedimentation rate reconstructions. In summary, we show that the spectral moments approach is suitable to obtain first‐order insights in variable components embedded in a depth series.

The Wetteldorf section is the Global Stratotype Section and Point (GSSP) locality of the Lower–Middle Devonian boundary. The section is characterized by an alternation of marls, limestones, and siltstones, with the base of the Eifelian stage defined by the first occurrence of the conodont Polygnathus costatus partitus. The well-established conodont biostratigraphy for this section permits correlations with other time-equivalent sections; however, a cyclostratigraphic framework does not yet exist for this section. Here, we construct such a framework that provides an additional tool for long-distance correlation, and at the same time constrains the amount of geologic time represented by this stratigraphic interval. Microfacies analysis indicates that the entire section has been deposited in an open-marine shelf paleoenvironment, below the storm wave base and under oxygen-rich conditions. We collected high-resolution (5-cm spaced) geochemical data, using a handheld X-ray fluorescence (XRF) spectrometer, and quantified net intensity results for 23 chemical elements. By using factor analysis, we explain 52% of the observed variance in those elements by a model that consists of four factors. The most important factor, factor 1, is driven by geochemical elements related to detrital input (e. g. K, Rb, Si, Ti). Factor 2 is interpreted as a pyrite-indicator, while Factor 3 is driven by Ca and Sr and is thus tied to lithological variations in carbonate content. Factor 1 and 3 are used for cyclo-stratigraphy, as their respective power spectra show similar characteristics. Both factors exhibit a strong 58-cm period, which can be interpreted, according to the available biostratigraphic constraints, either as the imprint of the 18-kyr precession cycle or as the imprint of the 33-kyr obliquity cycle. However, the ratio between different spectral peaks in the different proxies’ power spectra strongly advocates the first option. Our interpretation in terms of astronomical climate forcing implies that the 8.65 meter thick Wetteldorf section (within the Happel Hut) is characterized by an average sedimentation rate of ∼ 3.2 cm/kyr and represents ∼ 250 kyr. From a paleoclimatological point of view, our results indicate that detrital input into the Rhenohercynian Basin was spurred under high-eccentricity configurations, which allowed for strongly enhanced seasonal rainfall. Under astronomical configurations that avoided extremes in seasonality for a prolonged period of time (simultaneously low eccentricity and obliquity), Factor 2 indicates the occurrence of pyrite, which we associate with temporary deteriorations of generally well-oxygenated conditions.

International audience; Until now, the middle Eocene has remained a poorly constrained interval of efforts to produce an astrochronological timescale for the entire Cenozoic. This has given rise to a so-called "Eocene astronomical timescale gap" (Vandenberghe et al., 2012). A high-resolution astrochronological calibration for this interval has proven to be difficult to realize, mainly because carbonate-rich deep-marine sequences of this age are scarce. In this paper, we present records from middle Eocene carbonate-rich sequences from the North Atlantic Southeast Newfoundland Ridge (IODP Exp. 342, Sites U1408 and U1410), of which the cyclical sedimentary patterns allow for an orbital calibration of the geologic timescale between ∼38 and ∼48 Ma. These carbonate-rich cyclic sediments at Sites U1408 and U1410 were deposited as drift deposits and exhibit prominent lithological alternations (couplets) between greenish nannofossil-rich clay and white nannofossil ooze. The principal lithological couplet is driven by the obliquity of Earth's axial tilt, and the intensity of their expression is modulated by a cyclicity of about 173 kyr. This cyclicity corresponds to the interference of secular frequencies s3 and s6 (related to the precession of nodes of the Earth and Saturn, respectively). This 173-kyr obliquity amplitude modulation cycle is exceptionally well recorded in the XRF (X-ray fluorescence)-derived Ca/Fe ratio. In this work, we first demonstrate the stability of the (s3-s6) cycles using the latest astronomical solutions. Results show that this orbital component is stable back to at least 50 Ma, and can thus serve as a powerful geochronometer in the mid-Eocene portion of the Cenozoic timescale. We then exploit this potential by calibrating the geochronology of the recovered middle Eocene timescale between magnetic polarity Chrons C18n.1n and C21n. Comparison with previous timescales shows similarities, but also notable differences in durations of certain magnetic polarity chrons. We present a revision of previous astronomical timescales from the Equatorial and South Atlantic, to overcome the differences between different mid-Eocene astrochronologies. Using our new records from the North Atlantic, combined with existing records from the South Atlantic (ODP Site 1263 and Hole 702B) and Equatorial Atlantic (ODP Site 1260), we revise the durations of magnetic polarity Chrons C18n.1n to C21n, thereby arriving at a robust and self-consistent closure of the middle Eocene astronomical timescale gap.

International audience; The Japan Sea is directly influenced by the Asian monsoon, a system that transports moisture and heat across southeast Asia during the boreal summer, and is a major driver of the Earth's ocean-atmospheric circulation. Foraminiferal and facies analyses of a 460-kyr record from Integrated Ocean Drilling Program Expedition 346 Site U1427 in the Japan Sea reveal a record of nutrient flux and oxygenation that varied due to sea level and East Asian monsoon intensity. The East Asian summer monsoon (EASM) was most intense during marine isotope stage (MIS) 5e, MIS 7e, MIS 9e, and MIS 11c when the Tsushima Warm Current flowed into an unrestricted well-mixed normal salinity Japan Sea, whereas East Asian winter monsoon (EAWM) conditions dominated MIS 2, MIS 4, MIS 6, and MIS 8 when sea level minima restricted the Japan Sea resulting in low-salinity and oxygen conditions in the absence of Tsushima flow. Reduced oxygen stratified, low-salinity, and higher productivity oceanic conditions characterize Terminations TV, TIII, TII, and TI when East China Sea coastal waters breached the Tsushima Strait. Chinese loess, cave, and Lake Biwa (Japan) and U1427 proxy records suggest EASM intensification during low to high insolation transitions, whereas the strongest EAWM prevailed during lowest insolation periods or high to low insolation transitions. Ice sheet/CO 2 forcing leads to the strongest EAWM events in glacials and enhanced EASM in interglacials. Mismatches between proxy patterns suggest that latitudinal and land/sea thermal contrasts played a role in East Asian monsoon variability, suggesting that a complex interplay between ice sheet dynamics, insolation, and thermal gradients controls monsoonal intensity.

The Late Devonian envelops one of Earth’s big five mass extinction events at the Frasnian–Famennian boundary (374 Ma). Environmental change across the extinction severely affected Devonian reef-builders, besides many other forms of marine life. Yet, cause-and-effect chains leading to the extinction remain poorly constrained as Late Devonian stratigraphy is poorly resolved, compared to younger cataclysmic intervals. In this study we present a global orbitally calibrated chronology across this momentous interval, applying cyclostratigraphic techniques. Our timescale stipulates that 600 kyr separate the lower and upper Kellwasser positive δ13C excursions. The latter excursion is paced by obliquity and is therein similar to Mesozoic intervals of environmental upheaval, like the Cretaceous Ocean-Anoxic-Event-2 (OAE-2). This obliquity signature implies coincidence with a minimum of the 2.4 Myr eccentricity cycle, during which obliquity prevails over precession, and highlights the decisive role of astronomically forced “Milankovitch” climate change in timing and pacing the Late Devonian mass extinction., Understanding of Late Devonian mass extinction mechanisms is poor due to imprecise stratigraphies. Here, using cyclostratigraphic techniques, the authors present a global orbitally-calibrated chronology and reveal the key role of astronomically-forced Milankovitch climate change.

The discovery of elevated concentrations of the cosmogenic radionuclide 3He in deep- sea sediments from Ocean Drilling Program (ODP) Site 926 (Atlantic Ocean) and ODP Site 757 (Indian Ocean) points toward ac- cretion of extraterrestrial matter, probably as a result of the catastrophic disruption of a large asteroid that produced the Veritas fam- ily of asteroids at ca. 8.3 ± 0.5 Ma, and which may have had important effects on the global climatic and ecologic systems. Here, we in- vestigated the signatures possibly related to the Veritas event by performing a high- resolution multiproxy stratigraphic analysis through the late Tortonian–early Messinian Monte dei Corvi section near Ancona, Italy. Closely spaced bulk-rock samples through a 36-m-thick section, approximately spanning from ca. 9.9 Ma to ca. 6.4 Ma, show an ~5- fold 3He anomaly starting at ca. 8.5 Ma and returning to background values at ca. 6.9 Ma, con rming the global nature of the event. We then analyzed, at 5 cm intervals, bulk-rock samples for sedimentary and environmen- tal proxies such as magnetic susceptibility, calcium carbonate content, total organic carbon, and bulk carbonate d18O and d13C, through a 21-m-thick section encompassing the 3He anomaly. Available high-resolution sea-surface temperature data (via alkenone analyses) for this site show a temperature decrease starting exactly at the inception of the 3He anomaly. Cyclostratigraphic fast- Fourier-transform spectral analyses of the proxies indicate an age of 8.47 ± 0.05 Ma for the inception of the 3He anomaly. A search for impact ejecta (analogous to what is pres- ent in the late Eocene, where both a 3He anomaly and large-scale impact events are recorded) was not successful. Detailed cy- clostratigraphic analyses of our data suggest that the changes in the stable isotope series and environmental proxy series through this late Tortonian time interval had a common forcing agent, and that perturbations of or- bitally forced climate cycles are present ex- actly through the interval with the enhanced in ux of extraterrestrial 3He. Thus, the che- mostratigraphic evidence for a collisional event that created the Veritas family of aster- oids, coinciding with climate perturbations on Earth, suggests yet another form of inter- action between Earth and the solar system.

Geophysical research letters 44(13), 6914-6925(2017). doi:10.1002/2017GL072977, Published by Wiley, Hoboken, NJ

Earth’s climate has undergone different intervals of gradual change as well as abrupt shifts between climate states. Here we aim to characterize the corresponding changes in climate response to astronomical forcing in the icehouse portion of the Cenozoic, from the latest Eocene to the present. As a tool, we use a 35-m.y.-long d18Obenthic record compiled from different high-resolution benthic isotope records spliced together (what we refer to as a megasplice).We analyze the climate response to astronomical forcing during four 800-k.y.- long time windows. During the mid-Miocene Climatic Optimum (ca. 15.5 Ma), global climate variability was mainly dependent on Southern Hemisphere summer insolation, ampli ed by a dynamic Antarctic ice sheet; 2.5 m.y. later, relatively warm global climate states occurred during maxima in both Southern Hemisphere and Northern Hemisphere summer insolation. At that point, the Antarctic ice sheet grew too big to pulse on the beat of precession, and the Southern Hemisphere lost its overwhelming in uence on the global climate state. Likewise, we juxtapose response regimes of the Miocene (ca. 19 Ma) and Oligocene (ca. 25.5 Ma) warming periods. Despite the similarity in d18Obenthic values and variability, we nd different responses to precession forcing. While Miocene warmth occurs during summer insolation maxima in both hemispheres, Oligocene global warmth is consistently triggered when Earth reaches perihelion in the Northern Hemisphere summer. This pattern is in accordance with previously published paleoclimate modeling results, and suggests an amplifying role for Northern Hemisphere sea ice.

The geological record of the Paleozoic often exhibits cyclic features, in many cases the result of changes in paleoclimate. However, a thorough understanding of the processes that were driving Paleozoic climate change has not yet been reached. The main reason is relatively poor time-control on Paleozoic paleoclimate proxy records. This problem can be overcome by the identification of cyclic features resulting from astronomical climate forcing in the stratigraphic record. To correctly identify these cyclic features, it is necessary to quantify the effects of astronomical climate forcing under conditions different from today. In this work, we apply Late Devonian (375 Ma) boundary conditions to the Hadley Centre general circulation model (HadSM3). We estimate the response of Late Devonian climate to astronomical forcing by keeping all other forcing factors (e.g. paleogeography, p CO2, vegetation distribution) fixed. Thirty-one different "snapshots" of Late Devonian climate are simulated, by running the model with different combinations of eccentricity (e), obliquity (?) and precession. From the comparison of these 31 simulations, it appears that feedback mechanisms play an important role in the conversion of astronomically driven insolation variations into climate change, such as the formation of sea-ice and the development of an extensive snow cover on Gondwana. We compare the "median orbit" simulation to lithic indicators of paleoclimate to evaluate whether or not HadSM3 validly simulates Late Devonian climates. This comparison suggests that the model correctly locates the major climate zones. This study also tests the proposed link between the formation of ocean anoxia and high eccentricity (De Vleeschouwer et al., 2013) by comparing the ?18Ocarb record of the Frasnian-Famennian boundary interval from the Kowala section (Poland) with a simulated time series of astronomically forced changes in mean annual temperature at the paleolocation of Poland. The amplitude of climate change suggested by the isotope record is greater than that of the simulated climate. Hence, astronomically forced climate change may have been further amplified by other feedback mechanisms not considered here (e.g. CO2 and vegetation). Finally, the geologic and simulated time series correlate best when the Frasnian-Famennian negative isotope excursion aligns with maximum mean annual temperature in Poland, which is obtained when eccentricity and obliquity are simultaneously high. This finding supports a connection between Devonian ocean anoxic events and astronomical climate forcing.

The International Ocean Discovery Program (IODP) Expedition 356 Site U1461 cored a Miocene to Holocene sedimentary sequence in the upper bathyal carbonate offshore northwestern Australia (NWA). The siliciclastic component of these strata is primarily derived from the Australian continent. Radiocarbon dating on macrofossils and planktonic foraminifera shows that the upper 14 m section at Site U1461 preserves Holocene sediments, recording regional climate variability. K/Ca ratios determined by X-ray fluorescence elemental analyses and %K determined by shipboard natural gamma ray analysis are interpreted as indicators of riverine run-off from the Australian continent. We document the consequences of the variability of the Australian Summer Monsoon (ASM) on the continental shelf of NWA. We report an increase in terrigenous input due to a riverine run-off after 11.5 ka, which reaches a maximum at similar to 8.5 ka. This maximum is the result of the enhanced ASM-derived precipitation in response to the southern migration of the Intertropical Convergence Zone (ITCZ). A decrease in riverine run-off due to a weakening of precipitation in the NWA region after 8.5 ka was caused by the northern migration of the ITCZ. We conclude that the ITCZ reached its southernmost position at 8.5 ka and enhanced precipitation in the NWA region. This Holocene record shows that even during interglacial periods, monsoonal variability was primarily controlled by the position of the ITCZ.

Accurately deciphering periodic variations in paleoclimate proxy signals is essential for cyclostratigraphy. Classical spectral analysis often relies on methods based on (fast) Fourier transformation. This technique has no unique solution separating variations in amplitude and frequency. This characteristic can make it difficult to correctly interpret a proxy's power spectrum or to accurately evaluate simultaneous changes in amplitude and frequency in evolutionary analyses. This drawback is circumvented by using a polynomial approach to estimate instantaneous amplitude and frequency in orbital components. This approach was proven useful to characterize audio signals (music and speech), which are non-stationary in nature. Paleoclimate proxy signals and audio signals share similar dynamics; the only difference is the frequency relationship between the different components. A harmonic-frequency relationship exists in audio signals, whereas this relation is non-harmonic in paleoclimate signals. However, this difference is irrelevant for the problem of separating simultaneous changes in amplitude and frequency. Using an approach with overlapping analysis frames, the model (Astronomical Component Estimation, version 1: ACE v.1) captures time variations of an orbital component by modulating a stationary sinusoid centered at its mean frequency, with a single polynomial. Hence, the parameters that determine the model are the mean frequency of the orbital component and the polynomial coefficients. The first parameter depends on geologic interpretations, whereas the latter are estimated by means of linear least-squares. As output, the model provides the orbital component waveform, either in the depth or time domain. Uncertainty analyses of the model estimates are performed using Monte Carlo simulations. Furthermore, it allows for a unique decomposition of the signal into its instantaneous amplitude and frequency. Frequency modulation patterns reconstruct changes in accumulation rate, whereas amplitude modulation identifies eccentricity-modulated precession. The functioning of the time-variant sinusoidal model is illustrated and validated using a synthetic insolation signal. The new modeling approach is tested on two case studies: (1) a Pliocene–Pleistocene benthic δ18O record from Ocean Drilling Program (ODP) Site 846 and (2) a Danian magnetic susceptibility record from the Contessa Highway section, Gubbio, Italy. Matthias Sinnesael thanks the Research Foundation – Flanders (FWO) for the awarded PhD fellowship (FWOTM782). David De Vleeschouwer was funded through European Research Council (ERC) Consolidator Grant “Earthsequencing” (grant agreement no. 617462). This work was supported in part by the Fund for Scientific Research (FWO-Vlaanderen grant G009113N and support to Johan Schoukens), by the Flemish Government (Methusalem METH1 to Johan Schoukens), and by the Belgian Government through the Inter University Pole of Attraction (IUAP VII) Program (P7/15 PLANET TOPERS and P7/19 DYSCO).

The regional patterns and timing of the Younger Dryas cooling in the North Atlantic realm were complex and are mechanistically incompletely understood. To enhance understanding of regional climate patterns, we present molecular biomarker records at subannual to annual resolution by mass spectrometry imaging (MSI) of sediments from the Lake Meerfelder Maar covering the Allerød-Younger Dryas transition. These analyses are supported by conventional extraction-based molecular-isotopic analyses, which both validate the imaging results and constrain the sources of the target compounds. The targeted fatty acid biomarkers serve as a gauge of the response of the local aquatic and terrestrial ecosystem to climate change. Based on the comparison of our data with existing data from Meerfelder Maar, we analyse the short-term environmental evolution in Western Europe during the studied time interval and confirm the previously reported delayed hydrological response to Greenland cooling. However, despite a detected delay of Western European environmental change of ∼135 years, our biomarker data show statistically significant correlation with deuterium excess in Greenland ice core at ∼ annual resolution during this time-transgressive cooling. This suggests a coherent atmospheric forcing across the North Atlantic realm during this transition. We propose that Western European cooling was postponed due to major reorganization of the westerlies that were intermittently forcing warmer and wetter air masses from lower latitudes to Western Europe and thus resulted in delayed cooling relative to Greenland.