The kinetics of iron colloid aggregation in estuaries have been simulated with a model estuary in which seawater is continuously pumped into a reservoir initially containing river water. Profiles of colloidal Fe concentration versus salinity produced in this apparatus closely resembled field data for actual estuaries. Synthetic Fe colloids prepared by peptising Fe(OH)₃ with humic acid and phosphate showed very similar kinetic behaviour. Aggregation rate was found to be almost independent of velocity shear rate, implying that most aggregations are induced by brownian interparticle collisions. A heterogeneous kinetic model is proposed to explain the kinetic behaviour of Fe colloids during seawater-induced aggregation. This model describes kinetic behaviour in terms of a log-normal distribution of rate constants characterised by a mean value 〈k〉 and a standard deviation γ. Experiments showed that 〈k〉 is linearly related to the rate constant for salinity increase during mixing. This coupling of salinity changes and aggregation rate leads to Fe-salinity profiles that are nearly independent of the rate of salinity change, but which are dependent on γ.