Grassland Hillslope Research Articles

Phosphorus and iron-oxide transport from a hydrologically isolated grassland hillslope

Dissolved reactive phosphorus (DRP) loss from agricultural soils can negatively affect water quality. Shallow subsurface pathways can dominate P losses in grassland soils, especially in wetter months when waterlogging is common. This study investigated the processes controlling intra- and inter-event and seasonal DRP losses from poorly drained permanent grassland hillslope plots. Temporal flow related water samples were taken from surface runoff and subsurface (in-field pipe) discharge, analysed, and related to the likelihood of anaerobic conditions and redoximorphic species including nitrate (NO3−) over time. Subsurface drainage accounted for 89% of total losses. Simple linear regression and correlation matrices showed positive relationships between DRP and iron and soil moisture deficit; and negative relationships between these three factors and NO3− concentrations in drainage. These data indicate that waterlogging and low NO3− concentrations control the release of P in drainage, potentially via reductive dissolution. The relationship between DRP and metal release was less obvious in surface runoff, as nutrients gathered from P-rich topsoil camoflaged redox reactions. The data suggest a threshold in NO3− concentrations that could exacerbate P losses, even in low P soils. Knowledge of how nutrients interact with soil drainage throughout the year can be used to better time soil N and P inputs via, for example, fertiliser or grazing to avoid to excessive P loss that could harm water quality.

Decoupled cycling of carbon, nitrogen, and phosphorus in a grassland soil along a hillslope mediated by clay and soil moisture

Grasslands extend across a variety of topographies including non-flat hilly areas with varied soil texture and moisture that can mediate soil biogeochemical cycling of carbon (C), nitrogen (N) and phosphorus (P). In this study we examined soil organic C, total N and P pools (both inorganic and organic), as well as gross N mineralisation (GNM, as a measure of microbial activity and N dynamics) and microbial C, N and P in winter and spring along a grassland hillslope. Our results showed that variation in soil clay and associated soil moisture were the most important mediators affecting C and N dynamics along the grassland elevational gradient. Total organic C pools were highest where the clay content was highest, likely because of the sorption capacity of clay stabilising organic matter against microbial decomposition. Likewise, the variation in soil organic N pools along the hillslope was mostly controlled by clay and soil moisture effects on microbial stabilisation of organic N. Furthermore, microbial C, N and GNM increased with increased soil moisture, suggesting that microbial activity was limited by soil moisture. In contrast, soil P pools were not related to clay and soil moisture, while soil organic P was decoupled from soil organic C and N pools along the hillslope. Possibly, differences in microbial stabilisation and mineralisation pathways (oxidation for C and N, hydrolysis for P), microbial homeostatic regulation of C, N and P, and P fixation resulted in a decoupling of P from C and N cycling in this grassland. We conclude that variation in clay and associated soil moisture along grassland hillslopes can have important but non-unidirectional effects on soil C, N and P cycling.

Stable isotope profiles of soil organic carbon in forested and grassland landscapes in the Lake Alaotra basin (Madagascar): insights in past vegetation changes

Abstract. The extent to which the central highlands of Madagascar were once covered by forests is still a matter of debate: while reconstructing past environments is inherently difficult, the debate is further hampered by the fact that the evidence documenting land cover changes and their effects on carbon and sediment dynamics in Madagascar has hitherto mainly been derived from lake coring studies. Such studies provide an integrated view over relatively large areas but do not provide information on how land-use change affects hillslopes in terms of carbon and sediment dynamics. Such information would not only be complementary to lake inventories but may also help to correctly interpret lake sediment data. Carbon stable isotope ratios (δ13C) are particularly useful tracers to study the past dynamics of soil carbon over time spans ranging from years to millennia and thus to understand the consequences of land-use change over such time spans. We analysed soil profiles down to a depth of 2 m from pristine forests and grasslands in the Lake Alaotra region in central Madagascar. Along grassland hillslopes, soil organic carbon (SOC) content was low, from 0.4 % to 1.7 % in the top layer, and decreased rapidly to ca. 0.2 % below 100 cm depth. The current vegetation predominantly consists of C4 grasses (δ13C ∼ −13 ‰), yet topsoil δ13C-OC ranges between −23.0 ‰ and −15.8 ‰, and most profiles show a decrease in δ13C-OC with depth. This contrasts with our observations in the C3-dominated forest profiles, which show a typical profile whereby δ13C values increase slightly with depth. Moreover, the SOC stock of grasslands was ∼ 55.6 % lower than along the forested hillslopes for the upper 0–30 cm layer. δ13C values in grassland and forest profiles converge to similar values (within 2.0 ± 1.8 ‰) at depths below ∼ 80 cm, suggesting that the grasslands in the Lake Alaotra region have indeed developed on soils formerly covered by a tree vegetation dominated by C3 plants. We also observed that the percent of modern carbon (pMC) of the bulk OC in the top, middle and lower middle positions of grasslands was less than 85 % near the surface. This could reflect a combination of (i) the long residence time of forest OC in the soil, (ii) the slow replacement rate of grassland-derived OC (iii) and the substantial erosion of the top positions towards the valley position of grasslands. At the valley positions under grassland, the upper 80 cm contains higher amounts of recent grass-derived OC in comparison to the hillslope positions. This is likely to be related to the higher productivity of the grassland valleys (due to higher moisture and nutrient availability), and the deposition of OC that was eroded further upslope may also have contributed. The method we applied, which is based on the large difference in δ13C values between the two major photosynthetic pathways (C3 and C4) in (sub-)tropical terrestrial environments, provides a relatively straightforward approach to quantitatively determine changing vegetation cover, and we advocate for its broader application across Madagascar to better understand the island's vegetation history.

Experimental analysis of soil moisture response to rainfall in a typical grassland hillslope under different vegetation treatments

The responses of soil moisture to rainfall are of great significance for watershed hydrological modeling. However, few studies have been done to investigate these responds on hillslope in a typical semi-arid grassland region. This study used high temporal resolution soil moisture data to explore the soil moisture dynamics, response conditions and its controls of 0–40 cm soil profile in the upslope (14°), midslope (9°), and downslope (4°) of a typical grassland inland river basin under bare ground (BG), stubble (SG), and natural grassland (CK) treatments. The results showed that soil water content and water storage increased in the downslope direction, and all showed as BG > SG > CK. The dry and wet changes in fast-changing layer (5 cm) and active layer (10 cm) were rapid, while soil moisture below 20 cm was relatively stable and fluctuated only in heavy or continuous rainfalls. The soil moisture response process varied greatly under different rainfall, rainfall intensity and antecedent soil moisture conditions, which explained 41.1% of the total difference. The rainfall replenishment threshold and the required initial soil profile water content of soil moisture response in 5 cm, 10 cm and 20 cm soil layers were 5.8 mm, 8.0 mm, 11.4 mm and 8.7 vol%, 9.4 vol%, 10.8 vol%, respectively. Soil properties, vegetation characteristics and topography could explain 38.8%, 14.5% and 5.6% of the soil moisture variation on the hillslope. In addition, under the comprehensive influence of environmental factors, changes in soil moisture of the upslope were significantly affected by soil sand content, the differences in the midslope were mainly due to soil clay content and belowground biomass, whereas the vegetation characteristics were the main factors in the downslope. This study can contribute to the further understanding of slope-scale ecohydrological processes and hydrological simulation of semi-arid grassland watersheds.

Effect of rainfall gradient and vegetation restoration on gully initiation under a large-scale extreme rainfall event on the hilly Loess Plateau: A case study from the Wuding River basin, China

As the most serious form of soil erosion, gully erosion can be triggered by individual high-intensity rainfall events. In this study, a total of 369 small catchments in 24 sites were sampled to investigate the relationship between rainfall and gully erosion on hillslopes and to study the impacts of vegetation restoration following heavy rainstorms in the central Loess Plateau, China. A total of 280 newly formed gullies on hillslopes were identified by comparing pre-storm Google Earth images and post-storm unmanned aerial vehicle (UAV) images. The results showed that the dimensions and density of gullies increased significantly with rainfall gradient increasing from the periphery to the storm center. When the rainfall amount exceeded 200 mm, gully volumetric density reached up to 928.39 m3/km2 and the mean gully volume was 15.74 m3, 12.8 times and 2.3 times the mean gully volume for rainfall amounts of 106 and 150 mm, respectively. In the sampled small catchments, where cropland was dominant, the relationships between the gully densities and rainfall amount could be fitted with exponential functions. Vegetation restoration was found to reduce the densities and dimensions of gullies on hillslopes. Compared to those in cropland-dominated catchments, the density of gullies in grassland-dominated catchments was found to be lower by > 60%, while the individual gully volume was found to be 1.6 times higher. In small catchments, no new hillslope gullies were observed when the rainfall amount fell below 106.7 mm. Therefore, the rainfall thresholds for (1) ephemeral-gully initiation on grassland hillslopes, (2) permanent-gully initiation on grassland hillslopes, and (3) permanent-gully initiation on cropland hillslopes are concluded to be not >106.7 mm, not >136.1 mm, and not >110.2 mm, respectively. This suggests that the restoration of cropland to grassland would reduce the rainfall threshold for gully initiation.

Spatial and Temporal Dynamics of Hillslope‐Scale Soil Moisture Patterns: Characteristic States and Transition Mechanisms

Recent advances in wireless sensor technology allow monitoring of soil moisture dynamics with high temporal resolution at varying spatial scales. The objectives of this study were to: (i) develop an efficient strategy for monitoring soil moisture dynamics at the hillslope scale using a wireless sensor network; and (ii) characterize spatial patterns of soil moisture and infer hydrological processes controlling the dynamics of such patterns, using a method of analysis that allows the identification of the relevant hydrological dynamics within large data sets. We combined soil hydrological and pedological expertise with geophysical measurements and methods from digital soil mapping for designing the monitoring setup for a grassland hillslope in the Schäfertal catchment, central Germany. Hypothesizing a wet and a dry soil moisture state to be characteristic of the spatial pattern of soil moisture, we described the spatial and temporal evolution of such patterns using a method of analysis based on the Spearman rank correlation coefficient. We described the persistence and switching mechanisms of the two characteristic states, inferring the local properties that control the observed spatial patterns and the hydrological processes driving the transitions. The spatial organization of soil moisture appears to be controlled by different processes in different soil horizons, and the topsoil's moisture does not mirror processes that take place within the soil profile. The results will help to improve conceptual understanding for hydrological model studies at similar or smaller scales and to transfer observation concepts and process understanding to larger or less instrumented areas.

Stable water isotope tracing through hydrological models for disentangling runoff generation processes at the hillslope scale

Abstract. Hillslopes are the dominant landscape components where incoming precipitation becomes groundwater, streamflow or atmospheric water vapor. However, directly observing flux partitioning in the soil is almost impossible. Hydrological hillslope models are therefore being used to investigate the processes involved. Here we report on a modeling experiment using the Catchment Modeling Framework (CMF) where measured stable water isotopes in vertical soil profiles along a tropical mountainous grassland hillslope transect are traced through the model to resolve potential mixing processes. CMF simulates advective transport of stable water isotopes 18O and 2H based on the Richards equation within a fully distributed 2-D representation of the hillslope. The model successfully replicates the observed temporal pattern of soil water isotope profiles (R2 0.84 and Nash–Sutcliffe efficiency (NSE) 0.42). Predicted flows are in good agreement with previous studies. We highlight the importance of groundwater recharge and shallow lateral subsurface flow, accounting for 50 and 16% of the total flow leaving the system, respectively. Surface runoff is negligible despite the steep slopes in the Ecuadorian study region.

Interactions and connectivity between runoff generation processes of different spatial scales

AbstractMonitoring runoff generation processes in the field is a prerequisite for developing conceptual hydrological models and theories. At the same time, our perception of hydrological processes strongly depends on the spatial and temporal scale of observation. Therefore, the aim of this study is to investigate interactions between runoff generation processes of different spatial scales (plot scale, hillslope scale, and headwater scale). Different runoff generation processes of three hillslopes with similar topography, geology and soil properties, but differences in vegetation cover (grassland, coniferous forest, and mixed forest) within a small v‐shaped headwater were measured: water table dynamics in wells with high spatial and temporal resolution, subsurface flow (SSF) of three 10 m wide trenches at the bottom of the hillslopes subdivided into two trench sections each, overland flow at the plot scale, and catchment runoff. Bachmair et al. () found a high spatial variability of water table dynamics at the plot scale. In this study, we investigate the representativity of SSF observations at the plot scale versus the hillslope scale and vice versa, and the linkage between hillslope dynamics (SSF and overland flow) and streamflow. Distinct differences in total SSF within each 10 m wide trench confirm the high spatial variability of the water table dynamics. The representativity of plot scale observations for hillslope scale SSF strongly depends on whether or not wells capture spatially variable flowpaths. At the grassland hillslope, subsurface flowpaths are not captured by our relatively densely spaced wells (3 m), despite a similar trench flow response to the coniferous forest hillslope. Regarding the linkage between hillslope dynamics and catchment runoff, we found an intermediate to high correlation between streamflow and hillslope hydrological dynamics (trench flow and overland flow), which highlights the importance of hillslope processes in this small watershed. Although the total contribution of SSF to total event catchment runoff is rather small, the contribution during peak flow is moderate to substantial. Additionally, there is process synchronicity between spatially discontiguous measurement points across scales, potentially indicating subsurface flowpath connectivity. Our findings stress the need for (i) a combination of observations at different spatial scales, and (ii) a consideration of the high spatial variability of SSF at the plot and hillslope scale when designing monitoring networks and assessing hydrological connectivity. Copyright © 2013 John Wiley & Sons, Ltd.

Hillslope characteristics as controls of subsurface flow variability

Abstract. Hillslope hydrological dynamics, particularly subsurface flow (SSF), are highly variable and complex. A profound understanding of factors controlling this variability is needed. Therefore we investigated the relationship between variability of shallow water table dynamics and various hillslope characteristics. We ask whether measurable hillslope properties explain patterns of subsurface flow variability. To approach this question, shallow water table dynamics of three adjacent large-scale hillslopes were monitored with high spatial and temporal resolution over 18 months. The hillslopes are similar in terms of topography and parent material, but different in vegetation cover (grassland, coniferous forest, and mixed forest). We expect vegetation to be an important driver of water table dynamics at our study site, especially given the minor differences in topography. Various hillslope properties were determined in the field and via GIS analysis: common topography descriptors, well depth, soil properties via slug tests, and several vegetation parameters. Response variables characterizing the water table response per well were calculated for different temporal scales (entire time series, seasonal scale, event scale). Partial correlation analysis and a Random Forest machine learning approach were carried out to assess the explainability of SSF variability by measurable hillslope characteristics. We found a complex interplay of predictors, yet soil properties and topography showed the highest single explanatory power. Surprisingly, vegetation characteristics played a minor role. Solely throughfall and canopy cover exerted a slightly stronger control, especially in summer. Most importantly, the examined hillslope characteristics explained only a small proportion of the observed SSF variability. Consequently there must be additional important drivers not represented by current measurement techniques of the hillslope configuration (e.g. bedrock properties, preferential pathways). We also found interesting differences in explainability of SSF variability among temporal scales and between both forested hillslopes and the grassland hillslope.

Dynamics of critical source areas: Does connectivity explain chemistry?

Critical source area approaches to catchment management are increasingly being recognised as effective tools to mitigate sediment and nutrient transfers. These approaches often assume hydrological connectivity as a driver for environmental risk, however this assumption has rarely been tested. Using high resolution monitoring, 14 rainfall events of contrasting intensity were examined in detail for spatial and temporal dynamics of overland flow generation at a hydrologically isolated grassland hillslope in Co. Down, Northern Ireland. Interactions between overland flow connectivity and nutrient transfers were studied to test the critical source area hypothesis. While total and soluble phosphorus loads were found to be representative of the size of the overland flow contributing area (P=<0.05), the dynamics of concentrations throughout storm hydrographs were found to be complex and storm dependant. Near linear relationships were observed between the contributing area and total overland flow volumes (R2=0.86). Export coefficients (kgha−1) calculated using plot size were found to under estimate annual losses of total phosphorus by a factor of 17, when compared to those calculated using the contributing area. This study shows that current critical source area definitions for implementing mitigation measures may be overlooking the importance of storm characteristics in determining nutrient transfers and hence may be insufficient in determining catchment scale risk.

Intercomparing hillslope hydrological dynamics: Spatio‐temporal variability and vegetation cover effects

Generalizable process knowledge on hillslope hydrological dynamics is still very poor, yet indispensable for numerous theoretical and practical applications. To gain insight into the organization of hillslope hydrological dynamics we intercompared 90 observations of shallow water table dynamics at three neighboring large‐scale (33 × 75 m) hillslopes with similar slope, aspect, curvature, geologic, and pedologic properties but differences in vegetation cover (grassland, coniferous forest, and mixed forest) over a time period of 9 months. High‐resolution measurements of water table fluctuations, rainfall, and discharge in the creek at the foot of all hillslopes allowed a good system characterization. The aim of this study was to explore the spatio‐temporal variability of water table fluctuations within and between hillslopes, the effect of event and antecedent characteristics on the observed dynamics, and how the hillslope subsurface flow (SSF) response is reflected in the runoff response. To intercompare the SSF behavior we conducted an event‐based analysis of the percentage of well activation, several metrics characterizing the shape and timing of the water table response curves, rainfall characteristics, antecedent wetness conditions, and several runoff response metrics. The analysis reveals that there are distinct differences in SSF response between the grassland hillslope and the forested hillslopes, with a lower frequency of well activation and absolute water table rise at the grassland hillslope. Second, spatial patterns of water table dynamics differ between wet fall/winter/spring (predominantly saturation of the lower part of the hillslope, weaker water table response, and slower response times) and dry summer conditions (whole‐hillslope activation but higher spatial variability, generally stronger water table dynamics, and quicker response times). The observed seasonally changing water table dynamics suggest the development of a preferential flow network during high‐intensity rainstorms under dry summer conditions. Third, catchment runoff is strongly driven by hillslope dynamics, yet contrasting hydrographs during events with similar hillslope dynamics indicate the influence of additional processes. Overall, the observed high spatio‐temporal variability of seemingly homogeneous hillslopes calls for rethinking of current monitoring strategies and developing and testing new conceptual models of hillslope hydrologic processes.

Preferential Flow Effects on Infiltration and Runoff in Grassland and Forest Soils

Due to the heterogeneity of soil hydraulic properties, there are many ways in which natural soils respond to rainfall, making upscaling of flow processes from plot to catchment scale difficult. The objectives of this study were to qualitatively characterize the flow pathways on forest and grassland hillslopes and to quantitatively define the relevant parameters controlling surface runoff generation by using infiltration and dye tracer experiments supplemented with measurements of the saturated hydraulic conductivity Ksat, the structural porosity nSP, and the pore‐size distribution PSD. While infiltration excess overland flow dominates in grassland, forest soil structure characterized by relatively high values of Ksat and nSP enhances the infiltrability of the soil and consequently prevents or at least reduces surface runoff. The dye patterns suggest that macropores are more efficient in forest than in grassland soil. The low efficiency of grassland soil macropores in transporting all water vertically downward can be explained by (i) the fine and dense few topsoil layers caused by the land use that limit water flux into the underlying macropores and (ii) their restricted number, their tortuosity, and the restricted interaction between macropores and the matrix below the topsoil layer. The larger root water uptake of forest soil as compared to grassland soil can be viewed as an additional factor enhancing its storage capacity and, consequently, may reduce the generation of surface runoff. It remains unclear, however, what effect the low interaction between macropores and soil matrix in the upper part of the subsoil has on surface runoff in grassland soil; this should be investigated in future studies.

Using sediment travel distance to estimate medium‐term erosion rates: a 16‐year record

AbstractFine grained (80 µm) magnetite was introduced onto a semi‐arid grassland hillslope in 1992, as part of a set of rainfall‐simulation experiments. Using measurements of magnetic susceptibility, the median distance travelled by these magnetite grains during subsequent natural runoff events in the 16‐year period up to 2008 was estimated. Coupling this estimate to direct measurements of sediment flux obtained during the rainfall‐simulation experiments has enabled estimation of the erosion rate over this period. The estimated average erosion rate of between 2·61 × 10−2 and 4·36 × 10−2 kg m−1 year−1, is equivalent to a rate of ground lowering between 0·020 and 0·033 mm year−1. This estimate is consistent with (in the sense of being less than) an estimate of total sediment detachment over the same period. The rate of erosion measured using this travel‐distance approach is an order of magnitude less that obtained from a study based on 137Cs in a nearby catchment, and compatible with the longevity of continents. Copyright © 2010 John Wiley &amp; Sons, Ltd.

Overland flow initiation from a drumlin grassland hillslope

AbstractAccurate prediction of overland flow requires an understanding of the hydrological processes controlling its occurrence over a range of spatial and temporal scales. This study investigated the factors controlling the initiation of overland flow from grassland on a drumlin. From January 2003 to December 2006, fine scale monitoring of soil moisture, rainfall and overland flow were carried out at a grazed grassland site in Northern Ireland. Logistic regression analyses were used to identify the factors controlling the initiation of overland flow. Relationships between the observed volumetric soil moisture (VSM) and soil moisture deficit (SMD) values predicted from regional meteorological data were also compared. Results demonstrated that although saturation excess overland flow occurs at this site, 59% of overland flow events occurred on days when soil moisture was below field capacity. The logistic regression analysis confirmed that when soils were below field capacity, average rainfall intensity was the most important variable in determining the probability of overland flow occurring followed by the SMD. Although VSM values at this site were correlated with two sets of modelled values of SMD, one based on the meteorological office rainfall and evaporation calculation system model of the UK Meteorological Office and the other based on a study from Ireland, neither model provided accurate indicators of the risk of overland flow. Each significantly overestimated the number of days when soils were at or above field capacity. This uncertainty and the predicted increase in high‐intensity rainfall events as a result of climate change pose challenges to the use of SMD as an indicator of the risk of overland flow.

The characteristics of soil water cycle and water balance on steep grassland under natural and simulated rainfall conditions in the Loess Plateau of China

Summary Large-scale vegetation restoration has been helpful to prevent serious soil erosion, but also has aggravated water scarcity and resulted in soil desiccation below a depth of 200 cm in the Loess Plateau of China. To understand the dynamic mechanism of soil desiccation formation is very important for sustainable development of agriculture in the Loess Plateau. Based on natural and simulated rainfall, the characteristics of soil water cycle and water balance in the 0–400 cm soil layer on a steep grassland hillslope in Changwu County of Shaanxi Loess Plateau were investigated from June to November in 2002, a drought year with annual rainfall of 460 mm. It was similarly considered to represent a rainy year with annual rainfall of 850 mm under simulated rainfall conditions. The results showed that the temporal variability of water contents would decrease in the upper 0–200 cm soil layer with the increase in rainfall. The depth of soil affected by rainfall infiltration was 0–200 cm in the drought year and 0–300 cm in the rainy year. During the period of water consumption under natural conditions, the deepest layer of soil influenced by evapotranspiration (ET) rapidly reached a depth of 200 cm on July 21, 2002, and soil water storage decreased by 48 mm from the whole 0–200 cm soil layer. However, during the same investigation period under simulated rainfall conditions, soil water storage in the 0–400 cm soil layer increased by only 71 mm, although the corresponding rainfall was about 640 mm. The extra-simulated rainfall of 458 mm from May 29 to August 10 did not result in the disappearance of soil desiccation in the 200–400 cm deep soil layer. Most infiltrated rainwater retained in the top 0–200 cm soil layer, and it was subsequently depleted by ET in the rainy season. Because very little water moved below the 200 cm depth, there was desiccation in the deep soil layer in drought and normal rainfall years.

Rapid immobilisation and leaching of wet-deposited nitrate in upland organic soils

Nitrate (NO3-) is often observed in surface waters draining terrestrial ecosystems that remain strongly nitrogen (N) limited. It has been suggested that this occurs due to hydrological bypassing of soil or vegetation N retention, particularly during high flows. To test this hypothesis, artificial rain events were applied to 12 replicate soil blocks on a Welsh podzolic acid grassland hillslope, labelled with 15N-enriched NO3- and a conservative bromide (Br-) tracer. On average, 31% of tracer-labelled water was recovered within 4 h, mostly as mineral horizon lateral flow, indicating rapid vertical water transfer through the organic horizon via preferential flowpaths. However, on average only 6% of 15N-labelled NO3- was recovered. Around 80% of added NO3- was thus rapidly immobilised, probably by microbial communities present on the surfaces of preferential flowpaths. Transitory exceedance of microbial N-uptake capacity during periods of high water and N flux may therefore provide a mechanism for NO3- leaching.

Depth distribution of interrill overland flow and the formation of rills

AbstractA Gumbel distribution for maxima is proposed as a model for the depths of interrill overland flow. The model is tested against three sets of field measurements of interrill overland flow depths obtained on shrubland and grassland hillslopes at Walnut Gulch Experimental Watershed, southern Arizona. The model is found to be a satisfactory fit to 81 of the 90 measured distributions. The shape δ and location λ parameters of all fitted distributions are strongly correlated with discharge. However, whereas a common relationship exists between discharge and δ for all depth distributions, the relationships with λ vary systematically downslope. Using the Gumbel distribution as a model for the distribution of overland flow depths, a probabilistic model for the initiation of rills is developed, drawing upon the previous work of Nearing. As an illustration of this approach, we apply this model to the shrubland and grassland hillslopes at Walnut Gulch. It is concluded that the presence of rills on the shrubland, but not on the grassland, is due to the greater runoff coefficient for the shrubland and/or the greater propensity of the shrubland for soil disturbance compared with the grassland. Finally, a generalized conceptual model for rill initiation is proposed. This model takes account of the depth distribution of overland flow, the probability of flow shear stress in excess of local soil shear strength, the spatial variability in soil shear strength and the diffusive effect of soil detachment by raindrops. Copyright © 2005 John Wiley &amp; Sons, Ltd.

Prediction of sediment‐bound nutrient delivery from semi‐arid California watersheds

Soil carbon (C), nitrogen (N), and phosphorus (P) are lost from hillslopes in particulate forms through soil erosion. The fate of the eroded C (e.g., sequestration or oxidation) may affect the global C budget, and delivery of N and P to waterbodies can lead to eutrophication. Whereas the magnitude of particulate nutrient losses may be similar to or greater than dissolved losses, it is rarely estimated. We couple a sediment delivery model with measurements of C, N, and P in soil to account explicitly for hillslope sediment transport processes that yield sediment‐bound nutrients to fluvial networks. The model is applied to a site in California dominated by coastal sage scrub and gopher‐rich grasslands. Although the magnitude of sediment delivery predicted by the model has been tested with reservoir sedimentation records, no data exist to test the predicted rates of nutrient delivery. Nevertheless, the model results are provocative; it predicts that losses of particulate C from sage covered hillslopes (23 kg/ha/yr) are nearly double that from grassland hillslopes (13 kg/ha/yr), despite a lower annual sediment yield from the sage hillslopes. The model predicts similar average annual N and P losses for sage and grasslands but dramatic differences in the frequency and magnitude of delivery events. Nutrient delivery from grasslands is chronic whereas delivery from the coastal sage is highly episodic, with large pulses driven by fire frequency. These results suggest that changes in the vegetation community can alter the delivery regime of sediment‐bound C, N, and P.

Production of CO2 in Soil Profiles of a California Annual Grassland

Soils play a key role in the global cycling of carbon (C), storing organic C, and releasing CO2 to the atmosphere. Although a large number of studies have focused on the CO2 flux at the soil–air interface, relatively few studies have examined the rates of CO2 production in individual layers of a soil profile. Deeper soil horizons often have high concentrations of CO2 in the soil air, but the sources of this CO2 and the spatiotemporal dynamics of CO2 production throughout the soil profile are poorly understood. We studied CO2 dynamics in six soil profiles arrayed across a grassland hillslope in coastal southern California. Gas probes were installed in each profile and gas samples were collected weekly or biweekly over a three-year period. Using soil air CO2 concentration data and a model based on Fick’s law of diffusion, we modeled the rates of CO2 production with soil profile depth. The CO2 diffusion constants were checked for accuracy using measured soil air 222Rn activities. The modeled net CO2 production rates were compared with CO2 fluxes measured at the soil surface. In general, the modeled and measured net CO2 fluxes were very similar although the model consistently underestimated CO2 production rates in the surficial soil horizons when the soils were moist. Profile CO2 production rates were strongly affected by the inter- and intra-annual variability in rainfall; rates were generally 2–10 times higher in the wet season (December to May) than in the dry season (June to November). The El Niño event of 1997–1998, which brought above-average levels of rainfall to the study site, significantly increased CO2 production in both the surface and subsurface soil horizons. Whole profile CO2 production rates were approximately three times higher during the El Niño year than in the following years of near-average rainfall. During the dry season, when the net rates of CO2 flux from the soil profiles are relatively low (4–11 mg C– CO2 m−2 h−1), 20%–50% of the CO2 diffusing out of the profiles appears to originate in the relatively moist soil subsurface (defined here as those horizons below 40 cm in depth). The natural abundance 14C signatures of the CO2 and soil organic C suggest that the subsurface CO2 is derived from the microbial mineralization of recent organic C, possibly dissolved organic C transported to the subsurface horizons during the wet season.

Landslides on coastal sage-scrub and grassland hillslopes in a severe El Niño winter: The effects of vegetation conversion on sediment delivery

During the 1997-1998 El Nino, record rainfall triggered.150 shallow landslides within a 9.5 km2 area near Santa Barbara, Cal- ifornia. They were studied to analyze the sediment delivery to val- ley floors from landslides in coastal sage scrub and converted grasslands. The conversion of coastal sage to grasslands, primarily to provide pasturage for cattle, is common in the region, and the landscape's response may affect water quality, reservoir infilling, and debris flow hazards. We explore the relationship between lateral- root reinforcement and landslide volume by developing a slope- stability analysis that incorporates root cohesion along the sides of the failure. The stability analysis correctly predicts an inverse re- lationship between landslide volume and hillslope angle in the sage. The volumes of failures in the grasslands do not vary systemati- cally with slope and are generally smaller than those in the sage. From aerial-photograph analysis and field mapping, we find that there are 22.9 failures per square kilometer in the grasslands com- pared to 13.2 failures per square kilometer in the sage. Despite the lower failure density in the coastal sage, greater failure volumes and longer transport distances delivered more sediment to valley floors, with a specific volumetric flux of 2.8 3 10 -2 m 3 ·m -1 for this El Nino compared to 1.73 10 -2 m 3 ·m -1 in the grasslands. We con- clude that the conversion from vegetation with stronger and deeper roots (coastal sage) to vegetation with weaker and shallower roots (grass) has caused a pulse of increased landsliding in the grasslands because the soils are currently too thick for the prevailing root reinforcement. We suggest that, over time, soils in the grassland hollows will become thinner as the evacuation by landslides is re- peated until the landsliding rate declines to balance the soil sup- plied from local colluvium production and diffusive processes upslope.