We have developed a new, simple Eulerian–Lagrangian approach to solve equations for sinking particulate organic matter in the ocean. We rely on the known parameterizations, but our approach to solving the problem differs, allowing the algorithm to be incorporated into biogeochemical global ocean models with relative ease. New analytical and numerical solutions have confirmed that feedback between the degradation rate and sinking velocity significantly alters particulate matter fluxes.