Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Environmetrics”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

17 recordsLinked to original sources
USGS1996

A model for field toxicity tests

Toxicity tests conducted under field conditions present an interesting challenge for statistical modelling. In contrast to laboratory tests, the concentrations of potential toxicants are not held constant over the test. In addition, the number and identity of toxicants that belong in a model as explanatory factors are not known and must be determined through a model selection process. We present one model to deal with these needs. This model takes the record of mortalities to form a multinomial distribution in which parameters are modelled as products of conditional daily survival probabilities. These conditional probabilities are in turn modelled as logistic functions of the explanatory factors. The model incorporates lagged values of the explanatory factors to deal with changes in the pattern of mortalities over time. The issue of model selection and assessment is approached through the use of generalized information criteria and power divergence goodness-of-fit tests. These model selection criteria are applied in a cross-validation scheme designed to assess the ability of a model to both fit data used in estimation and predict data deleted from the estimation data set. The example presented demonstrates the need for inclusion of lagged values of the explanatory factors and suggests that penalized likelihood criteria may not provide adequate protection against overparameterized models in model selection.

Environmetrics↗
USGS2011

Semiparametric bivariate zero-inflated Poisson models with application to studies of abundance for multiple species

Ecological studies involving counts of abundance, presence–absence or occupancy rates often produce data having a substantial proportion of zeros. Furthermore, these types of processes are typically multivariate and only adequately described by complex nonlinear relationships involving externally measured covariates. Ignoring these aspects of the data and implementing standard approaches can lead to models that fail to provide adequate scientific understanding of the underlying ecological processes, possibly resulting in a loss of inferential power. One method of dealing with data having excess zeros is to consider the class of univariate zero-inflated generalized linear models. However, this class of models fails to address the multivariate and nonlinear aspects associated with the data usually encountered in practice. Therefore, we propose a semiparametric bivariate zero-inflated Poisson model that takes into account both of these data attributes. The general modeling framework is hierarchical Bayes and is suitable for a broad range of applications. We demonstrate the effectiveness of our model through a motivating example on modeling catch per unit area for multiple species using data from the Missouri River Benthic Fishes Study, implemented by the United States Geological Survey.

Environmetrics↗
USGS2013

A comparison of adaptive sampling designs and binary spatial models: A simulation study using a census of Bromus inermis

Commonly in environmental and ecological studies, species distribution data are recorded as presence or absence throughout a spatial domain of interest. Field based studies typically collect observations by sampling a subset of the spatial domain. We consider the effects of six different adaptive and two non-adaptive sampling designs and choice of three binary models on both predictions to unsampled locations and parameter estimation of the regression coefficients (species–environment relationships). Our simulation study is unique compared to others to date in that we virtually sample a true known spatial distribution of a nonindigenous plant species, Bromus inermis . The census of B. inermis provides a good example of a species distribution that is both sparsely (1.9 % prevalence) and patchily distributed. We find that modeling the spatial correlation using a random effect with an intrinsic Gaussian conditionally autoregressive prior distribution was equivalent or superior to Bayesian autologistic regression in terms of predicting to un-sampled areas when strip adaptive cluster sampling was used to survey B. inermis . However, inferences about the relationships between B. inermis presence and environmental predictors differed between the two spatial binary models. The strip adaptive cluster designs we investigate provided a significant advantage in terms of Markov chain Monte Carlo chain convergence when trying to model a sparsely distributed species across a large area. In general, there was little difference in the choice of neighborhood, although the adaptive king was preferred when transects were randomly placed throughout the spatial domain.

Environmetrics↗
USGS2014

County-level analysis of the impact of temperature and population increases on California wildfire data

The extent to which the apparent increase in wildfire incidence and burn area in California from 1990 to 2006 is affected by population and temperature increases is examined. Using generalized linear models with random effects, we focus on the estimated impacts of increases in mean daily temperatures and populations in different counties on wildfire in those counties, after essentially controlling for the overall differences between counties in their overall mean temperatures and populations. We find that temperature increase appears to have a significant positive impact on both total burn area and number of observed wildfires. Population growth appears to have a much less pronounced impact on total burn area than do annual temperature increases, and population growth appears to be negatively correlated with the total number of observed wildfires. These effects are especially pronounced in the winter season and in Southern California counties.

California↗
USGS2015

Reconstruction of late Holocene climate based on tree growth and mechanistic hierarchical models

Reconstruction of pre-instrumental, late Holocene climate is important for understanding how climate has changed in the past and how climate might change in the future. Statistical prediction of paleoclimate from tree ring widths is challenging because tree ring widths are a one-dimensional summary of annual growth that represents a multi-dimensional set of climatic and biotic influences. We develop a Bayesian hierarchical framework using a nonlinear, biologically motivated tree ring growth model to jointly reconstruct temperature and precipitation in the Hudson Valley, New York. Using a common growth function to describe the response of a tree to climate, we allow for species-specific parameterizations of the growth response. To enable predictive backcasts, we model the climate variables with a vector autoregressive process on an annual timescale coupled with a multivariate conditional autoregressive process that accounts for temporal correlation and cross-correlation between temperature and precipitation on a monthly scale. Our multi-scale temporal model allows for flexibility in the climate response through time at different temporal scales and predicts reasonable climate scenarios given tree ring width data.

Environmetrics↗
USGS2016

Hierarchical animal movement models for population-level inference

New methods for modeling animal movement based on telemetry data are developed regularly. With advances in telemetry capabilities, animal movement models are becoming increasingly sophisticated. Despite a need for population-level inference, animal movement models are still predominantly developed for individual-level inference. Most efforts to upscale the inference to the population level are either post hoc or complicated enough that only the developer can implement the model. Hierarchical Bayesian models provide an ideal platform for the development of population-level animal movement models but can be challenging to fit due to computational limitations or extensive tuning required. We propose a two-stage procedure for fitting hierarchical animal movement models to telemetry data. The two-stage approach is statistically rigorous and allows one to fit individual-level movement models separately, then resample them using a secondary MCMC algorithm. The primary advantages of the two-stage approach are that the first stage is easily parallelizable and the second stage is completely unsupervised, allowing for an automated fitting procedure in many cases. We demonstrate the two-stage procedure with two applications of animal movement models. The first application involves a spatial point process approach to modeling telemetry data, and the second involves a more complicated continuous-time discrete-space animal movement model. We fit these models to simulated data and real telemetry data arising from a population of monitored Canada lynx in Colorado, USA.

Environmetrics↗
USGS2017

A model-based approach to wildland fire reconstruction using sediment charcoal records

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 the 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 paleofire reconstruction and extends existing methods to infer regional fire history from multiple lake records with uncertainty following directly from posterior distributions.

Environmetrics↗
USGS2017

Group inverse sampling: An economical approach to inverse sampling

Inverse sampling is an adaptive design in the sense that the final sampling effort during a search for rare events will depend on what is found during the survey. Conventional inverse sampling (CIS) designs successively select individual sampling units to find, for example, the k th rare event. In real sampling situations, use of successive one‐by‐one sampling can be cost prohibitive. Here, we introduce an inverse sampling design that uses successive selection of groups instead of individuals, named group inverse sampling (GIS). An unbiased estimator and its variance estimator of the population mean are derived based on the Murthy estimator. CIS is a special case of the generalized design with group size equal to one. We simulate the GIS design to evaluate its efficiency using populations of rare freshwater mussels in West Virginia, USA. For cost consideration, we calculate distance traveled among the sampling units. Results show that GIS was more cost efficient than CIS in all cases. The group size for successive sampling (d ) was the most influential design parameter for reducing cost and increasing precision. Also, GIS found more rare units with greater consistency compared to simple random sampling without replacement (SRS). An important characteristic of the GIS design is that sampling stops when the target number of rare units is found, which prevents unnecessary sampling and contrasts favorably with other adaptive designs such as adaptive cluster sampling.

Environmetrics↗
USGS2018

Process convolution approaches for modeling interacting trajectories

Gaussian processes are a fundamental statistical tool used in a wide range of applications. In the spatiotemporal setting, several families of covariance functions exist to accommodate a wide variety of dependence structures arising in different applications. These parametric families can be restrictive and are insufficient in some situations. In contrast, process convolutions represent a flexible, interpretable approach to defining the covariance of a Gaussian process and have modest requirements to ensure validity. We introduce a generalization of the process convolution approach that employs multiple convolutions sequentially to form a “process convolution chain”. In our proposed multistage framework, complex dependencies that arise from a combination of different interacting mechanisms are decomposed into a series of interpretable kernel smoothers. We demonstrate an application of process convolution chains to model killer whale movement, in which the paths taken by multiple individuals are not independent but reflect dynamic social interactions within the population. Our proposed model for dependent movement provides inference for the latent dynamic social structure in the study population. Additionally, by leveraging the positive dependence among individual paths, we achieve a reduction in uncertainty for the estimated locations of the killer whales, compared with a model that treats paths as independent.

Environmetrics↗
USGS2019

Estimating population size with imperfect detection using a parametric bootstrap

We develop a novel method of estimating population size from imperfectly detected counts of individuals and a separate estimate of detection probability. Observed counts are separated into classes within which detection probability is assumed constant. Within a detection class, counts are modeled as a single binomial observation X with success probability p where the goal is to estimate index N . We use a Horvitz–Thompson‐like estimator for N and account for uncertainty in both sample data and estimated success probability via a parametric bootstrap. Unlike capture–recapture methods, our model does not require repeated sampling of the population. Our method is able to achieve good results, even with small X . We show in a factorial simulation study that the median of the bootstrapped sample has small bias relative to N and that coverage probabilities of confidence intervals for N are near nominal under a wide array of scenarios. Our methodology begins to break down when P ( X =0)>0.1 but is still capable of obtaining reasonable confidence coverage. We illustrate the proposed technique by estimating (1) the size of a moose population in Alaska and (2) the number of bat fatalities at a wind power facility, both from samples with imperfect detection probabilities, estimated independently.

Environmetrics↗
USGS2019

Nonlinear reaction–diffusion process models improve inference for population dynamics

Partial differential equations (PDEs) are a useful tool for modeling spatiotemporal dynamics of ecological processes. However, as an ecological process evolves, we need statistical models that can adapt to changing dynamics as new data are collected. We developed a model that combines an ecological diffusion equation and logistic growth to characterize colonization processes of a population that establishes long-term equilibrium over a heterogeneous environment. We also developed a homogenization strategy to statistically upscale the PDE for faster computation and adopted a hierarchical framework to accommodate multiple data sources collected at different spatial scales. We highlighted the advantages of using a logistic reaction component instead of a Malthusian component when population growth demonstrates asymptotic behavior. As a case study, we demonstrated that our model improves spatiotemporal abundance forecasts of sea otters in Glacier Bay, Alaska. Furthermore, we predicted spatially varying local equilibrium abundances as a result of environmentally driven diffusion and density-regulated growth. Integrating equilibrium abundances over the study area in our application enabled us to infer the overall carrying capacity of sea otters in Glacier Bay, Alaska.

Alaska↗
USGS2023

Bayesian spatio-temporal survival analysis for all types of censoring with application to a wildlife disease study

In this article, we consider modeling arbitrarily censored survival data with spatio-temporal covariates. We demonstrate that under the piecewise constant hazard function, the likelihood for uncensored or right-censored subjects is proportional to the likelihood of multiple conditionally independent Poisson random variables. To address left- or interval-censored subjects, we propose to impute the exact event times and convert them into uncensored subjects, enabling the application of the integrated nested Laplace approximation to update model parameters using the imputed data. We introduce an iterative algorithm that alternates between imputing event times for left- and interval-censored subjects and re-estimating model parameters. The proposed method is assessed through a simulation study and applied to analyze a spatio-temporal survival dataset in a wildlife disease study investigating bovine tuberculosis in white-tailed deer in Michigan.

Environmetrics↗
USGS2025

Spatiotemporal causal inference with mechanistic ecological models: Evaluating targeted culling on chronic wasting disease dynamics in cervids

Spatiotemporal causal inference methods are needed to detect the effect of interventions on indirectly measured epidemiological outcomes that go beyond studying spatiotemporal correlations. Chronic wasting disease (CWD) causes neurological degeneration and eventual death to white-tailed deer ( Odocoileus virginianus ) in Wisconsin. Targeted culling involves removing deer after traditional hunting seasons in areas with high CWD prevalence. The evaluation of the causal effects of targeted culling in the spread and growth of CWD is an important unresolved research and CWD management question that can guide surveillance efforts. Reaction–diffusion partial differential equations (PDEs) can be used to mechanistically model the underlying spatiotemporal dynamics of wildlife diseases, like CWD, allowing researchers to make inference about unobserved epidemiological quantities. These models indirectly regress spatiotemporal covariates on diffusion and growth rates parameterizing such PDEs, obtaining associational conclusions. In this work we develop an innovative method to obtain causal estimators for the effect of targeted culling interventions on CWD epidemiological processes using an inverse-probability-of-treatment-weighted technique by means of marginal structural models embedded in the PDE fitting process. Additionally we establish a novel scheme for sensitivity analysis under unmeasured confounder for testing the hypothesis of a significant causal effect in the indirectly measured epidemiological outcomes. Our methods can be broadly used to study the impact of spatiotemporal interventions and treatment exposures in the epidemiological evolution of infectious diseases that can help to inform future efforts to mitigate public health implications and wildlife disease burden.

Wisconsin↗
USGS2002

Large scale wildlife monitoring studies: Statistical methods for design and analysis

Techniques for estimation of absolute abundance of wildlife populations have received a lot of attention in recent years. The statistical research has been focused on intensive small-scale studies. Recently, however, wildlife biologists have desired to study populations of animals at very large scales for monitoring purposes. Population indices are widely used in these extensive monitoring programs because they are inexpensive compared to estimates of absolute abundance. A crucial underlying assumption is that the population index ( C ) is directly proportional to the population density ( D ). The proportionality constant, β , is simply the probability of 'detection' for animals in the survey. As spatial and temporal comparisons of indices are crucial, it is necessary to also assume that the probability of detection is constant over space and time. Biologists intuitively recognize this when they design rigid protocols for the studies where the indices are collected. Unfortunately, however in many field studios the assumption is clearly invalid. We believe that the estimation of detection probability should be built into the monitoring design through a double sampling approach. A large sample of points provides an abundance index, and a smaller sub-sample of the same points is used to estimate detection probability. There is an important need for statistical research on the design and analysis of these complex studies. Some basic concepts based on actual avian, amphibian, and fish monitoring studies are presented in this article.

Environmetrics↗
USGS2002

Defining surfaces for skewed, highly variable data

Skewness of environmental data is often caused by more than simply a handful of outliers in an otherwise normal distribution. Statistical procedures for such datasets must be sufficiently robust to deal with distributions that are strongly non-normal, containing both a large proportion of outliers and a skewed main body of data. In the field of water quality, skewness is commonly associated with large variation over short distances. Spatial analysis of such data generally requires either considerable effort at modeling or the use of robust procedures not strongly affected by skewness and local variability. Using a skewed dataset of 675 nitrate measurements in ground water, commonly used methods for defining a surface (least-squares regression and kriging) are compared to a more robust method (loess). Three choices are critical in defining a surface: (i) is the surface to be a central mean or median surface? (ii) is either a well-fitting transformation or a robust and scale-independent measure of center used? (iii) does local spatial autocorrelation assist in or detract from addressing objectives? Published in 2002 by John Wiley & Sons, Ltd.

Environmetrics↗
USGS2025

Animal trajectory imputation and uncertainty quantification via deep learning

Imputing missing data in animal trajectories is crucial for understanding animal movements during unobserved periods. However, the traditional methods, such as linear interpolation and the continuous-time correlated random walk model, are often inadequate to capture the complexity of animal movements. Here, we develop a deep learning approach to animal trajectory imputation by a conditional diffusion model. Unlike the traditional methods, our deep learning method uses observed data and external covariates to impute missing positions along an animal trajectory, capturing periodic patterns and the influence of covariates, which leads to more accurate imputations. In a case study of imputing deer trajectories, our method not only provides more accurate deterministic imputations than existing approaches but also achieves uncertainty quantification through probabilistic imputation.

Environmetrics↗
USGS2005

Dynamic design of ecological monitoring networks for non-Gaussian spatio-temporal data

Many ecological processes exhibit spatial structure that changes over time in a coherent, dynamical fashion. This dynamical component is often ignored in the design of spatial monitoring networks. Furthermore, ecological variables related to processes such as habitat are often non-Gaussian (e.g. Poisson or log-normal). We demonstrate that a simulation-based design approach can be used in settings where the data distribution is from a spatio-temporal exponential family. The key random component in the conditional mean function from this distribution is then a spatio-temporal dynamic process. Given the computational burden of estimating the expected utility of various designs in this setting, we utilize an extended Kalman filter approximation to facilitate implementation. The approach is motivated by, and demonstrated on, the problem of selecting sampling locations to estimate July brood counts in the prairie pothole region of the U.S.

Environmetrics↗