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:

Abstract. Sea-level observations provide information on a variety of processes occurring over different temporal and spatial scales that may contribute to coastal flooding and hazards. However, global research on sea-level extremes is restricted to hourly datasets, which prevent the quantification and analyses of processes occurring at timescales between a few minutes and a few hours. These shorter-period processes, like seiches, meteotsunamis, infragravity and coastal waves, may even dominate in low tidal basins. Therefore, a new global 1 min sea-level dataset – MISELA (Minute Sea-Level Analysis) – has been developed, encompassing quality-checked records of nonseismic sea-level oscillations at tsunami timescales (T<2 h) obtained from 331 tide-gauge sites (https://doi.org/10.14284/456, Zemunik et al., 2021b). This paper describes data quality control procedures applied to the MISELA dataset, world and regional coverage of tide-gauge sites, and lengths of time series. The dataset is appropriate for global, regional or local research of atmospherically induced high-frequency sea-level oscillations, which should be included in the overall sea-level extremes assessments.

Abstract: Aim Delimiting recently diverged species is challenging. During speciation, genetic differentiation may be distributed unevenly across the genome, as different genomic regions can be subject to different selective pressures and evolutionary histories. Reliance on limited numbers of genetic markers that may be underpowered can make species delimitation even more challenging, potentially resulting in taxonomic inconsistencies. Rockhopper penguins of the genus Eudyptes comprise three broadly recognized taxa: northern (E. moseleyi), southern (E. chrysocome) and eastern rockhopper (E. filholi). Their taxonomic status has been controversial for decades, with researchers disagreeing about whether E. chrysocome and E. filholi are distinct species or conspecific. Our goal is to evaluate genome-wide patterns of divergence to evaluate genetic differentiation and species delimitation in rockhopper penguins, and to assess which mechanisms may underlie previous discordance among nuclear versus mitochondrial analyses. Location Sub-Antarctic and temperate coastal regions of the Southern Hemisphere. Methods We generated reduced-representation genomic libraries using double digest restriction-site associated DNA (ddRAD) sequencing to evaluate genetic differentiation, contemporary migration rates and admixture among colonies of rockhopper penguins. Results The extent of genetic differentiation among the three taxa was consistently higher than population-level genetic differentiation found within these and other penguin species. There was no evidence of admixture among the three taxa, suggesting the absence of ongoing gene flow among them. Species delimitation analyses based on molecular data, along with other lines of evidence, provide strong support for the taxonomic distinction of three species of rockhopper penguins. Main conclusions Our results provide strong support for the existence of three distinct species of rockhopper penguins. The recognition of this taxonomic diversity is crucial for the management and conservation of this widely distributed species group. This study illustrates that widespread dispersive seabird lineages lacking obvious morphological differences may nevertheless have complex evolutionary histories and comprise cryptic species diversity.

Abstract: Each winter, the North Atlantic Ocean is the stage for numerous cyclones, the most severe ones leading to seabird mass-mortality events called “winter wrecks.”1, 2, 3 During these, thousands of emaciated seabird carcasses are washed ashore along European and North American coasts. Winter cyclones can therefore shape seabird population dynamics4,5 by affecting survival rates as well as the body condition of surviving individuals and thus their future reproduction. However, most often the geographic origins of impacted seabirds and the causes of their deaths remain unclear.6 We performed the first ocean-basin scale assessment of cyclone exposure in a seabird community by coupling winter tracking data for ∼1,500 individuals of five key North Atlantic seabird species (Alle alle, Fratercula arctica, Uria aalge, Uria lomvia, and Rissa tridactyla) and cyclone locations. We then explored the energetic consequences of different cyclonic conditions using a mechanistic bioenergetics model7 and tested the hypothesis that cyclones dramatically increase seabird energy requirements. We demonstrated that cyclones of high intensity impacted birds from all studied species and breeding colonies during winter but especially those aggregating in the Labrador Sea, the Davis Strait, the surroundings of Iceland, and the Barents Sea. Our broad-scale analyses suggested that cyclonic conditions do not increase seabird energy requirements, implying that they die because of the unavailability of their prey and/or their inability to feed during cyclones. Our study provides essential information on seabird cyclone exposure in a context of marked cyclone regime changes due to global warming.8

Abstract: Tides are one of the basic types of ocean water motions. Previous studies have reported that the M2 constituent exhibits seasonal variations (annual cycles) in some regions. However, based on the newly proposed method of modified two-step harmonic analysis (HA) and its application at 240 global tide gauges, we find that the M2 constituent as well as the S2 and K1 do not have significant seasonal variations at these tide gauges. The seasonal variations of the M2 constituent reported in previous studies are caused by its satellites, the H1 and H2 constituents, which are not resolved in these studies due to the short time window (one month or three months) used in HA. Because the frequency of the H1 (H2) constituent is equal to that of the M2 minus (plus) the frequency of annual cycles, the superposition of the M2, H1 and H2 constituents with constant amplitudes is equivalent to the M2 constituent with seasonally varying amplitudes. Compared with the new method, some adaptations to traditional HA aiming to capture variations in amplitudes and phase lags of constituents have some limitations, because they either neglect some satellites of the major constitutes or introduce spurious fluctuations resulting from an unreasonably large number of independent points. The nodal modulations of predominant constituents are also explored in this study. On the global scale, the nodal modulations of the M2, K1 and O1 constituents agree with the theoretical predictions, except a cold spot region with reduced nodal modulation in the Gulf of Maine and a hot spot region with enhanced nodal modulation in the South China Sea for the M2. Nodal modulation is also found for the S2 constituent (in theory, the S2 has no nodal modulation), which is 0.8% averaged at 164 tide gauges where the S2 is not too weak.