nueva página del texto (beta)
Inglés (pdf)
Artículo en XML
Referencias del artículo
Traducción automática
Enviar artículo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
INTRODUCTION
The central part of Mexico is regularly affected by earthquakes related to subduction processes, volcanic activity, and crustal faults. In the past 500 years, severe crustal earthquakes with epicenters located on faults within the Trans-Mexican Volcanic Belt (TMVB), particularly on faults of the Chapala-Tula fault zone and the Taxco-San Miguel de Allende fault system, caused considerable human losses and material damage. Some examples are the following: in 1567, in Ameca, Mw 7.2 (Suter, 2015) or Mw 7.5-7.8 (Suárez et al., 1994); in 1875, in San Cristóbal de la Barranca, MI VII (Bárcena, 1875; Iglesias et al., 1877); in 1912 in Acambay, Mw 6.8-7 (Urbina and Camacho 1913; and Langridge et al., 2000); in 1920, in Xalapa, Ms 6.2-6.4 (Flores et al., 1922; Suter et al., 1996). In spite of these studies, crustal earthquakes have not been the objective of comprehensive researches (Zúñiga et al., 2003).
The TMVB is a volcanic arc located in central Mexico, between 19°N and 21°N, formed by the interaction of the Cocos, Rivera, and North-American plates since the Miocene (Mooser, 1972; Demant, 1978; Ferrari et al., 2012). The central part ofthe TMVB is deformed by the extensional Morelia-Acambay fault system (MAFS, Figure 1). This system belongs to the Chapala-Tula fault zone and includes some of the most active crustal faults in Mexico. This fault system is also influenced by the Quaternary activity of the Michoacán-Guanajuato Volcanic Field (Johnson and Harrison, 1990; Hasenaka and Carmichael, 1985).
Figure 1 Location of the Trans-Mexican Volcanic Belt (TMVB), the Michoacán-Guanajuato volcanic field (MGVF; yellow, dotted line) and tectonic structures of the Chapala-Tula fault zone and the Morelia-Acambay fault system (MAFS). Stars indicate sites with published paleoseismic studies and numbers indicate historical earthquakes: 1 = Pátzcuaro, 1845 and 1858; 2 = Morelia, 2007; 3 = Araró, 1845; 4 = Ucareo-Zinapécuaro-Acámbaro, 1872 and 1874; 5 = Maravatío, 1979; 6 = Venta de Bravo, 1734-1755 and 1853-1854; 7 = Acambay, 1912. Abbreviations: TVSF, Tzitzio-Valle de Santiago inferred fault (from Garduño-Monroy et al., 2009 and Arce et al., 2012); MCHF, Maravatío-Ciudad Hidalgo inferred fault (from Garduño-Monroy et al., 2009 and Arce et al., 2012); LATF = La Alberca-Teremendo Fault; CDMX = Mexico City (Images obtained from ArcMap 10.2.2 and modified from Gómez-Tuena et al., 2005).
Despite the seismic and paleoseismic studies carried out along the MAFS, most of the faults remain unreported and their seismogenic potential is not considered in urban planning. This omission is a serious problem considering the population increase in central Mexico. So far, paleoseismic studies have been concentrated in the eastern part of the MAFS, around the epicenter of the 1912 Acambay earthquake (Lacan et al., 2018; Velázquez-Bucio and Garduño-Monroy, 2018; Ortuño et al., 2019), leaving the western portion of the system with scarce information about recurrence intervals and potential magnitudes (Garduño-Monroy et al., 2009; Singh et al., 2012; Suter, 2016). This information is required for a reliable seismic hazard assessment throughout the system.
This work focuses on the characterization of the La Alberca-Teremendo fault activity (LATF, Figure 1) from a paleoseismic and neotectonic perspective. The seismogenic potential of the LATF has never been assessed, despite geomorphological evidence of recent activity, such as a linear escarpment inside a well-defined Quaternary graben, and displacement of soils and phreatomagmatic sequences of La Alberca de Guadalupe maar (LAG), dated as ~21 ka (Siebe et al., 2014; Kshirsagar et al., 2015).
The main purpose of this work is to identify and characterize the rupture events of the LATF, their potential magnitude Mw, and the associated seismic hazard. These parameters were identified via a paleoseismic analysis, emphasizing soils pedogenesis, combined with a geomorphologic characterization. We propose a work scheme that can be used for the study of other faults that displace recent soils.
BACKGROUND
The MAFS developed in a left-lateral, transtensional, stress field that accommodates deformation in the center of the TMVB, according to the oblique subduction model with slip partitioning proposed by Ego and Ansan (2002). The MAFS presents geomorphological features and superficial scarps related to historical seismicity (Suter et al., 2001). Aside from the 1912 Acambay earthquake, other historical crustal earthquakes in the MAFS occurred in Pátzcuaro, in 1845 with MI VIII (Sánchez-Garcilazo, 2000), Maravatío in 1979 with Mb = 5.3 (Astíz-Delgado, 1980; Garduño-Monroy and Gutiérrez-Negrín, 1992) and Morelia in 2007 with Mw = 4.3 (Garduño-Monroy et al., 2009; Singh et al., 2012). There are also some poorly-recorded seismic crises between 1750 and 1890 in the Venta de Bravo fault, and the regions of Araró and Ucareo-Zinapécuaro-Acámbaro (Arreola, 1985; Dobson and Mahood, 1985; Suter et al., 1992, 1996), as well as shallow crustal seismicity in march and october, 1976 in the Mezquital graben-Los Aljibes semi-graben (Quintanar et al., 2004). Additionally, an intermediate-depth, in-slab earthquake occurred in 1858 in Pátzcuaro with Mw = 7.7 that caused severe damage in the area (Singh et al., 1996) (Figure 1).
The paleoseismic studies performed in the Acambay graben show that most of the faults were active during Late Pleistocene to Holocene times, with slip rates of 0.03-0.37 mm/yr., recurrence intervals of 2000 to 15000 years, and potential seismic magnitudes (Mw) of 6.4-7.0 (Langridge et al., 2000, 2013; Ortuño et al., 2015; Sunye-Puchol et al., 2015; Lacan et al., 2018). In the Pátzcuaro-Morelia-Cuitzeo region, paleoseismic studies have been performed on fault traces and lacustrine sequences, reporting slip rates of 0.009-2.78 mm/yr., recurrence intervals of 1200 to 20000 years, and potential seismic magnitude (Mw) of 5.8-7.1 (Garduño-Monroy et al., 2009; Suter, 2016).
STUDY AREA
The LATF is located in the northern part of the state of Michoacán, 30 km northwest of the capital, Morelia, within the Michoacán-Guanajuato Volcanic Field (MGVF) and the western portion of the MAFS (Figure 2). The area is a transition zone between Miocene volcanic activity and the Quaternary volcanism of the TMVB which renders the faulting poorly detectable around Cuitzeo. The oldest rocks exposed in the area are Miocene andesite (~19 Ma) and dacitic ignimbrite of the Mil Cumbres sequence, followed by a Miocene-Pliocene fluvial-lacustrine sequence of the Cuitzeo paleo-basin, Pleistocene andesite and dacite from the El Tzirate Volcanic Complex and Holocene basaltic andesite of the MGVF (Figure 2b) (Pasquarè et al., 1991; Garduño-Monroy, 1999; Israde-Alcántara et al., 2010; Gómez-Vasconcelos et al., 2015; Pérez-Orozco et al., 2018).
Figure 2 a) Location of the La Alberca-Teremendo fault (LATF) showing the complex geometric array. b) Local geology of the study area, main structures and stress field. c) Aerial view of the trenched area. The dotted line indicates the limit of the basaltic-andesitic lava flows from the Cerro Pelón volcano which underlay the La Alberca de Guadalupe maar (LAG). View of the collapsed summit of El Picacho volcano.
The La Alberca-Teremendo fault displays a complex geometry. The fault is composed of eleven segments with an en échelon array, separated by less than one kilometer, and individual lengths that range between 2 and 10 km; the total length of the combined segments is 26 km (Figure 2a). It could be interpreted as an immature normal fault whose segments are still not merged, but under tectonic stress could be connected; however, we need more geophysical data in order to clarify the connection between segments. The LATF trends N85°E and dips 75°NW. This structure shows predominantly dip slip, but features such as the displacement of the LAG crater and deviation of the drainage network to the left indicate a minor left-lateral displacement.
The morphological expressions of the LATF are semi-graben geometries in Cuitzeo, formed by uplifted, rotated, and tilting blocks of Miocene ignimbrite and fluvial-lacustrine sequences; a counter slope in the southern flank of El Picacho volcano; and the cut-in-half LAG maar (Figure 2c). The fault generates a maximum height throw of ~ 50 m visible on the basaltic-andesite lavas of the LAG maar basement exposed in the crater; the minimum height throw is ~ 1-1.4 m, located at the western end of the fault, where the scarp puts the maar phreatomagmatic sequences in contact with current soil.
The LATF also displays characteristics similar to magma-induced extensional structures, like a geometrical array composed by subparallel normal faults with lengths < 14 km and dips between 70°-80°, but originates fault scarps of > 2 m, and 8°-15° tilting in the hanging wall. The western limit of the fault is a NNE-SSW eruptive fissure that distributes felsic material from the El Tzirate Volcanic Complex (Pérez-Orozco et al., 2018); the eastern limit is El Picacho semi-shield volcano and the faulting zone of the southern part of Cuitzeo lake, where large fault planes are observed without related dikes. The observation of dike-induced structures is more conspicuous around Queréndaro, east of the study area, where ten monogenetic cones are aligned on an E-W eruptive fracture, and internally sustained by a ~ 700 ka-old dike; here the electrical tomography displays intrusive bodies related to NNW-SSE structures (Gómez-Vasconcelos et al., 2018).
The study of the LATF area must take into consideration of the influence of both, the magma-induced structures of the MGVF and the tectonic structures of the MAFS, and its extensional regime. In the area, the stress could be accommodated either by magmatism, normal faulting, or both, and the evaluation of the activity could be problematic. We consider that, in the study area, magmatic processes influence local seismicity and that volcano-tectonic earthquakes have the capacity to change stresses within the shallow and middle crust. These events are large, purely tectonic earthquakes rather than destructive ones; however, their study is important to understand the seismic hazard (Hackett et al., 1996).
STUDY METHODS
Conventional morphologic and structural analyses were performed to identify seismic-related features, assess the influence of tectonically active faults in the landscape, and locate sites for paleoseismic analysis. All these phases were performed with caution because of the ambiguous characteristics of the LATF. We worked with topographic maps (1:50000 scale), aerial photos, and LIDAR images with 5 m-resolution in the ArcMap platform.
Bathymetric, magnetometric and seismo-acoustic surveys were performed inside the LAG crater. We obtained data from twenty-seven bathymetric profiles, seven seismo-acoustic reflection profiles and six magnetometric profiles (lines with NE-SW, N-S and NW-SE orientations). Bathymetry was performed using a Furuno 808 grapher ecosound and a Garmin ecosound GPSMap 325 with a 200 KHz transductor; the seismo-acoustic survey was performed using a sidescan sonar with integrated GPS, fish tow, and 200 KHz transductor; the magnetometry was performed using a portable cesium-vapor magnetometer with integrated Garmin GPS, and Winglink v. 10 software for data corrections.
The paleoseismic study focused on the western segment of the LATF and consisted of a microtopographic survey performed using a Leica precision GPS device and the analysis of three trenches. The microtopography produced a topographic map with 20-cm-interval contour lines used to choose the trench sites. The study of the trenches was performed according to the methodological guidelines detailed in McCalpin (1996).
The first trench or T1 (19°48'25.2" N, 101°27'03.6" W, 2164 m a.s.l.) is a NE-SW-oriented dig in the LAG crater, that exposes the phreato-magmatic sequence deposited in the maar and in the related basement formed by the andesitic lava flow of Cerro Pelón volcano (Figure 3). Trenches T2 and T3 were excavated in the fault scarp and were north-south-oriented, but only T2 (19°48'14.6" N, 101°27'50.4" W, 2118 m a.s.l.) allowed to observe the fault and affected soils. Details oftrenches T1 and T2 will be discussed in the results section.
Figure 3 Left: Simplified stratigraphic sequence of the volcanic deposits in the area of the La Alberca de Guadalupe (LAG) maar (Not to scale). The inferred sequence and the dates are by Kshirsagar et al. (2015). Right: Details of the faulted phreatomagmatic sequence and photographic record of its structures: folding, andesitic lithic and impact sag, small scale graben-horst geometry, cross-stratification, accretionary lapilli with diameters between 0.2 and 1.5 cm, micro-faulting and ash deformation, and fault affecting material of varied granulometry (thick lapilli to fine ashes). Abbreviations: macrA = masive layer of ash with accretionary lapilli; sA = Stratified, fine ash layer; mLA = massive, fine ash layer.
A pedological analysis was carried out at the Laboratorio de Suelos of the Universidad Michoacana de San Nicolás de Hidalgo since the LATF displaces recent soils. Three sites were chosen for soil profiling: one inside trench T2, and the other two 10 m away from the trench, in the hanging wall and in the footwall, respectively. The materials were characterized physically and chemically, and classified according to the Food and Agriculture Organization (FAO) guidelines for genetic soil horizon interpretation and diagnosis (Jahn et al., 2006). Only three samples from sediments and soil material exposed in T2 were selected and dated since most of the area lacks of datable material. The samples ages were obtained with accelerator mass spectrometry (AMS) radiocarbon by BETA Analytic Laboratories. The ages were calibrated with the IntCal13 curves (Table 1).
Table 1 Radiocarbon ages measured for organic sediments from trench T2 following AMS technique. Results were provided by Beta Analytic Inc. The calibrated ages have been obtained using 2σ uncertainty, and the IntCal13 curves (Reimer et al., 2013).
| Sample name (unit) | Type of sample | Measured age (years BP) |
C13/C12 | Conventional age | Calibrated age (years BC/BP) |
|---|---|---|---|---|---|
| C1 | Soil (bulk) | 19160 ± 70 | -17.2 0/00 | 19290 ± 70 BP | 21485-21085 BC |
| (Unit 4c) | 23435-23035 BP | ||||
|
| |||||
| S1 | Soil (bulk) | 17310 ± 60 | -16.6 0/00 | 17450 ± 60 BP | 19290-18950 BC |
| (Primary soil: Unit 4b | 21240-20900 BP | ||||
| filling A1 fault) | |||||
|
| |||||
| A1 | Soil (bulk) | 10480 ± 40 | -18.5 0/00 | 10590 ± 40 BP | 10710-10575 BC |
| (Secundary soil: U3) | 12660-12525 BP | ||||
|
| |||||
| 10515-10485 BC | |||||
| 12465-12435 BP | |||||
RESULTS AND DISCUSSIONS
Geomorphology
The morphological analysis of the study area shows volcano-tectonic and seismotectonic landforms that match with the Dramis and Blumetti classification (2005). The formation of straight and steep scarps, counter-slopes in El Picacho volcano and small graben-horst features were related to tectonic stress (seismic/volcanic). According to Blumetti et al. (2002), these landforms are indicators of extensional environments. Other landforms in the area such as deviation of the drainage network, decapitated channels, alignment of monogenetic cones, and en échelon faults indicate a minor left-lateral displacement, according to Keller and Pinter (1996).
Similar seismotectonic landforms are reported in other areas of the MAFS (Ramírez-Herrera, 1998), such as the displacement of the southern flank of the Amealco Caldera (Aguirre-Díaz and Mc Dowell, 2000), the graben in the Temascalcingo volcano (Roldán-Quintana et al., 2011; Sunye-Puchol et al., 2015), the en échelon geometry of Los Azufres faults, the collapse of El Estribo volcano (Pola et al., 2014) and the uplift of lacustrine sequences in Pátzcuaro lake (Israde-Alcántara et al., 2005; Garduño-Monroy et al., 2009).
Bathymetry and magnetometry
The bathymetry shows a maximum depth of 11 m in the lake of the LAG maar; the profiles allow to trace the LATF underwater and visualize its vertical displacement on the lake floor (Figure 4 profile A, B). The vertical displacement was quantified from the sonograms, which record a vertical displacement of ~50 cm on the fault trace (Figure 4 profiles C, D). These profiles also differentiate four sedimentary units (U1-U4, Figure 4), and allow to delimit the diatreme geometry. Finally, the magnetometric survey shows a N84°oriented diatreme (volcanic neck), matching the superficial LATF trace (Figure 4).
Figure 4 a) Bathymetric map and localization of acoustic profiles through the lake of the LAG maar. Profiles A and B record displacement in the lacustrine floor, the arrows indicate anomalies. b) Sonograms of profiles C and D showing different sedimentary units and displacements in the lacustrine floor. c) Topographic map combined with the reduction-to-the-pole regional magnetic data map (colors). The dotted line indicates the La Alberca - Teremendo fault trace.
Structural analysis
Three fault systems with ENE-WSW, NW-SE and NE-SW trends are found in the study area. The interaction between those systems creates conspicuous en échelon geometries and pull-apart structures. We identified eight regional faults with ENE-WSW trend, belonging to the MAFS. All of the faults develop scarps with heights between 20 and 100 m, but only 25 % of them cause deviations in the drainage network (Supplementary Table 1, Figure 2).
In this study, a minor, left-lateral displacement of the ENE-WSW-trending faults was defined, with a striation measurement at the LATF plane (N80°E/76°NW dip, 60° rake) and the geomorphological features previously described. According to Pérez-Orozco et al. (2018), in this area, the ENE-WSW-trending faults present normal displacements with a minor left-lateral component, as the result of a NW-SE extension. In Mennella (2011), the NE-SW and NW-SE-trending faults exhibit normal displacements with minor left-, and right-lateral components, and seem to be reactivated by the movement of the ENE-WSW-trending faults.
The stress distribution in the zone shows that σ3 is oriented NW-SE, σ2 is oriented in NE-SW direction, and σ1 is in vertical position, slightly displaced towards the horizontal (Figure 2b). Similar stress distribution was found in the Acambay-Tixmadejé fault (Martínez-Reyes and Nieto-Samaniego, 1990; García-Palomo et al., 2000; Suter et al., 2001; Quintero-Legorreta, 2002).
Paleoseismic analysis
This section details the characteristics of the first (T1) and second (T2) trenches, where the faulting seems better preserved. Both sites expose the phreatomagmatic sequence deposited during the emplacement of the LAG maar. The sequence is composed of an alternation of four types of deposits: 1) a stratified ash fall layer, 2) a fine ash layer with formation of accretionary lapilli, 3) a discontinuous, wedged gravel layer, and 4) a massive, poorly sorted, grain-supported deposit, composed of sub-angular boulder-to gravel-size lithics. Ballistic impact structures and crossed stratification structures were also observed (Figure 3).
T1 trench
T1 (21 m long and 3 m high) is located in the eastern flank of the LAG crater (Figure 2c), and exposes the phreatomagmatic sequence affected by normal faults and conjugated fractures with NNW-SSE, NW-SE and E-W trends (Figure 5). These faults were classified in main faults (first-order structures) and small-scale faults (second-, third- and fourth-order structures) according to their size and the number of stratigraphic units affected (Supplementary Table 2). T1 shows a graben-horst-graben geometry, limited by the first-order structures, identified as F1, F2, F3, F4, F5, and F6 (Figure 5b, 5c). F1 (N67°W-53°SW) is the main fault exposed in the trench; we measured the vertical displacement on its plane. However, the total displacement in the trench could not be calculated because of the lack of correlative beds across the outcrop and the geometry of the wedged strata.
Figure 5 a) Composed photography of the trench T1, an artificial outcrop exposed along the path that connects the crater border with the lake, inside the LAG maar. In the photo, the main faults are delineated. b) Simplified diagram of T1 with the main faults identified (first-, second-, third- and fourth-order). c) Rose diagram of the measured faults indicating a preferential NW-SE strike. d) Sketch map showing the orientation of these outcrop-scale faults with respect to the La Alberca-Teremendo fault and their interpretation as Riedel faults.
The first graben (graben a) is delimited by F1, F2 (N23°W-83°SW), and F3 (N10°W-84°SW). In F1 we measured two displacements: at the top, a 92 cm-displacement of an alternated sequence of coarse ash layers with accretionary lapilli and gravel strata; at the bottom, a 107 cm-displacement of a lithified, fine ash layer (green and yellow layers in Figure 5a, respectively). The differential displacement indicates that, at least, two rupture events occurred on the F1 fault plane. F1 and the space between F2 and F3 are filled with a mix of phreatomagmatic materials, soil, and imbricated rocks with diameters between 10 and 20 cm.
The horst is delimited by F3 and F4 (N83°W-62°SW). The latter is filled with a mix of phreatomagmatic materials and soil; the vertical displacement measured on F4 was 123 cm of a layer of fine ash with accretionary lapilli (brown layer in Figure 5a). Finally, the second graben (graben c) is delimited by F4, F5 (N08°W-75°E), and F6 (N39°W-61°NE); F5 is filled with phreatomagmatic materials mixed with soil; the displacement measured on F6 was 20 cm of a gravel layer overlaid by lithified ash strata (purple layer in Figure 5a).
The small-scale faults (second-, third-, and fourth-order structures) are more abundant in the horst, followed by graben c, and less abundant in graben a; they create small-scale graben, semi-graben and horst geometries, and also merge forming faults nets. These structures do not show vertical displacement across the grabens, and the third-, and fourth-order structures disappear, both, upwards and downwards in the sequence. These characteristics suggest accommodation of materials due to horizontal extension of non-tectonic nature; accordingly, they cannot be used in the paleoseismic estimations.
T2 Trench
T2 trench (10 m long and 4.5 m deep) is located at the western end of the LATF (Figure 2c and Figure 6). The trench exposes the phreatomagmatic sequence of the LAG maar covered by two dark brown, superposed soils with thicknesses between 0.5 and 1.5 m. The soils characteristics are described in detail in the next section.
Figure 6 a) Microtopography of the La Alberca-Teremendo fault scarp at the western flank of the LAG maar crater (Equidistance of contour lines: 1 m), and location of trenches T2 and T3. b) Photography of the scarp at the T2 location. c) Photography of the contact between the phreatomagmatic sequence and the soil along the LATF plane. d) Scheme of the photography c). e) Opening of trench T2 with the fault plane exposed.
T2 shows the LATF trace and three smaller antithetic faults (A1-2-3, Figure 7), as well as second- and third-order faults and associated fractures, which displace the phreatomagmatic sequence. The LATF is a normal fault trending N260° and 76°NW dip; at the bottom, it changes to vertical, and at the top has an 89° dip, and a 60° rake. It has a total vertical displacement of335 cm and a net vertical displacement of262 cm. The antithetic faults that limit the blocks were identified as A1 (N60°E-80°SE), A2 (N80°E-60°SE) and A3 (N85°W-80°SW). These faults accumulate a vertical displacement of 73 cm (10 cm for A1, 57 cm for A2, and 6 cm for A3). The hanging wall shows a domino-style, normal faulting, integrated by three blocks which rotate 8° to 15° away from the foot wall; the spaces between blocks (4 to 61 cm) are filled with a mix of phreatomagmatic materials and soil (Figure 7).
Figure 7 Log of trench T2, rose diagram of the fault planes exposed, directions of the stress field, and pedological analysis and classification of the soil affected by the LATF. Abbreviations: A1-3 = Antithetic faults; s1-5 = second order faults; LATF = La Alberca-Teremendo fault.
Use of soil analysis as a tool for the identification of paleoearthquakes
The soils analyzed in the study area (hanging wall and footwall) were classified as a Luvisol with well-developed A and B horizons, formed from the phreatomagmatic materials with andesitic composition. The on-fault exposed soil (T2 trench) shows variations in its physicochemical properties: the bottom was classified as Luvisol, but at the top, the vertic properties become important, suggesting an evolution towards a Vertisol (Supplementary Table 3). The main variations in the soil exposed on T2 trench are:
A noticeable textural change caused by the increase in percentage of clay at 100 cm- depth (terminus Abruptic, Jahn et al., 2006).
The presence of a lithological discontinuity at 80 cm-depth caused by an anomalous increase in the percentage of organic material (terminus Ruptic, Jahn et al., 2006).
The presence of recent sediments, 50-cm thick, deposited over an old soil (terminus Novic, Jahn et al., 2006).
The previous characteristics are not related to the andic properties of the soil, but rather suggest a disordered soil, with high rates of skeleton and expansive clay.
The soil profile exposed in T2 can be interpreted as two superimposed soils. At the base is an older, primary soil composed by two horizons 2C, characterized by vertic properties, a clayey texture and the presence of slickensides. This soil lacks horizons A and B, possibly due to increased erosion (Figure 7). According to the radiocarbon dates, the lower limit of this soil has an age of 21240-20900 cal yr. BP.
The primary soil is covered by a secondary soil, made of horizons A1-2, B and C. A1-2 have humified organic matter mixed with altered minerals; horizon B is characterized by argillic properties, accumulation of clay, high expansibility, and slickensides; finally, horizon C has vertic properties, slickensides, and a sandy texture (Figure 7). The radiocarbon age at the base of the soil is 12660-12435 cal yr. BP, and continues its development until today.
The differences between the superimposed soils that evolved at different times, suggest a disruption in the pedogenic evolution of the buried soil. This fact may be related to the vertical displacement of the fault and the subsequent formation of a surface soil with unique characteristics derived of the previous soil as parental material. This new soil was displaced again in recent times. The soils analysis suggests that, at least, two ancient earthquakes in the LATF occurred after the emplacement of the LAG maar.
Number of seismic events
The net vertical displacement measured for the LATF in the fault plane in T2 was of 262 cm (335 cm of total synthetic displacement minus 73 cm of total antithetic displacement). We suggest that the observed displacement could be the result of three events: the first could be a volcano-tectonic event related to the LAG maar emplacement; the second event was identified from a net vertical displacement of 59 cm of the buried soil; finally, the third earthquake was identified from a net vertical displacement of 115 cm of the surface soil that likely occurred at prehistoric times. It is clear that the surface soil is displaced and the scarp shows little erosion, but neither historic nor instrumental record exists about this event; therefore, we assigned a maximum age of 1000 yr. BP for this earthquake.
However, estimations should be carefully considered because of the nature of the LATF, and the coexistence of magmatic and tectonic activity in the area.
Deformation model
According to the scenario previously proposed, we defined a model for the reconstruction of the deformation in trench T2 (Figure 8). The first phase is the deposit of the non-deformed, phreatomagmatic sequence during the formation of the La Alberca de Guadalupe maar around 23434-21000 yr. cal BP. The second phase is the occurrence of a volcano-tectonic earthquake related to the formation of the LAG maar around 23434-21000 years ago, causing an 88-cm vertical displacement of the phreatomagmatic sequence and generating a fracture filling of pyroclastic material that resembles a colluvial wedge. This event is followed by a quiescence period (third phase) characterized by the erosion of the upper layers of the phreatomagmatic sequence in the hanging wall (Units 5a-d), and followed by the development of a primary soil (phases fourth and fifth, Units 4a-b, Figure 8). The sixth phase is the occurrence of a tectonic event 12660-12435 years ago, that caused a vertical displacement of 59 cm of the primary soil and the phreatomagmatic sequence. The seventh phase is a quiescence period characterized by the erosion of unit 4a in the hanging wall, followed by the development of a secondary soil from material of the buried soil (eighth phase). Finally, the ninth phase is a third event occurred in recent times, causing a vertical displacement of115 cm of the newly formed soil. The scarp is still visible today.
Abbreviations: A1-3 = Antithetic faults; VD = Net vertical displacement; LAG = La Alberca de Guadalupe.
Figure 8 Retrodeformation model for the activity of La Alberca-Teremendo fault (LATF) based on the record from trench T2, indicating the main events (earthquakes) and the current aspect of the trench site.
Slip Rates (SR)
To measure the slip rate (SR) we followed the methodology used by McCalpin (1996) that considers:
The average slip rate was estimated for the net vertical displacement of the LATF (262 cm) measured in T2 and for the displacements of F1 measured in T1 (92 and 107 cm). The slip rates for individual events were estimated but considered with caution because of the lack of datable material in the area (Table 2).
Table 2 Slip rates (SR) obtained for partial displacements measured for NW-SE faults in trench T1 and for total and partial displacements of LATF measured in trench T2. Recurrence intervals (RI) obtained for the measured displacements. Ev = Displacement event. (*) = age of the LAG maar from Kshirsagar et al., 2015.
| T1 Trench: F1 | T2 Trench: LATF | Retrodeformation scenario | |||||
|---|---|---|---|---|---|---|---|
| Ev 1 | Ev2 | Complete fault | Ev 1 | Ev 2 | Ev 3 | ||
| Total displacement (mm) | 920 | 1070 | 3350 | 1130 | 820 | 1400 | |
| Net displacement (mm) | — | — | 2620 | 880 | 590 | 1150 | |
| Average vertical displacement (mm) | — | — | Total | 1117 | |||
| Net | 873 | ||||||
|
| |||||||
| Time elapsed (years) | *21000 | 23434-23035 | 23434-23035 | 10999 - 10,375 | 12660 - 12,435 | ||
| SR (mm/year) | 0.046 | Total | 0.142 - 0.145 | 0.048 - 0.049 | 0.074 - 0.079 | 0.11 - 0.112 | |
| Net | 0.112 - 0.114 | 0.037 - 0.038 | 0.053 - 0.056 | 0.09 - 0.092 | |||
| RI (years) | ~11630 | Total | 7784 ± 81 | ||||
| Net | 7726 ± 68 | ||||||
If we assumed a purely tectonic behavior, the average SR obtained for the net vertical displacement of the LATF was 0.114 mm/year, this value represents 23 % of the SR estimated by Suter et al. (2001) using geomorphological indicators (0.5 mm/year). The slip rates of individual events show great variation (from 0.037 to 0.141 mm/year) indicating a high degree of uncertainty, consequently, they cannot be used in the evaluation of seismic cycles. The slip rate obtained for the vertical displacements measured in F1 was 0.046 mm/year, but F1 is a secondary fault that could have moved with the volcano-tectonic events of the LATF.
The slip rate values for the LATF are similar to the SR obtained for other faults in the MAFS, especially in the Morelia-Cuitzeo region (0.05 to 0.18 mm/year) (Suter et al., 2001); whereas in the Acambay graben and the Pátzcuaro region, the SR is more variable: 0.03-0.37 mm/year and 0.009-2.78 mm/year, respectively (Suter et al., 1992; Langridge et al., 2000; Garduño-Monroy et al., 2009).
Recurrence intervals
According to McCalpin (1996), the recurrence interval can be assessed with two methods: the direct method and the geological method. The first one produces an average recurrence interval; the second method produces individual recurrence intervals, better reflecting the dynamics of individual faults and their seismic cycles (Scholz, 1990).
For the LATF, we used the direct method to estimate the average recurrence interval for the average vertical displacement considering three displacement events (total = 111.7 and net = 87.3 cm), with the following relation:
Where VDaver is the average vertical displacement (mm); SR is coseismic slip rate (mm/yr) and C is slip rate per creep. Recurrence intervals for individual events in T2 trench cannot be calculated because of the uncertainty of the associated slip rates. For all the calculations, we assumed that in the LATF, C = 0.
We obtained average recurrence intervals of 7784 ± 81 years for total displacement and 7726 ± 68 for net displacement for the LATF in T2; and an average of 11630 years recurrence interval for F1. These values are closer to the recurrence intervals calculated for other segments of the MAFS such as the Pastores fault (10-15 Ka according to Langridge et al., 2013; Ortuño et al., 2015), and the San Mateo fault (11570 ± 5320 years according to Sunye-Puchol et al., 2015). In the faults of the Pátzcuaro-Morelia region the values are more variable: our RI is similar to the Huiramba fault, the La Paloma fault, the Charo fault and the Queréndaro-Indaparapeo fault (4.5 to 20 Ka; Garduño-Monroy et al., 2009). These recurrence intervals are a first approach to evaluate the fault activity and should be taken cautiously, as occurrence of strong earthquakes could be clustered in time, with few years and even days of separation, as has been observed elsewhere in the world (Iezzi et al., 2018). We must also consider the frequency of events in the LATF that is influenced by the magmatic activity in the area, which likely changes the fault stresses.
Probable magnitude of individual earthquakes and maximum potential magnitude expected for a seismic rupture in the LATF
The probable magnitude for fault ruptures can be assessed using empirical relations that use fault parameters such as: the surface rupture length, vertical displacement per event, slip rate, and recurrence intervals. Published empirical relations were constructed with global data bases and can be used as reference for sites with scarce or null information. The most used empirical relationships apply the surface rupture length (SRL), but its main uncertainty resides in the lack of information about the geometry and prehistoric length of the faults (Lafuente-Tomás, 2011).
The LATF has characteristics that indicate it could move under the influence of both, an extensional volcanic environment and the actual tectonic regime. So, in order to obtain a reliable evaluation of the deformation for this fault, a specific paleoseismic relation should be constructed; however, more information is needed from others structures in this area in order to achieve a reliable relation. We used published relations for the LATF, although we took into consideration that they could produce an overestimation of magnitudes of the paleoseismic events, and magnify the fault seismic potential. We applied relations that include the average vertical displacement (VD) to assess the probable magnitude associated to individual events, and relations that use SRL and SR to estimate the maximum magnitude for an earthquake that displaces the entire fault trace; finally, we made a comparison between relations to establish the best estimation (Supplementary Table 4).
For the LATF (Table 3), if we assume that all the deformation is of tectonic nature and that the segments are connected in depth, the magnitude Mw values based on the total surface rupture length of 26.1 km, ranges between 6.7 and 7.3, and between 5.9 and 6.7 for superficial magnitude (MS). These values are slightly higher than those obtained for other segments of the MAFS. For example, for the Pastores, Acambay-Tixmadejé, Maravatío, La Paloma, Venta de Bravo, San Mateo faults, values range between 6.4 and 7.0, and for the Pátzcuaro-Morelia faults, the Mw values are 5.8 to 7.1. Meanwhile, for the individual events, the probable magnitude Mw, based on the average slip, ranges between 6.6 and 7.0. Finally, for the slip measured in F1 of T1, the estimated magnitudes range Mw = 6.5-6.6. The similar displacements values of trenches T1 and T2 suggest that the inferred events on the LATF trace influence the fault F1 movement, and the faulting could be simultaneous. The Mw figures obtained using the relations of Anderson et al., (1996) and Stirling et al., (2002) are higher than those obtained using the relations of Wells and Coppersmith (1994) and Wesnousky (2008), but the minimum magnitude remains around 6.7 for a rupture of 26 km of the LATF.
Table 3 Resolved scaling relations for assessment of potential Mw magnitude in the displacement events identified in T1 and T2 trenches. The values obtained for the Mohammadioun and Serva (2001) relation are referred as surface wave magnitude (Ms). VD = Vertical displacement. SRL = Surface rupture length. Data marked with (*) were obtained from Clemente-Chávez et al. (2013) and Singh et al. (2011, 2012).
| Event | Used relations | Assessed magnitudes (Mw) | ||||
|---|---|---|---|---|---|---|
| Maximum rupture event for the La Alberca - Teremendo fault (SLR) |
Wells and Coppersmith, (1994) | 6.72 | ||||
| Anderson et al., (1996) | 6.77-7.29 | |||||
| Stirling et al., (2002) | 6.70-7.32 | |||||
| Wesnousky (2008) | 6.78 | |||||
| Mohammadioun and Serva (2001) | 5.88 (Ms) | using Δσ = 11.1 bar (*average stress drop for the TMVB) |
||||
| 6.69 (Ms) | using Δσ = 50 bar (*minimum stress drop for Morelia region) |
|||||
| T1 | T2 | |||||
| F1 Average VD | Average VD | Events | ||||
| Ev 1 | Ev 2 | Ev 3 | ||||
| Individual net displacements measured for T1 and T2 trenches (VD) |
Hanks y Kanamori (1979) | 6.5 | 6.6 - 6.9 | 6.7 | 6.6 | 6.8 |
| Wells and Coppersmith (1994) | 6.6 | 6.8 - 7.0 | 6.8 | 6.6 | 7.0 | |
If we assume a purely tectonic activity, the magnitudes predicted for the LATF are greater than expected, suggesting the possibility of multi-fault ruptures (primary and secondary displacements triggered by other fault ruptures in the area), in a similar scenario to that of the Acambay graben during the 1912 earthquake (Arzate, et al., 2018; Ortuño et al., 2019). We must consider that secondary behavior on a fault system can also be triggered by large volcanic eruptions (as mentioned in Ortuño et al., 2019). For the LATF in particular both triggers are equally probable.
In a different scenario, if we assume that the movement of the LATF is the result of volcano-tectonic events produced by dike-induced structures under an extensional environment, the estimated magnitudes Mw should be in a range of 4.0 to 5.5 (exceptionally up to 6.5, Hackett et al., 1996), as has been observed in others areas of the world (e.g. Africa, Iceland and the USA, Hackett et al., 1996). Moderate magnitude events could be the most probable scenario for a displacement of the LATF or the NW-SE associated faults. However, moderate magnitudes events are not harmless; for example, in the Chapala graben, on October 2, 1847, an earthquake (MI 5.7) caused several human losses and serious damages to the villages of Poncitlán and Ocotlán (Suter, 2017).
Seismic hazard
The favorable orientation of the LATF in relation with the current stress field in central México enhances the probability of a future rupture in the fault. Considering the minimum potential magnitude for a tectonic event on the fault (Mw 5.9), the earthquake would be perceived in a radius of ~ 200 km, where some of the largest cities in the country concentrate almost 43 % of the population. For a volcano-tectonic earthquake (Mw 4-5) the area would be reduced, but it could affect parts of the states of Michoacán, Guanajuato and México (Figure 9a).
Figure 9 a) Perception area for a moderate event with magnitude Mw 5.9 and epicenter in the La Alberca-Teremendo fault. b) Areas prone to develop environmental effects during an earthquake with intensities between VIII - X (1000 km2-5000 km2). Numbers indicate Mexico states: 1 = Michoacán; 2 = Estado de México; 3 = México City; 4 = Querétaro; 5 = Guanajuato; 6 = Jalisco.
According to the Environmental Seismic Intensity scale (ESI-2007, Michetti et al., 2007), the LATF could generate events with intensities from IV to X, characterized by easily observable environmental effects, often permanent and diagnostic, starting from intensity VII. The analysis of environmental effects is recommended in sparsely populated areas, but they are less suitable for intensity assessment than the effects on humans and manmade structures. For such events, the damaged areas could be as large as 1000 km2 (IX) and 5000 km2 (X) (Figure 9b). Permanent and temporary effects could appear within these areas. Primary effects could be surface ruptures of kilometric scale, and local uplift and subsidence. Secondary effects could be ground cracks of metric scale, slope movements with bulk volumes of 103-106 m3, liquefaction processes, anomalous waves with over 1 m run-up in close basins, temporary spring drying, and streams changes. All these alterations would be particularly important in cities north of Michoacán, such as Morelia, Pátzcuaro, Zacapu, and Cuitzeo, where the granular basement modifies the accelerations of the terrain and can generate liquefaction and infrastructural damage.
CONCLUSIONS
The La Alberca-Teremendo fault (LATF) is a 26 km long, complex, en échelon structure, which has normal displacement with a minor left-lateral component. The activity of the LATF is closely linked to the emplacement of the La Alberca de Guadalupe maar and triggered the collapse of the Picacho volcano summit. These volcano-tectonic relationships reflect recurrent fault activity during the Quaternary.
The LATF activity was evaluated in the volcanic phreatomagmatic sequences of the LAG maar. The slip rate of this fault was estimated at 0.114 mm/year with an average recurrence interval of 7726 ± 68 years and an average net vertical displacement of 87.3 cm. Table 4 shows a compilation of the measured parameters in this study.
Table 4 Summary of parameters for the paleoseismic analysis in trenches T1 (NW-SE faults) and T2 (La Alberca - Teremendo fault).
| Parameter | |
|---|---|
| Thickness of seismogenic layer | 12-20 km (Ortuño et al., 2015) |
|
| |
| Crust rigidity | 3 x 1011 dyn/cm2 |
|
| |
| Maximal surface rupture length (SRLmax) for LATF | 26.176 km |
|
| |
| Rupture area (A) for LATF | 8.22 km2 |
|
| |
| Accumulated maximum throw for LATF | 50 m |
|
| |
| Partial vertical displacement of F1 fault in T1 (VDT1) | 92 cm (top) 107 cm (bottom) |
|
| |
| Total vertical displacement of LATF in T2 (VDT) | 335 cm |
|
| |
| Net vertical displacement of LATF in T2 (VDN) | 262 cm |
|
| |
| Total average displacement of LATF in T2 (VDaverT) | 111.7 cm |
|
| |
| Net average displacement of LATF in T2 (VDaverN) | 87.3 cm |
|
| |
| Individual displacement event 1 | 113 cm (Total) 88 cm (Net) |
|
| |
| Individual displacement event 2 | 82 cm (Total) 59 cm (Net) |
|
| |
| Individual displacement event 3 | 140 cm (Total) 115 cm (Net) |
|
| |
| Elapsed time for displacement of F1 fault in T1 | 21000 years (Kshirsagar et al., 2015) |
|
| |
| Elapsed time for accumulated displacement of LATF in T2 | 23035-23434 years |
|
| |
| Slip rate for F1 fault in T1, (SRf1) | 0.046 mm/yr. |
|
| |
| Average slip rate for VD (SRaverT) in T2 | 0.142-0.145 mm/yr. (Total) 0.112-0.114 mm/yr. (Net) |
|
| |
| Recurrence interval obtained for F1 fault in T1 | ~11630 years |
|
| |
| Recurrence interval obtained for LATF in T2 | 7784 ± 81 years (Total) 7726 ± 68 years (Net) |
Paleoseismic analysis suggested that the fault has moved three times in the past 23000 years: the first was a volcano-tectonic event that could be cogenetic to the LAG maar emplacement; we assume a moderate to low magnitude for this event. The other two events were observed more clearly and were related to earthquakes with magnitudes Mw of 6.6-7 if we assumed that they were the result only of tectonic stress.
The two earthquakes were identified through the superposition of two soils observed during the pedological analysis. This study revealed a secondary soil developed from an eroded buried soil, whose evolution was interrupted by the fault activity. The soil analysis provides crucial information for the assessment of the seismogenic potential of the LATF and could be helpful in similar scenarios where faults affect recent soils.
Although the LATF is a complex structure within the MAFS, its seismogenic potential is important because the fault is favorably oriented in relation to the current stress field in central Mexico. Empirical relationships suggest that the fault may generate earthquakes with moderate magnitude (4-5) during a volcano-tectonic event, or even with large magnitude (Mw 6-7) during a purely tectonic event.
In Mexico, historical intraplate earthquakes with similar magnitude to that estimated for the LATF have caused severe damages across vast areas, e.g. Chapala Graben 1847, MI 5.7; Acambay 1912, Mw 7; and Axochiapan 2017, Mw 7.1. The effects of these earthquakes were comparable to those produced by the largest subduction seisms in the country, e.g. Jalisco 1932, Mw 8.1-8.4; Lázaro Cárdenas 1985, Mw 8 (Singh and Suárez, 1987), and Colima 1995, Mw 8, among others.
Several fault segments with characteristics similar to those of the LATF exist in the western sector of the MAFS, and remain unstudied (for example, Pajacuarán or Ixtlán faults). The proximity of the LATF to some of the largest cities in the country turns it into a seismic source that should be considered in hazard assessment studies and urban planning of Michoacán and central Mexico.
