Numerical Feasibility Study for Treated Wastewater Recharge as a Tool to Impede Saltwater Intrusion in the Coastal Aquifer of Gaza – Palestine

The ongoing depletion of the coastal aquifer in the Gaza strip due to groundwater overexploitation has led to the process of seawater intrusion, which is continually becoming a serious problem in Gaza, as the seawater has further invaded into many sections along the coastal shoreline. As a first step to get a hold on the problem, the artificial neural network (ANN)-model has been applied as a new approach and an attractive tool to study and predict groundwater levels without applying physically based hydrologic parameters, and also for the purpose to improve the understanding of complex groundwater systems and which is able to show the effects of hydrologic, meteorological and anthropogenic impacts on the groundwater conditions. Prediction of the future behaviour of the seawater intrusion process in the Gaza aquifer is thus of crucial importance to safeguard the already scarce groundwater resources in the region. In this study the coupled three-dimensional groundwater flow and density-dependent solute transport model SEAWAT, as implemented in Visual MODFLOW, is applied to the Gaza coastal aquifer system to simulate the location and the dynamics of the saltwater–freshwater interface in the aquifer in the time period 2000-2010. A very good agreement between simulated and observed TDS salinities with a correlation coefficient of 0.902 and 0.883 for both steady-state and transient calibration is obtained. After successful calibration of the solute transport model, simulation of future management scenarios for the Gaza aquifer have been carried out, in order to get a more comprehensive view of the effects of the artificial recharge planned in the Gaza strip for some time on forestall, or even to remedy, the presently existing adverse aquifer conditions, namely, low groundwater heads and high salinity by the end of the target simulation period, year 2040. To that avail, numerous management scenarios schemes are examined to maintain the ground water system and to control the salinity distributions within the target period 2011-2040. In the first, pessimistic scenario, it is assumed that pumping from the aquifer continues to increase in the near future to meet the rising water demand, and that there is not further recharge to the aquifer than what is provided by natural precipitation. The second, optimistic scenario assumes that treated surficial wastewater can be used as a source of additional artificial recharge to the aquifer which, in principle, should not only lead to an increased sustainable yield of the latter, but could, in the best of all cases, revert even some of the adverse present-day conditions in the aquifer, i.e., seawater intrusion. This scenario has been done with three different cases which differ by the locations and the extensions of the injection-fields for the treated wastewater. The results obtained with the first (do-nothing) scenario indicate that there will be ongoing negative impacts on the aquifer, such as a higher propensity for strong seawater intrusion into the Gaza aquifer. This scenario illustrates that, compared with 2010 situation of the baseline model, at the end of simulation period, year 2040, the amount of saltwater intrusion into the coastal aquifer will be increased by about 35 %, whereas the salinity will be increased by 34 %. In contrast, all three cases of the second (artificial recharge) scenario group can partly revert the present seawater intrusion. From the water budget point of view, compared with the first (do nothing) scenario, for year 2040, the water added to the aquifer by artificial recharge will reduces the amount of water entering the aquifer by seawater intrusion by 81, 77and 72 %, for the three recharge cases, respectively. Meanwhile, the salinity in the Gaza aquifer will be decreased by 15, 32 and 26% for the three cases, respectively.

Analysis of water-level response to rainfall and implications for recharge pathways in the Chalk aquifer, SE England

Water table response to rainfall was investigated at six sites in the Upper, Middle and Lower Chalk of southern England. Daily time series of rainfall and borehole water level were cross-corretated to investigate seasonal variations in groundwater-level response times, based on periods of 3-month duration. The time tags (in days) yielding significant correlations were compared with the average unsaturated zone thickness during each 3-month period. In general, for cases when the unsaturated zone was greater than 18 m thick, the time tag for a significant water-level response increased rapidly once the depth to the water table exceeded a critical value, which varied from site to site. For shallower water tables, a linear relationship between the depth to the water table and the water-level response time was evident. The observed variations in response time can only be partially accounted for using a diffusive model for propagation through the unsaturated matrix, suggesting that some fissure flow was occurring. The majority of rapid responses were observed during the winter/spring recharge period, when the unsaturated zone is thinnest and the unsaturated zone moisture content is highest, and were more likely to occur when the rainfall intensity exceeded 5 mm/day. At some sites, a very rapid response within 24 h of rainfall was observed in addition to the longer term responses even when the unsaturated zone was up to 64 m thick. This response was generally associated with the autumn period. The results of the cross-correlation analysis provide statistical support for the presence of fissure flow and for the contribution of multiple pathways through the unsaturated zone to groundwater recharge. (c) 2006 Elsevier B.V. All rights reserved.

Multivariate relationships between groundwater chemistry and toxicity in an urban aquifer

Multivariate statistical methods were used to investigate file Causes of toxicity and controls on groundwater chemistry from 274 boreholes in an Urban area (London) of the United Kingdom. The groundwater was alkaline to neutral, and chemistry was dominated by calcium, sodium, and Sulfate. Contaminants included fuels, solvents, and organic compounds derived from landfill material. The presence of organic material in the aquifer caused decreases in dissolved oxygen, sulfate and nitrate concentrations. and increases in ferrous iron and ammoniacal nitrogen concentrations. Pearson correlations between toxicity results and the concentration of individual analytes indicated that concentrations of ammoinacal nitrogen, dissolved oxygen, ferrous iron, and hydrocarbons were important where present. However, principal component and regression analysis suggested no significant correlation between toxicity and chemistry over the whole area. Multidimensional Scaling was used to investigate differences in sites caused by historical use, landfill gas status, or position within the sample area. Significant differences were observed between sites with different historical land use and those with different gas status. Examination of the principal component matrix revealed that these differences are related to changes in the importance of reduced chemical species.

Geophysical characterization of the complex dynamics of groundwater and seawater exchange in a highly stressed aquifer system linked to a coastal lagoon (SE Spain)

Anthropogenic pressure influences the two-way interactions between shallow aquifers and coastal lagoons. Aquifer overexploitation may lead to seawater intrusion, and aquifer recharge from rainfall plus irrigation may, in turn, increase the groundwater discharge into the lagoon. We analyse the evolution, since the 1950s up to the present, of the interactions between the Campo de Cartagena Quaternary aquifer and the Mar Menor coastal lagoon (SE Spain). This is a very heterogeneous and anisotropic detrital aquifer, where aquifer–lagoon interface has a very irregular geometry. Using electrical resistivity tomography, we clearly identified the freshwater–saltwater transition zone and detected areas affected by seawater intrusion. Severity of the intrusion was spatially variable and significantly related to the density of irrigation wells in 1950s–1960s, suggesting the role of groundwater overexploitation. We distinguish two different mechanisms by which water from the sea invades the land: (a) horizontal advance of the interface due to a wide exploitation area and (b) vertical rise (upconing) caused by local intensive pumping. In general, shallow parts of the geophysical profiles show higher electrical resistivity associated with freshwater mainly coming from irrigation return flows, with water resources mostly from deep confined aquifers and imported from Tagus river, 400 km north. This indicates a likely reversal of the former seawater intrusion process.

Uranium series disequilibria in ground waters from a fractured bedrock aquifer (Morungaba Granitoids-Southern Brazil): Implications to the hydrochemical behavior of dissolved U and Ra

Activity concentrations of dissolved U-234, U-238, Ra-226 and Ra-228 were determined in ground waters fromtwo deep wells drilled in Morungaba Granitoids (Southern Brazil). Sampling was done monthly for little longer than 1 year. Significant disequilibrium between U-238, U-234 and Ra-226 were observed in all samples. The variation of U-238 and U-234 activity concentrations and U-234/U-238 activity ratios is related to seasonal changes. Although the distance between the two wells is short (about 900m), systematic differences of activity concentrations of U isotopes, as well as of U-234/U-238, Ra-226/U-234 and Ra-228/Ra-226 activity ratios were noticed, indicating distinct host rock-water interactions. Slightly acidic ground water percolation through heterogeneous host rock, associated with different recharge processes, may explain uranium and radium isotope behavior. (c) 2008 Elsevier Ltd. All rights reserved.

Geochemistry of natural chromium occurrence in a sandstone aquifer in Bauru Basin, Sao Paulo State, Brazil

Anomalous concentrations of Cr(VI) occur in groundwaters of the Adamantina Aquifer, in a large region in the western state of Sao Paulo, sometimes exceeding the potability limit (0.05 mg L(-1)). To identify the possible geochemical reactions responsible for the occurrence of Cr in groundwater in Urania, borehole rock samples were collected in order to carry out mineralogical and chemical analyses. In addition, multilevel monitoring wells were installed and groundwater samples were analyzed. Analyses of the borehole rock samples show the occurrence of a geochemical anomaly of Cr in the quartzose sandstones (average concentrations of 221 ppm). Chrome-diopside is one of the main minerals contributing to this anomaly, having an average Cr content of 1505 ppm. Sequential extraction experiments indicated weakly adsorbed Cr in the order of 0.54 ppm, and this quantity is enough to provide the Cr concentrations observed in groundwater. Groundwaters from the monitoring wells proved to be stratified, with the highest concentrations of Cr(VI) (0.13 mg L(-1)) being associated with high redox and pH values (over 10) and high concentrations of Na. Geochemical reactions that may explain the release of Cr from the solid phase to groundwater involve the release of Cr(III) from minerals (like chrome-diopside and Cr-Fe hydroxide), followed by oxidation of Cr(III) to Cr(VI), probably related to the reduction of Mn oxides present in the aquifer. Then cation exchange occurs and dissolution of carbonates which increases the pH of groundwater, resulting in the desorption and mobilization of Cr(VI) into groundwater. (C) 2011 Elsevier Ltd. All rights reserved.

Relation between sedimentary framework and hydrogeology in the Guarani Aquifer System in Sao Paulo state, Brazil

This paper presents the results of a new investigation of the Guarani Aquifer System (SAG) in Sao Paulo state. New data were acquired about sedimentary framework, flow pattern, and hydrogeochemistry. The flow direction in the north of the state is towards the southwest and not towards the west as expected previously. This is linked to the absence of SAG outcrop in the northeast of Sao Paulo state. Both the underlying Piramboia Formation and the overlying Botucatu Formation possess high porosity (18.9% and 19.5%, respectively), which was not modified significantly by diagenetic changes. Investigation of sediments confirmed a zone of chalcedony cement close to the SAG outcrop and a zone of calcite cement in the deep confined zone. The main events in the SAG post-sedimentary history were: (1) adhesion of ferrugineous coatings on grains, (2) infiltration of clays in eodiagenetic stage, (3) regeneration of coatings with formation of smectites, (4) authigenic overgrowth of quartz and K-feldspar in advanced eodiagenetic stage, (5) bitumen cementation of Piramboia Formation in mesodiagenetic stage, (6) cementation by calcite in mesodiagenetic and telodiagenetic stages in Piramboia Formation, (7) formation of secondary porosity by dissolution of unstable minerals after appearance of hydraulic gradient and penetration of the meteoric water caused by the uplift of the Serra do Mar coastal range in the Late Cretaceous, (8) authigenesis of kaolinite and amorphous silica in unconfined zone of the SAG and cation exchange coupled with the dissolution of calcite at the transition between unconfined and confined zone, and (9) authigenesis of analcime in the confined SAG zone. The last two processes are still under operation. The deep zone of the SAG comprises an alkaline pH, Na-HCO(3) groundwater type with old water and enriched delta(13)C values (<-3.9), which evolved from a neutral pH, Ca-HCO(3) groundwater type with young water and depleted delta(13)C values (>-18.8) close to the SAG outcrop. This is consistent with a conceptual geochemical model of the SAG, suggesting dissolution of calcite driven by cation exchange, which occurs at a relatively narrow front recently moving downgradient at much slower rate compared to groundwater flow. More depleted values of delta(18)O in the deep confined zone close to the Parana River compared to values of relative recent recharged water indicate recharge occur during a period of cold climate. The SAG is a ""storage-dominated"" type of aquifer which has to be managed properly to avoid its overexploitation. (C) 2011 Elsevier Ltd. All rights reserved.

Anomalous content of chromium in a Cretaceous sandstone aquifer of the Bauru Basin, state of Sao Paulo, Brazil

Anomalous and natural concentrations of Cr(6+), occasionally exceeding the permitted limit for human consumption (0.05 mg/L), have been detected in groundwater in the northwestern region of the state of Sao Paulo. As part of a water-rock interaction investigation, this article describes the chemical and mineralogical characterization of rock samples taken from boreholes in the municipality of Urania, with the objective of identifying Cr-bearing minerals and determining how chromium is associated with these minerals. Rock sample analysis were performed using X-ray Fluorescence, X-ray Diffraction, Scanning Electron Microscopy, electron microprobe and sequential extraction techniques. Chemical analyses indicated that the quartzose sandstones show a geochemical anomaly of chromium, with an average content of 221 ppm, which is higher than the reported chromium content of generic sandstones (35 ppm). Diopside was identified as the primary Cr-bearing mineral potentially subject to weathering processes, with a chromium content of up to 1.2% as Cr(2)O(3). Many of the diopside grains showed dissolution features, confirming the occurrence of weathering. Sequential extraction experiments indicated that 99.3% of the chromium in samples is tightly bonded to minerals, whereas 0.24% is weakly bonded via adsorption. Assuming hypothetically that all adsorbed chromium is released via desorption, the theoretical Cr concentration in water would be one order of magnitude higher than the concentrations of Cr(6+) detected in groundwater. (C) 2010 Elsevier Ltd. All rights reserved.

Stochastic simulation of time-series models combined with geostatistics to predict water-table scenarios in a Guarani Aquifer System outcrop area, Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Using neural networks for estimation of aquifer dynamical behavior

The systems of water distribution from groundwater wells can be monitored using the changes observed on its dynamical behavior. In this paper, artificial neural networks are used to estimate the depth of the dynamical water level of groundwater wells in relation to water flow, operation time and rest time. Simulation results are presented to demonstrate the validity of the proposed approach. These results have shown that artificial neural networks can be effectively used for the identification and estimation of parameters related to systems of water distribution.

Stable isotopes (H-2, O-18 and C-13) in groundwaters from the northwestern portion of the Guarani Aquifer System (Brazil)

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

A combined method for evaluating radon and progeny in waters and its use at Guarani aquifer, São Paulo State, Brazil

A combined method for evaluating radon (Rn-222) and progeny (Pb-214 and Bi-214) in water was developed by using inexpensive alpha scintillation Counting and gamma ray spectrometry through NaI(Tl) scintillation detectors. A groundwater sample collected at the Pocos de Caldas alkaline massif in Brazil was submitted to the technique in order to assure its applicability by comparing the volumetric activities by different methods. Similar volumetric activity was determined for Pb-214 and Bi-214 in the sample analyzed that is compatible with the expected condition of radioactive equilibrium between these nuclides. The combined method was successfully used to analyze groundwater samples from Guarani aquifer in S (a) over tildeo Paulo State, Brazil, and the results of the measurements indicated that Pb-214 and Bi-214 provide useful information concerning the evaluation of the drinking water quality in terms of radiological aspects. This is because they are directly identified in the water samples, without the need of requiring the assumption of the establishment of the transient equilibrium condition with its parent Rn-222. (c) 2005 Elsevier Ltd. All rights reserved.

Hydro(radio)chemical relationships in the giant Guarani aquifer, Brazil

This investigation was carried out within the Parana sedimentary basin and involved the sampling of 78 pumped tubular wells for evaluating the hydrochemistry and radioactivity due to the nuclides (238)U, (234)U, (222)Ra, (226)Ra, and (228)Ra in the Brazilian part of Guarani aquifer. Several significant correlations were found involving the geostatic pressure, for instance, specific flow rate, CO(3)(2-), SO(4)(2-) temperature, dissolved O(2), free CO(2), pH, redox potential Eh, conductivity, Na, HCO(3)-, CO(3)(2-) , SI(calcite), Cl(-), F(-), SO(4)(2-), and B. Carbonates precipitation was evidenced by inverse correlation between CO(3)(2-) and Ca, Mg, Sr, and Ba, whereas Na exhibited an opposite trend, dissolving rather than precipitating with increasing CO(3)(2-) concentration. An inverse correlation between 3 and K was found, possibly related to the increasing tendency of K to recombine with the thickness of the clayey layers. HCO(3)-played an important role on Na, Ca, Mg, and Sr dissolution. The dissolved U content and (234)U/(238)U activity ratio data were plotted on a two-dimensional diagram that was successfully utilized on identifying an unreported zone of U accumulation, though not necessarily of economic size and grade. The variability in chemical and radionuclides data indicated an important influence of the underlying Paleozoic sediments in the composition of waters from Guarani aquifer. The available data allowed estimate the groundwater residence time by two U-isotopes disequilibrium methods. Values of 45-61 ka were initially calculated, depending on the adopted porosity (15-20%), but a longer residence time (- 640 ka) was also estimated, which is more compatible with the hydraulic conductivity data in Guarani aquifer and groundwater flow velocity occurring at Milk River aquifer, Alberta, Canada. Such time range agrees with previously reported (14)C ages exceeding 30 ka BP at the more central parts of the Parana sedimentary basin. (c) 2005 Elsevier B.V. All rights reserved.

Some possible evolutionary scenarios suggested by Cl-36 measurements in Guarani aquifer groundwaters

The Guarani aquifer underlies 1.2 M km(2) in the Parana sedimentary basin of South America and is an important source of water for industry, agriculture, and domestic supplies. To determine the sustainability of this aquifer we need to understand the dynamics of the groundwater system. This paper describes the first Cl-36 measurements on aquifer groundwaters and some measurements on South American rainwaters, thought to be indicative of the recharge water. The results are compared to previous work in the region, including other radioisotope analyses. A simple model is developed, incorporating radioactive decay, allowing scenarios to be developed for mixing different waters at different mixing rates. Thus, mixing scenarios consistent with other hydrogeological and hydrogeochemical data could be assessed. A model that mixes fresh recharging waters with formational waters, that contain elevated chloride levels, but low (in situ) Cl-36 levels, can explain most of the results presented here. The expectation that rainwater samples would provide a good end-member for modelling recharge proved problematic, however. As a consequence, it is suggested that either: the recharge waters are not sourced from the same locations as the rains; that the current rainfall and fallout conditions were significantly different in the past; or that the low levels of chloride in rainfall may have allowed some contamination of the samples by old (Cl-36-free) chloride during the recharge process. (c) 2008 Elsevier Ltd. All rights reserved.

Dissolved Po-210 and Pb-210 in Guarani aquifer groundwater, Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)