diff --git a/src/biogeophys/TotalWaterAndHeatMod.F90 b/src/biogeophys/TotalWaterAndHeatMod.F90
index b774fe9761..79ee0a77d6 100644
--- a/src/biogeophys/TotalWaterAndHeatMod.F90
+++ b/src/biogeophys/TotalWaterAndHeatMod.F90
@@ -518,7 +518,7 @@ subroutine AccumulateLiqIceMassLake(bounds, num_c, filter_c, &
           c = filter_c(fc)
           ! calculate lake liq and ice content per lake layer first
           h2olak_liq = dz_lake(c,j) * denh2o * (1 - lake_icefrac(c,j))
-          h2olak_ice = dz_lake(c,j) * denice * lake_icefrac(c,j)
+          h2olak_ice = dz_lake(c,j) * denh2o * lake_icefrac(c,j) ! use water density of liquid water as layer depth is not adjusted
           
           liquid_mass(c) = liquid_mass(c) + h2olak_liq
           ice_mass(c) = ice_mass(c) + h2olak_ice
@@ -1086,7 +1086,8 @@ subroutine AccumulateHeatLake(bounds, num_c, filter_c, &
                 heat_liquid = lake_heat_liquid(c), &
                 latent_heat_liquid = lake_latent_heat_liquid(c))
             ! ice heat
-            h2olak_ice = dz_lake(c,j) * denice * lake_icefrac(c,j)      
+            ! use water density as lake layer does not adjust
+            h2olak_ice = dz_lake(c,j) * denh2o * lake_icefrac(c,j)      
             lake_heat_ice(c) = lake_heat_ice(c) + & 
                 TempToHeat(temp=t_lake(c,j), cv = (h2olak_ice * cpice))
         end do