Synthetic seismic anisotropy models within a slab impinging on the core-mantle boundary

The lowermost few hundreds of kilometres of the Earth's mantle are elastically anisotropic; seismic velocities vary with direction of propagation and polarization. Observations of strong seismic anisotropy correlate with regions where subducted slab material is expected. In this study, we evaluate the hypothesis that crystal preferred orientation (CPO) in a slab, as it impinges on the core-mantle boundary, is the cause of the observed anisotropy. Next, we determine if fast polarization directions seen by shear waves can be mapped to directions of geodynamic flow. This approach is similar to our previous study performed for a 2-D geodynamic model. In this study, we employ a 3-D geodynamic model with temperature-dependent viscosity and kinematic velocity boundary conditions defined at the surface of the Earth to create a broad downwelling slab. Tracers track the deformation that we assume to be accommodated by dislocation creep. We evaluate the models for the presence of perovskite or post-perovskite and for different main slip systems along which dislocation creep may occur in post-perovskite [(100),(010) and (001)]-resulting in four different mineralogical models of CPO. Combining the crystal pole orientations with single crystal elastic constants results in seismically distinguishable models of seismic anisotropy. The models are evaluated against published seismic observations by analysing different anisotropic components: the radial anisotropy, the splitting for (sub-)vertical phases (i.e. azimuthal anisotropy), and the splitting for subhorizontal phases. The patterns in radial anisotropy confirm our earlier results in 2-D. Observations of radial anisotropy and splitting in subhorizontal phases are mostly consistent with our models of post-perovskite with (010)-slip and (001)-slip. Our model of (001)-slip predicts stronger splitting than for (010)-slip for horizontally propagating phases in all directions. The strongest seismic anisotropy in this model occurs where the slab impinges on the core-mantle boundary. The azimuthal anisotropy pattern for (001)-slip shows fast axis directions at the edges of the slab (sub-)parallel to flow directions, suggesting horizontal flows may be mapped out in the lowermost mantle using seismic observations.