Characterization of seismogenic crustal faults in the Gulf of Guayaquil , Ecuador

Few moderate-to-strong earthquakes associated with active and capable geological faults have been documented for the southern coastal region of Ecuador. The seismic record of Ecuador initiates with the Guayaquil earthquake in 1787 (Mw 6.5), followed by the 1943 (Mw 6.2), and the most recent one in 1980 (Mw 6.1). The available data is insufficient to evaluate the seismic hazards associated with faults capable of generating seismic events of magnitude Mw≥6.0 in the region. Also, earthquakes of minor magnitudes can be disregarded as they do not induce significant ground coseismic effects. In this context, this study presents a catalog of geological faults, delineating 40 segments of capable and active faults on the sea floor of the Gulf of Guayaquil and inland segments of Guayas, Santa Elena, and El Oro provinces. This methodological approach estimates a variety of seismicity levels ranging between Mw 6.2 and Mw 7.2, with rock peak ground accelerations between 0.24 g and 0.41 g. These values have been obtained from empirical regression equations applied to the length of capable geological faults. The F-40 seismogenic structure located in the accretionary prism, close to the subduction zone, is capable of generating Mw 8.2 earthquakes and potentially causing coseismic ground damage to the city of Guayaquil located ca. 177 km to the NE. Furthermore, local tsunami hazards may affect severely areas that are densely populated, with developing industrial areas, on the coast of the Gulf of Guayaquil. This structural geological analysis provides useful new data for seismic hazard assessment.

The main seismogenic structure of the subduction zone is located between 50 km and 140 km west of the limits of the Gulf of Guayaquil (Fig. 1).This structure is capable of generating earthquakes of Mw 8.0 and higher, as documented in the past (Kanamori, 1977;Abe, 1979;Mendoza and Dewey, 1984;White et al., 2003).The Gulf of Guayaquil itself represents a zone with a seismic gap (Contreras, 2013), and therefore, it is feasible to hypothesize that this area probably stores seismic energy, enough to generate an earthquake of magnitude comprise between Mw 6 and Mw 8 (Ioualalen et al., 2014).In this context, this study aims at recognizing and evaluating the seismic potential of crustal faults situated in the offshore of the Gulf of Guayaquil.
The methodology considers the estimation of the lengths of the known active geological faults that can generate earthquakes with magnitudes Mw≥6.0, based on works such as those from Well and Coppersmith (1994) and Michetti et al. (2007).Smaller geological faults which may be able to generate earthquakes with a magnitude below Mw 5.5 have not been considered in this study since these earthquakes cause smaller side-effects or coseismic geological effects.As historical earthquakes related with crustal faults are poorly recorded on the mainland at the southern coast of Ecuador (Fig. 1), this methodology will provide a better understanding about the maximum magnitudes of future potential earthquakes, and their associated PGA values.Also, a resulting data set may contribute to determining the seismic hazard level for densely populated and developing industrial areas on the southern coast of Ecuador.

Introduction
Regression equations applied to kinematic and geometric parameters of both active and capable geological faults have been used in the estimation of the maximum magnitude of earthquakes (Wells and Coppermisth, 1994;Wesnousky, 2008;Leonard, 2010;Stirling et al., 2013), and peak ground acceleration on rocks, PGA (Fukushima and Tanaka, 1990).Empirical relationships between isoseismal area and major axis radius of intensities and values of magnitudes have also been used to determine the seismic epicenter (Marín et al., 2008).Subsequently, this kind of methodology indeed provides estimated values of the levels of seismicity, in particular in areas where the previous seismic history is reduced, or even absent.
The Gulf of Guayaquil in Ecuador is situated on the eastern side of the Pacific Ring of Fire.This is part of an active continental margin where the oceanic Nazca Plate subducts below the South American plate approximately 47 mm/yr (Barazangi and Isacks, 1976;De Mets et al., 1990;Trenkamp et al., 2002;Eguez et al., 2003;Veloza et al., 2012).The Nazca Plate includes the aseismic Carnegie Ridge, associated with the Galapagos hot spot magmatism (Shumway, 1954;Lonsdale, 1978;Freymuller et al., 1993;Jaillard et al., 1995;Toulkeridis, 2011).The deformation of continental crust seems to be influenced also by the collision, and mechanical coupling, of the Carnegie Ridge against the South American plate (Gutscher et al., 1999;Dumont et al., 2014).Deformation is partitioned in an NNE-SSW trending strike-slip and reverse faults (Hughes and Pilatasig, 2002).The continental domain is formed of two entities: the Caribbean, and the South-American continental plates.The contact among them corresponds to a deformation zone named Guayaquil-Caracas Mega-Fault, which extends from the Gulf of Guayaquil in Ecuador until Venezuela (Kellogg and Vega, 1995; de Guayas, Santa Elena, y El Oro.Este procedimiento metodológico ha permitido estimar niveles de sismicidad con magnitudes que varían entre Mw 6,2 a Mw 7,2, y aceleraciones en roca, PGA, en el orden de los 0,24 g a los 0,41 g.Estos valores fueron obtenidos con ecuaciones de ajuste aplicadas a fallas geológicas capaces.La estructura sismogénica F-40 localizada en el prisma de acreción (junto a la zona de subducción) es capaz de generar un sismo de magnitud aproximada a Mw 8.2, pudiendo causar daños cosísmicos a la ciudad de Guayaquil, ubicada a aproximadamente a 177 km de distancia hacia el NE, y también generar tsunamis sobre áreas densamente pobladas, y sectores industriales, asentados en el borde costero del Golfo de Guayaquil.Toda esta información geológica estructural proporciona nuevos datos para considerar en los estudios de peligro sísmico.

Seismic history of Ecuador
From 1653 until May 2016, the territory of Ecuador recorded 28 earthquakes related with crustal faults and the subduction process, with magnitudes between Mw 5.5 to Mw 7.8 and depths of hypocenters shallower than 20 km (Table 1).These data document a short record of earthquakes in the Gulf of Guayaquil (Chunga and Quiñonez, 2013).The first earthquake was recorded on July 9, 1653, with a magnitude of Mw 5.7 as converted from macroseismic values VI (CERESIS, 19851 ).On July 11, 1787, the first macroseismic earthquake hit the coastal sector of Ecuador with a possible epicenter in the Guayas province (MHQ, 1879).This event generated building damages and significant coseismic effects on the ground, mainly in the city of Guayaquil, while in nearby towns within a 50 km radius such damages were not observed (Chunga et al., 2005a).This suggests a shallow hypocenter with the epicenter in or near the city of Guayaquil.Both events (1653 and 1787) are probably associated with crustal  Moderate earthquakes higher than 5 in the Gulf of Guayaquil and in the peninsula of Santa Elena as well as Guayas provinces are recorded from the local network of the IGEPN.Gray shaded boxes represent the earthquakes near the Gulf of Guayaquil that originated tsunamis with wave heights run-up between 0.2 and 1.8 m.In the event of 1901 the wave heights reached above 2 m (CERESIS).
faults, while the following documented registered seismic events are classified as subduction-related earthquakes.On January 7, 1901, an earthquake with a magnitude of Mw 7.8 occurred in front of the Santa Elena Peninsula (CERESIS, 1985).The degree of the macroseismic intensities and the height of a tsunami run-up in the coastal plains of Salinas and La Libertad remained undetermined.The lack of such data did not allow to localize the epicentral area and the corresponding seismogenic features.
Another seismic event with a magnitude of about Mw 6.9, probably linked with the same subduction interface seismic source, occurred on October 2, 1933.Tsunamis waves, as high as 2-2.5 m, flooded without generating turbulent waves in the coastal edge of the Santa Elena peninsula (Silgado, 1957;Espinoza, 1992;Arreaga and Ortiz, 2002).In 1953, between the Ecuadorian-Peruvian border, an earthquake with a magnitude of Mw 7.8 and an epicenter ca.142 km to the south of Salinas triggered a tsunamis with a wave height of 20 cm that impacted the population of La Libertad allowing rapid flooding on the coastal strip of the city of Salinas (Chunga et al., 2005b).Silgado (1957) estimated for this event a macroseismic intensity of VI for the city of Guayaquil.
The earthquake of February 7, 1959 (Mw 7.2) does not provide relevant ground coseismic effects.Four earthquakes with moderate magnitudes (Mw between 6 and 6.2) were reported in the Gulf of Guayaquil from April 1961 to March 1962 (Table 1).The most recent earthquake being an Mw 6.3 occurred on May 21, 2005 and had a focal length (H) of 39 km (NEIC, 2017).These data suggest that the recurrence of earthquakes with magnitudes greater than Mw 6 is low.
Based on the documented seismic events in the Gulf of Guayaquil, it is necessary to understand what method allows characterizing active or capable faults.This characterization would provide useful insights regarding the seismic potential of this area, based on each of the seismogenic sources, and regarding maximum magnitudes, and Peak Ground Acceleration (PGA).

Ground rock and soil foundation of the Guayaquil area
In Guayaquil area, the features of rocks and soils change locally: a) the urban and commercial centers are founded on clay soils interbedded with unconsolidated sandy to clayey alluvial Holocene levels; b) the southern and southwestern areas are directly lying on muddy estuarine sediments, where the filling layer has a thickness that varies between 1 m and 2 m; c) the north of the city is built on Cretaceous hyaloclastite, sandstones, shales, siltstones, and competent greywacke (Cayo and Guayaquil formations); d) the marginal urban areas (Sergio Toral, Monte Sinai, Trinidad de Dios and Horizontes del Fortín) lie mostly on unstable hills slopes; e) the Kennedy and Urdesa sectors are built along two estuarine branches composed of soft unconsolidated sediments.These described lithological features of rocks and soils certainly can amplify or attenuate seismic waves in the ground during an earthquake of magnitude higher than Mw 6.0 (Chunga et al., 2005a;Aguiar et al., 2016).

Morphology of the Gulf of Guayaquil
The Gulf of Guayaquil is an extensional structure with dextral and sinistral shear faults related with the opening of a fore-arc basin as highlighted by seismic reflection sections (CEPE, 19862 ) and by re-evaluations of offshore lithologic logs (Witt et al., 2006;Cobos and Montenegro, 2010).Cobos and Montenegro (2010) define the Gulf of Guayaquil as a "pull-apart" type basin.Regionally, the primary stress regime is associated with the tectonic thrust of the oceanic plate beneath the Gulf of Guayaquil.This tectonic setting favored the uplift of the Paleocene to Holocene lithologic units in coastal margins (Costa et al., 2006).Submarine morphology of the Gulf of Guayaquil is relatively uneven, ranging from the continental shelf to the continental slope.The multibeam-highresolution sonar, acoustic reflectivity, and acoustics for sediments penetration (about 150 m strata) data illustrate continental deformations related with the 70 km-long active faulting landslide (Megasplay in figure 1) located in the accretionary prism in front of the subduction zone graben (Convemar, 2010).Hazards of submarine mega-landslides have been delineated from bathymetric data, allowing to reconstruct their sedimentary paleo-environment.This reconstruction has been performed assuming that these processes have been potentially active in past glacial periods when the sea level stood 120 meters below the current sea level (Bintanja et al., 2005;Rohling et al., 2009;Elderfield et al., 2012;Chunga and Quiñónez, 2013).Near the graben areas, continental slopes are highly unstable and able to produce submarine landslides, causing re-deposition into depressed areas (Fig. 2).In the intermediate slopes are well-structured and noticeable submarine canyons (e.g., Santa Elena Canyon), where high erosion rates demonstrate the geometry of structures with strong incisions in the seafloor.The Guayaquil canyon, located southward, also indicates considerable circular submarine landslides scarps (Fig. 2).The identification of a thrust fault with a length of approximately 150 km has also been documented in the accretionary prism of the Gulf of Guayaquil (Convemar, 2010; Table 2).

Capable faults and shallow earthquakes
Strong-to-moderate crustal earthquakes in the southern coast of Ecuador have been poorly documented.The historical accounts obtained from the Cuenca and Quito Culture House archives (MHQ, 1879), evidence of the earthquake that hit the city of Guayaquil on June 11 th , 1787.This earthquake has been associated with a seismic source in the Carrizal faults (F-18) or Estero Salado fault (F-17) located 35 and 16 kilometers, respectively, far from the urban area (Table 2).The macroseismic intensity estimated for this earthquake is between VIII to IX degrees of the ESI-07 scale, while the calculated magnitude FIG. 2. Submarine landslides have the potential to generate devastating tsunamis next to the Guayaquil canyon, within the Guayaquil Gulf.ranges between Mw 6 to Mw 6.5 (Chunga, 2010;Chunga et al., 2013).Another Mw 6.2 earthquake, (i.e., Ruffilli, 19483 ) occurred in 1943, in the Guayaquil area.Historical records suggest that this earthquake was accompanied by a strong underground noise related to a massive landslide in a nearby quarry.
The coseismic effects on the ground allowed to assume that the earthquake epicenter was located quite close to the urban area.The subjective perception of the Guayaquil inhabitants is that this earthquake was more violent and shorter than the one occurred on May 13, 1942.The earthquake of August 18, 1980, of magnitude Mw 6.1, had its epicenter located 28 kilometers northwest of the city of Guayaquil (Fig. 3).The macroseismic intensity values varied from VII to VIII degrees on the Mercalli scale (Egred, 1975;Mera, 1999;Chunga et al., 2005a).Multiple houses collapses, as well as telephone and electrical services, were interrupted.The earthquake was felt strongly in the southern, western, and central area of the city.Argudo et al. (1993) described the damages of old buildings where 29 walls of houses collapsed with the partial or destruction of seven of them (Fig. 4).Further evaluations determined that 49 houses needed later demolition.In the historical accounts, the housing damages and the levels of human perception of quakes were described.However, these seismic events certainly generated earthquake environmental effects on the ground that probably were not described.The homes grounded in soft clayey soils in the southwest part of the city may have suffered accelerations of the natural subsidence rate, sinkholes, and soil liquefactions processes.Guayaquil also has hills with unstable slopes where rock falls or landslides may have occurred but remained undocumented.Egred (1975) for the earthquake of August 18, 1980 (Mw 6.1).

Estimation of the maximum magnitude and PGA
Instrumental data of earthquakes with magnitudes higher than Mw 6 and their association with crustal capable faults are poorly documented in the Gulf of Guayaquil.In this context, it is defined as a "capable fault" the type of fault that may be a potential source of future earthquakes.It presents evidence of dislocations or surface displacements occurred in the last 30,000 years, and if the historical seismicity can be linked with a precise fault (IAEA, 2002;Chunga, 2010).Moreover, a fault is considered "potentially active," as a potential source for future earthquakes, if the surface dislocations occurred at least in the last 50,000 years (IAEA, 2002;Robert and Michetti, 2004;Michetti et al., 2007).
To estimate the seismicity associated with a crustal fault in the southern coast of Ecuador, a catalog with 40 faults capable of deforming the ground surface and to generate modest-to-strong earthquakes (6.2<Mw <7.2) was compiled.The database includes faults recognized at the seafloor and on the continental area of the Santa Elena, El Oro, and The Guayas provinces (Table 3).Relevant information about the geometry and kinematics movements of each fault, as well as evidence of vertical displacement from late Pleistocene to the Holocene has been provided by numerous studies (CEPE, 1986;Eguez et al., 2003;Dumont et al., 2005;Witt et al., 2006;Calahorrano et al., 2008;Witt and Bourgois, 2010;Cobos and Montenegro, 2010).Seismic profiles available from EP-Petroecuador were interpreted by CEPE (1986), and they documented that many of the faults do not reach the ground surface because they are covered by late Pleistocene to Holocene sediments.In this study, this kind of faults was classified as "active" and unable.Knowing the fault features allows evaluating possible earthquakes and their features as magnitude, macroseismic intensity, and PGA.
The geometric parameters considered for each selected fault includes: 1) the fault spatial projection on the surface, 2) the fault geometry and kinematics, 3) the structural dip and fault displacement angle estimated (rake), and 4) seismogenic structure width.Based on the described parameters, the faults have been modeled with several short segments instead of long segments (Fig. 5; Chunga, 2010).The historical seismicity of the Guayaquil gulf show seism with Mw comprise between 6 and 7.1 and a recurrence time lower than 15 years.These data can be justified only with fault with longitude lower than 70 kilometers.For this reason, the maximum estimated magnitude will be lower, and the fault slip rate requires of many smaller earthquakes to accommodate a cumulative seismic moment (Well and Coppersmith, 1994).

TABLE 3. CATALOG OF THE CAPABLE FAULTS THAT ARE ABLE TO GENERATE EARTHQUAKES WITH MAGNITUDE
The most appropriate approach to estimate the maximum magnitude of an earthquake is a relationship with the fault rupture length.The maximum magnitudes for each of the crustal faults identified, and the maximum vertical displacement, are estimated with the empirical regression relationships (magnitude earthquake/displacement of geological fault) given by Well and Coppersmith (1994).In equations 1 and 2, L f is the capable fault length.Leonard (2010) proposed changes and corrections on the relation used to estimate the maximum magnitudes, providing the equations, Mw=a*log (L f )+b (3) In equation 3, a=1.52 and b=4.33.The estimated magnitude is given by the equation 4: (M e )=1,52*Log 10 (L f )+4,33 Similarly, Wesnousky (2008), proposed other relations for each type of capable faults: Shear faults: Mw=5.56+0.87*Log(Lf ) (5) Normal faults: Mw=6.12+0.47*Log(Lf ) Reverse faults: Mw=4.11+1.88*Log(Lf ) (7) These regression equations show that different kind of faults, with the same dimension, may generate earthquakes with different magnitudes (Wesnousky, 2008;Stirling et al., 2013).This information was applied on the capable faults outlined in the Gulf of Guayaquil (Fig. 6).The reverse faults are considered capable of generating earthquakes higher than those of equal length but related with shears and normal faults.Another value that we have estimated is the Peak Ground Acceleration (PGA) and used the equation proposed by Fukushima and Tanaka (1990).
The PGA estimated in the analysis ranges from 0.24g to 0.41g (Table 3).These PGA values may be compared with the seismic zonation map of Ecuador (NEC, 2011; Fig. 7).The reliability levels for each of the capable fault have been implemented from seismological (instrumental earthquake records outlined along the structure) and morphological (fault scarps or delineation of reliefs associated with lifting or tectonic subsidence) analysis.Three main categories of capable faults are recognized, being true "I," deductible "II," and uncertain or hypothetical "III" (Table 3).For a reliability level of being "true," it has been required that the geological faults present evidence of dislocation on the ground, seismicity, and lateral displacement during the Quaternary.For a "deduced" level the fault must present displacement FIG. 5. Maximum Credible Earthquakes (MCEs) generated by fault segments capable in Ecuador (Chunga, 2010), confirms the relationship between length and vertical displacement of crustal faults in the Gulf of Guayaquil.
or dislocation evidence that may be distinguished on the morphometry relief.Finally, the hypothetical level corresponds when the structural guidelines may be associated with an active fault, but the displacement cannot be recognized, or where the epicenter of the earthquakes are aligned, with the hypocenter being less than 20 km in depth.

Discussion and conclusions
Forty crustal faults capable of generating earthquakes in the Gulf of Guayaquil have been identified.The relation proposed by Wesnousky (2008) allowed to estimate the associated maximum magnitudes and PGA values in the Gulf of Guayaquil and the coastal segment of the southern part of continental Ecuador.
These faults may generate earthquakes with magnitudes ranging from Mw 6.2 to Mw 7.2.The F-40 reverse fault, located 177 SW km west of the city of Guayaquil, on the accretionary prism area (opposite respect to the subduction zone), has a length of approximately 150 km and can generate a potential earthquake of magnitude Mw 8.2, and a PGA of 0.50 g.The southern coastal area of Ecuador, including Guayaquil, is a densely-populated area with an important strategic industrial sector as well as a petroleum refinery.In this area, earthquakes with a magnitude between Mw 7 and Mw 7.2, and PGA values from 0.34 g to 0.41 g may be generated in a range of 23 and 80 km with respect to the city of Guayaquil.The crustal faults and analyzed, which could generate these strong earthquakes are: F-01, F-08, F-09, F-10, F-11, F-21, F-25, and F-31.They are dextral shear faults, located to the north of the Island of Puná (Fig. 6).
The normal fault F-05 (Fig. 8), located 3 km east of Guayaquil, has a fault striking N300° and dipping N25° S-SW.Based on analysis of morphotectonic features, an earthquake of magnitude Mw 6.7 and a PGA value of 0.35g is estimated.This fault is evident in a quarry having a parallel trend with the strata direction (210/25) of limestones and shales of the San Eduardo and Guayaquil formations (Fig. 8).The Carrizal fault (F-18) or the Estero Salado fault (F-17) has been the responsible of the seismic event passed on June 11 th , 1787 (macroseismic intensity VIII; Chunga, 2010).Both structures have well-defined FIG. 7. Estimation of maximum PGA (Peak Ground Acceleration) obtained from the analysis of capable faults using the method proposed by Fukushima and Tanaka (1990).F-01 to F-40 indicate the geological faults outlined in the southern coast of Ecuador.The graphical representation of the analytically determined values were plotted using the Arcgis and Surfer software, by using the interpolation by "gridding" creating an iso-values map.
morphological and seismological features and an "educated guess" of the PGA of 0.33 g.This value has been estimated from the hypocentral distance of the faults being located between 10 and 16 km depth, and the maximum magnitudes have been estimated to be Mw 6.7, approximately.
For the sub-basin of Esperanza, located in the seafloor of the Gulf of Guayaquil, faults F-24 to F38 (Table 2) are identified as capable faults with estimated PGA values on the rock between 0.33 g and 0.41 g.In this segment of the continental shelf, the maximum estimated magnitudes have values between Mw 6.4 and Mw 7.0.The maximum vertical displacements have been estimated between 0.8 to 1.6 m.
The historical records illustrate that Guayaquil, the most populated city in Ecuador, may be prone to strong earthquake-induced damage, with coseismic effects on the surface, by moderate local earthquakes and very strong of, even, distant earthquakes of magnitude higher than Mw 7.8.

TABLE 2 . CHARACTERIZATION OF SEISMIC SOURCE FOR THE GUAYAQUIL EARTHQUAKE OF 11.06.1787 AND ITS POSSIBLE TECTONIC ASSOCIATION WITH THE ESTERO SALADO FAULT (F17) AND/OR THE CARRIZAL FAULT (F18).
Calculated magnitude based onWell and Coppermisth (1994); Estimated Magnitude from morpho-structural Zonation based on Chunga (2010); Fault length estimated from morphological analysis; Reliability levels for faults are from seismological and morphological analyzes.