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 829 records · Page 46Linked to original sources

Spatial and temporal snowpack variation in the crown of the continent ecosystem

Snowpack related ecosystem changes such as glacier recession and alpine treeline advance have been documented in the Crown of the Continent Ecosystem (CCE) over the course of the previous 150 years. Using data from the Natural Resource Conservation Service's SNOTEL sites and snow course surveys, we examined the spatial and temporal variation in snowpack in the region. SNOTEL data suggest CCE snowpacks are larger and more persistent than in most regions of the Western U.S., and that water year precipitation, rather than mean temperature, is the primary control on April 1 snow water equivalent (SWE). Snow course data indicate a statistically significant downward trend in mean April 1 SWE for the period 1950-2001 but no statistically significant trend in mean May 1 SWE for the longer period 1922-2001. Further analysis reveals that variations in both April 1 and May 1 mean SWE are closely tied to the Pacific Decadal Oscillation, an ENSO-like interdecadal pattern of Pacific Ocean climate variability. Despite no significant trend in mean May 1 SWE between 1922-2001, glaciers in Glacier National Park receded steadily during this period, implying changing climatic conditions crossed a threshold for glacier mass balance maintenace sometime between the Little Ice Age glacial maxima and 1922.

Conference Paper↗

A global empirical model for near real-time assessment of seismically induced landslides

Earthquake-triggered landslides are a significant hazard in seismically active regions, but our ability to assess the hazard they pose in near real-time is limited. In this study, we present a new globally applicable model for seismically induced landslides based on the most comprehensive global dataset available; we use 23 landslide inventories that span a range of earthquake magnitudes and climatic and tectonic settings. We use logistic regression to relate the presence and distribution of earthquake-triggered landslides with spatially distributed estimates of ground shaking, topographic slope, lithology, land-cover type, and a topographic index designed to estimate variability in soil wetness to provide an empirical model of landslide distribution. We tested over 100 combinations of independent predictor variables to find the best-fitting model, using a diverse set of statistical tests. Blind validation tests show the model accurately estimates the distribution of available landslide inventories. The results indicate that the model is reliable and stable, with high “balanced accuracy” (correctly vs. incorrectly classified pixels) for the majority of test events. A cross validation analysis shows high balanced accuracy for a majority of events as well. By combining near-real time estimates of ground shaking with globally available landslide susceptibility data, this model provides a tool to estimate the distribution of co-seismic landslide hazard within minutes of the occurrence of any earthquake worldwide for which a USGS ShakeMap is available.

Journal of Geophysical Research: Earth Surface↗

Simultaneous Middle Pleistocene eruption of three widespread tholeiitic basalts in northern California (USA): Insights into crustal magma transport in an actively extending back arc

Mapping and chronology are central to understanding spatiotemporal volcanic trends in diverse tectonic settings. The Cascades back arc in northern California (USA) hosts abundant lava flows and normal faults, but tholeiitic basalts older than 200 ka are difficult to discriminate by classic mapping methods. Paleomagnetism and chemistry offer independent means of correlating basalts, including the Tennant, Dry Lake, and Hammond Crossing basalt fields. Paleomagnetic analysis of these chemically similar basalts yield notable overlap, with statistical analysis yielding 7 chances in 1,000,000 that their similar mean remanent directions are random. These basalts also have overlapping 40 Ar/ 39 Ar ages of 272.5 ± 30.6 ka (Tennant), 305.8 ± 23.9 ka (Dry Lake), and 300.4 ± 15.2 and 322.6 ± 17.4 ka (Hammond Crossing). Chemical and paleomagnetic analyses indicate that these spatially distributed basalts represent simultaneous (<100 yr uncertainty) eruptions, and thus we use 305.5 ± 9.8 ka (weighted mean) as the eruption age. Their vents align on a N25°W trend over a distance of 39 km. Tennant erupted the largest volume (3.55 ± 0.75 km 3 ) at the highest elevation; both factors decay to the south-southeast at Dry Lake (0.75 ± 0.15 km 3 ) and Hammond Crossing (0.15 ± 0.05 km 3 ). We propose vertical magma ascent beneath the Tennant vent area, where the most evolved, high-SiO 2 magma erupted, with lateral dike propagation in the brittle crust. Propagation was near orthogonal to east-west extension (0.3–0.6 mm/yr) along north-northwest–trending normal faults.

California↗

Analysis and machine mapping of the distribution of band recoveries

A method of calculating distance and bearing from banding site to recovery location based on the solution of a spherical triangle is presented. X and Y distances on an ordinate grid were applied to computer plotting of recoveries on a map. The advantages and disadvantages of tables of recoveries by State or degree block, axial lines, and distance of recovery from banding site for presentation and comparison of the spatial distribution of band recoveries are discussed. A special web-shaped partition formed by concentric circles about the point of banding and great circles at 30-degree intervals through the point of banding has certain advantages over other methods. Comparison of distributions by means of a X? contingency test is illustrated. The statistic V = X?/N can be used as a measure of difference between two distributions of band recoveries and its possible use is illustrated as a measure of the degree of migrational homing.

Special Scientific Report - Wildlife↗

Artificial intelligence based decision support for trumpeter swan management

The number of trumpeter swans (Cygnus buccinator) breeding in the Tri-State area where Montana, Idaho, and Wyoming come together has declined to just a few hundred pairs. However, these birds are part of the Rocky Mountain Population which additionally has over 3,500 birds breeding in Alberta, British Columbia, Northwest Territories, and Yukon Territory. To a large degree, these birds seem to have abandoned traditional migratory pathways in the flyway. Waterfowl managers have been interested in decision support tools that would help them explore simulated management scenarios in their quest towards reaching population recovery and the reestablishment of traditional migratory pathways. I have developed a decision support system to assist biologists with such management, especially related to wetland ecology. Decision support systems use a combination of models, analytical techniques, and information retrieval to help develop and evaluate appropriate alternatives. Swan management is a domain that is ecologically complex, and this complexity is compounded by spatial and temporal issues. As such, swan management is an inherently distributed problem. Therefore, the ecological context for modeling swan movements in response to management actions was built as a multiagent system of interacting intelligent agents that implements a queuing model representing swan migration. These agents accessed ecological knowledge about swans, their habitats, and flyway management principles from three independent expert systems. The agents were autonomous, had some sensory capability, and could respond to changing conditions. A key problem when developing ecological decision support systems is empirically determining that the recommendations provided are valid. Because Rocky Mountain trumpeter swans have been surveyed for a long period of time, I was able to compare simulated distributions provided by the system with actual field observations across 20 areas for the period 1988-2000. Applying the Matched Pairs Multivariate Permutation Test as a statistical tool was a new approach for comparing flyway distributions of waterfowl over time that seemed to work well. Based on this approach, the empirical evidence that I gathered led me to conclude that the base queuing model does accurately simulate swan distributions in the flyway. The system was insensitive to almost all model parameters tested. That remains perplexing, but might result from the base queuing model, itself, being particularly effective at representing the actual ecological diversity in the world of Rocky Mountain trumpeter swans, both spatial and temporally.

Idaho, Montana, Utah, Wyoming↗

Fire frequency in the Interior Columbia River Basin: Building regional models from fire history data

Fire frequency affects vegetation composition and successional pathways; thus it is essential to understand fire regimes in order to manage natural resources at broad spatial scales. Fire history data are lacking for many regions for which fire management decisions are being made, so models are needed to estimate past fire frequency where local data are not yet available. We developed multiple regression models and tree-based (classification and regression tree, or CART) models to predict fire return intervals across the interior Columbia River basin at 1-km resolution, using georeferenced fire history, potential vegetation, cover type, and precipitation databases. The models combined semiqualitative methods and rigorous statistics. The fire history data are of uneven quality; some estimates are based on only one tree, and many are not cross-dated. Therefore, we weighted the models based on data quality and performed a sensitivity analysis of the effects on the models of estimation errors that are due to lack of cross-dating. The regression models predict fire return intervals from 1 to 375 yr for forested areas, whereas the tree-based models predict a range of 8 to 150 yr. Both types of models predict latitudinal and elevational gradients of increasing fire return intervals. Examination of regional-scale output suggests that, although the tree-based models explain more of the variation in the original data, the regression models are less likely to produce extrapolation errors. Thus, the models serve complementary purposes in elucidating the relationships among fire frequency, the predictor variables, and spatial scale. The models can provide local managers with quantitative information and provide data to initialize coarse-scale fire-effects models, although predictions for individual sites should be treated with caution because of the varying quality and uneven spatial coverage of the fire history database. The models also demonstrate the integration of qualitative and quantitative methods when requisite data for fully quantitative models are unavailable. They can be tested by comparing new, independent fire history reconstructions against their predictions and can be continually updated, as better fire history data become available.

California, Idaho, Montana, Nevada, Oregon, Utah, ↗

HyFlood: A surrogate-model-based framework for compound coastal flooding

Compound coastal flooding is a major threat to low-lying coastal regions and is expected to intensify under future climate change projections. However, modeling the joint interaction of waves, storm surge, tides, and rainfall remains computationally demanding, limiting the development of fast and reliable forecast tools. Here we present HyFlood, a hybrid statistical-numerical downscaling framework capable of computing and mapping high-resolution compound flood hazards while substantially reducing the computational cost compared with fully process-based hydrodynamic modeling. HyFlood combines statistical sampling and selection algorithms with a cascade of reduced-complexity surrogate models that emulate nearshore wave transformation, surf-zone hydrodynamics, and coastal, fluvial, and pluvial flooding. The surrogate models employ machine-learning and regression algorithms applied to a low-dimensional representation of the flooding outputs, obtained through statistical dimensionality reduction. The framework is demonstrated in southern O'ahu, Hawai'i, a region exposed to elevated sea levels driven by tides, waves, and storm surge along with frequent precipitation-driven flash flooding. Validation of the surrogates against the physics-based model outputs demonstrates that HyFlood accurately reproduces daily maxima of spatially distributed flooding depths. This hybrid approach offers a scalable and efficient tool to better quantify how changes in flooding drivers translate into hazard and impact assessments, and to support compound-flood risk assessments and climate-change adaptation planning.

Hawaii↗

Stochastic analysis of solute arrival time in heterogeneous porous media

Longitudinal advective solute movement in heterogeneous porous media is investigated by considering the solute arrival time at a plane perpendicular to the mean fluid velocity. The moments of the solute arrival time are defined in terms of the stochastic properties of a statistically anisotropic hydraulic conductivity field. The flux‐averaged concentration is specified by introducing the moments of the arrival time into a probability density function for the arrival time. The quadratic dependence of the arrival time variance on position in the vicinity of an injection point is indicative of a nondiffusive process. For the assumed spatial correlation of the hydraulic conductivity, the variance of the arrival time asymptotically approaches a linear dependence on the position from the injection point similar to a diffusion process. The impact of assuming solute movement to be a diffusive process from the onset of the solute injection causes erroneous estimates of flux‐averaged concentrations at distances from the injection point that are of the order of the correlation length of the hydraulic conductivity. The arrival time analysis and the particle position analysis given in Dagan (1982, 1984) are complementary interpretations of advective solute movement that yield different definitions of the solute concentration; the position analysis intrinsically defines the resident or volume‐averaged concentration, while the flux‐averaged concentration is defined from the arrival time analysis. The temporal variation of the resident and flux‐averaged concentration are similar at a given position for small values of the variance of the hydraulic conductivity, or at distances from the solute injection point that are large relative to correlation length of the hydraulic conductivity.

Water Resources Research↗

Streamflow Simulations and Percolation Estimates Using the Soil and Water Assessment Tool for Selected Basins in North-Central Nebraska, 1940-2005

The U.S. Geological Survey, in cooperation with the Upper Elkhorn, Lower Elkhorn, Upper Loup, Lower Loup, Middle Niobrara, Lower Niobrara, Lewis and Clark, and Lower Platte North Natural Resources Districts, used the Soil and Water Assessment Tool to simulate streamflow and estimate percolation in north-central Nebraska to aid development of long-term strategies for management of hydrologically connected ground and surface water. Although groundwater models adequately simulate subsurface hydrologic processes, they often are not designed to simulate the hydrologically complex processes occurring at or near the land surface. The use of watershed models such as the Soil and Water Assessment Tool, which are designed specifically to simulate surface and near-subsurface processes, can provide helpful insight into the effects of surface-water hydrology on the groundwater system. The Soil and Water Assessment Tool was calibrated for five stream basins in the Elkhorn-Loup Groundwater Model study area in north-central Nebraska to obtain spatially variable estimates of percolation. Six watershed models were calibrated to recorded streamflow in each subbasin by modifying the adjustment parameters. The calibrated parameter sets were then used to simulate a validation period; the validation period was half of the total streamflow period of record with a minimum requirement of 10 years. If the statistical and water-balance results for the validation period were similar to those for the calibration period, a model was considered satisfactory. Statistical measures of each watershed model's performance were variable. These objective measures included the Nash-Sutcliffe measure of efficiency, the ratio of the root-mean-square error to the standard deviation of the measured data, and an estimate of bias. The model met performance criteria for the bias statistic, but failed to meet statistical adequacy criteria for the other two performance measures when evaluated at a monthly time step. A primary cause of the poor model validation results was the inability of the model to reproduce the sustained base flow and streamflow response to precipitation that was observed in the Sand Hills region. The watershed models also were evaluated based on how well they conformed to the annual mass balance (precipitation equals the sum of evapotranspiration, streamflow/runoff, and deep percolation). The model was able to adequately simulate annual values of evapotranspiration, runoff, and precipitation in comparison to reported values, which indicates the model may provide reasonable estimates of annual percolation. Mean annual percolation estimated by the model as basin averages varied within the study area from a maximum of 12.9 inches in the Loup River Basin to a minimum of 1.5 inches in the Shell Creek Basin. Percolation also varied within the studied basins; basin headwaters tended to have greater percolation rates than downstream areas. This variance in percolation rates was mainly was because of the predominance of sandy, highly permeable soils in the upstream areas of the modeled basins.

Scientific Investigations Report↗

Nutrient and sediment concentrations, yields, and loads in impaired streams and rivers in the Taunton River Basin, Massachusetts, 1997-2008

Rapid development, population growth, and the changes in land and water use accompanying development are placing increasing stress on water resources in the Taunton River Basin. An assessment by the Massachusetts Department of Environmental Protection determined that a number of tributary streams to the Taunton River are impaired for a variety of beneficial uses because of nutrient enrichment. Most of the impaired reaches are in the Matfield River drainage area in the vicinity of the City of Brockton. In addition to impairments of stream reaches in the basin, discharge of nutrient-rich water from the Taunton River contributes to eutrophication of Mount Hope and Narragansett Bays. To assess water quality and loading in the impaired tributary stream reaches in the basin, the U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection compiled existing water-quality data from previous studies for the period 1997-2006, developed and calibrated a Hydrological Simulation Program-FORTRAN (HSPF) precipitation-runoff model to simulate streamflow in areas of the basin that contain the impaired reaches for the same time period, and collected additional streamflow and water-quality data from sites on the Matfield and Taunton Rivers in 2008. A majority of the waterquality samples used in the study were collected between 1999 and 2006. Overall, the concentration, yield, and load data presented in this report represent water-quality conditions in the basin for the period 1997-2008. Water-quality data from 52 unique sites were used in the study. Most of the samples from previous studies were collected between June and September under dry weather conditions. Simulated or measured daily mean streamflow and water-quality data were used to estimate constituent yields and loads in the impaired tributary stream reaches and the main stem of the Taunton River and to develop yield-duration plots for reaches with sufficient water-quality data. Total phosphorus concentrations in the impaired-reach areas ranged from 0.0046 to 0.91 milligrams per liter (mg/L) in individual samples (number of samples (n)=331), with a median of 0.090 mg/L; total nitrogen concentrations ranged from 0.34 to 14 mg/L in individual samples (n=139), with a median of 1.35 mg/L; and total suspended solids concentrations ranged from 2/d) for total phosphorus and 100 lb/mi 2 /d for total nitrogen in these reaches. In most of the impaired reaches not affected by the Brockton Advanced Water Reclamation Facility outfall, yields were lower than in reaches downstream from the outfall, and the difference between measured and threshold yields was fairly uniform over a wide range of flows, suggesting that multiple processes contribute to nonpoint loading in these reaches. The Northeast and Mid-Atlantic SPAtially-Referenced Regression On Watershed (SPARROW) models for total phosphorus and total nitrogen also were used to estimate annual nutrient loads in the impaired tributary stream reaches and main stem of the Taunton River and predict the distribution of these loads among point and diffuse sources in reach drainage areas. SPARROW is a regional, statistical model that relates nutrient loads in streams to upstream sources and land-use characteristics and can be used to make predictions for streams that do not have nutrient-load data. The model predicts mean annual loads based on longterm streamflow and water-quality data and nutrient source conditions for the year 2002. Predicted mean annual nutrient loads from the SPARROW models were consistent with the measured yield and load data from sampling sites in the basin. For conditions in 2002, the Brockton Advanced Water Reclamation Facility outfall accounted for over 75 percent of the total nitrogen load and over 93 percent of the total phosphorus load in the Salisbury Plain and Matfield Rivers downstream from the outfall. Municipal point sources also accounted for most of the load in the main stem of the Taunton River. Multiple municipal wastewater discharges in the basin accounted for about 76 and 46 percent of the delivered loads of total phosphorus and total nitrogen, respectively, to Mount Hope Bay. For similarly sized watersheds, total delivered loads were lower in watersheds without point sources compared to those with point sources, and sources associated with developed land accounted for most of the delivered phosphorus and nitrogen loads to the impaired reaches. The concentration, yield, and load data evaluated in this study may not be representative of current (2012) point-source loading in the basin; in particular, most of the water-quality data used in the study (1999-2006) were collected prior to completion of upgrades to the Brockton Advanced Water Reclamation Facility that reduced total phosphorus and nitrogen concentrations in treated effluent. Effluent concentration data indicate that, for a given flow rate, effluent loads of total phosphorus and total nitrogen declined by about 80 and 30 percent, respectively, between the late 1990s and 2008 in response to plant upgrades. Consequently, current (2012) water-quality conditions in the impaired reaches downstream from the facility likely have improved compared to conditions described in the report.

Massachusetts↗

Exploring the regional dynamics of U.S. irrigated agriculture from 2002 to 2017

The United States has a geographically mature and stable land use and land cover system including land used as irrigated cropland; however, changes in irrigation land use frequently occur related to various drivers. We applied a consistent methodology at a 250 m spatial resolution across the lower 48 states to map and estimate irrigation dynamics for four map eras (2002, 2007, 2012, and 2017) and over four 5-year mapping intervals. The resulting geospatial maps (called the Moderate Resolution Imaging Spectroradiometer (MODIS) Irrigated Agriculture Dataset or MIrAD-US) involved inputs from county-level irrigated statistics from the U.S. Department of Agriculture, National Agricultural Statistics Service, agricultural land cover from the U.S. Geological Survey National Land Cover Database, and an annual peak vegetation index derived from expedited MODIS satellite imagery. This study investigated regional and periodic patterns in the amount of change in irrigated agriculture and linked gains and losses to proximal causes and consequences. While there was a 7% overall increase in irrigated area from 2002 to 2017, we found surprising variability by region and by 5-year map interval. Irrigation land use dynamics affect the environment, water use, and crop yields. Regionally, we found that the watersheds with the largest irrigation gains (based on percent of area) included the Missouri, Upper Mississippi, and Lower Mississippi watersheds. Conversely, the California and the Texas–Gulf watersheds experienced fairly consistent irrigation losses during these mapping intervals. Various drivers for irrigation dynamics included regional climate fluctuations and drought events, demand for certain crops, government land or water policies, and economic incentives like crop pricing and land values. The MIrAD-US (Version 4) was assessed for accuracy using a variety of existing regionally based reference data. Accuracy ranged between 70% and 95%, depending on the region.

Conterminous United States↗

Science summary in support of Manatee Protection Area (MPA) design in Puerto Rico

The Antillean manatee (Trichechus manatus manatus), a subspecies of the West Indian manatee, is listed as endangered by the US Department of Interior. In accordance with its listing, the U.S. Fish and Wildlife Service’s Caribbean Field Office (USFWS) is mandated to create one or more Manatee Protection Areas (MPAs) for Puerto Rico. Designation of these areas must comply with the legal definition of an MPA’s purpose: to prevent or reduce take of manatees (CFR 50: 44 FR 60964, Oct. 22, 1979). To meet this goal, we pursued two objectives: 1) identify areas which include the specific ecological attributes necessary to support manatee populations, and 2) identify areas where take can be reduced through approved MPA regulatory frameworks. We achieved these objectives through literature review, expert elicitation, and geospatial modeling. This report delivers to USFWS a set of nine potential MPA regions. These regions represent the spatial realization of experts’ hypotheses regarding manatee requirements and threats, and the potential to implement MPA strategies (e.g. watercraft access, speed regulations, signage and boater education). The nine regions are compared based on a number of factors, including their potential to reduce take, quality of the habitat encompassed, and total area. These maps and statistics serve as suitable starting points to select one or more MPA sites, but we recommend that the mapped attributes and threats (i.e., boating activity) of MPAs be ground-truthed to visually confirm the local presence of resources, threats, and manatees before any area is selected. Once established, the effectiveness of MPAs can be monitored and updated through processes of adaptive monitoring and management. Aerial surveys, radio tracking studies, and public surveys are all valuable tools to assess the success of an MPA. Establishing MPAs is a management action that, integrated within the species Recovery Plan, should enhance the conservation of manatees.

Cooperator Science Series↗

Characterization and relation of precipitation, streamflow, and water-quality data at the U.S. Army Garrison Fort Carson and Piñon Canyon Maneuver Site, Colorado, water years 2013–14

To evaluate the influence of military training activities on streamflow and water quality, the U.S. Geological Survey, in cooperation with the U.S. Department of the Army, began a hydrologic data collection network on the U.S. Army Garrison Fort Carson in 1978 and on the Piñon Canyon Maneuver Site in 1983. This report is a summary and characterization of the precipitation, streamflow, and water-quality data collected at 43 sites between October 1, 2012, and September 30, 2014 (water years 2013 and 2014). Variations in the frequency of daily precipitation, seasonal distribution, and seasonal and annual precipitation at 5 stations at the U.S. Army Garrison Fort Carson and 18 stations at or near the Piñon Canyon Maneuver Site were evaluated. Isohyetal diagrams indicated a general pattern of increase in total annual precipitation from east to west at the U.S. Army Garrison Fort Carson and the Piñon Canyon Maneuver Site. Between about 54 and 79 percent of daily precipitation was 0.1 inch or less in magnitude. Precipitation events were larger and more frequent between July and September. Daily streamflow data from 16 sites were used to evaluate temporal and spatial variations in streamflow for the water years 2013 and 2014. At all sites, median daily mean streamflow for the 2-year period ranged from 0.0 to 9.60 cubic feet per second. Daily mean streamflow hydrographs are included in this report. Five sites on the Piñon Canyon Maneuver Site were monitored for peak stage using crest-stage gages. At the Piñon Canyon Maneuver Site, five sites had a stage recorder and precipitation gage, providing a paired streamflow-precipitation dataset. There was a statistically significant correlation between precipitation and streamflow based on Spearman’s rho correlation (rho values ranged from 0.17 to 0.35). Suspended-sediment samples were collected in April through October for water years 2013–14 at one site at the U.S. Army Garrison Fort Carson and five sites at the Piñon Canyon Maneuver Site. Suspended-sediment-transport curves were used to illustrate the relation between streamflow and suspended-sediment concentration. All these sediment-transport curves showed a streamflow dependent suspended-sediment concentration relation except for the U.S. Geological Survey station Bent Canyon Creek at mouth near Timpas, CO. Water-quality data were collected and reported from seven sites on the U.S. Army Garrison Fort Carson and the Piñon Canyon Maneuver Site during water years 2013–14. Sample results exceeding an established water-quality standard were identified. Selected water-quality properties and constituents were stratified to compare spatial variation among selected characteristics using boxplots. Trilinear diagrams were used to classify water type based on ionic concentrations of water-quality samples collected during the study period. At the U.S. Army Garrison Fort Carson and the Piñon Canyon Maneuver Site, 27 samples were classified as very hard or brackish. Seven samples had a lower hardness character relative to the other samples. Four of those nine samples were collected at two U.S. Geological Survey stations (Turkey Creek near Fountain, CO, and Little Fountain Creek above Highway 115 at Fort Carson, CO), which have different geologic makeup. Three samples collected at the Piñon Canyon Maneuver Site had a markedly lower hardness likely because of dilution from an increase in streamflow.

Colorado↗

Methodology and application of combined watershed and ground-water models in Kansas

Increased irrigation in Kansas and other regions during the last several decades has caused serious water depletion, making the development of comprehensive strategies and tools to resolve such problems increasingly important. This paper makes the case for an intermediate complexity, quasi-distributed, comprehensive, large-watershed model, which falls between the fully distributed, physically based hydrological modeling system of the type of the SHE model and the lumped, conceptual rainfall-runoff modeling system of the type of the Stanford watershed model. This is achieved by integrating the quasi-distributed watershed model SWAT with the fully-distributed ground-water model MODFLOW. The advantage of this approach is the appreciably smaller input data requirements and the use of readily available data (compared to the fully distributed, physically based models), the statistical handling of watershed heterogeneities by employing the hydrologic-response-unit concept, and the significantly increased flexibility in handling stream-aquifer interactions, distributed well withdrawals, and multiple land uses. The mechanics of integrating the component watershed and ground-water models are outlined, and three real-world management applications of the integrated model from Kansas are briefly presented. Three different aspects of the integrated model are emphasized: (1) management applications of a Decision Support System for the integrated model (Rattlesnake Creek subbasin); (2) alternative conceptual models of spatial heterogeneity related to the presence or absence of an underlying aquifer with shallow or deep water table (Lower Republican River basin); and (3) the general nature of the integrated model linkage by employing a watershed simulator other than SWAT (Wet Walnut Creek basin). These applications demonstrate the practicality and versatility of this relatively simple and conceptually clear approach, making public acceptance of the integrated watershed modeling system much easier. This approach also enhances model calibration and thus the reliability of model results. (C) 2000 Elsevier Science B.V.Increased irrigation in Kansas and other regions during the last several decades has caused serious water depletion, making the development of comprehensive strategies and tools to resolve such problems increasingly important. This paper makes the case for an intermediate complexity, quasi-distributed, comprehensive, large-watershed model, which falls between the fully distributed, physically based hydrological modeling system of the type of the SHE model and the lumped, conceptual rainfall-runoff modeling system of the type of the Stanford watershed model. This is achieved by integrating the quasi-distributed watershed model SWAT with the fully-distributed ground-water model MODFLOW. The advantage of this approach is the appreciably smaller input data requirements and the use of readily available data (compared to the fully distributed, physically based models), the statistical handling of watershed heterogeneities by employing the hydrologic-response-unit concept, and the significantly increased flexibility in handling stream-aquifer interactions, distributed well withdrawals, and multiple land uses. The mechanics of integrating the component watershed and ground-water models are outlined, and three real-world management applications of the integrated model from Kansas are briefly presented. Three different aspects of the integrated model are emphasized: (1) management applications of a Decision Support System for the integrated model (Rattlesnake Creek subbasin); (2) alternative conceptual models of spatial heterogeneity related to the presence or absence of an underlying aquifer with shallow or deep water table (Lower Republican River basin); and (3) the general nature of the integrated model linkage by employing a watershed simulator other than SWAT (Wet Walnut Creek basin). These applications demonstrate the practicality and ve

Journal of Hydrology↗

Trends in publications in fluvial geomorphology over two decades: A truly new era in the discipline owing to recent technological revolution?

Trends in the field of fluvial geomorphology have been reviewed by a number of authors, who have emphasized the dramatic change occuring in the field in the last two decades of the twentieth century, largely as a result of technological advances. Nevertheless, no prior authors have systematically compiled data on publications in fluvial geomorphology over a long period and statistically analyzed the resulting data set. In this contribution we present a quantitative analysis of fluvial geomorphology papers published in the twenty-two-year period 1987–2009 in five journals of the discipline with a more specific focus on Geomorphology and Earth Surface Processes and Landforms ( ESPL ), identifying authorships, geographic origin of authors, and spatial and temporal scales covered. We also documented the tools employed, demonstrating the transformation of the field with the emergence of new tools over this period, and conducted a cluster to highlight links between tools and a set of factors (country of author's origin, journals, time, and spatial and temporal scales). Of the 1717 papers published in the five journals during this period, the results showed an increased diversity in the nationality of the first author, mainly when dealing with present time scale, and channel feature. Our data show a significant change in methods used in the field as a result of the increase in data availability and new sources of information from remote sensing (ground, airborne and, satellite). Clearly, a new era in knowledge production is observed since 2000, showing the emergence of a second period of active quantification and an internationalization of the fields.

Geomorphology↗

A universal approximation to grain size from images of non-cohesive sediment

The two-dimensional spectral decomposition of an image of sediment provides a direct statistical estimate, grid-by-number style, of the mean of all intermediate axes of all single particles within the image. We develop and test this new method which, unlike existing techniques, requires neither image processing algorithms for detection and measurement of individual grains, nor calibration. The only information required of the operator is the spatial resolution of the image. The method is tested with images of bed sediment from nine different sedimentary environments (five beaches, three rivers, and one continental shelf), across the range 0.1 mm to 150 mm, taken in air and underwater. Each population was photographed using a different camera and lighting conditions. We term it a “universal approximation” because it has produced accurate estimates for all populations we have tested it with, without calibration. We use three approaches (theory, computational experiments, and physical experiments) to both understand and explore the sensitivities and limits of this new method. Based on 443 samples, the root-mean-squared (RMS) error between size estimates from the new method and known mean grain size (obtained from point counts on the image) was found to be ±≈16%, with a 95% probability of estimates within ±31% of the true mean grain size (measured in a linear scale). The RMS error reduces to ≈11%, with a 95% probability of estimates within ±20% of the true mean grain size if point counts from a few images are used to correct bias for a specific population of sediment images. It thus appears it is transferable between sedimentary populations with different grain size, but factors such as particle shape and packing may introduce bias which may need to be calibrated for. For the first time, an attempt has been made to mathematically relate the spatial distribution of pixel intensity within the image of sediment to the grain size.

Journal of Geophysical Research F: Earth Surface↗

Comparison of temperature, specific conductance, pH, and dissolved oxygen at selected basic fixed sites in south-central Texas, 1996-98

One component of the surface-water part of the U.S. Geological Survey National Water-Quality Assessment Program is the use of continuous water-quality monitors to help characterize the spatial and temporal distribution of general water quality in relation to hydrologic conditions. During 1996-98, six continuous water-quality monitors in the South-Central Texas study unit collected water temperature, specific conductance, pH, and dissolved oxygen data. The data were compared among the six sites using boxplots of monthly mean values, summary statistics of monthly values, and hydrographs of daily mean values.

Texas↗

Correlation analysis of a ground-water level monitoring network, Miami-Dade County, Florida

The U.S. Geological Survey cooperative ground-water monitoring program in Miami-Dade County, Florida, expanded from 4 to 98 continuously recording water-level monitoring wells during the 1939-2001 period. Network design was based on area specific assessments; however, no countywide statistical assessments of network coverage had been performed for the purpose of assessing network redundancy. To aid in the assessment of network redundancy, correlation analyses were performed using S-PLUS 2000 statistical analysis software for daily maximum water-level data from 98 monitoring wells for the November 1, 1973, to October 31, 2000 period. Because of the complexities of the hydrologic, water-supply, and water-management systems in Miami-Dade County and the changes that have occurred to these systems through time, spatial and temporal variations in the degree of correlation had to be considered. To assess temporal variation in correlation, water-level data from each well were subdivided by year and by wet and dry seasons. For each well, year, and season, correlation analyses were performed on the data from those wells that had available data. For selected wells, the resulting correlation coefficients from each year and season were plotted with respect to time. To assess spatial variation in correlation, the coefficients determined from the correlation analysis were averaged. These average wet- and dry-season correlation coefficients were plotted spatially using geographic information system software. Wells with water-level data that correlated with a coefficient of 0.95 or greater were almost always located in relatively close proximity to each other. Five areas were identified where the water-level data from wells within the area remained correlated with that of other wells in the area during the wet and dry seasons. These areas are located in or near the C-1 and C-102 basins (2 wells), in or near the C-6 and C-7 basins (2 wells), near the Florida Keys Aqueduct Authority Well Field (2 wells), near the Hialeah-Miami Springs Well Field (6 wells), and near the West Well Field (21 wells). Data from the remaining 65 wells (most of the wells in the network) generally were not correlated with those of other wells during both the wet and dry seasons with an average coefficient of 0.95 or greater for the comparison. Because many of the wells near the West Well Field and some near the Hialeah-Miami Springs Well Field had not been in operation for very long (most having been installed in 1994), the averaged correlation coefficients for these wells were often determined using only a few seasons of data. For the few instances where water-level data were found to be well correlated on average for a lengthy period of record, short-term declines in correlation were often identified. In general, it would be beneficial to compare data for longer periods of record than currently available.

Open-File Report↗