Wavefield extrapolation in caustic-free normal ray coordinates

Normal ray coordinates are conventionally constructed from ray tracing, which inherently requires smooth velocity profiles. To use rays as coordinates, the velocities have to be smoothed further to avoid caustics, which is detrimental to the mapping process. Solving the eikonal equation numerically for a line source at the surface provides a platform to map normal rays in complex unsmoothed velocity models and avoid caustics. We implement reverse-time migration (RTM) and downward continuation in the new ray coordinate system, which allows us to obtain efficient images and avoid some of the dip limitations of downward continuation