Uncertainty in the simulation of runoff due to the parameterization of frozen soil moisture using the Global Soil Wetness Project methodology

Four simulations of the region 30 degrees N-90 degrees N are performed usin g the Global Soil Wetness Project methodology and a single land surface sch eme. Four methods are used to represent soil ice: an explicit representatio n of the thermal and hydrological effects of soil ice; two implicit methods (which only account for the hydrological effects); and finally the simples t approach where soil ice is not accounted for. Substantial impacts on tota l runoff, evaporation and temperature result from the choice of parameteriz ation. The partitioning of total runoff between drainage and surface runoff is also changed. The impacts on temperature are large enough to cause prob lems for "fingerprinting" of global change while the change in the runoff g eneration process, and the timing of maximum runoff are large enough to con cern ocean modelers. Evidence presented here and elsewhere indicates that l and surface schemes should include the thermal effects of soil ice melting and freezing. However, the hydrological effects of soil ice suppressing inf iltration and encouraging surface runoff may be based on observations taken at a scale inappropriate to climate model parameterization. We show that f or one basin, the Mackenzie, a land surface model which ignores soil ice en tirely simulates runoff better than the other methodologies tested here. We therefore hypothesize that it may be preferable to not include soil ice in the runoff formulations used in land surface models until we have more obs ervations at an appropriate spatial scale. Testing of frozen soil moisture parameterizations in other catchments with high-quality observed runoff dat a should be conducted to test this hypothesis.