This work examines elastic wave scattering around cavities embedded in a continuum with depth-dependent shear modulus and under conditions of plane strain. A restricted case of inhomogeneity is considered, where the Poisson ratio is fixed at 0.25 and where the density profile also varies, but proportionally to the shear modulus. For this specific case, the wave speeds remain macroscopically constant and it becomes possible to recover the exact Green functions by using an algebraic transformation method. These functions are subsequently used as kernels in a standard 2D boundary element formulation defined in the Laplace transform domain. The final step involves an inverse Laplace transformation, whereby the transient behavior of cavities in the aforementioned inhomogeneous continuum can be recovered. Two basic examples are solved, namely the circular cylindrical cavity under sudden internal explosion and under a pressure wave sweep. In the latter case, it is possible to investigate the effect that the angle of wave incidence has on the displacement and stress that develop along the cavity's perimeter, given the fact that the shear modulus is changing along the vertical direction. These examples serve to illustrate the present approach and to reveal some interesting differences that are observed in transient wave scattering phenomena between homogeneous and continuously inhomogeneous models, where the latter models yield a more realistic representation of geological formations.