USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geophysical Journal”

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,621 records · Page 90Linked to original sources
USGS2012

Interactive coupling of a lake thermal model with a regional climate model

A one‐dimensional model of lake temperature, evaporation, and ice has been coupled in an interactive mode with the climate version of the National Center for Atmospheric Research/Pennsylvania State University regional (mesoscale) atmospheric model (MM4). The coupled MM4‐lake model makes possible high‐resolution simulations of climate in the proximity of large water bodies. Atmospheric variables required as input for the lake model are supplied by MM4 and simulated values of lake temperature, evaporation, and ice cover axe supplied to MM4 by the lake model. We have tested the coupled model system with a 60‐day, summertime simulation at Pyramid Lake, Nevada, and with a 10‐day, wintertime simulation of the North American Great Lakes and vicinity. Both simulations were conducted at a 60‐km resolution. Results from these tests indicate the coupled model system produces realistic simulations of lake temperature, evaporation, and ice cover and that the coupled system is applicable to simulations of regional climate change.

Journal of Geophysical Research-Atmospheres↗
USGS2012

The geology and distribution of impact craters on Venus: What are they telling us?

Magellan has revealed an ensemble of impact craters on Venus that is unique in many important ways. We have compiled a data base describing the 842 craters on 89% of Venus' surface mapped through orbit 2578. (The craters range in diameter from 1.5 to 280 km.) We have studied the distribution, size‐density, morphology, geology, and associated surface properties of these craters both in the aggregate and, for some craters, in greater detail. We find that (1) the spatial distribution of craters is highly uniform; (2) the size‐density distribution of large craters (diameters ≥35 km) is similar to the young crater populations on other terrestrial planets but at a much lower density that indicates an average age of only about 0.5 Ga (based on the estimated population of Venus‐crossing asteroids); (3) unlike the case on other planets, the density of small craters (diameters ≤35 km) declines rapidly with decreasing diameters because of atmospheric filtering; (4) the spectrum of crater modification differs greatly from that on other planets: 62% of all craters are pristine, only 4% are embayed by lavas, and the remainder are affected by tectonism, but none are severely and progressively depleted (as extrapolated from the size‐density distribution of larger craters); (5) large craters have a progression of morphologies generally similar to those on other planets, but small craters are typically irregular or multiple rather than bowl shaped; (6) diffuse radar‐bright or ‐dark features surround some craters, and 367 similar diffuse “splotches” with no central crater are observed; and (7) other crater features unique to Venus include radar‐bright or ‐dark parabolic arcs opening westward and extensive outflows originating in crater ejecta. The first three of these observations are entirely unexpected. We interpret them as indicating that the planet's cratering record was erased by a global resurfacing event or events, the latest ending about 0.5 Ga, after which volcanic activity declined (but did not cease entirely). Since the last resurfacing event, a maximum of 10% of the planet has been resurfaced and only about 4% of the craters have been obliterated. Convective thermal evolution models support this interpretation (Arkani‐Hamed and Toksoz, 1984). Observations 3–7 confirm quantitatively the expectation that the dense atmosphere of Venus has strongly affected the production of craters. Large impactors have been relatively unaffected, intermediate‐sized ones have been fragmented and have produced overlapping or multiple craters, a narrow size range has produced shock‐induced “splotches” but no craters, and the smallest bodies have had no observable effect on the surface. The number of craters eliminated by the “atmospheric filter” is enormous, about 98% of the craters between 2 and 35 km in diameter that Magellan might have observed on a hypothetical airless Venus. Unique crater‐related features such as parabolas and outflow deposits demonstrate the roles of Venus' high atmospheric density and temperature in modifying the crater formation process. Finally, heavily fractured craters and lava‐embayed craters are found to have higher than average densities along the major fracture belts and rifted uplands connecting Aphrodite Terra and Atla, Beta, Themis, and Phoebe regiones. These craters thus provide physical evidence for recent volcanic and tectonic activity at a low level.

Journal of Geophysical Research E: Planets↗
USGS2012

Enhanced visualization for the interpretation of Magellan radar data: Supplement to the Magellan special issue

The differences of radar data from more familiar photographic images, coupled with the unique geologic characteristics of Venus's surface, can make interpretation of the standard Magellan data products difficult for both planetary scientists and nonspecialists. We describe a set of digital processing techniques for transforming individual remote‐sensing datasets in order to make the information they contain more apparent, and for combining multiple datasets of diverse resolution and information content into synthetic images that allow the viewer to explore correlations among the data. The creation of derived data products by techniques like those we describe is becoming increasingly widespread and important, both for scientific research and for use in educational and media presentations. We apply our processing techniques to Magellan synthetic aperture radar images, altimetry, and microwave emissivity measurements of four areas that exemplify some of the most important styles of geologic activity on Venus.

Journal of Geophysical Research E: Planets↗
USGS2012

Effects of hydrothermal alteration on the magnetization of the Oligocene Carpenter Ridge Tuff, Bachelor Caldera, San Juan Mountains, Colorado

Intracaldera Oligocene Carpenter Ridge Tuff fills the Bachelor caldera in the central San Juan caldera complex and hosts mineral deposits of the Creede mineral district. The Carpenter Ridge Tuff and unaltered portions of its intracaldera Bachelor Mountain Member, have strong, high‐coercivity, reverse magnetizations with average magnetic susceptibility (MS) and natural remanent magnetization (NRM) of 6 × 10 −3 volume SI and 8 A/m, respectively. Oxide phenocrysts in these rocks are titanomagnetite and ilmenite; however, magnetization appears to be controlled by microcrystic titanomaghemite based on thermal demagnetization unblocking temperatures and Curie temperatures between 580°C and 620°C. Much of the intracaldera tuff was affected by potassic metasomatism, a type of hydrothermal alteration characterized by addition of K and loss of Ca and Na, between 27.3 Ma and 25.1 Ma. Potassic metasomatism resulted in the replacement of original feldspars by potassium feldspar and quartz, oxidation of the original oxide phenocrysts to hematite and rutile, and consequent suppression of MS and NRM by a factor of 5. Also present, however, are metasomatized rocks that have high magnetizations even though their original oxide phenocrysts were destroyed; values of MS and NRM are similar to those of outflow tuff. Such rocks are suspected of containing secondary magnetite on the basis of (1) Curie and thermal unblocking temperatures at and below 580°C, indicating that magnetization is carried by magnetite but that primary microcrysts of maghemite were destroyed, and (2) observation in one sample of magnetite rimming cores of hematite+rutile. The presence of secondary magnetite in metasomatized rocks requires a local shift to more reducing fluid chemistry, although we are unable to determine whether this shift occurred during potassic metasomatism or during later alteration or ore deposition. These altered rocks may record a change from a pervasive alteration by alkaline, oxidizing fluids to a vein‐controlled alteration dominated by more reduced fluids in which magnetite was stable.

Colorado↗
USGS2012

Broad trends in geomagnetic paleointensity on Hawaii during Holocene time

Paleointensity determinations have been obtained from 22 basaltic lava flows on the island of Hawaii using the Thelliers' method. Radiocarbon dating indicates that these flows erupted at intervals ranging from about 200 to 1000 years, and results of the experiments provide an estimate of broad trends in geomagnetic paleointensity during Holocene time in the vicinity of Hawaii. Most of the samples were obtained from quickly cooled flow margins and, as a consequence, typically contain two titanomagnetite populations. Only two of the 79 samples analyzed failed to yield an estimate of the paleofield, demonstrating that paleointensities can be obtained from such samples if carefully selected. Virtual dipole moments calculated for the flows are compared with a published curve of dipole field intensity that was constructed using worldwide archeomagnetic data. The large nondipole fields that were previously postulated for the vicinity of Hawaii are confirmed, and the present data indicate that they were present from about 5000 years B.P. to perhaps as recently as the past 200 years. The data indicate, however, that these nondipole sources must have been virtually absent between about 12,000 and 5000 years B.P. as they are at the present time.

Hawaii↗
USGS2012

Magnetic grain-size variations through an ash flow sheet: influence on magnetic properties and implications for cooling history

Rock magnetic studies of tuffs are essential to the interpretation of paleomagnetic data derived from such rocks, provide a basis for interpretation of aeromagnetic data over volcanic terranes, and yield insights into the depositional and cooling histories of ash flow sheets. A rhyolitic ash flow sheet, the Miocene-aged Tiva Canyon Member of the Paintbrush Tuff, contains both titanomagnetite phenocrysts, present in the magma prior to eruption, and cubic Fe-oxide microcrystals that grew after emplacement. Systematic variations in the quantity and magnetic grain size of the microcrystals produce large variations in magnetic properties through a section of the ash flow sheet penetrated in a borehole on the Nevada Test Site. Natural remanent magnetization varies from less than 1 × 10 −4 to more than 8 × 10 −4 A m 3 kg −1 , and in-phase magnetic susceptibility varies from less than 1 × 10 −6 to more than 10 × 10 −6 m 3 kg −1 . The microcrystals, which include both magnetite and maghemite, have Curie points and maximum unblocking temperatures between 580°C and 640°C. Rock magnetic data, including in-phase and quadrature magnetic susceptibilities as well as hysteresis parameters, demonstrate that these microcrystals are of superparamagnetic and single-domain sizes. Titanomagnetite phenocrysts are the dominant remanence carriers in the central 50 m of the section, whereas microcrystals are important contributors to remanent magnetization and magnetic susceptibility in two 15-m-thick zones at the top and bottom. Within these zones the size of microcrystals decreases both toward the quenched margins and toward the interior of the sheet. The decrease in microcrystal size toward the interior of the sheet is interpreted to indicate the presence of a cooling break; possibly represented by a concentration of pumice. Laboratory heating produces changes in magnetic properties which vary complexly from sample to sample. Such changes include (1) decrease in saturation magnetization, (2) both increase and decrease in in-phase and quadrature magnetic susceptibilities, and (3) shifts in coercivity spectra. The complex sample-to-sample variation in response to heating is attributed to differences in initial grain-size distributions.

Journal of Geophysical Research↗
USGS2012

Predictive modeling of the seismic cycle of the greater San Francisco Bay region

The seismic cycle for the San Francisco Bay region is synthesized by a model combining the pre- and post-1906 seismic histories. The long-term acceleration of seismic release (seismic moment, Benioff strain release, or event count) in the seismic cycle and the shorter-term accelerations preceding the larger earthquakes within that cycle are modeled using an empirical predictive technique, called a time-to-failure analysis, in which rate of seismic release is proportional to an inverse power of the remainng time to failure. The exponent of time to failure in the accelerating sequences appears to be scale invariant, and the length of the full cycle is estimated at 269±50 yr. The 1989 Loma Prieta earthquake should have been predictable with an uncertainty of two years in time and 0.5 in magnitude, although the specific location (at Loma Prieta) was not predictable by this technique. If the model is correct, the San Francisco Bay region should be entering a relatively long (20-50 yr) period of seismic quiescence above magnitude 6.

California↗
USGS2012

An episode of reinflation of the Long Valley Caldera, eastern California: 1989-1991

Following the episodes of inflation of the resurgent dome associated with the May 1980 earthquake sequence (four M 6 earthquakes) and the January 1983 earthquake swarm (two M 5.2 events), 7 years of frequently repeated two-color geodimeter measurements spanning the Long Valley caldera document gradually decreasing extensional strain rates from 5 ppm/yr in mid-1983, when the measurements began, to near zero in mid-1989. The corresponding seismic activity within the caldera persisted at a low rate of fewer than 10 M ≥ 1.2 earthquakes per week from 1985 through November 1989 with no events exceeding M 3.0. Early October 1989 marked a change in activity when measurements of the two-color geodimeter network showed a significant increase in extensional strain rate (9 ppm/yr) across the caldera. The seismic activity began exceeding 10 M ≥ 1.2 per week in early December 1989 and rapidly increased to a sustained level of tens of M ≥ 1.2 per week with bursts having hundreds of events per day. Many events exceeded M 3.0 and the largest event was M ≈ 4. The 1989–1991 inflation episode is the first time that we have sufficient geodetic measurements in Long Valley to define the temporal relation between onset of an inflation episode and onset of brittle failure (earthquake swarm within the caldera). Here, the onset of deformation preceded the onset of increased earthquake activity by more than 2 months. The seismicity rate began to decrease in mid-July 1990, consistent with a gradually slowing of extension across the caldera as measured by the two-color geodimeter. The recent episode of inflation can be modeled by a single Mogi point source located about 7 km beneath the center of the resurgent dome. In contrast, the deformation pattern observed between mid-1983 and mid-1989 is best reproduced by fault slip in the south moat, inflation at 6.5 km depth near Casa Diablo Hot Springs and inflation beneath the resurgent dome. It appears that the 7-km source beneath the resurgent dome that was active for the earlier episodes is the primary source for the more recent episode. The model used to satisfy the line length observations predicts 7.5 cm of uplift along leveling route along highway 395 from mid-1983 to mid-1989 and an additional 11 cm through the end of 1991. To comp are with the energy release from seismicity, the modeled inflation from late 1989 through the end of 1991 has a moment that is a factor of 40 more than the cumulative seismic moment from earthquakes located within the caldera from the same period. Thus the recent inflation episode represents a significant portion of the observed geodetic deformation with only little seismic energy release.

Journal of Geophysical Research↗
USGS2012

Helium isotope and gas discharge variations associated with crustal unrest in Long Valley Caldera, California, 1989-1992

The onset of anomalous seismic activity in 1989 beneath Mammoth Mountain on the southwestern rim of the Long Valley caldera, California, was followed within ∼4 months by a large increase in 3 He/ 4 He in vapor discharged from a fumarole on the north side of the mountain. The helium isotopic ratio at this vent rose to a maximum of 6.7 RA in July 1990 and subsequently declined to values near 5 R A . Potential sources of the 3 He‐rich vapors include degassing of fresh magma, degassing from fresh surfaces generated in newly fractured igneous rocks, and volatile release from a 3 He‐rich gas chamber situated above previously emplaced intrusives. The magnitude of the increase in helium isotopic composition (from 3.8 to 6.7 R A ), the persistence of relatively high values (>5 R A ) over a period of 3 years, the increase in the flux of total He relative to gases in air‐saturated water, and the increases in the rates of discharge of steam and gas from this fumarole indicate that magmatic intrusion did in fact begin in 1989 beneath Mammoth Mountain. Seismic activity and limited measurements of extensional deformation at the surface suggest that the depth of intrusion may be as shallow as 2 km, consistent with the prompt appearance of increased 3 He/ 4 He ratios in the fumarolic gas, and that the intrusive process may have persisted for ∼1 year. In contrast, a similar combination of magmatic intrusion and anomalous seismic activity beneath the resurgent dome‐south moat region during the 1989–1991 period resulted in at most relatively small changes in 3 He/ 4 He in fumarolic discharge at the southern edge of the resurgent dome. The more subdued response may result from a combination of greater intrusive depths and greater dilution of 3 He‐rich inputs to thermal fluid reservoirs in the shallow hydrothermal system in this area compared with Mammoth Mountain.

California↗
USGS2012

Evolution of the caldera‐forming eruption at Crater Lake, Oregon, indicated by component analysis of lithic fragments

Crater Lake caldera (8 × 10 km), formed 6845 years B. P. ( 14 C age) during the climactic eruption of the volcanic edifice known as Mount Mazama, is intermediate in size between small calderas associated with central vent eruptions and large calderas that have ring fracture vent systems. Our quantitative study of lithic fragments in the ejecta confirms the existing model of changes in vent configuration during the climactic eruption of Mount Mazama. Initial activity was from a single vent that produced a rhyodacite pumice fall from a Plinian column. Altered preexisting volcanic rocks are the predominant lithic type in the Plinian deposit, and their extensive hydrothermal alteration is considered as evidence of their relatively deep origin. The Wineglass Welded Tuff lies atop the Plinian deposit and contains a higher proportion of fresh volcanic rocks, suggesting enlargement of the single vent by slumping of its walls. This same vent enlargement caused the Plinian eruption column to collapse and feed valley‐hugging pyroclastic flows that deposited the Wineglass Welded Tuff. When enough material was erupted from the shallow magma chamber that its roof was no longer adequately supported, Mount Mazama collapsed to form the caldera, while highly energetic pyroclastic flows produced the climactic ignimbrite. A lag breccia that represents the proximal facies of the compositionally zoned climactic ignimbrite lies atop the Wineglass Welded Tuff and contains predominantly altered volcanic rocks of deeper origin, accompanied by minor granitoids from the magma chamber walls. Azimuthal differences in lithic component proportions in the lag breccia correlate well with the geology of the caldera walls, indicating that the climactic ignimbrite was ejected by multiple vents along a ring fracture system. Systematic lithic component changes within the lag breccia suggest different quarrying levels that reflect waxing and waning of the discharge rate during the caldera collapse phase of the climactic eruption. Our lithic component analysis demonstrates that calderas that may be too small to experience structural resurgence, such as Crater Lake, nevertheless may form by syneruptive subsidence along ring fractures.

Oregon↗
USGS2012

Three-dimensional velocity structure, seismicity, and fault structure in the Parkfield region, central California

This study examines the three-dimensional velocity structure in a 60- by 80-km region containing the Parkfield segment of the San Andreas fault. We use local earthquake and shot P arrival times in an iterative simultaneous inversion for velocity and hypocentral parameters. Using the three-dimensional model, we relocated 5251 events that occurred from 1969 to 1991, as well as the 1966 aftershocks, and computed 664 fault plane solutions. The San Andreas fault (SAF), characterized by a sharp across-fault velocity gradient, is the primary feature in the velocity solution. There is a 5–20% lateral change in velocity over a 4-km width, the contrast being sharper where there is better resolution. The model also shows significant variations in the velocity and in the complexity of the velocity patterns along the SAF. The largest across fault velocity difference is below Middle Mountain, where a large volume of low-velocity material impinges on the SAF from the northeast. This material is inferred to be overpressured and may be key to understanding the unusual behavior in the Parkfield preparation zone. A 20-km-long high-velocity slice is imaged northeast of the SAF near Gold Hill. Its along-fault length corresponds to the length of the maximum slip in 1966. The relocated seismicity shows that the San Andreas fault is a planar vertical fault zone at seismogenic depths. Ninety percent of the fault plane solutions that are on, or near, the SAF were right-lateral strike-slip on subvertical fault planes that parallel the SAF. Thus the surface fault complexities do not appear to extend to depth and therefore do not explain the rupture character at Parkfield. At Parkfield, variations in material properties play a key role in fault segmentation and deformation style. Our observations suggest that there may be a general relation between increasing velocity and increasing ability of the rocks to store strain energy and release it as brittle failure.

Journal of Geophysical Research↗
USGS2012

Precise U‐Pb ages of Duluth Complex and related mafic intrusions, northeastern Minnesota: Geochronological insights to physical, petrogenetic, paleomagnetic, and tectonomagmatic processes associated with the 1.1 Ga Midcontinent Rift System

Precise resolution of the timing of igneous activity is crucial to understanding the dynamic processes associated with continental rifting. Mafic intrusive rocks constitute a major portion of the exposed 1100 Ma (Keweenawan) Midcontinent Rift system in northeastern Minnesota; however, prior to this study, geochronological data were insufficient to allow rigorous interpretation of intrusive histories and their relationships to extrusive suites. Eight anorthositic and gabbroic intrusives were chosen to represent both the temporal and spatial ranges of plutonic activity that formed the Duluth Complex and related mafic intrusions. U‐Pb isotopic analyses from zircons and baddeleyites result in U‐Pb concordant ages with little or no ambiguity introduced by inherited components, Pb loss or common Pb. The earliest Keweenawan plutonism exposed in Minnesota occurs along the northeastern flank of the Duluth Complex as a series of layered gabbros (Nathan's layered series) emplaced at 1106.9 ± 0.6 Ma. This sequence of gabbro sheets shares temporal, spatial, and compositional similarities with the nearby Logan sills in Ontario. Four Duluth Complex anorthositic and troctolitic series samples from widely separated areas have unresolvable ages between 1099.3 ± 0.3 and 1098.6 ± 0.5 Ma, indicating a very short duration for peak intrusive activity (0.5–1 m.y.). The unresolvable ages between anorthositic and troctolitic plutons suggest that these two magma series are more closely related than previously modeled and argue strongly for the need to reexamine their fundamental petrogenetic relationships. These dates also imply that the major reverse‐to‐normal magnetic polarity switch, used throughout the rift system as an important correlation tool, occurred prior to 1099 Ma. This date is several million years earlier than previously suspected and emphasizes the need for further paleomagnetic and geochronological data from the overlying volcanics. Much of the hypabyssal intrusive suite within the volcanic pile overlying Duluth Complex plutons may be significantly younger than the main pulse of plutonic activity. Two hypabyssal bodies, the Sonju Lake intrusion and gabbro at Silver Bay, were emplaced at 1096.1 ± 0.8 Ma and 1095.8 ± 1.2 Ma, respectively. Dates reported here and in previous studies support the concept of episodic tectonomagmatic rift development where magmatism was apparently concentrated in episodes of short duration (<1–3 m.y.) interspersed with longer hiatuses (2–8 m.y.).

Michigan, Minnesota, Ontario, Wisconsin↗
USGS2012

Seismic evidence for active magmatic underplating beneath the Basin and Range Province, western United States

Near-vertical and wide-angle seismic reflection data provide evidence for the presence of a magma body at the base of the crust beneath Buena Vista Valley in northwestern Nevada. The seismic response of this hypothesized magma body is characterized by high-amplitude, near-vertical P wave reflections and a comparably strong P -to- S converted phase. The magma body, referred to here as the Buena Vista Magma Body, is probably a single sill with thickness no greater than 200 m and length no greater than 1.8 km. The melt fraction in the sill is undoubtedly greater than 20–30%, and probably exceeds 50%. Melt composition is unconstrained. Although the age of the Buena Vista Magma Body is difficult to determine precisely, it is probably no older than 500,000 years. This suggests that magmatism in the Basin and Range Province is an ongoing process, despite the relative paucity of volcanic rocks erupted at the surface during the last 6 m.y.

Nevada, Idaho↗
USGS2012

Deformation from 1973 through 1991 in the epicentral area of the 1992 Landers, California, Earthquake (Ms = 7.5)

Deformation of a 50 × 60 km trilateration network that spans the epicenter of the 1992 Landers earthquake( M s = 7.5) was measured by seven surveys over the 19 years preceding the earthquake. Three moderate earthquakes (1979 Homestead Valley, M s = 5.6; 1986 North Palm Springs, M s = 6.0; and 1992 Joshua Tree, M s = 6.1) occurred within the network during those 19 years. Here we use geodetic and seismic data to construct a dislocation model for each of the three moderate earthquakes. Coseismic changes due to these three moderate earthquakes as predicted by the dislocation models are then removed from the trilateration data. The residual geodetic changes appear to be uniform in time. We take those changes to represent secular strain accumulation. This strain accumulation clearly shows right‐lateral shear across the San Andreas fault but also suggests a northwest‐southeast extension northeast of the eastern end of the “Big Bend” of the San Andreas fault. The Landers earthquake ruptured along a north‐northwest trend across this region of northwest‐southeast extension.

California↗
USGS2012

The 1992 Landers earthquake sequence; seismological observations

The ( M W 6.1, 7.3, 6.2) 1992 Landers earthquakes began on April 23 with the M W 6.1 1992 Joshua Tree preshock and form the most substantial earthquake sequence to occur in California in the last 40 years. This sequence ruptured almost 100 km of both surficial and concealed faults and caused aftershocks over an area 100 km wide by 180 km long. The faulting was predominantly strike slip and three main events in the sequence had unilateral rupture to the north away from the San Andreas fault. The M W 6.1 Joshua Tree preshock at 33°N58′ and 116°W19′ on 0451 UT April 23 was preceded by a tightly clustered foreshock sequence ( M ≤4.6) beginning 2 hours before the mainshock and followed by a large aftershock sequence with more than 6000 aftershocks. The aftershocks extended along a northerly trend from about 10 km north of the San Andreas fault, northwest of Indio, to the east-striking Pinto Mountain fault. The M w 7.3 Landers mainshock occurred at 34°N13′ and 116°W26′ at 1158 UT, June 28, 1992, and was preceded for 12 hours by 25 small M ≤3 earthquakes at the mainshock epicenter. The distribution of more than 20,000 aftershocks, analyzed in this study, and short-period focal mechanisms illuminate a complex sequence of faulting. The aftershocks extend 60 km to the north of the mainshock epicenter along a system of at least five different surficial faults, and 40 km to the south, crossing the Pinto Mountain fault through the Joshua Tree aftershock zone towards the San Andreas fault near Indio. The rupture initiated in the depth range of 3–6 km, similar to previous M ∼5 earthquakes in the region, although the maximum depth of aftershocks is about 15 km. The mainshock focal mechanism showed right-lateral strike-slip faulting with a strike of N10°W on an almost vertical fault. The rupture formed an arclike zone well defined by both surficial faulting and aftershocks, with more westerly faulting to the north. This change in strike is accomplished by jumping across dilational jogs connecting surficial faults with strikes rotated progressively to the west. A 20-km-long linear cluster of aftershocks occurred 10–20 km north of Barstow, or 30–40 km north of the end of the mainshock rupture. The most prominent off-fault aftershock cluster occurred 30 km to the west of the Landers mainshock. The largest aftershock was within this cluster, the M w 6.2 Big Bear aftershock occurring at 34°N10′ and 116°W49′ at 1505 UT June 28. It exhibited left-lateral strike-slip faulting on a northeast striking and steeply dipping plane. The Big Bear aftershocks form a linear trend extending 20 km to the northeast with a scattered distribution to the north. The Landers mainshock occurred near the southernmost extent of the Eastern California Shear Zone, an 80-km-wide, more than 400-km-long zone of deformation. This zone extends into the Death Valley region and accommodates about 10 to 20% of the plate motion between the Pacific and North American plates. The Joshua Tree preshock, its aftershocks, and Landers aftershocks form a previously missing link that connects the Eastern California Shear Zone to the southern San Andreas fault.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Formation of left-lateral fractures within the Summit Ridge Shear Zone, 1989 Loma Prieta, California, Earthquake

The 1989 Loma Prieta, California, earthquake is characterized by the lack of major, throughgoing, coseismic, right-lateral faulting along strands of the San Andreas fault zone in the epicentral area. Instead, throughout the Summit Ridge area there are zones of tension cracks and left-lateral fracture zones oriented about N45°W, that is, roughly parallel to the San Andreas fault in this area. The left-lateral fractures zones are enigmatic because their left-lateral slip is opposite to the right-lateral sense of the relative motion between the Pacific and North American plates. We suggest that the enigmatic fractures can be understood if we assume that coseismic deformation was by right-lateral shear across a broad zone, about 0.5 km wide and 4 km long, beneath Summit Ridge. The shear zone has an orientation of about N70° W and is oblique to the San Andreas fault zone, which has a trend of N45° to 50°W, so that right-lateral shearing, along with some dilation, can account for the orientation of the fracture zones. Using an analog, kinematic model of the Summit Ridge shear zone and the observable geometric parameters (a shear zone about 0.5 km wide, an acute angle of 25° between the tension cracks and the shear-zone walls, left-lateral slip of 5 cm, and spacing of about 100 m in the tension cracks) we estimate that the amount of right-lateral shift across the Summit Ridge shear zone was on the order of 1.4 m. This is comparable to the amount of slip for coseismic faulting at depth, 1.9 to 2.4 m, as computed by several investigators. Thus contrary to most previous reports on the Loma Prieta earthquake, which assert that coseismic, right-lateral ground rupture was restricted to considerable (greater than 4 km) depths in the epicentral area, we find that nearly all the right-lateral offset is represented at the ground surface by the Summit Ridge shear zone.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Importance of model parameterization in finite fault inversions: Application to the 1974 Mw 8.0 Peru Earthquake

The spatial and temporal slip distributions for the October 3, 1974 ( M w = 8.0), Peru subduction zone earthquake and its largest aftershock on November 9 ( M s = 7.1) are calculated and analyzed in terms of the inversion parameterization and tectonic significance. Teleseismic, long-period World-Wide Standard Seismograph Network, P and SH waveforms are inverted to obtain the rupture histories. We demonstrate that erroneous results are obtained if a parameterization is used that does not allow for a sufficiently complex source, involving spatial variation in slip amplitude, risetime, and rupture time. The inversion method utilizes a parameterization of the fault that allows for a discretized source risetime and rupture time. Well-located aftershocks recorded on a local network have the same general pattern as teleseismically determined hypocenters and help to constrain the geometry of the subduction zone. For the main shock a hinged fault is preferred having a shallow plane with a dip of 11° and a deeper, landward plane with a dip of 30°. The preferred nucleation depth lies between 11 and 15 km. A bilateral rupture is obtained with two major concentrations of slip, one 60 to 70 km to the northwest of the epicenter and a second 80 to 100 km to the south and southeast of the epicenter. For these source regions, risetimes vary from 6 to 18 s. Our estimates of risetimes are consistent with the time for the rupture to traverse the dominant local asperity. The slip distribution for the November 9 aftershock falls within a conspicuous hole in the main shock rupture pattern, near the hypocenter of the main shock. The November 9 event has a simple risetime function with a duration of 2 s. Aftershocks recorded by the local network are shown to cluster near the hypocenter of the impending November 9 event and downdip from the largest main shock source region. Slip during the main shock is concentrated at shallow depths above 15 km and extends updip from the hypocenter to near the plate boundary at the trench axis. The large amount of slip at shallow depths is attributed to the absence of any significant accretionary wedge of sediments, and the relatively young age and high convergence rate of the subducted plate, which results in good seismic coupling near the trench axis.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Variations in magma supply rate at Kilauea Volcano, Hawaii

When an eruption of Kilauea lasts more than 4 months, so that a well-defined conduit has time to develop, magma moves freely through the volcano from a deep source to the eruptive site at a constant rate of 0.09 km 3 /yr. At other times, the magma supply rate to Kilauea, estimated from geodetic measurements of surface displacements, may be different. For example, after a large withdrawal of magma from the summit reservoir, such as during a rift zone eruption, the magma supply rate is high initially but then lessens and exponentially decays as the reservoir refills. Different episodes of refilling may have different average rates of magma supply. During four year-long episodes in the 1960s, the annual rate of refilling varied from 0.02 to 0.18 km 3 /yr, bracketing the sustained eruptive rate of 0.09 km 3 /yr. For decade-long or longer periods, our estimate of magma supply rate is based on long-term changes in eruptive rate. We use eruptive rate because after a few dozen eruptions the volume of magma that passes through the summit reservoir is much larger than the net change of volume of magma stored within Kilauea. The low eruptive rate of 0.009 km 3 /yr between 1840 and 1950, compared to an average eruptive rate of 0.05 km 3 /yr since 1950, suggests that the magma supply rate was lower between 1840 and 1950 than it has been since 1950. An obvious difference in activity before and since 1950 was the frequency of rift zone eruptions: eight rift zone eruptions occurred between 1840 and 1950, but more than 20 rift zone eruptions have occurred since 1950. The frequency of rift zone eruptions influences magma supply rate by suddenly lowering pressure of the summit magma reservoir, which feeds magma to rift zone eruptions. A temporary drop of reservoir pressure means a larger-than-normal pressure difference between the reservoir and a deeper source, so magma is forced to move upward into Kilauea at a faster rate.

Journal of Geophysical Research B: Solid Earth↗