The sediment suspension and transport process under complex breaking wave situation is investigated using large eddy simulation(abbreviated as LES hereafter) method. The coupled level set(LS) and volume of fluid(VOF) method is used to accurately capture the evolution of air-water interface. The wall effect at the bottom is modeled based on the wave friction term while the complicate bottom boundary condition for sediment is tackled using Chou and Fringer’s sediment erosion and deposition flux method. A simulation is carried out to study the sediment suspension and transport process under periodic plunging breaking waves. The comparison between the results by CLSVOF method and those obtained by the LS method is given. It shows that the latter performs as well as the CLSVOF method in the pre-breaking weak-surface deformation situation. However, a serious mass conservation problem in the later stages of wave breaking makes it inappropriate for this study by use of the LS method and thus the CLSVOF method is suggested. The flow field and the distribution of suspended sediment concentration are then analyzed in detail. At the early stage of breaking, the sediment is mainly concentrated near the bottom area. During the wave breaking process, when the entrapped large-scale air bubble travels downward to approach the bottom, strong shear is induced and the sediment is highly entrained. The coupled suspension set-up and transport process under complex breaking wave situation is investigated using large eddy simulation (abbreviated as LES hereafter) method. The coupled level set (LS) and volume of fluid (VOF) method is used to accurately capture the evolution of air- water interface. The wall effect at the bottom is modeled based on the wave friction term while the complicate bottom boundary condition for sediment is tackled using Chou and Fringer’s sediment erosion and deposition flux method. A simulation is carried out to study the sediment suspension and transport process under periodic plunging breaking waves. The comparison between the results by CLSVOF method and those obtained by the LS method is given. It shows that the latter performs well as the CLSVOF method in the pre-breaking weak- surface deformation situation. However, a serious mass conservation problem in the later stages of wave breaking makes it inappropriate for this study by use of the LS method and thus The flow field and the distribution of suspended sediment concentration are then analyzed in detail. At the early stage of breaking, the sediment is mainly concentrated near the bottom area. During the wave breaking process, when the entrapped large- scale air bubble travels downward to approach the bottom, strong shear is induced and the sediment is highly entrained.