scholarly journals High-Resolution Virtual Catchment Simulations of the Subsurface-Land Surface-Atmosphere System

Author(s):  
Bernd Schalge ◽  
Jehan Rihani ◽  
Gabriele Baroni ◽  
Daniel Erdal ◽  
Gernot Geppert ◽  
...  

Abstract. Combining numerical models, which simulate water and energy fluxes in the subsurface-land surface-atmosphere system in a physically consistent way, becomes increasingly important to understand and study fluxes at compartmental boundaries and interdependencies of states across these boundaries. Complete state evolutions generated by such models, when run at highest possible resolutions while incorporating as many processes as attainable, may be regarded as a proxy of the real world – a virtual reality – which can be used to test hypotheses on functioning of the coupled terrestrial system and may serve as source for virtual measurements to develop data-assimilation methods. Such simulation systems, however, face severe problems caused by the vastly different scales of the processes acting in the compartments of the terrestrial system. The present study is motivated by the development of cross-compartmental data-assimilation methods, which face the difficulty of data scarcity in the subsurface when applied to real data. With appropriate and realistic measurement operators, the virtual reality not only allows taking virtual observations in any part of the terrestrial system at any density, thus overcoming data-scarcity problems of real-world applications, but also provides full information about true states and parameters aimed to be reconstructed from the measurements by data assimilation. In the present study, we have used the Terrestrial Systems Modeling Platform TerrSysMP, which couples the meteorological model COSMO, the land-surface model CLM, and the subsurface model ParFlow, to set up the virtual reality for a regional terrestrial system roughly oriented at the Neckar catchment in southwest Germany. We find that the virtual reality is in many aspects quite close to real observations of the catchment concerning, e.g., atmospheric boundary-layer height, precipitation, and runoff. But also discrepancies become apparent both in the ability of such models to correctly simulate some processes – which still need improvement – and the realism of the results of some observation operators like the SMOS and SMAP sensors, when faced with model states. In a succeeding step, we will use the virtual reality to generate observations in all compartments of the system for coupled data assimilation. The data assimilation will rely on a coarsened and simplified version of the model system.

2020 ◽  
Author(s):  
Bernd Schalge ◽  
Gabriele Baroni ◽  
Barbara Haese ◽  
Daniel Erdal ◽  
Gernot Geppert ◽  
...  

Abstract. Coupled numerical models, which simulate water and energy fluxes in the subsurface-land surface-atmosphere system in a physically consistent way are a prerequisite for the analysis and a better understanding of heat and matter exchange fluxes at compartmental boundaries and interdependencies of states across these boundaries. Complete state evolutions generated by such models may be regarded as a proxy of the real world, provided they are run at sufficiently high resolution and incorporate the most important processes. Such a virtual reality can be used to test hypotheses on the functioning of the coupled terrestrial system. Coupled simulation systems, however, face severe problems caused by the vastly different scales of the processes acting in and between the compartments of the terrestrial system, which also hinders comprehensive tests of their realism. We used the Terrestrial Systems Modeling Platform TerrSysMP, which couples the meteorological model COSMO, the land-surface model CLM, and the subsurface model ParFlow, to generate a virtual catchment for a regional terrestrial system mimicking the Neckar catchment in southwest Germany. Simulations for this catchment are made for the period 2007–2015, and at a spatial resolution of 400 m for the land surface and subsurface and 1.1 km for the atmosphere. Among a discussion of modelling challenges, the model performance is evaluated based on real observations covering several variables of the water cycle. We find that the simulated (virtual) catchment behaves in many aspects quite close to observations of the real Neckar catchment, e.g. concerning atmospheric boundary-layer height, precipitation, and runoff. But also discrepancies become apparent, both in the ability of the model to correctly simulate some processes which still need improvement such as overland flow, and in the realism of some observation operators like the satellite based soil moisture sensors. The whole raw dataset is available for interested users. The dataset described here is available via the CERA database (Schalge et al., 2020): https://doi.org/10.26050/WDCC/Neckar_VCS_v1.


2016 ◽  
Author(s):  
G. J. Schürmann ◽  
T. Kaminski ◽  
C. Köstler ◽  
N. Carvalhais ◽  
M. Voßbeck ◽  
...  

Abstract. We describe the Max Planck Institute Carbon Cycle Data Assimilation System (MPI-CCDAS) built around the tangent-linear version of the land surface scheme of the MPI-Earth System Model v1 (JSBACH). The simulated terrestrial biosphere processes (phenology and carbon balance) were constrained by observations of the fraction of photosynthetically active radiation (TIP-FAPAR product) and by observations of atmospheric CO2 at a global set of monitoring stations for the years 2005–2009. The system successfully, and computationally efficiently, improved average foliar area and northern extra-tropical seasonality of foliar area when constrained by TIP-FAPAR. Global net and gross carbon fluxes were improved when constrained by atmospheric CO2, although the system tended to underestimate tropical productivity. Assimilating both data streams jointly allowed the MPI-CCDAS to match both observations (TIP-FAPAR and atmospheric CO2) equally well as the single data stream assimilation cases, therefore overall increasing the appropriateness of the resultant parameter values and biosphere dynamics. Our study thus highlights the role of the TIP-FAPAR product in stabilising the underdetermined atmospheric inversion problem and demonstrates the value of multiple-data stream assimilation for the simulation of terrestrial biosphere dynamics. The constraint on regional gross and net CO2 flux patterns is limited through the parametrisation of the biosphere model. We expect improvement on that aspect through a refined initialisation strategy and inclusion of further biosphere observations as constraints.


2020 ◽  
Author(s):  
Elizabeth Cooper ◽  
Eleanor Blyth ◽  
Hollie Cooper ◽  
Rich Ellis ◽  
Ewan Pinnington ◽  
...  

Abstract. Soil moisture predictions from land surface models are important in hydrological, ecological and meteorological applications. In recent years the availability of wide-area soil-moisture measurements has increased, but few studies have combined model-based soil moisture predictions with in-situ observations beyond the point scale. Here we show that we can markedly improve soil moisture estimates from the JULES land surface model using field scale observations and data assimilation techniques. Rather than directly updating soil moisture estimates towards observed values, we optimize constants in the underlying pedotransfer functions, which relate soil texture to JULES soil physics parameters. In this way we generate a single set of newly calibrated pedotransfer functions based on observations from a number of UK sites with different soil textures. We demonstrate that calibrating a pedotransfer function in this way can improve the performance of land surface models, leading to the potential for better flood, drought and climate projections.


Author(s):  
Nemesio Rodriguez-Fernandez ◽  
Patricia de Rosnay ◽  
Clement Albergel ◽  
Philippe Richaume ◽  
Filipe Aires ◽  
...  

The assimilation of Soil Moisture and Ocean Salinity (SMOS) data into the ECMWF (European Centre for Medium Range Weather Forecasts) H-TESSEL (Hydrology revised - Tiled ECMWF Scheme for Surface Exchanges over Land) model is presented. SMOS soil moisture (SM) estimates have been produced specifically by training a neural network with SMOS brightness temperatures as input and H-TESSEL model SM simulations as reference. This can help the assimilation of SMOS information in several ways: (1) the neural network soil moisture (NNSM) data have a similar climatology to the model, (2) no global bias is present with respect to the model even if regional differences can exist. Experiments performing joint data assimilation (DA) of NNSM, 2 metre air temperature and relative humidity or NNSM-only DA are discussed. The resulting SM was evaluated against a large number of in situ measurements of SM obtaining similar results to those of the model with no assimilation, even if significant differences were found from site to site. In addition, atmospheric forecasts initialized with H-TESSEL runs (without DA) or with the analysed SM were compared to measure of the impact of the satellite information. Although, NNSM DA has an overall neutral impact in the forecast in the Tropics, a significant positive impact was found in other areas and periods, especially in regions with limited in situ information. The joint NNSM, T2m and RH2m DA improves the forecast for all the seasons in the Southern Hemisphere. The impact is mostly due to T2m and RH2m, but SMOS NN DA alone also improves the forecast in July- September. In the Northern Hemisphere, the joint NNSM, T2m and RH2m DA improves the forecast in April-September, while NNSM alone has a significant positive effect in July-September. Furthermore, forecasting skill maps show that SMOS NNSM improves the forecast in North America and in Northern Asia for up to 72 hours lead time.


2012 ◽  
Vol 9 (4) ◽  
pp. 4587-4631 ◽  
Author(s):  
W. B. Anderson ◽  
B. F. Zaitchik ◽  
C. R. Hain ◽  
M. C. Anderson ◽  
M. T. Yilmaz ◽  
...  

Abstract. Drought in East Africa is a recurring phenomenon with significant humanitarian impacts. Given the steep climatic gradients, topographic contrasts, general data scarcity, and, in places, political instability that characterize the region, there is a need for spatially distributed, remotely derived monitoring systems to inform national and international drought response. At the same time, the very diversity and data scarcity that necessitate remote monitoring also make it difficult to evaluate the reliability of these systems. Here we apply a suite of remote monitoring techniques to characterize the temporal and spatial evolution of the 2010–2011 Horn of Africa drought. Diverse satellite observations allow for evaluation of meteorological, agricultural, and hydrological aspects of drought, each of which is of interest to different stakeholders. Focusing on soil moisture, we apply triple collocation analysis (TCA) to three independent methods for estimating soil moisture anomalies to characterize relative error between products and to provide a basis for objective data merging. The three soil moisture methods evaluated include microwave remote sensing using the Advanced Microwave Scanning Radiometer – Earth Observing System (AMSR-E) sensor, thermal remote sensing using the Atmosphere-Land Exchange Inverse (ALEXI) surface energy balance algorithm, and physically-based land surface modeling using the Noah land surface model. It was found that the three soil moisture monitoring methods yield similar drought anomaly estimates in areas characterized by extremely low or by moderate vegetation cover, particularly during the below-average 2011 long rainy season. Systematic discrepancies were found, however, in regions of moderately low vegetation cover and high vegetation cover, especially during the failed 2010 short rains. The merged, TCA-weighted soil moisture composite product takes advantage of the relative strengths of each method, as judged by the consistency of anomaly estimates across independent methods. This approach holds potential as a remote soil moisture-based drought monitoring system that is robust across the diverse climatic and ecological zones of East Africa.


2021 ◽  
Author(s):  
Ann Scheliga ◽  
Manuela Girotto

<p>Sea level rise (SLR) projections rely on the accurate and precise closure of Earth’s water budget. The Gravity Recovery and Climate Experiment (GRACE) mission has provided global-coverage observations of terrestrial water storage (TWS) anomalies that improve accounting of ice and land hydrology changes and how these changes contribute to sea level rise. The contribution of land hydrology TWS changes to sea level rise is much smaller and less certain than contributions from glacial melt and thermal expansion. Although land hydrology TWS plays a smaller role, it is still important to investigate to improve the precision of the overall global water budget. This study analyzes how data assimilation techniques improve estimates of the land hydrology contribution to sea level rise. To achieve this, three global TWS datasets were analyzed: (1) GRACE TWS observations alone, (2) TWS estimates from the model-only simulation using Catchment Land Surface Model, and (3) TWS estimates from a data assimilation product of (1) and (2). We compared the data assimilation product with the GRACE observations alone and the model-only simulation to isolate the contribution to sea level rise from anthropogenic activities. We assumed a balanced water budget between land hydrology and the ocean, thus changes in global TWS are considered equal and opposite to sea level rise contribution.  Over the period of 2003-2016, we found sea level rise contributions from each dataset of +0.35 mm SLR eq/yr for GRACE, -0.34 mm SLR eq/yr for model-only, and a +0.09 mm SLR eq/yr for DA (reported as the mean linear trend). Our results indicate that the model-only simulation is not capturing important hydrologic processes. These are likely anthropogenic driven, indicating direct anthropogenic and climate-driven TWS changes play a substantial role in TWS contribution to SLR.</p>


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>


2006 ◽  
Vol 7 (3) ◽  
pp. 421-432 ◽  
Author(s):  
Wade T. Crow ◽  
Emiel Van Loon

Abstract Data assimilation approaches require some type of state forecast error covariance information in order to optimally merge model predictions with observations. The ensemble Kalman filter (EnKF) dynamically derives such information through a Monte Carlo approach and the introduction of random noise in model states, fluxes, and/or forcing data. However, in land data assimilation, relatively little guidance exists concerning strategies for selecting the appropriate magnitude and/or type of introduced model noise. In addition, little is known about the sensitivity of filter prediction accuracy to (potentially) inappropriate assumptions concerning the source and magnitude of modeling error. Using a series of synthetic identical twin experiments, this analysis explores the consequences of making incorrect assumptions concerning the source and magnitude of model error on the efficiency of assimilating surface soil moisture observations to constrain deeper root-zone soil moisture predictions made by a land surface model. Results suggest that inappropriate model error assumptions can lead to circumstances in which the assimilation of surface soil moisture observations actually degrades the performance of a land surface model (relative to open-loop assimilations that lack a data assimilation component). Prospects for diagnosing such circumstances and adaptively correcting the culpable model error assumptions using filter innovations are discussed. The dual assimilation of both runoff (from streamflow) and surface soil moisture observations appears to offer a more robust assimilation framework where incorrect model error assumptions are more readily diagnosed via filter innovations.


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.


Sign in / Sign up

Export Citation Format

Share Document