Seismicity patterns associated with the September 10th, 2008 Qeshm earthquake, South Iran

The b value of the Gutenberg-Richter relation and the standard deviate, Z, were calculated to investigate the temporal and spatial variations in seismicity patterns associated with the September 10th, 2008 (Mw?=?6.1) Qeshm earthquake. The temporal variations of b value illustrate a distinct dramatic drop preceding the Qeshm earthquake, and the spatial changes in b value highlight a zone with an abnormally low b value around the epicenter of this event. The cumulative number and Z value as a function of time show a precursory seismic quiescence preceding the 2008 Qeshm earthquake that observed for 1?year in a circle with R?=?50?km around its epicenter. The spatial distribution map of the standard deviate, Z, also exhibits an obvious precursory seismic quiescence region before the 2008 Qeshm event around the epicenter of this event. Interestingly, the precursory seismic quiescence region is approximately consistent with low b value anomaly region, and both have E–W to NE–SW trend. These two precursory anomalies took place in relatively large regions, which were possibly relevant to the preparation zone of the 2008 Qeshm event.     

Assessing Lateral Continuity within the Yammama Reservoir in the Foroozan Oilfield, Offshore Iran: An Integrated Study

Located in Iranian sector of the Persian Gulf, Foroozan Oilfield has been producing hydrocarbons via seven different reservoirs since the 1970s. However, understanding fluid interactions and horizontal continuity within each reservoir has proved complicated in this field. This study aims to determine the degree of intra-reservoir compartmentalization using gas geochemistry, light hydrocarbon components, and petroleum bulk properties, comparing the results with those obtained from reservoir engineering indicators. For this purpose, a total of 11 samples of oil and associated gas taken from different producing wells in from the Yammama Reservoir were selected. Clear distinctions, in terms of gas isotopic signature and composition, between the wells located in northern and southern parts of the reservoir (i.e. lighter δ13C1, lower methane concentration, and negative sulfur isotope in the southern part) and light hydrocarbon ratios (e.g. nC7/toluene, 2,6-dmC7/1,1,3-tmcyC5 and m-xylene/4-mC8) in different oil samples indicated two separate compartments. Gradual variations in a number of petroleum bulk properties (API gravity, V/Ni ratios and asphaltene concentration) provided additional evidence on the reservoir-filling direction, signifying that a horizontal equilibrium between reservoir fluids across the Yammama Reservoir is yet to be achieved. Finally, differences in water-oil contacts and reservoir types further confirmed the compartmentalization of the reservoir into two separate compartments.     

Fast hyperbolic deconvolutive Radon transform using generalized Fourier slice theorem

The hyperbolic Radon transform has a long history of applications in seismic data processing because of its ability to focus/sparsify the data in the transform domain. Recently, deconvolutive Radon transform has also been proposed with an improved time resolution which provides improved processing results. The basis functions of the (deconvolutive) Radon transform, however, are time-variant, making the classical Fourier based algorithms ineffective to carry out the required computations. A direct implementation of the associated summations in the time–space domain is also computationally expensive, thus limiting the application of the transform on large data sets. In this paper, we present a new method for fast computation of the hyperbolic (deconvolutive) Radon transform. The method is based on the recently proposed generalized Fourier slice theorem which establishes an analytic expression between the Fourier transforms associated with the data and Radon plane. This allows very fast computations of the forward and inverse transforms simply using fast Fourier transform and interpolation procedures. These canonical transforms are used within an efficient iterative method for sparse solution of (deconvolutive) Radon transform. Numerical examples from synthetic and field seismic data confirm high performance of the proposed fast algorithm for filling in the large gaps in seismic data, separating primaries from multiple reflections, and performing high-quality stretch-free stacking.     

Evaluating sampling locations in river water quality monitoring networks: application of dynamic factor analysis and discrete entropy theory

In this paper, a methodology is proposed for evaluating sampling locations in an existing river water quality monitoring network. The dynamic factor analysis is utilized to extract the independent dynamic factors from time series of water quality variables. Then, the entropy theory is applied to the independent dynamic factors to construct transinformation–distance (T–D) curves. The computation time in the case of using dynamic factors is significantly less than when the raw data is used because the number of independent dynamic factors is usually much less than the number of monitored water quality variables. In this paper, it is also shown that by clustering the study area to some homogenous zones and developing T–D curves for each zone, the accuracy of the results is significantly increased. To evaluate the applicability and efficiency of the proposed methodology, it is applied to the Karoon River which is the most important river system in Iran.     

Characterisation of the Sarcheshmeh copper mine tailings,Kerman province,southeast of Iran

Mineral processing operation at the Sarcheshmeh porphyry copper mine has produced huge quantities of tailings materials containing sulphide minerals in particular pyrite. These tailings materials were geochemically and mineralogically characterised to assess pyrite and chalcopyrite oxidation, acid mine drainage generation, and trace element mobility to lead development of a proper remediation plan. Five vertical trenches up to 4.2 m deep were excavated from the tailings surface, and 70 solid samples were taken in 0.3 m intervals. The samples were first mineralogically analysed. Pyrite was the main sulphide mineral found in the tailings. The gangue minerals include quartz ± muscovite–illite ± chlorite ± albite ± orthoclase ± halite. The samples were geochemically analysed for total concentrations of 62 elements, paste pH, SO₄ ^2?, CO₃ ^2?, and HCO₃ ^?. The maximum concentrations of SO₄ ^2? (1,300, 1,170, 1,852, 1,960 and 837 mg/L) were observed at a depth of 0.9 m in profiles A, B, C, D and E, respectively. The tailings have a high acid-producing potential and low acid-neutralising potential (pyrite 4–6 wt %, calcite 1 wt %). Fe₂(SO₄)₃, CuSO₄, MgSO₄ and MnSO₄ were the dominant secondary sulphate minerals in the tailings. The lowest pH values (2.9, 3 and 3) were measured at a depth of 0.3 m in the profiles A, B and C, 3.9 at a depth of 0.6 m in the profile D and 3 at a depth of 0.9 m in the profile E. The upper portions of the profiles C (1.8 m) and D (2.1 m) were moderately oxidised, while oxidation in the profiles A, B and E did not extend more than 1.2, 1.2 and 1.5 m beneath the tailings surface. Zn, Pb, Rb, U, Hf, Nd, Zr and Ga show almost a constant trend with depth. Cd, Sr, Th, La and Ce increased with increasing depth of the tailings materials while, Co, V, Ti, Cr, Cu, As, Mn, Ag, Mo and Ni exhibit initially a decreasing trend from tailings surface to the depths that vary between 0.9 and 1.2. They then remained constant with the depth. The results show pyrite and chalcopyrite oxidation at surface layers of the tailings and subsequent leaching of the oxidation products and trace elements by infiltrated atmospheric precipitation.     

Volume,gravitational potential energy reduction,and regional centroid position change in the wake of landslides triggered by the 14 April 2010 Yushu earthquake of China

In recent years, earthquake-triggered landslides have attracted much attention in the scientific community as a main form of seismic ground response. However, little work has been performed concerning the volume and gravitational potential energy reduction of earthquake-triggered landslides and their severe effect on landscape change. This paper presents a quantitative study on the volume, gravitational potential energy reduction, and change in landscape related to landslides triggered by the 14 April 2010 Yushu earthquake. At least 2,036 landslides were triggered by the earthquake. A total landslide scar area of 1.194 km² was delineated from the visual interpretation of aerial photographs and satellite images and was supported by selected field checking. In this paper, we focus on possible answers to the following five questions: (1) What is the total volume of the 2,036 landslides triggered by the earthquake, and what is the average landslide erosion thickness in the earthquake-stricken area? (2) What are the elevations of all landslide materials in relation to pre- and post-landsliding? (3) How much was the gravitational potential energy reduced due to the sliding of these landslide materials? (4) What is the average elevation change caused by these landslides in the study area? (5) What is the vertical change of the regional centroid position above sea level, as induced by these landslides? It is concluded that the total volume of the 2,036 landslides is 2.9399?×?10⁶ m³. The landslide erosion thickness throughout the study area is 2.02 mm. The materials of these landslides moved from an elevation of 4,145.243 to 4,104.697 m, resulting in a decreased distance of 40.546 m. The gravitational potential energy reduction related to the landslides triggered by the earthquake was 2.9213?×?10¹² J. The average regional elevation of the study area is 4,427.160 m, a value consistent with the assumption that the accumulated materials were remained in situ. This value changes from 4,427.160 to 4,427.158 m with all landslide materials moved out of the study area, resulting in a reduction in elevation of 2 mm. Based on the assumption that all landslide materials moved out of the study area, the elevations of the centroid of the study area’s crust changed from 2,222.45967 to 2,222.45867 m, which means the centroid value decreased by 1 mm. This value is 0.001 mm when assuming that the materials were remained in situ, which is almost negligible, compared with the situation of “all landslide materials moved out of the study area.”     

Non-parametric trend analysis of the aridity index for three large arid and semi-arid basins in Iran

Currently, an important scientific challenge that researchers are facing is to gain a better understanding of climate change at the regional scale, which can be especially challenging in an area with low and highly variable precipitation amounts such as Iran. Trend analysis of the medium-term change using ground station observations of meteorological variables can enhance our knowledge of the dominant processes in an area and contribute to the analysis of future climate projections. Generally, studies focus on the long-term variability of temperature and precipitation and to a lesser extent on other important parameters such as moisture indices. In this study the recent 50-year trends (1955–2005) of precipitation (P), potential evapotranspiration (PET), and aridity index (AI) in monthly time scale were studied over 14 synoptic stations in three large Iran basins using the Mann–Kendall non-parametric test. Additionally, an analysis of the monthly, seasonal and annual trend of each parameter was performed. Results showed no significant trends in the monthly time series. However, PET showed significant, mostly decreasing trends, for the seasonal values, which resulted in a significant negative trend in annual PET at five stations. Significant negative trends in seasonal P values were only found at a number of stations in spring and summer and no station showed significant negative trends in annual P. Due to the varied positive and negative trends in annual P and to a lesser extent PET, almost as many stations with negative as positive trends in annual AI were found, indicating that both drying and wetting trends occurred in Iran. Overall, the northern part of the study area showed an increasing trend in annual AI which meant that the region became wetter, while the south showed decreasing trends in AI.     

Application of singular value decomposition (SVD) and semi-discrete decomposition (SDD) techniques in clustering of geochemical data: an environmental study in central Iran

Common multivariate clustering techniques are ineffective in identifying subtle patterns of correlation, and clustering of variables or samples within complex geochemical datasets. This study compares the combination of singular value decomposition (SVD) and semi discrete decomposition (SDD), with that of hierarchical cluster analysis (HCA), to examine patterns within a multielement soil geochemical dataset from an agricultural area in the vicinity of Pb–Zn mining operations in central Iran. SVD was used to both identify patterns of correlation between variables and samples and to “denoise” the data, and SDD to simultaneously cluster the samples and variables. The results reveal various spatial associations of mining waste-associated metals As, Ba, Pb and Zn, and within the remaining elements whose distribution is largely controlled by the major oxides. SVD–SDD was found to be superior to HCA, in its ability to detect subtle clusters in soil geochemistry indicative of mine-related contamination in the study area.     

Preparation of frost atlas using different interpolation methods in a semiarid region of south of Iran

In this research, suitability of different kriging and inverse distance weighted (IDW) methods in estimating occurrence date of frost was evaluated. Data included minimum daily air temperature values from 27 meteorological stations of Fars province in southern Iran from 18 to 45 years. Data ranges of 0 to ?1.5, ?1.5 to ?3 and below ?3°C were considered as mild, moderate and severe frost intensities, respectively. Starting with the first day of autumn, iso-occurrence days for the frost intensities and occurrence probabilities (25%, 50%, 75% and 90%) were estimated using ordinary kriging, cokriging, residual kriging type 1 (RK1), residual kriging type 2 (RK2), universal kriging and IDW methods. In these models, the errors of estimated frost intensities at different probabilities were lowest in the RK2 model, but lack of establishment of spatial structure due to long distance between stations caused the predictions not to be acceptable in some cases. In a proposed method (modified inverse distance weighted, MIDW), the trend between the first and last days of frost occurrence with earth elevation was removed, and the reminder values were estimated by (IDW) method. Although, the errors for estimated frost dates by MIDW and RK2 methods were the same, but the MIDW method did not have the spatial establishment shortcoming. Furthermore, the simplicity and practicality of the MIDW method makes it a reasonable selection.