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.     

Potassium induced salinity tolerance in wheat (Triticum aestivum L.)

Water culture experiments were conducted to study the response of ten wheat genotypes to external K application (10 mmol KCI dm^?3) at seedling stage under saline condition (0 and 100 mmol NaCl dm^?3). The data showed that there was an increase in the shoot and root length with the application of external K. The increase was more pronounced under control than under saline conditions. The better performing genotypes under two treatments were Bhitai, NIAB-41, NIAB-I076 and Khirman. The enhanced growth of these genotypes under saline condition might be due to the quick response to external K application, resulting in high K/Na ratio. The results indicated that the genotypes, which have the ability of enhanced K/Na discrimination, might perform better under saline conditions when sufficient potassium is applied in the rooting medium.     

Urban development and flooding in Houston Texas, inferences from remote sensing data using neural network technique

The effects of urban development on the natural ecosystem and its link to the increased flooding in Houston, Texas were evaluated. Houston is suitable for this type of analysis due to its 1.95 million population, large geographic area and fast growth rate. Using neural network techniques, four Landsat Thematic Mapper images were grouped into five land use classes for the period 1984 to 2003: vegetation, bare ground, water, concrete and asphalt. Results show that asphalt and concrete increased 21% in the time period 1984–1994, 39% in 1994–2000 and 114%, from 2000 to 2003, while vegetation suffered an overall decrease. When change detection data are compared with runoff ratio data, a relationship between increased runoff and urban development is apparent, which indicates increased chances of flooding. Initial results of this work are made available to the public in GIS format via internet using Arc Internet Map Server (ArcIMS) at .     

Selection of hazardous waste dumpsites based on parameters effecting soil adsorption capacity – a case study

A careful selection of waste dumpsites, particularly hazardous ones, is very important for sustainable water resources management. Several laboratory experiments were carried out on the field samples to study adsorption capacity using p-dichlorobenzene (a solvent used in various industrial processes) as the test contaminant. The effect of parameters such as organic matter, clay, and iron and aluminium oxides, which are known to influence the soil adsorption capacity, are studied in the present work. Several soil samples from the Patancheru Industrial Area (Hyderabad, India) were collected and characterized. Only three soils, which had a comparatively high percentage of organic matter, clay, iron and aluminium oxide contents were used for the adsorption studies. The results clearly indicated a decrease in the adsorption capacity of the soils by as much as 75% when organic matter was removed. The other parameters such as clay and iron and aluminium oxides also play an important role in adsorption (57 and 39.8% reduction respectively). It was observed that out of the selected factors organic matter in the soils has the maximum effect regarding the adsorption of p-dichlorobenzene. Since the selected soils contain comparatively more organic matter, clay and iron and aluminium oxides in the selected industrial area, these can be used as sites for dumping hazardous waste, which can be further treated by methods like bioremediation.     

Dynamic analysis of strain localization in water‐saturated clay using a cyclic elasto‐viscoplastic model

A computational framework is presented for dynamic strain localization and deformation analyses of water‐saturated clay by using a cyclic elasto‐viscoplastic constitutive model. In the model, the nonlinear kinematic hardening rule and softening due to the structural degradation of soil particles are considered. In order to appropriately simulate the large deformation phenomenon in strain localization analysis, the dynamic finite element formulation for a two‐phase mixture is derived in the updated Lagrangian framework. The shear band development is shown through the distributions of viscoplastic shear strain, the axial strain, the mean effective stress, and the pore water pressure in a normally consolidated clay specimen. From the local stress–strain relations, more brittleness is found inside the shear bands than outside of them. The effects of partially drained conditions and mesh‐size dependency on the shear banding are also investigated. The effect of a partially drained boundary is found to be insignificant on the dynamic shear band propagation because of the rapid rate of applied loading and low permeability of the clay. Using the finer mesh results in slightly narrower shear bands; nonetheless, the results manifest convergency through the mesh refinement in terms of the overall shape of shear banding and stress–strain relations. Copyright © 2013 John Wiley & Sons, Ltd.     

Complexity measurement of regional groundwater resources system using improved Lempel-Ziv complexity algorithm

In this paper, the self-adaptive artificial fish swarm algorithm (SAAFSA) is used to optimize the coarse graining of segment numbers, which are used in the Lempel-Ziv complexity algorithm. This approach improves the Lempel-Ziv complexity (LZC) algorithm of equal probability coarse graining. As a case study, the complexities of monthly series of groundwater depth were analyzed at seven farms in the Hongxinglong Administration. GIS technology was used to create spatial distributions of monthly groundwater depths. A projection pursuit model based on the SAAFSA was established and used for complexity attribution analysis at selected farms with different degrees of complexity. The three selected farms, Hongqiling, 852, and Youyi, each represent a certain degree of complexity. Further analysis shows that precipitation, evaporation, temperature, and human activities are the primary factors that cause complexity variations in local groundwater depth. The results reveal the evolution of the complexity characteristics of local groundwater depth and provide scientific evidence for the need to effectively allocate regional water resources. Additionally, the proposed method can be applied in complexity analyses of other hydrologic features, as well as in research regarding nonlinear time series in economic, engineering, medical, and signal analyses.     

Sequence stratigraphic analysis of Eocene Rock Strata,Offshore Indus,southwest Pakistan

In this study, seismic data from two wells (Pak G2-1 and Indus Marine-1C) and age diagnostic larger benthic foraminifera (LBF) within drill cuttings has been used for the first time to identify depositional sequences within the carbonates in the Offshore Indus Basin, Pakistan. The Offshore Indus is tectonically categorized as a passive continental margin where carbonates occur as shelf carbonates in the near offshore and on volcanic seamounts in deeper waters. Seismic data analysis has indicated the presence of minor faults and carbonate buildups above the igneous basement in the south. Patterns of the seismic reflections enabled definition of three seismic facies units identified as: Unit 1 basement, represented by chaotic, moderate amplitude reflection configuration; while parallel bedding and the drape of overlying strata is typical character of Unit 2, carbonate mound facies. The younger Miocene channels represent Unit 3. The diagnosis of Alveolina vredenburgi/cucumiformis biozone confirmed the Ilerdian (55–52 Ma) stage constituting a second order cycle of deposition for the Eocene carbonates (identified as Unit 2). The carbonate succession has been mainly attributed to an early highstand system tract (HST). The environmental conditions remained favorable leading to the development of keep-up carbonates similar to pinnacle buildups as a result of aggradation during late transgressive system tract and an early HST. The carbonate sequence in the south (Pak G2-1) is thicker and fossiliferous representing inner to middle shelf depths based on fauna compared to the Indus Marine-1C in the north, which is devoid of fossils. Three biozones (SBZ 5, SBZ 6 and SBZ 8) were identified based on the occurrence of LBF. The base of the SBZ 5 zone marks the larger foraminifera turnover and the Paleocene–Eocene (P–E) boundary. The LBF encountered in this study coincides with earlier findings for the P–E boundary. Our findings indicate that the entire Ilerdian stage ranges from 55.5 to 52 Ma that was the episode of warmer water conditions on the carbonate shelves leading to the diversification of K-strategist larger foraminifera. The larger foraminiferal assemblage encountered in this study confirms the findings. The possible indication of stratigraphic-combination traps, revealed as reflection terminations, make carbonate mounds in the south a potential exploration target.     

Practical Incorporation of Multivariate Parameter Uncertainty in Geostatistical Resource Modeling

There is a need to bridge theory and practice for incorporating parameter uncertainty in geostatistical simulation modeling workflows. Simulation workflows are a standard practice in natural resource and recovery modeling, but the incorporation of multivariate parameter uncertainty into those workflows is challenging. However, the objectives can be met without considerable extra effort and programming. The sampling distributions of statistics comprise the core theoretical notion with the addition of the spatial degrees of freedom to account for the redundancy in the spatially correlated data. Prior parameter uncertainty is estimated from multivariate spatial resampling. Simulation-based transfer of prior parameter uncertainty results in posterior distributions which are updated by data conditioning and the model domain extents and configuration. The results are theoretically tractable and practical to achieve, providing realistic assessments of uncertainty by accounting for large-scale parameter uncertainty, which is often the most important component impacting a project. A simulation-based multivariate workflow demonstrates joint modeling of intrinsic shale properties and uncertainty in estimated ultimate recovery in a shale gas project. The multivariate workflow accounts for joint prior parameter uncertainty given the current well locations and results in posterior estimates on global distributions of all modeled properties. This is achieved by transferring the joint prior parameter uncertainty through conditional simulations.