Lithium has proven a powerful tracer of weathering processes and chemical seawater evolution. Skeletal components of marine calcifying organisms, and in particular brachiopods, present promising archives of Li signatures. However, Li incorporation mechanisms and potential influence from biological processes or environmental conditions require a careful assessment. In order to constrain Li systematics in brachiopod shells, we present Li concentrations and isotope compositions for 11 calcitic brachiopod species collected from six different geographic regions, paralleled with data from culturing experiments where brachiopods were grown under varying environmental conditions and seawater chemistry (pH–pCO2, temperature, Mg/Ca ratio). The recent brachiopod specimens collected across different temperate and polar environments showed broadly consistent δ7Li values ranging from 25.2 to 28.1‰ (with mean δ7Li of 26.9 ± 1.5‰), irrespective of taxonomic rank, indicating that incorporation of Li isotopes into brachiopod shells is not strongly affected by vital effects related to differences among species. This results in Δ7Licalcite–seawater values (per mil difference in 7Li/6Li between brachiopod calcite shell and seawater) from −2.9‰ to −5.8‰ (with mean Δ7Licalcite–seawater value of −3.6‰), which is larger than the Δ7Licalcite–seawater values calculated based on data from planktonic foraminifera (~0‰ to ~−4‰). This range of values is further supported by results from brachiopods cultured experimentally. Under controlled culturing conditions simulating the natural marine environment, the Δ7Licalcite–seawater for Magellania venosa was −2.5‰ and not affected by an increase in temperature from 10 to 16 °C. In contrast, a decrease in Mg/Ca (or Li/Ca) ratio of seawater by addition of CaCl2 as well as elevated pCO2, and hence low-pH conditions, resulted in an increased Δ7Licalcite-seawater up to −4.6‰. Collectively, our results indicate that brachiopods represent valuable archives and provide an envelope for robust Li-based reconstruction of seawater evolution over the Phanerozoic.

In order to allow full comparability with other ocean acidification data sets, the R package seacarb (Gattuso et al, 2021) was used to compute a complete and consistent set of carbonate system variables, as described by Nisumaa et al. (2010). In this dataset the original values were archived in addition with the recalculated parameters (see related PI). The date of carbonate chemistry calculation by seacarb is 2022-04-19.