85

Revista de Geografía Norte Grande, 53: 85-106 (2012)

Tema central

Riesgo de tsunami en caleta Tubul,

Región del Biobío: escenarios extremos y
transformaciones territoriales
posterremoto1
Carolina Martínez2, Octavio Rojas3, Rafael Aránguiz4, Arturo Belmonte5,
Álvaro Altamirano6 y Paulina Flores7

RESUMEN
Se analiza el riesgo de inundación por tsunami para un evento extremo en la
localidad de Tubul (37ºS), Región del Biobío. Para el análisis de peligrosidad se
determinaron tres escenarios asociados a los eventos tsunamigénicos locales de
1835 y 2010 a los cuales se aplicó modelamiento numérico. Se incluyó el evento
de 1877 para determinar los efectos de un tsunami lejano. Se analizó la vulnerabi-
lidad en sus dimensiones física, socioeconómica y educativa. El tsunami de 1835
se determinó como evento extremo, con alturas de inundación de 10 m y cota de
inundación de 10 m. Se establecieron niveles de vulnerabilidad altos y medios
explicados por la deficiente materialidad de la vivienda y perfil socioeconómico
vulnerable de la población, bajos niveles de escolaridad y reacción frente a estos
eventos. El riesgo natural obtenido se ubicó en niveles altos para toda la localidad.
Se contrastan los resultados con el Plan de Reconstrucción vigente.
Palabras clave: Desastres naturales, vulnerabilidad, resiliencia, planificación terri-
torial.

ABSTRACT
The risk of tsunami inundation is analyzed for an extreme event in the town of
Tubul (37°S), Bio-Bio Region. Three scenarios were identified for risk assessment,
in all of which numerical simulation was applied. Two of the three scenarios were
local tsunamis, the events of 1835 and 2010, and the event of 1877 was included
to determine the effects of a far field tsunami. The tsunami of 1835 was determi-
ned as an extreme event, which gave inundation heights of 10m and run-up of
10m. Vulnerability was analyzed from physical, socioeconomic and educational
points of views. Two levels of vulnerability were defined, namely medium and
high. These levels were selected based on the poor housing materials, the vulne-
rable socioeconomic profiles of the population, low educational levels and the
population´s reactions to these tsunami events. The results show that natural risk
is obtained at the high level in the whole town. These results were also compared
with the current Reconstruction Plan.
Key words: Natural Disasters, Vulnerability, Resilience, Territorial Planning

1 Proyecto FONDECYT Nº 11100379. Artículo reci- 4. Departamento de Ingeniería Civil, Universidad Cató-
bido el 30 de agosto de 2012, aceptado el 20 de lica de la Santísima Concepción (Chile).
septiembre de 2012 y corregido el 4 de octubre de E-mail: [email protected]
2012. 5 Departamento de Geofísica, Universidad de Con-
2 Departamento de Geografía, Universidad de Con- cepción (Chile). E-mail: [email protected]
cepción (Chile). E-mail: [email protected] 6 Departamento de Geografía, Universidad de Con-
3 Centro EULA-Chile y Departamento de Geogra- cepción (Chile). E-mail: [email protected]
fía, Universidad de Concepción (Chile). Becario 7 Departamento de Geografía, Universidad de Con-
CONICYT 2012. E-mail: [email protected] cepción (Chile). E-mail: [email protected]
86 REVISTA DE GEOGRAFÍA NORTE GRANDE

Los tsunamis son fenómenos recurrentes cutando un Plan de Reconstrucción que está
en la costa chilena debido a su contexto generando la consolidación residencial y el
tectónico. La mayoría de estos se generan paso rápido de área rural a urbana, el propó-
debido a terremotos de subducción, los cua- sito de esta investigación es analizar el riesgo
les presentan un notable registro histórico de de tsunami para un evento extremo, que per-
efectos devastadores (Lockridge, 1985; Lom- mita generar criterios para el manejo del ries-
nitz, 1970, 1971; Kausel, 1986; Díaz, 1992; go, así como al proceso de Reconstrucción
Urrutia y Lanza, 1993; Monge, 1993; Lorca que vive actualmente la región.
y Recabarren, 1994; Lagos, 2000; Palacios,
2012). Recientemente, el terremoto Mw = 8,8 Materiales y métodos
del 27/F de 2010 y posterior tsunami produ-
jeron efectos devastadores en la costa chilena
Área de estudio
a lo largo de unos 600 km, sin embargo, los
más destructivos ocurrieron en las pequeñas
La localidad de Tubul se localiza en el
localidades costeras de la Región del Biobío,
golfo de Arauco (37º18’S), y desde el punto
la mayoría de estas dedicadas a la actividad
de vista administrativo pertenece a la comuna
pesquera artesanal. Entre estas, destaca la
de Arauco; comprende una superficie apro-
localidad rural de Tubul, una pequeña caleta
ximada de 3 km2 y una población de unos
de unos 2.000 habitantes, que hasta antes
2.000 habitantes, principalmente concentra-
del tsunami tenía como principal actividad da en la ribera del río Tubul (Figura Nº1) y
la explotación del alga gracilaria (pelillo) que vinculada a la pesca, fundamental actividad
crecía a expensas del humedal Tubul-Raqui. económica de la localidad.
La localidad forma parte de un sistema inter-
conectado de cuencas costeras cuyo contacto El área se inserta en un sistema interco-
con el mar ha formado una extensa llanura nectado de cuencas costeras cuya superficie
litoral con comportamiento de marisma sa- conjunta es de 274 km2. El clima es de tipo
lobre, formando así uno de los humedales mediterráneo de influencia oceánica con
costeros más relevantes de la región y que estación seca prolongada y lluvias concentra-
antes del terremoto del 27/F fi guraba entre das en la estación invernal (Devynck, 1970).
los sitios prioritarios para la conservación y El sistema se desarrolla sobre la formación
se proyectaba como futuro sitio RAMSAR. El Tubul de edad Pliocena (Pineda, 1983), com-
terremoto del 27/F generó en esta localidad puesta por depósitos sedimentarios, princi-
un alzamiento costero de 1,4 m que produ- palmente areniscas. En el área se reconocen
jo cambios importantes en el nivel de base tres niveles de terrazas marinas que fluctúan
de los ríos Tubul y Raqui, desecamiento del entre 50 y 100 m. Los procesos eustáticos y
humedal y pérdida consiguiente de la biodi- tectónicos cuaternarios generaron áreas de
versidad, así como cambios morfológicos en relleno y solevantamientos que han influido
el sistema estuarial, el cual extendió su zona la morfología del área.
intermareal en unos 100 m con estados alter-
nantes entre estuario barrera y lagoon costero La interacción entre el régimen pluvial de
(Quezada et al., 2010, Farías et al., 2010; las cuencas costeras y la marea de régimen
Martínez et al., 2011; Vargas et al., 2011). micromareal ha conformado uno de los hu-
medales más importantes de la región debido
Actualmente en la Región del Biobío se a la alta diversidad biológica que lo caracte-
desarrolla un proceso de Reconstrucción riza, condición que se ha ido recuperando a
donde se registran profundas transformacio- dos años de ocurrido el terremoto y tsunami
nes territoriales y que hace prioritario contar del 27/F de 2010. La asociación Sacocornio-
con estudios científicos aplicables a la toma Spartinentum densiflorae son las especies
de decisiones, dado que muchas de las áreas más relevantes adaptadas a las variaciones
arrasadas por el tsunami están hoy en día de salinidad (Stuardo y Valdovinos, 1989;
siendo repobladas en los mismos lugares San Martin et al., 1992; Stuardo et al., 1993;
afectados e incluso con un aumento de la Constabel, 1993; Vásquez, 2009).
densificación del uso residencial. Dado que
esta localidad fue una de las más afectadas
por el tsunami del 27/F, y en ella se está eje-
R IESGO DE TSUNAMI EN CALETA T UBUL , R EGIÓN DEL B IOBÍO : ESCENARIOS 87
EXTREMOS Y TRANSFORMACIONES TERRITORIALES P O S T E R R E M O T O

Figura Nº 1
Localización de Tubul, comuna de Arauco, Región del Biobío, Chile.

Fuente: Elaboración propia.

Procedimientos modelamiento numérico requirió la siguiente

secuencia de actividades:
Para evaluar el riesgo de tsunami en la
localidad de Tubul, se definió como la pro- a) Determinación del escenario extremo y
babilidad de ocurrencia de un peligro na- de la deformación inicial que genera el
tural extremo (tsunami), que puede generar tsunami
potenciales daños y pérdidas en el medio
antrópico, y poner a prueba su capacidad A través del uso de registros históricos y
para responder y/o recuperarse (vulnerabili- catálogos sísmicos, se elaboró una base de
dad) una vez que el peligro se convierte en datos histórica con sismos tsunamigénicos
desastre. De este modo, el riesgo representa para la zona centro-sur de Chile. Los eventos
la plasmación territorial de una actuación considerados corresponden a sismos interpla-
humana poco acorde con las características ca con magnitudes superiores a M = 8,0 que
físico-naturales del territorio donde tiene lu- generaron tsunamis en la costa de la región.
gar (modificado de Rojas y Martínez, 2012, Los eventos seleccionados correspondieron a
en prensa). los de 1835 y 2010. Adicionalmente, se in-
cluyó el evento de 1877 ocurrido en el norte
El análisis del peligro de inundación por de Chile, con el propósito de establecer los
tsunami fue realizado a partir de los resulta- efectos de grandes eventos generados fuera
dos de modelamiento numérico aplicando el de la región centro-sur pero con impacto po-
código NEOWAVE desarrollado en la Uni- tencial en ella, ya que la secuencia histórica
versidad de Hawai (Yamazaki et al., 2010). El de grandes terremotos alojados en la zona
88 REVISTA DE GEOGRAFÍA NORTE GRANDE

de ruptura de dicho sismo, permite asignar sidad registrados en documentos históricos,

un periodo de recurrencia de 121 ± 33 años tales como las observaciones realizadas por
(Comte y Pardo, 1991), por lo tanto, un even- Charles Darwin (Gil, 1945). Se define un lar-
to de similares características podría repetirse go de ruptura de 400 km y un ancho de 150
en un futuro cercano. km, con un desplazamiento interplaca de 8
m. Del mismo modo, Lomnitz (1970) recoge
Los parámetros de falla determinados observaciones sobre el fenómeno que ayudan
para calcular la condición inicial para cada a localizar de manera aproximada el evento.
tsunami (deformación), fueron aquellos que Por ejemplo, menciona un levantamiento en
caracterizan el plano de ruptura y el tama- isla Mocha y en el sector de punta Lavapié,
ño del sismo: locación o hipocentro, Mag- de este modo, la zona de ruptura se puede
nitud Momento Mw, azimuth (az), manteo localizar lo suficientemente al sur como para
(dip), ángulo de deslizamiento (rake), largo generar levantamiento en estas zonas, pero
(L), ancho (W) y deslizamiento interplaca no una subsidencia en la desembocadura
(D). Estos parámetros fueron deducidos del río Maule, que también fue descrito. Los
del conocimiento de las características parámetros sísmicos utilizados para generar
sismo-tectónicas de la región, así como de el tsunami de 1835 se muestran en el Cuadro
la relación existente entre la magnitud-mo- N°1.
mento del sismo, el momento sísmico y las
variables área de ruptura y deslizamiento Para el caso del tsunami del 27/F de 2010,
interplaca. se consideraron estos parámetros a partir del
registro de alturas de inundación, run-up,
Para definir la condición inicial del tsu- tiempos de llegada y dirección de propaga-
nami de 1835 no se dispone de mucha in- ción de las ondas del tsunami, considerando
formación, sin embargo, es posible estimar la existencia de más de una fuente tsunami-
de manera tentativa la fuente sísmica. Las di- génica. En el caso del evento de 2010, se
mensiones de la fuente fueron estimadas por consideró el modelo de ruptura no homogé-
Nishenko (1985) a partir de datos de inten- neo propuesto por el USGS (Hayes, 2010),

Cuadro Nº 1
Parámetros de falla para eventos históricos

Parámetros Tsunami 1835 Tsunami 1877

Origen (borde norte sobre fosa oceánica) 35,2º S - 73,9º W 19,0º S - 72º W

Dislocación 8m 12 m

Largo 400 km 420 km

Ancho 150 km 130 km

Rumbo N15° E N

Buzamiento 18° 19°

Profundidad mínima 5 km 10 km
Angulo
90° 100°
desplazamiento
Magnitud sismo (Escala Richter) 8.8 8.8

Fuente: Elaboración propia a partir de información histórica.

R IESGO DE TSUNAMI EN CALETA T UBUL , R EGIÓN DEL B IOBÍO : ESCENARIOS 89
EXTREMOS Y TRANSFORMACIONES TERRITORIALES P O S T E R R E M O T O

el cual ha sido validado regionalmente por b) Elaboración de grillas de diferente resolu-

Yamazaki & Cheung (2011), usando registros ción espacial
de las boyas DART 32412, 51406, 32411 y
43412. Este modelo de ruptura utiliza 180 fa- Se generaron cuatro grillas anidadas de
llas finitas de largo y ancho fijos, pero con los diferente resolución espacial (2’, 30’’, 6’’ y
demás parámetros sísmicos variables. 1’’) con el propósito de aplicar el modelo
numérico. Las grillas de menor resolución (2’
Para la determinación de los parámetros y 30’’) utilizaron bases de datos batimétricas
sísmicos del evento de 1877 registrado en internacionales tales como General Bathy-
el norte de Chile, se utilizaron las siguientes metric Chart of the Ocean (GEBCO) y 5-Mi-
fuentes principales: mecanismos focales de nute Gridded Global Relief Data Collection
grandes terremotos como los obtenidos para (ETOPO-5), mientras que la grilla de re-
los eventos de 1995 (Antofagasta) y 2001 solución 6” se basó en las cartas náuticas
(Arequipa), a partir de los cuales es posible del Servicio Hidrográfico y Oceanográfico
estimar los parámetros focales de cada sismo, (SHOA), y la de mayor resolución utilizada en
lo que representa los parámetros del plano de la localidad de Tubul, fue creada utilizando
ruptura o plano de Wadati-Benioff (Delouis datos batimétricos, topobatimétricos y topo-
et al., 2010; Pinares, 2006), y los resultados gráficos de detalle obtenidos en terreno entre
del estudio asociado al contacto sismogéni- 2011 y 2012 mediante GPS TRIMBLE doble
co interplaca a lo largo de Chile a partir del frecuencia y ecosonda Garmin acoplado a
análisis del conjunto de mecanismos focales bote pesquero. Estos datos fueron vinculados
asociados a la sismicidad de la región para a vértices creados especialmente para la zona
el periodo 1960-1995 obtenidos de la tesis dado que el terremoto del 27/F de 2010 des-
de grado de Belmonte (1997). Los parámetros truyó la red geodésica regional. Los levanta-
sísmicos utilizados en los sismos de 1835 mientos batimétricos cubrieron los ríos Tubul
y 1877 se muestran en el Cuadro Nº 1. La y Raqui, además del sector marino adyacente
Figura Nº 2 presenta la condición inicial del a la caleta, mientras que los levantamientos
tsunami para cada uno de los escenarios ana- topográficos cubrieron la zona estuarial y el
lizados. humedal hasta la cota de los 20 m. Para co-
tas mayores se utilizaron datos topográficos

Figura Nº2
Condición inicial de los tsunamis 1835, 1877 y 2010

Fuente: Elaboración propia

90 REVISTA DE GEOGRAFÍA NORTE GRANDE

extraídos de imágenes AsterGDem, con reso- modelaciones en Tubul, mientras que la Fi-
lución de 30 m, disponibles en: http://www. gura Nº 3B muestra solo parte de la malla de
gdem.aster.ersdac.or.jp/, generándose contor- mayor resolución en la bahía de Concepción.
nos cada 5 m.
Tal como se mencionó anteriormente,
c) Aplicación del modelo numérico y valida- el modelo de ruptura propuesto por el U.S.
ción de resultados Geological Survey (USGS) para el sismo de
2010 ha sido validado regionalmente en el
La aplicación del modelo de propagación trabajo de Yamazaki & Cheung (2011), por lo
de tsunami consideró los siguientes aspectos: que el presente trabajo realiza una validación
las condiciones iniciales de generación de local utilizando la información de terreno re-
tsunami aplicadas para deformación estática cogida por los autores durante el postsunami
y dinámica; la batimetría y topografía en re- en las diferentes localidades afectadas de la
gión de interés plasmada en cuatro niveles región. Esta información ha quedado registra-
de grillas de tal forma que las menores están da en los trabajos de Quezada et al. (2010),
contenidas en las mayores; la definición de Aránguiz (2010), Martínez et al. (2011). Los
los tamaños de celda en cada grilla; las cons- resultados del modelo numérico en Tubul fue-
tantes físicas del medio y un periodo máximo ron comparados a su vez con las mediciones
de propagación e interacción con el entorno realizadas por Fritz et al. (2011) en trabajos
variable entre 6 horas para los eventos de de terreno efectuados en la zona. La Figura
campo cercano y 10 horas para el tsunami Nº 4 indica las mediciones en la zona en
de campo lejano. Se elaboraron además ma- estudio donde el círculo indica la altura de
reogramas sintéticos para dos sectores de la inundación y el triángulo el run-up. Adicio-
caleta (en la desembocadura del río y frente nalmente, se realizó una modelación numé-
a la playa de Tubul) con el propósito de rica en la bahía de Concepción utilizando la
comparar las alturas de ola y los tiempos de misma resolución de mallas y se comparó el
arribo del tsunami. Del mismo modo, se de- mareógrafo de Talcahuano con un mareógrafo
finió un mareógrafo sintético en Talcahuano, sintético en la misma ubicación (Figura Nº
el cual fue usado para validar los eventos de 3B).
1877 y 2010. La Figura Nº 3 muestra la ubi-
cación de los mareógrafos sintéticos en Tubul Por otro lado, para validar el tsunami de
y Talcahuano. La Figura Nº 3A muestra la 1877 se analizó la serie de tiempo de un ma-
malla de mayor resolución utilizada para las reógrafo virtual en Talcahuano, ya que existen

Figura Nº3
Localización de mareógrafos sintéticos

Fuente: Elaboración propia.

R IESGO DE TSUNAMI EN CALETA T UBUL , R EGIÓN DEL B IOBÍO : ESCENARIOS 91
EXTREMOS Y TRANSFORMACIONES TERRITORIALES P O S T E R R E M O T O

observaciones visuales del tsunami de 1877 ción máxima asociada al evento extremo. Los
en dicha localidad (Soloviev & Go, 1975). niveles de peligrosidad fueron determinados
Según estas observaciones, en Talcahuano se a partir de alturas de flujo de 0, 0.5 y 2.0 m
observó una oscilación importante del mar (Beyer y Aránguiz, 2010).
unas 3 horas después de ocurrido el sismo.
Las observaciones visuales de esa época in- Para determinar áreas vulnerables frente a
dican también que en torno a las 0:30 horas tsunami, se definió la vulnerabilidad como el
(cuatro horas y media después del sismo) el grado de pérdida o grupos de elementos bajo
nivel del mar alcanzó 1,2 metros por debajo riesgo, resultado de la probable ocurrencia de
de la bajamar regular, esto es, unos 1,8 me- un suceso desastroso (UNDRO, 1979; Cardo-
tros aproximadamente, y a las 01:00 horas na, 1993). Bajo este enfoque, la vulnerabili-
(5 horas después del sismo) alcanzó los 2 dad total frente a tsunami (Wilches-Chaux,
metros sobre la marea máxima regular, esto 1988, en Maskrey, 1993) fue determinada
es, unos 2,8 metros sobre el nivel medio del por el estudio de tres vulnerabilidades indi-
mar. Finalmente, el parámetro para medir la viduales: física, socioeconómica y educativa.
intensidad y altura de los tsunamis de 1835 Para su análisis, el poblado fue divido en
y 1877 se realizó a partir de run-up de áreas unidades de análisis (UA) definidas por man-
máximas de inundación en conjunto con va-
zanas derivadas del Plan Regulador Comunal
lores de series de tiempo.
(PRC) de Arauco (Ilustre Municipalidad de
Arauco, 2005); las manzanas con dificultades
A partir de los resultados del modela-
de definición fueron delimitadas con méto-
miento numérico, se elaboró una carta de
dos de terreno mediante GPS cartográfico.
peligro de inundación por tsunami para la
La división generó 29 UA, que incluyen las
localidad de Tubul, utilizando al SIG ArcGis
aldeas provisorias emplazadas postevento
9.3, en la cual se incluye el área de propaga-

Figura Nº 4
Mediciones en Tubul realizadas por Fritz et al. (2011)

Fuente: Elaboración propia.

92 REVISTA DE GEOGRAFÍA NORTE GRANDE

tsunamigénico del 27-F de 2010. La base car- ción; b) actividad económica; c) pobreza. La
tográfica se obtuvo de la Oficina Nacional de vulnerabilidad educativa fue determinada en
Emergencia (ONEMI) Región del Biobío. La función de: a) nivel de enseñanza; b) nivel de
población de interés fue definida por los jefes información sobre las causas del fenómeno; y
de hogares mayores de 18 años residentes en c) reacción frente al fenómeno. Para cada in-
la aldea de Tubul. El número de la población dicador se utilizaron tres niveles de vulnera-
se determinó en función de datos de la Ilus- bilidad: alto, medio y bajo, definidos en fun-
tre Municipalidad de Arauco y se actualizó ción de datos previos. Para determinar áreas
mediante levantamiento en terreno, que con- jerarquizadas en niveles de vulnerabilidad
sistió en censar el número de viviendas poste- total frente a tsunami, se evaluaron los resul-
rior al tsunami del 27-F, debido a la pérdida tados de las tres cartografías de vulnerabilida-
de viviendas y su posterior emplazamiento des específicas, según la siguiente ecuación:
provisorio. CVT=VF+VSE+VE, donde: CVT (cartografía
de vulnerabilidad total); VF (vulnerabilidad
El tamaño muestral se determinó según los física); VSE (vulnerabilidad socioeconómica);
siguientes factores: amplitud de universo (494 VE (vulnerabilidad educativa).
jefes de hogares); nivel de confianza 95%;
precisión (5%) y una proporción del 90%, El riesgo final se obtuvo en función de la
estimada según datos del Ministerio de Plani- ecuación R=P*V, donde R=Riesgo, P=Peligro,
ficación para la Región del Biobío (Larrañaga A=Amenaza. Los factores fueron integrados
y Herrera, 2010). Aplicada la fórmula para en una matriz (Eckert et al., 2012; Jalínek et
poblaciones finitas, se obtuvo un tamaño de al., 2012), de la multiplicación se obtuvie-
muestra de 109 personas, por razones de pér- ron tres niveles de riesgo: alto, medio y bajo
didas o rechazo, finalmente se encuestó a 176 (Cuadro Nº 2), con puntajes desde 1 a 9. El
sujetos. Para el levantamiento de las variables nivel de riesgo se aplicó a las unidades de
incluidas en cada vulnerabilidad, se elaboró análisis mencionadas en el apartado de vul-
una encuesta con preguntas cerradas, estruc- nerabilidad en un entorno de SIG.
turada según las variables de interés de las
vulnerabilidades específicas (Cuadro Nº 2). La
muestra se aplicó en forma aleatoria simple;
Vulnerabilidad de la caleta de
el análisis de los datos se realizó en el soft- Tubul
ware InfoStat (Di Rienzo et al., versión 2011).
Validación de resultados
Para analizar las vulnerabilidades específi-
cas, se utilizó análisis multicriterio (AMC). En a) Tsunami 1877
la vulnerabilidad física las variables de interés
fueron: a) ubicación de la vivienda; b) tipo La Figura Nº 5 indica la serie de tiempo
de vivienda; c) materialidad; d) número de en Talcahuano obtenida para el evento M=8.8
pisos. Para la vulnerabilidad socioeconómica de 1877. Se puede observar una buena co-
se consideró: a) nivel de ingreso de la pobla- rrespondencia con las observaciones descri-

Cuadro Nº 2
Matriz de riesgo de inundación por tsunami, localidad de Tubul

X Peligro
Nivel Bajo (1) Medio (2) Alto (3)
Bajo (1) B 1 X 1 =1 B 1 X 2 =1 M1 X 3 =3
Vulnerabilidad
Medio (2) B 2 X 1 =2 M 2 X 2 =4 A 2 X 3=6
Alto (3) M 3 X 1 =3 A 3 X 2 =6 A 3 X 3 =9

Rango de riesgo: Bajo (1-2), Medio (3-4), Alto (6-9).

Fuente: Elaboración propia.
R IESGO DE TSUNAMI EN CALETA T UBUL , R EGIÓN DEL B IOBÍO : ESCENARIOS 93
EXTREMOS Y TRANSFORMACIONES TERRITORIALES P O S T E R R E M O T O

tas en Soloviev y Go (1975), pues la primera var que tanto la llegada de las dos primeras
onda comienza a aumentar después de 3 ondas como la máxima inundación de la pri-
horas de ocurrido el sismo y luego el nivel mera son bien representadas. Incluso las dos
desciende a poco más de 2 m bajo el nivel oscilaciones menores ocurridas entre las dos
medio del mar. Además, el máximo nivel al- ondas más grandes son captadas por el mode-
canzado coincide tanto en la hora como en lo numérico. En la Figura Nº 7 se presenta el
el valor registrado, esto es, 3 metros, 5 horas área de inundación medida (izquierda) y mo-
después de ocurrido el sismo. delada (derecha) en Tubul para el evento de
2010, se puede observar una concordancia
b) Tsunami 2010 entre ambas áreas de inundación.

La Figura Nº 6 indica la variación del ni- El Cuadro N° 3 presenta una comparación

vel del mar en Talcahuano registrada durante de las mediciones realizadas por Fritz et al.
las dos primeras horas del evento y lo obteni- (2011) (indicadas anteriormente en la Figura
do por el modelo numérico. Se puede obser- Nº 4) y los valores obtenidos por el modelo

Figura Nº 5
Serie de tiempo de mareógrafo virtual en Talcahuano para un evento M=8.8 frente a Iquique

Fuente: Elaboración propia.

Figura Nº 6
Comparación de la variación del nivel del mar medido en Talcahuano y el obtenido en el modelo
numérico

Fuente: Elaboración propia.

94 REVISTA DE GEOGRAFÍA NORTE GRANDE

Figura Nº 7
Área de inundación en Tubul medida (A) y obtenida mediante el modelo numérico (B)

Fuente: A) geoportal.sernageomin.cl/geoportal/rest B) Elaboración propia.

Cuadro Nº 3
Comparación de resultados modelo numérico y Fritz et al. (2011)

E N Medición (m) Modelo numérico (m) Diferencia (m) Tipo

-37,22586 -73,43676 6,2 6,13 0,07 runup
-37,22441 -73,43876 5,4 5,49 -0,09 runup
-37,22556 -73,44132 5,4 4,47 0,93 runup
-37,22904 -73,45708 3,2 3,57 -0,37 runup
-37,24009 -73,42923 8,4 6,50 1,90 runup
-37,22588 -73,43798 5,1 5,29 -0,19 Inundación
-37,22678 -73,43808 6,2 6,22 -0,02 Inundación
-37,23165 -73,43953 5,9 5,22 0,68 Inundación
-37,23018 -73,45706 4,1 3,77 0,33 Inundación

Fuente: Elaboración propia.

numérico. En general, se observa una muy los demás datos se obtuvieron diferencias de
buena correspondencia entre ambos datos, centímetros solamente.
principalmente en la zona plana y desembo-
cadura del río Tubul. Sin embargo, la dife- Propagación del tsunami
rencia más grande se obtiene para el run-up
medido de 8.4 m en la zona este de la des- La Figura Nº 8 indica la serie de tiempo
de los mareógrafos virtuales frente a la playa
embocadura del río Raqui y no en la caleta
de Tubul para los eventos de 1835 y 2010.
de Tubul. Una explicación posible es que el
En general, se puede observar que ambos
sector donde se midió este valor, posee una eventos muestran un comportamiento simi-
topografía en forma de V que genera concen- lar, a pesar de utilizar fuentes tsunamigé-
traciones de energía que el modelo numérico nicas distintas, una con modelo de ruptura
y la resolución de la malla (~30 m) no son homogéneo (1835) y la otra con modelo no
capaces de representar adecuadamente. Para homogéneo con 180 fallas finitas (2010). Se
R IESGO DE TSUNAMI EN CALETA T UBUL , R EGIÓN DEL B IOBÍO : ESCENARIOS 95
EXTREMOS Y TRANSFORMACIONES TERRITORIALES P O S T E R R E M O T O

puede observar que si bien se registra un au- de observar que la primera onda alcanza un
mento del nivel del mar desde el comienzo nivel máximo de 200 minutos (3 horas 20 mi-
del sismo, esta primera onda no es la más nutos) después de ocurrido el sismo. Al igual
destructiva y solo genera un aumento gra- que en el caso anterior, la primera onda no es
dual del nivel del mar, durante los primeros la más grande, sino las ondas posteriores, las
50 minutos, y puede alcanzar los 3 metros cuales alcanzan un nivel máximo de 2 m.
en el caso del evento de 1835. Además, la
mayor inundación se produce con la tercera La Figura Nº 10 indica el área de inunda-
onda, que corresponde a una onda de orilla ción máxima y las velocidades de flujo obte-
que viaja de norte a sur, la cual ingresa por nidas del modelo numérico en la localidad
el extremo norte del golfo de Arauco, la que para los eventos de 1835, 1877 y 2010. Se
alcanza 9 metros para el evento de 1835 y 6 observó que para el evento de 1835, tanto
metros para el de 2010. el área de inundación como la velocidad de
flujo fueron máximas, alcanzando valores
Del mismo modo, la Figura Nº 9 muestra de hasta 10 m y 7m/s, respectivamente. El
la variación del nivel del mar en Tubul para evento de 2010 presentó características muy
el evento de 1877 durante 10 horas. Se pue- similares al de 1835, sin embargo, de efectos

Figura N° 8
Serie de tiempo en el mareógrafo virtual frente a Tubul para los eventos de 1835 y 2010

Fuente: Elaboración propia.

Figura N° 9
Serie de tiempo en el mareógrafo virtual frente a Tubul para el evento de 1877

Fuente: Elaboración propia.

96 REVISTA DE GEOGRAFÍA NORTE GRANDE

Figura Nº 10
Área de inundación por tsunami y velocidades máximas de flujo

Fuente: Elaboración propia.

R IESGO DE TSUNAMI EN CALETA T UBUL , R EGIÓN DEL B IOBÍO : ESCENARIOS 97
EXTREMOS Y TRANSFORMACIONES TERRITORIALES P O S T E R R E M O T O

menores dado que la altura máxima de inun- corresponde a la onda primaria que ingresa
dación alcanzó los 6 m en la zona, mientras al golfo entre la isla Santa María y la penín-
que las velocidades de flujo registraron has- sula de Arauco, con una dirección hacia el
ta 5 m/s. Para el evento lejano de 1877, la este. Además, se puede observar el efecto de
altura de ola máxima fue de 2 m, afectando refracción inducido por la presencia del ca-
principalmente al litoral arenoso del área ñón submarino del Biobío, el cual genera un
entre Tubul y playa Raqui, con velocidades de cambio de dirección en la propagación del
flujo de hasta 2 m/s. La cota de inundación frente principal, que a su vez es incrementa-
para este evento fue de 2,5 m. Hacia el inte- do debido al efecto de difracción de la isla
rior de la localidad, las alturas de ola máxi- Santa María (Aránguiz, 2012). Las Figuras Nº
mas fueron de 1,7 m. En todos los casos se 11C y 11D muestran la propagación del tsu-
observó un efecto amortiguador del cerro isla nami desde el tiempo T=110 minutos a 150
localizado entre los ríos Tubul y Raqui. minutos, tal que se captura la propagación de
la onda de orilla que genera la máxima inun-
La Figura Nº 11 muestra la evolución del dación en Tubul (ver Figura Nº 8). Se puede
tsunami de 1835 a medida que ingresa al observar que esta onda ingresa al golfo de
golfo de Arauco y se propaga hacia Tubul. La Arauco en dirección norte-sur y en dirección
Figura Nº 11A corresponde a los primeros 5 SW a Tubul, lo cual genera un escenario muy
min después de ocurrido el evento, y se pue- desfavorable.
de observar claramente que la primera onda
que se registra en Tubul (ver Figura Nº 8)

Figura Nº 11
Propagación del tsunami, evento 1835

Fuente: Elaboración propia.

98 REVISTA DE GEOGRAFÍA NORTE GRANDE

Figura Nº 12
Alturas de inundación y cortes transversales en la zona de Tubul (37° S), Región del Biobío

Fuente: Elaboración propia.

Alturas de flujo y niveles de peligrosidad Niveles de vulnerabilidad

La Figura Nº 12 presenta el área de inun- La vulnerabilidad física para la localidad

dación correspondiente al evento de 1835, se presentó solo niveles altos y medios. El nivel
indican también dos cortes transversales que medio representó un 58% de la población;
permiten visualizar las alturas de inundación las viviendas construidas se caracterizaron
y alturas de flujo. Se puede observar que las por ser viviendas tipo casa, construidas con
máximas alturas de flujo en la zona de Tubul material de albañilería y/o madera. En nivel
alto se ubicó el 42% de la población, que
superaron los 5 m, pues la elevación del te-
reside en viviendas construidas con madera,
rreno es de 3-4 m, mientras que las alturas de
internit o lata, de solo un piso (Figura Nº 14).
inundación alcanzaron los 10 m.
En ambos casos el 100% de las viviendas se
encuentra bajo la cota de los 10 m.s.n.m.
Niveles de peligrosidad de inundación
por tsunami La vulnerabilidad socioeconómica de la
población de Tubul se concentra en niveles
Considerando el evento de 1835 como medios y altos. El nivel alto alcanza al 33%
evento extremo, los niveles de peligrosidad de la población, caracterizándose por pre-
fueron obtenidos a partir de las alturas de sentar bajos salarios y niveles de pobreza-
flujo y campo de velocidades obtenidas indigencia elevada. La vulnerabilidad media
del modelo numérico. Dado que el área se abarca al 64% de la población. Solamente
comporta por sus características topográficas un 3% de la población presenta vulnerabi-
como una superficie relativamente uniforme, lidad socioeconómica baja. En los niveles
alto y medio predominaron como actividad
las alturas de inundación alcanzaron la cota
económica la extracción de productos del
de los 10 m, por lo cual toda el área registró
mar; mientras que el nivel bajo corresponde
una inundación a ese nivel (> 2 m) asociadas
a personas que laboran en el sector terciario,
a velocidades de hasta 8 m/s.
como profesionales y técnicos.

La superficie con alta peligrosidad corres- Para la vulnerabilidad educativa se distin-

pondió al 66,7% del área, dentro de la cual guieron niveles altos y medios. La vulnerabi-
se ubica la totalidad del asentamiento de Tu- lidad media quedó representada por el 79%
bul (Figura Nº 13). de la población; en su mayoría presentaron
un nivel de escolaridad que no superó la
R IESGO DE TSUNAMI EN CALETA T UBUL , R EGIÓN DEL B IOBÍO : ESCENARIOS 99
EXTREMOS Y TRANSFORMACIONES TERRITORIALES P O S T E R R E M O T O

Figura Nº 13
Niveles de peligrosidad por tsunami, localidad de Tubul

Fuente: Elaboración propia.

enseñanza primaria, su nivel de información a) Áreas de riesgo por tsunami

sobre las causas de un tsunami supera el
40%, y ante la posibilidad de un evento ellos Se determinaron dos niveles de riesgo de
evacuarán a zonas seguras que no han sido inundación por tsunami en la localidad. El
determinadas por instituciones respectivas. nivel alto se registró para el 1,64% del área
En vulnerabilidad alta se encontró el 21% de de estudio, mientras que el nivel medio para
la población; este nivel no difirió estadística- el 65,1% de la localidad. Este último nivel se
mente en cuanto a las características del nivel registró en las áreas actualmente no ocupadas
de instrucción y conocimiento de las causas pero que se localizan dentro del área de peli-
de un tsunami en comparación a la vulnera- grosidad de tsunami. Las áreas con alto riesgo
bilidad alta; sin embargo, no supieron cómo incluyen el asentamiento humano de Tubul y
reaccionar y enfrentar el evento. corresponden a sectores que registraron una
alta vulnerabilidad social sobre áreas de alta
Los niveles de vulnerabilidad total de la peligrosidad de inundación por tsunami (Fi-
localidad fueron medios y altos para el 100% gura Nº 16).
de las unidades de análisis. Estos niveles se
explicaron por la exposición de las viviendas, La localidad de Tubul actualmente expe-
materialidad, nivel de ingreso de la pobla- rimenta un fuerte proceso de urbanización,
ción, escolaridad y reacción frente al fenóme- que de acuerdo al Plan Maestro de Recons-
no; la totalidad de estas variables se situaron trucción considera ocupar las áreas de riesgo
en rangos de vulnerabilidad medios y altos, medio y reocupar a través de la entrega de
explicando la variabilidad de los resultados subsidios para la reconstrucción aquellas vi-
obtenidos (Figura Nº 15). viendas destruidas por el terremoto y tsunami
del 27/F de 2010. La Figura Nº18 presenta el
100 REVISTA DE GEOGRAFÍA NORTE GRANDE

Figura Nº 14
Materialidad de las viviendas, localidad de Tubul

Fuente: Colección personal de los autores.

tipo de viviendas emplazadas dentro de la de la costa con el propósito de inferir los po-
denominada “zona de viviendas antitsuna- tenciales efectos de un tsunami extremo.
mi”. Este Plan considera un parque de mitiga-
ción en el borde costero y la construcción de Algunos rasgos morfológicos mayores
una costanera con ciclovías. tales como cañones submarinos generan
cambios en la propagación del tsunami tal
Conclusiones como ha sido analizado por Aránguiz (2011
y 2012) para el caso del cañón submarino del
Los efectos diferenciales registrados en la Biobío, el cual evidencia un efecto mitigador,
costa del centro sur de Chile, generados por a diferencia de otros cañones submarinos en
el tsunami del 27F de 2010, establecieron el mundo que generan efectos amplificado-
que la configuración de la costa es un aspec- res y mayores áreas de inundación (caso Sri
to clave a considerar en los escenarios extre- Lanka y Blangadesh). Este proceso ya había
mos a modelar. sido advertido por Quezada (2000) a través
del análisis de refracción. De acuerdo con
El tsunami del 27F utilizó los cursos lo- Aránguiz (2011), la localización de la isla
cales de las pequeñas localidades costeras Santa María al interior del golfo de Arauco,
(principalmente paleobahías y estuarios) para generaría procesos de difracción, que en
propagar la energía al interior de estas, ge- conjunto con la refl exión-transmisión en el
nerando una mayor superficie inundada. Las cañón, tenderían a desviar el frente de onda
características geomorfológicas, en especial al lado sur del golfo de Arauco con impacto
de la Región del Biobío, generan condiciones directo en las localidades de Arauco, Tubul y
propicias para recibir los impactos de estos Llico, lo cual coincidió con las alturas de olas
eventos dada su herencia morfogenética don- mayores registradas, donde el run-up alcanzó
de los procesos eustáticos y tectónicos han entre 4 y 13 m de acuerdo con observacio-
construido planicies de playa con depósitos nes locales (Quezada et al., 2010; Fritz et
de cubierta eólica antigua sobre asentamien- al., 2011). Los resultados del modelamiento
tos humanos construidos sobre el antiguo numérico coincidieron con estas mediciones
delta del Biobío (Isla et al., 2012) o campos realizadas en terreno y se ajustaron mejor a
dunares que protegieron algunas ciudades la zona más expuesta de la bahía de Tubul
costeras de efectos aún mayores como lo fue donde se localiza el asentamiento, que en
en Cobquequra, Tirúa y Vegas del Itata. Resul- el sector protegido del río Raqui. Por lo tan-
ta indispensable conocer estas características to, independientemente de la dirección de
R IESGO DE TSUNAMI EN CALETA T UBUL , R EGIÓN DEL B IOBÍO : ESCENARIOS 101
EXTREMOS Y TRANSFORMACIONES TERRITORIALES P O S T E R R E M O T O

Figura Nº 15
Niveles de vulnerabilidad, localidad de Tubul

Fuente: Elaboración propia.

Figura Nº 16
Áreas de riesgo de inundación por tsunami, localidad de Tubul

Fuente: Elaboración propia.

102 REVISTA DE GEOGRAFÍA NORTE GRANDE

Figura Nº 17
Plan de Maestro de Reconstrucción posterremoto 27F 2010, localidad de Tubul

Fuente: Elaboración propia a partir de Plan Maestro de Reconstrucción-Tubul (GORE,

2010)

Figura Nº 18
Tipo de viviendas insertas en el proceso de Reconstrucción en la localidad de Tubul

Fuente: Colección personal de los autores.

procedencia del tren de ondas del tsunami, yores daños en bahías expuestas tales como
estas localidades estarán más expuestas a los Tirúa, Lebu, Tubul, Llico y Dichato, así tam-
efectos de un evento extremo, como ha sido bién las localidades de Tumbes y la isla Santa
demostrado con los resultados del modela- María experimentaron alzamientos costeros,
miento numérico de los tsunamis de 1835 en Tubul se determinó en 2 metros (Cisternas
y 2010, los cuales presentarían similitudes et al., 2010), muy cercano al establecido para
respecto a sus efectos ya que de acuerdo a el terremoto de 2010 de 1,4 metros (Quezada
crónicas históricas el primero produjo los ma- et al., 2010; Farías et al., 2010).
R IESGO DE TSUNAMI EN CALETA T UBUL , R EGIÓN DEL B IOBÍO : ESCENARIOS 103
EXTREMOS Y TRANSFORMACIONES TERRITORIALES P O S T E R R E M O T O

Por otro lado, se sabe que el terremoto del clusivamente por mujeres. Estas actividades,
27F de 2010 generó cambios morfológicos aunque no habían causado impacto en los
drásticos en la costa de la Región del Biobío, indicadores socioeconómicos, dado que en
con alzamientos costeros de 1,4 m en Tubul la localidad el nivel de pobreza ha sido histó-
y 1,7 m en Lebu, situación que afectó a lo- ricamente alto y vinculado a otros problemas
calidades ubicadas al sur del golfo de Arauco sociales, el terremoto y tsunami de 2010
(Quezada et al., 2010), lo que ha generado generó un aumento en las características vul-
una nueva realidad geográfica que debe in- nerables de la población (Altamirano, 2011),
cluirse al momento de efectuar modelamiento que hoy contrastan con una visión de ocupa-
numérico. Esta nueva realidad, en conjunto ción del territorio que está fuera de la manera
con las transformaciones territoriales que idónea de gestionar las áreas de riesgo ya
actualmente experimentan las localidades que se está priorizando el emplazamiento
afectadas por el sismo-tsunami del 27F, hace de infraestructura urbana por sobre las nece-
que permanentemente se deba modelar rea- sidades que vinculan a la población con su
lidades dinámicas que hacen que tanto la espacio natural (arraigo e identidad), que es
peligrosidad como la vulnerabilidad cambien la pesca y el uso de recursos naturales tales
en el tiempo. Tubul es, actualmente, una de como la recolección de algas. La caleta Tubul
las localidades más intervenidas producto en este sentido se localiza en uno de los prin-
de la Reconstrucción y se ha observado, cipales humedales de la región, hoy en día en
solo en el último año, un rápido proceso de lenta recuperación (Valdovinos 2012), y ocu-
urbanización en áreas de alta peligrosidad, pa la antigua planicie de playa construida por
así también se proyecta la construcción de procesos fluviomarinos y tectónicos sobre la
una costanera y ciclovía dentro de la zona cual se desarrolla la marisma salobre, por lo
intermareal, donde eventos comunes en el cual es una zona frágil y que requiere de ma-
área tales como marejadas y olas de tormen- nejo ambiental. El proceso de Reconstrucción
tas tienen la capacidad de producir cambios en este sentido, debe considerar la naturaleza
morfológicos. Uno de los principales factores y fragilidad de los sistemas costeros en donde
que elevan la vulnerabilidad por tsunami en se interviene, con el fin de evitar la degrada-
Tubul es el tipo de vivienda y su materialidad, ción ambiental de estos y consecuencias ne-
localizadas al borde de la zona intermareal gativas ante futuras situaciones de riesgo de
y dispuestas en parcelas irregulares, estas desastre, que como se ha demostrado los fe-
viviendas, construidas en madera, de un piso nómenos de tsunamis son eventos recurrentes
y materiales ligeros, pueden ser masivamente en la costa chilena y en especial en la región
arrasadas cuando el tsunami desarrolla velo- del Biobío (Lagos, 2000; Palacios, 2012).
cidades de corrientes del orden de los 2 m/s
y alturas de flujo o profundidad de la inunda- Los estudios de riesgo hoy en día enfati-
ción sobre los 2 m; así, experiencias mundia- zan en los conceptos de adaptación y resi-
les con eventos extremos han establecido una liencia a partir de las lecciones aprendidas
destrucción total de las viviendas con alturas de eventos extremos en diferentes partes del
de flujo cercanas a 5 m y velocidades de las mundo (Paton & Johnston, 2006; Cutter et al.,
corrientes de 4 m/s (Lagos et al., 2008). Los 2008; Jabareen, 2012). Se hace necesario,
resultados del modelamiento numérico esta- por lo tanto, incluir e incorporar a la forma
blecieron alturas de ola de 9 m para el evento de ocupar el territorio, caracterizado por
de 1835 y 6 m para el de 2010, mientras que una intensa antropización de la zona coste-
las velocidades de las corrientes fueron de ra, la extensa recurrencia histórica de estos
8m/s para el evento de 1835 y de 6m/s para fenómenos y la determinación de escenarios
2010, lo cual elevó el riesgo debido a que la extremos, con el fin de disminuir o controlar
localidad está emplazada en áreas de alta pe- la situación de desastres que estos provocan.
ligrosidad de inundación por tsunami.
Referencias bibliográficas
Por otro lado, la caleta, históricamente,
surge vinculada estrechamente a la actividad ALTAMIRANO, A. Análisis de vulnerabili-
pesquera y en especial a la recolección del dad por tsunami en la localidad de Tubul, Re-
alga gracilaria, esta última explotada casi ex- gión del Bio-Bío, Chile. Concepción: Tesis para
104 REVISTA DE GEOGRAFÍA NORTE GRANDE

optar al Título de Geógrafo, Universidad de le. Concepción: Tesis para optar al Grado de
Concepción, 2011 (Inédito). Doctor en Ciencias Ambientales, Universidad
de Concepción, 1993 (Inédito).
ARÁNGUIZ, R. Modelación numérica del
tsunami de 2010 en la Bahía de Concepción CUTTER, S.; BARNES, L.; BERRY, M.; BUR-
utilizando deformación estática no uniforme. TON, C.; EVANS, E.; TATE, E. & WEBB, J.
Obras y Proyectos, 2010, Nº 8, p. 12-18 A place-based model for a understanding
community resilience. Global Environmental
ARÁNGUIZ, R. Estudio preliminar sobre el Change, 2008, Nº 18, p. 598-606.
efecto del cañón submarino del Biobío en la
propagación de tsunamis. En: Sociedad Chile- DELOUIS, B.; NOCQUET, J. & VALEÉ,
na de Ingeniería Hidraúlica, XX Congreso Chi- M. Slip distribution of the February 27, 2010
leno de Hidráulica, Santiago de Chile, 2011. Mw=8.8 Maule Earthquake, central Chile.
Disponible en Internet: http://cimp.ucsc.cl/wp- From static and high-rate GPS, InSAR and te-
content/uploads/2011/10/Aranguiz.pdf leseismic data. Geophysical Research Letters,
2010, Nº 37, p. L17305.
ARÁNGUIZ, R. Effect of a submarine can-
yon on tsunami propagation in the Gulf of DEVYNCK, J. Contribución al estudio de la
Arauco, Chile. In: International Conference on circulación atmosférica en Chile y el clima de
Coastal Engineering, Santander, 2012. la región del BioBío. Concepción: Memoria de
Título para optar al título de Meteorólogo, Uni-
BELMONTE, A. Análisis del Contacto Sis- versidad de Concepción, 1970 (Inédito).
mogénico Inter-Placa a lo largo de Chile.
Santiago de Chile: Tesis para optar el Título de DÍAZ, J. Estudio de fuentes de tsunamis
Ingeniero Civil Mención Estructuras y Magíster y de terremotos: aplicación en el Norte de
en Ciencias Mención Geofísica, Universidad Chile y Sur de Perú. Valparaíso: Memoria para
de Chile, 1997 (Inédito). optar al Título de Oceanógrafo, Pontificia Uni-
versidad Católica de Valparaíso, Chile, 1992
BEYER, J. y ARÁNGUIZ, R. Estudios de (Inédito).
Riesgos de Sismo y Maremoto para las comu-
nas costeras de la región del Bio-Bío. Con- Di RIENZO, J.; CASANOVES, F.; BALZA-
cepción: Centro de Investigación Marítimo RINI, M.; GONZÁLEZ, L.; TABLADA, M. y
Portuario, Universidad Católica de la Santísima ROBLEDO, C. InfoStat versión 2011. Córdoba:
Concepción (UCSC), 2010 (Inédito). Grupo InfoStat, FCA, Universidad Nacional de
Córdoba, Argentina. Disponible en internet:
CARDONA, O. Evaluación de la amenaza, http://www.infostat.com.ar
la vulnerabilidad y el riesgo: Elementos para el
Ordenamiento y la Planeación del Desarrollo. ECKERT, S.; JELINEK, R.; ZEUG, G. &
La Red, 1993. Disponible en Internet: http:// KRAUSMANN, E. Remote sensing-based as-
www.desenredando.org/public/libros/1993/ sessment of tsunami vulnerability and risk in
ldnsn/html/cap3.htm Alexandria, Egypt. Applied Geography, 2012,
Vol. 32, Nº 2, p. 714-723.
CISTERNAS, M.; MELNICK, D.; ELY, L.;
WESSON, R. & NORAMBUENA, R. Similarities FARÍAS, M.; VARGAS, G.; TASSARA, A.;
between the great Chilean earthquakes of 1835 CARRETIER, S.; BAIZE, S.; MELNICK, D. &
and 2010. American Geophysical Union Chap- BATAILLE, K. Land-level changes produced by
man Conference, 2010, p. 19. the 2010 Mw 8.8 Chile earthquake. Science,
2010, Vol. 329, p. 916.
COMTE, D. & PARDO, M. Reappraisal of
great historical earthquakes in the northern FRITZ, H.; PETROFF, C.; CATALÁN, P.;
Chile and southern Peru seismic gap. Natural CIENFUEGOS, R.; WINCKLER, P.; KALLIGERIS,
Hazards, 1991, Vol. 4, p. 23–44, doi:10.1007/ N.; WEISS, R.; BARRIENTOS, S.; MENESES, G.;
BF00126557. VALDERAS-BERMEJO, C.; EBELIG, C.; PAPA-
CONSTABEL, S. Ambientes sedimentarios DOPOULOS, A.; CONTRERAS, M.; ALMAR,
de los estuarios Tubul y Raqui, VIII Región, Chi- R.; DOMÍNGUEZ, J. & SYNOLAKYS, C. Field
R IESGO DE TSUNAMI EN CALETA T UBUL , R EGIÓN DEL B IOBÍO : ESCENARIOS 105
EXTREMOS Y TRANSFORMACIONES TERRITORIALES P O S T E R R E M O T O

Survey on the 27 February 2010 Chilean Tsu- LAGOS, M. Tsunamis de origen cercano a
nami. Pure and Applied Geophysics, 2011, Nº las costas de Chile. Revista de Geografía Norte
168, p. 1989-2010. Grande, 2000, Nº 27, p. 93-102.

GIL, J. Charles Darwin, viaje de un natu- LAGOS, M.; CISTERNAS, M. y MARDO-

ralista alrededor del mundo. Buenos Aires: El NES, M. Construcción de viviendas sociales en
áreas de riesgo de tsunami. Revista de la Cons-
Ateneo, 1945.
trucción, 2008, Vol. 7, Nº 2, p. 4-16.

GOBIERNO REGIONAL, REGIÓN DEL LARRAÑAGA, O. y HERRERA, R. Encuesta

BÍO-BÍO. Plan de Reconstrucción del Borde Post Terremoto: Principales resultados. Efectos
Costero, PRBC-18. Concepción: Plan Maestro en la calidad de vida de la población afectada
Tubul, 2010. Disponible en Internet: http:// por el terremoto/tsunami. Santiago de Chile:
www.minvu.cl/opensite_20101207193158. MIDEPLAN, 2010.
aspx
LOCKRIDGE, P. Tsunamis in Chile-Perú.
HAYES, G. Finite fault model, updated re- Boulder: Report SE-39 Worl Data Centera for
sult of the feb 27, 2010 Mw 8,8 Maule Chile solid earth and geophysics, 1985.
Earthquake. National earthquake Information
LOMNITZ, C. Major earthquakes and tsu-
Center (NEIC) of United States Geological Sur-
namis in Chile during the period 1535 to 1955.
vey 2010. Disponible en Internet en:
Soderdruck aus der Geologischen Rundschau
http://earthquake.usgs.gov/earthquakes/ Band, 1970, Nº 59, p. 938-690.
eqinthenews/2010/us2010tfan/finite_fault1.php
LOMNITZ, C. Terremotos y Tsunamis en
ILUSTRE MUNICIPALIDAD DE ARAUCO. Chile. Revista Geofísica Panamericana, 1971,
Plan Regulador Comunal, 2005. Concepción. Nº 1, p. 151-178.
Disponible en Internet en: http://www.munia-
rauco.cl/tpma/plan_regulador.html LORCA, E. y RECABARREN, M. Terremotos
y Tsunamis o Maremotos. Santiago de Chile:
ISLA, F.; QUEZADA, J.; MARTÍNEZ, C.; Comisión Oceanográfica Intergubernamental,
FERNÁNDEZ, A. & JAQUE, E. The evolution Centro Internacional de Informaciones, 1994.
of the BíoBío delta and the coastal plains of
MARTÍNEZ, C.; ROJAS, O.; JAQUE, E.;
the Arauco Gulf, BíoBío Region: the Holocene
QUEZADA, J.; VÁSQUEZ, D. y BELMONTE, A.
sea-level curve of Chile. Journal of Coastal Re-
Efectos territoriales del tsunami del 27 de Febre-
search, 2012, Nº 28, p. 102-111. ro de 2010 en la costa de la Región del Bio-Bío.
Revista Geográfica de América Central, Número
JABAREEN, Y. Vulnerability of cities to ex- Especial EGAL, 2011 Costa Rica, p. 1-16.
treme space weather events: A new frontier of
a multidisciplinary urban research. Natural Sci- MASKREY, A. Vulnerabilidad y mitigación
ence, 2012, Vol. 4, Nº 6, p. 368-371. de desastres. En: MASKREY, A (Compilador).
Los Desastres No Son Naturales. La Red, 1993.
JELÍNEK, R.; KRAUSMANN, E.; Disponible en internet: http://www.desenre-
GONZÁLEZ, M.; ÁLVAREZ-GÓMEZ, J.; BIRK- dando.org/public/revistas/dys/rdys03/dys3-1.0-
MANN, J. & WELLE, T. Approaches for tsunami oct-24-2001-TODO.pdf
risk assessment and application to the city of
MONGE, J. Estudios de riesgo de tsunami
Cádiz, Spain. Natural Hazards, 2012, Vol. 60,
en costas chilenas. Jornadas Chilenas de Sis-
Nº 2, p. 273-293. mología e Ingeniería Antisísmica, 1993, Nº 2,
p. 3-22.
KAUSEL, E. Los terremotos de agosto de
1868 y mayo de 1877 que afectaron al sur de NISHENKO, S. Seismic potential for large
Perú y norte de Chile. Boletín de la Academia and great interplate earthquakes along the chi-
Chilena de Ciencias, 1986, Vol. 3, N° 1, p. 8-13. lean and southern peruvian margins of south
106 REVISTA DE GEOGRAFÍA NORTE GRANDE

america: a quantitative reappraisal. Journal STUARDO, J.; VALDOVINOS, C.; FIGUEROA,

of Geophysical Research, 1985, Nº 90, p. R. y OCHIPINTI, A. Los ambientes costeros del
3589–3615. Golfo de Arauco y área adyacente. Concepción
Serie Monografías Científicas, Centro EULA, 1993,
PALACIOS, A. Sismicidad histórica de la p. 31-43.
ciudad de Concepción desde su fundación en
1550 hasta su traslado en 1751. Santiago de STUARDO, J. y VALDOVINOS, C. Estuarios
Chile: Servicio Nacional de Geología y Mine- y lagunas costeras: ecosistemas importantes de
ría, Boletín Nº 64, 2012. Chile central. Revista Ambiente y Desarrollo,
1989, Nº 5, p. 107-115.
PATON, D. & JOHNSTON, D. Disaster Re-
silience, An Integrated Approach. Illinois: Char- SOLOVIEV, S. & GO, Ch. A catalogue of tsu-
les C. Thomas Publisher Ltda., 2006. namis on the Eastern shore of the Pacific Ocean.
Moscow: Nauka Publishing House, 1975.
PINARES, P. Estudio de la distribución es-
pacial y temporal de la sismicidad en el norte UNITED NATIONS DISASTER RELIEF OR-
de Chile. Santiago de Chile: Tesis para optar al GANIZATION (UNDRO). Natural Disasters
título de Ingeniero Civil, Facultad de Ciencias and Vulnerability Analysis, Report of Experts
Físicas y Matemáticas, Universidad de Chile, Group Meeting. Génova, 1979.
2006 (Inédito).
URRUTIA, R. y LANZA, C. Catástrofes en
PINEDA, V. Evolución paleogeográfica de la Chile 1541-1992. Santiago de Chile: Editorial
Península de Arauco durante el Cretácico Supe- La Noria, 1993.
rior-terciario. Concepción: Memoria de Título
para optar al Título de Geólogo, Universidad VARGAS, G.; FARÍAS, M.; CARRETIER, S.;
de Concepción, Chile, 1983 (Inédito). TASSARA, A.; BAIZE, S. & MELNICK, D. Coastal
uplift and tsunami effects associated to the 2010
QUEZADA, J. Peligrosidad de tsunamis en Mw8.8 Maule earthquake in Central Chile. An-
la zona de Concepción. Congreso Geológico dean Geology, 2011, Nº 38, p. 219 -238.
Chileno, N° 9, 2000, p. 92-96.
VALDOVINOS, C.; MUÑOZ, M.; VÁS-
QUEZADA. J.; JAQUE, E.; BELMONTE, A.; QUEZ, D. y OLMOS, V. Desastres naturales
FERNÁNDEZ, A.; VÁSQUEZ, D. y MARTÍNEZ, y biodiversidad: el caso del humedal costero
C. Movimientos cosísmicos verticales y cam- Tubul-Raqui. Revista Sociedad Hoy, 2012, Nº
bios geomorfológicos generados durante el 19, p. 33-52.
terremoto Mw=8,8 del 27 de Febrero de 2010
en el centro-sur de Chile. Revista Geográfica VÁSQUEZ, D. Delimitación y Zonificación
del Sur, 2010, Vol. 1, Nº 2, p. 11-45. de los Humedales Costeros de Lenga y Tubul-
Raqui: Orientaciones para la Planificación Te-
ROJAS, O. y MARTÍNEZ, C. Riesgos Natura- rritorial de la Región del Bío-Bío. Concepción:
les: evolución y modelos conceptuales. Revista Memoria para optar al Título de Geógrafo, Uni-
Universitaria de Geografía, 2012 (en prensa). versidad de Concepción, Chile, 2009 (Inédito).

RUEGG, J.; OLCAY, Y. & LAZO, D. Co, YAMAZAKI, Y.; CHEUNG, K. & KOW-
post and pre-seismic displacements associated ALIK, Z. Depth-integrated, non-hydrostatic
with 1995 July 30 Mw=8.1 Antofagasta, Chile, model with grid nesting for tsunami genera-
earthquake as constrained by InSAR and GPS tion, propagation and run-up. International
observations. Seismological Research Letters, Journal for Numerical Methos in Fluids, 2010,
2011, Nº 72, p. 673-678. DOI:10.1002/fld.2485.

SAN MARTÍN, C.; CONTRERAS, D.; SAN YAMAZAKI, Y. & CHEUNG, K. Shelf re-
MARTÍN, J. y RAMÍREZ, C. Vegetación de sonance and impact of near -field tsunami
las marismas del centro-sur de Chile. Revista generated by the 2010 Chile earthquake.
Chilena de Historia Natural, 1992, Nº 65, p. Geophysical Research Letters, 2011, Vol. 38,
327-342. p. L12605