USGS Science⌕ Search

SEARCH · USGS Science

Results for “Computers & Geosciences”

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

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

97 records · Page 6Linked to original sources
USGS2024

Effective site coefficients for the 2024 International Building Code (IBC)

The U.S. National Seismic Hazard Models (NSHMs), developed by the U.S. Geological Survey (USGS), have long been the scientific foundation of seismic design guidelines and have been used to compute design ground motions for construction of new buildings and retrofit of existing buildings in the United States and its territories. The 2018 NSHM is adopted by the 2024 International Building Code (IBC). Prior to the 2018 NSHM update, hazard calculations were required only at one reference site condition defined by V S30 =760 m/s (where V S30 is the time-averaged shear wave velocity from the surface to a depth of 30 m) and three periods (peak ground acceleration, PGA, and pseudo spectral accelerations at periods of 0.2 and 1 s, S s and S 1 ). Site coefficients, F PGA , F a , and F v , were then defined by the Building Seismic Safety Council (BSSC) Provisions Update Committee (PUC) in the site-specific procedures of National Earthquake Hazard Reduction Program (NEHRP) Recommended Seismic Provisions to calculate ground motions for other site classes with different V S30 values at the given periods. The design ground motions at other periods were then estimated using a generic spectral shape that was also defined by the BSSC PUC in NEHRP provisions. In recent years, the engineering community has realized there were deficiencies with the F a and F v site coefficients and the generic spectral shape. To avoid potentially dangerous underestimations of design ground motions for long period structures on soft site conditions in high seismicity regions, the BSSC PUC recommended the use of multi-period response spectra (MPRS) in 2017. As a result, the USGS produced multi-period and multi-V S30 response spectral values in the 2018 NSHM for calculations of design ground motions and the site coefficients F a and F v were eliminated from the 2020 NEHRP Provisions. As these site coefficients were widely used inside and outside of the United States, in this study we back-calculate the “effective” site coefficients F a,eff , and F v,eff by comparing MPRS for various site classes with the MPRS for the reference site condition, and discuss the changes that are observed in the 2024 IBC compared to its previous version in 2021. The effective site coefficients are presented for test site locations and their dependence on various factors including period, ground motion intensity, and regional models are discussed. Ratio maps between the new effective site coefficients and the old ones are then presented for soft site classes and for short and long periods. For soft site classes at short periods, the new effective site coefficients are lower than the old site coefficients for high seismicity regions and higher for low seismicity regions. As it was expected, for soft site classes at long periods and high seismicity regions, the new effective site coefficients are much larger than the old site coefficients without imposing the 50% increase as a penalty that was suggested in the 2021 IBC, whereas they could be much smaller if the 50% increase would have been imposed particularly around New Madrid and Charleston high seismicity regions. For low seismicity regions, the long period effective site coefficients can be smaller or larger by 20% compared to the 2021 IBC coefficients.

Conference Paper↗
USGS2025

Phase 1 technical implementation plan for the expansion of the ShakeAlert earthquake early warning system to Alaska

Executive Summary The conference report accompanying the fiscal year (FY) 2022 Consolidated Appropriations Act (Public Law 117–103) for the U.S. Department of the Interior and related agencies directed the U.S. Geological Survey (USGS) to “work with the State of Alaska to develop an implementation plan to be completed within two years in order to put ShakeAlert/Earthquake Early Warning in Alaska” (p. 29). Congress included $1 million in the FY 2022 appropriation to conduct this effort. The USGS Earthquake Hazards Program, along with partner organizations, has developed the ShakeAlert earthquake early warning (EEW) system for the West Coast, which currently operates in California, Oregon, and Washington. The purpose of the system and its alert delivery partners is to reduce the impact of earthquakes and save lives and property by delivering ShakeAlert-powered alerts that are transmitted to the public via mass notification technologies, and by providing more detailed data streams to institutional users and commercial service providers to trigger automated, user-specific, protective actions. ShakeAlert was designed in such a way that it could be expanded to other U.S. regions with high earthquake risk, after the build-out of seismic and geodetic networks to support ShakeAlert in a specified region is completed and the necessary funding is secured for long-term operation and maintenance. When an earthquake occurs, seismic waves radiate from the rupturing fault like waves on a pond. It is these waves that people feel as earthquake shaking and that can cause damage to structures. Using networks of ground-motion sensors and sophisticated computer algorithms, ShakeAlert can detect an earthquake seconds after it begins, calculate its location and magnitude, and estimate the resulting intensity of shaking. Early warnings of impending shaking are then sent to people and systems that may experience damaging shaking, allowing them to take appropriate protective actions. Depending on the user’s distance from the earthquake, alerts may be delivered before, during, or after the arrival of strong shaking. There will almost always be a region near the earthquake epicenter where alerts arrive after damaging shaking has begun. The ShakeAlert system updates its ground-motion estimates as an earthquake grows larger. In response to the FY 2022 congressional direction, the USGS worked with the State of Alaska to devise this implementation plan for ShakeAlert expansion to Alaska. The USGS engaged with the Alaska Division of Homeland Security and Emergency Management (DHS&EM) and the Alaska Division of Geological and Geophysical Surveys (DGGS). A cooperative agreement was awarded to the Alaska Earthquake Center (AEC) at the University of Alaska Fairbanks (UAF) for their contributions to the plan and their work coordinating with other networks in Alaska. The USGS engaged with the Alaska Seismic Hazards Safety Commission (ASHSC) throughout the process. The USGS also held a series of Alaska stakeholder engagements. The process of developing the implementation plan was facilitated by contracted staff from Corner Alliance, which is a government consulting firm. This implementation plan describes the details and estimates the costs for a Phase 1 expansion of the ShakeAlert system to Alaska. A geographically limited Phase 1 goal was chosen that covers the highest risk and most populated areas of Alaska. The areas proposed encompass the State’s main population centers and 90 percent of the State’s population. This Phase 1 design is considered very challenging and ambitious from the viewpoint of network operators. The lessons learned if this plan is implemented could be used to consider subsequent phases to expand EEW beyond Phase 1 in Alaska in the future. ShakeAlert is built on the foundation of the sensor networks and data processing infrastructure of the USGS-led Advanced National Seismic System (ANSS). This implementation plan calls for a total of 450 high-quality, real-time EEW-capable ANSS seismic stations in Alaska: 270 new stations, 160 upgraded stations, and 20 existing stations. These seismic station numbers are based on a station spacing of 10 kilometers (km) in urban areas, 20 km in seismic source areas that endanger population centers, and 40 km in other areas. The associated costs also include support for some EEW-capable global navigation satellite system (GNSS) stations, with a focus on improving warnings for large subduction zone earthquakes. For effective EEW, ShakeAlert requires low-latency, high-availability, robust telemetry links to deliver continuous, real-time data from field stations to the data centers. The Alaska data processing hardware infrastructure would follow the general design for fail-safe operation that is used for the ShakeAlert system on the West Coast. The ShakeAlert architecture uses two independent layers: the production layer for earthquake processing and the alert layer to make alerting decisions and serve alerts to users. This implementation plan includes two geographically separated data centers in Alaska, each with two fully independent production and alert layers using the same system design developed for the West Coast. As of March 2024, the ShakeAlert system is at version 3.0.1, with more advanced versions in the development and testing pipeline. ShakeAlert originally used two algorithms to determine the location and magnitude of earthquakes using seismic data. A third algorithm that can calculate very large magnitudes of very large earthquakes with geodetic data was added in March 2024. ShakeAlert publishes several data and alert products to meet the needs of different users. All messages include the location of the earthquake, either as a point or a line, and its magnitude. Ground-shaking estimates are published in two forms, as ground-motion contours and a map grid. Providing adequate warning time for strong shaking (the “target threshold”) requires sending alerts at a threshold lower than that strong shaking level (the “alert threshold”). The thresholds for public alerting in Alaska would be a joint USGS and State decision. To have the greatest benefit, ShakeAlert-powered alerts would be delivered to institutional users and individuals by all practical pathways. The USGS alert layer can support thousands of institutional users and alert redistributors, but the USGS does not have the mission nor the infrastructure and expertise to perform mass notifications to the public or implement automatic actions for end users of the alerts. To meet this need, ShakeAlert recruits private sector “technology enablers” that have the necessary expertise to develop end-user implementations using EEW alerts with the goal of stimulating an EEW industry. Earthquake early warning alerts are useless if people do not know how to respond to them. Although the alert messages include instructions about what to do (drop, cover, and hold on), alerts are more effective if people have been trained in advance. Messages about ShakeAlert’s capabilities, limitations, and benefits could be integrated with existing earthquake education programs, including State-run programs. Therefore, ShakeAlert would coordinate with both public and private partners and stakeholders through various partnerships and agreements to accomplish consistent and ongoing public earthquake hazard education. The estimated capital cost of completing the computing infrastructure and sensor networks for the Phase 1 ShakeAlert expansion to Alaska is approximately $66 million in 2024 dollars. The annual operation and maintenance cost of the completed system is estimated to be $12 million per year in 2024 dollars when fully built out.

Alaska↗
USGS1999

Comparison of two approaches for determining ground-water discharge and pumpage in the lower Arkansas River Basin, Colorado, 1997-98

In March 1994, the Colorado Division of Water Resources (CDWR) adopted “Rules Governing the Measurement of Tributary Ground Water Diversions Located in the Arkansas River Basin” (Office of the State Engineer, 1994); these initial rules were amended in February 1996 (Office of the State Engineer, 1996). The amended rules require users of wells that divert tributary ground water to annually report the water pumped monthly by each well. The rules allow a well owner to report the pumpage measured by a totalizing flowmeter (TFM) or pumpage determined from electrical power data and a power conversion coefficient (PCC) (Hurr and Litke, 1989). Opinions by representatives of the State of Kansas, presented before the Special Master hearing a court case [State of Kansas v. State of Colorado, No. 105 Original (1996)] concerning post-Compact well pumping, stated that the PCC approach does not provide the same level of accuracy and reliability as a TFM when used to determine pumpage. In 1997, the U.S. Geological Survey (USGS), in cooperation with the CDWR, began a 2-year study to compare ground-water pumpage estimates made using the TFM and the PCC approaches. The study area was along the Arkansas River between Pueblo, Colorado, and the Colorado-Kansas State line (fig. 1). The two approaches for estimating ground-water discharge and pumpage were compared for more than 100 wells completed in the alluvial aquifer of the Arkansas River Basin. The TFM approach uses an inline flowmeter to directly measure instantaneous discharge and the total volume of water pumped at a well. The PCC approach uses electrical power consumption records and a power conversion coefficient to estimate the pumpage at ground-water wells. This executive summary describes the results of the comparison of the two approaches. Specifically, (1) the differences in instantaneous discharge measured with three portable flowmeters and measured with an inline TFM are evaluated, and the statistical differences in paired instantaneous discharge between the two approaches are determined; (2) short- and long-term variations in the PCC’s are presented; (3) differences in pumpage between the two approaches are evaluated, and the statistical differences in pumpage between the two approaches are determined; (4) potential sources of discrepancy between pumpage estimates are discussed; and (5) differences in total network pumpage using the two approaches are presented. During the irrigation seasons of 1997 and 1998, instantaneous discharge and electrical power demand were measured at randomly selected wells to determine PCC’s. At more than 100 wells, the PCC’s determined during the 1998 season were applied to total electrical power consumption data that was recorded between the initial and final readings at each network well site in 1998 to estimate total ground-water pumpage. At each site, an inline TFM was installed in a full-flowing, acceptable test section of pipe on the discharge side of the pump where the measurement of discharge was made. Measurements of instantaneous ground-water discharge also were made using three different types of portable flowmeters. The average velocity multiplied by the cross-sectional area of the discharge pipe was used to compute the discharge in gallons per minute. Whenever possible, discharge measurements were made at each network site using all three types of portable flowmeters.

Colorado↗
USGS2000

Regional evaluation of evapotranspiration in the Everglades

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

Florida↗
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↗
USGS2026

Water-withdrawal and groundwater-level fluctuations, physical and infiltration characteristics of Ohio River sediment, and isotope-based estimates of water sources relative to a riverbank filtration system in an alluvial aquifer, Louisville, Kentucky, 2000–20

The U.S. Geological Survey and the Louisville Water Company, from 2000 to 2020, cooperatively investigated physical and hydrologic conditions in the Ohio River alluvial aquifer and adjacent riverbed sediments in the northeast part of Louisville, Kentucky, to increase knowledge of processes affecting water sources to collector wells during operation of a riverbank filtration system (RBF). The RBF, in 2020, included an initial collector well (CW1) completed in 1999 and a tunnel and collector well system (TCWS) composed of four collector wells (CW2–CW5) completed in late 2010 that were connected by a tunnel that routed water to a pumping station and a treatment plant. Daily mean water-withdrawal rates typically ranged from about 10.4 to 20.7 million gallons per day (Mgal/d) from 2005 to 2007, then were cyclic from 2008 to 2010 and ranged from as low as about 9 to 13 Mgal/d in winter-spring months up to about 19 to 20 Mgal/d in summer-autumn months. From 2011 through 2017, peak daily mean water-withdrawal rates from CW2–CW5 fluctuated between winter-spring lows of about 20–30 Mgal/d to summer-autumn peak rates of about 41.6 Mgal/d in 2013 to about 48.2 Mgal/d in 2012. From 2017 to 2020, TCWS peak daily mean water-withdrawal rates decreased from about 45.7 Mgal/d in 2017 to about 37.5 Mgal/d in 2020 and peak daily mean rates from CW1 ranged from about 20.7 Mgal/d in 2020 to about 21.6 Mgal/d in 2018. Periods in 2019–20 with tunnel pumping station daily mean water-withdrawal rates smaller than about 30 Mgal/d had generally stable water levels in CW2 and indicated the sustainable yield of the TCWS collector wells under the prevailing hydrologic conditions. Drawdown during July–August of 2008 and 2009 was associated with CW1, the sole operating collector well. Groundwater-flow directions outside the CW1 cone of depression during July–August of 2008 and 2009 were from areas near the alluvial valley margin toward the Ohio River and in 2009 toward Harrods Creek. Collector wells CW2 through CW5 derived water during 2011–19 surveys from surface-water infiltration, as indicated by drawdown that extended beneath the Ohio River and Harrods Creek. Increases in drawdown toward the alluvial valley margins in July–August of 2013 and 2015, compared with July 2011, indicated less available recharge from those flow directions and decreased yield of groundwater from the aquifer after progressively longer TCWS water withdrawals. Drawdown at CW2–CW5 and at observation wells between them lessened in October 2019 and indicated that onshore groundwater levels can recover when collector well withdrawals decrease. Overlapping cones of depression in 2011–15 potentiometric surfaces centered around CW2, CW3, CW4, and CW5, and drawdown along divides between those wells indicated well interference between collector wells. Well-interference effects increased during 2013 and 2015 survey periods, despite decreased water-withdrawal rates compared with July 2011, and expanded to include CW1 when its use resumed. Groundwater-level fluctuations in observation wells near the Ohio River were consistent with gradual depletion of groundwater yield during tunnel pumping station operation at river elevations of about 420 feet above the National Geodetic Vertical Datum of 1929 (NGVD 29) and enhanced recharge of the aquifer when river elevations were above flood stage of about 430 feet above NGVD 29. Groundwater-level declines in observation wells near the Ohio River were most apparent in summer and some autumn months between 2012 and 2017 during periods when tunnel pumping station mean daily water-withdrawal rates were about 35 Mgal/d or more. Periods of relatively stable groundwater levels during 2019–20 in observation wells between CW2–CW5, when tunnel pumping station water-withdrawal rates were about 30–33 Mgal/d or less, also indicated that the prevailing sustainable yield was about 30 Mgal/d or less. Wells associated with CW3 and CW4 had the largest increases in slopes of best-fit lines relating water-level difference to water-withdrawal rate between 2012 and 2020, an indication of decreased efficiencies of those wells to supply water or of changes in groundwater yield around the well. Increased slopes of best-fit lines relating water-level difference to water-withdrawal rate from 2019 to 2020, when groundwater levels were otherwise relatively stable, also indicated that the specific capacity of one or more TCWS collector wells may have declined during that period. Riverbed surface elevation, roughness, and hardness changes between hydroacoustic surveys indicated active accretion and erosion of bed sediment offshore from the riverbank filtration system in 2002–4. Riverbed hydraulic conductivity computed from 2013 to 2015 in tests at 11 sites ranged from 0.2 to 230 feet per day. Ten temporary piezometer sites with adjacent operating collector wells had drawdowns that ranged from about 4 to 11.7 feet and hammer blows greater than or equal to 49 counts, which indicated a relatively great density of alluvial sediments and the potential for diminished riverbed infiltration rates. Contributions of Ohio River infiltration to collector wells estimated from deuterium and oxygen-18 isotopic ratios in representative Ohio River and groundwater samples were considerably larger in warmer months than cooler months and ranged from 80–90 percent in CW1–CW5 in May 2020 to 100 percent in May 2019 at CW2, CW4, and CW5. Groundwater contributions to water withdrawals by collector wells increased in cooler autumn months and ranged from 40 percent at CW3 and CW5 in November 2020 to 70 percent at CW2, CW4, and CW5 in November–December 2019. November 2019 and November 2020 isotopic compositions of water from the collector wells had enriched deuterium and oxygen-18 isotopic ratios consistent with mixtures of evaporation-modified water from the Ohio River and meteoric-origin groundwater.

Indiana, Kentucky↗
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↗