USGS ScienceSearch

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 415 records · Page 23Linked to original sources

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

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

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

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

Crustal structure and composition of the southern Foothills Metamorphic Belt, Sierra Nevada, California, from seismic data

The Foothills Metamorphic Belt is an accreted terrane consisting of Paleozoic and Mesozoic metamorphic rocks that separates the Great Valley from the Sierra Nevada batholith in northern and central California. Until recently, the only available geophysical data for this area were reconnaissance refraction surveys, and gravity and magnetic data. New insights into the structure of the deep crust are provided by the interpretation of a seismic reflection profile (CC‐2), acquired in 1984 by the U.S. Geological Survey at the southern end of the Foothills Metamorphic Belt. Our interpretation is constrained by a new seismic velocity model derived from coincident microearthquake data. Earthquake hypocenters that occur at unusually great depths of 12 to 30 km make the data set particularly useful for obtaining deep crustal velocity information. The velocity model shows velocities of 5.2 to 6.3 km s −1 for the upper 12 km of the crust, and 6.7 to 6.8 km s −1 from 12 km to an estimated Moho at 32 km. The upper crustal velocities correspond to metamorphic rocks and serpentinites of the Foothills Metamorphic Belt as well as to diorites and granodiorites of the Sierra Nevada batholith, while the lower crustal velocities are interpreted to represent intermediate to mafic granulites. The majority of the earthquake hypocenters as well as a 6.7 km s −1 layer in the velocity model corresponds in depth to thick zones of west dipping midcrustal reflections that may represent major shear zones formed during the late Jurassic Nevadan orogeny or synbatholithic ductile shear zones that accommodated crustal extension associated with batholith intrusion. These reflections are truncated updip by an inferred subvertical contact that coincides with the western edge of the Sierra Nevada batholith and the southward trace of the Bear Mountains fault zone. The updip truncation of midcrustal shear zones and high lower crustal velocities indicate that strike‐slip faulting and magmatic underplating can be important processes during the docking and welding of an accreted terrane.

California

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

Fluid expulsion sites on the Cascadia accretionary prism: mapping diagenetic deposits with processed GLORIA imagery

Point-discharge fluid expulsion on accretionary prisms is commonly indicated by diagenetic deposition of calcium carbonate cements and gas hydrates in near-surface (<10 m below seafloor; mbsf) hemipelagic sediment. The contrasting clastic and diagenetic lithologies should be apparent in side scan images. However, sonar also responds to variations in bottom slope, so unprocessed images mix topographic and lithologic information. We have processed GLORIA imagery from the Oregon continental margin to remove topographic effects. A synthetic side scan image was created initially from Sea Beam bathymetric data and then was subtracted iteratively from the original GLORIA data until topographic features disappeared. The residual image contains high-amplitude backscattering that we attribute to diagenetic deposits associated with fluid discharge, based on submersible mapping, Ocean Drilling Program drilling, and collected samples. Diagenetic deposits are concentrated (1) near an out-of-sequence thrust fault on the second ridge landward of the base of the continental slope, (2) along zones characterized by deep-seated strikeslip faults that cut transversely across the margin, and (3) in undeformed Cascadia Basin deposits which overlie incipient thrust faults seaward of the toe of the prism. There is no evidence of diagenetic deposition associated with the frontal thrust that rises from the dècollement. If the dècollement is an important aquifer, apparently the fluids are passed either to the strike-slip faults which intersect the dècollement or to the incipient faults in Cascadia Basin for expulsion. Diagenetic deposits seaward of the prism toe probably consist dominantly of gas hydrates.

Cascadia Basin

Evolution of the Precambrian lithosphere: Seismological and geochemical constraints

Several recent models of crustal evolution are based on the belief that the thickness of the continental crust is proportional to its age, with ancient crust being the thickest. A worldwide review of seismic structure contradicts this belief and falsifies these models, at least for the Archean. Proterozoic crust has a thickness of 40–55 km and a substantial high‐velocity (>7 km/s) layer at its base, while Archean crust is only 27–40 km thick (except at the site of younger rifts and collisional boundaries) and lacks the basal high‐velocity layer. Seismology also provides evidence that the lithosphere is thickest beneath Archean cratons, while diamond ages show that this lithospheric keel must have already existed in the Archean. Geochemical data also indicate significant differences between Archean and Proterozoic lithosphere. Major and trace element studies of sediments show a change in upper crustal composition between the Archean and Proterozoic. Archean rocks are depleted in Si and K and enriched in Na, Ca, and Mg. There is also a marked change in the Eu/Eu * ratio. Mantle xenoliths and continental flood basalts show that the mantle lithosphere beneath Archean crust is ultradepleted in FeO compared to that beneath post‐Archean crust. The secular change in the crust‐forming process is attributed to a decline in mantle temperature, leading to a change in the composition of the lithospheric mantle. The higher temperature of the Archean mantle led to the eruption of komatiitic lavas, producing a refractory lithospheric mantle which is ultradepleted in FeO and volatiles. The resultant lithospheric keel is intrinsically less dense than the surrounding mantle and thus not susceptible to delamination. It was sufficiently thick and cool for diamonds to form during the Archean. In contrast, Proterozoic crust developed above fertile mantle. The eruption of continental flood basalts and underplating of basaltic sills is attributed to subsequent heating and partial melting of the lithospheric mantle. Consequently, Proterozoic crust is thickened and has a high‐velocity basal layer.

Journal of Geophysical Research B: Solid Earth

Seismic guided waves trapped in the fault zone of the Landers, California, earthquake of 1992

A mobile seismic array of seven stations was deployed at 11 sites along the fault trace of the M 7.4 Landers earthquake of June 28, 1992, with a maximum offset of 1 km from the trace. We found a distinct wave train with a relatively long period following the S waves that shows up only when both the stations and the events are close to the fault trace. This wave train is interpreted as a seismic guided wave trapped in a low‐velocity fault zone. To study the distribution of amplitude of the guided waves with distance from the fault trace and also their attenuation with travel distance along the fault zone, we eliminated source and recording site effects by the coda normalization method. The normalized amplitudes of guided waves show a spectral peak at 3–4 Hz, which decays sharply with distance from the fault trace. Spectral amplitudes at high frequencies (8–15 Hz) show an opposite trend, increasing with distance from the fault trace. The normalized amplitudes of guided waves at 3–4 Hz also show a systematic decrease with hypocentral distance along the fault zone, from which we infer an apparent Q of 50. In order to confirm the existence of the guided waves, a dense array of 31 stations was deployed at one of the 11 sites. The resultant records revealed unequivocal evidence for the existence of guided waves associated with the fault zone. By modeling the waveforms as S waves trapped in a low‐velocity waveguide sandwiched between two homogeneous half‐spaces with velocity V s = 3.0 km/s, we infer a waveguide width of about 180 m, a shear velocity of 2.0–2.2 km/s, and a Q of ∼50. Hypocenters of aftershocks with clear guided waves show a systematic distribution both laterally and with depth delineating the extent of the low‐velocity fault zone in three dimensions. We find that the zone extends to a depth of at least 10 km. This zone apparently continues to the south across the Pinto Mountain fault because guided waves are observed at stations north of the Pinto Mountain fault for earthquakes with epicenters south of it. On the other hand, the zone appears to be discontinuous at the fault bend located about 20 km north of the mainshock epicenter; guided waves were observed for stations and epicenters which are located on the same sides of the fault bend but not for those on the opposite sides.

California

Postseismic deformation following the 1989 (M = 7.1): Loma Prieta, California, earthquake

Postseismic deformation along a 90‐km profile bisecting the projected surface trace of the coseismic rupture of the 1989 Loma Prieta earthquake has been monitored by frequent GPS surveys for 3.3 years following the earthquake. In addition to the expected deformation associated with secular strain accumulation on the San Andreas and Calaveras faults, deformation associated with postseismic readjustment has been detected. Most of that deformation can be attributed to 1.5 m right‐lateral and 0.9 m reverse postseismic slip on a 5‐km‐wide downdip extension of the Loma Prieta rupture. In addition, there seems to be a 0.1 m postseismic collapse of the Loma Prieta rupture zone in the direction perpendicular to the plane of the rupture. The fault‐normal (N48°E) surface displacements plotted as a function of time exhibit a curvature suggesting a relaxation time of about 1.4 years. Similar plots of the fault‐parallel (N42°W) displacement components do not exhibit significant curvature. Presumably, the deformation shown in those plots is dominated by secular strain accumulation along the San Andreas and Calaveras faults rather than postseismic relaxation.

California

The seismic velocity structure of the Newfoundland Appalachian orogen

The deep structure of the Newfoundland Appalachian orogen is investigated by analyses of three intersecting seismic refraction/wide-angle reflection profiles which traverse the Gander and Dunnage zones or central mobile belt of Newfoundland. A simultaneous travel time inversion for velocity and interface was applied to the in-line seismic refraction/wide-angle reflection data and constrained by synthetic amplitude models. The results of the modeling procedure show a subhorizontally layered crust with upper crustal velocities ranging from 5.4 to 6.2 km/s, a midcrustal velocity of 6.25–6.35 km/s, and a lower crustal velocity of 6.7±0.2 km/s. The top of the lower crust is marked by a series of prominent reflections between 18 and 23 km depth which suggest a complex layered velocity interface. Strong laterally coherent Moho reflections indicate a sharp crust-mantle transition at 35 ± 3 km. The uppermost mantle has a velocity of 8.0±0.2 km/s, and a reflecting horizon at 55 km depth suggests an increase to velocities approaching 8.5 km/s. Normal moveout corrections applied to fan profiles provide constraining evidence for the reflecting horizon at the top of the lower crust and laterally continuous Moho reflections at 11–12 s two-way travel time. Comparisons with a coincident deep seismic reflection profile show that the refraction and reflection Mohos match to better than 2–3 km. Bulk Poisson's ratios of 0.23–0.24 for the whole crust calculated from P m P / S m S travel times suggest a crust dominated by quartzofeldspathic lithologies and a notable absence of voluminous mafic additions to the lower crust. The absence of a deep crustal root, coupled with the bulk intermediate composition inferred for the lower crust from the seismic refraction/wide-angle data, implies that the crust beneath central Newfoundland has undergone multiple periods of reactivation and equilibration following successive orogenic episodes.

Newfoundland and Labrador

Aeromagnetic evidence for a buried Early Cretaceous magmatic arc, northeast Japan

Positive aeromagnetic anomalies, recent drilling, and models constructed from these data delineate the plutonic roots of the Early Cretaceous Kitakami magmatic arc in northeast Japan. Buried plutons, mostly offshore, produce belts of positive magnetic anomalies. These anomalies and magnetotelluric data suggest that the plutons form a batholith 70–120 km wide, nearly 800 km long, and 10–15 km thick. The batholith may mark the location of the main Kitakami arc. Most of the exposed Kitakami plutons are 2–20 km in diameter; some are 3 km thick. The small plutons line up along NW trending faults; some may have been satellite vents that tapped into the magma supply of the main arc. The batholithic roots of the main arc now compose almost half of the modern Japan forearc basement. Steep magnetic gradients, offset anomalies, and basin stratigraphy portray extensive faulting of the Kitakami batholith during oblique subduction in the Late Cretaceous and rifting in the Miocene. The eastern boundary of the Kitakami batholith lies between 90 and 140 km west of the modern trench, much closer than the 300-km distance between the active arc and trench. The Early Cretaceous forearc basin and accretionary prism may underlie the modern forearc basin east of the batholith, but clear evidence is lacking. Much of the Early Cretaceous margin, including most of its forearc therefore is missing. How the material was removed is unknown: it could have been strike-slip faulted, eroded by subduction-related processes, or both.

Journal of Geophysical Research B: Solid Earth

Introduction to the special section northeast Japan: A case history of subduction

Subduction-related tectonic events such as strike-slip faulting, uplift of high-pressure metamorphic rocks in the forearc, back arc spreading, arc magmatism, and continental collisions have shaped northeast Japan, making it a case history for subduction. Much of the record of these events is preserved and can be used to reconstruct the tectonic history of the region. Although many geological, geophysical, and geochemical data have been collected in Japan for the last 30 years, interpretation of these data into a plate tectonic context did not really begin until about 10 years ago. Much of the data and interpretation on northeast Japan have been published in Japanese or as individual local studies. The primary goal of this special section is to present new investigations covering a wide range of topics and scales and place them in a tectonic framework. The papers in this special section focus on three main subjects: Mesozoic aspects of northeast Japan; the relation of Japan Sea opening with the tectonics and geochemistry of northeast Japan in the Tertiary; and the geophysical setting. The purpose of this overview is to link these subjects and provide a tectonic framework for north- east Japan and eastern Asia since the Mesozoic.

Journal of Geophysical Research B: Solid Earth

Strain accumulation along the Laguna Salada Fault, Baja California, Mexico

Strain accumulation observed over the 1978–1991 interval in a 30×100 km aperture trilateration network spanning the Laguna Salada fault is described by the principal strain rates 0.101±0.012 strain/yr N80°E±2° and −0.021±0.012 strain/yr N10°W±2°, extension reckoned positive. These strain accumulation rates have been corrected to remove coseismic effects of the nearby 1979 Imperial Valley ( = 6.5), 1980 Victoria (Baja California) ( = 6.4), 1987 Superstition Hills ( = 6.5), and 1987 Elmore Ranch ( = 5.9) earthquakes. The observed strain rates indicate extension at a rate of about 0.08 strain/yr perpendicular to the trend (N35°W) of the Salton trough as well as a right‐lateral tensor shear strain rate 0.05 strain/yr across it. The extension perpendicular to the trough is observed neither farther north near the Salton Sea nor farther south across the Gulf of California. However, Holocene slip on the Laguna Salada fault, about equal parts right‐lateral and normal slip, is consistent with the observed strain accumulation. A simple dislocation model intended to explain the observed strain accumulation as a product of slip at depth on the Laguna Salada fault would require that the fault be listric.

California

Three-dimensional modeling of pull-apart basins: implications for the tectonics of the Dead Sea Basin

We model the three-dimensional (3-D) crustal deformation in a deep pull-apart basin as a result of relative plate motion along a transform system and compare the results to the tectonics of the Dead Sea Basin. The brittle upper crust is modeled by a boundary element technique as an elastic block, broken by two en echelon semi-infinite vertical faults. The deformation is caused by a horizontal displacement that is imposed everywhere at the bottom of the block except in a stress-free &ldquo;shear zone&rdquo; in the vicinity of the fault zone. The bottom displacement represents the regional relative plate motion. Results show that the basin deformation depends critically on the width of the shear zone and on the amount of overlap between basin-bounding faults. As the width of the shear zone increases, the depth of the basin decreases, the rotation around a vertical axis near the fault tips decreases, and the basin shape (the distribution of subsidence normalized by the maximum subsidence) becomes broader. In contrast, two-dimensional plane stress modeling predicts a basin shape that is independent of the width of the shear zone. Our models also predict full-graben profiles within the overlapped region between bounding faults and half-graben shapes elsewhere. Increasing overlap also decreases uplift near the fault tips and rotation of blocks within the basin. We suggest that the observed structure of the Dead Sea Basin can be described by a 3-D model having a large overlap (more than 30 km) that probably increased as the basin evolved as a result of a stable shear motion that was distributed laterally over 20 to 40 km.

Dead Sea Basin

Interseismic uplift at the Nankai subduction zone, southwest Japan, 1951–1990

Uplift as a function of time from 1951 through 1990 has been deduced from annual mean sea level measured at 15 tide gages along the Nankai subduction zone, southwest Japan. The recurrence interval for rupture of the Nankai subduction zone is about 100 years, and the most recent rupture was in late 1946. Thus the 1951–1990 uplift record covers most of the first half of the earthquake cycle. The precision of the uplift record is better than could be obtained currently by annual Global Positioning System (GPS) measurements. The pre‐1960 uplift record shows rapid deformation that appears to terminate in 1959. After 1959 the uplift record is remarkably linear in time: Significant curvature in the uplift‐versus‐time plots is detected at only three of the 15 tide gages. The inferred uplift rates are not quantitatively consistent with the predictions of either the viscoelastic coupling or elastic half‐space models of subduction, but the agreement is sufficient to suggest that the causative physical mechanisms have been identified. The immediate postseismic response is controlled by the propagation of slip downward along the downdip extension of the coseismic rupture, and the interseismic response is caused by the accumulation of a slip deficit on the main thrust zone. Asthenosphere relaxation is not required to explain the observations.

Journal of Geophysical Research B: Solid Earth

Seismic velocity structure and composition of the continental crust: A global view

Seismic techniques provide the highest‐resolution measurements of the structure of the crust and have been conducted on a worldwide basis. We summarize the structure of the continental crust based on the results of seismic refraction profiles and infer crustal composition as a function of depth by comparing these results with high‐pressure laboratory measurements of seismic velocity for a wide range of rocks that are commonly found in the crust. The thickness and velocity structure of the crust are well correlated with tectonic province, with extended crust showing an average thickness of 30.5 km and orogens an average of 46.3 km. Shields and platforms have an average crustal thickness nearly equal to the global average. We have corrected for the nonuniform geographical distribution of seismic refraction profiles by estimating the global area of each major crustal type. The weighted average crustal thickness based on these values is 41.1 km. This value is 10% to 20% greater than previous estimates which underrepresented shields, platforms, and orogens. The average compressional wave velocity of the crust is 6.45 km/s, and the average velocity of the uppermost mantle ( Pn velocity) is 8.09 km/s. We summarize the velocity structure of the crust at 5‐km depth intervals, both in the form of histograms and as an average velocity‐depth curve, and compare these determinations with new measurements of compressional wave velocities and densities of over 3000 igneous and metamorphic rock cores made to confining pressures of 1 GPa. On the basis of petrographic studies and chemical analyses, the rocks have been classified into 29 groups. Average velocities, densities, and standard deviations are presented for each group at 5‐km depth intervals to crustal depths of 50 km along three different geotherms. This allows us to develop a model for the composition of the continental crust. Velocities in the upper continental crust are matched by velocities of a large number of lithologies, including many low‐grade metamorphic rocks and relatively silicic gneisses of amphibolite facies grade. In midcrustal regions, velocity gradients appear to originate from an increase in metamorphic grade, as well as a decrease in silica content. Tonalitic gneiss, granitic gneiss, and amphibolite are abundant midcrustal lithologies. Anisotropy due to preferred mineral orientation is likely to be significant in upper and midcrustal regions. The bulk of the lower continental crust is chemically equivalent to gabbro, with velocities in agreement with laboratory measurements of mafic granulite. Garnet becomes increasingly abundant with depth, and mafic garnet granulite is the dominant rock type immediately above the Mohorovicic discontinuity. Average compressional wave velocities of common crustal rock types show excellent correlations with density. The mean crustal density calculated from our model is 2830 kg/m 3 , and the average SiO 2 content is 61.8%.

Journal of Geophysical Research B: Solid Earth

Estimation of rod scale errors in geodetic leveling

Comparisons among repeated geodetic levelings have often been used for detecting and estimating residual rod scale errors in leveled heights. Individual rod-pair scale errors are estimated by a two-step procedure using a model based on either differences in heights, differences in section height differences, or differences in section tilts. It is shown that the estimated rod-pair scale errors derived from each model are identical only when the data are correctly weighted, and the mathematical correlations are accounted for in the model based on heights. Analyses based on simple regressions of changes in height versus height can easily lead to incorrect conclusions. We also show that the statistically estimated scale errors are not a simple function of height, height difference, or tilt. The models are valid only when terrain slope is constant over adjacent pairs of setups (i.e., smoothly varying terrain). In order to discriminate between rod scale errors and vertical displacements due to crustal motion, the individual rod-pairs should be used in more than one leveling, preferably in areas of contrasting tectonic activity. From an analysis of 37 separately calibrated rod-pairs used in 55 levelings in southern California, we found eight statistically significant coefficients that could be reasonably attributed to rod scale errors, only one of which was larger than the expected random error in the applied calibration-based scale correction. However, significant differences with other independent checks indicate that caution should be exercised before accepting these results as evidence of scale error. Further refinements of the technique are clearly needed if the results are to be routinely applied in practice.

Journal of Geophysical Research B: Solid Earth