Characterizing near-surface structures at the Ostia archaeological site based on instantaneous-phase coherency inversion

Abstract

Traditional least-squares full-waveform inversion (FWI) suffers from severe local minima problems in the case of the presence of strongly dispersive surface waves. In addition, recorded wavefields are often characterized by amplitude errors due to varying source coupling and incorrect 3D-to-2D geometric-spreading correction. Thus, least-squares FWI is considered less suitable for near-surface applications. We introduce an amplitude-unbiased coherency measure as a misfit function that can be incorporated into FWI. Such coherency has been used previously in phase-weighted stacking to enhance weak but coherent signals. The benefit of this amplitude-unbiased misfit function is that it can extract information uniformly for all seismic signals (surface waves, reflections, and scattered waves). Using the adjoint-state method, we show how to calculate the gradient of this new misfit function. We validate the robustness of the new approach using checkerboard tests and synthetic data contaminated by random noise. We then apply the new FWI approach to a field data set acquired at an archaeological site located in Ostia, Italy. The goal of this survey is to map the unexcavated archaeological remains with high resolution. We identify a known tumulus in the FWI results. The instantaneous-phase coherency FWI results also establish that the shallow subsurface under the survey lines is quite heterogeneous. The instantaneous-phase coherency FWI of near-surface data can be a promising tool to image shallow small-scale objects buried under shallow soil covers, as found at archaeological sites.