Understanding and modeling basin hydrology: interpreting the hydrogeological signature

Basin landscapes possess an identifiable spatial structure, fashioned by climate, geology and land use, that affects their hydrologic response. This structure defines a basin's hydrogeological signature and corresponding patterns of runoff and stream chemistry. Interpreting this signature expresses a fundamental understanding of basin hydrology in terms of the dominant hydrologic components: surface, interflow and groundwater runoff. Using spatial analysis techniques, spatially distributed watershed characteristics and measurements of rainfall and runoff, we present an approach for modelling basin hydrology that integrates hydrogeological interpretation and hydrologic response unit concepts, applicable to both new and existing rainfall‐runoff models. The benefits of our modelling approach are a clearly defined distribution of dominant runoff form and behaviour, which is useful for interpreting functions of runoff in the recruitment and transport of sediment and other contaminants, and limited over‐parameterization. Our methods are illustrated in a case study focused on four watersheds (24 to 50 km²) draining the southern coast of California for the period October 1988 though to September 2002. Based on our hydrogeological interpretation, we present a new rainfall‐runoff model developed to simulate both surface and subsurface runoff, where surface runoff is from either urban or rural surfaces and subsurface runoff is either interflow from steep shallow soils or groundwater from bedrock and coarse‐textured fan deposits. Our assertions and model results are supported using streamflow data from seven US Geological Survey stream gauges and measured stream silica concentrations from two Santa Barbara Channel–Long Term Ecological Research Project sampling sites. Copyright © 2004 John Wiley & Sons, Ltd.     

中国冰冻圈水文过程变化研究新进展

冰冻圈显著的变化已经对冰冻圈水文过程产生了一系列影响。本文重点梳理和分析了近20年，尤其是近10年以冰川融水、融雪径流、冻土水文等为主体的中国冰冻圈水文过程变化研究方面取得的新进展：①在冰川融水变化研究方面，对不同尺度的冰川融水开展了全面研究，发现冰川融水呈现全面增加之势；对冰川融水"拐点"是否出现进行了科学辨识，有了基本认识；对冰川融水过程进行了模型模拟，取得显著进展。②在融雪径流变化研究方面，通过对不同流域融雪径流估算，可基本掌握各河流的融雪贡献率；中国融雪径流变化差异较大，增减不一；融雪期变化具有普遍性，突出特点是峰值提前。③在冻土水文研究方面，通过对地表水-活动层壤中流-多年冻土层上水之间关系的研究，揭示了冻土区径流形成的重力和热力耦合机制；多年冻土变化对地表径流的影响已经显现，主要表现在冬季（枯水季）径流增加；已经发现多年冻土退化对径流有直接补给作用，在一些流域补给量可能达到一定量级。     

Areal extent of the 1988 flood in Bangladesh: How much did the satellite imagery show?

This paper presents the methods and results of visual interpretation of satellite imagery for estimating areal extent of the 1988 flood in Bangladesh. The main limitation of this method was that cloud-free images were not available for the peak flood period in the first week of September. Relatively cloud-free images on three selected dates in the last three weeks of September showed that areas of inundation on those dates ranged from 44000 to 60000 km² ± 5% (31 to 42% of Bangladesh). The Government of Bangladesh official estimates of the areal extent of flooding were tentative in nature, since these were based on a choropleth map of relative proportions (percentages) of flooding in different parts of the country. The official estimate of the maximum extent of flooding was 82000 km² (57% of Bangladesh), i.e. 15 to 26% larger than the area shown on satellite imagery. The actual extent of flooding was certainly larger than the estimates from satellite imagery but, perhaps, smaller than the official estimates, since a reduction of 15 to 26% of flood area in one to three weeks seemed unlikely. An alternative method of mapping flood-affected areas by using newspaper-interpreted data was attempted, but the method had limited values because of reporting bias.     

Effects of soil moisture dynamics on slope failure at Hyrum Reservior,Utah

Field observations of shoreline conditions at Hyrum Reservoir, Utah, were conducted during the summers of 1991 to 1993. A process of bluff retreat is described for a multiple-layered bluff environment of sand and clay layers. Failure is initiated by wetting and drying of clay sediments, which produces horizontal cracks within bluff material. These cracks appear to penetrate to a depth of approximately 100-150 mm before initiating vertical cracking in the sediments. The vertical cracks are propagated by continued drying of the surface sediment, ultimately leading to failure of the bluff material. The physical dimensions of sediment blocks succumbing to this mechanism range from a few hundred millimetres up to 3 m on a side, with a depth of approximately 100-150 mm. The mechanism described here appears to operate optimally when the supply of subsurface moisture is abundant and nearly continuous throughout the spring and early summer. Reservoir draw-down, large capillary fringe effects in the bluff and periodic wetting from upslope undrained hollows are the dominant moisture controls at this site. Moisture delivery to the face is strongly influenced by anisotropy of saturated hydraulic conductivity in the alternating clay and sand layers and related differences in sediment texture.     

Changing drainage patterns in a semi-arid area and their effects on groundwater resources

The magnitude of groundwater resources is poorly correlated with surface catchment area in a group of wadi bottom alluvial aquifers in Northern Oman. It is suggested that this is due to progressive water gap captures across a limestone ridge.     

Identification of buried sinkholes using refraction tomography at Ft. Campbell Army Airfield,Kentucky

Karst aquifers are highly susceptible to contamination, with numerous points of entry for contaminants through recharge features such as sinkholes, swallow holes and solutionally enlarged fractures. These recharge features may be filled or obscured at the surface, requiring the use of geophysical or remote sensing techniques for their identification. This study uses seismic refraction data collected at the Ft. Campbell Army Airfield (CAAF), Kentucky, USA, to test the hypothesis that refraction tomography is a useful tool for imaging bedrock depressions beneath thick overburden (greater than 20 m of unconsolidated sediment). Southeast of the main taxiway of CAAF seismic velocity tomograms imaged a bedrock low, possibly a closed depression, at a depth of 25 m that had been earlier identified through delay-time analysis of the same refraction data. Tomography suggests the bedrock low is about 250-m wide by 10-m deep at its widest point. High rates of contaminant vapor extraction over the western extension of this feature suggest a high concentration of contaminants above, and within, this filled bedrock low, the base of which may contain solutionally enlarged fractures (i.e. karst conduits) that could funnel these contaminants to the upper or lower bedrock aquifers. This study thus demonstrates the viability of seismic refraction tomography as a tool for identification of filled sinkholes and bedrock depressions in karst areas.     

Stochastic sketching of infiltration-advance isochrones

A technique has been developed for predicting the irregular advance pattern often observed as water spreads on the surface of the ground. The technique is a combination of stochastic sketching, potential theory, probability theory, and a mass balance equation in the form of an advance equation. The technique can be used on flat as well as sloping terrain and addresses any form of obstructions or constraints to the flow of the water. The stochastic sketching portion of the technique uses cellular automata with transition probability movement rules to sketch the dynamics of small volume water elements in the defined environment. Randomly selected small volume flow path segments are computed and plotted. The envelope of these segments defines the wetted area and the advance front. Several examples are presented showing the patterns produced for various situations.     

Stochastic perturbation analysis of groundwater flow. Spatially variable soils,semi-infinite domains and large fluctuations

As is well known, a complete stochastic solution of the stochastic differential equation governing saturated groundwater flow leads to an infinite hierarchy of equations in terms of higher-order moments. Perturbation techniques are commonly used to close this hierarchy, using power-series expansions. These methods are applied by truncating the series after a finite number of terms, and products of random gradients of conductivity and head potential are neglected. Uncertainty regarding the number or terms required to yield a sufficiently accurate result is a significant drawback with the application of power series-based perturbation methods for such problems. Low-order series truncation may be incapable of representing fundamental characteristics of flow and can lead to physically unreasonable and inaccurate solutions of the stochastic flow equation. To support this argument, one-dimensional, steady-state, saturated groundwater flow is examined, for the case of a spatially distributed hydraulic conductivity field. An ordinary power-series perturbation method is used to approximate the mean head, using second-order statistics to characterize the conductivity field. Then an interactive perturbation approach is introduced, which yields improved results compared to low-order, power-series perturbation methods for situations where strong interactions exist between terms in such approximations. The interactive perturbation concept is further developed using Feynman-type diagrams and graph theory, which reduce the original stochastic flow problem to a closed set of equations for the mean and the covariance functions. Both theoretical and practical advantages of diagrammatic solutions are discussed; these include the study of bounded domains and large fluctuations.