We produced carbon and nitrogen isoscapes across the entire Southern Ocean (>40°S) using surface particulate organic matter isotope data, collected over the past 50 years. We used Integrated Nested Laplace Approximation -based approaches to predict mean annual isoscapes and four seasonal isoscapes using a suite of environmental data as predictor variables.