A reactive-transport model for weathering rind formation on basalt

Alexis Navarre-Sitchler, Carl I. Steefel, Peter Benjamin Sak, Susan Louise Brantley

Research output: Contribution to journalArticle

67 Citations (Scopus)

Abstract

Saprolite formation rates influence many important geological and environmental issues ranging from agricultural productivity to landscape evolution. Here we investigate the chemical and physical transformations that occur during weathering by studying small-scale "saprolites" in the form of weathering rinds, which form on rock in soil or saprolite and grow in thickness without physical disturbance with time. We compare detailed observations of weathered basalt clasts from a chronosequence of alluvial terraces in Costa Rica to diffusion-reaction simulations of rind formation using the fully coupled reactive transport model CrunchFlow. The four characteristic features of the weathered basalts which were specifically used as criteria for model comparisons include (1) the mineralogy of weathering products, (2) weathering rind thickness, (3) the coincidence of plagioclase and augite reaction fronts, and (4) the thickness of the zones of mineral reaction, i.e. reaction fronts. Four model scenarios were completed with varying levels of complexity and degrees of success in matching the observations. To fit the model to all four criteria, however, it was necessary to (1) treat diffusivity using a threshold in which it increased once porosity exceeded a critical value of 9%, and (2) treat mineral surface area as a fitting factor. This latter approach was presumably necessary because the mineral-water surface area of the connected (accessible) porosity in the Costa Rica samples is much less than the total porosity (Navarre-Sitchler et al., 2009). The model-fit surface area, here termed reacting surface area, was much smaller than the BET-measured surface area determined for powdered basaltic material. In the parent basalt, reacting surface area and diffusivity are low due to low pore connectivity, and early weathering is therefore transport controlled. However, as pore connectivity increases as a result of weathering, the reacting surface area and diffusivity also increase and weathering becomes controlled by mineral reaction kinetics. The transition point between transport and kinetic control appears to be related to a critical porosity (9%) at which pore connectivity is high enough to allow rapid transport. Based on these simulations, we argue that the rate of weathering front advance is controlled by the rate at which porosity is created in the weathering interface, and that this porosity increases because of mineral dissolution following a rate that is largely surface-reaction controlled.

Original languageEnglish (US)
Pages (from-to)7644-7667
Number of pages24
JournalGeochimica et Cosmochimica Acta
Volume75
Issue number23
DOIs
StatePublished - Dec 1 2011

Fingerprint

reactive transport
Weathering
weathering
basalt
surface area
Porosity
porosity
Minerals
diffusivity
connectivity
saprolite
mineral
Mineral Waters
physical disturbance
Mineralogy
landscape evolution
augite
chronosequence
Surface reactions
reaction kinetics

All Science Journal Classification (ASJC) codes

  • Geochemistry and Petrology

Cite this

Navarre-Sitchler, Alexis ; Steefel, Carl I. ; Sak, Peter Benjamin ; Brantley, Susan Louise. / A reactive-transport model for weathering rind formation on basalt. In: Geochimica et Cosmochimica Acta. 2011 ; Vol. 75, No. 23. pp. 7644-7667.
@article{3fee543375ac4266a2daa4e730571d1b,
title = "A reactive-transport model for weathering rind formation on basalt",
abstract = "Saprolite formation rates influence many important geological and environmental issues ranging from agricultural productivity to landscape evolution. Here we investigate the chemical and physical transformations that occur during weathering by studying small-scale {"}saprolites{"} in the form of weathering rinds, which form on rock in soil or saprolite and grow in thickness without physical disturbance with time. We compare detailed observations of weathered basalt clasts from a chronosequence of alluvial terraces in Costa Rica to diffusion-reaction simulations of rind formation using the fully coupled reactive transport model CrunchFlow. The four characteristic features of the weathered basalts which were specifically used as criteria for model comparisons include (1) the mineralogy of weathering products, (2) weathering rind thickness, (3) the coincidence of plagioclase and augite reaction fronts, and (4) the thickness of the zones of mineral reaction, i.e. reaction fronts. Four model scenarios were completed with varying levels of complexity and degrees of success in matching the observations. To fit the model to all four criteria, however, it was necessary to (1) treat diffusivity using a threshold in which it increased once porosity exceeded a critical value of 9{\%}, and (2) treat mineral surface area as a fitting factor. This latter approach was presumably necessary because the mineral-water surface area of the connected (accessible) porosity in the Costa Rica samples is much less than the total porosity (Navarre-Sitchler et al., 2009). The model-fit surface area, here termed reacting surface area, was much smaller than the BET-measured surface area determined for powdered basaltic material. In the parent basalt, reacting surface area and diffusivity are low due to low pore connectivity, and early weathering is therefore transport controlled. However, as pore connectivity increases as a result of weathering, the reacting surface area and diffusivity also increase and weathering becomes controlled by mineral reaction kinetics. The transition point between transport and kinetic control appears to be related to a critical porosity (9{\%}) at which pore connectivity is high enough to allow rapid transport. Based on these simulations, we argue that the rate of weathering front advance is controlled by the rate at which porosity is created in the weathering interface, and that this porosity increases because of mineral dissolution following a rate that is largely surface-reaction controlled.",
author = "Alexis Navarre-Sitchler and Steefel, {Carl I.} and Sak, {Peter Benjamin} and Brantley, {Susan Louise}",
year = "2011",
month = "12",
day = "1",
doi = "10.1016/j.gca.2011.09.033",
language = "English (US)",
volume = "75",
pages = "7644--7667",
journal = "Geochmica et Cosmochimica Acta",
issn = "0016-7037",
publisher = "Elsevier Limited",
number = "23",

}

A reactive-transport model for weathering rind formation on basalt. / Navarre-Sitchler, Alexis; Steefel, Carl I.; Sak, Peter Benjamin; Brantley, Susan Louise.

In: Geochimica et Cosmochimica Acta, Vol. 75, No. 23, 01.12.2011, p. 7644-7667.

Research output: Contribution to journalArticle

TY - JOUR

T1 - A reactive-transport model for weathering rind formation on basalt

AU - Navarre-Sitchler, Alexis

AU - Steefel, Carl I.

AU - Sak, Peter Benjamin

AU - Brantley, Susan Louise

PY - 2011/12/1

Y1 - 2011/12/1

N2 - Saprolite formation rates influence many important geological and environmental issues ranging from agricultural productivity to landscape evolution. Here we investigate the chemical and physical transformations that occur during weathering by studying small-scale "saprolites" in the form of weathering rinds, which form on rock in soil or saprolite and grow in thickness without physical disturbance with time. We compare detailed observations of weathered basalt clasts from a chronosequence of alluvial terraces in Costa Rica to diffusion-reaction simulations of rind formation using the fully coupled reactive transport model CrunchFlow. The four characteristic features of the weathered basalts which were specifically used as criteria for model comparisons include (1) the mineralogy of weathering products, (2) weathering rind thickness, (3) the coincidence of plagioclase and augite reaction fronts, and (4) the thickness of the zones of mineral reaction, i.e. reaction fronts. Four model scenarios were completed with varying levels of complexity and degrees of success in matching the observations. To fit the model to all four criteria, however, it was necessary to (1) treat diffusivity using a threshold in which it increased once porosity exceeded a critical value of 9%, and (2) treat mineral surface area as a fitting factor. This latter approach was presumably necessary because the mineral-water surface area of the connected (accessible) porosity in the Costa Rica samples is much less than the total porosity (Navarre-Sitchler et al., 2009). The model-fit surface area, here termed reacting surface area, was much smaller than the BET-measured surface area determined for powdered basaltic material. In the parent basalt, reacting surface area and diffusivity are low due to low pore connectivity, and early weathering is therefore transport controlled. However, as pore connectivity increases as a result of weathering, the reacting surface area and diffusivity also increase and weathering becomes controlled by mineral reaction kinetics. The transition point between transport and kinetic control appears to be related to a critical porosity (9%) at which pore connectivity is high enough to allow rapid transport. Based on these simulations, we argue that the rate of weathering front advance is controlled by the rate at which porosity is created in the weathering interface, and that this porosity increases because of mineral dissolution following a rate that is largely surface-reaction controlled.

AB - Saprolite formation rates influence many important geological and environmental issues ranging from agricultural productivity to landscape evolution. Here we investigate the chemical and physical transformations that occur during weathering by studying small-scale "saprolites" in the form of weathering rinds, which form on rock in soil or saprolite and grow in thickness without physical disturbance with time. We compare detailed observations of weathered basalt clasts from a chronosequence of alluvial terraces in Costa Rica to diffusion-reaction simulations of rind formation using the fully coupled reactive transport model CrunchFlow. The four characteristic features of the weathered basalts which were specifically used as criteria for model comparisons include (1) the mineralogy of weathering products, (2) weathering rind thickness, (3) the coincidence of plagioclase and augite reaction fronts, and (4) the thickness of the zones of mineral reaction, i.e. reaction fronts. Four model scenarios were completed with varying levels of complexity and degrees of success in matching the observations. To fit the model to all four criteria, however, it was necessary to (1) treat diffusivity using a threshold in which it increased once porosity exceeded a critical value of 9%, and (2) treat mineral surface area as a fitting factor. This latter approach was presumably necessary because the mineral-water surface area of the connected (accessible) porosity in the Costa Rica samples is much less than the total porosity (Navarre-Sitchler et al., 2009). The model-fit surface area, here termed reacting surface area, was much smaller than the BET-measured surface area determined for powdered basaltic material. In the parent basalt, reacting surface area and diffusivity are low due to low pore connectivity, and early weathering is therefore transport controlled. However, as pore connectivity increases as a result of weathering, the reacting surface area and diffusivity also increase and weathering becomes controlled by mineral reaction kinetics. The transition point between transport and kinetic control appears to be related to a critical porosity (9%) at which pore connectivity is high enough to allow rapid transport. Based on these simulations, we argue that the rate of weathering front advance is controlled by the rate at which porosity is created in the weathering interface, and that this porosity increases because of mineral dissolution following a rate that is largely surface-reaction controlled.

UR - http://www.scopus.com/inward/record.url?scp=80355145701&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=80355145701&partnerID=8YFLogxK

U2 - 10.1016/j.gca.2011.09.033

DO - 10.1016/j.gca.2011.09.033

M3 - Article

AN - SCOPUS:80355145701

VL - 75

SP - 7644

EP - 7667

JO - Geochmica et Cosmochimica Acta

JF - Geochmica et Cosmochimica Acta

SN - 0016-7037

IS - 23

ER -