Home About us Contact
|
Home About us Contact | ||
|
|
||
Faulting
Kinds of Faulting Selected AbstractsDiscrete-element modelling of detachment foldingBASIN RESEARCH, Issue 4 2005Stuart Hardy ABSTRACT A two-dimensional, discrete-element modelling technique is used to investigate the initiation and growth of detachment folds in sedimentary rocks above a weak décollement level. The model depicts the sedimentary rocks as an assemblage of spheres that obey Newton's equations of motion and that interact with elastic forces under the influence of gravity. Faulting or fracturing between neighbouring elements is represented by a transition from repulsive,attractive forces to solely repulsive forces. The sedimentary sequence is mechanically heterogeneous, consisting of intercalated layers of markedly different strengths and thicknesses. The interlayering of weak and strong layers within the sedimentary rocks promotes the localization of flexural flow deformation within the weak layers. Even with simple displacement boundary conditions, and straightforward interlayering of weak and strong layers, the structural geometries that develop are complex, with a combination of box, lift-off and disharmonic detachment fold styles forming above the décollement. In detail, it is found that the modelled folds grow by both limb rotation and limb lengthening. The combination of these two mechanisms results in uplift patterns above the folds that are difficult, or misleading, to interpret in terms of simple kinematic models. Comparison of modelling results with natural examples and with kinematic models highlights the complexities of structural interpretation in such settings. [source] Dextral-Slip Thrust Faulting and Seismic Events of the Ms 8.0 Wenchuan Earthquake, Longmenshan Mountains, Eastern Margin of the Tibetan PlateauACTA GEOLOGICA SINICA (ENGLISH EDITION), Issue 4 2009Zhenhan WU Abstract: Dextral-slip thrust movement of the Songpan-Garzę terrain over the Sichuan block caused the Ms 8.0 Wenchuan earthquake of May 12, 2008 and offset the Central Longmenshan Fault (CLF) along a distance of ,250 km. Displacement along the CLF changes from Yingxiu to Qingchuan. The total oblique slip of up to 7.6 m in Yingxiu near the epicenter of the earthquake, decreases northeastward to 5.3 m, 6.6 m, 4.4 m, 2.5 m and 1.1 m in Hongkou, Beichuan, Pingtong, Nanba and Qingchuan, respectively. This offset apparently occurred during a sequence of four reported seismic events, EQ1,EQ4, which were identified by seismic inversion of the source mechanism. These events occurred in rapid succession as the fault break propagated northeastward during the earthquake. Variations in the plunge of slickensides along the CLF appear to match these events. The Mw 7.5 EQ1 event occurred during the first 0,10 s along the Yingxiu-Hongkou section of the CLF and is characterized by 1.7 m vertical slip and vertical slickensides. The Mw 8.0 EQ2 event, which occurred during the next 10,42 s along the Yingxiu-Yanziyan section of the CLF, is marked by major dextral-slip with minor thrust and slickensides plunging 25°,35° southwestward. The Mw 7.5 EQ3 event occurred during the following 42,60 s and resulted in dextral-slip and slickensides plunging 10° southwestward in Beichuan and plunging 73° southwestward in Hongkou. The Mw 7.7 EQ4 event, which occurred during the final 60,95 s along the Beichuan-Qingchuan section of the CLF, is characterized by nearly equal values of dextral and vertical slips with slickensides plunging 45°,50° southwestward. These seismic events match and evidently controlled the concentrations of landslide dams caused by the Wenchuan earthquake in Longmenshan Mountains. [source] The Ohalo II prehistoric camp (19.5 Ky): New evidence for environmental and tectonic changes at the Sea of GalileeGEOARCHAEOLOGY: AN INTERNATIONAL JOURNAL, Issue 5 2002Shmuel Belitzky Combined archaeological data, shore surveys, and aerial photos of submerged sediments in the Sea of Galilee provide new insights into environmental and tectonic events, their dating, and their impact on the Ohalo II prehistoric camp (ca. 19,500 yr B.P.) and its surroundings. The Ohalo II waterlogged campsite contains excellently preserved brush hut remains and other in situ features, all embedded in late Pleistocene lacustrine strata. The findings indicate relatively short occupation of the site, not more than months or several years at a time. The high quality in situ preservation of delicate organic materials, as well as the short occupation period, suggests a quick and gentle burial by fine sediments. The evident fast submergence (water level rise of the Sea of Galilee) could have been the result of climatic fluctuations towards the end of the last glaciation and/or small-scale tectonic subsidence. The site is located on a tectonic block formed in the western fault belt of the Dead Sea Rift. We present new evidence of post-occupational folding of the late Pleistocene strata and recent tilting and faulting. A westward tectonic tilt may have caused the blockage of the old Jordan River outlet after A. D. 1106. Excellent preservation of the fault traces to the east of the site is attributed to the young age of the displacement on the fault. The last displacement apparently post-dates the blockage of the old Jordan River. © 2002 Wiley Periodicals, Inc. [source] Permo-Triassic development from Ireland to Norway: basin architecture and regional controlsGEOLOGICAL JOURNAL, Issue 6 2009tolfová Abstract Extensive occurrences of Permo-Triassic strata are preserved along the Northwest European Atlantic margin. Seismic reflection and well data are used to describe large-scale Permo-Triassic basin geometries along a swath of the continental shelf more than 2000,km long extending from the Irish to the mid-Norwegian sectors. Successions in the Celtic Sea, the flanks of the Irish Rockall Basin, basins west and north of Scotland, and the Trřndelag and Horda platforms west of Norway are described. The large-scale Permo-Triassic depositional geometries commonly represent erosional remnants of larger basins modified by later rifting episodes, uplift, inversion and continental breakup. However, the interpreted geometries reveal spatial and temporal differences in rifting style. The basins developed above a complex mosaic of petrologically heterogeneous crustal terranes with inherited crustal fabrics, which had a significant impact on the depositional basin geometries. Small Permian basins with growth faulting developed in the southern Celtic Sea region. Extensive, uniformly thick Triassic strata are characteristic of the wide rift basins in the southeastern Rockall Basin and northwest of the Solan Bank High. Thick, fault-controlled basins developed in the Horda and Trřndelag platform regions. The main controls on Permo-Triassic basin architecture are (a) crustal thickness and composition, which determined the development of narrow or wide rift basin geometries, (b) inherited Variscan, Caledonian and Precambrian basement structures and (c) pre-rift palaeotopography. Copyright © 2009 John Wiley & Sons, Ltd. [source] Tectono-sedimentary evolution of the northernmost margin of the NE German Basin between uppermost Carboniferous and Late Permian (Rotliegend)GEOLOGICAL JOURNAL, Issue 1 2001H. Rieke Abstract The tectono-sedimentary evolution of the Rotliegend deposits of the northernmost margin of NE German Basin (NEGB) has been analysed on the basis of detailed sedimentary logs of 300,m of core material together with the re-evaluation of 600,km of seismic lines. Three distinct phases were recognized. During the initial Phase I, basin geometry was largely controlled by normal faulting related to deep-seated ductile shearing leading to a strong asymmetric shape, with a steep fault-controlled eastern margin and a gently, dipping western margin. The results of forward modelling along a cross-section fit the basin geometry in width and depth and reveal a footwall uplift of c. 1000,m. Adjacent to the steep faults, local sedimentation of Lithofacies Type I was confined to non-cohesive debris flow-dominated alluvial fans, whereas the gently dipping western margin was dominated by alluvial-cone sedimentation. During the post-extensional period (Phase II), cooling of the lithosphere generated additional accommodation space. The sediments of Lithofacies Type II, comprising mainly clast-supported conglomerates, are interpreted as braided ephemeral stream flow-surge deposits. Tectonic quiescence and an increase in flood events resulting from wetter climate led to progradation of this facies over the entire region. At the end of this period, the accommodation space was almost completely filled resulting in a level topography. Phase III was controlled by the thermal-induced subsidence of the southerly located NEGB in post-Illawarra times. The formerly isolated region tilted towards the SW, thus forming the northern margin of the NEGB during uppermost Havel and Elbe Subgroup times. The sediments of Lithofacies Type III were divided into a marginal sandstone-dominated environment and a finer-grained facies towards the SW. The former consists of poorly-sorted coarse-grained sandstones of a proximal and medial ephemeral stream floodplain facies. The latter comprise mud flat fines and fine-grained distal ephemeral stream deposits. The end of the tectono-sedimentary evolution is marked by the basinwide Zechstein transgression. Copyright © 2001 John Wiley & Sons, Ltd. [source] Seismic constraints on the three-dimensional geometry of low-angle intracrustal reflectors in the Southern Iberia Abyssal PlainGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 2 2008S. M. Dean SUMMARY Several lines of evidence suggest that simple shear rifting of the continental crust, in the form of low-angle detachment faulting, occurred during the final stages of continental breakup between West Iberia and the Grand Banks. The primary evidence for such faulting is the occurrence of low-angle, high amplitude reflectors within the basement adjacent to the ocean,continent transition zone. Here we present a series of intersecting, depth migrated seismic reflection profiles that image one such reflector, the H-reflector, located on the southern edge of Galicia Bank. ,H' lies beneath several boreholes drilled during ODP Legs 149 and 173, in a region where the oceanward extent of extended continental crust steps at least 150 km westward from its location in the southern Iberia Abyssal Plain to its location off the relatively shallow Galicia Bank. In our profiles ,H' appears to define a surface that extends over a region of at least 200 km2 and that dips down ,19° to the north, towards Galicia Bank. The profiles show that a close affinity exists between ,H' and the most seaward continental crust. Based on geophysical data and ODP drilling results, we infer that the basement above ,H' is composed of continental crust deformed by extensional faults into a series of wedge-shaped blocks and thin slivers. These basement wedges have a complex 3-D geometry. ,H' rises to the basement surface on a number of the seismic profiles and appears to define locally the oceanward extent of continental fault blocks. [source] A structural model for the seismicity of the Arudy (1980) epicentral area (Western Pyrenees, France)GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2007Noalwenn Dubos-Sallée SUMMARY The Western Pyrenees presents a diffuse and moderate (M, 5.7) instrumental seismicity. It nevertheless historically suffered from strong earthquakes (I = IX MSK). The seismic sources of these events are not yet clearly identified. We focus on the Arudy (1980) epicentral area (M= 5.1) and propose here the reactivation of early Cretaceous normal faults of the Iberian margin as a potential source. The late Cretaceous inversion of this basin, first in a left-lateral strike-slip mode and then in a more frontal convergence, resulted in a pop-up geometry. This flower structure attests of the presence of a deep crustal discontinuity. The present-day geodynamic arrangement suggests that this accident is reactivated in a right lateral mode. This reactivation leads to a strain partitioning between the deep discontinuity that accommodates the lateral component of the motion and shallow thrusts, rooted on this discontinuity. These thrusts accommodate the shortening component of the strain. The distribution of the instrumental seismicity fits well the structural model of the Arudy basin. Whatever the compressive regional context, the structural behaviour of the system explains too the extensive stress tensor determined for the Arudy crisis if we interpret it in terms of strain ellipsoid. Indeed numerical modelling has shown that this concomitant activity of strike-slip and thrust faulting results in an extensive component that can rise 50 per cent of the finite strain. We identify too a 25,30 km long potential seismic source for the Arudy area. The size of the structure and its potential reactivation in a strike-slip mode suggest that a maximum earthquake magnitude of ,6.5 could be expected. The extrapolation of this model at the scale of the Western Pyrenees allows to propose other potential sources for major regional historical earthquakes. [source] 2001 August earthquake swarm at Shadwan Island, Gulf of Suez, EgyptGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2006Ahmed Badawy SUMMARY The earthquake swarm that struck Shadwan Island at the entrance of the Gulf of Suez in 2001 August included 408 events. Almost all of these events (94 per cent) were microearthquakes and only 6 per cent had small measurable magnitudes (5.0 > ML, 3.0). Most of the earthquakes were weak and followed each other so closely in time that they could not be identified at more distant stations. The fault plane solutions of the strongest events of the swarm show almost identical focal mechanisms, predominately normal faulting with a significant sinistral strike-slip component for nodal planes trending NW,SE. A comparison with the mechanisms of the 1969 and 1972 events which took place 20 km north of the swarm epicentral region shows similarities in faulting type and orientation of nodal planes. The azimuths of T -axes determined from focal mechanisms in this study are oriented in the NNE,SSW direction. This direction is consistent with the present-day stress field derived from borehole breakouts in the southern Gulf of Suez and the last phase of stress field changes in the Late Pleistocene, as well as with recent GPS results. The source parameters of the largest (ML, 3.0) events of the 2001 August Shadwan swarm have been estimated from the P -wave spectra of the Egyptian National Seismograph Network (ENSN). Averaging of the values obtained at different stations shows relatively similar source parameters, including a fault length of 0.65 ,L, 2 km, a seismic moment of 7.1 × 1012,Mo, 3.0 × 1014 N m and a stress drop of 0.4 ,,,, 10 bar. [source] Crosswell seismic waveguide phenomenology of reservoir sands & shales at offsets >600 m, Liaohe Oil Field, NE ChinaGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2005P. C. Leary SUMMARY Crosswell seismic data recorded at 620,650 m offsets in an oil-bearing sand/shale reservoir formation at the Liaohe Oil Field, northeast China, provide robust evidence for waveguide action by low-velocity reservoir layers. Crosswell-section velocity models derived from survey-well sonic logs and further constrained by observed waveguide seismic wavegroup amplitudes and phases yield plausible evidence for interwell reservoir,sand continuity and discontinuity. A pair of back-to-back Liaohe crosswell vector-seismic surveys were conducted using a source well between two sensor wells at 650 and 620 m offsets along a 200-m-thick reservoir formation dipping 7° down-to-east between depths of 2.5 and 3 km. A downhole orbital vibrator generated seismic correlation wavelets with frequency range 50,350 Hz and signal/noise ratio up to 5:1 over local downhole ambient noise. The sensor wells were instrumented with a mobile 12- to 16-level string of clamped vector-motion sensor modules at 5 m intervals. Using 5 m source depth increments, crosswell Surveys 1 and 2 cover source/sensor well intervals above and through the reservoir of, respectively, 600 m/600 m (13 000 vector traces in 9 common sensor fans) and 300 m/560 m (7000 vector traces in 7 common sensor fans). Survey 1 common sensor gathers show clear, consistent high-amplitude 20 ms waveletgroup lags behind the first-arrival traveltime envelope. Such arrivals are diagnostic of seismic low-velocity waveguides connecting the source and sensor wells. Observed Survey 1 retarded wavegroup depths tally with source and sensor depths in low-velocity layers identified in sonic well logs. Finite-difference acoustic model wavefields computed for waveguide acoustic layers constrained by well-log sonic velocity data match the observed waveguide traveltime and amplitude systematics. Model waveforms duplicate the observed m-scale and ms-scale sensitivity of waveguide spatio-temporal energy localization. Survey 2 crosswell data, in contrast, provide no comparable evidence for waveguide action despite a sensor-well sonic well log similar to that of Survey 1. Instead, acoustic wavefield modelling of Survey 2 data clearly favours an interpreted waveguide model with 10° downdip interrupted by a 75,100 m throw down-fault near the sensor well. The absence of clear waveguide arrivals is adequately explained by dispersal of waveguide energy at the fault discontinuity. Auxiliary well sonic velocity and lithologic logs confirm the model-implied 75,100 m of down-throw faulting near the sensor well. [source] Contemporary kinematics of the southern Aegean and the Mediterranean RidgeGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2004Corné Kreemer SUMMARY This study focuses on the kinematics of the southern Aegean and the Mediterranean Ridge (MR). A quantification of the deformation of the MR is essential for both evaluating physical models of accretionary wedges in general and for obtaining a self-consistent model of the surface deformation over the entire Nubia,Eurasia (NU,EU) plate boundary zone in the eastern Mediterranean. Previous kinematic studies have not properly considered the deformation field south of the Hellenic arc. Although this study focuses on the deformation field of the MR, we also discuss the kinematics of the southern Aegean, because the geometry and movement of the Hellenic arc determine to a large extent the kinematic boundary conditions for kinematic studies of the MR. We calculate a continuous velocity and strain rate field by interpolating model velocities that are fitted in a least-squares sense to published Global Positioning System (GPS) velocities. In the interpolation, we use information from a detailed data set of onshore and offshore active faulting to place constraints on the expected style and direction of the model strain rate field. In addition, we use the orientations of tracks left by seamounts travelling into the wedge to further constrain the offshore deformation pattern. Our model results highlight the presence of active shear partitioning within the Mediterranean ridge. High compressional strain rates between the ridge crest and the deformation front accommodate approximately 60,70 per cent of the total motion over the wedge, and the outward growth rate of the frontal thrust is , 4 mm yr,1. Strain partitioning within the wedge leads to 19,23 mm yr,1 of dextral motion at the wedge,backstop contact of the western MR, whereas the Pliny and Strabo trenches in the eastern MR accommodate 21,23 mm yr,1 of sinistral motion. The backstop of the western MR is kinematically part of the southern Aegean, which moves as a single block [the Aegean block (AE)] at 33,34 mm yr,1 in the direction of S24°W ± 1° towards stable Nubia (NU). Our model confirms that there is a clear divergence between the western and eastern Hellenic arc and we argue for a causal relation between the outward motion of the arc and the gradient in the regional geoid anomaly. Our results suggest that a significant driving source of the surface velocity field lies south of the Hellenic arc and only for the southeastern Aegean could there be some effect as a result of gravitational collapse associated with density differences within the overriding plate. [source] Analysis of the 2002 May earthquake sequence in the central Pyrenees, consequences for the evaluation of the seismic risk at Lourdes, FranceGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2004Noalwenn Dubos SUMMARY Three earthquakes of magnitudes 4.6, 4.3 and 3.7 occurred in 2002 May at two locations 20 km from the pilgrimage city of Lourdes in the French Pyrenees. They were well recorded by the permanent Pyrenean seismic networks, by a temporary local network, as well as by accelerometric stations. In order to understand their tectonic contexts, and to come to a better evaluation of the seismic risk at Lourdes, a detailed analysis of these events is performed. The first two events are located south of Lourdes in an area where only a few earthquakes have occurred up to now. Their focal solutions derived from first-motion polarities indicate reverse faulting, with a N110°E strike consistent with the geological structures. 10 aftershocks were recorded and relocated with respect to the main events, benefiting from the waveform similarity of the various events. This analysis reveals that the two main events concern probably the same fault, the second rupture being in the prolongation of the first one, whereas the other small aftershocks are located on fault segments in the vicinity of the hypocentre of the second event. The third large event, located to the SE of Lourdes, involves a normal mechanism with a N120°E plane parallel to the main geological structures. It occurred in a region of intense activity, including in particular an event of maximum macroseismic intensity IX in 1660. The first two events are at the boundary of a large quiet zone. In order to understand the related structural context, a new crustal tomographic model has been computed. It reveals that this quiet zone coincides with a block of high P -velocity. In contrast, the seismicity appears to be stronger at the northern and eastern boundaries of this block. The accelerometric data of the three main events recorded at Lourdes have been used to estimate the maximum peak ground accelerations in this city if a large event occurred, similar to those which damaged the city in the seventeenth and eighteenth centuries. Horizontal accelerations of 0.25 ± 0.07 g are predicted in the frequency domain 1,5 Hz at the location of the Sanctuary for a magnitude 6 event occurring 10 km away from the city. Taking into account the error bars, these values could in some cases exceed those specified by the building codes in this region. [source] The source process of the 2001 July 26 Skyros Island (Greece) earthquakeGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2004Zafeiria Roumelioti SUMMARY The spatial and temporal distribution of slip during the 2001 July 26 Skyros (Greece) earthquake Moment magnitude (M 6.5) is investigated using broadband data recorded at regional distances. The applied method involves estimation of the source time functions of the examined event through an empirical Green's function approach and inversion of their shapes to estimate kinematic source parameters. Our test inversions to statistically identify the fault plane, together with the distribution of aftershocks clearly indicate sinistral strike-slip faulting. In view of the fact that the Skyros epicentre lies near the western termination of the dextral strike-slip North Anatolian Fault (NAF) into the Aegean Sea, this sinistral strike-slip motion, for the first time instrumentally identified, has great tectonic significance. The best values searched through the inversion are 0.7 s for the rise time, and 2.4 km s,1 for the rupture velocity. Most of the slip appears to be concentrated in a relatively small area around the hypocentre, while a smaller slip patch was found at relatively large depth (18,24 km). At least two of the large aftershocks following the main event also occurred at the deeper part of the fault. Smaller amounts of slip are distributed in a wider area with dimensions similar to those inferred from the aftershock distribution studies and the empirical relations applicable to Greece. [source] Fault interactions and subduction tectonics: a re-examination of the Weber, New Zealand, earthquake sequence of 1990GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2003Russell Robinson SUMMARY Two moderate magnitude (Mw= 6.2 and 6.4) earthquakes in the Hikurangi subduction margin, North Island, New Zealand, occurred 3 months apart in 1990. The epicentres are nearly coincident, but the first (Weber 1, primarily normal faulting) occurred within the subducting Pacific Plate (depth about 28 km) and the second (Weber 2, a mix of thrusting and right-lateral motion) occurred within the overlying Australian Plate (depth about 13 km), the plate interface in between. The plate interface is interpreted to be locked trenchward (SE) from about the position of these events, with a transition to aseismic slip further down-dip to the NW. Several stress interaction questions are examined. First, to see whether Weber 1 triggered Weber 2, a range of possible mainshock parameters are used to calculate induced changes in the static Coulomb failure stress (,CFS). In most cases the results are consistent with triggering. Secondly, previous work showed that the rate of aftershock occurrence for Weber 1 decreased markedly about 35 days before Weber 2, recovering afterwards. To see whether aseismic pre-slip on the Weber 2 fault, as predicted by rate and state friction, could be the cause of the decrease, the same fault parameters have been used in reverse. The results are ambiguous, some fault parameters giving results consistent with the hypothesis and others not. The amount of pre-slip required for significant inhibition, however, is about equal to that in the mainshock and distributed over the entire fault plane. Thirdly, observations of episodic, aseismic slip down-dip from locked sections of other plate interfaces are becoming more common. Could such slip have triggered both Weber events? The induced changes in CFS for such slip are uniformly positive on the Weber 1 fault plane, and mostly positive on the Weber 2 fault plane, so the answer is yes. Although there is no independent evidence for aseismic slip prior to the Weber sequence, this case shows that such slip may trigger events on other nearby faults, besides loading the locked section of the plate interface. Static stress triggering considerations are thus likely to be important in subduction environments. [source] The 1933 Ms= 7.3 Baffin Bay earthquake: strike-slip faulting along the northeastern Canadian passive marginGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2002Allison L. Bent Summary The 1933 November 20 (Ms= 7.3) Baffin Bay earthquake is one of the largest instrumentally recorded passive margin earthquakes. Analysis of seismograms of this earthquake shows strong evidence for strike-slip faulting, which contrasts with the generally accepted belief that Baffin Bay is dominated by thrust faulting. The best-fitting solution consists of a large strike-slip subevent (strike 172°, dip 82°, rake 6°) followed by two smaller oblique-thrust subevents (strike 190°, dip 30°, rake 62°). All subevents occur at a depth of about 10 km. An instrumental moment magnitude of 7.4 was determined. Preliminary analysis of subsequent large (magnitude , 6.0) earthquakes in Baffin Bay finds additional evidence for strike-slip faulting in the region. The results for Baffin Bay, together with those for other passive margin earthquakes, suggest strike-slip faulting may be more prevalent in these regions than was previously believed. [source] On strike-slip faulting in layered mediaGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2002Maurizio Bonafede Summary We study the effects of structural inhomogeneities on the stress and displacement fields induced by strike-slip faults in layered media. An elastic medium is considered, made up of an upper layer bounded by a free surface and welded to a lower half-space characterized by different elastic parameters. Shear cracks with assigned stress drop are employed as mathematical models of strike-slip faults, which are assumed to be vertical and planar. If the crack is entirely embedded within the lower medium (case A), a Cauchy-kernel integral equation is obtained, which is solved by employing an expansion of the dislocation density in Chebyshev polynomials. If the crack is within the lower medium but it terminates at the interface (case B), a generalized Cauchy singularity appears in the integral kernel. This singularity affects the singular behaviour of the dislocation density at the crack tip touching the interface. Finally, the case of a crack crossing the interface is considered (case C). The crack is split into two interacting sections, each placed in a homogeneous medium and both open at the interface. Two coupled generalized Cauchy equations are obtained and solved for the dislocation density distribution of each crack section. An asymptotic study near the intersection between the crack and the interface shows that the dislocation densities for each crack section are bounded at the interface, where a jump discontinuity is present. As a corollary, the stress drop must be discontinuous at the interface, with a jump proportional to the rigidity contrast between the adjoining media. This finding is shown to have important implications for the development of geometrical complexities within transform fault zones: planar strike-slip faults cutting across layer discontinuities with arbitrary stress drop values are shown to be admissible only if the interface between different layers becomes unwelded during the earthquake at the crack/interface junction. Planar strike-slip faulting may take place only in mature transform zones, where a repetitive earthquake cycle has already developed, if the rheology is perfectly elastic. Otherwise, the fault cannot be planar: we infer that strike-slip faulting at depth is plausibly accompanied by en-echelon surface breaks in a shallow sedimentary layer (where the stress drop is lower than prescribed by the discontinuity condition), while ductile deformation (or steady sliding) at depth may be accommodated by multiple fault branching or by antithetic faulting in the upper brittle layer (endowed with lower rigidity but higher stress). [source] Seismic activity triggered by the 1999 Izmit earthquake and its implications for the assessment of future seismic riskGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2001Ali Pinar Summary A serious question has remained as to the location of the western end of the main rupture zone associated with the 1999 Izmit, Turkey, earthquake. A clear answer to this question is extremely important for the assessment of future seismic risk in the eastern Marmara Sea region, Turkey. In this paper we show an effective approach to answering this important question, unifying different kinds of information such as seismic activity, focal mechanism solutions and stress changes caused by the main shock into a clear image. We first point out that the major moment release is 1.6 × 1020 N m and covered the area between 29.7°E and 30.5°E and we then claim that the enhanced seismic activity after the main shock in the eastern Marmara Sea region should be regarded as activity triggered by the increase of stress, rather than as aftershock activity along the ruptured zone. We propose three fault segments with an average stress increase on each in the western extension of the main-shock rupture zone as potential sites for future large earthquakes, namely (i) the 50 km long Yalova,Hersek segment (0.45 MPa), (ii) the NW,SE-trending right-lateral strike-slip fault known as the Princes Islands segment (0.18 MPa), and (iii) the Çinarcik,Yalova segment (0.09 MPa) characterized by normal faulting, which was subject to rupture in 1963. [source] Analyses of the stress field in southeastern France from earthquake focal mechanismsGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 2 2001Emmanuel Baroux Summary Owing to the apparent deformation field heterogeneity, the stress regimes around the Provence block, from the fronts of the Massif Central and Alpine range up to the Ligurian Sea, have not been well defined. To improve the understanding of the SE France stress field, we determine new earthquake focal mechanisms and compute the present-day stress states by inversion of the 89 available focal mechanisms around the Provence domain, including 17 new ones calculated in the current study. This study provides evidence of six distinct deformation domains around the Provence block, with different tectonic regimes. On a regional scale, we identify three zones characterized by significantly different stress regimes: a western one affected by an extensional stress (normal faulting) regime; a southeastern one characterized by a compressional stress (reverse to strike-slip faulting) regime with NNW- to WNW-trending ,1; and a northeastern one, namely the Digne nappe front, marked by a NE-trending compression. Note that the Digne nappe back domain is controlled by an extensional regime that is deforming the western Alpine core. This extensional regime could be a response to buoyancy forces related to the Alpine high topography. The stress regimes in the southeast of the Argentera Massif and around the Durance fault are consistent with a coherent NNW-trending ,1, implying a left-lateral component of the active reverse oblique slip of the Moyenne Durance Fault. In the Rhone Valley, an E-trending extension characterizes the tectonic regime, implying a normal component of the present-day N,^mes fault displacement. This study provides evidence for short-scale variation of the stress states, which arises from abrupt changes in the boundary force influences on upper crustal fragments (blocks). These spatial stress changes around the Provence block result from the coeval influence of forces applied at both its extremities, namely in the northeast the Alpine front push, and in the southeast the northward African plate drift. In addition to these boundary forces, the mantle plume under the Massif Central influences the western block boundary. [source] Hydraulic pathways in the crystalline rock of the KTBGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2000Günter Zimmermann Fracture systems and fluid pathways must be analysed in order to understand the dynamical processes in the upper crust. Various deterministic as well as stochastic fracture networks in the depth section of the Franconian Lineament (6900 to 7140 m), which appears as a brittle ductile shear zone and prominent seismic reflector, were modelled to simulate the hydraulic situation at the two boreholes of the Continental Deep Drilling Program (KTB). They led to estimations of the hydraulic permeability in crystalline rock. The geometrical parameters of the fractures, such as fracture locations and orientations, were determined from structural borehole measurements, which create an image of the borehole wall. The selection of potentially open fractures was decided according to the stress field. Only fractures with the dip direction (azimuth) of the fracture plane perpendicular to the maximum horizontal stress field were assumed to be open. The motivation for this assumption is the fact that the maximum horizontal stress is higher than the vertical stress from the formation, indicating that the state of stress is a strike-slip faulting. Therefore, the probability of open fractures due to this particular stress field at the KTB sites is enhanced. Length scales for fracture apertures and extensions were stochastically varied and calibrated by hydraulic experiments. The mean fracture aperture was estimated to be 25 ,m, assuming an exponential distribution, with corresponding permeability in the range of 10,16 m2. Similar results were also obtained for log-normal and normal distributions, with a variation of permeability of the order of a factor of 2. The influence of the fracture length on permeability of the stochastic networks was also studied. Decreasing the fracture length beyond a specific threshold of 10 m led to networks with vanishing connectivity and hence vanishing permeability. Therefore, we assume a mean fracture length exceeding the threshold of 10 m as a necessary assumption for a macroscopic hydraulically active fracture system at the KTB site. The calculated porosity due to the fracture network is of the order of 10,3 per cent, which at first sight contradicts the estimated matrix porosity of 1 to 2 per cent from borehole measurements and core measurements. It can be concluded from these results, however, that if the fluid transport is due to a macroscopic fracture system, only very low porosity is needed for hydraulic flow with permeabilities up to several 10,16 m2, and hence the contribution of matrix porosity to the hydraulic transport is of a subordinate nature. [source] Microseismicity and faulting geometry in the Gulf of Corinth (Greece)GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 2 2000Denis Hatzfeld During the summer of 1993, a network of seismological stations was installed over a period of 7 weeks around the eastern Gulf of Corinth where a sequence of strong earthquakes occurred during 1981. Seismicity lies between the Alepohori fault dipping north and the Kaparelli fault dipping south and is related to both of these antithetic faults. Focal mechanisms show normal faulting with the active fault plane dipping at about 45° for both faults. The aftershocks of the 1981 earthquake sequence recorded by King et al. (1985) were processed again and show similar results. In contrast, the observations collected near the western end of the Gulf of Corinth during an experiment conducted in 1991 (Rigo et al. 1996), and during the aftershock studies of the 1992 Galaxidi and the 1995 Aigion earthquakes (Hatzfeld et al. 1996; Bernard et al. 1997) show seismicity dipping at a very low angle (about 15°) northwards and normal faulting mechanisms with the active fault plane dipping northwards at about 30°. We suggest that the 8,12 km deep seismicity in the west is probably related to the seismic,aseismic transition and not to a possible almost horizontal active fault dipping north as previously proposed. The difference in the seismicity and focal mechanisms between east and west of the Gulf could be related to the difference in the recent extension rate between the western Gulf of Corinth and the eastern Gulf of Corinth, which rotated the faults dipping originally at 45° (as in the east of the Gulf) to 30° (as in the west of the Gulf). [source] An extended finite element framework for slow-rate frictional faulting with bulk plasticity and variable frictionINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 13 2009Fushen Liu Abstract We present an extended finite element (FE) approach for the simulation of slow-rate frictional faulting in geologic media incorporating bulk plasticity and variable friction. The method allows the fault to pass through the interior of FEs without remeshing. The extended FE algorithm for frictional faulting, advocated in two recent articles, emanates from a variational equation formulated in terms of the relative displacement on the fault. In the present paper we consider the combined effects of bulk plasticity and variable friction in a two-dimensional plane strain setting. Bulk plasticity is localized to the fault tip and could potentially be used as a predictor for the initiation and propagation of new faults. We utilize a variable velocity- and state-dependent friction, known as the Dieterich,Ruina or ,slowness' law, formulated in a slip-weakening format. The slip-weakening/variable friction model is then time-integrated according to the generalized trapezoidal rule. We present numerical examples demonstrating the convergence properties of a global Newton-based iterative scheme, as well as illustrate some interesting properties of the variable friction model. Copyright © 2009 John Wiley & Sons, Ltd. [source] Inelastic constitutive properties and shear localization in Tennessee marbleINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 2 2001D. J. Holcomb Abstract The inelastic response of Tennessee marble is modelled by an elastic plastic constitutive relation that includes pressure dependence of yield, strain-softening and inelastic volume strain (dilatancy). Data from 12 axisymmetric compression tests at confining pressures from 0 to 100 MPa are used to determine the dependence of the yield function and plastic potential, which are different, on the first and second stress invariants and the accumulated inelastic shear strain. Because the data requires that the strain at peak stress depends on the mean stress, the locus of peak stresses is neither a yield surface nor a failure envelope, as is often assumed. Based on the constitutive model and Rudnicki and Rice criterion, localization is not predicted to occur in axisymmetric compression although faulting is observed in the tests. The discrepancy is likely due to the overly stiff response of a smooth yield surface model to abrupt changes in the pattern of straining. The constitutive model determined from the axisymmetric compression data describes well the variation of the in-plane stress observed in a plane strain experiment. The out-of-plane stress is not modelled well, apparently because the inelastic normal strain in this direction is overpredicted. In plane strain, localization is predicted to occur close to peak stress, in good agreement with the experiment. Observation of localization on the rising portion of the stress,strain curve in plane strain does not, however, indicate prepeak localization. Because of the rapid increase of mean stress in plane strain, the stress,strain curve can be rising while the shear stress versus shear strain curve at constant mean stress is falling (negative hardening modulus). Copyright © 2001 John Wiley & Sons, Ltd. [source] Anisotropy of magnetic susceptibility and petrofabric studies in the Garhwal synform, Outer Lesser Himalaya: Evidence of pop-up klippenISLAND ARC, Issue 3 2009Upasana Devrani Abstract Geological field, petrographic, and anisotropy of magnetic susceptibility studies help in understanding the evolutionary history of the Garhwal synform that lies in the western Outer Lesser Himalaya. Orientations of the magnetic susceptibility axes reveal large variations at short distances as a result of superimposed deformation, and predominant stress conditions favorable for normal faulting. Rocks forming the outer limbs of the Garhwal Synform are metamorphosed up to the lower greenschist facies. The metamorphic grade increases to chlorite zone in the inner limb and the core is characterized by chlorite,biotite to garnet zones. The different grades of metamorphism are separated by thrusts and the structure is described as a pop-up klippen. [source] Structure of Sumatra and its implications for the tectonic assembly of Southeast Asia and the destruction of PaleotethysISLAND ARC, Issue 1 2009Anthony J. Barber Abstract It is now generally accepted that Southeast Asia is composed of continental blocks which separated from Gondwana with the formation of oceanic crust during the Paleozoic, and were accreted to Asia in the Late Paleozoic or Early Mesozoic, with the subduction of the intervening oceanic crust. From east to west the Malay peninsula and Sumatra are composed of three continental blocks: East Malaya with a Cathaysian Permian flora and fauna; Sibumasu, including the western part of the Malay peninsula and East Sumatra, with Late Carboniferous,Early Permian ,pebbly mudstones' interpreted as glaciogenic diamictites; and West Sumatra, again with Cathaysian fauna and flora. A further unit, the Woyla nappe, is interpreted as an intraoceanic arc thrust over the West Sumatra block in the mid Cretaceous. There are varied opinions concerning the age of collision of Sibumasu with East Malaya and the destruction of Paleotethys. In Thailand, radiolarites have been used as evidence that Paleotethys survived until after the Middle Triassic. In the Malay peninsula, structural evidence and the ages of granitic intrusions are used to support a Middle Permian to Early Triassic age for the destruction of Paleotethys. It is suggested that the West Sumatra block was derived from Cathaysia and emplaced against the western margin of Sibumasu by dextral transcurrent faulting along a zone of high deformation, the Medial Sumatra Tectonic Zone. These structural units can be traced northwards in Southeast Asia. The East Malaya block is considered to be part of the Indochina block, Sibumasu can be traced through Thailand into southern China, the Medial Sumatra Tectonic Zone is correlated with the Mogok Belt of Myanmar, the West Burma block is the extension of the West Sumatra block, from which it was separated by the formation of the Andaman Sea in the Miocene, and the Woyla nappe is correlated with the Mawgyi nappe of Myanmar. [source] Transpressional tectonics of the Mineoka Ophiolite Belt in a trench,trench,trench-type triple junction, Boso Peninsula, JapanISLAND ARC, Issue 4 2005Ryota Mori Abstract Structures developed in metamorphic and plutonic blocks that occur as knockers in the Mineoka Ophiolite Belt in the Boso Peninsula, central Japan, were analyzed. The aim was to understand the incorporation processes of blocks of metamorphic and plutonic rocks with an arc signature into the serpentinite mélange of the Mineoka Ophiolite Belt in relation to changes in metamorphic conditions during emplacement. Several stages of deformation during retrogressive metamorphism were identified: the first faulting stage had two substage shearing events (mylonitization) under ductile conditions inside the crystalline blocks in relatively deeper levels; and the second stage had brittle faulting and brecciation along the boundaries between the host serpentinite bodies in relatively shallower levels (zeolite facies). The first deformation occurred during uplift before emplacement. The blocks were intensively sheared by the first deformation event, and developed numerous shear planes with spacing of a few centimeters. The displacement and width of each shear plane were a few centimeters and a few millimeters, respectively, at most. In contrast, the fault zone of the second shearing stage reached a few meters in width and developed during emplacement of the Mineoka Ophiolite. Both stages occurred under a right-lateral transpressional regime, in which thrust-faulting was associated with strike-slip faulting. Such displacement on an outcrop scale is consistent with the present tectonics of the Mineoka Belt. This implies that the same tectonic stress has been operating in the Boso trench,trench,trench-type triple junction area in the northwest corner of the Pacific since the emplacement of the Mineoka Ophiolite. The Mineoka Ophiolite Belt must have worked as a forearc sliver fault during the formation of a Neogene accretionary prism further south. [source] Fault rock analysis of the northern part of the Chelungpu Fault and its relation to earthquake faulting of the 1999 Chi-Chi earthquake, TaiwanISLAND ARC, Issue 1 2005Kohtaro UjiieArticle first published online: 3 MAR 200 Abstract The 1999 Chi-Chi earthquake in Taiwan (Mw = 7.6) produced a surface rupture along the north,south-striking Chelungpu thrust fault with pure dip-slip (east side up) and left lateral strike-slip displacements. Near-field strong-motion data for the northern part of the fault illustrate a distinct lack of the high-frequency seismic radiation associated with a large slip (10,15 m) and a rapid slip velocity (2,4 m/s), suggesting a smooth seismic slip associated with low dynamic frictional resistance on the fault. A drillhole was constructed at shallow depths in the possible fault zones of the northern part of the Chelungpu Fault, which may have slipped during the 1999 earthquake. One of the zones consists of a 20-cm-thick, unconsolidated fault breccia with a chaotic texture lacking both discrete slip surfaces (e.g. Riedel shears) and grain crushing. Other possible fault zones are marked by the narrow (less than a few centimeters) gouge zone in which clayey material intrudes into the damaged zone outside of the gouge zone. These characteristic fault rock textures suggest that the slip mechanisms at shallow levels during the earthquake involved either granular flow of initially unconsolidated material or slip localization under elevated pore pressure along the narrow clayey gouge zone. Because both mechanisms lead to low dynamic frictional resistance on the fault, the rapid seismic slip in the deep portions of the fault (i.e. the source region of strong-motion radiation) could have been accommodated by frictionless slip on the shallow portions of the fault. The combination of strong-motion data and fault rock analysis suggests that smooth slip associated with low dynamic friction occurred on both the deep and shallow portions of the fault, resulting in a large slip between the source region and the surface in the northern region. [source] Monazite geochronology in central New England: evidence for a fundamental terrane boundaryJOURNAL OF METAMORPHIC GEOLOGY, Issue 3 2008F. S. SPEAR Abstract Monazite crystallization ages have been measured in situ using SIMS and EMP analysis of samples from the Bronson Hill anticlinorium in central New England. In west-central New Hampshire, each major tectonic unit (nappe) displays a distinctive P,T path and metamorphic history that requires significant post-metamorphic faulting to place them in their current juxtaposition, and monazite ages were determined to constrain the timing of metamorphism and nappe assembly. Monazite ages from the low-pressure, high-temperature Fall Mountain nappe range from c. 455 to 355 Ma, and Y zoning indicates that these ages comprise three to four distinct age domains, similar to that found in the overlying Chesham Pond nappe. The underlying Skitchewaug nappe contains monazite ages that range from c. 417 to 307 Ma. 40Ar/39Ar ages indicate rapid cooling of the Chesham Pond and Fall Mountain nappes after 350 Ma, which is believed to represent the time of emplacement of the high-level Chesham Pond and Fall Mountain nappes onto rocks of the underlying Skitchewaug nappe. Garnet zone rocks from western New Hampshire contain monazite that display a range of ages (c. 430,340 Ma). Both the metamorphic style and monazite ages suggest that the low-grade belt in western New Hampshire is continuous with the Vermont sequence to the west. Rocks of the Big Staurolite nappe in western New Hampshire contain monazite that crystallized between c. 370 and 290 Ma and the same unit along strike in northern New Hampshire and central Connecticut records ages of c. 257,300 Ma. Conspicuously absent from this nappe are the older age populations that are found in both the overlying nappes and underlying garnet zone rocks. These monazite ages confirm that the metamorphism observed in the Big Staurolite nappe occurred significantly later than that in the units structurally above and below. These data support the hypothesis that the Big Staurolite nappe represents a major tectonic boundary, along which rocks of the New Hampshire metamorphic series were juxtaposed against rocks of the Vermont series during the Alleghanian. [source] Controls on low-pressure anatexisJOURNAL OF METAMORPHIC GEOLOGY, Issue 2 2006C. C. GERBI Abstract Low-pressure anatexis, whereby rocks melt in place after passing through the andalusite stability field, develops under more restricted conditions than does low-pressure metamorphism. Our thermal modelling and review of published work indicate that the following mechanisms, operating alone, may induce anatexis in typical pelitic rocks without inducing wholesale melting in the lower crust: (i) magmatic advection by pervasive flow; (ii) crustal-scale detachment faulting; and (iii) the presence of a high heat-producing layer. Of these, only magmatic advection by pervasive flow and crustal-scale detachment faulting have been shown quantitatively to provide sufficient heat to cause widespread melting. Combinations of the above mechanisms with pluton-scale magmatic advection, shear heating, removal of the lithospheric mantle, or with each other provide additional means of developing suitable high temperatures at shallow crustal levels to generate low-pressure anatexis. [source] Pressure,temperature,time evolution of the Central East Greenland Caledonides: quantitative constraints on crustal thickening and synorogenic extensionJOURNAL OF METAMORPHIC GEOLOGY, Issue 9 2003A. P. White Abstract Whereas geologists have known for three-quarters of a century that there was significant crustal thickening in the central East Greenland Caledonides, the crucial role of extensional faulting during Caledonian orogenesis has only been recognized during the past decade. In this paper, new petrographic and thermobarometric observations are presented from migmatitic metasedimentary gneisses of the Forsblad Fjord region (c. 72.5°N). Samples of the Krummedal Sequence, collected from the footwall of the upper of two significant splays of the main extensional fault system in the region,the Fjord Region Detachment (FRD),enable us to establish a relative sequence of metamorphism. Our pressure (P),temperature (T) results imply a clockwise loop in P,T space. As recorded by mineral assemblages in the Krummedal gneisses, prograde metamorphism involved a net increase of c. 4 kbar and 250 °C, with peak conditions of c. 10.5 kbar at 785 °C. Early burial and heating was followed by near-isothermal decompression of 4.5 kbar, a process which is attributed to roughly 18 km of tectonostratigraphic throw on the upper splay of the FRD. Combining data reported here with the published data, it is estimated that the approximate tectonostratigraphic throw along the lower splay of the FRD was c. 16 km. In situ U,Th,Pb-monazite electron microprobe dating suggests that the earliest phase of metamorphism recorded in the Krummedal Sequence gneisses of Forsblad Fjord occurred during the Caledonian orogeny. Furthermore, the combination of our new data with existing conventional TIMS U-Pb and 40Ar/39Ar data imply that: (1) movement along the uppermost splay of the FRD (c. 425,423 Ma) occurred at maximum time-averaged slip-rates equivalent to c. 9 mm of vertical displacement per year; and (2) that the final stages of metamorphism occurred prior to c. 411 Ma, although part of this denudation was likely accommodated on overlying extensional structures that may have been active more recently. There is close agreement between our data and results from the Krummedal Sequence north of the field area (72.5°,74°N), and rocks of the Smallefjord Sequence (75°,76°N) that are suggested to correlate with the Krummedal Sequence. This leads us to infer that the events recorded in the Forsblad Fjord region are of orogen-scale significance. [source] Regional variation in exhumation and strain rate of the high-pressure Sambagawa metamorphic rocks in central Shikoku, south-west JapanJOURNAL OF METAMORPHIC GEOLOGY, Issue 7 2002K. Yagi Abstract Regional variation in the P,T path of the Sambagawa metamorphic rocks, central Shikoku, Japan has been inferred from compositional zoning of metamorphic amphibole. Rocks constituting the northern part (Saruta River area) exhibit a hairpin type P,T path, where winchite/actinolite grew at the prograde stage, the peak metamorphism was recorded by the growth of barroisite to hornblende and sodic amphibole to winchite/actinolite grew at the retrograde stage. In the southern part (Asemi River area), rocks exhibit a clockwise type P,T path, where barroisite to hornblende core is rimmed by winchite to actinolite. The difference in P,T path could suggest a faster exhumation rate (i.e. more rapid decompression) in the southern than in the northern part. On the other hand, physical conditions of deformation during the exhumation stage have been independently inferred from microstructures in deformed quartz. Recrystallized quartz grains in rocks from the low-grade (chlorite and garnet) zones are much more stretched in the southern part (aspect ratio , 4.0) than in the northern part (aspect ratio< 4.0), indicating a higher strain rate in the former than in the latter. These facts may indicate that the exhumation and strain rates are correlated (i.e. the exhumation rate increases with increasing the strain rate). The difference in the exhumation rate inferred from amphibole zoning between the northern and southern parts could be explained by an extensional model involving normal faulting, where the lower plate can be exhumed faster than the upper plate due to the displacement along the fault. Furthermore, the model may explain the positive correlation between the exhumation and strain rates, because the lower plate tended to support more stress than the upper plate. [source] P,T modelling of the andalusite,kyanite,andalusite sequence and related assemblages in high-Al graphitic pelites.JOURNAL OF METAMORPHIC GEOLOGY, Issue 6 2001Prograde, retrograde paths in a late kyanite belt in the Variscan Iberia Abstract The exceptional andalusite,kyanite,andalusite sequence occurs in Al-rich graphitic slates in a narrow pelite belt on the hangingwall of a ductile normal fault in NW Variscan Iberia. Early chiastolite is replaced by Ky,Ms,Pg aggregates, which are overgrown by pleochroic andalusite near granites intruded along the fault. Slates plot in AKFM above the chloritoid-chlorite tie-line. Their P,T grids are modelled with Thermocalc v2.7 and the 1998 databases in the NaKFMASH and KFMASH systems. The univariant reaction Ctd + And/Ky = St + Chl + Qtz + H2O ends at progressively lower pressure as F/FM increases and A/AFM decreases, shrinking the assemblage Cld,Ky,Chl, and opening a chlorite-free Cld,Ky trivariant field on the low temperature reaction side. This modelling matches the observed absence of chlorite in high F/FM rocks, which is restricted to low pressure in the andalusite stability field. The P,T path deduced from modelling shows a first prograde event in the andalusite field followed by retrogression into the kyanite field, most likely coupled with a slight pressure increase. The final prograde evolution into the andalusite field can be explained by two different prograde paths. Granite intrusion caused the first prograde part of the path with andalusite growth. The subsequent thermal relaxation, together with aH2O decrease, generated the retrograde andalusite,kyanite transformation, plus chlorite consumption and chloritoid growth. This transformation could have been related to folding in the beginning, and aided later by downthrowing due to normal faulting. Heat supplied by syntectonic granite intrusion explains the isobaric part of the path in the late stages of evolution, causing the prograde andalusite growth after the assemblage St,Ky,Chl. Near postectonic granites, a prograde path with pressure decrease originated the assemblage St,And,Chl. [source] | ||||||||||||||||||||||||||||