USGS Science⌕ Search

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

The crustal structure of the Wrangellia Terrane along the East Glenn Highway, eastern‐southern Alaska

Recently acquired seismic refraction data from eastern‐southern Alaska provide new information on the structure and composition of the Wrangellia and adjacent terranes. The data comprise a 160‐km‐long refraction profile along the East Glenn (Tok‐Cutoff) Highway that was collected as part of the U.S. Geological Survey's multidisciplinary Trans‐Alaska Crustal Transect program. The upper 3 km of the Wrangellia terrane and associated rocks is characterized by low compressional wave velocities (V p = 1.9, 3.3, 4.6, 5.6 km s −1 ) and high‐velocity gradients common to most onshore seismic refraction profiles. There is also clear seismic expression of the West Fork fault system as a steep, down‐to‐the‐southwest fault that separates the Peninsular terrane to the southwest and the metamorphic complex of Gulkana River to the northeast. In contrast, no seismic expression occurs for the Paxson Lake fault, which separates the Wrangellia terrane from the metamorphic complex of Gulkana River. Adjacent to the Denali fault, within the Wrangellia terrane, two high‐velocity bodies (V p = 6.6 km s −1 ) occur in the upper crust. One of these extends to ∼10‐km depth and correlates with a late Paleozoic dioritic complex, suggesting that the Wrangellia terrane is at least 10 km thick in this part of Alaska. From 5 to 23 km depth, the crust appears seismically homogeneous, with velocity increasing from V p = 6.2 to V p = 6.6 km s −1 . Beneath this level, the crust is less well resolved, although evidence exists for a low‐velocity zone between 23 and 26 km and a possible southwest dipping interface at 35 km. No identifiable mantle refraction or reflection is observed, possibly indicating a crust as thick as 55 km. The relatively low seismic velocities in the upper 23 km of the crust compare favorably with laboratory‐measured velocities on pelitic schists and intermediate‐composition plutonic rocks (granites and granodiorites), both of which are recognized in Wrangellia. We interpret the seismic velocities to indicate that silicic‐to‐intermediate‐composition rocks are important constituents of the basement of this part of Wrangellia. Geologic evidence indicates that the Alaskan part of the Wrangellia terrane is a Paleozoic and Mesozoic island arc: our seismic evidence indicates it may have been built mostly on continental crust as opposed to the fragment of Wrangellia from Vancouver Island which was probably built on oceanic crust.

Alaska↗
USGS2012

Tertiary basin development and tectonic implications, Whipple Detachment System, Colorado River Extensional Corridor, California and Arizona

Colorado River extensional corridor. In the Mohave Mountains and Aubrey Hills of Arizona and the eastern Whipple Mountains of California near Parker Dam, these deposits comprise four unconformity-bounded sequences composed of locally derived epiclastic and volcanic rocks and the Peach Springs Tuff. The three older sequences represent syntectonic units that were deposited coeval with detachment faulting, and the fourth is interpreted to be postextensional. The sequences are correlated between four fault-bounded regions, which are the remnants of four different depositional basins. Similar sequences can be correlated over broad areas of the extensional corridor despite the general lack of widespread units. The basins developed in about the same positions, relative to each other and to volcanic sources, as they occupy at present. This is shown by gradational changes of pre-Tertiary rock types between regions, systematic variations in the abundance of magmatic units, and correlative volcanic units that occur in two adjacent regions. The basins formed in the early Miocene from segmentation of the upper crust into blocks bounded by high-angle faults that trended both parallel and perpendicular to the direction of extension and which were terminated at middle crustal depths by a low-angle detachment fault. Extreme rotation of one large crustal block, which constitutes the central Mohave Mountains, is recorded by a major unconformity in the lower Miocene section of one basin. Because coeval sections that formed in adjoining basins do not record this rotation, the underlying crustal blocks must have been separated by transfer faults that allowed them to rotate independently. These proposed transfer faults are represented at present by major faults with trends that parallel the direction of extension on the Whipple detachment system.

Arizona, California↗
USGS2012

Lithology and evolution of the crust-mantle boundary region in the southwestern Basin and Range Province

Mantle and crustal xenoliths from volcanic rocks in the southwestern Basin and Range province and Colorado Plateau Transition Zone reveal histories of episodic magmatism and deformation that have profoundly influenced the crustal structure of this region. Seismic transects in this area show a strongly reflective Moho of generally low relief, which, in the area of modern transects, consists of a thin zone (<2 km thick) of short reflectors. The upper mantle is transparent and has a P n of 7.8–8.0 km/s similar to much of the western United States. A lower crustal zone, 2–13 km thick, has variable internal reflectivity and a relatively low velocity of 6.6–6.8 km/s. Upper mantle peridotite xenoliths show both ductile and brittle deformational features and have structures and compositions affected by magmatic intrusion; intrusions form complex dike systems and extensive zones of grain boundary infiltration in peridotite xenoliths. Whereas melt infiltration preceded and followed ductile deformation, brittle deformation, represented by closely spaced joint systems and faults, followed ductile deformation and is related to the youngest magmatic episodes. These structural characteristics and high uppermost mantle temperature (∼1000°C) may combine to explain the relatively low P n . Alternating layers of ductily deformed and undeformed peridotites, with or without igneous intrusions, may contribute to the reflectivity of the Moho. Lower crustal xenoliths are dominantly igneous-textured pyroxenites and mafic to intermediate gabbros identical to the dikes in peridotite xenoliths. The crustal xenoliths also commonly are jointed, and in addition many show partial melting and have abundant cavities that probably were filled with CO 2 -rich fluids. These rocks are interpreted as products of underplated magmas that were fed through the mantle dike systems and may represent the lowest crustal unit identified in the seismic records. The mafic compositions and high densities of the crustal xenoliths indicate that the low velocity of the lower crust may be caused in part by fracture systems, partial melts, and high temperatures. Garnet granulite xenoliths from a locality with no mantle peridotite xenoliths probably represent crust of the region before late Miocene extension. Felsic granulite xenoliths from two localities have velocities like those of the two lower crustal units identified seismically and could be present in the modern crust as unequilibrated remnants of old crust. The preferred model for the evolution of the lower lithosphere is one in which extension affects the upper mantle as well as the crust and is overlapped in time by multiple magmatic episodes. The earliest magmatic events preceded extension, and later events accompanied and followed extension.

Arizona, California, Nevada↗
USGS2012

Keweenaw hot spot: Geophysical evidence for a 1.1 Ga mantle plume beneath the Midcontinent Rift System

The Proterozoic Midcontinent Rift System of North America is remarkably similar to Phanerozoic rifted continental margins and flood basalt provinces. Like the younger analogues, the volcanism within this older rift can be explained by decompression melting and rapid extrusion of igneous material during lithospheric extension above a broad, asthenospheric, thermal anomaly which we call the Keweenaw hot spot. Great Lakes International Multidisciplinary Program on Crustal Evolution seismic reflection profiles constrain end-member models of melt thickness and stretching factors, which yield an inferred mantle potential temperature of 1500°–1570°C during rifting. Combined gravity modeling and subsidence calculations are consistent with stretching factors that reached 3 or 4 before rifting ceased, and much of the lower crust beneath the rift consists of relatively high density intruded or underplated synrift igneous material. The isotopic signature of Keweenawan volcanic rocks, presented in a companion paper by Nicholson and Shirey (this issue), is consistent with our model of passive rifting above an asthenospheric mantle plume.

Midcontinent Rift System↗
USGS2012

Flank tectonics of Martian volcanoes

On the flanks of Olympus Mons is a series of terraces, concentrically distributed around the caldera. Their morphology and location suggest that they could be thrust faults caused by compressional failure of the cone. In an attempt to understand the mechanism of faulting and the possible influences of the interior structure of Olympus Mons, we have constructed a numerical model for elastic stresses within a Martian volcano. In the absence of internal pressurization, the middle slopes of the cone are subjected to compressional stress, appropriate to the formation of thrust faults. These stresses for Olympus Mons are ∼250 MPa. If a vacant magma chamber is contained within the cone, the region of maximum compressional stress is extended toward the base of the cone. If the magma chamber is pressurized, extensional stresses occur at the summit and on the upper slopes of the cone. For a filled but unpressurized magma chamber, the observed positions of the faults agree well with the calculated region of high compressional stress. Three other volcanoes on Mars, Ascraeus Mons, Arsia Mons, and Pavonis Mons, possess similar terraces. Extending our analysis to other Martian volcanoes, we find that only these three and Olympus Mons have flank stresses that exceed the compressional failure strength of basalt, lending support to the view that the terraces on all four are thrust faults.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Dark materials in Valles Marineris: Indications of the style of volcanism and magmatism on Mars

Rifting on the equatorial canyon system of Valles Marineris provides a unique view of the interior of the Martian crust to depths reaching 7 km, exposing several in situ bedrock units which testify to past volcanic and magmatic processes on Mars. A thick, regionally extensive deposit observed in Coprates and Juventae chasmata is interpreted on the basis of spectral reflectance, erosional morphology, and tendency for eolian mobilization to be composed of mafic glass, possibly an ancient Martian analogue of the lunar terra mantling deposits. Spectral mapping suggests that the dark floor‐covering materials in the lower canyons are derived from this unit. A series of cliffs in the Ophir Chasma wall rock is interpreted to be exposures of resistant bedrock; the spectral signature of this massive and uniform unit most closely resembles that of terrestrial mafic rocks altered to or coated by crystalline hematite. Application of computer mapping techniques to probable young volcanic materials in the central troughs yields an inferred distribution of volcanic activity consistent with the interpretation of extrusion along faults near the margins of the canyon floors. This result supports the hypothesis that the valles originated through tectonic extension.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Basal slip and mechanical anisotropy of biotite

The basal slip systems of biotite and their mechanical expressions have been investigated by shortening single crystals oriented to maximize and minimize shear stresses on (001). Samples loaded at 45° to (001) exhibit gentle external rotations associated with dislocation glide. High‐angle kink bands in these samples, unlike those developed in micas loaded parallel to (001), are limited to sample corners. Samples shortened perpendicular to (001) show no evidence of nonbasal slip and fail by fracture over all conditions tested. The mechanical response of biotite shortened at 45° to (001) is nearly perfectly elastic‐plastic; stress‐strain curves are characterized by a steep elastic slope, a sharply defined yield point, and continued deformation at low (mostly 1%. Stresses measured beyond the yield point are insensitive to confining pressure over the range 200 to 500 MPa and exhibit weak dependencies upon strain rate and temperature. Assuming an exponential relationship between differential stress σ and strain rate of the form , the data collected over strain rates and temperatures of 10−7 to 10−4 s−1 and 20° to 400°C, respectively, are best fit by an exponential constant α of 0.41±0.08 MPa‐1 and an activation energy of 82±13 kJ/mol. A power law fits the data equally well with = 18±4 and = 51±9 kJ/mol. Samples oriented favorably for slip in directions [100] and [110] are measurably weaker than those shortened at 45° to [010] and [310], consistent with the reported Burgers vectors 〈100〉, 1/2 〈110〉, and 1/2 〈110〉. The anisotropy of biotite is further revealed by contrasting these plastic strengths with results of samples deformed parallel and perpendicular to (001). Previous studies have shown that biotite loaded in the (001) plane is strong prior to the nucleation of kink bands. The strength of biotite shortened perpendicular to (001) exceeds that measured parallel to (001) and is pressure dependent. Application of the results to deformation within the continental crust suggests that biotite oriented favorably for slip is much weaker than most other silicates over a wide range of geologic conditions. Its presence within foliated rocks and shear zones may limit locally the stresses that can be supported.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

The 40Ar/39Ar thermochronology of the eastern Mojave Desert, California, and adjacent western Arizona with implications for the evolution of metamorphic core complexes

Mesozoic thickening and Cenozoic extension resulted in the juxtaposition of upper and middle crustal rocks in the eastern Mojave Desert, southeastern California and western Arizona. The application of 40 Ar/ 39 Ar thermochronology to rocks in this region provides information about the timing and nature of thrusting, plutonism, metamorphism, denudation, and detachment faulting. The 40 Ar/ 39 Ar ages of 175 to 125 Ma from the Clipper, Piute, Turtle, Mohave, Bill Williams, and Hualapai Mountains are interpreted to be the result of a middle Mesozoic thermal event(s) caused by crustal thickening. The 40 Ar/ 39 Ar data from the Clipper and Piute Mountains suggest that this thermal event was followed by a period of cooling at rates of 1°–5°C/m.y. Orogenesis culminated during the Late Cretaceous when rocks exposed in the Old Woman‐Piute, Chemehuevi, and Sacramento Mountains attained temperatures >500°C which reset the K–Ar systems of minerals from Proterozoic rocks. High‐grade metamorphism in the Old Woman Mountains area was caused by the intrusion of the Old Woman‐Piute batholith at 73±1 Ma. Cooling rates following batholith emplacement in the Old Woman Mountains were ∼100°C/m.y. between 73 and 70 Ma and 5°–10°C/m.y. from 70 to ∼30 Ma. Between 65 and 25 Ma the entire eastern Mojave Desert underwent a period of cooling at a rate of 2°–10°C/m.y. By 30 Ma, rocks exposed in the Old Woman‐Piute, Marble, Ship, Clipper, and Turtle Mountains were below ∼100°C. The 40 Ar/ 39 Ar ages from the Sacramento Mountains suggest that mylonitization caused by the onset of regional extension occurred at 23±1 Ma. When extension started in the Chemehuevi Mountains, rocks exposed in the southwestern and northeastern portions of footwall to the Chemehuevi detachment fault were at ∼180°C and ∼350°C, respectively. This suggests that the exposed part of the Chemehuevi detachment fault initiated at a dip of 5°–30° or as a series of higher‐angle faults that cut to a depth of 10–12 km and were later rotated to their present dip. Unroofing of the footwalls to detachment faults in the Sacramento and Chemehuevi Mountains resulted in average cooling rates of 10°–50°C/m.y. between 22 and 15 Ma.

California↗
USGS1996

Crustal structure of a transform plate boundary: San Francisco Bay and the central California continental margin

Wide-angle seismic data collected during the Bay Area Seismic Imaging Experiment provide new glimpses of the deep structure of the San Francisco Bay Area Block and across the offshore continental margin. San Francisco Bay is underlain by a veneer (<300 m) of sediments, beneath which P wave velocities increase rapidly from 5.2 km/s to 6.0 km/s at 7 km depth, consistent with rocks of the Franciscan subduction assemblage. The base of the Franciscan at-15-18 km depth is marked by a strong wide-angle reflector, beneath which lies an 8- to 10-km-thick lower crust with an average velocity of 6.75??0.15 km/s. The lower crust of the Bay Area Block may be oceanic in origin, but its structure and reflectivity indicate that it has been modified by shearing and/or magmatic intrusion. Wide-angle reflections define two layers within the lower crust, with velocities of 6.4-6.6 km/s and 6.9-7.3 km/s. Prominent subhorizontal reflectivity observed at near-vertical incidence resides principally in the lowermost layer, the top of which corresponds to the "6-s reflector" of Brocher et al. [1994]. Rheological modeling suggests that the lower crust beneath the 6-s reflector is the weakest part of the lithosphere; the horizontal shear zone suggested by Furlong et al. [1989] to link the San Andreas and Hayward/Calaveras fault systems may actually be a broad zone of shear deformation occupying the lowermost crust. A transect across the continental margin from the paleotrench to the Hayward fault shows a deep crustal structure that is more complex than previously realized. Strong lateral variability in seismic velocity and wide-angle reflectivity suggests that crustal composition changes across major transcurrent fault systems. Pacific oceanic crust extends 40-50 km landward of the paleotrench but, contrary to prior models, probably does not continue beneath the Salinian Block, a Cretaceous arc complex that lies west of the San Andreas fault in the Bay Area. The thickness (10 km) and high lower-crustal velocity of Pacific oceanic crust suggest that it was underplated by magmatism associated with the nearby Pioneer seamount. The Salinian Block consists of a 15-km-thick layer of velocity 6.0-6.2 km/s overlying a 5-km-thick, high-velocity (7.0 km/s) lower crust that may be oceanic crust, Cretaceous arc-derived lower crust, or a magmatically underplated layer. The strong structural variability across the margin attests to the activity of strike-slip faulting prior to and during development of the transcurrent Pacific/North American plate boundary around 29 Ma. Copyright 1996 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

Isotopic and trace element compositions of upper mantle and lower crustal xenoliths, Cima volcanic field, California: Implications for evolution of the subcontinental lithospheric mantle

Ultramafic and mafic xenoliths from the Cima volcanic field, southern California, provide evidence of episodic modification of the upper mantle and underplating of the crust beneath a portion of the southern Basin and Range province. The upper mantle xenoliths include spinel peridotite and anhydrous and hydrous pyroxenite, some cut by igneous-textured pyroxenite-gabbro veins and dikes and some by veins of amphibole ?? plagioclase. Igneous-textured pyroxenites and gabbros like the dike rocks also occur abundantly as isolated xenoliths inferred to represent underplated crust. Mineral and whole rock trace element compositions among and within the different groups of xenoliths are highly variable, reflecting multiple processes that include magma-mantle wall rock reactions, episodic intrusion and it filtration of basaltic melts of varied sources into the mantle wall rock, and fractionation. Nd, Sr, and Pb isotopic compositions mostly of clinopyroxene and plagioclase mineral separates show distinct differences between mantle xenoliths (??Nd = -5.7 to +3.4; 87Sr/86Sr = 0.7051 - 0.7073; 206Pb/204Pb = 19.045 - 19.195) and the igneous-textured xenoliths (??Nd = +7.7 to +11.7; 87Sr/86Sr = 0.7027 - 0.7036 with one carbonate-affected outlier at 0.7054; and 206Pb/204Pb = 18.751 - 19.068), so that they cannot be related. The igneous-textured pyroxenites and gabbros are similar in their isotopic compositions to the host basaltic rocks, which have ??Nd of+5.1 to +9.3; 87Sr/86Sr of 0.7028 - 0.7050, and 206Pb/204Pb of 18.685 - 21.050. The igneous-textured pyroxenites and gabbros are therefore inferred to be related to the host rocks as earlier cogenetic intrusions in the mantle and in the lower crust. Two samples of peridotite, one modally metasomatized by amphibole and the other by plagioclase, have isotopic compositions intermediate between the igneous-textured xenoliths and the mantle rock, suggesting mixing, but also derivation of the metasomatizing magmas from two separate and distinct sources. Sm-Nd two-mineral "isochrons" yield apparent ages for petrographically identical rocks believed to be coeval ranging from -0 to 113 ?? 26 Ma, indicating the unreliability of dating these rocks with this method. Amphibole and plagioclase megacrysts are isotopically like the host basalts and probably originate by mechanical breakup of veins comagmatic with the host basaltic rocks. Unlike other Basin and Range localities, Cima Cr-diopside group isotopic compositions do not overlap with those of the host basalts. Copyright 1997 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Fault slip and seismic moment of the 1700 Cascadia earthquake inferred from Japanese tsunami descriptions

The 1700 Cascadia earthquake attained moment magnitude 9 according to new estimates based on effects of its tsunami in Japan, computed coseismic seafloor deformation for hypothetical ruptures in Cascadia, and tsunami modeling in the Pacific Ocean. Reports of damage and flooding show that the 1700 Casscadia tsunami reached 1-5 m heights at seven shoreline sites in Japan. Three sets of estimated heights express uncertainty about location and depth of reported flooding, landward decline in tsunami heights from shorelines, and post-1700 land-level changes. We compare each set with tsunami heights computed from six Cascadia sources. Each source is vertical seafloor displacement calculated with a three-dimensional elastic dislocation model, for three sources the rupture extends the 1100 km length of the subduction zone and differs in width and shallow dip; for the other sources, ruptures of ordinary width extend 360-670 km. To compute tsunami waveforms, we use a linear long-wave approximation with a finite difference method, and we employ modern bathymetry with nearshore grid spacing as small as 0.4 km. The various combinations of Japanese tsunami heights and Cascadia sources give seismic moment of 1-9 ?? 1022 N m, equivalent to moment magnitude 8.7-9.2. This range excludes several unquantified uncertainties. The most likely earthquake, of moment magnitude 9.0, has 19 m of coseismic slip on an offshore, full-slip zone 1100 km long with linearly decreasing slip on a downdip partial-slip zone. The shorter rupture models require up to 40 m offshore slip and predict land-level changes inconsistent with coastal paleoseismological evidence. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Toggling of seismicity by the 1997 Kagoshima earthquake couplet: A demonstration of time-dependent stress transfer

Two M ??? 6 well-recorded strike-slip earthquakes struck just 4 km and 48 days apart in Kagoshima prefecture, Japan, in 1997, providing an opportunity to study earthquake interaction. Aftershocks are abundant where the Coulomb stress is calculated to have been increased by the first event, and they abruptly stop where the stress is dropped by the second event. This ability of the main shocks to toggle seismicity on and off argues that static stress changes play a major role in exciting aftershocks, whereas the dynamic Coulomb stresses, which should only promote seismicity, appear to play a secondary role. If true, the net stress changes from a sequence of earthquakes might be expected to govern the subsequent seismicity distribution. However, adding the stress changes from the two Kagoshima events does not fully capture the ensuing seismicity, such as its rate change, temporal decay, or migration away from the ends of the ruptures. We therefore implement a stress transfer model that incorporates rate/state friction, in which seismicity is treated as a sequence of independent nucleation events that are dependent on the fault slip, slip rate, and elapsed time since the last event. The model reproduces the temporal response of seismicity to successive stress changes, including toggling, decay, and aftershock migration. Nevertheless, the match of observed to predicted seismicity is quite imperfect, due perhaps to inadequate knowledge of several model parameters. However, to demonstrate the potential of this approach, we build a probabilistic forecast of larger earthquakes on the expected rate of small aftershocks, taking advantage of the large statistical sample the small shocks afford. Not surprisingly, such probabilities are highly time- and location-dependent: During the first decade after the main shocks, the seismicity rate and the chance of successive large shocks are about an order of magnitude higher than the background rate and are concentrated exclusively in the stress triggering zones. Copyright 2003 by the American Geophysical Uion.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Location of long-period events below Kilauea Volcano using seismic amplitudes and accurate relative relocation

We present methods for improving the location of long-period (LP) events, deep and shallow, recorded below Kilauea Volcano by the permanent seismic network. LP events might be of particular interest to understanding eruptive processes as their source mechanism is assumed to directly involve fluid transport. However, it is usually difficult or impossible to locate their source using traditional arrival time methods because of emergent wave arrivals. At Kilauea, similar LP waveform signatures suggest the existence of LP multiplets. The waveform similarity suggests spatially close sources, while catalog solutions using arrival time estimates are widely scattered beneath Kilauea's summit caldera. In order to improve estimates of absolute LP location, we use the distribution of seismic amplitudes corrected for station site effects. The decay of the amplitude as a function of hypocentral distance is used for inferring LP location. In a second stage, we use the similarity of the events to calculate their relative positions. The analysis of the entire LP seismicity recorded between January 1997 and December 1999 suggests that a very large part of the LP event population, both deep and shallow, is generated by a small number of compact sources. Deep events are systematically composed of a weak high-frequency onset followed by a low-frequency wave train. Aligning the low-frequency wave trains does not lead to aligning the onsets indicating the two parts of the signal are dissociated. This observation favors an interpretation in terms of triggering and resonance of a magmatic conduit. Instead of defining fault planes, the precise relocation of similar LP events, based on the alignment of the high-energy low-frequency wave trains, defines limited size volumes. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Theoretical analysis of tsunami generation by pyroclastic flows

Pyroclastic flows are a common product of explosive volcanism and have the potential to initiate tsunamis whenever thick, dense flows encounter bodies of water. We evaluate the process of tsunami generation by pyroclastic flow by decomposing the pyroclastic flow into two components, the dense underflow portion, which we term the pyroclastic debris flow, and the plume, which includes the surge and coignimbrite ash cloud parts of the flow. We consider five possible wave generation mechanisms. These mechanisms consist of steam explosion, pyroclastic debris flow, plume pressure, plume shear, and pressure impulse wave generation. Our theoretical analysis of tsunami generation by these mechanisms provides an estimate of tsunami features such as a characteristic wave amplitude and wavelength. We find that in most situations, tsunami generation is dominated by the pyroclastic debris flow component of a pyroclastic flow. This work presents information sufficient to construct tsunami sources for an arbitrary pyroclastic flow interacting with most bodies of water. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Groundwater flow, heat transport, and water table position within volcanic edifices: Implications for volcanic processes in the Cascade Range

The position of the water table within a volcanic edifice has significant implications for volcano hazards, geothermal energy, and epithermal mineralization. We have modified the HYDROTHERM numerical simulator to allow for a free-surface (water table) upper boundary condition and a wide range of recharge rates, heat input rates, and thermodynamic conditions representative of continental volcano-hydrothermal systems. An extensive set of simulations was performed on a hypothetical stratovolcano system with unconfined groundwater flow. Simulation results suggest that the permeability structure of the volcanic edifice and underlying material is the dominant control on water table elevation and the distribution of pressures, temperatures, and fluid phases at depth. When permeabilities are isotropic, water table elevation decreases with increasing heat flux and increases with increasing recharge, but when permeabilities are anisotropic, these effects can be much less pronounced. Several conditions facilitate the ascent of a hydrothermal plume into a volcanic edifice: a sufficient source of heat and magmatic volatiles at depth, strong buoyancy forces, and a relatively weak topography-driven flow system. Further, the plume must be connected to a deep heat source through a pathway with a time-averaged effective permeability ???1 ?? 10-16 m2, which may be maintained by frequent seismicity. Topography-driven flow may be retarded by low permeability in the edifice and/or the lack of precipitation recharge; in the latter case, the water table may be relatively deep. Simulation results were compared with observations from the Quaternary stratovolcanoes along the Cascade Range of the western United States to infer hydrothermal processes within the edifices. Extensive ice caps on many Cascade Range stratovolcanoes may restrict recharge on the summits and uppermost flanks. Both the simulation results and limited observational data allow for the possibility that the water table beneath the stratovolcanoes is relatively deep.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Location of seismic events and eruptive fissures on the Piton de la Fournaise volcano using seismic amplitudes

We present a method for locating the source of seismic events on Piton de la Fournaise. The method is based on seismic amplitudes corrected for station site effects using coda site amplification factors. Once corrected, the spatial distribution of amplitudes shows smooth and simple contours for many types of events, including rockfalls, long-period events and eruption tremor. On the basis of the simplicity of these distributions we develop inversion methods for locating their origins. To achieve this, the decrease of the amplitude as a function of the distance to the source is approximated by the decay either of surface or body waves in a homogeneous medium. The method is effective for locating rockfalls, long-period events, and eruption tremor sources. The sources of eruption tremor are usually found to be located at shallow depth and close to the eruptive fissures. Because of this, our method is a useful tool for locating fissures at the beginning of eruptions.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Near-field postseismic deformation associated with the 1992 Landers and 1999 Hector Mine, California, earthquakes

After the Landers earthquake (Mw = 7.3, 1992.489) a linear array of 10 monuments extending about 30 km N50??E on either side of the earthquake rupture plus a nearby offtrend reference monument were surveyed frequently by GPS until 2003.2. The array also spans the rupture of the subsequent Hector Mine earthquake (Mw = 7.1, 1999.792 . The pre-Landers velocities of monuments in the array relative to interior North America were estimated from earlier trilateration and very long baseline interferometry measurements. Except at the reference monument, the post-Landers velocities of the individual monuments in the array relaxed to their preseismic values within 4 years. Following the Hector Mine earthquake the velocities of the monuments relaxed to steady rates within 1 year. Those steady rates for the east components are about equal to the pre-Landers rates as is the steady rate for the north component of the one monument east of the Hector Mine rupture. However, the steady rates for the north components of the 10 monuments west of the rupture are systematically ???10 mm yr1 larger than the pre-Landers rates. The relaxation to a steady rate is approximately exponential with decay times of 0.50 ?? 0.10 year following the Landers earthquake and 0.32 ?? 0.18 year following the Hector Mine earthquake. The postearthquake motions of the Landers array following the Landers earthquake are not well approximated by the viscoelastic-coupling model of Pollitz et al. [2000]. A similar viscoelastic-coupling model [Pollitz et al., 2001] is more successful in representing the deformation after the Hector Mine earthquake.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Imaging the complexity of an active normal fault system: The 1997 Colfiorito (central Italy) case study

Six moderate magnitude earthquakes (5 < Mw < 6) ruptured normal fault segments of the southern sector of the North Apennine belt (central Italy) in the 1997 Colfiorito earthquake sequence. We study the progressive activation of adjacent and nearby parallel faults of this complex normal fault system using ???1650 earthquake locations obtained by applying a double-difference location method, using travel time picks and waveform cross-correlation measurements. The lateral extent of the fault segments range from 5 to 10 km and make up a broad, ???45 km long, NW trending fault system. The geometry of each segment is quite simple and consists of planar faults gently dipping toward SW with an average dip of 40??-45??. The fault planes are not listric but maintain a constant dip through the entire seismogenic volume, down to 8 km depth. We observe the activation of faults on the hanging wall and the absence of seismicity in the footwall of the structure. The observed fault segmentation appears to be due to the lateral heterogeneity of the upper crust: preexisting thrusts inherited from Neogene's compressional tectonic intersect the active normal faults and control their maximum length. The stress tensor obtained by inverting the six main shock focal mechanisms of the sequence is in agreement with the tectonic stress active in the inner chain of the Apennine, revealing a clear NE trending extension direction. Aftershock focal mechanisms show a consistent extensional kinematics, 70% of which are mechanically consistent with the main shock stress field.

Journal of Geophysical Research B: Solid Earth↗