Abstract
Physically based equations for unsaturated groundwater flow and solute transport have been coupled with kinetic rate laws for mineral dissolution–precipitation, and mass balance/mass action equations for aqueous species, in a numerical model that is capable of simulating rock–water interactions in a weathering profile subjected to fluctuating boundary conditions. A numerical experiment was conducted to demonstrate how incipient soil development may proceed in a warm subhumid environment. The simulation involved a hypothetical coarse-textured parent material that was subjected to frequent wetting and drying during an annual water cycle. The hypothetical weathering profile evolved rapidly; dissolution of primary minerals (enstatite, forsterite, and diopside) and precipitation of secondary clay–minerals (kaolinite and Ca-montmorillonite) occurred monotonically despite the abrupt fluctuations in soil-moisture content. In contrast, the activities of aqueous species and dissolution–precipitation rates of calcite were very sensitive to the changing moisture conditions in the upper part of the profile. Although the simulation involved numerous simplifying assumptions, reasonable results were achieved and the calculated (from the model) rate of chemical denudation fell within the range of contemporary denudation rates determined from the dissolved loads of rivers.