New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology
Flourishing ocean drives the end-Permian marine mass extinction
Edited by Timothy W. Lyons, University of California, Riverside, CA, and accepted by the Editorial Board June 30, 2015 (received for review February 23, 2015)

Significance
This study provides geochemical evidence that links some of the most important finds associated with the end-Permian mass extinction, including climate warming, enhanced weathering, increased primary productivity, and widespread marine anoxia under a common denominator: the linked biogeochemical sulfur and carbon cycles. Lethal marine conditions are likely the result of climate feedback mechanisms acting to increase nutrient input to the ocean, thereby stimulating global organic carbon production. With future projected climate change in mind, such climate feedback mechanisms could induce widespread eutrophication and expansion of anoxic and sulfidic zones, thereby fundamentally altering marine ecosystems.
Abstract
The end-Permian mass extinction, the most severe biotic crisis in the Phanerozoic, was accompanied by climate change and expansion of oceanic anoxic zones. The partitioning of sulfur among different exogenic reservoirs by biological and physical processes was of importance for this biodiversity crisis, but the exact role of bioessential sulfur in the mass extinction is still unclear. Here we show that globally increased production of organic matter affected the seawater sulfate sulfur and oxygen isotope signature that has been recorded in carbonate rock spanning the Permian−Triassic boundary. A bifurcating temporal trend is observed for the strata spanning the marine mass extinction with carbonate-associated sulfate sulfur and oxygen isotope excursions toward decreased and increased values, respectively. By coupling these results to a box model, we show that increased marine productivity and successive enhanced microbial sulfate reduction is the most likely scenario to explain these temporal trends. The new data demonstrate that worldwide expansion of euxinic and anoxic zones are symptoms of increased biological carbon recycling in the marine realm initiated by global warming. The spatial distribution of sulfidic water column conditions in shallow seafloor environments is dictated by the severity and geographic patterns of nutrient fluxes and serves as an adequate model to explain the scale of the marine biodiversity crisis. Our results provide evidence that the major biodiversity crises in Earth’s history do not necessarily implicate an ocean stripped of (most) life but rather the demise of certain eukaryotic organisms, leading to a decline in species richness.
The end-Permian mass extinction (EPME) is marked by ∼80% marine biodiversity loss (1). This event is linked with turmoil in the biogeochemical carbon and sulfur cycles (2⇓–4), alongside evidence for abrupt climate change and widespread euxinic (free H2S) and anoxic water column conditions (5⇓–7). Climate feedback mechanisms might have affected the biogeochemical cycles and may have spawned large-scale adverse conditions detrimental for many organisms (8). Regionally increased primary productivity (9⇓–11) and successive enhanced organic matter degradation by microbial sulfate reduction (MSR), which produces H2S as a by-product, has been held responsible for euxinic conditions during this time interval (12, 13). Even at a micromolar level, hydrogen sulfide is poisonous for most eukaryotic organisms, as it obstructs mitochondrial energy metabolism (14). The sulfate-reducing prokaryotes that drive this H2S production are a key catalyst of the biogeochemical sulfur cycle (15) and are associated with a cell-internal reaction pathway connected with kinetic isotope fractionation for both sulfate oxygen and sulfur (16, 17). This biochemical transformation of sulfate linked with an additional oxidative sulfur subcycle, mainly as microbial sulfur disproportionation, and sedimentary sulfide burial can leave a distinct imprint on the isotope composition of the marine sulfate pool (15, 18). In addition, MSR is intimately tied to organic carbon (OC) availability and is the most important contributor to sedimentary OC remineralization on the continental shelf, after aerobic respiration, and thereby fuels the biogeochemical carbon cycle (15, 19, 20). Measuring sulfate sulfur and oxygen isotopes in marine sediments can provide important information about biologically induced sulfur transmissions among different reservoirs (21⇓–23) and can record a global signal (15, 23).
We produced a new high-resolution dataset of carbonate-associated sulfate (CAS) sulfur and CAS oxygen isotopes spanning the time interval of the EPME. Changes in the marine dissolved sulfate isotope composition are considered to be faithfully captured by CAS (4, 15, 23), and chromium-reducible sulfur (CRS) records the products of MSR (15, 24). Carbonate rock was collected at the Kuh-e-Ali Bashi 1 (38.940°N, 45.520°E) and Zal (38.733°N, 45.580°E) sections of the Julfa region (NW Iran, Fig. 1), which were situated on the north-northeastern margin of the NW Iranian Terrane, belonging to the Cimmerian Microcontinent during the latest Permian (ref. 25, Fig. 1). These sections represent an excellent stratigraphic correlation and age assignment (Timescale and Figs. S1 and S2). Carbonate rock successions spanning the Permian−Triassic (P−Tr) boundary interval in this area can be divided into two major lithological units: the Ali Bashi Formation made up of the Paratirolites Limestone Member marking the upper part and the Zal Member marking the lower part (Fig. S1). The Aras Member (26), or “Boundary Clay,” and the following limestone unit, belong to the Elikah Formation (Fig. S1). The Zal Member consists, predominantly, of shales with some intercalations of limestone and marls. The red nodular Paratirolites Limestone beds are classified as argillaceous limestone. Macrofossil assemblages consist of rare ammonoids, nautiloids, rugose corals, and brachiopods. This unit is conformably overlain by the Aras Member, which records a marked reduction of CaCO3 (26). In addition, the Permian faunal elements disappear in this unit; therefore the boundary between these units is considered the extinction horizon (ref. 26, Fig. S1). On top of the clay-rich Aras Member, the remainder of the Elikah Formation consists of limestone beds. These lithologies were deposited as a deep-water carbonate platform and did not undergo large-scale sea level changes (Facies Description). Therefore, this carbonate platform with an isolated location in the Paleotethys (Fig. 1) bears a large potential for obtaining secular changes of ocean chemistry across the EPME (Fig. 1). Stratigraphic patterns in these CAS sulfur and oxygen isotope records allow for a better understanding of the sulfur cycle across the EPME, thus enabling the assessment of its potential impact on biodiversity.
Stratigraphy and geography. (A) Isotopic composition of CAS (circle) and CRS (diamond) of carbonates from Kuh-e-Ali Bashi (white) and Zal (gray). Significant changes in these geochemical records, discernable in the blue lines (loess fit), mark the extinction horizon and suggest a substantial perturbation of the sulfur cycle. Error bars denote reproducibility (2σ; Tables S1 and S2). The temporal framework is constructed by ages for conodont biozones (Figs. S1 and S2 and Timescale): 1, C. Orientalis; 2, transition; 3, C. subcarinata; 4, C. chanxingensis; 5, C. bachmanni; 6, C. nodosa; 7, C. yini; 8, C abadehensis; 9, C. hauschkei; 10, H. praeparvus-M. ultima; 11, H. parvus; 12, H. lobota; 13, I staeschei; and 14, I. isarcica. The size of individual CAS data points stands for the confidence of representing true seawater chemistry: large, high confidence; small, low confidence (high CRS and non-CAS); medium, intermediate confidence (high CRS or non-CAS; Fig. S3, Tables S3 and S4, and Reliability of the CAS Proxy). Geographic location of the sections (B) today and (C) during the P−Tr (modified from ref. 25).
Stratigraphic column of the (A) Kuh-e-Ali Bashi 1 and (B) Zal sections. Lithological column and stratigraphic position of δ34SCRS (diamond), δ34SCAS (circle), and δ18OCAS (circle). Conodont zones: a, C. bachmanni; b, C. nodosa; other zones follow from C. (C) Characteristic Late Permian–Early Triassic conodonts from the Julfa region. (Scale bars, 100 μm.) All specimens are stored in the collection of the Ferdowsi University of Mashhad (overview in refs. 26 and 71); c, Clarkina yini, FUM# G258-6, Ali Bashi Formation, Main Valley section, upper view; d, Clarkina abadehensis abadehensis, FUM# 1J248-6, Ali Bashi Formation, Ali Bashi 1 section, upper view; e, Clarkina hauschkei, FUM# 1J249D-9, Ali Bashi Formation, Ali Bashi 1 section, upper view; f, Hindeodus praeparvus, FUM# G274-6, Elikah Formation, Main Valley section, lateral view; g, Hindeodus changxingensis, FUM# 1J252-11, Elikah Formation, Ali Bashi 1 section, lateral view; h, Stepanovites ?mostleri, M element, FUM# 1J253-16, Elikah Formation, Ali Bashi 1 section, lateral view; i, Merrillina ultima, Pa element, FUM# 4J203-7, Elikah Formation, Ali Bashi 4 section, lateral view’ j-1, Hindeodus parvus, FUM# G278-17, Elikah Formation, Main Valley section, lateral view; j-2, Hindeodus parvus, FUM# AJ206-2, Elikah Formation, Aras Valley section, lateral view; k-1, Hindeodus lobata, FUM# 4J217-5, Elikah Formation, Ali Bashi 4 section, upper view; k-2, Hindeodus lobata, FUM# 1J265-4, Elikah Formation, Ali Bashi 1 section, lateral view; i-1, Isarcicella staeschei, FUM# 1J270-4, Elikah Formation, Ali Bashi 1 section, lateral view; i-2, Isarcicella staeschei, FUM# 4J222-7, Elikah Formation, Ali Bashi 4 section, upper view; j-1, Isarcicella isarcica, FUM# 1J281-1, Elikah Formation, Ali Bashi 1 section, upper view; and j-2, Isarcicella isarcica, FUM# 4J225-4, Elikah Formation, Ali Bashi 4 section, upper view.
Age model for the P−Tr sections at Iran. Note the high sedimentation rates for the Paratirolites Limestone unit when relying on age-constrained conodont zones of the latest Permian (dotted line). An alternative age model for the latest Permian—used in this study to time-constrain the geochemical data—represents the sedimentation rate calculated for the Changhsingian nodular limestones of the Abadeh section (solid line). This latter sedimentation rate has been calculated by interpolated dates of the extinction horizon (EH) and the base of the Changhsingian system of the Abadeh section.
Facies Description
The successive lithological change from dominantly clastic to nodular Paratirolites Limestone has been related to drowning of the carbonate platform (52, 53). This notion is enforced by lithological features such as frequent hardgrounds as well as the distinctive red color of the latter unit—plausibly related to authigenic iron accumulation and suggesting a low rate of carbonate accumulation (53). The Zal Member in contrast to the Paratirolites Limestone Member probably had a higher siliciclastic input and resultantly higher sedimentation rates. Deposition of these units occurred at a depth below the storm wave base in a well-oxygenated environment by a productive skeletal carbonate factory (53, 54). The Aras Member has been assigned to increased terrigenous input, based on the thickness of the unit (52, 54), expressing high sediment accumulation rates (55). A concomitant reduction in carbonate production could have been equally important (52, 53). However, the water depth of deposition was still likely to have been below the storm wave base in an oxic environment, based on ostracod assemblages (55). Somewhat higher CRS values in this unit might relate to a higher local production of sulfide by MSR (Tables S3 and S4). Increased pyrite content, or other iron sulfides, in the sediment can be a function of lower oxygenation levels or higher OC availability, but the two parameters are likely linked (22). The sulfide oxygen front at these sites was most likely located within the sediments, as a high-resolution microfacies analysis (54) did not reveal framboidal pyrite, a sign of water column sulfide (31). A continuation of deep-water conditions, especially for the nodular lithology, which resembles the preextinction strata—deep water facies—has been implied for the first few meters above the Aras Member at the Kuh-e-Ali Bashi site (54). The platy limestone beds of the Elikah Formation might, however, indicate a somewhat shallower water depth (54) but no high-amplitude sea level changes (52, 54). Lithological indices for pervasive anoxia and euxinia, such as framboidal pyrite or high OC accumulation, are absent. However, assemblages of ostracods together with abundant bivalve specimens (Claraia) might indicate slight or episodic dysoxic bottom water conditions (52, 54, 55).
Timescale
U−Pb geochronology has provided absolute dates for the chemostratigraphically and biostratigraphically well-resolved Meishan P−Tr section (56). Recent advances in U−Pb dating with isotopic dilution-thermal ionization mass spectrometery and especially the use of a more precisely calibrated tracer solution resulted in dates with analytical uncertainties of less than 100 ky (57). We used these radiometric dates to make linear interpolations for conodont zones of the Meishan P−Tr section (58). Notably, the top of Hindeodus parvus is 251.873 ± 0.004 Mya, and the first occurrence date (FAD) of Clarkina yini is 252.046 Mya. In addition to these interpolated dates, the FAD of H. parvus (251.902 ± 0.024 Mya), the extinction horizon [base bed 25 at the Meishan P−Tr section, 251.941 ± 0.037 Mya (57)], and the FAD of Clarkina wangi [254.150 Mya (56)] are also used in the age model. The linear interpolations have some inherent uncertainties, which potentially could be large, as the sediment accumulation rates, which are used, may not account for hiatuses in the Meishan P−Tr section. However, these dates, in our opinion, provide the best possibility for constructing a timeframe for the study sites in NW Iran.
Recently developed Late Permian high-resolution conodont zonation constructed by the sample population taxonomic approach for the Kuh-e-Ali Bashi as well as for Meishan P−Tr sections, together with some new data for the Early Triassic in Iran (Fig. S1), enables the correlation of these sections. However, the correlation is restricted to the extinction horizon and H. parvus zone only (55). The preextinction strata cannot be correlated based on conodont zonation, possibly due to the endemism of conodonts from the south Chinese basin (55). This is further supported by the resulting sedimentation rates for the Paratirolites Limestone Member and the Zal Member based on a conodont correlation scheme, which would contradict the facies interpretation of these units (Fig. S2 and Facies Description).
We, therefore, consider a sedimentation rate of 0.3 cm/ky for the Paratirolites Limestone Member (cf. ref. 52). Such a sedimentation rate can be calculated for a nodular limestone unit that spans the complete Changhsingian of the Abadeh section (central Iran; Fig. S2). In this way, the preextinction samples, belonging to the Paratirolites Limestone Member, together with the postextinction samples, can be placed into a temporal framework for the sites in NW Iran with a relatively high level of confidence.
Results and Discussion
To obtain a primary seawater chemistry signal, samples likely to be contaminated during rock diagenesis by allochtonous sulfur sources, such as leachable sulfate (non-CAS) and sulfide minerals, are identified (Reliability of the CAS Proxy and Fig. S3). This enables the construction of the primary stratigraphic CAS oxygen and sulfur isotopic record (Fig. 1). In contrast to previous CAS studies of this time interval (3, 4), the preextinction CAS isotope record (∼2 My) is remarkably stable around 17‰ for both δ34SCAS and δ18OCAS values with a spread of 1‰ and 2–3‰, respectively (Fig. 1). A ∼3‰ negative δ34SCAS excursion and a 5–7‰ positive shift in δ18OCAS can be traced for both sites during a ∼40 ky period following the marine extinction horizon. A positive δ34SCAS excursion is found above the conodont-defined P−Tr boundary (∼3 m; Fig. S1) with values ∼5‰ heavier compared with preextinction strata.
Relative concentrations and isotope composition of the extracted sulfur species. CRS (diamond) and leachable non-CAS (triangle) against CAS (circle) sulfur isotope values plotted against CAS oxygen isotope values of (A) Kuh-e-Ali Bashi 1 and (B) Zal sections. Most δ34S from CRS, non-CAS, and CAS have distinctive values. However, some CAS are marked by δ34S that plot in the region of non-CAS δ34S values, and are often linked to comparatively high CRS/CAS or non-CAS/CAS ratios (red dots). These anomalous values accompanied by lower δ18OCAS might represent oxidation of sulfide during diagenesis. On the other hand, some δ34SCAS, mainly the Zal section, are enriched as well as paired with high CRS/CAS or non-CAS/CAS ratios and high δ34SCRS, δ34Snon-CAS, and δ18OCAS and could therefore constitute postdepositional incorporation of CAS, from an evolved sulfate pool in a closed system. (C) CRS/CAS or (D) non-CAS/CAS ratios plotted against CAS sulfur isotopes. The lines represent static cutoff boundary values; samples to the right are presumed vulnerable to postdepositional alteration because relative high concentrations of CRS or non-CAS and could cause mixing of isotope signals.
Microbial sulfate reduction is an important aspect in the sulfur cycle, and the associated isotopic fractionation leads to 34S and 18O enrichments of the residual sulfate pool (4, 27). Most of the MSR-produced H2S will be reoxidized, and only a fraction is buried as sedimentary sulfide serving as a sink for 32S (22); some H2S can escape the sediment, especially under low O2 concentrations (8). Hydrogen sulfide release from the sediment is, however, also a factor of sulfide formation and reactive iron availability (8, 22). Continental oxidative weathering and mantle-derived sulfur are fluxes that have a depleted sulfur and oxygen isotopic composition compared with seawater sulfate (4, 24, 27). Evaporite formation and successive burial (with a negligible fractionation) is a sink of sulfate, but its importance depends largely on the extent of suitable environments for deposition [with net evaporation rates (28)]. In contrast to the isotope value of sulfate sulfur, δ18OSO4 is considerably affected by MSR and microbial sulfur disproportionation (29). Sulfur disproportionation produces sulfate with an oxygen isotope fractionation of ∼18‰ over the oxygen isotope composition of ambient sea water (29). In contrast to the sulfur isotope part of the sulfur cycle, this microbial-dominated sulfur subcycle is effectively a sink for 16O (4, 29). These different reaction pathways for the sulfur and oxygen ionic compounds of sulfate might lead to a dissimilar temporal evolution of marine sulfate sulfur and oxygen isotope composition, such as recorded by the deviating CAS δ34S and δ18O at the extinction horizon in the dataset (Fig. 1).
A depleted sulfate reservoir, in comparison with modern values of 28 mM, has been implied for allowing rapid changes of seawater δ34SSO4 (21), such as observed here for the EPME (>0.1‰⋅ky−1). The offset in δ34S between CAS and CRS (Δ34S) was often attributed to the kinetic isotope fractionation of MSR (15). A closed system of diagenetic sedimentary sulfide formation can lead to a small Δ34S due to ongoing sulfate depletion by MSR (15, 30), whereas a sulfate unlimited system, such as water column pyrite formation (e.g., in the modern Black Sea), can lead to a Δ34S of 45 ± 15‰ (15, 30, 31). The observed constant Δ34S (15−16‰) is invariant to burial rate changes (Fig. S2), and might suggest low sulfate levels within sedimentary porewater and the water column (cf. refs. 24 and 30).
We applied a box model calculation (Methods) to delineate the fluxes that determine the isotopic composition of the marine sulfate across the P−Tr transition. Paleontological and sedimentological evidence from the Iranian sections show no indications of pervasive water column anoxia (Facies Description); this is in contrast to the inferred Late Permian marine anoxia at sites in South China and Australia (13, 32). This spatial heterogeneous pattern can be best explained by increased nutrient influx only into parts of the Tethys. This influx caused regional eutrophication associated with increasing extents of oxygen minimum zones (33) and euxinic conditions (12) rather than deep bottom water anoxia caused by a stagnant water column. Widespread ocean fertilization and drawdown of water column O2 levels have been related to an increased influx of terrestrial material (12, 33) from soil erosion (34) and enhanced weathering (35) initiated by acid precipitation and destabilization of vegetation and soils (34), global warming, and an associated intensified hydrological cycle at that time (7). The climax of this biogeochemical cascade correlates with the greenhouse warming often associated with coeval Siberian Trap volcanism (2, 5). The importance of organic substrates for MSR justifies the incorporation of the OC pool and the link to continental weathering in our box model. This simplified approach is further justified because the continental supply of nutrients fertilizes the oceanic shelf connected with highest MSR activity and OC decomposition rates, which are almost equal in oxic and euxinic sediments (19, 20).
A best model fit, based on sensitivity experiments (Sensitivity Experiments and Fig. S4), suggests a 16-fold increase in the weathering flux that stimulates a >16-fold global marine OC inventory, starting at the extinction horizon and lasting for a period of 40 ky (Fig. 2). This stimulates large MSR activity (>17-fold) but also sulfur disproportionation (Methods), although it has been suggested that the latter could not explain an oxygen isotope fractionation of >18‰ between water and sulfate (29). An average δ18Oseawater of ∼−1‰ (6) and δ18OCAS excursion to values of >20‰ make sulfur disproportionation an unlikely candidate. Increased OC degradation by MSR is a strong and robust argument for the reconstructed scenario and explains the observed δ18OCAS excursion. Independent numerical exercises with higher complexity models suggest that >4 times enhanced continental weathering considerably increases the primary productivity in shelf settings (12, 36), followed by anoxic and sulfidic conditions as a consequence of increased (an)aerobic OC remineralization.
Isotopic composition of CAS (circle) and numerical solution when perturbing the system with a 16-fold increase in the global weathering flux after the extinction horizon (black line). Gray area represents range of sensitivity experiments (Fig. S4). Time-equivalent seawater temperature rise, seen in δ18Ophos (7), can enforce an enhanced hydrological cycle stimulating larger terrestrial weathering fluxes.
Sensitivity experiment for a volcanic- and weathering-induced sulfur cycle perturbation. (A) This scenario encompasses an influx of simultaneous isotopically depleted sulfur (δ34SSO4 = 3‰) and isotopically enriched oxygen flux (δ18O SO4 = 34‰) from a volcanic source after the EH. Shaded time interval stands for the perturbation. (B) The fraction of sulfide buried as pyrite is lowered 10 times for an interval of 20 ky after the EH, after that, the fraction of pyrite burial is increased 2 times for an equivalent time period. Shaded time interval stands for the perturbation. Continental weathering is maintained to 7 times the preextinction value from the P−Tr boundary until the end of the model run. Red dotted line stands for the scenario closely correlating with CAS isotope data.
The initial negative δ34SCAS excursion also suggests that pyrite burial did not increase with enhanced MSR activity, as is often assumed (4, 27). The main controlling factors for pyrite formation are labile OC, sulfide formation, oxygenation conditions, and highly reactive iron [Fe(HR)] availability in the sediment and water column (22). In sulfidic water column conditions, such as the modern Black Sea, which represents only a fraction of the total modern ocean volume (27), pyrite formation is a syngenetic process (Fig. 3) rather than a diagenetic process (31, 37). This pyrite sedimentation style depletes Fe(HR) already in the water column and largely inhibits the pyrite formation within the sediment (31). In such euxinic settings, pyrite formation becomes iron-limited and pyrite burial becomes sensitive to the supply of Fe(HR) to the ocean (38). Mechanisms for Fe(HR) supply are hydrothermal input, windblown dust, riverine input, and diagenetic iron cycling on the shelf (39, 40). Although the iron flux from the continent is potentially large, 70–90% of the Fe(HR) is already lost as particles in estuaries (40). The iron reaching the ocean is sensitive to water turbulence and grain concentrations delivered by rivers (40). Relating all this to the end-Permian scenario, a situation under a high sea level (5) and a sudden influx of soils (34) accompanied by high physical erosion rates (35) due to an intensified hydrological cycle in a warmer ocean (7) would create an environment that already depletes Fe(HR) before entering the ocean (Fig. 3).
Conceptual model for environmental changes across the P−Tr boundary, depicting the shelf settings that are prone to be affected in a scenario of intensified continental influx. Preextinction Late Permian marine communities would consist of multiple trophic levels. Authigenic pyrite formation is supported by sulfide production during anaerobic OC remineralization via MSR and nonlimiting iron availability. The extinction event marks an increased terrigenous influx causing local eutrophication as well as the input of isotopically depleted sulfur, aided by an intensified hydrological cycle forced by global warming. Increased sediment load and turbulence of river discharge flocculates iron as Fe oxy(hydr)oxides in estuaries, reducing the ocean Fe(HR) inventory. Consequential OC remineralization by MSR and Fe(HR) limitation raises seawater H2S. Sulfide toxicity alongside marine anoxia sustains a low-species richness but a high-biomass marine community. Pyrite sedimentation is limited to small amounts of syngenetic pyrite precipitation. Relaxation of physical continental weathering reduces ocean input of nutrients and raises Fe(HR) delivery to the postextinction ocean, thereby reducing sulfidic water column conditions and relocating pyrite sedimentation toward the sediments, where decreased O2 irrigation due to a less active benthic community stimulates authigenic pyrite formation. The enlarged OC pool is sustained due to carbon recycling by an active microbial community.
We do not provide an estimate for the pyrite flux change. However, under the Late Permian oceanic conditions, it was not necessary that pyrite burial increased concomitantly with increasing MSR. The diverging CAS sulfur and oxygen isotope pattern can be best explained by increased levels of MSR but decreased pyrite burial, compatible with an enlarged global extent of iron-limited euxinic seawater (Figs. 2 and 3). The presence of green sulfur bacteria—which thrive in H2S-rich surface water—during the time interval of the sulfur cycle perturbation, evidenced by biomarker analyses (32), suggests widespread euxinia. Moreover, ocean fertilization is in concert with an independent model for the Late Permian ocean, suggesting 10 times higher nutrient levels to induce globally pervasive water column H2S (12). The new CAS isotope data and the related environmental conditions suggest that the euxinic zones expanded suddenly at the extinction horizon (Fig. 3). Synchronicity of the two observations support the notion that sulfidic toxicity drove the marine biodiversity loss at the EPME, as first hypothesized by Kump et al. (8). However, the entire continental shelf area [∼5.8e6 km2 for the latest Permian (28)] would be primarily affected in our conceptual model, in contrast to only upwelling regions with a comparatively smaller total dimension [∼0.4e6 km2 for the modern ocean (8)]. Predominant sulfate reduction on the shelf, rather than in the deep sea, agrees with observations in the modern ocean of decreasing MSR with increasing distance to land (19, 22). In addition, a marine region would be targeted that is the locus of Paleozoic biological activity (the neritic realm), and thereby this conceptual model fulfills all criteria as suggested by Wignall and Twitchett (5) to explain a wholesale collapse of marine invertebrate diversity. However, the geography and nutrient load of weathering fluxes would also be an important determinant of the spatial distribution and severity of productivity-driven redox changes (oxic, nonsulfidic anoxic, and euxinic) (41). Regionally biased ecosystem destruction would leave refugia with less adverse conditions and would allow rapid recolonization of vacant ecospaces through migration of survivors. Postextinction species are frequently found to occur along a spectrum of environments, and within-habitat diversity is typically higher than diversity between habitats, a pattern explained by species that exploit the full range of the available niches under low competition pressure (42). Also, refilling of vacant ecospace by an ecologically diverse group of survivors can explain why, on a global scale, ecosystems remained operating even though species diversity collapsed (43).
This scenario to explain the EPME enhances the notion that life can influence seawater chemistry (22) and the relationship to OC remineralization (13). CAS isotopes provide a more direct view on biological productivity in past oceans, suggesting a flourishing of life. It is in stark contrast with the concept of decreased OC production, or “Strangelove” ocean, often implied in extinction scenarios (44). This study also emphasizes that, besides the property of organisms to construct a habitable planet (21), they can also act as a catalyst for destruction. However, marine life might have been very different with mainly prokaryotes dominating over eukaryotes (11, 45). A return to nearly preextinction δ34SCAS, within 40 ky, suggests recommencing pyrite burial; the disturbance was therefore a rather transient event, which agrees with the postulated short duration of soil erosion and local euxinic surface water conditions (32, 34). It argues for an enhanced supply of iron, possibly under reduced detrital loading of river discharge, a waning of the global extent of euxinic conditions and consequential increase of diagenetic pyrite accumulation, or a combination of both (Fig. 3).
The subsequent Early Triassic positive δ34SCAS is apparently not directly related to the biodiversity crisis (Fig. 1). Rather it records increased sequestration of sulfur by authigenic pyrite caused by the lack of organisms that normally irrigated sediments with O2 by burrowing (22) (Fig. 3). Reduced sediment mixing has been documented by the decreased thickness of the sedimentary mixed layer and the absence of burrowing faunas at localities globally, a possible symptom of the widespread devastation (5, 46). The relatively stable and high δ18OCAS in this period indicates that MSR remained important during the Early Triassic. Persisting high levels of MSR suggests a prevailing large OC pool, together with a long-lasting (∼800 ky) increase in continental weathering (35), a perturbed carbon cycle, and high global temperatures, the latter two returning in intervals throughout a period of almost 4 My (47, 48). These Earth surface processes are regarded as an important negative feedback loop of the carbon cycle, where enhanced production and sequestration of OC is stimulated by global warming and subsequent chemical weathering rates (36). The prolonged disturbance after the EPME contradicts a fast return (<100 ky) to predisturbance climate and carbon cycle, enforced by carbon sequestration, as suggested for other events marked by greenhouse warming, such as the Cenozoic hyperthermals (49). Degradation of organic substrates on the seafloor by sulfate-reducing microbes, possibly in conjunction with more-efficient water column carbon cycling of labile bacterial OC (45) and low fecal pellet production by reduced numbers of zooplankton (10), might all facilitate a sustained large marine OC pool. The postextinction prokaryote-dominated marine communities, in particular sulfate-reducing microbes, might therefore have had a crucial control on the Early Triassic carbon cycle and would affect marine redox conditions and climate with an adverse effect on the biotic recovery.
Reliability of the CAS Proxy
Sulfate is incorporated into the calcite lattice as an anion substitution with no apparent fractionation of the sulfur isotopes and therefore records the δ34S composition of the ambient seawater sulfate (23). The reliability of CAS oxygen isotopes in recording certain aspects of the global sulfur cycle has been acknowledged in a recent study (59). These authors, however, point out that questions remain about the influence of possible vital fractionation effects and specific microhabitats upon δ18OCAS. In contrast, other studies (4, 23) have shown that ancient CAS oxygen isotopes agree well with time-equivalent δ18O from sulfate locked in evaporites, suggesting that δ18OCAS can be used as a proxy for seawater chemistry. Multiple high-amplitude variations in the isotopic composition of CAS—thought to represent original seawater chemistry—have been used to explain scenarios describing the EPME (3). This large δ34SCAS variability, however, may also reflect postdepositional alteration (3).
CRS, leachable non-CAS, and CAS are marked with distinct isotopic values (Fig. S3). Based on this observation, we infer that leachable non-CAS and CRS could be potential sources to alter the original isotopic composition of CAS. Hence, to identify possible alteration of the CAS isotope signature, the relative abundance of CRS, leachable non-CAS, and CAS for each sample is helpful in this respect, more so, considering the evidence that initially large CAS concentrations buffer against diagenetic overprinting of the primary isotope signal (60). To approximate the likelihood that our obtained isotope values reflect ambient seawater sulfate δ34S and δ18O, concentrations of CRS and non-CAS sulfate relative to CAS are compared with δ34SCAS of each sample (Fig. S3).
Samples with comparatively high concentrations of CRS seem to be predisposed toward lower concentrations of δ34SCAS, whereas both depleted and enriched values dominate the sample population marked by a relatively high amount of leachable non-CAS (Fig. S3). It is noteworthy that depleted values of δ34SCAS, accompanied by high CRS/CAS or leachable non-CAS/CAS, are mainly found at an interval with a low carbonate content around the extinction horizon, especially for the Kuh-e Ali Bashi section (Figs. S1 and S3). It is therefore perceivable that these depleted δ34SCAS from low-carbonate samples, linked with low CAS concentrations, are in part an effect of oxidation of sulfide species, such as H2S or sulfide minerals during rock diagenesis (23). These likely altered δ34SCAS values lie in the same range as the low δ34S values of the CRS and leachable non-CAS fractions (Fig. S3). In addition, these low-carbonate rocks are also marked by the presence of some depleted δ18OCAS relative to values bracketing those (Fig. S1). This observation suggests a common diagenetic mechanism for these samples with paired low δ34SCAS and low δ18OCAS, for which oxidation of sulfide is the most likely one. This chemical reaction mainly uses H2O as an oxygen source, with an oxygen isotope value close to that of the ambient seawater (23).
On the other hand, some of the samples, mainly from the Triassic part of the Zal section, with relatively high leachable non-CAS concentrations are marked by enriched δ34SCAS. These enriched values could be indicative of an isotopic evolution reflecting progressing MSR in an evolving sulfate pool approaching closed system conditions (23). This would also agree with the often relatively heavy concentrations of δ18OCAS and δ34S from leachable non-CAS and CRS, respectively (e.g., Zl 4.15 m and 10.05 m, Tables S3 and S4). Enriched sulfate δ34S produced by isotope fractionation associated with MSR might be incorporated during authigenic carbonate formation.
As such, we argue that the ratios of CRS to leachable non-CAS with CAS provide a good impression of the likelihood that an actual seawater sulfate isotope signal is recorded in CAS. Based on the density distribution (Fig. S3), CRS/CAS of 0.1 and non-CAS/CAS of 0.4 are used as static cutoff boundaries for the herein-presented CAS isotope values. The resulting time series for δ34SCAS and δ18OCAS, with the highest probability of capturing seawater sulfate, show nearly identical isotope patterns (Fig. 2).
Steady-State P−Tr Sulfur Cycle Model
Based on the measured preextinction δ34SCAS, the observed lower background value for Late Permian marine sulfate sulfur isotopic composition is 17‰ (Fig. 1). Low sulfur isotope marine sulfate reflects either a change in sulfur input (specifically the input of sulfur with a lower isotopic composition) or reduced global pyrite burial (61).We hypothesize that exposed Carboniferous and Permian strata on the continents were rich in isotopically depleted pyrite sulfur (62). Hence, to achieve steady state, we solved Eq. S1 for δw, assuming that ΔMSR is −35‰ (24, 63), and taking Holocene values for the other parameters (Table S5).
Sensitivity Experiments
By conducting sensitivity experiments, we determined if the solutions with a best fit to the data are unique or can be explained in alternative scenarios. A critical aspect of the Permian and Triassic is the lowered marine sulfate levels with respect to modern oceanic levels. A low oceanic sulfate reservoir can explain sulfate isotope excursions constrained to short time intervals—shorter than the residence time of marine sulfate (submillion year)—without unrealistically large changes in fluxes (30). This would be appropriate for the excursions in the records of sulfate sulfur and oxygen isotope presented here, which occur within an interval shorter than 40 ky. This would require that P−Tr marine sulfate levels ranging from 18 mM [based on fluid inclusion in marine Late Permian halite (64)] to ≤4 mM [based on rate variation of 100‰ My−1 in sulfur isotopes in the Early Triassic (30)] are valid and possibly a consequence of evaporite deposition during the Permian (30). We tested the influence of the initial seawater sulfate concentration by performing each experiment with a spectrum of sulfate concentrations between 2 mM and 18 mM (Fig. S4).
A plausible scenario to explain the observed sulfur cycle perturbation is volcanism (4), besides scenarios associated with the fluxes of MSR (3). The total volume of sulfur released by volcanic outgassing of Siberian Trap basalts is estimated at 6,300–7,800 Gt (65), with sulfur isotope values close to Holocene values for volcanism. Although crustal contamination of evaporite-derived sulfur has been suggested to enrich evolved sulfate with 34S, this mechanism does not, in general, drive sulfur isotope values heavier than 5‰ (66). The main eruption phase of the Siberian Trap basalts has been constrained to a short magmatic event of <2 My, with the main phase correlated to bed 24 at the Meishan P−Tr section (67). Based on conodont correlation, this is close to the extinction horizon found at the studied sites in NW Iran. Volcanic pulses of less than 10 ky may have contributed to concentrated fluxes of volatiles into the atmosphere, in maybe as much as 25 pulses (65). This would constitute a sulfur flux ∼30 times bigger than background volcanism based on the accumulative volatile release during these pulses. The sulfate oxygen isotope composition of this flux could near 34‰ when considering that equilibrium between water and sulfate is faster in an acidic medium (68). Acidic meteoric water associated with Siberian trap volcanism could be realistic as concomitant land plant destruction has been related to acid rain (34). When constructing a model scenario with these boundary conditions, an excursion as observed in the geochemical data requires an unrealistically large sulfate contribution from volcanism over an extended period (Fig. S4).
Based on a set of analyses with our sulfur cycle coupled to an OC pool, we consider an increased weathering flux as most probable for producing a deviating pattern of sulfate δ34S and δ18O. This scenario does not necessitate unrealistic changes in the isotopic composition or sizes of fluxes. The simultaneous reduction of sulfur burial as pyrite, and the increase of sulfur input via enhanced weathering, can explain the negative sulfate sulfur isotope excursion as observed in CAS, without temporary sequestration of microbial-produced H2S (3); the latter is unlikely, as preextinction large-scale bottom water anoxia in a stagnant ocean has been largely dismissed (13, 33). A low marine sulfate level (2 mM) for the Late Permian to Early Triassic ocean (30) is, however, necessary to explain the observed fluctuations in δ34SCAS.
Supplementary Notes
The R programming environment and the ggplot2 package have been used for graphical plots, statistics, and the box model (69, 70).
Methods
The applied CAS extraction method minimizes the potential alteration of the CAS sulfur or oxygen isotope signal during the analytical procedure by the use of subsequent NaCl leaches (50). The NaCl leaches are precipitated with a BaCl2 solution, and, subsequently, the concentration of the sulfate fraction not structurally substituted within the CaCO3 lattice, leachable non-CAS, is determined gravimetrically. These steps are repeated until no leachable non-CAS is observed. In follow-up, acid digestion of the sample is performed with 25% HCl under restricted atmospheric exchange and in a time window generally not more than 3 h to prevent oxidation of sulfide species. The sulfate of the digested carbonate rock is also precipitated as barium sulfate, and CAS concentrations are gravimetrically determined.
The sulfur isotopic composition of the leachable non-CAS and CAS (Fig. S3) was obtained by an Elemental Analyzer (EA) connected to a ThermoFinnigan Delta Plus mass spectrometer and is reported in the standard (δ34S) notation relative to Vienna Cañon Diablo Troilite (VCDT) in per mil. Precision was confirmed by repeated measurements on the international standards and internal laboratory standards. This yielded an accuracy of 0.09‰ and a reproducibility on average better than 0.40‰ (Table S1). The oxygen isotopic composition of the leachable non-CAS and CAS was measured with a high temperature conversion/EA interfaced to a ThermoFinnigan Delta Plus XL mass spectrometer and is reported in the standard (δ18O) notation relative to Vienna Standard Mean Ocean Water (VSMOW) in per mil. International standards and an in-house standard were used to determine the precision of the procedure. This resulted in an accuracy of 0.40‰ and a reproducibility of ±0.80‰ (Table S2).
Sulfur isotope composition of international and in-house standards
Oxygen isotope composition of international and in-house standards
The CRS content was extracted from the remaining residue after the CAS extraction by the chromium reduction method (51). Under continuous N2 flow, samples were reacted with HCl and chromium powder to produce acidic chromium chloride and were heated to a gentle boil for 1.5 h. Sulfur was trapped as zinc sulfide in a zinc acetate solution and converted to silver sulfide for filtrating and drying. The CRS concentrations of the residues were determined gravimetrically. For sulfur isotopic analyses (δ34S), extracted silver sulfide was weighed into tin capsules withV2O5 added to aid in combustion. Isotope measurements were performed on a Costech Elemental Combustion System 4010 connected to a Thermo Delta V+ Isotope Ratio mass spectrometer with results reported relative to VCDT in per mil. Precision was confirmed by repeated measurements on the international standards [International Atomic Energy Agency (IAEA)-S-1, IAEA-S-2, IAEA-S-3, National Bureau of Standards (NBS)-127] and resulted in an accuracy of 0.24‰ and a reproducibility on average better than 0.40‰.
A box model was constructed, to delineate which fluxes could be responsible for isotopic changes in marine sulfate and thus could explain the stratigraphic CAS isotope record. Sulfate mass (MSO4) and isotopic composition (δSO4) are estimated with the following equations, adapted from ref. 24:
where Qw represents the weathering input flux, δw is the isotopic composition of weathering input, Qv is the volcanic input flux, δv is the isotopic composition of volcanic input, Qm is the hydrothermal input flux, δm is the isotopic composition of the hydrothermal input flux, Fe is the evaporite burial flux, Fp is the pyrite burial flux, and ΔMSR is the average isotopic fractionation associated with MSR. The oxygen and sulfur isotopic compositions of marine sulfate are affected by different processes. Two aspects are mainly responsible for this difference, microbial sulfur cycling and the exchange and cycling through the oceanic crust. Microbial sulfur cycling results in the loss of oxygen during MSR (FMSR), whereas the product (sulfide) is oxidized and mainly via microbial sulfur disproportionation (29) transformed back to sulfate, Qox, both with an effect on the isotopic composition. The surplus of these bacterial processes will be buried as sedimentary sulfide, mainly pyrite, and can be envisioned as
whereas marine sulfate δ34S is only affected by hydrothermal S input in one direction, sulfate δ18O is modified by oxidation of sulfides at these hydrothermal vents (Qh) and acquires a δ18O similar to ocean water (29). On the other hand, anhydrite formation in the oceanic crust is an effective sink for sulfate oxygen (Fs), without fractionation (29) and can be described by
To adjust for these differences imposed upon the marine sulfate δ18O, Eqs. 3 and 4, are substituted into Eq. 1, yielding
where δox stands for the isotopic composition of sulfate originating from sulfur disproportionation. Steady state was achieved by modifying Holocene reservoir and flux size and isotopic composition, to approximate the sulfur cycle of the P−Tr interval (Steady-State Sulfur Cycle Model and Table S5). To exploit the important impact of MSR on OC remineralization (19), an OC pool (MOC) was introduced to the model, by linking Qw linearly to ocean primary productivity (Qpo). OC burial (Fbo) and OC remineralization (FOC) are determined by the size of the OC pool as first-order removal constants (k values); k values are determined by assuming that 99% of OC is remineralized of which part is remineralized by MSR (FMSR; Table S5), and according to the stoichiometry of the two biochemical reactions
Dataset Kuh-e-Ali Bashi 1
Dataset Zal
Boundary conditions for reservoirs and fluxes used in the box model
The size of the marine OC pool is described by
where (i) stands for the initial pool and (t) stands for the pool changing with time. The relation between the OC pool and MSR is described in the following equation, taken from ref. 22:
which considers that the size of the OC pool has a direct impact on the availability of OC for sulfate reduction. Factor α is calculated by solving the equation with the initial values, and y has been taken as 0.3, in accordance with ref. 22. The model was perturbed by varying the weathering flux concomitant with the extinction horizon for a duration of 40 ky. Various simulations were performed, with the weathering flux changing from 2 to 64 times larger than the background value. The Triassic weathering flux is taken to be 7 times higher than the Permian (35). To simulate the effect of a largely iron-limited euxinic ocean, the fraction of pyrite buried was set to 10 times reduced for 20 ky following the extinction horizon, followed by a twofold increase for another 20 ky compared with the preperturbed value.
Acknowledgments
The Aras Free Zone office and Adel Najafzadeh are thanked for permission to sample the locations in the Julfa region (NW Iran). We furthermore thank the following persons for technical assistance; Robert Schreiber, Melanie Rühl, Kathrin Krahn, and Artur Fugmann. Acknowledgments are also due to Robyn Hannigan and Sonny A. Walton for comments on the manuscript. R.H. is thanked for access to the Environmental Analytical Facility at the University of Massachusetts Boston [National Science Foundation (NSF) Award 09-42371 DBI, R.H. and Alan Christian]. This project was funded by the Deutsche Forschungsgemeinschaft (Projects KO2011/8-1 and KO1829/12-2 provided to C.K. and D.K.) and the NSF Graduate Research Fellowship Program (Grant DGE-1356104).
Footnotes
- ↵1To whom correspondence should be addressed. Email: schobbenmartin{at}gmail.com.
Author contributions: M.S., D.K., and C.K. designed research; M.S., A.S., and H.S. performed research; M.S. analyzed data; M.S., A.S., A.G., H.S., D.K., and C.K. wrote the paper; and A.G. performed micropalaeontological analysis.
The authors declare no conflict of interest.
This article is a PNAS Direct Submission. T.W.L. is a guest editor invited by the Editorial Board.
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1503755112/-/DCSupplemental.
References
- ↵.
- Alroy J, et al.
- ↵
- ↵
- ↵
- ↵
- ↵.
- Joachimski MM, et al.
- ↵
- ↵.
- Kump LR,
- Pavlov A,
- Arthur MA
- ↵.
- Algeo T, et al.
- ↵.
- Algeo TJ, et al.
- ↵
- ↵.
- Meyer KM,
- Kump LR,
- Ridgwell A
- ↵.
- Brennecka GA,
- Herrmann AD,
- Algeo TJ,
- Anbar AD
- ↵
- ↵
- ↵.
- Sim MS,
- Bosak T,
- Ono S
- ↵
- ↵.
- Johnston DT, et al.
- ↵
- ↵
- ↵
- ↵.
- Canfield DE,
- Farquhar J
- ↵
- ↵
- ↵
- ↵
- ↵.
- Owens JD, et al.
- ↵.
- Halevy I,
- Peters SE,
- Fischer WW
- ↵.
- Turchyn AV,
- Schrag DP
- ↵
- ↵
- ↵.
- Grice K, et al.
- ↵.
- Winguth C,
- Winguth AME
- ↵.
- Sephton MA, et al.
- ↵.
- Algeo TJ,
- Twitchett RJ
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Payne JL, et al.
- ↵.
- Sun Y, et al.
- ↵
- ↵.
- Wotte T,
- Shields-Zhou GA,
- Strauss H
- ↵
- ↵
- ↵.
- Aghai PM,
- Vachard D,
- Krainer K
- ↵.
- Leda L, et al.
- ↵
- ↵.
- Shen SZ, et al.
- ↵.
- Burgess SD,
- Bowring S,
- Shen SZ
- ↵
- ↵
- ↵.
- Lyons TW,
- Walter LM,
- Gellatly AM,
- Martini AM,
- Blake RE
- ↵.
- Paytan A,
- Kastner M,
- Campbell D,
- Thiemens MH
- ↵.
- Berner RA,
- Canfield DE
- ↵.
- Arthur MA
- ↵
- ↵.
- Black BA,
- Elkins-Tanton LT,
- Rowe MC,
- Peate IU
- ↵
- ↵
- ↵.
- Llyod RM
- ↵.
- RCoreTeam
- ↵.
- Wickham H
- ↵.
- Ghaderi A
Citation Manager Formats
Sign up for Article Alerts
Article Classifications
- Physical Sciences
- Earth, Atmospheric, and Planetary Sciences