Geochimica et Cosmochimica Acta, Sept 15, 2012, Vol.93, p.235(27)
To link to full-text access for this article, visit this link: http://dx.doi.org/10.1016/j.gca.2012.03.021 Byline: Joel Moore (a)(b), Peter C. Lichtner (c), Art F. White (d), Susan L. Brantley (a)(b) Abstract: The reactive transport model FLOTRAN was used to forward-model weathering profiles developed on granitic outwash alluvium over 40-3000ka from the Merced, California (USA) chronosequence as well as deep granitic regolith developed over 800ka near Davis Run, Virginia (USA). Baseline model predictions that used laboratory rate constants (k.sub.m), measured fluid flow velocities (v), and BET volumetric surface areas for the parent material (A.sub.B,m.sup.o) were not consistent with measured profiles of plagioclase, potassium feldspar, and quartz. Reaction fronts predicted by the baseline model are deeper and thinner than the observed, consistent with faster rates of reaction in the model. Reaction front depth in the model depended mostly upon saturated versus unsaturated hydrologic flow conditions, rate constants controlling precipitation of secondary minerals, and the average fluid flow velocity (v.sub.a). Unsaturated hydrologic flow conditions (relatively open with respect to CO.sub.2(g)) resulted in the prediction of deeper reaction fronts and significant differences in the separation between plagioclase and potassium feldspar reaction fronts compared to saturated hydrologic flow (relatively closed with respect to CO.sub.2(g)). Under saturated or unsaturated flow conditions, the rate constant that controls precipitation rates of secondary minerals must be reduced relative to laboratory rate constants to match observed reaction front depths and measured pore water chemistry. Additionally, to match the observed reaction front depths, v.sub.a was set lower than the measured value, v, for three of the four profiles. The reaction front gradients in mineralogy and pore fluid chemistry could only be modeled accurately by adjusting values of the product k.sub.m A.sub.B,m.sup.o. By assuming k.sub.m values were constrained by laboratory data, field observations were modeled successfully with TST-like rate equations by dividing measured values of A.sub.B,m.sup.o by factors from 50 to 1700. Alternately, with sigmoidal or Al-inhibition rate models, this adjustment factor ranges from 5 to 170. Best-fit models of the wetter, hydrologically saturated Davis Run profile required a smaller adjustment to A.sub.B,m.sup.o than the drier hydrologically unsaturated Merced profiles. We attributed the need for large adjustments in v.sub.a and A.sub.B,m.sup.o necessary for the Merced models to more complex hydrologic flow that decreased the reactive surface area in contact with bulk flow water, e.g., dead-end pore spaces containing fluids that are near or at chemical equilibrium. Thus, rate models from the laboratory can successfully predict weathering over millions of years, but work is needed to understand how to incorporate changes in what controls the relationship between reactive surface area and hydrologic flow. Author Affiliation: (a) Department of Geosciences, Penn State University, University Park, PA, United States (b) Center for Environmental Kinetics Analysis, Earth and Environmental Systems Institute, Penn State University, University Park, PA, United States (c) Los Alamos National Laboratory, Los Alamos, NM, United States (d) US Geological Survey, Menlo Park, CA, United States Article History: Received 18 October 2010; Accepted 14 March 2012 Article Note: (miscellaneous) Associate editor: Chen Zhu
Precipitation (Meteorology) -- Usage ; Precipitation (Meteorology) -- Analysis ; Flow (Dynamics) -- Usage ; Flow (Dynamics) -- Analysis
Cengage Learning, Inc.