An areal distribution of sensors can be used for estimating the direction of incoming waves through beamforming. Beamforming may be implemented as a phase-shifting and stacking of data recorded on the different sensors (i.e., conventional beamforming). Alternatively, beamforming can be applied to cross-correlations between the waveforms on the different sensors. We derive a kernel for beamforming cross-correlated data and call it cross-correlation beamforming (CCBF). We point out that CCBF has slightly better resolution and aliasing characteristics than conventional beamforming. When auto-correlations are added to CCBF, the array response functions are the same as for conventional beamforming. We show numerically that CCBF is more resilient to non-coherent noise. Furthermore, we illustrate that with CCBF individual receiver-pairs can be removed to improve mapping to the slowness domain. An additional flexibility of CCBF is that cross-correlations can be time-windowed prior to beamforming, e.g., to remove the directionality of a scattered wavefield. The observations on synthetic data are confirmed with field data from the SPITS array (Svalbard). Both when beamforming an earthquake arrival and when beamforming ambient noise, CCBF focuses more of the energy to a central beam. Overall, the main advantage of CCBF is noise suppression and its flexibility to remove station pairs that deteriorate the signal-related beampower.

We investigated the shear strain field ahead of a supershear rupture. The strain array data along the sliding fault surfaces were obtained during the large-scale biaxial friction experiments at the National Research Institute for Earth Science and Disaster Resilience. These friction experiments were done using a pair of meter-scale metagabbro rock specimens whose simulated fault area was 1.5 m × 0.1 m. A 2.6-MPa normal stress was applied with loading velocity of 0.1 mm/s. Near-fault strain was measured by 32 two-component semiconductor strain gauges installed at an interval of 50 mm and 10 mm off the fault and recorded at an interval of 1 MHz. Many stick-slip events were observed in the experiments. We chose ten unilateral rupture events that propagated with supershear rupture velocity without preceding foreshocks. Focusing on the rupture front, stress concentration was observed and sharp stress drop occurred immediately inside the ruptured area. The temporal variation of strain array data is converted to the spatial variation of strain assuming a constant rupture velocity. We picked up the peak strain and zero-crossing strain locations to measure the cohesive zone length. By compiling the stick-slip event data, the cohesive zone length is about 50 mm although it scattered among the events. We could not see any systematic variation at the location but some dependence on the rupture velocity. The cohesive zone length decreases as the rupture velocity increases, especially larger than [Formula: see text] times the shear wave velocity. This feature is consistent with the theoretical prediction.

The receiver function (RF) technique is a well-established method to investigate the crustal and upper mantle structures based on three-component seismograms of teleseismic events. In the present study, we propose a modified automatic procedure to determine the back azimuth and polarization angles of a teleseismic event based on the RF technique. The method is tested for the recording of 3 permanent and 3 temporary broadband seismic stations located in the vicinity of Poland. Additionally, the analysis of Rayleigh wave polarization is conducted to show that the new procedure is not sensitive to incorrect seismometer orientation. The synthetic modelling of RF by a modified ray-tracing method for 2.5D models beneath each seismic station down to a depth of 60 km is performed to show the effectiveness of the proposed method in the calculation of RF for a complex structure with dipping layers.

The M 6.2 (Mj 6.8) Nagano (Japan) earthquake of 22 November 2014 produced a 9.3-km long surface rupture zone with a thrust-dominated displacement of up to 1.5 m, which duplicated the pre-existing Kamishiro Fault along the Itoigawa-Shizuoka Tectonic Line (ISTL), the plate-boundary between the Eurasian and North American plates, northern Nagano Prefecture, central Japan. To characterize the activity of the seismogenic fault zone, we conducted a paleoseismic study of the Kamishiro Fault. Field investigations and trench excavations revealed that seven morphogenic paleohistorical earthquakes (E2-E8) prior to the 2014 M 6.2 Nagano earthquake (E1) have occurred on the Kamishiro Fault during the last ca. 6000 years. Three of these events (E2-E4) are well constrained and correspond to historical earthquakes occurring in the last ca. 1200 years. This suggests an average recurrence interval of ca. 300-400 years on the seismogenic fault of the 2014 Kamishiro earthquake in the past 1200 years. The most recent event prior to the 2014 earthquakes (E1) is E2 and the penultimate and antepenultimate faulting events are E3 and E4, respectively. The penultimate faulting event (E3) occurred during the period of AD 1800-1400 and is associated with the 1791 M 6.8 earthquake. The antepenultimate faulting event (E4) is inferred to have occurred during the period of ca. AD 1000-700, likely corresponding to the AD 841 M 6.5 earthquake. The oldest faulting event (E8) in the study area is thought to have occurred during the period of ca. 5600-6000 years. The throw rate during the early Holocene is estimated to be 1.2-3.3 mm/a (average, 2.2 mm/a) with an average amount of characteristic offset of 0.7-1.1 m produced by individual event. When compared with active intraplate faults on Honshu Island, Japan, these slip rates and recurrence interval estimated for morphogenic earthquakes on the Kamishiro Fault along the ISTL appear high and short, respectively. This indicates that present activity on this fault is closely related to seismic faulting along the plate boundary between the Eurasian and North American plates.

This study presents a plan for seismic monitoring of a region around a potential nuclear power plant. Seismic monitoring is needed to evaluate seismic risk. The International Atomic Energy Agency has set guidelines on seismic hazard evaluation and monitoring of such areas. According to these guidelines, we have made a plan for a local network of seismic stations to collect data for seismic source characterization and seismotectonic interpretations, as well as to monitor seismic activity and natural hazards. The detection and location capability of the network were simulated using different station configurations by computing spatial azimuthal coverages and detection threshold magnitudes. Background noise conditions around Pyhäjoki were analyzed by comparing data from different stations. The annual number of microearthquakes that should be detected with a dense local network centered around Pyhäjoki was estimated. The network should be dense enough to fulfill the requirements of azimuthal coverage better than 180° and automatic event location capability down to ML ∼ 0 within a distance of 25 km from the site. A network of 10 stations should be enough to reach these goals. With this setup, the detection threshold magnitudes are estimated to be ML = -0.1 and ML = 0.1 within a radius of 25 and 50 km from Pyhäjoki, respectively. The annual number of earthquakes detected by the network is estimated to be 2 (ML ≥ ∼ -0.1) within 25 km radius and 5 (ML ≥ ∼-0.1 to ∼0.1) within 50 km radius. The location accuracy within 25 km radius is estimated to be 1-2 and 4 km for horizontal coordinates and depth, respectively. Thus, the network is dense enough to map out capable faults with horizontal accuracy of 1-2 km within 25 km radius of the site. The estimation is based on the location accuracies of five existing networks in northern Europe. Local factors, such as seismic noise sources, geology and infrastructure might limit the station configuration and detection and location capability of the network.

A steadily increasing number of citizen seismological stations, often located in populated areas, record a plethora of man-made events. These events are especially of interest, when they are caused by criminal activity or man-made explosions. On 4 December 2021, during an attempted robbery of an automated teller machine (ATM) in Bonn-Röttgen, Germany, the burglars used explosions. The seismic effects of the explosions were recorded with a Raspberry Shake (RS) station at a distance of 580 m from the site. While working on the analysis of this signal, another attack on an ATM on 23 February 2022 in Kürten-Dürscheid was recorded by another RS station, this time at 830-m distance with an instrument that also included an air pressure channel. The seismic signatures of both events indicate similar procedures in both cases whereby a larger explosion was quickly followed a smaller explosion after 21 s and 49 s, respectively. An estimate of the charge weight of the explosions shows that ratios of the strength of the first to second explosion were 21:1 and 9.4:1 in the Röttgen and Dürscheid attacks, respectively.

This paper presents a seismic hazard map for the southern part of Ghana prepared using the probabilistic approach, and seismic hazard assessment results for six cities. The seismic hazard map was prepared for 10% probability of exceedance for peak ground acceleration in 50 years. The input parameters used for the computations of hazard were obtained using data from a catalogue that was compiled and homogenised to moment magnitude (Mw). The catalogue covered a period of over a century (1615-2009). The hazard assessment is based on the Poisson model for earthquake occurrence, and hence, dependent events were identified and removed from the catalogue. The following attenuation relations were adopted and used in this study-Allen (for south and eastern Australia), Silva et al. (for Central and eastern North America), Campbell and Bozorgnia (for worldwide active-shallow-crust regions) and Chiou and Youngs (for worldwide active-shallow-crust regions). Logic-tree formalism was used to account for possible uncertainties associated with the attenuation relationships. OpenQuake software package was used for the hazard calculation. The highest level of seismic hazard is found in the Accra and Tema seismic zones, with estimated peak ground acceleration close to 0.2 g. The level of the seismic hazard in the southern part of Ghana diminishes with distance away from the Accra/Tema region to a value of 0.05 g at a distance of about 140 km.

Ambient seismic noise is caused by a number of sources in specific frequency bands. The quantification of ambient noise makes it possible to evaluate station and network performance. We evaluate noise levels in Norway from the 2013 data set of the Norwegian National Seismic Network as well as two temporary deployments. Apart from the station performance, we studied the geographical and temporal variations, and developed a local noise model for Norway. The microseism peaks related to the ocean are significant in Norway. We, therefore, investigated the relationship between oceanic weather conditions and noise levels. We find a correlation of low-frequency noise (0.125-0.25 Hz) with wave heights up to 900 km offshore. High (2-10 Hz) and intermediate (0.5-5 Hz) frequency noise correlates only up to 450 km offshore with wave heights. From a geographic perspective, stations in southern Norway show lower noise levels for low frequencies due to a larger distance to the dominant noise sources in the North Atlantic. Finally, we studied the influence of high-frequency noise levels on earthquake detectability and found that a noise level increase of 10 dB decreases the detectability by 0.5 magnitude units. This method provides a practical way to consider noise variations in detection maps.

An intraplate earthquake doublet, with 11-min delay between the events, devastated the city of Varzeghan in northwestern Iran on August 11, 2012. The first 6.5 strike-slip earthquake, which occurred after more than 200 years of low seismicity, was followed by an 6.4 oblique thrust event at an epicentral separation of about 6 km. While the first event can be associated with a distinct surface rupture, the absence of a surface fault trace and no clear aftershock signature makes it challenging to identify the fault plane of the second event. We use teleseismic body wave inversion to deduce the slip distribution in the first event. Using both P and SH waves stabilize the inversion and we further constrain the result with the surface rupture extent and the aftershock distribution. The obtained slip pattern shows two distinct slip patches with dissimilar slip directions where aftershocks avoid high-slip areas. Using the estimated slip for the first event, we calculate the induced Coulomb stress change on the nodal planes of the second event and find a preference for higher Coulomb stress on the N-S nodal plane. Assuming a simple slip model for the second event, we estimate the combined Coulomb stress changes from the two events on the focal planes of the largest aftershocks. We find that 90% of the aftershocks show increased Coulomb stress on one of their nodal planes when the N-S plane of the second event is assumed to be the correct fault plane.

The seismic hazard and risk analysis for the onshore Groningen gas field requires information about local soil properties, in particular shear-wave velocity (). A fieldwork campaign was conducted at 18 surface accelerograph stations of the monitoring network. The subsurface in the region consists of unconsolidated sediments and is heterogeneous in composition and properties. A range of different methods was applied to acquire in situ values to a target depth of at least 30 m. The techniques include seismic cone penetration tests (SCPT) with varying source offsets, multichannel analysis of surface waves (MASW) on Rayleigh waves with different processing approaches, microtremor array, cross-hole tomography and suspension P-S logging. The offset SCPT, cross-hole tomography and common midpoint cross-correlation (CMPcc) processing of MASW data all revealed lateral variations on length scales of several to tens of metres in this geological setting. SCPTs resulted in very detailed profiles with depth, but represent point measurements in a heterogeneous environment. The MASW results represent information on a larger spatial scale and smooth some of the heterogeneity encountered at the sites. The combination of MASW and SCPT proved to be a powerful and cost-effective approach in determining representative profiles at the accelerograph station sites. The measured profiles correspond well with the modelled profiles and they significantly enhance the ground motion model derivation. The similarity between the theoretical transfer function from the profile and the observed amplification from vertical array stations is also excellent.

The earthquake (Mw 6 from the SHEEC defined by the MDPs) that occurred in the central part of Slovenia on 14 April, 1895, affected a broad region, causing deaths, injuries, and destruction. This event was much studied but not fully explained; in particular, its causative source model is still debated. The aim of this work is to contribute to the identification of the seismogenic source of this destructive event, calculating peak ground velocity values through the use of different ground motion prediction equations (GMPEs) and computing a series of ground motion scenarios based on the result of an inversion work proposed by Jukić in 2009 and on various fault models in the surroundings of Ljubljana: Vič, Želimlje, Borovnica, Vodice, Ortnek, Mišjedolski, and Dobrepolje faults. The synthetic seismograms, at the basis of our computations, are calculated using the multi-modal summation technique and a kinematic approach for extended sources, with a maximum peak ground velocity value of 1 Hz. The qualitative and quantitative comparison of these simulations with the macroseismic intensity database allows us to discriminate between various sources and configurations. The quantitative validation of the seismic source is done using ad hoc ground motion to intensity conversion equations (GMICEs), expressly calculated for this study. This study allows us to identify the most probable causative source model of this event, contributing to the improvement of the seismotectonic knowledge of this region. The candidate fault that has the lowest values of average differences between observed and calculated intensities and chi-squared is a strike slip fault with a toward-north rupture as the Ortnek fault.

Microzonation is one of the essential tools in seismology to mitigate earthquake damage by estimating the near-surface velocity structure and developing land usage plans and intelligent building design. The number of microzonation studies increased in the last few years as induced seismicity becomes more relevant, even in low-risk areas. While of vital importance, especially in densely populated cities, most of the traditional techniques suffer from different shortcomings. The microzonation technique presented here tries to reduce the existing ambiguity of the inversion results by the combination of single-station six-component (6C) measurements, including three translational and three rotational motions, and more traditional H/V techniques. By applying this new technique to a microzonation study in the downtown area of Munich (Germany) using an iXblue blueSeis-3A rotational motion sensor together with a Nanometrics Trillium Compact seismometer, we were able to estimate Love and Rayleigh wave dispersion curves. These curves together with H/V spectral ratios are then inverted to obtain P- and S-wave velocity profiles of the upper 100 m. In addition, there is a good correlation between the estimated velocity models and borehole-derived lithology, indicating the potential of this single-station microzonation approach.

The single-station microtremor horizontal-to-vertical spectral ratio (MHVSR) method was initially proposed to retrieve the site amplification function and its resonance frequencies produced by unconsolidated sediments overlying high-velocity bedrock. Presently, MHVSR measurements are predominantly conducted to obtain an estimate of the fundamental site frequency at sites where a strong subsurface impedance contrast exists. Of the earthquake site characterization methods presented in this special issue, the MHVSR method is the furthest behind in terms of consensus towards standardized guidelines and commercial use. The greatest challenges to an international standardization of MHVSR acquisition and analysis are (1) the - the underlying composition of the microtremor wavefield is site-dependent, and thus, the appropriate theoretical (forward) model for inversion is still debated; and (2) the - many factors and options are involved in the data acquisition, processing, and interpretation stages. This paper reviews briefly a historical development of the MHVSR technique and the physical basis of an MHVSR (the ). We then summarize recommendations for MHVSR acquisition and analysis (the ). Specific sections address MHVSR interpretation and uncertainty assessment.

Recent estimates of fracture energy in earthquakes show a power-law dependence with slip which can be summarized as ∝ where is a positive real slightly larger than one. For cracks with sliding friction, fracture energy can be equated to : the post-failure integral of the dynamic weakening curve. If the dominant dissipative process in earthquakes is friction, and should be comparable and show a similar scaling with slip. We test this hypothesis by analyzing experiments performed on various cohesive and non-cohesive rock types, under wet and dry conditions, with imposed deformation typical of seismic slip (normal stress of tens of MPa, target slip velocity > 1 m/s and fast accelerations ≈ 6.5 m/s). The resulting fracture energy is similar to the seismological estimates, with and being comparable over most of the slip range. However, appears to saturate after several meters of slip, while in most of the reported earthquake sequences, appears to increase further and surpasses at large magnitudes. We analyze several possible causes of such discrepancy, in particular, additional off-fault damage in large natural earthquakes.

Several recent studies have presented evidence that significant induced earthquakes occurred in a number of oil-producing regions during the early and mid-twentieth century related to either production or wastewater injection. We consider whether the 21 July 1952 M 7.5 Kern County earthquake might have been induced by production in the Wheeler Ridge oil field. The mainshock, which was not preceded by any significant foreshocks, occurred 98 days after the initial production of oil in Eocene strata at depths reaching 3 km, within ~1 km of the White Wolf fault (WWF). Based on this spatial and temporal proximity, we explore a potential causal relationship between the earthquake and oil production. While production would have normally be expected to have reduced pore pressure, inhibiting failure on the WWF, we present an analytical model based on industry stratigraphic data and best estimates of parameters whereby an impermeable splay fault adjacent to the main WWF could plausibly have blocked direct pore pressure effects, allowing the poroelastic stress change associated with production to destabilize the WWF, promoting initial failure. This proof-of-concept model can also account for the 98-day delay between the onset of production and the earthquake. While the earthquake clearly released stored tectonic stress, any initial perturbation on or near a major fault system can trigger a larger rupture. Our proposed mechanism provides an explanation for why significant earthquakes are not commonly induced by production in proximity to major faults.

We study the ground motion simulations based on three finite-source models for the 2007 Mw6.6 Niigata Chuetsu-oki, Japan, earthquake in order to discuss the performance of the input ground motion estimations for the near-field seismic hazard analysis. The three models include a kinematic source inverted from the regional accelerations, a dynamic source on a planar fault with three asperities inferred from the very-near-field ground motion particle motions, and another dynamic source model with conjugate fault segments. The ground motions are calculated for an available 3D geological model using a finite-difference method. For the comparison, we apply a goodness-of-fit score to the ground motion parameters at different stations, including the nearest one that is almost directly above the ruptured fault segments. The dynamic rupture models show good performance. We find that seismologically inferred earthquake asperities on a single fault plane can be expressed with two conjugate segments. The rupture transfer from one segment to another can generate a significant radiation; this could be interpreted as an asperity projected onto a single fault plane. This example illustrates the importance of the fault geometry that has to be taken into account when estimating the very-near-field ground motion.

Earthquakes hit urban centres in Europe infrequently, but occasionally with disastrous effects. Obtaining an unbiased view of seismic hazard (and risk) is therefore very important. In principle, the best way to test probabilistic seismic hazard assessments (PSHAs) is to compare them with observations that are entirely independent of the procedure used to produce PSHA models. Arguably, the most valuable information in this context should be information on long-term hazard, namely maximum intensities (or magnitudes) occurring over time intervals that are at least as long as a seismic cycle. The new observations can provide information of maximum intensity (or magnitude) for long timescale as an input data for PSHA studies as well. Long-term information can be gained from intact stalagmites in natural caves. These formations survived all earthquakes that have occurred over thousands of years, depending on the age of the stalagmite. Their 'survival' requires that the horizontal ground acceleration (HGA) has never exceeded a certain critical value within that time period. Here, we present such a stalagmite-based case study from the Little Carpathians of Slovakia. A specially shaped, intact and vulnerable stalagmite in the Plavecká priepast cave was examined in 2013. This stalagmite is suitable for estimating the upper limit of horizontal peak ground acceleration generated by prehistoric earthquakes. The critical HGA values as a function of time going back into the past determined from the stalagmite that we investigated are presented. For example, at the time of Jókő event (1906), the critical HGA value cannot have been higher than 1 and 1.3 m/s at the time of the assumed Carnuntum event (∼340 AD), and 3000 years ago, it must have been lower than 1.7 m/s. We claimed that the effect of Jókő earthquake (1906) on the location of the Plavecká priepast cave is consistent with the critical HGA value provided by the stalagmite we investigated. The approach used in this study yields significant new constraints on the seismic hazard, as tectonic structures close to Plavecká priepast cave did not generate strong earthquakes in the last few thousand years. The results of this study are highly relevant given that the two capitals, Vienna and Bratislava, are located within 40 and 70 km of the cave, respectively.

The study focuses on the spatial organization of seismicity and the relation between fracture pattern and earthquakes in the Friuli (north-eastern Italy) and western Slovenia seismic regions. The structural setting is characterized by a complex structure resulting from the superposition of several tectonic phases that generated NW-SE trending Dinaric faults and about E-W trending Alpine faults. The upper crust is characterized by lithological and mechanical heterogeneities. The fractal analysis shows that, in general, the seismicity only partially fills a plane. Only in a few cases, the earthquakes distribute on planar structures. The orientation of planes that fit through the hypocentres shows a different disposition at the two depth intervals analysed. The shallower interval (0-10 km) is characterized by planes with highly variable orientations. The spatial seismicity is investigated in the context of a general damage model, represented by the crack density distribution. The results evidence that the seismicity appears mostly located along sharp transition areas from low crack density to higher crack density, i.e., from zones of low damage to zones of intermediate damage. These zones are characterized by high heterogeneity due to the superposition of different tectonic phases and by the maximum interference between Dinaric and Alpine domains. The orientation of the planes fitting the seismicity at 10-20-km depth appears less dispersed, coinciding with the trend of Dinaric sub-vertical faults in the northern and eastern parts of the study area, and with Alpine low-angle faults in the western and southern parts.

Field investigations and analyses of satellite images and aerial photographs reveal that the 2016 7.1 (Mj 7.3) Kumamoto earthquake produced a ∼40-km surface rupture zone striking NE-SW on central Kyushu Island, Japan. Coseismic surface ruptures were characterized by shear faults, extensional cracks, and mole tracks, which mostly occurred along the pre-existing NE-SW-striking Hinagu-Futagawa fault zone in the southwest and central segments, and newly identified faults in the northeast segment. This study shows that (i) the Hinagu-Futagawa fault zone triggered the 2016 Kumamoto earthquake and controlled the spatial distribution of coseismic surface ruptures; (ii) the southwest and central segments were dominated by right-lateral strike-slip movement with a maximum in-site measured displacement of up to 2.5 m, accompanied by a minor vertical component. In contrast, the northeast segment was dominated by normal faulting with a maximum vertical offset of up to 1.75 m with a minor horizontal component that formed graben structures inside Aso caldera; (iii) coseismic rupturing initiated at the jog area between the Hinagu and Futagawa faults, then propagated northeastward into Aso caldera, where it terminated. The 2016 7.1 Kumamoto earthquake therefore offers a rare opportunity to study the relationships between coseismic rupture processes and pre-existing active faults, as well as the seismotectonics of Aso volcano.

The 2014 6.5 (6.1) Ludian earthquake occurred in the eastern Sichuan-Yunnan border region of western China. This earthquake caused much more severe engineering damage than the usual earthquakes with the same magnitude in China. The National Strong Motion Network obtained large set of ground motion recordings during the earthquake. To investigate the engineering interested characteristics of ground motion from Ludian earthquake and compare it with the 7.9 Wenchuan and the 6.6 Lushan earthquakes in western China, studies on the ground motion field, attenuation relationship, distance dependence of significant duration, and site amplification were carried out. Some conclusion is drawn. Specifically, the ground motion field reveals a directional feature, and the distribution characteristics of the two horizontal components are similar. The attenuation relationship for Ludian earthquake is basically consistent with the ground motion prediction equation (GMPE) for western China, except the slight smaller than the GMPE predicted at short periods. The distance dependences of ground motion duration are different in Sichuan and Yunnan regions due to the local physical dispersion and value. The site amplification factors are dominated by linear site response for lower reference ground motion, but the nonlinearity becomes notable for higher reference ground motion. This feature is basically consistent with the empirical model for western China. All the results indicate that the spatial distribution of ground motion, the attenuation characteristics, and the site amplification effect should be considered in characterization of near-field ground motion.