Abstract: In December 2019, the 13th revision of the International Geomagnetic Reference Field (IGRF) was released by the International Association of Geomagnetism and Aeronomy (IAGA) Division V Working Group V-MOD. This revision comprises two new spherical harmonic main field models for epochs 2015.0 (DGRF-2015) and 2020.0 (IGRF-2020) and a model of the predicted secular variation for the interval 2020.0 to 2025.0 (SV-2020-2025). The models were produced from candidates submitted by fifteen international teams. These teams were led by the British Geological Survey (UK), China Earthquake Administration (China), Universidad Complutense de Madrid (Spain), University of Colorado Boulder (USA), Technical University of Denmark (Denmark), GFZ German Research Centre for Geosciences (Germany), Institut de physique du globe de Paris (France), Institut des Sciences de la Terre (France), Pushkov Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation (Russia), Kyoto University (Japan), University of Leeds (UK), Max Planck Institute for Solar System Research (Germany), NASA Goddard Space Flight Center (USA), University of Potsdam (Germany), and Université de Strasbourg (France). The candidate models were evaluated individually and compared to all other candidates as well to the mean, median and a robust Huber-weighted model of all candidates. These analyses were used to identify, for example, the variation between the Gauss coefficients or the geographical regions where the candidate models strongly differed. The majority of candidates were sufficiently close that the differences can be explained primarily by individual modeling methodologies and data selection strategies. None of the candidates were so different as to warrant their exclusion from the final IGRF-13. The IAGA V-MOD task force thus voted for two approaches: the median of the Gauss coefficients of the candidates for the DGRF-2015 and IGRF-2020 models and the robust Huber-weighted model for the predictive SV-2020-2025. In this paper, we document the evaluation of the candidate models and provide details of the approach used to derive the final IGRF-13 products. We also perform a retrospective analysis of the IGRF-12 SV candidates over their performance period (2015–2020). Our findings suggest that forecasting secular variation can benefit from combining physics-based core modeling with satellite observations.

Abstract: There is a growing interest in studying consistency and site fidelity of individuals to assess, respectively, how individual behaviour shapes the population response to environmental changes, and to highlight the critical habitats needed by species. In Antarctica, the foraging activity of central place foragers like Adélie penguins (Pygoscelis adeliae) is constrained by the sea-ice cover during the breeding season. We estimated the population-level repeatability in foraging trip parameters and sea-ice conditions encountered by birds across successive trips over several years, and we examined their foraging site fidelity linked to sea-ice concentrations throughout the chick-rearing season. Penguins’ foraging activity was repeatable despite varying annual sea-ice conditions. Birds’ site fidelity is constrained by both sea-ice conditions around the colony that limit movements and resources availability, and also behavioural repeatability of individuals driven by phenological constraints. Adélie penguins favoured sea-ice concentrations between 20–30%, as these facilitate access to open water while opening multiple patches for exploration in restricted areas in case of prey depletion. When the sea-ice concentration became greater than 30%, foraging site fidelity decreased and showed higher variability, while it increased again after 60%. Between two trips, the foraging site fidelity remained high when sea-ice concentration changed by ± 10% but showed greater variability when sea-ice concentrations differed on a larger range. In summary, Adélie penguins specialize their foraging behaviour during chick-rearing according to sea-ice conditions to enhance their reproductive success. The balance between being consistent under favourable environmental conditions vs. being flexible under more challenging conditions may be key to improving foraging efficiency and reproductive success to face fast environmental changes.

Abstract: Sea stars (Echinodermata: Asteroidea) are a key component of Southern Ocean benthos, with 16% of the known sea star species living there. In temperate marine environments, sea stars commonly play an important role in food webs, acting as keystone species. However, trophic ecology and functional role of Southern Ocean sea stars are still poorly known, notably due to the scarcity of large-scale studies. Here, we report 24,332 trophic marker (stable isotopes and elemental contents of C, N, and S of tegument and/or tube feet) and biometric (arm length, disk radius, arm to disk ratio) measurements in 2,456 specimens of sea stars. Samples were collected between 12 January 1985 and 8 October 2017 in numerous locations along the Antarctic littoral and subantarctic islands. The spatial scope of the data set covers a significant portion of the Southern Ocean (47.717° S to 86.273° S; 127.767° W to 162.201° E; depth, 6–5,338 m). The data set contains 133 distinct taxa, including 72 currently accepted species spanning 51 genera, 20 families, and multiple feeding guilds/functional groups (suspension feeders, sediment feeders, omnivores, predators of mobile or sessile prey). For 505 specimens, mitochondrial CO1 genes were sequenced to confirm and/or refine taxonomic identifications, and those sequences are already publicly available through the Barcode of Life Data System. This number will grow in the future, as molecular analyses are still in progress. Overall, thanks to its large taxonomic, spatial, and temporal extent, as well as its integrative nature (combining genetic, morphological, and ecological data), this data set can be of wide interest to Southern Ocean ecologists, invertebrate zoologists, benthic ecologists, and environmental managers dealing with associated areas. Please cite this data paper in research products derived from the data set, which is freely available without copyright restrictions.

Abstract: The state-of-stress within subducting oceanic plates controls rupture processes of deep intraslab earthquakes. However, little is known about how the large-scale plate geometry and the stress regime relate to the physical nature of the deep intraslab earthquakes. Here we find, by using globally and locally observed seismic records, that the moment magnitude 7.3 2021 East Cape, New Zealand earthquake was driven by a combination of shallow trench-normal extension and unexpectedly, deep trench-parallel compression. We find multiple rupture episodes comprising a mixture of reverse, strike-slip, and normal faulting. Reverse faulting due to the trench-parallel compression is unexpected given the apparent subduction direction, so we require a differential buoyancy-driven stress rotation, which contorts the slab near the edge of the Hikurangi plateau. Our finding highlights that buoyant features in subducting plates may cause diverse rupture behavior of intraslab earthquakes due to the resulting heterogeneous stress state within slabs.

Abstract: Water isotope ratios of ice cores are a key source of information on past temperatures. Through fractionation within the hydrological cycle, temperature is imprinted in the water isotopic composition of snowfalls. However, this signal of climatic interest is modified after deposition when snow remains at the surface exposed to the atmosphere. Comparing time series of surface snow isotopic composition at Dome C with satellite observations of surface snow metamorphism, we found that long summer periods without precipitation favor surface snow metamorphism altering the surface snow isotopic composition. Using excess parameters (combining D,17O, and 18O fractions) allow the identification of this alteration caused by sublimation and condensation of surface hoar. The combined measurement of all three isotopic compositions could help identifying ice core sections influenced by snow metamorphism in sites with very low snow accumulation.