Next Article in Journal
Improving Jason-2 Sea Surface Heights within 10 km Offshore by Retracking Decontaminated Waveforms
Previous Article in Journal
Agricultural Soil Spectral Response and Properties Assessment: Effects of Measurement Protocol and Data Mining Technique
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Data Synergy between Altimetry and L-Band Passive Microwave Remote Sensing for the Retrieval of Sea Ice Parameters—A Theoretical Study of Methodology

1
Ministry of Education Key Laboratory for Earth System Modeling, Department of Earth System Science, Tsinghua University, Beijing 100084, China
2
Department of Atmospheric and Environmental Sciences, University at Albany, State University of New York, Albany, NY 12222, USA
3
State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics (LASG), Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2017, 9(10), 1079; https://doi.org/10.3390/rs9101079
Submission received: 6 August 2017 / Revised: 18 October 2017 / Accepted: 20 October 2017 / Published: 23 October 2017
(This article belongs to the Section Ocean Remote Sensing)

Abstract

:
Accurate knowledge of the sea ice parameters, including the thickness and the snow depth over sea ice, are key to both climate change studies and operational forecast in polar regions. The estimation of these parameters mainly relies on satellite based remote sensing, and current retrieval algorithms usually focus on the retrieval of a single parameter under simple assumptions over the other. In this article, we explore the potential of combined retrieval of both sea ice thickness and snow depth through the data synergy two types of concurrent observations of the sea ice cover: the active altimetry and the L-band passive remote sensing. The data synergy is based on two physical constrains: (1) L-band (1.4 GHz) radiation model for the sea ice cover, and (2) the hydrostatic equilibrium as used in satellite altimetry. Two schemes of data synergy are proposed: (1) the synergy between L-band brightness temperature ( T B ) from passive microwave remote sensing and sea ice freeboard ( F B i c e ) as measured by radar altimetry, and (2) the synergy between L-band T B and snow freeboard ( F B s n o w ) as measured by laser altimetry. Based on retrievability studies, we show that both parameters can be retrieved using the two sets of data. Specifically, we show that there is potential problem of ill-posedness for the synergy between L-band T B and F B s n o w , with two possible retrieval solutions for a small portion of the solution space. On the other hand, the synergy between L-band T B and F B i c e is always well-posed. In terms of sensitivity, lower uncertainty is witnessed for thin ice for the retrieval with F B i c e , while the retrieval with F B s n o w shows advantage for thick ice. Besides the input parameters of T B , F B i c e and F B s n o w , the uncertainty associated with certain model parameters such as snow and ice densities is not negligible for the uncertainty estimation of the retrieved parameters. Verification is carried out with observational data from Operation IceBridge (OIB) campaigns and SMOS satellite, showing that both sea ice thickness and snow depth can be attained by the proposed retrieval algorithms. These algorithms serve as the basis for large-scale retrieval with satellite remote sensing data, including concurrent observation of the Arctic Ocean by independent satellite campaigns such as SMOS, CryoSat-2 and ICESat.

Graphical Abstract

1. Introduction

The sea ice forms as the result of the cooling and the phase change of the sea surface [1]. As an important factor in the climate system, the sea ice influences the overall radiative budget of the earth through albedo effects. Besides, it actively modulates the exchanges of heat, momentum and gas between the polar atmosphere and ocean, and plays an important role in the propagation of salinity and freshwater [2,3,4,5]. The snow cover is an integrated component of the sea ice, which is the direct indicator of the polar atmospheric hydrological cycle. With higher albedo and better thermal insulation as compared with the sea ice, the snow cover is also of crucial importance to the associated atmosphere-ocean interactions. Therefore, the accurate estimation of the sea ice thickness and snow depth over sea ice are the key to the understanding of processes in polar regions and the projections of future changes [6], especially under the context of drastic shrinkage and thinning of the Arctic sea ice cover during the recent decades. Aside from scientific questions, it is also indispensable to applications such as the sea ice forecast operations and seasonal prediction [7].
Large-scale observation of sea ice thickness is challenging due to the scarcity of in-situ measurements. Active satellite remote sensing based on radar or laser altimetry is the major source of sea ice thickness estimations, which are available during the freeze-up period. By using freeboard measurements from satellite altimetry such as ICESat [8] or CryoSat-2 [9], the thickness of sea ice is retrieved under certain assumptions over the snow depth and snow loading [10,11,12,13]. However, large uncertainty is usually present for satellite altimetry when the sea ice thickness is small (less than 1 m) [10,11]. Besides active altimetry, sea ice thickness can also be retrieved through passive microwave satellite remote sensing. The thickness of thin sea ice can be estimated with thermal imagery [14,15], which might be susceptible to cloud contamination. Besides, the thin ice thickness (10 to 20 cm) can be also estimated with the brightness temperature ( T B ) measurements from various frequency bands of passive microwave radiometer such as Special Sensor Microwave Imager (SSM/I) and Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) sensors [16,17,18,19]. Recently, in Kaleschke et al. [20] and related works, the retrieval of thin ice thickness with L-band (1.4 GHz) T B measurements from Soil Moisture and Ocean Salinity (SMOS) satellite is carried out, based on a forward sea ice radiation model. It is shown that SMOS-derived ice thickness has a lower uncertainty as compared with active altimetry for thin ice, and can serve as a complementary data set [20,21]. However, for thicker ice (over 0.5 m), the uncertainty is very large due to the saturation of T B with respect to the thickness [22]. A few recent work [21,23] focused on the synergy of the retrieved sea ice thickness from SMOS and that based on CryoSat-2. Ricker et al. [24] created a merged weekly product based on optimal interpolation scheme. It is worth noting that current retrieval of sea ice thickness are all based on simple assumptions over the snow loading over the sea ice. For example, the retrieval based on CryoSat-2 (or the merged product) is based on (adapted) climatology snow depth in [25], and that based on ICESat resorts to the combination of: (1) the snowfall from reanalysis data throughout the freeze-up period, and (2) the climatological snow depth at the beginning of the freeze-up period [11]. On one hand, the validity of climatological snow depth may be compromised due to its limited efficacy [26] and large uncertainty that arise from interannual variability [13,25]. On the other hand, the snow depth as estimated from reanalysis usually has unconstrained uncertainty and is subjected to model biases.
As compared with sea ice thickness, the observations of the snow depth over the sea ice are also sparse [27,28]. The climatology of snow depth in the Arctic Ocean has been developed in Warren et al. [25] by compiling multi-decadal (1954 to 1991) in-situ measurements with drifting stations. Snow depth retrieval with satellite based passive microwave remote sensing are proposed in various studies [29,30]. Comiso et al. [29] obtained snow depth over Antarctic sea ice and first-year sea ice (FYI) in the Arctic based on difference between 19 and 37 GHz T B (vertical polarization) from AMSR-E. Maaß et al. [31] (and related works) adopted a radiation model that incorporates the radiative properties of both snow and sea ice, and achieved snow depth retrieval over thick Arctic multi-year ice (MYI) using SMOS T B . It is worth noting that the sea ice thickness is a prerequisite for the snow depth retrieval, and the retrievable depth is about 35 cm due to T B saturation.
In this study, we carry out theoretical study over the simultaneous retrieval of both sea ice thickness and snow depth, by data synergy of concurrent active and passive remote sensing of the sea ice cover. Two data synergy schemes are proposed: (1) the synergy of L-band T B and radar altimetry such as CryoSat-2 (for sea ice freeboard measurements), and (2) the synergy between L-band T B and laser altimetry such as ICESat and ICESat-2 (for snow or total freeboard measurements). Similar to current active and passive remote sensing algorithms [9,11,31], the proposed retrieval framework relies on physical constraints for the retrieval, including the L-band radiation model and the hydrostatic equilibrium relationship. However, the existing algorithms usually target at a single sea ice parameter by using a single observational dataset, while the proposed framework aims at the simultaneous retrieval of both parameters with two independent observational datasets. Theoretical studies, including the analysis of the well-posedness of the retrieval problem, show that there exists good retrievability for the synergy of two independent datasets. The overall validity of the proposed algorithms is further verified with successful retrieval for typical scenarios, as well as the synergy study based on SMOS and airborne data of OIB. Furthermore, through sensitivity studies we analyze the quantitative contribution of various parameters to the uncertainty of the retrieved parameters. The proposed data synergy serves as a theoretical basis for the basin-scale retrieval of sea ice parameters using concurrent active and passive satellite remote sensing data, including SMOS and CryoSat-2 or ICESat-2.

2. Retrievability Studies

The active and passive remote sensing of the sea ice cover can be summarized in Figure 1. Radar and laser altimetry detects the elevation of the main reflectance plane of the surface by active scanning. Either the sea ice freeboard or the snow freeboard is attained, which is converted to sea ice thickness following buoyancy relationships. On the other hand, passive remote sensing (in microwave bands) such as SMOS measures the radiative properties of the earth’s surface, which usually features lower spatial resolution but better spatial coverage than active altimetry. Active and passive remote sensing are usually carried out by independent campaigns, such as ICESat and CryoSat-2 for altimetry, and SMOS for L-band passive remote sensing.
For the synergy of L-band T B and freeboard measurements, we adopt two existing physical models that constrain the observational data with the sea ice parameters. The first model is the L-band radiation model [32], which simulates the L-band (1.4 GHz) brightness temperature of the sea ice cover based on its parameters. The second model is the hydrostatic equilibrium model based on buoyancy relationship, which is widely used in satellite altimetry. Section 2.1 gives a short introduction to these models, including the parameters as adopted by the retrieval algorithms. Based on these two forward models, the inverse problem of retrieval is studied in Section 2.2 for the synergy between T B and F B i c e and that between T B and F B s n o w . The analysis of the properties of the solution space and the potential problem of ill-posedness is the basis of the design of retrieval methods in Section 3.

2.1. Forward Models

The first physical model as adopted by the retrieval is the L-band (1.4 GHz) radiation model. The original model was developed for modeling the radiative transfer of soil moisture at L-band, as in Burke et al. [33]. The adapted version of the model was applied for the sea ice that are related to the retrieval with SMOS data [31]. Specifically, in Kaleschke et al. [21] introduced the detailed modeling of the dielectric properties of the sea ice, and in Maaß et al. [31] the model was extended to include a snow layer over the sea ice and further applied for the snow depth retrieval. In Zhou et al. [32], the model was further improved to include more realistic characterization of the small-scale sea ice variabilities. Due to the sensitivity of the radiative properties to the salinity in the sea ice [21], the vertical salinity profile is integrated in the multi-layer formulation of the model. Sea ice type dependent salinity profile was adopted to reflect the more thorough salinity drainage and flushing of MYI as compared with FYI. Besides, as discussed in Zhou et al. [32], the open water or (refrozen) sea ice leads could have profound impact on the overall L-band T B , we limit the discussion in this article to the theoretical studies which: (1) only apply to normal Arctic winter conditions, and (2) involve no mixture of open water or sea ice types. In Zhou et al. [32], the improved multi-layer radiation model was verified with observational data including SMOS and Operation IceBridge (OIB [34]). There was good agreement between the modeled and the observational T B ( R 2 as high as 0.81, with a correction factor of about 1.8 K).
Figure 2 shows the simulated T B with respect to sea ice parameters ( h i and h s ), under normal Arctic winter conditions (the surface temperature of −30 C). There is non-linear relationship between T B and h i or h s . With the deepening of either the sea ice or the snow cover (i.e., increasing h i or h s ), the value of T B gradually saturates for both FYI and MYI. For h i larger than about 2.5 m (1.2 m) for MYI (FYI), the value of T B saturates. For the retrieval of h i based on a prescribed h s , this implies that a small perturbation/error in the observed T B would result in a large difference in the retrieved h i , as in Tian-Kunze et al. [22]. Also due to the nonlinear relationship, the retrieval cannot be formulated in an explicit form. We use Equation (1) to represent the nonlinear relationship between T B and related parameters, including: surface temperature ( T s f c ), sea ice type ( i c e _ t y p e ), sea water temperature ( T w a t e r ) and sea water salinity ( S w a t e r ) at the bottom of the sea ice.
T B = TB ( h i , h s , T s f c , i c e _ t y p e , T w a t e r , S w a t e r )
The second physical model as adopted by the retrieval is the hydrostatic equilibrium model based on the buoyancy relationship between h i , h s and the density of sea ice, snow and sea water ( ρ i c e , ρ s n o w and ρ w a t e r ). This model is the basis of both radar and laser satellite altimetry [8,10,35]. For radar altimetry as in CryoSat-2 [13,36,37], the sea ice thickness can be estimated with F B i c e as derived from satellite data, according to Equation (2). For laser altimetry as in ICESat [8,38], the same model is applied with a different form involving F B s n o w , as shown in Equation (3).
h i = ρ w a t e r ρ w a t e r ρ i c e · F B i c e + ρ s n o w ρ w a t e r ρ i c e · h s
h i = ρ w a t e r ρ w a t e r ρ i c e · F B s n o w ρ w a t e r ρ s n o w ρ w a t e r ρ i c e · h s
Since the linear relationship between h i and F B i c e or F B s n o w , the retrieval of h i can be formulated in the explicit form above. The value of h s is treated as a priori knowledge which can be estimated from other sources (see Section 1).
For the retrieval, we treat T B , F B i c e and F B s n o w as input parameters, and h i and h s the retrieved parameters. Other parameters as used by the models are categorized as model parameters, which are either treated as constants or assumed to be derived from other sources. For the radiation model, T s f c is assumed to be derived from remote sensing or reanalysis data. In this article, we mainly consider the typical winter condition of the Arctic, and the surface temperature is assumed to be −30 C Besides, we also consider the relatively warmer condition ( T s f c = −15 C) in the sensitivity studies, which is becoming more common with global warming [39]. The values of T w a t e r and S w a t e r are considered as constants ( T w a t e r = −1.8 C and S w a t e r = 33 g/kg). Following common practice in satellite altimetry, density values are assumed constant: ρ w a t e r = 1024 kg/m 3 and ρ i c e = 915 kg/m 3 , as derived from field measurements [40], and ρ s n o w = 320 kg/m 3 following Warren et al. [25]. Although these model parameters are assumed to be constant, the uncertainty as caused by these parameters are all analyzed systematically for the sensitivity studies in Section 4.
With the transformed forms of Equations (2) and (3), we can also derive F B i c e and F B s n o w based on h i and h s . In Figure 2, the colored lines (with labels) represent constant F B i c e (blue lines) and constant F B s n o w (red lines). Since F B i c e and F B s n o w are input parameters, the solution to the retrieval based on T B and freeboard synergy resides on the corresponding constant freeboard lines. For the synergy between T B and F B i c e , T B increases monotonically on each constant F B i c e line, in spite of the saturation and minor decrease of T B with respect to h i . This implies good potential of retrieving h i and h s based on T B and F B i c e . However, for the synergy between T B and F B s n o w , under certain values of F B s n o w (e.g., 0.12 m for FYI and 0.30 m for MYI), there is no monotonic increase of T B with the increase of h i . This implies that for certain combinations of T B s n o w and T B , there may exist more than 1 solution to the retrieval problem. It is worth noting that for thick sea ice, there is potentially better sensitivity for both schemes of synergy, as indicated by the non-flat constant freeboard lines with respect to T B .

2.2. Retrievability Analysis

We construct the solution space for the retrieval problems by scanning the sea ice parameters. For each scanned values of h i and h s , we generate the values of T B , F B i c e and F B s n o w . We exclude the cases involving inundation (i.e., F B i c e < 0 ), which is uncommon for Arctic regions. The solution space for each of the data synergy scheme is plotted and examined with respect to the input parameters.
For the synergy between T B and F B i c e , Figure 3a,b show the solution space for h i and h s respectively for FYI, and Figure 3c,d show those for MYI. A reasonable range of the solution space for both h i and h s is included ( h i within about 10 m and h s within 1 m). With a certain value of F B i c e , there exists a range of possible T B values. For FYI, with the increase in F B i c e , the lowest possible value of T B also increases. For MYI, there is a slight decrease of the lowest bound of T B when F B i c e is high. Due to the saturation of T B with respect to either h i or h s , there also exists an upper bound for T B . Under certain value of T B (or F B i c e ), the solutions of both h i and h s become larger with the increase of F B i c e (or T B ).
Figure 4 shows the solution space for the synergy between T B and F B s n o w . As compared with the corresponding information in Figure 3, there exists evident differences. The most prominent difference from the synergy between T B and F B i c e is the presence of two potential solutions for certain part of the solution space. Therefore, the retrieval problem of the synergy between T B and F B s n o w is not well-posed. In each subfigure, the portion of the solution space that involves two potential solutions is outlined by dotted lines in the main part. The other solution of this portion of the solution space is also shown by embedded graph in each subfigure. The corresponding input parameters are as follows: (1) for FYI, F B s n o w around 0.1 m and T B around 237 K, (2) for MYI, F B s n o w around 0.37 m and T B around 247 K. For FYI, this portion consists of a small part (3.32%) of the solution space, and the corresponding two solutions for h i ( h s ) are about 0.4 m and 0.9 m (0.1 m and 0.07 m). For MYI, the proportion is higher (6.67%), and the two potential solutions for h i ( h s ) are about 1.3 m and 2.8 m (0.3 m and 0.12 m).
We also consider the relatively warmer Arctic condition with T s f c = −15 C, instead of the normal winter condition ( T s f c = −30 C). Examination of the corresponding solution space is carried out in a similar manner. The results (not shown) indicate that under the relatively warmer condition, the range of T B is from 247 K to 255 K (250 K to 257 K) when F B i c e is large in FYI (MYI) algorithm. This range is narrower as compared to winter condition (Figure 3). However, the well-posedness (ill-posedness) of the retrieval problem for the synergy between T B and F B i c e ( F B s n o w ) is not changed. Besides, the portion of the solution space that involves two solutions takes up about 9.28 % ( 5.49 % ) for FYI (MYI).
To summarize, there exists retrievability of sea ice parameters with the observational data of L-band T B and freeboard measurements. For typical parameters, both h i and h s can be attained given the combination of T B and F B i c e or F B s n o w . The synergy between T B and F B i c e is always well-posed, while for that between T B and F B s n o w there is potential of two solutions. These characteristics of the solution space should be accounted for by the retrieval algorithms.

3. Retrieval Algorithms

To accommodate the different characteristics of the data synergy schemes, we construct two retrieval algorithms with the forward physical models. For the synergy between T B and F B i c e , due to the retrieval problem is well-posed, we start with a minimum guess for snow depth ( h s = 0 ), and the corresponding h i can be computed according to Equation (2). At each iteration, the values of h i and h s are used to compute the TB using the forward radiation model. The modeled TB is compared against the observed TB. If the modeled TB is lower, h s is then increased for the next iteration. Otherwise, a small step in increasing h s is adopted. Therefore an iterative process (involving binary search) is carried out to locate the proper value of h s (as well as h i ) that generates the observed TB. The outline of the iterative algorithm is shown in Figure 5.
For the synergy between T B and F B s n o w , since there exists possibility of two potential solutions, instead of a simple iterative process in Figure 5, two steps are involved for the retrieval. Firstly, a scanning process within the possible range of h s is carried out. The scanning targets at individual values of h s by a sufficiently small step (e.g., 1 cm). Potential solutions are detected when the values of modeled T B from two adjacent values of h s reside on different side of the observed T B . For the second step, for each potential solution, a binary search process can then be carried out. For the case involving two potential solutions, both of them are reported by the retrieval algorithm. The outline of this algorithm is shown in Figure 6.
In order to verify the retrieval algorithms, we construct 8 idealized retrieval scenarios that are typical for the Arctic region. As shown in Table 1, these scenarios cover both sea ice types, as well as a range of sea ice thickness and snow depth. The derived parameters based on physical models in Section 2.1 are also shown, including F B i c e , F B s n o w and L-band T B . We denote these scenarios by number I to VIII, and their locations in the solution space for both types of data synergy are indicated in Figure 3 and Figure 4.
We carry out retrieval with the derived values of F B i c e , F B s n o w and L-band T B , using the corresponding retrieval algorithm. Table 2 shows the retrieval results. By comparison with the true values of the parameters, we demonstrate that the true values of h i and h s are attained with both types of data synergy. For the synergy between T B and F B i c e , the retrieval solution is unique for all scenarios. For the synergy between T B and F B s n o w , two potential solutions are possible for scenario II and VI. The true solution is included in the retrieval result, and the other solution is shown in red for each scenario. It is further verified this solution does produce the same F B s n o w and T B as the true solution. As compared with the true solution, the snow is thicker and the sea ice is thinner. Consistent with the retrievability study, it is indicated that scenario II and VI are included in the portion of the solution space that contain ambiguous solutions, as in Figure 4.

4. Sensitivity Studies

In this section we carry out sensitivity study over the retrieved parameters and the effect of both the input and the model parameters. We first examine that of input parameters in Section 4.1 and Section 4.2, and further examine that of model parameters in Section 4.3. The uncertainty for these parameters are listed in Table 3. The values are the in standard deviation ( σ ) as derived from typical observations or datasets. For example, the uncertainty of T B is 0.5 K, as specified by SMOS data [22], and the uncertainty of S w a t e r represents Arctic basin-wide variability of surface salinity. Therefore, these values serve as references to the uncertainty of the data during actual retrieval practice. Uncertainty in the retrieved parameters are quantified by either differentiation or Monte Carlo simulations with perturbations to both individual parameters and a set of parameters. In Monte-Carlo simulations, all the perturbations to all parameters except the freeboard values follow normal distribution with the prescribed standard deviations. The perturbation to freeboard values ( F B i c e and F B s n o w ) follow log-normal distributions to avoid negative freeboards. Results show that T B and freeboard plays different roles in the uncertainty of h i and that of h s . Furthermore, model parameters play a non-negligible role as compared with input parameters. For further study with realistic data, both types of parameters should be accounted for the uncertainty estimation of the retrieved parameters.

4.1. Effect of Input Parameters for the Synergy between T B and F B i c e

Based on the solution spaces of both types of data synergy, we compute the differences of the retrieved h i and h s as caused by perturbations to either T B or freeboard values. First, we compute the absolute uncertainty of h i and h s , denoted ϵ a b s o l u t e : ϵ a b s o l u t e | y = Δ y Δ x · σ x , where x is the input parameter ( T B or F B i c e or F B s n o w ), and y the retrieved parameter ( h i or h s ). The values of Δ x and Δ y are derived by differentiation within the solution space in Figure 3 and Figure 4, and the value of Δ y Δ x characterizes the rate of change in y as caused by a small perturbation in x. Therefore ϵ a b s o l u t e represents the uncertainty of y in terms of standard deviation. Second, we compute the relative uncertainty based on the absolute uncertainty and the actual value of the parameter. The relative uncertainty is defined as: ϵ r e l a t i v e | y = ϵ a b s o l u t e | y y × 100 % , which is the percentage of the uncertainty of y with respect to y.
For the data synergy between T B and F B i c e , Figure 7 and Figure 8 show the uncertainty analysis for perturbations to F B i c e and T B , respectively. As shown in Figure 7, the uncertainty in h i as caused by that of F B i c e is generally lower than 10% (subfigure a and c). For comparison, as in Figure 8, the uncertainty in h i as caused by T B is generally lower than that caused by F B i c e when T B is relatively low. However, when T B is high ( h i > 4 m and h s > 0.9 m for FYI, or h i > 7 m and h s > 70 cm for MYI), the uncertainty as caused by T B is larger than that of F B i c e . This change is attributed to the saturation feature of T B for thick ice. It is worth noting that these parameters are relatively rare during Arctic observations. For h s , when comparing Figure 7 and Figure 8, it is shown that different from the case of h i , T B always plays a dominant role in the uncertainty of h s . The uncertainty in h s as caused by F B i c e is only higher than 10% when T B is lower than about 238 K or higher than 250 K for FYI, and lower than 246 K or higher than 253% for MYI. However, T B consistently causes nearly 10% or higher uncertainty in h s for both FYI and MYI. As a result, T B causes larger uncertainty in h s for a larger portion of the solution space than F B i c e . Similar to h i , when both h i and h s are relatively large, there exists increase in the uncertainty of h s for both absolute and relative terms (>0.4 m or 20% for FYI, and >0.5 m or 30% for MYI).
We further carry out sensitivity studies for the 8 typical retrieval scenarios as introduced in Section 3. In order to evaluate the uncertainty in the retrieved h i and h s , we apply Monte Carlo simulations with the retrieval algorithms. For each type of data synergy, three simulations are carried out, containing random perturbations to the value of freeboard, to that of T B , and to both of them. Each simulation contains 5000 samples. Each simulation sample contains a perturbation (or a set of perturbations), and the retrieval is carried out with the perturbed parameter(s). Perturbations follow normal distribution for T B and log-normal distribution for freeboard (see Table 3). The uncertainty in the retrieved h i and h s is then estimated by computing the standard deviation among the samples.
Table 4 (first 3 columns of the results) shows the relative uncertainty as estimated from the simulations for the synergy between T B and F B i c e . The relative uncertainty of parameter y, denoted ϵ r e l a t i v e | y is computed as σ y y × 100 % , where σ y is the standard deviation as derived from Monte-Carlo samples and y the true value of the parameter. In order to contrast the relative importance between T B and F B i c e , the larger value is shown in bold. For both h i and h s , F B i c e plays an more dominant role or relative thin ice (scenario I, II and V), and T B plays an more important role for thick MYI (scenario VII and VIII). For other scenarios (III, IV, VI), F B i c e plays a dominant role for h i , but T B caused higher uncertainty for h s . Also the uncertainty in h i and h s as caused by both input parameters is lower than 11% and 28%, respectively. Besides, they are generally lower than the combined values as caused by each parameter under the assumption of independence. This implies that the uncertainty in h i and h s is correlated, although the uncertainty in input parameters are independent.
Aside from the normal winter condition of the Arctic, we also consider a relatively warm condition ( T s f c =−15 C ), to reflect the effect of global warming for the typical scenarios. The uncertainty is estimated with Monte-Carlo simulations, and the results are shown in Table 5. Comparison with Table 4 shows that there is little change or slight increase in the uncertainty caused by F B i c e , but the increase in that caused by T B is much larger. In relative terms, T B has become the dominant source of uncertainty of h i for scenario IV and VI, and that of h s for scenario II. With the warmer condition, the saturation of T B with respect to sea ice parameters is more premature, which is the reason for the increase of its role in the uncertainty of the retrieved parameters.

4.2. Effect of Input Parameters for the Synergy between T B and F B s n o w

We carry out sensitivity study for the data synergy between T B and F B s n o w as in Section 4.1. On the solution space, the most prominent difference from the synergy between T B and F B i c e is the region with two potential solutions. In Figure 9 and Figure 10, the sensitivity of the retrieved parameters to F B s n o w and T B in the solution space are shown, respectively. The sensitivity for the portion that contains two potential solutions are shown by embedded figures, similar to Figure 4. For scenario II and IV with two potentials solutions, the effect of double solution is accounted for during the uncertainty estimation. Both F B s n o w and T B have direct influence over the retrieved h i and h s . When the sea ice is thin, the relative uncertainty in h i as caused by F B s n o w grows to over 15% for both FYI and MYI. When the sea ice thickness is high, this value drops under 2%. Similar situation applies to h s , for shallower snow cover or thinner sea ice, the relative uncertainty as caused by F B s n o w is high. For comparison, the effects of T B is relatively more dominant that of F B s n o w than the data synergy of F B i c e and T B . The influence on the uncertainty of h s is more prominent. For FYI, when T B is low, the relative uncertainty of h s rises up to 20%, and for MYI, to over 15%. Table 6 (first 3 columns) lists the uncertainty estimation for representative scenarios, which is in direct comparison with Table 4. The relative uncertainty as caused by both input parameters drops with the increase in sea ice thickness and snow depth. As contrast from Table 4, T B causes a larger portion of the uncertainty in both h i and h s (6 for h i and 7 for h s , among 8 scenarios), while F B s n o w and T B cause uncertainty in h i and h s on equal terms for thin FYI and thin MYI.
When compared with the data synergy between T B and F B i c e , the major difference is that for relatively thinner sea ice, the synergy between T B and F B i c e yields lower uncertainty in sea ice parameters, and for thicker ice, that between T B and F B s n o w shows advantage. We also carry out the study for these scenarios under relatively warm Arctic conditions (shown in Table 7). For relatively thinner sea ice cover, there is no general trend of the uncertainty in the retrieved parameters. However, for scenarios with thicker ice, the relative uncertainty as caused by T B generally double or triple. Similar to the other type of synergy, the uncertainty caused by freeboard measurements stays unchanged, due to the irrelevance of the physical model of buoyancy to the change of temperature from −30 C to −15 C.

4.3. Uncertainty Analysis from Model Parameters

Apart from input parameters, the effect of model parameters on the retrieved sea ice parameters are evaluated in terms of uncertainty. We extend the Monte-Carlo simulations in Section 4.1 and Section 4.2 to include perturbations to model parameters. The analyzed parameters are listed in Table 3. Since the effects of T w a t e r and S w a t e r are very small, causing less than 0.01% change in the retrieved h i and h s . This is due to two reasons: (1) they have very small impact on the density of the water and the corresponding retrieval process of altimetry, and (2) the radiative properties of the sea ice covered ocean are dominated by that of the sea ice cover, and the small change in the radiative property of the sea water as caused by these parameters has little impact on the sea ice retrieval. Therefore, we focus on the other 3 parameters: ρ i c e , ρ s n o w and T s f c . The density values mainly affect the physical model for altimetry, influencing the quantitative relationship between h i and h s under a linear relationship (see Equations (2) and (3)). Due to the observed large uncertainty and seasonal evolution of ρ s n o w [25] as well as their important role in altimetry [41,42], their uncertainty values are chosen accordingly to reflect their change during the Arctic winter season, ranging from 220 kg/m 3 to 420 kg/m 3 (see also [43,44]). The uncertainty of ρ i c e is set to 10 kg/m 3 , after [41,43,44]. On the other hand, T s f c mainly affects the retrieval results through its influence over the vertical temperature structure of the sea ice cover and its radiative properties [32]. The uncertainty of T s f c of about 1 K is derived from OIB dataset [45], which reflects the variability of surface temperature within a range of 25 km which is commonly used for the production of polar remote sensing data products.
In Table 4, Table 5, Table 6 and Table 7, the uncertainty as caused by ρ i c e , ρ s n o w , T s f c , and their combination are shown in column 4 to 7 (of the results). When compared with input parameters, they show non-negligible influence over the uncertainty of retrieved parameters. Their combined effect is comparable to that of input parameters across all scenarios. Besides, for the synergy between T B and F B i c e , with the increase of h i and h s , the uncertainty caused by these model parameters increase, but for the synergy between T B and F B s n o w , the decrease of the uncertainty for h i and h s is witnessed. This is in general consistency with the analysis for input parameters, and coincides with the relative advantage of the two types of synergy over different range of sea ice parameter. Under warmer conditions, the difference as caused by density values generally remain unchanged, while that of T s f c increase by 100% or higher for both h i and h s .
For the uncertainty estimation of the typical scenarios, we also produce a set of Monte-Carlo simulation that contain perturbations to all the parameters (including both input and model parameters), as shown by the last column in each table. For the retrieval problem, the model parameters are mainly derived from other supportive data sources, and not directly provided by satellite altimetry or L-band remote sensing data. However, according to the analyses, they can have equally profound influence over the retrieved parameters, which is consistent with existing studies such as the analysis of satellite altimetry [41,42]. Their effect should also be accounted for during the actual retrieval problems for the data synergy proposed by this article.

5. Verification with OIB and SMOS Data

5.1. Data

For the verification, we use airborne measurements from OIB and SMOS satellite based TB data, by considering the co-registered scanning of the sea ice by the two independent data sets. OIB data are provided as 40 m resolution along-track scanning of the sea ice cover. The onboard snow radar and laser altimeter measure the nadir snow depth and total freeboard, and the sea ice thickness can be derived with these parameters. All the data from OIB campaigns during year 2012 and 2015 are used for the retrieval and verifications [34]. The level-4 SMOS TB data product [46] is adopted, which are regridded on the 12.5 km EASE grid.
SMOS TB field on the same date of each OIB campaign is attained for the data synergy. Due to the inherent resolution difference between OIB (limited swath) and SMOS, we treat OIB as samples of the underlying sea ice cover that corresponds to a single TB value in the daily SMOS TB field. Although the resolution of the SMOS TB product is 12.5 km, the native resolution of SMOS measurements is about 40 km. Therefore, all the measurements from OIB that are within the 9 EASE grid cells are treated as samples, covering an area of 37.5 km × 37.5 km which is consistent in general with the native resolution of SMOS measurements.
For the verification, we align the spatial resolution of OIB data to that of SMOS. Points with mixed sea ice type, open water or lead are excluded for further retrieval and verification, in order to exclude the uncertainty associated with the mixture of the sea ice cover. For each SMOS TB, the arithmetic mean values of F B s n o w , h s , F B i c e and h i among all the OIB samples are computed. These values are denoted F B s n o w ¯ , h s ¯ , F B i c e ¯ and h i ¯ , respectively. For the synergy between radar altimetry and L-band TB, we treat OIB measurements of F B i c e as that measured by radar altimetry. Specifically we use F B i c e ¯ and T B for the retrieval, and compare the retrieve value of h i and h s with h i ¯ and h s ¯ , respectively. For the synergy between laser altimetry and L-band TB, we use F B s n o w ¯ and T B for the retrieval, and the verification is also between the retrieved parameters and observed ones, including h i ¯ and h s ¯ .

5.2. Results and Discussion

The TB as modeled using h i ¯ and h s ¯ for available OIB points is first compared with the corresponding SMOS TB. The root-mean-square error (RMSE) of the TB is about 3.1 K. Since OIB campaigns contain certain areas with recursive fly-over and better spatial coverage, we further compute the RMSE of TB for these points (95 percentile in terms of OIB sample count). The RMSE drops from 3.1 K to 1.41 K, which indicates that the coverage of airborne (or any type of altimetry) may potentially play an important role in the data synergy.
Figure 11 shows the result of verification with all available data. For both types of data synergy, there exists good match between the retrieval result and the observation for both sea ice thickness ( h i ) and snow depth ( h s ). For data synergy between F B i c e and T B , the R 2 for h i is 0.87 (0.86) for the least squares fitting (under the slope = 1 constraint) between retrieved parameter and observation, and no significant bias is present for h i (Figure 11a). For h s , the R 2 is 0.26, which is lower than that for h i (Figure 11b). The least squares fit under the constraint of slope = 1 yields an underestimation of about 4 cm for h s .
In order to analyze the source of error in the retrieved parameters, we further limit the retrieval to the points with lower difference between the modeled TB and the observed TB (within 1.5 K), shown in Figure 12. Based on all the points within 1.5 K of the observed TB (accounting for half of all points), the R 2 for h i and h s for the synergy between F B i c e and T B increase to 0.90 and 0.49, respectively. Besides, two extra improvements are present: (1) a minor decrease in the bias of h s is also witnessed (about −3 cm), and (2) the least squares fitting is very close to that under the constraint of slope = 1. As a comparison, in Figure 11b, the under-estimation of h s for relatively thick snow cover is present, lowering the quality of matching of the retrieved h s to observation. Limiting the retrieval to points with lower difference in T B removes these points in effect. Also the quality in the retrieved h i is also improved (0.87 to 0.90 for R 2 ). This result indicates that the error in T B is a major source of the error in the retrieved parameters, especially h s .
For the data synergy between F B s n o w and T B , the R 2 for h i and h s are 0.78 and 0.45, respectively (Figure 11c,d). Also, there is very small difference between the least squares fit and that under the constraint of slope = 1. If the retrieval is limited to points within 1.5 K, there exists increase in R 2 for h s (0.45 to 0.50), while that for h i remains the same (Figure 12c,d). As discussed in Section 3, the retrieval with F B s n o w and T B is not always well-posed, as contrast from that with F B i c e and T B . Among all the data in the verification, there exists potential of two solutions for 13.6% of all points (mainly MYI). When two potential solutions are encountered during retrieval, the one with smaller h s is adopted by default which is used for the verification.
For the actual retrieval with satellite altimetry data, several factors should be considered due to the difference from the airborne data. First, each satellite altimetry campaign features specific resolution settings, with fine spatial scans for laser altimetry (about 70 m for ICESat), and relatively large footprint for radar altimetry (e.g., 300 m by 1.5 km for CryoSat-2). Since the spatial coverage of OIB plays an important role in the error for both TB and retrieval, the specific information for each satellite campaign should be taken into account for the synergy with TB measurements (e.g., 40 km for SMOS). Better spatial coverage of satellite altimetric scans is also present for cross-over points of the scanning tracks on the daily basis. Second, there exists spatial variability of the sea ice cover which can be represented by the altimetry scans due to their relatively higher spatial resolution. Since there exists nonlinear relationship of TB to sea ice parameters as indicated by the radiation model, the variability on the small scale for both sea ice thickness and snow depth should be accounted for the actual synergy of remote sensing data. Third, other factors such as the presence of open water and (refrozen) leads can play a profound role in modulating the large-scale TB, as studied in Zhou et al. [32]. In this study, points involving these factors are ignored for the verification. For the actual retrieval with satellite data, third-party data sources such as sea ice lead maps, can be integrated in the retrieval process. Their role in the uncertainty of the retrieved parameters can be accounted for in a systematic way, by adopting the approach of Monte-Carlo simulations as in Section 4.

6. Summary and Future Work

In this study we propose the data synergy between active altimetry and L-band passive microwave remote sensing of the sea ice cover, for the retrieval of both sea ice thickness and snow depth. The synergy is based on two physical models: (1) the hydrostatic equilibrium relationship which is widely adopted in satellite altimetry, and (2) the L-band radiation model [32]. Specifically, the L-band radiation model is an improved version based on Maaß et al. [30] and verified with observational datasets of OIB and SMOS. The Arctic winter condition is assumed, during which the satellite altimetry is carried out [8,36]. Theoretical study over the retrievability and the well-posedness of the retrieval problem are carried out, and it is demonstrated that both sea ice parameters can be retrieved with two observational data. Furthermore, during the data synergy between L-band T B and F B s n o w , we discover the potential of the existence of two solutions, due to the non-monotonous relationship between sea ice thickness and T B . Although it consists of a small portion of the solution space for retrieval (3.32% and 6.67% for FYI and MYI, respectively), this feature should be accounted for in the retrieval algorithm of the data synergy between T B and F B s n o w . On the contrary, the data synergy between T B and F B i c e is always well-posed. Based on solution space analysis and retrieval with typical scenarios, we further carry out sensitivity studies over the retrieved parameters. While both input parameters affect the uncertainty of the retrieved parameters, the two types of data synergy have advantage over each other over different range of sea ice parameters. For thin sea ice, the synergy between T B and F B i c e shows lower uncertainty, while for thick sea ice, that between T B and F B s n o w produces lower uncertainty. Besides, under warmer conditions in the Arctic during winter, the uncertainty of the retrieved parameters increases, mainly due to the saturation of L-band T B . Furthermore, it is demonstrated that model parameters such as density and temperature have comparable influence over the retrieved parameters as input parameters. This indicates that for the actual retrieval, the uncertainty of both types of parameters should be accounted for in a systematic way, such as Monte-Carlo simulation as introduced in Section 4.
The proposed data synergy and corresponding retrieval algorithms serve as a basis for the retrieval with satellite measurements of the Arctic sea ice. Traditional retrieval methods usually utilize a single type of observational data, such as sea ice (or snow) freeboard in altimetry or L-band T B for the studies involving SMOS, and the objective for the retrieval is a single type of sea ice parameters ( h i or h s ). As a comparison, the proposed data synergy schemes aim at: (1) the utilization of two independent but concurrent observations, and (2) the retrieval of both h i and h s . By using SMOS TB and the co-registering OIB data, we simulate the retrieval based on both types of altimetry and the passive L-band measurement, and demonstrate that on the spatial scale of SMOS TB, both h i and h s can be retrieved by combining active and passive remote sensing data. For the study with actual data such as CryoSat-2 and SMOS, extension to the theoretical framework should be made. As contrast from the idealized retrieval scenario in this study, the difference in the spatial and temporal coverage of different observations should be considered for the design of retrieval algorithms. Specifically, SMOS satellite typically produces daily Arctic wide observations, but satellite altimetry as in CryoSat-2 can only achieve basin coverage on a monthly basis. Besides, there exists inherent resolution difference among various data (e.g., 40 km for SMOS, about 1 km for CryoSat-2, 40 m for OIB). This resolution difference should be considered during retrieval, using specific spatial coverage information of both data during the synergy. The inherent spatial variability of the sea ice cover, such as the sea ice cover with relatively lower ice concentration or the presence of leads, should also be accounted for by treating it as a mixture of various types, similar to [32]. With respect to the uncertainty for the large-scale retrieval problem, the Monte-Carlo simulation based methods in Section 4 can be adopted, with a systematic treatment to both observational data and model parameters. These research directions serve as future directions for the retrieval of sea ice parameters on the basin-scale.

Acknowledgments

The authors would like to thank the editors and potential reviewers in their effort to help improve the manuscript. This work is partially support by the National Science Foundation of China under the grant number of 41575076 and National Key R&D Program of China under the grant number of 2017YFA0603902.

Author Contributions

S.X. conceived and designed the algorithms and experiments; L.Z. and S.X. performed the experiments; S.X., L.Z. and J.L. analyzed the data; L.Z., S.X., J.L., B.W. and H.L. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Thomas, D.N.; Dieckmann, G.S. Sea Ice, 2nd ed.; Wiley-Blackwell: Oxford, UK, 2010. [Google Scholar]
  2. McPhee, M.; Proshutinsky, A.; Morison, J.H.; Steele, M.; Alkire, M. Rapid change in freshwater content of the Arctic Ocean. Geophys. Res. Lett. 2009, 36. [Google Scholar] [CrossRef]
  3. Kurtz, N.; Markus, T.; Farrell, S.; Worthen, D.; Boisvert, L. Observations of recent Arctic sea ice volume loss and its impact on ocean-atmosphere energy exchange and ice production. J. Geophys. Res. Oceans 2011, 116. [Google Scholar] [CrossRef]
  4. Rabe, B.; Karcher, M.; Kauker, F.; Schauer, U.; Toole, J.M.; Krishfield, R.A.; Pisarev, S.; Kikuchi, T.; Su, J. Arctic Ocean basin liquid freshwater storage trend 1992–2012. Geophys. Res. Lett. 2014, 41, 961–968. [Google Scholar] [CrossRef]
  5. Blazey, B.; Holland, M.; Hunke, E. Arctic Ocean sea ice snow depth evaluation and bias sensitivity in CCSM. Cryosphere 2013, 7, 1887–1900. [Google Scholar] [CrossRef]
  6. Stroeve, J.; Barrett, A.; Serreze, M.; Schweiger, A. Using records from submarine, aircraft and satellites to evaluate climate model simulations of Arctic sea ice thickness. Cryosphere 2014, 8, 1839–1854. [Google Scholar] [CrossRef] [Green Version]
  7. Lindsay, R.; Zhang, J. Assimilation of ice concentration in an ice-ocean model. J. Atmos. Ocean. Technol. 2006, 23, 742–749. [Google Scholar] [CrossRef]
  8. Kwok, R.; Cunningham, G.; Zwally, H.; Yi, D. Ice, Cloud, and land Elevation Satellite (ICESat) over Arctic sea ice: Retrieval of freeboard. J. Geophys. Res. Oceans 2007, 112. [Google Scholar] [CrossRef] [Green Version]
  9. Laxon, S.W.; Giles, K.A.; Ridout, A.L.; Wingham, D.J.; Willatt, R.; Cullen, R.; Kwok, R.; Schweiger, A.; Zhang, J.; Haas, C.; et al. CryoSat-2 estimates of Arctic sea ice thickness and volume. Geophys. Res. Lett. 2013, 40, 732–737. [Google Scholar] [CrossRef]
  10. Laxon, S.; Peacock, N.; Smith, D. High interannual variability of sea ice thickness in the Arctic region. Nature 2003, 425, 947–950. [Google Scholar] [CrossRef] [PubMed]
  11. Kwok, R.; Cunningham, G. ICESat over Arctic sea ice: Estimation of snow depth and ice thickness. J. Geophys. Res. Oceans 2008, 113. [Google Scholar] [CrossRef]
  12. Ricker, R.; Hendricks, S.; Helm, V.; Skourup, H.; Davidson, M. Sensitivity of CryoSat-2 Arctic sea-ice freeboard and thickness on radar-waveform interpretation. Cryosphere 2014, 8, 1607–1622. [Google Scholar] [CrossRef] [Green Version]
  13. Tilling, R.L.; Ridout, A.; Shepherd, A. Near-real-time Arctic sea ice thickness and volume from CryoSat-2. Cryosphere 2016, 10, 2003–2012. [Google Scholar] [CrossRef]
  14. Yu, Y.; Rothrock, D.A. Thin ice thickness from satellite thermal imagery. J. Geophys. Res. Atmos. 1996, 1012, 25753–25766. [Google Scholar] [CrossRef]
  15. Mäkynen, M.; Cheng, B.; Similä, M. On the accuracy of thin-ice thickness retrieval using MODIS thermal imagery over Arctic first-year ice. Ann. Glaciol. 2013, 54, 87–96. [Google Scholar] [CrossRef]
  16. Martin, S.; Drucker, R.; Kwok, R.; Holt, B. Improvements in the estimates of ice thickness and production in the Chukchi Sea polynyas derived from AMSR-E. Geophys. Res. Lett. 2005, 32. [Google Scholar] [CrossRef]
  17. Tamura, T.; Ohshima, K.I.; Markus, T.; Cavalieri, D.J.; Nihashi, S.; Hirasawa, N. Estimation of thin ice thickness and detection of fast ice from SSM/I data in the Antarctic Ocean. J. Atmos. Ocean. Technol. 2007, 24, 1757–1772. [Google Scholar] [CrossRef]
  18. Nihashi, S.; Ohshima, K.I.; Tamura, T.; Fukamachi, Y.; Saitoh, S.I. Thickness and production of sea ice in the Okhotsk Sea coastal polynyas from AMSR-E. J. Geophys. Res. Atmos. 2009, 114, 157–165. [Google Scholar] [CrossRef]
  19. Singh, R.K.; Oza, S.R.; Vyas, N.K.; Sarkar, A. Estimation of thin ice thickness from the advanced microwave scanning radiometer-EOS for coastal polynyas in the Chukchi and Beaufort Seas. IEEE Trans. Geosci. Remote Sens. 2011, 49, 2993–2998. [Google Scholar] [CrossRef]
  20. Kaleschke, L.; Tian-Kunze, X.; Maaß, N.; Mäkynen, M.; Drusch, M. Sea ice thickness retrieval from SMOS brightness temperatures during the Arctic freeze-up period. Geophys. Res. Lett. 2012, 39. [Google Scholar] [CrossRef]
  21. Kaleschke, L.; Maaß, N.; Haas, C.; Hendricks, S.; Heygster, G.; Tonboe, R. A sea-ice thickness retrieval model for 1.4 GHz radiometry and application to airborne measurements over low salinity sea-ice. Cryosphere 2010, 4, 583–592. [Google Scholar] [CrossRef] [Green Version]
  22. Tian-Kunze, X.; Kaleschke, L.; Maaß, N.; Mäkynen, M.; Serra, N.; Drusch, M.; Krumpen, T. SMOS-derived thin sea ice thickness: Algorithm baseline, product specifications and initial verification. Cryosphere 2014, 8, 997–1018. [Google Scholar] [CrossRef] [Green Version]
  23. Kaleschke, L.; Tian-Kunze, X.; Maaß, N.; Ricker, R.; Hendricks, S.; Drusch, M. Improved retrieval of sea ice thickness from SMOS and CryoSat-2. In Proceedings of the 2015 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Milan, Italy, 26–31 July 2015; pp. 5232–5235. [Google Scholar]
  24. Ricker, R.; Hendricks, S.; Kaleschke, L.; Tian-Kunze, X.; King, J.; Haas, C. A weekly Arctic sea-ice thickness data record from merged CryoSat-2 and SMOS satellite data. Cryosphere 2017, 11, 1607–1623. [Google Scholar] [CrossRef]
  25. Warren, S.G.; Rigor, I.G.; Untersteiner, N.; Radionov, V.F.; Bryazgin, N.N.; Aleksandrov, Y.I.; Colony, R. Snow depth on Arctic sea ice. J. Clim. 1999, 12, 1814–1829. [Google Scholar] [CrossRef]
  26. Webster, M.A.; Rigor, I.G.; Nghiem, S.V.; Kurtz, N.T.; Farrell, S.L.; Perovich, D.K.; Sturm, M. Interdecadal changes in snow depth on Arctic sea ice. J. Geophys. Res. Oceans 2014, 119, 5395–5406. [Google Scholar] [CrossRef]
  27. Markus, T.; Cavalieri, D.J. Snow Depth Distribution Over Sea Ice in the Southern Ocean from Satellite Passive Microwave Data. In Antarctic Sea Ice: Physical Processes, Interactions and Variability; American Geophysical Union: Washington, DC, USA, 2013; pp. 19–39. [Google Scholar]
  28. Gerland, S.; Haas, C. Snow-depth observations by adventurers traveling on Arctic sea ice. Ann. Glaciol. 2011, 52, 369–376. [Google Scholar] [CrossRef]
  29. Comiso, J.C.; Cavalieri, D.J.; Markus, T. Sea ice concentration, ice temperature, and snow depth using AMSR-E data. IEEE Trans. Geosci. Remote Sens. 2003, 41, 243–252. [Google Scholar] [CrossRef]
  30. Maaß, N.; Kaleschke, L.; Tian-Kunze, X.; Tonboe, R.T. Snow thickness retrieval from L-band brightness temperatures: A model comparison. Ann. Glaciol. 2015, 56, 9–17. [Google Scholar] [CrossRef]
  31. Maaß, N.; Kaleschke, L.; Tian-Kunze, X.; Drusch, M. Snow thickness retrieval over thick Arctic sea ice using SMOS satellite data. Cryosphere 2013, 7, 1971–1989. [Google Scholar] [CrossRef]
  32. Zhou, L.; Xu, S.; Liu, J.; Lu, H.; Wang, B. Improving L-band radiation model and representation of small-scale variability to simulate brightness temperature of sea ice. Int. J. Remote Sens. 2017, 38, 7070–7084. [Google Scholar] [CrossRef]
  33. Burke, W.; Schmugge, T.; Paris, J. Comparison of 2.8-and 21-cm microwave radiometer observations over soils with emission model calculations. J. Geophys. Res. Oceans 1979, 84, 287–294. [Google Scholar] [CrossRef]
  34. Kurtz, N.; Farrell, S.; Studinger, M.; Galin, N.; Harbeck, J.; Lindsay, R.; Onana, V.; Panzer, B.; Sonntag, J. Sea ice thickness, freeboard, and snow depth products from Operation IceBridge airborne data. Cryosphere 2012, 6, 4771–4827. [Google Scholar] [CrossRef]
  35. Kwok, R. Annual cycles of multiyear sea ice coverage of the Arctic Ocean: 1999–2003. J. Geophys. Res. Oceans 2004, 109. [Google Scholar] [CrossRef]
  36. Hendricks, S.; Ricker, R.; Helm, V.; Haas, C.; Skourup, H.; Herber, A.; Schwegmann, S.; Gerdes, R.; Davidson, M. CryoSat-2 Sea-Ice Freeboard and Thickness. In Proceedings of the ESA Living Planet Symposium, Edinburgh, Scotland, 9–13 September 2013. [Google Scholar]
  37. Kwok, R.; Cunningham, G. Variability of Arctic sea ice thickness and volume from CryoSat-2. Philos. Trans. R. Soc. A 2015, 373. [Google Scholar] [CrossRef] [PubMed]
  38. Kwok, R. Satellite remote sensing of sea-ice thickness and kinematics: A review. J. Glaciol. 2010, 56, 1129–1140. [Google Scholar] [CrossRef]
  39. Ricker, R.; Hendricks, S.; Girard-Ardhuin, F.; Kaleschke, L.; Lique, C.; Tian-Kunze, X.; Nicolaus, M.; Krumpen, T. Satellite-observed drop of Arctic sea ice growth in winter 2015–2016. Geophys. Res. Lett. 2017, 44, 3236–3245. [Google Scholar] [CrossRef]
  40. Wadhams, P.; Tucker, W., III; Krabill, W.; Swift, R.; Comiso, J.; Davis, N. Relationship between sea ice freeboard and draft in the Arctic Basin, and implications for ice thickness monitoring. J. Geophys. Res. 1992, 97, 20325–20334. [Google Scholar] [CrossRef]
  41. Tilling, R.L.; Ridout, A.; Shepherd, A.; Wingham, D.J. Increased Arctic sea ice volume after anomalously low melting in 2013. Nat. Geosci. 2015, 8, 643–646. [Google Scholar] [CrossRef]
  42. Zygmuntowska, M.; Rampal, P.; Ivanova, N.; Smedsrud, L.H. Uncertainties in Arctic sea ice thickness and volume: New estimates and implications for trends. Cryosphere 2014, 8, 705–720. [Google Scholar] [CrossRef] [Green Version]
  43. Alexandrov, V.; Sandven, S.; Wahlin, J.; Johannessen, O. The relation between sea ice thickness and freeboard in the Arctic. Cryosphere 2010, 4, 373–380. [Google Scholar] [CrossRef] [Green Version]
  44. Kurtz, N.T.; Farrell, S.L.; Studinger, M.; Galin, N.; Harbeck, J.P.; Lindsay, R.; Onana, V.D.; Panzer, B.; Sonntag, J.G. Sea ice thickness, freeboard, and snow depth products from Operation IceBridge airborne data. Cryosphere 2013, 7, 1035–1056. [Google Scholar] [CrossRef]
  45. Kurtz, N.T.; Farrell, S.L. Large-scale surveys of snow depth on Arctic sea ice from Operation IceBridge. Geophys. Res. Lett. 2011, 38. [Google Scholar] [CrossRef]
  46. Tian-Kunze, X.; Kaleschke, L.; Maaß, N. SMOS Daily Polar Gridded Brightness Temperatures; ICDC, University of Hamburg, Digital Media: Hamburg, Germany, 2012. [Google Scholar]
Figure 1. Schematic diagram of active and passive remote sensing of the sea ice cover. Sea ice parameters ( h i and h s ), as well as the directly observed values of freeboard ( F B i c e amd F B s n o w ) and brightness temperature ( T B ) are also shown.
Figure 1. Schematic diagram of active and passive remote sensing of the sea ice cover. Sea ice parameters ( h i and h s ), as well as the directly observed values of freeboard ( F B i c e amd F B s n o w ) and brightness temperature ( T B ) are also shown.
Remotesensing 09 01079 g001
Figure 2. The relationship between T B and sea ice parameters ( h i and h s ) under typical Arctic winter conditions (surface temperature is −30 C) for FYI (a) and MYI (b). The corresponding constant freeboard lines are shown. Constant F B s n o w lines and corresponding values (in m) of F B s n o w are in red, and those for F B i c e in blue.
Figure 2. The relationship between T B and sea ice parameters ( h i and h s ) under typical Arctic winter conditions (surface temperature is −30 C) for FYI (a) and MYI (b). The corresponding constant freeboard lines are shown. Constant F B s n o w lines and corresponding values (in m) of F B s n o w are in red, and those for F B i c e in blue.
Remotesensing 09 01079 g002
Figure 3. The solution space of h i (a) and h s (b) based on synergy between T B and F B i c e for FYI. (c,d) show the solution spaces for MYI. The results correspond to typical Arctic winter conditions. Asterisks and Roman numerals (I to VIII) in the figures represent representative retrieval scenarios with typical Arctic sea ice parameters.
Figure 3. The solution space of h i (a) and h s (b) based on synergy between T B and F B i c e for FYI. (c,d) show the solution spaces for MYI. The results correspond to typical Arctic winter conditions. Asterisks and Roman numerals (I to VIII) in the figures represent representative retrieval scenarios with typical Arctic sea ice parameters.
Remotesensing 09 01079 g003
Figure 4. Same as Figure 3 but for the synergy between T B and F B s n o w . Subfigures are ordered as in Figure 3, with (a,b) showing the sol2zuole tion space of h i and h s for FYI respectively, and (c,d) showing those for MYI.
Figure 4. Same as Figure 3 but for the synergy between T B and F B s n o w . Subfigures are ordered as in Figure 3, with (a,b) showing the sol2zuole tion space of h i and h s for FYI respectively, and (c,d) showing those for MYI.
Remotesensing 09 01079 g004
Figure 5. Retrieval algorithm based on data synergy between L-band T B and F B i c e .
Figure 5. Retrieval algorithm based on data synergy between L-band T B and F B i c e .
Remotesensing 09 01079 g005
Figure 6. Retrieval algorithm based on data synergy between L-band T B and F B s n o w .
Figure 6. Retrieval algorithm based on data synergy between L-band T B and F B s n o w .
Remotesensing 09 01079 g006
Figure 7. Sensitivity study of retrieved parameters to freeboard measurements. Results for the synergy between T B and F B i c e are shown. (a,b) show the sensitivity of the retrieved h i and h s to F B i c e for FYI, respectively. (c,d) show the results for MYI. The filled color contour is the absolute uncertainty of h i or h s (in m e t e r s ), and the labeled lines are the relative uncertainty (units: %). The uncertainty for F B i c e is quantified in terms of standard deviations according to Table 3. The absolute uncertainty in h i or h s (denoted ϵ a b s o l u t e ) is the quantitative difference in h i or h s as caused by perturbation to F B i c e or T B by 1 standard deviation. The relative uncertainty of y is computed as: ϵ r e l a t i v e | y = ϵ a b s o l u t e | y y × 100 % , where y is the corresponding value of h i or h s .
Figure 7. Sensitivity study of retrieved parameters to freeboard measurements. Results for the synergy between T B and F B i c e are shown. (a,b) show the sensitivity of the retrieved h i and h s to F B i c e for FYI, respectively. (c,d) show the results for MYI. The filled color contour is the absolute uncertainty of h i or h s (in m e t e r s ), and the labeled lines are the relative uncertainty (units: %). The uncertainty for F B i c e is quantified in terms of standard deviations according to Table 3. The absolute uncertainty in h i or h s (denoted ϵ a b s o l u t e ) is the quantitative difference in h i or h s as caused by perturbation to F B i c e or T B by 1 standard deviation. The relative uncertainty of y is computed as: ϵ r e l a t i v e | y = ϵ a b s o l u t e | y y × 100 % , where y is the corresponding value of h i or h s .
Remotesensing 09 01079 g007
Figure 8. Same as Figure 7 but for the effect of the uncertainty in T B for the synergy between T B and F B i c e . The order of subfigures are also the same as Figure 7, with (a,b) showing the sensitivity of the retrieved h i and h s to T B for FYI, and (c,d) the results for MYI.
Figure 8. Same as Figure 7 but for the effect of the uncertainty in T B for the synergy between T B and F B i c e . The order of subfigures are also the same as Figure 7, with (a,b) showing the sensitivity of the retrieved h i and h s to T B for FYI, and (c,d) the results for MYI.
Remotesensing 09 01079 g008
Figure 9. Same as Figure 7 but for the effect of the uncertainty in F B s n o w for the synergy of T B and F B s n o w . The order of subfigures are also the same as Figure 7, with (a,b) showing the sensitivity of the retrieved h i and h s to F B s n o w for FYI, and (c,d) the results for MYI.
Figure 9. Same as Figure 7 but for the effect of the uncertainty in F B s n o w for the synergy of T B and F B s n o w . The order of subfigures are also the same as Figure 7, with (a,b) showing the sensitivity of the retrieved h i and h s to F B s n o w for FYI, and (c,d) the results for MYI.
Remotesensing 09 01079 g009
Figure 10. Same as Figure 7 but for the effect of the uncertainty in T B for the synergy of T B and F B s n o w . The order of subfigures are also the same as Figure 7, with (a,b) showing the sensitivity of the retrieved h i and h s to T B for FYI, and (c,d) the results for MYI.
Figure 10. Same as Figure 7 but for the effect of the uncertainty in T B for the synergy of T B and F B s n o w . The order of subfigures are also the same as Figure 7, with (a,b) showing the sensitivity of the retrieved h i and h s to T B for FYI, and (c,d) the results for MYI.
Remotesensing 09 01079 g010
Figure 11. Verification with co-registered OIB and SMOS TB data. (a,c) show the results of data synergy between F B i c e and T B for h i and h s , respectively. (b,d) show those between F B s n o w and T B . Blue triangles represent FYI, and blue squares MYI. Solid black line in each subfigure is the 1:1 line, the (dotted-)dashed line the least squares fit line (under the constraint that the slope is 1).
Figure 11. Verification with co-registered OIB and SMOS TB data. (a,c) show the results of data synergy between F B i c e and T B for h i and h s , respectively. (b,d) show those between F B s n o w and T B . Blue triangles represent FYI, and blue squares MYI. Solid black line in each subfigure is the 1:1 line, the (dotted-)dashed line the least squares fit line (under the constraint that the slope is 1).
Remotesensing 09 01079 g011
Figure 12. Similar to Figure 11, but limited to points of modeled T B within +/− 1.5 K of SMOS TB. The ordering of the subfigures are the same as in Figure 11, with (a,b) showing the results for h i and h s during the synergy between F B i c e and T B , and (c,d) the results during the synergy between F B s n o w and T B .
Figure 12. Similar to Figure 11, but limited to points of modeled T B within +/− 1.5 K of SMOS TB. The ordering of the subfigures are the same as in Figure 11, with (a,b) showing the results for h i and h s during the synergy between F B i c e and T B , and (c,d) the results during the synergy between F B s n o w and T B .
Remotesensing 09 01079 g012
Table 1. Sea ice parameters ( h i and h s ) and observed parameters ( F B i c e , F B s n o w and T B ) for 8 typical retrieval scenarios for the Arctic.
Table 1. Sea ice parameters ( h i and h s ) and observed parameters ( F B i c e , F B s n o w and T B ) for 8 typical retrieval scenarios for the Arctic.
Ice TypeScenariohi (m)hs (m)FBice (m)FBsnow (m)TB (K)
FYII0.50.050.03760.0876231.5
II1.00.030.09710.1271238.2
III1.50.100.12840.2284243.7
IV2.50.250.18800.4380246.5
MYIV1.50.150.11270.2627243.8
VI2.50.150.21920.3692248.9
VII3.00.350.21000.5600252.0
VIII5.00.400.40720.8072250.4
Table 2. Retrieved sea ice parameters ( h i and h s ) for typical scenarios in Table 1. For the case with two potential solutions, the bold texts show the true solution, while the red texts in parentheses the other potential solution. The case with double solutions is only possible for the synergy of T B and F B s n o w .
Table 2. Retrieved sea ice parameters ( h i and h s ) for typical scenarios in Table 1. For the case with two potential solutions, the bold texts show the true solution, while the red texts in parentheses the other potential solution. The case with double solutions is only possible for the synergy of T B and F B s n o w .
ResultsData SynergyScenario
IIIIIIIVVVIVIIVIII
Retrieved h i (m) T B & F B i c e 0.50001.00001.50002.50001.50002.50003.00005.0000
T B & F B s n o w 0.50001.00001.50002.50001.50002.50003.00005.0000
(0.5552)(1.5499)
True h i (m)0.501.001.502.501.502.503.005.00
Retrieved h i (m) T B & F B i c e 0.05000.03000.10000.25000.15000.15000.35000.4000
T B & F B s n o w 0.05000.03000.10000.25000.15000.15000.35000.4000
(0.0989)(0.2971)
True h s (m)0.050.030.100.250.150.150.350.40
Table 3. The uncertainty ( σ ) of input parameters and model parameters in sensitivity studies.
Table 3. The uncertainty ( σ ) of input parameters and model parameters in sensitivity studies.
Input Parameters
σ F B i c e (m) σ F B s n o w (m) σ T B (K)
0.010.010.5
Model Parameters
σ ρ i c e (kg/m 3 ) σ ρ s n o w (kg/m 3 ) σ T s f c (K) σ T w a t e r (K) σ S w a t e r (psu)
1033.31.0310.22
Table 4. Relative uncertainty (units: %) in the retrieved h i and h s for the synergy between T B and F B i c e , as estimated by Monte-Carlo simulations. Normal Arctic winter condition is assumed (surface temperature of −30 C). In each Monte-Carlo simulation sample, the retrieval process is carried out with a random perturbation to the designated parameter(s). The perturbations are Gaussian, with the standard deviation in Table 3. The relative uncertainty ϵ r e l a t i v e | y is given by: σ y y × 100 % , where σ y is the standard deviation of y estimated from samples of the corresponding Monte-Carlo simulation, and y the true value (y for h i or h s ). Sensitivity studies are carried out for each input or model parameter, each set of parameters, as well as all parameters (in different columns). The numbers in bold are the relatively larger one between the two values associated with the input parameters.
Table 4. Relative uncertainty (units: %) in the retrieved h i and h s for the synergy between T B and F B i c e , as estimated by Monte-Carlo simulations. Normal Arctic winter condition is assumed (surface temperature of −30 C). In each Monte-Carlo simulation sample, the retrieval process is carried out with a random perturbation to the designated parameter(s). The perturbations are Gaussian, with the standard deviation in Table 3. The relative uncertainty ϵ r e l a t i v e | y is given by: σ y y × 100 % , where σ y is the standard deviation of y estimated from samples of the corresponding Monte-Carlo simulation, and y the true value (y for h i or h s ). Sensitivity studies are carried out for each input or model parameter, each set of parameters, as well as all parameters (in different columns). The numbers in bold are the relatively larger one between the two values associated with the input parameters.
ScenarioResults FB ice TB Both Input Param. ρ ice ρ snow T sfc All Model Param.All Param.
I h i 11.18 % 1.15 % 10.69 % 3.93 % 0.72 % 1.52 % 4.39 % 11.50 %
h s 27.61 % 3.91 % 26.73 % 10.27 % 5.60 % 5.17 % 13.10 % 29.55 %
II h i 7.79 % 1.18 % 7.93 % 5.92 % 0.10 % 1.64 % 6.16 % 9.85 %
h s 18.45 % 13.42 % 23.24 % 13.95 % 8.28 % 18.61 % 25.16 % 33.80 %
III h i 6.74 % 2.39 % 7.08 % 7.71 % 0.53 % 2.54 % 8.24 % 11.24 %
h s 2.23 % 12.22 % 12.40 % 2.77 % 5.55 % 12.95 % 14.52 % 20.21 %
IV h i 4.75 % 3.92 % 6.27 % 9.29 % 1.31 % 3.54 % 10.04 % 11.73 %
h s 3.54 % 13.36 % 13.72 % 6.75 % 4.03 % 12.06 % 14.67 % 20.05 %
V h i 3.71 % 1.62 % 3.96 % 4.10 % 1.23 % 1.79 % 4.45 % 5.98 %
h s 9.67 % 5.41 % 10.55 % 10.69 % 5.12 % 5.99 % 12.93 % 16.49 %
VI h i 3.21 % 2.16 % 3.77 % 6.12 % 0.58 % 1.96 % 6.48 % 7.75 %
h s 3.85 % 12.28 % 12.52 % 6.38 % 5.49 % 11.10 % 13.77 % 18.95 %
VII h i 3.88 % 6.22 % 7.36 % 10.13 % 1.55 % 3.69 % 12.65 % 16.59 %
h s 2.66 % 18.18 % 18.54 % 8.54 % 3.96 % 10.78 % 18.97 % 33.30 %
VIII h i 3.21 % 6.40 % 7.13 % 12.83 % 1.35 % 4.38 % 15.32 % 18.03 %
h s 5.66 % 27.24 % 27.80 % 23.34 % 4.12 % 18.63 % 34.59 % 50.31 %
Table 5. Same as Table 4 but for the relatively warm Arctic conditions (surface temperature of −15 C).
Table 5. Same as Table 4 but for the relatively warm Arctic conditions (surface temperature of −15 C).
ScenarioResults FB ice TB Both Input Param. ρ ice ρ snow T sfc All Model Param.All Param.
I h i 11.88 % 1.93 % 11.48 % 4.23 % 0.20 % 3.53 % 5.66 % 12.80 %
h s 25.35 % 6.56 % 25.18 % 9.29 % 7.81 % 12.02 % 17.39 % 30.68 %
II h i 7.45 % 2.20 % 7.78 % 5.65 % 0.71 % 3.37 % 6.68 % 10.13 %
h s 22.14 % 24.92 % 33.74 % 16.77 % 16.92 % 38.29 % 44.08 % 54.91 %
III h i 6.07 % 5.24 % 7.97 % 6.91 % 1.41 % 4.92 % 8.71 % 12.52 %
h s 1.37 % 26.78 % 26.76 % 1.56 % 15.81 % 25.13 % 30.02 % 43.55 %
IV h i 4.16 % 8.81 % 9.80 % 8.21 % 1.85 % 6.68 % 11.10 % 15.36 %
h s 1.58 % 30.01 % 29.73 % 3.19 % 14.83 % 22.75 % 28.72 % 44.01 %
V h i 3.42 % 2.37 % 4.17 % 3.87 % 0.48 % 3.14 % 5.04 % 6.73 %
h s 9.84 % 8.11 % 12.65 % 11.06 % 6.65 % 10.70 % 16.96 % 21.27 %
VI h i 3.46 % 5.93 % 6.86 % 6.88 % 0.67 % 4.63 % 8.67 % 16.50 %
h s 1.59 % 33.68 % 33.49 % 3.25 % 11.64 % 26.30 % 30.62 % 71.31 %
VII h i 4.85 % 54.97 % 56.27 % 14.36 % 2.04 % 12.50 % 32.46 % 86.22 %
h s 5.17 % > 100 % >100% 19.32 % 13.48 % 36.50 % 76.16 % >100%
VIII h i 3.07 % 14.08 % 14.77 % 11.86 % 1.36 % 9.26 % 17.23 % 28.86 %
h s 5.04 % 59.97 % 61.86 % 19.64 % 13.64 % 39.41 % 51.66 % >100%
Table 6. Same as Table 4 but for the data synergy between T B and F B s n o w .
Table 6. Same as Table 4 but for the data synergy between T B and F B s n o w .
ScenarioResults FB ice TB Both Input Param. ρ ice ρ snow T sfc All Model Param.All Param.
I h i 32.96%35.71%33.14%32.90%37.12%34.51%33.24%33.61%
h s 61.79 % 60.40 % 61.92 % 54.80 % 62.52 % 57.54 % 57.13 % 62.04 %
II h i 28.11 % 24.30 % 28.79 % 28.10 % 24.81 % 24.95 % 28.48 % 29.83 %
h s 54.03 % 55.47 % 57.48 % 54.53 % 51.12 % 56.77 % 58.72 % 64.14 %
III h i 5.60 % 4.30 % 6.98 % 6.23 % 3.50 % 4.60 % 8.61 % 11.31 %
h s 1.75 % 9.99 % 10.08 % 2.08 % 4.59 % 10.69 % 11.84 % 16.09 %
IV h i 2.43 % 4.34 % 4.90 % 4.69 % 3.46 % 3.91 % 7.08 % 8.61 %
h s 1.96 % 6.71 % 6.83 % 3.78 % 1.86 % 6.06 % 7.28 % 10.28 %
V h i 18.67 % 22.35 % 22.43 % 20.13 % 13.16 % 22.92 % 23.19 % 24.56 %
h s 38.84 % 40.67 % 45.13 % 43.55 % 21.27 % 42.70 % 46.33 % 46.04 %
VI h i 25.31 % 26.60 % 27.49 % 27.18 % 25.27 % 26.33 % 29.40 % 30.88 %
h s 35.56 % 37.84 % 39.19 % 36.18 % 35.95 % 37.14 % 40.38 % 43.25 %
VII h i 2.26 % 7.01 % 7.26 % 4.97 % 4.32 % 4.60 % 8.02 % 10.54 %
h s 1.59 % 9.30 % 9.81 % 3.02 % 2.36 % 6.01 % 7.17 % 11.70 %
VIII h i 0.99 % 3.99 % 4.16 % 3.58 % 2.38 % 2.93 % 5.29 % 6.68 %
h s 1.74 % 7.73 % 7.98 % 6.31 % 1.25 % 5.67 % 8.75 % 11.63 %
Table 7. Same as Table 5 but for the data synergy between T B and F B s n o w .
Table 7. Same as Table 5 but for the data synergy between T B and F B s n o w .
ScenarioResults FB ice TB Both Input Param. ρ ice ρ snow T sfc All Model Param.All Param.
I h i 22.72 % 21.43 % 24.07 % 19.07 % 20.41 % 22.45 % 22.81 % 26.71 %
h s 33.59 % 39.53 % 38.73 % 38.04 % 38.80 % 33.12 % 34.10 % 40.09 %
II h i 35.11 % 31.87 % 36.07 % 32.50 % 30.59 % 33.18 % 36.70 % 38.70 %
h s 65.16 % 57.95 % 63.67 % 57.63 % 58.09 % 54.32 % 57.42 % 59.10 %
III h i 22.46 % 26.36 % 27.92 % 11.34 % 25.48 % 20.72 % 29.56 % 30.48 %
h s 33.11 % 43.44 % 42.57 % 16.17 % 38.12 % 41.13 % 44.20 % 47.93 %
IV h i 3.02 % 13.61 % 13.57 % 5.57 % 10.65 % 10.42 % 15.63 % 19.12 %
h s 1.14 % 20.55 % 20.21 % 2.12 % 12.17 % 16.45 % 19.73 % 26.72 %
V h i 18.29 % 22.08 % 22.64 % 18.82 % 20.58 % 22.93 % 24.64 % 25.65 %
h s 45.23 % 47.91 % 49.65 % 48.24 % 44.07 % 46.84 % 49.28 % 47.62 %
VI h i 25.35 % 29.36 % 29.88 % 26.29 % 25.77 % 28.88 % 30.74 % 33.50 %
h s 34.60 % 40.92 % 41.90 % 34.23 % 35.37 % 40.10 % 42.10 % 47.46 %
VII h i 1.74 % 12.16 % 12.29 % 3.94 % 6.80 % 8.03 % 12.49 % 16.36 %
h s 1.79 % 16.01 % 15.85 % 4.06 % 5.40 % 10.75 % 13.04 % 19.90 %
VIII h i 1.03 % 7.99 % 7.80 % 3.77 % 3.83 % 6.63 % 8.67 % 11.62 %
h s 1.69 % 15.37 % 14.99 % 6.19 % 3.80 % 13.04 % 14.58 % 21.29 %

Share and Cite

MDPI and ACS Style

Xu, S.; Zhou, L.; Liu, J.; Lu, H.; Wang, B. Data Synergy between Altimetry and L-Band Passive Microwave Remote Sensing for the Retrieval of Sea Ice Parameters—A Theoretical Study of Methodology. Remote Sens. 2017, 9, 1079. https://doi.org/10.3390/rs9101079

AMA Style

Xu S, Zhou L, Liu J, Lu H, Wang B. Data Synergy between Altimetry and L-Band Passive Microwave Remote Sensing for the Retrieval of Sea Ice Parameters—A Theoretical Study of Methodology. Remote Sensing. 2017; 9(10):1079. https://doi.org/10.3390/rs9101079

Chicago/Turabian Style

Xu, Shiming, Lu Zhou, Jiping Liu, Hui Lu, and Bin Wang. 2017. "Data Synergy between Altimetry and L-Band Passive Microwave Remote Sensing for the Retrieval of Sea Ice Parameters—A Theoretical Study of Methodology" Remote Sensing 9, no. 10: 1079. https://doi.org/10.3390/rs9101079

APA Style

Xu, S., Zhou, L., Liu, J., Lu, H., & Wang, B. (2017). Data Synergy between Altimetry and L-Band Passive Microwave Remote Sensing for the Retrieval of Sea Ice Parameters—A Theoretical Study of Methodology. Remote Sensing, 9(10), 1079. https://doi.org/10.3390/rs9101079

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