Locating the Nordstream explosions without a velocity model using polarization analysis

The seismic events that preceded the leaks in the Nordstream pipelines in the Baltic Sea have been interpreted as explosions on the seabed, most likely man-made. We use a polarization-based location method initially developed for marsquakes to locate the source region without a subsurface velocity model. We show that the 2 largest seismic events can be unambiguously attributed to the methane plumes observed on the sea surface. The two largest events can be located with this method, using 4 and 5 stations located around the source, with location uncertainties of 30km and 10x60km. We can further show that both events emitted seismic energy for at least ten minutes after the initial explosion, indicative of resonances in the water column or the depressurizing pipeline.Comment: 6 pages, 2 figures, submitted as fast report to Seismic

Expected seismicity and the seismic noise environment of Europa

Seismic data will be a vital geophysical constraint on internal structure of Europa if we land instruments on the surface. Quantifying expected seismic activity on Europa both in terms of large, recognizable signals and ambient background noise is important for understanding dynamics of the moon, as well as interpretation of potential future data. Seismic energy sources will likely include cracking in the ice shell and turbulent motion in the oceans. We define a range of models of seismic activity in Europa's ice shell by assuming each model follows a Gutenberg-Richter relationship with varying parameters. A range of cumulative seismic moment release between $10^{16}$ and $10^{18}$ Nm/yr is defined by scaling tidal dissipation energy to tectonic events on the Earth's moon. Random catalogs are generated and used to create synthetic continuous noise records through numerical wave propagation in thermodynamically self-consistent models of the interior structure of Europa. Spectral characteristics of the noise are calculated by determining probabilistic power spectral densities of the synthetic records. While the range of seismicity models predicts noise levels that vary by 80 dB, we show that most noise estimates are below the self-noise floor of high-frequency geophones, but may be recorded by more sensitive instruments. The largest expected signals exceed background noise by $\sim$50 dB. Noise records may allow for constraints on interior structure through autocorrelation. Models of seismic noise generated by pressure variations at the base of the ice shell due to turbulent motions in the subsurface ocean may also generate observable seismic noise.Comment: 24 pages, 11 figures, Added in supplementary information from revision submission, including 3 audio files with sonification of Europa noise records. To view attachments, please download and extract the gzipped tar source file listed under "Other formats

Bayesian parameter-estimation of Galactic binaries in LISA data with Gaussian Process Regression

The Laser Interferometer Space Antenna (LISA), which is currently under construction, is designed to measure gravitational wave signals in the milli-Hertz frequency band. It is expected that tens of millions of Galactic binaries will be the dominant sources of observed gravitational waves. The Galactic binaries producing signals at mHz frequency range emit quasi monochromatic gravitational waves, which will be constantly measured by LISA. To resolve as many Galactic binaries as possible is a central challenge of the upcoming LISA data set analysis. Although it is estimated that tens of thousands of these overlapping gravitational wave signals are resolvable, and the rest blurs into a galactic foreground noise; extracting tens of thousands of signals using Bayesian approaches is still computationally expensive. We developed a new end-to-end pipeline using Gaussian Process Regression to model the log-likelihood function in order to rapidly compute Bayesian posterior distributions. Using the pipeline we are able to solve the Lisa Data Challange (LDC) 1-3 consisting of noisy data as well as additional challenges with overlapping signals and particularly faint signals.Comment: 12 pages, 10 figure

Accelerating global parameter estimation of gravitational waves from Galactic binaries using a genetic algorithm and GPUs

The Laser Interferometer Space Antenna (LISA) is a planned space-based gravitational wave telescope with the goal of measuring gravitational waves in the milli-Hertz frequency band, which is dominated by millions of Galactic binaries. While some of these binaries produce signals that are loud enough to stand out and be extracted, most of them blur into a confusion foreground. Current methods for analyzing the full frequency band recorded by LISA to extract as many Galactic binaries as possible and to obtain Bayesian posterior distributions for each of the signals are computationally expensive. We introduce a new approach to accelerate the extraction of the best fitting solutions for Galactic binaries across the entire frequency band from data with multiple overlapping signals. Furthermore, we use these best fitting solutions to omit the burn-in stage of a Markov chain Monte Carlo method and to take full advantage of GPU-accelerated signal simulation, allowing us to compute posterior distributions in 2 seconds per signal on a laptop-grade GPU.Comment: 13 pages, 11 figure

Empirical H/V spectral ratios at the InSight landing site and implications for the martian subsurface structure

The horizontal-to-vertical (H/V) spectral ratio inversion is a traditional technique for deriving the local subsurface structure on Earth. We calculated the H/V from the ambient vibrations at different wind levels at the InSight landing site, on Mars, and also computed the H/V from the S-wave coda of the martian seismic events (marsquakes). Different H/V curves were obtained for different wind periods and from the marsquakes. From the ambient vibrations, the recordings during low-wind periods are close to the instrument self-noise level. During high-wind periods, the seismic recordings are highly contaminated by the interaction of the lander with the wind and the martian ground. Therefore, these recordings are less favourable for traditional H/V analysis. Instead, the recordings of the S-wave coda of marsquakes were preferred to derive the characteristic H/V curve of this site between 0.4 and 10 Hz. The final H/V curve presents a characteristic trough at 2.4 Hz and a strong peak at 8 Hz. Using a full diffuse wavefield approach as the forward computation and the Neighbourhood Algorithm as the sampling technique, we invert for the 1-D shear wave velocity structure at the InSight landing site. Based on our inversion results, we propose a strong site effect at the InSight site to be due to the presence of a shallow high-velocity layer (SHVL) over low-velocity units. The SHVL is likely placed below a layer of coarse blocky ejecta and can be associated with Early Amazonian basaltic lava flows. The units below the SHVL have lower velocities, possibly related to a Late Hesperian or Early Amazonian epoch with a different magmatic regime and/or a greater impact rate and more extensive weathering. An extremely weak buried low velocity layer (bLVL) between these lava flows explains the data around the 2.4 Hz trough, whereas a more competent bLVL would not generate this latter feature. These subsurface models are in good agreement with results from hammering experiment and compliance measurements at the InSight landing site. Finally, this site effect is revealed only by seismic events data and explains the larger horizontal than vertical ground motion recorded for certain type of marsquakes.The authors acknowledge National Aeronautics and Space Administration (NASA), Centre National D’études Spatiales (CNES), their partner agencies and institutions (United Kingdom Space Agency [UKSA], Swiss Space Office [SSO], Deutsches Zentrum für Luft-und Raumfahrt [DLR], Jet Propulsion Laboratory [JPL], Institut du Physique du Globe de Paris [IPGP]—Centre National de la Recherche Scientifique-École Normale Supérieure [CNRS], Eldgenössische Technische Hochschule Zürich [ETHZ], Imperial college [IC], Max Planck Institute for Solar System Research [MPS-MPG]), and the flight operations team at JPL, SEIS on Mars Operation Center (SISMOC), Mars SEIS Data Service (MSDS), Incorporated Research Institutions for Seismology–Data Management Center (IRIS-DMC) and Planetary Data System (PDS) for providing SEED Seismic Experiment for Interior Structure (SEIS) data. We acknowledge funding from (1) Swiss State Secretariat for Education, Research and Innovation (SEFRI project ‘Marsquake Service-Preparatory Phase’), (2) ETH Research grant ETH-0617–02, and (3) ETH + 02 19–1: Planet MARS. The research was carried out in part at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration (80NM0018D0004). This is the InSight contribution number 268

First Focal Mechanisms of Marsquakes

Since February 2019, NASA's InSight lander is recording seismic signals on the planet Mars, which, for the first time, allows to observe ongoing tectonic processes with geophysical methods. A number of Marsquakes have been located in the Cerberus Fossae graben system in Elysium Planitia and further west, in the Orcus Patera depression. We present a first study of the focal mechanisms of three well-recorded events (S0173a, S0183a, S0235b) to determine the processes dominating in the source region. We infer for all three events a predominantly extensional setting. Our method is adapted to the case of a single, multicomponent receiver and based on fitting waveforms of P and S waves against synthetic seismograms computed for the initial crustal velocity model derived by the InSight team. We explore the uncertainty due to the single-station limitation and find that even data recorded by one station constrains the mechanisms (reasonably) well. For the events in the Cerberus Fossae region (S0173a, S0235b) normal faulting with a relatively steep dipping fault plane is inferred, suggesting an extensional regime mainly oriented E-W to NE-SW. The fault regime in the Orcus Patera region is not determined uniquely because only the P wave can be used for the source inversion. However, we find that the P and weak S waves of the S0183a event show similar polarities to the event S0173, which indicates similar fault regimes

The Polarization of Ambient Noise on Mars

Seismic noise recorded at the surface of Mars has been monitored since February 2019, using the InSight seismometers. This noise can reach −200 dB. It is 500 times lower than on Earth at night and it increases of 30 dB during the day. We analyze its polarization as a function of time and frequency in the band 0.03–1 Hz. We use the degree of polarization to extract signals with stable polarization independent of their amplitude and type of polarization. We detect polarized signals at all frequencies and all times. Glitches correspond to linear polarized signals which are more abundant during the night. For signals with elliptical polarization, the ellipse is in the horizontal plane below 0.3 Hz. In the 0.3-1Hz high frequency band (HF) and except in the evening, the ellipse is in the vertical plane and the major axis is tilted. While polarization azimuths are different in the two frequency bands, they both vary as a function of local hour and season. They are also correlated with wind direction, particularly during the daytime. We investigate possible aseismic and seismic origins of the polarized signals. Lander or tether noise can be discarded. Pressure fluctuations transported by wind may explain part of the HF polarization but not the tilt of the ellipse. This tilt can be obtained if the source is an acoustic emission coming from high altitude at critical angle. Finally, in the evening when the wind is low, the measured polarized signals may correspond to the seismic wavefield of the Mars background noise