New insights into combined surfzone and estuarine bathing hazards

Stokes, Christopher; Poate, Timothy; Masselink, Gerd; Scott, Tim; Instance, Steve

doi:https://doi.org/10.5194/egusphere-2024-482

Preprints

https://doi.org/10.5194/egusphere-2024-482

Preprints

15 Mar 2024

| 15 Mar 2024

New insights into combined surfzone and estuarine bathing hazards

Christopher Stokes, Timothy Poate, Gerd Masselink, Tim Scott, and Steve Instance

Abstract. Rip currents are the single largest cause of beach safety incidents globally, but where an estuary mouth intersects a beach, additional flows are created that can exceed the speed of a typical rip current, significantly increasing the hazard level for bathers. However, there is a paucity of observations of surfzone currents at estuary mouth beaches, and our understanding and ability to predict how the bathing hazard varies under different wave and tide conditions is therefore limited. Using field observations and process-based XBeach modelling, we demonstrate how surfzone currents can be driven by combinations of estuary discharge and wave-driven rip currents under various combinations of wave and tide forcing. While previous studies have demonstrated the high hazard that rip currents pose, typically during lower stages of the tide, here we demonstrate that an estuary mouth beach can exhibit flows reaching 1.5 m/s – up to 50 % stronger than typical rip current flows – with a high proportion (>60 %) of simulated bathers exiting the surfzone during the upper half of the tidal cycle. The three dimensional ebb shoal delta was found to strongly control surfzone currents by providing a conduit for estuary flows and acting as a nearshore bar system to generate wave-driven ‘river channel rips’. Despite significant spatiotemporal variability in the position of the river channels on the beach face, it was found to be possible to hindcast the timing and severity of past bathing incidents from model simulations, providing a means to forewarn bathers of hazardous flows.

Received: 18 Feb 2024 – Discussion started: 15 Mar 2024

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 preprint. The responsibility to include appropriate place names lies with the authors.

Download & links

Preprint (PDF, 3493 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 (3493 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

26 Nov 2024

New insights into combined surfzone, embayment, and estuarine bathing hazards

Christopher Stokes, Timothy Poate, Gerd Masselink, Tim Scott, and Steve Instance

Nat. Hazards Earth Syst. Sci., 24, 4049–4074, https://doi.org/10.5194/nhess-24-4049-2024,https://doi.org/10.5194/nhess-24-4049-2024, 2024

Short summary

Christopher Stokes, Timothy Poate, Gerd Masselink, Tim Scott, and Steve Instance

Interactive discussion

Status: closed

RC1:
'Comment on egusphere-2024-482', Anonymous Referee #1, 09 May 2024

Review of “New insights into combined surfzone and estuarine bathing hazards” by Stokes et al.
This manuscript examines surfzone and estuarine flows at Crantock Beach with XBeach simulation. Modeled simulations were calibrated with and tested against observations at this field site. The flows were quantified in terms of bather hazard metrics (maximum offshore flows and particle exits) and used to inform a forecast of bather hazard. This manuscript is mostly well written and clearly articulates the gaps it seeks to fill. However, I have a few concerns, which are discussed below.
Hydrodynamics: The paper describes a range of simulations with varying tidal levels and wave conditions, primarily representative of relatively quiescent conditions (Hs < 1.4 m). However, only results for varying wave power and tide are presented here. I have two concerns regarding this: First, the time series shows several larger wave events and more oblique waves than represented in the selection of simulations. Why did the authors select the set of hydrodynamics to simulate? Are they representative of times when swimmers are at the beach? Considering there are many users in HS3 conditions, should larger wave events be considered? Second, wave power is the only wave property presented here (Hs * Tp). Wave direction is only mentioned as ‘unimportant’ during mid-tide but is not otherwise described in the paper. Presumably, wave direction is important if waves are adequately oblique to reduce the energy entering the embayment. Wave direction is also shown to be important for channel rip hazards (e.g., Dusek and Seim, 2013). Is there little dependence on wave direction here due to the embayment geometry? Additionally, previous work on channel rip currents (e.g., Moulton et al., 2017) suggests that wave height is important for rip speeds but does not incorporate wave period. Assuming wave breaking is triggered at a wave height to depth ratio, wave height defines surfzone width. By presenting these values solely as wave power (which combines both period and height), this paper potentially glosses over some of the dynamics relevant for ‘blocking exits’, etc.
Hazard quantities: The paper primarily relies on two metrics to define hazard: maximum offshore flows and percent surfzone exits. Both are helpful metrics to assess bather hazards. Percent surfzone exits represent a free-floating bather's likelihood of being ejected offshore (but does not represent speed). Maximum offshore flows target how fast a bather may be advected offshore and, therefore, the feasibility for a bather to react (swim) from an offshore flow (but does not represent the offshore flow distribution). These metrics do not represent the number of locations with sufficiently strong offshore flows to eject a swimmer offshore. Providing an additional metric could help. For example, this could be represented as the percent of the locations alongshore with Uoff exceeding a threshold value or, possibly, the Uavg and the rms alongshore of Uoff. This similarly ties into the section on how these flows change with morphology. While Uoff,max does not vary strongly, the position and possibly this distribution of these flows change. This could be explored with an additional hazard metric.
Hazard forecasts: While most methods and results are thorough, the forecasting bathing hazard section needs to be better described. I found the definitions of different terms and how the hazards were predicted and allocated difficult to follow, especially since some definitions are different in the figure caption. For example, the use of seemingly redundant terms (risk, hazard) and the overly brief explanation of the hazard scoring.
Limitations: The paper should describe the limitations of these hazard predictions. For example, the modeled velocities are often underestimated, resulting in non-conservative hazard estimates. The surf-beat depth-averaged model cannot represent all rip current types (e.g., flash rips). The findings are highly tuned to this specific estuary and may not represent other combined surfzone and estuary flows.
Figure quality: Figures 9 and similar layouts are very challenging to read. The quivers are barely visible, and the figures are grainy. Consider using plots similar to Figure 8c,f.

Line-by-line:
L126: I expect the SfM DEMS to perform well at GCP locations because those locations were input into the algorithm to resolve camera geometry. Thus, this may not represent the accuracy of the DEM well. Can the DEM accuracy be checked by comparing the regions overlapping with both surveys (intertidal zone)?
L133: How was this RMSE computed? By comparing with?
L143: List drifter dimensions. Is it truly a surface drifter or representative of a depth-avg current?
Figure 5: Define acronyms in Figure.
L264: How do these tuned values compare with previous studies?
L336: Report bias since the sign is important here (i.e., if the flows are over or underpredicted).
Figure 8: Why are there large drifter passes at single points in the inner surf zone? Do drifters stagnate there?
Section 4.4: Is this supposed to be a new section?
Figure 10: Error in the x-label -> HsTp/overbar{HsTp}.
Section 4.6: The introduction claimed that the increase in swimmer rescues in recent years may be due to changes in the river channel. Here, the authors could add a simulation with synthetic bathymetry representing the previous, potentially less hazardous river channel to see how it compares with these morphologies.
L555: The spatial variability described here should be incorporated more often throughout the paper.
L600: NOAA’s rip current hazard forecasts consider hydrodynamic conditions (https://oceanservice.noaa.gov/news/apr21/rip-current-forecast.html). Perhaps specify that this is the only forecast model considering estuarine flow.
Reference:
Dusek, G., and H. Seim. "Rip current intensity estimates from lifeguard observations." Journal of Coastal Research 29.3 (2013): 505-518.
Moulton, Melissa, et al. "Comparison of rip current hazard likelihood forecasts with observed rip current speeds." Weather and Forecasting 32.4 (2017): 1659-1666.

Citation: https://doi.org/10.5194/egusphere-2024-482-RC1
- AC1: 'Reply on RC1', Christopher Stokes, 30 Jun 2024
  
  The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-482/egusphere-2024-482-AC1-supplement.pdf
  
  Citation: https://doi.org/10.5194/egusphere-2024-482-AC1
RC2:
'Comment on egusphere-2024-482', Anonymous Referee #2, 10 May 2024

General comments:
This observational and modeling study explored swimmer hazards in an understudied setting, where estuary mouth flows encounter surfzone currents including bathymetric rip currents and boundary rip currents with large tidal variations. The authors found that river channel morphology can facilitate not only strong estuary flows, but also strong rip current flows when the river channel modifies wave breaking patterns, similar to what occurs in a surfzone bar-channel system. In addition, prior studies of surfzone hazards have typically found hazard to depend on the water level, which modifies wave breaking, with no dependence on the tidal phase (ebbing vs flooding); in contrast, this study found bathing hazard was different during rising versus falling tides when ebbing or flooding estuarine flows were interacting with the surf zone. I found these conclusions to be very interesting, novel, and supported by the analysis. I do have several concerns about (1) the framing of the paper, (2) the forecasting hazards analysis, and (3) the clarity of the text and figures.
(1) The paper emphasizes estuary mouth flows and bathymetric rip currents. Some discussion of headland/boundary rip currents is included but should be expanded given the clear importance of the boundary flows in this system. In addition, the paper lacks discussion of other rip-current types like flash rip currents, which I would expect to be present, as well as embayment rip currents, which would form in the center of an embayment rather than at the boundaries. Given the importance of the flows resulting from the embayment geometry, I wonder if the title and framing of the paper should be adjusted to “Combined surfzone, embayment, and estuarine bathing hazards.”
(2) The forecasting bathing hazard section isn’t clearly described and is lacking important details. In particular, the analysis relies on a look-up table of hazard statistics from a prior study, but the authors don’t provide a summary of this study or it’s applicability to their study. The authors emphasize how this study site and combination of processes is understudied, so it warrants some explanation why a hazard model developed for a different setting would be the right choice here. The relationship between risk, hazard, and exposure isn’t explicitly stated, and there is some redundancy in presenting all of these separately. Proportions of hazard scores are presented, but a forecast skill assessment should be performed.
(3) Prior to publication I think some improvements to the clarity of the text and figures are needed, particularly to emphasize limitations of the study and make sure that key results are discernable in the figures (see Specific comments).

Specific comments:
L38-43: This text suggests that rip-current patterns/behaviors, rather than speeds, classify the hazard level. This is then inconsistent with the next statement that places importance on speeds. I think this section would be clearer if the authors started with a statement that it is expected that a combination of factors, including pattern and speed, influence the hazard.
L46: Here it is concluded that flows in estuary channels may pose “an equal or potential even higher bathing hazard than rip currents” because the speeds are equal or greater than the speed of rip currents. I don’t think this is known.
L75: “embaymentisation ratio (length/depth)” Maybe “headland amplitude” and “embayment width” would be clearer terms than depth/length? Labeling these scales in Figure 2 could be helpful too.
L75-77: When do boundary rip currents occur, in which fast flows are along the headlands, versus a headland circulation, in which the fast flows are in the middle of the embayment? If a headland circulation is occurring some of the time, could this enhance the flows out of the estuary channel?
L82: It may be helpful to label features such as “ebb tide delta” on the figure.
L98: Did you introduce available data on water users?
L133: Give some information on how the echosounder and UAV datasets are merged. Is this the same as the process described later for the model bathymetry?
L136: Flow measurements were at 0.1 m above the seafloor. It would be worth describing why this is expected to be a good representation of swimmer hazard, and if there are times when it might not be.
L144: Are there concerns about when drifters may be poorly tracking the currents, e.g., if they are scraping the seafloor or surfing waves?
L175: can delete “respectively”
L200: “tidal variation was imposed uniformly on all four corners of the model domain” – what does this mean?
L206: Wave directions are mentioned here, but not wave directional spread, which seems to vary tidally in the observational record. How does this influence the results of this study? I don’t think spread was varied in the model runs? In addition, it would be helpful to know if the observed wave spectra are well described by a single peak period and direction, or if a wave systems approach would be more accurate. How does this affect the results? I suggest referencing the observations here to say how the model spans the observational conditions (shown for this year, is this similar for other years?).
L213-223: Does this method of forcing the estuary flows miss any river-estuary interactions, or does it include them because it’s based on a measurement near the estuary mouth? A small discharge is added to “conservatively account for fluvial flow” – can you elaborate and say if the results are sensitive to these choices?
L232: “Each virtual drifter was advected for 20 minutes, or until they had returned to a safe water depth (<0.7 m).” Could you elaborate on these choices and how they affect the results? Why 20 minutes? Why is a safe water depth 0.7 m? Why not keep running the drifters to see what happens next even if they enter safe water?
L246: should this be lowercase u_off?
L328: Mention what may be the cause of these fluctuations (e.g., related to infragravity pulsations, instabilities, or flash rip currents?).
L355: Given the spatial complexity of the flows, a point to point agreement may not be expected. You could plot the modeled maximum flows and flow range within a spatial region around the observational point for a more fair comparison?
L395-404: There is little or no discussion of rip currents here. How does the variation in rip current strength and characteristics with wave power or other factors influence the hazard metrics?
L427-429: The ebb shoal delta acting as a bar rip system is interesting. Is this mentioned elsewhere in the literature, e.g., papers on flows near a small river mouth encountering a surf zone (Kastner et al., Rodriguez et al., 2018).
L510: It seems the authors are following a prior method (Austin et al., 2013), but more explanation here of how these thresholds and scores were developed is needed to understand this section.
L512: Assuming I’m correct that Risk = Hazard x Exposure (which isn’t spelled out here), and Exposure is a measured quantity, it seems redundant to me to present accuracy results for both Risk and Hazard.

Figure comments:
Figure 1: x and y axes are not labeled in the top panel.
Figure 5: The panel b transect where the Gannel estuary enters the beach shows the topography over a long distance, which does not seem necessary, and makes it harder to see the channel. It may be more useful to see a transect across the shoals to show the scales of the estuary channel and other channels connected to the swimmer hazards.
Figure 6: Were other days similar? This could be interesting to look at as a composite, though maybe that would be difficult given the variation in the behavior with tidal range.
Figure 7: Do other time periods look similar? It could be interesting to show a scatter plot comparing measured and modeled flows.
Figure 8: This figure compares observed drift tracks with gridded results of model drift tracks. It may be helpful to also show example model drift tracks (not expected to reproduce the observed tracks, but presumably similar patterns), and gridded observational data, so that there’s a more one to one comparison. Does this figure only show two of the three regimes described in the text (L374-379)? I suggest using colormaps that are colorblind friendly and perceptually uniform (e.g., Thyng et al., 2016).
Figures 9&11&13: Vectors are small and very difficult to see. Colormap for the vectors includes dark blue, which is the same color as the background water color. This figure shows the variation in the wave breaking, water levels, and inundated bathymetry, but is not readable for information about flow patterns.
Figure 10: This is an interesting figure, though a bit difficult to interpret. It’s clear there is a difference between rising tide and falling tide. Showing line plots of Uoff vs the wave factor for three example tidal elevations for rising and falling tides may be helpful as a summary of this figure.
Figure 12: I’m not sure I find the change plots very helpful. It may be more useful to show how the flows differ for the different cases, to show different spatial patterns, but similar magnitudes?
Figure 14: Is proportion the most useful comparison? This doesn’t show anything about the timing of events. Why isn’t a forecast skill assessment shown? I’m not quite sure what to take from this figure.

References:
Kastner, S. E., A. R. Horner-Devine, and J. M. Thomson (2019), A Conceptual Model of a River Plume in the Surf Zone, J. Geophys. Res. Ocean., 124(11), 8060–8078, doi:10.1029/2019JC015510.
Rodriguez, A. R., S. N. Giddings, and N. Kumar (2018), Impacts of Nearshore Wave-Current Interaction on Transport and Mixing of Small-Scale Buoyant Plumes, Geophys. Res. Lett., 45(16), 8379–8389, doi:10.1029/2018GL078328.
Thyng, K. M., Greene, C. A., Hetland, R. D., Zimmerle, H. M., & DiMarco, S. F. (2016). True colors of oceanography. Oceanography, 29(3), 10. link: http://tos.org/oceanography/assets/docs/29-3_thyng.pdf, https://matplotlib.org/cmocean/

Citation: https://doi.org/10.5194/egusphere-2024-482-RC2
- AC2: 'Reply on RC2', Christopher Stokes, 30 Jun 2024
  
  The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-482/egusphere-2024-482-AC2-supplement.pdf
  
  Citation: https://doi.org/10.5194/egusphere-2024-482-AC2

Interactive discussion

Status: closed

RC1:
'Comment on egusphere-2024-482', Anonymous Referee #1, 09 May 2024

Review of “New insights into combined surfzone and estuarine bathing hazards” by Stokes et al.
This manuscript examines surfzone and estuarine flows at Crantock Beach with XBeach simulation. Modeled simulations were calibrated with and tested against observations at this field site. The flows were quantified in terms of bather hazard metrics (maximum offshore flows and particle exits) and used to inform a forecast of bather hazard. This manuscript is mostly well written and clearly articulates the gaps it seeks to fill. However, I have a few concerns, which are discussed below.
Hydrodynamics: The paper describes a range of simulations with varying tidal levels and wave conditions, primarily representative of relatively quiescent conditions (Hs < 1.4 m). However, only results for varying wave power and tide are presented here. I have two concerns regarding this: First, the time series shows several larger wave events and more oblique waves than represented in the selection of simulations. Why did the authors select the set of hydrodynamics to simulate? Are they representative of times when swimmers are at the beach? Considering there are many users in HS3 conditions, should larger wave events be considered? Second, wave power is the only wave property presented here (Hs * Tp). Wave direction is only mentioned as ‘unimportant’ during mid-tide but is not otherwise described in the paper. Presumably, wave direction is important if waves are adequately oblique to reduce the energy entering the embayment. Wave direction is also shown to be important for channel rip hazards (e.g., Dusek and Seim, 2013). Is there little dependence on wave direction here due to the embayment geometry? Additionally, previous work on channel rip currents (e.g., Moulton et al., 2017) suggests that wave height is important for rip speeds but does not incorporate wave period. Assuming wave breaking is triggered at a wave height to depth ratio, wave height defines surfzone width. By presenting these values solely as wave power (which combines both period and height), this paper potentially glosses over some of the dynamics relevant for ‘blocking exits’, etc.
Hazard quantities: The paper primarily relies on two metrics to define hazard: maximum offshore flows and percent surfzone exits. Both are helpful metrics to assess bather hazards. Percent surfzone exits represent a free-floating bather's likelihood of being ejected offshore (but does not represent speed). Maximum offshore flows target how fast a bather may be advected offshore and, therefore, the feasibility for a bather to react (swim) from an offshore flow (but does not represent the offshore flow distribution). These metrics do not represent the number of locations with sufficiently strong offshore flows to eject a swimmer offshore. Providing an additional metric could help. For example, this could be represented as the percent of the locations alongshore with Uoff exceeding a threshold value or, possibly, the Uavg and the rms alongshore of Uoff. This similarly ties into the section on how these flows change with morphology. While Uoff,max does not vary strongly, the position and possibly this distribution of these flows change. This could be explored with an additional hazard metric.
Hazard forecasts: While most methods and results are thorough, the forecasting bathing hazard section needs to be better described. I found the definitions of different terms and how the hazards were predicted and allocated difficult to follow, especially since some definitions are different in the figure caption. For example, the use of seemingly redundant terms (risk, hazard) and the overly brief explanation of the hazard scoring.
Limitations: The paper should describe the limitations of these hazard predictions. For example, the modeled velocities are often underestimated, resulting in non-conservative hazard estimates. The surf-beat depth-averaged model cannot represent all rip current types (e.g., flash rips). The findings are highly tuned to this specific estuary and may not represent other combined surfzone and estuary flows.
Figure quality: Figures 9 and similar layouts are very challenging to read. The quivers are barely visible, and the figures are grainy. Consider using plots similar to Figure 8c,f.

Line-by-line:
L126: I expect the SfM DEMS to perform well at GCP locations because those locations were input into the algorithm to resolve camera geometry. Thus, this may not represent the accuracy of the DEM well. Can the DEM accuracy be checked by comparing the regions overlapping with both surveys (intertidal zone)?
L133: How was this RMSE computed? By comparing with?
L143: List drifter dimensions. Is it truly a surface drifter or representative of a depth-avg current?
Figure 5: Define acronyms in Figure.
L264: How do these tuned values compare with previous studies?
L336: Report bias since the sign is important here (i.e., if the flows are over or underpredicted).
Figure 8: Why are there large drifter passes at single points in the inner surf zone? Do drifters stagnate there?
Section 4.4: Is this supposed to be a new section?
Figure 10: Error in the x-label -> HsTp/overbar{HsTp}.
Section 4.6: The introduction claimed that the increase in swimmer rescues in recent years may be due to changes in the river channel. Here, the authors could add a simulation with synthetic bathymetry representing the previous, potentially less hazardous river channel to see how it compares with these morphologies.
L555: The spatial variability described here should be incorporated more often throughout the paper.
L600: NOAA’s rip current hazard forecasts consider hydrodynamic conditions (https://oceanservice.noaa.gov/news/apr21/rip-current-forecast.html). Perhaps specify that this is the only forecast model considering estuarine flow.
Reference:
Dusek, G., and H. Seim. "Rip current intensity estimates from lifeguard observations." Journal of Coastal Research 29.3 (2013): 505-518.
Moulton, Melissa, et al. "Comparison of rip current hazard likelihood forecasts with observed rip current speeds." Weather and Forecasting 32.4 (2017): 1659-1666.

Citation: https://doi.org/10.5194/egusphere-2024-482-RC1
- AC1: 'Reply on RC1', Christopher Stokes, 30 Jun 2024
  
  The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-482/egusphere-2024-482-AC1-supplement.pdf
  
  Citation: https://doi.org/10.5194/egusphere-2024-482-AC1
RC2:
'Comment on egusphere-2024-482', Anonymous Referee #2, 10 May 2024

General comments:
This observational and modeling study explored swimmer hazards in an understudied setting, where estuary mouth flows encounter surfzone currents including bathymetric rip currents and boundary rip currents with large tidal variations. The authors found that river channel morphology can facilitate not only strong estuary flows, but also strong rip current flows when the river channel modifies wave breaking patterns, similar to what occurs in a surfzone bar-channel system. In addition, prior studies of surfzone hazards have typically found hazard to depend on the water level, which modifies wave breaking, with no dependence on the tidal phase (ebbing vs flooding); in contrast, this study found bathing hazard was different during rising versus falling tides when ebbing or flooding estuarine flows were interacting with the surf zone. I found these conclusions to be very interesting, novel, and supported by the analysis. I do have several concerns about (1) the framing of the paper, (2) the forecasting hazards analysis, and (3) the clarity of the text and figures.
(1) The paper emphasizes estuary mouth flows and bathymetric rip currents. Some discussion of headland/boundary rip currents is included but should be expanded given the clear importance of the boundary flows in this system. In addition, the paper lacks discussion of other rip-current types like flash rip currents, which I would expect to be present, as well as embayment rip currents, which would form in the center of an embayment rather than at the boundaries. Given the importance of the flows resulting from the embayment geometry, I wonder if the title and framing of the paper should be adjusted to “Combined surfzone, embayment, and estuarine bathing hazards.”
(2) The forecasting bathing hazard section isn’t clearly described and is lacking important details. In particular, the analysis relies on a look-up table of hazard statistics from a prior study, but the authors don’t provide a summary of this study or it’s applicability to their study. The authors emphasize how this study site and combination of processes is understudied, so it warrants some explanation why a hazard model developed for a different setting would be the right choice here. The relationship between risk, hazard, and exposure isn’t explicitly stated, and there is some redundancy in presenting all of these separately. Proportions of hazard scores are presented, but a forecast skill assessment should be performed.
(3) Prior to publication I think some improvements to the clarity of the text and figures are needed, particularly to emphasize limitations of the study and make sure that key results are discernable in the figures (see Specific comments).

Specific comments:
L38-43: This text suggests that rip-current patterns/behaviors, rather than speeds, classify the hazard level. This is then inconsistent with the next statement that places importance on speeds. I think this section would be clearer if the authors started with a statement that it is expected that a combination of factors, including pattern and speed, influence the hazard.
L46: Here it is concluded that flows in estuary channels may pose “an equal or potential even higher bathing hazard than rip currents” because the speeds are equal or greater than the speed of rip currents. I don’t think this is known.
L75: “embaymentisation ratio (length/depth)” Maybe “headland amplitude” and “embayment width” would be clearer terms than depth/length? Labeling these scales in Figure 2 could be helpful too.
L75-77: When do boundary rip currents occur, in which fast flows are along the headlands, versus a headland circulation, in which the fast flows are in the middle of the embayment? If a headland circulation is occurring some of the time, could this enhance the flows out of the estuary channel?
L82: It may be helpful to label features such as “ebb tide delta” on the figure.
L98: Did you introduce available data on water users?
L133: Give some information on how the echosounder and UAV datasets are merged. Is this the same as the process described later for the model bathymetry?
L136: Flow measurements were at 0.1 m above the seafloor. It would be worth describing why this is expected to be a good representation of swimmer hazard, and if there are times when it might not be.
L144: Are there concerns about when drifters may be poorly tracking the currents, e.g., if they are scraping the seafloor or surfing waves?
L175: can delete “respectively”
L200: “tidal variation was imposed uniformly on all four corners of the model domain” – what does this mean?
L206: Wave directions are mentioned here, but not wave directional spread, which seems to vary tidally in the observational record. How does this influence the results of this study? I don’t think spread was varied in the model runs? In addition, it would be helpful to know if the observed wave spectra are well described by a single peak period and direction, or if a wave systems approach would be more accurate. How does this affect the results? I suggest referencing the observations here to say how the model spans the observational conditions (shown for this year, is this similar for other years?).
L213-223: Does this method of forcing the estuary flows miss any river-estuary interactions, or does it include them because it’s based on a measurement near the estuary mouth? A small discharge is added to “conservatively account for fluvial flow” – can you elaborate and say if the results are sensitive to these choices?
L232: “Each virtual drifter was advected for 20 minutes, or until they had returned to a safe water depth (<0.7 m).” Could you elaborate on these choices and how they affect the results? Why 20 minutes? Why is a safe water depth 0.7 m? Why not keep running the drifters to see what happens next even if they enter safe water?
L246: should this be lowercase u_off?
L328: Mention what may be the cause of these fluctuations (e.g., related to infragravity pulsations, instabilities, or flash rip currents?).
L355: Given the spatial complexity of the flows, a point to point agreement may not be expected. You could plot the modeled maximum flows and flow range within a spatial region around the observational point for a more fair comparison?
L395-404: There is little or no discussion of rip currents here. How does the variation in rip current strength and characteristics with wave power or other factors influence the hazard metrics?
L427-429: The ebb shoal delta acting as a bar rip system is interesting. Is this mentioned elsewhere in the literature, e.g., papers on flows near a small river mouth encountering a surf zone (Kastner et al., Rodriguez et al., 2018).
L510: It seems the authors are following a prior method (Austin et al., 2013), but more explanation here of how these thresholds and scores were developed is needed to understand this section.
L512: Assuming I’m correct that Risk = Hazard x Exposure (which isn’t spelled out here), and Exposure is a measured quantity, it seems redundant to me to present accuracy results for both Risk and Hazard.

Figure comments:
Figure 1: x and y axes are not labeled in the top panel.
Figure 5: The panel b transect where the Gannel estuary enters the beach shows the topography over a long distance, which does not seem necessary, and makes it harder to see the channel. It may be more useful to see a transect across the shoals to show the scales of the estuary channel and other channels connected to the swimmer hazards.
Figure 6: Were other days similar? This could be interesting to look at as a composite, though maybe that would be difficult given the variation in the behavior with tidal range.
Figure 7: Do other time periods look similar? It could be interesting to show a scatter plot comparing measured and modeled flows.
Figure 8: This figure compares observed drift tracks with gridded results of model drift tracks. It may be helpful to also show example model drift tracks (not expected to reproduce the observed tracks, but presumably similar patterns), and gridded observational data, so that there’s a more one to one comparison. Does this figure only show two of the three regimes described in the text (L374-379)? I suggest using colormaps that are colorblind friendly and perceptually uniform (e.g., Thyng et al., 2016).
Figures 9&11&13: Vectors are small and very difficult to see. Colormap for the vectors includes dark blue, which is the same color as the background water color. This figure shows the variation in the wave breaking, water levels, and inundated bathymetry, but is not readable for information about flow patterns.
Figure 10: This is an interesting figure, though a bit difficult to interpret. It’s clear there is a difference between rising tide and falling tide. Showing line plots of Uoff vs the wave factor for three example tidal elevations for rising and falling tides may be helpful as a summary of this figure.
Figure 12: I’m not sure I find the change plots very helpful. It may be more useful to show how the flows differ for the different cases, to show different spatial patterns, but similar magnitudes?
Figure 14: Is proportion the most useful comparison? This doesn’t show anything about the timing of events. Why isn’t a forecast skill assessment shown? I’m not quite sure what to take from this figure.

References:
Kastner, S. E., A. R. Horner-Devine, and J. M. Thomson (2019), A Conceptual Model of a River Plume in the Surf Zone, J. Geophys. Res. Ocean., 124(11), 8060–8078, doi:10.1029/2019JC015510.
Rodriguez, A. R., S. N. Giddings, and N. Kumar (2018), Impacts of Nearshore Wave-Current Interaction on Transport and Mixing of Small-Scale Buoyant Plumes, Geophys. Res. Lett., 45(16), 8379–8389, doi:10.1029/2018GL078328.
Thyng, K. M., Greene, C. A., Hetland, R. D., Zimmerle, H. M., & DiMarco, S. F. (2016). True colors of oceanography. Oceanography, 29(3), 10. link: http://tos.org/oceanography/assets/docs/29-3_thyng.pdf, https://matplotlib.org/cmocean/

Citation: https://doi.org/10.5194/egusphere-2024-482-RC2
- AC2: 'Reply on RC2', Christopher Stokes, 30 Jun 2024
  
  The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-482/egusphere-2024-482-AC2-supplement.pdf
  
  Citation: https://doi.org/10.5194/egusphere-2024-482-AC2

Peer review completion

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

ED: Reconsider after major revisions (further review by editor and referees) (02 Jul 2024) by Andreas Günther

AR by Christopher Stokes on behalf of the Authors (07 Aug 2024) Author's response Author's tracked changes Manuscript

ED: Referee Nomination & Report Request started (08 Aug 2024) by Andreas Günther

RR by Anonymous Referee #2 (13 Aug 2024)

Suggestions for revision or reasons for rejection

This is the revised version of an interesting observational and modeling study exploring swimmer hazards in an understudied environment where estuary mouth flows encounter surfzone currents in an embayment. The authors addressed all of the minor line-by-line suggestions. My earlier concerns were (1) the framing of the paper, (2) the forecasting hazards analysis, and (3) the clarity of the text and figures.

For (1), I appreciate the mention of boundary rip currents in the intro and the added comment about cellular circulation is great. I still think the abstract (and possibly the title) should mention the embayed setting and boundary rip currents to accurately represent the paper and to highlight the unique setting that contributes to the paper’s interest and novelty. This system has an estuary entering a beach within an embayment, which has modified wave conditions and boundary rip currents not present on an open-coast. The authors note this in the Discussion: “The embaymentisation […] elicit[s] specific flow behaviours (shadow rips, for example) that won’t necessarily occur in the same way at other estuarine surf beaches.” In addition to many mentions of boundary rips in the results, the authors note in the hazard section that other forecasts have not “yet included dynamics from channel rips, boundary rips, and estuary flow,” in the Discussion there is a section titled “Embayment, estuary, and wave controls on surfzone exits,” and the Conclusions state that “The river channels act to constrain both estuarine and wave-driven currents, […] often connecting with boundary and channel rip currents lower on the beach face.” To me, the embayment setting and boundary rips are a significant part of the story of this paper, but if I’m misunderstanding, maybe the text needs to be revised some to reflect this.

For (2), the revised and expanded forecasting section is much more detailed and clear.

For (3), the authors improved the clarity of the figures and discussion of the study limitations. I like the added time series in Figure 10 and the change to Figure 12 showing the signed cross-shore velocity in color is an improvement over the quivers. In Figure 8, the zoom regions could be shown with boxes in a, d.

I have no further comments.

Hide

RR by Anonymous Referee #1 (31 Aug 2024)

ED: Publish subject to minor revisions (review by editor) (03 Sep 2024) by Andreas Günther

AR by Christopher Stokes on behalf of the Authors (09 Sep 2024) Author's response Author's tracked changes Manuscript

ED: Publish as is (11 Sep 2024) by Andreas Günther

AR by Christopher Stokes on behalf of the Authors (12 Sep 2024)

Journal article(s) based on this preprint

26 Nov 2024

New insights into combined surfzone, embayment, and estuarine bathing hazards

Christopher Stokes, Timothy Poate, Gerd Masselink, Tim Scott, and Steve Instance

Nat. Hazards Earth Syst. Sci., 24, 4049–4074, https://doi.org/10.5194/nhess-24-4049-2024,https://doi.org/10.5194/nhess-24-4049-2024, 2024

Short summary

Christopher Stokes, Timothy Poate, Gerd Masselink, Tim Scott, and Steve Instance

Viewed

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

HTML	PDF	XML	Total	BibTeX	EndNote
335	130	31	496	21	27

HTML: 335
PDF: 130
XML: 31
Total: 496
BibTeX: 21
EndNote: 27

Views and downloads (calculated since 15 Mar 2024)

Month	HTML	PDF	XML	Total
Mar 2024	65	20	6	91
Apr 2024	36	12	6	54
May 2024	51	13	7	71
Jun 2024	39	11	4	54
Jul 2024	91	59	6	156
Aug 2024	19	3	1	23
Sep 2024	17	6	0	23
Oct 2024	8	2	1	11
Nov 2024	9	4	0	13

Cumulative views and downloads (calculated since 15 Mar 2024)

Month	HTML	PDF	XML	Total
Mar 2024	65	20	6	91
Apr 2024	36	12	6	54
May 2024	51	13	7	71
Jun 2024	39	11	4	54
Jul 2024	91	59	6	156
Aug 2024	19	3	1	23
Sep 2024	17	6	0	23
Oct 2024	8	2	1	11
Nov 2024	9	4	0	13

Viewed (geographical distribution)

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

Country	#	Views	%

Latest update: 26 Nov 2024

Download

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

Preprint (3493 KB)
Metadata XML

Short summary

Currents at beaches with an estuary mouth have rarely been studied before. Using field measurements and computer modelling, we show that surfzone currents can be driven by both estuary flow and rip currents. We show that an estuary mouth beach can have flows reaching 1.5 m/s and have a high likelihood of taking bathers out of the surfzone. The river channels on the beach direct the flows and even though they change position over time, it was possible to predict when peak hazards would occur.


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