XIX CONGRESO PERUANO DE GEOLOGÍA

LIMA 2018

DEFORMACIÓN DE LA CORTEZA TERRESTRE EN EL BORDE OCCIDENTAL DEL PERÚ Y

DETERMINACIÓN DEL ACOPLAMIENTO SÍSMICO LA ZONA DE SUBDUCCIÓN DEL PERÚ A
PARTIR DE MEDICIONES GPS/GNSS

Juan Carlos Villegas Lanza1*, Oscar Balladares Purizaga1, Franco Saavedra Pari1, Eliana Mantilla Calisaya1, Wendy
Quiroz Sifuentes1, Rocio Parillo Saico1, Mohamed Chlieh2, Jean-Mathieu Nocquet2, y José Chire3
1
Instituto Geofísico del Perú, Calle Badajoz 169 Urb. Mayorazgo IV Etapa, Ate, Lima, Perú ([email protected])
*
2
Institut de Recherche pour le Developpement, Calle 17 Urb Corpac, San Isidro, Lima, Perú
3
Instituto Geográfico Nacional, Av. Aramburú 1184, Surquillo, Lima, Perú

Resumen cycle. Nowadays GPS geodetic measurements allow us

En las últimas dos décadas las herramientas de to quantify the displacement of a given site on the
geodesia espacial, en particular el Sistema Global de Earth’s surface with a milimetric level of precision.
Posicionamiento por Satélite (GNSS) ha contribuido a During the last decade we have installed some
mejorar el conocimiento de los procesos geofísicos permanent GPS stationn and performed several
asociados a la deformación de la corteza terrestre y al campaign GPS measurements with support of national
ciclo sísmico de los grandes terremotos. Hoy en día las funding and international cooperation. In this study we
mediciones con GPS permiten cuantificar el present the results of the analysis of GPS data in the
desplazamiento de un punto sobre la superficie de la Peruvian subduction zone. The GPS velocity field show
Tierra con un alto nivel de precisión (milimétrico). En different deformation patterns that are result of two
los últimos años mediante el financiamiento de main processes: i) elastic deformation induced by the
proyectos nacionales y de cooperación internacional se stress accumulation along the plate interface, and ii)
ha venido efectuando la instalación de instrumentos de long-term deformation induced by the movement of a
monitoreo permanente y la densificación y medición de new tectonic block located on the western margin of
puntos de control geodésicos temporales. En este Peru. The modeling of GPS displacements allowed us to
estudio se presentan resultados del análisis de datos estimate the inter-seismic coupling along the subduction
GPS adquiridos en los últimos años en el borde plate interface. Our results show that at least 3 areas
occidental del Perú. Los resultados del campo de along the subduction margin of Peru (Lima, Nasca,
velocidad GPS evidencian diferentes patrones de Moquegua-Tacna) have a high inter-seismic coupling,
deformación que principalmente son resultado de dos whose potential to generate large to great magnitude
procesos: i) deformación elástica inducida por la earthquakes (M> 8.0) in the next few years is very high.
acumulación de esfuerzos en la zona de subducción y ii) Keywords: Subduction zone, crustal deformation, GPS,
deformación de largo plazo inducida por el movimiento Earthquake cycle, Peru.
de un nuevo bloque tectónico en el margen occidental
del Perú. El modelado numérico de los desplazamientos 1. Introducción
GPS ha permitido estimar el acoplamiento sísmico en la El Perú se encuentra ubicado en una zona de
interfaz de subducción. Los resultados indican que al convergencia de placas, también conocida como zona de
menos tres áreas a lo largo del margen de subducción subducción, en donde la Placa oceánica de Nazca con
del Perú (Lima, Nasca, Moquegua-Tacna), presentan un frente de más de 6000 km de longitud desde
actualmente alto nivel de acoplamiento sísmico, cuyo Colombia (latitud 6ºN) hasta Chile (latitud 46ºS), se
potencial para generar terremotos de gran magnitud introduce por debajo de la Placa continental
(M>8.0) en los próximos años es muy alto. Sudamericana a una velocidad relativa de 6 a 7 cm/año
Palabras Clave: Zona de Subducción, deformación (Norabuena et al., 1998; Kendrick et al., 2003). Esta
cortical, GPS, ciclo sísmico, terremotos, Perú. zona de subducción es una de regiones sísmicas más
activas de la Tierra, grandes terremotos de magnitud
Abstract M≥8.0, generadores de tsunamis devastadores, ocurren
In the last two decades the space geodesy tools, in con relativa frecuencia cada ~15-20 años.
particular the Global Navigation Satellite System La zona de subducción en el Perú, en función de sus
(GNSS), contributed to improve the comprehension and características geomorfológicas y de la recurrencia de
knowledge of several geophysical processes including grandes terremotos, se divide en 3 segmentos: Sur,
the deformation of the Earth's crust and the earthquake Centro y Norte (delimitados por la dorsal de Nazca y la
XIX CONGRESO PERUANO DE GEOLOGÍA
LIMA 2018

fractura de Mendaña, ver Figura 2 y 3). La historia comprensión de los procesos físicos asociados a la
sísmica de los segmentos Sur y Centro indica que en generación de los grandes terremotos. Hoy en día la
estas regiones han ocurrido periódicamente terremotos comunidad científica, en particular los geofísicos,
de gran magnitud acompañados de tsunamis, con utilizamos las herramientas como el GPS para detectar y
periodos de recurrencia de 150 a 200 años (Dorbath et cuantificar con alto nivel de precisión (1 mm), cualquier
al., 1990), siendo los más recientes los ocurridos en cambio de posición de un punto sobre la superficie de la
Nasca 1996 (Mw7.7), Arequipa 2001 (Mw8.4) y Pisco Tierra y así estudiar cada una de las etapas del ciclo
2007 (Mw8.0). Por el contrario, el segmento Norte se sísmico. Para el estudio y monitoreo de puntos
caracteriza por la ausencia de sismos de gran magnitud geodésicos se utilizan dos tipos de mediciones: 1)
(M≥8.0), al menos así lo evidencian los reportes mediciones temporales o de campaña, es decir se
históricos desde el siglo XVI (Silgado, 1978; Seiner, instalan monumentos geodésicos que son medidos con
2009). No obstante, eventos de magnitud moderada instrumentos GPS cada 1 o 2 años durante sesiones de
como los ocurridos en Piura 1960 (Mw7.6) y Chimbote registro de 48 horas o más, y 2) mediciones continuas o
1996 (Mw7.5), que produjeron tsunamis locales parecen permanentes, para lo cual se establecen redes geodésicas
ser característicos de esta región (Villegas-Lanza, 2014, compuestas por estaciones de monitoreo continuo con
2016). Por otro lado, en el interior del continente tasas de registro de 30 a 1 segundo o menos. Esta
también ocurren sismos de magnitud moderada información resulta de mucha utilidad para cuantificar
(M<7.0), que son producto de la interacción de las fallas los desplazamientos e identificar las áreas de mayor
corticales, y que al ocurrir cerca de la superficie (< 30 deformación asociada a la acumulación de energía
km) suelen ser perjudiciales. sísmica, permitiendo así identificar las áreas de mayor
potencial sísmico en donde ocurrirán los próximos
2. El ciclo sísmico y medición de la deformación terremotos y tsunamis de gran magnitud. Por otro lado,
cortical utilizando datos GPS/GNSS otros países ubicados en zonas de subducción, como
2.1. El ciclo sísmico Japón, Chile Ecuador, vienen implementando este tipo
Los grandes terremotos de subducción en el Perú y el de instrumentación con la finalidad de de complementar
mundo se originan en la interfaz o zona de contacto de y mejorar los sistema de alerta temprana en caso de
placas. Los avances logrados por la sismología y la sismo y tsunami.
geodesia en diferentes regiones del mundo han
permitido caracterizar esta zona como un conjunto de AQUÍ INCLUIR FIGURA 1
áreas o asperezas distribuidas de manera heterogénea,
las cuales impiden el continuo deslizamiento de placas, 3. Datos y procesamiento
ocasionando bloqueo mecánico, acumulación de En este estudio se analizaron datos geodésicos de más
esfuerzos en la interfaz de contacto y deformación de 100 puntos de control geodésicos temporales y
elástica en la placa superior. A este proceso que se da permanentes ubicados en el Perú recolectados durante
con periodos de tiempo de decenas a cientos de años se un periodo de 8 años a partir del 2008. Los datos
le conoce como fase inter-sísmica (Figura 1a). Cuando corresponden a las estaciones de las redes geodésicas
los esfuerzos acumulados en el tiempo sobrepasan cierto instaladas por el IGP, IRD, IGN e IGS (para esta ultima
umbral, estos se liberan súbitamente dando lugar a un red se incluyeron estaciones ubicadas en la placa
sismo. A este proceso de ruptura se le denomina fase Sudamericana y de Nazca). Asimismo se utilizaron las
co-sísmica y suele durar desde algunos segundos hasta tablas de corrección de mareas y datos de efemérides y
varios minutos, dependiendo del tamaño del sismo de orbitas de los satélites proporcionados por el IGS
(Figura 1b). Luego de la ruptura le sigue un periodo de para efectuar las correcciones correspondientes. El
estabilización caracterizado por la ocurrencia de procesamiento y análisis de datos se llevó a cabo
réplicas, deslizamiento asísmico y la respuesta utilizando el software Gamit/Globk (Herring, et al
viscoelástica de la parte inferior de la litósfera. A este 2010), el cual permite obtener soluciones de posiciones
último proceso se le conoce como fase post-sísmica y diarias con un alto nivel de precisión dentro de un
puede durar desde algunos meses hasta varias decenas marco de referencia determinado. En este estudio se
de años dependiendo del tamaño del sismo. Estas tres utilizó el marco de referencia del ITRF2008.
fases principales, que se repiten de manera sistemática
entre un sismo y otro de similares características, son las 4. Resultados
que definen el ciclo de los grandes terremotos. 4.1. Campo de velocidad GPS y deformación actual de
la corteza terrestre en el Perú
2.2. Medición de la deformación de la corteza terrestre En la Figura 2 se presentan los resultados obtenidos
mediante técnicas de Geodesia Espacial GPS/GNSS del campo de velocidad GPS. Las velocidades,
En las últimas décadas las técnicas de geodesia expresadas en milímetros por año en el marco de
espacial han contribuido a un mejor conocimiento y referencia de Sudamérica estable, muestran distintos
XIX CONGRESO PERUANO DE GEOLOGÍA
LIMA 2018

patrones de deformación que varían a lo largo del de Perú, que abarca desde la provincia de Ilo, Tacna
margen peruano: (1) En el noroeste del Perú las hasta el Norte de Chile (>150 km).
velocidades GPS muestran un desplazamiento cuasi-
constante de 5-6 mm/año en dirección sureste; 5. Discusión de los resultados
comportamiento que es consistente con la hipótesis del Los resultados obtenidos son consistentes con la
movimiento rígido de un bloque tectónico que historia sísmica del margen peruano. En las áreas con
comprendería la zona de ante-arco desde la fosa hasta la mayor acoplamiento sísmico, desde hace varias décadas
cordillera occidental extendiéndose hacia el sur del a cientos de años no se ha producido un terremoto de
Perú. Dicha hipótesis ha sido puesta en evidencia y gran magnitud, por tanto, la acumulación de esfuerzos
discutida ampliamente en los estudios de Nocquet et al., en dichas áreas se sigue incrementando, lo que eleva el
(2014) y Villegas-Lanza, (2014, 2016). (2) En la región peligro sísmico en la costa peruana. Otro resultado
central de Perú se observan altas velocidades GPS con interesante de este estudio indica que el segmento de
dirección noreste-este, similar al sentido de la subducción del norte del país presenta un nivel de
convergencia. Estas velocidades que alcanzan tasas de acoplamiento sísmico débil a moderado, lo que sugiere
hasta ~21 mm/año en la costa, disminuyen que la convergencia de placas en este sector se daría de
sistemáticamente hacia los Andes y el interior del manera asísmica, es decir sin que se produzca
continente, observándose además variaciones laterales importante acumulación de energía como para generar
al norte y sur de Lima, lo cual refleja cambios laterales terremotos de gran magnitud (M>8.5). No obstante,
de un fuerte acoplamiento sísmico en la interfaz de dado que se han registrado sismos de magnitud
subducción de dicha región. (3) En el sur del Perú, a moderada como los ocurridos en 1996 y 1960 (M ≈7.5)
partir de la dorsal de Nazca, se observan un es posible que este tipo de eventos puedan ocurrir
comportamiento heterogéneo, con áreas que tienen altos nuevamente.
valores de velocidades GPS entre 12 y 18 mm/año en la
costa, lo cual sugiere que el acoplamiento sísmico en 6. Conclusiones
esta región es alto, mientras que otras áreas como la En las últimas décadas las observaciones geodésicas
ubicada frente al área de ruptura del sismo del 2001 en zonas de alta actividad sísmica han contribuido a
(Mw8.4) presentan un acoplamiento menor. En general. mejorar el conocimiento sobre los procesos geofísicos
el campo de velocidad observado refleja el efecto asociados a la deformación cortical y al ciclo sísmico de
superpuesto de dos contribuciones principales: (i) el los grandes terremotos.
fuerte acoplamiento heterogéneo a lo largo de la interfaz En este estudio se obtuvieron dos importantes
de subducción, y (ii) deformación tectónica de largo resultados: (1) El campo de velocidad GPS a escala
plazo inducida por el bloque tectónico y el acortamiento nacional, que muestra los actuales patrones de
cortical en la placa continental. deformación de la corteza terrestre, los cuales son
AQUÍ INCLUIR FIGURA 2 producto de la deformación elástica inducida por la
acumulación de esfuerzos en la interfaz de subducción y
4.2. Acoplamiento inster-sísmico en la interfaz de deformación tectónica de largo plazo; y (2) El
contacto de la zona de subducción del Perú acoplamiento sísmico a lo largo de la interfaz o zona de
Para reproducir el efecto de la convergencia en la subducción del Perú, el cual permite identificar tres
interfaz de subducción, es necesario corregir el campo importantes áreas en donde el acoplamiento sísmico es
de velocidad de la contribución del movimiento del alto, y en donde desde hace muchos años no ocurre
bloque tectónico (Nocquet et al 2014, Villegas-Lanza, sismos importantes, lo que incrementa el potencial para
2014, 2016), de modo que se pueda modelar únicamente generar terremotos de gran magnitud (M>8.0) en los
el acoplamiento sísmico. Para el modelado del próximos años. Estas áreas, de norte a sur, son: la región
acoplamiento sísmico se llevó a cabo la la inversión de central de Perú, que abarca el segmento desde Barranca
los desplazamientos GPS utilizando la metodología de hasta Pisco, el segmento frente a las ciudades de Nasca
Okada (1985). Los resultados del modelado han y Chala, y la región sur de Perú, que abarca desde la
permitido identificar tres importantes áreas que provincia de Ilo, Tacna hasta el Norte de Chile.
actualmente presentan alto acoplamiento sísmico. Estas Finalmente, estos resultados constituyen nuevo
áreas están fuertemente acopladas acumulando energía conocimiento científico que permite caracterizar mejor
desde hace muchos años e incrementando su potencial las áreas de mayor peligro sísmico a lo largo del margen
para generar sismos de gran magnitud (M>8.0). Estas de subducción del Perú, y también sirve de base para la
áreas, de norte a sur, son: (1) la región central de Perú, toma de acciones dentro del ámbito de la Gestión de
que abarca el segmento desde Barranca hasta Pisco Riesgo de Desastres, ya que permite priorizar zonas de
(>350 km de longitud), (2) el segmento frente a las intervención para estimar el riesgo y tomar acciones de
ciudades de Nasca y Chala (~150km), y (3) la región sur reducción de la vulnerabilidad y de preparación ante la
ocurrencia de sismos. Dada la importancia para nuestro
XIX CONGRESO PERUANO DE GEOLOGÍA
LIMA 2018

país, se recomienda implementar un mayor número de

instrumentos geodésicos y geofísicos de monitoreo
permanente a fin de mejorar el conocimiento y estar
mejor preparados.

Agradecimientos
A las autoridades y pobladores de las ciudades y
centros poblados donde de realizaron las mediciones de
los puntos de control geodésicos por brindar su apoyo
para la realización de los trabajos. Al Fondo Nacional
de Desarrollo Científico, Tecnológico y de Innovación
por brindar oportunidades para el desarrollo de
proyectos de investigación. Parte de los resultados de
este estudio Figura 1. Esquema de las principales fases del Ciclo
Sísmico en una zona de subducción. a) En la fase inter-
sísmica se produce acoplamiento sísmico y acumulación
Referencias de energía. b) En la fase co-sísmica los esfuerzos
acumulados se liberan y se produce el sismo y tsunami.
Dorbath, L., A. Cisternas, and C. Dorbath (1990),
Assessment of the size of large and great historical
earthquakes in Peru, Bull. Seismol. Soc. Am., 80(3),
551–576.
Kendrick et al., (2003). The Nazca-South America Euler
vector and its rate of change. J.S.A.E.S, 16, 125-131.
Nocquet, J.-M. et al. (2014), Motion of continental
slivers and creeping subduction in the northern Andes,
Nat. Geosci., 7(March), 287–291,
doi:10.1038/NGEO2099
Norabuena, E. et al. (1999), Space Geodetic
Observations of Nazca-South America Convergence
Across the Central Andes, Science, 279(5349), 358–
362, doi:10.1126/science.279.5349.358.
Okada, (1985). Surface deformation due to shear and
tensile faults in a halfspace. Bull. of the Seism. Soc.
Am., vol. 75, 1135-1154.
Seiner Lizárraga, L. (2009), Historia de los sismos en el
Perú. Catálogo: Siglos XV-XVII, edited by F.
Editorial, Universidad de Lima, Lima.
Tavera, H., & E. Buforn (1998), Sismicidad y
sismotectónica del Perú, Fis. Tierra, 10, 187–219.
Silgado, E. (1978), Historia de los sismos más notables
ocurridos en el Perú (1513-1974), Lima.
Villegas-Lanza, J. C. (2014), Earthquake cycle and
continental deformation along the Peruvian subduction
zone, PhD Diss. Univ. Nice. S. Antipolis, France. Figura 2. Se muestra el campo de velocidad GPS en el
Villegas-Lanza, J. C. et al. (2016), Active tectonics of borde occidental del Perú expresado en el marco de
Peru: Heterogeneous interseismic coupling along the referencia de Sudamérica estable. Los vectores
Nazca megathrust, rigid motion of the Peruvian Sliver, representan la velocidad de deformación o
and Subandean shortening accommodation, J. desplazamiento de la corteza terrestre en milímetros por
Geophys. Res. Solid Earth,121, 7371-7394, doi: año. En general se observa que los vectores presentan
10.1002/2016JB013080. desplazamientos hacia el NEE consistentes con la
dirección de convergencia de placas, lo cual refleja
Ilustraciones acoplamiento en la interfaz de placas. Se observa que
las tasas más altas de deformación están en las regiones
el centro y sur mientras que en el norte del Perú son mas
bajas.
XIX CONGRESO PERUANO DE GEOLOGÍA
LIMA 2018