Geometry and dynamics of the freshwater—seawater interface in a coastal aquifer in southeastern Spain

Abstract

The contact between freshwater and seawater in coastal aquifers is studied using a relatively simple model for homogeneous aquifers. However, for real aquifers it is not so simple. The desalination plant built to supply water to the city of Almería is situated over the aquifer in the southern part of the River Andarax Delta. Its design capacity is 1100 L s^?1, and it is supplied from boreholes pumping water from beneath the freshwater—seawater contact in this aquifer. Well logs kept over a period of two years have allowed us to accurately define the interface geometry of the freshwater—seawater contact. Lithological data collected from 31 boreholes have also indicated the existence of strata with low hydraulic conductivity, within others of high conductivity. During a simultaneous pumping test of six wells with 690 L s^?1 total discharge, electrical conductivity measurements showed the influx of seawater 6–10 m below sea level and a drawdown of the interface in the piezometers closest to the pumping wells.     

Estimation of anisotropic hydraulic conductivity using geophysical data in a coastal aquifer of Karnataka,India

Estimation of hydraulic parameters is essential to understand the interaction between groundwater flow and seawater intrusion. Though several studies have addressed hydraulic parameter estimation, based on pumping tests as well as geophysical methods, not many studies have addressed the problem with clayey formations being present. In this study, a methodology is proposed to estimate anisotropic hydraulic conductivity and porosity values for the coastal aquifer with unconsolidated formations. For this purpose, the one-dimensional resistivity of the aquifer and the groundwater conductivity data are used to estimate porosity at discrete points. The hydraulic conductivity values are estimated by its mutual dependence with porosity and petrophysical parameters. From these estimated values, the bilinear relationship between hydraulic conductivity and aquifer resistivity is established based on the clay content of the sampled formation. The methodology is applied on a coastal aquifer along with the coastal Karnataka, India, which has significant clayey formations embedded in unconsolidated rock. The estimation of hydraulic conductivity values from the established correlations has a correlation coefficient of 0.83 with pumping test data, indicating good reliability of the methodology. The established correlations also enable the estimation of horizontal hydraulic conductivity on two-dimensional resistivity sections, which was not addressed by earlier studies. The inventive approach of using the established bilinear correlations at one-dimensional to two-dimensional resistivity sections is verified by the comparison method. The horizontal hydraulic conductivity agrees with previous findings from inverse modelling. Additionally, this study provides critical insights into the estimation of vertical hydraulic conductivity and an equation is formulated which relates vertical hydraulic conductivity with horizontal. Based on the approach presented, the anisotropic hydraulic conductivity of any type aquifer with embedded clayey formations can be estimated. The anisotropic hydraulic conductivity has the potential to be used as an important input to the groundwater models.     

Influence of a coarse interlayer on seawater intrusion and contaminant migration in coastal aquifers

Vertical 2D slice laboratory experiments were carried out in homogenous and layered sand tanks to elucidate the effects of a highly permeable (coarse‐grained sand) interlayer on seawater intrusion and transport of contaminants to a coastal sea. Tidal fluctuations produced oscillations in the seawater–freshwater transition zone, fluctuations of the contaminant infiltration rate and a zigzag contaminant plume outline. The seawater wedge became discontinuous at the (vertical) edges of the interlayer because of increased lateral movement of the seawater–freshwater interface within the interlayer. The contaminant plume formed a tail within the interlayer depending on the tidal stage, and similar to the wedge, its movement was accentuated. A simple analytical model that neglected vertical flow reliably predicted steady‐state seawater intrusion into the coastal aquifer. Numerical modeling was used to gain insight into the groundwater hydrodynamics and contaminant migration. The numerical results confirmed the experimental findings, i.e. that a highly permeable interlayer can provide a rapid transit path for contaminants to reach the seaward boundary and that the interlayer amplifies the effects of tidal fluctuations, resulting in wider transition zones for the seawater wedge and contaminant plume. Numerical simulations further showed that, with increasing interlayer hydraulic conductivity, the maximum seawater intrusion distance inside the interlayer increases approximately linearly. For the fixed‐head contaminant injection condition used, the model showed that contaminant infiltration increases approximately logarithmically with increasing interlayer hydraulic conductivity (other factors held fixed). Copyright © 2013 John Wiley & Sons, Ltd.     

Adjustment of a Formulated Sediment for Sediment Testing with Caenorhabditis elegans (Nematoda)

The OECD soil, which is standardised for use in the sediment toxicity with the midge Chironomus riparius revealed inhibition of the nematode Caenorhabditis elegans in body growth length and egg production of 41.0% and 84%, respectively. 18 formulated sediments were tested to optimise the sediment composition for growth and reproduction of C. elegans. Their components were chosen to simulate native sediments, whereby different quantities and different clay minerals were tested. The mixture that we found to be optimal in our experiments consisted of 5% sphagnum peat, 70% calcitic sand, 0.5% dolomite limestone, 4.5% iron(III) oxide, and a clay combination of 1.5% chlorite and 18.5% aluminium(III) oxide. When applying this mixture to C. riparius, the new formulated sediment improved also growth in comparison to the OECD sediment.     

Modelling of colloidal particle and heavy metal transfer behaviours during seawater intrusion and refreshing processes

The proper management of coastal aquifers commonly requires an understanding of regional mass flow and complete seawater–freshwater circulation. In this study, time series observations of seawater intrusion and refreshing were conducted using a column experiment based on natural flow conditions in coastal groundwater and a sampled medium from a coastal sandy aquifer without chemical treatment. Ranges of hydrodynamic and hydrochemical variables were tested and analysed. The results showed that the zeta potential of suspended colloids in aqueous solution in an aquifer polluted with 0.5 g/kg of heavy metals exhibited an isoelectric point for pH values ranging from 5.70 to 6.07 when freshwater or seawater completely occupied the aquifer pores, which is representative of natural hydrochemical conditions. In this scenario, a high background concentration of heavy metals induced colloidal immobilization. Otherwise, seawater–freshwater circulation enabled colloid mobilization due to ionic strength and pH fluctuations. The migration of multiple heavy metals occurred at a characteristic time of approximately 1 pore volume after each intrusion stage began and when the peak rate of colloid release was reached. At these times, the colloid behaviour determined the quantity and pathway of heavy metal transport. On the basis of the influences of seawater and freshwater interactions, the quantity of mobilized particles generally decreased and was uniformly distributed in each fraction due to particle loss and decreased porous connectivity. We speculate that the decrease in the total surface area of the migratory colloids may cause colloid‐associated heavy metal transport to decrease. The experimental results provide a useful basis for testing coastal groundwater flow and mass transport models because these phenomena require full characterization to precisely evaluate the associated fluxes from the field scale to the microscopic dimension.     

Using a pore‐scale model to quantify the effect of particle re‐arrangement on pore structure and hydraulic properties

A pore‐scale model based on measured particle size distributions has been used to quantify the changes in pore space geometry of packed soil columns resulting from a dilution in electrolyte concentration from 500 to 1 mmol l^?1 NaCl during leaching. This was applied to examine the effects of particle release and re‐deposition on pore structure and hydraulic properties. Two different soils, an agricultural soil and a mining residue, were investigated with respect to the change in hydraulic properties. The mining residue was much more affected by this process with the water saturated hydraulic conductivity decreasing to 0·4% of the initial value and the air‐entry value changing from 20 to 50 cm. For agricultural soil, there was little detectable shift in the water retention curve but the saturated hydraulic conductivity decreased to 8·5% of the initial value. This was attributed to localized pore clogging (similar to a surface seal) affecting hydraulic conductivity, but not the microscopically measured pore‐size distribution or water retention. We modelled the soil structure at the pore scale to explain the different responses of the two soils to the experimental conditions. The size of the pores was determined as a function of deposited clay particles. The modal pore size of the agricultural soil as indicated by the constant water retention curve was 45 µm and was not affected by the leaching process. In the case of the mining residue, the mode changed from 75 to 45 µm. This reduction of pore size corresponds to an increase of capillary forces that is related to the measured shift of the water retention curve. Copyright © 2007 John Wiley & Sons, Ltd.     

Source speciation resolving hydrochemical complexity of coastal aquifers

There is a growing concern of seawater intrusion to freshwater aquifers due to groundwater overexploitation in the eastern coastal belt of Southern India. The problem becomes complex in the regions where industrial effluents are also contaminating the freshwater aquifers. In order to understand the hydrochemical complexity of the system, topographic elevation, static water level measurements, major ion chemistry, ionic cross plots, water type contours and factor analysis were applied for 144 groundwater samples of shallow and deep sources from Quaternary and Tertiary coastal aquifers, located within the industrial zone of 25 km² area near Cuddalore, Southern India. The ionic cross plots indicates dissolution of halite minerals from marine sources and seawater mixing into inland aquifers up to the level of 9.3%. The factor analysis explains three significant factors totaling 86.3% of cumulative sample variance which includes varying contribution from marine, industrial effluent and freshwater sources.     

Evolution of Heterogeneous Mixed Siliciclastic/Carbonate Aquifers Containing Metastable Sediments

Mixed carbonate and siliciclastic marine sediments commonly become freshwater aquifers in eastern coastal regions of the United States and many other global locations. As these deposits age, the carbonate fraction of the sediment is commonly removed by dissolution and the aquifer can become a solely siliciclastic system or contain zones or beds of pure quartz sand. During aquifer evolution, the sediment grain size characteristics, hydraulic conductivity, and porosity change. An investigation of these changes using mixed carbonate/siliciclastic sediment samples collected from a modern barrier island beach in southern Florida showed that the average mean grain diameter decreased with removal of the carbonate fraction, but the average hydraulic conductivity and porosity increased slightly, but not to statistical significance. This counterintuitive result occurs because of the change in the pore types from a combined shelter and intergranular pore system producing a dual porosity system in the mixed sediments to a single intergranular pore system in the siliciclastic sediment fraction. Within the mixed carbonate/siliciclastic sediment, in the pure carbonate fraction, large shell fractions form grain‐supported large pores, which become filled with sand‐sized quartz as the shell fragments decrease in size or as the sediment becomes compacted. The hydraulic conductivity increases because the shell fragments that were oriented perpendicular to flow caused an increase in the length of the flow path, or a larger scale tortuosity, compared with the flow through pure quartz sand.     

Vertical Hydraulic Gradient Estimation in Clay Till,Using MiHPT Advanced Direct-Push Technology

The vertical transport of contaminants from source areas is employed in many risk assessment models and screening tools in order to estimate the contaminant mass discharge (CMD) into underlying aquifers. The key parameters for estimating CMD are the contaminant source area and concentration, and the vertical water flux, the latter of which depends on the vertical hydraulic conductivity and the vertical hydraulic gradient in the subsurface. This study focuses on advancing the use of the combined membrane interface probe hydraulic profiling tool (MiHPT) to investigate the vertical hydraulic gradient across a clay till overlying a sandy aquifer at a contaminated site in Denmark. Only the HPT is necessary for the estimate of vertical hydraulic gradient. The hydraulic head, clay till thickness, and resulting vertical hydraulic gradients found using the MiHPT compared well with observations from nearby nested wells. The parameter with the largest discrepancy was the thickness of the clay till. The advantage of the MiHPT is its relatively quick depth discrete access to information regarding subsurface permeability, vertical hydraulic gradients and contaminant distribution across a site. In this case study, performance of additional dissipationtests during the HPT log to acquire determination of the vertical hydraulic gradient increased the cost by 3% compared to standard HPT logs.     

Quantifying transient post‐overwash aquifer recovery for atoll islands in the Western Pacific

Marine overwash events for atoll islands in the Pacific and Indian Oceans, which cause salinization of fresh groundwater because of infiltrating seawater, pose a significant challenge for island community sustainability in regard to water supply. Understanding transient fresh groundwater development during a post‐overwash period for a range of island sizes, geologic characteristics, and rainfall patterns is essential for water management. This paper presents a methodology for quantifying this development for an atoll nation, with methods applied to the 32 atolls of the Federated States of Micronesia (FSM) in the western Pacific. Using the numerical groundwater modelling code SUTRA, overwash events and post‐overwash freshwater–seawater dynamics are simulated for the range of island widths (200 to 1100 m), geologic characteristics (hydraulic conductivity corresponding to leeward and windward islands), and rainfall patterns (western, central, and eastern regions) present in the FSM, thereby providing results for each atoll island. Results show that 10–17, 8–12, and 6–12 months are required to achieve 60% freshwater lens recovery for leeward islands in the western, central, and eastern FSM, respectively, with variation due to rainfall rate and island width. In contrast, 4–9 months is required for 60% recovery for windward islands. However, the natural thinness of the lend on windward islands typically precludes extensive use of groundwater under average rainfall conditions. Overwash characteristics (depth, duration, and seasonal timing) did not significantly affect recovery times. For the region of lowest rainfall (western FSM), 6–10 months is required to achieve potable groundwater at the typical depth of hand‐dug wells. Results provide water resource managers and atoll island communities with important information regarding timing of potential fresh groundwater use following an overwash event. Copyright © 2015 John Wiley & Sons, Ltd.     

Groundwater,salt and contaminant transport in a coastal aquifer system with a low-permeability layer in the intertidal zone

Low-permeability layer (LPL), formed by natural deposit or artificial reclamation and commonly found below the intertidal zone of coastal groundwater system, can retard the ingress of seawater and contaminants, and shorten the travel time of the land-sourced contaminant to the marine environment compared with a homogenous sandy coastal aquifer. However, there is limited understanding on how an intertidal LPL, a condition occurred in a coastal aquifer at Moreton Bay, Australia, influences the groundwater and contaminant transport across the shallow beach aquifer system. We characterized the aquifer hydrological parameters, monitored the in situ groundwater heads, and constructed a 2-D numerical model to analyses the cross-shore hydrological processes in this stratified system. The calibrated model suggests that in the lower aquifer, the inland-source fresh groundwater flowed horizontally towards the sea, upwelled along the freshwater–saltwater interface, and exited the aquifer at the shore below the LPL. Whereas in the upper aquifer, the tidally driven seawater circulation formed a barrier that prevented fresh groundwater from horizontal transport and discharge to the beach above the LPL, thereby directing its leakage to the lower aquifer. A contaminant represented by a conservative tracer was ‘released’ the upper aquifer in the model and results showed that the spreading extent of the contaminant plume, the maximum rate of contaminant discharge to the ocean, and its plume length decreased compared with a simulation case in a homogenous sandy aquifer. Sensitivity analysis was also conducted to investigate the characteristics of the LPL, including its continuity and hydraulic conductivity, which were found to vary along the beach at Moreton Bay. The result shows that with a lower hydraulic conductivity and continuous layer of LPL reduced the groundwater exchange and contaminant transport between upper and lower aquifer. The findings from the combined field and modelling investigations on the impact of an intertidal LPL on coastal aquifer systems highlight its significant implications to alter the groundwater and mass transport across the land–ocean interface.     

Groundwater modelling and marine intrusion in the semi‐arid systems (Cap‐Bon,Tunisia)

This article studies the effect of drought and pumping discharge on groundwater supplies and marine intrusion in the Korba aquifer (Cap‐Bon peninsula, Tunisia). The Groundwater Modelling System has been used to model the groundwater flow and to simulate the seawater intrusion. The calibration is based on the groundwater levels in the steady state from 1963, and in the transient state from the groundwater levels from 1963 to 2005. The main objective is to quantify the components of the groundwater mass balance and to estimate the hydraulic conductivity distribution. The impact of pumping discharge on the groundwater level evolution has been examined by two pumping scenarios P1 (no. 8420) and P2 (no. 8862) wells. The hydrodynamic modelling shows the increasing drawdowns after 14 years of pumping: 4 m in P1 well and about 5 m in P2 well below sea level. The drawdowns are accompanied by the inverse hydraulic gradient. The numerical model was used to discuss the management of the groundwater resources of Cap‐Bon. As the population continues to grow and the demand for groundwater pumping intensifies beyond the 1963 level, it can be expected that the actual extent of seawater intrusion in the future would be more severe than the model prediction. Better strategies for groundwater development and management will be necessary to protect the freshwater aquifers to the marine intrusion. Copyright © 2011 John Wiley & Sons, Ltd.     

Hydraulic conductivity in upland blanket peat: measurement and variability

A key parameter used in wetland hydrological and landform development models is hydraulic conductivity. Head recovery tests are often used to measure hydraulic conductivity, but the calculation techniques are usually confined to rigid soil theory. This is despite reports demonstrating the misapplication of rigid soil theory to non‐rigid soils such as peats. Although values of hydraulic conductivity calculated using compressible techniques have been presented for fenland peats, these data have never, to the authors' knowledge, been compared with such calculations in other peat types. Head recovery tests (slug withdrawal) were performed on piezometers at depths ranging from 10 to 80 cm from the surface on north Pennines blanket peats. Results were obtained using both rigid and compressible soil theories, thus allowing comparison of the two techniques. Compressible soil theory gives values for hydraulic conductivity that are typically a factor of five times less than rigid soil calculations. Hydraulic conductivity is often assumed to decrease with depth in upland peats, but at the study site in the northern Pennines it was not found to vary significantly with depth within the range of peat depths sampled. The variance within depth categories was not significantly different to the variance between depth categories showing that individual peat layers did not have characteristic hydraulic conductivity values. Thus, large lateral and vertical differences in hydraulic conductivity over short distances create problems for modelling but may help account for the high frequency of preferential flow pathways within what is otherwise a low matrix hydraulic conductivity peat. Hydraulic conductivity was found to vary significantly between sampling sites, demonstrating that hillslope‐ or catchment‐scale variability may be more important than plot‐scale variability. Values for compressibility of the peats are also reported. These generally decline with depth, and they also vary significantly between sampling sites. There are implications for the way in which measurements of hydraulic conductivity and other properties of blanket peat are interpreted, as the effects of environmental change in one part of a peat catchment may be very different to those in another. Copyright © 2002 John Wiley & Sons, Ltd.     

Stochastic analysis of flux and head moments in a heterogeneous aquifer system

This study investigates the behavior of flux and head in a strongly heterogeneous three-dimensional aquifer system. The analyses relied on data from 520 slug tests together with 38,000 one-foot core intervals lithological data from the site of the General Separations Area in central Savannah River Site, South Carolina, USA. The skewness in the hydraulic conductivity histograms supported the geologic information for the top two aquifers, but revealed stronger clay content, than was reported for the bottom aquifer. The log-normal distribution model described adequately the hydraulic conductivity measurements for all three aquifers although, other distributions described equally well the bottom aquifer measurements. No apparent anisotropy on the horizontal plane was found for the three aquifers, but ratios of horizontal to vertical correlation lengths between 33 and 75 indicated a strong stratification at the site. Three-dimensional Monte Carlo stochastic simulations utilized a grid with larger elements than the support volume of measurements, but of sub-REV (representative elementary volume) dimensions. This necessitated, on one hand, the use of upscaled hydraulic conductivity expressions, but on the other hand did not allow for the use of anisotropic effective hydraulic conductivity expressions (Sarris and Paleologos in J Stoch Environ Res Risk Assess 18: 188–197, 2004). Flux mean and standard deviations components were evaluated on three vertical cross-sections. The mean and variance of the horizontal flux component normal to a no-flow boundary tended to zero at approximately two to three integral scales from that boundary. Close to a prescribed head boundary both the mean and variance of the horizontal flux component normal to the boundary increased from a stable value attained at a distance of about five integral scales from that boundary. The velocity field 〈q_x〉 was found to be mildly anisotropic in the top two aquifers, becoming highly anisotropic in the bottom aquifer; 〈q_y〉 was anisotropic in all three aquifers with directions of high continuity normal to those of the 〈q_x〉 field; finally, 〈q_z〉 was highly anisotropic in all three aquifers, with higher continuity along the east–west direction. The mean head field was found to be continuous, despite the high heterogeneity of the underlying hydraulic conductivity field. Directions of high continuity were in alignment with field boundaries and mean flow direction. Conditioning did not influence significantly the expected value of the flux terms, with more pronounced being the effect on the standard deviation of the flux vector components. Conditioning reduced the standard deviations of the horizontal flux components by as much as 50% in the bottom aquifer. Variability in the head cross-sections was affected only marginally, with an average 10% reduction in the respective standard deviation. Finally, the location of the conditioning data did not appear to have a significant effect on the surrounding area, with uniform reduction in standard deviations.     

Comparison of Delivery Methods for Enhanced In Situ Remediation in Clay Till

Three methods for enhanced delivery of in situ remediation amendments in low-permeability deposits have been tested at a site in Denmark: pneumatic fracturing, direct-push delivery, and hydraulic fracturing. The testing was carried out at an uncontaminated part of a farm site, previously used for storage of chlorinated solvents, underlain by basal clay till with hydraulic conductivity ranging from 7.1× 10^–11 to 3.5 × 10^–7 m/s at testing depths 2.5 to 9.5 m b.s. Fluorescent tracers fluorescein and rhodamine WT were delivered. Tests of all three delivery methods have not been carried out at a single site before, and thus, this study provides unique data for comparison of enhanced delivery methods in both the vadose and saturated zone. Results show that pneumatic fracturing with nitrogen gas and propagation pressures of 1 to 9 bar had a distribution radius of less than 2 m, and produced dense networks of tracer-filled natural fractures above the redox boundary (0 to 3 m b.s.) and widely spaced, discrete, induced, tracer-filled subhorizontal fractures at depth (>3 m b.s.). Direct-push delivery at pressures of 8 to 30 bar had a distribution radius of approximately 1 m, distributed tracer primarily in natural fractures above the redox boundary and in discrete, closely spaced (but not merging) induced fractures below the redox boundary. Hydraulic fracturing with a sand-guar mixture at pressures of 0 to 6 bar produced an elliptical, asymmetrical, bowl-shaped fracture with a physical radius of approximately 3.5 m at 3 m b.s. The geometry of hydraulic fractures attempted emplaced at 6.5 and 9.5 m b.s. is uncertain, but clearly not horizontal as desired. The direct-push delivery method is robust and efficient for enhanced delivery at the clay till site in question, which based on thorough geological characterization is deemed a geologically representative basal clay till site.     

EXPLORATION OF BURUNDI NICKELIFEROUS LATERITES BY ELECTRICAL METHODS*

The laterites in Burundi, which are formed by weathering of ultrabasic rocks, show a complete profile with the following horizons: canga, the ferruginous crust capping, ferralite, consisting essentially of iron hydroxides, and saprolite, which contains a large quantity of hydrosilicate minerals. Nickel bearing minerals occur in the saprolite and the lower portion of ferralite. Resistivity well-logging and resistivity sounding indicated that the electrical properties of rocks depend upon their composition: Canga and ferralite showed high resistivities of 6,500 Ωm and 800 Ωm, respectively. The resistivity of saprolite was found to be much lower, between 10 Ωm and 20 Ωm. The laterite is underlain by resistive peridotite. The chargeability of saprolite was found to be lower than that of the upper horizons and the bedrock. Electrolytic conductivity of laterite, which depends on the geometry of the deposit, was found to be low, because the laterite contains moisture and ground water, which are highly resistive. The relatively high conductivity of saprolite is caused by nickeliferous hydrosilicates, which exhibit the electrical properties of clay minerals, with an apparent maximum conductivity of 0.25 S/m. The conductivity of saprolite corresponds to a concentration between 30% and 50% of conductive silicate minerals distributed in the pore space of deposit. A nickel enrichment of up to 6% was estimated from the resistivity of the saprolite. Prospecting for laterites by electrical sounding showed that the development of laterite horizons in a nickel deposit correlates with the surface morphology of weathered ultrabasic massif. Thus the method can be used in preliminary exploration of such deposits.     

Investigation of Failure Mechanisms and Migration of Organic Chemicals at Wilsonville, Illinois

Ground water contamination was discovered in 1981 in a monitoring well at the Earthline disposal facility near Wilsonville, Illinois. Organic chemicals had migrated at a rate 100 to 1000 times greater than predicted when the site received its permit to operate in 1978. Postulated failure mechanisms included migration through previously unmapped permeable zones, subsidence of an underground mine, organic-chemical and clay-mineral interactions, acid-mine drainage and clay interactions, trench-cover settlement, and erosion.
In this investigation, the Illinois State Geological Survey found the primary reason for the rapid migration: the presence of previously undetermined fractures and joints in glacial till. The inaccurate predictions of hydraulic conductivity were based on laboratory-determined values that did not adequately measure the effects of fractures and joints on the transit time calculations. Field-measured hydraulic conductivity values were generally 10 to 1000 times greater than their laboratory-measured counterparts, thus largely accounting for the discrepancy between predicted and actual migration rates in the transit time calculations. The problem was compounded, however, by the burial of liquid wastes and by trench covers that allowed excess surface runoff to enter the trenches. Organic-chemical and clay-mineral interactions may also have exacerbated the problem in areas where liquid organic wastes were buried.     

Modelling study on the impact of deep building foundations on the groundwater system

Coastal areas are usually the preferred place of habitation for human beings. Anthropogenic activities such as the construction of high‐rise buildings and underground transport systems usually require extensive deep foundations and ground engineering works, which may unintentionally modify the coastal groundwater system because the construction materials of foundations are usually of low hydraulic conductivity. In this paper, the impact of these building foundations on the groundwater regime is studied using hypothetical flow and transport models. Various possible realizations of foundation distributions are generated using stochastic parameters derived from a topographical map of an actual coastal area in Hong Kong. The effective hydraulic conductivity is first calculated for different realizations and the results show that the effective hydraulic conductivity can be reduced significantly. Then a hypothetical numerical model based on FEFLOW is set up to study the change of hydraulic head, groundwater discharge, and saltwater‐fresh water interface. The groundwater level and flow are modified to various degrees, depending on the foundations percentage and the distribution pattern of the buildings. When the foundations percentage is high and the building foundations are aggregated, the hydraulic head is raised significantly and the originally one‐dimensional groundwater flow field becomes complicated. Seaward groundwater discharge will be reduced and some groundwater may become seepage through the ground surface. The transport model shows that, after foundations are added, overall the seawater and fresh groundwater interface moves landward, so extensive foundations may induce seawater intrusion. It is believed that the modification of the coastal groundwater system by building foundations may have engineering and environmental implications, such as submarine groundwater discharge, foundation corrosion, and slope stability. Copyright © 2007 John Wiley & Sons, Ltd.     

A Leaky Aquifer below Champlain Sea Clay: Closed‐Form Solutions for Natural Seepage

Closed‐form solutions are proposed for natural seepage in semiconfined (leaky) aquifers such as those existing below the massive Champlain Sea clay layers in the Saint‐Lawrence River Valley. The solutions are for an ideal horizontal leaky aquifer below an ideal aquitard that may have either a constant thickness and a constant hydraulic head at its surface, or a variable thickness and a variable hydraulic head at its surface. A few simplifying assumptions were needed to obtain the closed‐form solutions. These have been verified using a finite element method, which did not make any of the assumptions but gave an excellent agreement for hydraulic heads and groundwater velocities. For example, the difference between the two solutions was smaller than 1 mm for variations in the 5 to 8 m range for the hydraulic head in the semiconfined aquifer. Note that fitting the hydraulic head data of monitoring wells to the theoretical solutions gives only the ratio of the aquifer and aquitard hydraulic conductivities, a clear case of multiple solutions for an inverse problem. Consequently, field permeability tests in the aquitard and the aquifer, and pumping tests in the aquifer, are still needed to determine the hydraulic conductivity values.     

Modelling metolachlor exports in subsurface drainage water from two structured soils under maize (eastern France)

The study investigated the processes involved in metolachlor transport in two artificially drained, structured soils in eastern France. Measured losses of bromide and metolachlor in drainage water were compared with results simulated by the mechanistic, stochastic AgriFlux model. Simulated drainage water volumes were generally similar to the measured volumes when the spatial variability of the soil water properties was taken into account. When such variability was disregarded, cumulative water volumes of the clay soil were over- or underestimated by more than 20%. Two types of adsorption were tested. For instantaneous, reversible adsorption, using the partition coefficient Koc, metolachlor losses were underestimated in the first drainage water volumes and overestimated for the total study period. The use of slow adsorption and desorption kinetics (ADK) produced an export pattern similar to the observed one. A sensitivity analysis indicated that the simulated results are very sensitive to the values of the ADK rates, especially for the silty loam soil. The effect of ADK on the attenuation of metolachlor exports was more significant than the effect of degradation (2.3 and 6.7 times higher for the clay and silty loam soils, respectively). For the same four-month period, the bromide and metolachlor losses (using ADK) in the clay soil were 2.1 and 1.3 times greater, respectively, if the macroporosity was set at 10% than if it was not simulated. Conversely, macroporosity did not significantly affect these losses in the silty loam. The main factors involved in the metolachlor transport in the studied soils were: (i) the macroporosity, especially in the clay soil because of the low hydraulic conductivity of the matrix and (ii) the sorption kinetics rates which varied according to the soil physico-chemical characteristics.