RADARSAT Geophysical Processor System Research Articles

On the sensitivity of sea ice deformation statistics to plastic damage

Abstract. We implement a plastic damage parametrization, distinct from the elastic damage in the elasto-brittle framework, in the standard viscous–plastic (VP) sea ice model to disentangle its effect from resolved model physics (visco-plastic with and without damage) on its ability to reproduce observed scaling laws of deformation. To this end, we compare scaling properties and multifractality of simulated divergence and shear strain rate, as proposed in the Sea Ice Rheology Experiment (SIREx) studies, with those derived from the RADARSAT Geophysical Processor System (RGPS). Results show that including a plastic damage parametrization in the standard viscous–plastic model increases the spatial but decreases the temporal localization of simulated linear kinematic features (LKFs) and brings all spatial deformation rate statistics in line with observations from RGPS without the need to increase the mechanical shear strength of sea ice as recently proposed for lower-resolution viscous–plastic sea ice models. In fact, including damage with a healing timescale of th=30 d and an increased mechanical strength unveils multifractal behavior that does not fit the theory. Therefore, a plastic damage parametrization is a powerful tuning knob affecting the deformation statistics of viscous–plastic sea ice.

Deformation lines in Arctic sea ice: intersection angle distribution and mechanical properties

Abstract. Despite its relevance for the Arctic climate and ecosystem, modeling sea-ice deformation, i.e., the opening, shearing, and ridging of sea ice, along linear kinematic features (LKFs) remains challenging, as the mechanical properties of sea ice are not yet fully understood. The intersection angles between LKFs provide valuable information on the internal mechanical properties, as they are linked to them. Currently, the LKFs emerging from sea-ice rheological models do not reproduce the observed LKF intersection angles, pointing to a gap in the model physics. We aim to obtain an intersection angle distribution (IAD) from observational data to serve as a reference for high-resolution sea-ice models and to infer the mechanical properties of the sea-ice cover. We use the sea-ice vorticity to discriminate between acute and obtuse LKF intersection angles within two sea-ice deformation datasets: the RADARSAT Geophysical Processor System (RGPS) and a new dataset from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) drift experiment. Acute angles dominate the IAD, with single peaks at 48∘±2 and 45∘±7. The IAD agrees well between both datasets, despite the difference in scale, time period, and geographical location. The divergence and shear rates of the LKFs also have the same distribution. The dilatancy angle (the ratio of shear and divergence) is not correlated with the intersection angle. Using the IAD, we infer two important mechanical properties of the sea ice: we found an internal angle of friction in sea ice of μI=0.66±0.02 and μI=0.75±0.05. The shape of the yield curve or the plastic potential derived from the observed IAD resembles a teardrop or a Mohr–Coulomb shape. With these new insights, sea-ice rheologies used in models can be adapted or redesigned to improve the representation of sea-ice deformation.

Sea Ice Rheology Experiment (SIREx): 1. Scaling and Statistical Properties of Sea‐Ice Deformation Fields

AbstractAs the sea‐ice modeling community is shifting to advanced numerical frameworks, developing new sea‐ice rheologies, and increasing model spatial resolution, ubiquitous deformation features in the Arctic sea ice are now being resolved by sea‐ice models. Initiated at the Forum for Arctic Modeling and Observational Synthesis, the Sea Ice Rheology Experiment (SIREx) aims at evaluating state‐of‐the‐art sea‐ice models using existing and new metrics to understand how the simulated deformation fields are affected by different representations of sea‐ice physics (rheology) and by model configuration. Part 1 of the SIREx analysis is concerned with evaluation of the statistical distribution and scaling properties of sea‐ice deformation fields from 35 different simulations against those from the RADARSAT Geophysical Processor System (RGPS). For the first time, the viscous‐plastic (and the elastic‐viscous‐plastic variant), elastic‐anisotropic‐plastic, and Maxwell‐elasto‐brittle rheologies are compared in a single study. We find that both plastic and brittle sea‐ice rheologies have the potential to reproduce the observed RGPS deformation statistics, including multi‐fractality. Model configuration (e.g., numerical convergence, atmospheric representation, spatial resolution) and physical parameterizations (e.g., ice strength parameters and ice thickness distribution) both have effects as important as the choice of sea‐ice rheology on the deformation statistics. It is therefore not straightforward to attribute model performance to a specific rheological framework using current deformation metrics. In light of these results, we further evaluate the statistical properties of simulated Linear Kinematic Features in a SIREx Part 2 companion paper.

Sea Ice Rheology Experiment (SIREx): 2. Evaluating Linear Kinematic Features in High‐Resolution Sea Ice Simulations

AbstractSimulating sea ice drift and deformation in the Arctic Ocean is still a challenge because of the multiscale interaction of sea ice floes that compose the Arctic Sea ice cover. The Sea Ice Rheology Experiment (SIREx) is a model intercomparison project of the Forum of Arctic Modeling and Observational Synthesis (FAMOS). In SIREx, skill metrics are designed to evaluate different recently suggested approaches for modeling linear kinematic features (LKFs) to provide guidance for modeling small‐scale deformation. These LKFs are narrow bands of localized deformation that can be observed in satellite images and also form in high resolution sea ice simulations. In this contribution, spatial and temporal properties of LKFs are assessed in 36 simulations of state‐of‐the‐art sea ice models and compared to deformation features derived from the RADARSAT Geophysical Processor System. All simulations produce LKFs, but only very few models realistically simulate at least some statistics of LKF properties such as densities, lengths, or growth rates. All SIREx models overestimate the angle of fracture between conjugate pairs of LKFs and LKF lifetimes pointing to inaccurate model physics. The temporal and spatial resolution of a simulation and the spatial resolution of atmospheric boundary condition affect simulated LKFs as much as the model's sea ice rheology and numerics. Only in very high resolution simulations (≤2 km) the concentration and thickness anomalies along LKFs are large enough to affect air‐ice‐ocean interaction processes.

Reassessing the Quality of Sea‐Ice Deformation Estimates Derived From the RADARSAT Geophysical Processor System and Its Impact on the Spatiotemporal Scaling Statistics

AbstractWe reassess the trajectory errors inherent to sea‐ice deformation estimates with a new propagation of uncertainty derivation and show that previous formulations applied to deformation estimates from the RADARSAT Geophysical Processor System (RGPS) are either too high due to incorrect assumptions or too low due to neglected terms in certain cases. We show that when the resulting signal‐to‐noise ratios are used to discriminate the deformation estimates based on their quality, as done for buoy records, the spatiotemporal scaling exponents for the mean total deformation rate increase, especially at smaller scale, such that a space‐time coupling of the scaling—which is otherwise absent—emerges from the RGPS deformation data set, in accord with previous analyses performed with buoy observations. We also show that the preprocessing method used to reduce the effects of irregular sampling of the Lagrangian deformation fields can significantly impact the value of the deformation statistics and could possibly explain part of previous discrepancies between deformation statistics obtained with buoy records and large‐scale synthetic aperture radar (SAR) imagery. Specifically, we show that spurious lines of deformation appear when interpolating RGPS trajectories that presenttemporal sampling inconsistencies. In the context of using observed sea‐ice deformation statistics to constrain and improve the performance of sea‐ice models, high confidence in the observed deformation field statistics is necessary. Using appropriate, well‐documented, methods to derive the set of statistics to be reproduced by models therefore becomes crucial.

Feature-based comparison of sea ice deformation in lead-permitting sea ice simulations

Abstract. The sea ice modeling community is progressing towards pan-Arctic simulations that explicitly resolve leads in the simulated sea ice cover. Evaluating these simulations against observations poses new challenges. A new feature-based evaluation of simulated deformation fields is introduced, and the results are compared to a scaling analysis of sea ice deformation. Leads and pressure ridges – here combined into linear kinematic features (LKFs) – are detected and tracked automatically from deformation and drift data. LKFs in two pan-Arctic sea ice simulations with a horizontal grid spacing of 2 km are compared with an LKF dataset derived from the RADARSAT Geophysical Processor System (RGPS). One simulation uses a five-class ice thickness distribution (ITD). The simulated sea ice deformation follows a multi-fractal spatial and temporal scaling, as observed from RGPS. The heavy-tailed distribution of LKF lengths and the scale invariance of LKF curvature, which points to the self-similar nature of sea ice deformation fields, are reproduced by the model. Interannual and seasonal variations in the number of LKFs, LKF densities, and LKF orientations in the ITD simulation are found to be consistent with RGPS observations. The lifetimes and growth rates follow a distribution with an exponential tail. The model overestimates the intersection angle of LKFs, which is attributed to the model's viscous-plastic rheology with an elliptical yield curve. In conclusion, the new feature-based analysis of LKF statistics is found to be useful for a comprehensive evaluation of simulated deformation features, which is required before the simulated features can be used with confidence in the context of climate studies. As such, it complements the commonly used scaling analysis and provides new useful information for comparing deformation statistics. The ITD simulation is shown to reproduce LKFs sufficiently well for it to be used for studying the effect of directly resolved leads in climate simulations. The feature-based analysis of LKFs also identifies specific model deficits that may be addressed by specific parameterizations, for example, a damage parameter, a grounding scheme, and a Mohr–Coulombic yield curve.

Leads and ridges in Arctic sea ice from RGPS data and a new tracking algorithm

Abstract. Leads and pressure ridges are dominant features of the Arctic sea ice cover. Not only do they affect heat loss and surface drag, but they also provide insight into the underlying physics of sea ice deformation. Due to their elongated shape they are referred to as linear kinematic features (LKFs). This paper introduces two methods that detect and track LKFs in sea ice deformation data and establish an LKF data set for the entire observing period of the RADARSAT Geophysical Processor System (RGPS). Both algorithms are available as open-source code and applicable to any gridded sea ice drift and deformation data. The LKF detection algorithm classifies pixels with higher deformation rates compared to the immediate environment as LKF pixels, divides the binary LKF map into small segments, and reconnects multiple segments into individual LKFs based on their distance and orientation relative to each other. The tracking algorithm uses sea ice drift information to estimate a first guess of LKF distribution and identifies tracked features by the degree of overlap between detected features and the first guess. An optimization of the parameters of both algorithms, as well as an extensive evaluation of both algorithms against handpicked features in a reference data set, is presented. A LKF data set is derived from RGPS deformation data for the years from 1996 to 2008 that enables a comprehensive description of LKFs. LKF densities and LKF intersection angles derived from this data set agree well with previous estimates. Further, a stretched exponential distribution of LKF length, an exponential tail in the distribution of LKF lifetimes, and a strong link to atmospheric drivers, here Arctic cyclones, are derived from the data set. Both algorithms are applied to output of a numerical sea ice model to compare the LKF intersection angles in a high-resolution Arctic sea ice simulation with the LKF data set.

Sea-ice deformation in a coupled ocean–sea-ice model and in satellite remote sensing data

Abstract. A realistic representation of sea-ice deformation in models is important for accurate simulation of the sea-ice mass balance. Simulated sea-ice deformation from numerical simulations with 4.5, 9, and 18 km horizontal grid spacing and a viscous–plastic (VP) sea-ice rheology are compared with synthetic aperture radar (SAR) satellite observations (RGPS, RADARSAT Geophysical Processor System) for the time period 1996–2008. All three simulations can reproduce the large-scale ice deformation patterns, but small-scale sea-ice deformations and linear kinematic features (LKFs) are not adequately reproduced. The mean sea-ice total deformation rate is about 40 % lower in all model solutions than in the satellite observations, especially in the seasonal sea-ice zone. A decrease in model grid spacing, however, produces a higher density and more localized ice deformation features. The 4.5 km simulation produces some linear kinematic features, but not with the right frequency. The dependence on length scale and probability density functions (PDFs) of absolute divergence and shear for all three model solutions show a power-law scaling behavior similar to RGPS observations, contrary to what was found in some previous studies. Overall, the 4.5 km simulation produces the most realistic divergence, vorticity, and shear when compared with RGPS data. This study provides an evaluation of high and coarse-resolution viscous–plastic sea-ice simulations based on spatial distribution, time series, and power-law scaling metrics.

Using sea‐ice deformation fields to constrain the mechanical strength parameters of geophysical sea ice

AbstractWe investigate the ability of viscous‐plastic (VP) sea‐ice models with an elliptical yield curve and normal flow rule to reproduce the shear and divergence distributions derived from the RADARSAT Geophysical Processor System (RGPS). In particular, we reformulate the VP elliptical rheology to allow independent changes in the ice compressive, shear and isotropic tensile strength parameters ( , respectively) in order to study the sensitivity of the deformation distributions to changes in the ice mechanical strength parameters. Our 10 km VP simulation with standard ice mechanical strength parameters = 27.5 kN m−2, = 6.9 kN m−2, and = 0 kN m−2 (ellipse aspect ratio of e = 2) does not reproduce the large shear and divergence deformations observed in the RGPS deformation fields, and specifically lacks well‐defined, active linear kinematic features (LKFs). Probability density functions (PDFs) for the shear and divergence of are nonetheless not Gaussian. Reducing the ice compressive strength (with constant S∗ and T∗) or increasing the ice shear strength (with constant P∗ and T∗) both results in shear and divergence PDFs in better agreement with RGPS distributions. The isotropic tensile strength of sea ice does not significantly affect the shear and divergence distributions. When considering additional metrics such as the ice drift error, mean ice thickness fields, and spatial scaling of the total deformations, our results suggest that reducing the ice compressive strength P∗ (while keeping S∗ constant, i.e. reducing the ellipse aspect ratio) is a better solution than increasing the shear strength to improve simulations of the Arctic sea‐ice cover with the VP elliptical rheology.

On producing sea ice deformation data sets from SAR-derived sea ice motion

Abstract. We propose a method to reduce the error generated when computing sea ice deformation fields from synthetic aperture radar (SAR)-derived sea ice motion. The method is based on two steps. The first step consists of using a triangulation of the positions taken from the sea ice trajectories to define a mesh on which a first estimate of sea ice deformation is computed. The second step consists of applying a specific smoother to the deformation field to reduce the artificial noise that arises along discontinuities in the sea ice motion field. This method is here applied to RADARSAT Geophysical Processor System (RGPS) sea ice trajectories having a temporal and spatial resolution of about 3 days and 10 km, respectively. From the comparison between unfiltered and filtered fields, we estimate that the artificial noise causes an overestimation of about 60% of opening and closing. The artificial noise also has a strong impact on the statistical distribution of the deformation and on the scaling exponents estimated with multifractal analysis. We also show that a similar noise is present in the deformation fields provided in the widely used four-point deformation RGPS data set. These findings may have serious implications for previous studies as the constant overestimation of the opening and closing could lead to a large overestimation of freezing in leads, salt rejection and sea ice ridging.

Variability of sea ice deformation rates in the Arctic and their relationship with basin-scale wind forcing

Abstract. The temporal variability of the moments of probability distribution functions (pdfs) of total sea ice deformation rates in the Arctic is analyzed in the context of the basin-scale wind forcing acting on the ice. The pdfs are estimated for 594 satellite-derived sea ice deformation maps from 11 winter seasons between 1996/1997 and 2007/2008, provided by the RADARSAT Geophysical Processor System. The temporal scale analyzed equals 3 days. The moments of the pdfs, calculated for a range of spatial scales (12.5–900 km), have two dominating components of variability: a seasonal cycle, with deformation rates decreasing throughout winter towards a minimum in March; and a short-term, synoptic variability, strongly correlated with the area-averaged magnitude of the wind stress over the Arctic, estimated based on the NCEP-DOE Reanalysis-2 data (correlation coefficient of 0.71 for the mean deformation rate). Due to scaling properties of the moments, logarithms of higher moments are strongly correlated with the wind stress as well. Exceptions are observed only at small spatial scales, as a result of extreme deformation events, not directly associated with large-scale wind forcing. By repeating the analysis within regions of different sizes and locations, we show that the wind–ice deformation correlation is largest at the basin scale and decreases with decreasing size of the area of study. Finally, we suggest that a positive trend in seasonally averaged correlation between sea ice deformation rates and the wind forcing, present in the analyzed data, may be related to an observed decrease in the multi-year ice area in the Arctic, indicating possibly even stronger correlations in the future.

Spatial scaling of Arctic sea ice deformation

Arctic sea ice deformation arises from spatial gradients in the ice velocity field. This deformation occurs across a wide range of spatial scales, from meters to thousands of kilometers. We analyze 7 years of sea ice deformation data from the RADARSAT Geophysical Processor System (RGPS) covering the western Arctic Ocean. We find that the mean deformation rate is related to the spatial scale over which it is measured according to a power law with exponent ∼−0.2, over a scale range from 10 to 1000 km (e.g., deformation rate doubles for a 30‐fold reduction in scale). Both the exponent and the deformation rate have distinct annual cycles. The exponent becomes more negative in summer as the ice pack weakens and internal stresses are not as readily transmitted over long distances. The deformation rate reaches a minimum in late winter when the ice pack is strongest. The deformation also exhibits considerable localization, in which the largest deformation rates are confined to smaller and smaller areas as the scale of measurement decreases. This supports a model for sea ice based on granular or fracture mechanics. The scaling exponent in the power law relationship tends to be larger in magnitude where the concentration of multiyear ice is low, consistent with a thinner and weaker ice pack. With decreasing multiyear ice in the Arctic and a thinning ice pack, an increase in the deformation rate has already been documented (from buoy data). However, the net effect of several deformation/thickness feedbacks is still uncertain.

Evaluation of high‐resolution sea ice models on the basis of statistical and scaling properties of Arctic sea ice drift and deformation

Sea ice drift and deformation from models are evaluated on the basis of statistical and scaling properties. These properties are derived from two observation data sets: the RADARSAT Geophysical Processor System (RGPS) and buoy trajectories from the International Arctic Buoy Program (IABP). Two simulations obtained with the Louvain‐la‐Neuve Ice Model (LIM) coupled to a high‐resolution ocean model and a simulation obtained with the Los Alamos Sea Ice Model (CICE) were analyzed. Model ice drift compares well with observations in terms of large‐scale velocity field and distributions of velocity fluctuations although a significant bias on the mean ice speed is noted. On the other hand, the statistical properties of ice deformation are not well simulated by the models: (1) The distributions of strain rates are incorrect: RGPS distributions of strain rates are power law tailed, i.e., exhibit “wild randomness,” whereas models distributions remain in the Gaussian attraction basin, i.e., exhibit “mild randomness.” (2) The models are unable to reproduce the spatial and temporal correlations of the deformation fields: In the observations, ice deformation follows spatial and temporal scaling laws that express the heterogeneity and the intermittency of deformation. These relations do not appear in simulated ice deformation. Mean deformation in models is almost scale independent. The statistical properties of ice deformation are a signature of the ice mechanical behavior. The present work therefore suggests that the mechanical framework currently used by models is inappropriate. A different modeling framework based on elastic interactions could improve the representation of the statistical and scaling properties of ice deformation.

Variability of sea ice simulations assessed with RGPS kinematics

Sea ice drift and deformation from coupled ice‐ocean models are compared with high‐resolution ice motion from the RADARSAT Geophysical Processor System (RGPS). In contrast to buoy drift, the density and extent of the RGPS coverage allows a more extensive assessment and understanding of model simulations at spatial scales from ∼10 km to near basin scales and from days to seasonal timescales. This work illustrates the strengths of the RGPS data set as a basis for examining model ice drift and its gradients. As it is not our intent to assess relative performance, we have selected four models with a range of attributes and grid resolution. Model fields are examined in terms of ice drift, export, deformation, deformation‐related ice production, and spatial deformation patterns. Even though the models are capable of reproducing large‐scale drift patterns, variability among model behavior is high. When compared to the RGPS kinematics, the characteristics shared by the models are (1) ice drift along coastal Alaska and Siberia is slower, (2) the skill in explaining the time series of regional divergence of the ice cover is poor, and (3) the deformation‐related volume production is consistently lower. Attribution of some of these features to specific causes is beyond our current scope because of the complex interplay between model processes, parameters, and forcing. The present work suggests that high‐resolution ice drift observations, like those from the RGPS, would be essential for future model developments, improvements, intercomparisons, and especially for evaluation of the small‐scale behavior of models with finer grid spacing.

Intercomparison of passive microwave sea ice concentration retrievals over the high‐concentration Arctic sea ice

Measurements of sea ice concentration from the Special Sensor Microwave Imager (SSM/I) using seven different algorithms are compared to ship observations, sea ice divergence estimates from the Radarsat Geophysical Processor System, and ice and water surface type classification of 59 wide‐swath synthetic aperture radar (SAR) scenes. The analysis is confined to the high‐concentration Arctic sea ice, where the ice cover is near 100%. During winter the results indicate that the variability of the SSM/I concentration estimates is larger than the true variability of ice concentration. Results from a trusted subset of the SAR scenes across the central Arctic allow the separation of the ice concentration uncertainty due to emissivity variations and sensor noise from other error sources during the winter of 2003–2004. Depending on the algorithm, error standard deviations from 2.5 to 5.0% are found with sensor noise between 1.3 and 1.8%. This is in accord with variability estimated from analysis of SSM/I time series. Algorithms, which primarily use 85 GHz information, consistently give the best agreement with both SAR ice concentrations and ship observations. Although the 85 GHz information is more sensitive to atmospheric influences, it was found that the atmospheric contribution is secondary to the influence of the surface emissivity variability. Analysis of the entire SSM/I time series shows that there are significant differences in trend between sea ice extent and area, using different algorithms. This indicates that long‐term trends in surface and atmospheric properties, unrelated to sea ice concentration, influence the computed trends.

Contrasts in sea ice deformation and production in the Arctic seasonal and perennial ice zones

Four years (1997–2000) of RADARSAT Geophysical Processor System (RGPS) data are used to contrast the sea ice deformation and production regionally, and in the seasonal (SIZ) and perennial (PIZ) ice zones. Ice production is of seasonal ice in openings during the winter. Three‐day estimates of these quantities are provided within Lagrangian elements initially 10 km on a side. A distinct seasonal cycle is seen in both zones with these estimates highest in the late fall and with seasonal minimums in the midwinter. Regional divergence over the winter could be up to 30%. Spatially, the highest deformation is seen in the SIZ north of coastal Alaska. Both ice deformation and production are higher in the SIZ: deformation‐related ice production in the SIZ (∼0.5 m) is 1.5–2.3 times that of the PIZ (∼0.3 m): this is connected to ice strength and thickness. Atmospheric forcing and boundary layer structure contribute to only the seasonal and interannual variability. Seasonal ice growth in ice fractures accounts for ∼25–40% of the total ice production of the Arctic Ocean. Uncertainties in these estimates are discussed. By itself, this deformation‐ice production relationship could be considered a negative feedback when thickness is perturbed. However, the overall effect on ice production in the face of increasing seasonal and thinner/weaker ice coverage could be modified by local destabilization of the water column promoting overturning of warmer water due to increased brine rejection; and the upwelling of the pynocline associated with increased occurrence of large shear motion in sea ice. Divergence is shown to be negligibly correlated to cyclonic motion in summer and winter in both ice zones.

Comparison of thin ice thickness distributions derived from RADARSAT Geophysical Processor System and advanced very high resolution radiometer data sets

Thin ice thickness distributions estimated from advanced very high resolution radiometer (AVHRR) and RADARSAT Geophysical Processor System (RGPS) data sets were compared over the Beaufort Sea and the Canada Basin for the period December 1996 to February 1997. The comparisons show a compelling agreement. High correlations were found in cases where thin ice grew in large, wide leads extending several hundred kilometers. At these large scales, estimates from AVHRR images and RGPS showed similar amounts of thin ice in leads. However, when major surface deformation occurred on small scales (100 m to 10 km), the finer spatial resolution (100 m) of RADARSAT images enabled the RGPS algorithm to derive more thin ice than that of AVHRR. Under such conditions the correlation between the two dropped, and a small negative bias (about 1%) was observed in the estimates from AVHRR. This bias, mostly concentrated at the very thin end of the thickness distribution, caused a further deficit in the AVHRR‐derived thin ice growth, roughly 0.2 cm/d. Although the AVHRR and RGPS algorithms treat snowfall differently in the ice thickness calculations, both snow assumptions appear reasonable. However, RGPS may underestimate the thin ice production because of the 3‐day sampling interval. With a better understanding of the sources of uncertainty and improved satellite estimates, a combination of these two satellite data could offer a wider coverage of thin ice thickness observations in the Arctic Basin.

The RADARSAT Geophysical Processor System: Quality of Sea Ice Trajectory and Deformation Estimates

Abstract NASA's RADARSAT Geophysical Processor System (RGPS) uses sequential synthetic aperture radar (SAR) images to track the trajectories of some 30 000 points on the Arctic sea ice for periods of up to 6 months. Much of the Arctic basin is imaged and tracked every 3 days. The result is a highly detailed picture of how the sea ice moves and deforms. The points are initially spaced 10 km apart and are organized into four-cornered cells. The area and the strain rates are calculated for each cell for each new observation of its corners. The accuracy of the RGPS ice tracking, area changes, and deformation estimates is needed to make the dataset useful for analysis, model validation, and data assimilation. Two comparisons are made to assess the accuracy. The first compares the tracking performed at two different facilities (the Jet Propulsion Laboratory in Pasadena, California, and the Alaska SAR Facility in Fairbanks, Alaska), between which the primary difference is the operator intervention. The error sta...

Changes in the modeled ice thickness distribution near the Surface Heat Budget of the Arctic Ocean (SHEBA) drifting ice camp

In the polar oceans the ice thickness distribution controls the exchange of heat between the ocean and the atmosphere and determines the strength of the ice. The Surface Heat Budget of the Arctic Ocean (SHEBA) experiment included a year‐long field program centered on a drifting ice station in the Beaufort and Chukchi Seas in the Arctic Ocean from October 1997 through October 1998. Here we use camp observations and develop methods to assimilate ice thickness and open water observations into a model in order to estimate the evolution of the thickness distribution in the vicinity of the camp. A thermodynamic model is used to simulate the ice growth and melt, and an ice redistribution model is used to simulate the opening and ridging processes. Data assimilation procedures are developed and then used to assimilate observations of the thickness distribution. Assimilated observations include those of the thin end of the distribution determined by aircraft surveys of the surface temperature and helicopter photographic surveys and aircraft microwave estimates of the open water fraction. The deformation of the ice was determined primarily from buoy and RADARSAT Geophysical Processor System (RGPS) measurements of the ice velocity. Because of the substantial convergence and ridging observed in the spring and summer, the estimated mean ice thickness increases by 59%, from 1.53 to 2.44 m, over the year in spite of a net thermodynamic ice loss for most multiyear ice.

Sea‐ice deformation rates from satellite measurements and in a model

The deformation of sea ice is an important element of the Arctic climate system because of its influence on the ice‐thickness distribution and on the rates of ice production and melt. New data obtained from the Radarsat Geophysical Processor System (RGPS) using satellite synthetic aperture radar images of the ice offers an opportunity to compare observations of the ice deformation to estimates obtained from models. The RGPS tracks tens of thousands of points, spaced roughly at 10‐km intervals, for an entire season in a Lagrangian fashion. The deformation is computed from cells formed by the tracked points, typically at 3‐day intervals. We used a coupled ice/ocean model with ice‐thickness and enthalpy distributions that covers the entire Arctic Ocean with a 40‐km grid. Model‐only and model‐with‐data‐assimilation runs were analysed. The data assimilation runs were analysed in order to determine the validity of the comparison techniques and to find the comparisons under the best of circumstances, when many buoy measurements are available for assimilation. This step is necessary because the RGPS and model data differ in spatial and temporal sampling characteristics. The assimilated data included buoy motion and Special Sensor Microwave/Imager (SSM/I)‐derived ice motion. The Pacific half of the Arctic Basin was analysed for a 10‐month period in 1997 and 1998. Comparisons of ice velocity observations to the modelled velocities showed excellent agreement from the model‐with‐data‐assimilation run but poorer agreement for the model‐only run. At a scale of 320 km, the deformation from the data assimilation run was in modest agreement with the observations but where many buoys were available for assimilation the agreement was quite good. Both model runs showed poor agreement during summer. Comparisons of the deformation distribution functions suggest why the agreements were poor even though the velocity agreements were good. Decreasing the ice strength parameter in the model improved the deformation comparisons for the model‐only runs.