Presentación Tesis InSAR

UNIVERSIDAD TÉCNICA FEDERICO SANTA MARÍA
Academia de Ciencias Aeronáuticas
Interferometría SAR aplicada al estudio del movimiento en

laderas aledañas al volcán Calbuco con ayuda de imágenes
Sentinel-1A
Autor: David Hermosilla Diaz

Profesor/a: Roxana Trujillo
Rodrigo Suárez
Marcelo Durán
Estructura
Área de estudio
Objetivos de Investigación
Fundamento teórico
Desarrollo Metodológico
Resultados
Propuesta de Implementación
Conclusiones
Santiago
Universidad Técnica Federico Santa María 07-10-2016
Área de Estudio
Laderas y proximidades del volcán Calbuco.
Ubicado en provincia de Llanquihue en la X Región de Los Lagos.
• Altura: 2003 m s. n. m.
• Estratovolcán.
• Cercano a volcán Osorno,
Puntiagudo y Puyehue.
• 4ta posición en nivel de
peligrosidad.
• Cono truncado de perfil
irregular.
• Última erupción en 1961,
Sub-Plininana.
Fuente: Google Earth
Santiago
Objetivos de Investigación
Objetivo General:
Utilizar técnica InSAR en las proximidades del volcán Calbuco con la ayuda del procesamiento de
imágenes provenientes del satélite Sentinel-1A, para determinar la deformación asociada al
suceso eruptivo que tuvo lugar durante el mes de Abril del 2015 en la región de Los Lagos.
Objetivos Específicos:
• Estudiar teoría física-matemática que sustenta la técnica InSAR, para aplicarla al estudio y entendimiento
de la deformación en superficies terrestres.
• Emplear herramientas y softwares especializados con el propósito de poder estimar efectos de
subsidencia o levantamiento del área en estudio.
• Comparar los datos obtenidos mediante tecnología InSAR con datos provenientes de otras instituciones
como el SERNAGEOMIN con el fin de validar la veracidad de los resultados.
• Idear una propuesta de implementación que integre la técnica de InSAR a la red de monitoreo volcánico
del SERNAGEOMIN, reforzando así la observación de la actividad volcánica en Chile.
Santiago
Fundamento Teórico
Radar de Apertura Sintética (SAR)
• Radar
- Ofrece su propia iluminación para la toma de imágenes (Teledetección activa):
Funciona de día y de noche, de forma independiente de la luz solar.
- Radiación en el espectro de las microondas:
Es posible atravesar nubes, tempestades y dependiendo de la banda utilizada, áreas foliares o nieve.
a
Nombre Banda Ancho de banda 𝝀(Cm) Frecuencia f(GHz)
K 0,83-2,75 36-10,9
X 2,75-5,21 10,9-5,75
C 5,21-7,69 5,75-3,9
S 7,69-19,4 3,9-1,55
L 19,4-76,9 1,55-0,39
P 76,9-133 0,39-0,225
Fuente: DLR. Fuente: Elaboración propia.
Santiago
Fundamento Teórico
• Concepto básico Radar
1.- El Radar emite un pulso de energía
(microondas) hacia el suelo.
1.- Lapulso
2.- El antena emite en
es esparcido untodas
pulso
las de energía (microondas)
direcciones. hacia el suelo.
2.- Una vez que el pulso toca el suelo es esparcido en todas las direcciones.
3.- 3.- Parte de este pulso retorna en dirección de la
antena receptora del radar (Retro-dispersión)
4.- Este pulso es captado en función de su

amplitud, fase, tiempo de retorno y polarización.
5.- Estos parámetros son procesados en conjunto, Fuente: Centre for Remote Imagin, Sensing & Processing (CRISP).
para formar una imagen de la superficie capturada.
Santiago
Fundamento Teórico
• Funcionamiento básico SAR
- “SAR es un proceso en donde se almacenan todas las señales
retornadas, con sus respectivas amplitudes y fases durante el
tiempo T, desde la posición A a la D. Como el cono de
observación se desplaza a lo largo de la trayectoria del sensor,
una apertura sintética es generada por esta señal. Cuando T es
grande, la apertura de la antena también lo es, logrando una
resolución mayor” (Christian Wolff, 2008)
Fuente: Carlos López

Fuente: Christian Wolff.
Santiago
Fundamento Teórico
• Resolución de sensores SAR
-En Alcance o Rango:
𝑐
𝜌𝐴 ≈ 2𝐵
𝑐: Vel. Luz en 𝑚/𝑠
B: Ancho de banda en 𝑀𝐻𝑧
-En Azimut:
𝐿𝑎
𝜌𝐴𝑧 = 2
𝐿𝑎: Longitud física de Antena
Resolución en azimut para sensor RAR:

𝑅𝜆 Fuente: Traducido desde trabajo de Felipe Geremia Nievinski.
𝜌𝐴𝑧 =
𝐿𝑎
R: Distancia Sensor-Objeto
Santiago
Fundamento Teórico
La imagen SAR:
• Componente Amplitud
-Intensidad de la retro-dispersión de la superficie, se observa en escala de grises.
-Depende de la rugosidad, composición química y propiedades dieléctricas del terreno.
• Componente Fase
-Ondas electromagnéticas que describen un movimiento sinusoidal.
Cada vez que la partícula recorre una distancia igual

a la longitud de onda 𝜆, se dice que recorrió un ciclo
igual a 2𝜋.
2𝑅
𝜙=− 2𝜋
𝜆
Fuente: ESA.
Santiago
Amplitud en imagen SAR
Santiago
Fundamento Teórico
Interferometría con Radar de Apertura Sintética (InSAR)
• Interferometría: “Medición de la Interferencia” entre las señales de
Imagen 1 (Main) e Imagen 2 (Slave).
La diferencia de fase entre dos pixeles vecinos de dos imágenes es igual a:
4𝜋 𝐵⊥ 𝑠 4𝜋 𝐵⊥ ℎ 4𝜋
Δ𝜙 =− − + 𝑑 + 𝜙 𝑟𝑢𝑖𝑑𝑜+𝑎𝑡𝑚
𝜆 𝑅 𝑡𝑎𝑛𝜃 𝜆 𝑅 𝑠𝑒𝑛𝜃 𝜆
𝐷𝑖𝑓.𝑑𝑒 𝑓𝑎𝑠𝑒 𝑡𝑜𝑡𝑎𝑙
𝑐𝑜𝑚𝑝.𝑜𝑏𝑙𝑖𝑐𝑢𝑎 𝑐𝑜𝑚𝑝.𝑎𝑙𝑡𝑢𝑟𝑎 𝑚𝑜𝑣
𝐵⊥ : Línea de Base perpendicular entre ambas capturas

s: Resolución en rango
R: Distancia sensor-objeto
𝜃: Ángulo de incidencia
h: Diferencia de altura
d: Movimiento relativo en dirección Línea de Visada (LOS)
Fuente: InSAR Imaging of Aleutian Volcanoes from Space.
Santiago
Fundamento Teórico
• Fase Interferométrica está compuesta de tres factores elementales:
4𝜋 𝐵⊥ 𝑠 4𝜋 𝐵⊥ ℎ 4𝜋
Δ𝜙 =− − + 𝑑 + 𝜙 𝑟𝑢𝑖𝑑𝑜+𝑎𝑡𝑚
𝜆 𝑅 𝑡𝑎𝑛𝜃 𝜆 𝑅 𝑠𝑒𝑛𝜃 𝜆
𝐷𝑖𝑓.𝑑𝑒 𝑓𝑎𝑠𝑒 𝑡𝑜𝑡𝑎𝑙
𝑐𝑜𝑚𝑝.𝑜𝑏𝑙𝑖𝑐𝑢𝑎 𝑐𝑜𝑚𝑝.𝑎𝑙𝑡𝑢𝑟𝑎 𝑚𝑜𝑣
4𝜋 𝐵⊥ ℎ 4𝜋
1.- Δ𝜙 =− + 𝑑 + 𝜙 𝑟𝑢𝑖𝑑𝑜+𝑎𝑡𝑚 Creación DEM
𝜆 𝑅 𝑠𝑒𝑛𝜃 𝜆
𝐷𝑖𝑓.𝑑𝑒 𝑓𝑎𝑠𝑒 𝑡𝑜𝑡𝑎𝑙 𝑡𝑜𝑝𝑜𝑔𝑟𝑎𝑓í𝑎 𝑚𝑜𝑣
4𝜋 𝐵⊥ ℎ 4𝜋
2.- Δ𝜙 = + 𝑑 + 𝜙 𝑟𝑢𝑖𝑑𝑜+𝑎𝑡𝑚 Movimiento
𝜆 𝑅 𝑠𝑒𝑛𝜃 𝜆
𝐷𝑖𝑓.𝑑𝑒 𝑓𝑎𝑠𝑒 𝑡𝑜𝑡𝑎𝑙 𝐷𝐸𝑀 𝑚𝑜𝑣
Δ𝜙𝐷𝐼𝑛𝑆𝐴𝑅 = Δ𝜙𝑡𝑜𝑡𝑎𝑙 − 𝜙𝑡𝑜𝑝𝑜𝑔𝑟𝑎𝑓𝑖𝑎

𝑀𝑜𝑣 𝐷𝐸𝑀
Santiago
Fundamento Teórico
• Interferometría Diferencial SAR (DInSAR):
Δ𝜙𝐷𝑖𝑛𝑆𝐴𝑅
𝑀𝑜𝑣 Movimiento
= Δ𝜙𝑡𝑜𝑡𝑎𝑙 − 𝜙𝑡𝑜𝑝𝑜𝑔𝑟𝑎𝑓𝑖𝑎
𝐷𝐸𝑀
Al realizar la diferencia entre la fase

interferométrica total y la altura
topográfica proporcionada por un DEM,
se obtiene la fase interferométrica
diferencial equivalente al movimiento
superficial en LOS.
Fuente: NASA
Santiago
Fundamento Teórico
• Interpretación del Interferograma:
1.- Observar patrón de cambio en colores de
anillos desde afuera hacia adentro.
2.- Si el patrón se mueve en dirección opuesta al
cambio de fase, indica que la superficie se acerca
hacia el sensor, revelando un Levantamiento.
3.- Si los colores se mueven en la misma
dirección que el cambio de fase, indica que la
superficie se aleja del sensor, revelando un
hundimiento de la zona o subsidencia.
4.- Cada ciclo completado representa 𝜆 2 en LOS.
→
FUENTE: Zhong Lu, Daniel Dzurisin. “InSAR Imaging of Aleutian Volcanoes”.
Santiago
Fundamento Teórico
Naturaleza de los Volcanes
• Distribución de la actividad volcánica:
-Las erupciones volcánicas son maneras que
posee la tierra para enfriarse y transmitir la
energía y el calor interno a la superficie
terrestre.
-La mayor actividad volcánica del mundo se ha

registrado en el Cinturón de Fuego.
-Se debe al movimiento constante de las

placas tectónicas.
Fuente: Infobae.
Santiago
Fundamento Teórico
Índice de Explosividad Volcánica (VEI)
Volumen material
VEI Clasificación Descripción Altura Pluma Frecuencia Ejemplo
volcánico
0 Hawaiano No-Explosiva <100 mts. >1.000 m³ Diariamente 1983 Kilauea, HAW
Hawaiano/
1 Ligera 100-1000 mts. >10.000 m³ Diariamente 2002 Nyiarongongo, CON
Estromboliano
Estramboliano/
2 Explosiva 1 – 5 Km >1.000.000 m³ Semanalmente 2015 Villarrica, CHI
Vulcaniano
3 Vulcaniano Violenta 3 – 15 Km >10.000.000 m³ Anualmente 1985 Nev del Ruiz, COL
Vulcaniano/
4 Cataclísmica 10 – 25 Km >100.000.000 m³ Cada 10 años 1982 Galunggung, IND
Pliniano
5 Pliniano Paroxística >25 Km >1 km³ Cada 100 años 2015 Calbuco, CHI
Pliniano/
6 Colosal >25 Km >10 km³ Cada 100 años 1883 Krakatoa, IND
Ultra-Pliniano
7 Ultra-Pliniano Mega-Colosal >25 Km >100 km³ Cada 1.000 años 1815 Mont Tambora, IND
640.000 AP Yellowstone
8 Ultra-Pliniano Apocalíptica >25 Km >1000 km³ Cada 10.000 años
EEUU
Fuente: Elaboración propia
Santiago
Fundamento Teórico
• Peligros asociados a una erupción:
Tefra: Todo fragmento de lava o roca volcánica de
cualquier tamaño arrojado al aire por las explosiones de
gases calientes de una erupción.
Flujos piroclásticos: Colapso gravitacional de cenizas,
fragmentos de roca y gases volcánicos calientes que
descienden por las laderas de un volcán en forma de
avalanchas.
Lahares: Mezcla de materiales volcánicos sólidos, que
mezclados con agua caen por efecto de la gravedad,
descendiendo a gran velocidad por su naturaleza
líquida.
Lava: Masas de roca fundida que se desparraman por
la superficie terrestre.
Fuente: USGS
Santiago
Santiago
Fundamento Teórico
• Peligros asociados a la erupción de Calbuco en Abril 2015:
-Columna eruptiva que alcanzó los 15 Km de altura y 240 millones de metros cúbicos.
-Emisión de fragmentos de material incandescente.
-Emisión de flujos de lava.
-Dispersión de material piroclástico del tipo ceniza y pómez.
-Hasta 1 metro de cenizas se detectó en Ensenada.
-Lahares que viajaron hasta una distancia cercana a los 15 kilómetros.
-Cambios morfológicos en las laderas del volcán.
Santiago
Desarrollo Metodológico
1.- Descarga Par Interferométrico
Santiago
2.- Sentinel-1 Toolbox (S1TBX) y apertura de imágenes SAR
Santiago
3.- Corregistro
Santiago
4.- Deburts
Santiago
5.- Merge
Santiago
6.- Subset y Fase Interferométrica
Santiago
7.- Remoción de fase Topográfica (DEM SRTM3ArcSEC)
Santiago
8.- Filtro de Goldstein
Santiago
9.- Terrain Correction (UTM, WGS84 Huso 19)
Santiago
9.- Vizualización en Google Earth (Lat/Long WGS84)
Santiago
Resultados
En base a la regla interpretativa para Interferogramas ya descrito, se observa:
• Los colores observados describen un patrón Azul-Naranjo-Amarillo hacia el interior de los anillos.
• Existe alargamiento de Fase (-2,71 → 2,67).
• Avistamiento de 4 anillos hacia el interior.
Para Sentinel-1ª, cuya longitud de onda es λ~5,6:
• La erupción del Volcán Calbuco en Abril de 2015 provocó
una subsidencia de ~12cm en dirección Línea de Visada
(Line of Sight) en el flanco Sur-Oeste del cono.
• Resultados corroborados por
SERNAGEOMIN
Fuente: Elaboración propia desde S1TBX.

Fuente: CEOS Pilot Proyect.
Santiago
Propuesta de Implementación
Integración de tecnología InSAR a las redes de monitoreo volcánico
SERNAGEOMIN:
Santiago
Conclusiones
 Tecnología Radar óptima para obtener datos terrenos en independencia de luz
solar y condiciones meteorológicas.
 Tecnología SAR permite resoluciones útiles desde espacio en el marco del estudio
de movimientos superficiales.
 Con InSAR es posible obtener DEM’s en ausencia de movimiento superficial.
 Asertividad de técnica DInSAR. Permite obtener diferenciales de movimiento a
nivel sub-métrico en superficies observadas, con asistencia de DEM externo.
 Eficiencia de Software S1TBX para fines interferométricos con Sentinel 1A.
 La erupción del Volcán Calbuco en Abril de 2015 provocó una subsidencia de
~12cm en dirección Línea de Visada (Line of Sight) en el flanco Sur-Oeste del cono.
 Datos corroborados por SERNAGEOMIN.
 Factibilidad de implementación Tecnología InSAR en redes de monitoreo volcánico
de SERNAGEOMIN.
Santiago
Conclusiones
• Tabla comparativa entre sistema InSAR, GPS e Inclinómetros:
InSAR GPS Inclinómetros
a) Instalaciones en terreno Innecesaria Necesaria Necesaria
b) Información en masa Posible Imposible Incompleta
c) Observación continua Imposible Posible Posible
d) Medición de deformación acumulada Indirectamente Observable Observable
e) Detección del movimiento 1 Dimensión 3 Dimensión 1 Dimensión
f) Tiempo de observación Cada 12/6 días Continua Continua
Fuente: Elaboración Propia.
Santiago
Gracias por su atención
Santiago

Presentación Tesis InSAR

Cargado por

Copyright:

Formatos disponibles

Presentación Tesis InSAR

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Presentación Tesis InSAR

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD TÉCNICA FEDERICO SANTA MARÍA

Academia de Ciencias Aeronáuticas

Interferometría SAR aplicada al estudio del movimiento en

Autor: David Hermosilla Diaz

Fuente: Google Earth

4.- Este pulso es captado en función de su

para formar una imagen de la superficie capturada.

Fuente: Carlos López

Resolución en azimut para sensor RAR:

Cada vez que la partícula recorre una distancia igual

La diferencia de fase entre dos pixeles vecinos de dos imágenes es igual a:

𝐵⊥ : Línea de Base perpendicular entre ambas capturas

Fuente: InSAR Imaging of Aleutian Volcanoes from Space.

Δ𝜙𝐷𝐼𝑛𝑆𝐴𝑅 = Δ𝜙𝑡𝑜𝑡𝑎𝑙 − 𝜙𝑡𝑜𝑝𝑜𝑔𝑟𝑎𝑓𝑖𝑎

Al realizar la diferencia entre la fase

FUENTE: Zhong Lu, Daniel Dzurisin. “InSAR Imaging of Aleutian Volcanoes”.

-La mayor actividad volcánica del mundo se ha

-Se debe al movimiento constante de las

0 Hawaiano No-Explosiva <100 mts. >1.000 m³ Diariamente 1983 Kilauea, HAW

3 Vulcaniano Violenta 3 – 15 Km >10.000.000 m³ Anualmente 1985 Nev del Ruiz, COL

Fuente: Elaboración propia desde S1TBX.

InSAR GPS Inclinómetros

a) Instalaciones en terreno Innecesaria Necesaria Necesaria

b) Información en masa Posible Imposible Incompleta

c) Observación continua Imposible Posible Posible

d) Medición de deformación acumulada Indirectamente Observable Observable

e) Detección del movimiento 1 Dimensión 3 Dimensión 1 Dimensión

f) Tiempo de observación Cada 12/6 días Continua Continua

Fuente: Elaboración Propia.

También podría gustarte