UW-Hydro / VIC

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

VIC 4.2.d active layer depth computation #901

Open rlilhare opened 4 years ago

rlilhare commented 4 years ago

Dear all: I am a little bit confused about how VIC 4.2.d calculates the active layer/thaw depths?

I am analyzing TDEPTH (depth of thawing front) output variable from the VIC model. Frozen soil is true in all VIC simulations, 5 thermal nodes, 10 m damping depth, and my study area fully comes under discontinuous permafrost region in north-central Canada. The model continuously runs from 1981 to 2070, which includes 5 years (1981-1984) of initial state condition.

I do NO FLUX: TRUE in the VIC simulations (use no flux boundary) to see the changes in ALT/Thaw depths in the projected climate.
This setting seems to give me decreased ALD in the future. The thaw depth/ALT computes nearly zero in winter but I see a decrease (negative changes) (the 1990s-2050s) in the ALT for my study area, which is odd to me.

When I do NO FLUX: FALSE (constant bottom boundary), the thaw depth/ALT computes 1.0-1.5 m in the winter for most of the grids in higher latitudes, which is also not true for winter? However, these simulations show a future increase in ALT/thaw depths (the 1990s-2050s).

Any expert advice on these results would be really helpful.

Thank you so much! Raj