USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research B: 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 451 records · Page 25Linked to original sources
USGS1996

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

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

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

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

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

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

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

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

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

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

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

Shear wave velocity structure in North America from large-scale waveform inversions of surface waves

A two-step nonlinear and linear inversion is carried out to map the lateral heterogeneity beneath North America using surface wave data. The lateral resolution for most areas of the model is of the order of several hundred kilometers. The most obvious feature in the tomographic images is the rapid transition between low velocities in the tectonically active region west of the Rocky Mountains and high velocities in the stable central and eastern shield of North America. The model also reveals smaller-scale heterogeneous velocity structures. A high-velocity anomaly is imaged beneath the state of Washington that could be explained as the subducting Juan de Fuca plate beneath the Cascades. A large low-velocity structure extends along the coast from the Mendocino to the Rivera triple junction and to the continental interior across the southwestern United States and northwestern Mexico. Its shape changes notably with depth. This anomaly largely coincides with the part of the margin where no lithosphere is consumed since the subduction has been replaced by a transform fault. Evidence for a discontinuous subduction of the Cocos plate along the Middle American Trench is found. In central Mexico a transition is visible from low velocities across the Trans-Mexican Volcanic Belt (TMVB) to high velocities beneath the Yucatan Peninsula. Two elongated low-velocity anomalies beneath the Yellowstone Plateau and the eastern Snake River Plain volcanic system and beneath central Mexico and the TMVB seem to be associated with magmatism and partial melting. Another low-velocity feature is seen at depths of approximately 200 km beneath Florida and the Atlantic Coastal Plain. The inversion technique used is based on a linear surface wave scattering theory, which gives tomographic images of the relative phase velocity perturbations in four period bands ranging from 40 to 150 s. In order to find a smooth reference model a nonlinear inversion based on ray theory is first performed. After correcting for the crustal thickness the phase velocity perturbations obtained from the subsequent linear waveform inversion for the different period bands are converted to a three-layer model of S velocity perturbations (layer 1, 25–100 km; layer 2, 100–200 km; layer 3, 200–300 km). We have applied this method on 275 high-quality Rayleigh waves recorded by a variety of instruments in North America (IRIS/USGS, IRIS/IDA, TERRAscope, RSTN). Sensitivity tests indicate that the lateral resolution is especially good in the densely sampled western continental United States, Mexico, and the Gulf of Mexico.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Uplift of the Colorado Plateau due to lithosphere attenuation during Laramide low-angle subduction

The Colorado Plateau is blanketed by Phanerozoic marine and nonmarine strata as young as Cretaceous that are now exposed at elevations of about 2 km. Crustal thickening due to magmatism and horizontal crustal shortening was far less than necessary to cause this uplift, which is commonly attributed to the consequences of mantle lithosphere thinning and heating. The Colorado Plateau and the midcontinent region around Iowa consist of Precambrian bedrock overlain by a similar amount of Paleozoic platformal strata, and thus both regions once had similar lithospheric buoyancy. Mesozoic sedimentation increased the crustal thickness and lithospheric buoyancy of the Colorado Plateau relative to the midcontinent region. Backstripping calculations yield elevation without these sediments and lead to a calculated elevation difference between the two areas of about 1200 m, which represents unexplained plateau uplift. Review of constraints on uplift timing finds little support for a late Cenozoic uplift age and allows early to middle Cenozoic uplift, which is consistent with uplift mechanisms related to low-angle subduction that ended in the middle Cenozoic. Finite element heat flow calculations of low-angle subduction and lithosphere attenuation, using a range of initial lithosphere thicknesses and degree of attenuation, indicate that required uplift can result from tectonic removal of about 120 km of mantle lithosphere from an initially 200-km-thick lithosphere. This allows for partial preservation of North American mantle lithosphere with its distinctive isotopic signature in some late Cenozoic volcanic rocks and is consistent with normal Pn velocities in the uppermost mantle beneath the plateau.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Three-dimensional models of deformation near strike-slip faults

We use three-dimensional elastic models to help guide the kinematic interpretation of crustal deformation associated with strike-slip faults. Deformation of the brittle upper crust in the vicinity of strike-slip fault systems is modeled with the assumption that upper crustal deformation is driven by the relative plate motion in the upper mantle. The driving motion is represented by displacement that is specified on the bottom of a 15-km-thick elastic upper crust everywhere except in a zone of finite width in the vicinity of the faults, which we term the “shear zone.” Stress-free basal boundary conditions are specified within the shear zone. The basal driving displacement is either pure strike slip or strike slip with a small oblique component, and the geometry of the fault system includes a single fault, several parallel faults, and overlapping en echelon faults. We examine the variations in deformation due to changes in the width of the shear zone and due to changes in the shear strength of the faults. In models with weak faults the width of the shear zone has a considerable effect on the surficial extent and amplitude of the vertical and horizontal deformation and on the amount of rotation around horizontal and vertical axes. Strong fault models have more localized deformation at the tip of the faults, and the deformation is partly distributed outside the fault zone. The dimensions of large basins along strike-slip faults, such as the Rukwa and Dead Sea basins, and the absence of uplift around pull-apart basins fit models with weak faults better than models with strong faults. Our models also suggest that the length-to-width ratio of pull-apart basins depends on the width of the shear zone and the shear strength of the faults and is not constant as previously suggested. We show that pure strike-slip motion can produce tectonic features, such as elongate half grabens along a single fault, rotated blocks at the ends of parallel faults, or extension perpendicular to overlapping en echelon faults, which can be misinterpreted to indicate a regional component of extension. Zones of subsidence or uplift can become wider than expected for transform plate boundaries when a minor component of oblique motion is added to a system of parallel strike-slip faults.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

The 1954 Rainbow Mountain-Fairview Peak-Dixie Valley earthquakes: A triggered normal faulting sequence

In 1954, four earthquakes of M > 6.0 occurred within a 30 km radius in a period of six months. The Rainbow Mountain-Fairview Peak-Dixie Valley earthquakes are among the largest to have been recorded geodetically in the Basin and Range province. The Fairview Peak earthquake (M = 7.2, December 12, 1954) followed two events in the Rainbow Mountains ( M = 6.2, July 6, and M = 6.5, August 24, 1954) by 6 months. Four minutes later the Dixie Valley fault ruptured ( M = 6.7, December 12, 1954). The changes in static stresses caused by the events are calculated using the Coulomb-Navier failure criterion and assuming uniform slip on rectangular dislocations embedded in an elastic half-space. Coulomb stress changes are resolved on optimally oriented faults and on each of the faults that ruptured in the chain of events. These calculations show that each earthquake in the Rainbow Mountain-Fairview Peak-Dixie Valley sequence was preceded by a static stress change that encouraged failure. The magnitude of the stress increases transferred from one earthquake to another ranged from 0.01 MPa (0.1 bar) to over 0.1 MPa (1 bar). Stresses were reduced by up to 0.1 MPa over most of the Rainbow Mountain-Fairview Peak area as a result of the earthquake sequence.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Geometry of the 1954 Fairview Peak-Dixie Valley earthquake sequence from a joint inversion of leveling and triangulation data

In 1954, four earthquakes greater than M s = 6.0 occurred within a 30-km radius and in a period of 6 months. Elevation and angle changes calculated from repeated leveling and triangulation surveys which span the coseismic period provide constraints on the fault geometries and coseismic slip of the faults which were activated. The quality of the coseismic geodetic data is assessed. Corrections are applied to the leveling data for subsidence due to groundwater withdrawal in the Fallon area, and a rod miscalibration error of 150 ± 30 ppm is isolated in leveling surveys made in 1967. The leveling and triangulation observations are then simultaneously inverted using the single value decomposition (SVD) inversion method to determine fault geometries and coseismic slip. Using SVD, it is possible to determine on which faults slip is resolvable given the data distribution. The faults are found to dip between 50° and 80° and extend to depths of 5 to 14 km. The geodetically derived slip values are generally equal to, or greater than, the maximum observed displacement along the surface scarps. Where slip is resolvable the geodetic data indicates the 1954 sequence contained a significant component of right-lateral slip. This is consistent with the N15°W trending shear zone which geodetic surveys have detected in western Nevada.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

The 1994 Northridge, California, earthquake: Investigation of rupture velocity, risetime, and high-frequency radiation

A hybrid global search algorithm is used to solve the nonlinear problem of calculating slip amplitude, rake, risetime, and rupture time on a finite fault. Thirty-five strong motion velocity records are inverted by this method over the frequency band from 0.1 to 1.0 Hz for the Northridge earthquake. Four regions of larger-amplitude slip are identified: one near the hypocenter at a depth of 17 km, a second west of the hypocenter at about the same depth, a third updip from the hypocenter at a depth of 10 km, and a fourth updip from the hypocenter and to the northwest. The results further show an initial fast rupture with a velocity of 2.8 to 3.0 km/s followed by a slow termination of the rupture with velocities of 2.0 to 2.5 km/s. The initial energetic rupture phase lasts for 3 s, extending out 10 km from the hypocenter. Slip near the hypocenter has a short risetime of 0.5 s, which increases to 1.5 s for the major slip areas removed from the hypocentral region. The energetic rupture phase is also shown to be the primary source of high-frequency radiation (1–15 Hz) by an inversion of acceleration envelopes. The same global search algorithm is used in the envelope inversion to calculate high-frequency radiation intensity on the fault and rupture time. The rupture timing from the low- and high-frequency inversions is similar, indicating that the high frequencies are produced primarily at the mainshock rupture front. Two major sources of high-frequency radiation are identified within the energetic rupture phase, one at the hypocenter and another deep source to the west of the hypocenter. The source at the hypocenter is associated with the initiation of rupture and the breaking of a high-stress-drop asperity and the second is associated with stopping of the rupture in a westerly direction.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Seismic velocities for hydrate-bearing sediments using weighted equation

A weighted equation based on the three-phase time-average and Wood equations is applied to derive a relationship between the compressional wave (P wave) velocity and the amount of hydrates filling the pore space. The proposed theory predicts accurate P wave velocities of marine sediments in the porosity range of 40-80% and provides a practical means of estimating the amount of in situ hydrate using seismic velocity. The shear (S) wave velocity is derived under the assumption that the P to S wave velocity ratio of the hydrated sediments is proportional to the weighted average of the P to S wave velocity ratios of the constituent components of the sediment. In the case that all constituent components are known, a weighted equation using multiphase time-average and Wood equations is possible. However, this study showed that a three-phase equation with modified matrix velocity, compensated for the clay content, is sufficient to accurately predict the compressional wave velocities for the marine sediments. This theory was applied to the laboratory measurements of the P and S wave velocities in permafrost samples to infer the amount of ice in the unconsolidated sediment. The results are comparable to the results obtained by repeatedly applying the two-phase wave scattering theory. The theory predicts that the Poisson's ratio of the hydrated sediments decreases as the hydrate concentration increases and the porosity decreases. In consequence, the amplitude versus offset (AVO) data for the bottom-simulating reflections may reveal positive, negative, or no AVO anomalies depending on the concentration of hydrates in the sediments.

Journal of Geophysical Research B: Solid Earth↗