Ranking the rock slope instability potential using the Interaction Matrix (IM) technique; a case study in Iran

All the conventional techniques for the analysis of slope stability ranging from simple kinematic analysis using stereonets, to the various widely used limit equilibrium methods, to sophisticated numerical methods belong to a category that are generally known as the analytic approaches and thus are only able to consider a limited number of affecting factors and then solve the problem in details. In contrast, the systems approaches not only can examine the problem in its totality with a complete list of the components, but also can take the interactions between the factors into account. This paper presents a complete application of a well-known systems technique named the Interaction Matrix (IM) in ranking the instability potential of rock slopes of the Khosh-Yeylagh Main Road, Iran as the case study of the research. For this purpose, 15 stations have been selected and a relatively comprehensive database containing the fieldwork information has been constructed. Following the IM technique, the most important factors relating to the general environment and to the rock mass characteristics have been considered. Their reciprocal causes and effects have been analyzed in order to weight each parameter according to its degree of interactivity in the system. Then, the slope instability index has been calculated which refers to the inherent potential instability of each slope of the examined region. The final instability ranking has been presented for the investigated slopes in Khosh-Yeylagh Main Road based on a simple classification. The main aim of the study is to extend the use of systems approach and specifically the IM technique in slope stability analysis. Also, this research shows the importance of consideration of an approximately complete set of key parameters affecting the stability of rock slopes.     

Geodetic constraints on the period of episodic tremors and slips using least squares harmonic estimation method with application to Cascadia subduction zone

It is expected that episodic tremors and slips (ETS) as a type of slow earthquakes with periodic property will be detected using harmonic estimation techniques. The principal goal of this paper is the detection of these kinds of earthquakes using least squares harmonic estimation (LSHE). To accomplish this, the raw time series of 38 permanent GPS stations of the Pacific Northwest Geodetic Array have been analyzed. Previously, some independent techniques could confirm the occurrences of the aforementioned quakes at these stations. However, the current research intends to evaluate the spectrum of each of the de-trended time series using the LS-HE method. In each station, the period of the detected harmonic with the maximum power spectrum is equivalent to the average period previously reported for these events. According to the obtained results, the recurrence interval of these events ranges from 9 months to 3 years. In sum, the study confirms this method being efficient for investigating the occurrences of ETSs when the length of the GPS time series is sufficiently large.     

A Monte Carlo adapted finite element method for dislocation simulation of faults with uncertain geometry

Dislocation modelling of an earthquake fault is of great importance due to the fact that ground surface response may be predicted by the model. However, geological features of a fault cannot be measured exactly, and therefore these features and data involve uncertainties. This paper presents a Monte Carlo based random model of faults with finite element method incorporating split node technique to impose the effects of discontinuities. Length and orientation of the fault are selected as random parameters in the domain model, and hence geometrical uncertainties are encountered. Mean and standard deviation values, as well as probability density function of ground surface responses due to the dislocation are computed. Based on analytical and numerical calculation of dislocation, two approaches of Monte Carlo simulations are proposed. Various comparisons are examined to illustrate the capability of both methods for random simulation of faults.     

Identification of hydrogeochemical processes and pollution sources of groundwater resources in the Marand plain,northwest of Iran

The main aims of the present study are to identify the major factors affecting hydrogeochemistry of groundwater resources in the Marand plain, NW Iran and to evaluate the potential sources of major and trace elements using multivariate statistical analysis such as hierarchical clustering analysis (HCA) and factor analysis (FA). To achieve these goals, groundwater samples were collected in three sampling periods in September 2013, May 2014 and September 2014 and analyzed with regard to ions (e.g., Ca²⁺, Mg²⁺, Na⁺ and K⁺, HCO₃ ^?, SO₄ ^2?, Cl^?, F^? and NO₃ ^?) and trace metals (e.g., Cr, Pb, Cd, Mn, Fe, Al and As). The piper diagrams show that the majority of samples belong to Na–Cl water type and are followed by Ca–HCO₃ and mixed Ca–Na–HCO₃. Cross-plots show that weathering and dissolution of different rocks and minerals, ion exchange, reverse ion exchange and anthropogenic activities, especially agricultural activities, influence the hydrogeochemistry of the study area. The results of the FA demonstrate that 6 factors with 81.7% of total variance are effective in the overall hydrogeochemistry, which are attributed to geogenic and anthropogenic impacts. The HCA categorizes the samples into two clusters. Samples of cluster C1, which appear to have higher values of some trace metals like Pb and As, are spatially located at the eastern and central parts of the plain, while samples of cluster C2, which express the salinization of the groundwater, are situated mainly westward with few local exceptions.     

Effects of climatic factors and soil salinity on the distribution of vegetation types containing <Emphasis Type="Italic">Anabasis aphylla</Emphasis> in Iran: a multivariate factor analysis

In order to identify factors affecting the distribution of Anabasis aphylla in rangelands of Iran, 132 climatic parameters which contributed to the species distribution (based on the ecological conditions of the study area) were selected. Factor analysis was then applied to determine the most important factors affecting species distribution. In the next stage, the spatial distribution map of An. aphylla throughout Iran was prepared. According to the results of factor analysis, ten factors with Eigenvalues greater than one explained 92.96% of the total variance. These factors were temperature below zero, winter rainfall, summer rainfall, wind, sunlight, warm season dust, rainfall in fall, thunderstorms, relative humidity in September, and cloudy days in winter and explained 34.34, 9.71, 9.69, 8.85, 5.99, 5.35, 4.97, 4.78, 4.73, and 4.51% of the total variance, respectively. Moreover, six vegetation types containing An. aphylla were identified throughout the country. These types were distributed in the central areas of Iran and also in patches within the southeastern parts of the country. The results showed that in addition to the nine above climatic factors, soil variables particularly salinity influence establishing the vegetation types in their distribution areas and the vegetation types are located in low to relatively high salinity, and soils with these salinity levels (low to relatively high salinity) are suitable for these types of plants.     

Calcareous nannofossil biostratigraphy of the Paleocene at the Misheiti section,East Central Sinai,Egypt

Integration of calcareous nannofossil data, δ¹³C and δ¹⁸O values, and carbonate contents of the lower Paleocene–upper Paleocene sequence that crops out at the Misheiti section, East Central Sinai, Egypt, were used to denote the Danian/Selandian (D/S) and Selandian/Thanetian (S/T) stage boundaries. The study interval belongs to the Dakhla and Tarawan formations. Four calcareous nannofossil zones (NP4, NP5, NP6, and NP7/8) were recognized. The base of the Selandian Stage is tentatively placed at the lowest occurrences (LOs) of taxa ascribable to the second radiation of fasciculiths (i.e., Lithoptychius janii). This level is marked by a sudden drop of δ¹³C and δ¹⁸O values and carbonate content. No distinctive lithological changes were observed across the D/S boundary at the study section. A hiatus at the NP5/NP6 zonal boundary is indicated by the condensation of zones NP5 and NP6.The base of the Thanetian is placed at the base of Zone NP7/8 at the lithological change observed in correspondence to the boundary between the Dakhla and Tarawan formations. The δ¹³C and δ¹⁸O values abruptly decrease slightly above the base of Zone NP7/8. No consistent variations in the carbonate contents were recorded within Zone NP6 or across the NP6/NP7/8 zonal boundary.     

Scour and flow field around a spur dike in a 90° bend

Spur dike is an important element in fiver training that creates rapid variations in flow field, sediment transport and bed topography. The mechanism of flow and sediment transport in a channel bend is very complex, especially when a spur dike is constructed in a bend. Most of previous investigations on flow behavior and scour around spur dike were carried out in straight channels. In this paper results of experiments on flow field and scour around a spur dike in a 90 degree channel bend are presented, Sand with uniform grain size was used as the bed material, Experiments were conducted for different locations and different lengths of spur dikes at the bend with different values of discharge, The three dimensional flow fields around a spur dike were investigated, The maximum depth of scour was correlated to the Froude numbers, lengths and the locations of spur dike in the bend.     

Analytical‐numerical solution for stress distribution around tunnel reinforced by radial fully grouted rockbolts

This paper presents an analytical‐numerical approach to obtain the distribution of stresses and deformations around a reinforced tunnel. The increase in the radial stress of the reinforced tunnel, based on the performance of a bolt, is modeled by a function, which its maximum value is in the vicinity of the bolt periphery and it exponentially decreases in the far distance from the bolt. On the basis of this approach, the shear stiffness between the bolt and the rock mass and the shear stress distribution around the bolt within the rock mass are also analytically obtained. The results are compared with those obtained by the assumption of ‘uniform increase of radial stress’ method, which is made by the previous studies. The analyses show when the bolts' spacing is large, the safety factor must be increased if the ‘uniform increase of radial stress’ method is used for the design. Finally, a procedure is introduced to calculate the non‐equal deformation of the rock mass between the bolts at any radius that can be useful to compute the bending moment in shotcrete layer in New Austrian Tunnelling Method (NATM) approach. Copyright © 2016 John Wiley & Sons, Ltd.     

The Bam Earthquake of 26 December 2003

The devastating earthquake of 26 December 2003 claimed more than 26,000 lives in the city of Bam and surrounding towns and villages in Southeast Iran, and left the majority of the Bam population homeless. The reason for this tragedy was an unfortunate combination of geological, social and human circumstances. The causative fault practically traversed the city of Bam and the earthquake occurred at a shallow depth. The residential buildings were completely inappropriate for a seismic region, being extremely vulnerable to earthquake shaking, and the earthquake occurred early in the morning when most people were still sleeping. The damage pattern was nearly symmetric about a line 3 km to the west of the surface expression of the Bam fault, and the damage attenuated rapidly with distance from this line. The industrial facilities and the lifelines performed relatively well and experienced slight to moderate damage, but this might have been due to their distance from the earthquake epicentre. However, many of the qanat (traditional subterranean irrigation channels) chains that served the twin cities of Bam and Baravat collapsed. Emergency facilities (hospitals, police and fire stations), schools and the university were destroyed and/or heavily damaged during the earthquake. The geotechnical effects of the earthquake were not significant. There was little evidence that site response effects played a major role in the damage pattern in the city. There were no reports of liquefaction and only minor sliding activity took place during the event. A unique set of strong motion acceleration recordings were obtained at the Bam accelerograph station. The highest peak ground acceleration (nearly 1g) was recorded for the vertical component of the motion. However, the longitudinal component (fault-parallel motion in N–S direction) clearly had the largest energy flux, as well as the largest maximum velocity and displacement.