A mathematical model for steady‐state regolith production at constant erosion rate

It has been hypothesized that many soil profiles reach a steady‐state thickness. In this work, such profiles were simulated using a one‐dimensional model of reaction with advective and diffusive solute transport. A model ‘rock’ is considered, consisting of albite that weathers to kaolinite in the presence of chemically inert quartz. The model yields three different steady‐state regimes of weathering. At the lowest erosion rates, a local‐equilibrium regime is established where albite is completely depleted in the weathering zone. This regime is equivalent to the transport‐limited regime described in the literature. With an increase in erosion rate, transition and kinetic regimes are established. In the transition regime, both albite and kaolinite are present in the weathering zone, but albite does not persist to the soil–air interface. In the weathering‐limited regime, here called the kinetic regime, albite persists to the soil–air interface. The steady‐state thickness of regolith decreases with increasing erosion rate in the local equilibrium and transition regimes, but in the kinetic regime, this thickness is independent of erosion rate. Analytical expressions derived from the model are used to show that regolith production rates decrease exponentially with regolith thickness. The steady‐state regolith thickness increases with the Darcy velocity of the pore fluid, and in the local equilibrium regime may vary markedly with small variations in this velocity and erosion rate. In the weathering‐limited regime, the temperature dependences for chemical weathering rates are related to the activation energy for the rate constant for mineral reaction and to the ΔH of dissolution, while for local equilibrium regimes they are related to the ΔH only. The model illustrates how geochemical and geomorphological observations are related for a simple compositional system. The insights provided will be useful in interpreting natural regolith profiles. Copyright © 2010 John Wiley & Sons, Ltd.     

Climate controls on coupled processes of chemical weathering,bioturbation, and sediment transport across hillslopes

Most hillslope studies examining the interplay between climate and earth surface processes tend to be biased towards eroding parts of landscapes. This limitation makes it difficult to assess how entire upland landscapes, which are mosaics of eroding and depositional areas, evolve physio‐chemically as a function of climate. Here we combine new soil geochemical data and published ¹⁰Be‐derived soil production rates to estimate variations in chemical weathering across two eroding‐to‐depositional hillslopes spanning a climate gradient in southeastern Australia. At the warmer and wetter Nunnock River (NR) site, rates of total soil (–3 to –14 g m^‐2 yr^‐1; negative sign indicates mass loss) and saprolite (–18 to –32 g m^‐2 yr^‐1) chemical weathering are uniform across the hillslope transect. Alternatively, the drier hillslope at Frog's Hollow (FH) is characterized by contrasting weathering patterns in eroding soils (–30 to –53 g m^‐2 yr^‐1) vs. depositional soils (+91 g m^‐2 yr^‐1; positive sign indicates mass addition). This difference partly reflects mineral grain size sorting as a result of upslope bioturbation coupled with water‐driven soil erosion, as well as greater vegetative productivity in moister depositional soils. Both of these processes are magnified in the drier climate. The data reveal the importance of linking the erosion–deposition continuum in hillslope weathering studies in order to fully capture the coupled roles of biota and erosion in driving the physical and chemical evolution of hillslopes. Our findings also highlight the potential limitations of applying current weathering models to landscapes where particle‐sorting erosion processes are active. Copyright © 2018 John Wiley & Sons, Ltd.     

Rock damage and regolith transport by frost: an example of climate modulation of the geomorphology of the critical zone

In this article we craft process‐specific algorithms that capture climate control of hillslope evolution in order to elucidate the legacy of past climate on present critical zone architecture and topography. Models of hillslope evolution traditionally comprise rock detachment into the mobile layer, mobile regolith transport, and a channel incision or aggradation boundary condition. We extend this system into the deep critical zone by considering a weathering damage zone below the mobile regolith in which rock strength is diminished; the degree of damage conditions the rate of mobile regolith production. We first discuss generic damage profiles in which appropriate length and damage scales govern profile shapes, and examine their dependence upon exhumation rate. We then introduce climate control through the example of rock damage by frost‐generated crack growth. We augment existing frost cracking models by incorporating damage rate limitations for long transport distances for water to the freezing front. Finally we link the frost cracking damage model, a mobile regolith production rule in which rock entrainment is conditioned by the damage state of the rock, and a frost creep transport model, to examine the evolution of an interfluve under oscillating climate. Aspect‐related differences in mean annual surface temperatures result in differences in bedrock damage rate and mobile regolith transport efficiency, which in turn lead to asymmetries in critical zone architecture and hillslope form (divide migration). In a quasi‐steady state hillslope, the lowering rate is uniform, and the damage profile is better developed on north‐facing slopes where the frost damage process is most intense. Because the residence times of mobile regolith and weathered bedrock in such landscapes are on the order of 10 to 100 ka, climate cycles over similar timescales result in modulation of transport and damage efficiencies. These lead to temporal variation in mobile regolith thickness, and to corresponding changes in sediment delivery to bounding streams. Copyright © 2012 John Wiley & Sons, Ltd.     

Local topography and erosion rate control regolith thickness along a ridgeline in the Sierra Nevada,California

The ridgelines of mountain ranges are a source of geomorphic information unadulterated by the arrival of sediment from upslope. Studies along ridgecrests, therefore, can help identify and isolate the controls on important regolith properties such as thickness and texture. A 1.5 km section of ridgeline in the Sierra Nevada (CA) with a tenfold decrease in erosion rate (inferred from ridgetop convexity) provided an opportunity to conduct a high‐resolution survey of regolith properties and investigate their controls. We found that regolith along the most quickly eroding section of the ridge was the rockiest and had the lowest clay concentrations. Furthermore, a general increase in regolith thickness with a slowing of erosion rate was accompanied by an increase in biomass, changes in vegetation community, broader ridgeline profiles, and an apparent increase in total available moisture. The greatest source of variation in regolith thickness at the 10–100 m scale, however, was the local topography along the ridgeline, with the deepest regolith in the saddles and the thinnest on the knobs. Because regolith in the saddles had higher surface soil moisture than the knobs, we conclude that the hydrological conditions primarily driven by local topography (i.e. rapid vs. slow drainage and water‐storage potential) provide the fundamental controls on regolith thickness through feedbacks incorporating physical weathering by the biota and chemical weathering. Moreover, because the ridgeline saddles are the uppermost extensions of first‐order valleys, we propose that the fluvial network affects regolith properties in the furthest reaches of the watershed. Copyright © 2015 John Wiley & Sons, Ltd.     

Relating the depth of the water table to the depth of weathering

Weathering of bedrock creates and occludes permeability, affecting subsurface water flow. Often, weathering intensifies above the water table. On the contrary, 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 the eventual steady-state position of the water table with respect to the weathering reaction front. Depending on 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). 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. Published 2020. This article is a U.S. Government work and is in the public domain in the USA     

A low‐dimensional model of bedrock weathering and lateral flow coevolution in hillslopes: 2. Controls on weathering and permeability profiles,drainage hydraulics,and solute export pathways

The advance of a chemical weathering front into the bedrock of a hillslope is often limited by the rate weathering products that can be carried away, maintaining chemical disequilibrium. If the weathering front is within the saturated zone, groundwater flow downslope may affect the rate of transport and weathering—however, weathering also modifies the rock permeability and the subsurface potential gradient that drives lateral groundwater flow. This feedback may help explain why there tends to be neither “runaway weathering” to great depth nor exposed bedrock covering much of the earth and may provide a mechanism for weathering front advance to keep pace with incision of adjacent streams into bedrock. This is the second of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory. Here, we show how a simplified kinetic model of 1‐D rock weathering can be extended to consider lateral flow in a 2‐D hillslope. Exact and approximate analytical solutions for the location and thickness of weathering within the hillslope are obtained for a number of cases. A location for the weathering front can be found such that lateral flow is able to export weathering products at the rate required to keep pace with stream incision at steady state. Three pathways of solute export are identified: “diffusing up,” where solutes diffuse up and away from the weathering front into the laterally flowing aquifer; “draining down,” where solutes are advected primarily downward into the unweathered bedrock; and “draining along,” where solutes travel laterally within the weathering zone. For each pathway, a different subsurface topography and overall relief of unweathered bedrock within the hillslope is needed to remove solutes at steady state. The relief each pathway requires depends on the rate of stream incision raised to a different power, such that at a given incision rate, one pathway requires minimal relief and, therefore, likely determines the steady‐state hillslope profile.     

Weathering and erosion of fractured bedrock systems

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, 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. Published 2017. This article is a U.S. Government work and is in the public domain in the USA     

Space weathering processes on airless bodies: Fe isotope fractionation in the lunar regolith

Nanophase Fe metal grains (np-Fe°) are a product of space weathering, formed by processes related to meteorite impacts, and solar-wind sputtering on airless planetary bodies, such as the Moon. Iron isotopes of lunar soils are fractionated during these processes, and the np-Fe° in the finest (<10 μm), mature, size fractions of the soil become enriched in heavier isotopes by ∼0.3‰ in ⁵⁶Fe/⁵⁴Fe in comparison to the bulk rocks (0.03±0.05‰), from which the soil was formed. A positive correlation of δ⁵⁶Fe values with the soil maturity index, I_S/FeO, suggests that the high δ⁵⁶Fe values reflect production of nanophase Fe metal that is produced by space weathering that occurs on airless planetary bodies. Furthermore, the enrichment of δ⁵⁶Fe in the smallest size fraction of lunar soils supports a model for creation of np-Fe° through vapor deposition induced by micrometeorites, as well as that by solar-wind sputtering.     

The concept of effective length in hillslopes: assessing the influence of climate and topography on the contributing areas of catchments

This paper introduces a concept of ‘effective length’ in hillslopes to define the effective area influencing the runoff‐producing saturated zones of a hillslope or catchment. This effective area of a catchment usually is less than that given by its physical boundaries, particularly in regions where the total potential evaporation exceeds total rainfall on an annual basis. In this paper, expressions for effective lengths in hillslopes with different scale, shape and soil properties are derived for given climatic conditions. The influence of these variables on effective length is investigated. It is shown that, for a given rainfall frequency and soil parameters, the effective length changes with the planform geometry and profile shape of a hillslope; it is also a function of the ratio of available travel opportunity time to the hillslope's scale response time. The application of the concept to three natural catchments, subdivided into a number of simple hillslopes, is described. It is shown that, for these three test catchments and over 24 years of record, rarely would the entire catchment areas contribute to flow at the respective outlets. The implications of the concept of effective length for several land‐use practices, such as clearing for forest for greater water yield, and planting trees for salinity control, are discussed. Copyright © 2002 John Wiley & Sons, Ltd.     

Temporal variation of chemical and mechanical weathering in NE Iceland: Evaluation of a steady-state model of erosion

This study critically assesses the temporal sensitivity of the steady-state model of erosion that has been applied to chemical and mechanical weathering studies of volcanic islands and the continents, using only one sample from each catchment. The model assumes a geochemical mass balance between the initially unweathered rock of a drainage basin and the dissolved and solid loads of the river.Chemical composition of 178 samples of suspended and dissolved inorganic river constituents, collected in 1998–2002, were studied from five basaltic river catchments in NE Iceland. The Hydrological Service in Iceland has monitored the discharge and the total suspended inorganic matter concentration (SIM) of the glacial rivers for ~ four decades, making it possible to compare modelled and measured SIM fluxes.Concentration of SIM and grain size increased with discharge. As proportion of clay size particles in the SIM samples increased, concentrations of insoluble elements increased and of soluble decreased. The highest proportion of altered basaltic glass was in the clay size particles.The concentration ratio of insoluble elements in the SIM was used along with data on chemical composition of unweathered rocks (high-Mg basalts, tholeiites, rhyolites) to calculate the pristine composition of the original catchment rocks. The calculated rhyolite proportions compare nicely with area-weighted average proportions, from geological maps of these catchments.The calculated composition of the unweathered bedrock was used in the steady-state model, together with the chemical composition of the suspended and dissolved constituents of the river. Seasonal changes in dissolved constituent concentrations resulted in too low modelled concentrations of SIM_mod at high discharge (and too high SIM_mod at low discharge). Samples collected at annual average river dissolved load yielded SIM_mod concentrations close to the measured ones. According to the model, the studied rivers had specific mechanical denudation rates of 1.3–3.0 kg/m²/yr whereas the average measured rates were 0.8–3.5 kg/m²/yr which are among the highest on Earth.This study validates the use of a steady-state model of erosion to estimate mechanical weathering rates at the scale of a river catchment when the collected riverine dissolved load represents the average chemical composition over a mean hydrological year.     

A low‐dimensional model of bedrock weathering and lateral flow coevolution in hillslopes: 1. Hydraulic theory of reactive transport

This is the first of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory (also known as Dupuit or Boussinesq theory). Here, we examine the effect of lateral flow on the downward fluxes of water and solutes through perched groundwater at steady state. We derive analytical expressions describing the decline in the downward flux rate with depth. Using these, we obtain analytical expressions for water age in a number of cases. The results show that when the permeability field is homogeneous, the spatial structure of water age depends qualitatively on a single dimensionless number, Hi. This number captures the relative contributions to the lateral hydraulic potential gradient of the relief of the lower‐most impermeable boundary (which may be below the weathering front within permeable or incipiently weathered bedrock) and the water table. A “scaled lateral symmetry” exists when Hi is low: age varies primarily in the vertical dimension, and variations in the horizontal dimension x almost disappear when the vertical dimension z is expressed as a fraction z/H(x) of the laterally flowing system thickness H(x). Taking advantage of this symmetry, we show how the lateral dimension of the advection–diffusion‐reaction equation can be collapsed, yielding a 1‐D vertical equation in which the advective flux downward declines with depth. The equation holds even when the permeability field is not homogeneous, as long as the variations in permeability have the same scaled lateral symmetry structure. This new 1‐D approximation is used in the accompanying paper to extend chemical weathering models derived for 1‐D columns to hillslope domains.     

Quantifying weathering on variable rocks,an extension of geochemical mass balance: Critical zone and landscape evolution

We present a statistical model of soil and rock weathering in deep profiles to expand the capacity to assess weathering to heterogeneous bedrock types, which are common at the Earth's surface. We developed the Weathering Trends (WT) model by extending the fractional mass change calculation (tau) of the geochemical mass balance model in two important ways. First, WT log transforms the elemental ratio data, to discern the log‐linear patterns that naturally develop from thermodynamic and kinetic laws of chemistry. Second, WT statistically fits log‐transformed element concentration ratio data – log(c_j/c_i), the only depth‐varying term in tau – as a function of depth to determine characteristic depths of transitions in weathering processes, along with confidence intervals. With no prior assumptions, WT estimates average parent material composition, average composition of the upper weathered zone and mean fractional mass change of each element over the entire weathering profile. WT displays the mean shape of weathering profiles of log‐transformed geochemical data bounded by calculated confidence intervals. We share the WT model code as an open‐source R package ( https://github.com/fisherba/WeatheringTrends ). The WT model was designed to interpret two 21 m cores from the Laurels Schist bedrock in the Christina River Basin Critical Zone Observatory in the Pennsylvania Piedmont, where our morphological and elemental data provided inconclusive estimates of bedrock depth. The WT model differentiated between rock variability and weathering to delineate the maximum extent of weathering at 12.3 m (CI 95% [9.2, 21.3]) in Ridge Well 1 and 7.2 m (CI 95% [4.3, 13.0]) in Interfluve Well 2. The water table was 5–8 m below fresh rock at Ridge Well 1, but at the same depth as fresh rock at the lower elevation interfluve. We assess statistical approaches to identify the best immobile element for use in WT and tau calculations. Copyright © 2017 John Wiley & Sons, Ltd.     

Weathering of granitic rocks by chlorides: effect of the nature of the solution on weathering morphology

Five granites from NW Spain were subjected to two salt‐spray weathering tests under controlled atmospheric conditions. Granite samples were exposed to a sodium chloride solution in the first test and to a complex solution (sea water) in the second. Subsequent examination of the rocks by light and petrographic microscopy clearly demonstrated the development of different weathering morphologies in each test. The distribution of dissolved ions in samples taken at different depths from the weathered surfaces at the end of the experiments, and the changes in weight of the samples during the tests, also differed. In our opinion, these results were determined by the nature of the solution used; although sodium chloride is the predominant salt in sea water, the presence of other ions modifies its solubility, thereby varying the dynamics of mobility and precipitation and thus the weathering morphology that develops. Copyright © 2003 John Wiley & Sons, Ltd.     

A one‐dimensional model for simulating armouring and erosion on hillslopes: 2. Long term erosion and armouring predictions for two contrasting mine spoils

This paper investigates the dynamics of soil armouring as a result of fluvial erosion for a non‐cohesive sandy gravel spoil from the Ranger Mine, Australia, and a cohesive silt loam spoil from the Northparkes Mine, Australia, using a model for hillslope soil armouring. These long term predictions concentrate on the temporal and spatial changes of the spoil grading and erosion over 100–200 years for the flat cap regions (1–2%) and steep batter edges (10–30%) typically encountered on waste rock dumps. The existence of a significant rock fragment fraction in the Ranger spoil means that it armours readily, while Northparkes does not. For Ranger the waste rock showed reductions in (1) cumulative erosion of up to 81% from that obtained by extrapolating the initial erosion rate out 100 years and (2) the erosion/year by more than 10‐fold. For Northparkes reductions were less marked, with the maximum reduction in erosion/year being 37% after 200 years. For Ranger the reductions were greatest and fastest for intermediate gradient hillslopes. For the steepest hillslopes the armouring decreased because the flow shear stresses were large enough to mobilize all material in the armour layer. Model uncertainty was assessed with probabilistic confidence limits demonstrating that these erodibility reductions were statistically significant. A commonly used hillslope erosion model (sediment flux = β₁ discharge m₁ slope n₁) was fitted to these predictions. The erodibility, β₁, and m₁ decreased with time, which was consistent with our physical intuition about armouring. At Ranger the parameter m₁ asymptoted to 1·5–1·6 while at Northparkes it asymptoted to 1·2–1·3. At Ranger transient spatial trends in armouring led to a short term (50–200 years in the future) reduction in n₁, to below zero under certain circumstances, recovering to an asymptote of about 0·5–1. At Northparkes n₁ asymptoted to about 0·6, with no negative transients predicted. The m₁ and n₁ parameters predicted for Ranger were shown to be consistent with field data from a 10‐year‐old armoured hillslope and consistent with published relationships between erodibility and rock content for natural hillslopes. Copyright © 2007 John Wiley & Sons, Ltd.     

Field evidence for the influence of weathering on rock erodibility and channel form in bedrock rivers

Erosion processes in bedrock‐floored rivers shape channel cross‐sectional geometry and the broader landscape. However, the influence of weathering on channel slope and geometry is not well understood. Weathering can produce variation in rock erodibility within channel cross‐sections. Recent numerical modeling results suggest that weathering may preferentially weaken rock on channel banks relative to the thalweg, strongly influencing channel form. Here, we present the first quantitative field study of differential weathering across channel cross‐sections. We hypothesize that average cross‐section erosion rate controls the magnitude of this contrast in weathering between the banks and the thalweg. Erosion rate, in turn, is moderated by the extent to which weathering processes increase bedrock erodibility. We test these hypotheses on tributaries to the Potomac River, Virginia, with inferred erosion rates from ~0.1 m/kyr to >0.8 m/kyr, with higher rates in knickpoints spawned by the migratory Great Falls knickzone. We selected nine channel cross‐sections on three tributaries spanning the full range of erosion rates, and at multiple flow heights we measured (1) rock compressive strength using a Schmidt hammer, (2) rock surface roughness using a contour gage combined with automated photograph analysis, and (3) crack density (crack length/area) at three cross‐sections on one channel. All cross‐sections showed significant (p < 0.01 for strength, p < 0.05 for roughness) increases in weathering by at least one metric with height above the thalweg. These results, assuming that the weathered state of rock is a proxy for erodibility, indicate that rock erodibility varies inversely with bedrock inundation frequency. Differences in weathering between the thalweg and the channel margins tend to decrease as inferred erosion rates increase, leading to variations in channel form related to the interplay of weathering and erosion rate. This observation is consistent with numerical modeling that predicts a strong influence of weathering‐related erodibility on channel morphology. Copyright © 2017 John Wiley & Sons, Ltd.     

Laboratory simulation of the salt weathering of schist: 1. Weathering of schist blocks in a seasonally wet tropical environment

Data describing sediment generation focusing on the temporal evolution of size gradation are required for the prediction of long‐term landform evolution. This paper presents such data for the salt weathering of a quartz‐chlorite schist obtained from the Ranger Uranium Mine in northern Australia. Rock fragment samples are subjected to three different climate regimes: (1) a dry season climate; (2) a wet season climate (both based on observations at the Ranger site); and (3) an oven‐drying sequence designed to test the sensitivity of the weathering process by exposing the rocks to more extreme temperatures. Two MgSO₄ salt solutions are applied, one being typical of wet season runoff and the other a more concentrated solution. Salt solution is applied daily in the wet season experiments and once only at the beginning of the dry season experiments. Results of the experiments reveal four stages of weathering. The kinetics of each stage are described and related to the formation of sediment of different sizes. Wet season climate conditions are shown to produce greater moisture variability and lead to faster weathering rates. However, salt concentrations in the wet season are typically lower and so when climate is combined with observed salt concentrations, the dry and wet season experiments weather at approximately equal rates. Finally, small variations in rock properties were shown to have a large impact on weathering rates, leading to the conclusion that rock weathering experiments need to be carefully designed if results are to be used to help predict weathering behaviour at the landscape scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Use of geochemical mass balance modelling to evaluate the role of weathering in determining stream chemistry in five mid-atlantic watersheds on different lithologies

The importance of mineral weathering was assessed and compared for five mid-Atlantic watersheds receiving similar atmospheric inputs but underlain by differing bedrock. Annual solute mass balances and volume-weighted mean solute concentrations were calculated for each watershed for each year of record. In addition, primary and secondary mineralogy were determined for each of the watersheds through analysis of soil samples and thin sections using petrographic, scanning electron microscope, electron microprobe and X-ray diffraction techniques. Mineralogical data were also compiled from the literature. These data were input to NETPATH, a geochemical program that calculates the masses of minerals that react with precipitation to produce stream water chemistry. The feasibilities of the weathering scenarios calculated by NETPATH were evaluated based on relative abundances and reactivities of minerals in the watershed. In watersheds underlain by reactive bedrocks, weathering reactions explained the stream base cation loading. In the acid-sensitive watersheds on unreactive bedrock, calculated weathering scenarios were not consistent with the abundance of reactive minerals in the underlying bedrock, and alternative sources of base cations are discussed. © 1997 John Wiley & Sons, Ltd.     

Using Rock‐Eval 6 pyrolysis for tracking fossil organic carbon in modern environments: implications for the roles of erosion and weathering

This work relates to the debate on the fossil organic carbon (FOC) input in modern environments and its possible implication for the carbon cycle, and suggests the use of Rock‐Eval 6 pyrolysis as a relevant tool for tracking FOC in such environments. Considering that such a delivery is mainly due to supergene processes affecting the continental surface, we studied organic matter in different reservoirs such as bedrocks, alterites, soils and rivers in two experimental catchments at Draix (Alpes de Haute Provence, France). Samples were subjected to geochemical (Rock‐Eval 6 pyrolysis) investigations and artificial bacterial degradations. After comparing the geochemical fingerprint of samples, geochemical markers of FOC were defined and tracked in the different reservoirs. Our results confirm the contribution of FOC in modern soils and rivers and display the various influences of weathering and erosional processes on the fate of FOC during its exchange between these pools. In addition, the contrasting behaviour of these markers upon the supergene processes has also highlighted the refractory or labile characters of the fossil organic matter (FOM). Bedrock to river fluxes, controlled by gully erosion, are characterized by a qualitative and quantitative preservation of FOM. Bedrock to alterite fluxes, governed by chemical weathering, are characterized by FOC mineralization without qualitative changes in deeper alterites. Alterite to soils fluxes, controlled by (bio)chemical weathering, are characterized by strong FOC mineralization and qualitative changes of FOM. Thus weathering and erosional processes induce different FOM evolution and affect the fate of FOC towards the global carbon cycle. In this study, gully erosion would involve maintenance of an ancient sink for the global carbon cycle, while (bio)chemical processes provide a source of CO₂. Finally, this study suggests that Rock‐Eval 6 pyrolysis can be considered as a relevant tool for tracking FOC in modern environments. Copyright © 2006 John Wiley & Sons, Ltd.     

Rainfall runoff and erosion in Napa Valley vineyards: effects of slope,cover and surface roughness

The effects of slope, cover and surface roughness on rainfall runoff, infiltration and erosion were determined at two sites on a hillside vineyard in Napa County, California, using a portable rainfall simulator. Rainfall simulation experiments were carried out at two sites, with five replications of three slope treatments (5%, 10% and 15%) in a randomized block design at each site (0%bsol;64 m² plots). Prior to initiation of the rainfall simulations, detailed assessments, not considered in previous vineyard studies, of soil slope, cover and surface roughness were conducted. Significant correlations (at the 95% confidence level) between the physical characteristics of slope, cover and surface roughness, with total infiltration, runoff, sediment discharge and average sediment concentration were obtained. The extent of soil cracking, a physical characteristic not directly measured, also affected analysis of the rainfall–runoff–erosion process. Average cumulative runoff and cumulative sediment discharge from site A was 87% and 242% greater, respectively, than at site B. This difference was linked to the greater cover, extent of soil cracking and bulk density at site B than at site A. The extent of soil cover was the dominant factor limiting soil loss when soil cracking was not present. Field slopes within the range of 4–16%, although a statistically significant factor affecting soil losses, had only a minor impact on the amount of soil loss. The Horton infiltration equation fit field data better than the modified Philip's equation. Owing to the variability in the ‘treatment’ parameters affecting the rainfall–runoff–erosion process, use of ANOVA methods were found to be inappropriate; multiple‐factor regression analysis was more useful for identifying significant parameters. Overall, we obtained similar values for soil erosion parameters as those obtained from vineyard erosion studies in Europe. In addition, it appears that results from the small plot studies may be adequately scaled up one to two orders of magnitude in terms of land areas considered. Copyright © 2000 John Wiley & Sons, Ltd.