Tungurahua Cartesian Moment Tensors
Cartesian moment tensors from Tungurahua Volcano, Ecuador
data(egl)
A list of 84 moment tensors, each elelment consists of: lam1, lam2, lam3, vec1, vec2,vec3, ratio, force.
See below
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.
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)
}Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.