Deciphering heavy metal contamination zones in soils of a granitic terrain of southern India using factor analysis and GIS

Soil contamination by heavy metals has been a major concern for last few decades due to increase in urbanization and industrialization. The main objective of this research was to identify the heavy metal contaminated zones in the study area. Twenty five soil samples collected throughout the agriculture, residential and industrial areas were analysed by X-ray Fluorescence Spectrometer (XRF) for trace metals and major oxides. These metals can affect the quality of soil and infiltrate through the soil, thereby causing groundwater pollution. Based on the chemical analysis of major oxides (SiO_2, Al₂O_3, ?Fe₂O_3, MnO, MgO, CaO, Na₂O, K₂O, TiO₂, and P₂O₅) and their distribution; it is observed that these soils are predominantly siliceous type with slight enrichment of alumina component in the study area. Correlation matrix (CM) and factor analysis (FA) is employed to the heavy metal variables, viz., Ba, Cr, Cu, Ni, Pb, Rb, Sr, V, Y, Zn and Zr of the soil to determine the dominant factors contributing to the soil contamination in the area. In the analysis, five factors emerged as significant contributors to the soil quality. The total contribution of these five factors is about 90%. The contribution of the first factor is about 45% and has significant positive loadings of Co, Cr, Cu, Ni and Zn. The contribution of second factor is 22% and has significant positive loadings of Rb, Sr and Y. The contribution of third, fourth and fifth factors is 10, 8 and 5% and show positive loadings for lead, molybdenum and barium respectively to the soil contamination. The spatial variation maps deciphering different zones of heavy metal concentration in the soil were generated in a GIS (geographic information system) based environment using ArcGIS 9.3.1. The results reveal that heavy metal contamination in the area is mainly due to anthropogenic activities.     

Reservoir Quality Assessment of the Upper Permian Chhidru Formation,Salt and Surghar Ranges,Pakistan

The Late Permian succession of the Upper Indus Basin in northeastern Pakistan is represented by the carbonate-dominated Zaluch Group, which consists of the Amb, Wargal and Chhidru formations, which accumulated on the southwestern shelf of the Paleo-Tethys Ocean, north of the hydrocarbon-producing Permian strata of the Arabian Peninsula. The reservoir properties of the mixed clastic-carbonate Chhidru Formation (CFm) are evaluated based on petrography, using scanning electron microscopy (SEM), energy dispersive x-ray spectroscopy (EDX) and x-ray diffraction (XRD) techniques. The diagenetic features are recognized, ranging from marine (isopachous fibrous calcite, micrite), through meteoric (blocky calcite-I, neomorphism and dissolution) to burial (poikilotopic cement, blocky calcite-II-III, fractures, fracture-filling, and stylolites). Major porosity types include fracture and moldic, while inter- and intra-particle porosities also exist. Observed visual porosity ranges from 1.5%–7.14% with an average of 5.15%. The sandstone facies (CMF-4) has the highest average porosity of 10.7%, whereas the siliciclastic grainstone microfacies (CMF-3) shows an average porosity of 5.3%. The siliciclastic mudstone microfacies (CMF-1) and siliciclastic wacke-packestone microfacies (CMF-2) show the lowest porosities of 4.8% and 5.0%, respectively. Diagenetic processes like cementation, neomorphism, stylolitization and compaction have reduced the primary porosities; however, processes of dissolution and fracturing have produced secondary porosity. On average, the CFm in the Nammal Gorge, Salt Range shows promise and at Gula Khel Gorge, Trans-Indus, the lowest porosity.     

Facies Analysis, Sequence Stratigraphy and Diagenetic Studies of the Jurassic Carbonates of the Kohat Basin, Northwest Pakistan: Reservoir Implications

The present study deals with the depositional facies, diagenetic processes and sequence stratigraphy of the shallow marine carbonates of the Samana Suk Formation, Kohat Basin, in order to elucidate its reservoir quality. The Samana Suk Formation consists of thin to thick-bedded, oolitic, bioclastic, dolomitic and fractured limestone. Based on the integration of outcrop, petrographic and biofacies analyses, the unit is thought to have been deposited on a gentle homoclinal ramp in peritidal, lagoonal and carbonate shoal settings. Frequent variations in microfacies based sea-level curve have revealed seven Transgressive Systems Tracts (TSTs) and six Regressive Systems Tracts (RSTs). The unit has undergone various stages of diagenetic processes, including mechanical and chemical compaction, cementation, micritization, dissolution and dolomitization. The petrographic analyses show the evolution of porosity in various depositional and diagenetic phases. The fenestral porosity was mainly developed in peritidal carbonates during deposition, while the burial dissolution and diagenetic dolomitization have greatly enhanced the reservoir potential of the rock unit, as is further confirmed by the plug porosity and permeability analyses. The porosities and permeabilities were higher in shoal facies deposited in TSTs, as compared to lagoonal and peritidal facies, except for the dolomite in mudstone, deposited during RSTs. Hence good, moderate and poor reservoir potential is suggested for shoal, lagoonal and peritidal facies, respectively.     

A Periglacial Palaeoenvionment in the Upper Carboniferous–Lower Permian Tobra Formation of the Salt Range, Pakistan

The Upper Carboniferous—Lower Permian(Upper Pennsylvanian-Asselian) Tobra Formation is exposed in the Salt and Trans Indus ranges of Pakistan.The formation exhibits an alluvial plain(alluvial fan-piedmont alluvial plain) facies association in the Salt Range and Khisor Range.In addition,a stream flow facies association is restricted to the eastern Salt Range.The alluvial plain facies association is comprised of clast-supported massive conglomerate(Gmc),diamictite(Dm)facies,and massive sandstone(Sm) Hthofacies whereas the stream flow-dominated alluvial plain facies association includes fine-grained sandstone and siltstone(Fss),fining upwards pebbly sandstone(Sf),and massive mudstone(Fm) Hthofacies.The lack of glacial signatures(particularly glacial grooves and striatums) in the deposits in the Tobra Formation,which are,in contrast,present in their timeequivalent and palaeogeographically nearby strata of the Arabian peninsula,e.g.the AI Khlata Formation of Oman and Unayzah B member of the Saudi Arabia,suggests a pro-to periglacial,i.e.glaciofluvial depositional setting for the Tobra Formation.The sedimentology of the Tobra Formation attests that the Salt Range,Pakistan,occupied a palaeogeographic position just beyond the maximum glacial extent during Upper Pennsylvanian-Asselian time.     

Impact of grazing on soil,vegetation and ewe production performances in a semi-arid rangeland

Controlled grazing is considered a good management strategy to maintain or increase the live weight of livestock and to reduce vegetation degradation of rangelands. The present study investigated soil characteristics, aboveground vegetation biomass dynamics and controlled grazinginduced changes in the live weight of local ewes in the semi-arid rangeland of Ahmadun, Ziarat, Balochistan province of Pakistan. An area of 115 ha was protected from livestock grazing in April 2014. In June 2015, soil characteristics within 0-30 cm depth i.e. soil organic matter (SOM), mineral nitrogen, pH and texture in controlled and uncontrolled grazing sites were assessed. Aboveground vegetation biomass measured in early (June) and late summer (August) in 2015 and 2016. The nutritional value i.e. crude protein, phosphorus (P), neutral detergent fiber (NDF), acid detergent fiber (ADF), calcium (Ca), magnesium (Mg) and potassium (K) of dominant plant species were assessed at the beginning of experiment in 2015. Vegetation cover of controlled and uncontrolled grazing sites was also measured during the two years of the study period using the VegMeasure software. From June to November in 2015 and 2016, controlled and uncontrolled livestock grazing sites were grazed on a daily basis by local ewes with a stocking rate of 2 and 1 head ha-1 respectively. Results reveal that the organic matter contents of coarse-textured, slightly alkaline soil of the study site were in the range of 9.4 - 17.6 g kg-1 soil and showed a strong positive correlation with aboveground vegetation biomass. The biomass of plants was 56.5% and 33% greater at controlled than uncontrolled grazing site in 2015 and 2016 respectively and plant cover was also higher at controlled than uncontrolled grazing site in both years. The nutrient contents were significantly (P<0.05) lower in grasses than shrubs. In both years, the controlled grazing increased the weight gain of ewes about two folds compared to the uncontrolled grazing. The results indicate that controlled grazing improved the vegetation biomass production and small ruminant productivity.     

Application of GIS for the identification and demarcation of selective heavy metal concentrations in the urban groundwater

Groundwater is the most appropriate and widely used source of drinking water, which is increasingly threatened by pollution from industrial and agricultural activities. To check the severity of the problem, 156 groundwater samples were collected from various depths (60-110 ft) of 52 different localities in Faisalabad city, the third largest metropolis in Pakistan, and analyzed for the metals (Zn, Cu, Cd, Ni, Pb, Mn and Fe) concentration in 2009. Quantification was done by using Flame Atomic Absorption Spectrophotometer technique and the results were compared with WHO standards for drinking water quality. Results showed that the levels of Cu, Mn and Fe were below the WHO standards while the concentrations of Zn, Cd, Ni and Pb were above the recommended levels of safe drinking water. Correlation analysis among the occurrence of these heavy metals revealed a highly significant and positive correlation of Mn with Zn and Fe. A significant and positive correlation of Cd was also found with Cu and groundwater depth showing that there is strong association between Cu-Cd pair and that the Cd concentration varies with depth of groundwater in the study area. Regional patterns of heavy metals occurrence were mapped using Geographical Information System (GIS) for the identification and demarcation of risk areas. The concentration maps may be used by policymakers of the city to mitigate groundwater pollution.     

Velocity and Structural Modeling of Mesozoic Chiltan Limestone and Goru Formation for Hydrocarbon Evaluation in the Bitrisim Area, Lower Indus Basin, Pakistan

The present study focuses on building a workflow for structural interpretation and velocity modeling and implementing to Jurassic-Cretaceous succession (Chiltan Limestone and Massive sand of the Lower Goru Formation). 2D-Migrated seismic sections of the area are used as data set and in order to confirm the presence of hydrocarbons in the study area, P and S-wave seismic velocities are estimated from single-component seismic data. Some specific issues in the use of seismic data for modeling and hydrocarbon evaluation need to deal with including distinguishing the reservoir and cap rocks, and the effects of faults, folds and presence of hydrocarbons on these rocks. This study has carried out the structural interpretation and modeling of the seismic data for the identification of traps. The results demonstrate existence of appropriate structural traps in the form of horst and grabens in the area. 2D and 3D velocity modeling of the horizons indicates the presence of high velocity zones in the eastern half of the study while relatively low velocity zones are encountered in the western half of the area. Two wells were drilled in the study area (i.e. Fateh-01 and Ichhri-01) and both are dry. Immature hydrocarbons migration is considered as a failure reason for Fateh-01 and Ichhri-01 well.     

A Seismic Hazard Map Based on Geology and Shear-wave Velocity in Rawalpindi–Islamabad,Pakistan

The ‘twin cities' district of Rawalpindi–Islamabad is among the most endangered seismic regions in Pakistan,with the seismic hazard assessed(0.32 g) to intensity IX MMI for a 475-year return period. A seismic hazard map for Rawalpindi–Islamabad is presented herein, based on 85 shear-wave velocity(VS) profiles obtained through geophysical H/V measurements and from the geological map of the region. Relationships between the average top thirty-meter shear-wave velocities(VS30) and surficial geological units have been determined. The peak ground acceleration(PGA) maps for 150,475 and 2475-year return periods were converted into a seismic intensity map. Intensity increments for different soils were used to compute PGA values for 150, 475 and 2475-year return periods. Sites located on softer ground experienced a higher degree of damage from moderate earthquakes. Due to the presence of soft clay or lique?able soil and lateral spreading, a few locations may be classified as hybrid sites class C and D. This map is a critical step in facilitating code-based site classification and seismic design throughout Rawalpindi–Islamabad. Although the seismic hazard map based on seismic intensities is no longer used in engineering geology, it is still important in seismological analysis and for civil protection purposes.     

Timing of metamorphism,melting and exhumation of the Leo Pargil dome,northwest India

The Leo Pargil dome, northwest India, is a 30 km‐wide, northeast‐trending structure that is cored by gneiss and mantled by amphibolite facies metamorphic rocks that are intruded by a leucogranite injection complex. Oppositely dipping, normal‐sense shear zones that accommodated orogen‐parallel extension within a convergent orogen bound the dome. The broadly distributed Leo Pargil shear zone defines the southwest flank of the dome and separates the dome from the metasedimentary and sedimentary rocks in the hanging wall to the west and south. Thermobarometry and in‐situ U–Th–Pb monazite geochronology were conducted on metamorphic rocks from within the dome and in the hanging wall. These data were combined with U–Th–Pb monazite geochronology of leucogranites from the injection complex to evaluate the relationship between metamorphism, crustal melting, and the onset of exhumation. Rocks within the dome and in the hanging wall contain garnet, kyanite, and staurolite porphyroblasts that record prograde Barrovian metamorphism during crustal thickening that reached ～530–630 °C and ～7–8 kbar, ending by c. 30 Ma. Cordierite and sillimanite overgrowths on Barrovian assemblages within the dome record dominantly top‐down‐to‐the‐west shearing during near‐isothermal decompression of the footwall rocks to ～4 kbar by 23 Ma during an exhumation rate of 1.3 mm year^?1. Monazite growth accompanied Barrovian metamorphism and decompression. The leucogranite injection complex within the dome initiated at 23 Ma and continued to 18 Ma. These data show that orogen‐parallel extension in this part of the Himalaya occurred earlier than previously documented (>16 Ma). Contemporaneous onset of near‐isothermal decompression, top‐down‐to‐the‐west shearing, and injection of the decompression‐driven leucogranite complex suggests that early crustal melting may have created a weakened crust that was proceeded by localization of strain and shear zone development. Exhumation along the shear zone accommodated decompression by 23 Ma in a kinematic setting that favoured orogen‐parallel extension.     

Drought Forecasting in a Semi-arid Watershed Using Climate Signals:a Neuro-fuzzy Modeling Approach

Large-scale annual climate indices were used to forecast annual drought conditions in the Maharlu-Bakhtegan watershed,located in Iran,using a neuro-fuzzy model.The Standardized Precipitation Index（SPI） was used as a proxy for drought conditions.Among the 45 climate indices considered,eight identified as most relevant were the Atlantic Multidecadal Oscillation（AMO）,Atlantic Meridional Mode（AMM）,the Bivariate ENSO Time series（BEST）,the East Central Tropical Pacific Surface Temperature（NINO 3.4）,the Central Tropical Pacific Surface Temperature（NINO 4）,the North Tropical Atlantic Index（NTA）,the Southern Oscillation Index（SOI）,and the Tropical Northern Atlantic Index（TNA）.These indices accounted for 81% of the variance in the Principal Components Analysis（PCA） method.The Atlantic surface temperature（SST：Atlantic） had an inverse relationship with SPI,and the AMM index had the highest correlation.Drought forecasts of neuro-fuzzy model demonstrate better prediction at a two-year lag compared to a stepwise regression model.