Transport of the environmental radionuclides in an alpine watershed

The radionuclides¹³⁷Cs,²¹⁰Pb and⁷Be have been examined in the alpine Rhoˆne watershed (Switzerland) during a period of two years in order to evaluate their usefulness as tracers of the removal and transport rate of top-soil particles and particle-reactive contaminants of atmospheric origin. The specific activities of the radionuclides in fluvial suspension show a distinct seasonal pattern which depends on the hydrologic regime of the stream and the sources of the suspended matter. Input-output budgets based on the atmospheric deposition and fluvial removal of three radionuclides in the alpine Rhoˆne watershed are used to estimate their erosional residence times. The simplest one box model yields mean residence times of about 800 and 1400 years for¹³⁷Cs and²¹⁰Pb, respectively. The removal rate of short-lived⁷Be suggests that a part of the watershed (0.6–2.3% of the total surface) is exposed to a rapid erosion, in which the mean residence time of the radionuclides is in the range of 1–220 days. This has little influence on the calculated residence time of¹³⁷Cs but increases the estimated residence time of²¹⁰Pb in soil to over 1800 years. The use of²¹⁰Pb-⁷Be pair in fluvial output may be very helpful in the assessment of the impact of atmospheric pollutants on the water quality in rivers and lakes.     

Effects of climatic variability on the hydrologic response of a freshwater watershed

A generalized watershed model was used to evaluate the effects of global climate changes on the hydrologic responses of freshwater ecosystems. The Enhanced Trickle Down (ETD) model was applied to W-3 watershed located near Danville, Vermont. Eight years of field data was used to perform model calibration and verification and the results were presented in Nikolaidis et al., (1993). Results from the Goddard Institute for Space Studies (GISS) and the Geophysical Fluid Dynamics Laboratory (GFDL) general circulation models which simulated the doubling of present day atmospheric CO₂ scenarios were used to perform the hydrologic simulations for the W-3 watershed. The results indicate that the W-3 watershed will experience increases in annual evapotranspiration and decreases in annual outflow and soil moisture. Stochastic models that simulate collective statistical properties of meteorological time series were developed to generate data to drive the ETD model in a Monte-Carlo fashion for quantification of the uncertainty in the model predictions due to input time series. This coupled deterministic and stochastic model was used to generate probable scenarios of future hydrology of the W-3 watershed. The predicted evapotranspiration and soil moisture under doubling present day atmospheric CO₂ scenarios exceed the present day uncertainty due to input time series by a factor greater than 2. The results indicate that the hydrologic response of the W-3 watershed will be significantly different than its present day response. The Enhanced Trickle Down model can be used to evaluate land surface feedbacks and assessing water quantity management in the event of climate change.     

Estimation of baseflow residence times in watersheds from the runoff hydrograph recession: method and application in the Neversink watershed,Catskill Mountains,New York

A method for estimation of mean baseflow residence time in watersheds from hydrograph runoff recession characteristics was developed. Runoff recession characteristics were computed for the period 1993–96 in the 2 km² Winnisook watershed, Catskill Mountains, southeastern New York, and were used to derive mean values of subsurface hydraulic conductivity and the storage coefficient. These values were then used to estimate the mean baseflow residence time from an expression of the soil contact time, based on watershed soil and topographic characteristics. For comparison, mean baseflow residence times were calculated for the same period of time through the traditional convolution integral approach, which relates rainfall δ¹⁸O to δ¹⁸O values in streamflow. Our computed mean baseflow residence time was 9 months by both methods. These results indicate that baseflow residence time can be calculated accurately using recession analysis, and the method is less expensive than using environmental and/or artificial tracers. Published in 2002 by John Wiley & Sons, Ltd.     

Macroscopic representation of structural geometry for simulating water and solute movement in dual-porosity media

The structure of macroporous or aggregated soils and fractured rocks is generally so complex that it is impractical to measure the geometry at the microscale (i.e., the size and the shape of soil aggregates or rock matrix blocks, and the myriad of fissures or fractures), and use such data in geometry-dependent macroscale flow and transport models. This paper analyzes a first-order type dual-porosity model which contains a geometry-dependent coefficient, β, in the mass transfer term to macroscopically represent the size and shape of soil or rock matrix blocks. As a reference, one- and two-dimensional geometry-based diffusion models were used to simulate mass transport into and out of porous blocks of defined shapes. Estimates for β were obtained analytically for four different matrix block geometries. Values for β were also calculated by directly matching analytical solutions of the diffusion models for a number of selected matrix block geometries to results obtained with the first-order model assuming standard boundary conditions. Direct matching improved previous results for cylindrical macropore geometries, especially when relatively small ratios between the outer soil mantle and the radius of the inner cylinder were used. Results of our analysis show that β is closely related to the ratio of the effective surface area available for mass transfer, and the soil matrix volume normalized by the effective characteristic length of the matrix system. Using values of β obtained by direct matching, an empirical function is derived to estimate macroscopic geometry coefficients from medium properties which in principle are measurable. The method permits independent estimates of β, thus allowing the dual-porosity approach eventually to be applied to media with complex and mixed types of structural geometry.     

Comparative streamflow characteristics in urbanizing basins in the Portland Metropolitan Area,Oregon, USA

This study investigates changes in streamflow characteristics for urbanizing watersheds in the Portland Metropolitan Area of Oregon for the period from 1951 to 2000. The objective of this study was to assess how mean annual runoff ratio, mean seasonal runoff ratio, annual peak runoff ratio, changes in streamflow in response to storm amount, the fraction of time that the daily mean flow exceeds the annual mean flow, 3‐day recession constants, and dry/wet flow ratio vary among watersheds with different degrees of urban development. There were no statistically significant changes in annual runoff ratio and annual peak runoff ratio for the mixed land‐use watershed (Tualatin River watershed) and the urban watershed (Johnson Creek watershed) during the entire study period. The Tualatin River watershed, where most of the urban development occurred in a lower part of the watershed, showed a statistically significant increase in annual peak runoff ratio during the 1976 and 2000 period. The Upper Tualatin River watershed illustrated a significant decrease in annual peak runoff ratio for the entire study period. With significant differences in seasonal runoff ratio, only Johnson Creek exhibited a significant increase in both wet and dry season runoff ratios. Streamflow during storm events declined rapidly in the urban watershed, with a high 3‐day recession constant. At an event storm scale, streamflow in Fanno Creek, which is the most urbanized watershed, responded quickly to precipitation input. The fraction of time that the daily mean flow exceeded the annual mean flow and dry/wet flow ratio are all lower in Johnson Creek. This suggests a shorter duration of storm runoff and lower baseflow in the urbanized watershed when compared to the mixed land use watershed. The findings of this study demonstrate the importance of spatial and temporal scale, climate variability, and basin physiographic characteristics in detecting the hydrologic effects of urbanization in the Pacific Northwest of the USA. Copyright © 2006 John Wiley & Sons, Ltd.     

Controls on event runoff coefficients in the eastern Italian Alps

Analyses of event runoff coefficients provide essential insight on catchment response, particularly if a range of catchments and a range of events are compared by a single indicator. In this study we examine the effect of climate, geology, land use, flood types and initial soil moisture conditions on the distribution functions of the event runoff coefficients for a set of 14 mountainous catchments located in the eastern Italian Alps, ranging in size from 7.3 to 608.4 km². Runoff coefficients were computed from hourly precipitation, runoff data and estimates of snowmelt. A total of 535 events were analysed over the period 1989–2004. We classified each basin using a “permeability index” which was inferred from a geologic map and ranged from “low” to “high permeability”. A continuous soil moisture accounting model was applied to each catchment to classify ‘wet’ and ‘dry’ initial soil moisture conditions. The results indicate that the spatial distribution of runoff coefficients is highly correlated with mean annual precipitation, with the mean runoff coefficient increasing with mean annual precipitation. Geology, through the ‘permeability index’, is another important control on runoff coefficients for catchments with mean annual precipitation less than 1200 mm. Land use, as indexed by the SCS curve number, influences runoff coefficient distribution to a lesser degree. An analysis of the runoff coefficients by flood type indicates that runoff coefficients increase with event snowmelt. Results show that there exists an intermediate region of subsurface water storage capacity, as indexed by a flow–duration curve-based index, which maximises the impact of initial wetness conditions on the runoff coefficient. This means that the difference between runoff coefficients characterised by wet and dry initial conditions is negligible both for basins with very large storage capacity and for basins with small storage capacity. For basins with intermediate storage capacities, the impact of the initial wetness conditions may be relatively large.     

The influence of freeze–thaw cycles of active soil layer on surface runoff in a permafrost watershed

As a result of global warming, the discharges from rivers in permafrost regions have varied significantly. However, its mechanism remains unclear. One of possible factors is active soil freeze–thaw cycle, which may influence surface runoff in the variation of permafrost water cycle processes. In this study, a typical permafrost watershed in the Qinghai-Tibet plateau was selected, its hydrological processes were monitored from 2004 to 2007, and the effects of the freezing and thawing depth of the soil active layer on runoff processes were assessed. The runoff modulus, runoff coefficient, direct runoff ratio, recession gradient and their seasonal variations were estimated and analyzed. The active soil dynamics and water budget were analyzed to prove the features of the surface runoff and the influences of active soil freeze–thaw processes. The primary factors influencing surface runoff processes during different seasons were analyzed by Principal Component Analysis (PCA) and statistical regression methods. The results showed that the high runoff coefficient and low direct runoff ratio were the main characteristics during the spring flood period (May–June) and during the autumn recession period (September). The runoff modulus and its year-to-year variability were the greatest in the summer flood period. The direct runoff ratio decreased from 0.43 in May to 0.29 in September, with the exception of the highest ratio, which occurred during the summer recession period (July). The active soil thawing in the upper layer of depth of 60 cm had contributed to increase in discharge, but the increase in thawing depth deeper than 60 cm led to a decrease in surface runoff and slowness in the recession process. Precipitation played a small role in the spring flood runoff and the autumn runoff. The soil active layer freeze–thaw variation, which affected seasonal soil water dynamic and water budget and reformed seasonal runoff characteristics, along with vegetation cover changes, is considered the potential major factor in control of the hydrological processes in the permafrost region.     

Simulated water budget of a small forested watershed in the continental/maritime hydroclimatic region of the United States

Annual streamflows have decreased across mountain watersheds in the Pacific Northwest of the United States over the last ~70 years; however, in some watersheds, observed annual flows have increased. Physically based models are useful tools to reveal the combined effects of climate and vegetation on long‐term water balances by explicitly simulating the internal watershed hydrological fluxes that affect discharge. We used the physically based Simultaneous Heat and Water (SHAW) model to simulate the inter‐annual hydrological dynamics of a 4 km² watershed in northern Idaho. The model simulates seasonal and annual water balance components including evaporation, transpiration, storage changes, deep drainage, and trends in streamflow. Independent measurements were used to parameterize the model, including forest transpiration, stomatal feedback to vapour pressure, forest properties (height, leaf area index, and biomass), soil properties, soil moisture, snow depth, and snow water equivalent. No calibrations were applied to fit the simulated streamflow to observations. The model reasonably simulated the annual runoff variations during the evaluation period from water year 2004 to 2009, which verified the ability of SHAW to simulate the water budget in this small watershed. The simulations indicated that inter‐annual variations in streamflow were driven by variations in precipitation and soil water storage. One key parameterization issue was leaf area index, which strongly influenced interception across the catchment. This approach appears promising to help elucidate the mechanisms responsible for hydrological trends and variations resulting from climate and vegetation changes on small watersheds in the region. Copyright © 2015 John Wiley & Sons, Ltd.     

Water yield,annual peaks and exposure in mountainous terrain

In steep, mountainous terrain of the western United States, where summer rains are negligible and annual runoff is primarily dependent on melting snow, the runoff distribution is such that a high energy recipient (HE) watershed generally yields more water per unit area than a low energy recipient (LE) watershed except during the latter's flood peak period, at which time the ratio of yields is reversed. HE watersheds peak sooner than their corresponding LE watersheds, but also vary more in their peaking date. The ratio of maximum peak flows may be greater or less than 1.0 depending on environmental circumstances.     

A framework to estimate sediment loads using distributions with covariates: Beaurivage River watershed (Québec,Canada)

A framework to estimate sediment loads based on the statistical distribution of sediment concentrations and various functional forms relating distribution characteristics (e.g. mean and variance) to covariates is developed. The covariates are used as surrogates to represent the main processes involved in sediment generation and transport. Statistical models of increasing complexity are built and compared to assess their relative performance using available sediment concentration and covariate data. Application to the Beaurivage River watershed (Québec, Canada) is conducted using data for the 1989–2004 period. The covariates considered in this application are streamflow and calendar day. A comparison of different statistical models shows that, in this case, the log‐normal distribution with a mean value depending on streamflow (power law with an additive term) and calendar day (sinusoidal), a constant coefficient of variation for streamflow dependence and a constant standard deviation for calendar day dependence provide the best result. Model parameters are estimated using the maximum likelihood estimation technique. The selected model is then used to estimate the distribution of annual sediment loads for the Beaurivage River watershed for a selected period. A bootstrap parametric method is implemented to account for uncertainties in parameter values and to build the distributions of annual loads. Comparison of model results with estimates obtained using the empirical ratio estimator shows that the latter were rarely within the 0·1–0·9 quantile interval of the distributions obtained with the proposed approach. Copyright © 2008 John Wiley & Sons, Ltd.     

Spatial and temporal variability of 0‐ to 5‐m soil–water storage at the watershed scale

Dynamic relationships among rainfall patterns, soil water distribution, and plant growth are crucial for sustainable conservation of soil and water resources in water‐limited ecosystems. Spatial and temporal variation in deep soil water content at a watershed scale have not yet been characterized adequately due to the lack of deep soil water data. Deep soil–water storage (SWS) up to a depth of 5 m (n = 73) was measured at 19 sampling occasions at the LaoYeManQu watershed on the Chinese Loess Plateau (CLP). At a depth of 0–1.5 m, the annual mean SWS was highly correlated with rain intensity, and the correlation decreased with depth, but within the layers at 1.5–5.0 m, the changes in SWS indicated a lag between precipitation and the replenishment of soil water. Geostatistical parameters of SWS were also highly dependent on depth, and the mean SWS presented similar spatial structures in two adjacent layers. Temporal stability of SWS as indicated by mean relative difference, standard deviation of the relative difference (SDRD), and mean absolute bias error (MABE) was significantly weaker at the shallow than at deeper layers. Soil separates and organic carbon content controlled the spatial pattern of SWS at the watershed scale. One representative location (Site 57) was identified to estimate the mean SWS in the 1‐ to 5‐m layer of the watershed. Semivariograms of the SDRD and MABE were best fitted by an isotropic spherical model, and their spatial distributions were depth‐dependent. Both temporal stability and spatial variability of SWS increased over depth. This study is helpful for deep SWS estimation and sustainable management of soil and water on the CLP, and for other similar regions around the world.     

Assessment of transmission losses and groundwater recharge from runoff events in a wadi under shortage of data on lateral inflow,Negev, Israel

A hydrological–lithostratigraphical model was developed for assessment of transmission losses and groundwater recharge from runoff events in arid water courses where hydrological and meteorological records are incomplete. Water balance equations were established for reaches between hydrometric stations. Because rainfall and tributary flow data are scarce, lateral inflow, which is an essential component of the water balance equation, could not be estimated directly. The solution was obtained by developing a method which includes a hydrological–lithostratigraphical analogy. This is based on the following assumptions: (a) runoff resulting from a given rainfall event is related to the watershed surface lithology; (b) for a given event, the spatial distribution of runoff reflects the distribution of rainfall: and (c) transmission losses are uniquely related to the total inflow to the reach. The latter relationship, called the loss function, and the water balance equation comprise a model which simultaneously assesses lateral inflow and transmission losses for runoff events recorded at the terminal stations. The model was applied to three reaches of the arid Nahal Tsin in Israel. In this case study, the transmission losses were of the same order of magnitude as the flow at the major hydrometric stations. The losses were subdivided into channel moistening, which subsequently evaporates, and deep percolation, which recharges groundwater. For large runoff events, evaporation was substantially smaller than the losses. The mean annual recharge of groundwater from runoff events in the Tsin watershed was 4·1×10⁶ m³, while the mean annual flow volume at the major stations ranged from 0·6 to 1·5×10⁶ m³. Once in 100 years, the annual recharge may be seven times higher than the mean annual value, but the recharge during most years is very small. Copyright © 1999 John Wiley & Sons, Ltd.     

Simulation of spatially-distributed copper pollution in a large river basin using the ECOMAG-HM model

The semi-distributed physically-based model ECOMAG-HM was developed to simulate cycling of heavy metals in large river basins: on the surface, and in soil, groundwater and river water. The model was applied to study the spatial distribution and temporal dynamics of copper concentrations in watercourses of the Nizhnekamskoe Reservoir watershed in Russia. This watershed is characterized by high background concentrations of heavy metals due to wide occurrence of ore deposits and considerable concentrations of ore-parent elements in rocks. The model was found to adequately reproduce the spatial variation of the mean annual copper concentrations at different monitoring points of the river network. The mean annual specific copper washoff, with the surface and subsurface components of river runoff, and the total copper washoff from the watershed into the river network were calculated and mapped. The contributions of natural and anthropogenic factors to river water pollution by copper were evaluated.     

Effect of channel size on solute residence time distributions in rivers

The effect of channel size on residence time distributions (RTDs) of solute in rivers is investigated in this paper using tracer test data and the variable residence time (VART) model. Specifically, the investigation focuses on the influence of shear dispersion and hyporheic exchange on the shape of solute RTD, and how these two transport processes prevail in larger and smaller streams, respectively, leading to distinct tails of RTD. Simulation results show that (1) RTDs are dispersion-dependent and thereby channel-size (scale) dependent. RTDs increasing longitudinal dispersion coefficient. Small streams with negligible dispersion coefficient may display various types of RTD from upward curving patterns to a straight line (power-law distributions) and further to downward curving lognormal distributions when plotted in log–log coordinates. Moderate-sized rivers are transitional in terms of RTDs and commonly exhibit lognormal and power-law RTDs; (2) the incorporation of water and solute losses/gains in the VART model can improve simulation results and make parameter values more reasonable; (3) the ratio of time to peak concentration to the minimum mean residence time is equal to the recovery ratio of tracer. The relation provides a simple method for determining the minimum mean residence time; and (4) the VART model is able to reproduce various RTDs observed in rivers with 3–4 fitting parameters while no user-specified RTD functions are needed.     

Hydrologic functioning of the deep critical zone and contributions to streamflow in a high‐elevation catchment: Testing of multiple conceptual models

High‐elevation mountain catchments are often subject to large climatic and topographic gradients. Therefore, high‐density hydrogeochemical observations are needed to understand water sources to streamflow and the temporal and spatial behaviour of flow paths. These sources and flow paths vary seasonally, which dictates short‐term storage and the flux of water in the critical zone (CZ) and affect long‐term CZ evolution. This study utilizes multiyear observations of chemical compositions and water residence times from the Santa Catalina Mountains Critical Zone Observatory, Tucson, Arizona to develop and evaluate competing conceptual models of seasonal streamflow generation. These models were tested using endmember mixing analysis, baseflow recession analysis, and tritium model “ages” of various catchment water sources. A conceptual model involving four endmembers (precipitation, soil water, shallow, and deep groundwater) provided the best match to observations. On average, precipitation contributes 39–69% (55 ± 16%), soil water contributes 25–56% (41 ± 16%), shallow groundwater contributes 1–5% (3 ± 2%), and deep groundwater contributes ~0–3% (1 ± 1%) towards annual streamflow. The mixing space comprised two principal planes formed by (a) precipitation‐soil water‐deep groundwater (dry and summer monsoon season samples) and (b) precipitation‐soil water‐shallow groundwater (winter season samples). Groundwater contribution was most important during the wet winter season. During periods of high dynamic groundwater storage and increased hydrologic connectivity (i.e., spring snowmelt), stream water was more geochemically heterogeneous, that is, geochemical heterogeneity of stream water is storage‐dependent. Endmember mixing analysis and ³H model age results indicate that only 1.4 ± 0.3% of the long‐term annual precipitation becomes deep CZ groundwater flux that influences long‐term deep CZ development through both intercatchment and intracatchment deep groundwater flows.     

Estimation of potential evapotranspiration in the mountainous Panama Canal watershed

Spatially distributed hydrometeorological and plant information within the mountainous tropical Panama Canal watershed is used to estimate parameters of the Penman–Monteith evapotranspiration formulation. Hydrometeorological data from a few surface climate stations located at low elevations in the watershed are complemented by (a) typical wet‐ and dry‐season fields of temperature, wind, water vapour and pressure produced by a mesoscale atmospheric model with a 3 × 3 km² spatial and hourly temporal resolution, and (b) leaf area index fields estimated over the watershed during a few years using satellite data with two different spatial and temporal resolutions. The mesoscale model estimates of spatially distributed surface hydrometeorological variables provide the basis for the extrapolation of the surface climate station data to produce input for the Penman–Monteith equation. The satellite information and existing digital spatial databases of land use and land cover form the basis for the estimation of Penman–Monteith spatially distributed parameter values. Spatially distributed 3 × 3 km² potential evapotranspiration estimates are obtained for the 3300 km² Panama Canal watershed. Estimates for Gatun Lake within the watershed are found to reproduce well the monthly and annual lake evaporation obtained from submerged pans. Sensitivity analysis results of potential evapotranspiration estimates with respect to cloud cover, dew formation, leaf area index distribution and mesoscale model estimates of surface climate are presented and discussed. The main conclusion is that even the limited spatially distributed hydrometeorological and plant information used in this study contributes significantly toward explaining the substantial spatial variability of potential evapotranspiration in the watershed. These results also allow the determination of key locations within the watershed where additional surface stations may be profitably placed. Copyright © 2006 John Wiley & Sons, Ltd.     

Assessment of episodic acidification in Sierra Nevada,California

Monte-Carlo simulations were used to assess the extent of shortterm alkalinity depressions occuring in Sierra Nevada lakes due to acidic deposition events. The Episodic Event Model (EEM) was used to simulate spring snowmelt events. Snow course data, precipitation data and lake acidification surveys were used to derive values for the EEM parameters. Spring snowmelt events were shown to have great impacts on the water quality of Sierran lakes. Lakes are likely to be most affected by the early-spring snowmelt event because the epilimnion depth is at a minimum, which indicates minimum dilution. Under annual average loading conditions, no Sierran lake has been reported as acidic although 29% of the lakes have alkalinities less than 40 µeq/L indicating a sensitivity to acidification. In simulations of early-spring snowmelt events, using present-day acidic loading conditions, it was estimated 79% ± 9% of the lakes would experience shortterm lake alkalinity depressions to levels less than 40 µeq/L. The results provided by the model simulations are valuable in establishing upper and lower limits on the extent of possible episodic acidification to lake-resources-at-risk. The most critical parameters controlling the magnitude of lake alkalinity depressions during snowmelt episodic events are a) the lake area to watershed area ratio — a measure of input loading, and b) the epilimnion volume — a measure of dilution and mixing.     

Delineation of reservoirs using remote sensing and their storage estimate: an example of the Yellow River basin,China

Reservoirs are an integral component of water resources planning and management. Periodic and accurate assessment of the water storage change in reservoirs is an extraordinarily important aspect for better watershed management and water resources development. In view of the shortcomings of conventional approaches in locating reservoirs' spatial location and quantifying their storage, the remote sensing technique has several advantages, either for a single reservoir or for a group of reservoirs. The satellite‐based remote sensing data, encompassing spatial, spectral and temporal attributes, can provide high‐resolution synoptic and repetitive information with short time intervals on a large scale. Using remote sensing images in conjunction with Google Earth and field check of representative reservoirs, the spatial distribution of constructed reservoirs in the Yellow River basin was delineated, and their storage volume and the residence time of the stored water were estimated. The results showed that 2816 reservoirs were extracted from the images, accounting for 89·5% of the registered total. All large‐ and medium‐sized reservoirs were extracted while small reservoirs may not be extracted due to coarse resolution and cloud‐cover shadows. An empirical relationship between the extracted water surface area and the compiled storage capacity of representative reservoirs was developed. The water storage capacity was estimated to be 66·71 km³, about 92·7% of the total storage capacity reported by the authority. Furthermore, the basin was divided into 10 sub‐basins upon which the water's residence time was analysed. The water discharge in the basin has been greatly regulated. The residence time has surged to 3·97 years in recent years, ranking the Yellow River in the top three of the list in terms of residence time and flow regulation among large river systems in the world. It is expected that it will be further extended in future owing to decreasing water discharge and increasing reservoir storage capacity. Copyright © 2011 John Wiley & Sons, Ltd.     

南方红壤区长期水土流失综合治理显著减少河流输沙——以鄱阳湖流域赣江上游平江为例

长时间序列水沙数据分析有助于科学评价流域尺度水土保持减流减沙效应,为科学开展区域水土保持成效评估提供实现路径.本文选取南方红壤区水土保持工程最为集中、持续时间最长的典型流域——平江流域,基于1975-2014年的逐日降水量、流量、含沙量数据,综合应用预置白M-K检验法、Theil-Sen趋势度估计法、Pettitt检验法与累积量斜率变化率比较法,开展河流水沙变化特征及归因分析,揭示长期水土流失治理对河流水沙的影响.结果表明:(1)1975-2014年平江年输沙量、汛期输沙量和非汛期输沙量均存在显著减少趋势,年均减少率分别为2.38×10⁴、1.75×10⁴和0.44×10⁴t/a,汛期输沙量在流域开始实施系统水土流失治理时即发生突变,而年输沙量和非汛期输沙量在系列国家水土保持重点工程实施一段时间后方开始突变;(2)流域长期水土保持显著减少平江输沙,但对径流无显著影响,重点治理期、治理后期以水土保持为主的生态保护工程对输沙量减少的贡献均在94.0%以上,且生态保护工程对汛期输沙量减少的贡献较非汛期的更为突出;(3)以水土保持为主的长期生态保护工程导致平江全年、汛期、非汛期输沙量分别减少59.3%、60.7%和55.7%;(4)在剔除大型水利工程(尤其是大中型水库)建设影响后,南方红壤区可以基于长时间序列水沙数据分析科学评价流域尺度水土保持减流减沙效应.研究结果可以为南方红壤区流域尺度水土保持减流减沙效应评价提供参考,并有望为区域系列国家生态保护与建设工程的布局和管理提供科学依据.     

Effects of afforestation on groundwater recharge and water budgets in the western region of Uruguay

Uruguay has stimulated the development of its forest sector since the promulgation of Forest Law N° 15 939 in December of 1987. Nevertheless, the substitution of natural grasslands with forest plantations for industrial use has raised concerns regarding hydrological processes of groundwater recharge and water consumption involving evapotranspiration. The purpose of this study is to assess the effects of this substitution approach on water resources. Input data were collected from two small experimental watersheds of roughly 100–200 hectares located in western Uruguay. The watersheds are characterized by Eucalyptus Globulus ssp. Maidenni and natural grasslands for cattle use. Total rainfall, stream discharge, rainfall redistribution, soil water content and groundwater level data were collected. Groundwater recharge was estimated from water table fluctuations and from groundwater contributions to base flows. Seasonal and annual water budgets were computed from October of 2006 to September of 2014 to evaluate changes in the hydrological processes. The data show a decrease in annual specific discharge of roughly 17% for mean hydrological years and no conclusive effects on annual groundwater recharge in the forested watershed relative to the reference pasture watershed. Reduced annual specific discharge is equivalent to the mean annual interception. The computed actual annual evapotranspiration is consistent with international catchment measurements. Reduction rates vary seasonally and according to accumulated rainfall and its temporary distribution. The degree of specific discharge decline is particularly high for drier autumns and winters (32 to 28%) when the corresponding rainfall varies from 275 to 400 mm. These results are of relevance for water resources management efforts, as water uses downstream can be affected. These findings, based on a study period dominated by anomalous wet springs and summers and by dry autumns and winters, oppose earlier results based on 34 years of rainfall and discharge data drawn from Uruguayan large basins. Copyright © 2016 John Wiley & Sons, Ltd.