Showing papers in "Hydrological Processes in 1988"

Groundwater recharge in arid regions: Review and critique of estimation methods

Abstract: Arid-site recharge, while generally low, can be highly variable. Recharge under similar climate and soil conditions but with different plant cover and topography can vary from zero to more than the annual precipitation. Simple estimates of recharge based on fixed fractions of annual precipitation are misleading because they do not reflect the plant and soil factors controlling recharge. Detailed water balance models, successful for irrigated agriculture, fail to predict evapotranspiration accurately under conditions where plants suffer seasonal water stress and cover is sparse. Recharge, when estimated as a residual in water balance models, may be in error by as much as an order of magnitude. Similar errors can occur when soil water flow models are used with measured or estimated soil hydraulic conductivities and tension gradients. Lysimetry and tracer tests offer the best hope for evaluating recharge at arid sites, particularly in siting waste disposal facilities, where reliable recharge estimates are needed. Quantification of drainage using lysimetry over several years under a given set of soil, plant, and climate conditions for a specific site can provide a basis for calibrating models for recharge prediction. Tracer tests using such long-lived tracers as 36Cl or perhaps stable isotopes (180, deuterium) can provide qualitative estimates of recent recharge at a given site.

Chloride in precipitation and streamwater for the upland catchment of river severn, mid‐wales; some consequences for hydrochemical models

Abstract: Variations in the concentration of Cl in rainfall and stream runoff are presented for two catchments in the Hafren forest of mid-Wales, Great Britain. Despite the large fluctuations in rainfall concentrations, Cl in the streamwater remains relatively constant. Using the two-reservoir Birkenes model, an attempt was made to simulate observed Cl in streamwater. The original model was unable to reproduce the observations and several modifications are suggested to provide better simulations. The resulting model is not the only one capable of reproducing the observations; other hydrochemical models will most probably also achieve this although emphasis will in each case be placed on different aspects. In this paper, it is suggested that the stochastic properties of water movement and chemical processes can account for the streamwater chemistry responses observed. On the catchment scale these processes will lead to an apparently deterministic behaviour that may well be described by simple relationships.

Stream nitrate levels in a small catchment in South West England over a period of 15 years (1970-1985)

Abstract: Stream nitrate levels in a small catchment of mixed land use (the Slapton Wood catchment) have been studied since September 1970; a record of this length is possibly unique in the United Kingdom for such a small basin (94 ha). A sustained increase in nitrate concentration has been observed during the study period. In addition to this long-term trend, short-term changes in nitrate concentrations relate to stream discharge levels and to seasonal variations. Multivariate statistical analysis has been used to quantify these trends and to identify those factors controlling the production and loss of nitrate from the catchment system. The main period of nitrate removal occurs in winter when high concentrations coincide with the main period of throughflow generation. The influence of climatic variability is illustrated by reference to the 1975/76 drought and post-drought period.

Sources and sinks of nutrients in a New Zealand hill pasture catchment II. Phosphorus

Abstract: The management of the riparian zone has been suggested as a technique for controlling the amounts of phosphorus (P) entering watercourses draining pasture catchments. A study was therefore made of P entering a stream from various sources (rainfall, surface and subsurface derived runoff, direct fallout from aerial topdressing), with the object of providing a rational basis for the design of effective riparian management schemes. P entrained in surface runoff could account for virtually all of the P entering the stream during storms. Approximately 20 per cent of the annual P export from the catchment could be accounted for by direct aerial input of P to the stream during autumn fertilizer topdressing. More than 85 per cent of the P was exported from the catchment as particulate P. Stream sediment had higher P sorption capacities, and were enriched with P relative to the soils from which they were derived. There was a pronounced seasonal variation in sediment enrichment which could be predicted (r2 = 0.92) by the logarithm of the rainfall since fertilizer topdressing (LNFERT) and flood intensity. The amount of P lost in streamflow during any flood event was predicted (r2 = 0.94) by peak flow, seven day antecedent peak flow and LNFERT. Approximately 40 per cent of the 1.3 kg P ha−1 exported during 1981 occurred in four storms with recurrence intervals of more than three months. From a P budget compiled from nine events it was hypothesized that the stream acted as a net sink for P at baseflow and low-medium intensity floods but was a source of P at higher flood intensities. It was concluded that P losses from hill pasture catchments could be reduced by avoidance of direct application of P fertilizer to the stream channel, and by fencing out stock from seasonally saturated areas during periods of saturation. The ultimate success of the latter technique would depend on whether buffer vegetation could retain accumulated P during extreme storm events.

Interception loss from eucalypt forest: Lysimeter determination of hourly rates for long term evaluation

Abstract: Analyses of the response by a weighing lysimeter in Kioloa State Forest during and after rainfall provided values of interception loss rate. The derived rates for time scales between 0.1 and 1.0mm h−1 were generally similar throughout storm events to losses determined from throughfall and stemflow observations. During post-rainfall periods of canopy drying, enhanced rates of lysimeter evaporation were consistent with micrometeorological determinations of the partitioning of available radiant energy, based on atmospheric gradients of humidity and temperature. Interception losses from the eucalypt forest, deduced from the lysimeter response, varied between 10 and 15 per cent of gross rainfall in three consecutive 12 month periods whereas the corresponding rainfall ranged between 590 and 1530 mm yr−1. Daytime losses accounted for about two-thirds of total interception loss with a similar fraction occurring during rain periods. Storage capacity of the evergreen forest canopy was inferred to be 0.35 mm. Hourly loss rates during rainfall ranged up to 0.8 mm h−1 but with decreasing mean values and variability with increasing time scale resulting in a monthly mean value computed for the total number of hours of rain of approximately 0.1 mm h−1. A preliminary analysis of loss rate in terms of storm windspeed and rainfall intensity explained about half of its variation in statistically derived relationships. Improved time resolution of the order of seconds was considered a prerequisite to the physical understanding of turbulent transport from saturated canopies. The small value of interception storage capacity was considered in relation to that for pine forest as a basis for explaining observed differences in interception behaviour between eucalypt forest and coniferous plantations in the same area. Large differences in interception losses between the Kioloa site and evergreen forest in the South Island of New Zealand and also eucalypt forest in Western Australia were attributed to dissimilar meteorological conditions at the various sites.

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

Abstract: 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.

A model of river flow forecasting for a small forested mountain catchment

Abstract: A model for quantitatively expressing the hydrological cycle in a forested mountain catchment is proposed as a HYCYMODEL. HYCYMODEL is able to predict both short- and long-term hydrographs because the model parameters remain independent of time. It shows a good applicability for ten years of continuous data at both hourly and daily intervals for the Kiryu catchment—a forested mountain basin. Since HYCYMODEL does not need hydrograph separation between storm flow and base flow, it is a particularly attractive model.

Sources and sinks of nutrients in a New Zealand hill pasture catchment. III: Nitrogen

Abstract: A study was made of the nitrogen (N) inputs to, and exports from, a stream draining a pasture catchment near Hamilton, New Zealand, in order to plan measures for minimizing N losses to natural waters. An estimated 7 kg N ha-' was exported from the catchment during 1981 of which 86 per cent was in reduced forms (Kjeldahl-N, TKN) and the remainder as nitrate-N (N(3-N). Virtually all of the reduced N inputs came from saturated overland flow whereas N03-N inputs were dominantly subsurface derived. Thc TKN exported by individual storm events could be predicted (R' = 0.97) from peak flow and from the peak flow rate in the seven days preceding the storm. A TKN balance for eight events showed that except for large floods (return period approximately a year) the stream system was a net sink for TKN. During large floods, scouring of the organic rich seepage areas resulted in the stream system itself being a net source of TKN. Microbial assays for nitrification and denitrification activity indicated that the main nitrate source was the well-aerated greywacke and ash soils and that the permanently saturated seepage zones were a significant nitrate sink. An in-stream nitrate addition experiment showed that up to 30mg N m-: h-' was removed from the stream. Simultaneous measurements of in situ denitrification activity demonstrated that only about 1 per cent of this removal could be accounted for by denitrification. It was inferred that plant uptake was responsible for the remainder. Retention of near-stream seepage areas is suggested as a measure for minimizing N03-N export, whilst removal of stock from seasonally saturated areas during periods of saturation should reduce soil loss and hence TKN inputs to the stream.

Modelling long term stream acidification trends in upland wales at plynlimon

Abstract: Historical reconstructions and predictions of streamwater acidification are presented for moorland and afforested catchments in the Welsh Uplands at Plynlimon. MAGIC (Model of Acidification of Groundwater In Catchments) is validated by application to an afforested catchment. MAGIC is used to illustrate that atmospheric deposition is primarily responsible for stream acidification but that conifer afforestation can enhance stream acidity. The historical trends determined by the model illustrate that long-term acidification has been present since the turn of the century and will continue unless either deposition levels are reduced significantly or liming is undertaken on a major scale.

Modelling pipeflow contributions to stream runoff

Abstract: Runoff from natural soil pipes has been shown to be a significant contributor to stream discharge in parts of upland Wales. Attempts have been made to model pipeflow contributions using both theoretical and empirical approaches, but little progress has yet been made towards producing generally applicable models of complete pipeflow systems. The paper identifies some of the problems of devising general models with reference to data from the Maesnant Experimental Catchment. The data suggest that in-pipe travel times are more important here than inferred elsewhere, and that neither simple hydraulic models nor kinematic wave theory can adequately explain the patterns of response. In particular, stormflow tends to begin first at the network outfall and even peak discharge often occurs at the outfalls before the headward zones.

One year irrigation experiment to assess losses and runoff volume relationships for a residential road in hertfordshire, England

Abstract: Six sections of a residential road (75 mm bituminous macadam over 200mm lean mix concrete over 100 mm hoggin with a recent surface dressing of 10 mm granite chippings and K1-70 binder) that drain to individual instrumented gully pots were irrigated along the kerb and then over the whole road approximately monthly for a year. The aim was the determination of terminal infiltration losses, initial losses, percentage runoff, and infiltration curves for the kerb and road surface. The results were not as expected from the literature. There is an annual cycle of infiltration losses at the kerb with a winter peak caused by frost action that is 3.2 times greater than the terminal loss rate at kerbs in summer. The terminal loss rate for an ‘average catchment’ was 6.4251 min−1 from the road surface and 14.251 min−1 in summer and 46.281 min−1 in winter at the kerb. Evaporation was usually more than an order of magnitude less significant than infiltration. The wide variation in initial losses before runoff commenced was inexplicable. Two sections of road behaved in the classic manner with initial losses averaging 0.8 mm, two other catchments had highly variable initial losses in the range 1.2 to 8.8mm, and the last two pieces of road were even more erratic. The percentage runoff for those irrigations of over 15 mm of equivalent rainfall was never more than 10 per cent. The maximum per cent runoff was around 50 per cent following 10 mm of equivalent rainfall for kerbside irrigation and only 5 mm of irrigation over the whole road. There were no significant simple or multiple regression relationships between percentage runoff from the kerb or the whole road irrigations and irrigation amount, slope, UCWI, and SMD. Infiltration curves, for kerb and road irrigation, were so diverse that they do not represent the ‘simple impervious surfaces’ envisaged at the start of the experiment and described in the literature. Since the artificial irrigation of kerbs and roads has failed to substantiate existing theory, these experiments should be repeated at a variety of sites with a high rate of irrigation.

Soil saturation patterns in steep, convergent hill-slopes under forest and pasture vegetation

Abstract: The relationships between hillslope form and the development of zones of hillslope saturation are central to studies of process hydrology. These relationships are examined in two small, steep (30°) first-order basins under forest and pasture vegetation. Hillslope saturation patterns are described using a probability analysis of observations of peak water table elevations made at individual sites in each basin. In both first-order basins, the effects of topographic convergence on the spatial variability of hillslope saturation zones are most pronounced on shallow slopes and in steeper areas where hillslope form is strongly concave. The implications of vegetation type on the formation and variation of hillslope saturation are also examined. The differences in hillslope saturation patterns in each vegetation type are attributed primarily to differences in net rainfall and not soil physical properties or other features of topography.

Sources and sinks of nutrients in a New Zealand hill pasture catchment I. Stormflow generation

Abstract: The processes of stormflow generation were studied in a hill pasture catchment near Hamilton, New Zealand. Although rainfall was relatively evenly distributed throughout the year, stormflow was highly seasonal and over 65 per cent occurred during the winter. Three main processes contributing to stormflow were identified which could be related to soil type and physiographic position. On gleyed soils derived from rhyolitic colluvium, saturation overland flow was the dominant process. Hydrographs from ‘Whipkey’ throughflow troughs also indicated that there was a subsurface response (saturated wedge) from this soil type. On steeper convex slopes, more permeable soils were derived from weathered greywacke. The presence of ephemeral springs on the hillslopes and direct observation during storm events indicated that storm runoff was generated as return flow from this soil. It was noted that nitrate concentrations from subsurface sources were 5–10 times higher than surface runoff. This difference in concentration was utilized in a chemical mixing equation which partitioned stormflow sources. This was compared with the stormflow predicted from rain falling on to saturated areas. There was good agreement between the two models for winter-spring events with respect to the volumes of surface runoff predicted, however the saturated areas model underestimated total stormflow. The results of the study are briefly discussed in terms of the potential for water quality management.

Field application of a numerical method for the derivation of baseflow recession constant

Abstract: In an earlier paper (Bako and Hunt, 1988), a method for the derivation of the baseflow recession constant (K) using one-way analysis of variance was presented. This paper presents the results of the field application of this method. The K values obtained by using the numerical equation of Bako and Hunt (1988) were inserted in the exponential recession equation (Barnes, 1939) to generate a series of baseflows. The fit between the model and the historical flows was found to be greater than 99 per cent thus confirming the applicability of the numerical method under field situation. The main advantage of this technique is its amenability for computerized application thus making it relatively faster than any of the existing techniques of fitting the recession equation. For this reason, the subjectivity inherent in most of the existing techniques is eliminated and a measure of procedural consistency can be guaranteed. Consistency is necessary if intercatchment comparison or interpolation of K values is to be meaningful.

On numerical errors associated with the iterative alternating direction implicit (Iadi) finite difference solution of the two dimensional transient saturated-unsaturated flow (Richards) equation

Abstract: Numerical models for simulation of multidimensional unsaturated flow are becoming increasingly available, but relatively little has been reported on the detailed analysis of numerical errors associated with such schemes. For unsaturated-saturated flow, further complexity is introduced to the highly non-linear unsaturated problem as the form of governing equation changes within the flow field. In this paper, two-dimensional simulation of infiltration to a water table is considered. Numerical errors associated with selection of time step, grid geometry, convergence criteria, and the representation of internodal hydraulic conductivity are discussed with respect to moisture content profiles and mass balance error. Although solution sensitivity to numerical parameters is problem specific, the results presented indicate the nature and magnitude of numerical effects which should not be overlooked in model applications.

The physical framework of the dependence between channel flow hydrographs and drainage network morphometry

Abstract: The physical basis of the linkage between magnitude and timing of channel flow hydrographs and drainage network morphometry is reviewed. Small Hortonian and structurally Hortonian networks are analysed using numerical runoff simulation. For Hortonian networks the variability of the geometry of individual channels and subcatchments within each Strahler order has generally little effect upon the overall character of the hydrograph in channels of higher order. If the network is also structurally Hortonian, the analysis of the simultaneous formation, travel, and concentration of the hydrographs in all channels of the network can be simplified to a sequence of one representative hydrograph per channel order. This approach is used in this study. Three major runoff processes control the flow hydrograph characteristics: the overland flow process which determines the water supply to the drainage network; the channel flow process which translates the hydrograph in space and time; and the drainage network process which concentrates and magnifies the flow at the junctions of the drainage network. Functional relations for the hydrograph peak, timing, and flow velocity are presented. For a given uniform rainfall and infiltration rate, the peak of the channel flow hydrograph is shown to increase geometrically with channel order, and its magnitude is directly related to the bifurcation ratio. The travel time of the peak also increases geometrically with channel order, and it is directly related to the channel length ratio over velocity ratio. The flow velocity of the peak changes in a downstream direction as a function of the bifurcation and slope ratio. It was also found that for negligible channel storage the channel flow and drainage network processes do not contribute significantly to the observed nonlinear response of a watershed to precipitation.

The mechanism of unsaturated flow through a volcanic ash layer under humid climatic conditions

Abstract: It is very interesting and meaningful to investigate the rainfall-groundwater recharge process under the humid climatic condition of Japan, where mean annual precipitation is about 1600 mm. The present study has investigated soil water movement in the unsaturated zones of a volcanic ash layer, called the ‘Kanto Loam formation’, using environmental tritium as a tracer. The site selected is a flat ground surface on a terraced upland which has a deep unsaturated zone (about 20 m) with a relatively high water content (about 70 per cent) consisting of nearly uniform Kanto Loam formation. The tritium concentrations in groundwater, soil waters having different matric potentials, precipitation, and the seepage water moving through the formation into a man-made cave were measured to characterize the rainfall-groundwater recharge process and the effect of large pore spaces in the formation mentioned by previous studies. Because of the humid climate of Japan, there appears to be a unique soil water flow characteristic which may involve percolation through large pore spaces during heavy rainfall. However, in a fine grained and high water content soil like the Kanto Loam formation, the existence of this flow through large pore spaces does not have a significant effect upon the whole recharge process. The recharge model of displacement flow with dispersion is useful in estimating the tritium concentration profile of soil water. The calculated result shows a recharge rate of 2.5 mm/day. The value obtained reflects the hydrological characteristics of the uplands covered with volcanic ash.

Some aspects of the hydrology of one‐day annual minimum low flows of malawian rivers

Abstract: Five aspects of the hydrology of one-day annual minimum flows QIM, have been studied using data from twelve catchments in Malawi. Results indicate that the log-normal distribution can be fitted to all twelve catchments. Four of the rivers studied are intermittent. Application of statistical methods developed in meterology to the dichotomous-transformed data of these catchments revealed that two are ‘flow-dominant’ and the other two are ‘dry-dominant’. Another catchment is entirely dominated by a hydraulic gradient towards the Shire River and Elphant Marsh and so dries up every dry season for considerable periods of time despite the relatively high rainfall in the catchment. QIM, t-days after the date of occurrence of QIM(May), can be better estimated from simple regression than from an empirically determined recession constant.

Interactions between groundwater seepage and sediment porewater sulphate concentration profiles in lake anna, Virginia

Abstract: Most sulphur diagenesis models predict that SO42- concentrations decrease exponentially with increasing sediment depth and are lower than that of the overlying water throughout the sediments. Low SO42- concentrations (less than 0.2 mM) are common in the sediments of Lake Anna that receive acid mine drainage; however, sediment with as much as 20 mM SO42- at about 20cm below the sediment surface is also seen in this section of the lake. A decision tree was proposed to investigate the cause of the high SO42- concentrations at depth (HSD) in the sediment. The first possibility proposed was that an increase in the quantity of groundwater flowing through Lake Anna sediments may increase groundwater advection of SO42- or oxygen which would induce sulphide oxidation. This hypothesis was tested by measuring groundwater flow. HSD profiles were found in a discrete region of the lake; however, stations having these profiles did not have higher groundwater flow than other sites sampled. Alternate explanations for the HSD profiles were that the region in which they occurred had: (1) unusual sediment chemical compositions; (2) a different source of regional groundwater, or (3) a lateral intrusion of high SO42- groundwater. There were no differences in sulphide and organic matter concentrations between the two regions. The area which has HSD in the sediment covers a large area in the middle of the lake, so it is unlikely that it has a unique source of regional groundwater. The third alternative was supported by the fact that in all three sample years, HSD stations were located in the preimpoundment stream channel, which is a likely lateral flow path for groundwater containing high SO42- concentrations.

Effect of groundwater interception and irrigation on salinity and piezometric levels of an aquifer

Abstract: The major strategy used to prevent the discharge of highly saline groundwater to the River Murray in southeastern Australia is groundwater interception and disposal. The basic design principle assumes that the extraction of groundwater from an aquifer hydraulically connected to the river, using a line of pumps positioned close and roughly parallel to the river, will decrease piezometric heads thereby reducing the discharge of saline groundwater to the river. The paper considers one of these schemes which was designed for the Mildura area on the basis of a hydrogeological investigation. It analyses the effects on piezometric head and groundwater salinity due to the groundwater interception scheme and adjacent irrigation activity over a period of several years from January 1980. It is shown that piezometric heads have decreased significantly in the stretch close to the river. A slight reduction in groundwater salinity is also apparent in this stretch except for an area between the river and a holding basin used for disposal of the saline effluents emanating from the groundwater interception scheme. This general reduction in groundwater salinity is mainly caused by pumping from the groundwater interception scheme and recharge from irrigation. The exception in the trend in groundwater salinity is due to the movement of a highly saline body of groundwater from the holding basin towards the River Murray. Results of this Australian experience should be helpful to the designers of similar salinity mitigation schemes elsewhere.

Modelling rainfall/watertable relations using linear response functions

Abstract: Linear system theory can be used to model and predict watertable responses to precipitation inputs in an artificially drained field. The response function is mathematically equivalent to the Unit Hydrograph concept familiar to hydrologists. This paper shows that it is possible to derive such response functions, and comments on the problems encountered in their derivation. Response functions for two contrasted sites are presented, and the possibilities for their use discussed.

A weighing lysimeter in a regenerating eucalypt forest: Design, construction, and performance

Abstract: A 36 tonne monolithic weighing lysimeter (3.7 m diameter and 1.5 m deep) was installed in Kioloa State Forest near Batemans Bay, New South Wales, Australia, to provide a continuous record of water use from a regenerating natural eucaplyt community, five years old at installation. The resolution is equivalent to 0.05 mm of evaporation and provides hourly information for diurnal variation of forest evaporation. The sensitivity of 14 ppm compares favourably with those reported for other lysimeters, similar in magnitude, in the U.S.A. Design criteria for selecting, isolating and suspending a sample community ‘in situ’ with an electro-mechanical balance are outlined. Procedures adopted during installation are described with emphasis placed on precautions taken to minimize disturbance and to avert possible changes in character of the isolated natural community. Lysimeter performance was assessed both in terms of its accuracy and its operation as a remote facility with infrequent attention.

Drainage basin peak discharge rating curve

Abstract: The stream gauge rating curve for a drainage basin can be transformed into a drainage basin peak discharge rating curve that is more stable than the rating curve from which it is derived. The resulting drainage basin peak discharge rating curve can be used to predict peak discharge, identify anomalous discharges caused by channel obstructions or other causes, evaluate the effect of flood retarding structures, and evaluate historical records. The drainage basin peak discharge rating curve is valid for drainage basins of any size, for any discharge up to the time of concentration, and for snowmelt.

Three instruments used in rainfall—runoff simulation experiments

Abstract: Three instruments recently developed for use in the Runoff Formation Laboratory of the Institute of Geography, Academia Sinica are briefly described. The experimental data obtained indicate that all of them meet the design requirements and they are being successfully used in the Laboratory.

Finite element computation of shape factors for Kirkham's piezometer tube method for hydraulic conductivity of saturated soil

Abstract: The Kirkham piezometer tube method of measuring hydraulic conductivity of saturated soil uses a numerical constant (shape factor) to account for the relative geometrical parameters of the cavity and the soil domain. Published shape factors have mostly been determined using an electrical analogue. In this paper the application of the finite element method to shape factor determination is considered. Convergence is studied and methods of extrapolation evaluated. A determination is made of the size of finite region required to represent adequately the radially unbounded domain assumed for the piezometer method. Computed shape factors were found to be in good agreement with analogue data.

Exceedance probability distributions of steady conduit flow in karst aquifers

Abstract: Discharge exceedance probabilities are calculated for a simple model karst aquifer composed of a few multilevel conduits with recharge from a single sinking stream with an exponential flow exceedance distribution. It is assumed that outflow instantaneously matches inflow, so that the conduit volume is constant but discharge is governed by the head in a volumeless shaft at the top end of the system. It is shown that small single conduit aquifers will frequently overflow at the surface during floods and the exceedance probability of flow through the aquifer and over the surface can be defined as a function of the inflow distribution and the form of the aquifer. Systems with multiple conduits will overflow less frequently, but each conduit will exhibit a flow distribution characteristic of its form and position in the vertical hierarchy. Comparison of these findings with actual flow data from a conduit aquifer shows that the approach is valid, although imprecise. The model is unlikely to be applied directly, as it requires unusually detailed data. However, it provides for the first time an indication of typical flow durations for surface overflows and individual conduits in a karst aquifer. Contrasts in flow duration will have a profound influence on solutional and sedimentary processes in the karst system.

Solute concentrations in the lower niger river and the source rock contribution

Abstract: Water quality analyses for the Niger River for the 1980/81 hydrological year are presented. The samples were collected from the main river at Lokoja, and from two main tributaries, the Kaduna and the Benue Rivers. Different water types were distinguished by the concentrations of major ions. The type Ca > Na > Mg > K - HCO3 > SO4 > Cl was represented at all stations during at least part of the year. Chloride was found to dominate the sulphate ion in the Kaduna and Niger, while the Benue maintained a higher concentration of sulphate relative to chloride all year round. Distinct patterns of seasonal variation in the ion concentrations were observed, particularly for the samples collected at Lokoja. Low ion concentrations were prominent during periods of high discharge, while low flow periods coincided with high dissolved ion concentrations. The contribution of rainwater to the total dissolved solids in the river waters was assessed indirectly using rainwater chemistry data from the Gulf of Guinea. The estimated rainwater contribution to the Lower Niger amounts to 5.15 mg 1−1. Geochemical weathering calculations involving reactions of the four major minerals of granitic rocks - anorthite, biotite, albite, and K-feldspar - with carbon dioxide and water, can account for the average water composition of the Lower Niger. The proportion of the ionic components was also related to the occurrence of the respective element in the minerals.