USGS Science⌕ Search

Geology topics

Mevin Hooten

Publications and source records attributed to Mevin Hooten.

At least 37 records · Page 2Linked to original sources
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↗
USGS2019

Making Recursive Bayesian inference accessible

Bayesian models provide recursive inference naturally because they can formally reconcile new data and existing scientific information. However, popular use of Bayesian methods often avoids priors that are based on exact posterior distributions resulting from former studies. Two existing Recursive Bayesian methods are: Prior- and Proposal-Recursive Bayes. Prior-Recursive Bayes uses Bayesian updating, fitting models to partitions of data sequentially, and provides a way to accommodate new data as they become available using the posterior from the previous stage as the prior in the new stage based on the latest data. ProposalRecursive Bayes is intended for use with hierarchical Bayesian models and uses a set of transient priors in first stage independent analyses of the data partitions. The second stage of Proposal-Recursive Bayes uses the posteriors from the first stage as proposals in an MCMC algorithm to fit the full model. We combine Prior- and Proposal-Recursive concepts to fit any Bayesian model, and often with computational improvements. We demonstrate our method with two case studies. Our approach has implications for big data, streaming data, and optimal adaptive design situations.

American Statistician↗
USGS2019

Extreme site fidelity as an optimal strategy in an unpredictable and homogeneous environment

1.Animal site fidelity structures space use, population demography and ultimately gene flow. Understanding the adaptive selection for site fidelity patterns provides a mechanistic understanding to both spatial and population processes. This can be achieved by linking space use with environmental variability (spatial and temporal) and demographic parameters. However, rarely is the environmental context that drives the selection for site fidelity behaviour fully considered. 2. We use ecological theory to understand whether the spatial and temporal variability in breeding site quality can explain the site fidelity behaviour and demographic patterns of Gunnison sage-grouse (Centrocercus minimus). We examined female site fidelity patterns across multiple spatial scales: proximity of consecutive year nest locations, space-use overlap within and across the breeding and brooding seasons, and fidelity to a breeding patch. We also examined the spatial and temporal variability in nest, chick, juvenile and adult survival. 3. We found Gunnison sage-grouse to be site faithful to their breeding patch, area of use within the patch and generally where they nest, suggesting an “Always Stay” site fidelity strategy. This is an optimal evolutionary strategy when site quality is unpredictable. Further, we found limited spatial variability in survival within age groups, suggesting little demographic benefit to moving among patches. We suggest Gunnison sage-grouse site fidelity is driven by the unpredictability of predation in a relatively homogeneous environment, the lack of benefits and likely costs to moving across landscape patches and leaving known lek and breeding/brooding areas. 4. Space use and demography are commonly studied separately. More so, site fidelity patterns are rarely framed in the context of ecological theory, beyond questions related to the win-stay:lose-switch rule. To move beyond describing patterns and understand the adaptive selection driving species movements and their demographic consequences require integrating movement, demography and environmental variability in a synthetic framework. 5. Site fidelity theory provides a coherent framework to simultaneously investigate the spatial and population ecology of animal populations. Using it to frame ecological questions will lead to a more mechanistic understanding of animal movement, spatial population structuring and meta-population dynamics.

Colorado↗
USGS2019

Bringing Bayesian models to life

Bringing Bayesian Models to Life empowers the reader to extend, enhance, and implement statistical models for ecological and environmental data analysis. We open the black box and show the reader how to connect modern statistical models to computer algorithms. These algorithms allow the user to fit models that answer their scientific questions without needing to rely on automated Bayesian software. We show how to handcraft statistical models that are useful in ecological and environmental science including: linear and generalized linear models, spatial and time series models, occupancy and capture-recapture models, animal movement models, spatio-temporal models, and integrated population-models.

Book↗
USGS2019

Comparing and improving methods for reconstructing peatland water-table depth from testate amoebae

Proxies that use changes in the composition of ecological communities to reconstruct temporal changes in an environmental covariate are commonly used in paleoclimatology and paleolimnology. Existing methods, such as weighted averaging and modern analog technique, relate compositional data to the covariate in very simple ways, and different methods are seldom compared systematically. We present a new Bayesian model that better represents the underlying data and the complexity in the relationships between species’ abundances and a paleoenvironmental covariate. Using testate amoeba-based reconstructions of water-table depth as a test case, we systematically compare new and existing models in a cross-validation experiment on a large training dataset from North America. We then apply the different models to a new 7500-year record of testate amoeba assemblages from Caribou Bog in Maine and compare the resulting water-table depth reconstructions. We find that Bayesian models represent an improvement over existing methods in three key ways: more complete use of the underlying compositional data, full and meaningful treatment of uncertainty, and clear paths toward methodological improvements. Furthermore, we highlight how developing and systematically comparing methods leads to an improved understanding of the proxy system. This paper focuses on testate amoebae and water-table depth, but the framework and ideas are widely applicable to other proxies based on compositional data.

Continental United States and Alaska↗
USGS2019

The rise of an apex predator following deglaciation

Aim Sea otters ( Enhydra lutris ) are an apex predator of the nearshore marine community and nearly went extinct at the turn of the 20th century. Reintroductions and legal protection allowed sea otters to re‐colonize much of their former range. Our objective was to chronicle the colonization of this apex predator in Glacier Bay, Alaska, to help understand the mechanisms that governed their successful colonization. Location Glacier Bay is a tidewater glacier fjord in southeastern Alaska that was entirely covered by glaciers in the mid‐18th century. Since then, it has endured the fastest tidewater glacier retreat in recorded history. Methods We collected and analysed several data sets, spanning 20 years, to document the spatio‐temporal dynamics of an apex predator expanding into an area where they were formerly absent. We used novel quantitative tools to model the occupancy, abundance and colonization dynamics of sea otters, while accounting for uncertainty in the data collection process, the ecological process and model parameters. Results Twenty years after sea otters were first observed colonizing Glacier Bay, they became one of the most abundant and widely distributed marine mammal. The population grew exponentially at a rate of 20% per year. They colonized Glacier Bay at a maximum rate of 6 km per year, with faster colonization rates occurring early in the colonization process. During colonization, sea otters selected shallow areas, close to shore, with a steep bottom slope, and a relatively simple shoreline complexity index. Main conclusions The growth and expansion of sea otters in Glacier Bay demonstrate how legal protection and translocation of apex predators can facilitate their successful establishment into a community in which they were formerly absent. The success of sea otters was, in part, a consequence of habitat that was left largely unperturbed by humans for the past 250 years. Further, sea otters and other marine predators, whose distribution is limited by ice, have the potential to expand in distribution and abundance, reshaping future marine communities in the wake of deglaciation and global loss of sea ice.

Alaska↗
USGS2019

Accounting for phenology in the analysis of animal movement

The analysis of animal tracking data provides important scientific understanding and discovery in ecology. Observations of animal trajectories using telemetry devices provide researchers with information about the way animals interact with their environment and each other. For many species, specific geographical features in the landscape can have a strong effect on behavior. Such features may correspond to a single point (eg, dens or kill sites), or to higher dimensional subspaces (eg, rivers or lakes). Features may be relatively static in time (eg, coastlines or home‐range centers), or may be dynamic (eg, sea ice extent or areas of high‐quality forage for herbivores). We introduce a novel model for animal movement that incorporates active selection for dynamic features in a landscape. Our approach is motivated by the study of polar bear (Ursus maritimus) movement. During the sea ice melt season, polar bears spend much of their time on sea ice above shallow, biologically productive water where they hunt seals. The changing distribution and characteristics of sea ice throughout the year mean that the location of valuable habitat is constantly shifting. We develop a model for the movement of polar bears that accounts for the effect of this important landscape feature. We introduce a two‐stage procedure for approximate Bayesian inference that allows us to analyze over 300 000 observed locations of 186 polar bears from 2012 to 2016. We use our model to estimate a spatial boundary of interest to wildlife managers that separates two subpopulations of polar bears from the Beaufort and Chukchi seas.

Biometrics↗
USGS2019

A hierarchical Bayesian approach for handling missing classification data

Ecologists use classifications of individuals in categories to understand composition of populations and communities. These categories might be defined by demographics, functional traits, or species. Assignment of categories is often imperfect, but frequently treated as observations without error. When individuals are observed but not classified, these “partial” observations must be modified to include the missing data mechanism to avoid spurious inference. We developed two hierarchical Bayesian models to overcome the assumption of perfect assignment to mutually exclusive categories in the multinomial distribution of categorical counts, when classifications are missing. These models incorporate auxiliary information to adjust the posterior distributions of the proportions of membership in categories. In one model, we use an empirical Bayes approach, where a subset of data from one year serves as a prior for the missing data the next. In the other approach, we use a small random sample of data within a year to inform the distribution of the missing data. We performed a simulation to show the bias that occurs when partial observations were ignored and demonstrated the altered inference for the estimation of demographic ratios. We applied our models to demographic classifications of elk ( Cervus elaphus nelsoni ) to demonstrate improved inference for the proportions of sex and stage classes. We developed multiple modeling approaches using a generalizable nested multinomial structure to account for partially observed data that were missing not at random for classification counts. Accounting for classification uncertainty is important to accurately understand the composition of populations and communities in ecological studies.

Ecology and Evolution↗
USGS2019

Spatially-structured statistical network models for landscape genetics

A basic understanding of how the landscape impedes, or creates resistance to, the dispersal of organisms and hence gene flow is paramount for successful conservation science and management. Spatially structured ecological networks are often used to represent spatial landscape‐genetic relationships, where nodes represent individuals or populations and resistance to movement is represented using non‐binary edge weights. Weights are typically assigned or estimated by the user, rather than observed, and validating such weights is challenging. We provide a synthesis of current methods used to estimate edge weights and an overview of common model types, stressing the advantages and disadvantages of each approach and their ability to model landscape‐genetic data. We further explore a set of spatial‐statistical methods that provide ecologists with alternative approaches for modeling spatially explicit processes that may affect genetic structure. This includes an overview of spatial autoregressive models, with a particular focus on how correlation and partial correlation are used to represent neighborhood structure with the inverse of the covariance matrix (i.e., precision matrix). We then demonstrate how to model resistance by specifying an appropriate statistical model on the nodes, conditioned on the edge weights, through the precision matrix. This integration of network ecology and spatial statistics provides a practical analytical framework for landscape‐genetic studies. The results can be used to make statistical inferences about the relative importance of individual landscape characteristics, such as the vegetative cover, hillslope, or the presence of roads or rivers, on gene flow. In addition, the R code we include allows readers to explore landscape‐genetic structure in their own datasets, which will potentially provide new insights into the evolutionary processes that generated ecological networks, as well as valuable information about the optimal characteristics of conservation corridors.

Ecological Monographs↗
USGS2019

Estimating lake–climate responses from sparse data: An application to high elevation lakes

Although many studies demonstrate lake warming, few document trends from lakes with sparse data. Diel and seasonal variability of surface temperatures limit conventional trend analyses to datasets with frequent repeated observations. Thus, remote lakes, including many high elevation lakes, are underrepresented in trend analyses. We used a Bayesian technique to analyze sparse data that explicitly incorporated diel and seasonal variability. This approach allowed us to estimate lake warming in a region of limited knowledge: high elevation lakes (> 2100 m ASL) of the Southern Rocky Mountains, U.S.A. The analysis allowed for inclusion of lakes with few repeated measurements, and observations made before 1980 when more intensive lake monitoring began. We accumulated the largest dataset of high elevation lake temperatures analyzed to date. Data from 590 high elevation lakes in the Southern Rocky Mountains showed a 0.13°C decade −1 increase in surface temperatures and a 14% increase in seasonal degree days since 1955. This result is lower than other regional and global estimates of lake warming; however, it is similar to other high elevation lake studies. Our approach can be applied to other understudied regions, increasing our overall understanding of the effects of climate change on lakes and their temporal dynamics.

Colorado↗
USGS2019

Selecting ecological models using multi-objective optimization

Choices in ecological research and natural resource management require balancing multiple, often competing objectives. Examples include maximizing species persistence in a wildlife conservation context, while minimizing cost, or balancing opposing stakeholder objectives when managing wildlife populations. Multiple-objective optimization (MOO) provides a unifying framework for solving multiple objective problems. Model selection is a critical component of ecological inference and prediction and requires balancing the competing objectives of model fit and model complexity. The tradeoff between model fit and model complexity provides a basis for describing the model-selection problem within the MOO framework. We discuss MOO and two strategies for solving the MOO problem; modeling preferences pre-optimization and post-optimization. Most conventional model selection methods can be formulated as solutions of MOO problems via specification of pre-optimization preferences. We reconcile model selection within the MOO framework. We also consider model selection using post-optimization specification of preferences. That is, by first identifying Pareto optimal solutions, and then selecting among them. We demonstrate concepts with an ecological application of model selection using avian species richness data in the continental United States.

Ecological Modelling↗
USGS2019

On the relationship between conditional (CAR) and simultaneous (SAR) autoregressive models

We clarify relationships between conditional (CAR) and simultaneous (SAR) autoregressive models. We review the literature on this topic and find that it is mostly incomplete. Our main result is that a SAR model can be written as a unique CAR model, and while a CAR model can be written as a SAR model, it is not unique. In fact, we show how any multivariate Gaussian distribution on a finite set of points with a positive-definite covariance matrix can be written as either a CAR or a SAR model. We illustrate how to obtain any number of SAR covariance matrices from a single CAR covariance matrix by using Givens rotation matrices on a simulated example. We also discuss sparseness in the original CAR construction, and for the resulting SAR weights matrix. For a real example, we use crime data in 49 neighborhoods from Columbus, Ohio, and show that a geostatistical model optimizes the likelihood much better than typical first-order CAR models. We then use the implied weights from the geostatistical model to estimate CAR model parameters that provides the best overall optimization.

Spatial Statistics↗
USGS2019

Predicting paleoclimate from compositional data using multivariate Gaussian process inverse prediction

Multivariate compositional count data arise in many applications including ecology, microbiology, genetics and paleoclimate. A frequent question in the analysis of multivariate compositional count data is what underlying values of a covariate(s) give rise to the observed composition. Learning the relationship between covariates and the compositional count allows for inverse prediction of unobserved covariates given compositional count observations. Gaussian processes provide a flexible framework for modeling functional responses with respect to a covariate without assuming a functional form. Many scientific disciplines use Gaussian process approximations to improve prediction and make inference on latent processes and parameters. When prediction is desired on unobserved covariates given realizations of the response variable, this is called inverse prediction. Because inverse prediction is often mathematically and computationally challenging, predicting unobserved covariates often requires fitting models that are different from the hypothesized generative model. We present a novel computational framework that allows for efficient inverse prediction using a Gaussian process approximation to generative models. Our framework enables scientific learning about how the latent processes co-vary with respect to covariates while simultaneously providing predictions of missing covariates. The proposed framework is capable of efficiently exploring the high dimensional, multi-modal latent spaces that arise in the inverse problem. To demonstrate flexibility, we apply our method in a generalized linear model framework to predict latent climate states given multivariate count data. Based on cross-validation, our model has predictive skill competitive with current methods while simultaneously providing formal, statistical inference on the underlying community dynamics of the biological system previously not available.

Annals of Applied Statistics↗
USGS2018

Running on empty: Recharge dynamics from animal movement data

Vital rates such as survival and recruitment have always been important in the study of population and community ecology. At the individual level, physiological processes such as energetics are critical in understanding biomechanics and movement ecology and also scale up to influence food webs and trophic cascades. Although vital rates and population-level characteristics are tied with individual-level animal movement, most statistical models for telemetry data are not equipped to provide inference about these relationships because they lack the explicit, mechanistic connection to physiological dynamics. We present a framework for modelling telemetry data that explicitly includes an aggregated physiological process associated with decision making and movement in heterogeneous environments. Our framework accommodates a wide range of movement and physiological process specifications. We illustrate a specific model formulation in continuous-time to provide direct inference about gains and losses associated with physiological processes based on movement. Our approach can also be extended to accommodate auxiliary data when available. We demonstrate our model to infer mountain lion ( Puma concolor ; in Colorado, USA) and African buffalo ( Syncerus caffer ; in Kruger National Park, South Africa) recharge dynamics.

Ecology Letters↗
USGS2018

Time-varying predatory behavior is primary predictor of fine-scale movement of wildland-urban cougars

Background While many species have suffered from the detrimental impacts of increasing human population growth, some species, such as cougars ( Puma concolor ), have been observed using human-modified landscapes. However, human-modified habitat can be a source of both increased risk and increased food availability, particularly for large carnivores. Assessing preferential use of the landscape is important for managing wildlife and can be particularly useful in transitional habitats, such as at the wildland-urban interface. Preferential use is often evaluated using resource selection functions (RSFs), which are focused on quantifying habitat preference using either a temporally static framework or researcher-defined temporal delineations. Many applications of RSFs do not incorporate time-varying landscape availability or temporally-varying behavior, which may mask conflict and avoidance behavior. Methods Contemporary approaches to incorporate landscape availability into the assessment of habitat selection include spatio-temporal point process models, step selection functions, and continuous-time Markov chain (CTMC) models; in contrast with the other methods, the CTMC model allows for explicit inference on animal movement in continuous-time. We used a hierarchical version of the CTMC framework to model speed and directionality of fine-scale movement by a population of cougars inhabiting the Front Range of Colorado, U.S.A., an area exhibiting rapid population growth and increased recreational use, as a function of individual variation and time-varying responses to landscape covariates. Results We found evidence for individual- and daily temporal-variability in cougar response to landscape characteristics. Distance to nearest kill site emerged as the most important driver of movement at a population-level. We also detected seasonal differences in average response to elevation, heat loading, and distance to roads. Motility was also a function of amount of development, with cougars moving faster in developed areas than in undeveloped areas. Conclusions The time-varying framework allowed us to detect temporal variability that would be masked in a generalized linear model, and improved the within-sample predictive ability of the model. The high degree of individual variation suggests that, if agencies want to minimize human-wildlife conflict management options should be varied and flexible. However, due to the effect of recursive behavior on cougar movement, likely related to the location and timing of potential kill-sites, kill-site identification tools may be useful for identifying areas of potential conflict.

Colorado↗
USGS2018

Accounting for location uncertainty in azimuthaltelemetry data improves ecological inference

Background Characterizing animal space use is critical for understanding ecological relationships. Animal telemetry technology has revolutionized the fields of ecology and conservation biology by providing high quality spatial data on animal movement. Radio-telemetry with very high frequency (VHF) radio signals continues to be a useful technology because of its low cost, miniaturization, and low battery requirements. Despite a number of statistical developments synthetically integrating animal location estimation and uncertainty with spatial process models using satellite telemetry data, we are unaware of similar developments for azimuthal telemetry data. As such, there are few statistical options to handle these unique data and no synthetic framework for modeling animal location uncertainty and accounting for it in ecological models. We developed a hierarchical modeling framework to provide robust animal location estimates from one or more intersecting or non-intersecting azimuths. We used our azimuthal telemetry model (ATM) to account for azimuthal uncertainty with covariates and propagate location uncertainty into spatial ecological models. We evaluate the ATM with commonly used estimators (Lenth (1981) maximum likelihood and M-Estimators) using simulation. We also provide illustrative empirical examples, demonstrating the impact of ignoring location uncertainty within home range and resource selection analyses. We further use simulation to better understand the relationship among location uncertainty, spatial covariate autocorrelation, and resource selection inference. Results We found the ATM to have good performance in estimating locations and the only model that has appropriate measures of coverage. Ignoring animal location uncertainty when estimating resource selection or home ranges can have pernicious effects on ecological inference. Home range estimates can be overly confident and conservative when ignoring location uncertainty and resource selection coefficients can lead to incorrect inference and over confidence in the magnitude of selection. Furthermore, our simulation study clarified that incorporating location uncertainty helps reduce bias in resource selection coefficients across all levels of covariate spatial autocorrelation. Conclusion The ATM can accommodate one or more azimuths when estimating animal locations, regardless of how they intersect; this ensures that all data collected are used for ecological inference. Our findings and model development have important implications for interpreting historical analyses using this type of data and the future design of radio-telemetry studies.

Movement Ecology↗
USGS2018

A guide to Bayesian model checking for ecologists

Checking that models adequately represent data is an essential component of applied statistical inference. Ecologists increasingly use hierarchical Bayesian statistical models in their research. The appeal of this modeling paradigm is undeniable, as researchers can build and fit models that embody complex ecological processes while simultaneously accounting for observation error. However, ecologists tend to be less focused on checking model assumptions and assessing potential lack of fit when applying Bayesian methods than when applying more traditional modes of inference such as maximum likelihood. There are also multiple ways of assessing the fit of Bayesian models, each of which has strengths and weaknesses. For instance, Bayesian P values are relatively easy to compute, but are well known to be conservative, producing P values biased toward 0.5. Alternatively, lesser known approaches to model checking, such as prior predictive checks, cross‐validation probability integral transforms, and pivot discrepancy measures may produce more accurate characterizations of goodness‐of‐fit but are not as well known to ecologists. In addition, a suite of visual and targeted diagnostics can be used to examine violations of different model assumptions and lack of fit at different levels of the modeling hierarchy, and to check for residual temporal or spatial autocorrelation. In this review, we synthesize existing literature to guide ecologists through the many available options for Bayesian model checking. We illustrate methods and procedures with several ecological case studies including (1) analysis of simulated spatiotemporal count data, (2) N‐mixture models for estimating abundance of sea otters from an aircraft, and (3) hidden Markov modeling to describe attendance patterns of California sea lion mothers on a rookery. We find that commonly used procedures based on posterior predictive P values detect extreme model inadequacy, but often do not detect more subtle cases of lack of fit. Tests based on cross‐validation and pivot discrepancy measures (including the “sampled predictive P value”) appear to be better suited to model checking and to have better overall statistical performance. We conclude that model checking is necessary to ensure that scientific inference is well founded. As an essential component of scientific discovery, it should accompany most Bayesian analyses presented in the literature.

Ecological Monographs↗