SEARCH · USGS Science

Results for “Near Surface Geophysics”

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

SourceResearch connectionFrom yearThrough year

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

At least 919 records · Page 51Linked to original sources

USGS2011

Vulnerability of high-latitude soil organic carbon in North America to disturbance

This synthesis addresses the vulnerability of the North American high-latitude soil organic carbon (SOC) pool to climate change. Disturbances caused by climate warming in arctic, subarctic, and boreal environments can result in significant redistribution of C among major reservoirs with potential global impacts. We divide the current northern high-latitude SOC pools into (1) near-surface soils where SOC is affected by seasonal freeze-thaw processes and changes in moisture status, and (2) deeper permafrost and peatland strata down to several tens of meters depth where SOC is usually not affected by short-term changes. We address key factors (permafrost, vegetation, hydrology, paleoenvironmental history) and processes (C input, storage, decomposition, and output) responsible for the formation of the large high-latitude SOC pool in North America and highlight how climate-related disturbances could alter this pool's character and size. Press disturbances of relatively slow but persistent nature such as top-down thawing of permafrost, and changes in hydrology, microbiological communities, pedological processes, and vegetation types, as well as pulse disturbances of relatively rapid and local nature such as wildfires and thermokarst, could substantially impact SOC stocks. Ongoing climate warming in the North American high-latitude region could result in crossing environmental thresholds, thereby accelerating press disturbances and increasingly triggering pulse disturbances and eventually affecting the C source/sink net character of northern high-latitude soils. Finally, we assess postdisturbance feedbacks, models, and predictions for the northern high-latitude SOC pool, and discuss data and research gaps to be addressed by future research.

Journal of Geophysical Research G: Biogeosciences↗

USGS2008

Sediment deposition, erosion, and bathymetric change in central San Francisco Bay: 1855-1979

Central San Francisco Bay is the hub of a dynamic estuarine system connecting the San Joaquin and Sacramento River Deltas, Suisun Bay, and San Pablo Bay to the Pacific Ocean and South San Francisco Bay. To understand the role that Central San Francisco Bay plays in sediment transport throughout the system, it is necessary to first determine historical changes in patterns of sediment deposition and erosion from both natural and anthropogenic forces. The first extensive hydrographic survey of Central San Francisco Bay was conducted in 1853 by the National Ocean Service (NOS) (formerly the United States Coast and Geodetic Survey (USCGS)). From 1894 to 1979, four additional surveys, composed of a total of approximately 700,000 bathymetric soundings, were collected within Central San Francisco Bay. Converting these soundings into accurate bathymetric models involved many steps. The soundings were either hand digitized directly from the original USCGS and NOS hydrographic sheets (H-sheets) or obtained digitally from the National Geophysical Data Center's (NGDC) Geophysical Data System (GEODAS) (National Geophysical Data Center, 1996). Soundings were supplemented with contours that were either taken directly from the H-sheets or added in by hand. Shorelines and marsh areas were obtained from topographic sheets. The digitized soundings, depth contours, shorelines, and marsh areas were entered into a geographic information system (GIS) and georeferenced to a common horizontal datum. Using surface modeling software, bathymetric grids with a horizontal resolution of 25 m were developed for each of the five hydrographic surveys. Before analyses of sediment deposition and erosion were conducted, interpolation bias was removed and all of the grids were converted to a common vertical datum. These bathymetric grids were then used to develop bathymetric change maps for subsequent survey periods and to determine long-term changes in deposition and erosion by calculating volumes and rates of net sediment change. Central San Francisco Bay experienced periods of both deposition and erosion, but overall experienced a net gain in sediment from 1855 to 1979 of approximately 42x10 6 m 3 (0.33x10 6 m 3 / yr). Over this same time period, 92 percent of the tidal marsh and 69 percent of the intertidal mudflats were lost as human activity increased and the shorefront was developed. During the first time period, from 1855 to 1895, Central San Francisco Bay was erosional, losing roughly 2x10 6 m 3 / yr of sediment. The next time period was depositional, with a net gain of approximately 3x10 6 m 3 / yr of sediment from 1895 to 1947. The last time period, from 1947 to 1979, was erosional again, losing roughly 2x10 6 m 3 / yr of sediment. Sedimentation patterns also varied spatially. The northern part of Central San Francisco Bay was depositional during all change periods while the eastern region alternated between erosional and depositional. Central San Francisco Bay sedimentation patterns have also been strongly impacted by anthropogenic activities, such as dredging and dredge disposal, borrow pits, and sand mining. For example, bathymetric change at a borrow pit created near Bay Farm Island sometime between the 1947 and 1979 surveys indicates roughly 25x10 6 m 3 of sediment was removed from the system

California↗

USGS2008

First in situ investigation of a dark wind streak on Mars

Prominent low‐albedo wind streaks issue from embayments at the north end of Victoria crater in Meridiani Planum, the site of surface investigations by the Mars Exploration Rover Opportunity. Opportunity was sent to examine the darkest of these wind streaks and compare it to the adjacent soil in the brighter corridor between the streaks. Two nearby sites inside and outside of the streak were selected for visible and infrared remote sensing and for elemental abundance measurements and images of the microscopic soil morphology. The rover was next sent to study the interaction of the wind streak with an obstacle, a small rock near the source of the streak. Opportunity's observations suggest that the Victoria wind streaks are deposits of basaltic sand blown out of the crater from the dark dunes nestled below the crater rim, particularly at the base of the alcove leading up to the darkest streak. No local sources of sand have been identified within the Victoria crater wall rock, suggesting that the sand is supplied from outside the crater and is presently escaping its temporary topographic trap. This process presents a possible explanation for the serrated margin of Victoria crater through abrasion of the soft rock as trapped sands are blown out of the crater and carve alcoves under various seasonal winds. Some indication of the rate of activity of these wind streaks will be provided by continued monitoring of the erasure of the rover's tracks.

Journal of Geophysical Research E: Planets↗

USGS1996

Type of faulting and orientation of stress and strain as a function of space and time in Kilauea's south flank, Hawaii

Earthquake focal mechanisms of events occurring between 1972 and 1992 in the south flank of Kilauea volcano, Hawaii, are used to infer the state of stress and strain as a function of time and space. We have determined 870 fault plane solutions from P wave first motion polarities for events with magnitudes M L ≥ 2.5 and depth ranging between 6 and 12 km. Faulting is characterized by a mixture of decollement, reverse, and normal faults. Most large earthquakes with magnitude M < 7 slip on reverse faults striking NE at 40° and dipping SE between 60° and 70°. In Hawaii, the earthquakes with M > 7 rupture the decollement plane, since it is the only surface large enough to generate magnitude 7 or larger earthquakes. The percentage of reverse faulting events is high compared to the decollement and normal faulting mechanisms for the period 1972–1983. The percentage of decollement type focal mechanisms becomes dominant after 1983. This pattern of faulting activity suggests that pressure was building up within Kilauea's rift zone prior to the 1983 Puu'Oo eruption. Overall, a single stress orientation with the maximum compressive stress oriented SE perpendicular to the rift and dipping at 45° is compatible with the coeval existence of decollement, reverse, and normal faults. However, in a crustal volume east of longitude 155°10′W, we find a change of the orientation of σ 1 from nearly horizontal to plunging 45° SE occurring in 1979. This stress rotation suggests magma movements within the aseismic part of Kilauea's east rift zone. The strain and stress orientations are coaxial in the south flank except within the volume where the stress rotation is observed. We observe a change in the relationship between stress and strain directions caused either by the shifting of seismic activity from reverse faults to decollements, while stress stays constant, or by a rotation of stress, while strain remains constant. Assuming that the model of a noncohesive Coulomb wedge is appropriate for Kilauea's south flank, we find that high pore pressures are prevalent along the decollement and within the wedge for a coefficient of friction equal to 0.85.

Hawaii↗

USGS2025

A decade of precursory seismicity at Mauna Loa: 2013–2022

The 2022 eruption of Mauna Loa was preceded by nearly 10 years of seismic unrest. Using the earthquake catalog of the Hawaiian Volcano Observatory, we describe and explore many aspects of unrest between January 1, 2013, and the start of the eruption on November 28, 2022 (UTC). Earthquakes associated with unrest occurred above an inflating dike at shallow depths under and to the south of Moku‘āweoweo and in the cluster of seismicity to the west of Moku‘āweoweo associated with a buried rift zone. Earthquakes also occurred in a broad region at the base of the crust, likely associated with a zone of magma accumulation. Unrest is broken down into four phases of increasing intensity with time and into individual geographic regions. While activity occurring at shallow depths under the caldera often alternated with activity on the west flank in the first two phases of unrest, the third phase of unrest featured increased rates of earthquake activity in both areas acting concurrently. We also document a slow shift in seismicity across multiple years from the south of Moku‘āweoweo to an area under the caldera that likely represents a change in the location of the intrusion through time. Broadly speaking, unrest observed in advance of the 2022 eruption occurred in the same areas as the unrest preceding eruptions in 1975 and 1984. Furthermore, precursory unrest prior to eruptions in 1975, 1984 and 2022 lasted years and culminated in an intense phase of seismic activity lasting on the order of hours that represented the final dike propagation to the surface. Concurrent earthquake activity under Moku‘āweoweo and the west flank of Mauna Loa, paired with a shift in earthquakes under Moku‘āweoweo were the two most notable long-term precursors to an eruption; however, neither provided an actionable timeline of an impending eruption. The decade of precursory seismicity prior to the 2022 eruption is similar to long precursory phases prior to eruptions in 1975 and 1984 and underscores the challenges of informing on eruption hazards in a timely manner. By combining seismic activity with other geophysical fields such as geodesy and gravity, we present a conceptual model to aid in the interpretation of future unrest at Mauna Loa.

Hawaii↗

USGS2013

Summit crater lake observations, and the location, chemistry, and pH of water samples near Mount Chiginagak volcano, Alaska: 2004-2012

Mount Chiginagak is a hydrothermally active volcano on the Alaska Peninsula, approximately 170 km south–southwest of King Salmon, Alaska (fig. 1). This small stratovolcano, approximately 8 km in diameter, has erupted through Tertiary to Permian sedimentary and igneous rocks (Detterman and others, 1987). The highest peak is at an elevation of 2,135 m, and the upper ~1,000 m of the volcano are covered with snow and ice. Holocene activity consists of debris avalanches, lahars, and lava flows. Pleistocene pyroclastic flows and block-and-ash flows, interlayered with andesitic lava flows, dominate the edifice rocks on the northern and western flanks. Historical reports of activity are limited and generally describe “steaming” and “smoking” (Coats, 1950; Powers, 1958). Proximal tephra collected during recent fieldwork suggests there may have been limited Holocene explosive activity that resulted in localized ash fall. A cluster of fumaroles on the north flank, at an elevation of ~1,750 m, commonly referred to as the “north flank fumarole” have been emitting gas throughout historical time (location shown in fig. 2). The only other thermal feature at the volcano is the Mother Goose hot springs located at the base of the edifice on the northwestern flank in upper Volcano Creek, at an elevation of ~160 m (fig. 2, near sites H1, H3, and H4). Sometime between November 2004 and May 2005, a ~400-m-wide, 100-m-deep lake developed in the snow- and ice-filled summit crater of the volcano (Schaefer and others, 2008). In early May 2005, an estimated 3 million cubic meters (3×106 m3) of sulfurous, clay-rich debris and acidic water exited the crater through tunnels at the base of a glacier that breaches the south crater rim. More than 27 km downstream, these acidic flood waters reached approximately 1.3 m above normal water levels and inundated a fertile, salmon-spawning drainage, acidifying the entire water column of Mother Goose Lake from its surface waters to its maximum depth of 45 m (resulting pH ~2.9), and preventing the annual salmon run in the King Salmon River. A simultaneous release of gas and acidic aerosols from the crater caused widespread vegetation damage along the flow path. Since 2005, we have been monitoring the crater lake water that continues to flow into Mother Goose Lake by collecting surface water samples for major cation and anion analysis, measuring surface-water pH of affected drainages, and photo-documenting the condition of the summit crater lake. This report describes water sampling locations, provides a table of chemistry and pH measurements, and documents the condition of the summit crater between 2004 and 2011. In September 2013, the report was updated with results of water-chemistry samples collected in 2011 and 2012, which were added as an addendum.

Alaska↗

USGS2003

Hydrogeology of shallow basin-fill deposits in areas of Salt Lake Valley, Salt Lake County, Utah

A study of recently developed residential/commercial areas of Salt Lake Valley, Utah, was done from 1999 to 2001 in areas in which shallow ground water has the potential to move to a deeper aquifer that is used for public supply. Thirty monitoring wells were drilled and sampled in 1999 as part of the study. The ground water was either under unconfined or confined conditions, depending on depth to water and the presence or absence of fine-grained deposits. The wells were completed in the shallowest water-bearing zone capable of supplying water. Monitoring-well depths range from 23 to 154 feet. Lithologic, geophysical, hydraulic-conductivity, transmissivity, water-level, and water-temperature data were obtained for or collected from the wells. Silt and clay layers noted on lithologic logs correlate with increases in electrical conductivity and natural gamma radiation shown on many of the electromagnetic-induction and natural gamma logs. Relatively large increases in electrical conductivity, determined from the electromagnetic-induction logs, with no major changes in natural gamma radiation are likely caused by increased dissolved-solids content in the ground water. Some intervals with high electrical conductivity correspond to areas in which water was present during drilling. Unconfined conditions were present at 7 of 20 monitoring wells on the west side and at 2 of 10 wells on the east side of Salt Lake Valley. Fine-grained deposits confine the ground water. Anthropogenic compounds were detected in water sampled from most of the wells, indicating a connection with the land surface. Data were collected from 20 of the monitoring wells to estimate the hydraulic conductivity and transmissivity of the shallow ground-water system. Hydraulic-conductivity values of the shallow aquifer ranged from 30 to 540 feet per day. Transmissivity values of the shallow aquifer ranged from 3 to 1,070 feet squared per day. There is a close linear relation between transmissivity determined from slug-test analysis and transmissivity estimated from specific capacity. Water-level fluctuations were measured in the 30 monitoring wells from 1999 to July 2001. Generally, water-level changes measured in wells on the west side of the valley followed a seasonal trend and wells on the east side showed less fluctuation or a gradual decline during the 2-year period. This may indicate that a larger percentage of recharge to the shallow ground-water system on the west side is from somewhat consistent seasonal sources, such as canals and unconsumed irrigation water, as compared to sources on the east side. Water levels measured in monitoring wells completed in the shallow ground-water system near large-capacity public-supply wells varied in response to ground-water withdrawals from the deeper confined aquifer. Water temperature was monitored in 23 wells. Generally, little or no change in water temperature was measured in monitoring wells with a depth to water greater than about 40 feet. The shallower the water level in the well, the greater the water-temperature change measured during the study. Comparison of water levels measured in the monitoring wells and deeper wells in the same area indicate a downward gradient on the east side of the valley. Water levels in the shallow and deeper aquifers in the secondary recharge area on the west side of the valley were similar to those on the east side. Water levels measured in the monitoring wells and nearby wells completed in the deeper aquifer indicate that the vertical gradient can change with time and stresses on the system.

Utah↗

USGS2016

Insights into shallow magmatic processes at Kīlauea Volcano, Hawaiʻi, from a multiyear continuous gravity time series

Continuous gravity data collected near the summit eruptive vent at Kīlauea Volcano, Hawaiʻi, during 2011–2015 show a strong correlation with summit-area surface deformation and the level of the lava lake within the vent over periods of days to weeks, suggesting that changes in gravity reflect variations in volcanic activity. Joint analysis of gravity and lava level time series data indicates that over the entire time period studied, the average density of the lava within the upper tens to hundreds of meters of the summit eruptive vent remained low—approximately 1000–1500 kg/m 3 . The ratio of gravity change (adjusted for Earth tides and instrumental drift) to lava level change measured over 15 day windows rose gradually over the course of 2011–2015, probably reflecting either (1) a small increase in the density of lava within the eruptive vent or (2) an increase in the volume of lava within the vent due to gradual vent enlargement. Superimposed on the overall time series were transient spikes of mass change associated with inflation and deflation of Kīlauea's summit and coincident changes in lava level. The unexpectedly strong mass variations during these episodes suggest magma flux to and from the shallow magmatic system without commensurate deformation, perhaps indicating magma accumulation within, and withdrawal from, void space—a process that might not otherwise be apparent from lava level and deformation data alone. Continuous gravity data thus provide unique insights into magmatic processes, arguing for continued application of the method at other frequently active volcanoes.

Hawaii↗

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↗

USGS2000

Results of the Imager for Mars Pathfinder windsock experiment

The Imager for Mars Pathfinder (IMP) windsock experiment measured wind speeds at three heights within 1.2 m of the Martian surface during Pathfinder landed operations. These wind data allowed direct measurement of near-surface wind profiles on Mars for the first time, including determination of aerodynamic roughness length and wind friction speeds. Winds were light during periods of windsock imaging, but data from the strongest breezes indicate aerodynamic roughness length of 3 cm at the landing site, with wind friction speeds reaching 1 m/s. Maximum wind friction speeds were about half of the threshold-of-motion friction speeds predicted for loose, fine-grained materials on smooth Martian terrain and about one third of the threshold-of-motion friction speeds predicted for the same size particles over terrain with aerodynamic roughness of 3 cm. Consistent with this, and suggesting that low wind speeds prevailed when the windsock array was not imaged and/or no particles were available for aeolian transport, no wind-related changes to the surface during mission operations have been recognized. The aerodynamic roughness length reported here implies that proposed deflation of fine particles around the landing site, or activation of duneforms seen by IMP and Sojourner, would require wind speeds >28 m/s at the Pathfinder top windsock height (or >31 m/s at the equivalent Viking wind sensor height of 1.6 m) and wind speeds >45 m/s above 10 m. These wind speeds would cause rock abrasion if a supply of durable particles were available for saltation. Previous analyses indicate that the Pathfinder landing site probably is rockier and rougher than many other plains units on Mars, so aerodynamic roughness length elsewhere probably is less than the 3-cm value reported for the Pathfinder site.

Journal of Geophysical Research E: Planets↗

USGS1982

The Giles County, Virginia, seismogenic zone-- Seismological results and geological interpretations

This paper describes and interprets a newly-recognized 40-km-long seismogenic zone, which is inferred to have been the locus of a damaging earthquake in 1897. That shock was the second largest known to have occurred in the southeastern United States (MMI VIII, m b estimated at 5.8, felt over 725,000 km 2 ). It struck Giles County in southwestern Virginia, and a recurrence would affect populous regions on and near the central Atlantic seaboard. This paper attempts to aid in evaluating that hazard by presenting and synthesizing new seismological data with geological inferences and deductions. A five-station, 60-km aperture seismic network has been in operation in the Giles County locale since early 1978. For the subsequent 3-year monitoring period, 10 microearthquakes (M < 2) have been detected. Eight of those 10 events, plus an additional 4 relocated felt earthquakes (3.2 < M < 4.1; 1959-1976) have a tabular distribution centered at Pearisburg, Virginia. That distribution is about 40 km long, 10 km wide, strikes N. 43° E., and has a nearly vertical extent of from 5 km to 25 km in depth. Thus, a Giles County seismogenic zone is defined presently by 12 earthquakes that span 4 orders of magnitude (0 < M < 4) and 2 decades of time (1959-1980). We conclude that the 1897 earthquake occurred on that seismogenic zone. From the orientation of the tabular zone, from evidence that greatest horizontal compressive stress trends east-northeasterly at seismogenic depths in and near Giles County and from sparse P-wave first-motion data, we infer that the monitored microseismicity probably occurs by right-reverse motion on the seismogenic zone, with the southeast side dropping down with respect to the northwest side. In the Giles County locale, the upper 3-6 km of the crust are Paleozoic sedimentary rocks that have moved some tens of kilometers northwest on nearly horizontal detachment faults. The above-mentioned reliable hypocenters for the region lie below the deepest likely detachment, indicating that Giles County seismicity probably has no simple relationship to surface geology. Since Precambrian time, three deformational episodes could have formed steep faults under today's surface structures, at the observed hypocentral depths. These episodes were as follows: (1) As the Iapetus Ocean (Atlantic's predecessor) opened in late Precambrian or early Paleozoic time, northeast-striking normal faults formed, probably at the inferred Iapetan continental edge in central Virginia and at least as far northwest of that locus as Giles County. (2) In late Paleozoic time, detachment faults loaded the crust with several kilometers of overthrust sedimentary rocks, perhaps forming northeast-striking thrust-load faults in a brittle analogue of isostatic depression caused by thrust masses and much lighter continental glaciers. (3) As the Atlantic Ocean opened in Mesozoic time, other northeast-striking normal faults formed on the present continental margin and inland of it. The N. 43° E.-striking seismogenic zone seems most likely to have resulted from compressional reactivation of an Iapetan normal fault, which also may have been reactivated by late Paleozoic compression and Mesozoic extension. First, the seismogenic zone probably does not occur on a thrust-load fault. The zone underlies detached structures of southern Appalachian orientations (east-northeast), but those structures are not known to be displaced where they cross the zone. Thus, if the zone occurs on a thrust-load fault, the fault and its coeval causative central Appalachian detachments would pre-date the southern Appalachian structures. That deduction contradicts stratigraphic and structural estimates of relative ages of southern and central Appalachian detachments. Second, the zone probably does not result from a Mesozoic normal fault, because known locations of Mesozoic normal faults and grabens are well to the southeast of Giles County. Not yet known is where else in the East reactivated Iapetan normal faults might generate shocks similar to that of 1897. However, our analysis enables us to suggest specific geological and geophysical investigations that may produce results useful in answering that question. Such investigations can concentrate on defining the area of probable occurrence of other Iapetan normal faults, and on determining whether the one inferred to underlie Giles County is uniquely active or is typical of others that might exist elsewhere.

Virginia↗

USGS2015

Low resistivity and permeability in actively deforming shear zones on the San Andreas Fault at SAFOD

The San Andreas Fault Observatory at Depth (SAFOD) scientific drillhole near Parkfield, California crosses the San Andreas Fault at a depth of 2.7 km. Downhole measurements and analysis of core retrieved from Phase 3 drilling reveal two narrow, actively deforming zones of smectite-clay gouge within a roughly 200 m-wide fault damage zone of sandstones, siltstones and mudstones. Here we report electrical resistivity and permeability measurements on core samples from all of these structural units at effective confining pressures up to 120 MPa. Electrical resistivity (~10 ohm-m) and permeability (10-21 to 10-22 m2) in the actively deforming zones were one to two orders of magnitude lower than the surrounding damage zone material, consistent with broader-scale observations from the downhole resistivity and seismic velocity logs. The higher porosity of the clay gouge, 2 to 8 times greater than that in the damage zone rocks, along with surface conduction were the principal factors contributing to the observed low resistivities. The high percentage of fine-grained clay in the deforming zones also greatly reduced permeability to values low enough to create a barrier to fluid flow across the fault. Together, resistivity and permeability data can be used to assess the hydrogeologic characteristics of the fault, key to understanding fault structure and strength. The low resistivities and strength measurements of the SAFOD core are consistent with observations of low resistivity clays that are often found in the principal slip zones of other active faults making resistivity logs a valuable tool for identifying these zones.

Journal of Geophysical Research↗

USGS2022

Microbially induced anaerobic oxidation of magnetite to maghemite in a hydrocarbon-contaminated aquifer

Iron mineral transformations occurring in hydrocarbon-contaminated sites are linked to the biodegradation of the hydrocarbons. At a hydrocarbon-contaminated site near Bemidji, Minnesota, USA, measurements of magnetic susceptibility (MS) are useful for monitoring the natural attenuation of hydrocarbons related to iron cycling. However, a transient MS, previously observed at the site, remains poorly understood and the iron mineral phases acting as reactants and products associated with this MS perturbation remain largely unknown. To address these unknowns, we acquired mineral magnetism measurements, including hysteresis loops, backfield curves, and isothermal remanent magnetizations on sediment core samples retrieved from the site and magnetite-filled mineral packets installed within the aquifer. Our data show that the core samples and magnetite packs display decreasing magnetization with time and that this loss in magnetization is accompanied by increasing bulk coercivity consistent with decreased average grain size and/or partial oxidation. Low-temperature magnetometry on all samples displayed behavior consistent with magnetite, but samples within the plume also show evidence of maghemitization. This interpretation is supported by the occurrence of shrinkage cracks on the surface of the grains imaged via scanning electron microscopy. Magnetite transformation to maghemite typically occurs under oxic conditions, here, we propose that maghemitization occurs within the anoxic portions of the plume via microbially mediated anaerobic oxidation. Mineral dissolution also occurs within the plume. Microorganisms capable of such anaerobic oxidation have been identified within other areas at the Bemidji site, but additional microbiological studies are needed to link specific anaerobic iron oxidizers with this loss of magnetization.

Minnesota↗

USGS2007

Hydrogeology of the Piney Point-Nanjemoy, Aquia, and Upper Patapsco aquifers, Naval Air Station Patuxent River and Webster Outlying Field, St. Marys County, Maryland, 2000–06

Recent and projected population growth in southern Maryland continues to bring ground-water-quality and quantity issues to the forefront. Lithologic, borehole geophysical, water-level, and water-use data were compiled and interpreted to revise understanding of the hydrogeologic framework of the Piney Point-Nanjemoy, Aquia, and Upper Patapsco aquifers in southern Maryland, with emphasis on the Naval Air Station Patuxent River and Webster Outlying Field. Understanding of the hydrogeologic framework for the Upper Patapsco aquifer also has been revised based on the results of aquifer testing and water-quality sampling of two wells. The Piney Point-Nanjemoy aquifer is 50 to 70 feet thick, with a top altitude of 213 to 260 feet below the North American Vertical Datum of 1988 and a hydraulic conductivity of 2 feet per day at Naval Air Station Patuxent River and Webster Outlying Field. Ground-water withdrawals from the Piney Point-Nanjemoy aquifer have been minimal since 1999 and water levels in the aquifer have not changed substantially since the 1950s. An overall decline of about 2.5 feet has been observed since 1997, however. The Aquia aquifer is 100 to 145 feet thick, with a top altitude of approximately 450 feet below the North American Vertical Datum of 1988 and a hydraulic conductivity of 6 to 10 feet per day at Naval Air Station Patuxent River. The Aquia aquifer is approximately 50 feet thick, with a top altitude of 470 feet below sea level and a hydraulic conductivity of 6 to 10 feet per day at Webster Outlying Field. Water levels in the Aquia aquifer declined in response to increased withdrawals from the aquifer from the early 1940s through about 2000 at Naval Air Station Patuxent River and Webster Outlying Field, but have been generally stable from about 1999 through April 2006. The Upper Patapsco aquifer at the Naval Air Station Patuxent River and Webster Outlying Field consists of layers of sand interbedded with layers of clay that total over 200 feet in thickness. The top of the aquifer near Naval Air Station Patuxent River and Webster Outlying Field ranges from approximately 620 to 680 feet below the North American Vertical Datum of 1988. The potentiometric surface of the Upper Patapsco aquifer near Naval Air Station Patuxent River and Webster Outlying Field has been declining at a rate of about 2 feet per year for the past several years. Aquifer testing indicates the transmissivity of the aquifer is about 2,100 to 3,900 feet squared per day. Water samples collected from the Upper Patapsco aquifer at Naval Air Station Patuxent River and Webster Outlying Field in 2000 and 2001 met all Federal drinking-water standards.

Scientific Investigations Report↗

USGS2012

Shallow structure and deformation along the San Andreas Fault in Cholame Valley, California, based on high-resolution reflection profiling

The mapped active traces of the San Andreas fault are separated by a 1-km-wide right-stepping offset in Cholame Valley. The geometry of this offset, defined in other strike-slip systems as a releasing bend or a dilational jog, has resulted in the formation of a pullapart basin. Various researchers have inferred that this offset served as a rupture terminus for earthquakes on both strands of the San Andreas fault (1966 Parkfield and 1857 Fort Tejon); thus, this en echelon offset may represent a barrier to the propagation of rupture between two segments of the San Andreas fault. We collected 18 km of high-resolution seismic reflection data specifically designed to image the San Andreas fault zone in the shallow crust surrounding this offset. Four short profiles (≤ 3.3 km long) of Mini-Sosie reflection data (1 s of two-way travel time ≈ 1.5 km deep) were collected perpendicular to the San Andreas fault; three of these profiles were tied by a 7-km-long profile that trended northwest through Cholame Valley, subparallel to the San Andreas fault. A zone of incoherent energy, narrow at the surface but widening with depth, underlies the mapped active traces of the San Andreas fault and abruptly terminates shallow reflections on both sides. The reflection profiles and available well data indicate that west of the mapped active traces of the San Andreas fault the shallow subsurface structure of the crust consists of thin (≤400 m thick), offset packages of reflections, laterally coherent on the scale of tens of meters, overlying deformed clastic sedimentary rocks. East of the San Andreas fault, the structure of the shallow crust in southern Cholame Valley is characterized by thick packages of reflections, laterally coherent on the scale of kilometers, overlying the Franciscan complex. All of the strata east of the fault (within Cholame Valley) dip toward the San Andreas fault and the offset, into an approximately 1-km-deep sedimentary basin abutting the south strand of San Andreas fault. The offset in Cholame Valley is characterized by a gentle downwarping of sediments into the offset, the presence of many small faults and discontinuous reflections (consistent with areally distributed shear) between the traces of the main fault, localized subsidence abutting the main strike-slip fault, the formation of a basin, near the offset, that is about as deep as the jog is wide, and the southward propagation of the deformation associated with the offset. Strain field modeling based on simple geometries of the San Andreas and associated faults successfully predicts the general features of the observed topography and subsurface structure of southern Cholame Valley, including subsidence and basin formation near the offset, a discontinuous San Andreas fault plane, and at least one fault in southeastern Cholame Valley.

Journal of Geophysical Research↗

USGS2015

The Palos Verdes Fault offshore southern California: late Pleistocene to present tectonic geomorphology, seascape evolution and slip rate estimate based on AUV and ROV surveys

The Palos Verdes Fault (PVF) is one of few active faults in Southern California that crosses the shoreline and can be studied using both terrestrial and subaqueous methodologies. To characterize the near-seafloor fault morphology, tectonic influences on continental slope sedimentary processes and late Pleistocene to present slip rate, a grid of high-resolution multibeam bathymetric data, and chirp subbottom profiles were acquired with an autonomous underwater vehicle (AUV) along the main trace of PVF in water depths between 250 and 600 m. Radiocarbon dates were obtained from vibracores collected using a remotely operated vehicle (ROV) and ship-based gravity cores. The PVF is expressed as a well-defined seafloor lineation marked by subtle along-strike bends. Right-stepping transtensional bends exert first-order control on sediment flow dynamics and the spatial distribution of Holocene depocenters; deformed strata within a small pull-apart basin record punctuated growth faulting associated with at least three Holocene surface ruptures. An upper (shallower) landslide scarp, a buried sedimentary mound, and a deeper scarp have been right-laterally offset across the PVF by 55 ± 5, 52 ± 4 , and 39 ± 8 m, respectively. The ages of the upper scarp and buried mound are approximately 31 ka; the age of the deeper scarp is bracketed to 17–24 ka. These three piercing points bracket the late Pleistocene to present slip rate to 1.3–2.8 mm/yr and provide a best estimate of 1.6–1.9 mm/yr. The deformation observed along the PVF is characteristic of strike-slip faulting and accounts for 20–30% of the total right-lateral slip budget accommodated offshore Southern California.

California↗

USGS2011

Hydrogeologic setting and simulation of groundwater flow near the Canterbury and Leadville Mine Drainage Tunnels, Leadville, Colorado

The Leadville mining district is historically one of the most heavily mined regions in the world producing large quantities of gold, silver, lead, zinc, copper, and manganese since the 1860s. A multidisciplinary investigation was conducted by the U.S. Geological Survey, in cooperation with the Colorado Department of Public Health and Environment, to characterize large-scale groundwater flow in a 13 square-kilometer region encompassing the Canterbury Tunnel and the Leadville Mine Drainage Tunnel near Leadville, Colorado. The primary objective of the investigation was to evaluate whether a substantial hydraulic connection is present between the Canterbury Tunnel and Leadville Mine Drainage Tunnel for current (2008) hydrologic conditions. Altitude in the Leadville area ranges from about 3,018 m (9,900 ft) along the Arkansas River valley to about 4,270 m (14,000 ft) along the Continental Divide east of Leadville, and the high altitude of the area results in a moderate subpolar climate. Winter precipitation as snow was about three times greater than summer precipitation as rain, and in general, both winter and summer precipitation were greatest at higher altitudes. Winter and summer precipitation have increased since 2002 coinciding with the observed water-level rise near the Leadville Mine Drainage Tunnel that began in 2003. The weather patterns and hydrology exhibit strong seasonality with an annual cycle of cold winters with large snowfall, followed by spring snowmelt, runoff, and recharge (high-flow) conditions, and then base-flow (low-flow) conditions in the fall prior to the next winter. Groundwater occurs in the Paleozoic and Precambrian fractured-rock aquifers and in a Quaternary alluvial aquifer along the East Fork Arkansas River, and groundwater levels also exhibit seasonal, although delayed, patterns in response to the annual hydrologic cycle. A three-dimensional digital representation of the extensively faulted bedrock was developed and a geophysical direct-current resistivity field survey was performed to evaluate the geologic structure of the study area. The results show that the Canterbury Tunnel is located in a downthrown structural block that is not in direct physical connection with the Leadville Mine Drainage Tunnel. The presence of this structural discontinuity implies there is no direct groundwater pathway between the tunnels along a laterally continuous bedrock unit. Water-quality results for pH and major-ion concentrations near the Canterbury Tunnel showed that acid mine drainage has not affected groundwater quality. Stable-isotope ratios of hydrogen and oxygen in water indicate that snowmelt is the primary source of groundwater recharge. On the basis of chlorofluorocarbon and tritium concentrations and mixing ratios for groundwater samples, young groundwater (groundwater recharged after 1953) was indicated at well locations upgradient from and in a fault block separate from the Canterbury Tunnel. Samples from sites downgradient from the Canterbury Tunnel were mixtures of young and old (pre-1953) groundwater and likely represent snowmelt recharge mixed with older regional groundwater that discharges from the bedrock units to the Arkansas River valley. Discharge from the Canterbury Tunnel contained the greatest percentage of old (pre-1953) groundwater with a mixture of about 25 percent young water and about 75 percent old water. A calibrated three-dimensional groundwater model representing high-flow conditions was used to evaluate large-scale flow characteristics of the groundwater and to assess whether a substantial hydraulic connection was present between the Canterbury Tunnel and Leadville Mine Drainage Tunnel. As simulated, the faults restrict local flow in many areas, but the fracture-damage zones adjacent to the faults allow groundwater to move along faults. Water-budget results indicate that groundwater flow across the lateral edges of the model controlled the majority of flow in and out of the aquifer (79 percent and 63 percent of the total water budget, respectively). The largest contributions to the water budget were groundwater entering from the upper reaches of the watershed and the hydrologic interaction of the groundwater with the East Fork Arkansas River. Potentiometric surface maps of the simulated model results were generated for depths of 50, 100, and 250 m. The surfaces revealed a positive trend in hydraulic head with land-surface altitude and evidence of increased control on fluid movement by the fault network structure at progressively greater depths in the aquifer. Results of advective particle-tracking simulations indicate that the sets of simulated flow paths for the Canterbury Tunnel and the Leadville Mine Drainage Tunnel were mutually exclusive of one another, which also suggested that no major hydraulic connection was present between the tunnels. Particle-tracking simulations also revealed that although the fault network generally restricted groundwater movement locally, hydrologic conditions were such that groundwater did cross the fault network at many locations. This cross-fault movement indicates that the fault network controls regional groundwater flow to some degree but is not a complete barrier to flow. The cumulative distributions of adjusted age results for the watershed indicate that approximately 30 percent of the flow pathways transmit groundwater that was younger than 68 years old (post-1941) and that about 70 percent of the flow pathways transmit old groundwater. The particle-tracking results are consistent with the apparent ages and mixing ratios developed from the chlorofluorocarbon and tritium results. The model simulations also indicate that approximately 50 percent of the groundwater flowing through the study area was less than 200 years old and about 50 percent of the groundwater flowing through the study area is old water stored in low-permeability geologic units and fault blocks. As a final examination of model response, the conductance parameters of the Canterbury Tunnel and Leadville Mine Drainage Tunnel were manually adjusted from the calibrated values to determine if altering the flow discharge in one tunnel affects the hydraulic behavior in the other tunnel. The examination showed no substantial hydraulic connection. The multidisciplinary investigation yielded an improved understanding of groundwater characteristics near the Canterbury Tunnel and the Leadville Mine Drainage Tunnel. Movement of groundwater between the Canterbury Tunnel and Leadville Mine Drainage Tunnel that was central to this investigation could not be evaluated with strong certainty owing to the structural complexity of the region, study simplifications, and the absence of observation data within the upper sections of the Canterbury Tunnel and between the Canterbury Tunnel and the Leadville Mine Drainage Tunnel. There was, however, collaborative agreement between all of the analyses performed during this investigation that a substantial hydraulic connection did not exist between the Canterbury Tunnel and the Leadville Mine Drainage Tunnel under natural flow conditions near the time of this investigation.

Colorado↗

USGS2012

Hydrologic effects on water level changes associated with episodic fault creep near Parkfield, California

As part of the Parkfield, California, earthquake prediction experiment, water level is monitored in a well 460 m from the main trace of the San Andreas fault on Middle Mountain, in the preparation zone of the anticipated Parkfield earthquake. The well configuration allows water level to be monitored in two fluid reservoirs at depths of 85 and 250 m below land surface. During 1987, water level changes were recorded during 12 of the 18 episodes of accelerated fault creep detected by a creep meter spanning the fault trace 750 m northwest of the well. The creep-related water level changes in the shallow reservoir have durations of less than 1 day, whereas in the deeper reservoir the changes persist for as long as 2 months. These data suggest that the transient nature of the water level changes in the shallow interval is due to vertical flow to the water table and is not evidence that creep events propagate past the well. Phase leads of earth tidal constituents in the water level data from the shallow interval relative to the same constituents in the local volume strain tide support the interpretation of significant flow to the water table at periods of 1 day or less. The form of the water level changes in the deep interval is affected by horizontal flow to the well bore. This effect can be removed from the water level records using a theoretical response curve constrained by the phases of earth tidal constituents in the deep interval relative to the local volume strain tide. For the events where the signal in the shallow interval has been large enough to measure, the sizes of the simultaneous water level changes in the two reservoirs are consistent with the same amounts of volume strain occurring at both depths.

Journal of Geophysical Research Solid Earth↗