Search USGS⌕ Search

SEARCH · Search USGS

Results for “Ecological Modeling”

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 973 records · Page 54Linked to original sources

Modeled climate-induced glacier change in Glacier National Park, 1850-2100

The glaciers in the Blackfoot–Jackson Glacier Basin of Glacier National Park, Montana, decreased in area from 21.6 square kilometers (km 2 ) in 1850 to 7.4 km 2 in 1979. Over this same period global temperatures increased by 0.45°C (± 0.15°C). We analyzed the climatic causes and ecological consequences of glacier retreat by creating spatially explicit models of the creation and ablation of glaciers and of the response of vegetation to climate change. We determined the melt rate and spatial distribution of glaciers under two possible future climate scenarios, one based on carbon dioxide–induced global warming and the other on a linear temperature extrapolation. Under the former scenario, all glaciers in the basin will disappear by the year 2030, despite predicted increases in precipitation; under the latter, melting is slower. Using a second model, we analyzed vegetation responses to variations in soil moisture and increasing temperature in a complex alpine landscape and predicted where plant communities are likely to be located as conditions change.

Montana↗

A methodology for ecosystem-scale modeling of selenium

The main route of exposure for selenium (Se) is dietary, yet regulations lack biologically based protocols for evaluations of risk. We propose here an ecosystem-scale model that conceptualizes and quantifies the variables that determinehow Se is processed from water through diet to predators. This approach uses biogeochemical and physiological factors from laboratory and field studies and considers loading, speciation, transformation to particulate material, bioavailability, bioaccumulation in invertebrates, and trophic transfer to predators. Validation of the model is through data sets from 29 historic and recent field case studies of Se-exposed sites. The model links Se concentrations across media (water, particulate, tissue of different food web species). It can be used to forecast toxicity under different management or regulatory proposals or as a methodology for translating a fish-tissue (or other predator tissue) Se concentration guideline to a dissolved Se concentration. The model illustrates some critical aspects of implementing a tissue criterion: 1) the choice of fish species determines the food web through which Se should be modeled, 2) the choice of food web is critical because the particulate material to prey kinetics of bioaccumulation differs widely among invertebrates, 3) the characterization of the type and phase of particulate material is important to quantifying Se exposure to prey through the base of the food web, and 4) the metric describing partitioning between particulate material and dissolved Se concentrations allows determination of a site-specific dissolved Se concentration that would be responsible for that fish body burden in the specific environment. The linked approach illustrates that environmentally safe dissolved Se concentrations will differ among ecosystems depending on the ecological pathways and biogeochemical conditions in that system. Uncertainties and model sensitivities can be directly illustrated by varying exposure scenarios based on site-specific knowledge. The model can also be used to facilitate site-specific regulation and to present generic comparisons to illustrate limitations imposed by ecosystem setting and inhabitants. Used optimally, the model provides a tool for framing a site-specific ecological problem or occurrence of Se exposure, quantify exposure within that ecosystem, and narrow uncertainties abouthowto protect it by understanding the specifics of the underlying system ecology, biogeochemistry, and hydrology.

Integrated Environmental Assessment and Management↗

Models for inference in dynamic metacommunity systems

A variety of processes are thought to be involved in the formation and dynamics of species assemblages. For example, various metacommunity theories are based on differences in the relative contributions of dispersal of species among local communities and interactions of species within local communities. Interestingly, metacommunity theories continue to be advanced without much empirical validation. Part of the problem is that statistical models used to analyze typical survey data either fail to specify ecological processes with sufficient complexity or they fail to account for errors in detection of species during sampling. In this paper, we describe a statistical modeling framework for the analysis of metacommunity dynamics that is based on the idea of adopting a unified approach, multispecies occupancy modeling, for computing inferences about individual species, local communities of species, or the entire metacommunity of species. This approach accounts for errors in detection of species during sampling and also allows different metacommunity paradigms to be specified in terms of species‐ and location‐specific probabilities of occurrence, extinction, and colonization: all of which are estimable. In addition, this approach can be used to address inference problems that arise in conservation ecology, such as predicting temporal and spatial changes in biodiversity for use in making conservation decisions. To illustrate, we estimate changes in species composition associated with the species‐specific phenologies of flight patterns of butterflies in Switzerland for the purpose of estimating regional differences in biodiversity.

Ecology↗

Application of the Hydroecological Integrity Assessment Process for Missouri Streams

Natural flow regime concepts and theories have established the justification for maintaining or restoring the range of natural hydrologic variability so that physiochemical processes, native biodiversity, and the evolutionary potential of aquatic and riparian assemblages can be sustained. A synthesis of recent research advances in hydroecology, coupled with stream classification using hydroecologically relevant indices, has produced the Hydroecological Integrity Assessment Process (HIP). HIP consists of (1) a regional classification of streams into hydrologic stream types based on flow data from long-term gaging-station records for relatively unmodified streams, (2) an identification of stream-type specific indices that address 11 subcomponents of the flow regime, (3) an ability to establish environmental flow standards, (4) an evaluation of hydrologic alteration, and (5) a capacity to conduct alternative analyses. The process starts with the identification of a hydrologic baseline (reference condition) for selected locations, uses flow data from a stream-gage network, and proceeds to classify streams into hydrologic stream types. Concurrently, the analysis identifies a set of non-redundant and ecologically relevant hydrologic indices for 11 subcomponents of flow for each stream type. Furthermore, regional hydrologic models for synthesizing flow conditions across a region and the development of flow-ecology response relations for each stream type can be added to further enhance the process. The application of HIP to Missouri streams identified five stream types ((1) intermittent, (2) perennial runoff-flashy, (3) perennial runoff-moderate baseflow, (4) perennial groundwater-stable, and (5) perennial groundwater-super stable). Two Missouri-specific computer software programs were developed: (1) a Missouri Hydrologic Assessment Tool (MOHAT) which is used to establish a hydrologic baseline, provide options for setting environmental flow standards, and compare past and proposed hydrologic alterations; and (2) a Missouri Stream Classification Tool (MOSCT) designed for placing previously unclassified streams into one of the five pre-defined stream types.

Open-File Report↗

Estimating species misclassification with occupancy dynamics and encounter rates: A semi-supervised, individual-level approach

1. Large-scale, long-term biodiversity monitoring is essential to conservation, land management, and identifying threats to biodiversity. However, multispecies surveys are prone to various types of observation error, including false positive/negative detection, and misclassification, where a species is thought to have been encountered but not correctly identified. Previous methods assume an imperfect classifier produces species-level classifications, but in practice, particularly with human observers, we may end up with extraspecific classifications including `unknown', morphospecies designations, and taxonomic identifications coarser than species. Disregarding these types of species misclassification in biodiversity monitoring datasets can bias estimates of ecologically important quantities such as demographic ratess, occurrence, and species richness. 2. Here we present a joint classification-occupancy model that accounts for species non-detection and misclassification. Our framework accommodates extinction and colonization dynamics, allows for additional uncertain `morphospecies' designations, and makes use of individual specimens with known species identities in a semi-supervised setting. We compare the performance of our model to a classification-only model that discards information about occupancy and encounter rate. We illustrate our model with an empirical case study of the carabid beetle (Carabidae) community at the National Ecological Observatory Network Niwot Ridge Mountain Research Station, near Boulder, CO, USA. We also use simulations to evaluate model performance through validation metrics where varying fractions of the data are confirmed. 3. The model supported imperfect classifier accuracy and favored certain true species classifications strongly for some morphospecies. The model outperformed (e.g., precision) the reduced model that discarded occupancy information, and these differences were most pronounced for abundant species. 4. Spatial and temporal dynamics from modeled occupancy and encounter rates may inform species misclassification probability, but this idea has not yet been tested. Our statistical framework explores this opportunity, and can be applied to datasets with imperfect species detection and classification, limited verification data, and non-species classifications.

Methods in Ecology and Evolution↗

Abundance estimation and conservation biology

Abundance is the state variable of interest in most population–level ecological research and in most programs involving management and conservation of animal populations. Abundance is the single parameter of interest in capture–recapture models for closed populations (e.g., Darroch, 1958; Otis et al., 1978; Chao, 2001). The initial capture–recapture models developed for partially (Darroch, 1959) and completely (Jolly, 1965; Seber, 1965) open populations represented efforts to relax the restrictive assumption of population closure for the purpose of estimating abundance. Subsequent emphases in capture–recapture work were on survival rate estimation in the 1970’s and 1980’s (e.g., Burnham et al., 1987; Lebreton et al.,1992), and on movement estimation in the 1990’s (Brownie et al., 1993; Schwarz et al., 1993). However, from the mid–1990’s until the present time, capture–recapture investigators have expressed a renewed interest in abundance and related parameters (Pradel, 1996; Schwarz & Arnason, 1996; Schwarz, 2001). The focus of this session was abundance, and presentations covered topics ranging from estimation of abundance and rate of change in abundance, to inferences about the demographic processes underlying changes in abundance, to occupancy as a surrogate of abundance. The plenary paper by Link & Barker (2004) is provocative and very interesting, and it contains a number of important messages and suggestions. Link & Barker (2004) emphasize that the increasing complexity of capture–recapture models has resulted in large numbers of parameters and that a challenge to ecologists is to extract ecological signals from this complexity. They offer hierarchical models as a natural approach to inference in which traditional parameters are viewed as realizations of stochastic processes. These processes are governed by hyperparameters, and the inferential approach focuses on these hyperparameters. Link & Barker (2004) also suggest that our attention should be focused on relationships between demographic processes such as survival and recruitment, the two quantities responsible for changes in abundance, rather than simply on the magnitudes of these quantities. They describe a type of Jolly–Seber capture–recapture model that permits inference about the underlying relationship between per capita recruitment rates and survival rates (Link & Barker, this volume). Implementation used Bayesian Markov Chain Monte Carlo methods and appeared to work well, yielding inferences about the relationship between recruitment and survival that were robust to selection of prior distribution. We believe that readers will find their arguments compelling, and we expect to see increased use of hierarchical modeling approaches in capture–recapture and related fields. Otto (presentation without paper) also recommended use of hierarchical models in analysis of multiple data sources dealing with population dynamics of North American mallards. He integrated survival inferences from ringing data, abundance information from aerial survey data, and recruitment information based on age ratios from a harvest survey. He used a Leslie matrix population projection model as an integrating framework and obtained estimates of breeding population size using all data.Otto’s approach also permitted inference about biases in estimated quantities. As with the work of Link & Barker (2004), we find Otto’s recommendation to use hierarchical models to integrate data from multiple sources to be very compelling. Alisauskas et al. (2004) report results of an analysis of capture–recapture data for a askatchewan population of white–winged scoters. They used the approach of Pradel (1996) to estimate population growth rate (See the PDF) directly. Estimates for 1975–1985 were quite low, but estimates for the recent period, 2000–2003,increased to values > 1. Parameter estimates for seniority, survival and per capita recruitment (Pradel, 1996) led to the inference that increased recruitment was largely responsible for the improvements in population status and growth. However, various data sources also indicated that this increase in recruitment was likely a result of increased immigration rather than improved reproduction on the area. This latter inference is important from a conservation perspective in indicating the importance of birds in other locations to growth and health of the study population. Lukacs and Burnham presented material to be published elsewhere that dealt with the use of genetic markers in capture–recapture studies. The data sources for such studies are samples of hair or feces, which are then analyzed using molecular genetic techniques in order to determine individual genotypes with respect to a usually small number of loci. Two types of classification error can arise in such analyses. First, if only a small number of loci is examined, then there may be nonnegligible probabilities that multiple individual animals will have the same genotypes. The second type of error arises during the polymerase chain reaction (PCR) process and can result from failure of alleles to amplify (allelic dropout) or from PCR inhibitors in hair and feces that produce the appearance of false alleles or misprinting (Creel et al., 2003). Lukacs and Burnham developed models that formally incorporate possible misclassification of samples resulting from these errors. These models permit estimation of parameters such as abundance and survival in a manner that properly incorporates this uncertainty of individual identity. We anticipate that noninvasive sampling based on molecular genetic analyses of hair or feces will become extremely important for some species, and that the models of Lukacs and Burnham will become very popular for such analyses. MacKenzie & Nichols (2004) discuss the use of occupancy (proportion of patches or habitat area that is occupied) as a surrogate for abundance. In cases of territorial species and where birds occur at low densities, the number of occupied patches may provide a reasonable estimate of abundance. In other cases, occupancy can be viewed as providing information about one tail of the abundance distribution, P (N = 0). The motivation for considering occupancy as a surrogate for abundance is that occupancy is based on so–called presence–absence surveys that are frequently less expensive of time and effort than methods that estimate abundance directly. We describe one set of models that can be used to estimate occupancy for a single season and another that can be used to estimate parameters such as local probabilities of extinction and colonization that are associated with occupancy dynamics. We outline a possible hybrid approach that combines occupancy data with data on marked individuals in order to betterexplore the mechanisms underlying occupancy dynamics. These five presentations made for an interesting session containing useful information and recommendations for future work. A number of themes connecting these presentations could be emphasized. For example, two of the presentations considered alternatives to standard capture–recapture sampling that can be used to draw inferences about abundance, or a portion of the abundance distribution, with field methods that should be less expensive than usual capture–recapture approaches of handling animals. We believe that the most important theme of the session was the emphasis on the processes responsible for changes in abundance. In particular, we are excited by the potential for using hierarchical models as a means of investigating relationships among vital rates and as a means of combining multiple sources of data relevant to system dynamics. Indeed, we expect the importance of this session theme to be reflected in the content and presentations of the next EURING meeting.

Animal Biodiversity and Conservation↗

Tigers and their prey: Predicting carnivore densities from prey abundance

The goal of ecology is to understand interactions that determine the distribution and abundance of organisms. In principle, ecologists should be able to identify a small number of limiting resources for a species of interest, estimate densities of these resources at different locations across the landscape, and then use these estimates to predict the density of the focal species at these locations. In practice, however, development of functional relationships between abundances of species and their resources has proven extremely difficult, and examples of such predictive ability are very rare. Ecological studies of prey requirements of tigers Panthera tigris led us to develop a simple mechanistic model for predicting tiger density as a function of prey density. We tested our model using data from a landscape-scale long-term (1995-2003) field study that estimated tiger and prey densities in 11 ecologically diverse sites across India. We used field techniques and analytical methods that specifically addressed sampling and detectability, two issues that frequently present problems in macroecological studies of animal populations. Estimated densities of ungulate prey ranged between 5.3 and 63.8 animals per km 2 . Estimated tiger densities (3.2-16.8 tigers per 100 km 2 ) were reasonably consistent with model predictions. The results provide evidence of a functional relationship between abundances of large carnivores and their prey under a wide range of ecological conditions. In addition to generating important insights into carnivore ecology and conservation, the study provides a potentially useful model for the rigorous conduct of macroecological science.

Proceedings of the National Academy of Sciences↗

Modelling heterogeneity in the recoveries of marked animal populations with covariates of individual animals, groups of animals or recovery time

A general framework is developed for modelling rates of survival and recovery of marked animal populations in terms of auxiliary information collected at the time of marking. The framework may be used to estimate differences in survival or recovery among individual animals, groups of animals, and recovery times. Analyses of the recoveries of tagged fish and banded bird populations are used to illustrate the specification and selection of various models.

Environmental and Ecological Statistics↗

The southern megalopolis: using the past to predict the future of urban sprawl in the Southeast U.S.

The future health of ecosystems is arguably as dependent on urban sprawl as it is on human-caused climatic warming. Urban sprawl strongly impacts the urban ecosystems it creates and the natural and agro-ecosystems that it displaces and fragments. Here, we project urban sprawl changes for the next 50 years for the fast-growing Southeast U.S. Previous studies have focused on modeling population density, but the urban extent is arguably as important as population density per se in terms of its ecological and conservation impacts. We develop simulations using the SLEUTH urban growth model that complement population-driven models but focus on spatial pattern and extent. To better capture the reach of low-density suburban development, we extend the capabilities of SLEUTH by incorporating street-network information. Our simulations point to a future in which the extent of urbanization in the Southeast is projected to increase by 101% to 192%. Our results highlight areas where ecosystem fragmentation is likely, and serve as a benchmark to explore the challenging tradeoffs between ecosystem health, economic growth and cultural desires.

PLoS ONE↗

Estimating habitat carrying capacity for migrating and wintering waterfowl: Considerations, pitfalls and improvements

Population-based habitat conservation planning for migrating and wintering waterfowl in North America is carried out by habitat Joint Venture (JV) initiatives and is based on the premise that food can limit demography (i.e. food limitation hypothesis). Consequently, planners use bioenergetic models to estimate food (energy) availability and population-level energy demands at appropriate spatial and temporal scales, and translate these values into regional habitat objectives. While simple in principle, there are both empirical and theoretical challenges associated with calculating energy supply and demand including: 1) estimating food availability, 2) estimating the energy content of specific foods, 3) extrapolating site-specific estimates of food availability to landscapes for focal species, 4) applicability of estimates from a single species to other species, 5) estimating resting metabolic rate, 6) estimating cost of daily behaviours, and 7) estimating costs of thermoregulation or tissue synthesis. Most models being used are daily ration models (DRMs) whose set of simplifying assumptions are well established and whose use is widely accepted and feasible given the empirical data available to populate such models. However, DRMs do not link habitat objectives to metrics of ultimate ecological importance such as individual body condition or survival, and largely only consider food-producing habitats. Agent-based models (ABMs) provide a possible alternative for creating more biologically realistic models under some conditions; however, ABMs require different types of empirical inputs, many of which have yet to be estimated for key North American waterfowl. Decisions about how JVs can best proceed with habitat conservation would benefit from the use of sensitivity analyses that could identify the empirical and theoretical uncertainties that have the greatest influence on efforts to estimate habitat carrying capacity. Development of ABMs at restricted, yet biologically relevant spatial scales, followed by comparisons of their outputs to those generated from more simplistic, deterministic models can provide a means of assessing degrees of dissimilarity in how alternative models describe desired landscape conditions for migrating and wintering waterfowl.

Wildfowl↗

Unravelling the influence of landscape alteration from flow alteration on benthic macroinvertebrate assemblage response in the Delaware River Basin

Quantifying the effects of streamflow alteration on assemblage response is central to understanding the role humans play in shaping aquatic environments. These changes represent a level of complexity that impedes developing quantitative links between flow and ecological response because stream hydrology is strongly intertwined with natural and anthropogenic factors. Better management outcomes require disentangling these linkages. Benthic macroinvertebrate data were combined with GIS-derived natural and anthropogenic basin characteristics to identify factors associated with changes in flow processes and assemblage characteristics. Models linking streamflow metrics and macroinvertebrate response at basin and subregion scales were developed using boosted regression tree (BRT) analysis. Basin-scale BRT analyses revealed that links between macroinvertebrate response and flow metrics were often obscured, whereas more homogeneous subregions were better able to discern relations with flow. Urban land cover was the primary factor accounting for changes in flow characteristics. Elevation, land cover, and high flow frequency were the principal variables driving changes in assemblage structure within subregions. Assemblage metrics and traits were equally useful for building response models and were affected similarly by streamflow alteration. Results indicate that response models should be developed based on upland and coastal subregions. However, when defining subregions, care should be taken to maintain data sufficiency. Developing practical flow-protection standards that support a balance between human water requirements and ecological integrity requires models that reduce uncertainty and identify management-relevant drivers. However, effective management often differs by location and models developed at the subregion level may be more applicable to management and stakeholder interests.

Delaware, New Jersey, New York, Pennsylvania↗

Assessing conservation relevance of organism-environment relations using predicted changes in response variables

1. Organism–environment models are used widely in conservation. The degree to which they are useful for informing conservation decisions – the conservation relevance of these relations – is important because lack of relevance may lead to misapplication of scarce conservation resources or failure to resolve important conservation dilemmas. Even when models perform well based on model fit and predictive ability, conservation relevance of associations may not be clear without also knowing the magnitude and variability of predicted changes in response variables. 2. We introduce a method for evaluating the conservation relevance of organism–environment relations that employs confidence intervals for predicted changes in response variables. The confidence intervals are compared to a preselected magnitude of change that marks a threshold (trigger) for conservation action. To demonstrate the approach, we used a case study from the Chihuahuan Desert involving relations between avian richness and broad-scale patterns of shrubland. We considered relations for three winters and two spatial extents (1- and 2-km-radius areas) and compared predicted changes in richness to three thresholds (10%, 20% and 30% change). For each threshold, we examined 48 relations. 3. The method identified seven, four and zero conservation-relevant changes in mean richness for the 10%, 20% and 30% thresholds respectively. These changes were associated with major (20%) changes in shrubland cover, mean patch size, the coefficient of variation for patch size, or edge density but not with major changes in shrubland patch density. The relative rarity of conservation-relevant changes indicated that, overall, the relations had little practical value for informing conservation decisions about avian richness. 4. The approach we illustrate is appropriate for various response and predictor variables measured at any temporal or spatial scale. The method is broadly applicable across ecological environments, conservation objectives, types of statistical predictive models and levels of biological organization. By focusing on magnitudes of change that have practical significance, and by using the span of confidence intervals to incorporate uncertainty of predicted changes, the method can be used to help improve the effectiveness of conservation efforts.

Methods in Ecology and Evolution↗

Streamflow as a stressor: Disentangling hydrology and water quality impacts to characterize flow-ecology relationships for two stream assemblages across two southeastern landscapes

Disassociating the independent effects of flow and water quality on the ecology of flowing waters is an overarching goal in water resource science needed to improve the efficacy of watershed management. However, the interrelatedness of these gradients and their subsequent alteration due to land use change has constrained progress made on this front. The objective of this study was to use benthic macroinvertebrate and fish assemblage data to characterize flow-ecology relationships that were unchanged by water quality impacts across two southeastern landscapes in the USA to help detect ecological change driven by flow alteration. General linear latent models were used to identify taxa that were responsive to high or low flow metrics and water quality gradients. Bayesian hierarchical generalized additive models were then developed using these indicator taxa and three biological metrics to identify flow-specific relationships that were unaffected by water quality impacts. Three low flow-specific relationships were identified, illustrating how potential agricultural or urban impacts to hydrology reduced stream biological health. Importantly, flow-ecology relationships developed using indicator taxa in this study effectively captured hydrology-specific impacts while biological metrics typical of state monitoring and assessment programs did not. Therefore, developing flow-specific biological metrics is a critical step when developing management strategies targeting flow alteration. Implementing standardized frameworks such as the one characterized here can limit contradictory findings and improve streamflow enhancement and restoration project efficacy. These low flow-specific relationships will enhance managers' capacity to develop environmental flow standards, monitor their success, and better understand urban and agricultural impacts on stream assemblages.

Alabama↗

Tools for quantifying isotopic niche space and dietary variation at the individual and population level.

Ecologists are increasingly using stable isotope analysis to inform questions about variation in resource and habitat use from the individual to community level. In this study we investigate data sets from 2 California sea otter ( Enhydra lutris nereis ) populations to illustrate the advantages and potential pitfalls of applying various statistical and quantitative approaches to isotopic data. We have subdivided these tools, or metrics, into 3 categories: IsoSpace metrics, stable isotope mixing models, and DietSpace metrics. IsoSpace metrics are used to quantify the spatial attributes of isotopic data that are typically presented in bivariate (e.g., δ 13 C versus δ 15 N) 2-dimensional space. We review IsoSpace metrics currently in use and present a technique by which uncertainty can be included to calculate the convex hull area of consumers or prey, or both. We then apply a Bayesian-based mixing model to quantify the proportion of potential dietary sources to the diet of each sea otter population and compare this to observational foraging data. Finally, we assess individual dietary specialization by comparing a previously published technique, variance components analysis, to 2 novel DietSpace metrics that are based on mixing model output. As the use of stable isotope analysis in ecology continues to grow, the field will need a set of quantitative tools for assessing isotopic variance at the individual to community level. Along with recent advances in Bayesian-based mixing models, we hope that the IsoSpace and DietSpace metrics described here will provide another set of interpretive tools for ecologists.

Journal of Mammalogy↗

Cohesive framework for modeling plant cover class data

The study of plant distribution and abundance is a fundamental pursuit in ecology and conservation biology. Measuring plant abundance by visually assessing percent cover and recording a cover class is a common field method that yields ordinal data. Statistical models for ordinal data exist but entail cumbersome interpretations and sometimes restrictive assumptions. We propose a Bayesian hierarchical framework for analysing cover class data that allows for linking ordinal observations to a latent beta distribution and accounts for zero inflation. Harnessing a latent beta distribution supports interpreting changes in abundance in terms of mean percent cover rather than odds ratios of cumulative cover classes as for cumulative link models. The zero augmentation allows for simultaneous inferences on both occurrence (distribution) and abundance. We show how our model can account for true and false zeros, misclassification of cover classes, multiple species and hierarchical sampling designs, using empirical examples and simulations. Simulated observation errors, when ignored, led to models overestimating abundance and underestimating occurrence. Based on simulations, we found no substantial difference between mean percent cover estimates when analyzing ordinal cover classes versus continuous percent cover as the response. Our empirical datasets displayed high probability of detection (>0.85 on average for all species), likely due to the sampling design used and training of observers. Probability of occurrence was slightly underestimated for bare ground, Artemisia tridentata , Elycap medusae , and Poa secunda using a model that ignored imperfect detection. Estimated mean percent cover was not substantially impacted by ignoring measurement error for five plant species and bare ground. Our modelling framework for cover class data allows for an explicit separation of distribution from abundance and, importantly, allows for interpreting species–environment relationships in terms of variation in mean percent cover as compared to cumulative odds ratios. The beta distribution inherently accommodates heteroscedasticity and skewness, statistical properties that are a consequence of spatially aggregated patterns common to plant survey data. Recording cover classes provides a reliable, efficient way to measure plants and our simulations suggest little loss of information compared to assuming continuous percent cover. We provide JAGS and Stan model code for implementation.

Methods in Ecology and Evolution↗

Modeling habitat dynamics accounting for possible misclassification

Land cover data are widely used in ecology as land cover change is a major component of changes affecting ecological systems. Landscape change estimates are characterized by classification errors. Researchers have used error matrices to adjust estimates of areal extent, but estimation of land cover change is more difficult and more challenging, with error in classification being confused with change. We modeled land cover dynamics for a discrete set of habitat states. The approach accounts for state uncertainty to produce unbiased estimates of habitat transition probabilities using ground information to inform error rates. We consider the case when true and observed habitat states are available for the same geographic unit (pixel) and when true and observed states are obtained at one level of resolution, but transition probabilities estimated at a different level of resolution (aggregations of pixels). Simulation results showed a strong bias when estimating transition probabilities if misclassification was not accounted for. Scaling-up does not necessarily decrease the bias and can even increase it. Analyses of land cover data in the Southeast region of the USA showed that land change patterns appeared distorted if misclassification was not accounted for: rate of habitat turnover was artificially increased and habitat composition appeared more homogeneous. Not properly accounting for land cover misclassification can produce misleading inferences about habitat state and dynamics and also misleading predictions about species distributions based on habitat. Our models that explicitly account for state uncertainty should be useful in obtaining more accurate inferences about change from data that include errors.

Landscape Ecology↗

Compilation of geologic, hydrologic, and ground-water flow modeling information for the Spokane Valley-Rathdrum Prairie aquifer, Spokane County, Washington, and Bonner and Kootenai Counties, Idaho

The U.S. Geological Survey, in cooperation with the Idaho Department of Water Resources and Washington Department of Ecology compiled and described geologic, hydrologic, and ground-water flow modeling information about the Spokane Valley-Rathdrum Prairie (SVRP) aquifer in northern Idaho and northeastern Washington. Descriptions of the hydrogeologic framework, water-budget components, ground- and surface-water interactions, computer flow models, and further data needs are provided. The SVRP aquifer, which covers about 370 square miles including the Rathdrum Prairie, Idaho and the Spokane valley and Hillyard Trough, Washington, was designated a Sole Source Aquifer by the U.S. Environmental Protection Agency in 1978. Continued growth, water management issues, and potential effects on water availability and water quality in the aquifer and in the Spokane and Little Spokane Rivers have illustrated the need to better understand and manage the region's water resources. The SVRP aquifer is composed of sand, gravel, cobbles, and boulders primarily deposited by a series of catastrophic glacial outburst floods from ancient Glacial Lake Missoula. The material deposited in this high-energy environment is coarser-grained than is typical for most basin-fill deposits, resulting in an unusually productive aquifer with well yields as high as 40,000 gallons per minute. In most places, the aquifer is bounded laterally by bedrock composed of granite, metasedimentary rocks, or basalt. The lower boundary of the aquifer is largely unknown except along the margins or in shallower parts of the aquifer where wells have penetrated its entire thickness and reached bedrock or silt and clay deposits. Based on surface geophysics, the thickness of the aquifer is about 500 ft near the Washington-Idaho state line, but more than 600 feet within the Rathdrum Prairie and more than 700 feet in the Hillyard trough based on drilling records. Depth to water in the aquifer is greatest in the northern Rathdrum Prairie (about 500 feet) and least near the city of Spokane along the Spokane River (less than about 50 feet). Ground-water flow is south from near the southern end of Lake Pend Oreille and Hoodoo Valley, through the Rathdrum Prairie, then west toward Spokane. In Spokane, the aquifer splits and water moves north through the Hillyard Trough as well as west through the Trinity Trough. From the Trinity Trough water flows north along the western arm of the aquifer. The aquifer's discharge area is along the Little Spokane River and near Long Lake, Washington. A compilation of estimates of water-budget components, including recharge (precipitation, irrigation, canal leakage, septic tank effluent, inflow from tributary basins, and flow from the Spokane River) and discharge (withdrawals from wells, flow to the Spokane and Little Spokane Rivers, evapotranspiration, and underflow to Long Lake) illustrates that these estimated values should be compared with caution due to several variables including the area and time period of interest as well as methods employed in making the estimates. Numerous studies have documented the dynamic ground-water and surface-water interaction between the SVRP aquifer and the Spokane and Little Spokane Rivers. Gains and losses vary throughout the year, as well as the locations of gains and losses. September 2004 streamflow measurements indicated that the upper reach of the Spokane River between Post Falls and downstream at Flora Road lost 321 cubic feet per second. A gain of 736 cubic feet per second was measured between the Flora Road site and downstream at Green Street Bridge. A loss of 124 cubic feet per second was measured for the reach between the Green Street Bridge and the Spokane River at Spokane gaging station. The river gained about 87 cubic feet per second between the Spokane River at Spokane gaging station and the TJ Meenach Bridge. Overall, the Spokane River gained about 284 cubic feet per second between the Post Falls,

Scientific Investigations Report↗

Estimated probability of arsenic in groundwater from bedrock aquifers in New Hampshire, 2011

Probabilities of arsenic occurrence in groundwater from bedrock aquifers at concentrations of 1, 5, and 10 micrograms per liter (µg/L) were estimated during 2011 using multivariate logistic regression. These estimates were developed for use by the New Hampshire Environmental Public Health Tracking Program. About 39 percent of New Hampshire bedrock groundwater was identified as having at least a 50 percent chance of containing an arsenic concentration greater than or equal to 1 µg/L. This compares to about 7 percent of New Hampshire bedrock groundwater having at least a 50 percent chance of containing an arsenic concentration equaling or exceeding 5 µg/L and about 5 percent of the State having at least a 50 percent chance for its bedrock groundwater to contain concentrations at or above 10 µg/L. The southeastern counties of Merrimack, Strafford, Hillsborough, and Rockingham have the greatest potential for having arsenic concentrations above 5 and 10 µg/L in bedrock groundwater. Significant predictors of arsenic in groundwater from bedrock aquifers for all three thresholds analyzed included geologic, geochemical, land use, hydrologic, topographic, and demographic factors. Among the three thresholds evaluated, there were some differences in explanatory variables, but many variables were the same. More than 250 individual predictor variables were assembled for this study and tested as potential predictor variables for the models. More than 1,700 individual measurements of arsenic concentration from a combination of public and private water-supply wells served as the dependent (or predicted) variable in the models. The statewide maps generated by the probability models are not designed to predict arsenic concentration in any single well, but they are expected to provide useful information in areas of the State that currently contain little to no data on arsenic concentration. They also may aid in resource decision making, in determining potential risk for private wells, and in ecological-level analysis of disease outcomes. The approach for modeling arsenic in groundwater could also be applied to other environmental contaminants that have potential implications for human health, such as uranium, radon, fluoride, manganese, volatile organic compounds, nitrate, and bacteria.

New Hampshire↗