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 847 records · Page 47Linked to original sources
USGS1997

Uplift of the Transantarctic Mountains and the bedrock beneath the East Antarctic ice sheet

In recent years the Transantarctic Mountains (TAM), the largest noncontractional mountain belt in the world, have become the focus of modelers who explained their uplift by a variety of isostatic and thermal mechanisms. A problem with these models is a lack of available data to compare with model predictions. We report here the results of a 312-km-long geophysical traverse conducted in 1993/1994 in the hinterland of the TAM. Using detailed subglacial topography and gravity measurements, we confirm the origin of the TAM as a flexural uplift of the edge of East Antarctica. Using an elastic model with a free edge, we can jointly fit the topography and the gravity with a plate having an elastic thickness of 85±15 km and a preuplift elevation of 700±50 m for East Antarctica. Using a variety of evidence, we argue that the uplift is coincident with a relatively minor tectonic event of transtensional motion between East and West Antarctica during the Eocene rather than the Late Cretaceous rifting event that created the Ross Embayment. We suggest that this transtensional motion caused the continuous plate to break, which created an escarpment that significantly increased the rates of erosion and exhumation. Results from the geophysical traverse also extend our knowledge of the bedrock geology from the exposures within the TAM to the ice covered interior. Our interpretation suggests that the Ferrar flood basalts extend at least 100 km westward under the ice. The Beacon Supergroup of Paleozoic and Mesozoic sediments thins gradually under the ice and its reconstructed thickness is reminiscent of profiles of foreland basins. Finally, there is no indication in the gravity field for an incomplete rebound due to significant melting of the East Antarctic ice sheet since the last glacial period.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

Stressing of the New Madrid seismic zone by a lower crust detachment fault

A new mechanical model for the cause of the New Madrid seismic zone in the central United States is analyzed. The model contains a subhorizontal detachment fault which is assumed to be near the domed top surface of locally thickened anomalous lower crust (“rift pillow”). Regional horizontal compression induces slip on the fault, and the slip creates a stress concentration in the upper crust above the rift pillow dome. In the coseismic stage of the model earthquake cycle, where the three largest magnitude 7–8 earthquakes in 1811–1812 are represented by a single model mainshock on a vertical northeast trending fault, the model mainshock has a moment equivalent to a magnitude 8 event. During the interseismic stage, corresponding to the present time, slip on the detachment fault exerts a right-lateral shear stress on the locked vertical fault whose failure produces the model mainshock. The sense of shear is generally consistent with the overall sense of slip of 1811–1812 and later earthquakes. Predicted rates of horizontal strain at the ground surface are about 10 −7 year −1 and are comparable to some observed rates. The model implies that rift pillow geometry is a significant influence on the maximum possible earthquake magnitude.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Paleomagnetism, paleogeographic origins, and uplift history of the Coast Range ophiolite at Mount Diablo, California

Divergent paleogeographic origins have been proposed for the Coast Range ophiolite of western California which are testable using paleomagnetic methods. Paleomagnetic data for Middle Jurassic pillow lavas and diabase sills of the Coast Range ophiolite at Mount Diablo, northern California, indicate that they contain two components of remanent magnetization. The characteristic directions have normal and reversed polarities and apparently are carried by Ti-poor magnetite. This magnetization is inferred to have been acquired during emplacement and seafloor alteration at an ancient spreading ridge. The paleolatitude calculated from its structurally corrected mean direction is 20°N ± 9° and agrees with the expected direction for stable North America; this result is also consistent with the concordant paleolatitude (32°N ± 8°) recently determined for Upper Jurassic Coast Range ophiolite at Stanley Mountain in southern California. In addition, clockwise vertical axis rotation of Mount Diablo (143° ± 11°) is indicated by the characteristic magnetization direction. An overprint component is inferred to have been acquired during uplift of Mount Diablo since the Miocene.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

A generalized law for brittle deformation of Westerly granite

A semiempirical constitutive law is presented for the brittle deformation of intact Westerly granite. The law can be extended to larger displacements, dominated by localized deformation, by including a displacement-weakening break-down region terminating in a frictional sliding regime often described by a rate- and state-dependent constitutive law. The intact deformation law, based on an Arrhenius type rate equation, relates inelastic strain rate to confining pressure P c , differential stress σ Δ , inelastic strain ε i and temperature T . The basic form of the law for deformation prior to fault nucleation is where σ o and ε o are normalization constants (dependent on confining pressure), a is rate sensitivity of stress, and α is a shape parameter. At room temperature, eight experimentally determined coefficients are needed to fully describe the stress-strain-strain rate response for Westerly granite from initial loading to failure. Temperature dependence requires apparent activation energy ( E * ∼ 90 kJ/mol) and one additional experimentally determined coefficient. The similarity between the prefailure constitutive law for intact rock and the rate- and state-dependent friction laws for frictional sliding on fracture surfaces suggests a close connection between these brittle phenomena.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Detrital zircon U-Pb geochronology of Cambrian to Triassic miogeoclinal and eugeoclinal strata of Sonora, Mexico

One hundred and eighty two individual detrital zircon grains from Cambrian through Permian miogeoclinal strata, Ordovician eugeoclinal rocks, and Triassic post-orogenic sediments in northwestern Sonora have been analyzed. During Cambrian, Devonian, Permian, and Triassic time, most zircons accumulating along this part of the Cordilleran margin were shed from 1.40–1.45 and 1.62–1.78 Ga igneous rocks that are widespread in the southwestern United States and northwestern Mexico. Zircons with ages of approximately 1.11 Ga are common in Cambrian strata and were apparently shed from granite bodies near the sample site. The sources of 225–280 Ma zircons in our Triassic sample are more problematic, as few igneous rocks of these ages are recognized in northwestern Mexico. Such sources may be present but unrecognized, or the grains could have been derived from igneous rocks of the appropriate ages to the northwest in the Mojave Desert region, to the east in Chihuahua and Coahuila, or to the south in accreted(?) arc-type terranes. Because the zircon grains in our Cambrian and Devonian to Triassic samples could have accumulated in proximity to basement rocks near their present position or in the Death Valley region of southern California, our data do not support or refute the existence of the Mojave-Sonora megashear. Ordovician strata of both miogeoclinal and eugeoclinal affinity are dominated by >1.77 Ga detrital zircons, which are considerably older than most basement rocks in the region. Zircon grains in the miogeoclinal sample were apparently derived from the Peace River arch area of northwestern Canada and transported southward by longshore currents. The eugeoclinal grains may also have come from the Peace River arch region, with southward transport by either sedimentary or tectonic processes, or they may have been shed from off-shelf slivers of continents (perhaps Antarctica?) removed from the Cordilleran margin during Neoproterozoic rifting. It is also possible that the Ordovician eugeoclinal strata are far traveled and exotic to North America.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Observations on basaltic lava streams in tubes from Kilauea Volcano, island of Hawai'i

From 1986 to 1997, the Pu'u 'O'o-Kupaianaha eruption of Kilauea produced a vast pahoehoe flow field fed by lava tubes that extended 10–12 km from vents on the volcano's east rift zone to the ocean. Within a kilometer of the vent, tubes were as much as 20 m high and 10–25 m wide. On steep slopes (4–10°) a little farther away from the vent, some tubes formed by roofing over of lava channels. Lava streams were typically 1–2 m deep flowing within a tube that here was typically 5 m high and 3 m wide. On the coastal plain (<1°), tubes within inflated sheet flows were completely filled, typically 1–2 m high, and several tens of meters wide. Tubes develop as a flow's crust grows on the top, bottom, and sides of the tubes, restricting the size of the fluid core. The tubes start out with nearly elliptical cross-sectional shapes, many times wider than high. Broad, flat sheet flows evolve into elongate tumuli with an axial crack as the flanks of the original flow were progressively buried by breakouts. Temperature measurements and the presence of stalactites in active tubes confirmed that the tube walls were above the solidus and subject to melting. Sometimes, the tubes began downcutting. Progressive downcutting was frequently observed through skylights; a rate of 10 cm/d was measured at one skylight for nearly 2 months.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Viscoelastic coupling model of the San Andreas fault along the Big Bend, southern California

The big bend segment of the San Andreas fault is the 300-km-long segment in southern California that strikes about N65°W, roughly 25° counterclockwise from the local tangent to the small circle about the Pacific-North America pole of rotation. The broad distribution of deformation of trilateration networks along this segment implies a locking depth of at least 25 km as interpreted by the conventional model of strain accumulation (continuous slip on the fault below the locking depth at the rate of relative plate motion), whereas the observed seismicity and laboratory data on fault strength suggest that the locking depth should be no greater than 10 to 15 km. The discrepancy is explained by the viscoelastic coupling model which accounts for the viscoelastic response of the lower crust. Thus the broad distribution of deformation observed across the big bend segment can be largely associated with the San Andreas fault itself, not subsidiary faults distributed throughout the region. The Working Group on California Earthquake Probabilities [1995] in using geodetic data to estimate the seismic risk in southern California has assumed that strain accumulated off the San Andreas fault is released by earthquakes located off the San Andreas fault. Thus they count the San Andreas contribution to total seismic moment accumulation more than once, leading to an overestimate of the seismicity for magnitude 6 and greater earthquakes in their Type C zones.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Quantitative measure of the variation in fault rheology due to fluid-rock interactions

We analyze friction data from two published suites of laboratory tests on granite in order to explore and quantify the effects of temperature ( T ) and pore water pressure ( P p ) on the sliding behavior of faults. Rate-stepping sliding tests were performed on laboratory faults in granite containing “gouge” (granite powder), both dry at 23° to 845°C [ Lockner et al. , 1986], and wet ( P p = 100 MPa) at 23° to 600°C [ Blanpied et al ., 1991, 1995]. Imposed slip velocities ( V ) ranged from 0.01 to 5.5 μm/s, and effective normal stresses were near 400 MPa. For dried granite at all temperatures, and wet granite below ∼300°C, the coefficient of friction (μ) shows low sensitivity to V , T , and P p . For wet granite above ∼350°, μ drops rapidly with increasing T and shows a strong, positive rate dependence and protracted strength transients following steps in V , presumably reflecting the activity of a water-aided deformation process. By inverting strength data from velocity stepping tests we determined values for parameters in three formulations of a rate- and state-dependent constitutive law. One or two state variables were used to represent slip history effects. Each velocity step yielded an independent set of values for the nominal friction level, five constitutive parameters (transient parameters a , b 1 , and b 2 and characteristic displacements D c1 and D c2 ), and the velocity dependence of steady state friction ∂μ ss /∂ ln V = a - b 1 − b 2 . Below 250°, data from dry and most wet tests are adequately modeled by using the “slip law” [ Ruina , 1983] and one state variable ( a = 0.003 to 0.018, b = 0.001 to +0.018, D c ≈ 1 to 20 μm). Dried tests above 250° can also be fitted with one state variable. In contrast, wet tests above 350° require higher direct rate dependence ( a = 0.03 to 0.12), plus a second state variable with large, negative amplitude ( b 2 = −0.03 to −0.14) and large characteristic displacement ( D c2 = 300 to >4000 μm). Thus the parameters a , b 1 , and b 2 for wet granite show a pronounced change in their temperature dependence in the range 270° to 350°C, which may reflect a change in underlying deformation mechanism. We quantify the trends in parameter values from 25° to 600°C by piecewise linear regressions, which provide a straightforward means to incorporate the full constitutive response of granite into numerical models of fault slip. The modeling results suggest that the succeptibility for unstable (stick-slip) sliding is maximized between 90° and 360°C, in agreement with laboratory observations and consistent with the depth range of earthquakes on mature faults in the continental crust.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

The excitation and characteristic frequency of the long-period volcanic event: An approach based on an inhomogeneous autoregressive model of a linear dynamic system

We present a method to quantify the source excitation function and characteristic frequencies of long-period volcanic events. The method is based on an inhomogeneous autoregressive (AR) model of a linear dynamic system, in which the excitation is assumed to be a time-localized function applied at the beginning of the event. The tail of an exponentially decaying harmonic waveform is used to determine the characteristic complex frequencies of the event by the Sompi method. The excitation function is then derived by operating an AR filter constructed from the characteristic frequencies to the entire seismogram of the event, including the inhomogeneous part of the signal. We apply this method to three long-period events at Kusatsu-Shirane Volcano, central Japan, whose waveforms display simple decaying monochromatic oscillations except for the beginning of the events. We recover time-localized excitation functions lasting roughly 1 s at the start of each event and find that the estimated functions are very similar to each other at all the stations of the seismic network for each event. The phases of the characteristic oscillations referred to the estimated excitation function fall within a narrow range for almost all the stations. These results strongly suggest that the excitation and mode of oscillation are both dominated by volumetric change components. Each excitation function starts with a pronounced dilatation consistent with a sudden deflation of the volumetric source which may be interpreted in terms of a choked-flow transport mechanism. The frequency and Q of the characteristic oscillation both display a temporal evolution from event to event. Assuming a crack filled with bubbly water as seismic source for these events, we apply the Van Wijngaarden-Papanicolaou model to estimate the acoustic properties of the bubbly liquid and find that the observed changes in the frequencies and Q are consistently explained by a temporal change in the radii of the bubbles characterizing the bubbly water in the crack.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Sulfur in serpentinized oceanic peridotites: Serpentinization processes and microbial sulfate reduction

The mineralogy, contents, and isotopic compositions of sulfur in oceanic serpentinites reflect variations in temperatures and fluid fluxes. Serpentinization of <1 Ma peridotites at Hess Deep occurred at high temperatures (200°–400°C) and low water/rock ratios. Oxidation of ferrous iron to magnetite maintained low ƒO 2 and produced a reduced, low-sulfur assemblage including NiFe alloy. Small amounts of sulfate reduction by thermophilic microbes occurred as the system cooled, producing low-δ 34 S sulfide (1.5‰ to −23.7‰). In contrast, serpentinization of Iberian Margin peridotites occurred at low temperatures(∼20°–200°C) and high water/rock ratios. Complete serpentinization and consumption of ferrous iron allowed evolution to higher ƒO 2 . Microbial reduction of seawater sulfate resulted in addition of low-δ 34 S sulfide (∼15 to ∼43‰) and formation of higher-sulfur assemblages that include valleriite and pyrite. The high SO 4 /total S ratio of Hess Deep serpentinites (0.89) results in an increase of total sulfur and high δ 34 S of total sulfur (mean ∼8‰). In contrast, Iberian Margin serpentinites gained large amounts of 34 S-poor sulfide (mean total S = 3800 ppm), and the high sulfide/total S ratio (0.61) results in a net decrease in δ 34 S of total sulfur (mean ≈ −5‰). Thus serpentinization is a net sink for seawater sulfur, but the amount fixed and its isotopic composition vary significantly. Serpentinization may result in uptake of 0.4–14 × 10 12 g S yr −1 from the oceans, comparable to isotopic exchange in mafic rocks of seafloor hydrothermal systems and approaching global fluxes of riverine sulfate input and sedimentary sulfide output.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Alteration and mineralization of an oceanic forearc and the ophiolite-ocean crust analogy

Mineralogical, chemical, and isotopic (O, C, S, and Sr) analyses were performed on minerals and bulk rocks from a forearc basement section to understand alteration processes and compare with mid-ocean ridges (MOR) and ophiolites. Ocean Drilling Program Hole 786B in the Izu-Bonin forearc penetrates 103 m of sediment and 725 m into volcanic flows, breccias, and basal dikes. The rocks comprise boninites and andesites to rhyolites. Most of the section was affected by low-temperature (<100°C) seawater alteration, with temperatures increasing downward. The rocks are partly (5–25%) altered to smectite, Fe-oxyhydroxide, calcite, and phillipsite, and exhibit gains of K, Rb, and P, loss of Ca, variable changes in Si, Na, Mg, Fe, Sr, and Y, and elevated δ 18 O and 87 Sr/ 86 Sr. Higher temperatures (∼150°C) in the basal dikes below 750 m led to more intense alteration and formation of chlorite-smectite, corrensite, albite, K-feldspar, and quartz (±chlorite). A 5 m thick hydrothermally altered and pyritized zone at 815 m in the basal dikes reacted with mixtures of seawater and hydrothermal fluids to Mg-chlorite, albite, and pyrite, and gained Mg and S and lost Si and Ca. Focused flow of hydrothermal fluids produced serialization halos (Na-K sericite, quartz, pyrophyllite, K-feldspar, and pyrite) along quartz veins at temperatures of 200°–250°C. High 87 Sr/ 86 Sr ratios of chloritized (∼0.7055) and sericitized (∼0.7065) rocks indicate involvement of seawater via mixing with hydrothermal fluids. Low δ 34 S of sulfide (−2 to −5.5‰) and sulfate (12.5‰) are consistent with input of magmatic SO 2 into hydrothermal fluids and disproportionation to sulfide and sulfate. Alteration processes were generally similar to those at MORs, but the arc section is more intensively altered, in part because of the presence of abundant glassy rocks and mafic phases. The increase in alteration grade below 750 m and the mineralization in the basal dikes are analogous to changes that occur near the base of the volcanic section in MOR and the Troodos ophiolite.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Hydrothermal circulation at the Cleft-Vance overlapping spreading center: Results of a magnetometric resistivity survey

We report on a magnetometric resistivity sounding carried out in the overlapping spreading center between the Cleft and Vance segments of the Juan de Fuca Ridge. The data collected reveal a strong three dimensionality in the crustal electrical resistivity structure on wavelengths of a few kilometers. Areas of reduced crustal electrical resistivities, with values approaching that of seawater, are seen beneath the neovolcanic zones of both active spreading centers. We interpret these reduced resistivities as evidence of active hydrothermal circulation within the uppermost 1 km of hot, young oceanic crust.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Paleoseismic investigations in the Santa Cruz mountains, California: Implications for recurrence of large-magnitude earthquakes on the San Andreas fault

Trenching, microgeomorphic mapping, and tree ring analysis provide information on timing of paleoearthquakes and behavior of the San Andreas fault in the Santa Cruz mountains. At the Grizzly Flat site alluvial units dated at 1640–1659 A.D., 1679–1894 A.D., 1668–1893 A.D., and the present ground surface are displaced by a single event. This was the 1906 surface rupture. Combined trench dates and tree ring analysis suggest that the penultimate event occurred in the mid-1600 s, possibly in an interval as narrow as 1632–1659 A.D. There is no direct evidence in the trenches for the 1838 or 1865 earthquakes, which have been proposed as occurring on this part of the fault zone. In a minimum time of about 340 years only one large surface faulting event (1906) occurred at Grizzly Flat, in contrast to previous recurrence estimates of 95–110 years for the Santa Cruz mountains segment. Comparison with dates of the penultimate San Andreas earthquake at sites north of San Francisco suggests that the San Andreas fault between Point Arena and the Santa Cruz mountains may have failed either as a sequence of closely timed earthquakes on adjacent segments or as a single long rupture similar in length to the 1906 rupture around the mid-1600 s. The 1906 coseismic geodetic slip and the late Holocene geologic slip rate on the San Francisco peninsula and southward are about 50–70% and 70% of their values north of San Francisco, respectively. The slip gradient along the 1906 rupture section of the San Andreas reflects partitioning of plate boundary slip onto the San Gregorio, Sargent, and other faults south of the Golden Gate. If a mid-1600 s event ruptured the same section of the fault that failed in 1906, it supports the concept that long strike-slip faults can contain master rupture segments that repeat in both length and slip distribution. Recognition of a persistent slip rate gradient along the northern San Andreas fault and the concept of a master segment remove the requirement that lower slip sections of large events such as 1906 must fill in on a periodic basis with smaller and more frequent earthquakes.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Stress transferred by the 1995 Mw = 6.9 Kobe, Japan, shock: Effect on aftershocks and future earthquake probabilities

The Kobe earthquake struck at the edge of the densely populated Osaka-Kyoto corridor in southwest Japan. We investigate how the earthquake transferred stress to nearby faults, altering their proximity to failure and thus changing earthquake probabilities. We find that relative to the pre-Kobe seismicity, Kobe aftershocks were concentrated in regions of calculated Coulomb stress increase and less common in regions of stress decrease. We quantify this relationship by forming the spatial correlation between the seismicity rate change and the Coulomb stress change. The correlation is significant for stress changes greater than 0.2–1.0 bars (0.02–0.1 MPa), and the nonlinear dependence of seismicity rate change on stress change is compatible with a state- and rate-dependent formulation for earthquake occurrence. We extend this analysis to future mainshocks by resolving the stress changes on major faults within 100 km of Kobe and calculating the change in probability caused by these stress changes. Transient effects of the stress changes are incorporated by the state-dependent constitutive relation, which amplifies the permanent stress changes during the aftershock period. Earthquake probability framed in this manner is highly time-dependent, much more so than is assumed in current practice. Because the probabilities depend on several poorly known parameters of the major faults, we estimate uncertainties of the probabilities by Monte Carlo simulation. This enables us to include uncertainties on the elapsed time since the last earthquake, the repeat time and its variability, and the period of aftershock decay. We estimate that a calculated 3-bar (0.3-MPa) stress increase on the eastern section of the Arima-Takatsuki Tectonic Line (ATTL) near Kyoto causes fivefold increase in the 30-year probability of a subsequent large earthquake near Kyoto; a 2-bar (0.2-MPa) stress decrease on the western section of the ATTL results in a reduction in probability by a factor of 140 to 2000. The probability of a M w = 6.9 earthquake within 50 km of Osaka during 1997–2007 is estimated to have risen from 5–6% before the Kobe earthquake to 7–11% afterward; during 1997–2027, it is estimated to have risen from 14–16% before Kobe to 16–22%.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Two-way coupling between Vesuvius eruptions and southern Apennine earthquakes, Italy, by elastic stress transfer

During the past 1000 years, eruptions of Vesuvius have often been accompanied by large earthquakes in the Apennines 50–60 km to the northeast. Statistical investigations had shown that earthquakes often preceded eruptions, typically by less than a decade, but did not provide a physical explanation for the correlation. Here, we explore elastic stress interaction between earthquakes and eruptions under the hypothesis that small stress changes can promote events when the Apennine normal faults and the Vesuvius magma body are close to failure. We show that earthquakes can promote eruptions by compressing the magma body at depth and opening suitably oriented near-surface conduits. Voiding the magma body in turns brings these same normal faults closer to Coulomb failure, promoting earthquakes. Such a coupling is strongest if the magma reservoir is a dike oriented normal to the regional extension axis, parallel to the Apennines, and the near-surface conduits and fissures are oriented normal to the Apennines. This preferred orientation suggests that the eruptions issuing from such fissures should be most closely linked in time to Apennine earthquakes. Large Apennine earthquakes since 1400 are calculated to have transferred more stress to Vesuvius than all but the largest eruptions have transferred to Apennine faults, which may explain why earthquakes more commonly lead than follow eruptions. A two-way coupling may thus link earthquakes and Vesuvius eruptions along a 100-km-long set of faults. We test the statistical significance of the earthquake-eruption correlation in the two-way coupling zone, and find a correlation significant at the 95% confidence level.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Waveform inversion of very long period impulsive signals associated with magmatic injection beneath Kilauea Volcano, Hawaii

We use data from broadband seismometers deployed around the summit of Kilauea Volcano to quantify the mechanism associated with a transient in the flow of magma feeding the east rift eruption of the volcano. The transient is marked by rapid inflation of the Kilauea summit peaking at 22 μrad 4.5 hours after the event onset, followed by slow deflation over a period of 3 days. Superimposed on the summit inflation is a series of sawtooth displacement pulses, each characterized by a sudden drop in amplitude lasting 5–10 s followed by an exponential recovery lasting 1–3 min. The sawtooth waveforms display almost identical shapes, suggesting a process involving the repeated activation of a fixed source. The particle motion associated with each sawtooth is almost linear, and its major swing shows compressional motion at all stations. Analyses of semblance and particle motion are consistent with a point source located 1 km beneath the northeast edge of the Halemaumau pit crater. To estimate the source mechanism, we apply a moment tensor inversion to the waveform data, assuming a point source embedded in a homogeneous half-space with compressional and shear wave velocities representative of the average medium properties at shallow depth under Kilauea. Synthetic waveforms are constructed by a superposition of impulse responses for six moment tensor components and three single force components. The origin times of individual impulses are distributed along the time axis at appropriately small, equal intervals, and their amplitudes are determined by least squares. In this inversion, the source time functions of the six tensor and three force components are determined simultaneously. We confirm the accuracy of the inversion method through a series of numerical tests. The results from the inversion show that the waveform data are well explained by a pulsating transport mechanism operating on a subhorizontal crack linking the summit reservoir to the east rift of Kilauea. The crack acts like a buffer in which a batch of fluid (magma and/or gas) accumulates over a period of 1–3 min before being rapidly injected into a larger reservoir (possibly the east rift) over a timescale of 5–10 s. The seismic moment and volume change associated with a typical batch of fluid are approximately 10 14 N m and 3000 m 3 , respectively. Our results also point to the existence of a single force component with amplitude of 10 9 N, which may be explained as the drag force generated by the flow of viscous magma through a narrow constriction in the flow path. The total volume of magma associated with the 4.5-hour-long activation of the pulsating source is roughly 500,000 m 3 in good agreement with the integrated volume flow rate of magma estimated near the eruptive site.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Earthquake triggering by transient and static deformations

Observational evidence for both static and transient near-field and far-field triggered seismicity are explained in terms of a frictional instability model, based on a single degree of freedom spring-slider system and rate- and state-dependent frictional constitutive equations. In this study a triggered earthquake is one whose failure time has been advanced by Δ t (clock advance) due to a stress perturbation. Triggering stress perturbations considered include square-wave transients and step functions, analogous to seismic waves and coseismic static stress changes, respectively. Perturbations are superimposed on a constant background stressing rate which represents the tectonic stressing rate. The normal stress is assumed to be constant. Approximate, closed-form solutions of the rate-and-state equations are derived for these triggering and background loads, building on the work of Dieterich [1992, 1994]. These solutions can be used to simulate the effects of static and transient stresses as a function of amplitude, onset time t 0 , and in the case of square waves, duration. The accuracies of the approximate closed-form solutions are also evaluated with respect to the full numerical solution and t 0 . The approximate solutions underpredict the full solutions, although the difference decreases as t 0 approaches the end of the earthquake cycle. The relationship between Δ t and t 0 differs for transient and static loads: a static stress step imposed late in the cycle causes less clock advance than an equal step imposed earlier, whereas a later applied transient causes greater clock advance than an equal one imposed earlier. For equal Δ t , transient amplitudes must be greater than static loads by factors of several tens to hundreds depending on t 0 . We show that the rate-and-state model requires that the total slip at failure is a constant, regardless of the loading history. Thus a static load applied early in the cycle, or a transient applied at any time, reduces the stress at the initiation of failure, whereas static loads that are applied sufficiently late raise it. Rate-and-state friction predictions differ markedly from those based on Coulomb failure stress changes (ΔCFS) in which Δ t equals the amplitude of the static stress change divided by the background stressing rate. The ΔCFS model assumes a stress failure threshold, while the rate-and-state equations require a slip failure threshold. The complete rate-and-state equations predict larger Δ t than the ΔCFS model does for static stress steps at small t 0 , and smaller Δ t than the ΔCFS model for stress steps at large t 0 . The ΔCFS model predicts nonzero Δ t only for transient loads that raise the stress to failure stress levels during the transient. In contrast, the rate-and-state model predicts nonzero Δ t for smaller loads, and triggered failure may occur well after the transient is finished. We consider heuristically the effects of triggering on a population of faults, as these effects might be evident in seismicity data. Triggering is manifest as an initial increase in seismicity rate that may be followed by a quiescence or by a return to the background rate. Available seismicity data are insufficient to discriminate whether triggered earthquakes are “new” or clock advanced. However, if triggering indeed results from advancing the failure time of inevitable earthquakes, then our modeling suggests that a quiescence always follows transient triggering and that the duration of increased seismicity also cannot exceed the duration of a triggering transient load. Quiescence follows static triggering only if the population of available faults is finite.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Paleomagnetism of the Middle Proterozoic Electra Lake Gabbro, Needle Mountains, southwestern Colorado

The Electra Lake Gabbro is a small 1.435 Ga pluton that intrudes 1.7 to 1.6 Ga gneisses and schists of the Needle Mountains in southwestern Colorado. Paleomagnetic samples were collected from the main phases of the gabbro, diabase dikes, granite, and alaskite dikes that cut the gabbro and from a partially melted zone in gneiss along the southern margin of the pluton. Gabbro, diabase, and some melt zone samples have a single-polarity characteristic magnetization of northeast declination ( D ) and moderate negative inclination ( I ). Demagnetization behavior and rock magnetic characteristics indicate that the remanence is carried by nearly pure magnetite. After correction for the minor west dip of overlying Paleozoic strata, we obtain a mean direction of D = 32.1°, I = −41.9° ( k : = 94, α 95 = 3.3°, N = 21 sites) and a paleomagnetic pole at 21.1°S, 221.1°E, ( K = 89, A 95 = 3.4°). This pole is similar to poles from the Middle Proterozoic Belt Supergroup but is located at a higher southerly latitude than poles from other 1.47–1.44 Ga plutons from North America, most of which plot at equatorial latitudes. The reason for this discrepancy is not clear but may result from a combination of factors, including unrecognized tilting of the gabbro, the failure of this relatively small pluton to fully average paleosecular variation, and uncertainties in the overall reliability of other 1.5–1.4 Ga poles of the North American apparent polar wander path.

Journal of Geophysical Research B: Solid Earth↗