S Juglea

Y Kerr

A Mialon

J.-P Wigneron

E Lopez-Baeza

A Cano

A Albitar

C Millan-Scheiding

M Carmen

S Delwart

Modelling soil moisture at SMOS scale by use of a SVAT model over the Valencia Anchor Station

The main goal of the SMOS (Soil Moisture and Ocean Salinity) mission is to deliver global fields of surface soil moisture and sea surface salinity using L-band (1.4 GHz) radiometry. Within the context of the Science preparation for SMOS, the Valencia Anchor Station (VAS) experimental site, in Spain, was chosen to be one of the main test sites in Europe for Calibration/Validation (Cal/Val) activities. In this framework, the paper presents an approach consisting in accurately simulating a whole SMOS pixel by representing the spatial and temporal heterogeneity of the soil moisture fields over the wide VAS surface (50×50 km 2 ). Ground and meteorological measurements over the area are used as the input of a Soil-Vegetation-Atmosphere-Transfer (SVAT) model, SURFEX (Externalized Surface) -module ISBA (Interactions between Soil-Biosphere-Atmosphere) to simulate the spatial and temporal distribution of surface soil moisture. The calibration as well as the validation of the ISBA model are performed using in situ soil moisture measurements. It is shown that a good consistency is reached when point comparisons between simulated and in situ soil moisture measurements are made.

Actually, an important challenge in remote sensing approaches concerns product validation. In order to obtain an representative soil moisture mapping over the Valencia Anchor Station (50×50 km 2 area), a spatialization method is Correspondence to: S. Juglea (

Introduction

Soil moisture is a key variable controlling the exchanges of water and energy at the surface/atmosphere interface [START_REF] Betts | The land surface-atmosphere interaction : A review based on observational and global modeling perspectives[END_REF][START_REF] Entekhabi | Mutual Interaction of Soil Moisture state and atmospheric processes[END_REF]. It is highly variable both spatially and temporally as the result of the spatial heterogeneity of soil and vegetation properties, topography, land cover, rainfall and evapo-transpiration [START_REF] Bosch | Large scale measurements of soil moisture for validation of remotely sensed data: Georgia soil moisture experiment of 2003[END_REF][START_REF] Entekhabi | An analytic framework for the characterization of the space-time variability of soil moisture[END_REF]. Observing the spatial distribution of soil moisture at the catchment scale is a Published by Copernicus Publications on behalf of the European Geosciences Union.

difficult task requiring intensive field instrumentation for accurate spatial and temporal representation.

Nowadays, remote sensing technology has matured to the point that surface soil moisture can be estimated at global scale from space [START_REF] Wigneron | Retrieving near-surface soil moisture from microwave radiometric observations: current status and future plans[END_REF][START_REF] Wagner | Operational readiness of microwave remote sensing of soil moisture for hydrologic applications[END_REF]. Microwave remote sensing at low frequencies have been found to produce the best results [START_REF] Kerr | Soil Moisture from space: Where we are?[END_REF][START_REF] Wagner | Operational readiness of microwave remote sensing of soil moisture for hydrologic applications[END_REF][START_REF] Njoku | Passive microwave remote sensing of soil moisture[END_REF][START_REF] Jones | A microwave satellite observational operator for variational data assimilation of soil moisture[END_REF]. In spite of the importance of soil moisture observations, the instruments that have been or are currently operating are not adapted to soil moisture monitoring. Nevertheless, there are a number of soil moisture products available from different sensors. The Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) [START_REF] Njoku | Soil moisture retrieval from AMSR-E[END_REF] on board the National Aeronautics and Space Administration's (NASA) Aqua satellite and the scatterometers (SCAT) on board the European Remote Sensing Satellites 1 and 2 (ERS-1 and ERS-2) (Wagner et al., 1999a) provide soil moisture products. Both instruments use frequencies above 5 GHz.

The SMOS (Soil Moisture and Ocean Salinity) [START_REF] Kerr | Soil moisture retrieval from Space: The soil moisture and ocean salinity (SMOS) mission[END_REF] mission was designed to measure soil moisture over continental surfaces as well as ocean salinity using a low microwave frequency -L-band (1.4 GHz). At this frequency, microwave observations are sensitive to soil moisture through the effects of moisture (water) on the dielectric constant and hence on the emissivity of the soil. The soil emission is integrated over a soil depth of a few centimeters, giving a more representative measurement of soil moisture conditions over this layer. Consequently, the SMOS mission benchmark is to provide global maps of soil moisture with an accuracy better than 0.04 m 3 /m 3 [START_REF] Kerr | Soil moisture retrieval from Space: The soil moisture and ocean salinity (SMOS) mission[END_REF]. SMOS will achieve a maximum spatial resolution of 50 km over land (43 km on average over the field of view), providing multi-angular dual polarized (or fully polarized) brightness temperatures over the globe [START_REF] Kerr | Soil moisture retrieval from Space: The soil moisture and ocean salinity (SMOS) mission[END_REF]. Launched in November 2009, SMOS will deliver, for the first time, global surface soil moisture measurement twice a day (06:00 a.m. and 06:00 p.m. LT -local time) in less than 3 days.

L-band passive microwave radiometry is a very useful tool for soil moisture monitoring, allowing nearly all weather observation and surface vegetation cover information. Numerous field experiments using ground based and airborne Lband observations indicated a soil moisture retrieval capability of better than 0.04 m 3 /m 3 accuracy (Wang et al., 1990a;[START_REF] Schmugge | Passive microwave remote sensing of soil moisture: results from HAPEX, FIFE and MON-SOON'90[END_REF][START_REF] Jackson | Large area mapping of soil moisture using the ESTAR passive microwave radiometer[END_REF][START_REF] Jackson | Soil moisture Mapping at regional scales using microwave radiometry: The southern great plains hydrology experiment[END_REF]. In this context, the strategy adapted by ESA for its Soil Moisture and Ocean Salinity mission was to develop specific land product validation activities over well equipped monitoring sites. The Valencia Anchor Station (Lopez-Baeza et al., 2005a), in eastern Spain, and the Upper Danube Catchment [START_REF] Delwart | Overall SMOS Cal/Val Plan and Requirements, SMOS 7th workshop[END_REF], in southern Germany, are chosen as the two main test sites in Europe for the SMOS Calibration/Validation (Cal/Val) activities. This article will focus over the Valencia Anchor Station site which is a large reference area, equipped with ground soil moisture probes and fully characterized so as to contribute to SMOS land product validation.

Several papers evaluated the soil moisture remote sensing products [START_REF] Wagner | Soil moisture from operational meteorological satellites[END_REF][START_REF] Albergel | An evaluation of ASCAT surface soil moisture products with in-situ observations in Southwestern France[END_REF][START_REF] Draper | An evaluation of AMSR-E derived soil moisture over Australia[END_REF][START_REF] Rüdiger | An Intercomparison of ERS-SCAT and AMSR-E Soil Moisture Observations with Model Simulations over France[END_REF][START_REF] Gruhier | Soil moisture active and passive microwave products: intercomparison and evaluation over a Sahelian site[END_REF]. [START_REF] Draper | An evaluation of AMSR-E derived soil moisture over Australia[END_REF] provided a comparison of four soil moisture products all based on AMSR-E sensor over a temperate climate in Australia during 2006. [START_REF] Rüdiger | An Intercomparison of ERS-SCAT and AMSR-E Soil Moisture Observations with Model Simulations over France[END_REF] showed a comparison of several remotely sensed surface soil moisture products and one simulation (land surface model predictions) over the mainland of France from 2003 to 2005, in addition to a ground measurement comparison. [START_REF] Gruhier | Soil moisture active and passive microwave products: intercomparison and evaluation over a Sahelian site[END_REF] provided an inter-comparison and evaluation of five products derived from different active and passive microwaves sensors using local ground station measurements from three different ground sites over a Sahelian area (located in the Gourma-Mali region) during two consecutive years (2005)(2006).

Validating soil moisture products is a challenge and up to now, in most cases, papers describe how to relate one point measurement, or a value derived from a sparse network to a satellite product.

In the framework of SMOS Cal/Val activities, it was decided to select a set of areas scattered around the globe and representative of different types of ecoclimates. These sites are to deliver, continuously, a value representative of a whole pixel which can be compared to a satellite product at any overpass time for Cal/Val purposes. To achieve this goal it is necessary to characterise and monitor an area slightly larger than the actual pixel (3dB footprint) in terms of brightness temperature, so that it is possible to convolute the antenna pattern on it. To acquire such a large field of soil moisture, ground measurements are not tractable so we rely on a limited set of ground sites and spatialize the soil moisture information with use of a SVAT -ISBA (Interactions between Soil-Biosphere-Atmosphere; [START_REF] Noilhan | A simple parameterization of land surface fluxes processes for meteorological models[END_REF][START_REF] Noilhan | The ISBA land surface parameterization scheme[END_REF] coupled to a good set of forcings and a very good knowledge of soil types and land use. Once the soil moisture fields are known, it is possible to compute satellite level brightness temperatures (to check calibration for instance) or to compare to satellite products. As the model runs with a reasonably fine time step we can always have values at the time of overpass. To check the validity of the approach we did a test with existing sensors (AMSR-E, ERS-SCAT). The paper describes how such approach is validated over one such site, the Valencia Anchor Station (VAS). The idea is then to extend the approach on several other sites (arid, temperate, boreal etc.) as, such an approach being exhaustive, it can only be applied to a limited set of sites.

Characteristics of the area

The Valencia Anchor Station test site represents a reasonably homogeneous and mostly flat area of about 50×50 km 2 (Figs. 2 and3). The main cover type is vineyards, about 56%, followed by trees, shrubs, forest, industrial and urban. Beside the vineyard growing season, the area remains mostly under bare soil conditions. In spite of its relatively flat topography, the small altitude variations of the region clearly influence climate. It oscillates between semiarid in the areas of the towns of Utiel and Caudete de las Fuentes and dry-sub-humid towards Villagordo del Cabriel (about 16 km from Caudete de las Fuentes). Annual mean temperatures oscillate between 12

Available data over the area

To reproduce and compare the soil moisture fields over the VAS 50×50 km 2 area, in situ measurements and remotely sensed data products are used. The characteristics of these data are depicted next.

In situ measurements

Valencia Anchor Station is characterized by an extensive set of measurements at different levels (in the atmosphere and in the soil) in order to derive surface energy fluxes. Over the 50×50 km 2 area 22 meteorological stations are available (Table 1), 4 fully equipped and 18 rain gauges are not uniformly distributed (Fig. 1). Only the 4 fully equipped stations measures meteorological data: air temperature and humidity at screen level, atmospheric pressure, precipitation, wind speed and direction and solar and atmospheric radiation.

In the VAS area the soil texture is a parameter that depends mainly on lithology [START_REF] Lopez-Baeza | and team: Validation of SMOS Products over Mediterranean Ecosystem Vegetation at the Valencia Anchor Station Reference Area, Experimental Plan SMOS Validation Rehearsal Campaign[END_REF]. An accurate map representing the spatial distribution of clay and sand [START_REF] Millan-Scheiding | VAS Soil and Vegetation Characterization, SMOS meeting in Bordeaux[END_REF] at 10 m resolution covering all the 50×50 km 2 area is available (Fig. 3). The division of the texture is made in 23 main classes.

Leaf area index (LAI), roughness and fraction of vegetation are accessible during short time periods. However, as the period considered in this study is from 2004 to 2008, remote sensed LAI data are used (see Sect. 2.2.2.). The historical data for the roughness and the fraction of vegetation are www.hydrol-earth-syst-sci.net/14/831/2010/ Hydrol. Earth Syst. Sci., 14, 831-846, 2010 obtained from ECOCLIMAP (Fig. 2), a global land use maps database at 1 km resolution [START_REF] Masson | A global database of land surface parameters at 1km resolution in meteorological and climate models[END_REF].

Over the 50×50 km 2 area two major ground measurement campaigns took place. In order to characterize the shrubs and vineyards, in the framework of the Mediterranean Ecosystem L-Band characterization Experiment (Melbex 1 and Melbex 2), ground based L-band radiometry experiments have been developed to fully account for different soil moisture conditions and different vegetation growth development stages.

The first campaign, Melbex 1 (39.553 • N, 1.273 • W), was carried out between June 2005 and January 2006 to observe the surface emission of Mediterranean shrubs [START_REF] Cano | The SMOS Medierranean Ecosystem L-band experiment (Melbex-1) over natural shrubs, Remote Sens[END_REF]. The soil was characterized as sandy, with a soil texture composed of 47% sand, 38% silt and 15% clay. The vegetation is well adapted to dry conditions in summer and to freeze conditions in winter. The vegetation biomass is only subject to small variations throughout the year, it does not generally grow over a meter high and its distribution is random. Soil moisture measurements were carried out for the top first 5 cm of the soil, at 12 points every 10 min using capacitive probes. The ground soil moisture measurements were randomly scattered over the study area by placing probes both over bare soil and under shrubs. The probes were calibrated under laboratory conditions at the end of the experiment using the same soil type in order to correctly convert the raw voltage values into volumetric soil moisture content (m 3 /m 3 ).

The second soil campaign, Melbex 2 (39.526 • N, 1.288 • W), was carried out from April 2007 to December 2007 to observe the surface emission of vineyards [START_REF] Cano | The Mediterranean Ecosystem L-Band EXperiment over vineyards (Melbex-2), 10th Specialist Meeting on Microwave Radiometry and Remote Sensing for the Environment[END_REF]. The soil is characterized as sandy clay loam, with a texture composed of 45% sand, 29% silt and 26% clay. As in the previous experiment, soil moisture measurements were carried out at different representative points every 10 min using the same capacitive probes. In the area, the soil was ploughed at least 3 times during the growing period of vineyards.

Remote sensing data

Satellites data are used in this study. A short description of each of these data is given below.

AMSR-E data

The Advanced Microwave Scanning Radiometer (AMSR) of the Earth Observing System (EOS) is a multi-channel passive microwave instrument, launched on the Aqua satellite in May 2002. It operates in polar sun-synchronous orbit with equator crossing at 01:30 p.m. and 01:30 a.m. local solar time. Global coverage is achieved every two days or less depending on the latitude. The AMSR-E instrument www.hydrol-earth-syst-sci.net/14/831/2010/ Hydrol. Earth Syst. Sci., 14, 831-846, 2010 measures dual polarized radiation at six frequencies in the range of 6.9 to 89 GHz, with an incidence angle of 55 • . The mean spatial resolution at 6.9 GHz is about 56 km with a swath width of 1445 km.

In order to minimize the atmospheric effects and to maximize vegetation and soil penetration, the inversion algorithm for the retrieval of soil moisture was designed to use the Cband frequency in preference to the higher ones. However, due to the high level of RFI (Radio Frequency Interference) observed by AMSR-E at 6.9 GHz, the current AMSR-E soil moisture retrievals use only the 10.7 GHz and higher frequencies [START_REF] Njoku | Soil moisture retrieval from AMSR-E[END_REF].

The data used in this study are from the National Snow and Ice Data Center (NSIDC) Level 3 AMSR-E dataset [START_REF] Njoku | AMSR-E/AQUA daily L3 surface soil moisture, interpretive parms, & QC EASE-Grids[END_REF]. The daily averages of brightness temperature and soil moisture products are re-sampled to a global cylindrical 25 km Equal-Area Scalable Earth Grid (EASE-Grid) cell spacing [START_REF] Njoku | AMSR-E/AQUA daily L3 surface soil moisture, interpretive parms, & QC EASE-Grids[END_REF].

In addition to the soil moisture, the polarization ratio (PR) at 6.9 Ghz is used and is defined as:

PR = T b v -T b h T b v + T b h (1)
It normalizes out the surface temperature and leaves a quantity that depends primarily on soil moisture, vegetation and atmosphere (Kerr and Njoku, 1990;[START_REF] Njoku | Soil moisture retrieval from AMSR-E[END_REF]Owe et al., 2001). At low microwave frequencies, the polarization ratio has often been used to study soil moisture and vegetation effects. Its dynamic is well related to the soil moisture variations. At increasingly large angles (55 • in this case) there is a longer observation path through the vegetation layer, causing greater attenuation of the emission from the underlying soil and reducing the sensitivity to the soil moisture [START_REF] Njoku | Soil moisture retrieval from AMSR-E[END_REF].

Several studies investigated the validation and evaluation of AMSR-E soil moisture product [START_REF] Gruhier | Evaluation of AMSR-E soil moisture product based on ground measurements over temperate and semi-arid regions[END_REF][START_REF] Rüdiger | An Intercomparison of ERS-SCAT and AMSR-E Soil Moisture Observations with Model Simulations over France[END_REF][START_REF] Draper | An evaluation of AMSR-E derived soil moisture over Australia[END_REF]. As the AMSR-E soil moisture product shows biases and very small amplitude, a normalization between (0, 1) is done using:

y ′ = y -y min y max -y min (2)
where y ′ is the normalized curve and y is the input curve (in this case y is considered as the soil moisture product). Consequently, the discussion of this paper is focus on the normalized dataset.

ERS-SCAT data

The ERS (European Remote Sensing Satellites) scatterometer is an active low-resolution microwave sensor flown on the board of the ERS-1 and ERS-2 satellites. ERS-1 was launched in July 1991 followed by the identical ERS-2 in 1995. The first objective of this sensor is to measure wind over oceans, but its measurements have been shown to be highly suitable for surface soil moisture remote sensing [START_REF] Magagi | Retrieval of soil moisture and vegetation characteristics by use of ERS-1 wind scatterometer over arid and semi-arid areas[END_REF]Wagner et al., 1999a). The ERS scatterometer operates at 5.3 GHz (C-band), vertical polarization, collecting backscatter measurements over an incidence angle range from 18 • to 57 • . It operates in polar sun-synchronous orbit with equator crossing times at 10:30/22:30. The spatial resolution of the ERS-SCAT footprint is about 50 km with a 12.5 km spatial sampling interval.

The surface soil moisture data are retrieved from the radar backscattering coefficients, using the change detection method suggested by Dobson and Ulaby (1976). The methodology is described by Wagner et al. (1999a,b) which takes advantage of the information provided by the dual incidence angle measurements acquired by the ERS scatterometer. The backscattering coefficients are normalized to a reference incidence angle of 40 • . The relative soil moisture data ranging from 0% to 100% are derived by scaling the normalized backscattering coefficients between the lowest/highest values corresponding to the driest/wettest soil conditions.

MODIS data

The Moderate Resolution Imaging Spectroradiometer (MODIS; http://modis.gsfc.nasa.gov/) is an instrument on board of NASA's Terra and Aqua platforms. The MODIS LAI product is globally tiled and is projected on a sinusoidal grid which is an equivalent projection conserving the surface areas. It is at 1 km spatial resolution provided on a daily and 8-day basis and they are used as input to the SVAT model.

METEOSAT data

The METEOSAT radiometer is a geostationary weather satellite launched by the European Space Agency (ESA). In order to run the SVAT model and so to produce soil moisture maps, an atmospheric forcing is needed. One of the parameters of the atmospheric forcing is the shortwave radiation. As already mentioned, only 4 meteorological stations measuring shortwave radiation are available over the VAS area. For a better resolution over the entire area, the shortwave radiation flux products developed by Météo-France in the framework of the Satellite Application Facility on Land Surface Analysis (Land-SAF) are used (http://www.meteo.pt/landsaf/). The product is based on the 0.6 µm, 0.8 µm and 1.6 µm channels of the Meteosat/SEVIRI instrument and is calculated and distributed in near real time.

moisture. Spatially distributed fields and forcing enable to simulate soil moisture spatial and temporal behavior and thus averaged soil moisture at any moment for the whole pixel (50×50 km 2 ). The model used is SURFEX (stands for surface externalisée -Le [START_REF] Le Moigne | SURFEX scientific documentation[END_REF] and was developed at the National Center for Meteorological Research (CNRM) at Météo-France. It gathers all the developments and improvements made in surface schemes, containing four different modules: ISBA (Interactions between Soil-Biosphere-Atmosphere), Sea and ocean, TEB (Town Energy Balance) and Lake. In this article only the module for the soil and vegetation -ISBA [START_REF] Noilhan | A simple parameterization of land surface fluxes processes for meteorological models[END_REF]) is used. ISBA is a SVAT scheme which describes the exchanges of heat and water between the low-level atmosphere, the vegetation and the soil. It depends on the type of soil and of vegetation. It has been widely validated over vegetated and bare ground surfaces [START_REF] Mahfouf | Comparative study of various formulations of evaporation from bare soil using in situ data[END_REF][START_REF] Calvet | Retriving the rootzone soil moisture from surface soil moisture or temperature estimates : a feasibility study based on field measurements[END_REF].

The soil module in ISBA can be run in different configurations: 2-layers, 3-layers (ISBA-2L, ISBA-3L -with forcerestored discretization) and diffusive (ISBA-DIF). New possibilities comparing with ISBA-2 or 3 layers are available by using ISBA-DIF: the computation of a vertical profile of the temperature, liquid water and ice content over as many layers as needed. This scheme has already been applied successfully over a fallow site [START_REF] Boone | The influence of the inclusion of soil freezing on simulations by a soilvegetation-atmosphere transfer scheme[END_REF] and an agricultural site [START_REF] Boone | The inclusion of a third soil layer in a Land Surface Scheme using the Force-Restore method[END_REF]. In order to select the most appropriate configuration, different tests are made. A significant decrease in error is obtained in the case of a diffusive scheme so for our study the ISBA-DIF model is used.

Soil water transfer (infiltration, runoff, diffusion and drainage) in SVAT's is computed by equations which attempt to characterize the soil through a set of hydrological parameters. The ISBA scheme uses the [START_REF] Clapp | Empirical Equation for some soil hydraulic properties[END_REF] soil water model, which is common to a large number of surface parametrization schemes. The estimation of the diffusion of water in the soil is based on Darcy's law, where the water vertical flux is proportional to the gradient of the matric potential through the hydraulic conductivity.

F = -k ∂ ∂z ( + z) -D v ∂ ∂z -K d (3)
where D v is the vapor conductivity [START_REF] Braud | Bareground surface heat and water exchanges under dry conditions: Observations and parameterization[END_REF], is the soil water matric potential (m) , K d is an additional linear background drainage term (m s -1 ) and k is the hydraulic conductivity (m s -1 ). The hydraulic conductivity k (m s -1 ) and the soil water matric potential (m) are related to the liquid volumetric soil water content through [START_REF] Brooks | Properties of porous media affecting fluid flow[END_REF][START_REF] Clapp | Empirical Equation for some soil hydraulic properties[END_REF]:

k = k sat ( w l w sat ) 2b+3 (4) = sat ( w l w sat ) -b (5)
where b is the coefficient of the water retention curve.

SVAT configuration

In this section, the different sensitivity studies made as well as the parametrization chosen for the soil hydraulic functions are described. The characteristics of the data used for the calibration, validation and spatialization of SVAT model are also depicted.

The atmospheric forcing, needed to run the ISBA model, is composed of: air temperature and humidity at screen level, atmospheric pressure, precipitation, wind speed and direction and solar and atmospheric radiation. ISBA can have 12 patch types to characterize land use and related vegetation parametrization. For our case study, as the vegetation on the VAS site is mainly composed of vineyards, almonds trees (groves) and shrubs, the crops case is considered.

An important aspect is the soil layer discretization that enables one to compare realistic configurations as a function of the penetration depth, between ground measurements and/or the remote sensing data. A sensitivity study was conducted in order to test the influence of different parameters. The most representative configuration was chosen with 13 layers, with different thickness, from 1 cm at the surface down to 1.50 m of depth (1, 2, 3, 4, 5, 7, 9, 10, 30, 50, 80, 100, 150 cm).

Calibration of the SVAT

The calibration of the SVAT model is done to be applied over the entire test site for any season/year. In order to accurately perform the SVAT model calibration, ground measurements are needed. For this purpose, ground measurements from Melbex 1 campaign are used. The period considered goes from July to December 2005. Soil moisture measurements were carried out for the top first 5 cm of the soil every 10 min using capacitive probes. The atmospheric forcing is established as follows: the precipitation events are from a rain gauge (Caudete de las Fuentes) recording data at an hourly basis; the temperature, pression, wind speed, wind direction and relative humidity are from the nearest complete meteorological station; the shortwave is a Land-SAF product and the longwave is calculated using the formulation from [START_REF] Brutsaert | On a derivable formula for long-wave radiation from clear skies[END_REF]. In order to reproduce the exact condition from the Melbex 1 site, the same texture is considered, 47% sand, 38% silt and 15% clay. The LAI is from MODIS and both the fraction of vegetation and the roughness are from ECO-CLIMAP. The SVAT model is run on from 2004 until 2008 in a hourly basis.

The performance of the land-surface schemes and hence the soil moisture simulations are sensitive to the choice of soil hydraulic parameters [START_REF] Shao | Validation of soil moisture simulation in land surface parameterisation schemes with HAPEX data[END_REF]. Most of these hydrological parameters are site dependent. They are obtained from measurements or they are prescribed. It is difficult to prescribe a value for the wilting point (w wilt ), field capacity (w fc ), hydraulic conductivity at saturation (k sat ), saturated soil moisture (w sat ), the coefficient of the water retention curve (b) or for the matric www.hydrol-earth-syst-sci.net/14/831/2010/ Hydrol. Earth Syst. Sci., 14, 831-846, 2010 [START_REF] Giordani | Expériences de validation unidimensionnelles du schéma de surface NP89 aux normes Arpége sur trois sites de la campagne EFEDA 91[END_REF][START_REF] Noilhan | GCM Grid-Scale Evaporation from Mesoscale Modeling[END_REF] and calibrated (CALIBRATED ISBA (2) see [START_REF] Cosby | A Statistical exploration of the relationships of soil moisture characteristics to the physical properties of soils[END_REF] and [START_REF] Boone | The inclusion of a third soil layer in a Land Surface Scheme using the Force-Restore method[END_REF] potential at saturation ( sat ). To take into account the characteristics of the VAS area, the establishment of new databases for soil hydraulic parameters is necessary to improve the soil moisture simulations. The hydrological parameters of the soil (k sat , w sat , b, sat ) are calculated using empirical equations as a function of the percentages of sand and clay. Using ISBA's default relations [START_REF] Giordani | Expériences de validation unidimensionnelles du schéma de surface NP89 aux normes Arpége sur trois sites de la campagne EFEDA 91[END_REF][START_REF] Noilhan | GCM Grid-Scale Evaporation from Mesoscale Modeling[END_REF] to compute the soil parameters (see Table 2), the simulated soil moisture obtained is not in perfect accordance with the ground measurements recorded during the Melbex 1 campaign. In order to minimize this difference, a new set of equations (see Table 2) for the soil hydraulic parameters are established using [START_REF] Cosby | A Statistical exploration of the relationships of soil moisture characteristics to the physical properties of soils[END_REF]; [START_REF] Boone | The inclusion of a third soil layer in a Land Surface Scheme using the Force-Restore method[END_REF]. These calibrated equations are optimized inside the confidence interval defined in [START_REF] Cosby | A Statistical exploration of the relationships of soil moisture characteristics to the physical properties of soils[END_REF]; [START_REF] Boone | The inclusion of a third soil layer in a Land Surface Scheme using the Force-Restore method[END_REF]. Both sets of data, the one used by default by ISBA and the one from the calibrated version, are obtained from the same 11 textural classes and the same dataset. The results of the comparison between ground measurements and the simulated soil moisture using the new set of equations are given in Sect. 4.

Validation of the SVAT

In order to validate the calibration of the SVAT over another representative land use and other season/year, the Melbex 2 data are used. Carried out from April to December 2007 to observe the surface emission of vineyards, the soil moisture measurements were recorded for the first 5 cm of the soil every 10 min. Due the short distance between Melbex 1 and Melbex 2 campaigns sites (about 3 km), the same atmospheric forcing is used for both cases. The Caudete de las Fuentes rain gauge is situated about 3 km away from Melbex 2 site. However, the texture (45% sand, 29% silt and 26% clay), the LAI (MODIS), the fraction of vegetation and the roughness (ECOCLIMAP) are specific to the Melbex 2 site. In this case also the SVAT model was run on from 2004 until 2008 in a hourly basis (spin up).

SVAT distribution over the 50×50 km 2 area

The distribution of soil moisture patterns throughout a catchment plays a critical role in a variety of hydrological processes. Observing the spatial distribution of soil moisture at the catchment scale is a difficult task requiring intensive field instrumentation for an accurate spatial representation. A SVAT model driven with a fine resolution of meteorological forcing and land surface data can help understanding these processes. For this purpose, a trade off between the simulation time and the needs in spatial data was found by dividing the 50×50 km 2 area into 25 grid surfaces of 10×10 km 2 each (see Fig. 1). The available data sets over the area are in different formats and resolution so they had to be transformed so as to fit the gridded area. The way these scaling are performed is depicted in the next section.

Spatialization method

The 50×50 km 2 is divided into 25 areas of 10×10 km 2 each so as to better reproduce the high temporal and spatial heterogeneity of soil moisture fields over the entire VAS area. In situ measurement (soil moisture), detailed knowledge of the environment (land use, texture) and meteorological stations are used to characterize VAS. The format and the use of these data as an input to the SVAT model is presented in this section. The spatialized soil moisture obtained is compared with remote sensing data from AMSR-E and ERS-SCAT. To reduce the scaling issue, these remote sensing products are transformed as depicted next.

Interpolation

The atmospheric forcing is used as an input to the SVAT model to obtain the surface soil moisture. According to the dataset, in the 4 fully equipped meteorological stations located into the VAS 50×50 km 2 area, the measured data are registered on a 30/60 min basis. In addition, among the rain gauges, some of them are recording the weather information daily. In order to run the SVAT models with a suitable temporal resolution, standard diurnal cycles are reconstructed from the daily data.

Figure 1 presents the spatial distribution of the available meteorological station over the VAS 50×50 km 2 area. An irregular distribution of the stations can be noticed, especially in the center of the area where there is no data. So as to obtain a good representation of soil moisture over the entire area, an interpolation of all the available meteorological stations is necessary. In order to choose the most appropriate interpolation method between the inverse distance weighted (IDW) method and kriging, tests are done for different dates and for different meteorological stations/rain gauges. By doing a cross-validation analysis in general both techniques give the same behavior. The differences between the use of IDW or kriging are not significant so the choice of a sophisticated technique like the kriging is not justified. Moreover, the limited number of meteorological station/rain gauges as well as their location over the wide VAS area, was an additional reason for selecting the IDW method as the interpolation technique.

Inverse distance weighted methods are based on the assumption that the interpolated surface should be influenced mostly by the nearby points and less by the more distanced points. A general form of finding an interpolated value u for a given point x is an interpolated function:

u(x) = N k=0 w k (x) u k N k=0 w k (x) (6)
where the weight function is:

w k (x) = 1 d(x,x k ) (7)
defined by [START_REF] Shepard | A two-dimensional interpolation function for irregularly-spaced data[END_REF], x denotes an interpolated (arbitrary point), x k is the interpolated (known) point, d is a given distance from the known point x k to the unknown point x and N is the total number of known points. Using IDW method, the temperature, atmospheric pressure, wind speed, wind direction and the relative humidity are interpolated over the 10×10 km 2 grid by using the 4 complete meteorological stations. The shortwave fluxes are extracted over the same grid from the Land-SAF radiation product while the longwave fluxes are calculated using the interpolated data and the formulation from [START_REF] Brutsaert | On a derivable formula for long-wave radiation from clear skies[END_REF]. For the precipitation interpolation, all the 22 stations/rain gauges are considered. Following the interpolation, we have an optimal spatial and temporal distribution of the atmospheric forcing over the VAS 50×50 km 2 area.

Aggregation

The surface characteristics are also important to be considered in SVAT's input. The LAI (MODIS), the roughness and the fraction of vegetation (ECOCLIMAP) are 1 km resolution products. Due to their different spatial resolutions when compared to the 10×10 km 2 grid, these products are aggregated though a spatial mean. For the texture, as the maps available are at 10 m resolution, the aggregation to the 10×10 km 2 is done. In this case, the majority texture class is considered into the grid area.

The data obtained after spatialization thus the aggregated data used in order to simulate the spatialized soil moisture are depicted in Table 3. This allows to simulate the soil moisture over the chosen grid: in this case 25 points.

Mean

Once the soil moisture fields are known over the 10×10 km 2 grid, it is possible to compare to satellite products. To check both the approach presented in this paper and to validate all the aggregation techniques (several parameters are non linear but we rely on the fact that the overall variations are smooth and/or small enough to allow one to consider them as pseudo linear), we compared the spatialized soil moisture to existing products derived from either AMSR-E or from ERS-SCAT. The remote sensing products used are considered so as to correspond to the 50×50 km 2 area. The AMSR-E brightness temperature and soil moisture products are resampled to a global cylindrical 25 km Equal-Area Scalable Earth Grid (EASE-Grid) cell spacing [START_REF] Njoku | AMSR-E/AQUA daily L3 surface soil moisture, interpretive parms, & QC EASE-Grids[END_REF]. Two AMSR-E soil moisture sampled pixel are covering the VAS area. The average of these two pixels is considered to be representative for the 50×50 km 2 area. For the ERS-SCAT, the footprint is about 50 km with a 12.5 km spatial sampling interval. This gives 16 ERS-SCAT soil moisture products over the 50×50 km 2 area. In order to have a maximum temporal and spatial cover, the mean value of the 16 pixels is considered to be representative over the VAS area.

Results

An evaluation of the surface soil moisture obtained from ISBA was undertaken so as to quantify the improvement gained from the calibration. In a first step the differences when using the default or calibrated version of ISBA are shown by comparing with in situ measurements registered during the Melbex 1 campaign. Then, the calibrated version of ISBA is tested by comparing with data from Melbex 2 campaign (different place and different period). Point-like and 10×10 km 2 simulations are compared to ground measurements. In a second step, the calibrated version of ISBA as well as the spatialization method described in the previous section are used to obtain a spatial distribution of soil moisture over the entire area. The validity of the method was tested and the next paragraph presents a comparison between the spatialized soil moisture and remotely sensed data. As the area was divided into 25 pseudopixels (10×10 km 2 each), in order to have a representative value over the entire 50×50 km 2 , resulting values were averaged both spatially and over time. The soil moisture simulations were extracted for the time steps close to the overpass times of the satellites. The penetration depth was also taken into account, 2 cm for AMSR-E and for the ERS-SCAT product. Using the initial equations, the model tended to overestimate soil moisture in the dry season (from July to September) and to underestimate soil moisture for the rest of the period. In general a good agreement between the two datasets (correlation R 2 =0.793-/-, Nash efficiency Eff=0.619-/-) is observed, but the RMSE (Root Mean Squared Error) value equal to 0.042 m 3 /m 3 is higher than the SMOS requirements -0.04 m 3 /m 3 . In order to minimize this error, the SVAT model is calibrated (see Sect. 2.3). Table 4 presents the calculated soil hydraulic values using the default equations and the calibrated ones (Table 2).

Using the default equations (see Table 2), the SVAT model is not able to represent faithfully the dynamics of the first layer of soil during the dry season, keeping it at high soil moisture values (more than 0.10 m 3 /m 3 ). In order to enable lower values of soil moisture into the dry season, the wilting point equation is modified as follows. The value obtained for the wilting point using the default equation is about 0.140 m 3 /m 3 , whereas the minimum observed value of soil moisture during the campaign is 0.040 m 3 /m 3 . This behavior was also observed by [START_REF] Pellarin | Hydrological modelling and associated microwave emission of a semi-arid region in South-western Niger[END_REF], where a value of wilting point of 0.040 m 3 /m 3 was used for the simulation in Niger. The default SVAT equation is modified so as to retrieve a more representative value to simulate the dry period. The default equation for volumetric water content at saturation (w sat ) is also optimized. The goal is to reach a more representative value over a sandy soil as in Melbex 1 case (47% sand, 38% silt and 15% clay) at least for the first centimeters of soil.

The simulated soil moisture is driven mostly by the weather patterns and especially by the precipitation. The temporal evolution of surface soil moisture has a sharp increase and exponential decline which are caused, respectively by rainfall events and consecutive drying periods. By increasing the coefficient of the retention curve and reducing the hydraulic conductivity at saturation by the use of the calibrated equations, the drainage dynamics are slowed down in order to encounter the same behavior as for the measured soil moisture.

Using the calibrated version of ISBA a significant improvement is obtained for the modelled soil moisture at the first 5 cm (R 2 =0.908-/-, Eff=0.895-/--see Table 5). A good quantitative agreement is found (RMSE=0.022 m 3 /m 3 ) between the two soil moisture data: the same variability, the same drying slope, same low levels and amplitudes. At the beginning of November a higher level of modelled surface soil moisture dynamics is observed compared with in situ data. This can be due to the high value of w sat .

Validation of ISBA new parametrization using Melbex 2 campaign

The first step of our study was to find a parametrization of the surface model which minimizes the error compared with in situ measurements. In order to evaluate the validity of the chosen optimization, the same equations (see Table 2) as described in Sect. 2.3 are used for Melbex 2 area. Figure 5 presents a comparison between Melbex 2 data (in black) and simulated soil moisture (in red). A good agreement is retrieved between the two soil moisture data RMSE=0.024 m 3 /m 3 , R 2 =0.910-/-. Some differences can be observed on 10 August and also on 15 September. These differences can be mostly associated to the fact that the meteorological station used is situated at almost 3 km from the place were the campaign took place. The precipitations between the different location can easily differ showing, incidentally, that the spatial spatial distribution of rain is a key factor.

Ground measurements versus 10×10 km 2 soil moisture

In order to address the uncertainties associated with the spatial averaging, comparison between point like and 10×10 km 2 simulations against ground measurements are done (Fig. 6). The 10×10 km 2 area used is the one covering both Melbex 1 and Melbex 2 campaigns. or Melbex 2 period. When comparing in situ soil moisture with 10×10 km 2 simulations, similar results are observed.

All the scores obtained are in the required range (less than 0.04 m 3 /m 3 ). A slight overestimation of soil moisture within hight values is observed in all cases. This may be due to the volumetric water content at saturation which gives after calibration higher values than the recorded ones (see Sect. 4.1.1).

Using the SVAT for the spatialization as well as a linear interpolation method we can obtain for instance a soil moisture map over the entire 50×50 km 2 like in Fig. 7. The map represents the averaged soil moisture from January until the end of April 2005. A representation of soil moisture heterogeneity within the 50×50 km 2 is showed, with two marked zones of drier soil moisture. For this case a mean soil moisture value of 0.164 m 3 /m 3 is obtained.

Comparison with remote sensing data

Comparison with AMSR-E data

A comparison between spatialized soil moisture and the AMSR-E soil moisture product (Njoku L3) is made. The simulated soil moisture as well as the AMSR-E soil moisture product used are representative for the 50×50 km 2 area. The penetration depth of AMSR-E sensor is considered to be 2 cm so the soil moisture for the first two simulated layers is considered. The comparison is done from 2005 to 2007. Here are presented results from 2005 as the same evolution is observed for the other years. In a first step, the absolute values of the AMSR-E soil moisture product are compared with the simulated spatialized data. A severe lack of soil moisture dynamics and also a big difference between the absolute value of the two dataset are observed (RMSE=0.066 m 3 /m 3 , MBIAS=0.015 m 3 /m 3 , R 2 =0.051-/-). Because of the different soil moisture dynamics and biases, it is difficult to compare the various datasets in detail. Consequently, all next comparison are undertaken with normalized data, leading to the loss of the absolute aspects. The results of this normalization can be seen in Fig. 8. The dynamics of the soil moisture are very well captured at the beginning and also at the end of the year, from January to March we observed an RMSE=0.146-/-(R 2 =0.470-/-, MBIAS=-0.100-/-, Eff=-0.610-/-) and at the end of the year, from October to December RMSE=0.150-/-(R 2 =0.444-/-, MBIAS=0.041-/-, Eff=0.207-/-). In the winter season, the signal of AMSR-E soil moisture shows a small difference when compared to the spatialized soil moisture. This may be explained by the sensitivity of the microwave signal to soil freezing and by the reduced dynamics of the surface soil moisture during winter. In the middle of the year, from April to September, the opposite trend is observed between both datasets. From April to September no correlation is observed (R 2 =0.150-/-) and the RMSE is twice that of the rest of the year (0.348-/-). The inversion algorithm for the AMSR-E soil moisture uses the 10.7 GHz and 18.7 GHz brightness temperature data [START_REF] Njoku | Soil moisture retrieval from AMSR-E[END_REF]. The increased attenuation by vegetation and the superficial sensing depth (on the order of 1 cm for bare soil) for higher frequencies is a limit in the soil moisture retrieval from AMSR-E data. This can be seen by plotting the leaf area index (MODIS) corresponding to the site. When the growing season begins, the AMSR-E signal follows the vegetation signal (Fig. 8).

The polarization ratio provides a better agreement (than the soil moisture product from AMSR-E) with simulated soil moisture even in the vegetation growing period. This is shown by the scores obtained: at the beginning of the year, from January to March we observed an RMSE=0.165-/-(R 2 =0.655-/-, MBIAS=-0.142-/-, Eff=-0.872-/-), at the end of the year, from October to December RMSE=0.163-/-(R 2 =0.604-/-, MBIAS=-0.115-/-, Eff=-0.116-/-) and in the middle of the year, from April to September R 2 =0.466-/-Fig. 8. Comparison between surface soil moisture AMSR-E (black line), spatialized soil moisture from ISBA (red line) and the polarization ratio at 6.9 GHz AMSR-E (blue line). The soil moisture data are representative over 2 cm depth. All values are normalized between (0, 1). The leaf area index from MODIS is also represented here (green stars). and the RMSE=0.206-/-(MBIAS=-0.178-/-, Eff=-1.402-/-). If the entire year is considered, better scores are obtained: RMSE=0.186-/-, R 2 =0.627-/-, MBIAS=-0.153-/-, Eff=-0.553-/-. It shows that its dynamics is well correlated to the soil moisture variations.

Comparison with ERS-SCAT data

A comparison is also performed between spatialized soil moisture and the ERS-SCAT soil wetness product (Wagner SSM). The derived product represents the relative wetness of the first centimeters between totally dry conditions (0%) and total water capacity (100%).

In order to compare the two data sets, the absolute values of the simulated soil moisture are normalized between (0, 1). The penetration depth of ERS-SCAT sensor is considered to be at about 2 cm so the mean of the first two simulated layers from ISBA is considered. Figure 9 presents the comparison between simulated and observed surface soil moisture during a one year period. At the beginning of the year, from January to end of March RMSE=0.195-/-, R 2 =0.339-/-, MBIAS=-0.085-/-, Eff=-0.472-/-and at the end of the year, from the end of September to December RMSE=0.149-/-, R 2 =0.330-/-, MBIAS=0.036-/-, Eff=-0.735-/-. An underestimation of the soil moisture level by the ERS-SCAT product is observed. In the middle of the year, from April to September (R 2 =0.222-/-, RMSE=0.206-/-, MBIAS=0.102-/-, Eff=-1.014-/-), as for the AMSR-E soil moisture product, the vegetation influence the ERS-SCAT signal. This leads to an overestimation of the soil moisture estimates during the vegetation growing period.

For 2005 only 45 observations are available over the 50×50 km 2 area. This is a limit in completely understanding the soil moisture variability. A frequent revisit time is important for hydrologic applications, especially to obtain adequate sampling of surface wetting and drying between precipitation events.

Conclusions

In the framework of ESA's Soil Moisture and Ocean Salinity (SMOS) mission, this paper investigates the ability to reproduce the high temporal and spatial heterogeneity of soil moisture fields at SMOS pixel scale. This is the first step in the process of Calibration and Validation (Cal/Val) activities of the SMOS data. The study was performed for 2004-2008 over the Valencia Anchor Station (VAS), which was selected to be one of the main key test sites for the SMOS Cal/Val activities. In this purpose, in situ measurement including two main campaigns (Melbex 1 and Melbex 2), detailed knowledge of the environment (land use, texture) and meteorological stations are used to describe VAS. The hydrological process as well as the spatialization of the soil moisture fields are performed by the use of the SVAT model, SUR-FEX (Externalized Surface) -module ISBA (Interactions between Soil-Biosphere-Atmosphere) from Météo-France. In order to minimize the error with respect to SMOS mission requirements, both the calibration and the validation of the SVAT model were done using in situ soil moisture data from two campaigns: Melbex 1 and Melbex 2. Based on local atmospheric and surface observations from Valencia Anchor Station, it was found that the calibrated ISBA was able to faithfully reproduce the hydrological processes at the surface level. The RMSE decreases from 0.042 m 3 /m 3 when using the default ISBA to 0.022 m 3 /m 3 using the calibrated version. A new database for soil hydraulic parameters was established so as to improve the soil moisture simulations. The new parametrization was validated by the use of other dataset of soil moisture ground measurements recorded during Melbex 2 campaing. The value of 0.024 m 3 /m 3 obtained for the RMSE when comparing the simulated soil moisture with data from Melbex II campaign is adequate for assessing the SMOS validation with an accuracy better than 0.04 m 3 /m 3 . A comparison between point like and 10×10 km 2 simulations against ground measurements was also done to address the issue related to the spatial averaging. A slight overestimation of soil moisture is observed at high values but the scores obtained are within the required range (less than 0.04 m 3 /m 3 ).

The validation of soil moisture remote sensing products is an important issue and in most cases until now, papers describe how to associate point/network measurement to remote sensing products. None of these studies characterize the entire pixel as viewed by a satellite. The detailed knowledge of the area as well as the use of a interpolation method for the distribution of the atmospheric forcing allowed the obtaining of a spatial distribution of the soil moisture fields over the 50×50 km 2 area.

Comparisons of this spatialized soil moisture with three kinds of remote sensing information is done in order to test the approach. AMSR-E's soil moisture and polarization ratio as well as ERS-SCAT soil moisture products are used in this study. The penetration depth and the re-sampling grid of the soil moisture products used for each satellite are also considered. Although AMSR-E surface soil moisture product is not able to capture the absolute value, it provides reliable information on surface soil moisture temporal variability, at seasonal and rainy events scale. In the middle of the year, from April to September, the vegetation has an important influence on the measured signal. During the growing season the AMSR-E signal is very perturbed. The polarization ratio 6.9 GHz provides a better agreement with simulated soil moisture even in the vegetation growing period. A comparison is done between spatialized soil moisture and the ERS-SCAT soil moisture product (Wagner SSM). In this case also the seasons are well marked but the lack of a higher temporal resolution (45 observations were available over the area for 2005) prevented capturing the soil moisture variability over the VAS area.

Radiometer signals received in L and C band are susceptible to man-made radio frequency interferences (Le [START_REF] Vine | RFI at L-band in synthetic aperture radiometers[END_REF][START_REF] Njoku | Global survey and statistics of radio-frequency interference in AMSR-E land observations[END_REF]. SMOS, now launched, has detected a significant amount of RFI sources over the globe. Among these RFI sources, some of them are over Spain, affecting also the VAS test site. The Spanish authorities concentrate their efforts and managed already to stop most of these sources. The actual brightness temperature values over VAS are within the range of expected values. As the soil moisture simulation process is now validated, future works will imply simulating the SMOS brightness temperature as part of the Cal/Val activities.

Fig. 1 .

 1 Fig. 1. Location and distribution of the meteorological stations over the 50×50 km 2 VAS area.

Fig. 2 .Fig. 3 .

 23 Fig. 2. ECOCLIMAP land cover over the 50×50 km 2 .

  .92+0.049*CLAY*100+(100-CLAY*100-SAND*100)*0.014 b=b mean+b st dev sat mean =(1.54+(-0.010)*SAND*100+0.006*(100-CLAY*100-SAND*100)) sat (m) sat =-10 (1.85-0.88 * SAND) *0.01 sat st dev =(0.72+(-0.0026)*(100-CLAY*100-SAND*100)+0.001*CLAY*100) sat =-(10 sat mean -sat st dev /100) k sat mean =(-0.60+0.013*SAND*100+(-0.0064)*CLAY*100) k sat (m s -1 ) k sat =see ( * ) k sat st dev =(0.43+0.003*(100-CLAY*100-SAND*100)+0.001*CLAY*100) k sat =10 k sat mean -k sat st dev *(2.54/360000) w sat mean =(50.5+(-0.142) * SAND * 100+(-0.037) * CLAY * 100)/100 w sat (m 3 m -3 ) w sat =0.001* (-108*SAND+494.305) w sat st dev =(8.23+(-0.081)*CLAY*100+(-0.007)*SAND*100)/100 w 1sat =w sat mean +w sat st dev w wilt (m 3 m -3 ) w wilt =37.134 E-3*CLAY 0.5 w wilt =17.134 E-3*CLAY 0.5 w fc (m 3 m -3 ) w fc =89.047 E-3*CLAY 0.349 w fc =89.047 E-3*CLAY 0.349 * : 1.0e-6*10 (0.162E+01-0.582E+01*CLAY-0.907E-01*SAND+0.529E+01*CLAY 2 +0.120E+01*SAND 2 )

Fig. 4 .

 4 Fig. 4. Comparison between observed (black line) and simulated soil moisture using the default ISBA (clear grey line) and the calibrated ISBA (red line) at 5 cm depth. The precipitation corresponding to the meteorological station Caudete de las Fuentes are represented in blue.

4. 1

 1 Figure 4 compares the soil moisture from the Melbex 1 campaign and the point-like soil moisture simulated with ISBA using the parametrization described in Sect. 2.3 (default and the new set of equations). The precipitation recorded at the meteorological station Caudete de las Fuentes are represented in blue. The simulations are done for the period 2004-2008. Recorded soil moisture estimates are used as initial condition for the SVAT as well as a spin up of more than a year is done. For graphical convenience, only the period corresponding to Melbex 1 campaign (14 July-31 December 2005) is presented.Using the initial equations, the model tended to overestimate soil moisture in the dry season (from July to September) and to underestimate soil moisture for the rest of the period. In general a good agreement between the two datasets (correlation R 2 =0.793-/-, Nash efficiency Eff=0.619-/-) is observed, but the RMSE (Root Mean Squared Error) value equal to 0.042 m 3 /m 3 is higher than the SMOS requirements -0.04 m 3 /m 3 . In order to minimize this error, the SVAT model is calibrated (see Sect. 2.3). Table4presents the calculated soil hydraulic values using the default equations and the calibrated ones (Table2).Using the default equations (see Table2), the SVAT model is not able to represent faithfully the dynamics of the first layer of soil during the dry season, keeping it at high soil moisture values (more than 0.10 m 3 /m 3 ). In order to enable lower values of soil moisture into the dry season, the wilting point equation is modified as follows. The value obtained for the wilting point using the default equation is about 0.140 m 3 /m 3 , whereas the minimum observed value of soil moisture during the campaign is 0.040 m 3 /m 3 . This behavior

R

  2 (-/-) RMSE (m 3 m -3 ) MEAN BIAS (m 3 m -3 ) NASH EFFICIENCY (-

Fig. 5 .

 5 Fig. 5. Comparison between observed (black line) and simulated soil moisture using the calibrated ISBA (red line) at 5 cm depth. The precipitation corresponding to the meteorological station Caudete de las Fuentes are represented in blue.

Fig. 6 .

 6 Fig. 6. Comparison between point like (left side) and 10×10 km 2 simulations (right side) against ground measurements from Melbex 1 (up side) and Melbex 2 (down side) campaigns. The 10×10 km 2 area is the one covering both Melbex 1 and Melbex 2 campaigns and data are extracted to correspond to both campaign periods.

Fig. 7 .

 7 Fig. 7. Distributed soil moisture over the VAS 50×50 km 2 area. Data represents the mean of soil moisture from 1 January to 30 April 2005.

Fig. 9 .

 9 Fig. 9. Comparison between surface soil moisture ERS-SCAT (black stars) and spatialized soil moisture mean from ISBA (red line). The soil moisture data are representative over 2 cm depth.

Table 1 .

 1 Coordinates of the meteorological stations located in the 50×50 km 2 area.

	Station Name	Longitude Latitude	Characteristics
	VAS	1.288 • W	39.571 • N Fully equipped station
	CASAS DE VES	1.330 • W	39.262 • N Rain gauge
	CASAS IBANEZ	1.465 • W	39.288 • N Rain gauge
	VILLAMALEA	1.598 • W	39.363 • N Rain gauge
	REQUENA LA PORTERA COOP.	1.101 • W	39.405 • N Rain gauge
	REQUENA CAMPO ARCIS	1.165 • W	39.436 • N Rain gauge
	DEL MORO C. H. JUCAR	1.355 • W	39.484 • N Rain gauge
	REQUENA	1.096 • W	39.484 • N Rain gauge
	CAUDETE DE LAS FUENTES	1.317 • W	39.523 • N Rain gauge
	MINGLANILLA	1.595 • W	39.538 • N Rain gauge
	PRESA DE CONTRERAS	1.505 • W	39.542 • N Rain gauge
	UTIEL C. H. JUCAR	1.206 • W	39.568 • N Rain gauge
	UTIEL	1.205 • W	39.575 • N Rain gauge
	UTIEL (LA CUBERA -AUTOMATICA) 1.249 • W	39.580 • N Fully equipped station
	CAMPORROBLES COOPERATIVA	1.402 • W	39.649 • N Rain gauge
	CAMPO ARCIS	1.168 • W	39.433 • N Fully equipped station
	CERRITO REQUENA	1.107 • W	39.480 • N Fully equipped station
	GRAJA DE INIESTA	1.674 • W	39.516 • N Rain gauge
	CONTRERAS	1.498 • W	39.540 • N Rain gauge
	CAUDETE DE LAS FUENTES I	1.280 • W	39.547 • N Rain gauge
	VILLAMALEA I	1.602 • W	39.365 • N Rain gauge
	CERRO	1.512 • W	39.259 • N Rain gauge

Table 2 .

 2 Equation of hydrological parameters used in default (DEFAULT ISBA (1) see

Table 3 .

 3 Resolution and data used as input of the SVAT model in order to obtain the spatialized soil moisture.

	VARIABLE	INPUT sources	OUTPUT resolution
	LAI	MODIS -1 km resolution	aggregated to 10×10 km 2
	ROUGHNESS	ECOCLIMAP -1 km resolution	aggregated to 10×10 km 2
	FRACTION OF VEGETATION ECOCLIMAP -1 km resolution	aggregated to 10×10 km 2
	TEXTURE	clay and sand map at 10 m resolution	aggregated to 10×10 km 2
	TEMPERATURE	4 meteorological stations	interpolated at 10×10 km 2
	PRESSURE	4 meteorological stations	interpolated at 10×10 km 2
	WIND SPEED	4 meteorological stations	interpolated at 10×10 km 2
	WIND DIRECTION	4 meteorological stations	interpolated at 10×10 km 2
	RELATIVE HUMIDITY	4 meteorological stations	interpolated at 10×10 km 2
	SHORTWAVE fluxes	METEOSAT	extracted at 10×10 km 2
	LONGWAVE fluxes	calculated using interpolated atmospheric data
	PRECIPITATION	22 meteorological stations	interpolated at 10×10 km 2

Table 4 .

 4 Soil hydraulic values for Melbex 1 site calculated using the default and the calibrated set of equations. b (-/-) sat (m) k sat (m s -1 ) w sat (m 3 m -3 ) w wilt (m 3 m -3 ) w fc (m 3 m -3 )

	DEFAULT ISBA	5.556 -0.172	1.225E-05	0.444	0.144	0.230
	CALIBRATED ISBA 7.519 -0.049	1.502E-06	0.499	0.066	0.230

Table 5 .

 5 Statistics obtained by comparing the in situ measurements from Melbex 1 campaign with the default/calibrated simulated soil moisture using ISBA.

  The point like and 10×10 km 2 soil moisture data are extracted within Melbex 1 www.hydrol-earth-syst-sci.net/14/831/2010/ Hydrol. Earth Syst. Sci., 14, 831-846, 2010

Hydrol. Earth Syst. Sci., 14, 831-846, 2010www.hydrol-earth-syst-sci.net/14/831/2010/

The SVAT modelThe SVAT model is used to generate, from atmospheric forcing and initial conditions, the temporal behavior of the soil Hydrol. Earth Syst. Sci.,14,[831][832][833][834][835][836][837][838][839][840][841][842][843][844][845][846] 2010 www.hydrol-earth-syst-sci.net/14/831/2010/

www.hydrol-earth-syst-sci.net/14/831/2010/ Hydrol. Earth Syst. Sci., 14, 831-846, 2010

Acknowledgements. The authors wish to thank the European Space Agency (ESA), the Centre National d'Etudes Spatiales (CNES), the Centre National de la Recherche Scientifique -Institut National des Sciences de l'Univers (CNRS-INSU SIC) and the French National Programme TOSCA (Terre, Océans, Surfaces Continentales et Atmosphére) for supporting this work. We also wish to thank the NASA National Snow and Ice Data Center (NSIDC) for providing AMSR-E data as well as the Institute for Photogrammetry and Hydrol. Earth Syst. Sci., 14,[831][832][833][834][835][836][837][838][839][840][841][842][843][844][845][846] 2010 www.hydrol-earth-syst-sci.net/14/831/2010/

for providing the ERS-SCAT data. We thank also the Centre national de Recherches Météorologiques and Jean Christophe Calvet (CNRM) -Météo-France for the SURFEX model. The authors wish to thank also the Spanish Agency for Meteorology (AEMet) and to the Jucar River Basin Authority (CHJ) for the meteorological data.

Edited by: N. Verhoest

The publication of this article is financed by CNRS-INSU.