PickBlue: Seismic Phase Picking for Ocean Bottom Seismometers With Deep Learning

Detecting phase arrivals and pinpointing the arrival times of seismic phases in seismograms is crucial for many seismological analysis workflows. For land station data, machine learning methods have already found widespread adoption. However, deep learning approaches are not yet commonly applied to ocean bottom data due to a lack of appropriate training data and models. Here, we compiled an extensive and labeled ocean bottom seismometer (OBS) data set from 15 deployments in different tectonic settings, comprising ∼90,000 P and ∼63,000 S manual picks from 13,190 events and 355 stations. We propose PickBlue, an adaptation of the two popular deep learning networks EQTransformer and PhaseNet. PickBlue joint processes three seismometer recordings in conjunction with a hydrophone component and is trained with the waveforms in the new database. The performance is enhanced by employing transfer learning, where initial weights are derived from models trained with land earthquake data. PickBlue significantly outperforms neural networks trained with land stations and models trained without hydrophone data. The model achieves a mean absolute deviation of 0.05 s for P-waves and 0.12 s for S-waves, and we apply the picker on the Hikurangi Ocean Bottom Tremor and Slow Slip OBS deployment offshore New Zealand. We integrate our data set and trained models into SeisBench to enable an easy and direct application in future deployments

Deep lithospheric structures along the southern central Chile Margin from wide-angle P-wave modellilng

Crustal- and upper-mantle structures of the subduction zone in south central Chile, between 42 degrees S and 46 degrees S, are determined from seismic wide-angle reflection and refraction data, using the seismic ray tracing method to calculate minimum parameter models. Three profiles along differently aged segments of the subducting Nazca Plate were analysed in order to study subduction zone structure dependencies related to the age, that is, thermal state, of the incoming plate. The age of the oceanic crust at the trench ranges from 3 Ma on the southernmost profile, immediately north of the Chile triple junction, to 6.5 Ma old about 100 km to the north, and to 14.5 Ma old another 200 km further north, off the Island of Chiloe. Remarkable similarities appear in the structures of both the incoming as well as the overriding plate. The oceanic Nazca Plate is around 5 km thick, with a slightly increasing thickness northward, reflecting temperature changes at the time of crustal generation. The trench basin is about 2 km thick except in the south where the Chile Ridge is close to the deformation front and only a small, 800-m-thick trench infill could develop. In south central Chile, typically three quarters (1.5 km) of the trench sediments subduct below the decollement in the subduction channel. To the north and south of the study area, only about one quarter to one third of the sediments subducts, the rest is accreted above. Similarities in the overriding plate are the width of the active accretionary prism, 35-50 km, and a strong lateral crustal velocity gradient zone about 75-80 km landward from the deformation front, where landward upper-crustal velocities of over 5.0-5.4 km s<SU-1</SU decrease seaward to around 4.5 km s<SU-1</SU within about 10 km, which possibly represents a palaeo-backstop. This zone is also accompanied by strong intraplate seismicity. Differences in the subduction zone structures exist in the outer rise region, where the northern profile exhibits a clear bulge of uplifted oceanic lithosphere prior to subduction whereas the younger structures have a less developed outer rise. This plate bending is accompanied by strongly reduced rock velocities on the northern profile due to fracturing and possible hydration of the crust and upper mantle. The southern profiles do not exhibit such a strong alteration of the lithosphere, although this effect may be counteracted by plate cooling effects, which are reflected in increasing rock velocities away from the spreading centre. Overall there appears little influence of incoming plate age on the subduction zone structure which may explain why the M-w = 9.5 great Chile earthquake from 1960 ruptured through all these differing age segments. The rupture area, however, appears to coincide with a relatively thick subduction channel

Distribution of Temperature and Strength in the Central Andean Lithosphere and Its Relationship to Seismicity and Active Deformation

We present three-dimensional (3D) models of the present-day steady-state conductive thermal field and strength distribution in the lithosphere beneath the Central Andes. Our primary objective was to investigate the influence that the structure of the Central Andean lithosphere has on its thermal and rheological state, and the relationship between the latter and the active deformation in the region. We used our previous data-driven and gravity-constrained 3D density model as starting point for the calculations. We first assigned lithology-derived thermal and rheological properties to the different divisions of the density model and defined temperature boundary conditions. We then calculated the 3D steady-state conductive thermal field and the maximum differential stresses for both brittle and ductile behaviors. We find that the thickness and composition of the crust are the main factors affecting the modeled thermal field, and consequently also the strength distribution. The orogen is characterized by a thick felsic crust with elevated temperatures and a low integrated strength, whereas the foreland and forearc are underlain by a more mafic and thinner crust with lower temperatures and a higher integrated strength. We find that most of the intraplate deformation coincides spatially with the steepest strength gradients and suggest that the high potential energy of the orogen together with the presence of rheological lateral heterogeneities produce high compressional stresses and strong strain localization along the margins of the orogen. We interpret earthquakes within the modeled ductile field to be related to the weakening effect of long-lived faults and/or the presence of seismic asperities.Fil: Ibarra, Federico. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Geociencias Básicas, Aplicadas y Ambientales de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Geociencias Básicas, Aplicadas y Ambientales de Buenos Aires; ArgentinaFil: Prezzi, Claudia Beatriz. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Geociencias Básicas, Aplicadas y Ambientales de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Geociencias Básicas, Aplicadas y Ambientales de Buenos Aires; ArgentinaFil: Bott, Judith. German Research Centre for Geosciences; AlemaniaFil: Scheck Wenderoth, Magdalena. German Research Centre for Geosciences; AlemaniaFil: Strecker, Manfred. Universitat Potsdam; Alemani

Skeletal muscle properties and fatigue resistance in relation to smoking history

Although smoking-related diseases, such as chronic obstructive pulmonary disease (COPD), are often accompanied by increased peripheral muscle fatigability, the extent to which this is a feature of the disease or a direct effect of smoking per se is not known. Skeletal muscle function was investigated in terms of maximal voluntary isometric torque, activation, contractile properties and fatigability, using electrically evoked contractions of the quadriceps muscle of 40 smokers [19 men and 21 women; mean (SD) cigarette pack years: 9.9 (10.7)] and age- and physical activity level matched non-smokers (22 men and 23 women). Maximal strength and isometric contractile speed did not differ significantly between smokers and non-smokers. Muscle fatigue (measured as torque decline during a series of repetitive contractions) was greater in smokers (P = 0.014), but did not correlate with cigarette pack years (r = 0.094, P = 0.615), cigarettes smoked per day (r = 10.092, P = 0.628), respiratory function (%FEV1pred) (r = −0.187, P = 0.416), or physical activity level (r = −0.029, P = 0.877). While muscle mass and contractile properties are similar in smokers and non-smokers, smokers do suffer from greater peripheral muscle fatigue. The observation that the cigarette smoking history did not correlate with fatigability suggests that the effect is either acute and/or reaches a ceiling, rather than being cumulative. An acute and reversible effect of smoking could be caused by carbon monoxide and/or other substances in smoke hampering oxygen delivery and mitochondrial function

Aftershock seismicity of the 2010 Maule Mw=8.8, Chile, earthquake: Correlation between co-seismic slip models and aftershock distribution?

International audienceThe 27 February 2010 Maule, Chile (Mw=8.8) earthquake is one of the best instrumentally observed subduction zone megathrust events. Here we present locations, magnitudes and cumulative equivalent moment of the first $2 months of aftershocks, recorded on a temporary network deployed within 2 weeks of the occurrence of the main-shock. Using automatically-determined onset times and a back projection approach for event association, we are able to detect over 30,000 events in the time period analyzed. To further increase the location accuracy, we systematically searched for potential S-wave arrivals and events were located in a regional 2D velocity model. Additionally, we calculated regional moment tensors to gain insight into the deformation history of the aftershock sequence. We find that the aftershock seismicity is concentrated between 40 and 140 km distance from the trench over a depth range of 10 to 35 km. Focal mechanisms indicate a predominance of thrust faulting, with occasional normal faulting events. Increased activity is seen in the outer-rise region of the Nazca plate, predominantly in the northern part of the rupture area. Further down-dip, a second band of clustered seismicity, showing mainly thrust motion, is located at depths of 40-45 km. By comparing recent published mainshock source inversions with our aftershock distribution, we discriminate slip models based on the assumption that aftershocks occur in areas of rapid transition between high and low slip, surrounding high-slip regions of the mainshock