Operational Mapping of Soil Moisture Using Synthetic Aperture Radar Data: Application to the Touch Basin (France).

Soil moisture is a key parameter in different environmental applications, suchas hydrology and natural risk assessment. In this paper, surface soil moisture mappingwas carried out over a basin in France using satellite synthetic aperture radar (SAR)images acquired in 2006 and 2007 by C-band (5.3 GHz) sensors. The comparisonbetween soil moisture estimated from SAR data and in situ measurements shows goodagreement, with a mapping accuracy better than 3%. This result shows that themonitoring of soil moisture from SAR images is possible in operational phase. Moreover,moistures simulated by the operational Météo-France ISBA soil-vegetation-atmospheretransfer model in the SIM-Safran-ISBA-Modcou chain were compared to radar moistureestimates to validate its pertinence. The difference between ISBA simulations and radarestimates fluctuates between 0.4 and 10% (RMSE). The comparison between ISBA andgravimetric measurements of the 12 March 2007 shows a RMSE of about 6%. Generally,these results are very encouraging. Results show also that the soil moisture estimatedfrom SAR images is not correlated with the textural units defined in the European Soil Geographical Database (SGDBE) at 1:1000000 scale. However, dependence was observed between texture maps and ISBA moisture. This dependence is induced by the use of the texture map as an input parameter in the ISBA model. Even if this parameter is very important for soil moisture estimations, radar results shown that the textural map scale at 1:1000000 is not appropriate to differentiate moistures zones.


Introduction
Soil moisture plays a crucial role in the continental water cycle, more specifically in the distribution of precipitation between surface runoff and infiltration, which is the main driver behind most hydrological and geomorphologic processes [1][2][3]. The retrieval of soil moisture data is of considerable importance in various applications, such as hydrology, risk prediction, agriculture, and meteorology. It is a key indicator for constraining the initial conditions of infiltration/runoff rates when modeling overland flow. The knowledge of the initial water content at the beginning of the rainfall event is a prerequisite because all infiltration models integrate this information. In the context of implementing the Water Directive and the forthcoming Soil Framework Directive in Europe, there is a need for operational tools to evaluate land management scenarios and to provide sound references for targeting land use planning and the protection of natural resources. More specifically, operational services in charge of flood prevention and forecasting, such as the SCHAPI (National Flood Forecasting Agency) in France are investing in the continuous distributed modeling of soil surface moisture.
Only spatial remote sensing allows monitoring of environmental problems over large areas at regular intervals. Moreover, radar sensors allow mapping whatever the meteorological conditions (clouds, fog, etc.), both day and night. This is not the case with optical sensors, which are difficult to operate if there is cloud cover, a frequent situation in periods where flood events occur. The Soil Surface Characteristics (SSC) can be estimated from microwave remote sensing sensors due to the sensitivity of radar signal to soil characteristics such as the soil's roughness and dielectric constant [4][5][6]. In addition, the radar signal depends on various radar parameters such as the polarization, incidence angle and frequency.
The SAR sensors currently operational are PALSAR/ALOS, ASAR/ENVISAT, RADARSAT-1 and ERS-2. ERS-2 and RADARSAT-1 provide data based on a single polarization (VV for ERS-2, and HH for RADARSAT-1) whereas ASAR allows the scientific community to acquire images in dualpolarization mode (HH/HV, HH/VV, VV/VH). The radar incidence angle is 23° for ERS-2, between 20° and 50° for RADARSAT-1, and ranges from 14° to 45° for ASAR. The nominal swath width for ASAR and ERS-2 is greater or equal to 100 km, with a spatial resolution of, at best, 25 m (12.5 m pixel size). For RADARSAT-1, the nominal swath width is greater than or equal to 50 km, with a spatial resolution of, at best, 10 m (6.25 m pixel size). These three radar sensors operating in C-band (~5.3 GHz), have been widely used for retrieving soil moisture [e.g. 4,[7][8][9][10][11][12][13][14][15]. The SAR data are more and more used to derive soil parameters. Although many improvements have been achieved, this field is still not fully operational.
When using ERS-1/2, RADARSAT-1 and ASAR, a better estimate of soil moisture is obtained with a radar configuration that minimizes the effects of other soil surface characteristics (mainly surface roughness). Several studies have shown that the best C-band estimates of soil moisture when using only one radar channel are obtained with SAR images acquired at low and medium incidence angles [13,15,16]. These optimal configurations of incidence angle (≤37°) can be attained with current SAR sensors (ERS-2, RADARSAT-1, ASAR, and ALOS) about 10 times per month on a study site in Europe. The arrival of new sensors (TerraSAR-X, COSMO-SkyMed, and RADARSAT-2) should improve this temporal frequency.
The use of two SAR images acquired at both low and high incidence angles markedly improves the precision of the moisture estimate because the backscattering coefficients ratio eliminates the effects of roughness and thus permits linking of the radar backscattering coefficients to only the moisture [13,15,16]. This solution is not possible with current SAR sensors and will not be possible either with the sensors planned for the near future (ALOS, RADARSAT-2 and Terra SAR-X). Indeed, the time separating two SAR images acquired by the same sensor at two selected incidence angles (20° and 40°) is of the order of a few days (about 7 days with ASAR over the Touch basin), which limits the use of this possibility because the soil surface characteristics could change rapidly between the two acquisitions. Only the single incidence case is thus possible. The use of SAR images with two polarizations (case of ASAR sensor) does not provide a significant improvement in estimating soil moisture (improvement <1%, [15,17]) when compared with the results obtained with a single polarization.
Soil moisture can be measured (in situ or using remote sensors) and also modeled. In situ soil moisture measurements are performed using two techniques: TDR (Time Domain Reflectometry) probes and gravimetry. TDR measurements are quicker to carry out, but less precise than those stemming from the gravimetric method. A drift in TDR measurements is often observed, and it is for this reason that the calibration of probes by several gravimetric measurements is often recommended. In practice, the volumetric water content on a field scale is assumed to be equal to the mean value estimated from several representative samples collected from the top layer of soil.
Estimation of soil moisture by inversion of SAR data can be performed in using physical or semiempirical approaches. The physical approach uses backscattering models capable of reproducing the radar backscattering coefficient from the sensor configuration (wavelength, polarization, and incidence angle) and soil parameters (soil moisture and surface roughness for bare soils). However, several studies have reported a poor agreement between measured radar signals and those predicted by the models [18][19][20][21]. Differences between simulations and measurements may reach several decibels, rendering the inversion results inaccurate. The second approach consists of establishing experimental calibration relationships linking the signal radar to soil parameters (surface roughness and soil moisture of bare soils) and to sensor parameters (frequency, polarization, and incidence). This approach requires an important experimental database, acquired from several study sites and representative of possible physical soil conditions. The validity of these relationships for use in operational mode depends on the quality of the database and its representativeness. With current SAR sensors, the cartography in operational mode of the surface moisture of bare soils is carried out on images with one channel (one incidence and one polarization), which leads to the elaboration of experimental relationships between the signal radar and the soil moisture whatever the surface roughness [11,15,16]. Using this technology the precision on the prediction of soil surface moisture varies between 3 and 10 % of moisture depending on the incidence angle. Uncertainty analyses carried out with the Green and Ampt infiltration/runoff model showed that this degree of accuracy is adapted to the use of this kind of hydrological models [22].
Another approach for the estimation of soil moisture is to use a Soil-Vegetation-Atmosphere-Transfer model (SVAT). In this case, the SVAT is forced by analyzed or observed meteorological data to estimate mesoscale soil moisture over large areas. Such an approach is used within the SAFRAN-ISBA-MODCOU chain in France, which is composed of a meteorological analysis system (SAFRAN), a SVAT (ISBA) and a hydrological model (MODCOU, [23]). This system is able to simulate the water and energy budget at the surface, the soil moisture and the discharge of the main French rivers. It is run on an operational mode to evaluate the soil moisture at the scale of France. The typical scale of application of the SAFRAN-ISBA-MODCOU model is 8 km x 8 km over all of France. A comparison with satellite data is of great interest to validate the variables simulated by the system.
The objective of the present paper is to test the applicability of experimental relationships established between the C-band radar backscattering coefficient (ERS-2, RADARSAT-1, and ASAR sensors) and the soil moisture from a great experimental database collected between 1995 and 2005 at several study sites. Next, in the soil moisture mapping process, these relationships will be inverted and applied to ERS-2 and ASAR images, acquired in 2006 and 2007 over the Touch catchment basin (located near Toulouse in France) to provide additional key indicators to improve flood forecasting over the Touch catchment basin. The correlation between textural units and soil moisture will also be studied. Moreover, the soil moisture maps estimated from SAR images will be compared with simulations obtained by the ISBA soil-vegetation-atmosphere transfer model. This model, developed by Météo-France models the main surface processes and is currently coupled to atmospheric prediction model. As the soil moisture is prognostic variable of ISBA, it is thus necessary to quantify its quality. Soil moisture in situ data, measured by gravimetry, will be used to validate obtained results. Given that the presence of dense and high vegetation cover prevents C-band radar signal (wavelength ~ 6 cm) from reaching the ground, the soil moisture mapping from ERS and ASAR sensors is focused on bare soils or zones with little vegetation cover.

Database 1995-2005
An extensive experimental database collected at seven study sites was used for building relationships between radar backscattering coefficient and soil moisture. This database is composed of SAR images acquired between 1995 and 2005, and in situ soil moisture measurements (bare soils only) carried out jointly with radar acquisitions: •  1994, 1998, and 1999 to measure soil moisture and surface roughness. These sites are characterized by agricultural fields intended for growing wheat and corn. The measurements of volumetric soil moisture have been carried out on the first top 10 cm using two techniques: a TDR probe and gravimetry. The soil moisture ranges from 5.4% to 47.3%. SAR images used were acquired by ERS-2, RADARSAT-1 and ASAR, with a great range of incidence angles (20° to 48°), and at HH/HV/VV polarizations.
The database available for this study consists of about 400 doublets of backscattering coefficients and soil moisture. Each point of this database corresponds to the mean radar signal over a reference plot and the corresponding in situ soil moisture. The reference plots were chosen with a low local topography (generally flat).
Soil roughness measurements were also carried out, using a 1 or 2 m long profilometer with a 1 or 2 cm sampling interval (pinmeter or laser). For radar applications, the surface roughness of a given bare soil is defined statistically by two variables, determined from the surface height profiles: the standard deviation of surface height (rms), which specifies the vertical scale of the roughness, and the correlation length (L), representing the horizontal scale. The rms values vary in our database between 0.25 cm and 5.29 cm; the lowest ones correspond mainly to sown fields and the highest ones to ploughed fields. The correlation length ranges from 2 to 20 cm for 95% of measurements. Recent investigations have indicated that roughness variables estimated on the basis of field measurement data and simulations are very sensitive to profile length [e.g. [24][25][26]. The limited length of conventional profilers (one or two metres long) leads to large uncertainties in the estimated roughness variables. In a theoretical study based on simulated profiles, Oh and Kay [24] demonstrated that in order to measure rms surface height and correlation length with a precision of ±5% of their mean values, the length of the roughness profiles had to be at least 100L and 270L respectively. Measurement precision can be improved by averaging multiple profiles. For a correlation length that ranges from 2 to 20 cm, and for 10 averaged profiles, we deduce that the 2 m profiles provide a precision better than ±5% for rms surface height and between ±5% and ±15% for correlation length. At 1 m profile length, the precision can reach 10% for rms height and 20% for correlation length.

Study site
Our study site is the Touch catchment basin, located at the south-west of Toulouse in France (latitude: 43°17' N to 43°37' N, longitude: 0°48' E to 01°24' E; Figure 1). It covers a total area of about 500 km². The site is composed mainly of agricultural fields intended for growing wheat and corn. It is flat or slightly undulated with altitudes varying between 135 and 377 m. The Touch basin was selected as it is a site representative of the region where several catchments suffer from excessive runoff due to the conjunction of fragile soils and intensive agricultural practices, responsible for many damages. It is now selected as reference basin for several official authorities of flood forecasting and prevention, such as DIREN (Regional Direction of the Environment) and SCHAPI (National Flood Forecasting Agency).

Radar and optical images
Images used in this study were acquired by the ERS-2 and ASAR SARs between 08 March 2006 and 12 March 2007. Fifteen ASAR images and three ERS-2 images were obtained with incidence angles between 18° and 44°, and at HH/VV/VH polarizations (Table 1). Optical images obtained from SPOT and ASTER sensors were also used to allow a reliable mapping of bare soils. Before analyzing the SAR images, the data were radiometrically calibrated, which enables extraction of the backscattering coefficient (σ°) from the signal intensity of each pixel (BestW software, [27]). The radiometric calibration of SAR images allows converting the raw digital numbers (DN) to backscattering coefficients (σ°): where θ is the radar incidence angle, and K the calibration constant given in the image header. The backscattering coefficients are then expressed in decibels (dB): This calibration makes it possible to carry out multi-temporal analysis of different images (with same incidence and polarization). All SAR and optical images were then georeferenced using topographic maps with a root mean square error of the control points of about 30 m. The registration error of SAR images is taken into account by selecting areas of interest (AOI) within each training site and the border pixels were removed (defined by GPS control points).
Speckle noise, due to the coherent interference of waves reflected from many elementary scatterers, is present on SAR images and makes the pixel-by-pixel interpretation of SAR images extremely difficult. This explains why the estimation of soil characteristics is generally carried out on homogeneous sectors with several pixels or at field scale (which helps reducing speckle). In practice, the mean backscattering coefficients are calculated from calibrated SAR images by averaging the linear σ° values of all pixels within the field (or sub-field). A reduction in speckle and an improvement in the quality of our estimations are highly dependent on the size of the fields [e.g. 28,29]. Consequently, parcels with a homogeneous surface of 200 pixels or more are often used to study the behaviour of the signal as a function of SSC. In the case of ERS and ASAR sensors, this corresponds to fields of around two hectares or more.

ISBA data
The soil moistures simulated by the ISBA (Interactions between Soil Biosphere and Atmosphere) SVAT [30,31] were also used. The simulations were compared to the soil moistures estimated from SAR images and measured by gravimetry on March 12, 2007. ISBA is a SVAT designed to run in stand alone mode or in atmospheric models. The version used in the SAFRAN-ISBA-MODCOU is the force restore version of ISBA. The temperature evolution is described using two temperatures (surface temperature and deep soil). Three layers are used to describe the hydrological processes: the surface layer (one centimeter), the root zone layer (~1.50 m, depending of the vegetation type in the grid), the sub-root zone (usually down to 2 m below the surface). The hydrological processes which influence the liquid water content of the surface layer are precipitation (except the part intercepted by the vegetation) and evaporation from bare ground. In addition, it is relaxed to an equilibrium value which account for the balance between gravity and capillary forces (this value depends on the soil texture). The second layer accounts for the root depth zone. It is in relation with the surface layer via the relaxation term and with the sub-root zone via the gravitational drainage and the diffusion term (capillary versus gravitational drainage). The surface runoff occurs on the saturated part of the ISBA grid, which depends on the root zone wetness. It is calculated using the Variable infiltration capacity scheme [23]. Since the purpose of this study was to depict the state of surface moisture, only the surface level soil moisture was used in this study. For temperature, a force restore approach is used. The soil and vegetation parameters are derived from the ECOCLIMAP database [32] and from textural information (clay and sand fractions).
The surface water and energy budget were calculated by ISBA using the SAFRAN [33,34] meteorological data as input. SAFRAN used meteorological model outputs and meteorological observations to derive the variables used to force ISBA (near surface air temperature, wind and humidity, precipitation and incoming long wave and short wave radiation). The analysis is produced for climatologically homogeneous regions (approximately 30 km x 30 km), then interpolated on the 8 km x 8 km grid used to run ISBA. The data used in this paper are the data produced on an operational mode by Météo-France and used for soil moisture and drought monitoring.
It is important to underline that in this study the radar soil moisture estimates are obtained exclusively for bare soils or for soils with little vegetation, whereas ISBA simulations are carried out whatever the landuse (bare soils, forest, cultures …). A preliminary study has shown the absence of correlation between the ISBA moisture values and the percentage of bare soil in ISBA grids. Comparisons between radar moisture estimated on only bare soil and ISBA moisture calculated on 8 km x 8 km grids regardless of land use can therefore be carried out. In situ measurements Simultaneously with the acquisition of the SAR image of 12 March 2007, true ground measurements of soil moisture were performed on several test fields. Gravimetric moisture measurements on 33 test fields of bare soil well distributed on the whole of the catchment basin within ± 3 h of the SAR overpass were carried out. Gravimetric measurements were carried out on the upper 0-10cm soil layer because the radar signal penetration depth is only of few centimetres at C-band [e.g. 5,35]. The soil moisture on a field scale is assumed to be equal to the mean value estimated from 5 samples collected from the top layer of soil. The gravimetric soil moisture content is then transformed into volumetric moisture content by multiplying it by the bulk density of the soil. The volumetric soil moisture measurements range from 19.4% to 32.6% (standard deviation on the field scale about ±3%), and the soil bulk density ranges from 0.9 g/cm 3 to 1.6 g/cm 3 (standard deviation on the field scale about ±0.15 g/cm 3 ). During the ground campaign of 12 March 2007, the surface roughness of test fields varied from smooth (wheat and corn sown fields) to rough (ploughed fields). No roughness measurement was taken on March 12, 2007.
The lack of in situ moisture data encourages the analysis of radar and ISBA moistures using climatic and pedologic factors. For climatic factors, daily precipitation and temperature data recorded on six meteorological stations installed on the basin (Toulouse Blagnac, Toulouse Francazal, Fabas, La Bastidette, Lherm, and Rieumes) were used. Concerning pedologic factors, soil texture was used as its influence on soil moisture has already been demonstrated in numerous studies [36]. Other pedological factors are usually used to determine moisture, as organic carbon content, density, and depth of root zone [37,38]. These factors were not considered in this study because they show little variation within the basin.

Textural map
The textural map of the Touch basin, extracted from the Europe map at 1:1000000, has been thus used (Soil Geographical Database of Europe, SGBDE, [39]). This database contains information on soil texture ( Figure 2) and it was possible to define four textural units on the Touch basin. The Southwestern part of the basin is mainly composed of clayey material. In the center of the basin, the majority formation is the clayey-sandy soils. Silty soils are primarily present in the North-East of the catchment.

Relationship between radar signal and soil moisture
For bare soils, the radar backscattering coefficient in decibels can be written as the sum of two functions : the first one, f (linear), describes the dependence of radar signal on volumetric surface soil moisture, while the second, g (exponential), illustrates the dependence of σ° on surface roughness [13,15]: where k is the wave number (≈1.11 cm -1 for C-band), and rms is the root mean square surface height (characteristic parameter of surface roughness). For a given radar wavelength, the coefficients a, b, and c are observed to be dependent on both radar incidence angle and polarization [15,16,17]. In this equation, the correlation length is not taken into account. To retrieve soil moisture (m v ) from a single radar configuration, it is necessary to establish a relationship between the radar backscattering coefficient (σ°) and m v alone, without having any knowledge of the rms surface height. As a first approximation, the radar backscattering coefficient (in dB) may be expressed as follows [15,40]: This simplified relationship ignores the surface roughness. The coefficient a is dependent on both incidence angle and polarization. The coefficient d is primarily controlled by incidence angle, polarization and surface roughness. However, the coefficients describing the linear relationship between radar signal and soil moisture can be different from one catchment to another and also from one year to the next, and sometimes need to be calibrated. This difference is mainly due to the effects of roughness variation [41].
The coefficients a and d of the equation (4)  These coefficients were established for radar configurations in C-band, as a function of polarization and incidence angle, regardless of soil roughness. They are valid only for the same ranges of soil texture and surface roughness. The large database acquired over the last twelve years from study sites in France and Canada shows that the sensitivity of the radar signal to soil moisture is important for low and medium incidence angles (≤37°). The sensitivity of radar signal to soil moisture is slightly different with polarization (for low incidences, 0.23 dB/% in VV compared to 0.19dB/% in HH and 0.15 dB/% in HV). It is of the same order of magnitude for incidences between 20° and 37° (between 0.19 dB/% for HH20°-28° and 0.23 dB/% for HH 34°-37°) and seems to decrease for incidences greater than 39° (0.11 dB/%). These results therefore show that moisture mapping is optimal at low and medium incidences (≤37°). These results are comparable with the observed sensitivities by various studies [11,16,17,42]. Next, the inversion procedure was applied to the database of 2006-2007 for the estimation of soil moisture. The validity of this procedure was verified by comparing the inversion output with experimental data of soil moisture.

Methodology
Based on single SAR images (one polarization and one incidence), a simple procedure was used for mapping the surface soil moisture over bare soils [equation (4)]. The mapping of soil moisture was carried out on the Touch catchment using SAR data acquired in 2006 and 2007. For each SAR image, the relationship σ°d B (m v ) corresponding to the same radar configuration (incidence and polarization) was used (equation (4), Figure 3). All images were used, even those which correspond to low sensitivities of radar signal to soil moisture (18 March 2006 -IS7), except the VV-33° and VV-44° configurations because the relationships between σ° and m v could not be established for lack of data. Three soil moisture mapping scales estimates were carried out: an estimate for each ISBA cell (8 km x 8 km), an estimate for each textural unit, and an average estimate on the basin. For each SAR image, the bare soils map used corresponds to the one obtaining from the nearest optical image (SPOT and ASTER). The surface of bare soil covered between 8.7% (30 March 2006) and 24.1% (08 to 18 March 2006) of the basin's total area (cf. Table 1).
For mapping soil moisture on the ISBA grid (8 km x 8 km), the basin was divided into contiguous windows before calculating the mean backscattering coefficient on the whole of bare soil pixels of each cell. For the mapping according to textural units, the mean radar signal is calculated for each textural unit (Figure 2). The estimation of mean soil moisture on the scale of catchment basin uses the mean backscattering coefficient of all bare soils pixels present on the basin. Finally, soil moisture estimates are obtained by inverting the relationship between radar signal and the soil moisture [equation (4)]. The gravimetric soil moisture on ISBA grid scale is assumed to be equal to the mean value of all gravimetric measurements carried out in each grid. On the scale of catchment basin, the gravimetric soil moisture is obtained in averaging all in situ measurements.
The inversion performance of radar signal for retrieving soil moisture is evaluated to show the robustness of the presented technique. The ground truth measurements of soil moisture carried out the March 12, 2007 allow comparing effectively estimated soil moisture by SAR image and observed soil moisture. As suggested by Willmott [43], three statistical indexes were used: Mean Absolute Error: Root Mean Square Error: Index of agreement: Oi Pi (7) where P is the field or simulated variable, O the estimated variable, and N the data number. MAE and RMSE are expressed in %, while d is dimensionless.

Results and discussion
Comparison between radar estimates and gravimetric data The comparison between soil moistures estimated from SAR images of March 12, 2007 and gravimetric moistures shows a good agreement between estimated and measured soil moistures, with a bias about 2% and a RMSE of the order of 3% ( Figure 4, Table 2). VV and VH polarizations give similar performances. According to the radiometric accuracy of ASAR and ERS sensors, typically of ±1dB, the accuracy on the soil moisture estimation will be about ±5% (for a sensitivity of radar signal to soil moisture of 0.20 dB/%). This value is consistent with results presented in Table 2. The relationships between radar signal and soil moisture are then robust and can be used in the mapping process of soil moisture. Figure 5 shows the procedure applied for soil moisture mapping on the SAR image acquired the 12 March 2007 (VV polarization; incidence angle ~ 23°). The map of bare soils was recognized by SPOT 5 image acquired the 15 February 2007. The soil moisture estimated was next compared to in situ soil moisture measurements.

Analysis of ISBA simulations
ISBA moisture values were analyzed spatially according to the ISBA grid (Figures 6a and 6b). For very wet soil conditions (approximately 35%), the soil moistures simulated by ISBA show a maximum variation between ISBA cells of about 8%. For moistures about 25%, the ISBA moistures vary strongly between grids to reach 15%. For low moistures (approximately 15%), the maximum variation between grids is about 5% except for January 04, 2007 where a strong difference between some grids is observed (about 23%). Indeed, for January 04, 2007, ISBA simulates curiously low moistures on grids 20419, 20420 and 20421 (moistures 5% approximately whereas the other grids are about 20%). These comparisons are made at radar acquisition time (T_radar). Only ISBA moistures obtained by the average of ISBA simulations over all the day of January 04 showed moisture values of the same order for all grids. The three grids having an atypical behavior (20419 to 20421) are located in the SW part of the catchment basin.  (20419 to 20421). The analysis of ISBA input parameters for January 04, 2007 at 10 h (T_radar) shows that the soil was frozen. Freezing was much more marked on the three grids 20419 to 20421 than on the other grids because they are in another meteorological zone. By adding the frozen bulk soil water content, the total soil water content in the superficial reservoir becomes about 35% on the three grids 20419 to 20421. With the use of the total soil water content instead of the liquid soil water content only, the difference between radar and ISBA moistures becomes lower than 5%. (d) Figure 6b shows also that moistures simulated by ISBA are distributed according to two groups. Indeed, cells located at the North-East of the basin correspond to ISBA moistures lower than those of other cells. Moreover, moistures of each group are similar for a given date (variation inter group inferior to 5%). The comparison between ISBA simulations and textural map (Figure 2) shows that ISBA moistures are correlated with the soils texture, due to the fact that infiltration and runoff processes are driven in ISBA by the texture characteristics of the soil.

Comparison between ISBA and radar moistures
The radar moisture estimation was aggregated to produce moisture estimation for each ISBA cell. It shows that the maximum variations between grids are of 7% (18 March 2006, Figure 6c) for the radar, while they are around 15% with the ISBA model. For radar acquisitions dates (March 2006 to March 2007), the moisture estimated by radar on ISBA cells shows values ranging from 17% to 33%, whereas ISBA simulations varied between 10 and 45%. The comparison carried out between radar and ISBA soil moistures in using the same ISBA cells show that the difference between radar and ISBA moistures vary from -15% to +12%, except for January 04, 2007 where the difference reaches 23% (Figure 6d). Results show also that ISBA simulations averaged on the basin are well comparables to radar estimates with a mean difference lower than 5% for 12 dates among the 18 examined dates. For the six remaining dates, differences between 5% and 10% are observed (Figure 7).

Correlation analysis between precipitations and soil moistures (ISBA and radar)
The lack of in situ moisture data on the basin at the time of radar acquisitions (except for March 12, 2007) encouraged the analysis of radar moistures according to daily precipitations. Precipitations were given by six stations installed in the basin. Figure 7 shows the temporal variation of soil moisture estimated by SAR sensors between March 08, 2006 and March 12, 2007. This variation correlated with precipitations, which confirms that radar estimates are qualitatively coherent. For example, after important precipitations on March 10, 2006 (45 mm) the radar soil moisture increased from 28% on March 08, 2006 to 31% on March 11, 2006. From March 11, 2006 the radar soil moisture starts to decrease to reach 23% on March 18, 2006. This tendency is completely coherent because there no were precipitations between these two dates and the soil has thus dried. A good correlation between radar moisture and precipitations is also observed on the database of 2007. The radar moisture increased from 26% on January 01 to 28% on January 04, and decreases then to 27% on January 07, 2007. This behavior is coherent with precipitations about 13 mm between 01 and 04 January, and 2 mm between 04 and 07 January. The mean temperature during the first fortnight of January was of approximately 10°C.
ISBA simulations follow well the distribution of precipitations. Figure 7 shows ISBA moistures varying between 35 and 40% for the period of 08 to 11     Soil moisture mapping according to soil surface texture of the SGBDE Spatial variations of soil moisture obtained by gravimetry, radar images, and ISBA simulations were studied at textural zones scale. The textural zones ( Figure 2) were defined according to the textural information contained in the 1:1000000 Soil Geographical database. The gravimetric data, measured on 33 test fields, were classified according to their textural zones. Concerning radar data, the mean soil moisture was calculated for each date and each textural zone in using the mean radar signal of all bare soil pixels presents inside each zone (Figure 8a). For ISBA moistures, moisture values corresponding to bare soils were averaged on each textural zone (Figure 8b).
The behavior of the soil moisture calculated for textural zones is different for ISBA, radar, and gravimetric data. Figure 8a shows that there is no dependence between soil moistures estimated by radar images and the textural zones. However, an important dependence is obtained between ISBA soil moisture and textural zones (Figure 8b). Indeed, the ISBA soil moisture for zone 4, mainly composed of clayey material, is stronger than that for the other textural zones. The mean difference between radar estimates and ISBA simulations at textural class scale varies from -3.27% for zone 4 to +4.58 for zone 1 (bias=3.27% for zone 2 and 4.16 for zone 3). The standard deviation varies from 5.57% to 5.98% except for zone 4 where the value is 7.45% (Figure 8c).
Moreover, spatial distribution of moisture heterogeneities was displayed at textural zone scale, using box plots (Figure 9). For gravimetric database, 33 data points were used. For both radar and ISBA databases, 72 moisture values were used (18 dates x 4 textural zones). Results obtained between gravimetric or radar moistures do not show dependence with soil texture. On the contrary, ISBA data are well correlated to textural units. Figure 8. Correlation between textural zones and surface soil moisture simulated from ISBA model, and estimated from ASAR images. Each symbol corresponds to the mean soil moisture on each textural zone. Zone 1 is mainly sandy and zone 2 is mainly composed of clay (cf. Figure 2).   Concerning gravimetry and radar data, the differences between the main quartiles (median, first and third quartiles) of 4 textural zones do not exceed 2% (Figures 9a and 9b). Indeed, moisture would have had to be higher for zone 4 (mainly clay) and to decrease while going towards zone 1 (rather sandy).
As the ISBA model uses the pedological characteristics among input parameters, it would be evident that the soil moisture is dependent of soil texture. In fact, ISBA database shows that zone 4, mainly composed of clay, is well differentiable of other zones (1, 2, and 3). Indeed, the difference between the quartiles of zone 4 and other zones varies between 6% and 10% (Figure 9c). Results show also that pedological zones 1, 2 and 3 differ only of about 2%.
In order to consolidate these results, a nonparametric Kruskal Wallis test was carried out with a confidence threshold of 5%. This test confirms that for radar and gravimetric data, the difference between soil moistures of textural zones is not significant. For ISBA moistures, a significant dependence was observed between the textural zones and soil moisture. Radar estimations and gravimetric measurements of March 12, 2007 have shown that the differences of moisture observed between textural zones are not significant. The low spatial resolution of the textural map (1:1000000 scale) is not adapted to illustrate the moisture variations inside the basin. Indeed, in situ gravimetric data of 12 March 2007 show a moisture variation between the 33 test fields around 12% (moisture values between 21% and 33%), whereas at the textural zones scale the moisture variation is very weak (2%). It is then probable that the textural zones defined by the 1:1000000 map are not homogeneous and thus contain a mixture of textures. However, following the comparison between gravimetric moistures and the soil texture at the SGBDE textural zones scale, it was concluded that this scale is not adapted to explain the moisture variations on the basin. It would be thus more judicious to use in ISBA model a finer textural grid, at least on the scale of ISBA grids. Nevertheless, many studies concerning pedotransfer functions showed the importance of texture on the moisture estimation [36]. The hypothesis of dependence between moisture and soil texture is not excluded, but with the low spatial resolution of SGBDE map, this is not possible to quantify. In fact, textural spatial heterogeneities determining moisture values not appear at this scale. Thus, ISBA simulations based on textural criteria is legitimate but texture variations must be more precisely determined for obtain good moisture simulation.

Conclusions
This study examined the potential of ASAR and ERS-2 data for estimating volumetric soil moisture (m v ) over bare soils. Soil moisture estimates obtained from SAR acquisitions were compared to those obtained by the operational Météo-France ISBA model and by in situ measurements. These results appear promising for the development of simplified algorithms for retrieving soil moisture from SAR data, and for monitoring multi-temporal moisture changes. The mapping algorithm of soil moisture from SAR data performed very well. The results on the soil moisture estimates may be summarized as follows: • The comparison between radar soil moisture and in situ measurements shows a RMSE for the soil moisture estimate of the order of 3% (in comparison with gravimetric moistures of March 12, 2007). Other SAR images with a wide range of soil moisture (<20% and >35%) are required to better validate the linear regressions established between radar signal and soil moisture.
• The mean difference between ISBA simulations and radar estimates is lower than 5% for 12 dates among the 18 examined dates. The 6 remaining dates show a mean difference between 5% and 10%. The comparison between ISBA and gravimetric measurements of the 12 March 2007 shows a RMSE of about 6%. Generally, these results are very encouraging. • An important difference between ISBA and radar moistures is observed for some grids when the soil is frozen. The use of both liquid and frozen bulk soil water contents improves the results of the comparison between radar and ISBA model. • No dependence was observed between gravimetric or radar moistures and textural zones according to 1:1000000 textural map. The soil moisture mapping from SAR images at the textural scale shown a weak variation between textural units. The correlation observed between texture map and ISBA moisture is induced by the use of the texture map as an input parameter in the ISBA model. Even if this parameter is very important for soil moisture estimation, radar and gravimetric results shown that the scale at 1:1000000 of the used textural map is not appropriate to differentiate moistures zones. Finally, additional studies on the sensibility of soil moisture to each pedological parameter (texture, density, …) would be necessary in using higher resolution texture data.
Moreover, radar measurements are acquired with a given periodicity (some days) and cannot describe the continuous variations of moisture, whereas ISBA model provides moisture simulations every hour. Thus, it would be important to associate the radar remote sensing in the hydrological modelling. Assimilation of soil moisture data estimated from radar images could be used as calibration data to improve the forecast models (ISBA model for example).
The recent radar sensor (PALSAR/ALOS) and those that will come into operation early (TerraSAR-X, COSMO-SkyMed, RADARSAT-2) should provide us with polarimetric data (all polarizations), with the possibility of better characterising the soil surface. Polarimetry plays an important role as it allows a separation of soil moisture and surface roughness effects. Several inversion models are based on the use of fully polarimetric SAR images [14]. Polarimetric parameters such as, for example, the entropy, the α angle and the anisotropy should allow us to map two soil surface characteristics simultaneously [44]. Finally, the very high spatial resolution (metric) of future SAR sensors offers great potential in terms of improving the quality of SSC mapping. These new SAR sensors will provide a diagnosis suited to catchment areas where the parcels are of small size.
Currently it is difficult to map the moisture of soils covered by vegetation using C-band data. The arrival of L-band (ALOS) will enable to extend the soil moisture mapping from SAR data of bare soils only to soils with vegetation (depending on the density and height of the vegetation). Indeed, the penetration depth of radar wave in vegetable cover is more important at L-band than at C-band. Moreover, the use of L-band enables to map soils with vegetation without using optical images, which is not currently the case with C-band because an optical image is always used for this purpose. Dubois et al. [9] showed that the ratio of two polarizations HV/VV in L-band is strongly correlated with the NDVI (Normalized Difference Vegetation Index). If this information is confirmed, the cartography of soil moisture could be carried out directly from radar images without the use of optical data.