Reply to "Comment on 'A model of earthquake triggering probabilities and application to dynamic deformations constrained by ground motion observations' by Ross Stein"
Explore the source record for details and available documents.
SEARCH · Search USGS
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.
Explore the source record for details and available documents.
We use three-dimensional elastic models to help guide the kinematic interpretation of crustal deformation associated with strike-slip faults. Deformation of the brittle upper crust in the vicinity of strike-slip fault systems is modeled with the assumption that upper crustal deformation is driven by the relative plate motion in the upper mantle. The driving motion is represented by displacement that is specified on the bottom of a 15-km-thick elastic upper crust everywhere except in a zone of finite width in the vicinity of the faults, which we term the “shear zone.” Stress-free basal boundary conditions are specified within the shear zone. The basal driving displacement is either pure strike slip or strike slip with a small oblique component, and the geometry of the fault system includes a single fault, several parallel faults, and overlapping en echelon faults. We examine the variations in deformation due to changes in the width of the shear zone and due to changes in the shear strength of the faults. In models with weak faults the width of the shear zone has a considerable effect on the surficial extent and amplitude of the vertical and horizontal deformation and on the amount of rotation around horizontal and vertical axes. Strong fault models have more localized deformation at the tip of the faults, and the deformation is partly distributed outside the fault zone. The dimensions of large basins along strike-slip faults, such as the Rukwa and Dead Sea basins, and the absence of uplift around pull-apart basins fit models with weak faults better than models with strong faults. Our models also suggest that the length-to-width ratio of pull-apart basins depends on the width of the shear zone and the shear strength of the faults and is not constant as previously suggested. We show that pure strike-slip motion can produce tectonic features, such as elongate half grabens along a single fault, rotated blocks at the ends of parallel faults, or extension perpendicular to overlapping en echelon faults, which can be misinterpreted to indicate a regional component of extension. Zones of subsidence or uplift can become wider than expected for transform plate boundaries when a minor component of oblique motion is added to a system of parallel strike-slip faults.
The emission rate of carbon dioxide escaping from the summit of Kīlauea Volcano, Hawaiʻi, proved highly variable, averaging 4900 ± 2000 metric tons per day (t/d) in June–July 2003 during a period of summit inflation. These results were obtained by combining over 90 measurements of COSPEC-derived SO 2 emission rates with synchronous CO 2 /SO 2 ratios of the volcanic gas plume along the summit COSPEC traverse. The results are lower than the CO 2 emission rate of 8500 ± 300 t/d measured by the same method in 1995–1999 during a period of long-term summit deflation [Gerlach, T.M., McGee, K.A., Elias, T., Sutton, A.J. and Doukas, M.P., 2002. Carbon dioxide emission rate of Kīlauea Volcano: Implications for primary magma and the summit reservoir. Journal of Geophysical Research-Solid Earth, 107(B9): art. no.-2189.]. Analysis of the data indicates that the emission rates of the present study likely reflect changes in the magma supply rate and residence time in the summit reservoir. It is also likely that emission rates during the inflation period were heavily influenced by SO 2 pulses emitted adjacent to the COSPEC traverse, which biased CO 2 /SO 2 ratios towards low values that may be unrepresentative of the global summit gas plume. We conclude that the SO 2 pulses are consequences of summit re-inflation under way since 2003 and that CO 2 emission rates remain comparable to, but more variable than, those measured prior to re-inflation.
Geologists are in a quandary over the correct interpretation of paleomagnetic data for Jurassic rocks of the North American plate. Conflicting reference paleopoles and alternate methods of constructing apparent polar wander (APW) paths have led to a controversy regarding the configuration of Jurassic APW for cratonic North America. These differences have been recently disputed at meetings (see Eos , Spring Meeting Supplement, April 7, 1992, p. 94) and in an exchange of letters in the Journal of Geophysical Research-Solid Earth [ Butler et al ., 1992; Van Fossen and Kent , 1992a]. At stake is important information concerning the Jurassic paleogeography of North America and the whole of Pangea, as well as the nature and driving mechanism of plate movements. In addition, cratonic paleopoles provide a reference frame for measuring relative displacements of tectonostratigraphic terranes.
A technique to locate automatically boundaries between crustal blocks of disparate densities was applied to upward continued isostatic residual gravity data. The boundary analysis delineates a narrow gravitational trough that extends the length of the Pliocene and Quaternary volcanic arc from Mount Baker in northern Washington to Lassen Peak in California. Gravitational highs interrupt the trough at two localities: a northwest trending high in southern Washington and a northeast trending high between Mount Shasta and Lassen Peak. The latter anomaly is one of a set of northeast trending anomalies that, within the Quaternary arc, appear related to volcanic segmentation proposed previously on the basis of spatial and compositional distributions of volcanoes. These northeast trending anomalies extend hundreds of kilometers northeast of the arc, are caused by sources in the upper crust, and in some cases are related to exposed pre-Tertiary rocks. Segmentation models invoke geometric characteristics of the subducting plate as the primary factor controlling location and chemistry of volcanism, and these northeast trending gravity sources also may be a product of disturbance of the upper crust by the subduction process. More likely, the gravity sources may reflect upper crustal structures older than the High Cascades, possibly relicts from earlier accretionary events or more recent crustal deformation, that have actively influenced the spatial location of more recent volcanism. Much of the Pliocene and Quaternary volcanism of the Cascade arc has concentrated on or near contacts between crustal blocks of disparate density. These contacts may promote the ascension of magma to the Earth's surface.
Rapid changes in ground tilt and GPS positions on Kīlauea Volcano, Hawai'i, are interpreted as resulting from a shallow, two‐segment dike intrusion into the east rift zone that began at 1217 UTC (0217 HST) on 17 June 2007 and lasted almost 3 days. As a result of the intrusion, a very small volume of basalt (about 1500 m 3 ) erupted on 19 June. Northward tilt at a coastal tiltmeter, subsidence of south flank GPS sites, southeastward displacements at southwestern flank GPS sites, and a swarm of flank earthquakes suggest that a slow slip event occurred on the décollement beneath Kīlauea's south flank concurrent with the rift intrusion. We use 4 min GPS positions that include estimates of time‐dependent tropospheric gradients and ground tilt data to study the spatial and temporal relationships between the two inferred shallow, steeply dipping dike segments extending from the surface to about 2 km depth and décollement slip at 8 km depth. We invert for the temporal evolution of distributed dike opening and décollement slip in independent inversions at each time step using a nonnegative least squares algorithm. On the basis of these inversions, the intrusion occurred in two stages that correspond spatially and temporally with concentrated rift zone seismicity. The dike opening began on the western of the two segments before jumping to the eastern segment, where the majority of opening accumulated. Dike opening preceded the start of décollement slip at an 84% confidence level; the latter is indicated by the onset of northward tilt of a coastal tiltmeter. Displacements at southwest flank GPS sites began about 18 h later and are interpreted as resulting from slow slip on the southwestern flank. Additional constraints on the evolution of the intrusion and décollement slip come from inversion of an Envisat interferogram that spans the intrusion until 0822 UTC on 18 June 2007, combined with GPS and tilt data. This inversion shows that up to 0822 UTC on 18 June, décollement slip is only required in a limited region offshore of Ka'ena Point. A similar inversion of the complete event, which includes GPS and tilt data up to 21 June and a second Envisat interferogram spanning the complete intrusion until 21 June, shows décollement slip spread westward across the south flank. This may suggest westward migration of the décollement slip as the event progressed.
A major obstacle in lunar U-Pb chronology is the elusiveness of the primary Pb isotopic composition and U/Pb ratio and therefore the Pb evolution for the early history of the moon. In an attempt to seek the primitive lunar Pb isotopic composition, green glass from lunar clod 15426,49 was studied for U-Th-Pb systematics because it is extremely Mg-rich and known to be the most primitive among sampled lunar volcanic rocks. Green glass and yellow impact glass were handpicked, washed with water, then washed in an acid mixture (0.3 N HBr + 0.2 N HF), and the washes and residues were analyzed. In contrast to the Pb from Apollo 17 orange glass, Pb from the surface of the green glass spherules is totally different from that of the spherule interior, but is similar to lunar surface lead. Data for the interior Pb yield a 207 Pb/ 206 Pb age of 3.41 Ga. Because of the low Pb concentration and high U/Pb ratio observed for the interior, the initial Pb was poorly defined. Nevertheless, the data indicate that lunar Pb evolved in an environment with 238 U/ 204 Pb = 19–55, which is considerably lower than those for mare basalts (around 300) but higher than values for the earth (6–8).
The water level in an open well can change in response to deformation of the surrounding material, either because of applied strains (tidal or tectonic) or surface loading by atmospheric pressure changes. Under conditions of no vertical fluid flow and negligible well bore storage (static-confined conditions), the sensitivities to these effects depend on the elastic properties and porosity which characterize the surrounding medium. For a poroelastic medium, high sensitivity to applied areal strains occurs for low porosity, while high sensitivity to atmospheric loading occurs for high porosity; both increase with decreasing compressibility of the solid matrix. These material properties also influence vertical fluid flow induced by areally extensive deformation and can be used to define two types of hydraulic diffusivity which govern pressure diffusion, one for applied strain and one for surface loading. The hydraulic diffusivity which governs pressure diffusion in response to surface loading is slightly smaller than that which governs fluid flow in response to applied strain. Given the static-confined response of a water well to atmospheric loading and Earth tides, the in situ drained matrix compressibility and porosity (and hence the one-dimensional specific storage) can be estimated. Analysis of the static-confined response of five wells to atmospheric loading and Earth tides gives generally reasonable estimates for material properties.
A combined aeromagnetic and radio echo ice-sounding survey (4200 km of traverse) made in 1978 in Antarctica over the Dufek layered mafic intrusion of Jurassic age suggests a minimum area of about 50,000 km 2 , making it comparable in size with the Bushveld Complex of Africa. Comparisons of the magnetic and subglacial topographic profiles illustrate the usefulness of this combination of methods in studying bedrock geology beneath ice-covered areas. Rocks are exposed in only 3% of the inferred area of the intrusion. Magnetic anomalies measured a few hundred meters above outcrops of the intrusion range in peak-to-trough amplitude from ∼50 nT over the lowermost exposed portion of the section in the Dufek Massif to ∼3600 nT over the uppermost part of the section in the Forrestal Range. Theoretical magnetic anomalies, computed from models based on the subice topography fitted to the highest-amplitude observed magnetic anomalies, required normal and reversed magnetizations ranging from 10 −4 to 10 −2 emu/cm 3 having directions and magnetizations consistent with measurements previously made on oriented samples. This result is interpreted as indicating that the Dufek intrusion cooled through the Curie isotherm during one or more reversals of the earth's magnetic field.
K-Ar ages and paleomagnetic data for basalt samples from a new core hole (site E) at the Idaho National Engineering Laboratory (INEL) indicate that the age of the reversed polarity event recorded in Snake River Plain lavas is older than 465±50 ka (1000 years before present) reported previously by Champion et al. (1981). Nine basalt flows, eight with normal polarity and one with reversed polarity, were recognized in the site E core hole. The flows above and below the reversed flow have ages of 491±80 ka and 580±93 ka, respectively. The inclination of the paleomagnetic field direction of the reversed flow at site E agrees with the inclination of reversed flows elsewhere at INEL which have an age of 565±14 ka. These reversed flows were previously thought to be correlative with the Emperor event. We suggest that this polarity event is an older event which we name the Big Lost Reversed Polarity Subchronozone and Subchron. A review of data documenting short reversal records from volcanic and sedimentary rocks shows that there is evidence for eight polarity subchrons in the Brunhes and two besides the Jaramillo in the late Matuyama. These 10 short subchrons begin to indicate the many short events that Cox (1968) hypothesized must exist if polarity interval lengths have a Poisson distribution. These events are true subchrons, not excursions, and may or may not have low associated paleointensities, although low field strengths might explain why the reversal process aborts. The mean sustained polarity interval length since late Matuyama Chron time is 90,000 years. The similarity of this number with the 10 5 -year period of the Earth's orbital eccentricity suggests anew that linkage between geomagnetic, paleoclimatic, and possible underlying Earth orbital parameters should be evaluated.
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.
New seismic and geodetic data from Costa Rica provide insight into seismogenic zone processes in Central America, where the Cocos and Caribbean plates converge. Seismic data are from combined land and ocean bottom deployments in the Nicoya peninsula in northern Costa Rica and near the Osa peninsula in southern Costa Rica. In Nicoya, inversion of GPS data suggests two locked patches centered at 14 ± 2 and 39 ± 6 km depth. Interplate microseismicity is concentrated in the more freely slipping intermediate zone, suggesting that small interseismic earthquakes may not accurately outline the updip limit of the seismogenic zone, the rupture zone for future large earthquakes, at least over the short (∼1 year) observation period. We also estimate northwest motion of a coastal “sliver block” at 8 ± 3 mm/yr, probably related to oblique convergence. In the Osa region to the south, convergence is orthogonal to the trench. Cocos-Caribbean relative motion is partitioned here, with ∼8 cm/yr on the Cocos-Panama block boundary (including a component of permanent shortening across the Fila Costeña fold and thrust belt) and ∼1 cm/yr on the Panama block–Caribbean boundary. The GPS data suggest that the Cocos plate–Panama block boundary is completely locked from ∼10–50 km depth. This large locked zone, as well as associated forearc and back-arc deformation, may be related to subduction of the shallow Cocos Ridge and/or younger lithosphere compared to Nicoya, with consequent higher coupling and compressive stress in the direction of plate convergence.
No abstract available.
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.
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.
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.
As part of the Parkfield, California, earthquake prediction experiment, water level is monitored in a well 460 m from the main trace of the San Andreas fault on Middle Mountain, in the preparation zone of the anticipated Parkfield earthquake. The well configuration allows water level to be monitored in two fluid reservoirs at depths of 85 and 250 m below land surface. During 1987, water level changes were recorded during 12 of the 18 episodes of accelerated fault creep detected by a creep meter spanning the fault trace 750 m northwest of the well. The creep-related water level changes in the shallow reservoir have durations of less than 1 day, whereas in the deeper reservoir the changes persist for as long as 2 months. These data suggest that the transient nature of the water level changes in the shallow interval is due to vertical flow to the water table and is not evidence that creep events propagate past the well. Phase leads of earth tidal constituents in the water level data from the shallow interval relative to the same constituents in the local volume strain tide support the interpretation of significant flow to the water table at periods of 1 day or less. The form of the water level changes in the deep interval is affected by horizontal flow to the well bore. This effect can be removed from the water level records using a theoretical response curve constrained by the phases of earth tidal constituents in the deep interval relative to the local volume strain tide. For the events where the signal in the shallow interval has been large enough to measure, the sizes of the simultaneous water level changes in the two reservoirs are consistent with the same amounts of volume strain occurring at both depths.
Large landslides in the Martian equatorial troughs have been investigated with respect to morphology, geologic structure of the troughs, time of emplacement, similarity to terrestrial landslides, and origin and mechanism of transport. Viking Orbiter images provided a base for morphologic analyses, mapping, crater counts, and volumetric considerations. Results showed that the morphologic variations of the landslides can be attributed mainly to their degree of confinement on trough floors. Many prominent landslides appear to be of similar age and were emplaced after a major episode of faulting that created scarps of great relief. Most sliding apparently occurred after these scarps were dissected into spurs, gullies, and tributary canyons. Much of the slide activity coincided with a major late eruptive phase of the Tharsis volcanoes, as deduced from the crater density, over the total deposit area, of 570±130 craters larger than 1 km in diameter per 10 6 km 2 . The huge size of many landslides is due to their occurrence on fault scarps that may have attained several kilometers in height in the absence of vigorous fluvial erosion on Mars. The mechanical efficiency of the Martian landslides is high but in accord with predictions from large landslides on earth. The slides may have originated as gigantic mudflows with slump blocks at their heads; water that may have existed in aquifers in the wallrock below and behind a lid of ice may have been liberated through minor disturbances from Marsquakes. Liquefaction and collapse of extensive wall sections may then have ensued, and the debris rushed out onto the trough floors at speeds exceeding 100 km/h.