Premium
Models of transport and reaction describing weathering of fractured rock with mobile and immobile water
Author(s) -
Reis Fábio D. A. Aarão,
Brantley Susan L.
Publication year - 2017
Publication title -
journal of geophysical research: earth surface
Language(s) - English
Resource type - Journals
eISSN - 2169-9011
pISSN - 2169-9003
DOI - 10.1002/2016jf004118
Subject(s) - regolith , weathering , bedrock , geology , dissolution , mineralogy , exponential decay , clastic rock , albite , diffusion , soil production function , macropore , soil science , geomorphology , chemistry , quartz , structural basin , soil water , thermodynamics , physics , pedogenesis , paleontology , astrobiology , nuclear physics , mesoporous material , biochemistry , catalysis
We propose models for calculating the rate of chemical weathering of minerals as a function of depth in the weathering zone. A macropore network is assumed to be responsible for the transport of mobile water, which removes soluble weathering products at the interface of that network and the matrix. Conditions of infrequent rainfall (A) and of very frequent rainfall (B) are separately modeled, but both lead to a volumetric weathering rate with the general formν 0 exp − z / ξ , where the amplitude ν 0 and the equilibration length ξ depend on the pore geometry of regolith and on parameters describing transport across the macropore‐matrix interface (A) or mineral dissolution (B). This is obtained with no assumption of steady state for regolith evolution. Extrapolating these end‐members into intermediate conditions, the model is consistent with the exponential decay of regolith production rate versus depth reported by several authors and yields ξ in a variety of regolith types in the range from centimeters to meters. The velocity of the regolith‐bedrock interface also shows an exponential decay with weathering zone thickness and is enhanced by bedrock fractures when compared to models of unfractured bedrock. When the residence time of fluid in the weathering zone, τ f , is large, the bulk weathering rate R b is inversely proportional to this time, and for small τ f , R b reverts to the laboratory rate. Application to compiled data from several sites leads to estimates of a dissolution rate constant and specific surface area consistent with those of albite.