SEARCH · Search USGS

Results for “Computer contribution”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

465 records · Page 26Linked to original sources

USGS2026

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of Virginia’s Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis conducted on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore shows that groundwater conditions vary substantially throughout the Coastal Plain aquifer system due to local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have fundamentally altered groundwater flow from pre-development conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers supplying groundwater is a relatively small portion of the total groundwater water budget for the aquifer system. Analyses of groundwater budgets from the Virginia Coastal Plain model show that groundwater flow is generally outward from the surficial aquifer to rivers and coastal water bodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer currently is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from areas outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2025, rates of withdrawal substantially lower than the historical maximum have resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or in adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. Downward groundwater flow continues under current conditions, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, those rates are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Estimated groundwater inflow from coastal regions could be contributing to saltwater intrusion, though that was not measured directly in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, show that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up over 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. Currently, estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The current withdrawal rate from the confined Yorktown-Eastover system is near the highest reported rate for the Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and small rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

EarthArXiv↗

USGS2026

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of the Virginia Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis based on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore indicates that groundwater conditions vary substantially throughout the Coastal Plain aquifer system because of local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have altered groundwater flow from predevelopment conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers is a relatively small part of the total groundwater budget for the aquifer system compared to the rate of recharge at the land surface. Analyses of groundwater budgets from the Virginia Coastal Plain model indicate that groundwater flow is generally outward from the surficial aquifer to rivers and coastal waterbodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2023, rates of withdrawal substantially lower than the historical maximum resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. For the most recent simulated conditions (2023) downward groundwater flow continues, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, rates of downward flow are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Groundwater inflow from coastal regions could be contributing to saltwater intrusion, even though that was not measured in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, indicate that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up more than 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. For the most recent simulated period (2023), estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The reported withdrawal rate for this period from the confined Yorktown-Eastover system is near the highest reported rate for the Virginia Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and low rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

Virginia↗

USGS2014

Simulation of the Lower Walker River Basin hydrologic system, west-central Nevada, using PRMS and MODFLOW models

Walker Lake is a terminal lake in west-central Nevada with almost all outflow occurring through evaporation. Diversions from Walker River since the early 1900s have contributed to a substantial reduction in flow entering Walker Lake. As a result, the lake is receding, and salt concentrations have increased to a level in which Oncorhynchus clarkii henshawi (Lahontan Cutthroat trout) are no longer present, and the lake ecosystem is threatened. Consequently, there is a concerted effort to restore the Walker Lake ecosystem and fishery to a level that is more sustainable. However, Walker Lake is interlinked with the lower Walker River and adjacent groundwater system which makes it difficult to understand the full effect of upstream water-management actions on the overall hydrologic system including the lake level, volume, and dissolved-solids concentrations of Walker Lake. To understand the effects of water-management actions on the lower Walker River Basin hydrologic system, a watershed model and groundwater flow model have been developed by the U.S. Geological Survey in cooperation with the Bureau of Reclamation and the National Fish and Wildlife Foundation. The watershed model was developed using the precipitation runoff modeling system (PRMS) and the groundwater flow model was constructed using the MODular groundwater FLOW model (MODFLOW) and both were calibrated for the lower Walker River Basin. These models can be incorporated in an integrated Groundwater and Surface-water FLOW (GSFLOW) model of the lower Walker River Basin. Additionally, the MODFLOW model developed for this study is useful for efficiently simulating long-term and large-scale effects of water-management actions on groundwater hydrology, streamflow, and Walker Lake level, volume, and dissolved-solids concentrations. The lower Walker River Basin PRMS model (LWR_PRMS) was constructed using a subbasin approach to aid in development and calibration, and simulates a 30-year period from 1978 to 2007 using daily time steps. The LWR_PRMS was used to estimate the distribution of groundwater recharge specified in the MODFLOW model. The highest rates of groundwater recharge occur in the Wassuk Range beneath perennial and ephemeral stream channels, whereas lower rates of recharge occur beneath alluvial fans along mountain fronts. The total groundwater recharge estimated using PRMS was about 25,000 acre-feet per year. The lower Walker River Basin MODFLOW (LWR_MF) model simulates an 89-year period using monthly time steps. The LWR_MF was constructed with an initial steady-state simulation to represent dynamic equilibrium conditions from 1908 to 1918 and then a transient simulation representing the period 1919–2007. The model was calibrated using a combination of manual and automated methods of adjusting model parameters to minimize errors between model simulated results and weighted observations of groundwater levels, streamflows, and lake level. Hydrologic conditions simulated with the LWR_MF include the movement and change in storage of groundwater, and the water budgets for Walker River, Walker Lake, and the groundwater system. The LWR_MF computed dissolved-solids concentrations for Walker Lake using simulated lake volume and an assumed constant internal salt mass of 37.2 million tons. Effects of potential changes in water management on future conditions (scenarios) of the lower Walker River Basin hydrologic system and Walker Lake from 2011 to 2070 were evaluated. Several water-management scenarios were considered, including a baseline scenario that represents no changes in system management, improved irrigation efficiencies for the Walker River Indian Irrigation Project (WRIIP), a range of increased streamflows entering the lower Walker River Basin, and, the fallowing of fields on the WRIIP. For the baseline scenario, it was assumed that streamflow conditions from 1981 to 2010 will be repeated in the future. Results indicate that Walker Lake level and volume continue to decline but at a slower rate as the surface area of the lake becomes smaller and lake evaporation decreases. Dissolved-solids concentrations in Walker Lake continue to increase and increase much more rapidly during periods when minimal flows reach the lake due to a diminished lake volume. Alternatively, in years with high runoff, lake level increases are greater and dissolved-solids decreases are greater, compared with equivalent runoffs experienced during 1981–2010. The simulated effects of improving WRIIP efficiencies on Walker River streamflows, Walker Lake inflow, level, and dissolved-solids concentrations, and crop consumptive use, are compared with the baseline reference scenario for a range of irrigation efficiency improvements from 0 to 25 percent over 60 years. Results indicate that water is conserved through a reduction in irrigation-induced groundwater recharge and subsequent groundwater discharge through evapotranspiration. The conserved water mostly goes to increased streamflow to Walker Lake, followed by increased crop consumptive use, then increased evaporation from Weber Reservoir. The simulated effects of increased streamflows at Walker River at Wabuska streamgage (10301500) on Walker Lake inflow, level, and dissolved-solids concentrations, and crop consumptive use, are compared with the baseline scenario after 60 years under two different management methods for Weber Reservoir. Results indicate Walker Lake level and dissolved-solids concentrations stabilized with increased irrigation-season streamflow of about 40,000 acre-feet per year at the Walker River at Wabuska streamgage. Walker Lake level increased, and dissolved-solids concentration decreased, with increased flows of 50,000 acre-feet per year or more. After 60 years with additional irrigation-season streamflows of 50,000 acre-feet per year, Walker Lake level increased by about 48 feet, and lake dissolved-solids concentrations decreased by about 3,000 milligrams per liter (mg/L). With 75,000 acre-feet per year of additional streamflow, Walker Lake level increased by 70 feet, and dissolved-solids concentration decreased by 7,600 milligrams per liter. The effects of fallowing of Walker River Indian Irrigation Project fields from 2007 to 2010 on Walker Lake inflow, level, and dissolved solids were evaluated. Fallowing resulted in a near doubling of Walker River inflow to Walker Lake during this period, an increase in Walker Lake level of about 1.4 feet, and a decrease in dissolved-solids concentration of about 540 mg/L.

Nevada↗

USGS2024

U.S. Geological Survey Earthquake Hazards Program decadal science strategy, 2024–33

Executive Summary Earthquakes represent one of our Nation’s most significant and costly natural hazards, with estimated annual loses from earthquakes close to $15 billion in 2023. Over the past two centuries, 37 U.S. States have experienced an earthquake exceeding a magnitude of 5, and 50 percent of States have a significant potential for future damaging shaking; these statistics speak to the need for nationwide interest and investment in earthquake hazard characterization and risk reduction. Authorized under the Earthquake Hazards Reduction Authorization Act, the U.S. Geological Survey (USGS) Earthquake Hazards Program (EHP) provides the scientific information, situational awareness, and knowledge necessary to reduce deaths, injuries, and economic losses from earthquakes and earthquake-induced tsunamis, landslides, and soil liquefaction. The EHP supports activities in three focused topical areas: (1) earthquake monitoring, (2) hazard assessment, and (3) applied research, using the results of each—and the coordination among them—to further support risk translation and communication in regions at risk nationwide. For earthquake monitoring, the Advanced National Seismic System (ANSS), a cooperative effort of USGS networks, university partner regional seismic networks, and real-time geodetic networks, collects and analyzes data on earthquakes; issues timely, reliable notifications of their occurrence and impacts; and provides data for earthquake research, hazard, and risk assessment as a foundation for building an earthquake-resilient Nation. The USGS-operated ShakeAlert Earthquake Early Warning system is a recent addition to EHP’s ANSS infrastructure. In the realm of earthquake hazard assessment, the EHP contributes to earthquake risk mitigation strategies by developing the National Seismic Hazard Model and maps, and other related products, that describe the likelihood and potential effects of earthquakes nationwide, especially in the urban areas of highest risk. The EHP also conducts research on the causes, characteristics, and effects of earthquakes and prioritizes work that directly increases the accuracy and precision of earthquake hazards assessments, earthquake forecasts, and earthquake monitoring and situational-awareness products and that supports the Nation’s earthquake mitigation practices. Bridging the EHP’s efforts across research, hazard assessments, and earthquake monitoring is a broad and comprehensive collection of earthquake information products, including the National Seismic Hazard Model, ShakeAlert, and other products describing impact, such as ShakeMap and PAGER (Prompt Assessment of Global Earthquakes for Response), which have been developed and integrated into EHP’s real-time monitoring systems. EHP funds external partners to carry out many important collaborative activities through an active external grants program—one of the largest in the USGS—and through cooperative agreements with other partners such as the university-operated regional seismic networks, funded as part of the ANSS. To continue its support of earthquake hazard characterization and risk reduction, the EHP aims to strengthen its foundational products and practices while positioning itself to respond to the evolving needs of the Nation and follow best practices of the scientific community. This document describes a strategy for the program to ensure it can meet these demands. The foundational priorities outlined in this strategy represent those activities that remain critical to the core functionality of the program and those that can be supported under current fiscal year 2024-level appropriations. Priorities described as aspirational are important for future growth, and to maintain the program’s position as a leading global resource in earthquake science, but would require increases in appropriated funding to be fully realized. Across the program’s portfolio of activities, several major themes have been identified as the most critical activities to advance EHP science over the coming decade. Together, these activities provide the framework necessary to integrate critical hazard characterization and risk reduction activities across the program. They provide the structure for research to advance the understanding of where, when, and why earthquakes occur and how we can use improved knowledge to drive short-term and actionable forecasts of seismic activity. They expand the usefulness of critical earthquake products and advance the sophistication of those products to keep pace with the rapidly evolving needs of an ever-expanding user base while maintaining the position of the USGS as a global leader in earthquake science. Focus on system-level science. Establish an automated earthquake-processing pipeline. Enhance the accuracy and reliability of the ShakeAlert Earthquake Early Warning system and plan for extension to other regions. Implement time-dependent earthquake forecasting. Develop physically realistic models. Expand computational capacity. This science strategy is organized into three primary sections. The first section provides an overview of the EHP and its budget, governance, and program council. Readers familiar with the program may wish to focus on the second section, which describes the core of the science strategy, including priorities across each of the EHP’s major program activities in monitoring, hazard assessment, and targeted research. The third section outlines science priorities that cut across program activities, including those involving collaborations external to the EHP.

American Samoa, Guam, Northern Marianas Islands, U↗

USGS1997

Influence of evaporation, ground water, and uncertainty in the hydrologic budget of Lake Lucerne, a seepage lake in Polk County, Florida

Evaporation losses and the interaction of ground water with Lake Lucerne were studied to determine the influence of these two processes on the hydrologic budget of a seepage lake. Lake Lucerne is representative of the numerous seepage lakes of sinkhole origin in the karst terrain of central Florida. Because of permeable surficial deposits, ground-water inflow is the only significant contribution from the surrounding watershed. The lake recharges the underlying Upper Floridan aquifer and, as a result, is susceptible to increased leakage induced by pumping from this aquifer. Ground-water fluxes determined in the study were analyzed to define the proportion of the total lake leakage induced by pumping from the Upper Floridan aquifer. A hydrologic budget is analyzed for the 1-year period from October 1985 to September 1986. Ground-water inflow and leakage are significant components of the hydrologic budget. Changes in the quantity of either of these fluxes can substantially alter lake stage. Ground-water inflow contributed from 20 to 37 percent of the total annual inflow to the lake. Leakage from the lake accounted for 18 to 23 percent of the total annual outflow. Water withdrawals from the Upper Floridan aquifer increased annual lake leakage by 22 percent over nonpumping conditions. Most of the increase (92 percent) in leakage occurred during April, May, and June 1986, when local citrus irrigation was highest. For the study year, ground-water inflow and leakage volumes were calculated by flow-net analysis to be equal to 10.5 and 12.6 inches, respectively, of water depth above (or below) the lake surface. These estimates were revised upward on the basis of an analysis of the error in the hydrologic-budget equation. Revised ground water inflow exceeded annual leakage from the lake. Ground-water inflow rates were increased by 120 percent to 23.6 inches, and leakage was increased by 40 percent to 17.5 inches. Differences between the two estimates probably reflect the uncertainty in the hydraulic conductivity estimates of the porous media around the lake and the unaccounted effect of transient ground-water inflow. The geometry of the sinkhole complex beneath Lake Lucerne and pumping in the Upper Floridan aquifer are primary controls on ground water interactions with the lake and, in particular, lake leakage. A numerical ground-water model was used to test the effects of these two factors on ground-water interactions with the lake. Results indicate that the intermediate confining unit below Lake Lucerne has been breached and replaced by materials about two orders of magnitude more conductive. Anisotropy in the surficial aquifer is approximately 100 and controls the depth of the ground-water flow intercepted by the lake. Lake sediments having low permeability may control the distribution of leakage through the lakebed but did not appreciably reduce total leakage rates in these simulations. Evaporation loss was the major outflow component of the hydrologic budget. Annual lake evaporation determined by the energy-budget method was 57.9 inches, about 8 inches greater than long-term estimates for the region. The greater rate was attributed to drier than normal conditions: rainfall totaled 40.9 inches during the study year, about 10 inches less than the long-term average. Similar annual evaporation rates were determined by the energy-budget method, the simpler mass-transfer method, and by corrected pan evaporation from an onsite pan; however, the probable errors associated with these other two methods were greater than for the energy budget. Weekly energy-budget evaporation rates ranged from 0.04 inch per day in early January 1986 to 0.26 inch per day in early May 1986. The largest monthly energy-budget evaporation rates occurred in April and May 1986, 7.16 and 7.12 inches per month, respectively. Monthly evaporation estimated from corrected pan evaporation generally was within 10 percent of the energy-budget estimate but differed by as much as 35 percent. Daily energy-budget evaporation also was computed for 321 days during the year.

Florida↗

USGS1947

Age and growth of the kiyi, Leucichthys kiyi Koelz, in Lake Michigan

Ages were determined and individual growth histories were calculated from the examination and measurement of the scales of 1,649 kiyis captured at seven localities in Lake Michigan in 1931 and 1932. The numbers of individuals employed for the investigation of other phases of the life history (such as length-frequency distributions, length-weight relationship, and sex ratio) varied according to the amount of data available or required. Age-group IV was dominant in the 1931 collections from Racine, Port Washington, and Kewaunee, Wisconsin, and age-group V dominated the 1932 samples from the Fox Islands and from three localities southward of Manistique, Michigan. A trend was noticeable toward an increase in average age from south to north. Among the explanations suggested for the observed differences in age composition were: Variation with latitude in the natural span of life; differences in fishing intensity; fluctuations in the strength of year classes (to account possibly for the shift in the dominant age group from 1931 to 1932). The oldest male kiyi belonged to the VII group and the oldest female was a member of the X group. The possible distorting effects of such factors as gear selection traceable to differences in the mesh sizes of nets fished in 1931 and 1932, selection by nets on the basis of the condition (K) of the fish, and local variations in fishing intensity and hence in the selective destruction of rapidly growing individuals in the fishery were held to be sufficiently great to render doubtful the significance of most of the observed local differences in growth rate. Kiyis from all samples were combined to determine the general growth in length. The growth in weight of the Fox Islands fish, however, was considered separately as these fish were consistently lighter than kiyis of corresponding length from other localities. The Lake Michigan kiyi grows slowly, with the females growing slightly more rapidly than the males. The grand average calculated lengths indicated, for example, that the females did not attain a total length of 10 inches until the fifth year of life or the males until the sixth. Similarly, the calculated weight of 4 ounces was not reached until the fifth or sixth year (with the actual time varying with sex and locality). The season's growth of the kiyi probably begins sometime in May and most or all of the growth is completed by the end of August. The calculated lengths of the age groups exhibited large discrepancies that differed from “Lee's phenomenon” as ordinarily observed in that the data for the later rather than the earlier years of life were affected most severely. Chief among the factors held responsible for these discrepancies were gear selection and the selective destruction of the more rapidly growing individuals in the fishery. Errors inherent in the (direct-proportion) method of computing growth from scale measurements were considered to have been unimportant. The Lake Michigan kiyi exhibits growth compensation–the tendency for the smaller of the young fish to have the more rapid growth in the later years of life. Comparisons with the average lengths and weights of the age groups of the Lake Ontario kiyi given by Pritchard (1931) indicated the Lake Michigan fish to be the larger at the earlier ages (age-groups II and III) and the smaller at the later ages (age-groups IV to VI). The length-frequency distributions of the age groups exhibited extensive overlap. As many as eight age groups were represented in a single centimeter interval of length. The length frequencies and average lengths of all fish collected, arranged according to the mesh sizes of the gill nets by which they were captured, revealed that the selective action of these nets in the taking of kiyis was much more obvious in the numbers of fish in the catch than in their average size. As an illustration, in 1930–1931, the 2 3/4-inch mesh nets took fish that were only 0.1 inch longer than those in 2 1/2-inch meshes but captured less than one fourth as many. Gill nets fished in northern Lake Michigan in 1932 captured kiyis that averaged 0.2 to 0.4 inch longer than those taken in the same meshes in southern Lake Michigan in 1930–1931. Because of the more slender form of kiyis from the northeastern island region of Lake Michigan, data on the general length weight relationship were compiled separately for fish of that area and for those of the great central basins of the lake. In both regions the weight increased to a power slightly greater than the cube of the length. Available information on condition indicated that the coefficient (K) was higher in August and early September than in May, June, and July. Condition declined from early September to October and early November–the latter period the time of most active spawning. Spawning itself was accompanied by an additional loss of about 12 per cent of the body weight of females and of somewhat less than 2 per cent of the weight of males. Analysis of the variations of K within a group that was homogeneous with respect to age, sex, maturity, and time of collection revealed that a net of a particular mesh size tends to take the heavier of the shorter fish and the lighter of the longer fish within its range of effectiveness. Among fish of the same length the values of K tended to increase with increase in the mesh size of the nets employed for their capture. Practically all fish in the samples were mature (only 11 immature in more than 6,000). These “immature” fish were probably “non-functional” since all of them approached or exceeded the average length of the mature kiyis. Females were strongly predominant in the collections at all seasons but were relatively more plentiful during the summer (90 per cent of the total) than during the spawning period (75 per cent). Possible factors contributing to this predominance of females and to the change in the sex ratio at the spawning season were discussed. A decrease in the relative abundance of males with increase in age appears to be characteristic of the kiyi. This decrease indicates a differential mortality of th sexes (greater relative destruction of males in the spawning period when they are unusually abundant or a greater natural mortality rate for the males). Current fishery regulations on mesh size and closed seasons afford the kiyi good protection but offer no guarantee against depletion from too intensive fishing.

Transactions of the American Fisheries Society↗

USGS2004

Geohydrology of the French Creek Basin and simulated effects of drought and ground-water withdrawals, Chester County, Pennsylvania

This report describes the results of a study by the U.S. Geological Survey, in cooperation with the Delaware River Basin Commission, to develop a regional ground-water-flow model of the French Creek Basin in Chester County, Pa. The model was used to assist water-resource managers by illustrating the interconnection between ground-water and surface-water systems. The 70.7-square mile French Creek Basin is in the Piedmont Physiographic Province and is underlain by crystalline and sedimentary fractured-rock aquifers. Annual water budgets were calculated for 1969-2001 for the French Creek Basin upstream of streamflow-measurement station French Creek near Phoenixville (01472157). Average annual precipitation was 46.28 in. (inches), average annual streamflow was 20.29 in., average annual base flow determined by hydrograph separation was 12.42 in., and estimated average annual ET (evapotranspiration) was 26.10 in. Estimated average annual recharge was 14.32 in. and is equal to 31 percent of the average annual precipitation. Base flow made up an average of 61 percent of streamflow. Ground-water flow in the French Creek Basin was simulated using the finite-difference MODFLOW-96 computer program. The model structure is based on a simplified two-dimensional conceptualization of the ground-water-flow system. The modeled area was extended outside the French Creek Basin to natural hydrologic boundaries; the modeled area includes 40 square miles of adjacent areas outside the basin. The hydraulic conductivity for each geologic unit was calculated from reported specific-capacity data determined from aquifer tests and was adjusted during model calibration. The model was calibrated for above-average conditions by simulating base-flow and water-level measurements made on May 1, 2001, using a recharge rate of 20 in/yr (inches per year). The model was calibrated for below-average conditions by simulating base-flow and water-level measurements made on September 11 and 17, 2001, using a recharge rate of 6.2 in/yr. Average conditions were simulated by adjusting the recharge rate until simulated streamflow at streamflow-measurement station 01472157 matched the long-term (1968-2001) average base flow of 54.1 cubic feet per second. The recharge rate used for average conditions was 15.7 in/yr. The effect of drought in the French Creek Basin was simulated using a drought year recharge rate of 8 in/yr for 3 months. After 3 months of drought, the simulated streamflow of French Creek at streamflow-measurement station 01472157 decreased 34 percent. The simulations show that after 6 months of average recharge (15.7 in/yr) following drought, streamflow and water levels recovered almost to pre-drought conditions. The effect of increased ground-water withdrawals on stream base flow in the South Branch French Creek Subbasin was simulated under average and drought conditions with pumping rates equal to 50, 75, and 100 percent of the Delaware River Basin Commission Ground Water Protected Area (GWPA) withdrawal limit (1,393 million gallons per year) with all pumped water removed from the basin. For average recharge conditions, the simulated streamflow of South Branch French Creek at the mouth decreased 18, 28, and 37 percent at a withdrawal rate equal to 50, 75, and 100 percent of the GWPA limit, respectively. After 3 months of drought recharge conditions, the simulated streamflow of South Branch French Creek at the mouth decreased 27, 40, and 52 percent at a withdrawal rate equal to 50, 75, and 100 percent of the GWPA limit, respectively. The effect of well location on base flow, water levels, and the sources of water to the well was simulated by locating a hypothetical well pumping 200 gallons per minute in different places in the Beaver Run Subbasin with all pumped water removed from the basin. The smallest reduction in the base flow of Beaver Run was from a well on the drainage divide between the French Creek Basin and the Marsh Creek Basin to the south; the simulated base flow of Beaver Run at the mouth was reduced 1 percent. The greatest reduction in the base flow of Beaver Creek was from a well close to Beaver Run; the simulated base flow of Beaver Run at the mouth was reduced 8 percent. The simulations showed that (1) if the contributing area of a well is in a basin, pumping will affect stream base flow and water levels in that basin whether the well is inside or outside that basin; (2) wells in different areas of a basin away from a divide produce a similar reduction in base flow; (3) a well within a basin will derive more water from diverted base flow and less water from storage than a well on or near a basin divide; and (4) the reduction in base flow at the mouth of the stream is the same for a well in the headwaters and a well downstream near the confluence. Model simulations illustrate some of the typical analyses and results that can be produced. The model was calibrated using annual values for recharge and ground-water ET and then was run using the annual values in a seasonally independent transient mode to show changes with time. The timing and relative magnitude of some of the changes simulated with the model when viewed in terms of a normal climatic year may be subject to considerable uncertainty because of the variability in seasonal recharge and ground-water ET rates. Transient model simulations for short-term periods are indicative of possible hydrologic system response and are considered an approximation.

Pennsyvlania↗

USGS2024

Status and understanding of groundwater quality in the Mojave Basin Domestic-Supply Aquifer study unit, 2018—California GAMA Priority Basin Project

Groundwater quality in the western part of the Mojave Desert in San Bernardino County, California, was investigated in 2018 as part of the California State Water Resources Control Board Groundwater Ambient Monitoring and Assessment Program Priority Basin Project. The Mojave Basin Domestic-Supply Aquifer study unit (MOBS) region was divided into two study areas—floodplain and regional—to assess differences between the two major aquifers used for drinking water supply in the area. This assessment characterized the quality of ambient groundwater and not the quality of treated drinking water. The study included three components: (1) a status assessment, which characterized the quality of groundwater resources used for domestic drinking-water supply in the floodplain and regional study areas; (2) a brief understanding assessment, which evaluated factors that could potentially affect the quality of groundwater used by domestic wells in the region; and (3) a comparative assessment between the groundwater resources used by domestic wells and public-supply wells in the two study areas. The domestic-well assessment was based on data collected by the U.S. Geological Survey from 48 domestic wells in January–May 2018. The public-supply assessment was based on data for samples from 322 public-supply wells in 2008–18, either collected by the U.S. Geological Survey or compiled from the California State Water Resources Control Boards Division of Drinking Water publicly available database. Concentrations of water-quality constituents in ambient groundwater were compared to regulatory and non-regulatory benchmarks typically used by the State of California and Federal agencies as health-based or aesthetic standards for public drinking water. Relative concentrations, defined as the measured concentration divided by the benchmark concentration, were classified as high (greater than 1.0), moderate (greater than 0.5 for inorganic constituents or 0.1 for organic and special-interest constituents, and not high), or low (concentrations lower than moderate). The floodplain and regional study areas were divided into 15 and 35 grid cells, respectively, and grid-based methods were used to compute the areal proportions of the two study areas with high, moderate, or low relative concentrations of individual constituents and classes of constituents. For the domestic-supply assessment, one or more inorganic constituents with health-based benchmarks were detected at high relative concentrations in 58 percent of the regional study area and 13 percent of the floodplain study area. The inorganic constituents with health-based benchmarks detected at high relative concentrations in the regional study area were arsenic, chromium and hexavalent chromium, fluoride, adjusted gross alpha particle activity, uranium, molybdenum, strontium, and nitrate; only arsenic was detected at high relative concentrations in the floodplain study area. One or more inorganic constituents with secondary maximum contaminant level benchmarks were detected at high concentrations in 15 and 6.7 percent of the regional and floodplain study areas, respectively. The constituents detected at high relative concentrations in the regional study area were total dissolved solids, chloride, sulfate, and iron; only total dissolved solids and sulfate were detected at high relative concentrations in the floodplain study area. Organic constituents were not detected at moderate or high relative concentrations in either the regional or floodplain study areas. Volatile organic compounds were detected at low relative concentrations in 21 and 27 percent of the regional and floodplain study areas, respectively, and pesticides were detected at low relative concentrations in 9.1 and 20 percent of the regional and floodplain study areas, respectively. The only individual organic constituent detected in more than 10 percent of either study area was the trihalomethane trichloromethane. Total coliform bacteria were detected in 15 and 27 percent of the grid wells in the regional and floodplain study areas, respectively. The greater prevalence of high relative concentrations of many inorganic constituents in the regional study area compared to the floodplain area likely indicates the greater diversity of geologic material at depth in aquifer material and generally finer-grained alluvium compared to the floodplain study area combined with generally older groundwater that has had more contact time with aquifer materials. In general, trace element concentrations (1) increased with increasing groundwater age, (2) increased with distance from recharge sources in the mountains, and (3) increased with closer proximity to some types of geological units. In general, groundwater from domestic wells in the floodplain study area is young, with most samples containing a component of modern groundwater based on tritium and unadjusted carbon-14 activities, whereas groundwater from domestic wells in the regional study area generally is old, with most samples having unadjusted carbon-14 ages of 5,000–40,000 years. Public-supply wells in MOBS generally were deeper than domestic wells and presumably are in contact with older, more weathered alluvium that may have more mobile trace elements, such as arsenic or uranium. However, only 26 percent of the public-supply regional study area had high relative concentrations of inorganic constituents, compared to 58 percent for the domestic regional study area. The percentages of the public-supply and domestic floodplain study areas with high relative concentrations of inorganic constituents were 11 and 13 percent, respectively. The ages of groundwater used by public-supply and domestic wells in each study area were similar, which was not expected given the greater depth of the public-supply wells. Three potential factors may contribute to these results: (1) greater spatial footprint of domestic well network, which may result in domestic wells pumping groundwater from fractured bedrock or mineralized areas not used by public-supply wells; (2) greater pumping rates in public-supply wells, resulting in more water being withdrawn from coarse-grained, heterogeneous alluvium than finer-grained layers, which may have higher concentrations of (or more mobile) inorganic constituents; and (3) a greater degree of well management with public-supply wells, which may include pausing use of or decommissioning wells if treating or blending water is not feasible to lower constituent concentrations.

California↗

USGS2007

Mercury in precipitation in Indiana, January 2001–December 2003

Mercury in precipitation was monitored during 2001 through 2003 at four locations in Indiana as part of the National Atmospheric Deposition Program-Mercury Deposition Network (NADP-MDN). Monitoring stations were operated at Roush Lake near Huntington, Clifty Falls State Park near Madison, Monroe County Regional Airport near Bloomington, and Indiana Dunes National Lakeshore near Porter. At these monitoring stations, precipitation amounts were measured and weekly samples were collected for analysis of total mercury and methylmercury by low-level methods. Wet deposition was computed with the total mercury and methylmercury concentrations and the precipitation amounts. In 3 years of weekly samples collected at the four monitoring stations, the volume-weighted total mercury concentration was 11.5 ng/L (nanograms per liter). As a reference for comparison, the total mercury concentration in 47 percent of the samples analyzed was greater than the Indiana water-quality standard for mercury (12 ng/L, protecting aquatic life) and nearly all of the concentrations exceeded the Indiana water-quality standards for mercury in the Great Lakes system (1.8 ng/L, protecting human health, and 1.3 ng/L, protecting wild mammals and birds). The precipitation-weighted concentrations at three of the monitoring stations in Indiana in 2003 were in the top 40 percent of all monitoring stations in the NADP-MDN and the concentration at Indiana Dunes was the eighth highest in the NADP-MDN for 2003. At the four monitoring stations during the study period, the mean weekly total mercury deposition was 243 ng/m 2 (nanograms per square meter) and mean annual total mercury deposition was 12,623 ng/m 2 . The annual mercury deposition at the four monitoring stations in Indiana in 2003 was in the top 40 percent of all monitoring stations in the NADP-MDN and the annual mercury deposition at the Clifty Falls station was the tenth highest in the NADP-MDN for 2003. For the 3-year period, the median methylmercury concentration in weekly samples was 0.058 ng/L with a maximum of 5.77 ng/L. Normalized methylmercury deposition was 2.09 ng/m 2 per inch of precipitation and methylmercury deposition was 0.7 percent of the total mercury deposition. The annual and mean weekly methylmercury deposition was highest at the Roush Lake station. Among the monitoring stations in the NADP-MDN with methylmercury data, methylmercury deposition at the monitoring stations in Indiana appeared to be higher than at eight stations in Wisconsin and Minnesota for that same time period, although methylmercury concentrations in Indiana were similar to or lower than those in Wisconsin and Minnesota. Geographically, the weekly total mercury concentrations at Indiana Dunes and Clifty Falls were statistically higher than concentrations at Bloomington, although a statistical difference in weekly total mercury deposition was not found among the four monitoring stations. Annual mercury emissions from sources in the vicinity of Indiana Dunes and Clifty Falls in 2001 were more than 10 times those at Bloomington, although other factors may help explain the differences in total mercury concentrations, such as the types of mercury emissions, mercury transport from sources outside Indiana, and meteorological conditions. Mercury concentrations and deposition varied at the four monitoring stations during the 3-year period. Total mercury concentrations in weekly samples ranged from 1.54 to 77 ng/L and weekly mercury deposition ranged from 0.8 to 2,456 ng/m 2 . Data from weekly samples exhibited seasonal patterns. Total mercury concentrations and deposition were highest in spring and summer and lowest in winter. Methylmercury concentrations were highest in winter and methylmercury deposition was highest in spring. Annual precipitation at the four monitoring stations was highest in 2003, exceeding the precipitation normals in spring and summer 2003. Annual mercury deposition in 2003 at the Roush Lake, Clifty Falls, and Indiana Dunes was as much as 41 to 67 percent higher in 2003 than in 2001 or in 2002 at those stations. Total mercury deposition that was more than 10 percent of the mean annual deposition (1,262 ng/m 2 ) was recorded in 11 of 551 weekly samples from the study period. These samples contained approximately 3 inches or more of rain and most were collected in spring and summer 2003. The highest deposition (2,456 ng/m 2 in a sample from Roush Lake) was 15.7 percent of the annual deposition at that station and approximately 10 times the mean weekly deposition for Indiana. High deposition recorded in three weekly samples at Clifty Falls contributed 31 percent of the annual deposition at that station in 2003. Weekly samples with high mercury deposition may help to explain the differences in annual mercury deposition among the four monitoring stations in Indiana.

Indiana↗

USGS2022

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↗

USGS2009

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↗

USGS2021

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↗

USGS2012

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↗

USGS2013

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↗

USGS2000

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↗