Search

Groundwater-derived solute fluxes to the ocean have long been assumed static and subordinate to riverine fluxes, if not neglected entirely, in marine isotope budgets. Here we present concentration and isotope data for Li, Mg, Ca, Sr, and Ba in coastal groundwaters to constrain the importance of groundwater discharge in mediating the magnitude and isotopic composition of terrestrially derived solute fluxes to the ocean. Data were extrapolated globally using three independent volumetric estimates of groundwater discharge to coastal waters, from which we estimate that groundwater-derived solute fluxes represent, at a minimum, 5% of riverine fluxes for Li, Mg, Ca, Sr, and Ba. The isotopic compositions of the groundwater-derived Mg, Ca, and Sr fluxes are distinct from global riverine averages, while Li and Ba fluxes are isotopically indistinguishable from rivers. These differences reflect a strong dependence on coastal lithology that should be considered a priority for parameterization in Earth-system models.

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.

The geological calcium cycle is linked to the geological carbon cycle through the weathering and burial of carbonate rocks. As a result, calcium (Ca) isotope ratios (44Ca/40Ca, expressed as δ44/40Ca) can help to constrain ancient carbon cycle dynamics if Ca cycle behavior can be reconstructed. However, the δ44/40Ca of carbonate rocks is influenced not only by the δ44/40Ca of seawater but also by diagenetic processes and fractionation associated with carbonate precipitation. In this study, we investigate the dominant controls on carbonate δ44/40Ca in Upper Permian to Middle Triassic limestones (ca. 253 to 244 Ma) from south China and Turkey. This time interval is ideal for assessing controls on Ca isotope ratios in carbonate rocks because fluctuations in seawater δ44/40Ca may be expected based on several large carbon isotope (δ13C) excursions ranging from − 2 to + 8‰. Parallel negative δ13C and δ44/40Ca excursions were previously identified across the end-Permian extinction horizon. Here, we find a second negative excursion in δ44/40Ca of ~ 0.2‰ within Lower Triassic strata in both south China and Turkey; however, this excursion is not synchronous between regions and thus cannot be interpreted to reflect secular change in the δ44/40Ca of global seawater. Additionally, δ44/40Ca values from Turkey are consistently 0.3‰ lower than contemporaneous samples from south China, providing further support for local or regional influences. By measuring δ44/40Ca and Sr concentrations ([Sr]) in two stratigraphic sections located at opposite margins of the Paleo-Tethys Ocean, we can determine whether the data represent global conditions (e.g., secular variations in the δ44/40Ca of seawater) versus local controls (e.g., original mineralogy or diagenetic alteration). The [Sr] and δ44/40Ca data from this study are best described statistically by a log-linear correlation that also exists in many previously published datasets of various geological ages. Using a model of early marine diagenetic water-rock interaction, we illustrate that this general correlation can be explained by the chemical evolution of bulk carbonate sediment samples with different initial mineralogical compositions that subsequently underwent recrystallization. Although early diagenetic resetting and carbonate mineralogy strongly influence the carbonate δ44/40Ca values, the relationship between [Sr] and δ44/40Ca holds potential for reconstructing first-order secular changes in seawater δ44/40Ca composition.