scholarly journals The Land Variational Ensemble Data Assimilation fRamework: LaVEnDAR

2019 ◽  
Author(s):  
Ewan Pinnington ◽  
Tristan Quaife ◽  
Amos Lawless ◽  
Karina Williams ◽  
Tim Arkebauer ◽  
...  

Abstract. The Land Variational Ensemble Data Assimilation fRamework (LaVEnDAR) implements the method of Four-Dimensional Ensemble Variational data assimilation for land surface models. Four-Dimensional Ensemble Variational data assimilation negates the often costly calculation of a model adjoint required by traditional variational techniques (such as 4DVar) for optimising parameters/state variables over a time window of observations. In this paper we implement LaVEnDAR with the JULES land surface model. We show the system can recover seven parameters controlling crop behaviour in a set of twin experiments. We run the same experiments at the Mead continuous maize FLUXNET site in Nebraska, USA to show the technique working with real data. We find that the system accurately captures observations of leaf area index, canopy height and gross primary productivity after assimilation and improves posterior estimates of the amount of harvestable material from the maize crop by 74 %. LaVEnDAR requires no modification to the model that it is being used with and is hence able to keep up to date with model releases more easily than other data assimilation methods.

2020 ◽  
Vol 13 (1) ◽  
pp. 55-69 ◽  
Author(s):  
Ewan Pinnington ◽  
Tristan Quaife ◽  
Amos Lawless ◽  
Karina Williams ◽  
Tim Arkebauer ◽  
...  

Abstract. The Land Variational Ensemble Data Assimilation Framework (LAVENDAR) implements the method of four-dimensional ensemble variational (4D-En-Var) data assimilation (DA) for land surface models. Four-dimensional ensemble variational data assimilation negates the often costly calculation of a model adjoint required by traditional variational techniques (such as 4D-Var) for optimizing parameters or state variables over a time window of observations. In this paper we present the first application of LAVENDAR, implementing the framework with the Joint UK Land Environment Simulator (JULES) land surface model. We show that the system can recover seven parameters controlling crop behaviour in a set of twin experiments. We run the same experiments at the Mead continuous maize FLUXNET site in Nebraska, USA, to show the technique working with real data. We find that the system accurately captures observations of leaf area index, canopy height and gross primary productivity after assimilation and improves posterior estimates of the amount of harvestable material from the maize crop by 74 %. LAVENDAR requires no modification to the model that it is being used with and is hence able to keep up to date with model releases more easily than other DA methods.


2021 ◽  
Author(s):  
Eduardo Emilio Sanchez-Leon ◽  
Natascha Brandhorst ◽  
Bastian Waldowski ◽  
Ching Pui Hung ◽  
Insa Neuweiler ◽  
...  

<p>The success of data assimilation systems strongly depends on the suitability of the generated ensembles. While in theory data assimilation should correct the states of an ensemble of models, especially if model parameters are included in the update, its effectiveness will depend on many factors, such as ensemble size, ensemble spread, and the proximity of the prior ensemble simulations to the data. In a previous study, we generated an ensemble-based data-assimilation framework to update model states and parameters of a coupled land surface-subsurface model. As simulation system we used the Terrestrial Systems Modeling Platform TerrSysMP, with the community land-surface model (CLM) coupled to the subsurface model Parflow. In this work, we used the previously generated ensemble to assess the effect of uncertain input forcings (i.e. precipitation), unknown subsurface parameterization, and/or plant physiology in data assimilation. The model domain covers a rectangular area of 1×5km<sup>2</sup>, with a uniform depth of 50m. The subsurface material is divided into four units, and the top soil layers consist of three different soil types with different vegetation. Streams are defined along three of the four boundaries of the domain. For data assimilation, we used the TerrsysMP PDAF framework. We defined a series of data assimilation experiments in which sources of uncertainty were considered individually, and all additional settings of the ensemble members matched those of the reference. To evaluate the effect of all sources of uncertainty combined, we designed an additional test in which the input forcings, subsurface parameters, and the leaf area index of the ensemble were all perturbed. In all these tests, the reference model had homogenous subsurface units and the same grid resolution as all models of the ensemble. We used point measurements of soil moisture in all data assimilation experiments. We concluded that precipitation dominates the dynamics of the simulations, and perturbing the precipitation fields for the ensemble have a major impact in the performance of the assimilation. Still, considerable improvements are observed compared to open-loop simulations. In contrast, the effect of variable plant physiology was minimal, with no visible improvement in relevant fluxes such as evapotranspiration. As expected, improved ensemble predictions are propagated longer in time when parameters are included in the update.</p>


2017 ◽  
Vol 21 (4) ◽  
pp. 2015-2033 ◽  
Author(s):  
David Fairbairn ◽  
Alina Lavinia Barbu ◽  
Adrien Napoly ◽  
Clément Albergel ◽  
Jean-François Mahfouf ◽  
...  

Abstract. This study evaluates the impact of assimilating surface soil moisture (SSM) and leaf area index (LAI) observations into a land surface model using the SAFRAN–ISBA–MODCOU (SIM) hydrological suite. SIM consists of three stages: (1) an atmospheric reanalysis (SAFRAN) over France, which forces (2) the three-layer ISBA land surface model, which then provides drainage and runoff inputs to (3) the MODCOU hydro-geological model. The drainage and runoff outputs from ISBA are validated by comparing the simulated river discharge from MODCOU with over 500 river-gauge observations over France and with a subset of stations with low-anthropogenic influence, over several years. This study makes use of the A-gs version of ISBA that allows for physiological processes. The atmospheric forcing for the ISBA-A-gs model underestimates direct shortwave and long-wave radiation by approximately 5 % averaged over France. The ISBA-A-gs model also substantially underestimates the grassland LAI compared with satellite retrievals during winter dormancy. These differences result in an underestimation (overestimation) of evapotranspiration (drainage and runoff). The excess runoff flowing into the rivers and aquifers contributes to an overestimation of the SIM river discharge. Two experiments attempted to resolve these problems: (i) a correction of the minimum LAI model parameter for grasslands and (ii) a bias-correction of the model radiative forcing. Two data assimilation experiments were also performed, which are designed to correct random errors in the initial conditions: (iii) the assimilation of LAI observations and (iv) the assimilation of SSM and LAI observations. The data assimilation for (iii) and (iv) was done with a simplified extended Kalman filter (SEKF), which uses finite differences in the observation operator Jacobians to relate the observations to the model variables. Experiments (i) and (ii) improved the median SIM Nash scores by about 9 % and 18 % respectively. Experiment (iii) reduced the LAI phase errors in ISBA-A-gs but had little impact on the discharge Nash efficiency of SIM. In contrast, experiment (iv) resulted in spurious increases in drainage and runoff, which degraded the median discharge Nash efficiency by about 7 %. The poor performance of the SEKF originates from the observation operator Jacobians. These Jacobians are dampened when the soil is saturated and when the vegetation is dormant, which leads to positive biases in drainage and/or runoff and to insufficient corrections during winter, respectively. Possible ways to improve the model are discussed, including a new multi-layer diffusion model and a more realistic response of photosynthesis to temperature in mountainous regions. The data assimilation should be advanced by accounting for model and forcing uncertainties.


2012 ◽  
Vol 25 (18) ◽  
pp. 6304-6317 ◽  
Author(s):  
Kevin Raeder ◽  
Jeffrey L. Anderson ◽  
Nancy Collins ◽  
Timothy J. Hoar ◽  
Jennifer E. Kay ◽  
...  

Abstract The Community Atmosphere Model (CAM) has been interfaced to the Data Assimilation Research Testbed (DART), a community facility for ensemble data assimilation. This provides a large set of data assimilation tools for climate model research and development. Aspects of the interface to the Community Earth System Model (CESM) software are discussed and a variety of applications are illustrated, ranging from model development to the production of long series of analyses. CAM output is compared directly to real observations from platforms ranging from radiosondes to global positioning system satellites. Such comparisons use the temporally and spatially heterogeneous analysis error estimates available from the ensemble to provide very specific forecast quality evaluations. The ability to start forecasts from analyses, which were generated by CAM on its native grid and have no foreign model bias, contributed to the detection of a code error involving Arctic sea ice and cloud cover. The potential of parameter estimation is discussed. A CAM ensemble reanalysis has been generated for more than 15 yr. Atmospheric forcings from the reanalysis were required as input to generate an ocean ensemble reanalysis that provided initial conditions for decadal prediction experiments. The software enables rapid experimentation with differing sets of observations and state variables, and the comparison of different models against identical real observations, as illustrated by a comparison of forecasts initialized by interpolated ECMWF analyses and by DART/CAM analyses.


2019 ◽  
Author(s):  
Xinxuan Zhang ◽  
Viviana Maggioni ◽  
Azbina Rahman ◽  
Paul Houser ◽  
Yuan Xue ◽  
...  

Abstract. Vegetation plays a fundamental role not only in the energy and carbon cycle, but also the global water balance by controlling surface evapotranspiration. Thus, accurately estimating vegetation-related variables has the potential to improve our understanding and estimation of the dynamic interactions between the water and carbon cycles. This study aims to assess to what extent a land surface model can be optimized through the assimilation of leaf area index (LAI) observations at the global scale. Two observing system simulation experiments (OSSEs) are performed to evaluate the efficiency of assimilating LAI through an Ensemble Kalman Filter (EnKF) to estimate LAI, evapotranspiration (ET), interception evaporation (CIE), canopy water storage (CWS), surface soil moisture (SSM), and terrestrial water storage (TWS). Results show that the LAI data assimilation framework effectively reduces errors in LAI simulations. LAI assimilation also improves the model estimates of all the water flux and storage variables considered in this study (ET, CIE, CWS, SSM, and TWS), even when the forcing precipitation is strongly positively biased (extremely wet condition). However, it tends to worsen some of the model estimated water-related variables (SSM and TWS) when the forcing precipitation is affected by a dry bias. This is attributed to the fact that the amount of water in the land surface model is conservative and the LAI assimilation introduces more vegetation, which requires more water than what available within the soil. Future work should investigate a multi-variate data assimilation system that concurrently merges both LAI and soil moisture (or TWS) observations.


2011 ◽  
Vol 8 (1) ◽  
pp. 1831-1877 ◽  
Author(s):  
A. L. Barbu ◽  
J.-C. Calvet ◽  
J.-F. Mahfouf ◽  
C. Albergel ◽  
S. Lafont

Abstract. The performance of the joint assimilation in a land surface model of a Soil Wetness Index (SWI) product provided by an exponential filter together with Leaf Area Index (LAI) is investigated. The data assimilation is evaluated with different setups using the SURFEX modeling platform, for a period of seven years (2001–2007), at the SMOSREX grassland site in southwestern France. The results obtained with a Simplified Extended Kalman Filter demonstrate the effectiveness of a joint data assimilation scheme when both SWI and Leaf Area Index are merged into the ISBA-A-gs land surface model. The assimilation of a retrieved Soil Wetness Index product presents several challenges that are investigated in this study. A significant improvement of around 13% of the root-zone soil water content is obtained by assimilating dimensionless root-zone SWI data. For comparison, the assimilation of in situ surface soil moisture is considered as well. A lower impact on the root zone is noticed. Under specific conditions, the transfer of the information from the surface to the root zone was found not accurate. Also, our results indicate that the assimilation of in situ LAI data may correct a number of deficiencies in the model, such as low LAI values in the senescence phase by using a seasonal-dependent error definition for background and observations. In order to verify the specification of the errors for SWI and LAI products, a posteriori diagnostics are employed. This approach highlights the importance of the assimilation design on the quality of the analysis. The impact of data assimilation scheme on CO2 fluxes is also quantified by using measurements of net CO2 fluxes gathered at the SMOSREX site from 2005 to 2007. An improvement of about 5% in terms of rms error is obtained.


2021 ◽  
Author(s):  
Sara Modanesi ◽  
Christian Massari ◽  
Alexander Gruber ◽  
Luca Brocca ◽  
Hans Lievens ◽  
...  

<p>Worldwide, the amount of water used for agricultural purposes is rising because of an increasing food demand. In this context, the detection and quantification of irrigation is crucial, but the availability of ground observations is limited. Therefore, an increasing number of studies are focusing on the use of models and satellite data to detect and quantify irrigation. For instance, the parameterization of irrigation in large scale Land Surface Models (LSM) is improving, but it is still characterized by simplifying assumptions, such as the lack of dynamic crop information, the extent of irrigated areas, and the mostly unknown timing and amount of irrigation. Remote sensing observations offer an opportunity to fill this gap as they are directly affected by, and hence potentially able to detect, irrigation. Therefore, combining models and satellite information through data assimilation can offer a viable way to quantify the water used for irrigation.</p><p>The aim of this study is to test how well modelled soil moisture and vegetation estimates from the Noah-MP LSM, with or without irrigation parameterization in the NASA Land Information System (LIS), are able to mimic in situ observations or to capture the signal of high-resolution Sentinel-1 backscatter observations in an irrigated area. The experiments were carried out over select sites in the Po river Valley, an important agricultural area in Northern Italy. To prepare for a data assimilation system, Level-1 Sentinel-1 backscatter observations, aggregated and sampled onto the 1 km EASE-v2 grid, were used to calibrate a Water Cloud Model (WCM) using simulated soil moisture and Leaf Area Index estimates. The WCM was calibrated with and without activating an irrigation scheme in Noah-MP. Results demonstrate that the use of the irrigation scheme provides the optimal calibration of the WCM, confirming the ability of Sentinel-1 to track the impact of human activities on the water cycle. Additionally, a first data assimilation experiment demonstrates the potential of Sentinel-1 backscatter observations to correct errors in Land Surface Model (LSM) simulations that are caused by unmodelled or wrongly modelled irrigation.</p>


2011 ◽  
Vol 8 (7) ◽  
pp. 1971-1986 ◽  
Author(s):  
A. L. Barbu ◽  
J.-C. Calvet ◽  
J.-F. Mahfouf ◽  
C. Albergel ◽  
S. Lafont

Abstract. The performance of the joint assimilation in a land surface model of a Soil Wetness Index (SWI) product provided by an exponential filter together with Leaf Area Index (LAI) is investigated. The data assimilation is evaluated with different setups using the SURFEX modeling platform, for a period of seven years (2001–2007), at the SMOSREX grassland site in southwestern France. The results obtained with a Simplified Extended Kalman Filter demonstrate the effectiveness of a joint data assimilation scheme when both SWI and Leaf Area Index are merged into the ISBA-A-gs land surface model. The assimilation of a retrieved Soil Wetness Index product presents several challenges that are investigated in this study. A significant improvement of around 13 % of the root-zone soil water content is obtained by assimilating dimensionless root-zone SWI data. For comparison, the assimilation of in situ surface soil moisture is considered as well. A lower impact on the root zone is noticed. Under specific conditions, the transfer of the information from the surface to the root zone was found not accurate. Also, our results indicate that the assimilation of in situ LAI data may correct a number of deficiencies in the model, such as low LAI values in the senescence phase by using a seasonal-dependent error definition for background and observations. In order to verify the specification of the errors for SWI and LAI products, a posteriori diagnostics are employed. This approach highlights the importance of the assimilation design on the quality of the analysis. The impact of data assimilation scheme on CO2 fluxes is also quantified by using measurements of net CO2 fluxes gathered at the SMOSREX site from 2005 to 2007. An improvement of about 5 % in terms of rms error is obtained.


2014 ◽  
Vol 21 (1) ◽  
pp. 303-323 ◽  
Author(s):  
N. Gustafsson ◽  
J. Bojarova ◽  
O. Vignes

Abstract. A hybrid variational ensemble data assimilation has been developed on top of the HIRLAM variational data assimilation. It provides the possibility of applying a flow-dependent background error covariance model during the data assimilation at the same time as full rank characteristics of the variational data assimilation are preserved. The hybrid formulation is based on an augmentation of the assimilation control variable with localised weights to be assigned to a set of ensemble member perturbations (deviations from the ensemble mean). The flow-dependency of the hybrid assimilation is demonstrated in single simulated observation impact studies and the improved performance of the hybrid assimilation in comparison with pure 3-dimensional variational as well as pure ensemble assimilation is also proven in real observation assimilation experiments. The performance of the hybrid assimilation is comparable to the performance of the 4-dimensional variational data assimilation. The sensitivity to various parameters of the hybrid assimilation scheme and the sensitivity to the applied ensemble generation techniques are also examined. In particular, the inclusion of ensemble perturbations with a lagged validity time has been examined with encouraging results.


Sign in / Sign up

Export Citation Format

Share Document