GSIF (version 0.5-3)

edgeroi: The Edgeroi Data Set

Description

Soil samples and covariate layers for the Edgeroi area in NSW, Australia (ca 1500 square-km).

Usage

data(edgeroi)

Format

The edgeroi data set contains two data frames --- sites and horizons. Sites table contains the following columns:

Horizons table contains the following columns:

The edgeroi.grids data frame contains a list of covariates at 250 m resolution:

The edgeroi.grids100 data frame contains a list of covariates at 100 m resolution prepared for the study area:

Details

The Edgeroi is one of the standard soil data sets used to test soil mapping methods in Australia. Out of 359 profiles, 210 sites were sampled on a systematic, equilateral triangular grid with a spacing of 2.8 km between sites, the other sites are distributed more irregularly or on transects. The data set is described in detail in Malone et al. (2010) and McGarry et al. (1989). The edgeroi contains only a subset of the original NatSoil records. Observed soil classes for TAXGAUC are (alphabetically): Alluvial soil ("A"), Brown clay ("BC"), Black earth ("BE"), Earthy sand ("ES"), Grey clay ("GC"), Grey earth ("GE"), No suitable group ("NSG"), Prairie soil ("PS"), Rendzina ("R"), Red-brown earth ("RBE"), Red clay ("RC"), Red earth ("RE"), Red podzolic soil ("RP"), Solodic soil ("SC"), Soloth ("SH"), Solonchak ("SK"), Siliceous sand ("SS"), and Solonetz ("SZ").

References

• Malone, B.P., McBratney, A.B., Minasny, B. (2010) Mapping continuous depth functions of soil carbon storage and available water capacity. Geoderma 154, 138-152.
• McGarry, D., Ward, W.T., McBratney, A.B. (1989) Soil Studies in the Lower Namoi Valley: Methods and Data. The Edgeroi Data Set. (2 vols) (CSIRO Division of Soils: Adelaide).
• Minty, B., Franklin, R., Milligan, P., Richardson, L.M., and Wilford, J., (2009) The Radiometric Map of Australia. Exploration Geophysics, 40(4), 325-333.

Examples

library(rgdal)
library(aqp)
library(sp)

data(edgeroi)
edgeroi$sites[edgeroi$sites$SOURCEID=="399_EDGEROI_ed095_1",]
edgeroi$horizons[edgeroi$horizons$SOURCEID=="399_EDGEROI_ed095_1",]
## spPoints:
sites <- edgeroi$sites
coordinates(sites) <- ~ LONGDA94 + LATGDA94
proj4string(sites) <- CRS("+proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs")
sites <- spTransform(sites, CRS("+init=epsg:28355"))

## Not run: 
# ## plot points and grids:
# pnts <- list("sp.points", sites, pch="+", col="black")
# ## load the 250 m grids:
# con <- url("http://gsif.isric.org/lib/exe/fetch.php?media=edgeroi.grids.rda")
# load(con)
# str(edgeroi.grids)
# gridded(edgeroi.grids) <- ~x+y
# proj4string(edgeroi.grids) <- CRS("+init=epsg:28355")
# spplot(edgeroi.grids[1], sp.layout=pnts)
# ## load the 100 m grids:
# con2 <- url("http://gsif.isric.org/lib/exe/fetch.php?media=edgeroi.grids100.rda")
# load(con2) 
# str(edgeroi.grids100)
# gridded(edgeroi.grids100) <- ~x+y
# proj4string(edgeroi.grids100) <- CRS("+init=epsg:28355")
# spplot(edgeroi.grids100["TI1LAN6"], sp.layout=pnts)
# ## End(Not run)