USGS Home Page

USGS Open-File Report 02-362


Changes from the original module are written in bold text

NAME

precip_prms_dist2.f

NAME-OLD

precip_prms.f

MODULE PROCESS (TYPE)

Precipitation form and distribution

DEFINITION

This module determines whether measured precipitation is rain or snow and distributes it to the HRUs. The daily precipitation values for the weather sites are interpolated to each HRU by an inverse distance weighting scheme and then averaged.

KEYWORDS

CREATION DATE

March 2001

PARAMETERS DECLARED

adjmix_rain

Monthly factor to adjust rain proportion in a mixed rain/snow event.

dist2

The square of the distance between an HRU and each weather station. Does not need to be input because it is calculated. Identified as parameter because not a state variable. Units are the square of the input unit for hru_xlong, hru_ylat, psta_xlong, and psta_ylat parameters (see below).

hru_xlong

The x or longitude of the centroid of an HRU that is input in projection coordinates such as State Plane, Lambert, or UTM. Used to calculate distance between HRU and a weather site. Units can be feet, miles, meters, or kilometers but must be consistent with hru_ylat, psta_xlong, and psta_ylat.

hru_ylat

The y or latitude of the centroid of an HRU that is input in projection coordinates such as State Plane, Lambert, or UTM. Used to calculate distance between HRU and a weather site. Units can be feet, miles, meters, or kilometers but must be consistent with hru_xlong, psta_xlong, and psta_ylat.

psta_mon

Mean monthly precipitation for each weather station, in inches. Mean monthly values should be consistent with period of record used for rain_mon and snow_mon described below.

psta_xlong

The x or longitude of each weather station that is input in projection coordinates such as State Plane, Lambert, or UTM. Used to calculate distance between HRU and a weather site. Units can be feet, miles, meters, or kilometers but must be consistent with hru_ylat, hru_xlong, and psta_ylat.

psta_ylat

The y or latitude of each weather station that is input in projection coordinates such as State Plane, Lambert, or UTM. Used to calculate distance between HRU and a weather site. Units can be feet, miles, meters, or kilometers but must be consistent with hru_xlong, psta_xlong, and hru_ylat.

maxmon_prec

Maximum daily precipitation that can occur in each month for each weather station. If precipitation at a station is greater than this value for the month, it is assumed to be in error. Value is used to constrain bad or missing values, in inches.

rain_mon

Mean monthly HRU rain to account for differences in elevation, spatial variation, and topography, in inches. Can be obtained from the National Weather Service's spatial distribution of mean annual precipitation for the 1961-90 climate normal period.

snow_mon

Mean monthly HRU snow to account for differences in elevation, spatial variation, and topography, in inches. Can be obtained from the National Weather Service's spatial distribution of mean annual precipitation for the 1961-90 climate normal period.

tmax_allrain

If HRU maximum temperature is greater than or equal to this value precipitation is assumed to be rain, by month, in ° C or ° F, depending on units of data.

tmax_allsnow

If HRU maximum temperature is less than or equal to this value, precipitation is assumed to be snow, by month, in ° C or ° F, depending on units of data.

tstorm_mo

Monthly indicator for prevalent storm type, 0 indicates frontal storms prevalent, 1 indicates convective storms prevalent during month.

basin_area

Mean elevation for each HRU, in feet. [basin]

hru_area

Area for each HRU, in acres. [basin]

temp_units

Indicator for units for temperature data, 0= ° F and 1= ° C. [basin]

VARIABLES DECLARED

basin_ppt

Average basin precipitation, in inches.

hru_ppt

Precipitation on HRU, rain and snow, in inches.

hru_rain

Rain on HRU, in inches.

hru_snow

Snow on HRU, in inches.

newsnow

Indicator for new snow during time step.

pptmix

Indicator for mixed rain and snow during time step.

prmx

The proportion of rain in a mixture of rain and snow.

EXTERNAL VARIABLES USED

precip

Observed precipitation at each measurement station, in inches. [obs]

form_data

Indicator to override computed determination of the precipitation form, rain or snow. [obs]

tmax

Observed maximum temperature at each temperature measurement station, ° F or ° C, depending on units of data. [obs]

tmaxf or tmaxc

Maximum HRU temperature, ° F or ° C, depending on units of data. [temp]

tminf or tminc

Minimum HRU temperature, ° F or ° C, depending on units of data. [temp]

DESCRIPTION

Total daily precipitation depth (hru_ppt) received on an HRU is computed by

 

hruppt = (Summation((pcor multiplied by precip) divided by (dist2))) divided by (Summation dist2)  

 

where

pcor = rain_mon/psta_mon if precipitation is rain, or

pcor = snow_mon/pstat_mon if precipitation is snow,

precip = observed precipitation at a measurement station corresponding to the HRU, and

dist2 = the square of the distance between the station and the centroid of the HRU.

 

The above method of interpolating precipitation from a weather site to a HRU using the inverse of the square of the distance between the site and an HRU is from Dean and Snyder (1977). The correction factor (pcor) is based on the method of Bauer and Vaccaro (1987), which used mean annual values, and has been modified to use mean monthly values.

Precipitation form (rain, snow, or a mixture of both) on each HRU is estimated from the HRU maximum and minimum daily air temperatures and their relationship to a base temperature (tmax_allsnow). Precipitation is all snow if the maximum temperature is less than or equal to the tmax_allsnow and all rain if the minimum temperature is greater than or equal to tmax_allsnow . If the maximum temperature is greater than tmax_allsnow and the minimum temperature is below tmax_allsnow, then the precipitation is considered a mixture, and the rain is assumed to occur first. The portion of the total precipitation occurring as rain (prmx) is computed by

 

prmx = [(tmax minus tmax_allsnow) divided by (tmax minus tmin)] multiplied by adjmix_rain  

 

where

tmax is the maximum HRU temperature,

tmin is the minimum HRU temperature, and

adjmix_rain is a monthly factor to adjust the rain proportion in a mixed rain/snow event.

 

This mixture algorithm can be overridden in two ways. One is the use of the parameter tmax_allrain, which is an air temperature value that, when exceeded by tmax, forces the precipitation to be considered all rain. This parameter is useful for periods, such as in the spring, when the minimum daily temperatures may be below tmax_allsnow but precipitation is predominantly convective afternoon storms. The form of the precipitation may also be explicitly specified by including the variable form_data in the observed data file.

REFERENCES

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.

Dean, J.D., and Snyder, W.M., 1977, Temporally and areally distributed rainfall: Journal of Irrigation Division, American Society of Civil Engineers, TA 103, No. IR2, p. 221-224.

Leavesley, G.H., Lichty, R.W., Troutman, B.M., and Saindon, L.G., 1983, Precipitation-runoff modeling system--User's manual: U. S. Geological Survey Water-Resources Investigations Report 83-4238, 207 p.

Willen, D.W., Shumway, C.A., and Reid, J.E., 1971, Simulation of daily snow water equivalent and melt, in Western Snow Conference, Billings, Mont., 1971, Proceedings: v. 39, p. 1-8.

DEVELOPER NAME AND ADDRESS

 

George H. Leavesley

U.S. Geological Survey, WRD

Box 25046, MS 412, DFC

Denver, CO 80225

 

Telephone: 303-236-5026

Fax: 303-236-5034

Email: george@usgs.gov

 

Modified by:

John J. Vaccaro

U.S. Geological Survey, WRD

1201 Pacific Ave., Suite 600

Tacoma, WA 98402

 

Telephone: 253-428-3600 x2620

Fax: 253-428-3614

Email: jvaccaro@usgs.gov

 


URL for this page is http://pubsdata.usgs.gov/pubs/of/2002/ofr02362/htdocs/precip/precip_prms_dist2.htm
Page contact: Mark Mastin (mcmastin@usgs.gov), 253-428-3600, ext. 2609
Last modified: Friday, 11-Jan-2013 03:19:50 EST

Bobby Approved Valid HTML 4.01!