We develop a single-class ice and snow model embedded inside a 3D hydrodynamic model on unstructured grids and apply it to lake studies using highly variable mesh resolution. The model is able to reasonably capture the ice fields observed in both small and large lakes. For the first time, we attempt simulation of ice processes on very small scales (~ 1 m). Physically sound results are obtained at the expense of moderately increased computational cost, although more rigorous validation nearshore is needed due to lack of observation. We also outline challenges on developing new process-based capabilities for accurately simulating nearshore ice.