R: Tungurahua Cartesian Moment Tensors

egl {RFOC} | R Documentation

## Tungurahua Cartesian Moment Tensors

### Description

Cartesian moment tensors from Tungurahua Volcano, Ecuador

### Usage

    data(egl)

### Format

A list of 84 moment tensors, each elelment consists of: lam1, lam2, lam3, vec1, vec2,vec3, ratio, force.

### Source

See below

### References

Kim, K., Lees, J.M. and Ruiz, M., (2014) Source mechanism of Vulcanian eruption at Tungurahua Volcano, Ecuador, derived from seismic moment tensor inversions, J. Geophys. Res., February, 2014. Vol. 119(2): pp. 1145-1164.

### Examples

    data(egl)

    typl1=c(2,4,7,12,13,16,17,18,19,20,24,25,26,27,
    28,29,30,31,33,34,35,36,37,38,40,43,50,
    59,62,73,74,77,8,79,80,81,83,84)
    typl2=c(5,6,8,9,10,11,14,15,22,42,46,47,48,49,
    51,52,53,54,55,56,57,58,60,61,63,72,82)

    evtns=1:84

    par(mfrow=c(1,2))
    T1 = TapeBase()
    TapePlot(T1)

    for(i in 1:length(egl))
    {
    i1 = egl[[i]]

    E1 = list(values=c(i1$lam1, i1$lam2, i1$lam3),
    vectors = cbind(i1$vec1, i1$vec2, i1$vec3))

    testrightHAND(E1$vectors)

    E1$vectors = forcerighthand(E1$vectors)

    mo=sort(E1$values,decreasing=TRUE)
    # M=sum(mo)/3
    # Md=mo-M
    h = SourceType(mo)
    h$dip = 90-h$phi

    h1 = HAMMERprojXY(h$dip*pi/180, h$lam*pi/180)

    if(i %in% typl1) { col="red" }else{col="blue" }
    points(h1$x, h1$y,  pch=21, bg=col )

    }

    par(mai=c(0,0,0,0))
    hudson.net()
    for(i in 1:length(typl1))
    {
    egv=egl[[typl1[i]]]
    m=c(egv$lam1,egv$lam2,egv$lam3)
    col='red'
    hudson.plot(m=m,col=col)
    }

    for(i in 1:length(typl2))
    {
    egv=egl[[typl2[i]]]
    m=c(egv$lam1,egv$lam2,egv$lam3)
    col='blue'
    hudson.plot(m=m,col=col,lwd=2)
    }

[Package RFOC version 3.4-3 Index]