Using hysteretic behaviour and hydrograph classification to identify hydrological function across the “hillslope–depression–stream” continuum in a karst catchment

In cockpit karst landscapes, fluxes from upland areas contribute large volumes of water to low-lying depressions and stream flow. Hydrograph hysteresis and similarity between monitoring sites is important for understanding the space–time variability of hydrologic responses across the “hillslope–depression–stream” continuum. In this study, the hysteretic feature of hydrographs was assessed by characterizing the loop-like relationships between responses at upstream sites relative to subsurface discharge at the outlet of a small karst catchment. A classification of hydrograph responses based on the multi-scale smoothing Kernel -derived distance classifies the hydrograph responses on the basis of similarities between hillslope and depression sites, and those at the catchment outlet. Results demonstrate that the temporal and spatial variability of hydrograph hysteresis and similarity between hillslope flow and outlet stream flow can be explained by the local heterogeneity of depression aquifer. Large depression storage deficits emerging in the highly heterogeneous aquifer produce strong hysteresis and multiple relationships of upstream hydrographs relative to the outlet subsurface discharge. In contrast, when depression storage deficits are filled during consecutive rainfall events, depression hydrographs at the high permeability sites are almost synchronous or exhibit a monotonous function with the hydrographs at the outlet. This reduced hydrograph hysteresis enhances preferential flow paths in fractured rocks and conduits that can accelerate the hillslope flow to the outlet. Therefore, classification of hydrograph similarities between any upstream sites and the catchment outlet can help to identify the dominant hydrological functions in the heterogeneous karst catchment.     

The generation and redistribution of overland flow on a massive oxic soil in a eucalypt woodland within the semi-arid tropics of North Australia

There is a dearth of knowledge on the runoff processes of eucalypt woodland communities in the semi-arid tropics of Australia. The work was undertaken on a 100 m transect of a 0·8 degree hillslope typical of the ‘smooth plainlands’ of central-north Queensland. This paper introduces a new experimental design for measuring overland flow in such areas by way of a cascade system of unbounded runoff plots which allow the inputs and outputs between troughs to be calculated. Most storms generate overland flow. Time to overland flow ranges between 1 and 18 min where rain intensities are above 10mm hr⁻¹ and when the average detention storage of 3·6 mm is exceeded. The bare soil surfaces within the scattered grass understory control the runoff generation process through the temporal variability of field saturated hydraulic conductivity. The study demonstrated that overland flow is mainly redistributed over the freely-draining oxic soil. Some areas export more overland flow than they gain from upslope (runoff), others gain more overland flow than they export (runon). Over the study period only 2 per cent of total rain is transferred out of this 100 m transect as overland flow due to the short duration of storms, the relatively high soil permeability, and the low slope angle. The remainder adds to the large soil water store or deep drainage. The variability of runoff–runon over these ‘smooth plainlands’ highlights how results from bounded plots would be misleading in such areas.     

PLAN FORM AND PROFILE FORM OF HILLSLOPES

For 27 hillslope profiles located within four first-order drainage basins on the dip slope of the South Downs in East Sussex, contour curvature (P), measured in degrees per 100 m, is highly variable within any one hillslope. Hillslopes may be described as wholly concave in plan where P is less than ?40, wholly convex in plan where P is greater than +40 and convexo-concave in plan where P lies between ?40 and +40. One measure of P taken at the profile sampling line is found to give a better estimate of P than measuring P at the steepest point on the hillslope. P is found to be significantly negatively correlated with slope in 38 per cent of cases, correlated with distance from the divide in 55 per cent of cases (positively where P <0 and negatively where P <0) and positively correlated with profile curvature in 70 per cent of cases. Implications of these correlations for hillslope process-response models are discussed.     

Infiltration into frozen soil: From core‐scale dynamics to hillslope‐scale connectivity

Infiltration into frozen soil is a key hydrological process in cold regions. Although the mechanisms behind point‐scale infiltration into frozen soil are relatively well understood, questions remain about upscaling point‐scale results to estimate hillslope‐scale run‐off generation. Here, we tackle this question by combining laboratory, field, and modelling experiments. Six large (0.30‐m diameter by 0.35‐m deep) soil cores were extracted from an experimental hillslope on the Canadian Prairies. In the laboratory, we measured run‐off and infiltration rates of the cores for two antecedent moisture conditions under snowmelt rates and diurnal freeze–thaw conditions observed on the same hillslope. We combined the infiltration data with spatially variable data from the hillslope, to parameterise a surface run‐off redistribution model. We used the model to determine how spatial patterns of soil water content, snowpack water equivalent (SWE), and snowmelt rates affect the spatial variability of infiltration and hydrological connectivity over frozen soil. Our experiments showed that antecedent moisture conditions of the frozen soil affected infiltration rates by limiting the initial soil storage capacity and infiltration front penetration depth. However, shallow depths of infiltration and refreezing created saturated conditions at the surface for dry and wet antecedent conditions, resulting in similar final infiltration rates (0.3 mm hr^?1). On the hillslope‐scale, the spatial variability of snowmelt rates controlled the development of hydrological connectivity during the 2014 spring melt, whereas SWE and antecedent soil moisture were unimportant. Geostatistical analysis showed that this was because SWE variability and antecedent moisture variability occurred at distances shorter than that of topographic variability, whereas melt variability occurred at distances longer than that of topographic variability. The importance of spatial controls will shift for differing locations and winter conditions. Overall, our results suggest that run‐off connectivity is determined by (a) a pre‐fill phase, during which a thin surface soil layer wets up, refreezes, and saturates, before infiltration excess run‐off is generated and (b) a subsequent fill‐and‐spill phase on the surface that drives hillslope‐scale run‐off.     

Interactions between perched and saprolite aquifers on a small,salt-affected and deeply weathered hillslope

A small hillslope was chosen to investigate the role of throughflow as a mechanism responsible for the movement of soil water and solutes towards a saline seep and as a source of recharge to a permanent, regional aquifer at depth. The hydraulic properties, chemical characteristics and physical responses of both systems were studied on a deeply weathered, salt-affected hillslope. Additional data were also obtained from other sites in south-western Australia. Regional groundwater flow occurred in a variably textured, deeply weathered material in which the hydraulic conductivity varied from < 0·001 to 0·14m day^?1. Perched groundwater flow (throughflow) occurred in the higher permeability (? 1·5 m day^?1), near-surface soil materials. Throughflow occurred throughout winter, contributing approximately 530 m³ of fresh (? 160 mg l^?1 Cl) water to a saline seep. By contrast, the deep aquifer discharged approximately 1100 m³ of waters with salt concentrations of 2000–6000 mg l^?1 Cl. Recharge and discharge rates to and from the deep aquifer, were estimated to be of the order of 5–20 mm a^?1 and 50–300 mm a^?1 respectively. Saturated conditions existed throughout winter within the seep and the immediately adjacent non-saline area, with up to 60 per cent of the hillslope soils becoming saturated after major rainfall events ( > 20 mm day^?1). In the mid-slopes, in particular along a central depression, saturation of the shallow soils caused macropore channel recharge to take place through the clay-textured subsoils. Water-level responses suggest that approximately 25–30 per cent of annual recharge occurred from one storm studied in September 1984. Recharge through macropore channels is a significant mechanism in the concave slope areas on the hillslope. Throughflow was found to be a major source of water, but not salt, contributing to the saline seep. In general, the contribution of throughflow was found to decrease further inland at other sites studied. However, at inland sites where perennial, perched aquifers have developed in deep sands, saline areas have been caused by throughflow and not by deep aquifer discharge.     

Seasonal soil water dynamics in the jarrah forest,Western Australia. II: Results from a site with fine-textured soil profiles

Seasonal soil water dynamics were measured at a fine-textured, upslope site within the jarrah forest of southwest Western Australia and compared to the results from a coarse-textured hillslope transect. Gravity drainage dominated during winter and early spring. This reversed in early summer and an upward potential gradient was observed to 7 m depth. A shallow ephemeral saturation zone was observed above a clay pan at 1.5 m depth. This saturation zone persisted through late winter and early spring, contrasting with the short-lived saturation in the duricrust on the hillslope transect. The annual maximum to minimum unsaturated soil water storage was about 530 mm, 50 mm greater than the hillslope transect and higher than most values reported elsewhere in Australia. Significant soil water content changes following winter rain were generally restricted to 6 m but at one site occurred to 9 m. These depths were significantly less than the coarser-textured hillslope transect. Soil water drying rates averaged 5 mm day^?1 during extended dry periods compared to 3.5 mm day^?1 on the hillslope transect. The drying rate occurred uniformly through the profile until late summer when a significant decrease in the upper 3 m was observed.     

The quantity of stormwater runoff from ten stretches of road,a car park and eight roofs in Hertfordshire,England during 1983

Rainfall and runoff were monitored simultaneously for one year from a residential road, a car park, nine sections of road draining to individual gullies, two house roofs, two garage roofs, and three types of factory roof. The sites, which included an automatic weather station, were in Redbourn, Hertfordshire on Flood Studies Report Soil Type 1. The 2906 quality controlled ‘station-storms’ represented 193 rain storms and involved 57.2 per cent of the annual rainfall. 1732 storms were of less than 1.4mm of rain, whilst 77 had over 10mm. The percentage runoff averaged 11.4 per cent for roads and 56.9 per cent for roofs (28.3 per cent and 90.4 per cent respectively for rainfalls >5mm). Percentage runoff from the roads was cyclic with a peak during the summer months but there was a marked variation in monthly percentage runoff within and between sites. Regression analysis to explain percentage runoff was undertaken with various subsets of data for: each site; roads; and roofs. The regression analysis considered all storms; >1 percent runoff events; >5mm rainfalls; and events with > = 4 mm rain and > = 5 per cent runoff. The variable values in percentage runoff could not be explained satisfactorily with statistical methods. Only eight of the 72 equations explained more than 57 per cent of the variance. The most important explanatory variables for roads were short term rainfall intensity and rainfall amount, the former was the most important for roofs. ‘Seasonal’ variables had a positive relation ship for roads which shows that the percentage runoff from roads is higher in summer than winter. The antecedent variables showed that percentage runoff from roads and roofs is increased by antecedent rainfall. Seasonal factors and evaporation were unimportant for the percentage runoff from roofs. Depression storage, assessed by examining rainfalls that did and did not produce runoff, showed a diversity of monthly values. The depression storages derived by the regression intercept method were usually smaller. There were no relationships between depression storage and catchment or roof slope. The mean values for roofs and roads respectively were 0.52 mm and 1.23 mm for the classification method and 0.42 mm and 0.6mm with the regression approach. Peak runoff from the roads showed an attenuation to 12.8 per cent for 1 minute rainfall intensities and 24.2 per cent for 5 minute intensities. For roofs the attenuation averaged 36.8 per cent for 1 minute intensities and 92.6 for 5 minute intensities. Regression for peak runoff coefficients from roofs and roads explained negligible amounts of the variance except when events with 1 minute rainfall intensities of over 30 mm hr^?1 over the roads were analysed. Total rainfall was an important explanatory variable as was the slope of the road. There was evidence that peak coefficients for roads are greater during the summer.     

Hydrogeological response to climate change in alpine hillslopes

Climate change threatens water resources in snowmelt‐dependent regions by altering the fraction of snow and rain and spurring an earlier snowmelt season. The bulk of hydrological research has focused on forecasting response in streamflow volumes and timing to a shrinking snowpack; however, the degree to which subsurface storage offsets the loss of snow storage in various alpine geologic settings, i.e. the hydrogeologic buffering capacity, is still largely unknown. We address this research need by assessing the affects of climate change on storage and runoff generation for two distinct hydrogeologic settings present in alpine systems: a low storage granitic and a greater storage volcanic hillslope. We use a physically based integrated hydrologic model fully coupled to a land surface model to run a base scenario and then three progressive warming scenarios, and account for the shifts in each component of the water budget. For hillslopes with greater water retention, the larger storage volcanic hillslope buffered streamflow volumes and timing, but at the cost of greater reductions in groundwater storage relative to the low storage granite hillslope. We found that the results were highly sensitive to the unsaturated zone retention parameters, which in the case of alpine systems can be a mix of matrix or fracture flow. The presence of fractures and thus less retention in the unsaturated zone significantly decreased the reduction in recharge and runoff for the volcanic hillslope in climate warming scenarios. This approach highlights the importance of incorporating physically based subsurface flow in to alpine hydrology models, and our findings provide ways forward to arrive at a conceptual model that is both consistent with geology and hydrologic principles. Copyright © 2016 John Wiley & Sons, Ltd.     

Ten year remeasurement of chemical denudation on a magnesian limestone hillslope

Results from the measurement of microweight loss of Magnesian Limestone rock tablets placed at the soil–bedrock interface on a hillslope over 10 years (1982–1992) gave the same relative pattern of upslope increase in weight loss as did short-term measurements (1979–1981), but the absolute values were an order of magnitude lower (0·01–0·03 per cent per year long term and 0·1–0·6 per cent per year short term). Microweighed rock tablets may therefore give unreliable absolute rates and can only be used to indicate relative spatial differences rather than to give reliable data on temporal changes.     

Seasonal soil water dynamics in the jarrah forest,Western Australia. I: Results from a hillslope transect with coarse-textured soil profiles

Seasonal soil water dynamics were measured on a hillslope transect in the jarrah forest of southwest Western Australia over the period 1984-86 using mercury manometer tensiometers, gypsum blocks, and a neutron moisture meter. The soil water potential gradients indicated downward vertical drainage flux through winter and spring. There was generally a change to an upwards flux in early summer which was sustained through to autumn. A shallow ephemeral saturation zone was identified in and above a duricrust layer, lasting up to three days after heavy, late winter rainfall. The annual maximum to minimum unsaturated soil water storage on the hillslope was approximately 400 mm to 6 m depth and 480 mm to 15 m depth. This did not change significantly in years of substantially different winter rainfall. The magnitude of seasonal soil water storage was similar to other forested areas with deep soil profiles. The depth of observable infiltration was dependent on annual rainfall. This was consistent with the observation that groundwater levels responded to rainfall over the whole hillslope in wet years but only responded on the lower slopes in dry years. The average summer drying rate of the soil profile to 6 m depth of 3.5 mm day^?1 was within the range of values reported for forests elsewhere. In late summer, following an extended drought period, the drying rate decreased downslope but increased midslope.     

EROSION PROCESSES AT THE HILLSLOPE SCALE

The results are presented of an intensive monitoring programme to determine eroded volumes due to observed rainfall events on a hillslope surface. The surface investigated has been reproduced through a digital elevation model and analysed in terms of drainage network, contributing area and slope by performing planimetric and altimetric analysis. Erosion maps were derived from a comparison of consecutive digital elevation models relative to time. These maps reveal the spatial and temporal evolution of the erosion process at the hillslope scale. The erosion was uniform across the surface, supporting the assumption of randomness in the erosion process commonly used in surface drainage development models. The observed value of erosion has been estimated at approximately 0.11 m/m² per year, with almost 500 mm of total annual rainfall.     

Sediment budgeting: A case study in the Katiorin drainage basin,Kenya

The primary objective of this study was to compute a detailed budget for a small semiarid tropical drainage basin in Kenya. Results indicated that transfer of sediments (‘inputs’) from primary source areas was minor in comparison to changes in storage. The major sediment source area within the Katiorin drainage basin was the colluvial hillslope zone. The net change in storage within this zone was approximately 2100 Mg yr^?1. Surface wash and rilling were the dominant transport processes responsible for the remobilization of colluvial sediments. Sediment storage within the in-channel reservoir increased by 60 Mg yr^?1, which was minor when compared to the total store of sediment in this reservoir. During 1986, the channel network stored only a small fraction ( < 3 per cent) of the sediment delivered from the hillslope subsystem. Therefore, the in-channel reservoir had limited influence on sediment conveyance to the basin outlet. These data indicate that a static equilibrium condition cannot be assumed within the Katiorin drainage basin. Such an assumption would result in erosion estimates of approximately 5.5 mm yr^?1 for the entire basin (based on a sediment output of 7430 Mg km^?2 yr^?1 and a measured bulk density of 1.35 Mg m^?3). However, this masked the actual rates of 1.2 to 7.1 mm yr^?1 in subbasin primary source areas, and rates of 0.6 to 17 mm yr^?1 for colluvial material in the various subbasins. The extreme accelerated erosion rates resulted from minimal ground vegetation, steep slopes, soil crust formation, an erodible substrate, and a well-integrated drainage network for rapid conveyance of sediments from the hillslope subsystem to the basin outlet.     

III: Methods of correcting for systematic errors in atmospheric precipitation measurements in Czechoslovakia

The methods of correcting for systematic errors in precipitation measurements using the Czechoslovak gauge METRA 886 are presented. This gauge has an orifice area of 500 cm² and is elevated 1 m above the ground. The wetting correction amounts to 0.1–0.2 mm per measurement. The evaporation correction ranges from 3 per cent in frost-free periods to 10 per cent in frost periods. The wind-induced correction amounts to 5 per cent for rain and 45 per cent for snow. The total sum of corrections on average 10–15 per cent per year in lower localities.     

Distribution of organic carbon in bed sediments of Manoa Stream,Oahu, Hawaii

Organic carbon (OC) associated with fluvial bed sediment plays an important role in biotic and abiotic processes operating within drainage basins. Increasingly, there is a need to characterize storage and spatial distributions of OC in aquatic sediments, particularly under-sampled areas like tropical streams. The objectives of this study were to examine in detail the variation of OC concentration with bed sediment grain size, to characterize the influence of grain size variation on relative OC mass storage, and to compare weighted OC values to those in other aquatic sediments worldwide. The study area selected was a third-order dendritic drainage basin developed in a basaltic complex. Bed sediments along a 6 km section of Manoa Stream were systematically sampled every 50 m for a total of 113 sample site locations. Sediments were partitioned into six size fractions (< 2·0 mm) and OC was determined by dry combustion. Data indicate that the OC concentration increases with decreasing grain size, with the greatest values in the < 0·063 mm (silt + clay) fraction, approximately 4·6 times greater than the very coarse sand fraction (1·00–2·00 mm). Robust smoothing techniques illustrated a general decrease in OC concentration downstream for the size fractions < 0·25 mm. Bed sediments were dominated by size fractions coarser than 0·5 mm (80 per cent of the total distribution) and only about 2 per cent in the fractions less than 0·13 mm. Combining information on OC concentration per size fraction and the mass contribution of each fraction to the whole sample, it was observed that fractions coarser than 0·5 mm had eight to 12 times the storage of OC per kilogram of bed sediments than the fractions finer than 0·13 mm. Weighted OC values for Manoa Stream were on average 6·7 g-OC kg−1, and these were similar to those reported in the literature for a variety of sediments in aquatic environments, both freshwater and marine. These data provide important information on the relative mass storage of OC in bed sediments and their longitudinal patterns in a tropical fluvial environment. Copyright © 1999 John Wiley & Sons, Ltd.     

A portable experimental hillslope for frozen ground studies

Frozen ground hydrological effects on runoff, storage, and release have been observed in the field and tested in numerical models, but few physical models of frozen slopes (at scales from 1 to 15 m) exist partly because the design of such an experiment requires new engineering design for realistic whole‐slope freezing and physical model innovation. Here, we present a new freezable tilting hillslope physical model for hydrological system testing under a variety of climate conditions with the ability to perform multiple (up to 20 per year) freeze–thaw cycles. The 4 × 2 m hillslope is mobile and tiltable on the basis of a modified tri‐axle 4.88‐m (16′) dump trailer to facilitate testing multiple configurations. The system includes controllable boundary conditions on all surfaces; examples of side and baseflow boundary conditions include permeable membranes, impermeable barriers, semipermeable configurations, and constant head conditions. To simulate cold regions and to freeze the hillslope in a realistic and controlled manner, insulation and a removable freezer system are incorporated onto the top boundary of the hillslope. The freezing system is designed to expedite the freezing process by the addition of a 10,130‐KJ (9,600‐BTU) refrigeration coil to the top‐centre of the insulated ceiling. Centre placement provides radial freezing of the hillslope in a top‐down fashion, similar to what natural systems encounter in the environment. The perimeter walls are insulated with 100 mm of spray foam insulation, whereas the base of the hillslope is not insulated to simulate natural heat fluxes beneath the frozen layer of soil. Our preliminary testing shows that covers can be frozen down to ?10 °C in approximately 7 days, with subsequent thaw on a similar time frame.     

Modelling stream flow for use in ecological studies in a large,arid zone river,central Australia

Australian arid zone ephemeral rivers are typically unregulated and maintain a high level of biodiversity and ecological health. Understanding the ecosystem functions of these rivers requires an understanding of their hydrology. These rivers are typified by highly variable hydrological regimes and a paucity, often a complete absence, of hydrological data to describe these flow regimes. A daily time‐step, grid‐based, conceptual rainfall–runoff model was developed for the previously uninstrumented Neales River in the arid zone of northern South Australia. Hourly, logged stage data provided a record of stream‐flow events in the river system. In conjunction with opportunistic gaugings of stream‐flow events, these data were used in the calibration of the model. The poorly constrained spatial variability of rainfall distribution and catchment characteristics (e.g. storage depths) limited the accuracy of the model in replicating the absolute magnitudes and volumes of stream‐flow events. In particular, small but ecologically important flow events were poorly modelled. Model performance was improved by the application of catchment‐wide processes replicating quick runoff from high intensity rainfall and improving the area inundated versus discharge relationship in the channel sections of the model. Representing areas of high and low soil moisture storage depths in the hillslope areas of the catchment also improved the model performance. The need for some explicit representation of the spatial variability of catchment characteristics (e.g. channel/floodplain, low storage hillslope and high storage hillslope) to effectively model the range of stream‐flow events makes the development of relatively complex rainfall–runoff models necessary for multisite ecological studies in large, ungauged arid zone catchments. Grid‐based conceptual models provide a good balance between providing the capacity to easily define land types with differing rainfall–runoff responses, flexibility in defining data output points and a parsimonious water‐balance–routing model. Copyright © 2004 John Wiley & Sons, Ltd.     

Interactions and connectivity between runoff generation processes of different spatial scales

Monitoring 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. ( 2012 ) 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.     

Landscape position,surface hydraulic gradients and erosion processes

Different hydraulic gradients, especially due to seepage or drainage, at different locations on a hillslope profile may have a profound effect on the dominant erosion processes. A laboratory study was designed to simulate hillslope processes and quantify effects of surface hydraulic gradients on erosion for a Glynwood clay loam soil (fine, illitic, mesic Aquic Hapludalf). A 5 m long, 1·2 m wide soil pan was used at 5 and 10 per cent slopes with an external watering tube to vary the soil bed's hydrological conditions. Different combinations of slope steepness with seepage or drainage gradients were used to simulate the hydrologic conditions on a 5 m segment of a hillslope profile. Runoff samples were taken during rainfall-only and rainfall with added inflow. Results showed that, under drainage conditions, interrill processes dominated and rilling was limited. The surface contained scattered crescent-shaped pits after the run. Under seepage conditions, rilling processes dominated and the inflow introduced at the top of the soil pan further accelerated the headward erosion of the rills. Erosion rates increased by as much as 60 times under seepage conditions representative of the lower backslope when compared to drainage conditions that generally occur at the upper backslope. This indicated that rills and gullies on backslopes and footslopes may be catalysed or enhanced by seepage conditions rather than form from flow hydraulic shear stress alone. An understanding of spatial and temporal changes that affect both hillslope hydrology and erosional processes is needed to develop accurate process-based erosion prediction models. This knowledge may lead to different management practices on landscape positions where seepage occurs. © 1998 John Wiley & Sons, Ltd.     

Subsurface flow responses of a small forested catchment in the Ouachita mountains

Despite considerable research performed on forested catchments in the Ouachita Mountains of Oklahoma and Arkansas, little information on hydrological processes in operation is available. Based on catchment physical characteristics, subsurface flow was thought to be an important hydrological process in the region. Therefore, this study was undertaken to determine the occurrence, rates, timing and volumes of subsurface flow, and to estimate the importance of subsurface flow as a streamflow generating process. Subsurface flow was collected from three hillslope sites on a 7.7 ha forested catchment. Hillslope sites drained through natural seepage faces located near stream channels. Subsurface flow was collected from three depths at each hillslope site, below the litter layer, below the a horizon, and within the B horizon (Bt21). Subsurface flow occurred and was measured during 11 of 31 rainfall events. Subsurface flow responded rapidly to the initiation of and to changes in intensity of rainfall at all depths. the rapid response was indicative of flow through soil macropores. B horizon subsurface flow commenced within 10 to 180 min of the initiation of rainfall. Multiple linear regression showed that the volume of subsurface flow generated during a given storm was directly related to rainfall depth and a 7-day antecedent precipitation index used to represent antecedent water content. About 67 per cent of the total subsurface flow collected during the study was produced in one large storm under wet antecedent conditions. the storm was equal to the 2-year, 24-hour storm for the region. Measured subsurface flow volumes were extended to the watershed scale to provide estimates of catchment-wide contributions to streamflow. It was estimated that subsurface flow contributed from 1 to 48 per cent of total quickflow measured at the catchment outlet. Based on the timing of subsurface flow, it was estimated that subsurface flow May, contribute up to 70 per cent of quickflow before and soon after peak flow.