skip to main content


Title: A new state-dependent parameterization for the free drift of sea ice
Abstract. Free-drift estimates of sea ice motion are necessary to produce a seamless observational record combining buoy and satellite-derived sea ice motionvectors. We develop a new parameterization for the free drift of sea ice based on wind forcing, wind turning angle, sea ice state variables(thickness and concentration), and estimates of the ocean currents. Given the fact that the spatial distribution of the wind–ice–ocean transfercoefficient has a similar structure to that of the spatial distribution of sea ice thickness, we take the standard free-drift equation and introducea wind–ice–ocean transfer coefficient that scales linearly with ice thickness. Results show a mean bias error of −0.5 cm s−1(low-speed bias) and a root-mean-square error of 5.1 cm s−1, considering daily buoy drift data as truth. This represents a 35 %reduction of the error on drift speed compared to the free-drift estimates used in the Polar Pathfinder dataset (Tschudi et al., 2019b). Thethickness-dependent transfer coefficient provides an improved seasonality and long-term trend of the sea ice drift speed, with a minimum (maximum)drift speed in May (October), compared to July (January) for the constant transfer coefficient parameterizations which simply follow the peak inmean surface wind stresses. Over the 1979–2019 period, the trend in sea ice drift in this new model is +0.45 cm s−1 per decadecompared with +0.39 cm s−1 per decade from the buoy observations, whereas there is essentially no trend in a free-driftparameterization with a constant transfer coefficient (−0.09 cm s−1 per decade) or the Polar Pathfinder free-drift input data(−0.01 cm s−1 per decade). The optimal wind turning angle obtained from a least-squares fitting is 25∘, resulting in a meanerror and a root-mean-square error of +3 and 42∘ on the direction of the drift, respectively. The ocean current estimates obtained from theminimization procedure resolve key large-scale features such as the Beaufort Gyre and Transpolar Drift Stream and are in good agreement with oceanstate estimates from the ECCO, GLORYS, and PIOMAS ice–ocean reanalyses, as well as geostrophic currents from dynamical ocean topography, with aroot-mean-square difference of 2.4, 2.9, 2.6, and 3.8 cm s−1, respectively. Finally, a repeat of the analysis on two sub-sections of thetime series (pre- and post-2000) clearly shows the acceleration of the Beaufort Gyre (particularly along the Alaskan coastline) and an expansion ofthe gyre in the post-2000s, concurrent with a thinning of the sea ice cover and the observed acceleration of the ice drift speed and oceancurrents. This new dataset is publicly available for complementing merged observation-based sea ice drift datasets that include satellite and buoydrift records.  more » « less
Award ID(s):
1928126
NSF-PAR ID:
10377946
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
The Cryosphere
Volume:
16
Issue:
2
ISSN:
1994-0424
Page Range / eLocation ID:
533 to 557
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Seasonal predictability of the minimum sea ice extent (SIE) in the Laptev Sea is investigated using winter coastal divergence as a predictor. From February to May, the new ice forming in wind-driven coastal polynyas grows to a thickness approximately equal to the climatological thickness loss due to summer thermodynamic processes. Estimating the area of sea ice that is preconditioned to melt enables seasonal predictability of the minimum SIE. Wintertime ice motion is quantified by seeding passive tracers along the coastlines and advecting them with the Lagrangian Ice Tracking System (LITS) forced with sea ice drifts from the Polar Pathfinder dataset for years 1992–2016. LITS-derived landfast ice estimates are comparable to those of the Russian Arctic and Antarctic Research Institute ice charts. Time series of the minimum SIE and coastal divergence show trends of −24.2% and +31.3% per decade, respectively. Statistically significant correlation ( r = −0.63) between anomalies of coastal divergence and the following September SIE occurs for coastal divergence integrated from February to the beginning of May. Using the coastal divergence anomaly to predict the minimum SIE departure from the trend improves the explained variance by 21% compared to hindcasts based on persistence of the linear trend. Coastal divergence anomalies correlate with the winter mean Arctic Oscillation index ( r = 0.69). LITS-derived areas of coastal divergence tend to underestimate the total area covered by thin ice in the CryoSat-2/SMOS (Soil Moisture and Ocean Salinity) thickness dataset, as suggested by a thermodynamic sea ice growth model.

     
    more » « less
  2. Abstract

    Data collected by two buoy arrays that operated during the ice seasons of 2014/2015 and 2016/2017 were used to characterize annual cycles of ice motion and deformation in the western Arctic Ocean. An anomalously strong and weak Beaufort Gyre in 2014/2015 and 2016/2017 induced generally anticyclonic and cyclonic sea ice drift during 2014/2015 and 2016/2017, respectively. Cyclonic ice motion resulted in higher contributions of ice divergence to total ice deformation in 2016/2017 than in 2014/2015. In 2014, the autumn ice concentration and multiyear ice coverage were higher than in 2016; consequently, the response of ice motion to wind forcing was weak, and less ice deformation was observed in autumn 2014. During the autumn‐winter transition, the ice‐wind speed ratio, ice deformation rate and its spatial and temporal scaling exponents, and localization of ice deformation decreased markedly in both 2014/2015 and 2016/2017 as a result of freeze‐up and consolidation of ice floes. Such dynamic behavior was maintained through to spring with the further thickening of ice cover. Ice deformation increased due to weakened ice strength as summer approached. The amplitude of the annual cycle of ice deformation rate in the western Arctic Ocean in 2014/2015 and especially in 2016/2017 was larger than that observed during the Surface Heat Budget of the Arctic Ocean (SHEBA) program in 1997/1998. We attribute this phenomenon to ice loss during the recent summers, especially of thick multiyear ice.

     
    more » « less
  3. Abstract

    A theory for the mean ice thickness and the Transpolar Drift in the Arctic Ocean is developed. Asymptotic expansions of the ice momentum and thickness equations are used to derive analytic expressions for the leading-order ice thickness and velocity fields subject to wind stress forcing and heat loss to the atmosphere. The theory is most appropriate for the eastern and central Arctic, but not for the region of the Beaufort Gyre subject to anticyclonic wind stress curl. The scale analysis reveals two distinct regimes: a thin ice regime in the eastern Arctic and a thick ice regime in the western Arctic. In the eastern Arctic, the ice drift is controlled by a balance between wind and ocean drag, while the ice thickness is controlled by heat loss to the atmosphere. In contrast, in the western Arctic, the ice thickness is determined by a balance between wind and internal ice stress, while the drift is indirectly controlled by heat loss to the atmosphere. The southward flow toward Fram Strait is forced by the across-wind gradient in ice thickness. The basic predictions for ice thickness, heat loss, ice volume, and ice export from the theory compare well with an idealized, coupled ocean–ice numerical model over a wide range of parameter space. The theory indicates that increasing atmospheric temperatures or wind speed result in a decrease in maximum ice thickness and ice volume. Increasing temperatures also result in a decrease in heat loss to the atmosphere and ice export through Fram Strait, while increasing winds drive increased heat loss and ice export.

     
    more » « less
  4. Abstract

    To track sea ice motion, four ice‐tethered buoys were deployed at 84.6°N and 144.3°W, 87.3°N and 172.3°W, 81.1°N and 157.4°W, and 82.8°N and 166.5°W in summers of 2008, 2010, 2014, and 2016, respectively. In addition, the remote sensed ice motion product provided by National Snow and Ice Data Center was used to reconstruct backward and forward ice drifting trajectories from the buoy deployment sites during 1979–2016. Sea ice in the central Arctic Ocean in late summer is trending to have travelled from lower latitudes, and to be advected to the region more involved in the Transpolar Drift Stream (TDS) during 1979–2016. The strengthened TDS has played a crucial role in Arctic sea ice loss from a dynamic perspective. The trajectory of ice is found to be significantly related to atmosphere circulation indices. The Central Arctic Index (CAI), defined as the difference in sea level pressure between 84°N, 90°W and 84°N, 90°E, can explain 34–40% of the meridional displacement along the backward trajectories, and it can explain 27–40% of the zonal displacement along the forward trajectories. The winter Beaufort High (BH) anomaly can explain 18–27% of the zonal displacement. Under high positive CAI values or high negative winter BH anomalies, floes from the central Arctic tended to be advected out of the Arctic Ocean through Fram Strait or other marginal gateways. Conversely, under high negative CAI values or high positive winter BH anomalies, ice tended to become trapped within a region close to the North Pole or it drifted into the Beaufort Gyre region. The long‐term trend and spatial change in Arctic surface air temperature were more remarkable during the freezing season than the melt season because most energy from the lower troposphere is used to melt sea ice and warm the upper ocean during summer.

     
    more » « less
  5. null (Ed.)
    Abstract. Ocean–sea-ice coupled models constrained by various observations provide different ice thickness estimates in the Antarctic. We evaluatecontemporary monthly ice thickness from four reanalyses in the Weddell Sea: the German contribution of the project Estimating the Circulation and Climate ofthe Ocean Version 2 (GECCO2), the Southern Ocean State Estimate (SOSE), the Ensemble Kalman Filter system based on the Nucleus for European Modelling of the Ocean (NEMO-EnKF) and the Global Ice–Ocean Modeling and Assimilation System (GIOMAS). The evaluation is performed againstreference satellite and in situ observations from ICESat-1, Envisat, upward-looking sonars and visual ship-based sea-ice observations. Compared withICESat-1, NEMO-EnKF has the highest correlation coefficient (CC) of 0.54 and lowest root mean square error (RMSE) of 0.44 m. Compared within situ observations, SOSE has the highest CC of 0.77 and lowest RMSE of 0.72 m. All reanalyses underestimate ice thickness near the coast ofthe western Weddell Sea with respect to ICESat-1 and in situ observations even though these observational estimates may be biased low. GECCO2 andNEMO-EnKF reproduce the seasonal variation in first-year ice thickness reasonably well in the eastern Weddell Sea. In contrast, GIOMAS ice thicknessperforms best in the central Weddell Sea, while SOSE ice thickness agrees most with the observations from the southern coast of the Weddell Sea. Inaddition, only NEMO-EnKF can reproduce the seasonal evolution of the large-scale spatial distribution of ice thickness, characterized by the thickice shifting from the southwestern and western Weddell Sea in summer to the western and northwestern Weddell Sea in spring. We infer that the thickice distribution is correlated with its better simulation of northward ice motion in the western Weddell Sea. These results demonstrate thepossibilities and limitations of using current sea-ice reanalysis for understanding the recent variability of sea-ice volume in the Antarctic. 
    more » « less