Tracking Marine Debris in Northwest Spain: Assessing Wind Influence with a Lagrangian Transport Model

Rial-Osorio, Martiño; Pérez-Muñuzuri, Vicente; Cloux, Sara

doi:10.5194/egusphere-2025-2057

Preprints

https://doi.org/10.5194/egusphere-2025-2057

Preprints

15 May 2025

| 15 May 2025

Tracking Marine Debris in Northwest Spain: Assessing Wind Influence with a Lagrangian Transport Model

Martiño Rial-Osorio, Vicente Pérez-Muñuzuri, and Sara Cloux

Abstract. Marine debris is responsible for major problems in our oceans, causing serious environmental degradation, detrimental health effects and economic losses in sectors related to the marine environment. In this study we investigate the influence of wind forcing on the transport, accumulation and beaching of floating particles in the Ría de Arousa, an estuary on the northwest coast of the Iberian Peninsula. Using Lagrangian simulations of particle tracking under different wind drag coefficients (1 %, 3 % and 5 %), we evaluate the spatial and seasonal patterns of particle concentration, residence time and deposition on the coast. Our results show that wind plays a crucial role in modulating particle behavior. Low wind-driven conditions favor greater near-shore accumulation and longer residence times, especially in the northern and inner regions of the estuary. As wind influence increases, particle dispersion intensifies, leading to lower overall accumulation and weakening of correlations between river discharge and coastal deposition. Seasonal differences are also studied, with higher concentrations observed in the north during winter and in the south during summer.

Received: 02 May 2025 – Discussion started: 15 May 2025

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.

Download & links

Preprint (PDF, 2621 KB)

Notice on discussion status
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
Preprint (2621 KB)

Download & links

The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.

Journal article(s) based on this preprint

06 Mar 2026

Tracking marine debris in Northwest Spain: assessing wind influence with a Lagrangian transport model

Martiño Rial-Osorio, Vicente Pérez-Muñuzuri, and Sara Cloux

Ocean Sci., 22, 777–790, https://doi.org/10.5194/os-22-777-2026,https://doi.org/10.5194/os-22-777-2026, 2026

Short summary

Martiño Rial-Osorio, Vicente Pérez-Muñuzuri, and Sara Cloux

Interactive discussion

Status: closed

RC1:
'Comment on egusphere-2025-2057', Anonymous Referee #1, 09 Jul 2025

This paper investigates how the distribution of floating plastic debris in the Ria de Arausa that originates from Ulla River is affected by the windage term used in Lagrangian simulations. Using this model, they find that a lower wind drag results in higher nearshore accumulation and larger retention rates. Furthermore, they find that the beaching pattern is highly heterogenous, where over half of the beached particles stranded in only 10 of the 242 coastal segments of the Ria. The number of beached particles at these segments is correlated with the river flow rate, where they find highest correlation for the sections closest to the release point of the particles.
I think in general this set up of this study allows for interesting analysis, however I feel the analysis presented here is not complete. In addition, the way model settings are chosen and why there are valid is not sufficiently explained currently. There is no discussion on the limitations of the model used here and what this implies for the obtained results. Due to these points, I think mayor revisions are needed before this work is published. I provide the following comments and questions about the manuscript, which hopefully the authors will find useful.
Comments/questions per section:
Numerical model:
Fig 1: How were the inner/intermediate and outer region chosen? From the plot, it seems that they are not equally spaced relative to the river mouth.
Fig 2: Why only show the wind roses of DJF and JJA? In the result section also the periods MAM and SON are considered and based on the results of Fig 6 it seems that the SON period has interesting wind characteristics. In addition, it would be good to not only show the wind but also the (average) Eulerian currents, as the particle distribution is strongly driven by this current. This would help with understanding the results.
Is a homogenous beaching of 80% everywhere realistic? One would expect that the beaching is a function of the specific coastal features. I think the choices that went into setting this beaching parameter have to be part of the paper.
The effect of stokes drift is not separately included. Is the assumption here that this is captured by windage? How big is the effect of swell in this region? I feel this has to be discussed in the paper.
Why were windage factors of 1%, 3% and 5% chosen, and how realistic is this range compared to the windage of real buoyant plastics? I think the work would benefit from comparing these values by typical windage factors measured in experimental settings (for example for surface drifters) to show the relevance of this range.
Why was the river input represented as a point quite far into the Ria, and not further back where the actual river mouth is located? Why were the particles not released over a line that spans the entire width of the estuary at the release point? With the current release point, there are coastal segments located behind the river input location, including segments 6 and 10. Which raises whether beaching of particles measured at these segments is realistic for plastic that originates from the river.
I think the section on how particles were simulated should contain more information. Were particles simulated for a specific time (1 season/1 year) or were they simulated until they were beached or left the domain? And if particles were simulated for the entire 5 years, how does this affect the analysis of the accumulation pattern done? As in the 5^th year, the concentration of particles then will have increased.
Results:
Fig 3: the way the data is plotted makes it hard to compare, maybe a different color bar with more shades (like magma or viridis) and/or using a discrete colormap can help to see differences. In addition, it might be useful to plot the relative concentration change for the 5% relative to the 1% instead. It is not clear whether the concentration shown here is for particles that might be in a beach cell but are yet definitely beached? Or are particles that are permanently beached also considered?
You mention that the concentration of the northern part of the estuary is higher due to the Coriolis force. I think it would be better to say due to the currents/ or explain how the Coriolis force affects the currents in the estuary. For this it would be helpful to plot the Eulerian currents (and the seasonal variations therein if present) as well (see comment for fig 2). For plastics there can also be a direct Coriolis force working on the particle because of their inertia (see the Coriolis force in the Maxey-Riley equation for example on the work of Beron-Verra et al. 2019) thus the wording Coriolis force might be confusing.
How are the northern and southern regions defined? Are they right, resp. left compared as viewed from the river mouth in direction of the Ria? And is the Island then part of the northern or southern region? Or is it an actual north/south distribution, which seems a bit arbitrary as it does not align with the shape and orientation of the Ria? In general, the south has more area located further from the river mouth, thus that would automatically result in less beaching as particles already beached close to the river mouth cannot beach there.
Fig 4 (and the discussion thereof in line 120-128) The discussion states that “during warmer months (MAM and JJA) beaching along the southern coast increases”. However, the absolute value of particles beached decreases compared to the DFJ period as well for the south? Maybe of the total particles beached, the relative amount of particle beached in the south increases, but this is not a relevant quantity to highlight? Also, how does the number of beached particles relate to the total input for each period? I think this relative number would give more insight on if particles are more likely to beach in summer or winter or that there are just more particles released into the system in winter but beaching stays constant.
Fig 5 and table 1 and discussion on the correlation: Is the correlation the best quantity to measure? If you put a lot of particles in the system, you expect that a lot of particles beach a period T after you put them into the system. This period T will be a function of the total time a particle needs to reach the area. Thus, I do not think figure 5 is the most interesting figure to show. Also in table 1, it is expected that the time delay grows with distance and thus the direct correlation will become less. In addition, the regions selected are not suitable to reveal any non-trivial spatial pattern (especially as 6 and 10 are sort of located at/in the river mouth for the location of particle release here). I think a more useful quantity to study is the cross correlation and then measuring the time lag. Or since you do Lagrangian analysis you can calculate the age of the particle (i.e. how long ago they were released) at the point in time that they reach land/beach. This would be an analysis that can be done for all coastal segments, revealing if there are any patterns different from what you expect (further away from the river mouth the particles are older). You could consider normalizing the age with the distance from the coastal segment to the particle release location.
Figure 6. It seems that the inner, intermediate and outer regions are not equally spaced and/or of the same size. This would affect the results. I would expect that for the outer region there are particles that do not reenter the Ria and thus have a residence time of infinity. Or are particles removed from the simulation if they leave the simulation domain? Which is not similarly shaped to the inner and intermediate region set here. Are beached particles considered in the calculation of the residence time?
Conclusions:
You state that wind also disperses particles more widely, but you do not show/highlight this in your result section.
Maybe highlight that the release location strongly affects the patterns of beaching in the Ria as your results are very different from Cloux et al. 2022 fig 4.
There should be a discussion on the limitations of the model used here and how these might affect the results presented here.

Citation: https://doi.org/10.5194/egusphere-2025-2057-RC1
- CC1:
  'Reply on RC1', Martiño Rial-Osorio, 05 Sep 2025
  
  Thank you very much for your comments. Your contributions have helped to improve the paper. Please find attached a PDF document with our detailed responses to each of them.
  
  Citation: https://doi.org/10.5194/egusphere-2025-2057-CC1
  - AC1: 'Reply on RC2', Martiño Rial-Osorio, 13 Oct 2025
    
    Thank you for your comments.
    Answers can be found in the attached document.
    
    Citation: https://doi.org/10.5194/egusphere-2025-2057-AC1
- AC2: 'Final reply on RC1', Martiño Rial-Osorio, 27 Oct 2025
  
  I accidentally answered as CC, and then pasted the answer to RC2 as a response to your comment. I am subsequently pasting the same comment I initially posted as CC, which contains the intended answer to your comment:
  Thank you very much for your comments. Your contributions have helped to improve the paper. Please find attached a PDF document with our detailed responses to each of them.
  
  Citation: https://doi.org/10.5194/egusphere-2025-2057-AC2
RC2:
'Comment on egusphere-2025-2057', Anonymous Referee #2, 02 Oct 2025

Preprint: https://doi.org/10.5194/egusphere-2025-2057

Review:
Tracking Marine Debris in Northwest Spain: Assessing Wind Influence with a Lagrangian Transport Model
Martiño Rial-Osorio,Vicente Pérez-Muñuzuri, and Sara Cloux
The article addresses a topic of considerable environmental and scientific relevance: the transport and accumulation of plastics in estuarine systems under wind influence. The study combines Lagrangian simulations with varying windage coefficients to evaluate seasonal accumulation patterns, residence times, and correlations with river discharge in the Ría de Arousa (NW Spain). The work is well-structured, and the results are consistent with previous literature, providing valuable insights into the processes governing the retention of marine debris in estuarine environments.
Nonetheless, before publication, several aspects require clarification and further detail, particularly with respect to methodological details and contextual framing. Overall, I recommend the manuscript for publication, subject to minor revisions.

Some points that I believe require clarification or improvement:
1. Introduction
In the manuscript, the distinction between macro- and microplastics is presented in centimetres. However, the standard unit, used in the literature, is typically millimetres, as indicated in the references cited by the authors themselves (Min et al., 2020; Wayman & Niemann, 2021). I recommend aligning this definition with the units in the literature, to avoid terminological inconsistencies.
2.1 Study Area
The criteria used to divide the study area into the sections shown in Figure 1 are not clearly explained. Does this segmentation arise from hydrodynamic differences, socio-economic characteristics, or another specific parameter? The authors should clarify the rationale behind this segmentation to justify its relevance.
While I understand the choice to present simulations for winter (DJF) and summer (JJA), as these represent extreme conditions, it would be advisable to also include wind roses for the intermediate seasons (MAM and SON) in Figure 2, in order to provide a more comprehensive view of the seasonal variability of the wind patterns analysed in the study. Additionally, since the results also refer to these intermediate seasons, their inclusion would be pertinent.
Considering that the entire study is based on particle discharge proportional to the flow of the Ulla River, it would be important to provide, in this section, information on the average, minimum, and maximum river discharge values.
2.2 Input Data
Although the Numerical Methods section specifies that simulations were performed using 2D surface fields, I consider that this information should also appear in the Input Data section. This would prevent the reader, upon first encountering the description of the hydrodynamic model, from assuming that vertical processes or particle sinking were also considered. A simple statement indicating that only 2D surface fields were used would make the methodology clearer from the start.
2.3 Numerical Methods
The adoption of an 80% beaching probability is a central aspect of the study, yet the manuscript does not provide the rationale for this choice. It would strengthen the work if the authors explained this choice in more detail. I understand that a high retention rate may be associated with rocky coasts and intertidal zones, where stranded plastics tend to remain longer due to lower remobilization. However, in the correlation section (Table 1), the analysis focuses mainly on sandy beach segments, where higher remobilization would be expected. Therefore, it would be useful to clarify this apparent discrepancy and discuss the choice of the 80% value for the different types of coast included in the study.
A key element of the work is the analysis using different windage values (cw = 1%, 3%, and 5%). However, this choice could be better supported with additional references or, ideally, observational data. Providing such justification would increase the robustness and relevance of the adopted parameterization (e.g., https://doi.org/10.1016/j.envint.2020.105655 and https://doi.org/10.1038/s41598-018-22939-w).
Results and Conclusion:
Figure 3 is not entirely informative. The current colour palette limits the ability to discern differences between areas of higher particle accumulation, which reduces the figure’s interpretative value.
The authors state that particles tend to concentrate upstream due to the Coriolis force. While this is correct, it would be helpful to explicitly relate it to the resulting residual circulation. Including a figure illustrating the residual circulation would significantly enhance the discussion and provide clearer mechanistic insight.
There are also concerns regarding the discussion of the results, which primarily focuses on a single winter and a single summer scenario with different windage values. Given that, particle release is substantially higher in winter, Figure 5 appears largely redundant. If the goal is to assess whether beaching is more frequent in summer or winter and under varying windage conditions, it would be more scientifically robust to compare scenarios with similar particle release rates across seasons. Alternatively, considering particle age—the time from release to beaching—across all scenarios would provide a more meaningful and grounded comparison.
Addressing these points would greatly strengthen the manuscript by clarifying the drivers of particle accumulation and providing a more rigorous basis for seasonal and windage-related comparisons.
The definition and discussion of residence time appears to be somewhat unclear. From my perspective, this parameter seems strongly influenced by the size of the areas considered and the number of particles released. In principle, a larger number of particles or a more extensive area would tend to increase residence time, but this relationship is not necessarily true from a hydrodynamic standpoint.
Moreover, the calculations also appear to depend on the initial time of the simulation and the number of particles released. For example, if the simulation begins in January, the residence time estimated for this first month will inevitably be shorter than in March—not due to hydrodynamic conditions or wind, but simply because fewer particles are present to be account for.
The interpretation of Figure 6 is also not entirely clear. In the innermost zone of the estuary, residence time may be higher due to proximity to the source. In the intermediate zone, however, since particles that return to this area are counted again, the probability of this occurring is naturally higher, particularly, given the way the estuary has been subdivided. This hypothesis could be further supported by including a figure showing the Eulerian residual circulation of velocities, which would provide better context for the results.
Additionally, according to the methodological criterion that particles leaving the study area are considered “dead,” the downstream region is excluded from the residence time calculation, which limits interpretation for this outer area.
Therefore, it would be important to clarify how residence time was computed or, alternatively, to explicitly highlight these limitations in the discussion to avoid potentially biased conclusions.
I would encourage the authors to explicitly discuss the limitations of the implemented methodology. For instance, it is not clear whether a horizontal resolution of 300 m is sufficient, considering the dimensions of the zones chosen as accumulation areas. Similarly, the potential effects of three-dimensional processes should be addressed: to what extent might a 3D simulation alter the results compared to the current approach?
The role of waves in particle distribution within the estuary, as well as in the beaching process, also warrants consideration, since they may significantly influence particle trajectories and residence times.
Finally, it would be valuable to discuss the generalizability of the methodology: how could it be applied to other study areas, and are there comparable studies in different regions that support or challenge the current approach? Addressing these points would strengthen the manuscript by providing a clearer context and by acknowledging the methodological constraints.

Citation: https://doi.org/10.5194/egusphere-2025-2057-RC2
- AC1: 'Reply on RC2', Martiño Rial-Osorio, 13 Oct 2025
  
  Thank you for your comments.
  Answers can be found in the attached document.
  
  Citation: https://doi.org/10.5194/egusphere-2025-2057-AC1

Interactive discussion

Status: closed

RC1:
'Comment on egusphere-2025-2057', Anonymous Referee #1, 09 Jul 2025

This paper investigates how the distribution of floating plastic debris in the Ria de Arausa that originates from Ulla River is affected by the windage term used in Lagrangian simulations. Using this model, they find that a lower wind drag results in higher nearshore accumulation and larger retention rates. Furthermore, they find that the beaching pattern is highly heterogenous, where over half of the beached particles stranded in only 10 of the 242 coastal segments of the Ria. The number of beached particles at these segments is correlated with the river flow rate, where they find highest correlation for the sections closest to the release point of the particles.
I think in general this set up of this study allows for interesting analysis, however I feel the analysis presented here is not complete. In addition, the way model settings are chosen and why there are valid is not sufficiently explained currently. There is no discussion on the limitations of the model used here and what this implies for the obtained results. Due to these points, I think mayor revisions are needed before this work is published. I provide the following comments and questions about the manuscript, which hopefully the authors will find useful.
Comments/questions per section:
Numerical model:
Fig 1: How were the inner/intermediate and outer region chosen? From the plot, it seems that they are not equally spaced relative to the river mouth.
Fig 2: Why only show the wind roses of DJF and JJA? In the result section also the periods MAM and SON are considered and based on the results of Fig 6 it seems that the SON period has interesting wind characteristics. In addition, it would be good to not only show the wind but also the (average) Eulerian currents, as the particle distribution is strongly driven by this current. This would help with understanding the results.
Is a homogenous beaching of 80% everywhere realistic? One would expect that the beaching is a function of the specific coastal features. I think the choices that went into setting this beaching parameter have to be part of the paper.
The effect of stokes drift is not separately included. Is the assumption here that this is captured by windage? How big is the effect of swell in this region? I feel this has to be discussed in the paper.
Why were windage factors of 1%, 3% and 5% chosen, and how realistic is this range compared to the windage of real buoyant plastics? I think the work would benefit from comparing these values by typical windage factors measured in experimental settings (for example for surface drifters) to show the relevance of this range.
Why was the river input represented as a point quite far into the Ria, and not further back where the actual river mouth is located? Why were the particles not released over a line that spans the entire width of the estuary at the release point? With the current release point, there are coastal segments located behind the river input location, including segments 6 and 10. Which raises whether beaching of particles measured at these segments is realistic for plastic that originates from the river.
I think the section on how particles were simulated should contain more information. Were particles simulated for a specific time (1 season/1 year) or were they simulated until they were beached or left the domain? And if particles were simulated for the entire 5 years, how does this affect the analysis of the accumulation pattern done? As in the 5^th year, the concentration of particles then will have increased.
Results:
Fig 3: the way the data is plotted makes it hard to compare, maybe a different color bar with more shades (like magma or viridis) and/or using a discrete colormap can help to see differences. In addition, it might be useful to plot the relative concentration change for the 5% relative to the 1% instead. It is not clear whether the concentration shown here is for particles that might be in a beach cell but are yet definitely beached? Or are particles that are permanently beached also considered?
You mention that the concentration of the northern part of the estuary is higher due to the Coriolis force. I think it would be better to say due to the currents/ or explain how the Coriolis force affects the currents in the estuary. For this it would be helpful to plot the Eulerian currents (and the seasonal variations therein if present) as well (see comment for fig 2). For plastics there can also be a direct Coriolis force working on the particle because of their inertia (see the Coriolis force in the Maxey-Riley equation for example on the work of Beron-Verra et al. 2019) thus the wording Coriolis force might be confusing.
How are the northern and southern regions defined? Are they right, resp. left compared as viewed from the river mouth in direction of the Ria? And is the Island then part of the northern or southern region? Or is it an actual north/south distribution, which seems a bit arbitrary as it does not align with the shape and orientation of the Ria? In general, the south has more area located further from the river mouth, thus that would automatically result in less beaching as particles already beached close to the river mouth cannot beach there.
Fig 4 (and the discussion thereof in line 120-128) The discussion states that “during warmer months (MAM and JJA) beaching along the southern coast increases”. However, the absolute value of particles beached decreases compared to the DFJ period as well for the south? Maybe of the total particles beached, the relative amount of particle beached in the south increases, but this is not a relevant quantity to highlight? Also, how does the number of beached particles relate to the total input for each period? I think this relative number would give more insight on if particles are more likely to beach in summer or winter or that there are just more particles released into the system in winter but beaching stays constant.
Fig 5 and table 1 and discussion on the correlation: Is the correlation the best quantity to measure? If you put a lot of particles in the system, you expect that a lot of particles beach a period T after you put them into the system. This period T will be a function of the total time a particle needs to reach the area. Thus, I do not think figure 5 is the most interesting figure to show. Also in table 1, it is expected that the time delay grows with distance and thus the direct correlation will become less. In addition, the regions selected are not suitable to reveal any non-trivial spatial pattern (especially as 6 and 10 are sort of located at/in the river mouth for the location of particle release here). I think a more useful quantity to study is the cross correlation and then measuring the time lag. Or since you do Lagrangian analysis you can calculate the age of the particle (i.e. how long ago they were released) at the point in time that they reach land/beach. This would be an analysis that can be done for all coastal segments, revealing if there are any patterns different from what you expect (further away from the river mouth the particles are older). You could consider normalizing the age with the distance from the coastal segment to the particle release location.
Figure 6. It seems that the inner, intermediate and outer regions are not equally spaced and/or of the same size. This would affect the results. I would expect that for the outer region there are particles that do not reenter the Ria and thus have a residence time of infinity. Or are particles removed from the simulation if they leave the simulation domain? Which is not similarly shaped to the inner and intermediate region set here. Are beached particles considered in the calculation of the residence time?
Conclusions:
You state that wind also disperses particles more widely, but you do not show/highlight this in your result section.
Maybe highlight that the release location strongly affects the patterns of beaching in the Ria as your results are very different from Cloux et al. 2022 fig 4.
There should be a discussion on the limitations of the model used here and how these might affect the results presented here.

Citation: https://doi.org/10.5194/egusphere-2025-2057-RC1
- CC1:
  'Reply on RC1', Martiño Rial-Osorio, 05 Sep 2025
  
  Thank you very much for your comments. Your contributions have helped to improve the paper. Please find attached a PDF document with our detailed responses to each of them.
  
  Citation: https://doi.org/10.5194/egusphere-2025-2057-CC1
  - AC1: 'Reply on RC2', Martiño Rial-Osorio, 13 Oct 2025
    
    Thank you for your comments.
    Answers can be found in the attached document.
    
    Citation: https://doi.org/10.5194/egusphere-2025-2057-AC1
- AC2: 'Final reply on RC1', Martiño Rial-Osorio, 27 Oct 2025
  
  I accidentally answered as CC, and then pasted the answer to RC2 as a response to your comment. I am subsequently pasting the same comment I initially posted as CC, which contains the intended answer to your comment:
  Thank you very much for your comments. Your contributions have helped to improve the paper. Please find attached a PDF document with our detailed responses to each of them.
  
  Citation: https://doi.org/10.5194/egusphere-2025-2057-AC2
RC2:
'Comment on egusphere-2025-2057', Anonymous Referee #2, 02 Oct 2025

Preprint: https://doi.org/10.5194/egusphere-2025-2057

Review:
Tracking Marine Debris in Northwest Spain: Assessing Wind Influence with a Lagrangian Transport Model
Martiño Rial-Osorio,Vicente Pérez-Muñuzuri, and Sara Cloux
The article addresses a topic of considerable environmental and scientific relevance: the transport and accumulation of plastics in estuarine systems under wind influence. The study combines Lagrangian simulations with varying windage coefficients to evaluate seasonal accumulation patterns, residence times, and correlations with river discharge in the Ría de Arousa (NW Spain). The work is well-structured, and the results are consistent with previous literature, providing valuable insights into the processes governing the retention of marine debris in estuarine environments.
Nonetheless, before publication, several aspects require clarification and further detail, particularly with respect to methodological details and contextual framing. Overall, I recommend the manuscript for publication, subject to minor revisions.

Some points that I believe require clarification or improvement:
1. Introduction
In the manuscript, the distinction between macro- and microplastics is presented in centimetres. However, the standard unit, used in the literature, is typically millimetres, as indicated in the references cited by the authors themselves (Min et al., 2020; Wayman & Niemann, 2021). I recommend aligning this definition with the units in the literature, to avoid terminological inconsistencies.
2.1 Study Area
The criteria used to divide the study area into the sections shown in Figure 1 are not clearly explained. Does this segmentation arise from hydrodynamic differences, socio-economic characteristics, or another specific parameter? The authors should clarify the rationale behind this segmentation to justify its relevance.
While I understand the choice to present simulations for winter (DJF) and summer (JJA), as these represent extreme conditions, it would be advisable to also include wind roses for the intermediate seasons (MAM and SON) in Figure 2, in order to provide a more comprehensive view of the seasonal variability of the wind patterns analysed in the study. Additionally, since the results also refer to these intermediate seasons, their inclusion would be pertinent.
Considering that the entire study is based on particle discharge proportional to the flow of the Ulla River, it would be important to provide, in this section, information on the average, minimum, and maximum river discharge values.
2.2 Input Data
Although the Numerical Methods section specifies that simulations were performed using 2D surface fields, I consider that this information should also appear in the Input Data section. This would prevent the reader, upon first encountering the description of the hydrodynamic model, from assuming that vertical processes or particle sinking were also considered. A simple statement indicating that only 2D surface fields were used would make the methodology clearer from the start.
2.3 Numerical Methods
The adoption of an 80% beaching probability is a central aspect of the study, yet the manuscript does not provide the rationale for this choice. It would strengthen the work if the authors explained this choice in more detail. I understand that a high retention rate may be associated with rocky coasts and intertidal zones, where stranded plastics tend to remain longer due to lower remobilization. However, in the correlation section (Table 1), the analysis focuses mainly on sandy beach segments, where higher remobilization would be expected. Therefore, it would be useful to clarify this apparent discrepancy and discuss the choice of the 80% value for the different types of coast included in the study.
A key element of the work is the analysis using different windage values (cw = 1%, 3%, and 5%). However, this choice could be better supported with additional references or, ideally, observational data. Providing such justification would increase the robustness and relevance of the adopted parameterization (e.g., https://doi.org/10.1016/j.envint.2020.105655 and https://doi.org/10.1038/s41598-018-22939-w).
Results and Conclusion:
Figure 3 is not entirely informative. The current colour palette limits the ability to discern differences between areas of higher particle accumulation, which reduces the figure’s interpretative value.
The authors state that particles tend to concentrate upstream due to the Coriolis force. While this is correct, it would be helpful to explicitly relate it to the resulting residual circulation. Including a figure illustrating the residual circulation would significantly enhance the discussion and provide clearer mechanistic insight.
There are also concerns regarding the discussion of the results, which primarily focuses on a single winter and a single summer scenario with different windage values. Given that, particle release is substantially higher in winter, Figure 5 appears largely redundant. If the goal is to assess whether beaching is more frequent in summer or winter and under varying windage conditions, it would be more scientifically robust to compare scenarios with similar particle release rates across seasons. Alternatively, considering particle age—the time from release to beaching—across all scenarios would provide a more meaningful and grounded comparison.
Addressing these points would greatly strengthen the manuscript by clarifying the drivers of particle accumulation and providing a more rigorous basis for seasonal and windage-related comparisons.
The definition and discussion of residence time appears to be somewhat unclear. From my perspective, this parameter seems strongly influenced by the size of the areas considered and the number of particles released. In principle, a larger number of particles or a more extensive area would tend to increase residence time, but this relationship is not necessarily true from a hydrodynamic standpoint.
Moreover, the calculations also appear to depend on the initial time of the simulation and the number of particles released. For example, if the simulation begins in January, the residence time estimated for this first month will inevitably be shorter than in March—not due to hydrodynamic conditions or wind, but simply because fewer particles are present to be account for.
The interpretation of Figure 6 is also not entirely clear. In the innermost zone of the estuary, residence time may be higher due to proximity to the source. In the intermediate zone, however, since particles that return to this area are counted again, the probability of this occurring is naturally higher, particularly, given the way the estuary has been subdivided. This hypothesis could be further supported by including a figure showing the Eulerian residual circulation of velocities, which would provide better context for the results.
Additionally, according to the methodological criterion that particles leaving the study area are considered “dead,” the downstream region is excluded from the residence time calculation, which limits interpretation for this outer area.
Therefore, it would be important to clarify how residence time was computed or, alternatively, to explicitly highlight these limitations in the discussion to avoid potentially biased conclusions.
I would encourage the authors to explicitly discuss the limitations of the implemented methodology. For instance, it is not clear whether a horizontal resolution of 300 m is sufficient, considering the dimensions of the zones chosen as accumulation areas. Similarly, the potential effects of three-dimensional processes should be addressed: to what extent might a 3D simulation alter the results compared to the current approach?
The role of waves in particle distribution within the estuary, as well as in the beaching process, also warrants consideration, since they may significantly influence particle trajectories and residence times.
Finally, it would be valuable to discuss the generalizability of the methodology: how could it be applied to other study areas, and are there comparable studies in different regions that support or challenge the current approach? Addressing these points would strengthen the manuscript by providing a clearer context and by acknowledging the methodological constraints.

Citation: https://doi.org/10.5194/egusphere-2025-2057-RC2
- AC1: 'Reply on RC2', Martiño Rial-Osorio, 13 Oct 2025
  
  Thank you for your comments.
  Answers can be found in the attached document.
  
  Citation: https://doi.org/10.5194/egusphere-2025-2057-AC1

Peer review completion

AR – Author's response | RR – Referee report | ED – Editor decision | EF – Editorial file upload

AR by Martiño Rial-Osorio on behalf of the Authors (29 Oct 2025) Author's response Author's tracked changes Manuscript

ED: Referee Nomination & Report Request started (30 Oct 2025) by Damian Leonardo Arévalo-Martínez

RR by Anonymous Referee #3 (16 Jan 2026)

Suggestions for revision or reasons for rejection

This paper investigates the effects of different windage factors on floating particle transport in the Ría de Arousa. The results could be useful to interpreting model results that incorporated windage, predicting transport differences between particles with different windage factors, and understanding particle transport in the Ría de Arousa. Overall, the authors have responded satisfactorily to the previous reviewer comments. However, additional clarification on a few points would be worthwhile, as detailed below.

1. The beaching parametrization could be further clarified.
a. What shoreline interaction threshold was used?
b. The justification for using a high beaching factor is reasonable, but it would be helpful to state the typical range of beaching factors to give readers without prior context a basis for comparison.

2. The discussion of residence time could also be further clarified.
a. My understanding is that the boundaries of the outer region go to the edge of the model domain. The revised manuscript says that particles that exit the domain are not included in the residence time calculation. If so, only particles that exit the outer region into the intermediate region are included in the residence time calculation, which could lead to a significant bias in the residence time results for the outer region.
b. From the drawing in figure 1, the water-covered area in the outer region appears to be substantially larger than in the inner region. The different areas of the regions could impact the residence times. It would be helpful for the authors to briefly comment on the impact of region size on the residence times.
c. Error bars on figure 7 would be helpful to understanding the variability in the residence times.

3. Along with the previous reviewers, I think that a brief comment on how Coriolis affects the ría with references would help readers unfamiliar with the area to understand the publication.

4. The discussion of figure 5 should remind readers that the number of particles released is correlated with the river discharge.

5. The time lag discussion is very interesting. This suggestion is very much up to the authors, but it may be helpful to put the results in a map so that readers can quickly understand where the time lags are the highest. Alternatively, the beaches could be renumbered based on proximity to the river mouth.

6. The discussion of the limitations of the study is very useful. One additional consideration that may be worth mentioning is wind mixing of microplastics, since high wind conditions would lead to strong windage effects in the model but may mix microplastics down into the water column (e.g., Kukulka et al, 2012).

7. The authors note that their findings significantly deviate from Cloux et. al. 2022 due to a difference in particle release position. It would be helpful if the authors explicitly stated how their findings differed from Cloux et. al. 2022, so that readers unfamiliar with the prior work can understand the comment.

Hide

ED: Publish subject to minor revisions (review by editor) (19 Jan 2026) by Damian Leonardo Arévalo-Martínez

AR by Martiño Rial-Osorio on behalf of the Authors (27 Jan 2026) Author's response Author's tracked changes Manuscript

ED: Publish as is (27 Jan 2026) by Damian Leonardo Arévalo-Martínez

AR by Martiño Rial-Osorio on behalf of the Authors (30 Jan 2026) Author's response Manuscript

Journal article(s) based on this preprint

06 Mar 2026

Tracking marine debris in Northwest Spain: assessing wind influence with a Lagrangian transport model

Martiño Rial-Osorio, Vicente Pérez-Muñuzuri, and Sara Cloux

Ocean Sci., 22, 777–790, https://doi.org/10.5194/os-22-777-2026,https://doi.org/10.5194/os-22-777-2026, 2026

Short summary

Martiño Rial-Osorio, Vicente Pérez-Muñuzuri, and Sara Cloux

Viewed

Total article views: 1,059 (including HTML, PDF, and XML)

HTML	PDF	XML	Total	BibTeX	EndNote
827	203	29	1,059	21	42

HTML: 827
PDF: 203
XML: 29
Total: 1,059
BibTeX: 21
EndNote: 42

Views and downloads (calculated since 15 May 2025)

Month	HTML	PDF	XML	Total
May 2025	42	17	5	64
Jun 2025	53	12	4	69
Jul 2025	46	16	3	65
Aug 2025	83	3	0	86
Sep 2025	422	5	1	428
Oct 2025	54	19	6	79
Nov 2025	18	11	2	31
Dec 2025	34	32	0	66
Jan 2026	35	42	5	82
Feb 2026	32	42	3	77
Mar 2026	8	4	0	12

Cumulative views and downloads (calculated since 15 May 2025)

Month	HTML	PDF	XML	Total
May 2025	42	17	5	64
Jun 2025	53	12	4	69
Jul 2025	46	16	3	65
Aug 2025	83	3	0	86
Sep 2025	422	5	1	428
Oct 2025	54	19	6	79
Nov 2025	18	11	2	31
Dec 2025	34	32	0	66
Jan 2026	35	42	5	82
Feb 2026	32	42	3	77
Mar 2026	8	4	0	12

Viewed (geographical distribution)

Total article views: 1,007 (including HTML, PDF, and XML) Thereof 1,007 with geography defined and 0 with unknown origin.

Country	#	Views	%

Latest update: 06 Mar 2026

Download

The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.

Preprint (2621 KB)
Metadata XML

Short summary

We used the MOHID-Lagrangian model to simulate the transport and accumulation of floating debris in the Ría de Arousa (2018–2023). Two-dimensional surface simulations with windage coefficients of 1 %, 3 %, and 5 % represented different debris types. Seasonal and spatial patterns were analyzed using a high-resolution grid and 242 shoreline segments to assess wind and seasonal effects on debris distribution and retention.


Total:	0
HTML:	0
PDF:	0
XML:	0