Abstract. A one-dimensional model of the biogeochemistry and speciation of iron is coupled with the General Ocean Turbulence Model (GOTM) and a NPZD-type ecosystem model. The model is able to simulate the temporal patterns and vertical profiles of dissolved iron (dFe) in the upper ocean at the Bermuda Atlantic Time-series Study site reasonably well. Subsurface model profiles strongly depend on the parameter values chosen for the loss processes for iron, colloidal aggregation and scavenging onto particles. Estimates for these parameters based on observations in particle-rich waters result in depletion of dFe. A high stability constant of iron-binding organic ligands is required to reproduce the observed degree of organic complexation below the mixed layer. The scavenging residence time for iron in the model is shortest in spring and summer, because of the larger abundance of particles, and increases with depth towards values on the order of a hundred years. A solubility of atmospherically deposited iron higher than 2% lead to dFe concentrations incompatible with observations. Despite neglecting ultraviolet radiation, the model produces diurnal variations and mean vertical profiles of H₂O₂ and iron species that are in good agreement with observations.