Search

There have been considerable efforts to understand the hydrography of Storfjorden (Svalbard). A recurring winter polynya with large sea ice production makes it an important region of dense water formation at the scale of the Arctic Ocean. In addition, this fjord is seasonally influenced by freshwater inputs from sea ice melt and the surrounding islands of the Svalbard archipelago, which impacts the hydrography. However, the understanding of factors controlling the optical properties of the waters in Storfjorden are lacking and are crucial for the development of more accurate regional bio-optical models. Here, we present results from the first detailed optical field survey of Storfjorden conducted in early summer of 2020. Our observations are based on spectrometric analysis of water samples and in situ vertical profiles with an absorption and attenuation meter, a fluorometer, and a conductivity, temperature, and depth (CTD) sensor. In addition to the expected seasonal contribution from phytoplankton, we find that in early summer waters in Storfjorden are optically complex with a significant contribution from coloured dissolved organic matter (CDOM, 33 %–64 % of the non-water absorption at 443 nm) despite relatively low CDOM concentrations and in the nearshore or near the seabed from non-algal particles (up to 61 % of the non-water absorption at 550 nm). In surface waters, the spatial variability of light attenuation was mainly controlled by inorganic suspended matter originating from river runoff. A distinct subsurface maximum of light attenuation was largely driven by a subsurface phytoplankton bloom, controlled by stratification resulting from sea ice melt. Lastly, the cold dense bottom waters of Storfjorden from winter sea ice production periodically overflows the sill at the mouth of the fjord and can thus reach the Fram Strait. It contained elevated levels of both non-algal particles and dissolved organic matter, which are likely caused by the dense flows of the nepheloid layer interacting with the sea bed.

The Arctic marine climate system is changing rapidly, which is seen in the warming of the ocean and atmosphere, decline of sea ice cover, increase in river discharge, acidification of the ocean, and changes in marine ecosystems. Socio-economic activities in the coastal and marine Arctic are simultaneously changing. This calls for the establishment of a marine Arctic component of the Pan-Eurasian Experiment (MA-PEEX). There is a need for more in situ observations on the marine atmosphere, sea ice, and ocean, but increasing the amount of such observations is a pronounced technological and logistical challenge. The SMEAR (Station for Measuring Ecosystem–Atmosphere Relations) concept can be applied in coastal and archipelago stations, but in the Arctic Ocean it will probably be more cost-effective to further develop a strongly distributed marine observation network based on autonomous buoys, moorings, autonomous underwater vehicles (AUVs), and unmanned aerial vehicles (UAVs). These have to be supported by research vessel and aircraft campaigns, as well as various coastal observations, including community-based ones. Major manned drifting stations may occasionally be comparable to terrestrial SMEAR flagship stations. To best utilize the observations, atmosphere–ocean reanalyses need to be further developed. To well integrate MA-PEEX with the existing terrestrial–atmospheric PEEX, focus is needed on the river discharge and associated fluxes, coastal processes, and atmospheric transports in and out of the marine Arctic. More observations and research are also needed on the specific socio-economic challenges and opportunities in the marine and coastal Arctic, and on their interaction with changes in the climate and environmental system. MA-PEEX will promote international collaboration; sustainable marine meteorological, sea ice, and oceanographic observations; advanced data management; and multidisciplinary research on the marine Arctic and its interaction with the Eurasian continent.

We introduce the OSI-450, the SICCI-25km and the SICCI-50km climate data records of gridded global sea-ice concentration. These three records are derived from passive microwave satellite data and offer three distinct advantages compared to existing records: first, all three records provide quantitative information on uncertainty and possibly applied filtering at every grid point and every time step. Second, they are based on dynamic tie points, which capture the time evolution of surface characteristics of the ice cover and accommodate potential calibration differences between satellite missions. Third, they are produced in the context of sustained services offering committed extension, documentation, traceability, and user support. The three records differ in the underlying satellite data (SMMR & SSM/I & SSMIS or AMSR-E & AMSR2), in the imaging frequency channels (37 GHz and either 6 or 19 GHz), in their horizontal resolution (25 or 50 km), and in the time period they cover. We introduce the underlying algorithms and provide an evaluation. We find that all three records compare well with independent estimates of sea-ice concentration both in regions with very high sea-ice concentration and in regions with very low sea-ice concentration. We hence trust that these records will prove helpful for a better understanding of the evolution of the Earth's sea-ice cover.

The snow data from the Soviet airborne expeditions Sever in the Arctic collected over several decades in March, April and May have been analyzed in this study. The Sever data included more measurements and covered a much wider area, particularly in the Eurasian marginal seas (Kara Sea, Laptev Sea, East Siberian Sea and Chukchi Sea), compared to the Soviet North Pole drifting stations. The latter collected data mainly in the central part of the Arctic Basin. The following snow parameters have been analyzed: average snow depth on the level ice (undisturbed snow) height and area of sastrugi, depth of snow dunes attached to ice ridges and depth of snow on hummocks. In the 1970s–1980s, in the central Arctic, the average depth of undisturbed snow was 21.2 cm, the depth of sastrugi (that occupied about 30 % of the ice surface) was 36.2 cm and the average depth of snow near hummocks and ridges was about 65 cm. For the marginal seas, the average depth of undisturbed snow on the level ice varied from 9.8 cm in the Laptev Sea to 15.3 cm in the East Siberian Sea, which had a larger fraction of multiyear ice. In the marginal seas the spatial variability of snow depth was characterized by standard deviation varying between 66 and 100 %. The average height of sastrugi varied from 23 cm to about 32 cm with standard deviation between 50 and 56 %. The average area covered by sastrugi in the marginal seas was estimated to be 36.5 % of the total ice area where sastrugi were observed. The main result of the study is a new snow depth climatology for the late winter using data from both the Sever expeditions and the North Pole drifting stations. The snow load on the ice observed by Sever expeditions has been described as a combination of the depth of undisturbed snow on the level ice and snow depth of sastrugi weighted in proportion to the sastrugi area. The height of snow accumulated near the ice ridges was not included in the calculations because there are no estimates of the area covered by those features from the Sever expeditions. The effect of not including that data can lead to some underestimation of the average snow depth. The new climatology refines the description of snow depth in the central Arctic compared to the results by Warren et al. (1999) and provides additional detailed data in the marginal seas. The snow depth climatology is based on 94 % Sever data and 6 % North Pole data. The new climatology shows lower snow depth in the central Arctic comparing to Warren climatology and more detailed data in the Eurasian seas.

An open-source sea ice drift algorithm for Sentinel-1 SAR imagery is introduced based on the combination of feature tracking and pattern matching. Feature tracking produces an initial drift estimate and limits the search area for the consecutive pattern matching, which provides small- to medium-scale drift adjustments and normalised cross-correlation values. The algorithm is designed to combine the two approaches in order to benefit from the respective advantages. The considered feature-tracking method allows for an efficient computation of the drift field and the resulting vectors show a high degree of independence in terms of position, length, direction and rotation. The considered pattern-matching method, on the other hand, allows better control over vector positioning and resolution. The preprocessing of the Sentinel-1 data has been adjusted to retrieve a feature distribution that depends less on SAR backscatter peak values. Applying the algorithm with the recommended parameter setting, sea ice drift retrieval with a vector spacing of 4 km on Sentinel-1 images covering 400 km × 400 km, takes about 4 min on a standard 2.7 GHz processor with 8 GB memory. The corresponding recommended patch size for the pattern-matching step that defines the final resolution of each drift vector is 34 × 34 pixels (2.7 × 2.7 km). To assess the potential performance after finding suitable search restrictions, calculated drift results from 246 Sentinel-1 image pairs have been compared to buoy GPS data, collected in 2015 between 15 January and 22 April and covering an area from 80.5 to 83.5° N and 12 to 27° E. We found a logarithmic normal distribution of the displacement difference with a median at 352.9 m using HV polarisation and 535.7 m using HH polarisation. All software requirements necessary for applying the presented sea ice drift algorithm are open-source to ensure free implementation and easy distribution.

The question of how to derive and present uncertainty information in climate data records (CDRs) has received sustained attention within the European Space Agency Climate Change Initiative (CCI), a programme to generate CDRs addressing a range of essential climate variables (ECVs) from satellite data. Here, we review the nature, mathematics, practicalities, and communication of uncertainty information in CDRs from Earth observations. This review paper argues that CDRs derived from satellite-based Earth observation (EO) should include rigorous uncertainty information to support the application of the data in contexts such as policy, climate modelling, and numerical weather prediction reanalysis. Uncertainty, error, and quality are distinct concepts, and the case is made that CDR products should follow international metrological norms for presenting quantified uncertainty. As a baseline for good practice, total standard uncertainty should be quantified per datum in a CDR, meaning that uncertainty estimates should clearly discriminate more and less certain data. In this case, flags for data quality should not duplicate uncertainty information, but instead describe complementary information (such as the confidence in the uncertainty estimate provided or indicators of conditions violating the retrieval assumptions). The paper discusses the many sources of error in CDRs, noting that different errors may be correlated across a wide range of timescales and space scales. Error effects that contribute negligibly to the total uncertainty in a single-satellite measurement can be the dominant sources of uncertainty in a CDR on the large space scales and long timescales that are highly relevant for some climate applications. For this reason, identifying and characterizing the relevant sources of uncertainty for CDRs is particularly challenging. The characterization of uncertainty caused by a given error effect involves assessing the magnitude of the effect, the shape of the error distribution, and the propagation of the uncertainty to the geophysical variable in the CDR accounting for its error correlation properties. Uncertainty estimates can and should be validated as part of CDR validation when possible. These principles are quite general, but the approach to providing uncertainty information appropriate to different ECVs is varied, as confirmed by a brief review across different ECVs in the CCI. User requirements for uncertainty information can conflict with each other, and a variety of solutions and compromises are possible. The concept of an ensemble CDR as a simple means of communicating rigorous uncertainty information to users is discussed. Our review concludes by providing eight concrete recommendations for good practice in providing and communicating uncertainty in EO-based climate data records.

Synthetic Aperture Radar (SAR) data from RADARSAT-2 (RS2) in dual-polarization mode provide additional information for discriminating sea ice and open water compared to single-polarization data. We have developed an automatic algorithm based on dual-polarized RS2 SAR images to distinguish open water (rough and calm) and sea ice. Several technical issues inherent in RS2 data were solved in the pre-processing stage, including thermal noise reduction in HV polarization and correction of angular backscatter dependency in HH polarization. Texture features were explored and used in addition to supervised image classification based on the support vector machines (SVM) approach. The study was conducted in the ice-covered area between Greenland and Franz Josef Land. The algorithm has been trained using 24 RS2 scenes acquired in winter months in 2011 and 2012, and the results were validated against manually derived ice charts of the Norwegian Meteorological Institute. The algorithm was applied on a total of 2705 RS2 scenes obtained from 2013 to 2015, and the validation results showed that the average classification accuracy was 91 ± 4 %.

A computationally efficient, open-source feature-tracking algorithm, called ORB, is adopted and tuned for sea ice drift retrieval from Sentinel-1 SAR (Synthetic Aperture Radar) images. The most suitable setting and parameter values have been found using four Sentinel-1 image pairs representative of sea ice conditions between Greenland and Severnaya Zemlya during winter and spring. The performance of the algorithm is compared to two other feature-tracking algorithms, namely SIFT (Scale-Invariant Feature Transform) and SURF (Speeded-Up Robust Features). Having been applied to 43 test image pairs acquired over Fram Strait and the north-east of Greenland, the tuned ORB (Oriented FAST and Rotated BRIEF) algorithm produces the highest number of vectors (177 513, SIFT: 43 260 and SURF: 25 113), while being computationally most efficient (66 s, SIFT: 182 s and SURF: 99 s per image pair using a 2.7 GHz processor with 8 GB memory). For validation purposes, 314 manually drawn vectors have been compared with the closest calculated vectors, and the resulting root mean square error of ice drift is 563 m. All test image pairs show a significantly better performance of the HV (horizontal transmit, vertical receive) channel due to higher informativeness. On average, around four times as many vectors have been found using HV polarization. All software requirements necessary for applying the presented feature-tracking algorithm are open source to ensure a free and easy implementation.

A satellite database including 16 555 satellite images and ice charts displaying the area of Isfjorden, Hornsund, and the Svalbard region has been established with focus on the time period 2000–2014. 3319 manual interpretations of sea ice conditions have been conducted, resulting in two time series dividing the area of Isfjorden and Hornsund into "fast ice" (sea ice attached to the coastline), "drift ice", and "open water". The maximum fast ice coverage of Isfjorden is > 40 % in the periods 2000–2005 and 2009–2011 and stays 40 % in all considered years, except for 2012 and 2014, where the maximum stays

We assess different methods and input parameters, namely snow depth, snow density and ice density, used in freeboard-to-thickness conversion of Arctic sea ice. This conversion is an important part of sea ice thickness retrieval from spaceborne altimetry. A data base is created comprising sea ice freeboard derived from satellite radar altimetry between 1993 and 2012 and co-locate observations of total (sea ice + snow) and sea ice freeboard from the Operation Ice Bridge (OIB) and CryoSat Validation Experiment (CryoVEx) airborne campaigns, of sea ice draft from moored and submarine upward looking sonar (ULS), and of snow depth from OIB campaigns, Advanced Microwave Scanning Radiometer (AMSR-E) and the Warren climatology (Warren et al., 1999). We compare the different data sets in spatiotemporal scales where satellite radar altimetry yields meaningful results. An inter-comparison of the snow depth data sets emphasizes the limited usefulness of Warren climatology snow depth for freeboard-to-thickness conversion under current Arctic Ocean conditions reported in other studies. We test different freeboard-to-thickness and freeboard-to-draft conversion approaches. The mean observed ULS sea ice draft agrees with the mean sea ice draft derived from radar altimetry within the uncertainty bounds of the data sets involved. However, none of the approaches are able to reproduce the seasonal cycle in sea ice draft observed by moored ULS. A sensitivity analysis of the freeboard-to-thickness conversion suggests that sea ice density is as important as snow depth.

Sea ice thickness is one of the most sensitive variables in the Arctic climate system. In order to quantify changes in sea ice thickness, CryoSat-2 was launched in 2010 carrying a Ku-band radar altimeter (SIRAL) designed to measure sea ice freeboard with a few centimeters accuracy. The instrument uses the synthetic aperture radar technique providing signals with a resolution of about 300 m along track. In this study, airborne Ku-band radar altimeter data over different sea ice types have been analyzed. A set of parameters has been defined to characterize the differences in strength and width of the returned power waveforms. With a Bayesian-based method, it is possible to classify about 80% of the waveforms from three parameters: maximum of the returned power waveform, the trailing edge width and pulse peakiness. Furthermore, the maximum of the power waveform can be used to reduce the number of false detections of leads, compared to the widely used pulse peakiness parameter. For the pulse peakiness the false classification rate is 12.6% while for the power maximum it is reduced to 6.5%. The ability to distinguish between different ice types and leads allows us to improve the freeboard retrieval and the conversion from freeboard into sea ice thickness, where surface type dependent values for the sea ice density and snow load can be used.

Arctic sea ice area has been decreasing for the past two decades. Apart from melting, the southward drift through Fram Strait is the main ice loss mechanism. We present high resolution sea ice drift data across 79° N from 2004 to 2010. Ice drift has been derived from radar satellite data and corresponds well with variability in local geostrophic wind. The underlying East Greenland current contributes with a constant southward speed close to 5 cm s−1, and drives around a third of the ice export. We use geostrophic winds derived from reanalysis data to calculate the Fram Strait ice area export back to 1957, finding that the sea ice area export recently is about 25% larger than during the 1960's. The increase in ice export occurred mostly during winter and is directly connected to higher southward ice drift velocities, due to stronger geostrophic winds. The increase in ice drift is large enough to counteract a decrease in ice concentration of the exported sea ice. Using storm tracking we link changes in geostrophic winds to more intense Nordic Sea low pressure systems. Annual sea ice area export likely has a significant influence on the summer sea ice variability and we find low values in the 1960's, the late 1980's and 1990's, and particularly high values during 2005–2008. The study highlights the possible role of variability in ice export as an explanatory factor for understanding the dramatic loss of Arctic sea ice during the last decades.

Retrieval of Arctic sea ice thickness from CryoSat-2 radar altimeter freeboard data requires observational data to verify the relation between these two variables. In this study in-situ ice and snow data from 689 observation sites, obtained during the Sever expeditions in the 1980s, have been used to establish an empirical relation between thickness and freeboard of FY ice in late winter. Estimates of mean and variability of snow depth, snow density and ice density were produced on the basis of many field observations. These estimates have been used in the hydrostatic equilibrium equation to retrieve ice thickness as a function of ice freeboard, snow depth and snow/ice density. The accuracy of the ice thickness retrieval has been calculated from the estimated variability in ice and snow parameters and error of ice freeboard measurements. It is found that uncertainties of ice density and freeboard are the major sources of error in ice thickness calculation. For FY ice, retrieval of ≈ 1.0 m (2.0 m) thickness has an uncertainty of 46% (37%), and for MY ice, retrieval of 2.4 m (3.0 m) thickness has an uncertainty of 20% (18%), assuming that the freeboard error is ± 0.03 m for both ice types. For MY ice the main uncertainty is ice density error, since the freeboard error is relatively smaller than that for FY ice. If the freeboard error can be reduced to 0.01 m by averaging measurements from CryoSat-2, the error in thickness retrieval is reduced to about 32% for a 1.0 m thick FY floe and to about 18% for a 2.4 m thick MY floe. The remaining error is dominated by uncertainty in ice density. Provision of improved ice density data is therefore important for accurate retrieval of ice thickness from CryoSat-2 data.