A mass balance based soil heavy metal accumulation model was developed.
Trends of soil Cd content under five remediation scenarios was predicted.
Liming reduces Cd content in rice produce but increases soil Cd in the long term.
Higher inputs of Cd in the past were estimated based on retrospective simulations.
Integrating remediation measures reduced Cd in both soils and rice grains.