A critical look at the calculation of the binding characteristics and concentration of iron complexing ligands in seawater with suggested improvements

Environmental context The low concentration of iron in the oceans limits growth of phytoplankton. Dissolved organic molecules, called ligands, naturally present in seawater, bind iron thereby increasing its solubility and, consequently, its availability for biological uptake by phytoplankton. The characteristics of these ligands are determined indirectly with various mathematical solutions; we critically evaluate the underlying method and calculations used in these determinations. The determination of the thermodynamic characteristics of organic Fe binding ligands, total ligand concentration ([L-t]) and conditional binding constant (K '), by means of titration of natural ligands with Fe in the presence of an added known competing ligand, is an indirect method. The analysis of the titration data including the determination of the sensitivity (S) and underlying model of ligand exchange is discussed and subjected to a critical evaluation of its underlying assumptions. Large datasets collected during the International Polar Year, were used to quantify the error propagation along the determination procedure. A new and easy to handle non-linear model written in R to calculate the ligand characteristics is used. The quality of the results strongly depends on the amount of titration points or Fe additions in a titration. At least four titration points per distinguished ligand group, together with a minimum of four titration points where the ligands are saturated, are necessary to obtain statistically reliable estimates of S, K ' and [L-t]. As a result estimating the individual concentration of two ligands, although perhaps present, might not always be justified.