UW-Hydro / VIC

The Variable Infiltration Capacity (VIC) Macroscale Hydrologic Model
http://vic.readthedocs.io
MIT License
259 stars 385 forks source link

VIC 5.0.1: Negative OUT_EVAP spikes in winter, and snow in canopy in summer #886

Open BenjaminFTurner opened 5 years ago

BenjaminFTurner commented 5 years ago

This is a moderate temperate climate in the northern hemisphere (West Virginia) where temperatures never get below freezing in the summer months. However, the model is showing continuously positive values of OUT_SNOW_CANOPY from the middle of May to the end of September, which obviously is not physically realistic. Also there are negative spikes in OUT_EVAP as large as -70 mm in winter months, which seem to screw up the water balance substantially.

Using the more computationally expensive energy balance options decreases the values of OUT_SNOW_CANOPY in summer months, but it doesn't go away.

A subset of my current settings is as follows. FULL_ENERGY TRUE CLOSE_ENERGY TRUE QUICK_FLUX FALSE FROZEN_SOIL TRUE TFALLBACK TRUE AERO_RESIST_CANSNOW AR_406_FULL NODES 10 SNOW_BAND FALSE

MAX_SNOW_TEMP and MIN_RAIN_TEMP use the defaults.

I'm using VIC 5.0.1 with the Image driver, gcc, Ubuntu 18.04.

I'm also attaching time series of basin-averaged OUT_SNOW_CANOPY and OUT_EVAP, as well relative humidity calculated from input VP and temperature (which never goes above 100%), and a complete set of inputs.

I'm wondering if this is a quirk/bug of VIC, or could any of my inputs be causing this strange behavior? I've been racking my brain on this for a while, but haven't figured anything out.

Ben OUT_SNOW_CANOPY OUT_EVAP RH wf_vic_inputs.zip