The remote monitoring of plant canopies is critically needed for understanding of terrestrial ecosystem mechanics and biodiversity as well as capturing the short- to long-term responses of vegetation ...to disturbance and climate change. A variety of orbital, sub-orbital, and field instruments have been used to retrieve optical spectral signals and to study different vegetation properties such as plant biochemistry, nutrient cycling, physiology, water status, and stress. Radiative transfer models (RTMs) provide a mechanistic link between vegetation properties and observed spectral features, and RTM spectral inversion is a useful framework for estimating these properties from spectral data. However, existing approaches to RTM spectral inversion are typically limited by the inability to characterize uncertainty in parameter estimates. Here, we introduce a Bayesian algorithm for the spectral inversion of the PROSPECT 5 leaf RTM that is distinct from past approaches in two important ways: First, the algorithm only uses reflectance and does not require transmittance observations, which have been plagued by a variety of measurement and equipment challenges. Second, the output is not a point estimate for each parameter but rather the joint probability distribution that includes estimates of parameter uncertainties and covariance structure. We validated our inversion approach using a database of leaf spectra together with measurements of equivalent water thickness (EWT) and leaf dry mass per unit area (LMA). The parameters estimated by our inversion were able to accurately reproduce the observed reflectance (RMSEVIS=0.0063, RMSENIR-SWIR=0.0098) and transmittance (RMSEVIS=0.0404, RMSENIR-SWIR=0.0551) for both broadleaved and conifer species. Inversion estimates of EWT and LMA for broadleaved species agreed well with direct measurements (CVEWT=18.8%, CVLMA=24.5%), while estimates for conifer species were less accurate (CVEWT=53.2%, CVLMA=63.3%). To examine the influence of spectral resolution on parameter uncertainty, we simulated leaf reflectance as observed by ten common remote sensing platforms with varying spectral configurations and performed a Bayesian inversion on the resulting spectra. We found that full-range hyperspectral platforms were able to retrieve all parameters accurately and precisely, while the parameter estimates of multispectral platforms were much less precise and prone to bias at high and low values. We also observed that variations in the width and location of spectral bands influenced the shape of the covariance structure of parameter estimates. Our Bayesian spectral inversion provides a powerful and versatile framework for future RTM development and single- and multi-instrumental remote sensing of vegetation.
•Novel Bayesian algorithm for RTM inversion is developed.•Accuracy comparable to past studies despite reflectance as sole input•PROSPECT can be inverted on idealized data from all sensors.•Inversion precision increases with spectral resolution.
Ecologists increasingly rely on complex computer simulations to forecast ecological systems. To make such forecasts precise, uncertainties in model parameters and structure must be reduced and ...correctly propagated to model outputs. Naively using standard statistical techniques for this task, however, can lead to bias and underestimation of uncertainties in parameters and predictions. Here, we explain why these problems occur and propose a framework for robust inference with complex computer simulations. After having identified that model error is more consequential in complex computer simulations, due to their more pronounced nonlinearity and interconnectedness, we discuss as possible solutions data rebalancing and adding bias corrections on model outputs or processes during or after the calibration procedure. We illustrate the methods in a case study, using a dynamic vegetation model. We conclude that developing better methods for robust inference of complex computer simulations is vital for generating reliable predictions of ecosystem responses.
Model error is a major problem for statistical inference with complex computer simulations due to their more pronounced nonlinearity and interconnectedness. Here, we propose a framework for robust inference including rebalancing of data and adding bias corrections on model outputs or processes during or after calibration. We conclude that methods for robust inference of complex computer simulations are vital for generating useful predictions of ecosystems responses.
We analysed the responses of 11 ecosystem models to elevated atmospheric CO₂ (eCO₂) at two temperate forest ecosystems (Duke and Oak Ridge National Laboratory (ORNL) Free‐Air CO₂ Enrichment (FACE) ...experiments) to test alternative representations of carbon (C)–nitrogen (N) cycle processes. We decomposed the model responses into component processes affecting the response to eCO₂ and confronted these with observations from the FACE experiments. Most of the models reproduced the observed initial enhancement of net primary production (NPP) at both sites, but none was able to simulate both the sustained 10‐yr enhancement at Duke and the declining response at ORNL: models generally showed signs of progressive N limitation as a result of lower than observed plant N uptake. Nonetheless, many models showed qualitative agreement with observed component processes. The results suggest that improved representation of above‐ground–below‐ground interactions and better constraints on plant stoichiometry are important for a predictive understanding of eCO₂ effects. Improved accuracy of soil organic matter inventories is pivotal to reduce uncertainty in the observed C–N budgets. The two FACE experiments are insufficient to fully constrain terrestrial responses to eCO₂, given the complexity of factors leading to the observed diverging trends, and the consequential inability of the models to explain these trends. Nevertheless, the ecosystem models were able to capture important features of the experiments, lending some support to their projections.
Earth system models (ESMs) have been developed to represent the role of terrestrial ecosystems on the energy, water, and carbon cycles. However, many ESMs still lack representation of ...within-ecosystem heterogeneity and diversity. In this paper, we present the Ecosystem Demography model version 2.2 (ED-2.2). In ED-2.2, the biophysical and physiological processes account for the horizontal and vertical heterogeneity of the ecosystem: the energy, water, and carbon cycles are solved separately for a series of vegetation cohorts (groups of individual plants of similar size and plant functional type) distributed across a series of spatially implicit patches (representing collections of micro-environments that have a similar disturbance history). We define the equations that describe the energy, water, and carbon cycles in terms of total energy, water, and carbon, which simplifies the differential equations and guarantees excellent conservation of these quantities in long-term simulation (< 0.1 % error over 50 years). We also show examples of ED-2.2 simulation results at single sites and across tropical South America. These results demonstrate the model's ability to characterize the variability of ecosystem structure, composition, and functioning both at stand and continental scales. A detailed model evaluation was conducted and is presented in a companion paper (Longo et al., 2019a). Finally, we highlight some of the ongoing model developments designed to improve the model's accuracy and performance and to include processes hitherto not represented in the model.
Understanding the manner in which changes in disturbance regimes will affect forest biodiversity is an important goal of global change research. Prevailing theories of recruitment after disturbance ...center on the role of pioneer species; predictions of forest biodiversity focus almost exclusively on dispersal and shade tolerance while vegetative reproduction is virtually omitted from models and serious discussions of the topic. However, the persistence of live damaged trees increases understory shade, generates fine-scale environmental heterogeneity, and moderates ecosystem responses to damage, while the sprouting of damaged trees offers a shortcut to reestablishment of the canopy. While a number of studies document snapshots of post-disturbance vegetative reproduction, we lack an understanding of the underlying demographic processes needed in order to both comprehend and predict observed patterns. In this study we quantify the abundance, competitive ability, and interspecific variability of vegetative reproduction in 18 replicated experimental gaps in the southern Appalachians and Carolina Piedmont, USA, in order to assess the potential role of sprouting in driving gap dynamics. Annual rates of damaged adult survival, sprout initiation, growth, and mortality were monitored over four years and compared to the performance of gap-regenerating saplings. Recruitment from sprouts was found to constitute 26—87% of early gap regeneration and forms the dominant pathway of regeneration for some species. Sprouts from recently damaged trees also grow significantly faster than the saplings with which they compete. For all measured demographic rates (damaged tree survival, sprout initiation, number, growth, and survival) differences among species are large and consistent across sites, suggesting that vegetative reproduction is an important and non-neutral process in shaping community composition. Sprouting ability does not correlate strongly with other life-history trade-offs, thus sprouting potentially provides an alternate trait axis in promoting diversity.
Phenology, by controlling the seasonal activity of vegetation on the land surface, plays a fundamental role in regulating photosynthesis and other ecosystem processes, as well as competitive ...interactions and feedbacks to the climate system. We conducted an analysis to evaluate the representation of phenology, and the associated seasonality of ecosystem‐scale CO2 exchange, in 14 models participating in the North American Carbon Program Site Synthesis. Model predictions were evaluated using long‐term measurements (emphasizing the period 2000–2006) from 10 forested sites within the AmeriFlux and Fluxnet‐Canada networks. In deciduous forests, almost all models consistently predicted that the growing season started earlier, and ended later, than was actually observed; biases of 2 weeks or more were typical. For these sites, most models were also unable to explain more than a small fraction of the observed interannual variability in phenological transition dates. Finally, for deciduous forests, misrepresentation of the seasonal cycle resulted in over‐prediction of gross ecosystem photosynthesis by +160 ± 145 g C m−2 yr−1 during the spring transition period and +75 ± 130 g C m−2 yr−1 during the autumn transition period (13% and 8% annual productivity, respectively) compensating for the tendency of most models to under‐predict the magnitude of peak summertime photosynthetic rates. Models did a better job of predicting the seasonality of CO2 exchange for evergreen forests. These results highlight the need for improved understanding of the environmental controls on vegetation phenology and incorporation of this knowledge into better phenological models. Existing models are unlikely to predict future responses of phenology to climate change accurately and therefore will misrepresent the seasonality and interannual variability of key biosphere–atmosphere feedbacks and interactions in coupled global climate models.
EuroAmerican land-use and its legacies have transformed forest structure and composition across the United States (US). More accurate reconstructions of historical states are critical to ...understanding the processes governing past, current, and future forest dynamics. Here we present new gridded (8x8km) reconstructions of pre-settlement (1800s) forest composition and structure from the upper Midwestern US (Minnesota, Wisconsin, and most of Michigan), using 19th Century Public Land Survey System (PLSS), with estimates of relative composition, above-ground biomass, stem density, and basal area for 28 tree types. This mapping is more robust than past efforts, using spatially varying correction factors to accommodate sampling design, azimuthal censoring, and biases in tree selection.
We compare pre-settlement to modern forests using US Forest Service Forest Inventory and Analysis (FIA) data to show the prevalence of lost forests (pre-settlement forests with no current analog), and novel forests (modern forests with no past analogs). Differences between pre-settlement and modern forests are spatially structured owing to differences in land-use impacts and accompanying ecological responses. Modern forests are more homogeneous, and ecotonal gradients are more diffuse today than in the past. Novel forest assemblages represent 28% of all FIA cells, and 28% of pre-settlement forests no longer exist in a modern context. Lost forests include tamarack forests in northeastern Minnesota, hemlock and cedar dominated forests in north-central Wisconsin and along the Upper Peninsula of Michigan, and elm, oak, basswood and ironwood forests along the forest-prairie boundary in south central Minnesota and eastern Wisconsin. Novel FIA forest assemblages are distributed evenly across the region, but novelty shows a strong relationship to spatial distance from remnant forests in the upper Midwest, with novelty predicted at between 20 to 60km from remnants, depending on historical forest type. The spatial relationships between remnant and novel forests, shifts in ecotone structure and the loss of historic forest types point to significant challenges for land managers if landscape restoration is a priority. The spatial signals of novelty and ecological change also point to potential challenges in using modern spatial distributions of species and communities and their relationship to underlying geophysical and climatic attributes in understanding potential responses to changing climate. The signal of human settlement on modern forests is broad, spatially varying and acts to homogenize modern forests relative to their historic counterparts, with significant implications for future management.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Plant functional traits determine vegetation responses to environmental variation, but variation in trait values is large, even within a single site. Likewise, uncertainty in how these traits map to ...Earth system feedbacks is large. We use a vegetation demographic model (VDM), the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), to explore parameter sensitivity of model predictions, and comparison to observations, at a tropical forest site: Barro Colorado Island in Panama. We define a single 12-dimensional distribution of plant trait variation, derived primarily from observations in Panama, and define plant functional types (PFTs) as random draws from this distribution. We compare several model ensembles, where individual ensemble members vary only in the plant traits that define PFTs, and separate ensembles differ from each other based on either model structural assumptions or non-trait, ecosystem-level parameters, which include (a) the number of competing PFTs present in any simulation and (b) parameters that govern disturbance and height-based light competition. While single-PFT simulations are roughly consistent with observations of productivity at Barro Colorado Island, increasing the number of competing PFTs strongly shifts model predictions towards higher productivity and biomass forests. Different ecosystem variables show greater sensitivity than others to the number of competing PFTs, with the predictions that are most dominated by large trees, such as biomass, being the most sensitive. Changing disturbance and height-sorting parameters, i.e., the rules of competitive trait filtering, shifts regimes of dominance or coexistence between early- and late-successional PFTs in the model. Increases to the extent or severity of disturbance, or to the degree of determinism in height-based light competition, all act to shift the community towards early-successional PFTs. In turn, these shifts in competitive outcomes alter predictions of ecosystem states and fluxes, with more early-successional-dominated forests having lower biomass. It is thus crucial to differentiate between plant traits, which are under competitive pressure in VDMs, from those model parameters that are not and to better understand the relationships between these two types of model parameters to quantify sources of uncertainty in VDMs.
The newest version of the Geostationary Operational Environmental Satellite series (GOES-16 and GOES-17) includes a near infrared band that allows for the calculation of normalized difference ...vegetation index (NDVI) at a 1 km at nadir spatial resolution every five minutes throughout the continental United States and every ten minutes for much of the western hemisphere. The usefulness of individual NDVI observations is limited due to the noise that remains even after cloud masks and data quality flags are applied, as much of this noise is negatively biased due to scattering within the atmosphere. Fortunately, high temporal resolution NDVI allows for the identification of consistent diurnal patterns. Here, we present a novel statistical model that utilizes this pattern, by fitting double exponential curves to the diurnal NDVI data, to provide a daily estimate of NDVI over forests that is less sensitive to noise by accounting for both random observation errors and atmospheric scattering biases. We fit this statistical model to 350 days of observations for fifteen deciduous broadleaf sites in the United States and compared the method to several simpler potential methods. Of the days 60% had more than ten observations and were able to be modeled via our methodology. Of the modeled days 72% produced daily NDVI estimates with <0.1 wide 95% confidence intervals. Of the modeled days 13% were able to provide a confident NDVI value even if there were less than five observations between 10:00–14:00. This methodology provides estimates for daily midday NDVI values with robust uncertainty estimates, even in the face of biased errors and missing midday observations.