Abstract. The transport mechanisms of floating marine debris in coastal zones remain poorly understood, primarily due to the complex geometry and influence of coastal processes which pose difficulties in incorporating them into Lagrangian numerical models at coastal scales. The numerical model presented in this study addresses these challenges using nested grids of varying hydrodynamic resolutions (2.5 km, 350 m, and 70 m). The model couples Eulerian hydrodynamic data with a Lagrangian particle solver to accurately simulate the motion of plastic particles in coastal regions. A Lagrangian validation of the model was conducted using drifter data to assess the model’s skill and confidence in simulation predictions. The results demonstrated high skill and stability in 72-hour forecast horizons and high skill score values close to the shoreline. Additionally, a beaching sensitivity analysis was performed to determine particle beaching parameterisations suitable for coastal zones. This analysis demonstrated that the real-time distance of particles to the shore during simulations was the most accurate method for detecting the land-water boundary at coastal scales, accurately representing residence times and spatiotemporal beaching patterns. Simulations were conducted using nested grids and a single coarse grid, both using the same plastic input from river outflows in a pilot test at the Barcelona coastline with a complex morphology due to coastal structures, two major river mouths, and a large harbour. The simulations revealed marked differences in beaching amounts and residence times. Both simulations exhibited high levels of beaching, with the nested grid simulation registering 91.5 % beaching and the low-resolution grid simulation showing 95.8 % beaching, surpassing other studies conducted at larger scales. Beaching amounts exhibited high levels of variability within demarcated areas between simulations, with the highest flux observed near the Llobregat River mouth release point. The inner port area showed that particles have 18 times longer residence times when using high-resolution data, demonstrating the model’s ability to solve complex coastal geometrical structures. By using nested grids, the model can resolve coastal processes with high-resolution hydrodynamic data and accurately simulate potential accumulation zones and litter hot-spots at coastal scales.