The transient behavior of the immiscible two liquids interface, which is initiated by the rising gas bubble was investigated using Smoothed Particle Hydrodynamics (SPH) model. This developed numerical method is using fully Lagrangian particle-based model, which can track the movement of both the gas and the liquid phase directly. Numerical simulations have been performed for different conditions corresponding to different values of surface tension, viscosity and density, and the predicted topological changes as well as the theoretical pressure and interfacial shape of bubbles are validated. In the case of immiscible two liquids, the column of a lower liquid phase penetrating into the upper liquid phase influences interface area, whose shape strongly depends on the wake flow pattern of a bubble. Thus, the dynamic balance between the buoyancy and the liquid-liquid interfacial tension determines an interface area. Under higher surface tension condition, such as molten metal-slag system, the liquid-liquid interface shape is greatly influenced by the fluctuation of a bubble. Then, nonlinearly changed interface shape can be observed, indicating that this shape becomes easily unstable by slight change of the curvature. (C) 2015 Elsevier Ltd. All rights reserved.