Evaluation of EEMT as calculated from MAT and MAP
Submitted by dylan on Fri, 2007-11-23 03:37.
EEMT Calculation from MAP and MAAT [See Rasmussen & Tabor; SSSAJ 71: 1719 - 1729]
## EEMT function based on MAP and MAT
f <- function(mat, map) { e <- 347134 * exp( -0.5*( ( (mat - 21.5)/-10.1 )^2 + ( (map - 4412)/1704 )^2 ) ) }
##
## generate sequence of possible mat and map values
mat.seq <- seq(-10, 30, by=2)
map.seq <- seq(100, 3500, by=50)
## create DF of all possible combinations
d <- expand.grid(mat=mat.seq, map=map.seq, KEEP.OUT.ATTRS=FALSE\)
## predict EEMT from mat and map
d$eemt <- f(d$mat, d$map)
## 3D viz: wireframe() works better than persp()
library(lattice)
wireframe(eemt ~ map * mat, data=d)
eemt_wireframe: EEMT vs. MAP vs. MAT
Great Group vs. EEMT: US Soil Taxonomy Great Group level classification vs. EEMT
EEMT for Central CA: EEMT as calculated from PRISM data. Points are locations of Soil-Veg records.
EEMT Critical Values (?): An approach for identifying the 'red soil zone' ???
|