We derive two reaction–diffusion systems modeling the impact of a nuclear waste contamination on food chain population sets. We exhibit the numerical difficulties the problem features and provide second order accurate nonnegative approximate solutions in two space dimensions.