Convolutional Neural Network for Seismic Phase Classification, Performance Demonstration over a Local Seismic Network

Over the past two decades, the amount of available seismic data has increased significantly, fueling the need for automatic processing to use the vast amount of information contained in such data sets. Detecting seismicity in temporary aftershock networks is one important example that has become a huge challenge because of the high seismicity rate and dense station coverage. Additionally, the need for highly accurate earthquake locations to distinguish between different competing physical processes during the postseismic period demands even more accurate arrival‐time estimates of seismic phase. Here, we present a convolutional neural network (CNN) for classifying seismic phase onsets for local seismic networks. The CNN is trained on a small dataset for deep‐learning purposes (411 events) detected throughout northern Chile, typical for a temporary aftershock network. In the absence of extensive training data, we demonstrate that a CNN‐based automatic phase picker can still improve performance in classifying seismic phases, which matches or exceeds that of historic methods. The trained network is tested against an optimized short‐term average/long‐term average (STA/LTA) based method (Rietbrock et al., 2012) in classifying phase onsets for a separate dataset of 3878 events throughout the same region. Based on station travel‐time residuals, the CNN outperforms the STA/LTA approach and achieves location residual distribution close to the ones obtained by manual inspection

Seismic anisotropy in the Sumatra subduction zone

An important tool for understanding deformation occurring within a subduction zone is the measurement of seismic anisotropy through observations of shear wave splitting (SWS). In Sumatra, two temporary seismic networks were deployed between December 2007 and February 2009, covering the fore arc between the fore-arc islands to the back arc. We use SKS and local SWS measurements to determine the type, amount, and location of anisotropy. Local SWS measurements from the fore-arc islands exhibit trench-parallel fast directions which can be attributed to shape preferred orientation of cracks/fractures in the overriding sediments. In the Sumatran Fault region, the predominant fast direction is fault/trench parallel, while in the back-arc region it is trench perpendicular. The trench-perpendicular measurements exhibit a positive correlation between delay time and raypath length in the mantle wedge, while the fault-parallel measurements are similar to the fault-parallel fast directions observed for two crustal events at the Sumatran Fault. This suggests that there are two layers of anisotropy: one due to entrained flow within the mantle wedge and a second layer within the overriding crust due to the shear strain caused by the Sumatran Fault. SKS splitting results show a NNW-SSE fast direction with delay times of 0.8–3.0 s. The fast directions are approximately parallel to the absolute plate motion of the subducting Indo-Australian Plate. The small delay times exhibited by the local SWS (0.05–0.45 s), in combination with the large SKS delay times, suggest that the anisotropy generating the teleseismic SWS is dominated by entrained flow in the asthenosphere below the slab

The Devastating 2022 M6.2 Afghanistan Earthquake: Challenges, Processes, and Implications

On June 21st, a Mw6.2 earthquake struck the Afghan-Pakistan-border-region, situated within the India-Asia collision. Thousand thirty-nine deaths were reported, making the earthquake the deadliest of 2022. We investigate the event\u27s rupture processes by combining seismological and geodetic observations, aiming to understand what made it that fatal. Our Interferometric Synthetic Aperture Radar-constrained slip-model and regional moment-tensor inversion, confirmed through field observations, reveal a sinistral rupture with maximum slip of 1.8 m at 5 km depth on a N20°E striking, sub-vertical fault. We suggest that not only external factors (event-time, building stock) but fault-specific factors made the event excessively destructive. Surface rupture was favored by the rock foliation, coinciding with the fault strike. The distribution of Peak-Ground-Velocity was governed by the sub-vertical fault. Maximum slip was large compared to other events globally and might have resulted in peak-frequencies coinciding with resonance-frequencies of the local buildings and demonstrates the devastating impact of moderate-size earthquakes

Imaging slab-transported fluids and their deep dehydration from seismic velocity tomography in the Lesser Antilles subduction zone

Volatiles play a pivotal role in subduction zone evolution, yet their pathways remain poorly constrained. Studying the Lesser Antilles subduction zone can yield new constraints, where old oceanic lithosphere formed by slow-spreading subducts slowly. Here we use local earthquakes recorded by the temporary VoiLA (Volatile recycling in the Lesser Antilles) deployment of ocean-bottom seismometers in the fore- and back-arc to characterize the 3-D seismic structure of the north-central Lesser Antilles subduction zone. Along the slab top, mapped based on seismicity, we find low Vp extending to 130–150 km depth, deeper than expected for magmatic oceanic crust. The slab\u27s most prominent, elevated Vp/Vs anomalies are beneath the fore- and back-arc offshore Guadeloupe and Dominica, where two subducted fracture zones lie with the obliquely subducting boundary between Proto-Caribbean and Equatorial Atlantic lithosphere. These structures, therefore, enhance hydration of the oceanic lithosphere as it forms and evolves and the subsequent dehydration of mantle serpentinite when subducted. Above the slab, we image the asthenosphere wedge as a high Vp/Vs and moderate Vp feature, indicating slab-dehydrated fluids rising through the overlying cold boundary layer that might induce melting further to the west. Our results provide new evidence for the impact of spatially-variable oceanic plate formation processes on slab dehydration and mantle wedge volatile transfer that ultimately impact volcanic processes at the surface, such as the relatively high magmatic output observed on the north-central islands in the Lesser Antilles

The devastating 2022 M6.2 Afghanistan earthquake: Challenges, processes, and implications

On June 21st, a Mw6.2 earthquake struck the Afghan-Pakistan-border-region, situated within the India-Asia collision. Thousand thirty-nine deaths were reported, making the earthquake the deadliest of 2022. We investigate the event's rupture processes by combining seismological and geodetic observations, aiming to understand what made it that fatal. Our Interferometric Synthetic Aperture Radar-constrained slip-model and regional moment-tensor inversion, confirmed through field observations, reveal a sinistral rupture with maximum slip of 1.8 m at 5 km depth on a N20°E striking, sub-vertical fault. We suggest that not only external factors (event-time, building stock) but fault-specific factors made the event excessively destructive. Surface rupture was favored by the rock foliation, coinciding with the fault strike. The distribution of Peak-Ground-Velocity was governed by the sub-vertical fault. Maximum slip was large compared to other events globally and might have resulted in peak-frequencies coinciding with resonance-frequencies of the local buildings and demonstrates the devastating impact of moderate-size earthquakes

Triggered crustal earthquake swarm across subduction segment boundary after the 2016 Pedernales, Ecuador megathrust earthquake

Megathrust ruptures and the ensuing postseismic deformation cause stress changes that may induce seismicity on upper plate crustal faults far from the coseismic rupture area. In this study, we analyze seismic swarms that occurred in the north Ecuador area of Esmeraldas, beginning two months after the 2016 M$_{w}$ 7.8 Pedernales, Ecuador megathrust earthquake. The Esmeraldas region is 70 km from the Pedernales rupture area in a separate segment of the subduction zone. We characterize the Esmeraldas sequence, relocating the events using manual arrival time picks and a local a-priori 3D velocity model. The earthquake locations from the Esmeraldas sequence outline an upper plate fault or shear zone. The sequence contains one major swarm and several smaller swarms. Moment tensor solutions of several events include normal and strike-slip motion and non-double-couple components. During the main swarm, earthquake hypocenters increase in distance from the first event over time, at a rate of a few hundred meters per day, consistent with fluid diffusion. Events with similar waveforms occur within the sequence, and a transient is seen in time series of nearby GPS stations concurrent with the seismicity. The events with similar waveforms and the transient in GPS time series suggest that slow aseismic slip took place along a crustal normal fault during the sequence. Coulomb stress calculations show a positive Coulomb stress change in the Esmeraldas region, consistent with seismicity being triggered by the Pedernales mainshock and large aftershocks. The characteristics of the seismicity indicate that postseismic deformation involving fluid flow and slow slip activated upper plate faults in the Esmeraldas area. These findings suggest the need for further investigation into the seismic hazard potential of shallow upper plate faults and the potential for megathrust earthquakes to trigger slow-slip and shallow seismicity across separate segments of subduction zones

Do distribution volumes and clearances relate to tissue volumes and blood flows? A computer simulation

BACKGROUND: Kinetics of inhaled agents are often described by physiological models. However, many pharmacokinetic concepts, such as context-sensitive half-times, have been developed for drugs described by classical compartmental models. We derived classical compartmental models that describe the course of the alveolar concentrations (F(A)) generated by the physiological uptake and distribution models used by the Gas Man(® )program, and describe how distribution volumes and clearances relate to tissue volumes and blood flows. METHODS: Gas Man(® )was used to generate F(A )vs. time curves during the wash-in and wash-out period of 115 min each with a high fresh gas flow (8 L.min(-1)), a constant alveolar minute ventilation (4 L.min(-1)), and a constant inspired concentration (F(I)) of halothane (0.75%), isoflurane (1.15%), sevoflurane (2%), or desflurane (6%). With each of these F(I), simulations were ran for a 70 kg patient with 5 different cardiac outputs (CO) (2, 3, 5, 8 and 10 L.min(-1)) and for 5 patients with different weights (40, 55, 70, 85, and 100 kg) but the same CO (5 L.min(-1)). Two and three compartmental models were fitted to F(A )of the individual 9 runs using NONMEM. After testing for parsimony, goodness of fit was evaluated using median prediction error (MDPE) and median absolute prediction error (MDAPE). The model was tested prospectively for a virtual 62 kg patient with a cardiac output of 4.5 L.min(-1 )for three different durations (wash-in and wash-out period of 10, 60, and 180 min each) with an F(I )of 1.5% halothane, 1.5% isoflurane, sevoflurane 4%, or desflurane 12%. RESULTS: A three-compartment model fitted the data best (MDPE = 0% and MDAPE ≤ 0.074%) and performed equally well when tested prospectively (MDPE ≤ 0.51% and MDAPE ≤ 1.51%). The relationship between CO and body weight and the distribution volumes and clearances is complex. CONCLUSION: The kinetics of anesthetic gases can be adequately described e by a mammilary compartmental model. Therefore, concepts that are traditionally thought of as being applicable to the kinetics of intravenous agents can be equally well applied to anesthetic gases. Distribution volumes and clearances cannot be equated to tissue volumes and blood flows respectively

Selection of Medical Diagnostic Codes for Analysis of Electronic Patient Records. Application to Stroke in a Primary Care Database

BACKGROUND: Electronic patient records from primary care databases are increasingly used in public health and health services research but methods used to identify cases with disease are not well described. This study aimed to evaluate the relevance of different codes for the identification of acute stroke in a primary care database, and to evaluate trends in the use of different codes over time.METHODS: Data were obtained from the General Practice Research Database from 1997 to 2006. All subjects had a minimum of 24 months of up-to-standard record before the first recorded stroke diagnosis. Initially, we identified stroke cases using a supplemented version of the set of codes for prevalent stroke used by the Office for National Statistics in Key health statistics from general practice 1998 (ONS codes). The ONS codes were then independently reviewed by four raters and a restricted set of 121 codes for 'acute stroke' was identified but the kappa statistic was low at 0.23.RESULTS: Initial extraction of data using the ONS codes gave 48,239 cases of stroke from 1997 to 2006. Application of the restricted set of codes reduced this to 39,424 cases. There were 2,288 cases whose index medical codes were for 'stroke annual review' and 3,112 for 'stroke monitoring'. The frequency of stroke review and monitoring codes as index codes increased from 9 per year in 1997 to 1,612 in 2004, 1,530 in 2005 and 1,424 in 2006. The one year mortality of cases with the restricted set of codes was 29.1% but for 'stroke annual review,' 4.6% and for 'stroke monitoring codes', 5.7%.CONCLUSION: In the analysis of electronic patient records, different medical codes for a single condition may have varying clinical and prognostic significance; utilisation of different medical codes may change over time; researchers with differing clinical or epidemiological experience may have differing interpretations of the relevance of particular codes. There is a need for greater transparency in the selection of sets of codes for different conditions, for the reporting of sensitivity analyses using different sets of codes, as well as sharing of code sets among researchers