Showing papers on "Fault (geology) published in 2010"

Geologically current plate motions

Abstract: SUMMARY We describe best-fitting angular velocities and MORVEL, a new closure-enforced set of angular velocities for the geologically current motions of 25 tectonic plates that collectively occupy 97 per cent of Earth's surface. Seafloor spreading rates and fault azimuths are used to determine the motions of 19 plates bordered by mid-ocean ridges, including all the major plates. Six smaller plates with little or no connection to the mid-ocean ridges are linked to MORVEL with GPS station velocities and azimuthal data. By design, almost no kinematic information is exchanged between the geologically determined and geodetically constrained subsets of the global circuit—MORVEL thus averages motion over geological intervals for all the major plates. Plate geometry changes relative to NUVEL-1A include the incorporation of Nubia, Lwandle and Somalia plates for the former Africa plate, Capricorn, Australia and Macquarie plates for the former Australia plate, and Sur and South America plates for the former South America plate. MORVEL also includes Amur, Philippine Sea, Sundaland and Yangtze plates, making it more useful than NUVEL-1A for studies of deformation in Asia and the western Pacific. Seafloor spreading rates are estimated over the past 0.78 Myr for intermediate and fast spreading centres and since 3.16 Ma for slow and ultraslow spreading centres. Rates are adjusted downward by 0.6–2.6 mm yr−1 to compensate for the several kilometre width of magnetic reversal zones. Nearly all the NUVEL-1A angular velocities differ significantly from the MORVEL angular velocities. The many new data, revised plate geometries, and correction for outward displacement thus significantly modify our knowledge of geologically current plate motions. MORVEL indicates significantly slower 0.78-Myr-average motion across the Nazca–Antarctic and Nazca–Pacific boundaries than does NUVEL-1A, consistent with a progressive slowdown in the eastward component of Nazca plate motion since 3.16 Ma. It also indicates that motions across the Caribbean–North America and Caribbean–South America plate boundaries are twice as fast as given by NUVEL-1A. Summed, least-squares differences between angular velocities estimated from GPS and those for MORVEL, NUVEL-1 and NUVEL-1A are, respectively, 260 per cent larger for NUVEL-1 and 50 per cent larger for NUVEL-1A than for MORVEL, suggesting that MORVEL more accurately describes historically current plate motions. Significant differences between geological and GPS estimates of Nazca plate motion and Arabia–Eurasia and India–Eurasia motion are reduced but not eliminated when using MORVEL instead of NUVEL-1A, possibly indicating that changes have occurred in those plate motions since 3.16 Ma. The MORVEL and GPS estimates of Pacific–North America plate motion in western North America differ by only 2.6 ± 1.7 mm yr−1, ≈25 per cent smaller than for NUVEL-1A. The remaining difference for this plate pair, assuming there are no unrecognized systematic errors and no measurable change in Pacific–North America motion over the past 1–3 Myr, indicates deformation of one or more plates in the global circuit. Tests for closure of six three-plate circuits indicate that two, Pacific–Cocos–Nazca and Sur–Nubia–Antarctic, fail closure, with respective linear velocities of non-closure of 14 ± 5 and 3 ± 1 mm yr−1 (95 per cent confidence limits) at their triple junctions. We conclude that the rigid plate approximation continues to be tremendously useful, but—absent any unrecognized systematic errors—the plates deform measurably, possibly by thermal contraction and wide plate boundaries with deformation rates near or beneath the level of noise in plate kinematic data.

Earthquake Fault Scaling: Self-Consistent Relating of Rupture Length, Width, Average Displacement, and Moment ReleaseEarthquake Fault Scaling: Self-Consistent Relating of Rupture Length, Width, Average Displacement

Abstract: In this paper, I propose the scaling relation W = C 1 L β (where β ≈2/3) to describe the scaling of rupture width with rupture length. I also propose a new displacement relation ![Graphic][1] , where A is the area ( LW ). By substituting these equations into the definition of seismic moment (![Graphic][2] ), I have developed a series of self-consistent equations that describe the scaling between seismic moment, rupture area, length, width, and average displacement. In addition to β , the equations have only two variables, C 1 and C 2, which have been estimated empirically for different tectonic settings. The relations predict linear log–log relationships, the slope of which depends only on β . These new scaling relations, unlike previous relations, are self-consistent, in that they enable moment, rupture length, width, area, and displacement to be estimated from each other and with these estimates all being consistent with the definition of seismic moment. I interpret C 1 as depending on the size at which a rupture transitions from having a constant aspect ratio to following a power law and C 2 as depending on the displacement per unit area of fault rupture and so static stress drop. It is likely that these variables differ between tectonic environments; this might explain much of the scatter in the empirical data. I suggest that these relations apply to all faults. For small earthquakes ( M ∼5) earthquakes β =2/3, so L 2.5 applies. For dip-slip earthquakes this scaling applies up to the largest events. For very large ( M >∼7.2) strike-slip earthquakes, which are fault width-limited, β =0 and assuming ![Graphic][3] , then L 1.5 scaling applies. In all cases, M ∝ A 1.5 fault scaling applies. [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif

Earthquake Monitoring in Southern California for Seventy-Seven Years (1932-2008)

Abstract: The Southern California Seismic Network (SCSN) has produced the SCSN earthquake catalog from 1932 to the present, a period of more than 77 yrs. This catalog consists of phase picks, hypocenters, and magnitudes. We present the history of the SCSN and the evolution of the catalog, to facilitate user understanding of its limitations and strengths. Hypocenters and magnitudes have improved in quality with time, as the number of stations has increased gradually from 7 to ~400 and the data acquisition and measuring procedures have become more sophisticated. The magnitude of completeness (M_c) of the network has improved from M_c ~3.25 in the early years to M_c ~1.8 at present, or better in the most densely instrumented areas. Mainshock–aftershock and swarm sequences and scattered individual background earthquakes characterize the seismicity of more than 470,000 events. The earthquake frequency-size distribution has an average b-value of ~1.0, with M≥6.0 events occurring approximately every 3 yrs. The three largest earthquakes recorded were 1952 M_w 7.5 Kern County, 1992 M_w 7.3 Landers, and 1999 M_w 7.1 Hector Mine sequences, and the three most damaging earthquakes were the 1933 M_w 6.4 Long Beach, 1971 M_w 6.7 San Fernando, and 1994 M_w 6.7 Northridge earthquakes. All of these events ruptured slow-slipping faults, located away from the main plate boundary fault, the San Andreas fault. Their aftershock sequences constitute about a third of the events in the catalog. The fast slipping southern San Andreas fault is relatively quiet at the microseismic level and has not had an M>6 earthquake since 1932. In contrast, the slower San Jacinto fault has the highest level of seismicity, including several M>6 events. Thus, the spatial and temporal seismicity patterns exhibit a complex relationship with the plate tectonic crustal deformation.

The application of failure mode diagrams for exploring the roles of fluid pressure and stress states in controlling styles of fracture‐controlled permeability enhancement in faults and shear zones

Abstract: Geofluids (2010) 10, 217–233 Abstract Permeability enhancement associated with deformation processes in faults and shear zones plays a key role in facilitating fluid redistribution between fluid reservoirs in the crust. Especially in high fluid flux hydrothermal systems, fracture-controlled permeability can be relatively short-lived, unless it is repeatedly regenerated by ongoing deformation. Failure mode diagrams in pore fluid factor and differential stress space, here termed λ–σ failure mode diagrams, provide a powerful tool for analysing how fluid pressure and stress states drive failure, associated permeability enhancement and vein styles during deformation in faults and shear zones. During fault-valve behaviour in the seismogenic regime, relative rates of recovery of pore fluid factor, differential stress and fault cohesive strength between rupture events impact on styles of veining and associated, fracture-controlled permeability enhancement in faults and shear zones. Examples of vein-rich fault zones are used to illustrate how constraints can be placed, not just on fluid pressure and stress states at failure, but also on the fluid pressurization and loading paths associated with failure and transitory permeability enhancement in faults and shear zones. This provides insights about when, during the fault-valve cycle, various types of veins can form. The use of failure mode diagrams also provides insights about the relative roles of optimally oriented faults and misoriented faults as hydraulically conductive structures. The analysis highlights the dynamics of competition between fluid pressures and loading rates in driving failure and repeated permeability regeneration in fracture-controlled, hydrothermal systems.

Complex rupture during the 12 January 2010 Haiti earthquake

Abstract: Initially, the devastating Mw 7.0, 12 January 2010 Haiti earthquake seemed to involve straightforward accommodation of oblique relative motion between the Caribbean and North American plates along the Enriquillo‐Plantain Garden fault zone. Here, we combine seismological observations, geologic field data and space geodetic measurements to show that, instead, the rupture process involved slip on multiple faults. Primary surface deformation was driven by rupture on blind thrust faults with only minor, deep, lateral slip along or near the main Enriquillo‐Plantain Garden fault zone; thus the event only partially relieved centuries of accumulated left-lateral strain on a small part of the plate-boundary system. Together with the predominance of shallow off-fault thrusting, the lack of surface deformation implies that remaining shallow shear strain will be released in future surface-rupturing earthquakes on the Enriquillo‐Plantain Garden fault zone, as occurred in inferred Holocene and probable historic events. We suggest that the geological signature of this earthquake—broad warping and coastal deformation rather than surface rupture along the main fault zone—will not be easily recognized by standard palaeoseismic studies. We conclude that similarly complex earthquakes in tectonic environments that accommodate both translation and convergence—such as the San Andreas fault through the Transverse Ranges of California—may be missing from the prehistoric earthquake record.

Seismic imaging of the lithospheric structure of the Zagros mountain belt (Iran)

Abstract: We present a synthesis and a comparison of the results of two temporary passive seismic experiments installed for a few months across Central and Northern Zagros. The receiver function analysis of teleseismic earthquake records gives a high-resolution image of the Moho beneath the seismic transects. On both cross-sections, the crust has an average thickness of 43±2 km beneath the Zagros fold-and-thrust belt and the Central domain. The crust is thicker beneath the hanging wall of the Main Zagros Reverse Fault (MZRF), with a larger maximum Moho depth in Central (69±2 km) than in Northern Zagros (56±2 km). The thickening affects a narrower region (170 km) beneath the Sanandaj-Sirjan zone of Central Zagros and a wider region (320 km) in Northern Zagros. We propose that this thickening is related to overthrusting of the crust of the Arabian margin by the crust of Central Iran along the MZRF considered as a major thrust fault cross-cutting the whole crust. The fault is imaged as a lowvelocity layer in the receiver function data of the Northern Zagros profile. Moreover, the crustal-scale thrust model reconciles the imaged seismic Moho with the Bouguer anomaly data measured on the Central Zagros transect. At upper mantle depth, P-wave tomography confirms the previously observed strong contrast between the faster velocities of the Arabian margin and the lower velocities of the Iranian micro-blocks. Our higher-resolution tomography combined with surface-wave analysis locates the suture in the shallow mantle of the Sanandaj-Sirjan zone beneath Central Zagros. The Arabian upper-mantle has shield-like shear-wave velocities, while the lower velocities of the Iranian upper mantle are likely due to higher temperature. But these velocities are not low enough and the low-velocity layer not thick enough to conclude on a delamination of the lithospheric mantle lid beneath Iran. The lack of a high-velocity anomaly in the mantle beneath Central Iran suggests that the Neotethyan oceanic lithosphere is probably detached from the Arabian margin.

Cenozoic right-slip faulting along the eastern margin of the Pamir salient, northwestern China

Abstract: The Pamir salient defines the western end of the Himalayan-Tibetan orogen and has overthrust the Tarim-Tajik basin to the north by ∼300 km along a late Cenozoic, south-dipping intracontinental subduction zone. Field mapping, structural measurements, and analysis of mesoscale structures along a 32-km-long reach of the Yarkand River document the tectonic evolution of the east flank of this salient, between the North Pamir to the west and the Western Kunlun Shan to the east. The study area is cut by a set of four, north-northwest–striking, steeply dipping brittle faults. Microstructures and asymmetric outcrop- to map-scale folds indicate right slip along these faults. Between these structures, fault-bounded panels of Phanerozoic strata are deformed by en echelon folds with axes that trend more westerly than the adjacent faults, consistent with dextral transpression. The fault system described here extends for ∼350 km along the eastern flank of the Pamir salient. Transpressional right slip along this set of faults, here called the Kashgar-Yecheng transfer system, appears to have accommodated late Cenozoic separation of the North Pamir from the Western Kunlun Shan during south-directed intracontinental subduction beneath the leading edge of the Pamir salient. Correlation of major faults suggests total slip along the Kashgar-Yecheng transfer system is likely on the order of ∼280 km. This offset estimate implies long-term slip rates of 7–15 mm/a along the Kashgar-Yecheng transfer system when combined with previous sedimentologic, stratigraphic, and thermochronologic data that indicate deformation along the east flank of the Pamir started between the late Eocene and early Miocene. These results imply that the first-order structures on the western and eastern flanks of the Pamir are asymmetric: previous work has shown that deformation in the west was accommodated by anticlockwise vertical axis rotation of the Pamir over the eastern margin of the Tajik basin. This rotation is generally interpreted to reflect northwest-directed radial thrusting, in contrast to the transpressional right-slip transfer faulting on the east side reported here.

Extension on the Tibetan plateau: recent normal faulting measured by InSAR and body wave seismology

Abstract: We use InSAR and body wave modelling to determine the faulting parameters for a series of five Mw 5.9–7.1 normal faulting earthquakes that occurred during 2008, including the March 20 Yutian earthquake (Mw 7.1), one of the largest normal faulting events to have occurred recently on the continents. We also study three earlier normal faulting earthquakes that occurred in southern Tibet between 1992 and 2005. Coseismic deformation for each of these eight events is measured with ascending and descending interferograms from ENVISAT, ERS and ALOS SAR data. Elastic dislocation modelling of the line-of-sight InSAR displacements and body wave seismological modelling of P and SH waves are used to estimate fault parameters and are found to be in good agreement for all the events studied. The use of InSAR to measure deformation allows a relatively precise determination of the fault location in addition to resolving the focal plane ambiguity. Only five of the eight events are associated with a clear surface topographic feature, suggesting that an underestimation of the amount of extension would result from using the surface expressions of normal faulting alone. The observations, in all cases, are consistent with slip on planar surfaces, with dips in the range 40–50°, that penetrates the uppermost crust to a depth of 10–15 km. We find no evidence for active low-angle (dip less than 30°) normal faulting. The contribution of the normal faulting to overall extension estimated by summing seismic moments over earthquakes for the past 43 yr is 3–4 mm yr−1, or 15–20 per cent of the rates of extension measured across the plateau using GPS. 85 per cent of the moment release in normal faulting over the past 43 yr has occurred in regions whose surface height exceeds 5 km. This observation adds weight to the suggestions that the widespread normal faulting on the plateau is the result of variations in the gravitational potential energy of the lithosphere.

Faults and fault properties in hydrocarbon flow models

Abstract: Geofluids (2010) 10, 94–113 Abstract The petroleum industry uses subsurface flow models for two principal purposes: to model the flow of hydrocarbons into traps over geological time, and to simulate the production of hydrocarbon from reservoirs over periods of decades or less. Faults, which are three-dimensional volumes, are approximated in both modelling applications as planar membranes onto which predictions of the most important fault-related flow properties are mapped. Faults in porous clastic reservoirs are generally baffles or barriers to flow and the relevant flow properties are therefore very different to those which are important in conductive fracture flow systems. A critical review and discussion is offered on the work-flows used to predict and model capillary threshold pressure for exploration fault seal analysis and fault transmissibility multipliers for production simulation, and of the data from which the predictions derive. New flow simulation models confirm that failure of intra-reservoir sealing faults can occur during a reservoir depressurization via a water-drive mechanism, but contrary to anecdotal reports, published examples of production-induced seal failure are elusive. Ignoring the three-dimensional structure of fault zones can sometimes have a significant influence on production-related flow, and a series of models illustrating flow associated with relay zones are discussed.

Relation between continental strike-slip earthquake segmentation and thickness of the crust

Abstract: [1] High-resolution maps of large continental strike-slip earthquake surface ruptures show that they are formed of fault segments. These segments are bounded by fault bends, step overs, or combinations of the two. The lowest limit in size for such segments may not be relevant in the understanding of earthquake mechanics, as it pertains to the granular properties of fault zones. The maximum limit in segment length, however, is important as it is directly relates to the maximum extent of seismic rupture. To measure the length of the segments, a new quantitative method based on piecewise linear fitting is developed and is used to automatically retrieve segments from earthquake rupture maps. Next, this approach is tested against a set of ten continental strike-slip earthquake ruptures derived from similar, high quality maps. The test suggests that segments have a maximum length of ∼18 km, independent of regional tectonic setting. Slip-inversions for earthquakes, based on seismological and/or geodetic data, most often are not unique and can show some variability even for one particular event. Some basic characteristics, however, such as total moment release or general source geometry, seem to persist that are relevant to earthquake mechanics. Measurements of the maximum horizontal extent of individual slip-patches derived from seismic source inversion for strike-slip ruptures show that their strike dimension does not increase infinitely with magnitude, but instead reaches a maximum value of ∼25 km. These two independent lines of observations, complemented by earlier data and analog experiments, suggest that it is the thickness of the seismogenic crust that controls the structural scaling of the length of seismic segments, and that it is independent of the ultimate size of individual earthquakes.

The Yakutat terrane: Dramatic change in crustal thickness across the Transition fault, Alaska

Abstract: We present new constraints on the crustal structure of the Yakutat terrane and evidence of the role of the Transition fault in southern Alaska. The Yakutat terrane south of Yakutat Bay includes crystalline crust that is 24–27 km thick overlain by sedimentary units that are 4.5–7.5 km thick. The Yakutat terrane crustal thickness and velocity structure are consistent with an oceanic plateau origin. The southern edge of the Yakutat terrane is bounded by the Transition fault, which is imaged as a near-vertical fault zone ∼1 km wide. The Transition fault is coincident with a dramatic change in Moho depth from 32 km for Yakutat oceanic plateau crust to 11.5 km for Pacific Ocean crust occurring over a horizontal distance of 0–5 km. There is no evidence for underthrusting of the Pacific Ocean crust beneath the Yakutat terrane at the Transition fault. We argue that the Yakutat terrane formed on the Kula or Farallon plate and was later juxtaposed next to younger Pacific Ocean crust by the Transition fault.

Triggering of New Madrid seismicity by late-Pleistocene erosion

Abstract: The New Madrid seismic zone, in a now heavily populated area of the central United States, was responsible for the 1811–12 New Madrid earthquakes of magnitudes of 7 or greater. The extent of the current seismic hazard in the region is hotly debated. Eric Calais and colleagues present evidence that the geologically recent sequence of large earthquakes in this region was triggered by the rapid removal of sediments by the rivers of the northern Mississippi embayment at the end of the last ice age. Models indicate that fault segments that have already ruptured are unlikely to fail again soon, but stress changes from sediment unloading and previous earthquakes may eventually be sufficient to bring to failure other nearby segments that have not yet ruptured, indicating that the hazard may be more widespread than previously thought. These authors argue that the concentration of magnitude-7 or larger earthquakes in the New Madrid seismic zone since the end of the last ice age results from the recent, climate-controlled, erosional history of the northern Mississippi embayment. They show that the upward flexure of the lithosphere caused a reduction of normal stresses in the upper crust sufficient to unclamp pre-existing faults close to failure equilibrium. The spatiotemporal behaviour of earthquakes within continental plate interiors is different from that at plate boundaries. At plate margins, tectonic motions quickly reload earthquake ruptures, making the location of recent earthquakes and the average time between them consistent with the faults’ geological, palaeoseismic and seismic histories. In contrast, what determines the activation of a particular mid-continental fault and controls the duration of its seismic activity remains poorly understood1. Here we argue that the concentration of magnitude-7 or larger earthquakes in the New Madrid seismic zone of the central United States2,3 since the end of the last ice age results from the recent, climate-controlled, erosional history of the northern Mississippi embayment. We show that the upward flexure of the lithosphere caused by unloading from river incision between 16,000 and 10,000 years ago caused a reduction of normal stresses in the upper crust sufficient to unclamp pre-existing faults close to failure equilibrium. Models indicate that fault segments that have already ruptured are unlikely to fail again soon, but stress changes from sediment unloading and previous earthquakes may eventually be sufficient to bring to failure other nearby segments that have not yet ruptured.

Coseismic ground deformation of the 6 April 2009 L’Aquila earthquake (central Italy, Mw6.3)

Abstract: [1] We provide field data of coseismic ground deformation related to the 6 April Mw 63 L'Aquila normal faulting earthquake Three narrow fracture zones were mapped: Paganica-Colle Enzano (P-E), Mt Castellano-Mt Stabiata (C-S) and San Gregorio (SG) These zones define 13 km of surface ruptures that strike at 130–140° We mapped four main types of ground deformation (free faces on bedrock fault scarps, faulting along synthetic splays and fissures with or without slip) that are probably due to the near-surface lithology of the fault walls and the amount of slip that approached the surface coseismically The P-E and C-S zones are characterized by downthrow to the SW (up to 10 cm) and opening (up to 12 cm), while the SG zone is characterized only by opening Afterslip throw rates of 05–06 mm/day were measured along the Paganica fault, where paleoseismic evidence reveals recurring paleo-earthquakes and post-248 kyr slip-rate ≥ 024 mm/yr

Nanocoatings of clay and creep of the San Andreas fault at Parkfield, California

Abstract: Mudrock samples were investigated from two fault zones at ∼3066 m and ∼3296 m measured depth (MD) located outside and within the main damage zone of the San Andreas Fault Observatory at Depth (SAFOD) drillhole at Parkfield, California. All studied fault rocks show features typical of those reported across creep zones with variably spaced and interconnected networks of polished displacement surfaces coated by abundant polished films and occasional striations. Electron microscopy and X-ray diffraction study of the surfaces reveal the occurrence of neocrystallized thin film clay coatings containing illite-smectite (I-S) and chlorite-smectite (C-S) minerals. 40 Ar/ 39 Ar dating of the illitic mix-layered coatings demonstrated Miocene to Pliocene crystallization and revealed an older fault strand (8 ± 1.3 Ma) at 3066 m MD, and a probably younger fault strand (4 ± 4.9 Ma) at 3296 m MD. Today, the younger strand is the site of active creep behavior, reflecting a possible (re)activation of these clay-weakened zones. We propose that the majority of slow fault creep is controlled by the high density of thin (

The West Andean Thrust, the San Ramón Fault, and the seismic hazard for Santiago, Chile

Abstract: [1] The importance of west verging structures at the western flank of the Andes, parallel to the subduction zone, appears currently minimized. This hampers our understanding of the Andes-Altiplano, one of the most significant mountain belts on Earth. We analyze a key tectonic section of the Andes at latitude 33.5°S, where the belt is in an early stage of its evolution, with the aim of resolving the primary architecture of the orogen. We focus on the active fault propagation–fold system in the Andean cover behind the San Ramon Fault, which is critical for the seismic hazard in the city of Santiago and crucial to decipher the structure of the West Andean Thrust (WAT). The San Ramon Fault is a thrust ramp at the front of a basal detachment with average slip rate of ∼0.4 mm/yr. Young scarps at various scales imply plausible seismic events up to Mw 7.4. The WAT steps down eastward from the San Ramon Fault, crossing 12 km of Andean cover to root beneath the Frontal Cordillera basement anticline, a range ∼5 km high and >700 km long. We propose a first-order tectonic model of the Andes involving an embryonic intracontinental subduction consistent with geological and geophysical observations. The stage of primary westward vergence with dominance of the WAT at 33.5°S is evolving into a doubly vergent configuration. A growth model for the WAT-Altiplano similar to the Himalaya-Tibet is deduced. We suggest that the intracontinental subduction at the WAT is a mechanical substitute of a collision zone, rendering the Andean orogeny paradigm obsolete.

First global positioning system results in northern Myanmar: Constant and localized slip rate along the Sagaing fault

Abstract: We present the first recorded global positioning system (GPS) data from Myanmar measured at the northern tip of the Sagaing fault. This area is in a very complex geodynamic context, where rigid and semirigid plates interact. The 12 GPS sites measured in 2005 and 2008 in northern Myanmar show that the slip rate is 18 mm/yr and is localized along a single active narrow fault trace (<20 km wide). The same rate was previously demonstrated and remeasured, 500 km southward, in central Myanmar. Despite the geodynamic regional complexity induced by interaction between the Sunda and India plates, the Burma microplate, and the highly deformable eastern Himalayan syntaxis, the slip rate remains surprisingly constant along this fault. However, the modeled locking depth varies from 20 km in central Myanmar to 5 km in the north.