Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Interpretation”

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

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

At least 1,117 records · Page 62Linked to original sources

Structure and Paleozoic stratigraphy of a complex of thrust plates in the Fish Creek Reservoir area, south-central Idaho

Permian, Pennsylvanian, Mississippian, Devonian, and Silurian marine rocks of diverse facies are brought together in a complex of six thrust sheets in the Fish Creek Reservoir area on the north edge of the Snake River Plain, Idaho. The lowest structural element, the parautochthon, is made of more than 450 m (1,500 ft) of folded and faulted Devonian miogeosynclinal carbonate rocks present in a 6.5-km 2 (2.5-mi 2 ) window. Along the east margin of the window, a sliver of continental margin transitional carbonate rocks of Early Devonian and Late Silurian age assigned to the Roberts Mountains Formation is thrust over the miogeosynclinal rocks. The window of middle Paleozoic rocks is overridden along the Fish Creek thrust fault by the flysch facies of the Copper Basin Formation, a turbidite-submarine-fan sequence more than 1,000 m (3,300 ft) thick, of Mississippian age. About 4.8 km (3 mi) southwest of the window, about 100 m (300 ft) of deepwater siliceous oceanic facies clastic rocks are exposed, which are assigned with question to the Milligen(?) Formation of Devonian age. These clastic rocks are interpreted to be thrust over the Copper Basin Formation. The highest structural elements are sequences more than 610 m (2,000 ft) thick of interbedded sandy and conglomeratic limestones, quartzites, and conglomerates and interbedded siltstones and argillites of the Wood River Formation of Middle Pennsylvanian to Early Permian age. The Wood River Formation is in thrust contact with the Milligen(?) Formation in the southwest part of the mapped area and with Copper Basin Formation along the west side of Fish Creek Reservoir. All the thrust sheets have moved eastward. The minimum distance moved is estimated from sedimentation models and facies reconstructions to range from perhaps several kilometres for the allochthon of the Roberts Mountains Formation to 48 km (30 mi) for the Milligen(?) Formation allochthon. The principal period of thrusting was post-Early Permian (post-Wood River Formation) and preEocene (pre-Challis Volcanics) and is of probable Sevier age. Middle Paleozoic rocks of the Milligen and Roberts Mountains Formations, however, also may have been involved in an earlier period of thrusting of latest Devonian to earliest Mississippian age related to the Antler orogeny. The thrust sheets were deformed into a northwest-trending dome in late Mesozoic time and were broken by basin-range faults during the Tertiary.

Idaho↗

Recalculated probability of M ≥ 7 earthquakes beneath the Sea of Marmara, Turkey

New earthquake probability calculations are made for the Sea of Marmara region and the city of Istanbul, providing a revised forecast and an evaluation of time-dependent interaction techniques. Calculations incorporate newly obtained bathymetric images of the North Anatolian fault beneath the Sea of Marmara [Le Pichon et al., 2001; Armijo et al., 2002]. Newly interpreted fault segmentation enables an improved regional A.D. 1500-2000 earthquake catalog and interevent model, which form the basis for time-dependent probability estimates. Calculations presented here also employ detailed models of coseismic and postseismic slip associated with the 17 August 1999 M = 7.4 Izmit earthquake to investigate effects of stress transfer on seismic hazard. Probability changes caused by the 1999 shock depend on Marmara Sea fault-stressing rates, which are calculated with a new finite element model. The combined 2004-2034 regional Poisson probability of M≥7 earthquakes is ~38%, the regional time-dependent probability is 44 ± 18%, and incorporation of stress transfer raises it to 53 ± 18%. The most important effect of adding time dependence and stress transfer to the calculations is an increase in the 30 year probability of a M ??? 7 earthquake affecting Istanbul. The 30 year Poisson probability at Istanbul is 21%, and the addition of time dependence and stress transfer raises it to 41 ± 14%. The ranges given on probability values are sensitivities of the calculations to input parameters determined by Monte Carlo analysis; 1000 calculations are made using parameters drawn at random from distributions. Sensitivities are large relative to mean probability values and enhancements caused by stress transfer, reflecting a poor understanding of large-earthquake aperiodicity.

Sea Of Marmara↗

Modeling haul-out behavior of walruses in Bering Sea ice

Understanding haul-out behavior of ice-associated pinnipeds is essential for designing and interpreting popula-tion surveys and for assessing effects of potential changes in their ice environments. We used satellite-linked transmitters to obtain sequential information about location and haul-out state for Pacific walruses, Odobenus rosmarus divergens (Il-liger, 1815), in the Bering Sea during April of 2004, 2005, and 2006. We used these data in a generalized mixed model of haul-out bout durations and a hierarchical Bayesian model of haul-out probabilities to assess factors related to walrus haul-out behavior, and provide the first predictive model of walrus haul-out behavior in sea ice habitat. Average haul-out bout duration was 9 h, but durations of haul-out bouts tended to increase with durations of preceding in-water bouts. On aver-age, tagged walruses spent only about 17% of their time hauled out on sea ice. Probability of being hauled out decreased with wind speed, increased with temperature, and followed a diurnal cycle with the highest values in the evening. Our haul-out probability model can be used to estimate the proportion of the population that is unavailable for detection in spring surveys of Pacific walruses on sea ice.

Canadian Journal of Zoology↗

New method for quantification of vuggy porosity from digital optical borehole images as applied to the karstic Pleistocene limestone of the Biscayne aquifer, southeastern Florida

Vuggy porosity is gas- or fluid-filled openings in rock matrix that are large enough to be seen with the unaided eye. Well-connected vugs can form major conduits for flow of ground water, especially in carbonate rocks. This paper presents a new method for quantification of vuggy porosity calculated from digital borehole images collected from 47 test coreholes that penetrate the karstic Pleistocene limestone of the Biscayne aquifer, southeastern Florida. Basically, the method interprets vugs and background based on the grayscale color of each in digital borehole images and calculates a percentage of vuggy porosity. Development of the method was complicated because environmental conditions created an uneven grayscale contrast in the borehole images that makes it difficult to distinguish vugs from background. The irregular contrast was produced by unbalanced illumination of the borehole wall, which was a result of eccentering of the borehole-image logging tool. Experimentation showed that a simple, single grayscale threshold would not realistically differentiate between the grayscale contrast of vugs and background. Therefore, an equation was developed for an effective subtraction of the changing grayscale contrast, due to uneven illumination, to produce a grayscale threshold that successfully identifies vugs. In the equation, a moving average calculated around the circumference of the borehole and expressed as the background grayscale intensity is defined as a baseline from which to identify a grayscale threshold for vugs. A constant was derived empirically by calibration with vuggy porosity values derived from digital images of slabbed-core samples and used to make the subtraction from the background baseline to derive the vug grayscale threshold as a function of azimuth. The method should be effective in estimating vuggy porosity in any carbonate aquifer. ?? 2003 Published by Elsevier B.V.

Journal of Applied Geophysics↗

Discrimination of carbon and nitrogen isotopes from milk to serum and vibrissae in Alaska Steller sea lions (Eumetopias jubatus)

Knowledge of diet-tissue stable isotope discrimination is required to properly interpret stable isotope values and to identify possible diet shifts, such as might be expected from nursing through weaning. This study compared ??13C and ??15N of paired serum and vibrissal roots with those of ingested milk (n = 52) from free-ranging Steller sea lion (Eumetopias jubatus (Schreber, 1776)) pups (1-11 months) and juveniles (14-27 months) to estimate diet-tissue discrimination. Mean 15N enrichment from ingested milk to serum was 2.1??? ?? 0.6%??? and ??15N at the root of the vibrissae (representing current growth) were not significantly different from serum values. Milk was enriched for mean 13C by 5.0??? ?? 1.0%??? and 7.3??? ?? 1.2??? relative to serum and vibrissal roots, respectively, which was due to the presence of 13C-depleted lipids in milk. This was confirmed by lipid extraction from a subset of milk and serum samples, resulting in a 5.8??? ?? 1.0??? change only in milk. This study established that vibrissal roots and serum are reflective of a milk diet with approximately 2.0??? 15N enrichment, and vibrissal roots reflect serum and lipid-extracted milk values with approximately 2.0??? 13C enrichment. These discrimination factors are important to establish for stable isotope studies assessing diet shifts. ?? 2008 NRC.

Canadian Journal of Zoology↗

Interpretation of hydraulic conductivity in a fractured-rock aquifer over increasingly larger length dimensions

A comparison of the hydraulic conductivity over increasingly larger volumes of crystalline rock was conducted in the Piedmont physiographic region near Bethesda, Maryland, USA. Fluid-injection tests were conducted on intervals of boreholes isolating closely spaced fractures. Single-hole tests were conducted by pumping in open boreholes for approximately 30 min, and an interference test was conducted by pumping a single borehole over 3 days while monitoring nearby boreholes. An estimate of the hydraulic conductivity of the rock over hundreds of meters was inferred from simulating groundwater inflow into a kilometer-long section of a Washington Metropolitan Area Transit Authority tunnel in the study area, and a groundwater modeling investigation over the Rock Creek watershed provided an estimate of the hydraulic conductivity over kilometers. The majority of groundwater flow is confined to relatively few fractures at a given location. Boreholes installed to depths of approximately 50 m have one or two highly transmissive fractures; the transmissivity of the remaining fractures ranges over five orders of magnitude. Estimates of hydraulic conductivity over increasingly larger rock volumes varied by less than half an order of magnitude. While many investigations point to increasing hydraulic conductivity as a function of the measurement scale, a comparison with selected investigations shows that the effective hydraulic conductivity estimated over larger volumes of rock can either increase, decrease, or remain stable as a function of the measurement scale. Caution needs to be exhibited in characterizing effective hydraulic properties in fractured rock for the purposes of groundwater management.

Hydrogeology Journal↗

Time series of high-resolution images enhances efforts to monitor post-fire condition and recovery, Waldo Canyon fire, Colorado, USA

Interpretations of post-fire condition and rates of vegetation recovery can influence management priorities, actions and perception of latent risks from landslides and floods. In this study, we used the Waldo Canyon fire (2012, Colorado Springs, Colorado, USA) as a case study to explore how a time series (2011–2016) of high-resolution images can be used to delineate burn extent and severity, as well as quantify post-fire vegetation recovery. We applied an object-based approach to map burn severity and vegetation recovery using Worldview-2, Worldview-3 and QuickBird-2 imagery. The burned area was classified as 51% high, 20% moderate and 29% low burn-severity. Across the burn extent, the shrub cover class showed a rapid recovery, resprouting vigorously within 1 year, whereas 4 years post-fire, areas previously dominated by conifers were divided approximately equally between being classified as dominated by quaking aspen saplings with herbaceous species in the understorey or minimally recovered. Relative to using a pixel-based Normalised Difference Vegetation Index (NDVI), our object-based approach showed higher rates of revegetation. High-resolution imagery can provide an effective means to monitor post-fire site conditions and complement more prevalent efforts with moderate- and coarse-resolution sensors.

Colorado↗

Semiempirical model of soil water hysteresis

In order to represent hysteretic soil water retention curves accurately using as few measurements as possible, a new semiempirical model has been developed. It has two postulates related to physical characteristics of the medium, and two parameters, each with a definite physical interpretation, whose values are determined empirically for a given porous medium. One parameter represents the fraction of the pore space that is not subject to Haines jump hysteresis. Its associated postulate is that a single value of this parameter characterizes a given medium, or, equivalently, that the medium is self-similar with respect to the division between hysteretic and nonhysteretic pore space. The second parameter is the effective body-to-neck size ratio of the medium's largest pore. The second postulate specifies a particular relation between the size distributions of pore bodies and of pore necks. Tests of the model show that it provides high-quality optimized fits to measured water content vs. matric pressure wetting curves for a wide variety of media. The parameter values obtained through optimization correlate plausibly with such media properties as uniformity of particles, complexity of structure, and degree of compaction. A practical use of this model is to provide a complete simulated main wetting curve for a medium where only a main drying curve and two points on the wetting curve have been measured. With additional development, it may be possible to do without the two measured wetting points if parameters can be evaluated from other measurements or known properties of the medium.

Soil Science Society of America Journal↗

Combined interpretation of radar, hydraulic, and tracer data from a fractured-rock aquifer near Mirror Lake, New Hampshire, USA

An integrated interpretation of field experimental cross-hole radar, tracer, and hydraulic data demonstrates the value of combining time-lapse geophysical monitoring with conventional hydrologic measurements for improved characterization of a fractured-rock aquifer. Time-lapse difference-attenuation radar tomography was conducted during saline tracer experiments at the US Geological Survey Fractured Rock Hydrology Research Site near Mirror Lake, Grafton County, New Hampshire, USA. The presence of electrically conductive saline tracer effectively illuminates permeable fractures or pathways for geophysical imaging. The geophysical results guide the construction of three-dimensional numerical models of ground-water flow and solute transport. In an effort to explore alternative explanations for the tracer and tomographic data, a suite of conceptual models involving heterogeneous hydraulic conductivity fields and rate-limited mass transfer are considered. Calibration data include tracer concentrations, the arrival time of peak concentration at the outlet, and steady-state hydraulic head. Results from the coupled inversion procedure suggest that much of the tracer mass migrated outside the three tomographic image planes, and that solute is likely transported by two pathways through the system. This work provides basic and site-specific insights into the control of permeability heterogeneity on ground-water flow and solute transport in fractured rock.

New Hampshire↗

Implications of tagging effects for interpreting the performance of sea lamprey traps in a large river

Abundance estimates can be crucial for managing species of economic concern. The accuracy of these estimates can depend on the methods used to track animals and to estimate abundance from tracking data. We tested experimentally if disparate estimates of trapping efficiency calculated for sea lamprey ( Petromyzon marinus ) in the St. Marys River near Sault Ste. Marie, Canada could be explained by effects related to the invasiveness and handling involved in tagging or the tag size used in the marking procedures. Trapping is used to gauge adult abundance, trapping efficiency, and success of a binational sea lamprey control program in the Laurentian Great Lakes, North America. Our experiment compared nightly catches of sea lamprey marked with external fin clips, surgically-implanted passive integrated transponder tags (PIT-only), and surgically-implanted PIT and acoustic tags (PIT+acoustic). We found no evidence that the probability of being trapped was affected by the added invasiveness and handling of internal tagging. Nightly recaptures of PIT-only tagged sea lamprey, relative to fin-clipped sea lamprey, were not different from expectations based on the numbers of individuals released from each treatment group. Conversely, there was evidence of effects related to tag size. Nightly recaptures of PIT+acoustic tagged sea lamprey, relative to PIT-only tagged sea lamprey, were lower than expected based on numbers of individuals released from each treatment group. Effects related to tag size partially explain the disparate estimates in trapping efficiency observed for sea lamprey.

Ontario↗

Topographic stress perturbations in southern Davis Mountains, west Texas 1. Polarity reversal of principal stresses

We have applied a previously developed analytical stress model to interpret subsurface stress conditions inferred from acoustic televiewer logs obtained in two municipal water wells located in a valley in the southern Davis Mountains near Alpine, Texas. The appearance of stress-induced breakouts with orientations that shift by 90° at two different depths in one of the wells is explained by results from exact solutions for the effects of valleys on gravity and tectonically induced subsurface stresses. The theoretical results demonstrate that above a reference depth termed the hinge point, a location that is dependent on Poisson's ratio, valley shape, and magnitude of the maximum horizontal tectonic stress normal to the long axis of the valley, horizontal stresses parallel to the valley axis are greater than those normal to it. At depths below this hinge point the situation reverses and horizontal stresses normal to the valley axis are greater than those parallel to it. Application of the theoretical model at Alpine is accommodated by the fact that nearby earthquake focal mechanisms establish an extensional stress regime with the regional maximum horizontal principal stress aligned perpendicular to the valley axis. We conclude that the localized stress field associated with a valley setting can be highly variable and that breakouts need to be examined in this context when estimating the orientations and magnitudes of regional principal stresses.

Texas↗

High-rate very-long-period seismicity at Yasur volcano, Vanuatu: source mechanism and decoupling from surficial explosions and infrasound

Yasur volcano, Vanuatu is a continuously active open-vent basaltic-andesite stratocone with persistent and long-lived eruptive activity. We present results from a seismo-acoustic field experiment at Yasur, providing locally dense broad-band seismic and infrasonic network coverage from 2016 July 27 to August 3. We corroborate our seismo-acoustic observations with coincident video data from cameras deployed at the crater and on an unoccupied aircraft system (UAS). The waveforms contain a profusion of signals reflecting Yasur’s rapidly occurring and persistent explosive activity. The typical infrasonic signature of Yasur explosions is a classic short-duration and often asymmetric explosion waveform characterized by a sharp compressive onset and wideband frequency content. The dominant seismic signals are numerous repetitive very-long-period (VLP) signals with periods of ∼2–10 s. The VLP seismic events are ‘high-rate’, reoccurring near-continuously throughout the data set with short interevent times (∼20–60 s). We observe variability in the synchronization of seismic VLP and acoustic sources. Explosion events clearly delineated by infrasonic waveforms are underlain by seismic VLPs. However, strong seismic VLPs also occur with only a weak infrasonic expression. Multiplet analysis of the seismic VLPs reveals a systematic progression in the seismo-acoustic source decoupling. The same dominant seismic VLP multiplet occurs with and without surficial explosions and infrasound, and these transitions occur over a timescale of a few days during our field campaign. We subsequently employ template matching, stacking, and full-waveform inversion to image the source mechanism of the dominant VLP multiplet. Inversion of the dominant VLP multiplet stack points to a composite source consisting of either a dual-crack (plus forces) or pipe-crack (plus forces) mechanism. The derived mechanisms correspond to a point-source directly beneath the summit vents with centroid depths in the range ∼900–1000 m below topography. All mechanisms suggest a northeast trending crack dipping relatively shallowly to the northwest and indicate a VLP source centroid and mechanism controlled by a stable structural geologic feature beneath Yasur. We interpret the results in the framework of gas slug ascent through the conduit responsible for Yasur explosions. The VLP mechanism and timing with infrasound (when present) are explained by a shallow-buffered top-down model in which slug ascent is relatively aseismic until reaching the base of a shallow section. Slug disruption in this shallow zone triggers a pressure disturbance that propagates downward and couples at the conduit base (VLP centroid). If the shallow section is open, an explosion propagates to the surface, producing infrasound. In the case of (the same multiplet) VLPs occurring without surficial explosions and weak or no infrasound, the decoupling of the dominant VLPs at ∼900–1000 m depth from surficial explosions and infrasound strongly indicates buffering of the terminal slug ascent. This buffering could be achieved by a variety of conditions at or directly beneath the vents, such as a high-viscosity layer of crystal-rich magma, a debris cap from backfill, a foam layer, or a combination of these. The dominant VLP at Yasur captured by our experiment has a source depth and mechanism separated from surface processes and is stable over time.

Geophysical Journal International↗

Mapping saline groundwater beneath the Sea Galilee and its vicinity using time domain electromagnetic (TDEM) geophysical technique

An extensive time domain electromagnetic (TDEM) survey covering the Sea of Galilee with a dense grid of points has been recently carried out. A total of 269 offshore and 33 supplementary onshore TDEM soundings were performed along six N-S and ten W-E profiles and at selected points both offshore and onshore along the whole coastal line. The interpreted resistivities were calibrated with the direct salinity measurements in the Haon-2 borehole and relatively deep (5 m) cores taken from the lake bottom. It was found that resistivities below 1 ohm-m are solely indicative of groundwater salinity exceeding 10,000 mg Cl/l. Such low resistivities (high salinities) were detected at depths greater than 15 m below almost the entire bottom of the lake. At some parts of the lake, particularly in the south, the saline water was detected at shallower depths, sometimes at a few meters below the bottom. Relatively high resistivity (fresh groundwater) was found along the margins of the lake down to roughly 100 m, the maximum exploration depth of the system. The detected sharp lateral contrasts at the lake margin between high and low resistivities coincide with the faults separating the carbonate and clastic units, respectively. The geometry of the fresh/saline groundwater interface below the central part of the lake is very similar to the shape of the lake bottom, probably due to the diffusive salt transport from the bottom sediments to the lake water. The above geophysical observations suggest differentsalt transport mechanisms from the sediments to the central part of the lake (diffusion) and from regional aquifers to the margins of the lake (advection). ?? 2004 Science From Israel/LPPLtd.

Israel Journal of Earth Sciences↗

Since “Groundwater and surface water–A single resource”: some U.S. Geological Survey advances in modeling groundwater/surface-water interactions

Field and interpretive studies conducted by T.C. Winter and U.S. Geological Survey colleagues, and summarized in the 1998 publication “Groundwater and Surface Water – A Single Resource”, inspired a new generation of research centered on extensions of the groundwater-flow code MODFLOW to more sophisticated simulation of coupled groundwater and surface-water systems. Guided by emerging concerns with water availability, safe yields from wells, health of aquatic habitat, and climate change, the changes to MODFLOW involve: 1) the ability to more precisely and accurately represent the interface between surface and subsurface flows and 2) the consideration of a variety of mechanisms that influence their interaction. A review of the most important changes to the code is supplemented in this article by selected case studies in an effort to show the scope of the advances. The updates discussed include the Streamflow Routing (SFR), Lake (LAK), and Unsaturated-Zone Flow (UZF) Packages in MODFLOW-2005 and the Groundwater Management (GWM), Local Grid Refinement (LGR), and Newton (NWT) formulation versions of MODFLOW-2005. New developments feature the integration of rainfall-runoff modeling with MODFLOW in GSFLOW, coupling of GFLOW and MODFLOW in a hybrid code, and the forthcoming unstructured grid version of MODFLOW. They promise continued advances in the ability to use science to protect the single water resource.

Acque Sotterranee: Italian Journal of Groundwater↗

Testing the suitability of geologic frameworks for extrapolating hydraulic properties across regional scales

The suitability of geologic frameworks for extrapolating hydraulic conductivity ( K ) to length scales commensurate with hydraulic data is difficult to assess. A novel method is presented for evaluating assumed relations between K and geologic interpretations for regional-scale groundwater modeling. The approach relies on simultaneous interpretation of multiple aquifer tests using alternative geologic frameworks of variable complexity, where each framework is incorporated as prior information that assumes homogeneous K within each model unit. This approach is tested at Pahute Mesa within the Nevada National Security Site (USA), where observed drawdowns from eight aquifer tests in complex, highly faulted volcanic rocks provide the necessary hydraulic constraints. The investigated volume encompasses 40 mi 3 (167 km 3 ) where drawdowns traversed major fault structures and were detected more than 2 mi (3.2 km) from pumping wells. Complexity of the five frameworks assessed ranges from an undifferentiated mass of rock with a single unit to 14 distinct geologic units. Results show that only four geologic units can be justified as hydraulically unique for this location. The approach qualitatively evaluates the consistency of hydraulic property estimates within extents of investigation and effects of geologic frameworks on extrapolation. Distributions of transmissivity are similar within the investigated extents irrespective of the geologic framework. In contrast, the extrapolation of hydraulic properties beyond the volume investigated with interfering aquifer tests is strongly affected by the complexity of a given framework. Testing at Pahute Mesa illustrates how this method can be employed to determine the appropriate level of geologic complexity for large-scale groundwater modeling.

Nevada↗

Parametric study of the physical properties of hydrate-bearing sand, silt, and clay sediments: 1. Electromagnetic properties

The marked decrease in bulk electrical conductivity of sediments in the presence of gas hydrates has been used to interpret borehole electrical resistivity logs and, to a lesser extent, the results of controlled source electromagnetic surveys to constrain the spatial distribution and predicted concentration of gas hydrate in natural settings. Until now, an exhaustive laboratory data set that could be used to assess the impact of gas hydrate on the electromagnetic properties of different soils (sand, silt, and clay) at different effective stress and with different saturations of hydrate has been lacking. The laboratory results reported here are obtained using a standard geotechnical cell and the hydrate-formed tetrahydrofuran (THF), a liquid that is fully miscible in water and able to produce closely controlled saturations of hydrate from dissolved phase. Both permittivity and electrical conductivity are good indicators of the volume fraction of free water in the sediment, which is in turn dependent on hydrate saturation. Permittivity in the microwave frequency range is particularly predictive of free water content since it is barely affected by ionic concentration, pore structure, and surface conduction. Electrical conductivity (or resistivity) is less reliable for constraining water content or hydrate saturation: In addition to fluid-filled porosity, other factors, such as the ionic concentration of the pore fluid and possibly other conduction effects (e.g., surface conduction in high specific surface soils having low conductivity pore fluid), also influence electrical conductivity.

Journal of Geophysical Research↗

Seismicity and shear strain in the southern Great Basin of Nevada and California

This study examines the relationship between the distribution of small earthquakes (M L ≤4.3) and mechanisms of strain accumulation and relaxation in an area with long repeat times between large events, the Southern Great Basin Seismic Network (SGBSN) region. The Great Basin is a unique continental extensional province characterized by normal and strike-slip faulting, high heat flow, crust of thin to normal thickness, and high elevations. The SGBSN is operated to provide data to address suitability issues pertaining to Yucca Mountain, Nevada which is being evaluated as a potential site for a national mined geologic nuclear waste repository. Suitability issues include estimation of the probability of occurrence of future damaging earthquakes, the characterization of the mechanisms that drive hydrologic flow, and the identification of fractures (faults) that might act as flow conduits or barriers. This study attempts to explain the distribution of small earthquakes in terms of spatial variations in the shear strain field; where strain concentrates there should be a greater number of small earthquakes. Strain field models are constructed under the assumption that long term fault behavior perturbs an otherwise uniform strain field. These strain field models are then interpreted with regard to the regional tectonics and site suitability issues. Modeling results provide one possible explanation of why earthquake clusters cover regions much larger than the surface projections of any of mapped major faults; clusters in a wide band along and extending northeast of the northern half of the Furnace Creek fault may correspond to elevated shear strains along the fault and a broad cluster in the Pahranagat Shear Zone may be associated with shear strain arising from a distribution of smaller localized faults. The relatively large number of small earthquakes in the southern and eastern portions of the Nevada Test Site is consistent with the strain field models. A minimum in shear strain at Yucca Mountain is predicted by all models consistent with an almost total lack of earthquakes observed there. The region to the west of the Death Valley/Furnace Creek fault system, the portion of the study area with the most active deformation but few small earthquakes, is an area of low shear strain. A possible reason for this is that the fault configuration in the area is optimal for accommodating regional deformation via large earthquakes or creep. While there is also a relative lack of earthquakes at Yucca Mountain, this may be indicative of a lack of accumulating strain energy and thus, a lower potential for a large earthquake.

Journal of Geophysical Research↗

Alteration in the Madera limestone and Sandia formation from core hole VC-1, Valles Caldera, New Mexico

Core hole VC-1 penetrated the southwestern ring fracture zone of the 1.1 Ma Valles caldera and at a depth of 333 m intersected the top of the Paleozoic section including the Abo Formation, Madera Limestone, and Sandia Formation, reaching a total depth of 856 m. The Paleozoic rocks, which consist of thin-bedded limestone, siltstone, mudstone, sandstone, and local conglomerate, are overlain by volcanic rocks of the caldera moat that are less than 0.6 Ma. Diagenetic and at least three hydrothermal alteration stages were identified in the Madera Limestone and Sandia Formation. Diagenetic clay alteration was pervasive throughout the sedimentary rocks. Volcanic activity at 16.5 Ma and continuing through the formation of the Valles caldera resulted in high thermal gradients, which caused recrystallization of diagenetic clay minerals. Interstratified smectite-illite is the most diagnostic clay mineral throughout the section; structurally, the illite component in the ordered interstratified illite-smectite changes gradationally from 70% at the top of the Madera Limestone to 95% at the base of the section in the Sandia Formation. Pyrite that occurs as small clots and lenses as well as finely disseminated is interpreted as being of diagenetic origin, especially in organic-rich beds. Low permeability of much of the Paleozoic section precluded the deposition of hydrothermal minerals except in fractures and intergranular space in some of the more permeable sandstone and brecciated horizons. Three stages of hydrothermal mineral deposition are defined. Stage I is widespread and includes mainly chlorite, calcite, pyrite, and interstratified smectite-illite that was formed prior to caldera development at temperatures approximating 200°C. Stage II is characterized by quartz, sericite, and scarce sulfides deposited locally by fluids at approximately 275°C shortly before or at the time of early caldera-related volcanism. The Stage III hydrothermal event, associated with formation of the Valles caldera, caused complex brecciation accompanied by deposition of phengite and molybdenite at temperatures locally approaching 300°C mostly in rocks underlying the Madera Limestone. During the present drilling program, thermal aquifers were encountered at 483 m (113°C), 533 m (125°C), and 845 m (184°C), very close to the Stage II and III mineralized zones; however, measured temperatures of the present hydrothermal system are much lower than during the main hydrothermal events. Good permeability has been maintained in these three levels for a long time while the rest of the Madera Limestone and Sandia Formation was sealed by hydrothermal mineral precipitation.

Journal of Geophysical Research Solid Earth↗