Source parameters of March 10 and September 13, 2007, United Arab Emirates earthquakes

On March 10 and September 13, 2007 two earthquakes with moment magnitudes 3.66 and 3.94, respectively, occurred in the eastern part of the United Arab Emirates (UAE). The two events were widely felt in the northern Emirates and Oman and were accompanied by a few aftershocks. Ground motions from these events were well recorded by the broadband stations of Dubai (UAE) and Oman seismological networks and provide an excellent opportunity to study the tectonic process and present day stress field acting in this area. In this study, we report the focal mechanisms of the two main shocks by two methods: first motion polarities and regional waveform moment tensor inversion. Our results indicate nearly pure normal faulting mechanisms with a slight strike slip component. We associated the fault plane trending NNE–SSW with a suggested fault along the extension of the faults bounded Bani Hamid area. The seismicity distribution between two earthquake sequences reveals a noticeable gap that may be a site of a future event. The source parameters (seismic moment, moment magnitude, fault radius, stress drop and displacement across the fault) were also estimated from displacement spectra. The moment magnitudes were very consistent with waveform inversion. The recent deployment of seismic networks in Dubai and Oman reveals tectonic activity in the northern Oman Mountains that was previously unknown. Continued observation and analysis will allow for characterization of seismicity and assessment of seismic hazard in the region.     

Aftershock ground motion characteristics during the 2012 Varzaghan–Ahar doublet events,northwest of Iran

The characteristics of aftershocks can be quite different from those of mainshocks. However, most of the works in the past have generally focused on the aftershocks of strong earthquakes, while some moderate earthquakes can cause aftershocks with magnitudes close to the mainshocks. This paper focuses on the aftershocks characteristics of the 2012 Varzaghan–Ahar Doublet event (Mw?=?6.4, 6.2), northwest of Iran. Numerous aftershocks with magnitudes 3.7–5.5 recorded by the local seismic stations provide useful data to study the spectral characteristics at various periods. The predictive capability of the 2008 Next Generation Attenuation along with two regional models has been examined by several aftershock records obtained during these events. This paper is the first attempt made to compare the predictive capability of modern models based on significant aftershock data obtained for these two moderate events. Results confirmed that the ground motion models which have been derived based on larger-magnitude events should not directly be extrapolated to predict aftershock ground motions with magnitude smaller than 5.5, especially when we use relationships that developed without considering aftershock data. The residual analysis showed that the model of Abrahamson and Silva (Earthq Spectra, 24(1): 67–97, 2008) performed the best overall predictions in this study. However, this model performs poorly in the short period of 0.2 s at distances greater than 50 km which underestimates the spectral values for this case.

     

Neotectonics of Boconó fault, Western Venezuela

Neotectonic morphologic evidence along the Boconó fault (with a mapped length of 500 km) consists of the typical features found along strike-slip faults; offset alluvium and drainage, shutterridges, closed depressions, sag ponds and marshes, fault scarps and trenches, triangular facets, and zones of mylonite and fault gouge. Evidence on fault planes, such as slickensides, suggests a predominant strike-slip displacement, and morphologic evidence suggests that this offset is right-lateral, with a magnitude of 60–100 m during the Holocene, and of several kilometers during the Quaternary. Calculations based on different empirical relationships suggest maximum expected Richter magnitudes of 7.2–7.9 for earthquakes along the fault (using rupture length estimates) and probable intervals of less than 200 years for events of magnitude 8 (using observed total displacement during the Holocene).     

The rupture process during the 1999 Düzce, Turkey, earthquake from joint inversion of teleseismic and strong-motion data

Teleseismic and strong-motion data are inverted to determine the rupture process during the November 1999 Düzce earthquake in NW Turkey. The fault geometry, rise time and rupture velocity are determined from the aftershock distribution and preliminary inversions of the teleseismic data. Joint inversion of the teleseismic and strong-motion data is then carried out for the slip distribution. We obtain the strike 264°, dip 64°, rake −172°, seismic moment 5.0×10¹⁹ N m (M_w 7.1), and average stress drop 7 MPa. This earthquake was characterized by bilateral fault rupture and asymmetric slip distribution. Two asperities (areas of large slip) are identified, the eastern one being 1.5 times larger than the western one. The derived slip distribution is consistent with the aftershock distribution, surface rupture and damage. The point of rupture initiation in this Düzce earthquake coincided with the eastern tip of the aftershock distribution of the August 1999 Izmit earthquake.     

The Erzincan earthquake of 13 March 1992 in Eastern Turkey

The 13 March 1992 Erzincan earthquake, M=6.8, occurred in the eastern half of the Erzincan basin. The largest aftershock took place near Pülümür on 15 March 1992. No clear surface breaks were observed, although teleseismic studies suggested that it was a strike-slip earthquake striking parallel to the North Anatolian fault, with a focus of approximately 10±2 km depth, 30 km rupture length, 95 cm of slip, and a 1.16×10²⁶ dyn.cm seismic moment. The aftershock distribution concentrated at an area of the intersection between the North Anatolian fault and the Ovacik fault. These results indicate that the previously suggested seismic gap along the North Anatolian fault, east of Erzincan, still remains unruptured.     

The 1988 Tennant Creek,northern territory,earthquakes: A synthesis

Three large earthquakes with surface‐wave magnitudes 6.3–6.7 on 22 January 1988 were associated with 32 km of surface faulting on two main scarps 30 km southwest of Tennant Creek in the Northern Territory. These events provide an excellent opportunity to study the mechanics of midplate earthquakes because of the abundance of geological and geophysical data in the area, the proximity of the Warramunga seismic array and the ease of access to the fault zone. The 1988 earthquakes were located in the North Australian Craton in an area that had no history of moderate or large earthquakes before 1986. Additionally, no smaller earthquakes from the fault zone were identified at the Warramunga array, which is situated only 30 km from the nearest scarp, between the 1965 installation of the array and 1986. The main shocks were preceded by a swarm of moderatesized (magnitude 4–5) earthquakes in January 1987 and many smaller aftershocks throughout 1987. Careful relocation of all teleseismically recorded earthquakes from the fault zone shows that the 1987 activity was concentrated in an area only 6 km across in the gap between the two main fault scarps. The main shocks also nucleated in the centre of the fault zone near the 1987 activity. Field observations of scarp morphology indicate that the scarp is divided into three segments, each showing primarily reverse faulting. However, whereas the western and eastern segments show movement of the southern block over the northern, the central scarp segment shows the opposite, with the northern block thrust over the southern block.

Analysis of the first arrival times at Warramunga suggests that the three main shocks were associated with the western, central and eastern scarp segments, respectively. The locations of aftershocks determined using data from temporary seismograph arrays in the epicentral area define three inclined zones of activity that are interpreted as fault planes. In the western and eastern portions of the aftershock zone, these concentrations of activity dip to the south at 45° and 35°, respectively, but in the central section the aftershock zone dips to the north at 55°. Focal mechanisms derived from modelling broadband teleseismic data show thrust and oblique thrust faulting for the three main shocks. The first event ruptured unilaterally up and to the northwest on the westernmost fault segment, while the third main shock ruptured horizontally to the southeast. Modelling of repeat levelling data from the epicentral area requires at least three distinct fault planes, with the eastern and western planes dipping to the south and the central plane dipping to the north. The combination of scarp morphology, aftershock distribution and elevation data makes a strong case for rupture of fault planes in conjugate orientation during the 22 January 1988 Tennant Creek earthquakes. More than 20000 aftershocks have been recorded at Warramunga and activity continues to the present‐day with occasional shocks felt in the town of Tennant Creek and some recent off‐fault aftershocks located directly under the Warramunga seismic array. Stratigraphic relationships exposed in trenches excavated across the scarps suggest that during the Quaternary, a large earthquake ruptured the surface along one segment of the 1988 scarps.     

Evidence for surface rupture associated with the Mw 6.3 L'Aquila earthquake sequence of April 2009 (central Italy)

An earthquake of Mw = 6.3 struck L'Aquila town (central Italy) on 6 April 2009 rupturing an ∼18-km-long SW-dipping normal fault. The aftershock area extended for a length of more than 35 km and included major aftershocks on 7 and 9 April and thousands of minor events. Surface faulting occurred along the SW-dipping Paganica fault with a continuous extent of ∼2.5 km. Ruptures consist of open cracks and vertical dislocations or warps (0.1m maximum throw) with an orientation of N130°–140°. Small triggered slip and shaking effects also took place along nearby synthetic and antithetic normal faults. The observed limited extent and small surface displacement of the Paganica ruptures with respect to the height of the fault scarps and vertical throws of palaeo-earthquakes along faults in the area put the faulting associated with the L'Aquila earthquake in perspective with respect to the maximum expected magnitude and the regional seismic hazard.     

The 2002 Molise earthquake sequence: What can we learn about the tectonics of southern Italy?

On October 31, 2002 a M_L = 5.4 earthquake occurred in southern Italy, at the margin between the Apenninic thrust belt (to the west) and the Adriatic plate (to the east). In this area, neither historical event nor seismogenic fault is reported in the literature. In spite of its moderate magnitude, the earthquake caused severe damage in cities close to the epicenter and 27 people, out of a total of 29 casualties, were killed by the collapse of a primary school in S. Giuliano di Puglia. By inverting broadband regional waveforms, we computed moment tensor solutions for 15 events, as small as M_L = 3.5 (M_w = 3.7). The obtained focal mechanisms show pure strike-slip geometry, mainly with focal planes oriented to NS (sinistral) and EW (dextral). In several solutions focal planes are rotated counterclockwise, in particular for later events, occurring west of the mainshock. From the relocated aftershock distribution, we found that the mainshock ruptured along an EW plane, and the fault mechanisms of some aftershocks were not consistent with the mainshock fault plane. The observed stress field, resulting from the stress tensor inversion, shows a maximum principal stress axis with an east–west trend (N83°W), whereas the minimum stress direction is almost N–S. Considering both the aftershock distribution and moment tensor solutions, it appears that several pre-existing faults were activated rather than a single planar fault associated with the mainshock. The finite fault analysis shows a very simple slip distribution with a slow rupture velocity of 1.1 km/s, that could explain the occurrence of a second mainshock about 30 h after. Finally, we attempt to interpret how the Molise sequence is related to the normal faulting system to the west (along the Apennines) and the dextral strike-slip Mattinata fault to the east.     

Spectral P-wave magnitudes, Aki's ω-square model and source parameters of earthquakes

The earthquake magnitude is a quantity sampling the spectrum of the far-field radiation.With a suit of properly defined magnitudes in a sufficiently broad range of frequencies, the radiated spectrum can be restored and analyzed.A method is proposed for the extraction of stress drop, fault length and seismic moment from magnitudes on a routine basis. Thereby, the theoretical spectrum as predicted by the ω-square model of Aki is utilized.In applying the method to earthquakes which occurred in several parts of Asia over a time-span of 3 years, it is shown that in most cases earthquakes in a given region are characterized by the same stress drop, varying however from region to region. In one region a change of the stress drop with time is found, eventually indicating a change in the state of stress in the particular region during the time interval investigated.     

Fractal analysis of major faults in India on a regional scale

The seismicity of a region is implicit of the causal faulting mechanisms and geodynamic diversity of the subsurface regime nucleating earthquakes of different magnitudes, several of which may be as devastating as ones historically reported in global perspective of tectonic complexity as in the case of India. Fractal analysis using box-counting method for the major fault networks across the country estimates fractal dimension, D_f, values to be varying between 0.88 and 1.36. The fault segments in parts of northwest Himalayas, northeast India and Indo-Gangetic plains, are observed to be associated with higher D_f values implicating high seismicity rates. On the other hand, low D_f values in the peninsular India indicate isolated pattern of the underlying faults. The fractal dimension is observed to be indicative of predominant faulting types — higher values conforming to thrust faulting mechanism while lower to strike slip tectonism.     

Depth variation of seismic source scaling relations: implications for earthquake hazard in southeastern Australia

Advances in earthquake data acquisition and processing techniques have allowed for improved quantification of source parameters for local Australian earthquakes. Until recently, only hypocentral locations and local magnitudes (M_L) had been determined routinely, with little attention given to the inversion of additional source parameters. The present study uses these new source data (e.g. seismic moment, stress drop, source dimensions) to further extend our understanding of seismicity and the continental stress regime of the Australian landmass and its peripheral regions.

Earthquake activity within Australia is typically low, and the proportion of small to large events (i.e. the b value) is also low. It is observed that average stress drops for southeastern Australian earthquakes appear to increase with seismic moment to relatively high levels, up to approximately 10 MPa for M_L 5.0 earthquakes. This is thought to be indicative of high compressive crustal stress, coupled with strong rocks and fault asperities. Furthermore, the data indicates that shallow focus earthquakes (shallower than 6 km) appear to produce lower than average stress drops than deeper earthquakes (between 6 and 20 km) with similar moment.

Recurrence estimates were obtained for a discrete seismogenic zone in southeastern Australia. Decreasing b values with increasing focal depth for this zone indicate that larger earthquakes (with high stress drops) tend to occur deeper in the crust. This may offer an explanation for the apparent increase of stress drop with hypocentral depth. Consequently, earthquake hazard estimates that assume a uniform Gutenburg–Richter distribution with depth (i.e. constant b value) may be too conservative and therefore slightly overestimate seismic hazard for surface sites in southeastern Australia.     

Earthquake activity related to seismic cycles in a model for a heterogeneous strike-slip fault

We investigate the evolution of seismicity within large earthquake cycles in a model of a discrete strike-slip fault in elastic solid. The model dynamics is governed by realistic boundary conditions consisting of constant velocity motion of regions around the fault, static/kinetic friction and dislocation creep along the fault, and 3D elastic stress transfer. The fault consists of brittle parts which fail during earthquakes and undergo small creep deformation between events, and aseismic creep cells which are characterized by high ongoing creep motion. This mixture of brittle and creep cells is found to generate realistic aftershock sequences which follow the modified Omori law and scale with the mainshock size. Furthermore, we find that the distribution of interevent times of the simulated earthquakes is in good agreement with observations. The temporal occurrence, however, is magnitude-dependent; in particular, the small events are clustered in time, whereas the largest earthquakes occur quasiperiodically. Averaging the seismicity before several large earthquakes, we observe an increase of activity and a broadening scaling range of magnitudes when the time of the next mainshock is approached. These results are characteristics of a critical point behavior. The presence of critical point dynamics is further supported by the evolution of the stress field in the model, which is compatible with the observation of accelerating moment release in natural fault systems.     

Active faulting and seismicity along the Siculo–Calabrian Rift Zone (Southern Italy)

Southern Italy is dominated by extensional tectonics that in the Calabrian arc and Eastern Sicily produced the development of the Siculo–Calabrian Rift Zone (SCRZ). This zone is represented by a ≈ 370 km-long fault belt consisting of 10 to 50 km long distinct fault segments which extend both offshore and on land being also responsible of the crustal seismicity of this region. The geological and morphological observations indicate that the active normal faults of the SCRZ are characterized by throw-rates ranging from 0.7 to 3.1 mm/a. They accommodate an almost uniform horizontal extension-rate of about 3.0 mm/a along a WNW–ESE regional extension direction. Based on our field observations and following empirical relationships between magnitude and surface rupture length connections between large crustal earthquakes and distinct fault segments of the SCRZ have been also tentatively tested. Our data indicate moreover that the magnitudes (M) of the historical and instrumental earthquakes are consistent with the estimated values and that the geometry and kinematics of the fault segments and the related different crustal features of the SCRZ control the different seismic behaviours of adjacent portions of the active rift zone.     

Seismic Hazard Analysis for the Bangalore Region

Indian peninsular shield, which was once considered to be seismically stable, is experiencing many earthquakes recently. As part of the national level microzonation programme, Department of Science and Technology, Govt. of India has initiated microzonation of greater Bangalore region. The seismic hazard analysis of Bangalore region is carried out as part of this project. The paper presents the determination of maximum credible earthquake (MCE) and generation of synthetic acceleration time history plot for the Bangalore region. MCE has been determined by considering the regional seismotectonic activity in about 350 km radius around Bangalore city. The seismotectonic map has been prepared by considering the faults, lineaments, shear zones in the area and historic earthquake events of more than 150 events. Shortest distance from the Bangalore to the different sources is measured and then peak ground acceleration (PGA) is calculated for the different source and moment magnitude. Maximum credible earthquake found in terms of moment magnitude is 5.1 with PGA value of 0.146 g at city centre with assuming the hypo central distance of 15.88 km from the focal point. Also, correlations for the fault length with historic earthquake in terms of moment magnitude, yields (taking the rupture fault length as 5% of the total fault length) a PGA value of 0.159 g. Acceleration time history (ground motion) and a response acceleration spectrum for the corresponding magnitude has been generated using synthetic earthquake model considering the regional seismotectonic parameters. The maximum spectral acceleration obtained is 0.332 g for predominant period of 0.06 s. The PGA value and synthetic earthquake ground motion data from the identified vulnerable source using seismotectonic map will be useful for the PGA mapping and microzonation of the area.     

Multiseismic cycle deformation rates from Holocene normal fault scarps on Crete (Greece)

Morphotectonic features reveal recent seismic activity on normal faults on Crete allowing slip rates, palaeoearthquake magnitudes and earthquake recurrence intervals to be inferred. The studied faults show major escarpments (100s m high, 10s km long), separating uplifted Mesozoic rocks from Quaternary deposits. During the Holocene, slip rates outpaced erosion/sedimentation rates and 5–15 m high fresh scarps formed at the base of the major escarpments. Based on our field observations and following empirical relationships between magnitude, surface rupture length and maximum co‐seismic vertical displacement, it is possible to infer some important seismotectonic parameter relative to the investigated faults. Taking into account the possible sources of error and the corresponding uncertainties, our preferred long‐term slip rates range between 0.5 and 1.3 mm yr^?1, the maximum expected magnitudes between 6.3 and 6.6 (or 6.5 and 6.8, considering worst case scenarios), while the mean recurrence intervals range between 260 and c. 840 years. These estimates are in good agreement with the shallow seismicity thus allowing to improve our knowledge on the seismic potential in Crete. Horizontal extension rates could be also inferred from our geological data and compared with the few available GPS data.     

断裂长度与最大位移的关系及其影响因素

断裂最大位移与断裂迹长遵循幂律关系 :D =cLn,但幂指数n的大小有很大的变化范围。为探索幂指数n的大小和断裂机制 ,从已发表的文献中收集了 18组数据 ,这些数据的断裂长度具有 8个数量级的跨度。经相关分析 ,我们得到n值的大小变化于 0 5 5和 1 6 5之间 ,平均值为 1 0 839。由于走滑断裂的最大长度在其倾向方向 ,不宜与倾滑断裂一起统计 ,我们去掉一组走向滑动断裂的数据 ,幂指数平均值为 1 10 6 6。用双回归方法得到的幂指数峰值 (nd)是 1 0～ 1 1。这些结果表明断裂最大位移与断裂迹长应该是非常接近线性关系。这种线性关系可以用Dugdale模型加以解释。该模型认为弹塑性物质拉张裂缝端点的变形是非弹性变形。模型的适用范围是单一岩性 ,一次构造力作用。我们认为n值的大小之所以有很大的变化范围 ,有可能受到断裂迹线长度偏差的影响 ,造成长度偏差的因素包括 :不同的观察平面 ,断裂端点的分辨率 ,断裂连接作用 ,岩石力学性质变化 ,断裂多期活动等。     

Radiative forcing and climate impact resulting from SO2 injections based on a 200,000-year record of Plinian eruptions along the Central American Volcanic Arc

We present for the first time a self-consistent methodology connecting volcanological field data to global climate model estimates for a regional time series of explosive volcanic events. Using the petrologic method, we estimated SO₂ emissions from 36 detected Plinian volcanic eruptions occurring at the Central American Volcanic Arc (CAVA) during the past 200,000 years. Together with simple parametrized relationships collected from past studies, we derive estimates of global maximum volcanic aerosol optical depth (AOD) and radiative forcing (RF) describing the effect of each eruption on radiation reaching the Earth’s surface. In parallel, AOD and RF time series for selected CAVA eruptions are simulated with the global aerosol model MAECHAM5-HAM, which shows a relationship between stratospheric SO₂ injection and maximum global mean AOD that is linear for smaller volcanic eruptions (<5 Mt SO₂) and nonlinear for larger ones (≥5 Mt SO₂) and is qualitatively and quantitatively consistent with the relationship used in the simple parametrized approximation. Potential climate impacts of the selected CAVA eruptions are estimated using an earth system model of intermediate complexity by RF time series derived by (1) directly from the global aerosol model and (2) from the simple parametrized approximation assuming a 12-month exponential decay of global AOD. We find that while the maximum AOD and RF values are consistent between the two methods, their temporal evolutions are significantly different. As a result, simulated global maximum temperature anomalies and the duration of the temperature response depend on which RF time series is used, varying between 2 and 3 K and 60 and 90 years for the largest eruption of the CAVA dataset. Comparing the recurrence time of eruptions, based on the CAVA dataset, with the duration of climate impacts, based on the model results, we conclude that cumulative impacts due to successive eruptions are unlikely. The methodology and results presented here can be used to calculate approximate volcanic forcings and potential climate impacts from sulfur emissions, sulfate aerosol or AOD data for any eruption that injects sulfur into the tropical stratosphere.     

Synthetic seismic intensity for historic earthquakes in the North China Plain: implications for the regional seismic hazard

Multiple approaches are used to study the potential seismic hazard in the North China Craton (NCC, or North China Plain), where approximately 15 % of the Chinese population resides and under which active faults are located. In this study, we develop a new modified Mercalli intensity (MMI) attenuation relationship for the NCC using intensity data from 10 instrumentally recorded events. We then utilize this relationship to infer the magnitude and epicentral location of historic events based on the method proposed by Bakun and Wentworth (Bull Seismol Soc Am 87(6):1502–1521, 1997). In addition, a modified stochastic finite fault model is employed to simulate the strong ground motions caused by these historic events. The simulated peak ground accelerations and velocities are then converted into regional MMI distributions through empirical relationships, and these synthetic MMI maps are compared to field observations. The resultant MMI attenuation versus distance models of the 1976 M _w 7.6 Tangshan event and the 1679 M 8.0 Sanhe-Pinggu event are consistent with the empirical attenuation relationships, and the location and size of the meizoseismal area (>VIII) are consistent with observations. The successful modeling of these historic events indicates that a stochastic finite fault model constrained by the regional MMI attenuation relationship can be used to evaluate a wide range of scenarios based on modern computational simulations. These findings may also provide useful information for the estimation and mitigation of potential seismic hazards in this region.     

Effect of sampling and linkage on fault length and length–displacement relationship

The power-law exponent (n) in the equation: D=cL ⁿ , with D = maximum displacement and L = fault length, would be affected by deviations of fault trace length. (1) Assuming n=1, numerical simulations on the effect of sampling and linkage on fault length and length–displacement relationship are done in this paper. The results show that: (a) uniform relative deviations, which means all faults within a dataset have the same relative deviation, do not affect the value of n; (b) deviations of the fault length due to unresolved fault tip decrease the values of n and the deviations of n increase with the increasing length deviations; (c) fault linkage and observed dimensions either increase or decrease the value of n depending on the distribution of deviations within a dataset; (d) mixed deviations of the fault lengths are either negative or positive and cause the values of n to either decrease or increase; (e) a dataset combined from two or more datasets with different values of c and orders of magnitude also cause the values of n to deviate. (2) Data including 19 datasets and spanning more than eight orders of fault length magnitudes (10⁻²–10⁵ m) collected from the published literature indicate that the values of n range from 0.55 to 1.5, the average value being 1.0813, and the peak value of n _d (double regression) is 1.0–1.1. Based on above results from the simulations and published data, we propose that the relationship between the maximum displacement and fault length in a single tectonic environment with uniform mechanical properties is linear, and the value of n deviated from 1 is mainly caused by the sampling and linkage effects.     

Source parameters of earthquakes in the reservoir-triggered seismic (RTS) zone of Koyna–Warna, Western India

New empirical relations are derived for source parameters of the Koyna–Warna reservoir-triggered seismic zone in Western India using spectral analysis of 38 local earthquakes in the magnitude range M _L 3.5–5.2. The data come from a seismic network operated by the CSIR-National Geophysical Research Institute, India, during March 2005 to April 2012 in this region. The source parameters viz. seismic moment, source radius, corner frequency and stress drop for the various events lie in the range of 10¹³–10¹⁶ Nm, 0.1–0.4 km, 2.9–9.4 Hz and 3–26 MPa, respectively. Linear relationships are obtained among the seismic moment (M ₀), local magnitude (M _L), moment magnitude (M _w), corner frequency (fc) and stress drop (?σ). The stress drops in the Koyna–Warna region are found to increase with magnitude as well as focal depths of earthquakes. Interestingly, accurate depths derived from moment tensor inversion of earthquake waveforms show a strong correlation with the stress drops, seemingly characteristic of the Koyna–Warna region.