InTech uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Earth and Planetary Sciences » Geology and Geophysics » "Advances in Geoscience and Remote Sensing", book edited by Gary Jedlovec, ISBN 978-953-307-005-6, Published: October 1, 2009 under CC BY-NC-SA 3.0 license. © The Author(s).

Chapter 6

Some Features of The Volume Component of Radar Backscatter From Thick and Dry Snow Cover

By Boris Yurchak
DOI: 10.5772/8339

Article top

Overview

Scheme of sounding for flat surface and plane wavefront
Figure 1. Scheme of sounding for flat surface and plane wavefront
To an assessment of sphericity impact
Figure 2. To an assessment of sphericity impact
Schematic of sounding in nadir direction
Figure 3. Schematic of sounding in nadir direction
Estimation of complete burial of an increment of pulse due to the sphericity of the Earth
Figure 4. Estimation of complete burial of an increment of pulse due to the sphericity of the Earth
Illustration of the flattening of a spherical wavefront within a snow slab
Figure 5. Illustration of the flattening of a spherical wavefront within a snow slab
Modes of the volume scattering; a) depth-limited, b) beamwidth-limited
Figure 6. Modes of the volume scattering; a) depth-limited, b) beamwidth-limited
Assessment of the complete burial of radial size of the probing pulse
Figure 7. Assessment of the complete burial of radial size of the probing pulse
Relationships between the penetration depth (D
							p), penetration path (L
							
								p
							) and propagation depth (
								
									
										
											
												D
												p
												∗
											
										
									
								
								
							) in a snow slab
Figure 8. Relationships between the penetration depth (D p), penetration path (L p ) and propagation depth ( D p ∗ ) in a snow slab
An assessment of the scattering volume for a fully scattering snow slab
Figure 9. An assessment of the scattering volume for a fully scattering snow slab
Cross-section of the scattering volume for a spherical surface and wavefront. The position of the flattened wavefront at the moment when its central point (a’) reaches the penetration depth is depicted with a convex dashed line
Figure 10. Cross-section of the scattering volume for a spherical surface and wavefront. The position of the flattened wavefront at the moment when its central point (a’) reaches the penetration depth is depicted with a convex dashed line
Plot of the normalized effective depth of snow sounding versus the snow optical thickness
Figure 11. Plot of the normalized effective depth of snow sounding versus the snow optical thickness
The angular dependence of the normalized effective depth of snow sounding for the “thick” (nSD>>1) and “intermediate” regimes
Figure 12. The angular dependence of the normalized effective depth of snow sounding for the “thick” (nSD>>1) and “intermediate” regimes
Illustration of the expected relative spatial behavior of the backscatter coefficient under a linear dependence of the snow depth on a distance
Figure 13. Illustration of the expected relative spatial behavior of the backscatter coefficient under a linear dependence of the snow depth on a distance
Transect across hill (straight line crossing the dark spot on the SAR signature) in the north-eastern region of Greenland (left), and the corresponding backscatter (solid line, convex downwards curve) and elevation (dotted line, convex upwards curve) profiles (right)
Figure 14. Transect across hill (straight line crossing the dark spot on the SAR signature) in the north-eastern region of Greenland (left), and the corresponding backscatter (solid line, convex downwards curve) and elevation (dotted line, convex upwards curve) profiles (right)
Illustration of the backscatter coefficient values (left) range for an arbitrary transect (straight line on RadarSAT-1 SAR signature, right) from the percolation to the dry snow areas within the east margin of the Greenland Ice Sheet (directions are from left to right for the plot of the backscatter profile and contra versa for the SAR signature; date of data acquisition: November 2, 2007; data granule ID: R1_62609_SWB_271)
Figure 15. Illustration of the backscatter coefficient values (left) range for an arbitrary transect (straight line on RadarSAT-1 SAR signature, right) from the percolation to the dry snow areas within the east margin of the Greenland Ice Sheet (directions are from left to right for the plot of the backscatter profile and contra versa for the SAR signature; date of data acquisition: November 2, 2007; data granule ID: R1_62609_SWB_271)
The dependence of the backscatter coefficient on particle size in accordance with the A-U model for C-band radar and snow depth greater than the penetration depth (“thick” snow regime)
Figure 16. The dependence of the backscatter coefficient on particle size in accordance with the A-U model for C-band radar and snow depth greater than the penetration depth (“thick” snow regime)
Simplified slice approach scheme. The arbitrary position of a single slice is shown on the cross section of the scattering volume
Figure 17. Simplified slice approach scheme. The arbitrary position of a single slice is shown on the cross section of the scattering volume
Deviation factor versus particle radius variation coefficient for different deviations of Poisson index χ from the Poisson law
Figure 18. Deviation factor versus particle radius variation coefficient for different deviations of Poisson index χ from the Poisson law
Typical backscatter coefficient profile across a transect within the margin of the Greenland Ice Sheet (in the vicinity of the Swiss camp)
Figure 19. Typical backscatter coefficient profile across a transect within the margin of the Greenland Ice Sheet (in the vicinity of the Swiss camp)
The normalized snow depth profile derived from the backscatter coefficient profile in Figure 19.
Figure 20. The normalized snow depth profile derived from the backscatter coefficient profile in Figure 19.
Study area #3 (box with center: 75.070N and 54.300W) within the Greenland ice sheet margin with several transects normal to the coastline (left). Normalized snow depth profiles (right) for transect 3a (lower line in the study area) derived from RadarSAT-1 data. Distance scale: 1 sample=50 m. Solid line is nSD profile on 11/16/1997, dotted line is nSD profile on 11/12/2007
Figure 21. Study area #3 (box with center: 75.070N and 54.300W) within the Greenland ice sheet margin with several transects normal to the coastline (left). Normalized snow depth profiles (right) for transect 3a (lower line in the study area) derived from RadarSAT-1 data. Distance scale: 1 sample=50 m. Solid line is nSD profile on 11/16/1997, dotted line is nSD profile on 11/12/2007
Dimensionless spectral characteristics of the backscatter for several values of the slant snow depth based on the empirical spectral dependence of the penetration depth
Figure 22. Dimensionless spectral characteristics of the backscatter for several values of the slant snow depth based on the empirical spectral dependence of the penetration depth
Slant snow depth as a function of the magnitude difference between the dimensionless spectral components of the backscatter coefficients at 5 and 14 GHz
Figure 23. Slant snow depth as a function of the magnitude difference between the dimensionless spectral components of the backscatter coefficients at 5 and 14 GHz
Spectral pattern of the experimental data of (Baumgartner et al., 1999) for an incidence angle of 400 for two test sites in Greenland
Figure 24. Spectral pattern of the experimental data of (Baumgartner et al., 1999) for an incidence angle of 400 for two test sites in Greenland
Theoretical spectral dependence under the Rayleigh approach of the penetration depth for a temperature of -150C (our calculations) with comparison of the empirical relationship (7.4)
Figure 25. Theoretical spectral dependence under the Rayleigh approach of the penetration depth for a temperature of -150C (our calculations) with comparison of the empirical relationship (7.4)
Comparison of the spectral dependence of the penetration depth calculated with Mie formulas for a temperature of -10С (Ulaby et al., 1986) with our calculations under the Rayleigh approach for a temperature of -150С
Figure 26. Comparison of the spectral dependence of the penetration depth calculated with Mie formulas for a temperature of -10С (Ulaby et al., 1986) with our calculations under the Rayleigh approach for a temperature of -150С
Illustration of the reproducibility of the snow backscatter uniform profile taking into account only attenuation (bold dash line) and finite probing pulse length normalized to the penetration depth (parameter
								
									
										
											
												
													
														(
														
															
																S
																γ
															
														
														)
													
												
												
													γ
													
													
													1
												
											
											=
											
												R
												Δ
											
											
											
											1
										
									
								
								
							).
Figure 27. Illustration of the reproducibility of the snow backscatter uniform profile taking into account only attenuation (bold dash line) and finite probing pulse length normalized to the penetration depth (parameter ( S γ ) γ 1 = R Δ 1 ).

Some Features of The Volume Component of Radar Backscatter From Thick and Dry Snow Cover

Boris Yurchak1

1. Introduction

Radar monitoring of thick snow cover in polar regions with optical thickness (a product of the depth of the snow and the extinction coefficient) of order 1 or more from elevated above-ground and space platforms is of great importance for registration and for understanding glaciology processes caused by climatic change. Applicable to this issue, the volume component of the backscatter coefficient has a notable contribution to the total backscatter (Noveltis, 2005). Although the idea of the radar-cross section (RCS) term and its derivatives, like the backscatter coefficient, aims to separate as far as possible the sensor (radar) and target parameters, this distinction can rarely be fulfilled in the case of spatially extended geophysical targets (SEGT), such as atmospheric clouds and rain, as well as the thick snow cover that is the focus of this paper. Due to this feature, the analysis of backscatter from SEGT strongly depends on the relationship between the technical and physical-geometrical properties of the radar and target, respectively. The main parameters that govern the radar-target configuration for snow sounding are wavelength, antenna characteristics, pulse duration, sounding direction, extinction coefficient and the geometrical depth of the snow. A correct assessment of the volume component of the backscatter coefficient and an understanding of the realm of applicability of any backscattering model is possible only when the size of the scattering volume within a snow slab is known. This parameter depends on the factors mentioned above and should be distinguished for different situations in the practice of radar sounding of snow-covered terrain by Synthetic Aperture Radar (SAR), a scatterometer or an altimeter. This work attempts, at first, to determine the radar-target configurations inherent to volume scattering estimations of the thick snow cover under the different radar applications mentioned above. Next, we analyze the range of applications of the incoherent approach for backscatter magnitude estimation, currently one of the main techniques for snow characteristics assessment. The simplest incoherent approach, based on the so-called “particle” or “discrete” approximation, leads to the dependence of backscatter on the sixth moment of the particle size distribution function (PSDF) and the mean amount of particles in the scattering volume (Siegert & Goldstein, 1951, Battan, 1959, Ulaby et al., 1982). For a medium with losses (such as a thick snow slab) the modification of this approach is referred to as the semi-empirical model (Attema & Ulaby, 1978, , Ulaby et al.1982, Ulaby et al., 1996). There is some evidence of less backscatter occurring than expected by virtue of the conventional reflectivity factor and the backscatter coefficient for incoherent scatter from homogeneous thick snow slab (e.g., Rott et al., 1993). For dense media (mean distance between particles is less than the wavelength) there are many references in the literature that the backscatter is determined by media inhomogeneities (Naito & Atlas, 1967, Gossard & Srauch, 1983, Fung, 1994). To attempt to evaluate the contributions of these inhomogeneities, the radiative transfer (RT) and the dense media radiative transfer (DMRT) models were developed (e.g., Ulaby et al., 1982, Fung, 1994, Tsang et al., 2007). A review of these models and their modifications provided by Noveltis (2005) stated that in certain aspects, these approaches showed some successes. However, these models have not the close-form solutions that make difficult to use them and analyze results obtained. Kendra et al. (1998), based on experiments with artificial snow of varying depths, concluded that both conventional and dense-medium radiative transfer models fail to adequately explain the observed results.

Finding an appropriate and relatively simple approach to calculate the volume component of backscatter to explain the observable deviations from the classical (incoherent) model is therefore a relevant task. The majority of previous studies, focusing on measurements of a dry snow at temperate latitudes, where the snow depth seldom exceeds 1 m and the corresponding optical thickness is much less than 1, have led to the conclusion that such an approach is not feasible because of the much weaker interaction of the electromagnetic radiation of radar wavebands with ice particles within the snow compared with the backscatter from the soil beneath. The current study focuses on the case of the Greenland ice sheet, however, where the depth of snow significantly exceeds that found in temperate latitudes. This favorable condition, in conjunction with the recent finding that the incoherent approach is only a specific case of a more comprehensive description of electromagnetic wave interactions with spatially-extended individual scatterers (Yurchak, 2009), provides a reason for more detailed investigations of the possibilities of the semi-empirical model for interpreting the observable features of radar backscatter from thick snow cover.

2. Condition of complete burial of the probing pulse into snow medium

To estimate the volume component of backscatter from a snow slab, it is first necessary to understand whether the probing pulse is completely buried within the snow slab or whether the illuminated volume is only part of the probing pulse volume. In weather radar meteorology, applicable to rain and clouds, this problem is known as “partial (or incomplete) beam filling,” with a corresponding factor included in the weather radar equation (e.g., Clift, 1985). For sounding thick snow cover, this issue is practically not discussed and is different compared with the sounding of meteorological targets. Complete burial depends on the sounding configuration and snow slab depth. The main criterion for complete burial is the location of the pulse scattering volume with angular size equal to the angular antenna beam width and with spatial length equal to one half of the actual spatial duration of the transmitted pulse into a snow layer. Thus, for complete burial, the angular and radial sizes of the backscatter volume should be matched with the snow slab depth, which is usually known only roughly for a particular geographical region. This condition poses problems for practical applications. Nevertheless, to better understand possible situations where complete burial is feasible, the requirements for the angular antenna beam size and the probing pulse duration will be analyzed separately for different major sounding configurations.

2.1. Condition of complete burial of the angular (transverse) size of the pulse volume

2.1.1. Oblique sounding, flat surface, plane wavefront

The scheme of this configuration is shown in Figure 1 Here, and everywhere below, the snow slab is assumed to have limited depth and an unbounded horizontal extent.

media/image1.jpg

Figure 1.

Scheme of sounding for flat surface and plane wavefront

Mathematically, the condition for complete burial in this case is:

where d b is the burial depth of the advance point (p) of pulse edge (leading or trailing) when the lagging point (c) reaches the surface and h is the snow depth. As follows from the geometry of Figure 2.1:

db=cg=pc×sincpg cpg=θ pc2R0cosθtgϕ0.52 , and hence:
db=2R0cosθtgϕ0.52sinθ
(2)

For φ0.5<<1 one can state tgϕ0.52ϕ0.52 , and condition (2.1) has the form:

ϕ0.5hR01tgθ
(3)

If h~10 m and R 0 =800 km, then ϕ0.51.25tgθ105 and for θ~400, for example, the burial condition (2.1) is satisfied at roughly φ 0.5<<0.0010. Obviously, this condition can not be fulfilled for any space-based radar system with a real aperture, but it is possible for the SAR if one assumes that the effective synthetic aperture radar beam illuminating an element of spatial resolution Δx on a flat surface can be equal to an extremely narrow, pencil-like beam with ϕ0.5~ΔxR . For example, RadarSAT-1 has Δx ~12.5 m with an orbit height of approximately 800 km and thus, φ 0.5~0.9*10-3 degrees. If one assumes that the typical size of the main lobe of the conventional antenna pattern is of order ~10, the complete burial of the transverse size of the pulse volume is possible for hR0ϕ0.5tgθ~2102tgθ . This condition can definitely be fulfilled for airborne and above-surface elevated radars. The estimate obtained above should be considered only a rough approximation, because at nadir sounding (θ=0) and unbounded horizontal extent of snow slab, the wavefront can not be considered planar for the assessment provided. A more precise estimate is given below.

2.1.2. Oblique sounding, flat surface, spherical wavefront

A real wavefront within a snow slab has a spherical shape. The impact of this shape on the estimation of d b as the look angle θ (Figure 1.) decreases begins at the moment when the line tangent to the spherical front at point p coincides with the horizontal line pg, Figure 2. This situation takes place when θ becomes equal to the angle between pc and the tangent:

media/image13.jpg

Figure 2.

To an assessment of sphericity impact

As follows from geometrical sketch, the auxiliary angle β is equal to: β=πϕ0.52 and, thus, γ=π2β=ϕ0.52 . Therefore, the impact of sphericity should be taken into account only for small look angles, when

Thus, the estimate (2.3), provided in the previous subsection, is valid for θϕ0.52 . An estimate for configurations close to nadir sounding will be carried out in the next subsection.

2.1.3. Nadir sounding, flat surface spherical wavefront

media/image18.jpg

Figure 3.

Schematic of sounding in nadir direction

This configuration is shown schematically in Figure 3 and the corresponding condition for complete burial is:

From the geometrical configuration and φ 0.5<<1, it follows that:

db(flt)=R0(1cosϕ0.521)=R02sin2ϕ0.54cosϕ0.5218R0ϕ0.52
(7)

This estimate can be accepted for elevated platforms (for example, airborne radars and those mounted above the surface) but needs to be analyzed further for space platforms due to the sphericity of the surface, as the corresponding area illuminated on the Earth’s surface by conventional radar is large.

2.1.4. Nadir sounding, spherical surface, spherical wavefront

In the spherical surface approach, the depth of complete burial in altimeter mode is more than that for a flat surface described above by an increment Δz in the center of the beam. This effect is illustrated in Figure 4, and can be written as:

db(sp)=db+Δzh
(8)
media/image22.jpg

Figure 4.

Estimation of complete burial of an increment of pulse due to the sphericity of the Earth

Based on the geometrical relations of the sketch, the following system of equations may be composed

{yR0+q=tgϕ0.52RE2=(REq)2+y2
(9)

where q=Δzcosϕ0.52 , and y is the distance between the vertical axis and the beam periphery touching the spherical surface. Solving this system, the increment due to the sphericity is equal to:

q=RER0tg2ϕ0.52(1+tg2ϕ0.52)[1±1(1+tg2ϕ0.52)(R0tgϕ0.52)2(RER0tg2ϕ0.52)2]
(10)

where R E is the Earth’s radius. Taking into account that ϕ0.52 << 1, RER01 , and that the minus sign should be taken before the square root for physical reasons, the above formulae can be simplified significantly, yielding an increment equal to:

Δz12R0RER0(ϕ0.52)2(cosϕ0.52)1
(11)

A numerical assessment under the parameters values listed above gives Δz~5 m. The estimation of the parameter db(sp) can be represented in the form:

db(sp)=R0(1cosϕ0.521)+12R0RER0(ϕ0.52)21cosϕ0.5218R0ϕ0.52(1+R0RE)
(12)

This relationship coincides with the estimate obtained by Barrick (1972) for a spherical wave increment over the spherical mean sea surface. For R 0=800 km, R E=6400 km and ϕ0.5=10 , one calculates db(sp) ≈45 m. This result means that the Earth’s sphericity increases the complete burial parameter compared with the flat surface case (see 2.7) by a factor of ~( 1+R0RE ). For R0RE18 the complete burial parameter equals to 1.125, i.e., a 12.5% increase.

Based on known assessments of penetration depth of 2-5 m for Ku-band (e.g., Davis, 1996), the wave sphericity in practice does not allow for consideration of the complete burial of the transversal size of the probing pulse for a space-based altimeter and scatterometer. This sphericity causes the area illuminated on the surface to change as the probing pulse is buried into a snow slab. It can be easily shown that the radius of the illuminated area changes with the burial increment ΔR as rill2R0ΔR .

2.1.5. Wavefront flattening

The speed of an electromagnetic wave within a snow slab is lower than that in an air by ε , where ε’ is the real part of the dielectric permittivity of snow. For dry snow within a density range of 0.2-0.5 g/cm3, parameter ε’ changes roughly from 1.35 to 1.95 (Tiuri et al., 1984). Recently, a decrease in wave speed was also confirmed by direct measurements in snow (Scott et al., 2006). Due to the sphericities of the Earth and the wavefront, the paths of the wave for different rays into a snow slab for equal time are not the same that results from the distortion of the initial spherical wave front, as demonstrated in Figure 5.

media/image36.jpg

Figure 5.

Illustration of the flattening of a spherical wavefront within a snow slab

The peripheral ray (ray Ob, for example) has a different path into a snow slab as compared with the central ray (Og), which that has a portion ag extending into the snow. This feature causes the flattening of the spherical wavefront. Thus, the snow slab works like dielectric lens antenna (e.g., Lo & Lee, 1993). Let us make a rough estimate of the impact of this phenomenon on the increment Δz. For a spherical wavefront (dashed curve in Figure 2.5), the following equation is valid: Ob=Oa+au+ug. Dividing both sides of this equation by speed of the wave propagation, c, yields the time of the wavefront arrival at point g: tg(sp)=Obc=Oa+au+ugc .

Due to delays in snow in part of ag=au+ug the equation for the wavefront arrival at point g is: tg(flt)=Oac+au+ugv , where v is the wave speed in snow. Because v<c, tg(flt)tg(sp) , providing proof of this flattening. To rate this phenomenon, let us suppose that the flattening front crosses the vertical at point u. This condition implies that the wave arrival times at points u and b are the same. This yields the equations: Obc=Oac+auv and au=vc(ObOa)=vcag . Hence, the relative shortage (flattening) of the wavefront path in the center of the beam compared with propagation in the free space (air) is equal to: auag=vc=1ε , and the corrected increment Δz (2.11) should be decreased by ε ~1.2…1.4 times.

Moreover, due to the increase in snow density with depth, the dielectric permittivity also increases, additionally impacting the rays’ path configuration as takes place in a Luneburg lens. Thus, the actual meaning of the burial parameter lies between the two estimates carried out above:

db(flt)dbdb(sp)
(13)

Based on the above estimates, one can say that the surface sphericity practically contributes a relatively small amount to the estimate of the complete burial condition of the transverse dimension of the probing pulse, and consequently, the surface can be roughly considered as flat for any radar configurations regarding volume component assessment.

2.1.6. Modes of volume scattering for spherical wave

By analogy to spherical wave scattering from a surface (e.g., Moore & Williams, 1957) it is also reasonable to distinguish two modes of the volume component in the case of spherical wave scattering, as shown in Figure 6. For simplicity, the bounds of the snow slab are considered flat. Let us call the mode depth-limited (Figure 6a.) if at the moment when the central point of the spherical wavefront reaches the bottom of a snow slab or a penetration depth (D p ), the peripheral point of the wavefront, crossing the surface, is still within the beamwidth footprint on the surface. This condition takes place when:

hR012(ϕ0.52)2 if, or
DpR012(ϕ0.52)2
(15)
otherwise
media/image49.jpg

Figure 6.

Modes of the volume scattering; a) depth-limited, b) beamwidth-limited

Let us call the mode beamwidth-limited (Figure 6b.) if the peripheral points of the spherical wavefronts, crossing the surface, are beyond the beamwidth footprint on the surface for a significant part of the scattering volume. This situation takes place when:

ϕ0.522hR0 if or otherwise

Obviously, that the depth-limited mode is inherent to space based altimeter and scatterometr and the beamwidth-limited mode is realized for SAR. If the “significant part” in the definition of this mode is replaced with a more exact term, the “determined part,” the condition (2.15) can be modified and written in the form:

ϕ0.522hR0(1Δhh) if or
ϕ0.522DpR0(1ΔhDp)
(19)
otherwise

where Δhh or ΔhDp is the relative part of the scattering volume in which the beamwidth-limited mode is fulfilled. For example, if 90% of the scattering volume ( Δhh =0.9) is under the beamwidth-limited mode, the condition (2.15a) is ϕ0.520.2hR0 for

2.2. Condition of complete burial of radial size of the pulse volume

Suppose that some area of snow cover is illuminated by a radar located in orbit as shown in Figure 7. For simplicity, this scenario is depicted with a plane wavefront and flat surface.

media/image57.jpg

Figure 7.

Assessment of the complete burial of radial size of the probing pulse

Let us denote the height of the orbit as R 0 , the distance from the radar to the center of the area as R, the anglular width of the main lobe of the antenna pattern as φ 0.5 and, the look angle as α. In addition, we assume that the cross-section of the main lobe is the circular. The illuminated area is assumed to be flat horizontally, and the bottom boundary surface of the snow slab is also flat and parallel to the top surface. Thus, the incidence angle, θ is equal to the look angle, α (θ=α).

The main scaling parameters determining radial propagation are:

1) The radial length of the pulse scattering volume, equal to half of its spatial extent in the snow medium:

hsct(p)=vτ02
(21)

where v=cε is the wave propagation speed in snow, c is the wave propagation speed in the air, ε is the real part of the snow permittivity and τ 0 is the duration of the probing pulse;

2) A one-way path in the snow, where the incidence power is decreased by “e” times, usually called the “penetration depth” in the literature

where k e is the extinction coefficient, which characterizes the attenuation properties of the medium due to scattering and absorption. One important remark is necessary. Since radar sounding of snow is often performed in off-nadir mode, and the main lobe of the antenna pattern has a finite angular size, it should be

media/image62.jpg

Figure 8.

Relationships between the penetration depth (D p), penetration path (L p ) and propagation depth ( Dp ) in a snow slab

underlined that the extinction coefficient and the “penetration depth” are measured along the direction of the wave (ray) propagation and not only to the vertical. Therefore, calling the term ke1 as the “…depth” should be considered a little bit confusing. A more appropriate name for this term, from a physical point of view, is the “penetration path” (L p ), keeping in mind that the penetration depth (D p) is its value in the vertical direction. The depth reached by an electromagnetic wave propagating at an angle θ to the vertical and attenuating by “e” times represents the vertical component of the penetration path, and can be called the propagation depth (Dp) , as illustrated in Figure 8. Thus, from this point onward in this paper, the following definitions and relationships are used:

Penetration path: Lp=ke1

Propagation depth: Dp=Lpcosθ

Penetration depth: Dp=ke,1

where ke, is the extinction coefficient in the vertical direction. In general, in inhomogeneous medium, the penetration path is a function of distance and direction Lp=f(R,θ,ψ) , where ψ is the azimuthal angle. For homogenous medium Dp=Lp . We consider the probing pulse to be short if hsct(p)Lp ; otherwise it is considered long. Obviously, the mode of sounding with current wave (CW) always belongs to the long probing pulse configuration. Also one should distinguish a case with fully scattering snow slabs, when the geometrical snow depth is equal to or less than the propagation depth ( hDp ). Otherwise, one has a case of partially (because the wave does not penetrate to the bottom) scattering slow slab. The condition for complete burial, as it follows from Figure 7, is

hsct(p)hθΔh
(23)

where hθ=hcosθ is the slant snow depth and Δh=2R0tgϕ0.52tgθcosθR0ϕ0.5tgθcosθ .

All of the above conditions are summarized in Table 1.

#ConditionDescriptionComments
1Short/long probing pulse (P.P.)Complete burial is satisfied for:
short P.P. and fully scattering snow slab;
long P.P. and partially scattering snow slab
2or
Fully/partially scattering snow slab
3Complete burial of the radial size of the probing pulse
Table 2.1. Conditions for complete burial of the radial size of the probing pulse

Table 1.

Conditions for complete burial of the radial size of the probing pulse

Concluding remarks:

If the conditions for complete burial of the radial size of the probing pulse are completed, it can at least be said that one sample of the return signal is formed by the scattering volume

where hDp is the illuminated base of the scattering volume. If, in addition, the conditions for complete burial of the transverse size of the probing pulse (see previous subsection) are also satisfied, that the scattering volume is determined by the entire pulse volume:

Otherwise, under hsct(p)hθΔh ,

Vsct=hsct(p)Aill
(26)

2.3. An assessment of the scattering volume under incomplete burial condition

If the parameters of the probing pulse and of the snow slab do not satisfy the conditions presented in the table above, the radial size of the scattering volume is determined by the geometry of the snow slab. Several practically important cases for practical application are discussed below.

2.3.1. Flat surface, plane wavefront, long probing pulse and fully scattering snow slab

In this scenario, the illuminated area, A ill on the snow cover (Figure 1.) is an ellipse, with the minor semi-axis equal to the radar beam cross-section radius Aill and major semi-axis the same divided by the cosine of the incidence angle: Vsct=Vp=hsct(p)π(R0ϕ0.52)2 . Thus:

ΔR12R0(ϕ0,52)2
(27)

If one uses the height of the satellite orbit Vsct(ΔR)hsct(p)2πR0ΔR , (2.23) transforms to:

rmin=(Rϕ0.52)
(28)

The pattern of the scattering volume for this case is depicted in Figure 9, where the scattering volume is bounded by the surface and bottom planes.

media/image90.jpg

Figure 9.

An assessment of the scattering volume for a fully scattering snow slab

In this case, the size of the scattering volume is the sum of the volume of a circular cylinder V 0 with cylindrical element d 0 and two similar volumes of the truncated cylinders V 1 that mutually add up to a completed circular cylinder with cylindrical element d 1 (the bases of all volumes are the same and equal to rmjr=rmincosθ ):

Aill=π(Rϕ0.52)21cosθ
(29)

Because R0=Rcosθ , we have, taking expression (2.25) into account:

Aill=π(R0ϕ0.52)21cos3θ
(30)

Thus, the equivalent scattering volume is an elliptic cylinder with a base equal to the illuminated area on a surface and with a height element equal to the slab depth (under the assumption that the slab depth is less than the half of the spatial duration of the probe radar pulse). This result was obtained by Matzler (1987).

The duration of the SAR probing pulse duration is equal to several tens of microseconds (e.g., 42 μs for RadarSAT-1 and 37.1 μs for ERS-1). Due to frequency chirp, the compressed probing pulse duration decreases by many times, resulting in a volume radial size equal to only several meters (5-13 m for RadarSAT-1 and 9.7 m for ERS-1, for example). Numerical data are provided based on Alaska Satellite Facility documents (“RadarSAT-1 Standard Beam SAR Images”, 1999; and “ERS-1 and ERS-2 SAR Images”, 1996).

media/image95.jpg

Figure 10.

Cross-section of the scattering volume for a spherical surface and wavefront. The position of the flattened wavefront at the moment when its central point (a’) reaches the penetration depth is depicted with a convex dashed line

Because SAR images represent the backscatter pulse train from the entire path of a wave into a snow slab for any look angle as a point on a surface of some mean brightness (return power), the radial size of the scattering volume can be considered equal to either the slant size of the snow layer or the penetration path (whichever is smaller). That is, although the SAR pulse has a finite compressed spatial duration; it works like a long pulse due to the absence of radial discrimination in the sense of conventional radar terminology.

2.3.2. Long probing pulse, spherical wave, spherical surface, partially scattering slab

For this case, the bottom of the scattering volume is the geometrical placement of points located at a distance L p from the snow surface along the family of rays within the solid angle of the main antenna lobe, as shown in Figure 10.. The defining feature of this configuration is that the central and peripheral points of the spherical wavefront pass the length L p at different times, resulting in different radii of curvature for the top and bottom bounded surfaces.

3. Volume component of the backscatter coefficient

3.1. Semi-empirical model for flat surface, plane wave and long probing pulse

Based on the considerations provided above, one can formulate an estimation of the volume component of the backscatter within the framework of the so-called semi-empirical model (Attema & Ulaby, 1978), as for the case of incomplete burial of the probing pulse with a flat surface, plane wave and a long probing pulse, which better fits the conditions for sounding of a thick snow slab with SAR.

The total radar backscatter from the illuminated area is composed of four components (Fung, 1994):

where σ as is the radar cross-section (RCS) due to backscatter from the air/snow (top) interface, σ s is the RCS due to backscatter from the snow volume, σ g is the RCS due to backscatter from the snow/ground (bottom) interface and σ gv is the RCS due to rescattering between ground and snow volume irregularities. As was summarized by Koskinen (2001), only the volume component is essential for dry snow. Therefore the essential portion of the total backscatter consists of the volume component and the contribution of the ground beneath:

Vsct=V0+2V1=Aillcosθ(d0+d1)
(32)

Obviously, the radar becomes sensitive to the properties of the snow only when d0+d1=hθ . In accordance with Ulaby et al. (1982), the general form for the backscatter coefficient from a surface is:

Vsct=Aillh=πhθ(Rϕ0.52)2
(33)

This notation presupposes that the illuminated areas on the slab top and on the ground (bottom of the slab) are the same. This assumption applies for small angle divergences of the radar beam, as is valid for strong directed antennae, low refraction on the air-snow interface and a ground surface that is flat and parallel to the surface of the snow slab. Accordingly, the volume component of the backscatter coefficient is:

σt=σas+σs+σg+σgv
(34)

In the incoherent approach, σt=σs+σg ; i.e., σ s represents the total radar cross-section (RCS) of N scatterers contributing backscatter from a volume σsσg of a snow slab (see subsection 2.3.1). Due to attenuation of the electromagnetic wave upon propagation within a snow mass, the total RCS of the snow slab is equal to:

where the attenuation coefficient, inherent to the i-th particle with a distance of ξ i from the coordinate origin O, is equal to:

where k e is the extinction coefficient. The factor α takes into account the two-way distance of forth and back wave propagation. The assessment of the scattering volume described above relates only to the absolute value of the volume. In the case of wave directed propagation and, consequently, directed attenuation, the summation should be performed along the propagation axis, ξ. Because N>>1, the summation in (3.5) can be replaced by integration. As was mentioned above, this integration should be carried out along the direction of wave propagation, i.e., along the ξ axes. Small changes in the propagation direction at the air/snow interface are ignored due to minor differences in the corresponding refraction coefficients. Due to the random spatial distribution of scatterers within the scattering volume, this discussion considers only the mean backscatter characteristics averaged over several illuminated areas, with particles having independent spatial locations and RCSs.

σs=i=1Nσi
(37)

where σ(ξ) is the running RCS (by the unit of a distance along axis ξ, dimension is Vsct=Aillh ) of the scattering volume. Assuming that the attenuation properties of snow remain the same along the propagation path σs=i=1Nσiαi , one can write:

αi=exp(2keξi)
(38)

If the scatterers’ RCSs are independent of their locations within the scattering volume (condition of incoherent approach), we get:

σs=i=1Nσiαi=0hθσ(ξ)α(ξ)dξ=0hθσ(ξ)α(ξ)dξ
(39)

where <σ> is the mean running RCS of the scattering volume.

Based on (3.4) and taking into account (3.9), the mean backscatter coefficient can be written as:

Taking into account that (kef(ξ)) is the radial size of the scattering volume (along axis ξ ), α(ξ)=exp(2keξ) is the mean total RCS of snow (ice) particles within the scattering volume (an oblique cylinder) while ignoring the attenuation and σs=σ0hθα(ξ)dξ is the corresponding mean volume specific backscatter coefficient, we can transform the ratio σs0=σsAill=σAill0hθα(ξ)dξ to the form:

Therefore, substituting (3.11) into (3.10), we can state:

σs=σhθ
(42)

Taking into account (3.8) and assuming the homogeneity of snow slab σv=σs*V , we can conduct the integration of (3.12) and finally arrive at:

This expression is the mathematical formulation of the semi-empirical model and has been obtained by Attema & Ulaby (1978) and Ulaby et al. (1982). In this paper, this model will be cited as the “A-U model.” For this case, the backscatter coefficient (3.13) depends only on the specific volume backscatter of the snow medium and not on the pulse volume sizes. This remarkable feature is due to (1) “overcomplete” burial of the pulse length into the snow medium, and (2) the backscattering normalization factor A ill (see 3.4) is the basis of the probing pulse.

The obvious imperfection of the model, as applied to sounding of thick snow, is the assumption of a constant extinction coefficient within the snow slab. The problem can be solved by designing an appropriate stratification model for the selected study area and modeling the spatial distribution of the extinction coefficient (e.g., Drinkwater et al., 2001).

The next limitation of the model is due to wave sphericity. Under depth-limited mode conditions (see section 2), the illuminated area changes, and the running RCS in (3.7) can not be assumed to be statistically homogeneous within the scattering volume. The same is true regarding the extinction coefficient in (3.8) as well. Thus, for this case, the A-U model in form (3.13) should be used with care. For the beamwidth-limited mode, the form in (3.13) can be used, taking in mind the “determined part” of the scattering volume (see comments to 2.15a).

The ground component of the backscatter coefficient is localized by the bottom location and can thus be expressed by its backscatter coefficient, taking into account the attenuation:

σAill=σAillhθhθ=σhθAillhcosθ=σsVcosθ=σvcosθ
(44)

where σs0=σvcosθ0hθα(ξ)dξ is the backscatter coefficient of the ground bottom surface governed only by the surface properties and its orientation in regards to the incidence of radar illumination..

3.2. Effective snow depth

Given the expression

(kef(ξ,θ,ψ))
(45)

(3.13) can be presented in the simple form:

σs0=σv12ke[1exp(2kehθ)]cosθ
(46)

where H eff is the effective depth of snow sounding (EDS); i.e., the depth such that the backscatter from which occurs as if without attenuation, with backscatter equal to that which would occur from a slab of larger real depth with attenuation due to absorption and scattering.

Parameter σg0=σg,max0exp(2kehθ) is the optical thickness of snow (SOT) along the wave propagation direction. It is useful to express the backscatter coefficient through the dimensionless EDS using the geometrical snow depth, h: σg,max0 . Given that expression, equation (3.16) can be written in the form:

12Lpcosθ[1exp(2hLpcosθ)]=Heff
(47)

where

σs0=σvHeff
(48)

is the maximal value of the backscatter coefficient for a given snow depth, h, ignoring both energy losses due to attenuation and scattering and the angular dependence of the scattering volume. The normalized effective depth of snow sounding (nEDS≡ kehθ=hθLp=τθ ) plays the role of a correction factor, and is equal to:

And the main formula of the model (3.2) can be written as:

σs0=σvhHeff0=σs0maxHeff0
(50)

It is useful to find the dependence of the correction factor on the SOT magnitude and the incidence angle. Let us consider two extreme cases:

a) “shallow” snow: σs0max=σvh

Given this relation,

Heff0=12τθ[1exp(2τθ)]
(52)

meaning that the EDS is approximately equal to the geometrical snow depth. Nevertheless, the small magnitude of the volume component compared with the backscatter from the bottom layer ( σs0=σs0maxHeff0+σg0 ) makes this dependence difficult to use in practice. For example, a snow layer with h<~1 m is practically transparent to electromagnetic irradiance of the C-band.

b) “thick” snow: τθ1

Heff012τθ[1(12τθ)]=1
(53)

Given this relation, and in accordance with (3.18) and (3.20), one gets:

σt0=σvh+σg0
(54)

where σvhσg0 . Equation (3.24) shows that the backscatter depends only on the penetration path (depth) under a constant incidence angle. In the “thick” snow regime, no additional snow accumulation contributes to the total backscatter due to the saturation effect. On the other hand, in this regime, the backscatter coefficient demonstrates an angular dependence, as the angle of incidence affects the propagation depth ( τθ1 ). The sensitivity of the backscatter coefficient to the changes in snow thickness takes place in the so-called “intermediate” regime, when Heff012τθ=Lpcosθ2h=12Dph Field data gives 20-30 meters of the penetration depth for the C-band (e.g., Hoen & Zebker, 2000), and yield values that can be considered the “working” range for probable snow measurements in this wave band. Values of the nEDS (3.19) can be assessed from the plot of the nEDS as a function of SOT, as shown in Figure 11.

media/image140.jpg

Figure 11.

Plot of the normalized effective depth of snow sounding versus the snow optical thickness

To evaluate the angular dependence of the nEDS, we consider the normalized Snow Depth (nSD), which is the depth of snow normalized to the penetration path:

σt0=σs0max+σg0
(55)

Comparing nSD with the definition of the SOT yields:

σs0max=12σvDp(θ)
(56)

Substituting (3.26) into (3.19), we have

A plot of (3.27) reduced to zero dB at θ=0 for different nSD is shown in Figure 12.. As follows from this plot, the angular dependence of the volume component of the backscatter coefficient in a practically appropriate range of angles, 0-500 is rather weak and equal to ~2dB for the thick snow regime. For the intermediate regime, with nSD~1, its range of variation is about 1 dB. Because this range is of the same order as the errors, it is difficult to expect a notable angular dependence in practical measurements inherent to the intermediate regime.

media/image144.jpg

Figure 12.

The angular dependence of the normalized effective depth of snow sounding for the “thick” (nSD>>1) and “intermediate” regimes

On the other hand, a notable angular dependence of the backscatter coefficient indicates that snow depth is greater than the penetration path. The intermediate regime represents practical interest for snow measurements. Suppose, for instance, the linear dependence of the geometrical snow depth on a horizontal distance S over the flat surface:

where k is the proportionality factor. Given that relation and taking into account (3.18) - (3.20), the normalized backscatter coefficient as a function of a distance is:

To specify the proportionality factor, let us assume that the geometrical snow depth reaches the propagation depth value at a distance, for example, equal to S 0 , i.e., nSD=τθcosθ . Thus,

Heff0=12cosθnSD[1exp(2nSDcosθ)]
(60)
media/image149.jpg

Figure 13.

Illustration of the expected relative spatial behavior of the backscatter coefficient under a linear dependence of the snow depth on a distance

where h(s)=kS is the normalized distance. For the intermediate regime, we can ignore the angular dependence of the backscatter coefficient. Let us also assume that the specific volume scattering is approximately the same for all points along that distance. The corresponding plot of expression (3.27), additionally normalized by a factor of σs0(S)σs0max=12cosθ[1exp(2kSDp(S,θ))] , is depicted in Figure 13.:

Thus, in accordance with the semi-empirical model, the expected profile of the backscatter coefficient should have an exponential pattern for a transect taken over terrain with a linear increase in snow depth within the intermediate regime condition.

4. Comparison the A-U model estimates and the experimental radar data on the Greenland ice sheet

For accurate comparison of the model with any experimental data, the study area must satisfy the “thick” snow condition to exclude the impact of the ground. Second, the snow pack must contain primarily Rayleigh ice particles with sizes less than approximately one-tenth of the wavelength. The most appropriate area matching this condition is the dry snow zone of the Greenland ice sheet. This region is not impacted by melt-freezing processes due to permanent sub-zero temperatures year-round. Therefore, the auspicious conditions for creation the ice lenses, tubes and large size aggregates, which inherent, for instance, to the percolation zone, are absent here. Additionally, the snow (firn) is up to several tens of meters deep (e.g., Zwally & Brenner, 2001), corresponding to the intermediate or even the thick snow regime. Due to known difficulties regarding in situ measurements in Greenland, there are no examples of the comparison of backscatter with snow depth in the literature for the intermediate or thick snow regimes.

media/image153.jpg

Figure 14.

Transect across hill (straight line crossing the dark spot on the SAR signature) in the north-eastern region of Greenland (left), and the corresponding backscatter (solid line, convex downwards curve) and elevation (dotted line, convex upwards curve) profiles (right)

Nevertheless, a certain qualitative example of the sensitivity of backscatter to the snow depth relative to the intermediate regime can be provided based on a topographic feature (a hill) in the northeast region of Greenland between 810 and 820 N, as shown in Figure 14..

The counter-phase behavior of the backscatter coefficient and the elevation profiles plotted above can be explained supposing that the snow depth is minimal on the hilltop and increases proportional to the distance from the top of the hill down to its base. This assumption closely matches the field and modeled data provided by Jaedicke et al. (2000).

For quantitative assessment of the semi-empirical model more regular terrain should be chosen. A typical dry snow area is presented on the RadarSAT-1 SAR image of eastern Greenland in Figure 15..

As is clearly seen, the magnitude of the backscatter coefficient is about -20 dB in the dry snow zone. This assessment coincides with data in the known literature (e.g., Drinkwater et al., 2001, Partington, 1998; Baumgartner et al., 1999, Forster et al., 1999) which also gives values of -12…-20 dB. We can estimate the corresponding value using the A-U model. We assume that the thick snow regime is valid, and that the snow slab consists primary of particles of roughly the same sizes; i.e., it is a monodisperse medium.

For these conditions (see 3.24), ignoring backscatter from the bottom surface, we get:

media/image154.jpg

Figure 15.

Illustration of the backscatter coefficient values (left) range for an arbitrary transect (straight line on RadarSAT-1 SAR signature, right) from the percolation to the dry snow areas within the east margin of the Greenland Ice Sheet (directions are from left to right for the plot of the backscatter profile and contra versa for the SAR signature; date of data acquisition: November 2, 2007; data granule ID: R1_62609_SWB_271)

σs0(S)σs0max=12cosθ[1exp(2S)]
(62)

Further calculations below aim to discover the specific volume backscatter, σ v and the extinction coefficient, k e . The estimates listed below will be conducted with the Rayleigh approach. In accordance with Ulaby et al. (1986), this approach bounds the upper electric size of particle, S=SS0 , with the inequality 2cosθ , where m is the complex refractive index of the particle matter and r u is the corresponding upper size of a particle. For ice particles with 2cosθσs0(S)σs0max=1exp(2S) observed with C-band radar (λ=5.6 cm), the maximal size satisfying the Rayleigh approach is 2.6 mm. Following known relationships (e.g., Ulaby et al, 1986), the specific volume scattering is:

σs0=σvLpcosθ2=12σvcosθke
(63)

where

In this formulae

ζu|m˙|0.5
(65)

where |m˙|1.7776 is the refractive index of the particle matter and n and η are the real and imaginary parts of the refractive index, respectively. The value of these components can be derived from the value of the complex dielectric permittivity of a material σv=Zs1λ4 , where ε’ and ε” are the real and imaginary parts of the dielectric permittivity.

Zs=26Vπ5|K|2i=1Nri6
(66)
|K|2=|m˙21m˙2+2|2
(67)

For monodisperse media ( r i =r ) equation (3.3) has the form:

where n 0 is the particle concentration.

The extinction coefficient can be calculated as a sum of the absorption (k a ) and scattering (k s ) coefficients (Ulaby et al., 1986):

ε˙=ε+jε
(69)

where Q a and Q s are the absorption and scattering cross-sections respectively. These parameters, under the Rayleigh approach, can be calculated by the following formulas:

n=[0.5((ε)2+(ε)2+ε)]0.5
(70)
and
η=[0.5((ε)2+(ε)2ε)]0.5
(71)

In accordance with (Ulaby et al., 1982)

Zs=26π5|K|2n0r6
(72)

The real part of the dielectric permittivity for ice, contained in the formulas above, is equal to 3.15 and in practice does not depend on temperature or wavelength. The imaginary part of the dielectric permittivity for ice is not constant with changes in temperature or the illumination frequency. In accordance with Matzler (1987) it can be expressed through the following empirical formula:

ke=ka+ks=n0[Qa+Qs]
(73)

where F is frequency in GHz and A, B and C are the empirical coefficients. For a temperature of -150C, which is more appropriate for the dry snow case, the coefficients are equal to: A=3.5*10-4, B=3.6*10-5, C=1.2. Now, the expressions obtained for the backscatter coefficient calculation may be combined. Substituting (4.2) in (4.1), we get:

Qa=λ2πζ3Im{K}
(74)

Taking into account (4.7) and (4.8)-(4.10), one can get:

Qs=2λ23πζ6|K|2
(75)

Simplifying this expression, we have:

Im{K}3ε(ε+2)2
(76)

An important feature of this expression is its independence with regards to the particle concentration, a reflection of the property inherent to the saturation regime mentioned above in section III. Since the scattering volume is less than the spatial duration of the probing pulse an increasing particle number is equivalent to an increasing particle concentration. Thus, no additional amount of the snow over the snow pack with a thickness greater than the penetration depth can cause a notable increase in the backscatter coefficient. This circumstance is also very useful from a simulation point of view, since there is no microstructure parameter that needs to be assumed. For the numerical calculation listed below we assumed that the snow within the dry snow zone consists of snow grains that are actually grains of ice that look like grains of rice, with a typical size of ~1 mm. Taking into account the fact for ice: ε’=3.15 and ε’>> ε” (Tiuri et al., 1984, Matzler, 1987), ε=AF+BFC . Given that information, equation (4.15) for an intermediate incidence angle of ~400 (cosθ=0.77) is equal to:

σs0=(12Zsλ4ke)cosθ
(77)

Taking into account (4.11) and (4.12), we get:

σs0=12λ426π5|K|2n0r6n0[λ2πζ3Im{K}+2λ23πζ6|K|2]cosθ
(78)

The plot of (4.17) is depicted in Figure 16. for C-band radar with F=5.3 GHz (λ≈ 5.6 cm).

media/image178.jpg

Figure 16.

The dependence of the backscatter coefficient on particle size in accordance with the A-U model for C-band radar and snow depth greater than the penetration depth (“thick” snow regime)

It follows from this plot that the observable backscatter coefficient values of -12…-20 dB can be caused by particles with radii of less than 0.35-0.17 mm. The current estimates are slightly lower compared with the same made by Partington (1998), who found that the backscatter coefficient would range from -20 dB for a mono-distribution of grain sizes with a mean radius of 0.25 mm to -2 dB for a mean grain radius of 1 mm. Our results showed a 0.17 mm mean radius for -20 dB and a 1 mm mean radius for -3 dB. The small discrepancy is probably due to the temperature dependence of the imaginary part of the ice’s dielectric permittivity. The grain sizes values assessed do not coincide with the field data. In accordance with Table 2. the mean particle radius within a ~1 m depth surface layer of dry snow is 0.5…2.0 mm. As seen in Figure 16.., these values should produce a backscatter coefficient of approximately from -8 dB to 0 dB. Thus, there is a discrepancy of ~12 dB between the observable data and the values forecast by the semi-empirical model. Therefore, the A-U model overestimates the backscatter coefficient significantly. It should be noted, in addition, that the particle radius can only increases with the snow depth due to the depth-dependent grain radius model (Alley et al., 1982):

σs0=12cosθ1ζ3|K|2Im{K}+23
(79)

where r 0 is the mean radius at the surface, C is the crystal growth rate and H a is the mean annual layer thickness. Given this model, various simulation experiments (e.g., Forster et al., 1999, Drinkwater et al., 2001) and field measurements (e.g., Jezek et al., 1994, Woods, 1994; Lytle & Jezek, 1994), the mean particle radius for the entire scattering volume is even greater than the surface values and the aforementioned discrepancy becomes even more firmly grounded.

#Mean ice particle radius, mmSourceLocation
10.5-2.0Benson, 1996Greenland
21.0-1.5Schytt, 1964Spitsbergen
3Less than 0.75 mmHigham and Craven, 1997Antarctic
40.2-0.6Woods, 1994Greenland
50.1-0.7Lytle and Jezek, 1994Greenland
Table 4.1. Literature data on the mean ice particle size near the surface of dry snow cover

Table 2.

Thus, although the A-U model gives a qualitative assessment of the main features of backscatter behavior from a snow slab, the result of quantitative comparisons with field data does not match the theoretical predictions.

5. Enhanced semi-empirical model of the volume component of the backscatter coefficient

The problem highlighted in the previous section can be resolved by considering the statistical properties of the small-scale fluctuations in the particle concentration and its scattering properties within the scattering volume. The appropriate method for doing so that is the so-called “slice” approach, which was primarily suggested in weather radar meteorology (Marshal & Hitchfeld, 1953; Smith, 1964). This approach exploits the known radar feature in accordance with which the particles of a “cloud” located close to the front of the incident radar wave are considered to be approximately at the same distance from the radar and reflect incident electromagnetic wave almost coherently. One can consider that these particles are embedded in a fictitious thin cylindrical volume (“slice”), whose base coincides with the surface of spherical wave front and side-bounded by the main lobe of the antenna pattern. Thus, the scattering volume into a snow slab can be represented as an adjoining series of these slices, as illustrated in Figure 17.. Each slice is much narrower than the radar wavelength in the wave propagation direction (Δs<<λ).

media/image180.jpg

Figure 17.

Simplified slice approach scheme. The arbitrary position of a single slice is shown on the cross section of the scattering volume

A slice’s radial size, Δs can actually be considered as a minimal spatial scale of the backscatter property fluctuation, if this scale is much smaller than the wavelength. Using this approach, it was shown by the author (Yurchak, 2009), that the specific volume component of the backscatter from any spatial extended geophysical target (SEGT), included the snow as well, can be presented as the incoherent sum of the radar cross-sections of individual scatterers only if the number of particles in the slices (n) are distributed in accordance with the Poisson law; i.e., the variance of the number of particles (Var(n) ) is equal to the mean number ( |K|2=0.16 ). Otherwise, the “classical” specific volume component should be corrected by a so-called deviation factor Y(ξ a , χ):

σs0=0.3856.25ζ3Im{K}+23
(80)

where

σs0=0.57751.06ζ3(AF+BFC)+1
(81)
and r2(h)=r02+ChHa is the Poisson index, n is the variation coefficient of the particle radar equivalent length (PREL≡ σv=σvclassY(ξa,χ) , σ p is the radar cross-section of an individual particle). Formulas (5.1) and (5.2) reflect the fundamental physical principal that the fluctuations of the medium parameters (inhomogeneities) are the cause of the electromagnetic wave scattering (e.g., Atlas, 1964, Fabelinskii, 1968). The parameter ξ a might be expressed through the measurable parameters of the snow particle size distribution function (PSDF), such as the variation coefficient of particle radius, Y(ξa,χ)=ξa2+χξa2+1 , and the skewness coefficient, Sk:
χ=Var(n)n
(82)

Following Shi et al. (1993), if the PSDF for snow can be presented by the lognormal distribution, the parameters of relationship (5.3) are: ξa=Stdev(a)a , and apσp . Given these relations, ξr=Stdev(r)r . Thus, the PREL variation coefficient (5.3) depends, via ξ r only on the geometric standard deviation, S r . The corresponding plot of the deviation factor (5.2) is depicted in Figure 18. for different values of the Poisson index χ.

media/image192.jpg

Figure 18.

Deviation factor versus particle radius variation coefficient for different deviations of Poisson index χ from the Poisson law

As can be concluded from analysis of the plot, the deviation of ~12 dB emphasized in the previous section can be explained by the homogeneous microstructure of the dry snow zone, with approximately identical particle sizes (ξ r <0.1) and negligible variations of the particle concentration (χ<<1) within the scattering volume. As an example of data supporting the current considerations, the results of comparative measurements of C-band backscatter from firn and in situ measurements of statistical parameters of firn dielectric permittivity for a study area in Antarctic (Zahnen et al., 2002) can be taken into account. In this study, dielectric permittivity was considered to be a measure of the heterogeneity of the snow/firn pack. A notable positive correlation was observed between backscattering signal strength and the standard deviation of the dielectric constant collected along 4 and 15 meters depth interval of several drills within a dry firn area. At the same time, a negative correlation takes place for a mean dielectric constant. Because the dielectric constant is approximately linearly proportional to the density of dry snow (e.g., Tiuri et al., 1984), and therefore linearly proportional to the number of ice particles per unit volume as well, these correlations match the behavior of the deviation factor (5.2) on the Poisson index. The index and, hence, scattering increase with increasing variations in particle number and decrease with increasing ξa=3ξr32ξr1+2Skξr3+3ξr1+Sk . Of course, more comprehensive field data are needed to confirm this hypothesis.

The considerations provided above are related to the simplest case of complete burial of the probing pulse into a snow slab in the beamwidth-limited mode. Obviously, that the result is still the same for incomplete burial of the radial size of the scattering volume. The transverse size change of a slice within the scattering volume in the depth-limited mode results in an additional gradient in particle number among slices, and the contribution of this additional gradient to the backscatter should be analyzed in the future.

6. Normalized snow depth

6.1. Derivation of the normalized snow depth from the backscatter coefficient based on the enhanced semi-empirical model

After enhancing the semi-empirical model by the “slice” approach, one resume assessment of the snow depth, the key parameter in all mass-balance studies. To derive this parameter, let us to write down the enhanced semi-empirical model (ESEM) in a form more fitting for the Greenland ice sheet:

ξr=exp(Sr2)1
(83)

where Sk=exp(Sr2)1[2+exp(Sr2)] is the backscatter coefficient from the ice layer, L p is the penetration path (2.18), Y is the deviation factor (5.2), and Sk=ξr3+3ξr is the slant snow depth. Let us also denote: n as the saturation backscatter coefficient. Given these expressions, the normalized Snow Depth ( σt0=σs0+σi0=σv12Lpcosθ[1exp(2hθLp)]Y+σi,max0exp(2hθLp) , see (3.25)) as a function of distance S along snow cover with a variable snow depth, h, can be derived from (6.1) to take the form:

Here, hθ=hcosθ , σvLp2cosθ=σ0 and nSD=hLp to simplify the notation. If nSD(S)=h(S)Lp=12cosθ(S)ln[1σ(S)Yσ][1σiYσ] then

σ(S)σt0(S)
(85)

The equation above indicates that the estimated normalized snow depth depends on number and kind of scatterers (σ ) and their statistical characteristics (Y) at small (less then wavelength) scales. Thus, to extract the nSD from measured radar data, σ0σ , one needs to know the product of the deviation factor and the saturation coefficient. To find the physical snow depth, the penetration path must also be known. It should

media/image206.jpg

Figure 19.

Typical backscatter coefficient profile across a transect within the margin of the Greenland Ice Sheet (in the vicinity of the Swiss camp)

be underlined, additionally, that the nSD σi0σi is not quite the same as the snow optical thickness Yσσi . Obviously, nSD(S)=12cosθ(S)ln[1σ(S)Yσ] . nSD is a combination of the electromagnetic wave direction (L p ) and vertical geometrical (h) scaling factors of the snow depth. These two terms coincide only for sounding in the nadir direction.

From the point of view of practical measurements, detailed analysis recommends stopping the nSD estimate at a distance when σ(S) . It should also be noted also that the deviation factor (Y) might be less or more than unit (Figure 18.). Because the product (=hLp) is generally unknown, it is reasonable to establish its value as ~2 dB higher than the value of the backscatter profile plateau shown, for instance, in Figure 19. The backscatter profile has the main features inherent to the theoretical profile in the linear approach, as shown in Figure 13..

media/image212.jpg

Figure 20.

The normalized snow depth profile derived from the backscatter coefficient profile in Figure 19.

The profile has an exponential-like increase at the beginning of transect that transitions into a plateau-like, approximately constant level. Most probably, this level represents the saturation level caused by thick snow of depth greater than the propagation depth. The corresponding nSD is shown in Figure 20.

The nSD profile plotted in Figure 20. is obtained under the assumption that the deviation factor and the extinction coefficient are constant across the entire profile distance. Although the approach demonstrated here does not give the absolute value of the snow depth without certain assumptions, its advantage is that the assumptions are clearly delineated and put forth. In turn, if one has a calibration point (spot), the technique allows extension of the remote sensing methodology on vast remote areas with similar electromagnetic properties. In addition, the nSD is a relative indicator of the snow environment stability and redistribution.

The nSD in form (6.3) is derived from the ESEM, which in turn is based on consideration of a flat surface, plane wavefront and long probing pulse. Hence, the result obtained is most appropriate for SAR sounding of a thick snow slab.

6.2. An example of time transformation of nSD profiles within a selected Greenland marginal area

The normalized snow depth is sensitive to the geometrical snow depth and to the extinction coefficient. Assuming the snow microstructure to be homogeneous within a limited area, one can suppose that nSD depends only on the snow depth and thus apply nSD for snow mass-balance assessment. An illustrative example of an nSD profile change after a 10-year period is shown in Figure 21..

media/image213.jpg

Figure 21.

Study area #3 (box with center: 75.070N and 54.300W) within the Greenland ice sheet margin with several transects normal to the coastline (left). Normalized snow depth profiles (right) for transect 3a (lower line in the study area) derived from RadarSAT-1 data. Distance scale: 1 sample=50 m. Solid line is nSD profile on 11/16/1997, dotted line is nSD profile on 11/12/2007

The nSD profiles were derived from the RadarSAT-1 data archive from the Alaska Satellite Facility in Fairbanks (granules R1_10617_SWB_261 from 1997 and R1_62753_SWB_261 from 2007 with coordinates of the center scene 75.230N/54.290W and 75.220N/54.370W respectively, descending modes for both). As follows from the center locations of the data, the SAR shots were performed from practically the same point in the orbit and thus can be analyzed in parallel without additional corrections.

A comparative assessment of the nSD profiles indicates notable differences between the profiles. Generally, these changes can be caused snow mass redistribution within the 10-year period, variation of the extinction coefficient due to possible snow metamorphism, or both factors simultaneously. The detailed analysis of these factors and their “weights” in the deviations of these profiles is beyond the scope of this work but is planning for further consideration.

7. Model spectral dependence of the backscatter

7.1. Spectral dependence of the backscatter in accordance with the A-U model

The spectral (frequency) dependence of the backscatter coefficient within the framework of the semi-empirical model (A-U and ESEM) is the topic of the current section. The aim of this discussion is to find out how spectral features depend on the thick snow properties. Particularly, Davis & Moore (1993), namely assumed that a two-frequency radar system can provide quantitative estimates of snow physical properties.

The model consists of two components that are definitely frequency-dependent: the specific volume component (4.2) and the normalized effective depth of sounding (3.19). Transforming the wavelength (λ) into the frequency (=1Lphcosθ=hθLp) , the combination of these two components gives, for the “intermediate” snow regime:

The frequency-dependent core of (7.1) is

(σ(S)Yσ)dB2dB
(87)

where Yσ is the slant snow depth; (F=cλ) . The dimension of function (7.2) is [L-3].

7.2. Estimation of spectral differences based on field penetration depth data

For further analysis, it is convenient to deal with the dimensionless spectral function, which can be obtained by multiplying the spectral core (7.2) by the cube of the wavelength:

σs0={YZsc4F4h}{Lp(F)cosθ2h[1exp(2hLp(F)cosθ)]}
(88)

To evaluate (7.3) the spectral dependence of the penetration path (depth) should be estimated. One can state that the high frequency irradiance interacts primarily with the surface layer of snow with the relative small particles. And the low frequency irradiance penetrates deeper and, thus, is affected by particles which sizes are larger in accordance with the grain growth model (e.g., 4.18). For example, the grain radius changes from 0.25 to 0.6 mm within the depth of 10 m is in model of Drinkwater et al., (2001). In model of Forster et al., (1999) the range of the particles’ changes is even more (up to 1 mm) and additionally depends on the accumulation rate. Thus, to get more realistic spectral dependence, we should take into account the depth of penetration of the electromagnetic wave into the snow and the corresponding grain size change within this depth. Therefore, there is an obvious paradox: to calculate frequency dependence of the penetration depth we should know the latter a priory. Moreover, due to chain: “frequency-penetration depth-particle size” the mean particle radius impacted the radar irradiance turns out to be dependent on the irradiance frequency. In additional, the A-U model was derived assuming the constant extinction coefficient, as it was mentioned in section 3. The solution of these questions should be evaluated in further research. Taking in mind the restrictions listed above, the penetration depth is assumed to be dependent only on frequency in the following estimation. Since the significance of the particles’ size change compared with the fixed mean size in the extinction coefficient calculation is currently unclear, it is reasonable to take the field data available in the literature (Table 3.) for rough assessment the spectral dependence of the penetration depth.

#Frequency, GHzPenetration depth, mLiterature source
15
10
20
35
18
9
4.3
1.5
Hofer & Matzler, 1980
25
10
20
35
20
10
4.8
1.4
Matzler, 1987
35.2
6.6
10.3
10.7
18.0
37
21.7
17.8
8.1
10.4
3.2
0.85
Rott et al., 1993
4103.4±1.3Devis & Poznyak, 1993
513.53.5±1.5Davis, 1996
65.327±4Hoen & Zebker, 2000

Table 3.

Literature data on the experimental penetration path (depth)These field data can be approximated by a simple empirical power function

A(F,hθ)=21Zsσs0Ycosθ={(Fc)4}{Lp(F)[1exp(2hθLp(F))]}
(89)

with the coefficient of determination equals to 0.9. In (7.4) the penetration path is in meters and frequency is in GHz. It should be noted that field data, provided in the first three rows of Table 3., are obtained based on measurements within the surface layer of the snow (up to 3 m). It means that these data also do not take into account adequately the actual stratification of the thick dry snow cover of the Greenland ice sheet and, therefore, the further estimation should be considered only as a first approximation. A plot of (7.3) is shown in Figure 22. This plot demonstrates a notable discrimination in behavior of the spectral function B(F,h θ ) for different snow depths. Particularly, the magnitudes of the curves at a frequency of 5 GHz are notably different. Differences of the backscatter magnitude for frequencies of 5 and 14 GHz as a function of the slant snow depth are depicted in Figure 23.

media/image221.jpg

Figure 22.

Dimensionless spectral characteristics of the backscatter for several values of the slant snow depth based on the empirical spectral dependence of the penetration depth

media/image222.jpg

Figure 23.

Slant snow depth as a function of the magnitude difference between the dimensionless spectral components of the backscatter coefficients at 5 and 14 GHz

As follows from the approximation results,

media/image224.jpg

Figure 24.

Spectral pattern of the experimental data of (Baumgartner et al., 1999) for an incidence angle of 400 for two test sites in Greenland

where the slant snow depth is in meters and the spectral difference Lp(F)=ke(F)1 is in dB. An illustrative example of the spectral difference observed in Greenland can be provided from the known literature. Baumgartner at al. (1999) conducted field measurements of the backscatter coefficient from the dry snow site GITS (77.10N, 610W) and from the site NASA-U (730N, 50.50W) located between the dry snow and percolation zone for four wavelengths, 5.3, 10, 13.5 and 17 GHz, at temperature from -210C to -150C, with the ground-based radar. The primary goal of this experiment was to study the angular dependence of the backscatter. For our purposes, the spectral dependence related to the incidence angle of 400 (where the backscatter is expected to be due mainly to the volume scattering) can be extracted for the four points and is shown in Figure 24. Taking into account that B(F,hθ)=λ3A(F,hθ)={Fc}{Lp(F)[1exp(2hθLp(F))]} , where μ is an unknown coefficient that we assume to be the same for all wavelengths, one can obtain the spectral ratio or the spectral difference in the logarithmic units for a pair of frequencies F1 and F2:

Lp223.44F1.44
(91)

If, for example, F 1 =5 GHz and F 2 =14 GHz, we get:

hθ=12.14exp(0.6ΔB514)
(92)

Differences in the backscatter coefficients of these frequencies for both test sites are approximately the same and equal roughly to -13 dB, as shown in Figure 24.. Assuming the backscatter error of ~±0.5 dB, the corresponding spectral difference (from 7.7) is approximately equal to 0.4±0.5 dB. In accordance with (7.5), the slant snow depth is ~6.7-20.8 m. Since these data are related to an incidence angle of 400, the estimate of actual snow depth is 5.1-15.9 m. This value is in accordance with a note in the paper that the depth of snow was ~7-8 meters for the NASA-U site, at least.

Thus, a rough assessment of the snow depth within several selected dry snow areas, based on the spectral dependence of the backscatter in the frames of the SEM (ESEM), gives results comparable with in-situ measurements.

7.3. Spectral characteristics of the penetration depth for model monodisperse cold snow

It is interesting to compare the empirical (7.4) and theoretical frequency dependence of the penetration path (depth). Calculation of the extinction coefficient has been completed in the same manner as in section 3, with the additional estimation of particle concentration n 0 (see (4.8)):

where B(F,hθ)=μσs0F3 is the fractional volume of ice in snow (0.2-0.4), ρ s and ρ i are density of snow and ice, respectively. We also should demarcate the frequency and particle size inherent to the Rayleigh approach ( ΔBdB,F1F2=B(F1)dBB(F2)dB=[σs0(F1)dBσs0(F2)dB]30[lgF1lgF2] ). The penetration depth (path) in snow is calculated for several combinations of snow density and particle size for a temperature of -150C. The result of this calculation is shown in Figure 25. together with the empirical dependence (7.4). The gently sloping of the empirical curve compared with the theoretical ones indirectly support the idea, mentioned above, about the contribution of the particle size in spectral behavior of the penetration depth (see the corresponded comments in the previous subsection).

In Figure 7.5, the results obtained from these calculations for T= -150C and the snow density of 0.3 g/m3 are compared with the known values (Ulaby et al., 1986, Fig. 19.64, p.1608) for T= -10C and density of 0.24 g/m3. The plot indicates the increasing of the penetration depth at low temperatures for the low frequency region, a trend most obvious for small particles of 0.5 mm radius. This behavior is explainable because the losses in snow caused by absorption decrease as the temperature decreases (Ulaby et al., 1986, Fig.Е.3, p. 2027), and the absorption mechanism prevails over scattering in the low frequency region. Pointed changes become more apparent for small particles due to more significant role absorption plays for them as compared with large particles.

media/image232.jpg

Figure 25.

Theoretical spectral dependence under the Rayleigh approach of the penetration depth for a temperature of -150C (our calculations) with comparison of the empirical relationship (7.4)

media/image233.jpg

Figure 26.

Comparison of the spectral dependence of the penetration depth calculated with Mie formulas for a temperature of -10С (Ulaby et al., 1986) with our calculations under the Rayleigh approach for a temperature of -150С

8. Backscatter depth profile model for snow sounding with broadband pulses of short effective (compressed) length

8.1. Phenomenological approach in the case of very short probing pulse

The above considerations apply when the scattering volume was less then the probing pulse volume. In the case of altimeter sounding of thick snow cover with short broadband chirp pulses this condition is not satisfied, the backscatter coefficient is represented not by the point of a radar signature but in the form of a backscatter profile along the wave propagation path. When ΔBdB,514=[σs0(5)dBσs0(14)dB]+13.41 and n0=fv(43πr3)1 , the volume component of the total backscatter from a snow area with the center located at a slant distance of R beneath the surface should simply mirror the backscatter coefficient depth distribution taking attenuation into account:

where 2πrFc|m˙|0.5 is the total radar cross-section (RCS) of N scatterers contributing backscatter from a scattering volume V sct from a distance R. Since the scattering volume is hsct(p)hθ , and taking the correction on the deviation factor into account (see section 5), the equation (8.1) can be presented in the form:

where σs0(R)=σs(R)Aillexp(2keR) is the specific volume backscatter, and the return signal is proportional to the cosine of the incidence angle and has the envelope of the exponential form distorted by the variability of the statistical, scattering and absorption properties of the snow medium along the wave propagation path. The illuminated area A ill (transverse size of the scattering volume) changes during the burial of the probing pulse into the snow slab due (mainly) to the wave sphericity in the depth-limited mode (see section 2). However, this effect does not impact the accuracy of the specific volume backscatter coefficient presented here, due to its the mutual cancellation in the nominator (via V sct) and the denominator of equation (8.2). As a result of this cancellation, for very short probing pulses, the backscatter envelope depends only on changes in the specific volume backscatter with distance.

8.2. Analytical derivation for commensurate length of the scattering volume and snow depth

In general case, when σs(R)=i=1Nσi and Vsct=hsct(p)Aillcosθ , the interactions between a probing pulse with a normalized form U(R) and an extended target such as snow cover is described by the convolution integral (Moore & Williams, 1957):

σs0(R)=Y(R)σv(R)AillVsctexp(2keR)=Y(R)σv(R)hsct(p)cosθexp(2RLp(R))
(96)

where P r (t) is the return signal, P T (t) is the transmitted pulse profile and P s (t) is a term that includes the distribution range of the scattering facets, their scattering properties and the effects of the antenna gain. As applied our consideration of volume backscatter, one can write

If, for example, the probing pulse has a Gaussian form, which is usually an adequate approximation for short pulse systems (e.g., Brown, 1977), the weighting function of the scattering volume is hsct(p)~hcosθ , where Δ q is the spatial duration of the scattering volume weight function at level hsct(p)~Lp . For level q=0.5 (-3 dB), β=-4ln0.5=2.8. If, in addition, the microstructural, statistical and attenuation properties of the snow slab do not change along the propagation path within the snow pack, one can write:

Pr(t)=PT(t)Ps(t)
(98)

For short pulses, if their envelope can be represented by the delta-function, i.e., σs0(R)=cosθ0dxU(xR)Y(x)σv(x)exp(2xLp(x)) , equation (8.5) yields the corresponding result obtained earlier (8.3) for the very short pulses configuration. For the long pulses, where U(R)=exp(βR2Δq2) and q=exp(β4) , which yields to the standard form of the A-U model (3.13). For the general case, a closed-form solution of the equation (8.5) can be obtained.

Taking into account the tabulated integral (Abramowitz & Stegun, 1972):

σs0(R)=cosθYσv0dxU(xR)exp(2xLp)
(99)

where U(xR)~δ(xR) is the error function, one can get the integral of equation (8.5) for a Gaussian probing pulse:

If one denotes σs0~0Ldxexp(2xLp) , then (8.5) has the form:

0dxexp(ax2bx)=eb24aπ2aΦ(b+2ax2a)|0=12eb24aπa[1Φ(b2a)]
(101)

Finally, if one denotes: Φ(x)=2π0xexp(t2)dt (see 3.24), 0dxU(xR)exp(2xLp)=12Δqπβexp[1β(ΔqLp)22RDp][1Φ(1β(ΔqLpβΔqR))] and Δqβ=Δ , the closed-form solution of equation (8.6) using these dimensionless parameters can be written as:

σs0(R)=cosθYσvΔπ2exp[(ΔLp)22LpR][1Φ(ΔLp1ΔR)]
(102)

Distance normalization in (8.9) is chosen through the parameter S for the purpose of comparing the backscatter profile calculated here with the modeled “clear” backscatter profile derived for a very short pulse (8.2) and governed only by attenuation. The equation obtained is close to the volume component of the SV model of Davis & Moore (1993). The distinctions of the current approach are that it uses a slightly different description of the pulse-snow interaction, the end result is written in terms of the backscatter coefficient depth profile, and the state of homogeneity of the snow mass is taken into account through the deviation factor, Y. It is interesting, that equation (8.9) coincides with the backscatter model from the sea surface, obtained by Barrick (1972), exactly by the form with, of course, other meaning of model parameters.

8.3. Comparison of the results obtained with the phenomenological approach

8.3.1. Very short probing pulse

For the limit case when σmax=12σvLpcosθ , and taking into account that (ΔLp)=γ and S=RLp (Abramowitz & Stegun, 1972), the corresponding backscatter coefficient is:

σs0(S)=Yσmaxπγexp(γ22S)[1Φ(γSγ)]
(103)

This result confirms relationship (8.2) obtained above based on phenomenological considerations.

8.3.2. Long probing pulse

Taking into account the asymptotic value of the error function γ1 (Abramowitz & Stegun, 1972), and taking only the first term of the extension, we get the value of (8.9) for γ>>1:

Φ(x)=Φ(x)
(104)

For long pulses, the available distance variation is much less than the pulse length, i.e., limΦ(x)x=1 . Hence,

σs0(S)γ1=Yσmaxπγexp[2S][1Φ(1γS)]=Yσmaxπγexp[2S][1+Φ(1γS)]2Yσmaxπγexp(2S)=cosθYσvΔπexp(2S)
(105)

This expression also confirms the result of the A-U model (enhanced) under thick snow condition (3.24). Thus, an important positive feature of the model (8.9) is that its results coincide with those obtained with phenomenological approaches for extreme values of the scattering volume (much less and much more than the penetration path respectively). Nevertheless, it should be noted that ignoring the change in the illumination area during the burial of the probing pulse into snow media is, strictly speaking, valid only for a short probing pulse or for the beamwidth-limited mode. For long and intermediate-length probing pulses with wavefront sphericity in the commonly-existing depth-limited configuration, this effect should be taken into account because there are slices (see sections 2 and 5) of different transverse sizes within Φ(x)x11ex21πx(112x2...) .

8.4. Numerical calculation of the return signal for intermediate length of the probing pulse

A plot of the normalized backscatter coefficient σs0(S)γ1Yσmaxπγexp(γ22S)exp[γ2+2S(Sγ)2]1π(γSγ)1 (S max is the normalized distance where function (8.7) has the maximum value at a given γ) is shown in Figure 27.

media/image270.jpg

Figure 27.

Illustration of the reproducibility of the snow backscatter uniform profile taking into account only attenuation (bold dash line) and finite probing pulse length normalized to the penetration depth (parameter (Sγ)γ1=RΔ1 ).

Legend: solid line, γ=0.01; dot line, γ=0.035; dashed line, γ=0.06; dashed-dot line, γ=0.085; dashed-dot-dot line, γ=0.11; long dashed line, γ=0.135

This figure illustrates the contribution of the pulse length to the return signal envelope pattern. The shorter pulse clearly gives a better approximation of the real backscatter profile. Contemporary altimeters use a probing pulse with a physical duration of 20-50 microseconds and a corresponding effective duration (due to intrapulse frequency chirping) equal to 3.125 nanoseconds (e.g., Rosmorduc et al., 2006). In spatial units, these values are equivalent to σs0(S)γ1Yσmaxexp[(Sγ)2](1Sγ2)γ11Yσmax=12YσvDp ≈0.5 m. For L p ~10-50 m, the expected values of the parameter γ range from 0.05 to 0.01. An important feature of this plot is that the value of the maximum of the return signal for any maximum position (S max) and values of parameter γ are the same:

9. Summary

The semi-empirical model of the volumetric component of the backscatter coefficient, proposed by Attema & Ulaby (1978) (SEM or A-U model) is considered in more detail. The basic advantage of the SEM is its mathematical simplicity and clear physical sense but the SEM, derived using the incoherent approach, does not explain some remote sensing experimental data collected over the dry snow zone of the Greenland ice sheet. We discuss the main configurations of radar sounding of snow, including the relationships between the technical and positional parameters of radar on one hand and the geometrical and electromagnetic properties of snow cover on the other. Based on this consideration, the field of applicability of the SEM is figured out for the further assessments. Some qualitative assessments of the configuration of the scattering volume were carried out using plane and spherical surface approaches. More detailed study of the diffraction of space-bounded spherical electromagnetic wave on a large and strongly absorbing dielectric sphere with a rough surface is needed to evaluate this problem further. Detailed SEM estimates and experimental radar data over the dry snow zone of the Greenland ice sheet are provided for the negative temperatures inherent to the winter seasons. Data analysis shows a significant discrepancy between the expected and actually measured values of the backscatter coefficient. Possible reasons for the observed deviations of the SEM results are discussed based on the enhanced semi-empirical model (ESEM) and analysis of the dependence of the backscatter coefficient on the snow optical thickness and its spectral function. The ESEM is obtained on the basis of the so-called "slice" approach, taking into account partially-coherent backscattering due to finite wave front size and considering the contribution of the statistical characteristics of small-scale fluctuations of concentration and of sizes of the ice particles in a snow pack. The known semi-empirical model is only a special case of the ESEM. Using the ESEM, some observable radar signatures of the Greenland ice sheet, including the dry snow zone, can be interpreted in terms of the statistical characteristics of the snow microstructure. The model obtained retains the mathematical simplicity and physical clarity of the known A-U model and significantly expands its application for the interpretation of experimental data. Available field data, which confirm the consistency of the ESEM with observed values of the backscatter coefficient, are also provided. As a consequence of this analysis, the concept of the normalized snow depth (nSD) is formulated, and its practical application for monitoring snow mass displacement is described. The frequency-dependent core of the backscatter has been analyzed as well. Empirical and theoretical spectral characteristics of the penetration depth are also presented. The result allow for improvements in our understanding of opportunities for a two-frequency method of sounding thick snow covers and promote more accurate interpretations of the known experimental data. The extension of the ESEM for sounding thick snow cover with broadband pulses of short effective length is also discussed, resulting in an exact description of the return signal envelope as a function of the electromagnetic and microstructural characteristics of the snow cover. The validity of the model is supported by the coincidence of the model estimates with assessments obtained for several specific cases using phenomenological approach. In conclusion, this work demonstrates that the enhanced version of the semi-empirical model can be applied successfully in the analysis of radar backscatter from thick snow cover.

10. Acknowledgements

This work was supported by NASA’s Cryospheric Sciences Program. The author would like to thank Dr. W. Abdalati for his permanent attention to this work and useful discussions.

References

1 - M. Abramowitz, I. A. Stegun, Eds 1972 Hanbook of mathematical functions with formulas, graphs, and mathematical tables. New York: Wiley.
2 - R. B. Alley, J. F. Bolzan, I. M. Whillans, 1982 Polar firn densification and grain growth. Ann Glaciol., 3 7 11 .
3 - D. Atlas, 1964 Advances in radar meteorology. Advances in Geophysics, 10 317 479 .
4 - E. Attema, F. Ulaby, 1978 Vegetation modeled as water cloud. Radio Science, 13 357 364 .
5 - D. E. Barrick, 1972 Remote sensing of sea by radar. In Derr, V. ed. Remote sensing of the troposphere. Washington, DC, U.S. Government Printing Office, 12 1-12.46.
6 - L. J. Battan, 1959 Radar Meteorology, Univ. of Chicago Press, Chicago, Ill.
7 - F. Baumgarnter, K. Jezek, R. R. Forster, S. P. Gogineni, I. H. H. Zabel, 1999 Spectral and angular ground-based radar backscatter measurements of Greenland Snow facies. IGARSS’99 Proceedings, 28 June-2 July 1999, Humburg, 2 1053 1055 .
8 - C. S. Benson, 1996 Stratigraphic studies in the snow and firn of the Greenland ice sheet. U.S. Army Cold Reg. Res. and Eng. Lab (CRREL), Rep.70, 120 pp. Hanover, NH. (Revised edition of 1962 report).
9 - G. S. Brown, 1977 The average impulse response of a rough surface and its applications. IEEE Transactions on Antennas and Propagation, AP-25 , 1 January 1977, 67 74 .
10 - G. A. Clift, 1985 Use of radar in meteorology. Secretariat of the World Meteorological Organization, Geneva.
11 - C. H. Davis, 1996 Temporal change in the extinction coefficient of snow on the Greenland Ice sheet from an analysis of Seasat and Geosat altimeter data. IEEE Trans. Geosci. Remote Sens., 34 5 September 1996, 1066 1073 .
12 - C. H. Davis, V. I. Pozdnyak, 1993 The depth of penetration in Antarctic firn at 10 GHz. IEEE Trans. Geosci. Remote Sens., 31 5 September 1977, 1107 1111 .
13 - C. H. Davis, R. K. Moore, 1993 A combined surface- and volume scattering model for ice-sheet radar altimeter. Journal of Glaciology, 39 133 675 686 .
14 - M. D. Drinkwater, D. G. Long, A. W. Bingham, 2001 Greenland snow accumulation estimates from satellite radar scatterometer data. Journal of Geophysical Research, 106 D24 December 2001, 33 935-33,950.
15 - ERS-1 and 2 SAR Images. 1996 ASF document. September 18, 1996, 21 p. Available online: http://www.asf/alaska.edu/content/reference/pdf/ERSproducts.pdf
16 - D. M. Etheridge, L. P. Steele, R. J. Francey, R. L. Langenfelds, 2001 Historic CH4 records from Antarctic and Greenland ice cores, Antarctic firn data, and archived air samples from Cape Grim, Tasmania. Available online: http://ediac.ornl.gov/trends/atm_meth/lawdome_meth.html
17 - I. L. Fabelinskii, 1968 Molecular scattering of light. Translated from Russian, New York, Plenum Press.
18 - R. R. Forster, K. C. Jezek, J. Bolzan, F. Baumgartner, S. P. Gogineni, 1999 Relationship between radar backscatter and accumulation rates on the Greenland ice sheet. Int. J. Remote Sensing, 20 15 16, 3131 3147 .
19 - A.K. Fung, 1994 Microwave Scattering and Emission Models and Their Applications. Boston, Artech House, Inc.
20 - E. E. Gossard, R. G. Srauch, 1983 Radar observation of clear air and clouds. Elsevier, New York.
21 - M. Higham, M. Craven, 1997 Surface mass balance and snow surface properties from the Lambert Glacier Basin traverses 1990 94 . Tech. Rep. Res Rep. 9, Antarctic CRC, Hobard, Australia.
22 - E. W. Hoen, H. A. Zebker, 2000 Penetration depths inferred from interferometric volume decorrelation observed over the Greenland Ice Sheet. IEEE Trans. Geosci. Remote Sens., 38 6 November 2000, 2571 2583 .
23 - R. Hofer, C. Matzler, 1980 Investigations of snow parameters by radiometry in the 3 to 60 mm wavelength region. Journal of Geophysical Research, 85 453 460 .
24 - C. Jaedicke, T. K. Thiis, B. Bang, 2000 The snowdrift around a small hill in the High Arctic. Proceedings of the 4th Int. conference on Snow Engineering. Trondheim, Norway, 19-21 June 2000. Hjoth-Hansen, Holand, Loset & Norem (eds). Balkema, Rotterdam, 75 80 .
25 - K. C. Jezek, P. Gogineni, M. Shanableh, 1994 Radar measurements of melt zones on the Greenland Ice Sheet. Geophysical Research Letters, 21 1 33 36 .
26 - J. R. Kendra, K. Sarabandi, F. T. Ulaby, 1998 Radar measurements of snow: experiment and analysis. IEEE Trans. Geosci. Remote Sens., 36 3 May 1998, 864 879 .
27 - J. Koskinen, 2001 Snow monitoring using microwave radars. Thesis for the degree of Doctor of Technology. Helsinki University of Technology, Espoo, January 2001, Report 44, 31p.
28 - T. W. Liu, 2002 Progress in scatterometer application. Journal of Oceanography, 58 121 136 .
29 - Y. T. Lo, S. W. Lee, Eds, 1993 The antenna handbook. Volume II, Antenna theory. New York, Van Nostrand Reinhold.
30 - V. I. Lytle, K. Jezek, 1994 Dielectric permittivity and scattering measurements of Greenland firn at 26.5-40 GHz. IEEE Trans. Geosci. Remote Sens., 32 March 1994, 290 295 .
31 - C. Matzler, 1987 Application of Interaction of Microwaves with the natural snow cover. Remote Sensing Reviews, 2 259 387 .
32 - R. K. Moore, C. S. Williams, 1957 Radar terrain return at near vertical incidence. Proceedings Institute of Radio Engineers, 45 228 238 .
33 - K. Naito, D. Atlas, 1966 On microwave scatter by partially coherent clouds, In Proc. 12th Weather Radar Conf., AMS, Norman, OK, October 17-20, 1966, 7 12 .
34 - M. H. Newkirk, G. S. Brown, 1996 A waveform model for surface and volume scattering from ice and snow. IEEE Trans. Geosci. Remote Sens., 34 2 March 1996, 444 456 .
35 - M. H. Newkirk, G. S. Brown, 1992 Issues related to waveform computations for radar altimeter applications. IEEE Trans. Geosci. Remote Sens., 40 12 December 1992, 1478 1488 .
36 - C. R. Noveltis, 2005 Feasibility study of imaging the Antarctic ice using a spaceborne P-band radar. WP 200: Electromagnetic Model Review Document, ESA Contract 18195 04/NL/CB, 91 2005.
37 - RadarSAT-1 Standard Beam SAR Images. 1999 ASF document. June 21, 20 p. Available online: http://www.asf/alaska.edu/content/reference/pdf/R1stanprod.pdf
38 - V. Rosmorduc, J. Benveniste, O. Lauret, M. Milagro, N. Picot, 2006 Radar Altimetry Tutorial, J. Benveniste and N. Picot (Eds)., 301p. Available online: http//www.altimetry.info.
39 - H. Rott, K. Sturm, H. Miller, 1993 Active and passive microwave signatures of Antarctic firn by means of field measurements and satellite data. Annals of Glaciology, 17 337 343 .
40 - V. Schytt, 1964 Scientific results of the Swedish glaciological expedition to Nordaustlandet, Sptisbergen, 1957 and 1958. Geografiska Annaler, 44 243 281 .
41 - J. B. T. Scott, D. Mair, P. Nienow, V. Parry, E. Morris, 2006 A groud-based radar backscatter onvestigation in the percolation zone of the Greenland ice sheet. Remote Sensing of Environment, 104 361 373 .
42 - J. Shi, R. E. Davis, J. Dozier, 1993 Stereological determination of dry-snow parameters for discrete-scatterer microwave modeling. Ann. Glaciol., 17 295 299 .
43 - A. J. F. Siegert, H. Goldstein, 1951 Coherent and incoherent scattering from assemblies of scatterers, /In: Propagation of short radio waves. D.E. Kerr, (Ed), McGraw-Hill, New York.
44 - M. E. Tiuri, A. H. Sihvola, E. G. Nyfors, M. T. Hallikainen, 1984 The complex dielectric constant of snow at microwave frequencies. IEEE Journal of Oceanic Engineering, OE-9 , 5 377 382 , December 1984.
45 - L. Tsang, J. Pan, D. Liang, Z. Li, D. W. Cline, Y. Tan, 2007 Modeling active microwave remote sensing of snow using dense media radiative transfer (DRMT) theory with multiple-scattering effects. IEEE Trans. Geosci. Remote Sens.. 45 4 April 2007, 990 1004 .
46 - F. T. Ulaby, R. K. Moore, A. D. Fung, 1982 Microwave Remote Sensing, Active and Passive. Vol. II. Radar Remote Sensing and Surface Scattering and Emission Theory. Artech House, Inc., Norwood.
47 - F. T. Ulaby, R. K. Moore, A. D. Fung, 1986 Microwave Remote Sensing, Active and Passive. Vol. III. From Theory to Applications. Artech House, Inc., Norwood.
48 - F. T. Ulaby, P. Siqueira, A. Nashashibi, K. Sarabandi, 1996 Semi-empirical model for radar backscatter from snow at 35 and 95 GHz. IEEE Trans. Geosci. Remote Sens., 34 5 1059 1065 , September 1996.
49 - G. A. Woods, 1994 Grain growth behavior of the GISP2 ice core from central Greenland. Technical Report Series 94 002 , Earth System Science Center, Penn State University, USA.
50 - B. S. Yurchak, 2009 Radar volume backscatter from spatially extended geophysical targets in a “slice” approach. IEEE Trans. Geosci. Remote Sens., 47 10 (in press).
51 - N. Zahnen, F. Jung-Rothenhausler, H. Oerter, F. Wilhelms, H. Miller, 2002 Correlation between Antarctic dry snow properties and backscattering characteristics in Radarsat SAR imagery. Proceedings of EARSeL-LISSIG_Workshop Observing our Cryosphere from Space, Bern, March 11-13, 2002, 140 148 .
52 - H. J. Zwally, A. C. Brenner, 2001 Ice sheet dynamics and mass balance. In book: Satellite altimetry and earth sciences. L.L. Fu and A. Cazenave, (Eds). International Geophysics series, 69 351 369 .