Stress fields and fault reactivation angles of the 2000 western Tottori aftershocks and the 2001 northern Hyogo swarm in southwest Japan

We estimated the stress fields of the aftershocks of the 2000 western Tottori earthquake (M_w 6.6) and the northern Hyogo swarm (max M_w 5.2) by a stress tensor inversion of moment tensor solutions reported from the National Research Institute for Earth Science and Disaster Prevention (Japan). The maximum principal stress direction of the western Tottori sequence was estimated as N107°E with a strike–slip regime. In the northern Hyogo swarm, the orientations of the principal stress directions could not be well constrained by the observed data, but after examining the detailed characteristics of the solution, we obtained a most probable solution of N113°E for the σ₁ direction. These solutions are consistent with the maximum horizontal directions roughly estimated from the strike directions of large earthquakes occurring geographically between these two seismic activities. We measured the angle between each fault–slip direction and maximum principal stress direction to investigate the frictional properties of earthquakes. The distribution of the angles was forward modeled to estimate the coefficient of friction and the stress ratio, assuming uniformly distributed fault orientations. For the western Tottori sequence, a homogeneous stress field with a coefficient of friction less than 0.4 was estimated. A high stress level was also suggested because very little change occurred in the stress field during the mainshock. For the northern Hyogo sequence, the coefficient of friction was estimated to be between 0.5 and 1.0.     

Ground motion estimation during the Kashmir earthquake of 8th October 2005

In this article, analytical methods have been used to estimate ground motion during the 8 October 2005, Kashmir earthquake. Peak ground acceleration (PGA) values at several stations in the epicentral region have been estimated by empirical analytical source mechanism models. As an alternate analysis, PGA estimates have also been obtained using the stochastic finite fault seismological model. The estimated PGAs are compared with that obtained from damage values. A PGA contour map in the near-source region is provided. It is found that very near to the epicenter, PGA would have reached more than 1 g. It is demonstrated that empirical analytical models can be effectively used to estimate ground motion due to rupture of active faults.     

Earthquake source model using strong motion displacement as response of finite elastic media

The strong motion displacement records available during an earthquake can be treated as the response of the earth as the a structural system to unknown forces acting at unknown locations. Thus, if the part of the earth participating in ground motion is modelled as a known finite elastic medium, one can attempt to model the source location and forces generated during an earthquake as an inverse problem in structural dynamics. Based on this analogy, a simple model for the basic earthquake source is proposed. The unknown source is assumed to be a sequence of impulses acting at locations yet to be found. These unknown impulses and their locations are found using the normal mode expansion along with a minimization of mean square error. The medium is assumed to be finite, elastic, homogeneous, layered and horizontal with a specific set of boundary conditions. Detailed results are obtained for Uttarkashi earthquake. The impulse locations exhibit a linear structure closely associated with the causative fault. The results obtained are shown to be in good agreement with reported values. The proposed engineering model is then used to simulate the acceleration time histories at a few recording stations. The earthquake source in terms of a sequence of impulses acting at different locations is applied on a 2D finite elastic medium and acceleration time histories are found using finite element methods. The synthetic accelerations obtained are in close match with the recorded accelerations.     

Strong ground motion estimation in the Sea of Marmara region (Turkey) based on a scenario earthquake

We perform a broadband frequency bedrock strong ground motion simulation in the Marmara Sea region (Turkey), based on several fault rupture scenarios and a source asperity model. The technique combines a deterministic simulation of seismic wave propagation at low frequencies with a semi-stochastic procedure for the high frequencies. To model the high frequencies, we applied a frequency-dependent radiation pattern model, which efficiently removes the effective dependence of the pattern coefficient on the azimuth and take-off angle as the frequency increases. The earthquake scenarios considered consist of the rupture of the closest segments of the North Anatolian Fault System to the city of Istanbul. Our scenario earthquakes involve the rupture of the entire North Anatolian Fault beneath the Sea of Marmara, namely the combined rupture of the Central Marmara Fault and North Boundary Fault segments. We defined three fault rupture scenarios based on the location of the hypocenter, selecting a preferred hypocentral location near a fault bend for each case. We analysed the effect of location of the asperity, within the Central Marmara Fault, on the subsequent ground motion, as well as the influence of anelasticity on the high-frequency attenuation characteristics. The fault and asperity parameters for each scenario were determined from empirical scalings and from results of kinematic and dynamic models of fault rupture. We calculated the resulting time series and spectra for ground motion at Istanbul and evaluated the sensitivity of the predictions to choice of model parameters. The location of the hypocenter is thus shown to be a critical parameter for determining the worst scenario earthquake at Istanbul. We also found that anelasticity has a significant effect on the regional attenuation of peak ground accelerations. Our simulated ground motions result in large values of acceleration response spectra at long periods, which could be critical for building damage at Istanbul during an actual earthquake.     

Estimation of the spatial distribution of ground motion parameters for two recent earthquakes in Japan

A recent development in strong motion instrumentation in Japan provides an opportunity to collect valuable data sets, especially after moderate and large magnitude events. Gathering and modeling these data is a necessity for better understanding of regional ground motion characteristics. Estimations of the spatial distribution of earthquake ground motion plays an important role in early-stage damage assessments for both rescue operations by disaster management agencies as well as damage studies of urban structures. Subsurface geology layers and local soil conditions lead to soil amplification that contributes to the estimated ground motion parameters of the surface. We present a case study of the applicability of the nationally proposed GIS-based soil amplification ratios [J. Soil Dyn. Earthqu. Eng. 19 (2000) 41–53] to the October 6, 2000 Tottori-ken Seibu (western Tottori Prefecture) and the March 24, 2001 Geiyo earthquakes in Japan. First, ground motion values were converted to those at a hypothetical ground base-rock level (outcrop) using an amplification ratio for each 1×1 km area, based on geomorphological and subsurface geology information. Then a Kriging method, assuming an attenuation relationship at the base-rock as a trend component, is applied. Finally, the spatial distribution of ground motion at ground surface is obtained by applying GIS-based amplification factors for the entire region. The correlation between the observed and estimated ground motion values is reasonable for both earthquakes. Thus, the proposed method is applicable in near real-time early-damage assessments and seismic hazard studies in Japan.     

Analysis of soil effects and distribution of damage from the Pyrgos 1993 (Greece) earthquake

Summary On 26 March, 1993, a moderate magnitude earthquake (M _s=5.5) occurred at 3 km epicentral distance from the town of Pyrgos, in Southern Greece, causing extensive damage to masonry houses. To explain the variability of seismic intensity over the town and to propose measures against future seismic activity, a microzonation study was undertaken which combined geological, geophysical and geotechnical investigations, site specific analyses of seismic ground response and detailed recording of structural damage. The analytical predictions of ground response are correlated to soil conditions and then used to identify (micro-)zones of sites with similar seismic response. Furthermore, they are compared to quantitative estimates of damage distribution over the town. It is concluded that the peak ground acceleration, normalized against the input peak seismic acceleration, is a function of the local soil conditions as well as the seismic excitation characteristics. Hence, it cannot be defined uniquely at a site, without reference to the seismic excitation. However, the normalized peak ground velocity and the acceleration response spectra are mainly functions of the soil conditions and can be used as criteria for the practical definition of (micro-)zones. The distribution of damage in various parts of the town is at least partially attributable to local soil effects. The small epicentral distance of the earthquake, connected with the direction of the fault rupture, as well as the quality and techniques of construction, are additional factors that may have influenced the extent and distribution of damage.     

Strong ground motion during the largest aftershock (Mw=5.8) of the 1999 Izmit earthquake, Turkey

We present a preliminary study of strong ground motion during the largest aftershock (M_w 5.8) of the 1999 Izmit earthquake (M_w 7.4), Turkey, at 11:55 on 13 September 1999. The peak ground acceleration observed near the epicentre of this aftershock was in agreement with that predicted by standard empirical prediction equations. Its spectral source parameters of the largest aftershock are also typical for a M_w 5.8 earthquake. At greater epicentral distances, there is an order-of-magnitude in scatter in peak ground acceleration values for this aftershock, which is attributed to site effects. The presence of thick layers of low-velocity sediments caused significant amplification of S-waves in the Avcılar district of Istanbul, at frequencies of 1 Hz, explaining the observed concentration of damage there as a result of the Izmit mainshock.     

First-order rupture features of the 2011 MW 9.0 Tohoku (Japan) earthquake from surface waves

This study analyzed the rupture directivity of the 2011 Tohoku earthquake by using 100-s Rayleigh-wave travel-times, influenced by the finite source, to derive the fault parameters of the earthquake. The results demonstrated that the earthquake exhibited a slow rupture propagation with a rupture velocity of approximately 1.5–2.0 km/s and asymmetric bilateral faulting. The two rupture directions were N60°E and N127°E, with rupture lengths of approximately 276 km and 231 km, respectively. The rupture toward N60°E had a source duration of approximately 183 s, longer than that toward N127°E (approximately 156 s). Overall, the entire source duration of the earthquake faulting lasted approximately 183 s. Regarding historical seismicity in eastern Japan, the 2011 Tohoku earthquake not only ruptured a locked area in which large earthquakes have rarely occurred, but also ruptured the source regions of several historical earthquakes. With the exception of its slow rupture velocity and generation of a tsunami, the rupture features of the 2011 Tohoku earthquake were inconsistent with those of typical tsunami earthquakes.     

Slip distribution and source parameters of the 20 July 2017 Bodrum-Kos earthquake (Mw6.6) from GPS observations

Greek-Turkish boundary near the cities Kos and Bodrum has been shaken on July 20, 2017 by a Mw6.6 earthquake. The mainshock is located offshore and did not generate an on-land surface rupture. Analyzing pre- and post-earthquake continuous/survey-type static GPS observations, we investigated co-seismic surface displacements at 20 sites to characterize source parameters and slip-distribution of the mainshock. Fault plane solutions as well as co-seismic slip distribution have been acquired through the inversion of co-seismic GPS displacements modeling the event as elastic dislocations in a half space. Fault plane solution shows a southward dipping normal-type fault segment extending a depth down to ~12 km, which remains within the brittle upper crust. Results from the distributed slip inversion show that the mainshock activated a ~65 km fault section, which has three high slip patches, namely western, central and eastern patches, where the coseismic slips reach up to 13, 26, and 5 cm, respectively. This slip pattern indicates that the pre-earthquake coupling, which is storing the slip deficit, occurred on these three patches.     

Siderite mineralization of the Gemericum superunit (Western Carpathians, Slovakia): review and a revised genetic model

The Gemericum is a segment of the Variscan orogen subsequently deformed by the Alpine–Carpathian orogeny. The unit contains abundant siderite–sulphide and quartz–antimony veins together with stratabound siderite replacement deposits in limestones and stratiform sulphide mineralization in volcano-sedimentary sequences. The siderite–sulphide veins and siderite replacement deposits of the Gemericum represent one of the largest accumulations of siderite in the world, with about 160 million tonnes of mineable FeCO₃. More than 1200 steeply dipping hydrothermal veins are arranged in a regional tectonic and compositional pattern, reflecting the distribution of regional metamorphic zones. Siderite–sulphide veins are typically contained in low-grade (chlorite zone) sedimentary, volcano-sedimentary or volcanic Lower and Upper Paleozoic rocks. Quartz–antimony veins are hosted by higher-grade units (biotite zone). Siderite–sulphide veins are dominated by early siderite followed by a complex set of stages, including quartz–sulphide (chalcopyrite, tetrahedrite), barite, tourmaline–quartz, and sulphide-remobilization stages. The temporal evolution of these stages is difficult to study because of the widespread and repeated tectonic processes, within-vein replacement and recrystallization. Siderite–sulphide veins show considerable vertical (up to 1200 m) and lateral (up to 15 km) extent, and a thickness typically reaching several metres. Carbonate-replacement siderite deposits of the Gemericum are hosted by a Silurian limestone belt and are similar to stratabound siderite deposits of the Eastern Alps (e.g., Erzberg, Austria).Based on a review of geological, petrological and geochronological data for the Gemericum, and extensive stable and radiogenic isotope data and fluid inclusion data on hydrothermal minerals, the siderite–sulphide veins and siderite replacement deposits are classified as metamorphogenic in a broad sense. The deposits were formed during several stages of regional crustal-scale fluid flow. Isotope (S, C, Sr, Pb) fingerprinting identifies the metamorphosed rock complexes of the Gemericum as a source of most components of hydrothermal fluids. Fluid inclusion and stable isotope data evidence the participation of several contrasting fluid types, and the existence of contrasting P–T conditions during vein evolution. A high-δ¹⁸O, medium- to high-salinity, H₂O-type fluid is the most important component during siderite deposition, whereas H₂O–CO₂-type fluid inclusion containing dense liquid CO₂ and corresponding to minimal pressures between 1 and 3 kbar were found in a younger tourmaline–quartz stage. Younger quartz–ankerite(±siderite)–sulphide stages are characterized by high-salinity (17 to 35 wt.% NaCl equivalent) and low-temperature (T_h=90 to 180 °C) H₂O-type fluids.The vein deposits are interpreted as a result of multistage hydrothermal circulation, with Variscan and Alpine mineralization phases. Based on available indirect data, the most important mineralization phase was related to regional fluid flow during the uplift of a Variscan metamorphic core complex, producing siderite–sulphide (±barite) mineralization, while tourmaline–quartz stage and sulphide remobilization stages are related to Alpine processes. Two phases of vein evolution are evident from two groups of ⁸⁷Sr/⁸⁶Sr isotope ratios of Sr-rich, Rb-poor hydrothermal minerals: 0.71042–0.71541 in older barite and 0.7190–0.7220 in late-stage celestine and strontianite.