<abstract><p>This study investigates solutions for the time-fractional coupled system of the shallow-water equations. The shallow-water equations are employed for the purpose of elucidating the dynamics of water motion in oceanic or sea environments. Also, the aforementioned system characterizes a thin fluid layer that maintains a hydrostatic equilibrium while exhibiting uniform density. Shallow water flows have a vertical dimension that is considerably smaller in magnitude than the typical horizontal dimension. In the current work, we employ an innovative and effective technique, known as the natural transform decomposition method, to obtain the solutions for these fractional systems. The present methodology entails the utilization of both singular and non-singular kernels for the purpose of handling fractional derivatives. The Banach fixed point theorem is employed to demonstrate the uniqueness and convergence of the obtained solution. The outcomes obtained from the application of the suggested methodology are compared to the exact solution and the results of other numerical methods found in the literature, including the modified homotopy analysis transform method, the residual power series method and the new iterative method. The results obtained from the proposed methodology are presented through the use of tabular and graphical simulations. The current framework effectively captures the behavior exhibited by different fractional orders. The findings illustrate the efficacy of the proposed method.</p></abstract>