USGS ScienceSearch

SEARCH · USGS Science

Results for “Climate of the Past”

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.

1,687 records · Page 29Linked to original sources

Low-sulfidation epithermal deposits of the central Basin and Range Province, USA

The Basin and Range Province is host to many important low-sulfidation epithermal deposits. Within this broad zone of extension, epithermal deposits are hosted by specific areas of Miocene and younger bimodal volcanism. In northern Nevada, rifting and related volcanic activity occurred in response to thermal bulging during the development of the Yellowstone hotspot. The Colorado River Extensional Corridor in southern Nevada, adjoining eastern California and northwestern Arizona, is a major zone of Miocene crustal extension formed during the transformation of the western margin of North America from a convergent to a transform plate boundary. The style of mineralization in the low-sulfidation epithermal deposits in these areas is strongly controlled by the nature of the volcanic successions. Highgrade, low-tonnage deposits are commonly found in flow-dominated volcanic successions where ore deposition occurred through short-lived periods of fluid flashing. In contrast, low-grade, large-tonnage deposits are more commonly located in clastic-dominated successions where fluid infiltration of the permeable hosts and cooling predominated.

Arizona, California, Nevada, Oregon

Hydrogeologic conceptual model of groundwater occurrence and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River is a major concern in the southwestern United States where the river provides water to about 40 million people for municipal and industrial use and is used to irrigate about 5.5 million acres of land. Much of the salinity in the Colorado River Basin is derived from natural interactions of surface water and groundwater with various geologic materials (rocks, soils, and alluvial deposits). The Dolores River in southwest Colorado is a major tributary of the Colorado River that historically accounts for about 6 percent of the salinity load to the Upper Colorado River Basin with the Paradox Valley being the primary source of salinity to the Dolores River. The Paradox Valley, one of several salt-anticline valleys in the region, is a fault-bounded topographic basin aligned with and exposing an underlying salt-anticline core. Salt deposits in the Pennsylvanian Paradox Formation of the Hermosa Group form an elongated salt diapir oriented northwest to southeast that is up to 12,000 feet (ft) thick beneath the present valley floor. Surface erosion, groundwater circulation, and weathering during Tertiary and Quaternary valley formation contributed to development of a cap rock, collapse features, breccia, and brine at the top of the exposed salt diapir. Today (2023), brine occurring in the brecciated cap rock and underlying salt deposits is in hydraulic connection with an overlying freshwater alluvial aquifer, and depending on seasonal river stage and hydrologic conditions, the brine discharges to the Dolores River causing the observed increase in salinity as the river crosses the Paradox Valley. To reduce salinity concentrations in the Dolores River, the Bureau of Reclamation (Reclamation) operates the Paradox Valley Unit (PVU). The PVU project consists of nine shallow brine pumping wells near the Dolores River and one deep disposal well where the brine is injected for disposal. When operational, the PVU pumping wells extract brine from the base of the alluvial aquifer that is piped and injected into a deep disposal well about 3 miles southwest of the PVU. The PVU became fully operational July 1, 1996, and by 2015, operation of the PVU had reduced salinity concentrations in the Dolores River by as much as 70 percent compared to pre-PVU conditions. In response to a 4.5 magnitude earthquake, injection operations, and thus PVU pumping, were ceased from March 2019 to June 2022. A trial period of PVU operation began in June 2022 with a reduced injection rate, and thus PVU pumping rate, of about two-thirds capacity to gather additional information and guide future operational decisions. In cooperation with Reclamation, the U.S. Geological Survey (USGS) developed this report to present the current (2023) understanding of groundwater and brine occurrence and discharge to the Dolores River in the Paradox Valley. Results from the compilation of spatial datasets, groundwater sampling and age dating, and aquifer tests are presented to provide improved understanding of the Paradox Valley hydrogeology, to supply datasets for a numerical groundwater-flow and brine-transport model, and to support future operations of the PVU. The hydrogeologic data provided herein, along with the most recent loading analysis for the Dolores River in the Paradox Valley, and a previous conceptual model for brine discharge to the river are used to present a conceptual understanding of groundwater occurrence in the Paradox Valley.

Colorado

Ground motions from finite faults with unknown geometry

Earthquake emergency response and loss modeling require accurate models of ground‐motion intensity metrics (IMs), such as those produced by ShakeMap. In the critical minutes and hours after an earthquake, however, knowledge of the earthquake source geometry limits the accuracy of these models, especially in areas where seismic instrumentation is insufficient to constrain the shaking intensity pattern. Previous attempts to overcome this limitation have centered on computing the mean distance metrics from a suite of assumed source geometries. That approach is inadequate with modern ground‐motion models (GMMs) because they typically include terms for the hanging wall effect, depth to top of rupture, and other nonlinear terms. Such approaches also propagate the distance uncertainty through the GMMs and may not properly allocate the additional uncertainty into the between‐ and within‐event components. Here we present Finite‐Fault SIMulation with Mixed Effects Regression (FFSIMMER), a method to address these limitations by computing IM distributions from a suite of rupture realizations that fully accounts for all geometry‐related GMM terms and allows for the separation of the between‐ and within‐event components of the total uncertainty. The accuracy of this method is limited by the available constraints on the source parameters. To address this, FFSIMMER can incorporate constraints based on knowledge of the source, such as the approximate strike and dip based on knowledge of the plate tectonics and focal mechanism of the earthquake. We demonstrate the improvement in accuracy that can be achieved with the 2025 M 7.7 Mandalay, Burma (Myanmar), earthquake by incorporating constraints on the range of plausible source parameters.

Bulletin of the Seismological Society of America

Three-dimensional hydrogeologic framework model of the Russian River watershed, California

The Russian River watershed is in northern Sonoma County and southern Mendocino County, California, in the northern part of the California Coast Ranges. The Russian River serves as a supply for agricultural irrigation and for municipal, domestic, and commercial uses. Through a cooperative agreement with the California State Water Resources Control Board and Sonoma County Water Agency, the U.S. Geological Survey has completed studies to better understand the hydrogeologic system and develop numerical hydrologic modeling tools to evaluate and aid in managing groundwater resources. This report focuses on the development of a digital three-dimensional hydrogeologic framework model of the Russian River watershed for use in groundwater resource assessment and numerical models. The digital three-dimensional hydrogeologic framework model of the Russian River watershed portrays the altitude, thickness, and extent of five hydrogeologic units. These five hydrogeologic units include (1) a basement unit, (2) the Sonoma Volcanics, (3) a consolidated sedimentary rock unit, (4) an unconsolidated sediment unit, and (5) channel alluvium. Model input data were compiled from published geologic maps, interpreted well data, and a model of the top of basement derived from gravity data. These data were used to construct surfaces that represent the upper and lower subsurface boundaries of each hydrogeologic unit. Top surfaces were created for the five hydrogeologic units and then stacked in three dimensions to create a solid-volume digital model. The digital three-dimensional hydrogeologic framework model described in this report and the corresponding data represent the generalized geometry of the subsurface geologic units; the model reproduces the input geologic data with reasonable accuracy and is consistent with previously published subsurface conceptualizations of the region. The model indicates the overall geometry of the basement within the watershed and the spatial extent, altitude, and thickness of the basin-filling units. The hydrogeologic framework model is at a scale and resolution appropriate for use as the foundation for a numerical hydrologic model of the study area.

California

High‐resolution surface deformation and slip distribution observations for the 2023 Kahramanmaraş, Türkiye, earthquake sequence help constrain the rupture process

Splay, or branch, faults are a common geometric feature of earthquake surface ruptures and may provide constraints on the rupture behavior of an earthquake. The 2023 M w 7.8 Pazarcık and M w 7.5 Elbistan, Türkiye, earthquakes are examples of ruptures with multiple small splays, and the Pazarcık earthquake nucleated on a splay fault, the Narlı fault, before rupturing bilaterally on the East Anatolian fault (EAF). Here, we present 3‐m‐resolution surface displacement from subpixel correlation of Planet Dove optical images for the entirety of both ruptures with corresponding surface slip distributions. For a 30‐km‐long study region spanning the Narlı‐EAF intersection, we compare surface slip derived from five data sets with different resolutions (on‐the‐ground, WorldView, Planet Dove, Sentinel‐2, and Sentinel‐1) to elucidate complementary information. In addition, we integrate information from the surface expression of faulting with published dynamic rupture simulations and rupture process studies to constrain a rupture evolution for the Pazarcık earthquake that is consistent across data sets. This work highlights the complementary nature of disparate surface slip data sets and the role that high‐resolution surface displacement information, including from fault splays, can play in constraining nonunique rupture models and refining understanding of the earthquake rupture process.

Bulletin of the Seismological Society of America

Estimation of reservoir storage capacity and geomorphic change detection analysis from a multibeam bathymetric survey of Randy Poynter Lake, Rockdale County, Georgia

Rockdale County Department of Water Resources has a directive to update estimates of the reservoir storage capacity of Randy Poynter Lake, located in northern Georgia, and to assess recent sedimentation and associated storage capacity loss. In 2022, the U.S. Geological Survey completed a multibeam bathymetric survey of Randy Poynter Lake to update storage capacity estimates and to quantify storage capacity change since the first multibeam bathymetric survey in 2012 in consideration of estimated errors inherent to bathymetric surveys. Data from the 2022 survey were used to generate contours of the reservoir as well as compute storage capacity at regular increments of water-surface elevation. Storage capacity comparisons between 2012 and 2022 at Randy Poynter Lake show minimal changes that are within the estimated uncertainties, with consistent or slightly increased storage capacities observed at most water-surface elevations and reductions observed at the remaining few elevations. Comparison of the multibeam bathymetric data collected in 2012 with data collected in 2022 further allowed for a formal geomorphic change detection analysis to map, quantify, and infer causation of morphological change over time with respect to a level of detectable change. The volume change in Randy Poynter Lake for the decade between 2012 and 2022 was slightly net-depositional and within the estimated uncertainty. The spatial distribution of sediment deposition was primarily concentrated in the northern portion of the lake, where the principal tributary flows into Randy Poynter Lake. The results of the geomorphic change analysis were used to further understand the future implications to storage capacity change. Despite the challenges of confirming systematic biases because of uncertainties exceeding the observed changes, insights from the study help predict long-term reservoir sediment accumulation, indicating a reservoir half-life extending about 650 years from 2022 on the basis of the current sediment yield estimates.

Georgia

Natural source zone depletion of crude oil in the subsurface: Processes controlling mass losses of individual compounds

At many petroleum hydrocarbon spill sites, residual spilled product forms a long-term source of groundwater contamination. The phrase source zone natural depletion is used to refer to the mass loss rates. Overall mass lost under environmental conditions was analyzed using conservative biomarker concentrations for a 1979 oil spill in northern Minnesota, USA. After 40–41 years, an average of 50% of the mass was lost with values ranging from 22% to 57% depending on location. It is also important to understand the composition changes in the source. To understand controls on the losses of individual compounds, concentrations of volatile hydrocarbons in oil samples were compared with aqueous solubilities, and pore-space oil saturations. The results of the comparison show that losses of the oil compounds were controlled by pore-space oil saturations, solubility, and susceptibility to degradation under methanogenic conditions. Compounds that degrade under methanogenic conditions, including toluene, o -xylene, and n -alkanes are more depleted compared to benzene, ethylbenzene, and m - and p -xylene for which losses are dominated by dissolution. These rates and compound-specific behaviors form a foundation for improved modeling approaches and risk analyses.

Minnesota

Water use permits as of July 2024 and reported water use near the North Unit of Theodore Roosevelt National Park, North Dakota, 1980–2023

Starting in the early 2000s, increasing oil and gas development in western North Dakota created a need for additional water resources from surface-water and groundwater sources near the North Unit of Theodore Roosevelt National Park. To summarize the use of water in that area, the U.S. Geological Survey, in cooperation with the National Park Service, developed a map of surface-water and groundwater resources, aquifers, and water-use diversions, and plotted water-use trends from 1980 to 2023. Reported water used from permits in the map area has more than doubled since 2020, increasing from about 750 acre-feet in 2020 to about 2,300 acre-feet in 2022 and 2,000 acre-feet in 2023. Surface water provided the primary source of reported water used for the study period with an average of about 410 acre-feet per year from 1980 through 2017 and about 1,330 acre-feet per year from 2018 through 2023. After 2011, groundwater sourced from the Little Missouri River, Tobacco Garden Creek, Fox Hills, Fort Union, and Dakota aquifers became a larger portion of total annual reported water use from permits in the map area. From 1980 through 2015, water use for irrigation averaged 86 percent of the total annual reported surface-water and groundwater use in the map area. Starting in 2011, however, industrial uses became a proportionally larger total use of water, and in 2015, became the highest reported volume of water use in the map area. From 2011 to 2023, industrial use designated for water depots increased from 50 acre-feet to about 1,370 acre-feet, accounting for about 70 percent of total reported water use in the map area in 2023.

North Dakota

Miocene evolution of the Humboldt Current

Diatom records from the East Pisco Basin (EPB) of southern Peru and offshore Ocean Drilling Program (ODP) Hole 682 A reveal stepwise increases in the primary productivity of the Humboldt Current during the middle and late Miocene. Although diatoms are present back through the late middle Eocene, successively enhanced diatom production occurs during the Miocene in four steps. The first step between 14.2 and 13.8 Ma marks the onset of diatom deposition in the Pisco-0 sequence. A second step at ∼12.9 Ma coincides with a major drop in global sea level during which diatom deposition ceased in the EPB but continued in ODP 682 A as well as in Chile, Ecuador, and Colombia. Beginning at 10.4 to 10.0 Ma, a major increase in diatom sedimentation rates in both the EPB and in ODP 682 A signals an abrupt intensification of the Humboldt Current productivity, which may be related to the closure of the Central American Seaway to deeper water circulation. A fourth step coincides with the transgressive deposition of the diatom-rich Pisco-2 sequence at 8.4 Ma and is enhanced by the onset of the global Late Miocene Biogenic Bloom (8.0–4.5 Ma). Between 7 and 6 Ma, common subtropical diatoms in ODP 682 A contrast with high diatom deposition rates in the EPB, suggesting variable El Niño-La Niña conditions. During this time, fossil vegetation evidence from southern Peru supports the presence of seasonal periods of enhanced rainfall typical of El Niño conditions. The Miocene trend of increasing sea mammal diversity in the EPB parallels these steps, further supporting stepwise enhancement of primary productivity in the Humboldt Current during the Miocene.

Humboldt Current

Three-dimensional geologic framework model of the Rio San Jose groundwater basin and adjacent areas, New Mexico

As part of a U.S. Geological Survey study in cooperation with the Bureau of Reclamation and the Pueblo of Acoma, New Mexico, and the Pueblo of Laguna, New Mexico, a digital three-dimensional geologic framework model was constructed for the Rio San Jose and its surface-water drainage basin in west-central New Mexico. This three-dimensional model defines the altitude, thickness, and extent of 18 geologic units for use in a regional numerical hydrologic model. The model included an undifferentiated Proterozoic basement layer, 13 consolidated Paleozoic and Mesozoic rock units, and 4 Cenozoic units. Model input data were compiled from published cross sections, well data, structure contour maps, selected geophysical data, and data derived from geologic maps and structural features in the study area. These data were used to construct faulted surfaces that represent the upper and lower subsurface geologic unit boundaries. The digital three-dimensional geologic framework model combines faults, the altitude of the tops of each geologic unit, and boundary lines depicting the subsurface extent of each geologic unit. The digital three-dimensional geologic model described in this report and the corresponding data release represents the generalized geometry of the subsurface geologic units; it reproduces with reasonable accuracy the input geologic data and is consistent with previously published subsurface conceptualizations of the region. The geologic framework model is at a scale and resolution appropriate for use as the foundation for a numerical hydrologic model of the study area.

New Mexico

Getting to the surface of the problem: A dynamic rupture benchmark for shallowly dipping faults near Earth’s surface

The Statewide California Earthquake Center (SCEC) Dynamic Rupture Group workshop “Getting to the Surface of the Problem: A Dynamic Rupture Benchmark for Shallowly-Dipping Faults Near Earth’s Surface” was convened virtually, in an online meeting on November 4, 2024. 52 people participated, including two who sent pre-recorded lightning talks earlier, but were unable to join us that day due to the international time-zone difference. One-half of our participants were students or postdocs, and the group included scientists from 11 countries (USA, Canada, China, Czech Republic, France, Germany, Iceland, Japan, New Zealand, Saudi Arabia, Switzerland).

Final Report

Status and performance of the ShakeAlert® earthquake early warning system: 2019-2023

The U.S. Geological Survey (USGS)‐operated ShakeAlert® system is the United States West Coast earthquake early warning system ( Given et al. , 2018 ). In this study we detail ShakeAlert’s performance during some of the largest events seen by the system thus far. Statewide public alerting using ShakeAlert messages was authorized in California in October 2019. Over the next few years, public alerts were expanded into Oregon and then into Washington ( U.S. Geological Survey, 2024 ). ShakeAlert source results are routinely compared to the USGS Comprehensive Catalog (ComCat; Guy et al. , 2015 ; U.S. Geological Survey, Earthquake Hazards Program, 2017 ), which contains the earthquake location and magnitude determined using complete waveform data. M 4.5 and larger is the threshold used for public alerting and was deliberately set below the level where damage is likely to compensate for cases where the system underestimates the magnitude. Between 17 October 2019 and 1 September 2023, the ShakeAlert system created 95 events with maximum magnitude estimates of M ≥4.5, the public alerting threshold. 94 of the 95 events were due to real earthquakes. Seven were categorized “false” per ShakeAlert’s internal definition that there was no matching catalog event within 100 km and 30 s of origin time; however, all but one of these were real earthquakes that were poorly located, primarily because they were at the edges of the seismic network. Three detected events were labeled “missed” because they were very poorly located (>100 km location error). In addition, the system did not produce solutions for four ComCat events M ≥4.5 ( U.S. Geological Survey, Earthquake Hazards Program, 2017 ), which were all at the edge of the alerting and network boundaries. The ShakeAlert system has accurately detected the majority of earthquakes that have occurred within the operational region since completing the public rollout, and alerts from the system have been delivered to millions of cell phone users throughout the West Coast.

California, Oregon, Washington

Geophysical constraints on continental rejuvenation in central China: Implications for outward growth of the Tibetan Plateau

Continental rejuvenation results from the tectonic reactivation of crustal structures and lithospheric reworking by mantle flow. Geochemical observations and field mapping have traditionally provided the primary evidence for the secular evolution of crustal composition and tectonic processes during continental rejuvenation. Nonetheless, the impact of continental rejuvenation on the observed present-day strain rate and orogenic-scale lithospheric structure has not been well constrained. The pre-existing E-W–trending Central China Orogenic Belt has been overprinted by the N-S–trending Central Longitudinal Seismic Belt and constitutes the intracontinental West Qinling Syntaxis in central China, where the tectonic setting changes eastward from contraction to extension. Combining updated global positioning system data and high-resolution crustal seismic tomography, we reveal a modern continental rejuvenation process within the West Qinling Syntaxis in central China. The northward extrusion of the Tibetan Plateau's weak lithospheric layer (middle-lower crust and lithospheric mantle) of southwestern China relative to the rigid Sichuan Basin/Ordos Block of the eastern West Qinling Syntaxis results in regional dextral shearing that shapes the Central Longitudinal Seismic Belt and defines the eastern Tibetan Plateau margin. The pre-existing E-W–trending Central China Orogenic Belt has been preserved above the brittle-ductile transition zone, and the northward movement of the deep lithospheric layer drives the deformation of the upper crust in the West Qinling Syntaxis. Our results, along with previous studies, suggest the presence of an intracontinental lithospheric interchange structure in central China. The continental rejuvenation of the West Qinling Syntaxis results from a combination of fault reactivation in the upper crust (Stage I, Eocene–Oligocene) and reworking of the deep lithosphere (Stage II, middle–late Miocene) related to the plateau-wide shift in stress accommodation ultimately driven by the redistribution of mass outward from the central Tibetan Plateau. At present, the transition zone between the high- and low-velocity anomalies along the Central Longitudinal Seismic Belt not only shapes the landscape boundary but controls the size and recurrence interval of earthquakes within the West Qinling Syntaxis in central China.

GSA Bulletin

Estimation of dynamic geologic CO2 storage resources in the Illinois Basin, including effects of brine extraction, anisotropy, and hydrogeologic heterogeneity

Since the vast majority of carbon dioxide (CO 2 ) storage resources in the United States are in deep saline aquifers, optimizing the use of these saline storage resources could be crucial for efficient development of geologic CO 2 storage (GCS) resources and basin- or larger-scale deployment of GCS in the country. Maximum CO 2 injection rates can be enhanced by extracting brine from the CO 2 storage unit. However, disposal of the extracted brine is both a technological and economic challenge. The lowest-cost option would likely be reinjection of the extracted brine into another formation above or below the CO 2 storage unit. Therefore, it is important to estimate brine injectivity as it will constrain the potential to increase CO 2 injectivity at an injection site that has access to multiple geologic storage units where either CO 2 or brine can be injected. Using a simulation-optimization framework, coupled with a non-isothermal, multiphase CO 2 -water-salt equation-of-state module, we developed a computationally efficient method for evaluating optimization of simultaneous CO 2 injection, brine extraction, and brine (re)injection at hypothetical injection sites deployed across a geologic basin. The Illinois basin is ideal for testing our methodology because it contains multiple geologic storage units with seals in between them to isolate injection of CO 2 in one unit from interfering with the injection of either brine or CO 2 in another unit above or below it. In addition, we investigated the relative effects of variation in key geologic parameters as well as two reservoir structures (hydrogeologic heterogeneity/anisotropy and homogeneity/isotropy) on CO 2 injectivities and enhancement of CO 2 injectivity through extracting brine. Results suggest that permeability, depth, and especially thickness of the storage unit could be the most influential parameters determining CO 2 injectivity. They also suggest that only injecting CO 2 into the storage unit with the greatest injectivity, enhancing that unit’s injectivity by extracting brine, and disposing of the produced brine in other suitable units could maximize total CO 2 injectivity in limited regions of the basin. At the majority of simulated injection sites, however, we found that injecting CO 2 into all of the accessible and suitable storage units was more likely to maximize the CO 2 storage resource.

Illinois, Kentucky, Indiana

Geochemical evidence for the origin of late Quaternary loess, Seward Peninsula, Alaska

Loess is the most widespread surficial deposit in the state of Alaska. Although loess of last glacial age is common in mid-continental North America, records of last glacial loess in Alaska have been elusive. Here we report a record of last glacial loess on the Seward Peninsula, along with studies of a transect of loess localities also thought to date to the last glacial period. Highly diverse hypotheses on the origin of loess on the Seward Peninsula include local volcanic rocks and mountain ranges to the north and south, as well as offshore deposits. We present particle size and mineralogical and geochemical data that indicate that the most likely sources for loess in this part of eastern Beringia are glaciogenic silt particles of the Kobuk River and Noatak River, both of which drain the Brooks Range to the northeast. Northeasterly paleowinds were therefore responsible for loess on the Seward Peninsula. Such winds probably occurred during late summer or early autumn, when river discharge was at a minimum but temperatures were still above freezing. Northeasterly winds during the last glacial period are in good agreement with a recent paleoclimate model.

Alaska

Estimating the importance of floating surface material to the total phosphorus transport in Silver Creek, Wisconsin using Particle Image Velocimetry

Various techniques are used to estimate nutrient delivery in streams that combine flow and water-quality data. However, the transport of surface floating material is difficult to measure, and is therefore typically neglected when stream sampling and in the estimated nutrient delivery. Here, we describe an approach to estimate the amount of material (duckweed ( Lemna genus), filamentous algae, and other macrophyte fragments) and associated nutrients (in this case, phosphorus, P) transported on the surface of Silver Creek, Wisconsin, to determine if this material is an important transport mechanism and if historical P loads were underestimated. This approach includes estimating the transport of surface material using 10 s videos collected every 15 min from a downward-looking camera installed beneath a bridge. The average velocity of the surface material was first determined using Large-Scale Particle Image Velocimetry (LSPIV), which uses short videos to analyze surface particle movement. The amount of surface material in each video was then computed using computer-vision techniques. The P load associated with the transported surface material was then estimated by combining surface velocities, coverage of floating material, and laboratory-measured P content. Surface material transported ~9–11% of the total summer P load and ~4–7% of the annual load in Silver Creek.

Wisconsin

Groundwater and surface-water interactions in the Waihe‘e-Kahalu‘u watershed, O‘ahu, Hawai‘i—Analysis of historical data and numerical groundwater-model simulations

The Waihe‘e-Kahaluʻu watershed and surrounding areas lie in a dike-impounded groundwater setting of the Ko‘olau Range, O‘ahu, Hawai‘i, where groundwater withdrawal from wells and tunnels have resulted in reductions in streamflow. Resource managers and surface-water users seek information that can be used to balance the needs of competing uses of groundwater and streamflow in the watershed. This study assesses the effects of four groundwater developments in the Waihe‘e-Kahalu‘u watershed—Kahaluu Tunnel, Waihee Tunnel, Kahaluu Well, and Waihee Incline Wells—on Kahalu‘u and Waihe‘e Streams. Analyses of historical streamflow, withdrawal, and rainfall data indicate that base flow (the component of streamflow that is supplied from groundwater discharge) in Kahalu‘u and Waihe‘e Streams decreased in response to withdrawals from wells and tunnels within and outside the watershed. Groundwater-model simulations indicate that shutdown of selected tunnels and wells in the Waihe‘e-Kahaluʻu watershed will result in increases of groundwater discharge, but not all of the increase will occur in streams within the watershed; increases will also occur in streams outside the basin, the ocean, and tunnels other than those that were shut down. The analysis of historical data and model simulations is consistent in showing that the effects of withdrawal changes on streams are not constrained by the topographic boundaries of watersheds. Stream base flows in the Waihe‘e-Kahalu‘u watershed are subject to a complex interaction of multiple withdrawals from wells and tunnels installed at various times within and outside the watershed. The correspondence between withdrawals and base-flow changes is consistent with the conceptual model and understanding of the connection between groundwater and streamflow in the dike-impounded groundwater setting of the Waihe‘e-Kahalu‘u watershed and the Ko‘olau Range.

Hawaii

A newly identified creeping strand of the Concord fault, San Francisco Bay Area

The Concord fault constitutes a major branch of the Pacific–North America transform plate boundary in Northern California, bridging the strike‐slip Bartlett Springs ‐ Green Valley Fault system to the north with the Greenville and Calaveras Faults to the south. Like many faults in the San Francisco Bay Area its long‐term slip is partially accommodated by aseismic slip (creep). Although creep has been recognized and monitored on the northern half of the fault for decades, the precise location of the southern half of the fault and its slip rate—whether accommodated seismically or aseismically—has remained enigmatic. How slip transfers between the Concord and Greenville or Calaveras faults to the south remains an outstanding question. New field observations presented here indicate that the active trace of the fault south of downtown Concord is not where previously interpreted and is indeed actively creeping. We report observations of shallow creep continuing >7 km farther south along the Concord fault than previously reported, along a fault strand not previously recognized for most of its length. This is evident as right‐laterally deflected concrete curbs and sidewalk slabs on both sides of every street that crosses the fault at a high angle in southeast Concord and northeast Walnut Creek. We document the magnitude and location of these deflections to estimate accumulated right‐lateral aseismic slip expressed in engineered structures. Offsets of these piercing lines range from 8 to 18 cm, over widths varying from narrow breaks along centimeter‐scale concrete joints to 10‐m‐wide zones of deflection. Significantly, this active trace is ∼400 m west of where the Quaternary active trace has previously been inferred, placing it within—rather than bounding—the built area of suburban Concord. Slip along the fault has already caused infrastructure damage. These results revise our understanding of the southern Concord fault and help constrain its seismic potential.

California