Starting on the day after the mainshock, we mapped in the field and compiled all the available observations on earthquake environmental effects (EEEs) caused by the $Mw$ 7.8 subduction earthquake that hit the coastal region of Ecuador on 16 April 2016 (the Pedernales earthquake). These effects include: (a) permanent ground deformation, (b) open cracks, (c) liquefaction, (d) landslides, and (e) tsunami waves. We use these observations to evaluate the macroseismic field of the Environmental Seismic Intensity (ESI‐07) scale and compare our results with published macroseismic data collected using traditional damage‐based intensity scales. We found systematic difference in the assessed earthquake intensity and in the intensity attenuation with distance. A comparison of our dataset with the macroseismic data of six instrumental large‐subduction earthquakes that occurred between Ecuador and Chile suggests that, for the Pedernales earthquake, the ESI‐07 works better in the near field, whereas damage‐based scales are more reliable in the far field. Therefore, we generate a single integrated intensity dataset for the Pedernales earthquake, including data measured with the European Macroseismic Scale 1998 (EMS‐98), the modified Mercalli (MM) scale, and the ESI‐07 scale. This new integrated macroseismic field provides reliable and comparable information on the seismogenic source, suggesting that a similar approach could be successfully applied to refine the source model of past subduction earthquakes in Ecuador and elsewhere.

It has long been recognized that earthquakes change the stress in the upper crust around the fault rupture and can influence the behaviour of neighbouring faults and volcanoes. Rapid estimates of these stress changes can provide the authorities managing the post-disaster situation with valuable data to identify and monitor potential threads and to update the estimates of seismic and volcanic hazard in a region. Here we propose a methodology to evaluate the potential influence of an earthquake on nearby faults and volcanoes and create easy-to-understand maps for decision-making support after large earthquakes. We apply this methodology to the Mw 7.8, 2016 Ecuador earthquake. Using Sentinel-1 Interferometric Synthetic Aperture Radar (InSAR) and continuous GPS data, we measure the coseismic ground deformation and estimate the distribution of slip over the fault rupture. We also build an alternative source model using the Global Centroid Moment Tensor (CMT) solution. Then we use these models to evaluate changes of static stress on the surrounding faults and volcanoes and produce maps of potentially activated faults and volcanoes. We found, in general, good agreement between our maps and the seismic and volcanic events that occurred after the Pedernales earthquake. We discuss the potential and limitations of the methodology.

We provide a dataset of the landslides induced by the 2016 Pedernales megathrust earthquake, Ecuador (Mw 7.8, focal depth of 20 km) and compare their spatial distribution with mapped bedrock lithology, horizontal peak ground acceleration (PGA-h) and the macroseismic intensity based on earthquake-induced environmental effects (ESI-07). We studied 192 coseismic landslides (classified as coherent, disrupted and lateral spreads) located in the epicentral area, defined by the VII to IX_ESI-07 isoseismals. Based on our findings, lahar deposits, tuffs and volcanoclastic units are the most susceptible to landslides occurrence. Alluvial plains with fluvial loose fine sand are the most susceptible setting for lateral spreading, with a maximum intensity of IX_ESI-07. The coherent landslides are frequently found in altered shale and siltstone geological units with moderate slopes (8°–16°), with typical intensity ranging between VII and VIII_ESI-07. Our analysis draws a typical framework for slope movements triggered by subduction earthquakes in Ecuador. The most dangerous setting is the coastal region, a relatively highly urbanized area located near the epicenter and where liquefaction can trigger massive lateral spreading events. Coherent and disrupted landslides, dominating the more internal hilly region, can be triggered also in moderate slope settings (i.e., less than 10°). Indeed, the regression analysis between seismic intensity, PGA-h and landslide occurrence shows that most of the events occurred at PGA-h values between 0.4 g and 1.2 g, at a distance of 30 to 50 km from the rupture plane. Our database suggests that lithology and hillslope geometry are the main geological/geomorphological factors controlling coseismic landslides occurrence; while the distance from the rupture plane plays a significant role on determining the landslide size. Finally, we underline that coseismically-triggered landslides are among the most common environmental effects occurring during large subduction events that can be effectively used to properly evaluate the earthquake macroseismic field. The landslide inventory we compiled is suitable for assessing the vulnerability of physical environment from subduction earthquakes in Ecuador, and offers a primary data source for future worldwide analysis.

Seismically induced soil liquefaction has been documented after the M7.8, 2016 Pedernales earthquake. In the city of Jama, the acceleration recorded by soil amplification 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 floodplain deposits, for which the liquefaction probability profiles have been additionally obtained, may serve as a useful tool for edifices 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 amplification.

The city of Portoviejo in coastal Ecuador was severely affected during the 16 April 2016, Pedernales earthquake (Mw 7.8). Various coseismic liquefaction phenomena occurred, inducing lateral spreading, sand boils, ground subsidence, and sinkholes in soils with poor geotechnical quality in the alluvial and alluvial–colluvial sedimentary environment. Therefore, the main aim of this study was to collect data from standard penetration tests (SPT) and shear velocity and exploratory trenches and to calculate the liquefaction potential index (LPI) by considering a corresponding seismic hazard scenario with an a_max = 0.5 g. From these data, a liquefaction hazard map was constructed for the city of Portoviejo, wherein an Fs of 1.169 was obtained. It was determined that strata at a depth of between 8 and 12 m are potentially liquefiable. Our quantitative results demonstrate that the city of Portoviejo’s urban area has a high probability of liquefaction, whereas the area to the southeast of the city is less sensitive to liquefaction phenomena, due to the presence of older sediments. Our results are in accordance with the environmental effects reported in the aftermath of the 2016 earthquake.