Forward Computation
Radial and vertical components of synthetic seismograms are calculated using the Propagator Matrix method of Thomson and Haskell. In the case of offshore setting, we adopt a boundary condition that requires zero traction on the seafloor. Also, we use equations of acoustic wave propagation rather than imitating the acoustic waves by setting Vs ~ 0 on the equations of elastic waves.
Water-level deconvolution is used.
A Gaussian low-pass filter is applied to the deconvolved seismograms:
where a is a tuning parameter.