Systematic deflection and offset of the Yangtze River drainage system along the strike-slip Ganzi-Yushu-Xianshuihe Fault Zone,Tibetan Plateau

The Ganzi-Yushu-Xianshuihe Fault Zone (GYXFZ) is a typical active strike-slip fault that has triggered many large historic earthquakes, including the 2010 M_w 6.9 Yushu earthquake in the central Tibetan Plateau. This fault zone extends for ca. 800 km from the central Tibetan Plateau to its southeastern margin and varies in trend from WNW-ESE in the northwestern segment of the fault zone to NNW-SSE in the southeastern segment, having the geometry of an arc projecting northeastwards. In this study, we present evidence for the systematical sinistral deflection and/or offset of the Yangtze River and its branch stream channels and valleys along the GYXFZ. Topographic analysis of three-dimensional (3D) perspective images constructed using digital elevation model (DEM) data, 0.5 m-resolution WorldView and GeoEye images, and 15 m-resolution Landsat-Enhanced Thematic Mapper (ETM+) images, together with analysis of geological structures, reveals the following: (i) the main river channels and valleys of the Yangtze River drainage system show systematic sinistral deflections and/or offsets along the GYXFZ; (ii) various amounts of sinistral offset have accumulated on the tributary stream channels, valleys, and gullies of the Yangtze River along the fault, with a linear relation, D = aL, between the upstream length L from the deflected point and the offset amount D with a certain coefficient a; (iii) the maximum amount of sinistral offset is up to ca. 60 km, which was accumulated in the past 13–5 Ma; and (iv) the long-term average strike-slip rate is ca. 4.6–12 mm/year. Geological and geomorphic evidence, combined with geophysical data, demonstrates that the GYXFZ is currently active as one of the major seismogenic faults in the Tibetan Plateau, dominated by left-lateral strike-slip motion. Our findings supply important evidence for the tectonic evolution of strike-slip faults in the Tibetan Plateau since the Eurasia-India continental collision.     

Neotectonics of the Dinarides–Pannonian Basin transition and possible earthquake sources in the Banja Luka epicentral area

This study provides evidence for post-5 Ma shortening in the transition area between the Dinarides fold-and-thrust belt and the Pannonian Basin and reviews possible earthquake sources for the Banja Luka epicentral area (northern Bosnia and Herzegovina) where the strongest instrumentally recorded earthquake (M_L 6.4) occurred on 27 October 1969. Geological, geomorphological and reflection seismic data provide evidence for a contractional reactivation of Late Palaeogene to Middle Miocene normal faults at slip rates below 0.1 mm/a. This reactivation postdates deposition of the youngest sediments in the Pannonian Basin of Pontian age (c. 5 Ma). Fault plane solutions for the main 1969 Banja Luka earthquake (M_L 6.4) and its largest foreshock (M_L 6.0) indicate reverse faulting along ESE–WNW-striking nodal planes and generally N–S trending pressure axes. The spatial distribution of epicentres and focal depths, analyses of the macroseismic field and fault-plane solutions for several smaller events suggest on-going shortening in the internal Dinarides. Seismic deformation of the upper crust is also associated with strike-slip faults, likely related to the NE–SW trending, sinistral Banja Luka fault. Possibly, this fault transfers contraction between adjacent segments of the Dinarides thrust system. The study area represents the seismically most active region of the Dinarides apart from the Adriatic Sea coast and the bend zone around Zagreb. We propose that on-going thrusting in the internal Dinarides thrust system takes up a portion of the current Adria–Europe convergence.     

Quaternary fault kinematics and chronology in intraplate northeastern Brazil

Northeastern Brazil has experienced earthquake swarms, which include events up to m_b = 5.2 in the last 30 years. Ground-rupturing events, however, have not been reported in the last 200 years. We have used a multidisciplinary approach to describe the geometry, kinematics, age, paleostress field, and paleoseismological significance of ground ruptures filled by clastic sediments in the region. Methods included GIS- and remote sensing-assisted structural mapping, borehole and resistivity sounding, grain size and X-ray analysis, and Optically Stimulated Luminescence (OSL) and Single-Aliquot Regenerative-Dose (SAR) chronology. Sediment-filled faults concentrate on the hanging wall of the Jundiaí fault, a 35 km long, NE-trending, late Cretaceous to Quaternary fault. These sediment-filled faults cut across Precambrian granites and overlying alluvial and colluvial deposits and exhibit steeply planar geometry. The largest faults exhibit sediment fills that extend downward to at least 20 m. The faults alternate between normal and right-oblique normal slip, which allowed the gravitational influx of unlithified gravels to gravelly sediments into the granites. Field evidence indicates seismogenic origin associated with multistage opening and infilling. On the basis of OSL and SAR chronology, we distinguish six generations of infills, which represent the main periods of fault activity: 8.0–9.0, 11.0–15.0, 16.0–24.0, 37.0–45.5, 65.0–67.9, and 84.5–93.5 ka. These dates suggest a 15.8 ka recurrence period. A few OSL dates may represent maximum ages due to poor bleaching of sediments, which implies that the recurrence period may be underestimated. We conclude that the Jundiaí fault has been continuously active in respect of surface-rupturing during at least the last ～100 ka. The size of the surface ruptures is consistent with earthquake magnitude M ～ 5.3, close to magnitudes observed in northeastern Brazil in the short instrumental and historical record. We suggest that seismogenic-fault recurrence in intraplate settings is longer than human settlement and should be assessed by multidisciplinary methods, mainly those that provide subsurface data and chronology of clastic sediments. Faults that trap sediments during movements are ideal places for these studies. The structures we describe here may have analogs in both modern and ancient intraplate settings.     

New observations on the 1939 Erzincan Earthquake surface rupture on the Kelkit Valley segment of the North Anatolian Fault Zone,Turkey

The 1939 Erzincan Earthquake (M = 7.8), occurred on the North Anatolian Fault Zone (NAFZ), was one of the most active strike-slip faults in the world, and created a 360-km-long surface rupture. Traces of this surface rupture are still prominently observed. In the absence of detailed mapping to resolve the fault characteristics, detailed observations have been conducted at 20 different points on the 70-km-long Kelkit Valley Segment (KVS) of the NAFZ's between Niksar and Koyulhisar. Field data defining fault character and slip amounts were found at eight points and show right-lateral slip varying between 1.8 and 4.25 m and the vertical slip varying between 0.5 and 2.0 m.The KVS developed in the most morphologically prominent and narrowest part of the NAFZ. Therefore, the chances of finding evidence of more than one historical earthquake in trenches opened to investigate palaeoseismological aspects are higher. Faults observed in foundation and channel excavations opened for energy purposes in the Reşadiye region show this clearly and evidence for up to four seismic events including the 1939 Erzincan Earthquake have been discovered. Further studies are required to discover whether right-lateral deformation on at some locations on this segment is surface ruptures associated with the 1939 earthquake or later creep.     

Inferred fault geometry and slip distribution of the 2010 Jiashian,Taiwan, earthquake is consistent with a thick-skinned deformation model

We invert measurements of coseismic displacements from 139 continuously recorded GPS sites from the 2010, Jiashian, Taiwan earthquake to solve for fault geometry and slip distribution using an elastic uniform stress drop inversion. The earthquake occurred at a depth of ~ 23 km in an area between the Western Foothills fold-and-thrust belt and the crystalline high mountains of the Central Range, providing an opportunity to examine the deep fault structure under Taiwan. The inferred rupture plane is oblique to the prominent orientation of thrust faults and parallel to several previously recognized NW-striking transfer zones that appear to connect stepping thrusts. We find that a fault striking 318°–344° with dip of 26°–41° fits the observations well with oblique reverse-sinistral slip under a low stress drop of about 0.5 MPa. The derived geodetic moment of 2.92 × 10¹⁸ N-m is equivalent to a M_w = 6.24 earthquake. Coseismic slip is largely concentrated within a circular patch with a 10-km radius at the depth between 10 and 24 km and maximum slip of 190 mm. We suggest this earthquake ruptured the NW-striking Chishan transfer fault zone, which we interpret as a listric NE-dipping lateral ramp with oblique slip connecting stepping thrust faults (ramps). The inferred slip on the lateral ramp is considerably deeper than the 7–15 km deep detachment identified in previous studies of western Taiwan. We infer an active basal detachment under western Taiwan at a depth of at least ~ 20–23 km based on these inversion results. The earthquake may have nucleated at the base of the lateral ramp near the intersection with the basal detachment. Coulomb stress change calculations suggest that this earthquake moved several NE-striking active thrust faults in western Taiwan nearer to failure.     

Active tectonics of the Northern Rif (Morocco) from geomorphic and geochronological data

We present results of a geomorphological and morphotectonic analysis of the northeastern part of the Rif. We show that the present day kinematics of the Rif is characterized by active deformation along the Trougout and Nekor faults in the North-East. Digital Elevation Models of offset drainage features (streams, fluvial terraces) allow determining a normal-left-lateral motion along the Trougout fault and a left-lateral strike-slip motion along the Nekor fault. Preliminary ³He cosmogenic dates of tectonic markers yield vertical and horizontal slip rates of ∼0.9 mm/yr and ∼0.5 mm/yr, respectively along the Trougout fault. The present-day localized transtension seen in the north-eastern Rif morphology (Ras Tarf) is coeval with uplifted marine terraces near the Al Hoceima Bay. U/Th dating of shells yield an average uplift rate of ∼0.2 mm/yr during the past 500 ka. These data show that active transtension in the northeastern Rif is also associated with uplift. These new morphotectonic constraints are consistent with the GPS measurements showing southwestward overall motion of most of the Rif belt with respect to stable Africa.     

Spatiotemporal patterns of fault slip rates across the Central Sierra Nevada frontal fault zone

Patterns in fault slip rates through time and space are examined across the transition from the Sierra Nevada to the Eastern California Shear Zone–Walker Lane belt. At each of four sites along the eastern Sierra Nevada frontal fault zone between 38 and 39° N latitude, geomorphic markers, such as glacial moraines and outwash terraces, are displaced by a suite of range-front normal faults. Using geomorphic mapping, surveying, and ¹⁰Be surface exposure dating, mean fault slip rates are defined, and by utilizing markers of different ages (generally, ~ 20 ka and ~ 150 ka), rates through time and interactions among multiple faults are examined over 10⁴–10⁵ year timescales.At each site for which data are available for the last ~ 150 ky, mean slip rates across the Sierra Nevada frontal fault zone have probably not varied by more than a factor of two over time spans equal to half of the total time interval (~ 20 ky and ~ 150 ky timescales): 0.3 ± 0.1 mm year^? 1 (mode and 95% CI) at both Buckeye Creek in the Bridgeport basin and Sonora Junction; and 0.4 + 0.3/?0.1 mm year^? 1 along the West Fork of the Carson River at Woodfords. Data permit rates that are relatively constant over the time scales examined. In contrast, slip rates are highly variable in space over the last ~ 20 ky. Slip rates decrease by a factor of 3–5 northward over a distance of ~ 20 km between the northern Mono Basin (1.3 + 0.6/?0.3 mm year^? 1 at Lundy Canyon site) to the Bridgeport Basin (0.3 ± 0.1 mm year^? 1). The 3-fold decrease in the slip rate on the Sierra Nevada frontal fault zone northward from Mono Basin is indicative of a change in the character of faulting north of the Mina Deflection as extension is transferred eastward onto normal faults between the Sierra Nevada and Walker Lane belt.A compilation of regional deformation rates reveals that the spatial pattern of extension rates changes along strike of the Eastern California Shear Zone-Walker Lane belt. South of the Mina Deflection, extension is accommodated within a diffuse zone of normal and oblique faults, with extension rates increasing northward on the Fish Lake Valley fault. Where faults of the Eastern California Shear Zone terminate northward into the Mina Deflection, extension rates increase northward along the Sierra Nevada frontal fault zone to ~ 0.7 mm year^? 1 in northern Mono Basin. This spatial pattern suggests that extension is transferred from more easterly fault systems, e.g., Fish Lake Valley fault, and localized on the Sierra Nevada frontal fault zone as the Eastern California Shear Zone–Walker Lane belt faulting is transferred through the Mina Deflection.     

Active tectonics of the outer northern Apennines: Adriatic vs. Po Plain seismicity and stress fields

The recent earthquake sequences of 2012 (northern Italy) and 2013 (Marche offshore) provided new, fundamental constraints to the active tectonic setting of the outer northern Apennines. In contrast to the Po Plain area, where the 2012 northern Italy earthquakes confirmed active frontal thrusting, the new focal mechanisms obtained in this study for the 2013 Marche offshore earthquakes indicate that only minor thrust fault reactivation occurs in the Adriatic domain, even for a theoretically favourably oriented maximum horizontal compression. Recent seismicity in this domain appears to be mainly controlled by transcurrent crustal faults dissecting the Apennine thrust belt. The along-strike stress field variation from the Po Plain to the Adriatic area has been quantitatively investigated by applying the multiple inverse method (MIM) to the analysis of the entire seismicity recorded from January 1976 to August 2014, from the top 12 km of the crust (fault plane solutions from 127 earthquakes with M_W ≥ 4), allowing us to obtain a comprehensive picture of the state of stress over the outer zone of the fold and thrust belt. The present-day stress field has been defined for 39 cells of 1.5° × 1.5° surface area and 12 km depth. The obtained stress field maps point out that, although the entire outer northern Apennines belt is characterized by a sub-horizontal maximum compressive axis (σ₁), the minimum compression (σ₃) is sub-vertical only in the Po Plain area, becoming sub-horizontal in the Adriatic sector, thus confirming that the latter region is dominated by an active tectonic regime of strike-slip type.     

The Al Hoceima Mw 6.4 earthquake of 24 February 2004 and its aftershocks sequence

The Al Hoceima Mw 6.4 earthquake of 24 February 2004 that occurred in the eastern Rif region of Morocco already hit by a large event in May 1994 (Mw 5.9) has been followed by numerous aftershocks in the months following the event. The aftershock sequence has been monitored by a temporary network of 17 autonomous seismic stations during 15 days (28 March–10 April) in addition to 5 permanent stations of the Moroccan seismic network (CNRST, SPG, Rabat). This network allowed locating accurately about 650 aftershocks that are aligned in two directions, about N10-20E and N110-120E, in rough agreement with the two nodal planes of the focal mechanism (Harvard). The aftershock alignments are long enough, about 20 km or more, to correspond both to the main rupture plane. To further constrain the source of the earthquake main shock and aftershocks (mb > 3.5) have been relocated thanks to regional seismic data from Morocco and Spain. While the main shock is located at the intersection of the aftershock clouds, most of the aftershocks are aligned along the N10-20E direction. This direction together with normal sinistral slip implied by the focal mechanism is similar with the direction and mechanisms of active faults in the region, particularly the N10E Trougout oblique normal fault. Indeed, the Al Hoceima region is dominated by an approximate ENE-SSW direction of extension, with oblique normal faults. Three major 10–30 km-long faults, oriented NNE-SSW to NW-SE are particularly clear in the morphology, the Ajdir and Trougout faults, west and east of the Al Hoceima basin, respectively, and the NS Rouadi fault 20 km to the west. These faults show clear evidence of recent vertical displacements during the late Quaternary such as uplifted alluvial terraces along Oued Rihs, offset fan surfaces by the Rouadi fault and also uplifted and tilted abandoned marine terraces on both sides of the Al Hoceima bay.However, the N20E direction is in contrast with seismic sources identified from geodetic inversions, which favour but not exclusively the N110-120E rupture directions, suggesting that the 1994 and 2004 events occurred on conjugate faults. In any event, the recent seismicity is thus concentrated on sinistral N10-20E or N110-120E dextral strike-slip faults, which surface expressions remain hidden below the 3–5 km-thick Rif nappes, as shown by the tomographic images build from the aftershock sequence and the concentration of the seismicity below 3 km. These observations may suggest that strain decoupling between the thrusted cover and the underlying bedrock and highlights the difficulty to determine the source properties of moderate events with blind faults even in the case of good quality recorded data.     

Main crustal discontinuities of Morocco derived from gravity data

Sharp linear gradients in maps of potential field data are generally assumed to result from sharp discontinuities or boundaries between rocks having different densities or magnetic susceptibilities and are usually associated with faults or other geological contacts. The computation of the horizontal gradients of the gravity field permits us to localize the limits of such blocks and then the fault locations. The horizontal derivative maxima of the Bouguer anomaly and its upward continuation at several heights show lineaments that could reflect the layout of faults and/or contacts and their dip directions. The application of this method to the Bouguer anomaly map of Morocco (with 19,571 points, using an average crustal density ρ = 2.67 g/cm³) allowed us to perform a multiscale analysis of the gravimetric lineaments of the country. The obtained structural map is consistent with several faults already identified in previous studies, and highlights five new major subsurface faults systems with location and dip: the Saghro fault system; Bou-Arfa Midelt fault system; Sidi Slimane Mezquitem fault; Ksar El Kebir–Chefchaouen fault and the Rifan West Mediterranean fault. In addition, this study suggests a new shape and localization for the Agadir-Oujda trans-Moroccan major fault with a NE-SW direction and 900 km length, subdividing Morocco into two main domains. The results of this study contribute to the improvement of the regional structural map of the north western part of Africa, which is situated within the convergence zone between Africa and Eurasia.     

Bayesian analysis on earthquake magnitude related to an active fault in Taiwan

It is understood that sample size could be an issue in earthquake statistical studies, causing the best estimate being too deterministic or less representative derived from limited statistics from observation. Like many Bayesian analyses and estimates, this study shows another novel application of the Bayesian approach to earthquake engineering, using prior data to help compensate the limited observation for the target problem to estimate the magnitude of the recurring Meishan earthquake in central Taiwan. With the Bayesian algorithms developed, the Bayesian analysis suggests that the next major event induced by the Meishan fault in central Taiwan should be in M_w 6.44±0.33, based on one magnitude observation of M_w 6.4 from the last event, along with the prior data including fault length of 14 km, rupture width of 15 km, rupture area of 216 km², average displacement of 0.7 m, slip rate of 6 mm/yr, and five earthquake empirical models.     

The April–June 2007 Trichonis Lake earthquake swarm (W. Greece): New implications toward the causative fault zone

On 10 April 2007, three moderate earthquakes with M_w = 4.9–5.1 occurred in the vicinity of Trichonis Lake (W. Greece). A local network composed of 12 three-component digital seismographs was installed in the epicentral area and recorded more than 1600 events. The double-difference algorithm HYPODD, incorporating both catalog and waveform cross-correlation differential travel-time data, was applied for the successful relocation of 1490 earthquakes. The latter led to the distinction of a main NW-SE trending and NE-dipping zone, as well as of three neighboring faults; a conjugate NW-SE striking and SW-dipping marginal fault mapped along the northeastern flanks of the lake; a E-W trending and south-dipping low-angle normal fault, possibly related to the major Agrinio Fault Zone (AFZ), parallel to the northern bank of the lake; a NE-SW striking and NW-dipping normal fault, likely related to a segment of the active Evinos fault, located south of the lake. Calculation of the Coulomb stress induced by the combination of the 1975 M_w = 6.0 event and the three largest events of 10 April 2007 on the inferred structures, reveals that most of the seismicity lies within the “stress-loaded” region, except for the westernmost activity, which probably belongs to the deep part of the AFZ. A total of 178 reliable focal mechanisms were determined by regional and local body-wave modeling (5 largest events) and P-wave first motion polarity data. The types of the obtained focal mechanisms are predominantly normal and strike-slip, however, numerous earthquakes were found to exhibit reverse faulting. Inversion of focal mechanism data showed that the prevailing principal horizontal component σ₃ is quite homogeneous throughout the activated area with a roughly NW-SE trend, parallel to the strike of the Hellenides. On the contrary, the compressional field σ₁ appears in two patterns: NE-SW trending onshore and NW-SE trending beneath the lake. This apparent rotation of σ₁ by 90° reveals a complex system enclosed by the suggested NW-SE trending antithetic faults in depths between 7 and 9 km. The calculated stress ratios beneath the lake imply that vertical forces are close to the overburden pressure. The overall inferred stress pattern is rather linked to topographic variations, locally imposing increase or decrease of the vertical forces. The presence of the water in the lake possibly plays an additional important role, penetrating through the bedrock, reducing the friction coefficient, while the pore pressure and, consequently, the effective stress increase. Thus, shearing along mature fractures is enhanced, likely yielding the observed diversity.     

Reflection seismic imaging of the end-glacial Pärvie Fault system,northern Sweden

Reflection seismic data were acquired along a c. 23 km long profile over the Pärvie Fault system with a nominal receiver and source spacing of 20 m. An hydraulic breaking hammer was used as a source, generating signals with a penetration depth of about 5–6 km. Steeply dipping reflections from the end-glacial faults are observed, as well as sub-horizontal reflections. The location and orientation of the reflections from the faults agree well with surface geological observations of fault geometries. Reflections from a potential fourth end-glacial fault is observed further to the east along the profile. The more sub-horizontal reflections may originate from gabbroic bodies within the granitic basement or from deeper lying greenstones. Our results indicate that the end-glacial faults dip at moderate to steep dips down to at least 2–3 km depth, and possibly continue at this dip to depths of 6 km. This result has significant implications for determining the state of stress required to activate the faults in the past and in the future.     

Active segment of the 12 November 2003 Mw 5.6 earthquake at Salsipuedes oceanic basin,Gulf of California. Mexico

We analyzed records of eight seismic stations of the autonomous broadband seismograph network of a joint project between Utrecht University (the Netherlands), California Institute of Technology, and Centro de Investigación Científica y de Estudios Superiores de Ensenada (CICESE). These stations recorded the Mw 5.6 earthquake that occurred on 12 November 2003 at Salsipuedes basin in the middle of the Gulf of California 2 km west of the island Angel de la Guarda. This event was located at 29.16º N and 113.37º W, 30 km northeast of Bahia de los Angeles. A foreshock and hundreds of aftershocks were recorded in the 48 hours after its origin time. With the location of 29 earthquakes we identified the active segment, perpendicular to the main transform fault NW–SE of Canal de Ballenas, representing the transtensional boundary between the Pacific and North American plates. The direction of the active fault described is consistent with the normal fault mechanism reported by the National Earthquake Information Center (strike=39º, dip=34º, slip=–44º).From the duration magnitude of 456 aftershocks, we calculated a b-value of 1.14±0.28; furthermore, we calculated a seismic moment of (3.5 ±3.3) X10¹⁷Nm, a source radius of 3.7 ± 2.63 km, and a static stress drop of 3.94 ± 1.15 MPa (39.4 ± 11.5 bar.).     

Temporal and spatial variations of seismicity scaling behavior in Southern México

R/S analysis is used in this work to investigate the fractal correlations in terms of the Hurst exponent for the 1998–2011 seismicity data in Southern Mexico. This region is the most seismically active area in Mexico, where epicenters for severe earthquakes (e.g., September 19, 1985, M_w = 8.1) causing extensive damage in highly populated areas have been located. By only considering the seismic events that meet the Gutenberg–Ritcher law completeness requirement (b = 0.97, M_GR = 3.6), we found time clustering for scales of about 100 and 135 events. In both cases, a cyclic behavior with dominant spectral components at about one cycle per year is revealed. It is argued that such a one-year cycle could be related to tidal effects in the Pacific coast. Interestingly, it is also found that high-magnitude events (M_w ≥ 6.0) are more likely to occur under increased interevent correlations with Hurst exponent values H > 0.65. This suggests that major earthquakes can occur when the tectonic stress accumulates in preferential directions. In contrast, the high-magnitude seismic risk is reduced when stresses are uniformly distributed in the tectonic shell. Such cointegration between correlations (i.e., Hurst exponent) and macroseismicity is confirmed for spatial variations of the Hurst exponent. In this way, we found that, using the Hurst exponent standpoint, the former presumed Michoacan and the Guerrero seismic gaps are the riskiest seismic zones. To test this empirical finding, two Southern Mexico local regions with large earthquakes were considered. These are the Atoyac de Alvarez, Guerrero (M_w = 6.3), and Union Hidalgo, Oaxaca (M_w = 6.6), events. In addition, we used the Loma Prieta, California, earthquake (October 17, 1989, M_w = 6.9) to show that the high-magnitude earthquakes in the San Andreas Fault region can also be linked to the increments of determinism (quantified in terms of the Hurst exponent) displayed by the stochastic dynamics of the interevent period time series. The results revealed that the analysis of seismic activity by means of R/S analysis could provide further insights in the advent of major earthquakes.     

Earthquake scaling,fault segmentation,and structural maturity

Slip and length measurements on earthquakes suggest large stress drop variability. We analyze an extended set of slip-length measurements for large earthquakes (M ≥ 6) to seek for the possible origin(s) of this apparent variability. We propose that such variability arises from earthquakes breaking a variable number of major fault segments. That number depends on the strength of the inter-segment zones, which itself depends on the structural maturity of the faults. We propose new D_max–L parameterizations based on that idea of multiple segment-ruptures. In such parameterizations, each broken segment roughly scales as a crack, while the total multi-segment rupture does not. Stress drop on individual segments is roughly constant, only varying between 3.5 to 9 MPa. The slight variation that is still observed depends on fault structural maturity; more mature faults have lower stress drops than immature ones. The new D_max–L functions that we propose reduce uncertainties with respect to available relationships. They thus provide a more solid basis to estimate seismic hazard by integrating fault properties revealed by geological studies.     

Time-Domain Moment Tensors for shallow (h ≤ 40 km) earthquakes in the broader Aegean Sea for the years 2006 and 2007: The database of the Aristotle University of Thessaloniki

We present a catalog of moment tensor (MT) solutions and moment magnitudes, M_w, for 119 shallow (h ≤ 40 km) earthquakes in Greece and its surrounding lands (34°N–42°N, 19°E–30°E) for the years 2006 and 2007, computed with the 1D Time-Domain Moment Tensor inversion method (TDMT_INV code of Dreger, 2003). Magnitudes range from 3.2 ≤ M_w ≤ 5.7. Green's functions (GF) have been pre-computed to build a library, for a number of velocity profiles applicable to the broader Aegean Sea region, to be used in the inversion of observed broad band waveforms (10–50 s). All MT solutions are the outcome of a long series of tests of different reported source locations and hypocenter depths. Quality factors have been assigned to each MT solution based on the number of stations used in the inversion and the goodness of fit between observed and synthetic waveforms. In general, the focal mechanisms are compatible with previous knowledge on the seismotectonics of the Aegean area. The new data provide evidence for strike-slip faulting along NW–SE trending structures at the lower part of Axios basin, close to the heavily industrialized, and presently subsiding, region of the city of Thessaloniki. Normal faulting along E–W trending planes is observed at the Strimon basin, and in Orfanou Gulf in northern Greece. A sequence of events in the east Aegean Sea close to the coastline with western Anatolia sheds light on an active structure bounding the north coastline of Psara–Chios Islands about 20–25 km in length exhibiting right lateral strike-slip faulting.     

A new fault rupture scenario for the 2003 Mw 6.6 Bam earthquake,SE Iran: Insights from the high-resolution QuickBird imagery and field observations

The December 26, 2003 M_w 6.6 Bam earthquake is one of the most disastrous earthquakes in Iran. QuickBird panchromatic and multispectral satellite imagery with 61 cm and 2.4 m ground resolution, respectively provide new insights into the surface rupturing process associated with this earthquake. The results indicate that this earthquake produced a 2–5 km-wide surface rupture zone with a complex geometric pattern. A 10-km-long surface rupture zone developed along the pre-existing Bam fault trace. Two additional surface rupture zones, each 2–5 km long, are oblique to the pre-existing Bam fault in angles of 20–35°. An analysis of geometric and geomorphic features also shows that movement on the Bam fault is mainly right-lateral motion with some compressional component. This interpretation is consistent with field investigations, analysis of aftershocks as well as teleseismic inversion. Therefore, we suggest that the 2003 Bam earthquake occurred on the Bam fault, and that the surface ruptures oblique to the Bam fault are caused by secondary faulting such as synthetic shears (Reidel shears). Our fault model for the Bam earthquake provides a new tectonic scenario for explaining complex surface deformations associated with the Bam earthquake.     

Deriving fault-slip histories to test for secular variation in slip,with examples from the Kunlun and Awatere faults

Although offset and age data from displaced landforms are essential for identifying earthquake clusters and thus testing whether faults slip at uniform or secularly varying rates, it is not clear how the uncertainties in such measurements should be propagated so as to yield a robust fault-slip history (i.e., record of fault displacement over time). Here we develop a Monte Carlo approach for estimating the distribution of geologically reasonable fault-slip histories that fit the offset and age data from a population of dated and displaced landforms. The model assumes that the landforms share common faulting histories, the offset and age constraints are correct, and the fault has not reversed shear sense. Analysis of the model results yields both a precise average slip rate, in the case where a linear fit is applied to the data, and a best-fit fault-slip history, in the case where the linear constraint is removed. The method can be used to test for secular variation in slip because the uncertainty on this best-fit history is quantified. By applying the method to previously published morphochronologic data from faulted late Quaternary terrace risers along the Kunlun fault in China and the Awatere fault in New Zealand, we have assessed the extent to which our modeled average slip rates match previously reported values and the data support previous interpretations of uniform slip rate. The Kunlun data set yields average slip rates of 8.7 + 3.6/?2.1 mm/yr and 5.1 + 1.6/?1.2 mm/yr (68.27% confidence), for the central and eastern reaches of the fault, respectively, both of which match previously published slip rates. Our analysis further indicates that these fault reaches have both slipped uniformly over the latest Quaternary. In contrast, analysis of data from the Saxton River site along the Awatere fault reveals a mid-Holocene deceleration in slip rate from 6.2 + 1.6/?1.4 mm/yr to 2.8 + 1.0/?0.6 mm/yr. This result contradicts previous interpretations of uniform slip along the Awatere fault. The Monte Carlo method we present here for quantifying fault-slip histories using the offset and age data from a population of faulted landforms provides an important tool for distinguishing temporally uniform from secularly varying fault slip.