Microtidal river mouths are dynamic environments that evolve as a consequence of many forcing actions. Under the hydrodynamic viewpoint, river currents, sea waves and tides strongly interact, and their interplay determines specific sediment transport and morphological patterns. Beyond literature evidence, information comes from field observations made at the Misa River study site, a microtidal river along the Adriatic Sea (Italy), object of a long-going monitoring. The river runs for 48 km in a watershed of 383 km2, providing a discharge of about 400 m3/s for return periods of 100 years. The overall hydrodynamics, sediment transport and morphological evolution at the estuary are analyzed with particular attention to specific issues like: the generation of vortical flows at the river mouth, the influence of various wave modes (infragravity to tidal) propagating upriver, the role of sediment flocculation, the generation and evolution of bed features (river-mouth bars and longitudinal nearshore bars). Numerical simulations are also used to clarify specific mechanisms of interest.
- estuarine dynamics
- river-sea interaction
- river current
- infragravity waves
Estuaries are dynamic and complex environments. Estuarine hydrodynamics are the product of nonlinear interactions between freshwater flow, tidal, wave, and wind forcing, and bathymetric and topographic changes [1, 2]. Such processes directly affect sediment transport and the morphodynamics of the estuary, which may lead to the formation of complex morphological patterns like river mouth bars  and submerged sandbars .
Although the influence of classical riverine and marine hydrodynamic forcing mechanisms, such as gravity waves and tides, are typically accounted for to describe the estuarine dynamics, an increasing number of studies is focusing on infragravity (IG hereafter) waves, which are seen to play a non-negligible role in estuary evolution . IG waves are a specific type of low-frequency waves with periods between 20 and 30 seconds to 5 minutes, larger than those of sea/swell waves [6, 7]. Typically, IG waves are generated either (i) as long waves bound to short wave groups, or (ii) by a temporal variation of the breakpoint, or (iii) from swell–swell interactions [8, 9, 10].
IG waves are regarded as an important trigger of sandbar generation in the coastal area, especially in the short term, although their role has not been properly understood so far [11, 12]. IG waves are also thought to be of some importance for sediment transport and nearshore morphodynamics, though not a primary forcing . For instance, during storms, breakers mobilize a large quantity of sediment, whose transport is modulated by IG waves .
In general, subtidal bars generating and evolving along sandy coasts are typical of wave-dominated environments as a whole and are the result of complex hydrodynamic and sediment transport patterns taking place within the surf zone [11, 15]. Sandbars exist both close to and far from river mouths, due to their strong link to the marine forcing, although their behavior is largely affected by river current and sediment transport, as well as by artificial structures existing at the mouth [4, 16, 17]. Since sandbars promote wave breaking and energy dissipation, they represent a natural solution for beach protection and are fundamental for coastal stability in both short term (storm scale) and long term (scale of years/decades) timescales [18, 19]. Several studies focused on the definition of evolutionary patterns for sand bars. It is commonly acknowledged that waves break over the bars during sea storms and generate undertow profiles, which lead to sediment being advected seaward and sandbars migrating offshore. Conversely, although with some notable exceptions , onshore bar motion is typically observed under non-breaking waves and relatively mild wave climate, like that occurring in summertime [12, 21, 22, 23].
Evidence shows that IG waves can easily enter a river mouth and propagate upriver for long distances, even during fast river flow conditions. The presence of IG waves in the estuarine regions has been put in relation to edge waves entering the estuary and producing resonance, while IG wave modulation by tides has been observed along the Pescadero River (North California), Ría de Santiuste (Spain), the Albufeira Lagoon (Portugal) and the Misa River (Italy) [1, 24, 25, 26].
River mouth bars are morphological features that generate and evolve due to the direct effect of riverine currents and marine actions (waves, tides). Typically, the sediment transported via river flow deposits out of the mouth, at a distance of about twice the river width, due to flow expansion [27, 28, 29]. Nevertheless, the location where the bar is formed depends on the amount of sediment transport induced by the river discharge and the net residuals of marine actions. Therefore, in a wave-dominated environment, sediment deposits can also be present inside the river mouth and lead to mouth bars. Such bars are extremely relevant for the overall estuarine dynamics, potentially obstructing the river cross-section to the point of contributing to flooding or river overflows. Downriver migration of the bar occurs when the river discharge dominates the estuarine dynamics and the increased flow velocity on the bar crest is able to erode the sediment on the top, which then deposits seaward because of flow deceleration . On the other hand, weak river flows associated with comparatively more intense sea actions induce sediment accumulation and upriver migration of the bar.
Due to higher concentration values with respect to both upriver and river-mouth regions, a “sediment trap” can be generated in the upper estuary due to both suspended matter and residual circulation. A clay content of >5–10% can make fine-grained sediments behave in a cohesive manner [31, 32]. The sediment grains, and often “sticky” natural organic matter that is present, cause sediment particles to cease acting independently. The behavior is typically characterized by flocculation (i.e., individual sediment particles interact and “bond”), a process whereby cohesive and fine-grained mixed sediment particles have the potential to aggregate into flocs [33, 34], which contribute to the formation of near-bed suspension layers in estuaries and may alter water column turbulent mixing, rheology, and sedimentary particle residence (retention) times.
Experiments have shown that when the fine fraction and the larger non-cohesive sediments coexist as a single mixture in an estuary , there exists the potential for the two fractions to combine and exhibit some degree of interactive flocculation [35, 36]. Further, cohesive sediments that are mixed into a predominately cohesionless sandy region can create a ‘cage-like’ structure that can fully encompass the sand grains, thereby trapping the sand within a clay floc envelope . The degree of cohesion between the various sediment fractions tends to increase with the content of fine clay minerals within the sediment, especially for clay contents larger than 5–10% [37, 38, 39, 40], with the biological activity playing an important role in the flocculation of mixed suspensions . In terms of properties, the floc size ranges from microns to centimeters and their effective density generally decreases with size [42, 43, 44], while their settling velocity follows a size-dependent Stokes law settling relationship [45, 46]. Due to their fast settling velocities [47, 48], macroflocs tend to have the most influence on the mass settling flux .
To investigate the interplay between the characteristic processes of microtidal estuarine environments and their effects in the nearby river reach and coastal area, field observations have been performed along the Misa River (MR hereafter) and its estuary, with dedicated field campaigns [50, 51] and long-term measurements collected using both on-site instrumentation and remote sensors [4, 26, 52]. The present work illustrates a comprehensive overview of recent observations and analyses performed to achieve a better understanding of estuarine dynamics at a microtidal river mouth, from the hydrological, morphodynamic, and sedimentological viewpoints.
The chapter is divided as follows. Section 2 describes the field campaigns and the instrumentation deployed along the MR and at the estuary, while the main results are reported in Section 3. An overall discussion, followed by some concluding remarks, is presented in Section 4.
2. Materials and methods
Both short-term and long-term monitoring was carried out at the MR estuary. Specifically, two field campaigns were carried out in September 2013 and January 2014 in the most downstream part of the MR and within the nearby coastal area [50, 51]. Additionally, long-term measurements are being collected since 2015, with a set of instrumentation installed in river and sea, as well as using remote sensors. Numerical simulations are also used to support both experimental findings and speculations, and to better understand the complex dynamics at the MR estuary.
2.1 Study site
The MR runs from the Apennine Mountains (central Italy) to Senigallia (Marche Region), where it flows into the Middle Adriatic Sea after about 48 km. The MR is characterized by a 383 km2 watershed and its flow rate is around 400, 450, and 600 m3/s for return periods of 100, 200, and 500 years, respectively. While a low-flow regime is promoted by a relatively small amount of precipitation in the MR watershed during the summer, significantly higher flow regime characterizes the wintertime [50, 51].
An illustration of the coastal and offshore area is provided in Figure 1a, with a close-up view of the nearshore and estuarine regions shown in Figure 1b. The most downstream portion of the MR is also shown in Figure 1c. The final reach features a heavily engineered river mouth, characterized by cement walls in place of classical riverbanks, which allows one to easily collect and analyze the hydro-morphodynamic data of the microtidal environment. The beach located north of the estuary and harbor is protected by emerged rubble-mound breakwaters, while the southern natural beach is protected by multiple alignments of longshore sandbars (Figure 1b).
From a hydrological viewpoint, the MR estuary is classified as a salt-wedge estuary, where the water column is split into two parts: the outgoing river current flows in the upper portion, while the entering seawater flows in the lower portion . The salt-wedge behavior, i.e. the simultaneous existence of saltwater and freshwater within the water column extends some kilometers seaward, as confirmed by salinity measurements collected up to 2 km from the estuary .
Soil samples were collected within the lower MR reach and confirmed the presence of clay, characterized by high percentages of montmorillonite minerals. Such fine-grained clay sediments promote formation of flocs, i.e. aggregations of individual clay and other fine particles, organic matter, micro-algae and bacteria, and other reactive constituents. Based on both salinity measurements and the sampled clayey material at the riverbed, it was speculated the existence of an upriver flocculation zone (between 2 km and 700 m from the mouth), with salinity S < 10 psu, and a downstream deposition zone (between 700 m from the mouth and the mouth itself), with S > 10 psu [50, 51].
The MR is characterized by large sediment transport rates despite the moderate flow rate, similar to many rivers originating within the Apennine Mountains. Such intense sediment outflow supplies a large amount of material to the coastal areas around the estuary, this being true especially for the natural beach located south of the estuary, which presents an array of longshore sandbars with a long-term morphological evolution [4, 16].
2.2 Short-term and long-term measurements
Within the framework of two international projects funded by the Office of Naval Research Global (UK), named EsCoSed (“EStuarine COhesive SEDiments”) and MORSE (“Modeling and Observation of River-Sea Exchanges at a microtidal estuary”), a series of activities have been planned in the last decade within the microtidal estuary of the MR.
2.2.1 The EsCoSed project
During the EsCoSed project, two field campaigns were carried out in September 2013 and January 2014 with the purpose of investigating, respectively, the summertime regime, mainly characterized by low-flow conditions, and the wintertime regime, where low-flow conditions alternate to high-flow conditions [50, 51]. Observations of meteorology, hydrodynamics and morphodynamics were performed using instrumentation deployed for some days during both investigated periods. Furthermore, in the wintertime experiment, water and sediment were sampled within the river and estuary, as well as in the sediment plume generated during high-flow/stormy conditions.
The hydro-morphodynamic parameters were recorded using small quadpods (an overall height of ∼1 m and a roughly square base of ∼1 m2) both in the final MR reach and in the nearby sea, up to 7-m depth. The quadpods were specifically devised to accommodate a set of instruments for the measurement of water velocity profiles along the water column (six velocity profilers), seabed variations (two pencil beams), the turbidity along the lower portion of the water column (two CT probes), and wave characteristics (one Sentinel 1200KHz Acoustic Doppler Current Profiler - ADCP, from Teledyne RDI®, deployed at a depth of ∼7 m). The quadpod locations in the 2013–2014 experiments are shown in Figure 1b as red and blue triangles, respectively indicating quadpods deployed in the river (QR1, QR2, QR3) and sea (QS1, QS2, QS3). Additionally, information on both water surface level at further locations and surface current were collected using, respectively, two tide gauges (TGup, TGdown) and Lagrangian drifters launched into the final reach of the MR.
Finally, a video-monitoring station named “Sena Gallica Speculator” (SGS) was installed at the Senigallia harbor in 2015 and is currently operating. The station is composed of four cameras located on top of a tower and is oriented to encompass the mouth of the MR and a coastal area located between the MR estuary and a pier, called “Rotonda a Mare” and located 500 m south of the estuary .
2.2.2 The MORSE project
First, a Workhorse Sentinel 300 kHz ADCP (Teledyne RDI®) for the measurements of offshore wave conditions (Figure 1e) was installed within the MEDA station, a meteo-marine station located ∼1.5 nm north of Senigallia and 1.5 nm from the coast (purple triangle in Figure 1a). The station is property of the Italian National Research Council (CNR) and is additionally equipped with several sensors for the monitoring of both atmosphere and sea water.
A tide gauge was deployed within a protected area at the entrance of the Senigallia harbor to record tidal excursions and storm surges (green triangle in Figure 1a).
Finally, a stream gauge was installed within the MR for flow-rate measurements, which is fundamental for the understanding of river-sea interactions. The chosen site is located about 1.2 km upriver of the mouth, near the bridge known as “Ponte Garibaldi”, and very close to a pre-existent hydrometer, property of the Civil Protection (Marche Region), installed for the measurement of the river stages. Both the hydrometer and the stream gauge, the location of which is collectively indicated here as RG1 (Figure 1c), are employed to observe upriver wave propagation from the MR mouth. The RG1 hydrometer complements another hydrometer located at “Bettolelle” (here referred as RG2), about 10 km from the estuary (Figure 1c).
The data collected by some of the above-mentioned instruments are represented in Figure 2. These can be profitably used as initial or boundary conditions for the initialization and validation of the numerical modeling of the MR estuary, as illustrated in the following.
2.3 Numerical modeling
The large amount of observed data allowed us to set up a numerical model of the MR estuarine area using the Delft3D software suite [55, 56]. A two-dimensional, depth-averaged model was used to reproduce the hydro-morphodynamics occurring along the final stretch of the MR estuarine channel . Specifically, a coupled WAVE-FLOW simulation was performed to investigate the deformation and displacement of the river mouth bar under the actions of river current and waves. The WAVE model was forced using time series of wave height, period and direction recorded at the offshore MEDA station. The time series of the water level recorded at the tide gauge was instead used for a first validation of the model.
Observations and results coming from both short-term and long-term recordings are reported in the following. The local dynamics observed during the 2013–2014 experiments are first described, with focus on the main seasonal differences in the hydrodynamics and the interplay between river and marine forcing actions in winter, as well as the saltwater-freshwater interaction across the water column (Section 3.1). Insight from a flocculation model supports the observed local dynamics (Section 3.2). The wave entrance in the MR estuary and the upriver propagation of long-wave modes are detailed in Section 3.3. The morphological processes occurring within the lower river reach (Section 3.4) and in the nearshore area south of the MR estuary (Section 3.5) are finally described, with special focus on the bar evolution.
3.1 Local dynamics
The estuarine area of the MR is subjected to sea storms mainly coming from two directions, as typically observed in coastal regions of the Middle Adriatic Sea. Short, steeper waves generated by WNW, N or NNE (Bora) winds typically enter the MR mouth, since their incoming direction is almost perfectly aligned to the estuary orientation. Also relevant in the MR environment are the ESE-approaching waves, induced by Sirocco winds. Such waves are significantly angled and cannot easily enter the MR. Consequently, Sirocco waves are reflected by the river walls and strongly affect the morphology around the estuary, thus impacting on the evolution of nearby sandbars .
Significant differences exist between summertime and wintertime conditions in terms of wind and rainfall, these directly affecting the wave action and the river current interplaying at the MR estuary. Specifically, the wind blowing during the investigated periods promoted the generation of waves of different heights, which mainly depended on the wind direction rather than on its velocity. The wind direction was frequently changing during the summertime experiment, whereas two intense storms were observed in wintertime and were characterized by almost constant wind directions.
During mild/quiescent conditions in both summer and winter (i.e. time intervals before and after storms), an upriver flow propagation was observed very close to the MR estuary (at QR3, i.e. 290 m upriver), with farther inland locations (QR2, i.e. around 400 m upriver) presenting a significant tidal modulation of the water column. An enhanced salt-wedge behavior was also observed during the wintertime experiment . However, some differences arise between quiescent conditions occurring in summer and winter. Specifically, the marine action in summertime was comparable to the river forcing, as confirmed by the observed tracks of surface drifters deployed within the final reach of the MR. Recorded surface speeds in summer are generally smaller than those recorded in winter, with the surface flow being slowed down and sometimes reversed due to upriver-propagating waves and tide. Conversely, drifter deployed before and after the two winter storms showed an increasing surface velocity while moving downstream. Such behavior is further supported by recorded velocity profiles across the water column at QR2 and QR3 locations, which followed a marked salt-wedge pattern, additionally modulated by the tidal motion .
Hydraulic data recorded during the January 2014 experiment are illustrated in Figure 3. The stage at RG2 and the mean precipitation within the MR watershed are shown in Figure 3a. The two floods are here highlighted by the stage peaks occurred in the mornings of 25 and 28 January, almost simultaneous to high-tide conditions recorded at Ancona harbor (Figure 3b) and to the sea storms recorded at QS3 (significant height and peak periods are illustrated in Figure 3c). The instantaneous velocity magnitude recorded along the water column at QR2 are plotted in Figure 3d.
The increased river outflow during the January 2014 storms generated a large river plume with a considerable amount of sediment extending up to around 1.3 km offshore of the mouth. The occurrence of strong river outflow was coupled with a negligible modulation provided by tidal oscillation on the velocity distribution along the water column, especially at a relatively far distance from the mouth (e.g., at QR2). Additionally, during the winter storms a large sediment deposition was observed closer to the mouth (i.e., near the bend, at QR3), due to the convergence of hydrodynamic fluxes and suspended sediments from both river and sea. This also suggests the existence of a turbidity maximum zone (TMZ), typically observed in macro- and meso-tidal estuaries, but rarely in microtidal environments [58, 59, 60]. Although TMZs were observed during both winter storms, their vertical structure was different and depended on the energetic nature of the storms. The storm occurred during 28–29 January 2014 was much less intense than that occurred during 24–25 January 2014, especially in terms of incoming wave height and outgoing river flow. A relatively larger turbidity was observed during the smaller storm .
While the summer period is supposed to be characterized by a flocculation region at distances larger than 600 m upstream of the MR mouth and by large sediment deposition throughout the final river reach, wintertime stormy conditions enhanced the sediment transport and provided different morphological patterns due to the strong interplay between incoming and outgoing fluxes: 1) a relevant erosion upstream of QR3, where the river discharge dominates on the marine forcing; 2) an erosion/deposition pattern downstream of QR3 and at the mouth; 3) a modification of the rheological properties of the soil both at the MR mouth and within the river-plume area [50, 51].
3.2 Floc dynamics
To investigate the potential depositional effects within the Misa system, a rigorously proven flocculation model (FM) suite was used [47, 48, 61, 62]. The FM is built on a series of empirical-derived algorithms that can quantitatively assess the relative settling and mass flux dynamics for both pure mud flocs and floc populations derived from a range of mud-sand mixtures. It is based on flocculation concept of macroflocs (D > 160 μm) and microflocs (D < 160 μm), as outlined in the classic ‘order of aggregation’ [63, 64]. FM outputs include macrofloc/microfloc settling velocity, ratio of floc mass between the two size fractions and the total mass settling flux (MSF).
The MR estuary has been subdivided into three regions with different representative sediment compositions, two located 600 m and 100 m upriver of the estuary, one located 400 m off the estuary. A nominal representative mud:sand (M:S) composition of 100 M:0S is chosen at the inland site, 50 M:50S equal mud/sand mixture at the mid-zone, and pure sand (0M:100S) in the seaward region, based on the collected field data [50, 51]. The level of flocculation primarily depends upon the combined effects of Suspended Sediment Concentration (SSC) and turbulent mixing. Hence, the FM indicated a fast macrofloc settling (3.5 mm/s at shear stress of 0.35 Pa) and a relatively quick microfloc settling at the pure cohesive inland generated during small storm conditions (28–29 January 2014). Further, both the lower river and seaward zones were much less turbid. In the lower river, the less cohesive equally mixed sediment composition produced slower macroflocs and quicker microflocs than within the inland region, with an equal apportioning of floc mass between macro- and microfloc fractions. The seaward (pure-inert/non-flocculating) sediment settled fast at ∼6.8 mm/s.
Throughout the Misa system, the MSF was only ∼40% of that occurring during small storm conditions. Quiescent conditions (27 January 2014) saw a further 20–25% slowing in the floc settling velocities, with SSC being only 30–40% of that found during the small storm at each site. Specifically, the FM indicated that quiescent conditions favored smaller microfloc fraction dynamics, a much tighter spatial grouping and significantly smaller settling flux values.
3.3 Upriver propagation of long waves
IG waves are long-wave modes that are seen to easily propagate upriver in meso- and macro-tidal estuaries [1, 24, 25], while their impact on microtidal estuaries has been rarely investigated. Upriver propagation of IG waves were observed in the microtidal environment of the MR during the January 2014 experiment (Figure 4a–c) . The dominant generation mechanism of IG waves in this case is probably the bound wave mechanism , as suggested by a normalized slope parameter calculated off the mouth
where Hz is the wave angular frequency during the storm climax,
It was observed that, while the IG contribution in the open sea was significantly smaller than that provided by sea-swell waves, it gained a more relevant role within the final reach of the MR. This is exemplified in Figure 4c and d, which shows the time evolution of band-specific significant wave heights during the storm event occurred on 25–26 January 2014, at two locations representative of offshore (QS2) and riverine (QR2) areas. While the total wave height in the sea was mainly due to sea-swell throughout the storm (Figure 4c), the IG contribution was seen to be much more important in the MR, especially during the storm climax (Figure 4d). This is due to the intense wave breaking affecting the shortest wave components at the mouth, and further enhanced by strong opposing currents. Inspection of wave energy levels revealed a large energy decay across a wide region off the MR mouth (between −400 m and 0 m, at (2–5) m depths, as shown in Figure 4a). This decay especially involved the short-wave components, while the normalized IG component seemed to peak just past the MR mouth and values larger than the reference offshore value characterize the region between 0 and 200 m upriver (Figure 4b). Such behavior is connected to the large damping and deviation of storm-driven sea currents promoted by the river discharge , not enough to block IG components. Hence, the MR estuarine area removes higher-frequency waves and retains low-frequency energy, thus operating as a low-pass filter. This process has already been observed in field experiments carried out in other environments, like microtidal estuaries or energetic oceanic inlets [24, 25, 65].
Although IG waves were observed to propagate upriver for hundreds of meters and tide propagate upriver for kilometers, the interplay between such low-frequency modes is of importance, with the tide controlling the upriver propagation of IG components . In addition, despite the low tidal range at the MR estuary (less than 0.6 m during the January 2014 experiment), tidal oscillations reach distances larger than 1.2 km within the MR (e.g., at RG1), although tidal currents are low and their effect on the river current is negligible. While low-flow river conditions imposes very small resistance to the tide upriver propagation, the tidal forcing was easily suppressed by river floods and high-flow conditions [26, 51]. Furthermore, the tidal effect is important just off the estuary, in connection to a persistence of wave-induced setup during the storm decay.
On the other hand, no tidal influence was observed at RG2, as confirmed by more recent observations performed in 2019 (Figure 5a–c) . These observations also confirmed the presence of tide at RG1, as well as the upriver propagation of additional low-frequency modes. Specifically, the stream-gauge recordings of some flood events occurred in 2019 allowed for the observation of modes not detected by the hydrometer, due to the difference in the sampling rate of the two instruments (2′ for the stream gauge versus 30′ for the hydrometer). The spectral content of two storm events recorded by the stream gauge and lasting a bit more than one day each (12–13 and 15–16 November 2019) has been analyzed, in conjunction with the spectral content of the signals recorded by both hydrometers (RG1 and RG2) and the tide gauge (Figure 5d). A long-time range has been used for hydrometer and tide signals, i.e. between 10 September and 17 December 2019. The hydrometers at RG1 and RG2 showed a similar spectral pattern, especially for frequencies Hz, while significant differences occurred for Hz, due to diurnal (∼25.6 h) and semi-diurnal (∼12.8 h) tidal constituents observed at RG1, but not at RG2.
Furthermore, the analysis of the stream-gauge signals during the recorded events at RG1 does not show relevant peaks referring to tidal constituents, due to the reduced event duration. However, local peaks exist at ∼1.42 hours during both events, which can be probably ascribed to other long-wave modes generating in the Adriatic Sea. Enclosed and semi-enclosed basins, like the Adriatic Sea, show low-frequency oscillations like the seiche motion, which can be described by the natural period at a specific mode
Selecting the mode-2 oscillation (
3.4 River mouth bar
The dynamics of river mouth bars is strongly correlated with the mutual interaction of the river discharge and wave actions. The long-term monitoring of the MR estuary allowed us to correlate the behavior of the emerged mouth bar with the river and sea forcing. Hourly videos recorded by the SGS video-monitoring station were post-processed to create a 10-minute time-exposure image (timex). For each hour, the four timex images were stabilized, geo-rectified and merged to obtain the plan-view images. The ortho-rectified timex images, from 2017 to 2019, were analyzed using a semi-automatic procedure1 to detect the presence of the emerged bar and its geometric features. The evolution of the area and center of mass of the bar was correlated with the sea state (wave height and direction) time-series recorded by the offshore ADCP, the tidal levels, as well as flow rate and water level recorded along the MR.
The results showed a net downriver migration due to relevant flood events (Figure 6e–h) and a persistence of sediment accumulation during periods of weak river action. Moreover, the accurate monitoring made it possible to observe a slow upriver bar migration under wave action (Figure 6a–d). Numerical modeling of the MR estuary, performed by Delft3D software, also highlighted the link between river discharge and downriver bar migration, and between wave action and upriver bar migration .
3.5 Coastal impact
The evolution of submerged sand bars may well be crucially altered (intentionally or not) by the presence of man-made structures, like breakwaters or jetties, the latter of which are commonly present in riverine environments [66, 67]. Melito et al.  discussed sandbar dynamics and their correlations with incident wave climate and morphological constraints at a portion of sheltered beach south of the MR estuary (Figures 1b and 7). The beach object of the study is delimited, at its northern side, by the southern concrete jetty delimiting the final reach of the MR. The jetty provides a partial sheltering from wave attack coming from northern quadrants (mainly originated by Bora storms), but leaves the coastline exposed to incident waves from eastern directions.
The submerged beach is characterized by an array of three shallow bars, whose displacement in response to seasonal climate and storm events is monitored since 2015 with the aid of remote sensing products from the SGS station. The bars at the sheltered beach show a response dominated by seasonal oscillations in wave climate, presenting occasionally consistent onshore displacements in milder climates typical of summer months, and offshore migration in winter months. The bar array is generally poorly responsive to single storms; a circumstance shared with other portions of unprotected coastline far from the influence of the river jetty. The overall behavior of the bar structure can be therefore assimilated to a pattern of bars oscillating around a more or less well-defined point of equilibrium (OPE pattern) .
Two short storm events from SSE occurred in February 2016 and October 2018, however, imposed a remarkable change in this established pattern by generating storm-scale displacements in the order of 30–50 m to all submerged bars. Migrations of such entity and with such short response times are not repeated anytime during the investigation period, spanning from 2015 to 2019, even during much more intense wave attacks. This exceptional event is likely connected to the peculiar interaction between incoming waves from eastern directions and the presence of the concrete jetty. Storm waves approach almost normally to the river jetty, generating enhanced reflection and intensified return currents, ultimately leading to enhanced offshore displacements of bars.
The eventuality of different beach response to storm waves with different incidence was explored with a campaign of numerical simulations run with FUNWAVE  and focused on the interaction of the man-made landmark with the dominant wave field at the MR estuary . Two simulations were devised to represent wave attacks from a typical Bora storm (NNE direction) and a typical Sirocco storm (ESE). The two prototypical storms used as wave input are rather different in terms of wave period: while the Bora storm is characterized by steeper waves (
30-minute-averaged velocity fields and relative wave heights in the nearshore area delimited by the MR jetty are presented in Figure 7 for the Bora (NNE) storm and the Sirocco (ESE) storm. The FUNWAVE model predicts longshore currents on top of the bar in both cases; however, while the current is directed towards the jetty during the Bora storm (Figure 7c), Sirocco waves force a stronger current, directed away from the jetty, and escaping the nearshore giving birth to a marked circulation cell (Figure 7a). Stronger currents, along with higher relative wave heights over the bars for the Sirocco (ESE) event (Figure 7b) in comparison to those modeled during the Bora (NNE) event (Figure 7d), can thus be linked to a greater potential for sediment stirring and motion by Sirocco waves and, ultimately, bar migration during ESE storms.
4. Discussion and concluding remarks
The overall dynamics of the salt-wedge estuary of the MR have been observed exploiting short-term measurements in both river and sea during September 2013 and January 2014, as well as long-term monitoring at different and farther locations. The short-term experiments revealed that the marine forcing is able at propagating upriver for long distances, during sea storms but also during relatively mild/quiescent conditions.
In detail, the large waves observed during both big and small storms propagated upriver for some hundreds of meters, also suggesting the existence of a TMZ developing at different locations within the MR, similar to what observed in meso- and macro-tidal estuaries [57, 59]. Further, results of the FM in the final reach of the MR suggest a fast macrofloc settling and high flocculation occurring about 600 m upriver of the mouth during the small storm, and an efficient flocculation throughout the range of turbulent stimulation, while both lower-river and seaward zones were characterized by a much smaller MSF compared to that estimated inland . The impact of waves in the estuarine area is also suggested by the relevant seabed variations obtained comparing the bathymetric surveys of September 2013 and January 2014, which show significant erosion/deposition patterns just off the jetty (bed variations of ±1.5 m), as well as important changes in the final river reach (erosion up to 1 m) .
However, the wave forcing is not the only action that pushes upriver marine waters along the MR. Evidence of the upriver propagation of marine fluxes was given in the summertime experiment (September 2013), when the salinity level in the MR, at about 1.8 km from the mouth, was larger than zero and suggested the existence of a flocculation zone at about (1.8–0.6) km from the mouth. Further, a floc deposition was supposed to occur in the final 600 m of the MR, as also demonstrated by the seabed increase and sediment deposition occurred in the 2013 summer, between May and September [50, 51]. In addition, more recent observations confirmed that the marine forcing can propagate upriver for kilometers, although not as IG waves, which mainly affect the final reach of the MR like the sea-swell waves do, as observed during the January 2014 storms . Differently, very low-frequency waves like tides or seiches (
In terms of sediment transport and bed morphology, river mouth bars and coastal sandbars are directly linked to the forcing actions existing in the investigated microtidal area. The interplay among such actions is summarized by the evolution of a river mouth bar. While significant flood events promote a net downriver migration, a slow upriver bar migration is observed under wave action.
Linked to the main processes occurring within river and estuary are also the dynamics characterizing the coastal area just south of the MR jetty, where submerged sandbars evolve while protecting the natural beach. While the sediment supply from the MR significantly affects the sediment transport in such area, the main incoming direction of sea storms is thought to largely affect the SSC just off the estuary. Numerical hydrodynamic simulations suggested that during NNE storms a recirculation cell is generated south of the jetty and over the sandbars, while the velocity field at the estuary is relatively small . Hence, the river-induced plume propagating far offshore from the estuary seems to be not significantly affected by NNE waves, whereas a different velocity field seems to generate during ESE storms, providing a seaward-directed fast flow in correspondence of the jetty. Hence, on the one side, the MR plume seems to be enhanced by the wave-induced cross-shore current. On the other side, the sandbar array is supposed to evolve and migrate, due to a larger relative height over the bars themselves and possibly related sediment mobilization.
Many of the above-described aspects, especially the main processes related to the evolution of river mouth bars and river plume, as well as the accurate spatial analysis of the sea-forcing propagation within the MR during mild and stormy conditions, are currently under investigation and will be detailed in future contributions.
The financial support from the MORSE Project (Office of Naval Research Global - UK, Research Grant Number N62909-17-1-2148) and the FUNBREAK Project (MIUR PRIN 2017 - Italy, Grant Number 20172B7MY9) is gratefully acknowledged. AJM’s contribution towards this research was partly supported by the US National Science Foundation under grants OCE-1736668 and OCE-1924532, and HR Wallingford company research FineScale project (ACK3013_62). The authors would like to thank all colleagues who made significant contributions during the planning and execution of the field experiments, as well as for the following activities, including Edward F. Braithwaite III, Sara Corvaro, Giovanna Darvini, Michael Fuller, Kevin Lois, Carlo Lorenzoni, Alessandro Mancinelli, Pierluigi Penna, Allen Reed, Aniello Russo, Alex Sheremet, Luciano Soldini, Tracy Staples, Gianluca Zitti. The following authorities and companies are also acknowledged: the Municipality of Senigallia, the Capitaneria di Porto of Senigallia and of Ancona, MARIDIPART La Spezia and MARIFARI Venezia, GESTIPORT (Senigallia), Club Nautico (Senigallia), NOTA srl (Senigallia), Carmar Sub (Ancona), Sena Gallica (Senigallia), METIS S.R.L. (Senigallia). The data used in this paper are available at this link.
Conflict of interest
The authors declare no conflict of interest.