# C-Band Polarimetric Coherences and Ratios for Discriminating Sea Ice Roughness

Authors: Mukesh Gupta; Randall K. Scharien; David G. Barber
Journal: International Journal of Oceanography
 (2013)

---

## Abstract

The rapid decline of sea ice in the Arctic has resulted in a variable sea ice roughness that necessitates improved methods for efficient observation using high-resolution spaceborne radar. The utility of C-band polarimetric backscatter, coherences, and ratios as a discriminator of ice surface roughness is evaluated. An existing one-dimensional backscatter model has been modified to two-dimensions (2D) by considering deviation in the orientation (i.e., the slopes) in azimuth and range direction of surface roughness simultaneously as an improvement in the model. It is shown theoretically that the circular coherence (ρRRLL) decreases exponentially with increasing surface roughness. The crosspolarized coherence (ρHHVH) is found to be less sensitive to surface roughness, whereas the copolarized coherence (ρVVHH) decreases at far-range incidence angles for all ice types. A complete validation of the adapted 2D model using direct measurements of surface roughness is suggested as an avenue for further research.


---

## Body

## 1. Introduction

Arctic sea ice is going through a rapid decline [1, 2]. Thinner first-year ice (FYI) is replacing multiyear ice, leaving an ice cover, which is more sensitive to deformation and changes in atmospheric and ocean forcing. Increased open water and marginal ice zones (MIZs), due to the enhanced mobility of a relatively thinned pack ice, are further susceptible to increases in surface roughness and greater surface roughness variability [3]. Greater surface roughness in the MIZ is of importance due to higher rates of heat flux [4] and momentum [5] exchanges occurring across the ocean-sea ice-atmosphere interface, greater biological productivity [6], and potential limitations imposed on ship navigation. Although the literature contains information on how the MIZ responds to wind and wave forces, it is necessary to investigate the electromagnetic (EM) response of the MIZ to facilitate satellite-based observations. Satellite-based observation is necessary due to the scarcity of surface observations in a MIZ, as well as the difficulties in collecting physical measurements due to the instability and roughness of the ice floes.The use of polarimetric synthetic aperture radar (pol-SAR) represents a promising approach for satellite-based monitoring of surface roughness and, concurrently, discriminating sea ice types within a MIZ. A pol-SAR records the amplitude and phase information of backscattered energy for four transmit-receive polarizations (HH, HV, VH, and VV), thereby facilitating the derivation of the full polarimetric response of the target. It is recognizable that the diversity in polarization achievable by pol-SARs or even by dual-polarization SAR systems provides more complete inference of target features (e.g., sea ice) than conventional, single channel SARs. Furthermore, recently launched pol-SARs are capable of higher spatial resolution (<10 m) imaging, leading to enhanced potential for monitoring complex ice environments.Discrimination of ice types using SAR has been conventionally achieved by utilizing different combinations of linearly polarized backscattering coefficients [7–9]. Multiyear ice, smooth FYI, rough FYI, and new ice/open water in the Beaufort and Chukchi Seas during March have been identified using a single-polarization SAR image intensity-based classification scheme [9] while others used single-polarization SAR image texture analysis to discriminate new ice, FYI, and multiyear ice during the month of March in the Beaufort Sea and the Mould Bay, respectively [10, 11]. Dual copolarized backscattering coefficient differences in HH and VV have been used to discriminate FYI, multiyear ice, and lead areas in the Beaufort Sea during March [12]. However, the complexities in polarimetric signatures associated with the dynamic mixture of surface roughness and ice type conditions in an MIZ during fall freeze-up remain to be examined. Such an examination requires utilizing polarimetric radar backscatter, so that the material (dielectric) and geometrical properties of the surface, which influence backscatter, may be individually assessed.In this study, ship-based observations of co- (linear) and crosspolarized backscatter, circular polarimetric coherences (ρVVHH, ρHHVH, and ρRRLL, resp.), and copolarized and crosspolarized polarization ratios (γco and γcross, resp.), are used to evaluate their utility for ice surface discrimination capabilities using a polarimetric radar operating in C-band (5.5 GHz). Characteristics of these polarimetric parameters for a variety of ice types in an MIZ during fall freeze-up are assessed with the following objectives: (1)
to investigate the performance of polarimetricρRRLL for sea ice surface roughness discrimination by adapting the one-dimensional backscatter model of [13] to two dimensions and introducing roughness as deviations in range and azimuth directions,(2)
to evaluate the utility of C-band polarimetric backscatter, coherences, and polarization ratios as a discriminator of surface roughness or ice type in a MIZ during fall freeze-up.


## 2. Methodology

### 2.1. Study Area

The study area is located in the southeastern Beaufort Sea and Amundsen Gulf regions in the western Canadian Arctic (Figure1). The seasonal Cape Bathurst Polynya forms in the region and hosts a number of flaw leads during the winter [14]. During fall freeze-up, this area contains a variable mix of ice types under various stages of formation, for example, new ice, pancake ice, frost flowers, deformed ice, gray ice, and nilas (Figure 2). The photographs in Figures 2(a), 2(b), 2(c), and 2(e) were taken at an oblique angle from the port side of the Canadian Research Icebreaker Amundsen at approximately eight meters height using a handheld digital camera after a given scatterometer scan; Figure 2(d) was taken at nadir angle on the ice floe at about one meter height. In the present study, thin FYI types are considered (first stage: 30–50 cm—as per World Meteorological Organization nomenclature), which include snow-covered FYI, pancake ice, frost flowers, and deformed FYI located within the MIZ. Data described in the following sections were acquired as part of the Circumpolar Flaw Lead (CFL) System Study project of the International Polar Year (IPY) 2007-08 over the period October 2007–August 2008 [15]. Ancillary meteorological data were collected through a ship-based AXYS Technologies Inc. (Sydney, BC, Canada) Automatic Voluntary Observing Ships (AVOS) system. This system was mounted approximately 20 m above sea level on the wheelhouse to minimize the ship’s influence and could measure air temperature and wind speed.Figure 1
Geographic map of study area showing sampling locations.Photographs of ice types used in the study. (a) snow-covered first-year ice (SCFYI), (b) deformed first-year ice (DFYI), (c) consolidated pancake ice (PI), (d) snow-covered frost flower (SCFF), and (e) dense frost flower (DFF).
(a)
(b)
(c)
(d)
(e)


### 2.2. Theoretical Formulation

Sea ice is a distributed radar target, and the conditions of stationarity and homogeneity seldom hold for dynamically changing ice in a MIZ. The radar backscattering is therefore analyzed using temporally and spatially varying stochastic processes. Backscatter from sea ice is incoherent and either partially or completely polarized, as described by the polarimetric covariance matrix. The electric field vector of an incident (i) and scattered (s) EM wave can be given by
(1)Ei=EHih^+EViv^,Es=EHsh^+EVsv^,
where H and V represent horizontal and vertical polarizations respectively. h^ and v^ are the unit vectors in the horizontal and vertical directions of polarization, respectively. The incident (i) and scattered field (s) can be either H or V. The scattered electric field is related to the incident electric field by the scattering matrix, S defined as
(2)[EHsEVs]=e-jkrr[SHHSHVSVHSVV][EHiEVi],[SHHSHVSVHSVV]=S,
where e-jkr/r term accounts for wave propagation effects in amplitude and phase. If the orientation of a surface such as sea ice in azimuth direction is rotated by an angle, the corresponding new backscatter matrix can be constructed as provided by Lee et al. [16].The coherency matrices can be derived as copolarized (3), crosspolarized (4) and circular (RRLL: right-right left-left rotation of the electric field vector about the line of sight) (5) coherences in magnitude form [13, 16] as (for derivation of ρRRLL, see Appendix A),
(3)ρVVHH=〈|SVVSHH*|〉〈|SVV|2〉〈|SHH|2〉,(4)ρHHVH=〈|SHHSVH*|〉〈|SHH|2〉〈|SVH|2〉,(5)ρRRLL=〈|SHH-SVV|2〉-4〈|SHV|2〉〈|SHH-SVV|2〉+4〈|SHV|2〉,(6)γco=SVVSHH,(7)γcross=SHVSHH,
where S is the complex scattering matrix; an asterisk (*) represents the complex conjugate. The brackets 〈·〉 represent ensemble averages of the observed data. There were approximately 34 pulses sent per incidence angle. An ensemble average was performed on those 34 pulses. Raw data were processed into range profiles and were averaged in the azimuth for each measured incidence angle. Polarimetric ratios γco and γcross are simply power ratios of backscattered energy. Polarimetric coherences and polarization ratios have demonstrated utility in reducing the ambiguities caused by the nonlinearity between system response and target properties. Regarding Arctic sea ice, some literature is available on the use of ρVVHH, ρRRLL, and γco at different EM frequencies. C-band backscatter coefficients (HH, HV, and VV) and ρVVHH have been used to characterize various FYI types (compressed, rubble and ridge, and smooth) and multiyear ice [17]. Thin sea ice has been effectively discriminated from FYI using C-band γco ratio [18]. ρVVHH and γco have been used to discriminate Arctic leads using L-band radar signatures [19]. In a similar study, Wakabayashi et al. [20] described polarimetric characteristics of different FYI types (thin ice, smooth, and rough) using L-band ρRRLL and γco and showed the utility of coherences and ratios in discriminating ice types. Nakamura et al. [21] discriminated ice surface using γco ratio in an observational study of lake ice using airborne L- and X-band SAR. These studies lack a holistic overview of the utility of different polarimetric coherences and ratios to discriminate thin FYI types in a MIZ.


### 2.3. Active Microwave Backscattering Data

C-band polarimetric backscattering data were collected using a completely stationary ship-mounted scatterometer system developed by ProSensing Inc., (Amherst, MA, USA) and mounted 7.56 m above the mean sea level on the port side of theAmundsen (Table 1). The system acquires backscatter and phase data in terms of the combinations of linear transmit-receive polarization combinations, HH, HV, VH, and VV at incidence angles 20–60° (5° increments) over a 60° azimuth range. The calibration of the instrument was performed through the methods given elsewhere [22, 23]. Polarimetric backscattering data were collected from homogeneous samples of snow-covered (dry and fresh) first-year ice (SCFYI), deformed FYI (DFYI), consolidated pancake ice (PI), snow-covered frost flowers (SCFF), and dense frost flowers (DFF) on different dates during November 2007. Data from each ice type sample comprised three to four contiguous scatterometer scans, which took up to 35 minutes to complete. The scatterometer had a footprint of 1.1 m2 in the range direction at a 45° incidence angle [23] with the footprint increasing in size with incidence angle [22].Table 1
Technical properties and specifications of C-band scatterometer.
RF output frequency
5.25–5.75 GHz
Transmit power at bulkhead connector
12 dBm
Antenna diameter
0.61 m
Transmit bandwidth
500 MHz
Antenna beamwidth
5.5°
Antenna gain
28 dB, nominal
Crosspolarization isolation
>30 dB, measured at the peak of the beam
Transmit-receive polarizations
Linear, vertical, and horizontal
Sensitivity, minimum NRCS at 15 m range
–40 dB m2/m2Towards objective 2, scan data for each ice type were grouped by incidence angle representing near (20–25°), mid (35–40°), and far (55–60°) range groupings. These groupings best represent the diversity of scattering mechanisms available across the acquired incidence angle range. In the near range, surface scattering is expected to dominate the measured C-band backscatter, while surface-volume scattering is increasingly expected to influence C-band backscatter beyond approximately 30°, that is, mid to far ranges [24]. Furthermore, combining data from adjacent incidence angles doubled the number of samples from 8 to 18 depending on ice type, although at the expense of range resolution. Scatterometer data had unequal number of data points in each range group, which does not fulfil parametric ANOVA requirements for statistical significance testing. Polarimetric coherences and ratios of ice types were tested for independence from each other for each incidence angle grouping. Testing was done using the nonparametric Kruskal-Wallis statistic, with α=0.01 significance level (one tailed) used as the threshold for statistical independence.


### 2.4. Surface Roughness and Circular Coherence

In pursuit of objective 1, a polarimetric backscattering model is used which is mainly a Bragg backscattering (coherent scattering) model modified for surface roughness considering the surface slope by slightly changing the tilt of the surface from the horizontal. Microwave measurements of surface roughness using co- or crosspolarization backscattered power are most successful in flat areas. In sea ice microwave remote sensing, the dielectric constant and topography (slope in range and azimuth) are important. According to(22) in the one-dimensional scattering model of [13], the circular coherence is only sensitive to surface roughness. Surface roughness has been considered as a change in the slope of ice in azimuth and ground range directions [13, 16, 25]. This is implemented mathematically in the Bragg backscattering model by considering roughness as a depolarizer which conforms to reflection symmetry; that is, the backscattering properties are identical on either side of the plane of incidence and HV = VH [25, 26]. The distribution of azimuth slope angles θ1 is considered as one-dimensional Gaussian distributed [13].The rotation matrix [16] and the coherency matrix [13] are calculated after introducing the rotation in azimuth anticlockwise about range direction. In this case, ρRRLL is derived as [13]
(8)ρRRLL=e-8σθ12,
where σθ1 is the standard deviation of the orientation angle distribution in azimuth direction and θ1 is slope angle in azimuth direction. From (8), the ρRRLL is only dependent on the orientation of ice surface in the range direction, or the standard deviation of the orientation angle distribution (i.e., surface roughness). Here, the surface roughness is introduced through rotation by angle θ2 in the range direction anti-clockwise about azimuth direction (Figure 3). Angle θ2 is not shown in Figure 3 due to complexity of the geometry. In this case also, the corresponding distribution of shift in orientation angle is Gaussian distributed.Figure 3
Illustration of scattering plane geometry with slight deviations in the orientation angles in azimuth (θ1) and range directions (θ2: not shown), respectively, as means of two-dimensional surface roughness.The new rotation matrixU2 is given by
(9)U2=[cos2θ20-sin2θ2010sin2θ20cos2θ2].
The new averaged coherency matrix over the Gaussian distribution p(θ2) can be calculated as
(10)〈T〉θ2=[ζAμB0μB*2C000(1-ζ)A],
where μ(θ2)=∫cos2θ2p(θ2)dθ2 and ζ(θ2)=∫cos22θ2p(θ2)dθ2.  B, a part of an element of coherency matrix, is defined according to scattering matrix, S [13]. B* is the conjugate of B. Both B* and B are not used in the computation of ρRRLL.TheρRRLL can be computed as
(11)ρRRLL=T22-T33T22+T33=2C-(1-ζ(θ2))A2C+(1-ζ(θ2))A,(12)ρRRLL=4C-(1-e-8σθ22)A4C+(1-e-8σθ22)A,
where A=|SHH+SVV|2, C=(1/2)|SHH-SVV|2, and ζ(θ2)=(1/2)(1+e-8σθ22). Tij represents (i,j)th element of the matrix, 〈T〉θ2 given in (10).Given the above, theρRRLL is dependent on the standard deviation of the orientation angle distribution in range and the dielectric constant of the surface. Thus, it is shown that the new ρRRLL is exponentially changing with the change in orientation angle in the azimuth direction, but it behaves in a way given by (12) and is dependent on both surface roughness (standard deviation) and the dielectric constant (scattering matrix) of the surface when roughness in two directions is considered. In our model, when two-dimensional roughness is considered, circular coherence is observed to be sensitive to both surface roughness and dielectric constant, thus, making it difficult to differentiate roughness. 2D model being more realistic requires further considerations of separating dielectrics from roughness.Now, the slope-induced roughness is examined in the range direction only. Lee et al. [16] gave a relationship between slope in azimuth, slope in ground range, radar look angle (ϕ), and rotation in azimuth. Schuler et al. [13] expressed this relationship in terms of root mean square (rms) surface height (s) and correlation length (l), assuming that the range slope and orientation in azimuth are small perturbations around their means,
(13)ρRRLL=e-16(s2/l2sin2ϕ).
Figure 4 shows the incidence angle dependence of ρRRLL by varying the s2/l2 ratio. As the roughness increases, ρRRLL decreases. For l≫s, that is, the surface is very smooth, the maximum value of ρRRLL approaches unity. ρRRLL decreases exponentially from unity to a fixed value of s2/l2 ratio at a particular incidence angle. A rough surface yields a smaller ρRRLL, which increases with increasing radar look angle. The range of s2/l2 for the presented ice classes is expected to lie between 0.001 and 0.1 [27].Figure 4
ρ
RRLL varying with squared ratio of rms surface height and surface correlation length; ρRRLL decays exponentially; however, it decays faster at steep incidence angles.The relationship between slopes in azimuth and range direction is further demonstrated. Corresponding shifts and radar incidence angle are given by (see AppendixB),
(14)tanωtanθ1=sinθ2(tanγ·sinϕ+cosϕ)+cosθ2(-tanγ·cosϕ+sinϕ),
where tanω is azimuth slope, tanγ is range slope, θ1 and θ2 are the perturbations in orientation in azimuth and range directions, respectively, and ϕ is radar look angle. Figure 4 represents the case when orientation shift in the range direction is observed. In a sea ice remote sensing context, both surface roughness and the dielectric constant of ice affect ρRRLL when slope is changed in azimuth direction, whereas only surface roughness affects ρRRLL when slope is changed in range direction.


## 2.1. Study Area

The study area is located in the southeastern Beaufort Sea and Amundsen Gulf regions in the western Canadian Arctic (Figure1). The seasonal Cape Bathurst Polynya forms in the region and hosts a number of flaw leads during the winter [14]. During fall freeze-up, this area contains a variable mix of ice types under various stages of formation, for example, new ice, pancake ice, frost flowers, deformed ice, gray ice, and nilas (Figure 2). The photographs in Figures 2(a), 2(b), 2(c), and 2(e) were taken at an oblique angle from the port side of the Canadian Research Icebreaker Amundsen at approximately eight meters height using a handheld digital camera after a given scatterometer scan; Figure 2(d) was taken at nadir angle on the ice floe at about one meter height. In the present study, thin FYI types are considered (first stage: 30–50 cm—as per World Meteorological Organization nomenclature), which include snow-covered FYI, pancake ice, frost flowers, and deformed FYI located within the MIZ. Data described in the following sections were acquired as part of the Circumpolar Flaw Lead (CFL) System Study project of the International Polar Year (IPY) 2007-08 over the period October 2007–August 2008 [15]. Ancillary meteorological data were collected through a ship-based AXYS Technologies Inc. (Sydney, BC, Canada) Automatic Voluntary Observing Ships (AVOS) system. This system was mounted approximately 20 m above sea level on the wheelhouse to minimize the ship’s influence and could measure air temperature and wind speed.Figure 1
Geographic map of study area showing sampling locations.Photographs of ice types used in the study. (a) snow-covered first-year ice (SCFYI), (b) deformed first-year ice (DFYI), (c) consolidated pancake ice (PI), (d) snow-covered frost flower (SCFF), and (e) dense frost flower (DFF).
(a)
(b)
(c)
(d)
(e)


## 2.2. Theoretical Formulation

Sea ice is a distributed radar target, and the conditions of stationarity and homogeneity seldom hold for dynamically changing ice in a MIZ. The radar backscattering is therefore analyzed using temporally and spatially varying stochastic processes. Backscatter from sea ice is incoherent and either partially or completely polarized, as described by the polarimetric covariance matrix. The electric field vector of an incident (i) and scattered (s) EM wave can be given by
(1)Ei=EHih^+EViv^,Es=EHsh^+EVsv^,
where H and V represent horizontal and vertical polarizations respectively. h^ and v^ are the unit vectors in the horizontal and vertical directions of polarization, respectively. The incident (i) and scattered field (s) can be either H or V. The scattered electric field is related to the incident electric field by the scattering matrix, S defined as
(2)[EHsEVs]=e-jkrr[SHHSHVSVHSVV][EHiEVi],[SHHSHVSVHSVV]=S,
where e-jkr/r term accounts for wave propagation effects in amplitude and phase. If the orientation of a surface such as sea ice in azimuth direction is rotated by an angle, the corresponding new backscatter matrix can be constructed as provided by Lee et al. [16].The coherency matrices can be derived as copolarized (3), crosspolarized (4) and circular (RRLL: right-right left-left rotation of the electric field vector about the line of sight) (5) coherences in magnitude form [13, 16] as (for derivation of ρRRLL, see Appendix A),
(3)ρVVHH=〈|SVVSHH*|〉〈|SVV|2〉〈|SHH|2〉,(4)ρHHVH=〈|SHHSVH*|〉〈|SHH|2〉〈|SVH|2〉,(5)ρRRLL=〈|SHH-SVV|2〉-4〈|SHV|2〉〈|SHH-SVV|2〉+4〈|SHV|2〉,(6)γco=SVVSHH,(7)γcross=SHVSHH,
where S is the complex scattering matrix; an asterisk (*) represents the complex conjugate. The brackets 〈·〉 represent ensemble averages of the observed data. There were approximately 34 pulses sent per incidence angle. An ensemble average was performed on those 34 pulses. Raw data were processed into range profiles and were averaged in the azimuth for each measured incidence angle. Polarimetric ratios γco and γcross are simply power ratios of backscattered energy. Polarimetric coherences and polarization ratios have demonstrated utility in reducing the ambiguities caused by the nonlinearity between system response and target properties. Regarding Arctic sea ice, some literature is available on the use of ρVVHH, ρRRLL, and γco at different EM frequencies. C-band backscatter coefficients (HH, HV, and VV) and ρVVHH have been used to characterize various FYI types (compressed, rubble and ridge, and smooth) and multiyear ice [17]. Thin sea ice has been effectively discriminated from FYI using C-band γco ratio [18]. ρVVHH and γco have been used to discriminate Arctic leads using L-band radar signatures [19]. In a similar study, Wakabayashi et al. [20] described polarimetric characteristics of different FYI types (thin ice, smooth, and rough) using L-band ρRRLL and γco and showed the utility of coherences and ratios in discriminating ice types. Nakamura et al. [21] discriminated ice surface using γco ratio in an observational study of lake ice using airborne L- and X-band SAR. These studies lack a holistic overview of the utility of different polarimetric coherences and ratios to discriminate thin FYI types in a MIZ.


## 2.3. Active Microwave Backscattering Data

C-band polarimetric backscattering data were collected using a completely stationary ship-mounted scatterometer system developed by ProSensing Inc., (Amherst, MA, USA) and mounted 7.56 m above the mean sea level on the port side of theAmundsen (Table 1). The system acquires backscatter and phase data in terms of the combinations of linear transmit-receive polarization combinations, HH, HV, VH, and VV at incidence angles 20–60° (5° increments) over a 60° azimuth range. The calibration of the instrument was performed through the methods given elsewhere [22, 23]. Polarimetric backscattering data were collected from homogeneous samples of snow-covered (dry and fresh) first-year ice (SCFYI), deformed FYI (DFYI), consolidated pancake ice (PI), snow-covered frost flowers (SCFF), and dense frost flowers (DFF) on different dates during November 2007. Data from each ice type sample comprised three to four contiguous scatterometer scans, which took up to 35 minutes to complete. The scatterometer had a footprint of 1.1 m2 in the range direction at a 45° incidence angle [23] with the footprint increasing in size with incidence angle [22].Table 1
Technical properties and specifications of C-band scatterometer.
RF output frequency
5.25–5.75 GHz
Transmit power at bulkhead connector
12 dBm
Antenna diameter
0.61 m
Transmit bandwidth
500 MHz
Antenna beamwidth
5.5°
Antenna gain
28 dB, nominal
Crosspolarization isolation
>30 dB, measured at the peak of the beam
Transmit-receive polarizations
Linear, vertical, and horizontal
Sensitivity, minimum NRCS at 15 m range
–40 dB m2/m2Towards objective 2, scan data for each ice type were grouped by incidence angle representing near (20–25°), mid (35–40°), and far (55–60°) range groupings. These groupings best represent the diversity of scattering mechanisms available across the acquired incidence angle range. In the near range, surface scattering is expected to dominate the measured C-band backscatter, while surface-volume scattering is increasingly expected to influence C-band backscatter beyond approximately 30°, that is, mid to far ranges [24]. Furthermore, combining data from adjacent incidence angles doubled the number of samples from 8 to 18 depending on ice type, although at the expense of range resolution. Scatterometer data had unequal number of data points in each range group, which does not fulfil parametric ANOVA requirements for statistical significance testing. Polarimetric coherences and ratios of ice types were tested for independence from each other for each incidence angle grouping. Testing was done using the nonparametric Kruskal-Wallis statistic, with α=0.01 significance level (one tailed) used as the threshold for statistical independence.


## 2.4. Surface Roughness and Circular Coherence

In pursuit of objective 1, a polarimetric backscattering model is used which is mainly a Bragg backscattering (coherent scattering) model modified for surface roughness considering the surface slope by slightly changing the tilt of the surface from the horizontal. Microwave measurements of surface roughness using co- or crosspolarization backscattered power are most successful in flat areas. In sea ice microwave remote sensing, the dielectric constant and topography (slope in range and azimuth) are important. According to(22) in the one-dimensional scattering model of [13], the circular coherence is only sensitive to surface roughness. Surface roughness has been considered as a change in the slope of ice in azimuth and ground range directions [13, 16, 25]. This is implemented mathematically in the Bragg backscattering model by considering roughness as a depolarizer which conforms to reflection symmetry; that is, the backscattering properties are identical on either side of the plane of incidence and HV = VH [25, 26]. The distribution of azimuth slope angles θ1 is considered as one-dimensional Gaussian distributed [13].The rotation matrix [16] and the coherency matrix [13] are calculated after introducing the rotation in azimuth anticlockwise about range direction. In this case, ρRRLL is derived as [13]
(8)ρRRLL=e-8σθ12,
where σθ1 is the standard deviation of the orientation angle distribution in azimuth direction and θ1 is slope angle in azimuth direction. From (8), the ρRRLL is only dependent on the orientation of ice surface in the range direction, or the standard deviation of the orientation angle distribution (i.e., surface roughness). Here, the surface roughness is introduced through rotation by angle θ2 in the range direction anti-clockwise about azimuth direction (Figure 3). Angle θ2 is not shown in Figure 3 due to complexity of the geometry. In this case also, the corresponding distribution of shift in orientation angle is Gaussian distributed.Figure 3
Illustration of scattering plane geometry with slight deviations in the orientation angles in azimuth (θ1) and range directions (θ2: not shown), respectively, as means of two-dimensional surface roughness.The new rotation matrixU2 is given by
(9)U2=[cos2θ20-sin2θ2010sin2θ20cos2θ2].
The new averaged coherency matrix over the Gaussian distribution p(θ2) can be calculated as
(10)〈T〉θ2=[ζAμB0μB*2C000(1-ζ)A],
where μ(θ2)=∫cos2θ2p(θ2)dθ2 and ζ(θ2)=∫cos22θ2p(θ2)dθ2.  B, a part of an element of coherency matrix, is defined according to scattering matrix, S [13]. B* is the conjugate of B. Both B* and B are not used in the computation of ρRRLL.TheρRRLL can be computed as
(11)ρRRLL=T22-T33T22+T33=2C-(1-ζ(θ2))A2C+(1-ζ(θ2))A,(12)ρRRLL=4C-(1-e-8σθ22)A4C+(1-e-8σθ22)A,
where A=|SHH+SVV|2, C=(1/2)|SHH-SVV|2, and ζ(θ2)=(1/2)(1+e-8σθ22). Tij represents (i,j)th element of the matrix, 〈T〉θ2 given in (10).Given the above, theρRRLL is dependent on the standard deviation of the orientation angle distribution in range and the dielectric constant of the surface. Thus, it is shown that the new ρRRLL is exponentially changing with the change in orientation angle in the azimuth direction, but it behaves in a way given by (12) and is dependent on both surface roughness (standard deviation) and the dielectric constant (scattering matrix) of the surface when roughness in two directions is considered. In our model, when two-dimensional roughness is considered, circular coherence is observed to be sensitive to both surface roughness and dielectric constant, thus, making it difficult to differentiate roughness. 2D model being more realistic requires further considerations of separating dielectrics from roughness.Now, the slope-induced roughness is examined in the range direction only. Lee et al. [16] gave a relationship between slope in azimuth, slope in ground range, radar look angle (ϕ), and rotation in azimuth. Schuler et al. [13] expressed this relationship in terms of root mean square (rms) surface height (s) and correlation length (l), assuming that the range slope and orientation in azimuth are small perturbations around their means,
(13)ρRRLL=e-16(s2/l2sin2ϕ).
Figure 4 shows the incidence angle dependence of ρRRLL by varying the s2/l2 ratio. As the roughness increases, ρRRLL decreases. For l≫s, that is, the surface is very smooth, the maximum value of ρRRLL approaches unity. ρRRLL decreases exponentially from unity to a fixed value of s2/l2 ratio at a particular incidence angle. A rough surface yields a smaller ρRRLL, which increases with increasing radar look angle. The range of s2/l2 for the presented ice classes is expected to lie between 0.001 and 0.1 [27].Figure 4
ρ
RRLL varying with squared ratio of rms surface height and surface correlation length; ρRRLL decays exponentially; however, it decays faster at steep incidence angles.The relationship between slopes in azimuth and range direction is further demonstrated. Corresponding shifts and radar incidence angle are given by (see AppendixB),
(14)tanωtanθ1=sinθ2(tanγ·sinϕ+cosϕ)+cosθ2(-tanγ·cosϕ+sinϕ),
where tanω is azimuth slope, tanγ is range slope, θ1 and θ2 are the perturbations in orientation in azimuth and range directions, respectively, and ϕ is radar look angle. Figure 4 represents the case when orientation shift in the range direction is observed. In a sea ice remote sensing context, both surface roughness and the dielectric constant of ice affect ρRRLL when slope is changed in azimuth direction, whereas only surface roughness affects ρRRLL when slope is changed in range direction.


## 3. Field Results

### 3.1. Sea Ice Type Discrimination (Coherences and Ratios)

The date and hour of scatterometer data acquisitions corresponding to each sea ice type, as well as coincident meteorological parameters, namely, wind speed, air temperature, and relative humidity, are provided in Table2. The photographs of the selected ice samples are shown in Figure 2. With the exception of wind speed, there is negligible variation in meteorological conditions between ice type scans. As such, it is expected that between-scan, temperature-induced effects on the dielectric properties, and backscattering intensities from the different ice types are negligible.Table 2
Meteorological parameters associated with each ice type on different dates.
Sea ice type
Wind speed (m/s)
Air temperature (°C)
Relative humidity %
Nov. 15, 2007(Stn. 1117, 1400 hrs)
SCFYI
14.4
-
16.0
85
Nov. 19, 2007(Stn. 1100, 0030 hrs)
DFYI
13.9
-
16.2
73
Nov. 20, 2007(Stn. 1910, 0300 hrs)
PI
2.6
-
13.5
79
Nov. 21, 2007(Stn. 437, 1630 hrs)
SCFF
5.1
-
16.2
82
Nov. 25, 2007(Stn. 1812, 2100 hrs)
DFF
3.6
-
16.6
86Figure5 shows backscattering coefficients for co- (HH and VV) and crosspolarization (HV) configurations of each ice type. The two frost flower cases (DFF and SCFF) are plotted separately to exemplify differences in backscattering behavior on the basis of their different frost flower concentrations. The DFF and SCFF have a visually measured concentration of approximately >95% and 20%, respectively. While SCFYI is visually separable using HH, HV, and VV polarizations at all incidence angles (low backscatter), PI and DFYI signatures overlap and are difficult to separate from each other. This may be indicative of PI geometry within the scatterometer footprint, as PI comprises a series of upturned edges and flat areas of ice (see Figure 2). The curvature of upturned PI edges causes a backscatter response similar to that caused by the deformations (upturned ice) in the DFYI. DFF and SCFF are differentiable at HV and VV polarization at mid to far incidence angles.Co- (HH and VV) and cross- (HV) polarization backscatter intensities of snow-covered first-year sea ice (SCFYI), deformed first-year sea ice (DFYI), consolidated pancake ice (PI), snow-covered frost flowers (SCFF), and dense frost flowers (DFF).
(a)
(b)
(c)
(d)
(e)
(f)Mean coherences and polarization ratios for each ice type as a function of incidence angle grouping are documented in Table3. All sea ice types show high ρVVHH, indicating low depolarisation and primarily single (surface) backscattering. The ρHHVH for DFF is notably higher than that from the other ice types, which points to strong depolarisation caused by the frost flower structures. As shown in the previous section, a low value of ρRRLL indicates a rougher surface. At mid to far ranges in Table 3, the ρRRLL for DFYI is the lowest while for SCFYI it is the highest, which is consistent with the roughest and smoothest ice types, respectively. Furthermore, for frost flower-covered surfaces, that is, SCFF and DFF, the lower magnitude of ρRRLL is consistent with the higher concentration of frost flowers. At near-incidence angle range, the SCFYI shows higher roughness (i.e., lower ρRRLL=0.47, Table 3) compared to that of PI (0.61). This may be due to the fact that the snow is dry and has low salinity, which allows EM waves to penetrate through the snow. This is likely to provide roughness of snow-ice interface rather than air-snow interface. At mid incidence angle range, as expected, SCFYI shows lower roughness (i.e., higher ρRRLL=0.65, Table 3) compared to that of PI (0.61). Mid incidence angles are well suited for differentiating ice roughness/types using ρRRLL.Table 3
Mean C-band polarimetric coherences and ratios of selected ice types, for near (N), middle (M), and far (F) range incidence angle groupings (also shown graphically in Figure6). The number of data samples is: SCFYI, N=14; DFYI, N=8; PI, N=14; SCFF, N=18; DFF, N=10.
ρ
VVHH
ρ
HHVH
ρ
RRLL
γ
co (dB)
γ
cross (dB)
SCFYI
N
0.95
0.09
0.47
0.40
-
16.44
M
0.81
0.09
0.65
3.64
-
13.72
F
0.55
0.08
0.81
4.54
-
12.25
DFYI
N
0.88
0.10
0.42
1.36
-
12.38
M
0.91
0.08
0.52
2.47
-
12.94
F
0.58
0.05
0.58
1.86
-
11.40
PI
N
0.96
0.04
0.61
1.78
-
16.89
M
0.84
0.04
0.61
1.89
-
13.73
F
0.80
0.05
0.59
2.10
-
12.70
SCFF
N
0.96
0.07
0.73
0.89
-
17.88
M
0.89
0.07
0.70
-
0.59
-
15.83
F
0.66
0.05
0.77
0.46
-
14.08
DFF
N
0.84
0.15
0.35
2.41
-
09.71
M
0.91
0.15
0.59
3.65
-
12.06
F
0.79
0.12
0.72
3.23
-
12.23Looking at polarization ratios in Table3, the γco increases rapidly with incidence angle and is the highest at the far range for SCFYI. The γco ratio is also high for DFF, but it remains fairly constant across all incidence angles. High γco is also representative of saline ice surface (FYI in this case) or surface scattering. The presence of dry snow (~1-2 cm) allows the EM waves to penetrate through snow, which causes reflection from the ice-snow interface. The γco behavior of SCFYI is consistent with that of a surface, which is very smooth (i.e., a Bragg surface), where the ratio between backscattered H and V is only dependent on incidence angle and dielectric constant [28]. On the other hand, the γco behavior for DFF is consistent with that of a rough surface exhibiting backscatter from features with preferential vertical orientation [22]. Including the γcross ratio in this comparison further supports the distinction in backscattering mechanisms. The near range γcross ratio is much smaller for SCFYI than DFF, indicating it to be much smoother. The DFYI and DFF show the highest overall γcross, due to multiple scattering within deformities for DFYI and depolarisation caused by frost flowers for DFF.Figure6 shows box plots of coherences and polarization ratios of each ice type. Table 4 provides the significance values resulting from statistical tests for independence between each ice type based on a given coherence or ratio. All data in Figure 6 and Table 4 are based on the aforementioned incidence angle groupings from near to far range and, together, facilitate a conceptual approach to assessing the utility of each parameter for distinguishing ice types within an MIZ. Summarizing Figure 6 and Table 4, the near range ρHHVH and γcross provide the greatest separation between classes, while the far range ρVVHH and γco provide the greatest separation. By combining γco (far) with either of ρHHVH or γcross (near), all ice types are independent of each other. From Figure 4 it is known that a lower ρRRLL is associated with a rougher ice surface.Table 4
Matrix of significance values from non-parametric Kruskal-Wallis tests for independence between ice types based on polarimetric parameters and near (N), middle (M), and far (F) range groupings. The number of data samples is: SCFYI,N=14; DFYI, N=8; PI, N=14; SCFF, N=18; DFF, N=10.
Near range
Middle range
Far range
SCFYI
DFYI
PI
SCFF
DFF
SCFYI
DFYI
PI
SCFF
DFF
SCFYI
DFYI
PI
SCFF
DFF
ρ
VVHH
SCFYI
DFYI
.
285
.
015
.
094
PI
.
028
.
094
.
509
.
000
.
000
.
001
SCFF
.
463
.
119
.
002
.
011
.
322
.
001
.
014
.
199
.
000
DFF
.
004
.
013
.
028
.
000
.
005
.
483
.
000
.000
.
000
.
004
.
420
.000
ρ
HHVH
SCFYI
DFYI
.
201
.
201
.
035
PI
.
006
.
000
.
005
.
005
.
018
.
221
SCFF
.
429
.
008
.
000
.
125
.
156
.
001
.
046
.
184
.
277
DFF
.
002
.
002
.
000
.
000
.
001
.
000
.
000
.
000
.
001
.
000
.
000
.000
ρ
RRLL
SCFYI
DFYI
.
308
.
048
.
030
PI
.
164
.
106
.
082
.
513
.
000
.
357
SCFF
.
006
.
001
.
039
.
418
.
005
.
002
.
015
.
023
.
000
DFF
.
023
.
197
.
000
.
000
.
069
.
483
.
170
.
003
.
000
.
042
.
117
.002
γ
co
SCFYI
DFYI
.
285
.
018
.
001
PI
.
028
.
094
.
001
.
183
.
000
.
308
SCFF
.
463
.
119
.
002
.
000
.
000
.
000
.
000
.
001
.
000
DFF
.
004
.
013
.
028
.
000
.
466
.
017
.
001
.
000
.
000
.
003
.
005
.000
γ
cross
SCFYI
DFYI
.
005
.
149
.
041
PI
.
437
.
001
.
214
.
357
.
164
.
041
SCFF
.
118
.
000
.
179
.
006
.
005
.
001
.
000
.
013
.
001
DFF
.
000
.
010
.
000
.
000
.
004
.
204
.
006
.
000
.
420
.
042
.
354
.
00
0
#
#Bold numbers indicate important significant values.Figure 6
Box plots of coherences and polarization ratios of ice types based on near, middle, and far range incidence angle groupings. Significance values are provided in Table4.It is demonstrated using theory that lower values ofρRRLL indicate a rougher ice surface. Referring to Figures 4 and 6, ρRRLL is high for increasing incidence angles and for low surface roughness. This is only true for SCFYI and SCFF. In the presence of dry and fresh snow the volume contribution from FYI can be ignored, in which case ρRRLL dictates surface roughness of the snow-ice interface rather than air-snow interface. The coherence estimates are negligibly affected by the signal-to-noise ratio (typically >10 dB) during the processing of scatterometer data. These coherences can also be computed using polarimetric observations from space-based platforms.


## 3.1. Sea Ice Type Discrimination (Coherences and Ratios)

The date and hour of scatterometer data acquisitions corresponding to each sea ice type, as well as coincident meteorological parameters, namely, wind speed, air temperature, and relative humidity, are provided in Table2. The photographs of the selected ice samples are shown in Figure 2. With the exception of wind speed, there is negligible variation in meteorological conditions between ice type scans. As such, it is expected that between-scan, temperature-induced effects on the dielectric properties, and backscattering intensities from the different ice types are negligible.Table 2
Meteorological parameters associated with each ice type on different dates.
Sea ice type
Wind speed (m/s)
Air temperature (°C)
Relative humidity %
Nov. 15, 2007(Stn. 1117, 1400 hrs)
SCFYI
14.4
-
16.0
85
Nov. 19, 2007(Stn. 1100, 0030 hrs)
DFYI
13.9
-
16.2
73
Nov. 20, 2007(Stn. 1910, 0300 hrs)
PI
2.6
-
13.5
79
Nov. 21, 2007(Stn. 437, 1630 hrs)
SCFF
5.1
-
16.2
82
Nov. 25, 2007(Stn. 1812, 2100 hrs)
DFF
3.6
-
16.6
86Figure5 shows backscattering coefficients for co- (HH and VV) and crosspolarization (HV) configurations of each ice type. The two frost flower cases (DFF and SCFF) are plotted separately to exemplify differences in backscattering behavior on the basis of their different frost flower concentrations. The DFF and SCFF have a visually measured concentration of approximately >95% and 20%, respectively. While SCFYI is visually separable using HH, HV, and VV polarizations at all incidence angles (low backscatter), PI and DFYI signatures overlap and are difficult to separate from each other. This may be indicative of PI geometry within the scatterometer footprint, as PI comprises a series of upturned edges and flat areas of ice (see Figure 2). The curvature of upturned PI edges causes a backscatter response similar to that caused by the deformations (upturned ice) in the DFYI. DFF and SCFF are differentiable at HV and VV polarization at mid to far incidence angles.Co- (HH and VV) and cross- (HV) polarization backscatter intensities of snow-covered first-year sea ice (SCFYI), deformed first-year sea ice (DFYI), consolidated pancake ice (PI), snow-covered frost flowers (SCFF), and dense frost flowers (DFF).
(a)
(b)
(c)
(d)
(e)
(f)Mean coherences and polarization ratios for each ice type as a function of incidence angle grouping are documented in Table3. All sea ice types show high ρVVHH, indicating low depolarisation and primarily single (surface) backscattering. The ρHHVH for DFF is notably higher than that from the other ice types, which points to strong depolarisation caused by the frost flower structures. As shown in the previous section, a low value of ρRRLL indicates a rougher surface. At mid to far ranges in Table 3, the ρRRLL for DFYI is the lowest while for SCFYI it is the highest, which is consistent with the roughest and smoothest ice types, respectively. Furthermore, for frost flower-covered surfaces, that is, SCFF and DFF, the lower magnitude of ρRRLL is consistent with the higher concentration of frost flowers. At near-incidence angle range, the SCFYI shows higher roughness (i.e., lower ρRRLL=0.47, Table 3) compared to that of PI (0.61). This may be due to the fact that the snow is dry and has low salinity, which allows EM waves to penetrate through the snow. This is likely to provide roughness of snow-ice interface rather than air-snow interface. At mid incidence angle range, as expected, SCFYI shows lower roughness (i.e., higher ρRRLL=0.65, Table 3) compared to that of PI (0.61). Mid incidence angles are well suited for differentiating ice roughness/types using ρRRLL.Table 3
Mean C-band polarimetric coherences and ratios of selected ice types, for near (N), middle (M), and far (F) range incidence angle groupings (also shown graphically in Figure6). The number of data samples is: SCFYI, N=14; DFYI, N=8; PI, N=14; SCFF, N=18; DFF, N=10.
ρ
VVHH
ρ
HHVH
ρ
RRLL
γ
co (dB)
γ
cross (dB)
SCFYI
N
0.95
0.09
0.47
0.40
-
16.44
M
0.81
0.09
0.65
3.64
-
13.72
F
0.55
0.08
0.81
4.54
-
12.25
DFYI
N
0.88
0.10
0.42
1.36
-
12.38
M
0.91
0.08
0.52
2.47
-
12.94
F
0.58
0.05
0.58
1.86
-
11.40
PI
N
0.96
0.04
0.61
1.78
-
16.89
M
0.84
0.04
0.61
1.89
-
13.73
F
0.80
0.05
0.59
2.10
-
12.70
SCFF
N
0.96
0.07
0.73
0.89
-
17.88
M
0.89
0.07
0.70
-
0.59
-
15.83
F
0.66
0.05
0.77
0.46
-
14.08
DFF
N
0.84
0.15
0.35
2.41
-
09.71
M
0.91
0.15
0.59
3.65
-
12.06
F
0.79
0.12
0.72
3.23
-
12.23Looking at polarization ratios in Table3, the γco increases rapidly with incidence angle and is the highest at the far range for SCFYI. The γco ratio is also high for DFF, but it remains fairly constant across all incidence angles. High γco is also representative of saline ice surface (FYI in this case) or surface scattering. The presence of dry snow (~1-2 cm) allows the EM waves to penetrate through snow, which causes reflection from the ice-snow interface. The γco behavior of SCFYI is consistent with that of a surface, which is very smooth (i.e., a Bragg surface), where the ratio between backscattered H and V is only dependent on incidence angle and dielectric constant [28]. On the other hand, the γco behavior for DFF is consistent with that of a rough surface exhibiting backscatter from features with preferential vertical orientation [22]. Including the γcross ratio in this comparison further supports the distinction in backscattering mechanisms. The near range γcross ratio is much smaller for SCFYI than DFF, indicating it to be much smoother. The DFYI and DFF show the highest overall γcross, due to multiple scattering within deformities for DFYI and depolarisation caused by frost flowers for DFF.Figure6 shows box plots of coherences and polarization ratios of each ice type. Table 4 provides the significance values resulting from statistical tests for independence between each ice type based on a given coherence or ratio. All data in Figure 6 and Table 4 are based on the aforementioned incidence angle groupings from near to far range and, together, facilitate a conceptual approach to assessing the utility of each parameter for distinguishing ice types within an MIZ. Summarizing Figure 6 and Table 4, the near range ρHHVH and γcross provide the greatest separation between classes, while the far range ρVVHH and γco provide the greatest separation. By combining γco (far) with either of ρHHVH or γcross (near), all ice types are independent of each other. From Figure 4 it is known that a lower ρRRLL is associated with a rougher ice surface.Table 4
Matrix of significance values from non-parametric Kruskal-Wallis tests for independence between ice types based on polarimetric parameters and near (N), middle (M), and far (F) range groupings. The number of data samples is: SCFYI,N=14; DFYI, N=8; PI, N=14; SCFF, N=18; DFF, N=10.
Near range
Middle range
Far range
SCFYI
DFYI
PI
SCFF
DFF
SCFYI
DFYI
PI
SCFF
DFF
SCFYI
DFYI
PI
SCFF
DFF
ρ
VVHH
SCFYI
DFYI
.
285
.
015
.
094
PI
.
028
.
094
.
509
.
000
.
000
.
001
SCFF
.
463
.
119
.
002
.
011
.
322
.
001
.
014
.
199
.
000
DFF
.
004
.
013
.
028
.
000
.
005
.
483
.
000
.000
.
000
.
004
.
420
.000
ρ
HHVH
SCFYI
DFYI
.
201
.
201
.
035
PI
.
006
.
000
.
005
.
005
.
018
.
221
SCFF
.
429
.
008
.
000
.
125
.
156
.
001
.
046
.
184
.
277
DFF
.
002
.
002
.
000
.
000
.
001
.
000
.
000
.
000
.
001
.
000
.
000
.000
ρ
RRLL
SCFYI
DFYI
.
308
.
048
.
030
PI
.
164
.
106
.
082
.
513
.
000
.
357
SCFF
.
006
.
001
.
039
.
418
.
005
.
002
.
015
.
023
.
000
DFF
.
023
.
197
.
000
.
000
.
069
.
483
.
170
.
003
.
000
.
042
.
117
.002
γ
co
SCFYI
DFYI
.
285
.
018
.
001
PI
.
028
.
094
.
001
.
183
.
000
.
308
SCFF
.
463
.
119
.
002
.
000
.
000
.
000
.
000
.
001
.
000
DFF
.
004
.
013
.
028
.
000
.
466
.
017
.
001
.
000
.
000
.
003
.
005
.000
γ
cross
SCFYI
DFYI
.
005
.
149
.
041
PI
.
437
.
001
.
214
.
357
.
164
.
041
SCFF
.
118
.
000
.
179
.
006
.
005
.
001
.
000
.
013
.
001
DFF
.
000
.
010
.
000
.
000
.
004
.
204
.
006
.
000
.
420
.
042
.
354
.
00
0
#
#Bold numbers indicate important significant values.Figure 6
Box plots of coherences and polarization ratios of ice types based on near, middle, and far range incidence angle groupings. Significance values are provided in Table4.It is demonstrated using theory that lower values ofρRRLL indicate a rougher ice surface. Referring to Figures 4 and 6, ρRRLL is high for increasing incidence angles and for low surface roughness. This is only true for SCFYI and SCFF. In the presence of dry and fresh snow the volume contribution from FYI can be ignored, in which case ρRRLL dictates surface roughness of the snow-ice interface rather than air-snow interface. The coherence estimates are negligibly affected by the signal-to-noise ratio (typically >10 dB) during the processing of scatterometer data. These coherences can also be computed using polarimetric observations from space-based platforms.


## 4. Summary and Conclusions

The one-dimensional backscatter model of Schuler et al. [13] was modified to two dimensions of surface roughness by considering deviation in the orientation angles (i.e., the slopes) in azimuth and range direction simultaneously as an improvement in the model. Parameters derived from the fully polarimetric C-band microwave backscatter response from sea ice targets were demonstrated to have utility for small-scale (cm level) sea ice roughness identification. Circular coherence has been investigated for its usefulness in discriminating surface roughness among other polarimetric parameters. Circular coherence is theoretically shown to detect measurement sensitivity to surface roughness.The conclusions with reference to objective 1 are as follows. It was shown theoretically that theρRRLL decreases exponentially with increasing surface roughness. However, ρRRLL responds to both roughness (standard deviation) and dielectric constant (scattering matrix) of the surface in the case when the orientations of the ice target in azimuth direction are changed. It remains challenging to separate roughness effects from the dielectric effects using C-band backscatter measurements. ρRRLL independently does not provide a robust sea ice roughness discrimination scheme. However, ρRRLL provides an improved insight of sea ice surface roughness combined with other polarimetric coherences and channel ratios in the chosen samples. The experimental data also show that rougher ice surface exhibits lower mean value of ρRRLL (Table 3, Figure 4), though a complete validation of the effect of changing orientations of ice floe on ρRRLL is required. This would require polarimetric backscattering data and surface roughness information to be acquired at different lines of sight (i.e., orientation of ice floes). Unfortunately, difficulties associated with extreme weather conditions and limitations to navigation in the Arctic restrict such detailed data acquisition; however, a tank experiment could be a useful alternative.The utility of C-band polarimetric coherences and ratios is addressed in the light of objective 2 as follows: for coherences,ρVVHH is smaller at far range incidence angles for all ice types. ρHHVH is less sensitive to roughness and is not a good discriminator of roughness. Regarding channel ratios, based on Kruskal-Wallis test, γco is more sensitive to increasing surface roughness compared to γcross and demonstrates utility for separating ice types compared to the other observed parameters.The knowledge obtained through surface-based polarimetric coherences and ratios can readily be extended to discriminate sea ice roughness on small scales using C-band microwave satellites (currently in orbit RADARSAT-2, RISAT-1). Future work will be to develop an algorithm combining all polarimetric coherences and ratios to discriminate individual ice type in a MIZ. These observations may become particularly useful for satellite measurements once planned SAR constellations (Sentinel series) systems are available, as currently planned by National Aeronautics and Space Administration and European Space Agency.



---

*Source: 567182-2013-05-22.xml*