Sea ice emissivity modelling
With increased interest in sea ice and its effects on the global climate, efficient methods are required to monitor both its extent and exchange processes.
Satellite-mounted, microwave radiometers, such SSMI, AMSR and AMSU, are an ideal tool for the task because they can see through cloud cover, and they have frequent, global coverage.
A passive microwave instrument detects objects through emitted radiation since different substance have different emission spectra.
To detect sea ice more efficiently, there is a need to model these emission processes.
The interaction of sea ice with electromagnetic radiation in the microwave range is still not well understood.
[1][2][3] In general is collected information limited because of the large-scale variability due to the emissivity of sea ice.
[4] Satellite microwave data (and visible, infrared data depending on the conditions) collected from sensors assumes that ocean surface is a binary (ice covered or ice free) and observations are used to quantify the radiative flux.
During the melt seasons in spring and summer, sea ice surface temperature goes above freezing.
Thus, passive microwave measurements are able to detect rising brightness temperatures, as the emissivity increases to almost that of a blackbody, and as liquid starts to form around the ice crystals, but when melting continues, slush forms and then melt ponds and the brightness temperature goes down to that of ice free water.
[5] As established in the previous section, the most important quantity in radiative transfer calculations of sea ice is the relative permittivity.
Since it is not just the relative composition that is important, but also the geometry, the calculation of effective permittivities introduces a high level of uncertainty.
Vant et al. [6] have performed actual measurements of sea ice relative permittivities at frequencies between 0.1 and 4.0 GHz which they have encapsulated in the following formula:
This empirical model shows some agreement with dielectric mixture models based on Maxwell's equations in the low frequency limit, such as this formula from Sihvola and Kong
The two formulas, while they correlate strongly, disagree in both relative and absolute magnitudes.
[2] Pure ice is an almost perfect dielectric with a real permittivity of roughly 3.15 in the microwave range which is fairly independent of frequency while the imaginary component is negligible, especially in comparison with the brine which is extremely lossy.
[8] Meanwhile, the permittivity of the brine, which has both a large real part and a large imaginary part, is normally calculated with a complex formula based on Debye relaxation curves.
[8] When scattering is neglected, sea ice emissivity can be modelled through radiative transfer.
The diagram to the right shows a ray passing through an ice sheet with several layers.
Each layer is characterized by its physical properties: temperature, Ti, complex permittivity,
Since we assume plane-parallel geometry, all reflected rays will be at the same angle and we need only account for radiation along a single line-of-sight.
The most important quantity in this calculation, and also the most difficult to establish with any certainty, is the complex refractive index, ni.
Emissivity calculations based strictly on radiative transfer tend to underestimate the brightness temperatures of sea ice, especially in the higher frequencies, because both included brine and air pockets within the ice will tend to scatter the radiation.
[9] Indeed, as ice becomes more opaque with higher frequency, radiative transfer becomes less important while scattering processes begin to dominate.
Scattering in sea ice is frequently modelled with a Born approximation [10] such as in strong fluctuation theory.
The Microwave Emission Model of Layered Snowpack (MEMLS) [13] uses a six-flux radiative transfer model to integrate both the scattering coefficients and the effective permittivities with scattering coefficients calculated either empirically or with a distorted Born approximation.
Mills and Heygster,[2] for instance, show that sea ice ridging may have a significant effect on the signal.
In addition to ridging, surface scattering from smaller-scale roughness must also be considered.
Since the microstructural properties of sea ice tend to be anisotropic, permittivity is ideally modelled as a tensor.
This anisotropy will also affect the signal in the higher Stokes components, relevant for polarimetric radiometers such as WINDSAT.
Both a sloping ice surface, as in the case of ridging—see polarization mixing, [1] as well as scattering, especially from non-symmetric scatterers, [14] will cause a transfer of intensity between the different Stokes components—see vector radiative transfer.
