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.     

Modeling Local Tidal Constituents Using TOPEX/POSEIDON SSH for the Malaysian Marine Region

A new Local Ocean Tide Model, has been produced for the Exclusive Economic Zone (EEZ) of Malaysia, which incorporates some of the latest TOPEX/POSEIDON data for the years 1992 to 1998. Local tide gauge data are used as a comparison, along with another leading Global Ocean Tide Model, Ori96. The leading diurnal and semidiurnal constituents M₂, S₂, N₂, K₁, O₁, P₁ and Q₁ are reproduced using TOPEX/POSEIDON Sea Surface Heights (SSH) in a response analysis type least squares derivation following Munk and Cartwright (1966).     

Multiple-point geostatistical modeling based on the cross-correlation functions

An important issue in reservoir modeling is accurate generation of complex structures. The problem is difficult because the connectivity of the flow paths must be preserved. Multiple-point geostatistics is one of the most effective methods that can model the spatial patterns of geological structures, which is based on an informative geological training image that contains the variability, connectivity, and structural properties of a reservoir. Several pixel- and pattern-based methods have been developed in the past. In particular, pattern-based algorithms have become popular due to their ability for honoring the connectivity and geological features of a reservoir. But a shortcoming of such methods is that they require a massive data base, which make them highly memory- and CPU-intensive. In this paper, we propose a novel methodology for which there is no need to construct pattern data base and small data event. A new function for the similarity of the generated pattern and the training image, based on a cross-correlation (CC) function, is proposed that can be used with both categorical and continuous training images. We combine the CC function with an overlap strategy and a new approach, adaptive recursive template splitting along a raster path, in order to develop an algorithm, which we call cross-correlation simulation (CCSIM), for generation of the realizations of a reservoir with accurate conditioning and continuity. The performance of CCSIM is tested for a variety of training images. The results, when compared with those of the previous methods, indicate significant improvement in the CPU and memory requirements.     

Numerical model of antecedent rainfall effect on slope stability at a hillslope of weathered granitic soil formation

Modeling rainwater infiltration in slopes is vital to the analysis of slope failure induced by heavy rainfall. Although the significance of rainwater infiltration in causing landslides is widely recognized, there have been different conclusions as to the relative roles of antecedent rainfall to slope failure. In this study, a numerical model was developed to estimate the effect of antecedent rainfall on an unsaturated slope, the formation of a saturated zone, and the change in slope stability under weak rainfall and rainstorm event. Results showed that under a rainstorm event, slope failure occurred at comparably similar time although the antecedent rainfall drainage periods prior to the major rainfall were different (i.e., 24-h, 48-h and 96-h). However, under weak rainfall condition, differences of the antecedent rainfall drainage periods have significant effect on development of pore-water pressure. A higher initial soil moisture conditions caused faster increase in pore water pressure and thus decreasing the safety factor of the slope eventually increasing likelihood of slope failure.     

Elemental characterization of black shales of Khyber Pakthunkhawa (KPK) region of Pakistan using AAS

Black shale samples were collected from Chimiari Khyber Pakthunkhawa region of Pakistan and were analyzed for elemental compositions. Atomic Absorption Spectrophotometry (AAS) was utilized for the determination of elements in the digested solutions. The analysis of black shale was performed precisely with relative standard deviation (RSD) lower than 2%. Results showed that the samples contained high concentrations of Ca (11.98 %), Al (7.09%), Fe (3.03%), Mg (0.59%) and Ti (0.58%).     

Determining the origin of volcanic rocks in the mélange complex of Karangsambung based on the electrical resistivity imaging

An ENE-WSW-trending localized basalt-diabase outcrop along the SE margin of Luk Ulo Mélange Complex has been suggested as intrusive rocks cut through the Paleogene Totogan and Karangsambung formations. However, the absolute dating of the volcanics is older than the inferred relative age of the sedimentary formations, hence the in-situ intrusion theory is less likely. A subsurface imaging should delineate the possibility of the in-situ nature of volcanic rock by looking at the continuation of the rocks to the depth. In this study, we did a subsurface imaging by electrical resistivity method. The electrical resistivity surveys were conducted at 3 (three) lines across the ENE-WSW trend of the volcanic distribution. From those three measurements, we obtained three inversion models that present the distribution of the resistivity. We could differentiate between the high resistivity of volcanic rocks and the low resistivity of the clay-dominated sediments. Instead of the deep-rooted intrusions, the geometry of the volcanic rocks is concordant with the sedimentary strata. Since we do not observe any spatial continuity of the bodies, both laterally and vertically, the volcanic rocks might be part of broken intrusive rocks. Furthermore, the size and the sporadically distributed of the rocks also indicated that they are more likely as fragments during the olistostrome deposition, transported from its original location.     

Risk-based flood zoning employing expected annual damages: the Chenab River case study

Flooding proves to be the most devastating and annihilating natural hazard in Pakistan. Existing flood management strategies are riveted primarily to the structural measures that contribute limited loss reduction capability at the national level. Non-structural measures are not part of regular practices, as the adopted design standards, which are probabilistic in nature, are unable to assess their feasibilities. An improved risk-based assessment using expected annual damages (EAD) is introduced in this article for the evaluation of combined impacts. EAD treat the probabilistic nature of losses and provide an extended visualization of risk distributions in the form of damage curves and expected annual damages distribution maps. The Chenab River floodplain was selected to study the coalesced response of embankments and flood zoning, preliminary in economic terms. In this regard, the impacts of all likely floods are considered instead of the traditional focus on a single design flood. Damage curves and maps are compiled using estimated losses and probabilities of all floods. Flood zoning for agricultural land is performed. The results support choosing a multidirectional conjunctive approach that considers multiple measures to reduce flood losses. These results can be used as a vital input for the decision-making process.     

Screening Indian Mustard Genotypes for Phytoremediating Arsenic‐Contaminated Soils

Ten Indian mustard (Brassica juncea L.) genotypes were screened for their phytoremediation potential for arsenic (As) contaminated water under laboratory‐controlled conditions. The genotypes were grown in a hydroponic chamber for 20 days in 250‐mL beakers containing As‐contaminated water. During plant development, changes in plant growth, biomass, and total As were evaluated. Of the 10 genotypes (Pusa Agrani, BTO, Pusa Kranti, Pusa Bahar, Pusa Bold, Pusa Basant, Pusa Jai Kisan, Arka Vardhan, Varuna, and Vaibhav) Pusa Jai Kisan was the most effective in phytoremediating As‐contaminated water under hydroponic conditions. This will provide new information for Indian mustard genotypes for phytoremediating As‐contaminated soils.     

Hot plasma waves in Schwarzschild magnetosphere

In this paper we examine the wave properties of a hot plasma living in a Schwarzschild magnetosphere. The 3+1 GRMHD perturbation equations are formulated for this scenario. These equations are Fourier analyzed and then solved numerically to obtain the dispersion relations for a non-rotating, rotating non-magnetized and rotating magnetized plasma. The wave vector is evaluated, which is used to calculate the refractive index. These quantities are shown in graphs which are helpful to discuss the dispersive properties of the medium near the event horizon.     

Modelling interannual variations in catchment evapotranspiration considering vegetation and climate seasonality using the Budyko framework

The Budyko framework is an efficient tool for investigating catchment water balance, focusing on the effects of seasonal changes in climate (S) and vegetation cover (M) on catchment evapotranspiration (ET). However, the effects of vegetation seasonality on ET remain largely unknown. The present study explored these effects by modelling interannual variations in ET considering vegetation and climate seasonality using the Budyko framework. Reconstructed 15-day GIMMS NDVI_3g timeseries data from 1982 to 2015 were used to estimate M and extract the relative duration of the vegetation growing season (GL) in the Yellow River Basin (YRB). To characterize S, seasonal variations in precipitation and potential ET were extracted using a Gaussian algorithm. Analysis of the observed datasets for 19 catchments revealed that interannual variation in the catchment parameter ϖ (in Fuh's equation) was significantly and positively correlated with M and GL. Conversely, ϖ was significantly but negatively correlated with S. Furthermore, stepwise linear regression was used to calibrate the empirical formula of ϖ for these three dimensionless parameters. Following validation, based on observations in the remaining 11 catchments, ϖ was integrated into Fuh's equation to accurately estimate annual ET. Over 79% subcatchments showed an upward trend (0.9 mm yr⁻¹), whereas fewer than 21% subcatchments showed a downward trend (−0.5 mm yr⁻¹) across YRB. In the central region of the middle reach, ET increased with increased M, prolonged GL, and decreased S, whereas in the source region of YRB, ET decreased with decreased M and shortened GL. Our study provides an alternative method to estimate interannual ET in ungauged catchments and offers a novel perspective to investigate hydrological responses to vegetation and climate seasonality in the long-term.