Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysics Research”

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,387 records · Page 77Linked to original sources
USGS1997

Comparing nocturnal eddy covariance measurements to estimates of ecosystem respiration made by scaling chamber measurements at six coniferous boreal sites

During the growing season, nighttime ecosystem respiration emits 30–100% of the daytime net photosynthetic uptake of carbon, and therefore measurements of rates and understanding of its control by the environment are important for understanding net ecosystem exchange. Ecosystem respiration can be measured at night by eddy covariance methods, but the data may not be reliable because of low turbulence or other methodological problems. We used relationships between woody tissue, foliage, and soil respiration rates and temperature, with temperature records collected on site to estimate ecosystem respiration rates at six coniferous BOREAS sites at half-hour or 1-hour intervals, and then compared these estimates to nocturnal measurements of CO 2 exchange by eddy covariance. Soil surface respiration was the largest source of CO 2 at all sites (48–71%), and foliar respiration made a large contribution to ecosystem respiration at all sites (25–43%). Woody tissue respiration contributed only 5–15% to ecosystem respiration. We estimated error for the scaled chamber predictions of ecosystem respiration by using the uncertainty associated with each respiration parameter and respiring biomass value. There was substantial uncertainty in estimates of foliar and soil respiration because of the spatial variability of specific respiration rates. In addition, more attention needs to be paid to estimating foliar respiration during the early part of the growing season, when new foliage is growing, and to determining seasonal trends of soil surface respiration. Nocturnal eddy covariance measurements were poorly correlated to scaled chamber estimates of ecosystem respiration ( r 2 =0.06–0.27) and were consistently lower than scaled chamber predictions (by 27% on average for the six sites). The bias in eddy covariance estimates of ecosystem respiration will alter estimates of gross assimilation in the light and of net ecosystem exchange rates over extended periods.

Journal of Geophysical Research D: Atmospheres↗
USGS1997

A comparison of six methods for measuring soil-surface carbon dioxide fluxes

Measurements of soil-surface CO 2 fluxes are important for characterizing the carbon budget of boreal forests because these fluxes can be the second largest component of the budget. Several methods for measuring soil-surface CO 2 fluxes are available: (1) closed-dynamic-chamber systems, (2) closed-static-chamber systems, (3) open-chamber systems, and (4) eddy covariance systems. This paper presents a field comparison of six individual systems for measuring soil-surface CO 2 fluxes with each of the four basic system types represented. A single system is used as a reference and compared to each of the other systems individually in black spruce (Picea mariana), jack pine (Pinus banksiana), or aspen (Populus tremuloides) forests. Fluxes vary from 1 to 10 μmol CO 2 m −2 s −1 . Adjustment factors to bring all of the systems into agreement vary from 0.93 to 1.45 with an uncertainty of about 10–15%.

Journal of Geophysical Research D: Atmospheres↗
USGS1997

Dust deposition downwind of Owens (dry) Lake, 1991–1994: Preliminary findings

Salt‐rich dust derived from the Owens Lake playa is deposited in significant quantities to distances of at least 40 km north and south of the playa. Semiannual measurements from 1991 to 1994 of dust deposition rates (dust flux) and composition 2 m above the ground at seven sites in Owens Valley show that (1) dust in Owens Valley is derived mainly from the playa, although areas closer to the sites can also be sources; (2) south of the playa, dust flux is higher in the winter than in the summer, but north of the playa, dust flux is about the same or slightly lower in the winter; (3) on the playa, interannual variation in dust flux is large (factors of 5–10 during the 3 years), but at downwind sites, the variation is much smaller; (4) the dust typically has total soluble salt content as high as 30%, generally much higher than that of dust elsewhere in southern Nevada and California; and (5) to a distance of at least 40 km south of the playa, soluble salt flux is significantly higher than regional rates. The dust flux measurements indicate that significant quantities of salt‐rich dust are probably being added to the soils in the region around Owens Lake playa; these dust additions may affect soil pH and vegetation.

California↗
USGS1997

Accumulation and turnover of carbon in organic and mineral soils of the BOREAS northern study area

Rates of input, accumulation, and turnover of C differ markedly within soil profiles and in soils with different drainage in the BOREAS northern study area. Soil C storage increases from ∼3 kg C m −2 in well-drained, sandy soils to greater than 100 kg C m −2 in wetlands. Two modes of C accumulation were observed in upland soil profiles. Large annual C inputs (0.06–0.1 kg C m −2 yr −1 ) and slow decomposition (turnover times of 6–250 years) lead to rapid C accumulation in regrowing surface moss and detrital layers following fire. Deep organic layers that have accumulated over the millennia since the initiation of soil development, and are located below the most recent charred horizon, show slower rates of input (0.015–0.03 kg C m −2 yr −1 ) and turnover (100–1600 years) and accumulate C about 10 times slower than surface detrital layers. Rates of C input to soils derived from C and 14 C data were in accord with net primary production estimates, with highest rates of input (0.14–0.6 kg C m −2 yr −1 ) in wetlands. Turnover times for C in surface detrital layers were 6–15 years for well-drained sand soils that showed highest soil temperatures in summer, 30–40 years for wetlands, and 36–250 years for uplands with thick moss cover and black spruce trees. Long (>100 years) turnover times in upland black spruce/clay soils most likely reflect the influence of woody debris incorporated into detrital layers. Turnover times for deep organic and mineral layer C were controlled by drainage, with fastest turnover (80–130 years) in well-drained sand soils and slowest turnover (>3000 years) in wetlands. Total C accumulation rates, which account for C losses from both deep organic and surface detrital layers, are close to zero for sand/jack pine soils, 0.003–0.01 kg C m −2 yr −1 for moderately to poorly drained sites in mature forest stands, and 0.03 kg C m −2 yr −1 for a productive fen. Decomposition of organic matter more than several decades old accounts for 9–22% of total heterotrophic respiration at these sites. The rates of C accumulation derived here are decadal averages for specific stands and will vary as stands age or undergo disturbance. Extrapolation to larger regions and longer timescales, where burning offsets C gains in moss layers, will yield smaller rates of C storage.

Journal of Geophysical Research D: Atmospheres↗
USGS1997

Case study modeling of turbulent and mesoscale fluxes over the BOREAS region

Results from aircraft and surface observations provided evidence for the existence of mesoscale circulations over the Boreal Ecosystem-Atmosphere Study (BOREAS) domain. Using an integrated approach that included the use of analytical modeling, numerical modeling, and data analysis, we have found that there are substantial contributions to the total budgets of heat over the BOREAS domain generated by mesoscale circulations. This effect is largest when the synoptic flow is relatively weak, yet it is present under less favorable conditions, as shown by the case study presented here. While further analysis is warranted to document this effect, the existence of mesoscale flow is not surprising, since it is related to the presence of landscape patches, including lakes, which are of a size on the order of the local Rossby radius and which have spatial differences in maximum sensible heat flux of about 300 W m-2. We have also analyzed the vertical temperature profile simulated in our case study as well as high-resolution soundings and we have found vertical profiles of temperature change above the boundary layer height, which we attribute in part to mesoscale contributions. Our conclusion is that in regions with organized landscapes, such as BOREAS, even with relatively strong synoptic winds, dynamical scaling criteria should be used to assess whether mesoscale effects should be parameterized or explicitly resolved in numerical models of the atmosphere.

Journal of Geophysical Research D: Atmospheres↗
USGS1997

Martian drainage densities

Drainage densities on Mars range from zero over large areas of volcanic plains to 0.3–0.5 km −1 locally on some volcanoes. These values refer to geologic units, not to drainage basins, as is normal for terrestrial drainage densities. The highest values are close to the lowest terrestrial values derived by similar techniques. Drainage densities were determined for every geologic unit portrayed on the 1:15,000,000 geologic map of Mars. Except for volcanoes the geologic unit with the highest drainage density is the dissected Noachian plains with a drainage density of 0.0074 km −1 . The average drainage density for Noachian units is 0.0032 km −1 , for Hesperian units is 0.00047 km −1 , and for Amazonian units is 0.00007 km −1 , excluding the volcanoes. These values are 2–3 orders of magnitude lower than typical terrestrial densities as determined by similar techniques from Landsat images. The low drainage densities, despite a cumulative record that spans billions of years, indicate that compared with the Earth, the channel-forming processes have been very inefficient or have operated only rarely or that the surface is extremely permeable. The high drainage density on volcanoes is attributed to a local cause, such as hydrothermal activity, rather than to a global cause such as climate change.

Journal of Geophysical Research E: Planets↗
USGS1997

The Phanerozoic impact cratering rate: Evidence from the farside of the Moon

The relatively recent (<1 b.y.) flux of asteroids and comets forming large craters on the Earth and Moon may be accurately recorded by craters with bright rays on the Moon's farside. Many previously unknown farside rayed craters are clearly distinguished in the low-phase-angle images returned by the Clementine spacecraft. Some large rayed craters on the lunar nearside are probably significantly older than 1 Ga; rays remain visible over the maria due to compositional contrasts long after soils have reached optical maturity. Most of the farside crust has a more homogeneous composition and only immature rays are visible. The size-frequency distribution of farside rayed craters is similar to that measured for Eratosthenian craters (up to 3.2 b.y.) at diameters larger than 15 km. The areal density of farside rayed craters matches that of a corrected tabulation of nearside Copernican craters. Hence the presence of bright rays due to immature soils around large craters provides a consistent time-stratigraphic basis for defining the base of the Copernican System. The density of large craters less than ∼3.2 b.y. old is ∼3.2 times higher than that of large farside rayed craters alone. This observation can be interpreted in two ways: (1) the average cratering rate has been constant over the past 3.2 b.y. and the base of the Copernican is ∼1 Ga, or (2) the cratering rate has increased in recent geologic time and the base of the Copernican is less than 1 Ga. We favor the latter interpretation because the rays of Copernicus (800–850 m.y. old) appear to be very close to optical maturity, suggesting that the average Copernican cratering rate was ∼35% higher than the average Eratosthenian rate. Other lines of evidence for an increase in the Phanerozoic (545 Ga) cratering rate are (1) the densities of small craters superimposed on Copernicus and Apollo landing sites, (2) the rates estimated from well-dated terrestrial craters (≤120 m.y.) and from present-day astronomical observations, and (3) the Proterozoic rate suggested by the crater record of Australia. The hypothesis most consistent with several key observations is that the cratering rate has increased by ∼2x during the past ∼300 m.y‥

Journal of Geophysical Research E: Planets↗
USGS1997

Radar and photoclinometric studies of wrinkle ridges on Mars

Earth-based radar altimetry and image derived photoclinometric profiles were analyzed to examine both the long- and short-wavelength topography associated with wrinkle ridges on Mars. Photoclinometrically derived elevation data across wrinkle ridges were evaluated to determine the sensitivity of profiles to two empirical photoclinometric parameters, the horizontal digital number (HDN) and the scattered light value (SLV). The photoclinometric profiles are extremely sensitive to small variations in HDN. The sense of slope of a profile can be completely reversed over a range in HDN of as little as ±1. Comparably small variations in the SLV have relatively minor effects on the photoclinometrically derived elevations. The existence of elevation offsets from one side of the ridge to the other, reported in previous photoclinometric studies of martian wrinkles, were not confirmed through photoclinometry. In addition, no evidence of elevation offsets were found in Earth-based radar altimetry profiles across wrinkle ridges. In order to more accurately model wrinkle ridge topography, we controlled photoclinometrically derived elevations with long-wavelength topography obtained from the radar altimetry. The results of this study do not support kinematic models for the origin of planetary wrinkle ridges that involve deeply rooted thrust faults which separate crustal blocks at different elevations. A kinematic model involving buckling of shallow crustal layers into concentric folds that close, leading to the development of thrust faults, is consistent with wrinkle ridge morphology and terrestrial analogs. Recent geophysical studies of terrestrial analogs and the influence of shallow subsurface structures, particularly buried craters, on the localization of many wrinkle ridges on Mars suggest that thrust faults associated with the ridges are confined to the ridged plains material and do not extend into the lithosphere.

Journal of Geophysical Research E: Planets↗
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

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

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

California↗
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

Suppression of large earthquakes by stress shadows: A comparison of Coulomb and rate-and-state failure

Stress shadows generated by California's two most recent great earthquakes (1857 Fort Tejon and 1906 San Francisco) substantially modified 19th and 20th century earthquake history in the Los Angeles basin and in the San Francisco Bay area. Simple Coulomb failure calculations, which assume that earthquakes can be modeled as static dislocations in an elastic half-space, have done quite well at approximating how long the stress shadows, or relaxing effects, should last and at predicting where subsequent large earthquakes will not occur. There has, however, been at least one apparent exception to the predictions of such simple models. The 1911 M >6.0 earthquake near Morgan Hill, California, occurred at a relaxed site on the Calaveras fault. We examine how the more complex rate-and-state friction formalism based on laboratory experiments might have allowed the 1911 earthquake. Rate-and-state time-to-failure calculations are consistent with the occurrence of the 1911 event just 5 years after 1906 if the Calaveras fault was already close to failure before the effects of 1906. We also examine the likelihood that the entire 78 years of relative quiet (only four M ≥6 earthquakes) in the bay area after 1906 is consistent with rate-and-state assumptions, given that the previous 7 decades produced 18 M ≥6 earthquakes. Combinations of rate-and-state variables can be found that are consistent with this pattern of large bay area earthquakes, assuming that the rate of earthquakes in the 7 decades before 1906 would have continued had 1906 not occurred. These results demonstrate that rate-and-state offers a consistent explanation for the 78-year quiescence and the 1911 anomaly, although they do not rule out several alternate explanations.

California↗
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↗