Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Water Surface Techniques”

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 505 records · Page 28Linked to original sources

The surface elevation table and marker horizon technique: A protocol for measuring wetland elevation dynamics, narrative (Version 2.0)

The National Park Service (NPS), in response to the growing evidence and awareness of the effects of climate change on federal lands, determined that monitoring wetland elevation change is a top priority in North Atlantic Coastal parks (Stevens et al. 2010). As a result, the NPS Northeast Coastal and Barrier Network (NCBN) in collaboration with colleagues from the U.S. Geological Survey (USGS) and the National Oceanic and Atmospheric Administration (NOAA) have developed a protocol for monitoring wetland elevation change and other processes important for determining the viability of wetland communities. Although focused on North Atlantic Coastal parks, this document is applicable to all coastal and inland wetland regions. Wetlands exist within a narrow range of elevation that is influenced by local hydrologic conditions. For tidally influenced coastal wetlands, local hydrologic conditions may be changing as sea levels continue to rise. As sea level rises, coastal wetland systems may either a) build elevation to maintain favorable hydrologic conditions for their survival or b) become increasingly inundated beyond their physiological tolerance resulting in an eventual conversion to open water. A better understanding of these processes will help to determine the present and future viability of coastal wetlands managed by the NPS and can help address measures to ensure these communities exist into the future. This protocol provides the reader with instructions and guidelines on designing a monitoring plan or study to: Quantify elevation change in wetlands with the Surface Elevation Table (SET); understand the processes that influence elevation change, including vertical accretion (SET and Marker Horizon methods); survey the wetland surface and SET mark to a common reference datum to allow for comparing sample stations to each other and to local tidal datums; survey the SET mark to monitor its relative stability; and establish water level equipment to characterize the hydrology and tidal datums for a particular wetland site or sites. The protocol is divided into two documents: The narrative (this document) presents an overview of all aspects of monitoring wetland elevation dynamics, and a document of Standard Operating Procedures (SOPs) with detailed instructions on the design, installation, data collection, and data management in support of monitoring changes in wetland elevation.

Science Report↗

Parametric study of the physical properties of hydrate‐bearing sand, silt, and clay sediments: 2. Small‐strain mechanical properties

The small‐strain mechanical properties (e.g., seismic velocities) of hydrate‐bearing sediments measured under laboratory conditions provide reference values for calibration of logging and seismic exploration results acquired in hydrate‐bearing formations. Instrumented cells were designed for measuring the compressional (P) and shear (S) velocities of sand, silts, and clay with and without hydrate and subject to vertical effective stresses of 0.01 to 2 MPa. Tetrahydrofuran (THF), which is fully miscible in water, was used as the hydrate former to permit close control over the hydrate saturation S hyd and to produce hydrate from dissolved phase, as methane hydrate forms in most natural marine settings. The results demonstrate that laboratory hydrate formation technique controls the pattern of P and S velocity changes with increasing S hyd and that the small‐strain properties of hydrate‐bearing sediments are governed by effective stress, σ ′ v and sediment specific surface. The S velocity increases with hydrate saturation owing to an increase in skeletal shear stiffness, particularly when hydrate saturation exceeds S hyd ≈ 0.4. At very high hydrate saturations, the small strain shear stiffness is determined by the presence of hydrates and becomes insensitive to changes in effective stress. The P velocity increases with hydrate saturation due to the increases in both the shear modulus of the skeleton and the bulk modulus of pore‐filling phases during fluid‐to‐hydrate conversion. Small‐strain Poisson's ratio varies from 0.5 in soft sediments lacking hydrates to 0.25 in stiff sediments (i.e., subject to high vertical effective stress or having high S hyd ). At S hyd ≥ 0.5, hydrate hinders expansion and the loss of sediment stiffness during reduction of vertical effective stress, meaning that hydrate‐rich natural sediments obtained through pressure coring should retain their in situ fabric for some time after core retrieval if the cores are maintained within the hydrate stability field.

Journal of Geophysical Research↗

Groundwater geochemistry in the Seminole Well Field, Cedar Rapids, Iowa

The City of Cedar Rapids obtains its municipal water supply from four well fields in an alluvial aquifer along the Cedar River in east-central Iowa. Since 1992, the City and the U.S. Geological Survey have cooperatively studied the groundwater-flow system and water chemistry near the well fields. The geochemistry in the alluvial aquifer near the Seminole Well Field was assessed to identify potentially reactive minerals and possible chemical reactions that produce observed changes in water chemistry. Calcite, dolomite, ferrihydrite, quartz, rhodochrosite, and siderite were identified as potentially reactive minerals by calculating saturation indexes. Aluminosiicate minerals including albite, Ca-montmorillonite, gibbsite, illite, K-feldspar, and kaolinite were identified as potentially reactive minerals using hypothetical saturation indexes calculated with an assumed dissolved aluminum concentration of 1 microgram per liter. Balanced chemical equations derived from inverse-modeling techniques were used to assess chemical reactions as precipitation percolates to the water table. Calcite dissolution was predominate, but aluminosilicate weathering, cation exchange, and redox reactions also likely occurred. Microbial-catalyzed redox reactions altered the chemical composition of water infiltrating from the Cedar River into the alluvial aquifer by consuming dissolved oxygen, reducing nitrate, and increasing dissolved iron and manganese concentrations. Nitrate reduction only occurred in relatively shallow (3 to 7 meters below land surface) groundwater near the Cedar River and did not occur in water infiltrating to deeper zones of the alluvial aquifer.

Iowa↗

Introduction to Field Water-Quality Methods for the Collection of Metals - 2007 Project Summary

The U.S. Geological Survey (USGS), Region VI of the U.S. Environmental Protection Agency (USEPA), and the Osage Nation presented three 3-day workshops, in June-August 2007, entitled ?Introduction to Field Water-Quality Methods for the Collection of Metals.? The purpose of the workshops was to provide instruction to tribes within USEPA Region VI on various USGS surface-water measurement methods and water-quality sampling protocols for the collection of surface-water samples for metals analysis. Workshop attendees included members from over 22 tribes and pueblos. USGS instructors came from Oklahoma, New Mexico, and Georgia. Workshops were held in eastern and south-central Oklahoma and New Mexico and covered many topics including presampling preparation, water-quality monitors, and sampling for metals in surface water. Attendees spent one full classroom day learning the field methods used by the USGS Water Resources Discipline and learning about the complexity of obtaining valid water-quality and quality-assurance data. Lectures included (1) a description of metal contamination sources in surface water; (2) introduction on how to select field sites, equipment, and laboratories for sample analysis; (3) collection of sediment in surface water; and (4) utilization of proper protocol and methodology for sampling metals in surface water. Attendees also were provided USGS sampling equipment for use during the field portion of the class so they had actual ?hands-on? experience to take back to their own organizations. The final 2 days of the workshop consisted of field demonstrations of current USGS water-quality sample-collection methods. The hands-on training ensured that attendees were exposed to and experienced proper sampling procedures. Attendees learned integrated-flow techniques during sample collection, field-property documentation, and discharge measurements and calculations. They also used enclosed chambers for sample processing and collected quality-assurance samples to verify their techniques. Benefits of integrated water-quality sample-collection methods are varied. Tribal environmental programs now have the ability to collect data that are comparable across watersheds. The use of consistent sample collection, manipulation, and storage techniques will provide consistent quality data that will enhance the understanding of local water resources. The improved data quality also will help the USEPA better document the condition of the region?s water. Ultimately, these workshops equipped tribes to use uniform sampling methods and to provide consistent quality data that are comparable across the region.

Fact Sheet↗

Simulation of daily streamflows at gaged and ungaged locations within the Cedar River Basin, Iowa, using a Precipitation-Runoff Modeling System model

The U.S. Geological Survey, in cooperation with the Iowa Department of Natural Resources, conducted a study to examine techniques for estimation of daily streamflows using hydrological models and statistical methods. This report focuses on the use of a hydrologic model, the U.S. Geological Survey's Precipitation-Runoff Modeling System, to estimate daily streamflows at gaged and ungaged locations. The Precipitation-Runoff Modeling System is a modular, physically based, distributed-parameter modeling system developed to evaluate the impacts of various combinations of precipitation, climate, and land use on surface-water runoff and general basin hydrology. The Cedar River Basin was selected to construct a Precipitation-Runoff Modeling System model that simulates the period from January 1, 2000, to December 31, 2010. The calibration period was from January 1, 2000, to December 31, 2004, and the validation periods were from January 1, 2005, to December 31, 2010 and January 1, 2000 to December 31, 2010. A Geographic Information System tool was used to delineate the Cedar River Basin and subbasins for the Precipitation-Runoff Modeling System model and to derive parameters based on the physical geographical features. Calibration of the Precipitation-Runoff Modeling System model was completed using a U.S. Geological Survey calibration software tool. The main objective of the calibration was to match the daily streamflow simulated by the Precipitation-Runoff Modeling System model with streamflow measured at U.S. Geological Survey streamflow gages. The Cedar River Basin daily streamflow model performed with a Nash-Sutcliffe efficiency ranged from 0.82 to 0.33 during the calibration period, and a Nash-Sutcliffe efficiency ranged from 0.77 to -0.04 during the validation period. The Cedar River Basin model is meeting the criteria of greater than 0.50 Nash-Sutcliffe and is a good fit for streamflow conditions for the calibration period at all but one location, Austin, Minnesota. The Precipitation-Runoff Modeling System model accurately simulated streamflow at four of six uncalibrated sites within the basin. Overall, there was good agreement between simulated and measured seasonal and annual volumes throughout the basin for calibration and validation sites. The calibration period ranged from 0.2 to 20.8 percent difference, and the validation period ranged from 0.0 to 19.5 percent difference across all seasons and total annual runoff. The Precipitation-Runoff Modeling System model tended to underestimate lower streamflows compared to the observed streamflow values. This is an indication that the Precipitation-Runoff Modeling model needs more detailed groundwater and storage information to properly model the low-flow conditions in the Cedar River Basin.

Iowa↗

Occurrence, fate, and transport of aerially applied herbicides to control invasive buffelgrass within Saguaro National Park Rincon Mountain District, Arizona, 2015–18

The spread of the invasive and fire-adapted buffelgrass ( Cenchrus ciliaris L.) threatens desert ecosystems by competing for resources, increasing fuel loads, and creating wildfire connectivity. The Rincon Mountain District of Saguaro National Park addressed this natural resource threat with the use of glyphosate-based herbicides (GBHs). In 2010, the Rincon Mountain District initiated an aerial restoration plan to control dense buffelgrass patches in remote areas and implemented a trial project to evaluate the effects of aerial restoration techniques that included the helicopter application of GBHs. In 2014, more than 250 acres of buffelgrass in the Rincon Mountain District were treated with the aerial application of GBHs. This widespread aerial application of GBHs continued through 2018, but the potential transport and effects to aquatic ecosystems were unknown. In 2015–18, the U.S. Geological Survey, in cooperation with the National Park Service, studied the occurrence, distribution, fate, and transport of glyphosate in surface water and sediments derived from areas that were treated during past and current aerial herbicide applications. Three watersheds, treated with different regimens of GBHs, were sampled for glyphosate and the primary metabolite of glyphosate, aminomethylphosphonic acid (AMPA), during various hydrologic flow conditions. Water and aquatic sediment were collected from three watersheds, each in a different stage of application during the U.S. Geological Survey study. The unnamed watershed above the Loma Verde Trailhead referred to by the National Park Service as “Loma Verde canyon” had received no aerial treatment since 2014, whereas the Box Canyon watershed was aerially treated every year beginning in 2014. The Madrona Canyon watershed was first sprayed in 2016 and aerial application continued once a year though the entirety of the study. In addition, terrestrial soil samples were sampled from areas sprayed to understand dissipation rates and herbicide transport via sediments washing away during rainfall runoff. The concentrations present in water and sediment samples were compared to ecological benchmarks and characterized within the context of the environmental conditions of the park setting. Of the 48 water samples collected and analyzed for glyphosate and AMPA, 10.4 percent and 14.6 percent were detected above the laboratory minimum detection limit, respectively. Mean water concentrations, calculated using specific statistical methods for non-detects, were equal to the laboratory minimum detection limit of 0.02 microgram per liter for samples collected in all the watersheds. In aquatic sediments, glyphosate and AMPA were detected in 10.7 and 25.0 percent of the samples, whereas 89.5 and 100 percent of the terrestrial soil samples had detections for glyphosate and AMPA, respectively. Mean aquatic sediment concentrations were 1.13 and 4.42 micrograms per kilogram (μg/kg) for glyphosate and AMPA, respectively. Mean terrestrial soil concentrations were orders of magnitude greater than water and aquatic sediment with concentrations of 678 μg/kg for AMPA and 1,240 μg/kg for glyphosate. Hours after glyphosate-based herbicide was applied, the concentrations of glyphosate and AMPA were present in terrestrial soil samples near or above the laboratory maximum detection limit of 5,000 μg/kg. The Box Canyon watershed was the most intensively treated watershed in terms of total land area treated, total amount of GBH applied, and number of years treated. The frequent and large volume of treatment resulted in the highest number of detections of glyphosate and AMPA in water (3 and 7 detections, respectively) and in aquatic sediment (2 and 6 detections, respectively) samples. In comparison, the other two watersheds had two or fewer detections for glyphosate and AMPA in water and aquatic sediment. Glyphosate detected in pools was associated with increased rainfall closer in time to the last herbicide treatment. Glyphosate and AMPA concentration ratios above one, along with stable-isotope and tritium results, indicated that runoff processes were the primary transport mechanism for the two compounds when found in streams and pools rather than subsurface recharge or deeper flow paths. One pool in a small tributary of Box Canyon consistently had detections of glyphosate and AMPA in aquatic sediments, but these frequent concentrations were likely related to the intensive application upstream, near the steep terrain above the head of the channel that supplies the downstream pool. Intense flows during summer rainfall events move treated sediments into this channel where vegetation and the incised bedrock banks of the pool retained those sediments and ultimately led to frequent detections of both compounds. Isotope results in most of the pools and tinajas indicated that the water source had residence time representative of recently recharged waters, on the order of years. No water concentrations exceeded published criteria for human health or aquatic life. Median and maximum glyphosate and AMPA water concentrations were lower than those reported in other national assessments, but maximum concentrations observed in individual runoff samples were higher than median concentrations measured in the national assessments. A similar finding was observed with aquatic sediment concentrations measured in the Rincon Mountain District. Results from the study were compared and assessed in the context of other studies examining GBHs and their effects on amphibians, fish, and macroinvertebrates. This comparison was used to generalize the potential risk to aquatic species similar to those species in the Rincon Mountain District. Concentrations of published effect levels were several orders of magnitude greater than the highest concentration detected in water at the Rincon Mountain District. Most published studies evaluate acute and chronic toxicity for glyphosate and GBHs, and these criteria may not be representative of environmental conditions in the Rincon Mountain District. The classic lethal dose studies conducted in a controlled laboratory setting may not be suitable for comparison to the longer, variable, low-dose exposure conditions in the pools and tinajas in the Rincon Mountain District. However, this study determined that the fate of GBHs transported from treated areas to potential aquatic habitat was highly variable in occurrence, timing, and concentrations. This variability in glyphosate concentrations was too high, and the potential exposure was determined to be far too complex to directly compare with the results from controlled studies. This study provides the first information collected on GBHs used to control invasive buffelgrass in a remote, mountainous, and semiarid setting. The information about the transport and fate of herbicide application near aquatic habitat will help to inform managers about the broader ecosystem implications and provide useful information to other agencies implementing buffelgrass remediation strategies near aquatic habitat.

Arizona↗

Characterizing land surface phenology and exotic annual grasses in dryland ecosystems using Landsat and Sentinel-2 data in harmony

Invasive annual grasses, such as cheatgrass (Bromus tectorum L.), have proliferated in dryland ecosystems of the western United States, promoting increased fire activity and reduced biodiversity that can be detrimental to socio-environmental systems. Monitoring exotic annual grass cover and dynamics over large areas requires the use of remote sensing that can support early detection and rapid response initiatives. However, few studies have leveraged remote sensing technologies and computing frameworks capable of providing rangeland managers with maps of exotic annual grass cover at relatively high spatiotemporal resolutions and near real-time latencies. Here, we developed a system for automated mapping of invasive annual grass (%) cover using in situ observations, harmonized Landsat and Sentinel-2 (HLS) data, maps of biophysical variables, and machine learning techniques. A robust and automated cloud, cloud shadow, water, and snow/ice masking procedure (mean overall accuracy >81%) was implemented using time-series outlier detection and data mining techniques prior to spatiotemporal interpolation of HLS data via regression tree models (r = 0.94; mean absolute error (MAE) = 0.02). Weekly, cloud-free normalized difference vegetation index (NDVI) image composites (2016–2018) were used to construct a suite of spectral and phenological metrics (e.g., start and end of season dates), consistent with information derived from Moderate Resolution Image Spectroradiometer (MODIS) data. These metrics were incorporated into a data mining framework that accurately (r = 0.83; MAE = 11) modeled and mapped exotic annual grass (%) cover throughout dryland ecosystems in the western United States at a native, 30-m spatial resolution. Our results show that inclusion of weekly HLS time-series data and derived indicators improves our ability to map exotic annual grass cover, as compared to distribution models that use MODIS products or monthly, seasonal, or annual HLS composites as primary inputs. This research fills a critical gap in our ability to effectively assess, manage, and monitor drylands by providing a framework that allows for an accurate and timely depiction of land surface phenology and exotic annual grass cover at spatial and temporal resolutions that are meaningful to local resource managers.

California, Idaho, Nevada, Oregon↗

Chromium in minerals and selected aquifer materials

Between 1952 and 1964, hexavalent chromium, Cr(VI), was released into groundwater from a Pacific Gas and Electric Company (PG&E) compressor station in Hinkley, California, in the western Mojave Desert 80 miles northeast of Los Angeles, California. In 2015, the extent of anthropogenic Cr(VI) in groundwater in Hinkley and Water Valleys was uncertain, but some Cr(VI) in groundwater may be naturally occurring from rock and aquifer material. To evaluate potential sources of natural Cr(VI), chromium and other selected trace-element concentrations were measured by inductively coupled plasma-mass spectrometry (ICP-MS), with multi-acid digestion, on 34 samples of surficial alluvium and core material from Hinkley and Water Valleys, California, and on 2 samples of alluvium from the mafic Sheep Creek fan to the southwest. Chromium concentrations in Hinkley and Water Valleys ranged from 2 to 110 milligrams per kilogram (mg/kg), with a median concentration of 14 mg/kg; concentrations were highest in weathered mafic hornblende diorite associated with Iron Mountain. High chromium concentrations also were present within fine-textured materials and visually abundant iron- and manganese-oxide coatings on the surfaces of mineral grains. For comparison, chromium concentrations as high as 170 mg/kg were measured in mafic alluvium from the Sheep Creek fan. In contrast, chromium concentrations were lowest in Mojave-type deposits (Mojave River stream and lake margin deposits), with a median of 6 mg/kg. Chromium concentrations measured by ICP-MS compared favorably with concentrations measured by portable (handheld) X-ray fluorescence (pXRF; chapter B), on the basis of least-squares regression results and a coefficient of determination (R 2 ) of 0.97. Minerals in bulk samples and the heavy (dense) mineral fractions isolated from those samples were identified using optical techniques, X-ray diffraction (XRD), and scanning electron microscopy (SEM). Quartz and feldspar were the most abundant minerals, especially within recent and older Mojave River deposits. Chromium concentrations were as high as 1,250 mg/kg in the heavy-mineral fraction, with specific gravity greater than 3.32. Chromium was not commonly detected in the light-mineral fraction, with specific gravity less than 2.85. Most chromium within the heavy-mineral fraction was substituted within magnetite mineral grains less than 100 micrometers (μm) in diameter, and almost no chromite was present within the heavy-mineral fraction. Although magnetite is resistive to weathering, weathering of magnetite to hematite was identified (1) in Miocene materials underlying unconsolidated deposits in the western subarea of Hinkley Valley and (2) in alluvium within Water Valley that contains weathered minerals eroded from Miocene rock. Less-dense, more easily weathered chromium-containing amphiboles, such as actinolite in older Mojave River alluvium and hornblende in locally derived alluvium from Iron Mountain, were identified optically. Magnetite was not identified in weathered hornblende diorite and was less abundant in locally derived materials and in Miocene materials than in Mojave-type deposits. A comparison of ICP-MS data and sequential extraction data shows that approximately 90 percent of chromium in aquifer material within Hinkley and Water Valleys was not extractable and was interpreted to reside within unweathered mineral grains. Most extractable chromium was within the strong acid extractable fraction. Chromium within the weakly sorbed, and specifically sorbed extractable fractions in oxide accumulations within the regulatory Cr(VI) plume is potentially mobile into groundwater with changes in ionic strength or pH. Although Hinkley and Water Valleys are regionally low in chromium, natural geologic sources of chromium may be present in aquifer materials penetrated by wells completed in (1) weathered hornblende diorite bedrock underlying the western subarea; (2) Miocene deposits underlying the western subarea and unconsolidated material in the northern subarea and Water Valley containing basalt or weathered minerals eroded from Miocene deposits; (3) unconsolidated material containing visually abundant iron- and manganese-oxide coatings on the surfaces of mineral grains that are present near the water table and near lithologic or geologic contacts; and (4) brown clay and mudflat/playa deposits in the northern subarea. Brown clay and mudflat/playa deposits in the eastern subarea near Mount General have a low-chromium, felsic mineralogy similar to Mojave River deposits and do not contain high concentrations of chromium; however, manganese(IV) oxides within these materials may facilitate oxidation of trivalent chromium, Cr(III), to Cr(VI).

California↗

Inland water greenhouse gas budgets for RECCAP2: 1. State-of-the-art of global scale assessments

Inland waters are important sources of the greenhouse gasses (GHGs) carbon dioxide (CO 2 ), methane (CH 4 ) and nitrous oxide (N 2 O) to the atmosphere. In the framework of the 2 nd phase of the REgional Carbon Cycle Assessment and Processes (RECCAP-2) initiative, we review the state of the art in estimating inland water GHG budgets at global scale, which has substantially advanced since the first phase of RECCAP nearly ten years ago. The development of increasingly sophisticated upscaling techniques, including statistical prediction and process based models, allows for spatially explicit estimates which are needed for regionalized assessments of continental GHG budgets such as those established for RECCAP. A few recent estimates also resolve the seasonal and/or interannual variability in inland water GHG emissions. Nonetheless, the global-scale assessment of inland water emissions remains challenging because of limited spatial and temporal coverage of observations and persisting uncertainties in the abundance and distribution of inland water surface areas. To decrease these uncertainties, more empirical work on the contributions of hot-spots and hot-moments to overall inland water GHG emissions is particularly needed.

Global Biogeochemical Cycles↗

Radar image and data fusion for natural hazards characterisation

Fusion of synthetic aperture radar (SAR) images through interferometric, polarimetric and tomographic processing provides an all - weather imaging capability to characterise and monitor various natural hazards. This article outlines interferometric synthetic aperture radar (InSAR) processing and products and their utility for natural hazards characterisation, provides an overview of the techniques and applications related to fusion of SAR/InSAR images with optical and other images and highlights the emerging SAR fusion technologies. In addition to providing precise land - surface digital elevation maps, SAR - derived imaging products can map millimetre - scale elevation changes driven by volcanic, seismic and hydrogeologic processes, by landslides and wildfires and other natural hazards. With products derived from the fusion of SAR and other images, scientists can monitor the progress of flooding, estimate water storage changes in wetlands for improved hydrological modelling predictions and assessments of future flood impacts and map vegetation structure on a global scale and monitor its changes due to such processes as fire, volcanic eruption and deforestation. With the availability of SAR images in near real - time from multiple satellites in the near future, the fusion of SAR images with other images and data is playing an increasingly important role in understanding and forecasting natural hazards.

International Journal of Image and Data Fusion↗

A methodology to assess the historical environmental footprint of in-situ recovery (ISR) of uranium: A demonstration in the Goliad Sand in the Texas Coastal Plain, USA

In-situ recovery (ISR) has been the only technique used to extract uranium from sandstone-hosted uranium deposits in the Pliocene Goliad Sand in the Texas Coastal Plain. Water plays a crucial role throughout the ISR lifecycle of production and groundwater restoration yet neither the water use nor other environmental footprints have been well documented. The goal of this study is to examine historical records for all six ISR operations completed in the Goliad Sand to identify and quantify parameters that indicate the surface and aquifer disturbances, water use, and radon emissions. Overall, the average mine area was 0.00023 ± 0.00006 acres per pound (ac/lb) U 3 O 8 . The average mine pore volume was 48.9 ± 50 gal/lb U 3 O 8 with a minimum affected aquifer volume of 0.51 ± 0.08 cubic feet per pound (cu ft/lb) U 3 O 8. An average of 258 ± 40 gallons (gal) of fluid were disposed per pound (lb) U 3 O 8 , with an average of 169 ± 26 gal/lb U 3 O 8 attributed to restoration and 89 ± 36 gal/lb U 3 O 8 attributed to the uranium production phase. The average radon emitted was 1.06 × 10 −3 ± 7.4 × 10 −4 curies per pound (Ci/lb) U 3 O 8 . Goodness-of-fit ( R 2 ) values are ≥0.79 for linear regressions of the amount of uranium produced versus mine area, mine pore volumes, mine aquifer volumes, water pumped, and total water disposed. The R 2 value for radon emitted was 0.68. However, the water disposed only during the uranium production phase is more strongly correlated to the number of production days ( R 2 = 0.96) than to uranium production ( R 2 = 0.84), whereas the volume of water disposed during restoration is more strongly correlated to the “pore volume” ( R 2 = 0.97) than to uranium production ( R 2 = 0.90). Pore volume is an industry term used to describe the amount of fluid circulated through the aquifer during the uranium production period and stipulated in bond agreements in order to satisfy groundwater restoration requirements. Models constructed in this study can be used to estimate probable water use and the extent of surface and aquifer disturbances associated with ISR-amenable undiscovered uranium resources in the Goliad Sand. The historical perspective offered by the data compiled and correlations may prove useful to both industry and regulators.

Texas↗

Experimental abrasion of detrital gold

The physical breakdown and abrasion rates of gold were studied using a tumbler to simulate natural high-energy environments. The gold fragments were tumbled for periods ranging from 30 to 240 h with different combinations of sand, cobbles, and water at velocities of 0.5 and 2.0 mi/h (0.85 and 3.22 km/h). With sand and gravel, the common bedload of the rivers that deposited the gold-bearing Tertiary sedimentary rocks of the Sierra Nevada, gold is abraded at rates of 0.015 to 0.007 percent (by weight) per hour of travel (at 0.5 mi/h or 0.845 km/h). Cobbles, rather than sand, are responsible for most of the physical changes and abrasion of the gold. Ten gold fragments tumbled for 120 h with cobbles and water (no sand) were broken down to 68 recoverable fragments and lost about 25 percent of their weight to particles smaller than could be recovered using conventional panning techniques. Gold tumbled for 120 h with sand and water lost less than 1 percent of its weight. Gold was abraded faster by wet sand than by dry sand. Velocity appears to be more important as a factor in abrasion of gold than travel distance a fourfold increase in velocity produced a tenfold increase in hourly abrasion rates of gold. Scanning electron microscope examination of the gold fragments after the tumbling experiments revealed differences in surface texture between fragments tumbled with (1) sand, (2) sand and cobbles, and (3) cobbles only.

Journal of Research of the U.S. Geological Survey↗

Water-balance techniques for determining available soil-water storage for selected sandy and clay soil study sites in Cass County, North Dakota, 2016–17

The U.S. Geological Survey, in cooperation with the U.S. Department of Agriculture Natural Resources Conservation Service, collected field and remotely sensed data on precipitation, evapotranspiration (ET), and soil-water content to determine available soil-water storage (AWS) at six study sites on sandy and clay soils in Cass County, North Dakota. Data were collected at all the study sites from May 1–October 31, 2016, and from May 1–October 24, 2017. Estimated daily AWS was determined using daily meteorological and potential evapotranspiration (PET) data obtained from various climate stations, and estimated monthly AWS was determined using monthly meteorological and PET data and monthly ET data determined using the Operational Simplified Surface Energy Balance model. AWS during 2016 and 2017 was determined at daily and monthly time steps because of data availability and to assess results using varying time steps. Comparisons of measured and estimated daily values of AWS at the Brewer Lake site indicated poor agreement during May–October 2016 and May–October 2017. Comparisons of measured and estimated daily values of AWS at the Embden East and Embden West sites indicated poor and fair agreement respectively. At the Lynchburg Crop and Lynchburg Grass sites, comparisons of measured and estimated daily values of AWS indicated fair and good relations, respectively, even with the possible effects of soil cracks. Mean estimated values of daily runoff plus soil percolation for the four sandy soil sites indicated that a maximum of about 19 percent of the estimated runoff plus soil percolation could be considered runoff and that the remaining 81 percent could be considered soil percolation, and for the two clay soil sites about 13 percent of the runoff plus soil percolation could have been considered runoff and about 87 percent could have been considered soil percolation. Results indicated little difference between using monthly PET or monthly ET in water-balance equations to estimate monthly AWS for the grouped sandy soil sites, and only slightly better results were obtained using monthly PET than monthly ET to estimate monthly AWS for the grouped clay soil study sites. Overall, the monthly water-balance models did not perform as well as the daily water-balance models for determining AWS at the six study sites. Additional data collection from a longer-period study and adjustments to the models may improve results from the monthly water-balance techniques.

North Dakota↗

Determination of the delta(15N/14N)of Ammonium (NH4+) in Water: RSIL Lab Code 2898

The purpose of the technique described by Reston Stable Isotope Laboratory (RSIL) lab code 2898 is to determine the N isotopic composition, delta(15N/14N), abbreviated as d15N, of ammonium (NH4+) in water (freshwater and saline water). The procedure involves converting dissolved NH4+ into NH3 gas by raising the pH of the sample to above 9 with MgO and subsequently trapping the gas quantitatively as (NH4)2SO4 on a glass fiber (GF) filter. The GF filter is saturated with NaHSO4 and pressure sealed between two gas-permeable polypropylene filters. The GF filter 'sandwich' floats on the surface of the water sample in a closed bottle. NH3 diffuses from the water through the polypropylene filter and reacts with NaHSO4, forming (NH4)2SO4 on the GF filter. The GF filter containing (NH4)2SO4 is dried and then combusted with a Carlo Erba NC 2500 elemental analyzer (EA), which is used to convert total nitrogen in a solid sample into N2 gas. The EA is connected to a continuous-flow isotope-ratio mass spectrometer (CF-IRMS), which determines the relative difference in ratios of the amounts of the stable isotopes of nitrogen (15N and 14N) of the product N2 gas and a reference N2 gas. The filters containing the samples are compressed in tin capsules and loaded into a Costech Zero-Blank Autosampler on the EA. Under computer control, samples then are dropped into a heated reaction tube that contains an oxidant, where combustion takes place in a He atmosphere containing an excess of O2 gas. To remove S-O gases produced from the NaHSO4, a plug of Ag-coated Cu wool is inserted at the bottom of the reaction tube. Combustion products are transported by a He carrier through a reduction furnace to remove excess O2, toconvert all nitrogen oxides to N2, and to remove any remaining S-O gases. The gases then pass through a drying tube to remove water. The gas-phase products, mainly N2 and a small amount of background CO2, are separated by a gas chromatograph (GC). The gas is then introduced into the IRMS through a Finnigan ConFlo II interface. The ConFlo II interface is used to introduce not only sample into the IRMS but also N2 reference gas and He for sample dilution. The flash combustion is quantitative, so no isotopic fractionation is involved. The IRMS is a Finnigan Delta V CF-IRMS with 10 cups and is capable of detecting ion beams with mass/charge (m/z) 28, 29, 30. The ion beams from N2 are as follows: m/z 28 = 14N14N, m/z 29 = 14N15N, and m/z 30 = 15N15N. The ion beam with m/z 30 also represents 14N16O, which may indicate contamination or incomplete reduction.

Techniques and Methods↗

Groundwater-flow model of the Treasure Valley, southwestern Idaho, 1986–2015

Most of the population of the Treasure Valley and the surrounding area of southwestern Idaho and easternmost Oregon depends on groundwater for domestic supply, either from domestic or municipal-supply wells. Current and projected rapid population growth in the area has caused concern about the long-term sustainability of the groundwater resource. In 2016, the U.S. Geological Survey, in cooperation with the Idaho Water Resource Board and the Idaho Department of Water Resources, began a project to construct a numerical groundwater-flow model of the westernmost portion of the western Snake River Plain aquifer system, called the Treasure Valley. The development of the model was guided by several objectives, including: to improve the understanding of groundwater and surface water interactions; to facilitate conjunctive water management; to provide a tool for water resources planning; and to provide a tool for water allocation. The model was constructed with a spatial scale and level of detail that aimed to meet these objectives while balancing the sometimes-competing goals of fast runtimes, numerical stability, usability, and parsimony. The Treasure Valley Groundwater Flow Model (TVGWFM) is a three-dimensional finite-difference numerical model constructed using MODFLOW 6 (Langevin and others, 2017, Documentation for the MODFLOW 6 Groundwater Flow Model: U.S. Geological Survey Techniques and Methods, book 6, chap. A55, 197 p., https://doi.org/10.3133/tm6A55 ). The model covers the westernmost portion of the western Snake River Plain and is discretized into a regular grid of 64 by 65 cells with a side length of 1 mile and 6 layers of varying depth and active area. A historical model period was developed consisting of 360 month-long stress periods for 1986–2015. The model builds upon previous modeling efforts by adding a transient period, incorporating new head and discharge observations to constrain parameters, incorporating information from the hydrogeologic framework model (HFM) of Bartolino (2019, Hydrogeologic framework of the Treasure Valley and surrounding area, Idaho and Oregon: U.S. Geological Survey Scientific Investigations Report 2019–5138, https://doi.org/10.3133/sir20195138 ) and incorporating refined estimates of evapotranspiration and irrigation classification of lands in the study area. The TVGWFM includes all significant components of recharge to and discharge from the aquifer. Inflows include canal seepage, irrigation and precipitation recharge, mountain-front recharge, rivers and stream seepage, and seepage from Lake Lowell. Outflows include discharge to agricultural drainage ditches, discharge to rivers and streams, pumping, and discharge to Lake Lowell. Each recharge or discharge component is represented separately using individual MODFLOW 6 packages. Parameter values were derived with a combination of trial-and-error steps and automated parameter estimation using PEST software (Doherty, J.E., 2005, PEST, model-independent parameter estimation–User manual: Watermark Numerical Computing, https://pesthomepage.org/documentation ). Parameter estimates were constrained with several types of observation data, including water levels, water level changes, vertical water level differences, drain discharges, change in drain discharges, river seepage, seepage from Lake Lowell, and change in seepage from Lake Lowell. Material properties from the hydrogeologic framework were also used to assign the minimum and maximum values of some parameters. A final parameter realization was reached that minimized residuals between the observed and modelled values for the various observation groups. Mean residuals for the observation groups were 15.4 feet (ft) for water levels, 0.2 ft for water level changes, 19.4 ft for vertical water level differences, −3.9 cubic feet per second (ft 3 /s) for drain discharges, 0.0 ft 3 /s for changes in drain discharge, 45.0 ft 3 /s for river seepage, −40.1 ft 3 /s for Lake Lowell seepage, and 126.3 ft 3 /s for changes in Lake Lowell seepage. The quality of the model’s fit to observations varied spatially, with notable areas of under- or over-simulation of water levels present to the northwest and southwest of Lake Lowell, in the foothills along the eastern model boundary, and near the City of Eagle. Trends were observed in the residuals of many of the observation groups, indicating that the model is missing or not fully reproducing some phenomena that are observed in the system. The TVGWFM can be used as a tool for water resource planning, for understanding the interactions of groundwater and surface water at a basin scale, and for facilitating conjunctive management, but may lack the precision needed for water rights administration at a local scale. Additional sources of uncertainty or limitations of the model are noted. The quantity and spatial distribution of canal seepage and infiltration of irrigation water recharge, the largest sources of recharge to the system, are unknown and approximated indirectly. There is poor understanding of how canal seepage and incidental recharge change as land is converted from agricultural (irrigated) to suburban (semi-irrigated). These uncertainties will affect any scenarios that investigate changes to land use or irrigation practices. Finally, the model has relatively high water-level residuals around and to the southwest of Lake Lowell and should not be used to estimate water level effects in that region. The model was built with multiple, broadly expressed objectives and did not optimize performance for specific uses. However, the model and the tools included in an associated data release provide ample flexibility to improve the model for future uses. Adjustments and improvements could be made by refining the model in an area of interest, collecting additional calibration data, applying more rigorous boundary conditions, or re-estimating model parameters to optimize model performance for a specific model forecast.

Idaho↗

Revision of the documentation for a model for calculating effects of liquid waste disposal in deep saline aquifers

The model developed under this contract is a modified version of the deep well disposal model developed by INTERCOMP Resource Development and Engineering, Inc., for the U.S. Geological Survey (A model for calculating effects of liquid waste disposal in deep saline aquifers). The model is a finite-difference numerical solution of the partial differential equations describing (1) single phase fluid flow in aquifer, (2) energy transport by convection and conduction, and (3) contaminant transport dissolved in the fluid by convection and dispersion. Both the energy and the contaminant transports include molecular diffusion and hydrodynamic dispersion. The objective of the contract was to modify the original version of the model for more general applications. Some of the major features are as follows: (1) Fluid density is permitted to be a function of fluid pressure, temperature and contaminant concentration, the viscosity can be described as a function of pressure and temperature. (2) Aquifer heterogeneities in the hydrological properties can be described on a numerical grid block basis. (3) Free water surface is permitted to exist in the aquifer for shallow ground water applications. (4) Contaminant may absorb on rock surface or decay according to a first order reaction. The absorption is described by a linear adsorption isotherm function of rock type in the aquifer. (5) The model is extremely flexible in providing a wide choice of boundary conditions. These include natural flow in the aquifer, aquifer influence functions around the perimeter of the grid in recognition that the gridded region does not have no-flow boundaries, heat losses into the overlying impermeable strata, and the wellbore heat and pressure drop calculations coupled to the aquifer flow equations. (6) The model offers the option of selecting as iterative or direct solution technique, and selecting central or backward finite-difference approximations in both time and space. (7) Virtually any aquifer can be modeled by proper grid block description in three dimensions. In addition, the model is fully transient. The major limitation of the model arises using second-order correct (central-difference) finite-difference approximation in space. To avoid numerical oscillations in the solution, the user must restrict grid block and time step sizes depending upon the magnitude of the dispersivity.

Water-Resources Investigations Report↗

Analyzing spatial distributions and alignments of pitted cone features in Utopia Planitia on Mars

Martian geomorphology and surface features provide links to understanding past geologic processes such as fluid movement, local and regional tectonics, and feature formation mechanisms. Pitted cones are common features in the northern plains basins of Mars. They have been proposed to have formed from upwelling volatile-rich fluids, such as magma or water-sediment slurries. In this study, we map the spatial distributions of pitted cone features across the Utopia Planitia (UP) basin. Using the Average Nearest Neighbor technique, we find that pitted cone features appear to be clustered across the basin, occurring in a narrow band around the circumferential rim of UP. Additionally, we find that pitted cone features also appear to be aligned in chains that are sub-parallel with the basin margin. Parallel bands of cones generally follow elevation contours of the UP basin, which suggests elevation or a correlated factor plays a major role in pitted cone and cone chain formation. We propose that pitted cones and cone chains may be related to vertical fractures formed around the UP basin rim from subsidence of infilling basin material.

Icarus↗

Creep of water ices at planetary conditions: A compilation

Many constitutive laws for the flow of ice have been published since the advent of the Voyager explorations of the outer solar system. Conflicting data have occasionally come from different laboratories, and refinement of experimental techniques has led to the publication of laws that supersede earlier ones. In addition, there are unpublished data from ongoing research that also amend the constitutive laws. Here we compile the most current laboratory-derived flow laws for water ice phases I, II, III, V, and VI, and ice I mixtures with hard particulates. The rheology of interest is mainly that of steady state, and the conditions reviewed are the pressures and temperatures applicable to the surfaces and interiors of icy moons of the outer solar system. Advances in grain-size-dependent creep in ices I and II as well as in phase transformations and metastability under differential stress are also included in this compilation. At laboratory strain rates the several ice polymorphs are rheologically distinct in terms of their stress, temperature, and pressure dependencies but, with the exception of ice III, have fairly similar strengths. Hard particulates strengthen ice I significantly only at high particulate volume fractions. Ice III has the potential for significantly affecting mantle dynamics because it is much weaker than the other polymorphs and its region of stability, which may extend metastably well into what is nominally the ice II field, is located near likely geotherms of large icy moons.

Journal of Geophysical Research E: Planets↗