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

Age of the Peach Springs Tuff, southeastern California and western Arizona

Sanidine separates from pumice of the early Miocene Peach Springs Tuff are concordantly dated at 18.5±0.2 Ma by two isotopic techniques. The Peach Springs Tuff is the only known unit that can be correlated between isolated outcrops of Miocene strata from the central Mojave Desert of southeastern California to the western Colorado Plateau in Arizona, across five structural provinces, a distance of 350 km. Thus the age of the Peach Springs Tuff is important to structural and paleogeographic reconstructions of a large region. Biotite and sanidine separates from bulk samples of the Peach Springs Tuff from zones of welding and vapor-phase alteration have not produced consistent ages by the K-Ar method. Published ages of mineral separates from 17 localities ranged from 16.2 to 20.5 Ma. Discordant 40 Ar/ 39 Ar incremental release spectra were obtained for one biotite and two of the sanidine separates. Ages that correspond to the last gas increments are as old as 27 Ma. The 40 Ar/ 39 Ar incremental release determinations on sanidine separated from blocks of Peach Springs Tuff pumice yield ages of 18.3±0.3 and 18.6±0.4 Ma. Laser fusion measurements yield a mean age of 18.51±0.10. The results suggest that sanidine and biotite K-Ar ages older than about 18.5 Ma are due to inherited Ar from pre-Tertiary contaminants, which likely were incorporated into the tuff during deposition. Sanidine K-Ar ages younger than 18 Ma probably indicate incomplete extraction of radiogenic 40 Ar, whereas laser fusion dates of biotite and hornblende younger than 18 Ma likely are due to postdepositional alteration. Laser fusion ages as high as 19.01 Ma on biotite grains from pumice suggest that minerals from pre-Tertiary country rocks also were incorporated in the magma chamber.

Journal of Geophysical Research Solid Earth↗
USGS2012

Crustal strain near the Big Bend of the San Andreas Fault: Analysis of the Los Padres-Tehachapi Trilateration Networks, California

In the region of the Los Padres-Tehachapi geodetic network, the San Andreas fault (SAF) changes its orientation by over 30° from N40°W, close to that predicted by plate motion for a transform boundary, to N73°W. The strain orientation near the SAF is consistent with right-lateral shear along the fault, with maximum shear rate of 0.38±0.01 μrad/yr at N63°W. In contrast, away from the SAF the strain orientations on both sides of the fault are consistent with the plate motion direction, with maximum shear rate of 0.19±0.01 μrad/yr at N44°W. The strain rate does not drop off rapidly away from the fault, and thus the area is fit by either a broad shear zone below the SAF or a single fault with a relatively deep locking depth. The fit to the line length data is poor for locking depth d less than 25 km. For d of 25 km a buried slip rate of 30 ± 6 mm/yr is estimated. We also estimated buried slip for models that included the Garlock and Big Pine faults, in addition to the SAF. Slip rates on other faults are poorly constrained by the Los Padres-Tehachapi network. The best fitting Garlock fault model had computed left-lateral slip rate of 11±2 mm/yr below 10 km. Buried left-lateral slip of 15±6 mm/yr on the Big Pine fault, within the Western Transverse Ranges, provides significant reduction in line length residuals; however, deformation there may be more complicated than a single vertical fault. A subhorizontal detachment on the southern side of the SAF cannot be well constrained by these data. We investigated the location of the SAF and found that a vertical fault below the surface trace fits the data much better than either a dipping fault or a fault zone located south of the surface trace.

Journal of Geophysical Research Solid Earth↗
USGS2012

Variations in fault slip and strain accumulation at Parkfield, California: Initial results using two-color geodimeter measurements 1984-1988

Repeated length measurements of several geodetic baselines near Parkfield, California, have revealed significant variations in the local rates of shallow fault slip and strain. This network of baselines, surveyed several times each week, straddles the San Andreas fault in the transition zone between the creeping section to the northwest and the locked section to the southeast. The length measurements, characterized by a precision approaching 0.1 ppm, reveal large fluctuations in the rates of baseline extension. Principal mode analysis of the length change data indicates that the two largest components of the signal are (1) secular extensions and contractions consistent with surface slip on the main strand of the San Andreas fault, and (2) a large seasonal oscillation with no obvious spatial coherence. On most of these baselines, the second component appears to be in phase with seasonal rainfall. When data from the baselines with the largest amplitude of the seasonal signal are excluded, the remaining data can be modeled in terms of both spatial and temporal variations in surface slip, variations in the components of the changes in uniform strain, and the possible displacement of the central monument in this radial network. In parameterizing this model, the spatial variation of slip beneath the near surface is reflected by changes in shear strain. Although the computed secular shear is highly dependent upon the specified parameterization of surface slip, the data are consistent with the hypothesis that slip at intermediate depths lags behind the surface slip rate. However, the range in models that fit the data does not necessarily imply that there is a deficit in slip at depth relative to the surface. Comparison of the inferred values of surface slip from the model with the observed fault slip measured by very short baseline creep meters indicates close agreement in secular rates, but the short-term variations observed with the creep meters are either highly attenuated or nonexistent in the modeled slip since the modeled slip is a spatial average which smooths out possible short-wavelength variations in the surface slip for which the creep instruments are most sensitive. An interesting conclusion from the two-color data is that surface slip on the San Andreas fault appears to be spread over a 2-km-wide zone on the south flank of Middle Mountain but is confined to a very narrow zone to the south as the fault passes through the center of the network. This conclusion is dependent upon the assumption that a few critical monuments are stable and track tectonic displacements in the long term. Finally, the largest observed strain change is an extensional strain coincident with the Kettleman Hills earthquake M 5.5 in August 1985.

Journal of Geophysical Research Solid Earth↗
USGS2012

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

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

Rupture process of a multiple main shock sequence: analysis of teleseismic, local and field observations of the Tennant Creek, Australia, earthquakes of January 22, 1988

On January 22, 1988, three large intraplate earthquakes (with M S 6.3, 6.4, and 6.7) occurred within a 12-hour period near Tennant Creek, Australia. These earthquakes, which occurred over a small interval of time and within a small volume of space, present a unique opportunity to study the rupture process of the class of intraplate earthquakes that occur as multiple main shocks. Broadband displacement and velocity records of body waves from teleseismically recorded data are analyzed to determine source mechanisms, depths, and complexity of rupture of each of the three main shocks. Hypocenters of an additional 150 foreshocks and aftershocks constrained by local arrival time data and field observations of surface rupture are used to complement the source characteristics of the main shocks in order to derive as complete a description of the rupture process as possible. The interpretation of the combined data sets suggests that the overall rupture process involved unusually complicated stress release. As locations of the main shock hypocenters progressively moved from west to east, we infer that the first and third main shocks, denoted as MS 1 and MS 3 , produced the southeast-northwest trending scarps observed at the western end (the Kunayungku fault) and at the eastern end (the east end of the Lake Surprise fault), respectively, of the rupture zone. The epicenter of the only immediate foreshock was located in the gap between these two fault scarps. MS 1 nucleated near this epicenter and ruptured upward and to the northwest from a depth of 6.5 km. MS 3 ruptured predominantly to the SE at a depth of 4.5 km. The second main shock, MS 2 , is inferred to have produced the deformation of the southwest trending central scarp segment (the western end of the Lake Surprise fault). From the sense of thrusting seen at the surface and from the distribution of aftershock hypocenters, the south dipping nodal planes derived from waveform modeling are identified as the fault planes for earthquakes MS1 and MS 3 . In contrast, the dip of the central fault scarp is reversed relative to the dips of the western and eastern fault scarps. The rupture process Of MS 2 turns out to be commensurately complex and sufficiently explains the geological complexity. MS 2 consisted of three subevents. The southeast dipping nodal plane of the first two subevents is coplanar with a southeast dipping plane implied by locations of aftershocks which did not break the surface. Choice of the north dipping plane as the rupture plane of the third subevent, consistent with the surface deformation and coplanar with a second plane delineated by aftershocks, would imply conjugate faulting. The majority of the aftershocks are concentrated near the edges of the fault planes, and there is an absence of activity in the center of the planes. The areas of absent activity may represent the failed asperities of the main shocks in which substantial stress relief occurred. The rupture process of each main shock is characterized by the rapid release of energy followed by a much slower release of moment and by aftershock zones whose dimensions exceed the inferred dimensions of the rupture. These characteristics suggest that substantial slow slip occurred on each of the three fault interfaces that was not accompanied by major energy release. The first main shock nucleated at the deepest part of its aftershock zone and ruptured upward. In contrast, MS 2 and MS 3 nucleated at depths that were one third to one half of the maximum depth of their aftershock zones. This variation of focal depth and the strong increase of moment and radiated energy with each main shock imply that lateral variations of strength were more important than vertical gradients of shear stress in controlling the progression of rupture.

Journal of Geophysical Research Solid Earth↗
USGS2012

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

Slip partitioning of the Calaveras Fault, California, and prospects for future earthquakes

Examination of main shock and microearthquake data from the Calaveras Fault during the last 20 years reveals that main shock hypocenters occur at depths of 8-9 km near the base of the zone of microearthquakes. Microseismicity extends between depths of 4 and 10 km and defines zones of concentrated microseismicity and aseismic zones. Estimates of the fault regions which slipped during the Coyote Lake and Morgan Hill earthquakes as derived from seismic radiation coincide with zones which are otherwise aseismic. We propose that these persistent aseismic zones represent stuck patches which slip only during moderate earthquakes. From the pattern of microearthquake locations we recognize six aseismic zones where we expect future main shocks will rupture the Calaveras Fault. -from Authors

Journal of Geophysical Research Solid Earth↗
USGS2012

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

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

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

Teleseismic tomography of the compressional wave velocity structure beneath the Long Valley region, California

In 1982 and 1984 the U.S. Geological Survey used several seismic networks, totaling over 90 stations, to record teleseismic P waves and measure travel time residuals in an area centered on the Long Valley caldera. We inverted the travel time residuals to obtain a three-dimensional image of the velocity structure with resolution of 5–6 km to depths of 70 km beneath the array. Direct inversion of these data indicates that the 2- to 4-km-thick low-velocity caldera fill contaminates the signal from any midcrustal velocity anomalies beneath the caldera. Thus two methods were used to strip the effects of the upper crust from the travel time residuals: (1) ray tracing through upper crustal velocity models provided by seismic refraction experiments and gravity surveys, and (2) an iterative stripping scheme using the inversion itself. The methods produce essentially identical results and adequately remove the effects of the shallowest crustal structures, including the caldera fill and hydrothermal alteration effects. The resulting “stripped” models show two well-resolved midcrustal low-velocity bodies in the Long Valley region. The first body is centered between 7 and 20 km depth beneath the resurgent dome of the Long Valley caldera and has a volume of 150–600 km 3 . The second, with a similar volume, is centered between 10 and 20 km depth beneath the Mono Craters, about 10 km north of Long Valley. Velocity contrasts in both of these bodies are about 6–10%, and the features are interpreted as silicic magma chambers. This experiment does not preclude the presence of additional pockets of magma smaller than 5 km across in the upper crust, particularly beneath the resurgent dome of the caldera (which would be removed with the stripping methods). It is likely that the midcrust anomaly beneath the resurgent dome is a remnant of the caldera forming magma chamber. Upper mantle velocities are lower than average beneath Mono Craters and higher than average beneath Long Valley. We tentatively interpret these anomalies as mantle partial melt beneath the Mono Craters and refractory upper mantle residuum beneath Long Valley. The high eruptive rate of the Mono Craters and these upper mantle structures suggest that the focus of volcanism is shifting north from Long Valley to the Mono Craters.

Journal of Geophysical Research Solid Earth↗
USGS2012

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

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

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

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

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

Paleomagnetic rotations and the Cenozoic tectonics of the Cascade Arc, Washington, Oregon, and California

Paleomagnetic results from Cenozoic (62–12 Ma) volcanic rocks of the Cascade arc and adjacent areas indicate that moderate to large clockwise rotations are an important component of the tectonic history of the arc. Two mechanisms of rotation are suggested by the regional pattern of paleomagnetic rotations. The progressive increase in rotation toward the coast in arc and forearc rocks results from distributed dextral shear, which is likely driven by oblique subduction of oceanic plates to the west. Simple shear rotation is accommodated in the upper crust by strike-slip faulting. The right-lateral Mount St. Helens seismic zone may be an active manifestation of this process. Dextral shear probably obscures a subequal contribution to arc and forearc rotation that is driven by intraarc or backarc extension. This rotation is suggested by the average southward increase in continental margin rotations into the region outboard of the Basin and Range. The southward increase in rotation parallels a change in the arc tectonic regime from largely compressional in northern Washington to extensional in Oregon. Concomitant with this change is a southward increase in the volume of eruptive rocks and the number of basaltic vents in the arc. A progressive eastward shift of the arc volcanic front with time in the rotated arc terrane is the result of the westward pivoting of the arc block in front of a zone of extension since Eocene time. Westward migration of bimodal Basin and Range volcanism since at least 16 Ma is tracking westward rotation of the frontal arc block and growth of the Basin and Range in its wake.

Journal of Geophysical Research Solid Earth↗