USGS Science⌕ Search

SEARCH · USGS Science

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 541 records · Page 30Linked to original sources

Space-based monitoring of penguin diet links sea ice, food webs, and population change

Rapid environmental change is reshaping Antarctic ecosystems through shifts in temperature and sea ice dynamics, with implications for species such as the iconic Adélie penguin ( Pygoscelis adeliae ). Because Antarctic predators rely on a relatively small number of key prey resources in the region 1 whose life histories are intrinsically linked to sea ice, 2 , 3 these abiotic changes are expected to cascade through food webs. However, the scale and remoteness of Antarctica have limited efforts to link environmental change to ecological responses across the continent. Combining tools from imaging spectroscopy, stable isotope analysis, and hierarchical statistical modeling, we reconstructed Adélie penguin diet across the entirety of the species’ global range over a three-decade period (1984–2013) using satellite imagery from the Landsat program. This approach leveraged the distinct spectral properties (a generalized measure of color) of penguin guano as observed by satellites, 4 , 5 which vary according to penguin diet. 6 , 7 We found pronounced differences in diet across the continent and that year-to-year dietary shifts were strongly associated with sea ice dynamics. Higher sea ice corresponded to more fish-based diets, whereas lower sea ice was associated with a greater reliance on krill. Additionally, spatial differences in penguin diet were associated with long-term trends in penguin abundance, linking food web processes to large-scale population dynamics. This study represents the first use of satellite observations to capture trophic dynamics at continental and decadal scales and highlights how environmental change might restructure Antarctic food webs and influence the future of a key sentinel species.

Current Biology↗

Tracing metal sources and groundwater flow paths in the Upper Animas River watershed using rare earth elements and stable isotopes

Groundwater flow paths and processes that govern metal mobility and transport are difficult to characterize in mountainous bedrock watersheds. Despite the difficulty in holistic characterization, conceptual understanding of subsurface hydrologic and geochemical processes is key to developing remediation plans for locations affected by acid mine drainage, such as the Upper Animas River watershed in southwestern Colorado, USA. Stable isotopes of water and rare earth elements were utilized to evaluate groundwater flow and metal sources within this complex catchment. Stable isotope samples collected from draining mine adits and springs display systematic spatial variation wherein sample sites at higher elevations have greater seasonal variability than sites at lower elevations. The Upper Cement Creek watershed, where multiple draining mines are present, displays the lowest seasonal variation in stable isotopic signatures, potentially indicating the presence of a large, well-mixed volume of groundwater storage or interbasin groundwater flow. Rare earth elements display statistically significant variation between different alteration styles in the catchment. Overprinting of regional propylitic alteration is evident based on enrichment of middle rare earth elements in acidic springs and mines that are not spatially associated with surficial exposures of acid generating alteration styles. Europium anomaly and middle rare earth enrichment signatures from two flooded mine tunnels on opposite sides of a watershed divide indicate connections to the same subsurface flooded mine workings.

Colorado↗

Habitat classification modeling with incomplete data: Pushing the habitat envelope

Habitat classification models (HCMs) are invaluable tools for species conservation, land-use planning, reserve design, and metapopulation assessments, particularly at broad spatial scales. However, species occurrence data are often lacking and typically limited to presence points at broad scales. This lack of absence data precludes the use of many statistical techniques for HCMs. One option is to generate pseudo-absence points so that the many available statistical modeling tools can be used. Traditional techniques generate pseudoabsence points at random across broadly defined species ranges, often failing to include biological knowledge concerning the species-habitat relationship. We incorporated biological knowledge of the species-habitat relationship into pseudo-absence points by creating habitat envelopes that constrain the region from which points were randomly selected. We define a habitat envelope as an ecological representation of a species, or species feature's (e.g., nest) observed distribution (i.e., realized niche) based on a single attribute, or the spatial intersection of multiple attributes. We created HCMs for Northern Goshawk (Accipiter gentilis atricapillus) nest habitat during the breeding season across Utah forests with extant nest presence points and ecologically based pseudo-absence points using logistic regression. Predictor variables were derived from 30-m USDA Landfire and 250-m Forest Inventory and Analysis (FIA) map products. These habitat-envelope-based models were then compared to null envelope models which use traditional practices for generating pseudo-absences. Models were assessed for fit and predictive capability using metrics such as kappa, thresholdindependent receiver operating characteristic (ROC) plots, adjusted deviance (Dadj2), and cross-validation, and were also assessed for ecological relevance. For all cases, habitat envelope-based models outperformed null envelope models and were more ecologically relevant, suggesting that incorporating biological knowledge into pseudo-absence point generation is a powerful tool for species habitat assessments. Furthermore, given some a priori knowledge of the species-habitat relationship, ecologically based pseudo-absence points can be applied to any species, ecosystem, data resolution, and spatial extent. ?? 2007 by the Ecological Society of America.

Ecological Applications↗

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↗