Next Article in Journal
Ground Deformation and Seismic Fault Model of the M6.4 Durres (Albania) Nov. 26, 2019 Earthquake, Based on GNSS/INSAR Observations
Next Article in Special Issue
A Computer Aided Approach for River Styles—Inspired Characterization of Large Basins: A Structured Procedure and Support Tools
Previous Article in Journal
Dating the Sedimentary Protolith of the Daldyn Group Quartzite, Anabar Shield, Russia: New Detrital Zircon Constraints
Previous Article in Special Issue
2D Runout Modelling of Hillslope Debris Flows, Based on Well-Documented Events in Switzerland
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Modeling of the Hydro-Morphodynamics of a Distributary Channel of the Po River Delta (Italy) during the Spring 2009 Flood Event

1
Department of Hydrology and Hydrodynamics, Polish Academy of Sciences, Ksiecia Janusza 64, 01-452 Warsaw, Poland
2
Department of Earth and Environmental Sciences, Dalhousie University, Halifax, Nova Scotia B3H 4R2, Canada
*
Author to whom correspondence should be addressed.
Geosciences 2020, 10(6), 209; https://doi.org/10.3390/geosciences10060209
Submission received: 17 April 2020 / Revised: 22 May 2020 / Accepted: 28 May 2020 / Published: 30 May 2020
(This article belongs to the Special Issue Numerical Modeling of Surface Processes)

Abstract

:
One-dimensional (1D) numerical models generally provide reliable results when applied to simulate river hydraulics and morphodynamics upstream of the tidal influence, given the predominantly unidirectional flow conditions. Such models, however, can also be used to reproduce river hydraulics across the fluvial to marine transition zone when specific conditions occur, as during high discharge events, and the results obtained via these simple modeling tools can provide indicative trends that may guide more structured and detailed modeling of a particularly critical area. In this study, the application of a 1D model setup with hydrologic engineering centers river analysis system (HEC-RAS) for simulating the hydro-morphodynamic conditions of a distributary channel of the Po River Delta (Italy) during a flooding event that occurred in Spring 2009 is presented. The channel bathymetry and the grainsize composition was taken from field measurements, while the dimension of the plume offshore the delta was derived from a MODIS image acquired at the peak of the flood. The comparison between the numerical outcomes and the field evidence shows the reliability of the proposed 1D modeling approach in representing the delta dynamics at a large scale, as well as in showing locations where more spatially detailed studies are needed. The code was also able to adequately reproduce the channel hydro-morphodynamics and the sediment data as derived from a core sample taken a few km offshore during the flooding event of April–May 2009. Through a sensitivity analysis, it is also proven that the dimension of the river plume can influence the evolution of the prodelta, while having a rather negligible effect inland, because of the major stresses induced by the high river discharge during the flood event.

1. Introduction

Nearly 10% of the world’s population lives in coastal areas that are less than ten meters above sea level, along river deltas and adjacent coastlines [1]. Such transitional environments are under threat due to relative sea-level rise and increasing anthropogenic pressure [2] and it is thus fundamental to understand how they will respond to external perturbations [3,4]. To reach this goal, and to adequately simulate future delta dynamics resulting from natural processes and human activities, it is fundamental to quantify their present state and the fundamental mechanisms that control their growth. Likewise, reconstructions of the past evolution of river deltas rely on the knowledge of deltas’ morphodynamics that is gained by studying modern systems, as well as by the modeler´s ability to numerically reproduce them using adequate tools [5,6]. At the same time, improved insights on the interactions between river discharge, sediment transport, and tidal conditions are paramount for tackling future challenges such as delta erosion due to sediment trapped in upstream reservoirs, urbanization, or subsidence [7].
In the last decades, commercial and open-source 1D numerical models have been used for simulating fluvial conditions by many researchers and agencies dealing with water resources management and flood control [8]. Generally, these models allow for solving both the equations of water flow and the sediment continuity equation in the streamwise direction. Among others are hydrologic engineering centers river analysis system (HEC-RAS), developed by the US Army Corps of Engineers, Telemac Mascaret 1D by a European consortium based in France, Germany, and the UK, Mike11 by the Danish Hydraulic Institute, and SOBEK (https://www.deltares.nl/en/software/sobek/)1D by Deltares (NL), all of them widely applied and, consequently, well-tested in different conditions. It is worth noting that, while open-source and freeware codes have been optimized for specific applications, commercial software packages have greatly improved to easily integrate geographic information systems (GIS) with data derived from climatic and hydrological models, aiming to allow for a downscaling of the problem via more detailed (2D or even 3D) descriptions.
Albeit the present 1D models have been highly developed and are now able to simulate a variety of hydro-morphological conditions, a rather reduced literature exists about the application of very simplified models for reproducing tidal-driven environments like deltas, estuaries, and lagoons [9,10,11]. Indeed, most frequently, such transitional environments are reproduced using 2D or 3D solvers [5,12,13,14], given their embedded complexity that cannot be adequately schematized by simplistic 1D models [15]. As an example, a 2D model was applied for simulating the evolution of the Po River delta, which has a rather complex geometry with narrow bends and strong gradients in the riverbed profile [16,17,18]. To account for the presence of tidal effects, the authors applied the classical 2D Saint-Venant equations, written in conservation form, looking for a nonstationary solution. These partial differential equations were discretized by a finite-element scheme, aiming to reduce the complexity of the system but, at the same time, to adequately represent the phenomena involved.
Notwithstanding, river hydraulics in the fluvial to marine transition zone can also be simulated by 1D models, on the occurrence of particular conditions, for example during a river flood, when the river flow overwhelms tidal currents, generating one-directional flow. Aiming to provide evidence on that, in this study, a 1D HEC-RAS model was implemented for reproducing the hydro-morphodynamics of a distributary channel of the Po River (Italy), the Pila channel, during a flood event that occurred in the Spring 2009, following previous application of this code to the same river, even if more upstream [19,20]. This application pointed out that 1D models can be helpful in detecting trends and locations where more detailed modeling is necessary, definitely saving computational time and resources with respects to more thorough simulations with more complex codes. After presenting the case study and the data used in the model, the main numerical characteristics are described, and the results of the calibration-validation procedure are shown, highlighting the need to critically design the model setup. The results demonstrate that the model is able indeed to realistically reproduce both the water stages and the bed composition measured along the final part of the river and in the estuary.

2. Materials and Methods

2.1. Case Study

The Po River, the longest watercourse in Italy, flows eastward across northern Italy for about 660 km, from the Alps to the Adriatic Sea (Figure 1). With a drainage area of about 74,000 km2, its basin can be divided into three sectors based on the lithology and the maximum elevation: the Alpine sector of crystalline and carbonate rocks (maximum relief ∼4500 m asl), the Apennine sector, mostly composed of sedimentary rocks with high clay content (maximum relief ∼2000 m asl), and the central alluvial area that includes the Po plain and the Po River delta [17,18,19,20]. The middle and lower parts of the Po River are subjected to high flood hazard and, consequently, have been heavily embanked to protect the intensively anthropized surroundings [21], with important consequences for the environment [22]. Indeed, the present state of the Po River and its distributary channels results from the interaction between natural processes (acting on time scales from seasonal to millennial) and anthropic modifications [3,19].
The annual hydrograph at the gauging station of Pontelagoscuro shows two peaks in discharge, normally in autumn and spring, generated by rainfall and snowmelt, respectively [21]. The mean annual discharge measured at this gauging station during the period 1984–2018 was around 1400–1500 m3/s, with a peak discharge higher than 12,000 m3/s. The total annual sediment and freshwater yield delivered to the Northern Adriatic Sea were about 13 × 109 kg and 40–50 km3, respectively [23,24,25]. Before reaching the Adriatic Sea, the main stem of the Po River forms five major distributary channels within the delta; from north to south, these are the Maistra, Pila, Tolle, Gnocca (or Donzella), and Goro. Among these five, the Pila carries the majority of water and sediment (about 60%), with secondary contributions from the other distributary channels [26,27]. Spring and neap tides in the northern Adriatic Sea have amplitudes of ~86 cm and 20 cm, respectively. At the Pila mouth, the river flow is reversed during flood tide, when the river discharge is less than 1000 m3/s [12].

2.2. The Spring 2009 Flood

The April–May 2009 river flood was triggered by the combination of spring snowmelt and a significant atmospheric perturbation coming from the Atlantic Ocean [28]. In fact, the winter 2008–2009 was particularly wet, allowing for significant accumulation of snow on the Alps (Figure 1b), which was then released starting from mid-April, because of an increment of the temperature in the basin. At the end of April, heavy precipitation occurred in the entire watershed, but mainly concentrated in the western part of the Po River catchment. The co-occurrence of rainfall and snowmelt resulted in increasing the river flow, and the peak discharge recorded at the beginning of May 2009 was significantly over 8000 m3/s at Pontelagoscuro (Figure 2).
The most notable Po River flood event before the one that occurred in Spring 2009 was the flood of October 2000, the second-largest flood of the last 150 years [27]. As pointed out by Tesi et al. [29], these two major flooding events exhibited interesting contrasts and similarities. In detail, the two floods occurred during different seasons and were triggered by different events: exclusively exceptional precipitations caused the October 2000 flood, while a combination of snowmelt and rainfall drove the May 2009 event. In addition, whereas the 2000 Po River flood was a single event that reached a maximum discharge of 9650 m3/s at Pontelagoscuro, the 2009 flood occurred at the end of a long rainy period which led to a series of smaller events exceeding 4000 m3/s prior to the Spring flood. Furthermore, the 2009 river flood occurred during marine fair weather conditions (wave height < 1 m) and relatively small tidal oscillations (∼50 cm) (Figure 2d,e), suggesting possible deposition in the shallow part of the prodelta, as observed for the October 2000 flood [27].

2.3. Numerical Model Setup

To simulate the changes of the water profile and the bed elevation in the period of 10 April–10 May 2009, the well-known code HEC-RAS version 5.0.7 [28] was used, imposing the input data reported in Section 2.3.1. The calibration and validation procedures show that the code reproduces in a reliable manner the water levels measured in an intermediate section.
HEC-RAS (www.hec.usace.army.mil/software/hec-ras) allows the user to perform 1D steady flow, 1D and 2D unsteady flow calculations, sediment transport/mobile bed computations, and water temperature/water quality modeling, along a network of natural and man-made channels. In the present application, 1D unsteady flow calculations were run for computing the sediment transport and the morphological changes due to a flooding event, using the data reported below, and following the calibration and validation procedure described in Section 2.3.2. It is worth to note that, being HEC-RAS a 1D model, the bed shear stress values are cross-section averages, computed as a function of the average depth (total flow area divided by the top width) and the slope of the computed local energy grade line.

2.3.1. Input Data: Channel Geometry, Flow Parameters, and Grainsize Composition

The geometry of the Po River reach downstream of Pontelagoscuro was described by combining 57 cross-sections measured in 2005 during a LiDAR campaign (cross-sections available at geoportale.agenziapo.it) with cross-sections of the delta front and prodelta derived from bathymetric data [30] and from satellite images to reproduce the whole plume up to the location of the sediment core E16 (Figure 1c,d). For the offshore part, three sections spacing ca. 1000 m and having a rather regular shape were assumed to adequately represent the plume. Aiming to allow for a better spatial resolution, all the cross-sections were linearly interpolated imposing a maximum separation of 50 m. Different interpolation distances were tested to find a compromise between the spatial accuracy and the computational effort.
The hydrological boundary conditions were defined through the 30 min averaged discharge at the upstream end of the simulated reach (Pontelagoscuro), as derived from the discharge-stage relationship provided by the AIPO (Agenzia Interregionale per il fiume Po)—Po Water Authority. Following a previous report [26], the river discharge was reduced to 60% to account for the reduced transport along the Pila distributary channel. The downstream boundary conditions were set using the water levels measured at the river mouth and supposing that these levels do not change on the prodelta plume. All the water levels were derived from the AIPO website (www.agenziapo.it/content/monitoraggio-idrografico-0). To adequately reproduce the hydrological variations, unsteady flow simulations were performed, imposing an integration time-step of 30 s, which was chosen after a sensitivity analysis for code stability.
The bed grain size composition was derived from field samples, collected and analyzed by AIPO in the early 2000s, along the main channel of the Po River and the sediment sample in the prodelta (location E16 in Figure 1c,d) [31]. For computing the sediment transport, the Meyer Peter Müller formula was adopted, imposing equilibrium conditions between sediments and water flow at the upstream end. The reliability of the proposed approach was tested in previous applications made on the same river, also using other solvers [20,32,33].

2.3.2. Model Calibration and Validation

The model was calibrated and validated, in terms of bed roughness, for the period 10 April–10 May 2009, using the first two weeks as a calibration period and the rest for the validation. The 30-minute-averaged computed values were compared to the ones measured every 30 min at the gauging station of Polesella, located around 20 km downstream of Pontelagoscuro (Figure 1c).
As visible from Figure 3a, the calibration period covered relatively low water levels, while the flooding event characterized the validation timespan. Despite this difference, a roughness Manning coefficient n of 0.020 s/m1/3 in the main channel provides reliable values for all the water levels observed in April–May 2009, with a slight underestimation of the high water values during the rising limb of both floods (Figure 2b). Such values of roughness agree rather well with literature evidence [1,34].
Being offshore and representing the river plume over the prodelta, for representing the last to km, a very high Manning value of 0.10 s/m1/3 was assumed. Between the river mouth and the prodelta (around 500 m), the roughness values were linearly interpolated, as made also for the bed geometry.

3. Results

3.1. Longitudinal Profile

The numerical model shows that a drawdown M2 profile characterizes the water surface for all discharges during the flood and that water profiles with gentler slopes occur during high tides (Figure 4). Regardless of the tide, the drawdown profile starts from the Pila mouth and extends for around 15–20 km in the case of high discharges, while during low flow the influence is reduced to only 8–10 km, and the tide is not changing the water surface in a significant manner.

3.2. Cross-Section Profiles

3.2.1. Gauging Station at Pila

Generally, at the gauging station of Pila (Figure 1a), the tidal oscillations have a very low influence on the water levels, as reported in Figure 5a. Indeed, here, the system is controlled by the river flow, therefore the discharge is the main driver, as also recognizable looking at the average flow velocity (Figure 5b). Despite the distance between the gauging station and the sea, however, tides can have an influence in the case of relatively low river discharges, as visible from the high water levels measured around April 27 (Figure 5a,b). On the other part, such an oscillation does not force significant changes in terms of bed shear stress and median diameter, given that these quantities are driven by the flow velocity (i.e., by the river discharge). A similar behavior is observable also at the E16 core location, as shown in Figure 6.
Being correlated to the flow velocity, the bed shear stress has a similar trend (Figure 5c,d), while the active mean diameter has a rather dissimilar behavior. In this case, the discharges observed in April can mobilize sediments having a diameter around 140–150 μm, while, during the flooding event at the beginning of May, the model simulates the transport of particles with smaller diameters. This could be correlated to the presence of finer sediments along the banks, which will be moved in the case of floods, occupying the floodplains.
To summarize, the model points out that, at the Pila gauging station, the hydro-morphodynamics are controlled by the river flow in the case of discharges higher than 1000 m3/s, in accordance with Maicu et al. [12].

3.2.2. Location of the Core Sample, Cross-Section E16

As reported in Figure 6a, the effects of the tidal oscillations are recognizable in the river discharge computed at the cross-section E16: generally, the lower the discharge, the higher the influence of the tide on the prodelta, even if during flooding conditions. However, given that the water level is mostly controlled by the tidal variations, oscillations of the river discharge have a lower impact on it. The proposed 1D model computes the depth-averaged velocity along the longitudinal direction (Figure 6b), therefore it is not possible to derive a velocity field from this simple model. Nonetheless, it is possible to recognize the tidal oscillations also from the computed velocities, even if with lower resolution.
Lower discharges are associated with a lower stream power and a smaller bed shear stress, resulting in a lower capacity for moving sediments (Figure 6c). Furthermore, the shear stress, and consequently the sediment transport and the erosion/deposition processes, are clearly controlled by the tidal oscillations, even during high flow conditions. For high discharges (flooding conditions), the influence of the downstream boundary conditions (tide) is less significant, with respect to low discharges, therefore these oscillations are smoothed.
The modest shear stress results in low transport capacity, and, therefore, only the finest fractions of sediments can be moved offshore (Figure 6d). In particular, the flooding event contributed to transport of sediments having a median diameter >30 μm, while during low flow conditions, only very fine material can be moved [31].

3.3. Variation of the Plume Dimensions

The extension of the plume was derived from a satellite image (250 m resolution) acquired on 2 May 2009 (Figure 1) by MODIS and available through the NASA Earth Observing System Data and Information System (EOSDIS). To test the influence of such dimensions on the final results, two additional simulations were performed assuming smaller or larger plume cross-section, but maintaining the depths as derived from the bathymetry [30]. In designing the new geometries, the section widths were changed of around 30% with respect to the values adopted in [31].
The computed values at the Pila mouth (Figure 7) show that the dimensions of the plume have practically no influence in the inland. There are only small changes in terms of median diameter (Figure 7b) and bed shear stress (Figure 7c) during the flooding event, while no variations are visible for lower discharges during April 2009.
On the other part, a higher influence of the plume dimensions is recognizable looking at the prodelta. Figure 8 reports four parameters computed at the core location E16, ca. 2 km downstream of the Pila channel mouth (see Figure 1). Comparing the calibrated simulation (black line) with the one performed assuming narrower sections (blue line), it is possible to see that smaller dimensions drive to higher velocities (Figure 8d) and consequently higher bed shear stresses (Figure 8c), which are more prone to move coarser sediments (Figure 8b), especially during flooding conditions. On the contrary, by imposing larger sections (red line), an opposite behavior can be observed.

4. Discussion

Being a major source of freshwater, the Po River plays a fundamental role in driving the physical and biogeochemical processes of the Adriatic Sea. For this reason, in the past, several researchers focused on representing the river plume variability, accounting for both the inland drivers (sediments, water discharge) and the sea-related variables (wind, currents, and wave resuspension).
Even if derived from a relatively simple 1D modeling approach, the results reported here are in line with what was computed by Bever et al. [16] using a more complex 3D model. Indeed, by combining the river hydrodynamics with sea currents, waves, and wind, they have shown that the sediment deposition occurred directly offshore of the distributary mouths, with the highest accumulation near the largest source, the Pila distributary, as a combination of the buoyant river plume, wave resuspension, and suspended transport due to the currents of the Adriatic Sea. With respect to previous works [17,35,36], in the present study, the focus was on understanding the general dynamics of a Po River distributary channel. Therefore, gravity flows were not considered here for representing the channel hydro-morphodynamics during the Spring 2009 flood event. Indeed, as pointed out in a previous report [16], except on the steepest portions of the delta front, suspended transport and plume settling are likely to dominate sediment flux during floods, and therefore representing high-density gravity flows is not strictly needed.
Falcieri et al. [37] applied a 3D hydrodynamic model for reproducing the Po River plume pattern variability at different temporal scales, showing that, on long time scales, river discharges represent the dominant forcing in defining the plume size and surface pattern, while on time scales of few days the plume dynamics are modulated mostly by the wind structure. In particular, they demonstrated that, during the period May 2008–July 2009, the river discharges were the dominant driver in shaping the Po River plume, mostly because the weak wind forcing, according to the results reported here.
During high flow conditions, the Po River is a river-dominated system, and the resulting high-water profile can be adequately simulated via the classic backwater profile [4,38], similarly to other fluvial systems worldwide. As an example, Leonardi et al. [39] monitored the flow velocity and the water levels of a distributary channel of the Apalachicola Delta (Florida, USA) by using an acoustic Doppler current profiler (ADCP), pointing out an M2 profile for relatively high tides. Based on such measures, the authors showed that, in this micro-tidal system, a drop in sea level during low tides forces the drawdown of the backwater profile of the river, thus increasing the river flow velocity. The HEC-RAS model applied in the present study shows a similar result, with steeper backwater profiles during low tides (Figure 4b, red line), corresponding to peaking flow velocity.
On the contrary, the code implemented in this application cannot be applied to deltas having very complex geometries or tidal-dominated estuaries like the Fly River [21,40,41], where neap tides generate a bed shear stress higher than the critical one, therefore moving large quantities of sediments that redeposit, preventing the erosion of islands and bars [22].
Focusing on the location of the core sample, as the computed sediment transport is correlated with the velocity, in particular the suspended load, it is possible to conclude that the tidal oscillations have a major influence on the sediment transport during low flow conditions, while, in the case of high flow conditions, it is the river discharge that drives the overall process. Despite small, the oscillations of the river profile due to the tides during the peak of the flood can still affect the sediment transport capacity of the river and the sediment in suspension, leading to the formation of a tidal-modulated flood deposit in the prodelta [31].

5. Conclusions

This study shows that 1D numerical models can be adopted for describing the hydro-morpho dynamics of a distributary channel during a river flood event, in particular in the case of a fluvial-dominated deltaic system, as the Italian Po River. Using a combination of data acquired at different times and with multiple techniques, this study highlights the reliability of a 1D HEC-RAS model in reproducing the channel hydro-morphodynamics and the sediment data derived from a core sample taken a few km offshore during the flooding event of April–May 2009.
The use of simple codes with a reduced mathematical complexity, combined with data having a rather coarse scale, allows for quick reproduction of observed phenomena, especially in the case of relatively simple geometry, and for describing the effects of small tidal changes on a fluvial-dominated environment, like the Po River mouth. The results obtained with this simplistic approach can be used for implementing more complete models, in particular in the case of environments characterized by complex 2D or even 3D dynamics. In essence, to adequately describe riverine morphodynamics in tidal environments, 1D models can be used to provide the general picture, while providing more complex codes to allow for a detailed description of the involved mechanisms.

Author Contributions

Conceptualization, M.N. and V.M.; numerical modeling, M.N. and A.V.; data acquisition, V.M.; formal analysis, M.N., A.V., and V.M.; writing, M.N., V.M., and A.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The work done by Michael Nones and Arianna Varrani was supported within statutory activities No. 3841/E-232 41/S/2020, Ministry of Science and Higher Education of Poland. We greatly acknowledge the two reviewers and the Editor for their comments, which were paramount in improving the quality of our manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vousdoukas, M.; Ranasinghe, R.; Mentaschi, L.; Plomaritis, T.A.; Athanasiou, P.; Luijendijk, A.; Feyen, L. Sandy coastlines under threat of erosion. Nat. Clim. Chang. 2020, 10, 260–263. [Google Scholar] [CrossRef]
  2. Syvitski, J.P.; Kettner, A.J.; Overeem, I.; Hutton, E.W.; Hannon, M.T.; Brakenridge, G.R.; Day, J.; Vörösmarty, C.; Saito, Y.; Giosan, L.; et al. Sinking deltas due to human activities. Nat. Geosci. 2009, 2, 681–686. [Google Scholar] [CrossRef]
  3. Maselli, V.; Trincardi, F. Man made deltas. Sci. Rep. 2013, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Nienhuis, J.H.; Ashton, A.D.; Edmonds, D.A.; Hoitink, A.J.F.; Kettner, A.J.; Rowland, J.C.; Törnqvist, T.E. Global-Scale human impact on delta morphology has led to net land area gain. Nature 2020, 577, 514–518. [Google Scholar] [CrossRef]
  5. Geleynse, N.; Storms, J.E.; Walstra, D.-J.R.; Jagers, H.A.; Wang, Z.-B.; Stive, M. Controls on river delta formation; insights from numerical modelling. Earth Planet. Sci. Lett. 2011, 302, 217–226. [Google Scholar] [CrossRef]
  6. Baar, A.; Boechat-Albernaz, M.; Van Dijk, W.M.; Kleinhans, M.G. Critical dependence of morphodynamic models of fluvial and tidal systems on empirical downslope sediment transport. Nat. Commun. 2019, 10, 1–12. [Google Scholar] [CrossRef] [Green Version]
  7. Luan, H.L.; Ding, P.; Wang, Z.-B.; Yang, S.L.; Lu, J.Y. Morphodynamic impacts of Large-Scale engineering projects in the Yangtze River delta. Coast. Eng. 2018, 141, 1–11. [Google Scholar] [CrossRef]
  8. Nones, M.; Guerrero, M.; Ronco, P. Opportunities from low-resolution modelling of river morphology in remote parts of the world. Earth Surf. Dyn. 2014, 2, 9–19. [Google Scholar] [CrossRef] [Green Version]
  9. Walton, R.; Jonas, M.; Stroppini, A.; Martin, B. Hydraulic model study of Delta-Mendota canal. Impacts Glob. Clim. Chang. 2005, 1–10. [Google Scholar] [CrossRef]
  10. Beltaos, S. The 2014 Ice–Jam flood of the Peace-Athabasca delta: Insights from numerical modelling. Cold Reg. Sci. Technol. 2018, 155, 367–380. [Google Scholar] [CrossRef]
  11. Rehman, M.A.; Naeem, U.A.; Hashmi, H.N.; Shakir, A.S. Possible options to slow down the advancement rate of Tarbela delta. Environ. Monit. Assess. 2017, 190, 39. [Google Scholar] [CrossRef] [PubMed]
  12. Hanegan, K.; Georgiou, I. Tidal Modulated Flow and Sediment Flux through Wax Lake Delta Distributary Channels: Implications for Delta Development. In Proceedings of the International Association of Hydrological Sciences, New Orleans, LA, USA, 11–14 December 2014; pp. 391–398. [Google Scholar]
  13. Maicu, F.; De Pascalis, F.; Ferrarin, C.; Umgiesser, G. Hydrodynamics of the Po River-Delta-Sea System. J. Geophys. Res. Oceans 2018, 123, 6349–6372. [Google Scholar] [CrossRef]
  14. Wester, S.J.; Grimson, R.; Minotti, P.; Booij, M.J.; Brugnach, M. Hydrodynamic modelling of a tidal delta wetland using an enhanced quasi-2D model. J. Hydrol. 2018, 559, 315–326. [Google Scholar] [CrossRef]
  15. Gugliotta, M.; Saito, Y. Matching trends in channel width, sinuosity, and depth along the fluvial to marine transition zone of Tide-Dominated river deltas: The need for a revision of depositional and hydraulic models. Earth-Sci. Rev. 2019, 191, 93–113. [Google Scholar] [CrossRef]
  16. Bever, A.; Harris, C.; Sherwood, C.R.; Signell, R. Deposition and flux of sediment from the Po River, Italy: An idealized and wintertime numerical modeling study. Mar. Geol. 2009, 260, 69–80. [Google Scholar] [CrossRef]
  17. Ambrosi, D.; Corti, S.; Pennati, V.; Saleri, F. Numerical simulation of unsteady flow at Po River delta. J. Hydraul. Eng. 1996, 122, 735–743. [Google Scholar] [CrossRef]
  18. Traykovski, P.; Wiberg, P.; Geyer, W. Observations and modeling of wave-supported sediment gravity flows on the Po prodelta and comparison to prior observations from the Eel shelf. Cont. Shelf Res. 2007, 27, 375–399. [Google Scholar] [CrossRef] [Green Version]
  19. Amorosi, A.; Maselli, V.; Trincardi, F. Onshore to offshore anatomy of a late quaternary Source-To-Sink system (Po Plain–Adriatic Sea, Italy). Earth-Sci. Rev. 2016, 153, 212–237. [Google Scholar] [CrossRef]
  20. Maselli, V.; Pellegrini, C.; Del Bianco, F.; Mercorella, A.; Nones, M.; Crose, L.; Guerrero, M.; Nittrouer, J.A. River morphodynamic evolution under Dam-Induced backwater: An example from the Po River (Italy). J. Sediment. Res. 2018, 88, 1190–1204. [Google Scholar] [CrossRef] [Green Version]
  21. Domeneghetti, A.; Carisi, F.; Castellarin, A.; Brath, A. Evolution of flood risk over large areas: Quantitative assessment for the Po River. J. Hydrol. 2015, 527, 809–823. [Google Scholar] [CrossRef]
  22. Nones, M. Numerical modelling as a support tool for river habitat studies: An Italian case study. Water 2019, 11, 482. [Google Scholar] [CrossRef] [Green Version]
  23. Syvitski, J.; Kettner, A.J. On the flux of water and sediment into the Northern Adriatic Sea. Cont. Shelf Res. 2007, 27, 296–308. [Google Scholar] [CrossRef]
  24. Cozzi, S.; Giani, M. River water and nutrient discharge in the Northern Adriatic Sea: Cur-Rent importance and long term changes. Cont. Shelf Res. 2011, 31, 1881–1893. [Google Scholar] [CrossRef]
  25. Montanari, A. Hydrology of the Po River: Looking for changing patterns in river discharge. Hydrol. Earth Syst. Sci. 2012, 16, 3739–3747. [Google Scholar] [CrossRef] [Green Version]
  26. Syvitski, J.; Kettner, A.J.; Correggiari, A.; Nelson, B.W.; Correggiari, A. Distributary channels and their impact on sediment dispersal. Mar. Geol. 2005, 222, 75–94. [Google Scholar] [CrossRef]
  27. Wheatcroft, R.A.; Stevens, A.W.; Hunt, L.M.; Milligan, T.G. The Large-Scale distribution and internal geometry of the fall 2000 Po River flood deposit: Evidence from digital X-radiography. Cont. Shelf Res. 2006, 26, 499–516. [Google Scholar] [CrossRef]
  28. Brunner, G.W. HEC-RAS, River Analysis System: Hydraulic Reference Manual, Version 5.0. 2016. Available online: hec.usace.army.mil/software/hec-ras/documentation/HEC-RAS%205.0%20Reference%20Manual.pdf (accessed on 10 January 2020).
  29. Tesi, T.; Miserocchi, S.; Goni, M.; Turchetto, M.; Langone, L.; De Lazzari, A.; Albertazzi, S.; Correggiari, A. Influence of distributary channels on sediment and organic matter supply in event-dominated coastal margins: The Po prodelta as a study case. Biogeosciences 2011, 8, 365–385. [Google Scholar] [CrossRef] [Green Version]
  30. Trincardi, F.; Campiani, E.; Correggiari, A.; Foglini, F.; Maselli, V.; Remia, A. Bathymetry of the Adriatic Sea: The legacy of the last eustatic cycle and the impact of modern sediment dispersal. J. Maps 2013, 10, 151–158. [Google Scholar] [CrossRef] [Green Version]
  31. Maselli, V.; Normandeau, A.; Nones, M.; Tesi, T.; Langone, L.; Trincardi, F.; Bohacs, K.M. Tidal modulation of River-Flood deposits: How low can you go? Geology 2020. [Google Scholar] [CrossRef]
  32. Nones, M.; Archetti, R.; Guerrero, M. Time-Lapse Photography of the Edge-of-Water Line Displacements of a Sandbar as a Proxy of Riverine Morphodynamics. Water 2018, 10, 617. [Google Scholar] [CrossRef] [Green Version]
  33. Nones, M.; Pugliese, A.; Domeneghetti, A.; Guerrero, M. Po River Morphodynamics Modelled with the Open-source Code iRIC. In Recent Trends in Environmental Hydraulics; Kalinowska, M.B., Mrokowska, M.M., Rowiński, P.M., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 335–346. [Google Scholar]
  34. SIMPO. Study and Preliminary Design of the Hydraulic Regulation of the River Po Main Channel for Civil and Environmental Purposes; Magistrato per il Po—SIMPO Spa; Technical Report; SIMPO: Paulo, Brazil, 1980. (In Italian) [Google Scholar]
  35. Zanichelli, G.; Caroni, E.; Fiorotto, V. River bifurcation analysis by physical and numerical modeling. J. Hydraul. Eng. 2004, 130, 237–242. [Google Scholar] [CrossRef]
  36. Friedrichs, C.T.; Scully, M.E. Modeling deposition by Wave-Supported gravity flows on the Po River prodelta: From seasonal floods to prograding clinoforms. Cont. Shelf Res. 2007, 27, 322–337. [Google Scholar] [CrossRef]
  37. Falcieri, F.M.; Benetazzo, A.; Sclavo, M.; Russo, A.; Carniel, S. Po River plume pattern variability investigated from model data. Cont. Shelf Res. 2014, 87, 84–95. [Google Scholar] [CrossRef]
  38. Lamb, M.P.; Nittrouer, J.A.; Mohrig, D.; Shaw, J. Backwater and river plume controls on scour upstream of river mouths: Implications for Fluvio-Deltaic morphodynamics. J. Geophys. Res. Space Phys. 2012, 117. [Google Scholar] [CrossRef]
  39. Leonardi, N.; Kolker, A.S.; Fagherazzi, S. Interplay between river discharge and tides in a delta distributary. Adv. Water Resour. 2015, 80, 69–78. [Google Scholar] [CrossRef]
  40. Dalrymple, R.W.; Baker, E.K.; Harris, P.T.; Hughes, M.G.; Sidi, F.H.; Nummedal, D.; Posamentier, H.W. Sedimentology and stratigraphy of a tide-dominated, foreland-basin delta (Fly River, Papua New Guinea). In Tropical Deltas of Southeast Asia–Sedimentology, Stratigraphy, and Petroleum Geology; Sidi, F.H., Nummedal, D., Imbert, P., Darman, H., Posamentier, H.W., Eds.; SEPM Special Publication: Tulsa, OK, USA, 2003; Volume 76, pp. 147–173. [Google Scholar]
  41. Canestrelli, A.; Defina, A.; Lanzoni, S.; Fagherazzi, S. Erosional power and morphodynamic tendencies in the Fly River Delta, Papua New Guinea. In Proceedings of the River, Coastal and Estuarine Morphodynamics RCEM Conference 2009, Santa Fe, Argentina, 21–25 September 2009; Volume 2, p. 355. [Google Scholar]
Figure 1. (a) Location of the Po River catchment in Italy; (b) elevation map; (c) location of the gauging stations downstream of Pontelagoscuro, in the lowland part of the Po River; (d) locations of the Pila gauging station and the core sample E16.
Figure 1. (a) Location of the Po River catchment in Italy; (b) elevation map; (c) location of the gauging stations downstream of Pontelagoscuro, in the lowland part of the Po River; (d) locations of the Pila gauging station and the core sample E16.
Geosciences 10 00209 g001
Figure 2. (a) Measured water level and river discharge at the gauging stations of (a) Pila), (b) Cavanella, and (c) Polesella; (d) predicted astronomical tidal amplitude and (e) wave heights off the Pila distributary mouth.
Figure 2. (a) Measured water level and river discharge at the gauging stations of (a) Pila), (b) Cavanella, and (c) Polesella; (d) predicted astronomical tidal amplitude and (e) wave heights off the Pila distributary mouth.
Geosciences 10 00209 g002
Figure 3. Calibration and validation of the model at the Polesella gauging station: (a) comparison between measured (blue line) and computed (red line) water levels; (b) computed (y-axis) and measured (x-axis) water levels during the study period.
Figure 3. Calibration and validation of the model at the Polesella gauging station: (a) comparison between measured (blue line) and computed (red line) water levels; (b) computed (y-axis) and measured (x-axis) water levels during the study period.
Geosciences 10 00209 g003
Figure 4. Simulated water levels along the Po River, with respect to the position of the core samples (0 km): (a) water profiles for maximum discharges (Q = 8000 m3/s) at high tide (red line) and low tide (blue line), and for low flow conditions (Q = 1000 m3/s, green line). The thalweg is indicated with the black line, while the ground elevation is reported on the secondary y-axis; (b) zoom of the last 10 km, focusing on the water levels.
Figure 4. Simulated water levels along the Po River, with respect to the position of the core samples (0 km): (a) water profiles for maximum discharges (Q = 8000 m3/s) at high tide (red line) and low tide (blue line), and for low flow conditions (Q = 1000 m3/s, green line). The thalweg is indicated with the black line, while the ground elevation is reported on the secondary y-axis; (b) zoom of the last 10 km, focusing on the water levels.
Geosciences 10 00209 g004
Figure 5. Simulated parameters at the gauging station of Pila: (a) river discharge (black line) vs. water level (grey line); (b) averaged flow velocity (red line) vs. water level (grey line); (c) bed shear stress (green line) vs. active median diameter d50 (blue line); (d) bed shear stress (green line) vs. averaged flow velocity (red line).
Figure 5. Simulated parameters at the gauging station of Pila: (a) river discharge (black line) vs. water level (grey line); (b) averaged flow velocity (red line) vs. water level (grey line); (c) bed shear stress (green line) vs. active median diameter d50 (blue line); (d) bed shear stress (green line) vs. averaged flow velocity (red line).
Geosciences 10 00209 g005
Figure 6. Simulated parameters at the E16 core location: (a) river discharge (black line) vs. water level (grey line); (b) averaged flow velocity (red line) vs. water level (grey line); (c) bed shear stress (green line) vs. active median diameter d50 (blue line); (d) bed shear stress (green line) vs. averaged flow velocity (red line).
Figure 6. Simulated parameters at the E16 core location: (a) river discharge (black line) vs. water level (grey line); (b) averaged flow velocity (red line) vs. water level (grey line); (c) bed shear stress (green line) vs. active median diameter d50 (blue line); (d) bed shear stress (green line) vs. averaged flow velocity (red line).
Geosciences 10 00209 g006
Figure 7. Comparison between calibrated (black line), narrow (blue line) and large (red line) sections at the gauging station of Pila: (a) effective width W; (b) median diameter d50; (c) bed shear stress τ; (d) averaged flow velocity v.
Figure 7. Comparison between calibrated (black line), narrow (blue line) and large (red line) sections at the gauging station of Pila: (a) effective width W; (b) median diameter d50; (c) bed shear stress τ; (d) averaged flow velocity v.
Geosciences 10 00209 g007
Figure 8. Comparison between calibrated (black line), narrow (blue line) and large (red line) sections at the E16 core location: (a) effective width W; (b) median diameter d50; (c) bed shear stress τ; (d) averaged flow velocity v.
Figure 8. Comparison between calibrated (black line), narrow (blue line) and large (red line) sections at the E16 core location: (a) effective width W; (b) median diameter d50; (c) bed shear stress τ; (d) averaged flow velocity v.
Geosciences 10 00209 g008

Share and Cite

MDPI and ACS Style

Nones, M.; Maselli, V.; Varrani, A. Numerical Modeling of the Hydro-Morphodynamics of a Distributary Channel of the Po River Delta (Italy) during the Spring 2009 Flood Event. Geosciences 2020, 10, 209. https://doi.org/10.3390/geosciences10060209

AMA Style

Nones M, Maselli V, Varrani A. Numerical Modeling of the Hydro-Morphodynamics of a Distributary Channel of the Po River Delta (Italy) during the Spring 2009 Flood Event. Geosciences. 2020; 10(6):209. https://doi.org/10.3390/geosciences10060209

Chicago/Turabian Style

Nones, Michael, Vittorio Maselli, and Arianna Varrani. 2020. "Numerical Modeling of the Hydro-Morphodynamics of a Distributary Channel of the Po River Delta (Italy) during the Spring 2009 Flood Event" Geosciences 10, no. 6: 209. https://doi.org/10.3390/geosciences10060209

APA Style

Nones, M., Maselli, V., & Varrani, A. (2020). Numerical Modeling of the Hydro-Morphodynamics of a Distributary Channel of the Po River Delta (Italy) during the Spring 2009 Flood Event. Geosciences, 10(6), 209. https://doi.org/10.3390/geosciences10060209

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop