Prepared by David R. Maidment
This exercise is intended to introduce you to the use of soils data from the STATSGO soils coverage. It is divided into three main parts:
This exercise uses workstation Arc/Info version 7.0 and also Arcview2 or Arcview3. It uses a watersheds coverage Covshd10 which was derived from 3 arc-second digital elevation data for the West Austin map sheet using the procedure described in Exercise 2. It also uses the Statsgo coverage for Texas, and a clipped version of this called Soilshed which has been clipped to conform to the watershed boundaries of Covshd10, and a similarly clipped version of the Rf1 rivers coverage called Rf1shed, so that you can see the association of rivers and soils in the study area. Tabular data on the soil properties are provided by the Info tables Mapunit, Comp, and Layer, which are held separately from the soil coverages in the Info directory. To save space in your file area, you will just view the Statsgo coverage for Texas from my directory using Arcview, rather than copying it over into your own file directory.
Make a new directory for this assignment, say ex5 and change directories to this directory.
Download the data via anonymous ftp from ftp.crwr.utexas.edu/pub/gisclass/ex5. If necessary, see instructions on how to use anonymous ftp. You need the following coverages:
These files are in export format so you will need to import them using the arc import command. For the coverages, use:
Arc: import cover covshd10 covshd10
Arc: import cover rf1shed rf1shed
Arc: import cover soilshed soilshed
and for the info files use:
Arc: import info comp comp
Arc: import info layer layer
Arc: import info mapunit mapunit
Once you've imported the coverages and info files, you can use the Unix command rm to remove the .e00 files that you used to acquire the data. You can remove several files at once by listing them one after another with a space in between. Don't put a comma after the file names or the rm command interprets the comma as being a part of the file name:
$ rm covshd10.e00 rf1shed.e00 soilshed.e00 comp.e00 layer.e00 mapunit.e00
Start up Arcview (you need to be working on the alpha system to access the Statsgo. Open a new View and add the theme Statsgo which is resident in the directory /home1/alpha62/maidment. If you have trouble accessing this file, I've put an export file statsgo.e00 on the anonymous ftp site with the other files for this exercise, but beware, as this file is 17MB in size! Add the theme Covshd10 which is in your working directory. You can see the large number of Statsgo polygons which are needed to describe the soils of Texas (actually there are 4031 polygons in the Statsgo coverage for Texas). By superimposing the watersheds coverage on top of the soils you can see the region of analysis in Central Texas that we will work with for the remainder of the exercise.
The properties of the Statsgo coverage for Texas are:
Clearly this is an immense data base, but the large size of its polygons makes it suitable for analysis of fairly large areas which contain several Statsgo polygons. There is another database called SSURGO, still under development by US Dept of Agriculture, which details soil polygons at the component level rather than at the map unit level.
In order to create a study region that has reasonably sized files, the Statsgo coverage of Texas has been clipped using a set of watersheds Covshd10 derived for the West Austin region of Central Texas in Exercise 2.
Use Edit/Delete Themes to delete the Statsgo theme. Add themes from your working directory for Soilshed, Rf1shed. Highlight Soilshed and using the legend editor, use the field Muid> to label the mapunits. Muid is the mapunit identification number for each unique combination soil components depicted in Statsgo.
Notice how the properties of the soils are closely associated with the location of rivers, particularly in the North-East part of the Central Texas area that you now have displayed. Use the tool in the View window to determine the map unit that underlies many of the rivers in that area (Mapunit Muid = Tx560). Lets take a look at the properties of this map unit.
Click on Tables in the Project window, select Add, go to the work area where you have copied the files at the beginning of this exercise and click on the Info directory. In the left bottom box instead of "List Files of Type " .dbf (dbase files), select INFO files. From the Info directory, add the mapunit and comp tables to the Project. Each of these tables has a lot of fields that are not needed in this exercise so lets delete these unwanted fields from the display. Highlight the Mapunit table, and use Table/Properties to click off the check boxes Visible in all Fields except for Muid and Muname. Drag the end of the Muname field title bar to the right so you can read all of the Muname. You'll see that it is a combination of several different soil names.
Similarly, for the Comp table, click off all fields except for Muid, Seqnum, Compname, Comppct, Slopel, Slopeh, Surftex, and Hydgrp which is located near the end of the table.
These fields represent:
Linking the Tables:
Now we are going to do something cool. The Comp table is going to be called the Source table and the Mapunit table the Destination table, and we are going to establish a linkage between these two tables: Click on the "Source" table (Comp) and select the field to be used for the linkage (Muid) by clicking on its field name (you will see it gets depressed in the table). Click on the "Destination" table (Mapunit), and select the key field to make the linkage (again Muid). Go to the Tables pulldown menu and near the bottom you will see Link you should see the blue bar at the bottom of the Arcview window, rev up indicating that some action is taking place. Now, go to the Mapunit table and select mapunit TX560 by clicking on it (it goes yellow). Bring the selected feature up to the top of table display using the Uparrow button on the top line of the tool bar. If necessary, click on the Comp table and likewise hit the button and you should see the component records for the selected mapunit highlight at the top of the table. Pretty neat!! When you select a new mapunit, the new components are chosen and highlighted at the top of the component table.
This is a one-to-many relate where one record in the mapunit table is related to many records in the component table.
The average properties of soils in a particular map unit can be determined by using a weighted average of the component properties where the weights are equal to the percentages contained in the comppct field. For example, map unit TX560 contains soils with three surface textures: C (Clay) for component 1, CL (Clay Loam) for components 2 and 3, and FSL (Fine Sandy Loam) for components 4 and 5. The percentage of the total area in the map unit occupied by these components is 69%, 23%, 3%, 3%, and 2%, respectively. These values can be summed to give 69% surface texture Clay, 26% surface texture Clay Loam, and 5% surface texture Fine Sandy Loam.
To be turned in: Describe the TX560 soil in terms of its properties in the Mapunit and Component Tables. How many components does it have? What are their names? What percentage of the map unit does each component comprise? What is the predominant surface slope where this soil unit is found? What is the dominant soil texture? What percentage of the soil is in hydrologic soil groups A, B, C, D? Do these soil properties make sense considering where this soil is usually located?
The component properties describe the soil as a whole. Soils are usually divided into soil horizons beginning at the surface, and each horizon or layer has different physical properties. These are contained in the Layer table.
In the Project window, click on Tables and Add the Layer table to the display in the same manner that you added the mapunit and component tables. Click off all the layer properties except the following: Muid, Seqnum, Layernum, Laydepl, Laydeph, Awcl and Awch (Awc is located near the end of the table). These properties represent:
Now lets link the Layer and Comp tables together. Use the same procedure
as before: click on the Layer table and highlight the
Now lets see what we have got here. The mapunit table is linked to the comp table and the comp table to the layer table. Click on Muid TX560 in the mapunit table and once it is highlighted, you should see the corresponding components and layers highlighted. Use the Uparrow icon to put the selected records at the top of the table if necessary. Pretty cool!
Soil Depth The total depth of the soil layer is equal to the value of laydeph for its deepest component. For example, for map unit TX560, Component 1 has depth 0 - 6 inches for layer 1 and 6 - 62 inches for layer 2 so the total depth is 62 inches (or slightly more than 5 feet) for this component.
The average depth can be found as:
Average Depth = Sum over all components of (comppct/100) * component depth
Available Water Capacity
In the layer table, the Available Water Capacity is given in dimensionless units of inches of water / inch of soil. In soil water analysis, we need to know the Water Holding Capacity of the soil, that is, the total depth of water that the soil could store within its whole vertical profile. To calculate this, we first find an average value of available water capacity for each soil layer by averaging the low and high values for that layer, (awcl + awch)/2, multiply the result by the layer thickness (laydeph - laydepl), and sum over all the layers:
Water Holding Capacity of each Component = sum over its layers of (awcl + awch)/2 * (laydeph - laydepl)
This gives a result in inches of water. For example, for mapunit TX560, Component 1 is a Tinn Clay soil having two layers, each of which has awcl = 0.15, and awch = 0.20, so the average awc = (0.15 + 0.20)/2 = 0.175. Layer 1 is 6 inches thick so its water holding capacity = 0.175 * 6 = 1.05 inches, and similarly, Layer 2 being 62 - 6 = 56 inches thick, has a water holding capacity of 56 * 0.175 = 9.80 inches of water. The total water holding capacity of the Tinn Clay component soil is 1.05 + 9.80 = 10.85 inches of water, whose storage is distributed throughout the pore spaces of the 62 inches of soil.
The average value of the water holding capacity of the soils in the TX560 map unit is found as a weighted average of the values in each component using the comppct percentages as weights:
Average Water Holding Capacity for the Map Unit = Sum over all its components of (comppct/100) * Water Holding Capacity of each ComponentComputing the average water holding capacity for a map unit is a fairly involved calculation so it is probably best to transfer the data to Excel and do the calculations there. If you click on a particular Table to highlight it, you can go to the File/Export Menu option and export the selected records of that Table (the ones in yellow) into a new Table. If you choose the .dbf in the Export dialog you will produce a dBase file that can be read directly by Excel. Exporting and transferring the required records from the comp and layer tables to a spreadsheet may save you time in transcribing the data values manually.
To be turned in: For mapunit TX560, how many layers does each component have? What is the total soil depth (inches) for each layer and the average depth (inches) for the map unit? What is the total water holding capacity (inches of water) over the full soil depth for each component? What is the average water holding capacity (inches of water)for soils in this map unit?
Ok, lets clean up now.
Enter Command: select mapunit
Enter Command: erase mapunit
Enter Command: select comp
Enter Command: erase comp
Enter Command: select layer
Enter Command: erase layer
Enter Command: quit
Arc: kill covshd10 all
Arc: kill soilshed all
Arc: kill rf1shed all