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.     

Kinematics of the eastern part of the North Anatolian Fault Zone

The North Anatolian Fault Zone (NAFZ), which marks the boundary between Anatolia and the Eurasian plate, is one of the world's most seismically active structures. Although the eastern part of NAFZ has high seismic hazard, there is a lack of geodetic information about the present tectonics of this region. Even though many scientists would like to study this area, geographical and logistical problems make performing scientific research difficult. In order to investigate contemporary neotectonic deformation on the eastern NAFZ and in its neighborhood, a relatively dense Global Positioning System (GPS) monitoring network was established in 2003. Geodetic observations were performed in three GPS campaigns in an area of 350 km × 200 km with 12-month intervals. In addition, 14 new GPS stations were measured far from the deforming area. Since this region includes the intersection of the NAFZ and the East Anatolian Fault Zone (EAFZ), deformation is complex and estimating seismic hazard is difficult. One important segment is the Yedisu segment and it has not broken since the 1784 earthquake. After the 1992 Erzincan and 2003 Pulumur earthquakes, the Coulomb stress loading on the Yedisu segment of the NAFZ has increased significantly, emphasizing the need to monitor this region. We computed the horizontal velocity field with respect to Eurasia and strain rates field as well. GPS-derived velocities relative to Eurasia are in the range of 16–24 mm/year, which are consistent with the regional tectonics. The principal strain rates were derived from the velocity field. Results show that strain is accumulating between the NAFZ and EAFZ along small secondary fault branches such as the Ovacik Fault (OF).     

Modification of fault slip models of the Mw 9.0 Tohoku Earthquake by far field GPS observations

GPS data from Crustal Movement Observation Network of China (CMONOC) are used to derive far-field co-seismic displacements induced by the M_w 9.0 Tohoku Earthquake. Significant horizontal displacements about 30 mm, 10 mm, and 20 mm were caused by this large event in northeast China, north China, and on the Korean peninsula respectively. Vectors of relatively large horizontal displacements with dominant east components pointed to the epicenter of this earthquake. The east components show an exponential decay with the longitude, which is characteristic of the decay of the co-seismic horizontal displacements associated with earthquakes of thrust rupture. The exponential fit of the east components shows that the influence of the co-seismic displacements can be detected by GPS at a distance of about 3200 km from the epicenter of the earthquake. By considering the capability of the far field displacements for constraining the inversion of the fault slip model of the earthquake, we use spherically stratified Earth models to simulate the co-seismic displacements induced by this event. Using computations and comparisons, we discuss the effects of parameters of layered Earth models on the results of dislocation modeling. Comparisons of the modeled and observed displacements show that far field GPS observations are effective for constraining the fault slip model. The far field horizontal displacements observed by GPS are used to modify the slips and seismic moments of fault slip models. The result of this work is applicable as a reference for other researchers to study seismic source rupture and crustal deformation.     

GPS-derived velocity field of the Iznik-Mekece segment of the North Anatolian Fault Zone

Space-based tectonic studies on the western part of the North Anatolian Fault Zone (NAFZ) have been conducted over two decades. After the August 17, 1999, Izmit earthquake (M_w = 7.4), this region attracted greater scientific interest, and the collected data became more valuable. The Geodesy Department of the Kandilli Observatory and Earthquake Research Institute (KOERI) at Bogazici University established three micro-geodetic networks to the east of Akyazi, east of Iznik, and west of Lake Sapanca in the eastern part of the Marmara region; GPS data have been continually collected at these locations since 1994. The NAFZ branches out in the western part of the Marmara region and extends up to the Aegean Sea. Segments of the fault passing through the Marmara Sea are considered active, and this has increased concern regarding imminent earthquakes. Conventional geodetic measurements made between 1990 and 1994 are not sufficient for monitoring small movements. However, GPS has played a very important role in detecting such deformations in the area after 1994. The Iznik network, with 10 points, is bilaterally located on the Iznik-Mekece fault. Six years of GPS data for 2004–2010 collected for the monitoring of crustal deformation showed that the Iznik-Mekece fault segment moves westward at about 22 ± 1 mm/yr with respect to the Eurasia fixed reference frame. The GPS observations show that there is no strain accumulation in the area.     

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.     

Determining fault slip distribution of the Chi-Chi Taiwan earthquake with GPS and InSAR data using triangular dislocation elements

The reactivation of the Chelungpu fault triggered the 20 September 1999 Chi-Chi Taiwan earthquake (M_w = 7.6) which caused a 100-km long surface rupture that trends north–south. We reconstruct the fault geometry using 1068 planar triangular dislocation elements that approximate more realistically the curved three-dimensional fault surface. The fault slip distribution is then determined with the observed GPS coseismic displacements as well as interferometric synthetic aperture radar (InSAR) data. The results show that our smooth 3D fault slip model has improved the fit to the geodetic data by 44% compared with the previously published inversions. The slip distribution obtained both by inversion of GPS data only and by joint inversion of GPS and InSAR data indicates that notable slips occur on the sub-horizontal décollement at the depth of 6.1–8.9 km.     

Surface rupture zone of the 1303 Hongtong M=8 earthquake, Shanxi Province

Based on the latest displacement of Huoshan piedmont fault, Mianshan west-side fault and Taigu fault obtained from the beginning of 1990‘s up to the present, the characteristics of distribution and displacement of surface rupture zone of the 1303 Hongtong M = 8 earthquake, Shanxi Province are synthesized and discussed in the paper. If Taigu fault, Mianshan west-side fault and Huoshan piedmont fault were contemporarily active during the 1303 Hongtong M = 8 earthquake, the surface rupture zone would be 160 km long and could be divided into 3 segments, that is, the 50-km-long Huoshan piedmont fault segment, 35-km-long Mianshan west-side fault segment and 70-km-long Taigu fault segment, respectively. Among them, there exist 4 km and 8 km step regions. The surface rupture zone exhibits right-lateral features. The displacements of northern and central segments are respectively 6～7 m and the southern segment has the maximum displacement of 10 m. The single basin-boundary fault of Shanxi fault-depression system usually corresponds to M ≈ 7 earthquake, while this great earthquake (M = 8) broke through the obstacle between two basins. It shows that the surface rupture scale of great earthquake is changeable.     

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.     

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.     

Earthquakes produce carbon dioxide in crustal faults

Fourier transform infrared (FTIR) microanalysis of pseudotachylytes (i.e. friction-induced melts produced by seismic slip) from the Nojima fault (Japan) reveals that earthquakes almost instantaneously expel 99 wt.% of the wall rock CO₂ content. Carbon is exsolved because it is supersaturated in the friction melts. By extrapolation to a crustal-scale fault rupture, large events such as the M7.2 Kobe earthquake (1995) may yield a total production of 1.8 to 3.4 × 10³ tons CO₂ within a few seconds. This extraordinary release of CO₂ can cause a flash fluid pressure increase in the fault plane, and therefore enhance earthquake slip or trigger aftershocks; it may also explain the anomalous discharge of carbon monitored in nearby fault springs after large earthquakes. Because carbon saturation in silicate melts is pressure-dependent, FTIR can be used as a new tool to constrain the maximum depth of pseudotachylyte formation in exhumed faults.     

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.     

Sensitivity of ground-motion scenarios to earthquake source parameters in the Tehran metropolitan area,Iran

The northern Tehran fault (NTF) is a principal active fault of the Alborz mountain belt in the northern Iran. The fault is located north of the highly populated Metropolitan Area of Tehran. Historical records and paleoseismological studies have shown that the NTF poses a high seismic risk for the Tehran region and the surrounding cities (e.g. Karaj). A series of ground-motion simulations are carried out using a hybrid kinematic-stochastic model to calculate broadband (0.1–20 Hz) ground-motion time histories for deterministic earthquake scenarios (M7.2) on the NTF. We will describe the source characteristics of the target event to develop a list of scenario earthquakes that are probably similar to a large earthquake on the NTF. The effect of varying different rupture parameters such as rupture velocity and rise time on the resulting broadband strong motions has been investigated to evaluate the range of uncertainty in seismic scenarios. The most significant parameters in terms of ground-shaking level are the rise time and the value of the rupture velocity. For the worst-case scenario, the maximum expected horizontal acceleration, and velocity at rock sites in Tehran range between 128 and 1315 cm/s/s and 11–191 cm/s, respectively. For the lowest scenario, the corresponding values range between 102 and 776 cm/s/s and 12 to 81 cm/s. Nonlinear soil effects may change these results but are not accounted for in this study. The largest variability of ground motion is observed in neighborhood of asperity and also in the direction of rupture propagation. The calculated standard deviation of all ground-motion scenarios is less than 30% of the mean. The capability of the simulation method to synthesize expected ground motions and the appropriateness of the key parameters used in the simulations are confirmed by comparing the synthetic peak ground motions (PGA, PGV and response spectra) with empirical ground-motion prediction equations.     

Crustal deformation before the 2008 Wenchuan MS8.0 earthquake studied using GPS data

We used GPS velocities from approximately 700 stations in western China to study the crustal deformation before the Wenchuan M_S8.0 earthquake. The processing methods included analyses of the strain rate field, inversion of fault locking and the GPS velocity profiles. The GPS strain rate in the E-W direction in the Qinghai-Tibet block shows that extensional deformation was dominant in the western region of the block (west of 92.5° E), while compressive deformation predominated in the eastern region of the block (from 92.5° E to 100° E). On a regional scale, the hypocentral region of the Wenchuan earthquake was located at the edge of an intense compression deformation zone of about 1.9 × 10⁻⁸/a in an east-west direction. The characteristic deformation in the seismogenic fault was compressive with a dextral component. The compression deformation rate was greater in the fault's western region than in its eastern region, and the strain accumulation was very slow on the fault scale. The results of a fault locking inversion show that the locking fraction and slip deficit was greater in the middle-northern section of the seismogenic fault than in the southern section. The GPS velocity profile before the Wenchuan earthquake shows that the compression deformation was smaller than the dextral deformation, which is asymmetrical with respect to the distribution of co-seismic displacement. These deformation characteristics should provide some clues to the Wenchuan earthquake which occurred in the later period of the earthquake cycle.     

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.).     

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.     

Surface rupture zone of the 1303 Hongtong M=8 earthquake, Shanxi Province

Introduction The 1303 Shanxi Hongtong M=8 earthquake is the earliest M=8 event determined in histori-cal records in China and the largest recorded in Shanxi fault-depression system in history. Some researchers have discussed the tectonic environment of this earthquake (DENG, et al, 1973; DENG, 1984; DENG, XU, 1994, 1995; Seismo-geological Brigade, State Seismological Bureau, Depart-ment of Geology and Geography, Peking University, 1979; LIU, XIAO, 1982; ZHANG, JIA, 1986; SU, …     

GPS constraints on the deformation of Azerbaijan and surrounding regions

In this study, we present new GPS observations in Azerbaijan to provide an improved basis for determining the distribution of crustal deformation throughout the country and surrounding areas. The deformation field in the region has been analyzed with a dense GPS network configuration and a reliable quantification of the ongoing deformation was achieved. Results show that while contraction is dominant over the whole region, it is mostly concentrated on the middle and eastern parts of Caucasus Thrust Fault reaching up to 6.4 ± 0.2 mm/yr and Lesser Caucasus Fault does no accommodate more than 1–2 mm/yr of contraction. New network also clearly substantiates that the West Caspian Fault, which is a continuation of Caucasus Thrust Fault in the south, accommodates right-lateral slip rates of 7.1 ± 0.3 mm/yr in addition to 5.5 ± 0.3 mm/yr contraction rates.     

Near-realtime source analysis of the 20 March 2012 Ometepec-Pinotepa Nacional,Mexico earthquake

We apply a single-step, finite-fault analysis procedure to derive a coseismic slip model for the large M_W 7.4 Ometepec-Pinotepa Nacional, Mexico earthquake of 20 March 2012, using teleseismic P waveforms recorded by the Global Seismographic Network. The inversion is conducted in near-realtime using source parameters available from the USGS/NEIC and the Global Centroid Moment Tensor (gCMT) project. The fault orientation and slip angle are obtained from the gCMT mechanism assuming that the fault coincides with the shallow-dipping nodal plane. The fault dimensions and maximum rise time are based on the magnitude reported for the event. Teleseismic data from the USGS/NEIC Continuous Waveform Buffer database are used in the inversion with record start times set to the P-wave arrivals used to compute the earthquake hypocenter. The inversion is stabilized by requiring a smooth transition of slip across the fault while minimizing the seismic moment. These constraints are applied using a smoothing weight that is estimated from the inverse problem, allowing the recovery of the least-complicated rupture history in a single step. Inversion of the deconvolved, ground-displacement waveforms reveals a simple, circular rupture similar in extent to the source identified by the USGS/NEIC using body-and surface-wave data, indicating that the teleseismic P waves can provide a first-order source model for the event in near-realtime. Additional inversions conducted using velocity records identify a more-detailed rupture model characterized by an elliptical 2500 km² source region extending updip and downdip from the hypocenter. This elliptical source preserves the orientation and overall dimensions of a dual-source slip model obtained recently by other investigators using local strong motions and global seismic waveforms. The results indicate that velocity waveforms could provide additional details of the earthquake rupture in near-realtime, finite-fault inversions using teleseismic P waves.     

Refining the OSL age of the last earthquake on the Dheshir fault,Central Iran

In Central Iran there are several cities along the Dehshir fault, which have similar geological conditions to that of the city of Bam prior to the 2003 earthquake (Mw 6.5), during which more than 30,000 lives were lost. Optical stimulated luminescence (OSL) samples were collected from the Dehshir fault in order to place constraints on its seismic history. The single aliquot regenerative (SAR) dose measurement protocol on coarse grained quartz extracts was used for this study. This SAR protocol had to be optimized for the low OSL sensitivity by varying both the preheat temperatures and test doses used. Dose recovery tests showed that given laboratory dose could be successfully recovered. However, replicate palaeodose (D_e) data were scattered and consequently ages based on mean D_e's had large uncertainties. As this is thought to largely reflect poor bleaching conditions prior to sediment burial at the site, various statistical procedures were employed in conjunction with the stratigraphic knowledge of the site to try and extract more refined burial ages from the samples. From this the timing of the last earthquake was estimated around 2.0 ± 0.2 kyr. This refined age suggests that the earthquake catalogue of Iran is incomplete and more paleoseismological investigation is required to recognize and date the previous events of Dheshir fault.     

Constraining the long-term evolution of the slip rate for a major extensional fault system in the central Aegean,Greece, using thermochronology

The brittle/ductile transition is a major rheologic boundary in the crust yet little is known about how or if rates of tectonic processes are influenced by this boundary. In this study we examine the slip history of the large-scale Naxos/Paros extensional fault system (NPEFS), Cyclades, Greece, by comparing published slip rates for the ductile crust with new thermochronological constraints on slip rates in the brittle regime. Based on apatite and zircon fission-track (AFT and ZFT) and (U–Th)/He dating we observe variable slip rates across the brittle/ductile transition on Naxos. ZFT and AFT ages range from 11.8 ± 0.8 to 9.7 ± 0.8 Ma and 11.2 ± 1.6 to 8.2 ± 1.2 Ma and (U–Th)/He zircon and apatite ages are between 10.4 ± 0.4 to 9.2 ± 0.3 Ma and 10.7 ± 1.0 to 8.9 ± 0.6 Ma, respectively. On Paros, ZFT and AFT ages range from 13.1 ± 1.4 Ma to 11.1 ± 1.0 Ma and 12.7 ± 2.8 Ma to 10.5 ± 2.0 Ma while the (U–Th)/He zircon ages are slightly younger between 8.3 ± 0.4 Ma and 9.8 ± 0.3 Ma. All ages consistently decrease northwards in the direction of hanging wall transport. Most of our new thermochronological results and associated thermal modeling more strongly support the scenario of an identical fault dip and a constant or slightly accelerating slip rate of ∼ 6–8 km Myr^− 1 on the NPEFS across the brittle/ductile transition. Even the intrusion of a large granodiorite body into the narrowing fault zone at ∼ 12 Ma on Naxos does not seem to have affected the thermal structure of the area in a way that would significantly disturb the slip rate. The data also show that the NPEFS accomplished a minimum total offset of ∼ 50 km between ∼ 16 and 8 Ma.