1. Introduction
Wind-driven flaw polynyas are nonlinear-shaped regions of open water and thin ice within a closed ice cover, formed by offshore winds advecting the pack ice away from fast-ice edges (Reference Smith, Muench and PeaseSmith and others, 1990). In flaw polynyas, large amounts of ice are generated and salt is expelled into the underlying waters (Reference Renfrew, King and MarkusRenfrew and others, 2002; Reference Morales Maqueda, Willmott and BiggsMorales Maqueda and others, 2004). As Arctic sea-ice volume shrinks rapidly (Reference Haas, Pfaffling, Rabenstein, Etienne and RigorHaas and others, 2008; Reference Kwok, Cunningham, Wensnahan, Rigor, Zwally and YiKwok and others, 2009), the monitoring of ice production within Arctic polynyas with satellites and models is a crucial objective of ongoing polar research.
Passive microwave satellites allow for a direct deduction of polynya area and, together with atmospheric data, a retrieval of thin-ice thickness and ice production. With their global coverage and their ability to penetrate cloud cover, passive microwave sensors can be used for long-term and large-scale monitoring of polynya dynamics and associated ice formation (Reference Tamura, Ohshima, Markus, Cavalieri, Nihashi and HirasawaTamura and others, 2007; Reference KernKern, 2009; Reference Willmes, Adams, Schröder and HeinemannWillmes and others, in press). Unfortunately, the accuracy of the retrieval method depends on regional particularities of polynya size and shape. the coarse spatial resolution of passive microwave satellite data, if applied to narrow polynyas, generates errors through mixed signals at the fast-and pack-ice edges (Willmes and others, 2010), and results in an underestimation of ice production.
Thus, monitoring of ice production in narrow leads and polynyas remains difficult with satellites alone. This gap may be filled by flux models capable of simulating polynya evolution and ice fluxes. Flux balance models were first formulated by Reference PeasePease (1987), embracing the idea of Reference LebedevLebedev (1968) that wind-generated polynyas attain a maximum size resulting from a balance between the flux of frazil ice produced in the open-water area and the wind-driven offshore divergence of consolidated new ice. If fluxes are not in balance, the polynya area is evolving.
Questions still remain as to the consistency and accuracy of flux model assessments. Reference Markus and BurnsMarkus and Burns (1995) found that the area of a polynya retrieved from passive microwave data agreed reasonably well with area estimates made by a one-dimensional (1-D) polynya flux model. These good agreements were confirmed by Reference Haarpaintner, Haugan and GascardHaarpaintner and others (2001) and Reference Skogseth, Haugan and HaarpaintnerSkogseth and others (2004), who compared flux simulations in the Storfjorden (Svalbard) polynya with ice-area and ice-drift information obtained from synthetic aperture radar (SAR) images. the estimation of ice and salt fluxes in polynyas requires, besides correct representation of the polynya extent, an accurate determination of ice growth in both open-water and thin-ice zones. However, a direct comparison of ice-production estimates made by flux models with satellite-based approximations is still missing.
The aim of this paper is therefore to test whether flux models provide realistic ice-production estimates and consequently are capable of filling the gap in narrow-lead and polynya ice-production monitoring.
Below we apply a two-dimensional (2-D) flux model developed by Reference Morales Maqueda and WillmottMorales Maqueda and Willmott (2000) to simulate ice production of an 11 day polynya event that took place in the southern Laptev Sea, Russia, in late December 2007. Model results are compared to ice-thickness and ice-production estimates calculated using high-resolution thermal infrared satellite data obtained from the Moderate Resolution Imaging Spectroradiometer (MODIS) in conjunction with an atmospheric dataset (Reference Riggs, Hall and SalomonsonRiggs and others, 2003; Reference Yu and LindsayYu and Lindsay, 2003). the approach is limited to clear-sky conditions, but yields good results for ice thicknesses less than 0.5 m (Willmes and others, 2010). Based on the comparison, we evaluate model applicability and identify mechanisms in the model that are fairly well represented or require further adjustments.
2. Material and Methods
Below, we briefly introduce the concept of the 2-D time-dependent polynya flux model, the methodology to derive ice production from thermal infrared satellite data and the applied atmospheric dataset. the polynya event was chosen based on MODIS and Envisat SAR observations. In late December 2007, 11 days of offshore winds and low temperatures led to the formation of a wide coastal polynya,
located directly north of the Lena Delta, Russia: the Anabar–Lena polynya. Fast-ice and pack-ice edges and the polynya area are easily identifiable in SAR images in Figure 1. the Envisat C-band wide-swath data are VV-polarized and cover an area of ~400×800km2 with a spatial resolution of 150×150m2. According to SAR observations, the thin ice formed during 11 days of offshore winds covering ~3.5% of the area of the Laptev Sea. Because relatively calm wind conditions, combined with low air temperatures, limited the development of large open-water areas, passive microwave sea-ice concentration data reveal no or only minor polynya activity.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20171010073821-40448-mediumThumb-S0260305500252270_fig1g.jpg?pub-status=live)
Fig. 1. The upper left panel shows the Laptev Sea and mean recurrent coastal polynya locations (orange line). the red box represents the footprint of three Envisat SAR images, shown in the upper right, lower left and lower right panels. the scenes were acquired between 21 and 28 December 2007, and cover the position of the Anabar–Lena polynya, showing the fast-ice belt (south of white line), the active polynya zone and a region of freely floating pack ice (north of dashed white line).
2.1. Model description
The applied 2-D flux model was developed by Reference Morales Maqueda and WillmottMorales Maqueda and Willmott (2000) to simulate the location and temporal evolution of wind-driven polynyas. This model builds upon earlier 1-D unsteady and 2-D steady theories (Reference PeasePease, 1987; Reference Darby, Willmott and MysakDarby and others, 1994, Reference Darby, Wilmott and Somerville1995; Reference Willmott and MoraleWillmott and Morales Maqueda, 1997).
The model distinguishes between two regions in a wind-driven polynya (Fig. 2). Region 1 is the area of open water where frazil ice is formed. Region 2 is the thin-ice area. Region 2 starts where the frazil ice, arriving from region 1, consolidates to a new ice layer, and ends at the first-year ice pack. the term polynya includes both the open-water area and the thin-ice zone. the boundary between regions 1 and 2 is termed the ‘open-water edge’.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20171010073821-15887-mediumThumb-S0260305500252270_fig2g.jpg?pub-status=live)
Fig. 2. Schematic drawing illustrating the polynya model in the 1-D (a) and 2-D (b) cases (adapted from Reference Willmott and MoraleWillmott and Morales Maqueda, 1997). In the open-water area, frazil ice grows with rate ∂h/∂t, and is then herded downwind with speed u, until it arrives with thickness hc at the open-water edge, C(R,t) = constant. Here it piles up to a thickness H, and then drifts as consolidated ice away from the edge with speed U.
The open-water edge is represented by a curve C(R,t) = 0, where R is the position vector of a point on the open-water edge. From Reference Morales Maqueda and WillmottMorales Maqueda and Willmott, (2000) the evolution of point R is determined by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20171010071609300-0538:S0260305500252270:S0260305500252270_eqn1.gif?pub-status=live)
where hc and u are the thickness and velocity of frazil ice arriving at the edge, C. H is the thickness of the consolidated ice formed when the frazil ice ‘piles up’ against C, while U is the drift velocity of the consolidated new ice.
According to Reference Morales Maqueda and WillmottMorales Maqueda and Willmott (2000), the evolution of H, U, hc and u is determined as follows. When the wind blows offshore, the pack ice drifts away from the fast-ice edge and open water is formed. Within the open-water zone, frazil ice formation takes place. Following Reference Martin and KauffmanMartin and Kauffman (1981), the frazil ice production in the open-water area is calculated by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20171010071609300-0538:S0260305500252270:S0260305500252270_eqn2.gif?pub-status=live)
where ρ f is the frazil ice density of 950 kgm–3 (Reference Martin and KauffmanMartin and Kauffman, 1981) and L s = 234.14 kJ kg–1 is the latent heat of fusion for sea ice (Reference Yen, Cheng, Fukusako, Zarling and FaussettYen and others, 1991). the surface heat balance, Q net, can be decomposed as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20171010071609300-0538:S0260305500252270:S0260305500252270_eqn3.gif?pub-status=live)
where Q s and Q l are the turbulent sensible and latent heat fluxes calculated using a bulk parameterization with a transfer coefficient of 2.0×10–3 (Reference Cavalieri and MartinCavalieri and Martin, 1994; Reference Morales Maqueda, Willmott and BiggsMorales Maqueda and others, 2004), Q lw is the net longwave radiation of the sea surface, Q sw is the shortwave solar radiation and α is the surface albedo (Reference Cavalieri and MartinCavalieri and Martin, 1994; Reference Haarpaintner, Haugan and GascardHaarpaintner and others, 2001; Reference Morales Maqueda, Willmott and BiggsMorales Maqueda and others, 2004). If Q net is negative, the water body emits heat to the atmosphere and frazil ice is produced. In late December, the incoming solar radiation and the heat flux due to precipitation can be neglected. the atmospheric forcing is assumed to be uniform over regions 1 and 2.
After formation, Langmuir circulation, created by the interaction of the waves with the wind stress, herds frazil ice into slurries taking the form of long bands or plumes. Wind pushes the long bands of grease ice downwind (Reference PeasePease, 1987; Reference Martin, Steele, Turekian and ThorpeMartin, 2001). the velocity, u, of the frazil ice drift is determined following Reference Martin and KauffmanMartin and Kauffman (1981). the effect of surface currents on the frazil ice drift trajectories is neglected.
At the open-water edge, hc consolidates to a thin-ice layer with thickness H. We apply the parameterization for H developed by Reference Biggs, Morales and WillmottBiggs and others (2000) to avoid situations where hc becomes larger than H. Following Reference Biggs, Morales and WillmottBiggs and others (2000), H is derived in terms of the depth of frazil ice arriving at C and the component, normal to C, of the frazil ice velocity relative to the velocity of the consolidated new ice.
After consolidation, the continuous thermodynamic growth of new ice formed is calculated by Stefan’s law (Reference Maykut and UntersteinerMaykut, 1986; Reference Petrich, Eicken, Thomas and DieckmannPetrich and Eicken, 2009). Through some simplifications, the so-called degree-day model can make a fairly accurate prediction of sea-ice growth. the drift of the new ice is determined by Zubov’s law (Reference LeppärantaLeppäranta, 2005). the occurrence of rafting inside the thin-ice area is not taken into account.
For a comprehensive description of the model and numerical methods employed, we refer the reader to Reference Morales Maqueda and WillmottMorales Maqueda and Willmott (2000) and Reference Willmott, Holland, Morale, Smith and BarberWillmott and others (2007).
2.2. Thermal infrared satellite data
We use level 1B calibrated radiances (thermal infrared) from the MODIS sensor, provided by the NASA Level 1 and Atmosphere Archive and Distribution System (LAADS), to derive ice surface temperatures from brightness temperatures at 11 and 12 μm (MODIS channels 31 and 32, respectively; Reference Riggs, Hall and SalomonsonRiggs and others, 2003). the spatial resolution of the sensor is approximately 1×1 km2. the MODIS cloud mask (MOD35) product is used at the ‘confident clear’ level to determine clear-sky regions and to avoid a potential temperature bias due to sea smoke.
Following Reference Yu and LindsayYu and Lindsay (2003), the thermal thin-ice thickness, h th, is inferred from the obtained ice surface temperatures, T s, by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20171010071609300-0538:S0260305500252270:S0260305500252270_eqn4.gif?pub-status=live)
where T 0 = –1.86˚C is the sea-water temperature at freezing point and k i = 2.03Wm–1 K–1 is the thermal conductivity of sea ice (Reference Drucker, Martin and MoritzDrucker and others, 2003). Q net is calculated from Equation (3), with the same parameterizations as used for the model. the thickness retrieval is based on the assumption that the heat flux through the ice equals the total atmospheric heat flux. the method implies that vertical temperature profiles within the ice are linear and no snow is present on top. It yields good results for ice thicknesses less than 0.5 m, for which its accuracy is estimated to be ±20% (Reference Drucker, Martin and MoritzDrucker and others, 2003).
In total, we identified 67 MODIS scenes, covering the Anabar–Lena polynya between 20 and 30 December 2007. the study period is characterized by very low cloud coverage and the absence of snowfall, such that surface temperature composites could be assembled on a daily basis, except for 29 December.
Assuming that surface heat loss is entirely used for ice formation, the ice production for ice thinner than 0.5 m is determined by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20171010071609300-0538:S0260305500252270:S0260305500252270_eqn5.gif?pub-status=live)
where ρ i = 920 kgm–3 is the density of thin ice. Willmes and others (2010) provide a detailed description of ice-thickness retrieval from surface temperature information.
2.3. Atmospheric data
Model and MODIS computations are based on 6 hourly analysis data from the global numerical weather prediction model, GME, of the German Weather Service (Reference MajewskiMajewski and others, 2002). Mean sea-level pressure, 2m air temperature and humidity and 10m wind vector were extracted from a single gridpoint in the centre of the polynya (Fig. 3). the applicability of the analyzed data is shown by a comparison with automatic weather stations (Reference Schröder, Heinemann and WillmesSchröder and others, in press).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20171010073821-75560-mediumThumb-S0260305500252270_fig3g.jpg?pub-status=live)
Fig. 3. 2 m air temperatures and 10 m wind vectors during the study period, taken from atmospheric dataset. Wind vectors show the direction of airflow, with their lengths representing wind speed.
3. Results
The evolution of open water and thin ice, and associated ice production was simulated with the polynya flux model. Model results were interpolated to a grid with a spacing equal to the spatial resolution of the MODIS sensor. the drift algorithm for consolidated ice (Zubov’s law) requires the ice-drift deviation angle from the wind direction and ice-drift velocity in per cent of wind velocity to be known. These constants were taken from SAR images (Reference Haarpaintner, Haugan and GascardHaarpaintner and others, 2001). the model boundary is the fast-ice edge. Broken fast-ice pieces at the polynya edge hamper comparison of model results with MODIS thermal estimates. Hence, we restrict comparison to the 100 km wide polynya centre covered solely by thin ice. the area is located between 124˚ E and 127˚ E. the method to infer thermal infrared thickness and production is restricted to ice thinner than 0.5 m. We therefore limit the model evaluation to this thermal range (0–0.5 m). Model thicknesses and resultant ice production exceeding 0.5 m are excluded from comparison. Note that for 29 December a MODIS composite is not available. For this particular date, we interpolate between satellite approximations made on 28 and 30 December.
Figure 4 shows flux model results and satellite estimates for 25 December. Polynya ice-thickness distribution and ice production were computed from MODIS data and model results. A SAR image taken on 25 December at 0228 UTC (Fig. 1) shows the areal extent of the polynya. the banded zones of high and low backscatter between the fast-ice edge and pack ice correspond to young ice, formed between 20 and 25 December. These patterns are partially preserved in thermal MODIS ice-thickness information. the mean thickness for the polynya area below 0.5 m is 0.32 m, with minimum ice thicknesses of 0.04 m. the ice area adjacent to the pack-ice edge (cf. SAR scene) is thicker than 0.5 m, and consequently not included in the thermal approach. the mean model polynya thickness is lower (0.27 m), and the maximum thickness does not exceed 0.4 m. According to Figure 4, the agreement between MODIS and model ice production is high and the shape of the polynya is well resolved by the model. Note that the opening event is not visible in sea-ice concentration data derived from the Advanced Microwave Scanning Radiometer (AMSR-E).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20171010073821-58342-mediumThumb-S0260305500252270_fig4g.jpg?pub-status=live)
Fig. 4. Polynya flux model results and satellite estimates for 25 December. (a, b) Ice-thickness distribution inside the thin-ice area, from MODIS approximations (a) and model (b). (c) A SAR scene, taken at 0002 UTC. (d, e) Satellite-based (d) and modelled (e) daily ice volume production.
Below, we compare satellite and model polynya area evolution, thickness distribution and ice-production rates for ice thinner than 0.5 m.
The polynya area evolution, as retrieved from model, MODIS and Envisat SAR estimates, is presented in Figure 5. the thermal and model area estimates show the evolution of ice areas thinner than 0.5 m. In contrast, SAR area estimates include, based on a visual image segmentation, the entire thin-ice area formed between fast ice and the pack-ice edge during 11 days of polynya activity (see Fig. 1).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20171010073821-71888-mediumThumb-S0260305500252270_fig5g.jpg?pub-status=live)
Fig. 5. Total polynya area (light-blue shaded) and open-water fraction (blue shaded) calculated with the flux model. MODIS-based polynya area estimates are presented in the red shaded area. Diamonds show area information taken from Envisat SAR images.
According to the last SAR image available before the end of the observation period, the areal extent of the polynya reaches ~7.9×103 km2. This corresponds to an average opening width of 87 km. During the first 8 days of offshore winds (see Fig. 3), polynya area continuously increases regardless of the methodology used for area calculation. MODIS area approximations are slightly lower than SAR area estimates, since the thickness of thin ice close to the pack-ice edge is beyond the limit of the thermal ice-thickness retrieval. Ice growth beyond 0.5 m also causes the sharp drop in MODIS area from 4.8×103 km2 down to 1.1×103 km2 between 28 and 30 December. Model thin-ice areas exceed ice thicknesses of 0.5 m only at the very end of the observation period. Hence, the model area closely follows SAR area estimates.
A direct comparison of the computed open-water area with satellite approximations is not possible, since reliable satellite-based open-water estimates are missing. In SAR images, often no clear limit exists between open water and consolidated thin ice, due to the smooth transition from low to high frazil ice concentration. MODIS-based open-water estimates suffer from mixed pixel signatures. Following Equation (4), open-water areas are characterized by a surface temperature close to the freezing point of sea water. Unfortunately, the presence of frazil ice, Langmuir streaks or ice floes in the water column lowers the averaged surface temperature within the sensor footprint beneath freezing point, such that potential open-water areas are classified as very thin ice.
Figure 6 shows the zonally averaged ice-thickness profile across the polynya and its standard deviation, computed from MODIS composites and model results for 9 days of polynya activity. Both MODIS and model show an increase in ice thickness with increasing distance from the fast-ice edge. the simulated ice thickness tends to be generally lower than thermal ice thickness. Standard deviations are higher for thermal ice thicknesses than for model thickness.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20171010073821-97754-mediumThumb-S0260305500252270_fig6g.jpg?pub-status=live)
Fig. 6. Zonally averaged polynya ice-thickness distribution as a function of distance from fast-ice edge for 9 days of polynya activity. the red dotted curve represents MODIS thermal ice thickness, and the blue solid curve flux model estimations. Corresponding shaded areas indicate one standard deviation from mean thickness.
Near the fast-ice edge, the model shows a small strip of open water with ice thickness equal to 0. Within this zone, ice growth in the model takes place as frazil ice formation, which is herded offshore until it consolidates at the open-water edge. Following the parameterization developed by Reference Biggs, Morales and WillmottBiggs and others (2000), the consolidation thickness varies between 6 cm (27 December) and 15cm (25 December). After consolidation, the ice thickness in the model is governed by thermodynamic processes as described by Stefan’s law (Reference Petrich, Eicken, Thomas and DieckmannPetrich and Eicken, 2009).
The MODIS sensor does not show open-water zones. Instead, the zone close to the fast-ice edge is characterized by very low thicknesses varying between 6 cm (30 December) and 20 cm (22 December). With increasing distance from the fast-ice edge, the difference in MODIS and model estimates decreases, except for 21 and 22 December.
Figure 7 shows the zonally averaged ice production across the polynya computed from MODIS and model results. Modelled ice production agrees well with MODIS estimates. At the fast-ice edge, highest ice-production rates are achieved, either in the form of frazil ice production (model), or by means of conductive heat flux through very thin ice (MODIS; Fig. 6). Ice production then decreases with distance from the fast-ice edge, owing to the increasing ice thickness. Largest deviations in MODIS and model ice production estimates occur 0–10 km offshore of the fast-ice edge.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20171010073821-10595-mediumThumb-S0260305500252270_fig7g.jpg?pub-status=live)
Fig. 7. Zonally averaged 24 hour ice volume production as a function of distance from fast-ice edge for 9 days of polynya activity. the red dotted curve corresponds to ice production as derived from MODIS. the blue solid curve represents flux model estimations. Corresponding shaded areas indicate one standard deviation from mean thickness.
The accumulated total ice production over 11 days in the entire polynya area is presented in Figure 8. Because modelled ice thickness exceeds 0.5m only at the end of the observation period, its ice-production area extends far north (see also Fig. 5). MODIS and the flux model calculate the accumulated ice production per km2 in the area near the fast-ice edge to be 1.3 and 1.2 m, respectively. Beyond the fast-ice edge region, ice formation rates are higher in the model, because of lower ice-thickness estimates (Fig. 6). the mean accumulated ice production per km2 is 0.4 and 0.37 m, for MODIS and model respectively. the satellite-observed characteristics of ice production are well captured by the model.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20171010073821-22314-mediumThumb-S0260305500252270_fig8g.jpg?pub-status=live)
Fig. 8. Total amount of ice produced in 11 days of polynya activity (20–30 December 2007). Left panel: ice production based on MODIS thermal infrared data. Right panel: model-simulated ice production. In the area near the fast-ice edge the accumulated ice production amounts to 1.3 m (MODIS) and 1.2 m (flux model). the dashed black line shows the extent of the mask area applied in Figure 9.
Figure 9 presents the MODIS and model mean daily and accumulated total ice production. Calculations in the upper panel are limited to an area of 900 km2 adjacent to the fast-ice edge. This area includes open water and parts of the thin-ice zone. the extent of the mask area is drawn in Figure 8 and is equivalent to the polynya area 1.5 days after opening. During the first 8 days, the modelled mean daily and accumulated total production slightly exceeds MODIS estimates. However, an increase in MODIS ice production on 28 and 30 December causes the accumulated ice volumes to converge towards the end of the observation period. the total ice production inside the mask area is 0.9 km3.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20171010073821-39261-mediumThumb-S0260305500252270_fig9g.jpg?pub-status=live)
Fig. 9. MODIS and model mean daily and accumulated total ice production between 20 and 30 December 2007. Upper panel: Mean daily and accumulated total ice production within a predefined 900 km2 area offshore the fast-ice edge (mask area is shown in Fig. 8). Modelled and thermal mean daily ice production are represented by the solid red and dashed red curves, respectively (left axis). Modelled and thermal accumulated total ice production are represented by the solid black and dashed black curves, respectively (right axis). Lower panel: Mean and accumulated total ice production for the entire polynya area. the color and line code are the same as in the upper panel.
In the lower panel, all ice areas thinner than 0.5 m were taken into account for volume approximations. Again, at the beginning the model shows a higher daily mean production, until on 28 December the MODIS polynya area decreases owing to an increase in ice thickness beyond 0.5 m (Fig. 5). As a consequence, modelled daily ice production becomes twice as high as MODIS estimates, resulting in a 30% higher total ice volume.
4. Discussion
The determination of the open-water edge is a central element in the model concept. Unfortunately, the reliability of the modelled open-water width is difficult to assess directly, since for the prescribed situation satellite-based open-water width estimates are missing. This is because the identification of open water in SAR and MODIS scenes is hampered by the narrowness of the open-water zone and the presence of frazil ice in the water column.
Nevertheless, since ice-production rates are highest inside the open-water area, the spatial distribution of highly productive zones indicates the extent of open water. Consequently, the achieved accuracy of the open-water simulation can be determined through a comparison of satellite-based estimates with modelled zones of high ice-growth rates. Moreover, good agreement between modelled and thermal ice-thickness distribution is a second indication for model reliability.
An initial comparison of thermal infrared ice-thickness and -production estimates with model computations has shown that the flux model reproduces quite well the shape of the polynya and salient features of the landfast ice boundary.
Note that the equations to derive ice thickness and ice production from MODIS are, by definition, very close to the equations used in the model. In addition, the applied atmospheric forcing is identical. Hence, the accuracy of the simulation primarily depends on whether a realistic fast-ice boundary is used and whether the movement of the pack-ice edge is prescribed correctly (by estimating appropriate parameters for Zubov’s law from SAR images). Another well-known uncertainty in polynya simulations is the parameterization of the collection depth, H, to which the evolution of the open-water edge is extremely sensitive. Early 1-D and 2-D polynya flux models used constant values for H, typically 0.02–0.3m (Reference Martin and KauffmanMartin and Kauffman, 1981; Reference PeasePease, 1987; Reference Morales Maqueda and WillmottMorales Maqueda and Willmott, 2000; Reference Haarpaintner, Haugan and GascardHaarpaintner and others, 2001), chosen by fitting flux model results to the observed widths of the open-water zone. Reference Winsor and BjörkWinsor and Björk (2000) assume H to be a linear function of the 10m wind speed, increasing from 10 to 30 cm as the wind speed increases from 5 to 35 ms–1. Reference Biggs, Morales and WillmottBiggs and others (2000) developed a parameterization to avoid situations where the thickness of frazil ice arriving at the polynya edge exceeds H. All different parameterizations for H were tested in the model and results were compared to satellite approximations (not shown here). For the event investigated in this paper, we found that the parameterization developed by Reference Biggs, Morales and WillmottBiggs and others (2000) yields the most realistic open-water extent. Parameterizations using constant H are either numerically unstable or lead to a crude overestimation of open-water area.
Furthermore, regional differences between satellite estimates and model results are related to the representation of physical processes in the model and possible deficiencies in the atmospheric datasets. Below, we discuss parameterizations and mechanisms that are fairly well represented or require further adjustments.
Because the GME dataset is spatially quite uniform, we use atmospheric forcing fields extracted from a single gridpoint in the centre of the polynya. the mean deviation in alongshore wind velocity from the polynya centre is 0.4 ±0.21 ms–1. the average alongshore differences in wind direction and temperature are 4.4 ±3.1˚C and 0.34±0.32˚C, respectively. the resultant uncertainty in heat-flux computations is 35 ±30Wm–2, which is equivalent to a frazil ice production error of 1.1 ±0.9cmm–2 d–1. In reality, however, strong gradients of air temperature and wind speed across and along the polynya may exist which are not captured by the GME dataset and which may lead to uncertainties in ice-production estimates. Furthermore, in the atmospheric GME model, polynyas are not present, so the 2m air temperatures are valid over thin ice only. Following Reference Ebner, Schreoder and HeinemannEbner and others (in press), 2 m air temperatures over open water can be significantly higher than over the ice. This may result in an overestimation of both modelled and satellite-based ice production.
Some physical processes that do have an impact on the formation and maintenance of the polynya are not considered in the model. There is evidence that frazil ice trajectories are affected by ocean currents (Willmott and Morales Maqueda, 1997). In the flux model, we assume frazil ice motion to be purely wind-driven, since the strength and direction of near-surface currents are unknown. In addition, inertial motions and tidal currents might advect warmer water from outside into the polynya area, thus hindering frazil ice production (Reference Willmott, Holland, Morale, Smith and BarberWillmott and Morales Maqueda, 2007). the error associated with this uncertainty is unknown.
The simplifications made in the model become increasingly important as the open-water region grows in size. In our case, the simulated polynya develops under moderate wind speeds and low air temperatures, such that the open-water area stays relatively small. If the open-water width increases, which is the case for polynya events in the Laptev Sea taking place in late April (Willmes and others, 2010), spatial variations in atmospheric and oceanic forcing will significantly alter frazil ice motion.
Besides the parameterization of H, discrepancies between satellite estimates and model can also be a consequence of the incomplete, or missing, representation of dynamic processes taking place inside the thin-ice area. According to Reference Skogseth, Fer, Haughan, Drange, Dokken, Furevik, Gerdes and BergerSkogseth and others (2005), rafting and ridging in polynyas can play a significant role in ice thickening, the creation of open-water area and ice production. Rafting and ridging primarily result from a compression of the thin-ice zone, owing to differences between the drift of the outer pack ice and the evolution of the open-water edge.
During the opening of the Anabar–Lena polynya, we assume that compression results exclusively in rafting. This assumption is supported by aerial photographs taken on winter expeditions in the Laptev Sea during different flights across active polynyas (Transdrift XIII, Transdrift XV). the aerial images show that the occurrence of rafting increases with distance from the fast ice. This is because, at the onset of an opening event, when the thin-ice zone is rather small, compression operates on a relatively small thin-ice area, leading to high rafting rates. With increasing thin-ice width, the rafting probability weakens. This possibly explains the decrease in discrepancies between satellite approximations and model results towards the end of the observation period. To account for rafting in the model, one could simulate the evolution of the outer pack ice and open-water edge separately and take the difference between drifts as a measure of compression acting on the thin-ice zone. Nevertheless, an implementation of a dynamic thickening component is difficult, since studies on rafting probabilities in thin-ice zones are limited to characterizations of processes such as energy consumption (Reference Worby, Jeffries, Weeks, Morris and JañaWorby and others, 1996; Reference Babko, Rothrock and MaykutBabko and others, 2002).
The spatial and quantitative agreement between modelled and satellite-based total ice volume production is high (Fig. 8). the largest deviations in ice-production estimates coincide with areas of high disparity between simulated and thermal infrared thickness (Figs 6, 7 and 9). If we limit satellite and model comparison to the highly productive zone close to the fast-ice edge (masked area in Fig. 8 and upper panel of Fig. 9), differences in ice volume calculations are even smaller. Taking into account the entire polynya area (lower panel of Fig. 9), the model overestimates ice production by 30%. This is due to the differences in simulated and satellite-based polynya area estimates (Fig. 5), which is, again, a consequence of the underestimated model ice thickness.
5. Conclusion
We have compared two approaches commonly used to estimate ice production by applying them to an observed Laptev Sea polynya event. A 2-D flux model developed by Reference Morales Maqueda and WillmottMorales Maqueda and Willmott (2000) is used to calculate the location of the transition between open water and thin ice, the drift of consolidated new ice and associated ice formation. Results from this model are in agreement with ice-thickness and -production estimates obtained from MODIS thermal infrared data in conjunction with an atmospheric dataset.
Given a polynya boundary (coastline or fast-ice edge) and an atmospheric dataset, the model is capable of describing the shape of the polynya realistically using only a few simple conservation laws. the mean accumulated ice production per unit area for 11 days of polynya activity is 0.4 and 0.37 m for thermal and model estimates, respectively. the accumulated ice production in the area near the fast-ice edge amounts to 1.3 and 1.2 m, respectively. Further offshore, ice-formation rates are overestimated by the model because it underestimates the ice thickness.
We conclude that this simplified physical polynya model is a valuable tool for studying polynya dynamics and estimating associated fluxes of single polynya events. It provides an alternative to passive microwave polynya-monitoring techniques, which are often not capable of resolving narrow polynyas formed under moderate or low wind forcing.
The best model results are achieved with the consolidation thickness parameterization developed by Reference Biggs, Morales and WillmottBiggs and others (2000). the implementation of dynamic thin-ice thickening could eventually improve observed regional discrepancies between model and satellite. When simulating polynya events of larger alongshore dimension, the use of a spatially non-uniform forcing dataset would be required. However, with increasing complexity, a purely numerical approach is needed.
Acknowledgements
This work was carried out as part of the German–Russian cooperation ‘System Laptev Sea’, funded by the Federal Ministry of Education and Research (BMBF), Germany, under grant 03G0639A and the Alfred Wegener Institute. Satellite data were obtained through European Space Agency project EO-500 ‘Formation, transport and distribution of sediment-laden sea ice in the Arctic Shelf seas’. MODIS images were provided by LAADS/NASA, and AMSR-E ARTIST Sea Ice (ASI) sea concentrations by the University of Bremen, Germany. the GME analyses were made available by the German Weather Service.