SEARCH · USGS Science

Results for “Spatial Statistics”

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

SourceResearch connectionFrom yearThrough year

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

At least 469 records · Page 26Linked to original sources

USGS2013

The use of process models to inform and improve statistical models of nitrate occurrence, Great Miami River Basin, southwestern Ohio

Statistical models of nitrate occurrence in the glacial aquifer system of the northern United States, developed by the U.S. Geological Survey, use observed relations between nitrate concentrations and sets of explanatory variables—representing well-construction, environmental, and source characteristics— to predict the probability that nitrate, as nitrogen, will exceed a threshold concentration. However, the models do not explicitly account for the processes that control the transport of nitrogen from surface sources to a pumped well and use area-weighted mean spatial variables computed from within a circular buffer around the well as a simplified source-area conceptualization. The use of models that explicitly represent physical-transport processes can inform and, potentially, improve these statistical models. Specifically, groundwater-flow models simulate advective transport—predominant in many surficial aquifers— and can contribute to the refinement of the statistical models by (1) providing for improved, physically based representations of a source area to a well, and (2) allowing for more detailed estimates of environmental variables. A source area to a well, known as a contributing recharge area, represents the area at the water table that contributes recharge to a pumped well; a well pumped at a volumetric rate equal to the amount of recharge through a circular buffer will result in a contributing recharge area that is the same size as the buffer but has a shape that is a function of the hydrologic setting. These volume-equivalent contributing recharge areas will approximate circular buffers in areas of relatively flat hydraulic gradients, such as near groundwater divides, but in areas with steep hydraulic gradients will be elongated in the upgradient direction and agree less with the corresponding circular buffers. The degree to which process-model-estimated contributing recharge areas, which simulate advective transport and therefore account for local hydrologic settings, would inform and improve the development of statistical models can be implicitly estimated by evaluating the differences between explanatory variables estimated from the contributing recharge areas and the circular buffers used to develop existing statistical models. The larger the difference in estimated variables, the more likely that statistical models would be changed, and presumably improved, if explanatory variables estimated from contributing recharge areas were used in model development. Comparing model predictions from the two sets of estimated variables would further quantify—albeit implicitly—how an improved, physically based estimate of explanatory variables would be reflected in model predictions. Differences between the two sets of estimated explanatory variables and resultant model predictions vary spatially; greater differences are associated with areas of steep hydraulic gradients. A direct comparison, however, would require the development of a separate set of statistical models using explanatory variables from contributing recharge areas. Area-weighted means of three environmental variables—silt content, alfisol content, and depth to water from the U.S. Department of Agriculture State Soil Geographic (STATSGO) data—and one nitrogen-source variable (fertilizer-application rate from county data mapped to Enhanced National Land Cover Data 1992 (NLCDe 92) agricultural land use) can vary substantially between circular buffers and volume-equivalent contributing recharge areas and among contributing recharge areas for different sets of well variables. The differences in estimated explanatory variables are a function of the same factors affecting the contributing recharge areas as well as the spatial resolution and local distribution of the underlying spatial data. As a result, differences in estimated variables between circular buffers and contributing recharge areas are complex and site specific as evidenced by differences in estimated variables for circular buffers and contributing recharge areas of existing public-supply and network wells in the Great Miami River Basin. Large differences in areaweighted mean environmental variables are observed at the basin scale, determined by using the network of uniformly spaced hypothetical wells; the differences have a spatial pattern that generally is similar to spatial patterns in the underlying STATSGO data. Generally, the largest differences were observed for area-weighted nitrogen-application rate from county and national land-use data; the basin-scale differences ranged from -1,600 (indicating a larger value from within the volume-equivalent contributing recharge area) to 1,900 kilograms per year (kg/yr); the range in the underlying spatial data was from 0 to 2,200 kg/yr. Silt content, alfisol content, and nitrogen-application rate are defined by the underlying spatial data and are external to the groundwater system; however, depth to water is an environmental variable that can be estimated in more detail and, presumably, in a more physically based manner using a groundwater-flow model than using the spatial data. Model-calculated depths to water within circular buffers in the Great Miami River Basin differed substantially from values derived from the spatial data and had a much larger range. Differences in estimates of area-weighted spatial variables result in corresponding differences in predictions of nitrate occurrence in the aquifer. In addition to the factors affecting contributing recharge areas and estimated explanatory variables, differences in predictions also are a function of the specific set of explanatory variables used and the fitted slope coefficients in a given model. For models that predicted the probability of exceeding 1 and 4 milligrams per liter as nitrogen (mg/L as N), predicted probabilities using variables estimated from circular buffers and contributing recharge areas generally were correlated but differed significantly at the local and basin scale. The scale and distribution of prediction differences can be explained by the underlying differences in the estimated variables and the relative weight of the variables in the statistical models. Differences in predictions of exceeding 1 mg/L as N, which only includes environmental variables, generally correlated with the underlying differences in STATSGO data, whereas differences in exceeding 4 mg/L as N were more spatially extensive because that model included environmental and nitrogen-source variables. Using depths to water from within circular buffers derived from the spatial data and depths to water within the circular buffers calculated from the groundwater-flow model, restricted to the same range, resulted in large differences in predicted probabilities. The differences in estimated explanatory variables between contributing recharge areas and circular buffers indicate incorporation of physically based contributing recharge area likely would result in a different set of explanatory variables and an improved set of statistical models. The use of a groundwater-flow model to improve representations of source areas or to provide more-detailed estimates of specific explanatory variables includes a number of limitations and technical considerations. An assumption in these analyses is that (1) there is a state of mass balance between recharge and pumping, and (2) transport to a pumped well is under a steady state flow field. Comparison of volumeequivalent contributing recharge areas under steady-state and transient transport conditions at a location in the southeastern part of the basin shows the steady-state contributing recharge area is a reasonable approximation of the transient contributing recharge area after between 10 and 20 years of pumping. The first assumption is a more important consideration for this analysis. A gradient effect refers to a condition where simulated pumping from a well is less than recharge through the corresponding contributing recharge area. This generally takes place in areas with steep hydraulic gradients, such as near discharge locations, and can be mitigated using a finer model discretization. A boundary effect refers to a condition where recharge through the contributing recharge area is less than pumping. This indicates other sources of water to the simulated well and could reflect a real hydrologic process. In the Great Miami River Basin, large gradient and boundary effects—defined as the balance between pumping and recharge being less than half—occurred in 5 and 14 percent of the basin, respectively. The agreement between circular buffers and volume-equivalent contributing recharge areas, differences in estimated variables, and the effect on statisticalmodel predictions between the population of wells with a balance between pumping and recharge within 10 percent and the population of all wells were similar. This indicated process-model limitations did not affect the overall findings in the Great Miami River Basin; however, this would be model specific, and prudent use of a process model needs to entail a limitations analysis and, if necessary, alterations to the model.

Ohio↗

USGS2021

Development of regression equations for the estimation of the magnitude and frequency of floods at rural, unregulated gaged and ungaged streams in Puerto Rico through water year 2017

The methods of computation and estimates of the magnitude of flood flows were updated for the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent chance exceedance levels for 91 streamgages on the main island of Puerto Rico by using annual peak-flow data through 2017. Since the previous flood frequency study in 1994, the U.S. Geological Survey has collected additional peak flows at additional streamgages, and Puerto Rico has experienced numerous flood events. This updated study was performed using longer annual peak-flow datasets from more stations to provide more representative equations to predict flood flows. Screening criteria for these streamgages included 10 or more years of annual peak-flow data, unregulated flow, and less than 10 percent impervious drainage area. The magnitude and frequency of floods at selected streamgages in Puerto Rico were estimated using updated methods outlined in Bulletin 17C. The new procedures include a regional skew analysis that incorporates Bayesian regression techniques, the Expected Moments Algorithm to better represent missing record and estimate parameters of the log-Pearson Type III distribution, and the Multiple Grubbs-Beck test for low outlier detection. Regional regression equations were developed to estimate peak-flow statistics at ungaged locations by using selected basin and climatic characteristics as explanatory variables. These variables were determined from digital spatial datasets and geographic information systems by using the most recent data available. Ordinary least-squares regression techniques were used to filter the basin characteristics and determine two separate regions, region 1 (west) and region 2 (east), based on residuals. A generalized least-squares procedure was used to account for cross-correlation of sites and develop the final set of equations that have drainage area as the only explanatory variable. The average standard errors of prediction ranged from 18.7 to 46.7 percent in region 1 and 33.4 to 57.6 percent in region 2 for all annual exceedance probabilities (AEPs) examined. The updated statistics showed a greater accuracy of prediction when compared to those from the previous study using drainage area as the only explanatory variable for all AEPs examined in region 1 and the 0.01 and 0.002 AEP flows for region 2. When compared to equations developed in the previous study that have drainage area, mean annual rainfall, and (or) depth-to-rock as explanatory variables, the updated statistics show a greater accuracy of prediction in region 1 at AEP flows of 0.02 and lower (that is, higher flows). Those developed for region 2 do not show a greater accuracy of prediction for any AEP flows when compared to the equations having multiple explanatory variables in the previous study. The calculated regression equations, basin characteristics, and at-site statistics will be incorporated into the U.S. Geological Survey web application, StreamStats ( https://streamstats.usgs.gov/ss/ ). This application allows users to select a location on a stream, whether gaged or ungaged, to obtain estimates of basin characteristics and flow statistics.

Puerto Rico↗

USGS1986

Predictive spatial modeling of narcotic crop growth patterns

Spatial models for predicting the geographic distribution of marijuana crops have been developed and are being evaluated for use in law enforcement programs. The models are based on growing condition preferences and on psychological inferences regarding grower behavior. Experiences of local law officials were used to derive the initial model, which was updated and improved as data from crop finds were archived and statistically analyzed. The predictive models are changed as crop locations are moved in response to the pressures of law enforcement. The models use spatial data in a raster geographic information system. The spatial data are derived from the U.S. Geological Survey's US GeoData, standard 7.5-minute topographic quadrangle maps, interpretations of aerial photographs, and thematic maps. Updating of cultural patterns, canopy closure, and other dynamic features is conducted through interpretation of aerial photographs registered to the 7.5-minute quadrangle base. The model is used to numerically weight various data layers that have been processed using spread functions, edge definition, and categorization. The building of the spatial data base, model development, model application, product generation, and use are collectively referred to as the Area Reduction Program (ARP). The goal of ARP is to provide law enforcement officials with tactical maps that show the most likely locations for narcotic crops.

Open-File Report↗

USGS2012

Crop classification modelling using remote sensing and environmental data in the Greater Platte River Basin, USA

With an ever expanding population, potential climate variability and an increasing demand for agriculture-based alternative fuels, accurate agricultural land-cover classification for specific crops and their spatial distributions are becoming critical to researchers, policymakers, land managers and farmers. It is important to ensure the sustainability of these and other land uses and to quantify the net impacts that certain management practices have on the environment. Although other quality crop classification products are often available, temporal and spatial coverage gaps can create complications for certain regional or time-specific applications. Our goal was to develop a model capable of classifying major crops in the Greater Platte River Basin (GPRB) for the post-2000 era to supplement existing crop classification products. This study identifies annual spatial distributions and area totals of corn, soybeans, wheat and other crops across the GPRB from 2000 to 2009. We developed a regression tree classification model based on 2.5 million training data points derived from the National Agricultural Statistics Service (NASS) Cropland Data Layer (CDL) in relation to a variety of other relevant input environmental variables. The primary input variables included the weekly 250 m US Geological Survey Earth Observing System Moderate Resolution Imaging Spectroradiometer normalized differential vegetation index, average long-term growing season temperature, average long-term growing season precipitation and yearly start of growing season. An overall model accuracy rating of 78% was achieved for a test sample of roughly 215 000 independent points that were withheld from model training. Ten 250 m resolution annual crop classification maps were produced and evaluated for the GPRB region, one for each year from 2000 to 2009. In addition to the model accuracy assessment, our validation focused on spatial distribution and county-level crop area totals in comparison with the NASS CDL and county statistics from the US Department of Agriculture (USDA) Census of Agriculture. The results showed that our model produced crop classification maps that closely resembled the spatial distribution trends observed in the NASS CDL and exhibited a close linear agreement with county-by-county crop area totals from USDA census data ( R 2 = 0.90).

Greater Platte River Basin↗

USGS2001

Spatial paradigms of lotic diatom distribution: A landscape ecology perspective

Spatial distributional patterns of benthic diatoms and their relation to current velocity were investigated in an unshaded cobble-bottom reach of White Creek (Washington County, NY). On 27 August 1999, diatoms were sampled and current velocity and depth were measured on a regular square sampling grid with a grain size of 0.01 m2, interval of 0.5 m, and extent of 16 m2. The relative abundance of the 18 common diatom species enumerated in the 81 samples was subjected to detrended correspondence analysis (DCA). The first axis (DCA1) explained 51% of the variance in diatom data and separated the samples according to current regimes. The spatial autocorrelation of DCA1 sample scores in deposition and erosion regions of White Creek was determined by Moran's I statistic to indicate patch size. In White Creek the patch length of all diatom communities was more than 3.1 m, whereas the patch width was 1 m in the deposition region and 0.5 m in the erosion region. There were 5 dominant diatom taxa, Achnanthes minutissima Ku??tz. et vars, Fragilaria capucina Dezmazie??res et vars, F. crotonensis Kitt., Diatoma vulgaris Bory, and Synedra ulna (Nitz.) Ehr. et vars. The patch length of the dominant species varied from 1 to more than 4.1 m, whereas the patch width, if defined, was 0.5 m. Achnanthes minutissima and F. capucina, the two diatom species with the highest relative abundance, displayed spatially structured patches of low abundance and comparatively random patches of high abundance, suggesting broad scale abiotic control of species performance in low abundance regions and finer scale biotic control of high abundance areas. Another objective of this study was to test the hypothesis that higher current velocities, which generally impede immigration, would increase randomness and complexity (i.e. homogeneity of diatom distributional patterns). The spatial complexity in low versus high velocity transects was determined by calculating the respective fractal dimension (D) of DCA1 scores. D of DCA1 was higher in the higher current velocity transects, suggesting that spatial complexity and homogeneity of diatom communities increased in faster currents. Partial canonical correspondence analysis was conducted on diatom, environmental, and spatial data to assess how much of the variance in species distribution could be attributed to environmental (current velocity and depth) versus spatial factors. The variance of species data, explained by the environment (exclusively current velocity), was 38%; whereas space alone contributed only 10%, indicating that 1) current velocity was the major factor that controlled diatom distribution in streams and 2) there were other spatially dependent variables, most likely biotic, but their role in shaping diatom communities was minor.

Journal of Phycology↗

USGS2012

A high-resolution bioclimate map of the world: a unifying framework for global biodiversity research and monitoring

Aim: To develop a novel global spatial framework for the integration and analysis of ecological and environmental data. Location: The global land surface excluding Antarctica. Methods: A broad set of climate-related variables were considered for inclusion in a quantitative model, which partitions geographic space into bioclimate regions. Statistical screening produced a subset of relevant bioclimate variables, which were further compacted into fewer independent dimensions using principal components analysis (PCA). An ISODATA clustering routine was then used to classify the principal components into relatively homogeneous environmental strata. The strata were aggregated into global environmental zones based on the attribute distances between strata to provide structure and support a consistent nomenclature. Results: The global environmental stratification (GEnS) consists of 125 strata, which have been aggregated into 18 global environmental zones. The stratification has a 30 arcsec resolution (equivalent to 0.86 km2 at the equator). Aggregations of the strata were compared with nine existing global, continental and national bioclimate and ecosystem classifications using the Kappa statistic. Values range between 0.54 and 0.72, indicating good agreement in bioclimate and ecosystem patterns between existing maps and the GEnS. Main conclusions: The GEnS provides a robust spatial analytical framework for the aggregation of local observations, identification of gaps in current monitoring efforts and systematic design of complementary and new monitoring and research. The dataset is available for non-commercial use through the GEO portal (http://www.geoportal.org).

Global Ecology and Biogeography↗

USGS2015

Assessing juvenile salmon rearing habitat and associated predation risk in a lower Snake River reservoir

Subyearling fall Chinook salmon (Oncorhynchus tshawytscha) in the Columbia River basin exhibit a transient rearing strategy and depend on connected shoreline habitats during freshwater rearing. Impoundment has greatly reduced the amount of shallow-water rearing habitat that is exacerbated by the steep topography of reservoirs. Periodic dredging creates opportunities to strategically place spoils to increase the amount of shallow-water habitat for subyearlings while at the same time reducing the amount of unsuitable area that is often preferred by predators. We assessed the amount and spatial arrangement of subyearling rearing habitat in Lower Granite Reservoir on the Snake River to guide future habitat improvement efforts. A spatially explicit habitat assessment was conducted using physical habitat data, two-dimensional hydrodynamic modelling and a statistical habitat model in a geographic information system framework. We used field collections of subyearlings and a common predator [smallmouth bass (Micropterus dolomieu)] to draw inferences about predation risk within specific habitat types. Most of the high-probability rearing habitat was located in the upper half of the reservoir where gently sloping landforms created low lateral bed slopes and shallow-water habitats. Only 29% of shorelines were predicted to be suitable (probability >0.5) for subyearlings, and the occurrence of these shorelines decreased in a downstream direction. The remaining, less suitable areas were composed of low-probability habitats in unmodified (25%) and riprapped shorelines (46%). As expected, most subyearlings were found in high-probability habitat, while most smallmouth bass were found in low-probability locations. However, some subyearlings were found in low-probability habitats, such as riprap, where predation risk could be high. Given their transient rearing strategy and dependence on shoreline habitats, subyearlings could benefit from habitat creation efforts in the lower reservoir where high-probability habitat is generally lacking. Published 2015. This article is a U.S. Government work and is in the public domain in the USA.

Lower Snake river reservior↗

USGS2020

Bioclimatic modeling of potential vegetation types as an alternative to species distribution models for projecting plant species shifts under changing climates

Land managers need new tools for planning novel futures due to climate change. Species distribution modeling (SDM) has been used extensively to predict future distributions of species under different climates, but their map products are often too coarse for fine-scale operational use. In this study we developed a flexible, efficient, and robust method for mapping current and future distributions and abundances of vegetation species and communities at the fine spatial resolutions that are germane to land management. First, we mapped Potential Vegetation Types (PVTs) using conventional statistical modeling techniques (Random Forests) that used bioclimatic ecosystem process and climate variables as predictors. We obtained over 50% accuracy across 13 mapped PVTs for our study area. We then applied future climate projections as climate input to the Random Forest model to generate future PVT maps, and used field data describing the occurrence of tree and non-tree species in each PVT category to model and map species distribution for current and future climate. These maps were then compared to two previous SDM mapping efforts with over 80% agreement and equivalent accuracy. Because PVTs represent the biophysical potential of the landscape to support vegetation communities as opposed to the vegetation that currently exists, they can be readily linked to climate forecasts and correlated with other, climate-sensitive ecological processes significant in land management, such as fire regimes and site productivity.

Montana↗

USGS2018

Classifying physiographic regimes on terrain and hydrologic factors for adaptive generalization of stream networks

Automated generalization software must accommodate multi-scale representations of hydrographic networks across a variety of geographic landscapes, because scale-related hydrography differences are known to vary in different physical conditions. While generalization algorithms have been tailored to specific regions and landscape conditions by several researchers in recent years, the selection and characterization of regional conditions have not been formally defined nor statistically validated. This paper undertakes a systematic classification of landscape types in the conterminous United States to spatially subset the country into workable units, in preparation for systematic tailoring of generalization workflows that preserve hydrographic characteristics. The classification is based upon elevation, standard deviation of elevation, slope, runoff, drainage and bedrock density, soil and bedrock permeability, area of inland surface water, infiltration-excess of overland flow, and a base flow index. A seven class solution shows low misclassification rates except in areas of high landscape diversity such as the Appalachians, Rocky Mountains, and Western coastal regions.

International Journal of Cartography↗

USGS2004

Variation in the population structure of Yukon River chum and coho salmon: Evaluating the potential impact of localized habitat degradation

We used microsatellite and mitochondrial DNA-restriction fragment length polymorphism (mtDNA-RFLP) analyses to test the hypothesis that chum salmon Oncorhynchus keta and coho salmon O. kisutch in the Yukon River, Alaska, exhibit population structure at differing spatial scales. If the hypothesis is true, then the risk of losing genetic diversity because of habitat degradation from a gold mine near a Yukon River tributary could differ between the two species. For each species, collections were made from two tributaries in both the Innoko and Tanana rivers, which are tributaries to the lower and middle Yukon River. The results revealed a large difference in the degree and spatial distribution of population structure between the two species. For chum salmon, the microsatellite loci (F-statistic [FST] = 0.021) and mtDNA (F ST = -0.008) revealed a low degree of interpopulation genetic diversity on a relatively large geographic scale. This large-scale population structure should minimize, although not eliminate, the risk of genetic diversity loss due to localized habitat degradation. For coho salmon, the microsatellites (FST = 0.091) and mtDNA (FST = 0.586) revealed a high degree of interpopulation genetic diversity on a relatively small geographic scale. This small-scale population structure suggests that coho salmon are at a relatively high risk of losing genetic diversity due to lo-calized habitat degradation. Our study underscores the importance of a multispecies approach for evaluating the potential impact of land-use activities on the genetic diversity of Pacific salmon.

Transactions of the American Fisheries Society↗

USGS2010

A climate trend analysis of Kenya-August 2010

Introduction This brief report draws from a multi-year effort by the United States Agency for International Development's Famine Early Warning System Network (FEWS NET) to monitor and map rainfall and temperature trends over the last 50 years (1960-2009) in Kenya. Observations from seventy rainfall gauges and seventeen air temperature stations were analyzed for the long rains period, corresponding to March through June (MAMJ). The data were quality controlled, converted into 1960-2009 trend estimates, and interpolated using a rigorous geo-statistical technique (kriging). Kriging produces standard error estimates, and these can be used to assess the relative spatial accuracy of the identified trends. Dividing the trends by the associated errors allows us to identify the relative certainty of our estimates (Funk and others, 2005; Verdin and others, 2005; Brown and Funk, 2008; Funk and Verdin, 2009). Assuming that the same observed trends persist, regardless of whether or not these changes are due to anthropogenic or natural cyclical causes, these results can be extended to 2025, providing critical, and heretofore missing information about the types and locations of adaptation efforts that may be required to improve food security.

Fact Sheet↗

USGS1981

Pattern-recognition program modified and applied to southeastern United States seismicity

The pattern-recognition technique of Bongard and others (1966), Gelfand and others (1976), and Briggs and others (1977) was modified and applied in Southeastern United States to try to detect combinations or interactions of geologic, geophysical, and topographic parameters characteristic of areas with moderate and large earthquakes. The modification involved replacing the operator-defined criteria for significance of association of parameters and earthquakes with statistical estimates of significance. This was done in order to decrease the probability of random associations and to increase the flexibility of the pattern-recognition technique. The binomial distribution was used to determine the probability that each particular association could have occurred by chance. Bonferroni's inequality was invoked to determine a conservative significance level for each of the binomial tests. However, without being able to quantify the dependence among the associations, it is not possible to statistically evaluate the group of binomial tests as a whole. Preliminary results identify four independent associations characterizing spatially independent earthquake epicenters, or "dangerous" points, and seven independent associations characterizing the intervening "safe" points. When all points were reclassified using these traits, all epicenters were correctly classified as dangerous. In addition, three "safe" points were reclassified as dangerous and appear to be potential sites for future earthquakes. However, it is not possible to determine the degree of dangerousness of the points using this reclassification scheme. Two dangerous zones emerge in the northeast and southwest portions of the study area. Epicenters outside of these zones, including those of the large historic earthquakes (MMI > VIII) in Giles County, Virginia, and Charleston, South Carolina, have not been well characterized by the significant associations and may have fundamentally different causes.

southeastern United States↗

USGS2014

Assessing potential effects of highway runoff on receiving-water quality at selected sites in Oregon with the Stochastic Empirical Loading and Dilution Model (SELDM)

In 2012, the U.S. Geological Survey and the Oregon Department of Transportation began a cooperative study to demonstrate use of the Stochastic Empirical Loading and Dilution Model (SELDM) for runoff-quality analyses in Oregon. SELDM can be used to estimate stormflows, constituent concentrations, and loads from the area upstream of a stormflow discharge site, from the site of interest and in the receiving waters downstream of the discharge. SELDM also can be used to assess the potential effectiveness of best management practices (BMP) for mitigating potential effects of runoff in receiving waters. Nominally, SELDM is a highway-runoff model, but it is well suited for analysis of runoff from other land uses as well. This report provides case studies and examples to demonstrate stochastic-runoff modeling concepts and to demonstrate application of the model. Basin characteristics from six Oregon highway study sites were used to demonstrate various applications of the model. The highway catchment and upstream basin drainage areas of these study sites ranged from 3.85 to 11.83 acres and from 0.16 to 6.56 square miles, respectively. The upstream basins of two sites are urbanized, and the remaining four sites are less than 5 percent impervious. SELDM facilitates analysis by providing precipitation, pre-storm streamflow, and other variables by region or from hydrologically similar sites. In Oregon, there can be large variations in precipitation and streamflow among nearby sites. Therefore, spatially interpolated geographic information system data layers containing storm-event precipitation and pre-storm streamflow statistics specific to Oregon were created for the study using Kriging techniques. Concentrations and loads of cadmium, chloride, chromium, copper, iron, lead, nickel, phosphorus, and zinc were simulated at the six Oregon highway study sites by using statistics from sites in other areas of the country. Water‑quality datasets measured at hydrologically similar basins in the vicinity of the study sites in Oregon were selected and compiled to estimate stormflow-quality statistics for the upstream basins. The quality of highway runoff and some upstream stormflow constituents were simulated by using statistical moments (average, standard deviation, and skew) of the logarithms of data. Some upstream stormflow constituents were simulated by using transport curves, which are relations between stormflow and constituent concentrations. Stochastic analyses were done by using SELDM to demonstrate use of the model and to illustrate the types of information that stochastic analyses may provide: 1. An analysis was done to demonstrate use of dilution factors as an initial reconnaissance tool for comparing relative risk among sites. 2. An analysis of hardness-dependent, water-quality criteria was done to illustrate the effects of variations in hardness and flow on the application and interpretation of such criteria. This analysis shows that hardness-dependent criteria can vary by an order of magnitude among storm events because hardness is diluted by stormflows. 3. An analysis of uncertainties in input and output values was done to demonstrate that properly selected robust datasets are needed to represent conditions at a site of interest. This analysis shows that the rate of water-quality exceedances that are measured or simulated may depend on sample size and the luck of the draw. 4. An analysis was done to demonstrate that SELDM and other Monte Carlo models may generate extreme values from input statistics, which may or may not be feasible based on physicochemical or hydrological limits. 5. An analysis of BMP modeling methods was done to demonstrate use of the model for estimating treatment requirements for meeting water-quality objectives. 6. An analysis of the use of grab sampling and nonstochastic upstream modeling methods was done to evaluate the potential effects on modeling outcomes. Additional analyses using surrogate water-quality datasets for the upstream basin and highway catchment were provided for six Oregon study sites to illustrate the risk-based information that SELDM will produce. These analyses show that the potential effects of highway runoff on receiving-water quality downstream of the outfall depends on the ratio of drainage areas (dilution), the quality of the receiving water upstream of the highway, and the concentration of the criteria of the constituent of interest. These analyses also show that the probability of exceeding a water-quality criterion may depend on the input statistics used, thus careful selection of representative values is important.

Oregon↗

USGS1996

Water-quality assessment of the Connecticut, Housatonic, and Thames River Basins study unit: Analysis of available data on nutrients, suspended sediments, and pesticides, 1972-92

This retrospective report examines available nutrient, suspended sediment, and pesticide data in surface and ground water in the Connecticut, Housatonic and Thames Rivers Study Unit of the National Water-Quality Assessment Program. The purpose of this study is to improve the understanding of natural and anthropogenic factors affecting water quality in the study unit. Waterquality data were acquired from various sources, primarily, the U.S. Geological Survey and the U.S. Environmental Protection Agency. The report examines data for water years 1972-92, focusing on 1980-92, although it also includes additional data from as early as 1905. The study unit lies within the New England Physiographic Province and altitudes range from sea level in coastal Connecticut to 6,288 feet above sea level at Mount Washington, New Hampshire. Two major aquifer types underlie the study unit unconsolidated glacial deposits and fractured bedrock. The climate generally is temperate and humid, with four distinct seasons. Average annual precipitation ranges from 34 to 65 inches. The study unit has a population of about 4.5 million, which is most highly concentrated in southwestern Connecticut and along the south-central region of the Connecticut River Valley. Surface-water-quality data were screened to provide information about sites with adequate numbers of analyses (50) over sufficiently long periods (1980-90) to enable valid statistical analyses. In order to compare effects of different types of land use on surface-water quality, examination of data required application of several statistical and graphical techniques, including mapping, histograms, boxplots, concentration-discharge plots, trend analysis, and load estimation. Spatial and temporal analysis of surface-water-quality data indicated that, with a single exception, only stations in the Connecticut water-quality network had sufficient data collected over adequately long time periods to use in detailed analyses. Ground-water nutrient and pesticide data were compiled from several Federal and State agencies, primarily the U.S. Geological Survey, U.S. Environmental Protection Agency, and Connecticut Department of Health Services. Nutrient data were available for several thousand wells; nitrite plus nitrate as nitrogen was the most commonly reported constituent. Most wells with nutrient data are in Massachusetts and Connecticut. Relative to nutrient data in ground and surface water, pesticide data are less common. Pesticide data were available for slightly more than 200 surface-water sites and less than 500 wells; about 95 percent of the wells are completed in stratified-drift or till aquifers. Data for 81 pesticide compounds were available in various data bases. 2,4-D and silvex were the most commonly detected herbicides in surface water and dieldrin and diazinon were the most commonly detected insecticides. Most surface-water pesticide samples and detections are from bed sediment, but much of the data are not recent. Ethylene dibromide (EDB), a soil fumigant used in tobacco farming was detected in 268 wells in a 50 square-mile area of north-central Connecticut; EDB contamination also was detected in wells in Massachusetts. Atrazine, an herbicide commonly used in corn farming, commonly was detected in wells installed in tilled agricultural fields. Corn herbicides were commonly detected in the northern part of the study unit, although the sampling has been less frequent than in the southern part of the study unit. Pesticides were seldom detected in public-supply wells in Connecticut. Urban sites with relatively high population densities and high concentrations of dischargers were characterized by having the highest nutrient concentrations and loads when adjusted for differences in drainage area or population. Particularly, the Pequabuck, Naugatuck, and Quinnipiac River Basins were characterized by high nutrient concentrations median total nitrogen concentrations ranged from 3.3 to 4.2 mg/L (milligrams per liter) and median total phosphorus concentrations ranged from 0.42 to 0.8 mg/L. In contrast, the predominantly forested and low density residential land-use sites, such as Saugatuck and Salmon River Basins, were characterized by low nutrient concentrations median total nitrogen ranged from 0.50 to 0.60 mg/L and median total phosphorus concentrations ranged from 0.01 to 0.02 mg/L. Estimated total nitrogen loadings in median discharge years ranged from 940 kilograms per square mile at the Salmon River near East Hampton, Conn., to 5,800 kilograms per square mile at the Naugatuck River at Beacon Falls, Conn. Water quality, in terms of nutrient concentrations and areally adjusted loadings, for sites with large drainage basins integrating a wide variety of land-use categories fell between the extremes of the urban and forested sites total nitrogen was 1,400 kilograms per square mile per year at the Connecticut River at Thompsonville, Conn. Nitrate concentrations in ground water occasionally exceeded the safe drinking-water standard of 10 mg/L as nitrogen. The greatest number of detections exceeding the standard, however, were not in public-water supplies but in shallow observation wells in agricultural settings (the most frequently sampled type of well). None of the public-supply wells in Massachusetts exceeded the standard. Although nitrate concentrations for Vermont and New Hampshire generally were low, few data were available and those were seldom reported on the basis of drainage basin, making analysis difficult. Trend analysis indicated that flow-adjusted concentrations of total and dissolved phosphorus generally decreased during the period of analysis, however, total nitrogen did not change substantially. Decreases in ammonia concentrations with time were usually accompanied by increases in nitrate, suggesting improvements in sewage treatment. The lack of adequate data from more or less exclusively agricultural areas points to the need for further study of the effects of fanning on surface-water quality in the study unit. Furthermore, additional information is needed on the rates, transformations, and movements of nutrients and other materials through and between the aquatic and terrestrial components of the study unit.

Connecticut, Massachusetts, New Hampshire, Quebec,↗

USGS2017

A three-dimensional mapping of the ocean based on environmental data

The existence, sources, distribution, circulation, and physicochemical nature of macroscale oceanic water bodies have long been a focus of oceanographic inquiry. Building on that work, this paper describes an objectively derived and globally comprehensive set of 37 distinct volumetric region units, called ecological marine units (EMUs). They are constructed on a regularly spaced ocean point-mesh grid, from sea surface to seafloor, and attributed with data from the 2013 World Ocean Atlas version 2. The point attribute data are the means of the decadal averages from a 57-year climatology of six physical and chemical environment parameters (temperature, salinity, dissolved oxygen, nitrate, phosphate, and silicate). The database includes over 52 million points that depict the global ocean in x, y, and z dimensions. The point data were statistically clustered to define the 37 EMUs, which represent physically and chemically distinct water volumes based on spatial variation in the six marine environmental characteristics used. The aspatial clustering to produce the 37 EMUs did not include point location or depth as a determinant, yet strong geographic and vertical separation was observed. Twenty-two of the 37 EMUs are globally or regionally extensive, and account for 99% of the ocean volume, while the remaining 15 are smaller and shallower, and occur around coastal features. We assessed the vertical distribution of EMUs in the water column and placed them into classical depth zones representing epipelagic (0 m to 200 m), mesopelagic (200 m to 1,000 m), bathypelagic (1,000 m to 4,000 m) and abyssopelagic (>4,000 m) layers. The mapping and characterization of the EMUs represent a new spatial framework for organizing and understanding the physical, chemical, and ultimately biological properties and processes of oceanic water bodies. The EMUs are an initial objective partitioning of the ocean using long-term historical average data, and could be extended in the future by adding new classification variables and by introducing functionality to develop time-specific EMU distribution maps. The EMUs are an open-access resource, and as both a standardized geographic framework and a baseline physicochemical characterization of the oceanic environment, they are intended to be useful for disturbance assessments, ecosystem accounting exercises, conservation priority setting, and marine protected area network design, along with other research and management applications.

Oceanography↗

USGS2008

Landscape complexity and soil moisture variation in south Georgia, USA, for remote sensing applications

This research addressed the temporal and spatial variation of soil moisture (SM) in a heterogeneous landscape. The research objective was to investigate soil moisture variation in eight homogeneous 30 by 30 m plots, similar to the pixel size of a Landsat Thematic Mapper (TM) or Enhanced Thematic Mapper plus (ETM+) image. The plots were adjacent to eight stations of an in situ soil moisture network operated by the United States Department of Agriculture-Agriculture Research Service USDA-ARS in Tifton, GA. We also studied five adjacent agricultural fields to examine the effect of different landuses/land covers (LULC) (grass, orchard, peanuts, cotton and bare soil) on the temporal and spatial variation of soil moisture. Soil moisture field data were collected on eight occasions throughout 2005 and January 2006 to establish comparisons within and among eight homogeneous plots. Consistently throughout time, analysis of variance (ANOVA) showed high variation in the soil moisture behavior among the plots and high homogeneity in the soil moisture behavior within them. A precipitation analysis for the eight sampling dates throughout the year 2005 showed similar rainfall conditions for the eight study plots. Therefore, soil moisture variation among locations was explained by in situ local conditions. Temporal stability geostatistical analysis showed that soil moisture has high temporal stability within the small plots and that a single point reading can be used to monitor soil moisture status for the plot within a maximum 3% volume/volume (v/v) soil moisture variation. Similarly, t-statistic analysis showed that soil moisture status in the upper soil layer changes within 24 h. We found statistical differences in the soil moisture between the different LULC in the agricultural fields as well as statistical differences between these fields and the adjacent 30 by 30 m plots. From this analysis, it was demonstrated that spatial proximity is not enough to produce similar soil moisture, since t-test's among adjacent plots with different LULCs showed significant differences. These results confirm that a remote sensing approach that considers homogeneous LULC landscape fragments can be used to identify landscape units of similar soil moisture behavior under heterogeneous landscapes. In addition, the in situ USDA-ARS network will serve better in remote sensing studies in which sensors with fine spatial resolution are evaluated. This study is a first step towards identifying landscape units that can be monitored using the single point reading of the USDA-ARS stations network. ?? 2008 Elsevier B.V.

Journal of Hydrology↗

USGS2015

An empirical evaluation of landscape energetic models: Mallard and American black duck space use during the non-breeding period

Bird conservation Joint Ventures are collaborative partnerships between public agencies and private organizations that facilitate habitat management to support waterfowl and other bird populations. A subset of Joint Ventures has developed energetic carrying capacity models (ECCs) to translate regional waterfowl population goals into habitat objectives during the non-breeding period. Energetic carrying capacity models consider food biomass, metabolism, and available habitat to estimate waterfowl carrying capacity within an area. To evaluate Joint Venture ECCs in the context of waterfowl space use, we monitored 33 female mallards ( Anas platyrhynchos ) and 55 female American black ducks ( A. rubripes ) using global positioning system satellite telemetry in the central and eastern United States. To quantify space use, we measured first-passage time (FPT: time required for an individual to transit across a circle of a given radius) at biologically relevant spatial scales for mallards (3.46 km) and American black ducks (2.30 km) during the non-breeding period, which included autumn migration, winter, and spring migration. We developed a series of models to predict FPT using Joint Venture ECCs and compared them to a biological null model that quantified habitat composition and a statistical null model, which included intercept and random terms. Energetic carrying capacity models predicted mallard space use more efficiently during autumn and spring migrations, but the statistical null was the top model for winter. For American black ducks, ECCs did not improve predictions of space use; the biological null was top ranked for winter and the statistical null was top ranked for spring migration. Thus, ECCs provided limited insight into predicting waterfowl space use during the non-breeding season. Refined estimates of spatial and temporal variation in food abundance, habitat conditions, and anthropogenic disturbance will likely improve ECCs and benefit conservation planners in linking non-breeding waterfowl habitat objectives with distribution and population parameters. Published 2015. This article is a U.S. Government work and is in the public domain in the USA.

Arkansas, Delaware, Louisiana, Michigan, New Jerse↗

USGS2012

Downscaling future climate scenarios to fine scales for hydrologic and ecological modeling and analysis

Introduction Evaluating the environmental impacts of climate change on water resources and biological components of the landscape is an integral part of hydrologic and ecological investigations, and the resultant land and resource management in the twenty-first century. Impacts of both climate and simulated hydrologic parameters on ecological processes are relevant at scales that reflect the heterogeneity and complexity of landscapes. At present, simulations of climate change available from global climate models [GCMs] require downscaling for hydrologic or ecological applications. Methods Using statistically downscaled future climate projections developed using constructed analogues, a methodology was developed to further downscale the projections spatially using a gradient-inverse-distance-squared approach for application to hydrologic modeling at 270-m spatial resolution. Results This paper illustrates a methodology to downscale and bias-correct national GCMs to subkilometer scales that are applicable to fine-scale environmental processes. Four scenarios were chosen to bracket the range of future emissions put forth by the Intergovernmental Panel on Climate Change. Fine-scale applications of downscaled datasets of ecological and hydrologic correlations to variation in climate are illustrated. Conclusions The methodology, which includes a sequence of rigorous analyses and calculations, is intended to reduce the addition of uncertainty to the climate data as a result of the downscaling while providing the fine-scale climate information necessary for ecological analyses. It results in new but consistent data sets for the US at 4 km, the southwest US at 270 m, and California at 90 m and illustrates the utility of fine-scale downscaling to analyses of ecological processes influenced by topographic complexity.

Ecological Processes↗