Theseaicecoverdisplaysvariousdynamicalcharacteristicssuchasbreakup,rafting,andridgingunderexternalforces.Tomodeltheicedynamicprocessaccurately,theeffectivenumericalmodelingmethodshouldbeestablished.Inthispaper,amodifiedparticle-in-cell(PIC)methodforseaicedynamicsisdevelopedcouplingthefinitedifference(FD)methodandsmoothedparticlehydrodynamics(SPH).Inthismethod,theicecoverisfirstdiscretizedintoaseriesofLagrangianiceparticleswhichhavetheirownsizes,thicknesses,concentrationsandvelocities.TheicethicknessandconcentrationatEuleriangridpositionsareobtainedbyinterpolationwiththeGaussianfunctionfromtheirsurroundingiceparticles.ThemomentumoficecoverissolvedwithFDapproachtoobtaintheEuleriancellvelocity,whichisusedtoestimatetheiceparticlevelocitywiththeGaussianfunctionalso.Thethicknessandconcentrationoficeparticlesareadjustedwithparticlemassdensityandsmoothlength,whichareadjustedwiththeredistributionoficeparticles.WiththeabovemodifiedPICmethod,numericalsimulationsforicemotioninanidealizedrectangularbasinandtheicedynamicsintheBohaiSeaarecarriedout.ThesesimulationsshowthatthismodifiedPICmethodisapplicabletoseaicedynamicssimulation.