Introduction
Arctic sea ice has undergone a dramatic decline in recent years, with a well-documented retreat of the summer ice cover (Reference Serreze, Holland and StroeveSerreze and others, 2007; Reference Comiso, Parkinson, Gersten and StockComiso and others, 2008), a general thinning (Reference Giles, Laxon and RidoutGiles and others, 2008; Reference Haas, Pfaffling, Rabenstein, Etienne and RigorHaas and others, 2008; Reference Rothrock, Percival and WensnahanRothrock and others, 2008; Reference Kwok, Cunningham, Wensnahan, Zwally and YiKwok and others, 2009) and a transition to a younger, more vulnerable ice pack (Reference Maslanik, Fowler, Drobot and ZwallyMaslanik and others, 2007; Reference Nghiem, Rigor, Perovich, Weatherly and NeumannNghiem and others, 2007). A record September minimum ice extent was observed in 2007 (Reference StroeveStroeve and others, 2008) and was followed by the second smallest observed ice extent in 2008 (US National Snow and Ice Data Center, http://nside.org/arcticseaicenews/). Climate models project that these changes may continue, with a possible transition to ice-free summers later this century (Reference Wang and OverlandWang and Overland, 2009).
These observed changes are intricately linked to sea-ice dynamics and thermodynamics and are driven by atmosphere and ocean forcing. Several factors have been established as contributors to the decline in the ice cover: a general warming (Reference Richter-MengeRichter-Menge, 2009), changes in atmospheric circulation patterns (Reference Rigor and WallaceRigor and Wallace, 2004), changes in cloudiness (Reference Francis and HunterFrancis and Hunter, 2006; Reference Kay, Gettelman, Stephens and O’DellKay and others, 2008), advected ocean heat from lower latitudes (Reference Polyakov, Walsh, Colony and TimokhovPolyakov and others, 2003; Reference ShimadaShimada and others, 2006; Reference Woodgate, Aagaard and WeingartnerWoodgate and others, 2006), increased ice export from the Fram Strait (Reference Nghiem, Rigor, Perovich, Weatherly and NeumannNghiem and others, 2007) and increased solar heating of the upper ocean (Reference Perovich, Light, Jones, Runciman and NghiemPerovich and others, 2007, Reference Perovich, Jones and Light2008; Reference Steele, Ermold and ZhangSteele and others, 2008).
Sea-ice mass balance can provide insights into the thermodynamic influences on the declining ice cover. the mass balance is the difference between the amount of ice growth during the winter and the amount of surface and bottom melt during the summer. the mass balance is an integrator of both the surface heat budget and the ocean heat flux and thus provides a way to attribute observed changes in ice extent and ice thickness. In this paper, we examine observations of surface and bottom melting made from seven autonomous ice mass-balance buoys during the summer of 2008. Observed melting is compared to calculated estimates of solar radiation incident on the ice and transmitted directly into the ocean.
Approach
The ice mass balance was measured using autonomous ice mass-balance buoys (IMBs) (Reference Richter-Menge, Perovich, Elder, Claffey, Rigor and OrtmeyerRichter-Menge and others, 2006; Reference Perovich, Jones and LightPerovich and others, 2008). These buoys monitor changes in snow deposition and ablation, ice growth, and ice surface and bottom melt. Two acoustic rangefinders monitor the position of the ice bottom and the snow/ice surface. A string of thermistors measures vertical profiles of ice temperature. Ancillary data from the IMBs include barometric pressure, air temperature (initially 2m above the surface), water temperature (initially ~0.5m below the ice bottom), and position. Data are recorded using a Campbell Scientific data logger and are transmitted using the ARGOS (Advanced Research and Global Observation Satellite) system.
These buoys have routinely been deployed on undeformed multi-year sea ice in the Arctic since 2000, typically in conjunction with other sensor packages studying the atmosphere and ocean. As part of the 2007–09 International Polar Year (IPY), we were able to expand the number of buoys deployed. the IPY deployments were done in collaboration with the North Pole Environmental Observatory (Reference MorisonMorison and others, 2002), the Beaufort Gyre Environmental Observatory and the DAMOCLES program (Developing Arctic Modeling and Observing Capabilities for Long-term Environmental Studies) (Reference GascardGascard and others, 2008).
Results and Discussion
Seven IMBs operated during the summer of 2008. One was installed in 2006, another in 2007 and five in 2008. These buoys were located primarily in the western Arctic. Figure 1 shows the approximate summer position (white dot), the total amount of surface ablation (red bar) and the total bottom melt (yellow bar) for the summer of 2008. All of the buoys were installed in undeformed multi-year ice, with end-of-winter ice thicknesses ranging from 1.83 to 3.17m (Table 1). There was considerable regional variability in the amount of surface, bottom and total melting. the smallest amount of melting was at a location north of Greenland, where 0.3 m of surface ablation and 0.1 m of bottom melt resulted in a total thinning of only 0.4 m. the maximum melting was in the Beaufort Sea, where 3.2 m of thick ice completely melted by 23 August 2008, with 0.87m of surface ablation and 1.77 m of bottom melt measured before the buoy failed. on average, there was 0.47m of surface melting and 0.66m of bottom melting.
The background map in Figure 1 displays ice concentration in September 2008. A casual examination shows a tendency towards greater bottom ablation in regions of lower ice concentration. the largest observed bottom melting was at the ice-edge site in the Beaufort Sea, while the smallest was north of Greenland where the ice concentration was large throughout the summer. There were also large amounts of bottom melting in the Beaufort Sea in the summer of 2007, due in part to enhanced solar heating of the upper ocean (Reference Perovich, Jones and LightPerovich and others, 2008).
We can examine ocean solar heating in 2008 in more detail by calculating estimates of local solar heating at each of the seven buoy locations. the flux of solar heat input directly to the upper ocean (F rw) through leads can be estimated using the relationship
where F r is the incident solar irradiance, α w is the albedo of the ocean and I c is the ice concentration. F r and I c are determined for a cell containing an IMB using a 25 km×25km Equal-Area Scalable Earth Grid (EASE-Grid). the albedo of water was set to 0.07 (Reference Pegau and PaulsonPegau and Paulson, 2001). Values of F r were obtained from adjusted US National Centers for Environmental Prediction (NCEP) values. the area of open water was determined from passive microwave observations by using the NASA Team 2 algorithm (Reference Markus and CavalieriMarkus and Cavalieri, 2000).
Mean daily downwelling shortwave fluxes from the NCEP Reanalysis are known to have large biases due to incorrect cloud fractions, in contrast to those of the European Centre for Medium-Range Weather Forecasts reanalysis (ERA-40; Reference Serreze, Key, Box, Maslanik and SteffenSerreze and others, 1998; Reference Walsh, Chapman and PortisWalsh and others, 2009) which has a more realistic seasonal cycle for the clouds but ends in 2002. However, the fluxes in the two reanalyses are well correlated if seasonal and geographic variations are accounted for (Reference Liu, Curry, Rossow, Key and WangLiu and others, 2005). Seasonal and regional regression equations were determined between the two datasets so that the more recent NCEP Reanalysis data could be adjusted to match the more precise ERA-40 values. the NCEP predictor variables are the downwelling solar flux and the sea-level pressure. Five years of data were used (evenly spaced from 1979 to 1999), and separate regression equations were determined every 5 days of the seasonal cycle for every second gridpoint of the NCEP data. the mean root-mean-square error of the regression fits is 29Wm–2 for the daily means.
Additional solar energy is also transmitted to the upper ocean through the ice, but that amount is neglected in this analysis. the amount of solar heat transmitted through the ice depends primarily on ice thickness and pond fraction, parameters that are not known for our study. Reference PerovichPerovich (2005), using field observations from the Surface Heat Budget of the Arctic Ocean (SHEBA) program and a radiative transfer model, determined that leads were the major source of solar energy deposited in the ocean (61%), but that bare ice (23%) and ponded ice (16%) also contributed significantly.
Equation (1) was evaluated every day from 1 January 2008 to 16 October (end of melt) for each of the seven 2008 IMB sites and integrated over time to obtain the heat input to the upper ocean (Q w). the heat absorbed in the ocean was averaged over the entire ice–ocean area of the gridcell area.
The solar heat input to the ice (F ri) is
where α i is the albedo of the ice. This heat would be available for melting at the ice surface and interior. As stated above, the small fraction that would be transmitted through the ice into the ocean is neglected. the key to evaluating Equation (2) is determining the time-dependent evolution of albedo. We assumed that the ice albedo evolved from cold snow (0.85) to wet snow (0.81) to melting snow (decrease from 0.81 to 0.65) to bare ice (0.65), and then back to snow-covered ice after fall freeze-up (Reference Perovich, Grenfell, Light and HobbsPerovich and others, 2002). the timing of the albedo transitions was determined from satellite estimates of the onset of melt and freeze-up (Reference Markus, Stroeve and MillerMarkus and others, 2009). the method distinguishes between early melt (the first day of the year when melt is detected) and the first day of continuous melt (when melt conditions prevail for the rest of the summer). A similar distinction is made for freeze-up. Fall freeze-up is taken as the day that freeze-up is first detected. Using this approach for albedo, Equation (2) was evaluated daily and then integrated over the duration of the surface melt period for each of the buoys to obtain the heat input to the ice (Q i).
The role of solar heating is illustrated in Figure 2 by comparing results from the IMBs with the least amount of melting (buoy 5) and the most (buoy 2). the Beaufort site had 0.87 m of surface ablation and 1.77m of bottom melt, compared with only 0.24 m of surface ablation and 0.16m of bottom melt at the site north of Greenland. Figure 2 compares time-series results of incident shortwave irradiance, ice concentration, heat input to the ocean, ice albedo, and heat input to the ice for these two IMBs. the incident solar irradiance for the two cases (Fig. 2a) is quite similar, with the curves crossing multiple times. the difference, integrated over the summer, is only a few percent, not nearly enough to explain the large differences in observed melt. There was an enormous difference in ice concentration between the two sites, with the Beaufort Sea site always having a smaller ice concentration and reaching open water in late August. This resulted in the solar heat input to the ocean always being larger at the Beaufort site (Fig. 2c). the total solar heat input to the ocean at the Beaufort site was three times as large as at the site north of Greenland (390 versus 131 MJm–2). Estimated albedos for the two sites are presented in Figure 2d. the main difference is that melting at the Beaufort site started 21 days earlier. This led to a lower albedo for June and greater solar heat input to the ice (Fig. 2e). the ice albedo in June is critical to the surface heat budget, since June is typically the month of peak solar input at the surface.
We can further explore the relationship between solar heat input and melt by comparing the solar heat input to the amount of heat used in melting per square meter on the surface (Q ms) and bottom (Q mb) of the ice. Snowmelt is included in the surface melting calculation:
where ρ s is the density of Arctic snow (330 kgm–3), ρ i is the density of ice (900 kgm–3), L f is the latent heat of fusion of ice (334 kJ kg–1), ΔH s is the amount of snowmelt, ΔH i is the amount of surface ice melt and ΔH b is the amount of bottom ice melt.
Using the observed surface and bottom melting, Equations (3) and (4) were evaluated for each of the buoys. the heat used in surface melting ranged from 66 to 269 MJm–2, while the bottom melting heat varied from 48 to 532 MJm–2 (Table 1). the relationship between the heat used in melting and solar heat input is shown in Figure 3. Examining surface melting, there is a general weak trend relating increased
surface melting with increased solar heat deposited in the ice (Fig. 3a). the slope of the linear relationship is 0.42, with a correlation coefficient of 0.38 and considerable scatter. For example, ~600Wm–2 of solar heat absorbed in the ice resulted in surface melting ranging from 66 to 198Wm–2. This indicates that there are additional factors that influence surface ablation, which is not surprising as sensible heat, latent heat and longwave radiation also impact the surface heat budget. the net longwave radiation, in particular, has been established as a major factor in the surface heat budget and has a cooling effect on the surface in winter (Reference Persson, Fairall, Andreas, Guest and PerovichPersson and others, 2002; Reference Perovich, Grenfell, Light, Tucker and EickenPerovich and others, 2003). In summer the impact of net longwave radiation is more complex, with low clouds usually leading to enhanced surface melt.
The correlation is much stronger (correlation coefficient of 0.94) between the heat used in bottom melting and the solar heat deposited in the ocean (Fig. 3b). the slope of the line is 0.89, indicating an almost one-to-one increase in bottom melting with solar heat input to the ocean. the relationship holds for observations that vary widely in geographic location, ice concentration and bottom melting. This argues that the primary source of heat for bottom melting of the ice is solar radiation absorbed in areas of open water (Reference Maykut and McPheeMaykut and McPhee, 1995; Reference Perovich, Jones and LightPerovich and others, 2008). For all seven of the buoys, there was more than enough solar energy deposited in the ice and in the ocean to provide the heat necessary for surface and bottom melting.
Conclusions
Surface-based observations of sea-ice mass balance at seven locations in the summer of 2008 showed considerable regional variability in surface and bottom melting. Total summer melt ranged from 0.37 to 3.17 m. There was a strong linkage between solar heat input to the ocean and bottom melt. Variations in the solar heat input to the ocean were caused primarily by differences in ice concentration, with the largest amounts of bottom melting occurring in regions of reduced ice concentration.
These results have implications for a declining Arctic sea-ice cover. They indicate that decreases in ice concentration will lead directly to enhanced solar heat input to the ocean and increased bottom melting. Thinning of the ice will increase light transmission through the ice into the ocean, increasing melting in an ice-albedo feedback. In addition, earlier dates of melt onset will result in more solar heat input into the ice and more surface melting. Later dates of freezeup will have a similar result, but to a lesser extent.
Acknowledgements
We appreciate the support of the North Pole Environmental Observatory, the Beaufort Gyre Environmental Observatory and the DAMOCLES program for their efforts in deploying ice mass-balance buoys. This work has been funded by the US National Science Foundation, the National Aeronautics and Space Administration, and the National Oceanographic and Atmospheric Administration.