21.1. Zona Wadati-Benioff I. Perfil Topografico + Sismos/Mecanismos Focales.

Graficar sismos y mecanismos focales en un perfil E-W.

 

 

Introducción:En las zonas de convergencia activa, se produce la subducción de una placa tectónica por debajo de otra. En la placa que desciende se producen gran cantidad de sismos (y mecanismos focales) que son detectados en la superficie. Al graficarlos en perfiles perpendiculares a la zona de subducción, estos se alinean y conforman una zona de sismicidad denominada zona de Wadati-Benioff. A continuación se muestra un ejemplo de Sudamérica a los 32 ºS.

0. Datos:

• DEM: Grilla GMRT 3.6 (ver anexo 4).
• Sismos USGS: Archivos query_*.csv (USGS.zip; ver ANSS en base de datos).
• Mecanismos Focales: Archivos CMT_*.txt (comprimidos en CMT_1976-2019.zip; ver Global CMT en base de datos).

Script Mapa
Script Perfil

0. Mapa Base: con los sismos y mecanismos focales ploteados y perfil dibujado (ver 4.1 y 17.1 para detalles de los comandos).

1. Mapa base 1: Se filtran que se localizan en la región que se extiende 1 grado hacia el norte y el sur del perfil (entre -33 y -31 ) que van a ser graficados en el perfil perpendicular A-A'. Se utiliza select con el argumento -R donde se indica la región de interés.

SET    REGION2=-74/-64/-33/-31
gmt select -R%REGION2% "query_*.csv" -h1 -i2,1,3 > "temp_sismos"
gmt select -R%REGION2% "CMT_*.txt" > "temp_CMT"

gmt psxy -R%REGION% -J -T -O -K >> %OUT%
gmt psxy   -R -J -O -K "temp_sismos" >> %OUT% -Sc0.05   -C"temp_seis"
gmt psmeca -R -J -O -K "temp_CMT"    >> %OUT% -Sd0.12/0 -Gorange 

2. Perfil topográfico (superior): Se utilizan los comandos explicados en 11.1. Perfil Topográfico para crear el perfil, extraer los datos del dem, calcular cuanto mide (en km), graficar el perfil con la exageración vertical y su orientación.

3. Perfil en profundidad (inferior): Se adjunta un perfil en profundidad hasta los 300 km de profundidad.

SET    DepthMin=0
SET    DepthMax=300

gmt psbasemap -R -J -O -K -Bxaf+l"Longitud (\260)"  >> %OUT% --MAP_FRAME_AXES=Sn
gmt psbasemap -R -J -O -K -Byaf+l"Profundidad (km)" >> %OUT% --FONT_ANNOT_PRIMARY=8,Helvetica,blue --MAP_FRAME_AXES=wE

4. Filtrar datos: Con los mismos comandos del punto 2, se filtran los datos de los sismos y de los mecanismos focales según la región geográfica entre los 31 y 33 latitud sur.

SET    REGION2=%LongMin%/%LongMax%/-33/-31
gmt select -R%REGION2% "query_*.csv"  -h1 -i2,1,3 > "temp_sismos"
gmt select -R%REGION2% "CMT_*.txt" > "temp_CMT"

5. Plotear sismos en profundidad versus longitud: con psxy se grafican los sismos previamente filtrados (temp_sismos). Es importante indicar las columnas con los datos de longitud y profundidad (-i0,2). Se utiliza una CPT con valores entre 0 y 300 para pintar los sismos según su profundidad.

gmt makecpt > "temp_seis.cpt" -Crainbow -T0/300 -Z -I 
gmt psxy  -R -J -O -K "temp_sismos" >> %OUT% -W.1 -C"temp_seis" -Sc0.05c -i0,2,2
 

5. Mecanismos focales: pscoupe permite graficar los mecanimos focales previamente filtrados (temp_CMT) en un perfil con los siguientes argumentos:

• -Aa: Perfil (longitud y latitud inicial y final), la inclinación del perfil (90 = vertical), extensión hacia los lados del perfil (10 mil km, para incluir a todos los mecanismos) y profundidad mínima y máxima.
• -Sd0.15/0: beachballs de un tamaño de 0,15 (para los eventos de magnitud 5 en pulgadas por default). El segundo valor corresponde al tamaño de la fuente usado para el nombre del evento. Con 0 no se escribe.
• -Q: No crea los archivos de información de los mecanismos focales.
• -G: Color para los cuadrantes compresivos.
• -M: Tamaño fijo.

gmt pscoupe -R -J -O -K "temp_CMT"    >> %OUT% -Sd0.15/0 -Gred -M -Aa%LongMin%/%LatMin%/%LongMax%/%LatMax%/90/10000/%DepthMin%/%DepthMax%f -Q

6. Otros Perfiles Longitudinales:

Perfil a los 28 S entre 74 y 62 L y hasta 700 km de profundidad.

Perfil a los 28 S entre 74 y 42 L y hasta 300 km de profundidad.

Perfil a los 28 S entre 74 y 42 L y hasta 300 km de profundidad.

Perfil a los 28 S entre 74 y 42 L y hasta 300 km de profundidad.