Nonlinear evolution equations are typically used to model and analyse the majority of natural processes. These equations are extremely important for comprehending many intricate situations involving wave propagation, including heat transfer, fluid dynamics, optical fibres, electrodynamics, physics, chemistry, biology, condensed matter physics, ocean engineering, and many other nonlinear science fields. In this study, the anticipation of a bilinear Kdv model known as the (1+1)-dimensional integro-differential Ito equation is considered which depicts the oceanic shallow water wave dynamics. Firstly, we determine the bilinear form of the model and then we develop lump waves and analyse their interactions with periodic waves using symbolic computation approaches. Later, the interactions between single, double-kink and lump waves solutions, as well as interactions between periodic, lump and soliton solutions is considered. Furthermore, the ( G ′ / G 2 ) expansion approach and the Adomian technique are implemented to acquire a variety of novel configurations for the governing dynamical model. Later, the corresponding absolute accuracy is evaluated through a comparison of several types of observed solutions. Additionally, the stability of the solitary travelling wave solutions to the specified model is also established. The 3D, 2D, and contour plots have been employed to demonstrate the characteristics of the nonlinear model and aid in determine the proper set of parameters. As a result, a collection of options for bright, dark, periodic, rational, and elliptic functions is established. To the best of our knowledge, the current work presents a novel case study that has not been previously studied in order to generate several new solutions to the integro-differential equation appearing in shallow water waves. The results show that the strategies that have been employed are more effective and capable than the traditional methods found in previous research. • The interactions between single, double-kink and lump waves solutions is considered. • The (G'/G2) expansion and the Adomian techniques are implemented for the dynamical model. • We have established a comparison between analytical and semi-analytical solutions. • The 3D, 2D, and contour plots have been employed in determine the proper set of parameters. [ABSTRACT FROM AUTHOR]