No abstract

The large low shear velocity provinces observed in the lowermost mantle are widely accepted as chemically distinct thermochemical “piles,” but their origin and long‐term evolution remain poorly understood. The survival time and shape of the large low shear velocity provinces are thought to be mainly controlled by their compositional density, while their viscosity has been considered less important. Based on recent constraints on chemical reactions between mantle and core, a more complex viscosity structure of the lowermost mantle, possibly including high viscosity thermochemical pile material, seems reasonable. In this study, we use numerical models to identify a trade‐off between compositional viscosity and density contrasts required for long‐term stability of thermochemical piles, which permits lower‐density and higher‐viscosity piles. Moreover, our results indicate more restrictive stability conditions during periods of strong deformation‐induced entrainment, for example, during initial pile formation, which suggests long‐term pile survival.

We investigate the upper mantle seismic discontinuities at 410 and 660 km depth beneath the Indian Ocean Geoid Low (IOGL). To map the discontinuities' topography, we use differential travel times of PP and SS waves and their precursors. Our final data set consists of 37 events with Mw ≥ 5.8, which densely cover our investigation area, also with crossing ray paths. We use array methods to detect the low‐amplitude precursor signals. The best quality data show a deepened 410 km discontinuity in the center of the IOGL as well as a mostly elevated 660 km discontinuity beneath the northern Indian Ocean, which we interpret as a hot anomaly currently residing in the mantle transition zone. We conclude that the largest negative geoid anomaly might be caused by a combined effect of hot material in the midmantle below the innermost IOGL and cold material below 660 km farther south.

While hotspot tracks beneath thin oceanic lithosphere are visible as volcanic island chains, the plume‐lithosphere interaction for thick continental or cratonic lithosphere often remains hidden due to the lack of volcanism. To identify plume tracks with missing volcanism, we characterize the amplitude and timing of surface heat flux anomalies following a plume‐lithosphere interaction using mantle convection models. Our numerical results confirm an analytical relationship in which surface heat flux increases with the extent of lithosphere thinning, which is primarily controlled by the viscosity structure of the lower lithosphere and the asthenosphere. We find that lithosphere thinning is greatest when the plate is above the plume conduit, while the maximum heat flux anomaly occurs about 40–140 Myr later. Therefore, younger continental and cratonic plume tracks can be identified by observed lithosphere thinning, and older tracks by an increased surface heat flux, even if they lack extrusive magmatism.