A numerical study of caustics and polarizations of seismic waves interacting with three-dimensional geological structures