martes, 17 de septiembre de 2013

Generic Mapping Tools (GMT): Proyecciones

Proyección- lineal

Vamos a comenzar haciendo un mapa básico para los ejes x-y. Queremos que este vaya de 10 a 70 en el eje x, anotando cada 10 espacios,y de -3 a 8 en y, anotando cada 1. El mapa final debería ser de 4 X 3 pies de tamaño. Una forma de hacerlo es:

psbasemap -R10/70/-3/8 -JX4i/3i -B10/1:."My first plot": -P > plot.ps


Proyección- logarítmica

Vamos a construir un basemap para un gráfico log-log. Vamos a asumir que los datos raw en "x" van desde 3 a 9613 y que "y" va desde 3.2 10^20  a 6.8 10^24. Una manera de hacerlo es:

psbasemap -R1/10000/1e20/1e25 -JX9il/6il \-B2:"Wavelength (m)":/a1pf3:"Power (W)":WS > plot.ps



Proyección- Mercator

A pesar de los problemas en las altas latitudes, la proyección Mercator (-JM) sigue siendo la más usada por los científicos.Esta es una de más múltiples proyecciones cilíndricas ofrecidas por GMT; Nos enfocaremos en esta proyección. La sintaxis es simplemente

-JM(ancho)

Para hacer mapas de costas usamos pscoast el cual automáticamente nos da acceso a las costas, ríos y base de datos de fronteras derivadas de GSHHG. Además de los cambios comunes, vamos a necesitar algunas especificaciones adicionales para pscoast 



Uno de -W,-G,-S debe ser elegido. Nuestro primer ejemplo de Coastline es de Centro américa

pscoast -R-90/-70/0/20 -JM6i -P -B5g5 -Gchocolate > map.ps


Proyección-Eckert VI y Proyección VI


Proyección-Ortográfica







lunes, 16 de septiembre de 2013

Mecanismos focales

GEOMETRIA DE LA FALLA

Para describir la geometría de una falla, asumimos que esta es una superficie plana a través de la cual ha ocurrido movimiento relativo durante un terremoto.  Las observaciones geológicas de las fallas que alcanzan la superficie nos muestran que a menudo este es el caso, aunque las son complejidades son comunes. Así mismo veremos que este supuesto es usualmente (pero no siempre) consistente con los datos sísmicos. Por tanto la geometría de la falla es descrita en términos de la orientación del plano de falla y la dirección del deslizamiento a lo largo del plano.

La geometría de este modelo se muestra en la siguiente figura. El plano de falla se caracteriza por  n, este es el vector normal.  La dirección del movimiento está dado por d, el vector deslizamiento del plano de falla. El vector deslizamiento indica que la dirección en el cual el lado superior de la falla, conocido como hanging wall block  (Bloque de techo), se mueve con respecto al lado inferior, el foot wall block (Bloque de piso). Ya que el vector deslizamiento está en el plano de falla, este es perpendicular al vector normal.

Diferentes sistemas coordenados son útiles para el estudio de las fallas. Uno de ellos puede ser:
Alineado de forma que el eje x1 está en la dirección del strike, la intersección del plano de falla con la superficie terrestre. El eje x3 apunta hacia arriba y el eje x2 es perpendicular a los otros dos. El ángulo dip (delta) entrega la orientación del plano de falla respecto a la superficie. Como el eje x1 puede estar definido en dos direcciones, 180º , este se escoge de modo que el dip es medido en el eje  -x2  es menor que 90º.
La dirección del movimiento es representada por por el ángulo dip, (lambda) medido en sentido antihorario en el plano de falla desde la dirección x1, la que nos da el movimiento del bloque de techo respecto al bloque de piso. Para orientar de este sistema relativo a uno es geográfico, la fault strike (phi) es definido como el ángulo en el que el plano de la superficie terrestre es medido en sentido horario desde el norte al eje x1.

Alternativamente, la orientación de la falla y el slip pueden ser descritos dando vectores normales y de slip en un sistema de coordenadas geográficas con x' apuntando al norte, y' apuntando al oeste y z' hacia arriba. En ese  sistema coordenada, el vector unitario normal al plano de falla es:


y el vector slip, un vector unitario a la dirección del slip es,


Estos dos diferentes sistemas coordenados (phi_f, delta, lambda) y (n,d) son útiles para diferentes propósitos. Algunos cálculos son fáciles hacerlos con respecto a la falla, mientras otros son más fáciles hacerlos con respecto a las direcciones geográficas.

Aunque la dirección del slip varía de modo que el ángulo de slip tiene un rango entre 0º y 360º, hay diferentes geometrías básicas de fallas, descritas por valores especiales del angulo slip, son útiles y buenas tenerlas en mente. 


 Cuando los dos lados de la falla se deslizan horizontalmente uno con el otro, ocurre el movimiento pure strike-slip.

- Cuando lambda = 0º , el bloque de techo se mueve a la derecha y el movimiento se llama left-lateral 

- Similarmente, cuando lambda = 180º ocurre el movimiento right-lateral 


Para decir cual es cual, miramos a través de la falla y vemos la manera en que el otro lado se mueve.
Las otras geometrias básicas describen el movimiento dip-slip 

- Cuando lambda = 270º, el bloque de techo se desliza hacia abajo, causando una falla normal
Cuando lambda = 90º ocurre el caso opuesto, el bloque de techo se desliza hacia arriba originando una falla inversa o thrust faulting.


Muchos terremotos consisten en una combinación de estos movimientos y tienen ángulos dip entre estos valores. Es por ello útil, cuando pensamos en los mecanismos de los terremotos recordar estos 3 movimientos básicos. Los tipos básicos de falla están relacionados con las orientaciones de los esfuerzos principales. 

Las fallas reales, obviamente, tienen dimensiones finitas y geometrias complicadas. Si tratamos a las fallas como rectángulos, la dimensión a través del strike es llamado largo de la falla, y la dimensión en la dirección del dip es conocido como ancho de la falla. Las geometrías de los terremotos pueden ser mucho más complicadas que un rectángulo. La falla puede curvarse y requerir una descripción en tres dimensiones. La ruptura puede ocurrir sobre un gran periodo de tiempo y consistir en una serie de sub-eventos en diferentes partes de la falla con diferentes orientaciones.  Estos eventos complejos, sin embargo pueden ser tratados como una superposición de eventos simples. 

Si entendemos que las ondas sísmicas  generadas por una simple falla en dos dimensiones, podemos modelar este resultado a partir de un set más complicado de rupturas. Esta aplicación del principio de superposición se basa en el supuesto de la elasticidad lineal y es análoga a la manera en que construimos ondas sísmicas  sumando modos normales. 


sábado, 14 de septiembre de 2013

Construyendo paletas con GMT

GMT tiene 20 tablas de colores (master cpt files). El siguiente es un esquema de cada una de ellas:


Los comandos makecpt y grd2cpt son usados para acceder a estas tablas maestras y re-escalarlas para ajustarse al rango de valores z que el usuario desea. La tabla cpt final puede ser discreta (como la que está encima) o bien continua (parte baja).

El capítulo 4 del tutorial oficial de GMT  nos enseña  a diseñar nuestras propias paletas (ya sea discretas a continuas), con ejemplo bastante sencillos (que también voy a mostrar en esta entrada, pero con diferentes archivos claro, para ser un poquito original :P ) y posteriormente les voy a presentar un pequeño tutorial (bastante básico debo admitir) de cómo hacerla a partir de un archivo topográfico (.grd).

Primero, vamos a acceder a cpt-city, que es una página con un gran lista de paletas que puedes descargar según el tema que desees graficar: sismos, topografías, batimetría, clima, etc etc.

Para hacer archivos cpt discretos y continuos para datos que van de -20 a 60, con cambios de color cada 10, intenta con estas dos variantes (nota: la paleta reds_01.cpt yo la bajé de la página antes mencionada y la guardé en GMT4.5.8/share/cpt, elige una paleta cualquiera y escribes su nombre luego del -C  y  ve que pasa), el -Z es para que nuestra paleta sea continua:


makecpt -Creds_01 -T-20/60/10 > disc.cpt

makecpt -Creds_01 -T-20/60/10 -Z > cont.cpt




Opción
Propósito


-CcptfileEl archivo cpt requerido


-Dxpos/ypos/length/width[h]Fija la posición del centro/izquierda y dimensiones de la escala


Añadir para obtener una barra horizontal to get horizontal bar. 


-Imax_intensityAgregar efectos de iluminación




grdraster 9 -R-108/-103/35/40 -Gus.nc

Este archivo está ya viene en el GMT (lo buscan por us.nc)
Usando grdinfo encontramos que los datos van de ~1000m a los ~4300m, de modo que hacemos una paleta cpt acorde:

makecpt -Creds_01 -T1000/5000/500 -Z > topo.cpt




Opción
Propósito


-EdpiFija la resolución deseada de la imagen
[Default is data resolution]


-IintenfileUsa una iluminación artificla  usando las intensidad de un archivo intensfile


-MFuerza de las sombras grises usando la conversión (television) YIQ 



grdimage us.nc -JM6i -P -B2 -Ctopo.cpt -V -K > topo.ps  




Opción
Propósito


-AazimuthDirección azimutal para los gradientes


-MIndica que esta es una grilla geográfica


-N[t|e][norm[/offset]]Normaliza los gradientes por norm/offset [= 1/0 by default].
Insertar  para normalizar por la transformación  tan-1.
Insertar  para normalizar por la distrubución de Laplace acumulativa.

psscale -D3i/8.5i/5i/0.25ih -Ctopo.cpt -I0.4 -B/:m: -O >> topo.ps



Ahora vamos a plotear nuestras nuevas tablas de colores con psscale;  veamos las siguientes características listadas en la tabla:
La opción, -B puede ser usada para fijar el titulo y una etiqueta de unidad.

psbasemap -R0/8.5/0/11 -Jx1i -P -B0 -K > bar.ps  
psscale -D3i/3i/4i/0.5ih -Cdisc.cpt -B:discrete: -O -K >> bar.ps
psscale -D3i/5i/4i/0.5ih -Ccont.cpt -B:continuous: -O -K >> bar.ps
psscale -D3i/7i/4i/0.5ih -Cdisc.cpt -B:discrete: -I0.5 -O -K >> bar.ps
psscale -D3i/9i/4i/0.5ih -Ccont.cpt -B:continuous: -I0.5 -O >> bar.ps




Una vez hecho un archivo cpt es casi directo generar una imagen a color a partir de una grilla. Aquí se va a extraer un subset del   global 30" DEM (data id 9) del USGS:

grdraster 9 -R-108/-103/35/40 -Gus.nc

La grilla us.nc puedes encontrarla en el directorio de instalación de GMT . Usando grdinfo vemos que los datos rondan entre los ~1000m a ~4300m de modo que hacemos una nueva paleta acorde a estos rangos:

makecpt -Creds_01 -T1000/5000/500 -Z > topo.cpt


 Las imagenes a color se hacen con  grdimage  el cual toma los opciones de los comandos usuales  (por defecto el -R  se toma del archivo de datos ) y un cptfile; otras opciones son:


Opción
Propósito


-EdpiFija la resolución deseada de la imagen


-IintenfileUsa una iluminación artificial a partir de un archivo
 
intensfile


-MFuerza de las sombras grises usando la conversión
YIQ 


El mapa resultante es la región de Rocky Mountain. Lo que necesita es una iluminación artificial. Queremos simular un sombreado por una fuente solar en el este, por lo tanto derivamos las intensidades requeridas a través de gradientes topográficos en dirección N90ºE usando grdgradient . Otras opciones disponibles para los archivos de salida son: 



Opción
Propósito


-AazimuthDirección azimutal para los gradientes


-MIndica que es una grilla geográfica


-N[t|e][norm[/offset]]Normaliza los gradientes por  norm/offset [= 1/0 por default].
Insertar  para normalizar por la transformación tan-1 .
Insertar para normalizar por la distribución acumulativa de  Laplace.



Tanto -Ne como -Nt otorgan gradientes bien comportados.  Personalmente,  preferimos usar la opción -Ne; el valor de la norma  norm  es subjetivo y puedes experimentar entre un rango de 0.5 - 5. Para nuestro caso escogimos:

grdgradient us.nc -Ne0.8 -A100 -M -Gus_i.nc

Dado un archivo  cpt  y dos grillas podemos crear una imagen con relieve sombreado:

grdimage us.nc -Ius_i.nc -JM6i -P -B2 -Ctopo.cpt -K > topo.ps  


jueves, 12 de septiembre de 2013

Primeros mapas con GMT (Generic Mapping Tools)

Sin lugar a duda,  para las personas que nos dedicamos a las ciencias de la Tierra tener a disposición  un software para crear nuestros mapas con total libertad modificar, añadir datos, fácil de usar etc etc y sobretodo que sea gratuito sea hace imprescindible. En la Universidad durante del primer semestre dictan un curso introductorio a  Linux , donde precisamente nos enseñan a ocupar el Generic Mapping Tools más conocido por sus siglas GMT.

 Obviamente no los voy a aburrir con "la historia tras GMT" pero sí voy a dar algunos tipos que siempre son útiles al momento de querer llevar de manera gráfica nuestro trabajo.

La instalación de GMT pueden seguirla paso a paso a través del manual que deja Matt cada año:
http://mttmllr.com/GMT/contenidos/node49.html

Los requerimientos son:

  • - Tener un ordenador con sistema operativo Linux.
Ejemplo mapa hecho en GMT: Réplicas terremoto del Maule de 2010, Mw 8.8