Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Interpretation Research”

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

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

At least 289 records · Page 16Linked to original sources
USGS2008

Introduction to the special section northeast Japan: A case history of subduction

Subduction-related tectonic events such as strike-slip faulting, uplift of high-pressure metamorphic rocks in the forearc, back arc spreading, arc magmatism, and continental collisions have shaped northeast Japan, making it a case history for subduction. Much of the record of these events is preserved and can be used to reconstruct the tectonic history of the region. Although many geological, geophysical, and geochemical data have been collected in Japan for the last 30 years, interpretation of these data into a plate tectonic context did not really begin until about 10 years ago. Much of the data and interpretation on northeast Japan have been published in Japanese or as individual local studies. The primary goal of this special section is to present new investigations covering a wide range of topics and scales and place them in a tectonic framework. The papers in this special section focus on three main subjects: Mesozoic aspects of northeast Japan; the relation of Japan Sea opening with the tectonics and geochemistry of northeast Japan in the Tertiary; and the geophysical setting. The purpose of this overview is to link these subjects and provide a tectonic framework for north- east Japan and eastern Asia since the Mesozoic.

Journal of Geophysical Research B: Solid Earth↗
USGS2016

Lithospheric flexure under the Hawaiian volcanic load: Internal stresses and a broken plate revealed by earthquakes

Several lines of earthquake evidence indicate that the lithospheric plate is broken under the load of the island of Hawai`i, where the geometry of the lithosphere is circular with a central depression. The plate bends concave downward surrounding a stress-free hole, rather than bending concave upward as with past assumptions. Earthquake focal mechanisms show that the center of load stress and the weak hole is between the summits of Mauna Loa and Mauna Kea where the load is greatest. The earthquake gap at 21 km depth coincides with the predicted neutral plane of flexure where horizontal stress changes sign. Focal mechanism P axes below the neutral plane display a striking radial pattern pointing to the stress center. Earthquakes above the neutral plane in the north part of the island have opposite stress patterns; T axes tend to be radial. The M 6.2 Honomu and M 6.7 Kiholo main shocks (both at 39 km depth) are below the neutral plane and show radial compression, and the M 6.0 Kiholo aftershock above the neutral plane has tangential compression. Earthquakes deeper than 20 km define a donut of seismicity around the stress center where flexural bending is a maximum. The hole is interpreted as the soft center where the lithospheric plate is broken. Kilauea's deep conduit is seismically active because it is in the ring of maximum bending. A simplified two-dimensional stress model for a bending slab with a load at one end yields stress orientations that agree with earthquake stress axes and radial P axes below the neutral plane. A previous inversion of deep Hawaiian focal mechanisms found a circular solution around the stress center that agrees with the model. For horizontal faults, the shear stress within the bending slab matches the slip in the deep Kilauea seismic zone and enhances outward slip of active flanks.

Hawaii↗
USGS2018

Seismicity induced by wastewater injection in Washington County, Ohio: Influence of preexisting structure, regional stress regime, and well operations

Recent seismicity in Washington County, Ohio, has been suggested to be induced by wastewater disposal operations despite injection ~2 km above the Precambrian basement. We investigated the relationships between disposal well locations and operational histories, spatiotemporal patterns of seismicity enhanced by waveform correlation, and mapped subsurface structures. We also analyzed proxies for in situ stress conditions to evaluate existing principal stress orientations and magnitudes. Double‐difference relocations of small ( M L ≤ 2.1) seismicity revealed linear clusters defining a NE‐SW orientation that parallels the trend of small‐amplitude folds in Upper Paleozoic rocks within 1–2 km of the injection well. Hypocenter depths of ~3.1 to 4.3 km indicate seismicity occurred along faults below the Silurian injection interval in lower Paleozoic strata and Precambrian crystalline basement. Similar folds/faults imaged in nearby seismic reflection lines were associated with basement‐involved fault systems that cut the injection interval targeted by the injection well, providing permeability pathways for fluid pressure increases that initiated slip. Interpreted faults were also optimally oriented in the regional stress field for reactivation due to injection‐related decreases in effective stress. Well completion records suggest that during the recorded seismicity, the Long Run‐1 well was injecting disposal fluids close to or above the fracture pressure of the targeted Silurian injection interval. A comparison of observed seismicity rates and monthly injection records from the disposal well also indicate a close correspondence between earthquake frequency and injection rate. Taken together, these results provide a physical framework for understanding why wastewater disposal operations induced the Washington County seismic sequence.

Ohio↗
USGS2023

Ten best practices for effective phenological research

The number and diversity of phenological studies has increased rapidly in recent years. Innovative experiments, field studies, citizen science projects, and analyses of newly available historical data are contributing insights that advance our understanding of ecological and evolutionary responses to the environment, particularly climate change. However, many phenological data sets have peculiarities that are not immediately obvious and can lead to mistakes in analyses and interpretation of results. This paper aims to help researchers, especially those new to the field of phenology, understand challenges and practices that are crucial for effective studies. For example, researchers may fail to account for sampling biases in phenological data, struggle to choose or design a volunteer data collection strategy that adequately fits their project’s needs, or combine data sets in inappropriate ways. We describe ten best practices for designing studies of plant and animal phenology, evaluating data quality, and analyzing data. Practices include accounting for common biases in data, using effective citizen or community science methods, and employing appropriate data when investigating phenological mismatches. We present these best practices to help researchers entering the field take full advantage of the wealth of available data and approaches to advance our understanding of phenology and its implications for ecology.

International Journal of Biometeorology↗
USGS2002

Prehistoric ruptures of the Gurvan Bulag fault, Gobi Altay, Mongolia

The 1957 Gobi Altay M8.3 earthquake in southern Mongolia was associated with the simultaneous rupture of several faults, including the Gurvan Bulag reverse fault, which is located about 25 km south of the main strike-slip Bogd fault. Our study of paleoseismic excavations across the Gurvan Bulag fault suggests that the penultimate surface rupture occurred after 6.0 ka, most likely between 2.6 and 4.4 ka, and a possible earlier rupture occurred after 7.3 ka. Our interpretation of the stratigraphic relations in one of the exposures suggests that at least five earthquakes have generated surface rupture of the Gurvan Bulag fault since the abandonment of an ancient alluvial fan surface. Luminescence dating of sediment associated with this surface indicates that it formed either 26.6 ± 2.1 ka or 16.1 ± 2.0 ka. These data imply that the recurrence intervals for surface faulting on the Gurvan Bulag and Bogd faults are similar, on the order of several thousands of years, but that the penultimate surface ruptures of the two faults did not occur during the same earthquake.

Gobi Altay↗
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↗
USGS2012

Limitations of quantitative analysis of deep crustal seismic reflection data: Examples from GLIMPCE

Amplitude preservation in seismic reflection data can be obtained by a relative true amplitude (RTA) processing technique in which the relative strength of reflection amplitudes is preserved vertically as well as horizontally, after compensating for amplitude distortion by near-surface effects and propagation effects. Quantitative analysis of relative true amplitudes of the Great Lakes International Multidisciplinary Program on Crustal Evolution seismic data is hampered by large uncertainties in estimates of the water bottom reflection coefficient and the vertical amplitude correction and by inadequate noise suppression. Processing techniques such as deconvolution, F-K filtering, and migration significantly change the overall shape of amplitude curves and hence calculation of reflection coefficients and average reflectance. Thus lithological interpretation of deep crustal seismic data based on the absolute value of estimated reflection strength alone is meaningless. The relative strength of individual events, however, is preserved on curves generated at different stages in the processing. We suggest that qualitative comparisons of relative strength, if used carefully, provide a meaningful measure of variations in reflectivity. Simple theoretical models indicate that peg-leg multiples rather than water bottom multiples are the most severe source of noise contamination. These multiples are extremely difficult to remove when the water bottom reflection coefficient is large (>0.6), a condition that exists beneath parts of Lake Superior and most of Lake Huron.

Great Lakes↗
USGS2012

Reservoir analysis of the Denver earthquakes: A case of induced seismicity

Injection of fluid wastes into the fractured Precambrian crystalline bedrock beneath the Rocky Mountain Arsenal near Denver triggered earthquakes in the 1960's. An analysis, based on the assumption that fluid flow in the fractured reservoir can be approximated by flow in a porous medium, is presented. The configuration and hydrologic properties of the reservoir are determined from two lines of evidence: (1) locations of earthquake hypocenters determined by seismic arrays installed at the Arsenal and (2) observed long-term decline in fluid levels in the injection well. Together these two sets of data indicate that a long, narrow reservoir, aligned in the direction N 60°W, exists. The reservoir is 3.35 km in width, extends 30.5 km to the northwest and infinitely to the southeast, and spans a depth interval from 3.7 to 7.0 km below land surface. It has a transmissivity of 1.08×10-5 m2/s and a storage coefficient of 1.0×10-5. Computed pressure buildup along the length of the reservoir is compared with the spatial distribution of earthquake epicenters. The comparison shows that earthquakes are confined to that part of the reservoir where the pressure buildup exceeds 32 bars. This critical value is interpreted as the pressure buildup above which earthquakes occur. The migration of earthquake epicenters away from the injection well, a phenomenon noted by previous investigators, can be accounted for by the outward propagation of the critical pressure buildup. The analysis is extended to examining the effects of rapid flow in fractures opened by high injection pressure. The results show that the effect is confined to a small region within 1 km of the injection well. The existence of a critical pressure buildup above which earthquakes occur is completely consistent with the theory on the role of fluid pressure in fault movement as presented by Hubbert and Rubey.

Colorado↗
USGS1975

Preliminary results of a gravity survey of the Henrys Lake quadrangle, Idaho and Montana

A gravity survey of the Henrys Lake quadrangle shows that a gravity low with about 10 milligals of closure coincides with the Henrys Lake basin. The low is interpreted to reflect a basin fill of 1,100 m or more of Cenozoic sediments and volcanic rock. The data indicate that on the east and probably on the west the basin is bounded by northwest- and north-trending faults, and near its waist by the east-trending Centennial fault. A model constructed across the Henrys Lake basin suggests a northwest-trending fault, located near the southeast corner of Henrys Lake and concealed beneath the basin fill. The sparse gravity data east of the Centennial Mountains are not sufficient to provide any evidence on the extension of the Centennial fault across the Henrys Lake basin.

Idaho, Montana, Wyoming↗
USGS2001

Upper crustal structure in Puget Lowland, Washington: Results from the 1998 Seismic Hazards Investigation in Puget Sound

A new three-dimensional (3-D) model shows seismic velocities beneath the Puget Lowland to a depth of 11 km. The model is based on a tomographic inversion of nearly one million first-arrival travel times recorded during the 1998 Seismic Hazards Investigation in Puget Sound (SHIPS), allowing higher-resolution mapping of subsurface structures than previously possible. The model allows us to refine the subsurface geometry of previously proposed faults (e.g., Seattle, Hood Canal, southern Whidbey Island, and Devils Mountain fault zones) as well as to identify structures (Tacoma, Lofall, and Sequim fault zones) that warrant additional study. The largest and most important of these newly identified structures lies along the northern boundary of the Tacoma basin; we informally refer to this structure here as the Tacoma fault zone. Although tomography cannot provide information on the recency of motion on any structure, Holocene earthquake activity on the Tacoma fault zone is suggested by seismicity along it and paleoseismic evidence for abrupt uplift of tidal marsh deposits to its north. The tomography reveals four large, west to northwest trending low-velocity basins (Tacoma, Seattle, Everett, and Port Townsend) separated by regions of higher velocity ridges that are coincident with fault-bounded uplifts of Eocene Crescent Formation basalt and pre-Tertiary basement. The shapes of the basins and uplifts are similar to those observed in gravity data; gravity anomalies calculated from the 3-D tomography model are in close agreement with the observed anomalies. In velocity cross sections the Tacoma and Seattle basins are asymmetric: the basin floor dips gently toward a steep boundary with the adjacent high-velocity uplift, locally with a velocity "overhang" that suggests a basin vergent thrust fault boundary. Crustal fault zones grow from minor folds into much larger structures along strike. Inferred structural relief across the Tacoma fault zone increases by several kilometers westward along the fault zone to Lynch Cove, where we interpret it as a zone of south vergent faulting overthrusting Tacoma basin. In contrast, structural relief along the Seattle fault zone decreases west of Seattle, which we interpret as evidence that the N-S directed compression is being accommodated by slip transfer between the Seattle and Tacoma fault zones. Together, the Tacoma and Seattle fault zones raise the Seattle uplift, one of a series of east-west trending, pop-up structures underlying Puget Lowland from the Black Hills to the San Juan Islands.

Washington↗
USGS2019

Heat and mass transport in a vapor-dominated hydrothermal area in Yellowstone National Park, USA: Inferences from magnetic, electrical, electromagnetic, subsurface temperature and diffuse CO2 flux measurements

Vapor‐dominated hydrothermal systems are characterized by localized and elevated heat and gas flux. In these systems, steam and gas ascend from a boiling water reservoir, steam condenses beneath a low‐permeability cap layer, and liquid water descends, driven by gravity (“heat pipe” model). We combine magnetic, electromagnetic, and geoelectrical methods and CO 2 flux and subsurface temperature measurements in the Solfatara Plateau Thermal Area in the Yellowstone Caldera to address several fundamental questions: (1) What are the structural and/or lithological controls on heat and mass transport in vapor‐dominated areas? (2) What is the geometry and size of convecting multiphase thermal plumes? (3) Are thermal plumes associated with subsurface rock alteration and demagnetization? Magnetic and electromagnetic data inversions suggest an asymmetric 50‐ to 100‐m thick basin of glacial deposits with the thickest part adjacent to the margin of a rhyolite flow. The 3‐D electrical conductivity model in the glacial basin reveals a narrow vertical conductor interpreted as a focused multiphase plume, which coincides at the ground surface with the heat and CO 2 flux maxima. The magnetic data suggest that destruction of magnetic minerals due to rock alteration associated with the hydrothermal plume occurs mainly near the ground surface. We propose a model where the buoyant multiphase plume forms in response to decompression, boiling, and phase separation of pressurized thermal groundwater that discharges from the brecciated base of a rhyolite flow into the basin of glacial deposits. Results from multiphase groundwater flow and heat transport numerical simulations corroborate the first‐order characteristics of this model.

Wyoming↗
USGS2014

Fine scale habitat use by age-1 stocked muskellunge and wild northern pike in an upper St. Lawrence River bay

Radio telemetry of stocked muskellunge ( n = 6) and wild northern pike ( n = 6) was used to track late summer and fall movements from a common release point in a known shared nursery bay to test the hypothesis that age-1 northern pike and stocked muskellunge segregate and have different habitat affinities. Water depth, temperature, substrate and aquatic vegetation variables were estimated for each muskellunge ( n = 103) and northern pike ( n = 131) position and nested ANOVA comparisons by species indicated differences in habitat use. Muskellunge exhibited a greater displacement from the release point and used habitat in shallower water depths (mean = 0.85 m, SE = 0.10) than northern pike (mean = 1.45 m, SE = 0.08). Both principal components analysis (PCA) and principal components ordination (PCO) were used to interpret underlying gradients relative to fish positions in two-dimensional space. Our analysis indicated that a separation of age-1 northern pike and muskellunge occurred 7 d post-release. This first principal component explained 48% of the variation in habitat use. Northern pike locations were associated with deeper habitats that generally had softer silt substrates and dense submersed vegetation. Muskellunge locations post-acclimation showed greater association with shallower habitats containing firmer sandy and clay substrates and emergent vegetation. The observed differences in habitat use suggest that fine-scale ecological separation occurred between these stocked muskellunge and wild northern pike, but small sample sizes and potential for individual variation limit extension of these conclusions. Further research is needed to determine if these patterns exist between larger samples of fishes over a greater range of habitats.

Journal of Great Lakes Research↗
USGS2024

Capturing potential: Leveraging grass carp behavior Ctenopharyngodon idella for enhanced removal

Effective management of invasive species benefits from a comprehensive understanding of the species’ behavior and interactions with the invaded system. We investigated temporal dynamics of telemetry detections and the potential utility of a traitor approach for informing response efforts to the invasive grass carp ( Ctenopharyngodon idella ) population in the Sandusky River, a major tributary to Lake Erie. Telemetered grass carp exhibited heightened activity at night and early morning, suggesting that capture and removal be more effective during these time periods. Analysis of catch per unit effort (CPUE) across different removal methods, trammel nets, electrofishing, and hoop nets. suggested that incorporating the traitor approach could improve capture. Low catchability values (<0.026), based on the number of telemetered grass carp present in the river on a weekly basis and the number of those telemetered fish captured, suggest the species is difficult to capture. Optimizing response effort efficiency is important and refining catchability estimates will lessen errors in population models and improve interpretation of low CPUE data. Results from generalized additive models suggest capture could be improved using telemetry data, night removals, and by attempting exploratory removal efforts in fall and winter months. By incorporating telemetry data and acknowledging the complexities of grass carp behavior and ecology, we found that a multifaceted and data-driven approach to invasive species control could be beneficial, ultimately promoting conservation and sustainability in dynamic ecosystems like Lake Erie.

Ohio↗
USGS2001

Asian dust events of April 1998

On April 15 and 19, 1998, two intense dust storms were generated over the Gobi desert by springtime low-pressure systems descending from the northwest. The windblown dust was detected and its evolution followed by its yellow color on SeaWiFS satellite images, routine surface-based monitoring, and through serendipitous observations. The April 15 dust cloud was recirculating, and it was removed by a precipitating weather system over east Asia. The April 19 dust cloud crossed the Pacific Ocean in 5 days, subsided to the surface along the mountain ranges between British Columbia and California, and impacted severely the optical and the concentration environments of the region. In east Asia the dust clouds increased the albedo over the cloudless ocean and land by up to 10–20%, but it reduced the near-UV cloud reflectance, causing a yellow coloration of all surfaces. The yellow colored backscattering by the dust eludes a plausible explanation using simple Mie theory with constant refractive index. Over the West Coast the dust layer has increased the spectrally uniform optical depth to about 0.4, reduced the direct solar radiation by 30–40%, doubled the diffuse radiation, and caused a whitish discoloration of the blue sky. On April 29 the average excess surface-level dust aerosol concentration over the valleys of the West Coast was about 20–50 μg/m 3 with local peaks >100 μg/m 3 . The dust mass mean diameter was 2–3 μm, and the dust chemical fingerprints were evident throughout the West Coast and extended to Minnesota. The April 1998 dust event has impacted the surface aerosol concentration 2–4 times more than any other dust event since 1988. The dust events were observed and interpreted by an ad hoc international web-based virtual community. It would be useful to set up a community-supported web-based infrastructure to monitor the global aerosol pattern for such extreme aerosol events, to alert and to inform the interested communities, and to facilitate collaborative analysis for improved air quality and disaster management.

Journal of Geophysical Research D: Atmospheres↗
USGS1996

Movement patterns of wintering lesser scaup in Grand Calumet River - Indiana Harbor Canal, Indiana

Lesser scaup ( Aythya affinis ) winter in the Grand Calumet River—Indiana Harbor Canal (GCR-IHC) drainage into southern Lake Michigan, a polluted oasis in the highly urban and industrial corridor of northern Illinois and Indiana. The GCR-IHC is an important wintering area for lesser scaup in northwestern Indiana, especially after Lake Michigan freezes, because of the lack of other wildlife habitat in the area. We investigated the feeding and movement patterns of these wintering lesser scaup to better understand the ecological significance of this area to wintering waterfowl and to interpret contaminant levels in these lesser scaup. We implanted radio transmitters in 20 individuals in January 1994 and tracked them until early March 1994. Four lesser scaup disappeared and four died within the first 2 weeks after implantation; 12 lesser scaup remained in the GCR-IHC for the next 6–7 weeks. Individual radio-marked scaup (n = 12) were located in the GCR-IHC on an average of 62.2% (extreme values = 25–87%) of the searches. No radio-marked lesser scaup were located outside the GCR-IHC. Lesser scaup were feeding during 23% of the radio fixes.

Journal of Great Lakes Research↗
USGS2012

Compositional variations on the Moon: Recalibration of Galileo solid‐state imaging data for the Orientale region and farside

Updated radiometric calibration and systematic processing procedures for Galileo solid‐state imaging (SSI) data from the first (1990) Earth‐Moon encounter are presented. These procedures were applied to a whole‐disk imaging sequence of the Moon centered near Mare Orientale, called Lunmap 14 (L14). Processing of L14 data included radiometric calibration, subpixel coregistration, scattered light removal, geometric control and reprojection, photometric normalization, and calibration to Earth‐based spectra. Coregistration and scattered‐light removal procedures are improvements over the initial calibration of the SSI mosaics. The effects of scattered light correction are best seen using a whole‐Moon view such as L14; resolution of the debated amounts of light scattering from within or outside the camera field of view is not necessary. Scattered light removal particularly affects the 1‐μm spectral region and has implications for interpretation of mafic mineral signatures in mare deposits of the lunar limb and farside. Recalibrated spectra indicate that mare ponds of the limb show moderately deep 1‐μm absorptions, and thus mafic mineral contents, comparable to those of other nearside basalts. Mafic mineral contents of Schiller‐Shickard cryptomaria are higher than previously thought and are similar to some low‐Ti nearside basalts. Many of the recalibrated spectra from South Pole/Aitken are similar to those of Schiller‐Schickard cryptomaria, suggesting that many of these soils represent a mixed mare/highland lithology. The hypothesis that there is an olivine enrichment in southern South Pole/Aitken basin is not supported strongly by spectra shown here.

Journal of Geophysical Research: Planets↗
USGS2007

Three-dimensional P wave velocity model for the San Francisco Bay region, California

[1] A new three-dimensional P wave velocity model for the greater San Francisco Bay region has been derived using the double-difference seismic tomography method, using data from about 5,500 chemical explosions or air gun blasts and approximately 6,000 earthquakes. The model region covers 140 km NE-SW by 240 km NW-SE, extending from 20 km south of Monterey to Santa Rosa and reaching from the Pacific coast to the edge of the Great Valley. Our model provides the first regional view of a number of basement highs that are imaged in the uppermost few kilometers of the model, and images a number of velocity anomaly lows associated with known Mesozoic and Cenozoic basins in the study area. High velocity ( V p > 6.5 km/s) features at ∼15-km depth beneath part of the edge of the Great Valley and along the San Francisco peninsula are interpreted as ophiolite bodies. The relocated earthquakes provide a clear picture of the geometry of the major faults in the region, illuminating fault dips that are generally consistent with previous studies. Ninety-five percent of the earthquakes have depths between 2.3 and 15.2 km, and the corresponding seismic velocities at the hypocenters range from 4.8 km/s (presumably corresponding to Franciscan basement or Mesozoic sedimentary rocks of the Great Valley Sequence) to 6.8 km/s. The top of the seismogenic zone is thus largely controlled by basement depth, but the base of the seismogenic zone is not restricted to seismic velocities of ≤6.3 km/s in this region, as had been previously proposed.

California↗
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↗