Oxygenation of the Baltoscandian shelf linked to Ordovician biodiversification

Marine biodiversity increased markedly during the Ordovician Period (~487–443 million years ago). Some intervals within the Ordovician were associated with unusually rapid and prominent rises in taxonomic richness, the reasons for which remain debated. Links between increased oxygenation and biodiversification have been proposed, although supporting marine oxygen proxy data are limited. Here we present an expansive multi-site iodine-to-calcium (I/Ca) record from Lower–Middle Ordovician marine carbonates in Baltoscandia that provides a detailed account of the spatio-temporal development of oxygen conditions across this palaeoshelf. The data document progressive oxygenation of regional seafloor environments, with well-oxygenated waters sourced from the palaeoequatorward Iapetus Ocean and peak I/Ca values (that is, dissolved oxygen concentrations) coinciding with the most pronounced biodiversity increases and ecosystem reorganizations during this time interval. This occurred while the climate cooled, global sea level dropped and carbonate deposits became regionally dominant. The results suggest that ventilation of shelves played a critical role in regulating early Palaeozoic marine biodiversity via development of ecospace and aerobic–metabolic conditions. Animal diversification coincided with increasing oxygenation of the Baltoscandian continental shelf from the Early to Middle Ordovician, according to iodine and calcium records.

(atmospheric) oxygen levels and long-term biodiversity patterns throughout the Phanerozoic [11][12][13] .It has been postulated that rising oxygen concentrations had a fundamental effect on marine biodiversity during the Ordovician [14][15][16] , but crude spatio-temporal resolution and scarce direct marine redox proxy investigations have hampered detailed interpretations.Oxygen levels in the atmosphere were likely relatively low during much of the early Palaeozoic, although there is discord between empirical data and geochemical models (and between models) 17,18 .The patchy and scant Ordovician geochemical data suggest generally lower oxygen availability in the shallow marine realm when compared to the Mesozoic-Cenozoic 18,19 .
In recent years, advances have been made in the development of geochemical proxies that help to elucidate ancient environmental conditions at various spatio-temporal scales.With this, reconstruction of the redox evolution of marine environments through geologic time has

Iodine-to-calcium variations across the Baltoscandian palaeoshelf
The study sections represent an Early-Middle Ordovician regional shelf setting on the mid-latitude (~45° S) palaeocontinent of Baltica (Fig. 1).The sections form a biostratigraphically well-constrained transect across this shelf, going from shore proximal to distal roughly towards the present-day south and west (Figs. 1 and 2 and Supplementary Tables 1-12) 25,26 .Generally, low I/Ca values characterize the outer shelf.The most consistently high I/Ca values occur in the ~mid-shelf area.Eastward, I/Ca values are also relatively low.Long-term stratigraphic trends and shorter-term 'excursions' in I/Ca are consistent across the studied region-thus, the first-order stratigraphic patterns in I/Ca provide a faithful record of oxygen availability on the Baltoscandian shelf during the Early-Middle Ordovician.Differences between sections can largely be explained by varying proximity to a local-regional oxycline, oceanic connection/restriction, water depth and current regimes across the shelf, but diagenetic alteration may locally play a role in some stratigraphic horizons (Supplementary Information sections 1-3).
The basal Floian (Early Ordovician) I/Ca values span a relatively large range (~2-5 μmol mol −1 ), indicating spatially variable oxygen conditions across the palaeoshelf.Higher values may be from very shallow-water environments related to the sedimentary break that commonly separates Floian strata from underlying rocks regionally (Fig. 2) [26][27][28] .The Floian thereafter is characterized by low I/Ca (~1-2), but values gradually increase towards a peak (~2.5-14.5) in the Floian-Dapingian (Early-Middle Ordovician) boundary interval.This suggests that water columns overlying Baltica transitioned from poorly to (near) well-oxygenated conditions through the Floian.Notably, the Floian-Dapingian transition, which appears to coincide with relatively low sea level, was associated with the expansion of carbonate sedimentation into many areas previously dominated by mudstone and shale [25][26][27] .
A pronounced decline in I/Ca occurs in the lowermost Dapingian strata, suggesting a region-wide decrease in water-column oxygen contents (Fig. 2).I/Ca values then mainly increase throughout the Dapingian.Consistently high I/Ca values (maximum ~5-7.5)occur become increasingly detailed.Iodine is redox-sensitive and among the earliest elements to be reduced under low-oxygen conditions in marine settings [19][20][21] .The dominant species of iodine in seawater are iodate (IO 3 − , oxidized) and iodide (I − , reduced).Briefly and simplified, iodate can be structurally incorporated into carbonates whereas iodide cannot.Therefore, iodine concentrations can be used to assess marine oxygen conditions with the proxy expressed as iodine-to-calcium ratios (I/Ca) in carbonates.Such ratios are commonly reported as I/(Ca+Mg) to accommodate varying carbonate mineralogy; this is also done herein, but for brevity we use I/Ca notation.Studies in modern marine environments have established that I/Ca ratios >2.6 μmol mol −1 indicate locally well-oxygenated conditions [19][20][21] .Lower I/Ca ratios suggest at least some oxygen present, and <1.5 μmol mol −1 characterizes oxygen-limited conditions 22 .Here we apply a 1.5-2.6 μmol mol −1 transitional redox interval as a general reference.Although the long residence time of iodine entails similar concentrations in seas globally 20 , in practice, I/Ca constitutes a local-regional proxy for the presence and proximity of oxygen-depleted waters.Hence, care should be taken before attempting to extrapolate data at a global scale 19 .The geographic limitation however has clear strengths in that I/Ca records conditions at the sampling site and the oxygen sensitivity of iodine makes it useful for establishing initial local redox responses.Other geochemical proxies, applicable to various lithologies, record global redox signals but are not necessarily reflective of local conditions 23 .Furthermore, such proxies typically respond to highly reducing conditions, limiting their capability to reconstruct subtle initial changes and variations close to the oxic-anoxic boundary across relatively small shelf areas, especially in carbonate settings 24 .
Regional-scale systematic analyses of I/Ca in shelf carbonates have the potential to provide a detailed and nuanced spatiotemporal understanding of water-column redox conditions in relation to the timing and rates of Ordovician biodiversification.Such studies may also shed light on the connections among climate, oceanographic conditions and the habitability of shelf environments in deep time.

Article
https://doi.org/10.1038/s41561-023-01287-zregionally in upper Dapingian through middle Darriwilian strata, signalling a protracted time interval (~5 Myr) of well-oxygenated conditions across the shelf.This interval of sustained high I/Ca, which is most pronounced in the mid-shelf area, coincides with evidence of a marked drop in sea level 8,25,29 .A temporary episode of declining I/Ca (minimum ~0.7) occurs near the Lenodus variabilis-Yangtzeplacognathus crassus conodont zonal boundary (Fig. 2).The upper half of the Darriwilian is characterized by a return to relatively low I/Ca values and purportedly oxygen-limited conditions (Fig. 2).The outer shelf saw a return to (black) shale deposition as sea level apparently rose [25][26][27]29 .
Map projections of our I/Ca datasets indicate a gradual spreading of well-oxygenated waters that moved into the interior of Baltica throughout the Early-Middle Ordovician (Fig. 3a-h and Supplementary Table 13).Oxygen-rich waters appear to have spread from the part of Baltica facing the Iapetus Ocean, roughly equatorward (Fig. 1) 30 .The regional I/Ca trends indicate that the mid-shelf area hosted the most consistently well-oxygenated environments, but the locus of apparent peak marine oxygen concentrations migrated through time.This may reflect prevailing shelf currents (Supplementary Table 14) and bathymetric conditions, which shifted with the paleogeographic configuration of Baltica and/or possibly with early stages of the Caledonian Orogeny 25,26,30 (Supplementary Information section 2).

Shelf oxygenation, biotic and environmental changes
The increase in oxygen availability indicated by I/Ca from the latest Floian through middle Darriwilian coincides with a successive series

Article
https://doi.org/10.1038/s41561-023-01287-z of biodiversity 'bursts' and behavioural innovations among benthic, nektic and planktic organisms (Fig. 2).Dominant regional patterns of marine taxonomic richness closely resemble the trends within our I/Ca datasets.Notably, global graptolite diversity is remarkably similar in stratigraphic patterns compared to I/Ca throughout the entire study interval, suggesting that well-oxygenated upper ocean water columns along continental margins played an important role in driving biodiversification among early Palaeozoic zooplankton.The coincidence of major oceanic currents 31 (Fig. 1a) with the spatio-temporal patterns of I/Ca (Fig. 3) suggests that the source area of oxygen-rich waters was outside of Baltica.Thus, the Early-Middle Ordovician redox development on the Baltoscandian palaeoshelf was linked with events and processes at a broader, perhaps global, scale.Considering the relatively slow oxidation of iodide 19,22 , the I/Ca data indicate that widespread oxidizing conditions were prevalent in water masses in the vicinity of Baltica.There are at present no directly comparable I/Ca datasets from

Article
https://doi.org/10.1038/s41561-023-01287-zother parts of the world, but we note that our results agree with broadly correlative global redox proxy datasets 16,32 and models of atmospheric oxygen 14,17 that collectively suggest increasing oxygen during the late Floian-early Darriwilian and a decrease in the late Darriwilian (Fig. 2).This does however not imply that all parts of the global ocean experienced the same and synchronous development 33 .Peak I/Ca values at the Floian-Dapingian transition coincide with a regionally traceable hardground complex, colloquially termed 'Blommiga Bladet' (Sweden), 'Püstakkiht' (Estonia) and 'Steklo' (Russia) 26 .Innumerable bioerosional structures permeate the hardgrounds, with macroborings reaching unprecedented areal and volumetric density.The Floian-Dapingian boundary marks an important step in the global ichnofossil record as hardground-dwelling infaunal organisms became frequent at the macroscopic scale 34,35 .The pitted hardgrounds of Baltoscandia form a key example of this development.It is plausible that the combination of relatively slow sedimentation and elevated oxygen, which penetrated deeper into the seafloor, provided optimal conditions for behavioural experimentation (Fig. 3b,c).Potentially, an ambient marine oxygen threshold was passed, as younger strata are characterized by extensive bioturbation 25,28 .
The long-duration interval of high I/Ca in the upper Dapingianmiddle Darriwilian coincided with remarkable changes in the biotic and sedimentary records (Fig. 2).The corresponding time interval was characterized by a marked increase in taxonomic richness among marine organisms in the Baltoscandian region and high-resolution diversity patterns co-vary with trends in our I/Ca datasets (Fig. 2 and Supplementary Information section 2) 8,10,36,37 .As Ordovician biodiversification took off in earnest, the marine fauna of Baltoscandia flourished overall and regional seascapes became the scenes of richer and more varied ecosystems.Whereas preceding times were characterized by low-diversity, trilobite-dominated faunas (Fig. 3a), the Middle Ordovician saw a prominent influx of brachiopods, bryozoans, echinoderms and other sessile benthos and conodonts and molluscs 8,10,25,[27][28][29][37][38][39][40] (notably, large cephalopods 41 ; Fig. 3b-h). The furishing of organisms is reflected in the composition of the strata via a marked increase in the proportion of skeletal grains as rock-forming components 25,27,28,37 (Supplementary Fig. 1), also apparent at the global scale 42 .Arguably, the expansion of oxygenated waters across the shelf (Fig. 3b-h) played a fundamental role in this biotic development (below).
The distinct drop in I/Ca near the L. variabilis-Y.crassus conodont zonal boundary coincides with the so-called Täljsten interval in Sweden, which is characterized by unusual litho-and biofacies (Fig. 2) 25 .Microbial fabrics and structures (oncoids, stromatolites) occur locally in Baltoscandia and possibly also globally 43,44 .Mass occurrences of Sphaeronites sp.cystoids (extinct echinoderms) are common, and this near-monospecific macrofauna has been interpreted as opportunistic colonization during adverse environmental conditions 25 (Supplementary Information section 2).The overall characteristics of the 'Täljsten' are consistent with dominantly reducing conditions at the seafloor and suggest temporary restriction of local depositional environments.The re-establishment of relatively low I/Ca values in the upper Darriwilian may reflect the global environmental development, with widespread formation of black shales 32 .

Climate and oxygen catalysed Ordovician biodiversification
As sea level dropped globally through the late Floian-early Darriwilian, large areas of the Baltoscandian shelf probably became more influenced by hydrodynamic activity 25,28,39,44 that entailed more consistent and efficient mixing of atmospheric oxygen into marine environments.However, sea-level change alone does not suffice to explain the observed shift to well-oxygenated conditions in the study interval.This is because oxygen-rich waters appear to have encroached upon the shelf from the Iapetus Ocean (Fig. 3), but an off-shelf direction would be expected if oxygenation was a simple result of a shallowing water column and migration of a regional oxygen minimum zone 15 .Hence, it is more likely that the progressive oxygenation reflects global-scale oceanographic changes 8,33 .It has been suggested that cooling sea surfaces and sea-level fluctuations during the Middle Ordovician were related, with the global climate state transitioning into 'icehouse'-like conditions 7,8,10 .Climatic cooling provides a plausible background to the oxygenation history of the Baltoscandian shelf, as changes in marine current regimes and strengthening of thermohaline circulation would be expected from global ocean temperature decline 8,31 .The palaeogeographic position of Baltica probably played an important role 33 , as pervasive winds at mid-latitudes promoted mixing of atmospheric gasses into the regional shelf environment (Fig. 1a).
Cooling climate has repeatedly been promoted as a central driver behind Ordovician biodiversification 7,8,10 , as have increasing oxygen levels 14 .These two interpretations are not mutually exclusive but rather complementary, as cooling would have allowed more dissolved oxygen in marine water columns.Indeed, the main trends in our I/Ca data are inverse to temperature estimates (Fig. 2 and Supplementary Information section 2).Climate and oxygen are intimately related through feedback processes at various spatio-temporal scales, and both strongly influence the distribution and organization of organisms [45][46][47] -thus, the combined effects of oceanic cooling and oxygenation provide a more holistic explanatory scenario for large-scale Ordovician biotic changes.A synergistic effect can be envisioned, wherein decreasing temperatures and rising oxygen levels expanded and stabilized the area and volume of habitable space via reduction of environmental and physiological hazards and stressors and enabled relatively energetic lifestyles 48 .Fundamentally, the net outcome was more balanced in terms of aerobic-metabolic requirements among the biota (Supplementary Information section 2).
Our study suggests that the Ordovician marine redox landscape was dynamic and evolving via multiple interacting mechanisms.It is likely that a combination of factors led to progressive radiation of marine life throughout the Ordovician, and spatio-temporal differences in biodiversity development at the global scale 1,4,5 paint a complex picture wherein dominant influences apparently varied between places and times.Nevertheless, at least on Baltica, several major biotic 'events' clearly coincide with I/Ca peaks and long-term trends, suggesting that oxygen indeed played an important role in evolutionary history during the Ordovician.Our I/Ca data are congruent with global models of oceanic oxygen conditions for the studied time interval 33 , indicating clear linkages to Ordovician climate 9 and large-scale oceanic circulation 31 .Although seascapes differ, the inferred connections between climate and marine oxygenation in the Ordovician provide a glimpse into the future and the fate of modern shelf ecosystems as global temperatures rise.

Online content
Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at https://doi.org/10.1038/s41561-023-01287-z.
Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this license, visit http://creativecommons. org/licenses/by/4.0/.© The Author(s) 2023

Fig. 1 |
Fig. 1 | Palaeogeographic characteristics of the study area.a, Palaeogeographic reconstruction for the Middle Ordovician (early Darriwilian, ~466 Ma), with names of notable continental masses and intercontinental seas and modelled upper ocean circulation 31 around Baltica indicated.Rectangle indicates area figured in b.The approximate distribution of land areas (tan-coloured shading) is generalized and simplified.b, Map of the Baltoscandian region, indicating the study localities (red dots) and notable features of the palaeoshelf as determined from relatively stable facies differentiation 8,25 .The areal distribution and extent of shelf features varied

Fig. 3 |
Fig. 3 | Spatio-temporal development of relative oxygen conditions in Baltoscandia.a-h, Time series illustrating the spatio-temporal development of relative oxygen conditions (I/Ca) across the Baltoscandian shelf.a, Stippled lines indicate generalized shelf features (Fig. 1b) for reference.Opaque red circles indicate localities included in data processing and map projections.Arrows indicate dominant water current directions, based on cephalopod
17,37versity.All datasets scaled vertically according to biozonation (Supplementary Information section 3) and relative to the Geologic Time Scale 2020 (ref.49),regionalstages(R.S.) after Nielsen et al.26.Ma stands for million years ago; Ep. stands for epoch.Horizontal dashed lines indicate global Age (Stage) boundaries, for lateral reference; vertical dashed lines in a-d indicate transitional redox interval (I/Ca 1.5-2.6 μmol mol −1 ; see main text).More positive I/Ca values indicate more dissolved oxygen in the marine environment.The data indicate that oxygen increased through the late Floian and into the early mid-Darriwilian, concurrently with major increases in marine biodiversity.All I/(Ca+Mg) data reported in μmol mol −1 ; errors for individual measurements are smaller than the data symbols.a,I/(Ca+Mg),outershelf(Fig.1b).b,I/(Ca+Mg),middleshelf.A generalized interpretation of the main I/Ca trends indicated by shaded curve (tan-coloured field); the precision of this qualitative approximation is limited by variations in data density.c,I/(Ca+Mg),middle-innershelf.d,I/(Ca+Mg),innershelf.One 'outlier' datapoint (14.5 μmol mol −1 ) was omitted.e,Generalizedsea-levelcurve(redline)showingpresumedeustaticvariations;compositeafterNielsen29,Rasmussenet al.8, Lindskog and Eriksson 25 , Rasmussen et al.10and Lindskog et al.28,37, excluding high-resolution local variability.f,Globalsea-surfacetemperature(SST)curve(coloured solid field) based on oxygen isotopes (δ 18 O), modified from lowest-temperature decile (that is, highest δ18O decile, presumed to be the diagenetically best preserved signal) reconstruction of Goldberg et al.9and atmospheric oxygen (O 2 ) modified from Krause et al.17.g,