Geology ReportsSearch

SEARCH · Geology Reports

Results for “Spatial Statistics”

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 505 records · Page 28Linked to original sources

Scale-dependent approaches to modeling spatial epidemiology of chronic wasting disease.

This e-book is the product of a second workshop that was funded and promoted by the United States Geological Survey to enhance cooperation between states for the management of chronic wasting disease (CWD). The first workshop addressed issues surrounding the statistical design and collection of surveillance data for CWD. The second workshop, from which this document arose, followed logically from the first workshop and focused on appropriate methods for analysis, interpretation, and use of CWD surveillance and related epidemiology data. Consequently, the emphasis of this e-book is on modeling approaches to describe and gain insight of the spatial epidemiology of CWD. We designed this e-book for wildlife managers and biologists who are responsible for the surveillance of CWD in their state or agency. We chose spatial methods that are popular or common in the spatial epidemiology literature and evaluated them for their relevance to modeling CWD. Our opinion of the usefulness and relevance of each method was based on the type of field data commonly collected as part of CWD surveillance programs and what we know about CWD biology, ecology, and epidemiology. Specifically, we expected the field data to consist primarily of the infection status of a harvested or culled sample along with its date of collection (not date of infection), location, and demographic status. We evaluated methods in light of the fact that CWD does not appear to spread rapidly through wild populations, relative to more highly contagious viruses, and can be spread directly from animal to animal or indirectly through environmental contamination. We discovered that many of the wellpublished methods were developed for fast-spreading human diseases, such as influenza and measles. While these methods are applicable to fast spreading wildlife diseases, such as foot-and-mouth disease or West Nile virus, many are not likely to work well for CWD. Only limited data exist to evaluate geographic and spatial spread because many locations where we find CWD tend to be locations where samples have just been taken or sample sizes have just become large enough to have a high probability of detecting a low prevalence. Consequently, methods that work well to describe or predict the spread of foot-and-mouth disease throughout England, which occurred within a year, do not work well for describing or predicting CWD spread. We did not exclude methods that we regarded as inappropriate; rather, we included methods that are commonly used for disease epidemiology and then discussed their applicability for modeling the spatial epidemiology of CWD. We hope including inappropriate methods with an explanation of why they are ill-suited for CWD will make it easier to drop them from consideration and explain to others why they were not recommended for spatial modeling of CWD. We organized the three chapters by scale and extent for which each method was developed or best suited. The first chapter covers methods appropriate to multi-jurisdictional or multi-state modeling, which we call “regional” scale. The second chapter covers methods appropriate for within state areas such as wildlife management units or metapopulations, which we call “landscape” scale. The third chapter covers methods appropriate for population or individual-based modeling, which we call “fine” scale. We know this rubric is somewhat artificial because many methods work at multiple scales. We hope, however, that this structure addresses some of the challenges faced by managers that work at local, regional, state, and national scales. Further, the resolution of empirical data often changes with spatial scale, which affects the utility of different modeling approaches. For example, individual-based models work best at modeling spread within populations, while risk analysis is most useful for summarizing data over larger scales such as a region. Because some methods are applicable at several scales, however, we included a graphic at the beginning of each method that indicates the range of scales for which it applies. For example, the graphic to the right indicates that the method is most applicable for regional-scale modeling. There is also a question of resolution as well as scale and extent for each method. CWD surveillance data have been collected over large areas, such as a wildlife management unit or state, but the resolution of the data may be fine scale with GPS locations for many samples. For each method, we described the required resolution of the data and describe the type of data required, as well as what questions the method could answer and how useful the method is, given typical CWD data. For each scale, we presented a focal approach that would be useful for understanding the spatial pattern and epidemiology of CWD, as well as being a useful tool for CWD management. The focal approaches include risk analysis and micromaps for the regional scale, cluster analysis for the landscape scale, and individual based modeling for the fine scale of within population. For each of these methods, we used simulated data and walked through the method step by step to fully illustrate the “how to”, with specifics about what is input and output, as well as what questions the method addresses. We also provided a summary table to, at a glance, describe the scale, questions that can be addressed, and general data required for each method described in this e-book. We hope that this review will be helpful to biologists and managers by increasing the utility of their surveillance data, and ultimately be useful for increasing our understanding of CWD and allowing wildlife biologists and managers to move beyond retroactive fire-fighting to proactive preventative action.

Book

Millennial-scale variability in the local radiocarbon reservoir age of south Florida during the Holocene

A growing body of research suggests that the marine environments of south Florida provide a critical link between the tropical and high-latitude Atlantic. Changes in the characteristics of water masses off south Florida may therefore have important implications for our understanding of climatic and oceanographic variability over a broad spatial scale; however, the sources of variability within this oceanic corridor remain poorly understood. Measurements of ΔR, the local offset of the radiocarbon reservoir age, from shallow-water marine environments can serve as a powerful tracer of water-mass sources that can be used to reconstruct variability in local-to regional-scale oceanography and hydrology. We combined radiocarbon and U-series measurements of Holocene-aged corals from the shallow-water environments of the Florida Keys reef tract (FKRT) with robust statistical modeling to quantify the millennial-scale variability in ΔR at locations with (“nearshore”) and without (“open ocean”) substantial terrestrial influence. Our reconstructions demonstrate that there was significant spatial and temporal variability in ΔR on the FKRT during the Holocene. Whereas ΔR was similar throughout the region after ∼4000 years ago, nearshore ΔR was significantly higher than in the open ocean during the middle Holocene. We suggest that the elevated nearshore ΔR from ∼8000 to 5000 years ago was most likely the result of greater groundwater influence associated with lower sea level at this time. In the open ocean, which would have been isolated from the influence of groundwater, ΔR was lowest ∼7000 years ago, and was highest ∼3000 years ago. We evaluated our open-ocean model of ΔR variability against records of local-to regional-scale oceanography and conclude that local upwelling was not a significant driver of open-ocean radiocarbon variability in this region. Instead, the millennial-scale trends in open-ocean ΔR were more likely a result of broader-scale changes in western Atlantic circulation associated with an increase in the supply of equatorial South Atlantic water to the Caribbean and shifts in the character of South Atlantic waters resulting from variation in the intensity of upwelling off the southwest coast of Africa. Because accurate estimates of ΔR are critical to precise calibrations of radiocarbon dates from marine samples, we also developed models of nearshore and open-ocean ΔR versus conventional 14 C ages that can be used for regional radiocarbon calibrations for the Holocene. Our study provides new insights into the patterns and drivers of oceanographic and hydrologic variability in the Straits of Florida and highlights the value of the paleoceanographic records from south Florida to our understanding of Holocene changes in climate and ocean circulation throughout the Atlantic.

Florida

Revealing the extent of sea otter impacts on bivalve prey through multi-trophic monitoring and mechanistic models

Sea otters are apex predators that can exert considerable influence over the nearshore communities they occupy. Since facing near extinction in the early 1900s, sea otters are making a remarkable recovery in Southeast Alaska, particularly in Glacier Bay, the largest protected tidewater glacier fjord in the world. The expansion of sea otters across Glacier Bay offers both a challenge to monitoring and stewardship and an unprecedented opportunity to study the top-down effect of a novel apex predator across a diverse and productive ecosystem. Our goal was to integrate monitoring data across trophic levels, space, and time to quantify and map the predator–prey interaction between sea otters and butter clams Saxidomus gigantea , one of the dominant large bivalves in Glacier Bay and a favoured prey of sea otters. We developed a spatially-referenced mechanistic differential equation model of butter clam dynamics that combined both environmental drivers of local population growth and estimates of otter abundance from aerial survey data. We embedded this model in a Bayesian statistical framework and fit it to clam survey data from 43 intertidal and subtidal sites across Glacier Bay. Prior to substantial sea otter expansion, we found that butter clam density was structured by an environmental gradient driven by distance from glacier (represented by latitude) and a quadratic effect of current speed. Estimates of sea otter attack rate revealed spatial heterogeneity in sea otter impacts and a negative relationship with local shoreline complexity. Sea otter exploitation of productive butter clam habitat substantially reduced the abundance and altered the distribution of butter clams across Glacier Bay, with potential cascading consequences for nearshore community structure and function. Spatial variation in estimated sea otter predation processes further suggests that community context and local environmental conditions mediate the top-down influence of sea otters on a given prey. Overall, our framework provides high-resolution insights about the interaction among components of this food web and could be applied to a variety of other systems involving invasive species, epidemiology or migration.

Journal of Animal Ecology

Regression equations for estimation of annual peak-streamflow frequency for undeveloped watersheds in Texas using an L-moment-based, PRESS-minimized, residual-adjusted approach

Annual peak-streamflow frequency estimates are needed for flood-plain management; for objective assessment of flood risk; for cost-effective design of dams, levees, and other flood-control structures; and for design of roads, bridges, and culverts. Annual peak-streamflow frequency represents the peak streamflow for nine recurrence intervals of 2, 5, 10, 25, 50, 100, 200, 250, and 500 years. Common methods for estimation of peak-streamflow frequency for ungaged or unmonitored watersheds are regression equations for each recurrence interval developed for one or more regions; such regional equations are the subject of this report. The method is based on analysis of annual peak-streamflow data from U.S. Geological Survey streamflow-gaging stations (stations). Beginning in 2007, the U.S. Geological Survey, in cooperation with the Texas Department of Transportation and in partnership with Texas Tech University, began a 3-year investigation concerning the development of regional equations to estimate annual peak-streamflow frequency for undeveloped watersheds in Texas. The investigation focuses primarily on 638 stations with 8 or more years of data from undeveloped watersheds and other criteria. The general approach is explicitly limited to the use of L-moment statistics, which are used in conjunction with a technique of multi-linear regression referred to as PRESS minimization. The approach used to develop the regional equations, which was refined during the investigation, is referred to as the 'L-moment-based, PRESS-minimized, residual-adjusted approach'. For the approach, seven unique distributions are fit to the sample L-moments of the data for each of 638 stations and trimmed means of the seven results of the distributions for each recurrence interval are used to define the station specific, peak-streamflow frequency. As a first iteration of regression, nine weighted-least-squares, PRESS-minimized, multi-linear regression equations are computed using the watershed characteristics of drainage area, dimensionless main-channel slope, and mean annual precipitation. The residuals of the nine equations are spatially mapped, and residuals for the 10-year recurrence interval are selected for generalization to 1-degree latitude and longitude quadrangles. The generalized residual is referred to as the OmegaEM parameter and represents a generalized terrain and climate index that expresses peak-streamflow potential not otherwise represented in the three watershed characteristics. The OmegaEM parameter was assigned to each station, and using OmegaEM, nine additional regression equations are computed. Because of favorable diagnostics, the OmegaEM equations are expected to be generally reliable estimators of peak-streamflow frequency for undeveloped and ungaged stream locations in Texas. The mean residual standard error, adjusted R-squared, and percentage reduction of PRESS by use of OmegaEM are 0.30log 10 , 0.86, and -21 percent, respectively. Inclusion of the OmegaEM parameter provides a substantial reduction in the PRESS statistic of the regression equations and removes considerable spatial dependency in regression residuals. Although the OmegaEM parameter requires interpretation on the part of analysts and the potential exists that different analysts could estimate different values for a given watershed, the authors suggest that typical uncertainty in the OmegaEM estimate might be about +or-0.10 10 . Finally, given the two ensembles of equations reported herein and those in previous reports, hydrologic design engineers and other analysts have several different methods, which represent different analytical tracks, to make comparisons of peak-streamflow frequency estimates for ungaged watersheds in the study area.

Texas

A comparison of chlorophyll a values obtained from an autonomous underwater vehicle to satellite-based measures for Lake Michigan

Accurate methods to track changes in lake productivity through time and space are critical to fisheries management. Chlorophyll a is the most widely studied proxy for ecosystem primary production, and has been the topic of many studies. The main sources of chlorophyll a measurements are ship-based measures or multi-spectral satellite data. Autonomous underwater vehicles can survey large spatial extents approaching the scale of satellite data, but with the accuracy of ship-based water sampling methods. We use several statistical measures to compare measures of chlorophyll a collected in Lake Michigan with spatiotemporally matched satellite-derived measures of chlorophyll a from the MODIS Aqua multi-spectral sensor using NASA’s OC3 and the Great Lakes Fit algorithms. Our findings show a near one to one relationship between AUV data and both satellite-derived data sets when the AUV data are coarsened to the resolution of the satellite data. A comparison of satellite-based chlorophyll a to AUV-derived chlorophyll summarized in discrete water depth bins suggested that, based on decreasing coefficients of determination, satellite estimates of chlorophyll accounted for the most variability in chlorophyll a concentrations in the upper 10 m of the water column, even though satellite sensors may detect past this depth.

Lake Michigan

Effect of different sampling schemes on the spatial placement of conservation reserves in Utah, USA

We evaluated the effect of three different sampling schemes used to organize spatially explicit biological information had on the spatial placement of conservation reserves in Utah, USA. The three sampling schemes consisted of a hexagon representation developed by the EPA/EMAP program (statistical basis), watershed boundaries (ecological), and the current county boundaries of Utah (socio-political). Four decision criteria were used to estimate effects, including amount of area, length of edge, lowest number of contiguous reserves, and greatest number of terrestrial vertebrate species covered. A fifth evaluation criterion was the effect each sampling scheme had on the ability of the modeled conservation reserves to cover the six major ecoregions found in Utah. Of the three sampling schemes, county boundaries covered the greatest number of species, but also created the longest length of edge and greatest number of reserves. Watersheds maximized species coverage using the least amount of area. Hexagons and watersheds provide the least amount of edge and fewest number of reserves. Although there were differences in area, edge and number of reserves among the sampling schemes, all three schemes covered all the major ecoregions in Utah and their inclusive biodiversity. ?? 2003 Elsevier Science Ltd. All rights reserved.

Biological Conservation

Meta-analysis of a megafish: Assessing patterns and predictors of Alligator Gar movement across multiple populations

Background Freshwater megafishes are among some of the most commercially and ecologically important aquatic organisms yet are disproportionately threatened with range and population reduction. Anthropogenic alterations of rivers influencing migrations are among the most significant causes for these declines. However, migratory fishes do not always respond similarly to movement barriers and thus it is necessary to develop models to predict movements of freshwater migratory fishes in the face of anthropogenic alteration. Predicting movement of freshwater fishes is often investigated using statistical packages. However, empirical studies assessing these packages have led to mixed results, questioning its applicability to all taxa. We argue that spatial, temporal, and environmental attributes are more influential for movement of a migratory megafish, the Alligator Gar ( Atractosteus spatula ), than the current parameters explored in a globally relevant fish movement model. Methods This study explored two independent mobile telemetry datasets investigating Alligator Gar movement on the Brazos and Trinity rivers. Environmental associations were investigated to predict Alligator Gar displacement and dispersal using generalized additive models, generalized linear models, and model selection. Leptokurtosis of Alligator Gar populations was also assessed. Predictability of the movement model was tested by comparing observed to model derived stationary and mobile components making up a leptokurtic movement distribution. Results Our study suggests that current and antecedent measures of discharge and water temperature are positively correlated with Alligator Gar displacement and dispersal. However, these patterns are only detectable when monthly relocation intervals are explored rather than seasonal scales. Leptokurtosis was observed in both Alligator Gar populations. However, movement was normally distributed (i.e., mesokurtic) under tracking events following high flood pulses. Additionally, predicted Alligator Gar movement was significantly farther under modeled values compared to observed values, in part because the species undergoes cyclical migrations for reproduction that are sensitive to water temperature and discharge. Discussion In conclusion, this study provides an alternative framework to assess the movement patterns of migratory fishes, which could be tested on additional freshwater fishes, and suggests that assessing spatial, environmental, and temporal processes simultaneously are necessary to capture the complexities of fish movement which currently are unavailable for the movement model we investigated.

Texas

Comparison of statistical and theoretical habitat models for conservation planning: The benefit of ensemble prediction

Selection of a modeling approach is an important step in the conservation planning process, but little guidance is available. We compared two statistical and three theoretical habitat modeling approaches representing those currently being used for avian conservation planning at landscape and regional scales: hierarchical spatial count (HSC), classification and regression tree (CRT), habitat suitability index (HSI), forest structure database (FS), and habitat association database (HA). We focused our comparison on models for five priority forest-breeding species in the Central Hardwoods Bird Conservation Region: Acadian Flycatcher, Cerulean Warbler, Prairie Warbler, Red-headed Woodpecker, and Worm-eating Warbler. Lacking complete knowledge on the distribution and abundance of each species with which we could illuminate differences between approaches and provide strong grounds for recommending one approach over another, we used two approaches to compare models: rank correlations among model outputs and comparison of spatial correspondence. In general, rank correlations were significantly positive among models for each species, indicating general agreement among the models. Worm-eating Warblers had the highest pairwise correlations, all of which were significant ( P , 0.05). Red-headed Woodpeckers had the lowest agreement among models, suggesting greater uncertainty in the relative conservation value of areas within the region. We assessed model uncertainty by mapping the spatial congruence in priorities (i.e., top ranks) resulting from each model for each species and calculating the coefficient of variation across model ranks for each location. This allowed identification of areas more likely to be good targets of conservation effort for a species, those areas that were least likely, and those in between where uncertainty is higher and thus conservation action incorporates more risk. Based on our results, models developed independently for the same purpose (conservation planning for a particular species in a particular geography) yield different answers and thus different conservation strategies. We assert that using only one habitat model (even if validated) as the foundation of a conservation plan is risky. Using multiple models (i.e., ensemble prediction) can reduce uncertainty and increase efficacy of conservation action when models corroborate one another and increase understanding of the system when they do not.

Ecological Applications

Preliminary data used to assess the accuracy of estimating female white-tailed deer diel birthing-season home ranges using only daytime locations

Because many white-tailed deer (Odocoileus virginianus) home-range and habitat-use studies rely only on daytime radio-tracking data, we were interested in whether diurnal data sufficiently represented diel home ranges. We analyzed home-range and core-use size and overlap of 8 adult-female Global-Positioning-System-collared deer during May and June 2001 and 2002 in the Superior National Forest, Minnesota, USA. We used 2 traditional means of analysis: minimum-convex polygons (MCP) and fixed kernels (95% FK, home range and 50% FK, core use) and two methods to partition day and night location data: (1) daytime = 0800-2000 h versus nighttime = 2000-0800 h and (2) sunup versus sundown. We found no statistical difference in size of home-range and core-use areas across day and night comparisons; however, in terms of spatial overlap, approximately 30% of night-range areas on average were not accounted for using daytime locations, with even greater differences between core-use areas (on average approximately 50%). We conclude that diurnal data do not adequately describe diel adult-female-deer, May-June home-ranges due to differences in spatial overlap (location). We suggest research to determine (1) if our findings hold in other circumstances (e.g., exclusive of the parturition period, other age classes, etc.), (2) if our conclusions generalize under other conditions (e.g., across deer range, varying seasons, etc.), (3) if habitat-use conclusions are affected by the incomplete overlap between diurnal and diel data, (4) how many nocturnal locations must be included to generate sufficient overlap, and (5) the influence of using other kernel sizes (e.g., 75%, 90%).

Wildlife Biology in Practice

Data on selected herbicides and two triazine metabolites in precipitation of the Midwestern and Northeastern United States, 1990-91

Weekly precipitation (rain and snow) samples were collected from 81 National Atmospheric Deposition Program/National Trends Network sites in the Midwestern and Northeastern United States for the analysis of herbicides. In addition, five high- elevation background sites along the Rocky Mountains and in Alaska were sampled to provide data on herbicides in precipitation at sites far from the study area. The study began in March 1990 and continued through September 1991. The precipitation samples were shipped to the Central Analytical Laboratory operated by the Illinois State Water Survey for analyses of inorganic compounds. Subsamples of the precipitation were shipped to the U.S. Geological Survey laboratory in Lawrence, Kansas, for the analysis of 11 herbicides and 2 triazines metabolites. This report provides descriptions of the study area, sample-collection methods, laboratory methods, and quality assurance. The report also includes a compilation of herbicide concentration data from both enzyme-linked immunosorbent assay and gas chromatography/mass spectrometry methods. Laboratory analyses consisted of 6,230 samples that were analyzed by enzyme-linked immunosorbent assay using alachlor and atrazine microtiter plates and 2,341 samples that were confirmed by gas chromatography/mass spectrometry analysis. Graphical and statistical comparisons of the two analytical methods are given in this report. Data from this study have been useful in determining the spatial and temporal distribution of herbicide concentrations and deposition in precipitation of over a 26-State area of the Midwestern and Northeastern United States. The data also provide evidence of long-range atmospheric transport of herbicides and triazine metabolites.

Open-File Report

Changes in water quality of Michigan streams near urban areas, 1973-84

Monthly water-quality monitoring of streams was begun by Michigan Department of Natural Resources in 1973 to (1) determine temporal and spatial variability, (2) detect long-term trends, and (3) describe changes in water quality near urban areas. This report provides a statistical analysis and summary of data collected from 1973 through 1984. Concentrations and discharges of nine commonly measured water-quality constituents and specific conductance are examined. Twenty-three sites on inland streams (streams draining basins wholly within Michigan) and 20 sites on Detroit River are discussed. The changes in water quality in 9 rivers near 12 urban areas in Michigan's southern Lower Peninsula and the relation between streamflow and selected waterquality characteristics, including phosphorus, chloride, sulfate, nitrogen, specific conductance, and solid residues are described. Results show that the median dissolved-solids concentration in Clinton River downstream from Pontiac exceed Michigan's 1986 stream water-quality standard. Among inland streams, constituent concentrations and discharges generally were greatest in Saginaw River and least in Grand River upstream from Jackson. Upstream from Detroit, constituent concentrations in Detroit River did not differ appreciably across the Windmill Point Transect; downstream from Detroit, at the Fermi Transect across Detroit River, most constituent concentrations were higher near the American and Canadian shorelines. Among urban areas, greatest changes in constituent concentrations occurred in the Grand River near Jackson, in the Clinton River near Pontiac, and in the Tittabawassee River near Midland; the least changes in constituent concentrations occurred in the Saginaw River near Saginaw, in Detroit River near Detroit and the Kalamazoo River near Battle Creek. Greatest changes in constituent discharges occurred in the Detroit River near the Detroit area; the least occurred in the Chippewa River near Mount Pleasant. Of the 230 regressions between streamflow and constituent concentrations in inland streams, about 73 percent were significant at the 5-percent level. The degree of the correlation and nature of the relation varied among sites and constituents. Generally, higher streamflows were associated with lower concentrations. Changes in streamflow and changes in constituent concentrations near urban areas were correlated in 57 percent of the 120 analyses. Generally, higher changes in streamflow were associated with lower changes in concentrations.

Michigan

Ecological diffusion models are still useful in ecology

Context Connecting animal movement across scales—both among spatial scales and between individual and population level processes—is a central theme in ecology, including the sub-fields of movement and landscape ecology. However, most modeling frameworks operate at either the individual (Lagrangian) or population (Eulerian) level. Objectives We aimed to review how a partial differential equation (PDE) known as the ecological diffusion equation (EDE) offers connections among the scales of animal movement and how it can be embedded in commonly used statistical frameworks (e.g., hierarchical models) and fit to multiple types of data, such as count, presence-only, presence-absence, and telemetry data. In doing so, we also highlight emerging and promising avenues of research related to statistical EDE models. Methods We present how the EDE is derived from a simple set of first principles, can be used to represent both individual animal movement and spatiotemporal population dynamics, and can be fit to data by being applied in a hierarchical statistical framework. Finally, we review literature on recent applications and advances related to EDE models. Results In presenting the derivation of the EDE and reviewing recent advances and applications of EDE models, we established that the EDE can enable a mechanistic understanding of animal movement across scales, as well as the integration of multiple types of data. We also identified a number of topics where EDE models can be advanced to gain additional ecological insight. Conclusions Ecological diffusion is a process that emerges naturally from first principles governing animal movement and can be used to model both individual and population-level spatiotemporal dynamics. Statistical EDE models are related to other modeling frameworks, such as species distribution models, occupancy models, and step-selection analysis. The flexibility of hierarchical statistical modeling allows combining different types of data collected on movement and population dynamics, and the theoretical underpinnings of the EDE allow it to connect ecological processes occurring at fine and coarse spatial scales, as well as individual animal movement and population dynamics.

Landscape Ecology

Sampling design for spatially distributed hydrogeologic and environmental processes

A methodology for the design of sampling networks over space is proposed. The methodology is based on spatial random field representations of nonhomogeneous natural processes, and on optimal spatial estimation techniques. One of the most important results of random field theory for physical sciences is its rationalization of correlations in spatial variability of natural processes. This correlation is extremely important both for interpreting spatially distributed observations and for predictive performance. The extent of site sampling and the types of data to be collected will depend on the relationship of subsurface variability to predictive uncertainty. While hypothesis formulation and initial identification of spatial variability characteristics are based on scientific understanding (such as knowledge of the physics of the underlying phenomena, geological interpretations, intuition and experience), the support offered by field data is statistically modelled. This model is not limited by the geometric nature of sampling and covers a wide range in subsurface uncertainties. A factorization scheme of the sampling error variance is derived, which possesses certain atttactive properties allowing significant savings in computations. By means of this scheme, a practical sampling design procedure providing suitable indices of the sampling error variance is established. These indices can be used by way of multiobjective decision criteria to obtain the best sampling strategy. Neither the actual implementation of the in-situ sampling nor the solution of the large spatial estimation systems of equations are necessary. The required values of the accuracy parameters involved in the network design are derived using reference charts (readily available for various combinations of data configurations and spatial variability parameters) and certain simple yet accurate analytical formulas. Insight is gained by applying the proposed sampling procedure to realistic examples related to sampling problems in two dimensions.

Advances in Water Resources

Occupancy models for citizen-science data

Large‐scale citizen‐science projects, such as atlases of species distribution, are an important source of data for macroecological research, for understanding the effects of climate change and other drivers on biodiversity, and for more applied conservation tasks, such as early‐warning systems for biodiversity loss. However, citizen‐science data are challenging to analyse because the observation process has to be taken into account. Typically, the observation process leads to heterogeneous and non‐random sampling, false absences, false detections, and spatial correlations in the data. Increasingly, occupancy models are being used to analyse atlas data. We advocate a dual approach to strengthen inference from citizen science data for the questions the programme is intended to address: (a) the survey design should be chosen with a particular set of questions and associated analysis strategy in mind and (b) the statistical methods should be tailored not only to those questions but also to the specific characteristics of the data. We review the consequences of particular survey design choices that typically need to be made in atlas‐style citizen‐science projects. These include spatial resolution of the sampling units, allocation of effort in space, and collection of information about the observation process. On the analysis side, we review extensions of the basic occupancy models that are frequently necessary with atlas data, including methods for dealing with heterogeneity, non‐independent detections, false detections, and violation of the closure assumption. New technologies, such as cell‐phone apps and fixed remote detection devices, are revolutionizing citizen‐science projects. There is an opportunity to maximize the usefulness of the resulting datasets if the protocols are rooted in robust statistical designs and data analysis issues are being considered. Our review provides guidelines for designing new projects and an overview of the current methods that can be used to analyse data from such projects.

Methods in Ecology and Evolution

Ichthyoplankton abundance and variance in a large river system concerns for long-term monitoring

System-wide spatial patterns of ichthyoplankton abundance and variability were assessed in the upper Mississippi and lower Illinois rivers to address the experimental design and statistical confidence in density estimates. Ichthyoplankton was sampled from June to August 1989 in primary milieus (vegetated and non-vegated backwaters and impounded areas, main channels and main channel borders) in three navigation pools (8, 13 and 26) of the upper Mississippi River and in a downstream reach of the Illinois River. Ichthyoplankton densities varied among stations of similar aquatic landscapes (milieus) more than among subsamples within a station. An analysis of sampling effort indicated that the collection of single samples at many stations in a given milieu type is statistically and economically preferable to the collection of multiple subsamples at fewer stations. Cluster analyses also revealed that stations only generally grouped by their preassigned milieu types. Pilot studies such as this can define station groupings and sources of variation beyond an a priori habitat classification. Thus the minimum intensity of sampling required to achieve a desired statistical confidence can be identified before implementing monitoring efforts.

Illinois River, Mississippi River

Calibration of Precipitation-Runoff Modeling System (PRMS) to simulate prefire and postfire hydrologic response in the upper Rio Hondo Basin, New Mexico

The Precipitation-Runoff Modeling System (PRMS) is widely used to simulate the effects of climate, topography, land cover, and soils on landscape-level hydrologic responses and streamflow. The U.S. Geological Survey (USGS), in cooperation with the New Mexico Department of Homeland Security and Emergency Management, developed procedures to apply the PRMS model to simulate the effects of fire on hydrologic responses. A PRMS model was built of the upper Rio Hondo Basin from the headwaters to approximately 19 miles downstream from the USGS streamgage Rio Hondo above Chavez Canyon near Hondo, New Mexico, by using 24 hydrologic response units (HRUs), or hydrologically similar subareas, from the National Hydrologic Model. A quasi-graphical user interface was created to easily query and analyze published PRMS sensitivity-analysis data. Simulation of mean daily streamflow was most sensitive to parameters related to snowmelt or infiltration throughout the upper Rio Hondo Basin. In the basin’s eastern and northern HRUs, flashiness and timing of streamflow were most sensitive to interflow; in many western-basin HRUs (higher elevations), flashiness of streamflow was most sensitive to soil moisture parameters, and timing of streamflow was most sensitive to infiltration and evapotranspiration parameters. The PRMS model was calibrated for the fire-affected North Fork Eagle Creek subwatershed by comparing modeled to observed daily streamflow for the nonfrozen (May through October) period for a prefire and postfire time period. The prefire model was calibrated for the period 2007–12 before the 2012 fire, and the postfire model was calibrated for a 2-year (2014–15) period after the fire. Model parameterization combined manual adjustment of 8 parameters on the basis of prior knowledge and automated adjustment of the most sensitive parameters by using the Let Us Calibrate interface. A gridded, daily precipitation dataset that captured the spatial heterogeneity across the study watershed was used as the precipitation input for calibration. Model performance was assessed as satisfactory by using standard statistical measures for prefire and postfire periods. The calibrated model was run by using data from a single precipitation gage to better represent the effect of localized, extreme storms on postfire hydrologic response. The calibrated models for prefire and postfire conditions simulated streamflows with greater consistency than the uncalibrated model for the corresponding (prefire or postfire) period of hydrographic record. The effect of fire on streamflow was found to be primarily a shift from streamflow dominated by base flow prior to fire to streamflow dominated by surface runoff after fire.

New Mexico

Modeling landslide recurrence in Seattle, Washington, USA

To manage the hazard associated with shallow landslides, decision makers need an understanding of where and when landslides may occur. A variety of approaches have been used to estimate the hazard from shallow, rainfall-triggered landslides, such as empirical rainfall threshold methods or probabilistic methods based on historical records. The wide availability of Geographic Information Systems (GIS) and digital topographic data has led to the development of analytic methods for landslide hazard estimation that couple steady-state hydrological models with slope stability calculations. Because these methods typically neglect the transient effects of infiltration on slope stability, results cannot be linked with historical or forecasted rainfall sequences. Estimates of the frequency of conditions likely to cause landslides are critical for quantitative risk and hazard assessments. We present results to demonstrate how a transient infiltration model coupled with an infinite slope stability calculation may be used to assess shallow landslide frequency in the City of Seattle, Washington, USA. A module called CRF (Critical RainFall) for estimating deterministic rainfall thresholds has been integrated in the TRIGRS (Transient Rainfall Infiltration and Grid-based Slope-Stability) model that combines a transient, one-dimensional analytic solution for pore-pressure response to rainfall infiltration with an infinite slope stability calculation. Input data for the extended model include topographic slope, colluvial thickness, initial water-table depth, material properties, and rainfall durations. This approach is combined with a statistical treatment of rainfall using a GEV (General Extreme Value) probabilistic distribution to produce maps showing the shallow landslide recurrence induced, on a spatially distributed basis, as a function of rainfall duration and hillslope characteristics.

Washington

Disease‐structured N‐mixture models: A practical guide to model disease dynamics using count data

Obtaining inferences on disease dynamics (e.g., host population size, pathogen prevalence, transmission rate, host survival probability) typically requires marking and tracking individuals over time. While multistate mark–recapture models can produce high‐quality inference, these techniques are difficult to employ at large spatial and long temporal scales or in small remnant host populations decimated by virulent pathogens, where low recapture rates may preclude the use of mark–recapture techniques. Recently developed N ‐mixture models offer a statistical framework for estimating wildlife disease dynamics from count data. N ‐mixture models are a type of state‐space model in which observation error is attributed to failing to detect some individuals when they are present (i.e., false negatives). The analysis approach uses repeated surveys of sites over a period of population closure to estimate detection probability. We review the challenges of modeling disease dynamics and describe how N ‐mixture models can be used to estimate common metrics, including pathogen prevalence, transmission, and recovery rates while accounting for imperfect host and pathogen detection. We also offer a perspective on future research directions at the intersection of quantitative and disease ecology, including the estimation of false positives in pathogen presence, spatially explicit disease‐structured N ‐mixture models, and the integration of other data types with count data to inform disease dynamics. Managers rely on accurate and precise estimates of disease dynamics to develop strategies to mitigate pathogen impacts on host populations. At a time when pathogens pose one of the greatest threats to biodiversity, statistical methods that lead to robust inferences on host populations are critically needed for rapid, rather than incremental, assessments of the impacts of emerging infectious diseases.

Ecology and Evolution