USGS Science⌕ Search

SEARCH · USGS Science

Results for “The Earth Observer”

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,189 records · Page 66Linked to original sources

Slip history and dynamic implications of the 1999 Chi-Chi, Taiwan, earthquake

[1] We investigate the rupture process of the 1999 Chi‐Chi, Taiwan, earthquake using extensive near‐source observations, including three‐component velocity waveforms at 36 strong motion stations and 119 GPS measurements. A three‐plane fault geometry derived from our previous inversion using only static data [ Ji et al. , 2001 ] is applied. The slip amplitude, rake angle, rupture initiation time, and risetime function are inverted simultaneously with a recently developed finite fault inverse method that combines a wavelet transform approach with a simulated annealing algorithm [ Ji et al. , 2002b ]. The inversion results are validated by the forward prediction of an independent data set, the teleseismic P and SH ground velocities, with notable agreement. The results show that the total seismic moment release of this earthquake is 2.7 × 10 20 N m and that most of the slip occurred in a triangular‐shaped asperity involving two fault segments, which is consistent with our previous static inversion. The rupture front propagates with an average rupture velocity of ∼2.0 km s −1 , and the average slip duration (risetime) is 7.2 s. Several interesting observations related to the temporal evolution of the Chi‐Chi earthquake are also investigated, including (1) the strong effect of the sinuous fault plane of the Chelungpu fault on spatial and temporal variations in slip history, (2) the intersection of fault 1 and fault 2 not being a strong impediment to the rupture propagation, and (3) the observation that the peak slip velocity near the surface is, in general, higher than on the deeper portion of the fault plane, as predicted by dynamic modeling.

Journal of Geophysical Research B: Solid Earth↗

Inducing in situ, nonlinear soil response applying an active source

[1] It is well known that soil sites have a profound effect on ground motion during large earthquakes. The complex structure of soil deposits and the highly nonlinear constitutive behavior of soils largely control nonlinear site response at soil sites. Measurements of nonlinear soil response under natural conditions are critical to advancing our understanding of soil behavior during earthquakes. Many factors limit the use of earthquake observations to estimate nonlinear site response such that quantitative characterization of nonlinear behavior relies almost exclusively on laboratory experiments and modeling of wave propagation. Here we introduce a new method for in situ characterization of the nonlinear behavior of a natural soil formation using measurements obtained immediately adjacent to a large vibrator source. To our knowledge, we are the first group to propose and test such an approach. Employing a large, surface vibrator as a source, we measure the nonlinear behavior of the soil by incrementally increasing the source amplitude over a range of frequencies and monitoring changes in the output spectra. We apply a homodyne algorithm for measuring spectral amplitudes, which provides robust signal-to-noise ratios at the frequencies of interest. Spectral ratios are computed between the receivers and the source as well as receiver pairs located in an array adjacent to the source, providing the means to separate source and near-source nonlinearity from pervasive nonlinearity in the soil column. We find clear evidence of nonlinearity in significant decreases in the frequency of peak spectral ratios, corresponding to material softening with amplitude, observed across the array as the source amplitude is increased. The observed peak shifts are consistent with laboratory measurements of soil nonlinearity. Our results provide constraints for future numerical modeling studies of strong ground motion during earthquakes.

Journal of Geophysical Research B: Solid Earth↗

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

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

Journal of Geophysical Research Solid Earth↗

Relaxation response of critically stressed macroscale surficial rock sheets

Rock environments both underground and on Earth’s surface show indications of energetic macroscale fracture. In tunnels and excavations, these manifest as rockbursts—energetic explosions of rock that can damage engineering projects, and may pose ongoing financial and safety risk as rock stresses adjust during post-failure relaxation. In natural settings at the surface, evidence for rockbursts exist in the form of tent-like structures of ruptured exfoliation sheets, but few direct observations of such events exist, precluding the analysis of how natural rock formations may evolve after rupture. Here we investigate the post-failure evolution of a granitic rock dome following rapid fracture events (i.e., surficial rockbursts) that occurred in California, USA during 2014–2016. Building upon previous work that showed a thermal stress origin for the observed fracturing, we investigate the return to background stress conditions (i.e., stress relaxation) observed in both short- (week, month) and long-term (multi-year) rock deformation trends. Acoustic emissions, deformation, and environmental monitoring data indicate that partially detached rock sheets forming the surface of the dome undergo fracture aperture closing during cooling periods, concurrent with reduction of rock stress by the source of forcing (i.e., thermal stress). However, with sufficient critical and/or subcritical fracture, our observations also show that rock sheets can become decoupled from the source of stress, resulting in a long-term return to background stress conditions. Our results provide insight into the cyclic and likely ephemeral nature of rock fracture in surficial rock domes, as well as in underground rockburst environments.

Rock Mechanics and Rock Engineering↗

Mars and Earth: Comparison of cold-climate features

On Earth, glacial and periglacial features are common in areas of cold climate. On Mars, the temperature of the present-day surface is appropriate for permafrost, and the presence of water is suspected from data relating to the outgassing of the planet, from remote-sensing measurements over the polar caps and elsewhere on the Martian surface, and from recognition of fluvial morphological features such as channels. These observations and the possibility that ice could be in equilibrium with the atmosphere in the high latitudes north and south of ±40° latitude suggest that glacial and periglacial features should exist on the planet. Morphological studies based mainly on Viking pictures indicate many features that can be attributed to the action of ice. Among these features are extensive talus aprons; debris avalanches; flows that resemble glaciers or rock glaciers; ridges that look like moraines; various types of patterned ground, scalloped scarps, and chaotically collapsed terrain that could be attributed to thermokarst processes; and landforms that may reflect the interaction of volcanism and ice.

Icarus↗

Rare earth, major, and trace element composition of Monterey and DSDP chert and associated host sediment: Assessing the influence of chemical fractionation during diagenesis

Chert and associated host sediments from Monterey Formation and Deep Sea Drilling Project (DSDP) sequences were analyzed in order to assess chemical behavior during diagenesis of biogenic sediments. The primary compositional contrast between chert and host sediment is a greater absolute SiO 2 concentration in chert, often with final SiO 2 ≥ 98 wt%. This contrast in SiO 2 (and Si Al "> SiAl ) potentially reflects precursor sediment heterogeneity, diagenetic chemical fractionation, or both. SiO 2 concentrations and Si Al "> SiAl ratios in chert are far greater than in modern siliceous oozes, however and often exceed values in acid-cleaned diatom tests. Compositional contrasts between chert and host sediment are also orders-of-magnitude greater than between multiple samples of the host sediment. Calculations based on the initial composition of adjacent host, observed porosity reductions from host to chert and a postulated influx of pure SiO 2 , construct a chert composition which is essentially identical to observed SiO 2 values in chert. Thus, precursor heterogeneity does not seem to be the dominant factor influencing the current chert composition for the key elements of interest. In order to assess the extent of chemical fractionation during diagenesis, we approximate the precursor composition by analyzing host sediments adjacent to the chert. The SiO 2 concentration contrast seems caused by biogenic SiO 2 dissolution and transport from the local adjacent host sediment and subsequent SiO 2 reprecipitation in the chert. Along with SiO 2 , other elements are often added (with respect to Al) to Monterey and DSDP chert during silicification, although absolute concentrations decrease. The two Monterey quartz chert nodules investigated, in contrast to the opal-CT and quartz chert lenses, formed primarily by extreme removal of carbonate and phosphate, thereby increasing relative SiO 2 concentrations. DSDP chert formed by both carbonate/phosphate dissolution and SiO 2 addition from the host. Manganese is fractionated during chert formation, resulting in MnO Al 2 O 3 "> MnOAl2O3 ratios that no longer record the depositional signal of the precursor sediment. REE data indicate only subtle diagenetic fractionation across the rare earth series. Ce Ce ∗ "> CeCe* values do not change significantly during diagenesis of either Monterey or DSDP chert. Eu Eu ∗ "> EuEu* decreases slightly during formation of DSDP chert. La n Yb n "> LanYbn is affected only minimally as well. During formation of one Monterey opal-CT chert lens, REE Al "> REEAl ratios show subtle distribution changes at Gd and to a lesser extent near Nd and Ho. REE compositional contrasts between diagenetic states of siliceous sediment and chert are of a vastly smaller scale than has been noted between different depositional environments of marine sediment, indicating that the paleoenvironmental REE signature is not obscured by diagenetic overprinting.

Geochimica et Cosmochimica Acta↗

Self-similar rupture implied by scaling properties of volcanic earthquakes occurring during the 2004-2008 eruption of Mount St. Helens, Washington

We analyze a group of 6073 low-frequency earthquakes recorded during a week-long temporary deployment of broadband seismometers at distances of less than 3 km from the crater at Mount St. Helens in September of 2006. We estimate the seismic moment ( M 0 ) and spectral corner frequency ( f 0 ) using a spectral ratio approach for events with a high signal-to-noise (SNR) ratio that have a cross-correlation coefficient of 0.8 or greater with at least five other events. A cluster analysis of cross-correlation values indicates that the group of 421 events meeting the SNR and cross-correlation criteria forms eight event families that exhibit largely self-similar scaling. We estimate the M 0 and f 0 values of the 421 events and calculate their static stress drop and scaled energy ( E R / M 0 ) values. The estimated values suggest self-similar scaling within families, as well as between five of eight families (i.e., and constant). We speculate that differences in scaled energy values for the two families with variable scaling may result from a lack of resolution in the velocity model. The observation of self-similar scaling is the first of its kind for such a large group of low-frequency volcanic tectonic events occurring during a single active dome extrusion eruption.

Journal of Geophysical Research B: Solid Earth↗

Thresholds for post-wildfire debris flows: Insights from the Pinal Fire, Arizona, USA

Wildfire significantly alters the hydrologic properties of a burned area, leading to increases in overland flow, erosion, and the potential for runoff-generated debris flows. The initiation of debris flows in recently burned areas is well-characterized by rainfall intensity-duration (ID) thresholds. However, there is currently a paucity of data quantifying the rainfall intensities required to trigger post-wildfire debris flows, which limits our understanding of how and why rainfall ID thresholds vary in different climatic and geologic settings. In this study, we monitored debris-flow activity following the Pinal Fire in central Arizona, which differs from both a climatic and hydrogeomorphic perspective from other regions in the western U.S. where ID thresholds for post-wildfire debris flows are well-established, namely the Transverse Ranges of southern CA. Since the peak rainfall intensity within a rainstorm may exceed the rainfall intensity required to trigger a debris flow, the development of robust rainfall ID thresholds requires knowledge of the timing of debris flows within rainstorms. Existing post-wildfire debris-flow studies in Arizona only constrain the peak rainfall intensity within debris-flow-producing storms, which may far exceed the intensity that actually triggered the observed debris flow. In this study, we used pressure transducers within 5 burned drainage basins to constrain the timing of debris flows within rainstorms. Rainfall ID thresholds derived here from triggering rainfall intensities are, on average, 22 mm/h lower than ID thresholds derived under the assumption that the triggering intensity is equal to the maximum rainfall intensity recorded during a rainstorm. We then use a hydrologic model to demonstrate that the magnitude of the 15-minute rainfall ID threshold at the Pinal Fire site is associated with the rainfall intensity required to exceed a recently proposed dimensionless discharge threshold for debris-flow initiation. Model results further suggest that previously observed differences in regional ID thresholds between Arizona and the San Gabriel Mountains of southern CA may be attributed, in large part, to differences in the hydraulic properties of burned soils.

Arizona↗

Runoff-generated debris flows: observations and modeling of surge initiation, magnitude, and frequency

Runoff during intense rainstorms plays a major role in generating debris flows in many alpine areas and burned steeplands. Yet compared to debris flow initiation from shallow landslides, the mechanics by which runoff generates a debris flow are less understood. To better understand debris flow initiation by surface water runoff, we monitored flow stage and rainfall associated with debris flows in the headwaters of two small catchments: a bedrock-dominated alpine basin in central Colorado (0.06 km 2 ) and a recently burned area in southern California (0.01 km 2 ). We also obtained video footage of debris flow initiation and flow dynamics from three cameras at the Colorado site. Stage observations at both sites display distinct patterns in debris flow surge characteristics relative to rainfall intensity (I). We observe small, quasiperiodic surges at low I; large, quasiperiodic surges at intermediate I; and a single large surge followed by small-amplitude fluctuations about a more steady high flow at high I. Video observations of surge formation lead us to the hypothesis that these flow patterns are controlled by upstream variations in channel slope, in which low-gradient sections act as “sediment capacitors,” temporarily storing incoming bed load transported by water flow and periodically releasing the accumulated sediment as a debris flow surge. To explore this hypothesis, we develop a simple one-dimensional morphodynamic model of a sediment capacitor that consists of a system of coupled equations for water flow, bed load transport, slope stability, and mass flow. This model reproduces the essential patterns in surge magnitude and frequency with rainfall intensity observed at the two field sites and provides a new framework for predicting the runoff threshold for debris flow initiation in a burned or alpine setting.

California;Colorado↗

Identifying calcium sources at an acid deposition-impacted spruce forest: A strontium isotope, alkaline earth element multi-tracer approach

Depletion of calcium from forest soils has important implications for forest productivity and health. Ca is available to fine feeder roots from a number of soil organic and mineral sources, but identifying the primary source or changes of sources in response to environmental change is problematic. We used strontium isotope and alkaline earth element concentration ratios of trees and soils to discern the record of Ca sources for red spruce at a base-poor, acid deposition-impacted watershed. We measured 87 Sr/ 86 Sr and chemical compositions of cross-sectional stemwood cores of red spruce, other spruce tissues and sequential extracts of co-located soil samples. 87 Sr/ 86 Sr and Sr/Ba ratios together provide a tracer of alkaline earth element sources that distinguishes the plant-available fraction of the shallow organic soils from those of deeper organic and mineral soils. Ca/Sr ratios proved less diagnostic, due to within-tree processes that fractionate these elements from each other. Over the growth period from 1870 to 1960, 87 Sr/ 86 Sr and Sr/Ba ratios of stemwood samples became progressively more variable and on average trended toward values that considered together are characteristic of the uppermost forest floor. In detail the stemwood chemistry revealed an episode of simultaneous enhanced uptake of all alkaline earth elements during the growth period from 1930 to 1960, coincident with reported local and regional increases in atmospheric inputs of inorganic acidity. We attribute the temporal trends in stemwood chemistry to progressive shallowing of the effective depth of alkaline earth element uptake by fine roots over this growth period, due to preferential concentration of fine roots in the upper forest floor coupled with reduced nutrient uptake by roots in the lower organic and upper mineral soils in response to acid-induced aluminum toxicity. Although both increased atmospheric deposition and selective weathering of Ca-rich minerals such as apatite provide possible alternative explanations of aspects of the observed trends, the chemical buffering capacity of the forest floor-biomass pool limits their effectiveness as causal mechanisms.

Biogeochemistry↗

Modeling the height of young forests regenerating from recent disturbances in Mississippi using Landsat and ICESat data

Many forestry and earth science applications require spatially detailed forest height data sets. Among the various remote sensing technologies, lidar offers the most potential for obtaining reliable height measurement. However, existing and planned spaceborne lidar systems do not have the capability to produce spatially contiguous, fine resolution forest height maps over large areas. This paper describes a Landsat–lidar fusion approach for modeling the height of young forests by integrating historical Landsat observations with lidar data acquired by the Geoscience Laser Altimeter System (GLAS) instrument onboard the Ice, Cloud, and land Elevation (ICESat) satellite. In this approach, “young” forests refer to forests reestablished following recent disturbances mapped using Landsat time-series stacks (LTSS) and a vegetation change tracker (VCT) algorithm. The GLAS lidar data is used to retrieve forest height at sample locations represented by the footprints of the lidar data. These samples are used to establish relationships between lidar-based forest height measurements and LTSS–VCT disturbance products. The height of “young” forest is then mapped based on the derived relationships and the LTSS–VCT disturbance products. This approach was developed and tested over the state of Mississippi. Of the various models evaluated, a regression tree model predicting forest height from age since disturbance and three cumulative indices produced by the LTSS–VCT method yielded the lowest cross validation error. The R 2 and root mean square difference (RMSD) between predicted and GLAS-based height measurements were 0.91 and 1.97 m, respectively. Predictions of this model had much higher errors than indicated by cross validation analysis when evaluated using field plot data collected through the Forest Inventory and Analysis Program of USDA Forest Service. Much of these errors were due to a lack of separation between stand clearing and non-stand clearing disturbances in current LTSS–VCT products and difficulty in deriving reliable forest height measurements using GLAS samples when terrain relief was present within their footprints. In addition, a systematic underestimation of about 5 m by the developed model was also observed, half of which could be explained by forest growth that occurred between field measurement year and model target year. The remaining difference suggests that tree height measurements derived using waveform lidar data could be significantly underestimated, especially for young pine forests. Options for improving the height modeling approach developed in this study were discussed.

Remote Sensing of Environment↗

Reconstructing the velocity and deformation of a rapid landslide using multiview video

Noncontact measurements of spatially varied ground surface deformation during landslide motion can provide important constraints on landslide mechanics. Here, we present and test a new method for extracting measurements of rapid landslide surface displacement and velocity (accelerations of approximately 1 m/s 2 ) using sequences of stereo images obtained from a pair of inexpensive, stationary 4K video cameras with nominal frame rates of 29.97 Hz. The method combines elements of Structure from Motion with those of optical flow to extract data on 3‐D evolution of the ground surface during slope failure. We apply the method to an experiment at the U.S. Geological Survey debris‐flow flume in which a high‐speed, liquefying landslide was triggered by gradually adding water to a 6‐m 3 prism of loosely packed sediment on a 31° slope. Strip‐scanning lidar measurements made during the experiment corroborate our video‐based measurements, but the latter encompassed the entire landslide surface and were much lower in cost. Our video‐based measurements enabled computation of depth‐integrated landslide dilation/contraction rates. The range of computed rates was within the ranges inferred from independent measurements of evolving pore water pressures and reasonable estimates of the hydraulic permeability of the sediment. Dilation and contraction rates play a crucial role in landslide mechanics. The dilation and contraction we observe contradict the incompressible flow assumption used in many studies that have employed noncontact methods to infer landslide properties.

Oregon↗

Interpreting and reporting 40Ar/39Ar geochronologic data

The 40 Ar/ 39 Ar dating method is among the most versatile of geochronometers, having the potential to date a broad variety of K-bearing materials spanning from the time of Earth’s formation into the historical realm. Measurements using modern noble-gas mass spectrometers are now producing 40 Ar/ 39 Ar dates with analytical uncertainties of ∼0.1%, thereby providing precise time constraints for a wide range of geologic and extraterrestrial processes. Analyses of increasingly smaller subsamples have revealed age dispersion in many materials, including some minerals used as neutron fluence monitors. Accordingly, interpretive strategies are evolving to address observed dispersion in dates from a single sample. Moreover, inferring a geologically meaningful “age” from a measured “date” or set of dates is dependent on the geological problem being addressed and the salient assumptions associated with each set of data. We highlight requirements for collateral information that will better constrain the interpretation of 40 Ar/ 39 Ar data sets, including those associated with single-crystal fusion analyses, incremental heating experiments, and in situ analyses of microsampled domains. To ensure the utility and viability of published results, we emphasize previous recommendations for reporting 40 Ar/ 39 Ar data and the related essential metadata, with the amendment that data conform to evolving standards of being findable, accessible, interoperable, and reusable (FAIR) by both humans and computers. Our examples provide guidance for the presentation and interpretation of 40 Ar/ 39 Ar dates to maximize their interdisciplinary usage, reproducibility, and longevity.

GSA Bulletin↗

Sedimentary response to orogenic exhumation in the northern Rocky Mountain Basin and Range province, Flint Creek basin, west-central Montana

Middle Eocene through Upper Miocene sedimentary and volcanic rocks of the Flint Creek basin in western Montana accumulated during a period of significant paleoclimatic change and extension across the northern Rocky Mountain Basin and Range province. Gravity modelling, borehole data, and geologic mapping from the Flint Creek basin indicate that subsidence was focused along an extensionally reactivated Sevier thrust fault, which accommodated up to 800 m of basin fill while relaying stress between the dextral transtensional Lewis and Clark lineament to the north and the Anaconda core complex to the south. Northwesterly paleocurrent indicators, foliated metamorphic lithics, 64 Ma (40Ar/39Ar) muscovite grains, and 76 Ma (U-Pb) zircons in a ca. 27 Ma arkosic sandstone are consistent with Oligocene exhumation and erosion of the Anaconda core complex. The core complex and volcanic and magmatic rocks in its hangingwall created an important drainage divide during the Paleogene shedding detritus to the NNW and ESE. Following a major period of Early Miocene tectonism and erosion, regional drainage networks were reorganized such that paleoflow in the Flint Creek basin flowed east into an internally drained saline lake system. Renewed tectonism during Middle to Late Miocene time reestablished a west-directed drainage that is recorded by fluvial strata within a Late Miocene paleovalley. These tectonic reorganizations and associated drainage divide explain observed discrepancies in provenance studies across the province. Regional correlation of unconformities and lithofacies mapping in the Flint Creek basin suggest that localized tectonism and relative base level fluctuations controlled lithostratigraphic architecture.

Canadian Journal of Earth Sciences↗

Rayleigh wave amplitude uncertainty across the Global Seismographic Network and potential implications for global tomography

The Global Seismographic Network (GSN) is a multiuse, globally distributed seismic network used by seismologists, to both characterize earthquakes and study the Earth’s interior. Most stations in the network have two collocated broadband seismometers, which enable network operators to identify potential metadata and sensor issues. In this study, we investigate the accuracy with which surface waves can be measured across the GSN, by comparing waveforms of vertical‐component Rayleigh waves from M w "> M w Mw 6 and larger events between collocated sensor pairs. We calculate both the amplitude deviation and correlation coefficient between waveforms at sensor pairs. In total, we make measurements on over 670,000 event–station pairs from events that occurred from 1 January 2010 to 1 January 2020. We find that the average sensor‐pair amplitude deviation, and, therefore, GSN calibration level, is, approximately, 4% in the 25–250 s period band. Although, we find little difference in sensor‐pair amplitude deviations as a function of period across the entire network, the amount of useable data decreases rapidly as a function of increasing period. For instance, we determined that just over 12% of records at 250 s period provided useable recordings (e.g., sensor‐pair amplitude deviations of less than 20% and sensor‐pair correlation greater than 0.95). We then use these amplitude‐estimate deviations to identify how data coverage and quality could be limiting our ability to invert for whole Earth 3D attenuation models. We find an increase in the variance of our attenuation models with increasing period. For example, our degree 12 attenuation inversion at 250 s period shows 32% more variance than our degree 12 attenuation model at 25 s. This indicates that discrepancies of deep‐mantle tomography between studies could be the result of these large uncertainties. Because these high uncertainties arise from limited, high‐quality observations of long‐period ( ⁠ > 100    s ">

Bulletin of the Seismological Society of America↗

Teleseismic waves reveal anisotropic poroelastic response of wastewater disposal reservoir

Connecting earthquake nucleation in basement rock to fluid injection in basal, sedimentary reservoirs, depends heavily on choices related to the poroelastic properties of the fluid-rock system, thermo-chemical effects notwithstanding. Direct constraints on these parameters outside of laboratory settings are rare, and it is commonly assumed that the rock layers are isotropic. With the Arbuckle wastewater disposal reservoir in Osage County, Oklahoma, high-frequency formation pressure changes and collocated broadband ground velocities measured during the passing of large teleseismic waves show a poroelastic response of the reservoir that is both azimuthally variable and anisotropic; this includes evidence of static shifts in pressure that presumably relate to changes in local permeability. The azimuthal dependence in both the static response and shear coupling appears related to tectonic stress and strain indicators such as the orientations of the maximum horizontal stress and faults and fractures. Using dynamic strains from a nearby borehole strainmeter, we show that the ratio of shear to volumetric strain coupling is ∼ 0.41 "> ∼ 0.41 ∼0.41 which implies a mean Skempton's coefficient of A = 0.24 "> A = 0.24 A=0.24 over the plausible range of the undrained Poisson's ratio. Since these observations are made at relatively low confining pressure and differential stress, we suggest that the hydraulically conductive fracture network is a primary control on the coupling between pore pressure diffusion and elastic stresses in response to natural or anthropogenic sources.

Earth and Planetary Physics↗

New insights into Kilauea's volcano dynamics brought by large-scale relative relocation of microearthquakes

We investigated the microseismicity recorded in an active volcano to infer information concerning the volcano structure and long-term dynamics, by using relative relocations and focal mechanisms of microearthquakes. There were 32,000 earthquakes of the Mauna Loa and Kilauea volcanoes recorded by more than eight stations of the Hawaiian Volcano Observatory seismic network between 1988 and 1999. We studied 17,000 of these events and relocated more than 70%, with an accuracy ranging from 10 to 500 m. About 75% of these relocated events are located in the vicinity of subhorizontal decollement planes, at a depth of 8-11 km. However, the striking features revealed by these relocation results are steep southeast dipping fault planes working as reverse faults, clearly located below the decollement plane and which intersect it. If this decollement plane coincides with the pre-Mauna Loa seafloor, as hypothesized by numerous authors, such reverse faults rupture the pre-Mauna Loa oceanic crust. The weight of the volcano and pressure in the magma storage system are possible causes of these ruptures, fully compatible with the local stress tensor computed by Gillard et al. [1996]. Reverse faults are suspected of producing scarps revealed by kilometer-long horizontal slip-perpendicular lineations along the decollement surface and therefore large-scale roughness, asperities, and normal stress variations. These are capable of generating stick-slip, large-magnitude earthquakes, the spatial microseismic pattern observed in the south flank of Kilauea volcano, and Hilina-type instabilities. Rupture intersecting the decollement surface, causing its large-scale roughness, may be an important parameter controlling the growth of Hawaiian volcanoes.

Journal of Geophysical Research B: Solid Earth↗

Rise of a variable-viscosity fluid in a steadily spreading wedge-shaped conduit with accreting walls

Relatively rigid plates making up the outer 50 to 100 km of the Earth are steadily separating from one another along narrow globe-circling zones of submarine volcanism, the oceanic spreading centers. Continuity requires that the viscous underlying material rise beneath spreading centers and accrete onto the steadily diverging plates. It is likely that during the rise the viscosity changes systematically and that the viscous tractions exerted on the plates contribute to the unique pattern of submarine mountains and earthquake faults observed at spreading centers. The process is modeled by viscous creep in a wedge-shaped conduit (with apex at the sea floor) in which the viscosity varies as r m where r is distance from the apex and m is a parameter. For these conditions, the governing differential equations take a simple form. The solution for the velocity is independent of r and of the sign of m . As viscous stresses vary as r m-1 , the pattern of stress on the conduit wall is sensitive to viscosity variation. For negative m , the viscous pressure along the base of the conduit is quite uniform; for positive m , it falls toward zero in the axial region as the conduit base widens. For small opening angles, viscous forces push the plates apart, and for large ones, they oppose plate separation. Though highly idealized, the solution provides a tool for investigating tectonic processes at spreading centers.

Journal of Research of the U.S. Geological Survey↗