USGS ScienceSearch

SEARCH · USGS Science

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 181 records · Page 10Linked to original sources

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

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

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

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

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

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

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

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

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

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

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

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

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

Absence of strain accumulation in the Shumagin seismic gap, Alaska, 1980-1987

Measurements of the deformation of a trilateration network in the Shumagin seismic gap in the interval 1980–1987 failed to detect any significant strain accumulation (observed extension rate in the direction of plate convergence 0.00±0.03 μstrain/yr). Dislocation models of the subduction process and measurements at a comparable network at a known seismic subduction zone (Nankai Trough, Japan) suggest that a rate of the order of −0.2 μstrain/yr should have been observed if the main thrust zone beneath the Shumagin Islands were locked. The simplest explanation of the observed absence of strain accumulation in the Shumagin seismic gap is that the main thrust zone beneath the Shumagin Islands is not presently locked. Other possible explanations depend upon very particular circumstances.

Journal of Geophysical Research Solid Earth

Three-dimensional P velocity structure of the crust below Newberry Volcano, Oregon

Utilizing teleseismic P residuals, we have detected a column of high P velocity material extending from within 10 km of the surface below the summit of Newberry Volcano, Oregon, to midcrustal depths near 25 km. We interpret this column to be the expression of a swarm of predominantly subsolidus gabbroic sills and dikes which were intruded as the volcano was built. The high P velocities observed below the volcano severely limit the size of magma chambers which could presently exist in the crust below Newberry Volcano. Those possible include a few percent of partial melt distributed through large volumes of a mafic intrusion zone in the midcrust; a few, smaller, higher melt fraction zones in the midcrust with dimensions less than 6 km and whose aggregate volume is only a few percent of enclosing volumes of 200 km 3 ; small magma bodies with dimensions of a few kilometers located within the upper 10 km of the crust; or a mafic, crystal-rich magma of arbitrary dimensions located in the upper few km. The low P velocities detected in the upper 4 km beneath the center of the summit caldera may be partially caused by a magma chamber in the second of these catagories.

Journal of Geophysical Research Solid Earth

Electrical structure of Newberry Volcano, Oregon

From the interpretation of magnetotelluric, transient electromagnetic, and Schlumberger resistivity soundings, the electrical structure of Newberry Volcano in central Oregon is found to consist of four units. From the surface downward, the geoelectrical units are (1) very resistive, young, unaltered volcanic rock, (2) a conductive layer of older volcanic material composed of altered tuffs, (3) a thick resistive layer thought to be in part intrusive rocks, and (4) a lower-crustal conductor. This model is similar to the regional geoelectrical structure found throughout the Cascade Range. Inside the caldera, the conductive second layer corresponds to the steep temperature gradient and alteration minerals observed in the USGS Newberry 2 test hole. Drill hole information on the south and north flanks of the volcano (test holes GEO N-1 and GEO N-3, respectively) indicates that outside the caldera the conductor is due to alteration minerals (primarily smectite) and not high-temperature pore fluids. On the flanks of Newberry the conductor is generally deeper than inside the caldera, and it deepens with distance from the summit. A notable exception to this pattern is seen just west of the caldera rim, where the conductive zone is shallower than at other flank locations. The volcano sits atop a rise in the resistive layer, interpreted to be due to intrusive rocks. The intrusive material has served as a heat source to produce enhanced hydrothermal alteration and, perhaps in the case of the west-flank anomaly, elevated fluid temperatures. While no public drill hole information is available to confirm this hypothesis, the west-flank anomaly appears to be a good geothermal target. In addition to the possibility that a region on the west side of the volcano could be favorable for prospecting, part of the resistive structure under the center of the volcano could be due to a vapor-dominated environment with temperatures above 300°C. In other parts of the Cascades, pervasive alteration has produced mixed layer clays and zeolites, resulting in low-resistivity anomalies. Low resistivities cannot be assumed to indicate high-temperature pore fluids. The use of electrical methods that measure resistivity as a function of excitation frequency, such as spectral induced polarization, may provide a way of obtaining information about the type and extent of alteration.

Journal of Geophysical Research Solid Earth

The hydrothermal system at Newberry Volcano, Oregon

Results of recent geological and geophysical studies at Newberry Volcano have been incorporated into conceptual and numerical models of a magma-based hydrothermal system. Numerical simulations begin with emplacement of a small magma body, the presumed source of silicic eruptions at Newberry that began about 10,000 B.P., into a thermal regime representing 100,000 years of cooling of a large underlying intrusion. Simulated flow patterns and thermal histories for three sets of hypothetical permeability values are compatible with data from four geothermal drill holes on the volcano. Meteoric recharge cools the caldera-fill deposits, but thermal water moving up a central conduit representing a permeable volcanic vent produces temperatures close to those observed in drill holes within the caldera. Meteoric recharge from the caldera moves down the flanks and creates a near-isothermal zone that extends several hundred meters below the water table, producing temperature profiles similar to those observed in drill holes on the flanks. The temperatures observed in drill holes on the flanks are not influenced by the postulated Holocene magma body. The elevated temperature gradients measured in the lower portions of these holes may be related to the cumulative effect of older intrusions. The models also indicate that meteoric recharge to the deep hydrothermal system probably originates within or near the caldera. Relatively low fluid velocities at depth suggest that at least a significant fraction of the thermal fluid may be very old.

Journal of Geophysical Research Solid Earth

Paleomagnetic investigation of some volcanic rocks from the McMurdo volcanic province, Antarctica

Paleomagnetic data for lava flows from sporadic but long-lived eruptions in the McMurdo Sound region are combined with previously published geologic and geochronologic data to determine the general eruptive sequence of the area. Lava flows in the Walcott Bay area were erupted during the Gauss Normal, Matuyama Reversed, and Brunhes Normal Polarity Chrons. The youngest flows on Black Island probably erupted near the boundary between the Gilbert and Gauss chrons. The most recent activity was concentrated on the volcanic edifices of Mounts Morning and Discovery and on Ross Island where Mount Erebus is currently active. Combining virtual geomagnetic poles calculated for 39 units sampled during this study with those of eight flows that were published previously yields a mean paleomagnetic pole at 87.3°N, 317.3°E (α 95 = 6.3°). The ancient geomagnetic field dispersion about this mean pole is 23.5°, with upper and lower limits of 95% confidence equal to 27.4° and 20.5°, respectively. This value probably is a reasonable estimate of secular variation for the Antarctic continent during Pliocene and Pleistocene time.

Journal of Geophysical Research Solid Earth