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.

SourceResearch connectionFrom yearThrough year

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

At least 1,135 records · Page 63Linked to original sources

USGS1999

Element concentrations in bed sediment of the Yellowstone River basin, Montana, North Dakota, and Wyoming — A retrospective analysis

Chemical data for bed sediment were analyzed as part of the U.S. Geological Survey National Water-Quality Assessment Program investigation of the Yellowstone River Basin in parts of Montana, North Dakota, and Wyoming. The primary data set consisted of about 13,000 samples collected during 1974-79 for the National Uranium Resource Evaluation program. Data were available for 50 elements, although not all samples were analyzed for all elements. Element concentrations varied spatially and were associated with geologic settings or ecoregions. Factor analysis indicated three groups of associated elements: factor 1 elements were strongly correlated with basaltic rocks, factor 2 elements were strongly correlated with granitic rocks, and factor 3 elements were strongly correlated with carbonate rocks. Scores for factor 1 were highest for bed-sediment samples associated with volcanic rocks of Tertiary and Cretaceous age in the Absaroka volcanic field and crystalline rocks of Precambrian age in the Beartooth Mountains. Scores for factor 2 were highest for samples associated with volcanic rocks of Quaternary age on the Yellowstone Plateau, crystalline rocks of Precambrian age, and sedimentary rocks of Tertiary age in the Wyoming Basin ecoregion. Scores for factor 3 were highest in samples associated with sedimentary rocks of Paleozoic age and volcanic rocks of Cretaceous and Tertiary age. Descriptive statistics are presented to serve as a baseline for element concentrations in bed sediment associated with eight geologic settings or ecoregions in the study unit. Some of the concentrations of chromium, copper, lead, nickel, and zinc in bed-sediment samples from areas of crystalline rocks in the Beartooth Mountains and other formations in the western part of the study unit exceeded sediment-quality assessment values associated with toxic effects to aquatic life.

Idaho, Montana, North Dakota, South Dakota, Wyomin↗

USGS2025

A monitoring framework to assess forest bird population response to landscape scale mosquito suppression using the Incompatible Insect Technique

The Birds, Not Mosquitoes Monitoring and Support Science Working Group detailed methods for monitoring the population response of Hawaiian forest birds during implementation of the Incompatible Insect Technique (IIT) on the islands of Maui and Kauaʻi. The group prioritized methods for measuring the influence of mosquito suppression on populations within IIT treatment and control areas and identified focal species for IIT efficacy monitoring in birds. Three primary metrics were established to assess the impact of IIT on vulnerable species: population demography, density, and geographic range. Each metric can be evaluated using multiple methods. This report reviews those methods, with emphasis on approaches supported by pre-IIT baseline data and compatible with a before-after control-impact (BACI) study design for evaluating population responses over time. Focal avian species were selected based on population size estimates, fecundity, and disease susceptibility. We identified ʻākohekohe (Palmeria dolei), ʻiʻiwi (Drepanis coccinea), Maui ʻalauahio (Paroreomyza montana), Hawaiʻi ʻamakihi (Chlorodrepanis virens), Kauaʻi ʻamakihi (Chlorodrepanis stejnegeri), Kauaʻi ʻelepaio (Chasiempis sclateri), and ʻanianiau (Magumma parva) as focal species for monitoring population level response to disease suppression. Populations of kiwikiu (Pseudonestor xanthophrys), ʻakikiki (Oreomystis bairdi), akekeʻe (Loxops caeruleirostris), and the ʻiʻiwi population on Kauaʻi may be too small (e.g., <100 individuals) to effectively monitor, and it is unlikely that sufficient data can be collected from these birds to show IIT efficacy in a relatively short time frame (i.e., 5–10 years). Despite the logistical challenges to IIT implementation, there is potential to maintain disease-free status in individual populations of birds. Indeed, the continued existence of these critically endangered species in the wild within or near IIT treatment areas could be considered an accomplishment of IIT, given the current predictions for their extinction in the wild within 5–10 years. Demographic monitoring methods, including territory mapping, nest monitoring, mist-netting, and mark-recapture studies, provide direct evidence of survivorship and reproductive output. When combined with disease surveillance, these approaches could provide the most robust evidence of increased survivorship and productivity resulting from avian malaria suppression via IIT. However, demographic studies require several years of monitoring to achieve statistically robust BACI comparisons of survivorship and are more difficult to implement relative to other approaches. Given that these field efforts are labor-intensive and heavily reliant on personnel availability and funding, demographic monitoring could be conducted when adequate resources permit. On both Maui and Kauaʻi, passive acoustic monitoring (PAM) was identified as a priority method for monitoring the range, occupancy, and relative abundance of focal species. Autonomous recording units (ARUs) can record bird vocalizations in remote areas for several months. Innovative machine learning techniques permit rapid and semi-autonomous identification of most endemic honeycreepers on each island, maximizing sampling efficiencies and minimizing data processing costs. We predict mosquito suppression could support expansion of focal species into areas where disease transmission is currently excluding these species and expect acoustic monitoring data of focal species to reflect these spatial patterns. Additionally, the relative occupancy and call densities can be monitored temporally and spatially to assess the efficacy of IIT for supporting positive growth in vulnerable bird species. It is not yet clear if PAM is more effective than other methods, such as distance sampling, for detecting trends in the densities of rare species. However, the increased detections resulting from the larger sample size per observation point using ARUs will likely improve accuracy in detecting changes in species’ ranges. Collection of during and after treatment data within the BACI design could help to provide critical information to track avian population response, recovery, and potential range expansion related to IIT efforts. Point-transect distance sampling (point-counts) was prioritized as a method for monitoring population densities of focal species. Extensive historical sampling across focal species’ ranges provides a robust baseline for detecting change. These counts provide updated population densities and can be used to assess the distribution of focal species within IIT treatment areas. However, detecting subtle population changes with traditional distance sampling requires intensive spatial and temporal effort and may be less effective for rare species. To improve resolution, density surface modeling can integrate multiple data sources (e.g., point-counts, PAM, spot-mapping, and resightings) to estimate species-specific densities at finer spatial scales, including within and outside IIT treatment areas. This integrated modeling approach allows for detailed comparisons and may reveal early signs of recovery, including recolonization of formerly occupied sites. A coordinated monitoring strategy can allow managers to evaluate the success of mosquito suppression as a conservation intervention and support adaptive management in the face of emerging challenges.

Hawaii↗

USGS2016

Detecting significant change in stream benthic macroinvertebrate communities in wilderness areas

A major challenge in the biological monitoring of stream ecosystems in protected wilderness areas is discerning whether temporal changes in community structure are significantly outside of a reference condition that represents natural or acceptable annual variation in population cycles. Otherwise sites could erroneously be classified as impaired. Long-term datasets are essential for understanding these trends, to ascertain whether any changes in community structure significantly beyond the reference condition are permanent shifts or with time move back to within previous limits. To this end, we searched for long-term (>8 years) quantitative data sets of macroinvertebrate communities in wadeable rivers collected by similar methods and time of year in protected wilderness areas with minimal anthropogenic disturbance. Four geographic areas with datasets that met these criteria in the USA were identified, namely: McLaughlin Nature Reserve in California (1 stream), Great Smoky Mountains National Park in Tennesse-North Carolina (14 streams), Wind River Wilderness Areas in Wyoming (3 streams) and Denali National Park and Preserve in Alaska (6 streams). Two statistical approaches were applied: Taxonomic Distinctness (TD) to describe changes in diversity over time and non-metric multidimensional scaling (MDS) to describe changes over time in community persistence (Jaccards Index) and community stability (Bray–Curtis Index). Control charts were used to determine if years in MDS plots were significantly outside a reference condition. For Hunting Creek, TD showed three years outside natural variation which could be attributed to severe hydrological events but years outside the natural-variation funnel at sites in other geographical areas were inconsistent and could not be explained by environmental variables. TD identified simulated severe pollutant events which caused the removal of entire invertebrate assemblages but not simulated water temperature shifts. Within a region, both MDS analyses typically identified similar years as exceeding reference condition variation, illustrating the utility of the approach for identifying wider spatial scale effects that influence more than one stream. MDS responded to both simulated water temperature stress and a pollutant event, and generally outlying years on MDS plots could be explained by environmental variables, particularly higher precipitation. Multivariate control charts successfully identified whether shifts in community structure identified by MDS were significant and whether the shift represented a press disturbance (long-term change) or a pulse disturbance. We consider a combination of TD and MDS with control charts to be a potentially powerful tool for determining years significantly outside of a reference condition variation.

Alaska, California, North Carolina, Tennesse, Wyom↗

USGS2015

Evaluating a satellite-based seasonal evapotranspiration product and identifying its relationship with other satellite-derived products and crop yield: A case study for Ethiopia

Satellite-derived evapotranspiration anomalies and normalized difference vegetation index (NDVI) products from Moderate Resolution Imaging Spectroradiometer (MODIS) data are currently used for African agricultural drought monitoring and food security status assessment. In this study, a process to evaluate satellite-derived evapotranspiration (ETa) products with a geospatial statistical exploratory technique that uses NDVI, satellite-derived rainfall estimate (RFE), and crop yield data has been developed. The main goal of this study was to evaluate the ETa using the NDVI and RFE, and identify a relationship between the ETa and Ethiopia’s cereal crop (i.e., teff, sorghum, corn/maize, barley, and wheat) yields during the main rainy season. Since crop production is one of the main factors affecting food security, the evaluation of remote sensing-based seasonal ETa was done to identify the appropriateness of this tool as a proxy for monitoring vegetation condition in drought vulnerable and food insecure areas to support decision makers. The results of this study showed that the comparison between seasonal ETa and RFE produced strong correlation ( R 2 > 0.99) for all 41 crop growing zones in Ethiopia. The results of the spatial regression analyses of seasonal ETa and NDVI using Ordinary Least Squares and Geographically Weighted Regression showed relatively weak yearly spatial relationships ( R 2 < 0.7) for all cropping zones. However, for each individual crop zones, the correlation between NDVI and ETa ranged between 0.3 and 0.84 for about 44% of the cropping zones. Similarly, for each individual crop zones, the correlation ( R 2 ) between the seasonal ETa anomaly and de-trended cereal crop yield was between 0.4 and 0.82 for 76% (31 out of 41) of the crop growing zones. The preliminary results indicated that the ETa products have a good predictive potential for these 31 identified zones in Ethiopia. Decision makers may potentially use ETa products for monitoring cereal crop yields and early warning of food insecurity during drought years for these identified zones.

International Journal of Applied Earth Observation↗

USGS2021

A framework for assessing the ability to detect macroscale effects on fish growth

Various abiotic and biotic factors affect fish and their habitats at macroscales. For example, changes in global temperatures will likely alter demographic rates, including growth. However, to date, there is no statistical framework for assessing the ability to detect macroscale effects on fish growth under different sampling scenarios. We provide a generalized framework for calculating the frequentist and Bayesian power of detecting macroscale effects on fish growth. We illustrate this framework for a range of sampling scenarios that varied in the number of fish sampled per lake, the number of lakes sampled, and the magnitude of the temperature effect on growth for two case study species. However, the framework can be adapted to investigate other species, sampling scenarios, and environmental drivers. The ability to detect macroscale effects was more affected by the number of lakes sampled rather than the number of fish sampled from each lake. Confidently detecting macroscale effects likely requires sampling hundreds of lakes. This was true for both case study species, despite different life histories and extents of spatial variability in growth.

Canadian Journal of Fisheries and Aquatic Sciences↗

USGS2020

Spatially referenced models of streamflow and nitrogen, phosphorus, and suspended-sediment loads in streams of the midwestern United States

In this report, SPAtially Referenced Regression On Watershed attributes (SPARROW) models developed to describe long-term (2000–14) mean-annual streamflow, total nitrogen (TN), total phosphorus (TP), and suspended-sediment (SS) transport in streams of the Midwestern part of the United States (the Mississippi River, Great Lakes, and Red River of the North Basins) are described. The nutrient and suspended-sediment models have a base year of 2012, which means they were developed based on source inputs and management practices similar to those existing during or near 2012 and average hydrological conditions detrended to 2012 (2000–14), whereas the streamflow model has base years of 2000–14, which means it was developed based on the average input precipitation minus actual evapotranspiration from 2000 to 2014. In developing the models, several updates and improvements were made to the data inputs and statistical approaches used to calibrate/develop the models from those used in the previous 2002 SPARROW models. The 2012 SPARROW models were constructed using a higher resolution stream network, which resulted in a mean catchment size of 2.7 square kilometers compared to 480 square kilometers in the 2002 models; more detailed and updated wastewater treatment plant contribution estimates; inputs from background phosphorus sources that were not included in the 2002 model; and more accurate loads for calibration that were computed using a modified Beale ratio-estimator technique whenever no trend in load was determined. Statistical approaches were added to compensate for the unequal effect of each monitoring site during the calibration process by adjusting for the fraction of the basin included in other upstream monitored sites (nested share) and thinning the calibration sites if a negative statistical correlation between nearby sites was determined. Results from 2012 SPARROW models describe how much of each water, TN, TP, and SS source was delivered to the stream network, and the major landscape factors that affected their delivery. Atmospheric deposition and natural (background) sources of TN and TP, respectively, were the dominant sources in anthropogenically unaffected areas (especially in the Rocky Mountains and north-central areas of the Midwest), whereas fertilizers, manure, and fixation were dominant sources in agricultural areas, especially in the Corn Belt and near the Mississippi River. Urban sources of TN and TP were typically localized, but they were still important for some large areas, especially the Lake Erie Basin. All of the land-to-water delivery variables in the nutrient and sediment SPARROW models, such as runoff, soil erodibility, basin slope, and the amount of tile drains, are commonly included in process-driven models. In the SPARROW TN and TP models, best management practices (BMPs) reduced the delivery of these nutrients to streams. Long-term mean-annual flows and nutrient and sediment loads were simulated in streams throughout the Midwest. The simulated flows from the SPARROW flow model were used in the SPARROW TN, TP, and SS models to help describe nutrient and sediment transport from the watershed and through the stream network. Outputs from the TN, TP, and SS models describe loads and yields of these constituents throughout the Midwest, and from major drainage basins throughout the Midwest. Highest TN, TP, and SS yields and delivered yields were from the Lake Erie, Ohio River, Upper Mississippi River, and Lower Mississippi River Basins, whereas lowest yields were spread over most other areas. Losses during downstream delivery resulted in part of the TN, TP, and SS that reach the stream network not reaching the downstream receiving bodies: 14, 15, and 28 percent of the TN, TP, and SS, respectively, are lost during delivery to the Great Lakes and 19, 23, and 52 percent of the TN, TP, and SS, respectively, are lost during delivery to the Gulf of Mexico. The largest losses of nutrients and sediments during transport were in the Missouri and Arkansas River Basins. Information from these SPARROW models can help guide nutrient and sediment reduction strategies throughout the Midwest. Model results provide information on what may be the most appropriate general type of actions to reduce total loading by describing the relative importance of each source, and where to most efficiently place the efforts to reduce loading by describing the distribution of nutrient and sediment loading. By implementing management efforts addressing the major sources of the loads in areas contributing the highest loads, it may be possible to reduce nutrient loading throughout the Mississippi River Basin and thus reduce the size of the hypoxic zone in the Gulf of Mexico; reduce nutrient loading into lakes, and thus reduce the occurrence of harmful algal blooms; and reduce sediment losses, and thus improve the benthic habitat in streams and rivers throughout the Midwest.

Scientific Investigations Report↗

USGS2007

Spatial variability in water-balance model performance in the conterminous United States

A monthly water-balance (WB) model was tested in 44 river basins from diverse physiographic and climatic regions across the conterminous United States (U.S.). The WB model includes the concepts of climatic water supply and climatic water demand, seasonality in climatic water supply and demand, and soil-moisture storage. Exhaustive search techniques were employed to determine the optimal set of precipitation and temperature stations, and the optimal set of WB model parameters to use for each basin. It was found that the WB model worked best for basins with: (1) a mean elevation less than 450 meters or greater than 2000 meters, and/or (2) monthly runoff that is greater than 5 millimeters (mm) more than 80 percent of the time. In a separate analysis, a multiple linear regression (MLR) was computed using the adjusted R-square values obtained by comparing measured and estimated monthly runoff of the original 44 river basins as the dependent variable, and combinations of various independent variables [streamflow gauge latitude, longitude, and elevation; basin area, the long-term mean and standard deviation of annual precipitation; temperature and runoff; and low-flow statistics (i.e., the percentage of months with monthly runoff that is less than 5 mm)]. Results from the MLR study showed that the reliability of a WB model for application in a specific region can be estimated from mean basin elevation and the percentage of months with gauged runoff less than 5 mm. The MLR equations were subsequently used to estimate adjusted R-square values for 1,646 gauging stations across the conterminous U.S. Results of this study indicate that WB models can be used reliably to estimate monthly runoff in the eastern U.S., mountainous areas of the western U.S., and the Pacific Northwest. Applications of monthly WB models in the central U.S. can lead to uncertain estimates of runoff.

Conterminous United States↗

USGS2011

Long-term post-fire effects on spatial ecology and reproductive output of female Agassiz’s desert tortoises ( Gopherus agassizii ) at a wind energy facility near Palm Springs, California, USA

We studied the long-term response of a cohort of eight female Agassiz’s desert tortoises ( Gopherus agassizii ) during the first 15 years following a large fire at a wind energy generation facility near Palm Springs, California, USA. The fire burned a significant portion of the study site in 1995. Tortoise activity areas were mapped using minimum convex polygons for a proximate post-fire interval from 1997 to 2000, and a long-term post-fire interval from 2009 to 2010. In addition, we measured the annual reproductive output of eggs each year and monitored the body condition of tortoises over time. One adult female tortoise was killed by the fire and five tortoises bore exposure scars that were not fatal. Despite predictions that tortoises would make the short-distance movements from burned to nearby unburned habitats, most activity areas and their centroids remained in burned areas for the duration of the study. The percentage of activity area burned did not differ significantly between the two monitoring periods. Annual reproductive output and measures of body condition remained statistically similar throughout the monitoring period. Despite changes in plant composition, conditions at this site appeared to be suitable for survival of tortoises following a major fire. High productivity at the site may have buffered tortoises from the adverse impacts of fire if they were not killed outright. Tortoise populations at less productive desert sites may not have adequate resources to sustain normal activity areas, reproductive output, and body conditions following fire.

California↗

USGS2005

Distribution patterns of mercury in Lakes and Rivers of northeastern North America

We assembled 831 data points for total mercury (Hg t ) and 277 overlapping points for methyl mercury (CH 3 Hg + ) in surface waters from Massachussetts, USA to the Island of Newfoundland, Canada from State, Provincial, and Federal government databases. These geographically indexed values were used to determine: (a) if large-scale spatial distribution patterns existed and (b) whether there were significant relationships between the two main forms of aquatic Hg as well as with total organic carbon (TOC), a well know complexer of metals. We analyzed the catchments where samples were collected using a Geographical Information System (GIS) approach, calculating catchment sizes, mean slope, and mean wetness index. Our results show two main spatial distribution patterns. We detected loci of high Hg t values near urbanized regions of Boston MA and Portland ME. However, except for one unexplained exception, the highest Hg t and CH 3 Hg + concentrations were located in regions far from obvious point sources. These correlated to topographically flat (and thus wet) areas that we relate to wetland abundances. We show that aquatic Hg t and CH 3 Hg + concentrations are generally well correlated with TOC and with each other. Over the region, CH 3 Hg + concentrations are typically approximately 15% of Hg t . There is an exception in the Boston region where CH 3 Hg + is low compared to the high Hg t values. This is probably due to the proximity of point sources of inorganic Hg and a lack of wetlands. We also attempted to predict Hg concentrations in water with statistical models using catchment features as variables. We were only able to produce statistically significant predictive models in some parts of regions due to the lack of suitable digital information, and because data ranges in some regions were too narrow for meaningful regression analyses.

Maine, Massachusetts, New Brunswick, New Hampshire↗

USGS2018

Wetlands inform how climate extremes influence surface water expansion and contraction

Effective monitoring and prediction of flood and drought events requires an improved understanding of how and why surface water expansion and contraction in response to climate varies across space. This paper sought to (1) quantify how interannual patterns of surface water expansion and contraction vary spatially across the Prairie Pothole Region (PPR) and adjacent Northern Prairie (NP) in the United States, and (2) explore how landscape characteristics influence the relationship between climate inputs and surface water dynamics. Due to differences in glacial history, the PPR and NP show distinct patterns in regards to drainage development and wetland density, together providing a diversity of conditions to examine surface water dynamics. We used Landsat imagery to characterize variability in surface water extent across 11 Landsat path/rows representing the PPR and NP (images spanned 1985–2015). The PPR not only experienced a 2.6-fold greater surface water extent under median conditions relative to the NP, but also showed a 3.4-fold greater change in surface water extent between drought and deluge conditions. The relationship between surface water extent and accumulated water availability (precipitation minus potential evapotranspiration) was quantified per watershed and statistically related to variables representing hydrology-related landscape characteristics (e.g., infiltration capacity, surface storage capacity, stream density). To investigate the influence stream connectivity has on the rate at which surface water leaves a given location, we modeled stream-connected and stream-disconnected surface water separately. Stream-connected surface water showed a greater expansion with wetter climatic conditions in landscapes with greater total wetland area, but lower total wetland density. Disconnected surface water showed a greater expansion with wetter climatic conditions in landscapes with higher wetland density, lower infiltration and less anthropogenic drainage. From these findings, we can expect that shifts in precipitation and evaporative demand will have uneven effects on surface water quantity. Accurate predictions regarding the effect of climate change on surface water quantity will require consideration of hydrology-related landscape characteristics including wetland storage and arrangement.

Prairie Pothole Regino↗

USGS2003

Preliminary preview for a geographic and monitoring program project; a review of point source-nonpoint source effluent trading/offset systems in watersheds

Watershed-based trading and offset systems are being developed to improve policy-maker?s and regulator?s ability to assess nonpoint source impacts in watersheds and to evaluate the efficacy of using market-incentive programs for preserving environmental quality. An overview of the history of successful and failed trading programs throughout the United States suggests that certain political, economic, and scientific conditions within a temporal and spatial setting help meet water quality standards. The current lack of spontaneous trading among dischargers does not mean that a marketable permit trading system is an inherently inefficient regulatory approach. Rather, its infrequent use is the result of institutional and informational barriers. Improving and refining the earth science information and technologies may help determine whether trading is a suitable policy for improving water quality. However, it is debatable whether or not environmental information is the limiting factor. This paper reviews additional factors affecting the potential for instituting a trading policy. The motivation for investigating and reviewing the history of offsets and trading was inspired by a project in the preliminary stages being developed by U.S. Geological Survey Western Geographic Science Center and the Environmental Protection Agency Region IX. An offset feasibility study will be an integrated, map-based approach that incorporates environmental, economic, and statistical information to investigate the potential for using offsets to meet mercury Total Maximum Daily Loads in the Sacramento River watershed. A regional water-quality offset program is being studied that may help known point sources reduce mercury loading more cost effectively by the remediation of abandoned mines or other diffuse sources as opposed to more costly treatment at their own sites. An efficient offset program requires both a scientific basis and methods to translate that science into a regulatory decision framework.

Open-File Report↗

USGS2017

Characterization of water quality and suspended sediment during cold-season flows, warm-season flows, and stormflows in the Fountain and Monument Creek watersheds, Colorado, 2007–2015

From 2007 through 2015, the U.S. Geological Survey, in cooperation with Colorado Springs City Engineering, conducted a study in the Fountain and Monument Creek watersheds, Colorado, to characterize surface-water quality and suspended-sediment conditions for three different streamflow regimes with an emphasis on characterizing water quality during storm runoff. Data collected during this study were used to evaluate the effects of stormflows and wastewater-treatment effluent discharge on Fountain and Monument Creeks in the Colorado Springs, Colorado, area. Water-quality samples were collected at 2 sites on Upper Fountain Creek, 2 sites on Monument Creek, 3 sites on Lower Fountain Creek, and 13 tributary sites during 3 flow regimes: cold-season flow (November–April), warm-season flow (May–October), and stormflow from 2007 through 2015. During 2015, additional samples were collected and analyzed for Escherichia coli ( E. coli ) during dry weather conditions at 41 sites, located in E. coli impaired stream reaches, to help identify source areas and scope of the impairment. Concentrations of E. coli , total arsenic, and dissolved copper, selenium, and zinc in surface-water samples were compared to Colorado in-stream standards. Stormflow concentrations of E. coli frequently exceeded the recreational use standard of 126 colonies per 100 milliliters at main-stem and tributary sites by more than an order of magnitude. Even though median E. coli concentrations in warm-season flow samples were lower than median concentrations in storm-flow samples, the water quality standard for E. coli was still exceeded at most main-stem sites and many tributary sites during warm-season flows. Six samples (three warm-season flow and three stormflow samples) collected from Upper Fountain Creek, upstream from the confluence of Monument Creek, and two stormflow samples collected from Lower Fountain Creek, downstream from the confluence with Monument Creek, exceeded the acute water-quality standard for total arsenic of 50 micrograms per liter. All concentrations of dissolved copper, selenium, and zinc measured in samples were below the water-quality standard. Concentrations of dissolved nitrate plus nitrite generally increased from upstream to downstream during all flow periods. The largest downstream increase in dissolved nitrate plus nitrite concentration was measured between sites 07103970 and 07104905 on Monument Creek. All but one tributary that drain into Monument Creek between the two sites had higher median nitrate plus nitrite concentrations than the nearest upstream site on Monument Creek, site 07103970 (MoCr_Woodmen). Increases in the concentration of dissolved nitrate plus nitrite were also evident below wastewater treatment plants located on Fountain Creek. Most stormflow concentrations of dissolved trace elements were smaller than concentrations from cold-season flow or warm-season samples. However, median concentrations of total arsenic, lead, manganese, nickel, and zinc generally were much larger during periods of stormflow than during cold-season flow or warm-season fl. Median concentrations of total arsenic, total copper, total lead, dissolved and total manganese, total nickel, dissolved and total selenium, and dissolved and total zinc concentrations increased from 1.5 to 28.5 times from site 07103700 (FoCr_Manitou) to 07103707 (FoCr_8th) during cold-season and warm-season flows, indicating a large source of trace elements between these two sites. Both of these sites are located on Fountain Creek, upstream from the confluence with Monument Creek. Median suspended-sediment concentrations and median suspended-sediment loads increased in the downstream direction during all streamflow regimes between Monument Creek sites 07103970 (MoCr_Woodmen) and 07104905 (MoCr_Bijou); however, statistically significant increase (p-value less than 0.05) were only present during warm-season flow and stormflow. Significant increases in median suspended sediment concentrations were measured during cold-season flow and warm-season flow between Upper Fountain Creek site 07103707 (FoCr_8th) and Lower Fountain Creek site 07105500 (FoCr_Nevada) because of inflows from Monument Creek with higher suspended-sediment concentrations. Median suspended-sediment concentrations between sites 07104905 (MoCr_Bijou) and 07105500 (FoCr_Nevada) increased significantly during warm-season flow but showed no significant differences during cold-season flow and stormflow. Significant decreases in median suspended-sediment concentrations were measured between sites 07105500 (FoCr_Nevada) and 07105530 (FoCr_Janitell) during all flow regimes. Suspended-sediment concentrations, discharges, and yields associated with stormflow were significantly larger than those associated with warm-season flow. Although large spatial variations in suspended-sediment yields occurred during warm-season flows, the suspended-sediment yield associated with stormflow were as much as 1,000 times larger than the suspended-sediment yields that occurred during warm-season flow.

Colorado↗

USGS2009

Estimating Locations of Perennial Streams in Idaho Using a Generalized Least-Squares Regression Model of 7-Day, 2-Year Low Flows

Many State and Federal agencies use information regarding the locations of streams having intermittent or perennial flow when making management and regulatory decisions. For example, the application of some Idaho water quality standards depends on whether streams are intermittent. Idaho Administrative Code defines an intermittent stream as one having a 7-day, 2-year low flow (7Q2) less than 0.1 ft3/s. However, there is a general recognition that the cartographic representation of perennial/intermittent status of streams on U.S. Geological Survey (USGS) topographic maps is not as accurate or consistent as desirable from one map to another, which makes broad management and regulatory assessments difficult and inconsistent. To help resolve this problem, the USGS has developed a methodology for predicting the locations of perennial streams based on regional generalized least-squares (GLS) regression equations for Idaho streams for the 7Q2 low-flow statistic. Using these regression equations, the 7Q2 streamflow may be estimated for naturally flowing streams in most areas in Idaho. The use of these equations in conjunction with a geographic information system (GIS) technique known as weighted flow accumulation allows for an automated and continuous estimation of 7Q2 streamflow at all points along stream reaches. The USGS has developed a GIS-based map of the locations of streams in Idaho with perennial flow based on a 7Q2 of 0.1 ft3/s and a transition zone of plus or minus 1 standard error. Idaho State cooperators plan to use this information to make regulatory and water-quality management decisions. Originally, 7Q2 equations were developed for eight regions of similar hydrologic characteristics in the study area, using long-term data from 234 streamflow-gaging stations. Equations in five of the regions were revised based on spatial patterns observed in the initial perennial streams map and unrealistic behavior of the equations in extrapolation. The standard errors of prediction for the final equations ranged from a minimum of +75.0 to -42.9 percent in the central part of the study area to a maximum of +277 to -73.5 percent in the southern part of the study area. The equations are applicable only to unregulated, naturally-flowing streams and may produce unreliable results outside the range of explanatory variables used for equation development. Extrapolation outside the range of available data was necessary, however, to predict perennial flow initiation points and transition zones along stream reaches. The map of perennial streams was evaluated by comparing predicted stream classifications with four independent datasets, including field observations by other government agencies. Overall, 81 percent of the comparison data points agreed with the USGS perennial streams model. Regions with the highest number of disagreements had a high percentage of mountainous and forested area with potential mountain front recharge zones, and regions with the highest agreements had a high percentage of low gradient, low elevation area. As a whole, the USGS model predicted a higher number of perennial streams than predictions made with the independent datasets. Some disagreements were due to poor site location coordinates, timing of the comparison site visits during unusually wet or dry years, discrepancies in classification criteria, and variable ground water contributions to flow in some areas. The Idaho Department of Environmental Quality Beneficial Use Reconnaissance Program (BURP) dataset is considered the most representative dataset for comparison because it covered a range of climate conditions and the number of sites visited were consistent from year to year during the study period. Eighty-five percent of BURP comparison data points agreed with the USGS perennial streams model. Although site-specific flow data may be needed to correctly classify streams in some areas, this information rarely is available and is not always practical to o

Scientific Investigations Report↗

USGS2020

Comparing methods to estimate the proportion of turbine-induced bird and bat mortality in the search area under a road and pad search protocol

Estimating bird and bat mortality at wind facilities typically involves searching for carcasses on the ground near turbines. Some fraction of carcasses inevitably lie outside the search plots, and accurate mortality estimation requires accounting for those carcasses using models to extrapolate from searched to unsearched areas. Such models should account for variation in carcass density with distance, and ideally also for variation with direction (anisotropy). We compare five methods of accounting for carcasses that land outside the searched area (ratio, weighted distribution, non-parametric, and two generalized linear models ( glm )) by simulating spatial arrival patterns and the detection process to mimic observations which result from surveying only, or primarily, roads and pads (R&P) and applying the five methods. Simulations vary R&P configurations, spatial carcass distributions (isotropic and anisotropic), and per turbine fatality rates. Our results suggest that the ratio method is less accurate with higher variation relative to the other four methods which all perform similarly under isotropy. All methods were biased under anisotropy; however, including direction covariates in the glm method substantially reduced bias. In addition to comparing methods of accounting for unsearched areas, we suggest a semiparametric bootstrap to produce confidence-based bounds for the proportion of carcasses that land in the searched area.

Environmental and Ecological Statistics↗

USGS2021

Distribution, abundance and spatial variability of microplastic pollution on the surface of Lake Superior

In 2014, 94 paired neuston net samples (0.5 mm mesh) were collected from the surface waters of Lake Superior. These samples comprise the most comprehensive surface water survey for microplastics of any of the Great Lakes to date, and the first to employ double net trawls. Microplastic abundance estimates showed wide variability, ranging between 4000 to more than 100,000 particles/km 2 with most locations having abundances between 20,000 to 50,000 particles/km 2 . The average abundance in Lake Superior was ~30,000 particles/km 2 which was similar to previous estimates within this Laurentian Great Lake and suggests a total count of more than 2.4 billion (1.7 to 3.3 billion, 95% confidence interval) particles across the lake’s surface. Distributions of plastic particles, characterized by size fraction and type, differed between nearshore and offshore samples, and between samples collected in the eastern versus western portion of the lake. Most of the particles found were fibers (67%), and most (62%) were contained in the smallest classified size fraction (0.50–1 mm). The most common type of polymer found was polyethylene (51%), followed by polypropylene (19%). This is consistent with global plastics production and results obtained from other studies. No statistically significant difference was detected between the paired net samples, indicating that single net sampling should produce a representative estimate of microplastic particle abundance and distribution within a body of water.

Lake Superior↗

USGS1998

Analysis of the influence of spatial pattern in habitat selection studies

Design and analysis of wildlife habitat selection studies typically do not assess the effect of spatial pattern on the habitat selection process. Effects of landscape scale pattern on habitat selection cannot be accomplished without replicate study areas, because pattern is a single, albeit multifaceted, attribute of an area. For a single area, however, the influence of pattern-related characteristics, such as shape and edge shared with adjacent patches, can be estimated by using GLIM (McCullough and Neider 1983) procedures to model patch-specific frequency counts of animal use as a function of these parameters. This approach is evaluated and illustrated with simulated breeding-bird counts in a South Carolina study area for which a GIS land cover classification is available. A related technique for evaluating whether movement from patch to patch is selective is developed and illustrated for designs that involve collection of trajectory data from monitored individuals. These designs and analyses are feasible given current GIS and GPS technology. Statistical inferences from habitat selection studies should be interpreted within the context of a range of scales at which animals differentiate between patch attributes.

Journal of Agricultural, Biological, and Environme↗

USGS2021

Introduced mangroves along the coast of Moloka‘i, Hawai‘i may represent novel habitats for megafaunal communities

Mangrove forests are prevalent along tropical/subtropical coastlines and provide valuable ecosystem services including coastline stabilization, storm impact reduction, and enhanced coastal productivity. However, mangroves were absent from the Hawaiian Islands and their introduction to Moloka‘i in 1902 has provided an opportunity to examine their unique influence on coastal landscapes. Previous studies indicate an inability of native detritivores to utilize tannin-rich substrates, yielding poor cycling of mangrove-derived detritus in Hawaiian tidal zones. We hypothesize that in addition to altering detrital inputs, introduced mangroves facilitate the persistence of introduced species in the Hawaiian coastal zone by providing novel habitat for juvenile megafauna. To determine whether mangrove-dominated tidal zones harbor megafaunal assemblages distinct from open sandflats, we sampled in two mangrove (M1 and M2) and two adjacent sandflat (S1 and S2) sites along the southern coast of Moloka‘i, where the most mature mangrove forests occur in Hawai‘i. There were no statistical differences in total abundances between M1 and M2 or S1 and S2; therefore, results from individual deployments were pooled across the sites in order to conduct between-habitat (mangrove vs. sandflat) comparisons. Our mangrove study site had significantly higher abundances of megafauna, including several shrimp and crab species, compared to the sandflat site. The community composition within the mangrove site differed from the sandflat site, including higher abundances of non-native mangrove crabs ( Scylla serrata ), as well as native fish Bathygobius cocosensis and crustaceans ( Thalamita crenata , Palaemon pacificus , P. debilis ) than in the sandflat site, indicating that the mangrove site may provide niches for both invasive and native species. In addition, mean body length for several similar species was smaller in the mangrove site than in the sandflat site, suggesting that these mangroves may be providing a habitat for juvenile species. While our study was spatially limited to two mangrove and two adjacent sandflat sites, our results suggest that introduced mangroves in Moloka‘i may support small-bodied, native, and non-native megafauna, influencing coastal Hawaiian trophic dynamics. Our case study provides a baseline for megafaunal fish and invertebrate communities present prior to non-native mangrove removal as well as for monitoring potential community changes following expansion of mangrove habitats due to climate change.

Hawaii↗

USGS2025

Spatially explicit power analyses to inform occupancy‐based multi‐species wildlife monitoring programmes

1. Current and accurate information on wildlife populations is integral to successful biodiversity management and conservation globally. Nevertheless, many monitoring programs fail in their attempts to accurately monitor populations of interest due to interlinked issues including insufficient sample sizes, inappropriate duration, lack of reproducibility, and lack of clearly stated objectives. These common pitfalls could be avoided through the elicitation of explicit monitoring objectives and the a priori use of simulations to inform minimum sampling design requirements to meet said objectives. 2. Here, we provide a blueprint for using spatially explicit power analyses to inform the design and implementation of multi-species monitoring programs on landscape-scales. As a demonstration, we used spatially explicit simulations to devise a suitable sampling regime to meet clearly specified monitoring objectives in New York State: to use annual occupancy-based monitoring to be able to detect 25% and 50% changes in abundance of populations over five- and ten- year periods for all species of management interest in New York State, USA. We focused our simulation efforts on three challenging focal species (black bear, Ursus americanus, bobcat, Lynx rufus, and American marten, Martes americana) that differ notably in their morphology, life histories, space use, detection probability, habitat suitability, and population sizes/trajectories, and thus provide extremes in the challenges presented when it comes to sampling appropriately to detect changes in abundance. 3. Our simulations demonstrate variable context dependent trade-offs in sampling designs (i.e. number of sites [J] and number of sampling occasions [K]), and identify necessary minimum detection probabilities that must be attained to achieve statistical power to detect changes of varying magnitudes in populations of varying sizes in the three focal species. The simulations also highlight that monitoring population increases is likely beyond the reach of occupancy-based monitoring programs for wide-ranging or locally abundant species. 4. Synthesis and applications : We combine the results from the single-species simulations to produce a multi-species sampling design that meets the specified objectives for all three species. While the case study is centered on developing a multi-species sampling regime for New York State, it provides a reproducible step-by-step framework using established methods for wildlife managers and other practitioners to inform their own context- and objective- specific multi-species occupancy-based monitoring programs.

New York↗