ALBERT

Description: We investigated the infrasound signals from seismic ground motions induced by North Korea's underground nuclear explosions, including the recent third explosion on 2013 February 12. For the third explosion, the epicentral infrasound signals were detected not only by three infrasound network stations (KSGAR, ULDAR and YAGAR) in South Korea but also by two nearby International Monitoring System infrasound stations, IS45 and IS30. The detectability of the signals was limited at stations located on the relatively east side of the epicentre, with large azimuth deviations due to very favourable atmospheric conditions for eastward propagation at stratospheric height in 2013. The stratospheric wind direction was the reverse of that when the second explosion was conducted in 2009 May. The source location of the epicentral infrasound with wave parameters determined at the multiple stations has an offset by about 16.6 km from the reference seismic location. It was possible to determine the infrasonic location with moderate accuracy by the correction of the azimuth deviation due to the eastward winds in the stratosphere. In addition to the epicentral infrasonic signals, diffracted infrasound signals were observed from the second underground nuclear explosion in 2009. The exceptional detectability of the diffracted infrasound was a consequence of the temporal formation of a thin atmospheric inversion layer over the ocean surface when the event occurred.

Description: Mainland Portugal, on the southwestern edge of the European continent, is located directly north of the boundary between the Eurasian and Nubian plates. It lies in a region of slow lithospheric deformation (〈5 mm yr –1 ), which has generated some of the largest earthquakes in Europe, both intraplate (mainland) and interplate (offshore). Some offshore earthquakes are nucleated on old and cold lithospheric mantle, at depths down to 60 km. The seismicity of mainland Portugal and its adjacent offshore has been repeatedly classified as diffuse. In this paper, we analyse the instrumental earthquake catalogue for western Iberia, which covers the period between 1961 and 2013. Between 2010 and 2012, the catalogue was enriched with data from dense broad-band deployments. We show that although the plate boundary south of Portugal is diffuse, in that deformation is accommodated along several distributed faults rather than along one long linear plate boundary, the seismicity itself is not diffuse. Rather, when located using high-quality data, earthquakes collapse into well-defined clusters and lineations. We identify and characterize the most outstanding clusters and lineations of epicentres and correlate them with geophysical and tectonic features (historical seismicity, topography, geologically mapped faults, Moho depth, free-air gravity, magnetic anomalies and geotectonic units). Both onshore and offshore, clusters and lineations of earthquakes are aligned preferentially NNE–SSW and WNW–ESE. Cumulative seismic moment and epicentre density decrease from south to north, with increasing distance from the plate boundary. Only few earthquake lineations coincide with geologically mapped faults. Clusters and lineations that do not match geologically mapped faults may correspond to previously unmapped faults (e.g. blind faults), rheological boundaries or distributed fracturing inside blocks that are more brittle and therefore break more easily than neighbour blocks. The seismicity map of western Iberia presented in this article opens important questions concerning the regional seismotectonics. This work shows that the study of low-magnitude earthquakes using dense seismic deployments is a powerful tool to study lithospheric deformation in slowly deforming regions, such as western Iberia, where high-magnitude earthquakes occur with long recurrence intervals.

Description: The Piton de la Fournaise basaltic volcano, on La Réunion Island in the western Indian Ocean, is one of the most active volcanoes in the world. This volcano is classically considered as the surface expression of an upwelling mantle plume and its activity is continuously monitored, providing detailed information on its superficial dynamics and on the edifice structure. Deeper crustal and upper mantle structure under La Réunion Island is surprisingly poorly constrained, motivating this study. We used receiver function techniques to determine a shear wave velocity profile through the crust and uppermost mantle beneath La Réunion, but also at other seismic stations located on the hotspot track, to investigate the plume and lithosphere interaction and its evolution through time. Receiver functions (RFs) were computed at permanent broad-band seismic stations from the GEOSCOPE network (on La Réunion and Rodrigues), at IRIS stations MRIV and DGAR installed on Mauritius and Diego Garcia islands, and at the GEOFON stations KAAM and HMDM on the Maldives. We performed non-linear inversions of RFs through modelling of P -to- S conversions at various crustal and upper mantle interfaces. Joint inversion of RF and surface wave dispersion data suggests a much deeper Mohorovičić discontinuity (Moho) beneath Mauritius (~21 km) compared to La Réunion (~12 km). A magmatic underplated body may be present under La Réunion as a thin layer (≤3 km thick), as suggested by a previous seismic refraction study, and as a much thicker layer beneath other stations located on the hotspot track, suggesting that underplating is an important process resulting from the plume–lithosphere interaction. We find evidence for a strikingly low velocity layer starting at about 33 km depth beneath La Réunion that we interpret as a zone of partial melt beneath the active volcano. We finally observe low velocities below 70 km beneath La Réunion and below 50 km beneath Mauritius that could represent the base of the oceanic lithosphere.

Description: We analyse daily cross-correlation computed from continuous records by permanent stations operating in vicinity of the Klyuchevskoy group of volcanoes (Kamchatka). Seismic waves generated by volcanic tremors are clearly seen on the cross-correlations between some pairs of stations as strong signals at frequencies between 0.2 and 2 Hz and with traveltimes typically shorter than those corresponding to interstation propagation. First, we develop a 2-D source-scanning algorithm based on summation of the envelops of cross-correlations to detect seismic tremors and to determine locations from which the strong seismic energy is continuously emitted. In an alternative approach, we explore the distinctive character of the cross-correlation waveforms corresponding to tremors emitted by different volcanoes and develop a phase-matching method for detecting volcanic tremors. Application of these methods allows us to detect and to distinguish tremors generated by the Klyuchevskoy and the Tolbachik, volcanoes and to monitor evolution of their intensity in time.

Description: Powerful subduction zone earthquakes rupture thousands of square kilometres along continental margins but at certain locations earthquake rupture terminates. To date, detailed knowledge of the parameters that govern seismic rupture and aftershocks is still incomplete. On 2015 September 16, the M w 8.3 Illapel earthquake ruptured a 200 km long stretch of the Central Chilean subduction zone, triggering a tsunami and causing significant damage. Here, we analyse the temporal and spatial pattern of the coseismic rupture and aftershocks in relation to the tectonic setting in the earthquake area. Aftershocks cluster around the area of maximum coseismic slip, in particular in lateral and downdip direction. During the first 24 hr after the main shock, aftershocks migrated in both lateral directions with velocities of approximately 2.5 and 5 km hr –1 . At the southern rupture boundary, aftershocks cluster around individual subducted seamounts that are related to the downthrusting Juan Fernández Ridge. In the northern part of the rupture area, aftershocks separate into an upper cluster (above 25 km depth) and a lower cluster (below 35 km depth). This dual seismic–aseismic transition in downdip direction is also observed in the interseismic period suggesting that it may represent a persistent feature for the Central Chilean subduction zone.

Description: We propose a procedure for uncertainty quantification in Probabilistic Tsunami Hazard Analysis (PTHA), with a special emphasis on the uncertainty related to statistical modelling of the earthquake source in Seismic PTHA (SPTHA), and on the separate treatment of subduction and crustal earthquakes (treated as background seismicity). An event tree approach and ensemble modelling are used in spite of more classical approaches, such as the hazard integral and the logic tree. This procedure consists of four steps: (1) exploration of aleatory uncertainty through an event tree, with alternative implementations for exploring epistemic uncertainty; (2) numerical computation of tsunami generation and propagation up to a given offshore isobath; (3) (optional) site-specific quantification of inundation; (4) simultaneous quantification of aleatory and epistemic uncertainty through ensemble modelling. The proposed procedure is general and independent of the kind of tsunami source considered; however, we implement step 1, the event tree, specifically for SPTHA, focusing on seismic source uncertainty. To exemplify the procedure, we develop a case study considering seismic sources in the Ionian Sea (central-eastern Mediterranean Sea), using the coasts of Southern Italy as a target zone. The results show that an efficient and complete quantification of all the uncertainties is feasible even when treating a large number of potential sources and a large set of alternative model formulations. We also find that (i) treating separately subduction and background (crustal) earthquakes allows for optimal use of available information and for avoiding significant biases; (ii) both subduction interface and crustal faults contribute to the SPTHA, with different proportions that depend on source-target position and tsunami intensity; (iii) the proposed framework allows sensitivity and deaggregation analyses, demonstrating the applicability of the method for operational assessments.

Description: In order to improve our understanding of hazardous underground cavities, the development and collapse of a ~200 m wide salt solution mining cavity was seismically monitored in the Lorraine basin in northeastern France. The microseismic events show a swarm-like behaviour, with clustering sequences lasting from seconds to days, and distinct spatiotemporal migration. Observed microseismic signals are interpreted as the result of detachment and block breakage processes occurring at the cavity roof. Body wave amplitude patterns indicated the presence of relatively stable source mechanisms, either associated with dip-slip and/or tensile faulting. Signal overlaps during swarm activity due to short interevent times, the high-frequency geophone recordings and the limited network station coverage often limit the application of classical source analysis techniques. To overcome these shortcomings, we investigated the source mechanisms through different procedures including modelling of observed and synthetic waveforms and amplitude spectra of some well-located events, as well as modelling of peak-to-peak amplitude ratios for the majority of the detected events. We extended the latter approach to infer the average source mechanism of many swarming events at once, using multiple events recorded at a single three component station. This methodology is applied here for the first time and represents a useful tool for source studies of seismic swarms and seismicity clusters. The results obtained with different methods are consistent and indicate that the source mechanisms for at least 50 per cent of the microseismic events are remarkably stable, with a predominant thrust faulting regime with faults similarly oriented, striking NW–SE and dipping around 35°–55°. This dominance of consistent source mechanisms might be related to the presence of a preferential direction of pre-existing crack or fault structures. As an interesting byproduct, we demonstrate, for the first time directly on seismic data, that the source radiation pattern significantly controls the detection capability of a seismic station and network.

Description: The proposed study deals with synthesis of seismograms by BIEM (boundary integral equation method) taking into account all three base components-seismic source, wave path and local region of interest. Consider a laterally inhomogeneous geological profile situated in a half-plane with non-parallel layers. Seismic load is time-harmonic or transient in time. It is presented by incident SH wave or wave radiating from an embedded line seismic source. Two types of lateral inhomogeneities with arbitrary shape and located in the inhomogeneous half-plane are considered: (i) free-surface relief as a canyon or a hill; (ii) alluvial basin with properties different from those of the layered half-plane. The computational tool is BIEM based on the frequency-dependent elastodynamic fundamental solutions. A relation between displacements and tractions along the free surface and arbitrary interface of the soil stratum is derived, which is applicable for arbitrary geometry of the interfaces between soil layers. Validation and convergence study is presented. All simulations reveal the sensitivity of the synthetic seismic signals on the type and characteristics of the seismic time-harmonic or transient load, on the wave path inhomogeneity and on the specific geotechnical properties of the local geological region.

Description: Fractured rocks are known to exhibit seismic anisotropy and shear wave splitting (SWS). SWS is commonly used for fractured rock characterization and has been shown to be sensitive to fluid type. The presence of partial liquid/gas saturation is also known to affect the elastic properties of rocks. The combined effect of both fractures and partial liquid/gas saturation is still unknown. Using synthetic, silica-cemented sandstones with aligned penny-shaped voids, we conducted laboratory ultrasonic experiments to investigate the effect fractures aligned at an oblique angle to wave propagation would have on SWS under partial liquid/gas saturation conditions. The result for the fractured rock shows a saturation dependence which can be explained by combining a fractured rock model and a partial saturation model. At high to full water saturation values, SWS decreases as a result of the fluid bulk modulus effect on the quasi-shear wave. This bulk modulus effect is frequency dependent as a result of wave-induced fluid flow mechanisms, which would in turn lead to frequency dependent SWS. This result suggests the possible use of SWS for discriminating between full liquid saturation and partial liquid/gas saturation.

Description: The depth of an earthquake is difficult to estimate because of the trade-off between depth and origin time estimations, and because it can be biased by lateral Earth heterogeneities. To face this challenge, we have developed a new, blind and fully automatic teleseismic depth analysis. The results of this new method do not depend on epistemic uncertainties due to depth-phase picking and identification. The method consists of a modification of the cepstral analysis from Letort et al. and Bonner et al. , which aims to detect surface reflected (pP, sP) waves in a signal at teleseismic distances (30°–90°) through the study of the spectral holes in the shape of the signal spectrum. The ability of our automatic method to improve depth estimations is shown by relocation of the recent moderate seismicity of the Guerrero subduction area (Mexico). We have therefore estimated the depth of 152 events using teleseismic data from the IRIS stations and arrays. One advantage of this method is that it can be applied for single stations (from IRIS) as well as for classical arrays. In the Guerrero area, our new cepstral analysis efficiently clusters event locations and provides an improved view of the geometry of the subduction. Moreover, we have also validated our method through relocation of the same events using the new International Seismological Centre (ISC)-locator algorithm, as well as comparing our cepstral depths with the available Harvard–Centroid Moment Tensor (CMT) solutions and the three available ground thrust (GT5) events (where lateral localization is assumed to be well constrained with uncertainty 〈5 km) for this area. These comparisons indicate an overestimation of focal depths in the ISC catalogue for deeper parts of the subduction, and they show a systematic bias between the estimated cepstral depths and the ISC-locator depths. Using information from the CMT catalogue relating to the predominant focal mechanism for this area, this bias can be explained as a misidentification of sP phases by pP phases, which shows the greater interest for the use of this new automatic cepstral analysis, as it is less sensitive to phase identification.