For large-scale and long-term Arctic climate simulations appropriate parameterization of the surface albedo is required.
Therefore, the sea ice surface (SIS) albedo parameterization of the coupled ...regional climate model HIRHAM–NAOSIM was examined against broadband surface albedo measurements performed during the joint ACLOUD (Arctic CLoud Observations Using airborne measurements during polar Day) and PASCAL (Physical feedbacks of Arctic boundary layer, Sea ice, Cloud and AerosoL) campaigns, which were performed in May–June 2017 north of Svalbard. The SIS albedo parameterization was tested using measured quantities of the prognostic variables surface temperature and snow depth to calculate the surface albedo and the individual fractions of the ice surface subtypes (snow-covered ice, bare ice, and melt ponds) derived from digital camera images taken on board the Polar 5 and 6 aircraft.
The selected low-altitude (less than 100 m) flight sections of overall 12 flights were performed over surfaces dominated by snow-covered ice. It was found that the range of parameterized SIS albedo for individual days is smaller than that of the measurements. This was attributed to the biased functional dependence of the SIS albedo parameterization on temperature. Furthermore, a time-variable bias was observed with higher values compared to the modeled SIS albedo (0.88 compared to 0.84 for 29 May 2017) in the beginning of the campaign, and an opposite trend towards the end of the campaign (0.67 versus 0.83 for 25 June 2017). Furthermore, the surface type fraction parameterization was tested against the camera image product, which revealed an agreement within 1 %.
An adjustment of the variables, defining the parameterized SIS albedo, and additionally accounting for the cloud cover could reduce the root-mean-squared error from 0.14 to 0.04 for cloud free/broken cloud situations and from 0.06 to 0.05 for overcast conditions.
Melt ponds forming on Arctic sea ice in summer significantly reduce the surface albedo and impact the heat and mass balance of the sea ice. Therefore, their areal coverage, which can undergo rapid ...change, is crucial to monitor. We present a revised method to extract melt pond fraction (MPF) from Sentinel‐2 satellite imagery, which is evaluated by MPF products from higher‐resolution satellite and helicopter‐borne imagery. The analysis of melt pond evolution during the MOSAiC campaign in summer 2020, shows a split of the Central Observatory (CO) into a level ice and a highly deformed ice part, the latter of which exhibits exceptional early melt pond formation compared to the vicinity. Average CO MPFs are 17% before and 23% after the major drainage. Arctic‐wide analysis of MPF for years 2017–2021 shows a consistent seasonal cycle in all regions and years.
Plain Language Summary
In the Arctic summer, puddles of surface melt water, called melt ponds, form on the sea ice. These melt ponds reduce the ability of the surface to reflect the sunlight. Instead, they absorb more solar energy and pave the way into the ocean beneath where the energy is also absorbed. Thus, it is important to know where these melt ponds develop and what fraction of the surface they cover. To investigate this, we present a classification algorithm that is used to extract the areal fraction of melt ponds from satellite measurements. The special focus of this study is the MOSAiC campaign in summer 2020, where the research vessel Polarstern drifted with an ice floe for 1 year. We can see a separation of this floe into two parts. One of them shows melt pond formation much earlier than the other. This is because of different ice age and surface properties. Additionally, we use the classification algorithm to analyze the differences of melt pond fraction between different dates and regions in the Arctic.
Key Points
Algorithm to extract melt pond and open water areas from Sentinel‐2 imagery with maximum uncertainty of 6%
Exceptional early melt pond formation on MOSAiC Central Observatory, summer 2020, compared to broader vicinity
High spatial and temporal variability of melt pond fraction on local and regional scales
Surface reflectance (SR) estimation is the most critical preprocessing step for deriving geophysical parameters in multi-sensor remote sensing. Most state-of-the-art SR estimation methods, such as ...the vector version of the Second Simulation of the Satellite Signal in the Solar Spectrum (6SV) radiative transfer (RT) model, depend on accurate information on aerosol and atmospheric gases. In this study, a Simplified and Robust Surface Reflectance Estimation Method (SREM) based on the equations from 6SV RT model, without integrating information of aerosol particles and atmospheric gasses, is proposed and tested using Landsat 5 Thematic Mapper (TM), Landsat 7 Enhanced Thematic Mapper plus (ETM+), and Landsat 8 Operational Land Imager (OLI) data from 2000 to 2018. For evaluation purposes, (i) the SREM SR retrievals are validated against in situ SR measurements collected by Analytical Spectral Devices (ASD) from the South Dakota State University (SDSU) site, USA; (ii) cross-comparison between the SREM and Landsat spectral SR products, i.e., Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) and Landsat 8 Surface Reflectance Code (LaSRC), are conducted over 11 urban (2013–2018), 13 vegetated (2013–2018), and 11 desert/arid (2000 to 2018) sites located over different climatic zones at a global scale; (iii) the performance of the SREM spectral SR retrievals for low to high aerosol loadings is evaluated; (iv) spatio-temporal cross-comparison is conducted for six Landsat paths/rows located in Asia, Africa, Europe, and the United States of America from 2013 to 2018 to consider a large variety of land surfaces and atmospheric conditions; (v) cross-comparison is also performed for the Normalized Difference Vegetation Index (NDVI), the Enhanced Vegetation Index (EVI), and the Soil Adjusted Vegetation Index (SAVI) calculated from both the SREM and Landsat SR data; (vi) the SREM is also applied to the Sentinel-2A and Moderate Resolution Imaging Spectrometer (MODIS) data to explore its applicability; and (vii) errors in the SR retrievals are reported using the mean bias error (MBE), root mean squared deviation (RMSD), and mean systematic error (MSE). Results depict significant and strong positive Pearson’s correlation (r), small MBE, RMSD, and MSE for each spectral band against in situ ASD data and Landsat (LEDAPS and LaSRC) SR products. Consistency in SREM performance against Sentinel-2A (r = 0.994, MBE = −0.009, and RMSD = 0.014) and MODIS (r = 0.925, MBE = 0.007, and RMSD = 0.014) data suggests that SREM can be applied to other multispectral satellites data. Overall, the findings demonstrate the potential and promise of SREM for use over diverse surfaces and under varying atmospheric conditions using multi-sensor data on a global scale.
Arctic summer sea ice experiences rapid changes in its sea-ice concentration, surface albedo, and the melt pond fraction. This affects the energy balance of the region and demands an accurate ...knowledge of those surface characteristics in climate models. In this paper, the broadband albedo (300–3000 nm) of Arctic sea ice is derived from MEdium Resolution Imaging Spectrometer (MERIS) optical swath data by transforming the spectral albedo as an output from the Melt Pond Detector (MPD) algorithm with a newly developed spectral-to-broadband conversion (STBC). The new STBC replaces the previously applied spectral averaging method to provide a more accurate broadband albedo product, which approaches the accuracy of 0.02–0.05 required in climate simulations and allows a direct comparison to broadband albedo values from climate models. The STBC is derived empirically from spectral and broadband albedo measurements over landfast ice. It is validated on a variety of simultaneous spectral and broadband field measurements over Arctic sea ice, is compared to existing conversion techniques, and performs better than the currently published algorithms. The root-mean-square deviation (RMSD) between broadband albedo that was measured and converted by the STBC is 0.02. Other conversion techniques, the spectral averaging method and the linear combination of albedo values from four MERIS channels, result in higher RMSDs of 0.09 and 0.05, respectively. The improved MERIS-derived broadband albedo values are validated with airborne measurements. Results show a smaller RMSD of 0.04 for landfast ice than the RMSD of 0.07 for drifting ice. The MERIS-derived broadband albedo is compared to broadband albedo from ERA5 reanalysis to examine the albedo parameterization used in ERA5. Both albedo products agree over large scales and in temporal patterns. However, consistency in point-to-point comparison is rather poor, with differences up to 0.20, correlations between 0.69 and 0.79, and RMSDs in excess of 0.10. Differences in sea-ice concentration and cloud-masking uncertainties play a role, but most discrepancies can be attributed to climatological sea-ice albedo values used in ERA5. They are not adequate and need revising, in order to better simulate surface heat fluxes in the Arctic. The advantage of the resulting broadband albedo data set from MERIS over other published data sets is the accompanied data set of available melt pond fraction. Melt ponds are the main reason for the sea-ice albedo change in summer but are currently not represented in climate models and atmospheric reanalysis. Additional information about melt evolution, together with accurate albedo retrievals, can aid the challenging representation of sea-ice optical properties in those models in summer.
The directional reflection of solar radiation by the Arctic Ocean is mainly shaped by two dominating surface types: sea ice (often snow-covered) and open ocean (ice-free). In the transitional zone ...between them, the marginal sea ice zone (MIZ), the surface reflection properties are determined by a mixture of the reflectance of both surface types. Retrieval methods applied over the MIZ need to take into account the mixed directional reflectivity; otherwise uncertainties in the retrieved atmospheric parameters over the MIZ may occur. To quantify these uncertainties, respective measurements of reflection properties of the MIZ are needed. Therefore, in this case study, an averaged hemispherical–directional reflectance factor (HDRF) of the inhomogeneous surface (mixture of sea ice and open ocean) in the MIZ is derived using airborne measurements collected with a digital fish-eye camera during a 20 min low-level flight leg in cloud-free conditions. For this purpose, a sea ice mask was developed to separate the reflectivity measurements from sea ice and open ocean and to derive separate HDRFs of the individual surface types. The respective results were compared with simulations and independent measurements available from the literature.
It is shown that the open-ocean HDRF in the MIZ differs from homogeneous ocean surfaces due to wave attenuation. Using individual HDRFs of both surface types and the sea ice fraction, the mixed HDRF describing the directional reflectivity of the inhomogeneous surface of the MIZ was retrieved by a linear weighting procedure. Accounting for the wave attenuation, good agreement between the average measured HDRF and the constructed HDRF of the MIZ was found for the presented case study.
To evaluate the performance of the eXtensible Bremen Aerosol/cloud and surfacE
parameters Retrieval (XBAER) algorithm, presented in the Part 1 companion paper to this paper, we apply the XBAER ...algorithm to the Sea and Land
Surface Temperature Radiometer (SLSTR) instrument on board Sentinel-3. Snow
properties – snow grain size (SGS), snow particle shape (SPS) and specific
surface area (SSA) – are derived under cloud-free conditions. XBAER-derived
snow properties are compared to other existing satellite products and
validated by ground-based and aircraft measurements. The atmospheric
correction is performed on SLSTR for cloud-free scenarios using Modern-Era
Retrospective Analysis for Research and Applications (MERRA) aerosol optical
thickness (AOT) and the aerosol typing strategy according to the standard XBAER
algorithm. The optimal SGS and SPS are estimated iteratively utilizing a
look-up-table (LUT) approach, minimizing the difference between
SLSTR-observed and SCIATRAN-simulated surface directional reflectances at
0.55 and 1.6 µm. The SSA is derived for a retrieved SGS and SPS pair.
XBAER-derived SGS, SPS and SSA have been validated using in situ measurements from
the recent campaign SnowEx17 during February 2017. The comparison shows a
relative difference between the XBAER-derived SGS and SnowEx17-measured SGS of
less than 4 %. The difference between the XBAER-derived SSA and SnowEx17-measured SSA is 2.7 m2/kg. XBAER-derived SPS can be
reasonably explained by the SnowEx17-observed snow particle shapes. Intensive validation shows that (1) for SGS and SSA, XBAER-derived results
show high correlation with field-based measurements, with correlation
coefficients higher than 0.85. The root mean square errors (RMSEs) of SGS and
SSA are around 12 µm and 6 m2/kg. (2) For SPS, aggregate SPS
retrieved by XBAER algorithm is likely to be matched with rounded grains
while single SPS in XBAER is possibly linked to faceted crystals. The comparison with aircraft measurements, during the Polar Airborne
Measurements and Arctic Regional Climate Model Simulation Project
(PAMARCMiP) campaign held in March 2018, also shows good agreement (with
R=0.82 and R=0.81 for SGS and SSA, respectively). XBAER-derived SGS and
SSA reveal the variability in the aircraft track of the PAMARCMiP campaign. The
comparison between XBAER-derived SGS results and the Moderate Resolution Imaging Spectroradiometer (MODIS) Snow-Covered Area and
Grain size (MODSCAG) product over Greenland shows similar spatial
distributions. The geographic distribution of XBAER-derived SPS over
Greenland and the whole Arctic can be reasonably explained by campaign-based
and laboratory investigations, indicating a reasonable retrieval accuracy of
the retrieved SPS. The geographic variabilities in XBAER-derived SGS and SSA
both over Greenland and Arctic-wide agree with the snow metamorphism
process.