Search

When interpreting marine very high-resolution (VHR) single-channel seismic reflection data, the signal in the water column is generally considered as noise and is often eliminated by a water-mute application to focus on geological information under the seafloor. Alternatively, the signal in the water column can be used to study ocean currents or gas/fluid emissions. To provide images of the sedimentary formations and tectonic structures beneath the seafloor in shallow water regions, such as continental shelves and lakes, marine seismic reflection profiles are often acquired using a single-channel streamer and sparker-type source, providing VHR data, with limited penetration depth. To exploit the full potential of these single-channel data, we propose a simple algorithm, called REWARE (Recovery of Water-column Acoustic Reflectors). This algorithm allows to extract further geological information from the water-column data using open-source codes (Seismic Un*x), adding the coherent signal from the previous shots, recorded in the water column, to the previous traces. The record length becomes longer while maintaining a very high trace-to-trace consistency. To demonstrate its efficiency, we present two examples of the REWARE processing in two different geological contexts: the East Sardinia shelf (Italy) and the North Evia Gulf (Greece). This method provides deeper images than with original data for seismic data acquired across steep slopes, such as canyons or continental shelf breaks. Thus, depending on the seafloor geometry and subseafloor structures, it is possible to image or map sediment layers and tectonic structures at depth, keeping a very high structural resolution. [ABSTRACT FROM AUTHOR]

We present a computational technique to model hydroacoustic waveforms from teleseismic earthquakes recorded by mid-column Mermaid floats deployed in the Pacific, taking into consideration bathymetric effects that modify seismo-acoustic conversions at the ocean bottom and acoustic wave propagation in the ocean layer, including reverberations. Our approach couples axisymmetric spectral-element simulations performed for moment-tensor earthquakes in a 1-D solid Earth to a 2-D Cartesian fluid–solid coupled spectral-element simulation that captures the conversion from displacement to acoustic pressure at an ocean-bottom interface with accurate bathymetry. We applied our workflow to 1129 seismograms for 682 earthquakes from 16 Mermaid s (short for Mobile Earthquake Recording in Marine Areas by Independent Divers) owned by Princeton University that were deployed in the Southern Pacific as part of the South Pacific Plume Imaging and Modeling (SPPIM) project. We compare the modelled synthetic waveforms to the observed records in individually selected frequency bands aimed at reducing local noise levels while maximizing earthquake-generated signal content. The modelled waveforms match the observations very well, with a median correlation coefficient of 0.72, and some as high as 0.95. We compare our correlation-based traveltime measurements to measurements made on the same data set determined by automated arrival-time picking and ray- traced traveltime predictions, with the aim of opening up the use of Mermaid records for global seismic tomography via full-waveform inversion. [ABSTRACT FROM AUTHOR]

Accurate timing corrections for seismic data recorded by ocean bottom seismometers (OBSs) are essential for a wide range of applications. The synchronization of internal OBS clocks with Global Positioning System (GPS) is typically only possible prior to deployment on the seafloor and upon retrieval. Thus, untracked, clock errors in seismic data may accumulate over the deployment period. The measurement of the clock's offset from GPS at retrieval, referred to as 'skew', can be used to correct the data solely under the assumption of a uniform rate of clock drift throughout the whole deployment. However, clock errors can be non-linear. We, therefore, develop a new workflow along with an associated open-source, interactive graphical user interface to estimate clock drift of large aperture OBS arrays. We use the workflow to estimate OBS clock drift curves for 40 OBSs of the large-scale UPFLOW seafloor array in the Madeira–Azores–Canaries region deployed for ∼14 months in 2021–2022. We use the relative shift of daily empirical Green's functions obtained from seismic ambient noise recorded by all available data channels to track clock error. We find that 95 per cent of our OBS clock drift observations have a substantial non-linear component: most maximum deviations to linearity are ∼0.75–1 s (and up to 2 s) occurring mainly halfway through the deployment. We test our drift curves by using them to correct teleseismic earthquake recordings, which enables larger numbers of high-quality P -wave traveltime measurements than when using linear drift corrections. Our drift curves have on average an uncertainty of ∼0.11 s indicating the suitability of the corrected data for future seismological studies such as for seismic tomography, seismicity analysis and moment tensor inversions. [ABSTRACT FROM AUTHOR]

Accurate seismic phase detection and onset picking are fundamental to seismological studies. Supervised deep-learning phase pickers have shown promise with excellent performance on land seismic data. Although it may be acceptable to apply them to Ocean Bottom Seismometer (OBS) data that are indispensable for studying ocean regions, they suffer from a significant performance drop. In this study, we develop a generalized transfer-learned OBS phase picker—OBSTransformer, based on automated labelling and transfer learning. First, we compile a comprehensive data set of catalogued earthquakes recorded by 423 OBSs from 11 temporary deployments worldwide. Through automated processes, we label the P and S phases of these earthquakes by analysing the consistency of at least three arrivals from four widely used machine learning pickers (EQTransformer, PhaseNet, Generalized Phase Detection and PickNet), as well as the Akaike Information Criterion (AIC) picker. This results in an inclusive OBS data set containing ∼36 000 earthquake samples. Subsequently, we use this data set for transfer learning and utilize a well-trained land machine learning model—EQTransformer as our base model. Moreover, we extract 25 000 OBS noise samples from the same OBS networks using the Kurtosis method, which are then used for model training alongside the labelled earthquake samples. Using three groups of test data sets at subglobal, regional and local scales, we demonstrate that OBSTransformer outperforms EQTransformer. Particularly, the P and S recall rates at large distances (>200 km) are increased by 68 and 76 per cent, respectively. Our extensive tests and comparisons demonstrate that OBSTransformer is less dependent on the detection/picking thresholds and is more robust to noise levels. [ABSTRACT FROM AUTHOR]

Distributed acoustic sensing (DAS) is becoming increasingly popular in microseismic monitoring operations. This data acquisition technology converts fibre-optic cables into dense arrays of seismic sensors that can sample the seismic wavefield produced by active or passive sources with a high spatial density, over distances ranging from a few hundred metres to tens of kilometres. However, standard microseismic data analysis procedures have several limitations when dealing with the high spatial (intersensor spacing up to submetre scale) sampling rates of DAS systems. Here, we propose a semblance-based seismic event detection method that fully exploits the high spatial sampling of the DAS data. The detector identifies seismic events by computing waveform coherence of the seismic wavefield along geometrical hyperbolic trajectories for different curvatures and positions of the vertex, which are completely independent from external information (i.e. velocity models). The method detects a seismic event when the coherence values overcome a given threshold and satisfies our clustering criteria. We first validate our method on synthetic data and then apply it to real data from the FORGE geothermal experiment in Utah, USA. Our method detects about two times the number of events obtained with a standard method when applied to 24 hr of data. [ABSTRACT FROM AUTHOR]

Over the past several years, the use of optical fibre cables as ground motion sensors has become a central topic for seismologists, with successful applications of distributed acoustic sensing (DAS) in various key fields such as seismic monitoring, structural imaging and source characterization. DAS response is a combination of both instrument response and cable-ground coupling, with the latter having a strong, spatially variable, but yet largely unquantified effect. This limits the application of a large number of staple seismological techniques (e.g. earthquake magnitude estimation, waveform tomography) that can require accurate knowledge of a signal's amplitude and frequency content. Here we present a method for accurately simulating a DAS cable and its ground coupling. The scheme is based on molecular dynamic-like particle-based numerical modelling, allowing the investigation of the effect of varying DAS-ground coupling scenarios. We start by computing the full strain field directly, for each pair of neighbouring particles in the model. We then define a virtual DAS cable, embedded within the model and formed by a single string of interconnected particles. This allows us to control all aspects of the cable-ground coupling and their properties at an effective granular level through changing the bond stiffness and bond types (e.g. non-linearity) for both the cable and the surrounding medium. Arbitrary cable geometries and heterogeneous materials can be accommodated at the desired scale of investigation. We observe that at the meter scale, the cable-ground coupling and local site effects can substantially alter the recorded signal. We find that the stiffness of the thin layer of material to which the cable is coupled has the strongest effects, selectively amplifying portions of the wave train and contributing to substantial phase delays. These differences show that cable coupling and local site effects should be considered both when designing a DAS deployment and analysing its data when either true or along-cable relative amplitudes and/or frequencies are considered. The codes developed herein for calculating full waveform DAS responses and coupling are made publicly available. [ABSTRACT FROM AUTHOR]

Over a period of less than a decade, distributed acoustic sensing (DAS) has become a well-established technology in seismology. For historical and practical reasons, DAS manufacturers usually provide instruments that natively record strain (rate) as the principal measurement. While at first glance strain recordings seem related to ground motion waveforms (displacement, velocity and acceleration), not all the seismological tools developed over the past century (e.g. magnitude estimation, seismic beamforming, etc.) can be readily applied to strain data. Notably, the directional sensitivity of DAS is more limited than conventional particle motion sensors, and DAS experiences an increased sensitivity to slow waves, often highly scattered by the subsurface structure and challenging to analyse. To address these issues, several strategies have been already proposed to convert strain rate measurements to particle motion. In this study, we focus on strategies based on a quantity we refer to as 'deformation'. Deformation is defined as the change in length of the cable and is closely related to displacement, yet both quantities differ from one another: deformation is a relative displacement measurement along a curvilinear path. We show that if the geometry of the DAS deployment is made of sufficiently long rectilinear sections, deformation can be used to recover the displacement without the need of additional instruments. We validate this theoretical result using full-waveform simulations and by comparing, on a real data set, the seismic velocity recovered from DAS with that recorded by collocated seismometers. The limitations of this approach are discussed, and two applications are shown: enhancing direct P -wave arrivals and simplifying the magnitude estimation of seismic events. Converted displacement provides better sensitivity to high velocity phases, improves broadside response and permits the direct application of conventional seismological tools that are less effective when applied to strain (rate) data. [ABSTRACT FROM AUTHOR]

Extreme weather events threaten life and property in populated areas. Timely and precise weather event monitoring and risk assessment are critical, but can be hampered by limited meteorological stations in cities. Recent studies have shown that seismic stations are sensitive to storm-induced noise. This study aims to investigate the sensitivity of distributed acoustic sensing (DAS), a technology capable of turning existing optical fibres into dense seismic sensors, for recording storm-induced seismic noise. We analyse 4-month continuous DAS recordings (June–September 2021) from a 4.2-km-long underground fibre-optic array in State College, PA. We calibrate the DAS data by comparing it to various meteorological data (rainfall and wind speed) from nearby weather stations. We first show DAS-recorded low-frequency wind-induced noise (0.5–8 Hz) probably caused by light poles swaying in the wind, as observed resonant frequencies agree with theoretical natural frequencies of nearby light poles. We find a strong linear correlation between DAS energy and wind speed. We further characterize rain-induced noise. Detailed observations from two rain events: a moderate rain and a heavy rain from Hurricane Ida, suggest that rain-induced noise is not generated by direct raindrops hitting the ground. Instead, the low-frequency noise (2–8 Hz) is attributed to the acoustic noise generated by water flow in stormwater drainage systems. Strong high-frequency noise up to 125 Hz is likely related to the rapid rainwater filling from the surface to the drainage system during heavy rain. We show linear relations between rain-induced DAS energy and rainfall rate, where the slopes of relations are related to the volume of rainwater, suggesting the influence of surface water and rainwater flow in the drainage system on DAS signals. Our results show the possibility of using DAS-equipped pre-existing telecom fibre-optic cables for sensing windstorms and rainstorms in urban areas and their interactions with urban infrastructures. [ABSTRACT FROM AUTHOR]

Ocean-bottom seismometers (OBSs) are equipped with seismic sensors that record acoustic and seismic events at the seafloor, which makes them suitable for investigating tectonic structures capable of generating earthquakes offshore. One critical parameter to obtain accurate earthquake locations is the absolute time of the incoming seismic signals recorded by the OBSs. It is, however, not possible to synchronize the internal clocks of the OBSs with a known reference time, given that GNSS signals are unable to reach the instrument at the sea bottom. To address this issue, here we introduce a new method to synchronize the clocks of large-scale OBS deployments. Our approach relies on the theoretical time-symmetry of time-lapse (averaged) crosscorrelations of ambient seismic noise. Deviations from symmetry are attributed to clock errors. This implies that the recovered clock errors will be obscured by lapse crosscorrelations' deviations from symmetry that are not due to clock errors. Non-uniform surface wave illumination patterns are arguably the most notable source which breaks the time symmetry. Using field data, we demonstrate that the adverse effects of non-uniform illumination patterns on the recovered clock errors can be mitigated by means of a weighted least-squares inversion that is based on station-station distances. In addition, our methodology permits the recovery of timing errors at the time of deployment of the OBSs. This error can be attributed to either: i) a wrong initial time synchronization of the OBS or ii) a timing error induced by changing temperature and pressure conditions while the OBS is sunk to the ocean floor. The methodology is implemented in an open-source Python package named OCloC, and we applied it to the OBS recordings acquired in the context of the IMAGE project in and around Reykjanes, Iceland. As expected, most OBSs suffered from clock drift. Surprisingly, we found incurred timing errors at the time of deployment for most of the OBSs.

Instrumental earthquake monitoring in Greece started in 1899–1906 with the first five seismic stations of the National Observatory of Athens, Institute of Geodynamics (NOA). Subsequent upgrades and expansions led to today's Unified National Seismic Network, which includes almost all permanent seismic stations in Greece and provides waveform and parametric data. We examine the detection capabilities of the Greek seismic networks using phase, location and magnitude data from the International Seismological Centre (ISC) and the various National Institutes. We apply two methods to measure the network performance. In one, we form a grid, and find the 50 per cent probability detection threshold for each station-grid cell pair for different times. In the other, we find the probability-based magnitude of completeness grid for every three months from 1970 to 2014. Both methods show that in 1990 the detection threshold improved significantly in the north part of Greece. A much greater improvement took place in 2010 throughout Greece, yielding a completeness magnitude of |$M_p\, \sim$| 1. [ABSTRACT FROM AUTHOR]

Einstein Telescope (ET) is a proposed underground infrastructure in Europe to host future generations of gravitational-wave (GW) detectors. One of its design goals is to extend the observation band of terrestrial GW detectors from currently about 20 Hz down to 3 Hz. The coupling of a detector to its environment becomes stronger at lower frequencies, which makes it important to carefully analyse environmental disturbances at ET candidate sites. Seismic disturbances pose the greatest challenge since there are several important mechanisms for seismic vibrations to produce noise in ET, for example, through gravitational coupling, stray light, or through harmful constraints on the design of ET's control system. In this paper, we present an analysis of the time-variant properties of the seismic field at the Sardinia candidate site of ET connected to anthropogenic as well as natural phenomena. We find that temporal variations of source distributions and of the noise spectra generally follow predictable trends in the form of diurnal, weekly, or seasonal cycles. Specific seismic sources were identified such as road bridges, which produce observable disturbances underground. This information can be used to adapt a detector's seismic isolation and control system. [ABSTRACT FROM AUTHOR]

We present a basic algorithm for optimal experimental design in distributed fibre-optic sensing. It is based on the fast random generation of fibre-optic cable layouts that can be tested for their cost-benefit ratio. The algorithm accounts for the maximum available cable length, lets the cable pass through pre-defined points of interest, avoids obstacles that the cable must not traverse, permits the adaptation of geometric complexity of different cable segments and allows for the incorporation of topography. Furthermore, the algorithm can be combined with arbitrary measures of the cost-benefit ratio, and its simplicity enables easy adaptations to the needs of specific applications. In addition to a description of the basic concept, we provide examples that illustrate the circumnavigation of obstacles, the steering of geometric complexity and the cable layout optimization in the presence of topographic variations. [ABSTRACT FROM AUTHOR]

We present a new compilation and analysis of broad-band ocean bottom seismometer noise properties from 15 yr of seismic deployments. We compile a comprehensive data set of representative four-component (seismometer and pressure gauge) noise spectra and cross-spectral properties (coherence, phase and admittance) for 551 unique stations spanning 18 U.S.-led experiments. This is matched with a comprehensive compilation of metadata parameters related to instrumentation and environmental properties for each station. We systematically investigate the similarity of noise spectra by grouping them according to these metadata parameters to determine which factors are the most important in determining noise characteristics. We find evidence for improvements in similarity of noise properties when grouped across parameters, with groupings by seismometer type and deployment water depth yielding the most significant and interpretable results. Instrument design, that is the entire deployed package, also plays an important role, although it strongly covaries with seismometer and water depth. We assess the presence of traditional sources of tilt, compliance, and microseismic noise to characterize their relative role across a variety of commonly used seismic frequency bands. We find that the presence of tilt noise is primarily dependent on the type of seismometer used (covariant with a particular subset of instrument design), that compliance noise follows anticipated relationships with water depth, and that shallow, oceanic shelf environments have systematically different microseism noise properties (which are, in turn, different from instruments deployed in shallow lake environments). These observations have important implications for the viability of commonly used seismic analysis techniques. Finally, we compare spectra and coherences before and after vertical channel tilt and compliance noise removal to evaluate the efficacy and limitations of these now standard processing techniques. These findings may assist in future experiment planning and instrument development, and our newly compiled noise data set serves as a building block for more targeted future investigations by the marine seismology community. [ABSTRACT FROM AUTHOR]

The eruption of the submarine Hunga Tonga-Hunga Haʻapai (Hunga Tonga) volcano on 15 January 2022, was one of the largest volcanic explosions recorded by modern geophysical instrumentation. The eruption was notable for the broad range of atmospheric wave phenomena it generated and for their unusual coupling with the oceans and solid Earth. The event was recorded worldwide across the Global Seismographic Network (GSN) by seismometers, microbarographs and infrasound sensors. The broad-band instrumentation in the GSN allows us to make high fidelity observations of spheroidal solid Earth normal modes from this event at frequencies near 3.7 and 4.4 mHz. Similar normal mode excitations were reported following the 1991 Pinatubo (Volcanic Explosivity Index of 6) eruption and were predicted, by theory, to arise from the excitation of mesosphere-scale acoustic modes of the atmosphere coupling with the solid Earth. Here, we compare observations for the Hunga Tonga and Pinatubo eruptions and find that both strongly excited the solid Earth normal mode 0S29 (3.72 mHz). However, the mean modal amplitude was roughly 11 times larger for the 2022 Hunga Tonga eruption. Estimates of attenuation (Q) for 0S29 across the GSN from temporal modal decay give Q  = 332 ± 101, which is higher than estimates of Q for this mode using earthquake data (Q  = 186.9 ± 5). Two microbarographs located at regional distances (<1000 km) to the volcano provide direct observations of the fundamental acoustic mode of the atmosphere. These pressure oscillations, first observed approximately 40 min after the onset of the eruption, are in phase with the seismic Rayleigh wave excitation and are recorded only by microbarographs in proximity (<1500 km) to the eruption. We infer that excitation of fundamental atmospheric modes occurs within a limited area close to the site of the eruption, where they excite select solid Earth fundamental spheroidal modes of similar frequencies that are globally recorded and have a higher apparent Q due to the extended duration of atmospheric oscillations. [ABSTRACT FROM AUTHOR]

SUMMARY: Fibre-optic sensing based on transmission offer an alternative to scattering-based distributed acoustic sensing (DAS). The ability to interrogate fibres that are thousands of kilometres long opens opportunities for studies of remote regions, including ocean basins. However, by averaging deformation along the fibre, transmission systems produce integrated instead of distributed measurements. They defy traditional interpretations in terms of simple seismic phases, thereby inherently requiring a full-waveform approach. For this, we develop a formalism to calculate sensitivity kernels of transmitted optical phase changes with respect to (Earth) structure using optical phase delay measurements. We demonstrate that transmission-based sensing can effectively provide distributed measurements when optical phase delays are analysed in different time windows. The extent to which a potentially useful sensitivity coverage can be achieved depends on the fibre geometry, and specifically on its local curvature. This work establishes a theoretical foundation for tomographic inversions and experimental design using transmission-based optical sensing. [ABSTRACT FROM AUTHOR]

SUMMARY: We present a theory and conceptual examples for fibre-optic deformation sensing based on phase changes of transmitted light. As a first result, we establish an exact relation between observable phase changes and the deformation tensor along the fibre. This relation is nonlinear and includes effects related to both local changes in fibre length and deformation-induced changes of the local refractive index. In cases where the norm of the deformation tensor is much smaller than 1, a useful first-order relation can be derived. It connects phase changes to an integral over in-line strain along the fibre times the local refractive index. When spatial variations of the refractive index are negligible, this permits the calculation of phase change measurements from distributed strain measurements, for instance, from distributed acoustic sensing (DAS). An alternative form of the first-order relation reveals that a directional sensitivity determines the ability of a point along the fibre to measure deformation. This directional sensitivity is proportional to fibre curvature and spatial variability of the refractive index. In a series of simple conceptual examples, we illustrate how a seismic wavefield is represented in a phase change time-series and what the role of higher-order effects may be. Specifically, we demonstrate that variable curvature along the fibre may lead to a multiplication of seismic waves, meaning that a single seismic wave appears multiple times in a recording of optical phase changes. Furthermore, we show that higher-order effects may be observable in specific scenarios, including deformation exactly perpendicular to the fibre orientation. Though higher-order effects may be realized in controlled laboratory settings, they are unlikely to occur in seismic experiments where fibre geometries are irregular and waves asymptotically propagate in all directions with all possible polarizations as a consequence of 3-D heterogeneity. Our results provide the mathematical foundation for the analysis of emerging transmission-based fibre-optic sensing data, and their later use in seismic event characterization and studies of Earth structure. [ABSTRACT FROM AUTHOR]

Debris flows are hazardous phenomena occurring at volcanoes, and monitoring them has proved as challenging as imperative in several cases. The use of seismic instruments to measure and study the physical properties of debris flows has witnessed significant progress in the last years, with the use of improved sensors, innovative methodologies and high-resolution analysis. However, the application of such studies to the practical task of providing early warnings remains limited by the significant amount of infrastructural and technological resources commonly required for their deployment. In Ecuador, debris flows at volcanoes are detected by means of seismic instruments which are usually part of broader monitoring networks, thus requiring calibration to provide quantitative information about the flows and feed early-warning systems. In the present work, a theoretical approach based on the Buckingham Π-theorem is used to determine an expression that linearly correlates the seismic signal produced by a transiting debris flow with its discharge rate, for instruments installed in different substrata and at variable distances from the drainage. The expression is experimentally tested with Acoustic Flow Monitors and Broad-band seismometers installed in the vicinity of drainages at Tungurahua and Cotopaxi volcanoes, where actual debris flows occurred in relation to eruptive activity. The experiments consist in comparing the measured peak amplitude values of the seismic signal envelopes with the estimated peak discharge rates of several events. The results confirm the validity of the theoretical expression with linear correlations observed between the seismic amplitudes and the discharge rates, thus defining calibration expressions that can be generally applied to varied environments and instruments. The seismic instruments calibrated through this methodology can provide instantaneous and reliable predictions of debris flow discharge rates within less than an order of magnitude and only requiring limited data processing and storage. Such level of prediction could help to improve early warning systems based on seismic instruments installed in locations where more developed instrumental arrays are unavailable or unpractical. [ABSTRACT FROM AUTHOR]

Tilting of the ground due to loading by the variable atmosphere is known to corrupt very long period horizontal seismic records (below 10 mHz) even at the quietest stations. At BFO (Black Forest Observatory, SW-Germany), the opportunity arose to study these disturbances on a variety of simultaneously operated state-of-the-art broad-band sensors. A series of time windows with clear atmospherically caused effects was selected and attempts were made to model these 'signals' in a deterministic way. This was done by simultaneously least-squares fitting the locally recorded barometric pressure and its Hilbert transform to the ground accelerations in a bandpass between 100 and 3600 s periods. Variance reductions of up to 97 per cent were obtained. We show our results by combining the 'specific pressure induced accelerations' for the two horizontal components of the same sensor as vectors on a horizontal plane, one for direct pressure and one for its Hilbert transform. It turned out that at BFO the direct pressure effects are large, strongly position dependent and largely independent of atmospheric events for instruments installed on piers, while three post-hole sensors are only slightly affected. The infamous 'cavity effects' are invoked to be responsible for these large effects on the pier sensors. On the other hand, in the majority of cases all sensors showed very similar magnitudes and directions for the vectors obtained for the regression with the Hilbert transform, but highly variable from event to event especially in direction. Therefore, this direction most certainly has to do with the gradient of the pressure field moving over the station which causes a larger scale deformation of the crust. The observations are very consistent with these two fundamental mechanisms of how fluctuations of atmospheric surface pressure causes tilt noise. The results provide a sound basis for further improvements of the models for these mechanisms. The methods used here can already help to reduce atmospherically induced noise in long-period horizontal seismic records. [ABSTRACT FROM AUTHOR]

We study spatial and temporal characteristics of the microseismic noise field across Europe. Rather than focusing on the areas of noise generation, the scope of this work is to characterize, at the scale of Europe, the spatio-temporal evolution of the noise wavefield that results from the interplay of the seismic noise sources and the propagation effect. To that end, we perform single station analysis in three period bands (PB1: 2.5−5 s; PB2: 5−10 s and PB3: 10−20 s) using three-component seismic data recorded by ∼1000 broad-band stations in the time period 2011–2019. We calculate, for each period band, station and day, a set of parameters that are practically possible to apply to a large data set, yet yields insight into the spatio-temporal evolution of the wavefield. These parameters are: the total energy level, the dominant period of the Primary and Secondary microseismic peaks, the horizontal direction with the most energy, the horizontal direction of the dominant Rayleigh waves and the square root of the energy ratio between the horizontal and vertical components. The analysis of these parameters shows that the noise field in Europe is dominated by surface waves from the North Atlantic Ocean with, in PB1 and PB2, an additional and significant contribution from the eastern part of the Mediterranean Sea. The relative contribution of these two source regions depends on the season, the influence of the eastern Mediterranean Sea being strongest in summer. The map of the peak period of the Primary and Secondary microseismic peaks indicates that the relative contribution of these two source regions is frequency dependent: the period of the Primary microseismic peak exhibits an overall increase with distance to the North Atlantic Sources, because of stronger attenuation of high-frequency wave contents. By contrast, the period of the Secondary microseismic peak is simultaneously influenced by sources in both the North Atlantic Ocean and eastern Mediterranean Sea. We show that in both microseismic peaks (PB2 and PB3), the wavefield is dominated by Love waves, as the horizontal components have the highest energy at approximately 90° angle to the direction of elliptical polarization. Moreover, our results show that lateral heterogeneities in the crust have a major influence on the noise field. In particular, the propagation directions of Love and Rayleigh waves show strong dependency on location (but not on time of year), with very sharp boundaries for example at the edge of the Alps. Thus, the scattering that takes place in the heterogeneous Alpine crust partly randomizes the directions of the microseismic wavefield in particular in PB1 and PB2. Finally, we show that the temporal evolution of the amplitude ratio between the horizontal and vertical components reflects the relative amounts of surface waves from the North Atlantic Ocean with respect to body waves from sources in the Southern Hemisphere. Thus, this ratio can be used as a proxy to identify time periods where body waves are significant in the noise wavefield. [ABSTRACT FROM AUTHOR]

We present the first 16 months of data returned from a mobile array of 16 freely floating diving instruments, named mermaid for Mobile Earthquake Recording in Marine Areas by Independent Divers, launched in French Polynesia in late 2018. Our 16 are a subset of the 50  mermaid deployed over a number of cruises in this vast and understudied oceanic province as part of the collaborative South Pacific Plume Imaging and Modeling (SPPIM) project, under the aegis of the international EarthScope-Oceans consortium. Our objective is the hydroacoustic recording, from within the oceanic water column, of the seismic wavefield generated by earthquakes worldwide, and the nearly real-time transmission by satellite of these data, collected above and in the periphery of the South Pacific Superswell. This region, characterized by anomalously elevated oceanic crust and myriad seamounts, is believed to be the surface expression of deeply rooted mantle upwellings. Tomographically imaging Earth's mantle under the South Pacific with data from these novel instruments requires a careful examination of the earthquake-to- mermaid traveltimes of the high-frequency P -wave detections within the windows selected for reporting by the discrimination algorithms on board. We discuss a workflow suitable for a fast-growing mobile sensor database to pick the relevant arrivals, match them to known earthquakes in global earthquake catalogues, calculate their traveltime residuals with respect to global seismic reference models, characterize their quality and estimate their uncertainty. We detail seismicity rates as recorded by mermaid over 16 months, quantify the completeness of our catalogue and discuss magnitude–distance relations of detectability for our network. The projected lifespan of an individual mermaid is 5 yr, allowing us to estimate the final size of the data set that will be available for future study. To prove their utility for seismic tomography we compare mermaid data quality against 'traditional' land seismometers and their low-cost Raspberry Shake counterparts, using waveforms recovered from instrumented island stations in the geographic neighbourhood of our floats. Finally, we provide the first analyses of traveltime anomalies for the new ray paths sampling the mantle under the South Pacific. [ABSTRACT FROM AUTHOR]

Single-component ring lasers have provided high-resolution observations of Earth's rotation rate as well as local earthquake- or otherwise-induced rotational ground motions. Here, we present the design, construction and operational aspects of ROMY, a four-component, tetrahedral-shaped ring laser installed at the Geophysical Observatory Fürstenfeldbruck near Munich, Germany. Four equilateral, triangular-shaped ring lasers with 12 m side length provide rotational motions that can be combined to construct the complete vector of Earth's rotation from a point measurement with very high resolution. Combined with a classic broad-band seismometer, we obtain the most accurate 6 degree-of-freedom ground motion measurement system to date, enabling local and teleseismic observations as well as the analysis of ocean-generated Love and Rayleigh waves. The specific design and construction details are discussed as are the resulting consequences for permanent observations. We present seismic observations of local, regional and global earthquakes as well as seasonal variations of ocean-generated rotation noise. The current resolution of polar motion is discussed and strategies how to further improve long-term stability of the multicomponent ring laser system are presented. [ABSTRACT FROM AUTHOR]

The measurements of acoustic waves created by a quake are of great interest for planets with hot and dense atmospheres, like Venus, because surface deployments of seismometers will last only a few hours, whereas freeflying balloons could fly many days. Infrasound sensors can also be used to constrain subsurface properties during active seismic experiments. This study presents a controlled source seismo-acoustic experiment using infrasonic sensors and accelerometers mounted on a tethered helium balloon. Both the acoustic waves generated below the balloon by seismic surface waves, and the ones generated by strong ground motions above the seismic source are clearly observed and separated on the records of the various instruments. This data set allows various validations and investigations. First, it validates the ground to air coupling theory and our numerical modelling tools. Then, it allows us to demonstrate that antenna processing of infrasound sensors deployed below the balloon can estimate the arrival incidence angle of the acoustic waves within 10°. Finally, a polarization analysis of the accelerometers taped on the balloon envelope is presented. It demonstrates that accelerometer records are strongly dependent on their location on the balloon due to its deformations and rotations. However, the different acoustic signals can be distinguished through their polarization, and a best sensor location is estimated at the bottom of the balloon envelope. These results are a first step towards detecting and locating seismic activity using airborne acoustic sensors on Venus and elsewhere. However, some observations of earthquake signals in a more realistic geometry are still missing. [ABSTRACT FROM AUTHOR]

Instrumental timing and phase errors are a notorious problem in seismic data acquisition and processing. These can be frequency independent, for example due to clock drift, but may also be frequency dependent, for example due to imperfectly known instrument responses. A technique is presented that allows both types of errors to be recovered in a systematic fashion. The methodology relies on the time-symmetry usually inherent in time-averaged cross-correlations of ambient seismic noise: the difference between the arrival time of the direct surface-wave at positive time and the arrival time of the direct surface-wave at negative time is quantified. Doing this for all eligible receiver–receiver pairs of a large-N seismic array, including one or more receivers devoid of instrumental timing errors, the instrumental timing errors of all incorrectly timed receivers can be determined uniquely. Most notably, this is accomplished by means of a weighted least-squares inversion. The weights are based on the receiver–receiver distances and decrease the adverse effect of inhomogeneities in the noise illumination pattern on the recovered instrumental timing errors. Inversion results are furthermore optimized by limiting the inversion to receiver couples that (i) exceed a specific receiver–receiver distance threshold and (ii) whose time-averaged cross-correlations exceed a specific signal-to-noise ratio threshold. Potential frequency dependence of the timing errors is incorporated by means of an iterative, frequency-dependent approach. The proposed methodology is validated using synthetic recordings of ambient seismic surface-wave noise due to an arbitrary non-uniform illumination pattern. The methodology is successfully applied to time-averaged cross-correlations of field recordings of ambient seismic noise on and around the Reykjanes peninsula, SW Iceland. [ABSTRACT FROM AUTHOR]

The goodness of fit of first motion manually picked polarities from the reviewed Bulletin of the International Seismological Centre (ISC) is examined versus available earthquake mechanism solutions in order to study the behavior of seismic stations, and/or reporting agencies, over time. The probability of the reporting polarity data fitting (and not fitting) source mechanisms is quantified by calculating the probability distribution of several Bernoulli trials over a randomly perturbed set of hypocenters and velocity models for each earthquake mechanism–station polarity combination. The overall agreement of first motion polarities against source mechanisms is found to be good with a few cases of seismic stations showing indications of systematic phase reversals over certain time periods. Specifically, results were obtained for 50% of the registered stations at the ISC, and from these stations 70% show reliable operation during the operational time period under investigation, with only 3% showing the opposite, and 7% showing evidence of systematic changes in the quality of the reported first motion polarities. The rest showed great variability over short periods of time, which does not allow one to draw any conclusions. Comparing waveform data with the associated reported polarities revealed a mixture of cases of either questionable picking or true station phase reversals. The ISC is planning to make available on its website the detailed results for each station in order to make data users and seismic network operators aware of cases that need attention. [ABSTRACT FROM AUTHOR]

We tested the capability of seismic noise to monitor the stability conditions of a small rock block that we forced to fail in four following stages. Ambient vibrations were recorded with a broad-band 3C seismometer placed on top of the block and were processed to analyse their spectral and polarization characteristics with diverse algorithms. To analyse the spectral content of the records, we applied the multitaper method while seismic noise polarization features were investigated by means of the singular value decomposition of the Hermitian spectral density matrix. Numerical modelling was found to add limited value because of the uncertainty in estimating correctly spatial and mechanical features of the rock bridges between the block and the rock mass. Nevertheless, a modelling exercise we performed is in agreement with previous post-failure observations according to which unstable rocks may be coupled to the stable rock mass by rock bridges covering only a few per cent of the total surface of the fractures. Our analyses confirm that, when approaching final collapse, there is a trend of the block eigenmodes towards lower frequencies and show that polarized bands become narrower. [ABSTRACT FROM AUTHOR]

Since landing on Mars, the NASA InSight lander has witnessed eight Phobos and one Deimos transits. All transits could be observed by a drop in the solar array current and the surface temperature, but more surprisingly, for several ones, a clear signature was recorded with the seismic sensors and the magnetometer. We present a preliminary interpretation of the seismometer data as temperature‐induced local deformation of the ground, supported by terrestrial analog experiments and finite‐element modeling. The magnetic signature is most likely induced by changing currents from the solar arrays. While the observations are not fully understood yet, the recording of transit‐related phenomena with high sampling rate will allow more precise measurements of the transit times, thus providing additional constraints for the orbital parameters of Phobos. The response of the seismometer can potentially also be used to constrain the thermoelastic properties of the shallow regolith at the landing site. Plain Language Summary: The geophysical lander, InSight, has been operating on the surface of Mars since November 2018. Since then, the Martian moons Phobos and Deimos have been partially blocking the Sun, as seen from the InSight landing site, multiple times. Multiple InSight instruments have been measuring the effect of those transits; this surprisingly includes the seismometer and the magnetometer. We conclude that temperature‐induced deformation and tilt are responsible for the seismic measurements. The change observed in the magnetometer measurements are most likely the result of a drop in the solar array currents. We do not observe atmospheric modulations with InSight's weather station during the transit. These observations help constrain orbital parameters of the Martian moons, and the seismometer signal might allow investigating thermoelastic properties of the shallow Martian material. Key Points: Multiple geophysical InSight instruments observe unexpected signals during Phobos transitsLocal ground deformation due to surface temperature drops explain the tilt signals seen by the seismometerThe drop in the solar array currents results in a change of the magnetic field [ABSTRACT FROM AUTHOR]

A 2-D orthogonal distributed acoustic sensing (DAS) array designed for seismic experiments was buried horizontally beneath the Kafadar Commons Geophysical Laboratory on the Colorado School of Mines campus at Golden, Colorado. The DAS system using straight fibre-optic cables is a cost-efficient technology that enables dense seismic array deployment for long-term seismic monitoring, favouring both earthquake-based and ambient-noise-based surface wave analysis for subsurface characterization. In our study, the horizontally orthogonal DAS array records ambient noise data for a period of about two months from November 2018 to January 2019. During this time, the array also detected seismic signals from an M L3.6 earthquake at Glenwood Springs, Colorado, which exhibit opposite signal polarities in the orthogonal DAS section recordings. We derive the transformation matrix for DAS strain measurements in horizontally orthogonal cables to retrieve both Rayleigh and Love wave dispersion information from the single-component DAS signals using the 2-D multichannel analysis of surface waves method. In addition, ambient noise interferometry is applied to long-term DAS noise recordings. Our theoretical derivation demonstrates that Rayleigh and Love wave Green's functions are coupled in the noise cross-correlation functions (NCFs) of DAS receiver pairs. Stacking NCFs over the horizontally orthogonal DAS array can constructively recover the radial Rayleigh wave component but destructively suppress the Love wave component. The multimodal Monte Carlo inversion of the earthquake-based Rayleigh wave and Love wave dispersion measurements and the noise-based Rayleigh wave measurement reveals a 1-D layered structure that agrees qualitatively with geological surveys of the site. Our study demonstrates that although straight fibre-optic cables lack broadside sensitivity, using appropriate DAS array configuration and seismic array methods can extend the seismic acquisition ability of DAS and enable its application to a broad range of scenarios. [ABSTRACT FROM AUTHOR]

Breakthroughs in understanding the structure and dynamics of our planet will strongly depend upon instrumenting deep oceans. Progress has been made these last decades in ocean-bottom seismic observations, but ocean-bottom seismometer (OBS) temporary deployments are still challenging and face set-up limitations. Launched from oceanographic vessels, OBSs fall freely and may slightly drift laterally, dragged by currents. Therefore, their actual orientation and location on the landing sites are hard to assess precisely. Numerous techniques have been developed to retrieve this key information, but most of them are costly, time-consuming or inaccurate. In this work, we show how ship noise can be used as an acoustic source of opportunity to retrieve both the orientation and the location of OBSs on the ocean floor. To retrieve the OBS orientation, we developed a first method based on a combination of seismic and pressure data through the use of the acoustic intensity. This latter can be used to quantify the OBS orientation from the ship noise direction of arrival (DOA), which can then be compared with known ship trajectories obtained from the automatic identification system (AIS). To accurately relocate OBSs, we also developed a second method based on the hydrophone data which computes distances of acoustical sources by measuring time differences of arrival (TDOA) between direct and reverberated phases. The OBS location is then retrieved by fitting measured ship distances with known ship trajectories. In this study, a full network of OBSs deployed in the SW Indian Ocean was reoriented and a test station was relocated. We demonstrate that our new methods may quantify the OBS orientation with an accuracy of about one degree, and its location with an accuracy of a few tens of metres, depending on the number of ships used in the analysis. [ABSTRACT FROM AUTHOR]

Station noise levels play a fundamental limitation in our ability to detect seismic signals. These noise levels are frequency-dependent and arise from a number of physically different drivers. At periods greater than 100 s, station noise levels are often limited by the self-noise of the instrument as well as the sensitivity of the instrument to non-seismic noise sources. Recently, station operators in the Global Seismographic Network (GSN) have deployed several Streckeisen STS-6A very broad-band borehole seismometers. These sensors provide a potential replacement for the no-longer-produced Streckeisen STS-1 seismometer and the GeoTech KS-54 000 borehole seismometer. Along with showing some of the initial observational improvements from installing modern very broad-band seismometers at depth, we look at current limitations in the seismic resolution from earth tide periods 100 000 s (0.01 mHz) to Nyquist at most GSN sites (0.02 s or 50 Hz). Finally, we show the potential for improved observations of continuously excited horizontal Earth hum as well as the splitting of very long-period torsional modes. Both of these observations make use of the low horizontal noise levels which are obtained by installing very broad-band borehole seismometers at depth. [ABSTRACT FROM AUTHOR]

We use seismic noise cross-correlations to obtain a 3-D tomography model of SV -wave velocities beneath the western Indian Ocean, in the depth range of the oceanic crust and uppermost mantle. The study area covers 2000 × 2000 km2 between Madagascar and the three spreading ridges of the Indian Ocean, centred on the volcanic hotspot of La Réunion. We use seismograms from 38 ocean bottom seismometers (OBSs) deployed by the RHUM-RUM project and 10 island stations on La Réunion, Madagascar, Mauritius, Rodrigues, and Tromelin. Phase cross-correlations are calculated for 1119 OBS-to-OBS, land-to-OBS, and land-to-land station pairs, and a phase-weighted stacking algorithm yields robust group velocity measurements in the period range of 3–50 s. We demonstrate that OBS correlations across large interstation distances of >2000 km are of sufficiently high quality for large-scale tomography of ocean basins. Many OBSs yielded similarly good group velocity measurements as land stations. Besides Rayleigh waves, the noise correlations contain a low-velocity wave type propagating at 0.8–1.5 km s−1 over distances exceeding 1000 km, presumably Scholte waves travelling through seafloor sediments. The 100 highest-quality group velocity curves are selected for tomographic inversion at crustal and lithospheric depths. The inversion is executed jointly with a data set of longer-period, Rayleigh-wave phase and group velocity measurements from earthquakes, which had previously yielded a 3-D model of Indian Ocean lithosphere and asthenosphere. Robust resolution tests and plausible structural findings in the upper 30 km validate the use of noise-derived OBS correlations for adding crustal structure to earthquake-derived tomography of the oceanic mantle. Relative to crustal reference model CRUST1.0, our new shear-velocity model tends to enhance both slow and fast anomalies. It reveals slow anomalies at 20 km depth beneath La Réunion, Mauritius, Rodrigues Ridge, Madagascar Rise, and beneath the Central Indian spreading ridge. These structures can clearly be associated with increased crustal thickness and/or volcanic activity. Locally thickened crust beneath La Réunion and Mauritius is probably related to magmatic underplating by the hotspot. In addition, these islands are characterized by a thickened lithosphere that may reflect the depleted, dehydrated mantle regions from which the crustal melts where sourced. Our tomography model is available as electronic supplement. [ABSTRACT FROM AUTHOR]

In the last decade the use of passive methods has become appealing in reconstructing the properties of the propagation medium by seismic ambient noise data, without the use of localized natural or artificial sources. A temporary seismic network was installed in the urban area of Benevento (southern Italy) in order to characterize the shallow structure of the city using stable methods for the analysis of the seismic noise continuously acquired by stations. The city of Benevento is one of the italian areas with highest seismic hazard, and at present the region is affected by low energy swarms and sparse events (Ml ≤ 4.1). It has been struck by several destructive historical earthquakes, the strongest of which occurred in 1456, 1688, 1805 with associated MCS intensity up to X–XI. We used the sixteen seismic stations installed in Benevento to record ambient noise for about 1 month. The stations were equipped with different seismic instruments: (i) digitizers Quanterra Q330 connected to Le3d-5 s short-period sensors; (ii) Nanometrics Centaur digitizers coupled with Trillium Compact 120s broad-band velocimeters; (iii) one station with Episensor force balance accelerometer connected to a D6BB-DIN Staneo digitizer. Interstations Green's functions were reconstructed by the cross-correlation of continuous ambient noise data, and surface waves signals were extracted from Green's Functions (GFs) for investigating the elastic properties of the subsurface structure. In this regard, we performed the beamforming analysis to test the hypothesis of isotropy distribution of noise sources on which the cross-correlation method is based, and the particle motion analysis to confirm the presence of surface Rayleigh waves in the GFs. We analysed the temporal stability of the cross-correlated signals and the results show that 2 weeks of continuous measurements are sufficient to stabilize the surface waves signal extracted from the GFs. The phase velocity dispersion curves are computed for 115 station pairs through the use of a far-field representation of the surface-wave GFs and an image transformation technique. Our strategy based on cross-correlation analysis provides robust phase-velocity dispersion curves that vary approximately from 1.4 km s–1 at 0.7 Hz to 0.6 km s–1 at 5 Hz. Different pairs were selected for the inversion of phase-velocity dispersion curves aimed to derive 1-D shear-wave velocity (Vs) profiles (up to a maximum depth of about 500 m) representative of some areas of the city characterized by different soil deposits. [ABSTRACT FROM AUTHOR]

Long-period Rayleigh wave horizontal to vertical amplitude (H/V) ratios at a station provide information about local earth structure that is complementary to phase velocity. However, a number of studies have observed that significant scatter appears in these measurements making it difficult to use H/V ratio measurements to resolve earth structure. Some of the scatter in these measurements has been attributed to local geological structure while some has remained unaccounted for. Most Global Seismographic Network (GSN) stations contain two nearby high-quality broad-band seismometers (e.g. in the same vault, but on different piers or in different boreholes). For each broad-band sensor in the IRIS/USGS component of the GSN, we estimate H/V ratios of fundamental mode Rayleigh waves using M  > 6.5 earthquakes from 2001 to 2018 (around 19 000 measurements). We compute these ratios at a number of discrete periods (25, 50, 75, 100 and 150 s) and find that for well-isolated Rayleigh waves (windows where the correlation coefficients between radial and the phase-shifted vertical components are greater than 0.9) significant scatter in H/V ratios occurs between colocated sensors (greater than 25 per cent at 100 s period). This suggests the scatter in H/V ratio measurements can be at least partially attributed to extremely local phenomena such as sensor emplacement in the vault. We also find that H/V ratios can vary as a function of event backazimuth, indicating that care must be taken when computing average ratios for a station, as a large number of events from a given region could bias H/V ratio measurements at a station. [ABSTRACT FROM AUTHOR]

Recovery of ground displacements from strong motion accelerographs alone is still performed in an empirical way. This paper starts with a theoretical description of strong motion instruments. Then, based on observations from the non-tilt 2-D shake table tests, this paper proposes a baseline correction procedure with a cumulative correction pattern. The procedure is validated by recorded displacements of the shake tables. The paper also extends the baseline correction procedure to deal with 3-D strong motions by introducing a tilt correction pattern. The results from the proposed baseline correction procedure are compared with those from an existing empirical baseline correction procedure as well as nearby GPS measurements. The accuracy of the proposed procedure is studied considering 41 accelerographs and the corresponding GPS stations that recorded the 2011 M w 9.0 Tohoku earthquake in Japan. We find that the corrected displacements from recorded accelerations alone are not always consistent with the GPS measurements. In this specific Tohoku event of which the primary direction of motion was eastward, the ground displacements derived from strong motion stations using the proposed baseline correction procedure achieve the best agreement in the east–west direction. Meanwhile, the proposed procedure is computationally efficient. [ABSTRACT FROM AUTHOR]

We propose a theory for rotational and strain ambient noise interferometry, motivated by the recent development of rotational ground motion sensors and distributed acoustic sensing (DAS) technology. In this context, we demonstrate that displacement, strain and rotation interferograms can be generically written in the form of a representation theorem, that is, as a solution to the seismic wave equation that we refer to as the interferometric wavefield. The physical quantity (displacement, strain or rotation) determines the distributed source of the interferometric wavefield, as well as an observational operator that extracts the correct type of noise correlation function. The proposed interferometric equations are free of assumptions on the distribution of noise sources or the equipartitioning of the ambient field, typically required for Green's function retrieval. In addition to being valid for any kind of heterogeneous source and viscoelastic medium, they allow us to account for measurement details, such as the gauge length in DAS. We illustrate the practical feasibility of our approach with a series of numerical examples, based on regional-scale, spectral-element simulations of the interferometric wavefield. Specifically, we compare displacement and strain interferograms for homogeneous and heterogeneous earth models, and for homogeneous and heterogeneous noise sources. Ultimately, our developments are intended to enable adjoint-based waveform inversion with emerging measurement technologies that provide spatial gradient information in addition to conventional seismic displacement recordings. [ABSTRACT FROM AUTHOR]

Template matching has been widely applied in the detection of earthquakes and other seismic events due to its power in detecting weak signals. Recent studies using synthetics have shown that application of template matching to large-N arrays can potentially detect earthquakes substantially below the noise level. Here we apply template matching to the distributed acoustic sensing (DAS) data recorded in the Brady Hot Springs geothermal field, Nevada. Using 5 catalogued events, we detect 116 events and find 68 of them well below the noise level. We confirm 112 events are true earthquakes by examining the patterns of their sensor-to-sensor cross-correlation sections. This demonstrates that the combination of DAS and template matching has capability to detect microseismicity below the noise level, which is unusual for conventional seismic arrays and methods. With the updated catalogue, we observe a surge of earthquakes during the shutdown of a geothermal power plant nearby. In addition, the rapid increases in the downhole pressure record coincide with intense swarms of earthquakes. These observations show a strong correlation between the seismicity frequencies and the downhole pressure changes. Finally, we investigate several factors that may affect the detection performance and compare different strategies for spatial down-sampling, in order to provide helpful insights for future large-N design and data processing. [ABSTRACT FROM AUTHOR]

We present a simple and fast method to estimate preliminary earthquake location coordinates using Distributed Acoustic Sensing (DAS). Strain rate data is recorded on a 25-km long fibre-optic cable located in the northern suburbs of Athens (Greece). We apply frequency - wavenumber analysis on two segments of the cable forming an L-shaped seismic array, which offers good azimuthal sensitivity, whereas, the spatial resolution (spacing) and signal coherency control the ability to reliably resolve the apparent slowness of an incoming wavefield, and hence, to estimate the hypocentral distance to the array. We attempt to locate a local earthquake (2021/9/23, M L =3.4, Thiva) NW to the DAS array (approximately 50 km NW of Athens), and a regional earthquake (2021/10/12, M W =6.3, off the East coast of Crete) using S -wave onsets filtered within two frequency bands (0.5 Hz - 2.5 Hz and 1.0 Hz - 5.0 Hz). The obtained backazimuths agree with the observed backazimuths based on the locations reported by the Institute of Geodynamics, National Observatory of Athens (NOA) for both earthquakes, with our location for the local earthquake being roughly 10 km South of the NOA location, whereas, the location that we obtained for the regional earthquake suggested larger errors in distance, projected in our slowness estimation, possibly due to the array spatial resolution and the complex structure of the Hellenic subduction zone. [Display omitted] • An L-shaped DAS array was used in order to determine the epicentral coordinates of two earthquakes in Greece. • An array spatial extent and signal coherency trade-off controls the ability to resolve the apparent velocity of a wavefield. • The epicentre of a local earthquake was successfully calculated based on f-k analysis of DAS strain rate data. [ABSTRACT FROM AUTHOR]

To monitor the stability of a mountain slope in northern Italy, microseismic monitoring technique has been used since 2013. Locating microseismic events is a basic step of this technique. We performed a seismic tomographic survey on the mountain surface above the rock face to obtain a reliable velocity distribution in the rock mass for the localization procedure. Seismic travel-time inversion showed high heterogeneity of the rock mass with strong contrast in velocity distribution. Low velocities were found at shallow depth on the top of the rock cliff and intermediate velocities were observed in the most critical area of the rock face corresponding to a partially detached pillar. Using the 3D velocity model obtained from inversion, localization tests were performed based on the Equal Differential Time (EDT) localization method. The results showed hypocenter misfits to be around 15 m for the five geophones of the microseismic network and the error was significantly decreased compared to the results produced by a constant velocity model. Although the localization errors are relatively large, the accuracy is sufficient to distinguish microseismic events occurring in the most critical zone of the monitored rock mass from microseismic events generated far away. Thus, the 3D velocity model will be used in future studies to improve the classification of the recorded events.

Many applications in seismology rely on the accurate absolute timing of seismograms. However, both seismological land stations and ocean bottom seismometers (OBSs) can be affected by clock errors, which cause the absolute timing of seismograms to deviate from a highly accurate reference time signal, usually provided by GPS satellites. Timing problems can occur in land stations when synchronization with a GPS signal is temporarily or permanently lost. This can give rise to complicated, time-dependent clock drifts relative to GPS time, due to varying environmental conditions. Seismometers at the ocean bottom cannot receive GPS satellite signals, but operate in more stable ambient conditions than land stations. The standard protocol is to synchronize an OBS with a GPS signal immediately before deployment and after recovery. The measured timing deviation, called ‘skew’, is assumed to have accumulated linearly over the deployment interval, an assumption that is plausible but usually not verifiable. In recent years, cross-correlations of ambient microseismic noise have been put to use for correcting timing errors, but have been limited to interstation distances of at most a few tens of kilometres without reducing the temporal resolution. We apply noise cross-correlations to the evaluation of clock errors in four broad-band land stations and 53 wideband and broad-band OBSs, which were installed on and around the island of La Réunion in the western Indian Ocean during the RHUM-RUM (Réunion Hotspot and Upper Mantle-Réunions Unterer Mantel) experiment. We correlate all three seismic components, plus a hydrophone channel in OBS stations. Daily cross-correlation functions are derived for intermediate distances (∼20 km) for land-to-land station pairs; stable, 10 d stacks are obtained for very large interstation distances up to >300 km for land-to-OBS and OBS-to-OBS configurations. Averaging over multiple station pairs, and up to 16 component pairs per station, improves the accuracy of the method by a factor of four compared to the single-channel approaches of prior studies. The timing accuracy of our method is estimated to be ∼20 ms standard deviation or one sample at a sampling rate of 50 Hz. In land stations, nonlinear clock drifts and clock jumps of up to 6 min are detected and successfully corrected. For 52 out of 53 OBSs, we successfully obtain drift functions over time, which validate the common assumption of linear clock drift. Skew values that were available for 29 of these OBSs are consistent with our independent estimates within their observational error bars. For 23 OBSs that lacked skew measurements, linear OBS clock drifts range between 0.2 and 8.8 ms d−1. In addition to linear drift, three OBSs are affected by clock jumps of ∼1 s, probably indicating a missing sample problem that would otherwise have gone undetected. Thus we demonstrate the routine feasibility of high-accuracy clock corrections in land and OBSs over a wide range of interstation distances. [ABSTRACT FROM AUTHOR]

An earthquake rupture process can be kinematically described by rupture velocity, duration and spatial extent. These key kinematic source parameters provide important constraints on earthquake physics and rupture dynamics. In particular, core questions in earthquake science can be addressed once these properties of small earthquakes are well resolved. However, these parameters of small earthquakes are poorly understood, often limited by available data sets and methodologies. The Incorporated Research Institutions for Seismology Community Wavefield Experiment in Oklahoma deployed ∼350 three-component nodal stations within 40 km2 for a month, offering an unprecedented opportunity to test new methodologies for resolving small earthquake finite source properties in high resolution. In this study, we demonstrate the power of the nodal data set to resolve the variations in the seismic wavefield over the focal sphere due to the finite source attributes of an M2 earthquake within the array. The dense coverage allows us to tightly constrain rupture area using the second moment method even for such a small earthquake. The M2 earthquake was a strike-slip event and unilaterally propagated towards the surface at 90 per cent local S -wave speed (2.93 km s−1). The earthquake lasted ∼0.019 s and ruptured L c ∼70 m and W c ∼45 m. With the resolved rupture area, the stress-drop of the earthquake is estimated as 7.3 MPa for M w 2.3. We demonstrate that the maximum and minimum bounds on rupture area are within a factor of two, much lower than typical stress-drop uncertainty, despite a suboptimal station distribution. The rupture properties suggest that there is little difference between the M2 Oklahoma earthquake and typical large earthquakes. The new three-component nodal systems have great potential for improving the resolution of studies of earthquake source properties. [ABSTRACT FROM AUTHOR]

In academia and the industry, there is an increasing interest in generating and recording low seismic frequencies, which lead to better data quality, deeper signal penetration and can be important for full-waveform inversion. The common marine seismic source in acquisition is the air gun which is towed behind a vessel. The frequency content of the signal produced by the air gun mainly depends on its source depth as there are two effects that are presumed to counteract each other. First, there is the oscillating air bubble generated by the air gun that leads to more low frequencies for shallow source depths. Second, there is the interference of the downgoing wave with the first reflection from the sea surface, referred to as the ghost, which leads to more low frequencies for deeper source depths. It is still under debate whether it is beneficial to place the source shallow or deep to generate the strongest signal for frequencies below 5 Hz. Therefore, the ghost effect is studied in more detail by measuring the transmission at the water–air interface. We conduct experiments in a water tank where a small-volume seismic source is fired at different depths below the water surface to investigate how the ghost varies with frequency and depth. The signal from the seismic source is recorded with hydrophones inside water and air during the test to estimate the transmitted signal through the interface. In a second test, we perform experiments with an acoustic source located in air that is fired at different elevations above the water surface. The source in air is a starter gun and the signals are again recorded in water and air. The measured data indicate an increasing transmission of the signal through the water–air interface when the source is closer to the water surface which leads to a decreasing reflection for sources close to the surface. The measured results are compared with modelled data and the existing theory. The observed increase in transmission for shallow source depths could be explained by the theory of a spherical wave front striking the interface instead of assuming a plane wave front. The difference can be important for frequencies below 1 Hz. The results suggest that deploying a few sources very shallow during marine seismic acquisition could be beneficial for these very low frequencies. In addition, the effect of a spherical wave front might be considered for modelling far-field signatures of seismic sources for frequencies below 1 Hz. [ABSTRACT FROM AUTHOR]

In this work we present the analysis of more than 9000 signals collected from February 2013 to January 2016 by a microseismic monitoring network installed on a 300 m high limestone cliff in the Italian Prealps. The investigated area was affected by a major rockfall in 1969 and several other minor events up to nowadays. The network features five three-component geophones and a weather station and can be remotely accessed thanks to a dedicated radio link. We first manually classified all the recorded signals and found out that 95 per cent of them are impulsive broad-band disturbances, while about 2 per cent may be related to rockfalls or fracture propagation. Signal parameters in the time and frequency domains were computed during the classification procedure with the aim of developing an automatic classification routine based on linear discriminant analysis. The algorithm proved to have a hit rate higher than 95 per cent and a tolerable false alarm rate and it is now running on the field PC of the acquisition board to autonomously discard useless events. Analysis of lightning data sets provided by the Italian Lightning Detection Network revealed that the large majority of broad-band signals are caused by electromagnetic activity during thunderstorms. Crosscorrelation between microseismic signals and meteorological parameters suggests that rainfalls influence the hydrodynamic conditions of the rock mass and can trigger rockfalls and fracture propagation very quickly since the start of a rainfall event. On the other hand, temperature seems to have no influence on the stability conditions of the monitored cliff. The only sensor deployed on the rock pillar next to the 1969 rockfall scarp typically recorded events with higher amplitude as well as energy.We deem that this is due to seismic amplification phenomena and we performed ambient noise recording sessions to validate this hypothesis. Results confirm that seismic amplification occurs, although we were not able to identify any spectral peak with confidence because the sensors used are not suitable for this task. In addition, we found out that there is a preferential polarization of the wave field along the EW direction and this is in agreement with the geological analysis according to which the pillar is overhanging towards the 1969 rockfall scarp and may preferentially evolve in a wedge failure. Event location was not possible because of the lack of a velocity model of the rock mass. We tried to distinguish between near and far events by analysing the covariance matrix of the three-component recordings. Although the parameters and the outcomes of this analysis should be evaluated very carefully, it seems that about 90 per cent of the considered microseismic signals are related to the stability conditions of the monitored area. [ABSTRACT FROM AUTHOR]

Accurate determination of the seismometer orientation is a prerequisite for seismic studies including, but not limited to seismic anisotropy. While borehole seismometers on land produce seismic waveform data somewhat free of human-induced noise, they might have a drawback of an uncertain orientation. This study calculates a harmonic decomposition of teleseismic receiver functions from the P and PP phases and determines the orientation of a seismometer by minimizing a constant term in a harmonic expansion of tangential receiver functions in backazimuth near and at 0 s. This method normalizes the effect of seismic sources and determines the orientation of a seismometer without having to assume for an isotropic medium. Compared to the method of minimizing the amplitudes of a mean of the tangential receiver functions near and at 0 s, the method yields more accurate orientations in cases where the backazimuthal coverage of earthquake sources (even in the case of ocean bottom seismometers) is uneven and incomplete. We apply this method to data from the Korean seismic network (52 broad-band velocity seismometers, 30 of which are borehole sensors) to estimate the sensor orientation in the period of 2005-2016. We also track temporal changes in the sensor orientation through the change in the polarity and the amplitude of the tangential receiver function. Six borehole stations are confirmed to experience a significant orientation change (10°-180°) over the period of 10 yr. We demonstrate the usefulness of our method by estimating the orientation of ocean bottom sensors, which are known to have high noise level during the relatively short deployment period. [ABSTRACT FROM AUTHOR]

One of long-standing problems in underwater seismic studies is the inaccurate timing due to the fact the internal clock of Ocean Bottom Seismograph (OBS) is unable to synchronize with GPS. Here we present correcting large OBS clock errors and temporal drifts in a passivesource OBS array experiment in South China Sea by cross-correlating hydrophone recordings of OBS pairs. We show that, in this experiment, the noise cross-correlation function (NCCF) from hydrophone signals can retrieve higher signal-to-noise ratio (SNR) Scholte waves than the NCCF from seismometer. Because the hydrophone is positioned in the water above the seafloor, their NCCFs are thus less contaminated by more complicated solid-mode phases at the sediment-covered seafloor, leading to higher SNR for hydrophone's NCCFs. The relatively high SNR of Scholte waves enables us to use the daily NCCF (or stacked NCCFs of a few days) to constrain the temporal variations of the clock errors. A two-step approach is employed in this study to tackle large OBS clock errors: using predicted traveltimes of P phases from earthquakes to roughly correct the time, and then using the time asymmetry of two Scholte waves on NCCF to fine-tune the clock drifts. The uncertainty analysis indicates that the average error of our corrections is less than 0.2 s, suggesting the clock-corrected OBS data are valuable for seismic studies using surface waves and S waves. [ABSTRACT FROM AUTHOR]

This paper shows the results of the refined locations for underground Peaceful Nuclear Explosions (PNEs). Peaceful nuclear explosions (PNEs) were made for industrial applications in the Soviet Union. This study is based on a comparison of PNEs’ parameters. These explosions were recorded by seismographic stations in Kazakhstan from 1966 to 1988. The monitoring/verification community generally utilizes PNE locations from Sultanov et al. (1999). In reality, there are errors and some PNEs are poorly located. Our locations were determined using an integrated approach encompassing published open literature sources and archive seismogram analysis from Borovoye Geophysical Observatory. Treated PNEs seismograms have been available to researchers since 2001. They became available after the cooperation between Russian and U.S. organizations. The first one was the Institute of Geosphere Dynamics of the Russian Academy of Sciences (IDG RAS), the organization that operated the Observatory in the Soviet era. The second one was the National Nuclear Center of the Republic of Kazakhstan (NNC since 1992). The third one, from the U.S. side, was the Lamont‒Doherty Earth Observatory of Columbia University (LDEO). We present two digital seismograms of old-style seismograms from a digitized archive in ASCII format. We provide travel times for P-waves, some seismograms, and additional source parameters.

Located at the southernmost sector of the Western Branch of the East African Rift System, the Malawi Rift exemplifies an active, magma-poor, weakly extended continental rift. To investigate the controls on rifting, we image crustal and uppermost mantle structure beneath the region using ambient-noise and teleseismic Rayleigh-wave phase velocities between 9 and 100 s period. Our study includes six lake-bottom seismometers located in Lake Malawi (Nyasa), the first time seismometers have been deployed in any of the African rift lakes. Noise levels in the lake are lower than that of shallow oceanic environments and allow successful application of compliance corrections and instrument orientation determination. Resulting phase-velocity maps reveal slow velocities primarily confined to Lake Malawi at short periods (T <= 12 s), indicating thick sediments in the border-fault bounded rift basin. The slowest velocities occur within the Central Basin where Malawi Rift sedimentary strata may overlie older (Permo-Triassic) Karoo group sediments. At longer periods (T > 25 s), a prominent low-velocity anomaly exists beneath the Rungwe Volcanic Province at the northern terminus of the rift basin. Estimates of phase-velocity sensitivity indicates these low velocities occur within the lithospheric mantle and potentially uppermost asthenosphere, suggesting that mantle processes may control the association of volcanic centres and the localization of magmatism. Beneath the main portion of the Malawi Rift, a modest reduction in velocity is also observed at periods sensitive to the crust and upper mantle, but these velocities are much higher than those observed beneath Rungwe. [ABSTRACT FROM AUTHOR]

We present two independent, automated methods for estimating the absolute horizontal misorientation of seismic sensors. We apply both methods to 44 free-fall ocean-bottom seismometers (OBSs) of the RHUM-RUM experiment (http://www.rhum-rum.net/). The techniques measure the 3-D directions of particle motion of (1) P-waves and (2) Rayleigh waves of earthquake recordings. For P-waves, we used a principal component analysis to determine the directions of particle motions (polarizations) in multiple frequency passbands. We correct for polarization deviations due to seismic anisotropy and dipping discontinuities using a simple fit equation, which yields significantly more accurate OBS orientations. For Rayleigh waves, we evaluated the degree of elliptical polarization in the vertical plane in the time and frequency domain. The results obtained for the RHUM-RUM OBS stations differed, on average, by 3.1° and 3.7° between the methods, using circular mean and median statistics, which is within the methods' estimate uncertainties. Using P-waves, we obtained orientation estimates for 31 ocean-bottom seismometers with an average uncertainty (95 per cent confidence interval) of 11° per station. For 7 of these OBS, data coverage was sufficient to correct polarization measurements for underlying seismic anisotropy and dipping discontinuities, improving their average orientation uncertainty from 11° to 6° per station. Using Rayleigh waves, we obtained misorientation estimates for 40 OBS, with an average uncertainty of 16° per station. The good agreement of results obtained using the two methods indicates that they should also be useful for detecting misorientations of terrestrial seismic stations. [ABSTRACT FROM AUTHOR]

In the present paper we describe the on-land field operations integrated in the TOMO-ETNA experiment carried out in June-November 2014 at Mt. Etna volcano and surrounding areas. This terrestrial campaign consists in the deployment of 90 short-period portable three-component seismic stations, 17 Broadband seismometers and the coordination with 133 permanent seismic station belonging to Italy’s Istituto Nazionale di Geofisica e Vulcanologia (INGV). This temporary seismic network recorded active and passive seismic sources. Active seismic sources were generated by an array of air-guns mounted in the Spanish oceanographic vessel “Sarmiento de Gamboa” with a power capacity of up to 5200 cubic inches. In total more than 26,000 shots were fired and more than 450 local and regional earthquakes were recorded. We describe the whole technical procedure followed to guarantee the success of this complex seismic experiment. We started with the description of the location of the potential safety places to deploy the portable network and the products derived from this search (a large document including full characterization of the sites, owners and indication of how to arrive to them). A full technical description of the seismometers and seismic sources is presented. We show how the portable seismic network was deployed, maintained and recovered in different stages. The large international collaboration of this experiment is reflected in the participation of more than 75 researchers, technicians and students from different institutions and countries in the on-land activities. The main objectives of the experiment were achieved with great success.