The Carbon Capture Storage (CCS) is considered as a promising atmospheric CO₂ mitigation solution. There is some risk that the high concentrated carbon dioxide could be leaked accidentally from the CCS facilities. It requires making attention on the behavior of carbon dioxide (CO₂) in the ocean. In particular, how such substances behave in the different phases (such as gaseous, liquid, hydrate or non-hydrate) should be understood with the forms of existences as like as in bubble/droplet form or CO₂ lake. This study focused on the characteristics of liquid CO₂ jet/plume in a near-field area mainly and investigated the behavior of a liquid CO₂ jet/ plume released at 800 m depth using a series of numerical simulation. The numerical results were verified by comparing with the theoretical and previous experiments and showed a good agreement in appearance of bifurcation as well as the trajectory of the flow centerline. The present work described the development of CO₂ flow, a turbulent buoyant fluid driven by concentration difference, both with time and spatial aspects. The dependency of jet/plume cross sectional areas on the Reynolds number, the centerline density distribution, the completely dissolving distance, and the characteristic of turbulent intensity and turbulent dissipation rate were also investigated.