Difference between revisions of "STATSGO2-soil data"
m |
|||
Line 29: | Line 29: | ||
[http://soildataviewer.nrcs.usda.gov/userguide.aspx] | [http://soildataviewer.nrcs.usda.gov/userguide.aspx] | ||
− | == Step 5: | + | == Step 5: An Example of estimating soil properties using SDV in the Silver Fork of the American == |
Definitions: | Definitions: | ||
Revision as of 00:44, 4 February 2010
Contents
- 1 Step 1: Get soil data
- 2 Step 2: Install Soil Data Viewer
- 3 Step 3: Prepare Access database
- 4 Step 4: Working with Soil Data Viewer in ArcGIS
- 5 Step 5: An Example of estimating soil properties using SDV in the Silver Fork of the American
- 5.1 1) Soil Physical Properties: Surface Texture
- 5.2 2) Soil Qualities and Features: Depth to Any Restrictive Layer (RHESSys requires soil depth)
- 5.3 3) Soil Physical Properties: Saturated Hydraulic Conductivity (RHESSys requires Ksat at surface and decay with depth)
- 5.4 4) Soil Physical Properties: %Clay, %Sand, %Silt (used in RHESSys only for carbon cycling)
Step 1: Get soil data
Download state data from here: [1]
Download county or survey level data from here: [2]
Both sources should include a) tabular data b) spatial data c) Microsoft Access template
State data is spatially coarser than county or survey level data
Step 2: Install Soil Data Viewer
You need 1) Install ArcGIS 9.2 or later for soil data viewer 5.2 on a Windows XP system 2) Install soil data viewer version 5.2 [3]
Step 3: Prepare Access database
1) Open Access template 2) Enable macros 3) Type in full path name to tabular data 4) Import should run automatically
Step 4: Working with Soil Data Viewer in ArcGIS
1) Open ArcGIS - View > Toolbars > Soil Data Viewer Tools 2) Click on the Soil Data Viewer button 3) Browse to Access database 4) Choose the soil data layer to connect to the database
Soil data viewer online user guide is here: [4]
Step 5: An Example of estimating soil properties using SDV in the Silver Fork of the American
Definitions:
a) Dominant Condition (DCD)
b) Dominant Component (DCP)
c) Weighted Average (WA)
d) All Layers (ACP)
1) Soil Physical Properties: Surface Texture
Useful for a) naming the soil .def file b) understanding depth to soil restrictive layer c) understanding ksat at surface
2) Soil Qualities and Features: Depth to Any Restrictive Layer (RHESSys requires soil depth)
a) Generate Depth DCD, DCP, WA, ACP using both < and > for tie break rule with no cutoff (8 scenarios)
b) Compare these to summer Landsat imagery and use expert judgment to determine the most accurate layer
c) Depth To Any Restrictive Layer results for Silver Fork
- ACP, > 127,143,153 - not much spatial variability, bedrock not represented
- ACP, < 0 - basically junk
- WA - tie break doesn't matter, 17 for texture un-weathered bedrock
- DCP - tie break doesn't matter, very shallow in upper basin
- DCD, > zeros where surface texture is bedrock, 77,127,201
- DCD, < zeros where surface texture is bedrock, 77,127,201 but shallower than DCD, > (obviously)
3) Soil Physical Properties: Saturated Hydraulic Conductivity (RHESSys requires Ksat at surface and decay with depth)
a) Generate Depth DCD, DCP, WA, ACP for ALL LAYERS using both slowest and fastest for tie break rule with no cutoff (8 scenarios)
b) All Layers Ksat results for American Silver Fork
- DCD - tie break doesn't matter, portions with surface texture bedrock have no data. Range is 91 um/s, 29 um/s 18 um/s
- DCP - tie break doesn't matter, large portions of upper basin have no data. Range is 91 um/s to 28 um/s for lower
- WA - tie break doesn't matter, Ksat is 60 um/s in upper basin and 26 um/s in lower basin. More spatial variability than ACP >
- ACP < - Ksat is too slow
- ACP > - Ksat is 91 um/s in upper basin and 28 um/s in lower basin. Less spatial variability than WA
c) 50 cm layer increment - WA results for American Silver Fork (note 50cm increment chosen given depth to bedrock in step 1 - coarser might be appropriate for lower basin). Weighted Averaged used - this means lower/upper is the same so only 5 cases (0-50cm,51-100cm,101-150cm,151-200cm, 201-250cm).
- 0-50cm 64 um/s to 25 um/s
- 51-100cm: 64 um/s to 25 um/s
- 101-150cm: 69.4 um/s, 44.9 um/s, 27.5 um/s, 25.7 um/s, 16.0 um/s, 9.2 um/s -
- 151-200cm: 54.8 um/s, 32.8 um/s, 20.7 um/s, 14.5 um/s, 9.2 um/s -
- 201-250cm: 9.2 um/s - this is for the lowest elevation unit in the basin which is likely the deepest
d) 25 cm layer increment - MA results for American Silver Fork (note 50cm increment chosen given depth to bedrock in step 1 - coarser might be appropriate for lower basin).
- export attributes table, open with excel and save 2 columns (# and ksat) as csv
- csvread in matlab and plot Ksat for each soil unit (21 for Silver Fork, 3 are lakes) plot Ksat w/ depth
- determine m (how?) from Ksat in MATLAB
4) Soil Physical Properties: %Clay, %Sand, %Silt (used in RHESSys only for carbon cycling)
a) Focus is not on carbon cycling - simplify and use WA, lower and upper were the same for all 3 variable Clay - range is 2.7% to 12.7% Sand - range is 66.3% to 80.4% Silt - range is 15.9% to 24.1%