Summer flows in experimental catchments with different forest covers, Chile

Runoff and peak flows in four experimental catchments with different land uses are analyzed for summer periods. The catchments have a rainy temperate climate with annual precipitations between 2000 and 2500 mm, 70% of which is concentrated in the winter period between May and August. The final harvest of the forest plantation in one of these catchments generated increases in summer runoff. Also, differences between the maximum instantaneous discharge and the flow at the beginning of the storm then almost duplicated those registered in rainfall events of similar magnitude when the catchment was fully forested. Runoff analysis in this catchment is difficult because the two post-harvesting summer periods are much wetter than the two pre-harvesting ones but a double mass analysis shows the effect of harvesting clearly. In a paired catchment study, low cover in one of the two neighbour catchments explains higher direct runoff and base flows although lower maximum instantaneous specific discharge occurred in the less vegetated but larger catchment. Low vegetation cover explains increases in summer flows, although the size, topography, rainfall conditions, road density, extent of affected area and runoff generation processes play an important role in the hydrological effects of different land uses.     

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.     

A physically based distributed subsurface–surface flow dynamics model for forested mountainous catchments

This study was designed to develop a physically based hydrological model to describe the hydrological processes within forested mountainous river basins. The model describes the relationships between hydrological fluxes and catchment characteristics that are influenced by topography and land cover. Hydrological processes representative of temperate basins in steep terrain that are incorporated in the model include intercepted rainfall, evaporation, transpiration, infiltration into macropores, partitioning between preferential flow and soil matrix flow, percolation, capillary rise, surface flow (saturation‐excess and return flow), subsurface flow (preferential subsurface flow and baseflow) and spatial water‐table dynamics. The soil–vegetation–atmosphere transfer scheme used was the single‐layer Penman–Monteith model, although a two‐layer model was also provided. The catchment characteristics include topography (elevation, topographic indices), slope and contributing area, where a digital elevation model provided flow direction on the steepest gradient flow path. The hydrological fluxes and catchment characteristics are modelled based on the variable source‐area concept, which defines the dynamics of the watershed response. Flow generated on land for each sub‐basin is routed to the river channel by a kinematic wave model. In the river channel, the combined flows from sub‐basins are routed by the Muskingum–Cunge model to the river outlet; these comprise inputs to the river downstream. The model was applied to the Hikimi river basin in Japan. Spatial decadal values of the normalized difference vegetation index and leaf area index were used for the yearly simulations. Results were satisfactory, as indicated by model efficiency criteria, and analysis showed that the rainfall input is not representative of the orographic lifting induced rainfall in the mountainous Hikimi river basin. Also, a simple representation of the effects of preferential flow within the soil matrix flow has a slight significance for soil moisture status, but is insignificant for river flow estimations. Copyright © 2005 John Wiley & Sons, Ltd.     

Physically‐based modelling of double‐peak discharge responses at Slapton Wood catchment

Heavy winter rainfall produces double‐peak hydrographs at the Slapton Wood catchment, Devon, UK. The first peak is saturation‐excess overland flow in the hillslope hollows and the second (i.e. the delayed peak) is subsurface stormflow. The physically‐based spatially‐distributed model SHETRAN is used to try to improve the understanding of the processes that cause the double peaks. A three‐stage (multi‐scale) approach to calibration is used: (1) water balance validation for vertical one‐dimensional flow at arable, grassland and woodland plots; (2) two‐dimensional flow for cross‐sections cutting across the stream valley; and (3) three‐dimensional flow in the full catchment. The main data are for rainfall, stream discharge, evaporation, soil water potential and phreatic surface level. At each scale there was successful comparison with measured responses, using as far as possible parameter values from measurements. There was some calibration but all calibrated values at one scale were used at a larger scale. A large proportion of the subsurface runoff enters the stream from three dry valleys (hillslope hollows), and previous studies have suggested convergence of the water in the three large hollows as being the major mechanism for the production of the delayed peaks. The SHETRAN modelling suggests that the hillslopes that drain directly into the stream are also involved in producing the delayed discharges. The model shows how in the summer most of the catchment is hydraulically disconnected from the stream. In the autumn the catchment eventually ‘wets up’ and shallow subsurface flows are produced, with water deflected laterally along the soil‐bedrock interface producing the delayed peak in the stream hydrograph. Copyright © 2007 John Wiley & Sons, Ltd.     

Decoupling the effects of deforestation and climate variability in the Tapajós river basin in the Brazilian Amazon

Tropical river basins are experiencing major hydrological alterations as a result of climate variability and deforestation. These drivers of flow changes are often difficult to isolate in large basins based on either observations or experiments; however, combining these methods with numerical models can help identify the contribution of climate and deforestation to hydrological alterations. This paper presents a study carried out in the Tapaj?s River (Brazil), a 477,000 km² basin in South‐eastern Amazonia, in which we analysed the role of annual land cover change on daily river flows. Analysis of observed spatial and temporal trends in rainfall, forest cover, and river flow metrics for 1976 to 2008 indicates a significant shortening of the wet season and reduction in river flows through most of the basin despite no significant trend in annual precipitation. Coincident with seasonal trends over the past 4 decades, over 35% of the original forest (140,000 out of 400,000 km²) was cleared. In order to determine the effects of land clearing and rainfall variability to trends in river flows, we conducted hindcast simulations with ED2 + R, a terrestrial biosphere model incorporating fine scale ecosystem heterogeneity arising from annual land‐use change and linked to a flow routing scheme. The simulations indicated basin‐wide increases in dry season flows caused by land cover transitions beginning in the early 1990s when forest cover dropped to 80% of its original extent. Simulations of historical potential vegetation in the absence of land cover transitions indicate that reduction in rainfall during the dry season (mean of ?9 mm per month) would have had an opposite and larger magnitude effect than deforestation (maximum of +4 mm/month), leading to the overall net negative trend in river flows. In light of the expected increase in future climate variability and water infrastructure development in the Amazon and other tropical basins, this study presents an approach for analysing how multiple drivers of change are altering regional hydrology and water resources management.     

The influence of preferential flow on hillslope hydrology in a semi‐arid watershed (in the Spanish Dehesas)

Preferential flow is known to influence hillslope hydrology in many areas around the world. Most research on preferential flow has been performed in temperate regions. Preferential infiltration has also been found in semi‐arid regions, but its impact on the hydrology of these regions is poorly known. The aim of this study is to describe and quantify the influence of preferential flow on the hillslope hydrology from small scale (infiltration) to large scale (subsurface stormflow) in a semi‐arid Dehesa landscape. Precipitation, soil moisture content, piezometric water level and discharge data were used to analyse the hydrological functioning of a catchment in Spain. Variability of soil moisture content during the transition from dry to wet season (September to November) within horizontal soil layers leads to the conclusion that there is preferential infiltration into the soils. When the rainfall intensity is high, a water level rapidly builds up in the piezometer pipes in the area, sometimes even reaching soil surface. This water level also drops back to bedrock within a few hours (under dry catchment conditions) to days (under wet catchment conditions). As the soil matrix is not necessarily wet while this water layer is built up, it is thought to be a transient water table in large connected pores which drain partly to the matrix, partly fill up bedrock irregularities and partly drain through subsurface flow to the channels. When the soil matrix becomes wetter the loss of water from macropores to the matrix and bedrock decreases and subsurface stormflow increases. It may be concluded that the hillslope hydrological system consists of a fine matrix domain and a macropore domain, which have their own flow characteristics but which also interact, depending on the soil matrix and macropore moisture contents. The macropore flow can result in subsurface flow, ranging from 13% contribution to total discharge for a large event of high intensity rainfall or high discharge to 80% of total discharge for a small event with low intensity rainfall or low discharge. During large events the fraction of subsurface stormflow in the discharge is suppressed by the large amount of surface runoff. Copyright © 2008 John Wiley & Sons, Ltd.     

Isotope hydrology of a tropical coffee agroforestry watershed: Seasonal and event‐based analyses

Stable isotope variations are extremely useful for flow partitioning within the hydrologic cycle but remain poorly understood throughout the tropics, particularly in watersheds with rapidly infiltrating soils, such as Andisols in Central America. This study examines the fluctuations of stable isotope ratios (δ¹⁸O and δ²H) in the hydrologic components of a tropical coffee agroforestry watershed (~1 km²) with Andisol soils in Costa Rica. Samples were collected in precipitation, groundwater, springs, and stream water over 2 years. The local meteoric water line for the study site was δ²H = 8.5 δ¹⁸O + 18.02 (r² = 0.97, n = 198). The isotope ratios in precipitation exhibited an enriched trend during the dry season and a notable depletion at the beginning of the wet season. The δ¹⁸O compositions in groundwater (average = ?6.4‰, σ = 0.7) and stream water (average = ?6.7‰, σ = 0.6) were relatively stable over time, and both components exhibited more enriched values in 2013, which was the drier year. No strong correlation was observed between the isotope ratios and the precipitation amount at the event or daily time‐step, but a correlation was observed on a monthly scale. Stream water and base flow hydrograph separations based on isotope end‐member estimations showed that pre‐event water originating from base flow was prevalent. However, isotope data indicate that event water originating from springs appears to have been the primary driver of initial rises in stream flow and peak flows. These results indicate that isotope sampling improves the understanding of water balance components, even in a tropical humid location, where significant variations in rainfall challenge current modelling efforts. Further research using fine‐scale hydrometric and isotopic data would enhance understanding the processes driving spring flow generation in watersheds.     

Using hydrogeochemical signatures of stream water to assess pathways for rainfall events: towards a predictive model

Thanks to its simple division into agricultural and forestry land use, the Corbeira catchment (Galicia, Spain) is used as a case study to build a predictive model using hydrogeochemical signatures. Stream data acquired under recessional flow conditions over a one year period were obtained from a sampling station near the downstream end of the catchment, and using principal component analysis, it is shown that some of the analytical parameters are covariant, and some are negatively correlated. These findings support inferences about the pathways of rainfall in the catchment. Specific signatures may be associated with the dominant hydrological source, either surface runoff or subsurface waters: additionally, the dominant land use in that part of the catchment, where the flow originated, can also be predicted. The dominant runoff shows a strong covariance between suspended solids (SS) and particulate phosphorus (PP), with a clear negative correlation with pH. Dissolved organic carbon (DOC) data are associated with this covariant set when these compounds are available in the soils in question. Dissolved phosphorus, total organic nitrogen and dissolved nitrates are also associated with the same covariant set when the runoff flows through areas of extensive agricultural use. The SS ? PP covariance is less significant at lower flows. Typical base flow regimes show a significant covariance between salinity and pH, with a marked negative correlation with SS ? PP set, confirming the dominance of subsurface waters in the baseflow, as expected. Seasonally divergent DOC ? SS behaviour proves to be a useful tracer for rainfall regimes. The DOC trend shows a sinusoidal annual variation in amplitude, determined by the rainfall regime. As a result, flow from the catchment is dominated by surface water whenever there is synchronicity between the peaks of DOC and SS. Copyright © 2013 John Wiley & Sons, Ltd.     

Base flow separation for soil erosion simulation in a granitic forested headwater catchment using a process-based model,GeoWEPP

Distributed erosion models, which simulate the physical processes of water flow and soil erosion, are effective for predicting soil erosion in forested catchments. Although subsurface flow through multiple pathways is dominant for runoff generation in forested headwater catchments, the process-based erosion model, Geo-spatial interface for Water Erosion Prediction Project(Geo WEPP), does not have an adequate subsurface component for the simulation of hillslope water flow. In the current study, t...     

Runoff generation in zero-order precambrian shield catchments: The stormflow response of a heterogeneous landscape

Hydrometric and isotopic (oxygen-18) observations were used to delineate the runoff processes operating in several headwater catchments on the Precambrian Shield of Canada. The catchments comprise patches of conifer forest situated on thin soils among areas of lichen-covered granitic bedrock. Horton overland flow occurred from the lichen-bedrock areas in all precipitation events that exceeded 4–6 mm. Runoff from the forest stands occurred mainly as subsurface stormflow, but in some instances saturation overland flow was observed. The occurrence of saturation overland flow was controlled by the topography of the bedrock beneath the forest soils. The area contributing runoff and the pathway by which water was conveyed to the catchment outflow switched from the open lichen-bedrock areas producing overland flow on the rising limb of the storm hydrograph to the forest stands contributing subsurface stormflow on the recession limb of the hydrograph. The areal extent and position of the landscape units in the basin were important to the rate and magnitude of stormflow production. Runoff was generated from the catchments only during and immediately after snowmelt and/or rainfall events. The catchments were dry and/or frozen for about 70% of the year.     

Review of field observations of runoff generation in the tropics

Over the last thirty years, French hydrologists have observed flow characteristics at the outlets of small tropical watersheds, characteristics which vary between arid zones, humid zones and equatorial rainforest zones. The main physical phenomena related to the various types of runoff generation have been observed and analysed. The generation of runoff and subsurface runoff following tropical rainstorms depends on soil surface formations, texture and structure as well as vegetation cover, geological substrata, weathered material and climate.

Analysis is also made of experimental rainfall-runoff relationships, the shape of flood hydrographs and the conditions necessary for the occurrence of steady base flow. Scale effects between experimental plot and watershed have been studied to identify any effect this may have upon runoff generation and flow. Tables detail watershed characteristics and relationships between these and rainfall and runoff.     

Storm runoff processes and subcatchment characteristics in a new zealand hill country catchment

The quickflow responses of six subcatchment areas in a small hill country catchment in the Craigieburn Range, South Island, New Zealand, were compared for a range of storm sizes, rainfall intensities and antecedent wetness conditions. Topography and soil characteristics suggested that all subcatchments would receive subsurface stormflow input, but that some would receive larger saturation overland flow inputs than others. Quickflow yields and response ratios were positively correlated with storm size and antecedent wetness conditions in the subcatchment most suited to producing saturation overland flow. In subcatchments more likely to be dominated by subsurface flow, quickflow yields and response ratios were positively correlated with storm size, but were either not correlated, or negatively correlated, with antecedent wetness. Quickflow responses were either not significantly or negatively correlated with rainfall intensity variables. Quickflow from the subcatchment most suited to produce saturation overland flow providing an increasing proportion of total catchment quickflow in larger storms and as antecedent conditions became wetter. Subcatchment responses varied greatly in space and time and there was less pattern to the variation than had been expected. Where topographic and pedologic conditions permit substantial responses to storm rainfall by both saturation overland flow and subsurface stormflow, simple topographic and soil indicators may not be useful guides to the relative importance of runoff mechanisms, or to the identification of runoff-source areas.     

Field investigation and modelling of coupled stream discharge and shallow water‐table dynamics in a small Mediterranean catchment (Sardinia)

In the semi‐arid Mediterranean environment, the rainfall–runoff relationships are complex because of the markedly irregular patterns in rainfall, the seasonal mismatch between evaporation and rainfall, and the spatial heterogeneity in landscape properties. Watersheds often display considerable non‐linear threshold behavior, which still make runoff generation an open research question. Our objectives in this context were: to identify the primary processes of runoff generation in a small natural catchment; to test whether a physically based model, which takes into consideration only the primary processes, is able to predict spatially distributed water‐table and stream discharge dynamics; and to use the hydrological model to increase our understanding of runoff generation mechanisms. The observed seasonal dynamics of soil moisture, water‐table depth, and stream discharge indicated that Hortonian overland‐flow was negligible and the main mechanism of runoff generation was saturated subsurface‐flow. This gives rise to base‐flow, controls the formation of the saturated areas, and contributes to storm‐flow together with saturation overland‐flow. The distributed model, with a 1D scheme for the kinematic surface‐flow, a 2D sub‐horizontal scheme for the saturated subsurface‐flow, and ignoring the unsaturated flow, performed efficiently in years when runoff volume was high and medium, although there was a smoothing effect on the observed water‐table. In dry years, small errors greatly reduced the efficiency of the model. The hydrological model has allowed to relate the runoff generation mechanisms with the land‐use. The forested hillslopes, where the calibrated soil conductivity was high, were never saturated, except at the foot of the slopes, where exfiltration of saturated subsurface‐flow contributed to storm‐flow. Saturation overland‐flow was only found near the streams, except when there were storm‐flow peaks, when it also occurred on hillslopes used for pasture, where soil conductivity was low. The bedrock–soil percolation, simulated by a threshold mechanism, further increased the non‐linearity of the rainfall–runoff processes. Copyright © 2013 John Wiley & Sons, Ltd.     

Modelling impacts of agricultural practice on flood peaks in upland catchments: An application of the distributed TOPMODEL

Upland agricultural land management activities such as grazing, vegetation burning, and bare ground restoration impact hydrological elements of headwater catchments, many of which may be important for downstream flood peaks (e.g., overland flow and soil water storage). However, there is poor understanding of how these management practices affect river flow peaks during high magnitude rainfall events. Using the distributed TOPMODEL, spatial configurations of land management were modelled to predict flood response in an upland catchment, which contains different regions operating subsidized agricultural stewardship schemes. Heavy grazing leading to soil compaction and loss of vegetation cover in stewardship regions covering 79.8% of the catchment gave a 42‐min earlier flow peak, which was 82.2% higher (under a 1‐hr 15‐mm storm) than the current simulated hydrograph. Light grazing over the same regions of the catchment had much less influence on river flow peaks (18 min earlier and 32.9% increase). Rotational burning (covering 8.8% of the catchment), most of which is located in the headwater areas, increased the peak by 3.2% in the same rainfall event. Vegetation restoration with either Eriophorum or Sphagnum (higher density) in bare areas (5.8%) of the catchment provided a reduction of flood peak (3.9% and 5.2% in the 15‐mm storm event), whereas the same total area revegetated with Sphagnum in riparian regions delivered a much larger decrease (15.0%) in river flow peaks. We show that changes of vegetation cover in highly sensitive areas (e.g., near‐stream zones) generate large impacts on flood peaks. Thus, it is possible to design spatially distributed management systems for upland catchments, which reduce flood peaks while at the same time ensuring economic viability for upland farmers.     

Modelling the impacts of land-cover change on streamflow dynamics of a tropical rainforest headwater catchment

Abstract

A modelling experiment is used to examine different land-use scenarios ranging from extreme deforestation (31% forest cover) to pristine (95% forest cover) conditions and related Payment for Ecosystem Services (PES) schemes to assess whether a change in streamflow dynamics, discharge extremes and mean annual water balance of a 73.4-km² tropical headwater catchment in Costa Rica could be detected. A semi-distributed, conceptual rainfall–runoff model was adapted to conceptualize the empirically-based, dominant hydrological processes of the study area and was multi-criteria calibrated using different objective functions and empirical constraints on model simulations in a Monte Carlo framework to account for parameter uncertainty. The results suggest that land-use change had relatively little effect on the overall mean annual water yield (<3%). However, streamflow dynamics proved to be sensitive in terms of frequency, timing and magnitude of discharge extremes. For low flows and peak discharges of return periods greater than one year, land use had a minor influence on the runoff response. Below these thresholds (<1-year return period), forest cover potentially decreased runoff peaks and low flows by as much as 10%, and non-forest cover increased runoff peaks and low flows by up to 15%. The study demonstrated the potential for using hydrological modelling to help identify the impact of protection and reforestation efforts on ecosystem services.

Editor Z.W. Kundzewicz

Citation Birkel, C., Soulsby, C., and Tetzlaff, D., 2012. Modelling the impacts of land-cover change on streamflow dynamics of a tropical rainforest headwater catchment. Hydrological Sciences Journal, 57 (8), 1543–1561.     

Reduced channel morphological response to urbanization in a flood‐dominated humid tropical environment

Urbanization through the addition of impervious cover can alter catchment hydrology, often resulting in increased peak flows during floods. This phenomenon and the resulting impact on stream channel morphology is well documented in temperate climatic regions, but not well documented in the humid tropics where urbanization is rapidly occurring. This study investigates the long‐term effects of urbanization on channel morphology in the humid sub‐tropical region of Puerto Rico, an area characterized by frequent high‐magnitude flows, and steep coarse‐grained rivers. Grain size, low‐flow channel roughness, and the hydraulic geometry of streams across a land‐use gradient that ranges from pristine forest to high density urbanized catchments are compared. In areas that have been urbanized for several decades changes in channel features were measurable, but were smaller than those reported for comparable temperate streams. Decades of development has resulted in increased fine sediment and anthropogenic debris in urbanized catchments. Materials of anthropogenic origin comprise an average of 6% of the bed material in streams with catchments with 15% or greater impervious cover. At‐a‐station hydraulic geometry shows that velocity makes up a larger component of discharge for rural channels, while depth contributes a larger component of discharge in urban catchments. The average bank‐full cross‐sectional area of urbanized reaches was 1.5 times larger than comparable forested reaches, and less than the world average increase of 2.5. On average, stream width at bank‐full height did not change with urbanization while the world average increase is 1.5 times. Overall, this study indicates that the morphologic changes that occur in response to urban runoff are less in channels that are already subject to frequent large magnitude storms. Furthermore, this study suggests that developing regions in the humid tropics shouldn't rely on temperate analogues to determine the magnitude of impact of urbanization on stream morphology. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling the influence of hydrogeological parameters on low flow in hilly and mountainous regions of Serbia / Modélisation de l'influence de paramètres hydrogéologiques sur les basses eaux dans des régions de collines et de montagnes de Serbie

Abstract

During dry weather periods of the year with long rainless intervals, streams slacken to what is generally termed “low flow”. This work presents an analysis of the influence of hydrogeology on low flows, using multiple linear regression, in natural medium and small streams in hilly and mountainous regions of Serbia. The study cases encompass 61 gauged catchments south of the rivers Sava and Danube. Characteristic relevant minimum mean 30-day flows of 80- or 95-percentile exceedence (Q _80%, Q _95%) are taken as dependent variables. Independent variables are the observable hydrogeological quantities: catchment area upstream of a gauging station; surface area of a hydrogeological soil category in a catchment; number of perennial springs of minimum flow higher than or equal to 1 L/s in a gauged catchment; number of perennial springs, each of minimum flow higher than or equal to 1 L/s, in a given hydrogeological soil category of the catchment; cumulative perennial spring flow of minimum single flow higher than or equal to 1 L/s in a catchment; and cumulative perennial spring flow of minimum single flow higher than or equal to 1 L/s in a hydrogeological soil category of a catchment. Through multiple linear regression between the characteristic relevant low flow and the hydrogeological elements, 16 models are developed and analysed, each based on a different combination of hydrogeological elements and characteristic low flow. The regional relationships developed for the minimum mean 30-day flows of 80- and 95-percentile exceedences are evaluated. The statistical tests of the representation quality of each multiple regression relationship show that the models justify the use of hydrogeological elements.