Ice Thickness Research Articles

Predicting the current fishable habitat distribution of Antarctic toothfish (Dissostichus mawsoni) and its shift in the future under climate change in the Southern Ocean.

Global warming continues to exert unprecedented impacts on marine habitats. Species distribution models (SDMs) are proven powerful in predicting habitat distribution for marine demersal species under climate change impacts. The Antarctic toothfish, Dissostichus mawsoni (Norman 1937), an ecologically and commercially significant species, is endemic to the Southern Ocean. Utilizing occurrence records and environmental data, we developed an ensemble model that integrates various modelling techniques. This model characterizes species-environment relationships and predicts current and future fishable habitats of D. mawsoni under four climate change scenarios. Ice thickness, depth and mean water temperature were the top three important factors in affecting the distribution of D. mawsoni. The ensemble prediction suggests an overall expansion of fishable habitats, potentially due to the limited occurrence records from fishery-dependent surveys. Future projections indicate varying degrees of fishable habitat loss in large areas of the Amery Ice Shelf's eastern and western portions. Suitable fishable habitats, including the spawning grounds in the seamounts around the northern Ross Sea and the coastal waters of the Bellingshausen Sea and Amundsen Sea, were persistent under present and future environmental conditions, highlighting the importance to protect these climate refugia from anthropogenic disturbance. Though data deficiency existed in this study, our predictions can provide valuable information for designing climate-adaptive development and conservation strategies in maintaining the sustainability of this species.

Retrieval of sea ice thickness from FY-3E data using Random Forest method

In this study, we employ a Random Forest approach to estimate sea ice thickness (SIT) using Fengyun-3E (FY-3E) and Soil Moisture Ocean Salinity (SMOS) data. This method relies on four input parameters: incidence angle (θ), reflectivity (Γ), sea ice salinity (S), and sea ice temperature (T). In addition, FY-3E can receive both Global Positioning System (GPS) and Beidou Navigation Satellite System (BDS) reflected signals. Evaluation for the Arctic region based on data spanning from October 2022 to April 2023 reveals that the proposed models trained on GPS and BDS signals from FY-3E achieve high consistency and low error. Take GPS signals as an example, coefficients of determination are 0.97 and 0.91 and mean absolute errors are 0.019 m and 0.032 m for the training and test sets, respectively. In general, SIT inversion based on GPS signals slightly exhibits a higher accuracy than that based on BDS signals, but both approaches display high performances. The areas with the highest accuracy of SIT estimation based on GPS and BDS signals are the Shelikhov Bay and the Okhotsk Sea, followed by the Bering Sea and the Bering Strait. We conclude that machine learning and data fusion are effective for SIT estimation.

SWAT‐GL: A new glacier routine for the hydrological model SWAT

AbstractThe hydrological model Soil Water Assessment Tool (SWAT) is widely used in water resources management worldwide. It is also used to simulate catchment hydrology in high‐mountainous regions where glaciers play an important role. However, SWAT considers glaciers in a simplistic way. Although some efforts were done to overcome this limitation, there is no official version available that considers glaciers adequately. This strongly impairs its applicability in glacierized catchments. In this technical note, we propose a novel version of the traditional SWAT, called SWAT‐GL, which introduces (1) a mass balance module and (2) a glacier evolution routine to represent dynamic glacier changes. Mass balance calculations are based on a conceptual degree‐day approach, similar to the snow routine implemented in SWAT. Glacier evolution is realized using the delta‐h (∆h) parameterization, which requires a minimum of data and is thus suitable in data‐scarce regions. The approach allows users to simulate spatially distributed glacier changes. Annual mass balance changes are translated to distributed ice thickness changes depending on the glacier elevation. We demonstrate how SWAT‐GL is technically integrated into SWAT and how glaciers are merged with the existing spatial units. Model code and test data are freely accessible to promote further model development efforts and a wide application. Ultimately, SWAT‐GL aims to make SWAT easily applicable in glacierized catchments without the need of additional tools.

Mapping Ice Buried by the 1875 and 1961 Tephra of Askja Volcano, Northern Iceland Using Ground‐Penetrating Radar: Implications for Askja Caldera as a Geophysical Testbed for In Situ Resource Utilization

AbstractEruptions of the Askja Volcano in Northern Iceland in 1875 and 1961 blanketed the caldera with rhyolitic and basaltic tephra deposits, respectively, which preserved layers of seasonal snowpack as massive ice. Askja serves as an operational and geophysical analog to test ground‐penetrating radar field and analysis techniques for in situ resource utilization objectives relevant to the martian and lunar environments. We conducted ground‐penetrating radar surveys at center frequencies of 200, 400, and 900 MHz to map the thickness and extent of tephra deposits and underlying massive ice at three caldera sites. We identified up to 1 m of tephra preserving up to 4.4 m of massive ice. We measured the real dielectric permittivity of the overlying tephra and the total attenuation at each frequency of the tephra and ice. A key objective of our investigation was to determine if attenuation (or loss) could be used as an additional diagnostic signature of massive ice preserved at depth when compared to ice‐free stratigraphy. Loss rates of the ice‐rich subsurface decrease with increasing ice thickness relative to the overburden, which may constitute a possible signature. Attenuation also increased with increasing frequency. The tephra, ice, and other volcanic deposits at each of our three caldera sites and the ice‐free, pumice‐mantled 1961 Vikrahraun lava flow exhibited consistently low loss rates at all frequencies. This result highlights the ambiguity associated with identifying the unique signature of ice within low‐loss stratigraphies, a possible challenge for its identification in the martian or lunar subsurface using radar.

Glaciological studies in Mexico, 60 years of academic work: A summary

Glaciers have played a very important role in controlling the climate during most of the geologic history of our planet Earth. Of course, the glaciers have always been at higher latitudes (north and south), and some on high-altitude mountains. Current glaciers in Mexico are those inherited from the Last Glacial Maximum, (26000 - 19000 years before the present), increasing in size during the period of the Little Ice Age (1300 to 1850 common era), and they are unique in several ways; they are located at 19° north latitude, and they received snow precipitation from both the Pacific Ocean and the Gulf of Mexico (Atlantic Ocean). Research on glacial chronology, physical glaciology, and glacial geochemistry at Universidad Nacional Autónoma de México has provided valuable information on climate and environmental changes at different time scales, from millennial to decadal, and even annual. The first part of this work deals with the reconstruction of the glacial history in Mexico and establishing a glacial chronology from the Last Glacial Maximum to the Little Ice Age. The second part of this work focuses on monitoring recent changes (the last 60 years, 1960s to 2023) of the glacier extent on Iztaccíhuatl, Popocatépetl, and Citlaltépetl (the three highest mountains in Mexico), as well as having an updated inventory of all the glaciers at those mountains. Changes in glacier extent and thickness of ice are directly related to the increase in air temperature, variation in precipitation patterns, and glacier dynamics on some of the last glaciers of the northern tropics. The third part or this work focuses on a compilation of geochemical data from 17 years (from 2006 to 2013) of sampling ice (shallow ice cores) and snow at Iztaccíhuatl and Citlaltépetl glaciers. This database has the potential for providing interesting and useful information on natural and anthropogenic-induced changes related to the occurrence of heavy metals in tropical glaciers in the northern hemisphere of North America. Black carbon concentrations analyzed in snow and glacier ice, and preliminary data on stable isotopes of Zn, also add information on natural vs. anthropogenic sources of heavy metals in central Mexico.

Influence of New Parameterization Schemes on Arctic Sea Ice Simulation

Two coupled climate models that participated in the CMIP6 project (Coupled Model Intercomparison Project Phase 6), the Earth System Model of Chinese Academy of Sciences version 2 (CAS-ESM2-0), and the Nanjing University of Information Science and Technology Earth System Model version 3 (NESM3) were assessed in terms of the impact of four new sea ice parameterization schemes. These four new schemes are related to air–ice heat flux, radiation penetration and absorption, melt ponds, and ice–ocean flux, respectively. To evaluate the effectiveness of these schemes, key sea ice variables with and without these new schemes, such as sea ice concentration (SIC) and sea ice thickness (SIT), were compared against observation and reanalysis products from 1980 to 2014. The simulations followed the design of historical experiments within the CMIP6 framework. The results revealed that both models demonstrated improvements in simulating Arctic SIC and SIT when the new parameterization schemes were implemented. The model bias of SIC in some marginal sea ice zones of the Arctic was reduced, especially during March. The SIT was increased and the transpolar gradient of SIT was reproduced. The changes in spatial patterns of SIC and SIT after adding new schemes bear similarities between the two coupled models. This suggests that the new schemes have the potential for broad application in climate models for simulation and future climate scenario projection, especially for those with underestimated SIT.

An Improved NLCS Algorithm Based on Series Reversion and Elliptical Model Using Geosynchronous Spaceborne–Airborne UHF UWB Bistatic SAR for Oceanic Scene Imaging

Geosynchronous spaceborne–airborne (GEO-SA) ultra-high-frequency ultra-wideband bistatic synthetic aperture radar (UHF UWB BiSAR) provides high-precision images for marine and polar environments, which are pivotal in glacier monitoring and sea ice thickness measurement for polar ocean mapping and navigation. Contrasting with traditional high-frequency BiSAR, it faces unique challenges, such as the considerable spatial variability, significant range–azimuth coupling, and vast volumes of echo data, which impede high-resolution image reconstruction. This paper presents an improved bistatic nonlinear chirp scaling (NLCS) algorithm for imaging oceanic scenes with GEO-SA UHF UWB BiSAR. This methodology extends the two-dimensional (2-D) spectrum up to the sixth order via the method of series reversion (MSR) to meet accuracy demands and then employs an elliptical model to elucidate the alterations in the azimuth frequency modulation (FM) rate mismatch. Initially, the imaging geometry and signal model are introduced, and then a separation of bistatic slant ranges based on the configuration is proposed. In addition, during range processing, after eliminating linear range cell migration (RCM), the derivation process for the sixth-order 2-D spectrum is detailed and an improved filter is applied to correct the high-order RCM. Finally, during azimuth processing, the causes of the FM rate mismatch are analyzed, a cubic perturbation function derived from the elliptical model is used for FM rate equalization, and a unified sixth-order filter is applied to complete the azimuth compression. Experimental results with point targets and natural oceanic scenes validate the outstanding efficacy of the proposed NLCS algorithm, particularly in imaging quality enhancements for GEO-SA UHF UWB BiSAR.

The use of terrestrial monoscopic time-lapse cameras for surveying glacier flow velocity

Glacier flow velocity is crucial information to assess the impact of global warming on glaciers since it is related to the ice thickness, its variations provide information on mass balance and, in general, on the current state of glacier “health”. Moreover, flow anomalies are often an indicator of glacier instabilities. Therefore, more attention has been dedicated to surveying glacier kinematics. Historically, flow velocity was the first quantitative variable measured on glaciers since the 19th century. In the last decades, terrestrial monoscopic time-lapse digital cameras have permitted to conduct automatic surveying at high spatial and temporal resolutions using digital image correlation. Even though terrestrial time-lapse imagery is currently a consolidated technique in glacier monitoring, possible strategies, limitations and potentialities have never been systematically reviewed. This work aims to illustrate the typical procedures required to monitor glacier flow velocity using terrestrial monoscopic time-lapse cameras. We describe possible inconveniences that can arise during the survey and provide some guidelines to minimise such issues. We present six study cases in the European Alps (Mont Blanc and Bernina massifs) that feature different monitoring equipment, site geometry and glacier morphodynamics to illustrate possible solutions for terrestrial imagery monitoring. The results of this review highlight the high benefit-to-cost ratio of terrestrial time-lapse cameras in glacier flow surveying.

The foundations of the Patagonian icefields

The two vast Patagonian icefields are a global hotspot for ice-loss. However, not much is known about the total ice volume they store - let alone its spatial distribution. One reason is that the abundant record of direct thickness measurements has never been systematically exploited. Here, this record is combined with remotely-sensed information on past ice thickness mapped from glacier retreat. Both datasets are incorporated in a state-of-the-art, mass-conservation approach to produce a well-informed map of the basal topography beneath the icefields. Its major asset is the reliability increase of thicknesses values along the many marine- and lake-terminating glaciers. For these, frontal ice-discharge is notably lower than previously reported. This finding implies that direct climatic control was more influential for past ice loss. We redact a total volume for both icefields in 2000 of 5351 km3. Despite the wealth of observations used in this assessment, relative volume uncertainties remain elevated.

Ice volume and thickness of all Scandinavian glaciers and ice caps

Abstract We present a new map of bed topography and ice thickness together with a corresponding ice volume estimate representative of the years ~2010 for all Scandinavian ice caps and glaciers. Starting from surface observations, we invert for ice thickness by iteratively running an innovative ice dynamics model on a distributed grid and updating bed topography until modelled and observed glacier dynamics as represented by their rate of surface elevation change (dh/dt) fields align. The ice flow model used is the instructed glacier model (Jouvet and Cordonnier, 2023, Journal of Glaciology 1–15), a generic physics-informed deep-learning emulator that models higher-order ice flow with high-computational efficiency. We calibrate the modelled thicknesses against >11 000 ice thickness observations, resulting in a final ice volume estimate of 302.7 km3 for Norway, 18.4 km3 for Sweden and 321.1 km3 for the whole of Scandinavia with an error estimate of ~ $\pm 11\%$ . The validation statistics computed indicate good agreement between modelled and observed thicknesses (RMSE = 55 m, Pearson's r = 0.87, bias = 0.8 m), outperforming all other ice thickness maps available for the region. The modelled bed shapes thus provide unprecedented detail in the subglacial topography, especially for ice caps where we produce the first maps that show ice-dynamically realistic flow features.

CD-type discretization for sea ice dynamics in FESOM version 2

Abstract. Two recently proposed variants of CD-type discretizations of sea ice dynamics on triangular meshes are implemented in the Finite-VolumE Sea ice–Ocean Model (FESOM version 2). The implementations use the finite element method in spherical geometry with longitude–latitude coordinates. Both are based on the edge-based sea ice velocity vectors but differ in the basis functions used to represent the velocities. The first one uses nonconforming linear (Crouzeix–Raviart) basis functions, and the second one uses continuous linear basis functions on sub-triangles obtained by splitting parent triangles into four smaller triangles. Test simulations are run to show how the performance of the new discretizations compares with the A-grid discretization using linear basis functions. Both CD discretizations are found to simulate a finer structure of linear kinematic features (LKFs). Both show some sensitivity to the representation of scalar fields (sea ice concentration and thickness). Cell-based scalars lead to a finer LKF structure for the first CD discretization, but the vertex-based scalars may be advantageous in the second case.

The linkage between wintertime sea ice drift and atmospheric circulation in an Arctic ice-ocean coupled simulation

By analyzing an Arctic ice-ocean coupled simulation, we study the linkage between wintertime sea ice drift and atmospheric circulation, and interpret the driving force terms in the sea ice dynamic equation. Sea ice drift anomaly is featured by an anticyclonic (cyclonic) gyre when regulated by negative (positive) phase of Arctic Oscillation with positive (negative) phase of Arctic Dipole, and a quasi-meridional stream from Chukchi-Beaufort (Barents-Kara) Seas to Barents-Kara (Chukchi-Beaufort) Seas when regulated by positive (negative) phase of Arctic Oscillation with positive (negative) phase of Arctic Dipole. Sea ice drift anomaly, when regulated by the mode alone, resembles spatial pattern of leading atmospheric mode. Decomposing sea ice dynamical equation shows that wind-ice stress dominates sea ice drift in areas away from islands and continental coastlines, ocean-ice stress acts as a resistant power to partly cancel the wind-ice stress in these areas, while in the coastal areas such as the thick multiyear ice zone north of the Canadian Arctic Archipelago, the wind-ice and ocean-ice stresses are small, the balance exists between sea surface height potential gradient and internal ice stress divergence. Developing more sophisticated internal ice stress expression in ice model is of great important to correctly project future sea ice change for the ice modeling community.

Результаты определения профиля высот геоида и уклонения отвесной линии по сигналам ГНСС, отраженным от водной поверхности

The method for determining geoid elevations, based on the use of GNSS signals reflected from water surface, is called bistatic one. Its features are small dimensions and weight of the necessary on-board equipment, and that enables mounting them on a nano-satellite, as well as simultaneously determining several dozen geoid height profiles. It corresponds to the number of visible navigation satellites from among GLONASS, GPS, Galileo, BeiDou. Several space projects using the technology under description were implemented abroad, their analysis has shown the possibility of determining geoid heights with an error of 10–20 cm, as well as sea ice thickness. The authors provide a review of some similar projects, an assessment of the simultaneously determined number of geoid height profiles, and calculation of its elevations and deflection of vertical along a single profile based on primary measurements from the project of CYGNSS. To compensate for various interfering parameters when calculating geoid elevations, the following models were used

Limited increases in Arctic offshore oil and gas production with climate change and the implications for energy markets

Climate change impacts on sea ice thickness is opening access to offshore Arctic resources. The degree to which these resources are exploited will depend on sea-ice conditions, technology costs, international energy markets, and the regulatory environment. We use an integrated human-Earth system model, GCAM, to explore the effects of spatial–temporal patterns of sea-ice loss under climate change on future Arctic offshore oil and gas extraction, considering interactions with global energy markets and emission reduction scenarios. We find that under SSP5, a “fossil-fueled development” scenario, the effects of sea-ice loss are larger for Arctic offshore oil production than gas. Under SSP5, future extraction of Arctic offshore oil and gas through 2100 adds roughly 0.8–2.6 EJ/year to oil and gas markets but does not have large impacts on global oil and gas markets. Surprisingly, a low-carbon scenario results in greater Arctic offshore oil production to offset the more emissions-intensive unconventional oil production.

A Study on the Ice Resistance Characteristics of Ships in Rafted Ice Based on the Circumferential Crack Method

In previous studies of ship–ice interactions, most studies focused on ship–level ice interactions, overlooking potential rafted ice conditions in extreme ice conditions. The purpose of this study is to develop a numerical model for predicting ship resistance in rafted ice regions. Numerical modeling of rafted ice was carried out using preset grid cells. By comparing the model test results, the accuracy and reliability of the numerical model are verified. On this basis, we undertook the analysis of the impacts of different ice thicknesses, ship speeds, bending strengths, and crushing strengths on the ice resistance of ships under level and rafted ice conditions. The results show that the ice resistance of ships is significantly higher than that of rafted ice under the condition of level ice; however, level ice and rafted ice have different effects on ship ice resistance. Compared with level ice, the ice resistance of ships navigating in rafted ice is more concentrated. The findings of the present research can serve as a technical reference for studies focused on predicting ship resistance in rafted ice regions.

The impacts of low flow, ice‐cover and ice thickness on sediment load in a sub‐arctic river – Modelling sediment transport with particle image velocimetry calibration data sets

AbstractClimate change will be pronounced in sub‐arctic and arctic regions. Consequently, a shorter ice‐cover period, increased precipitation and changes in the timing and magnitude of discharge are expected. These hydroclimatic changes can have an impact on sediment transport in northern rivers. Studying spatial and temporal variation of the processes in low‐flow open‐water and ice‐covered conditions is crucial to improve the prediction of future changes. The investigation of under‐ice conditions has been challenging and new approaches for their measurement are much needed. We analyse the spatial sediment transport in a meandering subarctic river during different flow conditions, i.e., low flow open‐channel during autumn and ice‐covered during winter. We use one‐dimensional (1D) morphodynamic models with image‐based sediment transportation data sets. We also simulate sediment transport with different ice thicknesses to better understand the overall bedload transport in a subarctic river in mid‐winter conditions.Simulations for the studied meander bend showed a considerable flow‐driven decrease of sediment transport in ice‐covered conditions when compared to open‐water conditions. Erosion was more pronounced in open‐water conditions and deposition was the prevailing process in ice‐covered conditions. When an additional increase in ice thickness was simulated, bedload increased, because the thicker ice narrowed the river channel substantially and thus leading to an increase in flow velocity and shear forces. Instead of solely relying on the sediment samples, the additional consideration of an image‐based sediment data set enabled a more reliable model calibration in ice‐covered conditions. This encourages further usage of image‐based methods for sediment transport estimations, especially in ice‐covered conditions. In the future, changes in river‐ice can be expected to decrease under‐ice sediment load but increase annual open‐water sediment loads. The increasing summative sediment transport of small arctic and subarctic rivers can have significant consequences on the downstream waterbodies.

The roles of river discharge and sea ice melting in formation of freshened surface layers in the Kara, Laptev, and East Siberian seas

Wide areas of the Siberian Arctic shelf are covered by freshened surface water layers, which are among the largest in the World Ocean. River discharge is the main freshwater source for formation of these layers; therefore, they are commonly referred to as river plumes (the Ob-Yenisei plume in the Kara Sea and the Lena plume in the Laptev and East Siberian seas). The contribution of sea ice meltwater (SIM) to the Ob-Yenisei and Lena plumes is pointed out to be small, albeit its actual volume remains unknown. In this study, we use a novel dataset of satellite-derived sea ice thickness in the Arctic Ocean during the melt period to quantify the annual volume of SIM, which was received by the Ob-Yenisei and Lena plumes during 2012–2020. We reveal that SIM is a significant source for the Lena plume providing, on average, 20% of total annual freshwater content. Moreover, the share of SIM in the Lena plume shows large inter-annual (14%–29%) variability, i.e., during certain years, SIM provides almost one-third of freshwater volume of the Lena plume. This variability is governed by inter-annual variability of ice thickness, as well as seasonal variability of sea ice melting conditions. Conversely, the contribution of SIM to the Ob-Yenisei plume is relatively low (8% on average), and its total annual share varies from 6% to 11% during the study period. This difference is mainly caused by significantly smaller area of the Ob-Yenisei plume as compared with the Lena plume. The forecasted earlier onset of ice melting in the Arctic Ocean in future decades due to climate change could decrease the contribution of SIM to the Ob-Yenisei plume, whereas its influence on the Lena plume remains unclear.

Sea ice cover in the Copernicus Arctic Regional Reanalysis

Abstract. The Copernicus Arctic Regional Reanalysis (CARRA) is a novel regional high-resolution atmospheric reanalysis product that covers a considerable part of the European Arctic including substantial amounts of ice-covered areas. Sea ice in CARRA is modelled by means of a one-dimensional thermodynamic sea ice parameterisation scheme, which also explicitly resolves the evolution of the snow layer over sea ice. In the present study, we assess the representation of sea ice cover in CARRA and validate it against a wide set of satellite products and observations from ice mass balance buoys. We show that CARRA adequately represents general interannual trends towards thinner and warmer ice in the Arctic. Compared to ERA5, sea ice in CARRA shows a reduced warm bias in the ice surface temperature. The strongest improvement was observed for winter months over the central Arctic and the Greenland and Barents seas where a 4.91 °C median ice surface temperature error in ERA5 is reduced to 1.88 °C in CARRA on average. Over Baffin Bay, intercomparisons suggest the presence of a cold winter-time ice surface temperature bias in CARRA. No improvement over ERA5 was found in the ice surface albedo with spring-time errors in CARRA being up to 0.08 higher on average than those in ERA5 when computed against the CLARA-A2 satellite retrieval product. Summer-time ice surface albedos are comparable in CARRA and ERA5. Sea ice thickness and snow depth in CARRA adequately resolve the annual cycle of sea ice cover in the Arctic and bring added value compared to ERA5. However, limitations of CARRA indicate potential benefits of utilising more advanced approaches for representing sea ice cover in next-generation reanalyses.

Ice Thickness Assessment of Non-Freshwater Lakes in the Qinghai–Tibetan Plateau Based on Unmanned Aerial Vehicle-Borne Ice-Penetrating Radar: A Case Study of Qinghai Lake and Gahai Lake

Ice thickness has a significant effect on the physical and biogeochemical processes of a lake, and it is an integral focus of research in the field of ice engineering. The Qinghai–Tibetan Plateau, known as the Third Pole of the world, contains numerous lakes. Compared with some information, such as the area, water level, and ice phenology of its lakes, the ice thickness of these lakes remains poorly understood. In this study, we used an unmanned aerial vehicle (UAV) with a 400/900 MHz ice-penetrating radar to detect the ice thickness of Qinghai Lake and Gahai Lake. Two observation fields were established on the western side of Qinghai Lake and Gahai Lake in January 2019 and January 2021, respectively. Based on the in situ ice thickness and the propagation time of the radar, the accuracy of the ice thickness measurements of these two non-freshwater lakes was comprehensively assessed. The results indicate that pre-processed echo images from the UAV-borne ice-penetrating radar identified non-freshwater lake ice, and we were thus able to accurately calculate the propagation time of radar waves through the ice. The average dielectric constants of Qinghai Lake and Gahai Lake were 4.3 and 4.6, respectively. This means that the speed of the radar waves that propagated through the ice of the non-freshwater lake was lower than that of the radio waves that propagated through the freshwater lake. The antenna frequency of the radar also had an impact on the accuracy of ice thickness modeling. The RMSEs were 0.034 m using the 400 MHz radar and 0.010 m using the 900 MHz radar. The radar with a higher antenna frequency was shown to provide greater accuracy in ice thickness monitoring, but the control of the UAV’s altitude and speed should be addressed.

On the data-driven investigation of factors affecting the need for icebreaker assistance in ice-covered waters

Merchant vessels often require icebreaker (IB) assistance to create safe pathways and improve efficiency when navigating in the Baltic Sea. Since IB resources are limited, an accurate estimation on the need for IB assistance is important. Whether IB assistance is needed depends on multiple factors. While practical experience from captains is naturally a source of valuable information for the decision on the need for IB assistance, systematic analysis of the reasoning is limited. The primary aim of this paper is to holistically investigate the influencing factors and their effect on estimating the need for IB assistance through data-driven techniques. Based on a comprehensive list of potential factors, different of data such as traffic history, environmental conditions, and ship specifications are gathered to present complex navigational scenarios. Each scenario is labeled by different navigation modes (independent navigation or IB assistance), laying the foundation for influencing factor identification and effect quantification. Logistic regression is applied to evaluate the effect of the factors on the need for IB assistance. The results show that the impact of the factors is diverse, and ridged ice concentration has the most significant impact. The effectiveness of identified factors is measured by comparing it to that of the factors that have been implemented by the existing studies (e.g., the combination of ice concentration, thickness, and ship ice class, or only ship speed). By considering the factors in this study, the classification performance can be improved by at least 5.6%. The findings in this paper can provide insights for predicting IB workloads and optimizing IB resources and have the potential to support the development of an intelligent decision-support system for winter navigation.