scholarly journals Parameter optimisation for a better representation of drought by LSMs: inverse modelling vs. sequential data assimilation

2017 ◽  
Vol 21 (9) ◽  
pp. 4861-4878 ◽  
Author(s):  
Hélène Dewaele ◽  
Simon Munier ◽  
Clément Albergel ◽  
Carole Planque ◽  
Nabil Laanaia ◽  
...  
Keyword(s):  
Data Assimilation ◽  
Land Surface ◽  
Inverse Modelling ◽  
Complex Method ◽  
Sequential Data ◽  
Annual Maximum ◽  
Area Index ◽  
Available Water Content ◽  
Cereal Grain ◽  
Administrative Units

Abstract. Soil maximum available water content (MaxAWC) is a key parameter in land surface models (LSMs). However, being difficult to measure, this parameter is usually uncertain. This study assesses the feasibility of using a 15-year (1999–2013) time series of satellite-derived low-resolution observations of leaf area index (LAI) to estimate MaxAWC for rainfed croplands over France. LAI interannual variability is simulated using the CO2-responsive version of the Interactions between Soil, Biosphere and Atmosphere (ISBA) LSM for various values of MaxAWC. Optimal value is then selected by using (1) a simple inverse modelling technique, comparing simulated and observed LAI and (2) a more complex method consisting in integrating observed LAI in ISBA through a land data assimilation system (LDAS) and minimising LAI analysis increments. The evaluation of the MaxAWC estimates from both methods is done using simulated annual maximum above-ground biomass (Bag) and straw cereal grain yield (GY) values from the Agreste French agricultural statistics portal, for 45 administrative units presenting a high proportion of straw cereals. Significant correlations (p value  <  0.01) between Bag and GY are found for up to 36 and 53 % of the administrative units for the inverse modelling and LDAS tuning methods, respectively. It is found that the LDAS tuning experiment gives more realistic values of MaxAWC and maximum Bag than the inverse modelling experiment. Using undisaggregated LAI observations leads to an underestimation of MaxAWC and maximum Bag in both experiments. Median annual maximum values of disaggregated LAI observations are found to correlate very well with MaxAWC.

scholarly journals Parameter optimisation for a better representation of drought by LSMs: inverse modelling vs. sequential data assimilation

2017 ◽  
Author(s):  
Hélène Dewaele ◽  
Simon Munier ◽  
Clément Albergel ◽  
Carole Planque ◽  
Nabil Laanaia ◽  
...  
Keyword(s):  
Data Assimilation ◽  
Land Surface ◽  
Inverse Modelling ◽  
P Value ◽  
Complex Method ◽  
Sequential Data ◽  
Area Index ◽  
Available Water Content ◽  
Cereal Grain ◽  
Land Data Assimilation System

Abstract. Soil Maximum Available Water Content (MaxAWC) is a key parameter in Land Surface Models (LSMs). However, being difficult to measure, this parameter is usually unavailable. This study assesses the feasibility of using a fifteen-year (1999–2013) time-series of satellite-derived low resolution observations of Leaf Area Index (LAI) to retrieve MaxAWC for rainfed croplands over France. LAI inter-annual variability is simulated using the CO2-responsive version of the Interactions between Soil, Biosphere and Atmosphere (ISBA) LSM for various values of MaxAWC. Optimal value is then selected by using (1) a simple inverse modelling technique, comparing simulated and observed LAI, (2) a more complex method consisting in integrating observed LAI in ISBA through a Land Data Assimilation System (LDAS) and minimizing LAI analysis increments. The evaluation of the MaxAWC retrievals from both methods is done using simulated annual maximum above-ground biomass (Bag) and straw cereal grain yield (GY) values from the Agreste French agricultural statistics portal, for 45 administrative units presenting a high proportion of straw cereals. Significant correlations (p-value 

The Role of Different Sources of Uncertainty in Data Assimilation with Coupled Land Surface - Subsurface Models

2021 ◽  
Author(s):  
Eduardo Emilio Sanchez-Leon ◽  
Natascha Brandhorst ◽  
Bastian Waldowski ◽  
Ching Pui Hung ◽  
Insa Neuweiler ◽  
...  
Keyword(s):  
Plant Physiology ◽  
Data Assimilation ◽  
Land Surface ◽  
Reference Model ◽  
Land Surface Model ◽  
Open Loop ◽  
Surface Model ◽  
Model Parameters ◽  
Area Index ◽  
Sources Of Uncertainty

&lt;p&gt;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&amp;#215;5km&lt;sup&gt;2&lt;/sup&gt;, 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.&lt;/p&gt;

scholarly journals The effect of satellite-derived surface soil moisture and leaf area index land data assimilation on streamflow simulations over France

2017 ◽  
Vol 21 (4) ◽  
pp. 2015-2033 ◽  
Author(s):  
David Fairbairn ◽  
Alina Lavinia Barbu ◽  
Adrien Napoly ◽  
Clément Albergel ◽  
Jean-François Mahfouf ◽  
...  
Keyword(s):  
Soil Moisture ◽  
Data Assimilation ◽  
Leaf Area Index ◽  
Land Surface ◽  
River Discharge ◽  
Surface Soil ◽  
Land Surface Model ◽  
Surface Model ◽  
Surface Soil Moisture ◽  
Area Index

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.

Modeling of a forested study site with the Community Land Model version 5 using climate projections for the 21st century.

2020 ◽  
Author(s):  
Lukas Strebel ◽  
Klaus Goergen ◽  
Bibi S. Naz ◽  
Heye Bogena ◽  
Harry Vereecken ◽  
...  
Keyword(s):  
Data Assimilation ◽  
Water Content ◽  
Soil Water ◽  
Soil Water Content ◽  
Land Surface ◽  
Model Parameters ◽  
Climate Projections ◽  
Modeling Framework ◽  
Area Index ◽  
Community Land Model

&lt;p&gt;Modeling forest ecosystems is important to facilitate adaptations in forest management approaches necessary to address the challenges of climate change, particularly of interest are ecohydrological states and fluxes such as soil water content, biomass, leaf area index, and evapotranspiration.&lt;/p&gt;&lt;p&gt;The community land model in its current version 5 (CLM5) simulates a broad collection of important land-surface processes; from moisture and energy partitioning, through biogeophysical processes, to surface and subsurface runoff. Additionally, CLM5 contains a biogeochemistry model (CLM5-BGC) which includes prognostic computation of vegetation states and carbon and nitrogen pools. However, CLM5 predictions are affected by uncertainty related to uncertain model forcings and parameters. Here, we use data assimilation methods to improve model performance by assimilating soil water content observations into CLM5 using the parallel data assimilation framework (PDAF).&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;The coupled modeling framework was applied to the small (38.5 ha) forested catchment W&amp;#252;stebach located in the Eifel National Park near the German-Belgian border. As part of the terrestrial environmental observatories (TERENO) network, the SoilNet sensors at the study site provide soil water content and soil temperature measurements since 2009.&lt;/p&gt;&lt;p&gt;CLM5 simulations for the period 2009-2100 were made, using local atmospheric observations for the period of 2009-2018 and an ensemble of regional climate model projections for 2019-2100. Simulations illustrate that data assimilation of soil water content improves the characterization of past model states, and that estimated model parameters and default model parameters result in different trajectories of ecohydrological states for 2019-2100. The simulations also illustrate that this site is hardly affected by increased water stress in the future.&lt;/p&gt;&lt;p&gt;The developed framework will be extended and applied for both ecosystem reanalysis as well as further simulations using climate projections across forested sites over Europe.&lt;/p&gt;

scholarly journals The Influence of Assimilating Leaf Area Index in a Land Surface Model on Global Water Fluxes and Storages

2019 ◽  
Author(s):  
Xinxuan Zhang ◽  
Viviana Maggioni ◽  
Azbina Rahman ◽  
Paul Houser ◽  
Yuan Xue ◽  
...  
Keyword(s):  
Soil Moisture ◽  
Data Assimilation ◽  
Leaf Area Index ◽  
Leaf Area ◽  
Land Surface ◽  
Water Storage ◽  
Land Surface Model ◽  
Surface Model ◽  
Area Index ◽  
Global Water

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.

scholarly journals An Ensemble Square Root Filter for the joint assimilation of surface soil moiture and leaf area index within LDAS-Monde: application over the Euro-Mediterranean region

2019 ◽  
Author(s):  
Bertrand Bonan ◽  
Clément Albergel ◽  
Yongjun Zheng ◽  
Alina Lavinia Barbu ◽  
David Fairbairn ◽  
...  
Keyword(s):  
Soil Moisture ◽  
Data Assimilation ◽  
Leaf Area Index ◽  
Land Surface ◽  
Mediterranean Region ◽  
Surface Soil ◽  
Positive Impact ◽  
Square Root ◽  
Area Index ◽  
Square Root Filter

Abstract. This paper introduces an Ensemble Square Root Filter (EnSRF), a deterministic Ensemble Kalman Filter, to the context of assimilating jointly observations of surface soil moisture (SSM) and leaf area index (LAI) in the Land Data Assimilation System LDAS-Monde. By ingesting those satellite-derived products, LDAS-Monde constrains the Interaction between Soil, Biosphere and Atmosphere (ISBA) land surface model (LSM), coupled with the CNRM (Centre National de Recherches Météorologiques) version of the Total Runoff Integrating Pathways (CTRIP), to improve the reanalysis of land surface variables (LSVs). To evaluate its ability to produce improved LSVs reanalyses, the EnSRF is compared with the Simplified Extended Kalman, which has been routinely operated in LDAS-Monde, in a real case over the well-studied Euro-Mediterranean region at a 0.25° spatial resolution between 2008 and 2017. Both data assimilation approaches provide a positive impact on SSM and LAI estimates with respect to the model alone, putting them closer to assimilated observations. SEKF and EnSRF have a similar behaviour for LAI showing performances that are influenced by the vegetation type. For SSM, EnSRF estimates tend to be closer to observations than SEKF. The impact of assimilating SSM and LAI is also assessed on unobserved soil moisture in the other layers of soil. Unobserved control variables are updated in the EnSRF through covariances and correlations sampled from the ensemble linking them to observed control variables. In our context, a strong correlation between SSM and soil moisture in deeper soil layers is exhibited, as expected, showing seasonal patterns that vary geographically. Moderate correlation and anti-correlations are also noticed between LAI and soil moisture in spring, summer and autumn, their absolute value tending to be larger for soil moisture in root-zone areas, showing that assimilating LAI can have an influence on soil moisture. Finally an independent evaluation of both assimilation approaches is conducted using satellite estimates of evapotranspiration and gross primary production (GPP) as well as measures of river discharges from gauging stations. The EnSRF shows a systematic albeit moderate improvement for evapotranspiration and GPP and a highly positive impact on river discharges, while the SEKF exhibits a more contrasting performance.

scholarly journals The Land Variational Ensemble Data Assimilation Framework: LAVENDAR v1.0.0

2020 ◽  
Vol 13 (1) ◽  
pp. 55-69 ◽  
Author(s):  
Ewan Pinnington ◽  
Tristan Quaife ◽  
Amos Lawless ◽  
Karina Williams ◽  
Tim Arkebauer ◽  
...  
Keyword(s):  
Data Assimilation ◽  
Land Surface ◽  
Time Window ◽  
Real Data ◽  
Surface Model ◽  
State Variables ◽  
Maize Crop ◽  
Area Index ◽  
Ensemble Data Assimilation ◽  
Ensemble Data

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.

Is the annual maximum leaf area index an important driver of water fluxes simulated by a land surface model in temperate forests?

2021 ◽  
pp. 1-9
Author(s):  
Audrey Maheu ◽  
Cybèle Cholet ◽  
Rebeca Cordero Montoya ◽  
Louis Duchesne
Keyword(s):  
Leaf Area Index ◽  
Leaf Area ◽  
Land Surface ◽  
Regional Scale ◽  
Surface Model ◽  
Annual Maximum ◽  
Water Fluxes ◽  
Mixed Forests ◽  
Area Index ◽  
Study Region

In land surface models, vegetation is often described using plant functional types (PFTs), a classification that aggregates plant species into a few groups based on similar characteristics. Within-PFT variability of these characteristics can introduce considerable uncertainty in the simulation of water fluxes in forests. Our objectives were to (i) compare the variability of the annual maximum leaf area index (LAImax) within and between PFTs and (ii) assess whether this variability leads to significant differences in simulated water fluxes at a regional scale. We classified our study region in southwestern Quebec (Canada) into three PFTs (evergreen needleleaf, deciduous broadleaf, and mixed forests) and characterized LAImax using remotely sensed MODIS-LAI data. We simulated water fluxes with the Canadian Land Surface Scheme (CLASS) and performed a sensitivity analysis. We found that within-PFT variability of LAImax was 1.7 times more important than variability between PFTs, with similar mean values for the two dominant PFTs, deciduous broadleaf forests (6.6 m2·m−2) and mixed forests (6.3 m2·m−2). In CLASS, varying LAImax within the observed range of values (4.0–7.5 m2·m−2) led to changes of less than 2% in mean evapotranspiration. Overall, LAImax is likely not an important driver of the spatial variability of water fluxes at the regional level.

scholarly journals An ensemble square root filter for the joint assimilation of surface soil moisture and leaf area index within the Land Data Assimilation System LDAS-Monde: application over the Euro-Mediterranean region

2020 ◽  
Vol 24 (1) ◽  
pp. 325-347 ◽  
Author(s):  
Bertrand Bonan ◽  
Clément Albergel ◽  
Yongjun Zheng ◽  
Alina Lavinia Barbu ◽  
David Fairbairn ◽  
...  
Keyword(s):  
Soil Moisture ◽  
Data Assimilation ◽  
Land Surface ◽  
Surface Soil ◽  
Positive Impact ◽  
Surface Soil Moisture ◽  
Area Index ◽  
Land Data Assimilation ◽  
Land Data Assimilation System ◽  
Data Assimilation System

Abstract. This paper introduces an ensemble square root filter (EnSRF) in the context of jointly assimilating observations of surface soil moisture (SSM) and the leaf area index (LAI) in the Land Data Assimilation System LDAS-Monde. By ingesting those satellite-derived products, LDAS-Monde constrains the Interaction between Soil, Biosphere and Atmosphere (ISBA) land surface model (LSM), coupled with the CNRM (Centre National de Recherches Météorologiques) version of the Total Runoff Integrating Pathways (CTRIP) model to improve the reanalysis of land surface variables (LSVs). To evaluate its ability to produce improved LSVs reanalyses, the EnSRF is compared with the simplified extended Kalman filter (SEKF), which has been well studied within the LDAS-Monde framework. The comparison is carried out over the Euro-Mediterranean region at a 0.25∘ spatial resolution between 2008 and 2017. Both data assimilation approaches provide a positive impact on SSM and LAI estimates with respect to the model alone, putting them closer to assimilated observations. The SEKF and the EnSRF have a similar behaviour for LAI showing performance levels that are influenced by the vegetation type. For SSM, EnSRF estimates tend to be closer to observations than SEKF values. The comparison between the two data assimilation approaches is also carried out on unobserved soil moisture in the other layers of soil. Unobserved control variables are updated in the EnSRF through covariances and correlations sampled from the ensemble linking them to observed control variables. In our context, a strong correlation between SSM and soil moisture in deeper soil layers is found, as expected, showing seasonal patterns that vary geographically. Moderate correlation and anti-correlations are also noticed between LAI and soil moisture, varying in space and time. Their absolute value, reaching their maximum in summer and their minimum in winter, tends to be larger for soil moisture in root-zone areas, showing that assimilating LAI can have an influence on soil moisture. Finally an independent evaluation of both assimilation approaches is conducted using satellite estimates of evapotranspiration (ET) and gross primary production (GPP) as well as measures of river discharges from gauging stations. The EnSRF shows a systematic albeit moderate improvement of root mean square differences (RMSDs) and correlations for ET and GPP products, but its main improvement is observed on river discharges with a high positive impact on Nash–Sutcliffe efficiency scores. Compared to the EnSRF, the SEKF displays a more contrasting performance.

Sign in / Sign up

Export Citation Format

Share Document