Ice Roughness Research Articles

The Determination of Criticality for Ice Shapes Based on CCAR-25

Determining the criticality of ice shapes is a necessary condition for verifying compliance with icing airworthiness regulations. However, the clear, concise, and applicable criterion based on the geometric characteristics of ice shapes has not been clearly given out by current advisory circulars. To address this problem, this paper summarizes aerodynamic performance items and recommended ice shapes the latest version of CCAR-25 and corresponding advisory circulars for a variety of flight phases, including takeoff, holding, en route, DTO, etc., instead of the single phase of holding in the previous research. Based on the geometric classification of the ice shapes, the dominant parameters of various ice shapes are clarified by the correlation between the geometric parameters and aerodynamic effects. The geometric parameters to determine the criticality of specific ice shapes are defined as the roughness height and range for the roughness ice and the total projection height in the direction of lift for the horn ice. On this basis, the detailed determination criterion of critical ice shape geometries corresponding to different flight phases and aircraft components is formulated, which will provide an operational selection methodology for determining the geometries of critical ice shapes at the airworthiness certification stage.

River ice breakup classification using dual- (HH&HV) or compact-polarization RADARSAT Constellation Mission data

Ice jams and associated flooding during river ice breakup are seasonal hazards for many northern communities. Synthetic Aperture Radar (SAR) data enable timely monitoring of river ice conditions in inclement weather. River ice types and open water are discriminated during breakup using dual-polarized (HH&HV; DPH) and compact polarimetric (CP) C-band SAR imagery from the RADARSAT Constellation Mission (RCM). Study areas comprise nineteen locations on rivers in northern Ontario, the Northwest Territories, and Alberta, Canada. Rubble ice, sheet ice, and open water areas are sampled in RCM imagery using corroborating evidence from shoreline cameras and oblique airborne photography. Samples are obtained in the incidence angle range 19° to 48°, for open water with various wind speeds and flow conditions, and for rubble ice and sheet ice with varying surface roughness, snow cover, and wetness. Discrimination relies on a primary classification of rubble ice, sheet ice, and open water, and a secondary classification of ice roughness. The primary classification model development uses a recursive-partitioning machine-learning technique. Overall accuracies for the best DPH models are 83.8% to 89.6%, depending on image noise floor. DPH models use both HH and HV polarizations for low noise floor image modes, whereas only HH is used for higher noise floor image modes. The best CP model accuracy is 90.2%, using the RL, RR, and RVRH parameters. Secondary classification associates increasing ice roughness with increasing HH backscatter for DPH data, and with increasing RH backscatter for CP data. The angular dependencies of rubble ice or sheet ice are used to normalize backscatter, which is then divided into roughness categories. The DPH and CP classification models are named IceBC-DP and IceBC-CP, respectively. Qualitative analysis illustrates good overall ice type and open water classification. Confounding situations are mixed pixel issues, whitewater and wind roughening of water, moist snow microwave absorption, and superficial water on sheet ice. The development of robust methods for monitoring ice jam formation with RCM is an operational concern for departments within the Government of Canada and within Provincial and Territorial Governments.

Parametric investigation of aerodynamic performance degradation due to icing on a symmetrical airfoil

Ice accretion on lifting surfaces induces an aerodynamic penalty in lift and drag on an aircraft. This performance degradation depends on the geometric features, type, and surface characteristics of the accreted ice on the airfoil. In the present work, we propose a set of two-parameter, low-order models to represent some of the typical ice topologies: glaze, rime, and horn. The parametric space is swept for all types of ice to isolate the aerodynamic changes causing performance degradation on a canonical symmetrical airfoil, which is the representative airfoil used by the National Research Council of Canada's platform for ice accretion and coatings tests with ultrasonic readings platform for in-flight icing tests. The three ice topologies show a self-similar trend between the stall angle of attack and the ice thickness, with the horn-type of ice imparting the greatest drag and lift penalty due to strong boundary layer separation. The relative effect of ice roughness plays a secondary role in performance degradation, and in some cases, the roughness causes a thicker and more resilient boundary layer, which can, under very specific icing conditions, enhance the aerodynamic performance.

Rough ice cover modified hyporheic exchange: A numerical study on the mechanical effect propagating from the ice-water interface to the sediment–water interface

Many rivers are subject to ice-covered flow conditions in cold climate regions or during the winter months. The presence of ice cover can change stream flow hydrodynamics and thus affect the hyporheic exchange involved in many biogeochemical processes in aquatic environments. Among the various features of ice cover, its underside roughness is a main factor modifying the surface–subsurface flow hydrodynamics. To date, there has been little research on how ice cover affects the hyporheic exchange, even less on the effect of the ice cover roughness. This leaves open questions on the mechanical changes initiated at the ice-water interface and transferred to the water–sediment interface. In this study, a coupled surface–subsurface computational approach was used to model conjunctively the channel flow over dune bed and the underlying porous flow under different ice cover conditions. The results show that ice cover can enhance the hyporheic flux and reduce the hyporheic zone depth when compared to open channel flow at equivalent discharge. As the river flow depth increases to 4 times the height of the dune, the effect of the ice cover on the hyporheic flow becomes negligible. Within this effective range, the hyporheic flux via the channel water depth follows a power-law function. Ice cover roughness greatly augments the exchange rate while compresses the exchange space, therefore reduces flow residence time in the hyporheic zone. Based on the modeling results, prediction models were proposed for evaluating the impacts of ice cover roughness on the hyporheic exchange flux and depth. Our study indicates that hyporheic flow exchange is significantly impacted by the roughness of the ice cover which in turn the hyporheic exchange is likely to affect the river ice processes as well as the water quality and ecological functions throughout the river corridors.

Deformation of vegetated channel bed under ice-covered flow conditions

A common feature of cold regions is the presence of ice cover on water surface. In the presence of vegetation in the channel bed which is normally leafless in winter, interaction between river ice and vegetated channel bed becomes very complex. In the present study, the impacts of ice cover and submerged leafless vegetation on channel bed deformation and flow resistance have been investigated based on 144 experiments conducted in a large-scale outdoor flume. The independent variables associated with the maximum scour depth around vegetation elements have been assessed and equations have been developed. Results indicate that the most important variable affecting the maximum scour depth around vegetation under ice-covered flow conditions is the ratio of the ice cover roughness to the bed roughness (ni/nb). However, under open channel flow conditions, the flow Froude number is the most important variable influencing the maximum scour depth. The methods based on the law of the wall provide reliable Manning's coefficient for both smooth and rough ice covers in the presence of submerged vegetation in the channel bed. As the vegetation density increases, the size of scour holes becomes smaller. With the increase in vegetation density, the median grain size of the armour layer in scour holes decreases correspondingly. When vegetation elements are placed in a staggered configuration, the median grain size of the armour layer in scour holes is less than that of squared configuration of vegetation elements. Regardless of the roughness of ice cover on water surface and the grain size of sediment in the bed, the maximum depth of scour holes always occurs at the upstream front face of each vegetation element. Under the rough ice-covered flow condition, the maximum depth of scour holes is more than that under the smooth ice-covered and open flow conditions.

Characteristics of turbulent flow around bridge abutments in the presence of vegetation in channel bed under ice-covered flow conditions

During winter, ice cover frequently forms on the water surface of rivers with vegetated channel beds in cold regions. The investigation of the impacts of both ice cover and channel bed vegetation on flow structures around bridge abutments is essential for engineers to gain a comprehensive understanding of the complex interactions occurring in such a situation. In the present study, the flow structure around a rectangular bridge abutment in the presence of vegetation under ice-covered conditions has been studied. Considering different vegetation densities by arranging vegetation elements in square and staggered configurations, this study incorporates the influence of ice covers with different roughness, namely smooth and rough ice cover. Key turbulence parameters, including turbulence intensity, Reynolds shear stress (RSS), and turbulent kinetic energy (TKE), are also examined based on laboratory experiments. Results show that the shape of velocity profiles for flow over a vegetated channel bed changes from an S-shaped curve under an open flow condition to a convex shape under ice-covered conditions. The effects of an ice cover and vegetation on the flow around the bridge abutment create unpredictable turbulence intensity patterns. Under a rough-covered flow condition, there appears a larger area with negative Reynolds shear stress (RSS) downstream of the abutment. Turbulent kinetic energy (TKE) under ice-covered conditions has substantially lower magnitudes than in open flow conditions.

Arctic Sea Ice Topography Information From RADARSAT Constellation Mission (RCM) Synthetic Aperture Radar (SAR) Backscatter

AbstractSea ice topography information can be obtained from altimetry but these data are spatially and temporally limited compared to recent synthetic aperture radar (SAR) missions such as the RADARSAT Constellation Mission (RCM). We analyze the relationship between sea ice roughness and height obtained from three Ice, Cloud, and Land Elevation Satellite (ICESat)‐2 tracks on two dates in March 2022, with RCM backscatter from 17 images in the McClintock Channel, Canadian Arctic. We analyze how this relationship varies with ice type, polarization, and incidence angle. We find particularly notable relationships between sea ice roughness and horizontal‐transmit/vertical‐receive backscatter for first‐year ice, and sea ice height and backscatter for multi‐year ice. We develop a preliminary model for winter sea ice roughness retrieval using RCM. In comparison with independent ICESat‐2 data in our study region, we find the model performs effectively at estimating a roughness distribution and key roughness statistics, and characterizes spatial variations in roughness at a sub‐kilometer scale.

Influence of Surface Ice Roughness on the Aerodynamic Performance of Wind Turbines

Influence of Surface Ice Roughness on the Aerodynamic Performance of Wind Turbines

Simulating Radiative Heat Transfer in Multi‐Scattering Irregular Surfaces: Application to Snow and Ice Morphologies on Europa

AbstractWe developed a Monte‐Carlo‐based radiative heat transfer model capable of simulating solar exposure and subsequent warming of rough snow and ice surfaces on ice‐covered airless solar system bodies. The model accounts for wavelength‐dependent internal light scattering and heat conduction in the snow interior down to meter‐scale depths. We validated the model against analytical and experimental test cases with relevant applications to Europa, one of Jupiter's moons. We examined differential heating across the surface, from the centimeter to meter scale, to reveal potential patterns of preferential sublimation that could lead to rough ice morphologies, such as penitentes. An exploration of parameters such as penitente height‐width ratios, shape, size, snow grain size, and thermal properties revealed that taller, thinner, larger penitentes with sharper peaks, coarser snow grain sizes, and lower thermal inertias are more likely to grow in Europa's environment near the equator.

Investigating Future Arctic Sea Ice Loss and Near‐Surface Wind Speed Changes Related to Surface Roughness Using the Community Earth System Model

AbstractThe Arctic is undergoing a pronounced and rapid transformation in response to changing greenhouse gasses, including reduction in sea ice extent and thickness. There are also projected increases in near‐surface Arctic wind. This study addresses how the winds trends may be driven by changing surface roughness and/or stability in different Arctic regions and seasons, something that has not yet been thoroughly investigated. We analyze 50 experiments from the Community Earth System Model Version 2 (CESM2) Large Ensemble and five experiments using CESM2 with an artificially decreased sea ice roughness to match that of the open ocean. We find that with a smoother surface there are higher mean wind speeds and slower mean ice speeds in the autumn, winter, and spring. The artificially reduced surface roughness also strongly impacts the wind speed trends in autumn and winter, and we find that atmospheric stability changes are also important contributors to driving wind trends in both experiments. In contrast to the clear impacts on winds, the sea ice mean state and trends are statistically indistinguishable, suggesting that near‐surface winds are not major drivers of Arctic sea ice loss. Two major results of this work are: (a) the near‐surface wind trends are driven by changes in both surface roughness and near‐surface atmospheric stability that are themselves changing from sea ice loss, and (b) the sea ice mean state and trends are driven by the overall warming trend due to increasing greenhouse gas emissions and not significantly impacted by coupled feedbacks with the surface winds.

Shear speed in Arctic ice and underwater acoustic reflection

The underwater reflected sound wave under a sheet of Arctic ice contains information about the ice including thickness and mechanical properties. The ideal case of a perfectly flat ice floe floating above the water may be modeled using the OASES wavenumber integration code. A spectrogram of the acoustic response contains features related to the compressional and shear waves. The resonance associated with the shear wave speed in the ice is particularly distinctive. In the real world, the ice-water interface is not perfectly smooth. The finite element code SPECFEM2D is used to simulate the response of ice with a rough ice–water interface. It is particularly well suited to this problem since it accounts for all orders of multiple reflections. It shows the effects of fine-scale roughness on the acoustic response. Depending on the severity of the roughness, it may enhance or diminish the acoustic features related to the properties of the ice. [Work supported by the U. S. Office of Naval Research, Code 32 Grant N00014-20-1-2041 (N. Chotiros). G. Bayrakci and A. Best were funded by the UK Defence and Security Accelerator, Grant ACC2016927. The Texas Advanced Computing Center provided the high-performance computing resources.]

Turbulent structure of the Arctic boundary layer in early summer driven by stability, wind shear and cloud-top radiative cooling: ACLOUD airborne observations

Abstract. Clouds are assumed to play an important role in the Arctic amplification process. This motivated a detailed investigation of cloud processes, including radiative and turbulent fluxes. Data from the aircraft campaign ACLOUD were analyzed with a focus on the mean and turbulent structure of the cloudy boundary layer over the Fram Strait marginal sea ice zone in late spring and early summer 2017. Vertical profiles of turbulence moments are presented from contrasting atmospheric boundary layers (ABLs) from 4 d. They differ by the magnitude of wind speed, boundary-layer height, stability, the strength of the cloud-top radiative cooling and the number of cloud layers. Turbulence statistics up to third-order moments are presented, which were obtained from horizontal-level flights and from slanted profiles. It is shown that both of these flight patterns complement each other and form a data set that resolves the vertical structure of the ABL turbulence well. The comparison of the 4 d shows that especially during weak wind, even in shallow Arctic ABLs with mixing ratios below 3 g kg−1, cloud-top cooling can serve as a main source of turbulent kinetic energy (TKE). Well-mixed ABLs are generated where TKE is increased and vertical velocity variance shows pronounced maxima in the cloud layer. Negative vertical velocity skewness points then to upside-down convection. Turbulent heat fluxes are directed upward in the cloud layer as a result of cold downdrafts. In two cases with single-layer stratocumulus, turbulent transport of heat flux and of temperature variance are both negative in the cloud layer, suggesting an important role of large eddies. In contrast, in a case with weak cloud-top cooling, these quantities are positive in the ABL due to the heating from the surface. Based on observations and results of a mixed-layer model it is shown that the maxima of turbulent fluxes are, however, smaller than the jump of the net terrestrial radiation flux across the upper part of a cloud due to the (i) shallowness of the mixed layer and (ii) the presence of a downward entrainment heat flux. The mixed-layer model also shows that the buoyancy production of TKE is substantially smaller in stratocumulus over the Arctic sea ice compared to subtropics due to a smaller surface moisture flux and smaller decrease in specific humidity (or even humidity inversions) right above the cloud top. In a case of strong wind, wind shear shapes the ABL turbulent structure, especially over rough sea ice, despite the presence of a strong cloud-top cooling. In the presence of mid-level clouds, cloud-top radiative cooling and thus also TKE in the lowermost cloud layer are strongly reduced, and the ABL turbulent structure becomes governed by stability, i.e., by the surface–air temperature difference and wind speed. A comparison of slightly unstable and weakly stable cases shows a strong reduction of TKE due to increased stability even though the absolute value of wind speed was similar. In summary, the presented study documents vertical profiles of the ABL turbulence with a high resolution in a wide range of conditions. It can serve as a basis for turbulence closure evaluation and process studies in Arctic clouds.

Sea Ice Elevation Measurements Using 3-D Laser Scanner

This study aims to introduce a sea ice elevation dataset estimated by using a 3-D laser scanner during the ice camp of the 2022 Arctic summer field survey. The equipment used is FARO’s Laser scanner FOCUS 3D X 130 HDR. The observed sea ice floe is located in the Arctic Ocean (76°13' N, 174°35') and is a multi-year ice with several melt ponds and ice ridges. We scanned eight sections separately, considering the equipment’s maximum horizontal scan range and the ice surface’s topographic features. The raw data were co-registered based on the positions of reference spheres. The result indicated a significant level of accuracy with a target-based vertical mean error of 4.8 mm. The laser scanner data were merged into point clouds ranging from 160×210 m. As a result, sea ice elevation data were generated in 0 to 2.8 m based on the minimum vertical point in the observation range. Sea ice elevation data is an essential variable in estimating the various properties of sea ice, such as ice thickness and roughness. In addition, using climatic variables such as air temperature and energy budget observed simultaneously can help to understand the physical interaction between the sea ice surface and the atmosphere on a local scale.

Improved Simulation of the Polar Atmospheric Boundary Layer by Accounting for Aerodynamic Roughness in the Parameterization of Surface Scalar Exchange Over Sea Ice

AbstractA new, simple parameterization scheme for scalar (heat and moisture) exchange over sea ice and the marginal ice zone is tested in a numerical weather and climate prediction model. This new “Blended A87” scheme accounts for the influence of aerodynamic roughness on the relationship between momentum and scalar exchange over consolidated sea ice, in line with long‐standing theory and recent field observations, and in contrast to the crude schemes currently operational in most models. Using aircraft observations and Met Office Unified Model simulations of cold‐air outbreak (CAO) conditions over aerodynamically rough sea ice, we demonstrate striking improvements in model performance when the Blended A87 scheme replaces the model's operational treatment for surface scalar exchange, provided that the aerodynamic roughness over consolidated ice is appropriately prescribed. The mean biases in surface sensible heat flux, surface latent heat flux, near‐surface air temperature, and surface temperature reduce from 25 to 11 W m−2, 22 to 12 W m−2, 0.8 to 0.0 K, and 1.4 to 0.8 K, respectively. We demonstrate that such impacts on surface exchange over sea ice can have a marked impact on the evolution of the atmospheric boundary layer across hundreds of kilometres downwind of the sea ice during CAO conditions in the model. Our results highlight the importance of spatiotemporal variability in the topography of consolidated sea ice for both momentum and scalar exchange over sea ice; accounting for which remains a key challenge for modeling polar weather and climate.

Heterogeneous melting near the Thwaites Glacier grounding line

Thwaites Glacier represents 15% of the ice discharge from the West Antarctic Ice Sheet and influences a wider catchment1–3. Because it is grounded below sea level4,5, Thwaites Glacier is thought to be susceptible to runaway retreat triggered at the grounding line (GL) at which the glacier reaches the ocean6,7. Recent ice-flow acceleration2,8 and retreat of the ice front8–10 and GL11,12 indicate that ice loss will continue. The relative impacts of mechanisms underlying recent retreat are however uncertain. Here we show sustained GL retreat from at least 2011 to 2020 and resolve mechanisms of ice-shelf melt at the submetre scale. Our conclusions are based on observations of the Thwaites Eastern Ice Shelf (TEIS) from an underwater vehicle, extending from the GL to 3 km oceanward and from the ice–ocean interface to the sea floor. These observations show a rough ice base above a sea floor sloping upward towards the GL and an ocean cavity in which the warmest water exceeds 2 °C above freezing. Data closest to the ice base show that enhanced melting occurs along sloped surfaces that initiate near the GL and evolve into steep-sided terraces. This pronounced melting along steep ice faces, including in crevasses, produces stratification that suppresses melt along flat interfaces. These data imply that slope-dependent melting sculpts the ice base and acts as an important response to ocean warming.

Automatic roughness characterization of simulated ice shapes

During in-flight ice accretion, roughness plays an important role since it heavily influences the convective heat transfer and skin friction coefficients. This paper aims to assess the ability of existing ice accretion simulation tools to compute the growing ice’s roughness. To this purpose, a technique based on Self Organizing Maps is applied to numerical simulations of in-flight ice accretion to characterize the roughness. The numerical ice predictions are performed using a standard approach comprising RANS computations, Lagrangian particle tracking, the solution of the unsteady Stefan problem, and a morphogenetic model. Numerical simulations are performed on selected benchmark cases from the 1st AIAA Ice Prediction Workshop. Validation of roughness computation is performed on synthetic test cases, while ice roughness is compared directly to that extracted from ice scans. The results of simulated ice shapes compare reasonably well with experimental data. Computations can replicate the trend of the experimental mean ice shape and roughness distribution for both rime and glaze cases.

Fast Ice Thickness Distribution in the Western Ross Sea in Late Spring

AbstractWe present a 700 km airborne electromagnetic survey of late‐spring fast ice and sub‐ice platelet layer (SIPL) thickness distributions from McMurdo Sound to Cape Adare, providing a first‐time inventory of fast ice thickness close to its annual maximum. The overall mode of the consolidated ice (including snow) thickness was 1.9 m, less than its mean of 2.6 ± 1.0 m. Our survey was partitioned into level and rough ice, and SIPL thickness was estimated under level ice. Although level ice, with a mode of 2.0 m and mean of 2.0 ± 0.6 m, was prevalent, rough ice occupied 41% of the transect by length, 50% by volume, and had a mode of 3.3 m and mean of 3.2 ± 1.2 m. The thickest 10% of rough ice was almost 6 m on average, inclusive of a 2 km segment thicker than 8 m in Moubray Bay. The thickest ice occurred predominantly along the northwestern Ross Sea, due to compaction against the coast. The adjacent pack ice was thinner (by ∼1 m) than the first‐year fast ice. In Silverfish Bay, offshore Hells Gate Ice Shelf, New Harbor, and Granite Harbor, the SIPL transect volume was a significant fraction (0.30) of the consolidated ice volume. The thickest 10% of SIPLs averaged nearly 3 m thick, and near Hells Gate Ice Shelf the SIPL was almost 10 m thick, implying vigorous heat loss to the ocean (∼90 W m−2). We conclude that polynya‐induced ice deformation and interaction with continental ice influence fast ice thickness in the western Ross Sea.

Observed and Parameterized Roughness Lengths for Momentum and Heat Over Rough Ice Surfaces

AbstractTurbulent heat fluxes, that is, the sensible heat flux and latent heat flux, are important sources/sinks of energy for surface melt over glaciers and ice sheets. Therefore, credible simulations of for example, future Greenland Ice Sheet mass loss need an accurate description of these fluxes. However, the parameterization of surface turbulent heat fluxes in climate models requires knowledge about the surface roughness lengths for momentum, heat and moisture, which are currently either unknown or tuned to indirect observations. In this study we take advantage of a large data set of eddy covariance observations acquired during multiple years and at multiple sites over the Greenland Ice Sheet. These in‐situ observations are used to develop an improved parameterization for the roughness length for momentum, and update the parameterization for the roughness lengths for heat and moisture over rough ice surfaces. The newly derived parameterizations are implemented in a surface energy balance model that is used to compute surface melt. Sensitivity experiments confirm the high sensitivity of surface melt to the chosen roughness length models. The new parameterization models the sensible heat flux to within 10 W m−2, and the cumulative ice ablation within 10% at three out of four sites. Two case studies demonstrate the important contribution of the turbulent heat fluxes to surface ablation. The presented roughness parameterizations can be implemented in climate models to improve the physical representation of surface roughness over rough snow and ice surfaces, which is expected to improve the modeled turbulent heat fluxes and thus surface melt.

Contrast Icing Wind Tunnel Tests between Normal Droplets and Supercooled Large Droplets

In order to compare and analyze the similarities and differences between normal droplet icing shapes and supercooled large droplet icing shapes, SADRI carried out normal droplet and supercooled large droplet icing wind tunnel tests in the NRC−AIWT icing wind tunnel. Taking the typical glaze ice in normal droplet icing conditions as the reference, the freezing drizzle and freezing rain icing tests under the supercooled large droplet conditions were carried out. The test results show that compared with normal droplets, the ice horn height of supercooled large droplets decreases with the increase in droplet particle size, and even the ice horn characteristics are not obvious when the icing condition is freezing rain. At the same time, the range and height of rough element ice shape after the main ice horn of supercooled large droplets are significantly larger and higher than those of the normal droplets, while the difference in the rough element in different supercooled large droplet icing conditions is small.

Mapping Arctic Sea-Ice Surface Roughness with Multi-Angle Imaging SpectroRadiometer

Sea-ice surface roughness (SIR) is a crucial parameter in climate and oceanographic studies, constraining momentum transfer between the atmosphere and ocean, providing preconditioning for summer-melt pond extent, and being related to ice age and thickness. High-resolution roughness estimates from airborne laser measurements are limited in spatial and temporal coverage while pan-Arctic satellite roughness does not extend over multi-decadal timescales. Launched on the Terra satellite in 1999, the NASA Multi-angle Imaging SpectroRadiometer (MISR) instrument acquires optical imagery from nine near-simultaneous camera view zenith angles. Extending on previous work to model surface roughness from specular anisotropy, a training dataset of cloud-free angular reflectance signatures and surface roughness, defined as the standard deviation of the within-pixel lidar elevations, from near-coincident operation IceBridge (OIB) airborne laser data is generated and is modelled using support vector regression (SVR) with a radial basis function (RBF) kernel selected. Blocked k-fold cross-validation is implemented to tune hyperparameters using grid optimisation and to assess model performance, with an R2 (coefficient of determination) of 0.43 and MAE (mean absolute error) of 0.041 m. Product performance is assessed through independent validation by comparison with unseen similarly generated surface-roughness characterisations from pre-IceBridge missions (Pearson’s r averaged over six scenes, r = 0.58, p < 0.005), and with AWI CS2-SMOS sea-ice thickness (Spearman’s rank, rs = 0.66, p < 0.001), a known roughness proxy. We present a derived sea-ice roughness product at 1.1 km resolution (2000–2020) over the seasonal period of OIB operation and a corresponding time-series analysis. Both our instantaneous swaths and pan-Arctic monthly mosaics show considerable potential in detecting surface-ice characteristics such as deformed rough ice, thin refrozen leads, and polynyas.