Site U1313 Research Articles

Millennial‐Scale Climate Variability Potentially Shaped the Early Interglacial Optimum in Southern Europe

AbstractThe seasonal and latitudinal distribution of insolation is considered the main factor controlling the magnitude and timing of interglacial periods. However, despite small differences in insolation forcing, vegetation and hydrology in southern Europe during past interglacials are variable and the gradual change in insolation cannot explain the observed short‐lived forest optimum. Here we focus on vegetation and hydroclimatic changes at orbital‐ and suborbital‐scales in southwestern Europe during two past warm interglacial periods with reduced ice‐sheets, namely Marine Isotope Stages (MIS) 9e and 5e. We provide new pollen and sea surface temperatures records for MIS 9e from IODP Site U1385. This pollen record shows a forest expansion in southern Iberia over a 14 ky interval, bracketed by the millennial‐scale cooling events of Termination IV and MIS 9d. Between 334.5 and 332.5 ka, forest expansion reached a maximum, suggesting increased winter moisture during early MIS 9e. Model‐data comparison for MIS 9e and 5e shows that insolation is the main driver of the orbital‐scale vegetation and precipitation changes in Iberia, atmospheric CO2 forcing playing a secondary role. The high‐frequency component of the MIS 9e and 5e forest timeseries highlights the early interglacial forest and precipitation maxima as prominent suborbital events lasting ∼2 ky. We propose that the primarily insolation‐driven forest and precipitation optima were fostered by the non‐equilibrium conditions generated by the millennial‐scale deglacial variability during the early interglacials. Additionally, the early end of these optima may have been favored by a cooling and drying event that is part of the persistent intra‐interglacial variability.

Distinct seasonal changes and precession forcing of surface and subsurface temperatures in the mid-latitudinal North Atlantic during the onset of the Late Pliocene

Abstract. The Late Pliocene marks the intensification of Northern Hemisphere glaciation (iNHG), offering a unique opportunity to study climate evolution and ice-sheet-related feedback mechanisms. In this study, we present high-resolution Mg / Ca-based sea surface temperatures (SSTs) and subsurface temperatures (SubTs) derived from the foraminiferal species Globigerinoides ruber and Globorotalia hirsuta, respectively, at the Integrated Ocean Drilling Program (IODP) Expedition 306 Site U1313 in the mid-latitudinal North Atlantic during the early Late Pliocene, 3.65–3.37 million years ago (Ma). We find distinct differences between our new G. ruber Mg / Ca-based SST record and previously published alkenone-based SST records from the same location. These discrepancies in both absolute values and variations highlight distinctly different seasonal influences on the proxies. The G. ruber Mg / Ca-based SST data were primarily influenced by local summer insolation, showing a dominant precession cycle. Conversely, the variations in alkenone-based SST, dominated by the obliquity and lacking the precession cycle, are found to be more indicative of cold-season changes, despite previous interpretations of these records as reflecting annual mean temperatures. A simultaneous decline in Mg / Ca-based SST and SubT records from 3.65 to 3.5 Ma suggests a diminished poleward oceanic heat transport, implying a weakening of the North Atlantic Current (NAC). A comparison with Early Pleistocene G. ruber Mg / Ca-based SST records shows a shift in the dominant climatic cycle from precession to obliquity, alongside a marked increase in amplitude, indicating an enhanced influence of obliquity cycles correlated with the expansion of Northern Hemisphere ice sheets.

Vanadium isotope records of the transformation from carbonated melt to alkali basalt

Volcanic clasts from the International Ocean Discovery Program Site U1431 in the South China Sea (SCS) were analyzed for V isotopic compositions in this study. These volcanic clasts represent mantle-derived carbonated melts and record the transition from carbonated melts to alkali basalts. Based on the formation sequence, these clast samples were divided into early-stage (> 8.3 Ma) and late-stage (< 8.3 Ma) clasts. The early-stage volcanic clasts have δ51V values of −0.76‰ to −0.67‰, which provide an estimate for the V isotopic composition of primary carbonated melts. Their V isotopic compositions are slightly heavier than those of mid-ocean ridge basalts (−0.84‰ ± 0.10‰) and bulk silicate Earth (−0.86‰ to −0.91‰), reflecting the control of low-degree partial melting under conditions of high oxygen fugacity. The late-stage volcanic clasts have δ51V values ranging from −0.62‰ to 0.29‰, which are systematically heavier than those of early-stage volcanic clasts. The carbonated melts rose through and reacted with the lithospheric mantle and were transformed into alkali basalts by dissolving orthopyroxene and precipitating olivine and/or clinopyroxene. Because the reactants and reaction products do not lead to significant V isotope fractionation, the V isotopes show limited variations during transformation from carbonated melts to alkali basalts. The alkali basalts were emplaced into shallow crust by multiple pulses, in which the magma experienced variable degrees of magma differentiation, with δ51V values close to those of the initial alkali basalts or increasing due to fractional crystallization of FeTi oxides. The change in V isotopic compositions from early-stage volcanic clasts to late-stage volcanic clasts was driven by the varying magnitude of influence exerted by different magmatic processes en route from the deep mantle to the shallow crust.

Diagenetic and hydrothermal controls in a modern submarine rear-arc setting: Insights from Izu-Bonin volcanic arc (IODP Site U1437)

International Ocean Discovery Program (IODP) Expedition 350 drilled Site U1437 in a submarine volcano-bounded basin situated in the modern Izu-Bonin rear-arc, NW Pacific Ocean. Subaqueous volcaniclastic sediments and rocks with a maximum deposition age of 15.4 ± 0.8 Ma were recovered from 0 to 1,806.5 m below seafloor (mbsf). In order to document post-depositional processes in such geological setting, we describe variations in bulk and clay mineralogy over the entire volcaniclastic succession. Four alteration stages (1, 2, 3 and 4) were identified through the occurrence and development of diagenetic background mineral assemblages (smectite ± Na-Ca zeolites ± illite) that were further superimposed by hydrothermal alteration. Stages 1 and 2 are characterized by diagenetic reactions linked with low fluid/rock interactions that enabled glass devitrification and subsequent lithification under burial conditions. Stages 3 and 4 are characterized by moderate to pervasive alteration processes that are well developed in coarser-grained rocks, and that may be induced by thermal pulses associated with fluid inputs. Below 1,460 mbsf, infilling and replacement textures overprinted the background alteration and can be directly linked with the development of two hydrothermal mineral assemblages: (i) ordered C/S (chlorite-smectite mixed-layers) ± chlorite ± albite, and (ii) calcite ± chalcedony ± anhydrite ± laumontite. Both assemblages evidence relatively low-temperature (up to 225 °C) hydrothermal activity that affected subaqueous volcaniclastic rocks at Site U1437. These assemblages are comparable with propylitc alteration facies present in ore-bearing hydrothermal systems. The preferential development of alteration mineral assemblages in high-permeability, coarse-grained lithofacies, reflects the significant influence of the physical properties of volcaniclastic rocks with depth on chemical kinetics, in comparison with those imposed by the local geothermal gradient.

Large variability and 2H-depletion of Middle Miocene to Pleistocene alkenone hydrogen isotopes in the Equatorial Pacific reflect subsurface, low light haptophyte growth

Hydrogen isotope ratios of haptophyte derived long-chain alkenones (δ2HC37:2) have shown to be a useful tool for reconstructing past isotopic compositions of surface seawater (δ2HSSW). The δ2HSSW is related to global ice volume, sea surface salinity and the local hydrological cycle. Here, we present a hydrogen isotope record of alkenones spanning the last 14.5 Ma from IODP site U1338 in the east equatorial Pacific. The alkenone-based reconstructed δ2HSSW is substantially more negative and variable than reconstructed δ2HSSW based on published oxygen isotopes of coccolith carbonates. This suggests that factors other than the isotopic composition of seawater affect the hydrogen isotopic composition of alkenones. The relatively negative and highly variable δ2HC37:2 values are in line with published modern observations on alkenones from suspended particulate matter in the equatorial and north Pacific, with the highest values at relatively high light conditions at the surface and the lowest values at higher water depth and relatively low light conditions. This suggests that the relatively negative and highly variable δ2HC37:2 values in these Middle Miocene to Pleistocene sediments are likely derived from haptophytes growing below the sea surface under variable low light conditions. In regions where the contribution of alkenones from subsurface production, due to high subsurface nutrients, at low light intensities to the sediment is relatively high the δ2HC37:2 has to be interpreted with care as a proxy for δ2HSSW.

Oceanography of the Eastern Equatorial Pacific Ocean Across the Oligocene‐Miocene Transition

AbstractThe functioning of the Pacific Ocean—the world's largest ocean—during a warmer‐than‐present paleoclimate state remains underexplored. We present planktonic and benthic foraminiferal stable oxygen (δ18O) and carbon (δ13C) isotope records from Integrated Ocean Drilling Program (IODP) Site U1334 that span the Oligocene‐Miocene Transition (OMT) interval, from 24.15 to 21.95 million years ago (Ma). We reconstruct (sub‐)surface and deep‐water conditions and provide better constraints on the physical and chemical oceanography of the eastern equatorial Pacific Ocean (EEP). Positive trends in planktonic and benthic foraminiferal δ18O values, mark a largely uniform imprint of increased land‐ice volume/global cooling on surface‐ and deep‐waters. We document a delayed planktonic foraminiferal δ18O increase across the OMT as well as an increase in the amplitude variability of planktonic foraminiferal δ18O values on eccentricity timescales during the early Miocene. We interpret this as an enhanced glacioeustatic sea‐level control on Atlantic‐Pacific salinity exchange through the Central American Seaway (CAS) or as the onset of more variable surface currents and oceanic fronts in the EEP. Positive trends in planktonic and benthic foraminiferal δ13C values characterize the whole‐ocean depletion in 12C linked to organic carbon burial during the Oligocene‐Miocene carbon maximum (CM‐OM). However, this depletion is more pronounced in the planktonic foraminiferal δ13C record, especially during ∼400 Kyr eccentricity minima, reflecting an increase in nutrient upwelling and the efficacy of the biological carbon pump (BCP) when global temperatures decreased across the OMT and during the early Miocene. Our study highlights the dynamic behavior of the EEP in a warmer‐than‐present unipolar icehouse state.

New Insights Into the Surface‐Ocean Dynamics of the Northeastern Atlantic Ocean Across the Marine Isotope Stage 7

AbstractThe upper water‐column dynamics and surface productivity variability in the Northeastern Atlantic Ocean across the Marine Isotope Stage (MIS) 7 interglacial complex is not well understood. Here, we present high‐resolution planktic foraminiferal proxies combined with Artificial Neural Network based‐sea‐surface temperature (SST) and ice‐rafted detritus records from International Ocean Discovery Program Site U1385, SW Iberian Margin for the intervals representing MIS 8 deglaciation, MIS 7 interglacial complex and MIS 6 glacial inception. The long‐term SST pattern is modulated by insolation and precession parameters across the MIS 7 interglacial complex and is superimposed by the millennial‐scale variability (stadials at ∼250, ∼243, ∼230, ∼221, ∼203, ∼196 and ∼192 ka). The regional SST records indicate high temperature gradient (∼6°C) between the sub‐polar North Atlantic Ocean and the SW Iberian Margin during MIS 7d which enhanced the moisture transport from mid‐to‐high latitudes. Further, low obliquity with low insolation induced cooling at high latitudes and promoted the expansion of ice‐sheets during MIS 7d. Comparison of our faunal proxies with the published marine and terrestrial records from SW Europe and western Mediterranean Sea suggested a weakening and southward shift of Azores High (AH) pressure system, similar to the present‐day (−) NAO‐like atmospheric configuration during the early phases of MIS 7e, MIS 7c and MIS 7a, resulting reduced surface productivity of SW Iberian Margin. During the late phases of MIS 7e, MIS 7c and MIS 7a, the AH pressure system strengthened and shifted northward causing the intensification of the trade winds, a scenario similar to present‐day (+) NAO‐like atmospheric configuration, which resulted in high surface productivity of SW Iberian Margin.

Paleomagnetic secular variation and revised chronostratigraphy of Bering Sea (IODP expedition 323) deep-sea sediments (MIS 5)

This study presents new high-resolution, full-vector paleomagnetic secular variation (PSV) records from IODP Expedition 323 Sites U1339, U1343, U1344, and U1345 in the Bering Sea (51°N-60°N) during Marine Isotope Stage 5 (71–130 ka). The chronology of the records has been determined by oxygen isotope stratigraphy with an age uncertainty of ∼±2000 years. The data come from shipboard paleomagnetic measurements of deep-sea sediments with accumulation rates of 26–54 cm/ky. The sample interval in the composite records is 5 cm with a time resolution of ±100–200 years. The composite PSV records for each site were developed by correlation and analysis of 3–4 paleomagnetic records from individual holes at each site. The composite records are a combination of the best correlatable PSV data from each site. We have been able to significantly correlate the PSV records from all four sites. We have identified 56 inclination features (highs/lows) and 54 declination features (east/west extremes) in all the records. Relative paleointensity has been determined by normalizing the sediment natural remanence, after 20-mT af demagnetization, to magnetic susceptibility. This normalization process is complicated by environmental variability, thus our relative paleointensity records are limited to large-scale correlations among the four sites. We correlated our relative paleointensity records to the global PISO-1500 paleointensity record that is itself dated by oxygen isotope stratigraphy. We can identify 6 key intensity highs/lows in our composite records and the PISO-1500 record that provide additional dating isochrons for our records. There are three short magnetic field excursions recorded with replication at these four sites – The youngest excursion occurred at 100 ± 2 ka and is associated with the Fram Strait Excursion. Two older excursions occurred at 116 ± 2 ka and 119.5 ± 2 ka. We associate both of these with the Blake Event. All three of the excursions have very short time intervals (∼200–400 yrs) with somewhat out-of-phase inclination and declination variability. All three of these are Class I excursions. We have studied the statistical characteristics our composite records by calculating the 3ky and 9ky averaged inclinations, declinations, unit vectors, and relative paleointensity over 71–130 ka after removing all excursional directions. The averaged inclinations and declinations are consistent among the four sites. The averaged angular dispersion is bimodal in its overall pattern. Most of the time angular dispersion is low (∼±10°–12°) except for three shorter intervals that have significantly higher angular dispersion (∼±24–30°) and lower paleointensities. The high angular dispersion intervals are associated with the three short excursional intervals. The timing and pattern of interrelated high angular dispersion, low paleointensity, and excursions is synchronous with the PSV pattern recorded during MIS 5 in the central North Atlantic Ocean. We think this defines a distinctive low-intensity/low-energy state for the geomagnetic field.

Origins and paleoclimatic and paleoceanographic significance of laminated sediments of middle Pleistocene age from the southern Bering Sea

Laminated diatomaceous sediments occur intermittently in the Bering Sea over the past 5 million years. A varve (annually deposited) origin for the laminae has been suggested, but there is currently no consensus. Here, we report results of a study of two laminated intervals dating from ∼528 ka and ∼ 782 ka from IODP Site U1340 on the Bowers Ridge. We combine conventional micropaleontological methods with scanning electron microscope analysis that resolves the seasonal cycle of flux events recorded and demonstrates an annual origin for the laminae. Resting spores of Chaetoceros spp., and more rarely of Thalassiosira antarctica represent early spring blooms with the latter likely due to increased meltwater input. Diatom laminae containing concentrations of Shionodiscus trifultus, Thalassiosira spp., Actinocyclus curvulatus, and Neodenticula seminae, represent flux from the main spring-summer bloom episodes, with the latter indicating influx of seed populations from the North Pacific Alaskan Stream. Concentrations of Coscinodiscus spp. in the uppermost part of the diatom laminae represent “Fall Dump” sedimentation of these stratified-adapted diatoms in response to break down of summer stratification in autumn/early-winter storms. The lithogenic laminae represent mainly winter deposition and rare earth element analysis suggests provenance from the southern Bering Sea shelf and the Aleutian Arc. Productivity was high in the studied intervals with total mass fluxes around 5 times higher than modern values. Variation in lamina thickness and diatom composition contain periodicities of 2–8 years, as well as a bi-decadal variability likely related to influence of the Pacific Decadal Oscillation.

Revisiting the relationship between the pore water carbon isotope gradient and bottom water oxygen concentrations

Reconstructing the oxygen content at the ocean floor provides insight into ocean circulation, ventilation, and carbon storage in the deep sea. The microbial breakdown of organic carbon within marine sediment through aerobic respiration consumes oxygen in the pore fluid and releases dissolved inorganic carbon. The offset in the carbon isotopic composition, δ13C, of epifaunal and infaunal foraminifera, Δδ13C, is considered to reflect the aerobic respiration of organic carbon and can be used to reconstruct the oxygen content of the bottom water. Previous work provided an empirical calibration that was suggested to be valid for oxygen reconstructions between 55 and 235 μmol kg−1. In this study, we apply a biogeochemical reactive transport model (RTM) to extend and update this calibration, allowing for the reconstruction of oxygen concentrations ([O2]) up to ∼300 μmol kg−1. Using the RTM and new bottom water [O2] and pore fluid measurements from Iberian Margin sediment cores, we also demonstrate that the calibration between the Δδ13C and bottom water [O2] must account for the coupled changes in the carbon system due to the respiration of organic carbon in the overlying ocean including the concentration and carbon isotopic composition of the dissolved inorganic carbon, and the δ13C of the organic carbon within the sediment column. We apply the improved calibration to reconstruct the changes in oxygen content at International Ocean Discovery Program (IODP) Site U1385 in the deep North Atlantic over the past 1.5 Myr.

Surface water mass dynamics at IODP Site U1313 through Principal Component Analysis: Evidence from coccolith assemblages in the ∼25–7 kyr interval

Coccolith assemblages from Integrated Ocean Drilling Program Site U1313 (North Atlantic; 41°0′N, 32°57.4′W) are analysed to assess temporal changes in surface water conditions from the Last Glacial Maximum (LGM) to the Middle Holocene (∼25–7 kyr). We perform Principal Component Analysis (PCA), involving both calcareous nannofossil assemblage data and other paleoenvironmental proxies. We reconstruct sea-surface dynamics using calcareous nannofossil assemblage variations and comparison with other previously published proxy data from Site U1313 including C37:4 alkenone, dolomite/calcite and quartz/calcite ratios, Uk’37 SST, as well as the sedimentation rate.We use variations in paleoproductivity proxies, cold and warm water group abundance, and the PCA results to infer North Atlantic Current (NAC) and North Atlantic Transitional Waters (NATW) migration through the studied interval. These proxies are also used to highlight variations in the Subtropical Gyre dynamics. Warming episodes are recorded during the LGM and Heinrich Stadial (HS) 1 as a result of northward expansion of the Subtropical Gyre. High coccolithophore productivity episodes characterize HS 1 and 2, and are concomitant with increasing terrigenous input. Both processes are considered a consequence of iceberg influence, which transport nutrients and ice-rafted debris to the area. A southward shift of the Subpolar Front and NAC favour the expansion of cold water from 25.3 ka to the end of the Younger Dryas, whereas the early Holocene is characterized by NATW influence resulting from northward migration of the NAC and Azores Current.

Revisiting the South China Sea MORBs: Mg isotope and whole-rock geochemical constraints

Mid-Ocean Ridge Basalts (MORBs) from the South China Sea (SCS) have attracted great attention from geologists in the past decade, but their petrogenesis still remains controversial. Here, we show that the Mg isotopes of MORBs from Integrated Ocean Drilling Program (IODP) Sites U1431, U1433 and U1434 are heterogeneous and present the mechanisms causing such heterogeneity. The MORBs from Sites U1431, U1433 and U1434 have δ26Mg values ranging from −0.35‰ to −0.21‰, from −0.34‰ to −0.13‰, and from −0.30‰ to −0.20‰, respectively. Together with all the published data, we found the MORBs at Site U1431 with high MgO contents (> 9 wt%) usually have δ26Mg values that are generally lower than the normal mantle (−0.25‰ ± 0.04‰), whereas some of the Site U1433 MORBs are higher than the normal mantle. Shallow-level geological processes like seawater alteration have negligible effects on the Mg isotope fractionation of our samples. Elemental ratio indicative of the degree of partial melting (like Sm/Nd) from the same site also shows no relationship with the δ26Mg variations. Olivine is proposed to have light Mg isotopes. Considering the high MgO (> 9 wt%) contents of low-δ26Mg samples and the petrographic evidence, we propose that their slightly light Mg isotopic compositions of U1431 MORBs may result from the enrichment of olivine, rather than reflecting the mantle source. Moreover, low Ce/Pb ratios and most variable δ26Mg values of some U1433 MORBs suggest that lower continental crustal materials have contributed to their source. As all the previously published MORBs are evolved, we use the Revpet to estimate the mantle melting conditions of SCS MORBs, including the Site U1500 samples recovered from the northern margin of SCS. Compared with the global MORBs data, there is no thermal anomaly in SCS MORBs. Together with our new Mg isotope data, therefore, we propose that the role of the Hainan mantle plume in the opening and spreading of the SCS may be insignificant.

British–Irish Ice Sheet and polar front history of the Goban Spur, offshore southwest Ireland over the last 250 000 years

Deep Sea Drilling Program (DSDP) Site 548 was cored in 1984 at a water depth of 1256 m on the Goban Spur, offshore southwest Ireland. Coring retrieved a ~100‐m‐thick Pleistocene contourite sequence. This study uses planktonic foraminiferal assemblage and benthic foraminiferal oxygen isotope analyses to establish an age model for the upper 40 m of this core. This site's multidisciplinary analyses of planktonic foraminiferal assemblages, lithic grains, facies and calcium carbonate concentration reveal a 250 000‐year record of the North Atlantic polar front variability and British–Irish Ice Sheet (BIIS) history. The sequence is characterized by alternations of ice rafted debris (IRD) laden pelagic mud facies with calcium carbonate‐rich silty sand contourite facies that track glacial/interglacial cycles. The polar front migrated southward across the area several times during glacial maxima and stadial periods, while warmer Mediterranean Outflow Water (MOW) flowed northward across the region during interglacial and interstadial periods depositing contourites. Lithic analyses reveal a complex history of IRD deposition associated with iceberg calving from the Laurentide Ice Sheet and northwest European ice sheets, mainly the BIIS. Comparison between the Goban Spur (DSDP Site 548) and the Celtic Margin (MD03‐2692) and central North Atlantic Integrated Ocean Drilling Program (IODP) Site U1308 suggests differences between the ‘non‐Laurentide Ice Sheet’ Heinrich Events (HE) 6 and 3 at the Goban Spur, with IRD from the BIIS being prominent during HE 6 and IRD from other European ice sheets north of the BIIS likely being more dominant during HE 3. The nature of lithics in IRD‐rich horizons during Terminations 3, 3A, 2 and 1 suggests significant iceberg calving episodes preceding BIIS retreat during the onset of interstadial intervals.

Rapid strengthening of westerlies accompanied intensification of Northern Hemisphere glaciation

The trigger, pace, and nature of the intensification of Northern Hemisphere Glaciation (iNHG) are uncertain, but can be probed by study of ODP Site 1208 North Pacific marine sediments. Herein, we present magnetic proxy data that indicate a 4-fold increase of dust between ~ 2.73 and ~ 2.72 Ma, with subsequent increases at the start of glacials thereafter, indicating a strengthening of the mid-latitude westerlies. Moreover, a permanent shift in dust composition after 2.72 Ma is observed, consistent with drier conditions in the source region and/or the incorporation of material which could not have been transported via the weaker Pliocene winds. The sudden increase in our dust proxy data, a coeval rapid rise in dust recorded by proxy dust data in the North Atlantic (Site U1313), and the Site 1208 shift in dust composition, suggest that the iNHG represents a permanent crossing of a climate threshold toward global cooling and ice sheet growth, ultimately driven by lower atmospheric CO2.

Benthic foraminiferal survival through the early Paleocene (Danian) greenhouse climate interval based on analysis of IODP Site U1457 (Laxmi Basin, Northern Indian Ocean)

The early Paleocene (Danian) period included three major events: the Cretaceous-Paleogene (K-Pg) extinction at 66.0 ​Ma, the Dan-C2 hyperthermal at 65.2 ​Ma, and the latest Danian hyperthermal at ∼62.2 ​Ma. In this paper, we investigate benthic foraminiferal diversity patterns, morphotypes, and oxygen conditions along with the carbonates and magnetic susceptibility records at IODP Site U1457 (Laxmi Basin, Northern Indian Ocean) to understand the effects of these Danian events on the marine community in the Indian Ocean. Findings suggest that foraminifera persisted across these major events. Species belonging to Bolivina, Glandulina, Hoeglundina, Parrelloides and Quadrimorphina genus were dominant above the K-Pg boundary whereas Bolivina, Bulimina, Cassidulina, Cornuspira, Gyroidinoides, Melonis, Oolina, Pullenia, Reussoolina and Rutherfordoides dominated across subsequent hyperthermal events. We calculated the average oxygen content at 0.16 ​ml/L in accordance with oxyphilic species abundance, which shows that the Laxmi Basin comprised, mostly, a suboxic to dysoxic environment. The benthic foraminiferal diversity patterns, primary anomalies of calcium carbonates, and magnetic susceptibility, integrated with previously-published global datasets of carbon and oxygen isotopes, help to define the major geologic events at the study site and show how biotas responded to global change during the early Paleocene greenhouse mode.

Micropaleontological and Isotopic Perspective of Surface Water Mass Variability in the NE Atlantic from MIS 6 to 5e (188 to 115 ka)

The penultimate glaciation (marine isotope stage (MIS) 6) is considered regionally extreme compared to the last glacial maximum, in which the European ice sheets had a vast areal extent. In contrast to the last deglaciation (19–7 ka), the penultimate deglaciation (140–130 ka) hosts one of the most rapid oceanographic changes of the late Pleistocene. In this study, we reconstructed changes in the near-surface and thermocline in the central to northeast Atlantic by analyzing sediments from two Integrated Ocean Drilling Program Expedition 306 sites. Sites U1313 (41°00.6′ N, 32°57.4′ W) and U1314 (56°21.9′ N, 27°53.3′ W) were drilled on the eastern flank of the mid-Atlantic ridge and Gardar Drift of the eastern subpolar North Atlantic, respectively. We analyzed planktonic foraminiferal assemblages, ice-rafted debris (IRD), and oxygen isotopes in two planktonic foraminifers, Globigerina bulloides, and Globorotalia inflata, from MIS 6 to 5e (185–115 ka). Warmer and colder sea-surface conditions were marked by a change in the relative abundance of polar, subpolar, and transitional planktonic foraminifers. Oxygen isotopes in G. bulloides and G. inflata suggest that the thermocline deepened at the subtropical Site U1313 during MIS 6. The lack of Globorotalia inflata prevented us from profiling the mixed layer and thermocline at the subpolar Site U1314. In contrast to MIS 6, the mixed layer and thermocline were re-stratified during the last interglacial. The lack of major IRD events at both sites suggests the stability of the Laurentide ice sheet during MIS 6 compared to the subsequent glaciation. The presence of Heinrich event 11 indicates the discharge of freshwater that freshened the sea surface, resulting in mixing between the mixed layer and thermocline. Our results were placed into a broader context using published data that shed light on the sensitivity of freshwater discharge to the North Atlantic and the following changes with a transition from a penultimate glacial to an interglacial period in surface circulation.

Heavy magnesium isotopic compositions of basalts erupted during arc inception: Implications for the mantle source underlying the nascent Izu-Bonin-Mariana arc

Basalts formed during the early development of an arc are usually buried, but they record critical information relating to the mantle source of the overriding plate prior to the addition of subducting slab components. Basalts recovered at Site U1438 of International Ocean Discovery Program (IODP) Expedition 351, in the Amami Sankaku Basin (ASB), formed during the transition from forearc basalt (FAB) and boninite eruptions to stratovolcano developments in the Izu-Bonin-Mariana (IBM) arc. In this study, we present magnesium (Mg) isotopic data (δ26Mg) for the ASB basalt core samples, formed at arc inception in the absence of down-going slab components, and report the implications of these data for the characteristics of the mantle sources underlying the proto-IBM arc. The δ26Mg values of the ASB basalts range from −0.21‰ to +0.08‰, with an average value of −0.13 ± 0.07‰. These values are systematically higher than those of mid-ocean ridge basalts (MORBs) (−0.25 ± 0.06‰). No obvious effects of post-eruptive alteration, fractional crystallization, partial melting, or subduction component addition can be identified in the Mg isotopic compositions either of the ASB basalts or their mantle source, given the absence of correlations between the δ26Mg values and the proxies of these processes (e.g., K/Nb, Nb/Zr, Ba/La, Nb/Y, etc.). We conclude that the high δ26Mg values of the ASB basalts are inherited from their mantle source that was enriched in heavy Mg isotopes. Melting of an ultra-depleted, refractory spinel peridotite source containing talc-bearing serpentinite components can explain the origin of the heavy Mg isotopic compositions in the ASB basalts. The presence of water-rich serpentinite components in the mantle provides a new perspective on how a refractory mantle source can undergo partial melting at temperatures and pressures similar to the formation of MORBs without the influence of subducted materials, which has important implications for the initiation of plate subduction. Furthermore, the Mg isotopic contributions of subducted additions and mantle source materials should be carefully discussed in future studies due to the heterogeneity of Mg isotopes in the mantle.

Highly Unradiogenic Pb Isotopes Revealed by Plagioclase‐Hosted Melt Inclusions in E‐MORB From the South China Sea

AbstractMid‐ocean ridge basalts (MORBs) can record heterogeneity of the upper mantle. Early crystallized minerals and their melt inclusions commonly record greater isotopic variability than whole‐rocks and are powerful tools to investigate the mantle source. In this study, we develop a new in situ Pb analytical protocol for Pb‐poor materials (sub‐ppm level) and apply it to the plagioclase‐hosted melt‐inclusions from IODP Expedition 349 Site U1433, South China Sea (SCS). Together with element mapping, in situ major and trace element analyzes and in situ Sr isotopes of plagioclase phenocrysts, we show that unlike the low and uniform whole‐rock 87Sr/86Sr ratios, the plagioclase phenocrysts record extreme Sr isotopic heterogeneity. Strontium isotopic heterogeneity is observed between crystals even in a single thin section. Based on the high an contents of plagioclase phenocrysts and high Mg# values of their melt inclusions, we propose they crystallized early in the magma chamber and were most likely formed in different batches of mantle‐derived melts. The plagioclase‐hosted melt inclusions show highly variable Pb isotopes, exhibiting a typical linear variation in 207Pb/204Pb versus 206Pb/204Pb. Such a mixing trend indicates an end‐member with highly unradiogenic Pb isotopes. The LREE enrichment and negative correlations between La/Sm and Pb isotopes of the melt inclusions, indicate that ancient cratonic lower continent crustal materials, with low time‐integrated U/Pb ratios, are likely present in the sub‐ridge mantle of the southwest Sub‐basin, SCS.

A 1.5-million-year record of orbital and millennial climate variability in the North Atlantic

Abstract. Climate during the last glacial period was marked by abrupt instability on millennial timescales that included large swings of temperature in and around Greenland (Daansgard–Oeschger events) and smaller, more gradual changes in Antarctica (AIM events). Less is known about the existence and nature of similar variability during older glacial periods, especially during the early Pleistocene when glacial cycles were dominantly occurring at 41 kyr intervals compared to the much longer and deeper glaciations of the more recent period. Here, we report a continuous millennially resolved record of stable isotopes of planktic and benthic foraminifera at IODP Site U1385 (the “Shackleton Site”) from the southwestern Iberian margin for the last 1.5 million years, which includes the Middle Pleistocene Transition (MPT). Our results demonstrate that millennial climate variability (MCV) was a persistent feature of glacial climate, both before and after the MPT. Prior to 1.2 Ma in the early Pleistocene, the amplitude of MCV was modulated by the 41 kyr obliquity cycle and increased when axial tilt dropped below 23.5∘ and benthic δ18O exceeded ∼3.8 ‰ (corrected to Uvigerina), indicating a threshold response to orbital forcing. Afterwards, MCV became focused mainly on the transitions into and out of glacial states (i.e. inceptions and terminations) and during times of intermediate ice volume. After 1.2 Ma, obliquity continued to play a role in modulating the amplitude of MCV, especially during times of glacial inceptions, which are always associated with declining obliquity. A non-linear role for obliquity is also indicated by the appearance of multiples (82, 123 kyr) and combination tones (28 kyr) of the 41 kyr cycle. Near the end of the MPT (∼0.65 Ma), obliquity modulation of MCV amplitude wanes as quasi-periodic 100 kyr and precession power increase, coinciding with the growth of oversized ice sheets on North America and the appearance of Heinrich layers in North Atlantic sediments. Whereas the planktic δ18O of Site U1385 shows a strong resemblance to Greenland temperature and atmospheric methane (i.e. Northern Hemisphere climate), millennial changes in benthic δ18O closely follow the temperature history of Antarctica for the past 800 kyr. The phasing of millennial planktic and benthic δ18O variation is similar to that observed for MIS 3 throughout much of the record, which has been suggested to mimic the signature of the bipolar seesaw – i.e. an interhemispheric asymmetry between the timing of cooling in Antarctica and warming in Greenland. The Iberian margin isotopic record suggests that bipolar asymmetry was a robust feature of interhemispheric glacial climate variations for at least the past 1.5 Ma despite changing glacial boundary conditions. A strong correlation exists between millennial increases in planktic δ18O (cooling) and decreases in benthic δ13C, indicating that millennial variations in North Atlantic surface temperature are mirrored by changes in deep-water circulation and remineralization of carbon in the abyssal ocean. We find strong evidence that climate variability on millennial and orbital scales is coupled across different timescales and interacts in both directions, which may be important for linking internal climate dynamics and external astronomical forcing.

East Asian summer monsoon evolution recorded by the middle Miocene pelagic reddish clay, South China Sea

The East Asian summer monsoon (EASM) is essential for the livelihoods of billions of people in East and Southeast Asia. Its present features and evolution history of the Quaternary have been well studied, yet how did EASM operate in much older geological times, especially before the late Miocene, is relatively less constrained. Scientific drillings conducted in the South China Sea (SCS) offer excellent opportunities to resolve this issue to a certain extent considering the EASM has been prevailing in this area since at least the earliest Miocene. In this study, mineralogical and stratigraphic studies were carried out on the middle Miocene pelagic reddish clay recovered at International Ocean Discovery Program Site U1433 using clay mineralogical and natural gamma radiation data. Our results show that the smectite content and illite crystallinity of the deep SCS clay were direct consequences of the EASM intensity varying in precession rhythm and modulated by eccentricity cycles. A high-resolution temporal framework constructed with onboard age controls and cyclostratigraphic analyses demonstrates that the pelagic reddish clay from Site U1433 was accumulated from 17.17 to 12.69 Ma. During this period, deep-sea sedimentation of the SCS was controlled by the EASM and influenced by the sea-level change. Importantly, clay mineral is proved to be an efficient terrigenous standard of astronomical tuning and stratigraphic correlation for pelagic sediments. It is also verified that the EASM intensified during the Miocene Climatic Optimum and the following Middle Miocene Climate Transition, and induced higher humidity and seasonality, respectively. Then, the EASM gradually weakened since the expansion of the East Antarctic Ice Sheet.