Horizontal-to-Vertical Spectral Ratio of Ambient Vibration Obtained with Hilbert-Huang Transform.

The Horizontal-to-Vertical Spectral Ratio (HVSR) of ambient vibration measurements is a common tool to explore near surface shear wave velocity (Vs) structure. HVSR is often applied for earthquake risk assessments and civil engineering projects. Ambient vibration signal originates from the combination of a multitude of natural and man-made sources. Ambient vibration sources can be any ground motion inducing phenomena, e.g., ocean waves, wind, industrial activity or road traffic, where each source does not need to be strictly stationary even during short times. Typically, the Fast Fourier Transform (FFT) is applied to obtain spectral information from the measured time series in order to estimate the HVSR, even though possible non-stationarity may bias the spectra and HVSR estimates. This problem can be alleviated by employing the Hilbert-Huang Transform (HHT) instead of FFT. Comparing 1D inversion results for FFT and HHT-based HVSR estimates from data measured at a well studied, urban, permanent station, we find that HHT-based inversion models may yield a lower data misfit χ2 by up to a factor of 25, a more appropriate Vs model according to available well-log lithology, and higher confidence in the achieved model.

On the Utility of Horizontal-to-Vertical Spectral Ratios of Ambient Noise in Joint Inversion with Rayleigh Wave Dispersion Curves for the Large-N Maupasacq Experiment.

Horizontal-to-Vertical Spectral Ratios (HVSR) and Rayleigh group velocity dispersion curves (DC) can be used to estimate the shallow S-wave velocity (VS) structure. Knowing the VS structure is important for geophysical data interpretation either in order to better constrain data inversions for P-wave velocity (VP) structures such as travel time tomography or full waveform inversions or to directly study the VS structure for geo-engineering purposes (e.g., ground motion prediction). The joint inversion of HVSR and dispersion data for 1D VS structure allows characterising the uppermost crust and near surface, where the HVSR data (0.03 to 10s) are most sensitive while the dispersion data (1 to 30s) constrain the deeper model which would, otherwise, add complexity to the HVSR data inversion and adversely affect its convergence. During a large-scale experiment, 197 three-component short-period stations, 41 broad band instruments and 190 geophones were continuously operated for 6 months (April to October 2017) covering an area of approximately 1500km2 with a site spacing of approximately 1 to 3km. Joint inversion of HVSR and DC allowed estimating VS and, to some extent density, down to depths of around 1000m. Broadband and short period instruments performed statistically better than geophone nodes due to the latter's gap in sensitivity between HVSR and DC. It may be possible to use HVSR data in a joint inversion with DC, increasing resolution for the shallower layers and/or alleviating the absence of short period DC data, which may be harder to obtain. By including HVSR to DC inversions, confidence improvements of two to three times for layers above 300m were achieved. Furthermore, HVSR/DC joint inversion may be useful to generate initial models for 3D tomographic inversions in large scale deployments. Lastly, the joint inversion of HVSR and DC data can be sensitive to density but this sensitivity is situational and depends strongly on the other inversion parameters, namely VS and VP. Density estimates from a HVSR/DC joint inversion should be treated with care, while some subsurface structures may be sensitive, others are clearly not. Inclusion of gravity inversion to HVSR/DC joint inversion may be possible and prove useful.

The Collaborative Seismic Earth Model: Generation 1.

We present a general concept for evolutionary, collaborative, multiscale inversion of geophysical data, specifically applied to the construction of a first-generation Collaborative Seismic Earth Model. This is intended to address the limited resources of individual researchers and the often limited use of previously accumulated knowledge. Model evolution rests on a Bayesian updating scheme, simplified into a deterministic method that honors today's computational restrictions. The scheme is able to harness distributed human and computing power. It furthermore handles conflicting updates, as well as variable parameterizations of different model refinements or different inversion techniques. The first-generation Collaborative Seismic Earth Model comprises 12 refinements from full seismic waveform inversion, ranging from regional crustal- to continental-scale models. A global full-waveform inversion ensures that regional refinements translate into whole-Earth structure.

Earthquake crisis unveils the growth of an incipient continental fault system.

Large continental faults extend for thousands of kilometres to form boundaries between rigid tectonic blocks. These faults are associated with prominent topographic features and can produce large earthquakes. Here we show the first evidence of a major tectonic structure in its initial-stage, the Al-Idrissi Fault System (AIFS), in the Alboran Sea. Combining bathymetric and seismic reflection data, together with seismological analyses of the 2016 Mw 6.4 earthquake offshore Morocco - the largest event ever recorded in the area - we unveil a 3D geometry for the AIFS. We report evidence of left-lateral strike-slip displacement, characterise the fault segmentation and demonstrate that AIFS is the source of the 2016 events. The occurrence of the Mw 6.4 earthquake together with historical and instrumental events supports that the AIFS is currently growing through propagation and linkage of its segments. Thus, the AIFS provides a unique model of the inception and growth of a young plate boundary fault system.

Results from 3D seismic tomography in northeast Venezuela

A micro earthquake survey in 1990 provided high resolution digital recordings of 313 local events. An array of 58 portable seismological stations was deployed for 89 days covering an area of about 26000km2 in northeastern Venezuela including the El Pilar fault zone and the broad plate boundary between Caribbean and South American tectonic plates. Each recording instrument was equipped with a 3-component short period geophone and provided a direct and analogue recording of three seismic channels and a time-code on audio cassettes at three gain levels with a total dynamic range over 90 dB. The inversion employed 1382 P- and 780 S- arrival times with timing errors less than 0.1s. Although for the S-onsets only phases on the horizontal components were considered, the obtained results did not fit in with the well constrained P-velocity. The results presented in cross-sections reveal features like thrust faults with imbricate structure dipping to the north along the El Pilar fault zone. Basal decollements and deepseated thrusting could be the explanations for back thrust, ratated or tilted faults.(AU)