Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geochemistry, Geophysics, Geosystems”

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 109 records · Page 6Linked to original sources
USGS2011

Cooling rates and the depth of detachment faulting at oceanic core complexes: Evidence from zircon Pb/U and (U-Th)/He ages

Oceanic detachment faulting represents a distinct mode of seafloor spreading at slow spreading mid‐ocean ridges, but many questions persist about the thermal evolution and depth of faulting. We present new Pb/U and (U‐Th)/He zircon ages and combine them with magnetic anomaly ages to define the cooling histories of gabbroic crust exposed by oceanic detachment faults at three sites along the Mid‐Atlantic Ridge (Ocean Drilling Program (ODP) holes 1270D and 1275D near the 15°20′N Transform, and Atlantis Massif at 30°N). Closure temperatures for the Pb/U (∼800°C–850°C) and (U‐Th)/He (∼210°C) isotopic systems in zircon bracket acquisition of magnetic remanence, collectively providing a temperature‐time history during faulting. Results indicate cooling to ∼200°C in 0.3–0.5 Myr after zircon crystallization, recording time‐averaged cooling rates of ∼1000°C–2000°C/Myr. Assuming the footwalls were denuded along single continuous faults, differences in Pb/U and (U‐Th)/He zircon ages together with independently determined slip rates allow the distance between the ∼850°C and ∼200°C isotherms along the fault plane to be estimated. Calculated distances are 8.4 ± 4.2 km and 5.0 ± 2.1 km from holes 1275D and 1270D and 8.4 ± 1.4 km at Atlantis Massif. Estimating an initial subsurface fault dip of 50° and a depth of 1.5 km to the 200°C isotherm leads to the prediction that the ∼850°C isotherm lies ∼5–7 km below seafloor at the time of faulting. These depth estimates for active fault systems are consistent with depths of microseismicity observed beneath the hypothesized detachment fault at the TAG hydrothermal field and high‐temperature fault rocks recovered from many oceanic detachment faults.

Geochemistry, Geophysics, Geosystems↗
USGS2011

Modeling sulfate reduction in methane hydrate-bearing continental margin sediments: Does a sulfate-methane transition require anaerobic oxidation of methane?

The sulfate‐methane transition (SMT), a biogeochemical zone where sulfate and methane are metabolized, is commonly observed at shallow depths (1–30 mbsf) in methane‐bearing marine sediments. Two processes consume sulfate at and above the SMT, anaerobic oxidation of methane (AOM) and organoclastic sulfate reduction (OSR). Differentiating the relative contribution of each process is critical to estimate methane flux into the SMT, which, in turn, is necessary to predict deeper occurrences of gas hydrates in continental margin sediments. To evaluate the relative importance of these two sulfate reduction pathways, we developed a diagenetic model to compute the pore water concentrations of sulfate, methane, and dissolved inorganic carbon (DIC). By separately tracking DIC containing 12 C and 13 C, the model also computes δ 13 C‐DIC values. The model reproduces common observations from methane‐rich sediments: a well‐defined SMT with no methane above and no sulfate below and a δ 13 C‐DIC minimum at the SMT. The model also highlights the role of upward diffusing 13 C‐enriched DIC in contributing to the carbon isotope mass balance of DIC. A combination of OSR and AOM, each consuming similar amounts of sulfate, matches observations from Site U1325 (Integrated Ocean Drilling Program Expedition 311, northern Cascadia margin). Without AOM, methane diffuses above the SMT, which contradicts existing field data. The modeling results are generalized with a dimensional analysis to the range of SMT depths and sedimentation rates typical of continental margins. The modeling shows that AOM must be active to establish an SMT wherein methane is quantitatively consumed and the δ 13 C‐DIC minimum occurs. The presence of an SMT generally requires active AOM.

Geochemistry, Geophysics, Geosystems↗
USGS2011

Three-dimensional model of an ultramafic feeder system to the Nikolai Greenstone mafic large igneous province, central Alaska Range

The Amphitheater Mountains and southern central Alaska Range expose a thick sequence of Triassic Nikolai basalts that is underlain by several mafic‐ultramafic complexes, the largest and best exposed being the Fish Lake and Tangle (FL‐T) mafic‐ultramafic sills that flank the Amphitheater Mountains synform. Three‐dimensional (3‐D) modeling of gravity and magnetic data reveals details of the structure of the Amphitheater Mountains, such as the orientation and thickness of Nikolai basalts, and the geometry of the FL‐T intrusions. The 3‐D model (50 × 70 km) includes the full geographic extent of the FL‐T complexes and consists of 11 layers. Layer surfaces and properties (density and magnetic susceptibility) were modified by forward and inverse methods to reduce differences between the observed and calculated gravity and magnetic grids. The model suggests that the outcropping FL‐T sills are apparently connected and traceable at depth and reveals variations in thickness, shape, and orientation of the ultramafic bodies that may identify paths of magma flow. The model shows that a significant volume (2000 km 3 ) of ultramafic material occurs in the subsurface, gradually thickening and plunging westward to depths exceeding 4 km. This deep ultramafic material is interpreted as the top of a keel or root system that supplied magma to the Nikolai lavas and controlled emplacement of related magmatic intrusions. The presence of this deep, keel‐like structure, and asymmetry of the synform, supports a sag basin model for development of the Amphitheater Mountains structure and reveals that the feeders to the Nikolai are much more extensive than previously known.

Alaska↗
USGS2011

Tiny intraplate earthquakes triggered by nearby episodic tremor and slip in Cascadia

Episodic tremor and slip (ETS) has been observed in many subduction zones, but its mechanical underpinnings as well as its potential for triggering damaging earthquakes have proven difficult to assess. Here we use a seismic array in Cascadia of unprecedented density to monitor seismicity around a moderate 16 day ETS episode. In the 4 months of data we examine, we observe five tiny earthquakes within the subducting slab during the episode and only one more in the same area, which was just before and nearby the next ETS burst. These earthquakes concentrate along the sides and updip edge of the ETS region, consistent with greater stress concentration there than near the middle and downdip edge of the tremor area. Most of the seismicity is below the megathrust, with a similar depth extent to the background intraslab seismicity. The pattern of earthquakes that we find suggests slow slip has a more continuous temporal and spatial pattern than the tremor loci, which notoriously appear in bursts, jumps, and streaks. Copyright 2011 by the American Geophysical Union.

Geochemistry, Geophysics, Geosystems↗
USGS2011

Gravity lineaments of the Cocos Plate: Evidence for a thermal contraction crack origin

Lineaments in the gravity field with wavelengths of 100–200 km affect the south-central Pacific. Because they align with absolute plate motion, it has been proposed that they reflect small-scale convection cells beneath the lithosphere that become elongated by basal shear. Alternatively, it was suggested that they reflect channelized flow of low viscosity material following the base of the lithosphere toward the East Pacific Rise, or that they result from lithospheric-scale thermal contraction cracks. Here, we report about previously undetected gravity lineaments across the Cocos Plate. Similarly to the south-central Pacific lineaments, the Cocos lineaments affect a plate that is anomalously shallow, with seamounts aligning mostly within their troughs. However, the Cocos lineaments strike markedly oblique to absolute plate motion and follow instead trajectories that are perpendicular to seafloor isochrons, a characteristic best explained by the thermal contraction crack model. The presence of steep scarps at the base of seamounts and the seismic imaging of faults striking perpendicular to isochrons further support this interpretation. Assuming that the slow subsidence rates of the south-central Pacific and Cocos plates reflect a warmer upper mantle, we propose that the associated thinner elastic plates favor the formation of thermal contraction cracks. A thinner elastic plate may also explain the pattern of ridge propagation in both areas. At large ridge offsets with a history of steady migration, the propagating segments have been those cutting into the shallower flanks, consistent with the concept that a warmer, thinner plate is more easily cracked.

Geochemistry Geophysics Geosystems↗
USGS2011

Magnetic properties in an ash flow tuff with continuous grain size variation: a natural reference for magnetic particle granulometry

The Tiva Canyon Tuff contains dispersed nanoscale Fe-Ti-oxide grains with a narrow magnetic grain size distribution, making it an ideal material in which to identify and study grain-size-sensitive magnetic behavior in rocks. A detailed magnetic characterization was performed on samples from the basal 5 m of the tuff. The magnetic materials in this basal section consist primarily of (low-impurity) magnetite in the form of elongated submicron grains exsolved from volcanic glass. Magnetic properties studied include bulk magnetic susceptibility, frequency-dependent and temperature-dependent magnetic susceptibility, anhysteretic remanence acquisition, and hysteresis properties. The combined data constitute a distinct magnetic signature at each stratigraphic level in the section corresponding to different grain size distributions. The inferred magnetic domain state changes progressively upward from superparamagnetic grains near the base to particles with pseudo-single-domain or metastable single-domain characteristics near the top of the sampled section. Direct observations of magnetic grain size confirm that distinct transitions in room temperature magnetic susceptibility and remanence probably denote the limits of stable single-domain behavior in the section. These results provide a unique example of grain-size-dependent magnetic properties in noninteracting particle assemblages over three decades of grain size, including close approximations of ideal Stoner-Wohlfarth assemblages, and may be considered a useful reference for future rock magnetic studies involving grain-size-sensitive properties.

Nevada↗
USGS2011

Shallow lithological structure across the Dead Sea Transform derived from geophysical experiments

In the framework of the DEad SEa Rift Transect (DESERT) project a 150 km magnetotelluric profile consisting of 154 sites was carried out across the Dead Sea Transform. The resistivity model presented shows conductive structures in the western section of the study area terminating abruptly at the Arava Fault. For a more detailed analysis we performed a joint interpretation of the resistivity model with a P wave velocity model from a partially coincident seismic experiment. The technique used is a statistical correlation of resistivity and velocity values in parameter space. Regions of high probability of a coexisting pair of values for the two parameters are mapped back into the spatial domain, illustrating the geographical location of lithological classes. In this study, four regions of enhanced probability have been identified, and are remapped as four lithological classes. This technique confirms the Arava Fault marks the boundary of a highly conductive lithological class down to a depth of ∼3 km. That the fault acts as an impermeable barrier to fluid flow is unusual for large fault zone, which often exhibit a fault zone characterized by high conductivity and low seismic velocity. At greater depths it is possible to resolve the Precambrian basement into two classes characterized by vastly different resistivity values but similar seismic velocities. The boundary between these classes is approximately coincident with the Al Quweira Fault, with higher resistivities observed east of the fault. This is interpreted as evidence for the original deformation along the DST originally taking place at the Al Quweira Fault, before being shifted to the Arava Fault.

Dead Sea↗
USGS2011

Origin of minor and trace element compositional diversity in anorthitic feldspar phenocrysts and melt inclusions from the Juan de Fuca Ridge

Melt inclusions trapped in phenocryst phases are important primarily due to their potential of preserving a significant proportion of the diversity of magma composition prior to modification of the parent magma array during transport through the crust. The goal of this investigation was to evaluate the impact of formational and post-entrapment processes on the composition of melt inclusions hosted in high anorthite plagioclase in MORB. Our observations from three plagioclase ultra-phyric lavas from the Endeavor Segment of the Juan de Fuca Ridge document a narrow range of major elements and a dramatically greater range of minor and trace elements within most host plagioclase crystals. Observed host/inclusion partition coefficients for Ti are consistent with experimental determinations. In addition, observed values of D Ti are independent of inclusion size and inclusion TiO 2 content of the melt inclusion. These observations preclude significant effects from the re-homogenization process, entrapment of incompatible element boundary layers or dissolution/precipitation. The observed wide range of TiO 2 contents in the host feldspar, and between bands of melt inclusions within individual crystals rule out modification of TiO contents by diffusion, either pre-eruption or due to re-homogenization. However, we do observe comparatively small ranges for values of K 2 O and Sr compared to P 2 O 5 and TiO 2 in both inclusions and crystals that can be attributed to diffusive processes that occurred prior to eruption.

Geochemistry, Geophysics, Geosystems↗
USGS2012

Deep magmatic degassing versus scrubbing: Elevated CO2 emissions and C/S in the lead-up to the 2009 eruption of Redoubt Volcano, Alaska

We report CO 2, SO 2 , and H 2 S emission rates and C/S ratios during the five months leading up to the 2009 eruption of Redoubt Volcano, Alaska. CO 2 emission rates up to 9018 t/d and C/S ratios ≥30 measured in the months prior to the eruption were critical for fully informed forecasting efforts. Observations of ice-melt rates, meltwater discharge, and water chemistry suggest that surface waters represented drainage from surficial, perched reservoirs of condensed magmatic steam and glacial meltwater. These fluids scrubbed only a few hundred tonnes/day of SO 2 , not the >2100 t/d SO 2 expected from degassing of magma in the mid- to upper crust (3–6.5 km), where petrologic analysis shows the final magmatic equilibration occurred. All data are consistent with upflow of a CO 2 -rich magmatic gas for at least 5 months prior to eruption, and minimal scrubbing of SO 2 by near-surface groundwater. The high C/S ratios observed could reflect bulk degassing of mid-crustal magma followed by nearly complete loss of SO 2 in a deep magmatic-hydrothermal system. Alternatively, high C/S ratios could be attributed to decompressional degassing of low silica andesitic magma that intruded into the mid-crust in the 5 months prior to eruption, thereby mobilizing the pre-existing high silica andesite magma or mush in this region. The latter scenario is supported by several lines of evidence, including deep long-period earthquakes (−28 to −32 km) prior to and during the eruption, and far-field deformation following the onset of eruptive activity.

Alaska↗
USGS2012

Generation and evolution of hydrothermal fluids at Yellowstone: Insights from the Heart Lake Geyser Basin

We sampled fumaroles and hot springs from the Heart Lake Geyser Basin (HLGB), measured water and gas discharge, and estimated heat and mass flux from this geothermal area in 2009. The combined data set reveals that diverse fluids share an origin by mixing of deep solute-rich parent water with dilute heated meteoric water, accompanied by subsequent boiling. A variety of chemical and isotopic geothermometers are consistent with a parent water that equilibrates with rocks at 205°C ± 10°C and then undergoes 21% ± 2% adiabatic boiling. Measured diffuse CO 2 flux and fumarole compositions are consistent with an initial dissolved CO 2 concentration of 21 ± 7 mmol upon arrival at the caldera boundary and prior to southeast flow, boiling, and discharge along the Witch Creek drainage. The calculated advective flow from the basin is 78 ± 16 L s −1 of parent thermal water, corresponding to 68 ± 14 MW, or –1% of the estimated thermal flux from Yellowstone. Helium and carbon isotopes reveal minor addition of locally derived crustal, biogenic, and meteoric gases as this fluid boils and degasses, reducing the He isotope ratio (Rc/Ra) from 2.91 to 1.09. The HLGB is one of the few thermal areas at Yellowstone that approaches a closed system, where a series of progressively boiled waters can be sampled along with related steam and noncondensable gas. At other Yellowstone locations, steam and gas are found without associated neutral Cl waters (e.g., Hot Spring Basin) or Cl-rich waters emerge without significant associated steam and gas (Upper Geyser Basin).

Wyoming↗
USGS2012

Cenozoic tectono-thermal history of the Tordrillo Mountains, Alaska: Paleocene-Eocene ridge subduction, decreasing relief, and late Neogene faulting

Topographic development inboard of the continental margin is a predicted response to ridge subduction. New thermochronology results from the western Alaska Range document ridge subduction related orogenesis. K-feldspar thermochronology (KFAT) of bedrock samples from the Tordrillo Mountains in the western Alaska Range complement existing U-Pb, 40 Ar/ 39 Ar and AFT (apatite fission track) data to provide constraints on Paleocene pluton emplacement, and cooling as well as Late Eocene to Miocene vertical movements and exhumation along fault-bounded blocks. Based on the KFAT analysis we infer rapid exhumation-related cooling during the Eocene in the Tordrillo Mountains. Our KFAT cooling ages are coeval with deposition of clastic sediments in the Cook Inlet, Matanuska Valley and Tanana basins, which reflect high-energy depositional environments. The Tordrillo Mountains KFAT cooling ages are also the same as cooling ages in the Iliamna Lake region, the Kichatna Mountains of the western Alaska Range, and Mt. Logan in the Wrangell-St. Elias Mountains, thus rapid cooling at this time encompasses a broad region inboard of, and parallel to, the continental margin extending for several hundred kilometers. We infer these cooling events and deposition of clastic rocks are related to thermal effects that track the eastward passage of a slab window in Paleocene-Eocene time related to the subduction of the proposed Resurrection-Kula spreading ridge. In addition, we conclude that the reconstructed KFAT max negative age-elevation relationship is likely related to a long period of decreasing relief in the Tordrillo Mountains.

Alaska↗
USGS2012

Tracking lava flow emplacement on the east rift zone of Kilauea, Hawai’i with synthetic aperture radar (SAR) coherence

Lava flow mapping is both an essential component of volcano monitoring and a valuable tool for investigating lava flow behavior. Although maps are traditionally created through field surveys, remote sensing allows an extraordinary view of active lava flows while avoiding the difficulties of mapping on location. Synthetic aperture radar (SAR) imagery, in particular, can detect changes in a flow field by comparing two images collected at different times with SAR coherence. New lava flows radically alter the scattering properties of the surface, making the radar signal decorrelated in SAR coherence images. We describe a new technique, SAR Coherence Mapping (SCM), to map lava flows automatically from coherence images independent of look angle or satellite path. We use this approach to map lava flow emplacement during the Pu‘u ‘Ō‘ō-Kupaianaha eruption at Kīlauea, Hawai‘i. The resulting flow maps correspond well with field mapping and better resolve the internal structure of surface flows, as well as the locations of active flow paths. However, the SCM technique is only moderately successful at mapping flows that enter vegetation, which is also often decorrelated between successive SAR images. Along with measurements of planform morphology, we are able to show that the length of time a flow stays decorrelated after initial emplacement is linearly related to the flow thickness. Finally, we use interferograms obtained after flow surfaces become correlated to show that persistent decorrelation is caused by post-emplacement flow subsidence.

Hawai'i↗
USGS2012

Integrated geophysical and hydrothermal models of flank degassing and fluid flow at Masaya Volcano, Nicaragua

We investigate geologic controls on circulation in the shallow hydrothermal system of Masaya volcano, Nicaragua, and their relationship to surface diffuse degassing. On a local scale (~250 m), relatively impermeable normal faults dipping at ~60° control the flowpath of water vapor and other gases in the vadose zone. These shallow normal faults are identified by modeling of a NE-SW trending magnetic anomaly of up to 2300 nT that corresponds to a topographic offset. Elevated SP and CO2 to the NW of the faults and an absence of CO2 to the SE suggest that these faults are barriers to flow. TOUGH2 numerical models of fluid circulation show enhanced flow through the footwalls of the faults, and corresponding increased mass flow and temperature at the surface (diffuse degassing zones). On a larger scale, TOUGH2 modeling suggests that groundwater convection may be occurring in a 3-4 km radial fracture zone transecting the entire flank of the volcano. Hot water rising uniformly into the base of the model at 1 x 10-5 kg/m2s results in convection that focuses heat and fluid and can explain the three distinct diffuse degassing zones distributed along the fracture. Our data and models suggest that the unusually active surface degassing zones at Masaya volcano can result purely from uniform heat and fluid flux at depth that is complicated by groundwater convection and permeability variations in the upper few km. Therefore isolating the effects of subsurface geology is vital when trying to interpret diffuse degassing in light of volcanic activity.

Masaya Volcano↗
USGS2012

An apparent "vital effect" of calcification rate on the Sr/Ca temperature proxy in the reef coral Montipora capitata

Measuring the strontium to calcium ratio in coral skeletons reveals information on seawater temperatures during skeletal deposition, but studies have shown additional variables may affect the ratio. Here we measured Sr/Ca in the reef coral, Montipora capitata, grown in six mesocosms continuously supplied with seawater from the adjacent reef flat. Three mesocosms were ambient controls, and three had seawater chemistry simulating "ocean acidification" (OA). We found that Sr/Ca was not affected by the OA treatment, and neither was coral calcification for these small colonies (larger colonies did show an OA effect). The lack of OA effects allowed us to test the hypothesis that coral growth rate can affect Sr/Ca using the natural range in calcification rates of the corals grown at the same temperature. We found that Sr/Ca was inversely related to calcification rate (Sr/Ca = 9.39 - 0.00404 mmol/mol * mg day-1 cm-2, R2 = 0.32). Using a previously published calibration curve for this species, a 22 mg day-1 colony-1 increase in calcification rate introduced a 1°C warmer temperature estimate, with the 27 corals reporting "temperatures" ranging from 24.9 to 28.9, with mean 26.6 ± 0.9°C SD. Our results lend support to hypotheses invoking kinetic processes and growth rate to explain vital effects on Sr/Ca. However, uncertainty in the slope of the regression of Sr/Ca on calcification and a low R-squared value lead us to conclude that Sr/Ca could still be a useful proxy in this species given sufficient replication or by including growth rate in the calibration.

Geochemistry, Geophysics, Geosystems↗
USGS2012

Illumination of rheological mantle heterogeneity by the M7.2 2010 El Mayor-Cucapah earthquake

Major intracontinental strike-slip faults tend to mark boundaries between lithospheric blocks of contrasting mechanical properties along much of their length. Both crustal and mantle heterogeneities can form such boundaries, but the role of crustal versus mantle strength contrasts for localizing strain sufficiently to generate major faults remains unclear. Using the crustal velocity field observed through the Global Positioning System (GPS) in the epicentral area of the M7.2 2010 El Mayor-Cucapah earthquake, Baja California, we find that transient deformation observed after the event is anomalously small in areas of relatively high seismic velocity in the shallow upper mantle (∼50 km depth). This pattern is best explained with a laterally heterogeneous viscoelastic structure that mimics the seismic structure. The mantle of the Southern Colorado River Desert (SCRD) and Peninsular Ranges (PR), which bound the fault system to its east and west, respectively, have anomalously high viscosity and seismic velocity. We hypothesize that compared with the rest of the San Andreas fault (SAF) system to its north, the strike-slip fault system in northern Baja California is narrow because of the presence of the PR and SCRD high-viscosity regions which bound it.

Arizona, Baja California, California, Sonora↗
USGS2012

Copper-nickel-rich, amalgamated ferromanganese crust-nodule deposits from Shatsky Rise, NW Pacific

A unique set of ferromanganese crusts and nodules collected from Shatsky Rise (SR), NW Pacific, were analyzed for mineralogical and chemical compositions, and dated using Be isotopes and cobalt chronometry. The composition of these midlatitude, deep-water deposits is markedly different from northwest-equatorial Pacific (PCZ) crusts, where most studies have been conducted. Crusts and nodules on SR formed in close proximity and some nodule deposits were cemented and overgrown by crusts, forming amalgamated deposits. The deep-water SR crusts are high in Cu, Li, and Th and low in Co, Te, and Tl concentrations compared to PCZ crusts. Thorium concentrations (ppm) are especially striking with a high of 152 (mean 56), compared to PCZ crusts (mean 11). The deep-water SR crusts show a diagenetic chemical signal, but not a diagenetic mineralogy, which together constrain the redox conditions to early oxic diagenesis. Diagenetic input to crusts is rare, but unequivocal in these deep-water crusts. Copper, Ni, and Li are strongly enriched in SR deep-water deposits, but only in layers older than about 3.4 Ma. Diagenetic reactions in the sediment and dissolution of biogenic calcite in the water column are the likely sources of these metals. The highest concentrations of Li are in crust layers that formed near the calcite compensation depth. The onset of Ni, Cu, and Li enrichment in the middle Miocene and cessation at about 3.4 Ma were accompanied by changes in the deep-water environment, especially composition and flow rates of water masses, and location of the carbonate compensation depth.

Geochemistry, Geophysics, Geosystems↗
USGS2012

Temporal variations of geyser water chemistry in the Upper Geyser Basin, Yellowstone National Park, USA

Geysers are rare features that reflect a delicate balance between an abundant supply of water and heat and a unique geometry of fractures and porous rocks. Between April 2007 and September 2008, we sampled Old Faithful, Daisy, Grand, Oblong, and Aurum geysers in Yellowstone National Park's Upper Geyser Basin and characterized temporal variations in major element chemistry and water isotopes (δ 18 O, δD, 3 H). We compare these temporal variations with temporal trends of Geyser Eruption Intervals (GEI). SiO 2 concentrations and geothermometry indicate that the geysers are fed by waters ascending from a reservoir with temperatures of ∼190 to 210°C. The studied geysers display small and complex chemical and isotopic seasonal variations, and geysers with smaller volume display larger seasonal variations than geysers with larger volumes. Aurum and Oblong Geysers contain detectable tritium concentrations, suggesting that erupted water contains some modern meteoric water. We propose that seasonal GEI variations result from varying degrees of evaporation, meteoric water recharge, water table fluctuations, and possible hydraulic interaction with the adjacent Firehole River. We demonstrate that the concentrations of major dissolved species in Old Faithful Geyser have remained nearly constant since 1884 despite large changes in Old Faithful's eruption intervals, suggesting that no major changes have occurred in the hydrothermal system of the Upper Geyser Basin for >120 years. Our data set provides a baseline for monitoring future changes in geyser activity that might result from varying climate, earthquakes, and changes in heat flow from the underlying magmatic system.

Wyoming↗
USGS2012

Large-scale splay faults on a strike-slip fault system: The Yakima Folds, Washington State

The Yakima Folds (YF) comprise anticlines above reverse faults cutting flows of the Miocene Columbia River Basalt Group of central Washington State. The YF are bisected by the ~1100-km-long Olympic-Wallowa Lineament (OWL), which is an alignment of topographic features including known faults. There is considerable debate about the origin and earthquake potential of both the YF and OWL, which lie near six major dams and a large nuclear waste storage site. Here I show that the trends of the faults forming the YF relative to the OWL match remarkably well the trends of the principal stress directions at the end of a vertical strike-slip fault. This comparison and the termination of some YF against the OWL are consistent with the YF initially forming as splay faults caused by an along-strike decrease in the amount of strike-slip on the OWL. The hypothesis is that the YF faults initially developed as splay faults in the early to mid Miocene under NNW-oriented principal compressive stress, but the anticlines subsequently grew with thrust motion after the principal compressive stress direction rotated to N-S or NNE after the mid-Miocene. A seismic profile across one of the YF anticlines shows folding at about 7 km depth, indicating deformation of sub-basalt strata. The seismic profile and the hypothesized relationship between the YF and the OWL suggest that the structures are connected in the middle or lower crust, and that the faults forming the YF are large-scale splay faults associated with a major strike-slip fault system.

Washington↗