Modeling of lower mantle seismic anisotropy beneath subduction zones