Multivariate and multiscale data assimilation in terrestrial systems: a review.

More and more terrestrial observational networks are being established to monitor climatic, hydrological and land-use changes in different regions of the World. In these networks, time series of states and fluxes are recorded in an automated manner, often with a high temporal resolution. These data are important for the understanding of water, energy, and/or matter fluxes, as well as their biological and physical drivers and interactions with and within the terrestrial system. Similarly, the number and accuracy of variables, which can be observed by spaceborne sensors, are increasing. Data assimilation (DA) methods utilize these observations in terrestrial models in order to increase process knowledge as well as to improve forecasts for the system being studied. The widely implemented automation in observing environmental states and fluxes makes an operational computation more and more feasible, and it opens the perspective of short-time forecasts of the state of terrestrial systems. In this paper, we review the state of the art with respect to DA focusing on the joint assimilation of observational data precedents from different spatial scales and different data types. An introduction is given to different DA methods, such as the Ensemble Kalman Filter (EnKF), Particle Filter (PF) and variational methods (3/4D-VAR). In this review, we distinguish between four major DA approaches: (1) univariate single-scale DA (UVSS), which is the approach used in the majority of published DA applications, (2) univariate multiscale DA (UVMS) referring to a methodology which acknowledges that at least some of the assimilated data are measured at a different scale than the computational grid scale, (3) multivariate single-scale DA (MVSS) dealing with the assimilation of at least two different data types, and (4) combined multivariate multiscale DA (MVMS). Finally, we conclude with a discussion on the advantages and disadvantages of the assimilation of multiple data types in a simulation model. Existing approaches can be used to simultaneously update several model states and model parameters if applicable. In other words, the basic principles for multivariate data assimilation are already available. We argue that a better understanding of the measurement errors for different observation types, improved estimates of observation bias and improved multiscale assimilation methods for data which scale nonlinearly is important to properly weight them in multiscale multivariate data assimilation. In this context, improved cross-validation of different data types, and increased ground truth verification of remote sensing products are required.

Salix psammophila afforestations can cause a decline of the water table, prevent groundwater recharge and reduce effective infiltration.

Afforestation can reduce desertification and soil erosion. However, the hydrologic implications of afforestation are not well investigated, especially in arid and semi-arid regions. China has the largest area of afforestation in the world, with one-third of the world's total plantation forests. How the shrubs affect evapotranspiration, soil moisture dynamics, and groundwater recharge remains unclear. We designed two pairs of lysimeters, one being 1.2 m deep and the other one 4.2 m deep. Each pair consists of one lysimeter with bare soil, while on the other one a shrub is planted. The different water table depths were implemented to understand how depth to groundwater affects soil moisture and water table dynamics under different hydrological conditions. Soil moisture, water table depth, sap flow, and rainfall were measured concurrently. Our study confirms that for the current meteorological conditions in the Ordos plateau recharge is reduced or even prohibited through the large-scale plantation Salix psammophila. Shrubs also raise the threshold of precipitation required to increase soil moisture of the surface ground. For the conditions we analyzed, a minimum of 6 mm of precipitation was required for infiltration processes to commence. In addition to the hydrological analysis, the density of root distribution is assessed outside of the lysimeters for different water table depths. The results suggest that the root-density distribution is strongly affected by water table depth. Our results have important implications for the determination of the optimal shrub-density in future plantations, as well as for the conceptualization of plant roots in upcoming numerical models.

Real-time management of an urban groundwater well field threatened by pollution.

We present an optimal real-time control approach for the management of drinking water well fields. The methodology is applied to the Hardhof field in the city of Zurich, Switzerland, which is threatened by diffuse pollution. The risk of attracting pollutants is higher if the pumping rate is increased and can be reduced by increasing artificial recharge (AR) or by adaptive allocation of the AR. The method was first tested in offline simulations with a three-dimensional finite element variably saturated subsurface flow model for the period January 2004-August 2005. The simulations revealed that (1) optimal control results were more effective than the historical control results and (2) the spatial distribution of AR should be different from the historical one. Next, the methodology was extended to a real-time control method based on the Ensemble Kalman Filter method, using 87 online groundwater head measurements, and tested at the site. The real-time control of the well field resulted in a decrease of the electrical conductivity of the water at critical measurement points which indicates a reduced inflow of water originating from contaminated sites. It can be concluded that the simulation and the application confirm the feasibility of the real-time control concept.

A 3 km spatially and temporally consistent European daily soil moisture reanalysis from 2000 to 2015.

High-resolution soil moisture (SM) information is essential to many regional applications in hydrological and climate sciences. Many global estimates of surface SM are provided by satellite sensors, but at coarse spatial resolutions (lower than 25 km), which are not suitable for regional hydrologic and agriculture applications. Here we present a 16 years (2000-2015) high-resolution spatially and temporally consistent surface soil moisture reanalysis (ESSMRA) dataset (3 km, daily) over Europe from a land surface data assimilation system. Coarse-resolution satellite derived soil moisture data were assimilated into the community land model (CLM3.5) using an ensemble Kalman filter scheme, producing a 3 km daily soil moisture reanalysis dataset. Validation against 112 in-situ soil moisture observations over Europe shows that ESSMRA captures the daily, inter-annual, intra-seasonal patterns well with RMSE varying from 0.04 to 0.06 m3m-3 and correlation values above 0.5 over 70% of stations. The dataset presented here provides long-term daily surface soil moisture at a high spatiotemporal resolution and will be beneficial for many hydrological applications over regional and continental scales.