Search USGSSearch

SEARCH · Search USGS

Results for “Computer contribution”

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 307 records · Page 17Linked to original sources

Advective transport observations with MODPATH-OBS--documentation of the MODPATH observation process

The MODPATH-OBS computer program described in this report is designed to calculate simulated equivalents for observations related to advective groundwater transport that can be represented in a quantitative way by using simulated particle-tracking data. The simulated equivalents supported by MODPATH-OBS are (1) distance from a source location at a defined time, or proximity to an observed location; (2) time of travel from an initial location to defined locations, areas, or volumes of the simulated system; (3) concentrations used to simulate groundwater age; and (4) percentages of water derived from contributing source areas. Although particle tracking only simulates the advective component of conservative transport, effects of non-conservative processes such as retardation can be approximated through manipulation of the effective-porosity value used to calculate velocity based on the properties of selected conservative tracers. This program can also account for simple decay or production, but it cannot account for diffusion. Dispersion can be represented through direct simulation of subsurface heterogeneity and the use of many particles. MODPATH-OBS acts as a postprocessor to MODPATH, so that the sequence of model runs generally required is MODFLOW, MODPATH, and MODPATH-OBS. The version of MODFLOW and MODPATH that support the version of MODPATH-OBS presented in this report are MODFLOW-2005 or MODFLOW-LGR, and MODPATH-LGR. MODFLOW-LGR is derived from MODFLOW-2005, MODPATH 5, and MODPATH 6 and supports local grid refinement. MODPATH-LGR is derived from MODPATH 5. It supports the forward and backward tracking of particles through locally refined grids and provides the output needed for MODPATH_OBS. For a single grid and no observations, MODPATH-LGR results are equivalent to MODPATH 5. MODPATH-LGR and MODPATH-OBS simulations can use nearly all of the capabilities of MODFLOW-2005 and MODFLOW-LGR; for example, simulations may be steady-state, transient, or a combination. Though the program name MODPATH-OBS specifically refers to observations, the program also can be used to calculate model prediction of observations. MODPATH-OBS is primarily intended for use with separate programs that conduct sensitivity analysis, data needs assessment, parameter estimation, and uncertainty analysis, such as UCODE_2005, and PEST. In many circumstances, refined grids in selected parts of a model are important to simulated hydraulics, detailed inflows and outflows, or other system characteristics. MODFLOW-LGR and MODPATH-LGR support accurate local grid refinement in which both mass (flows) and energy (head) are conserved across the local grid boundary. MODPATH-OBS is designed to take advantage of these capabilities. For example, particles tracked between a pumping well and a nearby stream, which are simulated poorly if a river and well are located in a single large grid cell, can be simulated with improved accuracy using a locally refined grid in MODFLOW-LGR, MODPATH-LGR, and MODPATH-OBS. The locally-refined-grid approach can provide more accurate simulated equivalents to observed transport between the well and the river. The documentation presented here includes a brief discussion of previous work, description of the methods, and detailed descriptions of the required input files and how the output files are typically used.

Techniques and Methods

Models for inference in dynamic metacommunity systems

A variety of processes are thought to be involved in the formation and dynamics of species assemblages. For example, various metacommunity theories are based on differences in the relative contributions of dispersal of species among local communities and interactions of species within local communities. Interestingly, metacommunity theories continue to be advanced without much empirical validation. Part of the problem is that statistical models used to analyze typical survey data either fail to specify ecological processes with sufficient complexity or they fail to account for errors in detection of species during sampling. In this paper, we describe a statistical modeling framework for the analysis of metacommunity dynamics that is based on the idea of adopting a unified approach, multispecies occupancy modeling, for computing inferences about individual species, local communities of species, or the entire metacommunity of species. This approach accounts for errors in detection of species during sampling and also allows different metacommunity paradigms to be specified in terms of species‐ and location‐specific probabilities of occurrence, extinction, and colonization: all of which are estimable. In addition, this approach can be used to address inference problems that arise in conservation ecology, such as predicting temporal and spatial changes in biodiversity for use in making conservation decisions. To illustrate, we estimate changes in species composition associated with the species‐specific phenologies of flight patterns of butterflies in Switzerland for the purpose of estimating regional differences in biodiversity.

Ecology

Streamflow characteristics and trends at selected streamgages in southwest and south-central Kansas

Historical data for nine selected streamgages in southwest and south-central Kansas were used in an assessment of streamflow characteristics and trends. This information is required by the U.S. Fish and Wildlife Service and the Kansas Department of Wildlife, Parks and Tourism to assist with the effective management of Etheostoma cragini (Arkansas darter) habitats and populations in the State. Changing streamflow conditions, such as a reduction or elimination of streamflow, may adversely affect the Arkansas darter. Priority basins for the Arkansas darter represented by the selected streamgages include the Cimarron River, Rattlesnake Creek, the North Fork Ninnescah River, the South Fork Ninnescah River, the Medicine Lodge River, and the Chikaskia River. Streamflow conditions were assessed using annual streamflow characteristics computed for the period of record for each of the selected streamgages. Specific streamflow characteristics computed were mean discharge, mean base flow, 90th-percentile flow, 10th-percentile flow, minimum 7-day mean flow, minimum 28-day mean flow, number of days of flow less than 1 cubic foot per second, and number of zero-flow days. Two of the priority basins had statistically significant decreases in annual mean discharge during the period of record. In the Cimarron River Basin, there was a pronounced multidecadal decrease in the magnitude and variability of annual mean discharge. Concurrently, the percentage of the annual mean discharge that was contributed by base flow increased. In the Rattlesnake Creek Basin, there was a pre-1985 decrease in annual mean discharge. Typically, in these two basins, significant decreases were indicated for mean base flow, 90th-percentile flow, 10th-percentile flow, minimum 7-day mean flow, and minimum 28-day mean flow. No significant trend in annual mean discharge was indicated for the North Fork Ninnescah, South Fork Ninnescah, Medicine Lodge, and Chikaskia River Basins. For the Medicine Lodge and Chikaskia River Basins as well as the downstream part of the South Fork Ninnescah River Basin, a significant increase in mean base flow and 10th-percentile flow was indicated. Also, for the latter two basins, a significant increase was indicated for minimum 7-day mean flow. Factors investigated to explain long-term trends in annual mean discharge, or lack thereof, included precipitation and groundwater withdrawals. Annual precipitation in the study area varied substantially from 1951 to 2013 with no pronounced long-term trend. Thus, a precipitation-related explanation for the significant decrease in annual mean discharge in the Cimarron River and Rattlesnake Creek Basins was not supported. Because the most pronounced decreases in annual mean discharge were in the basin with the largest groundwater-level declines (that is, the Cimarron River Basin), both in terms of magnitude and areal extent, it is likely that groundwater withdrawals were a primary, if not dominant, causative factor. The occurrence of extremely low-flow (less than 1 cubic foot per second) and zero-flow days varied by basin and year. Typically, such days occurred in the summer and autumn for all basins.

Kansas

Simulation of fluid, heat transport to estimate desert stream infiltration

In semiarid regions, the contribution of infiltration from intermittent streamflow to ground water recharge may be quantified by comparing simulations of fluid and heat transport beneath stream channels to observed ground temperatures. In addition to quantifying natural recharge, streamflow infiltration estimates provide a means to characterize the physical properties of stream channel sediments and to identify suitable locations for artificial recharge sites. Rates of winter streamflow infiltration along stream channels are estimated based on the cooling effect of infiltrated water on streambed sediments, combined with the simulation of two-dimensional fluid and heat transport using the computer program VS2DH. The cooling effect of ground water is determined by measuring ground temperatures at regular intervals beneath stream channels and nearby channel banks in order to calculate temperature-depth profiles. Additional data inputs included the physical, hydraulic, and thermal properties of unsaturated alluvium, and monthly ground temperatures measurements over an annual cycle. Observed temperatures and simulation results can provide estimates of the minimum threshold for deep infiltration, the variability of infiltration along stream channels, and also the frequency of infiltration events.

Ground Water

South Cascade (USA/North Cascades)

The U.S. Geological Survey has closely monitored this temperate mountain glacier since the late 1950s. During 1958-2007, the glacier retreated about 0.7 km and shrank in area from 2.71 to 1.73 km 2 , although part of the area change was due to separation of contributing ice bodies from the main glacier. Maximum and average glacier thicknesses are about 170 and 80 m, respectively. Year-to-year variations of snow accumulation amounts on the glacier are largely attributable to the regional maritime climate and fluctuating climate conditions of the North Pacific Ocean. Long-term-average precipitation is about 4500 mm and most of that falls as snow during October through May. Average annual air temperature at 1,900 m altitude (the approximate ELA 0 ) was estimated to be 1.6°C during 2000-2009. Mass balances are computed yearly by the direct glaciological method. Mass balances measured at selected locations are used in an interpolation and extrapolation procedure that computes the mass balance at each point in the glacier surface altitude grid. The resulting mass balance grid is averaged to obtain glacier mass balances. Additionally, the geodetic method has been applied to compute glacier net balances in 1970, 1975, 1977, 1979-80, and 1985-97. Winter snow accumulation on the glacier during 2007/08 and 2008/09 was larger than the long-term (1959-2009) average. The 2007/08 preliminary summer balance (-3510 mm w.e.) was slightly more negative than the long-term average and this yielded a preliminary 2007/08 net balance (-290 mm w.e.), which was less negative than the average for the period of record (-600 mm w.e.). Summer 2009 was uncommonly warm and the preliminary 2008/09 summer balance (-4980 mm w.e.) was more negative than any on record for the glacier. The 2008/09 glacier net balance (-1860 mm w.e.) was among the 10 most negative for the period of net balance record (1953-2009). Material presented here is preliminary in nature and presented prior to final review. These data and information are provided with the understanding that they are not guaranteed to be correct or complete. Users are cautioned to consider carefully the provisional nature of these data and information before using them for decisions that concern personal or public safety or the conduct of business that involves substantial monetary or operational consequences. Conclusions drawn from, or actions undertaken on the basis of, such data and information are the sole responsibility of the user.

Washington

New composite spectra of Mars, 0.4-5.7 μm

About 15 areas were observed in the equatorial regions of Mars by the infrared spectrometers IRS (Mariner 6 and 7) and ISM (Phobos-2). The comparison between the spectra shows a remarkable consistency between two data sets acquired 20 years apart and calibrated independently. This similarity demonstrates the accuracy of ISM calibration above 2 μm, except for a possible stray light contribution above 2.6 μm, on the order of ∼1–2% of the solar flux at 2.7 μm. Most differences in spectral shapes are related to differences in spectral/spatial resolution and viewing geometries. No important variation in surface properties is detected, except for a spot in southern Arabia Terra which has a much deeper hydration feature in IRS spectra; differences in viewing geometries and spatial resolutions do not seem to account for this difference that could result from shifting or dehydration of surface materials. Composite spectra of several types of bright and dark materials are computed by modeling the thermal emission and are completed with telescopic spectra in the visible range. Modeled reflectance in the 3.0–5.7 μm range is consistent with basalts and palagonites. The bright regions and analog palagonite spectra are different from hematite in this range, but resemble several phyllosilicates. We infer that (1) although hematite dominates the spectra in the 0.4- to 2.5-μm range, the silicate-clay host is spectrally active beyond 3 μm and can be identified from this domain; (2) phyllosilicates such as montmorillonite or smectite may be abundant components of the martian soils, although the domain below 3 μm lacks the characteristic features of the most usual terrestrial clay minerals.

Icarus

Flow monitoring along the western Tamiami Trail between County Road 92 and State Road 29 in support of the Comprehensive Everglades Restoration Plan, 2007-2010

The construction of U.S. Highway 41 (Tamiami Trail), the Southern Golden Gate Estates development, and the Barron River Canal has altered the flow of freshwater to the Ten Thousand Islands estuary of Southwest Florida. Two restoration projects, the Picayune Strand Restoration Project and the Tamiami Trail Culverts Project, both associated with the Comprehensive Everglades Restoration Plan, were initiated to address this issue. Quantifying the flow of freshwater to the estuary is essential to assessing the effectiveness of these projects. The U.S. Geological Survey conducted a study between March 2006 and September 2010 to quantify the freshwater flowing under theTamiami Trail between County Road 92 and State Road 29 in southwest Florida, excluding the Faka Union Canal (which is monitored by South Florida Water Management District). The study period was after the completion of the Tamiami Trail Culverts Project and prior to most of the construction related to the Picayune Restoration Project. The section of the Tamiami Trail that was studied contains too many structures (35 bridges and 16 culverts) to cost-effectively measure each structure on a continuous basis, so the area was divided into seven subbasins. One bridge within each of the subbasins was instrumented with an acoustic Doppler velocity meter. The index velocity method was used to compute discharge at the seven instrumented bridges. Periodic discharge measurements were made at all structures, using acoustic Doppler current profilers at bridges and acoustic Doppler velocity meters at culverts. Continuous daily mean values of discharge for the uninstrumented structures were calculated on the basis of relations between the measured discharge at the uninstrumented stations and the discharge and stage at the instrumented bridge. Estimates of daily mean discharge are available beginning in 2006 or 2007 through September 2010 for all structures. Subbasin comparison is limited to water years 2008–2010. The Faka Union Canal contributed more than half (on average 60 percent) of the flow under the Tamiami Trail between State Road 29 and County Road 92 during water years 2008–2010. During water years 2008–2010, an average 9 percent of the flow through the study area came from west of the Faka Union Canal and an average 31 percent came from east of the Faka Union Canal. Flow data provided by this study serve as baseline information about the seasonal and spatial distribution of freshwater flow under the Tamiami Trail between County Road 92 and State Road 29, and study results provide data to evaluate restoration efforts.

Florida

Regional equations for estimation of peak-streamflow frequency for natural basins in Texas

Peak-streamflow frequency estimates are needed for flood-plain management; for objective assessment of flood risk; and for cost-effective design of dams, levees, other flood-control structures, roads, bridges, and culverts. Peak-streamflow frequency represents the peak discharges for recurrence intervals of 2, 5, 10, 25, 50, and 100 years. The U.S. Geological Survey, in cooperation with the Texas Department of Transportation, conducted an investigation to develop regional regression equations for the estimation of peak-streamflow frequency for ungaged sites in natural basins in Texas. Peak-streamflow frequency data from streamflow-gaging stations in natural drainage basins of Texas were used. Peak-streamflow frequency for 559 Texas stations with natural (unregulated and rural or nonurbanized) basins was estimated with annual peak-streamflow data through 1993. The peak-streamflow frequency and drainage-basin characteristics for the Texas stations were used to develop 16 sets of equations to estimate peak-streamflow frequency for ungaged natural stream sites in each of 11 regions in Texas. The relation between peak-streamflow frequency and contributing drainage area for 5 of the 11 regions is curvilinear, requiring that one set of equations be developed for drainage areas less than 32 square miles and another set be developed for drainage areas greater than 32 square miles. These equations, developed through multiple-regression analysis using weighted least squares, are based on the relation between peak-streamflow frequency and basin characteristics for streamflow-gaging stations. The regions represent areas with similar flood characteristics. The use and limitations of the regression equations also are discussed. Additionally, procedures are presented to compute the 50-, 67-, and 90-percent confidence limits for any estimation from the equations. Also, supplemental peak-streamflow frequency and basin characteristics for 105 selected stations bordering Texas are included in the report. This supplemental information will aid in interpretation of flood characteristics for sites near the state borders of Texas.

Texas

Water quality, hydrology, and phosphorus loading to Little St. Germain Lake, Wisconsin, with special emphasis on the effects of winter aeration and ground-water inputs

Little St. Germain Lake is a 978-acre, multibasin lake in Vilas County, Wisconsin. In the interest of protecting and improving the water quality of the lake, the Little St. Germain Lake District initiated several cooperative studies with the U.S. Geological Survey between 1991 and 2004 to (1) document the water quality and the extent of winter anoxia in the lake, (2) evaluate the success of aerators at eliminating winter anoxia, (3) develop water and nutrient budgets for the lake, and (4) assess how the water quality of the lake should respond to changes in phosphorus loading. This report presents the results of these cooperative studies with special emphasis on the water quality in the lake since 2000, including the effects of winter aeration and the importance of ground-water contributions of phosphorus to the productivity of the lake. Measurements collected during these studies indicate that the water quality in Little St. Germain Lake was consistently different among basins. The West Bay consistently had the best water quality, the South Bay had intermediate water quality, and the East and Upper East Bays consistently had the worst water quality. The water quality in each of the basins was relatively stable from 1991 to 2000; however, since 2001, the West Bay has changed from oligotrophic to mesotrophic, the South Bay has changed from mesotrophic to eutrophic, and the East and Upper East Bays have changed from eutrophic to eutrophic/hypereutrophic. Winter anoxia frequently occurred throughout most of the lake, except in the West Bay and just below the ice in the East Bay. To eliminate winter anoxia, coarse-bubble line aerators were installed and operated in the Upper East, East, and South Bays. The aerators in the Upper East and South Bays were very successful at eliminating winter anoxia; however, the aerator in the East Bay had little impact on the dissolved oxygen concentrations throughout its basin. Detailed water and phosphorus budgets computed for the lake indicated that inflow from Muskellunge Creek was the major source of phosphorus to the lake and that ground water was the secondary source. Results from a detailed ground-water-flow model indicated that ground water flows into the lake from all sides, except the south sides of the West and Second South Bays. Most of the phosphorus appears to come from natural sources, such as ground water and surface water flowing through relatively undeveloped areas surrounding Little St. Germain Lake and Muskellunge Lake. Several empirical water-quality models were used to simulate how the East and Upper East Bays of the lake should respond to reductions in phosphorus loading from Muskellunge Creek. Simulation results indicated that reductions in tributary loading could improve the water quality of the East and Upper East Bays. Improving the water quality of these bays would also improve the water quality of the South and Second South Bays because of the flow of water through the lake. However, even with phosphorus loading from Muskellunge Creek completely eliminated, most of the lake would remain borderline mesotrophic/eutrophic because of the contributions of phosphorus from ground water.

Wisconsin

Characterizing riverbed sediment using high-frequency acoustics 1: spectral properties of scattering

Bed-sediment classification using high-frequency hydro-acoustic instruments is challenging when sediments are spatially heterogeneous, which is often the case in rivers. The use of acoustic backscatter to classify sediments is an attractive alternative to analysis of topography because it is potentially sensitive to grain-scale roughness. Here, a new method is presented which uses high-frequency acoustic backscatter from multibeam sonar to classify heterogeneous riverbed sediments by type (sand, gravel,rock) continuously in space and at small spatial resolution. In this, the first of a pair of papers that examine the scattering signatures from a heterogeneous riverbed, methods are presented to construct spatially explicit maps of spectral properties from geo-referenced point clouds of geometrically and radiometrically corrected echoes. Backscatter power spectra are computed to produce scale and amplitude metrics that collectively characterize the length scales of stochastic measures of riverbed scattering, termed ‘stochastic geometries’. Backscatter aggregated over small spatial scales have spectra that obey a power-law. This apparently self-affine behavior could instead arise from morphological- and grain-scale roughnesses over multiple overlapping scales, or riverbed scattering being transitional between Rayleigh and geometric regimes. Relationships exist between stochastic geometries of backscatter and areas of rough and smooth sediments. However, no one parameter can uniquely characterize a particular substrate, nor definitively separate the relative contributions of roughness and acoustic impedance (hardness). Combinations of spectral quantities do, however, have the potential to delineate riverbed sediment patchiness, in a data-driven approach comparing backscatter with bed-sediment observations (which is the subject of part two of this manuscript).

Journal of Geophysical Research F: Earth Surface

Computer simulations of large asteroid impacts into oceanic and continental sites--Preliminary results on atmospheric, cratering and ejecta dynamics

Computer simulations have been completed that describe passage of a 10-km-diameter asteroid through the Earth's atmosphere and the subsequent cratering and ejecta dynamics caused by impact of the asteroid into both oceanic and continental sites. The asteroid was modeled as a spherical body moving vertically at 20 km/s with a kinetic energy of 2.6 × 10 30 ergs (6.2 × 10 7 Mt ). Detailed material modeling of the asteroid, ocean, crustal units, sedimentary unit, and mantle included effects of strength and fracturing, generic asteroid and rock properties, porosity, saturation, lithostatic stresses, and geothermal contributions, each selected to simulate impact and geologic conditions that were as realistic as possible. Calculation of the passage of the asteroid through a U.S. Standard Atmosphere showed development of a strong bow shock wave followed by a highly shock compressed and heated air mass. Rapid expansion of this shocked air created a large low-density region that also expanded away from the impact area. Shock temperatures in air reached ∼20, 000K near the surface of the uplifting crater rim and were as high as ∼2000K at more than 30 km range and 10 km altitude. Calculations to 30 s showed that the shock fronts in the air and in most of the expanding shocked air mass preceded the formation of the crater, ejecta, and rim uplift and did not interact with them. As cratering developed, uplifted rim and target material were ejected into the very low density, shock-heated air immediately above the forming crater, and complex interactions could be expected. Calculations of the impact events showed equally dramatic effects on the oceanic and continental targets through an interval of 120 s. Despite geologic differences in the targets, both cratering events developed comparable dynamic flow fields and by ∼29s had formed similar-sized transient craters ∼39km deep and ∼62km across. Transient-rim uplift of ocean and crust reached a maximum altitude of nearly 40 km at ∼30s and began to decay at velocities of 500 m/s to develop large-tsunami conditions. After ∼30s, strong gravitational rebound drove both craters toward broad flat-floored shapes. At 120 s, transient crater diameters were ∼80km (continental) and ∼105km (oceanic) and transient depths were ∼27km; crater floors consisting of melted and fragmented hot rock were rebounding rapidly upward. By 60 s, the continental crater had ejected ∼2 × 10 14 t, about twice the mass ejected from the oceanic crater. By 120 s, ∼70, 000km 3 (continental) and ∼90, 000km 3 (oceanic) target material were excavated (no mantle) and massive ejecta blankets were formed around the craters. We estimate that in excess of ∼70% of the ejecta would finally lie within ∼3 crater diameters of the impact, and the remaining ejecta (∼10 13 t), including the vaporized asteroid, would be ejected into the atmosphere to altitudes as high as the ionosphere. Effects of secondary volcanism and return of the ocean over hot oceanic crater floor could also be expected to contribute substantial material to the atmosphere.

International Journal of Impact Engineering

Framework for mapping liquefaction hazard–Targeted design ground motions

Liquefaction-induced ground failure poses substantial challenges to geotechnical earthquake engineering design. Current approaches for designing against liquefaction hazards, as specified in most seismic provisions, focus on estimating a liquefaction factor of safety ( 𝐹⁢𝑆𝐿 ) and typically characterize earthquake loading using design parameters based on probabilistic or deterministic ground motion levels. Because 𝐹⁢𝑆𝐿 is estimated deterministically, this basis of design neglects considerable uncertainties for estimating liquefaction triggering and its consequences and results in a lack of liquefaction-specific design criteria, particularly as structural design has advanced toward risk-targeted performance objectives. This study presents a framework for developing liquefaction-targeted design criteria based on a minimum acceptable return period of liquefaction, informed by probabilistic liquefaction hazard analysis (PLHA). PLHA quantifies annualized rates of liquefaction by considering contributions from (1) the full ground-motion probability space, and (2) uncertainties in liquefaction triggering using probabilistic models. PLHA is used in this study to characterize the current, effective return periods of 𝐹⁢𝑆𝐿 ( 𝑇𝑅,𝐹⁢𝑆 ) obtained from conventional liquefaction hazard analysis (CLHA) using uniform-hazard ground motions. 𝑇𝑅,𝐹⁢𝑆 is evaluated in a parametric study of nearly 100 sites throughout the conterminous United States. The results indicate large geographic variations in acceptable liquefaction hazard levels, with implied 𝑇𝑅,𝐹⁢𝑆 ranging between approximately 1,000 to 3,000 years. To address these inconsistencies without the computational demands of full PLHA, a framework is proposed for developing a liquefaction-targeted design peak ground acceleration, 𝑃⁢𝐺⁢𝐴𝐿 , for use in liquefaction models that result in consistent liquefaction design levels across all geographic locations. The mapped 𝑃⁢𝐺⁢𝐴𝐿 is shown to be somewhat sensitive to site-specific properties, and adjustment factors are developed and presented. The proposed 𝑃⁢𝐺⁢𝐴𝐿 mapping procedure produces 𝐹⁢𝑆𝐿 estimates that are consistent with those obtained from full PLHA at a target 𝑇𝑅,𝐹⁢𝑆 , providing a promising roadmap to incorporating PLHA concepts into current liquefaction design methods.

Journal of Geotechnical and Geoenvironmental Engin

Water-quality trends of urban streams in Independence, Missouri, 2005–18

The U.S. Geological Survey and the city of Independence, Missouri, Water Pollution Control Department has studied the water quality and ecological condition of urban streams within Independence since 2005. Selected physical properties, nutrients, chloride, fecal indicator bacteria (Escherichia coli and total coliform), total dissolved solids, and suspended-sediment concentration data for base-flow and stormflow samples were used to document temporal trends in concentrations and flow-weighted concentrations; and annual loads were computed and investigated for selected nutrients, chloride, and suspended sediment. The six study sites included in this report are located on five urban streams: Rock Creek, a tributary in the city that drains to the Missouri River; three tributaries of the Little Blue River within the city (East Fork Little Blue River, Adair Creek, and Spring Branch Creek); and two sites on the main stem of the Little Blue River (one upstream from the city and one downstream from the three tributaries). Many factors such as population, land use, and climate, and combinations of these factors contributed to the significant changes in the concentrations and transport of nutrients, chloride, fecal indicator bacteria, and suspended sediment in the urban streams within Independence. The population of Independence and the amount of developed land in the urban watersheds remained unchanged during the 2005–18 study. Differences were noted in precipitation and in streamflow during the study. Annual precipitation and streamflow were separated into two time periods within the study—period 1 (2006–10), having greater annual streamflow and precipitation, and period 2 (2011–18), having about 30 percent lower annual streamflow and less precipitation. Streamflow was an important factor in the transport of nitrogen, phosphorus, chloride, and suspended sediment from the urban watersheds. Changes in data collection methodology during the study period and improvements to the city stormwater and wastewater infrastructure also could have contributed to some of the trends. Between 2009 and 2015, more than 35 million dollars of improvements were made to stormwater and wastewater infrastructure within the city. These improvements, such as additional sewage overflow holding tanks, removal of septic tanks, and improved and expanded sanitary sewer lines and storm overflows, also could have affected the decreased nutrients and fecal indicator bacteria trends among the urban streams in the study area. Models were used for analyzing streamflow-related variability in constituent concentrations and loads to determine if the water quality changed significantly during the study period. Trends in concentration data at four sites were analyzed using a statistical package called R–QWTREND and trends in load data were analyzed at six sites using a statistical package called Weighted Regressions on Time, Discharge, and Season-Kalman filter (WRTDS–K); both developed by the U.S. Geological Survey and publicly available for use. Statistically significant trends in flow-weighted nutrient concentrations and loads generally were downward during the study period. The only nutrient compound with a statistically significant upward trend in flow-weighted concentration was dissolved orthophosphate as phosphorus at the Rock Creek site and the upstream site on the Little Blue River. A statistically significant downward trend in annual dissolved ammonia load was identified at the downstream Little Blue River site. A significant upward linear trend in annual orthophosphate as phosphorus load was identified on Adair Creek. A statistically significant upward trend in dissolved chloride concentrations was identified at the downstream Little Blue River site. Road salt application near the site during the winter could have resulted in higher concentrated runoff during wet weather conditions. Annual chloride loads significantly decreased in Adair Creek and Spring Branch Creek. The mean annual chloride load transported in the drier (2011–18) period 2 was significantly less than during the wetter (2006–10) period 1, indicating that trends in precipitation runoff are an important factor in trends in annual transport of chloride. Statistically significant downward trends in flow-weighted fecal indicator bacteria Escherichia coli (E. coli) population densities were noted for Rock Creek and the down-stream site on the Little Blue River. However, no trend was identified in E. coli population density at the upstream Little Blue River site. The downward trend in E. coli population density at the downstream site could be a result of decreased streamflow and precipitation over the study period, storage of fecal indicator bacteria in the Little Blue River streambed within the study area, die-off of fecal indicator bacteria during travel from upstream to downstream, changes in the sample collection methodology, improvements to the city’s storm-water and wastewater infrastructures, or a combination of these factors. The statistically significant downward trend in suspended-sediment concentration identified at the upstream Little Blue River site could be affected by the decreased streamflow and precipitation during the study period, by changes in sampling methods within the study period, and by the decrease in construction and urban land development upstream from the city. No statistically significant change was indicated in the annual suspended-sediment load transported from Independence to the Little Blue River during the study period. More than one-half the suspended sediment transported in the Little Blue River originated in the watershed upstream from Independence. The Little Blue River and many of its tributaries that drain Independence have been designated as recreational waters classified for whole-body contact class B and secondary contact recreation, and some have been listed as impaired for E. coli by the Missouri Department of Natural Resources from urban runoff and storm sewers. Observations were made among the available E. coli population density data for both Little Blue River sites to further understand water-quality conditions over the study period. Both Little Blue River sites had similar medians and geometric means for the recreational season (April through October) and during the full study period, both of which are greater than the regulatory population density for both recreational classes. The Little Blue River drainage area nearly doubles in size from the upstream to downstream site; therefore, the consistent geometric mean and median of E. coli population densities at the upstream and downstream Little Blue River sites could be primarily due to the larger volume of streamflow creating a dilution effect. Other possible factors could be storage of fecal indicator bacteria in stream bed sediments, die-off of fecal indicator bacteria during transport, improvements to the city’s wastewater and stormwater infrastructure, changes to sampling methodology, or a combination of these factors. Specific sources of the E. coli are currently (2019) unknown.

Missouri

Long-term video surveillance and automated analyses reveal arousal patterns in groups of hibernating bats

Understanding natural behaviours is essential to determining how animals deal with new threats (e.g. emerging diseases). However, natural behaviours of animals with cryptic lifestyles, like hibernating bats, are often poorly characterized. White-nose syndrome (WNS) is an unprecedented disease threatening multiple species of hibernating bats, and pathogen-induced changes to host behaviour may contribute to mortality. To better understand the behaviours of hibernating bats and how they might relate to WNS, we developed new ways of studying hibernation across entire seasons. We used thermal-imaging video surveillance cameras to observe little brown bats ( Myotis lucifugus ) and Indiana bats ( M. sodalis ) in two caves over multiple winters. We developed new, sharable software to test for autocorrelation and periodicity of arousal signals in recorded video. We processed 740 days (17,760 hr) of video at a rate of >1,000 hr of video imagery in less than 1 hr using a desktop computer with sufficient resolution to detect increases in arousals during midwinter in both species and clear signals of daily arousal periodicity in infected M. sodalis . Our unexpected finding of periodic synchronous group arousals in hibernating bats demonstrate the potential for video methods and suggest some bats may have innate behavioural strategies for coping with WNS. Surveillance video and accessible analysis software make it now practical to investigate long-term behaviours of hibernating bats and other hard-to-study animals.

Methods in Ecology and Evolution

Quantification and simulation of metal loading to the Upper Animas River, Eureka to Silverton, San Juan County, Colorado, September 1997 and August 1998

Drainage from abandoned and inactive mines and from naturally mineralized areas in the San Juan Mountains of southern Colorado contributes metals to the upper Animas River near Silverton, Colorado. Tracer-injection studies and associated synoptic sampling were performed along two reaches of the upper Animas River to develop detailed profiles of stream discharge and to locate and quantify sources of metal loading. One tracer-injection study was performed in September 1997 on the Animas River reach from Howardsville to Silverton, and a second study was performed in August 1998 on the stream reach from Eureka to Howardsville. Drainage in the upper Animas River study reaches contributed aluminum, calcium, copper, iron, magnesium, manganese, sulfate, and zinc to the surface-water system in 1997 and 1998. Colloidal aluminum, dissolved copper, and dissolved zinc were attenuated through a braided stream reach downstream from Eureka. Instream dissolved copper concentrations were lower than the State of Colorado acute and chronic toxicity standards downstream from the braided reach to Silverton. Dissolved iron load and concentrations increased downstream from Howardsville and Arrastra Gulch, and colloidal iron remained constant at low concentrations downstream from Howardsville. Instream sulfate concentrations were lower than the U.S. Environmental Protection Agency's secondary drinking-water standard of 250 milligrams per liter throughout the two study reaches. Elevated zinc concentrations are the primary concern for aquatic life in the upper Animas River. In the 1998 Eureka to Howardsville study, instream dissolved zinc load increased downstream from the Forest Queen mine, the Kittimack tailings, and Howardsville. In the 1997 Howardsville to Silverton study, there were four primary areas where zinc load increased. First, was the increase downstream from Howardsville and abandoned mining sites downstream from the Cunningham Gulch confluence, which also was measured during the 1998 study. The second affected reach was downstream from Arrastra Gulch, where the increase in zinc load seems related to a series of right-bank inflows with low pH Quantification and Simulation of Metal Loading to the Upper Animas River, Eureka to Silverton, San Juan County, Colorado, September 1997 and August 1998By Suzanne S. Paschke, Briant A. Kimball, and Robert L. Runkeland elevated dissolved zinc concentrations. A third increase in zinc load occurred 6,100 meters downstream from the 1997 injection site and may have been from ground-water discharge with elevated zinc concentrations based on mass-loading graphs and the lack of visible inflow in the reach. A fourth but lesser dissolved zinc load increase occurred downstream from tailings near the Lackawanna Mill. Results of the tracer-injection studies and the effects of potential remediation were analyzed using the one- dimensional stream-transport computer code OTIS. Based on simulation results, instream zinc concentrations downstream from the Kittimack tailings to upstream from Arrastra Gulch would approach 0.16 milligram per liter (the upper limit of acute toxicity for some sensitive aquatic species) if zinc inflow concentrations were reduced by 75 percent in the stream reaches receiving inflow from the Forest Queen mine, the Kittimack tailings, and downstream from Howardsville. However, simulated zinc concentrations downstream from Arrastra Gulch were higher than approximately 0.30 milligram per liter due to numerous visible inflows and assumed ground-water discharge with elevated zinc concentrations in the lower part of the study reach. Remediation of discrete visible inflows seems a viable approach to reducing zinc inflow loads to the upper Animas River. Remediation downstream from Arrastra Gulch is more complicated because ground-water discharge with elevated zinc concentrations seems to contribute to the instream zinc load.

Colorado

Resolution analysis of finite fault source inversion using one- and three-dimensional Green's functions 2. Combining seismic and geodetic data

Using numerical tests for a prescribed heterogeneous earthquake slip distribution, we examine the importance of accurate Green's functions (GF) for finite fault source inversions which rely on coseismic GPS displacements and leveling line uplift alone and in combination with near-source strong ground motions. The static displacements, while sensitive to the three-dimensional (3-D) structure, are less so than seismic waveforms and thus are an important contribution, particularly when used in conjunction with waveform inversions. For numerical tests of an earthquake source and data distribution modeled after the 1994 Northridge earthquake, a joint geodetic and seismic inversion allows for reasonable recovery of the heterogeneous slip distribution on the fault. In contrast, inaccurate 3-D GFs or multiple 1-D GFs allow only partial recovery of the slip distribution given strong motion data alone. Likewise, using just the GPS and leveling line data requires significant smoothing for inversion stability, and hence, only a blurred vision of the prescribed slip is recovered. Although the half-space approximation for computing the surface static deformation field is no longer justifiable based on the high level of accuracy for current GPS data acquisition and the computed differences between 3-D and half-space surface displacements, a layered 1-D approximation to 3-D Earth structure provides adequate representation of the surface displacement field. However, even with the half-space approximation, geodetic data can provide additional slip resolution in the joint seismic and geodetic inversion provided a priori fault location and geometry are correct. Nevertheless, the sensitivity of the static displacements to the Earth structure begs caution for interpretation of surface displacements, particularly those recorded at monuments located in or near basin environments.

Journal of Geophysical Research B: Solid Earth

Estimation of Sediment Sources Using Selected Chemical Tracers in the Perry Lake and Lake Wabaunsee Basins, Northeast Kansas

In Kansas and nationally, stream and lake sediment is a primary concern as related to several important issues including water quality and reservoir water-storage capacity. The ability to achieve meaningful decreases in sediment loads to reservoirs requires a determination of the relative importance of sediment sources within the contributing basins. To investigate sources of sediment within the Perry Lake and Lake Wabaunsee Basins of northeast Kansas, representative samples of channel-bank sources, surface-soil sources (cropland and grassland), and reservoir bottom sediment were collected, analyzed, and compared. Subbasins sampled within the Perry Lake Basin included Atchison County Lake, Banner Creek Reservoir, Gregg Creek, Mission Lake, and Walnut Creek. The samples were sieved to isolate the less than 63-micron fraction (that is, the silt and clay) and analyzed for selected nutrients (total nitrogen and total phosphorus), organic and total carbon, 25 trace elements, and the radionuclide cesium-137 (137Cs). To determine which of the 30 constituents provided the best ability to discriminate between channel-bank and surface-soil sources in the two basins, four selection criteria were used. To be selected, it was required that the candidate constituent (1) was detectable, (2) had concentrations or activities that varied substantially and consistently between the sources, (3) had concentration or activity ranges that did not overlap between the sources, and (4) had concentration or activity differences between the sources that were statistically significant. On the basis of the four selection criteria, total nitrogen (TN), total phosphorus (TP), total organic carbon (TOC), and 137Cs were selected. Of the four selected constituents, 137Cs likely is the most reliable indicator of sediment source because it is known to be conservative in the environment. Trace elements were not selected because concentrations in the channel-bank and surface-soil sources generally were similar or did not vary in a consistent manner. To further account for differences in particle-size composition between the sources and the reservoir bottom sediment prior to the sediment-source estimations, constituent ratio and clay-normalization techniques were used. Computed ratios included the ratio of TOC to TN, TOC to TP, and TN to TP. Constituent concentrations (TN, TP, TOC) and activities (137Cs) were normalized by dividing by the percentage of clay. Thus, the sediment-source estimations involved the use of seven sediment-source indicators (that is, three constituent ratios and the clay-normalized concentration or activity for four constituents). Sediment-source estimation for each reservoir was based on a comparison between the reservoir bottom sediment and the end member channel-bank and surface-soil sources. Within the Perry Lake Basin, the seven-indicator consensus indicated that both channel-bank and surface-soil sources were important contributors of the sediment deposited in Atchison County Lake and Banner Creek Reservoir, whereas channel-bank sources were the dominant source of sediment for Mission Lake. On the sole basis of 137Cs activity, surface-soil sources contributed the most sediment to Atchison County Lake, and channel-bank sources contributed the most sediment to Banner Creek Reservoir and Mission Lake. For Perry Lake, both the seven-indicator consensus and 137Cs indicated that channel-bank sources were dominant and that channel-bank sources increased in importance with distance downstream in the Perry Lake Basin. For Lake Wabaunsee, the seven-indicator consensus and 137Cs indicated that both channel-bank and surface-soil sources were important. Given that the relative contribution of sediment from channel-bank and surface-soil sources can vary within and between basins and over time, basin-specific strategies for sediment management and monitoring are appropriate.

Scientific Investigations Report

Sources, instream transport, and trends of nitrogen, phosphorus, and sediment in the lower Tennessee River basin, 1980-96

In 1997, the U.S. Geological Survey (USGS) began an assessment of the lower Tennessee River Basin as part of the National Water-Quality Assessment Program. Existing nutrient and sediment data from 1980 to 1996 were compiled, screened, and interpreted to estimate watershed inputs from nutrient sources, provide a general description of the distribution and transport of nutrients and sediments in surface water, and evaluate trends in nutrient and sediment concentrations in the lower Tennessee (LTEN) River Basin. Nitrogen inputs from major sources varied widely among tributary basins in the LTEN River Basin. Point source wastewater discharges contributed between 0 and 0.61 tons per square mile per year [(tons/mi 2 )/yr]. Of the nonpoint sources of nitrogen for which inputs were estimated (atmospheric deposition, nitrogen fixation, fertilizer application, and livestock waste) livestock waste contributed the largest input in about two-thirds (7 out of 11) of the tributary basins, and fertilizer application contributed the largest input in the remaining 4 basins. Nitrogen input from fertilizer application was the most variable spatially among the nonpoint sources of nitrogen, ranging from 1.5 to 23 (tons/mi 2 )/yr. Atmospheric deposition estimates varied the least from basin to basin, ranging from 1.6 to 2.0 (tons/mi 2 )/yr. Estimates of nitrogen input from livestock waste ranged between 2.0 to 13 (tons/mi 2 )/yr. The percentage of the input from each of these nonpoint sources that entered the surface-water system is not known. Wastewater discharge contributed between 0 and 0.14 (ton/mi 2 )/yr of phosphorus to tributary basins. Livestock waste contributed most of the input in 8 out of the 11 basins, and fertilizer application contributed the most in the remaining 3 basins. Estimates of phosphorus input for fertilizer application ranged from 0.35 to 5.1 (tons/mi 2 )/yr and from 0.62 to 4.3 (tons/mi 2 )/yr from livestock waste. Reservoirs on the main stem of the Tennessee River and on the Duck and Elk Rivers affect nutrient transport because hydrodynamic conditions in the reservoirs promote assimilation by aquatic plants and deposition of particulate matter. Observed decreases in total nitrite plus nitrate and dissolved-orthophosphorus concentrations in reservoirs or at sites downstream of reservoirs during summer months were probably related to seasonality of plant growth. Nutrient and sediment data used to estimate annual instream loads and yields were compiled from various water-quality monitoring programs and represent the best available data in the LTEN River Basin, but these data have several characteristics that limit accuracy of load estimates. Many of the monitoring programs were not designed with the objective of annual load estimation, and data representing storm transport are, therefore, sparse; sampling and analytical methods varied through time and among the monitoring programs, hampering spatial and temporal comparisons. The load estimates computed from these data are useful for evaluating broad spatial patterns of instream load, and comparisons of instream load to inputs, but may not be sufficiently accurate for local-scale evaluations of water quality. Estimates of the mean annual instream load of total nitrogen entering (Chattanooga, Tenn.) and leaving (Paducah, Ky.) the LTEN River Basin were 29,000 and 60,000 tons per year (tons/yr), respectively. These estimates represent a gain of 31,000 tons/yr, on average, across the area (18,930 mi 2 ) between these inlet and outlet sites. The sum of the mean annual instream load from gaged tributaries to the main stem within the study unit was 14,000 tons/yr; however, this number cannot be directly compared with the gain between the inlet and outlet sites because (1) the gaged area represents only 30 percent of the total area and (2) the period of record at many tributary sites did not correspond with the period of record at the inlet or outlet sites. Estimates of mean annual instream load of total phosphorus at the inlet and outlet sites of the LTEN River Basin were 1,300 and 5,000 tons/yr, respectively, representing a gain of 3,700 tons/yr, on average, across the study unit. The sum of the gaged tributary load, representing only 28 percent of the area contributing to the main stem, was 4,300 tons/yr. Although this number cannot be closely compared with the gain throughout the study unit, for the same reasons given for total nitrogen, a general comparison suggests that the main stem of the Tennessee River and the tributary embayments along the main stem function as a sink for total phosphorus, removing a substantial amount from the water column through deposition or assimilation. The estimates of inputs can be compared and correlated with yields (area-normalized instream loads); significant correlations between estimates of inputs and yields might be useful as predictive tools for instream water quality where monitoring data are not available. Yields of nitrogen correlated moderately well with inputs from nonpoint sources, based on 1992 estimates. Nitrogen yield was highest [3.5 (tons/mi 2 )/yr] for Town Creek, for which the balance of nonpoint-source inputs to agricultural lands (fertilizer application plus nitrogen fixation plus livestock waste minus harvest) was also the highest [15 (tons/mi 2 )/yr]. Nitrogen yield was low [1.0 (tons/mi 2 )/yr] for the Buffalo River, for which the balance of agricultural nonpoint-source input was correspondingly low [3.2 (tons/mi 2 )/yr, the second lowest]. Correlation of wastewater discharge with yield was poor, and contrasted with the significant correlation between wastewater discharge and median nitrogen concentration during low streamflow. The poor correlation between wastewater discharge and annual yield was expected, however, as wastewater discharge is a small fraction compared with annual yield. In contrast with nitrogen, phosphorus yield did not correlate well with any estimated inputs or land-use types for the tributary basins. Phosphorus yield was highest [1.1 and 0.93 (tons/mi 2 )/yr] at two sites along the Duck River and at Elk River near Prospect [0.89 (ton/mi 2 )/yr]; however, estimates of inputs at these sites were in the middle of their respective ranges. The influence of the outcrop of phosphatic limestone formations of the brown-phosphate districts in the lower Duck and lower Elk River Basins might be responsible for the poor correlation between estimated inputs and yields of phosphorus. The outcrop pattern of these phosphatic limestones are an important factor to consider as regional boundaries are established for attainable, region-specific water-quality criteria for total phosphorus. Estimates of sediment input from cropland soil erosion in 1992 ranged from 51 to 540 (tons/mi 2 )/yr among the major hydrologic units in the LTEN River Basin. Information was not available to estimate this input for individual tributaries. Sediment yield estimates ranged from 65 to 263 (tons/mi 2 )/yr for the three tributary monitoring basins for which instream data were available, and from 17 to 26 (tons/mi 2 )/yr for the Tennessee River at South Pittsburg and at Pickwick Landing Dam, respectively. Lower sediment yields for the main stem sites compared with the tributary sites is probably due to sediment deposition in the main stem of the Tennessee River and tributary embayments along the main stem. Most of the significant trends in nutrient concentrations from about 1985 to about 1995 were decreasing trends, except for total nitrite plus nitrate, which increased at one site on the Elk River. The spatial distribution of decreasing trends of total nitrogen and total ammonia corresponds with the spatial variation among basins in wastewater loading rate. The time period of observed trends corresponds to the period of improvements in municipal treatment, thus decreases in wastewater effluent concentrations of nitrogen might be responsible for the decreasing trend in instream concentrations at these sites. Concentrations of total phosphorus did not decrease during this period at these sites, as might have been expected considering the reductions in wastewater input of phosphorus during this period.

Alabama, Georgia, Kentucky, Mississippi, Tennessee