Seismically Induced Soil Liquefaction and Geological Conditions in the City of Jama due to the M7.8 Pedernales Earthquake in 2016, NW Ecuador

: Seismically induced soil liquefaction has been documented after the M7.8, 2016 Pedernales earthquake. In the city of Jama, the acceleration recorded by soil ampliﬁcation yielded 1.05 g with an intensity of VIII to IX ESI-07 . The current study combines geological, geophysical, and geotechnical data in order to establish a geological characterization of the subsoil of the city of Jama in the Manabi province of Ecuador. Then, the liquefaction potential index (LPI) has been evaluated considering the PGA-rock values calculated from deterministic methods applied to nearby geological faults, as well as the soil acceleration records for the city of Jama since the Pedernales megathrust earthquake. The importance of conducting geotechnical evaluations of particular colluvial, alluvial, and ﬂoodplain deposits, for which the liquefaction probability proﬁles have been additionally obtained, may serve as a useful tool for ediﬁces foundations or earthquake resistant designs. Finally, the site response analysis is presented using a linear equivalent analysis, where previously seismic records compatible with the target spectrum have been selected. Hereby, the results of ground surface effects have been compared with the spectra of the Ecuadorian Regulation of Construction (NEC) in the context of local seismic ampliﬁcation.


Introduction
There are a variety of geologic effects that may be triggered by strong earthquakes causing severe damage to human settlements as well as strategic infrastructure and even threatening the lives of inhabitants [1][2][3][4][5][6][7][8][9][10][11]. Such secondary coseismic effects may include the shaping and modification of a vast area and corresponding landscapes over a longer period of time due to the generation of landslides, ground shaking, and slope failure as well as tsunamis amongst many others, which may leave behind zones of destruction and or catastrophic consequences [12][13][14][15][16][17][18][19][20][21].
Within the earthquake induced geological and geomorphological hazards, the soil liquefaction is the main cause of destruction that is directly related to earthquakes and has the ability of subsidence, position variation or even collapse for buildings [2,[22][23][24][25][26].

Geodynamic Setting and Description of the Study Area
In the northwestern edge of South America, the coast of Ecuador has a short record of historical earthquakes, which demonstrates the importance of studying seismogenic structures of both, subduction and geological faults, in order to estimate the maximum probable levels of seismicity [3,47]. Similar investigations about high seismic potential and PGA-rock faults have been conducted on the southern, central and northern coastal segments of Ecuador [6,35,47,48].
The predominant seismic source of the Pacific coastline of Ecuador is the process of subduction in the interface zone, which results from the Nazca plate moving below the continental segment formed by the South American and Caribbean plates which occurs at a rate of approximately 47 mm/year [49][50][51][52][53]. The dip variation of the subducting slab, as well as fracture zones of the oceanic crust, define the segments and the corresponding seismic behaviors [47,51,[54][55][56][57].
For the continental coast of Ecuador, four main seismogenic segments have been recognized along the subduction convergent zone ( Figure 1; Table 1; IOC [57]). The first segment is the southern coast, which includes the Santa Elena Province and the Gulf of Guayaquil with a rupture length estimated of about 200 to 220 km. This corresponds to the coastal zone from Salango in southern Ecuador up to the offshore area of Guayaquil Gulf, with a maximum estimated magnitude of M8.1, having a much longer recurrence than the other tectonic segments.
Gulf, with a maximum estimated magnitude of M8.1, having a much longer recurrence than the other tectonic segments. Table 1. List of the parameters of the tectonic segments along the coast of Ecuador, as modified from the IOC [57].    The second segment is the Isla de la Plata around the Manta peninsula with a rupture area of 100-120 km reaching from Bahia de Caráquez to Machalilla, where the estimated magnitude is about M7.8. This is followed by the third segment in the central coast around Manabí Province with a rupture area of 110-120 km, extending from San Vicente to Muisne.

Segment Depht (km) DIP (°) Strike (°) Rake (°) Width (km) Length (km) Slip (km/m) M
The maximum magnitude of M7.9 was documented in Jama in 1942, and the recent M7.8 Pedernales earthquake in 2016 [1,18,33,34,58]. Finally, the fourth segment is located within the northern coast of the Province of Esmeraldas and southern Colombia with a rupture area over 450 km, reaching from the Galera peninsula to the Buenaventura Colombian villages, where the maximum estimated magnitude is about M 8.6 [18,35].
Several authors indicate that the third and fourth segments correspond to the same tectonic delineation (see Table 1, Galera I + II), where the maximum magnitude would reach M 8.9 [59][60][61].
The city of Jama is situated within the third tectonic segment, where the Jama 1942 (M7.9) and Pedernales 2016 (M7.8) earthquakes are the strongest thrust earthquakes recorded during the last 120 years of seismic history [18] (Figure 1). Other earthquakes of lesser magnitude were documented in the same segment, such as the M7.1 earthquake of 3 May 1896, followed by the M7.4 earthquake of 1 June 1907, the earthquake M 7.4 of 16 January 1956 and the earthquake of M7.1 of 4 August 1998 ( Figure 1). Historical seismicity suggests a recurrence period of 20 years for earthquakes M >7.0, and for larger earthquakes of M ≥7.6 a recurrence between 70 to 80 years [35].
Therefore, the NEC-11 seismic zonation map of Ecuador is comparable with the seismotectonic map of Ecuador, as both classify the coast of Ecuador as an area with high levels of seismicity associated with interplate earthquakes, indicating rock acceleration values between 0.45 and 0.55 g [46,47]. Due to its proximity to the subduction zone, the city of Jama is located in the province of Manabí within the central coast of continental Ecuador, whose territory is referred to the highest seismic zone with PGA-rock of ≥ 0.5 g [46].
The recent Pedernales earthquake reported in 2016, registered for the city of Jama a PGA-soil in the order of 1.0 to 1.05 g, while for the site of the epicenter in Pedernales a PGA-soil of 1.4 g was registered (IGEPN). Referring to the strongest aftershocks affecting the city of Jama, the most significant ones were registered on 18 May 2016 with M 6.7 and M 6.9 [18,58,71,72]. These seismic sequences allowed to establish that earthquakes from M6.9 can generate a liquefaction phenomenon in saturated soils where the water table is less than 10 m [18].
The second seismic source is represented by geological faults. In the current study, the seismic parameters of the faults have been consulted from the geometric parameters for each of the selected faults, including: (1) the spatial projection of the fault length on the ground, (2) type and kinematics of the fault, (3) the structural immersion and estimated angle of the fault displacement "in focal mechanism analysis", and (4) the width of the seismogenic structure. With these parameters, and applying the equations most adaptable to the tectonics of the region, the magnitude, macroseismic intensity, and acceleration (PGA-rock, peak ground acceleration-horizontal) have been determined. Segmented faults have a lower magnitude, and therefore the slip rate of the fault requires many smaller earthquakes in order to accommodate a cumulative seismic moment [73][74][75].
The current study comprises an area of 579 km 2 , which considers to be deterministic about the seismic hazards on the entire canton of Jama. Meanwhile, the detailed geological analysis and geodynamic conditions of the subsoil is performed over the urban area of Jama with covering 2.48 km 2 . The Jama canton registers 20,230 inhabitants and represents 1.7% of the total inhabitants of the province of Manabi, with an annual growth of 6.4% [76]. The cantonal capital is the city of Jama with 4719 inhabitants. Most of their houses are built in lowlands of alluvial plains and abandoned terraces of floodplains. The altitudes are between 5 and 10 m a.s.l. The highest reliefs do not exceed 25 m in altitude, being located southeast of the urban area.

Data and Analytical Procedures
For a better understanding of the seismic hazard of Jama, we prepared a geological faults catalog, which included 27 active segments capable of generating earthquakes of 5.3 to 6.5 of moderade magnitudes (M) and PGA-rock in the order of 0.17 g to 0.32 g. The faults are abbreviated from F01 to F27 (Table 2). Larger faults, reaching lengths between 16 and 18 km, such as F01, F12, and F24, have been classified as reverse faults. These faults are prone to generate earthquakes between 6.4 and 6.5 with a PGA of 0.30 g to 0.32 g. Seismicity levels of M6.2 and M6.3 are associated with shear-type faults (F06, F09, F16, F18 and F19), with lengths between 5.7 and 6.8 km, and the accelerations in rock are in the order of 0.29 g.
In this context, the city of Jama is susceptible to two types of seismic sources, being large subduction earthquakes with regional coseismic effects, and local earthquakes of moderade magnitudes associated with the activation by geological faults. Hereby, local earthquakes may reach magnitudes of 6-6.9 and strong earthquakes between 7 and 7.9. For the city of Jama, an acceleration in rock of 0.5 g is determined, which may be amplified in soils reaching values of up to 1.05 g, as estimated for the 2016 Pedernales earthquake (Table 2; [35]).
The evaluation of the liquefaction potential and the earthquake site response in the city of Jama includes three distinctive stages: (i) Field reconnaissance and outlining of Quaternary lithological units of the urban area, and their relationship with the field borehole. In-situ testing and laboratory experiments, allows the construction of a stratigraphic profile of the subsoil up to −45 m deep and determine the basement depth in Jama urban area using active seismic methods (Geotechnical and geophysical data available from Chunga [77] and Daza et al. [78]. (ii) Evaluation of the liquefaction potential index (LPI) through deterministic risk-based methods from SPT-based simplified procedure and shear wave velocity measurements, taking into consideration a seismic hazard scenario calculated for a return period of 475 years which is amax = 0.5 g. This technique supports the detection of the depth of the maximum deformations induced by liquefaction. The approach realized for the evaluation of the liquefaction potential index of soil units is based on results obtained from laboratory tests, has been based on the grain-size characteristics and Atterberg limits of the soil, considering geotechnical parameters proposed by several studies [29,30,41]. (iii) Determination of the geometric mean spectra in rock and soil using an one-dimensional model in the Deepsoil program, and has been compared with the seismic spectra of the NEC-11, for Jama D and E geotechnical soil profile types. Using the shear wave velocity values for the underlaying strata, the rock profile was found to be dense soil or soft rock (Supplementary Materials). Following the NEC-11 [46], for soil quality classification, this type of soil is categorized as C profile, meaning that the soil has bad quality making it prone to have considerable coseismic effects, so that the quality of the filtering is adjusted to the target spectrum, geotechnical parameters typical of the study scenarios were also entered into the PEER Ground Motion Database (Earthquake records from PEER NGA strong motion database, https://ngawest2.berkeley.edu/).
The faults delineation for Jama was generated considering the geomorphological features, direction, and hill reliefs intersections. After this, the dimension of the faults is approximated without oversizing the length of this seismogenic structure. The tectonic activity is associated to similar earthquakes of the upper plate geological faults with hypocentral distances less than 15 km [79].
The equations of Wells and Coppermish [73], used in the analysis of seismic hazard for the upper plate geological fault, are as follows: Estimated magnitude (M) = 5.08 + 1. 16 Fault displacement (in meters) = EXP (−1.38 + 1.02 × LOG (Lf)) (2) where Lf is the length of the capable geological fault. We consider the equations proposed by Wesnosuky [73] and Leonard [80], as these define that a fault of the same length may have different magnitudes depending on the type of faults, being inverse, normal or shear. The spatial delineation of the faults were developed on a GIS platform as illustrated in Figure 2. Wesnousky [73], proposed modifications and corrections to the previous formula in order to estimate maximum magnitudes, such as: Normal faults : M = 6.12 + 0.47 × Log (Lf) Reverse faults : M = 4.11 + 1.88 × Log (Lf) The equations of Wells and Coppermish [73], used in the analysis of seismic hazard for the upper plate geological fault, are as follows: where Lf is the length of the capable geological fault. We consider the equations proposed by Wesnosuky [73] and Leonard [80], as these define that a fault of the same length may have different magnitudes depending on the type of faults, being inverse, normal or shear. The spatial delineation of the faults were developed on a GIS platform as illustrated in Figure 2. Wesnousky [73], proposed modifications and corrections to the previous formula in order to estimate maximum magnitudes, such as:   These regression equations indicate, that not all types of faults of the same dimension are able to generate earthquakes of the same degree of magnitude [47,74,75]. Such a concept is applied in the geological faults delineated for the city of Jama. However, we consider that these faults are inverse as well as potentially capable of generating earthquakes greater than others of the same length, but with different tectonic stresses such as shears and traction.
Furthermore, the values of the peak ground acceleration (PGA-rock) have been estimated using the equation proposed by Fukushima and Tanaka [81], expressed as: where Hf is the depth of fault and M the estimated magnitude. The data obtained from magnitudes and PGA-rock have been interpolated using the ArcMap GIS and Surfer programs, applying the "gridding" method generating iso-value maps, while the PGA data interpolations in rock have been illustrated in Figure 2.

Geomorphological Conditions and Earthquake-Induced Environmental Effects
The Jama urban area has six well-delineated geomorphological features in the terrain such as: intertidal flat, alluvial plain, ancient floodplain, colluvial, abandoned floodplain and soft rocks. The soft rocks are composed of siltstones and claystones, from the Onzole and Borbon formations, in addition to other outcrops from the Punta Blanca, Tosagua, Dos Bocas, Villingota, Angostura formations, and some recently identified quaternary deposits. Their ages range from the Oligocene to the Pliocene, covered by unconsolidated sediments from the Pleistocene to the Holocene [82].
In the road cut of different hill slopes, the stratigraphic sequence of the colluvial soils with angular fragments of highly altered claystones is evidenced, interspersed at some levels by layers of distal rain triggered lahars [83]. The height of the slope reaches between 15 and 20 m. In the terraced areas, on the upper level there are sandy loamy soils, which are not very compact having desiccation cracks, there are some areas with less depressions and being flooded demonstrating anthropogenic subsidence (BH-03 in Figure 3). At the lower level, sandy soils of medium granulometry, somewhat silty, limit the floodplain soils of alluvial plains. In the riverbanks of the Jama River, appear in the alluvial plain areas silt layers with few intercalations of very fine sand, having moderate to high plasticity.
During the Pedernales earthquake of 16 April 2016, cracks in loose soil with a thickness of up to 20 cm are evidenced in the alluvial terrace area, while in other areas of the transitional limits of alluvial plain and terrace were reported lateral displacements with cracks that reached a 30 cm opening in loose soil. One of the most significant damages was the collapse of the CIBV, a nursery for children, located on sandy silty soils in an alluvial plain, where soil liquefaction and cracks of up to 25 cm opening in loose soils considerably deformed the access roads ( Figure 3d).
Shrimp farms and their transition to the alluvial plain, where these types of unstable soils presented soil liquefaction and lateral displacements, occupy the intertidal flats to the north of the urban area ( Figure 3). The areas most prone to coseismic deformations in the terrain are intertidal, alluvial plains, and ancient floodplains. Mixed construction and reinforced concrete buildings had partial and total damage, both located in areas of lower alluvial terrace in the urban area of Jama ( Figure 3f).
In the lateral steps of the bridge over the river Jama, vertical displacements of the metal and concrete structure were evidenced, where the streets were filled with stone materials in order to gain access to the town of Pedernales, due to the emergency of the earthquake and the shipment of humanitarian aid (Figure 3e). In addition, asphalt and reinforced concrete roads had transverse fractures to the axis of the lateral access road to the city of Jama. There, the macroseismic intensity assigned was IX, when applying the environmental intensity scale ESI-2007 [13,18]. sciences 2021, 11, x FOR PEER REVIEW 8 of 22 During the Pedernales earthquake of 16 April 2016, cracks in loose soil with a thickness of up to 20 cm are evidenced in the alluvial terrace area, while in other areas of the transitional limits of alluvial plain and terrace were reported lateral displacements with cracks that reached a 30 cm opening in loose soil. One of the most significant damages was the collapse of the CIBV, a nursery for children, located on sandy silty soils in an alluvial plain, where soil liquefaction and cracks of up to 25 cm opening in loose soils considerably deformed the access roads (Figure 3d).
Shrimp farms and their transition to the alluvial plain, where these types of unstable soils presented soil liquefaction and lateral displacements, occupy the intertidal flats to the north of the urban area ( Figure 3). The areas most prone to coseismic deformations in the terrain are intertidal, alluvial plains, and ancient floodplains. Mixed construction and reinforced concrete buildings had partial and total damage, both located in areas of lower alluvial terrace in the urban area of Jama (Figure 3f).
In the lateral steps of the bridge over the river Jama, vertical displacements of the metal and concrete structure were evidenced, where the streets were filled with stone ma-

Geological and Geophysical Characterization of the Subsoil
A geological profile (A-B) with a length of 2.8 km, trending NW-SE has been performed for the urban area of the city of Jama, where most of the damage to buildings occurred during the M7.8 Pedernales earthquake in 2016 ( Figure 4). The profile is located on alluvial plain deposits in the sedimentary deposition domain of an alluvial valley limited by hillside zones.
The first layer of artificial fill and waste materials is between 1 and 2 metershigh, being attributed to a modern anthropogenic soil. Different Quaternary geological units include the alluvium plain deposits (Qa) with thicknesses ranging between 2 to 4 m, conformed by MH and CL soils of Holocene age. The shear velocity from microtremor ranges from 110 to 150 m/s. The colluvium deposits (Qc) are constituted by ML soils, being of Late Pleistocene age, with thicknesses between 8 and 15 m, and a shear velocity from 200 to 260 m/s. This unit can be divided into two subunits, being more recent colluvial where the N60 reach values of 30, and older and more resistant sediments where the N60 reach values of 50.
The ancient alluvium plain deposits (Qaa) of Holocene to Late Pleistocene ages, are constituted by firm soft soils such as ML, MH and SM, having thicknesses between 6 to 8 m, while their shear velocities ranges from 160 to 210 m/s. In this unit, the rigid silt layer reaches 60 strokes, being the most resistant layer in the entire Middle-Upper Pleistocene sequence.
include the alluvium plain deposits (Qa) with thicknesses ranging between 2 to 4 m, conformed by MH and CL soils of Holocene age. The shear velocity from microtremor ranges from 110 to 150 m/s. The colluvium deposits (Qc) are constituted by ML soils, being of Late Pleistocene age, with thicknesses between 8 and 15 m, and a shear velocity from 200 to 260 m/s. This unit can be divided into two subunits, being more recent colluvial where the N60 reach values of 30, and older and more resistant sediments where the N60 reach values of 50.
The ancient alluvium plain deposits (Qaa) of Holocene to Late Pleistocene ages, are constituted by firm soft soils such as ML, MH and SM, having thicknesses between 6 to 8 m, while their shear velocities ranges from 160 to 210 m/s. In this unit, the rigid silt layer reaches 60 strokes, being the most resistant layer in the entire Middle-Upper Pleistocene sequence.  The rocky basement of Miocene age is constituted by siltstones and claystones of the Onzole and Borbon geological formation [82,84], where the shear velocity must be around 680 m/s. Results obtained from active seismic methods indicate that the Miocene rocks are between 40 and 100 m thick, delineating the lithological contact of the cretaceous basaltic basement between 90 and 100 m deep from the surface [77,78,85], and where the Vs reach 2800 m/s. Table 3, lists the characteristics of the geological materials of the six Quaternary units of Jama city and the Onzole, Bobon, and Piñon geological formations.

Geotechnical Soil Characterization and Liquefaction Susceptibility
The consequences of liquefaction may be very different according to site characteristics and being classified. Furthermore, liquefaction may occur due to cyclical mobility and flow [86]. Cyclic mobility is the phenomenon in which there is a progressive increase in pore pressure until it equals the effective restricting pressure, which means, that there is an accumulation of deformation but it does not fail and what exists is a loss of rigidity. Such a phenomenon is able to cause considerable damage to the structure [87].
The flow failure is generated by a sudden loss of resistance to an increase in pore pressure caused by the undrained stress. Furthermore, there are excessive deformations when considering flow failure, which refers to a loss of resistance in which the mass of soils is required to be subjected to a static shear stress greater than the undrained shear strength. Compared with that of cyclical mobility, the effects of flow failure can be classified as catastrophic [86].
In order for flow liquefaction to occur, the geotechnical criteria of having the soil saturated 100% need to be met. This tends to occur, if there is an undrained load which originates from different sources such as a seismic load from an earthquake of moderate to strong magnitude, or through waves that may be caused by construction systems such as an explosion or vibrations that are used for geotechnical improvement of soils or in the driving of piles that are also sources of waves. Other conditions are that the tendency of the soil during cutting is contractive and that the soil is susceptible to liquefy. These soils may be clean and uniform sands, usually non-plastic silts and, in rare cases, gravels and clays that are also able to liquefy but with little to none catastrophic effects [87].
In the current study, the liquefaction susceptibility was analyzed for soils that reported previously coseismic damage during the 2016 earthquake. The epicenter of the M7.8 earthquake was located some 50 km NE of the city of Jama. Several researchers of earthquake engineering have studied the liquefaction susceptibility of geologic units of recent environmental depositions, proposing criteria for their classification [23,31,37].
According to these guidelines, the Jama urban area is meeting the following criteria: (a) the urban area appears to have been constructed on the geological units such as alluvium plain and ancient floodplains (abandoned terrace) deposits, of Holocene to late Pleistocene age; (b) the groundwater is less than 5 m deep and the peak ground acceleration (PGA) has a 10% probability of being exceeded in 50 years in the range between 0.45 and 0.55 g [46,72,88,89]; (c) evidence of historical liquefaction surely were reported after the 13 May 1942 (M 7.9) and the 16 April 2016 (M7.8) megathrust earthquakes [1,34,35,58]. In liquefaction studies it is important to determine the depth of the potentially liquefiable strata, and of these, which layer may have greater probability of liquefaction [36]. In the current study, the maximum probability will be determined for liquefaction in deep strata.
As a first approach, particular attention is referred to the stratigraphy of the urban area, using a variety of study techniques [41,[90][91][92]. Then, the empirical approach of these stratigraphic units are classified as "probably liquefiable", when the liquid limit is less than 37% and its plasticity index is less than 12% [90], always that the water content is high (wc > 0.8 LL). Consequently, a re-evaluation of the geotechnical parameters from three boreholes was conducted (BH1, BH2 and BH3). The Municipality of Jama, and the University of Los Andes from Colombia and some private soil laboratories supplied the geotechnical data, which have been required to estimate the liquefaction potential according to the proposed methodologies [28,29,41,42] have provided these data [77,78].
Hereby, the majority of the drilling (boreholes) are delineated in the liquefiable area [90]. However, this criterion is not fully compatible with the low plasticity behavior of fine-sized soils of the Jama city, where historical evidence of liquefaction has been documented in the alluvium plain and ancient floodplain deposits. Figure 5 illustrates the results of the liquidity index and the water content of the soil of the city of Jama, used to determine the scale of the natural moisture content of the soil samples with reference to the liquid and plastic limits. Hereby, it provides an indicator of geological history presenting these soils with ultra-sensitive sand deposits. Geosciences 2021, 11, x FOR PEER REVIEW 13 of 22

Liquefaction Potential Mapping
Although the available stratigraphic information reaches 35 m of depth, the recorded coseismic deformations in the soils of Jama are less than 20 m deep. Therefore, a depth of 20 m has been considered for the liquefaction potential assessment. The liquefaction potential index (LPI) has been calculated in the urban area of Jama using the equation defined by Iwasaki [96] as: where z is the depth below the ground surface in meters and is calculated as w(z) = 10-0.5z; F(z) is a function of the safety factor and the liquefaction factor, SF, where F(z) = 1-SF when SF < 1 and if SF > 1 than F(z) = 0. Iwasaki [96] calibrated the severity of liquefaction induced damages as, LPI = 0 as not prone to liquefaction and 0 < LPI < 2, 2 < LPI < 5, 5 < LPI < 15, and LPI ˃ 15 as low, moderate, high, and very high susceptibility to liquefaction, respectively. The probability of liquefaction was calculated for all the boreholes using the logistic regression proposed by Chen and Juang [41] and by Juang et al. [92]: The results of the liquefaction probability for the different boreholes are listed in Table 4 and in Figure 5. The results indicated that in the urban area of Jama city there is a high probability that "almost certain" soil liquefaction will occur in the Holocene sedimentary layers according to the method proposed by Iwasaki [96]. Furthermore, there are very likely sectors where no soil liquefaction will occur.
The liquefaction susceptibility of the urban area of the city of Jama was calculated considering the data of the previous section. This result delineates the areas according to their probability of liquefaction. The liquefaction probability mapping of the area of study is obtained by interpolation of the liquefaction potential index (LPI) obtained for each borehole, and processed with the Arc-Map GIS and Surfer. The coseismic liquefaction evidences of the Pedernales earthquake (M7.8) are confronted with the results of the LPI data [96].
Regarding the geological section, the recent alluvium plain and ancient floodplain deposits have (referring to the location of BH-01) a very high probability of liquefaction up to 15 m depth (0.85 ≤ PL ≤ 1.00; 0.653 ≥ SF > 0.000). In the interface zone with the ancient alluvial plain deposits, the probability of liquefaction is moderate with up to 23 m depth The analysis for the urban area of the city of Jama demonstrates that there is a deep stratum of late Pleistocene age which may also be potentially liquefiable when the soils meet a saturation close to 100%. The statistical summary of all these geotechnical parameters for the city of Jama is given in Table 4. From the results listed in this table and the particle size distribution results (Figure 5), the following statements may be made: (a) the saturation in alluvium plain deposits is about 95% whilst for the ancient alluvium deposits it is about 85%; (b) the fraction of fines particles (Fc) ranges between 1% and 61% for alluvium plain deposits and from 75% to 97% for the ancient alluvium deposits; (c) the diameters with 50%, D50, is between 0.02 and 0.20 mm for alluvium plain and ancient alluvium; (d) the coefficient of uniformity (CU) lies between 1.7 and 2.1 and the degree of curvature between 0.8 and 2.5 for fine to coarse silt and fine sand according to the soil classification. The analysis of the previous data indicates that the alluvium plain and ancient floodplain deposits are the liquefiable soils of the area (Figure 4). The distribution of the SPT values in function of depth in the urban area of the city of Jama is illustrated in Figure 5.

Evaluation of the Liquefaction Potential
In order to evaluate the earthquake-induced liquefaction potential of the Late Pleistocene -Holocene geological units in the city of Jama, some risk-based methods from an SPT-based simplified procedure and shear wave velocity seismic measurements were considered. These methods were originally developed by Seed and Idriss [28] and updated and modified later [30,93,94]. These deterministic methods comprise the calculation of the safety factor per layer (SF), defined as the ratio between cyclic resistance ratio (CRR) and the cyclic stress ratio (CSR).
The CRR, according to Youd and Idriss [93] has been calculated by using the following equation: where (N1)60 is influenced by the measured standard penetration resistance N, relative to the overburden pressure factor (Cn), the correction for hammer energy ratio (ER) Ce, the correction for borehole diameter Cb, the correction factor for rod length Cr, and the correction for samplers with or without liners Cs. The Cn was calculated according to the equation proposed by Liao and Whitman [95], Cn = (Pa/σ'v)0.5 in function with (Pa) (atmospheric pressure) and the σ'v (effective vertical stress). This has been followed by the application of a "fine content" correction to the calculated N1(60) value in order to obtain a clean-sand-equivalent value N1(60) Cs given by the equations proposed by Youd et al. [94] as illustrated in Figure 5. The CSR defines the seismic demand and is expressed as: where σv is the total vertical stress at depth z, σ'v is the effective vertical stress at the same depth, amax is the peak horizontal ground acceleration, (g) is the acceleration due to gravity, and rd is the stress reduction factor. In this study, the term "rd" was estimated using the Liao and Whitman [95] (13) where FCI is an index of fines content (FC) defined as follows: FCI = 1 for FC <5%, FCI = 2 for 5% < FC < 12%, FCI = 3 for 12% < FC < 35%, and FCI = 4 for FC >35%. The use of the ordinal scale to characterize the effect of fines content is consistent with current geotechnical knowledge (Chen and Juang [41]). The values of CSR and CRR will result in the estimation of the safety factor (SF). The results of the safety factor (SF), for the different geological units in the city of Jama are illustrated in Table 4. All the soil layers with a safety factor higher than 1.25 are considered as non-liquefiable soil layer according to Seed and Idriss [91]. Below the suggested value of 1.25, the soil layer has been considered to be prone to experiment liquefaction, when complying the liquefaction susceptibility criteria suggested by Chen and Juang [41] and Juang et al. [92].

Liquefaction Potential Mapping
Although the available stratigraphic information reaches 35 m of depth, the recorded coseismic deformations in the soils of Jama are less than 20 m deep. Therefore, a depth of 20 m has been considered for the liquefaction potential assessment. The liquefaction potential index (LPI) has been calculated in the urban area of Jama using the equation defined by Iwasaki [96] as: (14) where z is the depth below the ground surface in meters and is calculated as w(z) = 10-0.5z; F(z) is a function of the safety factor and the liquefaction factor, SF, where F(z) = 1-SF when SF < 1 and if SF > 1 than F(z) = 0. Iwasaki [96] calibrated the severity of liquefaction induced damages as, LPI = 0 as not prone to liquefaction and 0 < LPI < 2, 2 < LPI < 5, 5 < LPI < 15, and LPI > 15 as low, moderate, high, and very high susceptibility to liquefaction, respectively. The probability of liquefaction was calculated for all the boreholes using the logistic regression proposed by Chen and Juang [41] and by Juang et al. [92]: The results of the liquefaction probability for the different boreholes are listed in Table 4 and in Figure 5. The results indicated that in the urban area of Jama city there is a high probability that "almost certain" soil liquefaction will occur in the Holocene sedimentary layers according to the method proposed by Iwasaki [96]. Furthermore, there are very likely sectors where no soil liquefaction will occur.
The liquefaction susceptibility of the urban area of the city of Jama was calculated considering the data of the previous section. This result delineates the areas according to their probability of liquefaction. The liquefaction probability mapping of the area of study is obtained by interpolation of the liquefaction potential index (LPI) obtained for each borehole, and processed with the Arc-Map GIS and Surfer. The coseismic liquefaction evidences of the Pedernales earthquake (M7.8) are confronted with the results of the LPI data [96].
Regarding the geological section, the recent alluvium plain and ancient floodplain deposits have (referring to the location of BH-01) a very high probability of liquefaction up to 15 m depth (0.85 ≤ PL ≤ 1.00; 0.653 ≥ SF > 0.000). In the interface zone with the ancient alluvial plain deposits, the probability of liquefaction is moderate with up to 23 m depth (0.35 ≤ PL ≤ 0.65; 1.102 ≥ SF > 0.837), while beyond that depth, the ancient colluvium alluvium deposits, the probability of liquefaction is zero (0.00 ≤ PL ≤ 0.15; ∞ ≥ SF > 1.411).
In the central park of Jama (referred to the location of drilling BH-02), the lithological units of ancient floodplain and abandoned floodplain (terrace) deposits have greater thicknesses, as paleo-channels are evidenced in seismic refraction profiles. Recent alluvium plain deposits have a very high probability of liquefaction up to 6 m depth, leveed channels of silty sand one meter high have a big probability of liquefaction (0.65 ≤ PL ≤ 0.85; 0.837 ≥ SF > 0.653). The sedimentary depositions of greater thickness in paleo-channels are filled by ancient floodplain deposits. In this unit the probability of liquefaction is moderated, up to 25 m of depth, while beyond that depth the sediments are classified as very rigid to hard clay sands, having zero probability of liquefaction.
The colluvial soils have thicknesses between 8 and 15 m, while underlying this lithological unit occur ancient alluvial plain deposits (referring to its location to well BH-03). In these two units the probability of liquefaction is high to very high, within a depth of 23 m depth. Beyond that depth of up to 30 m, ancient colluvium alluvium deposits, from firm to rigid silt, have a low probability of liquefaction (0.15 ≤ PL ≤ 0.35; 1.411 ≥ SF > 1.102). The probability of liquefaction of soil strata and safety factors have been interpreted based on the classification proposed by Chen and Juang [41].
These results revealed that the urban area of Jama presents soils with high liquefaction probability. The liquefaction probability (PL) are delineated for those recent sedimentary depositions like alluvium plain, ancient floodplain and abandoned floodplain deposits of late Pleistocene-Holocene age, where the geological and hydrological conditions are favorable for ground-coseismic liquefaction (Table 4; Figure 5).

Site Response Analysis
The city of Jama is located in zone VI, indicating a high seismic level with a factor of 0.5 g, with a 10% exceedance in 50 years and a return period of 475 years (NEC-11, [46]). The V S30 seismic data were obtained from downhole geophysical tests, indicating a better contrast in the stratigraphy of the alluvial and alluvial-colluvial soils of Jama. For colluvial soils, up to a depth of 15 m, shear wave velocities higher than 300 m/s were measured, while moving further below to older colluvium alluvium and alluvial valley deposits of medium consistency, where the soil velocity reached 440 m/s. Shear wave velocity is the most important variable when establishing how the ground will behave in the event of seismic excitation. Therefore, it is vital to have extensive information and knowledge of the shear wave velocity profiles for the different perforations.
With these geotechnical characteristics and corroborating the technical specifications of the NEC-11 [46], the soils of Jama are classified as geotechnical Type D. For recent alluvial and floodplain soils, there are granular levels between 10 and 15 m of depth, then up to 28 m their velocities are reduced to 120 m/s. In this interval, sediments of medium consistency begin with shear wave velocities up to 320 m/s. The geotechnical classification for this soil is type E (NEC-11 [46]).
Furthermore, the seismic zone factor Z has been obtained for more resistant soils (greater than 30 m in depth), designating the C profile for the geological units of ancient colluvium alluvium and alluvial valley fill deposits. These units meet the velocity criteria according to NEC-11, with values that fluctuate between 760 ≥ vs. ≥ 360 m/s, and the considerations of the soil amplification coefficients Fa, Fd , and SF reflect a layer with characteristics of dense to very dense non cohesive soil. With these, an elastic response spectrum of accelerations Sa is determined, expressed as a fraction of gravity for the different types of sedimentary units in the city of Jama (Figures S6-S10; Supplementary Materials).
The active seismic data allowed to estimate the soil/rock contact for Quaternary deposits, in a depth between 50 to 55 m, and for colluvial soils, in a depth between 8 and 15 m. The predominant rocks in the area are claystones from the Onzole formation of Tertiary age, below the volcanic basement of Cretaceous basalts of the Piñon formation [82,97].

Selection of Earthquakes with Compressive Tectonic Settings
For this earthquake analysis, ten earthquakes with similar tectonic characteristics have been selected for the Jama canton, central region of Manabí. Therefore, the PEER ground motion database was used, in which the period coordinates and spectral accelerations were entered as target spectrum, of the NEC-11 [46] spectrum ( Figure 6), for soil types C (dense soil or soft rock). The Chi Chi earthquake in Taiwan of M 7.6, associated with a seismogenic structure with inverse tectonic stresses, and a rupture surface of 80 km [98], has been selected for our site response analysis, considering ten seismic records with different horizontal components. The seismic traces for the vertical components have not been considered in this study, due to its low probability of significant earthquake environmental effects. Subsequently, a calibration of the scale factors was performed for the selected earthquakes, maintaining at all times a scale factor of 1 to 3, for a return period (Tr) of 475 years, indicating the seismological characteristics of the scaled movements (Figure 7). [98], has been selected for our site response analysis, considering ten seismic records with different horizontal components. The seismic traces for the vertical components have not been considered in this study, due to its low probability of significant earthquake environmental effects. Subsequently, a calibration of the scale factors was performed for the selected earthquakes, maintaining at all times a scale factor of 1 to 3, for a return period (Tr) of 475 years, indicating the seismological characteristics of the scaled movements (Figure 7). Figure 6. Elastic seismic spectra of NEC-11 [46] accelerations for soil and rock profiles.
The seismogenic structures allow an estimation of the maximum magnitudes and accelerations as expected in rocks [99][100][101]. Hence, the prediction of the movement in the ground under the action of a seismic event is typically obtained by a one-dimensional site response analysis. This analysis is usually performed using equivalent linear models since it requires direct properties of the soils and a simple computational calculation. One of the most well-known and calibrated linear equivalent models is implemented in the software of Deepsoil [44], which calculates the response of a system of homogeneous, viscoelastic layers of an infinite horizontal limit, subjected to a motion of the shear wave traveling vertically. Figure 6. Elastic seismic spectra of NEC-11 [46] accelerations for soil and rock profiles. Elastic seismic spectra Acceleration response spectrum (5% structural damping) for ten scaled seismic movements (gray lines) for inverse component earthquakes, with seismological characteristics similar to the expected uniform hazard spectrum, 10% probability of exceedance in 50 years, the median and median ± 1 standard deviation of the scaled movements for NEC-type rock profile.
The non-linearity of the soil is represented by the use of the equivalent linear soil properties through an iterative process. In the Linear Equivalent model process a rock outcrop motion, typically a recorded motion at the surface, is taken and converted into an "internal rock motion". This "internal" rock movement is then converted from the time domain that is, a temporal evolution of accelerations, to the frequency domain, while adding harmonic waves of different frequencies, where each has its own amplitude. However, under a compatible base rock condition, where the impedance ratio is not infinite, an "internal" rock movement can be affected by the impedance ratio of the rock-soil interface, the mass of the soil deposit and its response characteristics, such as the characteristic period [102].
The appropriate shear modulus and damping value, which are dependent on deformations, are set for each layer during each iteration. When convergence is reached, the problem turns to a viscoelastic material with constant properties. With these constant shear modulus, damping, and unit weight values for each layer, the wave equation can be solved using a complex stiffness response method. A transfer function based on this model is calculated to relate the movement at the base level, that is, internal rock movement, to the movement at the ground surface.
Regarding the stress-strain cyclical hysterical behavior of the models, these can be validated with laboratory curves available in the literature, in this specific case [103], were used for the fine fraction and [28] for the coarse fraction, which provides the degradation curves of the shear modulus G/Gmax (-) and the % damping ratio are obtained with the program of Deepsoil V7. Elastic seismic spectra Acceleration response spectrum (5% structural damping) for ten scaled seismic movements (gray lines) for inverse component earthquakes, with seismological characteristics similar to the expected uniform hazard spectrum, 10% probability of exceedance in 50 years, the median and median ± 1 standard deviation of the scaled movements for NEC-type rock profile.
The seismogenic structures allow an estimation of the maximum magnitudes and accelerations as expected in rocks [99][100][101]. Hence, the prediction of the movement in the ground under the action of a seismic event is typically obtained by a one-dimensional site response analysis. This analysis is usually performed using equivalent linear models since it requires direct properties of the soils and a simple computational calculation. One of the most well-known and calibrated linear equivalent models is implemented in the software of Deepsoil [44], which calculates the response of a system of homogeneous, viscoelastic layers of an infinite horizontal limit, subjected to a motion of the shear wave traveling vertically.
The non-linearity of the soil is represented by the use of the equivalent linear soil properties through an iterative process. In the Linear Equivalent model process a rock outcrop motion, typically a recorded motion at the surface, is taken and converted into an "internal rock motion". This "internal" rock movement is then converted from the time domain that is, a temporal evolution of accelerations, to the frequency domain, while adding harmonic waves of different frequencies, where each has its own amplitude. However, under a compatible base rock condition, where the impedance ratio is not infinite, an "internal" rock movement can be affected by the impedance ratio of the rocksoil interface, the mass of the soil deposit and its response characteristics, such as the characteristic period [102].
The appropriate shear modulus and damping value, which are dependent on deformations, are set for each layer during each iteration. When convergence is reached, the problem turns to a viscoelastic material with constant properties. With these constant shear modulus, damping, and unit weight values for each layer, the wave equation can be solved using a complex stiffness response method. A transfer function based on this model is calculated to relate the movement at the base level, that is, internal rock movement, to the movement at the ground surface.
Regarding the stress-strain cyclical hysterical behavior of the models, these can be validated with laboratory curves available in the literature, in this specific case [103], were used for the fine fraction and [28] for the coarse fraction, which provides the degradation curves of the shear modulus G/Gmax (-) and the % damping ratio are obtained with the program of Deepsoil V7.
The geometric mean of the spectra of elastic acceleration responses obtained for the horizontal components of the earthquakes compatible for the Profile D and E soils are presented. With return periods of 475 years for the Equivalent linear analysis for these spectra has been considered a structural damping equal to 5%. The inclusion of the geometric mean of the selected records allowed us to appreciate the amplification that occurs due to the diverse depositional environmental strata found at the site ( Figure 8). For the units of colluvial sediments (referred to BH-03), spectral acceleration of 1.68 g was obtained for a period between 0.5 s and 0.7 s, while for the recent alluvial and floodplain sediments (referred to BH-01), a peak acceleration of 1.04 g was obtained for periods between 0.87 s and 1.11 s. Finally, for the ancient alluvial and abandoned flooplain sediments (referred to BH-02, central park of Jama) yielded an acceleration of 1.03 g as obtained for periods between 1.3 and 1.5 s. In this context, the thickness of Quaternary deposits has a direct impact on the fundamental vibration period of the deposit, and therefore on the shape and the expected spectral amplification, being a variable that generates high sensitivity.
Therefore, we may conclude that the spectra obtained on the surface are greater than those of the acceleration plateau of the Ecuadorian Construction Regulations (NEC-11, [46]). Additional graphs were obtained, where the response spectra at the top of the rock is compared with the spectra of the soil profiles. Graphs were generated demonstrating the response comparison both on the surface and on the rock for the three boreholes.
The results highlighted the importance of geotechnical assessments of particular geological environments which have sediments prone to liquefaction. This tool provides a better understanding of how the constructions and buildings could be affected by an earthquake of high magnitude depending on the type of soil. This was a non-existing tool that could help engineers to improve urban area planification or buildings stabilization to lessen the severity on earthquake's effects in the city of Jama and other areas with similar geological conditions. It is important to mention that the areas displayed in the liquefaction hazard map were contrasted with the on-site data collected after the Pedernales earth- Finally, for the ancient alluvial and abandoned flooplain sediments (referred to BH-02, central park of Jama) yielded an acceleration of 1.03 g as obtained for periods between 1.3 and 1.5 s. In this context, the thickness of Quaternary deposits has a direct impact on the fundamental vibration period of the deposit, and therefore on the shape and the expected spectral amplification, being a variable that generates high sensitivity.
Therefore, we may conclude that the spectra obtained on the surface are greater than those of the acceleration plateau of the Ecuadorian Construction Regulations (NEC-11, [46]). Additional graphs were obtained, where the response spectra at the top of the rock is compared with the spectra of the soil profiles. Graphs were generated demonstrating the response comparison both on the surface and on the rock for the three boreholes.
The results highlighted the importance of geotechnical assessments of particular geological environments which have sediments prone to liquefaction. This tool provides a better understanding of how the constructions and buildings could be affected by an earthquake of high magnitude depending on the type of soil. This was a non-existing tool that could help engineers to improve urban area planification or buildings stabilization to lessen the severity on earthquake's effects in the city of Jama and other areas with similar geological conditions. It is important to mention that the areas displayed in the liquefaction hazard map were contrasted with the on-site data collected after the Pedernales earthquake.
Furthermore, the elastic design spectra were analyzed for a return period of 475 years, considering a level of structural damping equal to 5% and corresponding to an on-site response, that is, the interaction between soil and soil has not been considered foundation of the structure. For the structural calculation, it would be recommended to use an envelope of the spectra obtained in each of the analysis. Site response analysis was performed using the linear equivalent method and in terms of effective vertical stress.

Conclusions
The city of Jama is located in a high sesmicity area, having rock accelerations, PGArock of 0.5 g and its ground amplification can be in the range of 1.05 and 1.68 g.
The obtained results based on the liquefaction potential assessment method demonstrated that the main urban area of Jama exhibits a high liquefaction potential and thereby a very high probability of liquefaction up to 15 m depth, when considering a maximum acceleration of 0.5 g.
It was also evidenced that the shear wave velocity is the variable with the greatest weight when establishing how the soil is going to behave due to seismic excitation. Thus, for colluvial soils there are higher spectral amplifications of 1.68 g with periods between 0.5 and 0.7 s. Therefore, the obtained results indicate that soil deposits tend to have important degradations at high acceleration levels.
As the intensity level of the input signal increases, the maximum spectral amplification peaks will be at higher structural periods and will have amplitudes lower due to the increased hysteretic damping of the soil. Thus, by presenting low wave velocity and clay soils with high plasticity soil profiles, a very high loss of shear stiffness will be common due to the deformations imposed by the input signals.
Supplementary Materials: The following are available online at https://www.mdpi.com/2076 -3263/11/1/20/s1, Figure S1: SPT test results vs. depth. Figure S2: Results of the liquefaction potential index (LPI). Figure S3: Fine content. Figure S4: Plastic Index. Figure S5: LSN chart. Figure  S6: Records time-history and dynamic parameters of one of the earthquakes considered for the analysis corresponding to 475 years of return period (PEER, Center for the Development of Seismic Engineering of the Pacific). Figure S7: Dynamic parameters of one of the earthquakes considered for the analysis corresponding to 475 years of return period (PEER, Center for the Development of Seismic Engineering of the Pacific). Figure S8: Comparison of acceleration spectra for the 5% structural damping for the horizontal components of the 10 earthquakes for coluvial deposits (ie. BH-03). Figure  S9: Comparison of acceleration spectra for the 5% structural damping for the horizontal components of the 10 earthquakes for alluvial soils (ie, BH-01). Figure S10: Comparison of acceleration spectra for the 5% structural damping for the horizontal components of the 10 earthquakes for abandonated flooplain deposits (ie. BH-02). Figure S11: Rock response spectrum for JA-BH-03 and comparison with NEC C spectrum. Figure S12: Rock response spectrum for BH-01 and comparison with NEC C spectrum. Figure S13: Rock response spectrum for BH-02 and comparison with NEC C spectrum. Figure S14: Rock response spectrum for BH-03 and comparison with NEC D spectrum. Figure S15: Rock response spectrum for BH-01 and comparison with NEC E spectrum. Figure S16: Rock response spectrum for BH-02 and comparison with NEC E spectrum.