Next Article in Journal
A Modified kε Turbulence Model for a Wave Breaking Simulation
Next Article in Special Issue
Documenting a Century of Coastline Change along Central California and Associated Challenges: From the Qualitative to the Quantitative
Previous Article in Journal
Correction: Tapoglou, E., et al. Climate Change Impact on the Frequency of Hydrometeorological Extremes in the Island of Crete. Water 2019, 11, 587
Previous Article in Special Issue
Impacts of Sea Level Rise and River Discharge on the Hydrodynamics Characteristics of Jakarta Bay (Indonesia)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Reconstruction of Eroded Beach Ridges at the Ombrone River Delta

1
GNRAC—Gruppo Nazionale per la Ricerca in Ambiente Costiero, Corso Europa, 26-16132 Genova, Italy
2
GeoCoste s.n.c., Via Corsi, 19-50141 Florence, Italy
3
Department of Earth Science, University of Florence, Via Micheli, 4-50121 Florence, Italy
*
Author to whom correspondence should be addressed.
Water 2019, 11(11), 2281; https://doi.org/10.3390/w11112281
Submission received: 10 September 2019 / Revised: 25 October 2019 / Accepted: 28 October 2019 / Published: 31 October 2019
(This article belongs to the Special Issue Coastal Dynamic and Evolution)

Abstract

:
Several remotely sensed images, acquired by different sensors on satellite, airplane, and drone, were used to trace the beach ridges pattern present on the delta of the River Ombrone. A more detailed map of these morphologies, than those present in the literature, was obtained, especially at the delta apex, where beach ridges elevation in minor. Beach ridges crests, highlighted through image enhancement using ENVI 4.5 and a DTM based on LiDAR data, were then processed with ArcGIS 9.3 software. Starting from this map, a method to reconstruct beach ridges segments deleted by the transformations of the territory is proposed in this paper. The best crest-lines fitting functions were calculated through interpolation of their points with Curve Expert software, and further extrapolated to reconstruct the ridges morphology where human activity, riverbed migration, or coastal erosion eliminated them. This allowed to reconstruct the ridges pattern also offshore the present delta apex, where the shoreline retreated approximately 900 m in the last 150 years. Results can be further used to implement conceptual and numerical models of delta evolution.

1. Introduction

Deltas are dynamic depositional landforms, sensitive to changes in both the terrestrial and marine environment [1,2], constantly reshaped by river input, tides, and waves [3,4]. Wave-dominated deltas are characterized by a cuspate morphology, a steep offshore profile, and the presence of more or less curved beach ridges parallel to past shoreline positions [5]. This allowed to reconstruct Holocene delta formation phases for many rivers, and to date, the various morphologies when archaeological remnants were present or if radiometric dating were executed [6,7].
Fluctuations in fluvial sediment supply were also assessed for several deltas, being the most important factor controlling the spatial and temporal variation in beach ridges development [8,9]. Spits are frequently present near the river mouth, sometimes overpassed by the coastal progradation and detectable in the delta area, others eroded by waves. Beach ridges frequently evolve in coastal dunes which reflect a huge amount of sediments that are redistributed by waves and moved inland by winds forming elongated and parallel morphologies behind the coastline [10]. A fast accreting coast makes the formation of several low crests, whereas a slowly prograding or a stable one allows the growth of few high dunes [11,12]. Beach ridges morphology and related pattern could be used as markers of morphodynamic variations: they can give information on past wave regime, climate conditions, sediment supply, and sea level change [13,14,15,16].
In historical times, fast delta progradation is frequently induced by fluvial input increase, ascribed to the population development on the catchment area accompanied by intensive deforestation [17,18].
Just short (months) and medium time (years) changes in sediment input can induce delta apex reshaping, resulting in different delta lobes morphology and beach ridges patterns. Fast progradation periods, due to higher sediment input by the river, sees a sharp cusp and the rotation of the river mouth towards the dominant wave front direction, as described for the Ombrone and Arno River deltas [19]; on the contrary, a reduced fluvial input determines a flattening of the cusp up to a rectification of the coastline and downdrift shift of the river mouth [20,21].
Conceptual models of cuspate delta evolution have been proposed by ref. [19,22], whereas other authors ref. [2,21,23,24] used numerical models to explain delta asymmetry and wave reshaping of the river mouth during erosion and accretion phases. For deltas characterized by alternating phases of growth and erosion, not all the beach ridges at the delta apex can be identified, since their most offshore part has been beheaded during the erosion phases [12]; in addition, some segments can be leveled by human activity, especially agriculture, or cut by the river course migration.
A full reconstruction of the delta evolution needs to know the delta shape at its extreme expansions, even if later cut by erosion; this is of outmost importance when conceptual or numerical models based on wings surface are to be developed.
The possibility to observe the secular trends of beach ridges on the Ombrone delta plain represents a valuable opportunity to obtain useful information on its complex geomorphological evolution, where one of the main drawbacks is represented by the loss of geographical information due to natural and anthropic transformations of the territory over time and the lack, or inaccuracy, of ancient cartography. The problem was solved through a mathematical reconstruction of the missing portions of the beach ridges inspired form by the observation of their geometry and how their planar forms remind of specific mathematical functions.
On Earth, a number of geomorphic elements have a regular shape, including small-scale landforms such as cirques, drumlins, sand dunes, alluvial fans, dolines, polygonal soils and cracks, columnar jointing, and scoria cones, to large-scale landforms such as composite and shield volcanoes [25]. Geomorphometry is the science of quantitative land-surface analysis and a variety of landforms are fitted by mathematical functions, often by means of regression analysis [26].
Just to remain in coastal environment, zeta bays in equilibrium with approaching refracted waves are fitted by a logarithmic spiral as proposed by Silvester and Hsu [27], or the equilibrium nearshore profile is mathematically described by the Dean equation [28].
Following a similar approach, a method to reconstruct beach ridges, where they have been eroded or razed by human activity, is hereafter proposed using fitting equations.
Results can be further used for a detailed model on the delta evolution reconstruction, which is not the aim of this study.

2. Study Area

The Ombrone River delta is in Italy, located on the Southern Tuscany coast (Italy), between Castiglione della Pescaia and Monti dell’Uccellina, in the Maremma Regional Park (Figure 1). The dominant waves direction is from Souh-West [24], and for the present research, is assumed to be a constant wave climate over the historical times [29,30]. The regional longshore transport direction is from south to north [31,32]; subsidence in the alluvial plain was measured in 3 mm/year [33] and this is the only significant vertical movement in the considered time lapse. Sea level rise is about 1 m since the roman period [34].
The ancient gulf of Grosseto started to fill with the material transported by the Ombrone River during the Holocene transgression, but the most internal beach ridge is very likely associated to a sand bar closing a lagoon. When the lagoon was almost silted, the river could directly empty into the sea and started the delta formation. This is dated to the early Roman period [35] and explained with the intense deforestation carried out within the watershed by the Etruscan, first, and by the Romans, later.
During the last twenty centuries, the delta apex grew for approximately 6.5 km at the apex, with a very fast expansion in the 14th–18th century [36], when an intense deforestation occurred. Some authors also see in this delta expansion, the effect of the Little Ice Age [35,37]. Two main erosional phases have been identified, the first after the fall of the Roman Empire, and the second with the Black Death epidemic; in both cases a strong population reduction induced the abandonment of many cultivated areas and the growth of the forest [36].
Furthermore, at the end of the 19th century, these tendencies drastically reversed with the erosion of the coastline at the river mouth, caused by agriculture abandonment, riverbed quarrying, dams and weirs construction; whereas wetland reclamations, through river course diversion, started some centuries before depriving the river load mostly in its finer fraction. Results of these works were poor since sediment compaction was restoring the original wetlands [38].
Since the end of the 19th century, the delta apex started to retreat, and it is now approximately 900 m back from its most prominent position (depicted in a 1881 topographic map). Before shore protection works were carried out in 2013–2014 on the southern lobe (submerged groins) and milder in 2016 on the northern one (short wood permeable groins), the erosion rate was approximately 10 m/year [17].
Delta morphology is slightly different in the two lobes [6,36,37]. The northern side hosts inter-dune swales and, near the river course, some permanent ponds, possibly old channels to connect salterns to the sea. Dunes progressively merge to the north in few parallel ridges; the inner one is attributed to the Etruscan period, the next five to the Roman time, followed by seven Medieval and the last four belonging to the Modern Age [39].
Agriculture activity and urban development at Marina di Grosseto deleted some segments of these dunes. The southern side is almost uninhabited, dunes are closer, and their continuity is more preserved, but at the apex, where they are far lower, seasonal ponds form and in the wetter period most of the apex is flooded. The Monti dell’Uccellina, with their seaward projection, limit sediment dispersion and, even if in the past centuries the prevalent longshore transport was to the north, the delta could symmetrically grow. (Smaller volume but shorter base in the southern side). The River Ombrone delta, hosts many beheaded ridges; but the well-preserved morphology of those remained allows to test a methodology to mathematically reconstruct the lost segments.

3. Material and Methods

For this study, a detailed beach ridge map was obtained through the combination of several data sets, produced by active and passive remote sensors, and processed with various methodologies. A LiDAR (light detection and ranging) survey acquired in 2008 by the Ministry of the Environment, with 1 × 1 m ground resolution, was used to derive a 3D model of the beach ridges area. Where agriculture flattened the ridges, ad hoc topographic surveys were carried out, and a higher vertical resolution DTM (digital terrain model) was produced from low altitude drone imagery using the structure from motion (SfM) technique.
Satellite images used span from Landsat TM and Landsat ETM+ (from 1986 to 2002); Quickbird (2004); Landsat 8 (2015) and Sentinel 2a (2017); in addition, some Google Earth images were used. Pixel size goes from 2.44 m (Quickbird) to 30 m (Landsat TM).
Remotely sensed images helped in tracing ridge crests where ridges were flattened, and soil moistures helps to discriminate sand from silty/clayey soils present in the interdune swales. This was obtained with specific image processing, such as normalized difference of vegetation index (NDVI), tasseled cup transformation (brightness, greenness, and wetness) [40], principal components analysis (PCA), and RGB (red-green-blue) color composites of the results of the different elaborations. Recent topographic surveys, regional cartography, and old aerial photos acquired in 1944 by RAF (Royal Air Force ) were also used. Images were then wrapped on the Lidar 3D model to assist their photo interpretation.
All these data were processed using ENVI 4.5 (Harris Geospatial Solutions Inc., Broomfield, CO, USA) and ArcGIS 9.3 (ESRI) softwares in order to produce a detailed beach dunes/beach ridges map of the delta area in the WGS84 UTM32 reference system (Figure 2) and leading information was derived from LiDAR data, whose resolution was maintained also after the following image overlapping.
From LiDAR-DTM, several cross-shore profiles were extracted, and minor beach ridges sequence and evolution phases were recognized and relatively dated, starting from the main delta evolution periods described in the literature [36,37].
From the position, geometry, and height of some ridges it was possible to pair them across the gap formed by the river or by agricultural flattering. In Figure 3, profiles crossing the dunes/ridges system on the two delta wings are drawn and ridges associated and numbered (see Figure 1 for profiles location). In the northern profile, all the morphologies are compressed, both because the profile is located far from the delta apex, and because sediment is distributed along a wider area (the northern limit of the cell—Punta delle Rocchette—is 24 km from the river mouth, whereas the southern, one—Monti dell’Uccellina—is 7 km only).
Profiles show the differences in height and distance between ridges, and how some of the highest/largest ones on the delta wings and associated to apex erosion or stability periods, comprise more than one crest, due to an overlapping of more dunes.
Sixteen beach ridges from the early 15th to the middle-19th century have been selected representing different periods of the delta evolution, both when accreting and when eroding. In the first case, beach ridges plant geometry is more curved at the apex with ridges sharply converging (Figure 4); on the contrary, soon after a period of delta erosion their planform is more rectilinear since beach retreats at the apex but progradation goes on at the wings (actually it is a form of beach rotation) [12].
Along each beach ridge, several x,y points were digitized where each ridge crest was more evident; one point every 50 m on average. Along this distance, ridges are almost rectilinear and more points do not increase the accuracy. They were further rotated and shifted in a Cartesian coordinate system with the y-axis coincident with the direction of the terminal river course.
In order to mathematically reconstruct the missing parts of the beach ridges, 33 families of linear and not linear regression models were tested, both for progradation and for erosion beach ridge shapes. The calculation was done using a specific curve fitting software (CurveExpert Pro 1.4, Hyams Development, Chattanooga, TN, USA). The software provides a ranking of the best fitting functions and they are also displayed in the graphing window with statistical results, such as standard error and correlation coefficient. Values x,y for points of interpolated and extrapolated curves are also available and useful for the reconstruction of the missing parts of the delated beach ridges.

4. Results

Finding equations that well fit the ridges allows to reconstruct their original plan form in the stretches where they have been erased by human activities or by coastal erosion.
Among all the 33 regression models tested, two distinct families of functions, respectively Weibull (Equation (1)) and Morgan–Morgan–Finney (MMF) (Equation (2)), were found out with lowest values of standard error and highest correlation coefficient. The former function better approximates more rectilinear beach ridges of older coastlines or of erosional phases (more flat lines); whereas the latter function better fits curved beach ridges of fast accreting stages (Table 1).
To a lesser extent, Hoerl (Equation (3)) and 3th degree polynomial (Equation (4)) were the following two better equations in the fitting ranking.
Weibull :   y = a be cx d .
MMF :   y =   ab + cx d b + x d .
Hoerl :   y = a   b x x c .
Polynomial :   3 th - degree   y =   a + bx + cx 2 + dx 3 .
Obviously, the robustness of the method depends on the extension of interpolated beach ridges too, and is correlated to the segment form: linear or arcuate. For very short segments, results exponentially lose significance.
In the table are reported precise year, just to have a temporal reference to be attributed to every mapped beach ridge, but there may be indeed an uncertainty of several years. The correlation coefficients range are from 0.999 to 0.997.
In Figure 4 and Figure 5, the real 1625 (apex erosion phase) and 1750 (apex accretion phase) beach ridges on the south lobe of the delta are drawn together with the best 4 fitting models produced through extrapolation, near the river, of the first 4 km of each ridge. The graphs show Weibull (Equation (1)) better approximates more rectilinear erosional beach ridges and MMF (Equation (2)) curved beach ridges of progradation phases.
Starting from all the beach ridges mapped in Figure 2, the main evolution phases of the Ombrone River delta, for the considered period, are reported in Figure 6. The considered period is from the early 15th to the middle 19th century, after which reliable topographic surveys are available.
Beach ridges deleted near the river because of watercourse shifting or truncated by a subsequent erosive event have been mathematically reconstructed according to the methodology proposed.
Beach ridges beheaded at the river mouth are represented in cyan and those related to the cuspate progradation in red. In blue and purple, their reconstructed morphologies with mathematical functions.

5. Discussion

Quantitative approach in Geomorphology became established in the mid-20th century [41], but equations fitting coastal landforms arrived a few decades later, like those aforementioned, by Dean [28] and by Silvester and Hsu [27]; both, although developed on actual morphologies, are used to forecast the shape a landform (nearshore profile and zeta-bays) will reach at equilibrium. What is done in this paper is the reconstruction of ancient landforms in those parts which have been later eroded.
In a period in which approximately 31% of the world beaches are eroding [42], and considering that erosion generally starts at the river mouth and gradually expand laterally [22], many present cuspate river deltas lost their apex. Further, agricultural works and urbanization frequently flattened, or at least hid, original morphologies. This makes it impossible to fully reconstruct the natural delta morphology and study its historical evolution.
The proposed method needs high-resolution data, which was obtained from active and passive remote sensing systems. In this way, a detailed reconstruction of beach ridges morphologies was done, and several accretion and erosion phases mapped. Thirty-three different fitting curves were generated though interpolation of control points on the beach ridges/dune crests and their accuracy in describing those morphologies was assessed.
The further extrapolation of the last seaward points allowed to reconstruct the ridge segments which were eroded during the last 150 years, thus obtaining information on the past delta apex shape. On the other sides, interpolation made it possible to identify reaches flattened by agriculture or by other anthropogenic activities. For ridges formed during delta expansions phases (more curved at the tip), the MMF model proved to be more accurate than all the other tested; whereas ridges formed during delta tip erosion, more flattens, were better reconstructed by the Weibull equation.

6. Conclusions

The mathematical interpolation and extrapolation of the functions fitting the ridge crests at the Ombrone River delta allowed to reconstruct their entire morphologies, not only in the part eroded during the last 150 years, but also where similar processes worked in the Middle and Modern Age.
On these bases, it will be possible to develop more accurate models of delta evolution in order to better explain the different morphology characterizing the updrift and the downdrift delta sides, and to forecast future shoreline displacement in this coast rich in natural and economical elements.
The hoped application of this method to other deltas, using different equation coefficients, could help to solve similar problems and improve the method itself.

Author Contributions

Conceptualization, I.M., L.R., E.P.; methodology, I.M., L.R., E.P.; validation, L.R.; formal analysis, I.M.; investigation, I.M.; resources, E.P.; data curation, L.R.; writing—original draft preparation, I.M. and L.R.; writing—review and editing, L.R. and E.P.; supervision, E.P.

Funding

This work was supported by the Project Interreg Marittimo-IT-FR Maritime “Management des Risques de l’Erosion cotière et actions de Gouvernance Transfrontalière” (MAREGOT).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Syvitski, J.P.; Milliman, J.D. Geology, geography and humans battle for dominance over the delivery of sediment to the coastal ocean. J. Geol. 2007, 115, 1–19. [Google Scholar] [CrossRef]
  2. Nienhuis, J.H.; Ashton, A.D.; Roos, P.C.; Hulscher, S.J.M.H.; Giosan, L. Wave reworking of abandoned deltas. Geophys. Res. Lett. 2013, 40, 5899–5903. [Google Scholar] [CrossRef] [Green Version]
  3. Wright, L.D.; Coleman, J.M. Variations in morphology of major river deltas as functions of ocean wave and river discharge regimes. AAPG Bull. 1973, 57, 370–398. [Google Scholar]
  4. Galloway, W.E. Process framework for describing the morphologic and stratigraphic evolution of deltaic depositional systems. In Deltas, Models for Exploration; Broussard, M.L., Ed.; Houston Geological Society: Houston, TX, USA, 1975; pp. 87–98. [Google Scholar]
  5. Van Maren, D.S.; Hoekstra, P. Cyclic development of a wave dominated delta. In Proceedings of the Coastal Sediments ’03, Clearwater Beach, FL, USA, 18–23 May 2013; p. 12. [Google Scholar]
  6. Bellotti, P.; Caputo, C.; Davoli, L.; Evangelista, S.; Valeri, P. Lineamenti morfologici e sedimentologici della piana deltizia del Fiume Ombrone (Toscana Meridionale). Boll. Soc. Geol. It. 1999, 118, 141–148. [Google Scholar]
  7. Otvos, E.G. Beach ridges–definition and significance. Geomorphology 2000, 32, 83–103. [Google Scholar] [CrossRef]
  8. Anthony, E.J. Beach ridge progradation in response to sediment supply examples from West Africa. Marine Geol. 1995, 129, 175–186. [Google Scholar] [CrossRef]
  9. Anthony, E.J. Wave influence in the construction, shaping and destruction of river deltas: A review. Marine Geol. 2015, 361, 53–78. [Google Scholar] [CrossRef]
  10. Vespremeanu-Stroe, A.; Preoteasa, L.; Zainescu, F.; Rotaru, S.; Croitoru, L.; Timar-Gabor, A. Formation of Danube delta beach ridge plains and signatures in morphology. Quat. Int. 2016, 415, 268–285. [Google Scholar] [CrossRef]
  11. Psuty, N.P. Spatial variation in coastal foredune development. In Coastal Dunes: Geomorphology, Ecology and Management; Balkema: Rotterdam, The Netherlands, 1992; pp. 3–13. [Google Scholar]
  12. Pranzini, E. Airborne LIDAR survey applied to the analysis of the historical evolution of the Arno River delta (Italy). J. Coast. Res. 2007, 400–409. [Google Scholar]
  13. Curray, J.R.; Emmel, F.J.; Crampton, P.J.S. Holocene history of a strand plain, Lagoonal Coast, Nayarit, Mexico. In Lagunas Costeras; Costonares, A.A., Phelger, F.B., Eds.; Un Simposio; Universidad Nacional Autonoma de Mexico: Mexico City, Mexico, 1967; p. 686. [Google Scholar]
  14. Stapor, F.W. Beach ridge and beach ridge coast. In Encyclopedia of Beaches and Coastal Environments; Schwartz, M.L., Ed.; Hutchinson Ross: Stroudsburg, PA, USA, 1982; pp. 160–161. [Google Scholar]
  15. Fairbridge, R.W.; Hillaire-Marcel, C. An 8000-yr palaeoclimatic record of the double-hale 45-yr solar cycle. Nature 1977, 268, 413–416. [Google Scholar] [CrossRef]
  16. Taylor, M.; Stone, G.W. Beach-bridges: A review. J. Coast. Res. 1996, 12, 612–621. [Google Scholar]
  17. Pranzini, E. Bilancio sedimentario ed evoluzione storica delle spiagge. Il Quaternario 1994, 7, 197–202. [Google Scholar]
  18. Syvitski, J.P.M.; Vörösmarty, C.J.; Kettner, A.J.; Green, P. Impact of humans on the flux of terrestrial sediment to the global coastal ocean. Science 2005, 308, 376–380. [Google Scholar] [CrossRef] [PubMed]
  19. Pranzini, E. Updrift river mouth migration on cuspate deltas: Two examples from the coast of Tuscany (Italy). Geomorphology 2001, 38, 125–132. [Google Scholar] [CrossRef]
  20. Komar, P.D. Computer model of delta growth due to sediment input rivers and longshore transport. Geol. Soc. Am. Bull. 1973, 84, 2217–2219. [Google Scholar] [CrossRef]
  21. Nienhuis, J.H.; Ashton, A.D.; Giosan, L. Littoral steering of deltaic channels. Earth Planet. Sci. Lett. 2016, 453, 204–214. [Google Scholar] [CrossRef] [Green Version]
  22. Pranzini, E. A model for cuspate delta erosion. In Proceedings of the 6th Sympoium on Coastal and Ocean Management/ASCE, Charleston, SC, USA, 11–14 July 1989; pp. 4345–4357. [Google Scholar]
  23. Ashton, A.D.; Giosan, L. Wave-angle control of delta evolution. Hydrol. Land Surf. Stud. 2011, 38, 4051–4056. [Google Scholar] [CrossRef]
  24. Aminti, P.; Pranzini, E. Variations in longshore sediment transport rates as a consequence of beach erosion in a cuspate delta. Proc. EUROCOAST Assoc. 1990, 130–134. [Google Scholar]
  25. Pike, R.J.; Evans, I.S.; Hengl, T. Geomorphometry: A brief guide. In Developments in Soil Science; Hengl, T., Reuter, H.I., Eds.; Geomorphometry—Concepts, Software, Applications; Elsevier: Amsterdam, The Netherlands, 2009; Volume 33, pp. 3–30. [Google Scholar]
  26. King, C.A.M. Geometrical forms in geomorphology. Int. J. Math. Educ. Sci. Technol. 2006, 2, 153–169. [Google Scholar] [CrossRef]
  27. Silvester, R.; Hsu, J.R.C. Coastal Stabilization. Innovative Concepts; PRT Prentice Hall: Upper Saddle River, NJ, USA, 1993; p. 578. [Google Scholar]
  28. Dean, R.G. Equilibrium beach profiles: U.S. Atlantic and Gulf coasts. In Ocean Engineering Report, 12; Department of Civil Engineering, University of Delaware: Newark, NJ, USA, 1977. [Google Scholar]
  29. Li, L.; Walstra, D.R.; Storms, J. The impact of wave-induced longshore transport on a delta-shoreface system. J. Sediment. Res. 2015, 85, 6–20. [Google Scholar]
  30. Meini, L.; Mucci, G.; Vittorini, S. Ricerche meteo marine sul litorale toscano: Centoventi anni di osservazioni metereologiche a Livorno (1857–1976). Boll. Soc. Geogl. It. 1979, VIII, 449–474. [Google Scholar]
  31. Aiello, E.; Bartolini, C.; Caputo, C.; D’Alessandro, L.; Fanucci, F.; Fierro, G.; Gnaccolini, M.; La Monica, G.B.; Lupia Palmieri, E.; Piccazzo, M.; et al. Il trasporto litoraneo lungo la costa toscana tra la foce del Fiume Magra ed i Monti dell’Uccellina. Boll. Soc. Geol. It. 1976, 94, 1519–1571. [Google Scholar]
  32. Aminti, P.L. Ricostruzione del clima ondoso della Toscana Meridionale sulla base di misure anemometriche. Bolettino Degli Ingegneri 1983, 4, 34–56. [Google Scholar]
  33. Salvioni, G. I movimenti del suolo nell’Italia centro-settentrionale. Bollettino di Geodesia e Scienze Affini 1957, 16, 325–366. [Google Scholar]
  34. Antonioli, F.; Silenzi, S. Variazioni relative del livello del mare e vulnerabilità delle pianure costiere italiane. Quaderni della Società Geologica Italiana 2007, 2, 1–29. [Google Scholar]
  35. Citter, C.; Arnoldus-Huyzendveld, A. Archeologia urbana a Grosseto, la città nel contesto geografico della bassa valle dell’Ombrone. Origine e sviluppo di una città medievale nella “Toscana delle città deboli”. In Le Ricerche 1997–2005; Edizioni all’Insegna del Giglio s.a.s: Fiorentino, Italy, 2007; pp. 41–60. [Google Scholar]
  36. Innocenti, L.; Pranzini, E. Geomorphological evolution and sedimentology of the Ombrone River Delta, Italy. J. Coast. Res. 1993, 9, 481–493. [Google Scholar]
  37. Bellotti, P.; Caputo, C.; Davoli, L.; Evangelista, S.; Garzanti, E.; Pugliese, F.; Valeri, P. Morpho-sedimentary characteristics and holocene evolution of the emergent part of the Ombrone River delta (southern Tuscany). Geomorphology 2004, 61, 71–90. [Google Scholar] [CrossRef]
  38. Barsanti, D.; Rombai, L. La Guerra delle Acque in Toscana; Edizioni Medicea: Firenze, Italy, 1986; p. 169. [Google Scholar]
  39. Mori, A. L’ evoluzione della costa grossetana dal Pliocene ad oggi. In Annuario 1932–1935 del Regio Liceo di Grosseto; La Maremmana: Grosseto, Italy, 1935; p. 21. [Google Scholar]
  40. Crist, E.P.; Cicone, R.C. A physically based transformation of thematic mapper data. IEEE Trans. Geosci. Remote Sens. 1984, GE-22, 256–263. [Google Scholar] [CrossRef]
  41. Pitty, A.F. The Nature of Geomorphology; Methuen & Co.: London, UK, 1982; p. 161. [Google Scholar]
  42. Luijendijk, A.; Hagenaars, G.; Ranasinghe, R.; Baart, F.; Donchyts, G.; Aarninkhof, S. The state of the world’s beaches. Nat. Sci. Rep. 2018, 8, 6641. [Google Scholar] [CrossRef]
Figure 1. The Ombrone River delta and the position of the two topographic profiles drawn in Figure 3 (Google Earth image).
Figure 1. The Ombrone River delta and the position of the two topographic profiles drawn in Figure 3 (Google Earth image).
Water 11 02281 g001
Figure 2. Ombrone River delta, beach ridges pattern mapped through the combination of different data sets and elaboration methodologies.
Figure 2. Ombrone River delta, beach ridges pattern mapped through the combination of different data sets and elaboration methodologies.
Water 11 02281 g002
Figure 3. North (A) and south (B) profiles from the early 15th (1) century until the limit of the progradation in middle 19th century (9). Due to the difference in deposition rates, horizontal scale cannot be considered a time scale.
Figure 3. North (A) and south (B) profiles from the early 15th (1) century until the limit of the progradation in middle 19th century (9). Due to the difference in deposition rates, horizontal scale cannot be considered a time scale.
Water 11 02281 g003
Figure 4. Real beach ridge 1625 ca., correlated to an erosional period, with the four best fitting mathematical functions.
Figure 4. Real beach ridge 1625 ca., correlated to an erosional period, with the four best fitting mathematical functions.
Water 11 02281 g004
Figure 5. Real beach ridge 1750 ca., correlated to a progradation period, with the four best fitting mathematical functions.
Figure 5. Real beach ridge 1750 ca., correlated to a progradation period, with the four best fitting mathematical functions.
Water 11 02281 g005
Figure 6. Beach ridges reconstruction for main historical phases from the 15th to the 19th century.
Figure 6. Beach ridges reconstruction for main historical phases from the 15th to the 19th century.
Water 11 02281 g006
Table 1. Relation between beach ridge and best four fitting mathematical functions. ER = apex erosion, PR = apex progradation. In bold, the best standard error for each ridge.
Table 1. Relation between beach ridge and best four fitting mathematical functions. ER = apex erosion, PR = apex progradation. In bold, the best standard error for each ridge.
Year (ca.)Evolution TypeStandard Error [m] for Four of the Best Different Fitting Models
MMF Weibull.Hoerl3rd Degree Polynomial
1500ER5.13.24.03.3
1550PR14.416.616.422.2
1600ER16.612.614.519.2
1615PR7.011.415.827.5
1625ER14.312.816.823.9
1650ER12.58.27.216.3
1660PR5.59.312.318.9
1675ER16.610.011.912.6
1700ER15.27.19.48.4
1725ER13.48.510.218.2
1750PR9.312.121.216.2
1775ER18.814.719.632.7
1800ER16.411.513.19.5
1820ER16.55.618.721.1
1830PR6.414.915.523.8
1850ER29.713.220.711.8

Share and Cite

MDPI and ACS Style

Mammì, I.; Rossi, L.; Pranzini, E. Mathematical Reconstruction of Eroded Beach Ridges at the Ombrone River Delta. Water 2019, 11, 2281. https://doi.org/10.3390/w11112281

AMA Style

Mammì I, Rossi L, Pranzini E. Mathematical Reconstruction of Eroded Beach Ridges at the Ombrone River Delta. Water. 2019; 11(11):2281. https://doi.org/10.3390/w11112281

Chicago/Turabian Style

Mammì, Irene, Lorenzo Rossi, and Enzo Pranzini. 2019. "Mathematical Reconstruction of Eroded Beach Ridges at the Ombrone River Delta" Water 11, no. 11: 2281. https://doi.org/10.3390/w11112281

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