Ca Isotope Fractionation Research Articles

Understanding Water–Rock Interaction in Crystalline Shield Fluids Using Calcium Isotopes

Calcium isotopes provide a potentially robust tool for understanding the evolution of crystalline shield fluids, but previous applications have focused on near surface groundwaters. The tendency of Ca isotopes to be affected by mass-dependent fractionation during processes such as water–rock interaction provides a powerful tool for studying the evolution and origin of groundwaters hosted in crystalline rocks. We report Ca isotope ratios (δ44/40Ca) of deep fluids (>300m) from across the Canadian Shield and integrate these with Sr and Br isotope values to understand long-term fluid–rock interactions in crystalline shield environments. Ca isotope ratios have a wide range of values from 0.07‰ to 0.86‰. At individual sites, δ44/40Ca values are variable whereas 87Sr/86Sr ratios are relatively constant. Sr isotope ratios (87Sr/86Sr) have a negative relationship with the Ca vs. Na content of the fluid indicating different host-minerals contributing Sr to the fluid. Ca isotope fractionation was caused by metamorphic reactions and by the growth and dissolution of Ca-rich fracture-filling minerals. The δ44/40Ca signatures of these processes are overprinted by radiogenic ingrowth of 40Ca by decay of 40K, which is expected to affect older and more K-rich rocks. At one site, δ44/40Ca and δ81Br variability reflects gas-generating reactions in the fluid and/or water–rock interaction processes. These new results demonstrate the strength of combining multiple isotope analysis to elucidate the sources of groundwater salinity and decipher the complex long-term processes that occur in crystalline shield environments.

Ca–Zn isotope fractionation during fluid–rock interaction in subduction zones and its implication for deep carbon cycle

Subducting slabs transport Earth's surface carbon into the deep mantle, with a portion of this carbon potentially recycled through volcanic eruptions at mid-ocean ridges or ocean islands. Non-traditional stable isotopes such as calcium (Ca) and zinc (Zn) are increasingly used to trace this deep carbon cycle. However, the potential for Ca and Zn isotope fractionation in the subduction zone, critical for the successful application of these tracers, remains insufficiently explored. In this study, we investigate high-pressure carbonated metamafic rocks and enclosed vein systems from a Paleo-oceanic subduction complex (SW Tianshan) to explore the possible fractionation of Ca and Zn isotope during subduction-related metamorphism and fluid–rock interaction. The δ44/40Ca (0.77 ‰–0.90 ‰) and δ66Zn isotope compositions (0.20 ‰–0.39 ‰) of greenschist-, blueschist- and eclogite-facies metabasites are consistent with their protoliths (OIBs and MORBs). We observe no systematic fractionation of either Ca and Zn isotopes across prograde metamorphism and associated dehydration. However, the precipitation of carbonate in surrounding rock along fluid pathways leads to an increase in the δ44/40Ca of residual fluids (>0.90 ‰). More significantly, we find detectable Zn isotope fractionation (ΔA-H = ẟ66Znaltered - ẟ66Znhost of 0.06–0.10 ‰) during fluid-mediated metasomatism of some blueschists/eclogites. Light Zn isotopes are preferentially released into the infiltrating fluid, resulting in an elevated δ66Zn isotope composition in the immediate carbonated eclogites along the vein. The modified slab fluids during migration are likely charactered by heavy δ44/40Ca and light δ66Zn compositions, which may play a crucial role in determining the low δ66Zn values found in some arc lavas. However, the MORB-like δ44/40Ca compositions observed in global arc magmatic rocks suggest that the slab-released Ca (possibly along with some carbon) is probably sequestered and stored along the migration path of the fluid, thereby hindering the slab–arc carbon recycling. Our study emphasizes that intraslab fluid–rock interaction results in similar δ44/40Ca and δ66Zn signature in carbonate-bearing eclogites compared to sedimentary carbonates. Consequently, deeply subducted carbonate-bearing eclogites –not necessarily sedimentary carbonates– might carry light δ44/40Ca and heavy δ66Zn signatures into the deep mantle, potentially influencing the corresponding Ca-Zn isotope compositions of some mantle-derived basalts.

Experimental determination of Si, Mg, and Ca isotope fractionation during enstatite melt evaporation

Abstract Evaporation of silicate materials from Earth or its precursors may be important in shaping their primordial compositions represented by undifferentiated meteorites, e.g., enstatite chondrites; however, the conditions under which evaporation occurs and the extent of evaporation-induced elemental and isotope fractionation remain uncertain. Here, we experimentally determine the volatility and isotope fractionation of Si, Mg, Ca, Nb, and Ta during enstatite melt evaporation at 2423–2623 K using a high-temperature conical nozzle levitator. Homogenous glasses are recovered after experiments; then we use EPMA and LA-ICP-MS to measure the elemental compositions, MC-ICP-MS to measure the Si and Mg isotopes, and TIMS to measure the Ca isotopes. Our results show that the evaporation rates of Si are larger than Mg, and the mean vapor/melt isotope fractionation factors (α = Rvapor/Rmelt; R = isotope ratio) are 0.99585 ± 0.00002 for 29Si/28Si and 0.98942 ± 0.00130 for 25Mg/24Mg. However, neither evaporative loss of Ca, Nb, and Ta nor Ca isotope fractionation was observed within analytical uncertainty. In conjunction with previous studies, we find that in an evaporation experiment the saturation degree (partial vapor pressure/equilibrium vapor pressure) of Si (SSi) is larger than SMg when Si is more volatile than Mg, and vice versa. If the Mg/Ca and Si/Ca ratios and isotopes in the bulk silicate Earth are attributed to the evaporation of enstatite chondrite-like precursors, evaporation temperatures >5000 K and SSi < SMg are required.

Calcium isotope fractionation during weathering of argillaceous carbonates in a humid climate

The continental weathering is a key process that controls calcium (Ca) transportation from the continental crust to the waters. To elucidate the behavior of Ca isotopes during carbonate weathering, the concentrations and δ44/40Ca (relative to NIST SRM 915a) of bulk saprolites, exchangeable, acid-leachable and residual phases of a weathering profile developed on the marine carbonates, Guangdong province, South China, were investigated. Upwards the profile, δ44/40Ca values of the bulk saprolites systematically decrease from 0.77 ± 0.12 ‰ to −0.44 ± 0.12 ‰, suggesting that significant Ca isotope fractionation occurred during chemical weathering. The exchangeable fractions have δ44/40Ca values higher than those of the bulk saprolites with Δ44/40Caexchangeable-saprolite varying from −0.01 ‰ to 0.73 ‰, suggesting that heavy isotopes are preferentially adsorbed onto the clays. The acid-leachable phases display a relatively narrow δ44/40Ca range from 0.52 ‰ to 0.74 ‰ with Ca fractions varying from 7.4 % to 100.3 %, potentially indicating that limited Ca isotopic fractionation occurs during the dissolution of primary carbonates. The residual Ca pool is strongly fractionated with δ44/40Ca ranging from 0.64 ± 0.08 ‰ to −0.98 ± 0.02 ‰, systematically lower than their bulk saprolites, perhaps indicating light Ca isotopes are preferentially incorporated into the clay lattices. Altogether, it seems that the Ca isotopic fractionation directions are opposite between clay structural incorporation and adsorption. Our study provides important insight of Ca behavior and Ca isotopic fractionation during chemical weathering, which is critical to shape Ca isotopic compositions of the upper continental crust and trace the global biogeochemical cycle of Ca.

The Genesis of A‐Type Granites in Lower Yangtze River Belt: Evidence From Calcium Isotopes

AbstractA‐type granite is a favorable rock for the study of crustal evolution, crust‐mantle interaction and metallogenesis. However, the origin of A‐type granite remains controversial. In this study, we report high‐precision Ca isotopic compositions of a co‐genetic suite of A‐type granites from the Lower Yangtze River Belt (LYRB) in eastern China. Our results show that the δ44/40Ca of the A‐type granites ranges from 0.55 ± 0.11‰ to 1.44 ± 0.06‰, and is negatively correlated with CaO, Sr and Ba contents and Eu/Eu*, indicating that the Ca isotope fractionation of A‐type granites is dominated by plagioclase and potassium feldspar. There is kinetic fractionation of Ca isotopes during A‐type granite magma evolution. Sample BSL‐7, the least evolved sample, revealed that it had a low δ44/40Ca of 0.59‰, which is significantly lower than that of the upper continental crust (δ44/40Ca = 0.71–0.72‰) and the bulk silicate earth (δ44/40Ca = 0.94‰). The low δ44/40Ca can be best explained by partial melting of the enriched mantle metasomatized by subduction material. Combined with the geodynamic background of the LYRB, we propose that the formation of A‐type granites in the LYRB originates from partial melting of the enriched mantle, triggered by early Cretaceous ridge subduction of the Pacific and Izanagi plates.

Insights on the origin of oldhamite in enstatite meteorites from Ca stable isotopes

In order to understand the origin of oldhamite (CaS) in enstatite meteorites, we report Ca isotopic compositions (δ44/40Ca) of oldhamite (obtained from water leachate of bulk chondrites and aubrites and mineral separates from the Norton County aubrite) and silicate minerals from different types of enstatite chondrites and aubrite. The δ44/40Ca of the bulk enstatite chondrites range from 1.05 ‰ to 1.24 ‰, with an average of 1.13 ± 0.12 ‰, higher than that of the estimate of the bulk silicate earth (∼0.94 ‰). Major and trace element analyses show that the water leachates of enstatite chondrites are mainly composed of oldhamite, and they take over 20.6–68.5 % Ca of the bulk meteorite Ca budget. The Ca isotope fractionation between oldhamite and residual silicate (Δ44/40Caoldhamite-silicate) for the studied enstatite chondrites is minimum (−0.44 ‰) for Abee (impact-melt breccia) and maximum (+0.16) for St.Marks (EH5). The Ca isotope fractionation between oldhamite (individual mineral grains and leachate) and silicates in Norton County varies from −0.47 ‰ to −0.31 ‰ with an average of −0.41 ‰. These Δ44/40Caoldhamite-silicate correlate well with previous theoretical calculation and suggests that the oldhamites in Norton County are in isotopic equilibrium with co-existing silicates, and therefore were formed during magmatic processes. However, in enstatite chondrites, the large variation on Δ44/40Caoldhamite-silicate and its negative correlation with metamorphic temperature reflects the redistribution and equilibration of Ca isotopes during metamorphism. The variable Δ44/40Caoldhamite-silicate found in unequilibrated chondrites reflect kinetic Ca isotope fractionation between oldhamite and nebular gas and therefore is evidence for the formation of oldhamite by condensation in the solar nebula.

Ca isotopic gradients in epeiric marine carbonates: Diagenetic origins of and significance for Ca cycle reconstructions

Marine carbonates record geochemical information on the state of Earth’s surface environment in the geological past, but recrystallization during burial can make the records difficult to interpret. Here, we report shore to basin gradients in δ44Ca values in matrix limestone and burrow dolomite across a sequence of Late Ordovician (Katian Stage) burrow-mottled carbonates in the Williston Basin, an intracratonic epeiric marine basin in North America. Both gradients are interpreted to have diagenetic origins. The gradient in limestone δ44Ca values developed in a bathymetrically controlled sediment- to fluid-buffered diagenetic system operating between the center and the edges of the basin. Two dolomitizing events formed the gradient in dolomite δ44Ca values. The sediment-buffered end of the limestone gradient preserves the δ44Ca value of the original carbonate mud (∼–1.28 ± 0.10 ‰). The fluid-buffered end of the gradient conserves the δ44Ca value of contemporaneous seawater (–0.47 ± 0.10 ‰), because: (1) pore water exchange drivers operated with greater variety and intensity in shallow water settings, (2) shallow nearshore deposits exchanged more Ca with seawater during recrystallization to diagenetic calcite than deeper offshore deposits, and (3) diagenetic calcite formed with negligible fractionation under equilibrium conditions. The gradient conserves the local value of Δsed (–0.81 ± 0.20 ‰) despite its diagenetic origin, and the magnitude of Δsed is consistent with primary calcite precipitation in a ‘calcite sea’. Field assessment of the equilibrium Ca isotope fractionation factor between calcite and dolomite of –0.41 ± 0.17 ‰ at ambient temperature is applied to correct basin-centered dolomite for fractionation, yielding –1.62 ± 0.27 ‰, which is very different from the –0.25 ± 0.07 ‰ value of Ca in evaporatively concentrated seawater from which the dolomite originally formed. Applying the same fractionation factor to correct literature data on early diagenetic dolomite in the Great Basin, North America, yields –0.69 ± 0.22 ‰ for the δ44Ca value of seawater, which is ∼0.2 ‰ lower than ‘normal marine’ waters in stratigraphically equivalent deposits in the Willison Basin, and 0.44 ‰ lower than evaporatively concentrated seawater during a brief period of basin restriction and anhydrite deposition. The nearly 0.5 ‰ variation in the δ44Ca value of Katian seawater falls within the range of scatter in the brachiopod reconstruction of seawater Ca isotopes. Here, we propose that the true ocean record of secular change lies towards the bottom of the range of scatter, and that the higher values in the record reflect local Ca cycling and circulation restriction in epeiric seas. Epeiric seas cultivate the fluid-buffered diagenetic conditions required to promote Ca isotope exchange between reactive carbonate sediments and seawater through carbonate mineral dissolution/recrystallization, crystal coarsening by way of Ostwald ripening, and replacive dolomitization. It is not the apportioning of kinetically fractionated aragonite and calcite in the carbonate deposition flux that drove secular changes in global value of Δsed between ‘calcite’ and ‘aragonite’ seas, but rather the replacement of kinetic fractionations with equilibrium fractionations during fluid-buffered diagenesis with seawater.

Calcium isotope fractionation by intracellular amorphous calcium carbonate (ACC) forming cyanobacteria.

The formation of intracellular amorphous calcium carbonate (ACC) by various cyanobacteria is a widespread biomineralization process, yet its mechanism and importance in past and modern environments remain to be fully comprehended. This study explores whether calcium (Ca) isotope fractionation, linked to ACC-forming cyanobacteria, can serve as a reliable tracer for detecting these microorganisms in modern and ancient settings. Accordingly, we measured stable Ca isotope fractionation during Ca uptake by the intracellular ACC-forming cyanobacterium Cyanothece sp. PCC 7425. Our results show that Cyanothece sp. PCC 7425 cells are enriched in lighter Ca isotopes relative to the solution. This finding is consistent with the kinetic isotope effects observed in the Ca isotope fractionation during biogenic carbonate formation by marine calcifying organisms. The Ca isotope composition of Cyanothece sp. PCC 7425 was accurately modeled using a Rayleigh fractionation model, resulting in a Ca isotope fractionation factor (Δ44Ca) equal to -0.72 ± 0.05‰. Numerical modeling suggests that Ca uptake by these cyanobacteria is primarily unidirectional, with minimal back reaction observed over the duration of the experiment. Finally, we compared our Δ44Ca values with those of other biotic and abiotic carbonates, revealing similarities with organisms that form biogenic calcite. These similarities raise questions about the effectiveness of using the Ca isotope fractionation factor as a univocal tracer of ACC-forming cyanobacteria in the environment. We propose that the use of Δ44Ca in combination with other proposed tracers of ACC-forming cyanobacteria such as Ba and Sr isotope fractionation factors and/or elevated Ba/Ca and Sr/Ca ratios may provide a more reliable approach.

Calcium isotope compositions of subduction-related leucite-bearing rocks: Implications for the calcium isotope heterogeneity of the mantle and carbonate recycling in convergent margins

Large Ca isotope variations recorded in mantle xenoliths and mantle-derived igneous rocks have been widely used to decode the origin of mantle heterogeneity, especially to trace recycled sedimentary carbonates. However, it is still unclear whether and how the subduction of sedimentary carbonate causes significant Ca isotope fractionation in the mantle in correspondence of convergent margins. In this study, we present Ca isotope data for sixteen silica-undersaturated ultrapotassic leucite-bearing rocks and their clinopyroxene (Cpx) phenocrysts from the “Colli Albani” volcano (Italy), and for six carbonate-bearing sediments from nearby Apennine Chain. These subduction-related volcanic rocks have been widely accepted as deriving from a mantle source that involves recycled carbonated sediments indicated by their enriched incompatible trace elements and radiogenic isotopes, and olivine geochemistry. The δ44/40Ca values (normalized to NIST SRM 915a) of leucite-bearing rocks are uniform (from 0.64 to 0.79 ‰, with one exception of 0.86 ‰) with an average of 0.72 ± 0.03 ‰, 2SE, N = 16, which is 0.12 ± 0.03 ‰ (2SE) and 0.22 ± 0.03 ‰ (2SE) lower than MORBs (0.84 ± 0.02 ‰, 2SE, N = 25) and BSE (0.94 ± 0.01 ‰, 2SE, N = 14) respectively. In contrast, the carbonate-bearing sediments show variable and notably low δ44/40Ca values, ranging from 0.44 to 0.77 ‰. The δ44/40Ca of Cpx (0.70 to 0.76 ‰) is indistinguishable from their host rocks (0.73 to 0.76 ‰). In addition, the lack of correlations between δ44/40Ca and SiO2, CaO, Ni and P2O5, and the Ca isotopic simulation result indicate that Ca isotope fractionation caused by fractional crystallization is insignificant. Given that the Ca isotope fractionation is limited during the melting of spinel wehrlite-like mantle source (within ∼0.04 to 0.11 ‰), the slightly low δ44/40Ca values of the Italian leucite-bearing rocks require 10 % to 20 % addition of carbonated sediments in their mantle source. This indicates that recycled carbonated sediments can cause slight δ44/40Ca heterogeneity (∼0.1 to 0.2) of the mantle wedge in local/regional scales.

Calcium isotopes track volatile components in the mantle sources of alkaline rocks and associated carbonatites

The volatile components CO2 and H2O induce mantle melting and thus exert major controls on mantle heterogeneity. Primitive intraplate alkaline magmatic rocks are the closest analogues for incipient mantle melts and provide the most direct method to assess such mantle heterogeneity. Given the considerable Ca isotope differences among carbonate, clinopyroxene, garnet, and orthopyroxene in the mantle (up to 1 ‰ for δ44/40Ca), δ44/40Ca of alkaline rocks is a promising tracer of lithological heterogeneity. We present stable Ca isotope data for ca. 1.4 Ga lamproites, 590–555 Ma ultramafic lamprophyres and carbonatites, and 142 Ma nephelinites from Aillik Bay in Labrador, eastern Canada. These primitive alkaline rock suites are the products of three stages of magmatism that accompanied lithospheric thinning and rifting of the North Atlantic craton. The three discrete magmatic events formed by melting of different lithologies in a metasomatized lithospheric mantle column at various depths: (1) MARID-like components (mica-amphibole-rutile-ilmenite-diopside) in the source of the lamproites; (2) phlogopite-carbonate veins were an additional source component for ultramafic lamprophyres during the second event; and (3) wehrlites at shallower depths were an important source component for nephelinites during the final event.The Mesoproterozoic lamproites show lower δ44/40Ca values (0.58 to 0.66 ‰) than MORBs (0.84 ± 0.03 ‰, 2se). This cannot be explained by fractional crystallization or melting of the clinopyroxene-dominated source but can be attributed to a source enriched in the alkali amphibole K-richterite, which has characteristically low δ44/40Ca. The δ44/40Ca values of the ultramafic lamprophyre suite during the second rifting stage are remarkably uniform, with overlapping ranges for primary carbonated silicate melts (aillikite: 0.67 to 0.75 ‰), conjugate carbonatitic liquids (0.71 to 0.82 ‰) and silicate-dominated damtjernite liquid (primary damtjernite: 0.68 to 0.72 ‰). This suggests negligible Ca isotope fractionation during liquid immiscibility of carbonate-bearing magmas. Combined with previously reported δ44/40Ca values for carbonatites and kimberlites, our data suggest that carbonated silicate melts in Earth's mantle have δ44/40Ca compositions resolvably lower than those for MORBs (0.74 ± 0.02 ‰ versus 0.84 ± 0.03 ‰, 2se). The δ44/40Ca values of the Cretaceous nephelinites (0.72 to 0.78 ‰) are homogenous and similar to those of the 590–555 Ma ultramafic lamprophyres, suggesting that the wehrlitic source component for the nephelinites formed by mantle metasomatism during interaction with rising aillikite magmas during the second rifting stage. Our results highlight that both K-richterite and carbonate components in mantle sources can result in the systematically low δ44/40Ca values of alkaline magmas, which may explain previously reported low δ44/40Ca values of alkaline rocks and some carbonatites. Our study indicates that Ca isotopes are a robust tracer of lithological variation caused by volatiles in the Earth's upper mantle.

Calcium isotopic fractionation during aragonite and high-Mg calcite precipitation at methane seeps

The formation of authigenic carbonate in marine environments represents a process that buffers ocean chemistry and atmospheric carbon dioxide levels. The isotopic composition of calcium (δ44/40Ca) in authigenic carbonate can be used to investigate the calcium (Ca) cycle, seawater chemistry, and diagenesis of carbonate rocks, archiving key information on the evolution of the Earth's surface environments. Laboratory experiments have shown that Ca isotopic fractionation during carbonate precipitation – the basis of the δ44/40Ca proxy – is mainly determined by both mineral phase and precipitation rate. However, the precipitation rate of submarine authigenic carbonate is much lower than rates in laboratory experiments. This study investigated the Ca isotopic composition of pore water and carbonate nodules collected from two gravity piston cores and one push core from the Haima seeps in the northern South China Sea. Methane fluxes calculated for the three cores range from 12 to 134 µmol cm−2 yr−1. The core with intermediate methane flux showed evidence of aragonite precipitation, while high-Mg calcite was the dominant mineral phase in the other two cores as indicated by Mg/Ca and Sr/Ca ratios of pore water. Numerical modeling of pore water geochemistry revealed that: (1) precipitation rates of carbonate per unit reactive surface area range from 0.05 to 21 µmol m−2 h−1, which is one to three orders of magnitude lower than rates in precipitation experiments; (2) Ca isotopic fractionation for aragonite precipitation is –1.8 ± 0.3 ‰, close to the reported equilibrium isotopic fractionation in precipitation experiments; and (3) Ca isotopic fractionation for high-Mg calcite precipitation is similar for the two cores (−0.9 ± 0.3 ‰ and −0.8 ± 0.3 ‰), although the modeled precipitation rates differed by two orders of magnitude. These results indicate that mineralogy rather than precipitation rate is the primary control on Ca isotopic fractionation during carbonate precipitation at seeps, resulting in lower δ44/40Ca values of aragonite (0.91±0.06 ‰, N = 4) than of high-Mg calcite (1.41±0.04 ‰, N = 9) reported relative to NIST SRM915a. This study documents distinguishable Ca isotopic fractionation during the formation of methane-derived high-Mg calcite and aragonite, providing fundamental parameters for further exploration of the calcium isotopic composition of marine authigenic carbonates in the geological record.

Light calcium isotope anomaly observed in continental basaltic lavas: A mixed signal of recycled carbonate and fractionation during melting

We present high-precision Ca isotopic data for 21 Cenozoic basaltic lavas from eastern China to evaluate the potential of Ca isotopes as a tool to trace the recycling of surface carbonates. These samples were derived from a variety of peridotite to pyroxenite mantle sources hybridized by different slab-derived agents. The δ44/42Ca of eastern China Cenozoic basaltic lavas ranges from 0.24 ± 0.03‰ (2SE, if not specified) to 0.39 ± 0.03‰, with an average of 0.31 ± 0.09‰ (2SD, N = 21), systematically lower than mid-ocean ridge basalts (MORB, 0.39 ± 0.07‰, 2SD, N = 30). No convincing correlation has been observed with δ26Mg and δ66Zn, indicating that metasomatism by recycled carbonates may not always yield Ca isotopically distinct mantle sources. δ44/42Ca shows correlations with parameters reflecting degree of partial melting (e.g., K2O, Nb, LREE, U, Th, Rb, and Sr) and residual garnet signature (i.e., (Dy/Yb)N). Nevertheless, modeling results suggest that Ca isotope fractionation during mantle and/or slab melting with the presence of garnet can partially account for the observed variable δ44/42Ca. Recycled carbonates with low δ44/42Ca were likely involved in the origins of the lightest lavas. Particularly, the nephelinitic lightest lavas also have the lowest δ26Mg, Ti/Ti⁎ and Hf/Hf⁎ among all the studied samples, reflecting a role of carbonatitic metasomatism in their sources. This study illustrates that a strongly light Ca isotope anomaly (e.g., δ44/42Ca ∼ 0.24‰) observed in basaltic lavas is likely associated with the recycling of carbonates into the mantle.

Calcium isotopes reveal shelf acidification on southern Neotethyan margin during the Smithian-Spathian boundary cooling event

The Smithian-Spathian transition of the late Early Triassic was a critical period of environmental and biological upheavals, following the end-Permian mass extinction. Changes in carbonate deposition during this period have been attributed to intensified upwelling along shelf margins, but relevant studies are scarce. Here, we present calcium isotopes of bulk marine carbonate (δ44/40Cacarb) from a Smithian–Spathian boundary (SSB) succession (Guryul Ravine section, Kashmir) on the southern margin of the Neo-Tethys. Our smoothed δ44/40Cacarb curve reveals a ∼ 0.2‰ negative shift (from ∼ − 1.1‰ to ∼ − 1.3‰) across the SSB, concurrent with a ∼ +10‰ shift in δ13Ccarb. While increased Ca isotopic fractionation could play a role, we specifically examine potential impacts due to changes in marine Ca fluxes. Using a Ca-cycle mass balance model, we explore scenarios of decreased carbonate burial flux (Fcarb), decreased riverine flux (Friv), and a combination of these processes. The modeling suggest that a pulse decrease in Fcarb by 40% over ∼0.06 Myr match the negative shift in δ44/40Cacarb at Guryul Ravine. We infer that this decrease was likely related to intensified upwelling of acidic deep seawater due to invigorated global-oceanic circulation during the SSB cooling event. We suggest that the regionally diverse excursions in δ44/40Cacarb in the Tethyan region could be attributed to spatially varied upwellings in the shelf margin. The upwelling of acidic and anoxic deep seawater may have driven the second-order extinction of ammonoids and conodonts at the beginning of the SSB cooling event.

Subducted carbonates not required: Deep mantle melting explains stable Ca isotopes in kimberlite magmas

Ca isotope geochemistry has great potential for improving our understanding of magmatic systems and for tracing the deep Earth carbon cycle. There are still many open questions, however, regarding the proper application of this relatively novel proxy to the study of mantle-derived magmas, including (i) the possible effects of pressure on mineral-melt Ca isotope fractionation factors, and (ii) the potential for Ca isotopes to be used as tracers of recycled marine carbonates in mantle-derived magmas. Kimberlites are mantle-derived melts that are highly enriched in CO2 and are the deepest-sourced magmas (>200 km depth) known to erupt at Earth’s surface, providing an excellent opportunity to explore these questions. We present Ca isotope data combined with detailed petrographic observations, bulk-carbonate C-O isotope data, and bulk-rock major element analyses, for a suite of 23 well-characterized kimberlite samples from their type-locality (Kimberley, South Africa). These kimberlites have abundant previous evidence for recycled surface materials in their mantle source, including low S isotope and moderately radiogenic Sr isotope compositions, yet display only limited variations in their Ca isotope compositions (δ44CaBSE of −0.08‰ to −0.27‰), with an average of −0.17 ± 0.02‰ (2SE, n = 21). This composition is indistinguishable from average carbonatites [−0.19 ± 0.03‰ (2SE, n = 106)] and OIB from recent studies [−0.16 ± 0.01‰ (2SE, n = 41)], and slightly lower than average MORB [−0.11 ± 0.02‰ (2SE, n = 31)]. Although our samples display a wide range of emplacement styles, alteration conditions, extents of magmatic differentiation, and degrees of mantle-cargo entrainment (i.e., xenocryst accumulation), we find no correlations between Ca isotopes and any of these factors. Instead, we find that low-degree partial melting of the likely kimberlite source lithology (i.e., carbon-bearing garnet lherzolite) yields modelled melt δ44CaBSE values ranging between −0.12‰ and −0.16‰ (at 1400–1500 °C), in agreement with the measured Ca isotope compositions of the Kimberley kimberlites. This observation, and the lack of heavy carbon isotope signatures in the examined samples, indicates that kimberlites do not require subducted carbonates in their mantle sources, despite their very high CO2 contents. Although several recent studies have suggested that equilibrium mineral-melt Ca isotope fractionation factors (e.g., 1000lnαgrt-melt) could be significantly different at higher pressures (i.e., due to pressure-induced changes in CaO bond lengths and coordinations), our models successfully reproduce the kimberlite data using pressure-independent predictions for mineral-melt fractionations. It remains possible, however, that differences in isotopic fractionation due to the peculiar composition of kimberlite melts (e.g., high CO2, low SiO2) are effectively cancelled out by competing pressure effects, and future work independently targeting these factors will be especially important for our understanding of Ca isotope fractionation in mantle-derived melts and the Ca isotope systematics of Earth’s mantle.

Granite petrogenesis and the δ44Ca of continental crust

Stable Ca isotopes are an increasingly useful tool for understanding the sources and processes leading to the formation of magmatic rocks, yet Ca isotope fractionation during genesis of silicic continental crust is still poorly understood. Here, we present Ca, Sr, and Nd isotope, as well as major- and trace-element whole-rock geochemical data for A-, I-, and S-type granites (n = 30) from Australia/Tasmania, Canada, and France (δ44CaBSE of -0.6‰ to +0.2‰) and compare them to phase-equilibrium models for partial-melting (pelite, greywacke, MORB, enriched Archean tholeiite) and crystallization (hydrous arc basalt, A-type granite) that incorporate novel ab-initio predictions for Ca isotope fractionation in epidote and K-feldspar. The ab-initio calculations predict that epidote has similar δ44Ca to anorthite and that K-feldspar is the isotopically lightest known silicate mineral at equilibrium (Δ44Cakspar-melt of -0.4‰ at 1000 K). Our phase-equilibrium model results suggest that δ44Ca variations in all three granite types can be fully explained through magmatic processes, without necessarily requiring addition of isotopically exotic components (e.g., carbonate sediments). Heavy Ca isotope enrichments in A-type granites from the Lachlan Fold Belt, however, require isotopic disequilibrium between plagioclase and melt, which we use to constrain average plagioclase growth rates in these systems. This also serves to illustrate that whole-rock Ca isotope measurements can be used to estimate crystal growth rates, even in the absence of analyzable phenocrysts. In general, low Ca diffusivities and strong isotopic diffusivity ratios (D44/D40) in low-H2O granitic magmas should lead to resolvable isotopic disequilibrium effects in plagioclase, even at relatively slow growth rates (e.g., > 0.03 cm/yr). Combining our data with those from previous studies, we demonstrate that average granitoids and upper continental crust (with newly estimated δ44CaBSE of -0.25 ± 0.02‰, 2SE) have resolvable low δ44Ca compared to basalts and oceanic crust. Given that pressure has a major influence on Ca isotope fractionation across all of our models, this implies that melts feeding upper crustal granitoids dominantly evolve in the lower crust (10-14 kbar, through either partial-melting or fractional crystallization). This observation also suggests that heavier Ca isotopes are preferentially recycled back into the mantle through subduction and/or lower-crustal delamination events, but this is unlikely to have had a significant influence on the δ44Ca evolution of the upper mantle through geologic time.

Equilibrium Ca isotope fractionation and the rates of isotope exchange in the calcite-fluid and aragonite-fluid systems at 25 °C

The Ca isotope composition of calcite and aragonite can provide insights into their formation conditions. The accurate interpretation of the Ca isotope composition of natural samples, however, requires precise knowledge of the equilibrium Ca isotope fractionation between the solid and fluid phases. In this study, the three-isotope method with 42Ca, 43Ca and 44Ca has been used to estimate the equilibrium isotope fractionation of Ca between the CaCO3 minerals calcite and aragonite and the Ca2+ aquo ion. Reactive fluids were enriched with 43Ca and equilibrated with synthetic calcite and aragonite of natural Ca isotope distribution up to 3745 h at 25°C. The isotopic composition of solids and fluids was measured using MC-ICP-MS and the estimated equilibrium fractionation for calcite-Ca2+(aq) and aragonite-Ca2+(aq) was Δ44/42Casolid-fluid=−0.02±0.13‰ and −0.80±0.10‰, respectively. Textural observations of the reacted solids suggest that isotope equilibration in aragonite experiments occurs via extensive Ostwald ripening, yielding large crystalline needles at the end of the experiments. In contrast, calcite did not exhibit an observable increase in size during the course of the experiments. Isotope exchange rates in the case of calcite are similar to those reported in previous studies and ∼4 orders of magnitude lower than the far-from-equilibrium calcite dissolution rate. Calcium isotope exchange rates for aragonite are more rapid than calcite driven by a greater extent of Ostwald ripening occurring via dissolution/precipitation reactions. The results of this study suggest that the Ca isotope compositions of calcite and aragonite crystals in chemical equilibrium, but isotope disequilibrium, with natural fluids could be significantly altered without overt evidence of diagenesis, especially in the case of calcite. The extent to which Ca isotope compositions are altered, however, would depend strongly on the environmental conditions, such as fluid:solid ratio and permeability of the solid facilitating fluid transport.

Calcium isotope ratio in kidney stones: preliminary exploration of mechanism from the geochemical perspective.

Stable calcium (Ca) isotope ratios are sensitive and radiation-free biomarkers in monitoring biological processes in human bodies. Recently, the Ca isotope ratios of bone, blood, and urine have been widely reported to study bone mineral balance. However, as a pure Ca crystallization product, there is no report on the Ca isotope ratios of kidney stones, even though the prevalence of kidney stones is currently on the rise. Here, we measured Ca isotope data of 21 kidney stone samples collected in Beijing, China. The δ44/42CaNIST 915a values ranged from 0.25‰ to 2.85‰ for calcium oxalate, and from 0.38‰ to 3.00‰ and 0.61‰ to 0.69‰ for carbonate apatite and uric acid, respectively. Kidney stones have heavier Ca isotope ratios than bone or blood, which is probably because complexed Ca contains more heavy Ca isotopes than free Ca2+. Ca isotope evidence suggests that magnesium (Mg) affects kidney stone formation, as the δ44/42CaNIST 915a value is inversely correlated with the Ca/Mg ratio. This study provides important preliminary reference values on the Ca isotopic composition of kidney stones and proposes a factor influencing Ca isotope fractionation in biological processes for future research.

Empirical constraints on the mass dependence of isotope diffusion in minerals by modeling sub-solidus exchange: Calcium isotopes in the two-pyroxene system

Diffusion is a significant driver of isotope fractionation at high temperatures, but its precise role is difficult to evaluate in many situations due to the sparsity of data on the isotopic mass dependence of the diffusion coefficient. Such data are lacking particularly for diffusion in minerals. The ratio of diffusion coefficients of two isotopes of the same element in the same material can be characterized as the inverse mass ratio raised to an empirical exponent β, but experiments to determine β are challenging and few measurements have been reported. Here, a method is developed to empirically determine β based on numerical modelling of the diffusion-controlled elemental and isotopic redistribution between minerals during slow subsolidus cooling, and comparison of the results to data from natural mineral pairs. The method is applied to the Ca redistribution between orthopyroxene (opx) and clinopyroxene (cpx), where the elemental partitioning equilibrium as a function of temperature has been well characterized experimentally and Ca isotope fractionation data are available in natural samples that experienced slow cooling. We show that the Ca isotope fractionation during cooling is insensitive to the initial temperature, cooling rate or grain size over a broad range of conditions, provided that diffusion is sufficient to re-homogenize opx grains at near-peak temperatures during early cooling. Furthermore, we find that the Ca isotope fractionation established during cooling is much more strongly dependent on the value of β in opx than the value of β in cpx, as long as cpx dominates the Ca budget in the system. Transient heating events produce fractionation in the opposite sense to that developed during cooling, which can diminish or reverse the isotope fractionation for samples that experienced such events after cooling. Based on the largest Ca isotope fractionation that has been documented between opx and cpx in natural samples, it is inferred that for Ca, βopx = 0.04(1). Despite the relatively small inferred value of βopx, diffusive fractionation of Ca isotopes between opx and cpx is significant over a broad range of cooling rates. Similar behavior during sub-solidus cooling is anticipated for other elements with temperature-dependent inter-mineral partition coefficients, and the method developed here has the potential to be used to place constraints on β for many elements, even in minerals where experimental measurements are not feasible.

Calcium isotope fractionation during the partial dissolution of artificial calcite

We studied the calcium (Ca) isotope fractionation during the partial dissolution of artificial calcite in pure water and in artificial seawater at 20 °C. Calcite (1.00 g) was stirred in 500 ml solution while bubbling it with air rich in CO2 gas (400 and 500 ppm). Depending on the composition of the dissolving medium, the measured pH close to equilibrium decreased with increasing concentration of CO2 in air. The dissolution rate (R*; μmol/m2.h) increased with the concentration of CO2 in air. The Ca isotope composition (δ44/40Ca) of the dissolved Ca in pure water during different intervals of partial dissolution process was almost constant and very close to that of calcite (1.06‰ ± 0.11‰). This phenomenon indicated that partial dissolution occurred without Ca isotope fractionation. By contrast, the isotope composition of Ca in artificial seawater (0.92‰ ± 0.17‰) was almost constant during the partial dissolution of calcite (<10% of calcite dissolved) because the Ca isotope composition in both calcite and in the artificial seawater were close to each other and the change in [Ca2+] ions due to dissolution was not >9%.

Ca-isotopes as a robust tracer of magmatic differentiation

The large mass difference (∼10%) between the two most abundant isotopes of calcium, 40Ca and 44Ca, gives Ca great potential in tracking mass-dependent fractionation during magmatic processes. Resolvable Ca-isotope fractionation during fractional crystallization of magma, particularly by feldspar in evolved melts, has been theoretically inferred but not robustly tested in nature. To further explore the effects of magmatic differentiation on Ca-isotope systematics, we studied the late-Permian alkaline igneous suite of the Øyangen Caldera, Oslo Rift, Norway, consisting of volcanic and intrusive units ranging from basanitic to rhyolitic compositions. Major and trace element variations and modeling demonstrate that the main series of samples (N=21), including basanites, ring-dyke syenites, and central-dome syenites, likely documents a co-genetic and closed-system fractional crystallization sequence. Our data show minimal δ44/40Ca variation (<0.05‰) in the intermediate magma and a marked increase in δ44/40Ca in the felsic magma of the Øyangen Caldera (from 0.62±0.02‰ to 1.15±0.03‰ relative to Ca standard, SRM915a). The systematic increase is best explained by equilibrium isotopic fractionation dominated by alkali feldspar in the fractionating mineral assemblage. This is further supported by strong correlations between δ44/40Ca, CaO, and Eu/Eu* in the main-series samples. Implementing a Monte Carlo approach, isotopic modeling of the liquid line of descent using Rayleigh fractionation is highly consistent with the observed Ca-isotope evolution. For the first time, we confirm prominent Ca stable isotope fractionation in felsic-stage differentiation of alkaline magma and constrain the isotope fractionation factors of plagioclase and K-feldspar. Integrated with extant estimations on mineral fractionation factors from the literature, our results suggest increasing fractionation effects of rock-forming minerals with decreasing Ca content. The affirmation of significant Ca-isotope fractionation in alkaline magma by feldspar empowers the application of Ca as a versatile tracer of crustal evolution, allowing further tests in other magmatic conditions across various planetary objects.