Hydraulic properties of groundwater systems in the saprolite and sediments of the wheatbelt, Western Australia

Hydraulic properties of deeply weathered basement rocks and variably weathered sedimentary materials were measured by pumping and slug-test methods. Results from over 200 bores in 13 catchments, and eight pumping-test sites across the eastern and central wheatbelt of Western Australia were analysed. Measurements were made in each of the major lithological units, and emphasis placed on a ubiquitous basal saprolite aquifer. Comparisons were made between alternative drilling and analytical procedures to determine the most appropriate methods of investigation.

Aquifers with an average hydraulic conductivity of 0.55 m day⁻¹ occur in variably weathered Cainozoic sediments and poorly weathered saprolite grits (0.57 m day⁻¹). These aquifers are separated by an aquitard (0.065 m day⁻¹) comprising the mottled and pallid zones of the deeply weathered profile. Locally higher values of hydraulic conductivity occur in the saprolite aquifer, although after prolonged periods of pumping the values decrease until they are similar to those obtained from the slug-test methods. Hydraulic conductivities measured in bores drilled with rotary auger rigs were approximately an order of magnitude lower than those measured in the same material with bores drilled by the rotary air-blast method.

Wheatbelt aquifers range from predominantly unconfined (Cainozoic sediments), to confined (saprolite grit aquifer). The poorly weathered saprolite grit aquifer has moderate to high transmissivities (4–50 m² day⁻¹) and is capable of producing from less than 5 to over 230 kl day⁻¹ of ground water, which is often of a quality suitable for livestock. Yields are influenced by the variability in the permeability of isovolumetrically weathered materials from which the aquifer is derived.

The overlying aquitard has a low transmissivity (< 1 m² day⁻¹), especially when deeply weathered, indurated and silicified. The transmissivity of the variably weathered sedimentary materials ranges from less than 0.5 m² day⁻¹ to over 10 m² day⁻¹, depending on the texture of the materials and their position within the landscape. Higher transmissivity zones may occur as discrete layers of coarser textured materials. The salinity of the saprolite and sedimentary aquifers ranges from less than 2000 mgl⁻¹ to greater than 250000 mgl⁻¹ (total dissolved solids; TDS), depending on position within the landscape. Secondary soil salinization develops when groundwater discharge occurs from either saprolite or sedimentary aquifers.     

Water budget and vertical conductance for Lowry (Sand Hill) Lake in north-central Florida, USA

Water-budget components and the vertical conductance were determined for Lowry (Sand Hill) Lake in north-central Florida, USA. In this type of lake, which interacts with both the surface-water and groundwater systems, the inflow components are precipitation, surface-water inflow, groundwater inflow, and direct runoff (i.e. overland flow), and the outflow components are evaporation, groundwater outflow, and surface-water outflow. In a lake and groundwater system that is typical of many karst lakes in Florida, a large part of the groundwater outflow occurs by means of vertical leakage through an underlying confining unit to a deeper, highly transmissive aquifer called the upper Floridan aquifer. The water-budget component that represents vertical leakage to the upper Floridan aquifer was calculated as a residual using the water-budget equation. For the 13 month period from August 1994 to August 1995, relative to the surface area of the lake, rainfall at Lowry Lake was 1.55 m yr⁻¹, surficial aquifer inflow was 0.79 m yr⁻¹, surface-water inflow was 1.92 m yr⁻¹, and direct runoff was 0.01 m yr⁻¹. Lake evaporation was 1.11 m yr⁻¹, and surface-water outflow was 1.61 m yr⁻¹. The lake stage increased 0.07 m yr⁻¹, and the vertical leakage to the upper Floridan aquifer was 1.48 m yr⁻¹. Surficial aquifer outflow from the lake was negligible. At Lowry Lake, vertical leakage is a major component of the water budget, comprising about 35% of the outflow during the study period. The vertical conductance (K_V/b), a coefficient that represents the average of the vertical conductances of the hydrogeologic units between the bottom of a lake and the top of the upper Floridan aquifer, was determined to be 2.51 × 10⁻⁴ day⁻¹ for Lowry Lake.     

Hydrochemical and stable isotope evidence for the extent and nature of the effective Chalk aquifer of north Norfolk, UK

In eastern England the Chalk aquifer is covered by extensive Pleistocene deposits which influence the hydraulic conditions and hydrochemical nature of the underlying aquifer. In this study, the results of geophysical borehole logging of groundwater temperature and electrical conductivity and depth sampling for major ion concentrations and stable isotope compositions (δ¹⁸O and δ²H) are interpreted to reveal the extent and nature of the effective Chalk aquifer of north Norfolk. It is found that the Chalk aquifer can be divided into an upper region of fresh groundwater, with a Cl concentration of typically less than 100 mg l⁻¹, and a lower region of increasingly saline water. The transition between the two regions is approximately 50 m below sea-level, and results in an effective aquifer thickness of 50–60 m in the west of the area, but less than 25 m where the Eocene London Clay boundary is met in the east of the area. Hydrochemical variations in the effective aquifer are related to different hydraulic conditions developed in the Chalk. Where the Chalk is confined by low-permeability Chalky Boulder Clay, isotopically depleted groundwater (δ¹⁸O less than −7.5‰) is present, in contrast to those areas of unconfined Chalk where glacial deposits are thin or absent (δ¹⁸O about −7.0‰). The isotopically depleted groundwater is evidence for groundwater recharge during the late Pleistocene under conditions when mean surface air temperatures are estimated to have been 4.5°C cooler than at the present day, and suggests long groundwater residence times in the confined aquifer. Elevated molar Mg:Ca ratios of more than 0.2 resulting from progressive rock-water interaction in the confined aquifer also indicate long residence times. A conceptual hydrochemical model for the present situation proposes that isotopically depleted groundwater, occupying areas where confined groundwater dates from the late Pleistocene, is being slowly modified by both diffusion and downward infiltration of modem meteoric water and diffusive mixing from below with an old saline water body.     

Method for estimating available groundwater volume of small coral islands

We present a simple modelling method to estimate the volume of available groundwater in the freshwater lens of atoll islands under steady-state conditions. Model inputs include annual rainfall depth, island width for cross-sections along the length of the island, aquifer hydraulic conductivity, and depth to the contact between the upper sand aquifer and the lower limestone aquifer. The methodology is tested for nine islands of varying size in the Maldives and Micronesia. Sensitivity analysis indicates that lens volume on large islands typically is governed by the depth to the discontinuity, whereas lens volume for smaller islands is governed by rainfall rate and hydraulic conductivity. Volume curves, which relate lens volume to lens thickness, are developed for each of the nine islands and for three generic island shapes to allow rapid estimation of lens volume given field-estimated lens thickness. The methods presented in this study can be used for any small atoll island.     

Constraints on hydraulic parameters and implications for groundwater flux across the upland–estuary interface

Analyses of independent laboratory- and field-scale measurements from two sites on Sapelo Island, Georgia reveal heterogeneity in hydraulic parameters across the upland–estuary interface. Regardless of the method used (short-duration pumping tests, amplitude attenuation of tidal pumping data, sediment grain size distributions, and falling head permeameter tests), we obtain hydraulic conductivity of 10⁻⁴ m s⁻¹ for the fine-grained, well-sorted, clean sands that make up the upland areas. Proximal to the upland–estuary boundary, the tidal pumping analyses and permeameter tests suggest that hydraulic conductivities decrease by more than two orders of magnitude, a result consistent with the presence of a clogging layer. Such a clogging layer may arise due to a variety of physical, chemical, or biological processes. The extent and orientation of the layers of reduced hydraulic conductivity near the upland–estuary boundary influence the nature of the aquifer's response to tidal forcing. Where the lower conductivity layer forms a relatively flat creek bank, tidal pumping produces a primarily mechanical response in the adjacent aquifer. Where the creek bank is nearly vertical, there is a more direct hydraulic connection between the tidal creek and the adjacent aquifer. The clogging layer likely contributes to the development of complicated flow pathways across the upland–estuary boundary. Effective flow paths calculated from tidal pumping data terminate within the marsh, beyond the boundary of the upland aquifer, suggesting a diffuse regime of groundwater discharge in the marsh. We postulate that, in many settings, submarsh flow may be as important as seepage faces for groundwater discharge into the marsh–estuary complex.     

Groundwater recharge in the Victoria Nile basin of east Africa: support for the soil moisture balance approach using stable isotope tracers and flow modelling

Across equatorial Africa, increasing demand for groundwater has raised concerns about resource sustainability and has highlighted the need for reliable estimates of groundwater recharge. Recharge investigations in this environment are typically inhibited by a shortage of good quality meteorological and hydrogeological records. Moreover, when recharge studies are attempted they tend to rely on a single technique and frequently lack corroborating evidence to substantiate recharge predictions. In recent studies undertaken in the Aroca catchment of the Victoria Nile basin in central Uganda, the timing and magnitude of recharge determined by a soil moisture balance approach are supported by stable isotope data and groundwater flow modelling. The soil moisture balance study reveals that recharge averages in the order of 200 mm year⁻¹ and is more dependent on the number of heavy (more than 10 turn day⁻¹) rainfall events than the total annual volume of rainfall. Stable isotope data suggest independently that recharge occurs during the heaviest rains of the monsoons, and further establish that recharge stems entirely from the direct infiltration of rainfall, an assumption implicit in the soil moisture balance approach. Deforestation over the last 30 years is shown to have more than doubled the recharge estimate. Aquifer flow modelling supports the recharge estimates but demonstrates that the vast majority (over 99%) of recharging waters must be transmitted by the aquifer in the regolith rather the underlying bedrock fractures which have traditionally been developed for rural water supplies.     

Numerical Modeling of Atoll Island Hydrogeology

We implemented Ayers and Vachers' (1986) inclusive conceptual model for atoll island aquifers in a comprehensive numerical modeling study to evaluate the response of the fresh water lens to selected controlling climatic and geologic variables. Climatic factors include both constant and time-varying recharge rates, with particular attention paid to the effects of El Niño and the associated drought it brings to the western Pacific. Geologic factors include island width; hydraulic conductivity of the uppermost Holocene-age aquifer, which contains the fresh water lens; the depth to the contact with the underlying, and much more conductive, Pleistocene karst aquifer, which transmits tidal signals to the base of the lens; and the presence or absence of a semiconfining reef flat plate on the ocean side. Sensitivity analyses of steady-steady simulations show that lens thickness is most strongly sensitive to the depth to the Holocene-Pleistocene contact and to the hydraulic conductivity of the Holocene aquifer, respectively. Comparisons between modeling results and published observations of atoll island lens thicknesses suggest a hydraulic conductivity of approximately 50 m/d for leeward islands and approximately 400 m/d for windward islands. Results of transient simulations show that lens thickness fluctuations during average seasonal conditions and El Niño events are quite sensitive to island width, recharge rate, and hydraulic conductivity of the Holocene aquifer. In general, the depletion of the lens during drought conditions is most drastic for small, windward islands. Simulation results suggest that recovery from a 6-month drought requires about 1.5 years.     

Analysis of an unconfined aquifer subject to asynchronous dual-tide propagation

Most published solutions for aquifer responses to ocean tides focus on the one-sided attenuation of the signal as it propagates inland. However, island aquifers experience periodic forcing from the entire coast, which can lead to integrated effects of different tidal signals, especially on narrow high-permeability islands. In general, studies disregard a potential time lag as the tidal wave sweeps around the island. We present a one-dimensional analytical solution to the ground water flow equation subject to asynchronous and asymmetric oscillating head conditions on opposite boundaries and test it on data from an unconfined volcanic aquifer in Maui. The solution considers sediment-damping effects at the coastline. The response of Maui Aquifers indicate that water table elevations near the center of the aquifer are influenced by a combination of tides from opposite coasts. A better match between the observed ground water head and the theoretical response can be obtained with the proposed dual-tide solution than with single-sided solutions. Hydraulic diffusivity was estimated to be 2.3 × 10⁷ m²/d. This translates into a hydraulic conductivity of 500 m/d, assuming a specific yield of 0.04 and an aquifer thickness of 1.8 km. A numerical experiment confirmed the hydraulic diffusivity value and showed that the y -intercepts of the modal attenuation and phase differences estimated by regression can approximate damping factors caused by low-permeability units at the boundary.     

High concentration of solutes at the upper part of the saturated zone (water table) of a deep aquifer under sewage-irrigated land

The water-table region (upper 50 cm of the saturated zone) of a 25 m deep phreatic sandstone aquifer, lying under fields irrigated with sewage effluents for up to 22 yrs, was monitored in 1971 and 1984. Average concentrations of NO⁻₃, Cl⁻ and SO²⁻₄ of up to 225, 307 and 155 mg l⁻¹, respectively, were detected in the upper 50 cm of the saturated region in two research wells in 1984. These concentrations, which are related to effluent and fertilizer input to groundwater, were two to four times higher than those found deep (37–55 m) below the water table in nearby (1000 m distant) production wells. Nitrate data and the estimated transit time through the unsaturated zone (2 m yr⁻¹) support the model suggesting that the major source of nitrate pollution in the past should be related to the oxidation of soil organic matter. The SO²⁻₄/Cl⁻ ratio is found to be a useful indicator for the arrival of SO²⁻₄-fertilizers at the groundwater interface. The observations presented in this paper question the suitability of plans for using effluents as a water source for agriculture in regions which are the replenishment areas of phreatic aquifers.     

Tracing groundwater flow in the Borden aquifer using krypton-85

Krypton-85 was measured in air, soil gas, and ground water at the Borden aquifer in Ontario in October 1989. The measured specific activities in air and soil gas were 52.0 ± 2.0 and 53.6 ± 1.8 disintegrations per min (dpm) cm⁻³ krypton. These measurements are in excellent agreement with the global atmospheric trend and demonstrate that krypton-85 enters the water table at the Borden site without a lag in the soil gas reservoir. The krypton-85 specific activity in five groundwater samples ranged from 44.9 to 9.5 dpm cm⁻³ corresponding to groundwater ages of 2–17 years with a monotonic decrease in specific activity (increase in age) along the groundwater flow path. Travel times calculated from a two-dimensional steady-state model of groundwater flow agree well with the krypton-85 ages in the main recharge region of the aquifer where flow is predominantly vertical, but were 30–40% older than the krypton-85 age downstream of the main recharge area where the flow is mainly horizontal. The effect of dispersion on the distribution of krypton-85 was determined by modelling the transport of krypton-85 in the Borden aquifer with a two-dimensional time-dependent advection dispersion model using the steady-state flow field. Agreement between model specific activity and observed specific activity was excellent for samples in the main recharge region, but the model specific activities were 30–50% lower than observed specific activities in the region of horizontal flow. The differences in travel times and krypton-85 ages and in model krypton-85 and observed krypton-85 specific activities are considered to be small given the heterogeneities that exist in the hydraulic conductivity and aquifer geometry and hence in the groundwater flow field. The model simulated krypton-85 distribution was not sensitive to changes in longitudinal dispersivity and was only weakly sensitive to changes in transverse dispersivity. The geochemical inertness, well-defined source function, and insensitivity to dispersion of krypton-85 allow estimates of groundwater age to be made in a straightforward manner and measurement of krypton-85 can significantly enhance the characterization of groundwater flow in many shallow subsurface systems.     

Axisymmetric Modeling Using MODFLOW-USG

Axisymmetric groundwater models are used for simulating radially symmetric conditions. Groundwater simulators built specifically to model axisymmetric conditions are most commonly used for simulating aquifer tests. Although some numerical models capable of simulating flow and solute transport that are developed in the cartesian coordinate system framework offer flexibility to simulate axisymmetric conditions, most of the numerical groundwater models, such as the MODFLOW family of codes, are based on structured grids in which axisymmetric flows cannot be directly simulated. Researchers in the past have provided methods to manipulate aquifer properties to mimic axisymmetric conditions. This study presents a methodology that takes advantage of the unstructured grids of MODFLOW-USG to simulate axisymmetric models within the MODFLOW framework. To develop axisymmetric models, the intercell interface area arrays of MODFLOW-USG were calculated to accurately represent coaxial cylindrical model cells. Three examples are presented to demonstrate the application of MODFLOW-USG for axisymmetric modeling: a pumping well with delayed yield effects, a vadose zone flow model simulating an infiltration basin, and a density-dependent saltwater intrusion problem for a circular island. Results were verified against analytical solutions and published numerical codes.     

Combining measurements of transpiration and stable isotopes of water to determine groundwater discharge from forests

Discharge of saline groundwater from Eucalyptus forests on a semi-arid floodplain was directly determined by first measuring transpiration rates from the forests, and then partitioning the transpiration flux into groundwater discharge and soil water depletion. This partitioning was achieved by identifying the source of the transpired water with naturally occurring stable isotopes of water. Transpiration rates were low, being about 0.3 mm day⁻¹ from three E. largiflorens sites and up to 2 mm day⁻¹ from an E. camaldulensis site. There was no significant variation in transpiration across seasons, indicating that transpiration was limited by environmental factors other than evaporative demand. Despite its salinity (electrical conductivities of 11–33 dS m⁻¹), the groundwater was used by the forests at all sites and all times, and made up 100% of transpiration in more than half of the measurements, and 40–80% in the remainder. There was some consistency in water uptake patterns. E. camaldulensis tended to take up shallow soil water and groundwater simultaneously, as did trees at one of the E. largiflorens sites. At the driest sampling time, however, groundwater was the only source of water for trees at both of these sites. Trees at the remaining two E. largiflorens sites generally relied solely on the groundwater. The tree water source results indicate that groundwater discharge fluxes were between 40 and 100% of the transpiration fluxes at these sites. These groundwater discharge fluxes were small in terms of regional groundwater balances, but would be important in the salinisation of the soils. Additionally, uptake of water from the soil profile by the trees substantially increased groundwater discharge compared with discharge from the soils had they been bare of vegetation.     

Influence of variable salinity conditions in a tidal creek on riparian groundwater flow and salinity dynamics

While recent studies have revealed that tidal fluctuations in an estuary significantly affect groundwater flows and salt transport in the riparian zone, only seawater salinity in the estuary has been considered. A numerical study is conducted to investigate the influence of estuarine salinity variations on the groundwater flow and salt dynamics in the adjacent aquifer to extend our understanding of these complex and dynamic systems. Tidal salinity fluctuations (synchronous with estuary stage) were found to alter the magnitude and distribution of groundwater discharge to the estuary, which subsequently impacted on groundwater salinity patterns and residence times, especially in the riparian zone. The effects of salinity fluctuations were not fully captured by adopting a constant, time-averaged estuarine salinity. The modelling analysis also included an assessment of the impact of a seasonal freshwater flush in the estuary, similar to that expected in tropical climates (e.g. mean estuary level during flood significantly greater than average), on adjacent groundwater flow and salinity conditions. The three-month freshwater flushing event temporarily disrupted the salt distribution and re-circulation patterns predicted to occur under conditions of constant salinity and tidal water level fluctuations in the estuary. The results indicate that the salinity variations in tidal estuaries impact significantly on estuary–aquifer interaction and need to be accounted for to properly assess salinity and flow dynamics and groundwater residence times of riparian zones.     

Water flux estimates from a Belgian Scots pine stand: a comparison of different approaches

Four distinct approaches, that vary markedly in the spatial and temporal resolution of their measurement and process-level outputs, are used to investigate the daily and seasonal water vapour exchange in a 70-year-old Belgian Scots pine forest. Transpiration, canopy interception, soil evaporation and evapotranspiration are simulated, using a stand-level process model (SECRETS) and a soil water balance model (WAVE). Simulated transpiration was compared with up-scaled sap flow measurements and simulated evapotranspiration to eddy covariance measurements.

Reasonable agreement in the temporal trends and in the annual water balance between the two models was observed, however daily and weekly predictions often diverged. Most notably, WAVE estimated very low, to no transpiration during late autumn, winter and early spring when incident radiation fell below 50 W m⁻² while SECRETS simulated low (0.1–0.4 mm day⁻¹) fluxes during the same period. Both models exhibited similar daily trends in simulated transpiration when compared with sap flow estimates, although simulations from SECRETS were more closely aligned. In contrast, WAVE over-estimated transpiration during periods of no rainfall and under-estimated transpiration during rainfall. Yearly, total evapotranspiration simulated by the models were similar, i.e. 658 mm (1997) and 632 mm (1998) for WAVE and 567 mm (1997) and 619 mm (1998) for SECRETS.

Maximum weekly-average evapotranspiration for WAVE exceeded 5 mm day⁻¹, while SECRETS never exceeded 4 mm day⁻¹. Both models, in general, simulated higher evapotranspiration than that measured with the eddy covariance technique. An impact of the soil water content in the direct relationship between the models and the eddy covariance measurements was found.

The results suggest that: (1) different model formulations can reproduce similar results depending on the scale at which outputs are resolved, (2) SECRETS estimates of transpiration were well correlated with the empirical measurements, and (3) neither model fitted favourably to the eddy covariance technique.     

Isotope hydrological study of mean transit times in an alpine basin (Wimbachtal, Germany)

Measurements of tritium and ¹⁸O concentrations in precipitation and runoff were used to provide further insight into the groundwater storage properties of the Wimbachtal Valley, a catchment area of 33.4 km², extending between 636 and 2713 m a.s.l. in the Berchtesgaden Alps. The catchment includes three aquifer types: a dominant porous aquifer; a fractured dolomite; a karstic limestone aquifer. Employing a simple hydrological model, information about mean transit times of environmental tracers is derived for the groundwater runoff component and several karst springs from the application of the exponential and dispersion flow models to the isotopic input and output data. The mean transit times calculated from a dispersion model with transit times of 4.1 years for ¹⁸O and 4.2 years for tritium, which agree well, allow calculation of total (mobile + stagnant) groundwater storage volume, which is equivalent to 6.6 m of water depth. Direct runoff appears negligible as in many other cases.     

Groundwater dynamics of Fongafale Islet, Funafuti Atoll, Tuvalu

Geoelectric and hydrologic surveys during spring tides revealed the spatiotemporal distribution of groundwater quality produced by tidal forcing in Fongafale Islet, Funafuti Atoll, Tuvalu. The observed low resistivity showed that saline water largely immersed the surficial Holocene aquifer, indicating that there is no thick freshwater lens in Fongafale Islet, unlike in other atoll islands of comparable size. Half of the islet was constructed by reclaiming the original swamp with porous, highly permeable coral blocks; this reclaimed area should not be considered as part of the islet width for calculation of the expected thickness of the freshwater lens. The degree of aquifer salinization depends on the topographic characteristics and the hydrologic controls on the inland propagation of the tidal forcing. Large changes in bulk resistivity and the electrical conductivity of groundwater from wells indicate that periodic salinization in phase with the semidiurnal tides was occurring widely, especially in areas at lower elevation than the high-tide level and in reclaimed areas with high permeability. Thin sheets of nearly fresh and brackish water were observed in the surficial aquifer in areas above the high-tide level and in taro swamps, respectively. The thinness of the brackish and freshwater sheets suggests that the taro swamps and the fresh groundwater resources of the islet are highly vulnerable to salinization from anticipated sea-level rise. An understanding of the inherent geologic and topographic features of an atoll is necessary to evaluate the groundwater resources of the atoll and assess the vulnerability of its water resources to climate change.     

Groundwater level monitoring by precision measurements

The paper presents the results of processing the data on variations of groundwater level in an aquifer opened by a well in a depth interval of 76–115 m. The obtained series are compared with the volume deformation of the aquifer, calculated theoretically using the ETERNA 3.0 tidal analysis program. Two components are detected in the long-period area: the one coherent to the luni-solar tides and the one incoherent to them. It is established that the phase shift between the variations of volume deformation of the aquifer and the corresponding variations of water level in the well depends not only on the stressed-deformed state of the surrounding rock massif, but on the position of the piezometric surface of the aquifer as well. Some peculiarities of the effects of luni-solar tides on variations of groundwater level are considered. The amplitudes of groundwater level and volume deformation of the aquifer, as well as the phase shift between these parameters are estimated for the detected tidal waves. The value of the amplitude factor lies in the range of 0.04–0.058 cm/nanostrain.     

Dispersion of naturally occurring ions in groundwater from various rock types in a portion of the San Pedro river basin, Arizona

The groundwater in an alluvial basin in southern Arizona was analyzed for concentrations of Ca⁺⁺, Mg⁺⁺, Na⁺, and Cl⁻, ions.

The variety of rock types in the area, plus the undisturbed state of the groundwater basin, make comparative rock mineralization-groundwater ionization interpretations possible. Ionic dispersion in groundwater eminating from source areas composed of differing rock types is plotted as isogram maps. These isolated areas of differing mineral composition each exhibit a unique ionic contribution to groundwater. The ion concentrations in groundwater were then used as naturally occurring tracers to determine source areas of recharge and to delineate subsurface barriers to the normal basin flow net. Ion dispersion plots reveal the carbonates of the Dragoon Mountains to be a major contributor of Ca⁺⁺ and Mg⁺⁺ to the deep alluvial portion of the basin. Cl⁻ dispersion patterns show the granitic intrusives of the Tombstone Hills produce a barrier effect in the normal flow pattern of the basin as well as being a contributor of Cl⁻ to groundwater.     

Formation of a salt plume in the Coastal Plain aquifer of Israel: the Be''er Toviyya region

The formation and development of a salt plume (salinity up to 800 mg Cl 1⁻¹) in the inner part of the Coastal Plain aquifer of Israel is analyzed. Massive groundwater exploitation during the 1950s caused a large drop in the water level and formation of a hydrologic depression in the Be'er Toviyya-Kefar Warburg area. The depression reached a maximal depth during the late 1960s; thereafter a reduction in the rate of pumpage led to restoration of water levels and shallowing of the depression, until its complete disappearance towards the end of the 1980s. A spot of high salinity first appeared in 1956, following a deep drawdown in the water levels. This saline plume has been continuously expanding with increasing salinity concentrations (200–800 mg Cl 1⁻¹) in its center. The average rate of radial expansion was about 50 m year⁻¹. The expansion and salinization did not cease as the depression disappeared. Rather, equalization of water levels in wells situated within the plume area with those of situated along its margins resulted in the salinization of the latter within a period of 1 year.

Mass balances for water and chloride contents were made for the period 1967–1990. Taking into consideration the storage change, pumpage, natural replenishment and artificial recharge, the lateral inflow to the depression is estimated as 60 × 10⁶ m³. Upon addition of the chloride balance, and taking into consideration the chloride concentrations of the surrounding fresh water and the apparent possible end-member of the saline source (based on geochemical considerations), the saline inflow is estimated as (40–60) × 10⁶ m³. These estimates indicate that a large amount of saline water penetrated into the aquifer, of about half of the natural replenishment of the study area, with an estimated salinity of 1900–2700 mg Cl 1⁻¹.

It is suggested that the salt plume was formed as a result of a drop in water level combined with a flow of underlying saline water bodies from deeper strata. The chemical composition of the groundwater points to the existence of two saline water bodies of Ca-chloride composition and a marine Br/Cl ratio: (1) saline water with low Na/Cl (0.6), So₄/Cl, and B/Cl ratio; (2) saline water with higher Na/Cl (> 0.6), So₄/Cl, and B/Cl ratios. These chemical compositions resemble Ca-chloride saline waters found in other locations in the Coastal Plain aquifer and in underlying formations. The saline water bodies may occur in either pockets at the bottom of the aquifer or lumachelle and sandstone layers of high hydraulic conductivity in underlying sediments.     

Seasonal soil water storage changes beneath central Amazonian rainforest and pasture

Evaporation and infiltration were compared for tropical rainforest and pasture, near to Manaus, Brazil from October 1990 to February 1992 using measurements of soil water storage over a depth of 2 m. The soil is a clayey oxisol of low water available capacity. In both of the dry seasons studied, the maximum change in soil water storage in the forest was 154 mm and in the pasture it was 131 and 112 mm. Similar behaviour of the soil water reservoir below forest and pasture in the wet season implied that differences in evaporation and drainage were small. In the dry season, soil water storage behaviour in the upper metre of the soil was similar but there were marked differences in the second metre. The pasture took up little water from below 1.5 m but the forest appeared to utilise all of the available water in the 2 m profile in both seasons.

The water balance of the 2 m profile showed that the pasture evaporation rate was equal to that of the forest until storage had decreased 80 mm from the maximum. There was then a decline in pasture evaporation rate to 1.2 mm day⁻¹ as the storage decreased by a further 50 mm. In contrast, the forest uptake rate remained above 3.5 mm day⁻¹ until storage had decreased 140 mm from the maximum (within 15 mm of the extraction limit), before declining abruptly to less than 1.5 mm day⁻¹. There was strong evidence that the forest was able to abstract water from depths greater than 3.6 m.

Spatial variability of soil water storage was significantly greater beneath the pasture than beneath the forest, particularly following rainfall events in the dry season. This was largely the result of redistribution of rainfall as local surface runoff. There was no evidence of redistribution or runoff in the forest.