Crustal heterogeneities and complexities of fault processes

In this work we study the relation between crustal heterogeneities and complexities in fault processes. The first kind of heterogeneity considered involves the concept of asperity. The presence of an asperity in the hypocentral region of the M = 6.5 earthquake of June 17-th, 2000 in the South Iceland Seismic Zone was invoked to explain the change of seismicity pattern before and after the mainshock: in particular, the spatial distribution of foreshock epicentres trends NW while the strike of the main fault is N 7◦ E and aftershocks trend accordingly; the foreshock depths were typically deeper than average aftershock depths. A model is devised which simulates the presence of an asperity in terms of a spherical inclusion, within a softer elastic medium in a transform domain with a deviatoric stress field imposed at remote distances (compressive NE − SW, tensile NW − SE). An isotropic compressive stress component is induced outside the asperity, in the direction of the compressive stress axis, and a tensile component in the direction of the tensile axis; as a consequence, fluid flow is inhibited in the compressive quadrants while it is favoured in tensile quadrants. Within the asperity the isotropic stress vanishes but the deviatoric stress increases substantially, without any significant change in the principal stress directions. Hydrofracture processes in the tensile quadrants and viscoelastic relaxation at depth may contribute to lower the effective rigidity of the medium surrounding the asperity. According to the present model, foreshocks may be interpreted as induced, close to the brittle-ductile transition, by high pressure fluids migrating upwards within the tensile quadrants; this process increases the deviatoric stress within the asperity which eventually fails, becoming the hypocenter of the mainshock, on the optimally oriented fault plane. In the second part of our work we study the complexities induced in fault processes by the layered structure of the crust. In the first model proposed we study the case in which fault bending takes place in a shallow layer. The problem can be addressed in terms of a deep vertical planar crack, interacting with a shallower inclined planar crack. An asymptotic study of the singular behaviour of the dislocation density at the interface reveals that the density distribution has an algebraic singularity at the interface of degree ω between -1 and 0, depending on the dip angle of the upper crack section and on the rigidity contrast between the two media. From the welded boundary condition at the interface between medium 1 and 2, a stress drop discontinuity condition is obtained which can be fulfilled if the stress drop in the upper medium is lower than required for a planar trough-going surface: as a corollary, a vertically dipping strike-slip fault at depth may cross the interface with a sedimentary layer, provided that the shallower section is suitably inclined (fault "refraction"); this results has important implications for our understanding of the complexity of the fault system in the SISZ; in particular, we may understand the observed offset of secondary surface fractures with respect to the strike direction of the seismic fault. The results of this model also suggest that further fractures can develop in the opposite quadrant and so a second model describing fault branching in the upper layer is proposed. As the previous model, this model can be applied only when the stress drop in the shallow layer is lower than the value prescribed for a vertical planar crack surface. Alternative solutions must be considered if the stress drop in the upper layer is higher than in the other layer, which may be the case when anelastic processes relax deviatoric stress in layer 2. In such a case one through-going crack cannot fulfil the welded boundary conditions and unwelding of the interface may take place. We have solved this problem within the theory of fracture mechanics, employing the boundary element method. The fault terminates against the interface in a T-shaped configuration, whose segments interact among each other: the lateral extent of the unwelded surface can be computed in terms of the main fault parameters and the stress field resulting in the shallower layer can be modelled. A wide stripe of high and nearly uniform shear stress develops above the unwelded surface, whose width is controlled by the lateral extension of unwelding. Secondary shear fractures may then open within this stripe, according to the Coulomb failure criterion, and the depth of open fractures opening in mixed mode may be computed and compared with the well studied fault complexities observed in the field. In absence of the T-shaped decollement structure, stress concentration above the seismic fault would be difficult to reconcile with observations, being much higher and narrower.

Aplicación de la técnica de inyección del terreno por fracturación hidráulica a las cimentaciones de puentes en circunstancias de especial dificultad

La implantación de una gran obra de paso facilita en multitud de ocasiones la vida de aquellos que aprovechen las ventajas que ofrece esta estructura. Sin embargo, el camino que lleva a su construcción supone importantes y tortuosos retos cuando el terreno de cimentación presenta condiciones geotécnicas no compatibles con las grandes cargas que el viaducto necesita transmitirle para su buen funcionamiento. El trabajo que desarrolla esta Tesis Doctoral proporciona una herramienta eficaz y económica, por lo reducido de su extensión y medios, que permite allanar el camino que acomete la construcción de una estructura de tal envergadura. Mediante el análisis de la problemática del terreno y de las distintas soluciones de cimentación empleadas en la actualidad, se conduce al lector hacia una técnica de cimentación innovadora que combina la técnica del micropilotaje, para la canalización de las cargas estructurales, junto a la técnica de la inyección de fracturación hidráulica, que mejorará el terreno de implantación de los micropilotes, conformando así un medio capaz de recibir y transmitir grandes cargas en cualquier tipo de terreno. La técnica ya empleada del micropilotaje, por su trabajo esencialmente axil, requiere, aunque sin problema, la constitución de sistemas de fuerzas que equilibren el sistema de cargas provenientes del viaducto. Pero su capacidad resistente viene condicionada por el terreno circundante. Cuando la roca rodea el micropilote, su empleo se realiza sin problema y sin necesidad de mejorar el terreno. Pero sin terreno consistente, el empleo de inyecciones de fracturación hidráulica a través de los propios micropilotes, no sólo mejorará la capacidad resistente de los micropilotes, aumentando la inercia necesaria cuando las cargas son de origen sísmico o estructural ferroviario, sino que resolverá y eliminará los problemas de estabilidad que presentan las laderas que frecuentemente deben recibir las cargas de la obra de paso. Tras recoger el análisis ya realizado en el Trabajo de Investigación, donde se justificaba la alta capacidad resistente de un micropilote con terreno circundante mejorado por la inyección, en la actual Tesis se emplean modelos matemáticos sobre un caso real de viaducto sometido a un gran sismo e implantado en una enorme quebrada (500 m) de Colombia, sujeta a lluvias torrenciales y en la zona de mayor sismicidad del país. Con ello se comprueba la estabilización que se alcanza en el terreno de cimentación con el empleo de esta técnica de transmisión de carga y mejora del terreno. De esta forma se completa un ciclo que justifica las bondades de esta combinación de técnicas de cimentación, pero se abren las puertas a nuevos entornos de aplicación, como edificios antiguos de cualquier tipo que requieran recalces, y no sólo en la implantación de grandes obras de paso. ABSTRACT The establishment of a large bridge represents, in many cases, a better life for those who can take profit of the advantages provided by that structure. Nevertheless, the process of building this structure has to overcome important and difficult circumstances whenever the geotechnical conditions of the bridge site are not adequate to carry the large loads transmitted by the bridge structure. This study develops a method both effective and economical, due to the extension and means necessary for its application, which allows to solve properly the foundation of a structure of that importance. Considering the geotechnical problems inherent to the bridge site, along with the different foundation solutions that are presently used and their limitations, the study leads the reader to an innovative technique which combines the micropile system, for transmission of the structural loads of the bridge, with the technique of hydraulic fracture grouting for improvement of the ground around the micropiles, allowing to both stabilizing and transmitting large loads in any kind of ground. It is well known that the micropiles work axially, and this condition requires an adequate distribution of those units, in order to properly absorb the load system introduced by the viaduct. The resistance of the bridge foundation is, in any case, provided by the ground. When rock is encountered, the micropiles have been successfully used without improving the ground. However, as it is shown in this study, by using the micropiles as sleeve-pipes for hydraulic fracture grouting, not only the micropile resistance can be improved in any ground, but it is possible to develop grouted “solids” in the ground, whose inertia allows to absorb actions of structural and seismic origin. Additionally, as it is shown and analyzed in the Thesis, the ground improvement can give an adequate safety factor to the slopes frequently encountered in bridge sites. In order to properly justify those advantages of combining micropiles with ground improvement through fracture grouting around the micropiles, mathematical models have been developed and applied to a real case of a cable-stayed bridge installed on a very large ravine (500 m) in Colombia, located in the highest seismic zone of the country, and subject to torrential rains. The results of this numerical analysis show the high safety condition provided by the ground improvement to the viaduct site. In conclusion, the Thesis shows the important improvement that can be provided by the combination of micropiles and soil improvement, through fracture grouting, to the problem of founding bridges. However, it can be understood that this technique could be applied successfully to underpinning buildings, specially old buildings of any type, apart from its use in bridge foundations.

Late-Holocene Mass Movements in High Arctic East Lake, Melville Island (Western Canadian Arctic Archipelago)

East Lake, located at Cape Bounty (Melville Island, Canadian High Arctic), was mapped using a high-resolution swath bathymetric sonar and a 12 kHz sub-bottom profiler, allowing for the first time the imaging of widespread occurrence of mass movement deposits (MMDs) in a Canadian High Arctic Lake. Mass movements occurred mostly on steep slopes located away from deltaic sedimentation. The marine to lacustrine transition in the sediment favours the generation of mass movements where the underlying massive mud appears to act as a gliding surface for the overlying varved deposits. Based on acoustic stratigraphy, we have identified at least two distinct events that triggered failures in the lake during the last 2000 years. The synchronicity of multiple failures and their widespread distribution suggest a seismic origin that could be related to the nearby Gustaf-Lougheed Arch seismic zone. Further sedimentological investigations on the MMDs are however required to confirm their age and origin.

Late-Holocene Mass Movements in High Arctic East Lake, Melville Island (Western Canadian Arctic Archipelago)

East Lake, located at Cape Bounty (Melville Island, Canadian High Arctic), was mapped using a high-resolution swath bathymetric sonar and a 12 kHz sub-bottom profiler, allowing for the first time the imaging of widespread occurrence of mass movement deposits (MMDs) in a Canadian High Arctic Lake. Mass movements occurred mostly on steep slopes located away from deltaic sedimentation. The marine to lacustrine transition in the sediment favours the generation of mass movements where the underlying massive mud appears to act as a gliding surface for the overlying varved deposits. Based on acoustic stratigraphy, we have identified at least two distinct events that triggered failures in the lake during the last 2000 years. The synchronicity of multiple failures and their widespread distribution suggest a seismic origin that could be related to the nearby Gustaf-Lougheed Arch seismic zone. Further sedimentological investigations on the MMDs are however required to confirm their age and origin.

GPS inferred velocity and strain rate fields in eastern Canada

La vallée du fleuve Saint-Laurent, dans l’est du Canada, est l’une des régions sismiques les plus actives dans l’est de l’Amérique du Nord et est caractérisée par de nombreux tremblements de terre intraplaques. Après la rotation rigide de la plaque tectonique, l’ajustement isostatique glaciaire est de loin la plus grande source de signal géophysique dans l’est du Canada. Les déformations et les vitesses de déformation de la croûte terrestre de cette région ont été étudiées en utilisant plus de 14 ans d’observations (9 ans en moyenne) de 112 stations GPS fonctionnant en continu. Le champ de vitesse a été obtenu à partir de séries temporelles de coordonnées GPS quotidiennes nettoyées en appliquant un modèle combiné utilisant une pondération par moindres carrés. Les vitesses ont été estimées avec des modèles de bruit qui incluent les corrélations temporelles des séries temporelles des coordonnées tridimensionnelles. Le champ de vitesse horizontale montre la rotation antihoraire de la plaque nord-américaine avec une vitesse moyenne de 16,8±0,7 mm/an dans un modèle sans rotation nette (no-net-rotation) par rapport à l’ITRF2008. Le champ de vitesse verticale confirme un soulèvement dû à l’ajustement isostatique glaciaire partout dans l’est du Canada avec un taux maximal de 13,7±1,2 mm/an et un affaissement vers le sud, principalement au nord des États-Unis, avec un taux typique de −1 à −2 mm/an et un taux minimum de −2,7±1,4 mm/an. Le comportement du bruit des séries temporelles des coordonnées GPS tridimensionnelles a été analysé en utilisant une analyse spectrale et la méthode du maximum de vraisemblance pour tester cinq modèles de bruit: loi de puissance; bruit blanc; bruit blanc et bruit de scintillation; bruit blanc et marche aléatoire; bruit blanc, bruit de scintillation et marche aléatoire. Les résultats montrent que la combinaison bruit blanc et bruit de scintillation est le meilleur modèle pour décrire la partie stochastique des séries temporelles. Les amplitudes de tous les modèles de bruit sont plus faibles dans la direction nord et plus grandes dans la direction verticale. Les amplitudes du bruit blanc sont à peu près égales à travers la zone d’étude et sont donc surpassées, dans toutes les directions, par le bruit de scintillation et de marche aléatoire. Le modèle de bruit de scintillation augmente l’incertitude des vitesses estimées par un facteur de 5 à 38 par rapport au modèle de bruit blanc. Les vitesses estimées de tous les modèles de bruit sont statistiquement cohérentes. Les paramètres estimés du pôle eulérien de rotation pour cette région sont légèrement, mais significativement, différents de la rotation globale de la plaque nord-américaine. Cette différence reflète potentiellement les contraintes locales dans cette région sismique et les contraintes causées par la différence des vitesses intraplaques entre les deux rives du fleuve Saint-Laurent. La déformation de la croûte terrestre de la région a été étudiée en utilisant la méthode de collocation par moindres carrés. Les vitesses horizontales interpolées montrent un mouvement cohérent spatialement: soit un mouvement radial vers l’extérieur pour les centres de soulèvement maximal au nord et un mouvement radial vers l’intérieur pour les centres d’affaissement maximal au sud, avec une vitesse typique de 1 à 1,6±0,4 mm/an. Cependant, ce modèle devient plus complexe près des marges des anciennes zones glaciaires. Basées selon leurs directions, les vitesses horizontales intraplaques peuvent être divisées en trois zones distinctes. Cela confirme les conclusions d’autres chercheurs sur l’existence de trois dômes de glace dans la région d’étude avant le dernier maximum glaciaire. Une corrélation spatiale est observée entre les zones de vitesses horizontales intraplaques de magnitude plus élevée et les zones sismiques le long du fleuve Saint-Laurent. Les vitesses verticales ont ensuite été interpolées pour modéliser la déformation verticale. Le modèle montre un taux de soulèvement maximal de 15,6 mm/an au sud-est de la baie d’Hudson et un taux d’affaissement typique de 1 à 2 mm/an au sud, principalement dans le nord des États-Unis. Le long du fleuve Saint-Laurent, les mouvements horizontaux et verticaux sont cohérents spatialement. Il y a un déplacement vers le sud-est d’une magnitude d’environ 1,3 mm/an et un soulèvement moyen de 3,1 mm/an par rapport à la plaque l’Amérique du Nord. Le taux de déformation verticale est d’environ 2,4 fois plus grand que le taux de déformation horizontale intraplaque. Les résultats de l’analyse de déformation montrent l’état actuel de déformation dans l’est du Canada sous la forme d’une expansion dans la partie nord (la zone se soulève) et d’une compression dans la partie sud (la zone s’affaisse). Les taux de rotation sont en moyenne de 0,011°/Ma. Nous avons observé une compression NNO-SSE avec un taux de 3.6 à 8.1 nstrain/an dans la zone sismique du Bas-Saint-Laurent. Dans la zone sismique de Charlevoix, une expansion avec un taux de 3,0 à 7,1 nstrain/an est orientée ENE-OSO. Dans la zone sismique de l’Ouest du Québec, la déformation a un mécanisme de cisaillement avec un taux de compression de 1,0 à 5,1 nstrain/an et un taux d’expansion de 1.6 à 4.1 nstrain/an. Ces mesures sont conformes, au premier ordre, avec les modèles d’ajustement isostatique glaciaire et avec la contrainte de compression horizontale maximale du projet World Stress Map, obtenue à partir de la théorie des mécanismes focaux (focal mechanism method).

Weak Zone Related Seismic Cycles in Progressive Failure Leading to Collapse in Brittle Crust

Until quite recently our understanding of the basic mechanical process responsible for earthquakes and faulting was not well known. It can be argued that this was partly a consequence of the complex nature of fracture in crust and in part because evidence of brittle phenomena in the natural laboratory of the earth is often obliterated or obscured by other geological processes. While it is well understood that the spatial and temporal complexity of earthquakes and the fault structures emerge from geometrical and material built-in heterogeneities, one important open question is how the shearing becomes localized into a band of intense fractures. Here the authors address these questions through a numerical approach of a tectonic plate by considering rockmass heterogeneity both in microscopic scale and in mesoscopic scale. Numerical simulations of the progressive failure leading to collapse under long-range slow driving forces in the far-field show earthquake-like rupture behavior. $En Echelon$ crack-arrays are reproduced in the numerical simulation. It is demonstrated that the underlying fracturing induced acoustic emissions (or seismic events) display self-organized criticality------from disorder to order. The seismic cycles and the geometric structures of the fracture faces, which are found greatly depending on the material heterogeneity (especially on the macroscopic scale), agree with that observed experimentally in real brittle materials. It is concluded that in order to predict a main shock, one must have extremely detailed knowledge on very minor features of the earth's crust far from the place where the earthquake originated. If correct, the model proposed here seemingly provides an explanation as to why earthquakes to date are not predicted so successfully. The reason is not that the authors do not understand earthquake mechanisms very well but that they still know little about our earth's crust.

Analysis on deep crustal structure along the onshore-offshore seismic profile across the Binhai (littoral) fault zone in northeastern South China Sea

The onshore-offshore deep seismic experiment was carried out for the first time and filled the blankness of the seismic surveys in the transition area between South China and northeastern South China Sea. The seismic data were analyzed and processed. The different seismic phases were identified and their travel time arrivals were modeled by ray-tracing to study the P-wave velocity crustal structure of this area. The crustal structure of this area is the continental crust. The crust thickness is gradually decreasing southward along the on-shore-offshore seismic line. The low-velocity layer (5.5 similar to 5.9 km (.) s(-1)) exists generally in the middle crust (about 10.0 similar to 18.0km)with about 2.5 similar to 4.0 km thickness, which is also thinning seaward. No obvious high-velocity layer appears in the lower crust. The Binhai (littoral) fault zone is a low velocity zone, which is located about 35km southeast to the Nan'ao station and corresponding to the gradient belt of gravity & magnetism anomalies. The depth of the fault zone is close to the Moho discontinuity. The littoral fault zone is a boundary between the normal continental crust of South China and the thinned continental crust of the sea area.

Seismic profiles across the middle Tan-Lu fault zone in Laizhou Bay, Bohai Sea, eastern China

The stratigraphic architecture, structure and Cenozoic tectonic evolution of the Tan-Lu fault zone in Laizhou Bay, eastern China, are analyzed based on interpretations of 31 new 2D seismic lines across Laizhou Bay. Cenozoic strata in the study area are divided into two layers separated by a prominent and widespread unconformity. The upper sedimentary layer is made up of Neogene and Quaternary fluvial and marine sediments, while the lower layer consists of Paleogene lacustrine and fluvial facies. In terms of tectonics, the sediments beneath the unconformity can be divided into four main structural units: the west depression, central uplift, east depression and Ludong uplift. The two branches of the middle Tan-Lu fault zone differ in their geometry and offset: the east branch fault is a steeply dipping S-shaped strike-slip fault that cuts acoustic basement at depths greater than 8 km, whereas the west branch fault is a relatively shallow normal fault. The Tan-Lu fault zone is the key fault in the study area, having controlled its Cenozoic evolution. Based on balanced cross-sections constructed along transverse seismic line 99.8 and longitudinal seismic line 699.0, the Cenozoic evolution of the middle Tan-Lu fault zone is divided into three stages: Paleocene-Eocene transtension, Oligocene-Early Miocene transpression and Middle Miocene to present-day stable subsidence. The reasons for the contrasting tectonic features of the two branch faults and the timing of the change from transtension to transpression are discussed. Crown Copyright (C) 2008 Published by Elsevier Ltd. All rights reserved.

Seismic anisotropy in the Morcles nappe shear zone: Implications for seismic imaging of crustal scale shear zones

Microstructures and textures of calcite mylonites from the Morcles nappe large-scale shear zone in southwestern Switzerland develop principally as a function of 1) extrinsic physical parameters including temperature, stress, strain, strain rate and 2) intrinsic parameters, such as mineral composition. We collected rock samples at a single location from this shear zone, on which laboratory ultrasonic velocities, texture and microstructures were investigated and quantified. The samples had different concentration of secondary mineral phases (< 5 up to 40 vol.%). Measured seismic P wave anisotropy ranges from 6.5% for polyphase mylonites (~ 40 vol.%) to 18.4% in mylonites with < 5 vol.% secondary phases. Texture strength of calcite is the main factor governing the seismic P wave anisotropy. Measured S wave splitting is generally highest in the foliation plane, but its origin is more difficult to explain solely by calcite texture. Additional texture measurements were made on calcite mylonites with low concentration of secondary phases (≤ 10 vol.%) along the metamorphic gradient of the shear zone (15 km distance). A systematic increase in texture strength is observed moving from the frontal part of the shear zone (anchimetamorphism; 280 °C) to the higher temperature, basal part (greenschist facies; 350–400 °C). Calculated P wave velocities become increasingly anisotropic towards the high-strain part of the nappe, from an average of 5.8% in the frontal part to 13.2% in the root of the basal part. Secondary phases raise an additional complexity, and may act either to increase or decrease seismic anisotropy of shear zone mylonites. In light of our findings we reinterpret the origin of some seismically reflective layers in the Grône–Zweisimmen line in southwestern Switzerland (PNR20 Swiss National Research Program). We hypothesize that reflections originate in part from the lateral variation in textural and microstructural arrangement of calcite mylonites in shear zones.

Study of the evolution of the seismic cycle of stress and strain associated to the El Salvador Fault Zone

• Central America: – Regional studies in Central America (Seismic Hazard). – El Salvador Fault Zone (ESFZ). – Aguacaliente‐Navarro Fault Zone (ANFZ), Central Valley of Costa Rica. – Haiti (seismic hazard) • Spain: – Regional‐Nacional studies of seismic hazards (applications to building codes, eurocode, emergency plans, etc.) – Betic range zone, south of Spain. – Ibero‐Maghrebi region (collision zone)

New constraints on micro-seismicity and stress state in the western part of the North Anatolian Fault Zone : Observations from a dense seismic array

Major funding was provided by the UK Natural Environment Research Council (NERC) under grant NE/I028017/1 and partially supported by Boğaziçi University Research Fund (BAP) under grant 6922. We would like to thank all the project members from the University of Leeds, Boğaziçi University, Kandilli Observatory, Aberdeen University and Sakarya University. I would also like to thank Prof. Ali Pinar and Dr. Kıvanç Kekovalı for their valuable comments. Some of the figures were generated by GMT software (Wessel and Smith, 1995).

New constraints on micro-seismicity and stress state in the western part of the North Anatolian Fault Zone : Observations from a dense seismic array

Major funding was provided by the UK Natural Environment Research Council (NERC) under grant NE/I028017/1 and partially supported by Boğaziçi University Research Fund (BAP) under grant 6922. We would like to thank all the project members from the University of Leeds, Boğaziçi University, Kandilli Observatory, Aberdeen University and Sakarya University. I would also like to thank Prof. Ali Pinar and Dr. Kıvanç Kekovalı for their valuable comments. Some of the figures were generated by GMT software (Wessel and Smith, 1995).

Evidence of episodic brittle faulting in the cratonic part of the Peninsular India and its implications for seismic hazard in slow deforming regions

The region around Waclakkancheri, in the province of Kerala, India, which lies in the vicinity of Palghat-Cauvery ;hear zone (within the Precambrian crystalline terrain), has been a site of microseismic activity since 1989. Earlier studies had identified a prominent WNW-ESE structure overprinting on the E-W trending lineaments associated with Palghat-Cauvery shear zone. We have mapped this structure, located in a chamockite quarry near Desamangalam, Waclakkancheri, which we identify as a ca. 30 km-long south dipping reverse fault. This article presents the characteristics of this fault zone exposed on the exhumed crystalline basement and discusses its significance in understanding the earthquake potential of the region. This brittle deformation zone consists of fracture sets with small-scale displacement and slip planes with embedded fault gouges. The macroscopic as well as the microscopic studies of this fault zone indicate that it evolved through different episodes of faulting in the presence of fluids. The distinct zones within consolidated gouge and the cross cutting relationship of fractures indicate episodic fault activity. At least four faulting episodes can be recognized based on the sequential development of different structural elements in the fault rocks. The repeated ruptures are evident along this shear zone and the cyclic behavior of this fault consists of co-seismic ruptures alternating with inter-seismic periods, which is characterized by the sealed fractures and consolidated gouge. The fault zone shows a minimum accumulated dip/oblique slip of 2.1 m in the reverse direction with a possible characteristic slip of 52 cm (for each event). The ESR dating of fault gouge indicates that the deformation zone records a major event in the Middle Quaternary. The empirical relationships between fault length and slip show that this fault may generate events M >= 6. The above factors suggest that this fault may be characterized as potentially active. Our study offers some new pointers that can be used in other slow deforming cratonic hinterlands in exploring the discrete active faults.

Seismic hazard map of Coimbatore using subsurface fault rupture

This study presents the future seismic hazard map of Coimbatore city, India, by considering rupture phenomenon. Seismotectonic map for Coimbatore has been generated using past earthquakes and seismic sources within 300 km radius around the city. The region experienced a largest earthquake of moment magnitude 6.3 in 1900. Available earthquakes are divided into two categories: one includes events having moment magnitude of 5.0 and above, i.e., damaging earthquakes in the region and the other includes the remaining, i.e., minor earthquakes. Subsurface rupture character of the region has been established by considering the damaging earthquakes and total length of seismic source. Magnitudes of each source are estimated by assuming the subsurface rupture length in terms of percentage of total length of sources and matched with reported earthquake. Estimated magnitudes match well with the reported earthquakes for a RLD of 5.2% of the total length of source. Zone of influence circles is also marked in the seismotectonic map by considering subsurface rupture length of fault associated with these earthquakes. As earthquakes relive strain energy that builds up on faults, it is assumed that all the earthquakes close to damaging earthquake have released the entire strain energy and it would take some time for the rebuilding of strain energy to cause a similar earthquake in the same location/fault. Area free from influence circles has potential for future earthquake, if there is seismogenic source and minor earthquake in the last 20 years. Based on this rupture phenomenon, eight probable locations have been identified and these locations might have the potential for the future earthquakes. Characteristic earthquake moment magnitude (M-w) of 6.4 is estimated for the seismic study area considering seismic sources close to probable zones and 15% increased regional rupture character. The city is divided into several grid points at spacing of 0.01 degrees and the peak ground acceleration (PGA) due to each probable earthquake is calculated at every grid point in city by using the regional attenuation model. The maximum of all these eight PGAs is taken for each grid point and the final PGA map is arrived. This map is compared to the PGA map developed based on the conventional deterministic seismic hazard analysis (DSHA) approach. The probable future rupture earthquakes gave less PGA than that of DSHA approach. The occurrence of any earthquake may be expected in near future in these eight zones, as these eight places have been experiencing minor earthquakes and are located in well-defined seismogenic sources.

The hazard potential of the western segment of the Makran subduction zone, northern Arabian Sea

Evaluating the hazard potential of the Makran subduction zone requires understanding the previous records of the large earthquakes and tsunamis. We address this problem by searching for earthquake and tectonic proxies along the Makran Coast and linking those observations with the available constraints on historical seismicity and the tell-tale characteristics of sea floor morphology. The earthquake of Mw 8.1 of 1945 and the consequent tsunami that originated on the eastern part of the Makran are the only historically known hazardous events in this region. The seismic status of the western part of the subduction zone outside the rupture area of the 1945 earthquake remains an enigma. The near-shore shallow stratigraphy of the central part of Makran near Chabahar shows evidence of seismically induced liquefaction that we attribute to the distant effects of the 1945 earthquake. The coastal sites further westward around Jask are remarkable for the absence of liquefaction features, at least at the shallow level. Although a negative evidence, this possibly implies that the western part of Makran Coast region may not have been impacted by near-field large earthquakes in the recent past-a fact also supported by the analysis of historical data. On the other hand, the elevated marine terraces on the western Makran and their uplift rates are indicative of comparable degree of long-term tectonic activity, at least around Chabahar. The offshore data suggest occurrences of recently active submarine slumps on the eastern part of the Makran, reflective of shaking events, owing to the great 1945 earthquake. The ocean floor morphologic features on the western segment, on the contrary, are much subdued and the prograding delta lobes on the shelf edge also remain intact. The coast on the western Makran, in general, shows indications of progradation and uplift. The various lines of evidence thus suggest that although the western segment is potentially seismogenic, large earthquakes have not occurred there in the recent past, at least during the last 600 years. The recurrence period of earthquakes may range up to 1,000 years or more, an assessment based on the age of the youngest dated coastal ridge. The long elapsed time points to the fact that the western segment may have accumulated sufficient slip to produce a major earthquake.