In this article, we explore the submarine channel formation driven by the interaction of turbidity currents with an erodible bed. The theoretical analysis considers the three-dimensional continuity and momentum equations of the fluid phase, and the advection-diffusion and Exner equations of the solid phase. The governing equations are linearized by imposing periodic perturbations on the base flow. We study the response of both the base flow (profiles of velocity and suspended sediment concentration) and perturbations (growth rate and perturbation fields) to changes in key parameters related to the flow and sediment transport. The growth rate and the critical wavenumber are examined for a given quintet formed by the gravitational parameter, longitudinal bed slope, sediment concentration at the edge of the driving layer, Rouse number and erosion coefficient. The critical wavenumber reduces with an increase in gravitational parameter, longitudinal bed slope, sediment concentration at the edge of the driving layer and erosion coefficient, while it increases with the Rouse number. For the submarine channel formation, we identify the upper threshold values for the gravitational parameter, longitudinal bed slope, sediment concentration at the edge of the driving layer and erosion coefficient and the lower threshold value for the Rouse number. © 2022 The Authors.