In this study we investigate the variability in
near-bottom turbidity in Alfacs Bay (in the northwestern Mediterranean Sea). The bay is
characterized by a micro-tidal environment and seiching activity
that may lead to flow velocities of more than 50 cm s
Suspended sediment in the water column and subsequent deposition play a critical role in coastal ecosystems and the management of coastal environments. High levels of suspended sediment in the water column have substantial implications for aquatic ecosystems and natural habitats (Ellis et al., 2002), particularly over longer exposure periods (Newcombe and Macdonald, 1991). Furthermore, sediment supplied from rivers potentially transports significant amounts of organic matter, pollutants, and heavy metals that may be deposited on the seabed or even transported offshore (Palanques et al., 2017). The sediment dynamics are relevant in coastal bays and estuaries due to the large amount of sediment delivered by freshwater and the presence of potential fine-sediment trapping zones in these areas. In addition, sediment resuspension significantly contributes to the total nutrient load (Sondergaard et al., 1992) and prevents sunlight penetration (Mehta, 1989). The analysis and prevention of fine sediment within basins and channels also plays an important role in port engineering and is used to examine and monitor siltation processes (e.g. Ghosh et al., 2001; van Maren et al., 2015). Finally, the growth of harmful species, such as dinoflagellate cysts, may be related to significant local resuspension via the mixing of the upper layers, resulting in more homogenous cyst profiles in the sediment (Giannakourou et al., 2005).
In coastal areas, sediment transport is related to the hydrodynamic conditions. On long timescales, advection processes determine the final depositional pattern as a function of the sediment and water current variables (Ogston et al., 2000; Bever et al., 2009). Hydrodynamic processes driven by wind waves (Grifoll et al., 2013; Carlin et al., 2016), tides (Fan et al., 2004; Garel et al., 2009), winds (Sherwood et al., 1994; Hofmann et al., 2011), surface seiches (Jordi et al., 2008), or internal seiches (Shteinman et al., 1997) promote the resuspension, advection, and settling of fine sediment, conditioned by the continental sediment sources. Subsequent resuspension effects due to natural causes also contribute to the reworking and final deposition of the sediment load (Guillén et al., 2006; Grifoll et al., 2014a). Moreover, anthropogenic activities such as trawling, ship propellers, and waves generated by vessels, introduce additional energy into the water system and influence the resuspension, transport, and final sediment deposition, in particular in shallow waters (e.g. Garel et al., 2009; Hofmann et al., 2011).
This study focuses on Alfacs Bay (in the northwestern Mediterranean Sea – in the southern region of the Ebro Delta) which is a micro-tidal estuary. The bay area is intensively exploited by commercial activities, including tourism, fishing, and aquaculture; hence, the ecosystem has a significant economic importance. In the past, the bay has been extensively investigated in terms of its hydrodynamic response (Solé et al., 2009; Llebot et al., 2014; Cerralbo et al., 2015a, 2016, 2018), tidal wave propagation (Cerralbo et al., 2014), biochemical processes (Llebot et al., 2010, 2011), and optical water properties (Ramírez-Pérez et al., 2017). The estuary mainly receives freshwater discharge from the rice fields on the Ebro River. Several episodes of algal blooms have been reported and have been linked to increased nutrient concentrations, which were possibly triggered by resuspension mechanisms. Moreover, the presence of harmful bacteria has been found in bivalves in the region which has negative effects on aquaculture (Loureiro et al., 2009; Roque et al., 2009).
The goal of this study was to improve the knowledge regarding fine-sediment dynamics in coastal bays and provide insight into the controlling factors of sediment resuspension events observed within the Alfacs micro-tidal bay. Using sea-level heights, water currents, and wind speed measurements we investigated the driving mechanisms of the resuspension of fine bottom sediment within the bay. Subsequently, the spatial and temporal interpretation of the resuspension mechanisms were linked to hydrodynamic processes and analysed via the implementation of a coupled wave–current numerical model. This contribution aims to explain resuspension mechanisms which may have a positive impact with respect to management activities (e.g. harmful species resuspension or algal blooms with negative effects on aquaculture activities).
The water circulation in Alfacs Bay has been extensively analysed in previous studies, using observational data sets and numerical models (Camp and Delgado, 1987; Cerralbo et al., 2014, 2015a; Llebot et al., 2014). However, fine-sediment dynamics and their resuspension mechanisms have not been examined to date. Synchronous optical measurements, in combination with velocity and sea-level measurements, have facilitated an opportunity to investigate the resuspension mechanisms in this region. Considering that Alfacs Bay is a micro-tidal estuary, wind or wind waves are candidate mechanisms for the dispersal of fine sediment.
Alfacs Bay, located in the south of the Ebro Delta, is formed by the
prograding southern spit. The semi-enclosed bay is approximately 16 km long
and 4 km wide. The average depth is 4 m and the maximum depth is about 6.5 m
in the middle of the bay (Fig. 1). A central channel, which is 2.5 km long
and 6.5 m deep, connects the bay with the open ocean. Shallow edges that are around
1–2 m deep are found on both sides of the bay. To the north, the bay is surrounded by rice
fields. From April to December, these fields spill around 10 m
The bay is categorized as a salt-wedge estuary (Camp and Delgado, 1987) and has almost stable stratification all year. Tidal ranges during spring tides reach around 0.2 m, and the hydrodynamic fluctuations are controlled by the wind and modulated by the seiche activity during short periods (Cerralbo et al., 2015a). Both wind and salinity gradients due to freshwater discharge dominate the water circulation in the low-frequency band (Solé et al., 2009; Cerralbo et al., 2018). Intense regional winds from the north and northwest and orographic effects result in wind jets in the Ebro River valley (Grifoll et al., 2015, 2016). This offshore wind is characterized by noticeable spatial variability due to the surrounding topography (Cerralbo et al., 2015b). The water column within the bay is usually stratified due to the freshwater discharge, but well-mixed conditions are common during winter as a consequence of the hydrodynamic response to strong wind forcing (Llebot et al., 2014) and occasionally due to seiches (Cerralbo et al., 2015a). During the summer, the contribution of the temperature to stratification may also be substantial (Cerralbo et al., 2015a).
The bulk of the observational data were collected during a 2-month field campaign from July to mid-September 2013. Water currents were measured using two 2 MHz acoustic Doppler current profilers (ADCPs) that were moored in the mouth (Fig. 1 – A1) and inner bay (Fig. 1 – A2); these instruments were configured to record 10 min averaged data from 10 registers per minute with 25 cm vertical cells. Both devices were equipped with an optical backscatter sensor (Campbell Scientific OBS-3), a bottom pressure meter, and a temperature sensor. The instruments were mounted at the sea bottom at a depth of 6.5 m, whereas the sensors were 0.25 m above the seabed. The signals from the OBS instruments were transformed into nephelometric turbidity units (NTU) according to the device calibration report. In the past, a linear relation between the optical signal and the suspended sediment concentration has been observed in the study area (Guillén et al., 2000). The ADCP has a 20 cm blanking zone. Additional sea-level data were obtained using a sea-level gauge mounted in Sant Carles de la Ràpita harbour (Fig. 1) and bottom pressure systems from the ADCPs. Atmospheric data (wind, atmospheric pressure, solar radiation, and humidity) were obtained from a land-based station (M-Sc) located in Sant Carles de la Ràpita, which was mounted 10 m above the ground.
We use the coupled SWAN-ROMS models included in the COAWST system
in order to simulate the hydrodynamics within the bay. The COAWST system
(Warner et al., 2010) consists of several
state-of-the-art numerical models such as ROMS (Regional Ocean Modeling
System) for ocean and coastal circulation and SWAN (Simulating Waves
Nearshore) for surface wind-wave simulation. SWAN is a third-generation
numerical wave model that computes random, short-crested waves in coastal
regions with shallow water and ambient currents (Booij
et al., 1999). It is based on the wave action balance with sources and sinks
and incorporates state-of-the-art formulations of the processes of wave
generation, dissipation, and wave–wave interactions. ROMS is a
three-dimensional circulation model which solves the primitive variables on
a
The implementation of the COAWST system in Alfacs Bay consists of a regular
grid of
The bottom boundary layer was parameterized using the combined wave–current
(Styles and Glenn, 2000) that was adopted in ROMS and SWAN coupling (Warner et al., 2008). The input
parameters for the model are the velocity components and wave
characteristics near the bottom (wave period, wave direction, and wave
orbital direction). For each computational step, an initial assessment of
bed roughness length is estimated as a function of the grain size, ripples,
and sediment transport. Consequently, the pure current bottom stress (
In order to investigate the suspended sediment events within Alfacs Bay, we
used a subset of the total observations recorded at A2: from
2 to 8 August 2013. This is because the subset data selected
include the main hydrodynamic conditions likely to increase the
near-bottom turbidity. Figure 2 shows the time series recorded at A2 in
terms of NTU from the OBS, measured sea-level height (the sea-level
height measured at A1 is also shown), bottom current speed in m s
Time series of the variables measured during the field campaign.
Two typical wind conditions are considered (Fig. 2a, b): sea breeze
and the northwesterly winds (Cerralbo et al.,
2015a). The sea breeze is associated with an increase in wind speed during the
central hours of the day (approximately from 11:00 to 18:00 GMT with a
wind direction within the range of approximately 30 to
180
During the analysis period (on 3 August), a seiche event was also captured.
This seiche event was previously described (hydrodynamically) in
Cerralbo et al. (2015a) revealing an
oscillation (of 1 h periods) in sea level and currents. This oscillation is
characterized by a node (roughly located at A2) where the velocities
reach a maximum, and an antinode (approximately located at A1) where the
amplitude of the sea-level oscillation reaches a maximum (see sea-level height at A1 in
comparison to A2 in Fig. 2c). The homogeneous vertical profile in
velocities measured at A2 is shown in Fig. 3, with velocity peaks in the
order of 0.5 m s
Each panel shows (at the top) the wind speed measured at M-Sc in m s
The near-bottom turbidity shows a fluctuating behaviour with values ranging from almost 0 to over 10 NTU (Fig. 2e). Three distinct events with high turbidity are observed. These events are E1 (covering the period from 08:00 GMT on 3 August to 10:00 GMT on 5 August), E2 (from 03:00 to 12:00 GMT on 6 August), and E3 (from 08:00 GMT on 7 August to 15:00 GMT on 8 August). The maximum turbidity is measured during the E1 event (maximum turbidity of 41.1 NTU). This event also lasts longer than the E2 (maximum turbidity of 4.6 NTU) and E3 (maximum turbidity of 12.1 NTU) events.
The performance of the water circulation model used in this study was examined in terms of sea level, water currents, and temperature/salinity evolution in previous research (Cerralbo et al., 2014). However, in this work we pay attention to the near-bottom velocities due to their relevant role in the sediment resuspension and transport dynamics. Thus, the skill assessment of the near-bottom velocities at A1 and A2 is analysed using a Taylor diagram (Taylor, 2001). This diagram characterizes the similarity between numerical model and observations using their correlation, the root mean square difference (RMSD), and the amplitude of their variations (represented by their standard deviations). The skill of the model improves when the triangles are closer to the observation reference point in the diagram which denotes full agreement between the model and the observations (Fig. 4). In general, the model results showed good agreement with the observations in the prevalent along-shelf direction, with correlations larger than 0.5 and RMSD values below 1. In addition, the water current fluctuations are well represented in the model because the normalized standard deviation is close to 1 at both measuring points.
The bottom stress is obtained from the coupled numerical model implemented
in Alfacs Bay. Figures 5 and 6 show different snapshots of the modelling
results to examine the bottom stress pattern for two components
(i.e. wave-induced and current-induced bottom stress). These results correspond
to the different episodes identified from the previous observational analysis. The
plot scale of the bottom stress is transformed in
Taylor diagram comparing the error metrics between the observations and model results for the near-bottom currents. A1 and A2 correspond to the ADCP locations shown in Fig. 1.
After the seiche activity (second stage of E1), the wind speed increases due to the sea breeze and the current-induced bottom stress (5 August 2013, 08:00 GMT) decreases significantly, particularly in the centre of the bay (Fig. 5b, d, f). The bottom stress distribution shows how the maximum values are obtained near the shoreline (2.2 Pa) due to the contribution of the wave-induced bottom stress. At A2, the combined bottom stress is equal to 0.03 Pa (which is presumably too small to induce resuspension). For this event, the wave field during the sea breeze is shown in Fig. 7; it reveals that the maximum significant wave height (equal to 0.3 m) occurs near the northern and southern shallow edges which is consistent with the maximum wave-induced bottom stress.
Current-induced bottom stress (
The bottom-stress pattern during episode E2 (Fig. 6a, c, e) is similar to the second stage of episode E1. Both wave and current bottom stress (08:00 GMT, 7 August) tend to be low at A2 in comparison to the seiche event. Substantial bottom stress is only observed on the shallow edges of the bay due to the wave action originating from the sea breeze.
Current-induced bottom stress (
During the E3 episode (northwesterly wind; Fig. 6b, d, f), the combined bottom stress (23:00 GMT, 8 August) is dominated by both wave and current action. The southern part of the bay shows that the maximum wave-induced bottom stress is consistent with the wave climate (Fig. 7). Furthermore, the current-induced bottom stress presents non-negligible values within the bay. Focusing on A2, both mechanisms contribute in a similar manner (wave and current bottom stress are 0.09 and 0.06 Pa respectively) to the combined bottom stress.
Wave field for the E2 (sea breeze;
The synchronous time series of the meteo-oceanographic variables and
turbidity shown in Fig. 2, in combination with the bottom stress model provides a
good opportunity to characterize the turbidity peaks measured at A2. During
the first stage of the E1 episode, the bottom current speed responds to the
node–antinode pattern with velocities that increase to 0.4 m s
The turbidity still showed large values after the seiche had already
dissipated and the bottom current had decreased during the second stage of the
E1 event. Typical sea-breeze wind conditions were observed (gentle variation
of wind direction from 30 to 180
Episode E2 is attributed to a sea-breeze mechanism. This event is qualitatively less important in terms of the turbidity measured at A2. The comparison of the sea-breeze events on 4 and 6 August (both have similar wind and bottom current speed but different turbidity values) seems to indicate the relevance of the previous events and the subsequent advection of fine sediment, following the mechanism explained previously. Similar to the second stage of E1, the bottom stress is low (below 0.02 Pa) in the central basin of the bay, indicating that local resuspension is unlikely. Consequently, the turbidity measured at A2 is probably due to advection processes of suspended sediment from the shallowest areas (combined bottom stress more than 0.8 Pa) in the central basin.
Finally, episode E3 corresponds to a strong northwesterly wind event with wind speeds
in excess of 12 m s
E2 and E3 are examples of two mechanisms that may increase the turbidity: wind-driven current and wind waves. In Alfacs Bay, the roles of these mechanisms in sediment resuspension are less clear that the role of seiches, as the former are a function of wind speed without any clear correlation between the wind module and the observed turbidity. The resuspension of fine sediment due to wind and wind waves in shallow environments has been reported in literature (Luettich et al., 1990; Ogston et al., 2000; Guillén et al., 2006; Bever et al., 2011; Grifoll et al., 2014b; Hawley et al., 2014; Martyanov and Ryabchenko, 2016; López et al., 2017). Some of these studies highlight the complexity of the sediment processes due to the temporal and spatial variability of the importance of resuspension mechanisms and the presence of available material to be resuspended. Apparently, this is the case in our observations, because similar wind conditions do not imply the same turbidity measurements. A good example is the sea-breeze wind events on 4, 5 and 6 August during which different turbidity values were observed. As mentioned in the previous section, advective fluxes and the sequence of events may play an important role in the observed water turbidity. Many authors have reported an apparent influence from advective fluxes correlated with the suspended sediment concentration after an initial deposition of fine sediment (Sherwood et al., 1994; Ogston et al., 2000; Guillén et al., 2006; Harris et al., 2008; Bever et al., 2009; Grifoll et al., 2014b). This means that on longer timescales, advection of sediment by currents may redistribute sediment and determine final deposition patterns (Wright and Nittrouer, 1995). This may be the mechanism responsible for the high turbidity observed under relatively low hydrodynamic conditions. For instance, the sea-breeze event on 2 August does not cause high turbidity, in contrast to the event on 5 August (second stage of the E1 event); this may indicate that an energetic event (i.e. seiche) could mobilize sediment, which is then easily resuspended in subsequent events. The lack of proportionality of the resuspension related to hydrodynamics is also found in extended data time series, where divergences are mainly associated with sediment availability on the bottom, among other factors (e.g. in Wiberg et al., 1994 or López et al., 2017). In the case of Alfacs Bay, more extended observations may clarify the relation between wind intensity, wind waves, seiches, and the amount of suspended sediment and fluxes, taking the sequence of energetic events into account.
The sediment distribution in Alfacs Bay (high percentage of silt and clay in the central basin and sand prevalence on the southern, eastern, and western shores) is consistent with the modelling results shown in this study, where larger bottom stress values were obtained in the lateral shallow edges due to the contribution of the wave-induced bottom stress in shallow areas. However, the deposition mechanism is a complex process, composed of an initial settling and a subsequent dispersal, as described in Wright and Nittrouer (1995). Further sediment transport simulations, including those considering sediment classes and erosion and settling effects, would help to investigate the sediment settling dynamics and its final deposition. These processes should include the cohesive nature of the fine sediment or other phenomena, such as armouring or bioturbation, which may modify the physical properties of the sediment layers (van Ledden et al., 2004; Amoudry and Souza, 2011).
The characteristics of the bay, such as the relatively narrow and shallow entrance, favour the trapping of fine sediments, fed either by the freshwater outflow or the exchange between the open sea and the inner bay. The trapping effect of the bay may entail the presence of a thin surface bottom layer of fine sediment that is easily subject to resuspension. This behaviour is typical of shallow and sheltered environments such as lagoons or lakes. According to Luettich et al. (1990) and Hofmann et al. (2011), regular resuspension events in sheltered and shallow water bodies prevent sediment consolidation and the formation of a cohesive sediment layer. This could explain the high turbidity values observed in Alfacs Bay under relatively weak conditions, such as sea-breeze events, which would likely not occur if the sediment was cohesive.
The Sentinel-2 satellites provide imagery which allows for the further
identification of scenarios with resuspension linked to hydrodynamic
forcing. Figure 8 shows the total suspended matter (TSM in mgr L
Total suspended matter (TSM in milligrams per litre) obtained from Sentinel-2
imagery for Alfacs Bay for two different conditions: northwesterly winds
(
Also, the proximity of the Ebro River mouth (15 km at north) may increase the suspended sediment within the bay under specific circumstances. River discharge is the main driver of the Ebro River plume, followed by wind and regional oceanic circulation (Fernández-Nóvoa et al., 2015; Mestres et al., 2003). Analysis of the turbid plume by remote sensing products indicate that more than 70 % of the plume extension was located south of the river mouth, influenced by the regional oceanic circulation (Fernández-Nóvoa et al., 2015). Other external sediment sources may be associated with freshwater discharge from channels, overwash in the bar, flash floods from small creeks, or aeolian transport. A complete study of the suspended sediment dynamics will provide objective information to address the problem of degrading water quality within the bay, and will assess how to make use of natural mechanisms to limit undesirable concentrations of nutrients or pollutants. This applies in particular to harmful algae blooms that are prone to occur in the area under present and future conditions.
The observational set and the wave–current numerical results obtained for Alfacs Bay allowed for a thorough investigation of the resuspension mechanisms of fine sediment. The results indicate evidence of a clear mechanism of resuspension induced by eventual seiche events, which according to the bottom stress patterns may have a relevant spatial variability within the bay consistent with the node–antinode position. The wind and wind wave mechanisms are also responsible for fine-sediment resuspension during energetic wind events, especially in shallow areas of the bay. A similar contribution of current-induced bottom stress and wave-induced bottom stress in the assessment of the combined wave–current bottom stress is found for the northwesterly wind case. In the case of the sea breeze, the wave-induced bottom stress prevails. The lack of proportionality between the turbidity and the hydrodynamics highlights the relevance of the sequence of the events, suggesting an effect of advective sediment fluxes within the bay (from the lateral shallow edges to the middle of the bay) which increase the suspended sediment concentration. The trapping effect of the bay may entail the presence of a thin surface layer of fine sediment that is easily involved in resuspension, neglecting the expected cohesive effects. However, these points deserve further analysis using extended data sets and sediment transport modelling. The exchange of fine sediment between the bay and the open sea is also evident according to remote sensing images. However, these points require further analysis in the same fashion as the previously mentioned factors. As the Ebro Delta is a region with high-anthropogenic pressure, this research may contribute to furthering integrated development plans that consider sustainable aquaculture activities and climate change mitigation in the region.
Data from field campaigns and model outputs are available upon request from the first author.
MG led the research and the writing process. All co-authors contributed equally to this work.
The authors declare that they have no conflict of interest.
This article is part of the special issue “Coastal modelling and uncertainties based on CMEMS products”. It is not associated with a conference.
Data provided by Puertos del Estado and AEMET are greatly appreciated. This work received funding from the EU H2020 program under grant agreement no. 730030 (CEASELESS project). The authors also wish to thank the Secretaria d'Universitats i Recerca del Dpt. d'Economia i Coneixement de la Generalitat de Catalunya (Ref. 2014SGR1253) who supported the research group. Pablo Cerralbo acknowledges financial support from ECOSISTEMA (CTM2017-84275-R INECO/AEI/FEDER, UE). This paper contains modified Copernicus Sentinel data (2017/2018).
This paper was edited by Joanna Staneva and reviewed by three anonymous referees.