The removal mechanisms of contaminants in Fe0/H2O are still poorly understood, and characterized by contradictory findings. Therefore, the present study aims to improve the understanding of the processes involved in phosphate removal in Fe0/H2O system. Herein, the methylene blue method (MB method) is used to trace the dynamics within the investigated systems. The MB method utilizes the differential adsorptive affinity of MB onto sand and sand coated with iron corrosion products to evaluate the degree of Fe0 corrosion in Fe0/H2O systems. The extent of MB discoloration and phosphate removal in various Fe0-based systems was characterized in parallel quiescent batch experiments for two weeks and two months. Parallel experiments with Orange II (O-II) as a model contaminant allowed an improved discussion of the results. The investigated systems were: (i) Fe0 alone, (ii) MnO2 alone, (iii) sand alone, (iv) Fe0 + sand, (v) Fe0 + MnO2, and (vi) Fe0 + sand + MnO2. Additional experiments were conducted to test the influence of Fe0 type, Fe0 and sand mass loadings on dye discoloration and phosphate removal in Fe0/sand system. Each system was characterized by: (i) pH value, (ii) Fe concentration, (iii) dye discoloration (MB, O-II), and (iv) phosphate concentration. Results showed that the MB method was capable of tracing the extent of iron corrosion in the various investigated systems and clarified the role of in-situ generated iron corrosion products (FeCPs) in the removal of phosphate. The suitability of mixing sand aggregates with Fe0 for efficient phosphate removal is demonstrated through the MB method. Overall, the appropriateness of the MB method to characterize the dynamics of Fe0/H2O systems is validated.