On integral approach to regional gravity field modelling from satellite gradiometric data

Solution of the gradiometric boundary value problems leads to three integral formulas. If we are satisfied with obtaining a smooth solution for the Earth’s gravity field, we can use the formulas in regional gravity field modelling. In such a case, satellite gradiometric data are integrated on a sphere at satellite level and continued downward to the disturbing potential (geoid) at sea level simultaneously. This paper investigates the gravity field modelling from a full tensor of gravity at satellite level. It studies the truncation bias of the integrals as well as the filtering of noise of data. Numerical studies show that by integrating T _zz with 1 mE noise and in a cap size of 7°, the geoid can be recovered with an error of 12 cm after the filtering process. Similarly, the errors of the recovered geoids from T _xz,yz and T _{xx-yy, 2xy} are 13 and 21 cm, respectively.     

On a damaging earthquake-induced landslide in the Algerian Alps: the March 20, 2006 Laâlam landslide (Babors chain,northeast Algeria), triggered by the Kherrata earthquake (<Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> = 5.3)

On March 20, 2006, an earthquake (M _w = 5.3; SED) struck the mountainous region of the Babors chain (Wilaya of Bejaia, northeast Algeria). The seismic epicenter was located near the Kherrata village. This earthquake was felt on a large area of the northeastern part of Algeria. It reached an intensity of VII (EMS scale) at the Laalam village, situated at about 20 km northeast of Kherrata. Here, many old and recent houses were damaged or collapsed totally, four people died and 68 were injured. Field investigations revealed that these casualties were caused by a landslide triggered by the earthquake. Many fissures were visible on ground throughout the site. They were generated by both sliding and settling phenomena. The Laalam site is prone to landslide, as revealed by some evidences on old instabilities. This is due to two main factors: local geomorphology and geology. These factors intervene synchronously for reducing the slope instability at the Laalam village. The March 20, 2006 Kherrata earthquake was the trigger that released the Laalam landslide.     

The maximum expected earthquake magnitudes in different future time intervals of six seismotectonic zones of Iran and its surroundings

One of the crucial components in seismic hazard analysis is the estimation of the maximum earthquake magnitude and associated uncertainty. In the present study, the uncertainty related to the maximum expected magnitude μ is determined in terms of confidence intervals for an imposed level of confidence. Previous work by Salamat et al. (Pure Appl Geophys 174:763-777, 2017) shows the divergence of the confidence interval of the maximum possible magnitude m_max for high levels of confidence in six seismotectonic zones of Iran. In this work, the maximum expected earthquake magnitude μ is calculated in a predefined finite time interval and imposed level of confidence. For this, we use a conceptual model based on a doubly truncated Gutenberg-Richter law for magnitudes with constant b-value and calculate the posterior distribution of μ for the time interval T_f in future. We assume a stationary Poisson process in time and a Gutenberg-Richter relation for magnitudes. The upper bound of the magnitude confidence interval is calculated for different time intervals of 30, 50, and 100 years and imposed levels of confidence α?=?0.5, 0.1, 0.05, and 0.01. The posterior distribution of waiting times T_f to the next earthquake with a given magnitude equal to 6.5, 7.0, and 7.5 are calculated in each zone. In order to find the influence of declustering, we use the original and declustered version of the catalog. The earthquake catalog of the territory of Iran and surroundings are subdivided into six seismotectonic zones Alborz, Azerbaijan, Central Iran, Zagros, Kopet Dagh, and Makran. We assume the maximum possible magnitude m_max?=?8.5 and calculate the upper bound of the confidence interval of μ in each zone. The results indicate that for short time intervals equal to 30 and 50 years and imposed levels of confidence 1???α?=?0.95 and 0.90, the probability distribution of μ is around μ?=?7.16???8.23 in all seismic zones.     

Fracture characterization using converted waves

We present a method for inversion of fracture compliance matrix components from wide‐azimuth noisy synthetic PS reflection data and quantitatively show that reflection amplitude variations with offset and azimuth for converted PS‐waves are more informative than P‐waves for fracture characterization. We consider monoclinic symmetry for fractured reservoir (parameters chosen from Woodford Shale), which can be formed by two or more sets of vertical fractures embedded in a vertically transverse isotropic background. Components of effective fracture compliance matrices for a medium with monoclinic symmetry are related to the characteristics of the fractured medium. Monte Carlo simulation results show that inversion of PS reflection data is more robust than that of PP reflection data to uncertainties in our a priori knowledge (vertically transverse isotropic parameters of unfractured rock) than PP reflection data. We also show that, while inversion of PP reflections is sensitive to contrasts in elastic properties of upper and lower media, inversion of PS reflections is robust with respect to such contrasts.     

The 27 August 2010 Mw 5.7 Kuh-Zar earthquake (Iran): field investigation and strong-motion evidence

In this study, seismological aspects and field observation of the 2010 Kuh-Zar earthquake has been investigated. The Kuh-Zar earthquake, of magnitude 5.7 (Mw), occurred in northeastern Iran on August 27, 2010. The area is surrounded by branches of the active faults which are coated by the quaternary alluvium. During the past several decades, this area has been struck by a number of earthquakes. This earthquake with a moderate magnitude caused a higher rate of damage contrasted with previous earthquakes of the same magnitude range in Iran. Fortunately, the source of the Kuh-Zar earthquake was in a sparsely populated area, and therefore, it caused a few loss of life with the highest observed intensity of shaking VII (modified Mercalli intensity) in the Kuh-Zar village. The shock killed 4 people, injured 40, and destroyed more than 12 villages. According to the field observation, the mechanism of this shock is defined as a left-lateral strike slip. We also checked out the properties of strong ground motions in this earthquake using the records availed by Iranian strong motion network. At KUZ station, about 7 km east of the epicenter, the recorded PGA and PGV in both horizontal and vertical components were remarkably large for an event of this size, and visual inspection of the velocity time history reveals a pulse-like shape. Unfortunately no other recording stations were located close enough to the fault to capture a directivity pulse. Finally, according to the strong-motion properties and observed information, ShakeMaps of the earthquake have been generated by the native intensity observations and the recorded strong motions.     

Recent development of the earthquake strong motion-intensity catalog and intensity prediction equations for Iran

This study aims to develop a new earthquake strong motion-intensity catalog as well as intensity prediction equations for Iran based on the available data. For this purpose, all the sites which had both recorded strong motion and intensity values throughout the region were first searched. Then, the data belonging to the 306 identified sites were processed, and the results were compiled as a new strong motion-intensity catalog. Based on this new catalog, two empirical equations between the values of intensity and the ground motion parameters (GMPs) for the Iranian earthquakes were calculated. At the first step, earthquake “intensity” was considered as a function of five independent GMPs including “Log (PHA),” “moment magnitude (M_W),” “distance to epicenter,” “site type,” and “duration,” and a multiple stepwise regression was calculated. Regarding the correlations between the parameters and the effectiveness coefficients of the predictors, the Log (PHA) was recognized as the most effective parameter on the earthquake “intensity,” while the parameter “site type” was removed from the equations since it was determines as the least significant variable. Then, at the second step, a simple ordinary least squares (OLS) regression was fitted only between the parameters intensity and the Log (PHA) which resulted in more over/underestimated intensity values comparing to the results of the multiple intensity-GMPs regression. However, for rapid response purposes, the simple OLS regression may be more useful comparing to the multiple regression due to its data availability and simplicity. In addition, according to 50 selected earthquakes, an empirical relation between the macroseismic intensity (I₀) and M_W was developed.     

Composition and structure of species along altitude gradient in Moghan-Sabalan rangelands,Iran

This study provides a checklist of species distributed at the altitude gradient of Moghan-Sabalan rangelands in Ardabili province, Northwest Iran. We evaluated the changes in species composition, growth types of species, Raunkiaer’s life forms, geographical distribution, threat and endemicity status, and palatability of species along two altitudinal gradients in the sampling plots, which were conducted in eleven sites/habitats with 300 meters above sea level (masl) altitude intervals (from 100 to 3300 masl). We assessed the plant species composition with special reference to the gradient analysis, and identified overall 396 species, which was comprising 44 families and 194 genera. Results showed that Asteraceae family is by far the most species-rich family, followed by Poaceae, Fabaceae, Caryophyllaceae and Brassicaceae. Among the genera, Astragalus is the most diverse genus, followed by Allium, Veronica and Bromus, Galium, Silene and Ranunculus. Results indicated that the number of species increased as the altitude increased to 1200–1500 masl, but then starts to decline to 3300 masl. Family-to-genera ratio was 1:4.4, the family-tospecies ratio was 1:9, and the genera-to-species ratio was 1:2.04. Growth type of species analysis shows that the frequency of perennial plants was higher in the study area followed by annual species while the lower group was biennial species. The number of annuals showed a decreasing trend towards higher altitude. Hemicryptophytes and therophytes were the most frequent life forms constituted each with (41.9%). Hemicryptophytes showed an increasing trend with altitude, while therophytes showed a decreasing trend with altitude increase, followed by geophytes, chamaephytes, and phanerophytes. Results showed more than half of the species of the study area belonged to Iran-Turanian region and these species showed an increasing trend with altitude. In contrast, Sahara-Sindian species comprise a minor component of the spectrum, with decreasing trend with altitude. The rare and endangered species out of the surveyed taxonomic groups comprised 53 species in total which 29 of them are considered lower risk (LR), 13 data deficient (DD), 5 vulnerable (Vu) and with 3 rare (R) and identified endemic plants comprised 24 species. Some 56.6% species were identified as class III, 22.6% were class I and 20.8% were class II as the palatability variation. Moghan-Sabalan rangelands require strong conservation management policies in case of species loss and changing natural communities due to the occurrence of conversion into cropland, over-grazing and other anthropogenic effects.     

A three-dimensional drift pushover method for unsymmetrical plan buildings

The drift pushover analysis method for tall and regular buildings is extended in this paper to the third dimension. The focus of study is on the structures with important torsional response. For this purpose, 10, 15, 20 and 30-story steel moment frame buildings having unsymmetrical plans with 5–30% eccentricity ratios are studied. For evaluation of accuracy, nonlinear dynamic response of the buildings is determined under a consistent suit of earthquake ground motions. The maxima of the story drifts and shears and cumulative plastic hinge rotations of stories are calculated under the ground motions and their averages along with those of the modal pushover procedure are compared with the results of the presented method. The comparative analysis establishes the good accuracy of the three dimensional drift pushover method.     

Bearing capacity of non-associative coaxial granular materials by upper bound limit analysis and finite elements

In this paper, an upper bound estimate of the limit load on non-associative coaxial granular materials is presented. The kinematic approach of the upper bound limit analysis has been utilised. The failure mechanism is assumed to coincide with the direction of the shear bands at every point throughout the body. The shear band orientation in non-associative coaxial materials, i.e. those with the same major principal stress and major principal strain increment directions, can be found based on the angle of dilation and the major principal stress direction. Therefore, having known the stress field at limiting equilibrium, the orientation of the shear bands and hence, the failure mechanism can be obtained. In this study, the stress field is first determined by the method of stress characteristics. Then, the finite element interpolation technique is used to interpolate the stress field and to find the orientation of the shear bands at every point within the field. Once the failure mechanism and the stress state at every point along velocity discontinuities have been found, the upper bound limit analysis has been performed to estimate the limit load.     

Dynamic earth pressure on rigid retaining walls induced by a neighboring machine foundation,by the meshless local Petrov-Galerkin method

Dynamic earth pressure induced by machine foundations on a neighboring retaining wall is analyzed with emphasis on factors which control the intensity and location of the design forces. The meshless local Petrov-Galerkin(MLPG) method is used to analyze the problem for a variety of retaining wall and machine foundation geometries. The soil medium is assumed to be homogeneous and visco-elastic. The machine foundation is idealized as a harmonic sinusoidal dynamic force often encountered in practice. A number of analyses have been made to reveal the effect of the loading frequency, the location and size of the foundation and the soil shear wave velocity on the distribution and magnitude of the dynamic earth pressure. Results indicate that there is a critical frequency and a critical location for which the passive pressure takes the maxima in the entire duration of the dynamic load.