We present the solar-terrestrial transit process of three successive coronal mass ejections (CMEs) of November 4―5, 1998 originating from active region 8375 by using a time-dependent three-dimensional magnetohydrodynamics (MHD) simula-tion. These CMEs interacted with each other while they were propagating in inter-planetary space and finally formed a “complex ejecta”. A newly developed SIP-CESE MHD model was applied to solve MHD equations numerically. The quiet solar wind was started from Parker-like 1D solar wind solution and the magnetic field map was calculated from the solar photospheric magnetic field data. In our simulation, the ejections were initiated using pulse in the real active region 8375. The interplanetary disturbance parameters, such as speed, direction and angular size of the expanding CME, were determined from the SOHO/LASCO data with the cone-model. We discussed the three-dimensional aspects of the propagation, in-teraction and merging of the three ejections. The simulated interplanetary shocks were compared with the nearby-Earth measurement. The results showed that our simulation could reproduce and explain some of the general features observed by satellite for the “complex ejecta”. We present the solar-terrestrial transit process of three successive coronal mass ejections (CMEs) of November 4-5, 1998 originating from active region 8375 by using a time-dependent three-dimensional magnetohydrodynamics (MHD) simula- tion. These CMEs interacted with each other while they were propagating in inter-planetary space and finally formed a “complex ejecta”. A newly developed SIP-CESE MHD model was applied to solve MHD equations numerically. The quiet solar wind was started from Parker-like 1D solar In our simulation, the ejections were initiated using pulses in the real active region 8375. The interplanetary disturbance parameters, such as speed, direction and angular size of the expanding We discussed the three-dimensional aspects of the propagation, in-teraction and merging of the three ejections. The simulated interplanetary shocks were compared with the nearby-Earth measurement. The results showed that our simulation could reproduce and explain some of the general features observed by satellite for the “complex ejecta”.