Hydrological-hydraulic modelling of the bench terraces in hilly and mountain areas

Terrace building have been expanded in the 19th century because of the increased demographic pressure and the need to crop additional areas at steeper slopes. Terraces are also important to regulate the hydrological behavior of the hillslope. Bench terraces, reducing the terrain slope and the length of the overland flow, quantitatively control the runoff flow velocity, facilitating the drainage and thus leading to a reduction of soil erosion. The study of the hydrologic-hydraulic function of terraced slopes is essential in order to evaluate their possible use to cooperate for flood-risk mitigation also preserving the landscape value. Few studies in literature are available on rainfall-runoff transformation and flood risk mitigation in terrace areas. Then, research results in this field are still scarce. The goal of this work is to improve knowledge on hydrological processes affecting a terraced slope and their effect on flood control. Specific researches objectives are: • Studing the the reduction of peak runoff at the toe of a hillslope and the delay in the passage of peak flow, which are provided by sequence of dry-stone walls under different space arrangements along the hillslope; • Understanding the rainfall-runoff separation mechanism and the superficial and subsurface flow propagation in case of terraced slopes. In order to reach the above objectives the hydrological response of a bench terrace was investigated by using a research approache based on modelling and experimental activities. In the first part of the thesis the The FLO-2D model is used to analyse the runoff propagation mechanism of a terraced slope (sequence of dry-stone walls) by varying number and spacing of terraces and assuming two hydrological soil setting scenarios in terms of antecedent moisture conditions within the Soil Conservation Service-Curve Number method. The model analysis shows that the majority of runoff modifications at the outlet of a terraced system result from topographical modifications rather than local variations of the infiltration capacity at the dry-stone wall zone. Repeated modelling applications show that, given a quite-typical scenario of a 20°-sloped hillslope and a reference intense rainstorm, the peak discharge reduction at the hillslope outlet depends on the percentage of the area managed with terraces. The reduction can be calculated with a logarithmic-type function (for example, an increase of terraced area from 10% to 30% might bring to runoff peak reduction of almost 45%). This information can help determine where terrace additions are more effective in terms of hydrological benefit. The second part of the thesis focus on an experimental/modelling research that aims to better focus the times of the hydrological response, which are determined by a hillslope plot bounded by a dry-stone wall, considering both the overland flow and the groundwater. A physical model, characterized by a quasi-real scale, has been built to reproduce the behavior of a 3%, 6% and 9% outward sloped terrace at bare and vegetated soil condition.The model consists of a steel metal box (1 m large, 3.3 m long, 0.8 m high) containing the hillslope terrain. The terrain is equipped with two piezometers, 9 TDR sensors measuring the volumetric water content, a surface spillway at the head releasing the steady discharge under test, two scales one at the wall base to measure the groundwater discharge and another at the top of the wall to measure the surface runoff. The experiments deal with different initial moisture condition (high and low degree of saturation), and discharges of 19.5, 12.0 and 5.0 l/min. Each experiment has been replicated, conducting a total number of 35 tests. The volumetric water content analysis produced by the 9 TDR sensors was able to provide a quite satisfactory representation of the soil moisture during the runs. Then, different lag times at the outlet since the inflow initiation were measured both for runoff and groundwater. Moreover, the time of depletion and the piezometer response have been monitored and analyzed, well corroborating the findings on the kinematics of the terrace plot. Finally, the computation of the specific Curve Number (Soil Conservation Service) of the physical model has revealed values rather large if compared with those reported in the literature. This phenomenon was likely caused by the high values of the inflow discharge, the limited cross-width of the model (1 m), the increasing compactness of the soil owing to the experiment repetition and the confined waterproof box). The experimental results indicate that terrace soil was highly heterogeneous, including discontinuities and piping systems that facilitated a rapid infiltration and the development of fast subsurface flow. The Groundwater in general is a small part of the total outflow but in case the presence of pipe is important it is coupled with impulsesive infiltration rates. A conceptual hydrological model was implemented and calibrated based on the experimental data. The model results fit well the measurements even if the groundwater component is not properly modelled. This is due to the activation of important piping systems during some of the tests; the physical proces that describ this located losses were not studied and integrated in the model. These pioneering experiments have produced some remarkable outcomes on the important role of lag-times (runoff and groundwater) and provided new knowledgement on the hydrological functioning of bench terraced systems for addressing more efficient management and maintenance issues of this important agricoltural structures.

I sistemi terrazzati si sono diffusi nel diciannovesimo secolo a seguito della crescente pressione demografica e della conseguente necessità di estendere le coltivazioni anche su terreni ad elevata pendenza. Oltre che dal punto di vista agrario tali sistemi sono importanti ai fini della regolazione della risposta idrogeologica di un versante. Infatti essi riducono la pendenza e la lunghezza dello scorrimento superficiale, controllando quindi quantitativamente la velocità del deflusso superficiale, facilitando il drenaggio e contribuendo in questo modo alla riduzione dei fenomeni erosivi. Lo studio della funzione idrologico-idraulica dei versanti terrazzati è essenziale per valutarne il possibile utilizzo come misure di mitigazione del rischio idraulico capaci anche di preservare il valore paesaggistico dei territori su cui essi insistono. In letteratura sono disponibili pochi studi inerenti la risposta idrologica di versanti terrazzati; l’avanzamento della ricerca in tale ambito è l’obiettivo principale di questo lavoro. In particolare vengono affrontate le seguenti tematiche: - la valutazione degli effetti di mitigazione della pericolosità idraulica (riduzione del picco di piena e suo ritardo temporale) a scala di versante indotti dalla presenza di sistemi terrazzati; - lo studio dei meccanismi di trasformazione afflussi-deflussi e dei processi di propagazione degli stessi in sistemi terrazzati; Al fine di raggiungere tali obiettivi è stato implementato un approccio integrato basato su attività sperimentali e modellistiche. Nella prima parte del lavoro è stato utilizzato il modello idraulico FLO-2D per analizzare i processi di propagazione in atto in un versante terrazzato (composto da una sequenza di muri a secco) al variare del numero e della disposizione spaziale dei terrazzi e assumendo due diversi scenari di saturazione del suolo rappresentati da diversi valori di umidità iniziale antecedente l'evento, come previsto dal metodo Soil Conservation Service - Curve Number. L'analisi modellistica mostra che la riduzione del deflusso alla base del sistema terrazzato dipende maggiormente dalle modifiche topografiche piuttosto che dalle variazioni della capacità di infiltrazione del suolo adiacente il muro. Le simulazioni eseguite su di un versante con una pendenza di 20° e alimentato da un evento di precipitazione intensa, mostrano che il picco di piena alla sezione di chiusura si riduce in funzione della percentuale di area terrazzata. Tale riduzione può essere valutata attraverso una specifica funzione logaritmica (per esempio, al crescere dell'area terrazzata da 10% a 30% la riduzione del picco di piena può essere quasi del 45%). Questa informazione può aiutare a individuare il corretto inserimento dei terrazzi per una maggiore efficace in termini di benefici idrologici. La seconda parte del lavoro riguarda lo studio della risposta idrologica di un’unità terrazzata con un muro a secco attraverso attività sperimentali e modellistiche. In particolare è stato costruito un modello fisico a scala reale per riprodurre il comportamento di un terrazzo al variare della sua pendenza (3%, 6% e 9%) e del tipo di copertura del suolo (suolo nodo o inerbito). Il modello consiste in un box metallico (1 metro di larghezza, 3.3 metri di lunghezza e 0.8 m di altezza) che contiene al suo interno un terrazzo composto da un versante delimitato a valle da un muro a secco. Tale versante è stato strumentato con 9 sensori TDR per la misura del contenuto di umidità del suolo, uno sfioratore delle portate liquidi in ingresso al versante a monte dello stesso, due bilance per la misura del deflusso, una posizionata alla base del muro per la misura del deflusso sotterraneo e una in corrispondenza della parte superiore del muro per la misura del deflusso superficiale. Gli esperimenti sono stati caratterizzati da differenti condizioni di umidità iniziale (ad alto e basso grado di saturazione) e da portate liquide in ingresso costanti e pari a 19.5, 12 e 5 l/minuto. Ogni esperimento è stato replicato per un totale di 35 esperimenti eseguiti. L'esame delle misure dei 9 sensori TDR ha fornito una soddisfacente rappresentazione dell'andamento dell'umidità globale del suolo nel corso di ogni esperimento. Sono stati poi misurati diversi tempi caratteristici della risposta idrologica alla sezione di chiusura sia per il deflusso superficiale che per il deflusso sotterraneo. I risultati ottenuti aiutano a comprendere la cinematica dei processi idrologici che caratterizzano l’unità terrazzata. E’ stato calcolato uno specifico Curve Number (Soil Conservation Service) associato all’unità terrazzata che assume valori piuttosto alti se comparati a quelli riportati in letteratura. Questo comportamento è probabilmente legato alle alte portate in ingresso, alla limitata sezione idraulica (1m), alla crescente compattazione del suolo causata dal susseguirsi delle prove e al fatto che il terrazzo è confinato all'interno di una struttura metallica impermeabile. Un innovativo modello idrologico è stato implementato e calibrato sui dati sperimentali. I risultati modellistici riproducono in modo soddisfacente le misure soprattutto per quanto riguarda il deflusso superficiale che è la componente prevalente di deflusso. In generale il deflusso sotterraneo non risulta invece essere propriamente simulato in quanto il modello non tiene conto di particolari fenomeni di infiltrazione impulsiva presenti in alcune prove. Infatti, i risultati sperimentali indicano che il suolo all'interno del terrazzo è altamente eterogeneo, con la presenza di discontinuità e sistemi di canali sotterranei che facilitano una rapida infiltrazione e lo sviluppo di deflusso sub-superficiale impulsivo che va a sommarsi al deflusso profondo (generalmente di modesta entità) alimentato dall’infiltrazione connessa agli strati superficiali del suolo. La sperimentazione effettuata risulta innovativa e fornisce nuove conoscenze sulla funzione idrologica-idraulica di un sistema terrazzato che possono servire per indirizzare in modo più efficiente la gestione e la manutenzione di queste importanti sistemazioni agrarie.