Search USGS⌕ Search

SEARCH · Search USGS

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 811 records · Page 45Linked to original sources

Regional variations in water quality and relationships to soil and bedrock weathering in the southern Sacramento Valley, California, USA

Regional patterns in ground- and surface-water chemistry of the southern Sacramento Valley in California were evaluated using publicly available geochemical data from the US Geological Survey's National Water Information System (NWIS). Within the boundaries of the study area, more than 2300 ground-water analyses and more than 20,000 surface-water analyses were available. Ground-waters from the west side of the Sacramento Valley contain greater concentrations of Na, Ca, Mg, B, Cl and SO 4 , while the east-side ground-waters contain greater concentrations of silica and K. These differences result from variations in surface-water chemistry as well as from chemical reactions between water and aquifer materials. Sediments that fill the Sacramento Valley were derived from highlands to the west (the Coast Ranges) and east (the Sierra Nevada Mountains), the former having an oceanic provenance and the latter continental. These geologic differences are at least in part responsible for the observed patterns in ground-water chemistry. Thermal springs that are common along the west side of the Sacramento Valley appear to have an effect on surface-water chemistry, which in turn may affect the ground-water chemistry.

Applied Geochemistry↗

Indicators of biodiversity and ecosystem services: A synthesis across ecosystems and spatial scales

According to the Millennium Ecosystem Assessment, common indicators are needed to monitor the loss of biodiversity and the implications for the sustainable provision of ecosystem services. However, a variety of indicators are already being used resulting in many, mostly incompatible, monitoring systems. In order to synthesise the different indicator approaches and to detect gaps in the development of common indicator systems, we examined 531 indicators that have been reported in 617 peer-reviewed journal articles between 1997 and 2007. Special emphasis was placed on comparing indicators of biodiversity and ecosystem services across ecosystems (forests, grass- and shrublands, wetlands, rivers, lakes, soils and agro-ecosystems) and spatial scales (from patch to global scale). The application of biological indicators was found most often focused on regional and finer spatial scales with few indicators applied across ecosystem types. Abiotic indicators, such as physico-chemical parameters and measures of area and fragmentation, are most frequently used at broader (regional to continental) scales. Despite its multiple dimensions, biodiversity is usually equated with species richness only. The functional, structural and genetic components of biodiversity are poorly addressed despite their potential value across habitats and scales. Ecosystem service indicators are mostly used to estimate regulating and supporting services but generally differ between ecosystem types as they reflect ecosystem-specific services. Despite great effort to develop indicator systems over the past decade, there is still a considerable gap in the widespread use of indicators for many of the multiple components of biodiversity and ecosystem services, and a need to develop common monitoring schemes within and across habitats. Filling these gaps is a prerequisite for linking biodiversity dynamics with ecosystem service delivery and to achieving the goals of global and sub-global initiatives to halt the loss of biodiversity. ?? 2009 Oikos.

Oikos↗

Characterizing near-surface CO2 conditions before injection - Perspectives from a CCS project in the Illinois Basin, USA

The Midwest Geological Sequestration Consortium is conducting a large-scale carbon capture and storage (CCS) project in Decatur, Illinois, USA to demonstrate the ability of a deep saline formation to store one million tonnes of carbon dioxide (CO 2 ) from an ethanol facility. Beginning in early 2011, CO 2 will be injected at a rate of 1,000 tonnes/day for three years into the Mount Simon Sandstone at a depth of approximately 2,100 meters. An extensive Monitoring, Verification, and Accounting (MVA) program has been undertaken for the Illinois Basin Decatur Project (IBDP) and is focused on the 0.65 km2 project site. Goals include establishing baseline conditions to evaluate potential impacts from CO 2 injection, demonstrating that project activities are protective of human health and the environment, and providing an accurate accounting of stored CO 2 . MVA efforts are being conducted pre-, during, and post- CO 2 injection. Soil and net CO 2 flux monitoring has been conducted for more than one year to characterize near-surface CO 2 conditions. More than 2,200 soil CO 2 flux measurements have been manually collected from a network of 118 soil rings since June 2009. Three ring types have been evaluated to determine which type may be the most effective in detecting potential CO 2 leakage. Bare soil, shallow-depth rings were driven 8 cm into the ground and were prepared to minimize surface vegetation in and near the rings. Bare soil, deep-depth rings were prepared similarly, but were driven 46 cm. Natural-vegetation, shallow-depth rings were driven 8 cm and are most representative of typical vegetation conditions. Bare-soil, shallow-depth rings had the smallest observed mean flux (1.78 μmol m −2 s −1 ) versus natural-vegetation, shallow-depth rings (3.38 μmol m −2 s −1 ). Current data suggest bare ring types would be more sensitive to small CO 2 leak signatures than natural ring types because of higher signal to noise ratios. An eddy covariance (EC) system has been in use since June 2009. Baseline data from EC monitoring is being used to characterize pre-injection conditions, and may then be used to detect changes in net exchange CO 2 fluxes (Fc) that could be the result of CO 2 leakage into the near-surface environment during or following injection. When injection at IBDP begins, soil and net CO 2 monitoring efforts will have established a baseline of near-surface conditions that will be important to help demonstrate the effectiveness of storage activities.

Illinois↗

Numerical simulation of ground-water flow in the central part of the western San Joaquin Valley, California

The occurrence of selenium in agricultural drain water in the central part of the western San Joaquin Valley, California, has focused concern on strategies for managing shallow, saline ground water. To assess alternatives to agricultural drains, a three-dimensional, finite-difference numerical model of the regional groundwater flow system was developed. This report documents the mathematical approach used to model the flow system, the data base on which the model is based, and the methods used to calibrate the model. The 550-square-mile study area includes parts of the Panoche Creek alluvial fan and parts of the Little Panoche Creek and Cantua Creek alluvial fans. The model simulates transient flow in the semiconfined and confined zones above and below the Corcoran Clay Member of the Tulare Formation of Pleistocene age. The model incorporates areally distributed ground-water recharge, areally and vertically distributed pumping, regional-collector drains in the Wesdands Water District (operative from 1980 to 1985), on-farm drains in parts of the Panoche, Broadview, and Firebaugh Water Districts, and bare-soil evaporation (which occurs if the water table is within 7 feet of land surface). The model also incorporates texture-based estimates of hydraulic conductivity, where texture is defined as the fraction of coarse-grained deposits present in a given subsurface interval. The numerical model was developed using hydrologic data from 1972 to 1988. Most of the parameters incorporated into the model were evaluated independently of the model, including system geometry, the distribution of texture, the altitudes of the water table and potentiometric surface of the confined zone in 1972 (initial condition), the hydraulic conductivity of coarse-grained deposits derived from the Coast Ranges, the hydraulic conductivity of coarse-grained deposits derived from the Sierra Nevada, specific storage, recharge, pumping, and parameters needed to incorporate drains and bare-soil evaporation. Four parameters were calibration variables: the hydraulic conductivity of fine-grained deposits in the semiconfined zone, the hydraulic conductivity of the Corcoran Clay Member, specific yield, and the transmissivity of the confined zone. The model was calibrated in two phases. In the first phase, a steady-state model of the ground-water flow system in 1984 was used to constrain the relation between the hydraulic conductivity of fine-grained deposits in the semiconfined zone and the hydraulic conductivity of the Corcoran Clay Member, thus reducing the number of independent variables from four to three. In the second phase of calibration, the change in altitude of the water table from 1972 to 1984, the change in altitude of the potentiometric surface of the confined zone from 1972 to 1984, and the number of model cells subject to bare-soil evaporation from 1972 to 1988 were used to evaluate the remaining three variables. The calibrated model reproduces the average change in water-table altitude (1972-84) to within 0.4 foot (average measured change 11.5 feet) and the average change in confined zone head (1972- 84) to within 19 feet (average measured change 120 feet). The simulated time-series record of the total number of model cells subject to bare-soil evaporation (each cell is 1 mile square) is within the range of the measured data. The measured values are at a minimum in October and a maximum in July. The October values ranged from 103 in 1972 to 132 in 1984 (the drains were closed in 1985) to 151 in 1988. The July values ranged from 144 in 1973 to 198 in 1984, to 204 in 1988. The simulated values ranged from 103 in 1972 to 161 in 1984, to 208 in 1988.

Water Supply Paper↗

Regional evaluation of evapotranspiration in the Everglades

Nine sites in the Florida Everglades were selected and instrumented for collection of data necessary for evapotranspiration-determination using the Bowen-ratio energy-budget method. The sites were selected to represent the sawgrass or cattail marshes, wet prairie, and open-water areas that constitute most of the natural Everglades system. At each site, measurements necessary for evapotranspiration (ET) calculation and modeling were automatically made and stored on-site at 15- or 30-minute intervals. Data collected included air temperature and humidity at two heights, wind speed and direction, incoming solar radiation, net solar radiation, water level and temperature, soil moisture content, soil temperature, soil heat flux, and rainfall. Data summarized in this report were collected from January 1996 through December 1997, and the development of site-specific and regional models of ET for this period is described. Latent heat flux is the energy flux density equivalent of the ET rate. Modified Priestley-Taylor models of latent heat flux as a function of selected independent variables were developed at each site. These models were used to fill in periods of missing latent heat flux measurement, and to develop regional models of the entire Everglades region. The regional models may be used to estimate ET in wet prairie, sawgrass or cattail marsh, and open-water portions of the natural Everglades system. The models are not applicable to forested areas or to the brackish areas adjacent to Florida Bay. Two types of regional models were developed. One type of model uses measurements of available energy at a site, together with incoming solar energy and water depth, to estimate hourly ET. This available-energy model requires site data for net radiation, water heat storage, and soil heat flux, as well as data for incoming solar radiation and water depth. The other type of model requires only incoming solar energy, air temperature, and water depth data to provide estimates of hourly ET. The second model thus uses data that are more readily available than the data required for the available-energy model. Computed ET mean annual totals for all nine sites for the 1996-97 period ranged from 42.4 inches per year at a site where the water level is below land surface for several months each year to 57.4 inches per year at an open-water site with no emergent vegetation. Although the density of photosynthetically-active plant leaves has been shown to relate directly to ET in some studies, it does not appear to relate directly to ET in the Everglades, based on comparison of annual ET data with leaf-area index, defined as the Normalized Difference Vegetation Index (NDVI), data from satellite imagery. NDVI and ET appear to be inversely related in the Everglades. The greatest ET rates occurred at open-water sites where the NDVI data indicated the lowest leaf-area index. Among the remaining vegetated sites, there is no clear relation between ET and NDVI, though the highest ET rate corresponded to the lowest NDVI and one of the lowest ET rates corresponded to the highest NDVI value. The variation in ET follows a seasonal pattern, with lowest monthly ET totals occurring in December through February, and highest ET occurring in May through August. The greatest range in monthly ET among all nine sites for the 2-year period occurred at site 3: from 1.81 inches in December 1997 to 6.84 inches in July 1996. A study to compare the Bowen-ratio/energy balance method of ET measurement with the eddy-correlation method was done at one site from June 22, 1998, through September 28, 1998. This comparison indicated that both methods gave comparable values of the Bowen ratio, but there was a considerable difference in available energy measured by the two methods. The mean of all 30-minute measured turbulent heat fluxes from the eddy-correlation apparatus for June 22 through September 29, 1998, was 137.4 watts per square meter, and the mean of the corresponding measured energy was 163.6 watts per square meter, or about 20 percent greater. The disagreement in mean energy fluxes measured by the two methods is problematical and is not fully understood. Although the difference seems to be related to friction velocity, and is practically non-existent at values of friction velocity greater than 0.3 meter per second, the "correctness" of either method cannot be determined with the data available.

Florida↗

Managing vegetation in surface-flow wastewater-treatment wetlands for optimal treatment performance

Constructed wetlands that mimic natural marshes have been used as low-cost alternatives to conventional secondary or tertiary wastewater treatment in the U.S. for at least 30 years. However, the general level of understanding of internal treatment processes and their relation to vegetation and habitat quality has not grown in proportion to the popularity of these systems. We have studied internal processes in surface-flow constructed wastewater-treatment wetlands throughout the southwestern U.S. since 1990. At any given time, the water quality, hydraulics, water temperature, soil chemistry, available oxygen, microbial communities, macroinvertebrates, and vegetation each greatly affect the treatment capabilities of the wetland. Inside the wetland, each of these components plays a functional role and the treatment outcome depends upon how the various components interact. Vegetation plays a uniquely important role in water treatment due to the large number of functions it supports, particularly with regard to nitrogen transformations. However, it has been our experience that vegetation management is critical for achieving and sustaining optimal treatment function. Effective water treatment function and good wildlife quality within a surface-flow constructed wetland depend upon the health and sustainability of the vegetation. We suggest that an effective tool to manage and sustain healthy vegetation is the use of hummocks, which are shallow emergent plant beds within the wetland, positioned perpendicular to the water flow path and surrounded by water sufficiently deep to limit further emergent vegetation expansion. In this paper, we describe the use of a hummock configuration, in conjunction with seasonal water level fluctuations, to manage the vegetation and maintain the treatment function of wastewater-treatment wetlands on a sustainable basis.

Ecological Engineering↗

Quantification of natural vapor fluxes of trichloroethene in the unsaturated zone at Picatinny Arsenal, New Jersey

The upward flux of trichloroethene (TCE) vapor through the unsaturated zone above a contaminated, water-table aquifer at Picatinny Arsenal, New Jersey, has been studied under natural conditions over a 12-month period. Vertical gas-phase diffusion fluxes were estimated indirectly by measuring the TCE vapor concentration gradient in the unsaturated zone and using Fick's law to calculate the flux. The total gas-phase flux (e.g., the sum of diffusion and advection fluxes) was measured directly with a vertical flux chamber (VFC). In many cases, the upward TCE vapor flux was several orders of magnitude greater than the upward TCE diffusion flux, suggesting that mechanisms other than steady-state vapor diffusion are contributing to the vertical transport of TCE vapors through the unsaturated zone. The measured total flux of TCE vapor from the subsurface to the atmosphere is approximately 50 kg/yr and is comparable in magnitude to the removal rate of TCE from the aquifer by an existing pump-and-treat system and by discharge into a nearby stream. The net upward flux of TCE is reduced significantly during a storm event, presumably due to the mass transfer of TCE from the soil gas to the infiltrating rainwater and its subsequent downward advection. Several potential problems associated with the measurement of total gas-phase fluxes are discussed.

Environmental Science & Technology↗

Experimental warming and its legacy effects on root dynamics following two hurricane disturbances in a wet tropical forest

Tropical forests are expected to experience unprecedented warming and increases in hurricane disturbances in the coming decades; yet, our understanding of how these productive systems, especially their belowground component, will respond to the combined effects of varied environmental changes remains empirically limited. Here we evaluated the responses of root dynamics (production, mortality, and biomass) to soil and understory warming (+4°C) and after two consecutive tropical hurricanes in our in situ warming experiment in a tropical forest of Puerto Rico: Tropical Responses to Altered Climate Experiment (TRACE). We collected minirhizotron images from three warmed plots and three control plots of 12 m 2 . Following Hurricanes Irma and María in September 2017, the infrared heater warming treatment was suspended for repairs, which allowed us to explore potential legacy effects of prior warming on forest recovery. We found that warming significantly reduced root production and root biomass over time. Following hurricane disturbance, both root biomass and production increased substantially across all plots; the root biomass increased 2.8-fold in controls but only 1.6-fold in previously warmed plots. This pattern held true for both herbaceous and woody roots, suggesting that the consistent antecedent warming conditions reduced root capacity to recover following hurricane disturbance. Root production and mortality were both related to soil ammonium nitrogen and microbial biomass nitrogen before and after the hurricanes. This experiment has provided an unprecedented look at the complex interactive effects of disturbance and climate change on the root component of a tropical forested ecosystem. A decrease in root production in a warmer world and slower root recovery after a major hurricane disturbance, as observed here, are likely to have longer-term consequences for tropical forest responses to future global change.

Puerto Rico↗

CO 2 outgassing in a combined fracture and conduit karst aquifer near lititz spring, Pennsylvania

Lititz Spring in southeastern Pennsylvania and a nearby domestic well were sampled for 9 months. Although both locations are connected to conduits (as evidenced by a tracer test), most of the year they were saturated with respect to calcite, which is more typical of matrix flow. Geochemical modeling (PHREEQC) was used to explain this apparent paradox and to infer changes in matrix and conduit contribution to flow. The saturation index varied from 0.5 to 0 most of the year, with a few samples in springtime dropping below saturation. The log PCO 2 value varied from -2.5 to -1.7. Lower log PCO 2 values (closer to the atmospheric value of -3.5) were observed when the solutions were at or above saturation with respect to calcite. In contrast, samples collected in the springtime had high PCO 2 , low saturation indices, and high water levels. Geochemical modeling showed that when outgassing occurs from a water with initially high PCO 2 , the saturation index of calcite increases. In the Lititz Spring area, the recharge water travels through the soil zone, where it picks up CO 2 from soil gas, and excess CO 2 subsequently is outgassed when this recharge water reaches the conduit. At times of high water level (pipe full), recharge with excess CO 2 enters the system but the outgassing does not occur. Instead the recharge causes dilution, reducing the calcite saturation index. Understanding the temporal and spatial variation in matrix and conduit flow in karst aquifers benefited here by geochemical modeling and calculation of PCO 2 values. ?? 2006 Geological Society of America.

Special Paper of the Geological Society of America↗

Watershed characterization for precipitation-runoff modeling system, north fork, American River and east fork, Carson River watersheds, California

As part of its Global Change Hydrology Program, the U.S. Geological Survey (USGS) is investigating the potential effects of climate change on the water resources of several river basins in the United States. The American River Basin in California represents the windward slope of the north-central Sierra Nevada, and the California part of the Carson River Basin, most of which is in Nevada, represents the leeward slope. Parts of the American River and Carson River Basins—the North Fork American River and East Fork Carson River watersheds, both in California—were studied to determine the sensitivity of water resources to potential climate change. The water resources of both basins are derived primarily from snowmelt. A geographic information system (GIS) data base has been created to facilitate paired-basin analysis. The GIS data base incorporates (1) land-surface data, which include elevation, land use and land cover, soil type, and geology; (2) hydrologic data, such as stream networks and streamflow-gaging stations; and (3) climatic data, such as snow-course, snow-telemetry, radiosonde, and meteorological data. Precipitation-runoff models were developed and calibrated for the North Fork watershed within the American River Basin and for the East Fork watershed within the Carson River Basin. (These watersheds were selected to represent the climatic and physiographic variability of the two larger basins.) Synthesized climate scenarios then were used in the model to predict potential effects of climate change.

Hydrologic Atlas↗

Distributed routing rainfall-runoff model; version II

A computer program of a watershed model for routing storm runoff through a branched system of pipes and (or) natural channels using rainfall as input is described. The model provides detailed simulation of storm-runoff periods selected by the user and a daily soil-moisture accounting between storms. A drainage basin is represented as a set of overland-flow, channel, and reservoir seqments which ,jointly describe the drainage features of the basin. Kinematic wave theory is used for routing flows over contributing overland-flow areas and through the channel network. A set of model segments can be arranged into a network that will represent many complex drainage basins. The model is intended primarily for application to urban watersheds, but may have limited applications to rural watersheds.

Open-File Report↗

Water Budget of East Maui, Hawaii

Ground-water recharge is estimated from six monthly water budgets calculated using long-term average rainfall and streamflow data, estimated pan-evaporation and fog-drip data, and soil characteristics. The water-budget components are defined seasonally, through the use of monthly data, and spatially by broad climatic and geohydrologic areas, through the use of a geographic information system model. The long-term average water budget for east Maui was estimated for natural land-use conditions. The average rainfall, fog-drip, runoff, evapotranspiration, and ground-water recharge volumes for the east Maui study area are 2,246 Mgal/d, 323 Mgal/d, 771 Mgal/d, 735 Mgal/d, and 1,064 Mgal/d, respectively.

Water-Resources Investigations Report↗

Areas contributing recharge to production wells and effects of climate change on the groundwater system in the Chipuxet River and Chickasheen Brook Basins, Rhode Island

The Chipuxet River and Chickasheen Brook Basins in southern Rhode Island are an important water resource for public and domestic supply, irrigation, recreation, and aquatic habitat. The U.S. Geological Survey, in cooperation with the Rhode Island Department of Health, began a study in 2012 as part of an effort to protect the source of water to six large-capacity production wells that supply drinking water and to increase understanding of how climate change might affect the water resources in the basins. Soil-water-balance and groundwater-flow models were developed to delineate the areas contributing recharge to the wells and to quantify the hydrologic response to climate change. Surficial deposits of glacial origin ranging from a few feet to more than 200 feet thick overlie bedrock in the 24.4-square mile study area. These deposits comprise a complex and productive aquifer system. Simulated areas contributing recharge to the production wells covered a total area of 0.63 square miles for average well withdrawal rates from 2007 through 2011 (total rate of 583 gallons per minute). Simulated areas contributing recharge for the maximum well pumping capacities (total rate of 3,700 gallons per minute) covered a total area of 2.55 square miles. Most simulated areas contributing recharge extend upgradient of the wells to morainal and upland till deposits and to groundwater divides. Some simulated areas contributing recharge include small, isolated areas remote from the wells. Relatively short groundwater traveltimes from recharging locations to discharging wells indicated that the wells are vulnerable to contamination from land-surface activities; median traveltimes ranged from 3.5 to 8.6 years for the production wells examined, and 57 to 91 percent of the traveltimes were 10 years or less. Land cover in the areas contributing recharge includes a substantial amount of urban and agriculture land use for five wells adjacent to the Chipuxet River; for one well adjacent to a tributary stream, land use is less developed. The calibrated groundwater-flow model provided a single, best representation of the areas contributing recharge to a production well. The parameter variance-covariance matrix from model calibration was used to create parameter sets that reflect the uncertainty of the parameter estimates and the correlation among parameters to evaluate the uncertainty associated with the predicted contributing areas to the wells. A Monte Carlo analysis led to contributing areas expressed as a probability distribution that differed from a single deterministic contributing area. Because of the effects of parameter uncertainty, the size of the probabilistic contributing areas for both average and maximum pumping rates was larger than the size of the deterministic contributing areas for the wells. Thus, some areas not in the deterministic contributing area might actually be in the contributing area, including additional areas of urban and agricultural land use that has the potential to contaminate groundwater. Additional areas that might be in the contributing area included recharge originating near the pumping wells that have relatively short groundwater-flow paths and traveltimes. At the maximum pumping rates, areas associated with low probabilities extended long distances along groundwater divides in the uplands remote from the wells. Climate projections for the Chipuxet River and Chickasheen Brook Basins from downscaled output from general circulation models indicate that mean annual temperature might increase by 4.7 degrees Fahrenheit and 8.0 degrees Fahrenheit by the late 21st century (2070–99) compared with the late 20th century (1970–99) under scenarios of lower and higher emissions of greenhouse gases, respectively. By the late 21st century, winter and spring precipitation is projected to increase by 12 to 17 percent, summer precipitation to increase by about the same as mean annual precipitation (8 percent), and fall precipitation to decrease by 5 percent for both emission scenarios compared with the late 20th century. Soil-water-balance simulations indicate that, although precipitation is expected to increase in three seasons, only in winter do precipitation increases exceed actual evapotranspiration increases. Recharge is projected to decrease in fall and generally change little in spring and summer. By the late 21st century, winter recharge is expected to increase by 13 percent for the lower emissions scenario and by 15 percent for the higher emissions scenario. In fall, recharge is projected to diminish by 13 percent for the lower emissions scenario and by 24 percent for the higher emissions scenario. Although recharge is projected to change seasonally in the 21st century, mean annual recharge changes minimally. Soil moisture is projected to decrease in the 21st century from spring through fall because of increases in potential evapotranspiration, and in fall because of decreases in precipitation in addition to increases in potential evapotranspiration. By the late 21st century, soil moisture for the lower emissions scenario is expected to decrease by 11 percent in summer and 15 percent in fall, and for the higher emissions scenario, decrease by 23 percent for both seasons. These decreases in soil moisture during the growing season might have implications for agriculture in the study area. Predicted changes in the magnitude and seasonal distribution of recharge in the 21st century increase simulated base flows and groundwater levels in the winter months for both emission scenarios, but because of less recharge in the fall and less or about the same recharge in the preceding months of spring and summer, base flows and groundwater levels in the fall months decrease for both emission scenarios. October has the largest base flow and groundwater level decreases. By the late 21st century, base flows at the Chipuxet River in October are projected to decrease by 9 percent for the lower emissions scenario and 18 percent for the higher emissions scenario. For a headwater stream in the upland till with shorter groundwater-flow paths and lower storage properties in its drainage area, base flows in October are projected to diminish by 28 percent and 42 percent for the lower and higher emissions scenarios by the late 21st century. Groundwater level changes in the uplands show substantial decreases in fall, but because of the large storage capacity of stratified deposits, water levels change minimally in the valley. By the late 21st century, water levels in large areas of upland till deposits in October are projected to decrease by up to 2 feet for the lower emissions scenario, whereas large areas decrease by up to 5 feet, with small areas with decreases of as much as 10 feet, for the higher emissions scenario. For both emission scenarios, additional areas of till go dry in fall compared with the late 20th century. Thus projected changes in recharge in the 21st century might extend low flows and low water levels for the year later in fall and there might be more intermittent headwater streams compared with the late 20th century with corresponding implications to aquatic habitat. Finally, the size and location of the simulated areas contributing recharge to the production wells are minimally affected by climate change because mean annual recharge, which is used to determine the contributing areas to the production wells, is projected to change little in the 21st century.

Rhode Island↗

A prototype system for forecasting landslides in the Seattle, Washington, area

Empirical rainfall thresholds and related information form the basis of a prototype system for forecasting landslides in the Seattle area. The forecasts are tied to four alert levels, and a decision tree guides the use of thresholds to determine the appropriate level. From analysis of historical landslide data, we developed a formula for a cumulative rainfall threshold (CT), P 3 = 88.9 − 0.67 P 15 , defined by rainfall amounts in millimeters during consecutive 3 d (72 h) periods, P 3 , and the 15 d (360 h) period before P 3 , P 15 . The variable CT captures more than 90% of historical events of three or more landslides in 1 d and 3 d periods recorded from 1978 to 2003. However, the low probability of landslide occurrence on a day when the CT is exceeded at one or more rain gauges (8.4%) justifies a low-level of alert for possible landslide occurrence, but it does trigger more vigilant monitoring of rainfall and soil wetness. Exceedance of a rainfall intensity-duration threshold I = 82.73 D −1.13 , for intensity, I (mm/hr), and duration, D (hr), corresponds to a higher probability of landslide occurrence (30%) and forms the basis for issuing warnings of impending, widespread occurrence of landslides. Information about the area of exceedance and soil wetness can be used to increase the certainty of landslide forecasts (probabilities as great as 71%). Automated analysis of real-time rainfall and subsurface water data and digital quantitative precipitation forecasts are needed to fully implement a warning system based on the two thresholds.

Washington↗

Preliminary maps of Quaternary deposits and liquefaction susceptibility, nine-county San Francisco Bay region, California: a digital database

This report presents a preliminary map and database of Quaternary deposits and liquefaction susceptibility for the nine-county San Francisco Bay region, together with a digital compendium of ground effects associated with past earthquakes in the region. The report consists of (1) a spatial database of fivedata layers (Quaternary deposits, quadrangle index, and three ground effects layers) and two text layers (a labels and leaders layer for Quaternary deposits and for ground effects), (2) two small-scale colored maps (Quaternary deposits and liquefaction susceptibility), (3) a text describing the Quaternary map, liquefaction interpretation, and the ground effects compendium, and (4) the databse description pamphlet. The nine counties surrounding San Francisco Bay straddle the San Andreas fault system, which exposes the region to serious earthquake hazard (Working Group on California Earthquake Probabilities, 1999). Much of the land adjacent to the Bay and the major rivers and streams is underlain by unconsolidated deposits that are particularly vulnerable to earthquake shaking and liquefaction of water-saturated granular sediment. This new map provides a modern and regionally consistent treatment of Quaternary surficial deposits that builds on the pioneering mapping of Helley and Lajoie (Helley and others, 1979) and such intervening work as Atwater (1982), Helley and others (1994), and Helley and Graymer (1997a and b). Like these earlier studies, the current mapping uses geomorphic expression, pedogenic soils, and inferred depositional environments to define and distinguish the map units. In contrast to the twelve map units of Helley and Lajoie, however, this new map uses a complex stratigraphy of some forty units, which permits a more realistic portrayal of the Quaternary depositional system. The two colored maps provide a regional summary of the new mapping at a scale of 1:275,000, a scale that is sufficient to show the general distribution and relationships of the map units but cannot distinguish the more detailed elements that are present in the database. The report is the product of years of cooperative work by the USGS National Earthquake Hazards Reduction Program (NEHRP) and National Cooperative Geologic Mapping Program, William Lettis and & Associates, Inc. (WLA) and, more recently, by the California Division of Mines and Geology as well. An earlier version was submitted to the Geological Survey by WLA as a final report for a NEHRP grant (Knudsen and others, 2000). The mapping has been carried out by WLA geologists under contract to the NEHRP Earthquake Program (Grants #14-08-0001-G2129, 1434-94-G-2499, 1434-HQ-97-GR-03121, and 99-HQ-GR-0095) and with other limited support from the County of Napa, and recently also by the California Division of Mines and Geology. The current map consists of this new mapping and revisions of previous USGS mapping.

Open-File Report↗

Use of output from high‐resolution atmospheric models in landscape‐scale hydrologic models: An assessment

In this paper we investigate the feasibility of coupling regional climate models (RCMs) with landscape‐scale hydrologic models (LSHMs) for studies of the effects of climate on hydrologic systems. The RCM used is the National Center for Atmospheric Research/Pennsylvania State University mesoscale model (MM4). Output from two year‐round simulations (1983 and 1988) over the western United States is used to drive a lake model for Pyramid Lake in Nevada and a streamfiow model for Steamboat Creek in Oregon. Comparisons with observed data indicate that MM4 is able to produce meteorologic data sets that can be used to drive hydrologic models. Results from the lake model simulations indicate that the use of MM4 output produces reasonably good predictions of surface temperature and evaporation. Results from the streamflow simulations indicate that the use of MM4 output results in good simulations of the seasonal cycle of streamflow, but deficiencies in simulated wintertime precipitation resulted in underestimates of streamflow and soil moisture. Further work with climate (multiyear) simulations is necessary to achieve a complete analysis, but the results from this study indicate that coupling of LSHMs and RCMs may be a useful approach for evaluating the effects of climate change on hydrologic systems.

Water Resources Research↗

Sediment transport and capacity change in three reservoirs, Lower Susquehanna River Basin, Pennsylvania and Maryland, 1900-2012

The U.S. Geological Survey (USGS) has conducted numerous sediment transport studies in the Susquehanna River and in particular in three reservoirs in the Lower Susquehanna River Basin to determine sediment transport rates over the past century and to document changes in storage capacity. The Susquehanna River is the largest tributary to Chesapeake Bay and transports about one-half of the total freshwater input and substantial amounts of sediment and nutrients to the bay. The transported loads are affected by deposition in reservoirs (Lake Clarke, Lake Aldred, and Conowingo Reservoir) behind three hydropower dams. The geometry and texture of the deposited sediments in each reservoir upstream from the three dams has been a subject of research in recent decades. Particle size deposition and sediment scouring processes are part of the reservoir dynamics. A Total Maximum Daily Load (TMDL) for nitrogen, phosphorus, and sediment was established for Chesapeake Bay to attain water-quality standards. Six states and the District of Columbia agreed to reduce loads to the bay and to meet load allocation goals for the TMDL. The USGS has been estimating annual sediment loads at the Susquehanna River at Marietta, Pennsylvania (above Lake Clarke), and Susquehanna River at Conowingo, Maryland (below Conowingo Reservoir), since the mid-1980s to predict the mass balance of sediment transport through the reservoir system. Using streamflow and sediment data from the Susquehanna River at Harrisburg, Pennsylvania (upstream from the reservoirs), from 1900 to 1981, sediment loads were greatest in the early to mid-1900s when land disturbance activities from coal production and agriculture were at their peak. Sediment loads declined in the 1950s with the introduction of agricultural soil conservation practices. Loads were dominated by climatic factors in the 1960s (drought) and 1970s (very wet) and have been declining since the 1980s through 2012. The USGS developed a regression equation to predict the sediment scour load for daily mean streamflows greater than 300,000 cubic feet per second for the Lower Susquehanna River reservoirs. A compilation of data from various sources produced a range in total sediment transported through the reservoir system and allowed for apportioning to source (watershed or scour) for various streamflows. In 2011, Conowingo Reservoir was estimated to be about 92 percent of sediment storage capacity. Since construction of Conowingo Dam in 1929 through 2012, approximately 470 million tons of sediment was transported down the Susquehanna River into the reservoir system, approximately 290 million tons were trapped, and approximately 180 million tons were transported to Chesapeake Bay. Spatial and estimated total sand deposition in Conowingo Reservoir based on historical sediment cores indicated continued migration of sand downgradient toward the dam and the winnowing of silts and clays near the dam due to scour.

Maryland, Pennsylvania↗

Hydrogeologic characterization of Area B, Fort Detrick, Maryland

Groundwater in the karst groundwater system at Area B of Fort Detrick in Frederick County, Maryland, is contaminated with chlorinated solvents from the past disposal of laboratory wastes. In cooperation with U.S. Army Environmental Command and U.S. Army Garrison Fort Detrick, the U.S. Geological Survey performed a 3-year study to refine the conceptual model of groundwater flow in and around Area B of Fort Detrick at the site- to regional-scale. The investigation was designed to review the geologic setting, assess the temporal variability of the hydrologic system, evaluate the potential for interbasin groundwater flow, determine the degree of vertical connectivity of the aquifer, characterize the sources and timing of groundwater recharge, and identify if dyes from previous tracer tests continue to drain from the aquifer. This study established a continuous hydrologic monitoring network of 12 water level gages, 2 streamgages, a precipitation gage, and in situ fluorometric monitoring. A water budget analysis was performed using hydrologic monitoring data and a soil-water balance model constructed for the study. In this study each individual water budget term is calculated using available data or through modeling, and a water budget residual term is calculated. If the water budget residual term is small relative to the uncertainty of the underlying data, then an additional import or export of water (in other words, interbasin transfer) is not needed to fully describe the hydrologic system. Groundwater and spring samples from 20 locations were collected in a 2019 synoptic geochemical sampling event and analyzed for a suite of analytes that included groundwater age tracer constituents. The karst groundwater system was found to be highly responsive to hydrologic events, with strong water level and stream base flow responses to individual storm events and a historic wet period in 2017 and 2018. The water budget analysis included historic flooding in May 2018, though more typical hydrologic patterns were observed in 2019 and 2020. During most evaluated intervals, the water budget residual was less than the estimated uncertainty on the residual for the two Carroll Creek watersheds, which suggested no substantial net interbasin flow occurs from these watersheds. The watershed difference area, a region that includes Area B, had a significant negative water budget residual, which may be the result of a net interbasin import of groundwater or the result of focused groundwater recharge not simulated by the soil-water balance model. Geochemical analysis and groundwater age dating reveals shallow groundwater (approximately less than [<] 150 feet deep) appears to be relatively young (approximately <30 years) and to be recharged in the vicinity of Area B. In the deep groundwater sampled in this study (approximately greater than [>] 150 feet deep), older groundwater from a differing recharge source, based on stable isotopes and noble gas analyses, is observed and interpreted to represent less direct connectivity to the surface and increased proportions of water recharged to the north and (or) west of Area B. A clustering analysis to reveal groupings within the suite of geochemical data was used to define seven groups. The groupings generally show that wells in similar depths and lateral aquifer positions generally cluster together, with some exceptions. Although limited by suspended sediments, the in situ fluorometric monitoring at springs did not detect any dye leaving the system above the limit of detection for the method. Dye was only detected above the limit of detection in one well, which was used as an injection well during a previous dye tracer test. The results of this study support and refine the conceptual site model of groundwater hydrology at Area B. The geologic and geophysical log review in this study agrees with prior assessments of physical controls on groundwater flow. A literature review of mid-Atlantic karst studies identified similar controls reported in these environments. The additional characterization of hydrologic responsiveness in this study suggests that hydrologic conditions and events are important considerations when interpreting potentiometric surfaces and contaminant trends over time and highlights the importance of continuous hydrologic monitoring. There is evidence to suggest that either intense focused groundwater recharge occurs in the vicinity of Area B or net along-valley groundwater interbasin flow from the upper study watershed enters the lower watershed and discharges to Carroll Creek. Geochemical analyses also suggest that water recharged from Catoctin Mountain and the elevated areas to the north and (or) west of the site may be present in the older and deeper Area B groundwater.

Maryland↗