Wind-wave-current　interaction　is　a　phenomenon　which　is　present　in　many　practical　engineering　situations　and　vertical　earthquakes　may　have　significant　influence　on　the　hydrodynamic　pressure　of　seawater　in　regions　of　active　seismicity.　This　study　is　devoted　to　probe　the　wave　height　and　hydrodynamic　pressure　of　seawater　caused　by　the　joint　wind,　wave,　current　and　earthquake　action　based　on　computational　fluid　dynamics　(CFD)　method,　where　a　homogeneous　multiphase　model　based　on　the　noncompressible　Reynolds　Averaged　Navier-Stokes　(RANS)　equation　and　the　k-epsilon　turbulence　model　implemented　in　ANSYS-Fluent　code　was　used.　First,　a　two-dimensional　(2D)　numerical　flume　excited　by　wind,　wave,　current　and　earthquake　loadings　was　established　through　the　secondary　development　of　a　user-defined　function　(UDF).　Subsequent,　effects　of　wind-wave　interaction　and　wind-wave-current　interaction　on　the　wave　height　and　hydrodynamic　pressure　of　seawater　were　investigated.　Finally,　effects　of　wave-earthquake　interaction　and　wave-current-earthquake　interaction　on　the　wave　height　and　hydrodynamic　pressure　of　seawater　were　investigated.　The　nonlinearity　of　seawater　under　joint　wind,　wave,　current　and　earthquake　action　was　also　studied.