Scientific Investigations Report 2006–5318
U.S. GEOLOGICAL SURVEY
Scientific Investigations Report 2006–5318
basin_soilmoist_dpm.f
Calculates current soil moisture and soil saturation for the modeled area, minimum and maximum for the modeled area, and unsaturated and saturated deficits for each HRU and the modeled area.
Calculates current available water capacity and soil saturation for the modeled area, including the deficits for each soil moisture store and minimum and maximum HRU values for the modeled area. Module based on the soil moisture subroutine of Bauer and Vaccaro (1987) as modified to account for soil saturation by Bauer and Mastin (1997).
July, 2004
cov_type
HRU cover type: land use/cover type, from 1-31, no units. [basin]
hru_soil
HRU soil type, no units. [basin]
unsatmax (calculated)
Maximum available water capacity for a soil for a HRU, in inches.
stormax (calculated)
Maximum saturated water capacity for a soil for a HRU, in inches.
avlcap
Total available water capacity (field capacity) for a soil association, no units-as decimal fraction. [soilms]
spcyld
Specific yield for a soil association, no units-as decimal fraction. [soilms]
nlayer
Number of 6" layers for a soil association, no units. [soilms]
smsavg
Area-average of soil moisture capacity for all HRUs, in inches.
satavg
Area-average of soil saturation capacity for all HRUs, in inches.
smsdefct
Deficit of water in the unsaturated soil-moisture store for a HRU, in inches.
satdefct
Deficit of water in the saturated soil-moisture store for a HRU, in inches.
smscellmn
Current minimum soil moisture capacity for all HRUs, in inches.
smscellmx
Current maximum soil moisture capacity for all HRUs, in inches.
defctsms
Current deficit soil moisture capacity for all HRUs, in inches.
satcellmn
Current minimum saturated moisture capacity for all HRUs, in inches.
satcellmx
Current maximum saturated moisture capacity for all HRUs, in inches.
defctsat
Current deficit saturated moisture capacity for all HRUs, in inches.
wght
HRU area as fraction of total area, as a fraction. [basin]
unsatcur
Soil moisture in unsaturated soil moisture store for a HRU, in inches. [soilms]
storcur
Soil moisture in saturated soil moisture store for a HRU, in inches. [soilms]
The initialization part of the module determines the total (maximum) quantity of water that can be stored in the unsaturated store and the saturated store for each HRU; these calculations are for all land uses/covers except water bodies and impervious surfaces. For each HRU except those with covers of water or impervious surfaces, the current saturated soil moisture (storcur) is accumulated in the local variable satblk and the unsaturated soil moisture (unsatcur) is accumulated in the local variable smsblk. The area-weighted averages are then calculated using the variable wght and stored in the local variables satavg and smsavg. The maximum and minimum HRU values in the modeled area for saturation and available water capacity (the unsaturated store) are found and stored in the local variables satcellmx, satcellmn, smscellmx, and smscellmn, respectively. Concurrently, the HRU moisture deficit in the root zone is calculated for the saturated (satdefct) and unsaturated stores (smsdefct), and the basin averaged deficit for each store is calculated.
The deficit calculations are
smsblk = unsatcur(i)
satblk = storcur(i)
smsdefct(i) = smsblk - unsatmx(i)
satdefct(i) = satblk - stormax(i)
where
nhru is the number of HRUs,
i is the index for the HRU, from 1 to nhru,
nlayer is the number of layers for the association for this HRU,
smsblk is the current unsaturated water for a HRU,
satblk is the current saturated water for a HRU, and
other variables as previously defined.
Bauer, H.H., and Mastin, M.C., 1997, Recharge from precipitation in three small glacial-till mantled catchments in the Puget Sound Lowlands: U. S. Geological Survey Water-Resources Investigations Report 96-4219, 119 p.
Bauer, H.H., and Vaccaro, J.J., 1987, Documentation of a deep percolation model for estimating ground-water recharge: U. S. Geological Survey Open-File Report 86-536, 180 p.
Henry H. Bauer and John J. Vaccaro
U.S. Geological Survey
Washington Water Science Center
934 Broadway, Suite 300
Tacoma, WA 98402
Modified by:
John J. Vaccaro
U.S. Geological Survey
Washington Water Science Center
934 Broadway, Suite 300
Tacoma, WA 98402
Telephone: 253-552-1620
Fax: 253-552-1581
Email: jvaccaro@usgs.gov