Showing 1–41 of 41 results for author: Finley, A O

Spatio-temporal areal models to support small area estimation: An application to national-scale forest carbon monitoring

Authors: Elliot S. Shannon, Andrew O. Finley, Paul B. May, Grant M. Domke, Hans-Erik Andersen, George C. Gaines III, Sudipto Banerjee

Abstract: National Forest Inventory (NFI) programs can provide vital information on the status, trend, and change in forest parameters. These programs are being increasingly asked to provide forest parameter estimates for spatial and temporal extents smaller than their current design and accompanying design-based methods can deliver with desired levels of uncertainty. Many NFI designs and estimation methods… ▽ More National Forest Inventory (NFI) programs can provide vital information on the status, trend, and change in forest parameters. These programs are being increasingly asked to provide forest parameter estimates for spatial and temporal extents smaller than their current design and accompanying design-based methods can deliver with desired levels of uncertainty. Many NFI designs and estimation methods focus on status and are not well equipped to provide acceptable estimates for trend and change parameters, especially over small spatial domains and/or short time periods. Fine-scale space-time indexed estimates are critical to a variety of environmental, ecological, and economic monitoring efforts. Estimates for forest carbon status, trend, and change are of particular importance to international initiatives to track carbon dynamics. Model-based small area estimation (SAE) methods for NFI and similar ecological monitoring data typically pursue inference on status within small spatial domains, with few demonstrated methods that account for spatio-temporal dependence needed for trend and change estimation. We propose a spatio-temporal Bayesian model framework that delivers statistically valid estimates with full uncertainty quantification for status, trend, and change. The framework accommodates a variety of space and time dependency structures, and we detail model configurations for different settings. Through analysis of simulated datasets, we compare the relative performance of candidate models and a traditional direct estimator. We then apply candidate models to a large-scale NFI dataset to demonstrate the utility of the proposed framework for providing unique quantification of forest carbon dynamics in the contiguous United States. We also provide computationally efficient algorithms, software, and data to reproduce our results and for benchmarking. △ Less

Submitted 30 July, 2024; v1 submitted 13 July, 2024; originally announced July 2024.

Calibrating satellite maps with field data for improved predictions of forest biomass

Authors: Paul B. May, Andrew O. Finley

Abstract: Spatially explicit quantification of forest biomass is important for forest-health monitoring and carbon accounting. Direct field measurements of biomass are laborious and expensive, typically limiting their spatial and temporal sampling density and therefore the precision and resolution of the resulting inference. Satellites can provide biomass predictions at a far greater density, but these pred… ▽ More Spatially explicit quantification of forest biomass is important for forest-health monitoring and carbon accounting. Direct field measurements of biomass are laborious and expensive, typically limiting their spatial and temporal sampling density and therefore the precision and resolution of the resulting inference. Satellites can provide biomass predictions at a far greater density, but these predictions are often biased relative to field measurements and exhibit heterogeneous errors. We developed and implemented a coregionalization model between sparse field measurements and a predictive satellite map to deliver improved predictions of biomass density at a 1-by-1 km resolution throughout the Pacific states of California, Oregon and Washington. The model accounts for zero-inflation in the field measurements and the heterogeneous errors in the satellite predictions. A stochastic partial differential equation approach to spatial modeling is applied to handle the magnitude of the satellite data. The spatial detail rendered by the model is much finer than would be possible with the field measurements alone, and the model provides substantial noise-filtering and bias-correction to the satellite map. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Model-assisted estimation of domain totals, areas, and densities in two-stage sample survey designs

Authors: Hans-Erik Andersen, Göran Ståhl, Bruce D. Cook, Douglas C. Morton, Andrew O. Finley

Abstract: Model-assisted, two-stage forest survey sampling designs provide a means to combine airborne remote sensing data, collected in a sampling mode, with field plot data to increase the precision of national forest inventory estimates, while maintaining important properties of design-based inventories, such as unbiased estimation and quantification of uncertainty. In this study, we present a comprehens… ▽ More Model-assisted, two-stage forest survey sampling designs provide a means to combine airborne remote sensing data, collected in a sampling mode, with field plot data to increase the precision of national forest inventory estimates, while maintaining important properties of design-based inventories, such as unbiased estimation and quantification of uncertainty. In this study, we present a comprehensive set of model-assisted estimators for domain-level attributes in a two-stage sampling design, including new estimators for densities, and compare the performance of these estimators with standard poststratified estimators. Simulation was used to assess the statistical properties (bias, variability) of these estimators, with both simple random and systematic sampling configurations, and indicated that 1) all estimators were generally unbiased. and 2) the use of lidar in a sampling mode increased the precision of the estimators at all assessed field sampling intensities, with particularly marked increases in precision at lower field sampling intensities. Variance estimators are generally unbiased for model-assisted estimators without poststratification, while model-assisted estimators with poststratification were increasingly biased as field sampling intensity decreased. In general, these results indicate that airborne remote sensing, collected in a sampling mode, can be used to increase the efficiency of national forest inventories. △ Less

Submitted 16 February, 2024; originally announced February 2024.

A spatial mixture model for spaceborne lidar observations over mixed forest and non-forest land types

Authors: Paul B. May, Andrew O. Finley, Ralph O. Dubayah

Abstract: The Global Ecosystem Dynamics Investigation (GEDI) is a spaceborne lidar instrument that collects near-global measurements of forest structure. While expansive in scope, GEDI samples are spatially sparse and cover a small fraction of the land surface. Converting the sparse samples into spatially complete predictive maps is of practical importance for many ecological studies. A complicating factor… ▽ More The Global Ecosystem Dynamics Investigation (GEDI) is a spaceborne lidar instrument that collects near-global measurements of forest structure. While expansive in scope, GEDI samples are spatially sparse and cover a small fraction of the land surface. Converting the sparse samples into spatially complete predictive maps is of practical importance for many ecological studies. A complicating factor is that GEDI collects measurements over forested and non-forested land alike, with no automatic labeling of the land type. Such classification is important, as it categorically influences the probability distribution of the spatial process and the ecological interpretation of the observations/predictions. We implement a spatial mixture model, separating the spatial domain into two latent classes. The latent classes are governed by a Bernoulli spatial process and within each class the process is governed by a separate spatial model. Model predictions take the form of scalar predictions as well as discrete labeling of the class membership. Inference is conducted through a Bayesian paradigm, yielding rich quantification of prediction and uncertainty. We demonstrate the method using GEDI data over Wollemi National Park. When compared to a single spatial model, the mixture model achieves much higher posterior predictive densities on the true value. When compared to a random forest model, a common algorithmic approach in the remote sensing community, the random forest achieves better absolute prediction accuracy for prediction locations far from observed training data locations, but at the expense of location-specific assessments of uncertainty. The unsupervised binary classifications of the mixture model appear broadly ecologically interpretable as forest and non-forest when compared to optical imagery, but further comparison to ground-truth data is required. △ Less

Submitted 3 January, 2024; originally announced January 2024.

Bayesian Modeling of Incompatible Spatial Data: A Case Study Involving Post-Adrian Storm Forest Damage Assessment

Authors: Lu Zhang, Andrew O. Finley, Arne Nothdurft, Sudipto Banerjee

Abstract: Modeling incompatible spatial data, i.e., data with different spatial resolutions, is a pervasive challenge in remote sensing data analysis. Typical approaches to addressing this challenge aggregate information to a common coarse resolution, i.e., compatible resolutions, prior to modeling. Such pre-processing aggregation simplifies analysis, but potentially causes information loss and hence compro… ▽ More Modeling incompatible spatial data, i.e., data with different spatial resolutions, is a pervasive challenge in remote sensing data analysis. Typical approaches to addressing this challenge aggregate information to a common coarse resolution, i.e., compatible resolutions, prior to modeling. Such pre-processing aggregation simplifies analysis, but potentially causes information loss and hence compromised inference and predictive performance. To avoid losing potential information provided by finer spatial resolution data and improve predictive performance, we propose a new Bayesian method that constructs a latent spatial process model at the finest spatial resolution. This model is tailored to settings where the outcome variable is measured on a coarser spatial resolution than predictor variables -- a configuration seen increasingly when high spatial resolution remotely sensed predictors are used in analysis. A key contribution of this work is an efficient algorithm that enables full Bayesian inference using finer resolution data while optimizing computational and storage costs. The proposed method is applied to a forest damage assessment for the 2018 Adrian storm in Carinthia, Austria, that uses high-resolution laser imaging detection and ranging (LiDAR) measurements and relatively coarse resolution forest inventory measurements. Extensive simulation studies demonstrate the proposed approach substantially improves inference for small prediction units. △ Less

Submitted 15 August, 2024; v1 submitted 19 November, 2023; originally announced November 2023.

Comments: 27 pages, 8 figures

Spatial prediction of diameter distributions for the alpine protection forests in Ebensee, Austria, using ALS/PLS and spatial distributional regression models

Authors: Arne Nothdurft, Andreas Tockner, Sarah Witzmann, Christoph Gollob, Tim Ritter, Ralf Kraßnitzer, Karl Stampfer, Andrew O. Finley

Abstract: A spatial distributional regression model is presented to predict the forest structural diversity in terms of the distributions of the stem diameter at breast height (DBH) in the protection forests in Ebensee, Austria. In total 36,338 sample trees were measured via a handheld mobile personal laser scanning system (PLS) on 273 sample plots each having a 20 m radius. Recent airborne laser scanning (… ▽ More A spatial distributional regression model is presented to predict the forest structural diversity in terms of the distributions of the stem diameter at breast height (DBH) in the protection forests in Ebensee, Austria. In total 36,338 sample trees were measured via a handheld mobile personal laser scanning system (PLS) on 273 sample plots each having a 20 m radius. Recent airborne laser scanning (ALS) data was used to derive regression covariates from the normalized digital vegetation height model (DVHM) and the digital terrain model (DTM). Candidate models were constructed that differed in their linear predictors of the two gamma distribution parameters. Non-linear smoothing splines outperformed linear parametric slope coefficients, and the best implementation of spatial structured effects was achieved by a Gaussian process smooth. Model fitting and posterior parameter inference was achieved by using full Bayesian methodology and MCMC sampling algorithms implemented in the R-package BAMLSS. Spatial predictions of stem count proportions per DBH classes revealed that regeneration of smaller trees was lacking in certain areas of the protection forest landscape. △ Less

Submitted 3 November, 2023; originally announced November 2023.

spAbundance: An R package for single-species and multi-species spatially explicit abundance models

Authors: Jeffrey W. Doser, Andrew O. Finley, Marc Kéry, Elise F. Zipkin

Abstract: Numerous modeling techniques exist to estimate abundance of plant and wildlife species. These methods seek to estimate abundance while accounting for multiple complexities found in ecological data, such as observational biases, spatial autocorrelation, and species correlations. There is, however, a lack of user-friendly and computationally efficient software to implement the various models, partic… ▽ More Numerous modeling techniques exist to estimate abundance of plant and wildlife species. These methods seek to estimate abundance while accounting for multiple complexities found in ecological data, such as observational biases, spatial autocorrelation, and species correlations. There is, however, a lack of user-friendly and computationally efficient software to implement the various models, particularly for large data sets. We developed the spAbundance R package for fitting spatially-explicit Bayesian single-species and multi-species hierarchical distance sampling models, N-mixture models, and generalized linear mixed models. The models within the package can account for spatial autocorrelation using Nearest Neighbor Gaussian Processes and accommodate species correlations in multi-species models using a latent factor approach, which enables model fitting for data sets with large numbers of sites and/or species. We provide three vignettes and three case studies that highlight spAbundance functionality. We used spatially-explicit multi-species distance sampling models to estimate density of 16 bird species in Florida, USA, an N-mixture model to estimate Black-throated Blue Warbler (Setophaga caerulescens) abundance in New Hampshire, USA, and a spatial linear mixed model to estimate forest aboveground biomass across the continental USA. spAbundance provides a user-friendly, formula-based interface to fit a variety of univariate and multivariate spatially-explicit abundance models. The package serves as a useful tool for ecologists and conservation practitioners to generate improved inference and predictions on the spatial drivers of populations and communities. △ Less

Submitted 13 May, 2024; v1 submitted 30 October, 2023; originally announced October 2023.

Journal ref: Methods in Ecology and Evolution

Modeling complex species-environment relationships through spatially-varying coefficient occupancy models

Authors: Jeffrey W. Doser, Andrew O. Finley, Sarah P. Saunders, Marc Kery, Aaron S. Weed, Elise F. Zipkin

Abstract: Occupancy models are frequently used by ecologists to quantify spatial variation in species distributions while accounting for observational biases in the collection of detection-nondetection data. However, the common assumption that a single set of regression coefficients can adequately explain species-environment relationships is often unrealistic, especially across large spatial domains. Here w… ▽ More Occupancy models are frequently used by ecologists to quantify spatial variation in species distributions while accounting for observational biases in the collection of detection-nondetection data. However, the common assumption that a single set of regression coefficients can adequately explain species-environment relationships is often unrealistic, especially across large spatial domains. Here we develop single-species (i.e., univariate) and multi-species (i.e., multivariate) spatially-varying coefficient (SVC) occupancy models to account for spatially-varying species-environment relationships. We employ Nearest Neighbor Gaussian Processes and Polya-Gamma data augmentation in a hierarchical Bayesian framework to yield computationally efficient Gibbs samplers, which we implement in the spOccupancy R package. For multi-species models, we use spatial factor dimension reduction to efficiently model datasets with large numbers of species (e.g., > 10). The hierarchical Bayesian framework readily enables generation of posterior predictive maps of the SVCs, with fully propagated uncertainty. We apply our SVC models to quantify spatial variability in the relationships between maximum breeding season temperature and occurrence probability of 21 grassland bird species across the U.S. Jointly modeling species generally outperformed single-species models, which all revealed substantial spatial variability in species occurrence relationships with maximum temperatures. Our models are particularly relevant for quantifying species-environment relationships using detection-nondetection data from large-scale monitoring programs, which are becoming increasingly prevalent for answering macroscale ecological questions regarding wildlife responses to global change. △ Less

Submitted 4 August, 2023; originally announced August 2023.

Comments: JABES (2024)

Models to support forest inventory and small area estimation using sparsely sampled LiDAR: A case study involving G-LiHT LiDAR in Tanana, Alaska

Authors: Andrew O. Finley, Hans-Erik Andersen, Chad Babcock, Bruce D. Cook, Douglas C. Morton, Sudipto Banerjee

Abstract: A two-stage hierarchical Bayesian model is developed and implemented to estimate forest biomass density and total given sparsely sampled LiDAR and georeferenced forest inventory plot measurements. The model is motivated by the United States Department of Agriculture (USDA) Forest Service Forest Inventory and Analysis (FIA) objective to provide biomass estimates for the remote Tanana Inventory Unit… ▽ More A two-stage hierarchical Bayesian model is developed and implemented to estimate forest biomass density and total given sparsely sampled LiDAR and georeferenced forest inventory plot measurements. The model is motivated by the United States Department of Agriculture (USDA) Forest Service Forest Inventory and Analysis (FIA) objective to provide biomass estimates for the remote Tanana Inventory Unit (TIU) in interior Alaska. The proposed model yields stratum-level biomass estimates for arbitrarily sized areas. Model-based estimates are compared with the TIU FIA design-based post-stratified estimates. Model-based small area estimates (SAEs) for two experimental forests within the TIU are compared with each forest's design-based estimates generated using a dense network of independent inventory plots. Model parameter estimates and biomass predictions are informed using FIA plot measurements, LiDAR data that are spatially aligned with a subset of the FIA plots, and complete coverage remotely detected data used to define landuse/landcover stratum and percent forest canopy cover. Results support a model-based approach to estimating forest variables when inventory data are sparse or resources limit collection of enough data to achieve desired accuracy and precision using design-based methods. △ Less

Submitted 31 January, 2024; v1 submitted 13 February, 2023; originally announced February 2023.

Guidelines for the use of spatially-varying coefficients in species distribution models

Authors: Jeffrey W. Doser, Marc Kéry, Sarah P. Saunders, Andrew O. Finley, Brooke L. Bateman, Joanna Grand, Shannon Reault, Aaron S. Weed, Elise F. Zipkin

Abstract: Species distribution models (SDMs) are increasingly applied across macroscales. Such models typically assume that a single set of regression coefficients can adequately describe species-environment relationships and/or population trends. However, such relationships often show nonlinear and/or spatially-varying patterns that arise from complex interactions with abiotic and biotic processes that ope… ▽ More Species distribution models (SDMs) are increasingly applied across macroscales. Such models typically assume that a single set of regression coefficients can adequately describe species-environment relationships and/or population trends. However, such relationships often show nonlinear and/or spatially-varying patterns that arise from complex interactions with abiotic and biotic processes that operate at different scales. Spatially-varying coefficient (SVC) models can readily account for variability in the effects of environmental covariates. Yet, their use in ecology is relatively scarce due to gaps in understanding the inferential benefits that SVC models can provide compared to simpler frameworks. Here we demonstrate the inferential benefits of SVC SDMs, with a particular focus on how this approach can be used to generate and test ecological hypotheses regarding the drivers of spatial variability in population trends and species-environment relationships. We illustrate the inferential benefits of SVC SDMs with simulations and two case studies: one that assesses spatially-varying trends of 51 forest bird species in the eastern US over two decades and a second that evaluates spatial variability in the effects of five decades of land cover change on Grasshopper Sparrow occurrence across the continental US. We found strong support for SVC SDMs compared to simpler alternatives in both empirical case studies. These applications display the utility of SVC SDMs to help reveal the environmental factors that drive species distributions across both local and broad scales. We conclude by discussing the potential applications of SVC SDMs in ecology and conservation. △ Less

Submitted 4 August, 2023; v1 submitted 13 January, 2023; originally announced January 2023.

Journal ref: Global Ecology and Biogeography, 33, e13814 (2024)

Quantifying and correcting geolocation error in spaceborne LiDAR forest canopy observations using high spatial accuracy ALS: A Bayesian model approach

Authors: Elliot S. Shannon, Andrew O. Finley, Daniel J. Hayes, Sylvia N. Noralez, Aaron R. Weiskittel, Bruce D. Cook, Chad Babcock

Abstract: Geolocation error in spaceborne sampling light detection and ranging (LiDAR) measurements of forest structure can compromise forest attribute estimates and degrade integration with georeferenced field measurements or other remotely sensed data. Data integration is especially problematic when geolocation error is not well quantified. We propose a general model that uses airborne laser scanning (ALS… ▽ More Geolocation error in spaceborne sampling light detection and ranging (LiDAR) measurements of forest structure can compromise forest attribute estimates and degrade integration with georeferenced field measurements or other remotely sensed data. Data integration is especially problematic when geolocation error is not well quantified. We propose a general model that uses airborne laser scanning (ALS) data to quantify and correct geolocation error in spaceborne sampling LiDAR. To illustrate the model, LiDAR data from NASA Goddard's LiDAR Hyperspectral & Thermal Imager (G-LiHT) was used with a subset of LiDAR data from NASA's Global Ecosystem Dynamics Investigation (GEDI). The model accommodates multiple canopy height metrics derived from a simulated GEDI footprint kernel using spatially coincident G-LiHT, and incorporates both additive and multiplicative mapping between the canopy height metrics generated from both datasets. A Bayesian implementation provides probabilistic uncertainty quantification in both parameter and geolocation error estimates. Results show a systematic geolocation error of 9.62 m in the southwest direction. In addition, estimated geolocation errors within GEDI footprints were highly variable, with results showing a ~0.45 probability the true footprint center is within 20 m. Estimating and correcting geolocation error via the model outlined here can help inform subsequent efforts to integrate spaceborne LiDAR data, like GEDI, with other georeferenced data. △ Less

Submitted 23 August, 2023; v1 submitted 23 September, 2022; originally announced September 2022.

Joint species distribution models with imperfect detection for high-dimensional spatial data

Authors: Jeffrey W. Doser, Andrew O. Finley, Sudipto Banerjee

Abstract: Determining spatial distributions of species and communities are key objectives of ecology and conservation. Joint species distribution models use multi-species detection-nondetection data to estimate species and community distributions. The analysis of such data is complicated by residual correlations between species, imperfect detection, and spatial autocorrelation. While methods exist to accomm… ▽ More Determining spatial distributions of species and communities are key objectives of ecology and conservation. Joint species distribution models use multi-species detection-nondetection data to estimate species and community distributions. The analysis of such data is complicated by residual correlations between species, imperfect detection, and spatial autocorrelation. While methods exist to accommodate each of these complexities, there are few examples in the literature that address and explore all three complexities simultaneously. Here we developed a spatial factor multi-species occupancy model to explicitly account for species correlations, imperfect detection, and spatial autocorrelation. The proposed model uses a spatial factor dimension reduction approach and Nearest Neighbor Gaussian Processes to ensure computational efficiency for data sets with both a large number of species (e.g., > 100) and spatial locations (e.g., 100,000). We compare the proposed model performance to five candidate models, each addressing a subset of the three complexities. We implemented the proposed and competing models in the spOccupancy software, designed to facilitate application via an accessible, well-documented, and open-source R package. Using simulations, we found ignoring the three complexities when present leads to inferior model predictive performance. Using a case study on 98 bird species across the continental US, the spatial factor multi-species occupancy model had the highest predictive performance among the candidate models. Our proposed framework, together with its implementation in spOccupancy, serves as a user-friendly tool to understand spatial variation in species distributions and biodiversity metrics while addressing common complexities in multi-species detection-nondetection data. △ Less

Submitted 13 December, 2022; v1 submitted 6 April, 2022; originally announced April 2022.

Comments: 38 pages

spOccupancy: An R package for single-species, multi-species, and integrated spatial occupancy models

Authors: Jeffrey W. Doser, Andrew O. Finley, Marc Kéry, Elise F. Zipkin

Abstract: Occupancy modeling is a common approach to assess spatial and temporal species distribution patterns, while explicitly accounting for measurement errors common in detection-nondetection data. Numerous extensions of the basic single species occupancy model exist to address dynamics, multiple species or states, interactions, false positive errors, autocorrelation, and to integrate multiple data sour… ▽ More Occupancy modeling is a common approach to assess spatial and temporal species distribution patterns, while explicitly accounting for measurement errors common in detection-nondetection data. Numerous extensions of the basic single species occupancy model exist to address dynamics, multiple species or states, interactions, false positive errors, autocorrelation, and to integrate multiple data sources. However, development of specialized and computationally efficient software to fit spatial models to large data sets is scarce or absent. We introduce the spOccupancy R package designed to fit single-species, multi-species, and integrated spatially-explicit occupancy models. Using a Bayesian framework, we leverage Pólya-Gamma data augmentation and Nearest Neighbor Gaussian Processes to ensure models are computationally efficient for potentially massive data sets. spOccupancy provides user-friendly functions for data simulation, model fitting, model validation (by posterior predictive checks), model comparison (using information criteria and k-fold cross-validation), and out-of-sample prediction. We illustrate the package's functionality via a vignette, simulated data analysis, and two bird case studies, in which we estimate occurrence of the Black-throated Green Warbler (Setophaga virens) across the eastern USA and species richness of a foliage-gleaning bird community in the Hubbard Brook Experimental Forest in New Hampshire, USA. The spOccupancy package provides a user-friendly approach to fit a variety of single and multi-species occupancy models, making it straightforward to address detection biases and spatial autocorrelation in species distribution models even for large data sets. △ Less

Submitted 3 May, 2022; v1 submitted 23 November, 2021; originally announced November 2021.

Comments: 20 pages, 2 figures

Simplifying small area estimation with rFIA: a demonstration of tools and techniques

Authors: Hunter Stanke, Andrew O. Finley, Grant M. Domke

Abstract: The United States (US) Forest Service Forest Inventory and Analysis (FIA) program operates the national forest inventory of the US. Traditionally, the FIA program has relied on sample-based approaches -- permanent plot networks and associated design-based estimators -- to estimate forest variables across large geographic areas and long periods of time. These approaches generally offer unbiased inf… ▽ More The United States (US) Forest Service Forest Inventory and Analysis (FIA) program operates the national forest inventory of the US. Traditionally, the FIA program has relied on sample-based approaches -- permanent plot networks and associated design-based estimators -- to estimate forest variables across large geographic areas and long periods of time. These approaches generally offer unbiased inference on large domains but fail to provide reliable estimates for small domains due to low sample sizes. Rising demand for small domain estimates will thus require the FIA program to adopt non-traditional estimation approaches that are capable of delivering defensible estimates of forest variables at increased spatial and temporal resolution, without the expense of collecting additional field data. In light of this challenge, the development of small area estimation (SAE) methods for FIA data has become an active and highly productive area of research. Yet, SAE methods remain difficult to apply to FIA data, due in part to the complex data structures and inventory design used by the FIA program. Thus, we argue that a new suite of estimation tools (i.e., software) will be required to accommodate shifts in demand for inference on large geographic areas and long time periods to inference on small spatial and/or temporal domains. Herein, we present rFIA, an open-source R package designed to increase the accessibility of FIA data, as one such tool. Specifically, we present two case studies chosen to demonstrate rFIA's potential to simplify the application of a broad suite of SAE methods to FIA data: (1) estimation of contemporary county-level forest carbon stocks across the conterminous US using a spatial Fay-Herriot model; and (2) temporally-explicit estimation of multi-decadal trends in merchantable wood volume in Washington County, Maine using a Bayesian mixed-effects model. △ Less

Submitted 27 February, 2022; v1 submitted 8 September, 2021; originally announced September 2021.

Estimating timber volume loss due to storm damage in Carinthia, Austria, using ALS/TLS and spatial regression models

Authors: Arne Nothdurft, Christoph Gollob, Ralf Kraßnitzer, Gernot Erber, Tim Ritter, Karl Stampfer, Andrew O. Finley

Abstract: A spatial regression model framework is presented to predict growing stock volume loss due to storm Adrian which caused heavy forest damage in the upper Gail valley in Carinthia, Austria, in October 2018. Model parameters were estimated using growing stock volume measured with a terrestrial laser scanner on 62 sample plots distributed across five sub-regions. Predictor variables were derived from… ▽ More A spatial regression model framework is presented to predict growing stock volume loss due to storm Adrian which caused heavy forest damage in the upper Gail valley in Carinthia, Austria, in October 2018. Model parameters were estimated using growing stock volume measured with a terrestrial laser scanner on 62 sample plots distributed across five sub-regions. Predictor variables were derived from high resolution vegetation height measurements collected during an airborne laser scanning campaign. Non-spatial and spatial candidate models were proposed and assessed based on fit to observed data and out-of-sample prediction. Spatial Gaussian processes associated model intercepts and regression coefficients were used to capture spatial dependence. Results show a spatially-varying coefficient model, which allowed the intercept and regression coefficients to vary spatially, yielded the best fit and prediction. Two approaches were considered for prediction over blowdown areas: 1) an areal approach that viewed each blowdown as a single prediction unit indexed by its centroid; and 2) a block approach where each blowdown was partitioned into smaller prediction units to better align with sample plots' spatial support. Joint prediction was used to acknowledge spatial dependence among block units. Results demonstrated the block approach is preferable as it mitigated change-of-support issues encountered in the areal approach. Despite the small sample size, predictions for 55% of the total 564 blowdown areas, accounting for 93% of the total loss, had a coefficient of variation less than 25%. Key advantages of the proposed regression framework are the ability to quantify uncertainty in spatial covariance parameters, propagate parameter uncertainty through to prediction, and provide statistically valid prediction point and interval estimates for individual blowdowns and collections of blowdowns. △ Less

Submitted 16 June, 2021; originally announced June 2021.

Grid-Parametrize-Split (GriPS) for Improved Scalable Inference in Spatial Big Data Analysis

Authors: Michele Peruzzi, Sudipto Banerjee, David B. Dunson, Andrew O. Finley

Abstract: Rapid advancements in spatial technologies including Geographic Information Systems (GIS) and remote sensing have generated massive amounts of spatially referenced data in a variety of scientific and data-driven industrial applications. These advancements have led to a substantial, and still expanding, literature on the modeling and analysis of spatially oriented big data. In particular, Bayesian… ▽ More Rapid advancements in spatial technologies including Geographic Information Systems (GIS) and remote sensing have generated massive amounts of spatially referenced data in a variety of scientific and data-driven industrial applications. These advancements have led to a substantial, and still expanding, literature on the modeling and analysis of spatially oriented big data. In particular, Bayesian inferences for high-dimensional spatial processes are being sought in a variety of remote-sensing applications including, but not limited to, modeling next generation Light Detection and Ranging (LiDAR) systems and other remotely sensed data. Massively scalable spatial processes, in particular Gaussian processes (GPs), are being explored extensively for the increasingly encountered big data settings. Recent developments include GPs constructed from sparse Directed Acyclic Graphs (DAGs) with a limited number of neighbors (parents) to characterize dependence across the spatial domain. The DAG can be used to devise fast algorithms for posterior sampling of the latent process, but these may exhibit pathological behavior in estimating covariance parameters. While these issues are mitigated by considering marginalized samplers that exploit the underlying sparse precision matrix, these algorithms are slower, less flexible, and oblivious of structure in the data. The current article introduces the Grid-Parametrize-Split (GriPS) approach for conducting Bayesian inference in spatially oriented big data settings by a combination of careful model construction and algorithm design to effectuate substantial improvements in MCMC convergence. We demonstrate the effectiveness of our proposed methods through simulation experiments and subsequently undertake the modeling of LiDAR outcomes and production of their predictive maps using G-LiHT and other remotely sensed variables. △ Less

Submitted 10 January, 2021; originally announced January 2021.

Integrating automated acoustic vocalization data and point count surveys for estimation of bird abundance

Authors: Jeffrey W. Doser, Andrew O. Finley, Aaron S. Weed, Elise F. Zipkin

Abstract: Monitoring wildlife abundance across space and time is an essential task to study their population dynamics and inform effective management. Acoustic recording units are a promising technology for efficiently monitoring bird populations and communities. We present an integrated modeling framework that combines high-quality but temporally sparse bird point count survey data with acoustic recordings… ▽ More Monitoring wildlife abundance across space and time is an essential task to study their population dynamics and inform effective management. Acoustic recording units are a promising technology for efficiently monitoring bird populations and communities. We present an integrated modeling framework that combines high-quality but temporally sparse bird point count survey data with acoustic recordings. Using simulations, we compare the accuracy and precision of abundance estimates using differing amounts of acoustic vocalizations obtained from a clustering algorithm, point count data, and a subset of manually validated acoustic vocalizations. We also use our modeling framework in a case study to estimate abundance of the Eastern Wood-Pewee (Contopus virens) in Vermont, U.S.A. The simulation study reveals that combining acoustic and point count data via an integrated model improves accuracy and precision of abundance estimates compared with models informed by either acoustic or point count data alone. Combining acoustic data with only a small number of point count surveys yields estimates of abundance without the need for validating any of the identified vocalizations from the acoustic data. Within our case study, the integrated models provided moderate support for a decline of the Eastern Wood-Pewee in this region. Our integrated modeling approach combines dense acoustic data with few point count surveys to deliver reliable estimates of species abundance without the need for manual identification of acoustic vocalizations or a prohibitively expensive large number of repeated point count surveys. Our proposed approach offers an efficient monitoring alternative for large spatio-temporal regions when point count data are difficult to obtain or when monitoring is focused on rare species with low detection probability. △ Less

Submitted 6 May, 2021; v1 submitted 22 November, 2020; originally announced November 2020.

Comments: 25 pages, 4 figures

Journal ref: Methods in Ecology and Evolution (2021)

A Bayesian hierarchical model to estimate land surface phenology parameters with harmonized Landsat 8 and Sentinel-2 images

Authors: Chad Babcock, Andrew O. Finley, Nathaniel Looker

Abstract: We develop a Bayesian Land Surface Phenology (LSP) model and examine its performance using Enhanced Vegetation Index (EVI) observations derived from the Harmonized Landsat Sentinel-2 (HLS) dataset. Building on previous work, we propose a double logistic function that, once couched within a Bayesian model, yields posterior distributions for all LSP parameters. We assess the efficacy of the Normal,… ▽ More We develop a Bayesian Land Surface Phenology (LSP) model and examine its performance using Enhanced Vegetation Index (EVI) observations derived from the Harmonized Landsat Sentinel-2 (HLS) dataset. Building on previous work, we propose a double logistic function that, once couched within a Bayesian model, yields posterior distributions for all LSP parameters. We assess the efficacy of the Normal, Truncated Normal, and Beta likelihoods to deliver robust LSP parameter estimates. Two case studies are presented and used to explore aspects of the proposed model. The first, conducted over forested pixels within a HLS tile, explores choice of likelihood and space-time varying HLS data availability for long-term average LSP parameter point and uncertainty estimation. The second, conducted on a small area of interest within the HLS tile on an annual time-step, further examines the impact of sample size and choice of likelihood on LSP parameter estimates. Results indicate that while the Truncated Normal and Beta likelihoods are theoretically preferable when the vegetation index is bounded, all three likelihoods performed similarly when the number of index observations is sufficiently large and values are not near the index bounds. Both case studies demonstrate how pixel-level LSP parameter posterior distributions can be used to propagate uncertainty through subsequent analysis. As a companion to this article, we provide an open-source \R package \pkg{rsBayes} and supplementary data and code used to reproduce the analysis results. The proposed model specification and software implementation delivers computationally efficient, statistically robust, and inferentially rich LSP parameter posterior distributions at the pixel-level across massive raster time series datasets. △ Less

Submitted 10 September, 2020; originally announced September 2020.

Trends in bird abundance differ among protected forests but not bird guilds

Authors: Jeffrey W. Doser, Aaron S. Weed, Elise F. Zipkin, Kathryn M. Miller, Andrew O. Finley

Abstract: Improved monitoring and associated inferential tools to efficiently identify declining bird populations, particularly of rare or sparsely distributed species, is key to informed conservation and management across large spatio-temporal regions. We assess abundance trends for 106 bird species in a network of eight national park forests located within the northeast USA from 2006-2019 using a novel hi… ▽ More Improved monitoring and associated inferential tools to efficiently identify declining bird populations, particularly of rare or sparsely distributed species, is key to informed conservation and management across large spatio-temporal regions. We assess abundance trends for 106 bird species in a network of eight national park forests located within the northeast USA from 2006-2019 using a novel hierarchical model. We develop a multi-species, multi-region removal sampling model that shares information across species and parks to enable inference on rare species and sparsely sampled parks and to evaluate the effects of local forest structure. Trends in bird abundance over time varied widely across parks, but species showed similar trends within parks. Three parks (Acadia, Marsh-Billings-Rockefeller, and Morristown) decreased in bird abundance across all species, while three parks (Saratoga, Roosevelt-Vanderbilt, and Weir-Farm) increased in abundance. Bird abundance peaked at medium levels of basal area and high levels of percent forest and forest regeneration, with percent forest having the largest effect. Variation in these effects across parks could be a result of differences in forest structural stage and diversity. Our novel hierarchical model enables estimates of abundance at the network, park, guild, and species levels. We found large variation in abundance trends across parks but not across bird guilds, suggesting that local forest condition may have a broad and consistent effect on the entire bird community within a given park. Management should target the three parks with overall decreasing trends in bird abundance to further identify what specific factors are driving observed declines across the bird community. Understanding how bird communities respond to local forest structure and other stressors is crucial for informed and lasting management. △ Less

Submitted 6 May, 2021; v1 submitted 27 August, 2020; originally announced August 2020.

Comments: 39 pages, 6 figures

Journal ref: Ecological Applications (2021)

Highly Scalable Bayesian Geostatistical Modeling via Meshed Gaussian Processes on Partitioned Domains

Authors: Michele Peruzzi, Sudipto Banerjee, Andrew O. Finley

Abstract: We introduce a class of scalable Bayesian hierarchical models for the analysis of massive geostatistical datasets. The underlying idea combines ideas on high-dimensional geostatistics by partitioning the spatial domain and modeling the regions in the partition using a sparsity-inducing directed acyclic graph (DAG). We extend the model over the DAG to a well-defined spatial process, which we call t… ▽ More We introduce a class of scalable Bayesian hierarchical models for the analysis of massive geostatistical datasets. The underlying idea combines ideas on high-dimensional geostatistics by partitioning the spatial domain and modeling the regions in the partition using a sparsity-inducing directed acyclic graph (DAG). We extend the model over the DAG to a well-defined spatial process, which we call the Meshed Gaussian Process (MGP). A major contribution is the development of a MGPs on tessellated domains, accompanied by a Gibbs sampler for the efficient recovery of spatial random effects. In particular, the cubic MGP (Q-MGP) can harness high-performance computing resources by executing all large-scale operations in parallel within the Gibbs sampler, improving mixing and computing time compared to sequential updating schemes. Unlike some existing models for large spatial data, a Q-MGP facilitates massive caching of expensive matrix operations, making it particularly apt in dealing with spatiotemporal remote-sensing data. We compare Q-MGPs with large synthetic and real world data against state-of-the-art methods. We also illustrate using Normalized Difference Vegetation Index (NDVI) data from the Serengeti park region to recover latent multivariate spatiotemporal random effects at millions of locations. The source code is available at https://github.com/mkln/meshgp. △ Less

Submitted 5 October, 2020; v1 submitted 25 March, 2020; originally announced March 2020.

High-dimensional Multivariate Geostatistics: A Bayesian Matrix-Normal Approach

Authors: Lu Zhang, Sudipto Banerjee, Andrew O. Finley

Abstract: Joint modeling of spatially-oriented dependent variables is commonplace in the environmental sciences, where scientists seek to estimate the relationships among a set of environmental outcomes accounting for dependence among these outcomes and the spatial dependence for each outcome. Such modeling is now sought for massive data sets with variables measured at a very large number of locations. Baye… ▽ More Joint modeling of spatially-oriented dependent variables is commonplace in the environmental sciences, where scientists seek to estimate the relationships among a set of environmental outcomes accounting for dependence among these outcomes and the spatial dependence for each outcome. Such modeling is now sought for massive data sets with variables measured at a very large number of locations. Bayesian inference, while attractive for accommodating uncertainties through hierarchical structures, can become computationally onerous for modeling massive spatial data sets because of its reliance on iterative estimation algorithms. This manuscript develops a conjugate Bayesian framework for analyzing multivariate spatial data using analytically tractable posterior distributions that obviate iterative algorithms. We discuss differences between modeling the multivariate response itself as a spatial process and that of modeling a latent process in a hierarchical model. We illustrate the computational and inferential benefits of these models using simulation studies and analysis of a Vegetation Index data set with spatially dependent observations numbering in the millions. △ Less

Submitted 6 December, 2020; v1 submitted 22 March, 2020; originally announced March 2020.

Comments: 22 pages, 12 figures

spNNGP R package for Nearest Neighbor Gaussian Process models

Authors: Andrew O. Finley, Abhirup Datta, Sudipto Banerjee

Abstract: This paper describes and illustrates functionality of the spNNGP R package. The package provides a suite of spatial regression models for Gaussian and non-Gaussian point-referenced outcomes that are spatially indexed. The package implements several Markov chain Monte Carlo (MCMC) and MCMC-free Nearest Neighbor Gaussian Process (NNGP) models for inference about large spatial data. Non-Gaussian outc… ▽ More This paper describes and illustrates functionality of the spNNGP R package. The package provides a suite of spatial regression models for Gaussian and non-Gaussian point-referenced outcomes that are spatially indexed. The package implements several Markov chain Monte Carlo (MCMC) and MCMC-free Nearest Neighbor Gaussian Process (NNGP) models for inference about large spatial data. Non-Gaussian outcomes are modeled using a NNGP Polya-Gamma latent variable. OpenMP parallelization options are provided to take advantage of multiprocessor systems. Package features are illustrated using simulated and real data sets. △ Less

Submitted 15 April, 2021; v1 submitted 24 January, 2020; originally announced January 2020.

ForestFit : An R package for modeling tree diameter distributions

Authors: Mahdi Teimouri, Jeffrey W. Doser, Andrew O. Finley

Abstract: Modeling the diameter distribution of trees in forest stands is a common forestry task that supports key biologically and economically relevant management decisions. The choice of model used to represent the diameter distribution and how to estimate its parameters has received much attention in the forestry literature; however, accessible software that facilitates comprehensive comparison of the m… ▽ More Modeling the diameter distribution of trees in forest stands is a common forestry task that supports key biologically and economically relevant management decisions. The choice of model used to represent the diameter distribution and how to estimate its parameters has received much attention in the forestry literature; however, accessible software that facilitates comprehensive comparison of the myriad modeling approaches is not available. To this end, we developed an R package called ForestFit that simplifies estimation of common probability distributions used to model tree diameter distributions, including the two- and three-parameter Weibull distributions, Johnson's SB distribution, Birnbaum-Saunders distribution, and finite mixture distributions. Frequentist and Bayesian techniques are provided for individual tree diameter data, as well as grouped data. Additional functionality facilitates fitting growth curves to height-diameter data. The package also provides a set of functions for computing probability distributions and simulating random realizations from common finite mixture models. △ Less

Submitted 25 November, 2019; originally announced November 2019.

Comments: 10 pages, 7 figures

rFIA: An R package for estimation of forest attributes with the Forest Inventory and Analysis Database

Authors: Hunter Stanke, Andrew O. Finley, Aaron S. Weed, Brian F. Walters, Grant M. Domke

Abstract: Forest Inventory and Analysis (FIA) is a US Department of Agriculture Forest Service program that aims to monitor changes in forests across the US. FIA hosts one of the largest ecological datasets in the world, though its complexity limits access for many potential users. rFIA is an R package designed to simplify the estimation of forest attributes using data collected by the FIA Program. Specific… ▽ More Forest Inventory and Analysis (FIA) is a US Department of Agriculture Forest Service program that aims to monitor changes in forests across the US. FIA hosts one of the largest ecological datasets in the world, though its complexity limits access for many potential users. rFIA is an R package designed to simplify the estimation of forest attributes using data collected by the FIA Program. Specifically, rFIA improves access to the spatio-temporal estimation capacity of the FIA Database via space-time indexed summaries of forest variables within user-defined population boundaries (e.g., geographic, temporal, biophysical). The package implements multiple design-based estimators, and has been validated against official estimates and sampling errors produced by the FIA Program. We demonstrate the utility of rFIA by assessing changes in abundance and mortality rates of ash (Fraxinus spp.) populations in the Lower Peninsula of Michigan following the establishment of emerald ash borer (Agrilus planipennis). △ Less

Submitted 4 February, 2020; v1 submitted 14 November, 2019; originally announced November 2019.

Assessing soundscape disturbance through hierarchical models and acoustic indices: a case study on a shelterwood logged northern Michigan forest

Authors: Jeffrey W. Doser, Andrew O. Finley, Eric P. Kasten, Stuart H. Gage

Abstract: Assessing the effects of anthropogenic disturbances on wildlife is a necessary conservation task. The soundscape is a critical habitat component for acoustically communicating organisms, but the use of the soundscape as a tool for assessing disturbance impacts has been relatively unexplored until recently. Here we present a broad modeling framework for assessing disturbance impacts on soundscapes,… ▽ More Assessing the effects of anthropogenic disturbances on wildlife is a necessary conservation task. The soundscape is a critical habitat component for acoustically communicating organisms, but the use of the soundscape as a tool for assessing disturbance impacts has been relatively unexplored until recently. Here we present a broad modeling framework for assessing disturbance impacts on soundscapes, which we apply to quantify the influence of a shelterwood logging on soundscapes in northern Michigan. Our modeling approach can be broadly applied to assess anthropogenic disturbance impacts on soundscapes. The approach accommodates inherent differences in control and treatment sites to improve inference about treatment effects, while also accounting for extraneous variables (e.g., rain) that influence acoustic indices. Recordings were obtained at 13 sites before and after a shelterwood logging. Four sites were in the logging region and nine sites served as control recordings outside the logging region. We quantify the soundscapes using common acoustic indices (Normalized Difference Soundscape Index (NDSI), Acoustic Entropy (H), Acoustic Complexity Index (ACI), Acoustic Evenness Index (AEI), Welch Power Spectral Density (PSD)) and build two hierarchical Bayesian models to quantify the changes in the soundscape over the study period. Our analysis reveals no long-lasting effects of the shelterwood logging on the soundscape diversity as measured by the NDSI, but analysis of H, AEI, and PSD suggest changes in the evenness of sounds across the frequency spectrum, indicating a potential shift in the avian species communicating in the soundscapes as a result of the logging. Acoustic recordings, in conjunction with this modeling framework, can deliver cost efficient assessment of disturbance impacts on the landscape and underlying biodiversity as represented through the soundscape. △ Less

Submitted 8 November, 2019; originally announced November 2019.

Comments: 25 pages, 7 figures, 7 tables

Conjugate Nearest Neighbor Gaussian Process Models for Efficient Statistical Interpolation of Large Spatial Data

Authors: Shinichiro Shirota, Andrew O. Finley, Bruce D. Cook, Sudipto Banerjee

Abstract: A key challenge in spatial statistics is the analysis for massive spatially-referenced data sets. Such analyses often proceed from Gaussian process specifications that can produce rich and robust inference, but involve dense covariance matrices that lack computationally exploitable structures. The matrix computations required for fitting such models involve floating point operations in cubic order… ▽ More A key challenge in spatial statistics is the analysis for massive spatially-referenced data sets. Such analyses often proceed from Gaussian process specifications that can produce rich and robust inference, but involve dense covariance matrices that lack computationally exploitable structures. The matrix computations required for fitting such models involve floating point operations in cubic order of the number of spatial locations and dynamic memory storage in quadratic order. Recent developments in spatial statistics offer a variety of massively scalable approaches. Bayesian inference and hierarchical models, in particular, have gained popularity due to their richness and flexibility in accommodating spatial processes. Our current contribution is to provide computationally efficient exact algorithms for spatial interpolation of massive data sets using scalable spatial processes. We combine low-rank Gaussian processes with efficient sparse approximations. Following recent work by [1], we model the low-rank process using a Gaussian predictive process (GPP) and the residual process as a sparsity-inducing nearest-neighbor Gaussian process (NNGP). A key contribution here is to implement these models using exact conjugate Bayesian modeling to avoid expensive iterative algorithms. Through the simulation studies, we evaluate performance of the proposed approach and the robustness of our models, especially for long range prediction. We implement our approaches for remotely sensed light detection and ranging (LiDAR) data collected over the US Forest Service Tanana Inventory Unit (TIU) in a remote portion of Interior Alaska. △ Less

Submitted 23 July, 2019; originally announced July 2019.

Characterizing functional relationships between anthropogenic and biological sounds: A western New York state soundscape case study

Authors: Jeffrey W. Doser, Kristina M. Hannam, Andrew O. Finley

Abstract: Roads are a widespread feature of landscapes worldwide, and road traffic sound potentially makes nearby habitat unsuitable for acoustically communicating organisms. It is important to understand the influence of roads at the soundscape level to mitigate negative impacts of road sound on individual species as well as subsequent effects on the surrounding landscape. We seek to characterize the relat… ▽ More Roads are a widespread feature of landscapes worldwide, and road traffic sound potentially makes nearby habitat unsuitable for acoustically communicating organisms. It is important to understand the influence of roads at the soundscape level to mitigate negative impacts of road sound on individual species as well as subsequent effects on the surrounding landscape. We seek to characterize the relationship between anthropogenic and biological sounds in western New York and assess the extent to which available traffic data explains variability in anthropogenic noise. Recordings were obtained in the spring of 2016 at 18 sites throughout western New York. We used the Welch Power Spectral Density (PSD) at low frequencies (0.5-2 kHz) to represent anthropogenic noise and PSD values at higher frequencies (2-11 kHz) to represent biological sound. Relationships were modeled using a novel two-stage hierarchical Bayesian model utilizing beta regression and basis splines. Model results and map predictions illustrate that anthropogenic noise and biological sound have an inverse relationship, and anthropogenic noise is greatest in close proximity to high traffic volume roads. The predictions have large uncertainty, resulting from the temporal coarseness of public road data used as a proxy for traffic sound. Results suggest that finer temporal resolution traffic sound data, such as crowd-sourced time-indexed traffic data from geographic positioning systems, might better account for observed temporal changes in the soundscape. The use of such data, in combination with the proposed modeling framework, could have important implications for the development of sound management policies. △ Less

Submitted 20 January, 2020; v1 submitted 3 May, 2019; originally announced May 2019.

Comments: 35 pages, 6 figures

Bayesian spatially varying coefficient models in the spBayes R package

Authors: Andrew O. Finley, Sudipto Banerjee

Abstract: This paper describes and illustrates new functionality for fitting spatially varying coefficients models in the spBayes (version 0.4-2) R package. The new spSVC function uses a computationally efficient Markov chain Monte Carlo algorithm and extends current spBayes functions, that fit only space-varying intercept regression models, to fit independent or multivariate Gaussian process random effects… ▽ More This paper describes and illustrates new functionality for fitting spatially varying coefficients models in the spBayes (version 0.4-2) R package. The new spSVC function uses a computationally efficient Markov chain Monte Carlo algorithm and extends current spBayes functions, that fit only space-varying intercept regression models, to fit independent or multivariate Gaussian process random effects for any set of columns in the regression design matrix. Newly added OpenMP parallelization options for spSVC are discussed and illustrated, as well as helper functions for joint and point-wise prediction and model fit diagnostics. The utility of the proposed models is illustrated using a PM10 analysis over central Europe. △ Less

Submitted 30 December, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

EcoMem: An R package for quantifying ecological memory

Authors: Malcolm S. Itter, Jarno Vanhatalo, Andrew O. Finley

Abstract: Ecological processes may exhibit memory to past disturbances affecting the resilience of ecosystems to future disturbance. Understanding the role of ecological memory in shaping ecosystem responses to disturbance under global change is a critical step toward developing effective adaptive management strategies to maintain ecosystem function and biodiversity. We developed EcoMem, an R package for qu… ▽ More Ecological processes may exhibit memory to past disturbances affecting the resilience of ecosystems to future disturbance. Understanding the role of ecological memory in shaping ecosystem responses to disturbance under global change is a critical step toward developing effective adaptive management strategies to maintain ecosystem function and biodiversity. We developed EcoMem, an R package for quantifying ecological memory functions using common environmental time series data (continuous, count, proportional) applying a Bayesian hierarchical framework. The package estimates memory functions for continuous and binary (e.g., disturbance chronology) variables making no a priori assumption on the form of the functions. EcoMem allows users to quantify ecological memory for a wide range of ecosystem processes and responses. The utility of the package to advance understanding of the memory of ecosystems to environmental drivers is demonstrated using a simulated dataset and a case study assessing the memory of boreal tree growth to insect defoliation. △ Less

Submitted 20 February, 2019; originally announced February 2019.

Comments: 11 pages, 3 figures, 2 appendices

Remote sensing to reduce the effects of spatial autocorrelation on design-based inference for forest inventory using systematic samples

Authors: Chad Babcock, Andrew O. Finley, Timothy G. Gregoire, Hans-Erik Andersen

Abstract: Systematic sampling is often used to select plot locations for forest inventory estimation. However, it is not possible to derive a design-unbiased variance estimator for a systematic sample using one random start. As a result, many forest inventory analysts resort to applying variance estimators that are design-unbiased following simple random sampling to their systematic samples even though this… ▽ More Systematic sampling is often used to select plot locations for forest inventory estimation. However, it is not possible to derive a design-unbiased variance estimator for a systematic sample using one random start. As a result, many forest inventory analysts resort to applying variance estimators that are design-unbiased following simple random sampling to their systematic samples even though this typically leads to conservative estimates of error. We explore the influence of spatial autocorrelation on variance estimation when systematic sampling is employed using repeated sampling. We generate a sequence of 1000 synthetic populations with increasing spatial autocorrelation and repeatedly sample from each to examine how the performance of estimators change as spatial autocorrelation changes. We also repeatedly sample from a tree census plot in Harvard Forest, Massachusetts and examine the performance of similar estimators. Results indicate that applying variance estimators that are unbiased following simple random sampling to systematic samples from populations exhibiting stronger spatial autocorrelation tend to be more conservative. We also find that incorporating ancillary wall-to-wall covariates, e.g., remote sensing data, using generalized regression estimators can reduce variance over-estimation by explaining some or all of the spatial structure in the population. △ Less

Submitted 19 October, 2018; originally announced October 2018.

Spatial Factor Models for High-Dimensional and Large Spatial Data: An Application in Forest Variable Mapping

Authors: Daniel Taylor-Rodriguez, Andrew O. Finley, Abhirup Datta, Chad Babcock, Hans-Erik Andersen, Bruce D. Cook, Douglas C. Morton, Sudipto Banerjee

Abstract: Gathering information about forest variables is an expensive and arduous activity. As such, directly collecting the data required to produce high-resolution maps over large spatial domains is infeasible. Next generation collection initiatives of remotely sensed Light Detection and Ranging (LiDAR) data are specifically aimed at producing complete-coverage maps over large spatial domains. Given that… ▽ More Gathering information about forest variables is an expensive and arduous activity. As such, directly collecting the data required to produce high-resolution maps over large spatial domains is infeasible. Next generation collection initiatives of remotely sensed Light Detection and Ranging (LiDAR) data are specifically aimed at producing complete-coverage maps over large spatial domains. Given that LiDAR data and forest characteristics are often strongly correlated, it is possible to make use of the former to model, predict, and map forest variables over regions of interest. This entails dealing with the high-dimensional ($\sim$$10^2$) spatially dependent LiDAR outcomes over a large number of locations (~10^5-10^6). With this in mind, we develop the Spatial Factor Nearest Neighbor Gaussian Process (SF-NNGP) model, and embed it in a two-stage approach that connects the spatial structure found in LiDAR signals with forest variables. We provide a simulation experiment that demonstrates inferential and predictive performance of the SF-NNGP, and use the two-stage modeling strategy to generate complete-coverage maps of forest variables with associated uncertainty over a large region of boreal forests in interior Alaska. △ Less

Submitted 8 November, 2018; v1 submitted 6 January, 2018; originally announced January 2018.

Geostatistical estimation of forest biomass in interior Alaska combining Landsat-derived tree cover, sampled airborne lidar and field observations

Authors: Chad Babcock, Andrew O. Finley, Hans-Erik Andersen, Robert Pattison, Bruce D. Cook, Douglas C. Morton, Michael Alonzo, Ross Nelson, Timothy Gregoire, Liviu Ene, Terje Gobakken, Erik Næsset

Abstract: The goal of this research was to develop and examine the performance of a geostatistical coregionalization modeling approach for combining field inventory measurements, strip samples of airborne lidar and Landsat-based remote sensing data products to predict aboveground biomass (AGB) in interior Alaska's Tanana Valley. The proposed modeling strategy facilitates pixel-level mapping of AGB density p… ▽ More The goal of this research was to develop and examine the performance of a geostatistical coregionalization modeling approach for combining field inventory measurements, strip samples of airborne lidar and Landsat-based remote sensing data products to predict aboveground biomass (AGB) in interior Alaska's Tanana Valley. The proposed modeling strategy facilitates pixel-level mapping of AGB density predictions across the entire spatial domain. Additionally, the coregionalization framework allows for statistically sound estimation of total AGB for arbitrary areal units within the study area---a key advance to support diverse management objectives in interior Alaska. This research focuses on appropriate characterization of prediction uncertainty in the form of posterior predictive coverage intervals and standard deviations. Using the framework detailed here, it is possible to quantify estimation uncertainty for any spatial extent, ranging from pixel-level predictions of AGB density to estimates of AGB stocks for the full domain. The lidar-informed coregionalization models consistently outperformed their counterpart lidar-free models in terms of point-level predictive performance and total AGB precision. Additionally, the inclusion of Landsat-derived forest cover as a covariate further improved estimation precision in regions with lower lidar sampling intensity. Our findings also demonstrate that model-based approaches that do not explicitly account for residual spatial dependence can grossly underestimate uncertainty, resulting in falsely precise estimates of AGB. On the other hand, in a geostatistical setting, residual spatial structure can be modeled within a Bayesian hierarchical framework to obtain statistically defensible assessments of uncertainty for AGB estimates. △ Less

Submitted 20 December, 2017; v1 submitted 9 May, 2017; originally announced May 2017.

Efficient algorithms for Bayesian Nearest Neighbor Gaussian Processes

Authors: Andrew O. Finley, Abhirup Datta, Bruce C. Cook, Douglas C. Morton, Hans E. Andersen, Sudipto Banerjee

Abstract: We consider alternate formulations of recently proposed hierarchical Nearest Neighbor Gaussian Process (NNGP) models (Datta et al., 2016a) for improved convergence, faster computing time, and more robust and reproducible Bayesian inference. Algorithms are defined that improve CPU memory management and exploit existing high-performance numerical linear algebra libraries. Computational and inferenti… ▽ More We consider alternate formulations of recently proposed hierarchical Nearest Neighbor Gaussian Process (NNGP) models (Datta et al., 2016a) for improved convergence, faster computing time, and more robust and reproducible Bayesian inference. Algorithms are defined that improve CPU memory management and exploit existing high-performance numerical linear algebra libraries. Computational and inferential benefits are assessed for alternate NNGP specifications using simulated datasets and remotely sensed light detection and ranging (LiDAR) data collected over the US Forest Service Tanana Inventory Unit (TIU) in a remote portion of Interior Alaska. The resulting data product is the first statistically robust map of forest canopy for the TIU. △ Less

Submitted 2 March, 2018; v1 submitted 1 February, 2017; originally announced February 2017.

A Model-Based Approach to Wildland Fire Reconstruction Using Sediment Charcoal Records

Authors: Malcolm S. Itter, Andrew O. Finley, Mevin B. Hooten, Philip E. Higuera, Jennifer R. Marlon, Ryan Kelly, Jason S. McLachlan

Abstract: Lake sediment charcoal records are used in paleoecological analyses to reconstruct fire history including the identification of past wildland fires. One challenge of applying sediment charcoal records to infer fire history is the separation of charcoal associated with local fire occurrence and charcoal originating from regional fire activity. Despite a variety of methods to identify local fires fr… ▽ More Lake sediment charcoal records are used in paleoecological analyses to reconstruct fire history including the identification of past wildland fires. One challenge of applying sediment charcoal records to infer fire history is the separation of charcoal associated with local fire occurrence and charcoal originating from regional fire activity. Despite a variety of methods to identify local fires from sediment charcoal records, an integrated statistical framework for fire reconstruction is lacking. We develop a Bayesian point process model to estimate probability of fire associated with charcoal counts from individual-lake sediments and estimate mean fire return intervals. A multivariate extension of the model combines records from multiple lakes to reduce uncertainty in local fire identification and estimate a regional mean fire return interval. The univariate and multivariate models are applied to 13 lakes in the Yukon Flats region of Alaska. Both models resulted in similar mean fire return intervals (100-350 years) with reduced uncertainty under the multivariate model due to improved estimation of regional charcoal deposition. The point process model offers an integrated statistical framework for paleo-fire reconstruction and extends existing methods to infer regional fire history from multiple lake records with uncertainty following directly from posterior distributions. △ Less

Submitted 7 December, 2016; originally announced December 2016.

Comments: 26 pages, 6 figures, 1 table

Joint hierarchical models for sparsely sampled high-dimensional LiDAR and forest variables

Authors: Andrew O. Finley, Sudipto Banerjee, Yuzhen Zhou, Bruce D. Cook, Chad Babcock

Abstract: Recent advancements in remote sensing technology, specifically Light Detection and Ranging (LiDAR) sensors, provide the data needed to quantify forest characteristics at a fine spatial resolution over large geographic domains. From an inferential standpoint, there is interest in prediction and interpolation of the often sparsely sampled and spatially misaligned LiDAR signals and forest variables.… ▽ More Recent advancements in remote sensing technology, specifically Light Detection and Ranging (LiDAR) sensors, provide the data needed to quantify forest characteristics at a fine spatial resolution over large geographic domains. From an inferential standpoint, there is interest in prediction and interpolation of the often sparsely sampled and spatially misaligned LiDAR signals and forest variables. We propose a fully process-based Bayesian hierarchical model for above ground biomass (AGB) and LiDAR signals. The process-based framework offers richness in inferential capabilities, e.g., inference on the entire underlying processes instead of estimates only at pre-specified points. Key challenges we obviate include misalignment between the AGB observations and LiDAR signals and the high-dimensionality in the model emerging from LiDAR signals in conjunction with the large number of spatial locations. We offer simulation experiments to evaluate our proposed models and also apply them to a challenging dataset comprising LiDAR and spatially coinciding forest inventory variables collected on the Penobscot Experimental Forest (PEF), Maine. Our key substantive contributions include AGB data products with associated measures of uncertainty for the PEF and, more broadly, a methodology that should find use in a variety of current and upcoming forest variable mapping efforts using sparsely sampled remotely sensed high-dimensional data. △ Less

Submitted 5 December, 2016; v1 submitted 23 March, 2016; originally announced March 2016.

Variable Effects of Climate on Forest Growth in Relation to Climate Extremes, Disturbance, and Forest Stand Dynamics

Authors: Malcolm S. Itter, Andrew O. Finley, Anthony W. D'Amato, Jane R. Foster, John B. Bradford

Abstract: Changes in the frequency, duration, and severity of climate extremes are forecast to occur under global climate change. The impacts of climate extremes on forest productivity and health are complicated by potential interactions with disturbance events and stand dynamics. The effects of stand dynamics on forest responses to climate and disturbance are particularly important given forest characteris… ▽ More Changes in the frequency, duration, and severity of climate extremes are forecast to occur under global climate change. The impacts of climate extremes on forest productivity and health are complicated by potential interactions with disturbance events and stand dynamics. The effects of stand dynamics on forest responses to climate and disturbance are particularly important given forest characteristics driven by stand dynamics can be modified through forest management with the goal of increasing forest resistance and resilience to climate change. We develop a hierarchical Bayesian state-space model allowing climate effects on tree growth to vary over time and in relation to climate extremes, disturbance events, and stand dynamics. We apply the model to a dendrochronology dataset comprising measurements from forest stands of varying composition, structure, and development stage in northeastern Minnesota. Results indicate average forest growth was most sensitive to variables describing climatic water deficit. Forest growth responses to water deficit were partitioned into responses driven by climatic threshold exceedances and interactions with forest tent caterpillar defoliation. Forest growth was both resistant and resilient to climate extremes with the majority of forest growth responses occurring after multiple climatic threshold exceedances or insect defoliation events. Forest growth was most sensitive to water deficit during periods of high stem density following major regeneration events when average inter-tree competition was high. Results suggest that forest growth resistance and resilience to interactions between climate extremes and insect defoliation can be increased through management steps such as thinning to reduce competition during early stages of stand development and small-group selection harvests to maintain forest structures characteristic of older, mature stands. △ Less

Submitted 2 December, 2016; v1 submitted 23 February, 2016; originally announced February 2016.

Comments: 30 pages, 3 tables, 7 figures

Non-separable Dynamic Nearest-Neighbor Gaussian Process Models for Large spatio-temporal Data With an Application to Particulate Matter Analysis

Authors: Abhirup Datta, Sudipto Banerjee, Andrew O. Finley, Nicholas A. S. Hamm, Martijn Schaap

Abstract: Particulate matter (PM) is a class of malicious environmental pollutants known to be detrimental to human health. Regulatory efforts aimed at curbing PM levels in different countries often require high resolution space-time maps that can identify red-flag regions exceeding statutory concentration limits. Continuous spatio-temporal Gaussian Process (GP) models can deliver maps depicting predicted P… ▽ More Particulate matter (PM) is a class of malicious environmental pollutants known to be detrimental to human health. Regulatory efforts aimed at curbing PM levels in different countries often require high resolution space-time maps that can identify red-flag regions exceeding statutory concentration limits. Continuous spatio-temporal Gaussian Process (GP) models can deliver maps depicting predicted PM levels and quantify predictive uncertainty. However, GP based approaches are usually thwarted by computational challenges posed by large datasets. We construct a novel class of scalable Dynamic Nearest Neighbor Gaussian Process (DNNGP) models that can provide a sparse approximation to any spatio-temporal GP (e.g., with non-separable covariance structures). The DNNGP we develop here can be used as a sparsity-inducing prior for spatio-temporal random effects in any Bayesian hierarchical model to deliver full posterior inference. Storage and memory requirements for a DNNGP model are linear in the size of the dataset thereby delivering massive scalability without sacrificing inferential richness. Extensive numerical studies reveal that the DNNGP provides substantially superior approximations to the underlying process than low rank approximations. Finally, we use the DNNGP to analyze a massive air quality dataset to substantially improve predictions of PM levels across Europe in conjunction with the LOTOS-EUROS chemistry transport models (CTMs). △ Less

Submitted 13 April, 2016; v1 submitted 24 October, 2015; originally announced October 2015.

Dynamic spatial regression models for space-varying forest stand tables

Authors: Andrew O. Finley, Sudipto Banerjee, Aaron R. Weiskittel, Chad Babcock, Bruce D. Cook

Abstract: Many forest management planning decisions are based on information about the number of trees by species and diameter per unit area. This information is commonly summarized in a stand table, where a stand is defined as a group of forest trees of sufficiently uniform species composition, age, condition, or productivity to be considered a homogeneous unit for planning purposes. Typically information… ▽ More Many forest management planning decisions are based on information about the number of trees by species and diameter per unit area. This information is commonly summarized in a stand table, where a stand is defined as a group of forest trees of sufficiently uniform species composition, age, condition, or productivity to be considered a homogeneous unit for planning purposes. Typically information used to construct stand tables is gleaned from observed subsets of the forest selected using a probability-based sampling design. Such sampling campaigns are expensive and hence only a small number of sample units are typically observed. This data paucity means that stand tables can only be estimated for relatively large areal units. Contemporary forest management planning and spatially explicit ecosystem models require stand table input at higher spatial resolution than can be affordably provided using traditional approaches. We propose a dynamic multivariate Poisson spatial regression model that accommodates both spatial correlation between observed diameter distributions and also correlation between tree counts across diameter classes within each location. To improve fit and prediction at unobserved locations, diameter specific intensities can be estimated using auxiliary data such as management history or remotely sensed information. The proposed model is used to analyze a diverse forest inventory dataset collected on the United States Forest Service Penobscot Experimental Forest in Bradley, Maine. Results demonstrate that explicitly modeling the residual spatial structure via a multivariate Gaussian process and incorporating information about forest structure from LiDAR covariates improve model fit and can provide high spatial resolution stand table maps with associated estimates of uncertainty. △ Less

Submitted 3 November, 2014; originally announced November 2014.

Hierarchical Nearest-Neighbor Gaussian Process Models for Large Geostatistical Datasets

Authors: Abhirup Datta, Sudipto Banerjee, Andrew O. Finley, Alan E. Gelfand

Abstract: Spatial process models for analyzing geostatistical data entail computations that become prohibitive as the number of spatial locations become large. This manuscript develops a class of highly scalable Nearest Neighbor Gaussian Process (NNGP) models to provide fully model-based inference for large geostatistical datasets. We establish that the NNGP is a well-defined spatial process providing legit… ▽ More Spatial process models for analyzing geostatistical data entail computations that become prohibitive as the number of spatial locations become large. This manuscript develops a class of highly scalable Nearest Neighbor Gaussian Process (NNGP) models to provide fully model-based inference for large geostatistical datasets. We establish that the NNGP is a well-defined spatial process providing legitimate finite-dimensional Gaussian densities with sparse precision matrices. We embed the NNGP as a sparsity-inducing prior within a rich hierarchical modeling framework and outline how computationally efficient Markov chain Monte Carlo (MCMC) algorithms can be executed without storing or decomposing large matrices. The floating point operations (flops) per iteration of this algorithm is linear in the number of spatial locations, thereby rendering substantial scalability. We illustrate the computational and inferential benefits of the NNGP over competing methods using simulation studies and also analyze forest biomass from a massive United States Forest Inventory dataset at a scale that precludes alternative dimension-reducing methods. △ Less

Submitted 1 January, 2016; v1 submitted 27 June, 2014; originally announced June 2014.

spBayes for large univariate and multivariate point-referenced spatio-temporal data models

Authors: Andrew O. Finley, Sudipto Banerjee, Alan E. Gelfand

Abstract: In this paper we detail the reformulation and rewrite of core functions in the spBayes R package. These efforts have focused on improving computational efficiency, flexibility, and usability for point-referenced data models. Attention is given to algorithm and computing developments that result in improved sampler convergence rate and efficiency by reducing parameter space; decreased sampler run-t… ▽ More In this paper we detail the reformulation and rewrite of core functions in the spBayes R package. These efforts have focused on improving computational efficiency, flexibility, and usability for point-referenced data models. Attention is given to algorithm and computing developments that result in improved sampler convergence rate and efficiency by reducing parameter space; decreased sampler run-time by avoiding expensive matrix computations, and; increased scalability to large datasets by implementing a class of predictive process models that attempt to overcome computational hurdles by representing spatial processes in terms of lower-dimensional realizations. Beyond these general computational improvements for existing model functions, we detail new functions for modeling data indexed in both space and time. These new functions implement a class of dynamic spatio-temporal models for settings where space is viewed as continuous and time is taken as discrete. △ Less

Submitted 30 October, 2013; originally announced October 2013.

Hierarchical spatial models for predicting tree species assemblages across large domains

Authors: Andrew O. Finley, Sudipto Banerjee, Ronald E. McRoberts

Abstract: Spatially explicit data layers of tree species assemblages, referred to as forest types or forest type groups, are a key component in large-scale assessments of forest sustainability, biodiversity, timber biomass, carbon sinks and forest health monitoring. This paper explores the utility of coupling georeferenced national forest inventory (NFI) data with readily available and spatially complete… ▽ More Spatially explicit data layers of tree species assemblages, referred to as forest types or forest type groups, are a key component in large-scale assessments of forest sustainability, biodiversity, timber biomass, carbon sinks and forest health monitoring. This paper explores the utility of coupling georeferenced national forest inventory (NFI) data with readily available and spatially complete environmental predictor variables through spatially-varying multinomial logistic regression models to predict forest type groups across large forested landscapes. These models exploit underlying spatial associations within the NFI plot array and the spatially-varying impact of predictor variables to improve the accuracy of forest type group predictions. The richness of these models incurs onerous computational burdens and we discuss dimension reducing spatial processes that retain the richness in modeling. We illustrate using NFI data from Michigan, USA, where we provide a comprehensive analysis of this large study area and demonstrate improved prediction with associated measures of uncertainty. △ Less

Submitted 8 October, 2009; originally announced October 2009.

Comments: Published in at http://dx.doi.org/10.1214/09-AOAS250 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AOAS-AOAS250

Journal ref: Annals of Applied Statistics 2009, Vol. 3, No. 3, 1052-1079