Two mathematical models of Mg and Sr partitioning into solution during incongruent calcite dissolution: Implications for dripwater and speleothem studies

Abstract Positive correlation between Sr/Ca and Mg/Ca in cave dripwaters and speleothem is often attributed to ‘prior calcite precipitation’ (PCP), where exclusion of Mg and Sr from calcite precipitated upstream of the dripwater/speleothem site results in enrichment in the Mg/Ca and Sr/Ca ratios of residual effluent. Modelling predicts that PCP will produce a straight line correlation in graphs of ln(Sr/Ca) vs. ln(Mg/Ca) with a theoretical slope that is given by ( Kd Sr  − 1) / ( Kd Mg  − 1) (calculated to be ∼ 0.88 ± 0.13) regardless of rock or absolute dripwater composition ( Sinclair et al., submitted for publication ). While this makes the correlation slope potentially diagnostic for PCP, Mg and Sr are also preferentially released from calcite during incongruent calcite dissolution (ICD), and it is possible that ICD produces a similar correlation between ln(Sr/Ca) and ln(Mg/Ca). To examine this possibility, mathematical formulations for two different mechanisms of ICD are presented. Model 1 simulates net dissolution as a dynamic kinetic process where forward and reverse reactions compete: dissolution is congruent but ions are fractionated by a simultaneous reverse reaction where new calcite forms with Sr/Ca and Mg/Ca ratios that are determined by the solution composition and Kd Sr and Kd Mg . Model 2 simulates leaching of Mg and Sr from a fresh calcite surface by tracking the dissolution of a thin layer where high ion mobility allows the calcite to maintain a homogenous composition controlled by Kd Sr and Kd Mg . In both models, generalised kinetic equations are used in the derivations meaning that the results are not specific to a particular kinetic model. In both models of ICD, the mathematics predict that under conditions that may prevail during dissolution of limestone, the correlation between ln(Sr/Ca) and ln(Mg/Ca) is approximately linear with a slope given by ( Kd Sr  − 1) / ( Kd Mg  − 1). Thus, this property is not specific to just PCP. However, for both PCP and ICD, this slope is universal—expected to occur regardless of kinetics and the absolute composition of solution or host rock. A correlation slope of ∼ 0.88 in either karst waters or speleothem is therefore diagnostic for calcite–water interaction.

[1]  J. Tison,et al.  Solute generation and transfer from a chemically reactive alpine glacial-proglacial system , 1999 .

[2]  J. Schroeder.  Experimental dissolution of calcium, magnesium, and strontium from Recent biogenic carbonates; a model of diagenesis , 1969 .

[3]  S. Burns,et al.  Climate influence on geochemistry parameters of waters from Santana–Pérolas cave system, Brazil , 2007 .

[4]  L. Land.  Diagenesis of skeletal carbonates , 1967 .

[5]  D. Shaw.  Trace element fractionation during anatexis , 1970 .

[6]  E. Watson,et al.  A conceptual model for near-surface kinetic controls on the trace-element and stable isotope composition of abiogenic calcite crystals , 2004 .

[7]  D. Schrag.  Rapid analysis of high‐precision Sr/Ca ratios in corals and other marine carbonates , 1999 .

[8]  A. Kulik,et al.  Spontaneous movement of ions through calcite at standard temperature and pressure , 1998, Nature.

[9]  S. Stipp,et al.  The dynamic nature of calcite surfaces in air , 1996 .

[10]  P. W. Gast.  Trace element fractionation and the origin of tholeiitic and alkaline magma types , 1968 .

[11]  H. Elderfield,et al.  An intensity ratio calibration method for the accurate determination of Mg/Ca and Sr/Ca of marine carbonates by ICP‐AES , 2002 .

[12]  E. Bonatti,et al.  High-Magnesian Calcite: Leaching of Magnesium in the Deep Sea , 1970, Science.

[13]  A. Reeve,et al.  Carbonate geochemistry and the concentrations of aqueous Mg2+, Sr2+ and Ca2+: Western north coast of the Yucatan, Mexico , 1994 .

[14]  I. Fairchild,et al.  Partitioning of Sr2+ and Mg2+ into calcite under karst-analogue experimental conditions , 2001 .

[15]  B. Spiro,et al.  Controls on trace element Sr-Mg compositions of carbonate cave waters: implications for speleothem climatic records , 2000 .

[16]  I. Fairchild.  Balmy Shores and Icy Wastes: The Paradox of Carbonates Associated with Glacial Deposits in Neoproterozoic Times , 2009 .

[17]  I. Fairchild,et al.  An experimental study of incongruent dissolution of CaCO3 under analogue glacial conditions , 2005, Journal of Glaciology.

[18]  F. McDermott,et al.  Annual trace element cycles in calcite–aragonite speleothems: evidence of drought in the western Mediterranean 1200–1100 yr BP , 2005 .

[19]  J. Hellstrom,et al.  Multi-proxy constraints on the climatic significance of trace element records from a New Zealand speleothem , 2000 .

[20]  J. Tison,et al.  Hydrochemistry of carbonate terrains in alpine glacial settings , 1994 .

[21]  S. Burns,et al.  Evidence of rainfall variations in Southern Brazil from trace element ratios (Mg/Ca and Sr/Ca) in a Late Pleistocene stalagmite , 2007 .

[22]  Robert B. Lorens,et al.  Sr, Cd, Mn and Co distribution coefficients in calcite as a function of calcite precipitation rate , 1981 .

[23]  V. Wright.  Sedimentology Review 1 , 1993 .

[24]  A. Mucci,et al.  Proton/calcium ion exchange behavior of calcite. , 2009, Physical chemistry chemical physics : PCCP.

[25]  E. Busenberg,et al.  The kinetics of dissolution of dolomite in CO 2 -H 2 O systems at 1.5 to 65 degrees C and O to 1 atm PCO 2 , 1982 .

[26]  A. Starinsky,et al.  Behaviour of strontium in subsurface calcium chloride brines: Southern Israel and Dead Sea rift valley , 1979 .

[27]  J. Banner,et al.  Magnesium and strontium systematics in tropical speleothems from the Western Pacific , 2012 .