diff --git a/src/Physics/Surface/SEB/seb.jl b/src/Physics/Surface/SEB/seb.jl
index b5245cdb2e20edef8f3b91fd519ebd8f31f49e2f..0db86369721863149c4de5b689b254e3ee27fcc3 100644
--- a/src/Physics/Surface/SEB/seb.jl
+++ b/src/Physics/Surface/SEB/seb.jl
@@ -48,8 +48,19 @@ Utils.@properties SEBParams(
     # surface propertiess; includes soil and snow by default;
     # this can be easily extended by the user to include other layers by adding more
     # fields to `SEBParams` and then implementing `surfaceproperties` for those layers.
-    soil = SurfaceProperties(),
-    snow = SurfaceProperties(α=0.8, ϵ=0.99, z₀=5e-4),
+    soil = SurfaceProperties(
+        α = param(0.2, desc="Initial snow-free surface albedo"),
+        ϵ = param(0.97, desc="Initial snow-free surface emissivity"),
+        zâ‚€ = param(1e-3, units=u"m", desc="Initial snow-free surface roughness length"),
+        râ‚› = param(50.0, units=u"s/m", desc="Initial snow-free surface resistance against evapotranspiration and sublimation"),
+    ),
+    
+    snow = SurfaceProperties(
+        α = param(0.8, desc="Initial snow-covered surface albedo"),
+        ϵ = param(0.99, desc="Initial snow-covered surface emissivity"),
+        zâ‚€ = param(5e-4, desc="Initial snow-covered surface roughness length"),
+        râ‚› = param(50.0, units=u"s/m", desc="Initial snow-covered surface resistance against evapotranspiration and sublimation"),
+    ),
 
     # "natural" constant
     σ = 5.6704e-8u"J/(s*m^2*K^4)",   # Stefan-Boltzmann constant