A higher-order moving particle semi-implicit (MPS) method was developed to solve water entry problems. The Wendland kernel function was applied in the particle interaction model. Various models for pressure gradient were investigated. To overcome the inconsistency in the original MPS methods, a pressure gradient correction was implemented to guarantee the first-order consistency of gradient. The corrective matrix was modified by using the derivative of the kernel function. A particle shifting technique was also applied to improve the numerical stability. Validation studies were carried out for water entry of a rigid wedge with the tilting angles of 0°, 10° and 20°, and a rigid ship section. The solutions by the present method are generally in good agreement with experimental data and other published numerical results.