Perturbation kernels for generalized seismological data functionals (GSDF)

In seismic waveform analysis and inversion, data functionals can be used to quantify the misfit between observed and model-predicted (synthetic) seismograms. The generalized seismological data functionals (GSDF) of Gee & Jordan quantify waveform differences using frequency-dependent phase-delay times and amplitude-reduction times measured on time-localized arrivals and have been successfully applied to tomographic inversions at different geographic scales as well as to inversions for earthquake source parameters. The seismogram perturbation kernel is defined as the Fréchet kernel of the data functional with respect to the seismic waveform from which the data functional is derived. The data sensitivity kernel, which is the Fréchet kernel of the data functional with respect to structural model parameters, can be obtained by composing the seismogram perturbation kernel with the Born kernel through the chain rule. In this paper, we extend GSDF analysis to broad-band waveforms by removing constraints on two control parameters defined in Gee & Jordan and derive the seismogram perturbation kernels for the modified GSDF analysis. The modifications given in this paper are consistent with the original GSDF theory in Gee & Jordan around the centre frequency and improve the stability of GSDF analysis towards the edges of the passband. We also present numerical examples of perturbation kernels for the modified GSDF analysis and their data sensitivity kernels using a homogenous half-space structure model and a complex 3-D structure model.