USGS Science⌕ Search

SEARCH · USGS Science

Results for “The R Journal”

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

Fractal properties of tremor and gas piston events observed at Kilauea Volcano, Hawaii

Studies the fractal properties of shallow volcanic tremor and gas piston events associated with magma degassing at Kilauea Volcano, Hawaii, using data from two dense short-baseline arrays of seismographs deployed near the active crater of Puu Oo on the east rift of the volcano. The existence of a categorically stable attractor characterizing both types of activities strongly suggests that the excitation mechanism of tremor is similar to that of gas piston events, which in turn are correlated with visual observations at the volcanic vent. Low values derived for the dimension of the attractor in phase space points to significant self-organization in the process of generation of tremor and offers general constraints on the dimensionality of attractors derived from models of acoustic emission associated with magma flow, vesiculation, and degassing. -from Authors

Hawaii↗
USGS2012

Crustal subsidence and extension and Medicine Lake volcano, northern California

The pattern of historical ground deformation, seismicity, and crustal structure near Medicine Lake volcano illustrates a close relation between magmatism and tectonism near the margin of the Cascade volcanic chain and the Basin and Range tectonic province. Between leveling surveys in 1954 and 1989 the summit of Medicine Lake volcano subsided 389±43 mm with respect to a reference bench mark 40 km to the southwest (average rate = 11.1±1.2 mm/yr). A smaller survey across the summit caldera in 1988 suggests that the subsidence rate was 15–28 mm/yr during 1988–1989. Swarms of shallow earthquakes ( M ≤ 4.6) occurred in the region during August 1978, January–February 1981, and September 1988. Except for the 1988 swarm, which occurred beneath Medicine Lake caldera, most historical earthquakes were located at least 25 km from the summit. The spatial relation between subsidence and seismicity indicates (1) radially symmetric downwarping of the volcano's summit and flanks centered near the caldera and (2) downfaulting of the entire edifice along regional faults located 25–30 km from the summit. We propose that contemporary subsidence, seismicity, and faulting are caused by (1) loading of the crust by more than 600 km 3 of erupted products plus a large volume of mafic intrusives; (2) east‐west extension in the western Basin and Range province; and, to a lesser extent, (3) crystallization or withdrawal of magma beneath the volcano. Thermal weakening of the subvolcanic crust by mafic intrusions facilitates subsidence and influences the distribution of earthquakes. Subsidence occurs mainly by aseismic creep within 25 km of the summit, where the crust has been heated and weakened by intrusions, and by normal faulting during episodic earthquake swarms in surrounding, cooler terrain.

California↗
USGS2012

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

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

Great Lakes↗
USGS2012

Interseismic deformation at the Nankai Trough, Japan, subduction zone

Uplift along the coast of southwestern Japan following the 1944 Tonankai ( M s = 8.0) and 1946 Nankaido ( M s = 8.2) earthquakes has been inferred from the 1950–1985 tide gage records. Although uplift rates vary with geographic position, the temporal dependence at each site can be described as the superposition of an exponentially decaying (time constant ∼5 years) transient, significant only in the first decade postseismic, and a steady interseismic rate. The steady interseismic uplift rate is defined only over the 1956–1985 interval, but the linearity in time ofthat uplift record is such that uplift over the remainder of the interseismic cycle (∼80 years) might reasonably be expected to occur at the same rate. Moreover, the spatial pattern and magnitude of the steady uplift rate are reasonably close to those predicted by a simple dislocation model of subduction. However, the currently observed steady uplift rate is not consistent with the uplift rate inferred from geodetic leveling for the final half of the previous earthquake cycle. Thus, either the hypothesis of a steady interseismic uplift rate after the first decade postseismic is wrong or the strain accumulation pattern for one earthquake need not be the same as for the preceding earthquake on the same segment.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Paleomagnetism and geochronology of 23 Ma gabbroic intrusions in the Keku Strait, Alaska, and implications for the Alexander Terrane

Samples of Tertiary gabbro from 24 sites in the Keku Strait, Alaska, help constrain the displacement history of the Alexander terrane. Step heating experiments on a plagioclase separate from these previously undated intrusions indicate a discordant 40 Ar/ 39 Ar age of 23.1 ± 1.7 Ma. The characteristic magnetization resides in magnetite, is easily isolated by thermal and alternating field demagnetization, and has both normal and reversed polarities. The mean paleomagnetic pole, with no structural correction, is latitude 87.1°N, longitude 141.6°E, A 95 = 10.1°, with N = 20 sites. This pole indicates insignificant tectonic displacement (0.5° ± 8.2° southward) and rotation (0.6° ± 15.2° counterclockwise). We therefore conclude that any northward displacement or vertical axis rotation of the Alexander terrane, and/or tilting in the vicinity of the Keku Strait must have occurred before 23 Ma.

Alaska↗
USGS2012

Seasonal fluctuations in sea level on the South Carolina shelf and their relationship to the Gulf Stream

Near-bottom pressure measurements obtained over a 9-month period in 1978 on the outer continental shelf off South Carolina document that sea level rose 35 to 50 cm between July and October. Records of coastal sea level showed a similar rise. When the pressure and coastal sea level records were corrected for the effects of wind stress and temperature, the largest portion of the rise in the two variables remained. The amplitude of the change in residual sea level was largest at the shelf break, 28 cm, and decayed to 22 cm at the coast. The relatively small magnitude of the decay suggests that alongshelf pressure gradients imposed at the shelf break have minimum alongshelf wavelengths of 2000 km. The changes in residual sea level were seasonal and associated with seasonal changes in the transport of the Gulf Stream. The 3 years of residual sea level records examined here indicate that the amplitude of a typical decrease in transport of the Gulf Stream off South Carolina between July and October has a strong interannual variability. The change in transport can be 50% lower to over 100% higher than the historical average. The residual sea level records also suggest that a change in transport of the Gulf Stream is not strongly related to the onshore/offshore position of the Gulf Stream.

South Carolina↗
USGS2012

Evolution of the caldera‐forming eruption at Crater Lake, Oregon, indicated by component analysis of lithic fragments

Crater Lake caldera (8 × 10 km), formed 6845 years B. P. ( 14 C age) during the climactic eruption of the volcanic edifice known as Mount Mazama, is intermediate in size between small calderas associated with central vent eruptions and large calderas that have ring fracture vent systems. Our quantitative study of lithic fragments in the ejecta confirms the existing model of changes in vent configuration during the climactic eruption of Mount Mazama. Initial activity was from a single vent that produced a rhyodacite pumice fall from a Plinian column. Altered preexisting volcanic rocks are the predominant lithic type in the Plinian deposit, and their extensive hydrothermal alteration is considered as evidence of their relatively deep origin. The Wineglass Welded Tuff lies atop the Plinian deposit and contains a higher proportion of fresh volcanic rocks, suggesting enlargement of the single vent by slumping of its walls. This same vent enlargement caused the Plinian eruption column to collapse and feed valley‐hugging pyroclastic flows that deposited the Wineglass Welded Tuff. When enough material was erupted from the shallow magma chamber that its roof was no longer adequately supported, Mount Mazama collapsed to form the caldera, while highly energetic pyroclastic flows produced the climactic ignimbrite. A lag breccia that represents the proximal facies of the compositionally zoned climactic ignimbrite lies atop the Wineglass Welded Tuff and contains predominantly altered volcanic rocks of deeper origin, accompanied by minor granitoids from the magma chamber walls. Azimuthal differences in lithic component proportions in the lag breccia correlate well with the geology of the caldera walls, indicating that the climactic ignimbrite was ejected by multiple vents along a ring fracture system. Systematic lithic component changes within the lag breccia suggest different quarrying levels that reflect waxing and waning of the discharge rate during the caldera collapse phase of the climactic eruption. Our lithic component analysis demonstrates that calderas that may be too small to experience structural resurgence, such as Crater Lake, nevertheless may form by syneruptive subsidence along ring fractures.

Oregon↗
USGS1996

Cyclic stressing and seismicity at strongly coupled subduction zones

We use the finite element method to analyze stress variations in and near a strongly coupled subduction zone during an earthquake cycle. Deformation is assumed to be uniform along strike (plane strain on a cross section normal to the trench axis), and periodic earthquake slip is imposed consistent with the long-term rate of plate convergence and degree of coupling. Simulations of stress and displacement rate fields represent periodic fluctuations in time superimposed on an average field. The oceanic plate, descending slab, and continental lithosphere are assumed here to respond elastically to these fluctuations, and the remaining mantle under and between plates is assumed to respond as Maxwell viscoelastic. In the first part of the analysis we find that computed stress fluctuations in space and time are generally consistent with observed earthquake mechanism variations with time since a great thrust event. In particular, trench-normal extensional earthquakes tend to occur early in the earthquake cycle toward the outer rise but occur more abundantly late in the cycle in the subducting slab downdip of the main thrust zone. Compressional earthquakes, when they occur at all, have the opposite pattern. Our results suggest also that the actual timing of extensional outer rise events is controlled by the rheology of the shallow aseismic portion of the thrust interface. The second part of the analysis shows the effects of mantle relaxation on the rate of ground surface deformation during the earthquake cycle. Models without relaxation predict a strong overall compressional strain rate in the continental plate above the main thrust zone, with the strain rate constant between mainshocks. However with significant relaxation present, a localized region of unusually low compressional, or even slightly extensional, strain rate develops along the surface of the continental plate above and somewhat inland from the downdip edge of the locked main thrust zone. The low strain rate starts in the middle or late part of the cycle, depending on position. This result suggests that the negligible or small contraction measured on the Shumagin Islands, Alaska, during 1980 to 1991, may not invalidate an interpretation of that region as being a moderately coupled subduction zone. In contrast, mantle relaxation causes only modest temporal nonuniformity of uplift rates in the overriding plate and of extensional stress rates in the subducting plate, even when the Maxwell time is an order of magnitude less than the recurrence interval.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Modeling of fluidized ejecta emplacement over digital topography on Venus

The FLOW computer model of McEwen and Malin (1989) modified for application to the study of Venus fluidized ejecta blankets (FEBs) demonstrates that relatively low viscosities, yield strengths, and initial velocities are required to duplicate the observed flow paths of the outflow materials. The model calculates the velocities and simulated flow paths of gravity flows over Magellan topography. The model is formulated to determine flow movements from initial conditions, gravitational acceleration, and resistance to motion as described by Coulomb, viscous, and turbulent resistance forces. Successful duplication of observed FEB flow paths has been achieved for the FEB craters Addams, Isabella, and Cochran. When used as a simple energy‐line model, the model requires low coefficients of friction to extend FEBs to near their observed termini in the synthetic aperture radar (SAR) imagery, although the resulting straight flow lines do not follow the observed flow paths well. For Bingham flow, the model requires low values of viscosity and yield strength which are more similar to pyroclastic or debris flows than basaltic lavas. Flows of 100‐m depth require 1 to 2 orders of magnitude higher values of both viscosity and yield strength than 10‐m‐deep flows. The complicated nature of the flow lines for the low velocity model suggests that FEBs were probably emplaced under variably laminar and turbulent flow conditions, where underlying topography influenced both the direction and energy of flow materials.

Journal of Geophysical Research E: Planets↗
USGS1997

Land cover mapping, fire regeneration, and scaling studies in the Canadian boreal forest with 1 km AVHRR and Landsat TM data

A multitemporal 1 km advanced very high resolution radiometer (AVHRR) land cover analysis approach was used as the basis for regional land cover mapping, fire disturbance-regeneration, and multiresolution land cover scaling studies in the boreal forest ecosystem of central Canada. The land cover classification was developed by using regional field observations from ground and low-level aircraft transits to analyze spectral-temporal clusters that were derived from an unsupervised cluster analysis of monthly normalized difference vegetation index (NDVI) image composites (April-September 1992). Quantitative areal proportions of the major boreal forest components were determined for a 821 km ?? 619 km region, ranging from the southern grasslands-boreal forest ecotone to the northern boreal transitional forest. The boreal wetlands (mostly lowland black spruce, tamarack, mosses, fens, and bogs) occupied approximately 33% of the region, while lakes accounted for another 13%. Upland mixed coniferous-deciduous forests represented 23% of the ecosystem. A SW-NE productivity gradient across the region is manifested by three levels of tree stand density for both the boreal wetland conifer and the mixed forest classes, which are generally aligned with isopleths of regional growing degree days. Approximately 30% of the region was directly affected by fire disturbance within the preceding 30-35 years, especially in the Canadian Shield Zone where large fire-regeneration patterns contribute to the heterogeneous boreal landscape. Intercomparisons with land cover classifications derived from 30-m Landsat Thematic Mapper (TM) data provided important insights into the relative accuracy of the 1 km AVHRR land cover classification. Primarily due to the multitemporal NDVI image compositing process, the 1 km AVHRR land cover classes have an effective spatial resolution in the 3-4 km range; therefore fens, bogs, small water bodies, and small patches of dry jack pine cannot be resolved within the wet conifer mosaic. Major differences in the 1-km AVHRR and 30-m Landsat TM-derived land cover classes are most likely due to differences in the spatial resolution of the data sets. In general, the 1 km AVHRR land cover classes are vegetation mosaics consisting of mixed combinations of the Landsat classes. Detailed mapping of the global boreal forest with this approach will benefit from algorithms for cloud screening and to atmospherically correct reflectance data for both aerosol and water vapor effects. We believe that this 1 km AVHRR land cover analysis provides new and useful information for regional water, energy, carbon, and trace gases studies in BOREAS, especially given the significant spatial variability in land cover type and associated biophysical land cover parameters (e.g., albedo, leaf area index, FPAR, and surface roughness). Multiresolution land cover comparisons (30 m, 1 km, and 100 km grid cells) also illustrated how heterogeneous landscape patterns are represented in land cover maps with differing spatial scales and provided insights on the requirements and challenges for parameterizing landscape heterogeneity as part of land surface process research.

Journal of Geophysical Research D: Atmospheres↗
USGS1998

The coseismic slip distributions of the 1940 and 1979 Imperial Valley, California, earthquakes and their implications

Geodetic arrays observed by the U.S. Coast and Geodetic Survey span the Imperial fault in southern California. For the 1940 M 7.1 Imperial Valley earthquake, a 1934–1941 triangulation network has sufficient resolution to allow inversion for the coseismic slip distribution on fault segments 5 to 25 km long extending from the surface to a depth of 9 km. The estimated right-lateral slip is 0.8 to 1.7 m on the northern 30 km of the main trace of the Imperial fault, 4.8±0.2 m on a 10-km-long segment straddling the United States - Mexico border, and 1.3±0.4 m on a southern 25-km-long segment in Mexico. Fixing this strike-slip model and inverting 1940 leveling data only for dip slip yields 0.1 m of east-side-down dip slip. The seismic moment for this model is M 0 = (3.2±0.3) ×10 19 N m. The 1979 geodetic data set, mostly elevation changes from leveling routes, has insufficient resolution for inversion. However, it is possible to use this geodetic data set and results published by others to infer that the 1940 and 1979 earthquakes may be similar on the rupture zone common to both events. Our preferred 1940 model is similar to the 1979 geodetic results of Crook [1984] on the segments where both networks have good resolution. Elevation changes from 1940 and 1979 leveling data are very similar. Thus the geodetic data corroborate the surface slip evidence of Sharp [1982b] that the 1940 and 1979 slip distributions are examples of “characteristic slip” on the northern Imperial fault.

California↗
USGS2012

Elastic moduli of rock glasses under pressure to 8 kilobars and geophysical implications

Shear and longitudinal velocities were measured by the ultrasonic phase comparison method as a function of pressure to 8 kbar on synthetic glasses of basalt, andesite, rhyolite, and quartz composition and on natural obsidian. Velocities of most of the glasses decrease anomalously with pressure, but increasingly more-normal behavior occurs with decrease in SiO 2 content. The pressure derivatives of rigidity and bulk modulus increase linearly, from −3.39 to −0.26 and from −5.91 to +2.09, respectively, with decrease in SiO 2 content from 100 to 49%. The change from negative to positive in the pressure derivatives of both moduli and observed at Poisson's ratio of about 0.25 is consistent with the Smyth model for the anomalous elastic behavior of glass. If the temperature in the upper mantle is about 1500°C, tholeiitic basalt would be molten in accordance with the partial melt explanation for the low-velocity zone; at 1300°C and below, basalt would be in the glassy state, especially if more felsic than tholeiite. At a temperature of 1370°C and at 30 kbar, reasonable values for the upper mantle at 100 km depth, the basalt glass of this study would have a viscosity of about 10 13 P. On the basis of the theory of viscoelasticity the glass would support shear wave propagation at frequencies above 0.1 Hz. Under such conditions of PT, 10 to 30% basalt glass in a matrix either of eclogite or olivine would reduce the seismic velocities by 3 to 9% and could also account for the values observed in the low-velocity zone.

Journal of Geophysical Research Solid Earth↗
USGS2012

Terranes and suture zones in east central Alaska

East central Alaska, with its 17 terranes, forms a part of the mosaic of allochthonous terranes that join the North American and Siberian plates. These terranes range from continental and continental margin, such as the Tatonduk with its thick well-bedded succession of marine shelf rocks, to seamount, arc, and ocean floor terranes. The Yukon crystalline terrane, the largest described here, is a composite of at least four subterranes juxtaposed across the Tintina fault with the Tatonduk terrane, a northwestern extension of the North American plate in Alaska. Inboard of the Yukon crystalline terrane are packets of closely appressed microterranes separated from the Tatonduk and other terranes belonging to North America by major suture zones. These microterranes lie between North America and the mosaic of accretionary terranes that form the more southerly part of Alaska. The most obviously allochthonous microterranes within the suture zones are the Woodchopper Canyon, an Early Devonian basaltic seamount, and the White Mountains, an Ordovician volcanic arc terrane capped by Silurian and Devonian carbonate bank deposits. The nearest counterpart of these terranes is the Alexander terrane in southeastern Alaska. The Tintina fault of Mesozoic and Cenozoic age, like the Denali fault, primarily follows old suture zones that separate terranes. Strike slip faulting developed after collision in places where further convergence was oblique to the terrane margins. Where terranes met head-on, their leading edges lie along a multiple set of high-angle faults that outline microterranes in accretion zones.

Alaska↗
USGS2012

Recent movement on the Garlock Fault as suggested by water level fluctuations in a well in Fremont Valley, California

Water levels have been continuously recorded since March 1978 in a well in Fremont Valley, where several strands of the adjacent Garlock fault zone have exhibited both left-lateral displacement and components of normal displacement. Differences in water levels indicate that a fault segment lies between the observation well and a nearby irrigation well. During the 4-year recording period, six sharp fluctuations, or “spikes,” were noted. These fluctuations, occurring over 2- to 4-day periods, have amplitudes of 15–30 cm. They appear to be the result of creep events on a nearby fault. Two types of creep events are plausible: (1) normal slip on an en echelon trace of the Garlock fault less than 300 m south of the well, with the north side up relative to Fremont Valley or (2) left-lateral slip on the same fault. Because of the nature of the fluctuations we favor the latter interpretation. Dislocation models utilizing exponential, arc tangent, and skewed cosine functions were used to analyze the water level fluctuations, associated pressure distribution, and fault displacements. The results suggest that creep on the fault ranges from several millimeters to a centimeter for individual events. Estimates of cumulative creep for the period 1978–1982 range from 20 to 50 mm, depending on the particular model employed.

Journal of Geophysical Research Solid Earth↗
USGS2012

Implications of silicic vent patterns for the presence of large crustal magma chambers

On the basis of the distribution of silicic vents, many volcanic fields can be grouped with (1) igneous systems that may be small and whose vent locations are controlled by regional tectonics, (2) those that include sizable crustal magma bodies which erupt at sites determined by their anomalous local stress fields, or (3) relatively small volume systems that are transitional between categories 1 and 2. Linear vent patterns that are aligned normal to the regional least principal stress (σ 3 ) commonly are associated with absence of evidence for large, shallow magma bodies. The Coso volcanic field and the Inyo‐Mono domes in California and probably the South Sister area in Oregon are examples of such systems. The 1960 dacitic fissure eruption at Cordón Caulle in southern Chile evidently is linked to tectonic stress relaxation associated with the great earthquake that occurred 48 hours before the eruption began. Large, shallow magma chambers are thought to perturb the local stress field so that areal patterns of silicic vents are diffuse, radial, or arcuate. Such systems may erupt great volumes of pyroclastic material catastrophically and produce large calderas. Well‐preserved examples of late precaldera leaks of silicic magma occur at Long Valley, California, and Mount Mazama (Crater Lake), Oregon. Other possible examples are noted. Some of these, those which formed calderas smaller than Crater Lake, apparently were preceded by silicic eruptions from aligned vents. Bearing in mind that there exists a transitional group between tectonically controlled small systems and very large magma chambers, vent distributions can be useful in evaluating potential volcanic hazards for silicic volcanic fields.

California↗
USGS2012

Links between magma-tectonic rate balances, plutonism, and volcanism

With Smith's (1979) model of ash flow magmatism as the point of departure, quantitative regimes of crustal magmatic evolution are outlined. The average extrusive regime of ash flow magmatism operates at a rate of about 0.001 km 3 /yr; progressively larger episodes of increasing repose times characterize the evolution of caldera-forming eruptions. The intrusive regime operates at an inferred average rate of about 0.01 km 3 /yr expressed in terms of the volumetric influx of mafic magma from the mantle, based on the record of silicic rates of extrusion and a 10:1 intrusive to extrusive ratio derived from evidence of zoned silicic systems. This rate is almost identical with the mean cumulative rate for the Hawaiian-Emperor system of basaltic volcanism, which is determined to be 0.015 km 3 /yr. Fluctuations about a mean rate of 10 −2 km 3 /yr lead to a threefold classification of magmatic regimes in terms of source power, depending on whether they are persistently greater than (HHER regime), about equal to (MHER regime), or lower than (LHER regime) the mean Hawaiian-Emperor rate. Ash flow magmatism correlates with the MHER regime, averaged over a sufficiently long time. Thus the conditions of Hawaiian-Emperor magmatism are identified with a characteristic power of magma transfer rates that optimize mechanisms for creation of high-level chambers of derivative and fractionated silicic magma (terminology of Smith (1979)) within sufficiently localized regions of high residence times. These conditions represent a balance point, and a bifurcation, in dynamic paths of magmatic evolution leading to widely divergent types of plutonic-volcanic associations. Categorical examples of the principal types of igneous path are discussed relative to continental basalt, silicic ash flow, batholithic granite, and intrusive porphyry systems. The Henry Mountains, Utah, are thought to represent an example of an intrusive porphyry generated in the LHER regime of a mafic source system lacking sufficient power for volcanic expression or significant mineralization. Long Valley, California (and possibly other systems along the eastern front of the Sierra Nevada), is an example operating at the mean rate with transient excursions to higher and lower rates. In regions of rapid extension and/or volcanic propagation (whether linear or areal) these regimes give rise to basaltic eruptions ranging from isolated cinder cones to flood basalts.

Journal of Geophysical Research Solid Earth↗
USGS2012

Magmatic inclusions in silicic and intermediate volcanic rocks

Fine‐grained ellipsoidal inclusions from a few millimeters to over l m in size are present in many intermediate to silicic lava flows and domes. Only recently has it become widely accepted that such inclusions are chilled blobs of magma. Their magmatic origin is manifested by vesicularity and high groundmass porosity, by ellipsoidal shapes, by mingling at contacts with the host, and by textural evidence for groundmass crystallization in an undercooled state. Groundmass textures reflect the degree of undercooling, which is a function of thermal contrast (Δ T ) between inclusions and host magmas before mingling. Compositions of inclusions are invariably more mafic than their hosts and indicate that the magmatic systems were compositionally zoned. Many inclusions were formed from hybrid magmas. Four examples are given in order to document effects of decreasing Δ C (compositional contrast) and Δ T , as hosts for mafic andesitic to andesitic inclusions vary from high‐silica rhyolite (HSR) to silicic andesite. In the HSR example from the Coso volcanic field, California, inclusions track the evolution of zonation in a long‐lived magmatic system. At the other extreme, andesitic inclusions in silicic andesite from Crater Lake, Oregon, demonstrate zonation of a system of low Δ T and Δ C and show effects of in situ differentiation by gas filter‐pressing of residual liquid. Inclusions form when a relatively small proportion of mafic magma comes into contact with a larger proportion of silicic magma. Several mechanisms by which the magmas mix have been proposed: buoyant rise of vesiculated mafic magma or convective stirring at the interface between mafic and overlying silicic magma, forcible injection of mafic magma, diking of fractured silicic magma, and mixing during ascent in a conduit. No single mechanism explains all field observations. Many inclusions in plutons are analogous to magmatic inclusions in volcanic rocks but are likely to have suffered greater chemical and textural change during comparatively slow postentrapment cooling. Magmatic inclusions are evidence of open system behavior of the differentiated parts of magma bodies, and the inclusions themselves may behave as open systems where Δ T is low.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Morphometric variability within the axial zone of the southern Juan de Fuca Ridge: Interpretation from Sea MARC II, Sea MARC I, and deep-sea photography

The morphometric characteristics of the axial regions of oceanic spreading centers are determined by (1) the type of volcanic flows, (2) the relation between primary volcanic relief (on a scale of a few meters to tens of meters) and degree of sediment cover, and (3) the extent of surficial expression and timing of tectonic disruption of the young oceanic crust. Even within a single, continuous, linear spreading-ridge segment with relatively uniform axial valley dimensions over a distance of 50 or more kilometers, such as along the southern Juan de Fuca Ridge, the changes in morphometric characteristics along axis within the youngest crust indicate distinct variation in tectonic and volcanic activity over short distances within short time periods. An integrated analysis of Sea MARC I, Sea MARC II, and photographic data for the southernmost continuous segment of the Juan de Fuca Ridge shows that generalizations about tectonic and volcanic processes at spreading ridges must consider both the temporal scale of processes as well as the physical scales of observations if predictive models are to be successful. Comparison of the morphometric expression within the major hydrothermal vent area and the rest of the southernmost ridge segment suggests that the mapped distribution of hydrothermal vents may reflect the extent of survey effort rather than uniqueness of geologic setting.

Journal of Geophysical Research B: Solid Earth↗