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.
Large-scale environmental sequencing efforts have transformed our understanding of the spatial controls over soil microbial community composition and turnover. Yet, our knowledge of temporal controls ...is comparatively limited. This is a major uncertainty in microbial ecology, as there is increasing evidence that microbial community composition is important for predicting microbial community function in the future. Here, we use continental- and global-scale soil fungal community surveys, focused within northern temperate latitudes, to estimate the relative contribution of time and space to soil fungal community turnover. We detected large intra-annual temporal differences in soil fungal community similarity, where fungal communities differed most among seasons, equivalent to the community turnover observed over thousands of kilometers in space. inter-annual community turnover was comparatively smaller than intra-annual turnover. Certain environmental covariates, particularly climate covariates, explained some spatial-temporal effects, though it is unlikely the same mechanisms drive spatial vs. temporal turnover. However, these commonly measured environmental covariates could not fully explain relationships between space, time and community composition. These baseline estimates of fungal community turnover in time provide a starting point to estimate the potential duration of legacies in microbial community composition and function.
Predicted responses of transpiration to elevated atmospheric CO2 concentration (eCO2) are highly variable amongst process‐based models. To better understand and constrain this variability amongst ...models, we conducted an intercomparison of 11 ecosystem models applied to data from two forest free‐air CO2 enrichment (FACE) experiments at Duke University and Oak Ridge National Laboratory. We analysed model structures to identify the key underlying assumptions causing differences in model predictions of transpiration and canopy water use efficiency. We then compared the models against data to identify model assumptions that are incorrect or are large sources of uncertainty. We found that model‐to‐model and model‐to‐observations differences resulted from four key sets of assumptions, namely (i) the nature of the stomatal response to elevated CO2 (coupling between photosynthesis and stomata was supported by the data); (ii) the roles of the leaf and atmospheric boundary layer (models which assumed multiple conductance terms in series predicted more decoupled fluxes than observed at the broadleaf site); (iii) the treatment of canopy interception (large intermodel variability, 2–15%); and (iv) the impact of soil moisture stress (process uncertainty in how models limit carbon and water fluxes during moisture stress). Overall, model predictions of the CO2 effect on WUE were reasonable (intermodel μ = approximately 28% ± 10%) compared to the observations (μ = approximately 30% ± 13%) at the well‐coupled coniferous site (Duke), but poor (intermodel μ = approximately 24% ± 6%; observations μ = approximately 38% ± 7%) at the broadleaf site (Oak Ridge). The study yields a framework for analysing and interpreting model predictions of transpiration responses to eCO2, and highlights key improvements to these types of models.
The daunting complexity of ecosystems has led ecologists to use mathematical modelling to gain understanding of ecological relationships, processes and dynamics. In pursuit of mathematical ...tractability, these models use simplified descriptions of key patterns, processes and relationships observed in nature. In contrast, ecological data are often complex, scale‐dependent, space‐time correlated, and governed by nonlinear relations between organisms and their environment. This disparity in complexity between ecosystem models and data has created a large gap in ecology between model and data‐driven approaches. Here, we explore data assimilation (DA) with the Ensemble Kalman filter to fuse a two‐predator‐two‐prey model with abundance data from a 2600+ day experiment of a plankton community. We analyse how frequently we must assimilate measured abundances to predict accurately population dynamics, and benchmark our population model's forecast horizon against a simple null model. Results demonstrate that DA enhances the predictability and forecast horizon of complex community dynamics.
The power of forecasts to advance ecological theory Lewis, Abigail S. L.; Rollinson, Christine R.; Allyn, Andrew J. ...
Methods in ecology and evolution,
March 2023, 2023-03-00, 20230301, 2023-03-01, Letnik:
14, Številka:
3
Journal Article
Recenzirano
Odprti dostop
Ecological forecasting provides a powerful set of methods for predicting short‐ and long‐term change in living systems. Forecasts are now widely produced, enabling proactive management for many ...applied ecological problems. However, despite numerous calls for an increased emphasis on prediction in ecology, the potential for forecasting to accelerate ecological theory development remains underrealized.
Here, we provide a conceptual framework describing how ecological forecasts can energize and advance ecological theory. We emphasize the many opportunities for future progress in this area through increased forecast development, comparison and synthesis.
Our framework describes how a forecasting approach can shed new light on existing ecological theories while also allowing researchers to address novel questions. Through rigorous and repeated testing of hypotheses, forecasting can help to refine theories and understand their generality across systems. Meanwhile, synthesizing across forecasts allows for the development of novel theory about the relative predictability of ecological variables across forecast horizons and scales.
We envision a future where forecasting is integrated as part of the toolset used in fundamental ecology. By outlining the relevance of forecasting methods to ecological theory, we aim to decrease barriers to entry and broaden the community of researchers using forecasting for fundamental ecological insight.
Robust ecological forecasting of tree growth under future climate conditions is critical to anticipate future forest carbon storage and flux. Here, we apply three ingredients of ecological ...forecasting that are key to improving forecast skill: data fusion, confronting model predictions with new data, and partitioning forecast uncertainty. Specifically, we present the first fusion of tree‐ring and forest inventory data within a Bayesian state‐space model at a multi‐site, regional scale, focusing on Pinus ponderosa var. brachyptera in the southwestern US. Leveraging the complementarity of these two data sources, we parsed the ecological complexity of tree growth into the effects of climate, tree size, stand density, site quality, and their interactions, and quantified uncertainties associated with these effects. New measurements of trees, an ongoing process in forest inventories, were used to confront forecasts of tree diameter with observations, and evaluate alternative tree growth models. We forecasted tree diameter and increment in response to an ensemble of climate change projections, and separated forecast uncertainty into four different causes: initial conditions, parameters, climate drivers, and process error. We found a strong negative effect of fall–spring maximum temperature, and a positive effect of water‐year precipitation on tree growth. Furthermore, tree vulnerability to climate stress increases with greater competition, with tree size, and at poor sites. Under future climate scenarios, we forecast increment declines of 22%–117%, while the combined effect of climate and size‐related trends results in a 56%–91% decline. Partitioning of forecast uncertainty showed that diameter forecast uncertainty is primarily caused by parameter and initial conditions uncertainty, but increment forecast uncertainty is mostly caused by process error and climate driver uncertainty. This fusion of tree‐ring and forest inventory data lays the foundation for robust ecological forecasting of aboveground biomass and carbon accounting at tree, plot, and regional scales, including iterative improvement of model skill.
Great enthusiasm surrounds forest carbon sequestration as a way to offset carbon emissions and address the climate crisis, but ecological complexity contributes to scientific uncertainty about forest responses to climate change. Here, we fused together two important sources of information on tree growth – forest inventory and tree‐ring data – to quantify multiple drivers of Pinus ponderosa tree growth (climate sensitivity, competition, tree size, and site quality), and their interactions. Ecological forecasting approaches indicate that growth will decline by 56%–91% under future climate conditions and allow us to identify the causes of uncertainty surrounding future tree growth.
How, where, and why carbon (C) moves into and out of an ecosystem through time are long-standing questions in biogeochemistry. Here, we bring together hundreds of thousands of C-cycle observations at ...the Harvard Forest in central Massachusetts, USA, a mid-latitude landscape dominated by 80–120-yr-old closed-canopy forests. These data answered four questions: (1) where and how much C is presently stored in dominant forest types; (2) what are current rates of C accrual and loss; (3) what biotic and abiotic factors contribute to variability in these rates; and (4) how has climate change affected the forest's C cycle? Harvard Forest is an active C sink resulting from forest regrowth following land abandonment. Soil and tree biomass comprise nearly equal portions of existing C stocks. Net primary production (NPP) averaged 680–750 g C.m⁻².yr⁻¹; belowground NPP contributed 38–47% of the total, but with large uncertainty. Mineral soil C measured in the same inventory plots in 1992 and 2013 was too heterogeneous to detect change in soil-C pools; however, radiocarbon data suggest a small but persistent sink of 10–30 g C.m⁻².yr⁻¹. Net ecosystem production (NEP) in hardwood stands averaged ~300 g C.m⁻².yr⁻¹. NEP in hemlock-dominated forests averaged ~450 g C.m⁻².yr⁻¹ until infestation by the hemlock woolly adelgid turned these stands into a net C source. Since 2000, NPP has increased by 26%. For the period 1992–2015, NEP increased 93%. The increase in mean annual temperature and growing season length alone accounted for ~30% of the increase in productivity. Interannual variations in GPP and NEP were also correlated with increases in red oak biomass, forest leaf area, and canopy-scale light-use efficiency. Compared to long-term global change experiments at the Harvard Forest, the C sink in regrowing biomass equaled or exceeded C cycle modifications imposed by soil warming, N saturation, and hemlock removal. Results of this synthesis and comparison to simulation models suggest that forests across the region are likely to accrue C for decades to come but may be disrupted if the frequency or severity of biotic and abiotic disturbances increases.
This paper summarizes the open community conventions developed by the Ecological Forecasting Initiative (EFI) for the common formatting and archiving of ecological forecasts and the metadata ...associated with these forecasts. Such open standards are intended to promote interoperability and facilitate forecast communication, distribution, validation, and synthesis. For output files, we first describe the convention conceptually in terms of global attributes, forecast dimensions, forecasted variables, and ancillary indicator variables. We then illustrate the application of this convention to the two file formats that are currently preferred by the EFI, netCDF (network common data form), and comma-separated values (CSV), but note that the convention is extensible to future formats. For metadata, EFI's convention identifies a subset of conventional metadata variables that are required (e.g., temporal resolution and output variables) but focuses on developing a framework for storing information about forecast uncertainty propagation, data assimilation, and model complexity, which aims to facilitate cross-forecast synthesis. The initial application of this convention expands upon the Ecological Metadata Language (EML), a commonly used metadata standard in ecology. To facilitate community adoption, we also provide a Github repository containing a metadata validator tool and several vignettes in R and Python on how to both write and read in the EFI standard. Lastly, we provide guidance on forecast archiving, making an important distinction between short-term dissemination and long-term forecast archiving, while also touching on the archiving of code and workflows. Overall, the EFI convention is a living document that can continue to evolve over time through an open community process.
Our limited ability to accurately simulate leaf phenology is a leading source of uncertainty in models of ecosystem carbon cycling. We evaluate if continuously updating canopy state variables with ...observations is beneficial for predicting phenological events. We employed ensemble adjustment Kalman filter (EAKF) to update predictions of leaf area index (LAI) and leaf extension using tower-based photosynthetically active radiation (PAR) and moderate resolution imaging spectrometer (MODIS) data for 2002-2005 at Willow Creek, Wisconsin, USA, a mature, even-aged, northern hardwood, deciduous forest. The ecosystem demography model version 2 (ED2) was used as the prediction model, forced by offline climate data. EAKF successfully incorporated information from both the observations and model predictions weighted by their respective uncertainties. The resulting estimate reproduced the observed leaf phenological cycle in the spring and the fall better than a parametric model prediction. These results indicate that during spring the observations contribute most in determining the correct bud-burst date, after which the model performs well, but accurately modeling fall leaf senesce requires continuous model updating from observations. While the predicted net ecosystem exchange (NEE) of CO
2
precedes tower observations and unassimilated model predictions in the spring, overall the prediction follows observed NEE better than the model alone. Our results show state data assimilation successfully simulates the evolution of plant leaf phenology and improves model predictions of forest NEE.
Canopy radiative transfer is the primary mechanism by which models relate vegetation composition and state to the surface energy balance, which is important to light- and temperature-sensitive plant ...processes as well as understanding land–atmosphere feedbacks. In addition, certain parameters (e.g., specific leaf area, SLA) that have an outsized influence on vegetation model behavior can be constrained by observations of shortwave reflectance, thus reducing model predictive uncertainty. Importantly, calibrating against radiative transfer outputs allows models to directly use remote sensing reflectance products without relying on highly derived products (such as MODIS leaf area index) whose assumptions may be incompatible with the target vegetation model and whose uncertainties are usually not well quantified. Here, we created the EDR model by coupling the two-stream representation of canopy radiative transfer in the Ecosystem Demography model version 2 (ED2) with a leaf radiative transfer model (PROSPECT-5) and a simple soil reflectance model to predict full-range, high-spectral-resolution surface reflectance that is dependent on the underlying ED2 model state. We then calibrated this model against estimates of hemispherical reflectance (corrected for directional effects) from the NASA Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) and survey data from 54 temperate forest plots in the northeastern United States. The calibration significantly reduced uncertainty in model parameters related to leaf biochemistry and morphology and canopy structure for five plant functional types. Using a single common set of parameters across all sites, the calibrated model was able to accurately reproduce surface reflectance for sites with highly varied forest composition and structure. However, the calibrated model's predictions of leaf area index (LAI) were less robust, capturing only 46 % of the variability in the observations. Comparing the ED2 radiative transfer model with another two-stream soil–leaf–canopy radiative transfer model commonly used in remote sensing studies (PRO4SAIL) illustrated structural errors in the ED2 representation of direct radiation backscatter that resulted in systematic underestimation of reflectance. In addition, we also highlight that, to directly compare with a two-stream radiative transfer model like EDR, we had to perform an additional processing step to convert the directional reflectance estimates of AVIRIS to hemispherical reflectance (also known as “albedo”). In future work, we recommend that vegetation models add the capability to predict directional reflectance, to allow them to more directly assimilate a wide range of airborne and satellite reflectance products. We ultimately conclude that despite these challenges, using dynamic vegetation models to predict surface reflectance is a promising avenue for model calibration and validation using remote sensing data.