Ice Thickness Research Articles

Identifying key safety investments needed for arctic shipping via a fuzzy analytic hierarchy process (FAHP) approach

ABSTRACT The retreating ice coverage and reduced (ice) thickness scientifically recorded in the wider Arctic region, during the last few decades, are transforming it into a more accessible operating environment. Vessels engaged in maritime transport activities can now navigate the already established routes for more extended periods each year. With this background, Arctic shipping routes have attracted the interest of academics and maritime stakeholders because of the potential distance reduction compared to other existing (and often viewed as) conventional routes utilizing the Suez and Panama Canals today. However, the existing Arctic shipping environment requires specific future safety investments. This paper intends to provide input towards better-informed decision-making and identification of priority investments needed to improve the level of safety related to the Northern Sea Route (NSR) and Northwest Passage (NWP). A multi-attribute decision-making methodology, the Fuzzy Analytic Hierarchy Process, is applied to rank the various proposed safety investments. The results indicate that “infrastructure and facility” should come first, followed by “personnel,” “technology,” “measure,” and “management” investments. Ice monitoring and weather forecasting, strengthening Arctic seafaring expertise, and icebreakers stand out as the top three most important future investments required in the region.

A theoretical model of excess attenuation of acoustic signals propagating under cracked sea-ice landscapes.

The Arctic sheet is transitioning from a continuous cover of thick multi-year ice to a fragmented landscape of thin young ice. If the type of acoustic transmission allows repetitive interaction of rays with the sea surface, in the fragmented scenario acoustic rays will undergo a random sequence of reflections from water or sea-ice interfaces. Calm sea conditions in the water channels between the ice floes (leads) and the smooth, flat surface of the young ice bottom reduce scattering due to interface roughness, resulting only in scattering due to inhomogeneity in surface reflectivity. Using an idealized framework, this study investigates the extent to which the mid- to high-frequency underwater acoustic propagation is altered due to repetitive interactions of acoustic signals with a sea surface consisting of a random distribution of ice sheets and leads. An expression for the coherent field (the acoustic field averaged over an ensemble of realizations of sea-ice distributions) was derived from theory. Any deviation from a homogeneous surface condition (either by randomly adding ice slabs in a free ice surface or by including leads in a fully ice-covered sea surface) leads to an excess attenuation of the coherent field. Results are validated by numerical simulations.

Thermodynamic and Dynamic Variations in Sea Ice Thickness of the Ross Sea, Antarctica, Driven by Atmospheric Circulation

AbstractAtmospheric circulation has significant impacts on sea ice drifting patterns and mass balance, as wind drag induces pressure ridges and leads on the sea ice surface. In this study, the spatiotemporal distributions of these dynamic sea ice deformation features in the Ross Sea are examined using ICESat‐2 (IS2) ATL10 freeboard data (2019–2022). The temporal variation of the modal sea ice thickness (SIT), caused by thermodynamic ice growth and sea ice advection, varies from 0.7–1.0 m in April to 1.0–1.6 m in July–September and decreases thereafter in the northwest (NW) and northeast (NE) sectors. This temporal variation of modal SIT agrees with the air temperature (correlation coefficients >0.5). The southwest (SW) sector shows a consistently low modal SIT (<1.0 m) because of the production of new ice in polynyas and continuous northward sea ice drift. Meanwhile, the southeast (SE) sector shows the thickest ice in Octobers 2019 and 2020 because of the advection of thick ice from the Amundsen Sea, which was reduced in 2021 and 2022. In terms of dynamic sea ice deformation, the SE sector shows the largest deformation because of the wind‐driven convergence of sea ice movement. However, such intense deformation in the SE sector diminished in 2021 and 2022 due to the dominance of strong southerly wind associated with the Amundsen Sea Low (ASL). This study emphasizes the potential of IS2 sea ice products to assess the role of atmospheric driving forces on thermodynamic and dynamic sea ice changes.

A Simple and Robust CryoSat‐2 Radar Freeboard Correction Method Dedicated to TFMRA50 for the Arctic Winter Snow Depth and Sea Ice Thickness Retrieval

AbstractCryoSat‐2 has been successful in observing sea ice thickness from space by providing ice freeboard information. The initial estimate of the ice freeboard, called radar freeboard, is obtained by analyzing the observed waveform using a retracker. A series of corrections are needed to convert the radar freeboard to the ice freeboard. Those are the physical effects (e.g., changes in wave propagation speed and the distribution of scattering at snow and ice surfaces, etc.) and the bias of the retracker; however, traditionally, only the wave speed correction has been applied due to lack of enough information to perform the complete correction. Here, an alternative correction method for the CryoSat‐2 radar freeboard derived using the Threshold First‐Maximum Retracker Algorithm with a 50% threshold (TFMRA50) is proposed. Snow depth was used as a predictor for the correction, similar to the traditional wave speed correction, but the coefficients were empirically determined by performing a direct comparison of the radar freeboard from CryoSat‐2 and the ice freeboard from airborne observations. Consequently, this new empirical correction treats the physical effects and the retracker bias as a whole, which have been difficult to separate in the retrieval process. In this paper, we demonstrate that the retrieval accuracy of snow and ice variables and the consistency of the two independent retrieval methods are improved when the new correction is applied. The result of this study emphasizes the importance of compatibility between the retracker and the freeboard correction method.

Entrapment in CO2 and H2O Ices: Impact of Ice Matrix Thickness

The volatile inventories of planets depend on the precise partitioning of different volatiles between the gas and solid phases across the planet birth disk. For the case of hyper-volatiles (e.g., CO, N2, and noble gases), the partitioning also depends on how efficiently they are trapped into less volatile ice matrices. The thicknesses of these ice matrices can range from a few molecular layers to macroscopic bodies, and in this study we explore how entrapment efficiency depends on the ice thickness between tens of nanometers and a few micrometers (∼50–3000 monolayers, ML). We carry out a series of temperature-programmed desorption (TPD) experiments on H2O:CO and CO2:CO mixtures with 5:1 and 15:1 matrix-to-CO mixing ratios. Entrapment efficiencies range from 41% to 64%, with higher entrapment efficiencies for the more dilute ices. Surprisingly, we find no significant difference in entrapment across the studied ice thicknesses for either H2O or CO2 ice matrices. Complementary TPD experiments with the additional hyper-volatiles N2 and Ar see a similar trend with ice thickness. We speculate that these results may be due to surface topography such as cracks that lead to hyper-volatile escape from deep ice layers. In either case, these experiments show that entrapment in microscopic ices is relatively insensitive to ice thickness (above ∼50 ML). In protoplanetary disks we therefore expect efficient entrapment in icy grains of a range of grain sizes.

EDPNet: A Transmission Line Ice-Thickness Recognition End-Side Network Based on Efficient Dynamic Perception

Ice-covered transmission lines seriously affect the stable operation of the power system. Deploying a recognition network for measuring the ice thickness on transmission lines within a deicing robot, and controlling the robot to perform resonant deicing, is an effective solution. In order to solve the problem that the existing recognition network is not suitable for an edge device, an ice-thickness recognition network for transmission lines based on efficient dynamic perception (EDPNet) is proposed. Firstly, a lightweight multidimensional recombination convolution (LMRC) is designed to split the ordinary convolution for lightweight design and extract feature information of different scales for reorganization. Then, a lightweight deep fusion module (LDFM) is designed, which combines the attention mechanism with different features to enhance the information interaction between the encoder and decoder. Then, a new dynamic loss function is adopted in the training process to guide the model to perform refined detection of ice-covered boundaries. Finally, we count the ice pixels and calculate the ice thickness. The model is deployed on an OrangePi5 Plus edge computing board. Compared with the baseline model, the maximum ice-thickness detection error is 4.2%, the model parameters are reduced by 86.1%, and the detection speed is increased by 74.6%. Experimental results show that EDPNet can efficiently complete the task of identifying ice-covered transmission lines and has certain engineering application value.

Formation and evolution of thermokarst landslides in the Qinghai-Tibet Plateau, China

Thermokarst landslide (TL) activity in the Qinghai-Tibet Plateau (QTP) is intensifying due to climate warming-induced permafrost degradation. However, the mechanisms driving landslide formation and evolution remain poorly understood. This study investigates the spatial distribution, annual frequency, and monthly dynamics of TLs along the Qinghai-Tibet engineering corridor (QTEC), in conjunction with in-situ temperature and rainfall observations, to elucidate the interplay between warming, permafrost degradation, and landslide activity. Through the analysis of high-resolution satellite imagery and field surveys, we identified 1298 landslides along the QTEC between 2016 and 2022, with an additional 386 landslides recorded in a typical landslide-prone sub-area. In 2016, 621 new active-layer detachments (ALDs) were identified, 1.3 times the total historical record. This surge aligned with unprecedented mean annual and August temperatures. The ALDs emerged primarily between late August and early September, coinciding with maximum thaw depth. From 2016 to 2022, 97.8 % of these ALDs evolved into retrogressive thaw slumps (RTSs), identified as active landslides. Landslides typically occur in alpine meadows at moderate altitudes and on gentle northward slopes. The thick ice layer near the permafrost table serves as the material basis for ALD occurrence. Abnormally high temperature significantly increased the active layer thickness (ALT), resulting in melting of the ice layer and formation of a thawed interlayer, which was the direct causing factor for ALD. By altering the local material, micro-topography, and thermal conditions, ALD activity significantly increases RTS susceptibility. Understanding the mechanisms of ALD formation and evolution into RTS provides a theoretical foundation for infrastructure development and disaster mitigation in extreme environments.

Field trial of a seismoacoustic method for ice cover parameters monitoring on the Franz Josef Land archipelago

Among the physical parameters of the freezing seas ice cover, ice thickness is of key importance, and its measurement is one of the most important tasks. The increased interest in the state of the sea ice cover as an indicator of global climatic changes, as well as the growth of comprehensive development of the Arctic shelf has caused intensive development of technical and methodological bases for ice observations. Despite the great variety of approaches to ice thickness estimation, all of them are not without weaknesses. Thus, most contact methods imply direct human presence, which significantly complicates the procedure, taking into account, among other factors, the rough weather conditions of the Arctic. Remote methods depend on weather conditions and cannot always provide high spatial resolution. In this connection, it is promising to use satellite observations coupled with the results of autonomous “ground” measurements, which can be seismoacoustic data containing information on the characteristics of elastic waves propagating in the ice-covered sea, is promising. The purpose of this work is to experimentally test a new passive method for monitoring ice cover parameters along long profiles based on the analysis of natural seismoacoustic fields. The article analyzes the data of a full-scale seismoacoustic experiment with a multichannel group of geophones placed on the floating ice of Alexandra Island in the Franz Josef Land archipelago within the framework of a complex expedition of the Russian Geographical Society. The demonstrates that it is in principle possible to use flexural-gravity waves propagating in the floating ice to estimate its characteristics, both in the active mode and by analyzing the ambient noise, is demonstrated. The results of ice parameter reconstruction obtained in a nondestructive manner using seismoacoustic waves and averaged over long profiles are compared with the data of direct contact measurements. This can be further used for monitoring seasonal and multiyear variability of sea ice thickness of freezing seas, including shelf zones.

Ice-marginal volcanic sequence in Iceland found on a nondescript gradual hillslope: An unexpected record of ice thickness late in deglaciation

Volcanism increases when glaciers melt because isostatic rebound during deglaciation decreases the pressure on the mantle, which enhances decompression melting. Anthropogenic climate change is now causing ice sheets and valley glaciers to melt around the world and this deglaciation could stimulate volcanic activity and associated hazards in Iceland, Antarctica, Alaska, and Patagonia. However, current model predictions for volcanic activity associated with anthropogenic deglaciation in Iceland are poorly constrained, in part due to uncertainties in past volcanic output over time compared to ice sheet arrangements. Further work specifically characterizing glaciovolcanic and ice-marginal volcanoes in Iceland is needed to reconstruct volcanic output during time periods with changing ice cover. Here, we describe a previously unrecognized ice-marginal volcanic sequence on a broad, gradual hillslope southeast of Langjökull and the Jarlhettur volcanic chain in Iceland's Western Volcanic Zone. Although previously mapped as interglacial lavas, canyons in this area revealed two southwest-dipping sequences of pillow-bearing tuff-breccias between pāhoehoe lava flows above modern lake Sandvatn. These pillow-bearing tuff-breccias and the quenched meter-scale cavities in coherent lava and cube-jointed facies show lavas came into contact with ice and pockets of trapped meltwater. However, clasts within the tuff-breccias include a mixture of pillow lavas and pāhoehoe fragments, requiring that the subaqueous tuff-breccia facies were derived from subaerial flows. In addition, we observed interfingering of subaerial and transitional subaqueous-subaerial pāhoehoe lava flows with the pillow-bearing tuff-breccias. We propose that during a deglaciation, subaerial lavas sourced upslope from near Skálpanes flowed downslope to the south and came into contact with thin ice north of the modern lake Sandvatn. We constrain the local ice at this time to be ∼30–50 m thick. Importantly, this finding demonstrates that ice-marginal deposits that can provide paleo-environmental constraints may be hidden in terrains without morphologically distinct glaciovolcanic edifices.

Quantifying the influence of snow over sea ice morphology on L-band passive microwave satellite observations in the Southern Ocean

Abstract. Antarctic snow on sea ice can contain slush, snow ice, and stratified layers, complicating satellite retrieval processes for snow depth, ice thickness, and sea ice concentration. The presence of moist and brine-wetted snow alters microwave snow emissions and modifies the energy and mass balance of sea ice. This study assesses the impact of brine-wetted snow and slush layers on L-band surface brightness temperatures (TBs) by synergizing a snow stratigraphy model (SNOWPACK) driven by atmospheric reanalysis data and the RAdiative transfer model Developed for Ice and Snow in the L-band (RADIS-L) v1.0 The updated RADIS-L v1.1 further introduces parameterizations for brine-wetted snow and slush layers over Antarctic sea ice. Our findings highlight the importance of including both brine-wetted snow and slush layers in order to accurately simulate L-band brightness temperatures, laying the groundwork for improved satellite retrievals of snow depth and ice thickness using satellite sensors such as Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP). However, biases in modelled and observed L-band brightness temperatures persist, which we attribute to small-scale sea ice heterogeneity and snow stratigraphy. Given the scarcity of comprehensive in situ snow and ice data in the Southern Ocean, ramping up observational initiatives is imperative to not only provide satellite validation datasets but also improve process-level understanding that can scale up to improving the precision of satellite snow and ice thickness retrievals.

Full-parametric and joint inversion of multi-mode surface wave data for identifying glacial ice thickness and freezing extent in subglacial sediments via Hunger Games Search Algorithm

Full-parametric and joint inversion of multi-mode surface wave data for identifying glacial ice thickness and freezing extent in subglacial sediments via Hunger Games Search Algorithm

Aufeis thickness and volume estimations from stereo satellite imagery and terrestrial photographs: Evidence from Central Ladakh, India

Aufeis – a phenomenon associated with permafrost and the cold arid region of Ladakh – serves as a critical water resource for local communities. In several tributaries of the Indus River, aufeis accumulation is enhanced in ice reservoirs (commonly known as “artificial glaciers”) to store winter baseflow for crop irrigation during the water-scarce period in spring. This study investigates aufeis thickness and volume across four study sites in the Trans-Himalaya of Central Ladakh: the ice reservoirs of Phuktse and Igoo and the catchments of Gya and Sasoma, where natural aufeis fields occur. Aufeis thickness and volume estimates were derived via differencing of digital elevation models calculated from very high-resolution stereo Pléiades imagery and terrestrial photographs. The study revealed ice thickness up to a maximum of 3 m and could further demonstrate the amplification effect of the walls of the ice reservoirs on aufeis accumulation in both, differenced elevation from Pléiades and terrestrial photographs. Aufeis volumes across the four study sites range from 34,106 ± 13,440 m3 in Phuktse up to 105,790 ± 28,511 m3 in Sasoma, indicating substantial amounts that need to be considered in future hydrological studies in the region. While the approach with terrestrial photographs is especially suitable for local studies, the results from very high-resolution stereo satellite data are promising for aufeis studies on large spatial scales. It represents an important contribution for the understanding of spatial aufeis patterns in the Trans-Himalaya and affected regions worldwide.

Method for Producing Columnar Ice in Laboratory and Its Application

This study presents the design of a small open-circuit wind tunnel for laboratory use and a method for preparing columnar ice. The ice formation process was analyzed in terms of temperature and ice thickness variations under varying environmental temperatures and wind speeds. Observations revealed that as wind speed increased, the grain size of the columnar ice decreased. Key findings include the following: (1) the selection and validation of two cubic arcs for the wind tunnel contraction section, achieving an acceleration ratio of 6.7–6.8 and stable wind speeds of 1–10 m/s; (2) real-time temperature monitoring indicated rapid cooling before freezing and slower cooling post-freezing, with lower ambient temperatures and higher wind speeds accelerating the icing process; (3) the −1/2 power of grain size was found to be positively correlated with wind speed; and (4) the method’s feasibility for studying mechanical properties of polar columnar ice was confirmed. This technique offers a controlled approach for producing columnar ice in the laboratory, facilitating comprehensive research on ice properties and providing a foundation for future studies on the mechanical behavior of ice under windy polar conditions.

Calculation of Ice Thickness of 1% Exceedance Probability and Assessment of its Long-Term Variability at the Mouth of the Oka River

Calculation of Ice Thickness of 1% Exceedance Probability and Assessment of its Long-Term Variability at the Mouth of the Oka River

Winter arctic sea ice volume decline: uncertainties reduced using passive microwave-based sea ice thickness

Arctic sea ice volume (SIV) is a key climate indicator and memory source in sea ice predictions and projections, yet suffering from large observational and model uncertainty. Here, we test whether passive microwave (PMW) data constrain the long-term evolution of Arctic SIV, as recently hypothesized. We find many commonalities in Arctic SIV changes from a PMW sea ice thickness (SIT) 1992-2020 time series reconstructed with a neural network algorithm trained on lidar altimetry, and the reference PIOMAS reanalysis: relatively low differences in SIV mean (4615 km3, 37%), SIV trends (46 km3, 17%), and phased variability (r2=0.55). Key to reduced differences is the consistent evolution of many SIV contributors: seasonal and perennial ice coverage, their SIT contrast, whereas perennial SIT provides the largest remaining uncertainty source. We argue that PMW includes useful SIT information, reducing SIV uncertainty. We foresee progress from sea ice reanalyses combining dynamical models and data assimilation of PMW SIT estimates, in addition to the already assimilated PWM sea ice concentration.

Influence of the Atlantic and Pacific Multidecadal Variability on Arctic Sea Ice in Pacemaker Simulations during 1920–2013

Abstract The Atlantic multidecadal variability (AMV) and Pacific multidecadal variability (PMV) can influence Arctic sea ice and modulate its trend, but to what extent the AMV and PMV can affect Arctic sea ice and which processes are dominant are not well understood. Here, we analyze the Community Earth System Model, version 1, idealized and time-varying pacemaker ensemble simulations to investigate these issues. These experiments show that the sea ice concentration varies mainly over the marginal Arctic Ocean, while the sea ice thickness variations occur over the entire Arctic Ocean. The internal components of AMV and PMV can enhance or weaken the decadal sea ice loss rates over the marginal Arctic Ocean by more than 50%. The AMV- or PMV-induced anomalous atmospheric energy transport and downward longwave radiation related to low clouds (thermodynamical processes) and sea ice motion (dynamical processes) contribute to the Arctic surface air temperature and sea ice concentration and thickness changes. Anomalous oceanic heat flux is mainly a response to rather than a cause of sea ice variations. The dynamic processes contribute to the winter Arctic sea ice variations as much as the thermodynamic processes, but they contribute less (more) to the summer Arctic sea ice variability than the thermodynamic processes over the marginal Arctic Ocean (parts of the central Arctic Ocean). Sea ice loss enhances air–sea heat fluxes, which cause oceanic heat convergence and warm near-surface air and the lower troposphere, which in turn melt more sea ice.

The connection between sea ice dynamics enhancement and atmospheric circulation in the East Siberian Sea during spring

The change in sea ice thickness can be divided into dynamical and thermodynamical effect. In the last four decades, the drastic changes in multi-year ice thickness in East Siberian Sea during spring have made the sea ice more susceptible to dynamical effect. On one hand, the dynamical effect on sea ice in this region is the strongest, surpassing that of other Arctic marginal seas, and has been continuously strengthening since 1996. On the other hand, this region’s dynamical effect varies with the Central Arctic Index (CAI). An increase in CAI extreme positive phase years leads to more frequent occurrences of cyclonic circulation anomalies, favoring the development of dynamical effect in spring East Siberian Sea. Furthermore, the influence of 10 m wind field on dynamical effect has shifted from being dominated by the northward component (v) to a combined effect of both northward and eastward components (v and u). This change is a result of alteration in the response pattern of wind field to CAI across the interdecadal periods.

Research on icing model and calculation methods

Aircraft icing seriously threatens flight safety. This paper describes a modified shallow-water icing thermodynamic model that is applicable to unstructured grids and considers the effects of changes in water physical parameters and sublimation on icing. An icing calculation method enables the automatic determination of whether freezing occurs and the type of icing. An autonomous icing calculation program is developed to extract the geometric coordinates and mesh information of the icing surface and combine this information with the multiphase flow field of air-supercooled water droplets. The icing equations in the Godnov format are discretized to produce a large-scale sparse matrix that is solved iteratively using the biconjugate gradient stabilized method. The output includes parameter distributions for the ice thickness, water film thickness, and equilibrium temperature. Taking the National Advisory Committee for Aeronautic 0012 airfoil as the study object, the results of icing simulations are compared with experimental data from an ice wind tunnel and the ice shapes calculated by the FENSAP-ICE software. The results are found to be in good agreement with the experimental data, and the ice height errors at stagnation points are less than 15%. In the case of rime ice, single-horned ice shapes are simulated for both conventional and large supercooled droplets. For mixed ice and glaze ice, double-horned ice shapes are simulated for both droplet conditions. FENSAP-ICE fails to simulate the ice horns and produces large errors in ice thickness and ice range.

Temporal Behavior of Glacially Induced Stresses and Strains at Potential Sites for Long‐Term Storage of Used Nuclear Fuel in Canada

AbstractContinental‐scale glaciations cause deformation, geopotential, rotation and stress changes of the Earth. Subsurface stress changes have implications to future activities such as carbon capture and storage, enhanced oil recovery and deep geological disposal of nuclear waste. We model glacially induced stresses, strain changes and deformation for North America, with emphasis on the two potential sites for long‐term storage of used nuclear fuel in Canada (Saugeen Ojibway Nation (SON)‐South Bruce area in southwestern Ontario and Wabigoon Lake Ojibway Nation (WLON)‐Ignace area in northwestern Ontario). We apply a revised, high‐resolution ice history of the past glacial cycle from the University of Toronto Glacial Systems Model, assumed to be representative for future glacial cycles, together with a set of seven different one‐ and three‐dimensional earth structures. We find that glacially induced stresses and strains can vary strongly throughout a glacial cycle, whereas especially the horizontal components can change from tensional to compressive in nature. Such changes can happen within a few 1,000 years, caused by drastic and rapid ice thickness increase or decrease above the potential site. Despite SON‐South Bruce being located further away from the ice sheet center than WLON‐Ignace and temporarily in the forebulge of the developing ice sheet during glaciation, stresses and strains are very similar in magnitude and range at both sites. We also see the potential that the glacially induced stresses can alter the direction of the pre‐existing maximum horizontal stress at SON‐South Bruce. These results will be incorporated in the site safety and site selection process.

Chlorophyll Production in the Amundsen Sea Boosts Heat Flux to Atmosphere and Weakens Heat Flux to Ice Shelves

AbstractThe Amundsen Sea in West Antarctica features rapidly thinning ice shelves, large polynyas, and sizable spring phytoplankton blooms. Although considerable effort has gone into characterizing heat fluxes between the Amundsen Sea, its associated ice shelves, and the overlying atmosphere, the effect of the phytoplankton blooms on the distribution of heat remains poorly understood. In this modeling study, we implement a feedback from biogeochemistry onto physics into MITgcm‐BLING and use it to show that high levels of chlorophyll—concentrated in the Amundsen Sea Polynya and the Pine Island Polynya—have the potential to increase springtime surface warming in polynyas by steepening the attenuation profile of solar radiation with depth. The chlorophyll‐associated warm anomaly (on average between +0.2C and +0.3C) at the surface is quickly dissipated to the atmosphere, by increases in longwave, latent and sensible heat loss from open water areas. Outside of the coastal polynyas, the summertime warm anomaly leads to an average sea ice thinning of 1.7 cm across the region, and stimulates up to 20% additional seasonal melting near the fronts of ice shelves. The accompanying cold anomaly, caused by shading of deeper waters, persists year‐round and affects a decrease in the volume of Circumpolar Deep Water on the continental shelf. This cooling ultimately leads to an average sea ice thickening of 3.5 cm and, together with associated changes to circulation, reduces basal melting of Amundsen Sea ice shelves by approximately 7% relative to the model scenario with no phytoplankton bloom.