SEARCH · Geology Reports

Results for “Ecological Modeling”

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

SourceResearch connectionFrom yearThrough year

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

At least 865 records · Page 48Linked to original sources

USGS2017

Species distributions models in wildlife planning: agricultural policy and wildlife management in the great plains

We know economic and social policy has implications for ecosystems at large, but the consequences for a given geographic area or specific wildlife population are more difficult to conceptualize and communicate. Species distribution models, which extrapolate species-habitat relationships across ecological scales, are capable of predicting population changes in distribution and abundance in response to management and policy, and thus, are an ideal means for facilitating proactive management within a larger policy framework. To illustrate the capabilities of species distribution modeling in scenario planning for wildlife populations, we projected an existing distribution model for ring-necked pheasants ( Phasianus colchicus ) onto a series of alternative future landscape scenarios for Nebraska, USA. Based on our scenarios, we qualitatively and quantitatively estimated the effects of agricultural policy decisions on pheasant populations across Nebraska, in specific management regions, and at wildlife management areas.

Nebraska↗

USGS2018

Waterfowl spring migratory behavior and avian influenza transmission risk in the changing landscape of the East Asian-Australasian Flyway

Avian influenza has advanced from a regional concern to a global health issue with significant economic, trade, and public health implications. Wild birds, particularly waterfowl (Anseriformes), are known reservoirs for low-pathogenic avian influenza viruses (AIV) and recent studies have shown their potential in the spread of highly pathogenic forms of virus. East Asia remains an epicenter for the emergence of novel strains of AIV, however, information on movement ecology of waterfowl, and subsequently a clearer understanding of disease transmission risks in this region has been greatly lacking. To address this, we marked two species of wild waterfowl, northern pintail (Anas acuta) and Eurasian wigeon (Anas penelope) , with satellite transmitters on their wintering grounds in Hong Kong, China to study the northward spring migration in the East Asian-Australasian Flyway in relation to disease transmission factors. Northern pintail were found to initiate migration 42 days earlier, travel 2,150 km farther, and perform 4.4 more stopovers than Eurasian wigeon. We found both species used similar stopover locations including areas along the Yangtze River near Shanghai, Bohai Bay and Korea Bay in rapidly developing regions of the Yellow Sea, and the Sea of Okhotsk where the species appeared to funnel through a migratory bottleneck. Both species appeared to exhibit strong habitat selection for rice paddies during migration stopovers, a habitat preference which has the potential to influence risks of AIV outbreaks as rapid land use and land cover changes occur throughout China. Both species had greatest association with H5N1 outbreaks during the early stages of migration when they were at lower latitudes. While Eurasian wigeon were not associated with outbreaks after the mean date of wintering ground departures, northern pintail were associated with outbreaks until the majority of individuals departed from the Yellow Sea, a migratory stopover location. Our results show species-level differences in migration timing and behavior for these common and widespread species, demonstrating the need to consider their unique temporal and spatial movement ecology when incorporating wild birds into AIV risk modeling and management.

Frontiers in Ecology and Evolution↗

USGS2021

Instrumental variable methods in structural equation models

Instrumental variable regression (RegIV) provides a means for detecting and correcting parameter bias in causal models. Widely used in economics, recently several papers have highlighted its potential utility for ecological applications. Little attention has thus far been paid to the fact that IV methods can also be implemented within structural equation models (SEMIV). In this paper I present the motivations, requirements and basic procedures for using SEMIV. I first consider causal inference and IVs from the perspective of a randomized experiment with partial control of the cause of interest. I consider common sources of bias, the role of randomization and limits to its capacity to exclude bias. Sources of bias include omitted confounders, reciprocal causation, reverse causation and measurement error, all of which can all be seen as a single problem—endogeneity. The approach to estimating IV models most commonly used in econometric practice, two-stage least squares regression (2SLS), is explained, followed by a brief exposition of the covariance modelling approach used in SEM. Using data from an ecological field experiment, I illustrate the use of the treatment variable as an IV and then illustrate procedures for evaluating candidate variables that might serve as additional IVs. IV methods are shown to be useful for both detecting endogeneity and removing its influences. I illustrate some of the ways that bias can be generated, as well as diagnostic capabilities and means for remedy embedded within SEM. Procedures for screening and evaluating additional IVs reveal valuable lessons regarding the theoretical requirements and empirical standards for IVs. SEMIV provides a useful way to detect and control for bias. I suggest that the use of IVs within the SEM framework can support the simultaneous pursuit of causal inference and explanatory modelling, a common pair of aspirations for ecologists. Moving forward, there is a need for a better understanding of the capabilities of SEMIV and requirements for successful application.

Methods in Ecology and Evolution↗

USGS2013

A spatial mark–resight model augmented with telemetry data

Abundance and population density are fundamental pieces of information for population ecology and species conservation, but they are difficult to estimate for rare and elusive species. Mark-resight models are popular for estimating population abundance because they are less invasive and expensive than traditional mark-recapture. However, density estimation using mark-resight is difficult because the area sampled must be explicitly defined, historically using ad-hoc approaches. We develop a spatial mark-resight model for estimating population density that combines spatial resighting data and telemetry data. Incorporating telemetry data allows us to inform model parameters related to movement and individual location. Our model also allows <100% individual identification of marked individuals. We implemented the model in a Bayesian framework, using a custom-made Metropolis-within-Gibbs Markov chain Monte Carlo algorithm. As an example, we applied this model to a mark-resight study of raccoons on South Core Banks, a barrier island in Cape Lookout National Seashore, North Carolina. We estimated a population of 186.71 ± 14.81 individuals, which translated to a density of 8.29 ± 0.66 individuals per km 2 . The model presented here will have widespread utility in future applications, especially for species that are not naturally marked.

Ecology↗

USGS2021

Predicting non-native insect impact: Focusing on the trees to see the forest

Non-native organisms have invaded novel ecosystems for centuries, yet we have only a limited understanding of why their impacts vary widely from minor to severe. Predicting the impact of non-established or newly detected species could help focus biosecurity measures on species with the highest potential to cause widespread damage. However, predictive models require an understanding of potential drivers of impact and the appropriate level at which these drivers should be evaluated. Here, we used non-native, specialist herbivorous insects of forest ecosystems to test which factors drive impact and if there were differences based on whether they used woody angiosperms or conifers as hosts. We identified convergent and divergent patterns between the two host types indicating fundamental similarities and differences in their interactions with non-native insects. Evolutionary divergence time between native and novel hosts was a significant driver of insect impact for both host types but was modulated by different factors in the two systems. Beetles in the subfamily Scolytinae posed the highest risk to woody angiosperms, and different host traits influenced impact of specialists on conifers and woody angiosperms. Tree wood density was a significant predictor of host impact for woody angiosperms with intermediate densities (0.5–0.6 mg/mm 3 ) associated with highest risk, whereas risk of impact was highest for conifers that coupled shade tolerance with drought intolerance. These results underscore the importance of identifying the relevant levels of biological organization and ecological interactions needed to develop accurate risk models for species that may arrive in novel ecosystems.

Biological Invasions↗

USGS2006

Multistate survival models and their extensions in Program MARK

Program MARK provides .100 models for the estimation of population parameters from mark?encounter data. The multistate model of Brownie et al. (1993) and Hestbeck et al. (1991) allows animals to move between states with a probability of transition. The simplest multistate model is an extension of the Cormack?Jolly?Seber (CJS) live recapture model. arameters estimated are state-specific survival rates and encounter probabilities and transition probabilities between states. The multistate model provides a valuable framework to evaluate important ecological questions. For example, estimation of state-specific survival and transition probabilities between the biological states of breeders and nonbreeders allows estimation of the cost of reproduction. Transitions between physical states, such as spatial areas, provide estimates needed for meta-population models. The basic multistate model uses only live recaptures, but 3 extensions are included in MARK. A multistate model with live and dead encounters is available, although the dead encounters are not state specific. Robust-design multistate models are also included in MARK, with both open and closed robust designs. These models assume that animals move between states only between primary sessions of the robust design. For the closed robust design, we can specify 12 different data types for the modeling of encounter probabilities during the primary session, including 6 versions of the closed model likelihood incorporating population size (N) directly in the likelihood, and 6 versions of the Huggins model in which N is estimated as a derived parameter outside the likelihood. One assumption that is generally necessary to estimate state-specific survival rates in the multistate model is that transitions take place immediately before encounter occasions. Otherwise, survival rates over the interval between encounter occasions are a mix of survival rates over multiple states. Advantages of using MARK to estimate the parameters of the various multistate models include flexibility of model specification to include group, time, and individual covariates, estimation of variance components, model averaging of parameter estimates, and Bayesian parameter estimation using Markov chain Monte Carlo procedures on the logit scale.

Journal of Wildlife Management↗

USGS1996

Potential role of vegetation feedback in the climate sensitivity of high-latitude regions: A case study at 6000 years B.P.

Previous climate model simulations have shown that the configuration of the Earth's orbit during the early to mid-Holocene (approximately 10–5 kyr) can account for the generally warmer-than-present conditions experienced by the high latitudes of the northern hemisphere. New simulations for 6 kyr with two atmospheric/mixed-layer ocean models (Community Climate Model, version 1, CCMl, and Global ENvironmental and Ecological Simulation of Interactive Systems, version 2, GENESIS 2) are presented here and compared with results from two previous simulations with GENESIS 1 that were obtained with and without the albedo feedback due to climate-induced poleward expansion of the boreal forest. The climate model results are summarized in the form of potential vegetation maps obtained with the global BIOME model, which facilitates visual comparisons both among models and with pollen and plant macrofossil data recording shifts of the forest-tundra boundary. A preliminary synthesis shows that the forest limit was shifted 100–200 km north in most sectors. Both CCMl and GENESIS 2 produced a shift of this magnitude. GENESIS 1 however produced too small a shift, except when the boreal forest albedo feedback was included. The feedback in this case was estimated to have amplified forest expansion by approximately 50%. The forest limit changes also show meridional patterns (greatest expansion in central Siberia and little or none in Alaska and Labrador) which have yet to be reproduced by models. Further progress in understanding of the processes involved in the response of climate and vegetation to orbital forcing will require both the deployment of coupled atmosphere-biosphere-ocean models and the development of more comprehensive observational data sets.

Global Biogeochemical Cycles↗

USGS2025

Heterogeneity of locked‐pasture snow conditions modulate habitat and movement choices of a facultative migrant

Habitat selection and movement are key mechanisms by which animals can respond to and potentially cope with highly variable environmental conditions. Optimal responses likely vary, however, depending on the severity and scope of conditions. We tested this hypothesis using a facultative migrant species, the Great Gray Owl ( Strix nebulosa ), which exhibits high inter- and intra-individual variation in the timing, direction, and distance of winter movements. Specifically, we evaluated whether episodic, spatiotemporally variable “locked-pasture” snow conditions, which restrict access to subnivean food, prompted shifts in habitat selection or long-distance movements by owls. We quantified the movement of 42 owls using global positioning system (GPS) data within the Greater Yellowstone Ecosystem, USA, during 2017–2022. We used a novel ecological application of SnowModel, a snow evolution modeling system, to estimate fine-scale, physical snow properties likely to influence access to prey. Variables included snow depth, snow crusts produced by wind, and ice crusts produced by melt-freeze and rain-on-snow events. Owls avoided heterogeneously distributed wind crusts via local shifts in habitat selection. More homogenous ice crusts elicited long-distance movements away from affected home ranges. Finally, owls employed both proximate shifts in habitat selection and long-distance movements to avoid deeper snow. Ultimately, owls exhibited behavioral flexibility in response to limiting snow conditions that can vary in terms of severity, spatial extent, and duration. Such behavioral responses determine species distribution, with implications for population and community dynamics in spatiotemporally variable systems. Understanding the effects of, and responses to, environmental controls is increasingly important given the scope of on-going global change.

Idaho, Wyoming↗

USGS2023

Towards a better future for biodiversity and people: Modelling Nature Futures

The Nature Futures Framework (NFF) is a heuristic tool for co-creating positive futures for nature and people. It seeks to open up a diversity of futures through mainly three value perspectives on nature – Nature for Nature, Nature for Society, and Nature as Culture. This paper describes how the NFF can be applied in modelling to support decision-making. First, we describe key considerations for the NFF in developing qualitative and quantitative scenarios: i) multiple value perspectives on nature as a state space where pathways improving nature toward a frontier can be represented, ii) mutually reinforcing key feedbacks of social-ecological systems that are important for nature conservation and human wellbeing, iii) indicators of multiple knowledge systems describing the evolution of complex social-ecological dynamics. We then present three approaches to modelling Nature Futures scenarios in the review, screening, and design phases of policy processes. This paper seeks to facilitate the integration of relational values of nature in models and strengthen modelled linkages across biodiversity, nature’s contributions to people, and quality of life.

Global Environmental Change↗

USGS2024

A hierarchical model for eDNA fate and transport dynamics accommodating low concentration samples

Environmental DNA (eDNA) sampling is an increasingly important tool for answering ecological questions and informing aquatic species management; however, several factors currently limit the reliability of ecological inference from eDNA sampling. Two particular challenges are (1) determining species source location(s) and (2) accurately and precisely measuring low concentration eDNA samples in the presence of multiple sources of ecological and measurement variability. The recently introduced eDNA Integrating Transport and Hydrology (eDITH) model provides a framework for relating eDNA measurements to source locations in riverine networks, but little empirical work has been done to test and refine model assumptions or accommodate low concentration samples, that can be systematically undermeasured. To better understand eDNA fate and transport dynamics and our ability to reliably quantify low concentration samples, we developed a hierarchical model and used it to evaluate a fate and transport experiment. Our model addresses several low concentration challenges by modeling the number of copies in each PCR replicate as a latent variable with a count distribution and conditioning detection and quantification on replicate copy number. We provide evidence that the eDNA removal rate declined through time, estimating that over 80% of eDNA was removed over the first 10 m, traversed in 41 s. After this initial period of rapid decay, eDNA decayed slowly with consistent detection through our farthest site 1 km from the release location, traversed in 67.8 min. Our model further allowed us to detect extra-Poisson variation in the allocation of copies to replicates. We extended our hierarchical model to accommodate a continuous effect of inhibitors and used our model to provide evidence for the inhibitor hypothesis and explore the potential implications. While our model is not a panacea for all challenges faced when quantifying low-concentration eDNA samples, it provides a framework for a more complete accounting of uncertainty.

Environmental and Ecological Statistics↗

USGS2007

Invasive plants and their ecological strategies: Prediction and explanation of woody plant invasion in New England

Effective management of introduced species requires the early identification of species that pose a significant threat of becoming invasive. To better understand the invasive ecology of species in New England, USA, we compiled a character data set with which to compare non-native species that are known invaders to non-native species that are not currently known to be invasive. In contrast to previous biological trait-based models, we employed a Bayesian hierarchical analysis to identify sets of plant traits associated with invasiveness for each of three growth forms (vines, shrubs, and trees). The resulting models identify a suite of 'invasive traits' highlighting the ecology associated with invasiveness for each of three growth forms. The most effective predictors of invasiveness that emerged from our model were 'invasive elsewhere', 'fast growth rate', 'native latitudinal range', and 'growth form'. The contrast among growth forms was pronounced. For example, 'wind dispersal' was positively correlated with invasiveness in trees, but negatively correlated in shrubs and vines. The predictive model was able to correctly classify invasive plants 67% of the time (22/33), and non-invasive plants 95% of the time (204/215). A number of potential future invasive species in New England that deserve management consideration were identified. ?? 2007 The Authors.

Diversity and Distributions↗

USGS2021

Synthesis of data and studies relating to Delta Smelt biology in the San Francisco Estuary, emphasizing water year 2017

In the San Francisco Estuary (SFE), the effects of freshwater flow on the aquatic ecosystem have been studied extensively over the years and remains a contentious management issue. It is especially contentious with regards to the Delta Smelt (Hypomesus transpacificus), a species endemic to the SFE that has been listed as threatened under the Federal Endangered Species Act and endangered by the State of California. Early studies of Delta Smelt distribution within the SFE suggested that Delta Smelt habitat is determined largely by freshwater flow; however, the exact mechanisms and processes producing such benefits remained unclear. In the summer of 2017, the Flow Alteration Management, Analysis, and Synthesis Team (FLOAT-MAST) was established to analyze, synthesize, and summarize the data collected from the various flow-related monitoring and special studies occurring in 2017(see Table Intro 4). This report will focus on the 2017 summer-fall status of Delta Smelt and its habitat following a record wet year. There has been a long-term decline in the abundance of Delta Smelt associated with a decline in other pelagic fishes. Investigators concluded that the decline has likely been caused by the interactive effects of several causes, including changes in both physical and biotic habitats, many of which are tied to amount and timing of freshwater flow. For this report, we formulated a number of basic predictions about the likely effects of high flows in 2017 on Delta Smelt and their habitat (Table 3). We use a qualitative weight of evidence approach to evaluate whether these predictions were supported by available data. Data sources included a variety of long-term monitoring surveys conducted by Interagency Ecological Program (IEP) agencies, as well as model outputs. Delta Smelt population, health, and life history metrics rarely responded as predicted. Water temperature appears to have a stronger effect on Delta Smelt growth rate and some metrics of life history diversity than outflow or X2 position. Other life history diversity attributes varied but did not appear to be driven by outflow or temperature. Health status was difficult to interpret. Low prevalence of lesions and improved nutritional condition during the drought was contradicted by declining overall population levels. Because of the sparse catches of Delta Smelt in the post-POD years, we consider the data insufficient to reach firm conclusions about the predictions concerning range and distribution of Delta Smelt, especially in the fall. The prediction of high survival was not supported. The 2017 Delta Smelt year class began with poor recruitment in spring of 2017 and below average survival for spring to summer and summer to fall. Thus, low production and low survival led to low abundance of all life stages. During the fall to winter period survival improved, yet the resulting adults were low in number. Foraging success of the fish captured, as measured by stomach fullness, was high for juveniles and adults in 2017 relative to recent years associated with the higher densities of common zooplankton prey that occurred in 2017.

California↗

USGS2017

Linking spring phenology with mechanistic models of host movement to predict disease transmission risk

Disease models typically focus on temporal dynamics of infection, while often neglecting environmental processes that determine host movement. In many systems, however, temporal disease dynamics may be slow compared to the scale at which environmental conditions alter host space-use and accelerate disease transmission. Using a mechanistic movement modelling approach, we made space-use predictions of a mobile host (elk [ Cervus Canadensis ] carrying the bacterial disease brucellosis) under environmental conditions that change daily and annually (e.g., plant phenology, snow depth), and we used these predictions to infer how spring phenology influences the risk of brucellosis transmission from elk (through aborted foetuses) to livestock in the Greater Yellowstone Ecosystem. Using data from 288 female elk monitored with GPS collars, we fit step selection functions (SSFs) during the spring abortion season and then implemented a master equation approach to translate SSFs into predictions of daily elk distribution for five plausible winter weather scenarios (from a heavy snow, to an extreme winter drought year). We predicted abortion events by combining elk distributions with empirical estimates of daily abortion rates, spatially varying elk seroprevelance and elk population counts. Our results reveal strong spatial variation in disease transmission risk at daily and annual scales that is strongly governed by variation in host movement in response to spring phenology. For example, in comparison with an average snow year, years with early snowmelt are predicted to have 64% of the abortions occurring on feedgrounds shift to occurring on mainly public lands, and to a lesser extent on private lands. Synthesis and applications . Linking mechanistic models of host movement with disease dynamics leads to a novel bridge between movement and disease ecology. Our analysis framework offers new avenues for predicting disease spread, while providing managers tools to proactively mitigate risks posed by mobile disease hosts. More broadly, we demonstrate how mechanistic movement models can provide predictions of ecological conditions that are consistent with climate change but may be more extreme than has been observed historically.

Journal of Applied Ecology↗

USGS2022

A quantitative soil-geomorphic framework for developing and mapping ecological site groups

Land management decisions need context about how landscapes will respond to different circumstances or actions. As ecologists’ understanding of nonlinear ecological dynamics has evolved into state-and-transition models (STMs), they have put more emphasis on defining and mapping the soil, geomorphological, and climate parameters that mediate these dynamics. The US Department of Agriculture Natural Resources Conservation Service ecological site descriptions (ESDs) have become the foremost system in classifying lands into ecological units based on STMs. However, an exhaustive inventory of ESDs has proved challenging to complete in the United States, and there have been questions about the consistency of detail in areas completed and the ability to objectively support some assertions made in existing ESDs. To address these issues, this study examines ESDs in the diverse Upper Colorado River region, where ESDs are only partially complete, to look at quantitative approaches to generalizing ecological site concepts based on unifying underlying soil, geomorphology, and climate patterns. Using existing ESDs and vegetation monitoring plot data, results show that a simple hierarchical soil geomorphic unit (SGU) framework based on topographic mediation of moisture, soil salinity, soil depth, slope, rock content, and soil texture can represent much of the ecological dynamics cataloged in ESDs. Analyses of reference plant production data, ecological state attribution, and regional monitoring data show that the new SGUs represent more variation than common climate parameters. This study also included predictively mapping SGUs at 30-m resolution (Kappa of 0.53, 74% agreement with top two predictions in validation). An optimized combination of SGUs with climate zones derived from an aridity index and maximum temperature of the hottest month resulted in an ecological site group framework that condensed over 826 unique ecological site records at various stages of completeness in the regional soil survey down to 35 intuitive and mappable ecological site groups.

Rangeland Ecology and Management↗

USGS2017

A hierarchical model for estimating the spatial distribution and abundance of animals detected by continuous-time recorders

Motivation Several spatial capture-recapture (SCR) models have been developed to estimate animal abundance by analyzing the detections of individuals in a spatial array of traps. Most of these models do not use the actual dates and times of detection, even though this information is readily available when using continuous-time recorders, such as microphones or motion-activated cameras. Instead most SCR models either partition the period of trap operation into a set of subjectively chosen discrete intervals and ignore multiple detections of the same individual within each interval, or they simply use the frequency of detections during the period of trap operation and ignore the observed times of detection. Both practices make inefficient use of potentially important information in the data. Model and data analysis We developed a hierarchical SCR model to estimate the spatial distribution and abundance of animals detected with continuous-time recorders. Our model includes two kinds of point processes: a spatial process to specify the distribution of latent activity centers of individuals within the region of sampling and a temporal process to specify temporal patterns in the detections of individuals. We illustrated this SCR model by analyzing spatial and temporal patterns evident in the camera-trap detections of tigers living in and around the Nagarahole Tiger Reserve in India. We also conducted a simulation study to examine the performance of our model when analyzing data sets of greater complexity than the tiger data. Benefits Our approach provides three important benefits: First, it exploits all of the information in SCR data obtained using continuous-time recorders. Second, it is sufficiently versatile to allow the effects of both space use and behavior of animals to be specified as functions of covariates that vary over space and time. Third, it allows both the spatial distribution and abundance of individuals to be estimated, effectively providing a species distribution model, even in cases where spatial covariates of abundance are unknown or unavailable. We illustrated these benefits in the analysis of our data, which allowed us to quantify differences between nocturnal and diurnal activities of tigers and to estimate their spatial distribution and abundance across the study area. Our continuous-time SCR model allows an analyst to specify many of the ecological processes thought to be involved in the distribution, movement, and behavior of animals detected in a spatial trapping array of continuous-time recorders. We plan to extend this model to estimate the population dynamics of animals detected during multiple years of SCR surveys.

Western Ghats, Nagarahole Tiger Reserve↗

USGS2021

Bayesian modeling can facilitate adaptive management in restoration

There is an urgent need for near-term predictions of ecological restoration outcomes despite imperfect knowledge of ecosystems. Restoration outcomes are always uncertain but integrating Bayesian modeling into the process of adaptive management allows researchers and practitioners to explicitly incorporate prior knowledge of ecosystems into future predictions. Although barriers exist, employing qualitative expert knowledge and previous case studies can help narrow the range of uncertainty in forecasts. Software and processes that allow for repeatable methodologies can help bridge the existing gap between theory and application of Bayesian methods in adaptive management.

Restoration Ecology↗

USGS2012

An accessible method for implementing hierarchical models with spatio-temporal abundance data

A common goal in ecology and wildlife management is to determine the causes of variation in population dynamics over long periods of time and across large spatial scales. Many assumptions must nevertheless be overcome to make appropriate inference about spatio-temporal variation in population dynamics, such as autocorrelation among data points, excess zeros, and observation error in count data. To address these issues, many scientists and statisticians have recommended the use of Bayesian hierarchical models. Unfortunately, hierarchical statistical models remain somewhat difficult to use because of the necessary quantitative background needed to implement them, or because of the computational demands of using Markov Chain Monte Carlo algorithms to estimate parameters. Fortunately, new tools have recently been developed that make it more feasible for wildlife biologists to fit sophisticated hierarchical Bayesian models (i.e., Integrated Nested Laplace Approximation, ‘INLA’). We present a case study using two important game species in North America, the lesser and greater scaup, to demonstrate how INLA can be used to estimate the parameters in a hierarchical model that decouples observation error from process variation, and accounts for unknown sources of excess zeros as well as spatial and temporal dependence in the data. Ultimately, our goal was to make unbiased inference about spatial variation in population trends over time.

PLoS ONE↗

USGS2002

Hierarchical modeling of population stability and species group attributes from survey data

Many ecological studies require analysis of collections of estimates. For example, population change is routinely estimated for many species from surveys such as the North American Breeding Bird Survey (BBS), and the species are grouped and used in comparative analyses. We developed a hierarchical model for estimation of group attributes from a collection of estimates of population trend. The model uses information from predefined groups of species to provide a context and to supplement data for individual species; summaries of group attributes are improved by statistical methods that simultaneously analyze collections of trend estimates. The model is Bayesian; trends are treated as random variables rather than fixed parameters. We use Markov Chain Monte Carlo (MCMC) methods to fit the model. Standard assessments of population stability cannot distinguish magnitude of trend and statistical significance of trend estimates, but the hierarchical model allows us to legitimately describe the probability that a trend is within given bounds. Thus we define population stability in terms of the probability that the magnitude of population change for a species is less than or equal to a predefined threshold. We applied the model to estimates of trend for 399 species from the BBS to estimate the proportion of species with increasing populations and to identify species with unstable populations. Analyses are presented for the collection of all species and for 12 species groups commonly used in BBS summaries. Overall, we estimated that 49% of species in the BBS have positive trends and 33 species have unstable populations. However, the proportion of species with increasing trends differs among habitat groups, with grassland birds having only 19% of species with positive trend estimates and wetland birds having 68% of species with positive trend estimates.

Ecology↗