Response in a porous seabed under dynamic environmental loading is a vital engineering issue in marine geotechnics. Lots of investigations for seabed response under dynamic loading have been developed through mathematical, numerical and experimental approaches. Most previous numerical models for seabed response in marine environments were based on finite element models. In this paper, based on local radial basis function collection method (LRBFCM), a meshfree model is proposed for the seabed response in the marine environments. In the present model, partial dynamic approximation (u- p approximation) will be used, and three different types of natural loading will be considered, i.e., wave, current and earthquake loading.