Spatial analysis of the frequency–magnitude distribution of aftershock activity of 2003 Bam earthquake: southeast Iran

The Bam earthquake (2003 December 26, M _W = 6.6) was one of the largest earthquakes that occurred in southeast of Iran during last century. It took place along an N–S trending right-lateral strike-slip fault, almost near the southern end of Nyband–Gowk fault. In this study, we mapped the frequency–magnitude distribution of aftershock events spatially across the Bam aftershock zone. The b-value varies between 0.6 and 1.1 across the Bam rupture zone. The overall depth distribution of b-value in Bam aftershock zone reveals two distinct increases in b-value: (1) at depths of 8–10 km and (2) shallower than 4 km beneath the Bam city. There is no correlation between high b- value anomalies found in this study and the region of largest slip, whereas the spatial correlation between high b-value anomalies and the zone of low V _s and high σ (in earlier tomography study) is obvious. This correlation reveals that material properties and increasing heterogeneity are more important in controlling b-value distribution in Bam earthquake rupture zone. The high b-value anomaly near the surface of northern part of rupture zone may be related to unconsolidated and water-rich quaternary alluvial sediments and probable low-strength rocks beneath them. The high b-value anomaly at depth range 8–10 km can be correlated with fractured and fluid-filled mass, which may result from the movement of magma during Eocene volcanism in the Bam area. In this study, the induced changes in pore fluid pressure due to main shock are suggested as a mechanism for aftershock generation.     

Mass,Volume and Velocity of the Antarctic Ice Sheet: Present-Day Changes and Error Effects

This study examines present-day changes of the Antarctic ice sheet (AIS) by means of different data sets. We make use of monthly gravity field solutions acquired by the Gravity Recovery and Climate Experiment (GRACE) to study mass changes of the AIS for a 10-year period. In addition to ‘standard’ solutions of release 05, solutions based on radial base functions were used. Both solutions reveal an increased mass loss in recent years. For a 6-year period surface-height changes were inferred from laser altimetry data provided by the Ice, Cloud, and land Elevation Satellite (ICESat). The basin-scale volume trends were converted into mass changes and were compared with the GRACE estimates for the same period. Focussing on the Thwaites Glacier, Landsat optical imagery was utilised to determine ice-flow velocities for a period of more than two decades. This data set was extended by means of high-resolution synthetic aperture radar (SAR) data from the TerraSAR-X mission, revealing an accelerated ice flow of all parts of the glacier. ICESat data over the Thwaites Glacier were complemented by digital elevation models inferred from TanDEM-X data. This extended data set exhibits an increased surface lowering in recent times. Passive microwave remote sensing data prove the long-term stability of the accumulation rates in a low accumulation zone in East Antarctica over several decades. Finally, we discuss the main error sources of present-day mass-balance estimates: the glacial isostatic adjustment effect for GRACE as well as the biases between laser operational periods and the volume–mass conversion for ICESat.     

Coda <Emphasis Type="Italic">Q</Emphasis> in Qeshm Island,south of Iran,using aftershocks of the Qeshm earthquake of November 27, 2005

Coda wave attenuation is estimated for Qeshm Island which is located in the southeastern part of Zagros. For this purpose, the aftershocks of Qeshm earthquake in November 27, 2005, recorded within an epicentral distance less than 100 km, have been used. More than 829 earthquakes were recorded by a local temporary network consisting of 16 short period stations installed after a week after the main shock for ～10 weeks. The coda quality factor, Q _c, was estimated using the single-backscattering model in frequency bands of 0.5–24 Hz. In this research, lateral and vertical variations of coda Q in Qeshm Island are explored. In Qeshm Island, absence of significant lateral variation of coda Q is observed. To investigate the attenuation variation with depth, the coda Q value was calculated for coda time windows with different lengths (5, 10, 15, 20, 25, and 30 s). It is observed that coda Q increases with depth. However, in our study area, the rate of increase of coda Q with depth is not uniform. Beneath Qeshm Island, the rate of increase of coda Q is greater at depths less than ～40 km compared with those of larger depths. This is indicating the existence of a low attenuation anomalous structure under the ～40-km depth which may be correlated with the Moho depth in this region. The average frequency relation for this region is Q _c = 36 ± 1.2f ^{0.94 ± 0.039} at a 5 s-lapse time window length and Q _c = 110 ± 1.8f ^{0.88 ± 0.09} at a 30-s lapse time window length.     

Assessing and Modeling the Impacts of Wetland Land Cover Changes on Water Provision and Habitat Quality Ecosystem Services

Natural Resources Research - Understanding the spatial–temporal dynamics of wetland land cover (LC) changes and their impacts on ecosystem services (ESs) is essential for wetland conservation...     

Using combined AHP–genetic algorithm in artificial groundwater recharge site selection of Gareh Bygone Plain,Iran

Flood spreading is one of the suitable strategies to control and benefit from floods which in turn improve the groundwater recharge, makes soil more fertile, and increases nutrients in soil. It is also a method for reusing sediment, which is usually wasted. Thus, selection of suitable areas for flood spreading and directing the flood water into permeable formations are amongst the most effective strategies in flood spreading projects. Having combined analytic hierarchy process (AHP) of multi-criteria decision analysis and genetic algorithm (GA) of artificial intelligence approaches, this paper addresses the problem of finding the most suitable area location for flood spreading operation in the Gareh Bygone Plain of Iran. To this end, the nine effective geodata layers including slope, alluvium thickness, geology, morphology, electrical conductivity, land use, drainage density, aquifer transmissivity, and elevation were prepared in geographic information system environment. This stage was followed by elimination of the exclusionary areas for flood spreading while determining the potentially suitable ones. Having closely examined the potentially suitable areas using the proposed methodology, the land suitability map for flood spreading was produced. The AHP and GA were used for ranking all the alternatives and weighting the criteria involved, respectively. The results of the study showed that most suitable areas for the artificial groundwater recharge are located in Quaternary Qft ₂ and Qsf geologic units and in morphological units of pediment and Alluvial fans with slopes not exceeding 2 %. Finally, further evidence for the acceptable efficiency of the integrated AHP–GA method in locating most suitable flood spreading areas have been provided by such significant spatial coincidence between the produced map and the control areas located near Kowsar research station, where the earlier flood spreading projects were successfully performed.     

Modeling the strong ground motion and rupture characteristics of the March 31, 2006, Darb-e-Astane earthquake,Iran, using a hybrid of near-field SH-wave and empirical Green’s function method

We analyze the strong motion accelerograms of the moderate (M _w = 6.1), March 31, 2006, Darb-e-Astane earthquake of western Iran and also those of one of its prominently recorded, large (M _w = 5.1) foreshock and (M _w = 4.9) aftershock. (1) Using derived SH-wave spectral data, we first objectively estimate the parameters W _o\mathit{\Omega} _{\rm o} (long period spectral level), f _c (corner frequency) and Q(f) (frequency dependent, average shear wave quality factor), appropriate for the best-fit Brune ω ^− 2 spectrum of each of these three events. We then perform a non-linear least square analysis of the SH-wave spectral data to provide approximate near-field estimates of the strike, dip, and rake of the causative faults and also the seismic moment, moment magnitude, source size, and average stress drop of these three events. (2) In the next step, we use these approximate values and an empirical Green’s function approach, in an iterative manner, to optimally model the strong ground motion and rupture characteristics of the main event in terms of peak ground acceleration/velocity/displacement and duration of ground shaking and thereby provide improved, more reliable estimates of the causative fault parameters of the main event and its asperities. Our near-field estimates for both the main moderate event and the two smaller events are in good conformity with the corresponding far-field estimates reported by other studies.     

Estimation of coda and shear wave attenuation in the volcanic area in SE Sabalan Mountain,NW Iran

The quality factors of coda and shear waves have been estimated for the SE Sabalan Mountain, geothermal region in northwestern Iran. We have analyzed 65 local earthquakes with magnitude of 2.8 to 6.1 and 2.8 to 5 for shear and coda wave quality factor estimation, respectively. These events were recorded on five stations installed by Building and Housing Research Center Network. Coda normalization and Spectral decay methods have been used to estimate the frequency dependence attenuation relation for shear wave, and single back-scattering method for coda waves. We have observed that the coda normalization method has supplied significantly higher Q _S values as compared to the spectral method. The results show that, in general, Q values are significantly smaller for the entire frequency range as compared to tectonically active areas and are close to the values for volcanic areas.     

Hydraulic conductivity prediction based on grain-size distribution using M5 model tree

The hydraulic conductivity, K_s, is one of the most important hydraulic properties which controls the water and solute movement into the soil. It is measured on soil specimens in the laboratory. On the other hand, sometimes it is obtained by tests carried out in the field by a number of researchers. Therefore, several experimental formulas have developed to predict it. Recently, soft computing tools have been used to evaluate the hydraulic conductivity. However, these tools are not as transparent as empirical formulas. In this study, another soft computing approach, i.e. model trees, have been used for predicting the hydraulic conductivity. The main advantage of model trees is that, unlike the other data learning tools, they are easier to use and represent understandable mathematical rules more clearly. In this paper, a new formula that includes some parameters is derived to estimate the hydraulic conductivity. To develop the new formulas, experimental data sets of hydraulic conductivity were used. A comparison is made between the estimated hydraulic conductivity by this new formula and formulas given by other’s researches.     

Reliability-Based Design of Spatially Random Two-Layered Clayey Slopes

Geotechnical and Geological Engineering - Two-layered cohesive slopes are encountered in geotechnical applications involving embankments, dams, levees, and natural cut slopes. The reliability of...     

The First Light Curve Solutions and Period Study of BQ Ari

Astronomy Letters - The first analysis of the photometric observation in BVR filters of a W UMa type binary system BQ Ari was performed. Light curve analysis was performed using...