skip to main content

Title: The geodiv r package: Tools for calculating gradient surface metrics
1. The geodiv r package calculates gradient surface metrics from imagery and other gridded datasets to provide continuous measures of landscape heterogeneity for landscape pattern analysis. 2. geodiv is the first open-source, command line toolbox for calculating many gradient surface metrics and easily integrates parallel computing for applications with large images or rasters (e.g. remotely sensed data). All functions may be applied either globally to derive a single metric for an entire image or locally to create a texture image over moving windows of a user-defined extent. 3. We present a comprehensive description of the functions available through geodiv. A supplemental vignette provides an example application of geodiv to the fields of landscape ecology and biogeography. 4. geodiv allows users to easily retrieve estimates of spatial heterogeneity for a variety of purposes, enhancing our understanding of how environmental structure influences ecosystem processes. The package works with any continuous imagery and may be widely applied in many fields where estimates of surface complexity are useful.
; ; ; ; ;
Goslee, Sarah
Award ID(s):
1926568 1926567
Publication Date:
Journal Name:
Methods in Ecology and Evolution
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. The research data repository of the Environmental Data Initiative (EDI) is building on over 30 years of data curation research and experience in the National Science Foundation-funded US Long-Term Ecological Research (LTER) Network. It provides mature functionalities, well established workflows, and now publishes all ‘long-tail’ environmental data. High quality scientific metadata are enforced through automatic checks against community developed rules and the Ecological Metadata Language (EML) standard. Although the EDI repository is far along in making its data findable, accessible, interoperable, and reusable (FAIR), representatives from EDI and the LTER are developing best practices for the edge cases in environmental data publishing. One of these is the vast amount of imagery taken in the context of ecological research, ranging from wildlife camera traps to plankton imaging systems to aerial photography. Many images are used in biodiversity research for community analyses (e.g., individual counts, species cover, biovolume, productivity), while others are taken to study animal behavior and landscape-level change. Some examples from the LTER Network include: using photos of a heron colony to measure provisioning rates for chicks (Clarkson and Erwin 2018) or identifying changes in plant cover and functional type through time (Peters et al. 2020). Multi-spectral images are employedmore »to identify prairie species. Underwater photo quads are used to monitor changes in benthic biodiversity (Edmunds 2015). Sosik et al. (2020) used a continuous Imaging FlowCytobot to identify and measure phyto- and microzooplankton. Cameras at McMurdo Dry Valleys assess snow and ice cover on Antarctic lakes allowing estimation of primary production (Myers 2019). It has been standard practice to publish numerical data extracted from images in EDI; however, the supporting imagery generally has not been made publicly available. Our goal in developing best practices for documenting and archiving these images is for them to be discovered and re-used. Our examples demonstrate several issues. The research questions, and hence, the image subjects are variable. Images frequently come in logical sets of time series. The size of such sets can be large and only some images may be contributed to a dedicated specialized repository. Finally, these images are taken in a larger monitoring context where many other environmental data are collected at the same time and location. Currently, a typical approach to publishing image data in EDI are packages containing compressed (ZIP or tar) files with the images, a directory manifest with additional image-specific metadata, and a package-level EML metadata file. Images in the compressed archive may be organized within directories with filenames corresponding to treatments, locations, time periods, individuals, or other grouping attributes. Additionally, the directory manifest table has columns for each attribute. Package-level metadata include standard coverage elements (e.g., date, time, location) and sampling methods. This approach of archiving logical ‘sets’ of images reduces the effort of providing metadata for each image when most information would be repeated, but at the expense of not making every image individually searchable. The latter may be overcome if the provided manifest contains standard metadata that would allow searching and automatic integration with other images.« less
  2. Background. An assumption of Digital Image Correlation (DIC) is that the displacement field within each subset is relatively smooth, captured with reasonable accuracy by, for example, linear or quadratic shape functions. Although this assumption works well for many materials, it becomes problematic for heterogeneous materials, such as fiber networks, wherein the length scale of heterogeneity matches the size of a subset. Objective. Here we applied DIC to fibrous networks made of collagen, for which displacements at the scale of a subset are highly heterogeneous, but errors caused by the heterogeneity are difficult to quantify. We developed a method to quantify such errors. Methods. We began by generating a synthetic three-dimensional fiber network with structure matching that of gels made of fibrous collagen. We then formulated an algorithm to mimic the way in which a confocal microscope images the fibers at its focal plane, thereby generating synthetic images similar to those obtained in experiments. Displacement boundary conditions were applied to the synthetic fiber networks, and the resulting displacement fields were computed using a finite element solver. DIC was applied to the synthetic images, and displacements were compared to the data from the finite element method, enabling rigorous quantification of error. Results.more »Point-wise errors in the DIC-measured displacements were substantial, often exceeding 40%, but over regions far larger than the length scales of heterogeneity or the DIC subset size, errors were modest, e.g., ≤15%. Conclusions. Although DIC can accurately measure displacements of fiber networks at length scales larger than the subset window, quantification of mechanical behavior at the scale of material heterogeneity will require new methods to complement or replace the use of DIC.« less
  3. Abstract
    Excessive phosphorus (P) applications to croplands can contribute to eutrophication of surface waters through surface runoff and subsurface (leaching) losses. We analyzed leaching losses of total dissolved P (TDP) from no-till corn, hybrid poplar (Populus nigra X P. maximowiczii), switchgrass (Panicum virgatum), miscanthus (Miscanthus giganteus), native grasses, and restored prairie, all planted in 2008 on former cropland in Michigan, USA. All crops except corn (13 kg P ha−1 year−1) were grown without P fertilization. Biomass was harvested at the end of each growing season except for poplar. Soil water at 1.2 m depth was sampled weekly to biweekly for TDP determination during March–November 2009–2016 using tension lysimeters. Soil test P (0–25 cm depth) was measured every autumn. Soil water TDP concentrations were usually below levels where eutrophication of surface waters is frequently observed (> 0.02 mg L−1) but often higher than in deep groundwater or nearby streams and lakes. Rates of P leaching, estimated from measured concentrations and modeled drainage, did not differ statistically among cropping systems across years; 7-year cropping system means ranged from 0.035 to 0.072 kg P ha−1 year−1 with large interannual variation. Leached P was positively related to STP, which decreased over the 7 years in all systems. These results indicate that both P-fertilized and unfertilized cropping systems mayMore>>
  4. Abstract. Permafrost soils store between 1330 and 1580 Pg carbon (C), which is 3 times the amount of C in global vegetation, almost twice the amount of C in the atmosphere, and half of the global soil organic C pool. Despite the massive amount of C in permafrost, estimates of soil C storage in the high-latitude permafrost region are highly uncertain, primarily due to undersampling at all spatial scales; circumpolar soil C estimates lack sufficient continental spatial diversity, regional intensity, and replication at the field-site level. Siberian forests are particularly undersampled, yet the larch forests that dominate this region may store more than twice as much soil C as all other boreal forest types in the continuous permafrost zone combined. Here we present above- and belowground C stocks from 20 sites representing a gradient of stand age and structure in a larch watershed of the Kolyma River, near Chersky, Sakha Republic, Russia. We found that the majority of C stored in the top 1 m of the watershed was stored belowground (92 %), with 19 % in the top 10 cm of soil and 40 % in the top 30 cm. Carbon was more variable in surface soils (10 cm; coefficient of variation (CV)  =  0.35 between stands) than inmore »the top 30 cm (CV  =  0.14) or soil profile to 1 m (CV  =  0.20). Combined active-layer and deep frozen deposits (surface – 15 m) contained 205 kg C m−2 (yedoma, non-ice wedge) and 331 kg C m−2 (alas), which, even when accounting for landscape-level ice content, is an order of magnitude more C than that stored in the top meter of soil and 2 orders of magnitude more C than in aboveground biomass. Aboveground biomass was composed of primarily larch (53 %) but also included understory vegetation (30 %), woody debris (11 %) and snag (6 %) biomass. While aboveground biomass contained relatively little (8 %) of the C stocks in the watershed, aboveground processes were linked to thaw depth and belowground C storage. Thaw depth was negatively related to stand age, and soil C density (top 10 cm) was positively related to soil moisture and negatively related to moss and lichen cover. These results suggest that, as the climate warms, changes in stand age and structure may be as important as direct climate effects on belowground environmental conditions and permafrost C vulnerability.« less
  5. Record, Sydne (Ed.)
    1. LiDAR data are being increasingly used to provide a detailed characterization of the vertical profile of forests. This characterization enables the generation of new insights on the influence of environmental drivers and anthropogenic disturbances on forest structure as well as on how forest structure influences important ecosystem functions and services. Unfortunately, extracting information from LiDAR data in a way that enables the spatial visualization of forest structure, as well as its temporal changes, is challenging due to the high-dimensionality of these data. 2. We show how the Latent Dirichlet Allocation model applied to LiDAR data (LidarLDA) can be used to identify forest structural types and how the relative abundance of these forest types changes throughout the landscape. The code to fit this model is made available through the open-source R package LidarLDA in github. We illustrate the use of LidarLDA both with simulated data and data from a large-scale fire experiment in the Brazilian Amazon region. 3. Using simulated data, we demonstrate that LidarLDA accurately identifies the number of forest types as well as their spatial distribution and absorptance probabilities. For the empirical data, we found that LidarLDA detects both landscape-level patterns in forest structure as well as themore »strong interacting effect of fire and forest fragmentation on forest structure based on the experimental fire plots. More specifically, LidarLDA reveals that proximity to forest edge exacerbates the impact of fires, and that burned forests remain structurally different from unburned areas for at least seven years, even when burned only once. Importantly, LidarLDA generates insights on the 3D structure of forest that cannot be obtained using more standard approaches that just focus on top-of-the-canopy information (e.g., canopy height models based on LiDAR data). 4. By enabling the mapping of forest structure and its temporal changes, we believe that LidarLDA will be of broad utility to the ecological research community.« less