A Connection between Geometrical Spreading and the Adjoint Field in Travel Time Tomography
- 1 Lamont-Doherty Earth Observatory of Columbia University, Palisades, New York, USA
Abstract
The goal of tomography is to reconstruct a spatially-varying image function s ( x , m ), where x is position and m is a finite-length vector of parameters. Many reconstruction methods minimize the total L 2 error E ≡ e T e , where individual errors e i quantify misfit between predictions and observations, to quantify goodness of fit. So-called adjoint state methods allow the gradient <i>∂</i> E / <i>∂</i> m i to be computed extremely efficiently from an adjoint field, facilitating image reconstruction by gradient-descent methods. We examine the structure of the differential equation for the adjoint field under the ray approximation and find that it has the same form as the transport equation, whose solution involves the well-known geometrical spreading function R Consequently, as R is routinely tabulated as part of a ray calculation, no extra work is needed to compute the adjoint field, permitting a rapid calculation of the gradient <i>∂</i> E / <i>∂</i> m i .
- Aki, K., Christoffersson, A. and Husebye, E.S. (1976) Three-Dimensional Seismic Structure of the Lithosphere under Montana Lasa. Bulletin of the Seismological Society of America, 66, 501-524.
- Menke, W. (1977) Lateral in Homogeneities in P Velocity under the Tarbella Array of the Lesser Himalayas of Pakistan. Bulletin of the Seismological Society of America, 67, 725-734.
- Munk, W., Worcester, P. and Wunsch, C. (1995) Ocean Acoustic Tomography. Cambridge University Press, Cambridge, 433 p.
- Justice, J.H., Vassiliou, A.A., Singh, S., Logel, J.D., Hansen, P.A., Hall, B.R., Hurt, P.R. and Solanki, J.J. (1989) Acoustic Tomography for Monitoring Enhanced Oil Recovery. Leading Edge, 8, 12-19. https://doi.org/10.1190/1.1439605
- Santamarina, J.C. (1994) An Introduction to Geotomography. In: Woods, R.D., Ed., Geophysical Characterization of Sites, Oxford & IBH Publishing Co., New Delhi, 35-44.
- Müller, T. (2014) GeoViS—Relativistic Ray Tracing in Four-Dimensional Spacetimes. Computer Physics Communications, 185, 2301-2308. https://doi.org/10.1016/j.cpc.2014.04.013
- Nolet, G. (1987) Seismic Wave Propagation and Seismic Tomography. In: Nolet, G., Ed., Seismic Tomography with Applications in Global Seismology and Exploration Geophysics, Springer, New York, 1-23.
- Dahlen, F., Hung, S.-H. and Nolet, G. (2002) Fréchet Kernels for Finite-Frequency Travel Times—I. Theory. Geophysical Journal International, 141, 157-174. https://doi.org/10.1046/j.1365-246X.2000.00070.x
- Hall, M.C.G., Cacuci, D.G. and Schlesinger, M.E. (1982) Sensitivity Analysis of a Radiative-Convective Model by the Adjoint Method. Journal of the Atmospheric Sciences, 29, 2038-2050. https://doi.org/10.1175/1520-0469(1982)039 2.0.CO;2
- Bin Waheed, U., Flagg, G. and Yarman, C.E. (2016) First-Arrival Traveltime Tomography for Anisotropic Media Using the Adjoint-State Method. Geophysics, 81, R147-R155. https://doi.org/10.1190/geo2015-0463.1
- Levenberg, K. (1944) A Method for the Solution of Certain Non-Linear Problems in Least-Squares. Quarterly of Applied Mathematics, 2, 164-168. https://doi.org/10.1090/qam/10666
- Engl, H. (1993) Regularization Methods for the Stable Solution of Inverse Problems. Surveys in Mathematics for Industry, 3, 71-143.
- Menke, W. and Eilon, Z. (2015) Relationship between Data Smoothing and the Regularization of Inverse Problems. Pure and Applied Geophysics, 172, 2711-2726. https://doi.org/10.1007/s00024-015-1059-0