An official website of the United States government
Abstract Ecosystem models offer a rigorous way to formalize scientific theories and are critical to evaluating complex interactions among ecological and biogeochemical processes. In addition to simulation and prediction, ecosystem models are a valuable tool for testing hypotheses about mechanisms and empirical findings because they reveal critical internal processes that are difficult to observe directly.However, many ecosystem models are difficult to manage and apply by scientists who lack advanced computing skills due to complex model structures, lack of consistent documentation, and low-level programming implementation, which facilitates computing but reduces accessibility.Here, we present the ‘pnetr’ R package, which is designed to provide an easy-to-manage ecosystem modeling framework and detailed documentation in both model structure and programming. The framework implements a family of widely used PnET (net photosynthesis, evapotranspiration) ecosystem models, which are relatively parsimonious but capture essential biogeochemical cycles of water, carbon, and nutrients. We chose the R programming language since it is familiar to many ecologists and has abundant statistical modeling resources. We showcase examples of model simulations and test the effects of phenology on carbon assimilation and wood production using data measured by the Environmental Measurement Station (EMS) eddy-covariance flux tower at Harvard Forest, MA.We hope ‘pnetr’ can facilitate further development of ecological theory and increase the accessibility of ecosystem modeling and ecological forecasting.
more »
« less
pnetr : An R package for the PnET family of forest ecosystem models
Abstract Ecosystem models offer a rigorous way to formalize scientific theories and are critical to evaluating complex interactions among ecological and biogeochemical processes. In addition to simulation and prediction, ecosystem models are a valuable tool for testing hypotheses about mechanisms and empirical findings because they reveal critical internal processes that are difficult to observe directly.However, many ecosystem models are difficult to manage and apply by scientists because of complex model structures, lack of consistent documentation, and low‐level programming implementation.Here, we present the ‘pnetr’ R package, which is designed to provide an easy‐to‐manage ecosystem modelling framework and detailed documentation in both model structure and programming. The framework implements a family of widely used PnET (net photosynthesis, evapotranspiration) ecosystem models, which are relatively parsimonious but capture essential biogeochemical cycles of water, carbon and nitrogen. We chose the R programming language because it is familiar to many ecologists and has abundant statistical modelling resources. We showcase examples of model simulations and test the effects of phenology on carbon assimilation and wood production using data measured by the Environmental Measurement Station (EMS) eddy‐covariance flux tower at Harvard Forest, MA.We hope ‘pnetr’ can facilitate further development of ecological theory and increase the accessibility of ecosystem modelling and ecological forecasting.
more »
« less
- Award ID(s):
- 2205705
- PAR ID:
- 10599202
- Publisher / Repository:
- Wiley-Blackwell
- Date Published:
- Journal Name:
- Methods in Ecology and Evolution
- Volume:
- 16
- Issue:
- 7
- ISSN:
- 2041-210X
- Format(s):
- Medium: X Size: p. 1378-1388
- Size(s):
- p. 1378-1388
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Accurate quantification of soil carbon fluxes is essential to reduce uncertainty in estimates of the terrestrial carbon sink. However, these fluxes vary over time and across ecosystem types and so, it can be difficult to estimate them accurately across large scales. The flux‐gradient method estimates soil carbon fluxes using co‐located measurements of soil CO2concentration, soil temperature, soil moisture and other soil properties. The National Ecological Observatory Network (NEON) provides such data across 20 ecoclimatic domains spanning the continental U.S., Puerto Rico, Alaska and Hawai‘i.We present an R software package (neonSoilFlux) that acquires soil environmental data to compute half‐hourly soil carbon fluxes for each soil replicate plot at a given terrestrial NEON site. To assess the computed fluxes, we visited six focal NEON sites and measured soil carbon fluxes using a closed‐dynamic chamber approach.Outputs from theneonSoilFluxshowed agreement with measured fluxes (R2between measured andneonSoilFluxoutputs ranging from 0.12 to 0.77 depending on calculation method used); measured outputs generally fell within the range of calculated uncertainties from the gradient method. Calculated fluxes fromneonSoilFluxaggregated to the daily scale exhibited expected site‐specific seasonal patterns.While the flux‐gradient method is broadly effective, its accuracy is highly sensitive to site‐specific inputs, including the extent to which gap‐filing techniques are used to interpolate missing sensor data and to estimates of soil diffusivity and moisture content. Future refinement and validation ofneonSoilFluxoutputs can contribute to existing databases of soil carbon flux measurements, providing near real‐time estimates of a critical component of the terrestrial carbon cycle.more » « less
-
For everything there is a season: Analysing periodic mortality patterns with the cyclomort r packageAbstract Many important demographic processes are seasonal, including survival. For many species, mortality risk is significantly higher at certain times of the year than at others, whether because resources are scarce, susceptibility to predators or disease is high, or both. Despite the importance of survival modelling in wildlife sciences, no tools are available to estimate the peak, duration and relative importance of these ‘seasons of mortality’.We presentcyclomort, anrpackage that estimates the timing, duration and intensity of any number of mortality seasons with reliable confidence intervals. The package includes a model selection approach to determine the number of mortality seasons and to test whether seasons of mortality vary across discrete grouping factors.We illustrate the periodic hazard function model and workflow of cyclomort with simulated data. We then estimate mortality seasons of two caribouRangifer taranduspopulations that have strikingly different mortality patterns, including different numbers and timing of mortality peaks, and a marked change in one population over time.Thecyclomortpackage was developed to estimate mortality seasons for wildlife, but the package can model any time‐to‐event processes with a periodic component.more » « less
-
Abstract A prominent challenge for managing migratory species is the development of conservation plans that accommodate spatiotemporally varying distributions throughout the year. Migratory networks are spatially‐explicit models that incorporate migratory assignment and seasonal abundance data to define patterns of connectivity between stages of the annual cycle. These models are particularly useful for widespread application because different types of migratory data can be used to quantify individual and population‐level movement across the annual cycle of migratory species. While there are clear benefits of combining migratory assignment and abundance data for the development of conservation strategies, there is a concurrent need for corresponding user‐friendly software to facilitate the integration of these data for conservation.Here, we presentmignette(migratory network tools ensemble), an R package for developing migratory network models to estimate network connectivity among migratory populations. We demonstrate the functionality ofmignettewith three empirical examples that highlight the use of different types of tracking data for migratory assignment.mignettefacilitates the modelling of migratory networks by providing R functions to: (1) define breeding and nonbreeding nodes, (2) assemble abundance and assignment data and (3) model the migratory network. Additionally,mignetteprovides R functions to visualize modelled migratory networks.With increasing availability of migratory assignment and abundance data,mignetterepresents a valuable tool for developing effective conservation strategies for migratory species.more » « less
-
Abstract Given the current rates of climate change, with associated shifts in herbivore population densities, understanding the role of different herbivores in ecosystem functioning is critical for predicting ecosystem responses. Here, we examined how migratory geese and resident, non‐migratory reindeer—two dominating yet functionally contrasting herbivores—control vegetation and ecosystem processes in rapidly warming Arctic tundra.We collected vegetation and ecosystem carbon (C) flux data at peak plant growing season in the two longest running, fully replicated herbivore removal experiments found in high‐Arctic Svalbard. Experiments had been set up independently in wet habitat utilised by barnacle geeseBranta leucopsisin summer and in moist‐to‐dry habitat utilised by wild reindeerRangifer tarandus platyrhynchusyear‐round.Excluding geese induced vegetation state transitions from heavily grazed, moss‐dominated (only 4 g m−2of live above‐ground vascular plant biomass) to ungrazed, graminoid‐dominated (60 g m−2after 4‐year exclusion) and horsetail‐dominated (150 g m−2after 15‐year exclusion) tundra. This caused large increases in vegetation C and nitrogen (N) pools, dead biomass and moss‐layer depth. Alterations in plant N concentration and CN ratio suggest overall slower plant community nutrient dynamics in the short‐term (4‐year) absence of geese. Long‐term (15‐year) goose removal quadrupled net ecosystem C sequestration (NEE) by increasing ecosystem photosynthesis more than ecosystem respiration (ER).Excluding reindeer for 21 years also produced detectable increases in live above‐ground vascular plant biomass (from 50 to 80 g m−2; without promoting vegetation state shifts), as well as in vegetation C and N pools, dead biomass, moss‐layer depth and ER. Yet, reindeer removal did not alter the chemistry of plants and soil or NEE.Synthesis. Although both herbivores were key drivers of ecosystem structure and function, the control exerted by geese in their main habitat (wet tundra) was much more pronounced than that exerted by reindeer in their main habitat (moist‐to‐dry tundra). Importantly, these herbivore effects are scale dependent, because geese are more spatially concentrated and thereby affect a smaller portion of the tundra landscape compared to reindeer. Our results highlight the substantial heterogeneity in how herbivores shape tundra vegetation and ecosystem processes, with implications for ongoing environmental change.more » « less
