5 Search Results

Weathering of bedrock creates and occludes permeability, affecting subsurface water flow. Often, weathering intensifies above the water table. On the other hand, weathering can also commence below the water table. To explore relationships between weathering and the water table, a simplified weathering model for an eroding hillslope was formulated that takes into account both saturated and unsaturated subsurface water flow (but does not fully account for changes in dissolved gas chemistry). The phreatic line was calculated using solutions to mathematical treatments for both zones. In the model, the infiltration rate at the hill surface sets both the original and themore » eventual steady-state position of the water table with respect to the weathering reaction front. Depending upon parameters, the weathering front can locate either above or below the water table at steady state. Erosion also affects the water table position by changing porosity and permeability even when other hydrological conditions (e.g. hydraulic conductivity of parent material, infiltration rate at the surface) do not change. The total porosity in a hill (water storage capacity) was found to increase with infiltration rate (all else held constant). Furthermore, this effect was diminished by increasing the erosion rate. We also show examples of how the infiltration rate affects the position of the water table and how infiltration rate affects weathering advance.« less

We introduce a model of chemical reaction within hills to explore how evolving porosity (and by inference, permeability) affects flow pathways and weathering. The model consists of hydrologic and reactive-transport equations that describe alteration of ferrous minerals and feldspar. These reactions were chosen because previous work emphasized that oxygen- and acid-driven weathering affects porosity differently in felsic and mafic rocks. A parameter controlling the order of the fronts is presented. In the absence of erosion, the two reaction fronts move at different velocities and the relative depths depend on geochemical conditions and starting composition. In turn, the fronts and associatedmore » changes in porosity drastically affect both the vertical and lateral velocities of water flow. For these instances, estimates of weathering advance rates based on simple models that posit unidirectional constant-velocity advection do not apply. In the model hills, weathering advance rates diminish with time as the Darcy velocities decrease with depth. The system can thus attain a dynamical steady state at any erosion rate where the regolith thickness is constant in time and velocities of both fronts become equal to one another and to the erosion rate. The slower the advection velocities in a system, the faster it attains a steady state. For example, a massive rock with relatively fast-dissolving minerals such as diabase -- where solute transport across the reaction front mainly occurs by diffusion -- can reach a steady state more quickly than granitoid rocks in which advection contributes to solute transport. The attainment of a steady state is controlled by coupling between weathering and hydrologic processes that force water to flow horizontally above reaction fronts where permeability changes rapidly with depth and acts as a partial barrier to fluid flow.« less

In this paper we clarify that our weathering model from 2013 did not explicitly describe weathering of soil moving downhill along hillslopes. Additionally, we re-analyze the role of the term that we neglected that describes loss of regolith mass through mineral dissolution. We derive an equation for this term by including lateral flow of water inside the model hill. For the revised hill model, we define a dimensionless parameter that allows estimation of the effect of lateral flow on the steady-state hillslope. This parameter is equal to the ratio of averaged advective flux of dissolved species out of the hillmore » to the rate of total denudation. The parameter also yields a criterion for the existence of a steady state regolith thickness for systems experiencing unidirectional advection at a constant velocity: for a ridge, the rate of downward flow of water (qy) must be less than the rate of upward movement of rock (E) after normalization by a small parameter, α. This parameter is equal to the equilibrium aqueous concentration divided by the concentration of the reacting mineral in the rock. Alternately, a steady state may exist for the case of both vertical and lateral flow in a hill for any value of erosion rate if the Darcy velocities decrease with depth. Subsurface flow systems play an essential role in the existence of both steady-state hillslopes and steady-state regolith thicknesses.« less

We explore the contribution of fractures (joints) in controlling the rate of weathering advance for a low-porosity rock by using methods of homogenization to create averaged weathering equations. The rate of advance of the weathering front can be expressed as the same rate observed in non-fractured media (or in an individual block) divided by the volume fraction of non-fractured blocks in the fractured parent material. In the model, the parent has fractures that are filled with a more porous material that contains only inert or completely weathered material. The low-porosity rock weathers by reaction-transport processes. As observed in field systems,more » the model shows that the weathering advance rate is greater for the fractured as compared to the analogous non-fractured system because the volume fraction of blocks is <1. The increase in advance rate is attributed both to the increase in weathered material that accompanies higher fracture density, and to the increase in exposure of surface of low-porosity rock to reaction-transport. For constant fracture aperture, the weathering advance rate increases when the fracture spacing decreases. Equations describing weathering advance rate are summarized in the “List of selected equations”. If erosion is imposed at a constant rate, the weathering systems with fracture-bounded bedrock blocks attain a steady state. In the erosional transport-limited regime, bedrock blocks no longer emerge at the air-regolith boundary because they weather away. In the weathering-limited (or kinetic) regime, blocks of various size become exhumed at the surface and the average size of these exposed blocks increases with the erosion rate. For convex hillslopes, the block size exposed at the surface increases downslope. This model can explain observations of exhumed rocks weathering in the Luquillo mountains of Puerto Rico.« less