The global marine iron cycle has been modeled using REcoM2 (Regulated Ecosystem Model 2) coupled with the ocean circulation model MITgcm (MIT General Circulation Model). A general description of the model can be found in the supporting material for Hauck et al. 2013. In the simulations producing this data set, the organic complexation of iron was made dependent on seawater pH and concentration of dissolved organic matter. With this parameterization, the effect of pH change on the iron cycle and the biological feedback to the iron cycle have been examined through the pathway of ligand binding.In this data set, model output of dissolved iron (DFe), dissolved organic carbon (DOC), stability constants of ligands, free Fe concentration and net primary production (NPP) are included. Four different simulations are considered: R_std (standard), R_constL (with constant 1 nM ligand concentration), R_progL (with one prognostic ligand) and R_ph (with future pH field).