M 8.2 earthquake in northern Chile

Alright. Lets hope not many people are harmed as a result of this earthquake and likely large tsunami. We all know what to do. Drop, duck, cover, and hold on. Run to high ground. Stay there.

Check back here for updates throughout the night. I am posting information about the tsunami here.

UPDATE: 2023.04.01 See below the original report.

Here is an earthquake that ruptures today in the region of a recent swarm of earthquakes in the northern Chile subduction zone. This is the USGS web page for this earthquake.

Here is the USGS web page for a large aftershock (M 6.2).

Here is a great animation from IRIS that shows the tectonics and seismology associated with this earthquake.

This is a great animation of the seismic waves travelling through the transportable array.


This is the early moment tensor (Mww), which shows a pure thrust/reverse slip.

Here is a map showing the region with earthquakes from the last 30 days. The largest dot is today’s great earthquake.


This map is zoomed into the region of the swarm and mainshock. The earthquakes that ruptured prior to the megathrust earthquake may be now considered generally foreshocks. The colored lines represent the depth of the subduction zone fault shown in km. The two main red depth contours are 20 and 40 km respectively.

Here is a map put together by the Belfast Telegraph.

This is a map that shows instrumental and historic earthquake rupture regions. This figure from Chlieh et al. (2011) shows the slip models (earthquake slip in meters) and earthquake slip regions for pre-seimologic (prior to seismometers) earthquakes in grey. The swarm of earthquakes from this March are in the northern region of the 1877 M 8.8 subduction zone earthquake.


Here is an estimate of ground shaking intensity, with contours offshore and the fault slip region plotted at 21:00 PST:

This is the seismograph from the HSU seismometer, thanks to HSU Department of Geology facebook page.

This map shows a large M 6.2 aftershock right near the coastline. Lets hope people are safe from all these earthquakes. There is also a M 5.5 further south.

Here is the USGS PAGER page. These plots are generated by overlaying shaking intensity with estimates of population and infrastructure. These overlays are “multiplied” to predict what number of fatalities and amount of economic losses might be (in a probabilistic way). These are just models, but do help international and national efforts that may need to be arisen. Most earthquakes that I have shown these PAGER plots for have had low exposure, but this M 8.2 earthquake has possibly affected up to or over 100 fatalities and up to a hundred to hundreds of millions of USD.

This is an interesting plot, showing the coupling ratio. 0 means that the fault is not locked, or is slipping aseismically. 1 means that the fault is completely locked. The darker the map, the higher their estimate of % locking. This is from Metois et al., 2013.

This is the fault slip model just generated by the USGS. The red star represents the hypocenter and the red squares represent regions of higher slip. Notice how the region of highest slip is not collocated with the hypocenter.

This is the source time function, which tells us how much energy is released over time. This is generated by integrating the energy amplitudes from the seismic waves generated by the earthquake. If we were to go coring in a month, we would possibly find a turbidite with two main pulses.
3

UPDATE: 2023.04.01

In 2023, I prepared a new interpretive poster. I also add some supportive material. I will further update this report over the next year leading up to the 2024 decadal remembrance.

I remember the earthquake partly because there were two closely spaced earthquakes that slipped the megathrust during a 3 day period.

First was the 1 April 2014 M 8.2 earthquake. This had clearly defined area of aftershocks (I interpreted this area as the area that the fault slipped).

On 3 April there was an M 7.7 earthquake. This caused a sequence of aftershocks that were in a distinct area of the megathrust fault.

I interpreted these two earthquakes to be separate events (in other words, the M 7.7 was not an aftershock of the M 8.2, but a separate (triggered) earthquake (with its own aftershock sequence)).

I remember attending the Seismological Society of America Annual Meeting in Anchorage later in April 2014.

I debated with my esteemed colleague about this pair of earthquakes. They claimed that they thought this was just a single earthquake sequence (i.e., the M 7.7 was an aftershock of the M 8.2).

I discussed the aftershock spatio-temporal pattern, supporting my interpretation. They claimed that the In-SAR results showed that this was a single event.

However, as I was thinking about this later, I realized that the In-SAR data could not distinguish these events as there was no radar acquisition between the 8.2 and the 7.7. So, there was no measurement of the Earth’s surface between the earthquakes.

Needless to say, I tried not to lose sleep that I missed making an excellent point. (there is much more in life to be concerned about, like, is there clean water to drink when I am thirsty.)

I will fill in the gaps for this earthquake over the next year and attempt to do some of this today.

This report is (currently) Chile centric (I will fill in more to the north, in Peru/Ecuador, over the coming year). To learn more about the earthquake history to the north of Arica, Peru, check out the Earthquake Report from early 2023.

I just got back from a USGS Powell Center meeting where we developed the initial basis for Probabilistic Tsunami Hazard Assessment in the Pacific. So, this margin is on my mind right now.

Below is my interpretive poster for this earthquake

  • I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1922-2022 with magnitudes M ≥ 3.0 in one version.
  • I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
  • A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
  • Some basic fundamentals of earthquake geology and plate tectonics can be found on the Earthquake Plate Tectonic Fundamentals page.

    I include some inset figures. Some of the same figures are located in different places on the larger scale map below.

  • In the upper right corner is a map that shows the tectonic plates. I outline aftershock patches from most of the larger earthquakes along this margin over the past century. I include the mechanisms from these events.
  • To the left of the tectonic margin map is a pair of maps that show the landslide probability (left) and the liquefaction susceptibility (right) for this M 6.8 earthquake. I spend more time describing these types of data here. Read more about these maps here.
  • In the lower right corner I plot the USGS modeled intensity (Modified Mercalli Intensity scale, MMI) and the USGS “Did You Feel It?” observations (labeled in yellow). Above the map is a plot showing these same data plotted relative to distance from the earthquake. Read more about what these data sets are and what they represent in the report here.
  • In the upper center I include a figure from Schurr et al., 2014 that shows the earthquake slip areas for the M 8.2 and M 7.7 2014 Iquique earthquakes. They also include the slip area for the 1995 M 8.1. On the left is a time-space diagram that shows the latitudinal range of prehistoric and historic earthquakes, along with the year that those events happened.
  • In the upper left corner is a figure from Metious et al., 2016 that shows (from left to right) the number of earthquakes M>3 per year, the average earthquake slip for 2014, 2015, and 2015, and a map showing the coupling ratio for the megathrust fault.
  • Here is the map with 3 month’s seismicity plotted.

Shaking Intensity

  • Here is a figure that shows a more detailed comparison between the modeled intensity and the reported intensity. Both data use the same color scale, the Modified Mercalli Intensity Scale (MMI). More about this can be found here. The colors and contours on the map are results from the USGS modeled intensity. The DYFI data are plotted as colored dots (color = MMI, diameter = number of reports).
  • In the upper panel is the USGS Did You Feel It reports map, showing reports as colored dots using the MMI color scale. Underlain on this map are colored areas showing the USGS modeled estimate for shaking intensity (MMI scale).
  • In the lower panel is a plot showing MMI intensity (vertical axis) relative to distance from the earthquake (horizontal axis). The models are represented by the green and orange lines. The DYFI data are plotted as light blue dots. The mean and median (different types of “average”) are plotted as orange and purple dots. Note how well the reports fit the green line (the model that represents how MMI works based on quakes in California).
  • Below the lower plot is the USGS MMI Intensity scale, which lists the level of damage for each level of intensity, along with approximate measures of how strongly the ground shakes at these intensities, showing levels in acceleration (Peak Ground Acceleration, PGA) and velocity (Peak Ground Velocity, PGV).

USGS Historic Seismicity

  • Here is a poster that shows the significant earthquakes along this plate boundary. Note how there are earthquakes in the Nazca plate associated with the 2010 and 2015 megathrust subduction zone earthquakes. These are triggered earthquakes along the outer rise, not additional subduction zone earthquakes.
  • In the lower right corner is a figure from Beck (1998) that shows the spatial extent of the known earthquakes. I added the extent of the 2015 and 2010 earthquakes as green arrows.
  • In the upper right corner is an excellent figure from Horton (2018) that shows the plate tectonic setting for this area.

Tsunami

    Potential for Ground Failure

    • Below are a series of maps that show the potential for landslides and liquefaction. These are all USGS data products.
      There are many different ways in which a landslide can be triggered. The first order relations behind slope failure (landslides) is that the “resisting” forces that are preventing slope failure (e.g. the strength of the bedrock or soil) are overcome by the “driving” forces that are pushing this land downwards (e.g. gravity). The ratio of resisting forces to driving forces is called the Factor of Safety (FOS). We can write this ratio like this:

      FOS = Resisting Force / Driving Force

      When FOS > 1, the slope is stable and when FOS < 1, the slope fails and we get a landslide. The illustration below shows these relations. Note how the slope angle α can take part in this ratio (the steeper the slope, the greater impact of the mass of the slope can contribute to driving forces). The real world is more complicated than the simplified illustration below.


      Landslide ground shaking can change the Factor of Safety in several ways that might increase the driving force or decrease the resisting force. Keefer (1984) studied a global data set of earthquake triggered landslides and found that larger earthquakes trigger larger and more numerous landslides across a larger area than do smaller earthquakes. Earthquakes can cause landslides because the seismic waves can cause the driving force to increase (the earthquake motions can “push” the land downwards), leading to a landslide. In addition, ground shaking can change the strength of these earth materials (a form of resisting force) with a process called liquefaction.
      Sediment or soil strength is based upon the ability for sediment particles to push against each other without moving. This is a combination of friction and the forces exerted between these particles. This is loosely what we call the “angle of internal friction.” Liquefaction is a process by which pore pressure increases cause water to push out against the sediment particles so that they are no longer touching.

      An analogy that some may be familiar with relates to a visit to the beach. When one is walking on the wet sand near the shoreline, the sand may hold the weight of our body generally pretty well. However, if we stop and vibrate our feet back and forth, this causes pore pressure to increase and we sink into the sand as the sand liquefies. Or, at least our feet sink into the sand.

      Below is a diagram showing how an increase in pore pressure can push against the sediment particles so that they are not touching any more. This allows the particles to move around and this is why our feet sink in the sand in the analogy above. This is also what changes the strength of earth materials such that a landslide can be triggered.


      Below is a diagram based upon a publication designed to educate the public about landslides and the processes that trigger them (USGS, 2004). Additional background information about landslide types can be found in Highland et al. (2008). There was a variety of landslide types that can be observed surrounding the earthquake region. So, this illustration can help people when they observing the landscape response to the earthquake whether they are using aerial imagery, photos in newspaper or website articles, or videos on social media. Will you be able to locate a landslide scarp or the toe of a landslide? This figure shows a rotational landslide, one where the land rotates along a curvilinear failure surface.


    • Below is the liquefaction susceptibility and landslide probability map (Jessee et al., 2017; Zhu et al., 2017). Please head over to that report for more information about the USGS Ground Failure products (landslides and liquefaction). Basically, earthquakes shake the ground and this ground shaking can cause landslides.
    • I use the same color scheme that the USGS uses on their website. Note how the areas that are more likely to have experienced earthquake induced liquefaction are in the valleys. Learn more about how the USGS prepares these model results here.

      Some Relevant Discussion and Figures

      • Here is the cross section of the subduction zone just to the south of the Sept/Nov 2015 swarm (Melnick et al., 2006). Below I include the text from the Melnick et al. (2006) figure caption as block text.

      • (A) Seismotectonic segments, rupture zones of historical subduction earthquakes, and main tectonic features of the south-central Andean convergent margin. Earthquakes were compiled from Lomnitz (1970, 2004), Kelleher (1972), Comte et al. (1986), Cifuentes (1989), Beck et al. (1998 ), and Campos et al. (2002). Nazca plate and trench are from Bangs and Cande (1997) and Tebbens and Cande (1997). Maximum extension of glaciers is from Rabassa and Clapperton (1990). F.Z.—fracture zone. (B) Regional morphotectonic units, Quaternary faults, and location of the study area. Trench and slope have been interpreted from multibeam bathymetry and seismic-reflection profiles (Reichert et al., 2002). (C) Profile of the offshore Chile margin at ~37°S, indicated by thick stippled line on the map and based on seismic-reflection profiles SO161-24 and ENAP-017. Integrated Seismological experiment in the Southern Andes (ISSA) local network seismicity (Bohm et al., 2002) is shown by dots; focal mechanism is from Bruhn (2003). Updip limit of seismogenic coupling zone from heat-fl ow measurements (Grevemeyer et al., 2003). Basal accretion of trench sediments from sandbox models (Lohrmann, 2002; Glodny et al., 2005). Convergence parameters from Somoza (1998 ).

      • Here is the overview figure from Horton, 2018.

      • Maps of (A) tectonic framework, (B) topography, and (C) sedimentary basin configuration of South America. (A) Map of plate boundaries, Andean magmatic arc (including the northern, central, and southern volcanic zones), regions of flat slab subduction, modern stress orientations from earthquake focal mechanisms, eastern front of Andean fold-thrust belt, and key segments of the retroarc foreland basin system. Plate velocities are shown relative to stable South American plate (DeMets et al., 2010). (B) DEM topographic map showing the Andes mountains and adjacent foreland region, including the Amazon, Parana, Orinoco, and Magdalena (Mag) river systems. (C) Map of Andean retroarc basins, showing isopach thicknesses (in km) of Cretaceous-Cenozoic basin fill, forebulge axis (from Chase et al., 2009), and locations of 13 sites (8 foreland basins, 5 hinterland basins) considered in this synthesis

      • Here is the seismic hazard map is from Rhea et al. (2010).

      • Here is the seismicity map and space time diagram from Métois et al. (2016). The subduction zone fault in the region of Coquimbo, Chile changes geometry, probably because of the Juan Fernandez Ridge (this structure controls the shape of the subduction zone). This figure shows a map and cross section for two parts of the subduction zone (Marot et al., 2014). Note how the subduction zone flattens out with depth here. The M=6.7 quake was shallower than this, but the shape of the downgoing slab does affect the amount of slab pull (tension in the down-dip direction) is exerted along the plate.

      • Left estimated extent of large historical or instrumental ruptures along the Chilean margin adapted from ME´ TOIS et al. (2012). Gray stars mark major intra-slab events. The recent Mw[8 earthquakes are indicated in red. Gray shaded areas correspond to LCZs defined in Fig. 3. Right seismicity recorded by the Centro Sismologico Nacional (CSN) during
        interseismic period, color-coded depending on the event’s depth. Three zones have been defined to avoid including aftershocks and preshocks associated with major events: (1) in North Chile, we plot the seismicity from 2008 to january 2014, i.e., between the Tocopilla and Iquique earthquakes; (2) in Central Chile, we plot the seismicity on the entire 2000–2014 period; (3) in South-Central Chile, we selected events that occurred between 2000 and 2010, i.e., before the Maule earthquake.

      • This figure is the 3 panel figure in the interpretive poster showing how seismicity is distributed along the margin, how historic earthquake slip was distributed, and how the fault may be locked (or slipping) along the megathrust fault.

      • a Histogram depicts the rate of Mw>3 earthquakes registered by the CSN catalog during the interseismic period defined for each zone (see Fig. 2) on the subduction interface, on 0.2° of latitude sliding windows. Stars are swarm-like sequences detected by HOLTKAMP et al. (2011) depending on their occurrence date. Swarms located in the Iquique LCZ and Camarones segment are from RUIZ et al. (2014). Empty squares are significant intraplate earthquakes. b Red curve variations of the average coupling coefficient on the first 60 km of depth calculated on 0.2° of latitude sliding windows for our best model including an Andean sliver motion. Dashed pink curves are alternative models with different smoothing options that fit the data with nRMS better than 2 (see supplementary figure 6): the pink shaded envelope around our best model stands for the variability of the coupling along strike. Green curves coseismic distribution for Maule (VIGNY et al. 2011), Iquique (LAY et al. 2014) and Illapel earthquakes (RUIZ et al. 2016). Gray shaded areas stand for the identified low coupling zones (LCZs). LCZs and high coupling segments are named on the left. The apparent decrease in the average coupling North of 30°S is considered as an artifact of the Andean sliver motion (see Sect. 5.2). c Best coupling distribution obtained inverting for Andean sliver motion and coupling amount simultaneously. The rupture zones for the three major earthquakes are indicated as green ellipses. White shaded areas are zones where we lack resolution.

      • This is a figure that shows details about the coupling compared to some slip models for the 2010, 2014, and 2015 earthquakes.

      • Left coupling maps (color coded) versus coseismic slip distributions (gray shaded contours in cm) for the last three major Chilean earthquakes (epicenters are marked by white stars). From top to bottom Iquique area, white squares are pre-seismic swarm event in the month before the main shock, green star is the 2005, Tarapaca´ intraslab earthquake epicenter, blue star is the Mw 6.7 Iquique aftershock; Illapel area, green squares show the seismicity associated with the 1997 swarm following the Punitaqui intraslab earthquake (green star); Maule area, green star is the epicenter of the 1939 Chillan intraslab earthquake. Right interseismic background seismicity in the shallow part of the subduction zone (shallower than 60 km depth) for each region (red dots) together with 80 and 90 % coupling contours. White dots are events identified as mainshock after a declustering procedure following GARDNER and KNOPOFF (1974). Yellow areas extent of swarm sequences identified by HOLTKAMP et al. (2011) for South and Central Chile, and RUIZ et al. (2014) for North Chile.

      • This is the fault locking figure from Saillard et al. (2017), showing the percent coupling (how much of the plate convergence contributes to deformation of the plate boundary, which may tell us places on the fault that might slip during an earthquake. We are still learning about why this is important and what it means.

      • Comparison between the uplift rates, interseismic coupling, major bathymetric features, and peninsulas along the Andean margin (10°S–40°S). (a) Uplift rates of marine terraces reported in the literature (we present the average rate since terrace abandonment; Table S1 in the supporting information [Jara-Muñoz et al., 2015]). Each color corresponds to a marine terrace assigned to a marine isotopic stage (MIS). Gray dots are the uplift rates of the central Andean rasa estimated from a numerical model of landscape evolution [Melnick, 2016]. (b) Major bathymetric features and peninsulas and pattern of interseismic coupling of the Andean margin from GPS data inversion (this study). Gray shaded areas correspond to the areas where the spatial resolution of inversion is low due to the poor density of GPS observations (see text and supporting information for more details). The Peru-Chile trench (thick black line), the coastline (thin black line), and the convergence direction (black arrows) are indicated. We superimposed the curve obtained by shifting the trench geometry eastward by 110 km (trench-coast distance of 110 km; blue line) with the curve reflecting the 40 km isodepth of the subducting slab (red line; Slab1.0 from Hayes and Wald [2009]), a depth which corresponds approximately with the downdip end of the locked portion of the Andean seismogenic zone (±10 km) [Ruff and Tichelaar, 1996; Khazaradze and Klotz, 2003; Chlieh et al., 2011; Ruegg et al., 2009; Moreno et al., 2011; Métois et al., 2012]. The two curves are spatially similar in the erosive part of the Chile margin (north of 34°S), whereas they diverge along the shallower slab geometry in the accretionary part of the Chile margin (south of 34°S), where the downdip end of the locked zone may be shallower (Figure 4b). Red arrows indicate the low interseismic coupling associated with peninsulas and marine terraces and evidence of aseismic afterslip (after Perfettini et al. [2010] below the Pisco-Nazca Peninsula; Pritchard and Simons [2006], Victor et al. [2011], Shirzaei et al. [2012], Bejar-Pizarro et al. [2013], and Métois et al. [2013] for the Mejillones Peninsula; Métois et al. [2012, 2014] below the Tongoy Peninsula; and Métois et al. [2012] and Lin et al. [2013] for the Arauco Peninsula). FZ: Fracture zone. Horizontal blue bands are the areas where coastline is less than 110 km (light blue) or 90 km (dark blue) from the trench (see Figure 1).

      • Here is the slip patch figure from Schurr et al. (2014).
      • The colored areas show the amount that the fault slipped (0-5 meters) for each earthquake.
      • To the right of the map is a plot showing the density of earthquake energy (moment) for this region.
      • On the left is the historic record of earthquakes for this part of the subduction zone. The area that each earthquake slipped is represented by a vertical bar (labeled by the year of the earthquake).
      • The last time that this part of the subduction zone has had an earthquake that spanned the region of the 2014 event, south to the 1995 earthquake, was in 1877.
      • Over the past few years there have been several papers that document a prehistoric earthquake that generated a trans-pacific tsunami (and a very large local tsunami) in the same area. This earthquake/tsunami from about 3,800 years ago was on the order of the size of the 1960 M 9.5 earthquake.

      • Map of Northern Chile and Southern Peru showing historical earthquakes and instrumentally recorded megathrust ruptures. IPOC instruments used in the present study (BB, broadband; SM, strong motion) are shown as blue symbols. Left: historical1,2 and instrumental earthquake record. Centre: rupture length was calculated using the regression suggested in ref. 28, with grey lines for earthquakes M.7 and red lines for Mw.8. The slip distribution of the 2014 Iquique event and its largest aftershock derived in this study are colour coded, with contour intervals of 0.5 m. The green and black vectors are the observed and modelled horizontal surface displacements of the mainshock. The slip areas of the most recent other large ruptures 4,5,7 are also plotted. Right: moment deficit per kilometre along strike left along the plate boundary after the Iquique event for moment accumulated since 1877, assuming current locking (Fig. 3a). The total accumulated moment since 1877 from 17uS to 25u S (red solid line) is 8.97; the remaining moment after subtracting all earthquake events with Mw 7 (grey dotted line) is 8.91 for the entire northern Chile–southern Peru seismic gap.

      • Here are some figures from Salazar et al. (2022) where they document evidence for an earthquake 3800 years ago in the region of the 2014 earthquake.
      • These authors did extensive research to support their interpretations that included field work (both geological and archeological) and numerical modeling of tsunami waves.
      • Something that is truly remarkable about this earthquake and tsunami from 3800 years ago is the mere size! This 3800a earthquake spanned the megathrust where there were great earthquakes in 1868, 1877, and 1922. This 3.8ka earthquake was similar in size to the 1960 M 9.5 earthquake.
      • This is their map that incudes (from right to left) slip patches from historic earthquakes on the right, coupling ratio to the left of that map, the estimates for tsunami size from historical tsunami, and estimates for local tsunami runup from this 3800 year ago tsunami.

      • Tectonic plates setting, major historical tsunamigenic megathrusts, sites, and tsunami modeling results supporting a Mw ~9.5 earthquake scenario ~3800 cal yr B.P. in northern Chile. (A) Proposed ca. 1000-km-long megathrust rupture for a ~3800 cal yr B.P. tsunamigenic earthquake, encompassing three major segments in front of Iquique, Taltal, and Atacama regions. (B) Sites of geological and archaeological observations for coastal uplift and paleotsunami deposits dated to ~3800 cal yr B.P. and tsunami modeling results. 1: Pabellón de Pica; 2: Paquica Sur; 3: Cobija; 4: Michilla; 5: Playa Grande de Mejillones; 6: Mejillones Peninsula; 7: Zapatero; 8: San Ramón-15 mine; 9: Los Bronces; 10: Hornos de Cal; 11: Playa San Pedro; 12: Bahía Cisne, Caldera. (C) Tsunami run-ups associated with the major historical tsunamigenic earthquakes of 1868, 1877, and 1922 CE in northern Chile (from National Hydrographic and Oceanographic Service). (D) Map showing tectonic plates coupling indicating a high potential for large megathrust earthquakes along northern Chile (9). (E) Seismotectonic setting along northern Chile, with the subduction of the Nazca plate beneath the South American plate (8, 9) and historical megathrusts ruptures (11, 12).

      • This figure shows the numerical ages that form the basis for their interpretation of the timing of this earthquake/tsunami (Salazar et al., 2022).

      • Inferred ages for a tsunamigenic megathrust episode occurred during the Late Holocene, deduced from 12 sites located along northern Chile (see location in Fig. 1). Horizontal bars correspond to the minimum and maximum ages for each calibrated age range. PP, Pabellón de Pica; PS, Paquica Sur archaeological site; CO, Cobija; Mi, Michilla; PG, Playa Grande de Mejillones; MP, Mejillones Peninsula; Za, Zapatero archaeological site; SR, San Ramón archaeological site; LB, Los Bronces archaeological site; HC, Hornos de Cal archaeological site; SP, Playa San Pedro; BC, Bahía Cisnes.

      • This figure shows how changes in culture may have responded to this earthquake. The authors show how cemeteries were placed in different locations (distance from the coast and in elevation) following this period of time.

      • Changes in the archaeological record after the 3800 cal yr B.P. event at regional and local scales. (A) Artefactual and ecofactual density from occupations at the Zapatero site. (B) Artefactual and ecofactual density in residential sites of the Taltal-Paposo area. (C) Total site counts for residential occupations and cemeteries from the Taltal-Paposo area. (D) SPDs of available radiocarbon ages for coastal Antofagasta and Tarapacá regions. (E to G) Change in location of cemeteries between the Archaic IV and Archaic VI periods along the coast of the Antofagasta region.

      • Here is an animation showing the tsunami modeled for the 3.8ka tsunami

    Return to the Earthquake Reports page.

    Eeewww not good at all

  • Leave a Reply

    Your email address will not be published. Required fields are marked *