First-arrival traveltime tomography (FATT) has been one of the commanding tomographic approaches for long wavelength velocity model building. Nowadays, FATT models are often used as initial models for more resolving imaging methods such as full-waveform inversion (FWI). In this context, improving the accuracy of FATT model is crucial to mitigate the nonlinearity of FWI. In spite of the widespread use of FATT at all scales (from near surface to global scale), it suffers from ill-posedness in terms of non-uniqueness of the solution due to the deficient information carried out by the sole traveltime attribute. We promote the use of the traveltime perturbation with respect to the source/receiver positions (horizontal component of the slowness vector or slope) as a supplement to the first-arrival traveltime (Tavakoli et al., 2018). For multi-component controlled-source experiments, we can generally take advantage of the fine source sampling to measure slopes by tracking local coherency over few neighboring traces in addition to traveltimes in common-receiver gathers, while slopes at the receiver position may be estimated by polarization analysis. The information of differential nature carried out by slopes yields a higher order sensitivity with respect to the subsurface parameters. We implement our first-arrival slope+traveltime (FASTT) tomography with eikonal solver and the matrix-free adjoint-state method, as previously proposed by Taillandier et al.(2009) for FATT. We illustrate the resolution power of FASTT relatively to FATT for coarse ocean bottom seismometer (OBS) acquisitions using two benchmarks based upon the EAGE/SEG Overthrust model at the exploration scale and a synthetic model of the eastern-Nankai subduction zone at the deep crustal scale (Górszczyk et al., 2018), then assess the quality of the tomographic models as initial models for FWI. We finally present an application to real OBS data collected during the SFJ experiment in the eastern Nankai Trough (Tokai area).

Dense acquisition are more and more available in exploration and earthquake seismology. Tomographic approaches can now consider not only travel times but also the wavefront itself across the seismic network (Zhang and Thurber, 2003; Yuan et al., 2016). For dense controlled-source seismic experiments, double differences of travel times between receivers in a common-shot gather (resp between sources in a common-receiver gather) are estimated, namely the horizontal component of the slowness vector at source and receiver positions designed as slopes. These slopes associated with the two-way traveltimes are interpreted as a reflection/diffraction from a small reflector segment or diffractor are used in tomographic inversion (Lambaré, 2008; Tavakoli F. et al., 2017). Picking of locally-coherent events leads to dense volumetric dataset and hence higher-resolution tomographic results (Guillaume et al., 2008). The reflection setting introduces implicitly another class of unknowns which are scatterer positions. Resulting inverse problem is awkward due to the intrinsic coupling between velocities and scatterer positions. The first choice alternates positions and wavespeeds. The second performs the joint estimation of the two parameter classes. The third one relies on the projection of the scatterer positions subspace onto the wavespeed subspace leading to a reduced-space inversion. This reduced-space formulation can be implemented in the slope tomography using adjoint-state method. Two focusing equations, which depend on two observables among the three available ones (two-way traveltime and one slope in 2D), gives exact solutions of positions which are injected as constraints in the slope tomography (Chauris et al., 2002). These constraints explicitly enforce the positions in the velocity estimation problem, which reduces now to a mono-variate inverse problem by minimization of single-slope residuals, not yet used. 2D synthetic (see figure) and real data case studies show faster convergence toward more accurate minimizer achieved by this variable projection method compared to the alternated and joint strategies. This method, which can be extended to 3D configurations, draws also interesting perspective for the joint hypocenter-velocity inversion problem in earthquake seismology.