Adjoint Method and Its Application in Mantle Convection
In this paper, the authors try to introduce the adjoint method in a systematic way, from the theoretical basis to its practical implementation in numerical models of mantle convection, with several examples to help the understanding. The adjoint operator of a temporally evolving system can be derived based on the perturbation theory, where a mismatch in the model output against observation is attributed to an error in the input, with their relation approximated as a first-order derivative (gradient) of the least-squared mismatch with respect to the input. For a nonlinear system, iterative processing is inevitable; the efficiency and convergence rate depend on the amount of a-priori information about the input (e. g. the initial condition). Seismic tomography, which describes the present-day mantle structures, has seen steady progress on both regional and global scales, and allows inverting mantle convection to the past.