USGS ScienceSearch

SEARCH · USGS Science

Results for “Advances in Water Resources”

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.

496 records · Page 28Linked to original sources

Ground-water flow in low permeability environments

Certain geologic media are known to have small permeability; subsurface environments composed of these media and lacking well developed secondary permeability have groundwater flow sytems with many distinctive characteristics. Moreover, groundwater flow in these environments appears to influence the evolution of certain hydrologic, geologic, and geochemical systems, may affect the accumulation of pertroleum and ores, and probably has a role in the structural evolution of parts of the crust. Such environments are also important in the context of waste disposal. This review attempts to synthesize the diverse contributions of various disciplines to the problem of flow in low-permeability environments. Problems hindering analysis are enumerated together with suggested approaches to overcoming them. A common thread running through the discussion is the significance of size- and time-scale limitations of the ability to directly observe flow behavior and make measurements of parameters. These limitations have resulted in rather distinct small- and large-scale approaches to the problem. The first part of the review considers experimental investigations of low-permeability flow, including in situ testing; these are generally conducted on temporal and spatial scales which are relatively small compared with those of interest. Results from this work have provided increasingly detailed information about many aspects of the flow but leave certain questions unanswered. Recent advances in laboratory and in situ testing techniques have permitted measurements of permeability and storage properties in progressively “tighter” media and investigation of transient flow under these conditions. However, very large hydraulic gradients are still required for the tests; an observational gap exists for typical in situ gradients. The applicability of Darcy's law in this range is therefore untested, although claims of observed non-Darcian behavior appear flawed. Two important nonhydraulic flow phenomena, osmosis and ultrafiltration, are experimentally well established in prepared clays but have been incompletely investigated, particularly in undisturbed geologic media. Small-scale experimental results form much of the basis for analyses of flow in low-permeability environments which occurs on scales of time and size too large to permit direct observation. Such large-scale flow behavior is the focus of the second part of the review. Extrapolation of small-scale experimental experience becomes an important and sometimes controversial problem in this context. In large flow systems under steady state conditions the regional permeability can sometimes be determined, but systems with transient flow are more difficult to analyze. The complexity of the problem is enhanced by the sensitivity of large-scale flow to the effects of slow geologic processes. One-dimensional studies have begun to elucidate how simple burial or exhumation can generate transient flow conditions by changing the state of stress and temperature and by burial metamorphism. Investigation of the more complex problem of the interaction of geologic processes and flow in two and three dimensions is just beginning. Because these transient flow analyses have largely been based on flow in experimental scale systems or in relatively permeable systems, deformation in response to effective stress changes is generally treated as linearly elastic; however, this treatment creates difficulties for the long periods of interest because viscoelastic deformation is probably significant. Also, large-scale flow simulations in argillaceous environments generally have neglected osmosis and ultrafiltration, in part because extrapolation of laboratory experience with coupled flow to large scales under in situ conditions is controversial. Nevertheless, the effects are potentially quite important because the coupled flow might cause ultra long lived transient conditions. The difficulties associated with analysis are matched by those of characterizing hydrologic conditions in tight environments; measurements of hydraulic head and sampling of pore fluids have been done only rarely because of the practical difficulties involved. These problems are also discussed in the second part of this paper.

Water Resources Research

Characterization of springflow in the north coast limestone of Puerto Rico using physical, chemical, and stable isotopic methods

The carbonate sequence of middle Tertiary age of the north coast of Puerto Rico is characterized by the presence of numerous springs in the coastal areas. In order to advance the understanding of the hydrologic role of the springs in the north coast limestone aquifer system of Puerto Rico, a 4-year study was conducted by the U.S. Geological Survey in cooperation with the Puerto Rico Aqueduct and Sewer Authority. As part of this study, data were collected on the chemical, physical, bacteriological, oxygen-18, and deuterium composition of water from springs in the Dorado to Rincon area, in northwestern Puerto Rico. A group of springs in the Dorado to Arecibo area was selected for more detailed monitoring. Oxygen-18 and deuterium composition was also determined for water wells and monthly rainfall composites at a series of sites in the study area. Springs are associated with all the carbonate units of the middle Tertiary sequence of the northern karst belt of Puerto Rico, except the Camuy and San Sebastian Formations. These springs mostly drain the unconfined parts of the upper and lower aquifers in the north coast limestone aquifer system. There are no first and second order springs in the north coast limestone and of those present fifth and sixth order springs are the most numerous type. Springflow at the springs measured during the study ranged from less than 0.1 to 61 cubic feet per second. Springs in the north coast limestone can also be classified by their response to rainfall. There is little or no short-term response to rainfall at springs such as Ojo de Agua in Vega Baja, Mameyes in Manati, and Mackovic in Vega Alta. These springs are known as diffuse-type springs. Other springs such as Maguayo in Dorado, Ojo de Guillo in Manati, and San Pedro in Arecibo exhibit a strong short-term response to rainfall and are known as conduit-type springs. Spring water temperature, during the study, ranged from 22.5 to 28 °C and resembled air temperature. Specific conductance ranged from 289 to about 4,000 microsiemens per centimeter, and pH ranged from 6.9 to 7.8. Calcium, sodium, bicarbonate, and chloride are the main ionic species in water from the springs sampled during the study. The main water type is calcium-bicarbonate and secondary water types are calcium-bicarbonate chloride and sodium-bicarbonate-chloride. A seasonal and short-term transient relation exists, particularly in conduit-type springs, between springflow, physical properties, and water quality. Temporal and spatial variations in the oxygen-18 and deuterium composition of modern precipitation are significantly larger than those of springs and ground water in the study area. Regional flow in the upper aquifer appears to attenuate or average the variations in isotopic composition of rainfall. There is, however, a regional gradient in the deuterium composition of water from the upper aquifer in the north coast limestone, with isotopically heavier water occurring further north. It was possible to determine the source of water contributing to springs at some sites with more detailed data collection and analysis. A drainage basin of about 10 square kilometers was delineated for the Ojo de Agua spring in Vega Baja, for a base flow of 2 cubic feet per second and an estimated subregional recharge rate of around 20 inches per year. A delineation of drainage basins for conduit-type springs in the study area such as San Pedro, Ojo de Guillo, Maguayo, and others is very difficult because the boundaries of these systems are highly responsive to changing hydraulic conditions such as rapid and short-term variations in the hydraulic head distribution as a consequence of rainfall. However, a preliminary drainage basin of about 6 square kilometers was delineated for the San Pedro spring for a base flow of 2 cubic feet per second and an estimated subregional annual recharge of 12 inches.

Puerto Rico

Evaluation of candidate rain gages for upgrading precipitation measurement tools for the National Atmospheric Deposition Program

The National Atmospheric Deposition Program (NADP) was established in 1977 to investigate atmospheric deposition and its effects on the environment. Since its establishment, precipitation records have been obtained at all NADP sites using a gage developed approximately 50 years ago-the Belfort 5-780 mechanical rain gage. In 1998 and 1999, a study was done by the U.S. Geological Survey to evaluate four recently developed, technologically advanced rain gages as possible replacement candidates for the mechanical gage currently (2002) in use by the NADP. The gage types evaluated were the Belfort 3200, Geonor T-200, ETI Noah II, and the OTT PLUVIO. The Belfort 5-780 was included in the study to compare the performance of the rain gage currently (2002) used by NADP to the performance of the more recently developed gages. As a reference gage, the NovaLynx Model 260-2510 National Weather Service type stick gage also was included in the study. Two individual gages of each type were included in the study to evaluate precision between gages of the same type. A two-phase evaluation was completed. Phase I consisted of indoor bench tests with known amounts of simulated rainfall applied in 20 individual tests. Phase II consisted of outdoor testing by collecting precipitation during a 26-week period near Bay St. Louis, Mississippi. The ETI Noah II, OTT PLUVIO, and NovaLynx stick gages consistently recorded depths more commensurate with the amounts of applied simulated rainfall in Phase I testing than the Geonor T-200, Belfort 5-780, and Belfort 3200 gages. Gages where both the median difference between the measured and applied simulated rainfall and the interquartile range of all of their measured minus applied simulated rainfall differences were small (less than or equal to 0.01 inch) were judged to have performed very well in Phase I testing. The median and interquartile-range values were 0.01 inch or less for each of the ETI Noah II gages, OTT PLUVIO gages, and NovaLynx stick gages. The performance of the Geonor T-200 and Belfort 3200 gages was affected by technical problems during Phase I testing. As part of the evaluation of Phase II results, the average weekly precipitation totals obtained from the Belfort 5-780 gages and from each of the gages under consideration as possible replacements for the Belfort 5-780 gage were all compared with the average precipitation weekly totals obtained from two NovaLynx stick gages. The median absolute differences between a particular gage model and the NovaLynx stick reference gage for the 26 weeks of outdoor testing ranged from 0.04 inch for the ETI Noah II and OTT PLUVIO gages to 0.06 inch for the Geonor T-200. The total absolute difference between a particular gage type and the reference gage ranged from 1.23 inches for the Belfort 5-780 to 1.83 inches for the Geonor T-200 gages. Because the Belfort 3200 gages were inoperable for most of the Phase II testing, it is not meaningful to include the results from that gage type in a calculation of median or total absolute differences. The OTT PLUVIO proved to be the most reliable gage in Phase I and II testing, operating trouble free over the duration of the study.

Water-Resources Investigations Report

Historical methane hydrate project review

In 1995, U.S. Geological Survey made the first systematic assessment of the volume of natural gas stored in the hydrate accumulations of the United States. That study, along with numerous other studies, has shown that the amount of gas stored as methane hydrates in the world greatly exceeds the volume of known conventional gas resources. However, gas hydrates represent both a scientific and technical challenge and much remains to be learned about their characteristics and occurrence in nature. Methane hydrate research in recent years has mostly focused on: (1) documenting the geologic parameters that control the occurrence and stability of gas hydrates in nature, (2) assessing the volume of natural gas stored within various gas hydrate accumulations, (3) analyzing the production response and characteristics of methane hydrates, (4) identifying and predicting natural and induced environmental and climate impacts of natural gas hydrates, and (5) analyzing the effects of methane hydrate on drilling safety. Methane hydrates are naturally occurring crystalline substances composed of water and gas, in which a solid water-­‐lattice holds gas molecules in a cage-­‐like structure. The gas and water becomes a solid under specific temperature and pressure conditions within the Earth, called the hydrate stability zone. Other factors that control the presence of methane hydrate in nature include the source of the gas included within the hydrates, the physical and chemical controls on the migration of gas with a sedimentary basin containing methane hydrates, the availability of the water also included in the hydrate structure, and the presence of a suitable host sediment or “reservoir”. The geologic controls on the occurrence of gas hydrates have become collectively known as the “methane hydrate petroleum system”, which has become the focus of numerous hydrate research programs. Recognizing the importance of methane hydrate research and the need for a coordinated effort, the U.S. Congress enacted Public Law 106-­‐193, the Methane Hydrate Research and Development Act of 2000. This Act called for the Secretary of Energy to begin a methane hydrate research and development program in consultation with other U.S. federal agencies. At the same time a new methane hydrate research program had been launched in Japan by the Ministry of International Trade and Industry to develop plans for a methane hydrate exploratory drilling project in the Nankai Trough. Since this early start we have seen other countries including India, China, Canada, and the Republic of Korea establish large gas hydrate research and development programs. These national led efforts have also included the investment in a long list of important scientific research drilling expeditions and production test studies that have provided a wealth of information on the occurrence of methane hydrate in nature. The most notable expeditions and projects have including the following: -­‐Ocean Drilling Program Leg 164 (1995) -­‐Japan Nankai Trough Project (1999-­‐2000) -­‐Ocean Drilling Program Leg 204 (2004) -­‐Japan Tokai-­‐oki to Kumano-­‐nada Project (2004) -­‐Gulf of Mexico JIP Leg I (2005) -­‐Integrated Ocean Drilling Program Expedition 311 (2005) -­‐Malaysia Gumusut-­‐Kakap Project (2006) -­‐India NGHP Expedition 01 (2006) -­‐China GMGS Expedition 01 (2007) -­‐Republic of Korea UBGH Expedition 01 (2007) -­‐Gulf of Mexico JIP Leg II (2009) -­‐Republic of Korea UBGH Expedition 02 (2010) -­‐MH-­‐21 Nankai Trough Pre-­‐Production Expedition (2012-­‐2013) -­‐Mallik Gas Hydrate Testing Projects (1998/2002/2007-­‐2008) -­‐Alaska Mount Elbert Stratigraphic Test Well (2007) -­‐Alaska Iġnik Sikumi Methane Hydrate Production Test Well (2011-­‐2012) Research coring and seismic programs carried out by the Ocean Drilling Program (ODP) and Integrated Ocean Drilling Program (IODP), starting with the ODP Leg 164 drilling of the Blake Ridge in the Atlantic Ocean in 1995, have also contributed greatly to our understanding of the geologic controls on the formation, occurrence, and stability of gas hydrates in marine environments. For the most part methane hydrate research expeditions carried out by the ODP and IODP provided the foundation for our scientific understanding of gas hydrates. The methane hydrate research efforts under ODP-­‐IODP have mostly dealt with the assessment of the geologic controls on the occurrence of gas hydrate, with a specific goal to study the role methane hydrates may play in the global carbon cycle. Over the last 10 years, national led methane hydrate research programs, along with industry interest have led to the development and execution of major methane hydrate production field test programs. Two of the most important production field testing programs have been conducted at the Mallik site in the Mackenzie River Delta of Canada and in the Eileen methane hydrate accumulation on the North Slope of Alaska. Most recently we have also seen the completion of the world’s first marine methane hydrate production test in the Nankai Trough in the offshore of Japan. Industry interest in gas hydrates has also included important projects that have dealt with the assessment of geologic hazards associated with the presence of hydrates. The scientific drilling and associated coring, logging, and borehole monitoring technologies developed in the long list of methane hydrate related field studies are one of the most important developments and contributions associated with methane hydrate research and development activities. Methane hydrate drilling has been conducted from advanced scientific drilling platforms like the JOIDES Resolution and the D/V Chikyu, which feature highly advanced integrated core laboratories and borehole logging capabilities. Hydrate research drilling has also included the use of a wide array of industry, geotechnical and multi-­‐service ships. All of which have been effectively used to collect invaluable geologic and engineering data on the occurrence of methane hydrates throughout the world. Technologies designed specifically for the collection and analysis of undisturbed methane hydrate samples have included the development of a host of pressure core systems and associated specialty laboratory apparatus. The study and use of both wireline conveyed and logging-­‐while-­‐drilling technologies have also contributed greatly to our understanding of the in-­‐situ nature of hydrate-­‐bearing sediments. Recent developments in borehole instrumentation specifically designed to monitor changes associated with hydrates in nature through time or to evaluate the response of hydrate accumulations to production have also contributed greatly to our understanding of the complex nature and evolution of methane hydrate systems. Our understanding of how methane hydrates occur and behave in nature is still growing and evolving – we do not yet know if methane hydrates can be economically produced, nor do we know fully the role of hydrates as an agent of climate change or as a geologic hazard. But it is known for certain that scientific drilling has contributed greatly to our understanding of hydrates in nature and will continue to be a critical source of the information to advance our understanding of methane hydrates.

Report

Use case development for earth monitoring, analysis, and prediction (EarthMAP)—A road map for future integrated predictive science at the U.S. Geological Survey

Executive Summary The U.S. Geological Survey (USGS) 21st-century science strategy 2020–30 promotes a bureau-wide strategy to develop and deliver an integrated, predictive science capability that works at the scales and timelines needed to inform societally relevant resource management and protection and public safety and environmental health decisions (U.S. Geological Survey, 2021). This is the overarching goal of the USGS Earth Monitoring, Analysis, and Prediction (EarthMAP) vision, which consists of three components: (1) integrated data and information, (2) integrated predictive science, and (3) actionable information—all designed and delivered to respond to user needs. To launch this vision and help shape the design and implementation of integrated predictive science, the USGS Regional Offices each developed a set of use cases (hereafter Use Cases)—short descriptions of potential science applications that could clearly address high priority decision-making needs of our stakeholders and that align with an integrated science focus. Use Cases are not actionable science planning documents, nor stand-alone scholarly works, but should be considered as innovative, next-generation science ideas that can be considered as potential components of science plans still under development. The goal of Use Case development was to (1) identify and characterize existing USGS scientific capacities and expertise that can support science goals and products, (2) identify opportunities to leverage current capacities for next-generation science, and (3) foster engagement across the entire Bureau to further refine the USGS strategy for EarthMAP and integrated predictive science. The Use Case development effort documented in this report was coordinated by the Use Case Development Team (UCDT), consisting of representatives from each region. The UCDT undertook five tasks: (1) develop a unified approach to engage bureau scientists consistently across all regions in aspirational thinking about what can be accomplished; (2) work with the regions and their Science Centers to generate an initial set of Use Cases, authored directly by scientists; (3) characterize, summarize, and document the initial set of Use Case submissions from authors to illuminate bureau-level demand for integrated science; (4) compare existing and needed capacities from the Use Case descriptions with preliminary results of the EarthMAP Capacity Assessment (Keisman and others, 2021); and (5) describe lessons learned from the Use Case development process and provide recommendations to inform future efforts to generate integrated science activities. This report outlines the approach the UCDT developed to solicit Use Cases from the regions and summarizes the high-level qualitative findings from this first-round effort. The UCDT received 36 Use Cases from the regions and identified potential points of convergence and commonalities considered useful in making connections among the participating scientists. The Southwest (SW) Region and the Rocky Mountain (RM) Region asked scientists to give special consideration to Use Cases with applicability to the Colorado River Basin, and seven of the Use Cases specifically named that geographic area as a focus. Coastal hazards and coastal resilience were identified in Use Cases from the Alaska (AK), Northeast (NE), and Southeast (SE) Regions. Aspects of wildfire and post-wildfire response were part of Uses Cases from AK, RM, and SW Regions. The greatest convergence of Use Case themes was related to conservation of public lands and waters, which is a powerful linkage lending strength to future collaborative efforts. The most common type of stakeholder decisions that would be informed by the Use Case science applications were related to adaptation, mitigation, and response (for example, how to increase the resilience of coastal communities to climate-related stressors and how to prevent or respond to harmful algal blooms). Other common types of decisions included water and land management decisions (including operational water management decisions such as reservoir operations and land use planning in the sagebrush biome), decisions about how to manage and conserve habitats and species, and risk management decisions (such as managing the post-wildfire flood risks). These decision types are not exclusive because many Use Cases cross categories. Use Case authors identified existing and needed science and technology capabilities required for Use Case implementation, which were then aligned to capabilities assessed in the EarthMAP Capacity Assessment (Keisman and others, 2021). Strong alignment was found for data and information integration approaches, modeling and prediction approaches, and capabilities related to delivery of actionable information. A majority of Use Cases indicated insufficient current capacity for needed data collection methods, data integration, and modeling and prediction approaches, whereas only 25 percent indicated insufficient capacity for actionable information delivery. Overall, many Use Case capacity demand gaps could potentially be met by existing bureau-wide capacity. In addition, nearly half of the Use Cases could potentially be implemented within 3 years if funding, capabilities, and personnel impediments were removed and science priorities were realigned. Several challenges emerged during the Use Case development process. The first challenge was developing an approach that was flexible enough to accommodate regional differences in planning and implementation, while also ensuring enough guidance to promote meaningful summary analyses. The UCDT encountered a strong demand for continuous communication and education to improve overall understanding of the integrated predictive science strategy. Another challenge was managing expectations about EarthMAP activities as a design effort that was not aligned to an immediate funding opportunity. Connecting the Use Cases to stakeholder needs without the opportunity for direct stakeholder engagement was also challenging. The last notable challenge was in obtaining consistent interpretation and characterization of the qualitative data housed in the narrative descriptions of Use Cases, written in different styles. Overall, the 36 Use Cases can serve as components of a road map for advancing integrated monitoring and predictive science throughout the USGS by revealing opportunities to (1) encourage cross-region initiatives that address shared interests in common themes by integrating similar Use Cases and through direct involvement of stakeholders in identifying needs and designing effective responses, (2) leverage the Use Cases to target investments that are aligned with the Bureau and Department of the Interior (DOI) priorities, (3) connect Use Cases and the results of the companion EarthMAP Capacity Assessment (Keisman and others, 2021) to identify potential priorities for capacity building investments, and (4) raise awareness of common integrated and interdisciplinary science interests within and across the regions through Use Case and Capacity Assessment summary outreach activities.

Open-File Report

Salt-front movement in the Hudson River estuary, New York — Simulations by one-dimensional flow and solute-transport models

The Hudson River is being considered for use as a supplemental source of water supply for New York City during droughts. One proposal entails withdrawal of Hudson River water from locations near Newburgh, Chelsea, or Kingston, but the extent to which this could cause the salt front to advance upstream to points where it could adversely affect community water supplies is unknown. The U.S. Geological Survey (USGS) one-dimensional Branch-Network Dynamic Flow model (BRANCH) was used in conjunction with the USGS one-dimensional Branched Lagrangian Solute-Transport Model (BLTM) to simulate the effect of five water-withdrawal scenarios on the salt-front location. The modeled reach contains 132 miles of the lower Hudson River between the Federal Dam at Troy and Hastings-on-Hudson (near New York City). The BRANCH model was calibrated and verified to 19 tidal-cycle discharge measurements made at 11 locations by conventional and acoustic Doppler current-profiler methods. Maximum measured instantaneous tidal flow ranged from 20,000 ft 3 /s (cubic feet per second) at Albany to 368,000 ft 3 /s at Tellers Point; daily-mean flow at Green Island near Troy ranged from 3,030 ft 3 /s to 45,000 ft 3 /s during the flow measurements. Successive ebb- and flood-flow volumes were measured and compared with computed volumes; daily-mean bias was -1.6 percent (range from -21.0 to +23.7 percent; 13.5 percent mean absolute error). Daily-mean deviation between simulated and measured stage at eight locations (from Bowline Point to Albany) over the 19 tidal-cycle measurements averaged +0.06 ft (range from -0.31 to +0.40 ft; 0.21 ft root mean square error, RMSE). These results indicate that the model can accurately simulate flow in the Hudson River under a wide range of flow, tide, and meteorological conditions. The BLTM was used to simulate chloride transport in the 61-mi reach from Turkey Point to Bowline Point under two seasonal conditions in 1990.one representing spring conditions of high inflow and low salinity (April-June), the other representing typical summer conditions of low inflow and high salinity (July-August). Measured chloride concentrations at Bowline Point were used to drive the BLTM simulations, and data collected at West Point were used for calibration. Mean bias in simulated chloride concentration for the April-June 1990 (high flow) data (observed range from 12 to 201 mg/L [milligrams per liter]; 30 mg/L RMSE) was .16 mg/L, and mean bias for the July-August 1990 (low flow) data (observed range from 31 to 2,000 mg/L; 535 mg/ L RMSE) was +126 mg/L. The salt front (saltwater/ freshwater interface) on the Hudson River was defined as the furthest upstream location where the chloride concentration exceeded 100 mg/L. Data from August 1991 were used to evaluate solute transport between West Point and Poughkeepsie because a chloride concentration of 100 mg/L was not observed at Clinton Point in 1990. The BLTM then was used to simulate chloride concentrations at Chelsea Pump Station and Clinton Point. Regression equations, based on daily mean values of specific conductance measured at West Point, were used to estimate daily mean chloride concentrations at Chelsea Pump Station and Clinton Point for model analysis. Mean biases in BLTM-simulated daily mean chloride concentrations for August 1991 were .38 mg/L at Chelsea Pump Station (range from 189 to 551 mg/L; 103 mg/L RMSE) and .9 mg/L at Clinton Point (range from 53 to 264 mg/L; 62 mg/L RMSE). Hypothetical withdrawals at (1) Newburgh, (2) Chelsea, (3) Chelsea and Newburgh, (4) Chelsea and Kingston, and (5) Kingston and Newburgh, were simulated to compute the effects of withdrawals on salt-front movement. Withdrawals of 300 Mgal/d from any combination of Chelsea or Newburgh could result in upstream movement of the salt front of as much as 1.0 mi, given an initial salt-front location between West Point and Rogers Point. Scenarios that included withdrawals at Kingston caused the greatest upstream salt-front movement. Simulation of a 90-day April-June high-flow period during which discharges at Green Island averaged 25,200 ft 3 /s indicated that withdrawals of 1,939 Mgal/d (million gallons per day) at Chelsea Pump Station would not measureably increase chloride concentrations at Chelsea Pump Station under normal tidal and meteorological conditions, but withdrawals at twice that rate (3,878 Mgal/d) could increase the chloride concentration at Chelsea Pump Station to 250 mg/L.

New York

Santa Barbara and Foothill groundwater basins Geohydrology and optimal water resources management—Developed using density dependent solute transport and optimization models

Groundwater has been a part of the city of Santa Barbara’s water-supply portfolio since the 1800s; however, since the 1960s, the majority of the city’s water has come from local surface water, and the remainder has come from groundwater, State Water Project, recycled water, increased water conservation, and as needed, seawater desalination. Although groundwater from the Santa Barbara and Foothill groundwater basins only accounts for a small percentage of the long-term supply, it is an important source of supplemental water during times of surface-water shortages. During the late 1980s and early 1990s, production wells extracted additional groundwater to compensate for drought related water-delivery shortfalls from other sources; in response, water levels declined substantially in the Santa Barbara and Foothill groundwater basins (below sea level in the Santa Barbara groundwater basin). In coastal basins that have groundwater extraction near shore, seawater intrusion is often a problem. Seawater intrusion in the Santa Barbara groundwater basin is thought to be more limited than in other coastal basins because of an offshore fault that acts as a partial barrier to groundwater flow. During the late 1980s and early 1990s, seawater intrusion was observed in the Santa Barbara groundwater basin, as indicated by increased chloride concentrations at several monitoring wells that ranged from 200 ft to 1,300 ft from the ocean and as close as 2,900 ft to the nearest pumping well. This demonstrated that seawater can intrude into the Santa Barbara groundwater basin when groundwater levels fall below sea level near the coast. The city of Santa Barbara is interested in developing a better understanding of the sustainability of its groundwater supplies. In 2014, California adopted historic legislation to manage its groundwater: the Sustainable Groundwater Management Act (SGMA). The SGMA requires the development and implementation of “Groundwater Sustainability Plans” in 127 priority groundwater basins; although Santa Barbara was not a designated priority basin, the city is taking steps to achieve sustainability. Sustainability was defined in the SGMA in terms of avoiding undesirable results: significant and unreasonable groundwater-level declines, reduction in groundwater storage, seawater intrusion, water-quality degradation, land subsidence, and surface-water depletion. In this project, a cooperative study between the U.S. Geological Survey (USGS) and the city of Santa Barbara, sustainable yield is defined as the volume of groundwater that can be pumped from storage without causing water-level drawdowns and the associated increases in seawater intrusion (as indicated by increases in measured chloride concentrations) at selected wells. In order to estimate the sustainability of Santa Barbara’s groundwater basins, a three-dimensional density-dependent groundwater-flow and solute-transport model (the Santa Barbara Flow and Transport Model, or SBFTM) was developed on the basis of an existing groundwater-flow model. To simulate seawater intrusion to the Santa Barbara Basin under various management strategies, the SBFTM uses the USGS code SEAWAT to simulate salinity transport and variable-density flow. The completed SBFTM was coupled with a management optimization tool, in this case a multi-objective evolutionary algorithm, to determine optimal pumping strategies that maximize the sustainable yield and at the same time satisfy user-defined drawdown and chloride-concentration constraints. As part of this study, a three-dimensional hydrogeologic framework model was developed to quantify the extent and hydrogeologic characteristics of the Santa Barbara and Foothill groundwater basins and to help define the discretization and hydraulic properties used in the SBFTM. The development of the hydrogeologic framework model required the collection and reconciliation of geologic and geophysical data from existing maps, reports, and databases, along with geologic and hydrologic data from recently drilled wells. These data were integrated into a three-dimensional hydrogeologic framework model that defines the stratigraphy and geometry of the aquifer zones and the major geologic structures in the basin. The hydrogeologic framework model also quantifies the variation in sediment grain size within each aquifer zone as the percentage of coarse-grained sediment. Previous studies indicated that there are two principal water-producing zones in the Santa Barbara groundwater basin, the upper and lower producing zones; an additional thin, productive zone was identified as part of this study. This “middle producing zone” is not as areally extensive as the upper and lower producing zones and only exists in the coastal part of Storage Unit I. These producing zones are bounded at depth by less productive shallow, middle, and deep zones. Two versions of the SBFTM were constructed: an initial-condition model and a modern transient model. The initial-condition model is a long-term transient model that simulates flow and solute-transport conditions during a period with limited anthropogenic influences preceeding the modern transient model. The simulation-transient model simulates flow and transport conditions from 1929 through 2013; however, because of data availability, the focus of the model calibration was 1972–2013. The SBFTM was calibrated to measured groundwater levels and drawdown, as well as measured chloride concentrations and change in concentrations, using a combination of automated and trial-and-error parameter-estimation techniques. A sensitivity analysis indicated that, in general, the SBFTM was most sensitive to recharge- and pumping-distribution parameters, specifically those controlling the amount of small-catchment recharge and the distribution of water extraction by hydrogeologic layer for production wells. The model was also sensitive to parameters controlling stream-recharge rates, horizontal and vertical hydraulic conductivity, and porosity. From 1929 to 1971, most of the water entering the area represented by the SBFTM was from creek and small-catchment recharge, and the majority of water leaving the SBFTM area was from pumping, discharge to creeks, and drains. In addition, about 37 percent of the total pumpage came from a net reduction in groundwater storage. From 1972 to 2013, the amount of water entering and leaving the SBFTM was fairly similar as that from 1929 to 1971, except the reduction in pumpage added about 17,000 acre-ft of water to storage. During this later period, there were also times of storage loss. For example, during July 1990, a month when approximately 705 acre-ft of groundwater was pumped in the study area, the pumpage was much greater than all sources of recharge combined, and about 382 acre-ft of water was removed from groundwater storage. Simulated hydraulic heads replicated the observed data to an acceptable matching of the measured water-level, flow direction, and vertical gradients. Simulated hydrographs for selected wells were in good agreement with the measured data, with an average residual of -2.7 ft and a standard deviation of 14.5 ft, indicating that the simulated heads, on average, underestimated the observed water levels. An examination of the model fit indicated that most of the discrepancies were lower simulated heads at wells proximal to production well sites. The simulated chloride concentrations reasonably matched the rising limbs of the measured breakthrough curves in terms of timing and magnitude; however, the simulation overestimated the chloride concentrations on the falling limbs. The overestimation of low chloride concentrations was attributed to the model overestimating the advance of the chloride front during periods of heavy pumping and underestimating the retreat of the chloride front during periods of low pumping. These simulation errors would result in a conservative response by local water managers to seawater intrusion. The SBFTM was used to develop a collection of predictive simulations optimized to produce pumping schedules that maximize yield, subject to a set of constraints and competing objectives. The simulations were grouped as scenarios that differed in their time horizon, initial conditions for groundwater levels and chloride concentrations, as well as precipitation, which was incorporated into the model through simulated recharge. Overall, five scenarios were developed in a multi-objective framework to obtain optimal pumping rates for all of the wells managed by the city, while minimizing excessive drawdown and seawater intrusion. For the current study, complexities in the simulation model and the optimization formulation required additional considerations. Incorporating the solute-transport equations to simulate chloride transport added a highly nonlinear process that is solved iteratively in each time step of the groundwater-flow model. These nonlinearities, coupled with the highly refined grid in the current model, creates challenges for many traditional optimization methods. Therefore, an optimization method was needed that could address nonlinear relationships as well as a very large problem size. Lastly, the optimization problem was reformulated to include multiple objectives without requiring convergence to a single solution. This approach, guided by the city’s objectives, allowed the maximum extraction of information from the complex simulation. Borg, a multi-objective evolutionary algorithm, was chosen as the optimization algorithm for this study for several reasons: (1) it is very computationally efficient; (2) it can run in parallel; (3) it requires little user input; and (4) it can solve for multiple competing objectives. The first three points allow the algorithm to proceed toward the optimal solutions at the fastest possible rate. The fourth point is advantageous for large, complex optimization problems because it is difficult to formulate the optimization problem in a way that produces only one optimal solution. The problem formulation consisted of four competing objectives and a constraint set in accordance with the main concerns of the city. The objectives were maximizing total pumpage, minimizing seawater intrusion, minimizing total drawdown in production wells, and minimizing the maximum drawdown. The constraints were pump capacity, meeting drinking-water standards for chloride, maintaining a specified minimum flowrate to a groundwater treatment plant, and maintaining minimum water levels in pumping wells. The decision variables either were quarterly pumpage by well or total pumpage by basin. Five optimization scenarios were developed that allow the decision makers to evaluate a range of optimal solutions for a variety of water levels and chloride concentrations as well as potential future climatic conditions. Three scenarios (1, 2, and 5) were multi-objective optimization formulations that allowed for variations in management preferences and climatic conditions. The other two scenarios (3 and 4) were designed to examine the optimization results to answer specific questions. Scenario 1 described the best-case sustainable yield assuming a “full” basin (that is, high initial water levels) and typical climate conditions for 10 years. Scenario 2 also started with a “full” basin; however, this was followed by a 10-year drought. Scenario 3 determined if an “empty” basin (that is, low initial water levels) would recover to full conditions (1998 conditions) given climate assumptions and optimal pumping schedules from scenarios 1 and 2. Scenario 4 was designed to produce decision rules that can be used by water managers to help choose an optimal pumping schedule based on measured water-level or chloride data. Scenario 5 identified future pumping schedules based on short-term climate variations during a 2-year management horizon. The results from scenarios 1 and 2 described the differences in maximum pumpage in the basin under typical and dry long-term climate projections, respectively. The scenario 1 results indicated the maximum 10-year pumpage of the basin was about 31,300 acre-ft under typical conditions and controlling simulated seawater intrusion and drawdowns. For scenario 2, less recharge over the 10-year dry climate produced a maximum pumpage estimate of 30,000 acre-ft to control seawater intrusion and drawdowns. The larger pumpage for scenario 1 resulted in more seawater intrusion, but less total drawdown, compared to that of scenario 2. Results for scenarios 3 and 4 showed the basin’s response to management actions combined with climate projections. Both scenarios used the optimal pumping schedules and the 10-year climates from scenarios 1 and 2. The scenario 3 results showed that under minimal pumping, the basin did not fully recover to 1998 water levels within 10 years under either climate scenario. The relatively larger recharge from the typical climate resulted in less drawdown at coastal monitoring wells after the 10-year recovery period than that from the dry climate. The location of the seawater intrusion front was not appreciably different between the scenarios, however. Scenario 4 used the optimal results from scenarios 1 and 2 to produce decision-rule curves that illustrated the pumpage for each basin, given measured levels of chloride concentration or drawdown. This allowed the use of additional measurements at monitoring wells to assess future management decisions on the basis of the sensitivity of observations of drawdown and seawater intrusion to various pumping rates. Scenario 5 allowed managers to investigate the effects of short-term climate variations on optimal pumping schedules. Three specific 2-year simulations were optimized: typical-to-dry (scenario 5A), dry-to-typical (scenario 5B), and dry-to-dry (scenario 5C). The most noteable result from scenario 5 was the overall reduction in optimal pumpage for most schedules in scenario 5C, when the climate is simulated as dry-to-dry. There are also many optimal pumping schedules that produced an overall increase in waterlevels over the two-year simulation period, regardless of climatic condition. Similar to scenario 2, the scenario 5C results represents conservative yield estimates under a minimal-precipitation climatic condition.

California

Borehole-geophysical and hydraulic investigation of the fractured-rock aquifer near the University of Connecticut Landfill, Storrs, Connecticut, 2000 to 2001

An integrated borehole-geophysical and hydraulic investigation was conducted at the former landfill area near the University of Connecticut in Storrs, Connecticut, where solvents and landfill leachate have contaminated a fractured-bedrock aquifer. Borehole-geophysical techniques and hydraulic methods were used to characterize the site bedrock lithology and structure, fractures, and hydraulic properties. The geophysical and hydraulic methods included conventional logs, borehole imaging, borehole radar, flowmeter under ambient- and stressed hydraulic conditions, and discrete-zone hydraulic testing, sampling, and monitoring. The conventional geophysical-logging methods included caliper, deviation, electromagnetic induction, gamma, specific conductance, and fluid temperature. The advanced methods included optical and acoustic imaging of the borehole wall, heat-pulse flowmeter, and directional radar reflection. Borehole-geophysical methods were used to further define conductive features identified with surface-geophysical methods in the first phase of the investigation. The results of the surface- and borehole-geophysical logging were evaluated in an iterative and integrated manner to develop a conceptual model of ground-water flow at the site. The rock type, foliation, and fractures at the site were characterized from high-resolution optical televiewer (OTV) images of rocks penetrated by the boreholes and were compared to drilling logs and conventional geophysical logs. The rocks are interpreted as fine- to mediumgrained quartz-feldspar-biotite-garnet gneiss and schist with local intrusions of quartz diorite and pegmatite and minor concentrations of sulfide mineralization similar to rocks described as the Bigelow Brook Formation on regional geologic maps. Layers containing high concentrations of sulfide minerals appear as high electrical conductivity zones on electromagnetic-induction and borehole-radar logs. Foliation in the rocks generally strikes to the southwest and northeast, and dips to the northwest and southeast consistent with previous investigations in this area. The orientation of foliation, however, varies locally and with depth in some of the boreholes. These results are consistent with geologic mapping that has identified small-scale folding. The orientations of the transmissive fractures identified in the six boreholes logged for this investigation are similar to the fracture orientations mapped in a previous investigation. Many of these fractures are oriented with a north-northwest strike and have a shallow dip to the west. Other transmissive fractures have a southwest strike and dip at shallow angles to the northwest, and some strike roughly east-west and dip to the north and south. Flowmeter logging was used to identify transmissive fractures and to estimate the hydraulic properties in the boreholes. Ambient down flow was measured in one borehole, and ambient up flow and down flow were measured in another borehole. The other four bedrock boreholes did not have measurable vertical flow. Under low-rate pumping conditions (0.25 to 0.5 gallons per minute), one to three inflow zones were identified in each well. Commonly, fractures that are active under ambient conditions contribute to the well under pumping conditions. The ambient conditions were incorporated into the determination of the relative proportions of transmissivity. Specific capacity and transmissivity were determined for these open-hole low-rate pumping tests. Quasi-steady-state water levels were reached in four of the boreholes, including MW201R, MW204R, MW302R, and W202-NE. When pumped at low-rate conditions for 0.5 to 4 hours, the specific capacity ranged from 0.03 to 0.18 gallons per minute per foot. The open-hole transmissivity estimates ranged from 4.9 to 30 feet squared per day (ft2/d). Open-hole transmissivity was determined for boreholes that did not reach quasi-steady-state conditions under low-rate pumping conditions. Transmissivity was estimated for MW201R, MW202R, and MW203R using non-equilibrium methods, pumping rate, and the transient drawdown data to estimate the open-hole transmissivity. Transmissivity in these boreholes ranged from 0.98 to 3.2 ft 2 /d. The transmissivity and head of individual fractures or zones of fractures were estimated from heat-pulse flowmeter data acquired under ambient and stressed conditions. In the absence of ambient flow, data from two profiles of heat-pulse flowmeter data under two different stressed conditions were used to estimate the transmissivity and head of individual fracture zones. Only two boreholes, MW302R and W202-NE, had sufficient data for these analyses. The estimated transmissivity of individual transmissive zones ranged from 1.2 to 9.2 ft 2 /d. The transmissivity values determined by this numerical simulation method were less than the open-hole estimations, which were 15 and 30 ft 2 /d. Transmissivity also was measured directly over discrete intervals of the borehole using a straddle-packer apparatus and constant-rate pumping tests. Pumping rates were less than or equal to 0.25 gallons per minute. These discretezone single-hole pumping tests were conducted over a short period of time, usually about 30 minutes to 1 hour in duration. Pumping continued until the test zone reached a steady-state water level or until it was determined that the zone could not yield water at the pumped rate. The estimated transmissivity of individual transmissive zones ranged from about 0.21 to 11 ft 2 /d. The zone at a depth of 197 feet in W202-NE was the only zone that had discrete-interval testing with a straddle packer and sufficient heat-pulse flowmeter data for modeling the flow and estimating transmissivity and head. The two methods produced similar results. The straddle-packer method estimated a transmissivity of 4.7 ft 2 /d, and the heat-pulse flowmeter modeling results estimated a transmissivity of 6.9 ft 2 /d. A comparison of the transmissivity estimates indicate estimates typically are within an order of magnitude. The heat-pulse flowmeter methods used in this investigation to determine transmissivity of the boreholes and the individual fractures measure only the upper two or three orders of magnitude of transmissivity. Hence, other fractures in these boreholes permit the movement of water; their transmissivities, however, are lower than the detection limits of the methods that were used for this investigation and very small compared to the transmissive fractures that were studied. The data collected in this investigation were used to design discrete-zone monitoring systems for four of the boreholes used for monitoring. The results of the investigation are useful for refining the conceptual site model of ground-water flow, and for providing critical information for interpreting the results of water-quality sampling.

Connecticut

The geography and geology of Alaska; a summary of existing knowledge, with a section on climate, and a topographic map and description thereof

Alaska, the largest outlying possession of the United States, is that great land mass forming the northwestern extremity of the North American continent, whose western point is within 60 miles of the Asiatic coast (PI. II). About one-quarter of this area lies within the Arctic Circle, and from the standpoint of geographic position must be regarded as an arctic province; but the southern seaboard, exposed to the warm winds and waters of the Pacific, gives to the entire southern portion of the territory" a comparatively warm climate. It is not generally realized that the range of climate in Alaska is greater than that between Florida and Maine. At the southernmost point of the Pacific coast the mean annual temperature is not far from that of the city of Washington, the winters being warmer and characterized by less snowfall; the Yukon Valley on the other hand has a winter climate similar to that of northern Montana and Dakota; while in the extreme northern part of the territory the meteorologic conditions are invariably arctic. Though as yet only sparsely settled, Alaska's vast area and great resources make it one of the most important possessions of the United States and promise its rapid development. During the years 1890 to 1900 the population increased from 32,052* to 63,592. The mineral output, which in 1890 was valued at less than $800,000, exceeded $9,000,000 in 1904, and the fisheries show a corresponding growth. This rapid development has attracted public attention and led to urgent demand for explorations, surveys, and other investigations. So actively has this work been pushed, both by public and private enterprise, that exact knowledge of the geography, geology, and mineral resources of the interior has made greater strides within the last eight years than during the preceding thirty-one years since the acquisition of Alaska. The facts regarding the geography and geology, scattered as they are through the many books and reports of this period, are not always readily accessible, and the time seems ripe to present them in a summarized form. The topography of Alaska is varied and complex (see PI. I), and it is not easy to present briefly even the salient features. The limited number of pages here devoted to the subject precludes the possibility of detailed treatment, even if the facts were available. Much of the description has been taken from the results attained by other investigators, the writer being personally familiar with only a part of this large province. A list of the publications consulted is appended. The larger geographic features of Alaska are now fairly well known, though the detailed surveys which are demanded by the development of many localities have hardly been begun. Preliminary surveys have been completed of all but three 8 of the larger rivers. The most important mountain ranges have been at least outlined (fig. 3). Only three large areas remain almost entirely unmapped: One in southwestern Alaska, between Cook Inlet and the lower Kuskokwim, and the others in northern Alaska, embracing the Arctic watershed east and west of the Colville River. Nearly all the surveys of the interior, however, have been of a preliminary and exploratory character, and to meet the requirements of exact geography must be followed by more detailed mensuration. Though the coast line has been fairly well known for more than half a century, knowledge of the interior has been gained chiefly within the last two decades. This has not yet found its way into text-books and has too often been entirely ignored by cartographers. If facts are presented which may seem elementary, it is because even well-informed people have been known to harbor misconceptions in regard to the orographic features, climate, and general character of Alaska. Those who read of the perils and privations of winter travel and explorations are apt to picture a region of ice and snow; others, again, who have personal knowledge of the tourist route of southeastern Alaska, regard the whole district as one of rugged mountains and glaciers. In point of fact, glaciers are now nearly limited to the ranges bordering the Pacific and to the two slopes of the Alaska range; and even during the greatest development of glaciers but a small portion of Alaska was under ice (see map, PI. XXII). As a treatise on geography would hardly be complete without some discussion of the climate, meteorologic data have been compiled by Mr. Cleveland Abbe, jr., but the discussion of this does not pretend to be more than a cursory treatment of the subject. The scope of the paper seems to require also a brief summary of the development of geographic knowledge of Alaska. This subject, with its many ramifications, is of fascinating interest and offers a magnificent field for the trained historian. If the accompanying sketch of discovery and exploration awakens any measure of popular interest the writer will feel amply rewarded for having attacked a theme which hardly falls within the scope of his investigations. When this compilation was begun it was intended to be chiefly a description of the topography of Alaska, as illustrated by the accompanying map (PI. XXXIV, in pocket), which was compiled under the direction of the late R. U. Goode. In the course of the work there accumulated much geologic as well as geographic material which seemed worthy of inclusion in the report. As no comprehensive statement of the geology of Alaska has been made since the modern epoch of investigation was begun, an attempt will be made to give a summary of all results achieved. Since the writer has obtained much of his knowledge of the facts from the work of others, he disclaims any pretense of making an entirely original contribution to geologic science. He feels, however, that a personal familiarity with a considerable part of the province, gained during seven consecutive seasons of field work, will justify Mm in presenting conclusions which may in some cases be at variance with those in the reports on which he must draw for his facts. Throughout this report attempt will be made to credit borrowed material to the source from which it is drawn. Where such matter has been obtained entirely from published reports there is no difficulty in so doing; but as regards investigators of the Geological Survey, with whom the writer has collaborated both in field and in office, the case is somewhat different, for it is not always possible to know whether this or that theory originated with the writer or with one of his colleagues. It will, then, perhaps suffice to state that this report could not have been prepared without the explorations and researches of the geologists, F. C. Schrader, Walter C. Mendenhall, Arthur J. Collier, J. E. Spurr, and Arthur C. Spencer; and the surveys of the topographers, T. G. Gerdine, D. C. Witherspoon, D. L. Reaburn, W. J. Peters, and E. C. Barnard. Each of these men, in the course of from two to six years of field work, has made important contributions to the knowledge of the geography and geology of Alaska, and not all of these results have yet been put in print. In the last season (1903) L. M. Prindle, C. W. Wright, Arthur Hollick, G. C. Martin, F. L. Hess, and Fred H. Moffit have carried on geologic work in Alaska, and the writer has made use of their work now in course of publication. He has also been fortunate in having access to the manuscript reports of Walter C. Mendenhall and F. C. Schrader on the Copper River basin, to which references will be made. The matter here presented should be credited in a measure to all of these investigators, but for many of the theories advanced the writer alone is responsible. As this manuscript goes to press there has been opportunity to incorporate some of the results of the field work of 1904. As far as possible these have been embodied in the text, but in some instances it has been found advisable to add them only as footnotes. During the past summer F. E. and C. W. Wright extended the geologic reconnaissance in southeastern Alaska. In southwestern Alaska G. C. Martin and T. W. Stanton have determined the general Mesozoic section, while F. H. Moffit has made a reconnaissance of the northern part of the Kenai Peninsula. A. J. Collier has mapped the geology of the Cape Lisburne region, and L. M. Prindle and F. L. Hess have made contributions to the knowledge of the metamorphic terranes of the Yukon-Tanana district. It is the writer's purpose to describe in nontechnical language the larger geographic features and discuss their relation as far as the data available will permit. In the treatment of the geology, however, less effort will be made to make the matter acceptable to the lay reader. It is hoped, however, that a brief summary of the salient features of the geologic history' may be not without interest to the general public. If this paper serves in some measure to dispel the popular fallacies regarding Alaska and to disseminate more accurate knowledge of its geographic and geologic features, the purpose of its publication will be accomplished.

Alaska

Borehole-geophysical investigation of the University of Connecticut landfill, Storrs, Connecticut

A borehole-geophysical investigation was conducted to help characterize the hydrogeology of the fractured-rock aquifer and the distribution of unconsolidated glacial deposits near the former landfill and chemical waste-disposal pits at the University of Connecticut in Storrs, Connecticut. Eight bedrock boreholes near the landfill and three abandoned domestic wells located nearby were logged using conventional and advanced borehole-geophysical methods from June to October 1999. The conventional geophysical-logging methods included caliper, gamma, fluid temperature, fluid resistivity, and electromagnetic induction. The advanced methods included deviation, optical and acoustic imaging of the borehole wall, heat-pulse flowmeter, and directional radar reflection. Twenty-one shallow piezometers (less than 50-feet deep) were logged with gamma and electromagnetic induction tools to delineate unconsolidated glacial deposits. Five additional shallow bedrock wells were logged with conventional video camera, caliper, electromagnetic induction, and fluid resistivity and temperature tools. The rock type, foliation, and fracturing of the site were characterized from high-resolution optical-televiewer (OTV) images of rocks penetrated by the boreholes. The rocks are interpreted as fine- to medium-grained quartz-feldsparbiotite-garnet gneiss and schist with local intrusions of quartz diorite and pegmatite and minor concentrations of sulfide mineralization similar to rocks described as the Bigelow Brook Formation on regional geologic maps. Layers containing high concentrations of sulfide minerals appear as high electrical conductivity zones on electromagneticinduction and borehole-radar logs. Foliation in the rocks generally strikes to the northeast-southwest and dips to the west, consistent with local outcrop observations. The orientation of foliation and small-scale gneissic layering in the rocks, however, varies locally and with depth in some of the boreholes. In two of the boreholes, the foliation strikes predominantly to the northwest and dips to the northeast. Although small-scale faults and lithologic discontinuities were observed in the OTV data, no large-scale faults were observed that appear on regional geologic maps. Fractures were located and characterized through the use of conventional geophysical, OTV, acoustic-televiewer (ATV), and borehole-radar logs. The orientation of fractures varies considerably across the site; some fractures are parallel to the foliation, whereas others cross-cut the foliation. Many of the transmissive fractures in the bedrock boreholes strike about N170°E and N320°E with dips of less than 45°. Other transmissive fractures strike about N60°E with dips of more than 60°. Most of the transmissive fractures in the domestic wells strike about N60°E and N22°E with dips of more than 45°. The strike of N60°E is parallel to the trend of a thrust fault that appears on regional geologic maps. Vertical flow in the boreholes was measured with the heat-pulse flowmeter under ambient and (or) pumping conditions. Results of ATV, OTV, and conventional logs were used to locate specific zones for flowmeter testing. Ambient downflow was measured in three boreholes, ambient upflow was measured in two other boreholes, and both ambient downflow and upflow were measured in a sixth borehole. The other five bedrock boreholes and domestic wells did not have measurable vertical flow. The highest rate of ambient flow was measured in the background borehole in which upflow and downflow converged and exited the borehole at a fracture zone near a depth of 62 feet. Ambient flow of about 340 gallons per day was measured. In the other five wells, ambient flow of about 20 to 35 gallons per day was measured. Under low-rate pumping (0.25 to 1 gallon per minute), one to six inflow zones were identified in each well. Usually the fractures that are active under ambient conditions contribute to the well under pumping conditions. To prevent ambient vertical flow and the potential for cross-contamination, temporary borehole liners were installed in five of the boreholes. Specific-capacity and open-hole transmissivity values were determined in eight boreholes completed in bedrock. The specific capacity estimated for these boreholes ranges from 0.14 to 1.6 gallons per minute per foot. The values for open-hole transmissivity range over two orders of magnitude and when proportioned to individual fracture transmissivity, range from 23 to 340 feet squared per day. Two boreholes had been drilled to intersect electrically conductive zones identified by previous surface-geophysical investigations. The borehole-geophysical results indicate that the boreholes penetrate electrically conductive structures consistent with the anomalies interpreted from the surface-geophysical data. Borehole MW121R was located to intersect a dipping electrically conductive anomaly at about 60 feet, interpreted from the two-dimensional direct currentresistivity survey conducted on the western side of the landfill. The electromagnetic-conductivity log in the borehole contains a high electrical conductivity anomaly at a depth of 69 feet. The magnitude of this anomaly is nearly 10,000 millisiemens per meter and is coincident with a layer containing sulfide mineralization, rather than fractures. The other borehole, MW105R, was located to intersect another anomaly south of the landfill. This anomaly was interpreted as a north-south striking, westward dipping feature. In the borehole, two south-striking, westward dipping fractures were identified in the ATV, OTV, and radar logs. The specific conductance of the fluid measured near these fractures was as high as 1,250 microsiemens per centimeter. Water-quality samples collected in October 1999 from an isolated zone from 71.5 to 76.5 feet indicated high specific conductance (810 microsiemens per centimeter), high concentrations of iron and cadmium, negative oxidation-reduction potential, and chlorobenzene. Collectively, these parameters indicate that the high specific conductance in the borehole logs for MW105R was caused by landfill leachate. Therefore, the anomaly identified by boreholeand surface-geophysical surveys is interpreted as a conductive lithologic feature and a permeable fracture zone that contains landfill leachate.

Connecticut