Stiffeners play a vital role in strengthening thin panels in a wide range of engineering constructions by reducing additional structural weight. However, these structures are vulnerable to issues such as interlayer delamination or skin-stiffener interfacial debonding due to high stress levels developed from external environmental conditions and operational loadings. In contrast, ultrasonic-guided wave (UGW) techniques exhibit an efficient and precise approach for monitoring discontinuities or damages in composite structures. There is a lack of research on understanding the characteristics of the interaction between UGW and interfacial debonding when in-plane edge loading and environmental factors are simultaneously taken into account. Therefore, this study is motivated by the need to develop a multiphysics numerical model which employs a commercially available finite element software, COMSOL Multiphysics®, to simulate UGW propagation in a stiffened composite plate with debonding at the plate-stiffener interface through a piezoelectric transducer under the combined influence of in-plane edge load and hygrothermal environment. The stiffened plate and piezoelectric patches are modelled with the tetrahedral element, and the bottom surface of the attached stiffener has a through-width 0.1 mm deep groove simulated for debonding. The developed FE model is validated against the results of the conducted experiments and those found in the available literature through the correlation coefficient. Further, the study conducts a comprehensive parametric investigation on stiffened cross-ply (0/90/0) laminated plates, considering variations in debonding size, in-plane load, and hygrothermal load intensity through the excitation of A0 mode. The acquired response is processed to compare the peak amplitude of various modes and energy of the waveform. Additionally, statistical indices such as normalised correlation moment (NCM) and variance of the continuous wavelet transform (CWT) peak are estimated to understand the impact of various parameters on waveform. The results show that the presence of a 90° lamina in the cross-ply laminate generates a low amplitude S0 mode in the scattered response. Moreover, a mode conversion from A0 to S0 mode is observed due to perfect bonding between the plate and the stiffener, providing insights into the bonding state in the panel. Furthermore, it is found that the magnitude of the in-plane loading marginally affects the peak amplitude of various modes in the scattered response. Additionally, when temperature intensity rises, the energy and amplitude of the UGW signals acquired through piezoelectric patches positioned in a direct line with the actuator gradually increase. The NCM value enhances with debonding regardless of exposed hygrothermal condition and reduces with increasing temperature intensity. In addition, the variance of the CWT peak reduces with debonding. The findings of this research are expected to be helpful for the development of efficient algorithms for detecting damages for structural health monitoring of stiffened composite panels.