Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research, Solid Earth”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36Linked to original sources

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 “shear zone” 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↗

Rapid extension recorded by cooling‐age patterns and brittle deformation, Naxos, Greece

The metamorphic core complex exposed as the island of Naxos in the Aegean Sea records an unusually complete sequence of structures developed as a result of continental extension. The structures formed during Miocene rise and cooling from ductile, upper amphibolite facies and anatectic conditions to brittle near‐surface conditions beneath the Naxos detachment fault. Top‐to‐the‐north ductile fabrics in the footwall, which initially developed during amphibolite facies prograde metamorphism, were overprinted by a succession of north directed, normal sense lower‐temperature brittle structures as the footwall was tectonically unloaded and unroofed. Pseudotachylite and cataclasite formation, brittle faulting, alteration, and erosion of the footwall occurred during continued slip and tectonic denudation. Neogene conglomerate and megabreccia, in part derived from exhumation of this footwall, lie structurally above the peripheral Naxos fault. Published K‐Ar and 40 Ar/ 39 Ar ages for hornblende, white mica, and biotite in the footwall decrease northwestward; apparent ages 17–50 Ma in the southeast correspond to areas of low metamorphic grade where preextension argon was partially retained. Published ages 16–10 Ma in higher‐grade rocks of the domal core in the north are cooling ages that for each of the three minerals show a component of younging in the NNE direction of extension. Assuming this is the direction of unroofing, we interpret the rate of this younging as the fault slip rate as the footwall rocks moved >20 km SSW relative to their hanging wall along the base of the Naxos detachment fault. The calculated rates of slip average 5–8 mm/yr, comparable to maximum rates reported in the Basin and Range province.

Naxos↗

Effects of geothermal development on deformation in the Long Valley Caldera, eastern California, 1985-1994

Long Valley caldera in east central California has been the site of crustal unrest in the form of seismicity and ground deformation in response to magmatic inflation since 1980. Uplift of the resurgent dome has totaled ∼0.6 m over the 1975–1992 period. Within this region of uplift, and near the southwestern edge of the resurgent dome, is the Casa Diablo area, which experienced relative subsidence of 0.17 m between 1985 and 1992. Geothermal fluid production began at Casa Diablo in 1985 to supply a binary electric power plant; currently, three such plants generate a total of about 40 MW of electricity. The plants are supplied by 170°C water pumped from depths near 150m; all the produced water is reinjected at cooler temperatures at depths near 600 m. Analyses of data from (1) regional leveling lines and a network of bench marks in the geothermal well field area, (2) tilt observations from an L‐shaped array of bench marks, and (3) calculations relating pressure and temperature changes to subsidence indicate that relative subsidence at Casa Diablo results mainly from reductions in pressure in the shallow production reservoir and overlying formations and reductions in temperature in the underlying injection reservoir. Secondary effects on deformation are caused by offsets along the Taylor‐Bryant fault bounding the west side of the well field.

California↗

Identifying coseismic subsidence in tidal-wetland stratigraphic sequences at the Cascadia subduction zone of western North America

Tidal-wetland stratigraphy reveals that great plate boundary earthquakes have caused hundreds of kilometers of coast to subside at the Cascadia subduction zone. However, determining earthquake recurrence intervals and mapping the coastal extent of past great earthquake ruptures in this region are complicated by the effects of many sedimentologic, hydrographic, and oceanographic processes that occur on the coasts of tectonically passive as well as active continental margins. Tidal-wetland stratigraphy at many Cascadia estuaries differs little from that at similar sites on passive-margin coasts where stratigraphic sequences form through nonseismic processes unrelated to coseismic land level changes. Methods developed through study of similar stratigraphic sequences in Europe provide a framework for investigating the Cascadia estuarine record. Five kinds of criteria must be evaluated when inferring regional coastal subsidence due to great plate boundary earthquakes: the suddenness and amount of submergence, the lateral extent of submerged tidal-wetland soils, the coincidence of submergence with tsunami deposits, and the degree of synchroneity of submergence events at widely spaced sites. Evaluation of such criteria at the Cascadia subduction zone indicates regional coastal subsidence during at least two great earthquakes. Evidence for a coseismic origin remains equivocal, however, for the many peat-mud contacts in Cascadia stratigraphic sequences that lack (1) contrasts in lithology or fossils indicative of more than half a meter of submergence, (2) well-studied tsunami deposits, or (3) precise ages needed for regional correlation. Paleoecologic studies of fossil assemblages are particularly important in estimating the size of sudden sea level changes recorded by abrupt peat-mud contacts and in helping to distinguish erosional and gradually formed contacts from coseismic contacts. Reconstruction of a history of great earthquakes for the Cascadia subduction zone will require rigorous application of the above criteria and many detailed investigations.

Journal of Geophysical Research B: Solid Earth↗

Strain accumulation across the central Nevada seismic zone, 1973–1994

Five trilateration networks extending for 280 km along the central Nevada seismic zone (1915 Pleasant Valley, M = 7.3; 1954 Dixie Valley, M = 6.8; 1954 Stillwater, M = 6.8; 1954 Rainbow Mountain, M = 6.6; 1954 Fairview Peak, M = 7.1; and 1932 Cedar Mountain, M = 7.2) have been surveyed 6 times since 1973 to determine deformation along the zone. Within the precision of measurement the deformation appears uniform along the zone and is described by the principal strain rates 0.036±0.008 μstrain/yr N60°W±3° and −0.031±0.008 μstrain/yr N30°E±3°, extension reckoned positive. The observed strain rates are consistent with simple, right‐lateral, tensor shear at the rate of 0.033 μstrain/yr across a shear zone striking N15°W. This central Nevada shear zone appears to be the northward continuation of the eastern California shear zone. The orientation of the strike‐slip and normal‐slip ruptures within the central Nevada seismic zone are consistent with principal stress axes parallel to the measured principal strain rate axes. Space‐based geodetic measurements (very long baseline interferometry) indicate that the relative motion accommodated across the Basin and Range province west of Ely, Nevada, is about 9.1±1.5 mm/yr N16°W±8° (Dixon et al., 1995.) Notice that the right‐lateral shear zone postulated to explain deformation in the central Nevada seismic zone is properly oriented to accommodate that relative motion. However, a 135‐km effective width of the shear zone would be required to accommodate all of the 9.1 mm/yr relative motion at the strain rates observed in the Nevada seismic zone; only about 3 mm/yr of that relative motion is accommodated within the span of the trilateration networks.

California, Nevada↗

Tholeiitic‐alkalic transition at subglacial volcanoes, Tuya region, British Columbia, Canada

Ash Mountain, South Tuya, and Tuya Butte are three small basaltic volcanoes in the Stikine volcanic belt of northern British Columbia. The volcanoes rise 700, 500, and 400 m above their bases and are about 3.2, 1.6, and 2.6 km 3 in volume, respectively. They began eruptive activity under several hundred meters of overlying glacial ice, or water in an ice‐impounded lake, and undegassed pillow lava was erupted and forms the bases of all three. Later, as the vents grew into shallow water, explosive phreatomagmatic activity erupted partly degassed glassy tuffs. Finally, when the volcano emerged through the surface of the ice or water (or the water was drained), degassed subaerial lava flows were erupted and were converted to assemblages of foreset‐bedded pillow breccia and pillow lava when subaerial flows crossed a shoreline and flowed into meltwater lakes. The undegassed subglacial pillow base of Ash Mountain is overlain by partly degassed pillows and hyaloclastite tuff cut by dikes; at South Tuya the pillow base is overlain by hyaloclastite tuffs and lenses of pillow lava; at Tuya Butte the pillow base is overlain by foreset‐bedded pillow lava, pillow breccias, and hyaloclastite tuffs, which in turn are overlain by subaerial lava flows composing a small shield volcano. The undegassed basal subglacial pillow lava of the three volcanoes contain 0.10 ± 0.01 wt % sulfur and ∼0.5 wt % H 2 O. The overlying partly degassed assemblages contain 0.06 ± 0.02% sulfur and ∼0.2% H 2 O at Ash Mountain, 0.07±0.01% sulfur at South Tuya, and 0.03±0.01% sulfur at Tuya Butte. The differences in the degree of degassing can be related to the nature of eruption and quenching and the distance of flow of the subaerial lava. When the volcanoes switched from subglacial to shallow water or subaerial eruptions, as shown by change to more explosive activity and then to subaerial lava flows (and by a marked reduction of sulfur in volcanic glass), the magma shifted from tholeiitic to alkalic composition. This transition occurs at each of the three volcanoes. The tholeiitic and alkalic magmas cannot be related by shallow crystal fractionation and apparently originated by differing degrees of deep melting at a mantle source. Prior to eruption the tholeiitic melts overlay alkalic melts in shallow chambers underlying each of the volcanoes because of their lower density and were, therefore, the first to erupt under subglacial conditions. As the volcano grew through the ice (or ice‐impounded water), the volcanic conduit vented to the atmosphere, producing a partial depressurization of the conduit and the subsurface chamber. This sudden reduction in confining pressure caused enhanced vesiculation of volatile saturated melts, particularly of the more volatile‐rich alkalic melts, causing them to rise to the top of the chamber and erupt.

British Columbia↗

The 1992 M=7 Cape Mendocino, California, earthquake: Coseismic deformation at the south end of the Cascadia megathrust

We invert geodetic measurements of coseismic surface displacements to determine a dislocation model for the April 25, 1992, M = 7 Cape Mendocino, California, earthquake. The orientation of the model slip vector, which nearly parallels North America-Juan de Fuca relative plate convergence, and the location and orientation of the model fault relative to the offshore Cascadia megathrust, suggest that the 1992 Cape Mendocino earthquake is the first well-recorded event to relieve strain associated with the Cascadia subduction zone. We use data from three geodetic techniques: (1) the horizontal and vertical displacements of 13 monuments surveyed with the Global Positioning System, corrected for observed horizontal interseismic strain accumulation, (2) 88 section-elevation differences between leveling monuments, and (3) the uplift of 12 coastal sites observed from the die-off of intertidal marine organisms. Maximum observed displacements are 0.4 m of horizontal movement and 1.5 m of uplift along the coast. We use Monte Carlo techniques to estimate an optimal uniform slip rectangular fault geometry and its uncertainties. The optimal model using all the data resolves 4.9 m of slip on a 14 by 15 km fault that dips 28° SE. The fault extends from 1.5 to 8.7 km in depth and the main-shock hypocenter is close to the downdip projection of the fault. The shallowly dipping fault plane is consistent with the observed aftershock locations, and the estimated geodetic moment is 3.1 × 10 19 N m, 70% of the seismic moment. Other models that exclude leveling data collected in 1935 and 1942 are more consistent with seismological estimates of the fault geometry. If the earthquake is characteristic for this segment, the estimated horizontal slip vector compared with plate convergence rates suggests a recurrence interval of 140 years, with a 95% confidence range of 100–670 years. The coseismic uplift occurred in a region that also has high Quaternary uplift rates determined from marine terrace studies. If repeated ruptures of this southernmost segment of the Cascadia megathrust are responsible for the Quaternary uplift, a comparison of the coseismic uplift with coastal uplift rates suggests a recurrence interval of 200–400 years. Thus comparing horizontal and vertical coseismic to long-term deformation suggests a recurrence interval of about 100–300 years for M = 7 events at the south end of the Cascadia megathrust.

Journal of Geophysical Research B: Solid Earth↗

Equatorial origin for Lower Jurassic radiolarian chert in the Franciscan Complex, San Rafael Mountains, southern California

Lower Jurassic radiolarian chert sampled at two localities in the San Rafael Mountains of southern California (∼20 km north of Santa Barbara) contains four components of remanent magnetization. Components A , B ′, and B are inferred to represent uplift, Miocene volcanism, and subduction/accretion overprint magnetizations, respectively. The fourth component ( C ), isolated between 580° and 680°C, shows a magnetic polarity stratigraphy and is interpreted as a primary magnetization acquired by the chert during, or soon after, deposition. Both sequences are late Pliensbachian to middle Toarcian in age, and an average paleolatitude calculated from all tilt-corrected C components is 1° ± 3° north or south. This result is consistent with deposition of the cherts beneath the equatorial zone of high biologic productivity and is similar to initial paleolatitudes determined for chert blocks in northern California and Mexico. This result supports our model, in which deep-water Franciscan-type cherts were deposited on the Farallon plate as it moved eastward beneath the equatorial productivity high, were accreted to the continental margin at low paleolatitudes, and were subsequently distributed northward by strike-slip faulting associated with movements of the Kula, Farallon, and Pacific plates. Upper Cretaceous turbidites of the Cachuma Formation were sampled at Agua Caliente Canyon to determine a constraining paleolatitude for accretion of the Jurassic chert sequences. These apparently unaltered rocks, however, were found to be completely overprinted by the A component of magnetization. Similar in situ directions and demagnetization behaviors observed in samples of other Upper Cretaceous turbidite sequences in southern and Baja California imply that these rocks might also give unreliable results.

Journal of Geophysical Research B: Solid Earth↗

Recognizing and dating prehistoric liquefaction features: Lessons learned in the New Madrid seismic zone, central United States

The New Madrid seismic zone (NMSZ), which experienced severe liquefaction during the great New Madrid, Missouri, earthquakes of 1811 and 1812 as well as during several prehistoric earthquakes, is a superb laboratory for the study of world-class, earthquake-induced liquefaction features and their use in paleoseismology. In seismically active regions like the NMSZ, frequent large earthquakes can produce a complex record of liquefaction events that is difficult to interpret. Lessons learned studying liquefaction features in the NMSZ may help to unravel the paleoseismic record in other seismically active regions. Soil characteristics of liquefaction features, as well as their structural and stratigraphic relations to Native American occupation horizons and other cultural features, can help to distinguish prehistoric liquefaction features from historic features. In addition, analyses of artifact assemblages and botanical content of cultural horizons can help to narrow the age ranges of liquefaction features. Future research should focus on methods for defining source areas and estimating magnitudes of prehistoric earthquakes from liquefaction features. Also, new methods for dating liquefaction features are needed.

Journal of Geophysical Research B: Solid Earth↗

Thermal and unroofing history of a thick, tilted Basin-and-Range crustal section in the Tortilla Mountains, Arizona

We estimate here a geothermal gradient of only 17 ± 5°C km −1 for the tilted Grayback fault block in southeastern Arizona when extension began ∼25 Ma. This gradient is lower than preextension gradients estimated elsewhere in the Basin and Range, is only about 50% of typical gradients in the Basin and Range today, and needs to be accounted for in models of continental extension. The Grayback block exposes a 12-km-thick crustal section of Proterozoic and Cretaceous granitoids, which was tilted 90° during extension between 25 and 15 Ma. Zircon fission-track ages decrease structurally downward (westward) across the block and were all within a zone of partial track annealing prior to tilting and quenching. The zircon age gradient suggests that the 220°–240°C isotherm migrated downward 5–6 km during Paleogene erosion and regional cooling. Apatite fission-track ages decrease westward from ∼83 Ma in the structurally highest crystalline rocks to ∼24 Ma at ∼6-km paleodepth and then to ∼15 Ma another 6 km farther west. Track-length analysis confirms that apatites above the break in slope in age at ∼5.7-km paleodepth resided in a zone of partial annealing prior to tilting, and deeper apatites record rapid cooling upon tilting and unroofing beginning ∼25 Ma. At that time the 110 ± 10°C isotherm determined by the depth at which tracks in apatite were fully erased was at a basement paleodepth of ∼5.7 km, and the 220 ± 30°C isotherm as estimated from zircon data resided at a pretilting basement depth of ∼12.15 km. From consistent values of paleogeothermal gradient for two depth intervals we estimate the pretilt gradient was 17 ± 5°C km −1 . From 25 to 15 Ma the rotating Grayback block cooled rapidly as higher, westward moving blocks unroofed it tectonically at a rate of ∼1 km m.y. −1 .

Journal of Geophysical Research B: Solid Earth↗

Gravity field over the Sea of Galilee: Evidence for a composite basin along a transform fault

The Sea of Galilee (Lake Kinneret) is located at the northern portion of the Kinneret-Bet Shean basin, in the northern Dead Sea transform. Three hundred kilometers of continuous marine gravity data were collected in the lake and integrated with land gravity data to a distance of more than 20 km around the lake. Analyses of the gravity data resulted in a free-air anomaly map, a variable density Bouguer anomaly map, and a horizontal first derivative map of the Bouguer anomaly. These maps, together with gravity models of profiles across the lake and the area south of it, were used to infer the geometry of the basins in this region and the main faults of the transform system. The Sea of Galilee can be divided into two units. The southern half is a pull-apart that extends to the Kinarot Valley, south of the lake, whereas the northern half was formed by rotational opening and transverse normal faults. The deepest part of the basinal area is located well south of the deepest bathymetric depression. This implies that the northeastern part of the lake, where the bathymetry is the deepest, is a young feature that is actively subsiding now. The pull-apart basin is almost symmetrical in the southern part of the lake and in the Kinarot Valley south of the lake. This suggests that the basin here is bounded by strike-slip faults on both sides. The eastern boundary fault extends to the northern part of the lake, while the western fault does not cross the northern part. The main factor controlling the structural complexity of this area is the interaction of the Dead Sea transform with a subperpendicular fault system and rotated blocks.

Sea of Galilee↗

Three-dimensional P and S wave velocity structure of Redoubt Volcano, Alaska

The three‐dimensional P and S wave structure of Redoubt Volcano, Alaska, and the underlying crust to depths of 7–8 km is determined from 6219 P wave and 4008 S wave first‐arrival times recorded by a 30‐station seismograph network deployed on and around the volcano. First‐arrival times are calculated using a finite‐difference technique, which allows for flexible parameterization of the slowness model and easy inclusion of topography and source‐receiver geometry. The three‐dimensional P wave velocity structure and hypocenters are determined simultaneously, while the three‐dimensional S wave velocity model is determined using the relocated seismicity and an initial S wave velocity model derived from the P wave velocity model assuming an average Vp / Vs ratio of 1.78. Convergence is steady with approximately 73% and 52% reduction in P and S wave arrival time RMS, respectively, after 10 iterations. The most prominent feature observed in the three‐dimensional velocity models derived for both P and S waves is a relative low‐velocity, near‐vertical, pipelike structure approximately 1 km in diameter that extends from 1 to 6 km beneath sea level. This feature aligns axially with the bulk of seismicity and is interpreted as a highly fractured and altered zone encompassing a magma conduit. The velocity structure beneath the north flank of the volcano between depths of 1 and 6 km is characterized by large lateral velocity variations. High velocities within this region are interpreted as remnant dikes and sills and low velocities as regions along which magma migrates. No large low‐velocity body suggestive of a magma chamber is resolved in the upper 7–8 km of the crust.

Journal of Geophysical Research B: Solid Earth↗