A statistical learning approach to produce seasonal temperature forecasts in western Europe and Scandinavia was implemented and tested. The leading principal components (PCs) of sea surface temperature (SST) and the geopotential at the 150-hPa level (GPT) were derived from reanalysis datasets and used at different lags (from one to five seasons) as predictors. Random sampling of both the fitting years and the potential predictors together with the Least Absolute Shrinkage and Selection Operator regression (LASSO) was used to create a large ensemble of statistical models. Applying the models to independent test years shows that the ensemble performs well over the target areas and that the ensemble mean is more accurate than the best individu...