本報では,ため池堤体の越流侵食過程を安定的に解くことができる数値解析手法を紹介する。同手法では越流水の流れは浅水流方程式を有限体積法によって解き,それから得られる底面摩擦力から侵食速度を用いて侵食による堤体表面の形状変化を逐次的に解析する。しかし,堤体上を流れる越流水は,その堤体形状の変化に伴って急激に変化するため,安定的な侵食計算が困難な場合が多い。ここでは,侵食の安定化を図るための一手法として,侵食面のラプラシアン(凹凸)を考慮した解法を提案する。本報では,提案する解析手法を簡潔に説明し,その手法によって得られる解析結果を報告する。