Carbonate dissolution, the process by which carbonate minerals are stoichiometrically dissolved by an aqueous fluid, likely plays an important role in liberating carbon from subducting slabs and may thus make significant contributions to global carbon fluxes. It is therefore necessary to understand the fluid infiltration and flow geometries associated with this process. Despite its widespread use as a powerful tracer of fluid-rock interaction, strontium (Sr) has yet to be widely applied to carbonate dissolution scenarios. Previous work has identified a partially altered metacarbonate layer from Syros, Greece, that has undergone significant carbonate dissolution leading to ∼90% CO2 loss (Ague and Nicolescu, 2014). A bulk-rock compositional profile along this layer has a distinctive spike-shaped Sr profile. Sr concentrations are uniform in the altered region behind the reaction front, increase to form a spike just ahead of the front, and taper off to a baseline value moving farther into the unaltered metacarbonate rock. The origin of this type of spike shape and its potential utility for better understanding carbonate dissolution processes remain unaddressed. Here, we present a numerical diffusion-reaction model which demonstrates that both carbonate dissolution and multidimensional transport of Sr, in which the primary advective flow direction is perpendicular to the propagation of the diffusion front, are required to produce this shape. This model, as well as newly acquired Sr concentration data for calcite (former aragonite) and epidote, also indicate that the infiltration of a fluid with an elevated Sr concentration is necessary. We suggest that such spike-shaped features may aid in identifying carbonate dissolution in other settings, such as vein-selvage systems in metacarbonate rocks or reaction rinds on carbonate-bearing mélange blocks.
Read full abstract