{"id":859162,"date":"2026-06-11T05:23:15","date_gmt":"2026-06-11T05:23:15","guid":{"rendered":"https:\/\/www.europesays.com\/us\/859162\/"},"modified":"2026-06-11T05:23:15","modified_gmt":"2026-06-11T05:23:15","slug":"amplified-arctic-iceberg-traffic-reshapes-benthic-biodiversity","status":"publish","type":"post","link":"https:\/\/www.europesays.com\/us\/859162\/","title":{"rendered":"Amplified Arctic iceberg traffic reshapes benthic biodiversity"},"content":{"rendered":"<p>Iceberg imaging and sampling<\/p>\n<p>During expedition PS126 of RV Polarstern, the opportunity arose to visit a large iceberg in the vicinity of HAUSGARTEN carrying numerous dark stones. The iceberg was accessed by helicopter from the ship on 14 June 2021 at the coordinates 78\u00b0 35.66\u2032\u2009N, 3\u00b0 32.92\u2032\u2009W, allowing us to collect samples of the transported stone material.<\/p>\n<p>A series of overlapping, downward-facing images of the iceberg and its stone load were recorded from the helicopter at an altitude of about 100\u2009m using a Canon EOS digital camera with a lens operated at about 50\u2009mm. The images were stitched together using Agisoft Photoscan Pro. Stone piles placed approximately 5\u2009m apart by field personnel served as a size reference for both individual images and the final mosaic (Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">2a<\/a>). Additional close-range images of stone clusters were acquired on the iceberg using a Nikon D850 with a 50\u2009mm lens, with a ruler included in each frame for scale. The planar area of 20 randomly selected stones per image (n\u2009=\u2009500 in total) was measured using ImageJ. Furthermore, stones were sampled to determine their mineralogical composition. Representative specimens spanning the observed morphological range were collected for classification (Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">2c<\/a>).<\/p>\n<p>An ice core (9\u2009cm in diameter) was collected from the upper 2\u2009m of the iceberg (after snow removal) using a hand auger with an electric drill (Kovacs Enterprise). The coring location was about 50\u2009m away from the nearest stone pile to avoid stone interference during drilling. Three replicate measurements of electrical conductivity were made from sub-cores 27\u201357\u2009cm long (n\u2009=\u20096) to determine the salinity of the ice. Identical methods were used to assess the salinity of two sea ice floes nearby for reference (78\u00b0 54.11\u2032\u2009N, 3\u00b0 9.39\u2032\u2009W and 79\u00b0 01.44\u2032\u2009N, 5\u00b0 42.49\u2032\u2009W).<\/p>\n<p>Seafloor observations<\/p>\n<p>Images were recorded from the seafloor at HAUSGARTEN station EG-IV (2,500\u2009m depth) using either the OFOS or the Ocean Floor Observation and Bathymetry System between 2015 and 2021 (ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 49\" title=\"Purser, A. et al. Ocean floor observation and bathymetry system (OFOBS): a new towed camera\/sonar system for deep-sea habitat surveys. IEEE J. Ocean. Eng. 44, 87&#x2013;99 (2019).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR49\" id=\"ref-link-section-d58876969e1351\" rel=\"nofollow noopener\" target=\"_blank\">49<\/a>). The exact transect location varied by year due to ice conditions (Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig7\" rel=\"nofollow noopener\" target=\"_blank\">1b<\/a>). In 2015 and 2017, the same transect was surveyed, allowing for a comparison between years in an area of seafloor about 5,400\u2009m2. OFOS included a downward-facing camera (Canon EOS), strobes that flashed to illuminate each image and downward-facing laser points that served as a size scale for the image and stone plan area. Telemetry data showed distance to the seafloor, and all images were recorded at a target altitude of 1.5\u2009m. OFOS images were recorded automatically every 30\u2009s while the ship transited at 0.5\u2009knots. Images that were too bright, too dark or at an anomalous altitude (that is, outside the range of 1.3\u2013 1.6\u2009m) were considered ineligible for analysis.<\/p>\n<p>We observed glacial dropstones on the seafloor in OFOS images, ranging from 0.6\u2009cm to 73\u2009cm. The vast majority of stones on the deep seafloor north of 45\u00b0\u2009N are glacial in origin<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 50\" title=\"Kidd, R. B. &amp; Huggett, Q. J. Rock debris on abyssal plains in the Northeast Atlantic: a comparison of epibenthic sledge hauls and photographic surveys. Oceanol. Acta 4, 99&#x2013;104 (1981).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR50\" id=\"ref-link-section-d58876969e1363\" rel=\"nofollow noopener\" target=\"_blank\">50<\/a>. Near-bottom currents in the bathyal Fram Strait are &lt;10\u2009cm\u2009s\u22121, not fast enough to uncover buried stones or bedrock<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 14\" title=\"Soltwedel, T. et al. Natural variability or anthropogenically-induced variation? Insights from 15 years of multidisciplinary observations at the Arctic open-ocean LTER site HAUSGARTEN. Ecol. Indic. 65, 89&#x2013;102 (2016).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR14\" id=\"ref-link-section-d58876969e1369\" rel=\"nofollow noopener\" target=\"_blank\">14<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 21\" title=\"Hasemann, C. et al. Effects of dropstone-induced habitat heterogeneity on Arctic deep-sea benthos with special reference to nematode communities. Mar. Biol. Res. 9, 229&#x2013;245 (2012).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR21\" id=\"ref-link-section-d58876969e1372\" rel=\"nofollow noopener\" target=\"_blank\">21<\/a>.<\/p>\n<p>An exception occurs at a steep reef in the eastern Fram Strait, where stronger near-bottom currents expose bedrock and mobilize locally derived stones<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 51\" title=\"Meyer, K. S., Soltwedel, T. &amp; Bergmann, M. High biodiversity on a deep-water reef in the eastern Fram Strait. PLoS One 9, e105424 (2014).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR51\" id=\"ref-link-section-d58876969e1379\" rel=\"nofollow noopener\" target=\"_blank\">51<\/a>. However, this feature lies several hundred kilometres from station EG-IV, where our imagery was collected. The stones observed on both the seafloor and the visited iceberg consist of shale and quartz (Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>). Taken together, the geological, bathymetric and oceanographic context indicates that stones observed on the seafloor at station EG-IV represent glacial dropstones.<\/p>\n<p>The densities of stones in 30 randomly selected images from the beginning, middle and end of the transect were calculated each year (total of 90 images per year). Dropstones and dropstone-associated fauna were annotated in each image using the online tool BIIGLE (Bio-Image Indexing and Graphical Labelling Environment; <a href=\"http:\/\/biigle.de\" rel=\"nofollow noopener\" target=\"_blank\">biigle.de<\/a>; ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 52\" title=\"Langenk&#xE4;mper, D., Zurowietz, M., Schoening, T. &amp; Nattkemper, T. W. BIIGLE 2.0 - Browsing and annotating large marine image collections. Front. Mar. Sci. 4, 83 (2017).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR52\" id=\"ref-link-section-d58876969e1396\" rel=\"nofollow noopener\" target=\"_blank\">52<\/a>). OFOS images recorded at 3\u2009m altitude reliably show seafloor features and fauna &gt;1\u2009cm across<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 49\" title=\"Purser, A. et al. Ocean floor observation and bathymetry system (OFOBS): a new towed camera\/sonar system for deep-sea habitat surveys. IEEE J. Ocean. Eng. 44, 87&#x2013;99 (2019).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR49\" id=\"ref-link-section-d58876969e1400\" rel=\"nofollow noopener\" target=\"_blank\">49<\/a>. Images in this study were recorded at about 1.5\u2009m altitude, allowing for visualization of fauna and features to 0.6\u2009cm across. Fauna visible in our OFOS images included sponges, bryozoans, tunicates, anemones, soft corals, serpulid polychaetes, barnacles, sea stars and crinoids. Some taxa could be identified based on taxonomic voucher samples collected in previous studies throughout Fram Strait<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 53\" title=\"Bergmann, M., Dannheim, J., Bauerfeind, E. &amp; Klages, M. Trophic relationships along a bathymetric gradient at the deep-sea observatory HAUSGARTEN. Deep Sea Res. I Oceanogr. Res. Pap. 56, 408&#x2013;424 (2009).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR53\" id=\"ref-link-section-d58876969e1404\" rel=\"nofollow noopener\" target=\"_blank\">53<\/a>. Taxa for which no voucher had been identified were given morphotype descriptors. The densities of stones and of each dropstone-associated species in 2015 and 2017 were calculated by dividing the number in each image by the planar area of that image. Plan areas of 250 stones in randomly selected 2015 and 2017 images (n\u2009=\u2009250 per year) were measured using the rectangle select tool in BIIGLE. Plan areas of 500 stones from randomly selected 2021 images were measured using the straight line tool for measurement of length and width in ImageJ.<\/p>\n<p>Dropstones are typically deposited on the seafloor in clusters, and slight variations in the swing of the OFOS camera could cause these clusters to appear in one year but not another. To avoid biasing the results by a small number of images containing dense stone clusters, we removed outlier images with anomalously high stone densities (&gt;5 stones per\u00a0m2) before analysis (n\u2009=\u20091 image, 2015; n\u2009=\u20093 images, 2017). Differences in univariate metrics (that is, stone density, fauna density) were evaluated using parametric analysis of variance (ANOVA) tests or non-parametric Kruskal\u2013Wallis tests, when the assumptions of ANOVA were violated. Chi-square tests were used to test for differences in the size distributions of stones on the iceberg and on the seafloor and between years (2015 and 2017). Statistical analyses were conducted in the R environment using the packages vegan, pairwiseAdonis and ggplot2.<\/p>\n<p>Ship-based glacial ice observations<\/p>\n<p>Small icebergs, growlers and bergy bits generally evade detection by satellite imagery, necessitating visual observations for monitoring their presence. On board the research vessel Polarstern, iceberg observations have been routinely carried out every 3\u2009h by the German Weather Service (DWD). These observations are part of the Surface Synoptic Observations, a globally standardized meteorological data collection program maintained according to the guidelines set by the World Meteorological Organization (WMO). This continuous observational practice has been operational since the commissioning of Polarstern in 1984, providing essential information for onboard activities such as navigation, route planning and aircraft operations, while simultaneously contributing to global meteorological networks.<\/p>\n<p>Synoptic observations from ships include key atmospheric parameters, current weather conditions and basic oceanic parameters, such as air and water temperature, air pressure, wind speed and direction, humidity, precipitation, cloud cover and cloud types, visibility, wave height and periodicity, all encoded in a globally standardized reporting format \u2018FM 13 SHIP\u2019<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 54\" title=\"World Meteorological Organization Information System. Manual on Codes Vol 1.2 (WMO-No. 306) (WMO, 2019&#x2013;2025).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR54\" id=\"ref-link-section-d58876969e1440\" rel=\"nofollow noopener\" target=\"_blank\">54<\/a>. In ice-covered waters, observations include details on the ice-going ability of the vessel, the type and developmental stage of encountered sea ice, and sightings of glacial ice classified according to WMO<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 54\" title=\"World Meteorological Organization Information System. Manual on Codes Vol 1.2 (WMO-No. 306) (WMO, 2019&#x2013;2025).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR54\" id=\"ref-link-section-d58876969e1444\" rel=\"nofollow noopener\" target=\"_blank\">54<\/a> alphanumeric table 0439\/BUFR table 020035. For glacial ice, this code specifically describes the type (icebergs, bergy bits and growlers) and quantity of observed glacial ice, reported in increments of five (see legend in Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig11\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>). The absence of glacial ice and the inability to perform observations due to restricted visibility conditions are also explicitly noted.<\/p>\n<p>The processing and analysis of glacial ice observations were conducted as follows. Initially, all synoptic records collected between 1981 and 2024 were compiled separately by campaign (data for individual campaigns are accessible through PANGAEA<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 55\" title=\"Schmith&#xFC;sen, J. &amp; K&#xF6;nig-Lango, G. Meteorological Observations During POLARSTERN Cruises (ANT-1\/1 et seq). Dataset publication series (PANGAEA, 2025).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR55\" id=\"ref-link-section-d58876969e1454\" rel=\"nofollow noopener\" target=\"_blank\">55<\/a>). Subsequently, all valid glacial ice observations within the Fram Strait region (70\u00b0\u201385\u00b0\u2009N, 30\u00b0\u2009W\u201330\u00b0\u2009E) were filtered. Observations were defined as valid if visibility exceeded 1\u2009km. This resulted in approximately 8,200 valid observations, with glacial ice reported around the vessel in about 32% of these cases. The spatial distribution of observations is shown in Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3a<\/a>. Note that repeated sightings at the same position or within a radius smaller than the visibility at the time of observation were excluded to avoid double-counting.<\/p>\n<p>Following this filtering, campaign-specific and annual mean frequencies of iceberg sightings were calculated (Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3b<\/a>), as well as frequency distributions of individual glacial ice types. Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig11\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a> shows the frequencies of glacial ice observations, along with the corresponding classification codes. Most frequently observed were 1\u20135 icebergs accompanied by growlers and bergy bits (class 6), followed by occurrences of only icebergs (1\u20135, class 1) or only growlers and bergy bits (class 4).<\/p>\n<p>Visual glacial ice observations inherently involve uncertainties affecting their representativeness. First, sightings may vary among individual observers, particularly under challenging visibility conditions or low sun angles. Moreover, observer performance might evolve over time, for example, as observers gain experience. Although this subjective error is difficult to quantify precisely, it is considered small because of the consistent use of trained meteorological experts. Second, weather conditions such as precipitation or fog can limit the number of visual observations carried out during expeditions. On average, about 20.4% of planned observations are hindered by visibility conditions below the 1\u2009km threshold, although this percentage can be considerably higher for individual expeditions. However, analysis of the 40-year observational record shows no significant trend in visibility conditions, indicating that no systematic long-term biases are present in the observational time series (Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig12\" rel=\"nofollow noopener\" target=\"_blank\">6<\/a>). Third, larger and more prominent icebergs are proportionally detected more frequently than smaller growlers and bergy bits. This uncertainty also remains difficult to precisely quantify.<\/p>\n<p>Given these limitations, definitive quantification of absolute iceberg frequencies, interannual variability or spatial distribution patterns remains uncertain and should be interpreted cautiously. Nevertheless, Fram Strait, and specifically the HAUSGARTEN region, benefits from frequent observation efforts, on average traversed by the AWI 2.5 times per year, yielding a continuous and extensive dataset. The density and duration of this observational time series enable robust detection of relative changes in glacial ice occurrence despite the discussed observational constraints.<\/p>\n<p>Backtracking of glacial ice<\/p>\n<p>To identify the source regions (calving sites) of glacial ice and, consequently, the origin of stones deposited on the deep seafloor, we backtracked ship-based glacial ice observations using satellite-derived ice drift data. For this purpose, we used IceTrack<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 56\" title=\"Krumpen, T. et al. Smoother sea ice with fewer pressure ridges in a more dynamic Arctic. Nat. Clim. Chang. 15, 66&#x2013;72 (2025).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR56\" id=\"ref-link-section-d58876969e1488\" rel=\"nofollow noopener\" target=\"_blank\">56<\/a>, a tool originally developed for backtracking sea ice to the sites where it was formed.<\/p>\n<p>For the backtracking approach, ice is traced backwards on a daily basis from observed locations using three independent low-resolution sea ice drift products: (1) OSI-405-c ice-motion data from the Ocean and Sea Ice (OSI) Satellite Application Facility<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 57\" title=\"Lavergne, T. &amp; Down, E. Validation and Monitoring of the OSI SAF Low Resolution Sea Ice Drift Product (v.5). Technical Report No. SAF\/OSI\/CDOP\/Met.no\/T&amp;V\/RP\/131. The EUMETSAT Network of Satellite Application Facilities (OSI SAF, 2016).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR57\" id=\"ref-link-section-d58876969e1495\" rel=\"nofollow noopener\" target=\"_blank\">57<\/a> (SAF); (2) MERGED ice-motion vectors provided by the Center for Satellite Exploitation and Research<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 58\" title=\"Girard-Ardhuin, F. &amp; Ezraty, R. Enhanced Arctic Sea ice drift estimation merging radiometer and scatterometer data. IEEE Trans. Geosci. Remote Sens. 50, 2639&#x2013;2648 (2012).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR58\" id=\"ref-link-section-d58876969e1499\" rel=\"nofollow noopener\" target=\"_blank\">58<\/a> (CERSAT); and (3) Polar Pathfinder daily motion vectors (v.4.1) available from the National Snow and Ice Data Center<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 59\" title=\"Tschudi, M. A., Meier, W. N. &amp; Stewart, J. S. An enhancement to sea ice motion and age products at the National Snow and Ice Data Center (NSIDC). Cryosphere 14, 1519&#x2013;1536 (2020).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR59\" id=\"ref-link-section-d58876969e1503\" rel=\"nofollow noopener\" target=\"_blank\">59<\/a> (NSIDC). Backtracking is terminated whenever trajectories intersect coastlines, encounter fast-ice edges, or when the satellite-derived ice concentration (from OSI SAF products OSI-430-a and OSI-450-a) falls below 50%.<\/p>\n<p>This Lagrangian tracking method is applicable to glacial ice as long as icebergs, growlers and bergy bits are firmly embedded within the pack ice and carried along by it. Under looser ice-pack conditions or for deep-drafted icebergs, the validity of this assumption diminishes. Here, ocean currents begin to strongly influence the drift, causing icebergs to diverge from the general motion of the surrounding, mainly wind-driven, sea ice<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 60\" title=\"Marson, J. M., Myers, P. G., Garbo, A., Copland, L. &amp; Mueller, D. Sea ice-driven iceberg drift in Baffin Bay. J. Geophys. Res. Oceans 129, e2023JC020697 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR60\" id=\"ref-link-section-d58876969e1510\" rel=\"nofollow noopener\" target=\"_blank\">60<\/a>. To account for this limitation, we restricted backtracking to sightings of glacial ice firmly entrained within the pack ice at the time of observation and throughout their drift. Specifically, we imposed a criterion that satellite-derived ice concentration along each trajectory should not fall below 80% for more than 30 cumulative days. A total of 106 out of the 1,362 glacial ice sightings fulfilled this condition, with their derived drift tracks shown in Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3a<\/a>.<\/p>\n<p>The applied concentration criterion is grounded in an evaluation of the tracking approach against independent drift observations. First, trajectories of 50 GPS-equipped sea ice buoys deployed in the Arctic Ocean between 2012 and 2016 (Meereisportal.de data repository<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 61\" title=\"Grosfeld, K. et al. Online sea-ice knowledge and data platform. Polarforschung 85, 143&#x2013;155 (2016).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR61\" id=\"ref-link-section-d58876969e1520\" rel=\"nofollow noopener\" target=\"_blank\">61<\/a>) were reconstructed using IceTrack and directly compared with their observed GPS positions. The mean great-circle deviation between reconstructed and observed locations amounts to approximately 47\u2009\u00b1\u200948\u2009km after 100\u2009days and 109\u2009\u00b1\u200968\u2009km after 300\u2009days (Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig16\" rel=\"nofollow noopener\" target=\"_blank\">10a<\/a>), indicating reliable reproduction of large-scale drift pathways under compact pack-ice conditions. To assess tracking performance under looser ice conditions, iceberg drift was additionally evaluated using a small set of GPS-equipped icebergs that were partially exposed to open-water conditions. The applied dataset comprises 11 icebergs tracked between 2012 and 2025 in the Beaufort Gyre and north of the Canadian Arctic Archipelago (data sources: International Arctic Buoy Programme Ice Island Archive, <a href=\"https:\/\/iabp.apl.uw.edu\/Ice_Islands_2025.html\" rel=\"nofollow noopener\" target=\"_blank\">https:\/\/iabp.apl.uw.edu\/Ice_Islands_2025.html<\/a>; Carleton University Ice Island Drift Tracking Database, <a href=\"https:\/\/wirl.carleton.ca\/research\/ice\/ice-islands\/ibtd\/\" rel=\"nofollow noopener\" target=\"_blank\">https:\/\/wirl.carleton.ca\/research\/ice\/ice-islands\/ibtd\/<\/a>). The comparison between reconstructed and observed trajectories indicates larger deviations (55\u2009\u00b1\u200970\u2009km after 100\u2009days and 286\u2009\u00b1\u2009123\u2009km after 300\u2009days), particularly during periods when sea ice concentration falls below 80% (Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig16\" rel=\"nofollow noopener\" target=\"_blank\">10b<\/a>). Despite the higher deviations under loose ice conditions, positional deviations remain sufficiently small to distinguish iceberg source regions at basin and sector scale, although attribution to individual marine-terminating glaciers is not supported, particularly for trajectories involving extended drift durations.<\/p>\n<p>To assess whether the increasing number of icebergs observed in Fram Strait could, in part, be explained by an enhanced input of glacial ice from the Canadian Arctic, we applied the same Lagrangian tracking framework. Sea ice crossing a virtual gate at 83\u00b0\u2009N (25\u00b0\u2009W\u201325\u00b0\u2009E) was traced backwards in time using an identical Lagrangian backtracking configuration to that described in ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 45\" title=\"Krumpen, T. et al. Arctic warming interrupts the Transpolar Drift and affects long-range transport of sea ice and ice-rafted matter. Sci. Rep. 9, 5459 (2019).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR45\" id=\"ref-link-section-d58876969e1548\" rel=\"nofollow noopener\" target=\"_blank\">45<\/a>. The Beaufort fraction was defined as the proportion of tracked trajectories that intersected the Beaufort Gyre region prior to export through Fram Strait (Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig15\" rel=\"nofollow noopener\" target=\"_blank\">9<\/a>). Temporal changes in this fraction were analysed to evaluate whether the contribution of Beaufort-origin ice increased after 2000.<\/p>\n<p>Glacier changes<\/p>\n<p>In this study, we relate the observed increase in iceberg frequency to changes at the contributing marine-terminating glaciers. These changes can be quantified through several complementary approaches that differ in physical meaning, temporal resolution and data requirements. Although the main outlet glaciers in northeast Greenland are comparatively well observed, many marine-terminating glaciers in the Russian High Arctic lack the consistent, long-term datasets required for detailed dynamic assessments<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 62\" title=\"Moholdt, G., Wouters, B. &amp; Gardner, A. S. Recent mass changes of glaciers in the Russian high Arctic. Geophys. Res. Lett. 39, L10502 (2012).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR62\" id=\"ref-link-section-d58876969e1564\" rel=\"nofollow noopener\" target=\"_blank\">62<\/a>.<\/p>\n<p>Frontal ablation<\/p>\n<p>Calving is most directly reflected in changes in frontal ablation. Estimates of frontal ablation rely on high-resolution satellite imagery used to track terminus positions over time, typically through manual or semi-automated delineation. For northeast Greenland, these observations document the extensive retreat and fragmentation of the ZI and NG glaciers since the early 2000s (Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig10\" rel=\"nofollow noopener\" target=\"_blank\">4<\/a>, derived from Landsat 5\u20138 optical imagery), providing an indication of enhanced calving over the past two decades. By contrast, for many Russian High Arctic glaciers, terminus positions have been mapped only episodically or over relatively short time windows<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 35\" title=\"Zheng, W. et al. Accelerating glacier mass loss on Franz Josef Land, Russian Arctic. Remote Sens. Environ. 211, 357&#x2013;375 (2018).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR35\" id=\"ref-link-section-d58876969e1578\" rel=\"nofollow noopener\" target=\"_blank\">35<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 39\" title=\"Wytiahlowsky, H., Stokes, C. R. &amp; Evans, D. J. A. Remote sensing of glacier change (1965-2021) and identification of surge-type glaciers on Severnaya Zemlya, Russian high Arctic. J. Glaciol. 69, 1764&#x2013;1784 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR39\" id=\"ref-link-section-d58876969e1581\" rel=\"nofollow noopener\" target=\"_blank\">39<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 63\" title=\"Kochtitzky, W. et al. The unquantified mass loss of Northern Hemisphere marine-terminating glaciers from 2000-2020. Nat. Commun. 13, 5835 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR63\" id=\"ref-link-section-d58876969e1584\" rel=\"nofollow noopener\" target=\"_blank\">63<\/a>.<\/p>\n<p>Ice discharge<\/p>\n<p>A second approach quantifies ice discharge, defined as the flux of grounded ice across a gate near the grounding line. Discharge isolates the dynamic component of glacier change and, over multi-decadal timescales, typically correlates with changes in calving activity. Although discharge is not identical to frontal ablation, increased discharge generally reflects accelerating flow and dynamic thinning that favour enhanced calving. This method requires ice velocity, ice thickness and bedrock topography and is therefore available only on longer timescales for the major outlet glaciers in northeast Greenland. For ZI and NG, we use the discharge estimates in ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 64\" title=\"Mankoff, K. D. et al. Greenland Ice Sheet solid ice discharge from 1986 through March 2020. Earth Syst. Sci. Data 12, 1367&#x2013;1383 (2020).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR64\" id=\"ref-link-section-d58876969e1596\" rel=\"nofollow noopener\" target=\"_blank\">64<\/a>, which extend back to 1987 and allow quantification of long-term changes in dynamic mass loss. In Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3c<\/a>, cumulative discharge estimates relative to the stable period from 1985 to 2000 are shown, providing a robust baseline against which the post-2000 increase can be evaluated.<\/p>\n<p>Mass-change estimates<\/p>\n<p>Where neither frontal-ablation nor discharge records exist, glacier change can be assessed through mass-change estimates derived from remote-sensing products. These measurements reflect the combined effects of surface mass balance and dynamic thinning. For northeast Greenland, we calculate mass changes using NASA Operation IceBridge ATM data and satellite altimetry (ICESat, ICESat-2, CryoSat-2 and ENVISAT), referenced to a 1978 DEM<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 32\" title=\"Khan, S. A. et al. Extensive inland thinning and speed-up of Northeast Greenland ice stream. Nature 611, 727&#x2013;732 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR32\" id=\"ref-link-section-d58876969e1611\" rel=\"nofollow noopener\" target=\"_blank\">32<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 65\" title=\"Korsgaard, N. J. et al. Digital elevation model and orthophotographs of Greenland based on aerial photographs from 1978-1987. Sci. Data 3, 160032 (2016).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR65\" id=\"ref-link-section-d58876969e1614\" rel=\"nofollow noopener\" target=\"_blank\">65<\/a>. Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig13\" rel=\"nofollow noopener\" target=\"_blank\">7<\/a> presents the time series of mass loss for the lower sectors of ZI and NG<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 32\" title=\"Khan, S. A. et al. Extensive inland thinning and speed-up of Northeast Greenland ice stream. Nature 611, 727&#x2013;732 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR32\" id=\"ref-link-section-d58876969e1621\" rel=\"nofollow noopener\" target=\"_blank\">32<\/a>. The two mass-change time series indicate that there was no significant ice loss between 1978 and 2000, an assumption that is also applied when calculating the cumulative discharge shown in Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3c<\/a>. For SZ and FJL, we rely on published mass-change assessments based on DEM differencing, radar and optical altimetry<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 36\" title=\"Tepes, P. et al. Changes in elevation and mass of Arctic glaciers and ice caps, 2010&#x2013;2017. Remote Sens. Environ. 261, 112481 (2021).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR36\" id=\"ref-link-section-d58876969e1628\" rel=\"nofollow noopener\" target=\"_blank\">36<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 38\" title=\"Sommer, C., Seehaus, T., Glazovsky, A. &amp; Braun, M. H. Brief communication: increased glacier mass loss in the Russian High Arctic (2010-2017). Cryosphere 16, 35&#x2013;42 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR38\" id=\"ref-link-section-d58876969e1631\" rel=\"nofollow noopener\" target=\"_blank\">38<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 39\" title=\"Wytiahlowsky, H., Stokes, C. R. &amp; Evans, D. J. A. Remote sensing of glacier change (1965-2021) and identification of surge-type glaciers on Severnaya Zemlya, Russian high Arctic. J. Glaciol. 69, 1764&#x2013;1784 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR39\" id=\"ref-link-section-d58876969e1634\" rel=\"nofollow noopener\" target=\"_blank\">39<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 62\" title=\"Moholdt, G., Wouters, B. &amp; Gardner, A. S. Recent mass changes of glaciers in the Russian high Arctic. Geophys. Res. Lett. 39, L10502 (2012).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR62\" id=\"ref-link-section-d58876969e1637\" rel=\"nofollow noopener\" target=\"_blank\">62<\/a>. Although these studies differ in methodological foundations, resolution and temporal coverage, they consistently report sustained mass loss of marine-terminating glaciers since 2000.<\/p>\n<p>Surface velocities<\/p>\n<p>Where direct estimates of frontal ablation or ice discharge are unavailable, glacier dynamics can be assessed using satellite-derived surface velocity time series as indicators of changes in ice flow. We use glacier surface velocities for selected marine-terminating glaciers on FJL and SZ from the NASA MEaSUREs ITS_LIVE project, generated using the auto-RIFT feature-tracking algorithm<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 48\" title=\"Gardner, A. S., Fahnestock, M. A. &amp; Scambos, T. A. MEaSUREs ITS_LIVE Regional Glacier and Ice Sheet Surface Velocities (NSIDC-0776, Version 1) [Data Set]. (National Snow and Ice Data Center, 2019).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR48\" id=\"ref-link-section-d58876969e1650\" rel=\"nofollow noopener\" target=\"_blank\">48<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 66\" title=\"Gardner, A. S. et al. Increased West Antarctic and unchanged East Antarctic ice discharge over the last 7 years. Cryosphere 12, 521&#x2013;547 (2018).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR66\" id=\"ref-link-section-d58876969e1653\" rel=\"nofollow noopener\" target=\"_blank\">66<\/a>. Each velocity estimate shown in Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4<\/a> and Extended Data Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig14\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a> represents a temporally averaged surface speed over the image-pair separation interval. Seasonal mean velocities were computed using a moving temporal window centred on 1 July of each year (\u00b130\u2009days). Within each window, velocities were averaged using weights proportional to the image-pair temporal separation to account for unequal temporal coverage.<\/p>\n<p>FESOM2-iceberg drift model<\/p>\n<p>To quantify potential changes in the large-scale transport of icebergs towards Fram Strait, we conducted a controlled sensitivity experiment using the FESOM2<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 67\" title=\"Danilov, S., Sidorenko, D., Wang, Q. &amp; Jung, T. The Finite-volumE Sea ice-Ocean Model (FESOM2). Geosci. Model Dev. 10, 765&#x2013;789 (2017).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR67\" id=\"ref-link-section-d58876969e1672\" rel=\"nofollow noopener\" target=\"_blank\">67<\/a> equipped with an interactive iceberg module<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 68\" title=\"Rackow, T. et al. A simulation of small to giant Antarctic iceberg evolution: differential impact on climatology estimates. J. Geophys. Res. Oceans 122, 3170&#x2013;3190 (2017).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR68\" id=\"ref-link-section-d58876969e1676\" rel=\"nofollow noopener\" target=\"_blank\">68<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 69\" title=\"Ackermann, L. et al. A comprehensive Earth system model (AWI-ESM2.1) with interactive icebergs: effects on surface and deep-ocean characteristics. Geosci. Model Dev. 17, 3279&#x2013;3301 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR69\" id=\"ref-link-section-d58876969e1679\" rel=\"nofollow noopener\" target=\"_blank\">69<\/a>. The experiment prescribes temporally constant calving fluxes from predefined glacier sectors in NEGIS, FJL, and SZ, allowing us to examine changes in iceberg drift dynamics and associated processes along their trajectories over the simulation period.<\/p>\n<p>FESOM2 is a global, multi-resolution ocean-sea ice model that supports regional mesh refinement in areas of interest. In this study, we use a model configuration with a horizontal resolution of approximately 4.5\u2009km in the Arctic Ocean and Nordic Seas. The model is forced with the JRA-55 atmospheric reanalysis<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 70\" title=\"Kobayashi, S. et al. The JRA-55 reanalysis: general specifications and basic characteristics. J. Meteorol. Soc. Jpn. 93, 5&#x2013;48 (2015).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR70\" id=\"ref-link-section-d58876969e1686\" rel=\"nofollow noopener\" target=\"_blank\">70<\/a> and has previously been evaluated against observational datasets for Arctic Ocean circulation and sea ice dynamics<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" title=\"Heukamp, F. O., Wekerle, C., Kanzow, T., McPherson, R. &amp; Baumann, T. M. Atlantic water recirculation in the northern Barents Sea affects winter sea ice extent. Nat. Commun. 16, 5148 (2025).\" href=\"#ref-CR71\" id=\"ref-link-section-d58876969e1690\">71<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" title=\"Heukamp, F. O. et al. Cyclones modulate the control of the North Atlantic Oscillation on transports into the Barents Sea. Commun. Earth Environ. 4, 324 (2023).\" href=\"#ref-CR72\" id=\"ref-link-section-d58876969e1690_1\">72<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 73\" title=\"Wekerle, C. et al. Atlantic Water warming increases melt below Northeast Greenland&#x2019;s last floating ice tongue. Nat. Commun. 15, 1336 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR73\" id=\"ref-link-section-d58876969e1693\" rel=\"nofollow noopener\" target=\"_blank\">73<\/a>.<\/p>\n<p>The iceberg module represents icebergs as Lagrangian particles with prescribed geometric properties, including height, length and width. These dimensions evolve dynamically during the simulation through thermodynamic processes such as basal and lateral melting, wave-induced erosion and buoyant convection. The momentum balance of model icebergs includes Coriolis and gravitational forces as well as drag forces exerted by the atmosphere, ocean and sea ice. Sea ice drag follows the step-function parameterization described in ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 74\" title=\"Lichey, C. &amp; Hellmer, H. H. Modeling giant-iceberg drift under the influence of sea ice in the Weddell Sea, Antarctica. J. Glaciol. 47, 452&#x2013;460 (2001).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR74\" id=\"ref-link-section-d58876969e1700\" rel=\"nofollow noopener\" target=\"_blank\">74<\/a>. For sea ice concentrations below 15%, sea ice drag is set to zero. When sea ice concentration exceeds 90%, and sea ice strength is greater than 10,000\u2009N\u2009m\u22123, icebergs are assumed to be frozen into the pack ice and drift with the surrounding sea ice. A detailed description of the iceberg module is provided in refs.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 68\" title=\"Rackow, T. et al. A simulation of small to giant Antarctic iceberg evolution: differential impact on climatology estimates. J. Geophys. Res. Oceans 122, 3170&#x2013;3190 (2017).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR68\" id=\"ref-link-section-d58876969e1706\" rel=\"nofollow noopener\" target=\"_blank\">68<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 69\" title=\"Ackermann, L. et al. A comprehensive Earth system model (AWI-ESM2.1) with interactive icebergs: effects on surface and deep-ocean characteristics. Geosci. Model Dev. 17, 3279&#x2013;3301 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR69\" id=\"ref-link-section-d58876969e1709\" rel=\"nofollow noopener\" target=\"_blank\">69<\/a>.<\/p>\n<p>In the Southern Ocean, simulated iceberg drift, both in open water and while embedded in pack ice, has been evaluated against observed iceberg trajectories derived from satellite<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 75\" title=\"Budge, J. S. &amp; Long, D. G. A comprehensive database for Antarctic iceberg tracking using scatterometer data. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 11, 434&#x2013;442 (2018).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR75\" id=\"ref-link-section-d58876969e1716\" rel=\"nofollow noopener\" target=\"_blank\">75<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 76\" title=\"Stuart, K. M. &amp; Long, D. G. Iceberg size and orientation estimation using SeaWinds. Cold Reg. Sci. Technol. 69, 39&#x2013;51 (2011).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR76\" id=\"ref-link-section-d58876969e1719\" rel=\"nofollow noopener\" target=\"_blank\">76<\/a> and buoy datasets<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 77\" title=\"Schodlok, M. P., Hellmer, H. H., Rohardt, G. &amp; Fahrbach, E. Weddell Sea iceberg drift: five years of observations. J. Geophys. Res. 111, C06018 (2006).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR77\" id=\"ref-link-section-d58876969e1723\" rel=\"nofollow noopener\" target=\"_blank\">77<\/a> in refs. <a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 68\" title=\"Rackow, T. et al. A simulation of small to giant Antarctic iceberg evolution: differential impact on climatology estimates. J. Geophys. Res. Oceans 122, 3170&#x2013;3190 (2017).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR68\" id=\"ref-link-section-d58876969e1727\" rel=\"nofollow noopener\" target=\"_blank\">68<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 69\" title=\"Ackermann, L. et al. A comprehensive Earth system model (AWI-ESM2.1) with interactive icebergs: effects on surface and deep-ocean characteristics. Geosci. Model Dev. 17, 3279&#x2013;3301 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR69\" id=\"ref-link-section-d58876969e1730\" rel=\"nofollow noopener\" target=\"_blank\">69<\/a>. However, in the Arctic, direct validation is difficult because of the limited availability of observational data. Nevertheless, icebergs that become embedded in dense pack ice are expected to drift largely with the surrounding sea ice cover. Because FESOM2 realistically reproduces observed Arctic sea ice drift and variability<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 71\" title=\"Heukamp, F. O., Wekerle, C., Kanzow, T., McPherson, R. &amp; Baumann, T. M. Atlantic water recirculation in the northern Barents Sea affects winter sea ice extent. Nat. Commun. 16, 5148 (2025).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR71\" id=\"ref-link-section-d58876969e1734\" rel=\"nofollow noopener\" target=\"_blank\">71<\/a>, we assume that the drift of icebergs entrained in pack ice is captured realistically by the model.<\/p>\n<p>Whenever the simulated iceberg draft exceeds the local bathymetry, the iceberg is assumed to be grounded. Basal melting continues during grounding and can eventually reduce the draft sufficiently for the iceberg to refloat. Discrete icebergs are generated from an integrated calving flux. The prescribed ice discharge is converted into an equivalent calving area by dividing the flux by the assumed iceberg thickness (250\u2009m; refs.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 78\" title=\"Marson, J. M., Myers, P. G., Hu, X. &amp; Le Sommer, J. Using vertically integrated ocean fields to characterize Greenland icebergs&#x2019; distribution and lifetime. Geophys. Res. Lett. 45, 4208&#x2013;4217 (2018).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR78\" id=\"ref-link-section-d58876969e1742\" rel=\"nofollow noopener\" target=\"_blank\">78<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 79\" title=\"Marsh, R. et al. NEMO-ICB (v1.0): interactive icebergs in the NEMO ocean model globally configured at eddy-permitting resolution. Geosci. Model Dev. 8, 1547&#x2013;1562 (2015).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR79\" id=\"ref-link-section-d58876969e1745\" rel=\"nofollow noopener\" target=\"_blank\">79<\/a>). Individual iceberg sizes are sampled from a \u22122.2 power-law size distribution following ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 80\" title=\"Shiggins, C. J., Lea, J. M. &amp; Brough, S. Automated ArcticDEM iceberg detection tool: insights into area and volume distributions, and their potential application to satellite imagery and modelling of glacier-iceberg-ocean systems. Cryosphere 17, 15&#x2013;32 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR80\" id=\"ref-link-section-d58876969e1749\" rel=\"nofollow noopener\" target=\"_blank\">80<\/a>. Icebergs are assumed to have a quadratic surface and cuboidal geometry. For icebergs smaller than 250\u2009m in horizontal extent, height is set equal to the horizontal dimensions, whereas larger icebergs are assigned a maximum height of 250\u2009m to reduce the risk of instantaneous grounding after release. In an iterative adjustment procedure, iceberg dimensions are adjusted so that the total iceberg volume matches the prescribed integrated discharge.<\/p>\n<p>To isolate transport dynamics from variations in iceberg production, iceberg release was prescribed with temporally constant annual fluxes throughout the simulation period (1970\u20132023), thereby preventing changes in prescribed calving rates from influencing the simulated iceberg distribution. Icebergs were released continuously throughout the year rather than at a fixed time of release to avoid introducing artificial seasonal biases. Release locations were defined in sufficiently deep offshore waters rather than directly at glacier grounding lines, thereby minimizing artificial grounding effects and reducing the influence of poorly constrained near-coastal processes such as fast-ice storage or m\u00e9lange-controlled release.<\/p>\n<p>Icebergs were released in three source sectors identified by the backtracking analysis (Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3a<\/a>): NEG, FJL and SZ. These regions represent the primary sectors from which icebergs may enter Fram Strait. The simulations, therefore, provide a physically consistent large-scale representation of the dominant drift corridors linking these source regions Fram Strait. For the three source regions, the following integrated calving fluxes, averaged over the time period 2000\u20132010, were prescribed:<\/p>\n<ul class=\"u-list-style-bullet\">\n<li>\n<p>SZ: 3.15\u2009gigatonnes per annum (Gt\u2009a\u22121) (ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 63\" title=\"Kochtitzky, W. et al. The unquantified mass loss of Northern Hemisphere marine-terminating glaciers from 2000-2020. Nat. Commun. 13, 5835 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR63\" id=\"ref-link-section-d58876969e1770\" rel=\"nofollow noopener\" target=\"_blank\">63<\/a>)<\/p>\n<\/li>\n<li>\n<p>FJL: 8.78\u2009Gt\u2009a\u22121(ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 63\" title=\"Kochtitzky, W. et al. The unquantified mass loss of Northern Hemisphere marine-terminating glaciers from 2000-2020. Nat. Commun. 13, 5835 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR63\" id=\"ref-link-section-d58876969e1782\" rel=\"nofollow noopener\" target=\"_blank\">63<\/a>)<\/p>\n<\/li>\n<li>\n<p>NEG: 11.15\u2009Gt\u2009a\u22121 (ZI<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 64\" title=\"Mankoff, K. D. et al. Greenland Ice Sheet solid ice discharge from 1986 through March 2020. Earth Syst. Sci. Data 12, 1367&#x2013;1383 (2020).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR64\" id=\"ref-link-section-d58876969e1794\" rel=\"nofollow noopener\" target=\"_blank\">64<\/a>); 8.41\u2009Gt\u2009a\u22121 (NG<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 81\" title=\"Millan, R., Mouginot, J., Rabatel, A. &amp; Morlighem, M. Ice velocity and thickness of the world&#x2019;s glaciers. Nat. Geosci. 15, 124&#x2013;129 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR81\" id=\"ref-link-section-d58876969e1800\" rel=\"nofollow noopener\" target=\"_blank\">81<\/a>); 0.33\u2009Gt\u2009a\u22121 (NEG peripheral glaciers<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 63\" title=\"Kochtitzky, W. et al. The unquantified mass loss of Northern Hemisphere marine-terminating glaciers from 2000-2020. Nat. Commun. 13, 5835 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR63\" id=\"ref-link-section-d58876969e1807\" rel=\"nofollow noopener\" target=\"_blank\">63<\/a>).<\/p>\n<\/li>\n<\/ul>\n<p>The first 16 years (1970\u20131985) of the simulation are considered a spin-up period to reach an equilibrated number of active icebergs, as suggested in ref.\u2009<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 79\" title=\"Marsh, R. et al. NEMO-ICB (v1.0): interactive icebergs in the NEMO ocean model globally configured at eddy-permitting resolution. Geosci. Model Dev. 8, 1547&#x2013;1562 (2015).\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#ref-CR79\" id=\"ref-link-section-d58876969e1816\" rel=\"nofollow noopener\" target=\"_blank\">79<\/a>. To separate the effect of changing ocean and sea ice conditions on iceberg drift from the increase in active icebergs during the spin-up phase, the analysis of iceberg dynamics focuses on the period from 1986 to 2023. To quantify spatial patterns of iceberg transport (Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41586-026-10630-4#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>), we calculate iceberg transit frequency as the relative frequency of iceberg passages through each grid cell over the simulation period.<\/p>\n","protected":false},"excerpt":{"rendered":"Iceberg imaging and sampling During expedition PS126 of RV Polarstern, the opportunity arose to visit a large iceberg&hellip;\n","protected":false},"author":3,"featured_media":859163,"comment_status":"","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"footnotes":"","_share_on_mastodon":"0"},"categories":[23],"tags":[18844,52108,93872,93873,945,746,10046,10047,159,67,132,68],"class_list":["post-859162","post","type-post","status-publish","format-standard","has-post-thumbnail","category-environment","tag-biodiversity","tag-climate-change-ecology","tag-climate-change-impacts","tag-cryospheric-science","tag-ecology","tag-environment","tag-humanities-and-social-sciences","tag-multidisciplinary","tag-science","tag-united-states","tag-unitedstates","tag-us"],"share_on_mastodon":{"url":"","error":""},"_links":{"self":[{"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/posts\/859162","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/users\/3"}],"replies":[{"embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/comments?post=859162"}],"version-history":[{"count":0,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/posts\/859162\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/media\/859163"}],"wp:attachment":[{"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/media?parent=859162"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/categories?post=859162"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/tags?post=859162"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}