Valentijn R. N. Pauwels

The importance of the spatial patterns of remotely sensed soil moisture in the improvement of discharge predictions for small-scale basins through data assimilation

Abstract In this paper, we investigate to which degree information concerning the spatial patterns of remotely sensed soil moisture data are needed in order to improve discharge predictions from hydrological models. For this purpose, we use the TOPMODEL-based Land–Atmosphere Transfer Scheme (TOPLATS). The remotely sensed soil moisture values are determined using C-band backscatter data from the European Space Agency (ESA) European Remote Sensing (ERS) Satellites. A baseline run, without soil moisture assimilation, is established for both the distributed and lumped versions of the land–atmosphere scheme. The modeled discharge matches the observations slightly better for the distributed model than for the lumped model. The remotely sensed soil moisture data are assimilated into the distributed version of the model through the ‘nudging to individual observations’ method, and the ‘statistical correction assimilation’ method. The remotely sensed soil moisture data are also assimilated into the lumped version of the model through the ‘statistical correction assimilation’ method. The statistical correction assimilation method leads to similar, and improved, discharge predictions for both the distributed and lumped models. The nudging to individual observations method leads, for the distributed model, to only slightly better results than the statistical correction assimilation method. As a consequence, it is suggested that it is sufficient to assimilate the statistics (spatial mean and variance) of remotely sensed soil moisture data into lumped hydrological models when one wants to improve hydrological model-based discharge predictions.

Satellite-Scale Snow Water Equivalent Assimilation into a High-Resolution Land Surface Model

Four methods based on the ensemble Kalman filter (EnKF) are tested to assimilate coarse-scale (25 km) snow water equivalent (SWE) observations (typical of passive microwave satellite retrievals) into finescale (1 km) land model simulations. Synthetic coarse-scale observations are assimilated directly using an observation operator for mapping between the coarse and fine scales or, alternatively, after disaggregation (regridding) to the finescale model resolution prior to data assimilation. In either case, observations are assimilated either simultaneously or independently for each location. Results indicate that assimilating disaggregated finescale observations independently (method 1D-F1) is less efficient than assimilating a collection of neighboring disaggregated observations (method 3D-Fm). Direct assimilation of coarse-scale observations is superior to a priori disaggregation. Independent assimilation of individual coarse-scale observations (method 3D-C1) can bring the overall mean analyzed field close to the truth, but does not necessarily improve estimates of the finescale structure. There is a clear benefit to simultaneously assimilating multiple coarse-scale observations (method 3D-Cm) even as the entire domain is observed, indicating that underlying spatial error correlations can be exploited to improve SWE estimates. Method 3D-Cm avoids artificial transitions at the coarse observation pixel boundaries and can reduce the RMSE by 60% when compared to the open loop in this study.

State and bias estimation for soil moisture profiles by an ensemble Kalman filter: Effect of assimilation depth and frequency

[1] An ensemble Kalman filter for state estimation and a bias estimation algorithm were applied to estimate individual soil moisture profiles in a small corn field with the CLM2.0 model through the assimilation of measurements from capacitance probes. Both without and with inclusion of bias correction, the effect of the assimilation frequency, the assimilation depth, and the number of observations assimilated per profile were studied. Assimilation of complete profiles had the highest impact on deeper soil layers, and the optimal assimilation frequency was about 1–2 weeks, if bias correction was applied. The optimal assimilation depth depended on the calibration results. Assimilation in the surface layer had typically less impact than assimilation in other layers. Through bias correction the soil moisture estimate greatly improved. In general, the correct propagation of the innovations for both the bias-blind state and bias filtering from any layer to other layers was insufficient. The approximate estimation of the a priori (bias) error covariance and the choice of a zero-initialized persistent bias model made it impossible to estimate the bias in layers for which no observations were available.

Improvement of Modeled Soil Wetness Conditions and Turbulent Fluxes through the Assimilation of Observed Discharge

Abstract The objective of this paper is to improve the performance of a hydrologic model through the assimilation of observed discharge. Since an observation of discharge at a certain time is always influenced by the catchment wetness conditions and meteorology in the past, the assimilation method will have to modify both the past and present soil wetness conditions. For this purpose, a bias-corrected retrospective ensemble Kalman filter has been used as the assimilation algorithm. The assimilation methodology takes into account bias in the forecast state variables for the calculation of the optimal estimates. A set of twin experiments has been developed, in which it is attempted to correct the model results obtained with erroneous initial conditions and strongly over- and underestimated precipitation data. The results suggest that the assimilation of observed discharge can correct for erroneous model initial conditions. When the precipitation used to force the model is underestimated, the assimilation of...

Optimization of a coupled hydrology–crop growth model through the assimilation of observed soil moisture and leaf area index values using an ensemble Kalman filter

It is well known that the presence and development stage of vegetation largely influences the soil moisture content. In its turn, soil moisture availability is of major importance for the development of vegetation. The objective of this paper is to assess to what extent the results of a fully coupled hydrology-crop growth model can be optimized through the assimilation of observed leaf area index ( LAI) or soil moisture values. For this purpose the crop growth module of the World Food Studies ( WOFOST) model has been coupled to a fully process based water and energy balance model ( TOPMODEL-Based Land-Atmosphere Transfer Scheme ( TOPLATS)). LAI and soil moisture observations from 18 fields in the loamy region in the central part of Belgium have been used to thoroughly validate the coupled model. An observing system simulation experiment ( OSSE) has been performed in order to assess whether soil moisture and LAI observations with realistic uncertainties are useful for data assimilation purposes. Under realistic conditions ( biweekly observations with a noise level of 5 volumetric percent for soil moisture and 0.5 for LAI) an improvement in the model results can be expected. The results show that the modeled LAI values are not sensitive to the assimilation of soil moisture values before the initiation of crop growth. Also, the modeled soil moisture profile does not necessarily improve through the assimilation of LAI values during the growing season. In order to improve both the vegetation and soil moisture state of the model, observations of both variables need to be assimilated.

A soil‐vegetation‐atmosphere transfer scheme for the modeling of water and energy balance processes in high latitudes: 1. Model improvements

This paper describes the development of a process-based water and energy balance model for use in high latitudes under both summertime and wintertime conditions. The model is developed as a part of the Boreal Ecosystem-Atmosphere Study (BOREAS). The model differs from its original version in the representation of hydrological processes specific to climatic and ecologic conditions in high latitudes: the impact of an organic layer on the water and energy balance of the boreal forest; the parameterization of frozen ground and snow accumulation and ablation; and the effect of open water bodies on the calculation of the radiative energy and water budget. The model can be run either in a fully distributed or an aggregated statistical mode. The resulting macroscale hydrological model can be used, off line, to study the energy and water balance of the boreal forest and potentially can be used as a land-atmosphere parameterization in atmospheric models.

Global Calibration of the GEOS-5 L-Band Microwave Radiative Transfer Model over Nonfrozen Land Using SMOS Observations

AbstractA zero-order (tau-omega) microwave radiative transfer model (RTM) is coupled to the Goddard Earth Observing System, version 5 (GEOS-5) catchment land surface model in preparation for the future assimilation of global brightness temperatures (Tb) from the L-band (1.4 GHz) Soil Moisture Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) missions. Simulations using literature values for the RTM parameters result in Tb biases of 10–50 K against SMOS observations. Multiangular SMOS observations during nonfrozen conditions from 1 July 2011 to 1 July 2012 are used to calibrate parameters related to the microwave roughness h, vegetation opacity τ and/or scattering albedo ω separately for each observed 36-km land grid cell. A particle swarm optimization is used to minimize differences in the long-term (climatological) mean values and standard deviations between SMOS observations and simulations, without attempting to reduce the shorter-term (seasonal to daily) errors. After calibration, global T...

Crop Classification Using Short-Revisit Multitemporal SAR Data

Classification of crops and other land cover types is an important application of both optical/infrared and SAR satellite data. It is already an import application of present satellite systems, as it will be for planned missions, such as the Sentinels. An airborne SAR data set with a short revisit time acquired by the German ESAR system during the ESA-campaign, AgriSAR 2006, has been used to assess the performance of different polarization modes for crop classification. Both C-and L-band SAR data were acquired over the Demmin agricultural test site in North Eastern Germany on a weekly basis during the growing season. Single-and dual-polarization, and fully polarimetric data have been used in the analysis (fully polarimetric data were only available at L-band). The main results of the analysis are, that multitemporal acquisitions are very important for single-and dual-polarization modes, and that cross-polarized backscatter produces the best results, with errors down to 3%-6% at the two frequencies. There is a trade-off between the polarimetric information and the multitemporal information, where the best overall results are obtained using the multitemporal information. If only a few acquisitions are available, the polarimetric mode may perform better than the single-and dual polarization modes.

Optimization of Soil Hydraulic Model Parameters Using Synthetic Aperture Radar Data: An Integrated Multidisciplinary Approach

It is widely recognized that synthetic aperture radar (SAR) data are a very valuable source of information for the modeling of the interactions between the land surface and the atmosphere. During the last couple of decades, most of the research on the use of SAR data in hydrologic applications has been focused on the retrieval of land and biogeophysical parameters (e.g., soil moisture contents). One relatively unexplored issue consists of the optimization of soil hydraulic model parameters, such as, for example, hydraulic conductivity values, through remote sensing. This is due to the fact that no direct relationships between the remote-sensing observations, more specifically radar backscatter values, and the parameter values can be derived. However, land surface models can provide these relationships. The objective of this paper is to retrieve a number of soil physical model parameters through a combination of remote sensing and land surface modeling. Spatially distributed and multitemporal SAR-based soil moisture maps are the basis of the study. The surface soil moisture values are used in a parameter estimation procedure based on the extended Kalman filter equations. In fact, the land surface model is, thus, used to determine the relationship between the soil physical parameters and the remote-sensing data. An analysis is then performed, relating the retrieved soil parameters to the soil texture data available over the study area. The results of the study show that there is a potential to retrieve soil physical model parameters through a combination of land surface modeling and remote sensing.

A soil-vegetation-atmosphere transfer scheme for the modeling of water and energy balance processes in high latitudes: 2. Application and validation

In support of the overall scientific objective of the Boreal Ecosystem-Atmosphere Study (BOREAS), which encompasses the improved understanding of the interactions between the boreal forest and the atmosphere, a process-based water and energy balance model is applied to observed forcing data, and the results are presented and discussed. Observed tower forcing and validation data are analyzed. A consistent diurnal pattern in the energy balance closure of the validation data is obtained. Simulations are performed for a number of BOREAS flux tower sites. The model successfully simulates the temporally averaged Bowen ratio and the evaporative part of precipitation over the different BOREAS flux tower sites during the 1994 and 1996 intensive field campaigns. At finer temporal scales a small phase shift in sensible heat flux and net radiation exists between the observed and model-derived quantities. The ground heat flux is found to be slightly larger than the observations during the course of the day. It is suggested that the sensitivity of the model to parameters such as the moss thickness, thermal conductivity, and heat capacity is responsible for these differences. The moss moisture content and the different components of the energy balance were very well matched for a continuous simulation during 1996. Overall, the accuracy performance of the model is equivalent to the accuracy of the input forcing data.