Modeling organic iron-binding ligands in a three-dimensional biogeochemical ocean model

Modeling ligands prognostically (as opposed to assuming a uniform ligand concentration) leads to a more nutrient-like profile of iron that is more in accordance with data. It however, also leads to higher surface concentrations of dissolved iron and negative excess ligand L⁎ in some ocean regions. This is probably an indication that with more realistic and higher ligand concentrations near the surface, as opposed to the traditionally chosen low uniform concentration, iron modelers will have to re-evaluate their assumption of low scavenging rates for iron. Given their sensitivity to environmental conditions, spatio-temporal variations in ligand concentrations have the potential to impact primary production via changes in iron limitation.