USGS ScienceSearch

SEARCH · USGS Science

Results for “Nature Methods”

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 415 records · Page 23Linked to original sources

Methods of analysis by the U.S. Geological Survey National Water Quality Laboratory-A method supplement for the determination of Fipronil and degradates in water by gas chromatography/mass spectrometry

A method for the isolation and detemrination of fipronil and four of its degradates has been developed. This method adapts an analytical method created by the U.S. Geological Survey National Water Quality Laboratory in 1995 for the determination of a broad range of high-use pesticides typically found in filtered natural-water samples. In 2000, fipronil and four of its degradates were extracted, analyzed, and validated using this method. The recoveries for these five compounds in reagent-water samples fortified at 1 microgram per liter (ug/L) avereraged 98 percent. Initial method detection limits averaged 0.0029 ug/L. The performance of these five new compounds is consistent with the performance of the compounds in the initial method, making it possible to include them in addition to the other 41 pesticides and pesticide degradates in the original method.

Open-File Report

Hydrologic budget of the Beaverdam Creek basin, Maryland

A hydrologic budget is a statement accounting for the water gains and losses for selected periods in an area. Weekly measurements of precipitation streamflow, surface-water storage, ground-water stage, and soil resistivity were made during a 2year period, April 1, 1950, to March 28, 1952, in the Beaverdam Creek basin, Wicomico County, Md. The hydrologic measurements are summarized in two budgets, a total budget and a ground-water budget, and in supporting tables and graphs. The results of the investigation have some potentially significant applications because they describe a method for determining the annual replenishment of the water supply of a basin and the ways of water disposal under natural conditions. The information helps to determine the 'safe' yield of water in diversion from natural to artificial discharge. The drainage basin of Beaverdam Creek was selected because it appeared to have fewer hydrologic variables than are generally found. However, the methods may prove applicable in many places under a variety of conditions. The measurements are expressed in inches of water over the area of the basin. The equation of the hydrologic cycle is the budget balance: P= R+E+ASW+ delta SW + delta SM + delta GW where P is precipitation; R is runoff; ET is evapotranspiration; delta SW is change in surface-water storage; delta SM is change in soil moisture; and delta GW is change in ground-water storage. In this report 'change' is the final quantity minus the initial quantity and thus is synonymous with 'increase.' Further, ,delta GW= delta H .x Yg, in which delta H is the change in ground-water stage and Yg is the gravity yield, or the specific yield of the sediments as measured during the short periods of declining ground-water levels characteristic of the area. The complex sum of the revised equation P ? R - delta SW ? ET - delta SM, which is equal to delta H. x Yg, has been named the 'infiltration residual'; it is equivalent to ground-water recharge. Two unmeasured, but not entirely unknown, quantities, evapotranspiration, (ET) and gravity yield, (Yg), are included in the equation. They are derived statistically by a method of convergent approximations, one of the contributions of this investigation. On the basis of laboratory analysis, well-field tests, and general information on rates of drainage from saturated sediments, a gravity yield of 14 percent was assumed as a first approximation. The equation was then solved, by weeks, for evapotranspiration, ET. The evapotranspiration losses were plotted against the calendar week. Using the time of year as a control, a smooth curve was fitted to the evapotranspiration data, and modified values of ET were read from the curve. These were used to compute weekly values of the infiltration residual which were plotted against ground-water stage. The slope of the line of best fit gave a closer approximation of gravity yield, Yg. The process was repeated. The approximations converged, so that a fourth and final approximation resulted in a close grouping of all the points along a line whose slope indicated a Yg of 11.0 percent, and a slightly asymmetric bell-shaped curve of total evapotranspiration by weeks was obtained that is considered representative of this area. Check calculations of gravity yield were made during periods of low evapotranspiration and high infiltration, which substantiate the computed average of 11.0 percent. Refinements in the method of deriving the ground-water budget were introduced to supplement the techniques developed by Meinzer and Stearns in the study of the Pomperaug River basin in Connecticut in 1913 and 1916. The hydrologic equation for the ground-water cycle may be written Gr=D + delta H. x Yg + ETg, in which Gr is ground-water recharge (infiltration); D is ground-water drainage; delta H is the change in mean ground-water stage (final stage minus initial stage); Yg is gravity yield (taken as 11.0 percent in computations here); an

Water Supply Paper

A novel approach to assessing natural resource injury with Bayesian networks

Quantifying the effects of environmental stressors on natural resources is problematic because of complex interactions among environmental factors that influence endpoints of interest. This complexity, coupled with data limitations, propagates uncertainty that can make it difficult to causally associate specific environmental stressors with injury endpoints. The Natural Resource Damage Assessment and Restoration (NRDAR) regulations under the Comprehensive Environmental Response, Compensation, and Liability Act and Oil Pollution Act aim to restore natural resources injured by oil spills and hazardous substances released into the environment; exploration of alternative statistical methods to evaluate effects could help address NRDAR legal claims. Bayesian networks (BNs) are statistical tools that can be used to estimate the influence and interrelatedness of abiotic and biotic environmental variables on environmental endpoints of interest. We investigated the application of a BN for injury assessment using a hypothetical case study by simulating data of acid mine drainage (AMD) affecting a fictional stream-dwelling bird species. We compared the BN-generated probability estimates for injury with a more traditional approach using toxicity thresholds for water and sediment chemistry. Bayesian networks offered several distinct advantages over traditional approaches, including formalizing the use of expert knowledge, probabilistic estimates of injury using intermediate direct and indirect effects, and the incorporation of a more nuanced and ecologically relevant representation of effects. Given the potential that BNs have for natural resource injury assessment, more research and field-based application are needed to determine their efficacy in NRDAR. We expect the resulting methods will be of interest to many US federal, state, and tribal programs devoted to the evaluation, mitigation, remediation, and/or restoration of natural resources injured by releases or spills of contaminants

Integrated Environmental Assessment and Management

Digital representation of oil and natural gas well pad scars in southwest Wyoming

The recent proliferation of oil and natural gas energy development in southwest Wyoming has stimulated the need to understand wildlife responses to this development. Central to many wildlife assessments is the use of geospatial methods that rely on digital representation of energy infrastructure. Surface disturbance of the well pad scars associated with oil and natural gas extraction has been an important but unavailable infrastructure layer. To provide a digital baseline of this surface disturbance, we extracted visible oil and gas well pad scars from 1-meter National Agriculture Imagery Program imagery (NAIP) acquired in 2009 for a 7.7 million-hectare region of southwest Wyoming. Scars include the pad area where wellheads, pumps, and storage facilities reside, and the surrounding area that was scraped and denuded of vegetation during the establishment of the pad. Scars containing tanks, compressors, and the storage of oil and gas related equipment, and produced-water ponds were also collected on occasion. Our extraction method was a two-step process starting with automated extraction followed by manual inspection and clean up. We used available well-point information to guide manual clean up and to derive estimates of year of origin and duration of activity on a pad scar. We also derived estimates of the proportion of non-vegetated area on a scar using a Normalized Difference Vegetation Index derived using 1-meter NAIP imagery. We extracted 16,973 pad scars of which 15,318 were oil and gas well pads. Digital representation of pad scars along with time-stamps of activity and estimates of non-vegetated area provides important baseline (circa 2009) data for assessments of wildlife responses, land-use trends, and disturbance-mediated pattern assessments.

Wyoming

Spatial capture–recapture for categorically marked populations with an application to genetic capture–recapture

Recently introduced unmarked spatial capture–recapture (SCR), spatial mark–resight (SMR), and 2‐flank spatial partial identity models (SPIMs) extend the domain of SCR to populations or observation systems that do not always allow for individual identity to be determined with certainty. For example, some species do not have natural marks that can reliably produce individual identities from photographs, and some methods of observation produce partial identity samples as is the case with remote cameras that sometimes produce single‐flank photographs. Unmarked SCR, SMR, and SPIM share the feature that they probabilistically resolve the uncertainty in individual identity using the spatial location where samples were collected. Spatial location is informative of individual identity in spatially structured populations because a sample is more likely to have been produced by an individual living near the trap where it was recorded than an individual living further away from the trap. Further, the level of information about individual identity that a spatial location contains is related to two key ecological concepts, population density and home range size, which we quantify using a proposed Identity Diversity Index (IDI). We show that latent and partial identity SCR models produce imprecise and biased density estimates in many high IDI scenarios when data are sparse. We then extend the unmarked SCR model to incorporate categorical, partially identifying covariates, which reduce the level of uncertainty in individual identity, increasing the reliability and precision of density estimates, and allowing reliable density estimation in scenarios with higher IDI values and with more sparse data. We illustrate the performance of this “categorical SPIM” via simulations and by applying it to a black bear data set using microsatellite loci as categorical covariates, where we reproduce the full data set estimates with only slightly less precision using fewer loci than necessary for confident individual identification. We then discuss how the categorical SPIM can be applied to other wildlife sampling scenarios such as remote camera surveys, where natural or researcher‐applied partial marks can be observed in photographs. Finally, we discuss how the categorical SPIM can be added to SMR, 2‐flank SPIM, or other latent identity SCR models.

Ecosphere

Post-wildfire curve number estimates for the southern Rocky Mountains in Colorado, USA

The curve number method first developed by the US Department of Agriculture Soil Conservation Service (now the Natural Resources Conservation Service) is often used for post-wildfire runoff assessments. These assessments are critical for land and emergency managers making decisions on life and property risks following a wildfire event. Three approaches (i.e., historical event observations, linear regression model, and regression tree model) were used to help estimate a post-wildfire curve number from watershed and wildfire parameters. For the first method, we used runoff events from 102 burned watersheds in Colorado, southern Wyoming, northern New Mexico, and eastern Utah to quantify changes in curve number values from pre- to post-wildfire conditions. The curve number changes from the measured runoff events vary substantially between positive and negative values. The measured curve number changes were then associated with watershed characteristics (e.g., slope, elevation, northness, and eastness) and land cover type to develop prediction models that provide estimates of post-wildfire curve number changes. Finally, we used a regression tree method to demonstrate that accurate predications can be developed using the measured curve number changes from our study domain. These models can be used for future post-wildfire assessments within the region.

Colorado, New Mexico, Utah, Wyoming

A comparison of methods for streamflow uncertainty estimation

Streamflow time series are commonly derived from stage‐discharge rating curves, but the uncertainty of the rating curve and resulting streamflow series are poorly understood. While different methods to quantify uncertainty in the stage‐discharge relationship exist, there is limited understanding of how uncertainty estimates differ between methods due to different assumptions and methodological choices. We compared uncertainty estimates and stage‐discharge rating curves from seven methods at three river locations of varying hydraulic complexity. Comparison of the estimated uncertainties revealed a wide range of estimates, particularly for high and low flows. At the simplest site on the Isère River (France), full width 95% uncertainties for the different methods ranged from 3 to 17% for median flows. In contrast, uncertainties were much higher and ranged from 41 to 200% for high flows in an extrapolated section of the rating curve at the Mahurangi River (New Zealand) and 28 to 101% for low flows at the Taf River (United Kingdom), where the hydraulic control is unstable at low flows. Differences between methods result from differences in the sources of uncertainty considered, differences in the handling of the time‐varying nature of rating curves, differences in the extent of hydraulic knowledge assumed, and differences in assumptions when extrapolating rating curves above or below the observed gaugings. Ultimately, the selection of an uncertainty method requires a match between user requirements and the assumptions made by the uncertainty method. Given the significant differences in uncertainty estimates between methods, we suggest that a clear statement of uncertainty assumptions be presented alongside streamflow uncertainty estimates.

Water Resources Research

River meander modeling of the Wabash River near the Interstate 64 Bridge near Grayville, Illinois

Natural river channels continually evolve and change shape over time. As a result, channel evolution or migration can cause problems for bridge structures that are fixed in the flood plain. A once-stable bridge structure that was uninfluenced by a river’s shape could be encroached upon by a migrating river channel. The potential effect of the actively meandering Wabash River on the Interstate 64 Bridge at the border with Indiana near Grayville, Illinois, was studied using a river migration model called RVR Meander. RVR Meander is a toolbox that can be used to model river channel meander migration with physically based bank erosion methods. This study assesses the Wabash River meandering processes through predictive modeling of natural meandering over the next 100 years, climate change effects through increased river flows, and bank protection measures near the Interstate 64 Bridge.

Illinois

Isotopic tracing of perchlorate in the environment

Isotopic measurements can be used for tracing the sources and behavior of environmental contaminants. Perchlorate (ClO 4 − ) has been detected widely in groundwater, soils, fertilizers, plants, milk, and human urine since 1997, when improved analytical methods for analyzing ClO 4 − concentration became available for routine use. Perchlorate ingestion poses a risk to human health because of its interference with thyroidal hormone production. Consequently, methods for isotopic analysis of ClO 4 − have been developed and applied to assist evaluation of the origin and migration of this common contaminant. Isotopic data are now available for stable isotopes of oxygen and chlorine, as well as 36 Cl isotopic abundances, in ClO 4 − samples from a variety of natural and synthetic sources. These isotopic data provide a basis for distinguishing sources of ClO 4 − found in the environment, and for understanding the origin of natural ClO 4 − . In addition, the isotope effects of microbial ClO 4 − reduction have been measured in laboratory and field experiments, providing a tool for assessing ClO 4 − attenuation in the environment. Isotopic data have been used successfully in some areas for identifying major sources of ClO 4 − contamination in drinking water supplies. Questions about the origin and global biogeochemical cycle of natural ClO 4 − remain to be addressed; such work would benefit from the development of methods for preparation and isotopic analysis of ClO 4 − in samples with low concentrations and complex matrices.

Book chapter

Leveraging geodetic data to reduce losses from earthquakes

Seismic hazard assessments that are based on a variety of data and the best available science, coupled with rapid synthesis of real-time information from continuous monitoring networks to guide post-earthquake response, form a solid foundation for effective earthquake loss reduction. With this in mind, the Earthquake Hazards Program (EHP) of the U.S. Geological Survey (USGS) Natural Hazards Mission Area (NHMA) engages in a variety of undertakings, both established and emergent, in order to provide high quality products that enable stakeholders to take action in advance of and in response to earthquakes. Examples include the National Seismic Hazard Model (NSHM), development of tools for improved situational awareness such as earthquake early warning (EEW) and operational earthquake forecasting (OEF), research about induced seismicity, and new efforts to advance comprehensive subduction zone science and monitoring. Geodetic observations provide unique and complementary information directly relevant to advancing many aspects of these efforts (fig. 1). EHP scientists have long leveraged geodetic data for a range of influential studies, and they continue to develop innovative observation and analysis methods that push the boundaries of the field of geodesy as applied to natural hazards research. Given the ongoing, rapid improvement in availability, variety, and precision of geodetic measurements, considering ways to fully utilize this observational resource for earthquake loss reduction is timely and essential. This report presents strategies, and the underlying scientific rationale, by which the EHP could achieve the following outcomes: The EHP is an authoritative source for the interpretation of geodetic data and its use for earthquake loss reduction throughout the United States and its territories. The USGS consistently provides timely, high quality geodetic data to stakeholders. Significant earthquakes are better characterized by incorporating geodetic data into USGS event response products and by expanded use of geodetic imaging data to assess fault rupture and source parameters. Uncertainties in the NSHM, and in regional earthquake models, are reduced by fully incorporating geodetic data into earthquake probability calculations. Geodetic networks and data are integrated into the operations and earthquake information products of the Advanced National Seismic System (ANSS). Earthquake early warnings are improved by more rapidly assessing ground displacement and the dynamic faulting process for the largest earthquakes using real-time geodetic data. Methodology for probabilistic earthquake forecasting is refined by including geodetic data when calculating evolving moment release during aftershock sequences and by better understanding the implications of transient deformation for earthquake likelihood. A geodesy program that encompasses a balanced mix of activities to sustain missioncritical capabilities, grows new competencies through the continuum of fundamental to applied research, and ensures sufficient resources for these endeavors provides a foundation by which the EHP can be a leader in the application of geodesy to earthquake science. With this in mind the following objectives provide a framework to guide EHP efforts: Fully utilize geodetic information to improve key products, such as the NSHM and EEW, and to address new ventures like the USGS Subduction Zone Science Plan. Expand the variety, accuracy, and timeliness of post-earthquake information products, such as PAGER (Prompt Assessment of Global Earthquakes for Response), through incorporation of geodetic observations. Determine if geodetic measurements of transient deformation can significantly improve estimates of earthquake probability. Maintain an observational strategy aligned with the target outcomes of this document that includes continuous monitoring, recording of ephemeral observations, focused data collection for use in research, and application-driven data processing and analysis systems. Collaborate on research, development, and operation of affordable, high-precision seafloor geodetic methods that improve earthquake forecasting and event response. Advance computational techniques and instrumentation to enable use of strategies like repeat-pass imagery and low-cost geodetic sensors for earthquake response, monitoring, and research. Engage stakeholders and collaborate with partner institutions to foster operational and research objectives and to safeguard the continued health of geodetic infrastructure upon which we mutually depend. Maintaining a vibrant internal research program provides the foundation by which the EHP can remain an effective and trusted source for earthquake science. Exploiting abundant new data sources, evaluating and assimilating the latest science, and pursuing novel avenues of investigation are means to fulfilling the EHP’s core responsibilities and realizing the important scientific advances envisioned by its scientists. Central to the success of such a research program is engaging personnel with a breadth of competencies and a willingness and ability to adapt these to the program’s evolving priorities, enabling current staff to expand their skills and responsibilities, and planning holistically to meet shared workforce needs. In parallel, collaboration with external partners to support scientific investigations that complement ongoing internal research enables the EHP to strengthen earthquake information products by incorporating alternative perspectives and approaches and to study topics and geographic regions that cannot be adequately covered internally. With commensurate support from technical staff who possess diverse skills, including engineering, information technology, and proficiency in quantitative analysis combined with basic geophysical knowledge, the EHP can achieve the geodetic outcomes identified in this document.

Open-File Report

Connecting diverse disciplines to improve understanding of surface water-groundwater interactions

Laura K. Lautz is a premier mentor, collaborator, and researcher at the intersection of natural hydrologic systems and humans. Her research has shifted the paradigm around measuring and understanding the impacts of surface water and groundwater interactions across spatial and temporal scales. She has done this by testing and refining new methods and by collaborating with, training, supporting, and mentoring diverse scientists. Here, we review her research across five themes, summarizing the prior status of the field, what Lautz contributed, as well as new directions in the field inspired by her work. Lautz’s research expanded our understanding of the impacts of stream restoration on surface water-groundwater interactions, where she tested new field methods and showed that restoration structures increase hyporheic exchange, locally altering biogeochemical function of the streambed. She refined novel methods for measuring surface water-groundwater exchanges and worked to make these methods easily accessible through freely available software. Her research group greatly expanded the use of heat as a quantitative tracer of hydrologic processes via the well-used VFLUX and HFLUX programs. Her research evaluated the impacts of surface water-groundwater interactions in urban streams, showing the substantial fluxes of nutrients and chloride that can move through those exchanges and the potential for groundwater to help buffer contamination. To assess groundwater impacts on streamflow below tropical glaciers, she used a wide range of field methods to reveal the sensitivity of these systems to climate change. Finally, she built tools to quantify natural brine contamination of drinking water wells in areas that may later be subject to high-volume hydraulic fracturing, creating a needed ‘pre-fracking’ dataset. Through this process, she identified multiple sources of salinity that are already reaching wells in these systems. Overall, this research has been done with a focus on mentoring and training the next generation of hydrologists, including work to specifically train for careers beyond academia, and facilitating early career scientists to realize their innate potentials. With former trainees in careers across industry, government, and academia, Dr. Laura K. Lautz is now working to build cross-disciplinary research at even larger scales, across federal research units, guaranteeing that an even larger impact on hydrology is still to come.

Journal of Hydrology X

On the use of log-transformation vs. nonlinear regression for analyzing biological power laws

Power‐law relationships are among the most well‐studied functional relationships in biology. Recently the common practice of fitting power laws using linear regression (LR) on log‐transformed data has been criticized, calling into question the conclusions of hundreds of studies. It has been suggested that nonlinear regression (NLR) is preferable, but no rigorous comparison of these two methods has been conducted. Using Monte Carlo simulations, we demonstrate that the error distribution determines which method performs better, with NLR better characterizing data with additive, homoscedastic, normal error and LR better characterizing data with multiplicative, heteroscedastic, lognormal error. Analysis of 471 biological power laws shows that both forms of error occur in nature. While previous analyses based on log‐transformation appear to be generally valid, future analyses should choose methods based on a combination of biological plausibility and analysis of the error distribution. We provide detailed guidelines and associated computer code for doing so, including a model averaging approach for cases where the error structure is uncertain.

Ecology

Constructing an interdisciplinary flow regime recommendation

It is generally agreed that river rehabilitation most often relies on restoring a more natural flow regime, but credibly defining the desired regime can be problematic. I combined four distinct methods to develop and refine month-by-month and event-based flow recommendations to protect and partially restore the ecological integrity of the Cache la Poudre River through Fort Collins, Colorado. A statistical hydrologic approach was used to summarize the river's natural flow regime and set provisional monthly flow targets at levels that were historically exceeded 75% of the time. These preliminary monthly targets were supplemented using results from three Poudre-specific disciplinary studies. A substrate maintenance flow model was used to better define the high flows needed to flush accumulated sediment from the river's channel and help sustain the riparian zone in this snowmelt-dominated river. A hydraulic/habitat model and a water temperature model were both used to better define the minimum flows necessary to maintain a thriving cool water fishery. The result is a range of recommended monthly flows and daily flow guidance illustrating the advantage of combining a wide range of available disciplinary information, supplemented by judgment based on ecological principles and a general understanding of river ecosystems, in a highly altered, working river. ?? 2010 American Water Resources Association.

Journal of the American Water Resources Associatio

Methods to characterize environmental settings of stream and groundwater sampling sites for National Water-Quality Assessment

Characterization of natural and anthropogenic features that define the environmental settings of sampling sites for streams and groundwater, including drainage basins and groundwater study areas, is an essential component of water-quality and ecological investigations being conducted as part of the U.S. Geological Survey's National Water-Quality Assessment program. Quantitative characterization of environmental settings, combined with physical, chemical, and biological data collected at sampling sites, contributes to understanding the status of, and influences on, water-quality and ecological conditions. To support studies for the National Water-Quality Assessment program, a geographic information system (GIS) was used to develop a standard set of methods to consistently characterize the sites, drainage basins, and groundwater study areas across the nation. This report describes three methods used for characterization-simple overlay, area-weighted areal interpolation, and land-cover-weighted areal interpolation-and their appropriate applications to geographic analyses that have different objectives and data constraints. In addition, this document records the GIS thematic datasets that are used for the Program's national design and data analyses.

Scientific Investigations Report

Comparison of U.S. Geological Survey and Ohio Environmental Protection Agency fish-collection methods using the index of biotic integrity and modified index of well-being, 1996–97

The U.S. Geological Survey (USGS) and Ohio Environmental Protection Agency (OEPA) collected data on fish from 10 stream sites in 1996 and 3 stream sites in 1997 as part of a comparative study of fish community assessment methods. The sites sampled represent a wide range of basin sizes (ranging from 132–6,330 square kilometers) and surrounding land-use types (urban, agricultural, and mixed). Each agency used its own fish-sampling protocol. Using the Index of Biotic Integrity and Modified Index of Well-Being, differences between data sets were tested for significance by means of the Wilcoxon signed-ranks test ( α = 0.05). Results showed that the median of Index of Biotic Integrity differences between data sets was not significantly different from zero ( p = 0.2521); however, the same statistical test showed the median differences in the Modified Index of Well-Being scores to be significantly different from zero ( p = 0.0158). The differences observed in the Index of Biotic Integrity scores are likely due to natural variability, increased variability at sites with degraded water quality, differences in sampling methods, and low-end adjustments in the Index of Biotic Integrity calculation when fewer than 50 fish were collected. The Modified Index ofWell-Being scores calculated by OEPA were significantly higher than those calculated by the USGS. This finding was attributed to the comparatively large numbers and biomass of fish collected by the OEPA. By combining the two indices and viewing them in terms of the percentage attainment of Ohio Warmwater Habitat criteria, the two agencies’ data seemed comparable, although the Index of Biotic Integrity scores were more similar than the Modified Index of Well-Being scores.

Michigan, New York, Ohio, Pennsylvania

Methods of analysis by the U.S. Geological Survey Organic Geochemistry Research Group: Determination of acetamide herbicides and their degradation products in water using online solid-phase extraction and liquid chromatography/mass spectrometry

An analytical method for the determination of 6 acetamide herbicides (acetochlor, alachlor, dimethenamid, flufenacet, metolachlor, and propachlor) and 16 of their degradation products in natural water samples using solid-phase extraction and liquid chromatography/mass spectrometry is described in this report. Special consideration was given during the development of the method to prevent the formation of degradation products during the analysis. Filtered water samples were analyzed using octadecylsilane as the solid-phase extraction media on online automated equipment followed by liquid chromatography/mass spectrometry. The method uses only 10 milliliters of sample per injection. Three different water-sample matrices, a reagent-water, a ground-water, and a surface-water sample spiked at 0.10 and 1.0 microgram per liter, were analyzed to determine method performance. Method detection limits ranged from 0.004 to 0.051 microgram per liter for the parent acetamide herbicides and their degradation products. Mean recoveries for the acetamide compounds in the ground- and surface-water samples ranged from 62.3 to 117.4 percent. The secondary amide of acetochlor/metolachlor ethanesulfonic acid (ESA) was recovered at an average rate of 43.5 percent. The mean recoveries for propachlor and propachlor oxanilic acid (OXA) were next lowest, ranging from 62.3 to 95.5 percent. Mean recoveries from reagent-water samples ranged from 90.3 to 118.3 percent for all compounds. Overall the mean of the mean recoveries of all compounds in the three matrices spiked at 0.10 and 1.0 microgram per liter ranged from 89.9 to 100.7 percent, including the secondary amide of acetochlor/metolachlor ESA and the propachlor compounds. The acetamide herbicides and their degradation products are reported in concentrations ranging from 0.05 to 2.0 micrograms per liter. The upper concentration limit is 2.0 micrograms per liter for all compounds without dilution. With the exception of the secondary amide of acetochlor/metolachlor ESA, good precision and accuracy for the chloroacetanalide herbicides and their degradation compounds were demonstrated for the method in buffered reagent water, ground water, and surface water. The extraction method as used did not optimize the recovery of the secondary amide of acetochlor/metolachlor ESA.

Open-File Report

Groundwater-flow and land-subsidence model of Antelope Valley, California

Antelope Valley, California, is a topographically closed basin in the western part of the Mojave Desert, about 50 miles northeast of Los Angeles. The Antelope Valley groundwater basin is about 940 square miles and is separated from the northern part of Antelope Valley by faults and low-lying hills. Prior to 1972, groundwater provided more than 90 percent of the total water supply in the valley; since 1972, it has provided between 50 and 90 percent. Most groundwater pumping in the valley occurs in the Antelope Valley groundwater basin, which includes the rapidly growing cities of Lancaster and Palmdale. Groundwater-level declines of more than 270 feet in some parts of the groundwater basin have resulted in an increase in pumping lifts, reduced well efficiency, and land subsidence of more than 6 feet in some areas. Future urban growth and limits on the supply of imported water may increase reliance on groundwater. In 2011, the Los Angeles County Superior Court of California ruled that the Antelope Valley groundwater basin is in overdraft—groundwater extractions are in excess of the Court-defined safe yield of the groundwater basin. The Court determined that the safe yield of the adjudicated area of the basin was 110,000 acre-feet per year (acre-ft/yr). Natural recharge is an important component of total groundwater recharge in Antelope Valley; however, the exact quantity and distribution of natural recharge, primarily in the form of mountain-front recharge, is uncertain, with total estimates ranging from 30,000 to 160,000 acre-ft/yr. Technical experts, retained by parties to the adjudication, used 60,000 acre-ft/yr to estimate the sustainable yield of the basin, and this value was used in this study. In order to better understand the uncertainty associated with natural recharge and to provide a tool to aid in groundwater management, a numerical model of groundwater flow and land subsidence in the Antelope Valley groundwater basin was developed using old and new geohydrologic information. The groundwater-flow system consists of three aquifers: the upper, middle, and lower aquifers. The three aquifers, which were identified on the basis of the hydrologic properties, age, and depth of the unconsolidated deposits, consist of gravel, sand, silt, and clay alluvial deposits and clay and silty clay lacustrine deposits. Prior to groundwater development in the valley, recharge was primarily the infiltration of runoff from the surrounding mountains. Groundwater flowed from the recharge areas to discharge areas around the playas where it discharged from the aquifer system as either evapotranspiration or from springs. Partial barriers to horizontal groundwater flow, such as faults, have been identified in the groundwater basin. Water-level declines owing to groundwater development have eliminated the natural sources of discharge, and pumping for agricultural and urban uses have become the primary source of discharge from the groundwater system. Infiltration of return flow from agricultural irrigation has become an important source of recharge to the aquifer system. The groundwater-flow model of the basin was discretized horizontally into a grid of 130 rows and 118 columns of square cells 1 kilometer (0.621 mile) on a side, and vertically into four layers representing the upper (two layers), middle (one layer), and lower (one layer) aquifers. Faults that were thought to act as horizontal-flow barriers were simulated in the model. The model was calibrated to simulate steady-state conditions, represented by 1915 water levels and transient-state conditions during 1915–95, by using water-level and subsidence data. Initial estimates of the aquifer-system properties and stresses were obtained from a previously published numerical model of the Antelope Valley groundwater basin; estimates also were obtained from recently collected hydrologic data and from results of simulations of groundwater-flow and land-subsidence models of the Edwards Air Force Base area. Some of these initial estimates were modified during model calibration. Groundwater pumpage for agriculture was estimated on the basis of irrigated crop acreage and crop consumptive-use data. Pumpage for public supply, which is metered, was compiled and entered into a database used for this study. Estimated annual agricultural pumpage peaked at 395,000 acre-feet (acre-ft) in 1951 and then declined because of declining agricultural production. Recharge from irrigation return flows was assumed to be 30 percent of agricultural pumpage; delays associated with return flow moving through the unsaturated zone were also simulated. The annual quantity of mountain-front recharge initially was based on estimates from previous studies. The model was calibrated using the PEST software suite; prior information from the area was incorporated through the use of Tikhonov regularization. During model calibration, the estimated mountain-front recharge was reduced from the previous estimate of 30,300 acre-ft/yr to 29,150 acre-ft/yr. Results of the simulations using the calibrated model indicate that simulated groundwater pumpage exceeded recharge in most years, resulting in an estimated cumulative depletion in groundwater storage of 8,700,000 acre-ft during the transient-simulation period (1915–2005). About 15,000,000 acre-ft of cumulative groundwater pumpage was simulated during the transient-simulation period (1915–2005), reaching a maximum rate of about 400,000 acre-ft/yr in 1951. Groundwater pumpage resulted in simulated hydraulic heads declining by more than 150 feet (ft) compared to 1915 conditions in agricultural areas. The decline in hydraulic head in the groundwater basin is the result of this depletion of groundwater storage. In turn, the simulated decline in hydraulic head in the groundwater basin has resulted in the decrease in natural discharge from the basin and has caused compaction of aquitards, resulting in land subsidence. The areal distribution of total simulated land subsidence for 2005, after about 90 years of groundwater development, indicates that land subsidence occurred throughout almost the entire Lancaster subbasin, with a maximum of about 9.4 ft in the central and eastern parts of the subbasin. An important objective of this study was to systematically address the uncertainty in estimates of natural recharge and related aquifer parameters by using the groundwater-flow and land-subsidence model with observational data and expert knowledge. After the model was calibrated to the observations and a reasonable parameter set obtained, the parameter null space—parameter values that do not appreciably affect the model calibration but may have importance for prediction—was identified. The effect of parameter uncertainty on the estimation of mountain-front recharge was addressed using the Null-Space Monte Carlo method. The Pareto trade-off method of visualizing uncertainty was also used to portray the reasonableness of larger natural-recharge rates. Results indicate that the total mountain-front recharge likely ranges between 28,000 and 44,000 acre-ft/yr, which is appreciably less than published estimates of 60,000 acre-ft/yr. Additionally, expected errors associated with agricultural pumpage estimates used in this study were found to have relatively little effect on the estimates of mountain-front recharge, reflecting the difficulty in increasing recharge through manipulation of other components of the water budget. The calibrated model was used to simulate the response of the aquifer to potential future pumping scenarios: (1) no change in the distribution of pumpage, or status quo; (2) redistribution of pumpage; and (3) artificial recharge. All three of these scenarios specify a total pumpage throughout the Antelope Valley of 110,000 acre-ft/yr according to the safe yield value ruled by the Los Angeles County Superior Court of California. This reduction in groundwater pumpage is assumed uniform throughout the basin, based on a 10-percent reduction of the total pumpage in 2005 to achieve the 110,000 acre-ft/yr level. The calibrated Antelope Valley groundwater-flow and land-subsidence model was used to simulate the hydrologic effects of the three groundwater-management scenarios during a 50-year period by using the reduced, temporally constant, pumpage distribution. Results from the first scenario indicated that the total drawdown observed since predevelopment would continue, with values exceeding 325 ft near Palmdale; consequently, land subsidence would also continue, with additional subsidence (since 2005) exceeding 3 ft in the central part of the Lancaster subbasin. The second scenario evaluated redistributing pumpage from areas in the Lancaster subbasin where simulated hydraulic-head declines were the greatest to areas where declines were smallest. Neither a formal optimization algorithm nor water-rights allocations were considered when redistributing the pumpage. Results indicated that hydraulic heads near Palmdale, where the pumpage was reduced, would recover by about 200 ft compared to 2005 conditions, with only 30 ft of additional drawdown in the northwestern part of the Lancaster subbasin, where the pumpage was increased. The magnitude of the simulated additional land subsidence decreased slightly compared to the first, status quo, scenario but land subsidence continued to be simulated throughout most of the northern part of the Lancaster subbasin. The third scenario consisted of two artificial-recharge simulations along the Upper Amargosa Creek channel and at a site located north of Antelope Buttes. Results indicate that applying artificial recharge at these sites would yield continued drawdowns and associated land subsidence. However, the magnitudes of drawdown and subsidence would be smaller than those simulated in the status quo scenario, indicating that artificial-recharge operations in the Antelope Valley could be expected to reduce the magnitude and extent of continued water-level declines and associated land subsidence.

California

Increased scientific rigor will improve reliability of research and effectiveness of management

Rigorous science that produces reliable knowledge is critical to wildlife management because it increases accurate understanding of the natural world and informs management decisions effectively. Application of a rigorous scientific method based on hypothesis testing minimizes unreliable knowledge produced by research. To evaluate the prevalence of scientific rigor in wildlife research, we examined 24 issues of the Journal of Wildlife Management from August 2013 through July 2016. We found 43.9% of studies did not state or imply a priori hypotheses, which are necessary to produce reliable knowledge. We posit that this is due, at least in part, to a lack of common understanding of what rigorous science entails, how it produces more reliable knowledge than other forms of interpreting observations, and how research should be designed to maximize inferential strength and usefulness of application. Current primary literature does not provide succinct explanations of the logic behind a rigorous scientific method or readily applicable guidance for employing it, particularly in wildlife biology; we therefore synthesized an overview of the history, philosophy, and logic that define scientific rigor for biological studies. A rigorous scientific method includes 1) generating a research question from theory and prior observations, 2) developing hypotheses (i.e., plausible biological answers to the question), 3) formulating predictions (i.e., facts that must be true if the hypothesis is true), 4) designing and implementing research to collect data potentially consistent with predictions, 5) evaluating whether predictions are consistent with collected data, and 6) drawing inferences based on the evaluation. Explicitly testing a priori hypotheses reduces overall uncertainty by reducing the number of plausible biological explanations to only those that are logically well supported. Such research also draws inferences that are robust to idiosyncratic observations and unavoidable human biases. Offering only post hoc interpretations of statistical patterns (i.e., a posteriori hypotheses) adds to uncertainty because it increases the number of plausible biological explanations without determining which have the greatest support. Further, post hoc interpretations are strongly subject to human biases. Testing hypotheses maximizes the credibility of research findings, makes the strongest contributions to theory and management, and improves reproducibility of research. Management decisions based on rigorous research are most likely to result in effective conservation of wildlife resources.

Journal of Wildlife Management