Seafloor Spreading Rates Research Articles

Unexpectedly High Magma Productivity Inferred From Crustal Roughness and Residual Bathymetry on the Eastern Part of the Ultra‐Slow Spreading Gakkel Ridge Since ∼45 Ma, Eurasian Basin, Arctic Ocean

AbstractThe Gakkel Ridge in the Eurasian Basin has the slowest seafloor spreading worldwide. The western Gakkel Ridge (3°W–85°E; 14–11 mm/a) alternate between magmatic and sparsely magmatic zones, while the eastern Gakkel Ridge (85–126°E; 11–6 mm/a) appears to be dominated by magmatic zones despite ultraslow spreading. Little is known about the seafloor spreading conditions in the past along the entire ridge. Here, we exploit the residual bathymetry and basement roughness to assess the crustal accretion process of the Gakkel Ridge over time using 23 published regional multichannel seismic reflection profiles. Full seafloor spreading rates were faster (20–24 mm/a) up to ∼45 Ma, and residual bathymetry for the older crust is deeper than the world average in the entire Eurasian Basin. There is a sharp transition to 300–400 m shallower residual bathymetry for seafloor <45 Ma in the eastern Eurasian Basin. The crustal roughness versus spreading rate of the western Eurasian Basin is on the global trend, while that of the eastern is significantly below. Both low roughness and shallow residual bathymetry of the eastern Eurasian Basin is close to that of oceanic crust for spreading rates above 30 mm/a, demonstrating increased magmatic production of the eastern Gakkel Ridge since ∼45 Ma. A recent mantle tomography model predicts partial melting in the upper mantle based on the low Vs anomaly underneath. The sedimentary pattern toward the Lomonosov Ridge indicates that this hot mantle anomaly started to cause dynamic uplift of the area at ∼45 Ma.

Strain localization instabilities and the genesis of multiple axes of seafloor spreading in the Carmen basin, southern Gulf of California

We present new insight into the tectonic evolution of the Carmen basin (CB) in the southern Gulf of California (GC) from high-resolution bathymetry and two-dimensional seismic reflection data. Our goal is to document the seafloor morphology and sub-surface acoustic characteristics to understand the structure and crustal lithology across the CB. We identify three sub-basins with distinct geometries and evolutionary histories, with basement structures displaying a strong affinity with highly reflective, sigmoidal-shaped layers and the emplacement of high-amplitude tabular material underlying sediments with varying stratigraphic thicknesses. From the extent of new oceanic crust accreted along the CB, we estimate the age of the basin using a seafloor spreading rate of 52 mm/year, which is the average velocity of the relative plate motion between the Baja California microplate and the North American plate, as documented by previous authors. The southern and central sub-basins of the CB are mainly abandoned, while the northern one is currently the locus of seafloor spreading. This is evidenced by the juxtaposition of oceanic crust younger than ~ 1.9 Ma against older oceanic crust correlating in age with the adjacent Guaymas and Farallon basins to the northeast and southwest, respectively. We propose that mantle upwelling beneath the CB is a northward continuation of the East Pacific Rise, resulting in a fast-evolving system with sharp variations in strain localization within the seafloor spreading centers of the CB.Graphical

Seafloor spreading and the delivery of sulfur and metals to Earth’s oceans

Abstract Circulating fluids within Earth’s mid-ocean ridge system cool and alter the oceanic crust, contribute distinct chemistry to the ocean, and generate economically and geologically important metal-sulfide deposits at the seafloor. Yet, we have few constraints on the characteristics of these fluids at peak subseafloor pressure and temperature conditions or how the primary variable, seafloor spreading, affects these fluids’ delivery of metals and sulfur to seawater. Here, we develop a new, robust technique for estimating the peak endowment of heat and dissolved sulfur, iron, and copper in subseafloor hydrothermal fluids and determining their fate as these superheated fluids rise to the seafloor. Calculations using this technique indicate that >20%–70% of sulfur, iron, and copper dissolved at peak subseafloor conditions are lost during upflow due to cooling and concomitant decreases in sulfide mineral solubility. The interpretation of these estimates within the geologic context of vent fields allows us to demonstrate a strong inverse relationship between seafloor spreading rate and peak pressure-temperature conditions, subseafloor heat loss, and the magnitude of subseafloor sulfide mineralization. Our results demonstrate the extent to which the secular variation of Earth’s mid-ocean ridge system over geologic time has impacted sulfide deposition rates and hydrothermal fluxes of sulfur and metals to the ocean.

Revised onset age of magnetochron M0r: Chronostratigraphic and geologic implications

Abstract The timing of the onset of magnetochron M0r and its duration are disputed, reflecting both a limited set of radioisotopic dates and uncertain magnetostratigraphic correlations. We present a chronostratigraphic framework for a reversed polarity interval based on two chronometers (40Ar/39Ar, U-Pb) and newly published paleomagnetic data from the Qingshan Group, Jiaolai Basin, China. Bayesian modeling of U-Pb zircon and 40Ar/39Ar sanidine dates suggests a minimum duration of 540 ± 37 k.y. (95% credible interval) for the reversed polarity interval. These findings are compatible with an astrochronologic age model for M-sequence seafloor magnetic anomalies, indicating that the reversely magnetized sediments correspond to magnetochron M0r rather than the shorter chron “M-1r.” Integration of U-Pb and 40Ar/39Ar ages constrains the onset of M0r to 120.29 ± 0.09 Ma, which is ~1 m.y. younger than that inferred in the current geologic time scale (GTS 2020). This finding also implies that the Cretaceous normal superchron (CNS) began at 119.70 ± 0.12 Ma and that the average seafloor spreading rate during the CNS was ~3.5% higher than that inferred from GTS 2020. It also suggests that oceanic anoxic event 1a began at 119.40 ± 0.12 Ma, thereby providing an updated chronologic basis for exploring the primary trigger of this carbon cycle perturbation.

Deep time spatio‐temporal data analysis using pyGPlates with PlateTectonicTools and GPlately

AbstractPyGPlates is an open‐source Python library to visualize and edit plate tectonic reconstructions created using GPlates. The Python API affords a greater level of flexibility than GPlates to interrogate plate reconstructions and integrate with other Python workflows. GPlately was created to accelerate spatio‐temporal data analysis leveraging pyGPlates and PlateTectonicTools within a simplified Python interface. This object‐oriented package enables the reconstruction of data through deep geologic time (points, lines, polygons and rasters), the interrogation of plate kinematic information (plate velocities, rates of subduction and seafloor spreading), the rapid comparison between multiple plate motion models, and the plotting of reconstructed output data on maps. All tools are designed to be parallel‐safe to accelerate spatio‐temporal analysis over multiple CPU processors.

Planetary geodynamics and age constraints on circumstellar habitable zones around main sequence stars

AbstractPlanetary geodynamics may have an important influence over planetary habitability and the boundaries of the circumstellar habitable zone (CHZ) in space and time. To investigate this we use a minimal parameterized model of the co-evolution of the geosphere and atmosphere of Earth-like planets around F, G, K and M main sequence stars. We found the CHZ for the present Solar System located between 0.92 and 1.09 au for a 1.0 M$_{\oplus }$Earth-like planet, extendible to 1.36 au for a 4.0 M$_{\oplus }$planet. In the literature, the CHZ varies considerably in width and border location, but the outer edges tend to be more spread out than the inner edges, showing a higher difficulty in determining the outer edge. Planetary mass has a considerable effect on planetary geodynamics, with low-mass planets cooling down faster and being less capable of maintaining a rich carbon dioxide atmosphere for several billions of years. Age plays a particularly important role in the width of the CHZ as the CHZ contracts in both directions: from the inner edge (as stellar luminosity increases with time), and from the outer edge (as planetary heat flux and seafloor spreading rate decrease with time). This strongly affects long-lived habitability as the 5 Gyr continuous CHZ may be very narrow or even non-existent for low-mass planets (<0.5 M$_{\oplus }$) and fast-evolving high-mass stars (>1.1 M$_{\odot }$). Because of this, the mean age of habitable terrestrial planets in our Galaxy today may be younger than Earth's age. Our results suggest that the best targets for future surveys of biosphere signatures may be planets between 0.5 and 4.0 M$_{\oplus }$, in systems younger than the Solar System. These planets may present the widest and long-lived CHZ.

Diachronous basin evolution along northern South China Sea: Result of a migrating Hainan plume?

The diachronous basin evolution from west to east along the northern margin of the South China Sea (SCS) is closely related to the formation of the SCS, but in details remains ambiguous. We focused our study on the crustal structure and Cenozoic subsidence and sedimentary evolution of the Pearl River Mouth Basin (PRMB) and the Qiangdongnan Basin (QDNB) to reveal the divergences between the western and eastern segments of this margin. The back-stripping method and gravity modeling were employed based on drilling wells and geophysical data. We revealed that the crustal structure of the PRMB has a wedge shape, thinning monotonously from the coastline to the continent ocean transition (COT), while that of the QDNB has a symmetrical shape, thinning intensively in the Central Canyon. The characteristics of multi-episodic sedimentation and subsidence evolution were revealed in the PRMB, and its rapid subsidence rates were roughly synchronized with the decrease in the seafloor spreading rate of the SCS basin. However, the sedimentary and subsidence rate curves mainly featured accelerated sedimentation and subsidence in the western QDNB and an accelerated subsidence rate curve but a slow sedimentation rate curve in the eastern QDNB after 11.6 Ma. Some low-velocity conduits, which might be the crustal magmatic footprints of the Hainan mantle plume, were revealed in the middle and eastern QDNB. One possible dynamic model that can cause accelerated subsidence after 11.6 Ma in the eastern QDNB is that the Hainan mantle plume has migrated northwestward which may be derived from the contemporaneous southeastward movement of the South China Block. This conjecture coincides with magmatism and subsidence rate variations in the northwestern margin and is evidenced by the low-velocity conduits in the eastern QDNB.

New advances in the stratigraphy of Aptian oceanic anoxic events (Castro Urdiales, Basque-Cantabrian Basin, Spain)

The Aptian sedimentary deposits of the Castro Urdiales area in north Spain record an episode of disoxia/anoxia which postdates Oceanic Anoxic Event 1a (OAE 1a). Carbonate platform and overlying facies are analysed for stratigraphy, sedimentology, total organic carbon and carbon isotopes. Lower Aptian shallow-water limestones are covered by deeper-water marlstones and organic-carbon-rich black lutites of lower to upper Aptian, ranging from Dufrenoyia furcata Zone to Cheloniceras (Epicheloniceras) martinioides Zone. Biostratigraphy and chemostratigraphic calibration are based on ammonites, carbon isotopes and total organic carbon analyses, revealing consistency with the global Aptian reference framework. The Castro Urdiales black-shale levels of the lower to upper Aptian transition together with their encasing series are correlated with other Aptian sections of the Basque-Cantabrian Basin (northern Spain) and southeastern France. The Castro Urdiales black shales are correlated with the Aparein Level of the eastern Basque-Cantabrian Basin. Prior to this oxygen deficient episode, glauconite-and-ostreid facies were also deposited under oxic conditions during the Gutiolo volcanic event of the Basque-Cantabrian Basin. This volcanism correlates in part with the Cretaceous superplume eruptions of Ontong Java, Manihiki Plateaus and the North Atlantic opening as a result of increasing seafloor spreading rates.

Detailed reconstructions of India-Somalia plate motion, 60 Ma to present: Implications for Somalia plate absolute motion and India-Eurasia plate motion

Summary We estimate India-Somalia plate motion at 45 times since 60 Ma from ∼9,000 crossings of Carlsberg and northern Central Indian Ridge magnetic reversals C1 to C26 and numerous fracture zone crossings. The new rotations reveal at least seven significant spreading rate changes since ∼60 Ma, some previously unknown. The largest changes occurred before 46 Ma, when the forces acting on the India plate evolved rapidly due to the transition from subduction to continent-continent collision between India and Asia and the influence of the Reunion hotspot plume on the India plate. The new rotations reveal a gradual ∼50 percent decline in Carlsberg Ridge spreading rates from 57-52.7 Ma, but an end to the decline at ∼53 Ma, when spreading rates surged rapidly by up to 100 percent. From 52-46.7 Ma, Carlsberg Ridge spreading rates collapsed by ∼90 percent, possibly defining a protracted transition to continental collision between India and Asia. Significant kinematic events since 46.7 Ma have included a ∼25 percent spreading rate recovery from 42-40 Ma, ultraslow spreading from 38.6-33.2 Ma, a gradual rate doubling from 33-18 Ma, a ∼50 percent slowdown from 18-13 Ma, and apparently steady motion since 13 Ma. The new rotations successfully predict Carlsberg Ridge abyssal hill orientations for seafloor ages of 48-42 Ma and 20-0 Ma and Central Indian Ridge fracture zone traces for seafloor ages of 43 to 16 Ma, constituting useful tests of the rotation accuracies at these ages. When corrected for the movement of India relative to the Capricorn plate since 16 Ma, the new rotations also successfully restore magnetic lineations C13, C18, C20, and C21, and fracture zone segments from the Capricorn plate onto of their Somalia plate counterparts. This further confirms the accuracies of our new rotations back to C21n.1o (47.8 Ma), and validates a ∼16 Ma start date for India-Capricorn plate motion and published correction for India-Capricorn motion. Anticorrelated changes in India-Somalia and Antarctic-Somalia seafloor spreading rates from 37-18 Ma may be evidence that Somalia plate absolute motion changed during this period, possibly triggered by Somalia’s post-30-Ma detachment from the Arabian Peninsula or the kinematic effects of the Afar and/or Reunion mantle plumes on the Somalia plate. New India-Eurasia rotations that we estimate from an updated global plate circuit predict convergence rates from 53-47 Ma that are ∼30 percent faster than previous estimates and that decline ∼75 percent by ∼38 Ma. Changes in India-Somalia and India-Eurasia rates correspond closely with recently described Tibetan deformation pulses, consistent with linkages between all three. A joint inversion of Carlsberg and southern Central Indian ridge magnetic reversal and fracture zone data, including a correction for movement of the Capricorn plate relative to India, satisfactorily realigns the reconstructed magnetic lineations and fracture zones back to C23n.1n (50.7 Ma), but misfits some data by 100 km or more at earlier times. The misfits may be evidence for deformation within the IndoCapricorn and/or Somalia plates before 48 Ma or a misinterpretation of magnetic reversal and/or fracture zone data from times before 48 Ma.

A pulse of the Earth: A 27.5-Myr underlying cycle in coordinated geological events over the last 260 Myr

We performed spectral analyses on the ages of 89 well-dated major geological events of the last 260 Myr from the recent geologic literature. These events include times of marine and non-marine extinctions, major ocean-anoxic events, continental flood-basalt eruptions, sea-level fluctuations, global pulses of intraplate magmatism, and times of changes in seafloor-spreading rates and plate reorganizations. The aggregate of all 89 events shows ten clusters in the last 260 Myr, spaced at an average interval of ~ 26.9 Myr, and Fourier analysis of the data yields a spectral peak at 27.5 Myr at the ≥ 96% confidence level. A shorter period of ~ 8.9 Myr may also be significant in modulating the timing of geologic events. Our results suggest that global geologic events are generally correlated, and seem to come in pulses with an underlying ~ 27.5-Myr cycle. These cyclic pulses of tectonics and climate change may be the result of geophysical processes related to the dynamics of plate tectonics and mantle plumes, or might alternatively be paced by astronomical cycles associated with the Earth’s motions in the Solar System and the Galaxy.

Emplacement processes of proto‐arc basalt in the Izu–Bonin–Mariana arc system

AbstractAscertaining the emplacement mechanism of oceanic basaltic lavas is important in understanding how ocean floor topography is produced and oceanic plates evolve, particularly during the early stages of crustal development of a supra‐subduction zone. A detailed study of the volcanic stratigraphy at International Ocean Discovery Program (IODP) Site U1438 in the Amami Sankaku Basin, west of the Kyushu–Palau Ridge, has revealed the development of lava accretion and ridge topography on the Philippine Sea plate at about 49 Ma. Igneous basement rocks penetrated at Site U1438 are the uppermost 150 m of ~6 km‐thick oceanic crust, and comprise, in a downhole direction, sheet flows (12.6 m), lobate sheet flows (61.3 m), pillow lavas (50.7 m), and thin sheet flows (25.3 m). The lowermost sheet flows are intercalated with layers of limestone and epiclastic tuff. Lithofacies analysis reveals that the lowermost sheet flows, limestone, and tuff formed on an axial rise, the pillow lavas were emplaced on a ridge slope, and the lobate sheet flows formed off ridge on an abyssal plain. The lithofacies of the basement basalt corresponds to the upper portions of fast‐spreading oceanic crust, suggesting that subduction initiation was associated with intermediate to fast rates of seafloor spreading. The surface sheet flows are olivine–clinopyroxene‐phyric basalt and differ from the lower basalt flows that contain phenocrysts of olivine and plagioclase, with or without clinopyroxene. The depleted chrome‐spinel composition and olivine–clinopyroxene phenocryst assemblage in the surface sheet flows suggests a slight contribution of water for magma generation not present for the lower basalt flows. Considering the lithological difference between the backarc and forearc oceanic crust in the Izu–Bonin–Mariana arc, with sheet flow dominant in the former, seafloor spreading occurred faster in the later stage of subduction initiation.

A 9 million-year-long astrochronological record of the early–middle Eocene corroborated by seafloor spreading rates

The early-middle Eocene (ca. 56-41 Ma) is recorded in the pelagic Scaglia Rossa and Variegata Formations of the Umbria-Marche Basin (central Italy). Geochemical and magnetostratigraphic alignment between the Bottaccione section (Gubbio, central Italy) and the Smirra core (Cagli, central Italy) allows us to generate a continuous and wellpreserved new record that, combined with previously published data from the same area, creates a continuous high-resolution record from the Paleocene-Eocene thermal maximum (ca. 56 Ma) to the lower part of chron C21n. Comparison with carbon isotope records from Ocean Drilling Program Sites 1258 and 1263 reveals a satisfactory match, providing further evidence of the global significance of the long-term trend and superposed perturbations captured by the δ13C records. The identification of astronomically forced geochemical cycles allows us to develop a 405 k.y. tuned age model, thereby extending the astrochronology from ca. 56.0 to ca. 47.5 Ma. Marine magnetic anomaly profiles from major oceanic basins characterized by high seafloor spreading rates were used to independently test the astronomical polarity time scale associated with our tuning, as well as other polarity time scales. Our age model suggests the existence of periods of relatively constant seafloor spreading rates separated by rapid changes, while the other time scales generate more gradual variations and also include large and short-term deviations in spreading rates that occur simultaneously in different oceanic basins, implying errors in polarity reversal ages. The Umbria- Marche age model further contributes to the closure of the middle Eocene gap in the astronomical time scale.

Oxygen isotope composition of the Phanerozoic ocean and a possible solution to the dolomite problem

The 18O/16O of calcite fossils increased by ∼8‰ between the Cambrian and present. It has long been controversial whether this change reflects evolution in the δ18O of seawater, or a decrease in ocean temperatures, or greater extents of diagenesis of older strata. Here, we present measurements of the oxygen and ‟clumped" isotope compositions of Phanerozoic dolomites and compare these data with published oxygen isotope studies of carbonate rocks. We show that the δ18O values of dolomites and calcite fossils of similar age overlap one another, suggesting they are controlled by similar processes. Clumped isotope measurements of Cambrian to Pleistocene dolomites imply crystallization temperatures of 15-158 °C and parent waters having δ18OVSMOW values from -2 to +12‰. These data are consistent with dolomitization through sediment/rock reaction with seawater and diagenetically modified seawater, over timescales of 100 My, and suggest that, like dolomite, temporal variations of the calcite fossil δ18O record are largely driven by diagenetic alteration. We find no evidence that Phanerozoic seawater was significantly lower in δ18O than preglacial Cenozoic seawater. Thus, the fluxes of oxygen-isotope exchange associated with weathering and hydrothermal alteration reactions have remained stable throughout the Phanerozoic, despite major tectonic, climatic and biologic perturbations. This stability implies that a long-term feedback exists between the global rates of seafloor spreading and weathering. We note that massive dolomites have crystallized in pre-Cenozoic units at temperatures >40 °C. Since Cenozoic platforms generally have not reached such conditions, their thermal immaturity could explain their paucity of dolomites.

The crustal structure of the Enderby Basin, East Antarctica

The passive margin and ocean crust of the Enderby Basin, East Antarctica preserves a record of the breakup of East Gondwana. Using a suite of public domain geophysical data, we have examined and described the crustal morphology of the basin. Based on our geophysical observations, we divide the Enderby Basin into three distinct morphologic domains. The Eastern Domain demonstrates the most volcanic morphology of the basin, with abundant seaward dipping reflector packages and anomalously thick oceanic crust. These features suggest an early influence by the Kerguelen Hotspot on continental breakup within the domain. The Central Domain is characterized by two regions of oceanic crust of varying morphology segregated by a high amplitude magnetic anomaly. Geophysical observations suggest that the basement directly inboard of this magnetic anomaly is composed of thin, rugged, and poorly structured, proto-oceanic crust, similar in morphology to oceanic crust formed at ultraslow/slow mid-ocean ridged. Outboard of this anomaly, oceanic crust appears to be well-structured and of normal thickness. We offer three, non-exclusive, explanations for the observed change in ocean crustal structure: (1) melt production was initially low at the time of continental breakup, and the progressive decompression of the mantle led to a gradual increase in melt production and ocean crust thickness, (2) melt production was initially low to due lower extension rates and that melt production increased following a change in spreading rate, (3) a change in spreading ridge geometry led to more effective seafloor spreading rate and concurrent increase in melt production. The Western Domain of the Enderby Basin is characterized by abundant fracture zones and anomalously thin oceanic crust. We believe these features arose as a geometric consequence of the originally oblique orientation of continental rifting relative to the extension direction within the domain. Together these observations suggest that the breakup of East Gondwana was highly variable, with notable along-strike differences in crustal deformation and seafloor spreading processes.

The Impacts of Seawater Mg/Ca and Temperature on Element Incorporation in Benthic Foraminiferal Calcite

AbstractOn geological timescales, oceanic [Mg2+] and [Ca2+] vary with changing rates of weathering, seafloor spreading and dolomite formation. Accurate reconstruction of the ratio between [Mg2+] and [Ca2+] in seawater (Mg/Casw), may potentially be reconstructed using foraminiferal Mg/Ca ratios. Since both temperature and seawater Mg/Ca impact foraminiferal Mg/Ca, successful reconstruction of Mg/Casw requires quantification of both these parameters independently on foraminiferal Mg/Ca, as well as their combined effect on Mg‐incorporation. Here we present the combined and isolated impacts of temperature and Mg/Casw on Mg incorporation in two model species, the benthic hyaline (i.e., perforate) foraminifer Elphidium crispum and porcelaneous (i.e., miliolid) foraminifer Quinqueloculina sp. using controlled growth experiments. Specimens of these two species were kept at four different temperatures (ranging from 10 to 27°C) and three Mg/Casw's (3.4, 6.4 and 8.5 mol/mol), resulting in 12 experimental conditions. Newly grown calcite was analyzed for a number of elements (Na, Mg and Sr) by laser ablation‐ICP‐MS. Results show that although the Mg/Ca varied by more than an order of magnitude between species, the sensitivity of Mg incorporation with respect to temperature appeared not to be influenced by Mg/Casw. By extension, these results may also help improving accuracy in the reconstruction of past Mg/Casw based on foraminifera with contrasting Mg/Ca.

Astronomical calibration of the Ypresian timescale: implications for seafloor spreading rates and the chaotic behavior of the solar system?

Abstract. To fully understand the global climate dynamics of the warm early Eocene with its reoccurring hyperthermal events, an accurate high-fidelity age model is required. The Ypresian stage (56–47.8 Ma) covers a key interval within the Eocene as it ranges from the warmest marine temperatures in the early Eocene to the long-term cooling trends in the middle Eocene. Despite the recent development of detailed marine isotope records spanning portions of the Ypresian stage, key records to establish a complete astronomically calibrated age model for the Ypresian are still missing. Here we present new high-resolution X-ray fluorescence (XRF) core scanning iron intensity, bulk stable isotope, calcareous nannofossil, and magnetostratigraphic data generated on core material from ODP Sites 1258 (Leg 207, Demerara Rise), 1262, 1263, 1265, and 1267 (Leg 208, Walvis Ridge) recovered in the equatorial and South Atlantic Ocean. By combining new data with published records, a 405 kyr eccentricity cyclostratigraphic framework was established, revealing a 300–400 kyr long condensed interval for magnetochron C22n in the Leg 208 succession. Because the amplitudes are dominated by eccentricity, the XRF data help to identify the most suitable orbital solution for astronomical tuning of the Ypresian. Our new records fit best with the La2010b numerical solution for eccentricity, which was used as a target curve for compiling the Ypresian astronomical timescale (YATS). The consistent positions of the very long eccentricity minima in the geological data and the La2010b solution suggest that the macroscopic feature displaying the chaotic diffusion of the planetary orbits, the transition from libration to circulation in the combination of angles in the precession motion of the orbits of Earth and Mars, occurred ∼ 52 Ma. This adds to the geological evidence for the chaotic behavior of the solar system. Additionally, the new astrochronology and revised magnetostratigraphy provide robust ages and durations for Chrons C21n to C24n (47–54 Ma), revealing a major change in spreading rates in the interval from 51.0 to 52.5 Ma. This major change in spreading rates is synchronous with a global reorganization of the plate–mantle system and the chaotic diffusion of the planetary orbits. The newly provided YATS also includes new absolute ages for biostratigraphic events, magnetic polarity reversals, and early Eocene hyperthermal events. Our new bio- and magnetostratigraphically calibrated stable isotope compilation may act as a reference for further paleoclimate studies of the Ypresian, which is of special interest because of the outgoing warming and increasingly cooling phase. Finally, our approach of integrating the complex comprehensive data sets unearths some challenges and uncertainties but also validates the high potential of chemostratigraphy, magnetostratigraphy, and biostratigraphy in unprecedented detail being most significant for an accurate chronostratigraphy.

Structural variation of the oceanic Moho in the Pacific plate revealed by active-source seismic data

The characteristics of the oceanic Moho are known to depend on various factors, such as seafloor spreading rate, crustal age, and accretionary processes at a ridge. However, the effect of local magmatic activities on the seismic signature of the Moho is poorly understood. Here an active-source reflection and refraction survey is used to investigate crustal structure and Moho characteristics along a >1000-km-long profile southeast of the Shatsky Rise in a Pacific Ocean basin formed from the Late Jurassic to Early Cretaceous and spanning the onset of Shatsky Rise volcanism. Although the seismic velocity structure estimated from the refraction data showed typical characteristics of the oceanic crust of the old Pacific plate, the appearance of the Moho reflections was spatially variable. We observed clear Moho reflections such as those to be expected where the spreading rate is fast to intermediate only at the southwestern end of the profile, whereas Moho reflections were diffuse, weak, or absent along other parts of the profile. The poor Moho reflections can be explained by the presence of a thick crust–mantle transition layer, which is temporally coincident with the formation of the Shatsky Rise. We inferred that the crust–mantle transition layer was formed by changes in on-axis accretion process or modification of the primary Moho by off-axis magmatism, induced by magmatic activity of the Shatsky Rise.

Temporal record of osmium concentrations and 187Os/188Os in organic-rich mudrocks: Implications for the osmium geochemical cycle and the use of osmium as a paleoceanographic tracer

We present a compilation of 192Os concentrations (representing non-radiogenic Os) and initial 187Os/188Os isotope ratios from organic-rich mudrocks (ORM) to explore the evolution of the Os geochemical cycle during the past three billion years. The initial 187Os/188Os isotope ratio of a Re-Os isochron regression for ORM constrains the local paleo-seawater 187Os/188Os, which is governed by the relative magnitudes of radiogenic Os (old continental crust) and unradiogenic Os (mantle, extraterrestrial, and juvenile/mafic/ultramafic crust) fluxes to seawater. A first-order increase in seawater 187Os/188Os ratios occurs from the Archean to the Phanerozoic, and may reflect a combination of increasing atmosphere-ocean oxygenation and weathering of progressively more radiogenic continental crust due to in-growth of 187Os from radioactive decay of 187Re. Superimposed on this long-term trend are shorter-term fluctuations in seawater 187Os/188Os ratios as a result of climate change, emplacement of large igneous provinces, bolide impacts, tectonic events, changes in seafloor spreading rates, and lithological changes in crustal terranes proximal to sites of ORM deposition. Ediacaran-Phanerozoic ORM have mildly higher 192Os concentrations overall compared with pre-Ediacaran Proterozoic ORM based on the mean and 95% confidence interval of 10,000 median values derived using a bootstrap analysis for each time bin (insufficient Archean data exist for robust statistical comparisons). However, there are two groups with anomalously high 192Os concentrations that are distinguished by their initial 187Os/188Os isotope ratios. Ediacaran-Cambrian ORM from South China have radiogenic initial 187Os/188Os, suggesting their high 192Os concentrations reflect proximal Os-rich crustal source(s), ultraslow sedimentation rates, and/or other unusual depositional conditions. In contrast, the unradiogenic initial 187Os/188Os and high 192Os concentrations of some Mesozoic ORM can be tied to emplacement of large igneous provinces. Excluding these two anomalous groups and repeating the bootstrap analysis, we find that, overall, the 192Os concentrations for the Ediacaran-Phanerozoic and pre-Ediacaran Proterozoic time bins are not significantly different. An improved understanding of Os geochemical behavior in modern environments is required before our compilation can be fully used to constrain the temporal evolution of the seawater Os reservoir.

Statistics of the Earth’s Topography

The so-called “typical” values for the Earth’s topography are often used in the literature, such as the mean continental altitude (MCA), the Moho depth for “normal” continental crust, or “typical” depth of mid-oceanic ridges. However, the statistical relevance of those values is hardly discussed. Focussed on data for the global topography, this paper presents statistical analyses regarding various environments. It is shown in particular that the definition of the mid-oceanic ridge is not straightforward, and varies considerably according to what is actually considered: the ridge “inner-rift”, the ridge “crest”, or the “virtual ridge” at spreading centre. This definition is also a function of the spreading rate and has strong implications for the rationale on the age-depth relationship of the sea-floor. In addition, the latter relationship is highly dependent on how the topographic data are corrected from sediment load. The correction itself implies numerous aspects that relies on the precision and associated uncertainties of, in particular, the sediment thickness, sediment porosity, and the mantle, water and sediment densities. In this respect, the analysis carried out here favours a plate cooling model (PCM) for the age-depth dependence of the sea-floor. The topographic elevation at trench proves also to be related to the age of the sea-floor through a different PCM equation. Away from the trench, the oceanic lithosphere is affected by flexuration, for which equations can be defined assuming that the end-load position is not located at trench. On the other hand, the elevation of magmatic arc does not appear to be related to sea-floor age or spreading rate. However, the correlation between the arc-trench distance and the topographic elevation of arc for continental crust seems to be an indicator of slab dip and therefore the existence of slab roll-back processes. Along intra-oceanic magmatic arc, a periodicity in topographic elevation suggests a periodicity in the occurrence of magma chambers, and therefore magmatic processes that need to be further studied. At passive margin, the transition between continental and oceanic crust seems to be relatively sharp in average. Subdivision of the datasets according to the age of the continent-ocean boundary (COB) indicates that rift and passive margin shoulders are found within a couple of degree away from the COB and for ages younger than ca. 20 Ma. Finally, the statistical analysis of continental data assumed to be free of thinning or thickening effects suggests that the MCA should rather be considered in terms of “lowlands” and “highlands”. Relying on model of Moho depth, the “normal” crustal depth might be thinner than commonly accepted. In any case, the filtering of reduced topography can help to determine the impact of dynamic topography.

Objective estimates of mantle 3He in the ocean and implications for constraining the deep ocean circulation

Hydrothermal vents along the ocean's tectonic ridge systems inject superheated water and large amounts of dissolved metals that impact the deep ocean circulation and the oceanic cycling of trace metals. The hydrothermal fluid contains dissolved mantle helium that is enriched in 3He relative to the atmosphere, providing an isotopic tracer of the ocean's deep circulation and a marker of hydrothermal sources. This work investigates the potential for the 3He/4He isotope ratio to constrain the ocean's mantle 3He source and to provide constraints on the ocean's deep circulation. We use an ensemble of 11 data-assimilated steady-state ocean circulation models and a mantle helium source based on geographically varying sea-floor spreading rates. The global source distribution is partitioned into 6 regions, and the vertical profile and source amplitude of each region are varied independently to determine the optimal 3He source distribution that minimizes the mismatch between modeled and observed δ3He. In this way, we are able to fit the observed δ3He distribution to within a relative error of ∼15%, with a global 3He source that ranges from 640 to 850 mol yr−1, depending on circulation. The fit captures the vertical and interbasin gradients of the δ3He distribution very well and reproduces its jet-sheared saddle point in the deep equatorial Pacific. This demonstrates that the data-assimilated models have much greater fidelity to the deep ocean circulation than other coarse-resolution ocean models. Nonetheless, the modelled δ3He distributions still display some systematic biases, especially in the deep North Pacific where δ3He is overpredicted by our models, and in the southeastern tropical Pacific, where observed westward-spreading δ3He plumes are not well captured. Sources inferred by the data-assimilated transport with and without isopycnally aligned eddy diffusivity differ widely in the Southern Ocean, in spite of the ability to match the observed distributions of CFCs and radiocarbon for either eddy parameterization.