Adaptive phase space method with application to reflection traveltime tomography

In this work an adaptive strategy for the phase space method [5] for traveltime tomography is developed. The method uses those geodesics/rays that produces smaller mismatch with the measurements first and continue on in the spirit of layer stripping without the need to specifically defining the layers physically. The adaptive approach improves stability, efficiency and accuracy. We then extend our method to reflection traveltime tomography by incorporating broken geodesics/rays, for which a jump condition has to be imposed at the broken point for the geodesic flow. In particular we show that our method can distinguish non-broken and broken geodesics in the measurement and utilize accordingly in reflection traveltime tomography. We demonstrate that our method can recover the convex (with respect to the underlying metric) hull of unknown obstacles as well as the metric outside the convex hull.