Cutting out the middleman: calibrating and validating a dynamic vegetation model (ED2-PROSPECT5) using remotely sensed surface reflectance

Ecosystem models are often calibrated and/or validated against derived remote sensing data products, such as MODIS leaf area index. However, these data products are generally based on their own models, whose assumptions may not be compatible with those of the ecosystem model in question, and whose uncertainties are usually not well quantified. Here, we develop an alternative approach whereby we modify an ecosystem model to predict full-range, high spectral resolution surface reflectance, which can then be compared directly against airborne and satellite data. Specifically, we coupled the two-stream representation of canopy radiative transfer in the Ecosystem Demography model (ED2) with a leaf radiative transfer model (PROSPECT 5) and a simple soil reflectance model. We then calibrated this model against reflectance observations from the NASA Airborne VIsible/InfraRed Imaging Spectrometer (AVIRIS) and survey data from 54 temperate forest plots in the northeastern United States. The calibration successfully constrained the posterior distributions of model parameters related to leaf biochemistry and morphology and canopy structure for five plant functional types. The calibrated model was able to accurately reproduce surface reflectance and leaf area index for sites with highly varied forest composition and structure, using a single common set of parameters across all sites. We conclude that having dynamic vegetation models directly predict surface reflectance is a promising avenue for model calibration and validation using remote sensing data.https://gmd.copernicus.org/preprints/gmd-2020-324/gmd-2020-324.pdfFirst author draf

Towards mapping biodiversity from above: Can fusing lidar and hyperspectral remote sensing predict taxonomic, functional, and phylogenetic tree diversity in temperate forests?

Aim: Rapid global change is impacting the diversity of tree species and essential ecosystem functions and services of forests. It is therefore critical to understand and predict how the diversity of tree species is spatially distributed within and among forest biomes. Satellite remote sensing platforms have been used for decades to map forest structure and function but are limited in their capacity to monitor change by their relatively coarse spatial resolution and the complexity of scales at which different dimensions of biodiversity are observed in the field. Recently, airborne remote sensing platforms making use of passive high spectral resolution (i.e., hyperspectral) and active lidar data have been operationalized, providing an opportunity to disentangle how biodiversity patterns vary across space and time from field observations to larger scales. Most studies to date have focused on single sites and/or one sensor type; here we ask how multiple sensor types from the National Ecological Observatory Network’s Airborne Observation Platform (NEON AOP) perform across multiple sites in a single biome at the NEON field plot scale (i.e., 40 m × 40 m).Location: Eastern USA.Time period: 2017– 2018.Taxa studied: Trees.Methods: With a fusion of hyperspectral and lidar data from the NEON AOP, we as-sess the ability of high resolution remotely sensed metrics to measure biodiversity variation across eastern US temperate forests. We examine how taxonomic, functional, and phylogenetic measures of alpha diversity vary spatially and assess to what degree remotely sensed metrics correlate with in situ biodiversity metrics.Results: Models using estimates of forest function, canopy structure, and topographic diversity performed better than models containing each category alone. Our results show that canopy structural diversity, and not just spectral reflectance, is critical to predicting biodiversity.Main conclusions: We found that an approach that jointly leverages spectral properties related to leaf and canopy functional traits and forest health, lidar derived estimates of forest structure, fine-resolution topographic diversity, and careful consideration of biogeographical differences within and among biomes is needed to accurately map biodiversity variation from above

Hyperspectral reflectance as a tool to measure biochemical and physiological traits in wheat

Simultaneous indirect assessment of multiple and diverse plant parameters in an exact and expeditious manner is becoming imperative in irrigated arid regions, with a view toward creating drought-tolerant genotypes or for the management of precision irrigation. This study aimed to evaluate whether spectral reflectance indices (SRIs) in three parts of the electromagnetic spectrum ((visible-infrared (VIS), near-infrared (NIR)), and shortwave-infrared (SWIR)) could be used to track changes in morphophysiological parameters of wheat cultivars exposed to 1.00, 0.75, and 0.50 of the estimated evapotranspiration (ETc). Significant differences were found in the parameters of growth and photosynthetic efficiency, and canopy spectral reflectance among the three cultivars subjected to different irrigation rates. All parameters were highly and significantly correlated with each other particularly under the 0.50 ETc treatment. The VIS/VIS- and NIR/VIS-based indices were sufficient and suitable for assessing the growth and photosynthetic properties of wheat cultivars similar to those indices based on NIR/NIR, SWIR/NIR, or SWIR/SWIR. Almost all tested SRIs proved to assess growth and photosynthetic parameters, including transpiration rate, more efficiently when regressions were analyzed for each water irrigation rate individually. This study, the type of which has rarely been conducted in irrigated arid regions, indicates that spectral reflectance data can be used as a rapid and non-destructive alternative method for assessment of the growth and photosynthetic efficiency of wheat under a range of water irrigation rates.This work was supported by the Sustainable Modernization of Traditional Agriculture (MasAgro) initiative from the Secretariat of Agriculture, Livestock, Rural Development, Fisheries and Food (SAGARPA) and a collaboration between CIMMYT, CSIRO and The Australian National University. We acknowledge financial support for the PhD scholarship of VSP from CONACYT, Mexico number 207607, the financial support of the Australian Research Council Centre of Excellence for Translational Photosynthesis (CE140100015) and access to glasshouse and field experiments supported by Grains Research & Development Corporation funding (CSP00168). SPS was supported by the United States Department of Energy contract No. DE-SC0012704 to Brookhaven National Laborator

Identification of key parameters controlling demographically structured vegetation dynamics in a land surface model: CLM4.5(FATES)

Vegetation plays an important role in regulating global carbon cycles and is a key component of the Earth system models (ESMs) that aim to project Earth\u27s future climate. In the last decade, the vegetation component within ESMs has witnessed great progress from simple “big-leaf” approaches to demographically structured approaches, which have a better representation of plant size, canopy structure, and disturbances. These demographically structured vegetation models typically have a large number of input parameters, and sensitivity analysis is needed to quantify the impact of each parameter on the model outputs for a better understanding of model behavior. In this study, we conducted a comprehensive sensitivity analysis to diagnose the Community Land Model coupled to the Functionally Assembled Terrestrial Simulator, or CLM4.5(FATES). Specifically, we quantified the first- and second-order sensitivities of the model parameters to outputs that represent simulated growth and mortality as well as carbon fluxes and stocks for a tropical site with an extent of 1×1∘. While the photosynthetic capacity parameter (Vc,max25) is found to be important for simulated carbon stocks and fluxes, we also show the importance of carbon storage and allometry parameters, which determine survival and growth strategies within the model. The parameter sensitivity changes with different sizes of trees and climate conditions. The results of this study highlight the importance of understanding the dynamics of the next generation of demographically enabled vegetation models within ESMs to improve model parameterization and structure for better model fidelity

Hydraulic architecture explains species moisture dependency but not mortality rates across a tropical rainfall gradient

Intensified droughts are affecting tropical forests across the globe. However, the underlying mechanisms of tree drought response and mortality are poorly understood. Hydraulic traits and especially hydraulic safety margins (HSMs), that is, the extent to which plants buffer themselves from thresholds of water stress, provide insights into species-specific drought vulnerability. We investigated hydraulic traits during an intense drought triggered by the 2015–2016 El Niño on 27 canopy tree species across three tropical forest sites with differing precipitation. We capitalized on the drought event as a time when plant water status might approach or exceed thresholds of water stress. We investigated the degree to which these traits varied across the rainfall gradient, as well as relationships among hydraulic traits and species-specific optimal moisture and mortality rates. There were no differences among sites for any measured trait. There was strong coordination among traits, with a network analysis revealing two major groups of coordinated traits. In one group, there were water potentials, turgor loss point, sapwood capacitance and density, HSMs, and mortality rate. In the second group, there was leaf mass per area, leaf dry matter content, hydraulic architecture (leaf area to sapwood area ratio), and species-specific optimal moisture. These results demonstrated that while species with greater safety from turgor loss had lower mortality rates, hydraulic architecture was the only trait that explained species’ moisture dependency. Species with a greater leaf area to sapwood area ratio were associated with drier sites and reduced their transpirational demand during the dry season via deciduousness

Convergence in relationships between leaf traits, spectra and age across diverse canopy environments and two contrasting tropical forests

• Leaf age structures the phenology and development of plants, as well as the evolution of leaf traits over life histories. However, a general method for efficiently estimating leaf age across forests and canopy environments is lacking. • Here, we explored the potential for a statistical model, previously developed for Peruvian sunlit leaves, to consistently predict leaf ages from leaf reflectance spectra across two contrasting forests in Peru and Brazil and across diverse canopy environments. • The model performed well for independent Brazilian sunlit and shade canopy leaves (R2 = 0.75–0.78), suggesting that canopy leaves (and their associated spectra) follow constrained developmental trajectories even in contrasting forests. The model did not perform as well for mid-canopy and understory leaves (R2 = 0.27–0.29), because leaves in different environments have distinct traits and trait developmental trajectories. When we accounted for distinct environment–trait linkages – either by explicitly including traits and environments in the model, or, even better, by re-parameterizing the spectra-only model to implicitly capture distinct trait-trajectories in different environments – we achieved a more general model that well-predicted leaf age across forests and environments (R2 = 0.79). • Fundamental rules, linked to leaf environments, constrain the development of leaf traits and allow for general prediction of leaf age from spectra across species, sites and canopy environments

High-throughput characterization, correlation, and mapping of leaf photosynthetic and functional traits in the soybean (Glycine max) nested association mapping population

Photosynthesis is a key target to improve crop production in many species including soybean [Glycine max (L.) Merr.]. A challenge is that phenotyping photosynthetic traits by traditional approaches is slow and destructive. There is proof-of-concept for leaf hyperspectral reflectance as a rapid method to model photosynthetic traits. However, the crucial step of demonstrating that hyperspectral approaches can be used to advance understanding of the genetic architecture of photosynthetic traits is untested. To address this challenge, we used full-range (500-2,400 nm) leaf reflectance spectroscopy to build partial least squares regression models to estimate leaf traits, including the rate-limiting processes of photosynthesis, maximum Rubisco carboxylation rate, and maximum electron transport. In total, 11 models were produced from a diverse population of soybean sampled over multiple field seasons to estimate photosynthetic parameters, chlorophyll content, leaf carbon and leaf nitrogen percentage, and specific leaf area (with R2 from 0.56 to 0.96 and root mean square error approximately \u3c10% of the range of calibration data). We explore the utility of these models by applying them to the soybean nested association mapping population, which showed variability in photosynthetic and leaf traits. Genetic mapping provided insights into the underlying genetic architecture of photosynthetic traits and potential improvement in soybean. Notably, the maximum Rubisco carboxylation rate mapped to a region of chromosome 19 containing genes encoding multiple small subunits of Rubisco. We also mapped the maximum electron transport rate to a region of chromosome 10 containing a fructose 1,6-bisphosphatase gene, encoding an important enzyme in the regeneration of ribulose 1,5-bisphosphate and the sucrose biosynthetic pathway. The estimated rate-limiting steps of photosynthesis were low or negatively correlated with yield suggesting that these traits are not influenced by the same genetic mechanisms and are not limiting yield in the soybean NAM population. Leaf carbon percentage, leaf nitrogen percentage, and specific leaf area showed strong correlations with yield and may be of interest in breeding programs as a proxy for yield. This work is among the first to use hyperspectral reflectance to model and map the genetic architecture of the rate-limiting steps of photosynthesis

One Stomatal Model to Rule Them All?:Toward Improved Representation of Carbon and Water Exchange in Global Models

Stomatal conductance schemes that optimize with respect to photosynthetic and hydraulic functions have been proposed to address biases in land-surface model (LSM) simulations during drought. However, systematic evaluations of both optimality-based and alternative empirical formulations for coupling carbon and water fluxes are lacking. Here, we embed 12 empirical and optimization approaches within a LSM framework. We use theoretical model experiments to explore parameter identifiability and understand how model behaviors differ in response to abiotic changes. We also evaluate the models against leaf-level observations of gas-exchange and hydraulic variables, from xeric to wet forest/woody species spanning a mean annual precipitation range of 361–3,286 mm yr−1. We find that models differ in how easily parameterized they are, due to: (a) poorly constrained optimality criteria (i.e., resulting in multiple solutions), (b) low influence parameters, (c) sensitivities to environmental drivers. In both the idealized experiments and compared to observations, sensitivities to variability in environmental drivers do not agree among models. Marked differences arise in sensitivities to soil moisture (soil water potential) and vapor pressure deficit. For example, stomatal closure rates at high vapor pressure deficit range between −45% and +70% of those observed. Although over half the new generation of stomatal schemes perform to a similar standard compared to observations of leaf-gas exchange, two models do so through large biases in simulated leaf water potential (up to 11 MPa). Our results provide guidance for LSM development, by highlighting key areas in need for additional experimentation and theory, and by constraining currently viable stomatal hypotheses

A test of the 'one-point method' for estimating maximum carboxylation capacity from field-measured, light-saturated photosynthesis

Simulations of photosynthesis by terrestrial biosphere models typically need a specification of the maximum carboxylation rate (Vcmax). Estimating this parameter using A–Ci curves (net photosynthesis, A, vs intercellular CO2 concentration, Ci) is laborious, which limits availability of Vcmax data. However, many multispecies field datasets include net photosynthetic rate at saturating irradiance and at ambient atmospheric CO2 concentration (Asat) measurements, from which Vcmax can be extracted using a 'one-point method'.\ud \ud We used a global dataset of A–Ci curves (564 species from 46 field sites, covering a range of plant functional types) to test the validity of an alternative approach to estimate Vcmax from Asat via this 'one-point method'.\ud \ud If leaf respiration during the day (Rday) is known exactly, Vcmax can be estimated with an r2 value of 0.98 and a root-mean-squared error (RMSE) of 8.19 μmol m−2 s−1. However, Rday typically must be estimated. Estimating Rday as 1.5% of Vcmax, we found that Vcmax could be estimated with an r2 of 0.95 and an RMSE of 17.1 μmol m−2 s−1.\ud \ud The one-point method provides a robust means to expand current databases of field-measured Vcmax, giving new potential to improve vegetation models and quantify the environmental drivers of Vcmax variation