A mathematical model that comprehensively captures the real behavior of riverbed deformation, encompassing all pertinent effects, is developed. The underwater slope reformation process, with the generatrix aligned along the flow velocity in the model, is considered. A numerical model is introduced to calculate the flow involving a deformable bottom, and the model's validation is established through rigorous analysis of experimental findings. The research firmly confirms the suitability of the proposed mathematical and numerical model for describing deformations in uneven and unsteady river flows, including the movement of dredging slots and channel quarries. The model's minimal equation count and reliance on empirical constants demonstrate the model's efficiency. The model's predictions align strongly with experimental data, although optimal values of empirical coefficients vary slightly across different experiments. Hence, there is a call for further investigation to derive more universally applicable closure relationships for the model. The importance of validating the model with reliable field data and its potential extension to accommodate hydraulically diverse soils is emphasized. Such an extension is feasible due to the concentration transfer equation, enabling independent calculations for particle fractions of varying sizes as long as the total particle concentration in the stream remains within reasonable limits. This dedicated research contributes significantly to understanding riverbed deformations and advancing accurate modeling and management of riverine environments.