USGS Science⌕ Search

SEARCH · USGS Science

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 739 records · Page 41Linked to original sources
USGS2012

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

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

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

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

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

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

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

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

The July 12, 1993, Hokkaido-Nansei-Oki, Japan, earthquake: Coseismic slip pattern from strong-motion and teleseismic recordings

We employ a finite fault inversion scheme to infer the distribution of coseismic slip for the July 12, 1993 Hokkaido-Nansei-Oki earthquake using strong ground motions recorded by the Japan Meteorological Agency within 400 km of the epicenter and vertical P waveforms recorded by the Global Digital Seismograph Network at teleseismic distances. The assumed fault geometry is based on the location of the aftershock zone and comprises two fault segments with different orientations: a northern segment striking at N20°E with a 30° dip to the west and a southern segment with an N20°W strike. For the southern segment we use both westerly and easterly dip directions to test thrust orientations previously proposed for this portion of the fault. The variance reduction is greater using a shallow west dipping segment, suggesting that the direction of dip did not change as the rupture propagated south from the hypocenter. This indicates that the earthquake resulted from the shallow underthrusting of Hokkaido beneath the Sea of Japan. Static vertical movements predicted by the corresponding distribution of fault slip are consistent with the general pattern of surface deformation observed following the earthquake. Fault rupture in the northern segment accounts for about 60% of the total P wave seismic moment of 3.4 × 10 20 N m and includes a large circular slip zone (4-m peak) near the earthquake hypocenter at depths between 10 and 25 km. Slip in the southern segment is also predominantly shallower than 25 km, but the maximum coseismic displacements (2.0–2.5 m) are observed at a depth of about 5 km. This significant shallow slip in the southern portion of the rupture zone may have been responsible for the large tsunami that devastated the small offshore island of Okushiri. Localized shallow faulting near the island, however, may require a steep westerly dip to reconcile the measured values of ground subsidence.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Holocene paleoseismicity, temporal clustering, and probabilities of future large (M > 7) earthquakes on the Wasatch fault zone, Utah

The chronology of M > 7 paleoearthquakes on the central five segments of the Wasatch fault zone (WFZ) is one of the best dated in the world and contains 16 earthquakes in the past 5600 years with an average repeat time of 350 years. Repeat times for individual segments vary by a factor of 2, and range from about 1200 to 2600 years. Four of the central five segments ruptured between ∼620 ± 30 and 1230 ± 60 calendar years B.P. The remaining segment (Brigham City segment) has not ruptured in the past 2120 ± 100 years. Comparison of the WFZ space-time diagram of paleoearthquakes with synthetic paleoseismic histories indicates that the observed temporal clusters and gaps have about an equal probability (depending on model assumptions) of reflecting random coincidence as opposed to intersegment contagion. Regional seismicity suggests that for exposure times of 50 and 100 years, the probability for an earthquake of M > 7 anywhere within the Wasatch Front region, based on a Poisson model, is 0.16 and 0.30, respectively. A fault-specific WFZ model predicts 50 and 100 year probabilities for an M > 7 earthquake on the WFZ itself, based on a Poisson model, as 0.13 and 0.25, respectively. In contrast, segment-specific earthquake probabilities that assume quasi-periodic recurrence behavior on the Weber, Provo, and Nephi segments are less (0.01–0.07 in 100 years) than the regional or fault-specific estimates (0.25–0.30 in 100 years), due to the short elapsed times compared to average recurrence intervals on those segments. The Brigham City and Salt Lake City segments, however, have time-dependent probabilities that approach or exceed the regional and fault specific probabilities. For the Salt Lake City segment, these elevated probabilities are due to the elapsed time being approximately equal to the average late Holocene recurrence time. For the Brigham City segment, the elapsed time is significantly longer than the segment-specific late Holocene recurrence time.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

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

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

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

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

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

Mauna Loa lava accumulation rates at the Hilo drill site: Formation of lava deltas during a period of declining overall volcanic growth

Accumulation rates for lava flows erupted from Mauna Loa, as sampled in the uppermost 280 m of the Hilo drill hole, vary widely for short time intervals (several thousand years), but overall are broadly similar to those documented elsewhere on this volcano since 100 ka. Thickness variations and accumulation rates for Mauna Loa lavas at the Hilo drill site have been strongly affected by local paleotopography, including funneling and ponding between Mauna Kea and Kilauea. In addition, gentle submerged slopes of Mauna Kea in Hilo Bay have permitted large shoreline displacements by Mauna Loa flows. Ages of eruptive intervals have been determined from published isotopic data and from eustatic sea level curves modified to include the isostatic subsidence of the island of Hawaii at 2.2-2.6 mm/yr. Prior to 10 ka, rates of Mauna Loa lava accumulation at the drill site varied from 0.6 to 4.3 mm/yr for dateable intervals, with an overall rate of 1.8 mm/yr. Major eruptive pulses at about 1.3 and 10 ka, each probably representing a single long-lived eruption based on lack of weathering between flow units, increase the overall accumulation rate to 2.4 mm/yr. The higher rate since 10 ka reflects construction of thick near-shoreline lava deltas as postglacial sea levels rose rapidly. Large lava deltas form only along coastal segments where initially subaerial slopes have been submerged by the combined effects of eustatic sea level rise, isostatic subsidence, or spreading of volcano flanks. Overall accumulation of 239 m of lava at the drill site since 100-120 ka closely balances submergence of the Hilo area, suggesting that processes of coastal lava deposition have been modulated by rise in sea level. The Hilo accumulation rate is slightly higher than average rates of 1-2 mm/yr determined elsewhere along the Mauna Loa coast, based on rates of shoreline coverage and dated sea cliff and fault scarp exposures. Low rates of coastal lava accumulation since 100 ka, near or below the rate of island-wide isostatic subsidence, indicate that Mauna Loa is no longer growing vigorously or even maintaining its size above sea level.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Coral ages and island subsidence, Hilo drill hole

A 25.8-m-thick sedimentary section containing coral fragments occurs directly below a surface lava flow (the ∼1340 year old Panaewa lava flow) at the Hilo drill hole. Ten coral samples from this section dated by accelerator mass spectrometry (AMS) radiocarbon and five by thermal infrared multispectral scanner (TIMS) 230 Th/U methods show good agreement. The calcareous unit is 9790 years old at the bottom and 1690 years old at the top and was deposited in a shallow lagoon behind an actively growing reef. This sedimentary unit is underlain by a 34-m-thick lava flow which in turn overlies a thin volcaniclastic silt with coral fragments that yield a single 14 C date of 10,340 years. The age-depth relations of the dated samples can be compared with proposed eustatic sea level curves after allowance for island subsidence is taken. Island subsidence averages 2.2 mm/yr for the last 47 years based on measurements from a tide gage near the drill hole or 2.5–2.6 mm/yr for the last 500,000 years based on the ages and depths of a series of drowned coral reefs offshore from west Hawaii. The age-depth measurements of coral fragments are more consistent with eustatic sea levels as determined by coral dating at Barbados and Albrolhos Islands than those based on oxygen isotopic data from deep sea cores. The Panaewa lava flow entered a lagoon underlain by coral debris and covered the drill site with 30.9 m of lava of which 11 m was above sea level. This surface has now subsided to 4.2 m above sea level, but it demonstrates how a modern lava flow entering Hilo Bay would not only change the coastline but could extensively modify the offshore shelf.

Journal of Geophysical Research B: Solid Earth↗