Global mean surface air temperature (Tglobal) variability on subdecadal timescales can be of substantial magnitude relative to the long-term global warming signal, and such variability has been associated with considerable environmental and societal impacts. Therefore, probabilistic foreknowledge of short-term Tglobal evolution may be of value for anticipating and mitigating some course-resolution climate-related risks. Here we present a simple, empirically based methodology that utilizes only global spatial patterns of annual mean surface air temperature anomalies to predict subsequent annual Tglobal anomalies via partial least squares regression. The method\u27s skill is primarily achieved via information on the state of long-term global warming as well as the state and recent evolution of the El Niño–Southern Oscillation and the Interdecadal Pacific Oscillation. We test the out-of-sample skill of the methodology using cross validation and in a forecast mode where statistical predictions are made precisely as they would have been if the procedure had been operationalized starting in the year 2000. The average forecast errors for lead times of 1 to 4 years are smaller than naïve benchmarks on average, and they perform favorably relative to most dynamical Global Climate Models retrospectively initialized to the observed state of the climate system. Thus, this method can be used as a computationally efficient benchmark for dynamical model forecast systems
