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 more » 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. « less
Authors:
; ;
Award ID(s):
1928126
Publication Date:
NSF-PAR ID:
10377946
Journal Name:
The Cryosphere
Volume:
16
Issue:
2
Page Range or eLocation-ID:
533 to 557
ISSN:
1994-0424
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 divergencemore »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.

    « less
  2. 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, coupledmore »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.

    « less
  3. 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 coastmore »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.« less
  4. Abstract. The annual sea ice freeze–thaw cycle plays a crucial role in theArctic atmosphere—ice–ocean system, regulating the seasonal energy balanceof sea ice and the underlying upper-ocean. Previous studies of the sea icefreeze–thaw cycle were often based on limited accessible in situ or easilyavailable remotely sensed observations of the surface. To better understandthe responses of the sea ice to climate change and its coupling to the upperocean, we combine measurements of the ice surface and bottom usingmultisource data to investigate the temporal and spatial variations in thefreeze–thaw cycle of Arctic sea ice. Observations by 69 sea ice mass balancebuoys (IMBs) collected from 2001 to 2018 revealed that the average ice basalmelt onset in the Beaufort Gyre occurred on 23 May (±6 d),approximately 17 d earlier than the surface melt onset. The average icebasal melt onset in the central Arctic Ocean occurred on 17 June (±9 d), which was comparable with the surface melt onset. This difference wasmainly attributed to the distinct seasonal variations of oceanic heatavailable to sea ice melt between the two regions. The overall average onsetof basal ice growth of the pan Arctic Ocean occurred on 14 November (±21 d), lagging approximately 3 months behind the surface freezeonset. This temporal delay was caused by a combinationmore »of cooling the seaice, the ocean mixed layer, and the ocean subsurface layer, as well as thethermal buffering of snow atop the ice. In the Beaufort Gyre region, both(Lagrangian) IMB observations (2001–2018) and (Eulerian) moored upward-looking sonar (ULS) observations (2003–2018) revealed a trend towardsearlier basal melt onset, mainly linked to the earlier warming of thesurface ocean. A trend towards earlier onset of basal ice growth was alsoidentified from the IMB observations (multiyear ice), which we attributed tothe overall reduction of ice thickness. In contrast, a trend towards delayedonset of basal ice growth was identified from the ULS observations, whichwas explained by the fact that the ice cover melted almost entirely by theend of summer in recent years.« less
  5. Abstract
    These data were collected on the CCGS (Canadian Coast Guard) Louis St. Laurent during BGOS (Beaufort Gyre Observing System) research cruises in 2012-2014 and 2016-2021 in the Beaufort Sea area. They are underway pCO2 (Partial pressure of carbon dioxide) data collected using an equilibrator-infrared method (SUPER CO2 system from Sunburst Sensors). Ancillary data for calculation of air-sea CO fluxes include temperature, salinity, atmospheric CO2, wind speed, and gas transfer velocity (calculated from Wanninkhof et al. (2009). Fluxes are not corrected for fractional ice-coverage. The specific goal of the study is to continue to operate an Arctic Observing Network (AON) for the measurement of the partial pressure of CO2 (pCO2), pH, and dissolved O2 (DO) focused on the surface waters of the Arctic Ocean (specifically, the Canada Basin). These data were collected on the CCGS (Canadian Coast Guard) Louis St. Laurent during a BGOS (Beaufort Gyre Observing System) research cruise in the Beaufort Sea area. It is underway pCO2 (partial pressure of carbon dioxide) data collected using an equilibrator-infrared method (SUPER CO2 system from Sunburst Sensors). Ancillary data for calculation of air-sea CO2 fluxes include temperature, salinity, atmospheric CO2.