Strontium Isotope Stratigraphy

Strontium isotope stratigraphy (SIS) is a high-precision chemostratigraphic and geochronological method that determines the numerical age of marine sedimentary-rock strata and reconstructs the paleoecology and geographic mobility of fossil organisms by measuring the ratio of radiogenic strontium-87 (⁸⁷Sr) to stable non-radiogenic strontium-86 (⁸⁶Sr). Developed into a quantitative global dating tool through the foundational Phanerozoic seawater compilations of William H. Burke et al. (1982) and improved with fifth-decimal-place thermal ionization mass spectrometry (TIMS) and multi-collector inductively coupled plasma mass spectrometry (MC-ICP-MS) by John M. McArthur, Richard J. Howarth, and Graham A. Shields (the LOWESS calibration curve, 2001–2020), strontium isotope stratigraphy relies on two geochemical principles. First, ⁸⁷Sr is continuously generated over geological time by the β⁻ radioactive-decay of naturally occurring rubidium-87 (⁸⁷Rb) with a half-life of 4.88 × 10¹⁰ years (decay constant λ = 1.42 × 10⁻¹¹ yr⁻¹). Second, the residence time of dissolved strontium (Sr²⁺) in the global ocean (~2.0 to 2.5 × 10⁶ years) is more than three orders of magnitude longer than the global ocean thermohaline mixing time (~1.0 to 1.5 × 10³ years).

Because oceanic mixing is rapid compared to input and removal fluxes, open-marine seawater at any single instant in geological history possesses a globally uniform ⁸⁷Sr/⁸⁶Sr ratio across all latitudes and depths (modern seawater ⁸⁷Sr/⁸⁶Sr = 0.709175 ± 0.000002). Over millions of years, however, the marine ⁸⁷Sr/⁸⁶Sr ratio fluctuates between ~0.7068 and ~0.7092 as plate tectonics alters the balance between two primary strontium reservoirs: high-⁸⁷Sr/⁸⁶Sr riverine runoff from ancient continental crust (~0.711 to >0.730) and low-⁸⁷Sr/⁸⁶Sr hydrothermal exchange at mid-ocean ridge basalt (MORB) spreading centers (~0.7028 to 0.7035). When marine organisms precipitate calcium carbonate (CaCO₃) shells or calcium phosphate bioapatite [Ca₁₀(PO₄)₆(OH,F)₂] teeth and bone, divalent Sr²⁺ substitutes directly for Ca²⁺ in the crystal lattice without measurable biological or kinetic mass-dependent fractionation after instrumental normalization to ⁸⁶Sr/⁸⁸Sr = 0.1194, locking the ambient seawater (or local continental bedrock) ⁸⁷Sr/⁸⁶Sr signature into the fossil record.

Rubidium-87 Decay and Crust–Mantle Differentiation

Strontium (Z = 38, alkaline earth metal) has four naturally occurring stable isotopes: ⁸⁴Sr (0.56%), ⁸⁶Sr (9.86%), ⁸⁷Sr (7.00%), and ⁸⁸Sr (82.58%). Among these, only ⁸⁷Sr is radiogenic, produced by the β⁻ decay of the alkali metal ⁸⁷Rb (Z = 37):

⁸⁷Rb → ⁸⁷Sr + β⁻ + ν̄_e + Q (0.283 MeV)

Consequently, the present-day abundance of ⁸⁷Sr in a closed geological reservoir of age t depends on its initial ratio (⁸⁷Sr/⁸⁶Sr)₀ and its parent-to-daughter elemental ratio (⁸⁷Rb/⁸⁶Sr):

(⁸⁷Sr / ⁸⁶Sr)_t = (⁸⁷Sr / ⁸⁶Sr)₀ + (⁸⁷Rb / ⁸⁶Sr) · ( e^λ*t* − 1 ) ≈ (⁸⁷Sr / ⁸⁶Sr)₀ + (⁸⁷Rb / ⁸⁶Sr) · λ · t

During partial melting of the Earth's upper mantle, the larger ionic radius of Rb⁺ (1.52 Å, substituting for K⁺ in biotites, muscovites, and K-feldspars) makes rubidium more incompatible in mantle minerals than Sr²⁺ (1.18 Å, which substitutes for Ca²⁺ in plagioclase and clinopyroxene). As a result, granitic continental crust is enriched in rubidium (high Rb/Sr ≈ 0.15–0.50) and evolves high radiogenic ⁸⁷Sr/⁸⁶Sr ratios over billions of years (mean continental river flux ≈ 0.7119; Himalayan silicate rivers reaching 0.725–0.740), whereas the depleted upper mantle retains a low Rb/Sr ratio (~0.01) and preserves unradiogenic ⁸⁷Sr/⁸⁶Sr values (~0.7025–0.7035).

The Three-Component Seawater Mixing Equation

At steady state, the secular ⁸⁷Sr/⁸⁶Sr ratio of global seawater (R_SW) reflects the flux-weighted mass balance of three primary inputs—continental riverine weathering (J_riv ≈ 3.3 × 10¹⁰ mol Sr/yr, R_riv ≈ 0.7119), mid-ocean ridge hydrothermal circulation and basalt alteration (J_hyd ≈ 1.5 × 10¹⁰ mol Sr/yr, R_hyd ≈ 0.7035), and diagenetic carbonate recrystallization/benthic reflux (J_diag ≈ 0.3 × 10¹⁰ mol Sr/yr, R_diag ≈ 0.7084):

R_SW = [ J_riv · R_riv + J_hyd · R_hyd + J_diag · R_diag ] / [ J_riv + J_hyd + J_diag ]

During intervals of supercontinent assembly and continental collision—such as the Early Paleozoic Pan-African/Taconic orogenies and the Cenozoic Himalayan–Tibetan uplift (40–0 Ma)—enhanced chemical weathering of old, radiogenic felsic silicates elevates J_riv and R_riv, driving the seawater ⁸⁷Sr/⁸⁶Sr curve steeply upward. Conversely, during intervals of rapid seafloor spreading, flood-basalt volcanism (such as the Late Permian Siberian Traps or Jurassic–Cretaceous superplumes), and extensive epeiric marine submergence of continents, high hydrothermal fluxes J_hyd pull seawater ⁸⁷Sr/⁸⁶Sr down toward mantle minima (~0.7068 in the Late Permian and Oxfordian Late Jurassic).

Instrumental Mass Fractionation Correction (⁸⁶Sr/⁸⁸Sr = 0.1194)

A decisive advantage of ⁸⁷Sr/⁸⁶Sr over light stable isotopes (such as δ¹⁸O or δ¹³C) is that vital effects, water temperature, and trophic fractionation are eliminated during mass spectrometry. Although biological precipitation and thermal ionization inside a mass spectrometer induce slight mass-dependent fractionation between ⁸⁷Sr and ⁸⁶Sr, every measurement simultaneously records the invariant natural ratio of the two non-radiogenic isotopes ⁸⁶Sr/⁸⁸Sr = 0.1194. Applying Russell's power-law mass-fractionation equation (R_true = R_meas · (m₈₇/m₈₆)^f) during TIMS or MC-ICP-MS acquisition normalizes out both natural biological fractionation and instrumental mass bias to better than ±0.000010 (±10 ppm), calibrated against the international carbonate reference standard NIST SRM 987 (accepted value ⁸⁷Sr/⁸⁶Sr = 0.710248).

Substrate Selection and Diagenetic Alteration

Reliable strontium isotope stratigraphy requires fossil substrates that retain their primary marine ⁸⁷Sr/⁸⁶Sr composition without post-depositional exchange with meteoric pore fluids. Low-magnesium calcite (LMC) shells—specifically articulate brachiopods (secondary fibrous shell layer), belemnite rostra, inoceramid and ostreid bivalves, and planktic foraminifera—alongside conodont albid crowns (composed of densely crystalline francolite bioapatite), provide the highest resistance to diagenetic alteration. By contrast, primary aragonite (corals, ammonites, gastropods) and high-magnesium calcite (echinoderms) contain higher initial Sr²⁺ concentrations (1,000–8,000 ppm) but rapidly dissolve and recrystallize into secondary diagenetic calcite in the presence of crustal groundwater. Prior to isotopic dissolution, carbonate samples are screened using cathodoluminescence microscopy (non-luminescent primary calcite vs. bright orange Mn²⁺-activated diagenetic luminescence) and trace-element thresholds (typically requiring [Mn] < 100 ppm, [Fe] < 250 ppm, and high Sr/Mn ratios).

The temporal resolution of strontium isotope stratigraphy is inversely proportional to the absolute slope |*d*(⁸⁷Sr/⁸⁶Sr)/*dt*| of the seawater curve (LOWESS look-up table) and directly proportional to analytical precision. During intervals of steep, unidirectional change—such as the Neogene–Quaternary (35–0 Ma), where ⁸⁷Sr/⁸⁶Sr climbs monotonically from 0.7078 to 0.709175 at ~0.000040 per Myr, or the Late Cretaceous (Campanian–Maastrichtian) and Early Jurassic (Sinemurian–Pliensbachian)—SIS achieves numerical dating precisions of ±0.1 to ±0.5 Myr, correlating isolated shallow-water carbonate platforms, evaporite basins, and phosphorites that lack index ammonites or volcanic ash beds.

In vertebrate paleontology, archaeology, and ecology, because terrestrial soils and freshwater watersheds inherit the specific ⁸⁷Sr/⁸⁶Sr signature of their underlying bedrock geology (isoscapes), measuring ⁸⁷Sr/⁸⁶Sr via laser-ablation MC-ICP-MS across sequential growth increments of tooth enamel (which has <1% organic content, large hydroxyapatite crystals, and high resistance to diagenesis compared to porous dentine and bone) reconstructs seasonal migrations of Pleistocene mammoths and dinosaurs, hominid foraging territories, and the marine versus freshwater osmoregulation habitats of extinct temnospondyli, mosasaurs, and spinosaurids.

Why is the ⁸⁷Sr/⁸⁶Sr ratio of the ocean identical everywhere on Earth at any given time?

The global homogeneity of marine ⁸⁷Sr/⁸⁶Sr arises from the ratio between the oceanic residence time of strontium (~2.5 million years) and the global ocean mixing time (~1,500 years via thermohaline circulation). Because an atom of strontium circulates through the world's ocean basins more than 1,500 times before being removed by carbonate sedimentation, local inputs from rivers or hydrothermal vents are thoroughly homogenized across all open-ocean basins. Only severely restricted, brackish epicontinental seas or estuaries with <10% marine water exchange deviate from the global seawater curve.

Why don't temperature or metabolic "vital effects" alter the ⁸⁷Sr/⁸⁶Sr ratio of shells and teeth?

While biological metabolism and water temperature do cause tiny mass-dependent isotope fractionations (roughly 0.02% to 0.05% per atomic mass unit between ⁸⁶Sr, ⁸⁷Sr, and ⁸⁸Sr), mass spectrometers simultaneously measure the non-radiogenic pair ⁸⁶Sr/⁸⁸Sr. Because ⁸⁶Sr and ⁸⁸Sr are neither created nor destroyed by radioactive decay and exist in a fixed natural abundance ratio of 0.1194, normalizing every measured ⁸⁷Sr/⁸⁶Sr ratio to ⁸⁶Sr/⁸⁸Sr = 0.1194 using Russell's power-law mass-fractionation correction mathematically strips away all biological, temperature, and instrumental fractionation.

Why is tooth enamel preferred over bone or dentine when reconstructing the migration of fossil vertebrates?

Cortical bone and tooth dentine contain ~20–30% organic collagen and tiny, plate-like bioapatite crystallites (20–50 nm long) perforated by microscopic vascular canals and dentinal tubules. After burial, groundwater readily infiltrates bone and dentine, precipitating secondary minerals and resetting their ⁸⁷Sr/⁸⁶Sr ratio toward the burial soil. By contrast, tooth enamel is >96% mineral by weight, non-porous, avascular, and composed of large, tightly packed hydroxyapatite prisms up to 1,000 nm long that lock in the animal's juvenile dietary ⁸⁷Sr/⁸⁶Sr signature against diagenesis.

1. Burke, W. H., Denison, R. E., Hetherington, E. A., Koepnick, R. B., Nelson, H. F., & Otto, J. B. (1982). "Variation of seawater ⁸⁷Sr/⁸⁶Sr throughout Phanerozoic time." Geology, 10(10), 516–519. https://doi.org/10.1130/0091-7613(1982)10<516:VOSSTP>2.0.CO;210<516:VOSSTP>2.0.CO;2)
2. Veizer, J., Ala, D., Azmy, K., Bruckschen, P., Buhl, D., Bruhn, F., Carden, G. A. F., Diener, A., Ebneth, S., Godderis, Y., Jasper, T., Korte, C., Pawellek, F., Podlaha, O. G., & Strauss, H. (1999). "⁸⁷Sr/⁸⁶Sr, δ¹³C and δ¹⁸O evolution of Phanerozoic seawater." Chemical Geology, 161(1–3), 59–88. https://doi.org/10.1016/S0009-2541(99)00081-900081-9)
3. McArthur, J. M., Howarth, R. J., & Bailey, T. R. (2001). "Strontium isotope stratigraphy: LOWESS Version 3: Best fit to the marine Sr-isotope curve for 0–509 Ma and accompanying look-up table for deriving numerical age." The Journal of Geology, 109(2), 155–170. https://doi.org/10.1086/319243
4. Bentley, R. A. (2006). "Strontium isotopes from the Earth to the archaeological skeleton: a review." Journal of Archaeological Method and Theory, 13(3), 135–187. https://doi.org/10.1007/s10816-006-9009-x