Hydrogeologic assessment of escalating groundwater exploitation in the Indus Basin, Pakistan

Groundwater development has contributed significantly to food security and reduction in poverty in Pakistan. Due to rapid population growth there has been a dramatic increase in the intensity of groundwater exploitation leading to declining water tables and deteriorating groundwater quality. In such prevailing conditions, the hydrogeological appraisal of escalating groundwater exploitation has become of paramount importance. Keeping this in view, a surface water–groundwater quantity and quality model was developed to assess future groundwater trends in the Rechna Doab (RD), a sub-catchment of the Indus River Basin. Scenario analysis shows that if dry conditions persist, there will be an overall decline in groundwater levels of around 10 m for the whole of RD during the next 25 years. The lower parts of RD with limited surface water supplies will undergo the highest decline in groundwater levels (10 to 20 m), which will make groundwater pumping very expensive for farmers. There is a high risk of groundwater salinization due to vertical upconing and lateral movement of highly saline groundwater into the fresh shallow aquifers in the upper parts of RD. If groundwater pumping is allowed to increase at the current rate, there will be an overall decline in groundwater salinity for the lower and middle parts of RD because of enhanced river leakage.     

Study of geothermal water intrusion due to groundwater exploitation in the Puebla Valley aquifer system, Mexico

Significant intrusion of geothermal water into fresh groundwater takes place in the Puebla Valley aquifer system, Mexico. The decline in the potentiometric surface due to the overexploitation of the groundwater induces this intrusion. This hydrological system comprises three aquifers located in Plio-Quaternary volcanic sediments and Mesozoic calcareous rocks. The hydraulic balance of the aquifer shows that the annual output exceeds the natural inputs by 12 million m³. Between 1973 and 2002, a drop in the potentiometric surface, with an 80 m cone of depression, was identified in a 5-km-wide area located southwest of the city of Puebla. Chemical analyses performed on water samples since 1990 have shown an increase in total dissolved solids (TDS) of more than 500 mg/L, coinciding with the region showing a cone of depression in the potentiometric surface. A three-dimensional flow and transport model, based on the hydrogeological and geophysical studies, was computed by using the MODFLOW and MT3D software. This model reproduces the evolution of the aquifer system during the last 30 years and predicts for 2010 an additional drawdown in the potentiometric surface of 15 m, and an increase in the geothermal water intrusion.     

Quantifying Organization of Atmospheric Turbulent Eddy Motion Using Nonlinear Time Series Analysis

Using three methods from nonlinear dynamics, we contrast the level of organization inthe vertical wind velocity (w) time series collected in the atmospheric surface layer(ASL) and the canopy sublayer (CSL) for a wide range of atmospheric stability ()conditions. The nonlinear methods applied include a modified Shannon entropy, waveletthresholding, and mutual information content. Time series measurements collected overa pine forest, a hardwood forest, a grass-covered forest clearing, and a bare soil, desertsurface were used for this purpose. The results from applying all three nonlinear timeseries measures suggest that w in the CSL is more organized than that in the ASL, and that as the flows in both layers evolve from near-neutral to near-convective conditions, the level of organization increases. Furthermore, we found that the degree of organization in w associated with changes in is more significant than the transition from CSL to ASL.     

Foliation fields and 3D cartography in geology: Principles of a method based on potential interpolation

A modeling method that takes into account known points on a geological interface and plane orientation data such as stratification or foliation planes is described and tested. The orientations data do not necessarily belong to one of the interfaces but are assumed to sample the main anisotropy of a geological formation as in current geological situations. The problem is to determine the surfaces which pass through the known points on interfaces and which are compatible with the orientation data. The method is based on the interpolation of a scalar field defined in the space the gradient in which is orthogonal to the orientations, given that some points have the same but unknown scalar value (points of the same interface), and that scalar gradient is known on the other points (foliations). The modeled interfaces are represented as isovalues of the interpolated field. Preliminary two-dimensional tests carried-out with different covariance models demonstrate the validity of the method, which is easily transposable in three dimensions.     

Integrated resource management

Food, energy and mineral resources appeared to be available in limitless quantities. The accelerated growth in global population during recent centuries brought about a new situation. We reached the resource limits of Spaceship Earth. This, along with the dawning awareness of the extent of environmental damage to date, have made it necessary to turn to new attitudes as to the use of renewable and non-renewable resources.IRM (Integrated Resource Management) is proposed as a way to rationally balance the use and conservation of natural resources. A databank based on the Geographic Information System (GIS) concept, bringing together all relevant resources of a nation, is the centerpiece of IRM and will make it possible to establish strategies to optimize the long term use of resources, including all relevant environmental factors.     

http://www.sciencedirect.com/science/article/pii/S1674987111000405

This study has focused on the processes of soil degradation and chemical element concentration in tea-growing regions of Rwanda, Africa. Soil degradation accelerated by erosion is caused not only by topography but also by human activities. This soil degradation involves both the physical loss and reduction in the amount of topsoil associated with nutrient decline. Soil samples were collected from eleven tropical zones in Rwanda and from variable depth within each collecting site. Of these, Samples from three locations in each zone were analyzed in the laboratory, with the result that the pH of all soil samples is shown to be less than 5 (pH < 5) with a general average of 4.4. The elements such as iron (Fe), copper (Cu), manganese (Mn), and zinc (Zn) are present in high concentration levels. In contrast calcium (Ca) and sodium (Na) are present at low-level concentrations and carbon (C) was found in minimal concentrations. In addition, elements derived from fertilizers, such as nitrogen (N), phosphorous (P), and potassium (K) which is also from minerals such as feldspar, are also present in low-level concentrations. The results indicate that the soil in certain Rwandan tea plantations is acidic and that this level of pH may help explain, in addition to natural factors, the deficiency of some elements such as Ca, Mg, P and N. The use of chemical fertilizers, land use system and the location of fields relative to household plots are also considered to help explain why tea plantation soils are typically degraded.     

Modélisation géométrique 3D des granites stéphaniens du massif du Pelvoux (Alpes,France)

The structural analysis and the 3D modelling of Stephanian granites of the Pelvoux Massif characterize an emplacement along sinistral NW–SE- and dextral NE–SW-trending shear zones in the Pelvoux and in the Aiguilles Rouges–Mont Blanc Massifs, respectively. This Carboniferous shear system is consistent with a north–south extension direction known in the whole Variscan belt at this time. To cite this article: P. Strzerzynski et al., C. R. Geoscience 337 (2005).     

Exponential time integrators for stochastic partial differential equations in 3D reservoir simulation

The transport of chemically reactive solutes (e.g. surfactants, CO₂ or dissolved minerals) is of fundamental importance to a wide range of applications in oil and gas reservoirs such as enhanced oil recovery and mineral scale formation. In this work, we investigate exponential time integrators, in conjunction with an upwind weighted finite volume discretisation in space, for the efficient and accurate simulation of advection–dispersion processes including non-linear chemical reactions in highly heterogeneous 3D oil reservoirs. We model sub-grid fluctuations in transport velocities and uncertainty in the reaction term by writing the advection–dispersion–reaction equation as a stochastic partial differential equation with multiplicative noise. The exponential integrators are based on the variation of constants solution and solve the linear system exactly. While this is at the expense of computing the exponential of the stiff matrix representing the finite volume discretisation, the use of real Léja point or the Krylov subspace technique to approximate the exponential makes these methods competitive compared to standard finite difference-based time integrators. For the deterministic system, we investigate two exponential time integrators, the second-order accurate exponential Euler midpoint (EEM) scheme and exponential time differencing of order one (ETD1). All our numerical examples demonstrate that our methods can compete in terms of efficiency and accuracy compared with standard first-order semi-implicit time integrators when solving (stochastic) partial differential equations that model mixing and chemical reactions in 3D heterogeneous porous media. Our results suggest that exponential time integrators such as the ETD1 and EEM schemes could be applied to typical 3D reservoir models comprising tens to hundreds of thousands unknowns.     

Generation of 3D Spatially Variable Anisotropy for Groundwater Flow Simulations

Sedimentary units generally present anisotropy in their hydraulic properties, with higher hydraulic conductivity along bedding planes, rather than perpendicular to them. This common property leads to a modeling challenge if the sedimentary structure is folded. In this paper, we show that the gradient of the geological potential used by implicit geological modeling techniques can be used to compute full hydraulic conductivity tensors varying in space according to the geological orientation. For that purpose, the gradient of the potential, a vector normal to the bedding, is used to construct a rotation matrix that allows the estimation of the 3D hydraulic conductivity tensor in a single matrix operation. A synthetic 2D cross section example is used to illustrate the method and show that flow simulations performed in such a folded environment are highly influenced by this rotating anisotropy. When using the proposed method, the streamlines follow very closely the folded formation. This is not the case with an isotropic model.