Frozen ground can be important to flood production and is often heterogeneous within a watershed due to spatial 10 variations in the available energy, insulation by snowpack and ground cover, and the thermal and moisture properties of the soil. The widely-used Continuous Frozen Ground Index (CFGI) model is a degree-day approach and identifies frozen ground using a simple frost index, which varies mainly with elevation through a temperature-elevation relationship.Similarly, snow depth and its insulating effect are also estimated based on elevation. The objective of this work is to more accurately represent the spatial heterogeneity of frozen ground in a distributed hydrologic model, Gridded Surface 15 Subsurface Hydrologic Analysis (GSSHA), by modifying the CFGI method. Among the modifications, the snowpack and frost indices are simulated by replacing air temperature (a surrogate for the available energy) with a radiation-derived temperature that aims to better represent spatial variations in available energy. Ground cover is also included as an additional insulator of the soil. Furthermore, the modified Berggren Equation, which accounts for soil thermal conductivity and soil moisture, is used to convert the frost index into frost depth. The modified CFGI model is tested by application at six 20 test sites within the Sleepers River Experimental Watershed in Vermont. Compared to the CFGI model, the modified CFGI model more accurately captures the variations in frozen ground between the sites, inter-annual variations in frozen ground depths at a given site, and the occurrence of frozen ground.