The impact of wave-mean flow interaction on the Northern Hemisphere polar vortex after tropical volcanic eruptions

The current generation of Earth system models that participate in the Coupled Model Intercomparison Project phase 5 (CMIP5) does not, on average, produce a strengthened Northern Hemisphere (NH) polar vortex after large tropical volcanic eruptions as suggested by observational records. Here we investigate the impact of volcanic eruptions on the NH winter stratosphere with an ensemble of 20 model simulations of the Max Planck Institute Earth system model. We compare the dynamical impact in simulations of the very large 1815 Tambora eruption with the averaged dynamical response to the two largest eruptions of the CMIP5 historical simulations (the 1883 Krakatau and the 1991 Pinatubo eruptions). We find that for both the Tambora and the averaged Krakatau-Pinatubo eruptions the radiative perturbation only weakly affects the polar vortex directly. The position of the maximum temperature anomaly gradient is located at approximately 30 degrees N, where we obtain significant westerly zonal wind anomalies between 10 hPa and 30 hPa. Under the very strong forcing of the Tambora eruption, the NH polar vortex is significantly strengthened because the subtropical westerly wind anomalies are sufficiently strong to robustly alter the propagation of planetary waves. The average response to the eruptions of Krakatau and Pinatubo reveals a slight strengthening of the polar vortex, but individual ensemble members differ substantially, indicating that internal variability plays a dominant role. For the Tambora eruption the ensemble variability of the zonal mean temperature and zonal wind anomalies during midwinter and late winter is significantly reduced compared to the volcanically unperturbed period.