Modelling of the groundwater hydrological behaviour of the Langueyú creek basin by using N‐way multivariate methods

The groundwater hydrochemical behaviour of the Langueyú creek basin (Argentina) has been evaluated through a systematic survey, followed by application of hydrological and chemometric multivariate techniques. Ten physicochemical parameters were determined in groundwater samples collected from 26 wells during four sampling campaigns (June 2010; October 2010; February 2011 and June 2011), originating a tridimensional experimental dataset X . Univariate statistical and graphical hydrochemical tools (contour maps and Piper diagrams) applied to individual campaigns, allowed to reach some preliminary conclusions. However, a best visualization of the aquifer behaviour was achieved by applying Principal Component Analysis (MA‐PCA) and N‐way PCA procedures, Parallel Factor Analysis and Tucker3. Results were consistent with two‐term models, being Tucker3 [2 2 1] the most adequate, explaining a large amount of the dataset variance (50.7%) with a low complexity. The first Tucker3 [1 1 1] interaction (38.2% of variance) is related with (i) calcium/magnesium versus sodium/potassium ion exchange processes; (ii) an increase of ionic concentration and (iii) a decrease of nitrate pollution, all processes along the direction of the groundwater flow. The second [2 2 1] interaction (12.5% of variance), accounts for the predominant role played by conductivity, bicarbonate and magnesium in the dataset. The seasonal variations are closely related to concentration/dilution phenomena originated by the variations of the phreatic levels, although this point will require additional sampling to establish a definitive hydrochemical model. Copyright © 2013 John Wiley & Sons, Ltd.     

Influence of aquifer properties and the spatial and temporal distribution of recharge and abstraction on sustainable yields in semi‐arid regions

This paper aims to contribute to understanding the importance of four factors on the determination of sustainable yields: (i) aquifer properties; (ii) temporal distribution of recharge; (iii) temporal distribution of groundwater pumping; and (iv) spatial distribution of pumping wells. It is important to comprehend how the present‐day and future vulnerability of groundwater systems to pumping activities depend on these critical factors and what the risks are of considering sustainable yield as a fixed percentage of mean annual recharge (MAR). A numerical model of the Querença–Silves aquifer in Portugal is used to develop hypothetical scenarios with which these factors are studied. Results demonstrate the aquifer properties, particularly the storage coefficient, have an important role in determining the resilience of an aquifer and therefore to which degree it is dependent on the spatial and temporal distribution of abstraction and recharge, as well as the occurrence of extreme events. Sustainable yields are determined for the developed scenarios based on specific criteria rather than a fraction of MAR. Under simplified current recharge and abstraction conditions, the sustainable yield was determined at approximately 73% of MAR or 76 million m³. When considering a concentration of rainfall in time, as predicted by climate scenarios for the region, sustainable yield could drop to ca 70% of MAR. However, a more even distribution of pumping volumes throughout the year could increase this value. The location of the pumping wells is seen to affect the distribution of hydraulic heads in the aquifer, albeit without significant changes in sustainable yield. Copyright © 2011 John Wiley & Sons, Ltd.     

Geophysical surveys for identifying saline groundwater in the semi‐arid region of the central Altiplano,Bolivia

In the central part of the Bolivian Altiplano, the shallow groundwater presents electrical conductivities ranging from 0·1 to 20 mS/cm. In order to study the origin of this salinity pattern, a good knowledge is required of the geometry of the aquifer at depth. In this study, geophysics has been used to complement the sparse data available from drill holes. One hundred time‐domain electromagnetic (TDEM) soundings were carried out over an area of 1750 km². About 20 geological logs were available close to some of the TDEM soundings. Three intermediate results were obtained from the combined data: (i) the relationship between the electrical conductivity of the groundwater and the formation resistivity, (ii) geoelectrical cross‐sections and (iii) geoelectrical maps at various depths. The limited data set shows a relationship between resistivity and the nature of the rock. From the cross‐sections, a conductive substratum with a resistivity of less than 1 Ω·m was identified at most of the sites at depths ranging from 50 to 350 m. This substratum could be a clay‐rich formation containing brines. Using derived relationships, maps of the nature of the formation (sandy, intermediate and clayey sediments) were established at depths of 10 and 50 m. Discrimination between sand and clays was impossible where groundwater conductivity is high (>3 mS/cm). In the central part of the area, where the groundwater conductivity is low, sandy sediments are likely to be present from the surface to a depth of more than 200 m. Clayey sediments are more likely to be present in the south‐east and probably constitute a hydraulic barrier to groundwater flow. In conclusion, the study demonstrates the efficiency of the TDEM sounding method to map conductive zones. Copyright © 2001 John Wiley & Sons, Ltd.     

GIS based DRASTIC model for groundwater vulnerability assessment: Case study of the shallow mio-plio-quaternary aquifer (Southeastern Tunisia)

Groundwater is the main source of water in arid regions. Thus, groundwater pollution becomes a major issue due to the increasing contamination, which poses serious and harmful risk to the environment. Groundwater vulnerability maps can be used as a tool to help decision makers to protect groundwater resources from contamination. The vulnerability of the Mio-Plio-Quaternary shallow aquifer (Southeast Tunisia) has been assessed using a DRASTIC model based on Geographic Information System (GIS). The different parameters of the model were collected from several sources and converted into thematic maps using ArcGis©. Each DRASTIC parameter was assigned a weight and rating based on a range of information within the parameter. Groundwater vulnerability map shows a large area (48%) with high risk of pollution. It indicates that the Southern part of the aquifer and the wadi beds are the most susceptible to contamination. The measured nitrate concentration is coherent with the DRASTIC model results.     

Modelling of the groundwater flow and of tracer movement in the porous and fissured media: Chalk Aquifer (Northern part of Paris Basin,France)

A groundwater flow model has been developed in order to study the chalk aquifer of Paris Basin, based on most of the geological and hydrological available data. The numerical processes are intended to modelling the groundwater flow in the Senonian (Late Cretaceous) formations and to visualize the tracer movement in groundwater resources in the experimental site of LaSalle Beauvais (northern part Paris Basin). Both objectives were achieved as follows: (i) the comprehension of the spatial distribution of the hydraulic conductivity in the chalk aquifer taking into account the characteristics of the hydrogeological system and (ii) the use of the analytical solution for describing one‐dimensional to two‐dimensional solute transport in a unidirectional steady‐state flow tracer with scale‐dependent dispersion. Advection and diffusion mechanisms are taken into account. Comparison between the breakthrough curves of the analytical and the numerical solutions provided an excellent agreement for various ranges of scale‐related transport parameters of interest. The developed power series solution facilitates fast prediction of the breakthrough curves at each observation point. Thus, the derived new solutions are widely applicable and are very useful for the validation of numerical transport. The numerical approach is carried out by MT3DMS, a Modular 3‐D Multi‐Species Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems, and based on total variation‐diminishing method using the ULTIMATE algorithm. The estimation of the infected surface could constitute an approach in water management and allows to prevent the risks of pollution and to manage the groundwater resource from a durable development perspective. Copyright © 2016 John Wiley & Sons, Ltd.     

Monitoring of drawdown pattern during pumping in an unconfined physical aquifer model

In this study, we attempted to analyse a drawdown pattern around a pumping well in an unconfined sandy gravelly aquifer constructed in a laboratory tank by means of both experimental and numerical modelling of groundwater flow. The physical model consisted of recharge, aquifer and discharge zones. Permeability and specific yield of the aquifer material were determined by Dupuit approximation under steady‐state flow and stepwise gravitational drainage of groundwater, respectively. The drawdown of water table in pumping and neighbouring observation wells was monitored to investigate the effect of no‐flow boundary on the drawdown pattern during pumping for three different boundary conditions: (i) no recharge and no discharge with four no‐flow boundaries (Case 1); (ii) no recharge and reservoir with three no‐flow boundaries (Case 2); (iii) recharge and discharge with two no‐flow boundaries (Case 3). Based on the aquifer parameters, numerical modelling was also performed to compare the simulated drawdown with that observed. Results showed that a large difference existed between the simulated drawdown and that observed in wells for all cases. The reason for the difference could be explained by the formation of a curvilinear type water table between wells rather than a linear one due to a delayed response of water table in the capillary fringe. This phenomenon was also investigated from a mass balance study on the pumping volume. The curvilinear type of water table was further evidenced by measurement of water contents at several positions in the aquifer between wells using time domain reflectometry (TDR). This indicates that the existing groundwater flow model applicable to an unconfined aquifer lacks the capacity to describe a slow response of water table in the aquifer and care should be taken in the interpretation of water table formation in the aquifer during pumping. Copyright © 2001 John Wiley & Sons, Ltd.     

Evaluation of Bias Associated with Capture Maps Derived from Nonlinear Groundwater Flow Models

The impact of groundwater withdrawal on surface water is a concern of water users and water managers, particularly in the arid western United States. Capture maps are useful tools to spatially assess the impact of groundwater pumping on water sources (e.g., streamflow depletion) and are being used more frequently for conjunctive management of surface water and groundwater. Capture maps have been derived using linear groundwater flow models and rely on the principle of superposition to demonstrate the effects of pumping in various locations on resources of interest. However, nonlinear models are often necessary to simulate head‐dependent boundary conditions and unconfined aquifers. Capture maps developed using nonlinear models with the principle of superposition may over‐ or underestimate capture magnitude and spatial extent. This paper presents new methods for generating capture difference maps, which assess spatial effects of model nonlinearity on capture fraction sensitivity to pumping rate, and for calculating the bias associated with capture maps. The sensitivity of capture map bias to selected parameters related to model design and conceptualization for the arid western United States is explored. This study finds that the simulation of stream continuity, pumping rates, stream incision, well proximity to capture sources, aquifer hydraulic conductivity, and groundwater evapotranspiration extinction depth substantially affect capture map bias. Capture difference maps demonstrate that regions with large capture fraction differences are indicative of greater potential capture map bias. Understanding both spatial and temporal bias in capture maps derived from nonlinear groundwater flow models improves their utility and defensibility as conjunctive‐use management tools.     

Recharge and Groundwater Flow Within an Intracratonic Basin,Midwestern United States

The conservative nature of chloride (Cl^?) in groundwater and the abundance of geochemical data from various sources (both published and unpublished) provided a means of developing, for the first time, a representation of the hydrogeology of the Illinois Basin on a basin‐wide scale. The creation of Cl^? isocons superimposed on plan view maps of selected formations and on cross sections across the Illinois Basin yielded a conceptual model on a basin‐wide scale of recharge into, groundwater flow within and through the Illinois Basin. The maps and cross sections reveal the infiltration and movement of freshwater into the basin and dilution of brines within various geologic strata occurring at basin margins and along geologic structures. Cross‐formational movement of brines is also seen in the northern part of the basin. The maps and cross sections also show barriers to groundwater movement created by aquitards resulting in areas of apparent isolation/stagnation of concentrated brines within the basin. The distribution of Cl^? within the Illinois Basin suggests that the current chemical composition of groundwater and distribution of brines within the basin is dependent on five parameters: (1) presence of bedrock exposures along basin margins; (2) permeability of geologic strata and their distribution relative to one another; (3) presence or absence of major geologic structures; (4) intersection of major waterways with geologic structures, basin margins, and permeable bedrock exposures; and (5) isolation of brines within the basin due to aquitards, inhomogeneous permeability, and, in the case of the deepest part of the basin, brine density effects.     

A large weighing lysimeter for evapotranspiration and soil‐water–groundwater exchange studies

A large weighing lysimeter was installed at Yucheng Comprehensive Experimental Station, north China, for evapotranspiration and soil‐water–groundwater exchange studies. Features of the lysimeter include the following: (i) mass resolution equivalent to 0·016 mm of water to accurately and simultaneously determine hourly evapotranspiration, surface evaporation and groundwater recharge; (ii) a surface area of 3·14 m² and a soil profile depth of 5·0 m to permit normal plant development, soil‐water extraction, soil‐water–groundwater exchanges, and fluctuations of groundwater level; (iii) a special supply–drainage system to simulate field conditions of groundwater within the lysimeter; (iv) a soil mass of about 30 Mg, including both unsaturated and saturated loam. The soil consists mainly of mealy sand and light loam. Monitoring the vegetated lysimeter during the growing period of winter wheat, from October 1998 through to June 1999, indicated that during the period groundwater evaporation contributed 16·6% of total evapotranspiration for a water‐table depth from 1·6 m to 2·4 m below ground surface. Too much irrigation reduced the amount of upward water flow from the groundwater table, and caused deep percolation to the groundwater. Data from neutron probe and tensiometers suggest that soil‐water‐content profiles and soil‐water‐potential profiles were strongly affected by shallow groundwater. Copyright © 2000 John Wiley & Sons, Ltd.     

A supervised committee machine artificial intelligent for improving DRASTIC method to assess groundwater contamination risk: a case study from Tabriz plain aquifer,Iran

Vulnerability maps are designed to show areas of greatest potential for groundwater contamination on the basis of hydrogeological conditions and human impacts. The objective of this research is (1) to assess the groundwater vulnerability using DRASTIC method and (2) to improve the DRASTIC method for evaluation of groundwater contamination risk using AI methods, such as ANN, SFL, MFL, NF and SCMAI approaches. This optimization method is illustrated using a case study. For this purpose, DRASTIC model is developed using seven parameters. For validating the contamination risk assessment, a total of 243 groundwater samples were collected from different aquifer types of the study area to analyze \( {\text{NO}}_{ 3}^{ - } \) concentration. To develop AI and CMAI models, 243 data points are divided in two sets; training and validation based on cross validation approach. The calculated vulnerability indices from the DRASTIC method are corrected by the \( {\text{NO}}_{3}^{ - } \) data used in the training step. The input data of the AI models include seven parameters of DRASTIC method. However, the output is the corrected vulnerability index using \( {\text{NO}}_{3}^{ - } \) concentration data from the study area, which is called groundwater contamination risk. In other words, there is some target value (known output) which is estimated by some formula from DRASTIC vulnerability and \( {\text{NO}}_{3}^{ - } \) concentration values. After model training, the AI models are verified by the second \( {\text{NO}}_{3}^{ - } \) concentration dataset. The results revealed that NF and SFL produced acceptable performance while ANN and MFL had poor prediction. A supervised committee machine artificial intelligent (SCMAI), which combines the results of individual AI models using a supervised artificial neural network, was developed for better prediction of vulnerability. The performance of SCMAI was also compared to those of the simple averaging and weighted averaging committee machine intelligent (CMI) methods. As a result, the SCMAI model produced reliable estimates of groundwater contamination risk.     

Groundwater Deterioration of Shallow Groundwater Aquifers Due to Overexploitation in Northeast Jordan

Groundwater is the major water resource in Jordan and most of the groundwater basins are already exploited beyond their estimated safe yield. Azraq basin is one of the most important groundwater basins in Jordan, which supplies Amman with drinking water. However, due to overpumping from the shallow groundwater aquifers, the water level dropped dramatically and signs of salinization and depletion are starting to occur. The severe drawdown in the Azraq well‐field caused a reverse in the hydraulic gradient and consequently, the saltwater in the center of the basin (Qa‐Azraq) started to move in the direction of the well‐field. The salinization in the shallow aquifer (basalt/B5/B4) is believed to result from one of the following scenarios: (i) a reverse flow from Sabkha to the AWSA well field, (ii) an upward leakage from the middle aquifer system (B2/A7) and the combined B3 Aquitard‐B2/A7 aquifer, (iii) a dissolution process between the water and rock matrix due to lowering of the dynamic water levels during pumping which reached the mineralized formations underlying the Basalt. The salinization trend of some AWSA wells represented by the gradual increase of major ions is associated with rather constant stable isotopic contents. This indicates that these constituents originate from the main minerals existing in the matrix of the aquifers and thus this scenario is the most likely to occur.     

Modification of the DRASTIC Framework for Mapping Groundwater Vulnerability Zones

The DRASTIC technique is commonly used to assess groundwater vulnerability. The main disadvantage of the DRASTIC method is the difficulty associated with identifying appropriate ratings and weight assignments for each parameter. To mitigate this issue, ratings and weights can be approximated using different methods appropriate to the conditions of the study area. In this study, different linear (i.e., Wilcoxon test and statistical approaches) and nonlinear (Genetic algorithm [GA]) modifications for calibration of the DRASTIC framework using nitrate (NO₃) concentrations were compared through the preparation of groundwater vulnerability maps of the Meshqin-Shahr plain, Iran. Twenty-two groundwater samples were collected from wells in the study area, and their respective NO₃ concentrations were used to modify the ratings and weights of the DRASTIC parameters. The areas found to have the highest vulnerability were in the eastern, central, and western regions of the plain. Results showed that the modified DRASTIC frameworks performed well, compared to the unmodified DRASTIC. When measured NO₃ concentrations were correlated with the vulnerability indices produced by each method, the unmodified DRASTIC method performed most poorly, and the Wilcoxon–GA–DRASTIC method proved optimal. Compared to the unmodified DRASTIC method with an R² of 0.22, the Wilcoxon–GA–DRASTIC obtained a maximum R² value of 0.78. Modification of DRASTIC parameter ratings was found to be more efficient than the modification of the weights in establishing an accurately calibrated DRASTIC framework. However, modification of parameter ratings and weights together increased the R² value to the highest degree.     

Evaluating the use of a gridded climate surface for modelling groundwater recharge in a semi‐arid region (Okanagan Basin,Canada)

Spatially distributed groundwater recharge was simulated for a segment of a semi‐arid valley using three different treatments of meteorological input data and potential evapotranspiration (PET). For the same area, timeframe, land cover characteristics and soil properties, groundwater recharge was estimate using (i) single‐station climate data with monthly PET calculated by the Thornthwaite method; (ii) single‐station climate data with daily PET calculated by the Penman–Monteith method; and (iii) daily gridded climate data with spatially distributed PET calculated using the Penman–Monteith method. For each treatment, the magnitude and distribution of actual evapotranspiration (AET) for summer months compared well with those estimated for a 5‐year crop study, suggesting that the near‐surface hydrological processes were replicated and that subsequent groundwater recharge rates are realistic. However, for winter months, calculated AET was near zero when using the Thornthwaite PET method. Mean annual groundwater recharge varied from ～3·2 to 10·0 mm when PET was calculated by the Thornthwaite method, and from ～1·8 to 7·5 mm when PET was calculated by the Penman–Monteith method. Comparisons of bivariate plots of seasonal recharge rates estimated from single‐station versus gridded surface climate reveal that there is greater variability between the different methods for spring months, which is the season of greatest recharge. Furthermore, these seasonal differences are shown to provide different results when compared to the depth to water table, which could lead to different results of evaporative extinction depth. These findings illustrate potential consequences of using different approaches for representing spatial meteorological input data, which could provide conflicting predictions when modelling the influence of climate change on groundwater recharge. Copyright © 2010 John Wiley & Sons, Ltd.     

Stochastic environmental risk analysis: an integrated methodology for predicting cancer risk from contaminated groundwater

Stochastic environmental risk assessment considers the effects of numerous biological, chemical, physical, behavioral and physiological processes that involve elements of uncertainty and variability. A methodology for predicting health risks to individuals from contaminated groundwater is presented that incorporates the elements of uncertainty and variability in geological heterogeneity, physiological exposure parameters, and in cancer potency. An idealized groundwater basin is used to perform a parametric sensitivity study to assess the relative impact of (a) geologic uncertainty, (b) behavioral and physiological variability in human exposure and (c) uncertainty in cancer potency on the prediction of increased cancer risk to individuals in a population exposed to contaminants in household water supplied from groundwater. A two-dimensional distribution (or surface) of human health risk was generated as a result of the simulations. Cuts in this surface (fractiles of variability and percentiles of uncertainty) are then used as a measure of relative importance of various model components on total uncertainty and variability. A case study for perchloroethylene or PCE, shows that uncertainty and variability in hydraulic conductivity play an important role in predicting human health risk that is on the same order of influence as uncertainty of cancer potency.     

Characterization and evaluation of processes governing the groundwater quality in parts of the Sabarmati basin,Gujarat using hydrochemistry integrated with GIS

Intense agricultural and industrial activities in any area are likely to make groundwater vulnerable with respect to its quality. In one such area which is a part of Sabarmati river basin of Gujarat, factors influencing the groundwater hydrochemistry in pre‐ and post‐monsoon season were evaluated. Groundwater samples were collected from 5 km × 5 km grids on the basis of spectral signature of vegetation and soil, observed on satellite image. Integration of Conventional graphical plots, Piper plot, saturation index values (estimated using PHREEQC) and GIS was helpful not only to create the database for analysis of spatial variation in respective water quality parameters but also to decipher the hydrogeochemical process occurring in such a large area. USSL diagram and % sodium were used to characterise the suitability of groundwater for irrigation. It was observed that leaching of wastes disposed from anthropogenic activities and agrichemicals is the major factor influencing the groundwater quality, in addition to the natural processes such as weathering, dissolution and ion exchange. Sea water relics are also impacting the groundwater quality. Control of indiscriminate and unplanned exploitation of groundwater, application of fertilizers and disposal of industrial wastes in the affected areas can possibly ensure groundwater protection from further pollution and depletion. Copyright © 2011 John Wiley & Sons, Ltd.     

Surveying Upstate NY Well Water for Pesticide Contamination: Cayuga and Orange Counties

Private wells in Cayuga and Orange counties in New York were sampled to determine the occurrence of pesticide contamination of groundwater in areas where significant pesticide use coincides with shallow or otherwise vulnerable groundwater. Well selection was based on local groundwater knowledge, risk modeling, aerial photo assessments, and pesticide application database mapping. Single timepoint samples from 40 wells in each county were subjected to 93‐compound chromatographic scans. All samples were nondetects (reporting limits ≤1 μg/L), thus no wells from either county exceeded any of 15 state groundwater standards or guidance values. More sensitive enzyme‐linked immunosorbent assays (ELISA) found two wells with quantifiable atrazine in each county (0.1–0.3 μg/L), one well with quantifiable diazinon (0.1 μg/L) in Orange County, and one well with quantifiable alachlor (0.2 μg/L) in Cayuga County. Trace detections (<0.1 μg/L) in Cayuga County included atrazine (five wells), metolachlor (six wells), and alachlor (one well), including three wells with multiple detections. All 12 Cayuga County wells with ELISA detections had either corn/grain or corn/forage rotations as primary surrounding land uses (although 20 other wells with the same land uses had no detections) and all quantified detections and most trace detections occurred in wells up to 9‐m deep. Orange County trace (<0.1 μg/L) ELISA detections (atrazine three wells, diazinon one well, and metolachlor five wells) and quantified detections were only generally associated with agricultural land uses. Finding acceptable drinking water quality in areas of vulnerable groundwater suggests that water quality in less vulnerable areas will also be good.     

Use of Nested Flow Models and Interpolation Techniques for Science‐Based Management of the Sheyenne National Grassland,North Dakota,USA

Noxious weeds threaten the Sheyenne National Grassland (SNG) ecosystem and therefore herbicides have been used for control. To protect groundwater quality, the herbicide application is restricted to areas where the water table is less than 10 feet (3.05 m) below the ground surface in highly permeable soils, or less than 6 feet (1.83 m) below the ground surface in low permeable soils. A local MODFLOW model was extracted from a regional GFLOW analytic element model and used to develop depth‐to‐groundwater maps in the SNG that are representative for the particular time frame of herbicide applications. These maps are based on a modeled groundwater table and a digital elevation model (DEM). The accuracy of these depth‐to‐groundwater maps is enhanced by an artificial neural networks (ANNs) interpolation scheme that reduces residuals at 48 monitoring wells. The combination of groundwater modeling and ANN improved depth‐to‐groundwater maps, which in turn provided more informed decisions about where herbicides can or cannot be safely applied.     

Assessment of groundwater pollution risk as a characteristic of the stability of its quality

Definitions of hazard and risk of groundwater pollution are given. A deterministic method for the assessment of groundwater pollution risk using estimates of groundwater protection against and vulnerability to pollution and stability indicators of groundwater quality is considered. Also presented are the principal methodological approaches to the assessment of groundwater protection against pollution and the formation of the structure of indicators and indices characterizing the stability of groundwater quality. The structure of hazard, risks, and damages associated with groundwater pollution is shown. Expert appraisal method is used for the assessment of groundwater pollution risks.     

Geoelectrical structure and hydrogeological investigations of the southern Rif Cordillera (Morocco)

An electrical prospecting survey is conducted in the Rharb basin, a semi‐arid region in the southern part of the Rifean Cordillera (Morocco) to delineate characteristics of the aquifer and the groundwater affected by the marine intrusion related to Atlantic Ocean. Analysis and interpretations of electrical soundings, bi‐logarithmic diagrams and the geoelectrical sections highlight a monolayer aquifer in the southern part, a multilayer system in the northern part of the Rharb basin and lenticular semi‐permeable formations. Several electrical layers have been deduced from the analysis of bi‐logarithmic diagrams: resistant superficial level (R₀), conducting superficial level (C₀), resistant level (R), intermediary resistant level (R′), conducting level (Cp) and intermediary layer of resistivity (AT). Spatial distribution of the resistivity deduced from the interpretation of apparent resistivity maps (AB = 400 and 1000 m) and the decreasing of resistivity values (35–10 Ωm), in particular in the coastal zone show that this heterogeneity is related to several anomalies identified in the coastal area, which result from hydraulic and geological processes: (i) heterogeneous hydraulic conductivity in particular in the southern part of the Rharb; (ii) lateral facies and synsedimentary faulting and (iii) the relationship between the electrical conductivity and chloride concentration of groundwater shows that salinity is the most important factor controlling resistivity. The distribution of fresh/salt‐water zones and their variations in space along geoelectrical sections are established through converting subsurface depth‐resistivity models. Copyright © 2010 John Wiley & Sons, Ltd.     

Using chlorofluorocarbons (CFCs) and tritium to improve conceptual model of groundwater flow in the South Coast Aquifers of Laizhou Bay,China

The southern coastal plain of Laizhou Bay, which is the area most seriously affected by salt water intrusion in north China, is a large alluvial depression, which represents one of the most important hydrogeological units in the coastal region of northern China. Chlorofluorocarbons (CFCs, including CFC‐11, CFC‐12 and CFC‐113) and tritium were used together for dating groundwater up to 50 years old in the study area. There are two cones of depression, caused by intensive over‐exploitation of fresh groundwater in the south and brine water in the north. The assigned CFC apparent ages for shallow groundwater range from 8 a to >50 a. A binary mixing model based on CFC‐113 and CFC‐12 concentrations in groundwater was used to estimate fractions of young and pre‐modern water in shallow aquifers and to identify groundwater mixing processes during saltwater intrusion. Discordance between concentrations of different CFC compounds indicate that shallow groundwater around the Changyi cone of depression is vulnerable to contamination. Pumping activities, CFC contamination, mixing and/or a large unsaturated zone thickness (e.g. >20 m) may be reasons for some groundwater containing CFCs without tritium. Saline intrusion mainly occurs because of large head gradients between fresh groundwater in the south and saline water bodies in the north, forming a wedge of saline water below/within fresh aquifer layers. Both CFC and tritium dates indicate that the majority of the saline water is from >50 a, with little or no modern seawater component. Based on the distribution of CFC apparent ages, tritium contents plus chemical and physical data, a conceptual model of groundwater flow along the investigated Changyi‐Xiaying transect has been developed to describe the hydrogeological processes. Three regimes are identified from south to north: (i) fresh groundwater zone, with a mixing fraction of 0.80–0.65 ‘young’ water calculated with the CFC binary mixing model (groundwater ages <34 a) and 1.9–7.8TU of tritium; (ii) mixing zone characterized by a mixing fraction of 0.05–0.65 young groundwater (ages of 23–44 a), accompanied by local vertical recharge and upward leakage of older groundwater; and (iii) salt water zone, mostly comprising waters with ages beyond the dating range of both CFCs and tritium. Some shallow groundwater in the north of the Changyi groundwater depression belongs to the >50a water group (iii), indicating slow velocity of groundwater circulation and possible drawing in of saline or deep groundwater that is tracer‐free. Copyright © 2012 John Wiley & Sons, Ltd.