Next Article in Journal
Cross Deployment Networking and Systematic Performance Analysis of Underwater Wireless Sensor Networks
Next Article in Special Issue
Establishment of a Site-Specific Tropospheric Model Based on Ground Meteorological Parameters over the China Region
Previous Article in Journal
A Study on the Data Compression Technology-Based Intelligent Data Acquisition (IDAQ) System for Structural Health Monitoring of Civil Structures
Previous Article in Special Issue
GNSS-R Altimetry Performance Analysis for the GEROS Experiment on Board the International Space Station
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sea Ice Detection Based on Differential Delay-Doppler Maps from UK TechDemoSat-1

1
School of Geodesy and Geomatics and Collaborative Innovation Center for Geospatial Technology, Wuhan University, Wuhan 430079, China
2
Key Laboratory of Precise Engineering and Industry Surveying, National Administration of Surveying, Mapping and Geoinformation, Wuhan 430079, China
3
German Research Centre for Geosciences, GFZ Potsdam 14473, Germany
4
Institute of Geodesy and Geoinformation Science, Technische Universität Berlin, Berlin 10623, Germany
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(7), 1614; https://doi.org/10.3390/s17071614
Submission received: 21 June 2017 / Revised: 8 July 2017 / Accepted: 10 July 2017 / Published: 12 July 2017

Abstract

:
Global Navigation Satellite System (GNSS) signals can be exploited to remotely sense atmosphere and land and ocean surface to retrieve a range of geophysical parameters. This paper proposes two new methods, termed as power-summation of differential Delay-Doppler Maps (PS-D) and pixel-number of differential Delay-Doppler Maps (PN-D), to distinguish between sea ice and sea water using differential Delay-Doppler Maps (dDDMs). PS-D and PN-D make use of power-summation and pixel-number of dDDMs, respectively, to measure the degree of difference between two DDMs so as to determine the transition state (water-water, water-ice, ice-ice and ice-water) and hence ice and water are detected. Moreover, an adaptive incoherent averaging of DDMs is employed to improve the computational efficiency. A large number of DDMs recorded by UK TechDemoSat-1 (TDS-1) over the Arctic region are used to test the proposed sea ice detection methods. Through evaluating against ground-truth measurements from the Ocean Sea Ice SAF, the proposed PS-D and PN-D methods achieve a probability of detection of 99.72% and 99.69% respectively, while the probability of false detection is 0.28% and 0.31% respectively.

1. Introduction

Knowledge of sea ice coverage and thickness of the polar region is very important for ship routing [1], oil and gas exploration [2], and global weather studies [3]. Although in situ monitoring may provide continuous sea ice information, the spatial sampling is typically insufficient to provide information on global sea ice extent. Estimating global sea ice coverage with adequate temporal resolution can only be achieved using remote sensing technology.
Global Navigation Satellite System (GNSS) Reflectometry (GNSS-R) has been investigated for many years as an innovative remote sensing technique using GNSS signals reflected off the Earth’s surface [4]. GNSS-R was originally suggested and applied to ocean altimetry in 1993 as a multistatic radar technique [5]. With the availability of GNSS signals from GPS, GLONASS, Galileo and BeiDou satellite constellations, this technique shows great promising. In the following years, applications of GNSS-R were extended to the remote sensing of sea surface roughness [6,7], wind speed and direction retrieval [8,9,10,11], target detection and positioning [12], soil moisture determination [13,14,15,16,17], forest biomass [18,19], oil slick detection [20], tsunami parameter estimation and tsunami detection [21,22], sea surface altimetry [23,24], ocean topography [25,26], coastal and continental water altimetry [27] and snow depth estimation [28]. The potential of GNSS-R Delay-Doppler Maps (DDMs) for sea ice distribution detection is studied in this paper.
GNSS-R-based sea ice detection was first investigated by Komjathy et al. [29] by comparing the waveform peak power of reflected signals over ice obtained from an airborne GNSS receiver with RADARSAT back-scattering echoes. The feasibility of ice sheet sensing using GNSS-R was exploited in [30], which proposed a theoretical model of ice scattering. Then the shape of reflected GPS waveforms was studied in [31] to retrieve the permittivity and roughness of different types of sea ice. The United Kingdom-Disaster Monitoring Constellation (UK-DMC) launched in 2003 provided spaceborne GNSS-R data for sea ice detection for the first time, which captured the coherent reflection component of signals reflected from sea ice [32]. Unfortunately, the relationship between sea ice parameters and waveform peak power has not been established. In addition, airborne and ground-based studies using reflected GNSS signals were also performed to explore sea ice concentration [33]. After the launch of the UK TechDemoSat-1 (TDS-1), a huge number of DDMs are available from the Measurement of Earth Reflected Radio-navigation Signals by Satellite (MERRByS) website at www.merrbys.co.uk [34]. TDS-1 DDM data have been used by many researchers for different applications, including wind speed retrieval [35], ocean wave characteristic inversion [36], Tsunami detection and parameter estimation [37], soil moisture estimation [38], sea surface height estimation [23,39] and sea ice detection [40,41].
The first result of sea ice detection using TDS-1 GNSS-R DDMs was presented in [40], which proposed a pixel number based sea ice detection method. Another sea ice detection algorithm based on the similarity of received reflected waveform or DDM, which performs well on sea ice detection, was proposed by Alonso-Arroyo, et al. [42]. The goal of the paper presented here is to explore a new approach through estimating the degree of difference between adjacent DDMs with the capability of distinguishing between sea ice and water. A new GNSS-R differential DDM (dDDM) observable derived from adjacent DDMs is proposed firstly. Then, two different methods are proposed to characterize the dDDMs, which are power-summation (PS) and pixel-number (PN) based methods. The different transition types (water-water, water-ice, ice-ice and ice-water) can be classified according to the characteristics of dDDMs. Several data processing strategies such as dataset selection, noise floor subtraction, and normalization of DDMs are employed to reduce the impact of noise floor and speckle noise. In addition, an adaptive incoherent averaging technique is employed to improve the computational efficiency. Finally, the performance of the proposed methods is evaluated against the existing methods using a large number of TDS-1 DDMs and ground-truth sea ice edge maps.
The next section of this paper presents the details of the two proposed ice detection methods, including the definition of the new observation variable, the derivation of the formulas, the data processing techniques, and the detection procedures. Section 3 provides the experimental results obtained through processing the TDS-1 DDM dataset using the proposed methods and the existing methods. The final section summarizes the scientific results and indicates future directions for performance enhancement.

2. Theory and Method

2.1. GNSS-R Delay-Doppler Maps

The TDS-1, launched in July 2014, was placed into a quasi-Sun synchronous orbit with an altitude of ~635 km and an inclination of 98.4°. The Space GNSS Receiver Remote Sensing Instrument (SGR-ReSI) built by Surrey Satellite Technology Ltd. (SSTL, Guildford, UK) is one of the key payloads carried on TDS-1. The ReSI supports the GPS L1 frequency band using a downward pointing antenna with a gain of ~13 dBi. More details about TDS-1 and the SGR-ReSI are introduced in [43].
The TDS-1 receiver produces DDMs through cross correlating received scattered signals with locally generated and time-synchronized code replicas for different path delays and Doppler shifts [38]. Speckle noise has a significant influence on cross correlation, which is usually performed over 1 ms and less 20 ms of coherent integration time, so an incoherent summation of consecutive cross correlation values is employed to mitigate the noise. The typical incoherent summation time for TDS-1 DDMs is 1 s [43]. The spreading of DDM reveals the degree of roughness of the observing surface. It is well known that the surface of sea water is usually rougher than that of young ice, so the spreading in Delay and Doppler of DDMs collected over sea ice surface would be less than that of sea water, whose DDMs usually exhibit a horseshoe shape [7]. Typical DDMs from TDS-1 collected over sea ice and water are shown in Figure 1. Figure 1a depicts a DDM collected over an ice field, while a DDM of sea water is shown in Figure 1b, which clearly indicates that the DDM is heavily affected by noise and sea state.

2.2. GNSS-R Observables

GNSS-R DDMs are usually processed to produce observables that could describe the scattering strength of reflected surface. Different surface roughness related observables have been proposed in [44,45]. Those DDM observables have been employed in [40] with slight modification, which applied a pixel number based sea ice detection method with an incoherent summation.
The standard DDM observable is usually denoted as D D M ( τ , f ) , where τ and f represent delay and Doppler shift, respectively. Then the corresponding normalized DDM denoted as D D M ¯ ( τ , f ) is formulated as:
D D M ¯ ( τ , f ) = D D M ( τ , f ) / | D D M | max ,
where | D D M | max is the maximum absolute value of the standard DDM.

2.2.1. Forward GNSS-R Observables

One existing observable derived from D D M ¯ ( τ , f ) was defined in [45] and later was modified as weighted area in [44], which can be formulated as:
W e i g h t e d A r e a = D D M ¯ ( τ , f ) > thre s h o l d D D M ¯ ( τ , f ) δ τ δ f ,
where δ τ and δ f are delay and Doppler shift resolutions respectively. If δ τ δ f is set as one, this observable represents the normalized power summation of the grids whose D D M ¯ ( τ , f ) values are greater than a preset threshold, which is defined through specific processes described in Section 2.4.2.
A similar observable modified from an observable defined in [45] is defined in [40] as the number of grids (i.e., pixel number) whose values greater than a given threshold:
A r e a = DDM ¯ ( τ , f ) > t h r e s h o l d δ τ δ f .
If δ τ δ f is set as one, this observable represents the pixel number of D D M ¯ ( τ , f ) whose component values are greater than a given threshold, which is defined through specific processes described in Section 2.4.2.

2.2.2. Proposed GNSS-R Observables

In this paper, an observable called dDDM observable is proposed for distinguishing sea ice from water. Details of the proposed observable are provided as follows.
Two adjacent normalized DDM observables can be expressed as D D M i ¯ ( τ , f ) and D D M i + 1 ¯ ( τ , f ) respectively. Then the dDDM of these two adjacent normalized DDMs can be simply formulated as,
DDM D ( τ , f ) = D D M i ¯ ( τ , f ) D D M i + 1 ¯ ( τ , f ) ,
where DDM D ( τ , f ) denotes the dDDM. The dDDM can be normalized using Equation (1) and the normalized dDDM is denoted as DDM D ¯ ( τ , f ) .
Considering that the power value of DDM D ( τ , f ) can be either positive or negative, two new observables modified from the observable W e i g h t e d A r e a (Equation (2)) and A r e a (Equation (3)) ( W e i g h t e d A r e a and A r e a are math type variables) respectively are defined as,
P S = | D D M D ¯ ( τ , f ) | > D D M T D D M D ¯ ( τ , f ) d τ d f ,
P N = | D D M D ¯ ( τ , f ) | > D D M T D D M D ¯ ( τ , f ) | D D M D ¯ ( τ , f ) | d τ d f ,
where d τ and d f are delay and Doppler shift resolutions for DDM D ( τ , f ) respectively. If setting d τ d f = 1, PS represents the power summation and PN depicts abstract pixel number of the grids with absolute values ( | DDM D ¯ ( τ , f ) | ) greater than a preset threshold ( D D M T ).

2.3. Data Processing

2.3.1. Dataset Selection and Noise Floor Subtraction

As introduced above, a large number of GNSS-R DDMs have been collected over the Arctic region since the launch of TDS-1. The DDMs generated onboard the TDS-1 satellite are available on the data service platform MERRByS, which provides data track view to help select DDM dataset of the area of interest. Since the purpose is to explore methods for distinguishing between sea ice and water, so only DDMs collected over sea ice and water fields are selected for testing. Those datasets acquired over land are excluded with the geographic information obtained from the U.S. National Oceanic and Atmospheric Administration (NOAA).
Before noise floor subtraction, data alignment is employed to each DDM through adjusting the peak power point to the location without delay and Doppler shift ( τ = 0 , f = 0 ). Note that some DDMs may be seriously affected by noise and cannot demonstrate the characteristics of reflected surface correctly. Thus, DDMs with peak Signal-to-Noise Ratio (SNR) lower than 0 dB are not used. SNR can be calculated as follows:
SNR = ( D D M ( τ , f ) D D M n o i s e ) / D D M n o i s e ,
where D D M ( τ , f ) is the raw power value of DDM, D D M n o i s e is the noise floor, which can be calculated by:
D D M n o i s e = 1 N τ 1 τ m f 1 f m D D M ( τ , f ) ,
where D D M n o i s e represents the noise floor value; τ 1 and τ m are the pixel limits of delay rows of noise box; f 1 and f m are the pixel limits of Doppler shift of noise box; N represents the number of pixels of the noise box. The noise box chosen in this paper spans twenty Doppler bins and five delay chips starting from 1st to 20th in the standard DDM (shown in Figure 2). Then, the DDM values are respectively subtracted by the corresponding noise floor.

2.3.2. Incoherent Averaging and Normalization for Each DDM

Each dataset contains a large number of DDMs and a certain number of consecutive DDMs are associated with the same type of reflected surface. To improve computational efficiency and reduce the impact of speckle noise, incoherent averaging is employed to the DDMs. However, as pointed out in [40], GNSS-R specular point ground track traverses different kinds of surface at some time intervals, and hence one DDM can be related to a transition area with different surface types, such as ice and water. Similar to the idea of Yan [40], an adaptive incoherent averaging is implemented here and a shorter incoherent averaging interval is used for DDMs collected over transition areas. The TDS-1 data contain the specular point positions, so DDMs with specular points near coastline and ice-water transition areas will be processed with a shorter incoherent averaging duration. The sea ice-water transition area can be roughly detected through observing the changing trend of DDM peak power as abrupt power change occurs over sea ice-water transition areas.
After noise reduction and incoherent averaging, normalization is employed to the processed DDMs according to Equation (1). The examples of the resultant DDMs of sea ice and water with 5-s incoherent averaging are shown in Figure 3.

2.3.3. Differential DDM (dDDM) and overall normalization

Let D D M 1 i ¯ ( τ , f ) and D D M 2 i ¯ ( τ , f ) be two adjacent normalized DDMs from sea ice field, let D D M 1 w ¯ ( τ , f ) and D D M 2 w ¯ ( τ , f ) be two adjacent normalized DDMs from sea water field. Then dDDMs between ice and ice, ice and water, water and water can be obtained according to Equation (4). The corresponding equations are shown below:
{ DDM D i i ( τ , f )   =   D D M 1 i ¯ ( τ , f ) D D M 2 i ¯ ( τ , f ) DDM D i w ( τ , f )   =   D D M 2 i ¯ ( τ , f ) D D M 1 w ¯ ( τ , f ) DDM D w w ( τ , f )   =   D D M 1 w ¯ ( τ , f ) D D M 2 w ¯ ( τ , f ) DDM D w i ( τ , f )   =   D D M 2 w ¯ ( τ , f ) D D M 1 i ¯ ( τ , f ) .
Considering that the values of DDM D ( τ , f ) are relatively small, we normalized the various dDDMs by | D D M D ( τ , f ) | max , where | D D M D ( τ , f ) | max is the maximum absolute value of obtained dDDMs. Since it is useful to know what DDM D ( τ , f ) looks like, examples of DDM D ( τ , f ) associated with the four different scenarios are shown in Figure 4. Clearly, these four pictures are significantly different from each other, providing the basis to identify the transition type. To remain consistency and avoid performance degradation, the same normalization with the same maximum value is applied to all the dDDMs of the specific dataset.

2.4. Proposed Identification Approach

As illustrated in Figure 4, the differential power values of DDM D ¯ ( τ , f ) between the same reflected surfaces (Figure 4a,c) basically range from −0.2 to 0.3, while those between different reflected surfaces (Figure 4b,d) have a much wider range and show horseshoe shape. This significant difference provides an opportunity to distinguish dDDMs of different reflected surfaces from those of the same reflected surfaces. Moreover, dDDMs of water-ice transition show a very different power range from those of ice-water transition. This characteristic can be used to distinguish water-ice from ice-water transition. The dDDMs of water-water transition also show different spreading characteristics from those of ice-ice transition. The degree of difference between two adjacent DDMs reflected from ice field is quite small as the surface is suppressed by ice, while dDDMs of water-water transition show larger spreading characteristics than those of ice-ice transition due to the effect of sea surface. The different characteristics can be used to distinguish water-water from ice-ice transition as reflected signals are more sensitive to surface of open water than that of ice covered area. Unfortunately, it seems difficult to build the relationship between the characteristics of dDDMs and the roughness of reflected surface due to the nature of the TDS-1 data (without calibration) and the characteristics of the GPS-R instrument onboard the TDS-1 satellite was not optimized for altimetry purposes. Four scenarios for one dataset which contains a certain number of consecutive DDMs are shown in Figure 5.

2.4.1. Approach Description

The power values of DDM D ( τ , f ) range from −1 to 1, so a power value threshold ( D D M T , which ranges from 0 to 1) can be chosen to identify ice-water and water-ice transition. Another power value threshold D D M T ( D D M T D D M T ) is chosen to distinguish ice-ice transition from water-water transition. Two methods are proposed to identify the transition-case, which are presented as follows:
(a) Power Summation of dDDM (PS-D)
The power summation (PS) of the grids with absolute values ( | DDM D ¯ ( τ , f ) | ) greater than a preset threshold ( D D M T ) can be obtained using Equation (5). Similarly, the power summation (PS’) of the grids with absolute values ( | DDM D ¯ ( τ , f ) | ) greater than a threshold ( D D M T ) can also be obtained using Equation (5). If the power summation threshold for determining ice-water and water-ice transition is defined as P T , the threshold for distinguishing ice-ice transition from water-water transition is defined as P T , the identification procedure is described as follows:
  • If P S > P T , two adjacent DDMs are determined as ice-water transition;
  • If P S < P T , two adjacent DDMs are determined as water-ice transition;
  • If P T P S P T , two adjacent DDMs are determined as the same observing surface and further processing should be done to distinguish between ice-ice and water-water transition:
    If | P S | > P T , two adjacent DDMs are determined as water-water transition;
    If | P S | P T , two adjacent DDMs are determined as ice-ice transition.
(b) Pixel Number of dDDM (PN-D)
Similar to PS-D, the abstract pixel number (PN and PN’) of the grids with absolute values ( | DDM D ¯ ( τ , f ) | ) greater than a preset threshold ( D D M T and D D M T ) of each normalized dDDM can be obtained using Equation (6). If the pixel number threshold for determining ice-water and water-ice transition is defined as N T , the threshold for distinguishing between ice-ice and water-water transition is defined as N T , the corresponding sea ice detection scheme can be presented as follows:
  • If P N > N T , two adjacent DDMs are determined as ice-water transition;
  • If P N < N T , two adjacent DDMs are determined as water-ice transition;
  • If N T P N N T , two adjacent DDMs are determined as the same observing surface and further processing should be done to distinguish between ice-ice and water-water transition:
    If | P N | > N T , two adjacent DDMs are determined as water-water transition;
    If | P N | N T , two adjacent DDMs are determined as ice-ice transition.
As illustrated above, water-ice and ice-water transition can be figured out firstly using PS-D and PN-D while the discrimination between ice-ice and water-water transition need further processing. The observing surface types of DDMs from datasets collected over both ice and water field (Figure 5c) can be derived according to ice-water and water-ice transition. DDMs of datasets only collected over ice (Figure 5a) are distinguished from those of datasets only collected over water (Figure 5b) through using similar scheme for determining ice-water and water-ice transition.

2.4.2. Deriving Thresholds

As discussed above, the scenarios can be divided into different-surface (i.e., ice-water and water-ice transition) and same-surface (i.e., ice-ice and water-water transition). For PS-D based method, if the maximum absolute power summation value of same-surface ( | P S | max s ) is always smaller than the minimum absolute value of different-surface ( | P S | min d ), it will be error free to distinguish between different-surface and same-surface. Similarly, if the maximum absolute pixel number value of same-surface ( | P N | max s ) is smaller than the minimum absolute value of different-surface ( | P N | min d ), different-surface and same-surface can be correctly distinguished from each other.
After determining the thresholds for identifying ice-water and water-ice transition, datasets only collected over only ice or water are used to derive the thresholds for distinguishing between ice-ice and water-water transition. Theoretically, if the maximum absolute PS and PN values of ice-ice transition are always smaller than the minimum absolute values of water-water transition, it will be error free to distinguish ice-ice transition from water-water transition. While the difference between ice-ice transition and water-water transition is much smaller than that between different-surface and same-surface. Therefore, it may be difficult to obtain such error free thresholds. Actually, if more than eighty percent of dDDMs of one dataset collected over only one type of observing surfaces are determined as ice-ice transition, the type of surface of this dataset can be determined as ice.
In order to apply this sea ice detection scheme, six sets of DDM data [46] were used as training data to derive these thresholds so that they may be used as a guideline in the algorithm development. Sea ice edge maps [47] with a resolution of 10 km × 10 km provided by OSI SAF are used as reference data for deriving thresholds. The power value threshold ( D D M T ) for identifying ice-water and water-ice transition is set as 0.2 firstly and increased with a step size of 0.01. Then, different PS and PN results for each D D M T can be obtained according to PS-D and PN-D respectively. After deriving thresholds for identifying ice-water and water-ice transition, datasets only collected over ice or water field can be identified. Then, the power value threshold ( D D M T ) for distinguishing between ice-ice and water-water transition is set as D D M T ( D D M T = 0 . 2 ) firstly and decreased with a step size of 0.01. Then, different PS’ and PN’ results for each D D M T can be obtained according to PS-D and PN-D respectively. For PS-D based method, if the maximum absolute PS value of more than 85% of ice-ice transition dDDMs ( | P S | max i c e ) is smaller than the minimum absolute value of water-water transition ( | P S | min w a t e r ), it will provide the basis for distinguishing ice-ice transition from water-water transition. Similarly, if the maximum absolute PN value of more than 85% of ice-ice transition dDDMs ( | P N | max i c e ) is smaller than the minimum absolute value of water-water transition ( | P N | min w a t e r ), ice-ice transition and water-water transition can be correctly distinguished from each other. It was found that D D M T ranging from 0.22 to 0.53 and D D M T ranging from 0.12 to 0.26 can obtain the best detection results. The process of obtaining thresholds is illustrated in Figure 6.

3. Experimental Results

3.1. Experimental DDM Data Set

The proposed method is evaluated using GNSS-R DDM data produced on TDS-1 when flying over the Arctic region. The Arctic consists of the Arctic Ocean, adjacent seas, and islands. Land within the Arctic region has seasonally varying snow and ice cover, with predominantly treeless permafrost-containing tundra. Arctic seas contain seasonal sea ice in many places. Sea ice of Arctic has a significant effect on global climate and local development, so it is very important to study the ice coverage and its changing trend. The experiments include two parts. Firstly, a case study is implemented through analyzing 131 DDMs from one GNSS-R SP ground track collected on 26 October 2015 in Section 3.2. In Section 3.3, 17346 DDMs from 130 GNSS-R SP ground tracks collected in six separate days from 24 March 2015 to 26 March 2016 have been employed to detect sea ice and water. Meanwhile, the sea ice changing trend over one year can be derived from experimental results.

3.2. Case Study

In order to illustrate the proposed sea ice detection methods, a dataset collected over sea ice and water field is selected to depict the detailed process of distinguishing between sea ice and water. Some basic information about tested dataset is shown in Table 1 and the specular points (SP) ground-track is shown in Figure 7.
The sea ice detection scheme is applied to 131 GNSS-R DDMs and hence 130 dDDMs. For illustration purpose, 30 dDDMs (dDDM 45 to dDDM 74) are displayed in Figure 8.
The three rows of dDDMs presented in Figure 8 depict three different transition cases, respectively. However, in order to analyze a large amount of data automatically by software, PS-D and PN-D methods described in Section 2.4 are employed to obtained dDDMs. Results of PS-D and PN-D under five different power value thresholds are shown in Figure 9a,b, respectively. Both PS and PN of dDDMs show an abrupt change at dDDM 62 from variation to almost a constant value, so this location provides the basis of determining the boundary between water and sea ice. According to the proposed approach, the first 61 dDDMs are identified as water-water transition, while the 62th dDDM is determined as water-ice transition and the next 67 dDDMs are determined as ice-ice transition. Therefore, the first 62 DDMs were collected over sea water field, while the remaining 69 DDMs were collected over sea ice field.
There is a large variability in PS and PN from dDDM 1 to dDDM 37, then a relatively stable variability from dDDM 38 to dDDM 58, and again large change from dDDM 59 to dDDM 61. This phenomenon may be associated with the sea wave height (WH) over the studied area. The significant wave height data (shown in Figure 10) of 26 October 2015 was retrieved from the European Centre for Medium-Range Weather Forecasts (ECMWF, available online: https://www.ecmwf.int) for validating our results. The SP ground-track of 57 DDMs (DDM 1 to DDM 57) is marked in the WH map with a red line, while there is no valid WH data from DDM 58 to DDM 131. The WH of these 57 SP was retrieved through a 2-D linear interpolation according to the SP locations. The values of WH were scaled down as we only concentrate on the variability, which may explain the variability of dDDMs. The relationship between relative wave height and PS of dDDM is presented in Figure 11. The variability of WH can be classified into three stages: (1) the WH changes rapidly from dDDM 1 to dDDM 18 and (2) a relatively slow change from dDDM 19 to dDDM 38 and (3) a decrease from dDDM 38 to dDDM 56. Only the large variability in PS and PN from dDDM 1 to dDDM 18 shows agreement with the changing of WH. Therefore, the variability should be affected by the ocean surface (e.g., wave height, wind speed, sea state etc.), but it is difficult to make any specific conclusions or models about the ocean surface.

3.3. Overall Experimental Results

As introduced in Section 3.1, a number of datasets containing a large number of DDMs spanning over twelve months (March 2015, June 2015, October 2015, January 2016 and March 2016) were employed in this section to validate the proposed sea ice detection approach. Only DDMs collected over sea ice and water field are used and DDMs over land are excluded based on the position of specular point and geographic information. Moreover, datasets collected over sea water far from Arctic center region are also excluded.
In order to assess the performance of the proposed methods, another two observables (power summation of normalized DDM (PS-N) [30], pixel number of normalized DDM (PN-N) [31]) and a matcher filter (MF) approach proposed by Alonso-Arroyo et al. [42] were also used to detect sea ice. Thresholds for best detection results are used for each methods and the probability of false detection of the four methods are presented in Table 2. The sea ice edge data provided by OSI SAF are used as the ground-truth data to examine the detection results. The probability of false detection results for five methods are illustrated in Figure 12. The PS-D and PN-D methods achieve similar performance with a detection probability of 99.72% and 99.69%, respectively, while the detection probabilities of the PS-N, PN-N and MF methods are 96.38%, 96.45% and 98.16%, respectively. Compared to the proposed method PS-D and PN-D, the computation complexity of MF is about 10% lower than that of PS-D and PN-D, but the detection probability of MF is about 1.5% less than that of PS-D and PN-D. The results retrieved by MF show about 97.9% agreement with those obtained by PS-D and PN-D. Moreover, the performance of MF presented in this paper shows good agreement with that presented in [42].
It is useful and interesting to illustrate sea ice detection results directly on a geographical map. Figure 13 shows the sea ice detection results based on the PS-D method using data collected over six periods (24 March 2015, 4 June 2015, 23 August 2015, 26 October 2015, 21 January 2016 and 26 March 2016).
There was a significant decline in the coverage of sea ice over Arctic region from March to August 2015 due to the increase of temperature. On the other hand, as the temperature decreased from August 2015 to March 2016, the coverage of sea ice increased up to almost the same coverage of March 2015.
As shown in Figure 13a, the change of sea ice coverage mainly happened in the area of A, B and C, which are marked with blue rectangles. Sea ice in these areas can be recognized as young ice or first-year ice which includes both closed ice and open ice, while sea ice in area D marked with magenta rectangle can be regarded as multi-year ice (only include closed ice) [48]. The results show that the sea ice detection probability of closed ice and open ice in the area of A, B and C is 100% and 97.56%, respectively, while the sea ice detection probability of D (multi-year ice) is 98.74%. Generally speaking, first-year ice is relatively smoother [49] and DDMs collected over first-year ice show less spreading, which benefits for sea ice detection. The surface of multi-year ice may be rougher, likely increasing false detection probability.
Although the proposed methods hold great potential for monitoring sea ice, a small quantity of false detection still occurs around the open drift ice and multi-year ice areas. The received signals reflected from both sea ice and water may lead to false detection. Thus, more in-depth research into the relationship between reflected signals and sea ice concentration in needed in the future. Having a good knowledge about sea ice concentration will benefit sea ice detection a lot.

4. Conclusions

In this paper, GNSS-R DDM measurements collected by the SGR-ReSI on board the TDS-1 satellite are used for distinguishing between sea ice and water. A dDDM observable is proposed for the first time and two methods (PS-D and PN-D) based on power-summation and pixel-number of dDDMs are proposed to detect sea ice. Meanwhile, several data processing schemes including noise subtraction, incoherent averaging and overall normalization are employed to improve significantly the performance of proposed sea ice detection approaches. The feasibility of PS-D and PN-D methods is validated with ground-truth sea ice edge data provided by OSI SAF and the results indicate the good performance of PS-D and PN-D with a detection probability of 99.72% and 99.69% respectively and a probability of false detection of 0.28% and 0.31% respectively. The proposed methods considerably outperform the existing methods with an improved detection probability of more than 1.5% and significantly reduced probability of false detection.
The dDDMs of water-water transition show larger spreading characteristics than those of ice-ice transition due to the larger variability of reflected surface. This characteristic allows us to distinguish water-water from ice-ice transition, while it does not allow us to model the relationship between received signals and the roughness as the TDS-1 data was not calibrated and the instrument onboard the TDS-1 satellite was not optimized for altimetry purposes.
In the future, it is desirable to further reduce the sea ice false detection probability to obtain more accurate sea ice coverage. Open ice affects sea ice detection and new techniques are needed to handle this issue. Moreover, data coverage rate of TDS-1 is rather limited but CYGNSS launched in 2016 and the potential GEROS-ISS [50] mission will provide much wider coverage.

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China under Grant 41574031. The authors would like to thank the China Scholarship Council (CSC) for funding and UK TechDemoSat-1 team for providing TDS-1 DDM data.

Author Contributions

All authors have made significant contributions to this manuscript. Yongchao Zhu and Jingui Zou proposed the idea, analyzed the data and wrote the paper. Kegen Yu and Jens Wickert provided valuable advices and revised the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hamer, M. The polar road to riches. New Sci. 2006, 189, 24–25. [Google Scholar]
  2. Galley, R.J.; Else, B.G.T.; Prinsenberg, S.J.; Babb, D.; Barber, D.G. Summer Sea Ice Concentration, Motion, and Thickness near Areas of Proposed Offshore Oil and Gas Development in the Canadian Beaufort Sea–2009. Arctic 2013, 66, 105–116. [Google Scholar] [CrossRef]
  3. Francis, J.A.; Chan, W.; Leathers, D.J.; Miller, J.R.; Veron, D.E. Winter Northern Hemisphere weather patterns remember summer Arctic sea-ice extent. Geophys. Res. Lett. 2009, 36. [Google Scholar] [CrossRef]
  4. Martin-Neira, M.; Caparrini, M.; Font-Rosselo, J.; Lannelongue, S.; Vallmitjana, C.S. The PARIS concept: An experimental demonstration of sea surface altimetry using GPS reflected signals. IEEE Trans. Geosci. Remote Sens. 2001, 39, 142–150. [Google Scholar] [CrossRef]
  5. Maritin-Neira, M. A passive reflectometry and interferometry system (PARIS): Application to ocean altimetry. ESA J. 1993, 17, 331–355. [Google Scholar]
  6. Garrison, J.L.; Katzberg, S.J.; Hill, M.I. Effect of sea roughness on bistatically scattered range coded signals from the Global Positioning System. Geophys. Res. Lett. 1998, 25, 2257–2260. [Google Scholar] [CrossRef]
  7. Gleason, S.; Hodgart, S.; Sun, Y.P.; Gommenginger, C.; Mackin, S.; Adjrad, M.; Unwin, M. Detection and processing of bistatically reflected GPS signals from low earth orbit for the purpose of ocean remote sensing. IEEE Trans. Geosci. Remote Sens. 2005, 43, 1229–1241. [Google Scholar] [CrossRef]
  8. Komjathy, A.; Zavorotny, V.U.; Axelrad, P.; Born, G.H.; Garrison, J.L. GPS signal scattering from sea surface: Wind speed retrieval using experimental data and theoretical model. Remote Sens. Environ. 2000, 73, 162–174. [Google Scholar] [CrossRef]
  9. Zavorotny, V.U.; Voronovich, A.G. Scattering of GPS signals from the ocean with wind remote sensing application. IEEE Trans. Geosci. Remote Sens. 2000, 38, 951–964. [Google Scholar] [CrossRef]
  10. Garrison, J.L.; Komjathy, A.; Zavorotny, V.U.; Katzberg, S.J. Wind speed measurement using forward scattered GPS signals. IEEE Trans. Geosci. Remote Sens. 2002, 40, 50–65. [Google Scholar] [CrossRef]
  11. Komjathy, A.; Armatys, M.; Masters, D.; Axelrad, P.; Zavorotny, V.; Katzberg, S. Retrieval of ocean surface wind speed and wind direction using reflected GPS signals. J. Atmos. Ocean. Technol. 2004, 21, 515–526. [Google Scholar] [CrossRef]
  12. Ji, W.; Xiu, C.; Li, W.; Wang, L. Ocean surface target detection and positioning using the spaceborne GNSS-R Delay-Doppler maps. In Proceedings of the IEEE International Symposium on Geoscience and Remote Sensing (IGARSS), Quebec City, QC, Canada, 13–18 July 2014; pp. 3806–3809. [Google Scholar]
  13. Rodriguez-Alvarez, N.; Bosch-Lluis, X.; Camps, A.; Vall-llossera, M.; Valencia, E.; Fernando Marchan-Hernandez, J.; Ramos-Perez, I. Soil Moisture Retrieval Using GNSS-R Techniques: Experimental Results Over a Bare Soil Field. IEEE Trans. Geosci. Remote Sens. 2009, 47, 3616–3624. [Google Scholar] [CrossRef]
  14. Egido, A.; Paloscia, S.; Motte, E.; Guerriero, L.; Pierdicca, N.; Caparrini, M.; Santi, E.; Fontanelli, G.; Floury, N. Airborne GNSS-R Polarimetric Measurements for Soil Moisture and Above-Ground Biomass Estimation. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 1522–1532. [Google Scholar] [CrossRef]
  15. Alonso-Arroyo, A.; Camps, A.; Aguasca, A.; Forte, G.; Monerris, A.; Ruediger, C.; Walker, J.P.; Park, H.; Pascual, D.; Onrubia, R. Improving the Accuracy of Soil Moisture Retrievals Using the Phase Difference of the Dual-Polarization GNSS-R Interference Patterns. IEEE Geosci. Remote Sens. 2014, 11, 2090–2094. [Google Scholar] [CrossRef]
  16. Roussel, N.; Frappart, F.; Ramillien, G.; Darrozes, J.; Baup, F.; Lestarquit, L.; Ha, M.C. Detection of Soil Moisture Variations Using GPS and GLONASS SNR Data for Elevation Angles Ranging From 2 degrees to 70 degrees. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 4781–4794. [Google Scholar] [CrossRef]
  17. Small, E.E.; Larson, K.M.; Chew, C.C.; Dong, J.; Ochsner, T.E. Validation of GPS-IR Soil Moisture Retrievals: Comparison of Different Algorithms to Remove Vegetation Effects. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 4759–4770. [Google Scholar] [CrossRef]
  18. Ferrazzoli, P.; Guerriero, L.; Pierdicca, N.; Rahmoune, R. Forest biomass monitoring with GNSS-R: Theoretical simulations. Adv. Space Res. 2011, 47, 1823–1832. [Google Scholar] [CrossRef]
  19. Carreno-Luengo, H.; Amezaga, A.; Vidal, D.; Olive, R.; Fran Munoz, J.; Camps, A. First Polarimetric GNSS-R Measurements from a Stratospheric Flight over Boreal Forests. Remote Sens. 2015, 7, 13120–13138. [Google Scholar] [CrossRef]
  20. Valencia, E.; Camps, A.; Rodriguez-Alvarez, N.; Park, H.; Ramos-Perez, I. Using GNSS-R Imaging of the Ocean Surface for Oil Slick Detection. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2013, 6, 217–223. [Google Scholar] [CrossRef]
  21. Yu, K. Tsunami-Wave Parameter Estimation Using GNSS-Based Sea Surface Height Measurement. IEEE Trans. Geosci. Remote Sens. 2015, 53, 2603–2611. [Google Scholar] [CrossRef]
  22. Yu, K. Weak Tsunami Detection Using GNSS-R-Based Sea Surface Height Measurement. IEEE Trans. Geosci. Remote Sens. 2016, 54, 1363–1375. [Google Scholar] [CrossRef]
  23. Clarizia, M.P.; Ruf, C.; Cipollini, P.; Zuffada, C. First spaceborne observation of sea surface height using GPS-Reflectometry. Geophys. Res. Lett. 2016, 43, 767–774. [Google Scholar] [CrossRef]
  24. Roussel, N.; Ramillien, G.; Frappart, F.; Darrozes, J.; Gay, A.; Biancale, R.; Striebig, N.; Hanquiez, V.; Bertin, X.; Allain, D. Sea level monitoring and sea state estimate using a single geodetic receiver. Remote Sens. Environ. 2015, 171, 261–277. [Google Scholar] [CrossRef]
  25. Semmling, A.M.; Beckheinrich, J.; Wickert, J.; Beyerle, G.; Schoen, S.; Fabra, F.; Pflug, H.; He, K.; Schwabe, J.; Scheinert, M. Sea surface topography retrieved from GNSS reflectometry phase data of the GEOHALO flight mission. Geophys. Res. Lett. 2014, 41, 954–960. [Google Scholar] [CrossRef]
  26. Semmling, A.M.; Leister, V.; Saynisch, J.; Zus, F.; Heise, S.; Wickert, J. A Phase-Altimetric Simulator: Studying the Sensitivity of Earth-Reflected GNSS Signals to Ocean Topography. IEEE Trans. Geosci. Remote Sens. 2016, 54, 6791–6802. [Google Scholar] [CrossRef]
  27. Lestarquit, L.; Peyrezabes, M.; Darrozes, J.; Motte, E.; Roussel, N.; Wautelet, G.; Frappart, F.; Ramillien, G.; Biancale, R.; Zribi, M. Reflectometry With an Open-Source Software GNSS Receiver: Use Case With Carrier Phase Altimetry. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 4843–4853. [Google Scholar] [CrossRef]
  28. Yu, K.; Ban, W.; Zhang, X.; Yu, X. Snow Depth Estimation Based on Multipath Phase Combination of GPS Triple-Frequency Signals. IEEE Trans. Geosci. Remote Sens. 2015, 53, 5100–5109. [Google Scholar] [CrossRef]
  29. Zhang, Z.; Liu, S.; Wei, J.; Xu, J.; Guo, W.; Bao, W.; Jiang, Z. Mass Change of Glaciers in Muztag Ata-Kongur Tagh, Eastern Pamir, China from 1971/76 to 2013/14 as Derived from Remote Sensing Data. PLoS ONE 2016, 11, e0147327. [Google Scholar] [CrossRef] [PubMed]
  30. Wiehl, M.; Legresy, B.; Dietrich, R. Potential of reflected GNSS signals for ice sheet remote sensing. Prog. Electromagn. Res. 2003, 40, 177–205. [Google Scholar] [CrossRef]
  31. Rivas, M.B.; Maslanik, J.A.; Axelrad, P. Bistatic Scattering of GPS Signals Off Arctic Sea Ice. IEEE Trans. Geosci. Remote Sens. 2010, 48, 1548–1553. [Google Scholar] [CrossRef]
  32. Gleason, S. Towards Sea Ice Remote Sensing with Space Detected GPS Signals: Demonstration of Technical Feasibility and Initial Consistency Check Using Low Resolution Sea Ice Information. Remote Sens. 2010, 2, 2017–2039. [Google Scholar] [CrossRef]
  33. Zhang, Y.; Meng, W.; Gu, Q.; Han, Y.; Hong, Z.; Cao, Y.; Xia, Q.; Wang, W. Detection of Bohai Bay Sea Ice Using GPS-Reflected Signals. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 39–46. [Google Scholar]
  34. Jales, P.; Unwinn, M. MERRByS Product Manual: GNSS Reflectometry on TDS-1 with the SGR-ReSI. Available online: http://www.merrbys.co.uk:8080/CatalogueData/Documents/MERRByS%20Product%20Manual%20V2.pdf (assessed on 10 July 2017).
  35. Foti, G.; Gommenginger, C.; Jales, P.; Unwin, M.; Shaw, A.; Robertson, C.; Rosello, J. Spaceborne GNSS reflectometry for ocean winds: First results from the UK TechDemoSat-1 mission. Geophys. Res. Lett. 2015, 42, 5435–5441. [Google Scholar] [CrossRef]
  36. Soisuvarn, S.; Jelenak, Z.; Said, F.; Chang, P.S.; Egido, A. The GNSS Reflectometry Response to the Ocean Surface Winds and Waves. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 4678–4699. [Google Scholar] [CrossRef]
  37. Yan, Q.; Huang, W. Tsunami Detection and Parameter Estimation from GNSS-R Delay-Doppler Map. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 4650–4659. [Google Scholar] [CrossRef]
  38. Chew, C.; Shah, R.; Zuffada, C.; Hajj, G.; Masters, D.; Mannucci, A.J. Demonstrating soil moisture remote sensing with observations from the UK TechDemoSat-1 satellite mission. Geophys. Res. Lett. 2016, 43, 3317–3324. [Google Scholar] [CrossRef]
  39. Sarangi, R.K.; Devi, K.N. Space-based observation of chlorophyll, sea surface temperature, nitrate, and sea surface height anomaly over the Bay of Bengal and Arabian Sea. Adv. Space Res. 2017, 59, 33–44. [Google Scholar] [CrossRef]
  40. Yan, Q.; Huang, W. Spaceborne GNSS-R Sea Ice Detection Using Delay-Doppler Maps: First Results from the UK TechDemoSat-1 Mission. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 4795–4801. [Google Scholar] [CrossRef]
  41. Schiavulli, D.; Frappart, F.; Ramillien, G.; Darrozes, J.; Nunziata, F.; Migliaccio, M. Observing Sea/Ice Transition Using Radar Images Generated From TechDemoSat-1 Delay Doppler Maps. IEEE Geosci. Remote Sens. 2017, 14, 734–738. [Google Scholar] [CrossRef]
  42. Alonso-Arroyo, A.; Zavorotny, V.U.; Camps, A. Sea Ice Detection Using U.K. TDS-1 GNSS-R Data. IEEE Trans. Geosci. Remote Sens. 2017. [Google Scholar] [CrossRef]
  43. Jales, P.; Unwin, M. Mission Description—GNSS Reflectometry on TDS-1 with the SGR-ReSI. Available online: http://www.merrbys.co.uk:8080/CatalogueData/Documents/TDS-1%20GNSS-R%20Mission%20Description.pdf (assessed on 10 July 2017).
  44. Rodriguez-Alvarez, N.; Akos, D.M.; Zavorotny, V.U.; Smith, J.A.; Camps, A.; Fairall, C.W. Airborne GNSS-R Wind Retrievals Using Delay-Doppler Maps. IEEE Trans. Geosci. Remote Sens. 2013, 51, 626–641. [Google Scholar] [CrossRef]
  45. Marchan-Hernandez, J.F.; Rodriguez-Alvarez, N.; Camps, A.; Bosch-Lluis, X.; Ramos-Perez, I.; Valencia, E. Correction of the Sea State Impact in the L-Band Brightness Temperature by Means of Delay-Doppler Maps of Global Navigation Satellite Signals Reflected Over the Sea Surface. IEEE Trans. Geosci. Remote Sens. 2008, 46, 2914–2923. [Google Scholar] [CrossRef]
  46. Zhu, Y.; Yu, K.; Zou, J.; Wickert, J. Sea ice detection using GNSS-R Delay-Doppler Maps from UK TechDemoSat-1. In Proceedings of the IEEE International Symposium on Geoscience and Remote Sensing Symposium, Fort Worth, TX, USA, 23–28 July 2017. accepted. [Google Scholar]
  47. Aaboe, S.; Breivik, L.-A.; Sorensen, A.; Eastwood, S.; Lavergne, T. Global Sea Ice Edge and Type Product User’s Manual. Available online: http://osisaf.met.no/docs/osisaf_cdop2_ss2_pum_sea-ice-edge-type_v2p1.pdf (accessed on 10 July 2014).
  48. Aaboe, S.; Breivik, L.-A.; Sorensen, A.; Eastwood, S. Algorithm Theoretical Basis Document for the OSI SAF Global Sea Ice Edge and Type Product. Available online: http://osisaf.met.no/docs/osisaf_cdop2_ss2_atbd_sea-ice-edge_type_v1p2.pdf (accessed on 10 July 2017).
  49. Nolin, A.W.; Fetterer, F.M.; Scambos, T.A. Surface roughness characterizations of sea ice and ice sheets: Case studies with MISR data. IEEE Trans. Geosci. Remote Sens. 2002, 40, 1605–1615. [Google Scholar] [CrossRef]
  50. Wickert, J.; Cardellach, E.; Martin-Neira, M.; Bandeiras, J.; Bertino, L.; Andersen, O.B.; Camps, A.; Catarino, N.; Chapron, B.; Fabra, F.; et al. GEROS-ISS: GNSS REflectometry, Radio Occultation, and Scatterometry Onboard the International Space Station. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 4552–4581. [Google Scholar] [CrossRef]
Figure 1. Global Navigation Satellite System Reflectometry (GNSS-R) Delay-Doppler Maps (DDMs) collected over (a) sea ice field (close to a specular reflection point) and (b) sea water field obtained from the Space GNSS Receiver Remote Sensing Instrument (SGR-ReSI) receiver onboard the TDS-1 satellite on 21 January 2016.
Figure 1. Global Navigation Satellite System Reflectometry (GNSS-R) Delay-Doppler Maps (DDMs) collected over (a) sea ice field (close to a specular reflection point) and (b) sea water field obtained from the Space GNSS Receiver Remote Sensing Instrument (SGR-ReSI) receiver onboard the TDS-1 satellite on 21 January 2016.
Sensors 17 01614 g001
Figure 2. GNSS-R DDM collected over sea water surface and affected by noise. DDM n o i s e denoted as noise floor is the average value of grids within the white box.
Figure 2. GNSS-R DDM collected over sea water surface and affected by noise. DDM n o i s e denoted as noise floor is the average value of grids within the white box.
Sensors 17 01614 g002
Figure 3. Left plot: Processed GNSS-R DDM collected over sea ice field (middle specular point location: 5°40′23.14″ W, 77°52′1.62″ N, PRN 24) and right plot: Processed GNSS-R DDM collected over sea water field (middle specular point location: 17°44′41.03″ E, 73°42′20.78″ N, PRN 24) by the SGR-ReSI receiver onboard the TDS-1 satellite on 19 February 2015.
Figure 3. Left plot: Processed GNSS-R DDM collected over sea ice field (middle specular point location: 5°40′23.14″ W, 77°52′1.62″ N, PRN 24) and right plot: Processed GNSS-R DDM collected over sea water field (middle specular point location: 17°44′41.03″ E, 73°42′20.78″ N, PRN 24) by the SGR-ReSI receiver onboard the TDS-1 satellite on 19 February 2015.
Sensors 17 01614 g003
Figure 4. Normalized dDDMs for four different observation situations: (a) water-water transition, (b) water-ice transition, (c) ice-ice transition and (d) ice-water transition.
Figure 4. Normalized dDDMs for four different observation situations: (a) water-water transition, (b) water-ice transition, (c) ice-ice transition and (d) ice-water transition.
Sensors 17 01614 g004
Figure 5. Four scenarios for one dataset which includes: (a) DDMs only collected over sea ice field, (b) DDMs only collected over sea water field, (c,d) DDMs collected over both sea ice and water field.
Figure 5. Four scenarios for one dataset which includes: (a) DDMs only collected over sea ice field, (b) DDMs only collected over sea water field, (c,d) DDMs collected over both sea ice and water field.
Sensors 17 01614 g005
Figure 6. Flowchart of threshold determination for sea ice detection scheme. Six datasets over the Arctic region from TDS-1 are used as training data and sea ice edge data from OSI SAF is used as ground-truth data for determining thresholds.
Figure 6. Flowchart of threshold determination for sea ice detection scheme. Six datasets over the Arctic region from TDS-1 are used as training data and sea ice edge data from OSI SAF is used as ground-truth data for determining thresholds.
Sensors 17 01614 g006
Figure 7. GNSS-R Specular Points (SP) ground-track (green and red line) of TDS-1 signal related to the tested dataset collected over Greenland Sea on 26 October 2015 at 17:29 UTC to 17:31 UTC. The track goes from the East to the West. dDDMs corresponding to SP ground-track in red line are presented in Figure 8. The location of DDM 45 (3°23’35.44” E, 77°57’5.44” N) and DDM 75 (4°26’20.72” W, 78°14’11.77” N) are marked with a blue line and purple line, respectively.
Figure 7. GNSS-R Specular Points (SP) ground-track (green and red line) of TDS-1 signal related to the tested dataset collected over Greenland Sea on 26 October 2015 at 17:29 UTC to 17:31 UTC. The track goes from the East to the West. dDDMs corresponding to SP ground-track in red line are presented in Figure 8. The location of DDM 45 (3°23’35.44” E, 77°57’5.44” N) and DDM 75 (4°26’20.72” W, 78°14’11.77” N) are marked with a blue line and purple line, respectively.
Sensors 17 01614 g007
Figure 8. Examples of processed dDDMs (dDDM 45 to dDDM 74) obtained from DDM 45 to DDM 75 of water-water transition, water-ice transition and ice-ice transition.
Figure 8. Examples of processed dDDMs (dDDM 45 to dDDM 74) obtained from DDM 45 to DDM 75 of water-water transition, water-ice transition and ice-ice transition.
Sensors 17 01614 g008
Figure 9. (a) PS-D method results and (b) PN-D method results for five different thresholds ( D D M T = 0.36, D D M T = 0.38, D D M T = 0.40, D D M T = 0.42 and D D M T = 0.44).
Figure 9. (a) PS-D method results and (b) PN-D method results for five different thresholds ( D D M T = 0.36, D D M T = 0.38, D D M T = 0.40, D D M T = 0.42 and D D M T = 0.44).
Sensors 17 01614 g009
Figure 10. Significant Wave Height data from the European Centre for Medium-Range Weather Forecasts (HW-ECMWF) over high latitude area of the north hemisphere (Latitude: 56° N–84° N, Longitude: 44° W–40° E) on 26 October 2015. The red line is SP ground-track of DDMs (DDM 1 to DDM 57). There is no valid wave height from DDM 58 to DDM 131.
Figure 10. Significant Wave Height data from the European Centre for Medium-Range Weather Forecasts (HW-ECMWF) over high latitude area of the north hemisphere (Latitude: 56° N–84° N, Longitude: 44° W–40° E) on 26 October 2015. The red line is SP ground-track of DDMs (DDM 1 to DDM 57). There is no valid wave height from DDM 58 to DDM 131.
Sensors 17 01614 g010
Figure 11. The relationship between PS and WH-ECMWF of dDDM 1 to dDDM 56. Four groups of PS results ( D D M T = 0.36, D D M T = 0.38, D D M T = 0.40, D D M T = 0.42) is presented. The values of WH are scaled down as we just concentrate on the variability.
Figure 11. The relationship between PS and WH-ECMWF of dDDM 1 to dDDM 56. Four groups of PS results ( D D M T = 0.36, D D M T = 0.38, D D M T = 0.40, D D M T = 0.42) is presented. The values of WH are scaled down as we just concentrate on the variability.
Sensors 17 01614 g011
Figure 12. Sea ice detection probability of four different detection methods. PS-D (blue bar) and PN-D (cyan bar) based methods show better performance than PS-N (yellow bar), PN-N (magenta bar) and MF (red bar) based approaches.
Figure 12. Sea ice detection probability of four different detection methods. PS-D (blue bar) and PN-D (cyan bar) based methods show better performance than PS-N (yellow bar), PN-N (magenta bar) and MF (red bar) based approaches.
Sensors 17 01614 g012
Figure 13. Sea ice detection results based on PS-D using differential and normalized DDM observables in six different days: (a) 24 March 2015, (b) 4 June 2015, (c) 23 August 2015, (d) 26 October 2015, (e) 21 January 2016, (f) 26 March 2016. Lines in green and red represent detected sea ice and water respectively. Three areas marked with blue rectangles are defined as A, B and C respectively. The area marked with magenta rectangle is regarded as D. The white and gray area depicts the ground-truth sea ice coverage (sea ice edge product of the EUMETSAT Ocean and Sea Ice Satellite Application Facility, OSI SAF, www.osi-saf.org).
Figure 13. Sea ice detection results based on PS-D using differential and normalized DDM observables in six different days: (a) 24 March 2015, (b) 4 June 2015, (c) 23 August 2015, (d) 26 October 2015, (e) 21 January 2016, (f) 26 March 2016. Lines in green and red represent detected sea ice and water respectively. Three areas marked with blue rectangles are defined as A, B and C respectively. The area marked with magenta rectangle is regarded as D. The white and gray area depicts the ground-truth sea ice coverage (sea ice edge product of the EUMETSAT Ocean and Sea Ice Satellite Application Facility, OSI SAF, www.osi-saf.org).
Sensors 17 01614 g013
Table 1. Information for the tested dataset.
Table 1. Information for the tested dataset.
Start time17:29:01 26 October 2015
End time17:31:11 26 October 2015
Start point13°58′48.05″ E 77°09′55.70″ N
End point19°24′48.51″ W 78°10′31.35″ N
No. Of DDM131
Ground-track length793 km
Table 2. Probability of false detection for four observables using GNSS-R DDM.
Table 2. Probability of false detection for four observables using GNSS-R DDM.
Date of DatasetNumber of DatasetNumber of Tested DDMPS-DPN-DPS-NPN-NMF
24 March 20152228150.280.363.203.481.88
4 June 20152529060.340.413.993.921.72
23 August 20152026680.220.303.673.371.95
26 October 20152234110.290.233.283.341.91
21 January 20162127640.220.294.123.981.77
26 March 20162027820.290.293.523.241.83
Total130173460.280.313.623.551.84

Share and Cite

MDPI and ACS Style

Zhu, Y.; Yu, K.; Zou, J.; Wickert, J. Sea Ice Detection Based on Differential Delay-Doppler Maps from UK TechDemoSat-1. Sensors 2017, 17, 1614. https://doi.org/10.3390/s17071614

AMA Style

Zhu Y, Yu K, Zou J, Wickert J. Sea Ice Detection Based on Differential Delay-Doppler Maps from UK TechDemoSat-1. Sensors. 2017; 17(7):1614. https://doi.org/10.3390/s17071614

Chicago/Turabian Style

Zhu, Yongchao, Kegen Yu, Jingui Zou, and Jens Wickert. 2017. "Sea Ice Detection Based on Differential Delay-Doppler Maps from UK TechDemoSat-1" Sensors 17, no. 7: 1614. https://doi.org/10.3390/s17071614

APA Style

Zhu, Y., Yu, K., Zou, J., & Wickert, J. (2017). Sea Ice Detection Based on Differential Delay-Doppler Maps from UK TechDemoSat-1. Sensors, 17(7), 1614. https://doi.org/10.3390/s17071614

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