Estimating Arctic sea-ice freeze-up and break-up from the satellite record : A comparison of different approaches in the Chukchi and Beaufort Seas

The recognized importance of the annual cycle of sea ice in the Arctic to heat budgets, human behavior, and ecosystem functions, requires consistent definitions of such key events in the ice cycle as break-up and freeze-up. An internally consistent and reproducible approach to characterize the timing of these events in the annual sea-ice cycle is described. An algorithm was developed to calculate the start and end dates of freeze-up and break-up and applied to time series of satellite-derived sea-ice concentration from 1979 to 2013. Our approach builds from discussions with sea-ice experts having experience observing and working on the sea ice in the Bering, Chukchi and Beaufort Seas. Applying the algorithm to the 1979–2013 satellite data reveals that freeze-up is delayed by two weeks per decade for the Chukchi coast and one week per decade for the Beaufort coast. For both regions, break-up start is arriving earlier by 5–7 days per decade and break-up end is arriving earlier by 10–12 days per decade. In the Chukchi Sea, “early” break-up is arriving earlier by one month over the 34-year period and alternates with a “late” break-up. The calculated freeze-up and break-up dates provide information helpful to understanding the dynamics of the annual sea-ice cycle and identifying the drivers that modify this cycle. The algorithm presented here, and potential refinements, can help guide future work on changes in the seasonal cycle of sea ice. The sea-ice phenology of freeze-up and break-up that results from our approach is consistent with observations of sea-ice use. It may be applied to advancing our understanding and prediction of the timing of seasonal navigation, availability of ice as a biological habitat, and assessment of numerical models.


Introduction
The annual cycle of Arctic sea ice, from freeze-up to break-up, currently contracts and expands across an area ranging from 7 to 16 x 10 6 km 2 (Serreze et al., 2007).Sea-ice extent affects the fluxes of sensible and latent heat and modifies the regional surface albedo.Based on satellite data from 1979 to 2000, the perennial Arctic sea-ice cover is declining at −7% per decade.Including the years through 2013, the Arctic-wide trend is −14% per decade (Stroeve et al., 2014a,b).
Trends in satellite-derived surface brightness temperature suggest a lengthening of the melt season (Comiso, 2002;Parkinson, 2014).Melt onset is particularly affected by early season atmospheric storms that bring warm air at the same time that incident solar energy increases in May (Bitz et al., 1996;Perovich and Polashenski, 2012).Late spring atmospheric forcing, at the onset of melt, influences low-frequency variability sometimes for decades through nonlinearities in sea-ice volume (thickness and extent), based on energy balance models forced by anomalous melt onset dates (Bitz et al., 1996).Using an energy balance model, Bitz and Roe (2004) found that variations in season length had a "nontrivial" influence on sea-ice thickness even though the sea-ice growth-thickness relationship, where thin ice grows faster thermodynamically than thick ice, dominates the overall ice behavior.transitions from 1999 through 2012 poleward of 60°N at 4.5 km resolution.Similar to passive-microwave approaches, scatterometer data respond to changes in surface snow and ice dielectric properties which may or may not relate to geophysically or operationally relevant parameters.To eliminate unreasonable melt-freeze transitions, Mortin et al. (2014) relied on sea-ice concentration data from the Special Sensor Microwave/ Imager (SSM/I) and Special Sensor Microwave Imager/Sounder (SSMI/S) instruments.Using SSM/I and SSMI/S to constrain estimates of melt-freeze transitions suggests those data have additional utility in studies of the annual cycle of sea ice.
An analysis of abrupt sea-ice changes using ensemble numerical simulations by Holland et al. (2006) found that the dynamic effects of ice transport and ridging had "little direct role" in climate transitions and instead were thermodynamically driven.However, thermodynamic melt onset does not necessarily correspond to the first appearance of open ocean, although changes in microwave emissivity are sometimes erroneously taken for a break-up signal which is a process combining dynamics with thermodynamics.The emissivity spec trum is difficult to relate to specific ice surface type (Baordo and Geer, 2015).While some of the aforementioned studies validated or calibrated their approach based on surface air temperature or satellitederived ice concentration data, we are not aware of research that compares or validates remote-sensing derived data against ground-based observations of freeze-up or break-up.
In this study, we are motivated by three issues relevant to determining the timing of the annual sea-ice cycle.
(1) The character of sea ice is changing, including declining ice (Bjørgo et al., 1997;Parkinson et al., 1999) and thinning ice (Rothrock et al.,1999;Wadhams and Davis, 2000).The seasonal ice cycle is changing based on thermodynamic phase changes (Markus et al., 2009), and the length of the open water season is changing based on the 15% ice concentration threshold (Parkinson, 2014).In this paper we examine the timing and duration of the ice season by defining freeze-up and break-up having ecological significance, relating to human activities, and reflecting what users of the ice and ocean see.This approach draws on long-term observations by sea-ice experts in several Alaskan coastal communities (Eicken et al., 2014).Thus, when "freeze-up" occurs, persistent ice ends the "open water" season.When sea ice "breaks up" it permits seasonal travel on the ocean.Our work is particularly motivated by the need to arrive at an intercomparable, meaningful index of sea-ice phenology.Such phenology should be relevant to ice users, such as humans navigating and hunting on and among sea ice, useful to characterizing the ecological roles of the ice, and internally consistent to be meaningful to large-scale geophysical studies of sea ice.
(2) Analyzing passive-microwave derived ice concentrations for the Arctic from 1979 to 2013, we arrive at time series of freeze-up and break-up that allow us to examine and compare variability and trends.We rely on ice concentration from passive microwave measurements because that data set is well validated and well understood.Of particular interest is that the passive microwave data from 1979 onwards provide a longer context for changes in sea-ice phenology than recent studies (e.g., Brown et al., 2014) where the choice of sensors shifts the focus to the current decade.
(3) Using the passive microwave data from 1979 through 2013, we derive and analyze the timing of freeze-up and break-up in the coastal region of the Chukchi and Beaufort Seas in light of community-based observations of the ice cycle, sea-ice uses, and ice hazards.This region is of interest because it has experienced some of the largest changes in seasonal minimum ice extent (Parkinson, 2014) anywhere in the Arctic, with major increases in solar heating of the surface ice-ocean system (Perovich et al., 2011;Stroeve et al., 2014b) and implications for coastal dynamics (Barnhart et al., 2014).Comparing these data to other measures of freeze-up or break-up provides insight into potential biases and uncertainties of different approaches taken to determine key events in the ice cycle for operational or regulatory purposes.
Figure 1 illustrates the problem at hand by comparing five different measures of sea-ice seasonality and specifically freeze-up for the nearshore region off Barrow, Alaska.The methods address the last day of navigation, direct observations of freeze-up, presence of shore ice from Radarsat SAR, and the remote-sensing approach introduced here.There is a great deal of variation among the different methods.From 1990 to 2013, the dates shown in Figure 1 for "persistent ice", "freeze-up" and "last day of navigation" span more than two months, providing little clarity on the behavior of sea ice at that time, or from year to year.Depending on the specific use of the ice cover or the locally relevant definition, observations by local experts at the same site may not agree.Here we attempt to resolve the problem by developing a coherent framework for analyzing and characterizing the timing of freeze-up and break-up.
3. The sea-ice seasonal cycle of the coastal Chukchi and Beaufort Seas Barry et al. (1979) used LANDSAT imagery from 1973 through 1977 to characterize the extension/modification of fast ice for the central Chukchi coast from November through February (see their Table 1).They observed that "first openings and movement" of the ice were on June 10 and "new ice begins to form in October" with the actual freeze-up date varying by 10-15 days.Barry et al. (1979) noted that the fast ice off Icy Cape and Peard Bay (Chukchi coast) has melt and breakup characteristics similar to that along the Beaufort Sea coast.
Elementa: Science of the Anthropocene • 4: 000124 • doi: 10.12952/journal.elementa.000124Break-up along the Chukchi coast is driven by solar heating and atmospheric and oceanic dynamics (Barry et al., 1979;Petrich et al., 2012).The first open water is along the margin of the shorefast ice.The 20-m isobath, a proxy for grounded ridges stabilizing shorefast ice farthest away from the coast, is generally 10 to 25 km offshore in the Chukchi Sea and about twice that distance off the Beaufort Sea coast (Mahoney et al., 2014).
Ice conditions in the Beaufort Sea are different from conditions in the Chukchi Sea (Barry et al., 1979;Mahoney et al., 2014).The annual cycle for the Beaufort Sea begins with freeze-up over the low-salinity inner shelf and proceeds as follows (Barry et al., 1979;Mahoney et al., 2014) with approximate timing: 1) new ice formation near shore (October); 2) freezing in sheltered bays and lagoons; 3) seaward extension of the ice sheet in late fall and early winter (October-November) and disruption of the ice by late-winter storms (December-February) and pressure from the pack ice; 4) generally stable ice shoreward of the 15-m isobath in late winter and early spring (February-May); 5) over-ice flooding at the river mouths in spring (May); 6) surface melt and pond formation on the ice surface (May-June); 7) cracks or flaws in previously continuous ice sheets ( June); 8) movement of previously immobile ice; and finally 9) nearshore regions becoming largely free of ice ( July-August).
Along the Beaufort Sea coast, a stable and relatively smooth fast ice sheet is present between the coast and the barrier islands from January through May.Because the barrier islands extend along more than half of the 800 km coast, fast ice in the Beaufort Sea is well anchored and protected from pack ice seaward of the islands (Barry et al., 1979).The offshore region was characterized by Kovacs and Mellor (1974) and Reimnitz et al. (1978) as a zone of fast ice, grounded ridges, and seasonal ice pack.
Variations in break-up are due to the strength, thickness, and morphology of the ice and regional weather (Barry et al., 1979).Break-up in the Beaufort Sea, defined by "openings and movement", occurs after about 55-144 thawing degree days (TDDs, defined in °C) according to Barry et al. (1979).Clearing of the fast ice takes approximately 140-220 TDDs (°C).The first stage of break-up is flooding of the lagoons and estuaries in late May or early June (Barry et al., 1979)   Last day of navigation from the Barnett Ice Severity Index (blue), freeze-up as calculated here (red), freeze-up based on direct observations at Barrow by local observer J.C. George (green), date of first persistent ice noted by local observers J.C. George and J.M. Leavitt (black), and first shore ice from Radarsat SAR (white squares).Blue and red lines are linear fits to the data.doi: 10.12952/journal.elementa.000124.f001

Starting or ending date Criteria and algorithm
Start freeze-up Date for which the sea-ice concentration exceeds the average value for August and September ("summer") plus one standard deviation for the first time.A minimum threshold is set to 15%.
End freeze-up Date for which the sea-ice concentration exceeds for the first time the average ice concentration for January and February ("winter"), minus 10%.

Start break-up
Date for which the sea-ice concentration exceeds the threshold calculated as the last day for which the previous two weeks' ice concentration always exceeds the average value for winter ( Jan-Feb), minus two standard deviations.
End break-up Date for which the sea-ice concentration exceeds for the last time the average value for summer (Aug-Sep), plus one standard deviation.A minimum threshold is set to 15%.
doi: 10.12952/journal.elementa.000124.t001 Updating and extending the study of Barry et al. (1979), Mahoney et al. (2014) found that in recent years, landfast ice forms up to two weeks later and breaks up as much as two to four weeks earlier than in the mid-1970s based on Radarsat SAR data from 1996 through 2008.Note, however, that both datasets and definitions of freeze-up and break-up were not entirely consistent between the two studies, reinforcing the arguments made above for the need to develop a coherent analysis framework.
4. Methods: Approach for calculating freeze-up and break-up Indigenous experts primarily in the coastal communities of Wales, Gambell and Barrow, Alaska, with collective knowledge of the Bering, Chukchi and Beaufort Seas have been keeping records of sea-ice, ocean, and weather conditions as part of the Seasonal Ice Zone Observing Network (SIZONet).This communitybased observation component, described in detail by Eicken et al. (2014), has kept daily logs since 2006.Ice experts log weather, ocean and ice conditions relevant to activities and ice use by members of their respect ive communities.Log entries are transcribed and entered into a database for archival, sharing and analysis (Eicken et al., 2014).This work arrived at a representative seasonal ice cycle for 2006-2007 and then placed it in context of interannual variability from 2006 through 2011 of freeze-up and break-up dates.
Observations during the freeze-up season typically focused on ocean access for subsistence hunters, potential hazards due to the presence of slush ice, and the impact of freeze-up on animals and shoreline processes.Excerpts from two entries by Winton Weyapuk, Jr. from Wales during the 2007 freeze-up season illustrate the types of observations made (Apangalook et al., 2013): Skies are mostly cloudy, winds NE at 30,gusting to 35 [mph], temp about 19 F, visibility 7 miles.Light drifting snow.There are narrow strips of light slush in the ocean.Small pieces of pancake ice washed ashore along the tide line.Frozen slush on the beach at the high tide line.Frost smoke not very visible.There are more extensive areas of light slush towards the North.A few clams are washing in.Large flocks of seagulls feeding along the breakers near the shore.[…]. 7 Nov. 2007. Skies are partly cloudy, winds NE at 15 [mph], temp about 23 F, visibility 10 miles.There was occasional moderately heavy snow late in the morning.The slush berm created the day before has been knocked down to half its size or smaller.Thin slush in front of the village and to the North, thicker at the cove.The slush is still soft and not solid enough to walk on.Saw a few small seals just off the beach in front of the village.[…] Winds are higher offshore with white caps visible about 4 to 5 miles out.Small areas of frost smoke out there where the waves are breaking.9 Nov. 2007.
Selected log entries by J. M. Leavitt during the time of freeze-up at Barrow provide insight into the freezeup process, the criteria relevant to local observers in defining the nature of the ice present, and uncertainties associated with determining a freeze-up date.
1 November 2011: NW@16, 26°, Swell [s] are here at 5′. Slush ice goes about ½ mi out.Lots of seagulls along the beach, snowing 1.5 mi visibility [,]  As is evident from these observations, slush ice formed during a snow storm early in the season but did not persist, with open water observed until the first persistent ice set up on November 13, determined to be the date of freeze-up.Here, freeze-up was a result of a combination of advection of ice formed offshore and in situ freezing.For local ice use, the strength of the ice is particularly relevant, with a key criterion being sufficient thickness to walk on the ice.
Observations and analysis of freeze-up and break-up as described above were used in developing a statistical algorithm to be applied to the SSM/I and SSMI/S data to calculate freeze-up and break-up dates.The algorithm is described here and applied to the satellite record from October 28, 1978, through September 30, 2013.The algorithm was derived by comparing three different approaches to determine freeze-up and Elementa: Science of the Anthropocene • 4: 000124 • doi: 10.12952/journal.elementa.000124break-up dates based on time series of satellite ice concentration obtained from passive microwave data (M.-L.Kapsch and H. Eicken, unpublished) and selecting the algorithm that showed the smallest deviations between ground-based observations and remote sensing data.We focus in this paper on the northeastern Chukchi Sea and southern Beaufort Sea coasts near Barrow, Alaska.
The SSMR, SSM/I and SSMI/S sea-ice concentration data derived using the NASA Team algorithm were extracted from several platform frequencies and polarizations (Cavalieri et al., 1996).The SMMR data span 1978 to 1987 and the follow-up SSM/I onboard Defense Meteorological Satellite Program (DMSP) satellites F8, F11, and F13 cover 1987 through 2007.Beginning in January 2007 SSMI/S data collected aboard DMSP-F17 have become available.The SMMR data were collected every other day with gaps in August 1982 and 1984.The SSM/I and SSMI/S data were collected daily with a data gap from December 3, 1987 to January 13, 1988.The NASA Team algorithm uses the multi-frequency brightness temperature data to calculate total ice concentration, mapped to a 25-km polar stereographic grid, from which total ice area and total ice extent are derived.These data, readily available through the National Snow and Ice Data Center, and based on GHz frequency microwave (millimeter to centimeter wavelengths) satellite radiance measure ments, are among the best known, well validated, long-term, and pan-Arctic seaice records.
The sea-ice concentration grid points in this study (Figure 2) are located in the coastal Chukchi and Beaufort Seas.Grid points in proximity to land may be subject to errors from land-to-ocean spillover due to the coarseness of the satellite antennae pattern.Strong brightness temperature contrasts that are smeared across grid points can lead to falsely high estimates of sea-ice concentration.Spillover varies in time and is particularly important in late summer when there is open water near the coast.
To address this problem, the NASA Team algorithm applies the method detailed in Cavalieri et al. (1999) where "excess" concentration, the amount above a minimum threshold, is subtracted from concentration at coastal grid points.Monthly minimum concentration thresholds were established from 1984 and 1992 data.The spillover correction reduces values to the minimum threshold and was applied to daily sea-ice concentrations at grid points in the vicinity of land and open water.Figure 2 in Cavalieri et al. (1999) shows the 7 x 7 array defining the coastal grid point neighborhood, and their Figure 3 displays the results of applying the spillover correction.The correction scheme "reduces substantially the land spillover" (Cavalieri et al., 1999).
To evaluate the data used here for possible spillover, we visually inspected time series of summer seaice concentration in both the Chukchi and Beaufort regions.Ice concentration decreases in the offshore direction with grid points closest to shore generally having sea-ice concentrations that vary in time between 0% and ∼ 20%.For this study we calculated the average ice concentration from August through September from a set of grid points to establish a summer concentration threshold used in freeze-up.The two-month average is expected to reduce any remaining spillover effects that are not constant in time.
NSIDC overview documents (nsidc.org)state that these data are within ± 15% of the actual summer concentration in the presence of summer melt ponds and within ± 5% of the actual winter concentration.In winter when the ice concentration is high, we expect the break-up computations to be robust.In the following discussions the decreasing accuracy of the data as the proportion of thin ice increases, and potential issues with summer data when melt is present or when there is significant new sea ice forming, must be kept in mind.
Different approaches to defining freeze-up and break-up using time series of ice concentration derived from passive microwave data were compared (M.-L.Kapsch and H. Eicken, unpublished).The definitions and algorithm used in this paper and summarized in Table 1 provide the closest approximation to freeze-up and break-up dates as defined by local observers for the local area of coastal ocean near the Alaska coastal communities of Gambell, Wales and Barrow (Eicken et al., 2014).The 25-km NISDC grid points are marked by light tan dots.Green dots in the west denote "Chukchi", in the east "Beaufort", and red squares mark the subset for "Barrow".The Chukchi and Beaufort sets are each 24 points.Barrow is marked by the blue dot.Binary data were downloaded from the NSIDC website, linearly interpolated to daily values of seaice concentration centered on noon, and spatially averaged as described below.The time series data were Hann ing smoothed (weights of 0.25 to 0.5 to 0.25) to reduce short-term events lasting less than ∼ 3 days.This shifts some freeze-up and break-up dates toward climatological values.
For convenience we define the sea-ice year from September 1 to August 31 of the following year.Calculating freeze-up start from the criteria listed in Table 1 requires a summer average (August-September) sea-ice concentration.A threshold is set each year to the larger of the summer average sea-ice concentration plus one standard deviation of the August through September concentrations, or 15%.Once a threshold is set, it is compared with the concentrations from the daily time series beginning September 1. Freeze-up start is the date when the ice concentration first exceeds this threshold.To calculate freeze-up end, a threshold of the average winter ( January-February) concentration minus 10% is set.Freeze-up end is the first date when the daily time series concentration exceeds the threshold.
To calculate break-up start, a threshold is set to the winter ( January-February) average concentration minus two standard deviations of the January through February concentrations.Daily concentrations within a two-week window beginning January 1 are compared to the threshold value.If all daily concentrations exceed the threshold, the next day is evaluated.Break-up start is the last day in the two-week window when all values exceed the threshold.Break-up end sets a threshold to the average summer (August-September) concentration plus one standard deviation of the August through September concentrations, with a minimum threshold of 15%.Break-up end is the last day to exceed the threshold.Note that the algorithm requires the summer sea-ice concentration to calculate a break-up end date that generally occurs in the spring.When the August-September average sea-ice concentration exceeds 25%, the dates for freeze-up start and breakup end are not calculated.When the summer sea-ice concentration exceeds 40%, freeze-up end and breakup start dates are not calculated.
Daily time series of sea-ice concentration were created for the three sets of grid points shown in Figure 2 and are referred to in the following discussions as Chukchi, Beaufort, and Barrow.A daily climatology was computed for each set as the 34-year (1978 to 2013) mean at each Julian Day.Dates of freeze-up and break-up were calculated using the algorithm described above and in Table 1 from the daily data and from the climatology sea-ice concentration time series.

Results
The algorithm described here determines dates that capture sea-ice behavior important to indigenous experts and ice users who depend on detailed knowledge of the seasonal ice cycle.The algorithm uses thresholds that vary in time.Each year, the "summer" and "winter" ice conditions are used to determine when freeze-up and break-up occur, consistent with ice observers describing freeze-up and break-up as a relative annual measure.This approach is different from applying a fixed concentration threshold, such as 15%, to determine ice season length (Parkinson, 2014).
A schematic of freeze-up and break-up seasonality using this approach (Figure 3) shows the mean start and end dates of freeze-up and break-up calculated from the algorithm described above for the entire set of coastal grid points (Figure 2  The computed start date for freeze-up and break-up are both marked by solid green lines.The end dates for freeze-up and breakup are marked by solid red lines.Seasonal durations are shown by dark red arrows from "breakup start to freeze-up end" and "break-up end to freeze-up start".Dates from Markus et al. (2009) for early melt onset and early freeze onset from their Chukchi/ Beaufort Seas region are shown by the dotted magenta lines, and melt onset and freeze onset by dotted blue lines.Their "inner" and "outer" seasonal durations are marked by blue arrows.doi: 10.12952/journal.elementa.000124.f003melt season lengths (the "inner" and "outer" seasons as described previously) for the same period in the "Chukchi/ Beaufort seas" region examined by Markus et al. (2009).Their "Chukchi/Beaufort seas" region extends from Bering Strait to ∼ 81°N and from Wrangel Island to Banks Island (see their Figure 8).The lengths of the two seasons from our study differ more than the two seasons from Markus et al. (2009).Our approach separates well the annual cycle and marks the start and end of rapid sea-ice transitions in fall and spring.
Calculated freeze-up and break-up dates and linear trends from this study and Markus et al. (2009) are listed in Table 2. Compared to Markus et al. (2009), freeze-up start is 26 days later than early freeze onset (October 16 compared to September 19), freeze-up end is 46 days after freeze onset, break-up start is the same as early melt onset, and break-up end is 51 days after melt onset.While the mean dates are different, as expected, the trends from this study and from Markus et al. are all identical in sign.The freezing season trends have similar magnitude while the melt season trends from this study are more than twice that found by Markus et al. (2009).Figure 5 shows daily ice concentration for 1982 through 1986 over the (repeated) record-length climatology.September ice concentration exceeds climatology for summers of 1983 and 1985 that may be due to pack ice incursions that are well known to occur on the Chukchi coast.The most recent year for which local Iñupiaq ice experts reported such incursions was the summer of 2006, described by Richard Glenn from Barrow as similar to the types of ice conditions that were observed in the 1970s and 1980s, with ice linger ing through summer into the freeze-up period (R. Glenn, personal communication, 2006).
When the summer ice concentration remains high, as shown in Figure 5, a calculation of break-up end is undefined in the current algorithm.This occurred in seven of the 35 years evaluated.In these years, ice remained along the coast for extended periods in summer.In some cases the ice concentration increased enough to significantly change the threshold but lasted for only a short period.Ice incursions and similar events described by local experts are detailed in Johnson et al. (2014).
The algorithm successfully calculates freeze-up and break-up start and end dates for every year after 1994.Figure 6 shows ice concentration for 2007 through 2012 when there were no summer ice incursions and the summer ice concentration fell to minimum levels each summer.The length of the summer ice-free season is clearly longer from 2007 onward than in 1983 and 1985 as shown in Figure 5.
The dates for freeze-up and break-up start and end for the Chukchi Sea region are plotted in Figure 7.There is substantial inter-annual variability with standard deviations of the mean dates at around three weeks.The correlation coefficient for freeze-up start is 0.6, freeze-up end is 0.4, break-up start is −0.4 and break-up end is −0.6.The linear fit shows freeze-up arriving later and break-up arriving sooner (values in Table 4).We identified outliers and years when freeze-up or break-up could not be calculated using the definitions from Table 1.These are listed in Table 5 to show that freeze-up or break-up could not be calculated for seven of the 34 record-length years, all of which are before 1995.Most of these years had significant summer ice present.

Discussion
The algorithm that was applied to ice concentration time series arose following discussions with indigenous sea-ice experts who regularly observe and use sea ice.Based on those discussions, we constructed an algorithm to define the seasonal ice cycle by calculating dates of freeze-up and break-up start and end.
A comparison between freeze-up start and end dates to local observations (Figure 8) provides an indication of the uncertainty associated with the passive microwave data as to the onset and termination of freeze-up.This uncertainty is reflected for specific years in anomalies associated with a substantial discrepancy between passive microwave data and local observations or between local observations and webcam/photo records.
Most observations straddle the identity line where observations and the calculated dates agree, but data for 2006 and 2007 deviate significantly.For 2006, this deviation is explained by the delay in the formation of persistent, stable coastal and shorefast ice during a year with freeze-up dominated by advection into and export out of the area along with late local ice formation.As reported by J. M. Leavitt and A. Brower, Sr., ice kept forming and melting or being advected offshore with open water persisting over large areas from the beach outward through mid-December.On December 15 the first persistent young ice took hold alongshore; this ice was deemed safe to walk on the following day, but with large areas of open water remaining (Leavitt, in Apangalook et al., 2013).In 2007, a different mode of freeze-up occurred at Barrow, with young ice forming and persisting nearshore from November onwards.As confirmed by local experts, such ice can form in shallows or lagoons and persist alongshore despite the absence of extensive pack ice offshore.
The uncertainties evident in Figure 8 in defining a single freeze-up date reflect the recognition by local experts and other ice users (such as the maritime industry) that with a combination of in situ ice formation and advective processes at a particular location, freeze-up is an extended period whose bounds depend on the specific ice use or definition of freeze-up.The same holds for break-up in the spring.The definitions for both events in the seasonal cycle, and their start and end dates, reflect this fact and may hence be of value from an operational or ice use perspective.
Elementa: Science of the Anthropocene • 4: 000124 • doi: 10.12952/journal.elementa.000124From 34 years of satellite data the calculated trends of delayed freezing and earlier melting are consistent in sign and general magnitude with those of Markus et al. (2009), Mortin et al. (2014), andBrown et al. (2014).We can use the calculated long-term trends in the Chukchi Sea to extrapolate freeze-up and break-up dates for a future year.While the choice of a future year is somewhat arbitrary, we first choose 2030 as an example, in part because Overland and Wang (2013) label 2030 as a midpoint among time horizons for a "nearly sea-ice free summer" Arctic.The linear projection from data presented here is for freeze-up to start on December 6 and end on December 13, 2030.Break-up start projects to April 19 and break-up end projects to June 10.These results indicate that by around 2030, the open water season between break-up end and freeze-up start, is from early June until early December, a range nearly twice the 1978-2013 mean range of 97 days (see Table 4).
With freeze-up arriving later and break-up earlier, we can also determine from the trends the year when they coincide.At such a time we assume there would be no "winter" season because the ice would break-up as soon as it freezes-up.From this study's linear projections, freeze-up ends and break-up begins around the same time in 2100.In February-March 2100, there is no duration for the frozen season such that a "solid" winter ice cover might not form.While the trends are clear, this analysis is a local one focusing on the Chukchi and Beaufort coasts.We know of no reason to support the assumption that a linear trend would be sustained even as the ice cover undergoes major change.Rather, non-linear adjustments are likely, hence the need for long-term tracking of ice conditions.
The trends of later freezing and earlier melting raise questions about the role of sea-ice advection into the Chukchi Sea because advection of ice into a region affects the calculation of the date of break-up.There are seven years where break-up could not be calculated (Table 5), and they all fall before 1995, consistent with a possible retreat of the ice pack.We noted above the possible interpretation that by 2100 the sea ice may retreat poleward of Barrow.To examine this possibility, we analyzed where the 15% sea-ice concentration (also used by Mortin et al., 2014 to help constrain melt-freeze transition dates) fell on longitude 169°W (longitude of Bering Strait) from 1979 through 2013.The 15% isoline was found to retreat poleward at this longitude at a rate of ∼ 250 km per decade.If this trend persists, then the 15% isoline retreats to the pole by 2100.In this case, sea-ice incursions from the central Arctic pack may be the dominant source of summer season ice in the Chukchi Sea.
Inspection of the break-up start dates (Figure 7B) shows eleven years well below the regression line.To focus on these "early" break-up years, a Hovmöller diagram was constructed of sea-ice concentration anomalies (annual cycle removed) for the Chukchi points shown in Figure 2.This diagram (Figure 9) shows an oscillation in break-up start (green line) that is strongly bi-annual from 1982 through 2000 where early break-up years alternate with later break-up years that are generally later by about a month.Selecting only the "early break-up start" dates gives a linear correlation coefficient of -0.88, indicating a strong, linear trend for earlier "early" break-up start (black trend line in Figure 9).Since 1979, early break-up has become earlier by more than one month.
and at the mouths of the Colville, Sagavanirktok, and Kuparuk Rivers.(River deltas account for 16% of the total Beaufort Sea coast.)When rivers flood, nearby bottom fast ice floats up and cracks develop.Nearshore open water spreads laterally and seaward from mid-June through early July as melt ponds form and the ice sheet thins.Elementa: Science of the Anthropocene • 4: 000124 • doi: 10.12952/journal.elementa.000124 Figure 1Methods of measuring freeze-up.
doi: 10.12952/journal.elementa.000124.f002 ) for the period of 1979-2007.Also shown are the transition dates and Elementa: Science of the Anthropocene • 4: 000124 • doi: 10.12952/journal.elementa.000124 Figure 3Seasonal timing and the climatological sea-ice cycle.Climatology of the 1979-2007 seaice concentration for the combined Chukchi/Beaufort coastal region is shown by the solid black line from May through December.The computed start date for freeze-up and break-up are both marked by solid green lines.The end dates for freeze-up and breakup are marked by solid red lines.Seasonal durations are shown by dark red arrows from "breakup start to freeze-up end" and "break-up end to freeze-up start".Dates fromMarkus et al. (2009) for early melt onset and early freeze onset from their Chukchi/ Beaufort Seas region are shown by the dotted magenta lines, and melt onset and freeze onset by dotted blue lines.Their "inner" and "outer" seasonal durations are marked by blue arrows.
doi: 10.12952/journal.elementa.000124.t004Mean dates of freeze-up and break-up for the Barrow daily data and for climatology are shown in Figure 4 for two seasonal cycles, 1997-1998 and 2005-2006.These years were selected to illustrate freeze-up and break-up differences and the algorithm's sensitivity to variations in the ice concentration time series.Year 1997-1998 had highly variable daily ice concentration in winter fluctuating between 50 and 100%.In contrast, year 2005-2006 had persistent 100% ice concentration for periods up to a month.There is close agreement between the daily data and climatology freeze-up end dates for both 1997 and 2005, and similarly close dates for break-up end for 1998 and 2006.Freeze-up start from daily data is later than climatology for both 1997 and 2005.Break-up start is later than climatology for both 1998 and 2006.Elementa: Science of the Anthropocene • 4: 000124 • doi: 10.12952/journal.elementa.000124 Figure 4 Sea-ice concentration record off Barrow, Alaska.Calculated dates of freeze-up start and end (solid blue squares) and break-up start and end (solid red squares) are marked over the climatological (1978-2013) seaice concentration (thick black line) for Barrow.Calculated freeze-up and break-up dates are marked by open squares for the daily sea-ice concentration (thin blue line) for (A) 1997-1998 and (B) 2005-2006.doi: 10.12952/journal.elementa.000124.f004

Figure 5 Freeze
Figure 5 Freeze-up and break-up with two summer ice events.Calculated dates of freeze-up (open blue squares) and breakup (open red squares) from 1982 through 1986 over daily sea-ice concentration (blue line) and the 34-year climatology (black line).Dates for end of break-up could not be calculated for 1983 and 1985 when significant ice (red arrows) remained in the region or was advected into the region during summer.doi: 10.12952/journal.elementa.000124.f005

Figure 6 Freeze
Figure 6 Freeze-up and break-up from 2007 to 2012.Calculated dates of freeze-up (open blue squares) and break-up (open red squares) over daily sea-ice concentration (blue line) and the repeated 34-year climatology (black line).During summer in these years the sea-ice concentration is well below climatology and near zero.Extended winter periods of 100% ice concentration are absent.doi: 10.12952/journal.elementa.000124.f006

Figure 7 Chukchi
Figure 7 Chukchi Sea break-up and freezeup dates.Chukchi Sea (A) freeze-up start (green) and end (red) dates and (B) break-up start (green) and end (red) dates with linear trends.Freeze-up start delays by 14 days per decade and freeze-up end delays by 11 days per decade.The mean duration of freeze-up is 22 days.Break-up start begins earlier by 9 days per decade and breakup end begins earlier by 12 days per decade.The mean duration of break-up is 58 days.The x-axis is from 1978 to 2014.The y-axis limits are dates of the month.The standard deviations of the mean dates (horizontal lines) are shown at right (dotted vertical lines).doi: 10.12952/journal.elementa.000124.f007 Elementa: Science of the Anthropocene • 4: 000124 • doi: 10.12952/journal.elementa.000124 Figure 8 Comparison of freeze-up start and end dates.Comparison of freeze-up period derived from passive microwave satellite data with that derived from local observations at Barrow by observer Joe M. Leavitt and from webcam/ photographic images.Symbols mark mid-point of freeze-up interval, vertical lines indicate range from start to end of freezeup.The 2006 and 2007 outliers are discussed in the text.The dotted line marks the one-to-fit where local observations equal day of freeze-up from passive microwave data.doi: 10.12952/journal.elementa.000124.f008 doi: 10.12952/journal.elementa.000124.f009

Table 2 . Comparison of sea-ice cycle phenologies for the Chukchi/Beaufort Seas, 1979-2007
a Standard deviation in parentheses.For this study, n = 24 (freeze-up) and n = 20 (break-up).b Asterisk indicates significant trend at the 95% level; two asterisks, at the 99% level.c EFO indicates early freeze onset; FO, continuous freeze onset; EMO, early melt onset; MO, continuous melt onset.

Table 4 . Mean freeze-up and break-up dates and trends for the Chukchi and Beaufort Sea coastal regions, 1979-2013
bN values in parentheses.

Table 5 . Freeze-up and break-up outlier dates a Year Break-up Freeze-up Comments Start End Start End Early Late Early Late Early Late Early Late
a NC indicates date could not be calculated; x indicates feature described in comments column.doi: 10.12952/journal.elementa.000124.t005