The flow characteristics of a 2D slot jet vertically impinging on free surfaces are numerically investigated by the Large Eddy Simulation (LES) with a dynamic Sub-Grid Scale (SGS) model. The σ-coordinate transformation is introduced to map the depth-variable physical domain to a depth-uniform computational domain. The split-operator scheme, which splits the solution procedure into advection, diffusion and pressure propagation steps, is employed to solve the instantaneous velocity and pressure field. A fully nonlinear Lagrange-Euler method is used to compute the free surface elevation. The numerical results show that the jet retains good self-similarity in the Zone of Eestablished Flow (ZEF), while the turbulence of jet shifts gradually from isotropic to anisotropic in the Zone of Surface Impingement (ZSI). When the jet approaches the free surface, the centerline velocity decreases rapidly and the flow is deflected by the free surface. Two symmetrical surface jets are formed in the Zone of Horizontal Jets (ZHJ). An approximate Guassian distribution of super-elevation is also formed on the free surface. The computed results are in good agreement with the available experimental data. The flow characteristics of a 2D slot jet vertically impinging on free surfaces are numerically investigated by the Large Eddy Simulation (LES) with a dynamic Sub-Grid Scale (SGS) model. The σ-coordinate transformation is introduced to the map the depth-variable physical The split-operator scheme, which splits the solution procedure into advection, diffusion and pressure propagation steps, is employed to solve the instantaneous velocity and pressure field. A fully nonlinear Lagrange-Euler method is used to The numerical results show that the jet retains good self-similarity in the Zone of Eestablished Flow (ZEF), while the turbulence of jet shifts gradually from isotropic to anisotropic in the Zone of Surface Impingement (ZSI). When the jet approaches the free surface, the centerline velocity decreases rapidly and the flow is deflected by the free surface. Two symmetrical surface jets are formed in the Zon e of Horizontal Jets (ZHJ). An approximate Guassian distribution of super-elevation is also formed on the free surface. The computed results are in good agreement with the available experimental data.