Abstract. The neodymium (Nd) isotope composition (εNd) of seawater can be used to trace large-scale ocean circulation features. Yet, due to the elusive nature of marine Nd cycling, particularly in discerning non-conservative particle-seawater interactions, there remains considerable uncertainty surrounding a complete description of marine Nd budgets. Here, we present an optimisation of the Nd isotope scheme within the fast coupled atmosphere-ocean general circulation model (FAMOUS), using a statistical emulator to explore the parametric uncertainty and optimal combinations of three key model inputs relating to: (1) the efficiency of reversible scavenging, (2) the magnitude of the seafloor benthic flux, and (3) a riverine source scaling, accounting for release of Nd from river sourced particulate material. Furthermore, a suite of sensitivity tests provide insight on the regional mobilisation and spatial extent (i.e., testing a margin-constrained versus a seafloor-wide benthic flux) of certain reactive sediment components. In the calibrated scheme, the global marine Nd inventory totals 4.27 × 1012 g and has a mean residence time of 727 years. Atlantic Nd isotope distributions are represented well, and the weak sensitivity of North Atlantic Deep Water to highly unradiogenic sedimentary sources implies an abyssal benthic flux is of secondary importance in determining the water mass εNd properties under the modern vigorous circulation condition. On the other hand, Nd isotope distributions in the North Pacific are 3 to 4 εNd-units too unradiogenic compared to water measurements, and our simulations indicate that a spatially uniform flux of bulk sediment εNd does not sufficiently capture the mobile sediment components interacting with seawater. Our results of sensitivity tests suggest that there are distinct regional differences in how modern seawater acquires its εNd signal, in part relating to the complex interplay of Nd addition and water advection.