The Arctic icescape is rapidly transforming from a thicker multiyear ice cover to a thinner and largely seasonal first-year ice cover with significant consequences for Arctic primary production. One critical challenge is to understand how productivity will change within the next decades. Recent studies have reported extensive phytoplankton blooms beneath ponded sea ice during summer, indicating that satellite-based Arctic annual primary production estimates may be significantly underestimated. Here we present a unique time-series of a phytoplankton spring bloom observed beneath snow-covered Arctic pack ice. The bloom, dominated by the haptophyte algae Phaeocystis pouchetii, caused near depletion of the surface nitrate inventory and a decline in dissolved inorganic carbon by 16 ± 6 g C m−2. Ocean circulation characteristics in the area indicated that the bloom developed in situ despite the snow-covered sea ice. Leads in the dynamic ice cover provided added sunlight necessary to initiate and sustain the bloom. Phytoplankton blooms beneath snow-covered ice might become more common and widespread in the future Arctic Ocean with frequent lead formation due to thinner and more dynamic sea ice despite projected increases in high-Arctic snowfall. This could alter productivity, marine food webs and carbon sequestration in the Arctic Ocean.

Abstract. Melt ponds on sea ice strongly reduce the surface albedo and accelerate the decay of Arctic sea ice. Due to different spectral properties of snow, ice, and water, the fractional coverage of these distinct surface types can be derived from multispectral sensors like the Moderate Resolution Image Spectroradiometer (MODIS) using a spectral unmixing algorithm. The unmixing was implemented using a multilayer perceptron to reduce computational costs.Arctic-wide melt pond fractions and sea ice concentrations are derived from the level 3 MODIS surface reflectance product. The validation of the MODIS melt pond data set was conducted with aerial photos from the MELTEX campaign 2008 in the Beaufort Sea, data sets from the National Snow and Ice Data Center (NSIDC) for 2000 and 2001 from four sites spread over the entire Arctic, and with ship observations from the trans-Arctic HOTRAX cruise in 2005. The root-mean-square errors range from 3.8 % for the comparison with HOTRAX data, over 10.7 % for the comparison with NSIDC data, to 10.3 % and 11.4 % for the comparison with MELTEX data, with coefficient of determination ranging from R 2 = 0.28 to R 2 = 0.45. The mean annual cycle of the melt pond fraction per grid cell for the entire Arctic shows a strong increase in June, reaching a maximum of 15 % by the end of June. The zonal mean of melt pond fractions indicates a dependence of the temporal development of melt ponds on the geographical latitude, and has its maximum in mid-July at latitudes between 80 • and 88 • N.Furthermore, the MODIS results are used to estimate the influence of melt ponds on retrievals of sea ice concentrations from passive microwave data. Results from a case study comparing sea ice concentrations from ARTIST Sea Ice-, NASA Team 2-, and Bootstrap-algorithms with MODIS sea ice concentrations indicate an underestimation of around 40 % for sea ice concentrations retrieved with microwave algorithms.

The salinity and water oxygen isotope composition (δ18O) of 29 first‐year (FYI) and second‐year (SYI) Arctic sea ice cores (total length 32.0 m) from the drifting ice pack north of Svalbard were examined to quantify the contribution of snow to sea ice mass. Five cores (total length 6.4 m) were analyzed for their structural composition, showing variable contribution of 10–30% by granular ice. In these cores, snow had been entrained in 6–28% of the total ice thickness. We found evidence of snow contribution in about three quarters of the sea ice cores, when surface granular layers had very low δ18O values. Snow contributed 7.5–9.7% to sea ice mass balance on average (including also cores with no snow) based on δ18O mass balance calculations. In SYI cores, snow fraction by mass (12.7–16.3%) was much higher than in FYI cores (3.3–4.4%), while the bulk salinity of FYI (4.9) was distinctively higher than for SYI (2.7). We conclude that oxygen isotopes and salinity profiles can give information on the age of the ice and enables distinction between FYI and SYI (or older) ice in the area north of Svalbard.

Seven ice mass balance instruments deployed near 83°N on different first‐year and second‐year ice floes, representing variable snow and ice conditions, documented the evolution of snow and ice conditions in the Arctic Ocean north of Svalbard in January–March 2015. Frequent profiles of temperature and thermal diffusivity proxy were recorded to distinguish changes in snow depth and ice thickness with 2 cm vertical resolution. Four instruments documented flooding and snow‐ice formation. Flooding was clearly detectable in the simultaneous changes in thermal diffusivity proxy, increased temperature, and heat propagation through the underlying ice. Slush then progressively transformed into snow‐ice. Flooding resulted from two different processes: (i) after storm‐induced breakup of snow‐loaded floes and (ii) after loss of buoyancy due to basal ice melt. In the case of breakup, when the ice was cold and not permeable, rapid flooding, probably due to lateral intrusion of seawater, led to slush and snow‐ice layers at the ocean freezing temperature (−1.88°C). After the storm, the instruments documented basal sea‐ice melt over warm Atlantic waters and ocean‐to‐ice heat flux peaked at up to 400 W m−2. The warm ice was then permeable and flooding was more gradual probably involving vertical intrusion of brines and led to colder slush and snow‐ice (−3°C). The N‐ICE2015 campaign provided the first documentation of significant flooding and snow‐ice formation in the Arctic ice pack as the slush partially refroze. Snow‐ice formation may become a more frequently observed process in a thinner ice Arctic.

In recent years, sea‐ice conditions in the Arctic Ocean changed substantially toward a younger and thinner sea‐ice cover. To capture the scope of these changes and identify the differences between individual regions, in situ observations from expeditions are a valuable data source. We present a continuous time series of in situ measurements from the N‐ICE2015 expedition from January to June 2015 in the Arctic Basin north of Svalbard, comprising snow buoy and ice mass balance buoy data and local and regional data gained from electromagnetic induction (EM) surveys and snow probe measurements from four distinct drifts. The observed mean snow depth of 0.53 m for April to early June is 73% above the average value of 0.30 m from historical and recent observations in this region, covering the years 1955–2017. The modal total ice and snow thicknesses, of 1.6 and 1.7 m measured with ground‐based EM and airborne EM measurements in April, May, and June 2015, respectively, lie below the values ranging from 1.8 to 2.7 m, reported in historical observations from the same region and time of year. The thick snow cover slows thermodynamic growth of the underlying sea ice. In combination with a thin sea‐ice cover this leads to an imbalance between snow and ice thickness, which causes widespread negative freeboard with subsequent flooding and a potential for snow‐ice formation. With certainty, 29% of randomly located drill holes on level ice had negative freeboard.