SEARCH · Search USGS

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,063 records · Page 59Linked to original sources

USGS2026

Hydrologic variability drives environmental and geospatial relationships in Smallmouth Bass (Micropterus dolomieu) distribution

Hydrologic variation is a primary driver of stream ecosystems. Changing hydrology can lead to assemblage shifts and alterations in suitable habitat for freshwater species. As climate change is predicted to alter flow patterns in addition to increasing water temperatures, insight into relationships between species occupancy, hydrology, and temperature is critical for understanding current and future distributions. We examined how hydrologic variability, temperature, and other environmental variables interact to influence Micropterus dolomieu (Smallmouth Bass) occurrence. We used Spatial Stream Network models, allowing for the incorporation of spatial autocorrelation along streams' unique dendritic network, to examine Smallmouth Bass occupancy across a range of hydrologic variation in the Ozark-Ouachita Interior Highlands, USA. Hydrologic variation was the main driver of Smallmouth Bass occurrence, with occurrence more likely in groundwater streams with low hydrologic variation and high flow permanence. For groundwater streams, occurrence was positively associated with summer stream temperature and negatively associated with annual stream temperature. As variation increased, more variables showed significant relationships with occurrence. Distance metrics were important for all models, however as hydrologic disturbance increased, flow connected distance played a lesser role and stream distance played a greater role. Hydrologic variability was the overarching determinant of Smallmouth Bass occurrence and strongly influenced the predictive importance of environmental variables and geospatial relationships. Greater hydrologic variability resulted in stronger statistical relationships between occurrence and environmental variables and an increased importance of system connectivity. As climate change alters hydrologic processes and streams become more variable, understanding and accounting for these shifting relationships is essential.

Arkansas, Kansas, Missouri, Oklahoma↗

USGS2013

Effects of surface applications of biosolids on groundwater quality and trace-element concentrations in crops near Deer Trail, Colorado, 2004-2010

The U.S. Geological Survey (USGS), in cooperation with Metro Wastewater Reclamation District (Metro District), studied biosolids composition and the effects of biosolids applications on groundwater quality and trace-element concentrations in crops of the Metro District properties near Deer Trail, Colorado, during 2004 through 2010. Priority parameters for each monitoring component included the nine trace elements regulated by Colorado for biosolids (arsenic, cadmium, copper, lead, mercury, molybdenum, nickel, selenium, and zinc); other constituents also were analyzed. All concentrations for the priority parameters in monthly biosolids samples were less than Colorado regulatory limits, and the concentrations were relatively consistent. Biosolids likely were the largest source of nitrogen and phosphorus on the Metro District properties. Plutonium isotopes were not detected in the biosolids, but many organic wastewater compounds (organic wastewater compounds: wastewater indicators, pharmaceuticals, and hormones) were detected in substantial concentrations relative to minimum reporting levels and various surface-water concentrations. Bismuth, copper, mercury, nitrogen, phosphorus, silver, biogenic sterols, detergent degradates, disinfectants, fire retardants, fragrances, pharmaceuticals, and plasticizers would be the most likely biosolids signature to indicate the presence of Metro District biosolids in soil or streambed sediment from the study area. Antimony, cadmium, cobalt, copper, molybdenum, nickel, nitrogen, phosphorus, selenium, tungsten, vanadium, zinc, detergent degradates, disinfectants, fire retardants, fragrances, pharmaceuticals or their degradates, and plasticizers would be the most likely biosolids signature for groundwater and surface water in the study area. More biosolids-signature components detected and larger concentration differences from untreated materials, baseline, and blank samples indicate more evidence of biosolids presence or effects. Although the inorganic constituent concentrations were relatively large in samples from one monitoring well, the concentrations of organic wastewater compounds in groundwater samples were not correspondingly large. Concentrations of organic wastewater compounds in the groundwater samples from all five monitoring wells were less than the minimum reporting levels with only a few detections. Some of the organic wastewater compounds detected could have anthropogenic sources that are not biosolids. Concentrations of priority parameters in groundwater varied spatially and temporally but generally were less than Colorado regulatory limits. Concentrations of dissolved nitrate, arsenic, and selenium, in addition to chloride, sulfate, total dissolved solids, boron, iron, manganese, and uranium, in samples from some wells exceeded the Colorado standards. Concentrations of dissolved nitrate (three wells), molybdenum (one well), selenium (two wells), and uranium (one well) in shallow groundwater had significant (alpha = 0.05) upward trends in some parts of the study area. The biosolids-signature results indicate that the aquifers intercepted by the five routinely sampled wells likely have received some recharge through treated (biosolids-applied) fields or biosolids-affected ponds. Adverse effects from this biosolids-related recharge range from few (if any) at one well to large and significantly (alpha = 0.05) increasing nitrate concentrations at another well. A statistical evaluation of five paired wheat-grain samples from treated (biosolids-applied) fields and untreated (control) fields did not indicate any evidence that biosolids applications significantly (alpha = 0.05 or 0.10) increased concentration of any of these constituents in wheat grain. The wheat-grain concentrations from this study were similar to those from other studies for fields in North America where no biosolids were applied. The data for the limited crop samples indicate that biosolids applications are not increasing the concentrations of arsenic, cadmium, copper, lead, mercury, molybdenum, nickel, selenium, sulfur, and zinc in mature wheat grain from the study area.

Colorado↗

USGS2026

Groundwater-level elevations in the bedrock aquifers of the Denver Basin aquifer system, Elbert County, Colorado, 2015–23

Water users in Elbert County, Colorado, rely on groundwater from bedrock aquifers in the Denver Basin aquifer system (upper Dawson, lower Dawson, Denver, Arapahoe, and Laramie-Fox Hills aquifers) for approximately half of their water uses. Withdrawals from the bedrocks aquifers have increased to meet the water use needs of expanding regional population growth and development. The U.S. Geological Survey, in cooperation with the Elbert County Board of County Commissioners, began a study in 2015 to monitor groundwater levels within Elbert County. The primary purpose of this report is to present a summary of groundwater levels measured during the study period (2015–23) and present results from statistical analyses of changes in groundwater-level elevations through time. Discrete groundwater levels were measured at 36 wells within Elbert County. Seven of those wells contained equipment to make and record continuous groundwater-level measurements at hourly intervals. All aquifers, except the lower Dawson aquifer, had only declining groundwater-level elevations in discrete measurements for wells with statistically significant trends. Of the eight statistically significant trends in the lower Dawson aquifer, two wells indicated increasing groundwater-level elevation from discrete measurements. The groundwater-level elevation trend medians in the upper Dawson, lower Dawson, Denver, Arapahoe, and Laramie-Fox Hills aquifers were −0.23, −0.66, −0.64, −0.39, and −0.63 feet per year, respectively, for discrete groundwater-level elevation measurements. Trends in continuous groundwater-level elevations were in agreement with statistically significant trends in discrete groundwater-level elevations for all wells. The groundwater-level elevation trend medians in this study, compared to the overall trends in a 2015−2018 study, both indicated declining groundwater-level elevations except in the upper Dawson aquifer, where the trend direction was opposite, a positive trend from 2015 to 2018 and a negative trend (declining groundwater elevations) from 2015 to 2023. The change in trends within the upper Dawson aquifer may be affected by differences in the study period and the trend analysis applied. Trends during the 2015–23 study period were compared to departures from the median 2015 groundwater-level elevation for each site in each aquifer. In general, the departures from the 2015 median supported trends observed at each site and correlated spatially with greater departures near the western border of Elbert County. Additionally, 30-year precipitation data showing wet and dry periods were overlaid with the departure from the 2015 median to assess groundwater-level patterns in wells in the five aquifers. Departures from the 2015 median groundwater-level elevations appeared greatest during the dry period between 2020 and 2023. Potentiometric-surface maps of the upper and lower Dawson aquifers created from static April 2023 groundwater elevations indicated groundwater-flow direction is generally from the south to the north. Results of this study could be used to guide additional groundwater monitoring in Elbert County and could aid in long-term planning of water resources.

Colorado↗

USGS2001

Using resampling to assess reliability of audio-visual survey strategies for marbled murrelets at inland forest sites

Marbled Murrelets (Brachyramphus marmoratus) are threatened seabirds that nest in coastal old-growth coniferous forests throughout much of their breeding range. Currently, observer-based audio-visual surveys are conducted at inland forest sites during the breeding season primarily to determine nesting distribution and breeding status and are being used to estimate temporal or spatial trends in murrelet detections. Our goal was to assess the feasibility of using audio-visual survey data for such monitoring. We used an intensive field-based survey effort to record daily murrelet detections at seven survey stations in the Oregon Coast Range. We then used computer-aided resampling techniques to assess the effectiveness of twelve survey strategies with varying scheduling and a sampling intensity of 4-14 surveys per breeding season to estimate known means and SDs of murrelet detections. Most survey strategies we tested failed to provide estimates of detection means and SDs that were within A?20% of actual means and SDs. Estimates of daily detections were, however, frequently estimated to within A?50% of field data with sampling efforts of 14 days/breeding season. Additional resampling analyses with statistically generated detection data indicated that the temporal variability in detection data had a great effect on the reliability of the mean and SD estimates calculated from the twelve survey strategies, while the value of the mean had little effect. Effectiveness at estimating multi-year trends in detection data was similarly poor, indicating that audio-visual surveys might be reliably used to estimate annual declines in murrelet detections of the order of 50% per year.

Waterbirds↗

USGS2022

Data-driven prospectivity modelling of sediment-hosted Zn-Pb mineral systems and their critical raw materials

Demand for critical raw materials is expected to accelerate over the next few decades due to continued population growth and the shifting consumption patterns of the global economy. Sedimentary basins are important sources for critical raw materials and new discoveries of sediment–hosted Mississippi Valley–type (MVT) and/or clastic–dominated (CD) Zn–Pb deposits are likely required to mitigate future supply chain disruptions for Zn, Pb, Ag, Cd, Ga, Ge, Sb, and In. Herein we integrate public geoscience datasets using a discrete global grid to system to model the mineral potential for MVT and CD deposits across Canada, the United States of America, and Australia. Statistical analysis of the model results demonstrates that surface–wave tomography and derivative products from satellite gravity datasets can be used to map the most favourable paleo–tectonic settings of MVT and CD deposits inboard of orogenic belts and at the rifted edges of cratonic lithosphere, respectively. Basin development at pre–existing crustal boundaries was likely important for maintaining the low geothermal–gradients that are favourable for metal transport and generating the crustal fluid pathways that were reactivated during ore–formation, as suggested by the statistical association of both sediment–hosted mineral deposit types with the edges of upward–continued gravity and long–wavelength magnetic anomalies. Multivariate statistical analysis demonstrates that the most prospective combination of these geophysical datasets varies for each geological region and deposit type. We further demonstrate that maximum and minimum geological ages, coupled with Phanerozoic paleogeographic reconstructions, represent mappable proxies for the availability of oxidized, brine–generating regions that are the most likely source of ore–forming fluids (e.g., low– to mid–latitude carbonate platforms and evaporites). Ore deposition was likely controlled by interaction between oxidized, low–temperature brines and sulfidic and/or carbonaceous rocks, which, in some cases, can be mapped at the exposed surface or identified using the available rock descriptions. Baseline weights–of–evidence models are based on regional geophysics and are the least impacted by missing surface information but yield relatively poor results, as demonstrated by the low area–under–the–curve (AUC) for the spatially independent test set on the success–rate plot (AUC = 0.787 for MVT and AUC = 0.870 for CD). Model performance can be improved by: (1) using advanced methods that were trained and validated during a series of semi–automated machine learning competitions; and/or (2) incorporating geological and geophysical datasets that are proxies for each component of the mineral system. The best–performing gradient boosting machine models yield higher AUC for the test set (AUC = 0.983 for MVT and AUC = 0.991 for CD) and reduce the search space by >94%. The model results highlight the potential benefits of mapping sediment–hosted mineral systems at continental scale to improve mineral exploration targeting for critical raw materials.

Ore Geology Reviews↗

USGS2022

From drought to deluge: Spatiotemporal variation in migration routing, survival, travel time and floodplain use of an endangered migratory fish

We developed a novel statistical model to relate the daily survival and migration dynamics of an endangered anadromous fish to river flow and water temperature during both extreme drought and severe flooding in an intensively managed river system. Our Bayesian temporally stratified multistate mark recapture model integrates over unobserved travel times and route transitions to efficiently estimate covariate relationships and includes an adjustment for telemetry tag battery failure. We applied the model to acoustic-tagged juvenile Sacramento river winter-run Chinook salmon (Oncorhynchus tshawytscha) and found that survival decreased with decreasing river flows and increased water temperatures. We found that fish were likely to enter at a large floodplain during flood conditions and that survival in floodplain was comparable to the mainstem Sacramento river. Our study demonstrates the response of an endangered anadromous fish population to extreme spatial and temporal variability in habitat accessibility and quality. The general model framework we introduce here can be applied to telemetry of migratory fish through systems with multiple routes to efficiently estimate spatiotemporal variation in survival, travel time, and routing.

California↗

USGS2022

Estimating abundance, temporary emigration and the pattern of density dependence in a cyclic snowshoe hare population in Yukon, Canada

Estimates of demographic parameters based on capture-mark-recapture (CMR) methods may be biased when some individuals in the population are temporarily unavailable for capture (temporary emigration). We estimated snowshoe hare abundance, apparent survival, and probability of temporary emigration in a population of snowshoe hares (Lepus americanus Erxleben 1777) in the Yukon using Pollock’s robust design CMR model, and population density using spatially-explicit CMR models. Survival rates strongly varied among cyclic phases, seasons, and across five population cycles. We found strong evidence that temporary emigration was Markovian (i.e., non-random), suggesting that it varied among individuals that were temporary emigrant in the previous sampling period and those that were present in the sampled area. The probability of temporary emigration for individuals that were in the study area during the previous sampling occasion (γ´´) varied among cycles. Probability that individuals that were temporarily absent from the sampled area would remain temporary emigrants (γ´) showed strongly seasonal pattern, low in winter and high during summers. Snowshoe hare population density ranged from 0.017 (0.015–0.05) hares/ha to 4.43 (3.90–5.00) hares/ha and large-scale cyclical fluctuation. Autocorrelation functions and autoregressive analyses revealed that our study population exhibited statistically significant cyclic fluctuations, with a periodicity of 9-10 years.

Yukon↗

USGS2019

Balancing sampling intensity against spatial coverage for a community science monitoring programme

Community science is an increasingly integral part of biodiversity research and monitoring, often achieving broad spatial and temporal coverage but lower sampling intensity than studies conducted by professional scientists. When designing a community‐science monitoring programme, careful assessment of sampling designs that could be both feasible and successful at meeting programme goals is essential. Monarch butterflies ( Danaus plexippus ) are the focus of several successful community‐science projects in the U.S., but broader coverage is needed to monitor breeding areas and explain population declines observed in overwintering areas. The U.S. Monarch Conservation Science Partnership's Integrated Monarch Monitoring Program (IMMP) will representatively monitor monarchs and milkweed across North America. We performed a simulation‐based power analysis to predict trade‐offs between sampling breadth (number of sites and years) and sampling intensity (number of visits or subplots per site and year) for the IMMP. We evaluated whether each sampling design would produce sufficient statistical power to detect population trends and differences among land‐use sectors in densities of milkweed, monarch eggs, and adult monarchs. Sampling breadth had much stronger effects than sampling intensity on statistical power for all three monitoring targets. Depending on land‐use sector, monitoring 400–800 sites over 10–15 years would detect trends in densities of milkweed and adult monarchs, but no scenarios were successful for monarch eggs. Sampling 400–800 sites would also detect small (for adult monarchs) to large (for milkweed) differences among land‐use sectors in density of all three monitoring targets within the first 2–5 years. As more data become available from the IMMP, the sampling goals can be updated. Synthesis and applications . Careful sample design is an essential step in developing a successful monitoring programme. For monarchs and milkweed, we found that sampling breadth (number of sites and years) had a much stronger effect on statistical power than sampling intensity (number of visits or subsamples per site), suggesting field protocols could be tailored to maximize recruitment and retention of volunteers by minimizing the effort required to monitor each site. Many long‐term monitoring programmes might similarly benefit from evaluating the statistical trade‐offs between sampling breadth and intensity in their sampling designs.

Journal of Applied Ecology↗

USGS2009

Spatial and seasonal variations in mercury methylation and microbial community structure in a historic mercury mining area, Yolo County, California

The relationships between soil parent lithology, nutrient concentrations, microbial biomass and community structure were evaluated in soils from a small watershed impacted by historic Hg mining. Upland and wetland soils, stream sediments and tailings were collected and analyzed for nutrients (DOC, SO 4 = , NO 3 - ), Hg, MeHg, and phospholipid fatty acids (PLFA). Stream sediment was derived from serpentinite, siltstone, volcanic rocks and mineralized serpentine with cinnabar, metacinnabar and other Hg phases. Soils from different parent materials had distinct PLFA biomass and community structures that are related to nutrient concentrations and toxicity effects of trace metals including Hg. The formation of MeHg appears to be most strongly linked to soil moisture, which in turn has a correlative relationship with PLFA biomass in wetland soils. The greatest concentrations of MeHg (> 0.5??ng g - 1 MeHg) were measured in wetland soils and soil with a volcanic parent (9.5-37????g g - 1 Hg). Mercury methylation was associated with sulfate-reducing bacteria, including Desulfobacter sp. and Desulfovibrio sp., although these organisms are not exclusively responsible for Hg methylation. Statistical models of the data demonstrated that soil microbial communities varied more with soil type than with season.

Chemical Geology↗

USGS2013

A statistical analysis of the global historical volcanic fatalities record

A new database of volcanic fatalities is presented and analysed, covering the period 1600 to 2010 AD. Data are from four sources: the Smithsonian Institution, Witham (2005), CRED EM-DAT and Munich RE. The data were combined and formatted, with a weighted average fatality figure used where more than one source reports an event; the former two databases were weighted twice as strongly as the latter two. More fatal incidents are contained within our database than similar previous works; approximately 46% of the fatal incidents are listed in only one of the four sources, and fewer than 10% are in all four. 278,880 fatalities are recorded in the database, resultant from 533 fatal incidents. The fatality count is dominated by a handful of disasters, though the majority of fatal incidents have caused fewer than ten fatalities. Number and empirical probability of fatalities are broadly correlated with VEI, but are more strongly influenced by population density around volcanoes and the occurrence and extent of lahars (mudflows) and pyroclastic density currents, which have caused 50% of fatalities. Indonesia, the Philippines, and the West Indies dominate the spatial distribution of fatalities, and there is some negative correlation between regional development and number of fatalities. With the largest disasters removed, over 90% of fatalities occurred between 5 km and 30 km from volcanoes, though the most devastating eruptions impacted far beyond these distances. A new measure, the Volcano Fatality Index, is defined to explore temporal changes in societal vulnerability to volcanic hazards. The measure incorporates population growth and recording improvements with the fatality data, and shows prima facie evidence that vulnerability to volcanic hazards has fallen during the last two centuries. Results and interpretations are limited in scope by the underlying fatalities data, which are affected by under-recording, uncertainty, and bias. Attempts have been made to estimate the extent of these issues, and to remove their effects where possible. The data analysed here are provided as supplementary material. An updated version of the Smithsonian fatality database fully integrated with this database will be publicly available in the near future and subsequently incorporate new data.

Journal of Applied Volcanology↗

USGS2002

Spatial and temporal distribution of contaminated, effluent-affected sediment on the Palos Verdes margin, southern California

A sedimentary deposit on the continental margin near the Palos Verdes Peninsula, California is comprised of sewage effluent and geologic materials and is contaminated with metals, pesticides (including DDT and associated compounds), and PCBs. The deposit was mapped with subbottom acoustic profilers, and sediment cores were analyzed for geochemical and physical properties to determine the volume of the deposit and the distribution and mass of contaminants. Mapping showed that the deposit ranges up to 60-cm thick, has a total volume exceeding 9 million m3, and covers over 40 km2. Virtually the entire effluent-affected deposit is contaminated with DDT and PCBs. Nearly half of the area of the deposit lies on the continental slope, but 70-75% of the volume of the deposit and total mass of DDT reside on the continental shelf. Analysis of data collected biennially since 1981 by the Sanitation Districts of Los Angeles County show that the mass of DDT has apparently decreased at some stations but has remained essentially constant at others. Temporal changes m mass per unit area of DDT are not statistically significant (at the 90% confidence level) at the most contaminated locations over a 16-yr period. The results of this mapping effort were used as a basis for modeling efforts described elsewhere in this issue.

California↗

USGS2011

Variation in local abundance and species richness of stream fishes in relation to dispersal barriers: Implications for management and conservation

1.Barriers to immigration, all else being equal, should in principle depress local abundance and reduce local species richness. These issues are particularly relevant to stream-dwelling species when improperly designed road crossings act as barriers to migration with potential impacts on the viability of upstream populations. However, because abundance and richness are highly spatially and temporally heterogeneous and the relative importance of immigration on demography is uncertain, population- and community-level effects can be difficult to detect. 2.In this study, we tested the effects of potential barriers to upstream movements on the local abundance and species richness of a diverse assemblage of resident stream fishes in the Monongahela National Forest, West Virginia, U.S.A. Fishes were sampled using simple standard techniques above- and below road crossings that were either likely or unlikely to be barriers to upstream fish movements (based on physical dimensions of the crossing). We predicted that abundance of resident fishes would be lower in the upstream sections of streams with predicted impassable barriers, that the strength of the effect would vary among species and that variable effects on abundance would translate into lower species richness. 3.Supporting these predictions, the statistical model that best accounted for variation in abundance and species richness included a significant interaction between location (upstream or downstream of crossing) and type (passable or impassable crossing). Stream sections located above predicated impassable culverts had fewer than half the number of species and less than half the total fish abundance, while stream sections above and below passable culverts had essentially equivalent richness and abundance. 4.Our results are consistent with the importance of immigration and population connectivity to local abundance and species richness of stream fishes. In turn, these results suggest that when measured at appropriate scales (multiple streams within catchments), with simple protocols amenable to use by management agencies, differences in local abundance and species richness may serve as indicators of the extent to which road crossings are barriers to fish movement and help determine whether road-crossing improvements have restored connectivity to stream fish populations and communities. Published 2011. This article is a US Government work and is in the public domain in the USA.

Freshwater Biology↗

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↗