Modelling weathering processes at the catchment scale: The WITCH numerical model

A numerical model of chemical weathering in soil horizons and underlying bedrock (WITCH) has been coupled to a numerical model of water and carbon cycles in forest ecosystems (ASPECTS) to simulate the concentration of major species within the soil horizons and the stream of the Strengbach granitic watershed, located in the Vosges Mountains (France). For the first time, simulations of solute concentrations in soil layers and in the catchment river have been performed on a seasonal basis. The model is able to reproduce the concentrations of most major species within the soil horizons, as well as catching the first-order seasonal fluctuations of aqueous calcium, magnesium and silica concentrations. However, the WITCH model underestimates concentrations of Mg2+ and silica at the spring of the catchment stream, and significantly underestimates Ca2+ concentration. The deficit in calculated calcium can be compensated for by dissolution of trace apatite disseminated in the bedrock. However, the resulting increased Ca2+ release yields important smectite precipitation in the deepest model layer (in contact with the bedrock) and subsequent removal of large amount of silica and magnesium from solution. In contrast, the model accurately accounts for the concentrations of major species (Ca, Mg and silica) measured in the catchment stream when precipitation of clay minerals is not allowed. The model underestimation of Mg2+ and H4SiO4 concentrations when precipitation of well crystallized smectites is allowed strongly suggests that precipitation of well crystallized clay minerals is overestimated and that more soluble poorly crystallized and amorphous materials may be forming. In agreement with observations on other watersheds draining granitic rocks, this study indicates that highly soluble trace calcic phases control the aqueous calcium budget in the Strengbach watershed.