Geology ReportsSearch

SEARCH · Geology Reports

Results for “Materials Performance”

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 433 records · Page 24Linked to original sources

Surface-water nutrient conditions and sources in the United States Pacific Northwest

The SPAtially Referenced Regressions On Watershed attributes (SPARROW) model was used to perform an assessment of surface-water nutrient conditions and to identify important nutrient sources in watersheds of the Pacific Northwest region of the United States (U.S.) for the year 2002. Our models included variables representing nutrient sources as well as landscape characteristics that affect nutrient delivery to streams. Annual nutrient yields were higher in watersheds on the wetter, west side of the Cascade Range compared to watersheds on the drier, east side. High nutrient enrichment (relative to the U.S. Environmental Protection Agency's recommended nutrient criteria) was estimated in watersheds throughout the region. Forest land was generally the largest source of total nitrogen stream load and geologic material was generally the largest source of total phosphorus stream load generated within the 12,039 modeled watersheds. These results reflected the prevalence of these two natural sources and the low input from other nutrient sources across the region. However, the combined input from agriculture, point sources, and developed land, rather than natural nutrient sources, was responsible for most of the nutrient load discharged from many of the largest watersheds. Our results provided an understanding of the regional patterns in surface-water nutrient conditions and should be useful to environmental managers in future water-quality planning efforts.

Journal of the American Water Resources Associatio

Pilot studies for the North American Soil Geochemical Landscapes Project - Site selection, sampling protocols, analytical methods, and quality control protocols

In 2004, the US Geological Survey (USGS) and the Geological Survey of Canada sampled and chemically analyzed soils along two transects across Canada and the USA in preparation for a planned soil geochemical survey of North America. This effort was a pilot study to test and refine sampling protocols, analytical methods, quality control protocols, and field logistics for the continental survey. A total of 220 sample sites were selected at approximately 40-km intervals along the two transects. The ideal sampling protocol at each site called for a sample from a depth of 0-5 cm and a composite of each of the O, A, and C horizons. The <2-mm fraction of each sample was analyzed for Al, Ca, Fe, K, Mg, Na, S, Ti, Ag, As, Ba, Be, Bi, Cd, Ce, Co, Cr, Cs, Cu, Ga, In, La, Li, Mn, Mo, Nb, Ni, P, Pb, Rb, Sb, Sc, Sn, Sr, Te, Th, Tl, U, V, W, Y, and Zn by inductively coupled plasma-mass spectrometry and inductively coupled plasma-atomic emission spectrometry following a near-total digestion in a mixture of HCl, HNO 3 , HClO 4 , and HF. Separate methods were used for Hg, Se, total C, and carbonate-C on this same size fraction. Only Ag, In, and Te had a large percentage of concentrations below the detection limit. Quality control (QC) of the analyses was monitored at three levels: the laboratory performing the analysis, the USGS QC officer, and the principal investigator for the study. This level of review resulted in an average of one QC sample for every 20 field samples, which proved to be minimally adequate for such a large-scale survey. Additional QC samples should be added to monitor within-batch quality to the extent that no more than 10 samples are analyzed between a QC sample. Only Cr (77%), Y (82%), and Sb (80%) fell outside the acceptable limits of accuracy (% recovery between 85 and 115%) because of likely residence in mineral phases resistant to the acid digestion. A separate sample of 0-5-cm material was collected at each site for determination of organic compounds. A subset of 73 of these samples was analyzed for a suite of 19 organochlorine pesticides by gas chromatography. Only three of these samples had detectable pesticide concentrations. A separate sample of A-horizon soil was collected for microbial characterization by phospholipid fatty acid analysis (PLFA), soil enzyme assays, and determination of selected human and agricultural pathogens. Collection, preservation and analysis of samples for both organic compounds and microbial characterization add a great degree of complication to the sampling and preservation protocols and a significant increase to the cost for a continental-scale survey. Both these issues must be considered carefully prior to adopting these parameters as part of the soil geochemical survey of North America.

Applied Geochemistry

Exploring the potential of electrospray-Orbitrap for stable isotope analysis using nitrate as a model

Widely used isotope ratio mass spectrometers have limited capabilities to measure metabolites, drugs, or small polyatomic ions without the loss of structural isotopic information. A new approach has recently been introduced that uses electrospray ionization Orbitrap to measure multidimensional isotope signatures of intact polar compounds. Using nitrate as a model compound, this study aims to establish performance metrics for comparisons with conventional IRMS at the natural abundance level. We present a framework on how to convert isotopolog intensities to δ values that are commonly used in the isotope geochemistry community. The quantification of seven nitrate isotopologs provides multiple pathways for obtaining the primary N and O δ values including non-mass-dependent O isotope variations, as well as opportunities to explore nonrandom isotopic distributions (i.e., clumping effects) within molecular nitrate. Using automation and the adaptation of measurement principles that are specific to isotope ratio analysis, nitrate δ15NAIR, δ18OVSMOW, and δ17OVSMOW were measured with a long-term precision of 0.4‰ or better for isotopic reference materials and purified nitrate from environmental samples. In addition, we demonstrate promising results for unpurified environmental samples in liquid form. With these new developments, this study connects the two largely disparate mass spectrometry fields of bioanalytical MS and isotope ratio MS, thus providing a route to measure new isotopic signatures in diverse organic and inorganic solutes.

Analytical Chemistry

Implementation of the Next Generation Attenuation (NGA) ground-motion prediction equations in Fortran and R

This report presents two methods for implementing the earthquake ground-motion prediction equations released in 2008 as part of the Next Generation Attenuation of Ground Motions (NGA-West, or NGA) project coordinated by the Pacific Earthquake Engineering Research Center (PEER). These models were developed for predicting ground-motion parameters for shallow crustal earthquakes in active tectonic regions (such as California). Of the five ground-motion prediction equations (GMPEs) developed during the NGA project, four models are implemented: the GMPEs of Abrahamson and Silva (2008), Boore and Atkinson (2008), Campbell and Bozorgnia (2008), and Chiou and Youngs (2008a); these models are abbreviated as AS08, BA08, CB08, and CY08, respectively. Since site response is widely recognized as an important influence of ground motions, engineering applications typically require that such effects be modeled. The model of Idriss (2008) is not implemented in our programs because it does not explicitly include site response, whereas the other four models include site response and use the same variable to describe the site condition (VS30). We do not intend to discourage the use of the Idriss (2008) model, but we have chosen to implement the other four NGA models in our programs for those users who require ground-motion estimates for various site conditions. We have implemented the NGA models by using two separate programming languages: Fortran and R (R Development Core Team, 2010). Fortran, a compiled programming language, has been used in the scientific community for decades. R is an object-oriented language and environment for statistical computing that is gaining popularity in the statistical and scientific community. Derived from the S language and environment developed at Bell Laboratories, R is an open-source language that is freely available at http://www.r-project.org/ (last accessed 11 January 2011). In R, the functions for computing the NGA equations can be loaded as an add-on user-contributed code, which is referred to as a ?package? in R. The details of the nga package (Kaklamanos and Thompson, 2010) are presented in this report. In addition, differences between the R and Fortran implementations will be discussed later in this report. The NGA models have established a new baseline for seismic hazard assessments, and they have been incorporated into the most recent national seismic hazard maps published by the U.S. Geological Survey (Petersen and others, 2008). However, many of the new models are significantly more complicated than previous GMPEs and, therefore, require a substantial investment of time to implement and validate. We hope that the scientific and engineering communities find our implementations to be useful in research and practice. Our implementations may be considered as an alternate to the Microsoft Excel spreadsheet implementation available on the PEER NGA project Web site (http://peer.berkeley.edu/ngawest/, last accessed 11 January 2011). The implementations in Fortran and R are more appropriate for performing batch calculations than the implementation in Microsoft Excel. Spreadsheets and Fortran code for some of the individual models also are available on the PEER NGA project Web site; our programs implement the four GMPEs simultaneously. Our programs give the same results as the programs on the PEER NGA Web site, but we offer some additional flexibility of input, additional methods of estimating unknown input parameters, and additional options for output. Although these programs have been used by the U.S. Geological Survey (USGS), Tufts University, and others, no warranty, expressed or implied, is made by Tufts or the USGS as to the accuracy or functioning of the programs and related material, nor shall the fact of distribution constitute any such warranty, and no responsibility is assumed by Tufts or the USGS in connection therewith.

Open-File Report

Shear surface undulations modulate clayey gouge strength and contribute to divergent landslide acceleration

Landslides display a spectrum of speeds for incompletely known reasons. Sliding occurs along slickensided undulatory shear surfaces within boundary shear gouge. Laboratory tests reveal that gouge shear strength generally decreases with finite cumulative displacement during relatively rapid failure and may increase or decrease with increasing shear rate; these behaviors can result in accelerating or decelerating landslide motion, which significantly affects consequent hazards. However, mechanisms responsible for such behaviors are poorly understood. We performed advanced ring shear strength testing that revealed such variable strength of a landslide near Oso, Washington, USA. We hypothesized that millimeter-scale undulations along shear surfaces caused the strength variability by imparting shear strength but while also modifying stresses that locally increase and decrease the typically considered particle-scale shear strength. We tested our hypotheses in the laboratory and with finite element soil deformation modeling. Lab results suggest that undulations contribute strength that decays with finite cumulative displacement. Modeling similarly reveals this, and that rapid shearing across undulations locally reduces effective normal stress by persistently elevating pore-water pressure and causing dilation. Consequent effects on strength differ by material with granular-rich, high-friction gouge losing strength and clay-rich, low-friction gouge losing little or gaining strength as shear rates increase. Ample testing by others reveals similar patterns. Hence, the propensity for gouge-controlled accelerating or decelerating failure may be estimated from simple index tests. Our findings on the effects of undulations reveal previously unknown mechanisms that may help to explain why some landslides reactivate catastrophically while others do not.

Washington

Groundwater flow model for the Little Plover River basin in Wisconsin’s Central Sands

The Little Plover River is a groundwater-fed stream in the sand plains region of central Wisconsin. In this region, sandy sediment deposited during or soon after the last glaciation forms an important unconfined sand and gravel aquifer. This aquifer supplies water for numerous high-capacity irrigation, municipal, and industrial wells that support a thriving agricultural industry. In recent years the addition of many new wells, combined with observed diminished flows in the Little Plover and other nearby rivers, has raised concerns about the impacts of the wells on groundwater levels and on water levels and flows in nearby lakes, streams, and wetlands. Diverse stakeholder groups, including well operators, Growers, environmentalists, local land owners, and regulatory and government officials have sought a better understanding of the local groundwater-surface water system and have a shared desire to balance the water needs of the he liagricultural, industrial, and urban users with the maintenance and protection of groundwater-dependent natural resources. To help address these issues, the Wisconsin Department of Natural Resources requested that the Wisconsin Geological and Natural History Survey and U.S. Geological Survey cooperatively develop a groundwater flow model that could be used to demonstrate the relationships among groundwater, surface water, and well withdrawals and also be a tool for testing and evaluating alternative water management strategies for the central sands region. Because of an abundance of previous studies, data availability, local interest, and existing regulatory constraints the model focuses on the Little Plover River watershed, but the modeling methodology developed during this study can apply to much of the larger central sands of Wisconsin. The Little Plover River groundwater flow model simulates three-dimensional groundwater movement in and around the Little Plover River basin under steady-state and transient conditions. This model explicitly includes all high-capacity wells in the model domain and simulates seasonal variations in recharge and well pumping. The model represents the Little Plover River, and other significant streams and drainage ditches in the model domain, as fully connected to the groundwater system, computes stream base flow resulting from groundwater discharge, and routes the flow along the stream channel. A separate soil-water-balance (SWB) model was used to develop groundwater recharge arrays as input for the groundwater flow model. The SWB model uses topography, soils, land use, and climatic data to estimate recharge as deep drainage from the soil zone. The SWB model explicitly includes recharge originating as irrigation water, and computes irrigation using techniques similar to those used by local irrigation operators. The groundwater flow model uses the U.S. Geological Survey’s MODFLOW modeling code which is freely available, widely accepted, and commonly used by the groundwater community. The groundwater flow model and the SWB model use identical high-resolution numerical grids having model cells 100 feet on a side, with physical properties assigned to each grid cell. This grid allows accurate geographic placement of wells, streams, and other model features. The 3-dimensional grid has three layers; layers 1 and 2 represent the sand and gravel aquifer and layer 3 represents the underlying sandstone. The distribution of material properties in the model (hydraulic conductivity, aquifer thickness, etc.) comes from previous published geologic studies of the region, updated by calibration to recent streamflow and groundwater level data. The SWB model operates on a daily time step. The groundwater flow model was calibrated to monthly stress periods with time steps ranging from 1 to 16 days. More detailed time discretization is possible. The groundwater model was calibrated to water-level and streamflow data collected during 2013 and 2014 by adjusting model parameters (primarily hydraulic conductivity, storage, and recharge) until the model produced a conditionally optimal fit between field observations and model output, subject to consistency with previously published geologic studies. Calibration was performed under both steady and transient conditions, and used a sophisticated parameter-estimation procedure (PEST) for the calibration process and to identify important model parameters. For the Little Plover River, the two most important parameters are the global recharge multiplier and the hydraulic conductivity of the stream bed. The calibrated model produces water-level and mass-balance results that are consistent with field observations and previous studies of the area. The completed model is a powerful tool for testing and demonstrating alternative water-management scenarios. Example model applications described in this report include simulating how the cumulative impacts of pumping and land-use change have affected average baseflow in the Little Plover River. Depletion-potential mapping represents a method for predicting which wells and well locations have the greatest impact on nearby surface-water resources. The completed model is publicly available, along with a companion user’s guide to assist with its operation, at http://wgnhs.org/littleplover- river-groundwater-model.

Wisconsin

Utility of an instantaneous salt dilution method for measuring streamflow in headwater streams

Streamflow records are biased toward large streams and rivers, yet small headwater streams are often the focus of ecological research in response to climate change. Conventional flow measurement instruments such as acoustic Doppler velocimeters (ADVs) do not perform well during low-flow conditions in small streams, truncating the development of rating curves during critical baseflow conditions dominated by groundwater inflow. We revisited an instantaneous solute tracer injection method as an alternative to ADVs based on paired measurements to compare their precision, efficiency, and feasibility within headwater streams across a range of flow conditions. We show that the precision of discharge measurements using salt dilution by slug injection and ADV methods were comparable overall, but salt dilution was more precise during the lowest flows and required less time to implement. Often, headwater streams were at or below the depth threshold where ADV measurements could even be attempted and transects were complicated by coarse bed material and cobbles. We discuss the methodological benefits and limitations of salt dilution by slug injection and conclude that the method could facilitate a proliferation of streamflow observation across headwater stream networks that are highly undersampled compared to larger streams.

Virginia

Suspended particulate layers and internal waves over the southern Monterey Bay continental shelf: an important control on shelf mud belts?

Physical and optical measurements taken over the mud belt on the southern continental shelf of Monterey Bay, California documented the frequent occurrence of suspended particulate matter features, the majority of which were detached from the seafloor, centered 9–33 m above the bed. In fall 2011, an automated profiling mooring and fixed instrumentation, including a thermistor chain and upward-looking acoustic Doppler current profiler, were deployed at 70 m depth for 5 weeks, and from 12 to 16 October a long-range autonomous underwater vehicle performed across-shelf transects. Individual SPM events were uncorrelated with local bed shear stress caused by surface waves and bottom currents. Nearly half of all observed SPM layers occurred during 1 week of the study, 9–16 October 2011, and were advected past the fixed profiling mooring by the onshore phase of semidiurnal internal tide bottom currents. At the start of the 9–16 October period, we observed intense near-bed vertical velocities capable of lifting particulates into the middle of the water column. This “updraft” event appears to have been associated with nonlinear adjustment of high-amplitude internal tides over the mid and outer shelf. These findings suggest that nonlinear internal tidal motions can erode material over the outer shelf and that, once suspended, this SPM can then be transported shoreward to the middle and shallow sections of the mud belt. This represents a fundamental broadening of our understanding of how shelf mud belts may be built up and sustained.

California

Experimental deformation of polycrystalline H2O ice at high pressure and low temperature: Preliminary results

Interest in the mechanical properties of water ice under the conditions in which it exists in the outer solar system has motivated the development and use of a new high‐pressure, low‐temperature triaxial deformation apparatus. Constant displacement rate tests on 70 samples of pure polycrystalline water ice have been performed at temperatures 77≤≤258 K, confining pressures 0.1≤≤350 MPa, and strain rates 3.5×10≤ ≤3.5×10 s. In most cases, the ice polymorph tested was ice . Both brittle and ductile behavior have been observed. Brittle behavior of ice, promoted by lower pressure, lower temperature, and higher strain rate, is analogous to that in rocks, with the important exception that brittle fracture strength becomes independent of confining pressure above 50 MPa pressure and the fracture angle is approximately 45° to the loading direction (i.e., the coefficient of internal friction is approximately zero). Ductile flow, the predominant behavior in our tests at ≥195 K, follows a law of form = σ exp (−*/) (σ is stress; is the gas constant; , * are material constants). Three sets of material constants are required to fit the data, with changes in sets (or mechanisms) occurring near 243 K and 195 K. The value of remains near 4 throughout the measured ductile field, but * drops from 91 to 61 to 31 /mole as temperature decreases. The maximum brittle strength measured was 171 MPa; the maximum ductile strength measured was 91 MPa. At confining pressures near the phase transition pressure of ice → ice II, the ductile strength is observed to drop dramatically. Some overlap with previous work occurs at higher temperatures and lower pressures. Agreement with present work is generally good, both quantitatively in the values of and *, and qualitatively in the mechanism of deformation. Although the ductile strengths measured here are somewhat higher than expected on the basis of extrapolations of previous work, the low value of * at <195 K indicates that the ice layer on icy bodies in the solar system is much weaker than has generally been predicted.

Journal of Geophysical Research B: Solid Earth

Effective groundwater model calibration: With analysis of data, sensitivities, predictions, and uncertainty

Methods and guidelines for developing and using mathematical models Turn to Effective Groundwater Model Calibration for a set of methods and guidelines that can help produce more accurate and transparent mathematical models. The models can represent groundwater flow and transport and other natural and engineered systems. Use this book and its extensive exercises to learn methods to fully exploit the data on hand, maximize the model's potential, and troubleshoot any problems that arise. Use the methods to perform: Sensitivity analysis to evaluate the information content of data Data assessment to identify (a) existing measurements that dominate model development and predictions and (b) potential measurements likely to improve the reliability of predictions Calibration to develop models that are consistent with the data in an optimal manner Uncertainty evaluation to quantify and communicate errors in simulated results that are often used to make important societal decisions Most of the methods are based on linear and nonlinear regression theory. Fourteen guidelines show the reader how to use the methods advantageously in practical situations. Exercises focus on a groundwater flow system and management problem, enabling readers to apply all the methods presented in the text. The exercises can be completed using the material provided in the book, or as hands-on computer exercises using instructions and files available on the text's accompanying Web site. Throughout the book, the authors stress the need for valid statistical concepts and easily understood presentation methods required to achieve well-tested, transparent models. Most of the examples and all of the exercises focus on simulating groundwater systems; other examples come from surface-water hydrology and geophysics. The methods and guidelines in the text are broadly applicable and can be used by students, researchers, and engineers to simulate many kinds systems.

Book

The sedimentology and dynamics of crater-affiliated wind streaks in western Arabia Terra, Mars and Patagonia, Argentina

Wind streaks comprise recent aeolian deposits that have been extensively documented on Venus, Earth and Mars. Martian wind streaks are among the most abundant surface features on the planet and commonly extend from the downwind margins of impact craters. Previous studies of wind streaks emerging from crater interior deposits suggested that the mode of emplacement was primarily related to the deposition of silt-sized particles as these settled from plumes. We have performed geologic investigations of two wind streaks clusters; one situated in western Arabia Terra, a region in the northern hemisphere of Mars, and another in an analogous terrestrial site located in southern Patagonia, Argentina, where occurrences of wind streaks emanate from playas within maar craters. In both these regions we have identified bedforms in sedimentary deposits on crater floors, along wind-facing interior crater margins, and along wind streaks. These observations indicate that these deposits contain sand-sized particles and that sediment migration has occurred via saltation from crater interior deposits to wind streaks. In Arabia Terra and in Patagonia wind streaks initiate from crater floors that contain lithic and evaporitic sedimentary deposits, suggesting that the composition of wind streak source materials has played an important role in development. Spatial and topographic analyses suggest that regional clustering of wind streaks in the studied regions directly correlates to the areal density of craters with interior deposits, the degree of proximity of these deposits, and the craters' rim-to-floor depths. In addition, some (but not all) wind streaks within the studied clusters have propagated at comparable yearly (Earth years) rates. Extensive saltation is inferred to have been involved in its propagation based on the studied terrestrial wind streak that shows ripples and dunes on its surface and the Martian counterpart changes orientation toward the downslope direction where it extends into an impact crater. ?? 2009 Elsevier B.V.

Geomorphology

A terrain-based site characterization map of California with implications for the contiguous United States

We present an approach based on geomorphometry to predict material properties and characterize site conditions using the V S 30 parameter (time‐averaged shear‐wave velocity to a depth of 30 m). Our framework consists of an automated terrain classification scheme based on taxonomic criteria (slope gradient, local convexity, and surface texture) that systematically identifies 16 terrain types from 1‐km spatial resolution (30 arcsec) Shuttle Radar Topography Mission digital elevation models (SRTMDEMs). Using 853 V S 30 values from California, we apply a simulation‐based statistical method to determine the mean V S 30 for each terrain type in California. We then compare the V S 30 values with models based on individual proxies, such as mapped surface geology and topographic slope, and show that our systematic terrain‐based approach consistently performs better than semiempirical estimates based on individual proxies. To further evaluate our model, we apply our California‐based estimates to terrains of the contiguous United States. Comparisons of our estimates with 325 V S 30 measurements outside of California, as well as estimates based on the topographic slope model, indicate our method to be statistically robust and more accurate. Our approach thus provides an objective and robust method for extending estimates of V S 30 for regions where in situ measurements are sparse or not readily available.

California

Rotenone persistence in freshwater ponds: Effects of temperature and sediment adsorption

The persistence of rotenone was compared between a cement-lined pond (0.04 hectare) and an earthen-bottom pond (0.02 hectare) treated with 5 μL Noxfish/L (250 μg rotenone/L) during spring, summer, and fall. Water temperatures on the days of treatment in each season were 8, 22, and 15°C, respectively. Both ponds were filled with pond water from a common source 1 week before each of the three treatments. Water samples (filtered and unfiltered) and sediment samples were analyzed by high-performance liquid chromatography to monitor the decrease of rotenone until residues were at or below the detection limit (<2.0 μg/L for water and < 25 ng/g for sediments). The loss of rotenone from water generally followed a first-order rate ofdecay. Rotenone disappeared two to three times faster in the earthen pond than in the concrete pond. The rotenone half-life times in the spring, summer, and fall treatments were 3.7, 1.3, and 5.2 d, respectively, in the concrete pond, and 1.8, 0.7, and 1.8 d in the earthen pond. Rates of decay in both ponds were directly correlated with water temperature. Filtered water samples from both ponds contained less rotenone than unfiltered water, indicating that some rotenone was bound to suspended material. The highest concentration of rotenone in sediment samples was 102 ng/g; residues decreased to below the detection limit within 14 d in the spring treatment and within 3 d in the summer and fall treatments.

North American Journal of Fisheries Management

Modelling root reinforcement in shallow forest soils

A hypothesis used to explain the relationship between timber harvesting and landslides is that tree roots add mechanical support to soil, thus increasing soil strength. Upon harvest, the tree roots decay which reduces soil strength and increases the risk of management -induced landslides. The technical literature does not adequately support this hypothesis. Soil strength values attributed to root reinforcement that are in the technical literature are such that forested sites can't fail and all high risk, harvested sites must fail. Both unstable forested sites and stable harvested sites exist, in abundance, in the real world thus, the literature does not adequately describe the real world. An analytical model was developed to calculate soil strength increase due to root reinforcement. Conceptually, the model is composed of a reinforcing element with high tensile strength, i.e. a conifer root, embedded in a material with little tensile strength, i.e. a soil. As the soil fails and deforms, the reinforcing element also deforms and stretches. The lateral deformation of the reinforcing element is treated analytically as a laterally loaded pile in a flexible foundation and the axial deformation is treated as an axially loaded pile. The governing differential equations are solved using finite-difference approximation techniques. The root reinforcement model was tested by comparing the final shape of steel and aluminum rods, parachute cord, wooden dowels, and pine roots in direct shear with predicted shapes from the output of the root reinforcement model. The comparisons were generally satisfactory, were best for parachute cord and wooden dowels, and were poorest for steel and aluminum rods. A parameter study was performed on the root reinforcement model which showed reinforced soil strength increased with increasing root diameter and soil depth. Output from the root reinforcement model showed a strain incompatibility between large and small diameter roots. The peak increase in soil strength attributed to roots was controlled by the small (<4mm) diameter root fraction. These results were used to calculate the effect of timber harvesting on a small, approximately 7.6 m 3 (10 yd 3 ), hypothetical landslide in a shallow, cohesionless, forest soil. The root reinforcement model predicted a post-harvest reduction in soil strength of 14 and 19 percent for a soil with and without 5 kPa (105 lbs/ft 2 ) of cohesion, respectively.

Thesis

Vertical profiles of streambed hydraulic conductivity determined using slug tests in central and western Nebraska

Many issues of water-resources management rely on modeling of ground-water/surfacewater interactions, and streambed hydraulic conductivity is a key parameter controlling the water fluxes across the stream/aquifer interface. However, in central and western Nebraska, this parameter is generally undefined. The U.S. Geological Survey, in cooperation with the Nebraska Platte River Cooperative Hydrology Study Group, performed slug tests at 15 stream sites in the Platte, Republican, and Little Blue River watersheds to determine the hydraulic conductivity of streambeds in central and western Nebraska. Slug tests were completed at several discrete depth intervals using pneumatic or mechanical methods, and the water-level response was monitored on site using a pressure transducer and laptop computer. Responses were analyzed using either the Bouwer and Rice or Springer and Gelhar methods. Vertical profiles of hydraulic conductivity with depth were developed and were compared to available information on lithology. The profiles and corresponding lithology showed that different types of streambeds were tested and suggested that some streambeds display a large variability in hydraulic conductivity with depth. In some cases, hydraulic conductivity values associated with nonstreambed materials could be identified from nearby lithologic descriptions. Seven of 15 sites had streambed values that ranged over more than 3 orders of magnitude, and that variability increased significantly when the measurements considered to be from nonstreambed materials were included. Streambed profiles from the Platte and South Platte River sites generally were more homogeneous and of larger hydraulic conductivity than the other sites. No restrictive layers were detected at any of the streambed sites on the main stems or the flood plains of the main stems of their respective watersheds. Alternatively, the profiles characterized by a restrictive streambed layer at some depth below the streambed surface were all from tributary sites out of the main-stem flood plain. These profiles can be used to represent the streambed hydraulic conductivity in central and western Nebraska in various applications, including modeling ground-water/surfacewater interactions.

Nebraska

Evaluation of multiple laboratory performance and variability in analysis of recreational freshwaters by a rapid Escherichia coli qPCR method (Draft Method C)

There is interest in the application of rapid quantitative polymerase chain reaction (qPCR) methods for recreational freshwater quality monitoring of the fecal indicator bacteria Escherichia coli ( E. coli ). In this study we determined the performance of 21 laboratories in meeting proposed, standardized data quality acceptance (QA) criteria and the variability of target gene copy estimates from these laboratories in analyses of 18 shared surface water samples by a draft qPCR method developed by the U.S. Environmental Protection Agency (EPA) for E. coli . The participating laboratories ranged from academic and government laboratories with more extensive qPCR experience to “new” water quality and public health laboratories with relatively little previous experience in most cases. Failures to meet QA criteria for the method were observed in 24% of the total 376 test sample analyses. Of these failures, 39% came from two of the “new” laboratories. Likely factors contributing to QA failures included deviations in recommended procedures for the storage and preparation of reference and control materials. A master standard curve calibration model was also found to give lower overall variability in log 10 target gene copy estimates than the delta-delta Ct (ΔΔCt) calibration model used in previous EPA qPCR methods. However, differences between the mean estimates from the two models were not significant and variability between laboratories was the greatest contributor to overall method variability in either case. Study findings demonstrate the technical feasibility of multiple laboratories implementing this or other qPCR water quality monitoring methods with similar data quality acceptance criteria but suggest that additional practice and/or assistance may be valuable, even for some more generally experienced qPCR laboratories. Special attention should be placed on providing and following explicit guidance on the preparation, storage and handling of reference and control materials.

Water Research

Determination of steroid hormones and related compounds in filtered and unfiltered water by solid-phase extraction, derivatization, and gas chromatography with tandem mass spectrometry

A new analytical method has been developed and implemented at the U.S. Geological Survey National Water Quality Laboratory that determines a suite of 20 steroid hormones and related compounds in filtered water (using laboratory schedule 2434) and in unfiltered water (using laboratory schedule 4434). This report documents the procedures and initial performance data for the method and provides guidance on application of the method and considerations of data quality in relation to data interpretation. The analytical method determines 6 natural and 3 synthetic estrogen compounds, 6 natural androgens, 1 natural and 1 synthetic progestin compound, and 2 sterols: cholesterol and 3--coprostanol. These two sterols have limited biological activity but typically are abundant in wastewater effluents and serve as useful tracers. Bisphenol A, an industrial chemical used primarily to produce polycarbonate plastic and epoxy resins and that has been shown to have estrogenic activity, also is determined by the method. A technique referred to as isotope-dilution quantification is used to improve quantitative accuracy by accounting for sample-specific procedural losses in the determined analyte concentration. Briefly, deuterium- or carbon-13-labeled isotope-dilution standards (IDSs), all of which are direct or chemically similar isotopic analogs of the method analytes, are added to all environmental and quality-control and quality-assurance samples before extraction. Method analytes and IDS compounds are isolated from filtered or unfiltered water by solid-phase extraction onto an octadecylsilyl disk, overlain with a graded glass-fiber filter to facilitate extraction of unfiltered sample matrices. The disks are eluted with methanol, and the extract is evaporated to dryness, reconstituted in solvent, passed through a Florisil solid-phase extraction column to remove polar organic interferences, and again evaporated to dryness in a reaction vial. The method compounds are reacted with activated -methyl--trimethylsilyl trifluoroacetamide at 65 degrees Celsius for 1 hour to form trimethylsilyl or trimethylsilyl-enol ether derivatives that are more amenable to gas chromatographic separation than the underivatized compounds. Analysis is carried out by gas chromatography with tandem mass spectrometry using calibration standards that are derivatized concurrently with the sample extracts. Analyte concentrations are quantified relative to specific IDS compounds in the sample, which directly compensate for procedural losses (incomplete recovery) in the determined and reported analyte concentrations. Thus, reported analyte concentrations (or analyte recoveries for spiked samples) are corrected based on recovery of the corresponding IDS compound during the quantification process. Recovery for each IDS compound is reported for each sample and represents an absolute recovery in a manner comparable to surrogate recoveries for other organic methods used by the National Water Quality Laboratory. Thus, IDS recoveries provide a useful tool for evaluating sample-specific analytical performance from an absolute mass recovery standpoint. IDS absolute recovery will differ and typically be lower than the corresponding analyte’s method recovery in spiked samples. However, additional correction of reported analyte concentrations is unnecessary and inappropriate because the analyte concentration (or recovery) already is compensated for by the isotope-dilution quantification procedure. Method analytes were spiked at 10 and 100 nanograms per liter (ng/L) for most analytes (10 times greater spike levels were used for bisphenol A and 100 times greater spike levels were used for 3--coprostanol and cholesterol) into the following validation-sample matrices: reagent water, wastewater-affected surface water, a secondary-treated wastewater effluent, and a primary (no biological treatment) wastewater effluent. Overall method recovery for all analytes in these matrices averaged 100 percent, with overall relative standard deviation of 28 percent. Mean recoveries of the 20 individual analytes for spiked reagent-water samples prepared along with field samples and analyzed in 2009–2010 ranged from 84–104 percent, with relative standard deviations of 6–36 percent. Concentrations for two analytes, equilin and progesterone, are reported as estimated because these analytes had excessive bias or variability, or both. Additional database coding is applied to other reported analyte data as needed, based on sample-specific IDS recovery performance. Detection levels were derived statistically by fortifying reagent water at six different levels (0.1 to 4 ng/L) and range from about 0.4 to 4 ng/L for 16 analytes. Interim reporting levels applied to analytes in this report range from 0.8 to 8 ng/L. Bisphenol A and the sterols (cholesterol and 3-beta-coprostanol) were consistently detected in laboratory and field blanks. The minimum reporting levels were set at 100 ng/L for bisphenol A and at 200 ng/L for the two sterols to prevent any bias associated with the presence of these compounds in the blanks. A minimum reporting level of 2 ng/L was set for 11-ketotestosterone to minimize false positive risk from an interfering siloxane compound emanating as chromatographic-column bleed, from vial septum material, or from other sources at no more than 1 ng/L.

Techniques and Methods

Perturbational and nonperturbational inversion of Rayleigh-wave velocities

The inversion of Rayleigh-wave dispersion curves is a classic geophysical inverse problem. We have developed a set of MATLAB codes that performs forward modeling and inversion of Rayleigh-wave phase or group velocity measurements. We describe two different methods of inversion: a perturbational method based on finite elements and a nonperturbational method based on the recently developed Dix-type relation for Rayleigh waves. In practice, the nonperturbational method can be used to provide a good starting model that can be iteratively improved with the perturbational method. Although the perturbational method is well-known, we solve the forward problem using an eigenvalue/eigenvector solver instead of the conventional approach of root finding. Features of the codes include the ability to handle any mix of phase or group velocity measurements, combinations of modes of any order, the presence of a surface water layer, computation of partial derivatives due to changes in material properties and layer boundaries, and the implementation of an automatic grid of layers that is optimally suited for the depth sensitivity of Rayleigh waves.

Geophysics