Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research: Solid Earth”

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 217 records · Page 12Linked to original sources
USGS2012

Geochemistry, mineralogy, and petrogenesis of basalt from the Gorda Ridge

Basalt pillow lava with glassy rims was dredged from 17 sites along the Gorda Ridge between latitudes 43°N and 40.8°N. All samples are low-K 2 O mid-ocean ridge tholeiitic basalt having a narrow compositional range with Mg numbers clustered around 60–62; more primitive and evolved compositions are present but rare. None of the more fractionated samples is as evolved as ferrobasalt from the Juan de Fuca Ridge. Incompatible trace element ratios are typical of normal mid-ocean ridge basalt but show considerable scatter. The 87 Sr/ 86 Sr ratios, ranging from 0.70233 to 0.70267, are near the low end of the range for normal mid-ocean ridge basalt. The phenocrysts present are abundant plagioclase, less abundant olivine ± spinel, and traces of sulfide. Clinopyroxene is absent in the glasses. The mineral chemistry is typical of ocean floor basalt. However, unusually anorthitic plagioclase in some highly plagioclasephyric samples, complex zoning in plagioclase and olivine, and other disequilibrium features indicate complex petrogenetic processes involving magma mixing. Glass inclusions in olivine and plagioclase indicate primitive trapped melts. The chemical diversity of Gorda Ridge lava appears to have been primarily caused by variable degrees of partial melting, but crystal fractionation and magma mixing were also important processes. Compositional variations along the strike of the ridge indicate that lava from the segment south of the offset at latitude 41.6°N was generated by the smallest percentage of partial melting of a source less depleted than that of the two northern segments. Lava from the central and northern segments was generated by variable percentages of partial melting of a uniformly depleted source. Lava from the northern segment shows the greatest diversity with a large range in percentage of melting and both relatively primitive and evolved compositions in spatially and temporally closely associated lava. The diversity of lava from the northern segment may result from thermal perturbations caused by proximity to the Blanco Fracture Zone. Gorda Ridge lava is more similar to that from the Mid-Atlantic Ridge at 22°–25°N than to that from other spreading centers in the Pacific Ocean. Magma chambers under the Gorda Ridge appear to be small, discontinuous in time and space and possibly at greater depth than those beneath other Pacific spreading centers.

Journal of Geophysical Research Solid Earth↗
USGS2012

Uniformity and diversity in the composition of mineralizing fluids from hydrothermal vents on the southern Juan de Fuca Ridge

Abundances of Li, Na, K, Rb, Ca, Sr, Ba, Mn, Fe, Zn, and Si have been determined in fluid samples from seven vents located in three areas on the southern Juan de Fuca Ridge. The hydrothermal component estimated from the Mg contents of the samples ranges from 7% to 76%. Concentrations of Fe and Si, among other elements, in acid-stabilized solutions appear to be generally representative of the parental hydrothermal fluids, but some Zn determinations and most Ba values appear to be too low. Thermodynamic calculations indicate that the acidified samples remain supersaturated with respect to silica, barite, and pyrite; unacidified samples are supersaturated, in addition, with respect to ZnS, FeS, and many silicate phases. Within the constraints of limited sampling there appear to be differences in fluid compositions both within and between the three vent areas. Some uniform differences in the elemental abundances predicted for hydrothermal end-member fluids might be due to inmixing of fresh seawater at depth in the hydrothermal system. The Juan de Fuca hydrothermal fluids contain more Fe but otherwise have relative elemental abundances fairly similar to those in 13°N (East Pacific Rise) fluids, albeit at higher levels. In contrast, fluids from 21°N (East Pacific Rise) and Galapagos have lower K/Rb and much lower Sr and Na abundances; these compositional features probably result from interaction of these fluids with a different mineral assemblage, possibly more mature greenstone.

Journal of Geophysical Research Solid Earth↗
USGS2012

Mineralogical studies of sulfide samples and volatile concentrations of basalt glasses from the southern Juan de Fuca Ridge

Sulfide samples obtained from the U.S. Geological Survey's DSRV Alvin dives on the southern Juan de Fuca Ridge closely resemble those from the same area described by Koski et al. (1984). Major minerals include sphalerite, wurtzite, pyrite, marcasite, isocubanite, anhydrite, and chalcopyrite. Equilibrium, if attained at all, during deposition of most sulfides was a transient event over a few tens of micrometers at most and was perturbed by rapid temperature and compositional changes of the circulating fluid. Two new minerals were found: one, a hydrated Zn, Fe hydroxy-chlorosulfate, and the other, a (Mn, Mg, Fe) hydroxide or hydroxy-hydrate. Both were formed at relatively low temperatures. Lizardite, starkeyite, and anatase were found for the first time in such an environment. Sulfide geothermometry involving the system Cu-Fe-S indicates a vent temperature of <328°C for one sample. Fluid inclusion studies on crystals from the same vicinity of the same sample give pressure-corrected homogenization temperatures of 268° and 285°C. Ice-melting temperatures on inclusions from the same sample are about −2.8°C, indicating that the equivalent salinity of the trapped fluid is about 50% greater than that of seawater. Volatile concentrations from vesicle-free basalt glass from the vent field are about 0.013 wt % CO 2 and 0.16 wt % H 2 O. CO 2 contents in these samples yield an entrapment depth of 2200 m of seawater, which is the depth from which the samples were collected.

Journal of Geophysical Research Solid Earth↗
USGS2012

Relation of the spectroscopic reflectance of olivine to mineral chemistry and some remote sensing implications

Using high-resolution visible and near-infrared diffuse spectral reflectance, we have systematically investigated apparent wavelength shifts as a function of mineral chemistry in the Fe/Mg olivine series from Fo 11 to Fo 91 . The study also shows that trace amounts of nickel can be spectrally detected in the olivine structure. We show that significant compositional information can only be extracted at relatively high resolution, because the overall spectral characteristics of the olivines change only subtly as a function of the Fe/Mg ratio. Significant spectral variation as a function of grain size is also demonstrated, adding a further complication to the interpretation of remotely sensed data from olivine-rich surfaces. This laboratory study is expected to aid in the interpretation of remotely sensed data from both terrestrial and extraterrestrial bodies. Terrestrial applications may include the recognition of ultramafic, ultrabasic, and basaltic terrains which in themselves may have mineral potential. Among extraterrestrial applications, the asteroids are obvious candidates for further examination if instrumentation can provide the necessary wavelength coverage, resolution, and signal-to-noise ratio so that spectra can be compared to those laboratory data discussed here. Some permutations of Fe-Mg-Ni relations in olivines are discussed as they apply to the interpretation of asteroid surfaces and other extraterrestrial bodies.

Journal of Geophysical Research Solid Earth↗
USGS2012

Fitting degradation of shoreline scarps by a nonlinear diffusion model

The diffusion model of degradation of topographic features is a promising means by which vertical offsets on Holocene faults might be dated. In order to calibrate the method, we have examined present-day profiles of wave-cut shoreline scarps of late Pleistocene lakes Bonneville and Lahontan. It may be assumed that these scarps were initially at least as steep as the angle of repose. Offsets range from 1 to 12 m, and present slope angles range from 9° to 29°. A parameter called apparent diffusion age, defined as half the mean square horizontal extent of the slope function of each profile, is plotted as a function of scarp offset. The points show a clear trend of apparent age increasing nearly linearly with offset. If linear diffusion held and scarps were initially vertical, apparent diffusion age would be the same for all the shoreline profiles. The increasing trend can only partly be explained by nonvertical initial scarp slope, and therefore the rate of transport of material downslope must increase significantly faster than a linear law in the range of slopes spanned by the data. The transport law must become linear at small slope to reduce scatter between profiles with varying ambient slopes. The transport law adopted for the purpose of dating is k 0 (1 + 5 s 2 ), where s is local slope. The transport coefficient k 0 is correlated inversely with fan slope, suggesting that there is a dependence on the particle size distribution. A table is included that allows easy application of the model to scarps with simple initial shape.

Journal of Geophysical Research Solid Earth↗
USGS2012

Anomalous shear wave attenuation in the shallow crust beneath the Coso volcanic region, California

We use seismograms of local earthquakes to image relative shear wave attenuation structure in the shallow crust beneath the region containing the Coso volcanic-geothermal area of eastern California. SV and P wave amplitudes were measured from vertical component seismograms of earthquakes that occurred in the Cososouthem Sierra Nevada region from July 1983 to August 1985. Seismograms of 16 small earthquakes show SV amplitudes which are greatly diminished at some azimuths and takeoff angles, indicating strong lateral variations in S wave attenuation in the area. Three-dimensional images of the relative S wave attenuation structure are obtained from forward modeling and a back projection inversion of the amplitude data. The results indicate regions within a 20 by 30 by 10 km volume of the shallow crust (one shallower than 5 km) that severely attenuate SV waves passing through them. These anomalies lie beneath the Indian Wells Valley, 30 km south of the Coso volcanic field, and are coincident with the epicentral locations of recent earthquake swarms. No anomalous attenuation is seen beneath the Coso volcanic field above about 5 km depth. Geologic relations and the coincidence of anomalously slow P wave velocities suggest that the attenuation anomalies may be related to magmatism along the eastern Sierra front.

Journal of Geophysical Research Solid Earth↗
USGS2012

Volcanic hazards and public response

Although scientific understanding of volcanoes is advancing, eruptions continue to take a substantial toll of life and property. Some of these losses could be reduced by better advance preparation, more effective flow of information between scientists and public officials, and better understanding of volcanic behavior by all segments of the public. The greatest losses generally occur at volcanoes that erupt infrequently where people are not accustomed to dealing with them. Scientists sometimes tend to feel that the blame for poor decisions in emergency management lies chiefly with officials or journalists because of their failure to understand the threat. However, the underlying problem embraces a set of more complex issues comprising three pervasive factors. The first factor is the volcano: signals given by restless volcanoes are often ambiguous and difficult to interpret, especially at long-quiescent volcanoes. The second factor is people: people confront hazardous volcanoes in widely divergent ways, and many have difficulty in dealing with the uncertainties inherent in volcanic unrest. The third factor is the scientists: volcanologists correctly place their highest priority on monitoring and hazard assessment, but they sometimes fail to explain clearly their conclusions to responsible officials and the public, which may lead to inadequate public response. Of all groups in society, volcanologists have the clearest understanding of the hazards and vagaries of volcanic activity; they thereby assume an ethical obligation to convey effectively their knowledge to benefit all of society. If society resists, their obligation nevertheless remains. They must use the same ingenuity and creativity in dealing with information for the public that they use in solving scientific problems. When this falls short, even excellent scientific results may be nullified.

Journal of Geophysical Research Solid Earth↗
USGS2012

The mechanics and three-dimensional internal structure of active magmatic systems: Kilauea volcano, Hawaii

Interpretation of abundant seismic data suggests that Kilauea's primary conduit within the upper mantle is concentrically zoned to about 34-km depth. This zoned structure is inferred to contain a central core region of relatively higher permeability, surrounded by numerous dikes that are in intermittent hydraulic communication with each other and with the central core. During periods of relatively high magma transport, the entire cross section of the conduit is utilized. During periods of relatively low to moderate transport, however, only the central core is active. As the conduit penetrates the oceanic crust and enters the volcanic shield, it simultaneously supplies the deeper sections of the rift zones (6-to 10-km depth) and the roots of the summit reservoir with picritic magma. The rift zones at depth are inferred to be almost wholly molten and to possess a high degree of fluid continuity from Heiheiahulu in the East Rift Zone, 45 km westward through the roots of the summit magma reservoir, and well into the Southwest Rift Zone. Higher in the shield, the subcaldera magma reservoir and the shallow rift zones occupy the 2-to 4-km depth interval. Summit-differentiated olivine tholeiite (ρ ≈ 2.62 g cm −3 ) is periodically injected laterally along a horizon of neutral buoyancy within the rift zones, where the density of the magma is just balanced by the in situ density of the shield (Ryan, 1987 a , b ). Deep rift zone intrusions push seaward the deep tectonic blocks of the volcano's south flank. Shallow rift intrusions build a sheeted dike complex, inferred to be in isostatic equilibrium with the higher-density deep rift cores below. General finite element analyses are presented for the deformation and stress fields surrounding such dikes in the horizontal and vertical planes. The dike tip in two and three dimensions is surrounded by a tubular core of tensile (σ 1 , σ 2 ) and shear stress (τ max ). The displacement field is characterized by counterrotating cells on either side of the dike tip which, in vertical orientation, produce the characteristic subsidence above the dike complex, with uplift on either side, forming a ridge-trough-ridge structure. A finite element model of Kilauea's shield computes the displacement fields and principal stress (σ 1 ) distributions resulting from intrusive activity on each or both of the rift zones. Within the summit region, tensile stress lobes produced by the three-dimensional upward extension of the intrusions superpose constructively to produce calderawide regimes of tensile stress, conducive to caldera development. Parametric studies of (1) intrusion in the East Rift Zone only, (2) intrusion in the Southwest Rift Zone only, and (3) intrusion in both rift zones demonstrate their unique kinematic contributions. For case 1, the caldera undergoes a counterclockwise rotation (torque up state) conducive to the development of rightstepping en echelon eruptive fissures, as exemplified by the August 14, 1971, eruption. For case 2, the caldera undergoes a clockwise rotation (torque down state) conducive to the development of left-stepping eruptive fissures, as occurred during the December 31, 1974, eruption. For case 3, the caldera substructure is driven due southward, producing the southward migration of the upper portions of the summit magma reservoir.

Journal of Geophysical Research Solid Earth↗
USGS2012

The mechanics of ground deformation precursory to dome-building extrusions at Mount St. Helens 1981-1982

Detailed monitoring at Mount St. Helens since 1980 has enabled prediction of the intermittent eruptive activity (mostly dome growth) with unprecedented success. During 1981 and 1982, accelerating deformation of the crater floor around the vent (including radial cracks, thrust faults, and ground tilt) was the earliest indicator of impending activity. Numerical experiments using the finite element method to model the mechanics of the crater floor show that all deformation features can be explained by a uniform shear-stress boundary condition along the conduit wall. The magnitude of the shear stress required to match observed displacements (1–7 MPa) is inversely proportional to the conduit diameter (estimated to be 25–100 m). The most probable source of this shear stress is the flow of viscous magma up the conduit and into the lava dome. We propose a model in which the accelerating deformation, beginning as much as 4 weeks before extrusions, is caused by the increasing velocity of ascending magma in the conduit. This model is examined by using deformation data of the dome before four extrusions in 1981 and 1982 to estimate the volumetric flow rate through the conduit. This flow rate and an estimate of the effective viscosity of the magma enable calculation of an ascent velocity and an applied shear stress that, again, depend on the conduit diameter. The results of these calculations are consistent with the finite element experiments and show that the proposed model is feasible. In light of this model, events observed just before or near the time extrusions began, such as reversals of ground tilt direction from outward to inward and the sudden decrease in the number of shallow earthquakes, may indicate an abrupt decrease of shear stress in the conduit. This could be explained by a decrease in either the ascent velocity, or the effective viscosity of the magma ascending through the shallow conduit, or both, near the time of extrusion. Precursory deformation like that measured at Mount St. Helens should be observable at similar volcanoes elsewhere because it is caused by the fundamental process of magma ascent.

Journal of Geophysical Research Solid Earth↗
USGS2012

Resonance of a fluid-driven crack: Radiation properties and implications for the source of long-period events and harmonic tremor

A dynamic source model is presented, in which a three-dimensional crack containing a viscous compressible fluid is excited into resonance by an impulsive pressure transient applied over a small area ΔS of the crack surface. The crack excitation depends critically on two dimensionless parameters called the crack stiffness, C = ( b /μ)( L / d ), and viscous damping loss, F = (12η L )/(ρ ƒ d 2 α), where b is the bulk modulus, η is the viscosity, ρ ƒ is the density of the fluid, μ is the rigidity, α is the compressional velocity of the solid, L is the crack length, and d is the crack thickness. The first parameter characterizes the ability of the crack to vibrate and shapes the spectral signature of the source, and the second quantifies the effect of fluid viscosity on the duration of resonance. Resonance is sustained by a very slow wave trapped in the fluid-filled crack. This guided wave, called the crack wave, is similar to the tube wave propagating in a fluid-filled borehole; it is inversely dispersive, showing a phase velocity that decreases with increasing wavelength, and its wave speed is always lower than the acoustic velocity of the fluid, decreasing rapidly as the crack stiffness increases. The source spectrum shows many sharp peaks characterizing the individual modes of vibration of the crack; the variation of spectral shape, both in the number and width of peaks, is surprisingly complex, reflecting the interference between the lateral and longitudinal modes of resonance, as well as nodes for these modes. The far-field spectrum is marked by narrow-band dominant and subdominant peaks that reflect the interaction of the various source modes. The frequency of the dominant spectral peak radiated by the source is independent of the radiation direction. The frequency, bandwidth, and spacing of the resonant peaks are strongly dependent on the crack stiffness, larger values of the stiffness factor shifting these peaks to lower frequencies and decreasing their bandwidth. The excitation of a particular mode depends on the position of the trigger and on the extent of the crack surface affected by the pressure transient. Fluid viscosity decreases the amplitudes of the main spectral peaks, smears out the finer structure of the spectrum, and greatly reduces the duration of the radiated signal. The energy loss by radiation is stronger for high frequencies, producing a seismic signature that is marked by a high-frequency content near the onset of the signal and dominated by a longer-period component of much longer duration in the signal coda. Such signature is in harmony with those displayed by long-period events observed on active volcanoes and in hydrofracture experiments. The very low velocity which is possible in a crack with high stiffness ( C ≥ 100) also provides an attractive explanation for very long period tremor, such as type 2 tremor at Aso volcano, Japan, without the requirement of an unrealistically large magma container. The standing wave pattern set up on the crack surface by the sustained resonance in the fluid is observable in the near field of the crack, suggesting that the location and extent of the source may be estimated from the mapping of the pattern of nodes and antinodes seen in its vicinity. According to the model, the long-period event and harmonic tremor share the same source but differ in the boundary conditions for fluid flow and in the triggering mechanism setting up the resonance of the source, the former being viewed as the impulse response of the tremor generating system and the latter representing the excitation due to more complex forcing functions.

Journal of Geophysical Research Solid Earth↗
USGS2012

A magmatic model of Medicine Lake Volcano, California

Medicine Lake volcano is a Pleistocene and Holocene shield volcano of the southern Cascade Range. It is located behind the main Cascade arc in an extensional tectonic setting where high-alumina basalt is the most commonly erupted lava. This basalt is parental to the higher-silica calc-alkaline and tholeiitic lavas that make up the bulk of the shield. The presence of late Holocene, chemically identical rhyolites on opposite sides of the volcano led to hypotheses of a large shallow silicic magma chamber and of a small, deep chamber that fed rhyolites to the surface via cone sheets. Subsequent geophysical work has been unable to identify a large silicic magma body, and instead a small one has apparently been recognized. Some geologic data support the geophysical results. Tectonic control of vent alignments and the dominance of mafic eruptions both in number of events and volume throughout the history of the volcano indicate that no large silicic magma reservoir exists. Instead, a model is proposed that includes numerous dikes, sills, and small magma bodies, most of which are too small to be recognized by present geophysical methods.

Journal of Geophysical Research Solid Earth↗
USGS2012

Mineralization, alteration, and hydrothermal metamorphism of the ophiolite-hosted Turner-Albright sulfide deposit, southwestern Oregon

The Turner-Albright sulfide deposit, part of the Josephine ophiolite, formed on and below the seafloor during Late Jurassic volcanism at a back arc spreading center. Ore fluids were probably localized by faults which were active on the seafloor at the time of sulfide deposition. The uppermost massive sulfide formed on the seafloor at hydrothermal vents, similar to sulfide-depositing hot springs on modern spreading centers. The bulk of the sulfide mineralization formed below the seafloor within olivine basalt hyaloclastite erupted near the time of mineralization. Infiltration of hydrothermal fluid into the hyalo-clastite altered the rock to quartz + Fe-chlorite and quartz + sericite + Fe-chlorite. Intense alteration resulted in replacement of the protolith by quartz + sulfide. The fluid responsible for the hydrothermal alteration was evolved seawater with low p H and Mg and high Fe. Sulfide δ 534 S values average 4.7‰ and indicate contribution of sulfur from both basalt and seawater sulfate. The average value of sulfide δ 34 S and the difference between sulfide and contemporaneous seawater sulfate δ 34 S values are similar to ophiolite-hosted sulfide deposits in Cyprus. Chlorite-rich hydrothermally altered basalt has δ 18 O values of 3.0–6.8‰ due to high-temperature (250°–350°C) reaction with a seawater-derived ore fluid. Quartz in silicified basalt has higher δ 18 O values of 12.3–15.5‰, probably reflecting lower-temperature recrystallization of amorphous silica formed during initial alteration of basalt. Mudstone and clinopyroxene basalt above the sulfide horizons were not altered by the ore-transporting hydrothermal fluid, but these rocks were hydrothermally metamorphosed by altered seawater heated by deep circulation into hot oceanic crust. This subseafloor metamorphism produced a mineral assemblage typical of prehnite-pumpellyite facies metamorphism. Exchange with altered seawater increased the whole-rock δ 18 O of the basalts to values of 9.4–11.2‰.

Journal of Geophysical Research Solid Earth↗
USGS2012

Paleoseismicity and neotectonics of the Cordillera Blanca fault zone, northern Peruvian Andes

The Cordillera Bianca fault zone is a major west dipping normal fault that bounds the west side of a 120- to 170-km-wide zone of active extension along the crest of the northern Peruvian Andes. The fault is approximately 210 km long and exhibits continuous geomorphic evidence of repeated late Pleistocene and Holocene displacements but has not been the source of historical or teleseismically recorded earthquakes. Trenching and mapping of fault scarps provide new information on earthquake recurrence, slip rate, timing of the most recent events, and Andean neotectonics. At Quebrada Queroccocha, 55 km from the south end of the zone, an 11,000- to 14,000-year-old moraine is displaced 12–15 m, and younger valley fill lacustrine and fluvial deposits are displaced 7.5–8 m. Scarp profiles, tectonic terraces, and trench exposures indicate five to seven scarp-forming earthquakes of 2–3 m per event during the past 11,000–14,000 years at this location. Considering uncertainties in the number of events and in the age of the displaced moraine, the average recurrence interval is estimated to be 2440±1060 years. Detrital charcoal from colluviums that predate and postdate the most recent event yielded ages of 2480±65 and 750+80 14 C years B.P., respectively, and stratigraphic relationships suggest that the event may be closer to the older date. The late Quaternary slip rate at Quebrada Queroccocha is 0.86–1.36 mm/yr and appears to be representative of much of the fault zone. At Pachma Bajo, located 30 km from the north end of the zone on a separate segment, displacement during each of the past two events was 2 to ≥ 3 m. The morphology of scarps in alluvial fan and debris flow deposits indicates recurrence intervals of about 1000–3000 years. Relationships between pre-Inca walls and faulted deposits suggest that 1500–2000 years have elapsed since the most recent event. Fault behavior parameters, particularly displacement per event and fault segment length, show that the Cordillera Bianca fault zone is similar to seismogenic normal faults that nucleate M 7–7½ earthquakes at depths of approximately 15 km. The Quebrada Queroccocha slip rate and analysis of other normal faults in the region suggest cumulative late Quaternary extensional strain rates of about 1–4 mm/yr across the complete zone of Andean crestal extension.

Journal of Geophysical Research Solid Earth↗
USGS2012

Deformation in the Yakataga seismic gap, Southern Alaska, 1980-1986

A 60-by-40-km trilateration network in the Yakataga seismic gap was surveyed in 1980, 1982, 1984, and 1986 with precise electro-optical distance-measuring equipment to measure strain accumulation. The overall deformation is roughly approximated by a 0.24±0.03 µstrain/yr N32°W±2.4° uniaxial contraction that is uniform in time. However, the spatial distribution of deformation shows some concentration of convergence in the neighborhood of the Chugach-St. Elias fault and of right-lateral shear across the Contact fault. A simple dislocation model of the plate interaction in the Yakataga gap fits the observed deformation reasonably well but seems to require that the motion of the Pacific plate relative to the North American plate be directed more nearly N36°W than N15°W, the generally accepted direction of relative motion for this location. However, the direction of plate motion inferred from the dislocation model depends upon details of the interaction at the plate boundary that may not have been modeled accurately. A nearby but smaller trilateration network at Icy Bay was surveyed in 1982, 1984, and 1986. This network spans the southwest corner of the rupture zone of the 1979 St. Elias earthquake. The deformation at Icy Bay consists of left-lateral shear across a northeast trending zone. The relation of this deformation to strain accumulation in the Yakataga gap, postseismic relaxation associated with the 1979 earthquake, or rebound from the unloading associated with the rapid recession of the Guyot glacier is not understood.

Journal of Geophysical Research Solid Earth↗
USGS2012

Conductive heat flux in VC-1 and the thermal regime of Valles caldera, Jemez Mountains, New Mexico

Over 5% of heat in the western United States is lost through Quaternary silicic volcanic centers, including the Valles caldera in north central New Mexico. These centers are the sites of major hydrothermal activity and upper crustal metamorphism, metasomatism, and mineralization, producing associated geothermal resources. We present new heat flow data from Valles caldera core hole 1 (VC-1), drilled in the southwestern margin of the Valles caldera. Thermal conductivities were measured on 55 segments of core from VC-1, waxed and wrapped to preserve fluids. These values were combined with temperature gradient data to calculate heat flow. Above 335 m, which is probably unsaturated, heat flow is 247±16 mW m −2 . The only deep temperature information available is from an uncalibrated commercial log made 19 months after drilling. Gradients, derived from uncalibrated temperature logs, and conductivities are inversely correlated between 335 and 737 m, indicating a conductive thermal regime, and component heat fluxes over three depth intervals (335–539 m, 549–628 m, and 628–737 m) are in excellent agreement with each other with an average of 504±15 mW m −2 . Temperature logs to 518 m depth with well-calibrated temperature sensors result in a revised heat flow of 463±15 mW m. We use shallow thermal gradient data from 75 other sites in and around the caldera to interpret the thermal regime at the VC-1 site. A critical review of published thermal conductivity data from the Valles caldera yields an average thermal conductivity of ≥1 W m −1 K −1 for the near-surface tuffaceous material, and we assume that shallow gradient values (°C km −1 ) are approximately numerically equal to heat flow (mW m −2 ). Heat loss from the caldera is asymmetrically distributed, with higher values (400 mW m −2 or higher) concentrated in the west-southwestern quadrant of the caldera. This quadrant also contains the main drainage from the caldera and the youngest volcanism associated with the caldera. We interpret the shallow thermal gradient data and the thermal regime at VC-1 to indicate a long-lived hydrothermal (and magmatic) system in the southwestern Valles caldera that has been maintained through the generation of shallow magma bodies during the long postcollapse history of the caldera. High heat flow at the VC-1 site is interpreted to result from hot water circulating below the base of the core hole, and we attribute the lower heat flow in the unsaturated zone to hydrologic recharge.

Journal of Geophysical Research Solid Earth↗
USGS2012

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

Geothermal gradients in the conterminous United States

Geothermal gradients from published temperature/depth measurements in drill holes generally deeper than 600 m are used to construct a temperature gradient map of the conterminous United States. The broadly contoured map displays 284 temperature gradients that are applicable to a depth of 2 km. In terms of the number of contoured areas and the fraction of data points having a value not within a contour interval (outliers), the temperature gradient data and associated deep heat flow data have similar measures of contourability. Areally, most of the United States is contoured from 15° to 35°C/km. The eastern United States is generally cooler (average 25°C/km) than the western United States (average 34°C/km), in accordance with broad heat flow trends. Differences between the temperature gradient and heat flow maps are caused by areal differences in rock thermal conductivities. The effect of conductivity on gradients is particularly apparent in the eastern United States where heat flow is relatively constant over large areas. Gradients are elevated where thick, low-conductivity, sedimentary deposits occur such as in the Atlantic Coastal Plain province and in basins in the Allegheny Plateau and the Great Plains provinces. No clear gradient pattern emerges where both heat flow and conductivity vary widely, such as in the northern Basin and Range and Rocky Mountain provinces. Using the temperature gradients determined in this study and associated heat flow values, derived thermal conductivities are calculated for the depth range of a few hundred meters to 2 km; the average conductivity is 2.5±0.8 W/m °C. Some areas show little variation in derived thermal conductivity, while others show a wide range. The Atlantic Coastal Plain, Appalachian Plateaus, Superior Upland, Gulf Coastal Plain, Columbia Plateaus and Cascade Mountains, and Sierra Nevada provinces have restricted ranges of conductivities, which supports the concept of regional conductivities.

Journal of Geophysical Research Solid Earth↗
USGS2012

Distribution of late Cenozoic volcanic vents in the Cascade Range: Volcanic arc segmentation and regional tectonic considerations

Spatial, temporal, and compositional distributions of approximately 4000 volcanic vents formed since 16 Ma in Washington, Oregon, northern California, and northwestern Nevada illustrate the evolution of volcanism related to subduction of the Juan de Fuca plate system and extension of the Basin and Range province. Vent data were obtained from published map compilations and include monogenetic and small polygenetic volcanoes in addition to major composite centers. On the basis of the distribution of 2821 vents formed since 5 Ma, the Cascade Range is divided into five segments, with vents of the High Lava Plains along the northern margin of the Basin and Range province in Oregon forming a sixth segment. Some aspects of the Cascade Range segmentation can be related to gross structural features of the subducting Juan de Fuca plate. The orientation of the volcanic front of segments one and two changes from NW in northern Washington to NE in southern Washington, paralleling the strike of the subducting Juan de Fuca plate. Segments one and two are separated by a 90-km volcanic gap between Mount Rainier and Glacier Peak that is landward of the portion of the subducting plate having the least average dip to a depth of 60 km. A narrow, N-S trending belt of predominantly andesitic vents in Oregon constitutes a third segment, which is landward of the seismically quiet portion of the subduction zone. The narrowness of this segment may indicate steep dip of the subducting plate beneath the Cascade arc in Oregon. Vents are sparse between segment four (containing the Mount Shasta and Medicine Lake centers) and segment five (containing Lassen Peak), where the Juan de Fuca and Gorda North plates are characterized by differing age, amounts of subcrustal seismicity, and probably geometry. From the relation between seismicity at depth of 60 km and the position of the volcanic front of vents formed since 5 Ma, transitions between subducting-plate segments of varying geometry likely occur near boundaries between independently defined volcanic segments in northern Oregon and northern California. In the Basin and Range province east of the Cascade arc, volcanism migrated into the region adjacent to the Cascade Range during the interval 5–10 Ma. Since 5 Ma, the impingement of the two provinces is characterized by cessation of basin-range volcanism in southern Oregon, continuation of basaltic volcanism in northeastern California where the impingement process may not yet be complete, and contraction of the area of mafic volcanism around Mount Shasta, Medicine Lake, and Lassen Peak. In central Oregon where the northern margin of basin-range volcanism (the High Lava Plains) intersects the Cascade arc, impingement of basin-range extensional volcanism approximately coincides in time and space with the development of the High Cascade graben between Three Sisters and Mount Jefferson.

Journal of Geophysical Research Solid Earth↗