A document by Martijn van den Ende. Click on the document to view its contents.

The novel technique of distributed acoustic sensing (DAS) holds great potential for underwater seismology by transforming standard telecommunication cables, such as those currently traversing most of the world’s oceans, into dense arrays of seismo-acoustic sensors. To harness these measurements for seismic monitoring, the ability to record transient ground deformations using telecommunication fibers is investigated here by analyzing ambient noise, earthquake signals, and their associated phase velocities, on DAS records from three dark fibers in the Mediterranean Sea. The recording quality varies dramatically along the fibers and is strongly correlated with the bathymetry and the apparent phase velocities of the recorded waves. Apparent velocities are determined for several well-recorded earthquakes and used to convert DAS S-wave strain spectra to ground motion spectra. Excellent agreement is found between the spectra of nearby underwater and on-land seismometers and DAS converted spectra, when the latter are corrected for site effects. Apparent velocities greatly affect the ability to detect seismic deformations: for the same ground motions, slower waves induce higher strains and thus are more favorably detected than fast waves. The effect of apparent velocity on the ability to detect seismic phases, quantified by expected signal-to-noise ratios, is investigated by comparing signal amplitudes predicted by an earthquake ground motion model to recorded noise levels. DAS detection capabilities on underwater fibers are found to be similar to those of nearby broadband sensors, and superior to those of on-land fiber segments. The results demonstrate the great potential of underwater DAS for seismic monitoring and earthquake early warning.

Megathrust earthquakes are strongly influenced by the elastic properties of rocks surrounding the fault. However, these properties are often overestimated in numerical simulations, particularly in the shallow megathrust. Here we explore the influence that realistic depth-varying upper-plate elastic properties along the megathrust have on earthquake rupture dynamics and tsunamigenesis using 3D dynamic rupture and tsunami simulations. We compare results from three subduction zone scenarios with homogeneous and heterogeneous elastic media, and bimaterial fault. Elastic properties in the heterogeneous model follow a realistic depth-distribution derived from controlled-source tomography models of subduction zones. We assume the same friction properties for all scenarios. Simulations in the heterogeneous and homogeneous models show that rigidity variation of the country rock determines the depth-varying behavior of slip, slip rate, frequency content, and rupture time. Fault friction may provide additional constraints, but to a lesser extent. The depth-varying behavior of slip, frequency content, and rupture duration quantitatively agree with previous predictions based on worldwide data compilations, explaining the main depth-dependent traits of tsunami earthquakes and large shallow megathrust earthquakes. Large slip, slow rupture and slip rate amplification in bimaterial simulations are largely controlled by the elastic rock properties of the most compliant side of the fault, which in subduction zones is the upper plate. Large shallow slip and trenchward increasing upper-plate compliance of the heterogeneous model lead to the largest co-seismic seafloor deformation and tsunami amplitude. This highlights the importance of considering realistic variations in upper-plate rigidity to properly assess the tsunamigenic potential of megathrust earthquakes.

Pulse-like ruptures tend to be more sensitive to stress heterogeneity than crack-like ones. For instance, a stress-barrier can more easily stop the propagation of a pulse than that of a crack. While crack-like ruptures tend to homogenize the stress field within their rupture area, pulse-like ruptures develop heterogeneous stress fields. This feature of pulse-like ruptures can potentially lead to complex seismicity with a wide range of magnitudes akin to the Gutenberg-Richter law. Previous models required a friction law with severe velocity-weakening to develop pulses and complex seismicity. Recent dynamic rupture simulations show that the presence of a damaged zone around a fault can induce pulse-like rupture, even under a simple slip-weakening friction law, although the mechanism depends strongly on initial stress conditions. Here we aim at testing if fault zone damage is a sufficient ingredient to generate complex seismicity. In particular, we investigate the effects of damaged fault zones on the emergence and sustainability of pulse-like ruptures throughout multiple earthquake cycles, regardless of initial conditions. We consider a fault bisecting a homogeneous low-rigidity layer (the damaged zone) embedded in an intact medium. We conduct a series of earthquake cycle simulations to investigate the effects of two fault zone properties: damage level D and thickness H. The simulations are based on classical rate-and-state friction, the quasi-dynamic approximation and the software QDYN (https://github.com/ydluo/qdyn). Selected fully-dynamic simulations are also performed with a spectral element method. Our numerical results show the development of complex rupture patterns in some damaged fault configurations, including events of different sizes, as well as pulse-like, multi-pulse and hybrid pulse-crack ruptures. We further apply elasto-static theory to assess how D and H affect ruptures with constant stress drop, in particular the flatness of their slip profile, which is an indicator of pulse-like rupture. We find qualitative agreement between our theoretical and computational results regarding the range of damaged zone properties that enable pulse-like rupture.