Solute and isotope constraint of groundwater recharge simulation in an arid environment,Abu Dhabi Emirate,United Arab Emirates

A simple, well-constrained simulation of solute increases in a downgradient direction was used in the shallow unconfined aquifer of eastern Abu Dhabi Emirate (United Arab Emirates). The simulation indicates that the observed exponential increase in solute concentrations results from a combination of upward transport of solutes from underlying mudstones and evaporites, and groundwater losses by evaporation. Groundwater recharge and discharge flux in unconfined regional aquifers in arid regions commonly are difficult to estimate because there are few constraints on the flux of water lost or gained from the system. Total dissolved solids (TDS) and deuterium isotopes (δ²H) in groundwater are used to constrain estimated fluxes to the shallow aquifer of eastern Abu Dhabi Emirate. Vertical upward transport of solutes from underlying mudstones and evaporites accounts for solute increases along approximately the first 80 km of the simulated flow path, but a combination of upward solute transport and evaporation is necessary to explain observed solute concentrations beyond 80 km. Mobilization and transport of solutes in the unsaturated zone by recharging precipitation is not a significant factor.

---

Résumé Une simulation simple et bien contrainte de l’augmentation des solutés dans la direction de l’écoulement a été réalisée sur l’aquifère phréatique et libre de l’Est de l’émirat d’Abu Dhabi (émirats Arabes Unis). La simulation indique que l’augmentation observée de solutés, de forme exponentielle, est le produit de la combinaison du transport ascendant de solutés provenant des mudstones et des évaporites sous-jacents, et de l’évaporation des eaux souterraines. La recharge des eaux souterraines et les flux de vidange dans les aquifères régionaux libres des régions arides sont habituellement difficiles à estimer du fait du peu de contraintes agissant sur la perte ou le gain d’eau dans le système. La charge dissoute totale (TDS, en anglais) et les isotopes deutérium (δ²H) des eaux souterraines sont utilisés pour contraindre l’estimation des flux en direction de l’aquifère phréatique de l’Est de l’émirat d’Abu Dhabi. Le transport vertical ascendant des solutés à partir des mudstones et des évaporites fait parti des processus qui concentrent les solutés approximativement le long des 80 premiers kilomètres du trajet des écoulements simulés, mais une combinaison du transport ascendant des solutés et de l’évaporation est nécessaire pour expliquer les concentrations observées en soluté au delà des 80 km. La mobilisation et le transport des solutés par les précipitations efficaces dans la zone non saturée ne représentent pas des facteurs significatifs.

---

Resumen Se utilizó una simulación simple y muy restringida, de incrementos de soluto en la dirección del flujo, en el acuífero libre poco profundo del Emirato oriental de Abu Dhabi, (Emiratos árabes Unidos). La simulación indica que el aumento exponencial observado en las concentraciones del soluto, resulta de una combinación de transporte ascendente de solutos desde las evaporitas y lodolitas subyacentes, y por pérdidas del agua subterránea por evaporación. Los flujos de recarga y descarga de agua subterránea, en los acuíferos libres regionales en regiones áridas, normalmente son difíciles estimar, porque hay algunas restricciones en el flujo de agua perdido o ganado por el sistema. Los sólidos disueltos totales (TDS) y los isótopos de deuterio (δ²H) en el agua subterránea, se usan para forzar los flujos estimados hacia el acuífero poco profundo del Emirato oriental de Abu Dhabi. El transporte ascendente vertical de solutos desde las lodolitas y evaporitas subyacentes, involucra los aumentos del soluto a lo largo de aproximadamente los primeros 80 Km. de la dirección de flujo simulada, pero es necesaria una combinación de transporte ascendente del soluto y la evaporación para explicar las concentraciones del soluto observadas más allá de 80 Km. La movilización y transporte de solutos en la zona no saturada, por recarga debida a precipitación, no es un factor significante.

     

Archaean high-K granitoids produced by remelting of earlier Tonalite–Trondhjemite–Granodiorite (TTG) in the Sangmelima region of the Ntem complex of the Congo craton,southern Cameroon

We present a geochemical and isotopic study that, consistent with observed field relations, suggest Sangmelima late Archaean high-K granite was derived by partial melting of older Archaean TTG. The TTG formations are sodic-trondhjemitic, showing calcic and calc-alkalic trends and are metaluminous to peraluminous. High-K granites in contrast show a potassic calc-alkaline affinity that spans the calcic, calc-alkalic, alkali-calcic and alkalic compositions. The two rock groups (TTG and high-K granites) on the other hand are both ferroan and magnesian. They have a similar degree of fractionation for LREE but a different one for HREE. Nd model ages and Sr/Y ratios define Mesoarchaean and slab-mantle derived magma compositions respectively, with Nb and Ti anomalies indicating a subduction setting for the TTG. Major and trace element in addition to Sr and Nd isotopic compositions support field observations that indicate the derivation of the high-K granitic group from the partial melting of the older TTG equivalent at depth. Geochemical characteristics of the high-K granitic group are therefore inherited features from the TTG protolith and cannot be used for determining their tectonic setting. The heat budget required for TTG partial melting is ascribed to the upwelling of the mantle marked by a doleritic event of identical age as the generated high-K granite melts. The cause of this upwelling is related to linear delamination along mega-shear zones in an intracontinental setting.     

Environmental impact assessment of risk associated with groundwater overdraft remediation in cone of depression,Jining, China

In this paper, a simple scenario and probabilistic approach is used to assess the potential groundwater risk due to proposed overdraft remedial actions in cone of depression, Jining City, China. Focusing on the concentrations of Chloride ions (Cl⁻) and total hardness (TH), the impact of artificial recharge and reduced pumping on groundwater quality and quantity is analysed by using the three-dimensional finite difference groundwater flow and transport model, Visual MODFLOW, to simulate groundwater flow and transport within the study area based on scenarios, and utilizes SURFER software to map risk levels. Although 5, 10 or 15% reduced pumping with artificial recharge leads to more decrease in Cl⁻ and TH concentrations than the 25%, less volume increase is achieved for the remediation of land subsidence and other environmental problems in the cone of depression. The Cl⁻ concentrations in recovered groundwater are within the desired concentration of 200 mg/l; however, TH in some cases are above the maximum permissible limit of 500 mg/l, with an exceedence probability of about 0.67 for recharge and recharge with reduced pumping at 25%. The presence of fractures and hydrogeological complexity greatly determines impacts of remediation, and the 22% reduced pumping with artificial recharge offers an optimum strategy for overdraft remediation in the Jining cone of depression.     

Hybridization of the probability perturbation method with gradient information

Geostatistically based history-matching methods make it possible to devise history-matching strategies that will honor geologic knowledge about the reservoir. However, the performance of these methods is known to be impeded by slow convergence rates resulting from the stochastic nature of the algorithm. It is the purpose of this paper to introduce a method that integrates qualitative gradient information into the probability perturbation method to improve convergence. The potential of the proposed method is demonstrated on a synthetic history-matching example. The results indicate that inclusion of qualitative gradient information improves the performance of the probability perturbation method.     

A coupling of mixed and continuous Galerkin finite element methods for poroelasticity II: the discrete-in-time case

In this paper, we formulate a finite element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. Here, we approximate the pressure by a mixed finite element method and the displacements by a Galerkin method. Theoretical convergence error estimates are derived in a discrete-in-time setting. Of particular interest is the case when the lowest-order Raviart–Thomas approximating space or cell-centered finite differences are used in the mixed formulation and continuous piecewise linear approximations are used for displacements. This approach appears to be the one most frequently applied to existing reservoir engineering simulators.     

Estimating groundwater evapotranspiration rates using diurnal water-table fluctuations in a semi-arid riparian zone

In semi-arid climates, phreatophytes draw on shallow aquifers, and groundwater evapotranspiration (ET_G) is a principal component of groundwater budgets. Diurnal water table fluctuations, which often are a product of ET_G, were monitored in the riparian zone of Red Canyon Creek, Wyoming, USA. These fluctuations were higher in a riparian wetland (2–36 mm) than a grass-covered meadow (1–6 mm). The onset and cessation of water-table fluctuations correspond to daily temperatures relative to freezing. Spatial differences were due to vegetation type and specific yield, while temporal changes were due to vegetation dormancy. Ratios of ET_G to potential evapotranspiration (PET), K _c,GW, were similar to ratios of actual evapotranspiration (ET) to PET, K _c, in semi-arid rangelands. Before vegetation senescence, K _c,GW increased between precipitation events, suggesting phreatophytes pull more water from the saturated zone as soil moisture decreases. In contrast, K _c decreases with soil moisture following precipitation events as ET becomes increasingly water-limited. Error in ET_G is primarily from estimates of specific yield (S _y), which is difficult to quantify in heterogeneous sediments. ET_G values may be more reliable because the range of acceptable S _y is smaller than K _c and S _y does not change with vegetation type or soil moisture.     

Recent plumbing system of the Krakatau volcano revealed by teleseismic earthquake distribution

Spatial and temporal analysis of global seismological data 1964–2005 reveals a distinct teleseismic earthquake activity producing a columnar-like formation in the continental wedge between the Krakatau volcano at the surface and the subducting slab of the Indo-Australian plate. These earthquakes occur continuously in time, are in the body-wave (m _b) magnitude range 4.5–5.3 and in the depth range 1–100 km. The Krakatau earthquake cluster is vertical and elongated in the azimuth N30°E, suggesting existence of a deep-rooted fault zone cutting the Sunda Strait in the SSW-NNE direction. Possible continuation of the fault zone in the SW direction was activated by an intensive 2002/2003 aftershock sequence, elongated in the azimuth of N55°E. Beneath the Krakatau earthquake cluster, an aseismic gap exists in the Wadati-Benioff zone of the subducting plate at the depths 100–120 km. We interpret this aseismic gap as a consequence of partial melting inhibiting stress concentration necessary to generate stronger earthquakes, whereas the numerous earthquakes observed in the overlying lithospheric wedge beneath the volcano probably reflect magma ascent in the recent plumbing system of the Krakatau volcano. Focal depth of the deepest events (~100 km) of the Krakatau cluster constrains the location of the primary magma generation to greater depths. The ascending magmatic fluids stress fault segments within the Sunda Strait fault zone and change their friction parameters inducing the observed tectonic earthquakes beneath Krakatau.     

Determining the optimum cell size of digital elevation model for hydrologic application

Scale is one of the most important but unsolved issues in various scientific disciplines that deal with spatial data. The arbitrary choice of grid cell size for contour interpolated digital elevation models (DEM) is one of the major sources of uncertainty in the hydrologic modelling process. In this paper, an attempt was made to identify methods for determining an optimum cell size for a contour interpolated DEM in prior to hydrologic modelling. Twenty-meter interval contour lines were used to generate DEMs of five different resolutions, viz., 30, 45, 60, 75, and 90 m using TOPOGRID algorithm. The obtained DEMs were explored for their intrinsic quality using four different methods, i.e., sink analysis, fractal dimension of derived stream network, entropy measurement and semivariogram modelling. These methods were applied to determine the level artifacts (interpolation error) in DEM surface as well as derived stream network, spatial information content and spatial variability respectively. The results indicated that a 90 m cell size is sufficient to capture the terrain variability for subsequent hydrologic modelling in the study area. The significance of this research work is that it provides methods which DEM users can apply to select an appropriate DEM cell size in prior to detailed hydrologic modelling.     

Changes in nutrient uptake of phytoplankton under the interaction between sunlight and phosphate in the Changjiang （Yangtze） River Estuary

We conducted ship-board incubation experiments to investigate changes in nutrient uptake of phytoplankton under different phosphate concentrations and irradiances in the Changjiang River Estuary and its adjacent waters in China. Under 100% natural irradiance the uptake rates of phosphate, silicate, and nitrate were accelerated at high phosphate levels （1.84 μM）, while under low irradiance （about 50% natural irradiance） their uptake rates were restrained at the high but stimulated greatly at the intermediate phosphate concentrations （1.26 μM）, as the growth of phytoplankton, changes in nitrite and ammonium uptake didn＇t follow an obvious pattern. Our results also showed that there were linear relationships between nitrate, silicate and phosphate uptake at different phosphate concentrations under low and high irradiances, and the growth period of phytoplankton was prolonged both at the high phosphate concentrations under high irradiance and at the intermediate concentrations under low irradiance, suggesting that the limitation of phytoplankton growth mainly reflected changes in its growth period, and because no such environment （low irradiance and low phosphate concentrations） actually existed in a high turbidity zone, phytoplankton blooms hardly occurred there. In the absence of irradiance, denitrification occurred readily and phytoplankton was kept decreasing, which resulted in phosphate regeneration.