Next Article in Journal
Pulsed Orthogonal Time Frequency Space: A Fast Acquisition and High-Precision Measurement Signal for Low Earth Orbit Position, Navigation, and Timing
Previous Article in Journal
The Influence of Short-Term Water Level Fluctuations on the Habitat Response and Ecological Fragility of Siberian Cranes in Poyang Lake, China
Previous Article in Special Issue
Sea Ice Detection from GNSS-R Data Based on Local Linear Embedding
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Microwave Radiometer Sea Ice Concentration Products over the Baltic Sea

1
Finnish Meteorological Institute, FI-00101 Helsinki, Finland
2
Integrated Climate Data Center (ICDC), Center for Earth System Research and Sustainability (CEN), University of Hamburg, 20144 Hamburg, Germany
3
National Space Institute, Technical University of Denmark, 2800 Lyngby, Denmark
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(23), 4430; https://doi.org/10.3390/rs16234430
Submission received: 18 October 2024 / Revised: 14 November 2024 / Accepted: 25 November 2024 / Published: 27 November 2024

Abstract

:
Sea ice concentration (SIC) monitoring in the Arctic using microwave radiometer data is a well-established method with numerous published accuracy studies. For the Baltic Sea, accuracy studies have not yet been conducted. In this study, we evaluated five different SIC products over the Baltic Sea using MODIS (250 m) and Sentinel-2 (10 m) open water–sea ice classification charts. The selected SIC products represented different SIC algorithm types, e.g., climate data records and near-real-time products. The one-to-one linear agreement between the radiometer SIC dataset and the MODIS/Sentinel-2 SIC was always quite poor; the slope of the linear regression was from 0.40 to 0.77 and the coefficient of determination was from 0.26 to 0.80. The standard deviation of the difference was large and varied from 15.5% to 26.8%. A common feature was the typical underestimation of the MODIS/Sentinel-2 SIC at large SIC values (SIC > 60%) and overestimation at small SIC values (SIC < 40%). None of the SIC products performed well over the Baltic Sea ice, and they should be used with care in Baltic Sea ice monitoring and studies.

1. Introduction

In this study, we evaluated selected sea ice concentration (SIC) products derived from satellite microwave radiometer data over the Baltic Sea. This evaluation was conducted using MODIS 250 m pixel size and Sentinel-2 (S-2) 10 m pixel size open water–sea ice (OWSI) charts.
Over the Arctic, numerous SIC product validation studies have been conducted, e.g., [1,2,3,4,5], but none so far have been over the Baltic Sea. The usage of SIC products over the Baltic Sea has been limited; no SIC interannual trend studies have been published. This is due to the following reasons. First, over the Baltic Sea, radiometer SIC estimation is challenging due to its typical rugged coastline, with many islands which manifest as a land spillover effect on the measured brightness temperatures and further as erroneous SIC values if not properly corrected or filtered [6]. Global SIC estimation algorithms are tuned with tie points for the Arctic sea ice, and it is poorly quantified how they fit for low-salinity (0.2–2 psu) and thin (<0.8 m) Baltic Sea ice [7]. Another error source is variable weather in the Baltic Sea, with mid-winter air temperatures sometimes well above 0 °C making snow cover on the sea ice wet, which accelerates snow metamorphosis and creates surface and internal layering in the snow, and, therefore, a SIC algorithm tuned to the Arctic’s cold conditions may give false SIC values. Finally, the Baltic Sea is quite small compared to the resolution of some SIC products, e.g., those with a 25 km resolution. It is envisaged that validation data on global SIC products will show their usefulness for Baltic Sea ice monitoring.
The validation of SIC products follows methods used in [5] where Landsat 30 m imagery was used as reference data. Landsat OWSI charts were aggregated to SIC charts with the pixel sizes of radiometer SIC products by calculating the fraction of sea ice pixels to the total number of pixels in each radiometer grid cell. Next, various statistical parameters between the Landsat and radiometer SIC charts were calculated. Validations were also conducted separately for first-year ice (FYI)-, mixed FYI and multiyear ice (MYI), or MYI-dominated cases as well as for different ice regimes, e.g., ice-edge, freeze-up and melt conditions.
The validation study with the Landsat imagery in [5] had ten different SIC products including OSI SAF [8], ESA SICCI (Sea Ice Climate Change Initiative) [8,9], and NOAA [10,11] climate data records (CDR). A superior linear agreement between the radiometer SIC and Landsat SIC was found for the 25 and 50 km SICCI-2 products. We do not aim to conduct as comprehensive a validation study as that in [5]. We mainly include SIC products from the EUMETSAT OSI SAF (Ocean and Sea Ice Satellite Application Facility) and ESA SICCI (Sea Ice Climate Change Initiative) which represent both climate data records (CDR) and near-real-time products and also high–low resolution hybrid products by pansharpening.
In the following, Section 2 first briefly describes our study area, the Baltic Sea. Next, the chosen SIC products are introduced, followed by description of the MODIS and S-2 OWSI charts and methods used in their calculation. Section 4 presents our validation studies and results, and Section 5 presents our discussion and conclusions.

2. Study Area

The Baltic Sea is a semi-enclosed brackish sea water basin in northern Europe. The Baltic Sea ice cover usually begins to form in November and has its largest extent between January and March [12,13,14]. The normal ice break-up starts in April and the ice melts completely by the beginning of June. The maximum annual ice cover ranges from 9% to 100% of the whole Baltic Sea area, and the average is 50% [14,15]. The Baltic Sea ice extent and duration of the ice season depend on the indices of the North Atlantic Oscillation (NOA) and Arctic Oscillation (AO); for example, the annual maximum ice extent generally decreases with the increasing indices of the AO and NAO [13]. The ice in the Baltic Sea is categorized as drift ice and landfast ice (LFI). LFI occurs in the coastal and archipelago areas and usually extends to the 5–15 m isobath. Drift ice has a dynamic nature due to being forced to move by winds and currents. The motion of drift ice results in an uneven and broken ice field with distinct floes up to several kilometers in diameter, leads and cracks, brash ice barriers, rafted ice, and ice ridges. In the Bay of Bothnia, the annual maximum ice thickness is typically 0.65–0.80 m, and it reaches 0.3–0.5 m even in mild winters [12,13]. The measured all-time maximum is around 1.2 m. In the southern Baltic Sea, in the coastal areas of Germany and Poland and the Danish straits, the annual maximum ice thickness seldom exceeds 0.5 m [16]. The thickness of ice ridges (sail height plus keel depth) is typically 5 to 15 m [17]. The salinity of the Baltic Sea ice is typically from 0.2 to 2‰ depending on the water salinity, growth rate, and ice age [7].

3. Materials and Methods

3.1. Radiometer SIC Products

Five radiometer SIC products selected for this study are summarized in Table 1. Instead of including all available SIC products in our study, we selected products representing different SIC algorithm types. One of the products, the ARTIST (Arctic Radiation and Turbulence Interaction Study) Sea Ice (ASI) SIC algorithm using AMSR2 (ASI-AMSR2) 89 GHz data [18], represents high-resolution SIC products retrieved with 90 GHz brightness temperature ( T B ) channels in combination with low-resolution weather filters. The OSI-408-a is near-real-time SIC product based on the AMSR2 18.7 and 36.5 GHz T B data used in many SIC algorithms [19]. The OSI-450-a and OSI-458 are climate data records based on SSMIS and AMSR2 data, respectively [8,20]. The SICCI-HR-SIC (HR means high resolution) is a resolution-enhancing SSMIS product where SIC with 19.35 and 37.0 GHz T B data (the same SIC as in the OSI-450-a product) is pansharpened using a SIC with finer-resolution 91.7 GHz data [21]. These products have pixel sizes from 3.125 to 25 km. Some of the products are oversampled, like the ASI-AMSR2 3.125 km product. Products with larger pixel sizes are taken to be too coarse for the Baltic Sea, especially for the Bay of Bothnia and Gulf of Finland where sea ice forms every winter.
Descriptions of the SIC retrieval algorithms, e.g., tie point selection and applied weather filters, can be found in the references shown in Table 1. All products, except the ASI-AMSR2, use the daily dynamic tie points method from Maass and Kaleschke [6] to correct the land spillover already occurring at the T B level, a dynamic open water filter, and the atmospheric correction of T B s [8]. The ASI-AMSR2 is based on fixed tie points and fixed open water filters [18], and T B s are not corrected for atmospheric effects. All SIC products have SIC values between 0 and 100%, i.e., they are truncated to the range between 0% and 100%.

3.2. MODIS OWSI Chart

For the validation of the radiometer SIC products over the Baltic Sea SIC, we derived an open water—sea ice (OWSI) chart using MODIS band 1 reflectance data which have a 250 m spatial resolution. The total number of the MODIS OWSI charts is 62, and they cover February to April 2019. This OWSI chart derivation follows our earlier study for the Barents and Kara Seas [25]. In the following, the used MODIS datasets and their processing are first described, including the automatic and manual cloud masking of the MODIS reflectance data. Next, the OWSI chart derivation is described.
The MODIS spectrometer measures reflected solar and emitted thermal radiation in 36 spectral bands ranging from the visible to the thermal infrared (TIR) (0.415 to 14.235 μm) [26]. The MODIS instrument is aboard NASA Terra and Aqua satellites. In this study, we derived an OWSI chart for the Baltic Sea using band 1 reflectance ( R 1 ) data at a 250 m spatial resolution and the MODIS cloud mask product [27]. MODIS sea ice surface temperature (IST) data and an RGB image with a 2-1-1 band combination (RGB211) were used in finding cloud-free MODIS swath datasets. The following Terra and Aqua MODIS products were used: the MOD/MYD/02QKM-MODIS/Terra Calibrated Radiances 5-Min L1B Swath 250 m, MOD/MYD/03-MODIS Geolocation Fields 5-Min L1A Swath 1 km, MOD/MYD/35_L2-MODIS Cloud Mask and Spectral Test Results 5-Min L2 Swath 250 m and 1 km, and MOD/MYD/29-MODIS Sea Ice Extent 5-Min L2 Swath 1 km. The MOD/MYD29 product contained IST data, and the cloud mask was derived from the MOD/MYD35_L2 unobstructed field-of-view flag [28]. In the following, these products are only referred to with the ‘MOD’ prefix.
The MODIS products were rectified to the Mercator projection used for the Finnish Ice Service (FIS) chart using the MODIS Swath Reprojection Tool developed by the NASA LP DAAC (Land Processes Distributed Active Archive Center). In the rectification of the MODIS products, original pixel sizes (250 or 1000 m) were kept, and fixed Mercator grids were used. The size of the rectified data was 1330 km in the northing and 1120 km in the easting. The MODIS sensor scan angle ( θ s ) was calculated from the sensor zenith angle, Earth radius, and satellite orbit altitude. Raw data from the bands 1 and 2 were converted to top-of-the-atmosphere (TOA) sun angle corrected reflectances ( R 1 and R 2 ). The land mask for the Baltic Sea study area was retrieved from the MOD44W (Land Water Mask Derived L3 Global 250 m) product.
In all MODIS data analyses herein, the retrieval of the MODIS OWSI chart θ s was limited to a maximum of 40° (the overall maximum is 55°). At a θ s of 40° the MODIS spatial resolution in the across-track direction decreased by a factor of around two, e.g., from 250 m at the nadir to ~500 m at 40° [29]. In the MOD29 and MOD35 products, daylight data were defined as data collected with a θ z s less than 85° [27,28]. We used here a somewhat smaller limit of 80° as at high sun zenith angles the 1 / c o s ( θ z s ) term in the TOA reflectance calculation is highly variable, and, thus, the accuracy of the estimated θ z s is critical. This θ z s limitation excluded all the MOD02 data from December 2018 to January 2019 in the Bay of Bothnia. Therefore, we used here only MODIS data acquired over the Baltic Sea from 1 February 2019 onwards to end of April 2019.
The following procedure was used to select mostly cloud-free MOD02 data over the sea ice-covered part of the Baltic Sea. First, all available MODIS IST data were processed over the Baltic Sea from 1 February to 30 April 2019. Visual analysis of the IST charts with the MOD35 cloud mask was conducted to select times for the MOD02 data, which possibly had large cloud-free areas over the sea ice-covered part of the Baltic Sea. Next, RGB211 images with and without s cloud mask were used in the further screening of the MOD02 datasets for the OWSI chart derivation. We ended up with 20 MOD02 swath datasets for February, 21 for March, and 21 for April. Here, a MOD02 dataset was either one MOD02 swath or two swaths stitched together. For some days, there were two datasets (Terra and Aqua) and their temporal difference was mostly 10 or 15 min and at maximum 100 min.
Cloud masking was started by designating pixels that had ‘confident clear’, ‘probably clear’, or ‘probably cloudy’ values in the unobstructed field-of-view quality flag of the MOD35 product as cloud-free. This logic followed the MOD29 IST product generation with the target to increase the number of retrievals balanced against the cloud conservative nature of the cloud mask [28]. A visual analysis of the resulting cloud mask over the RGB211 image showed that sometimes cloud-free pixels of open water or the thin ice within led to larger open water/thin ice areas being classified as cloudy, and cloud-free sea ice, in general, was classified as cloudy. This is because NSIDC’s Near-real-time Ice and Snow Extent (NISE) product with a 25 km resolution was used to initialize the sea ice background flag for directing the cloud masking algorithm processing flow [28]. Open water or thin ice when the NISE sea ice flag is set may be classified as cloudy, and sea ice when the flag is not set may also be classified as cloudy. Methods have been developed for the correction of these cloud mask-induced errors in the OWSI classification, e.g., in [30], but we decided to correct these errors by the manual editing of the cloud mask. For the sea ice-covered Baltic Sea, which has a rather small area, the manual cloud mask editing for one MODIS dataset was not an overly laborious task, and it also served as a quality control step in the OWSI chart processing.
Before the manual cloud masking, all single or groups of two cloudy pixels within the NISE sea ice flag set area were first removed. These were assumed to be erroneous cloud detections over the open water or thin ice. Next, the cloud mask was aggregated into 5 km blocks (5 by 5 1 km pixels), and if there were more than or equal to three cloudy pixels (12% of all 25 pixels), then the whole block was flagged as cloudy. Using 5 by 5 km blocks for the cloud mask, we could better identify large cloud-free areas and discard areas where the cloud mask was tattered. Small cloud-free areas (‘holes’ in the mask) less than ten blocks in size were flagged as cloudy. Finally, the following manual editing procedures could be conducted: filling holes, removing erroneous cloud mask elements, and masking arbitrary polygonal areas as cloudy or clear. The manual editing was conducted using the RGB211 image. An example of the original and manually edited cloud masks is shown in Figure 1.
We aimed to develop a simple OWSI classification method based only on the threshold R 1 (bandwidth 620–670 nm, red). We only used the R 1 data as it had the best possible MODIS resolution, 250 m. The R 2 data (841–876 nm, near-infrared (NIR)) also had a 250 m resolution, but we evaluated that it was not needed in the OWSI classification. In addition, we have previously determined an R 1 -based OWSI classification for the Barents and Kara Seas [25]. In the Baltic Sea, MODIS R 1 data have been used to map the sea ice extent in the Gulf of Riga [31]. In this study, the threshold for the OWSI classification was determined automatically for each swath dataset based on a sampled R 1 histogram which was bimodal, showing OW and SI peaks. More accurate OWSI classification is also possible using other MODIS reflectance bands and shortwave and thermal infrared bands, e.g., as in [28,30,32], but this is at the cost of decreasing the resolution to 500 m or 1 km. In the MOD29 product, for daytime, there is sea ice extent data based on various reflectance tests, including the threshold R 1 , at 1 km resolution [28]. An algorithm called IceMap500 has been developed to generate 500 m sea ice extent maps with a systematic accuracy over 90% [30]. Machine learning methods have also been developed to classify OW, SI, and clouds; for example, Jiang et al. [33] constructed a training sample library and used a Multi-Feature Level Fusion Random Forest classification algorithm that integrated multiple features from the reflectance bands 1, 2, and 7. Their SI chart had a 250 m pixel size, but this required the interpolation of the band 7 data at a 500 m resolution to 250 m.
For the determination of the R 1 threshold for the OWSI classification in the Baltic Sea, we manually selected samples in the cloud-masked RGB211 images from February 2019 for OW, thin ice (THI), and snow-covered sea ice (SSI). These samples represented winter conditions in the Baltic Sea. The size of the sample windows for the OW and SSI was 5.25 by 5.25 km (21 by 21 pixels), and it was 2.75 by 2.75 km (11 by 11 pixels) for the THI. The number of sample windows was 210 for the OW, 278 for the THI, and 128 for the SSI. Another set of samples were selected in April 2019, which represented springtime melting sea ice conditions. Here, we selected windows for the OW (219 windows), SSI (35 windows), and melting sea ice (likely mostly snow-free) (MSI, 385 windows, 11 by 11 pixels window size).
The resulting R 1 probability density function (pdf) for the OW is very narrow (see Figure 2); for example, the 95th percentile of the R 1 is only 0.070 in the winter data. The SSI is well separated from the OW, and the 5th percentile of the R 1 is around 0.55. In the winter data, the THI has a large variation, from 0.16 to 0.51 by the 5th and 95th percentiles. The MSI likewise has a large variation, from 0.073 to 0.42. Setting the R 1 threshold to 0.10, only 0.12% of the OW samples in the winter data are classified as THI, and only 1.0% of the THI is classified as OW. In the spring data, the OW and MSI misclassifications are 0.46% and 4.5%, respectively. Unfortunately, it is not possible to separate fully the OW and very thin ice which both have small R 1 values. In addition, the R 1 data (TOA) here do not have atmospheric corrections, and in the springtime, water turbity close to the coast and river estuaries may raise R 1 . Therefore, it is better to use somewhat too high than too small of a threshold. We decided to use the threshold of 0.10 for all R 1 data as it gives very small OW misclassifications in our manually selected sample data:
sea   ice   if   R 1 > 0.10 .
This threshold was manually determined. In the resulting OWSI charts, very thin ice (thickness: a few cm) and very wet sea ice may be classified as OW. Fog and thin unmasked clouds over OW can also be misclassified as sea ice. In addition, the mixed-pixel effect affects the OWSI classification as even a small fraction of sea ice, especially that of snow-covered ice, within the 250 m pixel leads to assigning the SI class to the pixel. This would further lead to the overestimation of SIC when the MODIS data are used to calculate SIC at radiometer data resolutions. An example of the resulting OWSI charts in shown in Figure 3 for the same date and time as the RGB211 images in Figure 1.

3.3. Sentinel-2 OWSI Chart

Sentinel-2A and Sentinel-2B (S-2A/B) carry the MultiSpectral Instrument (MSI) with 13 different spectral bands from VNIR (visible–near-infrared) to SWIR (shortwave infrared) with band central wavelengths varying from 443 to 2202 (S-2A) / 2186 (S-2B) nm [34,35,36]. The spatial resolution for different bands is either 10, 20, or 60 m. The swath width of MSI is 290 km. The Sentinel-2 mission supports Copernicus land, marine, security, emergency, and climate change services and studies, including the monitoring of vegetation, soil, and water cover, as well as the observation of inland waterways and coastal areas. The MSI data can also be used for the remote sensing of sea ice, e.g., for open water–sea ice (OWSI) classification, but MSI lacks longwave infrared bands for good cloud masking over sea ice, e.g., 3.9 and 11 µm TIR bands, as in MODIS.
There are two end-user MSI products: Level-1C top-of-atmosphere (TOA) reflectance and Level-2A bottom-of-atmosphere (BOA) reflectance. We used here the L1C product with the processing baseline 02.07 version. The L1C product had 100 by 100 km fixed tiles (or granules) in the UTM/WGS84 projection. The L1C product had opaque (dense clouds) and cirrus cloud masks [37]. The cloud masks were processed with data sampled at a 60 m spatial resolution for all spectral bands. The land masking of the L1C tiles was conducted here with the ESA SNAP software (version 9.0) (fractional land/water mask). The fixed land mask was based on the ESA Global Land Cover Map data. It did not identify all the small islands in the Bay of Bothnia.
We aimed to use the 10 m resolution bands (best possible resolution)—B2 (490 nm; blue), B3 (560 nm; green), B4 (665 nm; red) and B8 (842 nm; NIR)—from the L1C product for the OWSI classification over the Bay of Bothnia. B4 was equal to MODIS band 1 (bandwidth of 620–670 nm) which was used here for the MODIS OWSI classification.
As S-2 cloud detection over sea ice may be inaccurate, we searched S-2 data over the Bay of Bothnia only for the dates when the MODIS data were mostly cloud-free and the sun zenith angle was mostly below 80°. The Bay of Bothnia was covered with five 100 km tiles in the UTM34 projection. Our data search for February-April 2019 resulted in S-2 imagery for 13 dates; see an example in Figure 4.
S-2 data have been used previously in sea ice remote sensing studies either as a validation dataset or to classify sea ice for further analyses, like for lead statistics, e.g., in [38,39,40]. Ludwig et al. [38] evaluated the accuracy of their MODIS and AMSR-E data-based SIC product with a 1 km pixel using 79 cloud-free (visually determined) S-2 L1C products acquired mainly over the Arctic FYI. They investigated the R s of the visible spectrum bands 2–4 for open water (OW), thin ice, and thick ice classification. These bands did not show significant differences in the classification, and there was a clear distinction between the thin and thick ice and a clear distinction between the OW and sea ice (SI) if only thick ice was present. They decided to use only the band 4 reflectance ( R 4 ) for the surface type classification and selected specific classification thresholds (OW vs. thin ice and thin vs. thick ice) for each R 4 image instead of using global thresholds for all images. Muchow et al. [39] investigated lead-width distribution for the Antarctic sea ice using 20 cloud-free L1C products. Lead detection was conducted by thresholding R 4 . For the threshold determination, they selected manually sample data for OW and the following four ice types: nilas (NI; thickness < 10 cm), gray sea ice (GI, 10–15 cm), gray–white ice (GWI; 15–30 cm) and sea ice covered with snow. The sample data were used to calculate R 4 histograms for each surface type, and then a summation of Gaussian functions was fitted to them. The threshold for each surface type was then determined from the intersection of two curves adjacent to each other. For the lead identification, two different thresholds were used: 0.10 for leads covered with OW and 0.17 for leads covered with OW and NI. The threshold of 0.10 classified 29% of NI as OW and the threshold of 0.17 classified 11% of GWI as NI. Wang et al. [40] developed a decision tree classification of OW and various ice types in the Liaodong Bay of the Bohai Sea using S-2A/2B data. Using one S-2B image, they selected manually sampled areas for different surface types: OW, NI, GWI, and landfast ice (LFI). Utilizing mean reflectance spectra ( R 2 to R 12 ), a decision tree classification was determined for these surface types. This classification tree included 10 m and 20 m bands. The total accuracy of the classification was 88% according to visual validation. These previous sea ice classification studies demonstrate that OWSI classification in the Baltic Sea is possible with the S-2 data.
We started the development of the OWSI classification method by first investigating the pdfs of the 10 m reflectances of R 2 (490 nm; blue), R 3 (560 nm; green), R 4 (665 nm; red), and R 8 (842 nm; NIR)), as well as the normalized difference water index (NDWI) in the Bay of Bothnia. The NDWI is calculated as follows [41]:
N D W I = R 3 R 8 R 3 + R 8 .
The NDWI detects open water features against soil and terrestrial vegetation features (they have a smaller NDWI than OW). It is noted that the normalized difference snow index (NDSI) could not be calculated from the S-2 10 m reflectance data. The pdfs were calculated using imagery of one tile in the northern Bay of Bothnia (tile 34WFT) acquired on 5 March and on 6 April 2019. The first date represented winter conditions (the air temperature T a at the Kemi 1 lighthouse was −10 °C), and there were large areas of thin ice (<15 cm) according to the Finnish Ice Service (FIS) ice chart; the second date had spring, early melting conditions (the T a was around 0 °C) when there was no thin ice in the Bay of Bothnia. For 5 March, all four 10 m R pdfs show a distinct peak at small R values, representing OW, and another peak at high R values, representing snow-covered SI, mostly LFI; see Figure 5. There are a few smaller peaks, although they are not very distinct, representing different types of level, deformed, and snow-free thin ice. There are no clear R thresholds between the OW and thin ice. For R 4 and R 8 , the OW peak is somewhat narrower than for R 2 and R 3 . The R pdfs on 6 April also have a OW peak, but they have much a smaller snow-covered SI peak as large areas of LFI snow had melted; see Figure 5. The snow melting resulted in large fractions of R values in the 0.2–0.3 to 0.6 range. There is now a clear R local minimum after the OW peak, as there was no new thin ice on 6 April. The NDWI pdf shows a clear LFI peak and smaller OW peak but no better discrimination of surface types than the single R pdfs.
Likely, there was not much difference that affected which R band we would use for the OWSI classification, but R 4 and R 8 could have been better than R 2 and R 3 due to their narrower OW peaks. Therefore, we selected R 4 for further analyses as it has also been used previously for lead detection and SIC estimation [38,39].
First, we investigated the R 4 statistics for the OW on 5 March and 6 April by visually selecting 195 OW rectangles from four S-2 tiles acquired on 5 March and 194 rectangles from five tiles on 6 April. The size of the rectangles was from 0.15 to 870 km2. On 5 March, the tile-based mean R 4 varied only from 0.047 to 0.056, the mode was either 0.05 or 0.06 when the pdf bin width was 0.01, and the 99th percentile was from 0.057 to 0.065. On 6 April, the tile based R 4 statistics were slightly larger: the mean R 4 varied from 0.048 to 0.059, and the 99th percentile varied from 0.056 to 0.078. The mode was again either 0.05 or 0.06.
These OW R 4 statistics suggest that fixed thresholds for the OWSI classification are not feasible. A tile-based threshold is a better approach. This tile-based OW threshold takes into account variation in the sun zenith and sensor zenith angles (TOA reflectance does not correct for atmospheric effects, which depend on the sensor zenith angle, i.e., the path length). In the above rectangle selection, we observed that it was sometimes difficult to visually discriminate between OW and very thin ice, i.e., dark nilas and new ice (including grease ice and slush), which had same kind of ‘dark’ appearance (in level and texture) in the R 4 image. There was also sometimes thin fog over the OW and thin ice, which made their identification more difficult and raised the R 4 . Most notably, there were areas where both the OW and very thin ice had a similar R 4 level (and no texture). The OW and new ice also had very similar S-2 spectra in [40]. These observations show that full discrimination between OW and very thin ice is not possible, and a R 4 threshold for the OWSI classification will classify some unknown part of thin ice to the OW class. This is likely also the case for OW and melting wet ice in late spring.
Next, we needed to find a way to determine the R 4 threshold for each L1C tile. One way to conduct this is to manually sample OW data from a tile and to calculate the 99th (or a lower 95th) percentile of R 4 . Unfortunately, this manual sampling is time consuming. We tested fitting a mixture of Gaussian pdfs to the R 4 OW and SI data in a tile and used the intersection between the pdf with the smallest mean (representing OW) and the pdf with the second smallest mean as the threshold. The fitting seemed to work fine when there was a local minimum after the OW peak but not when this was not the case. Therefore, we decided to set the threshold manually with the help of the R 4 pdf. The manually set threshold varied from 0.06 to 0.11, and it was larger in winter, from 0.09 to 0.11, than in the melting season (12 April onwards), from 0.06 to 0.08. The manual setting was easier in the melting season as new thin ice was not present.
Before the OWSI classification, the opaque (dense clouds) and cirrus cloud masks were combined, and the combined mask was manually edited. This editing was conducted with the help of an RGB image with R 4 - R 3 - R 2 bands (approximating a true color image); see an example in Figure 4. In the manual cloud masking editing we rather masked too much than too little to cover for certain all cloudy and cloud shadow areas. A visual comparison of the cloud masks of the R 4 and RGB images showed that the opaque cloud mask was sometimes showing erroneously cloudy areas over LFI. In some R 4 images from later in April, there were increased R 4 values over open water at river mouths due to water turbidity. These areas were included in the manual cloud mask as they would otherwise be classified as SI.
An example of the resulting OWSI charts is shown in Figure 6. Here, the S-2 tile data are the same as in Figure 4. All tile-based OWSI charts for the same S-2 overpass (max. five tiles) were combined to create an OWSI mosaic over the Bay of Bothnia.

3.4. Coastal Weather Station Data

Air temperature data from three coastal weather stations in the Bay of Bothnia (see Figure 7) were used to identify cold winter conditions (air temperature below 0 °C) and melting conditions for the SIC product evaluations. There were cold conditions roughly on 1–15 February, 17–23 February, 28 February–15 March, and 25–28 March; otherwise, there were warm melting conditions.

3.5. Co-Location and Comparison

For comparison with the MODIS OWSI chart, the SIC products were rectified to the Mercator projection of the MODIS chart with nearest neighbor sampling. Inside each SIC product pixel, the MODIS SIC was calculated as the fraction of sea ice pixels to the total number of pixels. To reduce the land spillover effect, only those radiometer SIC pixels which had, at maximum, 10% land in the MODIS land mask were included in the comparison.
In similar way, the SIC products were rectified to the UTM34 projection of the S-2 OWSI chart, and then the S-2 SIC was calculated within the SIC product grid cells which had, at maximum, 10% land in the S-2 land mask.
In selecting radiometer data pixels for comparison, the status flag of the SIC retrieval in the OSI-450, OSI-458, and SICCI-HR-SIC products was utilized. This flag showed when an open water filter and a land spill over filter had been applied. We only selected pixels with no filters applied in order to only have ‘true’ retrieved open water pixels in the comparison. There was also a status flag in the OSI-408-a product, but it did not have the needed information. The ASI-AMSR2 product did not have a status flag at all. For these two products, we excluded pixels with a radiometer SIC of exactly zero. These pixels were assumed to be mostly the result of applied filters. Keeping the rejected pixels would have had a large effect on the comparison statistics.
The comparison between the gridded MODIS or S-2 SIC and co-located radiometer SIC was conducted by calculating the mean difference of the radiometer SIC minus the MODIS/S-2 SIC, standard deviation of the difference, and linear regression line and coefficient of determination of the regression ( R 2 ).

4. Results

4.1. Comparison with the MODIS OWSI Chart

Statistics on the comparison between the MODIS SIC and different radiometer SIC products are shown in Table 2. Examples of the MODIS SIC and OSI-408-a SIC images, both with a 10 km pixel size, are shown in Figure 7. The images also show the coastline to illustrate how far the land mask extends from the coast. Here, the land mask extension is typically 2–3 pixels, at which distance there can be some leftover land spill effect after correcting or filtering the measured T B s. A scatterplot between the MODIS and OSI-408-a SIC datasets is shown in Figure 8.
For the OSI-408-a SIC, the overall mean difference is small, only +1.1%; however, at high SICs, the OSI-408-a SIC underestimates the MODIS SIC. For example, in the 85 to 95% MODIS SIC range, the average OSI-408-a SIC is only 76.0%. At small MODIS SICs, there is on the contrary overestimation; for example, in the 5 to 15% MODIS SIC range, the average OSI-408-a SIC is 37.0%. The standard deviation of the difference is large, at 19.7%. Dividing the data into winter and melting cases does not change statistics much. For example, the standard deviation is 20.5% for the winter case and 18.1% for the melting one. The slope of the linear regression is 0.71 and the R 2 is 0.80 shows a rather poor one-to-one linear agreement between the datasets.
The mean difference for the OSI-450-a SIC is 10.6%. At a high SIC, there is an underestimation of the MODIS SIC; in the 85 to 95% MODIS SIC range, the average OSI-450-a SIC is only 78.8%. Also, at a small SIC, we find overestimation; in the 5 to 15% MODIS SIC range, the average is 41.9%. The standard deviation of the difference is large here, at 21.0%, and it is almost the same for the winter and melting conditions. The slope is only 0.58 and the R 2 is 0.79.
For the OSI-458 SIC product, the mean difference is very small, 0.8%. There is a good match with the average when SIC ≥ 60%, but there is overestimation at SIC ≤ 50%. The standard deviation of the difference is modest, at 15.5%, ranging between 13.4% for the melting conditions and 16.2% for the winter. The slope is small, only 0.54, and the R 2 is 0.64.
For the SICCI-HR-SIC product, the scatterplot against the MODIS SIC also shows SIC underestimation at a high SIC; in the 85 to 95% MODIS SIC range, the average SICCI-HR-SIC is 77.5%. Also, there is overestimation at a small SIC; in the 5–15% MODIS SIC range, the mean SIC is 26.0%. The mean difference is 5.3% and the standard deviation of the difference is 17.8% overall, taking a larger value for the winter case, at 19.4%, and a smaller value for the melting case, at 16.4%. The slope is the largest of the products investigated, at 0.77, and the R 2 is 0.80.
The ASI-AMSR2 SIC has a very poor match with the MODIS SIC. The overall mean difference is −12.4%, the standard deviation of the difference is 26.8%, and there is typically underestimation of the MODIS SIC when the radiometer SIC is larger than 40; for example, in the 85 to 95% MODIS SIC range, the average ASI-AMSR2 SIC is only 54.6%. For the winter and melting conditions, the standard deviation of the difference is also large, at 28.2% and 24.0%, respectively. A poor match between the datasets is reflected in the particularly low value for the R 2 of only 0.50, while the slope of the linear regression is only 0.66.

4.2. Comparison with the Sentinel-2 OWSI Chart

Comparison statistics between the S-2 SIC and different SIC products are shown in Table 3. Figure 9 depicts a scatterplot between the S-2 SIC and OSI-408-a data. In the comparison, we again rejected samples where the radiometer SIC was exactly zero or status flags showed the application of open water and land spillover filters.
The mean difference, the OSI-408-a SIC minus the S-2 SIC, is small, only −0.6%, but there is a slight underestimation of the S-2 SIC when the S-2 SIC is larger than 60; for example, in the 85 to 95% S-2 SIC range, the average OSI-408-a SIC is 82.8%. At SICs below 50%, the OSI-408-a SIC overestimates compared to the S-2; in the 5 to 15% S-2 SIC range, the average OSI-408-a SIC is 31.0%. The standard deviation of the difference is modest, at 15.5%, ranging between 13.5% for winter conditions and 16.1% for melting conditions. The linear one-to-one agreement between the SIC datasets is not very high; the slope is 0.70 and the R 2 is 0.80, almost identical to the values obtained for the comparison to the MODIS SIC for this product.
For the OSI-450-a and OSI-458 products, the accuracy of the results of the comparison is limited by the small number of samples (max. 242). The mean difference is 1.3% for the OSI-450-a and 0.4% for the OSI-458. The standard deviation of the difference is large, at 22.0% and 20.2%, respectively. The slope is only 0.50 for the OSI-450-a and 0.40 for the OSI-458, which is considerably worse than that for the MODIS SIC.
For the SICCI-HR-SIC, the mean difference is small, at 2.5%, but the S-2 SIC is underestimated when SIC > 70% and overestimated when SIC ≤ 40%. The standard deviation of the difference is again large, at 22.1% for all data, with a substantial variation between the winter conditions (18.6%) and melting conditions (23.6%). The one-to-one linear agreement is quite poor, with a slope of 0.67 and R 2 of 0.68.
The ASI-AMSR2 underestimates the S-2 SIC when SIC ≥ 50%. The mean difference is large, at −12.0%. The standard deviation of the difference is again large, at 24.0% for all data, 26.6% for the winter, and 20.8% for the melting conditions. The slope of the linear regression is the highest of the products compared to the S-2 SIC, at 0.76, but the R 2 is only 0.48.

5. Discussion and Conclusions

Compared to the Arctic Sea ice statistics in [5], the standard deviations of the differences are here larger. For example, in [5], the largest observed standard deviation was 11.7% in the Arctic FYI case. Here, the smallest standard deviation was 15.5% and the largest was 26.8%. A part of the large differences in the standard deviation could be due to the usage of Arctic tie points for cold conditions under Baltic Sea melting conditions. Further, the Arctic tie points in the SIC algorithms may not fit with the low-salinity Baltic Sea ice, but this is likely a minor error source.
A common feature of the OSI SAF and SICCI-HR-SIC datasets was the typical underestimation of the MODIS SIC at large SIC values (SIC > 60%) and its overestimation at small SIC values (SIC < 40%). A similar underestimation of the S-2 SIC was also present when it was larger than 60 or 70%, and its overestimation typically happened when SIC ≤ 40 or 50%. In the case of the ASI-AMSR2 SIC, the underestimation of the MODIS and S-2 SIC occurred over a wider SIC range, SIC > 50%. This underestimation was partly due to a possible positive bias in the MODIS/S-2 SIC, because a pixel with even a tiny fraction of highly reflective sea ice, i.e., snow-covered sea ice, was classified as sea ice with the reflectance thresholding method. This underestimation effect was larger for finer-resolution SIC products, as the same number of MODIS/S-2 pixels classified as sea ice resulted in a larger aggregated MODIS/S-2 SIC. For example, 100 such MODIS pixels resulted in a SIC of 1% in the 25 km grid, but the 10 km grid gave a SIC of 6.25%. Another source for the underestimation is the occurrence of thin ice [2,42,43]. The SIC overestimation at a small SIC could be partly due to remnant land contamination, i.e., the land contamination correction was not perfect on the measured T B data.
The one-to-one linear agreement between the radiometer SIC dataset and the MODIS/S-2 SIC was always quite poor; the slope of the linear regression was from 0.40 to 0.77 and the R 2 ranged from 0.26 to 0.88.
Dividing the data into winter and melting conditions did not have a clear consistent effect on the statistics. In the case of the comparison against the MODIS SIC data, the standard deviation value of the difference typically slightly increased for the winter conditions compared to the value for all data, and it slightly decreased for the melting conditions. However, we did not always observe such a behavior in the comparison against the S-2 SIC data. This could be due to different data amounts in the MODIS and S-2 SIC cases and complex changes in snow and sea ice properties from winter to melting conditions and vice versa.
It is difficult to rank the SIC products based on the comparison results. For example, the OSI-458 had a smaller standard deviation of the difference (15.5%) than the MODIS SIC data, but against the S-2 SIC data, the standard deviation of the difference was considerably larger than the smaller OSI-408-a one (20.2% vs. 15.5%).
The two CDRs, the OSI-450 and OSI-458, were retrieved with the same algorithm. Only the input data were different, with SSMIS for the OSI-450 and AMSR2 for the OSI-458. The standard deviation of the difference was smaller for the OSI-458 product when against both the MODIS and S-2 data, as it should have been based on results in [5], because the OSI-458 SIC was computed from finer-resolution AMSR2 input data.
The worst overall performance was surprisingly observed for the ASI-AMSR2 product, which had the finest resolution of all the SIC products. Finer resolution meant larger spatial variation in the SIC and a larger effect from the spatial mismatch of an ice area between the radiometer SIC and MODIS/S-2 products (the radiometer SIC is a daily product whereas the MODIS/S-2 SIC is an instantaneous swath product). The ASI-AMSR2 product, using 89 GHz channels, was sensitive to atmospheric effects, which were not corrected for. However, cloud liquid water does not have role here, as clear-sky SIC data were used and, thus, the atmospheric effects were likely not that large. The other radiometer SIC products included the atmospheric correction of the T B data, and less sensitive lower-frequency channels were used in the SIC retrieval.
In summary, none of the SIC products performed well over the Baltic Sea ice, and they should be used with care in Baltic Sea ice monitoring and studies. The accuracy of the products could be improved by diminishing any further land spillover effects on the measured T B data and using Baltic Sea-specific tie points.

Author Contributions

Conceptualization, M.M., S.K., and R.T.; methodology, M.M., S.K., and R.T.; software, M.M.; validation, M.M.; formal analysis, M.M.; investigation, M.M.; resources, M.M.; data curation, M.M.; writing—original draft preparation, M.M.; writing—review and editing, M.M., S.K., and R.T.; visualization, M.M.; supervision, M.M.; project administration, M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Andersen, S.; Tonboe, R.; Kern, S.; Schyberg, H. Improved Retrieval of Sea Ice Total Concentration from Spaceborne Passive Microwave Observations Using Numerical Weather Prediction Model Fields: An Intercomparison of Nine Algorithms. Remote Sens. Environ. 2006, 104, 374–392. [Google Scholar] [CrossRef]
  2. Ivanova, N.; Pedersen, L.T.; Tonboe, R.T.; Kern, S.; Heygster, G.; Lavergne, T.; Sørensen, A.; Saldo, R.; Dybkjær, G.; Brucker, L.; et al. Inter-Comparison and Evaluation of Sea Ice Algorithms: Towards Further Identification of Challenges and Optimal Approach Using Passive Microwave Observations. Cryosphere 2015, 9, 1797–1817. [Google Scholar] [CrossRef]
  3. Kern, S.; Lavergne, T.; Notz, D.; Pedersen, L.T.; Tonboe, R.T.; Saldo, R.; Sørensen, A.M. Satellite Passive Microwave Sea-Ice Concentration Data Set Intercomparison: Closed Ice and Ship-Based Observations. Cryosphere 2019, 13, 3261–3307. [Google Scholar] [CrossRef]
  4. Kern, S.; Lavergne, T.; Notz, D.; Pedersen, L.T.; Tonboe, R. Satellite Passive Microwave Sea-Ice Concentration Data Set Inter-Comparison for Arctic Summer Conditions. Cryosphere 2020, 14, 2469–2493. [Google Scholar] [CrossRef]
  5. Kern, S.; Lavergne, T.; Pedersen, L.T.; Tonboe, R.T.; Bell, L.; Meyer, M.; Zeigermann, L. Satellite Passive Microwave Sea-Ice Concentration Data Set Intercomparison Using Landsat Data. Cryosphere 2022, 16, 349–378. [Google Scholar] [CrossRef]
  6. Maaß, N.; Kaleschke, L. Improving Passive Microwave Sea Ice Concentration Algorithms for Coastal Areas: Applications to the Baltic Sea. Tellus A 2010, 62, 393–410. [Google Scholar] [CrossRef]
  7. Hallikainen, M. Microwave Remote Sensing of Low-Salinity Sea Ice. In Geophysical Monograph Series; Carsey, F.D., Ed.; American Geophysical Union: Washington, DC, USA, 1992; Volume 68, pp. 361–373. ISBN 978-0-87590-033-9. [Google Scholar]
  8. Lavergne, T.; Sørensen, A.M.; Kern, S.; Tonboe, R.; Notz, D.; Aaboe, S.; Bell, L.; Dybkjær, G.; Eastwood, S.; Gabarro, C.; et al. Version 2 of the EUMETSAT OSI SAF and ESA CCI Sea-Ice Concentration Climate Data Records. Cryosphere 2019, 13, 49–78. [Google Scholar] [CrossRef]
  9. Toudal Pedersen, L.; Dybkjær, G.; Eastwood, S.; Heygster, G.; Ivanova, N.; Kern, S.; Lavergne, T.; Saldo, R.; Sandven, S.; Sørensen, A.; et al. ESA Sea Ice Climate Change Initiative (Sea_Ice_cci): Sea Ice Concentration Climate Data Record from the AMSR-E and AMSR-2 Instruments at 25 km Grid Spacing, Version 2.1; 2017. Available online: https://catalogue.ceda.ac.uk/uuid/f17f146a31b14dfd960cde0874236ee5 (accessed on 1 November 2024).
  10. Peng, G.; Meier, W.N.; Scott, D.J.; Savoie, M.H. A Long-Term and Reproducible Passive Microwave Sea Ice Concentration Data Record for Climate Studies and Monitoring. Earth Syst. Sci. Data 2013, 5, 311–318. [Google Scholar] [CrossRef]
  11. Meier, W.N.; Fetterer, F.; Savoie, M.; Mallory, S.; Duerr, R.; Stroeve, J. NOAA/NSIDC Climate Data Record of Passive Microwave Sea Ice Concentration, Version 3; 2017. Available online: https://nsidc.org/data/g02202/versions/3 (accessed on 1 November 2024).
  12. Seinä, A.; Peltola, J. Duration of Ice Season and Statistics of Fast Ice Thickness along the Finnish Coast 1961–1990; Finnish Marine Research Report No. 258; Finnish Institute of Marine Research: Helsinki, Finland, 1991; pp. 1–46. [Google Scholar]
  13. Vihma, T.; Haapala, J. Geophysics of Sea Ice in the Baltic Sea: A Review. Prog. Oceanogr. 2009, 80, 129–148. [Google Scholar] [CrossRef]
  14. Ice Season in the Baltic Sea—Finnish Meteorological Institute. Available online: https://en.ilmatieteenlaitos.fi/ice-season-in-the-baltic-sea (accessed on 8 October 2020).
  15. Seinä, A.; Palosuo, E. The Classification of the Maximum Annual Extent of Ice Cover in the Baltic Sea 1720–1995; Meri Report No. 27; Finnish Institute of Marine Research: Helsinki, Finland, 1996; pp. 79–91. [Google Scholar]
  16. Schmelzer, N.; Holfort, J. Climatological Ice Atlas for the Western and Southern Baltic Sea (1961–2010); Digital Supplement: Comparison of Ice Conditions in the 30-Year Periods 1961–1990, 1971–2000, 1981–2010; Bundesamt für Seeschifffahrt und Hydrographie (BSH): Hamburg/Rostock, Germany, 2012. [Google Scholar]
  17. Leppäranta, M.; Hakala, R. The Structure and Strength of First-Year Ice Ridges in the Baltic Sea. Cold Reg. Sci. Technol. 1992, 20, 295–311. [Google Scholar] [CrossRef]
  18. Spreen, G.; Kaleschke, L.; Heygster, G. Sea Ice Remote Sensing Using AMSR-E 89-GHz Channels. J. Geophys. Res. 2008, 113, 2005JC003384. [Google Scholar] [CrossRef]
  19. Boardo, F.; Tonboe, R.; Howe, E. Algorithm Theoretical Basis Document for Global Sea Ice Concentration Level 2 and Level 3; Danish Meteorological Institute: Copenhagen, Denmark, 2023. [Google Scholar]
  20. Lavergne, T.; Sorensen, A.; Tonboe, R.; Saldo, R.; Pedersen, L.T.; Strong, C.; Cherkaev, E.; Golden, K.M.; Eastwood, S. Global Sea Ice Concentration Climate Data Records, Algorithm Theoretical Basis Document; EUMETSAT OSI SAF: Oslo, Norway, 2022. [Google Scholar]
  21. Lavergne, T.; Sorensen, A.; Tonboe, R.; Pedersen, L.T. CCI+ Sea Ice ECV Sea Ice Concentration Algorithm Theoretical Basis Document; Norwegian Meteorological Institute: Oslo, Norway, 2021; p. 41. [Google Scholar]
  22. OSI SAF OSI-408-a Global AMSR Sea Ice Concentration—GCOM-W1 2017. Available online: https://navigator.eumetsat.int/product/EO:EUM:DAT:GCOM:OSI-408 (accessed on 1 November 2024).
  23. Lavergne, T.; Sørensen, A.; Tonboe, R.T.; Kreiner, M.; Saldo, R.; Birkedal, A.; Baordo, F.; Aspenes, T.; Eastwood, S. ESA Sea Ice Climate Change Initiative (Sea_Ice_cci): High(Er) Resolution Sea Ice Concentration Climate Data Record Version 3 (SSM/I and SSMIS) 2023, 21809 Files, 83566854591 B. Available online: https://catalogue.ceda.ac.uk/uuid/eade27004395466aaa006135e1b2ad1a (accessed on 1 November 2024).
  24. Kaleschke, L. AWI AMSR2 Ice Concentration Product V110; Alfred Wegener Institute: Bremen, Germany, 2022. [Google Scholar]
  25. Mäkynen, M.; Karvonen, J. MODIS Sea Ice Thickness and Open Water–Sea Ice Charts over the Barents and Kara Seas for Development and Validation of Sea Ice Products from Microwave Sensor Data. Remote Sens. 2017, 9, 1324. [Google Scholar] [CrossRef]
  26. Barnes, W.L.; Pagano, T.S.; Salomonson, V.V. Prelaunch Characteristics of the Moderate Resolution Imaging Spectroradiometer (MODIS) on EOS-AM1. IEEE Trans. Geosci. Remote Sens. 1998, 36, 1088–1100. [Google Scholar] [CrossRef]
  27. Ackerman, S.; Frey, R.; Strabala, K.; Liu, Y.; Gumley, L.; Baum, B.; Menzel, P. Discriminating Clear-Sky from Cloud with MODIS Algorithm Theoretical Basis Document (MOD35), Version 6.1; Cooperative Institute for Meteorological Satellite Studies, University of Wisconsin: Madison, WI, USA, 2010. [Google Scholar]
  28. Riggs, G.; Hall, D. MODIS Sea Ice Products User Guide to Collection 6; NSIDC: Boulder, CO, USA, 2015. [Google Scholar]
  29. Nishihama, M.; Wolfe, R.; Solomon, D.; Patt, F.; Blanchette, J.; Fleig, A.; Masuoka, E. MODIS Level 1a Earth Location: Algorithm Theoretical Basis Document Version 3.0; MODIS Science Data Support Team. NASA/Goddard Spaceflight Center: Greenbelt, MD, USA, 1997. [Google Scholar]
  30. Parera-Portell, J.A.; Ubach, R.; Gignac, C. An Improved Sea Ice Detection Algorithm Using MODIS: Application as a New European Sea Ice Extent Indicator. Cryosphere 2021, 15, 2803–2818. [Google Scholar] [CrossRef]
  31. Siitam, L.; Sipelgas, L.; Pärn, O.; Uiboupin, R. Statistical Characterization of the Sea Ice Extent during Different Winter Scenarios in the Gulf of Riga (Baltic Sea) Using Optical Remote-Sensing Imagery. Int. J. Remote Sens. 2017, 38, 617–638. [Google Scholar] [CrossRef]
  32. Riggs, G.A.; Hall, D.K.; Ackerman, S.A. Sea Ice Extent and Classification Mapping with the Moderate Resolution Imaging Spectroradiometer Airborne Simulator. Remote Sens. Environ. 1999, 68, 152–163. [Google Scholar] [CrossRef]
  33. Jiang, L.; Ma, Y.; Chen, F.; Liu, J.; Yao, W.; Shang, E. Automatic High-Accuracy Sea Ice Mapping in the Arctic Using MODIS Data. Remote Sens. 2021, 13, 550. [Google Scholar] [CrossRef]
  34. Drusch, M.; Del Bello, U.; Carlier, S.; Colin, O.; Fernandez, V.; Gascon, F.; Hoersch, B.; Isola, C.; Laberinti, P.; Martimort, P.; et al. Sentinel-2: ESA’s Optical High-Resolution Mission for GMES Operational Services. Remote Sens. Environ. 2012, 120, 25–36. [Google Scholar] [CrossRef]
  35. ESA. Sentinel-2 User Handbook; ESA Standard Document; ESA: Paris, France, 2015. [Google Scholar]
  36. Sentinel-2—Missions. Available online: https://sentinel.esa.int/web/sentinel/missions/sentinel-2 (accessed on 19 October 2023).
  37. Sentinel-2 Level-1C Cloud Masks. Available online: https://sentinel.esa.int/web/sentinel/technical-guides/sentinel-2-msi/level-1c/cloud-masks (accessed on 19 October 2023).
  38. Ludwig, V.; Spreen, G.; Pedersen, L.T. Evaluation of a New Merged Sea-Ice Concentration Dataset at 1 Km Resolution from Thermal Infrared and Passive Microwave Satellite Data in the Arctic. Remote Sens. 2020, 12, 3183. [Google Scholar] [CrossRef]
  39. Muchow, M.; Schmitt, A.U.; Kaleschke, L. A Lead-Width Distribution for Antarctic Sea Ice: A Case Study for the Weddell Sea with High-Resolution Sentinel-2 Images. Cryosphere 2021, 15, 4527–4537. [Google Scholar] [CrossRef]
  40. Wang, Z.; Sun, P.; Wang, L.; Zhang, M.; Wang, Z. Monitoring Sea Ice in Liaodong Bay of Bohai Sea during the Freezing Period of 2017/2018 Using Sentinel-2 Remote Sensing Data. J. Spectrosc. 2021, 2021, 9974845. [Google Scholar] [CrossRef]
  41. McFeeters, S.K. The Use of the Normalized Difference Water Index (NDWI) in the Delineation of Open Water Features. Int. J. Remote Sens. 1996, 17, 1425–1432. [Google Scholar] [CrossRef]
  42. Kwok, R.; Comiso, J.C.; Martin, S.; Drucker, R. Ross Sea Polynyas: Response of Ice Concentration Retrievals to Large Areas of Thin Ice. J. Geophys. Res. Ocean. 2007, 112, C03S21. [Google Scholar] [CrossRef]
  43. Shokr, M.; Kaleschke, L. Impact of Surface Conditions on Thin Sea Ice Concentration Estimate from Passive Microwave Observations. Remote Sens. Environ. 2012, 121, 36–50. [Google Scholar] [CrossRef]
Figure 1. An example of the MODIS RGB images with a band 2-1-1 combination over the Bay of Bothnia (a), the same RGB image with the MOD35_L2 cloud mask (confident cloudy in the unobstructed field-of-view quality flag) (b), and the RGB image with the manually edited cloud mask (c). In figure (b), yellow ovals highlight areas where clear sky lead pixels are classified as cloudy. The Aqua MODIS data were acquired on 11 February 2019, 11:00 UTC. The pixel size of the RGB images is 250 m. Clouds are presented in purple color.
Figure 1. An example of the MODIS RGB images with a band 2-1-1 combination over the Bay of Bothnia (a), the same RGB image with the MOD35_L2 cloud mask (confident cloudy in the unobstructed field-of-view quality flag) (b), and the RGB image with the manually edited cloud mask (c). In figure (b), yellow ovals highlight areas where clear sky lead pixels are classified as cloudy. The Aqua MODIS data were acquired on 11 February 2019, 11:00 UTC. The pixel size of the RGB images is 250 m. Clouds are presented in purple color.
Remotesensing 16 04430 g001
Figure 2. Probability density functions of MODIS band 1 reflectances (a) for open water, thin ice, and snow-covered ice in winter conditions and (b) for open water, melting ice, and snow-covered ice in spring conditions.
Figure 2. Probability density functions of MODIS band 1 reflectances (a) for open water, thin ice, and snow-covered ice in winter conditions and (b) for open water, melting ice, and snow-covered ice in spring conditions.
Remotesensing 16 04430 g002
Figure 3. MODIS open water–sea ice chart over the Bay of Bothnia on 11 February 2019, 11:00 UTC, based on the band 1 reflectance. Land is black, clouds are purple, open water is blue, and sea ice is white. The pixel size is 250 m.
Figure 3. MODIS open water–sea ice chart over the Bay of Bothnia on 11 February 2019, 11:00 UTC, based on the band 1 reflectance. Land is black, clouds are purple, open water is blue, and sea ice is white. The pixel size is 250 m.
Remotesensing 16 04430 g003
Figure 4. A sentinel-2 RGB image from bands 4-3-2 (approximating a true color image) with the opaque cloud mask (a) from the tile T34WES acquired on 5 March 2019, 10:10 UTC. The (b) image shows the cloud mask after manual editing. Clouds are in a purple color. The pixel size is 10 m.
Figure 4. A sentinel-2 RGB image from bands 4-3-2 (approximating a true color image) with the opaque cloud mask (a) from the tile T34WES acquired on 5 March 2019, 10:10 UTC. The (b) image shows the cloud mask after manual editing. Clouds are in a purple color. The pixel size is 10 m.
Remotesensing 16 04430 g004
Figure 5. Probability density functions of Sentinel-2 band 2, 3, 4, and 8 reflectances from imagery tile 34WFT acquired on (a) 5 March and (b) 6 April 2009.
Figure 5. Probability density functions of Sentinel-2 band 2, 3, 4, and 8 reflectances from imagery tile 34WFT acquired on (a) 5 March and (b) 6 April 2009.
Remotesensing 16 04430 g005
Figure 6. Sentinel-2 open water–sea ice chart for the tile T34WES acquired on 5 Mar 2019, 10:10 UTC, based on the band 4 reflectance. Land is black, clouds are purple, open water is blue, and sea ice is white. The pixel size is 10 m.
Figure 6. Sentinel-2 open water–sea ice chart for the tile T34WES acquired on 5 Mar 2019, 10:10 UTC, based on the band 4 reflectance. Land is black, clouds are purple, open water is blue, and sea ice is white. The pixel size is 10 m.
Remotesensing 16 04430 g006
Figure 7. (a) MODIS SIC chart of a 10 km pixel size on 12 February 2019, 09:50 UTC. Black dots show the locations of the coastal weather stations, from north to south: Kemi 1, Ulkokalla, and Valassaaret. Dark blue (SIC value: 10%) shows missing data due to cloud cover. (b) An OSI-408-a daily SIC chart on 12 February 2019. Both images have a 10 km pixel size.
Figure 7. (a) MODIS SIC chart of a 10 km pixel size on 12 February 2019, 09:50 UTC. Black dots show the locations of the coastal weather stations, from north to south: Kemi 1, Ulkokalla, and Valassaaret. Dark blue (SIC value: 10%) shows missing data due to cloud cover. (b) An OSI-408-a daily SIC chart on 12 February 2019. Both images have a 10 km pixel size.
Remotesensing 16 04430 g007
Figure 8. Scatterplot between MODIS SIC and OSI-408-a SIC data for February–April 2019. Blue line is linear regression line. Red triangles are mean OSI-408-a SIC calculated for MODIS SIC ranges of 0–5%, 5–15%, …, 85–95%, and 95–100%, and red vertical bars are ±1 standard deviations. Pixel size of both datasets is 10 km.
Figure 8. Scatterplot between MODIS SIC and OSI-408-a SIC data for February–April 2019. Blue line is linear regression line. Red triangles are mean OSI-408-a SIC calculated for MODIS SIC ranges of 0–5%, 5–15%, …, 85–95%, and 95–100%, and red vertical bars are ±1 standard deviations. Pixel size of both datasets is 10 km.
Remotesensing 16 04430 g008
Figure 9. Scatterplot between Sentinel-2 SIC and OSI-408-a SIC data for February–April 2019. Blue line is linear regression line. Red triangles are mean OSI-408-a SIC calculated for Sentinel-2 SIC ranges of 0–5%, 5–15%, …, 85–95%, and 95–100%, and red vertical bars are ±1 standard deviations. Pixel size of both datasets is 10 km.
Figure 9. Scatterplot between Sentinel-2 SIC and OSI-408-a SIC data for February–April 2019. Blue line is linear regression line. Red triangles are mean OSI-408-a SIC calculated for Sentinel-2 SIC ranges of 0–5%, 5–15%, …, 85–95%, and 95–100%, and red vertical bars are ±1 standard deviations. Pixel size of both datasets is 10 km.
Remotesensing 16 04430 g009
Table 1. Microwave radiometer SIC products evaluated over the Baltic Sea.
Table 1. Microwave radiometer SIC products evaluated over the Baltic Sea.
ProductInput dataGridReferences
OSI-408-aAMSR2: 18.7 GHz V-pol, 36.5 GHzPolarStereo 10 km[19,22]
OSI-450-aSSMIS: 19.35 GHz V-pol, 37 GHzEASE2.0 25 km[8,20]
OSI-458AMSR2: 18.7 GHz V-pol, 36.5 GHzEASE2.0 25 km[8,20]
SICCI-HR-SICSSMIS: 19.35 GHz V-pol, 37 GHz, 91.7 GHzEASE2.0 12.5 km[21,23]
ASI-AMSR2AMSR2: 89 GHzPolarStereo 3.125 km[18,24]
Table 2. Statistical parameters between the radiometer SIC and MODIS SIC using all data. The mean difference is the radiometer SIC minus the MODIS SIC. The slope and intercept are the coefficients of linear regression, and R2 is the coefficient of determination. N is the number of data pairs.
Table 2. Statistical parameters between the radiometer SIC and MODIS SIC using all data. The mean difference is the radiometer SIC minus the MODIS SIC. The slope and intercept are the coefficients of linear regression, and R2 is the coefficient of determination. N is the number of data pairs.
OSI-408-aOSI-450-aOSI-458SICCI-HR-SICASI-AMSR2
Mean difference [%]1.110.60.85.3−12.4
Std of difference [%]19.721.015.517.826.8
Slope0.710.580.540.770.66
Intercept17.029.537.011.614.2
R20.800.790.640.800.50
N8101224991218,291114,733
Table 3. Statistical parameters between the radiometer SIC and Sentinel-2 SIC using all data. The mean difference is the radiometer SIC minus the Sentinel-2 SIC. The slope and intercept are the coefficients of linear regression, and R2 is the coefficient of determination. N is the number of data pairs.
Table 3. Statistical parameters between the radiometer SIC and Sentinel-2 SIC using all data. The mean difference is the radiometer SIC minus the Sentinel-2 SIC. The slope and intercept are the coefficients of linear regression, and R2 is the coefficient of determination. N is the number of data pairs.
OSI-408-aOSI-450-aOSI-458SICCI-HR-SICASI-AMSR2
Mean difference [%]−0.61.30.42.5−12.0
Std of difference [%]15.522.020.222.124.0
Slope0.700.500.400.670.76
Intercept19.932.746.520.67.8
R280.00.590.260.680.48
N1293242139138618,755
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mäkynen, M.; Kern, S.; Tonboe, R. Evaluation of Microwave Radiometer Sea Ice Concentration Products over the Baltic Sea. Remote Sens. 2024, 16, 4430. https://doi.org/10.3390/rs16234430

AMA Style

Mäkynen M, Kern S, Tonboe R. Evaluation of Microwave Radiometer Sea Ice Concentration Products over the Baltic Sea. Remote Sensing. 2024; 16(23):4430. https://doi.org/10.3390/rs16234430

Chicago/Turabian Style

Mäkynen, Marko, Stefan Kern, and Rasmus Tonboe. 2024. "Evaluation of Microwave Radiometer Sea Ice Concentration Products over the Baltic Sea" Remote Sensing 16, no. 23: 4430. https://doi.org/10.3390/rs16234430

APA Style

Mäkynen, M., Kern, S., & Tonboe, R. (2024). Evaluation of Microwave Radiometer Sea Ice Concentration Products over the Baltic Sea. Remote Sensing, 16(23), 4430. https://doi.org/10.3390/rs16234430

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop