An evaporite-based high-resolution sulfur isotope record of Late Permian and Triassic seawater sulfate

Abstract Variations in the sulfur isotope composition of dissolved marine sulfate through time reflect changes in the global sulfur cycle and are intimately related to changes in the carbon and oxygen cycles. A large shift in the sulfur isotope composition of sulfate at the Permian/Triassic boundary has been recognized for long time and a number of studies were carried out to understand the causes and significance of this shift. However, data for the Middle and Late Triassic are very sparse and the stratigraphic evolution of the sulfur isotope composition of seawater is poorly constrained due to the small number of samples analyzed and/or due to the limited stratigraphic intervals studied. Moreover, in the last few years the Triassic timescale has significantly changed due to a wealth of new radiometric and stratigraphic data. In this study we show that for the Late Permian and the Triassic it is possible to obtain a precise reconstruction of the evolution of the sulfur cycle, for parts of it at sub-million year resolution, by analyzing exclusively gypsum and anhydrite deposits. We base our reconstruction on new data from the Middle and Late Triassic evaporites of Northern Switzerland and literature data from evaporites from Germany, Austria, Italy and the Middle East. We propose a revised correlation between the well-dated marine Tethyan sections in northern Italy and the evaporites from Northern Switzerland and from the Germanic Basin calibrated to the newest radiometric absolute age scale. This new correlation allows for a precise dating of the evaporites and constructing a composite sulfur isotope evolution of seawater sulfate from the latest Permian (Lopingian Epoch) to the Norian. We show that a rapid positive shift of approximately 24‰ at the Permian-Triassic boundary can be used to constrain seawater sulfate concentrations in the range of 2–6 mM, thus higher than previous estimates but with less rapid changes. Finally, we discuss two possible evolution scenarios for the Middle Triassic negative shift to values of 15‰ that subsequently remain constant during the Late Triassic.