A new back-and-forth iterative method for time-reversed convection modeling: Implications for the Cenozoic evolution of 3-D structure and dynamics of the mantle

The 3-D distribution of buoyancy in the convecting mantle drives a suite of convection-related manifestations. Although seismic tomography is providing increasingly resolved images of the present-day mantle heterogeneity, the distribution of mantle density variations in the geological past is unknown, and, by implication, this is true for the convection-related observables. The one major exception is tectonic plate motions, since geologic data are available to estimate their history and they currently provide the only available constraints on the evolution of 3-D mantle buoyancy in the past. We developed a new back-and-forth iterative method for time-reversed convection modeling with a procedure for matching plate velocity data at different instants in the past. The crucial aspect of this reconstruction methodology is to ensure that at all times plates are driven by buoyancy forces in the mantle and not vice versa. Employing tomography-based retrodictions over the Cenozoic, we estimate the global amplitude of the following observables: dynamic surface topography, the core-mantle boundary ellipticity, the free-air gravity anomalies, and the global divergence rates of tectonic plates. One of the major benefits of the new data assimilation method is the stable recovery of much shorter wavelength changes in heterogeneity than was possible in our previous work. We now resolve what appears to be two-stage subduction of the Farallon plate under the western U.S. and a deeply rooted East African Plume that is active under the Ethiopian volcanic fields during the Early Eocene.