Estimating causal networks in biosphere–atmosphere interaction with the PCMCI approach

Local meteorological conditions and biospheric activity are tightly coupled. Understanding these links is an essential prerequisite for predicting the Earth system under climate change conditions. However, many empirical studies on the interaction between the biosphere and the atmosphere are based on correlative approaches that are not able to deduce causal paths, and only very few studies apply causal discovery methods. Here, we use a recently proposed causal graph discovery algorithm, which aims to reconstruct the causal dependency structure underlying a set of time series. We explore the potential of this method to infer temporal dependencies in biosphere-atmosphere interactions. Specifically we address the following questions: How do periodicity and heteroscedasticity influence causal detection rates, i.e. the detection of existing and non-existing links? How consistent are results for noise-contaminated data? Do results exhibit an increased information content that justifies the use of this causal-inference method? We explore the first question using artificial time series with well known dependencies that mimic real-world biosphere-atmosphere interactions. The two remaining questions are addressed jointly in two case studies utilizing observational data. Firstly, we analyse three replicated eddy covariance datasets from a Mediterranean ecosystem at half hourly time resolution allowing us to understand the impact of measurement uncertainties. Secondly, we analyse global NDVI time series (GIMMS 3g) along with gridded climate data to study large-scale climatic drivers of vegetation greenness. Overall, the results confirm the capacity of the causal discovery method to extract time-lagged linear dependencies under realistic settings. The violation of the method's assumptions increases the likelihood to detect false links. Nevertheless, we consistently identify interaction patterns in observational data. Our findings suggest that estimating a directed biosphere-atmosphere network at the ecosystem level can offer novel possibilities to unravel complex multi-directional interactions. Other than classical correlative approaches, our findings are constrained to a few meaningful set of relations which can be powerful insights for the evaluation of terrestrial ecosystem models

Seasonal adaptation of the thermal‐based two‐source energy balance model for estimating evapotranspiration in a semiarid tree‐grass ecosystem

© 2020 by the authors.The thermal-based two-source energy balance (TSEB) model has accurately simulated energy fluxes in a wide range of landscapes with both remote and proximal sensing data. However, tree-grass ecosystems (TGE) have notably complex heterogeneous vegetation mixtures and dynamic phenological characteristics presenting clear challenges to earth observation and modeling methods. Particularly, the TSEB modeling structure assumes a single vegetation source, making it difficult to represent the multiple vegetation layers present in TGEs (i.e., trees and grasses) which have different phenological and structural characteristics. This study evaluates the implementation of TSEB in a TGE located in central Spain and proposes a new strategy to consider the spatial and temporal complexities observed. This was based on sensitivity analyses (SA) conducted on both primary remote sensing inputs (local SA) and model parameters (global SA). The model was subsequently modified considering phenological dynamics in semi-arid TGEs and assuming a dominant vegetation structure and cover (i.e., either grassland or broadleaved trees) for different seasons (TSEB-2S). The adaptation was compared against the default model and evaluated against eddy covariance (EC) flux measurements and lysimeters over the experimental site. TSEB-2S vastly improved over the default TSEB performance decreasing the mean bias and root-mean-square-deviation (RMSD) of latent heat (LE) from 40 and 82 W m−2 to −4 and 59 W m−2, respectively during 2015. TSEB-2S was further validated for two other EC towers and for different years (2015, 2016 and 2017) obtaining similar error statistics with RMSD of LE ranging between 57 and 63 W m−2. The results presented here demonstrate a relatively simple strategy to improve water and energy flux monitoring over a complex and vulnerable landscape, which are often poorly represented through remote sensing models.The research received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 721995. It was also funded by Ministerio de Economía y Competitividad through FLUXPEC CGL2012-34383 and SynerTGE CGL2015-G9095-R (MINECO/FEDER, UE) projects. The research infrastructure at the measurement site in Majadas de Tiétar was partly funded through the Alexander von Humboldt Foundation, ELEMENTAL (CGL 2017-83538-C3-3-R, MINECO-FEDER) and IMAGINA (PROMETEU 2019; Generalitat Valenciana).Peer reviewe

Multiple-constraint inversion of SCOPE. Evaluating the potential of GPP and SIF for the retrieval of plant functional traits

The most recent efforts to provide remote sensing (RS) estimates of plant function rely on the combination of Radiative Transfer Models (RTM) and Soil-Vegetation-Atmosphere Transfer (SVAT) models, such as the Soil-Canopy Observation Photosynthesis and Energy fluxes (SCOPE) model. In this work we used ground spectroradiometric and chamber-based CO2 flux measurements in a nutrient manipulated Mediterranean grassland in order to: 1) develop a multiple-constraint inversion approach of SCOPE able to retrieve vegetation biochemical, structural as well as key functional traits, such as chlorophyll concentration (Cab), leaf area index (LAI), maximum carboxylation rate (Vcmax) and the Ball-Berry sensitivity parameter (m); and 2) compare the potential of the of gross primary production (GPP) and sun-induced fluorescence (SIF), together with up-welling Thermal Infrared (TIR) radiance and optical reflectance factors (RF), to estimate such parameters. The performance of the proposed inversion method as well as of the different sets of constraints was assessed with contemporary measurements of water and heat fluxes and leaf nitrogen content, using pattern-oriented model evaluation. The multiple-constraint inversion approach proposed together with the combination of optical RF and diel GPP and TIR data provided reliable estimates of parameters, and improved predicted water and heat fluxes. The addition of SIF to this scheme slightly improved the estimation of m. Parameter estimates were coherent with the variability imposed by the fertilization and the seasonality of the grassland. Results revealed that fertilization had an impact on Vcmax, while no significant differences were found for m. The combination of RF, SIF and diel TIR data weakly constrained functional traits. Approaches not including GPP failed to estimate LAI; however GPP overestimated Cab in the dry period. These problems might be related to the presence of high fractions of senescent leaves in the grassland. The proposed inversion approach together with pattern-oriented model evaluation open new perspectives for the retrieval of plant functional traits relevant for land surface models, and can be utilized at various research sites where hyperspectral remote sensing imagery and eddy covariance flux measurements are simultaneously taken

How nitrogen and phosphorus availability change water use efficiency in a Mediterranean savanna ecosystem

Nutrient availability, especially of nitrogen (N) and phosphorus (P), is of major importance for every organism and at a larger scale for ecosystem functioning and productivity. Changes in nutrient availability and potential stoichiometric imbalance due to anthropogenic nitrogen deposition might lead to nutrient deficiency or alter ecosystem functioning in various ways. In this study, we present 6 years (2014–2020) of flux-, plant-, and remote sensing data from a large-scale nutrient manipulation experiment conducted in a Mediterranean savanna-type ecosystem with an emphasis on the effects of N and P treatments on ecosystem-scale water-use efficiency (WUE) and related mechanisms. Two plots were fertilized with N (NT, 16.9 Ha) and N + P (NPT, 21.5 Ha), and a third unfertilized plot served as a control (CT). Fertilization had a strong impact on leaf nutrient stoichiometry only within the herbaceous layer with increased leaf N in both fertilized treatments and increased leaf P in NPT. Following fertilization, WUE in NT and NPT increased during the peak of growing season. While gross primary productivity similarly increased in NT and NPT, transpiration and surface conductance increased more in NT than in NPT. The results show that the NPT plot with higher nutrient availability, but more balanced N:P leaf stoichiometry had the highest WUE. On average, higher N availability resulted in a 40% increased leaf area index (LAI) in both fertilized treatments in the spring. Increased LAI reduced aerodynamic conductance and thus evaporation at both fertilized plots in the spring. Despite reduced evaporation, annual evapotranspiration increased by 10% (48.6 ± 28.3 kg H2O m−2), in the NT plot, while NPT remained similar to CT (−1%, −6.7 ± 12.2 kgH2O m−2). Potential causes for increased transpiration at NT could be increased root biomass and thus higher water uptake or rhizosphere priming to increase P-mobilization through microbes. The annual net ecosystem exchange shifted from a carbon source in CT (75.0 ± 20.6 gC m−2) to carbon-neutral in both fertilized treatments [−7.0 ± 18.5 gC m−2 (NT) 0.4 ± 22.6 gC m−2 (NPT)]. Our results show, that the N:P stoichiometric imbalance, resulting from N addition (without P), increases the WUE less than the addition of N + P, due to the strong increase in transpiration at NT, which indicates the importance of a balanced N and P content for WUE

Resolving seasonal and diel dynamics of non-rainfall water inputs in a Mediterranean ecosystem using lysimeters

The input of liquid water to terrestrial ecosystems is composed of rain and non-rainfall water (NRW). The latter comprises dew, fog, and the adsorption of atmospheric vapor on soil particle surfaces. Although NRW inputs can be relevant to support ecosystem functioning in seasonally dry ecosystems, they are understudied, being relatively small, and therefore hard to measure. In this study, we apply a partitioning routine focusing on NRW inputs over 1 year of data from large, high-precision weighing lysimeters at a semi-arid Mediterranean site. NRW inputs occur for at least 3 h on 297 d (81 % of the year), with a mean diel duration of 6 h. They reflect a pronounced seasonality as modulated by environmental conditions (i.e., temperature and net radiation). During the wet season, both dew and fog dominate NRW, while during the dry season it is mostly the soil adsorption of atmospheric water vapor. Although NRW contributes only 7.4 % to the annual water input, NRW is the only water input to the ecosystem during 15 weeks, mainly in the dry season. Benefitting from the comprehensive set of measurements at our experimental site, we show that our findings are in line with (i) independent measurements and (ii) independent model simulations forced with (near-) surface energy and moisture measurements. Furthermore, we discuss the simultaneous occurrence of soil vapor adsorption and negative eddy-covariance-derived latent heat fluxes. This study shows that NRW inputs can be reliably detected through high-resolution weighing lysimeters and a few additional measurements. Their main occurrence during nighttime underlines the necessity to consider ecosystem water fluxes at a high temporal resolution and with 24 h coverage.</p

The three major axes of terrestrial ecosystem function.

The leaf economics spectrum1,2 and the global spectrum of plant forms and functions3 revealed fundamental axes of variation in plant traits, which represent different ecological strategies that are shaped by the evolutionary development of plant species2. Ecosystem functions depend on environmental conditions and the traits of species that comprise the ecological communities4. However, the axes of variation of ecosystem functions are largely unknown, which limits our understanding of how ecosystems respond as a whole to anthropogenic drivers, climate and environmental variability4,5. Here we derive a set of ecosystem functions6 from a dataset of surface gas exchange measurements across major terrestrial biomes. We find that most of the variability within ecosystem functions (71.8%) is captured by three key axes. The first axis reflects maximum ecosystem productivity and is mostly explained by vegetation structure. The second axis reflects ecosystem water-use strategies and is jointly explained by variation in vegetation height and climate. The third axis, which represents ecosystem carbon-use efficiency, features a gradient related to aridity, and is explained primarily by variation in vegetation structure. We show that two state-of-the-art land surface models reproduce the first and most important axis of ecosystem functions. However, the models tend to simulate more strongly correlated functions than those observed, which limits their ability to accurately predict the full range of responses to environmental changes in carbon, water and energy cycling in terrestrial ecosystems7,8

Global transpiration data from sap flow measurements: The SAPFLUXNET database

Plant transpiration links physiological responses of vegetation to water supply and demand with hydrological, energy, and carbon budgets at the land-atmosphere interface. However, despite being the main land evaporative flux at the global scale, transpiration and its response to environmental drivers are currently not well constrained by observations. Here we introduce the first global compilation of whole-plant transpiration data from sap flow measurements (SAPFLUXNET, https://sapfluxnet.creaf.cat/, last access: 8 June 2021). We harmonized and quality-controlled individual datasets supplied by contributors worldwide in a semi-automatic data workflow implemented in the R programming language. Datasets include sub-daily time series of sap flow and hydrometeorological drivers for one or more growing seasons, as well as metadata on the stand characteristics, plant attributes, and technical details of the measurements. SAPFLUXNET contains 202 globally distributed datasets with sap flow time series for 2714 plants, mostly trees, of 174 species. SAPFLUXNET has a broad bioclimatic coverage, with woodland/shrubland and temperate forest biomes especially well represented (80% of the datasets). The measurements cover a wide variety of stand structural characteristics and plant sizes. The datasets encompass the period between 1995 and 2018, with 50% of the datasets being at least 3 years long. Accompanying radiation and vapour pressure deficit data are available for most of the datasets, while on-site soil water content is available for 56% of the datasets. Many datasets contain data for species that make up 90% or more of the total stand basal area, allowing the estimation of stand transpiration in diverse ecological settings. SAPFLUXNET adds to existing plant trait datasets, ecosystem flux networks, and remote sensing products to help increase our understanding of plant water use, plant responses to drought, and ecohydrological processes. SAPFLUXNET version 0.1.5 is freely available from the Zenodo repository (10.5281/zenodo.3971689; Poyatos et al., 2020a). The "sapfluxnetr"R package-designed to access, visualize, and process SAPFLUXNET data-is available from CRAN. © 2021 Rafael Poyatos et al.This research was supported by the Minis-terio de Economía y Competitividad (grant no. CGL2014-55883-JIN), the Ministerio de Ciencia e Innovación (grant no. RTI2018-095297-J-I00), the Ministerio de Ciencia e Innovación (grant no. CAS16/00207), the Agència de Gestió d’Ajuts Universitaris i de Recerca (grant no. SGR1001), the Alexander von Humboldt-Stiftung (Humboldt Research Fellowship for Experienced Researchers (RP)), and the Institució Catalana de Recerca i Estudis Avançats (Academia Award (JMV)). Víctor Flo was supported by the doctoral fellowship FPU15/03939 (MECD, Spain)