Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research - Planets”

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

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

At least 1,621 records · Page 90Linked to original sources

Order and diversity in the modes of circum-Pacific earthquake recurrence

Recurrence characteristics of great circum-Pacific earthquakes and determinations of spatial distribution of seismic moment release are surveyed in order to delineate their general features and relate earthquake slip distribution to models of recurrent rupture. As noted by others, the pattern of moment release is typically very irregular, with strong concentrations in a few isolated regions of a much larger aftershock zone. Despite this complexity, rupture nucleation is notably systematic, with mainshock epicenter showing a strong tendency to locate in or immediately adjacent to identified regions of high moment release. This generalization suggests that earthquake recurrence is more likely to be controlled by maximum rather than average fault slip. Well-documented case histories from twelve plate boundary regions demonstrate that seismic strain release tends to be temporally well-ordered, while source dimensions, slip, and cumulative moment release vary considerably from cycle to cycle. On none of the segments studied do earthquakes consistently recur in nearly identical events. Instead, individual great earthquakes differ significantly from cycle to cycle or rerupture takes place in a sequence of two or more smaller events. Despite these differences the duration of the seismic cycle is approximately uniform. Furthermore, when strain release occurs in a sequence of large earthquakes, these events take place towards the end of the cycle and occur with increasing rupture length and magnitude. These ordered and irregular features of earthquake recurrence argue for the existence of corresponding elements on plate boundary faults. The ordered characteristics are identified with zones of concentrated moment release that slip comparable amounts in each cycle and have high shear strength. The irregular features are associated with intervening weaker regions that move in response to the stress concentration of dynamic rupture and slip by differing amounts in each event. Although these mechanistic associations are indirect and tentative, the features of recurrent behavior documented here have implications for long-term earthquake hazard assessment that are not dependent on the models proposed to explain them. Chief among these are the cycle-to-cycle differences among gap-filling events and the absence of shocks that fill major slip-deficient regions early in the seismic cycle.

Journal of Geophysical Research Solid Earth↗

Cooling rate and thermal structure determined from progressive magnetization of the dacite dome at Mount St. Helens, Washington

Our study of a magnetic anomaly associated with the recently active dacite dome at Mount St. Helens suggests that the dome consists of a hot, nonmagnetized core surrounded by a cool, magnetized carapace and flanking talus. The talus does not contribute to the anomaly because its constituent blocks are randomly oriented. Temporal changes in the magnetic anomaly indicate that the magnetized carapace thickened at an average rate of 0.03±0.01 m/d from 1984 to 1986. Petrographic and rock magnetic properties of dome samples indicate that the dominant process responsible for these changes is magnetization of extensively oxidized rock at progressively deeper levels within the dome as the rock cools through its blocking temperature, rather than subsequent changes in magnetization caused by further oxidation. Newly extruded material cools rapidly for a short period as heat is conducted outward in response to convective heat loss from its surface. The cooling rate gradually declines for several weeks, and thereafter the material cools at a relatively constant rate by convective heat loss from its interior along fractures that propagate inward. The rate of internal convective heat loss through fractures varies with rainfall, snowmelt, and large-scale fracturing during subsequent eruptive episodes. In accordance with a model for solidification of the 1959 lava lake at Kilauea Iki, Hawaii, we picture the dome's magnetized carapace as being a two-phase, porous, convective zone separated from the nonmagnetized core of the dome by a thin, single-phase conductive zone. As a consequence of the heat balance between the conductive and convective zones, the blocking-temperature isotherm migrates inward at a relatively constant rate. If the dome remains inactive, the time scale for its complete magnetization is estimated to be 18–36 years, a forecast which can be refined by shallow drilling into the dome and by continuing studies of its growing magnetic anomaly.

Journal of Geophysical Research Solid Earth↗

Earthquake depths and the relation to strain accumulation and stress near strike-slip faults in southern California

In this paper I report and attempt to interpret several observations about the depth distribution of earthquakes in the region of southern California containing the southern San Andreas and San Jacinto fault zones. These observations are as follows: (1) Earthquakes in the major fault zones are predominantly deep (maximum, 22 km; characteristic, 11–18 km). (2) Earthquakes in the crustal blocks bounding the fault zones are predominantly shallow (maximum, 17 km; characteristic, 1–6 km). (3) In the San Jacinto fault zone, maximum earthquake depths correlate with surface heat flow. These relations together with focal mechanisms, geodetic strain measurements, and fault zone models are consistent with the following ideas: (1) Interseismic plate motion is accommodated by aseismic slip along an extension of the major fault zone below a brittle zone that is locked between large earthquakes. (2) The aseismic slip in a narrow fault zone in the brittle-plastic transition region concentrates strain at the base of the brittle fault zone. (3) Deep earthquakes occur in the lower part of the brittle fault zone due to stick-slip failure of highly stressed patches. (4) Background earthquakes and aftershocks that occur several kilometers deeper than large earthquake hypocenters suggest that a zone of mixed slip behavior may exist between the stable sliding (deep) and stick-slip (shallow) regions of the fault zone. This zone of mixed slip may fail predominantly by stable sliding, analogous to the central San Andreas and Calaveras faults of California. Seismic precursors to large earthquakes may occur in this zone, and large earthquakes may nucleate near the upper edge of this zone. Furthermore, the difference in seismicity between the San Jacinto and southern San Andreas faults suggests that the nature of this mixed zone may evolve as total displacement in the fault zone increases. (5) Shear stress may be less in the crustal blocks than in the deep brittle fault zones and generally at a level sufficient to cause brittle failure only shallow in the crustal blocks. (6) In the stress field produced by plate motion and slip in the deep fault zone, the upper brittle fault zone is not oriented favorably for shear failure. Lack of shallow earthquakes in the fault zones and the predominance of shallow earthquakes on favorably oriented fractures in the adjacent crustal blocks suggest that either stress in the upper brittle fault zone is relatively low or the upper fault zone is effectively strong due to its orientation.

California↗

Age estimates and uplift rates for late Pleistocene marine terraces: Southern Oregon portion of the Cascadia forearc

Interest in the Cascadia subduction zone has increased because recent investigations have suggested that slip along plates at certain types of convergent margins is characteristically accompanied by large earthquakes. In addition, other investigations have suggested that convergent margins can be broadly classified by the magnitude of their uplift rates. The authors generated new uranium series, amino acid, and stable isotope data for southern Oregon marine terrace fossils. These data, along with terrace elevations and two alternative estimates of sea level at the time of terrace formation, allow one to determine terrace ages and uplift rates. Uranium series analysis of fossil coral yields an age of 83 {plus minus} 5 ka for the Whisky Run terrace at Coquille Point in Bandon, Oregon. A combination of amino acid and oxygen isotope data suggest ages of about 80 and 105 ka for the lowest two terraces at Cape Blanco. These ages indicate uplift rates of 0.45-1.05 and 0.81-1.49 m/kyr for Coquille Point and Cape Blanco, respectively. In order to assess the utility of the southern Oregon uplift rates for predicting the behavior of the Cascadia subduction zone, the authors compared late Quaternary uplift rates derived from terrace data from subduction zones around the world. On the basis of this comparison the southern Oregon rates of vertical deformation are not usually high or low. Furthermore, late Quaternary uplift rates show little relationship to the type of convergent margin. In the case of the southern Oregon coast, variability in uplift rate probably reflects local structures in the overriding plate, and the rate of uplift cannot be used as a simple index of the potential for great earthquakes along the southern Cascadia subduction zone.

California, Oregon↗

Reflected and mode-converted seismic waves within the shallow Aleutian subduction zone, southern Kenai Peninsula, Alaska

Pronounced secondary phases observed in local recordings of quarry shots and earthquakes on the southern Kenai Peninsula, Alaska, are identified as reflected P and S and converted S-to-P phases originating within four depth ranges: in the upper few kilometers of the Cook Inlet Tertiary basin, at midcrustal depths within the overthrust North American plate, at about 35 km depth near the top of the Wadati-Benioff seismic zone in proximity to the inferred interplate megathrust, and at about 5–10 km below the megathrust in the subducted Pacific plate. Beneath the landward margin of the accreted Chugach terrane, the mid-upper plate (MUP) discontinuity dips as steeply as 20°–30° to the west-northwest over the depth range 12–18 km. At shallower depths it flattens and possibly arches with the crest at about 10 km depth. Similar midcrustal reflectors have been imaged about 125 km to the southwest and about 350 km to the northeast along the structural trend of the Chugach terrane. The extensive reflectors may have a common origin, possible caused by the presence of underplated rocks. Relative amplitudes and polarities of the secondary phases originating at MUP depths provide few constraints on the nature of the discontinuity. It appears that the MUP discontinuity is seismically inactive and does not represent a brittle-ductile transition zone within the upper plate. The two converted S -to- P phases generated near the top of the subducted plate could indicate a low-velocity zone associated with subducted oceanic crust.

Journal of Geophysical Research Solid Earth↗

A tomographic glimpse of the upper mantle source of magmas of the Jemez lineament, New Mexico

The 800-km-long Jemez lineament is the most active volcanic feature in the southwestern United States. It is the southeastern tectonic boundary of the Colorado Plateau and crosses the Rio Grande rift at the Jemez Mountains. The primary volcanism of the lineament is basaltic and has occurred in the last 4.5 m.y. To infer spatial distributions of partial melt in the upper mantle source zones for the Rio Grande rift and the Jemez lineament, we investigated the lateral variations of P wave velocity in the upper mantle beneath these features. We used teleseismic P wave delays recorded at a 22-station network to perform a damped least squares, three-dimensional inversion for these lateral variations. Our technique employed velocity interpolation within a three-dimensional grid of points, rather than using blocks of constant P wave velocity. This method allows highly realistic computation of seismic ray paths as well as accurate computation of the matrix elements in our system of equations. Determinations of resolution of results were done in two independent ways, both of which gave consistent estimates of resolution. In our best resolved volume the inversion showed no significant concentration of relative low velocity for P waves beneath the Rio Grande rift. However, directly beneath the Jemez lineament there is a ∼100-km-wide, 1–2% low-velocity feature in the depth range of 50–160 km. Because of the association of the low P wave velocity with the Jemez volcanic lineament but not with the Rio Grande rift, because lowered P wave velocity can be associated with increased partial melt, and because the volume of recent volcanism at the lineament greatly exceeds that at the rift, we infer that a large magmatic source zone exists beneath the Jemez lineament but not beneath the Rio Grande rift. This implies that the volcanic potential of the Jemez lineaments continues to greatly exceed that of the Rio Grande rift. The mantle source zones for volcanics of the Jemez lineament are not overridden by, but rather track, the motion of the North American plate; this implies that these sources are within the lithospheric plate, as is clarified in the discussion. The magmatic source zones of the Jemez lineament are modeled as due to clockwise rotation of the Colorado Plateau about a pole in northeastern Colorado. This rotation caused extension of the lithosphere beneath the Jemez lineament, permitting concentration there of partially melted rock in the upper mantle.

Journal of Geophysical Research Solid Earth↗

Midcontinent rift volcanism in the Lake Superior region: Sr, Nd, and Pb isotopic evidence for a mantle plume origin

Between 1091 and 1098 Ma, most of a 15- to 20-km thickness of dominantly tholeiitic basalt erupted in the Midcontinent Rift System of the Lake Superior region, North America. The Portage Lake Volcanics in Michigan, which are the youngest MRS flood basalts, fall into distinctly high- and low-TiO 2 types having different liquid lines of descent. Incompatible trace elements in both types of tholeiites are enriched compared to depleted or primitive mantle (La/Yb = 4.3–5.3; Th/Ta = 2.12–2.16; Zr/Y = 4.3–4.4), and both basalt types are isotopically indistinguishable. Sr, Nd, and Pb isotopic compositions of the Portage Lake tholeiites have 87 Sr/ 86 Sr i ≈ 0.7038, ε Nd (1095 Ma) ≈ 0±2, and μ 1 ≈ 8.2. Model ages with respect to a depleted mantle source ( T DM ) average about 1950–2100 Ma. Portage Lake rhyolites fall into two groups. Type I rhyolites have Nd and Pb isotopic characteristics (ε Nd (1095 Ma) ≈ 0 to −4.7; μ 1 ≈ 8.2–7.8) consistent with contamination of tholeiitic rocks by 5–10% Archean crust. The one type II rhyolite analyzed has Nd and Pb isotopic compositions (ε Nd (1095 Ma) ≈ −13 to −16; μ 1 ≈ 7.6–7.7) which are consistent with partial melting of Archean crust. Early Proterozoic crust was not a major contaminant of MRS rocks in the Lake Superior region. Most reported Nd and Pb isotopic compositions of MRS tholeiites from the main stage of volcanism in the Lake Superior region and of the Duluth Complex are comparable to the Nd and Pb isotopic data for Portage Lake tholeiites. The isotopic enrichment of the MRS source compared to depleted mantle is striking and must have occurred at least 700 m.y. before 1100 Ma. There are two likely sources for such enriched MRS tholeiitic magmatism: subcontinental lithospheric mantle enriched during the early Proterozoic or enriched mantle derived from an upwelling plume. Subcontinental lithospheric mantle alone as a source cannot be ruled out strictly on isotopic grounds, but melting of enriched lithosphere potentially ranging in age from 3.6 to 1.8 Ga would have been less likely to generate either the homogeneity or the volume of MRS tholeiites. Decompression melting of an upwelling enriched mantle plume in a region of lithosphere thinned by extension could have successfully generated the enormous volume (850×10 3 km 3 ) of relatively homogeneous magma in a restricted time interval.

Journal of Geophysical Research Solid Earth↗

Keweenaw hot spot: Geophysical evidence for a 1.1 Ga mantle plume beneath the Midcontinent Rift System

The Proterozoic Midcontinent Rift System of North America is remarkably similar to Phanerozoic rifted continental margins and flood basalt provinces. Like the younger analogues, the volcanism within this older rift can be explained by decompression melting and rapid extrusion of igneous material during lithospheric extension above a broad, asthenospheric, thermal anomaly which we call the Keweenaw hot spot. Great Lakes International Multidisciplinary Program on Crustal Evolution seismic reflection profiles constrain end-member models of melt thickness and stretching factors, which yield an inferred mantle potential temperature of 1500°–1570°C during rifting. Combined gravity modeling and subsidence calculations are consistent with stretching factors that reached 3 or 4 before rifting ceased, and much of the lower crust beneath the rift consists of relatively high density intruded or underplated synrift igneous material. The isotopic signature of Keweenawan volcanic rocks, presented in a companion paper by Nicholson and Shirey (this issue), is consistent with our model of passive rifting above an asthenospheric mantle plume.

Midcontinent Rift System↗

Gabbroic xenoliths from the northern Gorda Ridge: Implications for magma chamber processes under slow spreading centers

Abundant gabbroic xenoliths in porphyritic pillow basalt were dredged from the northern Gorda Ridge. The host lava is a moderately fractionated, normal mid-ocean ridge basalt with a heterogeneous glass rind (Mg numbers 56-60). Other lavas in the vicinity range from near primary (Mg number 69) to fractionated (Mg number 56). On the basis of textures and mineral compositions, the xenoliths are divided into five types. The xenoliths are not cognate to the host lava, but they are genetically related. Chemistry of mineral phases in conjunction with textural features suggests that the xenoliths formed in different parts of a convecting magma chamber that underwent a period of closed system fractionation. The chamber was filled with a large proportion of crystalline mush when new, more primitive, and less dense magma was injected and mixed incompletely with the contents in the chamber, forming the hybrid host lava. -from Authors

Journal of Geophysical Research Solid Earth↗

Source complexity of the 1987 Whittier Narrows, California, earthquake from the inversion of strong motion records

Strong motion records for the Whittier Narrows earthquake are inverted to obtain the history of slip. Both constant rupture velocity models and variable rupture velocity models are considered. The results show a complex rupture process within a relatively small source volume, with at least four separate concentrations of slip. Two sources are associated with the hypocenter, the larger having a slip of 55–90 cm, depending on the rupture model. These sources have a radius of approximately 2–3 km and are ringed by a region of reduced slip. The aftershocks fall within this low slip annulus. Other sources with slips from 40 to 70 cm each ring the central source region and the aftershock pattern. All the sources are predominantly thrust, although some minor right-lateral strike-slip motion is seen. The overall dimensions of the Whittier earthquake from the strong motion inversions is 10 km long (along the strike) and 6 km wide (down the dip). The preferred dip is 30° and the preferred average rupture velocity is 2.5 km/s. Moment estimates range from 7.4 to 10.0×10 24 dyn cm, depending on the rupture model.

Journal of Geophysical Research Solid Earth↗

High spectral resolution reflectance spectroscopy of minerals

The reflectance spectra of minerals are studied as a function of spectral resolution in the range from 0.2 to 3.0 μm. Selected absorption bands were studied at resolving powers (λ/Δλ) as high as 2240. At resolving powers of approximately 1000, many OH‐bearing minerals show diagnostic sharp absorptions at the resolution limit. At low resolution, some minerals may not be distinguishable, but as the resolution is increased, most can be easily identified. As the resolution is increased, many minerals show fine structure, particularly in the OH‐stretching overtone region near 1.4 μm. The fine structure can enhance the ability to discriminate between minerals, and in some cases the fine structure can be used to determine elemental composition. For example, in amphiboles and talcs, four absorption bands are observed in the samples analyzed in this study that are due to hydroxyl linked to Mg 3 , Mg 2 Fe, MgFe 2 , and Fe 3 sites. The band intensities have been shown by other investigators to give the Fe:Fe+Mg ratio from transmission spectra. This study shows that the same equations can be used to obtain the ratio from reflectance spectra of unprepared samples. High‐resolution reflectance Spectroscopy of minerals may prove to be a very important tool in the laboratory, in the field using field‐portable spectrometers, from aircraft, and from satellites looking at Earth or other planetary surfaces.

Journal of Geophysical Research Solid Earth↗

Discontinuities in the shallow Martian crust at Lunae, Syria, and Sinai Plana

Detailed photoclinometric profiles across 125 erosional features and 141 grabens in the western equatorial region of Mars indicate the presence of three discontinuities within the shallow crust. Pits, troughs, and wall valleys (tributary canyons) within Noctis Labyrinthus and Valles Marineris and escarpments within the fretted terrain of Sacra Fossae and Kasei Valles show distinct erosional base levels at depths of 0.3–0.6 km, 1 km, and 2–3 km. The shallowest discontinuity corresponds to thickness estimates for the ridged plains unit in this region, and thus the discontinuity probably is the contact between a sequence of layered rock making up this unit and the underlying megaregolith. The 1-km discontinuity is reflected in the base levels of erosion of all the features studied, and it may correspond to the base of the proposed layer of ground ice. Model calculations using photoclinometric profiles of simple grabens (corrected for the effects of mass wasting) and the 60° dip of bounding faults (measured from the exposed traces of faults on trough walls) show that graben-bounding faults consistently intersect at the mechanical discontinuity at about 1 km depth. This discontinuity may represent an interface between ice-laden and dry regolith, ice-laden and water-laden regolith, or pristine and cemented regolith. A correlation between wall valley head depth and local thickness of the faulted layer suggests that the 1-km discontinuity also controlled the depth of the heads of sapping canyons. There is no apparent relation between the pit and trough depths and the local thickness of the faulted layer, which can be explained if pits and troughs developed by subsidence into underlying tension cracks that disrupted the 1-km discontinuity. The data do suggest, however, that erosion of shallow pits and troughs was influenced by the 1-km discontinuity. The third discontinuity, at a depth of 2–3 km, corresponds to the proposed base of the Martian megaregolith and is probably the interface between overlying, ejected breccia and in situ, fractured basement rocks.

Journal of Geophysical Research Solid Earth↗

The relation of catastrophic flooding of Mangala Valles, Mars, to faulting of Memnonia Fossae and Tharsis volcanism

Detailed stratigraphic relations indicate two coeval periods of catastrophic flooding and Tharsis centered faulting (producing Memnonia Fossae) in the Mangala Valles region of Mars. Major sequences of lava flows of the Tharsis Montes Formation and local, lobate plains flows were erupted during and between these channeling and faulting episodes. First, Late Hesperian channel development overlapped in time the Tharsis-centered faulting that trends north 75° to 90°E. Next, Late Hesperian/Early Amazonian flooding was coeval with faulting that trends north 55° to 70°E. In some reaches, resistant lava flows filled the early channels, resulting in inverted channel topography after the later flooding swept through. Both floods likely originated from the same graben, which probably was activated during each episode of faulting. Faulting broke through groundwater barriers and tapped confined aquifers in higher regions west and east of the point of discharge. The minimum volume of water required to erode Mangala Valles (about 5×10 12 m 3 ) may have been released through two floods that drained a few percent pore volume from a relatively permeable aquifer. The peak discharges of the floods may have lasted from days to weeks. The perched water discharged from the aquifer may have been produced by hydrothermal groundwater circulation induced by Tharsis magmatism, tectonic uplift centered at Tharsis Montes, and compaction of saturated crater ejecta due to loading by lava flows.

Journal of Geophysical Research Solid Earth↗

Flank tectonics of Martian volcanoes

On the flanks of Olympus Mons is a series of terraces, concentrically distributed around the caldera. Their morphology and location suggest that they could be thrust faults caused by compressional failure of the cone. In an attempt to understand the mechanism of faulting and the possible influences of the interior structure of Olympus Mons, we have constructed a numerical model for elastic stresses within a Martian volcano. In the absence of internal pressurization, the middle slopes of the cone are subjected to compressional stress, appropriate to the formation of thrust faults. These stresses for Olympus Mons are ∼250 MPa. If a vacant magma chamber is contained within the cone, the region of maximum compressional stress is extended toward the base of the cone. If the magma chamber is pressurized, extensional stresses occur at the summit and on the upper slopes of the cone. For a filled but unpressurized magma chamber, the observed positions of the faults agree well with the calculated region of high compressional stress. Three other volcanoes on Mars, Ascraeus Mons, Arsia Mons, and Pavonis Mons, possess similar terraces. Extending our analysis to other Martian volcanoes, we find that only these three and Olympus Mons have flank stresses that exceed the compressional failure strength of basalt, lending support to the view that the terraces on all four are thrust faults.

Journal of Geophysical Research B: Solid Earth↗

Dark materials in Valles Marineris: Indications of the style of volcanism and magmatism on Mars

Rifting on the equatorial canyon system of Valles Marineris provides a unique view of the interior of the Martian crust to depths reaching 7 km, exposing several in situ bedrock units which testify to past volcanic and magmatic processes on Mars. A thick, regionally extensive deposit observed in Coprates and Juventae chasmata is interpreted on the basis of spectral reflectance, erosional morphology, and tendency for eolian mobilization to be composed of mafic glass, possibly an ancient Martian analogue of the lunar terra mantling deposits. Spectral mapping suggests that the dark floor‐covering materials in the lower canyons are derived from this unit. A series of cliffs in the Ophir Chasma wall rock is interpreted to be exposures of resistant bedrock; the spectral signature of this massive and uniform unit most closely resembles that of terrestrial mafic rocks altered to or coated by crystalline hematite. Application of computer mapping techniques to probable young volcanic materials in the central troughs yields an inferred distribution of volcanic activity consistent with the interpretation of extrusion along faults near the margins of the canyon floors. This result supports the hypothesis that the valles originated through tectonic extension.

Journal of Geophysical Research B: Solid Earth↗

High-resolution reflectance spectra of Mars in the 2.3-μm region: evidence for the mineral scapolite

patially resolved reflectance spectra of Mars in the 2.2- to 2.4-μm spectral region were obtained in August 1988 using the NASA 3-m Infrared Telescope Facility. The spectra show weak absorption features due to Martian atmospheric carbon monoxide and a surface mineral. Both CO and the mineral absorptions are composed of overlapping narrow features, but in many locations, such as Hellas, Chryse, Eden, and Moab, the mineral absorptions are quite strong, at least 3 times stronger than at the most absorbing wavelengths of CO near 2.33 μm. Therefore CO complicates the analysis of the surface mineral but does not always overwhelm its signature. Model removal of the Martian atmospheric CO has been performed, and the remaining absorption bands are identified as scapolite. Relatively strong absorptions that match bands in the spectrum of scapolite and have little or no CO absorption interference are seen near 2.41, 2.39, and 2.29 μm. Absorption also occurs at the scapolite bands at 2.36 and 2.33 μm, but the analysis is complicated by uncertainty in the atmospheric CO removal at these wavelengths. Weaker scapolite bands are seen at 2.44 and 2.23 μm where there is virtually no atmospheric interference. The scapolite bands observed on Mars are due to HCO 3 − and HSO 4 − ions in the scapolite structure. The bicarbonate and bisulfate contents appear to vary with location: the scapolite in Hellas is more bisulfate-rich relative to that in the Chryse/Moab/Eden area. Other locations contain little (Arabia, Syrtis Major, Hellespontica, and Isidis) or no scapolite (e.g., Margaritifer, Ausonia, and Erythraeum). The calculated abundances are unconstrained because the amounts of HCO 3 − and HSO 4 − in the Martian scapolites as well as their grain sizes are not known. If the scapolites contain about 3 wt % of each, near the maximum possible, the scapolite abundances probably range from about 5 wt % scapolite at Eden and Hellas; 3–5% at Chryse, Moab, and Oxia Palus; 2–3% at Arabia, Syrtis Major, and Isidis; to less than 2% at Hellespontica, Syrtis Minor, and Margaritifer, assuming a relatively large grain size of 50–100 μm. If the characteristic grain sizes are smaller or the HCO 3 − and HSO 4 − contents are lower, the scapolite abundances required to match the observed band depths would be higher. The mineral bands are apparent in many of the Mars spectra measured, so it appears to be widely but not uniformly distributed. The newly observed fine structure also varies greatly in both depth and spectral detail with location on Mars. Thus there appears to be regional variations in composition. The mineral phases appear to reflect local or regional geology and are not primarily contained in the homogeneous, globally redistributed aeolian dust. Higher spectral resolution Martian spectra in the 2.3-μm region as well as at 3.9 μm are needed to confirm the scapolite identification and to constrain its abundance.

Journal of Geophysical Research↗

Three-dimensional P and S velocity structure in the Coalinga Region, California

The Coalinga earthquake sequence of 1983 provided a unique opportunity to perform a three-dimensional velocity and hypocenter inversion in an area of complex three-dimensional structure dominated by folding and blind thrusts. Additionally, since other varied geological and geophysical studies have been completed in this area, the three-dimensional inversion solution could be compared to other interpretations. Inversion of 7696 P and 1511 S first arrivals from earthquakes and 696 P first arrivals from refraction shots produced a three-dimensional velocity model with grid spacing of 1–2 km in the hypocentral area. The overall shape and location of velocity features correspond well to the mapped surface geology. The three-dimensional inversion yields details of folds where the resolution is good and the general shape where resolution is lower. The amounts of structural relief inferred for the local folds are similar to values inferred from geologic data. The three-dimensional velocity solution has several distinctive features. There is a linear high-velocity body (6.1–6.5 km/s), about 25 km long, from 6 to 8 km depth, that may represent a fragment of Coast Range ophiolite. A shallow low-velocity zone (LVZ), which extends for 20 km along the fold axis at about 6 km depth and correlates with LVZs observed in both the refraction and the reflection data, may indicate high pore pressure caused by lateral compressive strain. Deeper LVZs occur within inferred Franciscan material and are characterized by horizontal or southwest dipping zones of varied thickness, 4–8 km wide and 5–10 km long. These LVZs may represent multiple thrust faults. Their locations and geometry are consistent with thrust faults inferred with seismic reflection data. The three-dimensional velocity solution compares well to prior two-dimensional seismic reflection and refraction models and observed gravity. The shape of the inferred sedimentary section agrees well with the reflections from the Cenozoic strata. Compared to the refraction model, the three-dimensional solution has similar velocities and similar locations of velocity features but is more detailed in the hypocentral zone where it uses more data. The gravity computed from the three-dimensional velocities is similar to the observed gravity in both shape and amplitude. Both a simple one-dimensional initial model and a complex initial model derived from the refraction interpretation were tried. A simple starting model gave the best results. The S velocity solution has different resolution than the P velocity solution because it uses a different set of stations, and it has lower resolution because it uses fewer arrival times. While the general patterns of velocity variation are similar for both V p and V s , the V s solution tends to have more smearing of velocity features and can have somewhat different locations of velocity features.

Journal of Geophysical Research Solid Earth↗

Some additional observations on inclusions in the granitic rocks of the Sierra Nevada

Microgranular quartz diorite and diorite inclusions are widespread in central Sierra Nevada granitoid rocks and are almost exclusively restricted to hornblende-bearing rocks, most commonly felsic tonalites and mafic granodiorites. Strontium 87/strontium 86 values of the inclusions and host granitoids from individual plutons generally plot on single isochrons that agree closely with previously determined ages. Neodymium 143/neodymium 144 values of five inclusions ranged from 0.51229 to 0.51248 and yield ε Nd values that correlate closely with values of their respective hosts. The Nd-Sm and Rb-Sr systematics indicate that most inclusions were in isotopic equilibrium with enclosing materials at the time of formation. Silica contents of inclusions and granitoids are contiguous, but inclusions generally contain less than, and granitoids more than, 60% SiO 2 . Linear continuity of the femic oxide and H 2 O + trends relative to SiO 2 suggests many inclusions formed as concentrations of hydrous mafic minerals. Variation of other major element oxides and trace elements support this inference. The aforementioned data are generally inconsistent with a direct genetic link between the mafic inclusions and the mantle-derived basalt. However, many initial Sierra Nevada magmas were produced by mixing of mafic, mantle-derived melts and felsic crustal materials, resulting in isotopically varied tonalitic and granodioritic magmas. Most inclusions thus represent fragmented crystal accumulations of early-formed, near-liquidus minerals generated from these previously mixed magmas.

Journal of Geophysical Research Solid Earth↗