Propagation of the Hawaiian-Emperor volcano chain by Pacific plate cooling stress

The lithosphere crack model, the main alternative to the mantle plume model for age-progressive magma emplacement along the Hawaiian-Emperor volcano chain, requires the maximum horizontal tensile stress to be normal to the volcano chain. However, published stress fields calculated from Pacific lithosphere tractions and body forces (e.g., subduction pull, basal drag, lithosphere density) are not optimal for southeast propagation of a stress-free, vertical tensile crack coincident with the Hawaiian segment of the Hawaiian-Emperor chain. Here we calculate the thermoelastic stress rate for present-day cooling of the Pacific plate using a spherical shell finite element representation of the plate geometry. We use observed seafloor isochrons and a standard model for lithosphere cooling to specify the time dependence of vertical temperature profiles. The calculated stress rate multiplied by a time increment (e.g., 1 m.y.) then gives a thermoelastic stress increment for the evolving Pacific plate. Near the Hawaiian chain position, the calculated stress increment in the lower part of the shell is tensional, with maximum tension normal to the chain direction. Near the projection of the chain trend to the southeast beyond Hawaii, the stress increment is compressive. This incremental stress field has the form necessary to maintain and propagate a tensile crack or similar lithosphere flaw and is thus consistent with the crack model for the Hawaiian volcano chain.