Transient‐state groundwater flow in various geometries with wells at arbitrary positions: analytical methods

We have developed a method for analytically solving the porous medium flow equation in many different geometries for horizontal (two‐dimensional), homogeneous and isotropic aquifers containing impermeable boundaries and any number of pumping or injection wells located at arbitrary positions within the system. Solutions and results are presented for rectangular and circular aquifers but the method presented here is easily extendible to many geometries. Results are also presented for systems where constant head boundary conditions can be emulated internal to the aquifer boundary. Recommendations for extensions of the present work are briefly discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

Nitrogen isotope indicators of seasonal source variability to groundwater

 A nitrogen isotope study of soil water and groundwater in southern Indiana, USA, in 1991–1992 demonstrated considerable variations in nitrate degradation processes compared to an earlier investigation in 1986–1987. Although N-fertilizers were applied in May 1991, the δ¹⁵N values in soil water decreased in February 1992, indicating its delayed release into the system after substantial rainfall. The δ¹⁵N values of groundwater decreased from +12.3‰ in November 1991 to +11.3‰ in February 1992, and to +7.5‰ in March 1992. The increased residence time of nitrate in the soil resulted in increased denitrification, ammonia volatilization and plant uptake, and reduced threat to the groundwater quality. The 1986–1987 study in the area reported that excessive rainfall during the summer rapidly transported the nitrate to deeper horizons and drastically reduced volatilization and microbial reduction of nitrate, thus increasing the immediate threat to the groundwater quality in the area. The present study demonstrated that nitrogen isotopic signatures can be used to determine the effects of local soil type, rainfall, and land-use practices on the fate of nitrate in the subsurface. Received: 18 February 1997 · Accepted: 17 June 1997     

Seismic stratigraphy of the Mianwali and Bannu depressions,north-western Indus foreland basin

International Journal of Earth Sciences - Regional seismic reflection profiles, deep exploratory wells, and outcrop data have been used to study the structure and stratigraphic architecture of the...     

Imaging of Stratigraphic Pinch-Out Traps Within the Lower-Cretaceous Shaly-Sandstone System,Pakistan, Using 3D Quantitative Seismic Inverted Porosity–Velocity Modeling

Incised valleys form excellent stratigraphic pinch-out traps. Traditional seismic data analysis techniques fail to predict quantitatively the porous and low-velocity sand-fills for incised valleys. The 3D quantitative seismic inverted porosity–velocity (3DQSIPV) analysis was applied in the Indus Basin, SW Pakistan. The reflection strength attribute better portrayed the reservoir sandstone and faults compared to seismic amplitude attribute. The sweetness-based continuous wavelet transform authenticated the development of the stratigraphic play. The 17 Hz amplitude delineated the non-porous seal and porous reservoirs of sand-filled incised valley and strand plain, and faults. The integrated model of seismic attributes categorizes the reservoir and seal constituents. The petrophysical modeling corroborated the gas-bearing “sweet-spots” within the stratigraphic-based dynamical system. The facies modeling predicted the for coarse-grained sandstone and fine-grained shales, depositional environments, fluctuations of sea level and their impacts on the overall development of stratigraphic plays. The predicted density and P-wave velocity for the sandstone-filled incised valley of the lowstand system tract were?~?1.4–1.75 g/cc and?~?3217–3802 m/s, respectively. The predicted density and P-wave velocity for the sealing shales facies of strand plain of transgressive system tract were?~?1.9–2.1 g/cc and 2.55–2.7 g/cc and 3900–4700 m/s, respectively. The 3DQSIPV predicted?>?25% porosity and?~?3300 m/s velocity of reservoirs in the west. The eastern zones shows?<?12% porosity and high velocity of?~?4580 m/s. Cross-plots of porosity, velocity, and thickness showed correlation coefficients of R²?>?0.90 for inverted velocity. This workflow may serve as an analogue for the remaining oil and gas fields of the Indus Basins of Pakistan and similar geological settings of divergent plate margins.

     

Geospatial assessment of soil erosion intensity and sediment yield: a case study of Potohar Region,Pakistan

Estimation of spatial extent of soil erosion, one of the most serious forms of land degradation, is critical because soil erosion has serious implications on soil fertility, water ecosystem, crop productivity and landscape beauty. The primary objective of the current study was to assess and map the soil erosion intensity and sedimentation yield of Potohar region of Pakistan. Potohar is the rainfed region with truncated and complex topography lying at the top of the Indus Basin, the world’s largest irrigation networks of canals and barrages. Spatially explicit Revised Universal Soil Loss Equation (RUSLE) Model integrated with Remote Sensing-GIS techniques was used for detecting/mapping of erosion prone areas and quantification of soil losses. The results show that the Potohar region is highly susceptible to soil erosion with an average annual soil loss of 19 tons ha^?1 year^?1 of which the maximum erosion (70–208 tons ha^?1 year^?1) was near the river channels and hilly areas. The sediment yield due to the erosion is as high as 148 tons ha^?1 year^?1 with an average of 4.3 tons ha^?1 year^?1. It was found that 2.06% of the total area falls under severe soil erosion, 13.34% under high erosion, 15.35% under moderate soil erosion while 69.25% of the area lies in the low (tolerable) soil erosion. Chakwal and Jhelum districts of the region are seriously affected by erosion owing to their topography and soil properties. The information generated in this study is a step forward towards proper planning and implementation of strategies to control the erosion and for protection of natural resources. It is, hence, necessary that suitable water harvesting structures be made to control water to prevent soil erosion and provision of water in the lean season in this region. Tree plantation and other erosion control practices such as strip cropping can also minimize soil erosion in this region.     

ENSO and IOD analysis on the occurrence of floods in Pakistan

Different techniques have been used to discuss the existence of significant relation between the El Nino Southern Oscillation (ENSO) and Indian Ocean Dipole (IOD). Various studies present their interaction and influence on the natural disasters (i.e. drought, flood, etc.) over large parts of the globe. This study uses a Markov chain method to investigate the relation between the ENSO and IOD for the period of 62 years (1950–2011) and aggregates their influence on the occurrence of floods in Pakistan. Both data sets show similarities in the formation of transition matrices and expected number of visits from one state to another. The strong values of 2-dimensional correlation and high self-communication of the transition states confirm the existence of a possible relation between ENSO and IOD data. Moreover, significant values of dependency and stationary test endorse the applicability of the Markov chain analyses. The independent analysis shows that strong events of both data sets are co-occurred in the same flood years. During the study period maximum number of floods was observed during summer monsoon season. However, further analysis shows that after 1970, Pakistan observed the highest percentage of floods occurred per year during El Nino, Non-ENSO and positive IOD years. These observations and results demonstrate that climate variability especially ENSO and IOD should be incorporated into disaster risk analyses and policies in Pakistan.     

Peculiar motions of galaxy clusters: correlation function approach

The correlation function theory on the basis of prescribed boundary conditions provides a deeper understanding in studying the dynamical parameters of galaxy clusters. The approach approximates that the moderate dense systems discussed by a two point correlation function is helpful for describing the dynamical nature of galaxy clusters. The projected theory of two point correlation function for point mass and extended mass structures can be used an alternative tool in measuring the average peculiar motion and temperature profile of galaxy clusters.     

Gravitational clustering of galaxies in an expanding universe

We inquire the phenomena of clustering of galaxies in an expanding universe from a theoretical point of view on the basis of thermodynamics and correlation functions. The partial differential equation is developed both for the point mass and extended mass structures of a two-point correlation function by using thermodynamic equations in combination with the equation of state taking gravitational interaction between particles into consideration. The unique solution physically satisfies a set of boundary conditions for correlated systems and provides a new insight into the gravitational clustering problem.     

Long-term seasonal rainfall forecasting: efficiency of linear modelling technique

Using the lagged (past) climate indices, including El Nino–Southern Oscillation (ENSO) and Indian Ocean Dipole (IOD) as input parameters and long-term spring rainfall as outputs, calibration and validation of the linear multiple regression (MR) models have been performed. Since Australian rainfall varies both temporally and spatially, the analysis on the linear MR models was performed on regional scale. These models show the capability of linear MR technique for long-term predictions of Western Australian spring rainfall. The emphasis was given to assess the statistical correlations between Western Australian spring rainfall and dominating large-scale climate modes. The efficiency of linear modelling technique was evaluated to predict seasonal rainfall forecasting. At the same time, the Pearson correlation (R), mean absolute error, root-mean-square error and Willmott index agreement (d) were used to assess the capability of MR models. The models which fulfilled the limits of statistical significances were used for the prediction of future spring rainfall using independent data set. The results indicate that during calibration periods maximum achievable correlations varied from 0.47 to 0.53 for the selected stations. In regard to predict peaks and troughs of rainfall time series, it was found that correlations between predicted and actual peaks varied from 0.82 to 0.94 and between predicted and actual troughs varied from 0.53 to 0.91.