Search USGS⌕ Search

SEARCH · Search USGS

Results for “Statistical Methodology”

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.

At least 19 recordsLinked to original sources
USGS2010

A statistical methodology for estimating transport parameters: Theory and applications to one-dimensional advectivec-dispersive systems

A simulation nonlinear multiple-regression methodology for estimating parameters that characterize the transport of contaminants is developed and demonstrated. Finite difference contaminant transport simulation is combined with a nonlinear weighted least squares multiple-regression procedure. The technique provides optimal parameter estimates and gives statistics for assessing the reliability of these estimates under certain general assumptions about the distributions of the random measurement errors. Monte Carlo analysis is used to estimate parameter reliability for a hypothetical homogeneous soil column for which concentration data contain large random measurement errors. The value of data collected spatially versus data collected temporally was investigated for estimation of velocity, dispersion coefficient, effective porosity, first-order decay rate, and zero-order production. The use of spatial data gave estimates that were 2–3 times more reliable than estimates based on temporal data for all parameters except velocity. Comparison of estimated linear and nonlinear confidence intervals based upon Monte Carlo analysis showed that the linear approximation is poor for dispersion coefficient and zero-order production coefficient when data are collected over time. In addition, examples demonstrate transport parameter estimation for two real one-dimensional systems. First, the longitudinal dispersivity and effective porosity of an unsaturated soil are estimated using laboratory column data. We compare the reliability of estimates based upon data from individual laboratory experiments versus estimates based upon pooled data from several experiments. Second, the simulation nonlinear regression procedure is extended to include an additional governing equation that describes delayed storage during contaminant transport. The model is applied to analyze the trends, variability, and interrelationship of parameters in a mourtain stream in northern California.

Water Resources Research↗
USGS2003

Thematic accuracy of the 1992 National Land-Cover Data for the eastern United States: Statistical methodology and regional results

The accuracy of the 1992 National Land-Cover Data (NLCD) map is assessed via a probability sampling design incorporating three levels of stratification and two stages of selection. Agreement between the map and reference land-cover labels is defined as a match between the primary or alternate reference label determined for a sample pixel and a mode class of the mapped 3×3 block of pixels centered on the sample pixel. Results are reported for each of the four regions comprising the eastern United States for both Anderson Level I and II classifications. Overall accuracies for Levels I and II are 80% and 46% for New England, 82% and 62% for New York/New Jersey (NY/NJ), 70% and 43% for the Mid-Atlantic, and 83% and 66% for the Southeast.

Remote Sensing of Environment↗
USGS2024

Spatio-temporal ecological models via physics-informed neural networks for studying chronic wasting disease

To mitigate the negative effects of emerging wildlife diseases in biodiversity and public health it is critical to accurately forecast pathogen dissemination while incorporating relevant spatio-temporal covariates. Forecasting spatio-temporal processes can often be improved by incorporating scientific knowledge about the dynamics of the process using physical models. Ecological diffusion equations are often used to model epidemiological processes of wildlife diseases where environmental factors play a role in disease spread. Physics-informed neural networks (PINN) are deep learning algorithms that constrain neural network predictions based on physical laws and therefore are powerful forecasting models useful even in cases of limited and imperfect training data. In this paper, we develop a novel ecological modeling tool using PINNs, which fits a feedforward neural network and simultaneously performs parameter identification in a partial differential equation (PDE) with varying coefficients. We demonstrate the applicability of our model by comparing it with the commonly used Bayesian stochastic partial differential equation method and traditional machine learning approaches, showing that our proposed model exhibits superior prediction and forecasting performance when modeling chronic wasting disease in deer in Wisconsin. Furthermore, our model provides the opportunity to obtain scientific insights into spatiotemporal covariates affecting spread and growth of diseases. This work contributes to future machine learning and statistical methodology development by studying spatio-temporal processes enhanced by prior physical knowledge.

Spatial Statistics↗
USGS2017

Basis function models for animal movement

Advances in satellite-based data collection techniques have served as a catalyst for new statistical methodology to analyze these data. In wildlife ecological studies, satellite-based data and methodology have provided a wealth of information about animal space use and the investigation of individual-based animal–environment relationships. With the technology for data collection improving dramatically over time, we are left with massive archives of historical animal telemetry data of varying quality. While many contemporary statistical approaches for inferring movement behavior are specified in discrete time, we develop a flexible continuous-time stochastic integral equation framework that is amenable to reduced-rank second-order covariance parameterizations. We demonstrate how the associated first-order basis functions can be constructed to mimic behavioral characteristics in realistic trajectory processes using telemetry data from mule deer and mountain lion individuals in western North America. Our approach is parallelizable and provides inference for heterogenous trajectories using nonstationary spatial modeling techniques that are feasible for large telemetry datasets. Supplementary materials for this article are available online.

Journal of the American Statistical Association↗
USGS2010

A general science-based framework for dynamical spatio-temporal models

Spatio-temporal statistical models are increasingly being used across a wide variety of scientific disciplines to describe and predict spatially-explicit processes that evolve over time. Correspondingly, in recent years there has been a significant amount of research on new statistical methodology for such models. Although descriptive models that approach the problem from the second-order (covariance) perspective are important, and innovative work is being done in this regard, many real-world processes are dynamic, and it can be more efficient in some cases to characterize the associated spatio-temporal dependence by the use of dynamical models. The chief challenge with the specification of such dynamical models has been related to the curse of dimensionality. Even in fairly simple linear, first-order Markovian, Gaussian error settings, statistical models are often over parameterized. Hierarchical models have proven invaluable in their ability to deal to some extent with this issue by allowing dependency among groups of parameters. In addition, this framework has allowed for the specification of science based parameterizations (and associated prior distributions) in which classes of deterministic dynamical models (e. g., partial differential equations (PDEs), integro-difference equations (IDEs), matrix models, and agent-based models) are used to guide specific parameterizations. Most of the focus for the application of such models in statistics has been in the linear case. The problems mentioned above with linear dynamic models are compounded in the case of nonlinear models. In this sense, the need for coherent and sensible model parameterizations is not only helpful, it is essential. Here, we present an overview of a framework for incorporating scientific information to motivate dynamical spatio-temporal models. First, we illustrate the methodology with the linear case. We then develop a general nonlinear spatio-temporal framework that we call general quadratic nonlinearity and demonstrate that it accommodates many different classes of scientific-based parameterizations as special cases. The model is presented in a hierarchical Bayesian framework and is illustrated with examples from ecology and oceanography. ?? 2010 Sociedad de Estad??stica e Investigaci??n Operativa.

Test↗
USGS2000

Microsatellites: Evolutionary and methodological background and empirical applications at individual, population, and phylogenetic levels

The recent proliferation and greater accessibility of molecular genetic markers has led to a growing appreciation of the ecological and evolutionary inferences that can be drawn from molecular characterizations of individuals and populations (Burke et al. 1992, Avise 1994). Different techniques have the ability to target DNA sequences which have different patterns of inheritance, different modes and rates of evolution and, concomitantly, different levels of variation. In the quest for 'the right marker for the right job', microsatellites have been widely embraced as the marker of choice for many empirical genetic studies. The proliferation of microsatellite loci for various species and the voluminous literature compiled in very few years associated with their evolution and use in various research applications, exemplifies their growing importance as a research tool in the biological sciences. The ability to define allelic states based on variation at the nucleotide level has afforded unparalleled opportunities to document the actual mutational process and rates of evolution at individual microsatellite loci. The scrutiny to which these loci have been subjected has resulted in data that raise issues pertaining to assumptions formerly stated, but largely untestable for other marker classes. Indeed this is an active arena for theoretical and empirical work. Given the extensive and ever-increasing literature on various statistical methodologies and cautionary notes regarding the uses of microsatellites, some consideration should be given to the unique characteristics of these loci when determining how and under what conditions they can be employed.

Book chapter↗
USGS2021

Caution is warranted when using animal space-use and movement to infer behavioral states

Background Identifying the behavioral state for wild animals that can’t be directly observed is of growing interest to the ecological community. Advances in telemetry technology and statistical methodologies allow researchers to use space-use and movement metrics to infer the underlying, latent, behavioral state of an animal without direct observations. For example, researchers studying ungulate ecology have started using these methods to quantify behaviors related to mating strategies. However, little work has been done to determine if assumed behaviors inferred from movement and space-use patterns correspond to actual behaviors of individuals. Methods Using a dataset with male and female white-tailed deer location data, we evaluated the ability of these two methods to correctly identify male-female interaction events (MFIEs). We identified MFIEs using the proximity of their locations in space as indicators of when mating could have occurred. We then tested the ability of utilization distributions (UDs) and hidden Markov models (HMMs) rendered with single sex location data to identify these events. Results For white-tailed deer, male and female space-use and movement behavior did not vary consistently when with a potential mate. There was no evidence that a probability contour threshold based on UD volume applied to an individual’s UD could be used to identify MFIEs. Additionally, HMMs were unable to identify MFIEs, as single MFIEs were often split across multiple states and the primary state of each MFIE was not consistent across events. Conclusions Caution is warranted when interpreting behavioral insights rendered from statistical models applied to location data, particularly when there is no form of validation data. For these models to detect latent behaviors, the individual needs to exhibit a consistently different type of space-use and movement when engaged in the behavior. Unvalidated assumptions about that relationship may lead to incorrect inference about mating strategies or other behaviors.

Movement Ecology↗
USGS2018

Effects of environmental variables on invasive amphibian activity: Using model selection on quantiles for counts

Many different factors influence animal activity. Often, the value of an environmental variable may influence significantly the upper or lower tails of the activity distribution. For describing relationships with heterogeneous boundaries, quantile regressions predict a quantile of the conditional distribution of the dependent variable. A quantile count model extends linear quantile regression methods to discrete response variables, and is useful if activity is quantified by trapping, where there may be many tied (equal) values in the activity distribution, over a small range of discrete values. Additionally, different environmental variables in combination may have synergistic or antagonistic effects on activity, so examining their effects together, in a modeling framework, is a useful approach. Thus, model selection on quantile counts can be used to determine the relative importance of different variables in determining activity, across the entire distribution of capture results. We conducted model selection on quantile count models to describe the factors affecting activity (numbers of captures) of cane toads ( Rhinella marina ) in response to several environmental variables (humidity, temperature, rainfall, wind speed, and moon luminosity) over eleven months of trapping. Environmental effects on activity are understudied in this pest animal. In the dry season, model selection on quantile count models suggested that rainfall positively affected activity, especially near the lower tails of the activity distribution. In the wet season, wind speed limited activity near the maximum of the distribution, while minimum activity increased with minimum temperature. This statistical methodology allowed us to explore, in depth, how environmental factors influenced activity across the entire distribution, and is applicable to any survey or trapping regime, in which environmental variables affect activity.

Ecosphere↗
USGS2024

Assessing giant sequoia mortality and regeneration following high-severity wildfire

Fire is a critical driver of giant sequoia ( Sequoiadendron giganteum [Lindl.] Buchholz) regeneration. However, fire suppression combined with the effects of increased temperature and severe drought has resulted in fires of an intensity and size outside of the historical norm. As a result, recent mega-fires have killed a significant portion of the world's sequoia population (13%–19%), and uncertainty surrounds whether severely affected groves will be able to recover naturally, potentially leading to a loss of grove area. To assess the likelihood of natural recovery, we collected spatially explicit data assessing mortality, crown condition, and regeneration within four giant sequoia groves that were severely impacted by the SQF- (2020) and KNP-Complex (2021) wildfires within Sequoia and Kings Canyon National Parks. In total, we surveyed 5.9 ha for seedlings and assessed the crown condition of 1104 giant sequoias. To inform management, we used a statistical methodology that robustly quantifies the uncertainty in inherently “noisy” seedling data and takes advantage of readily available remote sensing metrics that would make our findings applicable to other recently burned groves. A loss of giant sequoia grove area would be a consequence of giant sequoia tree mortality followed by a failure of natural regeneration. We found that areas that experienced very high-severity fire (above ~800 RdNBR) are at substantial risk for the loss of grove area, with tree mortality rapidly increasing and giant sequoia seedling density simultaneously decreasing with fire severity. Such high-severity areas comprised 17.8, 142.0, 14.6, 1.6 ha and ~90%, ~14%, ~53%, and ~27% of Board Camp, Redwood Mountain, Suwanee, and New Oriole Lake groves, respectively. In all sampling areas, we found that seedling densities fell far below the average density measured after prescribed fires, where seedling numbers were almost certainly adequate to maintain giant sequoia populations and postfire conditions were more in keeping with historical norms. Importantly, spatial pattern is also important in assessing the risk of grove loss, and in two groves, Suwanee and New Oriole Lake, the high-severity patches were not always contiguous, potentially making some areas more resilient to regeneration failure due to the proximity of surviving trees.

Ecosphere↗
USGS2024

Quantifying the relative importance of biotic and abiotic factors in landscape-based models of stream fish distributions

Lotic fish species distributions are frequently predicted using remotely sensed habitat variables that characterize the adjacent landscape and serve as proxies for instream habitat. Recent advancements in statistical methodology, however, allow for leveraging fish assemblage data when predicting distributions. This is important because assemblage composition likely provides better information about instream habitat compared to landscape-derived metrics and therefore may improve predictions. To better understand the value of using multi-species fish data in species distribution modeling, we fit two conditional random fields (CRF) models to quantify the relative importance of fish assemblage co-occurrence, landscape-derived habitat variables, and interactions between these two predictor groups (i.e., effects of co-occurrence could be context-dependent) at over 1200 stream catchments in Pennsylvania, USA. We first compared predictive performance of CRF models against traditionally used single-species logistic regressions (generalized linear models; GLMs) and found that inclusion of fish assemblage data often improved predictive performance. The multi-species CRF models performed significantly better at predicting occurrence for 63% of species with an average percent increase in AUC of 25% compared to GLMs. Furthermore, the CRF identified species co-occurrences as more informative, and thus relatively more important, at predicting occurrence than the other effect types. The CRF also suggested that allowing these biotic effects to be context-dependent was important for predicting occurrence of many species. These findings illustrate the value of fish assemblage data for landscape-scale species distribution modeling and leveraging this information can improve predictions and inferences to help inform the management and conservation of freshwater fishes.

Pennsylvania↗
USGS2024

Combining terrestrial lidar with single line transects to investigate geomorphic change: A case study on the Upper Verde River, Arizona

The Upper Verde River in northern Arizona, USA is a vital resource for the wildlife and humans that rely on its waters. We characterize the riparian corridor topography using terrestrial laser scanner (TLS) data from 2021 to 2022. We also quantify geomorphic changes associated with human and climate-driven alterations in river flow and vegetation changes by combining the contemporary lidar surveys with legacy measurements from single line geomorphology transects measured by the United States Forest Service (USFS) in 2009. Seventeen plots along the Upper Verde River were surveyed with the TLS and the data were coregistered within individual plots with a Root Mean Square Error of <0.03 m among scan positions. Digital Elevation Models (DEM) were derived for each plot from the TLS data at 10 cm resolution and compared to the 2009 USFS cross-section data to quantify elevation changes. In areas with statistically significant change, we detected maximum changes in elevation due to erosion and deposition of −0.37 m and + 0.97 m, respectively. Topographic changes over the 13-year period were predominately aggradation and associated with sediment deposition, which we hypothesize might have resulted from altered river flow and vegetation encroachment. This study also demonstrates a quantitative and statistical methodology to fuse traditional single line cross-section data with contemporary lidar data to quantify geomorphic change. The novel approach demonstrated here is broadly applicable to natural resource managers for integrating and contextualizing legacy topographic data for understanding past, present, and future landscape and habitat changes.

Arizona↗
USGS2004

On the estimation of dispersal and movement of birds

The estimation of dispersal and movement is important to evolutionary and population ecologists, as well as to wildlife managers. We review statistical methodology available to estimate movement probabilities. We begin with cases where individual birds can be marked and their movements estimated with the use of multisite capture-recapture methods. Movements can be monitored either directly, using telemetry, or by accounting for detection probability when conventional marks are used. When one or more sites are unobservable, telemetry, band recoveries, incidental observations, a closed- or open-population robust design, or partial determinism in movements can be used to estimate movement. When individuals cannot be marked, presence-absence data can be used to model changes in occupancy over time, providing indirect inferences about movement. Where abundance estimates over time are available for multiple sites, potential coupling of their dynamics can be investigated using linear cross-correlation or nonlinear dynamic tools.

Condor↗
USGS2015

A functional model for characterizing long-distance movement behaviour

Advancements in wildlife telemetry techniques have made it possible to collect large data sets of highly accurate animal locations at a fine temporal resolution. These data sets have prompted the development of a number of statistical methodologies for modelling animal movement. Telemetry data sets are often collected for purposes other than fine-scale movement analysis. These data sets may differ substantially from those that are collected with technologies suitable for fine-scale movement modelling and may consist of locations that are irregular in time, are temporally coarse or have large measurement error. These data sets are time-consuming and costly to collect but may still provide valuable information about movement behaviour. We developed a Bayesian movement model that accounts for error from multiple data sources as well as movement behaviour at different temporal scales. The Bayesian framework allows us to calculate derived quantities that describe temporally varying movement behaviour, such as residence time, speed and persistence in direction. The model is flexible, easy to implement and computationally efficient. We apply this model to data from Colorado Canada lynx ( Lynx canadensis ) and use derived quantities to identify changes in movement behaviour.

Methods in Ecology and Evolution↗
USGS2009

Bird population density estimated from acoustic signals

Many animal species are detected primarily by sound. Although songs, calls and other sounds are often used for population assessment, as in bird point counts and hydrophone surveys of cetaceans, there are few rigorous methods for estimating population density from acoustic data. 2. The problem has several parts - distinguishing individuals, adjusting for individuals that are missed, and adjusting for the area sampled. Spatially explicit capture-recapture (SECR) is a statistical methodology that addresses jointly the second and third parts of the problem. We have extended SECR to use uncalibrated information from acoustic signals on the distance to each source. 3. We applied this extension of SECR to data from an acoustic survey of ovenbird Seiurus aurocapilla density in an eastern US deciduous forest with multiple four-microphone arrays. We modelled average power from spectrograms of ovenbird songs measured within a window of 0??7 s duration and frequencies between 4200 and 5200 Hz. 4. The resulting estimates of the density of singing males (0??19 ha -1 SE 0??03 ha -1 ) were consistent with estimates of the adult male population density from mist-netting (0??36 ha -1 SE 0??12 ha -1 ). The fitted model predicts sound attenuation of 0??11 dB m -1 (SE 0??01 dB m -1 ) in excess of losses from spherical spreading. 5.Synthesis and applications. Our method for estimating animal population density from acoustic signals fills a gap in the census methods available for visually cryptic but vocal taxa, including many species of bird and cetacean. The necessary equipment is simple and readily available; as few as two microphones may provide adequate estimates, given spatial replication. The method requires that individuals detected at the same place are acoustically distinguishable and all individuals vocalize during the recording interval, or that the per capita rate of vocalization is known. We believe these requirements can be met, with suitable field methods, for a significant number of songbird species. ?? 2009 British Ecological Society.

Journal of Applied Ecology↗
USGS2020

Time-series model, statistical methods, and software documentation for R–QWTREND—An R package for analyzing trends in stream-water quality

As part of a U.S. Geological Survey water-quality study started in 2018, in cooperation with the International Joint Commission, North Dakota Department of Environmental Quality, and Minnesota Pollution Control Agency, a publicly available software package called R–QWTREND was developed for analyzing trends in stream-water quality. The R–QWTREND package is a collection of functions written in R, an open source language and a general environment for statistical computing and graphics. The package uses a parametric time-series model to express logarithmically transformed concentration in terms of flow-related variability, trend, and serially correlated model errors. Flow-related variability captures natural variability in concentration on the basis of concurrent and antecedent streamflow. The trend identifies systematic changes in concentration in terms of potential step trends, piecewise monotonic trends, or user-specified trends. Maximum likelihood estimation is used to estimate model parameters and determine the best-fit trend model. This report describes the time-series model and statistical methodology behind R–QWTREND and provides formal documentation for installing and using the package.

Open-File Report↗