Comment on "The 1946 Hispaniola earthquake and the tectonics of the North America-Caribbean plate boundary zone, northeastern Hispaniola"
No abstract available.
SEARCH · USGS Science
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.
No abstract available.
Earthquake focal mechanisms of events occurring between 1972 and 1992 in the south flank of Kilauea volcano, Hawaii, are used to infer the state of stress and strain as a function of time and space. We have determined 870 fault plane solutions from P wave first motion polarities for events with magnitudes M L ≥ 2.5 and depth ranging between 6 and 12 km. Faulting is characterized by a mixture of decollement, reverse, and normal faults. Most large earthquakes with magnitude M < 7 slip on reverse faults striking NE at 40° and dipping SE between 60° and 70°. In Hawaii, the earthquakes with M > 7 rupture the decollement plane, since it is the only surface large enough to generate magnitude 7 or larger earthquakes. The percentage of reverse faulting events is high compared to the decollement and normal faulting mechanisms for the period 1972–1983. The percentage of decollement type focal mechanisms becomes dominant after 1983. This pattern of faulting activity suggests that pressure was building up within Kilauea's rift zone prior to the 1983 Puu'Oo eruption. Overall, a single stress orientation with the maximum compressive stress oriented SE perpendicular to the rift and dipping at 45° is compatible with the coeval existence of decollement, reverse, and normal faults. However, in a crustal volume east of longitude 155°10′W, we find a change of the orientation of σ 1 from nearly horizontal to plunging 45° SE occurring in 1979. This stress rotation suggests magma movements within the aseismic part of Kilauea's east rift zone. The strain and stress orientations are coaxial in the south flank except within the volume where the stress rotation is observed. We observe a change in the relationship between stress and strain directions caused either by the shifting of seismic activity from reverse faults to decollements, while stress stays constant, or by a rotation of stress, while strain remains constant. Assuming that the model of a noncohesive Coulomb wedge is appropriate for Kilauea's south flank, we find that high pore pressures are prevalent along the decollement and within the wedge for a coefficient of friction equal to 0.85.
A two-step nonlinear and linear inversion is carried out to map the lateral heterogeneity beneath North America using surface wave data. The lateral resolution for most areas of the model is of the order of several hundred kilometers. The most obvious feature in the tomographic images is the rapid transition between low velocities in the tectonically active region west of the Rocky Mountains and high velocities in the stable central and eastern shield of North America. The model also reveals smaller-scale heterogeneous velocity structures. A high-velocity anomaly is imaged beneath the state of Washington that could be explained as the subducting Juan de Fuca plate beneath the Cascades. A large low-velocity structure extends along the coast from the Mendocino to the Rivera triple junction and to the continental interior across the southwestern United States and northwestern Mexico. Its shape changes notably with depth. This anomaly largely coincides with the part of the margin where no lithosphere is consumed since the subduction has been replaced by a transform fault. Evidence for a discontinuous subduction of the Cocos plate along the Middle American Trench is found. In central Mexico a transition is visible from low velocities across the Trans-Mexican Volcanic Belt (TMVB) to high velocities beneath the Yucatan Peninsula. Two elongated low-velocity anomalies beneath the Yellowstone Plateau and the eastern Snake River Plain volcanic system and beneath central Mexico and the TMVB seem to be associated with magmatism and partial melting. Another low-velocity feature is seen at depths of approximately 200 km beneath Florida and the Atlantic Coastal Plain. The inversion technique used is based on a linear surface wave scattering theory, which gives tomographic images of the relative phase velocity perturbations in four period bands ranging from 40 to 150 s. In order to find a smooth reference model a nonlinear inversion based on ray theory is first performed. After correcting for the crustal thickness the phase velocity perturbations obtained from the subsequent linear waveform inversion for the different period bands are converted to a three-layer model of S velocity perturbations (layer 1, 25–100 km; layer 2, 100–200 km; layer 3, 200–300 km). We have applied this method on 275 high-quality Rayleigh waves recorded by a variety of instruments in North America (IRIS/USGS, IRIS/IDA, TERRAscope, RSTN). Sensitivity tests indicate that the lateral resolution is especially good in the densely sampled western continental United States, Mexico, and the Gulf of Mexico.
The Colorado Plateau is blanketed by Phanerozoic marine and nonmarine strata as young as Cretaceous that are now exposed at elevations of about 2 km. Crustal thickening due to magmatism and horizontal crustal shortening was far less than necessary to cause this uplift, which is commonly attributed to the consequences of mantle lithosphere thinning and heating. The Colorado Plateau and the midcontinent region around Iowa consist of Precambrian bedrock overlain by a similar amount of Paleozoic platformal strata, and thus both regions once had similar lithospheric buoyancy. Mesozoic sedimentation increased the crustal thickness and lithospheric buoyancy of the Colorado Plateau relative to the midcontinent region. Backstripping calculations yield elevation without these sediments and lead to a calculated elevation difference between the two areas of about 1200 m, which represents unexplained plateau uplift. Review of constraints on uplift timing finds little support for a late Cenozoic uplift age and allows early to middle Cenozoic uplift, which is consistent with uplift mechanisms related to low-angle subduction that ended in the middle Cenozoic. Finite element heat flow calculations of low-angle subduction and lithosphere attenuation, using a range of initial lithosphere thicknesses and degree of attenuation, indicate that required uplift can result from tectonic removal of about 120 km of mantle lithosphere from an initially 200-km-thick lithosphere. This allows for partial preservation of North American mantle lithosphere with its distinctive isotopic signature in some late Cenozoic volcanic rocks and is consistent with normal Pn velocities in the uppermost mantle beneath the plateau.
A spatially fixed or at least internally rigid hotspot reference frame has been assumed for determining relative plate motions on Earth. Recent 1:5,000,000 scale mapping of Venus, a planet without terrestrial-style plate tectonics and ocean cover, reveals a systematic age and dimensional progression of corona-like arachnoids occurring in an uncinate chain. The nonrandom associations between arachnoids indicate they likely formed from a deep-seated mantle plume in a manner similar to terrestrial hotspot features. However, absence of expected convergent "plate" margin deformation suggests that the arachnoids are the surface expression of a migratory mantle plume beneath a stationary surface. If mantle plumes are not stationary on Venus, what if any are the implications for Earth?
In 1954, four earthquakes of M > 6.0 occurred within a 30 km radius in a period of six months. The Rainbow Mountain-Fairview Peak-Dixie Valley earthquakes are among the largest to have been recorded geodetically in the Basin and Range province. The Fairview Peak earthquake (M = 7.2, December 12, 1954) followed two events in the Rainbow Mountains ( M = 6.2, July 6, and M = 6.5, August 24, 1954) by 6 months. Four minutes later the Dixie Valley fault ruptured ( M = 6.7, December 12, 1954). The changes in static stresses caused by the events are calculated using the Coulomb-Navier failure criterion and assuming uniform slip on rectangular dislocations embedded in an elastic half-space. Coulomb stress changes are resolved on optimally oriented faults and on each of the faults that ruptured in the chain of events. These calculations show that each earthquake in the Rainbow Mountain-Fairview Peak-Dixie Valley sequence was preceded by a static stress change that encouraged failure. The magnitude of the stress increases transferred from one earthquake to another ranged from 0.01 MPa (0.1 bar) to over 0.1 MPa (1 bar). Stresses were reduced by up to 0.1 MPa over most of the Rainbow Mountain-Fairview Peak area as a result of the earthquake sequence.
In 1954, four earthquakes greater than M s = 6.0 occurred within a 30-km radius and in a period of 6 months. Elevation and angle changes calculated from repeated leveling and triangulation surveys which span the coseismic period provide constraints on the fault geometries and coseismic slip of the faults which were activated. The quality of the coseismic geodetic data is assessed. Corrections are applied to the leveling data for subsidence due to groundwater withdrawal in the Fallon area, and a rod miscalibration error of 150 ± 30 ppm is isolated in leveling surveys made in 1967. The leveling and triangulation observations are then simultaneously inverted using the single value decomposition (SVD) inversion method to determine fault geometries and coseismic slip. Using SVD, it is possible to determine on which faults slip is resolvable given the data distribution. The faults are found to dip between 50° and 80° and extend to depths of 5 to 14 km. The geodetically derived slip values are generally equal to, or greater than, the maximum observed displacement along the surface scarps. Where slip is resolvable the geodetic data indicates the 1954 sequence contained a significant component of right-lateral slip. This is consistent with the N15°W trending shear zone which geodetic surveys have detected in western Nevada.
A hybrid global search algorithm is used to solve the nonlinear problem of calculating slip amplitude, rake, risetime, and rupture time on a finite fault. Thirty-five strong motion velocity records are inverted by this method over the frequency band from 0.1 to 1.0 Hz for the Northridge earthquake. Four regions of larger-amplitude slip are identified: one near the hypocenter at a depth of 17 km, a second west of the hypocenter at about the same depth, a third updip from the hypocenter at a depth of 10 km, and a fourth updip from the hypocenter and to the northwest. The results further show an initial fast rupture with a velocity of 2.8 to 3.0 km/s followed by a slow termination of the rupture with velocities of 2.0 to 2.5 km/s. The initial energetic rupture phase lasts for 3 s, extending out 10 km from the hypocenter. Slip near the hypocenter has a short risetime of 0.5 s, which increases to 1.5 s for the major slip areas removed from the hypocentral region. The energetic rupture phase is also shown to be the primary source of high-frequency radiation (1–15 Hz) by an inversion of acceleration envelopes. The same global search algorithm is used in the envelope inversion to calculate high-frequency radiation intensity on the fault and rupture time. The rupture timing from the low- and high-frequency inversions is similar, indicating that the high frequencies are produced primarily at the mainshock rupture front. Two major sources of high-frequency radiation are identified within the energetic rupture phase, one at the hypocenter and another deep source to the west of the hypocenter. The source at the hypocenter is associated with the initiation of rupture and the breaking of a high-stress-drop asperity and the second is associated with stopping of the rupture in a westerly direction.
A weighted equation based on the three-phase time-average and Wood equations is applied to derive a relationship between the compressional wave (P wave) velocity and the amount of hydrates filling the pore space. The proposed theory predicts accurate P wave velocities of marine sediments in the porosity range of 40-80% and provides a practical means of estimating the amount of in situ hydrate using seismic velocity. The shear (S) wave velocity is derived under the assumption that the P to S wave velocity ratio of the hydrated sediments is proportional to the weighted average of the P to S wave velocity ratios of the constituent components of the sediment. In the case that all constituent components are known, a weighted equation using multiphase time-average and Wood equations is possible. However, this study showed that a three-phase equation with modified matrix velocity, compensated for the clay content, is sufficient to accurately predict the compressional wave velocities for the marine sediments. This theory was applied to the laboratory measurements of the P and S wave velocities in permafrost samples to infer the amount of ice in the unconsolidated sediment. The results are comparable to the results obtained by repeatedly applying the two-phase wave scattering theory. The theory predicts that the Poisson's ratio of the hydrated sediments decreases as the hydrate concentration increases and the porosity decreases. In consequence, the amplitude versus offset (AVO) data for the bottom-simulating reflections may reveal positive, negative, or no AVO anomalies depending on the concentration of hydrates in the sediments.
Local earthquake P traveltime data is inverted to obtain a three-dimensional tomographic image of the region centered on the junction of the San Andreas and Calaveras faults. The resulting velocity model is then used to relocate more than 17,000 earthquakes and to produce a model of fault structure in the region. These faults serve as the basis for modeling the topography using elastic dislocation methods. The region is of interest because active faults join, it marks the transition zone from creeping to locked fault behavior on the San Andreas fault, it exhibits young topography, and it has a good spatial distribution of seismicity. The tomographic data set is extensive, consisting of 1445 events, 96 stations, and nearly 95,000 travel time readings. Tomographic images are resolvable to depths of 12 km and show significant velocity contrasts across the San Andreas and Calaveras faults, a low-velocity zone associated with the creeping section of the San Andreas fault, and shallow low-velocity sediments in the southern Santa Clara valley and northern Salinas valley. Relocated earthquakes only occur where v p > 5 km/s and indicate that portions of the San Andreas and Calaveras faults are non vertical, although we cannot completely exclude the possibility that all or part of this results from ray tracing problems. The new dips are more consistent with geological observations that dipping faults intersect the surface where surface traces have been mapped. The topographic modeling predicts extensive subsidence in regions characterized by shallow low-velocity material, presumably the result of recent sedimentation. Some details of the topography at the junction of the San Andreas and Calaveras faults are not consistent with the modeling results, suggesting that the current position of this “triple junction” has changed with time. The model also predicts those parts of the fault subject to contraction or extension perpendicular to the fault strike and hence the sense of any dip-slip component. In each locality the relative vertical motion across the fault is consistent with the fault dips found with the new hypocentral locations.
A model of a fault in the Earth is a sand-filled saw cut in a granite cylinder subjected to a triaxial test. The saw cut is inclined at an angle α to the cylinder axis, and the sand filling is intended to represent gouge. The triaxial test subjects the granite cylinder to a constant confining pressure and increasing axial stress to maintain a constant rate of shortening of the cylinder. The required axial stress increases at a decreasing rate to a maximum, beyond which a roughly constant axial stress is sufficient to maintain the constant rate of shortening. Such triaxial tests were run for saw cuts inclined at angles α of 20°, 25°, 30°, 35°, 40°, 45°, and 50° to the cylinder axis, and the apparent coefficient of friction μ a (ratio of the shear stress to the normal stress, both stresses resolved onto the saw cut) at failure was determined. Subject to the assumption that the observed failure involves slip on Coulomb shears (orientation unspecified), the orientation of the principal compression axis within the gouge can be calculated as a function of μ a for a given value of the coefficient of internal friction μ i . The rotation of the principal stress axes within the gouge in a triaxial test can then be followed as the shear strain across the gouge layer increases. For μ i ∼ 0.8, an appropriate value for highly sheared sand, the observed values μ a imply that the principal axis of compression within the gouge rotates so as to approach being parallel to the cylinder axis for all saw cut angles (20° < α < 50°). In the limiting state (principal compression axis parallel to cylinder axis) the stress state in the gouge layer would be the same as that in the granite cylinder, and the failure criterion would be independent of the saw cut angle.
We find that foreshocks provide clear evidence for an extended nucleation process before some earthquakes. In this study, we examine in detail the evolution of six California foreshock sequences, the 1986 Mount Lewis ( M L = 5.5), the 1986 Chalfant ( M L = 6.4), the 1986 Stone Canyon ( M L = 4.7), the 1990 Upland ( M L = 5.2), the 1992 Joshua Tree ( M W = 6.1), and the 1992 Landers ( M W = 7.3) sequence. Typically, uncertainties in hypocentral parameters are too large to establish the geometry of foreshock sequences and hence to understand their evolution. However, the similarity of location and focal mechanisms for the events in these sequences leads to similar foreshock waveforms that we cross correlate to obtain extremely accurate relative locations. We use these results to identify small-scale fault zone structures that could influence nucleation and to determine the stress evolution leading up to the mainshock. In general, these foreshock sequences are not compatible with a cascading failure nucleation model in which the foreshocks all occur on a single fault plane and trigger the mainshock by static stress transfer. Instead, the foreshocks seem to concentrate near structural discontinuities in the fault and may themselves be a product of an aseismic nucleation process. Fault zone heterogeneity may also be important in controlling the number of foreshocks, i.e., the stronger the heterogeneity, the greater the number of foreshocks. The size of the nucleation region, as measured by the extent of the foreshock sequence, appears to scale with mainshock moment in the same manner as determined independently by measurements of the seismic nucleation phase. We also find evidence for slip localization as predicted by some models of earthquake nucleation.
Numerical simulations of multiphase fluid and heat transport through a porous medium define combinations of rock properties and boundary conditions which lead to geyser‐like periodic discharge. Within the rather narrow range of conditions that allow geyser‐like behavior, eruption frequency and discharge are highly sensitive to the intrinsic permeabilities of the geyser conduit and the surrounding rock matrix, to the relative permeability functions assumed, and to pressure gradients in the matrix. In theory, heat pipes (concomitant upward flow of steam and downward flow of liquid) can exist under similar conditions, but our simulations suggest that the periodic solution is more stable. Simulated time series of geyser discharge are chaotic, but integrated quantities such as eruption frequency and mass discharge per eruption are free of chaos. These results may explain the observed sensitivity of natural geysers to small strains such as those caused by remote earthquakes, if ground motion is sufficient to induce permeability changes. Changes in geyser behavior caused by minor preseismic deformation, periodic surface loading, and Earth tides are more difficult to explain in the context of our current model.
Comparison of contemporaneous measurements of 84 distances in the range of 10 to 50 km by both Global Positioning System (GPS) and Geodolite (an electro-optical distance-measuring instrument) indicates that the Geodolite measurements are systematically longer by 0.283 ± 0.100 parts per million of the measured distance. Quoted uncertainty is 1 standard deviation. This amounts to 11.3 ± 4.0 mm at 40 km, which is near the maximum Geodolite range. The systematic difference is within the random uncertainty of an individual GPS-Geodolite comparison and was detected only from an analysis of a large number (84) of such comparisons. The source of the systematic difference has not been identified.
In thrusting and strike-slip situations, when the maximum principal horizontal stress S 1 acts nearly normal to a fault (a misoriented fault, such as the San Andreas), pore-fluid pressure > the lithostatic load, P f > S v , is required to reactivate movement on that fault. P f > S v may be achieved without causing hydraulic tensile fracturing if (1) previously existing cracks have regained cohesive strength by chemical processes, (2) subcritical crack growth has been blunted, and (3) the least principal horizontal stress S 3 nearly equals S v . Where P f > S v has been attained within a misaligned fault, increasing the stress difference ( S 1 – S 3 ) at constant P f > S v will not lead to shear failure, while a decrease in ( S 1 – S 3 ) can lead to shear failure of that fault. However, where the cohesive strength of material in a broad misaligned fault zone is less than that of the surrounding intact rock, increasing ( S 1 – S 3 ) while P f > S v can result in shear failure of fractures at near optimum angles to S 1 , but confined within this weak fault zone. If this faulting results in the local short-lived attainment of P f > S v (cataclastic deformation and frictional heating overcoming dilation) and a simultaneous decrease in ( S 1 – S 3 ), this combination of effects can trigger movement along the main trace of the misaligned fault. When increasing P f results in hydraulic failure, anisotropy in tensile strength or fracture toughness resulting from foliation within faults allows fractures to propagate along the planes of weakness rather than across the foliation perpendicular to S 3 .
New paleontologic and isotopic data from the Emma Creek and Marion Creek schists of the Coldfoot terrane, Arctic Alaska superterrane, central Brooks Range, suggest Devonian and possibly younger ages of deposition for their sedimentary protoliths. Conodonts from marble of the Emma Creek schist, intruded by a roughly 392 Ma orthogneiss, are late Lochkovian (early Early Devonian, between about 408 and 396 Ma) and Silurian to Devonian at two other locations. Spherical to oblong detrital zircons from quartz-mica schist of the overlying Marion Creek schist yield mostly discordant U—Pb data suggestive of provenance ages of 3.0, 2.0–1.8, and 1.5–1.4 Ga; however, several euhedral grains of zircon from Marion Creek quartz-mica schist have concordant U—Pb ages from 370 to 360 Ma. The Marion Creek schist in our study area therefore is at least 26 m.y. younger than the Emma Creek schist. The age data imply that the protolith of the Emma Creek schist is age correlative with Devonian carbonate rocks in the Hammond and North Slope terranes, whereas the Marion Creek schist is age correlative with Upper Devonian and Lower Mississippian clastic sedimentary rocks of the Endicott Group in the Endicott Mountains terrane and shale and carbonate units in the De Long Mountains and Sheenjek River terranes. Consequently, tectonic models restoring the entire Coldfoot terrane beneath partly or wholly coeval rocks of the Hammond, Endicott Mountains, De Long Mountains, and Sheenjek River terranes of the Arctic Alaska superterrane require revision. Alternative reconstructions, including restoration of the Coldfoot terrane inboard of the Endicott Mountains terrane or outboard of the De Long Mountains and Sheenjek River terranes are plausible but require either larger amounts of shortening than previously suggested or indicate problematic facies relations.
Glaciers, like rivers, have a branching structure which can be characterized by topological trees or networks. Probability distributions of various topological quantities in the networks are shown to satisfy the criterion for self-similarity, a symmetry structure which might be used to simplify future models of glacier dynamics. Two analytical methods of describing river networks, Shreve's random topology model and deterministic self-similar trees, are applied to the six glaciers of south central Alaska studied in this analysis. Self-similar trees capture the topological behavior observed for all of the glaciers, and most of the networks are also reasonably approximated by Shreve's theory.
The Endicott Mountains allochthon has been thrust far northward over the North Slope parautochthon in the northern Brooks Range. Progressively younger units are exposed northward within the allochthon. To the south, the incompetent Hunt Fork Shale has thickened internally by asymmetric folds and thrust faults. Northward, the competent Kanayut Conglomerate forms a duplex between a floor thrust in Hunt Fork and a roof thrust in the Kayak Shale. To the north, the competent Lisburne Group forms a duplex between a floor thrust in Kayak and a roof thrust in the Siksikpuk Formation. Both duplexes formed from north vergent detachment folds whose steep limbs were later truncated by south dipping thrust faults that only locally breach immediately overlying roof thrusts. Within the parautochthon, the Kayak, Lisburne, and Siksikpuk-equivalent Echooka Formation form a duplex identical to that in the allochthon. This duplex is succeeded abruptly northward by detachment folds in Lisburne. These folds are parasitic to an anticlinorium interpreted to reflect a fault-bend folded horse in North Slope “basement,” with a roof thrust in Kayak and a floor thrust at depth/These structures constitute two northward tapered, internally deformed wedges that are juxtaposed at the base of the allochthon. Within each wedge, competent units have been shortened independently between detachments, located mainly in incompetent units. The basal detachment of each wedge cuts upsection forward (northward) to define a wedge geometry within which units dip regionally forward. These dips reflect forward decrease in internal structural thickening by forward vergent folds and hindward dipping thrust faults.