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.

Terrain characteristics influencing aquifer salinisation of the Glenelg-Hopkins CMA region, using GIS

Salinisation of aquifers is an issue of great concern in the Glenelg-Hopkins region. The GlenelgHopkins region is located in south-west Victoria, south of the Great Dividing Range and covers 2.6 million hectares. The area receives an annual average rainfall of 500-910 mm and experiences a Mediterranean climate, with hot, dry summers and cold wet winters and has varied geology and soil types. Terrain characteristics, such as soil type, geology, depth-to-water table, land use and topography have been integrated into a Geographic Information System (GIS). A geostatistical approach, including the use of multiple linear regression is used to analyse the spatial variability of the relationships between aquifer salinity and terrain characteristics across the entire region. Results from this study should greatly improve knowledge of aquifer salinisation across the region. It is expected that this work will enable managers to determine the most appropriate mitigating measures for each specific area affected.

Sudden drawdown and drainage of a horizontal aquifer

Drainage of a saturated horizontal aquifer following a sudden drawdown is reanalyzed using the Boussinesq equation. The effect of the finite length of the aquifer is considered in detail. An analytical approximation based on a superposition principle yields a very good estimate of the outflow when compared to accurate numerical solutions. An illustration of the new analytical approach to analyze basin-scale field data is used to demonstrate possible field applications of the new solution.

Controlled field studies on soil aquifer treatment in a constructed coastal sandfill

A controlled artificial recharge experiment was conducted to investigate the effect of soil aquifer treatment during percolation of secondary and tertiary (ultrafiltered) treated wastewater through the shallow vadoze zone of a newly constructed coastal sandfill. The sandfill is a reclaimed land constructed from marine sand dredged from the seabed. To obtain 1-D flow, a stainless steel column was driven to a depth of 2.5 m, penetrating the phreatic surface. Wastewater was percolated through the column under fully-saturated and unsaturated conditions. Infiltration rates, dissolved organic carbon (DOC) and ultra-violet absorption (UVA) were monitored. The wastewaters were recharged at similar infiltration rates of approximately 5.5 m/day and 3.5 m/day under fully-saturated and unsaturated conditions, respectively. In both cases, clogging occurred 40 days after the start of recharge, under saturated conditions. For secondary treated wastewater, DOC concentration (mg/l) reduced by 28% and 13% under unsaturated and saturated conditions, respectively. The corresponding UVA reduction was 19.4% and 14.1%. Similar reductions in DOC were observed for the tertiary treated wastewater; however, the reduction in UVA was higher; 28% and 22% under unsaturated and saturated conditions, respectively. On an mass removal (mg/m(2) DOC) basis, DOC reduction appeared to be more significant under unsaturated conditions. This is attributed to the presence of interstitial oxygen.

DOC and UVA attenuation with soil aquifer treatment in the saturated zone of an artificial coastal sandfill

Results from a direct recharge experiment conducted in the field to investigate DOC and UVA(254) attenuation rates during the direct injection of UF treated wastewater into a artificial coastal sandfill are presented in this paper. Approximately 500 m(3) of ultra-filtered wastewater was injected into the saturated zone, over a period of 9 days. The movement of the plume was tracked over 80 days, during which time samples were obtained from multilevel samplers installed in transects across the drift axis of the plume. An analysis of fluorescein in the samples obtained during the drift of the UF plume showed that DOC and UVA were attenuated beyond rates predicted by conservative mixing, by up to a maximum of 45%. A degradation coefficient of 0.0175 day(-1) was found to be applicable for DOC degradation. After a drift period of 80 days, DOC and UVA reduced to approximately 4.5 mg/l and 0.100 cm(-1), respectively, from initial values of 8.06 mg/l and 0.199 cm(-1).

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)