Explaining the ~5-million-year delay in marine biotic recovery following the latest Permian mass extinction, the largest biotic crisis of the Phanerozoic, is a fundamental challenge for both geological and biological sciences. Ocean redox perturbations may have played a critical role in this delayed recovery. However, the lack of quantitative constraints on the details of Early Triassic oceanic anoxia (for example, time, duration, and extent) leaves the links between oceanic conditions and the delayed biotic recovery ambiguous. We report high-resolution U-isotope (δ²³⁸U) data from carbonates of the uppermost Permian to lowermost Middle Triassic Zal section (Iran) to characterize the timing and global extent of ocean redox variation during the Early Triassic. Our δ²³⁸U record reveals multiple negative shifts during the Early Triassic. Isotope mass-balance modeling suggests that the global area of anoxic seafloor expanded substantially in the Early Triassic, peaking during the latest Permian to mid-Griesbachian, the late Griesbachian to mid-Dienerian, the Smithian-Spathian transition, and the Early/Middle Triassic transition. Comparisons of the U-, C-, and Sr-isotope records with a modeled seawater PO₄^3- concentration curve for the Early Triassic suggest that elevated marine productivity and enhanced oceanic stratification were likely the immediate causes of expanded oceanic anoxia. The patterns of redox variation documented by the U-isotope record show a good first-order correspondence to peaks in ammonoid extinctions during the Early Triassic. Our results indicate that multiple oscillations in oceanic anoxia modulated the recovery of marine ecosystems following the latest Permian mass extinction.

A negative shift in the calcium isotopic composition of marine carbonate rocks spanning the end-Permian extinction horizon in South China has been used to argue for an ocean acidification event coincident with mass extinction. This interpretation has proven controversial, both because the excursion has not been demonstrated across multiple, widely separated localities, and because modeling results of coupled carbon and calcium isotope records illustrate that calcium cycle imbalances alone cannot account for the full magnitude of the isotope excursion. Here, we further test potential controls on the Permian-Triassic calcium isotope record by measuring calcium isotope ratios from shallow-marine carbonate successions spanning the Permian-Triassic boundary in Turkey, Italy, and Oman. All measured sections display negative shifts in δ44/40Ca of up to 0.6‰. Consistency in the direction, magnitude, and timing of the calcium isotope excursion across these widely separated localities implies a primary and global δ44/40Ca signature. Based on the results of a coupled box model of the geological carbon and calcium cycles, we interpret the excursion to reflect a series of consequences arising from volcanic CO2 release, including a temporary decrease in seawater δ44/40Ca due to short-lived ocean acidification and a more protracted increase in calcium isotope fractionation associated with a shift toward more primary aragonite in the sediment and, potentially, subsequently elevated carbonate saturation states caused by the persistence of elevated CO2 delivery from volcanism. Locally, changing balances between aragonite and calcite production are sufficient to account for the calcium isotope excursions, but this effect alone does not explain the globally observed negative excursion in the δ13C values of carbonate sediments and organic matter as well. Only a carbon release event and related geochemical consequences are consistent both with calcium and carbon isotope data. The carbon release scenario can also account for oxygen isotope evidence for dramatic and protracted global warming as well as paleontological evidence for the preferential extinction of marine animals most susceptible to acidification, warming, and anoxia.