We used Gaussian Process Empirical Dynamic Models (GP-EDM) to forecast sockeye run sizes. We considered a variety of predictors, including total returns, returns of each age class, and several environmental variables (NPGO, PDO, PNI, and CUTI). All models were fit in R using the GPEDM package (https://tanyalrogers.github.io/GPEDM). Data from prior to 2016 were used to train each model, and then one-year-ahead predictions (with sequential model updating) were made on the 2016-2025 data. Of the models that achieved similar median APEs on the ten year test set, a single model was selected for the Stellako, Raft, and Late Stuart populations. For the Chilko and Quesnel populations, different model classes performed similarly but produced notably different forecasts, so model averaging was performed.