USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Oceans”

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,459 records · Page 81Linked to original sources
USGS2012

Uranium-lead isotopic ages from the Sierra Nevada Batholith, California

This study provides new information on the timing and distribution of Mesozoic magmatic events in the Sierra Nevada batholithic complex chiefly between 36° and 37°N. latitude. U-Pb ages have been determined for 133 zircon and 7 sphene separates from 82 samples of granitoid rocks. Granitoid rocks in this area range in age from 217 to 80 m.y. Triassic intrusions are restricted to the east side of the batholith; Jurassic plutons occur south of the Triassic plutons east of the Sierra Nevada, as isolated masses within the Cretaceous batholith, and in the western foothills of the range; Cretaceous plutons form a continuous belt along the axis of the batholith and occur as isolated masses east of the Sierra Nevada. No granitic intrusions were emplaced for 37 m.y. east of the Sierra Nevada following the end of Jurassic plutonism. However, following emplacement of the eastern Jurassic granitoids, regional extension produced a fracture system at least 350 km long into which the dominantly mafic, calc-alkalic Independence dike swarm was intruded 148 m.y. ago. The dike fractures probably represents a period of regional crustal extension caused by a redistribution of the regional stress pattern accompanying the Nevadan orogeny. Intrusion of Cretaceous granitic plutons began in large volume about 120 m.y. ago in the western Sierra Nevada and migrated steadily eastward for 40 m.y. at a rate of 2.7 mm/y. This slow and constant migration indicates remarkably uniform conditions of subduction with perhaps downward migration of parent magma generation or a slight flattening of the subduction zone. Such steady conditions could be necessary for the production of large batholithic complexes such as the Sierra Nevada. The abrupt termination of plutonism 80 m.y. ago may have resulted from an increased rate of convergence of the American and eastern Pacific plates and dramatic flattening of the subduction zone. U-Pb ages of the Giant Forest-alaskite sequence in Sequoia National Park are all in the range 99±3 m.y., indicating a relatively short period of emplacement and cooling for this nested group of plutons. U-Pb ages of a mafic inclusion and its host granodiorite indicate that both were derived from a common source or that the mafic inclusion was totally equilibrated with the granodioritic magma. Comparison of isotopic ages determined by different methods such as zircon U-Pb, sphene U-Pb, hornblende K-Ar, and biotite K-Ar suggests that zircon U-Pb ages generally approximate the emplacement age of a pluton. However, some plutons probably contain inherited or entrained old zircons, and the zircons of some samples are disturbed by younger thermal and metamorphic events. The ages reported here are consistent with U-Pb age determinations previously made on granitic rocks to the north [Stern et al., 1981], The age distribution of granitic belts determined here is in general agreement with those established by K-Ar dating [Evernden and Kistler, 1970] but does not differentiate the five epochs of plutonism determined in their study.

California↗
USGS2012

Contemporary block tectonics: California and Nevada

Well-determined fault plane solution and the gross pattern of late-Cenozoic faulting in California and Nevada show a systematic relation between the orientation of fault planes and slip directions. In general, normal faults have northerly strikes, reverse faults have easterly strikes, and dextral and sinstral strike slip faults have northwesterly and northeasterly strikes, respectively. Kinematically, this relation is consistent with the response of clusters of fault-bounded crustal blocks to a regional stress field generated by the relative motion between the Pacific and North American plates. In this stress field, the greatest and least principal (compressive) stresses are restricted to northerly and easterly striking vertical planes, respectively. Clusters of crustal blocks bound by faults having the above attitudes form three basic regimes: (1) a spreading regime with a combination of normal and conjugate, strike slip faults, (2) a converging regime with a combination of reverse and conjugate, strike slip faults, and (3) a wrench regime with sets of subparallel, throughgoing, northwest striking (dextral) faults or northeast striking (sinstral) faults. These three regimes are typified by faulting patterns in the Basin and Range province (spreading), the Transverse Ranges (converging), and the San Andreas-Garlock fault systems (wrench), respectively. The gross deformation of each regime resulting from relative displacements between individual blocks is characterized by north-south shortening and east-west extension with the ratio of extensional to shortening strains (and the areal dilitation) decreasing systematically from spreading to wrench to compressional regimes. The wrench regime involves a component of net rotational deformation (clockwise for dextral slip and counter-clockwise for sinstral slip), while deformation of the spreading and converging regimes is irrotational. Local deviations from regional kinematic directions are concentrated along the boundaries between regimes reflecting the mismatch in gross deformation fields between regimes. Maximum principal and shear stress magnitudes will increase systematically from spreading to wrench to converging regimes provided that fault slip is controlled by frictional strength (Byerlee's law) along preexisting fractures and that pore pressure in the brittle crust is laterally uniform. A minimum strength difference between active, block-bounding faults and block interiors is 15–30%. Simple arrangements of such block clusters mimic the gross kinematic pattern of Quaternary faulting in California and Nevada. Some implications for contemporary tectonics emphasized by this model involve the westward displacement of the Sierra Nevada block.with respect to the stable interior of the North American plates, oblique thrusting of the Salinian block over the Pacific plate, and a progressive increase in the offset of the San Andreas fault represented by the ‘big bend’ through the Transverse Ranges.

Journal of Geophysical Research Solid Earth↗
USGS2012

The sea slope problem revisited

Discrepancies in sea surface topography based on comparisons between the results of steric leveling and repeated geodetic levelings have identified what is known as the ‘sea slope problem.’ This problem is actually twofold: (1) the sea surface relief based on steric leveling differs significantly from that based on geodetic leveling along several generally north-south coasts, and (2) successively propagated levelings between several widely separated tide stations indicate that the stationary sea slope seemingly has been changing with time, whereas differenced sea level means between these stations indicate that the sea surface relief has remained virtually invariant during the same intervals. Reexamination of the three reported discrepancies between geodetically and sterically determined sea slopes indicates that the Australian example is based on leveling of a quality inappropriate to the comparison. The discrepancy developed along the Atlantic coast of the United States is limited to the reach between Portsmouth–Hampton Roads, Virginia, and Key West, Florida, where the accuracy of steric leveling may be especially vulnerable owing to the dynamic effects of the Gulf Stream. Reconsideration of the example along the Pacific coast of the United States indicates that the various discrepancies are due chiefly to intrasurvey movement and resultant distortion of geodetically defined height differences between tide stations. Agreement between the results of steric and geodetic leveling along tectonically inactive north-south coasts devoid of strong boundary currents is generally good. This observation supports the conclusions (1) that any directionally dependent systematic error in geodetic leveling is measurably insignificant and (2) that where allowance is made for the possible effects of major boundary currents or intrasurvey movement during levelings between tide stations, the sea slope problem tends to vanish.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Earthquakes of Loihi submarine volcano and the Hawaiian hot spot

Loihi is an active submarine volcano located 35 km south of the island of Hawaii and may eventually grow to be the next and southernmost island in the Hawaiian chain. The Hawaiian Volcano Observatory recorded two major earthquake swarms located there in 1971–1972 and 1975 which were probably associated with submarine eruptions or intrusions. The swarms were located very close to Loihi's bathymetric summit, except for earthquakes during the second stage of the 1971–1972 swarm, which occurred well onto Loihi's southwest flank. The flank earthquakes appear to have been triggered by the preceding activity and possible rifting along Loihi's long axis, similar to the rift-flank relationship at Kilauea volcano. Other changes accompanied the shift in locations from Loihi's summit to its flank, including a shift from burst to continuous seismicity, a rise in maximum magnitude, a change from small earthquake clusters to a larger elongated zone, a drop in b value, and a presumed shift from concentrated volcanic stresses to a more diffuse tectonic stress on Loihi's flank. The 1971–1972 swarm began at depths of 20–50 km about 1 month before the shallow swarm started and suggests an upward migration of seismicity as the first stage of the 13-month swarm. The seismic ‘root’ of Kilauea volcano is well defined by earthquakes which plunge to the south and southwest to a depth of 50–60 km, terminating in a region in which earthquakes are associated with deep harmonic tremor and magma. The zone of these deeper earthquakes and tremor lies between Kilauea, Loihi, and Mauna Loa volcanoes, may feed magma to all three volcanoes, and probably locates the Hawaiian hot spot.

Journal of Geophysical Research Solid Earth↗
USGS2012

Rapid intrusion of magma into wet rock: Groundwater flow due to pore pressure increases

Analytical and numerical solutions are developed to simulate the pressurization, expansion, and flow of groundwater contained within saturated, intact host rocks subject to sudden heating from the planar surface of an igneous intrusion. For most rocks, water diffuses more rapidly than heat, assuring that groundwater is not heated along a constant-volume pressure path and that thermal expansion and pressurization adjacent to the intrusion drives a flow that extends well beyond the heated region. The forcing parameter for pressurization and flow is α Δ T , where α is a thermal expansion coefficient reflecting the overall expansion of water heated through the temperature difference Δ T between the initial ambient and intrusive values. Pore pressure increases due to heating are greatest when the intrusion is emplaced rapidly and where the intrusive contact is impervious to groundwater contained in stiff, impermeable rocks with high thermal diffusivities and porosities. The maximum velocity of water flowing in pores decays with the inverse square root of time and is insensitive to hydraulic properties of the host rocks. Pressures are lessened and flow directions are reversed with the onset of hydrothermal convection. This occurs at times ranging from hours to weeks after onset of intrusion. As magma rises into near-surface rocks, steam can be generated. Solutions indicate that pressure increases and velocities are sensitive to the overall amount of expansion rather than the behavior of the water-steam transition. Both the overall thermal expansion coefficient α and the temperature difference Δ T are greater in shallow (<1 km) environments than in deep (∼5 km) ones. Thus, for rocks with similar transport properties, pressure increases due to heating are greatest in shallow environments. Although solutions can be applied to rocks with a wide variety of properties, pressure increases are calculated for compliant quartz-rich sedimentary rocks with a porosities between 1 and 20% and permeabilities between 1 darcy and 1 μdarcy, subject to temperature increases of 500 and 1000 K at depths ranging from 0.1 to 5 km in a region of hydrostatic pressures and normal geothermal gradient. Such rocks, with porosities greater than 5%, permeabilities less than a 0.1 mdarcy, and drained hydrostatic compressibilities of 10 −4 /MPa, undergo pressure increases greater than 10 MPa (100 bars)for conditions typical of water table depths of 2.5 km and heating to 500 K above ambient. Similar rocks, but with permeabilities less than 1 mdarcy, undergo pressure increases of 10 MPa for conditions typical of 1 km water table depth. Rocks commonly considered to be good aquifers undergo pressure increases of less than 1 MPa, primarily because of their high permeability. Although these estimates neglect the effects of fracturing and brecciation that may accompany such pressure increases, calculations indicate that pressure increases due to heating of cool groundwater can lead to failure of host rocks by a phreatic mechanism.

Journal of Geophysical Research Solid Earth↗
USGS2012

Extension in the Rio Grande rift

A positive gravity anomaly along the axis of the Rio Grande rift reflects a volume of anomalous mass added at the base of the crust and intruded into the crust. Part of this volume can be associated with vertical uplift of the crust. The remainder of this anomalous volume, plus the volume of surficial graben fill, can be associated with horizontal crustal extension. The volume of crustal uplift in the Rio Grande rift is unknown, but this term can be eliminated by means of an independent equation provided by assumption of generalized isostatic equilibrium. The volume and mass equations combined provide a solution for extension of the crust in terms of the following parameters: total anomalous mass deficiency in the mantle lithosphere, total anomalous mass excess in the crust and its density contrast, total anomalous mass deficiency of surficial graben fill and its density contrast, and the volume of material eroded from the uplift. Using standard density estimates and masses determined by equivalent-source modeling of gravity profiles, I obtained 1-km extension at 37°N (Colorado-New Mexico border), 13-km extension at 35°N (Albuquerque, New Mexico), and 24-km extension at 33°N in southern New Mexico. These estimates are of necessity very poorly constrained and may be as much as ±60% in error. Taking results at face value, the best fitting Euler pole occurs at about 41°N, in north central Colorado. The calculated opening angle is 1.32°. For an approximately 30-m.y. duration of the extensional system (late Oligocene to the present), calculated average angular velocity is 7.8×10 −4 rad/m.y.; calculated average spreading half-rate at 33°N (for example) is a stately 0.04 cm/yr. The ratio of extension to uplift increases southward.

Journal of Geophysical Research Solid Earth↗
USGS2012

Permafrost, heat flow, and the geothermal regime at Prudhoe Bay, Alaska

Temperature measurements through permafrost in the oil field at Prudhoe Bay, Alaska, combined with laboratory measurements of the thermal conductivity of drill cuttings permit an evaluation of in situ thermal properties and an understanding of the general factors that control the geothermal regime. A sharp contrast in temperature gradient at ∼600 m represents a contrast in thermal conductivity caused by the downward change from interstitial ice to interstitial water at the base of permafrost under near steady state conditions. Interpretation of the gradient contrast in terms of a simple model for the conductivity of an aggregate yields the mean ice content (∼39%), and thermal conductivities for the frozen and thawed sections (8.1 and 4.7 mcal/cm s °C, respectively). These results yield a heat flow of ∼1.3 HFU, which is similar to other values on the Alaskan Arctic Coast; the anomalously deep permafrost is a result of the anomalously high conductivity of the siliceous ice-rich sediments. Curvature in the upper 160 m of the temperature profiles represents a warming of ∼1.8°C of the mean surface temperature and a net accumulation of 5–6 kcal/cm 2 by the solid earth surface during the last 100 years or so. Rising sea level and thawing of ice-rich sea cliffs probably caused the shoreline to retreat tens of kilometers in the last 20,000 years, inundating a portion of the continental shelf that is presently the target of intensive oil exploration. A simple conduction model suggests that this recently inundated region is underlain by near-melting ice-rich permafrost to depths of 300–500 m; its presence is important to seismic interpretations in oil exploration and to engineering considerations in oil production. With confirmation of the permafrost configuration by offshore drilling, heat conduction models can yield reliable new information on the chronology of arctic shorelines.

Journal of Geophysical Research Solid Earth↗
USGS2012

Ice sculpture in the Martian outflow channels

Many landforms in Martian outflow channels have characteristics that suggest sculpture by glaciers, ice streams, or ice sheets. Viking Orbiter and terrestrial satellite images were examined at similar resolution to compare features of the Martian outflow channels to features produced by the movement of ice on earth. Many resemblances were found. They include the anastomoses, sinuosities, and U‐shaped cross profiles of valleys; hanging valleys; and linear scour marks on valley walls, grooves and ridges on valley floors, and streamlining of bedrock highs. The question of whether ice could have moved in the Martian environment is investigated. Since gradients on Martian channel floors are extremely low or even reversed, a certain thickness of ice is required to initiate glacial motion. If ice filled the Martian channels to the level of interior plateaus, flow could have occurred, even under present climatic conditions, although warmer climates in the past would have been more favorable. Additionally, if the channel fluids were brines, ice might have been wet based and would have moved readily by slipping over the ground. It is envisioned that springs or small catastrophic outbursts discharged fluids from structural outlets or chaotic terrains. The fluids built icings that may have grown into substantial masses and eventually flowed like glaciers down preexisting valleys. Alternatively, the fluids may have formed rivers or floods that formed ice jams and consolidated into icy masses in places where obstacles blocked their flow. These masses of ice or slush were probably episodic and temporary, but where they occurred within the channels, they may have produced landforms that are characteristic of glacial sculpture on earth.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

A teleseismic analysis of the New Brunswick earthquake of January 9, 1982

The analysis of the New Brunswick earthquake of January 9, 1982, has important implications for the evaluation of seismic hazards in eastern North America. Although moderate in size ( m b 5.7), it was well-recorded teleseismically. Source characteristics of this earthquake have been determined from analysis of data that were digitally recorded by the Global Digital Seismograph Network. From broadband displacement and velocity records of P waves, we have obtained a dynamic description of the rupture process as well as conventional static properties of the source. The depth of the hypocenter is estimated to be 9 km from depth phases. The focal mechanism determined from the broadband data corresponds to predominantly thrust faulting. From the variation in the waveforms the direction of slip is inferred to be updip on a west dipping NNE striking fault plane. The steep dip of the inferred fault plane suggests that the earthquake occurred on a preexisting fault that was at one time a normal fault. From an inversion of bodywave pulse durations, the estimated rupture length is 5.5 km. Average properties of the rupture process were examined by a moment tensor analysis of long-period P and SH body waves. The long-period moment of this earthquake was 5.3 × 10 24 dyne cm. The static and dynamic stress drops are 41 and 65 bars, respectively, similar to those of many earthquakes with similar moment in regions that are more seismically active. The joint epicenter determination algorithm was used to locate, relative to the mainshock, the three teleseismically recorded aftershocks that occurred through March 31, 1982. The relocated hypocenters of the aftershocks are significantly different from each other and from that of the mainshock; they provide additional support for the source dimensions inferred from the waveform analysis.

Journal of Geophysical Research Solid Earth↗
USGS2012

High paleointensities of the geomagnetic field from thermomagnetic studies on rift valley pillow basalts from the Mid- Atlantic Ridge

Nineteen pillow basalts dredged within the rift valley of the Mid-Atlantic Ridge at36.8°N were studied by the Thellier stepwise heating method in order to determine the paleointensity of the geomagnetic field when they erupted on to the sea floor. Previously reported fission track ages are 2,000 to 6,000 years for the youngest rocks (mainly olivine basalts) and 10,000 to 100,000 years for the others (mainly plagioclase basalts and pyroxene basalts). All but three pillow basalts meet the conditions commonly considered as indicative of quite reliable paleointensity estimates: stability of the direction of NRM during its thermal demagnetization, constant ratio of NRM/TRM (natural remanent magnetization to thermoremanent magnetization) over 50% or more of the original NRM intensity (80 to 94% for 11 specimens), and reproducibility of low-temperature partial TRM (PTRM). However, strong field thermomagnetic measurements indicate that 11 of these 16 samples display a significant increase in Curie temperature (15 to 80°C) during the paleointensity experiments below 250°C, notwithstanding the linearity of the NRM-TRM plot in this temperature interval. This alteration, probably due to low-temperature oxidation of the specimens, seems typical of young pillow basalts and may result in paleointensity estimates which are too high. This result shows that excellence of remanence tests (NRM-TRM linearity and PTRM stability) does not ensure the absence of chemical changes during the Thellier experiments and therefore the validity of the paleointensity obtained. The assumption that reliable paleointensities are obtained when the Curie point increase is 10°C or less led to the selection of five of our specimens, all from the youngest group. Their mean paleointensity is 64.2±20.5 μT (standard deviation) and the corresponding virtual axial dipole moment (VADM) is 11.5±3.7× 10 22 A m 2 . Given the variations of the VADM of the Earth's magnetic field over the last 6,000 years, as established from archeomagnetic studies, our paleointensity results suggest that the latest eruptions on the inner floor of the Rift Valley at 36.8°N occurred 1,500±1,000 years ago.

Journal of Geophysical Research Solid Earth↗
USGS2012

In situ studies of velocity in fractured crystalline rocks

A study of the effects of macroscopic fractures on P and S wave velocities has been conducted in four wells drilled in granitic rock to depths between 0.6 and 1.2 km. The effect of macroscopic fractures is to decrease both V p and V s and increase V p / V s . In wells with a relatively low density of macroscopic fractures, the in situ velocity is similar to that of saturated core samples under confining pressure in the laboratory, and there is a clear correlation between zones with macroscopic fractures and anomalously low velocities. In wells with numerous macroscopic fractures, the in situ velocity is lower than that of intact samples under pressure, and there is a correlation between the rate at which in situ velocity increases with depth and the rate at which the velocity of laboratory samples increases with pressure. Differences in in situ P wave velocity between wells cannot be explained solely by differences in the degree of macroscopic fracturing, thus emphasizing the importance of composition and microcracks on velocity. In one highly fractured well the in situ P wave velocity is essentially the same for frequencies ranging from 10 Hz to 20 kHz; this suggests that the macrofractures affect velocity similarly over a broad frequency range. Chemical alteration of rock adjacent to macroscopic fractures appears to play an important role in reducing in situ velocities. Synthetic reflection seismograms generated from the velocity logs suggest that fracture zones are one possible source of deep-crustal reflectors observed on seismic reflection profiles.

Journal of Geophysical Research Solid Earth↗
USGS2012

Investigation of internal friction in fused quartz, steel, Plexiglass, and Westerly granite from 0.01 to 1.00 Hertz at 10-8 to 10-7 strain amplitude

A detailed evaluation on the method of internal friction measurement by the stress-strain hysteresis loop method from 0.01 to 1 Hz at 10 −8 to 10 −7 strain amplitude and 23.9°C is presented. Significant systematic errors in relative phase measurement can result from convex end surfaces of the sample and stress sensor and from end surface irregularities such as nicks and asperities. Preparation of concave end surfaces polished to optical smoothness having a radius of curvature >3.6×10 4 cm reduces the systematic error in relative phase measurements to <(5.5±2.2)×10 −4 radians. The values of Q E −1 (internal friction under uniaxial compression) determined from the relative phase measurements are | Q E −1 – Q s −1 |< 2.8×10 −3 for the tool steel sample and | Q E −1 – Q s −1 |< 2.2×10 −3 for the Westerly granite sample, where Q s −1 is the internal friction of the fused quartz stress sensor under uniaxial compression. These values are consistent with those inferred from the relative modulus dispersion data also presented in this paper. The polymethyl methacrylate (PMM, trade name Plexiglass) sample shows high values of internal friction ( Q E −1 ≅5×10 −2 ) with strong frequency dependence and with a maximum in Q E −1 at ≅0.4 Hz.

Journal of Geophysical Research Solid Earth↗
USGS2012

Wandering terranes in southern Alaska: The Aleutia Microplate and implications for the Bering Sea

Paleomagnetic and geological data suggest that much of southern Alaska is a collage of tectonostratigraphic terranes which originated in Mesozoic time at paleolatitudes far south of their present position. The time of ‘docking’ of the terranes against cratonic Alaska is critical to defining their amalgamated size and extent during their northward motion as well as their role in the evolution of the Bering Sea. One of the largest of the tectonostratigraphic terranes, the Peninsular terrane of south central and southwestern Alaska, extends offshore along the outer Bering Sea continental margin (Beringia). Paleomagnetic data suggest that this terrane has moved northward through all of Cenozoic time, but geologic data imply that the terrane had accreted to Alaska by the end of the Mesozoic. In early Cenozoic time the eastern part of the Aleutian arc appears to have been superimposed on the Peninsular terrane, and postulated northward Cenozoic motion of the terrane would therefore have required northward motion of the arc. Two accretion models, based on docking times for terranes in Alaska, are proposed, and they illustrate that large areas of the abyssal Bering Sea, the Alaska Peninsula, the Aleutian arc, and the Beringian continental margin may be part of a superterrane or microplate called Aleutia (microplate as defined by Beck et al. (1980), i.e., a microplate is a displaced segment of lithosphere that has crustal roots, whereas a superterrane is an amalgamation of terranes which may or may not be rootless). Model A implies that the Aleutian arc developed in situ on the southern edge of Aleutia after the microplate had docked. In model B, the final docking time of the Peninsular terrane is late Cenozoic, which implies that the Aleutia microplate encompasses a mammoth area that includes parts of southern Alaska, the Alaska Peninsula, the southern Beringian margin, the abyssal Bering Sea (Kula plate), and the Aleutian arc. If model A is correct, the docking time of the Peninsular terrane is late Mesozoic or earliest Tertiary. The Aleutia microplate in this model is made up solely of the abyssal Bering Sea (Kula plate), which presumably docked at the same time or slightly after the Peninsular terrane accreted against Alaska. If model B is correct, that is, if the Aleutia collided with nuclear Alaska during the Cenozoic, then a late Cenozoic suture zone, the vestige of a large open sea that must have closed between Aleutia and Alaska, must exist in south central and southwest Alaska. Either evidence for Cenozoic closure and suturing has been obliterated in Alaska or the inferences of Cenozoic terrane motion derived from paleomagnetic data are suspect.

Alaska↗
USGS2012

A dislocation model of strain accumulation and release at a subduction zone

Strain accumulation and release at a subduction zone are attributed to stick slip on the main thrust zone and steady aseismic slip on the remainder of the plate interface. This process can be described as a superposition of steady state subduction and a repetitive cycle of slip on the main thrust zone, consisting of steady normal slip at the plate convergence rate plus occasional thrust events that recover the accumulated normal slip. Because steady state subduction does not contribute to the deformation at the free surface, deformation observed there is completely equivalent to that produced by the slip cycle alone. The response to that slip is simply the response of a particular earth model to embedded dislocations. For a purely elastic earth model, the deformation cycle consists of a coseismic offset followed by a linear‐in‐time recovery to the initial value during the interval between earthquakes. For an elastic‐viscoelastic earth model (elastic lithosphere over a viscoelastic asthenosphere), the postearthquake recovery is not linear in time. Records of local uplift as a function of time indicate that the long‐term postseismic recovery is approximately linear, suggesting that elastic earth models are adequate to describe the deformation cycle. However, the deformation predicted for a simple elastic half‐space earth model does not reproduce the deformation observed along the subduction zones in Japan at all well if stick slip is restricted to the main thrust zone. As recognized earlier by Shimazaki, Seno, and Kato, the uplift profiles could be explained if stick slip were postulated to extend along the plate interface beyond the main thrust zone to a depth of perhaps 100 km, but independent evidence suggests that stick slip at such depths is unlikely.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Crustal and upper mantle structure of the northern and central Sierra Nevada

Teleseismic data were recorded within the Sierra Nevada to look for lateral variations in the upper mantle. The data were collected at both temporary and permanent stations, and P wave residuals were computed. After correcting the P residual data for crustal and topographic effects, there is still a variation of as much as 0.5-0.6 s from the north end of the Sierra Nevada to Mono Lake, located east of the central part of the range. In addition, there are significant variations in travel time patterns, depending on the azimuth of wave arrivals. Two simple modeling approaches have been used to infer the upper mantle velocity structure from the observed variations in travel time. In the first, it is assumed that the velocity variations are distributed throughout a depth range of about 60–160 km. One model that fits the data indicates a north-to-south trend in upper mantle P velocities from about 7.9 km/s at the north end of the Sierra to 7.6 km/s near Mono Lake. Superimposed is a west-to-east trend from 7.85 km/s under the Sierra Nevada crest to 7.7 km/s east of Lake Tahoe. In the second modeling approach, a fixed velocity contrast is assumed between the lithosphere and asthenosphere with variations in travel time resulting from variations in lithospheric thickness. One model, consistent with the observations, indicates a thinning of the lithosphere from 110 km at the north end of the Sierra to 60 km near Mono Lake.

Journal of Geophysical Research Solid Earth↗
USGS2012

Nonlinear strain buildup and the earthquake cycle on the San Andreas fault

Two contrasting models of the earthquake deformation cycle on strike slip faults predict significant temporal declines in shear strain rate near the fault, accompanied by a progressive broadening of the zone of deformation adjacent to it. In the thin lithosphere model, transient deformation results from flow in the asthenosphere due to stress relaxation following faulting through most or all of the lithosphere. For an earth model with a thick elastic lithosphere (plate thickness » depth of seismic slip), transient motions are due to postearthquake aseisrnic slip below the coseismic fault plane. Data from the San Andreas fault indicate a long-term temporal decrease in strain rate that persists for at least 30 years and may extend through the entire earthquake cycle. Observations support a cycle-long rate decrease and a temporal spreading of the deformation profile only if movement cycles on the northern and southern locked sections of the fault are basically similar. If so, the usually lower strain rates and broader deformation zone currently observed on the southern San Andreas represent a later evolutionary stage of the northern locked section, where a great earthquake is a more recent occurrence. Although the data allow some extreme models to be discarded, no sufficiently strong constraints exist to decide between the thin and thick lithosphere models. Regardless of the appropriate model the geodetic observations themselves indicate that strain buildup is sufficiently nonlinear to cause significant departures from recurrence estimates based on linear strain accumulation and the time-predictable model.

Journal of Geophysical Research Solid Earth↗
USGS2012

The use of MAGSAT data to determine secular variation

A combined spatial and secular variation model of the geomagnetic field, labeled M061581, is derived from a selection of Magsat data. Secular variation (SV) data computed from linear fits to midnight hourly values from 19 magnetic observatories were also included in the analysis but were seen to have little effect on the model. The SV patterns from this new model are compared with those from the 1980 IGRF and with those for 1970 computed by the DGRF and with the 1960 patterns computed using the GSFC(12/66) model. Most of the features of the M061581 are identical in location and level with those of the 1980 IGRF. Together they confirm that the reversals in sign of field change seen over Asia and North America between 1965 and 1975 are reverting to the pre-1965 states. The M061581 model gives −32 nT/yr for the dipole decay rate, larger than the 70% increase already reported since 1965. This abnormally high value is interpreted as being a defect of the model because it appears to result from a much larger (−100 nT/yr) drop in field over the polar regions not indicated by the 1980 IGRF. This north polar decrease is shown to be of external origin as the result of a combination of the seasonal effect of the north polar ionospheric (counterclockwise) afternoon Sq p cell increasing in intensity from the beginning (November 1979) to the end (June 1980) of the Magsat data collection period, coupled with an enhancement of its effect as the orbit lowered from the 350– 550 km initial altitudes to near 200 km just prior to burnup. This experiment indicates that secular variation can be obtained from satellite data for intervals of less than a full year if corrections can be made for seasonal effects and that ‘annual snapshots’ of the field by a satellite would allow easy and accurate models of secular change without the use of any surface data.

Journal of Geophysical Research Solid Earth↗
USGS2012

The 1979 Homestead Valley earthquake sequence, California: Control of aftershocks and postseismic deformation

The coseismic slip and geometry of the March 15, 1979, Homestead Valley, California, earthquake sequence are well constrained by precise horizontal and vertical geodetic observations and by data from a dense local seismic network. These observations indicate 0.52±0.10 m of right-lateral slip and 0.17±0.04 m of reverse slip on a buried vertical 6-km-long and 5-km-deep fault and yield a mean static stress drop of 7.2±1.3 MPa. The largest shock had M S = 5.6. Observations of the ground rupture revealed up to 0.1 m of right-lateral slip on two mapped faults that are subparallel to the modeled seismic slip plane. In the 1.9 years since the earthquakes, geodetic network displacements indicate that an additional 60±10 mm of postseismic creep took place. The rate of postseismic shear strain (0.53±0.13 μrad/yr) measured within a 30×30-km network centered on the principal events was anomalously high compared to its preearthquake value and the postseismic rate in the adjacent network. This transient cannot be explained by postseismic slip on the seismic fault but rather indicates that broadscale release of strain followed the earthquake sequence. We have calculated the postearthquake stress field caused by the modeled coseismic slip. We assume that failure is promoted when the sum of the shear stress plus 0.75 times the faultopening stress increases. Most aftershocks concentrate at points where the stresses are enhanced by 0.3 MPa (3 bars) or more; aftershocks are nearly absent where postearthquake stresses decrease by 0.3–0.5 MPa. Isolated off-fault clusters of aftershocks that locate at one fault length from the rupture plane are explainable by this hypothesis. We find that ground rupture and postseismic creep take place where near-surface stresses are calculated to increase within the preexisting fault zones. Two patches that extend 4 km from both ends of the seismic fault exhibited neither aftershocks nor measurable postseismic creep. The sensitivity of aftershocks and ground rupture to changes in stress that are less than 5% of the earthquake stress drop demonstrates that the region around the earthquakes was within a few percent of its failure threshold before the main shocks. The preearthquake stress field and the stress required for failure must also have been nearly uniform.

Journal of Geophysical Research Solid Earth↗