Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Complexity”

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

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

At least 1,477 records · Page 82Linked to original sources
USGS2017

Field-trip guide to Columbia River flood basalts, associated rhyolites, and diverse post-plume volcanism in eastern Oregon

The Miocene Columbia River Basalt Group (CRBG) is the youngest and best preserved continental flood basalt province on Earth, linked in space and time with a compositionally diverse succession of volcanic rocks that partially record the apparent emergence and passage of the Yellowstone plume head through eastern Oregon during the late Cenozoic. This compositionally diverse suite of volcanic rocks are considered part of the La Grande - Owyhee eruptive axis (LOEA), an approximately 300-kilometer-long (185 mile), north - northwest - trending, middle Miocene to Pliocene volcanic belt located along the eastern margin of the Columbia River flood basalt province. Volcanic rocks erupted from and preserved within the LOEA form an important regional stratigraphic link between the (1) flood basalt - dominated Columbia Plateau on the north, (2) bimodal basalt - rhyolite vent complexes of the Owyhee Plateau on the south, (3) bimodal basalt-rhyolite and time-transgressive rhyolitic volcanic fields of the Snake River Plain - Yellowstone Plateau, and (4) the High Lava Plains of central Oregon. This field-trip guide describes a 4 - day geologic excursion that will explore the stratigraphic and geochemical relationships among mafic rocks of the Columbia River Basalt Group and coeval and compositionally diverse volcanic rocks associated with the early “Yellowstone track” and High Lava Plains in eastern Oregon. Beginning in Portland, the Day 1 log traverses the Columbia River gorge eastward to Baker City, focusing on prominent outcrops that reveal a distal succession of laterally extensive, large-volume tholeiitic flood lavas of the Grande Ronde, Wanapum, and Saddle Mountains Basalt formations of the CRBG. These “great flows” are typical of the well-studied flood basalt-dominated Columbia Plateau, where interbedded silicic and calc-alkaline lavas are conspicuously absent. The latter part of Day 1 will highlight exposures of middle to late Miocene silicic ash-flow tuffs, rhyolite domes, and calc-alkaline lava flows overlying the CRBG across the northern and central parts of the LOEA. The Day 2 field route migrates to southern parts of the LOEA, where rocks of the CRBG are associated in space and time with lesser known and more complex silicic volcanic stratigraphy associated with middle Miocene, large-volume, bimodal basalt-rhyolite vent complexes. Key stops will provide a broad overview of the structure and stratigraphy of the middle Miocene Mahogany Mountain caldera and middle to late Miocene calc-alkaline lavas of the Owyhee basalt. Stops on Day 3 will progress westward from the eastern margin of the LOEA, examining a transition linking the Columbia River Basalt-Yellowstone province with a northwestward-younging magmatic trend of silicic volcanism that underlies the High Lava Plains of eastern Oregon. Initial field stops on Day 3 will examine key outcrops demonstrating the intercalated nature of middle Miocene tholeiitic CRBG flood basalts, prominent ash-flow tuffs, and “Snake River-type” large-volume rhyolite lava flows exposed along the Malheur River. Subsequent stops on Day 3 will focus upon the volcanic stratigraphy northeast of the town of Burns, which includes regional middle to late Miocene ash-flow tuffs, and lava flows assigned to the Strawberry Volcanics. The return route to Portland on Day 4 traverses across the western axis of the Blue Mountains, highlighting exposures of the widespread, middle Miocene Dinner Creek Tuff and aspects of Picture Gorge Basalt flows and northwest-trending feeder dikes situated in the central part of the CRBG province.

Oregon↗
USGS2017

U.S. Geological Survey Karst Interest Group Proceedings, San Antonio, Texas, May 16–18, 2017

Introduction and Acknowledgments Karst aquifer systems are present throughout parts of the United States and some of its territories, and have developed in carbonate rocks (primarily limestone and dolomite) and evaporites (gypsum, anhydrite, and halite) that span an interval of time encompassing more than 550 million years. The depositional environments, diagenetic processes, post-depositional tectonic events, and geochemical weathering processes that form karst aquifers are varied and complex. These factors involve biological, chemical, and physical changes that when combined with the diverse climatic regimes in which karst development has taken place, result in the unique dual- or triple-porosity nature of karst aquifers. These complex hydrogeologic systems typically represent challenging and unique conditions to scientists attempting to study groundwater flow and contaminant transport in these terrains. The dissolution of carbonate rocks and the subsequent development of distinct and beautiful landscapes, caverns, and springs have resulted in the most exceptional karst areas being designated as national or state parks. Tens of thousands of similar areas in the United States have been developed into commercial caverns and known privately owned caves. Both public and private properties provide access for scientists to study the flow of groundwater in situ . Likewise, the range and complexity of landforms and groundwater flow systems associated with karst terrains are enormous, perhaps more than for any other aquifer type. Karst aquifers and landscapes that form in tropical areas, such as the cockpit karst along the north coast of Puerto Rico, differ greatly from karst landforms in more arid climates, such as the Edwards Plateau in west-central Texas or the Guadalupe Mountains near Carlsbad, New Mexico, where hypogenic processes have played a major role in speleogenesis. Many of these public and private lands also contain unique flora and fauna associated with these karst hydrogeologic systems. As a result, numerous federal, state, and local agencies have a strong interest in the study of karst terrains. Many of the major springs and aquifers in the United States have developed in carbonate rocks, such as the Floridan aquifer system in Florida and parts of Alabama, Georgia, and South Carolina; the Ozark Plateaus aquifer system in parts of Arkansas, Kansas, Missouri, and Oklahoma; and the Edwards-Trinity aquifer system in west-central Texas. These aquifers, and the springs that discharge from them, serve as major water-supply sources and form unique ecological habitats. Competition for the water resources of karst aquifers is common, and urban development and the lack of attenuation of contaminants in karst areas due to dissolution features that form direct pathways into karst aquifers can impact the ecosystem and water quality associated with these aquifers. The concept for developing a platform for interaction among scientists within the U.S. Geological Survey (USGS) working on karst-related studies evolved from the November 1999 National Groundwater Meeting of the USGS. As a result, the Karst Interest Group (KIG) was formed in 2000. The KIG is a loose-knit, grass-roots organization of USGS and non-USGS scientists and researchers devoted to fostering better communication among scientists working on, or interested in, karst science. The primary mission of the KIG is to encourage and support interdisciplinary collaboration and technology transfer among scientists working in karst areas. Additionally, the KIG encourages collaborative studies between the different mission areas of the USGS as well as with other federal and state agencies, and with researchers from academia and institutes. To accomplish its mission, the KIG has organized a series of workshops that have been held near nationally important karst areas. To date (2017) seven KIG workshops, including the workshop documented in this report, have been held. The workshops typically include oral and poster sessions on selected karst-related topics and research, as well as field trips to local karst areas. To increase non-USGS participation an effort was made for the workshops to be held at a university or institute beginning with the fourth workshop. Proceedings of the workshops are published by the USGS and are available online at the USGS publications warehouse https://pubs.er.usgs.gov/ by using the search term “karst interest group.” The first KIG workshop was held in St. Petersburg, Florida, in 2001, in the vicinity of the large springs and other karst features of the Floridan aquifer system. The second KIG workshop was held in 2002, in Shepherdstown, West Virginia, in proximity to the carbonate aquifers of the northern Shenandoah Valley, and highlighted an invited presentation on karst literature by the late Barry F. Beck of P.E. LaMoreaux and Associates. The third KIG workshop was held in 2005, in Rapid City, South Dakota, near evaporite karst features in limestones of the Madison Group in the Black Hills of South Dakota. The Rapid City KIG workshop included field trips to Wind Cave National Park and Jewel Cave National Monument, and featured a presentation by Thomas Casadevall, then USGS Central Region Director, on the status of Earth science at the USGS. The fourth KIG workshop in 2008 was hosted by the Hoffman Environmental Research Institute and Center for Cave and Karst Studies at Western Kentucky University in Bowling Green, Kentucky, near Mammoth Cave National Park and karst features of the Chester Upland and Pennyroyal Plateau. The workshop featured a late-night field trip into Mammoth Cave led by Rickard Toomey and Rick Olsen, National Park Service. The fifth KIG workshop in 2011 was a joint meeting of the USGS KIG and University of Arkansas HydroDays, hosted by the Department of Geosciences at the University of Arkansas in Fayetteville. The workshop featured an outstanding field trip to the unique karst terrain along the Buffalo National River in the southern Ozarks, and a keynote presentation on paleokarst in the United States was delivered by Art and Peggy Palmer. The sixth KIG workshop was hosted by the National Cave and Karst Research Institute (NCKRI) in 2014, in Carlsbad, New Mexico. George Veni, Director of the NCKRI, served as a co-chair of the workshop with Eve Kuniansky of the USGS. The workshop featured speaker Dr. Penelope Boston, Director of Cave and Karst Studies at New Mexico Tech, Socorro, and Academic Director at the NCKRI, who addressed the future of karst research. The field trip on evaporite karst of the lower Pecos Valley was led by Lewis Land (NCKRI karst hydrologist), and the field trip on the geology of Carlsbad Caverns National Park was led by George Veni. This current seventh KIG workshop is being held in San Antonio at the University of Texas at San Antonio (UTSA). This 2017 workshop is being hosted by the Department of Geological Sciences’ Student Geological Society (SGS), and student chapters of the American Association of Petroleum Geologists (AAPG) and Association of Engineering Geologists (AEG), with support by the UTSA Department of Geological Sciences and Center for Water Research. The UTSA student chapter presidents, Jose Silvestre (SGS), John Cooper (AAPG), and Tyler Mead (AEG) serve as co-chairs of the 2017 workshop with Eve Kuniansky of the USGS. The technical session committee is chaired by Eve Kuniansky, USGS, and includes Michael Bradley, Tom Byl, Rebecca Lambert, John Lane, and James Kaufmann, all USGS, and Patrick Tucci, retired USGS. The logistics committee includes Amy Clark, Yongli Gao, and Lance Lambert (Department Chair), UTSA Department of Geological Sciences; and Ryan Banta and Allan Clark, USGS, San Antonio, Texas. The field trip committee is chaired by Allan Clark and includes Amy Clark, Yongli Gao, and Keith Muehlestein, UTSA; Marcus Gary, Edwards Aquifer Authority and University of Texas at Austin; Ron Green, Southwest Research Institute; Geary Schindel, Edwards Aquifer Authority; and George Veni, NCKRI. Additionally, two organizations have assisted the UTSA student chapters in hosting the meeting by donating funds to the chapters: the Edwards Aquifer Authority, San Antonio, Texas, and the Barton Springs Edwards Aquifer Authority, Austin, Texas. Additionally, Yongli Gao, Center for Water Research and Department of Geological Sciences, UTSA, helped develop sessions on cave and karst research in China for this workshop. These proceedings could not have been accomplished without the assistance of Lawrence E. Spangler as co-editor who not only has subject matter expertise, but also serves as an editor with the USGS Science Publishing Network. We sincerely hope that this workshop continues to promote future collaboration among scientists of varied and diverse backgrounds, and improves our understanding of karst aquifer systems in the United States and its territories. The extended abstracts of USGS authors were peer reviewed and approved for publication by the USGS. Articles submitted by university researchers and other federal and state agencies did not go through the formal USGS peer review and approval process, and therefore may not adhere to USGS editorial standards or stratigraphic nomenclature. However, all articles had a minimum of two peer reviews and were edited for consistency of appearance in the proceedings. The use of trade, firm or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government. The USGS Water Availability and Use Science Program funded the publication costs of the proceedings.

Scientific Investigations Report↗
USGS2017

A review of surface energy balance models for estimating actual evapotranspiration with remote sensing at high spatiotemporal resolution over large extents

Many approaches have been developed for measuring or estimating actual evapotranspiration ( ET a ), and research over many years has led to the development of remote sensing methods that are reliably reproducible and effective in estimating ET a . Several remote sensing methods can be used to estimate ET a at the high spatial resolution of agricultural fields and the large extent of river basins. More complex remote sensing methods apply an analytical approach to ET a estimation using physically based models of varied complexity that require a combination of ground-based and remote sensing data, and are grounded in the theory behind the surface energy balance model. This report, funded through cooperation with the International Joint Commission, provides an overview of selected remote sensing methods used for estimating water consumed through ET a and focuses on Mapping Evapotranspiration at High Resolution with Internalized Calibration (METRIC) and Operational Simplified Surface Energy Balance (SSEBop), two energy balance models for estimating ET a that are currently applied successfully in the United States. The METRIC model can produce maps of ET a at high spatial resolution (30 meters using Landsat data) for specific areas smaller than several hundred square kilometers in extent, an improvement in practice over methods used more generally at larger scales. Many studies validating METRIC estimates of ET a against measurements from lysimeters have shown model accuracies on daily to seasonal time scales ranging from 85 to 95 percent. The METRIC model is accurate, but the greater complexity of METRIC results in greater data requirements, and the internalized calibration of METRIC leads to greater skill required for implementation. In contrast, SSEBop is a simpler model, having reduced data requirements and greater ease of implementation without a substantial loss of accuracy in estimating ET a . The SSEBop model has been used to produce maps of ET a over very large extents (the conterminous United States) using lower spatial resolution (1 kilometer) Moderate Resolution Imaging Spectroradiometer (MODIS) data. Model accuracies ranging from 80 to 95 percent on daily to annual time scales have been shown in numerous studies that validated ET a estimates from SSEBop against eddy covariance measurements. The METRIC and SSEBop models can incorporate low and high spatial resolution data from MODIS and Landsat, but the high spatiotemporal resolution of ET a estimates using Landsat data over large extents takes immense computing power. Cloud computing is providing an opportunity for processing an increasing amount of geospatial “big data” in a decreasing period of time. For example, Google Earth Engine TM has been used to implement METRIC with automated calibration for regional-scale estimates of ET a using Landsat data. The U.S. Geological Survey also is using Google Earth Engine TM to implement SSEBop for estimating ET a in the United States at a continental scale using Landsat data.

Scientific Investigations Report↗
USGS2017

Correlation between basalt flows and radiochemical and chemical constituents in selected wells in the southwestern part of the Idaho National Laboratory, Idaho

Wastewater discharged to wells and ponds and wastes buried in shallow pits and trenches at facilities at the Idaho National Laboratory (INL) have contributed contaminants to the eastern Snake River Plain (ESRP) aquifer in the southwestern part of the INL. This report describes the correlation between subsurface stratigraphy in the southwestern part of the INL with information on the presence or absence of wastewater constituents to better understand how flow pathways in the aquifer control the movement of wastewater discharged at INL facilities. Paleomagnetic inclination was used to identify subsurface basalt flows based on similar inclination measurements, polarity, and stratigraphic position. Tritium concentrations, along with other chemical information for wells where tritium concentrations were lacking, were used as an indicator of which wells were influenced by wastewater disposal. The basalt lava flows in the upper 150 feet of the ESRP aquifer where wastewater was discharged at the Idaho Nuclear Technology and Engineering Center (INTEC) consisted of the Central Facilities Area (CFA) Buried Vent flow and the AEC Butte flow. At the Advanced Test Reactor (ATR) Complex, where wastewater would presumably pond on the surface of the water table, the CFA Buried Vent flow probably occurs as the primary stratigraphic unit present; however, AEC Butte flow also could be present at some of the locations. At the Radioactive Waste Management Complex (RWMC), where contamination from buried wastes would presumably move down through the unsaturated zone and pond on the surface of the water table, the CFA Buried Vent; Late Basal Brunhes; or Early Basal Brunhes basalt flows are the flow unit at or near the water table in different cores. In the wells closer to where wastewater disposal occurred at INTEC and the ATR-Complex, almost all the wells show wastewater influence in the upper part of the ESRP aquifer and wastewater is present in both the CFA Buried Vent flow and AEC Butte flow. The CFA Buried Vent flow and AEC Butte flow are also present in wells at and north of CFA and are all influenced by wastewater contamination. All wells with the AEC Butte flow present have wastewater influence and 83 percent of the wells with the more prevalent CFA Buried Vent flow have wastewater influence. South and southeast of CFA, most wells are not influenced by wastewater disposal and are completed in the Big Lost Flow and the CFA Buried Vent flow. Wells southwest of CFA are influenced by wastewater disposal and are completed in the Big Lost flow and CFA Buried Vent flow at the top of the aquifer. Basalt stratigraphy indicates that the CFA Buried Vent flow is the predominant flow in the upper part of the ESRP aquifer at and near the RWMC as it is present in all the wells in this area. The Late Basal Brunhes flow, Middle Basal Brunhes flow, Early Basal Brunhes flow, South Late Matuyama flow, and Matuyama flow are also present in various wells influenced by waste disposal. Some wells south of RWMC do not show wastewater influence, and the lack of wastewater influence could be due to low hydraulic conductivities. Several wells south and southeast of CFA also do not show wastewater influence. Low hydraulic conductivities or ESRP subsidence are possible causes for lack of wastewater south of CFA. Multilevel monitoring wells completed much deeper in the aquifer show influence of wastewater in numerous basalt flows. Well Middle 2051 (northwest of RWMC) does not show wastewater influence in its upper three basalt flows (CFA Buried Vent, Late Basal Brunhes, and Middle Basal Brunhes); however, wastewater is present in two deeper flows (the Matuyama and Jaramillo flows). Well USGS 131A (southwest of CFA) and USGS132 (south of RWMC) both show wastewater influence in all the basalt flows sampled in the upper 600 feet of the aquifer. Wells USGS 137A, 105, 108, and 103 completed along the southern boundary of the INL all show wastewater influence in several basalt flows including the G flow, Middle and Early Basal Brunhes flows, the South Late Matuyama flow and the Matuyama flow; however, the strongest wastewater influence appears to be in the South Late Matuyama flow. The concentrations of wastewater constituents in deeper parts of these wells support the concept of groundwater flow deepening in the southwestern part of the INL.

Idaho↗
USGS2018

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

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

California↗
USGS2018

Mercury methylation and bioaccumulation in Sinclair Inlet, Kitsap County, Washington

The U.S. Geological Survey evaluated the transformation of mercury to bioavailable methylmercury in Sinclair Inlet, Kitsap County, Washington, and assessed the effect of the transformation processes on the mercury burden in marine organisms and sediment. In August 2008, samples of sediment, water, and biota from six sites in Sinclair Inlet and three bays representative of Puget Sound embayments were collected. The extensive sediment sampling included analysis of methylmercury in sediment and porewater, estimates of methylation production potential, and analyses of ancillary constituents associated with organic carbon and reduction-oxidation (redox) conditions to assist in interpreting the mercury results. Analyses of methylmercury in water overlying incubated cores provided an estimate of the release of methylmercury to the water column. Collection of samples for mercury species in the aqueous, particulate (suspended solids), and biological phases, and for ancillary carbon and nitrogen constituents in surface water, continued, on about a monthly schedule, at four stations through August 2009. In February, June, and August 2009, seasonal sediment samples were collected at 20 stations distributed between greater Sinclair Inlet and Operable Unit B Marine of the Bremerton naval complex, Bremerton, Washington, to examine geographical and seasonal patterns of mercury biogeochemistry of sediment in Sinclair Inlet. At six of these seasonal sediment stations, porewater was collected and triplicate core incubation experiments were done. Median sediment-methylmercury concentrations were not statistically different between the representative bays and Sinclair Inlet. The percentage of sediment methylmercury (relative to total mercury) was actually lower in the Sinclair Inlet sites compared with the representative bays, reflecting the higher sediment total mercury concentration for the Sinclair Inlet stations compared with the representative bays. Likewise, median sediment methylmercury concentrations were not statistically different between the greater Sinclair Inlet stations and the Bremerton naval complex stations; whereas the percentage of sediment methylmercury to total mercury was lower in the Bremerton naval complex due to higher sediment total mercury concentrations than the greater Sinclair Inlet stations. The biogeochemical characteristics of each station, measured by redox, organic carbon, and the seasonal availability of nutrients controlled methylmercury biogeochemistry. Mercury methylation production potential was a function of temperature, concentration of total mercury in sediment, and the percentage of ferrous iron (relative to total measured iron) across all sites. Methylmercury porewater concentrations were best described by using concentrations of dissolved organic carbon and reduction-oxidation conditions. Likewise, the variable fluxes of methylmercury from incubated cores were best described using dissolved organic carbon and reduction-oxidation conditions. Sinclair Inlet exhibited the classic Puget Sound biological cycle, with spring and autumn phytoplankton blooms resulting in depletion of nitrate, orthophosphate, and silicate in the surface water. Although variable in timing between 2008 and 2009, a strong corresponding seasonal trend of increased availability, incorporation, and bioaccumulation of methylmercury into the food web of Sinclair Inlet occurred during the early spring and summer growing season.

Washington↗
USGS2019

Tritium as an indicator of modern, mixed, and premodern groundwater age

Categorical classification of groundwater age is often used for the assessment and understanding of groundwater resources. This report presents a tritium-based age classification system for the conterminous United States based on tritium ( 3 H) thresholds that vary in space and time: modern (recharged in 1953 or later), if the measured value is larger than an upper threshold; premodern (recharged prior to 1953) if the measured value is smaller than a lower threshold; or mixed if the measured value is between the two thresholds. Inclusion of spatially varying thresholds, rather than a single threshold, accounts for the observed systematic variation in 3 H deposition across the United States. Inclusion of time-varying thresholds, rather than a single threshold, accounts for the date of sampling given the radioactive decay of 3 H. The efficacy of the tritium-based age classification system was evaluated at national and regional scales. The system was evaluated at a national scale by classifying samples from 1,788 public-supply wells distributed across 19 principal aquifers and comparing those results with expectations based on hydrogeologic principles. The regional-scale data are from five paired networks of shallow and deep wells (287 wells). As expected, modern groundwater is more prevalent in shallow wells than in deeper wells, in fractured-rock and carbonate aquifers as compared to clastic aquifers, in unconfined areas as compared to confined areas, and in humid climates as compared to arid climates. The results from a tritium-based age classification system compared favorably with the results of 14 previous studies of groundwater ages that used different age tracers and analytical methods. The wells and samples from the Cambrian-Ordovician aquifer that had been analyzed using a more complex multi-tracer analysis were also analyzed using the tritium-based age classification system, and there was a close match between the two methods. The results from these various studies suggest that the tritium-based age classification system may be informative as a screening tool prior to selecting more expensive and complex age-dating tracers and methods, or to provide an explanatory variable for other water-quality data where more complex methods or tracers are not available. This work improves on previous groundwater age classification using 3 H by developing methods that (1) determine 3 H thresholds for groundwater recharged in 1953 or later that minimize the misclassification of modern samples as mixed; (2) determine a pre-1953 threshold to estimate premodern background concentrations; and (3) add a mixed category to classify samples that are clearly neither entirely modern nor entirely premodern. As with any tritium-based approach, it can fail when the 3 H record in precipitation does not accurately reflect the record of 3 H in recharge

Conterminous United States↗
USGS2022

Using microbial source tracking to identify contamination sources in Port Jefferson Harbor, Setauket Harbor, and Conscience Bay on Long Island, New York

The U.S. Geological Survey worked in cooperation with the New York State Department of Environmental Conservation to assess the potential sources of fecal contamination entering Port Jefferson Harbor, Setauket Harbor, and Conscience Bay, an embayment complex on the northern shore of Suffolk County, Long Island, New York. Water samples are routinely collected by the New York State Department of Environmental Conservation in the harbor and analyzed for fecal coliform bacteria, an indicator of fecal contamination, to determine the need for closure of shellfish beds for harvest and consumption. Fecal coliform and other bacteria are an indicator of the potential presence of pathogenic (disease-causing) bacteria. However, indicator bacteria alone cannot determine the biological or geographical sources of contamination; therefore, microbial source tracking was implemented to determine various biological sources of contamination. In addition, information such as the location, weather and season, and surrounding land use where a sample was collected help determine the geographical source and conveyance of land-based water to the embayment. Our analysis revealed that the most substantial source of fecal contamination to the Port Jefferson Harbor complex was discharge from sites draining ponds and wetlands, particularly during the summer months. Fecal coliform bacteria at sites where ponds and wetlands drain are increased by stormwater runoff, which is another substantial source of fecal contamination. Human markers, likely originating from the outfall of the Port Jefferson sewage treatment plant, were found at least once in every sample collected within Port Jefferson Harbor; however, fecal coliform concentrations were low, indicating that the sewage treatment plant is not a likely source of fecal contamination to the embayment. Canine markers detected in Conscience Bay were associated with high fecal coliform, particularly in wet summer samples. Waterfowl markers were most often detected in source water for Port Jefferson Harbor, but in Conscience Bay, waterfowl was frequently detected within the bay itself. Resuspension of bed sediment may contribute to fecal contamination in the harbor, but more targeted analyses are needed to support this finding. There was little evidence of groundwater-contributing fecal bacteria by direct discharge from the subsurface. A classification scheme was developed to convey the degree of fecal contamination to stakeholders and resource managers. Based on this classification scheme, the Culvert North of State Route 25A, Culvert North of Shore Road, Old Mill Creek Culvert, and Mill Pond Culvert to Conscience Bay sites were identified as locations that contribute substantial fecal contamination to the Port Jefferson Harbor complex.

New York↗
USGS2024

Evaluation of sensors for continuous monitoring of harmful algal blooms in the Finger Lakes region, New York, 2019 and 2020

In response to the increasing frequency of cyanobacterial harmful algal blooms (CyanoHABs) in the Finger Lakes region of New York State, a pilot study by the U.S. Geological Survey, in collaboration with the New York State Department of Environmental Conservation, was conducted to enhance CyanoHAB monitoring and understanding. High-frequency sensors were deployed on open water monitoring-station platforms at Seneca Lake in 2019–20, at Owasco Lake in 2019–20, and at Skaneateles Lake in 2019. One of the goals of this study was to evaluate the ability of in-place sensors to make representative measurements of dissolved organic matter, nutrients, and algal pigments (as indicators of phytoplankton biomass) while collecting routine field parameters (water temperature, specific conductance, pH, dissolved oxygen, turbidity, weather, and light) to provide additional information about environmental conditions. Despite challenges like power issues and sensor fouling, the sensors performed well overall. However, correlation analyses between sensor readings and laboratory measurements revealed variable performance. Results indicate the relation between the fluorescent dissolved organic matter sensor and laboratory-measured dissolved organic carbon was weak at all study lakes. The nitrate sensors can be sensitive to ambient temperature and have a substantial power requirement, and the relation between sensor- and laboratory-measured nitrate values differed among lakes. The orthophosphate sensors, which were complex and prone to data loss, yielded results that were difficult to interpret because orthophosphate detections are rare in the study lakes. The multichannel fluorometer was also complex to use and required several unique procedures for its operation. Chlorophyll measurements from the fluorometers correlated moderately well with laboratory-measured chlorophyll- a , although relations with total phytoplankton biovolume were weaker. Relations between phycocyanin concentration measurements from the dual-channel fluorometers and cyanobacterial biovolume were not significant; however, the cyanobacterial biovolume correlation was moderately strong with chlorophyll contribution from cyanobacteria measurements from the multichannel fluorometer. Of all collected parameters, water temperature was among the strongest correlated with chlorophyll- a , total phytoplankton biovolume, and cyanobacterial biovolume. Stepwise regression analysis was used to identify the best parameters for modeling variance in laboratory measures of phytoplankton biomass. This analysis included factors such as chlorophyll fluorescence, pH, water temperature, and others, which varied by lake. Overall, the models had limited explanatory power for chlorophyll- a and other biovolumes, possibly due to the absence of CyanoHABs at the open-water monitoring locations. Multivariate models did not outperform simple fluorescence-based models. Notably, turbidity was a more significant indicator of cyanobacterial biovolume variability than phycocyanin from dual-channel fluorometers. The study concludes that while single and multivariate models based on sensor data are useful, they did not explain any more variance than fluorescence-based models. Broader data collection, including more CyanoHAB events, is necessary to refine these models. Integrating machine learning could leverage large, complex datasets to improve CyanoHAB predictions, thereby enhancing the management and understanding of these blooms.

New York↗
USGS1956

A fluorimetric study of the thorium-morin system

Thorium reacts with morin to yield a yellow complex that fluoresces when irradiated with ultraviolet light. The effect on the fluorescence of such variable as concentration of acid, alcohol, thorium, morin, and complex; time, temperature, and wavelength of exciting light are studied to determine experimental conditions yielding maximum fluorescence. The effects of Zr 4+ , Al 3+ , Fe 3+ , Ca 2+ , and La 3+ are discussed. The fundamental relationships between light absorption and fluorescence are expressed in a general equation which applied to a three-component system when the fluorescence is measured in a transmission-type fluorimeter. This general equation is used to obtain an expression for the fluorescence of the thorium-morin system. Equations, derived from experimental data, related both the fraction of thorium reacted to form complex and the fraction of unquenched fluorescence to the concentration of uncombined morin. These functions, when combined with the general equation, give an expression which relates the total net fluorescence to the amount of uncombined morin in the solution. This last equation can be used to determine the one region for the concentration of uncombined morin that gives maximum sensitivity for the system. Calculated standard curves are in excellent agreement with experimental curves.

Trace Elements Investigations↗
USGS2003

Geochemistry of the Birch Creek Drainage Basin, Idaho

The U.S. Survey and Idaho State University, in cooperation with the U.S. Department of Energy, are conducting studies to describe the chemical character of ground water that moves as underflow from drainage basins into the eastern Snake River Plain aquifer (ESRPA) system at and near the Idaho National Engineering and Environmental Laboratory (INEEL) and the effects of these recharge waters on the geochemistry of the ESRPA system. Each of these recharge waters has a hydrochemical character related to geochemical processes, especially water-rock interactions, that occur during migration to the ESRPA. Results of these studies will benefit ongoing and planned geochemical modeling of the ESRPA at the INEEL by providing model input on the hydrochemical character of water from each drainage basin. During 2000, water samples were collected from five wells and one surface-water site in the Birch Creek drainage basin and analyzed for selected inorganic constituents, nutrients, dissolved organic carbon, tritium, measurements of gross alpha and beta radioactivity, and stable isotopes. Four duplicate samples also were collected for quality assurance. Results, which include analyses of samples previously collected from four other sites, in the basin, show that most water from the Birch Creek drainage basin has a calcium-magnesium bicarbonate character. The Birch Creek Valley can be divided roughly into three hydrologic areas. In the northern part, ground water is forced to the surface by a basalt barrier and the sampling sites were either surface water or shallow wells. Water chemistry in this area was characterized by simple evaporation models, simple calcite-carbon dioxide models, or complex models involving carbonate and silicate minerals. The central part of the valley is filled by sedimentary material and the sampling sites were wells that are deeper than those in the northern part. Water chemistry in this area was characterized by simple calcite-dolomite-carbon dioxide models. In the southern part, ground water enters the ESRPA. In this area, the sampling sites were wells with depths and water levels much deeper than those in the northern and central parts of the valley. The calcium and carbon water chemistry in this area was characterized by a simple calcite-carbon dioxide model, but complex calcite-silicate models more accurately accounted for mass transfer in these areas. Throughout the geochemical system, calcite precipitated if it was an active phase in the models. Carbon dioxide either precipitated (outgassed) or dissolved depending on the partial pressure of carbon dioxide in water from the modeled sites. Dolomite was an active phase only in models from the central part of the system. Generally the entire geochemical system could be modeled with either evaporative models, carbonate models, or carbonate-silicate models. In both of the latter types of models, a significant amount of calcite precipitated relative to the mass transfer to and from the other active phases. The amount of calcite precipitated in the more complex models was consistent with the amount of calcite precipitated in the simpler models. This consistency suggests that, although the simpler models can predict calcium and carbon concentrations in Birch Creek Valley ground and surface water, silicate-mineral-based models are required to account for the other constituents. The amount of mass transfer to and from the silicate mineral phases was generally small compared with that in the carbonate phases. It appears that the water chemistry of well USGS 126B represents the chemistry of water recharging the ESRPA by means of underflow from the Birch Creek Valley.

Idaho↗
USGS1983

Determination of dissolved aluminum in water samples

A technique has been modified for determination of a wide range of concentrations of dissolved aluminum (Al) in water and has been tested. In this technique, aluminum is complexed with 8-hydroxyquinoline at pH 8.3 to minimize interferences, then extracted with methyl isobutyl ketone (MIBK). The extract is analyzed colorimetrically at 395 nm. This technique is used to analyze two forms of monomeric Al, nonlabile (organic complexes) and labile (free, Al, Al sulfate, fluoride and hydroxide complexes). A detection limit 2 ug/L is possible with 25-ml samples and 10-ml extracts. The detection limit can be decreased by increasing the volume of the sample and (or) decreasing the volume of the methyl isobutyl ketone extract. The analytical uncertainty of this method is approximately + or - 5 percent. The standard addition technique provides a recovery test for this technique and ensures precision in samples of low Al concentrations. The average percentage recovery of the added Al plus the amount originally present was 99 percent. Data obtained from analyses of filtered standard solutions indicated that Al is adsorbed on various types of filters. However, the relationship between Al concentrations and adsorption remains linear. A test on standard solutions also indicated that Al is not adsorbed on nitric acid-washed polyethylene and polypropylene bottle wells. (USGS)

Water-Resources Investigations Report↗
USGS1987

Cenozoic stratigraphy and geologic history of the Tucson Basin, Pima County, Arizona

This report was prepared as part of a geohydrologic study of the Tucson basin conducted by the U.S. Geological Survey in cooperation with the city of Tucson. Geologic data from more than 500 water supply and test wells were analyzed to define characteristics of the basin sediments that may affect the potential for land subsidence induced by groundwater withdrawal. The Tucson basin is a structural depression within the Basin and Range physiographic province. The basin is 1,000 sq mi in units area and trends north to northwest. Three Cenozoic stratigraphic unit--the Pantano Formation of Oligocene age, the Tinaja beds (informal usage) of Miocene and Pliocene age, and the Fort Lowell Formation of Pleistocene age--fill the basin. The Tinaja beds include lower, middle, and upper unconformable units. A thin veneer of stream alluvium of late Quaternary age overlies the Fort Lowell Formation. The Pantano Formation and the lower Tinaja beds accumulated during a time of widespread continental sedimentation, volcanism, plutonism, uplift, and complex faulting and tilting of rock units that began during the Oligocene and continued until the middle Miocene. Overlying sediments of the middle and upper Tinaja beds were deposited in response to two subsequent episodes of post-12-million-year block faulting, the latter of which was accompanied by renewed uplift. The Fort Lowell Formation accumulated during the Quaternary development of modern through-flowing the maturation of the drainage. The composite Cenozoic stratigraphic section of the Tucson basin is at least 20,000 ft thick. The steeply tilted to flat-lying section is composed of indurated to unconsolidated clastic sediments, evaporites, and volcanic rocks that are lithologically and structurally complex. The lithology and structures of the section was greatly affected by the uplift and exhumation of adjacent metamorphic core-complex rocks. Similar Cenozoic geologic relations have been identified in other parts of southern Arizona. (Author 's abstract)

Water-Resources Investigations Report↗
USGS1995

Geohydrology and simulation of ground-water flow in the aquifer system near Calvert City, Kentucky

The U.S. Geological Survey, in cooperation with the Kentucky Natural Resources and Environmental Protection Cabinet, constructed a two-dimensional, steady-state ground-water-flow model to estimate hydraulic properties, contributing areas to discharge boundaries, and the average linear velocity at selected locations in an aquifer system near Calvert City, Ky. Nonlinear regression was used to estimate values of model parameters and the reliability of the parameter estimates. The regression minimizes the weighted difference between observed and calculated hydraulic heads and rates of flow. The calibrated model generally was better than alternative models considered, and although adding transmissive faults in the bedrock produced a slightly better model, fault transmissivity was not estimated reliably. The average transmissivity of the aquifer was 20,000 feet squared per day. Recharge to two outcrop areas, the McNairy Formation of Cretaceous age and the alluvium of Quaternary age, were 0.00269 feet per day (11.8 inches per year) and 0.000484 feet per day (2.1 inches per year), respectively. Contributing areas to wells at the Calvert City Water Company in 1992 did not include the Calvert City Industrial Complex. Since completing the fieldwork for this study in 1992, the Calvert City Water Company discontinued use of their wells and began withdrawing water from new wells that were located 4.5 miles east-southeast of the previous location; the contributing area moved farther from the industrial complex. The extent of the alluvium contributing water to wells was limited by the overlying lacustrine deposits. The average linear ground-water velocity at the industrial complex ranged from 0.90 feet per day to 4.47 feet per day with a mean of 1.98 feet per day.

Water-Resources Investigations Report↗
USGS1999

Geohydrology and numerical simulation of the ground-water flow system of Kona, Island of Hawaii

Prior to the early 1990's, ground-water in the Kona area, which is in the western part of the island of Hawaii, was withdrawn from wells located within about 3 mi from the coast where water levels were less than 10 feet above sea level. In 1990, exploratory drilling in the uplands east of the existing coastal wells first revealed the presence of high water levels (greater than 40 feet above sea level) in the Kona area. Measured water levels from 16 wells indicate that high water levels exist in a zone parallel to and inland of the Kona coast, between Kalaoa and Honaunau. Available hydrologic and geophysical evidence is generally consistent with the concept that the high ground-water levels are associated with a buried dike complex. A two-dimensional (areal), steady-state, freshwater-saltwater, sharp-interface ground-water flow model was developed for the Kona area of the island of Hawaii, to enhance the understanding of (1) the distribution of aquifer hydraulic properties, (2) the conceptual framework of the ground-water flow system, and (3) the regional effects of ground-water withdrawals on water levels and coastal discharge. The model uses the finite-difference code SHARP. To estimate the hydraulic characteristics, average recharge, withdrawals, and water-level conditions for the period 1991-93 were simulated. The following horizontal hydraulic-conductivity values were estimated: (1) 7,500 feet per day for the dike-free volcanic rocks of Hualalai and Mauna Loa, (2) 0.1 feet per day for the buried dike complex of Hualalai, (3) 10 feet per day for the northern marginal dike zone (north of Kalaoa), and (4) 0.5 feet per day for the southern marginal dike zone between Palani Junction and Holualoa. The coastal leakance was estimated to be 0.05 feet per day per foot. Measured water levels indicate that ground water generally flows from inland areas to the coast. Model results are in general agreement with the limited set of measured water levels in the Kona area. Model results indicate, however, that water levels do not strictly increase in an inland direction and that a ground-water divide exists within the buried dike complex. Data are not available, however, to verify model results in the area near and inland of the model-calculated ground-water divide. Three simulations to determine the effects of proposed withdrawals from the high water-level area on coastal discharge and water levels, relative to model-calculated, steady-state coastal discharge and water levels for 1997 withdrawal rates, show that the effects are widespread. During 1997, the total withdrawal of ground water from the high water-level area between Palani Junction and Holualoa was about 1 million gallons per day. Model results indicate that it may not be possible to withdraw 25.6 million gallons per day of freshwater from this area between Palani Junction and Holualoa, but that it may be possible to withdraw between 5 to 8 million gallons per day from the same area. For a proposed withdrawal rate of 5.0 million gallons per day uniformly distributed to 12 sites between Palani Junction and Holualoa, the model-calculated drawdown of 0.01 foot or more extends about 9 miles north-northwest and about 7 miles south of the proposed well sites. In all scenarios, freshwater coastal discharge is reduced by an amount equal to the additional freshwater withdrawal. Additional data needed to improve the understanding of the ground-water flow system in the Kona area include: (1) a wider spatial distribution and longer temporal distribution of water levels, (2) improved information about the subsurface geology, (3) independent estimates of hydraulic conductivity, (4) improved recharge estimates, and (5) information about the vertical distribution of salinity in ground water.

Hawaii↗
USGS2000

Site Selection for a Deep Monitor Well, Kualapuu, Molokai, Hawaii

Management of the ground-water resources near Kualapuu on the island of Molokai, Hawaii, is hindered by the uncertainty in the vertical salinity structure in the aquifer. In the State of Hawaii, vertical profiles of ground-water salinity are commonly obtained from deep monitor wells, and these profiles are used to estimate the thicknesses of the freshwater part of the ground-water flow system and the freshwater-saltwater transition zone. Information from a deep monitor well would improve the understanding of the ground-water flow system and the ability to effectively manage the ground-water resources near Kualapuu; however, as of mid-1999 no deep monitor wells had been drilled on the island of Molokai. Selection of an appropriate site for drilling a deep monitor well in the Kualapuu area depends partly on where future ground-water development may occur. Simulations using an areally two-dimensional, steady-state, sharp-interface ground-water flow model previously developed for the island of Molokai, Hawaii, indicate that the southeastern part of the Kualapuu area is a possible area of future ground-water development because (1) withdrawals from this area have a small effect on water levels at existing wells in the Kualapuu area (relative to effects from withdrawals in other parts of the Kualapuu area that are outside of the dike complex), and (2) model-calculated water levels in this part of the Kualapuu area are high relative to water levels in other parts of the Kualapuu area that are outside of the dike complex. Additional site-selection criteria include (1) ground-water level, (2) ground-surface altitude, (3) land classification, ownership, and accessibility, (4) geology, (5) locations of existing production wells, and (6) historical ground-water quality information. A deep monitor well in the Kualapuu area will likely be most useful for management purposes if it is located (1) in the vicinity of future ground-water development, (2) in an area where water levels are between about 8 and 12 feet above sea level, (3) at a ground-surface altitude that is between about 1,000 and 1,100 feet, (4) on government-owned land, (5) outside of the dike complex and as far from known volcanic vents as possible, (6) at least about 1,000 feet from, but within the same hydrogeologic setting as, existing or proposed production wells, and (7) east of well 0902-01. A viable area for drilling a deep monitor well is about a half mile southeast of existing wells 0801-01 to -03 and a half mile north of a known volcanic vent, Puu Luahine.

Water-Resources Investigations Report↗
USGS1969

Water resources of Windward Oahu, Hawaii

Windward Oahu lies in a large cavity--an erosional remnant of the Koolau volcanic dome at its greatest stage of growth. Outcrops include volcanic rocks associated with caldera collapse and the main fissure zone which is marked by a dike complex that extends along the main axis of the dome. The fissure zone intersects and underlies the Koolau Range north of Waiahole Valley. South of Waiahole Valley, the crest of the Koolau Range is in the marginal dike zone, an area of scattered dikes. The crest of the range forms the western boundary of windward Oahu. Dikes, mostly vertical and parallel or subparallel to the fissure zone, control movement and discharge of ground water because they are less permeable than the rocks they intrude. Dikes impound or partly impound ground water by preventing or retarding its movement toward discharge points. The top of this water, called high-level water in Hawaii, is at an altitude of about 1,000 feet in the north end of windward Oahu and 400 feet near the south end in Waimanalo Valley. It underlies most of the area and extends near or to the surface in poorly permeable rocks in low-lying areas. Permeability is high in less weathered mountain areas and is highest farthest away from the dike complex. Ground-water storage fluctuates to some degree owing to limited changes in the level of the ground-water reservoir--maximum storage is about 60,000 million gallons. The fluctuations control the rate at which ground water discharges. Even at its lowest recorded level, the reservoir contains a major part of the storage capacity because most of the area is perennially saturated to or near the surface. Tunnels have reduced storage by about 26,000 million gallons--only a fraction of the total storage--by breaching dike controls. Much of the reduction in storage can be restored if the .breached dike controls are replaced by flow-regulating bulkheads. Perennial streams intersect high-level water and collectively form its principal discharge. The larger streams are those that cut deepest into high-level reservoirs. Except near the coast in the northern end of the area, where dikes are absent, total base flow of streams equals total ground-water discharge. Development of high-level water by tunnels and wells diverts ground-water discharge from streams, decreasing the base flow of these streams. Construction of Haiku tunnel decreased the flow of Kahaluu Stream, 2 ? miles away, by about 26 percent. The dependable flow of water is estimated at 118 mgd (million gallons per day), of which 84 mgd is discharged by streams, tunnels, springs, and wells The remaining 34 mgd is underflow, most of it discharging into the sea near the northern end of ,the area. Average flow is estimated at 220 mgd, of which 159 mgd is. inventoried flow and 61 mgd is estimated underflow. Specific capacity of wells tapping lava flows of the Koolau Volcanic Series ranges from less than 1 to 11 gallons per minute per foot of drawdown in the dike-complex zone and from 2 to 100 in the marginal dike zone. A transmissivity of 4,000,000 gallons per day per foot was determined for the basal aquifer. Permeabilities of rocks in high mountainous areas penetrated by water-development tunnels were compared by recession constants determined from free-flow drainage. Evapotranspiration was estimated from regression curves obtained by correlating median annual rainfall and median annual pan evaporation. Evapotranspiration values from these curves compared favorably w4th values obtained from water-budget listings of rainfall and measured ground-water flow. The chemical quality of water in wells and tunnels tapping rocks of the Koolau and Honolulu Volcanic Series is excellent. Except in a few isolated areas near the shore, the chloride content of the water from these sources is generally less than 100 parts per million. Wells tapping calcareous materials are subject to sea-water contamination under heavy pumping.

Water Supply Paper↗
USGS2006

Costimulatory receptors in a teleost fish: Typical CD28, elusive CTLA4

T cell activation requires both specific recognition of the peptide-MHC complex by the TCR and additional signals delivered by costimulatory receptors. We have identified rainbow trout sequences similar to CD28 (rbtCD28) and CTLA4 (rbtCTLA4). rbtCD28 and rbtCTLA4 are composed of an extracellular Ig-superfamily V domain, a transmembrane region, and a cytoplasmic tail. The presence of a conserved ligand binding site within the V domain of both molecules suggests that these receptors likely recognize the fish homologues of the B7 family. The mRNA expression pattern of rbtCD28 and rbtCTLA4 in naive trout is reminiscent to that reported in humans and mice, because rbtCTLA4 expression within trout leukocytes was quickly up-regulated following PHA stimulation and virus infection. The cytoplasmic tail of rbtCD28 possesses a typical motif that is conserved in mammalian costimulatory receptors for signaling purposes. A chimeric receptor made of the extracellular domain of human CD28 fused to the cytoplasmic tail of rbtCD28 promoted TCR-induced IL-2 production in a human T cell line, indicating that rbtCD28 is indeed a positive costimulator. The cytoplasmic tail of rbtCTLA4 lacked obvious signaling motifs and accordingly failed to signal when fused to the huCD28 extracellular domain. Interestingly, rbtCTLA4 and rbtCD28 are not positioned on the same chromosome and thus do not belong to a unique costimulatory cluster as in mammals. Finally, our results raise questions about the origin and evolution of positive and negative costimulation in vertebrate immune systems. T cell activation is initiated through complex cell-to-cell interactions. TCRs expressed on the surface of T cells first recognize antigenic peptides presented by MHC molecules on the surface of APCs. TCRs and other surface receptors and ligands stabilize the contact between T cells and APCs, which triggers signal transduction pathways resulting in T cell activation. According to the two-signal model ( 1 , 2 ), T cell activation by Ag requires both specific recognition of the peptide by the TCR (signal 1) and additional signals delivered by other costimulatory receptors (signal 2). Among the known costimulatory receptors, CD28 and CTLA4 (also known as CD152), which are expressed on T cells and bind their corresponding ligands B7-1/B7-2 (CD80/CD86) on APCs, represent a well-studied system in mammals ( 3 ). CD28 is a glycosylated homodimeric protein expressed on the surface of double-positive thymocytes, mature CD4 + T cells, CD11b − CD8 + T cells, and γδ CD3 + T cells ( 4 , 5 ). CD28 is expressed at higher levels on activated T cells than on resting cells ( 6 ). In mammals, CD28 binds to B7-1 and B7-2, which are expressed on the surface of APCs, and delivers a critical costimulatory signal to T lymphocytes. In the absence of CD28 ligation, TCR binding either induces apoptosis or anergy in T cells ( 5 ). Engagement of CD28 alone cannot activate T cells, even if synergetic effects are provided via stimulation by T cell mitogens or anti-CD3 treatment ( 7 ). A tyrosine-based motif in cytoplasmic tail of CD28 functions as a binding site for the p85 PI3K subunit when it is phosphorylated ( 8 , 9 ). Other CD28 intracytoplasmic motifs, which associates with IL-2-inducible tyrosine kinase, lymphocyte-specific tyrosine kinase, and the adaptor growth-factor receptor bound protein-2 are responsible for signal transduction ( 10 , 11 , 12 ). The binding of growth-factor receptor bound protein-2 to CD28 also activates the GTPase RAS ( 13 ). CD28 signaling is also thought to contribute to the mobilization of lipid rafts to the immunological synapse ( 14 ), the region of contact between T cell and APC, which lowers the overall threshold of TCR engagement required for effective cytokine production and proliferation ( 15 ). In contrast, CTLA4 is a powerful negative regulator of T cell activation and was first cloned via differential screening of a cytotoxic T cell cDNA library ( 16 ). CD28 and CTLA4 belong to the same family and share high sequence similarity to each other ( 16 , 17 ). Both receptors interact with B7-1 and/or B7-2 but induce different signals: CD28 amplifies signaling triggered by the TCR-CD3 complex, whereas CTLA4 generates negative signals that inhibit T cell activation ( 18 ). CTLA4 has higher affinity for B7-1 and B7-2 compared with CD28, and its expression is also induced by TCR engagement ( 19 ). The mechanisms of CTLA4-mediated suppression likely involve both competition with CD28 for B7-1/B7-2 binding and potent inhibitory signals delivered by CTLA4 ( 20 ). The inhibitory signaling mechanisms triggered by CTLA4 are not clear. The CTLA4 cytoplasmic tail lacks typical ITIM motifs and shares with CD28 a p85 binding site ( 21 ). Thus, the extent of an immune response is likely controlled through the finely tuned expression of costimulatory receptors from the CD28/CTLA4 family on the surface of activated T cells ( 22 , 23 ). Very little is known about T cell responses and activation in bony fish. TCR-αβ cDNAs ( 24 , 25 , 26 , 27 ) and polymorphic class IA and class IIA/B MHC molecules ( 28 , 29 , 30 ) have been reported in several species, suggesting that the fish TCR recognizes antigenic peptides presented by MHC molecules similar to that of mammals. T cell-mediated responses in fish are supported by several lines of evidence. Allograft rejection provided the first experimental indications suggesting that bony fish possess a functional T cell-mediated immune response ( 31 , 32 ). In vitro assays for allospecific cytotoxicity have also been established using clonal and nonclonal catfish cell lines ( 33 , 34 ) or clonal rainbow trout ( 35 ). Autologous cell-mediated specific lysis of virus-infected syngenic target cells has also been reported in cloned goldfish ( 36 , 37 ). Finally, both public and private T cell-specific expansions have been observed in rainbow trout during secondary immune responses to viral hemorrhagic septicemia virus (VHSV), 3 a fish rhabdovirus, using TCR-β CDR3-length spectratyping ( 38 ). TCR-MHC-peptide interactions induce signaling events through CD3 in mammals and most likely in lower vertebrates because a complete set of CD3 genes, including TCR ζ, CD3 ε, and CD3 γδ, has been found in Pufferfish and Xenopus ( 39 , 40 ). CD8α has been identified in rainbow trout ( 41 ); however, the cytoplasmic tail lacks the consensus p56 lck motif, suggesting that TCR/CD3/CD8-mediated signaling events may be different in teleosts. In this study, we have identified and initially characterized two members of the CD28 family in rainbow trout, rbtCD28 and rbtCTLA4. Their sequence features and expression patterns suggest that they are the likely homologues of mammalian CD28 and CTLA4. The potential costimulatory capacities of these fish receptors were investigated using a human T cell line expressing chimeric receptors composed of the extracellular domain of human CD28 (hCD28) fused to the cytoplasmic tail of rbtCD28 or rbtCTLA4. The chimeric hCD28-rbtCD28 receptor mediated enhanced TCR-induced IL-2 production, suggesting a costimulatory function for rbtCD28. In contrast, the divergent cytoplasmic tail from rbtCTLA4 did not mediate similar signaling activities. This study therefore provides the first characterization of costimulatory receptors in lower vertebrates.

Journal of Immunology↗