USGS Science⌕ Search

SEARCH · USGS Science

Results for “Soil Systems”

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 739 records · Page 41Linked to original sources

Conceptual and numerical models of the glacial aquifer system north of Aberdeen, South Dakota

This U.S. Geological Survey report documents a conceptual and numerical model of the glacial aquifer system north of Aberdeen, South Dakota, that can be used to evaluate and manage the city of Aberdeen's water resources. The glacial aquifer system in the model area includes the Elm, Middle James, and Deep James aquifers, with intervening confining units composed of glacial till. The Elm aquifer ranged in thickness from less than 1 to about 95 feet (ft), with an average thickness of about 24 ft; the Middle James aquifer ranged in thickness from less than 1 to 91 ft, with an average thickness of 13 ft; and the Deep James aquifer ranged in thickness from less than 1 to 165 ft, with an average thickness of 23 ft. The confining units between the aquifers consisted of glacial till and ranged in thickness from 0 to 280 ft. The general direction of groundwater flow in the Elm aquifer in the model area was from northwest to southeast following the topography. Groundwater flow in the Middle James aquifer was to the southeast. Sparse data indicated a fairly flat potentiometric surface for the Deep James aquifer. Horizontal hydraulic conductivity for the Elm aquifer determined from aquifer tests ranged from 97 to 418 feet per day (ft/d), and a confined storage coefficient was determined to be 2.4x10 -5 . Estimates of the vertical hydraulic conductivity of the sediments separating the Elm River from the Elm aquifer, determined from the analysis of temperature gradients, ranged from 0.14 to 2.48 ft/d. Average annual precipitation in the model area was 19.6 inches per year (in/yr), and agriculture was the primary land use. Recharge to the Elm aquifer was by infiltration of precipitation through overlying outwash, lake sediments, and glacial till. The annual recharge for the model area, calculated by using a soil-water-balance method for water year (WY) 1975-2009, ranged from 0.028 inch in WY 1980 to 4.52 inches in WY 1986, with a mean of 1.56 inches. The annual potential evapotranspiration, calculated in soil-water-balance analysis, ranged from 21.8 inches in WY 1983 to 27.0 inches in WY 1985, with a mean of 24.6 inches. Water use from the glacial aquifer system primarily was from the Elm aquifer for irrigation, municipal, and suburban water supplies, and the annual rate ranged from 1.0 to 2.4 cubic feet per second (ft 3 /s). The MODFLOW-2005 numerical model represented the Elm aquifer, the Middle James aquifer, and the Deep James aquifer with model layers 1-3 respectively separated by confining layers 1-2 respectively. Groundwater flow was simulated with 75 stress periods beginning October 1, 1974, and ending September 30, 2009. Model grid spacing was 200 by 200 ft and boundaries were represented by specified-head boundaries and no-flow boundaries. The model used parameter estimation that focused on minimizing the difference between 954 observed and simulated hydraulic heads for 135 wells. Calibrated mean horizontal hydraulic conductivity values for model layers 1-3 were 94, 41, and 30 ft/d respectively. Vertical hydraulic conductivity values for confining layers 1 and 2 were 0.0002 and 0.0003 ft/d, respectively. Calibrated specific yield for model layer 1was 0.1 and specific storage ranged from 0.0003 to 0.0005 per foot. Calibrated mean recharge rates ranged from 2.5 in/yr where glacial till thickness was less than 10 ft to 0.8 in/yr where glacial till thickness was greater than 30 ft. Calibrated mean annual evapotranspiration rate was 8.8 in/yr. Simulated net streamflow gain from model layer 1 was 3.1 ft 3 /s.

South Dakota↗

Response of anaerobic carbon cycling to water table manipulation in an Alaskan rich fen

To test the effects of altered hydrology on organic soil decomposition, we investigated CO 2 and CH 4 production potential of rich-fen peat (mean surface pH = 6.3) collected from a field water table manipulation experiment including control, raised and lowered water table treatments. Mean anaerobic CO 2 production potential at 10 cm depth (14.1 ± 0.9 μmol C g −1 d −1 ) was as high as aerobic CO 2 production potential (10.6 ± 1.5 μmol C g −1 d −1 ), while CH4 production was low (mean of 7.8 ± 1.5 nmol C g −1 d −1 ). Denitrification enzyme activity indicated a very high denitrification potential (197 ± 23 μg N g −1 d −1 ), but net NO - 3 reduction suggested this was a relatively minor pathway for anaerobic CO 2 production. Abundances of denitrifier genes ( nirK and nosZ ) did not change across water table treatments. SO 2- 4 reduction also did not appear to be an important pathway for anaerobic CO 2 production. The net accumulation of acetate and formate as decomposition end products in the raised water table treatment suggested that fermentation was a significant pathway for carbon mineralization, even in the presence of NO - 3 . Dissolved organic carbon (DOC) concentrations were the strongest predictors of potential anaerobic and aerobic CO 2 production. Across all water table treatments, the CO 2 :CH 4 ratio increased with initial DOC leachate concentrations. While the field water table treatment did not have a significant effect on mean CO 2 or CH 4 production potential, the CO 2 :CH 4 ratio was highest in shallow peat incubations from the drained treatment. These data suggest that with continued drying or with a more variable water table, anaerobic CO 2 production may be favored over CH 4 production in this rich fen. Future research examining the potential for dissolved organic substances to facilitate anaerobic respiration, or alternative redox processes that limit the effectiveness of organic acids as substrates in anaerobic metabolism, would help explain additional uncertainty concerning carbon mineralization in this system.

Soil Biology and Biochemistry↗

Dust in the Critical Zone: North American case studies

The dust cycle facilitates the exchange of particles among Earth's major systems, enabling dust to traverse ecosystems, cross geographic boundaries, and even move uphill against the natural flow of gravity. Dust in the atmosphere is composed of a complex and ever-changing mixture that reflects the evolving human footprint on the landscape. The emission, transport, and deposition of dust interacts with and connects Critical Zone processes at all spatial and temporal scales. Landscape properties, land use, and climatic factors influence the wind erosion of soil and nutrient loss, which alters the long-term ecological dynamics at erosional locations. Once in the atmosphere, dust particles influence the amount of solar radiation reaching Earth, and interact with longwave (terrestrial) radiation, with cascading effects on the climate system. Finally, the wet and dry deposition of particles influences ecosystem structure, composition, and function over both short and long-term scales. Tracking dust particles from source to sink relies on monitoring and measurement of the geochemical composition and size distribution of the particles, space-borne and ground-based remote sensing, and dust modeling. Dust is linked to human systems via land use and policies that contribute to dust emissions and the health-related consequences of particulate loads and composition. Despite the significant influence dust has in shaping coupled natural-human systems, it has not been considered a key component of the Critical Zone. Here, we demonstrate that dust particles should be included as a key component of the Critical Zone by outlining how dust interacts with and shapes Earth System processes from generation, through transport, to deposition. We synthesize current understanding from global research and identify critical data and knowledge gaps while showcasing case studies from North America.

Earth-Science Reviews↗

Effects of tillage and application rate on atrazine transport to subsurface drainage: Evaluation of RZWQM using a six-year field study

Well tested agricultural system models can improve our understanding of the water quality effects of management practices under different conditions. The Root Zone Water Quality Model (RZWQM) has been tested under a variety of conditions. However, the current model's ability to simulate pesticide transport to subsurface drain flow over a long term period under different tillage systems and application rates is not clear. Therefore, we calibrated and tested RZWQM using six years of data from Nashua, Iowa. In this experiment, atrazine was spring applied at 2.8 (1990–1992) and 0.6 kg/ha/yr (1993–1995) to two 0.4 ha plots with different tillage (till and no-till). The observed and simulated average annual flow weighted atrazine concentrations (FWAC) in subsurface drain flow from the no-till plot were 3.7 and 3.2 μg/L, respectively for the period with high atrazine application rates, and 0.8 and 0.9 μg/L, respectively for the period with low application rates. The 1990–1992 observed average annual FWAC difference between the no-till and tilled plot was 2.4 μg/L while the simulated difference was 2.1 μg/L. These observed and simulated differences for 1993–1995 were 0.1 and 0.1 μg/L, respectively. The Nash–Sutcliffe model performance statistic (EF) for cumulative atrazine flux to subsurface drain flow was 0.93 for the no-till plot testing years (1993–1995), which is comparable to other recent model tests. The value of EF is 1.0 when simulated data perfectly match observed data. The order of selected parameter sensitivity for RZWQM simulated FWAC was atrazine partition coefficient > number of macropores > atrazine half life in soil > soil hydraulic conductivity. Simulations from 1990 to 1995 with four different atrazine application rates applied at a constant rate throughout the simulation period showed concentrations in drain flow for the no-till plot to be twice those of the tilled plot. The differences were more pronounced in the early simulation period (1990–1992), partly because of the characteristics of macropore flow during large storms. The results suggest that RZWQM is a promising tool to study pesticide transport to subsurface drain flow under different tillage systems and application rates over several years, the concentrations of atrazine in drain flow can be higher with no-till than tilled soil over a range of atrazine application rates, and atrazine concentrations in drain flow are sensitive to the macropore flow characteristics under different tillage systems and rainfall timing and intensity.

Agricultural Water Management↗

Ecosystem responses to elevated CO2 using airborne remote sensing at Mammoth Mountain, California

We present an exploratory study examining the use of airborne remote-sensing observations to detect ecological responses to elevated CO 2 emissions from active volcanic systems. To evaluate these ecosystem responses, existing spectroscopic, thermal, and lidar data acquired over forest ecosystems on Mammoth Mountain volcano, California, were exploited, along with in situ measurements of persistent volcanic soil CO 2 fluxes. The elevated CO 2 response was used to statistically model ecosystem structure, composition, and function, evaluated via data products including biomass, plant foliar traits and vegetation indices, and evapotranspiration (ET). Using regression ensemble models, we found that soil CO 2 flux was a significant predictor for ecological variables, including canopy greenness (normalized vegetation difference index, NDVI), canopy nitrogen, ET, and biomass. With increasing CO 2 , we found a decrease in ET and an increase in canopy nitrogen, both consistent with theory, suggesting more water- and nutrient-use-efficient canopies. However, we also observed a decrease in NDVI with increasing CO 2 (a mean NDVI of 0.27 at 200 g m −2 d −1 CO 2 reduced to a mean NDVI of 0.10 at 800 g m −2 d −1 CO 2 ). This is inconsistent with theory though consistent with increased efficiency of fewer leaves. We found a decrease in above-ground biomass with increasing CO 2 , also inconsistent with theory, but we did also find a decrease in biomass variance, pointing to a long-term homogenization of structure with elevated CO 2 . Additionally, the relationships between ecological variables changed with elevated CO 2 , suggesting a shift in coupling/decoupling among ecosystem structure, composition, and function synergies. For example, ET and biomass were significantly correlated for areas without elevated CO 2 flux but decoupled with elevated CO 2 flux. This study demonstrates that (a) volcanic systems show great potential as a means to study the properties of ecosystems and their responses to elevated CO 2 emissions and (b) these ecosystem responses are measurable using a suite of airborne remotely sensed data.

California↗

Delineation and evaluation of hydrologic-landscape regions in the United States using geographic information system tools and multivariate statistical analyses.

Hydrologic-landscape regions in the United States were delineated by using geographic information system (GIS) tools combined with principal components and cluster analyses. The GIS and statistical analyses were applied to land-surface form, geologic texture (permeability of the soil and bedrock), and climate variables that describe the physical and climatic setting of 43,931 small (approximately 200 km 2 ) watersheds in the United States. (The term “watersheds” is defined in this paper as the drainage areas of tributary streams, headwater streams, and stream segments lying between two confluences.) The analyses grouped the watersheds into 20 noncontiguous regions based on similarities in land-surface form, geologic texture, and climate characteristics. The percentage of explained variance (R-squared value) in an analysis of variance was used to compare the hydrologic-landscape regions to 19 square geometric regions and the 21 U.S. Environmental Protection Agency level-II ecoregions. Hydrologic-landscape regions generally were better than ecoregions at delineating regions of distinct land-surface form and geologic texture. Hydrologic-landscape regions and ecoregions were equally effective at defining regions in terms of climate, land cover, and water-quality characteristics. For about half of the landscape, climate, and water-quality characteristics, the R-squared values of square geometric regions were as high as hydrologic-landscape regions or ecoregions.

Environmental Management↗

Synopsis of strontium isotope variations in groundwater at Aspo, southern Sweden

Strontium isotope ratios are used to identify end-member ground-water compositions at Aspo in southeastern Sweden where the Hard Rock Laboratory (HRL) has been constructed to evaluate the suitability of crystalline rock for the geologic disposal of nuclear waste. The Hard Rock Laboratory is a decline (tunnel) constructed in 1.8 Ga-old granitic rock that forms islands in an archipelago along the Swedish coast. Ground-water samples were obtained for isotopic analyses from boreholes drilled from the surface and from side boreholes drilled within the HRL. Infiltration at Aspo occurs primarily through fractures zones in the granitic bedrock beneath thin soils throughout the area. Because of extremely low Sr concentrations, rain and snow are not important contributors to the Sr isotope budget of the ground-water system. At shallow levels, water percolating downward along fractures and fracture zones acquires a ??87Sr between +9.5 and +10.0??? and maintains this value downward while Sr concentrations increase by two orders of magnitude. Ground-water samples from both boreholes and from in the HRL show the effects of mixing with saline waters containing as much as 59 mg/L Sr and ??87Sr values as large as +13.92%, Baltic Sea water is a potential component of the groundwater system with ??87Sr values only slightly larger than modern marine values (+0.3???) but with much lower concentrations (1.5 mg/L) than ocean water (8 mg/L). However, because of large Sr concentration differences between the saline groundwater (59 mg/L) and Baltic Sea water (1.5 rag/L), ??87Sr values are not particularly sensitive indicators of sea-water intrusion even though their ??87Sr values differ substantially.

Applied Geochemistry↗

Rocky 7 prototype Mars rover field geology experiments: 1. Lavic Lake and Sunshine volcanic field, California

Experiments with the Rocky 7 rover were performed in the Mojave Desert to better understand how to conduct rover-based, long-distance (kilometers) geological traverses on Mars. The rover was equipped with stereo imaging systems for remote sensing science and hazard avoidance and 57 Fe Mössbauer and nuclear magnetic resonance spectrometers for in situ determination of mineralogy of unprepared rock and soil surfaces. Laboratory data were also obtained using the spectrometers and an X ray diffraction (XRD)/XRF instrument for unprepared samples collected from the rover sites. Simulated orbital and descent image data assembled for the test sites were found to be critical for assessing the geologic setting, formulating hypotheses to be tested with rover observations, planning traverses, locating the rover, and providing a regional context for interpretation of rover-based observations. Analyses of remote sensing and in situ observations acquired by the rover confirmed inferences made from orbital and simulated descent images that the Sunshine Volcanic Field is composed of basalt flows. Rover data confirmed the idea that Lavic Lake is a recharge playa and that an alluvial fan composed of sediments with felsic compositions has prograded onto the playa. Rover-based discoveries include the inference that the basalt flows are mantled with aeolian sediment and covered with a dense pavement of varnished basalt cobbles. Results demonstrate that the combination of rover remote sensing and in situ analytical observations will significantly increase our understanding of Mars and provide key connecting links between orbital and descent data and analyses of returned samples.

Journal of Geophysical Research E: Planets↗

Catchment disturbance and stream metabolism: Patterns in ecosystem respiration and gross primary production along a gradient of upland soil and vegetation disturbance

Catchment characteristics determine the inputs of sediments and nutrients to streams. As a result, natural or anthropogenic disturbance of upland soil and vegetation can affect instream processes. The Fort Benning Military Installation (near Columbus, Georgia) exhibits a wide range of upland disturbance levels because of spatial variability in the intensity of military training. This gradient of disturbance was used to investigate the effect of upland soil and vegetation disturbance on rates of stream metabolism (ecosystem respiration rate [ER] and gross primary production rate [GPP]). Stream metabolism was measured using an open-system, single-station approach. All streams were net heterotrophic during all seasons. ER was highest in winter and spring and lowest in summer and autumn. ER was negatively correlated with catchment disturbance level in winter, spring, and summer, but not in autumn. ER was positively correlated with abundance of coarse woody debris, but not significantly related to % benthic organic matter. GPP was low in all streams and generally not significantly correlated with disturbance level. Our results suggest that the generally intact riparian zones of these streams were not sufficient to protect them from the effect of upland disturbance, and they emphasize the role of the entire catchment in determining stream structure and function. ?? 2005 by The North American Benthological Society.

Journal of the North American Benthological Societ↗

Evaluating mass flow meter measurements from chambers for greenhouse gas emissions from orphan wells and other point sources

This study evaluates the performance of a rigid gas flux chamber equipped with a mass flow meter (MFM) for measuring gas emissions from leaking orphan wells and similar pressure-driven gas point sources. We conducted a series of laboratory and field experiments to evaluate the sensitivity, stability, and dynamic range of an MFM chamber system and found an optimal method for sealing the chamber to the ground to isolate the emission source. From these results, we estimate the effects of different soil gas permeabilities on measurements and identify the uncertainty of environmental processes that can impact measurements. Simulations of an MFM chamber are compared to those of a dynamic flux chamber to contrast the data derived with both methodologies and illustrate the potential for measuring high variability leaks with the MFM chamber. Using a low flow resistance MFM and a chamber well-sealed to the ground, it is possible to measure leaks down to 1.08 x 10 -3 cubic meters per hour (m 3 h −1 ) (refenced to 25°/1 atm), corresponding to 0.77 grams per hour (g h −1 ) methane or 2.11 g h −1 carbon dioxide, with a mean uncertainty of 0.89 % relative standard deviation. Environmental processes such as heated gas inside the chamber from solar gain, wind blowing across the chamber vent, and changing humidity in the chamber, can cause variation in MFM measurements. Over 11 d of continuous monitoring under varying weather conditions, the standard deviation of the environmentally sourced signals was found to be 7.40 x 10 -3 m 3 h −1 (equivalent to or 5.27 g h −1 methane or 14.45 g h −1 carbon dioxide). Strategies to obtain the highest quality data from MFM chambers include burying the edges of the chamber below the surface sufficiently deep to seal the chamber edges against gas flow and soaking the dirt with water to lower the chances of escaping gases, while monitoring the gas flow and adjusting the chamber seal to achieve a maximum flow rate.

Atmospheric Measurement Techniques↗

Fire in Mediterranean climate ecosystems-A comparative overview

Four regions of the world share a similar climate and structurally similar plant communities with the Mediterranean Basin. These five areas, known collectively as "mediterranean-type climate (MTC) regions", are dominated by evergreen sclerophyllous-leaved shrublands, semi-deciduous scrub, and woodlands, all of which are prone to widespread crown fires. Summer droughts produce an annual fire hazard that contributes to a highly predictable fire regime. Fire has been an important factor driving the convergence of these systems and is reflected in plant traits such as lignotubers in resprouting shrubs and delayed reproduction that restricts recruitment to a postfire pulse of seedlings. On fertile soils where postfire resprouting is very rapid, opportunities for postfire seedling recruitment are limited and thus these woody taxa have not opted for delaying reproduction. Such fire-independent recruitment is widespread in the floras of MTC regions of the Mediterranean Basin and California and postfire seeding tends to dominate at the more arid end of the gradient. Due to very different geological histories in South Africa and Western Australia, substrates are nutrient poor and thus postfire resprouters do not pose a similar competitive challenge to seedlings and thus postfire seeding is very widespread in these floras. Although circumstantial evidence suggests that the MTC region of Chile had fire-prone landscapes in the Tertiary, these were lost with the late Miocene completion of the Andean uplift, which now blocks summer lightning storms from moving into the region. Today these five regions pose a significant fire management challenge due to the annual fire hazard and metropolitan centers juxtaposed with highly flammable vegetation. This challenge varies across the five MTC landscapes as a function of differences in regional fuel loads and population density.

Israel Journal of Ecology and Evolution↗

Muted responses to chronic experimental nitrogen deposition on the Colorado Plateau

Anthropogenic nitrogen (N) deposition is significantly altering both community structure and ecosystem processes in terrestrial ecosystems across the globe. However, our understanding of the consequences of N deposition in dryland systems remains relatively poor, despite evidence that drylands may be particularly vulnerable to increasing N inputs. In this study, we investigated the influence of 7 years of multiple levels of simulated N deposition (0, 2, 5, and 8 kg N ha −1 year −1 ) on plant community structure and biological soil crust (biocrust) cover at three semi-arid grassland sites spanning a soil texture gradient. Biocrusts are a surface community of mosses, lichens, cyanobacteria, and/or algae, and have been shown to be sensitive to N inputs. We hypothesized that N additions would decrease plant diversity, increase abundance of the invasive annual grass Bromus tectorum, and decrease biocrust cover. Contrary to our expectations, we found that N additions did not affect plant diversity or B. tectorum abundance. In partial support of our hypotheses, N additions negatively affected biocrust cover in some years, perhaps driven in part by inter-annual differences in precipitation. Soil inorganic N concentrations showed rapid but ephemeral responses to N additions and plant foliar N concentrations showed no response, indicating that the magnitude of plant and biocrust responses to N fertilization may be buffered by endogenous N cycling. More work is needed to determine N critical load thresholds for plant community and biocrust dynamics in semi-arid systems and the factors that determine the fate of N inputs.

Arizona, Colorado, New Mexico, Utah↗

Excessive phosphorus loading contributes to future vulnerability of mangrove ecosystems by reducing net ecosystem exchange of carbon

J.N. “Ding” Darling National Wildlife Refuge (DDNWR) is located on Sanibel Island along the southwestern coast of Florida, USA. There, eutrophication attributed to agricultural discharge along the Caloosahatchee River has affected the area’s aquatic habitat. In anticipation of additional nutrient loading, we experimentally fertilized mangrove forests with nitrogen (+N; NH4) and phosphorus (+P; P2O5) for 3 years, and monitored soil and pneumatophore CO2 fluxes and tree sap flow from two mangrove species. Furthermore, we modeled individual tree and stand water use, from which we developed carbon (C) budgets for +N and + P vs. control simulations based on a novel application of water use efficiency conversion. Many of the measured response variables provided hints of subtle changes in response to +P rather than +N, which were enhanced when scaled. From this, we found that additional P loading is expected to reduce both gross and net primary productivity as well as CO2 uptake via net ecosystem exchange of C, likely pressing the system beyond metabolic capacity and leading to a 48–62% decrease in projected lateral C export. Greater eutrophication will likely compound vulnerabilities to sea-level rise submergence, especially where P concentrations are high and already reducing soil surface elevations.

Florida↗

Controls on selenium distribution and mobilization in an irrigated shallow groundwater system underlain by Mancos Shale, Uncompahgre River Basin, Colorado, USA

Elevated selenium (Se) concentrations in surface water and groundwater have become a concern in areas of the Western United States due to the deleterious effects of Se on aquatic ecosystems. Elevated Se concentrations are most prevalent in irrigated alluvial valleys underlain by Se-bearing marine shales where Se can be leached from geologic materials into the shallow groundwater and surface water systems. This study presents groundwater chemistry and solid-phase geochemical data from the Uncompahgre River Basin in Western Colorado, an irrigated alluvial landscape underlain by Se-rich Cretaceous marine shale. We analyzed Se species, major and trace elements, and stable nitrogen and oxygen isotopes of nitrate in groundwater and aquifer sediments to examine processes governing selenium release and transport in the shallow groundwater system. Groundwater Se concentrations ranged from below detection limit (< 0.5 μg L − 1 ) to 4070 μg L − 1 , and primarily are controlled by high groundwater nitrate concentrations that maintain oxidizing conditions in the aquifer despite low dissolved oxygen concentrations. High nitrate concentrations in non-irrigated soils and nitrate isotopes indicate nitrate is largely derived from natural sources in the Mancos Shale and alluvial material. Thus, in contrast to areas that receive substantial NO 3 inputs through inorganic fertilizer application, Se mitigation efforts that involve limiting NO 3 application might have little impact on groundwater Se concentrations in the study area. Soluble salts are the primary source of Se to the groundwater system in the study area at-present, but they constitute a small percentage of the total Se content of core material. Sequential extraction results indicate insoluble Se is likely composed of reduced Se in recalcitrant organic matter or discrete selenide phases. Oxidation of reduced Se species that constitute the majority of the Se pool in the study area could be a potential source of Se in the future as soluble salts are progressively depleted.

Colorado↗

Volcano-tectonic structures, gravity and helium in geothermal areas of Tuscany and Latium (Vulsini volcanic district), Italy

Since the early 1980s, geological and structural mapping, gravity, and helium soil-gas studies have been performed in the eastern sector of the Vulsini Volcanic District (Roman Magmatic Province) in an attempt to locate potential geothermal reservoirs. This area is characterised by an anomalous geothermal gradient of > 100??C/km, and by widespread hydrothermal mineralization, thermal springs, high gas fluxes, and fossil and current travertine deposits. The results of these surveys indicate the existence of a number of fault systems, with N-S and E-W structures that appear to be superimposed on older NW-SE and NE-SW features. Comparison of the results of the various studies also reveals differences in permeability and potential reservoir structures at depth.Since the early 1980s, geological and structural mapping, gravity, and helium soil-gas studies have been performed in the eastern sector of the Vulsini Volcanic District (Roman Magmatic Province) in an attempt to locate potential geothermal reservoirs. This area is characterised by an anomalous geothermal gradient of > 100??C/km, and by widespread hydrothermal mineralization, thermal springs, high gas fluxes, and fossil and current travertine deposits. The results of these surveys indicate the existence of a number of fault systems, with N-S and E-W structures that appear to be superimposed on older NW-SE and NE-SW features. Comparison of the results of the various studies also reveals differences in permeability and potential reservoir structures at depth.

Geothermics↗

Simulation of groundwater and surface-water resources of the San Antonio Creek Valley watershed, Santa Barbara County, California

In the San Antonio Creek Valley watershed (SACVW), western Santa Barbara County, California, groundwater is the primary source of water for agricultural irrigation, the town of Los Alamos, and supplemental water to Vandenberg Space Force Base (VSFB). Groundwater pumpage has increased since the 1970s as non-irrigated agricultural land has been converted to irrigated land and as local pumping for municipal use has increased. This increase in groundwater use has resulted in declining groundwater levels, adjustments in surface-water flows and species habitats, and changes in water quality. Water managers are addressing the challenges of meeting this increased demand while maintaining sustainable groundwater supplies. To address these challenges, Santa Barbara County Water Agency, Vandenberg Space Force Base (VSFB), and the U.S. Geological Survey (USGS) undertook a cooperative study to characterize the integrated hydrologic system of the SACVW and develop tools to better understand and manage the groundwater system. The objectives of this study were to improve the understanding of the integrated hydrologic system and incorporate the understanding into an integrated groundwater and surface-water flow model that can be used to help manage the water resources in the SACVW. The San Antonio Creek integrated model (SACIM) was developed using the USGS coupled groundwater and surface-water flow model to simulate the hydrologic system of the SACVW and provide annual and average water budgets for 1948–2018 water years. Results from the SACIM indicated that between 1948 and 2018, total groundwater from storage (storage depletion) for the period was 453,300 acre-feet (acre-ft). Agricultural pumpage was the largest discharge and accounted for a total of 1,020,000 acre-ft of groundwater. Increased pumpage since the mid-1980s (of which agricultural pumpage is the primary component) is tied to an increased rate of storage depletion and reduced rates of groundwater evapotranspiration and surface leakage (groundwater discharge to the surface and soil zone). The increased pumpage also reduced subsurface inflow to Barka Slough, resulting in a decline in upward flow through the underlying hydrogeologic units and surface leakage. In addition to quantifying historical changes in the integrated hydrologic system, the SACIM is a tool than can be used by water managers to evaluate the effects of different climatic and hydrologic conditions and management strategies.

California↗

Subsurface investigation for liquefaction analysis and piezometer calibration at Treasure Island Naval Station, California

Between January and March, 1994, a drilling program was conducted at the Treasure Island Naval Station to evaluate the liquefaction resistance of soils that did and did not liquefy during the Loma Prieta, California, earthquake of October 17, 1989. A second goal of this drilling program was to test and calibrate a retrievable piezometer system that is designed to monitor dynamic pore-water pressure during liquefaction. Retrievable Piezometer A retrievable piezometer can be used to replace failed transducers without redrilling, and the external casing can be installed without the piezometer itself. Many external casings can be installed throughout a region and used only when necessary. At two sites the USGS retrievable piezometer was placed at depths between 2.3 and 4.6m. The retrievable piezometer involves augering a hole to the testing depth and emplacing a 33-mm outside diameter pvc pipe with a porous stone. The hole is back filled and sealed with bentonite, the top of the boring is capped with a box flush to the ground. Later, a commercial transducer is connected to a 21-mm outside diameter pvc pipe and lowered down the 33 mm casing and screwed into the bottom porous-stone assembly. A calibrating transducer (the same type and model as in the USGS retrievable piezometer) was installed inside a penetrometer with a 60° conical tip and an external sleeve that protects the porous filter, located immediately behind the tip, during advancement through dry soil. After the instrument was advanced to the proper depth the tip with the porous filter was advanced past the protective sleeve. Pore pressure was elevated separately by dynamic impact and blasting. The first calibration tests were conducted within the U.S. Geotechnical Test Site established at the Treasure Island fire station (building 157) (de Alba and others, 1994). A 590 kg weight (diameter 72 cm) was dropped 0.69 to 1.63 m onto a steel plate (91 cm square, 0.6 cm thick) to elevate pore pressure, each test involved dropping the weight one time. The surficial distance from the energy source to the piezometers ranged from 1 to 3.6 m. At an empty field (bounded by 11th and 13th streets and H and I ave) pore pressure was elevated using the 590 kg weight and no. 8 blasting caps (50 grains, 3 grams) and primer cord. The explosives were placed 2 m from the piezometers at depths of 2.4 to 2.7 m. A USGS explosives expert handled the explosives under the supervision of Navy personnel. Liquefaction Gray and brown, fine to medium grained sand was hydraulically dredged from San Francisco Bay to create Treasure Island (fig. 1). During the Loma Prieta earthquake Treasure Island experience a peak ground acceleration of 0.16 g (Shakal and other, 1989) and portions of the soil beneath Treasure Island liquefied and were vented to the ground surface as sand boils. During November, 1989, a survey of Treasure Island was made to document ground effects such as sand boils, settlement, and ground cracking. During this survey samples of more than 30 sand boils were taken for grain size analysis (Bennett, in press). Although the soil beneath the fire station did not liquefy, surrounding areas did. The surrounding liquefaction may have affected the peak ground acceleration. Fifteen seconds into the acceleration record at the fire station there is a sudden drop in ground acceleration, and 16 seconds into the record there is practically no response (Idriss, 1991). De Alba and others (1994) ascribe the behavior of the acceleration record to the liquefaction of the underlying sand. Besides the generation of sand boils, Treasure Island experienced significant ground settlement and lateral deformation that damaged lifelines for water and gas (Seed and others, 1990). The primary objective of this report is to document the subsurface stratigraphy at the liquefaction and non-liquefaction sites (fig. 2), and to explore the relation between sand boils and subsurface sediment. This documentation adds to the geotechnical data base of liquefaction by clarifying which layers actually liquefy during earthquakes. Another objective is to briefly describe the piezometer calibration test in terms of what was done and where it was done, results of the calibration work will be reported later by the primary investigator, Behnam Hushmand of Hushmand Associates.

California↗

Simulation of groundwater flow in the regional aquifer system on Long Island, New York, for pumping and recharge conditions in 2005–15

A three-dimensional groundwater-flow model was developed for the aquifer system of Long Island, New York, to evaluate (1) responses of the hydrologic system to changes in natural and anthropogenic hydraulic stresses, (2) the subsurface distribution of groundwater age, and (3) the regional-scale distribution of groundwater travel times and the source of water to fresh surface waters and coastal receiving waters. The model also provides the groundwater flow components used to define model boundaries for possible inset models used for local-scale analyses. The three-dimensional, groundwater flow model developed for this investigation uses the numerical code MODFLOW–NWT to represent steady-state conditions for average groundwater pumping and aquifer recharge for 2005–15. The particle-tracking algorithm MODPATH, which simulates advective transport in the aquifer, was used to estimate groundwater age, delineate the areas at the water table that contribute recharge to coastal and freshwater bodies, and estimate total travel times of water from the water table to discharge locations. A three-dimensional, 1-meter (3.3-foot) topobathymetric model was used to determine land-surface altitudes for the island and seabed altitudes for the surrounding coastal waters. The mapped extents and surface altitudes of major geologic units were compiled and used to develop a three-dimensional hydrogeologic framework of the aquifer system, including aquifers and confining units. Lithologic data from deep boreholes and previous aquifer-test results were used to estimate the three-dimensional distribution of hydraulic conductivity in principal aquifers. Natural recharge from precipitation was estimated for 2005–15 using a modified Thornthwaite-Mather methodology as implemented in a soil-water balance model. Components of anthropogenic recharge—wastewater return flow, storm water inflow, and inflow from leaky infrastructure—also were estimated for 2005–15. Groundwater withdrawals for various sources, including public water supply, industrial, remediation, and agricultural, were compiled or estimated for the same period. These data were incorporated into the model development to represent the aquifer system geometry, boundaries, and initial hydraulic properties of the regional aquifers and confining units within the Long Island aquifer system. Average hydraulic conditions—water levels and streamflows—for 2005–15 were estimated using existing data from the U.S. Geological Survey National Water Information System database. Model inputs were adjusted to best match average hydrologic conditions using inverse methods as implemented in the parameter-estimating software PEST. The calibrated model was used to simulate average hydrologic conditions in the aquifer system for 2005–15. About 656 cubic feet per second of water was withdrawn on average annually for 2005–15 for water supply and an average of about 349 cubic feet per second of water recharged the aquifer annually from return flow and leaky infrastructure. Parts of New York City have drawdowns exceeding 25 feet, mostly because of urbanization and associated large decreases in recharge rates. Large areas in the western part of the island have drawdowns exceeding 10 feet, mostly from large groundwater withdrawals and sewering, which largely eliminates wastewater return flow. Water-table altitudes in eastern parts of the island increased by more than 2 feet in some areas as a result of wastewater return flow in unsewered areas and changes in land use. Changes in streamflows show a similar pattern as water-table altitudes. Streamflows generally decrease in western parts of the island where there are large drawdowns and increase in eastern parts of the island where water-table altitudes increase.

New York↗