Noroozi, M., Riahi, M.A. and Bostani, A., 2023. Numerical study of combined modulated electromagnetic wave parameters in forward modeling of reflected seismic waves from subsurface layers. Journal of Seismic Exploration, 32: 357-372. Nowadays, a combination of different methods are implemented for identifying subsurface reserves with the aim of increasing accuracy and reducing failure costs. In exploration geophysics, electromagnetic methods with a controlled source are used in and and sea environments to study subsurface structures. On the other hand, seismic methods are traditionally the most common method for hydrocarbon explorations. In this study, a combination of the two methods is used. As usual for solving electromagnetic wave equations, the dielectric coefficient parameter is assumed to be constant. But here the electromagnetic wave equations are solved by considering the dielectric coefficient parameter, resulting in a higher accuracy in identifying the subsurface layers. The importance of the dielectric coefficient parameter is that the value of this parameter for oil and gas with its value for the host rock surrounding the oil and gas trap is significantly different, therefore it facilitates the exploration process. Also in our numerical modeling, we used a wave packet which includes a sinusoidal wave and a square wave that is propagated into the Earth model in a modulated form. Those waves have different wavelengths because waves with long wavelengths can penetrate deeper while waves with short wavelengths result in more resolution. 0963-065 1/23/$5.00 © 2023 Geohysical Press Ltd. 358 We used the Marmousi2 model which is made up of several hydrocarbon sections including the Marmousi hydrocarbon section. Then, using the empirical relationship between seismic velocity and the resistivity parameter we extracted the dielectric coefficient parameter. According to the lithological information we produced a dielectric coefficient matrix. Finally, with this information, we made a numerical forward modeling code to generate the final output of this model which is the intensity of the ground- induced reflected fields received by the receivers. Then a forward model for the seismological method was prepared and the output related to the reflection coefficient was obtained for this model. The comparison of these two outputs with 3 methods, consisting of the mean squared error, structural similarity index, and histogram chart. we observe that when using the joint method, the results of the properties of the layers are obtained more accurately than when using the conventional seismic method alone.