When the empirical Green's function method is applied to seismograms recorded at a site which is located on soft soil layers, it is necessary to take into account the nonlinear behavior of the soft soil layers. In practice, it has often been assumed that the seismic wave is affected by soil nonlinearity only after it's incidence to the local soft soil layers. If we consider a seismic ray connecting the source and the site, however, it might be true that the ray crosses the soft soil layers several times. Therefore, it is reasonable to assume that the seismic wave is affected by soil nonlinearity several times during the transmission from the source to the receiver. This potential phenomena is referred to as “multiple nonlinear effects” in this article. Taking into account these effects in the empirical Green's function method, the authors propose “nonlinear parameters” which represent the deviation of material properties of the sediments due to soil nonlinearity from linear status and which are used to modify Green's functions. The parameters are applied to the simulation of the 2000 Tottori-ken Seibu and the 1995 Hyogo-ken Nanbu earthquake ground motions. It is found that, although the new parameters are quite simple, they greatly contribute to the improvement of the simulation results.