Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research - Planets”

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,099 records · Page 61Linked to original sources
USGS1996

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

Satellite imagery of the onset of streaming flow of ice streams C and D, West Antarctica

Five overlapping Landsat multispectral scanner satellite images of the interior of the West Antarctic ice sheet were enhanced with principal component analysis, high-pass filtering, and linear contrast stretching and merged into a mosaic by aligning surface features in the overlap areas. The mosaic was registered to geodetic coordinates, to an accuracy of about 1 km, using the five scene centers as control points. The onset of streaming flow of two tributaries of ice stream C and one tributary of ice stream D is visible in the mosaic. The onset appears to occur within a relatively short distance, less than the width of the ice stream, typically at a subglacial topographic feature, such as a step or ridge. The ice streams extend farther up into the interior than previously mapped. Ice stream D starts about 150 km from the ice divide, at an altitude of about 1500 m, approximately halfway up the convex-upward dome shape of the interior ice sheet. Ice stream D is relatively much longer than ice stream C, possibly because ice stream D is currently active, whereas ice stream C is currently inactive. The grounded portion of the West Antarctic ice sheet is perhaps best conceptualized as an ice sheet in which ice streams are embedded over most of its area, with slow-moving ice converging into fast-moving ice streams in a widely distributed pattern, much like that of streams and rivers in a hydrologic basin. A relic margin appears to parallel most of the south margin of the tributary of ice stream D, separated from the active shear margin by about 10 km or less for a distance of over 200 km. This means there is now evidence for recent changes having occurred in three of the five major ice streams which drain most of West Antarctica (B, C, and D), two of which (B and D) are currently active.

Journal of Geophysical Research C: Oceans↗
USGS1996

Southern Ocean monthly wave fields for austral winters 1985-1988 by Geosat radar altimeter

Four years of monthly averaged wave height fields for the austral winters 1985–1988 derived from the Geosat altimeter data show a spatial variability of the scale of 500–1000 km that varies monthly and annually. This variability is superimposed on the zonal patterns surrounding the Antarctic continent and characteristic of the climatology derived from the U.S. Navy [1992] Marine Climatic Atlas of the World. The location and the intensity of these large-scale features, which are not found in the climatological fields, exhibit strong monthly and yearly variations. A global underestimation of the climatological mean wave heights by more than 1 m is also found over large regions of the Southern Ocean. The largest monthly averaged significant wave heights are above 5 m and are found during August of every year in the Indian Ocean, south of 40°S. The monthly wave fields show more variability in the Atlantic and Pacific Oceans than in the Indian Ocean. The Seasat data from 1978 and the Geosat data from 1985 and 1988 show an eastward rotation of the largest wave heights. However, this rotation is absent in 1986 and 1987; the former was a year of unusually low sea states, and the latter was a year of unusually high sea states, which suggests a link to the El Niño-Southern Oscillation event of 1986.

Journal of Geophysical Research C: Oceans↗
USGS1996

Cadmium in the California Current system: Tracer of past and present upwelling

Over 100 samples were collected off the west coast of North America during 1991–1993 to determine the relation between wind‐driven upwelling and nearshore concentrations of dissolved silicate (Si), phosphate (P), and cadmium (Cd). Highly enriched in deep water offshore, these constituents are sensitive indicators of upwelling. Coastal water was sampled from the shore in January and June 1992 at 12 sites distributed between 36° and 48°N latitude. In January the composition of nearshore water along this transect was fairly uniform: 5–15 μmol/kg for Si, 0.5 to 1.0 μmol/kg for P, and 0.1–0.3 nmol/kg for Cd. In June, elevated concentrations of Si (30 μmol/kg), P (2.0 μmol/kg), and Cd (0.6 nmol/kg) revealed a region of intense upwelling between 38° and 40°N. The pattern is broadly consistent with meridional gradients in coastal upwelling calculated from the long‐term mean of alongshore winds compiled from ship reports. Nearshore water was also collected biweekly to monthly at two sites 3 km apart near San Francisco Bay (37.5°N) during 1991–1993. The variability seen in the time series suggests that the composition of nearshore water integrates the effect of alongshore winds over timescales of several weeks. Seasonal variations in Si (5–50 μmol/kg), P (0.5–2.5 μmol/kg), and Cd (0.1–0.8 nmol/kg) concentrations were consistent with upwelling during spring and summer. Maximum Si, P, and Cd concentrations reached in May 1991 were consistent with advection to the very nearshore region from a depth of about 300 m relative to a vertical profile at a distance of 200 km from the coast. Nearshore Si, P, and Cd concentrations were reduced relative to 1991 in 1992, and, to a lesser extent, in 1993 due to weaker upwelling linked to the warm phase of the El Niño‐Southern Oscillation. During periods of weaker upwelling or downwelling, variations in P, Si, and Cd concentrations became uncoupled. There is a good correlation between the coastal Cd time series near San Francisco Bay (37.5°N) and a second order polynomial function of the upwelling index of Bakun [1975] at 36°N, filtered with a 30‐day running mean ( r 2 = 0.71, n = 39). The index is a daily estimate of coastal upwelling calculated from 6‐hourly mean atmospheric pressure distributions at 36°N. From this function and a record of daily upwelling indices, we infer a range of annually averaged coastal Cd concentrations of at least 0.3–0.5 nmol/kg since 1967. Cd/Ca ratios in shells of foraminifera from San Francisco Bay suggest that average coastal Cd concentrations 3500–4500 years ago were at the upper end of this range.

Journal of Geophysical Research C: Oceans↗
USGS1996

Potential climatic impacts of vegetation change: A regional modeling study

The human species has been modifying the landscape long before the development of modern agrarian techniques. Much of the land area of the conterminous United States is currently used for agricultural production. In certain regions this change in vegetative cover from its natural state may have led to local climatic change. A regional climate version of the Colorado State University Regional Atmospheric Modeling System was used to assess the impact of a natural versus current vegetation distribution on the weather and climate of July 1989. The results indicate that coherent regions of substantial changes, of both positive and negative sign, in screen height temperature, humidity, wind speed, and precipitation are a possible consequence of land use change throughout the United States. The simulated changes in the screen height quantities were closely related to changes in the vegetation parameters of albedo, roughness length, leaf area index, and fractional coverage.

Journal of Geophysical Research D: Atmospheres↗
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↗