USGS Science⌕ Search

SEARCH · USGS Science

Results for “Spatial Statistics”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 901 records · Page 50Linked to original sources
USGS2022

Large-scale distribution models for optimal prediction of Eastern black rail habitat within tidal ecosystems

Eastern black rails ( Laterallus jamaicensis jamaicensis ) are among the rarest and least-studied birds in North America and were recently listed as threatened under the U.S. Endangered Species Act. Spatial models that predict habitat quality across the subspecies range are therefore needed to inform conservation, recovery, and monitoring efforts for this rare bird. We used data from 47,585 call-broadcast surveys collected at 7906 sites over a 3-decade period (1990s, 2000s, 2010s; 23 total years) to build species distribution models for eastern black rails. We used hierarchical Bayesian occupancy models and predictive model selection to develop multi-scale models that optimally predict habitat suitability for eastern black rails within tidal wetlands while also accounting for imperfect detection of these cryptic birds during field surveys. We also used raster regression techniques to translate model predictions into 30-m resolution maps of habitat suitability for eastern black rails within tidal wetlands along the eastern seaboard of the United States. The model predicted suitability of breeding habitat as a function of wetland attributes (e.g., cover of high marsh and terrestrial border), hydrologic modification, and disturbance from human development measured over multiple spatial scales. We also found differences in habitat relationships for eastern black rails when compared to models that included both North American subspecies of black rail. Important results included negative effects of shrub-scrub wetlands, and strong positive effects of high marsh, terrestrial border, and impoundments on breeding season occupancy. Our study provides an example of integrating detection-non-detection data and modern statistical methods to build predictive distribution models for an extremely rare species, while also providing rigorous predictions of breeding habitat quality for the eastern black rail within tidal wetlands. These models will facilitate optimal monitoring, habitat conservation, and recovery planning efforts for eastern black rails and provide a foundation for future research and conservation of this imperiled bird.

Global Ecology and Conservation↗
USGS2010

Geomorphic Framework to assess changes to aquatic habitat due to flow regulation and channel and floodplain alteration, Cedar River, Washington

Flow regulation, bank armoring, and floodplain alteration since the early 20th century have contributed to significant changes in the hydrologic regime and geomorphic processes of the Cedar River in Washington State. The Cedar River originates in the Cascade Range, provides drinking water to the Seattle metropolitan area, and supports several populations of anadromous salmonids. Flow regulation currently has limited influence on the magnitude, duration, and timing of high-flow events, which affect the incubation of salmonids as well as the production and maintenance of their habitat. Unlike structural changes to the channel and floodplain, flow regulation may be modified in the short-term to improve the viability of salmon populations. An understanding of the effects of flow regulation on those populations must be discerned over a range of scales from individual floods that affect the size of individual year classes to decadal high flow regime that influences the amount and quality of channel and off-channel habitat available for spawning and rearing. We present estimates of reach-scale sediment budgets and changes to channel morphology derived from historical orthoimagery, specific gage analyses at four long-term streamflow-gaging stations to quantify trends in aggradation, and hydrologic statistics of the magnitude and duration of peak streamflows. These data suggest a gradient of channel types from unconfined, sediment-rich segments to confined, sediment-poor segments that are likely to have distinct responses to high flows. Particle-size distribution data and longitudinal water surface and streambed profiles for the 56 km downstream of Chester Morse Lake measured in 2010 show the spatial extent of preferred salmonid habitat along the Cedar River. These historical and current data constitute a geomorphic framework to help assess different river management scenarios for salmonid habitat and population viability. PDF version of a presentation on changes to aquatic habitat at the Cedar River in Washington state. Presented at the American Geophysical Union Fall Meeting 2010.

Washington↗
USGS2016

Environmental drivers of differences in microbial community structure in crude oil reservoirs across a methanogenic gradient

Stimulating in situ microbial communities in oil reservoirs to produce natural gas is a potentially viable strategy for recovering additional fossil fuel resources following traditional recovery operations. Little is known about what geochemical parameters drive microbial population dynamics in biodegraded, methanogenic oil reservoirs. We investigated if microbial community structure was significantly impacted by the extent of crude oil biodegradation, extent of biogenic methane production, and formation water chemistry. Twenty-two oil production wells from north central Louisiana, USA, were sampled for analysis of microbial community structure and fluid geochemistry. Archaea were the dominant microbial community in the majority of the wells sampled. Methanogens, including hydrogenotrophic and methylotrophic organisms, were numerically dominant in every well, accounting for, on average, over 98% of the total Archaea present. The dominant Bacteria groups were Pseudomonas, Acinetobacter , Enterobacteriaceae, and Clostridiales, which have also been identified in other microbially-altered oil reservoirs. Comparing microbial community structure to fluid (gas, water, and oil) geochemistry revealed that the relative extent of biodegradation, salinity, and spatial location were the major drivers of microbial diversity. Archaeal relative abundance was independent of the extent of methanogenesis, but closely correlated to the extent of crude oil biodegradation; therefore, microbial community structure is likely not a good sole predictor of methanogenic activity, but may predict the extent of crude oil biodegradation. However, when the shallow, highly biodegraded, low salinity wells were excluded from the statistical analysis, no environmental parameters could explain the differences in microbial community structure. This suggests that the microbial community structure of the 5 shallow, up-dip wells was different than the 17 deeper, down-dip wells. Also, the 17 down-dip wells had statistically similar microbial communities despite significant changes in environmental parameters between oil fields. Together, this implies that no single microbial population is a reliable indicator of a reservoir's ability to degrade crude oil to methane, and that geochemistry may be a more important indicator for selecting a reservoir suitable for microbial enhancement of natural gas generation.

Louisiana↗
USGS2024

Increasing seasonal variation in the extent of rivers and lakes from 1984 to 2022

Knowledge of the spatial and temporal distribution of surface water is important for water resource management, flood risk assessment, monitoring ecosystem health, constraining estimates of biogeochemical cycles and understanding our climate. While global-scale spatiotemporal change detection of surface water has significantly improved in recent years due to planetary-scale remote sensing and computing, it has remained challenging to distinguish the changing characteristics of rivers and lakes. Here we analyze the spatial extent of permanent and seasonal rivers and lakes globally over the past 38 years based on new data of river system extents and surface water trends. Results show that while the total permanent surface area of both rivers and lakes has remained relatively constant, the areas with intermittent seasonal coverage have increased by 12 % and 27 % for rivers and lakes, respectively. The increase is statistically significant in over 84 % of global water catchments based on Spearman's rank correlations (rho) above 0.05 and p values less than 0.05. The seasonal river extent is nearly 32 % larger than the previously observed annual mean river extent, suggesting large seasonal variations that impact not only ecosystem health but also estimations of terrestrial biogeochemical cycles of carbon. The outcomes of our analysis are shared as the Surface Area of Rivers and Lakes (SARL) database, serving as a valuable resource for monitoring and research of hydrological cycles, ecosystem accounting, and water management.

Hydrology and Earth System Sciences↗
USGS2008

Data to support statistical modeling of instream nutrient load based on watershed attributes, southeastern United States, 2002

This report presents and describes the digital datasets that characterize nutrient source inputs, environmental characteristics, and instream nutrient loads for the purpose of calibrating and applying a nutrient water-quality model for the southeastern United States for 2002. The model area includes all of the river basins draining to the south Atlantic and the eastern Gulf of Mexico, as well as the Tennessee River basin (referred to collectively as the SAGT area). The water-quality model SPARROW (SPAtially-Referenced Regression On Watershed attributes), developed by the U.S. Geological Survey, uses a regression equation to describe the relation between watershed attributes (predictors) and measured instream loads (response). Watershed attributes that are considered to describe nutrient input conditions and are tested in the SPARROW model for the SAGT area as source variables include atmospheric deposition, fertilizer application to farmland, manure from livestock production, permitted wastewater discharge, and land cover. Watershed and channel attributes that are considered to affect rates of nutrient transport from land to water and are tested in the SAGT SPARROW model as nutrient-transport variables include characteristics of soil, landform, climate, reach time of travel, and reservoir hydraulic loading. Datasets with estimates of each of these attributes for each individual reach or catchment in the reach-catchment network are presented in this report, along with descriptions of methods used to produce them. Measurements of nutrient water quality at stream monitoring sites from a combination of monitoring programs were used to develop observations of the response variable - mean annual nitrogen or phosphorus load - in the SPARROW regression equation. Instream load of nitrogen and phosphorus was estimated using bias-corrected log-linear regression models using the program Fluxmaster, which provides temporally detrended estimates of long-term mean load well-suited for spatial comparisons. The detrended, or normalized, estimates of load are useful for regional-scale assessments but should be used with caution for local-scale interpretations, for which use of loads estimated for actual time periods and employing more detailed regression analysis is suggested. The mean value of the nitrogen yield estimates, normalized to 2002, for 637 stations in the SAGT area is 4.7 kilograms per hectare; the mean value of nitrogen flow-weighted mean concentration is 1.2 milligrams per liter. The mean value of the phosphorus yield estimates, normalized to 2002, for the 747 stations in the SAGT area is 0.66 kilogram per hectare; the mean value of phosphorus flow-weighted mean concentration is 0.17 milligram per liter. Nutrient conditions measured in streams affected by substantial influx or outflux of water and nutrient mass across surface-water basin divides do not reflect nutrient source and transport conditions in the topographic watershed; therefore, inclusion of such streams in the SPARROW modeling approach is considered inappropriate. River basins identified with this concern include south Florida (where surface-water flow paths have been extensively altered) and the Oklawaha, Crystal, Lower Sante Fe, Lower Suwanee, St. Marks, and Chipola River basins in central and northern Florida (where flow exchange with the underlying regional aquifer may represent substantial nitrogen influx to and outflux from the surface-water basins).

Open-File Report↗
USGS2026

ENSO and PDO drive shoreline position anomalies in the U.S. Pacific Northwest

Sandy beaches act as buffers against various coastal hazards but are vulnerable to episodic (seasonal) and chronic (interannual) erosion. Understanding the variation of shoreline position, a key metric in coastal morphology, over a spectrum of time scales is therefore crucial in assessing hazard vulnerability. Long-standing research has investigated the role of El Niño-Southern Oscillation (ENSO), the dominant mode of climate variability in the Pacific Basin, in seasonal shoreline variability. Yet, ENSO’s chronic influence—and that of another Pacific climate mode, the Pacific Decadal Oscillation (PDO)—on shoreline anomalies remains poorly understood. Here, we examine the variability of sandy beaches in the US Pacific Northwest, a ∼750 km long coastal region on the US West Coast. We leverage 40 years (1984–2024) of shoreline data from publicly available Earth-observing (Landsat) satellite imagery at a high spatial resolution (>10,000 shore-normal transects at 50-m alongshore spacing) and employ Convergent Cross Mapping (CCM), a methodology for inferring causality in dynamical systems. We discover that strong El Niño years are signified by erosion (75.1% of transects), and strong La Niña years exhibit accretional behavior (73.4% of transects). Furthermore, we establish, for the first time, that both ENSO and PDO exert a statistically significant control on interannual shoreline variability, particularly on the alongshore component (in 95 and 100% of littoral cells, respectively), with water level fluctuations playing a critical role. This effort advances our understanding of the seasonal-to-interannual interactions between Pacific Basin climate variability and the PNW’s coastal morphodynamics, with implications for sediment management and coastal adaptation.

California, Oregon, Washington↗
USGS2017

Translating statistical species-habitat models to interactive decision support tools

Understanding species-habitat relationships is vital to successful conservation, but the tools used to communicate species-habitat relationships are often poorly suited to the information needs of conservation practitioners. Here we present a novel method for translating a statistical species-habitat model, a regression analysis relating ring-necked pheasant abundance to landcover, into an interactive online tool. The Pheasant Habitat Simulator combines the analytical power of the R programming environment with the user-friendly Shiny web interface to create an online platform in which wildlife professionals can explore the effects of variation in local landcover on relative pheasant habitat suitability within spatial scales relevant to individual wildlife managers. Our tool allows users to virtually manipulate the landcover composition of a simulated space to explore how changes in landcover may affect pheasant relative habitat suitability, and guides users through the economic tradeoffs of landscape changes. We offer suggestions for development of similar interactive applications and demonstrate their potential as innovative science delivery tools for diverse professional and public audiences.

PLoS ONE↗
USGS2020

Potential duration of aftershocks of the 2020 southwestern Puerto Rico earthquake

Abstract Aftershocks (earthquakes clustered spatially and chronologically near the occurrence of a causative earthquake) are ongoing in southwestern Puerto Rico after a series of earthquakes, which include a magnitude 6.4 earthquake that occurred near Barrio Indios, Guayanilla, on January 7, 2020, and affected the surrounding area. This report estimates the expected duration of these aftershocks by incorporating observations of aftershocks as of January 17, 2020, into a well-established statistical model of how earthquake sequences behave. Aftershocks will persist for years to decades, although with decreasing frequency, and earthquakes will likely be felt on a daily basis for up to several months. These estimates have significant uncertainty owing to different scenarios of how the earthquake sequence may evolve over time and could also change if a new large aftershock occurs. This report also estimates the amount of time remaining until the annual probability of magnitude 5, 6, and 7 or greater aftershocks—which could cause additional damage—decreases to 50, 25, 10, 5, and 1 percent. As of this writing, the chance of having a magnitude 6 or greater earthquake within a given year, going forward, will not fall below 25 percent for another 3 months to 3 years. The chance of having a magnitude 5 or greater earthquake will not fall below 25 percent for a decade or more. The aftershocks discussed in this report would be located in the same general area as the aftershocks that have already occurred. Our results do not imply a change in the risk of earthquakes in other parts of Puerto Rico. Resumen Las réplicas (terremotos agrupados espacial y cronológicamente cerca de la ocurrencia de un terremoto causante) están en curso en el suroeste de Puerto Rico después de una serie de terremotos, que incluyen un terremoto de magnitud 6.4, ocurrido cerca del Barrio Indios, Guayanilla, el 7 de enero de 2020 y que afectaron las áreas circundantes. Este informe estima la duración esperada de las réplicas incorporando observaciones de réplicas a partir del 17 de enero de 2020 en un modelo estadístico bien establecido de cómo se comportan las secuencias de terremotos. Las réplicas persistirán durante años o décadas, aunque con una frecuencia decreciente, y los terremotos probablemente se sentirán a diario durante varios meses. Estos estimados tienen una incertidumbre significativa debido a diferentes escenarios de cómo la secuencia del terremoto puede evolucionar con el tiempo y también podrían cambiar si ocurre una nueva réplica grande. Este informe también estima la cantidad de tiempo restante hasta que la probabilidad anual de réplicas de magnitud 5, 6 y 7 o más - que podría causar daños adicionales- disminuya a un 50, 25, 10, 5 y 1 por ciento. Al momento de escribir este artículo, la posibilidad de tener un terremoto de magnitud 6 en un año determinado, en el futuro, no caerá por debajo del 25 por ciento durante otros 3 meses a 3 años. La probabilidad de tener un terremoto de magnitud 5 o mayor no será inferior al 25 por ciento durante una década o más. Las réplicas discutidas en este informe se ubicarían en la misma área general que las réplicas que ya han ocurrido. Nuestros resultados no implican un cambio en el riesgo de terremotos en otras partes de Puerto Rico.

Puerto Rico↗
USGS1994

Land use in, and water quality of, the Pea Hill Arm of Lake Gaston, Virginia and North Carolina, 1988-90

The City of Virginia Beach currently (1994) supplies water to about 400,000 people in southeastern Virginia. The city plans to withdraw water from the Pea Hill Arm of Lake Gaston to meet projected water needs of the population to the year 2030. The purpose of this report is to (1) describe the temporal and spatial distribution of selected water-quality constituents, (2) document current (1989) land use and land cover in the Pea Hill Arm drainage basin, and (3) discuss relations, if any, between the quality of water in the inlets within the Pea Hill Arm and land uses. The report focuses on water-quality problems in the basin, including changes in concentrations of major ions, nutrients, and algae associated with urban development adjacent to water bodies. The Pea Hill Arm was classified as mesotrophic on the basis of the range of concentrations of total phosphorus (0.001 to 0.61 milligrams per liter); the range of concentrations of total organic-plus-ammonia nitrogen (0.2 to 1.4 milligrams per liter); and the range of concentrations of chlorophyll a (1.4 to 56 micrograms per liter). These water-quality data were collected at 3 feet below the water surface during water years 1989-90. Thermal stratification in Pea Hill Arm generally began in April and ended in September. Water below a depth of about 25 feet generally became anoxic by June. Destratification generally began in late September and was completed by November. Lake Gaston followed the same general stratification and destratification pattern as Pea Hill Arm, except Lake Gaston was partially destratified during the summer when large amounts of water were released from John H. Kerr Reservoir and Lake Gaston Dams. During water year 1988, streamflows were 33 percent below the long-term mean-annual streamflows at one of the major streams to Lake Gaston. Low streamflows contributed to elevated specific conductances and concentrations of sodium, calcium, magnesium, and alkalinity from October 1988 to February 1989 at sampling stations in the Pea Hill Arm and Lake Gaston. About 75 percent of the land use in the Pea Hill Arm is forest land. The remaining 25 percent of the Pea Hill Arm drainage basin is 8 percent pasture/open land, 8 percent open water, 6 percent residential land, and 3 percent cropland. No statistical relations are present between water-quality constituents measured and developed land uses within 11 basins in the Pea Hill Arm Basin, except during periods of stormwater runoff. During a stormwater-runoff event, there was a relation between total nitrite plus nitrate and land use (Kendall's tau correlation coefficient of 0.69). The relation between the developed land use and total nitrite plus nitrate can also be related to the increased ground-water inputs during high base-flow periods. Spatial differences in water-quality constituents as determined by Wilcoxon (matched-pairs) signed-rank tests and cluster analyses were longitudinal and primarily grouped into riverine, transition, and lacustrine zones. These zones were grouped on the basis of flow characteristics and nutrient concentrations.

North Carolina, Virginia↗
USGS2009

Shifts in the trophic base of intermittent stream food webs

Understanding spatial and temporal variation in the trophic base of stream food webs is critical for predicting population and community stability, and ecosystem function. We used stable isotope ratios ( 13 C/ 12 C, and 15 N/ 14 N) to characterize the trophic base of two streams in the Ozark Mountains of northwest Arkansas, U.S.A. We predicted that autochthonous resources would be more important during the spring and summer and allochthonous resources would be more important in the winter due to increased detritus inputs from the riparian zone during autumn leaf drop. We predicted that stream communities would demonstrate increased reliance on autochthonous resources at sites with larger watersheds and greater canopy openness. The study was conducted at three low-order sites in the Mulberry River Drainage (watershed area range: 81-232 km 2 ) seasonally in 2006 and 2007. We used circular statistics to examine community-wide shifts in isotope space among fish and invertebrate consumers in relation to basal resources, including detritus and periphyton. Mixing models were used to quantify the relative contribution of autochthonous and allochthonous energy sources to individual invertebrate consumers. Significant isotopic shifts occurred but results varied by season and site indicating substantial variation in the trophic base of stream food webs. In terms of temporal variation, consumers shifted toward periphyton in the summer during periods of low discharge, but results varied during the interval between summer and winter. Our results did not demonstrate increased reliance on periphyton with increasing watershed area or canopy openness, and detritus was important at all the sites. In our study, riffle-pool geomorphology likely disrupted the expected spatial pattern and stream drying likely impacted the availability and distribution of basal resources.

Hydrobiologia↗
USGS2001

Bird-landscape relations in the Chihuahuan Desert: Coping with uncertainties about predictive models

During the springs of 1995–1997, we studied birds and landscapes in the Chihuahuan Desert along part of the Texas–Mexico border. Our objectives were to assess bird–landscape relations and their interannual consistency and to identify ways to cope with associated uncertainties that undermine confidence in using such relations in conservation decision processes. Bird distributions were often significantly associated with landscape features, and many bird–landscape models were valid and useful for predictive purposes. Differences in early spring rainfall appeared to influence bird abundance, but there was no evidence that annual differences in bird abundance affected model consistency. Model consistency for richness (42%) was higher than mean model consistency for 26 focal species (mean 30%, range 0–67%), suggesting that relations involving individual species are, on average, more subject to factors that cause variation than are richness–landscape relations. Consistency of bird–landscape relations may be influenced by such factors as plant succession, exotic species invasion, bird species' tolerances for environmental variation, habitat occupancy patterns, and variation in food density or weather. The low model consistency that we observed for most species indicates the high variation in bird–landscape relations that managers and other decision makers may encounter. The uncertainty of interannual variation in bird–landscape relations can be reduced by using projections of bird distributions from different annual models to determine the likely range of temporal and spatial variation in a species' distribution. Stochastic simulation models can be used to incorporate the uncertainty of random environmental variation into predictions of bird distributions based on bird–landscape relations and to provide probabilistic projections with which managers can weigh the costs and benefits of various decisions. Uncertainty about the true structure of bird–landscape relations (structural uncertainty) can be reduced by ensuring that models meet important statistical assumptions, designing studies with sufficient statistical power, validating the predictive ability of models, and improving model accuracy through continued field sampling and model fitting. Uncertainty associated with sampling variation (partial observability) can be reduced by ensuring that sample sizes are large enough to provide precise estimates of both bird and landscape parameters. By decreasing the uncertainty due to partial observability, managers will improve their ability to reduce structural uncertainty.

Chihuahua, Coahuila, Nuevo León, Tamaulipas, Texas↗
USGS2016

Mortality monitoring design for utility-scale solar power facilities

Introduction Solar power represents an important and rapidly expanding component of the renewable energy portfolio of the United States (Lovich and Ennen, 2011; Hernandez and others, 2014). Understanding the impacts of renewable energy development on wildlife is a priority for the U.S. Fish and Wildlife Service (FWS) in compliance with Department of Interior Order No. 3285 (U.S. Department of the Interior, 2009) to “develop best management practices for renewable energy and transmission projects on the public lands to ensure the most environmentally responsible development and delivery of renewable energy.” Recent studies examining effects of renewable energy development on mortality of migratory birds have primarily focused on wind energy (California Energy Commission and California Department of Fish and Game, 2007), and in 2012 the FWS published guidance for addressing wildlife conservation concerns at all stages of land-based wind energy development (U.S. Fish and Wildlife Service, 2012). As yet, no similar guidelines exist for solar development, and no published studies have directly addressed the methodology needed to accurately estimate mortality of birds and bats at solar facilities. In the absence of such guidelines, ad hoc methodologies applied to solar energy projects may lead to estimates of wildlife mortality rates that are insufficiently accurate and precise to meaningfully inform conversations regarding unintended consequences of this energy source and management decisions to mitigate impacts. Although significant advances in monitoring protocols for wind facilities have been made in recent years, there remains a need to provide consistent guidance and study design to quantify mortality of bats, and resident and migrating birds at solar power facilities (Walston and others, 2015). In this document, we suggest methods for mortality monitoring at solar facilities that are based on current methods used at wind power facilities but adapted for the unique conditions encountered at solar facilities. In particular, unlike at wind-power facilities, the unimpeded access to almost all areas within the facilities, the typically flat terrain, and general absence of thick vegetation allow distance-sampling techniques (Buckland and others, 2001, 2004) to be exploited to advantage at industrial solar sites. These protocols build on the work of Nicolai and others (2011), and as our understanding and techniques for monitoring improve, the methods may be further modified to incorporate improvements in the future. We present case studies based on monitoring methods currently implemented at different utility-scale solar facilities to illustrate how distance-sampling techniques may improve overall detectability without substantially increasing costs. Every facility is unique, and the protocols presented may be adapted based on specific monitoring objectives and conditions at each site. We provide guidance for designing monitoring programs whose objective it is to estimate the total number of bird and bat fatalities occurring at a facility over an extended period of time. We address spatial variation in causes of mortality, as well as potential sources of imperfect detection, for example, animals falling in or moving to unsearched areas, carcasses removed by predators, and carcasses missed by searchers. We suggest methods to estimate and account for each source of imperfect detection. This document focuses on monitoring design only and does not discuss approaches for estimating mortality from collected data. The development of statistically sound estimators relevant to the solar context is a current topic of research, although there are already strong foundations for estimation with distance-sampling methods in similar open, arid environments (Anderson and others, 2001; Freilich and others, 2005). Nonetheless, if protocols described in this document are followed, the resulting data will be adequate and sufficient for estimating mortality using newly formulated estimators.

Open-File Report↗
USGS2024

Peak streamflow trends in South Dakota and their relation to changes in climate, water years 1921–2020

Peak-flow (flood) frequency analysis is essential to water-resources management applications, including the design of critical infrastructure such as bridges and culverts, and floodplain mapping. Federal guidelines for performing peak-flow flood frequency analyses are presented in a U.S. Geological Survey Techniques and Methods Report known as Bulletin 17C. A basic assumption within Bulletin 17C, which documents the guidelines for determining annual peak streamflow frequency, is that, for basins without major hydrologic alterations (for example, regulation, diversion, and urbanization), statistical properties of the distribution of annual peak streamflows are stationary; that is, the mean, variance, and skew are constant through time. Nonstationarity is a statistical property of a peak-flow series such that the long-term (on the order of decades) distributional properties change one or more times either gradually or abruptly through time. Individual nonstationarities may be attributed to one source such as flow regulation, land-use change, or climate but are often the result of a combination of sources, making detection and attribution of nonstationarities challenging. In response to a growing concern regarding nonstationarity in peak streamflows in the region, the U.S. Geological Survey, in cooperation with the Departments of Transportation of Illinois, Iowa, Michigan, Minnesota, Missouri, South Dakota, and Wisconsin; the Montana Department of Natural Resources and Conservation; and the North Dakota Department of Water Resources, assessed the potential nonstationarity in peak streamflows in the north-central United States. This chapter characterizes the effects of natural hydroclimatic shifts and potential climate change on annual peak streamflows in the State of South Dakota. Annual peak and daily streamflow as well as model-simulated gridded climatic data were examined for temporal monotonic trends, change points, and other statistical properties indicative of changing climatic and environmental conditions. Changes in annual peak and daily flows were evaluated among 13, 35, and 81 qualifying U.S. Geological Survey streamgages for the 75-, 50-, and 30-year trend periods through water year 2020 (the period from October 1, 2019, to September 30, 2020) in South Dakota, respectively. No qualifying streamgages were in the 100-year trend period in the State. Statistical tests for autocorrelation (independent and identically distributed assumption), monotonic trends, and change points in the median and scale are analyzed to evaluate potential stationarity violations (nonstationarity) for performing at-site peak-flow flood-frequency analysis. The trends are reported using a likelihood approach as an alternative to simply reporting significant trends with an arbitrary p -value cutoff point. A distinct east-west spatial pattern of likely upward and downward monotonic trends and change points, respectively, was detected in 75- and 50-year trend periods, but an inconsistent spatial pattern was detected in the 30-year trend period. Additionally, change points in the median annual peak streamflows were detected in the late 1970s and early 1980s in the western part of the State, but in the east, the change point was more commonly detected in 1992–93. A similar east-west spatial pattern of likely upward and downward trends was detected in the annual peak-flow timing, the day of the year of the annal peak streamflow. In the western part of the State, the annual peak streamflows are arriving earlier, but in the east, the annual peak streamflows are arriving later. A peaks-over-threshold (POT) analysis where, on average, there are two events per year (POT2) and four events per year (POT4) was also used to evaluate changes in the frequency (count) of daily streamflows exceeding the threshold. Similar to detected changes in the annual peak streamflow, an east-west likely upward or downward change corresponding to an increase or decrease, respectively, in the frequency of daily streamflow greater than a POT2 and POT4 threshold was detected. A monthly water-balance model was used to evaluate hydroclimatic variation in annual and seasonal precipitation, snowfall, potential evapotranspiration, and soil moisture storage for all qualifying streamgages in the 75-, 50-, and 30-year trend periods. Detected trends in the annual hydroclimatic metrics for the 75- and 50-year trend periods indicate a spatially consistent statewide increase in precipitation, decrease in snowfall, increase in potential evapotranspiration, and increase in soil moisture storage. Furthermore, detected trends in seasonal precipitation in the 75- and 50-year trend periods highlight a pronounced change in precipitation in winter and later into the summer season, especially in the 50-year trend period in the eastern part of the State. Statewide increases in seasonal soil moisture storage were also detected, highlighting year-round increasing flood magnitudes, particularly in the eastern part of the State. Based on the results of these stationarity tests for the qualifying streamgages in South Dakota among the 75-, 50-, and 30-year trend periods, consistent temporal and spatial patterns of nonstationarity were detected among the 75- and 50-year trend periods. Furthermore, when nonstationarity is detected in daily streamflow, increased streamflow and volume (increasing frequency in POT), as well as potentially bridge scour, may have implications on culvert and highway design in the eastern part of South Dakota. Thus, when performing at-site peak-flow flood-frequency analyses in South Dakota, potential nonstationarities and alternative approaches are important considerations.

South Dakota↗
USGS2020

Social Values for Ecosystem Services, version 4.0 (SolVES 4.0)—Documentation and user manual

The geographic information system tool, Social Values for Ecosystem Services (SolVES), was developed to incorporate quantified and spatially explicit measures of social values into ecosystem service assessments. SolVES 4.0 provides an open-source version of SolVES, which was designed to assess, map, and quantify the social values of ecosystem services. Social values—the perceived, nonmarket values the public ascribes to ecosystem services, particularly cultural services, such as aesthetics and recreation—can be evaluated for various stakeholder groups. These groups are distinguishable by factors such as their attitudes and preferences regarding public uses (for example, motorized recreation and logging). As with previous versions, SolVES 4.0 derives a quantitative 10-point, social-values metric—the value index—from a combination of spatial and nonspatial responses to public value and preference surveys. The tool also calculates metrics characterizing the underlying environment, such as average distance to water and dominant landcover. SolVES 4.0 has been developed with Python using a QGIS user interface and a PostgreSQL database for required data. SolVES is integrated with Maxent maximum entropy modeling software to generate more complete social-value maps and offer robust statistical models describing the relation between the value index and explanatory environmental variables. A model’s goodness of fit to a primary study area and its potential performance in transferring social values to similar areas using value-transfer methods can be evaluated. SolVES 4.0 provides an improved open-source, public-domain tool for decision makers and researchers to evaluate the social values of ecosystem services and to facilitate discussions among diverse stakeholders regarding the tradeoffs among ecosystem services in a variety of biophysical and social contexts including mountain, forest, coastal, riparian, agricultural, and urban environments around the globe.

Techniques and Methods↗
USGS2011

Improving national-scale invasion maps: Tamarisk in the western United States

New invasions, better field data, and novel spatial-modeling techniques often drive the need to revisit previous maps and models of invasive species. Such is the case with the at least 10 species of Tamarix , which are invading riparian systems in the western United States and expanding their range throughout North America. In 2006, we developed a National Tamarisk Map by using a compilation of presence and absence locations with remotely sensed data and statistical modeling techniques. Since the publication of that work, our database of Tamarix distributions has grown significantly. Using the updated database of species occurrence, new predictor variables, and the maximum entropy (Maxent) model, we have revised our potential Tamarix distribution map for the western United States. Distance-to-water was the strongest predictor in the model (58.1%), while mean temperature of the warmest quarter was the second best predictor (18.4%). Model validation, averaged from 25 model iterations, indicated that our analysis had strong predictive performance (AUC = 0.93) and that the extent of Tamarix distributions is much greater than previously thought. The southwestern United States had the greatest suitable habitat, and this result differed from the 2006 model. Our work highlights the utility of iterative modeling for invasive species habitat modeling as new information becomes available.

Western United States↗
USGS2003

Incorporating precision, accuracy and alternative sampling designs into a continental monitoring program for colonial waterbirds

A comprehensive monitoring program for colonial waterbirds in North America has never existed. At smaller geographic scales, many states and provinces conduct surveys of colonial waterbird populations. Periodic regional surveys are conducted at varying times during the breeding season using a variety of survey methods, which complicates attempts to estimate population trends for most species. The US Geological Survey Patuxent Wildlife Research Center has recently started to coordinate colonial waterbird monitoring efforts throughout North America. A centralized database has been developed with an Internet-based data entry and retrieval page. The extent of existing colonial waterbird surveys has been defined, allowing gaps in coverage to be identified and basic inventories completed where desirable. To enable analyses of comparable data at regional or larger geographic scales, sampling populations through statistically sound sampling designs should supersede obtaining counts at every colony. Standardized breeding season survey techniques have been agreed upon and documented in a monitoring manual. Each survey in the manual has associated with it recommendations for bias estimation, and includes specific instructions on measuring detectability. The methods proposed in the manual are for developing reliable, comparable indices of population size to establish trend information at multiple spatial and temporal scales, but they will not result in robust estimates of total population numbers.

Ornis Hungarica↗
USGS2014

Monitoring conterminous United States (CONUS) land cover change with Web-Enabled Landsat Data (WELD)

Forest cover loss and bare ground gain from 2006 to 2010 for the conterminous United States (CONUS) were quantified at a 30 m spatial resolution using Web-Enabled Landsat Data available from the USGS Center for Earth Resources Observation and Science (EROS) (http://landsat.usgs.gov/WELD.php). The approach related multi-temporal WELD metrics and expert-derived training data for forest cover loss and bare ground gain through a decision tree classification algorithm. Forest cover loss was reported at state and ecoregional scales, and the identification of core forests' absent of change was made and verified using LiDAR data from the GLAS (Geoscience Laser Altimetry System) instrument. Bare ground gain correlated with population change for large metropolitan statistical areas (MSAs) outside of desert or semi-desert environments. GoogleEarth™ time-series images were used to validate the products. Mapped forest cover loss totaled 53,084 km2 and was found to be depicted conservatively, with a user's accuracy of 78% and a producer's accuracy of 68%. Excluding errors of adjacency, user's and producer's accuracies rose to 93% and 89%, respectively. Mapped bare ground gain equaled 5974 km2 and nearly matched the estimated area from the reference (GoogleEarth™) classification; however, user's (42%) and producer's (49%) accuracies were much less than those of the forest cover loss product. Excluding errors of adjacency, user's and producer's accuracies rose to 62% and 75%, respectively. Compared to recent 2001–2006 USGS National Land Cover Database validation data for forest loss (82% and 30% for respective user's and producer's accuracies) and urban gain (72% and 18% for respective user's and producer's accuracies), results using a single CONUS-scale model with WELD data are promising and point to the potential for national-scale operational mapping of key land cover transitions. However, validation results highlighted limitations, some of which can be addressed by improving training data, creating a more robust image feature space, adding contemporaneous Landsat 5 data to the inputs, and modifying definition sets to account for differences in temporal and spatial observational scales. The presented land cover extent and change data are available via the official WELD website (ftp://weldftp.cr.usgs.gov/CONUS_5Y_LandCover/ftp://weldftp.cr.usgs.gov/CONUS_5Y_LandCover/).

Remote Sensing of Environment↗
USGS2012

Fundamental questions of earthquake statistics, source behavior, and the estimation of earthquake probabilities from possible foreshocks

Estimates of the probability that an M L 4.8 earthquake, which occurred near the southern end of the San Andreas fault on 24 March 2009, would be followed by an M 7 mainshock over the following three days vary from 0.0009 using a Gutenberg–Richter model of aftershock statistics (Reasenberg and Jones, 1989) to 0.04 using a statistical model of foreshock behavior and long‐term estimates of large earthquake probabilities, including characteristic earthquakes (Agnew and Jones, 1991). I demonstrate that the disparity between the existing approaches depends on whether or not they conform to Gutenberg–Richter behavior. While Gutenberg–Richter behavior is well established over large regions, it could be violated on individual faults if they have characteristic earthquakes or over small areas if the spatial distribution of large‐event nucleations is disproportional to the rate of smaller events. I develop a new form of the aftershock model that includes characteristic behavior and combines the features of both models. This new model and the older foreshock model yield the same results when given the same inputs, but the new model has the advantage of producing probabilities for events of all magnitudes, rather than just for events larger than the initial one. Compared with the aftershock model, the new model has the advantage of taking into account long‐term earthquake probability models. Using consistent parameters, the probability of an M 7 mainshock on the southernmost San Andreas fault is 0.0001 for three days from long‐term models and the clustering probabilities following the M L 4.8 event are 0.00035 for a Gutenberg–Richter distribution and 0.013 for a characteristic‐earthquake magnitude–frequency distribution. Our decisions about the existence of characteristic earthquakes and how large earthquakes nucleate have a first‐order effect on the probabilities obtained from short‐term clustering models for these large events.

Bulletin of the Seismological Society of America↗