## 1 Introduction  
The electromagnetic fields excited by a dipole source in a layered dissipative medium structure have been investigated widely because of their applications in Through-the-Earth communication, underwater communication, ground penetrating radar technology, and antenna design [1; 2; 3; 4; 5; 6]. The developments in electromagnetic field propagation in layered medium structures have been analyzed by many investigators, including [PERSON] [7; 8; 9] and [PERSON] [7; 10; 11] who proposed the asymptotic methods, and the surface-impedance technique has been used for a two-layered region. Contour integrations and ranch cuts have been also used for the electromagnetic field in a layered region [8; 9]. [PERSON] et al. have obtained the complete formulas for electromagnetic fields excited by horizontal and vertical electric dipoles in two- and three-layered regions by using the proposed method in the recent studies. However, it is difficult to solve the arbitrary layered region problem. To overcome this drawback, [PERSON] et al. [12] presented a compact formulation of the electric-type and magnetic-type dyadic Green's functions for a plane-stratified, laterally unbounded multilayered structure, which is based on the tranmission-line network. [PERSON] et al. [13] has presented a near-fielding radiating dipole antenna next to a three-layered lossy sphere close to a human head model. Moreover, the numerical results have been verified via a unified method of moments (MoM) model based on the electric field integral equation (EFIE) given by [PERSON] [14]. These algorithms are suitable for analyzing the non-planar electromagnetic problem only when the structure is a concentric sphere, which limits their application in an arbitrary shaped non-planar layered region. [PERSON] et al. [15] studied the electromagnetic field propagation characteristic of shallow water in detail, and the influence of seabed layer is also considered.  
In most studies [16; 17; 18; 19], surface integral equation (SIE)-based methods have been adopted in the numerical analysis of the scattering of arbitrary shaped objects and the scattering of objects buried in layered structure, but few of them focus on the electromagnetic fields excited by dipoles in arbitrary layered dissipative media. In this paper, an SIE-based method to study the electromagnetic fields excited by electric dipole submerged in arbitrary layered homogeneous dissipative media is proposed. The total electromagnetic fields in each layer with three source components including the excitation source itself, equivalent electric/magnetic surface currents at the top of the layer, and equivalent electric/magnetic surface currents at the bottom of the layer are derived in detail. Then, the formulation proposed for the interface of each layer is discretized by using the Galerkin's method. After that, a matrix equation for a layered dissipative structure with a block tridiagonal matrix feature of a convenient solution is derived [20]. Three application scenarios are constructed and computed to analyze the electromagnetic characteristics with the vertical and horizontal dipole buried in wide, practically important media including sea, wet and other dissipative materials. Finally, the numerical results are presented to discuss the difference of the proposed method in comparison with the Computer Simulation Technology (CST) Studio Suite.