Gulraiz Akhter

Performance of CMORPH, TMPA, and PERSIANN rainfall datasets over plain, mountainous, and glacial regions of Pakistan

The present study aims at the assessment of six satellite rainfall estimates (SREs) in Pakistan. For each assessed products, both real-time (RT) and post adjusted (Adj) versions are considered to highlight their potential benefits in the rainfall estimation at annual, monthly, and daily temporal scales. Three geomorphological climatic zones, i.e., plain, mountainous, and glacial are taken under considerations for the determination of relative potentials of these SREs over Pakistan at global and regional scales. All SREs, in general, have well captured the annual north-south rainfall decreasing patterns and rainfall amounts over the typical arid regions of the country. Regarding the zonal approach, the performance of all SREs has remained good over mountainous region comparative to arid regions. This poor performance in accurate rainfall estimation of all the six SREs over arid regions has made their use questionable in these regions. Over glacier region, all SREs have highly overestimated the rainfall. One possible cause of this overestimation may be due to the low surface temperature and radiation absorption over snow and ice cover, resulting in their misidentification with rainy clouds as daily false alarm ratio has increased from mountainous to glacial regions. Among RT products, CMORPH-RT is the most biased product. The Bias was almost removed on CMORPH-Adj thanks to the gauge adjustment. On a general way, all Adj versions outperformed their respective RT versions at all considered temporal scales and have confirmed the positive effects of gauge adjustment. CMORPH-Adj and TMPA-Adj have shown the best agreement with in situ data in terms of Bias, RMSE, and CC over the entire study area.

Hydrostratigraphy and hydrogeology of the western part of Maira area, Khyber Pakhtunkhwa, Pakistan: a case study by using electrical resistivity

Hydrostratigraphy and hydrogeology of the Maira vicinity is important for the characterization of aquifer system and developing numerical groundwater flow models to predict the future availability of the water resource. Conventionally, the aquifer parameters are obtained by the analysis of pumping tests data which provide limited spatial information and turn out to be costly and time consuming. Vertical electrical soundings and pump testing of boreholes were conducted to delineate the aquifer system at the western part of the Maira area, Khyber Pakhtun Khwa, Pakistan. Aquifer lithology in the eastern part of the study area is dominated by coarse sand and gravel whereas the western part is characterized by fine sand. An attempt has been made to estimate the hydraulic conductivity of the aquifer system by establishing a relationship between the pumping test results and vertical electrical soundings by using regression technique. The relationship is applied to the area along the resistivity profiles where boreholes are not drilled. Our findings show a good match between pumped hydraulic conductivity and estimated hydraulic conductivity. In case of sparse borehole data, regression technique is useful in estimating hydraulic properties for aquifers with varying lithology.

Determining the depositional pattern by resistivity–seismic inversion for the aquifer system of Maira area, Pakistan

Velocity and density measured in a well are crucial for synthetic seismic generation which is, in turn, a key to interpreting real seismic amplitude in terms of lithology, porosity and fluid content. Investigations made in the water wells usually consist of spontaneous potential, resistivity long and short normal, point resistivity and gamma ray logs. The sonic logs are not available because these are usually run in the wells drilled for hydrocarbons. To generate the synthetic seismograms, sonic and density logs are required, which are useful to precisely mark the lithology contacts and formation tops. An attempt has been made to interpret the subsurface soil of the aquifer system by means of resistivity to seismic inversion. For this purpose, resistivity logs and surface resistivity sounding were used and the resistivity logs were converted to sonic logs whereas surface resistivity sounding data transformed into seismic curves. The converted sonic logs and the surface seismic curves were then used to generate synthetic seismograms. With the utilization of these synthetic seismograms, pseudo-seismic sections have been developed. Subsurface lithologies encountered in wells exhibit different velocities and densities. The reflection patterns were marked by using amplitude standout, character and coherence. These pseudo-seismic sections were later tied to well synthetics and lithologs. In this way, a lithology section was created for the alluvial fill. The cross-section suggested that the eastern portion of the studied area mainly consisted of sandy fill and the western portion constituted clayey part. This can be attributed to the depositional environment by the Indus and the Kabul Rivers.

OIL-Output input language for data connectivity between geoscientific software applications

Geoscientific computing has become so complex that no single software application can perform all the processing steps required to get the desired results. Thus for a given set of analyses, several specialized software applications are required, which must be interconnected for electronic flow of data. In this network of applications the outputs of one application become inputs of other applications. Each of these applications usually involve more than one data type and may have their own data formats, making them incompatible with other applications in terms of data connectivity. Consequently several data format conversion utilities are developed in-house to provide data connectivity between applications. Practically there is no end to this problem as each time a new application is added to the system, a set of new data conversion utilities need to be developed. This paper presents a flexible data format engine, programmable through a platform independent, interpreted language named; Output Input Language (OIL). Its unique architecture allows input and output formats to be defined independent of each other by two separate programs. Thus read and write for each format is coded only once and data connectivity link between two formats is established by a combination of their read and write programs. This results in fewer programs with no redundancy and maximum reuse, enabling rapid application development and easy maintenance of data connectivity links.

Composite use of numerical groundwater flow modeling and geoinformatics techniques for monitoring Indus Basin aquifer, Pakistan

The integration of the Geographic Information System (GIS) with groundwater modeling and satellite remote sensing capabilities has provided an efficient way of analyzing and monitoring groundwater behavior and its associated land conditions. A 3-dimensional finite element model (Feflow) has been used for regional groundwater flow modeling of Upper Chaj Doab in Indus Basin, Pakistan. The approach of using GIS techniques that partially fulfill the data requirements and define the parameters of existing hydrologic models was adopted. The numerical groundwater flow model is developed to configure the groundwater equipotential surface, hydraulic head gradient, and estimation of the groundwater budget of the aquifer. GIS is used for spatial database development, integration with a remote sensing, and numerical groundwater flow modeling capabilities. The thematic layers of soils, land use, hydrology, infrastructure, and climate were developed using GIS. The Arcview GIS software is used as additive tool to develop supportive data for numerical groundwater flow modeling and integration and presentation of image processing and modeling results. The groundwater flow model was calibrated to simulate future changes in piezometric heads from the period 2006 to 2020. Different scenarios were developed to study the impact of extreme climatic conditions (drought/flood) and variable groundwater abstraction on the regional groundwater system. The model results indicated a significant response in watertable due to external influential factors. The developed model provides an effective tool for evaluating better management options for monitoring future groundwater development in the study area.

Modelling the vulnerability of groundwater to contamination in an unconfined alluvial aquifer in Pakistan

Abstract The area of Thal Doab is located in the Indus Basin and is underlain by a thick alluvial aquifer called the Thal Doab aquifer (TDA). The TDA is undergone intense hydrological stress owing to rapid population growth and excessive groundwater use for livestock and irrigated agricultural land uses. The potential impact of these land uses on groundwater quality was assessed using a DRASTIC model in a Geographic Information System environment. Seven DRASTIC thematic maps were developed at fixed scale and then combined into a groundwater vulnerability map. The resultant vulnerability index values were grouped into four zones as low, moderate, high and very high. The study has established that 76% of the land area that is underlain by the TDA has a high to very high vulnerability to groundwater contamination mainly because of a thin soil profile, a shallow water table and the presence of soils and sediments with high hydraulic conductivity values. In addition, only 2 and 22% of the total area lie in low and moderate vulnerability zones, respectively. The outcomes of this study can be used to improve the sustainability of the groundwater resource through proper land-use management.

Integrated geoscience databanks for interactive analysis and visualization

There has been a worldwide revolution in geoscientific data availability and access. An effectively infinite and instantaneous free access to geoscientific data from the World Wide System of Geoscience Data Centers and Virtual Observatories is available. In addition, national databanks and commercially available large exploration data-sets also exist. These distributed data resources impose challenges for the future to move toward their objective integration and visualization to discover new knowledge. Such advancements can facilitate meaningful interpretations and decision-making for the benefit of society at global and local scales. This article presents the Digital Earth initiative at a national level to address multiple domains, such as effective management of natural resources, interactive planning of exploration activities and monitoring, mapping and mitigation of natural hazards. It discusses a distributed geospatial data infrastructure and its importance in geoscientific data integration for efficient and interactive data retrieval, analysis and visualization. Some examples are presented to demonstrate the advantages of integrated visualization in geoscientific analysis.

A study on the evaluation of the geoid-quasigeoid separation term over Pakistan with a solution of first and second order height terms

An attempt has been made to evaluate the geoid-quasigeoid separation term over Pakistan by using solutions of terms involving first and second order terrain heights. The first term, involving the Bouguer anomaly, has a significant value and requires being incorporated in any case for determination of the geoid from the quasigoidal solution. The results of the study show that the second term of separation, which involves the vertical gravity anomaly gradient, is significant only in areas with very high terrain elevations and reaches a maximum value of 2–3 cm. The integration radius of 18 km for the evaluation of the vertical gravity anomaly gradient was found to be adequate for the near zone contribution in the case of the vertical gravity anomaly gradient. The Earth Gravity Model EGM96 height anomaly gradient terms were evaluated to assess the magnitude of the model dependent part of the separation term. The density of the topographic masses was estimated with the linear operator of vertical gravity anomaly gradient using the complete Bouguer anomaly data with an initial arbitrary density of 2.67 g/cm3 to study the effect of variable Bouguer density on the geoid-quasigeoid separation. The density estimates seem to be reasonable except in the area of very high relief in the northern parts. The effect of variable density is significant in the value of the Bouguer anomaly-dependent geoid-quasigeoid separation and needs to be incorporated for its true applicability in the geoid-quasigeoid separation determination. The geoid height (N) was estimated from the geoid-quasigeoid separation term plus global part of height anomaly and terrain-dependant correction terms. The results were compared with the separation term computed from EGM96-derived gravity anomalies and terrain heights to estimate its magnitude and the possible amount of commission and omission effects.

Satellite Gravimetric Estimation of Groundwater Storage Variations Over Indus Basin in Pakistan

Like other agrarian countries, Pakistan is now heavily dependent on its groundwater resources to meet the irrigated agricultural water demand. Groundwater has emerged as a major source with more than 60% contribution in total water supplies. In the absence of groundwater regulation, the uneven and overexploitation of groundwater resource in Indus Basin has caused several problems of water table decline, groundwater mining, and deterioration of groundwater quality. This study evaluates the potential of Gravity Recovery and Climate Experiment Satellite (GRACE)-based estimation of changes in groundwater storage (GWS) as a cost-effective approach for groundwater monitoring and policy recommendations for sustainable water management in the Indus basin. The GRACE monthly gravity anomalies from 2003 to 2010 were analyzed as total water storage (TWS) variations. The variable infiltration capacity hydrological model-generated soil moisture and surface runoff were used for the separation of TWS into GWS anomalies. The GRACE-based GWS anomalies are found to favorably agree with trends inferred from in situ piezometric data. A general depletion trend is observed in Upper Indus Plain (UIP) where groundwater is found to be declining at a mean rate of about 13.5 mm per year in equivalent height of water during 2003-2010. A total loss of about 11.82 km3 per year fresh groundwater stock is inferred for UIP. Based on TWS variations and ground knowledge, the two southern river plains, Bari and Rechna are found to be under threat of extensive groundwater depletion. GRACE TWS data were also able to pick up signals from the large-scale flooding events observed in 2010 and 2014. These flooding events played a significant role in the replenishment of the groundwater system in Indus Basin. Our study indicates that the GRACE-based estimation of GWS changes is skillful enough to provide monthly updates on the trend of the GWS changes for resource managers and policy makers of Indus basin.

Geophysical Investigation of Fresh-Saline Water Interface: A Case Study from South Punjab, Pakistan

The importance of the study of fresh-saline water incursion cannot be over-emphasized. Borehole techniques have been widely used, but they are quite expensive, intrusive, and time consuming. The electrical resistivity method has proved very successful in groundwater assessment. This advanced technique uses the calculation of Dar-Zarrouk (D-Z) parameters, namely longitudinal unit conductance, transverse unit resistance, and longitudinal resistivity has been employed by using 50 vertical electrical sounding points to assess the groundwater and delineate the fresh-saline water interface over 1045 km2 area of Khanewal in Southern Punjab of Pakistan. The x-y plots and maps of D-Z parameters were produced to establish a decipherable vision for the occurrence and distribution of different water-bearing formations of fresh-saline water aquifers through a complicated situation of intermixing of different resistivity ranges for fresh-saline water bodies. This technique is useful to reduce the ambiguity produced by the process of equivalence and suppression which cause intermixing in differentiating fresh, brackish, and saline aquifers during interpretation. The fresh-saline water interface is correlated very well with the previous studies of water quality analysis carried out in Khanewal area. The results suggest that the D-Z parameters are useful for demarcating different aquifer zones. The behavior and pattern of D-Z parameters with respect to occurrence and distribution of different water-bearing formations were effectively identified and delineated in the study area.