Although not ideal, I understand that one option to include linear distributed springs (i.e., different linear p-y relationships along the pile) would be working with SubDyn. In this case, you should include in SubDyn the monopile (including the part below the seabed) and reaction joints with associated stiffness matrices.
Kind regards.