Numerical modelling of heavy metals for riverine systems using a new approach to the source term in the ADE