Quantifying soil and critical zone variability in a forested catchment through digital soil mapping

Abstract. Quantifying catchment-scale soil property variation yields insights into critical zone evolution and function. The objective of this study was to quantify and predict the spatial distribution of soil properties within a high-elevation forested catchment in southern Arizona, USA, using a combined set of digital soil mapping (DSM) and sampling design techniques to quantify catchment-scale soil spatial variability that would inform interpretation of soil-forming processes. The study focused on a 6 ha catchment on granitic parent materials under mixed-conifer forest, with a mean elevation of 2400 m a.s.l, mean annual temperature of 10 °C, and mean annual precipitation of ~ 85 cm yr−1. The sample design was developed using a unique combination of iterative principal component analysis (iPCA) of environmental covariates derived from remotely sensed imagery and topography, and a conditioned Latin hypercube sampling (cLHS) scheme. Samples were collected by genetic horizon from 24 soil profiles excavated to the depth of refusal and characterized for soil mineral assemblage, geochemical composition, and general soil physical and chemical properties. Soil properties were extrapolated across the entire catchment using a combination of least-squares linear regression between soil properties and selected environmental covariates, and spatial interpolation or regression residual using inverse distance weighting (IDW). Model results indicated that convergent portions of the landscape contained deeper soils, higher clay and carbon content, and greater Na mass loss relative to adjacent slopes and divergent ridgelines. The results of this study indicated that (i) the coupled application of iPCA and cLHS produced a sampling scheme that captured the greater part of catchment-scale soil variability; (ii) application of relatively simple regression models and IDW interpolation of residuals described well the variance in measured soil properties and predicted spatial correlation of soil properties to landscape structure; and (iii) at this scale of observation, 6 ha catchment, topographic covariates explained more variation in soil properties than vegetation covariates. The DSM techniques applied here provide a framework for interpreting catchment-scale variation in critical zone process and evolution. Future work will focus on coupling results from this coupled empirical–statistical approach to output from mechanistic, process-based numerical models of critical zone process and evolution.

[1]  N. Toomanian,et al.  Using Canonical Correspondence Analysis (CCA) to identify the most important DEM attributes for digi , 2011 .

[2]  W. Dietrich,et al.  Geomorphic transport laws for predicting landscape form and dynamics , 2013 .

[3]  Robert H. Whittaker,et al.  VEGETATION OF THE SANTA CATALINA MOUNTAINS, ARIZONA. V. BIOMASS, PRODUCTION, AND DIVERSITY ALONG THE ELEVATION GRADIENT' , 1975 .

[4]  Peter A. Troch,et al.  Separating physical and meteorological controls of variable transit times in zero‐order catchments , 2013 .

[5]  C. Lee Burras,et al.  Equations for Predicting Soil Organic Carbon Using Loss‐on‐Ignition for North Central U.S. Soils , 2002 .

[6]  Gary A. Peterson,et al.  Soil Attribute Prediction Using Terrain Analysis , 1993 .

[7]  G. Heuvelink,et al.  Optimization of sample patterns for universal kriging of environmental variables , 2007 .

[8]  Ramon C. Littell,et al.  Regression Using JMP , 2003 .

[9]  Craig Rasmussen,et al.  Geomorphically based predictive mapping of soil thickness in upland watersheds , 2009 .

[10]  K. Beven,et al.  A physically based, variable contributing area model of basin hydrology , 1979 .

[11]  S. Brantley,et al.  Approaches to Modeling Weathered Regolith , 2009 .

[12]  Budiman Minasny,et al.  A conditioned Latin hypercube method for sampling in the presence of ancillary information , 2006, Comput. Geosci..

[13]  P. Owens,et al.  Remembering Knowledge: An Expert Knowledge Based Approach to Digital Soil Mapping , 2013 .

[14]  Kate Maher,et al.  The dependence of chemical weathering rates on fluid residence time , 2009 .

[15]  T. G. Freeman,et al.  Calculating catchment area with divergent flow based on a regular grid , 1991 .

[16]  Igor V. Florinsky,et al.  Combined analysis of digital terrain models and remotely sensed data in landscape investigations , 1998 .

[17]  J. L. Ortiz,et al.  Comparação entre mapas de solos obtidos por sensoriamento remoto espectral e pelo método convencional , 2004 .

[18]  S. Miyamoto,et al.  Spectral properties of salt crusts formed on saline soils. , 2002, Journal of environmental quality.

[19]  Robert C. Reynolds,et al.  X-Ray Diffraction and the Identification and Analysis of Clay Minerals , 1989 .

[20]  N. McKenzie,et al.  Spatial prediction of soil properties using environmental correlation , 1999 .

[21]  D. Browning,et al.  Digital Soil Mapping in the Absence of Field Training Data: A Case Study Using Terrain Attributes and Semiautomated Soil Signature Derivation to Distinguish Ecological Potential , 2011 .

[22]  W. Dietrich,et al.  The soil production function and landscape equilibrium , 1997, Nature.

[23]  Cristiano Ballabio,et al.  A Plant ecology approach to digital soil mapping, improving the prediction of soil organic carbon content in alpine grasslands , 2012 .

[24]  W. Dietrich,et al.  Integration of geochemical mass balance with sediment transport to calculate rates of soil chemical weathering and transport on hillslopes , 2007 .

[25]  Dennis D. Eberl,et al.  User Guide to RockJock - A Program for Determining Quantitative Mineralogy from X-Ray Diffraction Data , 2003 .

[26]  Gerard B. M. Heuvelink,et al.  About regression-kriging: From equations to case studies , 2007, Comput. Geosci..

[27]  Feras M. Ziadat,et al.  Analyzing digital terrain attributes to predict soil attributes for a relatively large area , 2005 .

[28]  M. Jackson.  Soil Chemical Analysis - Advanced Course. , 1969 .

[29]  B. J,et al.  Soil regionalisation by means of terrain analysis and process parameterisation , 2002 .

[30]  M. Holleran.  Quantifying catchment scale soil variability in Marshall Gulch, Santa Catalina Mountains Critical Zone Observatory , 2013 .

[31]  A. Karnieli,et al.  Mapping of several soil properties using DAIS-7915 hyperspectral scanner data - a case study over clayey soils in Israel , 2002 .

[32]  O. Chadwick,et al.  From a black to a gray box ― a mass balance interpretation of pedogenesis , 1990 .

[33]  P. Troch,et al.  The effects of climate and landscape position on chemical denudation and mineral transformation in the Santa Catalina mountain critical zone observatory , 2011 .

[34]  J. Pelletier.  Quantitative Modeling of Earth Surface Processes , 2008 .

[35]  Pierre Goovaerts,et al.  Estimation or simulation of soil properties? An optimization problem with conflicting criteria , 2000 .

[36]  G. Milligan.  Ultrametric hierarchical clustering algorithms , 1979 .

[37]  John R. Jensen,et al.  Introductory Digital Image Processing: A Remote Sensing Perspective , 1986 .

[38]  J. Chorover,et al.  Fractionation of Dissolved Organic Matter by (Oxy)Hydroxide‐Coated Sands: Competitive Sorbate Displacement during Reactive Transport , 2014 .

[39]  D. Beaudette,et al.  Quantifying the Aspect Effect: An Application of Solar Radiation Modeling for Soil Survey , 2009 .

[40]  R. Jackson,et al.  Spectral response of a plant canopy with different soil backgrounds , 1985 .

[41]  E. Ben-Dor,et al.  The Spectral Reflectance Properties of Soil Structural Crusts in the 1.2‐ to 2.5‐μm Spectral Region , 2003 .

[42]  J. Harte,et al.  Soil Erosion: Data Say C Sink , 2008, Science.

[43]  William E. Dietrich,et al.  Constitutive mass balance relations between chemical composition, volume, density, porosity, and strain in metasomatic hydrochemical systems: Results on weathering and pedogenesis , 1987 .

[44]  Michael E. Schaepman,et al.  Representing major soil variability at regional scale by constrained Latin Hypercube Sampling of remote sensing data , 2013, Int. J. Appl. Earth Obs. Geoinformation.

[45]  Stephen J. Ventura,et al.  Fuzzy and isodata classification of landform elements from digital terrain data in Pleasant Valley, Wisconsin , 1997 .

[46]  Tomislav Hengl,et al.  Methods to interpolate soil categorical variables from profile observations: Lessons from Iran , 2007 .

[47]  J. Dozier,et al.  Estimating the spatial distribution of snow water equivalent in an alpine basin using binary regression tree models: the impact of digital elevation data and independent variable selection , 2005 .

[48]  E. Ben-Dor.  Quantitative remote sensing of soil properties , 2002 .

[49]  D. Sullivan,et al.  IKONOS Imagery to Estimate Surface Soil Property Variability in Two Alabama Physiographies , 2005 .

[50]  Paul E. Gessler,et al.  Modeling Soil–Landscape and Ecosystem Properties Using Terrain Attributes , 2000 .

[51]  Kristine Walraevens,et al.  Pumping test interpretation by combination of Latin hypercube parameter sampling and analytical models , 2009, Comput. Geosci..

[52]  P. Scull,et al.  A Comparison of Methods to Predict Soil Surface Texture in an Alluvial Basin , 2005 .

[53]  P. Fitze,et al.  Formulation of pedologic mass balance based on immobile elements: a revision. , 2000 .

[54]  C. Rasmussen,et al.  Lithologic controls on regolith weathering and mass flux in forested ecosystems of the southwestern USA , 2011 .

[55]  G. R. Mehuys,et al.  Landform classification from a digital elevation model and satellite imagery , 2008 .

[56]  John W. Salisbury,et al.  Infrared (8–14 μm) remote sensing of soil particle size , 1992 .

[57]  S. Mudd,et al.  Discrepancy between mineral residence time and soil age: Implications for the interpretation of chemical weathering rates , 2008 .

[58]  Gerard B. M. Heuvelink,et al.  Sampling design optimization for multivariate soil mapping , 2010 .

[59]  J. Demattê,et al.  Weathering sequence of soils developed from basalt as evaluated by laboratory (IRIS), airborne (AVIRIS) and orbital (TM) sensors , 2003 .

[60]  Gerard B. M. Heuvelink,et al.  Efficiency comparison of conventional and digital soil mapping for updating soil maps , 2012 .

[61]  Jon C. Helton,et al.  Use of Replicated Latin Hypercube Sampling to Estimate Sampling Variance in Uncertainty and Sensitivity Analysis Results for the Geologic Disposal of Radioactive Waste , 2010 .

[62]  Paul L. G. Vlek,et al.  Environmental correlation of three-dimensional soil spatial variability: a comparison of three adaptive techniques , 2002 .

[63]  F. Nachtergaele.  Soil taxonomy—a basic system of soil classification for making and interpreting soil surveys: Second edition, by Soil Survey Staff, 1999, USDA–NRCS, Agriculture Handbook number 436, Hardbound , 2001 .

[64]  S. Brantley,et al.  Chapter 1. CHEMICAL WEATHERING RATES OF SILICATE MINERALS: AN OVERVIEW , 1995 .

[65]  Lênio Soares Galvão,et al.  Relationships between the mineralogical and chemical composition of tropical soils and topography from hyperspectral remote sensing data , 2008 .

[66]  I. Moore,et al.  Digital terrain modelling: A review of hydrological, geomorphological, and biological applications , 1991 .

[67]  B. Minasny,et al.  On digital soil mapping , 2003 .

[68]  C. Rasmussen,et al.  Covariate selection with iterative principal component analysis for predicting physical soil properties , 2014 .

[69]  Luis A. Garcia,et al.  Comparison of Ordinary Kriging, Regression Kriging, and Cokriging Techniques to Estimate Soil Salinity Using LANDSAT Images , 2010 .

[70]  W. C. Krumbein.  : Factors of Soil Formation: A System of Quantitative Pedology , 1942 .

[71]  John Triantafilis,et al.  Digital soil mapping of compositional particle-size fractions using proximal and remotely sensed ancillary data , 2012 .

[72]  G. Heuvelink,et al.  A generic framework for spatial prediction of soil variables based on regression-kriging , 2004 .

[73]  Alex B. McBratney,et al.  Spatial prediction of soil properties from landform attributes derived from a digital elevation model , 1994 .

[74]  John P. Wilson,et al.  Terrain analysis : principles and applications , 2000 .

[75]  Alex B. McBratney,et al.  An overview of pedometric techniques for use in soil survey , 2000 .

[76]  T. Nauman.  Digital Soil-Landscape Classification for Soil Survey using ASTER Satellite and Digital Elevation Data in Organ Pipe Cactus National Monument, Arizona , 2009 .

[77]  M. Kleber,et al.  Erosion, deposition, and the persistence of soil organic matter: mechanistic considerations and problems with terminology , 2013 .