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 21

Soil moisture estimation using L-band radiometry

By Alessandra Monerris and Thomas Schmugge
DOI: 10.5772/8334

Article top

Overview

a) SMOS field of view. Blue lines indicate the observation angle, black ellipses indicate the pixel size depending on the position within the snap-shot, and red lines show how a single pixel can be observed under different viewing positions during a satellite overpass. (b) MIRAS during the IVT campaign at ESA-ESTEC Maxwell anechoic chamber premises (picture courtesy of EADS-CASA Espacio)
Figure 1. a) SMOS field of view. Blue lines indicate the observation angle, black ellipses indicate the pixel size depending on the position within the snap-shot, and red lines show how a single pixel can be observed under different viewing positions during a satellite overpass. (b) MIRAS during the IVT campaign at ESA-ESTEC Maxwell anechoic chamber premises (picture courtesy of EADS-CASA Espacio)
Picture of the MOUSE 2004, T-REX 2006, and SMOS REFLEX 2003 experiment sites.
Figure 2. Picture of the MOUSE 2004, T-REX 2006, and SMOS REFLEX 2003 experiment sites.
Variation of the emissivity with the field capacity at H-pol. Colours indicate the incidence angle: 25  (blue), 35  (red), and 45  (green).
Figure 3. Variation of the emissivity with the field capacity at H-pol. Colours indicate the incidence angle: 25 (blue), 35 (red), and 45 (green).
Minimum root mean square error of the soil moisture retrieval. Results for a soil dielectric constant model and different soil moisture ground-truth data are represented.
Figure 4. Minimum root mean square error of the soil moisture retrieval. Results for a soil dielectric constant model and different soil moisture ground-truth data are represented.
Variation of the emissivity measured during the T-REX 2006 as a function of the incidence angle for the four ploughs. Colours indicate the polarisation: H-pol (red) and V-pol (blue). Line style and icons distinguish the ploughing.
Figure 5. Variation of the emissivity measured during the T-REX 2006 as a function of the incidence angle for the four ploughs. Colours indicate the polarisation: H-pol (red) and V-pol (blue). Line style and icons distinguish the ploughing.
Variation with time of the emissivity measured during the T-REX 2006 experiment as a function of the incidence angle. The mean value of LAURA measurements (black icons) and their standard deviation during the whole experiment (bars) have been represented together with IEM simulations.
Figure 6. Variation with time of the emissivity measured during the T-REX 2006 experiment as a function of the incidence angle. The mean value of LAURA measurements (black icons) and their standard deviation during the whole experiment (bars) have been represented together with IEM simulations.
Comparison between T-REX 2006 measurements (black icons) and the models in Wang & Choudhury (1981), Mätzler (2006, section 4.7), Schneeberger et al. (2004), Mo & Schmugge (1987), and Escorihuela et al. (2007). Solid lines indicate results using the standard deviation of the surface height profile while dashed lines indicates results using the λ-average standard deviation of height
Figure 7. Comparison between T-REX 2006 measurements (black icons) and the models in Wang & Choudhury (1981), Mätzler (2006, section 4.7), Schneeberger et al. (2004), Mo & Schmugge (1987), and Escorihuela et al. (2007). Solid lines indicate results using the standard deviation of the surface height profile while dashed lines indicates results using the λ-average standard deviation of height
Mean value and standard deviation of the emissivity and of the ground-truth soil moisture as a function of (a) the incidence angle and for all azimuth angles, and (b) the azimuth angle and for all incidence angles. DoY 185 (solid line, wet soil) and DoY 181 (dashed-dotted line, dry soil) have been represented. Blue and red indicate vertical and horizontal polarisations, respectively.
Figure 8. Mean value and standard deviation of the emissivity and of the ground-truth soil moisture as a function of (a) the incidence angle and for all azimuth angles, and (b) the azimuth angle and for all incidence angles. DoY 185 (solid line, wet soil) and DoY 181 (dashed-dotted line, dry soil) have been represented. Blue and red indicate vertical and horizontal polarisations, respectively.
Dependence of the retrieved vegetation albedo (a), and opacity (b) on the incidence angle. Icons represent the value estimated from LAURA measurements for each day of experiment and all azimuth angles. The mean value is represented with a black solid line. (c) Ground truth soil moisture vs. estimates from radiometric measurements.
Figure 9. Dependence of the retrieved vegetation albedo (a), and opacity (b) on the incidence angle. Icons represent the value estimated from LAURA measurements for each day of experiment and all azimuth angles. The mean value is represented with a black solid line. (c) Ground truth soil moisture vs. estimates from radiometric measurements.

Soil moisture estimation usingL-band radiometry

Alessandra Monerris1 and Thomas Schmugge2

1. Introduction

About 70% of the Earth is covered with water, 97% of which is part of oceans. Heating of oceans by the Sun keeps the Earth's water in a continuous circulation from the atmosphere to the Earth and back to the atmosphere through condensation, evapotranspiration, and precipitation processes. This continuous water motion is called water cycle. The groundwater in the unsaturated or vadose zone of land is usually referred to as soil moisture. The thickness of this zone extends from soil surface to a few metres below the surface in humid regions, and to 300 m or more below surface in arid regions. Although the water held by soils is a small fraction of the Earth's water budget, soil moisture plays an important role in the water cycle since it controls the proportion of rainfall that percolates, runs off, or evaporates from the land, influences plants growth and transpiration, and is related to the precipitation variability within a region (Koster et al., 2004).

Although soil moisture is one of the main parameters used in climate models, because of the large temporal and spatial variations of soil moisture sparse in-situ measurements are inadequate to be of much use in these models. Remote sensing with sufficient accuracy would provide meaningful soil moisture data over large regions. The maximum sensitivity to soil moisture is achieved in the lower range of microwave frequencies, which goes from 1 to 5 GHz (wavelength from 30 cm to 5 cm, respectively). The justification to this behaviour can be found if the dielectric properties of a target are analysed, since they have a large influence on its microwave brightness temperature. Land surfaces can be considered as a mixture of soil, water, and air particles and thus soil moisture might be estimated from measurements if the contrast between the water and soil particles were large enough. At microwave frequencies the real part of dry soil and water dielectric constants are approximately 4 and 80, respectively (Jackson & Schmugge, 1989. The scientific literature concluded that microwave radiometry at L-band (1.4-1.427 GHz) is optimal to estimate soil moisture, not only because it is very sensitive to soil moisture, but also because provides all-weather coverage, since the atmosphere at microwave frequencies may be considered nearly transparent, and vegetation is semi-transparent, which allows observations of the underlying layers (Eagleman, 1976, Wang & Choudhury, 1981, Jackson & Schmugge 1991). This specific band was chosen because it is a radio astronomy quiet band.

1.1. Satellite missions for soil moisture estimation using passive microwave sensors

While L-band microwave sensors provide maximum sensitivity to soil moisture their long wavelength implies the need for large antennae to achieve useful spatial resolution. For example, to obtain a ground resolution of 50 km or less using classical solutions on low-orbit satellites implies an antenna size of up to 20 m. At present, different scientific groups are developing new techniques to face this problem, and two space missions have been proposed to measure soil moisture at global scale: ESA's SMOS and NASA's SMAP.

The Soil Moisture and Ocean Salinity mission

The Soil Moisture and Ocean Salinity (SMOS) mission will observe soil moisture over land and salinity over oceans, and will measure snow and ice areas contributing to studies of the cryosphere (Kerr et al., 2001, Barré et al., 2008). Scheduled for launch in autumn 2009, it is the second Earth Explorer mission and part of ESA's Living Planet Programme. SMOS was thought of as a cost-effective, demonstrator mission with a nominal (extended) lifetime of 3 (5) years. The orbit is quasi-circular, sun-synchronous and dawn-dusk, and will be in the low-Earth range, at 763 km.

The SMOS mission provides a completely new approach in the field of remote sensing at L-band. Its payload MIRAS (Microwave Radiometer by Aperture Synthesis) is the first-ever two-dimensional interferometric radiometer in space. MIRAS consists of a Y-shaped antenna array with three arms, each arm having an approximate length of 4.5 m and 21 dual-polarisation L-band antennae or LICEF (Light Cost Effective Front-end) spaced 0.875 times the wavelength. Each LICEF is a total power radiometer in its own. Nine redundant antenna and three full-polarimetric noise injection radiometers are located in a fixed structure in the centre of the array, while the arms are divided in three segments to be folded during the launch (McMullan et al., 2008). If classical solutions would have been used, with this antenna size and orbit altitude the field of view (FOV) would be of near 3000 km in diameter. However, because of the microwave interferometry technique, the instrument Y-shape and the antenna spacing, the resulting FOV is a hexagonal-like area of less than 1000 km (see Figure 1). The spatial resolution varies from 35 km at the FOV centre to 50 km at the border. MIRAS will provide observations of a single azimuth at various incidence angles (from 0 to 55 ) and radiometric resolutions depending on its position within the field of view during a satellite overpass. This fact will much improve the retrieval algorithms because a lot of independent information of each pixel will be registered, and will permit the estimation of soil roughness, vegetation opacity and albedo, etc. together with soil moisture.

To satisfy the scientific requirements, SMOS aims at providing global maps of soil moisture every 3 days with a ground resolution better than 50 km, and a 0.04 m3/m3 volumetric humidity. This soil moisture accuracy is referred to pixels outside mountainous, urban, and partially frozen or snow-covered areas. For sea salinity, maps with accuracy better than 1.2 psu and 200 km ground resolution will be acquired every 30 days.

media/image1.jpg

Figure 1.

a) SMOS field of view. Blue lines indicate the observation angle, black ellipses indicate the pixel size depending on the position within the snap-shot, and red lines show how a single pixel can be observed under different viewing positions during a satellite overpass. (b) MIRAS during the IVT campaign at ESA-ESTEC Maxwell anechoic chamber premises (picture courtesy of EADS-CASA Espacio)

The Soil Moisture Active and Passive mission

The Hydrosphere State (HYDROS) mission was proposed by NASA to enhance the understanding of the land hydrosphere state and improve the climate prediction models (Entekhabi et al., 2004). The mission had a circular, polar, sun-synchronous, 6 am/pm equator crossing orbit located at 670 km altitude. Hydros aimed at providing soil moisture estimates with a 4% volumetric accuracy in the top 2-5 cm and capturing freeze/thaw state transitions in integrated vegetation-soil continuum at the spatial scale of landscape variability (3 km). The payload consisted of both active and passive sensors. An L-band radiometer measured the first, second, and third Stokes parameters with a 40 km spatial resolution, and 1 K relative accuracy, while an L-band radar acquired VV, HH, and HV polarisations with a 10 km resolution, and 0.5 dB accuracy for VV and HH. Observations at a constant incidence angle between 35 and 50 . Both the radar and the radiometer share a 6 m diameter reflector antenna which rotates about its nadir axis. The HYDROS mission was cancelled in 2005, but many HYDROS science and technology issues are being reviewed to be implemented in the future Soil Moisture Active and Passive (SMAP) mission.

1.2. Scope of this chapter

Space missions for Earth observation require a huge amount of previous work to validate the existing science and technology. This is accomplished through airborne and ground-based field experiments. This chapter reviews the basics of microwave radiometry and the existing models of soil emissivity and soil dielectric constant. The field experiments over land carried out during the last three decades are summarised, paying special attention to the recent experiments carried out by the Remote Sensing Laboratory from the Universitat Politècnica de Catalunya (UPC) as part of the preparatory activities for the SMOS mission over land.

2. L-band emission of land covers

Microwave remote sensing is based on the measurement of the thermal radiation or brightness temperature of a target, which is determined by its physical temperature and emissivity. The emissivity of land covers depends on soil moisture, but also on soil temperature (Choudhury et al., 1982, Wigneron et al., 2001), soil surface roughness (Mo & Schmugee, 1987, Escorihuela et al., 2007), vegetation canopy (Jackson & Schmugee, 1991, Ferrazzoli et al., 2002), snow cover (Schwank et al, 2004), relief (Mätzler & Standley, 2000, Talone et al., 2007), etc.

2.1. Soil dielectric constant model

The dielectric constant determines the response of the soil to an incident electromagnetic wave. This response is composed of two parts, real and imaginary (ε = ε' + jε'') which determine the wave velocity and energy losses respectively. In a non-homogeneous medium such as soils, the dielectric properties have a strong impact on its microwave emission. However, the relationship between the soil dielectric constant and the soil physical properties is not straightforward. A large number of studies have been performed during the last decades to find out this relationship since it plays an important part in the soil moisture retrieval algorithms from remote sensing data (Hipp, 1974, Wang & Schmugge, 1980, Topp et al., 1980, Hallikainen et al., 1985, Dobson et al., 1985, Roth et al., 1992, Mironov et al., 2004). Some of these models are simple empirical models in which data is fitted by a curve unique for all soils; others propose semi-empirical approaches which take into account some soil physical properties. The dielectric constant of dry soils is almost independent of temperature (Topp et al., 1980) and frequency. On the contrary, wet soils show a complex behaviour depending on the interaction between soil, water, and air particles. Hallikainen et al. (1985) performed a series of dielectric constant measurements of five soils with different texture composition at frequencies between 1.4 and 18 GHz and found out that texture has a strong effect on the dielectric behaviour which is specially pronounced at frequencies below 5 GHz.

In the dielectric-mixing model by Roth et al. (1992), differences in soil texture and bound water and free water are ignored altogether. Other models use a semi-empirical approach that contains a model of the complex dielectric constant and the volume fraction of each of the soil components. This kind of approach was followed by two widely used models presented in Wang & Schmugge (1980) and Dobson et al. (1985). The starting point of both of them is the dielectric mixing model by (Birchak et al., 1974),

εα = Vsεα + Vaεaα + Vεfwfwα + Vbwεfwα
(1)

where Vs (ε s), Va (ε a), Vfw (ε fw) and Vbw (ε bw) are the volume fraction (dielectric constant) of solid phase, air, free water, and bound water in the soil, respectively. The expression in (1) can be rewritten as a function of the bulk density ρ b, the particle density ρ s, and the volumetric moisture mv as

εα = 1 + ρbρs(εαs - 1) + Vεfwfwα + Vbwεfwα - mv
(2)

The Wang & Schmugge (1980) model was proposed for 1.4 and 5 GHz frequencies and starts from (2) with α=1. This model provides separate dielectric constant equations for volumetric water content lower than or greater than the transition moisture wt. The transition moisture is the moisture content at which the free water phase begins to dominate the soil hydraulics, and it is strongly dependent on texture. The soil dielectric constant is then estimated as

ε = {mvε x+ (P - mv)ε a+ (1 - P)εr,mvwtwtε x+ (mv - wt)ε fw+ (P - mv)εa+(1 - P)εr,mvwt
(3)
with
εx = {ε i+ (ε fwε)imvwtγ,mvwtε i+ (ε fwε)iγ,mvwt
(4)

and γ = -0.57wt + 0.481, ε i, ε a, ε fw, and εr the dielectric constants of ice, air, free water, and rock, respectively, ε x the dielectric constant of the initially adsorbed water, and P the soil porosity. The variable ε fw is estimated using the Debye equation.

On the other hand, the Dobson et al. (1985) model starts from (2) and assumes that there is not distinction between bound and free water. Taking this into account, the new expression for the dielectric constant is

ε = (1 + ρbρs(εαs - 1) + mvβεfwα - mv)1/α
(5)

The real and imaginary parts of the dielectric constant of soils are obtained separately using the percentage of sand S and clay C in the soil, and β =(127.48-0.519S-0.152C)/100 for the real part, and β =(1.33979-0.603S-0.166C)/100 for the imaginary part. For further information on dielectric constant models refer to Behari (2005) and Chukhlantsev (2006).

2.2. Emission from bare soils

The emissivity of bare soils es is given by

Γsp=1esp
(6)

where Γs,p is the total reflectivity at p polarisation in a chosen observation direction, which can be expressed in terms of the bistatic scattering coefficients σpp and σpq as follows (Fung, 1994):

Γsp(θi,ϕi)=14π4π[σpp(θi,ϕi,θr,ϕr)+σpq(θi,ϕi,θr,ϕr)] dΩr
(7)

Subscripts i and r stand for the incident and scattered radiation, and θ and ϕ are the incidence and azimuth angles, respectively. The coefficient σ pp , often called coherent component, considers the scattering on the observation direction and polarisation, while σ pq , or incoherent component, considers the scattering in whatever other direction and polarisation. The integration over the entire half space and errors in the modelling of the bistatic scattering coefficients make the computation of (7) unfeasible or prone to errors. For this reason, the uniform half space model with smooth surface is commonly adopted. This model is simpler and thus suitable to be used in soil moisture retrieval algorithms from remotely sensed data. In this case, the total reflectivity in (7) equals the reflection coefficient given by the Fresnel formulation.

Soil brightness temperature is related to soil emission through the soil effective temperature, Teff, TBp = e p Teff. The theoretical effective temperature of a soil profile can be estimated as (Ulaby et al., 1986):

Teff=0T(z)α(z)exp[0xα(z')dz']dz
(8)

where T(z) is the thermodynamic temperature at a depth z, and α is the attenuation coefficient, which is a function of the soil dielectric constant and of the microwave emission wavelength λ. Several simple formulations have been developed to estimate the soil effective temperature from soil properties, and soil moisture and temperature profiles. Choudhury et al. (1982) proposed a parameterisation of Teff based on the soil temperature at deep soil (T) corresponding to a depth between 50 cm and 1m, and on a surface temperature (Tsurf) corresponding to a depth of 0 to 5 cm:

Teff=T+C(TsurfT)
(9)

The coefficient C was considered constant for a given frequency, and equal to 0.246 at L-band. On the other hand, Chanzy et al. (1997) presented a model for the soil effective temperature at L- and C-bands based on the air temperature, a deep soil temperature, and the brightness temperature measured at X-band and V-pol. Wigneron et al. (2001) proposed a parameterisation based on (9), but with a coefficient C dependent on the volumetric water content ws, and two semi-empirical parameters. Another formulation using the soil dielectric constant instead of the volumetric water content was proposed by Holmes et al. (2006). The performance of some of these approaches is analysed in Wigneron et al. (2008).

The effect of soil surface roughness on the brightness temperature has been an issue widely addressed in the literature (Choudhury et al., 1979, Wang & Choudhury, 1981, Mo & Schmugge, 1987, Wang, 1983, Wegmüller & Mätzler, 1999, Wigneron et al., 2001, Schneeberger et al., 2004, Escorihuela et al., 2007). Fung (1994) proposes a theoretical physical model based on surface characteristics derived from the measured soil height profile. A simple empirical roughness model which takes into account only the coherent term of the scattering was reported in Choudhury et al. (1979),

Γsp=Γs,p*exp[(2kσscosθ)2]
(10)

where Γs,p* is the Fresnel specular reflectivity of soil, k is the electromagnetic wave number, σ s is the standard deviation of the surface height, and θ is the incidence angle. This model was reviewed, and another formulation was proposed in Wang & Choudhury (1981):

Γsp=[(1Qs)Γs,p*+QsΓs,q*]exp(hscosθn)
(11)

In this case, two semi-empirical parameters (h s and n) were included to model the effects of the polarisation mixing Q s, and surface roughness. The dependence of these parameters on surface properties such as correlation length or standard deviation of height is not yet clear. Mo & Schmugge (1987) and Wigneron et al. (2001) conclude that the n=2 dependence in (12) is too strong for L-band. A value of n=0 at both polarisations was found to be consistent with measurements in Wigneron et al. (2001), while Escorihuela et al. (2007) distinguishes n values for both polarisations (1 at horizontal and -1 at vertical). Similarly, there are discrepancies on the value of roughness parameter h s. Some authors obtain h s from experimental data by best-fit (Wigneron et al., 2001) while others propose expressions for h s as a function of geophysical parameters. Mo & Schmugge (1987) obtain good results with two parameterisations of h s as a function of σ s and lc. Finally, there is a general agreement on the value of the cross-polarisation parameter Q s, which has been found to be very small (0 to 0.12) at L-band.

Apart from these considerations, the effects of frequency and incidence angle on the roughness parameters have not been studied thoroughly. Mo & Schmugge (1987) and Shi et al. (2002) pointed out that the roughness effects depend on both the frequency and the incidence angle. The later suggested a parameterisation of the surface reflectivity derived from data simulated for a wide range of soil water content and roughness properties using the integral equation model (Fung, 1994). The surface reflectivity model of Shi et al. (2002) was tested in Schneeberger et al. (2004) and found not to be capable of explaining discrepancies between the ground truth and remotely sensed data. As a consequence, a new model was developed for describing the influence of the topsoil structure on the L-band emission as an impedance matching between the dielectric constants of soil and air (Mätzler, 2006).

If the remote sensor is placed above a canopy looking downwards, the measured brightness temperature will contain not only information on the soil, but also on vegetation, since vegetation radiates its own energy and, moreover, attenuates and scatters the soil radiation. Chukhlantsev (2006) revises the theory and conducted experimental research over vegetated areas. Although the modelling of the land emission involves analytical solutions of the radiative transfer equation (Ferrazzoli et al., 2002), this approach is not easy to use with experimental data. Hence, the common practice is to use approximate formulas or semi-empirical models in which the different components of the brightness temperature could be differentiated. The brightness temperature of a soil covered by vegetation is usually estimated as the contribution of three terms: (i) the radiation from the soil that is attenuated by the overlying vegetation, (ii) the upward radiation from the vegetation, and (iii) the downward radiation from the vegetation, reflected by the soil, and attenuated by the canopy (Ulaby et al., 1986):

TB,pmodel=(1+1espLveg)(11Lveg)(1ω)Tveg+espLvegTs
(12)

where Tveg and Ts are the physical temperatures of the vegetation and soil, respectively, Lveg=exp(τ secθ ) is the attenuation due to the vegetation cover, τ=bVWC is the optical thickness, b is the b-factor (Van de Griend & Wigneron, 2004), VWC is the vegetation water content, and ω is the single scattering albedo. This formulation is known as the τ−ω model and it is based on the single scattering approach proposed in Kirdiashev et al. (1979). The optical depth is related to the vegetation density and the frequency. The single scattering albedo describes the scattering of the emitted radiation by the vegetation, and is a function of plant geometry.

3. Soil moisture retrieval algorithms

The brightness temperature of land covers is influenced by many variables, the most important being soil moisture and temperature, and vegetation characteristics. The challenge is to reconstruct the environmental parameters from the measured signal by using a minimum of ancillary data. To do this, different soil moisture retrieval algorithms have been developed. The first one is based on the experimental relationship between the geophysical variables and the radiative transfer equation using a regression technique. This approach has limited applicability, since often the regression is valid only for the test sites where they were obtained. The second approach is based on the use of neural networks. These algorithms have been used with satisfactory results in the retrieval of agricultural parameters from radiometric data (Del Frate et al., 2003), but need a training phase that is not always feasible. The third type of algorithms is widely used and is based on the inversion of radiative transfer models. Obviously, this approach has also disadvantages, since errors of the model lead to errors in the retrieval. The soil moisture models are used as forward models, and the geophysical variables are retrieved by minimisation of a cost function of the type (Wigneron et al., 2003, Pardé et al., 2004, Saleh et al., 2006)

=θ(TBhTBhmeas)2+(TBvTBvmeas)2σTB2+N(PNiniPN)2σPN2
(13)

The simulated brightness temperature TBp is computed using (13), TBpmeas is the measured brightness temperature, and P N is any of the parameters on which TBp depends. A first-guess value of the parameter P N ( PNini ) with associated standard deviation σ PN can be also considered in the cost function.

4. Soil moisture field experiments

Space missions for Earth observation require a huge amount of previous work to validate the existing science and technology. This is accomplished through airborne or ground-based field experiments. Table 1 summarises some of the campaigns that have been conducted in the last years over land surfaces using L-band radiometers. Most of these radiometers have been designed within the context of preparatory activities previous to a space mission launch.

In the early 70’s NASA was getting ready to launch its first microwave radiometers for Earth observations on the Nimbus 5 spacecraft and it needed to learn more about the capabilities of these sensors for land observations. To do this several field campaigns were conducted in the period 1971 to 1973. For the remote sensing of soil moisture irrigated agriculture test areas were selected in the southwestern part of the United States. These experiments demonstrated the sensitivity of microwave radiometers to surface soil moisture and that the longer wavelength sensors worked best (Schmugge et al. 1974, Burke et al., 1979). These data were used to study the effects of surface roughness (Choudhury et al., 1979), soil temperature (Choudhury et al., 1982) and soil texture (Schmugge, 1980). These experiments defined the basic sensitivity of microwave radiometers for soil moisture sensing and some of the factors affecting this sensitivity. Since that time much has been done to refine and improve this understanding. Some of these later experiments will be decribed in the following sections. A description of the MOUSE 2004, T-REX 2004/2006, and SMOS REFLEX 2003/2006 campaigns over land carried out by UPC and of their results is provided in the next sections. Radiometric measurements were acquired using the UPC L-band Automatic Radiometer (LAURA), which has a working frequency of 1.4135 GHz, the same as the SMOS payload MIRAS (Villarino, 2004). A view of the experiment sites and the LAURA radiometer is shown in Figure 2.

media/image19.jpg

Figure 2.

Picture of the MOUSE 2004, T-REX 2006, and SMOS REFLEX 2003 experiment sites.

YearCampaignSiteSurfaceInstrument
1971/1975Phoenix,AZ Imperial valley, CA (USA)Irrigated agricultureMulti-frequency Microwave Radiometers
1976/1978Joint US/USSR exchangeHand County, SD (USA)Dryland agriculturePMIS, MFMR
1985/1989FIFEKonza Prairie, Kansas (USA)Short grass, prairiePBMR
1990/1991Monsoon 90Walnut Gulch AZ (USA)DesertPBMR, ESTAR
1992HAPEX-SahelNiger, West AfricaSemi-arid savannahPBMR
1992/1994WashitaOklahoma (USA)AgriculturalESTAR, AIRSAR, PBMR
1997/1999SGPOklahoma (USA)AgriculturalESTAR, PSR, PAL
1999ESSCReading (UK)Bare soilSWaMP-L
2001INRAAvignon (FR)Bare, cornEMIRAD
2001EuroSTARSSToulouse (FR), Utiel (SP)VariousSTARSS
2002ETHEschikon (CH)BareELBARA
2002SMEXIowa (USA)VariousESTAR, PSR, PALS, AIRSAR
2003/2006SMOS REFLEXUtiel (SP)VinesLAURA
2003ELBARA 2003Wageningen (NL)Bare, grassELBARA
2003-2006SMOSREXToulouse (FR)Bare, fallowLEWIS
2003SMEXOK,AL,GA(USA)Various2DSTAR, PSR, AIRSAR
2004SMEXAR, NM (USA)Various2DSTAR, PSR
2004MOUSEIspra (IT)Bare soilLAURA
2004/2006T-REXAgramunt (SP)Bare soilLAURA
2004ETHEschikon (CH)GrassELBARA
2004BrayLes Landes (FR)ForestEMIRAD
2005ETHJulich (GE)ForestELBARA
2005SMEXIowa (USA)VariousAPMIR
2005NAFEUpper Hunter (AU)Grass, forestEMIRAD, PLMR
2006TuRTLEEl Brull (SP)MountainLAURA
2006/2007MELBEXUtiel (SP)Shrub, vinesEMIRAD
2006NAFEYanco (AU)VariousPSR, PALS
2007CLASICOklahoma (USA)VariousPSR, PALS
2008-2010GRAJOVadillo de la Guareña (SP)Grass, barley, bare soilLAURA,SMIGOL

Table 1.

Summary of L-band field campaigns over land performed since 1970. Acronyms stand for USA (United States of America), UK (United Kingdom), FR (France), SP (Spain), NL (The Netherlands), CH (Switzerland), GE (Germany), IT (Italy), and AU (Australia).

5. Texture effects on the L-band emissivity of bare soils: the MOUSE 2004 experiment

The Monitoring Underground Soil Experiment (MOUSE) 2004 was performed from day of year (DoY) 159 to 184, 2004, at the Joint Research Centre (JRC) complex of the European Union in Ispra, north of Italy (4548N, 837E; 213 m altitude). The test lane consisted of six 6 m long by 6 m wide by 1 m depth bare plots of different soil types which are distributed along a lane with West-East orientation. The LAURA radiometer, an infra-red video-camera, and the control rack were installed on a metallic gantry guided by a rail on one side and freely rolling on the other. A picture of the setup is shown in Figure 2(c). Soil temperature was measured using 107 probes from Campbell Scientifc installed just below the surface, and at 5 cm and 10 cm depth. An infra-red video camera was mounted by LAURA to acquire soil surface temperature images concurrently with the radiometric measurements. Most probably, satellite-borne radiometers will have to rely on this kind of auxiliary sensors to estimate soil moisture from radiometric observations. On the other hand, soil moisture was automatically acquired every 30 min by seven ML2x Thetaprobes per plot installed below the surface (x2), and at 5 cm (x2), 10 cm (x2), and 15 cm (x1). Plots were irrigated three times during the experiment, except RS which was kept dry and protected from water by a sliding tent. Soils at the experiment site were loam (LO), natural sand (SA), sieved sand used in the construction industry (RS), clay from rice fields (CL), a mixture of loam and of commercial products for gardening (OR), and ferromagnetic crushed volcanic soil from Naples (FE).

Figure 3 represents the sensitivity of the H-pol emissivity to the percentage of field capacity, which depends on soil texture. Measurements at 25 (blue), 35 (red), and 45 (green) incidence angles have been plotted. LO, RS, and FE emissivity has a linear dependence on field capacity. The R2 estimator at 35 for LO is 71%, for RS 70%, and for FE is 85%. Clay has the worst linear fit, with R2 lower than 20%. Soil samples of each plot were acquired and their dielectric constant measured at laboratory (Vall-llossera et al., 2005a).

Soil moisture was retrieved from radiometric measurements using the cost function in (14) with different constraints. The minimum soil moisture RMSE obtained considering different depths as ground-truth data and three different dielectric constant models (Wang & Schmugge, 1980, Dobson et al., 1985; Vall-llossera et al., 2005a) is represented in Figure 4. Best results for LO and FE are achieved for the soil moisture in the 0 to 5 cm layer, but few discrepancies between RMSE values at different depths were found out. This is mainly due to their almost linear soil moisture profile. Best results for SA were obtained for the soil moisture in the 0-15 cm layer. As sandy soils dry faster, the penetration depth is higher. On the other hand, as RS plot was kept dry throughout the experiment, no important differences in the RMSE have been found in the comparison. Clay had a high water content at 10 cm (>0.22 m3/m3) and hence the main signal contribution comes from the first 10 cm of soil. Similar RMSE values have been found for OR considering either the 0 to 5 cm layer or the 0 to 10 cm layer as ground-truth. Similarities are because the water content at 5 cm and 10 cm depth is almost the same during the experiment. On the other side, in Figure 4 it can be noted the performance of each dielectric constant models for each soil type. Results for Wang & Schmugge (1980) model are good for loamy and ferromagnetic soils, but their performance for sandy soils is worse than Dobson et al. (1985) and Vall-llossera et al. (2005a). The model in Dobson et al. (1985) offers best results for sandy soils (SA and RS) which registered lower soil moisture values and dried faster than the others. The expressions derived from UPC laboratory experiments show good performance for most of the fields. These facts confirm that the selected dielectric constant model has an impact on soil moisture estimation error and can increase it between 3 to 5% depending on the soil type.

media/image20.jpg

Figure 3.

Variation of the emissivity with the field capacity at H-pol. Colours indicate the incidence angle: 25 (blue), 35 (red), and 45 (green).

media/image21.jpg

Figure 4.

Minimum root mean square error of the soil moisture retrieval. Results for a soil dielectric constant model and different soil moisture ground-truth data are represented.

6. Emission of bare rough soils: the T-REX experiments

Several experiments have been performed in the past years to study the L-band emission from bare soils, both at small and large scales. Although small scale measurements are not representative of onboard sensor data, they can be useful to evaluate the various approaches to the land emission modelling. The Terrain-Roughness EXperiments (T-REX) were carried out in Agramunt, Spain (41 48N, 1 7W; 356 m altitude) from DoY 333 to DoY 337, 2004, and then again from DoY 143 to DoY 160, 2006. The sites consisted of bare soil plots with four different ploughing, but the same textural composition (21% clay, 21.3% sand). This section presents some results from the T-REX 2006 experiment. A picture of the site has been shown in Figure 2(a). Each parcel was 6 m x 15 m and the same ploughing was applied to two different plots to observe spatial diversity effects. LAURA was installed on a trailer towed by a car. Radiometric measurements were acquired at incidence angles from 40 to 65 in T-REX 2006, in 5 steps. Every plough type was measured twice each day of experiment.

Parameter (mm)Plot label
ABCD
lc12610019096
σs1716.6833
σλ76.82.517
hair-soil31261346

Table 2.

Roughness characteristics of the T-REX 2006 ploughs.

Roughness measurements were performed in the along and across directions of LAURA's field of view. The standard deviation of height (σ s) and the correlation length (l c) of the ploughs at the T-REX 2006 site are summarised in Table 2. Two additional parameters have also been included: (i) h air-soil, which is the transition layer thickness used to evaluate the small-scale soil roughness in the model proposed in Schneeberger et al. (2004), Mätlzer (2006, section 4.7) and (ii) σ λ , which is the standard deviation of height assuming a windowing of length λ. Soil temperature was measured at 0, 5, 10, 15, and 20 cm depth using thermometers. On the other hand, soil moisture was measured with ML2x ThetaProbe sensors at various positions within the radiometer's FOV. As the purpose of this experiment was to study the impact of soil roughness, the site was not irrigated to keep as constant a value of soil moisture as possible. The mean soil moisture during T-REX 2006 was of 4% to 6%.

Figure 5 shows the variation of the emissivity measured during the T-REX 2006 as a function of the incidence angle. Colours indicate the polarisation, H-pol (red) and V-pol (blue), while line styles distinguish the ploughing. While soil emission at H-pol decreases almost linearly with increasing incidence angle (and the trend is more significative at incidence angles above 55 ), the emission at V-pol slightly increases with incidence angle until an inflexion point at 55 , when it begins decreasing. Ploughs C and D, which correspond to the smoother and rougher plots in the site, have the highest (lowest) range of variation of the emissivity.

The comparison between radiometric measurements and numerical simulations obtained with a numerical approach by means of the integral equation method (IEM; Fung (1994)) are shown in Figure 6. Icons stand for the mean value of LAURA measurements while error bars indicate their standard deviation during the whole T-REX 2006 experiment. Although IEM predicts an increase of the vertical emission with incidence angle, dry soil measurements decrease as incidence angle increases. A good matching between simulations and measurements is observed at horizontal polarisation. The effects of taking into account the incoherent component are only relevant in the case of very rough ploughs (see Figure 6(c), plough D).

media/image22.jpeg

Figure 5.

Variation of the emissivity measured during the T-REX 2006 as a function of the incidence angle for the four ploughs. Colours indicate the polarisation: H-pol (red) and V-pol (blue). Line style and icons distinguish the ploughing.

media/image23.jpg

Figure 6.

Variation with time of the emissivity measured during the T-REX 2006 experiment as a function of the incidence angle. The mean value of LAURA measurements (black icons) and their standard deviation during the whole experiment (bars) have been represented together with IEM simulations.

media/image24.jpg

Figure 7.

Comparison between T-REX 2006 measurements (black icons) and the models in Wang & Choudhury (1981), Mätzler (2006, section 4.7), Schneeberger et al. (2004), Mo & Schmugge (1987), and Escorihuela et al. (2007). Solid lines indicate results using the standard deviation of the surface height profile while dashed lines indicates results using the λ-average standard deviation of height

T-REX 2006 measurements were also compared to the semi-empirical soil reflectivity models in Wang & Choudhury (1981), Mo & Schmugge (1987), Escorihuela et al. (2007), and the air-to-soil transition model proposed in Schneeberger et al. (2004). These models propose different expressions for the effective roughness parameter h s in (12). Results are shown in Figure 7. Upper plots correspond to vertical polarisation while bottom plots correspond to horizontal polarisation. Colours indicate the soil reflectivity model, while line types indicate whether l c and σ s of the height profile (solid lines) or their equivalents taking into account the λ-windowing (dashed lines) have been used. Both the model in Wang & Choudhury (1981) using as inputs the λ -windowing soil roughness characteristics and the air-to-soil transition model from Schneeberger et al. (2004) fit data. However, Schneeberger et al. (2004) does not fit T-REX 2006 measurements of wet soil, which suggests that inhomogeneities in soil moisture or the presence of dew make it difficult to properly characterise the upper soil layer, which is needed for the air-to-soil transition model.

7. Radiometric observations of fully developed vines

Several papers are available that present experimental studies of brightness temperatures from vegetation canopies (Jackson et al. 1982, Jackson et al., 1991, Ferrazzoli et al., 2002, Van de Griend & Wigneron, 2004, Saleh et al., 2006). A look at the field experiments conducted during the past decades over vegetated sites shows how they focused their attention on crops such as alfalfa, wheat corn, or soybean, on grass, on forests, and on bushes, but few studies had been conducted over vineyards. In this context, two ground-based field experiments were performed to study the variations on the emissivity due to vines development. The SMOS REFerence pixel L-band EXperiments (REFLEX) were carried out in July 2003 and then again from July to November 2006 at two vineyards within the València Anchor Station (VAS; Figure 2(b)), Spain (Vall-llossera et al., 2005b). The VAS has been selected as a SMOS Cal/Val site because of its almost homogeneous land cover (vineyards) and its size, which is approximately that of a SMOS pixel. The goal of the SMOS REFLEX experiments was to assess the impact of grapevines on the radiometric emission and on the sensitivity of the sensor to soil moisture. Results from the SMOS REFLEX 2003 experiment are presented hereafter.

media/image25.jpg

Figure 8.

Mean value and standard deviation of the emissivity and of the ground-truth soil moisture as a function of (a) the incidence angle and for all azimuth angles, and (b) the azimuth angle and for all incidence angles. DoY 185 (solid line, wet soil) and DoY 181 (dashed-dotted line, dry soil) have been represented. Blue and red indicate vertical and horizontal polarisations, respectively.

During SMOS REFLEX 2003 radiometric measurements were acquired from DoY 181 to 191 at nine incidence angles (from 25 to 65 in 5 steps) and seven azimuth angles. Pictures at every look direction were taken to estimate the percentage of FOV covered by plants. The site was irrigated until saturation on DoY 182 and 185, and then was let to dry out. Even though the field was irrigated until saturation, different soil moisture values were measured depending on the location of the test point within the vineyard. This is most probably because of surface inhomogeneities and variations on the compactness of the terrain due to plough and roots distribution. This large spatial variability in soil moisture, which in this case varies from 3% up to 35% in a 20 m x 20 m area, is one of the main problems soil moisture retrieval algorithms must face since it complicates the comparison between ground-truth and estimated ws at the plot scale. Volumetric soil moisture, soil temperature, and soil roughness were measured, and plants were fully characterised.

Figure 8(a) represents the mean value and standard deviation of the emissivity and of the ground-truth soil moisture as a function of the incidence angle and for all azimuth angles. Measurements from DoY 185 (solid line) and DoY 181 (dashed-dotted line) have been represented. The mean value of both polarisations seems independent of incidence angle for completely dry soils, which is the case of DoY 181, being the std. of measurements lower than 0.005 (0.01) in emissivity -1.5 K (3 K) in brightness temperature- for vertical (horizontal) polarisations. In contrast, for wet soil an ascending trend is observed for V-pol emissivity up to 45 , when V-pol emissivity seems to saturate, while H-pol emissivity decreases as incidence angle increases. The std. of measurements for wet soils is 0.025 (0.03) for vertical (horizontal) polarisations.

The mean value and standard deviation of the emissivity from DoY 181 and DoY 185 as a function of the azimuth angle and for all incidence angles is represented in Figure 8(b). For dry soils no dependence on the azimuth angle is observed, being the standard deviation of measurements lower than 1.5 K at both polarisations. In the case of wet soils, and since the contribution from vines is supposed to be the same for specular azimuth angles with respect to zero, the inhomogeneity in soil moisture content, which is higher for positive azimuth angles, leads to an asymmetry in the emission measurements of DoY 185 (solid lines), which as expected are higher for the lowest soil moisture. For a given azimuth angle, the standard deviation of the measured emissivity at V-pol goes from 0.014 (4 K) at 0 down to 0.008 (2 K) at 65 , whereas the largest standard deviation at H-pol goes from 0.007 (2.6 K) at 0 up to 0.022 (7 K) at 45 .

Figure 9(a) shows the variation of the retrieved albedo with respect to the incidence angle for all azimuth angles. Icons indicate the day of experiment, while the solid line is the mean value of the whole experiment at each incidence angle. The albedo is smaller than 0.15, having its maximum for an incidence angle of 35 and then decreasing as the incidence angle increases. On the other hand, the opacity tends to increase with incidence angle as shown in Figure 9(b) which is in accordance to the higher fraction of area covered by vegetation as the incidence angle. Good agreement between measured and computed values is appreciated for incidence angles below 55 . For 55 the algorithm does not always converge and at 60 and 65 , the convergence is rarely achieved, pointing out a problem in the emissivity model function for large incidence angles. This discrepancy can be due to the fact that as the incidence angle increases, the signal from the soil suffers from more attenuation and scattering in the canopy. The τ−ω model is a valid and simple forward model to be used in optimal estimation approaches, but improvements are required if it is to be used at high incidence angles.

media/image26.jpg

Figure 9.

Dependence of the retrieved vegetation albedo (a), and opacity (b) on the incidence angle. Icons represent the value estimated from LAURA measurements for each day of experiment and all azimuth angles. The mean value is represented with a black solid line. (c) Ground truth soil moisture vs. estimates from radiometric measurements.

In order to retrieve the soil water content from the LAURA's measurements, an iterative algorithm, based on the least squares minimisation procedure, was applied. A better performance was found out when the albedo, opacity and soil moisture were simultaneously retrieved using (14). The first guesses for the opacity and the albedo are the estimated values obtained in the iterative algorithm described before. The algorithm converges for incidence angles equal to or smaller than 50 , only for some day of experiments at 55 , and rarely for 60 or 65 due to the increasing discrepancy between the model and the data as the incidence angle increases. When the incidence angle increases, the canopy influence becomes more important (higher attenuation on the signal and higher scattering contribution), so a higher order approximation for the direct model becomes necessary. In Figure 9(c), the soil moisture estimated from LAURA data is compared to ground-truth soil moisture data. The slope of the computed linear regression is 0.9 (close to 1), which shows a good agreement between both values. The error between estimations and measurements is about 2.3%, which is smaller than the 4% accuracy required for SMOS.

8. Conclusions

The basics of land emission modelling have been summarised in this chapter, and an overview of the field experiments over land carried out during the last three decades has been given to present to the reader the state of the art. Some results of the recent SMOS REFLEX 2003, MOUSE 2004, and T-REX 2006 experiments have been reported here. These are, among others, part of the preparatory experiments for the SMOS mission, which in the near future will provide soil moisture maps at global scale.

Soil texture influences the electromagnetic behaviour of soils since it affects their dielectric constant. During MOUSE 2004, six types of bare at soils with different soil moisture contents

were measured at five incidence angles. As expected, variations on the moisture profile were noted depending on texture. The emission of loams and sands showed a linear dependence on the percentage of field capacity, especially at H-pol. In this case, the R-squared estimator at 35 was found to be between 70% and 85%. On the contrary, clay showed no linear trend which, however, can not be a general conclusion since clay data had a large scattering.

The modified Wang & Choudhury (1981) soil emission model, with Q s=0 and n=0, was used as forward model in the soil moisture retrieval algorithm. Best results were obtained when no constraints were applied to soil moisture and temperature. The estimated soil moisture was compared with the soil moisture measured at different depths. For loam and ferromagnetic soils best results were achieved considering as ground-truth data the soil moisture in the 0 to 5 cm layer, 0 to15 cm layer for sands, and 0 to 10 cm layer for clays. On the other hand, when results for every dielectric constant model were compared, it was found out that the Wang & Schmugge (1980) model had a better performance in the retrieval algorithm of loam and ferromagnetic soils, while the Dobson et al. (1985) model offered best results for sands.

Results from MOUSE 2004 suggest that the soil dielectric constant model should be selected as a function of texture, since the root mean squared error can vary from 2% to 7% depending on the soil texture, and that soil moisture estimates will be the averaged value in a soil layer which will depend on texture.

Soil roughness has a strong impact on land brightness temperature. This effect is more noticeable in the case of dry soils. Data from bare soils with a standard deviation of height varying from 8 mm to 33 mm was acquired during the T-REX experiments, and was compared to predictions from numerical and semi-empirical soil emission models. The goodness of these models is of key importance for an accurate soil moisture estimation from satellite data. As expected, dry soil emissivity at H-pol decreased almost linearly with increasing incidence angle, being the decrement inversely proportional to soil roughness. Although similar emission was measured at H-pol for all plough at low incidence angles, it was noted that emission decreases slowly for rough soils than for smooth soils. Above an incidence angle of 50 , the decreasing slope of H-pol emission was higher for all plough. The sensitivity of V-pol to roughness is lower than that of H-pol. The trend of V-pol emission was found to be different depending on whether the soil was wet or dry. When the soil is wet the emission increased with the incidence angle, which is in accordance with most model predictions. However, the trend for dry soils is decreasing with increasing incidence angle and roughness. This may suggest the existence of a relationship between soil moisture and the effective roughness.

The integral equation method (IEM) predicts an increase in V-pol emission with incidence angle, which is not in accordance with measurements of dry soils. On the other hand, the predicted descending trend at H-pol with increasing incidence angle is in accordance to measurements, although IEM underestimates the soil emission both for dry and wet soils. In a randomly ploughed field, without a significant tillage direction, the impact of choosing an exponential or gaussian height probability density function in the IEM model is minimum, whereas it was noted that the incoherent term of the reflectivity must be considered for rough soils.

In general, all semi-empirical land emission models follow the trend of dry soils measurements at H-pol, whereas discrepancies exist for wet soils. Neither the semi-empirical models nor the IEM describe the trend of dry soils V-pol measurements, being the error at this polarisation larger as the incidence angle increases. The lowest error between predictions and measurements for H-pol was obtained using the Wang & Choudhury (1981) model, but with roughness statistics averaged as a function of the measurements wavelength. The dependence on the incidence angle proposed by Wang & Choudhury (a squared cosine) had been discussed by other authors who considered it to be too much strong at L-band. However, this simple formulation has been tested with T-REX data with good results if the averaged standard deviation in transects equal to the wavelength is used.

The SMOS REFLEX 2003/2006 experiments site was in the Valencia Anchor Station, a selected area for the SMOS calibration and validation activities. No previous studies over vines were reported in the literature prior to these campaigns. In the first experiment, fully developed vines were characterised during two weeks, while controlled irrigations moistened the field. The second experiment was planned to monitor changes in the L-band emission of vineyards during different stages of plants development.

Since vines do not have a predominant vertical nor horizontal structure, the opacity and the albedo were found to be independent on the polarisation. The Wang & Choudhury (1981) model was used as the soil emission forward model in the retrieval algorithm, while the dielectric constant model was the one from Wang & Schmugge (1980). Good results were obtained for incidence angles up to 55 , but the convergence of the algorithm was rarely achieved above that value. The error between ground-truth and estimated soil moisture was 2.3%, better than the 4% required for SMOS. For incidence angles above 55 the convergence of the algorithm was rarely achieved, probably due to the larger effect of the vegetation at large incidence angles, not accurately described by the model. Higher order models should be accounted for at these angles.

Although the radiometric behaviour varies from one canopy to the other, at a larger scale such as that of SMOS it is very likely that the vegetation types can be averaged over the footprint and that it is not necessary to account for an accurate distinction between canopies.

. Acknowledgements

A. Monerris would like to thank Spanish Ministry of Science and Education projects MIDAS-2 ESP2002-11648-E, MIDAS-3 ESP2004-00671, and MIDAS-4 ESP2005-06823-C05-02, for funding the SMOS REFLEX, T-REX and MOUSE 2004 field experiments, and the FPU AP2003-1567 grant. Special thanks are also given to M. Vall-llossera and A. Camps, from UPC, for their valuable comments during the analysis of the experimental data.

References

1 - H. Barré, B. Duesmann, Y. H. Kerr, 2008 SMOS: The mission and the system. IEEE Trans. Geosci. and Remote Sens., 46 587 593
2 - J. Behari, 2005 Microwave dielectric behaviour of wet soils (Remote Sensing and Digital Image Processing), Springer, 101402032714 Netherlands
3 - N. C. Brady, R. R. Weil, 2003 Elements of the nature and properties of soils, 2nd ed., Prentice Hall, 10 NJ, USA
4 - W. J. Burke, T. J. Schmugge, J. F. Paris, (1979). Comparison of 2.8 and 21 cm microwave radiometer observations over soils with emission model calculations, J. Geophys. Res., 84 287 294 .
5 - A. Chanzy, S. Raju, J. P. Wigneron, 1997 Estimation of soil microwave effective temperature at L and C bands. IEEE Trans. Geosci. and Remote Sens., 35 3 570 580
6 - B. J. Choudhury, T. J. Schmugge, A. T. C. Chang, R. W. Newton, 1979 Effect of surface roughness on the microwave emission of soils. J. Geophys. Res., 84 5699 5705
7 - B. J. Choudhury, T. J. Schmugge, T. Mo, 1982 A simple parameterization of effective soil temperature for microwave emission. J. Geophys. Res., Vol., No., (month and year of the edition) 1301 1304
8 - A. A. Chukhlantsev, 2006 Microwave radiometry of vegetation canopies, Springer, 101402046812 Netherlands
9 - F. Del Frate, P. Ferrazzoli, G. Schiavon, 2003 Retrieving soil moisture and agriculturalvariables by microwave radiometry using neural networks. Rem. Sens. Env., 84 2 174 183
10 - M. C. Dobson, F. T. Ulaby, M. T. Hallikainen, M. A. Elrayes, 1985 Microwave dielectric behaviour of wet soils 2. Dielectric mixing models. IEEE Trans. Geosci. and Remote Sens., 23 35 46
11 - J. R. Eagleman, W. C. Lin, 1976 Remote sensing of soil moisture by a 21 cm passive radiometer.. J. Geophys. Res., 81 3660 3666
12 - D. Entekhabi, E. Njoku, P. Houser, M. Spencer, T. Doiron, J. Smith, R. Girard, S. Belair, W. Crow, T. Jackson, Y. Kerr, J. Kimball, R. Koster, K. Mc Donald, P. O’Neill, T. Pultz, S. Runnig, J. Shi, E. Wood, J. van Zyl, 2004 The Hydrosphere State (HYDROS) mission concept: An Earth system pathfinder for global mapping of soil moisture and land freeze/thaw. IEEE Trans. Geosci. and Remote Sens., 42 2184 2195
13 - M. J. Escorihuela, Y. H. Kerr, P. de Rosnay, J. P. Wigneron, J. C. Calvet, F. Lemaitre, 2007 A simple model of the bare soil microwave emission at L-band, IEEE Trans. Geosci. and Remote Sens., 45 1978 1987
14 - P. Ferrazzoli, L. Guerriero, J. P. Wigneron, 2002 Simulating L-band emission of forests in view of future satellite applications. IEEE Trans. Geosci. and Remote Sens., 40 2700 2708
15 - A. K. Fung, 1994 Microwave Scattering and Emission. Models and their applications, Artech House, 100890065233 MA.
16 - M. T. Hallikainen, F. T. Ulaby, M. C. Dobson, M. A. El -Rayes, L. Wu, 1985 Microwave dielectric behaviour of wet soil- Part 1: Empirical models and experimental observations. IEEE Trans. Geosci. and Remote Sens., GE-23 25 34
17 - J. E. Hipp, 1974 Soil electromagnetic parameters as a function of frequency, soil density and soil moisture, Proceedings of IEEE 98 103
18 - T. Holmes, P. de Rosnay, R. de Jeu, J. P. Wigneron, Y. Kerr, J. C. Calvet, M. J. Escorihuela, K. Saleh, F. Lemaitre, 2006 A new parameterization of the effective temperature for L- band radiometry.. Geophys. Res. Lett., 33, 7405
19 - T. J. Jackson, T. J. Schmugge, J. R. Wang, 1982 Passive microwave sensing of soil moisture under vegetation canopies.. Water Resources Research, 18 1137 1142
20 - T. J. Jackson, T. J. Schmugge, 1989 Passive microwave remote sensing system for soil moisture: Some supporting research.. IEEE Trans. Geosci. and Remote Sens., 27 2 225 235
21 - T. J. Jackson, T. J. Schmugge, 1991 Vegetation effects on the microwave emission from soils. Remote Sens. Environ., 36 203 212
22 - Y. H. Kerr, P. Waldteufel, J. P. Wigneron, J. M. Martinuzzi, J. Font, M. Berger, 2001 Soil moisture retrieval from space: the Soil Moisture and Ocean Salinity (SMOS) mission.. IEEE Trans. Geosci. and Remote Sens., 39 1729 1735
23 - K. P. Kirdiashev, A. A. Chukhlantsev, A. M. Shutko, 1979 Microwave radiation of the Earth’s surface in the presence of vegetation cover. Radio Eng. Electron. Phys., 2 37 56 .
24 - R. D. Koster, P. A. Dirmeyer, Z. Guo, G. Bonan, E. Chan, et al. 2004 Regions of strong coupling between soil moisture and precipitation.. Science, 305 5687 1138 1140
25 - C. Mätlzler, A. Standley, 2000 Relief effects for passive microwave remote sensing, J. Geophys. Res., 21 2403 2412
26 - C. Mätzler, 2006 Thermal Microwave Radiation: Applications for Remote Sensing, IET Electromagnetic Waves Series 52, 100863415733 UK
27 - K. D. Mc Mullan, M. A. Brown, M. Martín-Neira, W. Rits, J. Martí, J. Lemanczyk, 2008 SMOS: The payload. IEEE Trans. Geosci. and Remote Sens., 46 3 594 605
28 - V. L. Mironov, M. C. Dobson, V. H. Kaupp, S. A. Komarov, V. N. Kleshchenko, 2004 Generalized refractive mixing dielectric model for moist soils. IEEE Trans. Geosci. and Remote Sens., 42 4 773 785
29 - T. Mo, T. Schmugge, 1987 A parameterization ot the effect of surface roughness on microwave emission., IEEE Trans. Geosci. and Remote Sens., 25 47 54
30 - A. Monerris, P. Benedicto, M. Vall-llossera, E. Santanach, M. Piles, R. Prehn, 2008 Assessment of the topography impact on microwave radiometry at L-band. J. Geophys. Res., 113(B12202)
31 - M. Pardé, J. P. Wigneron, P. Waldteufel, Y. H. Kerr, A. Chanzy, S. S. Sobjaerg, N. Skou, 2004 N-parameter retrievals from L-band microwave observations acquired over a variety of crop fields. IEEE Trans. Geosci. and Remote Sens., 42 6 1168 1178
32 - C. H. Roth, M. A. Malicki, R. Plagge, 1992 Empirical evaluation of the relationship between soil dielectric constant and volumetric water content as the basis for calibrating soil moisture measurements by TDR. J. of Soil Sci., 43 1 3
33 - K. Saleh, J. P. Wigneron, P. de Rosnay, J. C. Calvet, M. J. Escorihuela, Y. Kerr, P. Waldteufel, 2006 Impact of rain interception by vegetation and mulch on the L-band emission of natural grass.. Remote Sens. Environ., 101 127 139
34 - T. Schmugge, P. Gloersen, T. T. Wilheit, F. Geiger, 1974Remote Sensing of Soil Moisture with Microwave Radiometers, NASA/GSFC X-652-72-305, Jour. Geophys. Res., 79(2), 317 323 .
35 - T. Schmugge, 1980Effect of Texture on the Microwave Emission from Soils, NASA Tech. Memo. 80632, IEEE Trans. on Geoscience and Remote Sensing, GE-18, 353 361 .
36 - K. Schneeberger, M. Schwank, Ch. Stamm, P. de Rosnay, Ch. Mätzler, H. Flühler, 2004 Topsoil structure influencing soil water retrieval by microwave radiometry.. Vadose Zone Journal, 3 1169 1179
37 - M. Schwank, M. Stähli, H. Wydler, J. Leuenberger, C. Mätzler, H. Flühler, 2004 Microwave L-band emission of freezing soil., IEEE Trans. Geosci. and Remote Sens., 42 1252 1261
38 - M. Talone, A. Camps, A. Monerris, M. Vall-llossera, P. Ferrazzoli, M. Piles, 2007 Surface topography and mixed pixel effects on the simulated L-band brightness temperatures. IEEE Trans. Geosci. and Remote Sens.,42 786 794
39 - G. C. Topp, J. L. Davies, A. P. Annan, 1980 Electromagnetic determination of soil water content: Measurements in coaxial transmission lines.. Water Resour. Res., 16 574 582
40 - F. T. Ulaby, R. K. Moore, A. K. Fung, 1986 Microwave Remote Sensing: Active and Passiv, vol. II, Radar remote sensing and surface scattering and emission theory, Addison-Wesley, 100890061912
41 - A. A. Van de Griend, J. P. Wigneron, 2004 The b-factor as a function of frequency andcanopy type at H-polarization. IEEE Trans. Geosci. and Remote Sens., 42 4 786 794
42 - M. Vall-llossera, M. Cardona, S. Blanch, A. Camps, A. Monerris, I. Corbella, F. Torres, N. Duffo, 2005a L-band dielectric properties of different soil types collected during the MOUSE 2004 field experiment. Proceedings IGARSS, 2 1109 1112
43 - M. Vall-llossera, A. Camps, I. Corbella, F. Torres, N. Duffo, A. Monerris, R. Sabia, D. Selva, C. Antolín, E. López-Baeza, J. F. Ferrer, K. Saleh, 2005b SMOS REFLEX 2003: L-band emissivity characterization of vineyards. IEEE Trans. Geosci. and Remote Sens., 43 5 973 982
44 - R. Villarino, 2004 Empirical Determination of the Sea Surface Emissivity at L-band: A contribution to ESA’s SMOS Earth Explorer Mission. Ph.D. thesis, RSLab, UPC. [Online] http://www.tesisenxarxa.net/TDX-1019104 170028 /index.html
45 - J. R. Wang, T. J. Schmugge, 1980 An empirical model for the complex dielectric permittivity of soils as a function of water content.. IEEE Trans. Geosci. and Remote Sens., 18 288 295
46 - J. R. Wang, 1983 Passive microwave sensing of soil moisture content: The effects of soil bulk density and surface roughness. . Remote Sens. Environ., 13 329 344
47 - J. R. Wang, B. J. Choudhury, 1981 Remote sensing of soil moisture content over bare field at 1.4 GHz frequency. J. Geophys. Res., 86 5277 5282
48 - U. Wegmüller, C. Mätzler, 1999 Rough bare soil reflectivity model. IEEE Trans. Geosci. and Remote Sens., 37 1391 1395
49 - J. Wigneron, L. Laguerre, Y. H. Kerr, 2001 A simple parameterization of the L-band microwave emission from rough agricultural soils. IEEE Trans. Geosci. and Remote Sens., 39 1697 1707
50 - J. P. Wigneron, J. . Calvet, T. Pellarin, A. A. de Griend, M. Van Berger, P. Ferrazzoli, 2003 Retrieving near-surface soil moisture from microwave radiometric observations: current status and future plans. Remote Sens. Environ., 85 489 506
51 - J.P. Wigneron, , A. Chanzy, P. de Rosnay, C. Räudiger, J. C. Calvet, 2008 Estimating the effective soil temperature at L-band as a function of soil properties. IEEE Trans. Geosci. And Remote Sens., 46 3 797 807