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.

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.

Modeling the two- and three-dimensional displacement field in Lorca, Spain, subsidence and the global implications.

Land subsidence associated with overexploitation of aquifers is a hazard that commonly affects large areas worldwide. The Lorca area, located in southeast Spain, has undergone one of the highest subsidence rates in Europe as a direct consequence of long-term aquifer exploitation. Previous studies carried out on the region assumed that the ground deformation retrieved from satellite radar interferometry corresponds only to vertical displacement. Here we report, for the first time, the two- and three-dimensional displacement field over the study area using synthetic aperture radar (SAR) data from Sentinel-1A images and Global Navigation Satellite System (GNSS) observations. By modeling this displacement, we provide new insights on the spatial and temporal evolution of the subsidence processes and on the main governing mechanisms. Additionally, we also demonstrate the importance of knowing both the vertical and horizontal components of the displacement to properly characterize similar hazards. Based on these results, we propose some general guidelines for the sustainable management and monitoring of land subsidence related to anthropogenic activities.