Sea Ice Research Articles

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.

Relative contribution of the ocean-air heat exchange and advective heat transport to the increase of the Barents Sea water temperature in the early 21st century

The annual water temperature in the major water masses of the Barents Sea (BS) has significantly increased since the early 2000s. Advective heat transport from the neighboring water areas and heat exchange through the sea surface are the major factors, which shape the hydrological conditions in the BS. The paper estimates the contributions of heat exchange at the sea-atmosphere boundary and advective heat transport to changes in the average water temperature of the BS for the entire sea area. The average annual heat balance of the BS is calculated using atmospheric and oceanic reanalysis data. The change in the average temperature of the BS water is estimated taking into account the heat consumption for ice melting. The average surface heat balance from 1993 to 2018 was negative throughout the entire sea area: –70…–100 W/m2 in the south and –10…–20 W/m2 in the north. The advective heat supply was calculated for 9 straits with neighboring water areas. The determining source of advective heat is the influx of Atlantic waters from the Norwegian Sea between Cape Nordkapp and Bear Island. An average of 40.8 TW of advective heat is supplied through this margin. The calculations showed the predominance of annual heat influx due to advection over heat loss from the sea surface. This excess heat influx resulted in an estimated increase in the water temperature of the BS from 1993 to 2018 at a rate of 0.28 °C per year (taking into account the heat consumption for ice melting). In conclusion, it can be argued that the analysis has validated the hypothesis proposed in the article about compensation of heat losses from the surface of the BS by advective heat flow. The hypothesis is quantitatively confirmed by calculations on a simple box model (with an accuracy of up to an order of magnitude) based on atmospheric and oceanic reanalysis data. The ERA5 and GLORYS12V1 reanalysis data reliably describe the basic patterns of observed variability of ocean, sea ice and atmospheric parameters in the Barents Sea.

Warming events and their causes at the Bering Sea section B in summer of 1999–2019

Based on hydrological data for the Bering Sea from 1999 to 2019 during the Chinese National Arctic Research Expeditions (CHINAREs), along with oceanographic and meteorological data from the National Centers for Environmental Information (NCEI) and the European Centre for Medium-Range Weather Forecasts (ECMWF), we investigated trends in summer at the Bering Sea section B over the past 20 years and their underlying causes. The results revealed a pronounced rise in temperature at the Bering Sea section B in 2003, 2014, 2016, 2018, and 2019, suggesting a shift towards a new phase of warmth starting in 2014. In warm years (2003, 2014, 2016, 2018, and 2019), the mean maximum temperature in the upper layer in the basin north of 57°N and on the continental shelf were roughly 2.7 °C, 2.6 °C higher than that in cold years (1999, 2008, 2010, and 2012). Additionally, the mean minimum temperature of the cold water in the middle layer in the basin was roughly 0.7 °C higher in the warm years compared to the cold years. This difference between the warm and cold years was evident in the heat content: the average heat content in the upper to middle layers in the basin north of 57°N and on the continental shelf were about 1.2 and 0.9 GJ/m2 higher, respectively, in the warm years than the cold years. However, there was a modest warming trend in the basin south of 57°N, and the temperature in the mid-layer on the continental shelf showed year-to-year stability, in contrast to the significant warming observed elsewhere. Further study suggested that temperature fluctuation in the basin was closely tied to seasonal 2-m air temperature. In addition, the southern basin, divided at 57°N, experienced larger positive or negative anomalous wind stress curl compared to the northern basin. This discrepancy resulted in differing causes of temperature variation between the two regions. Temperature change in the upper layer on the continental shelf was attributed to an increase in cumulative net heat flux resulting from reduced sea ice, while the temperature in the middle layer on the continental shelf was limited by the enhanced density gradient between the upper and middle layer caused by melt ice, preventing significant warming. Overall, the thermal state for the Bering Sea section B was influenced by a combination of factors, including net heat flux, sea ice, and atmospheric conditions, contributing to the year-to-year variation characterized by alternating cold and warm regimes.

Dense Water Production in Storfjorden, Svalbard, From a 1‐Year Time Series of Observations and a Simple Model: Are Polynyas in a Warming Arctic Exporting Heat to the Deep Ocean?

AbstractThe formation of dense Brine‐enriched Shelf Water (BSW) in Storfjorden is analyzed during Winter 2016–2017 from mooring observations, a polynya model nudged to satellite observations, and an original BSW production model. The ice season was two months shorter than average, yet 44.2 of sea ice were formed, in line with estimates for the period preceding the atlantification of the Barents Sea in the mid‐2000s: A thinner, more fragile ice may favor polynya openings and frazil ice production. A saline specimen of BSW was produced in large volumes, corresponding to an annual mean transport of 0.042 Sv, larger than previous estimates. The important production is due to the preconditioning of the polynya with a more saline source water, exceeding the pre‐2005 values by 0.37. The BSW overflow was observed on the West Spitsbergen shelf slope from hydrographic sections down to 750 m, thus entering the Norwegian Sea Deep Water layer. Its core temperature was about C warmer than the pre‐2005 values owing to the entrainment of a warmer water in Storfjordrenna, suggesting that a part of the excess surface heat of the Barents Sea could be exported into the deep ocean. Overall our results suggest that dense water formation in the Storfjorden polynya may not, at least for now, be hampered by the atlantification of the Barents Sea, and perhaps even temporarily favored by the more saline source water. Anomalous atmospheric warming during the Winter‐Spring may however disrupt the production, as was observed 1 year before.

Improved snow depth retrieval over Arctic sea ice from FY-3 satellites using machine learning and its future scenario projection

Snowpack atop sea ice plays a vital role in the Arctic heat and energy balance by reducing the energy transfer between the atmosphere and the ocean. This study presents an improved snow depth retrieval model over Arctic multiyear ice (MYI) from the Microwave Radiometer Imager (MWRI) data of the Fengyun-3 (FY-3) series satellites using the optimal machine learning approach identified. Validation against the Operation IceBridge (OIB) data revealed our snow depth estimates had a root mean square error (RMSE) of 7.88 cm and 6.83 cm for snow depth estimates over first-year ice (FYI) and MYI, respectively. FY-3 derived snow depth estimates from 2010 to 2014 were used to evaluate snow depth outputs from the Coupled Model Intercomparison Project Phase 6 (CMIP6). Based on evaluation results, seven CMIP6 models were selected for snow depth predictions, namely ACCESS-ESM1-5, AWI-CM-1-1-MR, CAMS-CSM1-0, CanESM5, FGOALS-f3-L, MRI-ESM2-0, and NorESM2-MM. Cold season snow depths from 2015 to 2100 were generated under four shared socioeconomic pathway-representative concentration pathway (SSP-RCP) scenarios (SSP126, SSP245, SSP370, and SSP585). The results indicated that short-term increases in radiation forcing have minimal impact on snow depth variability; however, as time goes by and radiative forcing accumulates, the rate of snow depth decrease accelerates. Relative to the average snow depth in the baseline period (2010–2014), the percentage change in snow depth during the late 21st century (2081–2100) under the SSP585 scenario (−63.53 %) is more than twice the percentage change observed under the SSP126 scenario (−29.63 %).These findings highlight the potential for substantial snow loss over Arctic sea ice under high greenhouse gas emission scenarios.

Links between atmospheric aerosols and sea state in the Arctic Ocean

Sea spray emission is the largest mass flux of aerosols to the atmosphere with important impact on atmospheric radiative transfer. However, large uncertainties still exit in constraining this mass flux and its climate forcing, in particular in the Arctic, where sea ice and relatively low wind speed in summer constitute a significantly different regime compared to the global ocean. Sea state conditions and marine boundary layer stability are also critical variables, but their contribution is often overlooked. Here we present concurrent observations of sea state using a novel stereo camera system, of sea spray through coarse mode aerosols, and of meteorological variables to determine boundary layer stability in the Barents and Kara Seas during the 2021 Arctic Century Expedition. Our findings reveal that aerosol concentrations were highest over open waters, closely correlating with wave height, followed by wind speed, wave steepness, and wave age. Notably, these correlations were stronger under unstable marine boundary layer conditions, reflecting immediate sea spray generation. By analysing various combinations of sea and atmospheric variables, we identified the wave height Reynolds number as the most effective indicator of atmospheric sea spray concentration, explaining 57% of its variability in unstable conditions. Our study underscores the need to consider sea state, wind, and boundary layer conditions together to accurately estimate atmospheric sea spray concentrations in the Arctic.

High‐Resolution Fracture Dynamics Simulation of Pack‐Ice and Drift‐Ice Formation During Sea Ice Break up Events Using the HiDEM2.0 Code

High‐Resolution Fracture Dynamics Simulation of Pack‐Ice and Drift‐Ice Formation During Sea Ice Break up Events Using the HiDEM2.0 Code

Arctic Amplification of marine heatwaves under global warming

Marine heatwaves (MHWs) and total heat exposures (THEs), extreme warming events occurring across the global oceans, seriously threaten marine ecosystems and coastal communities as the climate warms. However, future changes in MHWs and THEs in the Arctic Ocean, where unique marine ecosystems are present, are still unclear. Here, based on the latest CMIP6 climate simulations, we find that both MHWs and THEs in the Arctic Ocean are anticipated to intensify in a warming climate, mainly due to Arctic sea ice decline and long-term warming trend, respectively. Particularly striking is the projected rise in MHW mean intensity during the 21st century in the Arctic Ocean, surpassing the global average by more than sevenfold under the CMIP6 SSP585 scenario. This phenomenon, coined the ‘Arctic MHW Amplification’, underscores an impending and disproportionately elevated threat to the Arctic marine life, necessitating targeted conservation and adaptive strategies.

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.

Intersensor Calibration of Spaceborne Passive Microwave Radiometers and Algorithm Tuning for Long-Term Sea Ice Trend Analysis Based on AMSR-E Observations

Sea ice monitoring is key to analyzing the Earth’s climate system. Long-term sea ice extent (SIE) has been continuously monitored using various spaceborne passive microwave radiometers (PMRs) since November 1978. As the lifetime of a satellite is usually approximately 5 years, bias caused by differences in PMRs should be eliminated to obtain objective SIE trends. Most sea ice products have been analyzed for long-term trends with a bias adjustment based on the coarse resolution special sensor microwave imager (SSM/I) in operation for the longest period. However, since 2002, Japanese microwave radiometers of the Advanced Microwave Scanning Radiometer (AMSR) series, which have the highest spatial resolution in PMR, have been available. In this study, we developed standardization techniques for processing SIE including calibration of the brightness temperature (TB), tuning the sea ice concentration (SIC) algorithm, and adjusting the SIC threshold to retrieve a consistent SIE trend based on the AMSR for the Earth Observing System (AMSR-E, one of the AMSR that operated from May 2002 to October 2011). Analysis results showed that the root-mean-square error between AMSR-E SICs and those of moderate resolution imaging spectroradiometer (MODIS) was 15%. In this study, SIE was defined as the sum of the areas where the AMSR-E SIC was >15%. When retrieving SIE, we adjusted the SIC threshold for each PMR to be consistent with the SIE calculated based on the 15% SIC threshold for AMSR-E. We then calculated a time-series of the SIE trends over approximately 45 years using the adjusted SIE data. Therefore, we revealed the dramatic decrease in global sea ice extent since 1978. This technique enables retrieval of more accurate long-term sea ice trends for more than half a century in the future.

Regulation of Ocean Surface Currents and Seasonal Sea Ice Variations on the Occurrence and Transport of Organophosphate Esters in the Central Arctic Ocean.

Organophosphate esters (OPEs) have been observed in the remote Arctic Ocean, yet the influence of hydrodynamics and seasonal sea ice variations on the occurrence and transport of waterborne OPEs remains unclear. This study comprehensively examines OPEs in surface seawater of the central Arctic Ocean during the summer of 2020, integrating surface ocean current and sea ice concentration data. The results confirm significant spatiotemporal variations of the OPEs, with the total concentration of seven major OPEs averaging 780 ± 970 pg/L. Chlorinated OPEs, particularly tris(1-chloro-2-propyl) phosphate (TCPP), were dominant. The significant impact of hydrodynamics on the OPE transport is demonstrated by higher OPE concentrations in regions with strong surface currents, especially at the edge of the Beaufort Gyre and the confluence of the Beaufort Gyre and the Transpolar Drift. Furthermore, OPE levels were generally higher in drifting-ice-covered regions compared to ice-free regions, attributed to the volatilization of dissolved OPEs formerly trapped below the sea ice or newly released from melting snow and sea ice. Notably, TCPP decreased by only 19% in the ice-free area, while the more volatile triphenyl phosphate decreased by 63% compared with the partial ice region.

Sea ice dynamics influence movement patterns of adult female polar bears in southern Hudson Bay

Sea ice dynamics influence movement patterns of adult female polar bears in southern Hudson Bay

How many parameters are needed to represent polar sea ice surface patterns and heterogeneity?

Abstract. Sea ice surface patterns encode more information than can be represented solely by the ice fraction. The aim of this paper is thus to establish the importance of using a broader set of surface characterization metrics and to identify a minimal set of such metrics that may be useful for representing sea ice in Earth system models. Large-eddy simulations of the atmospheric boundary layer over various idealized sea ice patterns, with equivalent ice fractions and average floe areas, demonstrate that the spatial organization of ice and water can play a crucial role in determining boundary layer structures. Thus, various methods used to quantify heterogeneity in categorical lattice-based spatial data, such as those used in landscape ecology and Geographic Information System (GIS) studies, are employed here on a set of recently declassified high-resolution sea ice surface images. It is found that, in conjunction with ice fraction, patch density (representing the fragmentation of the surface), the splitting index (representing variability in patch size), and the perimeter–area fractal dimension (representing the tortuosity of the interface) are all required to describe the two-dimensional pattern exhibited by a sea ice surface. For surfaces with anisotropic patterns, the orientation of the surface relative to the mean wind is also needed. Finally, scaling laws are derived for these relevant landscape metrics, allowing for their estimation using aggregated spatial sea ice surface data at any resolution. The methods used in and the results gained from this study represent a first step toward developing further methods for quantifying variability in polar sea ice surfaces and for parameterizing mixed ice–water surfaces in coarse geophysical models.

Enhanced ocean heat storage efficiency during the last deglaciation

Proxy reconstructions suggest that increasing global mean sea surface temperature (GMSST) during the last deglaciation was accompanied by a comparable or greater increase in global mean ocean temperature (GMOT), corresponding to a large heat storage efficiency (HSE; ∆GMOT/∆GMSST). An increased GMOT is commonly attributed to surface warming at sites of deepwater formation, but winter sea ice covered much of these source areas during the last deglaciation, which would imply an HSE much less than 1. Here, we use climate model simulations and proxy-based reconstructions of ocean temperature changes to show that an increased deglacial HSE is achieved by warming of intermediate-depth waters forced by mid-latitude surface warming in response to greenhouse gas and ice sheet forcing as well as by reduced Atlantic meridional overturning circulation associated with meltwater forcing. These results, which highlight the role of surface warming and oceanic circulation changes, have implications for our understanding of long-term ocean heat storage change.

Development of a total variation diminishing (TVD) sea ice transport scheme and its application in an ocean (SCHISM v5.11) and sea ice (Icepack v1.3.4) coupled model on unstructured grids

Abstract. As the demand for increased resolution and complexity in unstructured sea ice models is growing, higher demands are also placed on the sea ice transport scheme. In this study, we couple the Semi-implicit Cross-scale Hydro-science Integrated System Model (SCHISM, v5.11) with Icepack (v1.3.4), the column physics package of the Los Alamos sea ice model (CICE); a key step is to implement a total variation diminishing (TVD) transport scheme for the multi-class sea ice module in the coupled model. Compared with the second-order upwind scheme and the finite-element flux-corrected transport (FEM-FCT) scheme, the TVD transport scheme is overall superior when evaluated based on conservation, accuracy, efficiency (even with very high resolution), and strict monotonicity. Although it is slightly weaker than FEM-FCT in terms of accuracy alone, the TVD scheme still outperforms the other two schemes in comprehensive performance. The new coupled model outperforms the existing single-class ice model of SCHISM in the case of Lake Superior. For the Arctic Ocean case, it successfully reproduces the long-term changes in the sea ice extent, sea ice boundary, concentration observations from satellites, and thickness from in situ measurement.

Dye tracing of upward brine migration in snow

Abstract Salt is often present in the snow overlying seasonal sea ice, and has profound thermodynamic and electromagnetic effects. However, its provenance and behaviour within the snow remain uncertain. We describe two investigations tracing upward brine movement in snow: one conducted in the laboratory and one in the field. The laboratory experiments involved the addition of dyed brine to the base of terrestrial snow samples, with subsequent wicking being measured. Our field experiment involved dye being added directly (without brine) to bare sea-ice and lake ice surfaces, with snow then accumulating on top over several days. On the sea ice, the dye migrated upwards into the snow by up to 5 cm as the snow's basal layer became more salty, whereas no migration occurred in our control experiment over non-saline lake ice. This occurred in relatively dry snowpacks where brine took up $< 6\%$ of the snow's calculated pore volume, suggesting pore saturation is not required for upward salt transport. Our results highlight the potential role of microstructural parameters beyond those currently retrievable with penetrometry, and the potential value of longitudinal, process-based field studies of young snowpacks.

Split-explicit external mode solver in the finite volume sea ice–ocean model FESOM2

Abstract. A novel split-explicit (SE) external mode solver for the Finite volumE Sea ice–Ocean Model (FESOM2) and its sub-versions is presented. It is compared with the semi-implicit (SI) solver currently used in FESOM2. The split-explicit solver for the external mode utilizes a dissipative asynchronous (forward–backward) time-stepping scheme that does not require additional filtering during the barotropic sub-cycling. Its implementation with arbitrary Lagrangian–Eulerian vertical coordinates like Z star (z*) and Z tilde (z̃) is explored. The comparisons are performed through multiple test cases involving idealized and realistic global simulations. The SE solver demonstrates lower phase errors and dissipation, while maintaining a simulated mean ocean state very similar to the SI solver. The SE solver is also shown to possess better run-time performance and parallel scalability across all workloads tested.

The Southern Ocean marine ice record of the early historical, circum-Antarctic voyages of Cook and Bellingshausen

Abstract. The circumnavigations of Cook (second voyage, 1772–1775) and Bellingshausen (1819–1821) were attempts to find any great southern land mass poleward of ∼ 50° S and consequently involved sailing for three or two summers, respectively, in polar latitudes around Antarctica. Extensive sea ice eventually blocked each voyage's southern probes, although Bellingshausen, unknowingly at the time, saw the Antarctic continent. However, these attempts meant sea ice and iceberg records from the early historical period were collected nearly simultaneously from around much of Antarctica. Here, these records are extracted from journals, analysed, and compared to each other and the modern satellite record of both forms of marine ice. They generally show an early historical period with a more northerly record of both forms of marine ice than normal for today, but to a geographically varying degree. However, the early historical period in the Pacific sector of the Southern Ocean saw marine ice generally within the range of modern observations for the same time of year, but the Weddell Sea and Indian Ocean marine ice, particularly on Cook's voyage, then extended several degrees further north than in today's extreme ice years.

Calibration of a hybrid sea ice model during an expedition to the Arctic

AbstractWe describe the acquisition and evaluation of data during the ArcWatch II expedition into the Arctic Ocean on board of the icebreaker Polarstern. The data is used to calibrate a hybrid sea ice model that focuses on modeling the marginal ice zone, which is the transition zone between pack ice and the open ocean. We describe the experimental setup used to track the motion of individual ice floes as well as the measurement of the wind and the ocean current. In the literature, the motion of an ice floe in the free‐drift regime is often estimated by Nansen's rule of thumb, which assumes the sea ice drifts at of the wind speed, deflected up to in clockwise direction. Analyzing the first datasets, we noticed that observed small ice floes in free drift often have a larger relative angle to the wind direction. Based on numerical simulation we quantify that these larger angles are related to the geometry of larger sea ice floes in the vicinity of the observed ice floe. These large icefloes influence the ocean current, which in turn, impacts the motion of the small ice floes in the surrounding.

A new high-resolution Coastal Ice-Ocean Prediction System for the East Coast of Canada

The Coastal Ice Ocean Prediction System for the East Coast of Canada (CIOPS-E) was developed and implemented operationally at Environment and Climate Change Canada (ECCC) to support a variety of critical marine applications. These include support for ice services, search and rescue, environmental emergency response and maritime safety. CIOPS-E uses a 1/36° horizontal grid (~ 2 km) to simulate sea ice and ocean conditions over the northwest Atlantic Ocean and the Gulf of St. Lawrence (GSL). Forcing at lateral open boundaries is taken from ECCC’s data assimilative Regional Ice-Ocean Prediction System (RIOPS). A spectral nudging method is applied offshore to keep mesoscale features consistent with RIOPS. Over the continental shelf and GSL, the CIOPS-E solution is free to evolve according to the model dynamics. Overall, CIOPS-E significantly improves the representation of tidal and sub-tidal water levels compared to ECCC’s lower resolution systems: RIOPS (~ 6 km) and the Regional Marine Prediction System – GSL (RMPS-GSL, 5 km). Improvements in the GSL are due to the higher resolution and a better representation of bathymetry, boundary forcing and dynamics in the upper St. Lawrence Estuary. Sea surface temperatures show persistent summertime cold bias, larger in CIOPS-E than in RIOPS, as the latter is constrained by observations. The seasonal cycle of sea ice extent and volume, unconstrained in CIOPS-E, compares well with observational estimates, RIOPS and RMPS-GSL. A greater number of fine-scale features are found in CIOPS-E with narrow leads and more intense ice convergence zones, compared to both RIOPS and RMPS-GSL.