Neodymium (Nd) isotopic composition (epsilon(Nd)) is an important tracer for water mass mixing and the reconstruction of past ocean circulation. To allow for a direct model-data comparison, we have implemented Nd isotopes in the ocean component of the Community Earth System Model (CESM1.3). The model is able to capture the major features of the observed modern distribution of both epsilon(Nd) and Nd concentrations. Our model provides a useful tool for the interpretation of epsilon(Nd) reconstructions. For example, we show that in an idealized North Atlantic freshwater hosing experiment, epsilon(Nd) changes in the Atlantic are documenting primarily the changes in water mass mixing and are hardly affected by the concomitant and large changes in the marine biological productivity and organic matter fluxes. However, the hosing experiment also shows that the end-member changes due to the change of ocean circulation can influence the interpretation of epsilon(Nd) in the Atlantic, depending on the location. The implementation of Nd, together with other existing tracers, such as delta O-18, Pa-231/Th-230, delta C-13, and radiocarbon in the same model, can improve our understanding of past ocean circulation significantly.