Hydrogeological framework and estimation of aquifer hydraulic parameters using geoelectrical data: a case study from West Iran

Detailed local geological, geophysical, and hydrogeological investigations were carried out for the alluvial aquifer in the Kangavar basin, West Iran to delineate the architecture of different subsurface geological horizons using lithologs and generated vertical electrical sounding (VES) data. An attempt has also been made to estimate aquifer transmissivity from resistivity data. Forty VESs were recorded with the Schlumberger electrode configuration in the study area; 28 of these were selected for evaluation. The maximum current electrode spacing was 400–500 m. The data obtained were interpreted by computer iterative modeling with curve matching for calibration purposes. In order to ascertain the subsurface geological framework, the general distribution of resistivity responses of the geological formations was obtained and geoelectrical sections along a number of lines were prepared. Probable aquifer horizons from these sections were identified. The transmissivity of the unconfined aquifer was computed by determining the Dar-Zarrouk parameters (longitudinal unit conductance and transverse unit resistance) and were compared with the actual field transmissivity. The results showed a direct relation between aquifer transmissivity and modified transverse resistance.     

Effects of shallow and deep geology on seismic hazard estimates: a case study of pseudo-acceleration response spectra for the northwestern Balkans

Northwestern Balkans represents one of the rare regions where data on both the shallow geology (i.e., local soil conditions) and the deep geology exist for stations that recorded hundreds of strong motion records. The strong motion database used in this study consists of 203 strong motion accelerograms (each with three orthogonal translational components), recorded in former Yugoslavia in the period 1976–1983, from 108 contributing earthquakes. In this paper, the results of a series of regression analyses are presented where empirical equations for scaling pseudo-acceleration response spectra were developed on the basis of 5 subsets of data and using three prediction models. Results of the regression analyses show that for ground motion in the horizontal direction, both the shallow and deep geology site conditions have to be taken into account or else the spectral estimates might be considerably biased. Results show that the shallow geology influences spectral amplitudes the most in the short period range and has much lesser effects for larger periods, while the influence of the deep geology is expressed in a wider range of intermediate to long periods. Results also show that if the prediction model that considers solely the shallow geology effects is used, the spectral peaks that have been related to the deep geology effects will completely vanish for larger source-to-site distances, while in the case when solely the deep geology is considered, the peaks related to the shallow geology will not be visible for any distance. As for the amplitudes in the vertical direction, although both the shallow and the deep geology effects are less expressed than in the horizontal one, they still cannot be neglected—the deep geology effects are visible for a wider range of vibration periods, while the shallow geology has some effects only for periods smaller than ~0.3 s.     

Regional modeling of geothermal energy systems in shallow aquifers: the Leibnitzer Feld case study (Austria)

Large thermal extractions and extensive implementation of groundwater heat pumps (GWHP) necessitate a validation of the sustainability of their use and possible detrimental effects on groundwater. The goal of this work is to develop a regional heat transport model (of ~13 km × 5 km) for real site conditions. This model should consider all relevant transport processes, despite the large area under investigation. The model is based on a two-dimensional, transient-calibrated groundwater flow model for the “Leibnitzer Feld” (Styria, Austria). The two-dimensional horizontal model is linked via the FEFLOW interface manager with a newly developed “Multi-Layer-Model”-tool, which reproduces thermal aquifer–atmosphere interaction. Based on the regional heat transport model, scenarios are delineated for heating and cooling purposes for large GWHPs, which are appropriate for a small manufacturing business, an administrative building and 10 family homes. First of all, these have large spacing and thereafter, effects of area-covering usage of geothermal systems are evaluated for five administrative buildings located in close proximity to one another (200–350 m) and also for a large number of smaller heat extractions (each representing a one family house system). Modeled spatial and temporal temperature effects on the shallow aquifer are discussed. It was possible to present a simulation of realistic heating and cooling scenarios. This simulation may be introduced into practice once some further simplifications to the system are made. Locally limited heat plumes (max. length: 625 m) were observed for the manufacturing business. Any thermal effects coming from the geothermal systems were shown to be temporally stable. As such, no distinct trend of reduced annual temperatures could be observed.     

Use of geochemical modeling to evaluate the hydraulic connection of aquifers: a case study from Chianan Plain, Taiwan

Aquifers are generally composed of highly permeable layers that can conduct a considerable amount of groundwater. Traditionally, aquifer units are correlated through the concept of lithostratigraphy. For low-permeable aquifers, it is difficult to define the spatial distribution of hydrogeological units, and this study attempts to use geochemical modeling to identify the groundwater flow paths in an area of Taiwan. Multiple geochemical analyses, including groundwater chemistry; stable isotopic compositions of hydrogen, oxygen and carbon; and radiocarbon contents were performed. Using these parameters as the constraints of geochemical models, the hydraulic connection was examined between pairs of possibly interlinked wells along four selected cross sections, and the conceptual groundwater model was accordingly established. The resultant model suggests that the hydraulic connection between aquifers should be correlated with the concept of chronological stratigraphy, especially for low-permeable, unconsolidated aquifers. Using Darcy’s law, the hydraulic conductivities of the fine-sand aquifers were estimated to be between 3.14 × 10⁻⁵ and 1.80 × 10⁻⁴ m/s, which are roughly one order of magnitude higher than those derived by in situ pumping tests. The substantial extraction of groundwater over a long period in the studied area could accelerate groundwater flow, leading to an overestimation of the aquifer permeability.

---

Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.

---

Résumé En règle générale, les aquifères sont constitués de couches très perméables capables de conduire des quantités d’eau considérables. Les unités aquifères sont traditionnellement corrélées d’un point de vue lithostratigraphique. Dans le cas des aquifères à faible perméabilité, il s’avère difficile de définir une distribution spatiale des unités hydrogéologiques, et cette étude tente d’utiliser la modélisation géochimique pour identifier les cheminements des eaux souterraines dans un secteur de Taiwan. Plusieurs analyses géochimiques ont été réalisées ; elles incluent la chimie des eaux souterraines, les isotopes stables de l’hydrogène, de l’oxygène et du carbone et l’abondance en carbone 14. En utilisant ces paramètres comme contraintes pour les modèles géochimiques, les liaisons hydrauliques entre les puits potentiellement interconnectés ont été étudiées deux par deux, le long de quatre coupes présélectionnées, et le modèle conceptuel a été établi en conséquence. Ce modèle résultant suggère que les connexions hydrauliques inter-aquifères devraient être corrélées avec les concepts stratigraphiques, surtout pour les aquifères non consolidés à faible perméabilité. Les perméabilités dans les aquifères de sables fins ont été estimées par la loi de Darcy entre 3.14 × 10⁻⁵ et 1.80 × 10⁻⁴ m/s, soit environ un ordre de grandeur au-dessus de celles issues des tests de pompage sur site. L’exploitation substantielle des eaux souterraines sur une longue période dans la zone d’étude peut accélérer les écoulements souterrains, menant à une surestimation de la perméabilité de l’aquifère.

---

Resumen Los acuíferos están generalmente compuestos de capas muy permeables que pueden conducir una cantidad considerable de agua subterránea. Tradicionalmente, se ponen en correlación las unidades acuíferas a través del concepto de litoestratigrafía. Para los acuíferos de baja-permeabilidad, es difícil de definir la distribución espacial de unidades hidrogeológicas y este estudio intenta usar modelamiento geoquímico para identificar las direcciones del flujo de agua subterránea en un área de Taiwán. Se realizó el análisis geoquímico múltiple, incluyendo la química del agua subterránea, la composición de isótopos estables de hidrógeno, oxígeno y carbono, y el contenido de radiocarbono. Usando estos parámetros como limitantes de modelos geoquímicos, la conexión hidráulica se examinó entre pares de pozos posiblemente interrelacionados, a lo largo de cuatro cortes transversales seleccionados y de acuerdo con esto se estableció el modelo conceptual del agua subterránea. El modelo resultante sugiere que la conexión hidráulica entre los acuíferos deba interrelacionarse con el concepto de estratigrafía cronológica, sobre todo para los acuíferos sin consolidar de baja permeabilidad. Usando la ley de Darcy, se estimaron las conductividades hidráulicas de los acuíferos de arena fina obteniendo valores entre 3.14 × 10⁻⁵ y 1.80 × 10⁻⁴ m/s, las cuales son aproximadamente un orden de magnitud mayor, que aquéllos obtenidos in situ por las pruebas de bombeo. La extracción considerable de agua subterránea durante un periodo largo en el área estudiada, podría acelerar el flujo de agua subterránea, llevando a una sobrestimación de la permeabilidad en los acuíferos.

     

黄土缓坡片蚀过程及其水力参数适宜性试验研究

为确定黄土缓坡片蚀过程中侵蚀限制条件的变化及水力参数的适宜性,采用模拟降雨试验,在60 mm/h雨强、3种坡度（7.5°、10°、15°）和2种坡长（5 m和10 m）条件下,验证含沙量作为表征参量判定片蚀过程侵蚀限制条件的可行性,提出检验水力参数与产沙量作用关系适宜性的评价方法。结果表明:①依据含沙量随径流的变化规律,片蚀过程依次经历输沙能力限制阶段、剥蚀能力限制阶段Ⅰ和剥蚀能力限制阶段Ⅱ共3个侵蚀阶段,其中剥蚀能力限制阶段Ⅰ的侵蚀量占绝对优势。②从水力参数的适宜性来看,径流功率和雷诺数均能很好地表征与产沙量的关系,径流剪切力和弗劳德数则不能有效反映与产沙量的关系。③在不同侵蚀限制阶段,径流功率和雷诺数与产沙量的作用关系是不同的。输沙能力限制阶段和剥蚀能力限制阶段Ⅱ时,径流功率和雷诺数与产沙量分别呈指数关系和线性正相关关系;剥蚀能力限制阶段Ⅰ时,该2种水力参数均与产沙量呈线性负相关关系,但作用关系不唯一,且存在雷诺数临界值,即在输沙能力限制阶段和剥蚀能力限制阶段Ⅱ,径流功率和雷诺数与产沙量的作用关系不受坡长效应影响;在剥蚀能力限制阶段Ⅰ,需要同时考虑坡长和雷诺数临界值的共同影响。该研究结果可为深入理解片蚀过程的复杂性提供参考依据。     

Determination of groundwater availability in shallow arid region aquifers utilizing GIS technology: a case study in Hada Al-Sham, Western Saudi Arabia

A simple, physically based method is developed in this paper to assist in the allocation of areas with high groundwater potential and for the determination of maximum allowed pumping rate to ensure proper groundwater management. This method utilizes the aquifer physical properties as well as GIS technology to accomplish this purpose. The design of this method was considered to be applicable in areas with little data, such as in most arid regions. This technique was applied to a catchment in an arid environment where qualitative as well as quantitative analyses of the results were undertaken. Locations of available groundwater and rates of maximum allowable pumping were compared with observations and experiments in the field and a good agreement was found. It was concluded that the best groundwater location was in the alluvial area, which represents only 16% of the total aquifer, which is a typical case in arid region catchments. The rate of maximum pumping was estimated to be 65 m³/h. However, to benefit 55% of the area, the maximum pumping rate should only be 40 m³/h with an average rate throughout the area (55%) of about 24 m³/h.This revised version was published online in December 2004 with corrections to the category.     

Modified DRASTIC assessment for intrinsic vulnerability mapping of karst aquifers: a case study

Groundwater in karstic aquifers can be dangerously sensitive to contamination. In this paper, DRASTIC assessment was modified and applied, for the first time, to address the intrinsic vulnerability for karst aquifers. The theoretical weights of two of DRASTIC’s parameters (aquifer media and hydraulic conductivity) were modified through sensitivity analysis. Two tests of sensitivity analyses were carried out: the map removal and the single parameter sensitivity analyses. The modified assessment was applied for the karst aquifers underlying Ramallah District (Palestine) as a case study. The aquifer vulnerability map indicated that the case study area is under low, moderate and high vulnerability of groundwater to contamination. The vulnerability index can assist in the implementation of groundwater management strategies to prevent degradation of groundwater quality. The modified DRASTIC assessment has proven to be effective because it is relatively straightforward, use data that are commonly available or estimated and produces an end product that is easily interpreted.     

Evaluation of hydrogeochemical parameters of groundwater for suitability of domestic and irrigational purposes: a case study from central Ganga Plain,India

The study area is one of the most productive part of central Ganga Plain and fertile tract for sugarcane cultivation. Hydrogeochemical parameters of groundwater of the study area was evaluated to know the suitability of groundwater for domestic and irrigational purposes. Fifty-five groundwater samples from hand pumps in pre-monsoon 2007 were analyzed for physical and chemical groundwater parameters (electrical conductivity, pH, total dissolved solid, Na, K, Ca, Mg, HCO₃, Cl, SO₄, and NO₃). Three major groups of groundwater with distinct chemical compositions had been identified on L–L diagram, i.e., Ca⁺Mg-HCO₃ type, mixed type, and alkali bicarbonate type. All possible species, such as Na-Cl, K-Cl, Na-HCO₃, Na-SO₄, Ca-HCO₃, Mg-HCO₃, Ca-SO₄, and Mg-SO₄ occur in the groundwater system. Groundwater comes under the category of moderately hard to very hard, mildly acidic to slightly alkaline in nature. Majority of the sample are within the permissible limit when compared with drinking water standards in terms of electrical conductivity, pH, total dissolved solid, Na, K, Ca, Mg, HCO₃, Cl, SO₄, and NO₃. According to Gibb’s ratio, most of groundwater samples fall in the rock dominance field. Assessing the domestic uses, all the samples are considered fit, as they are neither acidic nor strongly alkaline. Based on analytical results, irrigational quality parameters like sodium adsorption ratio, residual sodium carbonate, and permeability index were calculated which indicate that the groundwater is also suitable for irrigational uses.     

Multivariate statistical analysis of the groundwater flow in shallow aquifers: a case of the basins of northern Algeria

Water resources in Algeria are mainly controlled by climate change which creates enormous problems in its planning, management and distribution. While the surface water resources are perfectly managed and operated by means of dams and small dams built for several years, the groundwater resources remain long unknown and unusable because of the lack of relevant working tools (e.g., methods, formulas, maps, etc.) for planners and engineers working in the field of water resources exploration. To highlight the hydrodynamic processes of groundwater in shallow aquifers of the basins of northern Algeria, we conducted a study using 81 subwatersheds collected from different locations at the basins; taking into account the climatic and geomorphological factors, to understand water usage trends, analyse patterns, tap good shallow aquifers and ensure long lasting supplies of water through arid periods, mapping and modelling of groundwater are fundamental to problem resolution. Multivariate statistical techniques as well as cluster and principal component analysis were applied to the data on groundwater flow, with the objective of defining the main controls on the groundwater flows at the basins. These statistical techniques showed the presence of three groundwater flow groups with increasing importance according to precipitation. The first group was mainly influenced by climatic factors, the second was more controlled by the communication between the surface and underground flows and the third group revealed the influence of geomorphological factors on groundwater flows.     

Vertical geochemical variations and arsenic mobilization in the shallow alluvial aquifers of the Chapai-Nawabganj District, northwestern Bangladesh: implication of siderite precipitation

Core sediments from three disturbed boreholes (JOR, GHAT, and RAJ) and two undisturbed boreholes (DW1 and DW2) were collected in the study area of the Chapai-Nawabganj district of northwestern Bangladesh for geochemical analyses. In the study area, groundwater samples from fourteen As-contained private wells and five nested piezometers at both the DW1 and DW2 boreholes were also collected and analyzed. The groundwater arsenic concentrations in the uppermost aquifer (10–40 m of depth) range from 3 to 315 μg/L (mean 47.73 ± 73.41 μg/L), while the arsenic content in sediments range from 2 to 14 mg/kg (mean 4.36 ± 3.34 mg/kg). An environmental scanning electron microscope (ESEM) with an energy dispersive X-ray spectrometer was used to investigate the presence of major and trace elements in the sediments. Groundwaters in the study area are generally the Ca–HCO₃ type with high concentrations of As, but low levels of Fe, Mn, NO₃ ^? and SO ₄ ^?2 . The concentrations of As, Fe, Mn decrease with depth in the groundwater, showing vertical geochemical variations in the study area. Statistical analysis clearly shows that As is closely associated with Fe and Mn in the sediments of the JOR core (r = 0.87, p < 0.05 for Fe and r = 0.78, p < 0.05 for Mn) and GHAT core (r = 0.95, p < 0.05 for Fe and r = 0.93, p < 0.05 for Mn), while As is not correlated with Fe and Mn in groundwater. The comparatively low Fe and Mn concentrations in some groundwater and the ESEM image revealed that siderite precipitated as a secondary mineral on the surface of the sediment particles. The correlations along with results of sequential extraction experiments indicated that reductive dissolution of FeOOH and MnOOH represents a mechanism for releasing arsenic into the groundwater.     

Salinity mapping of coastal groundwater aquifers using hydrogeochemical and geophysical methods: a case study from north Kelantan,Malaysia

Integrated hydrogeochemical and geophysical methods were used to study the salinity of groundwater aquifers along the coastal area of north Kelantan. For the hydrogeochemical investigation, analysis of major ion contents of the groundwater was conducted, and other chemical parameters such as pH and total dissolved solids were also determined. For the geophysical study, both geoelectrical resistivity soundings and reflection seismic surveys were conducted to determine the characteristics of the subsurface and groundwater contained within the aquifers. The pH values range from 6.2 to 6.8, indicating that the groundwater in the study area is slightly acidic. Low content of chloride suggests that the groundwater in the first aquifer is fresh, with an average concentration of about 15.8 mg/l and high geoelectrical resistivity (>45 ohm m). On the other hand, the groundwater in the second aquifer is brackish, with chloride concentration ranging from 500 mg/l to 3,600 mg/l and very low geoelectrical resistivity (<45 ohm m) as well as high concentration of total dissolved solids (>1,000 mg/l). The groundwater in the third aquifer is fresh, with chloride concentrations generally ranging from 2 mg/l to 210 mg/l and geoelectrical resistivity of greater than 45 ohm m. Fresh and saltwater interface in the first aquifer is generally located directly in the area of the coast, but, for the second aquifer, both hydrogeochemical and geoelectrical resistivity results indicate that the fresh water and saltwater interface is located as far as 6 km from the beach. The considerable chloride ion content initially suggests that the salinity of the groundwater in the second aquifer is probably caused by the intrusion of seawater. However, continuous monitoring of the chloride content of the second aquifer indicated no significant changes with time, from which it can be inferred that the salinity of the groundwater is not affected by seasonal seawater intrusion. Schoeller diagrams illustrate that sulphate concentrations of the groundwater of the second aquifer are relatively low compared to those of the recent seawater. Therefore, this result suggests that the brackish water in the second aquifer is probably from ancient seawater that was trapped within the sediments for a long period of time, rather than due to direct seawater intrusion.     

Groundwater recharge zones mapping using GIS: a case study in Southern part of Jordan Valley,Jordan

This paper aims at mapping the potential groundwater recharge zones in the southern part of Jordan Valley (JV). This area is considered as the most important part for agricultural production in Jordan. The methodology adopted in this study is based on utilizing the open ended SLUGGER-DQL score model, which was developed by Raymond et al (2009). Geographic information systems were used in this study to build up the different layers of this model and to create the potential groundwater recharge zones. Based on the generated SLUGGER-DQL potential map, it was found that about 70.8 % of the investigated area was categorized as high potential for groundwater recharge, 18.7 % as moderate, and 10.5 % as low potential for groundwater recharge. To validate the model results, sensitivity analysis was carried out to assess the influence of each model parameter on the obtained results. Based on this analysis, it was found that the slope parameter (S) is the most sensitive parameter among SLUGGER-DQL model parameters, followed by water level in summer (L), well density (D), water quality (Q), runoff availability (R), land use/land cover, geology (GE), whereas the lowest sensitive parameter is the geology parameter (GE). Moreover, the parameters R, D, and Q show the lowest effective weights. The effective weight for each parameter was found to differ from the assigned theoretical weight by SLUGGER-DQL index model.     

Occurrence of minor toxic elements in volcanic-sedimentary aquifers: a case study in central Italy

The occurrence of minor toxic elements (arsenic, vanadium and fluoride among others) in the volcanic aquifers of central Italy, with concentrations often above the national standards for drinking waters, has been recognized since the 1970s. However, these groundwaters are still often used for human consumption. With the aim of providing insight into the water–rock interaction processes, as well as into the evolution of arsenic and other elements, water-sampling campaigns where undertaken within a 100 km² area at the eastern margin of the Viterbo region, where volcanites overlie an alluvial aquifer complex. The samples were analysed and geochemical modelling was applied to characterize the thermodynamic state of the waters. The results indicate the existence of direct relations among arsenic, vanadium and fluoride concentrations and of their inverse relationship with calcium concentration. An evolution scheme of groundwater composition from the upstream volcanic aquifer to the downstream sedimentary aquifer is also discussed. Two hypotheses are proposed which justify the observed decrease of the cited minor element concentrations, i.e. dilution processes and/or precipitation of fluorapatite, which can include in its structure small amounts of these elements. The precipitation hypothesis is supported by modelling results.     

跟管钻进在大口径页岩气井浅部复杂地层中的应用 ——以黔水地1井导管段为例

桂中-南盘江地区的浅部地层裂隙、溶洞发育,破碎带多,井漏、井垮严重,泥浆护壁钻进成孔难度较大。针对该类地层,首次将跟管钻进技术引入大口径页岩气井钻探——黔水地1井导管段钻探,有效地解决了该井导管段的成孔难等问题。并通过系统总结该技术在黔水地1井工程施工中的实际应用情况,为今后类似工程提供了宝贵经验和建议。     

Hydraulic parameters and solute velocities in triple-porosity karstic-fissured-porous carbonate aquifers: case studies in southern Poland

 Rock and flow parameters of three karstic-fissured-porous aquifers in the Krakow-Silesian Triassic formations were measured using various methods and compared. Though cavern and fissure porosities are shown to be very low (cavern porosity below 0.5% and fracture porosity below 0.2%), they contribute dominantly to the hydraulic conductivity (from about 1.3×10^–6 to about 11×10^–6 m/s). Matrix porosity (2–11%) is shown to be the main water reservoir for solute transport and the main or significant contributor to the specific yield (<2%). Though the matrix porosity is shown to be much larger than the sum of the cavern and fissure porosities, its contribution to the total hydraulic conductivity is practically negligible (hydraulic conductivity of the matrix is from about 5×10^–11 m/s to about 2×10^–8 m/s). On the other hand, the matrix porosity (for neglected cavern and fissure porosities) when combined with tracer ages (or mean travel times) is shown to yield proper values of the hydraulic conductivity (K) by applying the following formula: K≅(matrix porosity×mean travel distance)/(mean hydraulic gradient×mean tracer age). Confirming earlier findings of the authors, this equation is shown to be of great practical importance because matrix porosity is easily measured in the laboratory on rock samples, whereas cavern and fracture porosities usually remain unmeasurable. Received: 21 February 1997 · Accepted: 13 May 1997     

Qualitative risk assessment schemes using selected parameters for hillslope developments: a case study of Penang Island

Development of hilly areas requiring deforestation and increased human activities has raised many environmental issues in the world. This paper presents the results of investigation to assess the parameters that contribute to the risk to environment by reviewing available reports from Penang Island and surrounding areas. Nine parameters were identified that mostly contribute to the environmental risk. These are casualties of people, soil properties, earth coverage, soil grading characteristics, land use suitability, factor of safety, blasting area, distance between proposed structure to slope and geometry of slope. From these parameters, the rated risk approaches were proposed based on the prediction of occurrences of environmental issues. The rated risks were categorised from very high (5) to very low (1) with high, medium and low in between. The assessments were reviewed and evaluated based on rated risk for the chosen parameters. A map was developed identifying areas of environmental risk. The map can be used to quantify risk to the environment based on the parameters due to the proposed development. The overall result shows that areas PD3 of Tanjung Bungah and PD3, PD8 and PD9 of Paya Terubong have high-rated risk. No development area is classified as very high risk and very low-risk areas. These rating approaches can be used to develop potential risk zoning of proposed development areas. From this zoning guideline or standard, the viability of any proposed development can be measured.     

Management of saltwater intrusion in coastal aquifers using different wells systems: a case study of the Nile Delta aquifer in Egypt

Hydrogeology Journal - Saltwater intrusion (SWI) is a type of pollution that adversely affects the quality of groundwater in coastal aquifers. The Nile Delta aquifer (NDA) in Egypt contains a large...     

Groundwater recharge estimation using chloride,stable isotopes and tritium profiles in the sands of northwestern Senegal

A study of environmental chloride, deuterium, oxygen-18, and tritium in deep sand profiles (35 m) has been carried out in order to estimate their relative value for measuring average groundwater recharge. The investigation was located at a 0.1-km² site in Quaternary sands near the northwestern coast of Senegal in a zone of rainfed agriculture. By using a steady-state model for duplicate unsaturated zone chloride profiles, the long-term average recharge at the site was estimated to be 30 mm yr^–1 or around 10% of the average precipitation (290 mm). The chloride concentration of adjacent shallow groundwater was relatively uniform and comparable to the unsaturated zone average, while the spatial variability in the depth distribution of Cl^– in the unsaturated zone was considerable. Stable isotope (deuterium and oxygen-18) data show that there is some isotopic enrichment due to direct evaporation through the soil surface. The degree of heavy isotope enrichment is proportional to the extent of evaporative loss and there is good correspondance with the chloride enrichment. Nevertheless, stable isotopes cannot be used quantitatively to estimate the recharge. The excellent preservation of the peak in thermonuclear tritium in precipitation in the unsaturated zone at depths between 12 and 20 m enables an estimated annual recharge of 24 mm yr^–1 in this area to be calculated, using the piston flow model. Agreement therefore between Cl and³H as tools for recharge measurement is reasonable over the site.     

Arsenic mobilization in shallow aquifers of Datong Basin: Hydrochemical and mineralogical evidences

To better understand the sources and mobilization processes responsible for arsenic enrichment in groundwater in the central part of Datong Basin where serious arsenic poisoning cases have been reported, hydrochemical characteristics of the groundwater and the geochemical and mineralogical features of the aquifer sediments were studied. The aqueous arsenic levels are strongly depth-dependent in the study area and the high arsenic concentrations are found at depths between 15 m and 60 m, with a maximum up to 1820 μg/L. The hydrochemical characteristics of high arsenic groundwater from the Datong Basin indicate that the mobilization of arsenic is related to reductive dissolution of Fe oxides/oxyhydroxides and/or desorption from the Fe oxides/oxyhydroxides at high pH (above 8.0). The bulk chemical results of sediments show the arsenic and iron are moderately correlated, suggesting that arsenic is associated with iron-bearing minerals. Results of sequential-extraction experiment show that solid-phase arsenic is similarly distributed among the different pools of reservoir in the aquifer sediments. Strongly adsorbed arsenic and co-precipitated arsenic are its dominant species in the solid-phase. Geochemical studies using chemical analysis, X-ray diffraction and scanning electron microscopy on magnetically separated fractions demonstrate that iron oxides/oxyhydroxides with residual magnetite and chlorite, illite, iron oxides/oxyhydroxides-coated quartz and feldspar, and ankerite are the dominant carriers of arsenic in the sediments. The major processes of arsenic mobilization are probably linked to desorption of As from Fe oxides/oxyhydroxides and reductive dissolution of Fe-rich phases in the aquifer sediments under reducing and alkaline conditions.