The role of reaction affinity and secondary minerals in regulating chemical weathering rates at the Santa Cruz Soil Chronosequence, California

• Published in: Geochimica et cosmochimica acta 2009-05, Vol.73 (10), p.2804-2831
• Main Authors: Maher, Kate, Steefel, Carl I., White, Art F., Stonestrom, Dave A.
• Format: Article
• Language: English
• ISSN: 0016-7037; 1872-9533
• DOI: 10.1016/j.gca.2009.01.030

In order to explore the reasons for the apparent discrepancy between laboratory and field weathering rates and to determine the extent to which weathering rates are controlled by the approach to thermodynamic equilibrium, secondary mineral precipitation, and flow rates, a multicomponent reactive transport model (CrunchFlow) was used to interpret soil profile development and mineral precipitation and dissolution rates at the 226 ka Marine Terrace Chronosequence near Santa Cruz, CA. Aqueous compositions, fluid chemistry, transport, and mineral abundances are well characterized [White A. F., Schulz M. S., Vivit D. V., Blum A., Stonestrom D. A. and Anderson S. P. (2008) Chemical weathering of a Marine Terrace Chronosequence, Santa Cruz, California. I: interpreting the long-term controls on chemical weathering based on spatial and temporal element and mineral distributions. Geochim. Cosmochim. Acta 72 (1), 36–68] and were used to constrain the reaction rates for the weathering and precipitating minerals in the reactive transport modeling. When primary mineral weathering rates are calculated with either of two experimentally determined rate constants, the nonlinear, parallel rate law formulation of Hellmann and Tisserand [Hellmann R. and Tisserand D. (2006) Dissolution kinetics as a function of the Gibbs free energy of reaction: An experimental study based on albite feldspar. Geochim. Cosmochim. Acta 70 (2), 364–383] or the aluminum inhibition model proposed by Oelkers et al. [Oelkers E. H., Schott J. and Devidal J. L. (1994) The effect of aluminum, pH, and chemical affinity on the rates of aluminosilicate dissolution reactions. Geochim. Cosmochim. Acta 58 (9), 2011–2024], modeling results are consistent with field-scale observations when independently constrained clay precipitation rates are accounted for. Experimental and field rates, therefore, can be reconciled at the Santa Cruz site. Additionally, observed maximum clay abundances in the argillic horizons occur at the depth and time where the reaction fronts of the primary minerals overlap. The modeling indicates that the argillic horizon at Santa Cruz can be explained almost entirely by weathering of primary minerals and in situ clay precipitation accompanied by undersaturation of kaolinite at the top of the profile. The rate constant for kaolinite precipitation was also determined based on model simulations of mineral abundances and dissolved Al, SiO 2(aq) and pH in pore waters. Changes in the rate of kaol