USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Oceans”

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 1,657 records · Page 92Linked to original sources

Rhyolitic calderas of the Yukon-Tanana Terrane, east central Alaska: volcanic remnants of a mid-Cretaceous magmatic arc

Four large but poorly exposed rhyolitic calderas are present in the Yukon-Tanana terrane (YTT) in east central Alaska. At least two are mid-Cretaceous in age (~93 Ma). Similar volcanic rocks, the South Fork Volcanics, occur northeast of the Tintina fault in Yukon Territory. Evidence for the calderas consists of thick deposits of devitrified crystal- and lithic-rich densely welded tuff, interpreted as caldera fill, associated with lava domes or shallow intrusive rocks. Coeval outflow sheets have been largely stripped by erosion. The calderas are preserved within a northeast trending depression extending across the axis of the elongate mid-Cretaceous plutonic province. Trace element abundances in andesites and rhyolites associated with the caldera structures are similar to those of volcanic and plutonic rocks of subduction-related magmatic arcs developed on continental crust and thus are suggestive of formation in such an environment. Late Cretaceous and early Tertiary igneous rocks in the YTT near the calderas are interpreted to have been emplaced in a more extensional setting when the subduction-related magmatic front was farther oceanward.

Alaska↗

Comparison of ground motion from tremors and explosions in deep gold mines

Seismic body waves, from tamped chemical explosions, two with yields of 50 and one of 150 kg, were compared with corresponding data from three mining-induced tremors with a view to testing methods of discriminating between the two types of events. Detonated at depths of about 2 km, all three explosions generated P waves for which the low-frequency spectral asymptotes agree well with corresponding results scaled down from nuclear shots at the Nevada Test Site. For the two smaller explosions, recorded underground at hypocentral distances ranging from 234 to 871 m, the P wave corner frequencies, defined by the intersections of the low- and high-frequency asymptotes, are also in fine agreement with scaled-down results from the Nevada Test Site; for the 150-kg explosion, which was recorded only at the surface, it appears that 2.39 km of upward propagation caused the corner frequency to be at least a factor of 4 lower than anticipated due to attenuation. All three explosions generated S waves that appear to be a consequence of deviatoric stress release in the immediate environs of the explosions. The three tremors analyzed here generated P and S waves whose spectra and source parameters agree well with standard earthquake source models. With regard to discrimination based on P wave spectra, we conclude that for events of fixed low-frequency spectral asymptotes, the explosions typically have higher corner frequencies than tremors or earthquakes, although counterexamples certainly exist. Interestingly, the 150-kg explosion was identified as such on the basis of P and S wave polarities that are incompatible with the normally expected double-couple source model; instead, these initial motions are consistent with an explosion in conjunction with normal faulting. The body wave spectra of this explosion and those of a nearby tremor, however, were indistinguishable.

Journal of Geophysical Research Solid Earth↗

Crustal structure of the northwestern Basin and Range Province from the 1986 Program for Array Seismic Studies of the Continental Lithosphere Seismic Experiment

A portion of northwestern Nevada was imaged to determine the crustal structure and to assess reported differences between refraction versus reflection determinations of Moho depth and how the crustal composition and structure has been influenced by volcanic and extension mechanisms. Interpretation of the refraction/wide-angle reflection data suggests that the crust is fairly uniform in thickness and varies by less than 5 km over the 280 km east-west profile and 3 km over its 220 km north-south length. The velocity structure is characterized by five layers: 1) an uppermost crust, composed of sedimentary rocks and basement that has an average velocity of 5.7 km s-1; 2) a middle crust that extends to a depth of 18-22 with an average velocity of 6.1 km s-1; 3) a 10-12 km thick lower crust with an average velocity of 6.6 km s-1; 4) a 2-5 km thick transitional crust-mantle boundary defined by a 7.6 km s-1 velocity; and 5) an upper mantle with an average Pn velocity of 7.9-8.0 km s-1. A uniform upper mantle composition across the Basin and Range is suggested and the homogeneity of the velocity structure beneath the western Basin and Range argues for a youthful Moho and crust that has been reworked by province-wide late Cenozoic extension, episodic magmatism, and underplating. -from Authors

Journal of Geophysical Research↗

Infragravity waves over a natural barred profile

Measurements of cross-shore flow were made across the surf zone during a storm as a nearshore bar became better developed and migrated offshore. Measured infragravity band spectra were compared to synthetic spectra calculated numerically over the natural barred profile assuming a white run-up spectrum of leaky mode or high-mode edge waves. As in earlier studies, the spectra compared closely; however, for some frequencies the energy of the measured spectrum exceeded the energy of the synthetic spectrum, suggesting that the run-up spectrum was not white but had dominating frequencies. Utilizing cross-shore flow data and synthetic spectra from a number of cross-shore locations, an equivalent run-up spectrum was calculated for each day. On the first day of the storm, the equivalent run-up spectrum indicated a dominant wave that had a node in velocity reasonably close to the bar crest. Later during the storm, when the bar had migrated farther offshore, there was no evidence for a dominant motion having a velocity node at the bar crest. The structure of the equivalent run-up spectrum compared well with spectra of direct measurements of run-up obtained several hundred meters away. We have no clear evidence in support of the theory that infragravity waves might form or force the offshore migration of a bar. To confirm this finding, longer records obtained synoptically over a developing bar are required. The dominant wave observed early in the storm was consistent with Symond and Bowen's (1984) theoretical prediction of resonant amplification of discrete frequencies over a barred profile.

Journal of Geophysical Research Solid Earth↗

Sulfur dioxide and particles in quiescent volcanic plumes from Poás, Arenal, and Colima Volcanos, Costa Rica and Mexico

Measurements of SO 2 emission rates and concentrations and of particle distribution, size, shape, and composition were made in quiescent volcanic plumes emitted into the troposphere from Poás and Arenal volcanos, Costa Rica, and Colima volcano, Mexico. SO 2 emission rates were 700±180 metric tons per day (t/d) for Poás, 210±30 t/d for Arenal, and 320±50 t/d for Colima. The concentrations of SO 2 calculated from the COSPEC/lidar data were 5–380 ppb. Concentrations of SO 2 measured directly by flame photometry were 10–250 ppb. Particles collected in the plumes with a quartz crystal microbalance impactor were mostly less than 3 μm in diameter and consisted of droplets of dilute sulfur-bearing solutions and minor amounts of larger silicate particles coated with a sulfur-bearing film or crust. Total particle concentrations were 4.7 μg/m 3 for Poás and 18.8 μg/m 3 for Colima. Comparison of concentrations of SO 2 in the plumes with gas samples collected at fumaroles on the ground suggests that the plumes are diluted by the atmosphere by factors of up to 10 5 .

Journal of Geophysical Research Atmospheres↗

Influence of tidal range on the stability of coastal marshland

Early comparisons between rates of vertical accretion and sea level rise across marshes in different tidal ranges inspired a paradigm that marshes in high tidal range environments are more resilient to sea level rise than marshes in low tidal range environments. We use field‐based observations to propose a relationship between vegetation growth and tidal range and to adapt two numerical models of marsh evolution to explicitly consider the effect of tidal range on the response of the marsh platform channel network system to accelerating rates of sea level rise. We find that the stability of both the channel network and vegetated platform increases with increasing tidal range. Our results support earlier hypotheses that suggest enhanced stability can be directly attributable to a vegetation growth range that expands with tidal range. Accretion rates equilibrate to the rate of sea level rise in all experiments regardless of tidal range, suggesting that comparisons between accretion rate and tidal range will not likely produce a significant relationship. Therefore, our model results offer an explanation to widely inconsistent field‐based attempts to quantify this relationship while still supporting the long‐held paradigm that high tidal range marshes are indeed more stable.

Journal of Geophysical Research↗

Role of geophysics in identifying and characterizing sites for high-level nuclear waste repositories.

Evaluation of potential high-level nuclear waste repository sites is an area where geophysical capabilities and limitations may significantly impact a major governmental program. Since there is concern that extensive exploratory drilling might degrade most potential disposal sites, geophysical methods become crucial as the only nondestructive means to examine large volumes of rock in three dimensions. Characterization of potential sites requires geophysicists to alter their usual mode of thinking: no longer are anomalies being sought, as in mineral exploration, but rather their absence. Thus the size of features that might go undetected by a particular method take on new significance. Legal and regulatory considerations that stem from this different outlook, most notably the requirements of quality assurance (necessary for any data used in support of a repository license application), are forcing changes in the manner in which geophysicists collect and document their data. -Authors

Journal of Geophysical Research↗

Crustal structure of the western New England Appalachians and the Adirondack Mountains

Presents an interpretation of the crustal velocity structure based on a seismic refraction/wide-angle reflection experiment in eastern North America extending from the Adirondacks in New York State through the northern Appalachians in Vermont and New Hampshire to central Maine. Modeling of the eastern portion of the profile within the New England Appalachians shows a subhorizontal layered crust with upper crustal velocities ranging from 5.5 to 6.2 km/s, a midcrustal velocity of 6.4 km/s, and a lower crustal velocity of approximately 6.8 km/s. Crustal thickness increases from 36 km beneath Maine to 40 km in Vermont. This experiment provides further evidence for the obduction of the allochthonous western Appalachian units onto Grenvillian crust above the azone of detachment that penetrates at least to midcrustal depths and was the locus of successive Paleozoic thrusting. -from Authors

Journal of Geophysical Research↗

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↗

Flow of variably fluidized granular masses across three-dimensional terrain 2. Numerical predictions and experimental tests

Numerical solutions of the equations describing flow of variably fluidized Coulomb mixtures predict key features of dry granular avalanches and water-saturated debris flows measured in physical experiments. These features include time-dependent speeds, depths, and widths of flows as well as the geometry of resulting deposits. Threedimensional (3-D) boundary surfaces strongly influence flow dynamics because transverse shearing and cross-stream momentum transport occur where topography obstructs or redirects motion. Consequent energy dissipation can cause local deceleration and deposition, even on steep slopes. Velocities of surge fronts and other discontinuities that develop as flows cross 3-D terrain are predicted accurately by using a Riemann solution algorithm. The algorithm employs a gravity wave speed that accounts for different intensities of lateral stress transfer in regions of extending and compressing flow and in regions with different degrees of fluidization. Field observations and experiments indicate that flows in which fluid plays a significant role typically have high-friction margins with weaker interiors partly fluidized by pore pressure. Interaction of the strong perimeter and weak interior produces relatively steep-sided, flat-topped deposits. To simulate these effects, we compute pore pressure distributions using an advection-diffusion model with enhanced diffusivity near flow margins. Although challenges remain in evaluating pore pressure distributions in diverse geophysical flows, Riemann solutions of the depthaveraged 3-D Coulomb mixture equations provide a powerful tool for interpreting and predicting flow behavior. They provide a means of modeling debris flows, rock avalanches, pyroclastic flows, and related phenomena without invoking and calibrating Theological parameters that have questionable physical significance.

Journal of Geophysical Research B: Solid Earth↗

Mars Global Surveyor observations of Martian fretted terrain

The Martian fretted terrain between latitudes 30° and 50°N and between 315° and 360°W has been reexamined in light of new Mars Orbiter Camera (MOC) and Mars Orbiter Laser Altimeter (MOLA) data from Mars Global Surveyor. Much of the terrain in the 30°–50° latitude belt in both hemispheres has a characteristic stippled or pitted texture at MOC (1.5 m) scale. The texture appears to result from partial removal of a formerly smooth, thin deposit as a result of sublimation and deflation. A complex history of deposition and exhumation is indicated by remnants of a former, thicker cover of layered deposits. In some hollows and on some slopes, particularly those facing the pole, are smooth textured deposits outlined by an outward facing escarpment. Throughout the study area are numerous escarpments with debris flows at their base. The escarpments typically have slopes in the 20°–30° range. At the base of the escarpment is commonly a deposit with striae oriented at right angles to the escarpment. Outside this deposit is the main debris apron with a surface that typically slopes 2°–3° and complex surface textures suggestive of compression, sublimation, and deflation. The presence of undeformed impact craters indicates that the debris flows are no longer forming. Fretted valleys contain lineated fill and are poorly graded. They likely form from fluvial valleys that were initially like those elsewhere on the planet but were subsequently widened and filled by the same mass-wasting processes that formed the debris aprons. Slope reversals indicate that downvalley flow of the lineated fill is minor. The ubiquitous presence of breaks in slope formed by mass wasting and the complex surface textures that result from mass wasting, deflation, and sublimation decreases the recognizability of the shorelines formerly proposed for this area.

Journal of Geophysical Research E: Planets↗

Visible/near-infrared spectra of experimentally shocked plagioclase feldspars

High shock pressures cause structural changes in plagioclase feldspars such as mechanical fracturing and disaggregation of the crystal lattice at submicron scales, the formation of diaplectic glass (maskelynite), and genuine melting. Past studies of visible/ near-infrared spectra of shocked feldspars demonstrated few spectral variations with pressure except for a decrease in the depth of the absorption feature near 1250-1300 nm and an overall decrease in reflectance. New visible/near-infrared spectra (400-2500 nm) of experimentally shocked (17-56 GPa) albite- and anorthite-rich rock powders demonstrate similar trends, including the loss of minor hydrated mineral bands near 1410, 1930, 2250, and 2350 nm. However, the most interesting new observations are increases in reflectance at intermediate pressures, followed by subsequent decreases in reflectance at higher pressures. The amount of internal scattering and overall sample reflectance is controlled by the relative proportions of micro-fractures, submicron grains, diaplectic glass, and melts formed during shock metamorphism. We interpret the observed reflectance increases at intermediate pressures to result from progressively larger proportions of submicron feldspar grains and diaplectic glass. The ensuing decreases in reflectance occur after diaplectic glass formation is complete and the proportion of genuine melt inclusions increases. The pressure regimes over which these reflectance variations occur differ between albite and anorthite, consistent with thermal infrared spectra of these samples and previous studies of shocked feldspars. These types of spectral variations associated with different peak shock pressures should be considered during interpretation and modeling of visible/near-infrared remotely sensed spectra of planetary and asteroidal surfaces.

Journal of Geophysical Research E: Planets↗

Aeolian processes in Proctor Crater on Mars: Sedimentary history as analyzed from multiple data sets

Proctor Crater is a 150 km diameter crater in Noachis Terra, within the southern highlands of Mars. The analysis leading to the sedimentary history incorporates several data sets including imagery, elevation, composition, and thermal inertia, mostly from the Mars Global Surveyor mission. The resulting stratigraphy reveals that the sedimentary history of Proctor Crater has involved a complex interaction of accumulating and eroding sedimentation. Aeolian features spanning much of the history of the crater interior dominate its surface, including large erosional pits, stratified beds of aeolian sediment, sand dunes, erosional and depositional streaks, dust devil tracks, and small bright bed forms that are probably granule ripples. Long ago, up to 450 m of layered sediment filled the crater basin, now exposed in eroded pits on the crater floor. These sediments are probably part of an ancient deposit of aeolian volcaniclastic material. Since then, some quantity of this material has been eroded from the top layers of the strata. Small, bright dune forms lie stratigraphically beneath the large dark dune field. Relative to the large dark dunes, the bright bed forms are immobile, although in places, their orientations are clearly influenced by the presence of the larger dunes. Their prevalence in the crater and their lack of compositional and thermal distinctiveness relative to the crater floor suggests that these features were produced locally from the eroding basin fill. Dust devil tracks form during the spring and summer, following a west-southwesterly wind. Early in the spring the dust devils are largely restricted to dark patches of sand. As the summer approaches, dust devil tracks become more plentiful and spread to the rest of the crater floor, indicating that the entire region acquires an annual deposit of dust that is revealed by seasonal dust devils. The dark dunes contain few dust devil tracks, suggesting that accumulated dust is swept away directly by saltation, rather than by the passage of dust devils. Spectral deconvolution indicates that the dark dunes have infrared spectra consistent with basalt-like materials. The average thermal inertia calculated from Thermal Emission Spectrometer bolometric temperatures is 277 ?? 17 J m-2 s-0.5 K-1, leading to an effective grain size of 740 ?? 170 ??m, which is consistent with coarse sand and within the range expected for Martian sand. The coarse sand that composes the large dune field may have originated from outside the crater, saltating in from the southwest. Most of the transport pathway that delivered this sand to the dune field has since been eroded away or buried. The sand was transported to the east center of the crater floor, where beneath the present-day dunes a 50 m high mound of sand has accumulated. Dune slip faces indicate a wind regime consisting of three opposing winds. Some of these wind directions are correlated with the orientations of dust devil tracks and bright bed forms. The combination of a tall mound of sand and three opposing winds is consistent with a convergent wind regime, which produces the large reversing transverse and star dunes that dominate the dune field. The dark dunes have both active slip faces and seemingly inactive slip faces, suggesting that the dunes vary spatially in their relative activity. Nevertheless, the aeolian activity that has dominated the history of Proctor Crater still continues today. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

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↗

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↗

Athena Mars rover science investigation

Each Mars Exploration Rover carries an integrated suite of scientific instruments and tools called the Athena science payload. The primary objective of the Athena science investigation is to explore two sites on the Martian surface where water may once have been present, and to assess past environmental conditions at those sites and their suitability for life. The remote sensing portion of the payload uses a mast called the Pancam Mast Assembly (PMA) that provides pointing for two instruments: the Panoramic Camera (Pancam), and the Miniature Thermal Emission Spectrometer (Mini-TES). Pancam provides high-resolution, color, stereo imaging, while Mini-TES provides spectral cubes at mid-infrared wavelengths. For in-situ study, a five degree-of-freedom arm called the Instrument Deployment Device (IDD) carries four more tools: a Microscopic Imager (MI) for close-up imaging, an Alpha Particle X-Ray Spectrometer (APXS) for elemental chemistry, a Mo??ssbauer Spectrometer (MB) for the mineralogy of Fe-bearing materials, and a Rock Abrasion Tool (RAT) for removing dusty and weathered surfaces and exposing fresh rock underneath. The payload also includes magnets that allow the instruments to study the composition of magnetic Martian materials. All of the Athena instruments have undergone extensive calibration, both individually and using a set of geologic reference materials that are being measured with all the instruments. Using a MER-like rover and payload in a number of field settings, we have devised operations processes that will enable us to use the MER rovers to formulate and test scientific hypotheses concerning past environmental conditions and habitability at the landing sites. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

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↗

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↗