Hello everyone!
I recently tried implementing iron fertilization in the ecosys_forcing_mod.F90 module of the POP model, with the simulated fertilization region located in the eastern Pacific Ocean. However, I observed a surprising result: within just a few months, surface dissolved iron (Fe) concentrations changed significantly across the entire global ocean.
This seems physically implausible—given the timescales of large-scale ocean circulation, such a rapid and widespread response shouldn’t occur from a localized perturbation. The influence should be more regional and gradual.
Has anyone encountered a similar issue? Could this be related to model numerics, tracer advection schemes, boundary conditions, or perhaps an unintended global coupling in the biogeochemical module? I’d greatly appreciate any insights or suggestions to help diagnose this behavior.
Thank you in advance!
I recently tried implementing iron fertilization in the ecosys_forcing_mod.F90 module of the POP model, with the simulated fertilization region located in the eastern Pacific Ocean. However, I observed a surprising result: within just a few months, surface dissolved iron (Fe) concentrations changed significantly across the entire global ocean.
This seems physically implausible—given the timescales of large-scale ocean circulation, such a rapid and widespread response shouldn’t occur from a localized perturbation. The influence should be more regional and gradual.
Has anyone encountered a similar issue? Could this be related to model numerics, tracer advection schemes, boundary conditions, or perhaps an unintended global coupling in the biogeochemical module? I’d greatly appreciate any insights or suggestions to help diagnose this behavior.
Thank you in advance!