jats:titleSUMMARY</jats:title> jats:pWe apply the adjoint method to efficiently calculate the linearized sensitivity of body tide observations to perturbations in density, elastic/anelastic moduli and boundary topography. This theory is implemented practically within the context of normal mode coupling calculations, with an advantage of this approach being that much of the necessary technical machinery is present in existing coupling codes. A range of example sensitivity kernels are calculated relative to both spherically symmetric and laterally heterogeneous background models. These results reaffirm the conclusions of earlier studies that the M2 body tide is strongly sensitive to spherical harmonic degree-2 density variations at the base of the mantle. Moreover, it is found that the sensitivity kernels are only weakly dependent on the background model, and hence linearized methods are likely to be effective within inversions of body tide observations.</jats:p>