Introduction to special section: Paleoseismology
No abstract available.
Journal of Geophysical Research B: Solid Earth↗
SEARCH · Geology Reports
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.
Coral ages and island subsidence, Hilo drill hole
A 25.8-m-thick sedimentary section containing coral fragments occurs directly below a surface lava flow (the ∼1340 year old Panaewa lava flow) at the Hilo drill hole. Ten coral samples from this section dated by accelerator mass spectrometry (AMS) radiocarbon and five by thermal infrared multispectral scanner (TIMS) 230 Th/U methods show good agreement. The calcareous unit is 9790 years old at the bottom and 1690 years old at the top and was deposited in a shallow lagoon behind an actively growing reef. This sedimentary unit is underlain by a 34-m-thick lava flow which in turn overlies a thin volcaniclastic silt with coral fragments that yield a single 14 C date of 10,340 years. The age-depth relations of the dated samples can be compared with proposed eustatic sea level curves after allowance for island subsidence is taken. Island subsidence averages 2.2 mm/yr for the last 47 years based on measurements from a tide gage near the drill hole or 2.5–2.6 mm/yr for the last 500,000 years based on the ages and depths of a series of drowned coral reefs offshore from west Hawaii. The age-depth measurements of coral fragments are more consistent with eustatic sea levels as determined by coral dating at Barbados and Albrolhos Islands than those based on oxygen isotopic data from deep sea cores. The Panaewa lava flow entered a lagoon underlain by coral debris and covered the drill site with 30.9 m of lava of which 11 m was above sea level. This surface has now subsided to 4.2 m above sea level, but it demonstrates how a modern lava flow entering Hilo Bay would not only change the coastline but could extensively modify the offshore shelf.
Journal of Geophysical Research B: Solid Earth↗
Stress/strain changes and triggered seismicity at The Geysers, California
The principal results of this study of remotely triggered seismicity in The Geysers geothermal field are the demonstration that triggering (initiation of earthquake failure) depends on a critical strain threshold and that the threshold level increases with decreasing frequency, or, equivalently, depends on strain rate. This threshold function derives from (1) analyses of dynamic strains associated with surface waves of the triggering earthquakes, (2) statistically measured aftershock zone dimensions, and (3) analytic functional representations of strains associated with power production and tides. The threshold is also consistent with triggering by static strain changes and implies that both static and dynamic strains may cause aftershocks. The observation that triggered seismicity probably occurs in addition to background activity also provides an important constraint on the triggering process. Assuming the physical processes underlying earthquake nucleation to be the same, Gomberg [this issue] discusses seismicity triggered by the M W 7.3 Landers earthquake, its constraints on the variability of triggering thresholds with site, and the implications of time delays between triggering and triggered earthquakes. Our results enable us to reject the hypothesis that dynamic strains simply nudge prestressed faults over a Coulomb failure threshold sooner than they would have otherwise. We interpret the rate-dependent triggering threshold as evidence of several competing processes with different time constants, the faster one(s) facilitating failure and the other(s) inhibiting it. Such competition is a common feature of theories of slip instability. All these results, not surprisingly, imply that to understand earthquake triggering one must consider not only simple failure criteria requiring exceedence of some constant threshold but also the requirements for generating instabilities.
Journal of Geophysical Research B: Solid Earth↗
Hydrogeology of the Hawaii Scientific Drilling Project borehole KP-1 1. Hydraulic conditions adjacent to the well bore
Temperature and formation resistivity logs obtained in borehole KP-1 of the Hawaii Scientific Drilling Project indicate that the adjacent formation is characterized by several zones of distinctly different average temperature and water salinity. A series of hydraulic analyses and water sampling programs were conducted to rule out the possibility of local hydraulic effects associated with the presence of the borehole in the generation of these apparent groundwater zones. Hydraulic tests and sampling with the borehole cased to a depth of 710 m and open below that depth indicate that the deep aquifer contains seawater at a temperature nearly identical to that of the open ocean at the same depth. Various analyses give estimates of aquifer transmissivity of about 10-3 m2/s in the vicinity of the borehole. Isolation of this deeper aquifer from the overlying groundwater zones was investigated by perforating the casing at six locations and then measuring the changes in water level in the borehole, in the salinity of the fluid column, in the temperature profile of the fluid column, and in the rate of flow in the fluid column induced by the perforations. These results positively confirm that the zones of distinctly different formation properties indicated on the temperature and resistivity logs are not caused by flow in or around casing. Flow and fluid column salinity induced by the perforations also confirm significant differences between the hydraulic heads and geochemistry of the different groundwater zones inferred from the well logs.
Journal of Geophysical Research B: Solid Earth↗
Thermal and mass implications of magmatic evolution in the Lassen volcanic region, California, and minimum constraints on basalt influx to the lower crust
We have analyzed the heat and mass demands of a petrologic model of basalt-driven magmatic evolution in which variously fractionated mafic magmas mix with silicic partial melts of the lower crust. We have formulated steady state heat budgets for two volcanically distinct areas in the Lassen region: the large, late Quaternary, intermediate to silicic Lassen volcanic center and the nearby, coeval, less evolved Caribou volcanic field. At Caribou volcanic field, heat provided by cooling and fractional crystallization of 52 km 3 of basalt is more than sufficient to produce 10 km of rhyolitic melt by partial melting of lower crust. Net heat added by basalt intrusion at Caribou volcanic field is equivalent to an increase in lower crustal heat flow of ∼7 mW m −2 , indicating that the field is not a major crustal thermal anomaly. Addition of cumulates from fractionation is offset by removal of erupted partial melts. A minimum basalt influx of 0.3 km 3 (km 2 Ma) −1 is needed to supply Caribou volcanic field. Our methodology does not fully account for an influx of basalt that remains in the crust as derivative intrusives. On the basis of comparison to deep heat flow, the input of basalt could be ∼3 to 7 times the amount we calculate. At Lassen volcanic center, at least 203 km 3 of mantle-derived basalt is needed to produce 141 km 3 of partial melt and drive the volcanic system. Partial melting mobilizes lower crustal material, augmenting the magmatic volume available for eruption at Lassen volcanic center; thus the erupted volume of 215 km 3 exceeds the calculated basalt input of 203 km 3 . The minimum basalt input of 1.6 km 3 (km 2 Ma) −1 is >5 times the minimum influx to the Caribou volcanic field. Basalt influx high enough to sustain considerable partial melting, coupled with locally high extension rate, is a crucial factor in development of Lassen volcanic center; in contrast, Caribou volcanic field has failed to develop into a large silicic center primarily because basalt supply there has been insufficient.
Journal of Geophysical Research B: Solid Earth↗
Rheology of water ices V and VI
We have measured the mechanical strength (σ) of pure water ices V and VI under steady state deformation conditions. Constant displacement rate compressional tests were conducted in a gas apparatus at confining pressures from 400 < P < 800 MPa, temperatures from 209 < T < 270 K, and strain rates 7 × 10 −7 < < 7 × 10 −4 s −1 . Most of the results fit to an empirical flow law of the form = A σ n exp {−( E * + PV *)/ RT }, where the four material constants A , n , E *, and V * are (for in inverse seconds and P and σ in megapascals) 10 23.0 , 6.0, 136 kJ/mol, and 29 cm 3 /mol, respectively, for ice V; and 10 6.7 , 4.5, 66 kJ/mol, and 11 cm 3 /mol, respectively, for ice VI. Ice VI may weaken to a mechanism of higher E * at T > 250 K. Ices V and VI are thus rheologically distinct but by coincidence have approximately the same strength under the conditions chosen for these experiments. To avoid misidentification, these tests are therefore accompanied by careful observations of the occurrences and characteristics of phase changes. One sample each of ice V and VI was quenched at pressure to metastably retain the high-pressure phase and the acquired deformation microstructures; X ray diffraction analysis of these samples confirmed the phase identification. Surface replicas of the deformed and quenched samples suggest that ice V probably deforms largely by dislocation creep, while ice VI deforms by a more complicated process involving substantial grain size reduction through recrystallization.
Journal of Geophysical Research B: Solid Earth↗
Cyclic stressing and seismicity at strongly coupled subduction zones
We use the finite element method to analyze stress variations in and near a strongly coupled subduction zone during an earthquake cycle. Deformation is assumed to be uniform along strike (plane strain on a cross section normal to the trench axis), and periodic earthquake slip is imposed consistent with the long-term rate of plate convergence and degree of coupling. Simulations of stress and displacement rate fields represent periodic fluctuations in time superimposed on an average field. The oceanic plate, descending slab, and continental lithosphere are assumed here to respond elastically to these fluctuations, and the remaining mantle under and between plates is assumed to respond as Maxwell viscoelastic. In the first part of the analysis we find that computed stress fluctuations in space and time are generally consistent with observed earthquake mechanism variations with time since a great thrust event. In particular, trench-normal extensional earthquakes tend to occur early in the earthquake cycle toward the outer rise but occur more abundantly late in the cycle in the subducting slab downdip of the main thrust zone. Compressional earthquakes, when they occur at all, have the opposite pattern. Our results suggest also that the actual timing of extensional outer rise events is controlled by the rheology of the shallow aseismic portion of the thrust interface. The second part of the analysis shows the effects of mantle relaxation on the rate of ground surface deformation during the earthquake cycle. Models without relaxation predict a strong overall compressional strain rate in the continental plate above the main thrust zone, with the strain rate constant between mainshocks. However with significant relaxation present, a localized region of unusually low compressional, or even slightly extensional, strain rate develops along the surface of the continental plate above and somewhat inland from the downdip edge of the locked main thrust zone. The low strain rate starts in the middle or late part of the cycle, depending on position. This result suggests that the negligible or small contraction measured on the Shumagin Islands, Alaska, during 1980 to 1991, may not invalidate an interpretation of that region as being a moderately coupled subduction zone. In contrast, mantle relaxation causes only modest temporal nonuniformity of uplift rates in the overriding plate and of extensional stress rates in the subducting plate, even when the Maxwell time is an order of magnitude less than the recurrence interval.
Journal of Geophysical Research B: Solid Earth↗
The 40Ar/39Ar and K/Ar dating of lavas from the Hilo 1-km core hole, Hawaii Scientific Drilling Project
Mauna Kea lava flows cored in the Hilo hole range in age from <200 ka to about 400 ka based on 40 Ar/ 39 Ar incremental heating and K‐Ar analyses of 16 groundmass samples and one coexisting plagioclase. The lavas, all subaerially deposited, include a lower section consisting only of tholeiitic basalts and an upper section of interbedded alkalic, transitional tholeiitic, and tholeiitic basalts. The lower section has yielded predominantly complex, discordant 40 Ar/ 39 Ar age spectra that result from mobility of 40 Ar and perhaps K, the presence of excess 40 Ar, and redistribution of 39 Ar by recoil. Comparison of K‐Ar ages with 40 Ar/ 39 Ar integrated ages indicates that some of these samples have also lost 39 Ar. Nevertheless, two plateau ages of 391 ± 40 and 400 ± 26 ka from deep in the hole, combined with data from the upper section, show that the tholeiitic section accumulated at an average rate of about 7 to 8 m/kyr and has a mean recurrence interval of 0.5 kyr/flow unit. Samples from the upper section yield relatively precise 40 Ar/ 39 Ar plateau and isotope correlation ages of 326 ± 23, 241 ± 5, 232 ± 4, and 199 ± 9 ka for depths of −415.7 m to −299.2 m. Within their uncertainty, these ages define a linear relationship with depth, with an average accumulation rate of 0.9 m/kyr and an average recurrence interval of 4.8 kyr/flow unit. The top of the Mauna Kea sequence at −280 m must be older than the plateau age of 132 ± 32 ka, obtained for the basal Mauna Loa flow in the corehole. The upward decrease in lava accumulation rate is a consequence of the decreasing magma supply available to Mauna Kea as it rode the Pacific plate away from its magma source, the Hawaiian mantle plume. The age‐depth relation in the core hole may be used to test and refine models that relate the growth of Mauna Kea to the thermal and compositional structure of the mantle plume.
Journal of Geophysical Research B: Solid Earth↗
Origin and depositional environment of clastic deposits in the Hilo drill hole, Hawaii
Volcaniclastic units cored at depths of about 87, 164, 178, 226, and 246 m below sea level and carbonate units located between depths of 27 and 53 m below sea level in the Hilo drill core were found to be deposited at or near sea level. Four of these units are hydroclastic deposits, formed when subaerially erupted Mauna Loa lava flows entered the ocean and fragmented to produce quenched, glassy fragments during hydrovolcanic explosions. Ash units 24 and 26, at 178 m depth, accumulated at sea level in a freshwater bog. They contain pyroxenes crystallized from tholeiitic magma that we infer erupted explosively at the summit of Kilauea volcano. Two carbon‐rich layers from these ashes have a weighted average radiocarbon age of 38.6 ± 0.9 ka; the ashes probably correlate with the oldest and thickest part of the Pahala ash. Ash unit 44, at the transition from Mauna Kea to Mauna Loa lava flows, was probably nearly 3.2 m thick and is inferred to be equivalent to the lower thick part of the composite Homelani ash mapped in Hilo and on the flanks of Mauna Kea. The age of this part of Homelani ash is between 128 ± 33 and 200 ± 10 ka; it may have erupted subglacially during the Pohakuloa glacial maxima on Mauna Kea. Beach sand units 12 and 22 were derived from nearby Mauna Loa and Mauna Kea lava flows. The middle of beach sand unit 38 was derived mainly from lava erupted near the distal end of the subaerial east rift zone of Kilauea volcano; these sands were transported about 33 km northwest to Hilo Bay by prevailing longshore currents. Combined age, depth, and sea level markers in the core allow us to determine that lava flow recurrence intervals averaged one flow every 4 kyr during the past 86 kyr and one flow every 16 kyr between 86 and 200 ka at the drill site and that major explosive eruptions that deposit thick ash in Hilo have occurred only twice in the last 400 kyr. These recurrence intervals support the moderate lava flow hazard zonation (zone 3) for coastal Hilo previously determined from surficial mapping.
Journal of Geophysical Research B: Solid Earth↗
Crustal structure of the Colorado Plateau, Arizona: Application of new long-offset seismic data analysis techniques
The Colorado Plateau is a large crustal block in the southwestern United States that has been raised intact nearly 2 km above sea level since Cretaceous marine sediments were deposited on its surface. Controversy exists concerning the thickness of the plateau crust and the source of its buoyancy. Interpretations of seismic data collected on the plateau vary as to whether the crust is closer to 40 or 50 km thick. A thick crust could support the observed topography of the Colorado Plateau isostatically, while a thinner crust would indicate the presence of an underlying low-density mantle. This paper reports results on long-offset seismic data collected during the 1989 segment of the U.S. Geological Survey Pacific to Arizona Crustal Experiment that extended from the Transition Zone into the Colorado Plateau in northwest Arizona. We apply two new methods to analyze long-offset data that employ finite difference travel time calculations: (1) a first-arrival time inverter to find upper crustal velocity structure and (2) a forward-modeling technique that allows the direct use of the inverted upper crustal solution in modeling secondary reflected arrivals. We find that the crustal thickness increases from 30 km beneath the metamorphic core complexes in the southern Basin and Range province to about 42 km beneath the northern Transition Zone and southern Colorado Plateau margin. We observe some crustal thinning (to ∼37 km thick) and slightly higher lower crustal velocities farther inboard; beneath the Kaibab uplift on the north rim of the Grand Canyon the crust thickens to a maximum of 48 km. We observe a nonuniform crustal thickness beneath the Colorado Plateau that varies by ∼15% and corresponds approximately to variations in topography with the thickest crust underlying the highest elevations. Crustal compositions (as inferred from seismic velocities) appear to be the same beneath the Colorado Plateau as those in the Basin and Range province to the southwest, implying that the plateau crust represents an unextended version of the Basin and Range. Some of the variability in crustal structure appears to correspond to preserved lithospheric discontinuities that date back to the Proterozoic Era.
Journal of Geophysical Research B: Solid Earth↗
Analysis of fractures intersecting Kahi Puka Well 1 and its relation to the growth of the island of Hawaii
As part of the Hawaii Scientific Drilling Project, Kahi Puka Well 1 penetrated about 275 m of Mauna Loa basalts overlying a sequence of Mauna Kea flow units as it was drilled and cored to a total depth of 1053 m below land surface. A borehole televiewer (BHTV) was run in most of the well in successive stages prior to casing in order to obtain magnetically oriented acoustic images of the borehole wall. A total of 283 individual fractures were identified from this log and characterized in terms of strike and dip. These data are divided into three vertical sections based upon age and volcanic source, and lower hemisphere stereographic plots identify two predominant, subparallel fracture subsets common to each section. Assuming that most of the steeply dipping fractures observed in the BHTV log are tensile features generated within basalt flows during deposition and cooling, this fracture information can be combined with models of the evolution of the island of Hawaii to investigate the depositional history of these Mauna Loa and Mauna Kea basalts over the past 400 kyr. The directions of high-angle fractures appear to be generally parallel to topography or to the coastline at the time of deposition, as is supported by surface mapping of modern flows. Consequently, an overall counterclockwise rotation of about 75?? in the strike of these fractures from the bottom to the top of the well represents a systematic change in depositional slope direction over time. We attribute the observed rotation in the orientations of the two predominant fracture subsets over the past 400 kyr to changes in the configurations of volcanic sources during shield building and to the structural interference of adjacent volcanoes that produces shifts in topographic patterns.
Journal of Geophysical Research B: Solid Earth↗
Along-axis variability in crustal accretion at the Mid-Atlantic Ridge: Results from the OCEAN study
The OCEAN experiment is an integrated geophysical study of a region of the Cape Verde abyssal plain that formed at 140 Ma. Deep seismic reflection and ocean bottom hydrophone (OBH) refraction data were acquired along lines parallel and perpendicular to the paleoridge axis trend identified from a detailed magnetic anomaly survey. The igneous basement is overlain by about 1.3 km of sediment which enables improved imaging of intracrustal structure beyond that possible near the Mid-Atlantic Ridge axis. We describe the results of a 150-km long profile oriented parallel to magnetic anomalies M15 and M16, along which deep seismic reflection data collected by the British Institutions Reflection Profiling Syndicate are complemented by refraction data constrained by four OBHs. The line spans an entire spreading segment between two fracture zones; the northern of which has an offset of 40 km and the other (central) has an offset of only 10 km. Away from the fracture zones, the mean igneous crustal thickness is 7.2 km; near both fracture zones, thinning of up to 4 km is observed, giving a mean igneous crustal thickness over the whole segment of approximately 6.5 km. Differences are seen between the two fracture zones in their seismic velocity structure, in the associated basement topography, and in the presence of a strong reflection extending into the mantle beneath the northern fracture zone. The boundary between oceanic layers 2 and 3 correlates with variably coherent normal incidence reflections and a change in the character of the reflectivity. A number of planar reflections up to 10 km in length are present within the middle and lower crust, dipping outward from beneath low-amplitude basement highs at ∼15°; these appear to be present only within layer 3. The Moho has several expressions in the reflection data, including isolated reflection events, a local increase in reflected amplitudes, and a downward decrease in coherent reflections. At the center of the segment there is a zone at the base of the crust within which both high- and low-velocity materials are present. This zone shows an enhanced level of discontinuous normal incidence reflectivity and may represent an initial fractionation event as melt was emplaced at the spreading ridge.
Journal of Geophysical Research B: Solid Earth↗
The upper mantle structure of the central Rio Grande rift region from teleseismic P and S wave travel time delays and attenuation
The lithosphere beneath a continental rift should be significantly modified due to extension. To image the lithosphere beneath the Rio Grande rift (RGR), we analyzed teleseismic travel time delays of both P and S wave arrivals and solved for the attenuation of P and S waves for four seismic experiments spanning the Rio Grande rift. Two tomographic inversions of the P wave travel time data are given: an Aki-Christofferson-Husebye (ACH) block model inversion and a downward projection inversion. The tomographic inversions reveal a NE-SW to NNE-SSW trending feature at depths of 35 to 145 km with a velocity reduction of 7 to 8% relative to mantle velocities beneath the Great Plains. This region correlates with the transition zone between the Colorado Plateau and the Rio Grande rift and is bounded on the NW by the Jemez lineament, a N52°E trending zone of late Miocene to Holocene volcanism. S wave delays plotted against P wave delays are fit with a straight line giving a slope of 3.0 ± 0.4. This correlation and the absolute velocity reduction imply that temperatures in the lithosphere are close to the solidus, consistent with, but not requiring, the presence of partial melt in the mantle beneath the Rio Grande rift. The attenuation data could imply the presence of partial melt. We compare our results with other geophysical and geologic data. We propose that any north-south trending thermal (velocity) anomaly that may have existed in the upper mantle during earlier (Oligocene to late Miocene) phases of rifting and that may have correlated with the axis of the rift has diminished with time and has been overprinted with more recent structure. The anomalously low-velocity body presently underlying the transition zone between the core of the Colorado Plateau and the rift may reflect processes resulting from the modern (Pliocene to present) regional stress field (oriented WNW-ESE), possibly heralding future extension across the Jemez lineament and transition zone.
Journal of Geophysical Research B: Solid Earth↗
Volatile emissions from the crater and flank of Oldoinyo Lengai volcano, Tanzania
As a comparison to airborne infrared (IR) flux measurements, ground-based sampling of fumarole and soil gases was used to characterize the quiescent degassing of CO 2 from Oldoinyo Lengai volcano. Aerial and ground-based measurements are in good agreement: ∼75% of the aerially measured CO 2 flux at Lengai (0.05–0.06 × 10 12 mol yr −1 or 6000–7200 tonnes CO 2 d −1 ) can be attributed to seven large crater vents. In contrast to Etna and Vulcano Island, where 15–50% of the total CO 2 flux emanates diffusely through the volcanic flanks, diffuse emissions were measured only within 500 m of the crater rim at Lengai, contributing <2% of the total flux. The lack of extensive flank emissions may reflect the dimensions of the magma chamber and/or the lack of a shallow fluid flow system. Thermodynamic restoration of fumarole analyses shows that gases are the most CO 2 -rich and H 2 O-poor reported for any volcano, containing 64–74% CO 2 , 24–34% H 2 O, 0.88–1.0% H 2 , 0.1–0.4% CO and <0.1% H 2 S, HCl, HF, and CH 4 . Volatile emissions of S, Cl, and F at Oldoiyno Lengai are estimated as 4.5, 1.5, and 1.0 × 10 7 mol yr −1 , respectively. Accuracy of the airborne technique was also assessed by measuring the C emission rate from a coal-burning power plant. CO 2 fluxes were measured within ±10% near the plant; however, poor resolution at increased distances caused an underestimation of the flux by a factor of 2. The relatively large CO 2 fluxes measured for alkaline volcanoes such as Oldoinyo Lengai or Etna may indicate that midplate volcanoes represent a large, yet relatively unknown, natural source of CO 2 .
Journal of Geophysical Research B: Solid Earth↗
Isotopic evolution of Mauna Kea volcano: Results from the initial phase of the Hawaii Scientific Drilling Project
We have examined the Sr, Nd, and Pb isotopic compositions of Mauna Kea lavas recovered by the first drilling phase of the Hawaii Scientific Drilling Project. These lavas, which range in age from ∼200 to 400 ka, provide a detailed record of chemical and isotopic changes in basalt composition during the shield/postshield transition and extend our record of Mauna Kea volcanism to a late-shield period roughly equivalent to the last ∼100 ka of Mauna Loa activity. Stratigraphic variations in isotopic composition reveal a gradual shift over time toward a more depleted source composition (e.g., higher 143 Nd/ 144 Nd, lower 87 Sr/ 86 Sr, and lower 3 He/ 4 He). This gradual evolution is in sharp contrast with the abrupt appearance of alkalic lavas at ∼240 ka recorded by the upper 50 m of Mauna Kea lavas from the core. Intercalated tholeiitic and alkalic lavas from the uppermost Mauna Kea section are isotopically indistinguishable. Combined with major element evidence (e.g., decreasing SiO 2 and increasing FeO) that the depth of melt segregation increased during the transition from tholeiitic to alkalic volcanism, the isotopic similarity of tholeiitic and alkalic lavas argues against significant lithosphere involvement during melt generation. Instead, the depleted isotopic signatures found in late shield-stage lavas are best explained by increasing the proportion of melt generated from a depleted upper mantle component entrained and heated by the rising central plume. Direct comparison of Mauna Kea and Mauna Loa lavas erupted at equivalent stages in these volcanoes' life cycles reveals persistent chemical and isotopic differences independent of the temporal evolution of each volcano. The oldest lavas recovered from the drillcore are similar to modern Kilauea lavas, but are distinct from Mauna Loa lavas. Mauna Kea lavas have higher 143 Nd/ 144 Nd and 206 Pb/ 204 Pb and lower 87 Sr/ 86 Sr. Higher concentrations of incompatible trace elements in primary magmas, lower SiO 2 , and higher FeO also indicate that Mauna Kea lavas formed through smaller degrees of partial melting at greater depth than Mauna Loa lavas. These chemical and isotopic differences are consistently found between volcanoes along the western “Loa” and eastern “Kea” trends and reflect large-scale variations in source composition and melting environment. We propose a simple model of a radially zoned plume centered beneath the Loa trend. Loa trend lavas generated from the hot plume axis reflect high degrees of partial melting from a source containing a mixture of enriched plume-source material and entrained lower mantle. Kea trend lavas, in contrast, are generated from the cooler, peripheral portions of the plume, record lower degrees of partial melting, and tap a source containing a greater proportion of depleted upper mantle.
Journal of Geophysical Research B: Solid Earth↗
Frictional constraints on crustal faulting
We consider how variations in fault frictional properties affect the phenomenology of earthquake faulting. In particular, we propose that lateral variations in fault friction produce the marked heterogeneity of slip observed in large earthquakes. We model these variations using a rate- and state-dependent friction law, where we differentiate velocity-weakening behavior into two fields: the strong seismic field is very velocity weakening and the weak seismic field is slightly velocity weakening. Similarly, we differentiate velocity-strengthening behavior into two fields: the compliant field is slightly velocity strengthening and the viscous field is very velocity strengthening. The strong seismic field comprises the seismic slip concentrations, or asperities. The two “intermediate” fields, weak seismic and compliant, have frictional velocity dependences that are close to velocity neutral: these fields modulate both the tectonic loading and the dynamic rupture process. During the interseismic period, the weak seismic and compliant regions slip aseismically, while the strong seismic regions remain locked, evolving into stress concentrations that fail only in main shocks. The weak seismic areas exhibit most of the interseismic activity and aftershocks but can also creep seismically. This “mixed” frictional behavior can be obtained from a sufficiently heterogeneous distribution of the critical slip distance. The model also provides a mechanism for rupture arrest: dynamic rupture fronts decelerate as they penetrate into unloaded complaint or weak seismic areas, producing broad areas of accelerated afterslip. Aftershocks occur on both the weak seismic and compliant areas around a fault, but most of the stress is diffused through aseismic slip. Rapid afterslip on these peripheral areas can also produce aftershocks within the main shock rupture area by reloading weak fault areas that slipped in the main shock and then healed. We test this frictional model by comparing the seismicity and the coseismic slip for the 1966 Parkfield, 1979 Coyote Lake, and 1984 Morgan Hill earthquakes. The interevent seismicity and aftershocks appear to occur on fault areas outside the regions of significant slip: these regions are interpreted as either weak seismic or compliant, depending on whether or not they manifest interevent seismicity.
Journal of Geophysical Research B: Solid Earth↗
Frictional behavior of large displacement experimental faults
The coefficient of friction and velocity dependence of friction of initially bare surfaces and 1-mm-thick simulated fault gouges (<90 μm) of Westerly granite were determined as a function of displacement to >400 mm at 25°C and 25 MPa normal stress. Steady state negative friction velocity dependence and a steady state fault zone microstructure are achieved after ∼18 mm displacement, and an approximately constant strength is reached after a few tens of millimeters of sliding on initially bare surfaces. Simulated fault gouges show a large but systematic variation of friction, velocity dependence of friction, dilatancy, and degree of localization with displacement. At short displacement (<10 mm), simulated gouge is strong, velocity strengthening and changes in sliding velocity are accompanied by relatively large changes in dilatancy rate. With continued displacement, simulated gouges become progressively weaker and less velocity strengthening, the velocity dependence of dilatancy rate decreases, and deformation becomes localized into a narrow basal shear which at its most localized is observed to be velocity weakening. With subsequent displacement, the fault restrengthens, returns to velocity strengthening, or to velocity neutral, the velocity dependence of dilatancy rate becomes larger, and deformation becomes distributed. Correlation of friction, velocity dependence of friction and of dilatancy rate, and degree of localization at all displacements in simulated gouge suggest that all quantities are interrelated. The observations do not distinguish the independent variables but suggest that the degree of localization is controlled by the fault strength, not by the friction velocity dependence. The friction velocity dependence and velocity dependence of dilatancy rate can be used as qualitative measures of the degree of localization in simulated gouge, in agreement with previous studies. Theory equating the friction velocity dependence of simulated gouge to the sum of the friction velocity dependence of bare surfaces and the velocity dependence of dilatancy rate of simulated gouge fails to quantitatively account for the experimental observations.
Journal of Geophysical Research B: Solid Earth↗
Type of faulting and orientation of stress and strain as a function of space and time in Kilauea's south flank, Hawaii
Earthquake focal mechanisms of events occurring between 1972 and 1992 in the south flank of Kilauea volcano, Hawaii, are used to infer the state of stress and strain as a function of time and space. We have determined 870 fault plane solutions from P wave first motion polarities for events with magnitudes M L ≥ 2.5 and depth ranging between 6 and 12 km. Faulting is characterized by a mixture of decollement, reverse, and normal faults. Most large earthquakes with magnitude M < 7 slip on reverse faults striking NE at 40° and dipping SE between 60° and 70°. In Hawaii, the earthquakes with M > 7 rupture the decollement plane, since it is the only surface large enough to generate magnitude 7 or larger earthquakes. The percentage of reverse faulting events is high compared to the decollement and normal faulting mechanisms for the period 1972–1983. The percentage of decollement type focal mechanisms becomes dominant after 1983. This pattern of faulting activity suggests that pressure was building up within Kilauea's rift zone prior to the 1983 Puu'Oo eruption. Overall, a single stress orientation with the maximum compressive stress oriented SE perpendicular to the rift and dipping at 45° is compatible with the coeval existence of decollement, reverse, and normal faults. However, in a crustal volume east of longitude 155°10′W, we find a change of the orientation of σ 1 from nearly horizontal to plunging 45° SE occurring in 1979. This stress rotation suggests magma movements within the aseismic part of Kilauea's east rift zone. The strain and stress orientations are coaxial in the south flank except within the volume where the stress rotation is observed. We observe a change in the relationship between stress and strain directions caused either by the shifting of seismic activity from reverse faults to decollements, while stress stays constant, or by a rotation of stress, while strain remains constant. Assuming that the model of a noncohesive Coulomb wedge is appropriate for Kilauea's south flank, we find that high pore pressures are prevalent along the decollement and within the wedge for a coefficient of friction equal to 0.85.
Hawaii↗