USGS Science⌕ Search

SEARCH · USGS Science

Results for “Computer contribution”

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.

474 records · Page 27Linked to original sources

Science mission requirements for a globally ranging, riserless drilling vessel for U.S. Scientific Ocean Drilling

Through the collection and analysis of shallow and deep subseafloor sediments, rocks, fluids, and life, scientific ocean drilling has enriched our understanding of the complex Earth system. Among other achievements, scientific ocean drilling has documented the history of Earth’s climate, the waxing and waning of polar ice sheets, the past changes in ocean and atmospheric circulation, the existence and function of microbial life in the subseafloor, the compositional variations in Earth’s crust and underlying mantle, and the physical and chemical processes acting at subduction zones, including those associated with tsunamigenic earthquakes. Over the decades, more than 12,000 articles that depend on analyses of scientific ocean drilling samples and geophysical data have been published, many detailing breakthrough contributions to global knowledge about the Earth system. Approximately 45% of these publications were led by U.S.-affiliated authors (International Ocean Discovery Program Publication Services, 2021). Since the mid-1980s, the workhorse of this multidisciplinary, international research effort has been the riserless D/V JOIDES Resolution, operated by Texas A&M University with funding from the U.S. National Science Foundation (NSF). D/V JOIDES Resolution has conducted the vast majority of scientific ocean drilling expeditions and collected most of the scientific cores over that period, including 82% of the expeditions and 93% of the cores in the last decade alone, despite being one of three platforms that is operated within the International Ocean Discovery Program. However, D/V JOIDES Resolution is approaching the end of its useful life. With a strong commitment to continue scientific ocean drilling beyond the end of the current phase, the community developed a document outlining the research frontiers that should be pursued. Exploring Earth by Scientific Ocean Drilling: 2050 Science Framework (Koppers and Coggon, 2020) describes seven scientific strategic objectives that focus on understanding interconnections within the Earth system and five flagship initiatives that integrate these objectives into long-term research efforts that address issues facing society. Additional elements in the 2050 Science Framework, including STEM education, workforce development, technology development, and innovative applications of data analytics, will advance the goals of scientific ocean drilling. Addressing the 2050 Science Framework also requires building partnerships with allied U.S. and international science programs and strengthening existing ones. To implement a significant portion of the 2050 Science Framework, the U.S. scientific community seeks to lease or acquire a newly built, globally ranging, state-of-the art, riserless drilling vessel. The many and varied technical and human resources requirements for successful accomplishment of scientific and educational goals summarized in this document and described in detail in the 2050 Science Framework require broad community input and careful consideration. Following receipt of NSF’s formal Request for Assistance to the United States Science Support Program (USSSP), the U.S. scientific ocean drilling community conducted a one-year exercise to identify its national scientific needs and priorities in order to determine the Science Mission Requirements (SMRs) presented here. This community effort included: (1) a U.S. community-wide survey to identify the specific operational and technical capabilities critical to addressing science in the 2050 Science Framework; (2) a series of online workshops focusing on critical capabilities identified by the survey; and (3) a large in-person workshop to synthesize the results of the survey and the virtual workshops (Appendix 1). The approach was designed to reach as many participants as possible. Overall, 278 survey responses were received from U.S. community members, representing 104 unique institutions from 39 states and the District of Columbia, and 137 unique individuals participated in the workshops (Appendix 2). The results of this effort comprise two classes of SMRs: Foundational Science Mission Requirements and Primary Science Mission Requirements. Foundational SMRs define minimum criteria for a new riserless drilling vessel that can address significant portions of the 2050 Science Framework. Primary SMRs build upon the Foundational SMRs and will create more robust science opportunities and data collection capabilities, will increase progress in addressing the 2050 Science Framework objectives, and will provide more real-time ship-to-shore interaction to improve science productivity, engagement, and outreach. Modern safety and environmental standards, including meeting standards to access protected waters such as exclusive economic zones, extended continental shelves, or high latitudes, while being cognizant of the vessel’s environmental footprint. Safe and efficient operations in global locations and in water depths from 70 m to 6000+ m, with total drill string length of at least 7000 m. High-quality core and data collection from a range of key subseafloor environments. Advanced heave compensation, dynamic positioning, and drill pipe stability. Modern mud and cement/casing systems. Critical onboard measurements for safety, operational decision-making, documentation of ephemeral properties, mission-specific science, and long-term science goals that extend beyond a single expedition. Designated and appropriate space for sample and data preservation. Highly skilled onboard personnel, including technical staff for curation and core handling; support for safety, time-sensitive, and critical shipboard measurements; computer support; equipment and instrument repair; application support; and data assurance. Primary Science Mission Requirements include: Flexible shipboard space for laboratories and on deck to ensure safe, successful implementation of diverse science objectives and operations. Minimizing contamination of recovered samples. Over-the-side capabilities for science-supporting technology (e.g., remotely operated vehicles, water- column sampling, sediment-water interface sampling). Downhole logging and measurements. Expanded borehole observatory capabilities. Reliable and consistent ship-to-shore communications. NSF’s investment in a new globally ranging, riserless drilling vessel will have a powerful economic multiplier effect, including the infusion of additional science support funds in the United States for training and research, the development of new technologies and tools, and the associated scientific and technical workforce development. The skills and knowledge gained through scientific ocean drilling are translatable to careers in fields such as sustainable energy development (e.g., geothermal and offshore wind), carbon sequestration, data management and cyberinfrastructure, biotechnology, communications, science education, policy, hazard mitigation, and environmental management. The United States is a leader in a well-established and internationally collaborative scientific ocean drilling community. A modern, globally ranging, riserless drilling vessel will allow the United States to expand its leadership position, address broad scientific questions that current capabilities preclude, and cultivate equitable international, multidisciplinary collaborations that will ensure scientific ocean drilling’s future success.

Report↗

Hydrogeology and Simulation of Groundwater Flow in the Plymouth-Carver-Kingston-Duxbury Aquifer System, Southeastern Massachusetts

The glacial sediments that underlie the Plymouth-Carver-Kingston-Duxbury area of southeastern Massachusetts compose an important aquifer system that is the primary source of water for a region undergoing rapid development. Population increases and land-use changes in this area has led to two primary environmental effects that relate directly to groundwater resources: (1) increases in pumping that can adversely affect environmentally sensitive groundwater-fed surface waters, such as ponds, streams, and wetlands; and (2) adverse effects of land use on the quality of water in the aquifer. In response to these concerns, the U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, began an investigation in 2005 to improve the understanding of the hydrogeology in the area and to assess the effects of changing pumping and recharge conditions on groundwater flow in the Plymouth-Carver-Kingston-Duxbury aquifer system. A numerical flow model was developed based on the USGS computer program MODFLOW-2000 to assist in the analysis of groundwater flow. Model simulations were used to determine water budgets, flow directions, and the sources of water to pumping wells, ponds, streams, and coastal areas. Model-calculated water budgets indicate that approximately 298 million gallons per day (Mgal/d) of water recharges the Plymouth-Carver-Kingston-Duxbury aquifer system. Most of this water (about 70 percent) moves through the aquifer, discharges to streams, and then reaches the coast as surface-water discharge. Of the remaining 30 percent of flow, about 25 percent of the water that enters the aquifer as recharge discharges directly to coastal areas and 5 percent discharges to pumping wells. Groundwater withdrawals are anticipated to increase from the current (2005) rate of about 14 Mgal/d to about 21 Mgal/d by 2030. Pumping from large-capacity production wells decreases water levels and increases the potential for effects on surface-water bodies, which are affected by pumping and wastewater disposal locations and rates. Pumping wells that are upgradient of surface-water bodies potentially capture water that would otherwise discharge to these surface-water bodies, thereby reducing streamflow and pond levels. The areas most affected by proposed increases in groundwater withdrawals are in the Towns of Plymouth and Wareham where more than half of the proposed increase in pumping will occur. In response to an increase of about 7 Mgal/d of pumping, groundwater discharge to streams is reduced by about 6 cubic feet per second (ft3/s) (about 4 Mgal/d) from a total of about 325 ft3/s. Reduction in streamflow is moderated by an increase of artificial recharge from wastewater returned to the aquifer by onsite domestic septic systems and centralized wastewater treatment facilities. It is anticipated that about 3 Mgal/d of the 7 Mgal/d of increase in pumped water will be returned to the aquifer as wastewater by 2030. Currently (2005) about 3 percent of groundwater discharge to streams is from wastewater return flow to the aquifer during average conditions. During drought conditions, the component of streamflow augmented by wastewater return flow doubles as wastewater recharge remains constant and aquifer recharge rates decrease. Wastewater return flow, whether as direct groundwater discharge to streams or as an additional source of aquifer recharge, increases the height of the water table near streams, thereby moderating the effects of increased groundwater withdrawals on streamflow. An analysis of a simulated drought similar to the 1960s drought of record indicates that the presence of streams moderates the effects on water levels of reduced aquifer recharge. The area where water-table altitudes were least affected by drought was in the Weweantic River watershed in the Town of Carver. Water levels decreased by less than 2 feet from current average conditions compared to decreases of greater than 5

Scientific Investigations Report↗

Groundwater/surface-water interactions in the Partridge River Basin and evaluation of hypothetical future mine pits, Minnesota

The Partridge River Basin (PRB) covers 156 square miles in northeastern Minnesota with headwaters in the Mesabi Iron Range. The basin is characterized by extensive wetlands, lakes, and streams in poorly drained and often thin glacial material overlying Proterozoic bedrock. To better understand the interaction between these extensive surface water features and the groundwater system, a three-dimensional, steady-state, groundwater-flow model of the PRB was developed by the U.S. Geological Survey in cooperation with the Great Lakes Indian Fish & Wildlife Commission using the finite-difference computer code MODFLOW-NWT. The model simulates steady-state base flow in streams and groundwater interactions using the streamflow routing (SFR2) package. Existing mining features including tailings basins, stockpiles, pumped mine pits, and flooded mine pits were simulated using either high hydraulic conductivity zones or the drain (DRN) package. The unsaturated zone flow (UZF) package was used to better represent the groundwater system in areas with a high water table and for wetlands often associated with such areas. UZF typically is used to represent unsaturated zone processes but also can simulate the rejection of recharge and groundwater discharge to the land surface when the water table is near land surface. The steady-state model used data from the 2011 to 2013 period when 2011 high-resolution land surface (light detecting and ranging [lidar]) data were available that reflected land-surface and water elevations from mining activity in the basin. The parameter-estimation software suite PEST_HP was used to obtain a best fit of the modeled to measured groundwater levels, streamflow, pit inflow rates, and mapped peat deposits. The PEST calibration used the target residuals from two models with the same model parameters and targets from two separate periods: (1) a 1995–2015 calibration model, which provided a larger number of calibration targets, and (2) a 2011–2013 mining conditions model, which included calibration targets that reflected conditions consistent with the modeled mine-workings topography. Calibration of the PRB model resulted in ranges of glacial horizontal hydraulic conductivity parameters that generally agreed with literature values and other models of the region. Horizontal hydraulic conductivity of the bedrock was higher in the upper bedrock layers where numerous and continuous fractures have been observed and lower in the deeper bedrock layers. Average basin-wide calibrated infiltration was 5.3 inches per year. An average of 4.6 inches per year of infiltration crosses the water table and becomes recharge and 0.7 inch per year is rejected by UZF due to saturated conditions at the land surface. Simulated groundwater runoff (the sum of rejected recharge and groundwater seepage to the land surface) can either be routed to streams or removed from the model as evapotranspiration. The calibrated model indicates relatively shallow groundwater-flow paths dominating and approximately 50 percent of the stream base flow coming from groundwater runoff. The 2011–2013 mining conditions model was then used to develop five model scenarios simulating the response of the groundwater and surface-water system to potential hydrologic stress. The purpose of these mine pit scenarios is to present a possible workflow to quantify a model’s uncertainty for a given model forecast and serve as a possible guide for initial data collection that may improve a future model’s ability to make such a forecast. The scenarios included one scenario with the currently existing Peter Mitchell pit at final buildout and flooded to an elevation of 1,500 feet, and four scenarios with a hypothetical, new mine pit plus the flooded Peter Mitchell at final buildout. The five model scenarios were used to forecast streamflow at six locations in the PRB, pit inflow rates for the new mine pits and the flooded Peter Mitchell pit, and the average depth to water in 12 wetlands. A linear uncertainty analysis was performed using information from the PEST calibration and tools in the PyEMU python package to assess model uncertainty propagation to the model forecasts. Streamflows generally were reduced with future mining and the greatest streamflow reductions occurred from the flooded Peter Mitchell Pit, probably due to its large size. Average depth to groundwater in wetlands was most affected the closer the wetland was to a new mine pit. Linear uncertainty methods were also used to evaluate data worth, which is the ability for potential new groundwater elevation observations to reduce the uncertainty in scenario forecasts. Data worth was performed for a grid of new hydraulic head observations. Overall, areas with nonnegligible data worth generally corresponded to wetland areas with no groundwater seepage to land surface from UZF. These model behaviors indicated that the land-surface boundary condition simulated by the UZF package was pinning the groundwater elevations to the land surface in areas with groundwater seepage (33 percent of the 2011–2013 base conditions model) such that the sensitivity to new observations in these areas was minimal. Therefore, representing wetlands as boundary conditions minimized the usefulness of data worth calculations because wetland areas were present over a large part of the model domain. Probabilistic capture zones were estimated for each of the mines in the model scenarios. A capture zone represents the area contributing recharge to a model feature, like a well or a mine pit, and can be calculated by forward tracking particles from the water table. By using Monte Carlo techniques, it is possible to generate estimated capture zones that include the probability of recharge capture given the uncertainty present in the model. Monte Carlo techniques use randomly generated model parameter sets sampled from a plausible parameter range to create many possible realizations. The resulting capture zone arrays were calculated by tallying the total number of realizations in which a particle from a model cell was captured by the feature. Probabilities from the Monte Carlo runs ranged from 1 (captured in 100 percent of the runs) near the pits to 0 (captured in 0 percent of the runs) at the edges of the capture zone. Capture zones were not always spatially continuous; for example, the capture zone for the proposed mine pits south of the flooded Peter Mitchell pit was discontinuous with capture surrounding the proposed mine pit and north of the flooded Peter Mitchell pit. This northern section represents deeper groundwater flow paths that originate in the topographic high, move under the flooded pit, and discharge into the proposed pit. This pattern of capture indicates the possibility of some deeper flow through the upper fractured bedrock when the shallow groundwater flow system is modified. These results underscore that future site-specific applications of the base condition model require the input of site-specific data and recalibration to focus on the site of interest.

Minnesota↗

Predicting sea-level rise vulnerability of terrestrial habitat and wildlife of the Northwestern Hawaiian Islands

If current climate change trends continue, rising sea levels may inundate low-lying islands across the globe, placing island biodiversity at risk. Recent models predict a rise of approximately one meter (1 m) in global sea level by 2100, with larger increases possible in areas of the Pacific Ocean. Pacific Islands are unique ecosystems home to many endangered endemic plant and animal species. The Northwestern Hawaiian Islands (NWHI), which extend 1,930 kilometers (km) beyond the main Hawaiian Islands, are a World Heritage Site and part of the Papahanaumokuakea Marine National Monument. These NWHI support the largest tropical seabird rookery in the world, providing breeding habitat for 21 species of seabirds, 4 endemic land bird species and essential foraging, breeding, or haul-out habitat for other resident and migratory wildlife. In recent years, concern has grown about the increasing vulnerability of the NWHI and their wildlife populations to changing climatic patterns, particularly the uncertainty associated with potential impacts from global sea-level rise (SLR) and storms. In response to the need by managers to adapt future resource protection strategies to climate change variability and dynamic island ecosystems, we have synthesized and down scaled analyses for this important region. This report describes a 2-year study of a remote northwestern Pacific atoll ecosystem and identifies wildlife and habitat vulnerable to rising sea levels and changing climate conditions. A lack of high-resolution topographic data for low-lying islands of the NWHI had previously precluded an extensive quantitative model of the potential impacts of SLR on wildlife habitat. The first chapter (chapter 1) describes the vegetation and topography of 20 islands of Papahanaumokuakea Marine National Monument, the distribution and status of wildlife populations, and the predicted impacts for a range of SLR scenarios. Furthermore, this chapter explores the potential effects of SLR on wildlife breeding habitats for each island. The subsequent chapter (chapter 2) details a study of the Laysan Island ecosystem, describing a quantitative model that incorporates SLR, storm wave, and rising groundwater inundation. Wildlife, storm, and oceanographic data allowed for an assessment of the phenological and spatial vulnerability of Laysan Island's breeding bird species to SLR and storms. Using remote sensing and geospatial techniques, we estimated topography, classified vegetation, modeled SLR, and evaluated a range of climate change scenarios. On the basis of high-resolution airborne data collected during 2010-11 (root-mean-squared error = 0.05-0.18 m), we estimated the maximum elevation of 20 individual islands extending from Kure Atoll to French Frigate Shoals (range: 1.8-39.7 m) and computed the mean elevation (1.7 m, standard deviation 1.1 m) across all low-lying islands. We also analyzed general climate models to describe rainfall and temperature scenarios expected to influence adaptation of some plants and animals for this region. Outcomes for the NWHI predicted an increase in temperature of 1.8-2.6 degrees Celsius (°C) and an annual decrease in precipitation of 24.7-76.3 millimeters (mm) across the NWHI by 2100. Our models of passive SLR (excluding wave-driven effects, erosion, and accretion) showed that approximately 4 percent of the total land area in the NWHI will be lost with scenarios of +1.0 m of SLR and 26 percent will be lost with +2.0 m of SLR. Some atolls are especially vulnerable to SLR. For example, at Pearl and Hermes Atoll our analysis indicated substantial habitat losses with 43 percent of the land area inundated at +1.0 m SLR and 92 percent inundated at +2.0 m SLR. Across the NWHI, seven islands will be completely submerged with +2.0 m SLR. The limited global ranges of some tropical nesting birds make them particularly vulnerable to climate change impacts in the NWHI. Climate change scenarios and potential SLR impacts presented here emphasize the need for early climate change adaptation and mitigation planning, especially for species with limited breeding distributions and/or ranges restricted primarily to the low-lying NWHI: Cyperus pennatiformis var. bryanii , Black-footed Albatross ( Phoebastria nigripes ), Laysan Albatross ( P. immutabilis ), Bonin Petrel ( Pterodroma hypoleuca ), Gray-backed Tern ( Onychoprion lunatus ), Laysan Teal ( Anas laysanensis ), Laysan Finch ( Telespiza cantans ), and Hawaiian monk seal ( Monachus schauinslandi ). Furthermore, SLR scenarios that include the effects of wave dynamics and groundwater rise may indicate amplified vulnerability to climate change driven habitat loss on low-lying islands. In chapter 2, we incorporated the combined effects of SLR, dynamic wave-driven inundation, and rising groundwater in a quantitative study specifically for the Laysan Island ecosystem. This is the first hydrodynamic model to simulate the combined impacts of SLR and wave-driven inundation in the NWHI. We developed a high-resolution digital elevation model (mean vertical accuracy of 0.32 m) for the island. Then using recent satellite imagery, geospatial models, and historical oceanographic, storm, and biological data we estimated potential inundation extent, habitat loss, and wildlife population impacts for a range of potential SLR scenarios (0.00, +0.50, +1.00, +1.50, and +2.00 m) that may occur over the next century. Additionally, we estimated the carrying capacity of Laysan Island for five species based on the available population monitoring data and described how potential losses in nesting habitat could influence population dynamics for Black-footed Albatross, Laysan Albatross, Red-footed Booby (Sula sula), Laysan Teal, and Laysan Finch. For some other seabird populations (Masked Booby, S. dactylatra ; Brown Booby, S. leucogaster ; Great Frigatebird, Fregata minor ; and Sooty Tern, Onychoprion fuscata ), we used recent colony distribution data, land cover maps, and nesting behavior to estimate potential losses of nesting habitat from SLR and wave-driven inundation. We observed far greater potential impacts of SLR to wildlife with the dynamic wave-driven modeling approach than with the passive modeling approach. Depending on SLR scenario and coastal orientation, during storms under a +2.00 m SLR scenario, the wave-driven inundation model predicted three times more inundation than the passive model (17.2 percent of total terrestrial area versus 4.6 percent, respectively). Large-wave events generally added 1 m of water height to passive inundation surfaces, therefore our dynamic models (during storm events) forecasted comparable inundation extents earlier than passive models. Although wave-driven water levels were highest in the northwest quadrant of Laysan Island, the greatest extent of inundation occurred in the southeast where coastal dunes less than 3 m above mean sea level provide little protection from wave-driven inundation. When wave-driven inundation was included in the SLR model for Laysan Island greater nesting habitat loss and potential impacts on wildlife population dynamics were predicted. The consequences of habitat loss due to SLR may be worse for species with colonies in the wave-exposed coastal zones (for example, Black-footed Albatross) and for populations already near the island's carrying capacity (for example, Laysan Teal). Species whose peak incubation and chick-rearing periods coincide with seasonally high wave heights also will be increasingly vulnerable, especially those species nesting on the ground in areas vulnerable to inundation, such as Gray-backed Tern and Black-footed Albatross. Other species that have space for population growth, or are not restricted to a narrow range of habitat types on Laysan (for instance, Sooty Terns), may be less sensitive to habitat loss from SLR over the next century. Our assessments of inundation risk, habitat loss, and wildlife species vulnerability synthesize current knowledge about individual islands and contribute to a broader understanding of the impacts of inundation from SLR and storm-induced waves. Yet, most NWHI and their bird populations lack monitoring data to evaluate adaptations to and impacts of climate change. Exceptions include some data sets from long-term monitoring of wildlife populations, tides, or weather at French Frigate Shoals, Laysan Island, and Midway Atoll. These data sets are potentially valuable baselines, which could be informative for adaptive learning (integrating management and science) to predict, adapt, and mitigate the effects of climate change on NWHI wildlife in the future. This study provides the first quantitative vulnerability assessment for all of the low-lying NWHI, and results identify biological communities, locales, and resident endangered species of Papahanaumokuakea Marine National Monument expected to be at risk from SLR. This report is also intended as a reference for managers and conservation planners, a tool to identify and potentially reduce uncertainty, and a starting place for developing climate change monitoring priorities and future scientific studies.

Hawai'i↗

Investigations of groundwater system and simulation of regional groundwater flow for North Penn Area 7 Superfund site, Montgomery County, Pennsylvania

Groundwater in the vicinity of several industrial facilities in Upper Gwynedd Township and vicinity, Montgomery County, in southeast Pennsylvania has been shown to be contaminated with volatile organic compounds (VOCs), the most common of which is the solvent trichloroethylene (TCE). The 2-square-mile area was placed on the National Priorities List as the North Penn Area 7 Superfund site by the U.S. Environmental Protection Agency (USEPA) in 1989. The U.S. Geological Survey (USGS) conducted geophysical logging, aquifer testing, and water-level monitoring, and measured streamflows in and near North Penn Area 7 from fall 2000 through fall 2006 in a technical assistance study for the USEPA to develop an understanding of the hydrogeologic framework in the area as part of the USEPA Remedial Investigation. In addition, the USGS developed a groundwater-flow computer model based on the hydrogeologic framework to simulate regional groundwater flow and to estimate directions of groundwater flow and pathways of groundwater contaminants. The study area is underlain by Triassic- and Jurassic-age sandstones and shales of the Lockatong Formation and Brunswick Group in the Mesozoic Newark Basin. Regionally, these rocks strike northeast and dip to the northwest. The sequence of rocks form a fractured-sedimentary-rock aquifer that acts as a set of confined to partially confined layers of differing permeabilities. Depth to competent bedrock typically is less than 20 ft below land surface. The aquifer layers are recharged locally by precipitation and discharge locally to streams. The general configuration of the potentiometric surface in the aquifer is similar to topography, except in areas affected by pumping. The headwaters of Wissahickon Creek are nearby, and the stream flows southwest, parallel to strike, to bisect North Penn Area 7. Groundwater is pumped in the vicinity of North Penn Area 7 for industrial use, public supply, and residential supply. Results of field investigations by USGS at the site and results from other studies support, and are consistent with, a conceptual model of a layered leaky aquifer where the dip of the beds has a strong control on hydraulic connections in the groundwater system. Connections within and (or) parallel to bedding tend to be greater than across bedding. Transmissivities of aquifer intervals isolated by packers ranged over three orders of magnitude [from about 2.8 to 2,290 square feet per day (ft 2 /d) or 0.26 to 213 square meters per day (m 2 /d)], did not appear to differ much by mapped geologic unit, but showed some relation to depth being relatively smaller in the shallowest and deepest intervals (0 to 50 ft and more than 250 ft below land surface, respectively) compared to the intermediate depth intervals (50 to 250 ft below land surface) tested. Transmissivities estimated from multiple-observation well aquifer tests ranged from about 700 to 2,300 ft 2 /d (65 to 214 m 2 /d). Results of chemical analyses of water from isolated intervals or monitoring wells open to short sections of the aquifer show vertical differences in concentrations; chloride and silica concentrations generally were greater in shallow intervals than in deeper intervals. Chloride concentrations greater than 100 milligrams per liter (mg/L), combined with distinctive chloride/bromide ratios, indicate a different source of chloride in the western part of North Penn Area 7 than elsewhere in the site. Groundwater flow at a regional scale under steady-state conditions was simulated by use of a numerical model (MODFLOW-2000) for North Penn Area 7 with different layers representing saprolite/highly weathered rock near the surface and unweathered competent bedrock. The sedimentary formations that underlie the study area were modeled using dipping model layers for intermediate and deep zones of unweathered, fractured rock. Horizontal cell model size was 100 meters (m) by 100 meters (328 ft by 328 ft), and model layer thickness ranged from 6 m (19.7 ft) representing shallow weathered rock and saprolite up to 200 m (656 ft) representing deeper dipping bedrock. The model did not include detailed structure to account for local-scale differences in hydraulic properties, with the result that local-scale groundwater flow may not be well simulated. Additional detailed multi-well aquifer tests would be needed to establish the extent of interconnection between intervals at the local scale to address remediation of contamination at each source area. This regional groundwater-flow model was calibrated against measured groundwater levels (1996, 2000, and 2005) and base flow estimated from selected streamflow measurements by use of nonlinear-regression parameter-estimation algorithms to determine hydraulic conductivity and anisotropy of hydraulic conductivity, streambed hydraulic conductivity, and recharge during calibration periods. Results of the simulation using the calibrated regional model indicate that the aquifer appears to be anisotropic where hydraulic conductivity is greatest parallel to the orientation of bedding of the formations underlying the area and least in the cross-bed direction. The maximum hydraulic conductivity is aligned with the average regional strike of the formations, which is “subhorizontal” in the model because the altitudes of the beds and model cells vary in the strike, as well as dip, direction. Estimated subhorizontal hydraulic conductivities (in strike direction parallel to dipping beds) range from 0.001 to 1.67 meters per day (0.0032 to 5.5 feet per day). The ratio of minimum (dip direction) to maximum (strike direction) subhorizontal hydraulic conductivity ranges from 1/3.1 to 1/8.6, and the ratio of vertical to horizontal hydraulic conductivity ranges from 1/1 to 1/478. However, limited available field data precluded rigorous calibration of vertical anisotropy in the model. Estimated recharge rates corresponding to calibration periods in 1996, 2000, and 2005 are 150, 109, and 124 millimeters per year (5.9, 4.3, and 4.9 inches per year), respectively. The calibrated groundwater-flow model was used to simulate groundwater flow under steady-state conditions during periods of relatively high withdrawals (pumpage) (1990) and relatively low withdrawals (2000 and 2005). Groundwater-flow paths originating from recharge areas near known areas of soil contamination (sources) were simulated. Pumped industrial and production wells captured more groundwater from several of these sources during 1990 than after 1990 when pumping declined or ceased and greater amounts of contaminated groundwater moved away from North Penn Area 7 Superfund site to surrounding areas. Uncertainty in simulated groundwater-flow paths from contaminant sources and contributing areas, resulting from uncertainty in estimated hydraulic properties of the model, was illustrated through Monte Carlo simulations. The effect of uncertainty in the vertical anisotropy was not included in the Monte Carlo simulations. Contributing areas indicating the general configuration of groundwater flow towards production well MG-202 (L-22) in the study area also were simulated for the different time periods; as simulated, the flow paths do not pass through any identified contaminant source in North Penn Area 7. However, contributing areas to wells, such as MG-202, located near many pumped wells are particularly complex and, in some cases, include areas that contribute flow to streams that subsequently recharge the aquifer through stream loss. In these cases, water-quality constituents, including contaminants that are present in surface water may be drawn into the aquifer to nearby pumped wells. Results of a simulated shutdown of well MG-202 under steady-state 2005 conditions showed that the area contributing recharge for nearby production well MG-76 (L-17), when MG-202 is not pumping, shifts downstream and is similar to the area contributing recharge for MG-202 when both wells are pumping. Concentrations of constituents in groundwater samples collected in fall 2005 or spring 2006 were compared to simulated groundwater-flow paths for the year 2005 to provide a qualitative assessment of model results. The observed spatial distribution of selected constituents, including TCE, CFC-11, and CFC-113 in groundwater in 2005 and the chloride/bromide mass ratios in 2006, generally were consistent with the model results of the simulated 2005 groundwater-flow paths at North Penn Area 7, indicating the presence of several separate sources of contaminants within North Penn Area 7.

Pennsylvania↗

Hydrology and water and sediment quality at James Campbell National Wildlife Refuge near Kahuku, Island of Oahu, Hawaii

The James Campbell National Wildlife Refuge occupies two lowland marsh and pond complexes on the northern coastal plain of Oahu: the mostly natural ponds and wetlands of the Punamano Unit and the constructed ponds of the Kii Unit. The U.S. Fish and Wildlife Service manages the Refuge primarily to protect and enhance habitat for four endangered species of Hawaiian waterbirds. Kii Unit is fed by artesian wells and rainfall, whereas Punamano Unit is fed naturally by rainfall, runoff, and ground-water seepage. Streams drain from the uplands into lowland ditches that pass through Kii Unit on their way to the ocean. A high-capacity pump transfers water from the inner ditch terminus at Kii to the ocean outlet channel. Stormwaters also exit the inner ditch system over flood-relief swales near the outlet pump and through a culvert with a one-way valve. A hydrologic investigation was done from November 1996 through February 1998 to identify and quantify principal inflows and outflows of water to and from the Refuge, identify hydraulic factors affecting flooding, document ground-water/surface-water interactions, determine the adequacy of the current freshwater supply, and determine water and sediment quality. These goals were accomplished by installing and operating a network of stream-gaging stations, meteorology stations, and shallow ground-water piezometers, by computing water budgets for the two Refuge units, and by sampling and analyzing water and pond-bottom sediments for major ions, trace metals, and organic compounds. Streamflow during the study was dominated by winter stormflows, followed by a gradual recession of flow into summer 1997, as water that had been stored in alluvial fans drained to lowland ditches. Outflow at the ditch terminus in 1997 was 125 million gallons greater than measured inflow to the coastal plain, mainly reflecting gains from ground water along the ditches between outlying gages and the ditch terminus. Of the measured 1997 outflow, 98 percent was through the Kii outlet pump, with the outlet culvert valve only opening for brief periods during storms. Large volumes of stormflow overflowed the flood-relief swales unmeasured. The largest storm of the study, in November 1996, was estimated to have a flood frequency of about 3 to 4 years. Streamflow exceeded culvert capacity and overtopped Kamehameha Highway at Kalaeokahipa Stream and Hospital ditch. Slight overbank flooding in Kii ditch resulted strictly from high discharge. Minor overbank flooding farther out on the coastal plain probably was caused mainly by the small hydraulic gradients available to convey stormflows along the lowland ditches. Stormwaters flooded Kii ponds and flowed back upstream along Punamano ditch into Punamano marsh, introducing suspended sediment and possibly other contaminants to the Refuge. Two smaller storms in January 1997 resulted in smaller flows and no overbank flooding. The Kii outlet pump ran continuously for 7 days during the November 1996 storm and for 1 to 2 days during the January 1997 storms. During all three storms, the outlet culvert valve opened and the inner ditches overtopped the flood-relief swales, allowing free outflow of water from the inner ditch. Backwater effects hindered drainage during the January 1997 storms at Hospital ditch at Kamehameha Highway, and at Punamano ditch at Nudist Camp Road (where the backflow into Punamano marsh in November 1996 constituted an extreme backwater effect). A probable marine backwater effect was imposed at the ocean outlet ditch during the November 1996 storm through a combination of high spring tides and wave setup from large surf. Whether this backwater effect propagated upstream in the ditches to affect inland sites could not be determined conclusively. A sand plug may have built up in the ocean outlet channel before the November 1996 storm, but if so, it probably washed out prior to, or early in the storm, and was not present at the time of peak stage at inland sites. A season-long buildup of the sand plug in late 1997 was inferred from rising water levels in the outlet ditch. Seawater flows up the outlet channel or over the sand berm and into the outer ditch system on most high tides, and particularly during spring high tides. Ponds and ditches of the Refuge and surrounding lowlands have mud- and clay-lined bottoms that form an effective confining unit and inhibit interaction with an underlying shallow limestone aquifer. At Kii Unit, pond levels are higher than adjacent ditch levels and underlying ground-water levels, establishing lateral and downward head gradients that could foster seepage losses from the ponds. Regional ground-water discharge from the Koolau aquifer to the coastal-plain sediments is mostly diffuse, but is concentrated where ridges of Koolau Basalt plunge beneath coastal-plain sediments near Punamano Unit and at the head of Hospital ditch. Kii ditch gains brackish ground water downstream of Kamehameha Highway. Wastewater disposal from the sewage treatment plant adjoining Kii Unit poses little or no threat to Refuge habitat. Disposal is at six injection wells located 0.45 mi away at Kahuku, and the wells inject into confined limestone aquifers that do not extend to Kii Unit. The natural freshwater supply to Punamano Unit is adequate for maintaining the wildlife habitat, judging from stable pond levels and low salinities there. A monthly water budget for Punamano showed an apparent annual deficit in measurable flows in 1997, requiring unmeasured ground-water gains equalling 51 inches of water. The freshwater supply to Kii Unit is inadequate according to Refuge managers, because there is not enough water to manipulate levels adequately in the ponds during most of the year, and particularly during the driest months. This is confirmed by monthly deficits in the water budget for the Kii ponds during summer months. However, the Kii budget showed an annual surplus in measurable flows for 1997 equalling 24 inches of water. Unmeasured losses are required to explain the apparent annual surplus, such as discharge to the ditches through pond water-control structures and downward and/or lateral ground-water seepage. The apparent surplus at Kii is strictly hydrologic and is not a surplus in a management sense; it cannot be stored or used to supply the Refuge, but instead reflects losses from the system that render this amount of water unavailable for use. The budget year, 1997, was drier than normal (24 percent below long-term mean rainfall) and so the measured potential evaporation for 1997 was probably higher than the long-term mean. Few metals or organic compounds of potential concern were detected in pond and ditch waters and in pond-bottom sediments. Detected pesticides were at trace levels or just above minimum reporting limits. Exceptions that exceeded quality guidelines for freshwater sediment were copper and zinc in sediment from Kii ponds C and D, and copper in sediment from Punamano north pond. Therefore, urban and agricultural runoff probably have contributed little in the way of harmful metals or organic compounds to the Refuge, although the potential for such contribution remains from periodic flooding of the ponds by ditch stormflows. Salinity was low throughout most Refuge waters, qualifying as fresh to slightly brackish and suitable for the environmental needs of Refuge fauna. Higher salinities have been observed in ditches during past periods of sugarcane cultivation and saltwater aquaculture, however. Resumption of saltwater aquaculture could raise ditch salinities if saltwater effluents are disposed directly into the ditches, as they were in the past.

Hawaii↗