the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 04 Jul 2024
A wave-resolving 2DV Lagrangian approach to model microplastic transport in the nearshore
Abstract. Potentially acting as a source or a sink for plastic pollution to the open ocean, nearshore waters remain a challenging context for predicting the transport and deposition of plastic debris. In this study, we present an advanced modelling approach based on the SWASH wave model and the TrackMPD (v3.0) particle transport model to investigate the transport dynamics of floating and sinking microplastics in wave-dominated environments. This approach introduces novel features such as coupling with advanced turbulence models, simulating resuspension and bedload processes, implementing advanced settling and rising velocity formulations, and enabling parallel computation. The wave laboratory experiments conducted by Forsberg et al. (2020) were simulated to validate the model's ability to reproduce the transport of diverse microplastics (varying in density, shape, and size) along a comprehensive beach profile, capturing the whole water column. Our results underscore the robustness of the proposed model, showing good agreement with experimental data. High-density microplastics moved onshore near the bed accumulating in proximity to the wave-breaking zone, while the distribution of low-density particles varied along the coastal profile depending on the particle properties. The study also sheds light on the primary mechanisms driving microplastic transport, such as Stokes drift, wave asymmetry and settling/rising velocities. Sensitivity analyses on calibration parameters further confirm the robustness of the model results and the influence of these factors on transport patterns. This research establishes the SWASH-TrackMPD approach as a valuable tool, opening avenues for future studies to contextualize laboratory findings within the complexities of real-world nearshore environments and further refine our comprehension of microplastic dynamics across different beaches and wave-climate conditions.
- Preprint
(1731 KB) - Metadata XML
- BibTeX
- EndNote
Status: closed
-
CEC1: 'Comment on gmd-2024-100', Astrid Kerkweg, 15 Jul 2024
Dear authors,
in my role as Executive editor of GMD, I would like to bring to your attention our Editorial version 1.2: https://www.geosci-model-dev.net/12/2215/2019/
This highlights some requirements of papers published in GMD, which is also available on the GMD website in the ‘Manuscript Types’ section: http://www.geoscientific-model-development.net/submission/manuscript_types.html
In particular, please note that for your paper, the following requirements have not been met in the Discussions paper:
- "The main paper must give the model name and version number (or other unique identifier) in the title."
In order to simplify reference to your developments, please add a model name (and/or its acronym) and a version number in the title of your article in your revised submission to GMD.
Yours, Astrid Kerkweg
Citation: https://doi.org/10.5194/gmd-2024-100-CEC1 -
AC1: 'Reply on CEC1', Isabel Jalon-Rojas, 11 Sep 2024
Dear Editor,
Following the journal's requirements, we will modify the title to include the model name and version number: 'A wave-resolving 2DV Lagrangian approach to model microplastic transport in the nearshore based on TrackMPD v3.0.'"
Best regards,
Isabel
Citation: https://doi.org/10.5194/gmd-2024-100-AC1
-
RC1: 'Comment on gmd-2024-100', Anonymous Referee #1, 16 Jul 2024
This paper presents a novel approach to numerically study the transport of buoyant and non-buoyant microplastics in wave-dominated nearshore areas, a highly relevant and innovative topic. The Authors combine the SWASH numerical model to resolve surf zone hydrodynamics with the Lagrangian model TrakMPD to track microplastics. They enhanced the Lagrangian model to better represent turbulent mixing with spatial gradients from wave breaking and depth variations and included processes such as deposition, resuspension, and bedload transport.
The authors reproduced one of Forsberg's windless hydrodynamic laboratory experiments, parameterizing the model based on existing literature Under this condition, they evaluated the transport and dispersion of six types of microplastics (spheres, fibres, and sheets, both buoyant and non-buoyant) and found a good match between numerical results and laboratory findings in the final concentration distribution across different beach profile areas. The main discrepancy occurred with bottom-transported microplastics, which escaped the breaking zone into the adjacent shoaling zone. This is likely due to a slight misrepresentation of the experimental undertow current, the temporal evolution of vertical mixing, or inaccuracies in defining the compartment limits, as the authors suggest. Subsequently, the Authors conducted a sensitivity analysis of model parameters to explore the influence of turbulence from wave breaking, bedload transport, and resuspension on microplastic transport. They also investigated mechanisms driving microplastic transport, such as Stokes drift, wave asymmetry, and particle settling and rising velocities.
The topic addressed by this paper is highly relevant and potentially a valuable contribution to the Scientific Community. The paper is concise, well-organized, and well-written. To further enhance the manuscript, the following minor aspects and suggestions should be addressed:
Minor comments
L131-138: In Section 2.2.1, general aspects and model improvements are described, without addressing the specific configuration for the present study. Therefore, I would suggest relocating these lines (L131-137, “For the present simulation…. uniform throughout the domain.”) closer to line L223, where the chosen model parameters for the study are detailed.
L133: Please replace “surfzone” with “surf zone”.
L136: Here, a sentence could be added to justify the hypothesis that the vertical diffusion coefficient for microplastics equals the eddy diffusivity, given that the small size of microplastics ensures their behaviour as passive tracers predominantly governed by fluid turbulence.
L190: I would explicitly state that one of the two windless hydrodynamic conditions evaluated by Forsberg was replicated numerically.
L199: It is recommended to include a quantitative/statistical comparison between numerical and laboratory wave heights and an assessment of how much the model underestimates dissipation in the surf zone.
L207: I believe the Jalon-Rojas method (2022) estimates rising and settling velocities in calm water. De Leo et al. (2021; https://doi.org/10.3390/jmse9020142) found that settling velocities increase in the presence of waves. I would add a sentence noting this aspect.
L212 & L245: It should be clarified whether the results presented in the manuscript are the average of all runs or from a representative run.
P11, Section 3 (Results): Similarly to line L261, where it states “...54% and 67% of low-density fibres were gradually transported onshore...”, I would like to see more quantification in the description of spheres and sheets throughout this section.
Citation: https://doi.org/10.5194/gmd-2024-100-RC1 -
AC2: 'Reply on RC1', Isabel Jalon-Rojas, 11 Sep 2024
We thank the reviewer for carefully reading our manuscript, for positive comments, and useful suggestions. Below, we address each of the points and outline how we plan to revise the manuscript accordingly:
Minor comments
L131-138: In Section 2.2.1, general aspects and model improvements are described, without addressing the specific configuration for the present study. Therefore, I would suggest relocating these lines (L131-137, “For the present simulation…. uniform throughout the domain.”) closer to line L223, where the chosen model parameters for the study are detailed.
Yes, we will relocate these lines in Section 2.4 as suggested
L133: Please replace “surfzone” with “surf zone”.
Yes, we will replace it.
L136: Here, a sentence could be added to justify the hypothesis that the vertical diffusion coefficient for microplastics equals the eddy diffusivity, given that the small size of microplastics ensures their behaviour as passive tracers predominantly governed by fluid turbulence.
Good suggestion, we will include this explanation in the text as suggested
L190: I would explicitly state that one of the two windless hydrodynamic conditions evaluated by Forsberg was replicated numerically.
Yes, we will detail that we have replicated one of the two windless hydrodynamic conditions evaluated by Forsberg.
L199: It is recommended to include a quantitative/statistical comparison between numerical and laboratory wave heights and an assessment of how much the model underestimates dissipation in the surf zone.
We will modify the sentence to include a quantitative comparison: "Numerical wave heights are generally in good agreement with experimental measurements, as illustrated in Figure 1a. Specifically, the model tends to slightly underestimate surf zone dissipation, leading to a root mean square error (RMSE) of 2.2 cm and an average negative bias of 0.7 cm for wave heights in this region. Despite this minor discrepancy, the overall agreement in wave heights supports using SWASH simulations as a reliable hydrodynamic forcing for TrackMPD."
L207: I believe the Jalon-Rojas method (2022) estimates rising and settling velocities in calm water. De Leo et al. (2021; https://doi.org/10.3390/jmse9020142) found that settling velocities increase in the presence of waves. I would add a sentence noting this aspect.
De Leo et al. (2021) proposed a formulation for settling velocity to account for the influence of waves on vertical trajectories. In our study, the vertical trajectories of particles are already influenced by the wave-resolved hydrodynamics within the SWASH model. Since the effects of waves are inherently captured by the hydrodynamics, we do not need to directly adjust the settling velocity for wave-induced variations, so we use values in calm water. The adjustment proposed by Leo et al. (2021) could indeed be relevant for large-scale models that do not explicitly resolve wave dynamics but, in our case, the wave influences are accounted for through the hydrodynamic forcing.
L212 & L245: It should be clarified whether the results presented in the manuscript are the average of all runs or from a representative run.
Figures 2.b and 3.a-b illustrate the main results of the manuscript, showing the average values and standard deviations from five runs for each scenario. As noted in the text: “to ensure the robustness of our simulations, we conducted five simulations for each scenario, and the results consistently exhibited only minor variability in the number of particles within each region. This variability was within the same order of magnitude as that observed in the experiments as indicated by the error bars in Figure 2.” In Figures 2.a and 3.c, we present the trajectories and final positions from one of the runs, which, given the low variability between runs, is representative of the scenario. We will clarify this last point in the text.
P11, Section 3 (Results): Similarly to line L261, where it states “...54% and 67% of low-density fibres were gradually transported onshore...”, I would like to see more quantification in the description of spheres and sheets throughout this section.
We will include this quantification for low-density sheets as well. For low-density spheres and all high-density particles, nearly all particles ended up in a single compartment, as mentioned in the text.
Citation: https://doi.org/10.5194/gmd-2024-100-AC2
-
AC2: 'Reply on RC1', Isabel Jalon-Rojas, 11 Sep 2024
-
RC2: 'Comment on gmd-2024-100', Anonymous Referee #2, 02 Sep 2024
The manuscript describes the modifications proposed to an existing numerical modelling approach (TrackMPD) for plastic debris with the aim to improve the transport in coastal zones. The modified numerical approach was then benchmarked against experimental observations of the debris transport in the nearshore zone (breaking zone, swash zone) recently published by one of the Authors. Differently from the original formulations of TrackMPD, in the present study the hydrodynamics is solved used SWASH, a new formulation is proposed for the vertical diffusivity (linked to the wave breaking), and the computation of the transport at the bed. The topic is of interest from the scientific point of view and with possible high impact in terms of marine pollution. In fact, a reliable numerical model able to predict the plastic debris trajectories is absolutely need to contrast plastic pollution. The present contribution could be represent a good step forward. However, the Authors should address the following comments before reconsidering the manuscript for publications
General comments
Modelling the transport of plastic debris is known to be a challenging task for the complex mechanisms involved. The inertial character of the debris poses the major problems and several (simplified) solutions have been proposed in the last years. TrackMPD is one of the available suite for modelling the Lagrangian transport of particles in ocean and coastal environments. A major concern remains how the inertial transport is modelled.
Focusing the attention on the hydrodynamic context described in the manuscript, a series of papers have been recently published on the transport of inertial particles under the action of waves. For examples the following contributions:
DiBenedetto, M. H., Ouellette, N. T., & Koseff, J. R. (2018). Transport of anisotropic particles under waves. Journal of Fluid Mechanics, 837, 320-340.
DiBenedetto, M. H., & Ouellette, N. T. (2018). Preferential orientation of spheroidal particles in wavy flow. Journal of Fluid Mechanics, 856, 850-869.
DiBenedetto, M. H., Koseff, J. R., & Ouellette, N. T. (2019). Orientation dynamics of nonspherical particles under surface gravity waves. Physical Review Fluids, 4(3), 034301.
De Leo, A., & Stocchino, A. (2022). Dispersion of heavy particles under sea waves. Physics of Fluids, 34(1).
DiBenedetto, M. H., Clark, L. K., & Pujara, N. (2022). Enhanced settling and dispersion of inertial particles in surface waves. Journal of Fluid Mechanics, 936, A38.
Clark, L. K., DiBenedetto, M. H., Ouellette, N. T., & Koseff, J. R. (2023). Dispersion of finite-size, non-spherical particles by waves and currents. Journal of Fluid Mechanics, 954, A3.
The above list is not complete, but represent a good example on how the inertial behaviour of the plastic debris can be mathematically modelled. The main conclusions of the studies on inertial particles and waves are: settling/rising is strongly enhanced by the inertial effects and the spreading on inertial particles is almost never found to be fully diffusive.
It is clear that it is almost impossible to fully include the inertial effects in approaches like the ones used in TrackMPD. However the Authors should include at least part of the reference listed above and discuss the limitations of their modelling approach. For example, the assumption that the dispersion of the plastic debris is diffusive should be better explained. Closing the fluxes using a diffusion-type coefficient is a common approach even if it has been demonstrated that heavy particles don't follow a Brownian dispersion regime. However, it is convenient to introduce horizontal and vertical coefficient, as done in the present study. TrackMPD seems to include only gravitational effects through a settling velocity of the debris. No other inertial effect are included or modelled. The Authors should discuss these aspects in more details
Another major concern regards the treatment of the particle transport close to the bed (sections 2.2.2 and 2.2.3). To what extent is reasonable to apply concepts developed for sediment transport to the present case? the empirical formulations used in the case of bedload or suspended load are developed considering relatively high concentration of sediment, especially in the case of bedload. with this assumption, the emprical models are not lagrangian models, but describe the average behavior of a certain mass of sediment. Indeed, the sediment continuity equation (Exner equation) is written in Eulerian terms. Similarly, suspended sediment transport is described using the advection diffusion equation for the sediment concentration. Plastic debris concentrations are usually much less compared to suspended sediment concentration and, of course, at the bed. The Authors used the Soulsby formulation designed for sands. Is it reasonable to assume that the estimate of the critical shear stresses (bedload and suspension) is valid also in case of plastic debris? There are no evidence of this in the literature, even if this approach is commonly used. Moreover, it is not clear the formulation of the bedload transport and how, provided the mobility condition, the bedload is described in terms of an acceleration term.
Detailed comments
Introduction
Include at least part of the reference provided above and discuss in more details what have been done in order to describe the inertial particles transport under waves
Section 2 Methods
1. Please provide more details on the Lagrangian particle transport model equation implemented in TrackMPD
2. Please provide more details on the mesh used in SWASH and the integration time step used for both SWASH and TrackMPD.
3. What are the Stokes time of the simulated particles? is it comparable to the integration time step?
4. line 165, please provide the definition of D*
5. Section 2.2.3: what is the difference between Ch and Cd?
6. line 185. It is not clear the reason why the eddy viscosity is replaced the fluid viscosity
7. Section 3.1. The results shown in Figure 2 are not clear. How particles have been simulated to obtain the results? Are the grey lines the particle trajectories? If yes, why the black points represent the final positions? How long was the simulations in terms of wave periods? and what is the influence on the final results?
8. section 3.3. The discussion should be improved providing more details on the physical meaning of the transport mechanisms
9. line 360. it is not true that the interaction between plastic particles and waves received limited attention, see the reference listed above. Please rephrase
10 a section on the limitation of the approach should be added.
Citation: https://doi.org/10.5194/gmd-2024-100-RC2 -
AC3: 'Reply on RC2', Isabel Jalon-Rojas, 23 Sep 2024
We thank the reviewer for carefully reading our manuscript, for contrastive comments, and for useful suggestions. Below, we address each of the points and outline how we plan to revise the manuscript accordingly:
General comments
Modelling the transport of plastic debris is known to be a challenging task for the complex mechanisms involved. The inertial character of the debris poses the major problems and several (simplified) solutions have been proposed in the last years. TrackMPD is one of the available suite for modelling the Lagrangian transport of particles in ocean and coastal environments. A major concern remains how the inertial transport is modelled.
Focusing the attention on the hydrodynamic context described in the manuscript, a series of papers have been recently published on the transport of inertial particles under the action of waves. For examples the following contributions:
DiBenedetto, M. H., Ouellette, N. T., & Koseff, J. R. (2018). Transport of anisotropic particles under waves. Journal of Fluid Mechanics, 837, 320-340.
DiBenedetto, M. H., & Ouellette, N. T. (2018). Preferential orientation of spheroidal particles in wavy flow. Journal of Fluid Mechanics, 856, 850-869.
DiBenedetto, M. H., Koseff, J. R., & Ouellette, N. T. (2019). Orientation dynamics of nonspherical particles under surface gravity waves. Physical Review Fluids, 4(3), 034301.
De Leo, A., & Stocchino, A. (2022). Dispersion of heavy particles under sea waves. Physics of Fluids, 34(1).
DiBenedetto, M. H., Clark, L. K., & Pujara, N. (2022). Enhanced settling and dispersion of inertial particles in surface waves. Journal of Fluid Mechanics, 936, A38.
Clark, L. K., DiBenedetto, M. H., Ouellette, N. T., & Koseff, J. R. (2023). Dispersion of finite-size, non-spherical particles by waves and currents. Journal of Fluid Mechanics, 954, A3.
The above list is not complete, but represent a good example on how the inertial behaviour of the plastic debris can be mathematically modelled. The main conclusions of the studies on inertial particles and waves are: settling/rising is strongly enhanced by the inertial effects and the spreading on inertial particles is almost never found to be fully diffusive.
It is clear that it is almost impossible to fully include the inertial effects in approaches like the ones used in TrackMPD. However the Authors should include at least part of the reference listed above and discuss the limitations of their modelling approach. For example, the assumption that the dispersion of the plastic debris is diffusive should be better explained. Closing the fluxes using a diffusion-type coefficient is a common approach even if it has been demonstrated that heavy particles don't follow a Brownian dispersion regime. However, it is convenient to introduce horizontal and vertical coefficient, as done in the present study. TrackMPD seems to include only gravitational effects through a settling velocity of the debris. No other inertial effect are included or modelled. The Authors should discuss these aspects in more details.
We thank the reviewer for their valuable comment and the list of publications on the inertial behavior of particles. It is true that our approach neglects particle inertia, and this should indeed be more explicitly stated and discussed. Our current discussion touches on this point (lines 414–421), where we reference the experimental study by Alsina et al. (2020), which suggests that, in the nearshore shoaling region—our area of interest—factors aside from buoyancy, such as inertial effects, minimally influence the net drift of low-density particles, likely due to high turbulence and other complex processes in this region (air-water two-phase mixing, roller). The precise contribution of each of the involved processes remains to be established, following the growing research effort on the topic including the references provided by the reviewer. Our modelling approach should therefore be considered as a first bulk approach on the microplastic transport in the surf zone. We also noted that future model developments could consider incorporating the effects of particle drag on advection, as proposed by Stocchino et al. We agree that the studies cited by the reviewer provide valuable insights into particle inertia, and we will expand our discussion in the next version of the manuscript using these references to inform potential future developments.
Another major concern regards the treatment of the particle transport close to the bed (sections 2.2.2 and 2.2.3). To what extent is reasonable to apply concepts developed for sediment transport to the present case? the empirical formulations used in the case of bedload or suspended load are developed considering relatively high concentration of sediment, especially in the case of bedload. with this assumption, the emprical models are not lagrangian models, but describe the average behavior of a certain mass of sediment. Indeed, the sediment continuity equation (Exner equation) is written in Eulerian terms. Similarly, suspended sediment transport is described using the advection diffusion equation for the sediment concentration. Plastic debris concentrations are usually much less compared to suspended sediment concentration and, of course, at the bed. The Authors used the Soulsby formulation designed for sands. Is it reasonable to assume that the estimate of the critical shear stresses (bedload and suspension) is valid also in case of plastic debris? There are no evidence of this in the literature, even if this approach is commonly used. Moreover, it is not clear the formulation of the bedload transport and how, provided the mobility condition, the bedload is described in terms of an acceleration term.
We agree that there are some differences, but also similitudes, in the transport of sediment and microplastics as has been reviewed in the work by Waldschlager et al. (2022). We also concur with the reviewer that there is a lack of empirical work on microplastics, particularly for bed load transport. Therefore, we made the assumption of using a similar approach to that used for sediments as an initial step, which can be refined in the future as more empirical data become available. We will make this limitation clearer in the next version of the manuscript.
Regarding the calculation of the bed shear stress of microplastics, there are now several empirical studies available. For example, Waldschlager and Schttrumpf (2019) proposed a correction of Shield formulation to account for the sediment bed properties. Goral et al. (2023) proposed a new framework so that the incipient motion conditions for microplastic particles lying on a sediment bed are, for the first time, reconciled with the classical Shields diagram. One author of the present manuscript is also working on this topic (https://meetingorganizer.copernicus.org/EGU24/EGU24-20316.html). Given the current state of knowledge, TrackMPD include the Soulsby/Shield approach and the correcting proposed by Waldschlager and Schttrumpf (2019). Given that there is no sediment bed in the experimental case that we are reproducing, we have opted for using the Soulsby approach for this specific case. The approaches included in TrackMPD can be refined in future as more empirical data and unified formulations will be available. We will discuss these elements in more detail in the revised manuscript, addressing both the current limitations and the potential for future refinements based on emerging empirical data.
Waldschläger, K., Brückner, M. Z., Almroth, B. C., Hackney, C. R., Adyel, T. M., Alimi, O. S., ... & Wu, N. (2022). Learning from natural sediments to tackle microplastics challenges: A multidisciplinary perspective. Earth-Science Reviews, 228, 104021.
Waldschläger, K. and Schttrumpf, H.: Erosion behavior of different microplastic particles in comparison to natural sediments, Environmental science & technology, 53, 13 219–13 227, 2019b.
Goral, K. D., Guler, H. G., Larsen, B. E., Carstensen, S., Christensen, E. D., Kerpen, N. B., ... & Fuhrman, D. R. (2023). Shields diagram and the incipient motion of microplastic particles. Environmental Science & Technology, 57(25), 9362-9375.
Detailed comments
Introduction
Include at least part of the reference provided above and discuss in more details what have been done in order to describe the inertial particles transport under waves
Yes, we will complete the description of the state of the art using references provided above.
Section 2 Methods
- Please provide more details on the Lagrangian particle transport model equation implemented in TrackMPD
The advection-dispersion equations and numerical schemes implemented in TrackMPD were comprehensively detailed in the original publication (Jalón-Rojas et al., 2019). In this paper, we focus on presenting the new developments of the model rather than repeating previously published information. However, we do provide a summary of the key aspects in lines 89-91: “The model employs a 4th-order Runge-Kutta scheme to accurately advect virtual particles through a set of velocity fields. A random-walk approach is implemented to simulate the turbulent motion of particles in both the horizontal and vertical directions as a function of the horizontal and vertical diffusivity coefficients (see Section 2.2.1 for more details).” For more comprehensive information on the Lagrangian transport equations, we direct interested readers to the earlier publication (Jalón-Rojas et al., 2019).
- Please provide more details on the mesh used in SWASH and the integration time step used for both SWASH and TrackMPD.
The mesh and time steps used in SWASH and TrackMPD were provided in lines 193-195: The computational grid consisted of 175 points in the horizontal direction (resolution 3.45cm) and 15 sigma-layers in the vertical direction. A time step of 0.05 seconds was selected for both hydrodynamic and particle tracking simulations.
- What are the Stokes time of the simulated particles? is it comparable to the integration time step?
The integration time step (0.05 seconds) is much smaller than the Stokes time for all particles (ranging from 9.4 to 35.7 seconds, calculated as τ_Stokes=dp2/(12νβ), indicating that the model's time step is adequately small to resolve the dynamics of the particles accurately.
- line 165, please provide the definition of D*
We will do it.
- Section 2.2.3: what is the difference between Ch and Cd?
Ch refers to the particle drag coefficient, while Cd refers to the bottom drag coefficient. We will complete the description of Ch in line 178 to avoid any confusion between the two terms.
- line 185. It is not clear the reason why the eddy viscosity is replaced the fluid viscosity
The following explaination has been included in the Discussion section :
“Estimating the Stokes layer thickness as for the present case leads to a value about 0.6mm, which is larger than the thickness of tested fibers and sheets. For the spheres, the top of the particles is expected to rise above the laminar layer. However, as the exact distribution of small-scale viscous/turbulent shear remains unreachable to the present dataset, the fluid viscosity is used for each particle type to ensure robust comparisons.”
- Section 3.1. The results shown in Figure 2 are not clear. How particles have been simulated to obtain the results? Are the grey lines the particle trajectories? If yes, why the black points represent the final positions? How long was the simulations in terms of wave periods? and what is the influence on the final results?
We will clarify the color coding in the figures (gray lines for trajectories and black points for final positions). The simulations run for a sufficiently long duration (10 minutes, compared to the 1.2-second wave period) to ensure that the results are not influenced. We will include this information in the revised version. Sensitivity tests were even conducted to confirm this.
- section 3.3. The discussion should be improved providing more details on the physical meaning of the transport mechanisms
We will add a few more details on the physical meaning of the transport mechanisms. For example:
“superimposed to the oscillating backwards and forwards transport due to the wave motion, the plastic particles tend to spend more time in the faster onshore-moving layer underneath the crest than in the slower offshore-moving layer below the trough. This effect, which is the well-known Stokes drift, is enhanced by the increased non-linearity and asymmetry of shoaling and surf zone waves. As shown in Figure 5.a-f, these spheres predominantly traveled in the upper water layer, closely following the water surface in the vertical coordinate, as their high buoyancy prevents significant dispersion due to turbulence. Consequently, they followed the net drift velocity aligned with the wave propagation direction, i.e. the Stokes drift (van den Bremer and Breivik, 2018) (see residual velocity in Figure 1.c).”
“In surface layer, the driving mechanisms are the same that affected low-density particles described herebefore. When spreading throughout the water column due buoyancy-driven settling and vertical turbulent mixing, the particles are exposed to a offshore-directed return current, i.e. the classical undertow, which compensates for the wave-induced surface mass flux (Fig. 1.b-c).”
- line 360. it is not true that the interaction between plastic particles and waves received limited attention, see the reference listed above. Please rephrase
We agree that the interactions between plastic particles and waves have not received limited attention, as demonstrated by the studies listed in the Introduction and throughout the paper including the references raised by the reviewer. However, in this sentence in question, we specifically refer to wave non-linear processes (e.g., wave asymmetry and skewness) in shallow waters, which have received less attention. To better clarify it, we will add “in shallow waters” after “wave non-linear processes (e.g. Martins et al., 2020)”.
10 a section on the limitation of the approach should be added
The main limitations of our approach and potential future improvements were already included in the discussion section, lines 410-422. Nevertheless, we will try to complete this section including some points raised in general comments to ensure that the limitations (and future developments) are more explicitly detailed.
Citation: https://doi.org/10.5194/gmd-2024-100-AC3
-
AC3: 'Reply on RC2', Isabel Jalon-Rojas, 23 Sep 2024
Status: closed
-
CEC1: 'Comment on gmd-2024-100', Astrid Kerkweg, 15 Jul 2024
Dear authors,
in my role as Executive editor of GMD, I would like to bring to your attention our Editorial version 1.2: https://www.geosci-model-dev.net/12/2215/2019/
This highlights some requirements of papers published in GMD, which is also available on the GMD website in the ‘Manuscript Types’ section: http://www.geoscientific-model-development.net/submission/manuscript_types.html
In particular, please note that for your paper, the following requirements have not been met in the Discussions paper:
- "The main paper must give the model name and version number (or other unique identifier) in the title."
In order to simplify reference to your developments, please add a model name (and/or its acronym) and a version number in the title of your article in your revised submission to GMD.
Yours, Astrid Kerkweg
Citation: https://doi.org/10.5194/gmd-2024-100-CEC1 -
AC1: 'Reply on CEC1', Isabel Jalon-Rojas, 11 Sep 2024
Dear Editor,
Following the journal's requirements, we will modify the title to include the model name and version number: 'A wave-resolving 2DV Lagrangian approach to model microplastic transport in the nearshore based on TrackMPD v3.0.'"
Best regards,
Isabel
Citation: https://doi.org/10.5194/gmd-2024-100-AC1
-
RC1: 'Comment on gmd-2024-100', Anonymous Referee #1, 16 Jul 2024
This paper presents a novel approach to numerically study the transport of buoyant and non-buoyant microplastics in wave-dominated nearshore areas, a highly relevant and innovative topic. The Authors combine the SWASH numerical model to resolve surf zone hydrodynamics with the Lagrangian model TrakMPD to track microplastics. They enhanced the Lagrangian model to better represent turbulent mixing with spatial gradients from wave breaking and depth variations and included processes such as deposition, resuspension, and bedload transport.
The authors reproduced one of Forsberg's windless hydrodynamic laboratory experiments, parameterizing the model based on existing literature Under this condition, they evaluated the transport and dispersion of six types of microplastics (spheres, fibres, and sheets, both buoyant and non-buoyant) and found a good match between numerical results and laboratory findings in the final concentration distribution across different beach profile areas. The main discrepancy occurred with bottom-transported microplastics, which escaped the breaking zone into the adjacent shoaling zone. This is likely due to a slight misrepresentation of the experimental undertow current, the temporal evolution of vertical mixing, or inaccuracies in defining the compartment limits, as the authors suggest. Subsequently, the Authors conducted a sensitivity analysis of model parameters to explore the influence of turbulence from wave breaking, bedload transport, and resuspension on microplastic transport. They also investigated mechanisms driving microplastic transport, such as Stokes drift, wave asymmetry, and particle settling and rising velocities.
The topic addressed by this paper is highly relevant and potentially a valuable contribution to the Scientific Community. The paper is concise, well-organized, and well-written. To further enhance the manuscript, the following minor aspects and suggestions should be addressed:
Minor comments
L131-138: In Section 2.2.1, general aspects and model improvements are described, without addressing the specific configuration for the present study. Therefore, I would suggest relocating these lines (L131-137, “For the present simulation…. uniform throughout the domain.”) closer to line L223, where the chosen model parameters for the study are detailed.
L133: Please replace “surfzone” with “surf zone”.
L136: Here, a sentence could be added to justify the hypothesis that the vertical diffusion coefficient for microplastics equals the eddy diffusivity, given that the small size of microplastics ensures their behaviour as passive tracers predominantly governed by fluid turbulence.
L190: I would explicitly state that one of the two windless hydrodynamic conditions evaluated by Forsberg was replicated numerically.
L199: It is recommended to include a quantitative/statistical comparison between numerical and laboratory wave heights and an assessment of how much the model underestimates dissipation in the surf zone.
L207: I believe the Jalon-Rojas method (2022) estimates rising and settling velocities in calm water. De Leo et al. (2021; https://doi.org/10.3390/jmse9020142) found that settling velocities increase in the presence of waves. I would add a sentence noting this aspect.
L212 & L245: It should be clarified whether the results presented in the manuscript are the average of all runs or from a representative run.
P11, Section 3 (Results): Similarly to line L261, where it states “...54% and 67% of low-density fibres were gradually transported onshore...”, I would like to see more quantification in the description of spheres and sheets throughout this section.
Citation: https://doi.org/10.5194/gmd-2024-100-RC1 -
AC2: 'Reply on RC1', Isabel Jalon-Rojas, 11 Sep 2024
We thank the reviewer for carefully reading our manuscript, for positive comments, and useful suggestions. Below, we address each of the points and outline how we plan to revise the manuscript accordingly:
Minor comments
L131-138: In Section 2.2.1, general aspects and model improvements are described, without addressing the specific configuration for the present study. Therefore, I would suggest relocating these lines (L131-137, “For the present simulation…. uniform throughout the domain.”) closer to line L223, where the chosen model parameters for the study are detailed.
Yes, we will relocate these lines in Section 2.4 as suggested
L133: Please replace “surfzone” with “surf zone”.
Yes, we will replace it.
L136: Here, a sentence could be added to justify the hypothesis that the vertical diffusion coefficient for microplastics equals the eddy diffusivity, given that the small size of microplastics ensures their behaviour as passive tracers predominantly governed by fluid turbulence.
Good suggestion, we will include this explanation in the text as suggested
L190: I would explicitly state that one of the two windless hydrodynamic conditions evaluated by Forsberg was replicated numerically.
Yes, we will detail that we have replicated one of the two windless hydrodynamic conditions evaluated by Forsberg.
L199: It is recommended to include a quantitative/statistical comparison between numerical and laboratory wave heights and an assessment of how much the model underestimates dissipation in the surf zone.
We will modify the sentence to include a quantitative comparison: "Numerical wave heights are generally in good agreement with experimental measurements, as illustrated in Figure 1a. Specifically, the model tends to slightly underestimate surf zone dissipation, leading to a root mean square error (RMSE) of 2.2 cm and an average negative bias of 0.7 cm for wave heights in this region. Despite this minor discrepancy, the overall agreement in wave heights supports using SWASH simulations as a reliable hydrodynamic forcing for TrackMPD."
L207: I believe the Jalon-Rojas method (2022) estimates rising and settling velocities in calm water. De Leo et al. (2021; https://doi.org/10.3390/jmse9020142) found that settling velocities increase in the presence of waves. I would add a sentence noting this aspect.
De Leo et al. (2021) proposed a formulation for settling velocity to account for the influence of waves on vertical trajectories. In our study, the vertical trajectories of particles are already influenced by the wave-resolved hydrodynamics within the SWASH model. Since the effects of waves are inherently captured by the hydrodynamics, we do not need to directly adjust the settling velocity for wave-induced variations, so we use values in calm water. The adjustment proposed by Leo et al. (2021) could indeed be relevant for large-scale models that do not explicitly resolve wave dynamics but, in our case, the wave influences are accounted for through the hydrodynamic forcing.
L212 & L245: It should be clarified whether the results presented in the manuscript are the average of all runs or from a representative run.
Figures 2.b and 3.a-b illustrate the main results of the manuscript, showing the average values and standard deviations from five runs for each scenario. As noted in the text: “to ensure the robustness of our simulations, we conducted five simulations for each scenario, and the results consistently exhibited only minor variability in the number of particles within each region. This variability was within the same order of magnitude as that observed in the experiments as indicated by the error bars in Figure 2.” In Figures 2.a and 3.c, we present the trajectories and final positions from one of the runs, which, given the low variability between runs, is representative of the scenario. We will clarify this last point in the text.
P11, Section 3 (Results): Similarly to line L261, where it states “...54% and 67% of low-density fibres were gradually transported onshore...”, I would like to see more quantification in the description of spheres and sheets throughout this section.
We will include this quantification for low-density sheets as well. For low-density spheres and all high-density particles, nearly all particles ended up in a single compartment, as mentioned in the text.
Citation: https://doi.org/10.5194/gmd-2024-100-AC2
-
AC2: 'Reply on RC1', Isabel Jalon-Rojas, 11 Sep 2024
-
RC2: 'Comment on gmd-2024-100', Anonymous Referee #2, 02 Sep 2024
The manuscript describes the modifications proposed to an existing numerical modelling approach (TrackMPD) for plastic debris with the aim to improve the transport in coastal zones. The modified numerical approach was then benchmarked against experimental observations of the debris transport in the nearshore zone (breaking zone, swash zone) recently published by one of the Authors. Differently from the original formulations of TrackMPD, in the present study the hydrodynamics is solved used SWASH, a new formulation is proposed for the vertical diffusivity (linked to the wave breaking), and the computation of the transport at the bed. The topic is of interest from the scientific point of view and with possible high impact in terms of marine pollution. In fact, a reliable numerical model able to predict the plastic debris trajectories is absolutely need to contrast plastic pollution. The present contribution could be represent a good step forward. However, the Authors should address the following comments before reconsidering the manuscript for publications
General comments
Modelling the transport of plastic debris is known to be a challenging task for the complex mechanisms involved. The inertial character of the debris poses the major problems and several (simplified) solutions have been proposed in the last years. TrackMPD is one of the available suite for modelling the Lagrangian transport of particles in ocean and coastal environments. A major concern remains how the inertial transport is modelled.
Focusing the attention on the hydrodynamic context described in the manuscript, a series of papers have been recently published on the transport of inertial particles under the action of waves. For examples the following contributions:
DiBenedetto, M. H., Ouellette, N. T., & Koseff, J. R. (2018). Transport of anisotropic particles under waves. Journal of Fluid Mechanics, 837, 320-340.
DiBenedetto, M. H., & Ouellette, N. T. (2018). Preferential orientation of spheroidal particles in wavy flow. Journal of Fluid Mechanics, 856, 850-869.
DiBenedetto, M. H., Koseff, J. R., & Ouellette, N. T. (2019). Orientation dynamics of nonspherical particles under surface gravity waves. Physical Review Fluids, 4(3), 034301.
De Leo, A., & Stocchino, A. (2022). Dispersion of heavy particles under sea waves. Physics of Fluids, 34(1).
DiBenedetto, M. H., Clark, L. K., & Pujara, N. (2022). Enhanced settling and dispersion of inertial particles in surface waves. Journal of Fluid Mechanics, 936, A38.
Clark, L. K., DiBenedetto, M. H., Ouellette, N. T., & Koseff, J. R. (2023). Dispersion of finite-size, non-spherical particles by waves and currents. Journal of Fluid Mechanics, 954, A3.
The above list is not complete, but represent a good example on how the inertial behaviour of the plastic debris can be mathematically modelled. The main conclusions of the studies on inertial particles and waves are: settling/rising is strongly enhanced by the inertial effects and the spreading on inertial particles is almost never found to be fully diffusive.
It is clear that it is almost impossible to fully include the inertial effects in approaches like the ones used in TrackMPD. However the Authors should include at least part of the reference listed above and discuss the limitations of their modelling approach. For example, the assumption that the dispersion of the plastic debris is diffusive should be better explained. Closing the fluxes using a diffusion-type coefficient is a common approach even if it has been demonstrated that heavy particles don't follow a Brownian dispersion regime. However, it is convenient to introduce horizontal and vertical coefficient, as done in the present study. TrackMPD seems to include only gravitational effects through a settling velocity of the debris. No other inertial effect are included or modelled. The Authors should discuss these aspects in more details
Another major concern regards the treatment of the particle transport close to the bed (sections 2.2.2 and 2.2.3). To what extent is reasonable to apply concepts developed for sediment transport to the present case? the empirical formulations used in the case of bedload or suspended load are developed considering relatively high concentration of sediment, especially in the case of bedload. with this assumption, the emprical models are not lagrangian models, but describe the average behavior of a certain mass of sediment. Indeed, the sediment continuity equation (Exner equation) is written in Eulerian terms. Similarly, suspended sediment transport is described using the advection diffusion equation for the sediment concentration. Plastic debris concentrations are usually much less compared to suspended sediment concentration and, of course, at the bed. The Authors used the Soulsby formulation designed for sands. Is it reasonable to assume that the estimate of the critical shear stresses (bedload and suspension) is valid also in case of plastic debris? There are no evidence of this in the literature, even if this approach is commonly used. Moreover, it is not clear the formulation of the bedload transport and how, provided the mobility condition, the bedload is described in terms of an acceleration term.
Detailed comments
Introduction
Include at least part of the reference provided above and discuss in more details what have been done in order to describe the inertial particles transport under waves
Section 2 Methods
1. Please provide more details on the Lagrangian particle transport model equation implemented in TrackMPD
2. Please provide more details on the mesh used in SWASH and the integration time step used for both SWASH and TrackMPD.
3. What are the Stokes time of the simulated particles? is it comparable to the integration time step?
4. line 165, please provide the definition of D*
5. Section 2.2.3: what is the difference between Ch and Cd?
6. line 185. It is not clear the reason why the eddy viscosity is replaced the fluid viscosity
7. Section 3.1. The results shown in Figure 2 are not clear. How particles have been simulated to obtain the results? Are the grey lines the particle trajectories? If yes, why the black points represent the final positions? How long was the simulations in terms of wave periods? and what is the influence on the final results?
8. section 3.3. The discussion should be improved providing more details on the physical meaning of the transport mechanisms
9. line 360. it is not true that the interaction between plastic particles and waves received limited attention, see the reference listed above. Please rephrase
10 a section on the limitation of the approach should be added.
Citation: https://doi.org/10.5194/gmd-2024-100-RC2 -
AC3: 'Reply on RC2', Isabel Jalon-Rojas, 23 Sep 2024
We thank the reviewer for carefully reading our manuscript, for contrastive comments, and for useful suggestions. Below, we address each of the points and outline how we plan to revise the manuscript accordingly:
General comments
Modelling the transport of plastic debris is known to be a challenging task for the complex mechanisms involved. The inertial character of the debris poses the major problems and several (simplified) solutions have been proposed in the last years. TrackMPD is one of the available suite for modelling the Lagrangian transport of particles in ocean and coastal environments. A major concern remains how the inertial transport is modelled.
Focusing the attention on the hydrodynamic context described in the manuscript, a series of papers have been recently published on the transport of inertial particles under the action of waves. For examples the following contributions:
DiBenedetto, M. H., Ouellette, N. T., & Koseff, J. R. (2018). Transport of anisotropic particles under waves. Journal of Fluid Mechanics, 837, 320-340.
DiBenedetto, M. H., & Ouellette, N. T. (2018). Preferential orientation of spheroidal particles in wavy flow. Journal of Fluid Mechanics, 856, 850-869.
DiBenedetto, M. H., Koseff, J. R., & Ouellette, N. T. (2019). Orientation dynamics of nonspherical particles under surface gravity waves. Physical Review Fluids, 4(3), 034301.
De Leo, A., & Stocchino, A. (2022). Dispersion of heavy particles under sea waves. Physics of Fluids, 34(1).
DiBenedetto, M. H., Clark, L. K., & Pujara, N. (2022). Enhanced settling and dispersion of inertial particles in surface waves. Journal of Fluid Mechanics, 936, A38.
Clark, L. K., DiBenedetto, M. H., Ouellette, N. T., & Koseff, J. R. (2023). Dispersion of finite-size, non-spherical particles by waves and currents. Journal of Fluid Mechanics, 954, A3.
The above list is not complete, but represent a good example on how the inertial behaviour of the plastic debris can be mathematically modelled. The main conclusions of the studies on inertial particles and waves are: settling/rising is strongly enhanced by the inertial effects and the spreading on inertial particles is almost never found to be fully diffusive.
It is clear that it is almost impossible to fully include the inertial effects in approaches like the ones used in TrackMPD. However the Authors should include at least part of the reference listed above and discuss the limitations of their modelling approach. For example, the assumption that the dispersion of the plastic debris is diffusive should be better explained. Closing the fluxes using a diffusion-type coefficient is a common approach even if it has been demonstrated that heavy particles don't follow a Brownian dispersion regime. However, it is convenient to introduce horizontal and vertical coefficient, as done in the present study. TrackMPD seems to include only gravitational effects through a settling velocity of the debris. No other inertial effect are included or modelled. The Authors should discuss these aspects in more details.
We thank the reviewer for their valuable comment and the list of publications on the inertial behavior of particles. It is true that our approach neglects particle inertia, and this should indeed be more explicitly stated and discussed. Our current discussion touches on this point (lines 414–421), where we reference the experimental study by Alsina et al. (2020), which suggests that, in the nearshore shoaling region—our area of interest—factors aside from buoyancy, such as inertial effects, minimally influence the net drift of low-density particles, likely due to high turbulence and other complex processes in this region (air-water two-phase mixing, roller). The precise contribution of each of the involved processes remains to be established, following the growing research effort on the topic including the references provided by the reviewer. Our modelling approach should therefore be considered as a first bulk approach on the microplastic transport in the surf zone. We also noted that future model developments could consider incorporating the effects of particle drag on advection, as proposed by Stocchino et al. We agree that the studies cited by the reviewer provide valuable insights into particle inertia, and we will expand our discussion in the next version of the manuscript using these references to inform potential future developments.
Another major concern regards the treatment of the particle transport close to the bed (sections 2.2.2 and 2.2.3). To what extent is reasonable to apply concepts developed for sediment transport to the present case? the empirical formulations used in the case of bedload or suspended load are developed considering relatively high concentration of sediment, especially in the case of bedload. with this assumption, the emprical models are not lagrangian models, but describe the average behavior of a certain mass of sediment. Indeed, the sediment continuity equation (Exner equation) is written in Eulerian terms. Similarly, suspended sediment transport is described using the advection diffusion equation for the sediment concentration. Plastic debris concentrations are usually much less compared to suspended sediment concentration and, of course, at the bed. The Authors used the Soulsby formulation designed for sands. Is it reasonable to assume that the estimate of the critical shear stresses (bedload and suspension) is valid also in case of plastic debris? There are no evidence of this in the literature, even if this approach is commonly used. Moreover, it is not clear the formulation of the bedload transport and how, provided the mobility condition, the bedload is described in terms of an acceleration term.
We agree that there are some differences, but also similitudes, in the transport of sediment and microplastics as has been reviewed in the work by Waldschlager et al. (2022). We also concur with the reviewer that there is a lack of empirical work on microplastics, particularly for bed load transport. Therefore, we made the assumption of using a similar approach to that used for sediments as an initial step, which can be refined in the future as more empirical data become available. We will make this limitation clearer in the next version of the manuscript.
Regarding the calculation of the bed shear stress of microplastics, there are now several empirical studies available. For example, Waldschlager and Schttrumpf (2019) proposed a correction of Shield formulation to account for the sediment bed properties. Goral et al. (2023) proposed a new framework so that the incipient motion conditions for microplastic particles lying on a sediment bed are, for the first time, reconciled with the classical Shields diagram. One author of the present manuscript is also working on this topic (https://meetingorganizer.copernicus.org/EGU24/EGU24-20316.html). Given the current state of knowledge, TrackMPD include the Soulsby/Shield approach and the correcting proposed by Waldschlager and Schttrumpf (2019). Given that there is no sediment bed in the experimental case that we are reproducing, we have opted for using the Soulsby approach for this specific case. The approaches included in TrackMPD can be refined in future as more empirical data and unified formulations will be available. We will discuss these elements in more detail in the revised manuscript, addressing both the current limitations and the potential for future refinements based on emerging empirical data.
Waldschläger, K., Brückner, M. Z., Almroth, B. C., Hackney, C. R., Adyel, T. M., Alimi, O. S., ... & Wu, N. (2022). Learning from natural sediments to tackle microplastics challenges: A multidisciplinary perspective. Earth-Science Reviews, 228, 104021.
Waldschläger, K. and Schttrumpf, H.: Erosion behavior of different microplastic particles in comparison to natural sediments, Environmental science & technology, 53, 13 219–13 227, 2019b.
Goral, K. D., Guler, H. G., Larsen, B. E., Carstensen, S., Christensen, E. D., Kerpen, N. B., ... & Fuhrman, D. R. (2023). Shields diagram and the incipient motion of microplastic particles. Environmental Science & Technology, 57(25), 9362-9375.
Detailed comments
Introduction
Include at least part of the reference provided above and discuss in more details what have been done in order to describe the inertial particles transport under waves
Yes, we will complete the description of the state of the art using references provided above.
Section 2 Methods
- Please provide more details on the Lagrangian particle transport model equation implemented in TrackMPD
The advection-dispersion equations and numerical schemes implemented in TrackMPD were comprehensively detailed in the original publication (Jalón-Rojas et al., 2019). In this paper, we focus on presenting the new developments of the model rather than repeating previously published information. However, we do provide a summary of the key aspects in lines 89-91: “The model employs a 4th-order Runge-Kutta scheme to accurately advect virtual particles through a set of velocity fields. A random-walk approach is implemented to simulate the turbulent motion of particles in both the horizontal and vertical directions as a function of the horizontal and vertical diffusivity coefficients (see Section 2.2.1 for more details).” For more comprehensive information on the Lagrangian transport equations, we direct interested readers to the earlier publication (Jalón-Rojas et al., 2019).
- Please provide more details on the mesh used in SWASH and the integration time step used for both SWASH and TrackMPD.
The mesh and time steps used in SWASH and TrackMPD were provided in lines 193-195: The computational grid consisted of 175 points in the horizontal direction (resolution 3.45cm) and 15 sigma-layers in the vertical direction. A time step of 0.05 seconds was selected for both hydrodynamic and particle tracking simulations.
- What are the Stokes time of the simulated particles? is it comparable to the integration time step?
The integration time step (0.05 seconds) is much smaller than the Stokes time for all particles (ranging from 9.4 to 35.7 seconds, calculated as τ_Stokes=dp2/(12νβ), indicating that the model's time step is adequately small to resolve the dynamics of the particles accurately.
- line 165, please provide the definition of D*
We will do it.
- Section 2.2.3: what is the difference between Ch and Cd?
Ch refers to the particle drag coefficient, while Cd refers to the bottom drag coefficient. We will complete the description of Ch in line 178 to avoid any confusion between the two terms.
- line 185. It is not clear the reason why the eddy viscosity is replaced the fluid viscosity
The following explaination has been included in the Discussion section :
“Estimating the Stokes layer thickness as for the present case leads to a value about 0.6mm, which is larger than the thickness of tested fibers and sheets. For the spheres, the top of the particles is expected to rise above the laminar layer. However, as the exact distribution of small-scale viscous/turbulent shear remains unreachable to the present dataset, the fluid viscosity is used for each particle type to ensure robust comparisons.”
- Section 3.1. The results shown in Figure 2 are not clear. How particles have been simulated to obtain the results? Are the grey lines the particle trajectories? If yes, why the black points represent the final positions? How long was the simulations in terms of wave periods? and what is the influence on the final results?
We will clarify the color coding in the figures (gray lines for trajectories and black points for final positions). The simulations run for a sufficiently long duration (10 minutes, compared to the 1.2-second wave period) to ensure that the results are not influenced. We will include this information in the revised version. Sensitivity tests were even conducted to confirm this.
- section 3.3. The discussion should be improved providing more details on the physical meaning of the transport mechanisms
We will add a few more details on the physical meaning of the transport mechanisms. For example:
“superimposed to the oscillating backwards and forwards transport due to the wave motion, the plastic particles tend to spend more time in the faster onshore-moving layer underneath the crest than in the slower offshore-moving layer below the trough. This effect, which is the well-known Stokes drift, is enhanced by the increased non-linearity and asymmetry of shoaling and surf zone waves. As shown in Figure 5.a-f, these spheres predominantly traveled in the upper water layer, closely following the water surface in the vertical coordinate, as their high buoyancy prevents significant dispersion due to turbulence. Consequently, they followed the net drift velocity aligned with the wave propagation direction, i.e. the Stokes drift (van den Bremer and Breivik, 2018) (see residual velocity in Figure 1.c).”
“In surface layer, the driving mechanisms are the same that affected low-density particles described herebefore. When spreading throughout the water column due buoyancy-driven settling and vertical turbulent mixing, the particles are exposed to a offshore-directed return current, i.e. the classical undertow, which compensates for the wave-induced surface mass flux (Fig. 1.b-c).”
- line 360. it is not true that the interaction between plastic particles and waves received limited attention, see the reference listed above. Please rephrase
We agree that the interactions between plastic particles and waves have not received limited attention, as demonstrated by the studies listed in the Introduction and throughout the paper including the references raised by the reviewer. However, in this sentence in question, we specifically refer to wave non-linear processes (e.g., wave asymmetry and skewness) in shallow waters, which have received less attention. To better clarify it, we will add “in shallow waters” after “wave non-linear processes (e.g. Martins et al., 2020)”.
10 a section on the limitation of the approach should be added
The main limitations of our approach and potential future improvements were already included in the discussion section, lines 410-422. Nevertheless, we will try to complete this section including some points raised in general comments to ensure that the limitations (and future developments) are more explicitly detailed.
Citation: https://doi.org/10.5194/gmd-2024-100-AC3
-
AC3: 'Reply on RC2', Isabel Jalon-Rojas, 23 Sep 2024
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 315 | 87 | 147 | 549 | 7 | 11 |
- HTML: 315
- PDF: 87
- XML: 147
- Total: 549
- BibTeX: 7
- EndNote: 11
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1