Analisi ad alta risoluzione spazio-temporale di eventi grandinigeni nel bacino del Mediterraneo mediante una nuova tecnica satellitare ibrida multisensore

Negli ultimi anni si è sviluppato un crescente interesse per lo studio e il monitoraggio di fenomeni meteorologici estremi. In particolare, stanno aumentando osservazioni e testimonianze relative a manifestazioni temporalesche molto violente, spesso caratterizzate da grandine di notevole dimensione che rappresenta la causa principale di pericolo per le persone e di danni anche ingenti al comparto agricolo, infrastrutturale e alle automobili. Per questi motivi, in questo lavoro di tesi è stata presentata una nuova tecnica meteorologica denominata “Multi-sensor Approach for Satellite Hail Advection” (MASHA) che si basa sulla combinazione delle misure fornite dai satelliti geostazionari per trasportare l’informazione e da quelli in orbita polare per riconoscere le nubi grandinigene. I risultati della tecnica MASHA forniscono ogni cinque minuti su tutto il bacino del Mediterraneo la probabilità che la nube osservata abbia formato della grandine in un intervallo compreso tra 0 e 1. Le potenzialità del metodo MASHA riguardano soprattutto l’elevata risoluzione spaziale e temporale con la quale vengono riprodotti i sistemi grandinigeni, la capacità del modello di ricostruire i campi di grandine e l’evoluzione di questi sistemi anche in assenza di passaggi recenti dei satelliti polari. Per migliorare le performance della tecnica, è stato sviluppato un nuovo indice di severità chiamato Hail Severity Index (HSI) in grado di migliorare la localizzazione e la stima dell’intensità dei sistemi grandinigeni attraverso l’integrazione di nuove variabili: l’overshooting index, i fulmini e i dati radar. Le performance della tecnica MASHA e le potenzialità dell’indice HSI mostrate in questo lavoro di tesi, suggeriscono implementazioni operative a supporto di sistemi per il nowcasting di eventi temporaleschi particolarmente severi. Questo approccio, potrebbe essere molto efficace per fornire informazioni in tempo reale ai servizi meteorologici e nel contesto della protezione civile.

A comparison between vertical mixing parameterizations in the Levantine Sea

The goal of this thesis was the study of an optimal vertical mixing parameterization scheme in a mesoscale dominated field characterized from a strong vorticity and the presence of a layer of colder, less saline water at about 100 m depth (Atlantic Waters); in these conditions we compared six different experiments, that differ by the turbulent closure schemes, the presence or not of an enhanced diffusion parameterization and the presence or not of a double diffusion mixing parameterization. To evaluate the performance of the experiments and the model we compared the simulations with the ARGO observations of temperature and salinity available in our domain, in our period of interest. The conclusions were the following: • the increase of the resolution gives better results in terms of temperature in all the considered cases, and in terms of salinity. • The comparisons between the Pacanovski-Philander and the TKE turbulent closure schemes don’t show significant differences when the simulations are compared to the observations. • The removing of the enhanced diffusion parameterization in presence of the TKE turbulent closure submodel doesn’t give positive results, and show limitations in the resolving of gravitational instabilities near the surface • The k-ϵ turbulent closure model utilized in all the GLS experiments, is the best performing closure model among the three considered, with positive results in all the salinity comparison with the in situ observation and in most of the temperature comparisons. • The double mixing parameterization utilized in the k-ϵ closure submodel improves the results of the experiments improving both the temperature and salinity in comparison with the ARGO data.

Submesoscale features in regions with different baroclinic deformation radius

In questo studio, un modello nidificato "child" ad alta risoluzione (risoluzione 1/48°) è ottenuto attraverso la piattaforma SURF. Il modello "child" è ottenuto tramite "downscaling" dei campi medi giornalieri dal modello globale "parent" (1/12°). Questo permette di osservare le caratteristiche della sottomesoscala in due regioni dell'oceano Atlantico settentrionale, Azzorre e Bermuda, dal 4 al 12 gennaio 2021. Questa tesi si propone di condurre un'analisi preliminare della relazione fra il raggio di deformazione baroclino e l'esordio dell'attività di sottomesoscala, nelle regioni di interesse. A questo scopo, vengono effettuati molti confronti fra i campi risultanti dal "parent" e dal "child". In particolare, ci concentriamo sulla Mixed-Layer Instability (MLI) studiata attraverso variabili come la Mixed-Layer Depth (MLD), la vorticità relativa, le velocità orizzontali e verticali, l'energia cinetica e la frequenza di Brunt-Vaisala. Dai risultati, mentre il modello "parent" sembra inadeguato, quello "child" è in grado di rilevare la MLI e la presenza di filamenti e scie topografiche in entrambe le regioni, riproducendo meglio le correnti orizzontali e verticali alla sottomesoscala. Ciò fornisce una connessione tra mesoscala e sottomesoscala: mentre nelle Azzorre la MLI sembra svolgere un ruolo chiave nella ristratificazione della colonna d'acqua, lo stesso non sembra nelle Bermuda dove i vortici di mesoscala sono più ampi, influendo maggiormente sulla stratificazione verticale. Qui, i processi di ristratificazione sembrerebbero correlati alle mesoscale piuttosto che alle sottomesoscale. La MLI giocherebbe comunque un ruolo nella generazione di correnti alla sottomesoscala, insieme ad altri processi come la frontogenesi e le scie topografiche. In conclusione, il valore del raggio di deformazione baroclino non sembra influenzare l'attivazione dei processi di sottomesoscala, mentre sembrerebbe determinare l'importanza della MLI nel processo di ristratificazione.

Analisi dei tornado del 19 settembre 2021 in Pianura Padana tramite simulazioni numeriche ad altissima risoluzione

Il 19 settembre 2021 si sono formati in poche ore 7 tornado tra Lombardia ed Emilia Romagna, 4 dei quali di intensità F2. Una simile sequenza di tornado in Pianura Padana non ha precedenti. Il caso studio è stato analizzato sia dal punto di vista delle osservazioni che attraverso simulazioni numeriche ad altissima risoluzione con la catena modellistica BOLAM-MOLOCH, allo scopo di comprendere i meccanismi fisici che hanno portato alla genesi delle supercelle tornadiche. E' stata verificata la capacità del modello MOLOCH di simulare accuratamente lo sviluppo delle supercelle, alla luce dei modelli teorici in letteratura. Tramite le osservazioni al suolo, è stato osservato come durante l'evento sulla Pianura Padana si sia creato un "punto triplo", ossia una convergenza tra tre masse d'aria con caratteristiche diverse, nei pressi di dove si sono verificati i tornado. Per l'innesco della convezione è stato fondamentale il contributo di un cut-off low in quota, associato a un fronte freddo al suolo sulla pianura lombarda. Inoltre, l'avvezione di aria umida e instabile dall'Adriatico verso i settori di innesco, in prossimità della superficie, ha rivestito un ruolo rilevante. In Emilia invece è stata probabilmente determinante una convergenza tra il flusso instabile di origine adriatica e il flusso sud-occidentale sceso dagli Appennini. Per quanto riguarda l'intensificazione e il successivo sviluppo delle supercelle e dei tornado, sono state individuate due sorgenti principali di vorticità orizzontale: la vorticità baroclina che si è generata lungo il fronte freddo in Lombardia e la vorticità generata dallo shear verticale in presenza di un low-level jet. La circolazione nei bassi strati atmosferici ha trasportato la vorticità generata da entrambe queste sorgenti verso i sistemi convettivi in Lombardia, dove il processo di tilting ha trasformato la vorticità orizzontale in vorticità verticale ciclonica.

Geodetic strain rate and seismicity rates along the Apennines

Gli Appennini sono una regione tettonicamente attiva che è stata luogo di forti terremoti storici in passato, tra i quali il recente L’Aquila (2009) e le sequenze di Amatrice-Visso-Norcia (2016-2017). Una densità di stazioni GPS sempre maggiore permette di stimare il tasso di deformazione orizzontale con più risoluzione spaziale che in passato. In questa tesi studio la correlazione tra tassi di deformazione geodetica orizzontale e tassi di sismicità negli Appennini, esaminando possibili implicazioni per valutazioni del potenziale sismogenico. Uso un nuovo dataset di velocità GPS analizzando lo stato di deformazione attraverso il calcolo dello strain rate geodetico. Per questo calcolo utilizzo tre algoritmi, tra quelli presenti in letteratura scientifica, e discuto punti di forza e limitazioni di ciascuno. Seleziono poi una mappa di strain rate tra le altre e sviluppo un’analisi seguendo l’approccio di uno studio recente per investigare la relazione tra tassi di sismicità e tassi di deformazione lungo la catena Appenninica. Per farlo, verifico se le ipotesi che sostengono l’approccio sono verificate. Calcolo la sismicità di background a partire da un nuovo catalogo sismico per l’Italia e testo se gli eventi indipendenti seguono un processo di Poisson nel tempo. Infine, divido la regione di studio in aree a diversa deformazione sulla base dei tassi di deformazionee analizzado la distribuzione della sismicità di background in questi eventi e analizzando la relazione tra tassi di deformazione e tassi sismicità. Mostro che vi è una relazione è approssimativamente lineare ma con una dipendenza importante dalla risoluzione della mappa di strain rate. Descrivo poi alcune implicazioni che questa relazione avrebbe su stime di due quantità che concorrono a definire il potenziale sismogenico analizzando quali sono i parametri all’interno dell’approccio che possono influenzare maggiormente dette stime.

The routes of the Mediterranean Sea circulation

The Mediterranean Sea is a semi-enclosed basin connected to the Atlantic Ocean through the narrow and shallow Strait of Gibraltar and further subdivided in two different sub-basins, the Eastern Mediterranean and the Western Mediterranean, connected through the Stait of Sicily. On annual basis, a net heat budget of −7 W/m2, combined with exceeding evaporation over precipation and runoff together with wind stress, is responsible for the antiestuarine character of the zonal thermoaline circulation. The outflow at Gibraltar Strait is mainly composed of Levantine Intermediate Water (LIW) and deep water masses formed in the Western Mediterranean Sea. The aim of this thesis is to validate and quantitatively assess the main routes of water masses composing the ouflow at Gibraltar Strait, using for the first time in the Mediterranean Sea a lagrangian interpretation of the eulerian velocity field produced from an eddy-resolving reanalysis dataset, spanning from 2000 to 2012. A lagrangian model named Ariane is used to map out three-dimensional trajectories in order to describe the pathways of water mass transport from the Strait of Sicily, the Gulf of Lyon and the Northern Tyrrhenian Sea to the Gibraltar Strait. Numerical experiments were carried out by seeding millions of particles in the Strait of Gibraltar and following them backwards in time to track the origins of water masses and transport exchanged between the different sections of the Mediterranean. Finally, the main routes of the intermediate and deep water masses are reconstructed from virtual particles trajectories, which highlight the role of the Western Mediterranean Deep Water (WMDW) as the main contributor to the Gibraltar Strait outflow. For the first time, the quantitative description of the flow of water masses coming from the Eastern Mediterranean towards the Gibraltar Strait is provided and a new route that directly links the Northern Tyrrhenian Sea to Gibralatr Strait has been detected.

Analisi termodinamica di eventi meteorologici estremi nell'Artico

Il sistema terra è soggetto ad un aumento di temperatura causato in larga parte dall'intervento umano. Il riscaldamento nell'Artico avviene ad un velocità doppia rispetto a quella del resto del globo e induce lo scioglimento dei ghiacci polari oceanici e terrestri. Lo studio della variabilità meteorologica delle regioni polari e subpolari può aiutarci a comprendere il ruolo di eventi meteorologici estremi nel determinare forti variazioni di temperatura superficiale della calotta polare. Questo lavoro di tesi discute la fenomenologia degli eventi di blocco meteorologico, regime di circolazione anticiclonica persistente, nella troposfera delle alte latitudini ed il loro legame con estremi di temperatura superficiale nelle regioni artiche. Successivamente, una metodologia di analisi termodinamica per i suddetti è descritta ed applicata ad un caso studio descritto nella letteratura scientifica un blocco alle alte latitudini coincide con la propagazione di una perturbazione ciclonica dalle medie latitudini fino alla calotta artica. Questa particolare configurazione sinottica ha prodotto una variazione della temperatura della calotta artica di circa 7 K in 4 giorni. Sulla base dell’analisi termodinamica è possibile evidenziare una combinazione di processi fisici, diabatici e adiabatici, che concorrono nel determinare il repentino riscaldamento dell’Artico. L’applicazione di simili metodologie su ulteriori casi studio potrebbe contribuire alla definizione di un paradigma dinamico per il fenomeno delle ondate di calore artiche.

Characterization of ENSO changes under idealised radiative forcing in an intermediate complexity general circulation model

Intermediate-complexity general circulation models are a fundamental tool to investigate the role of internal and external variability within the general circulation of the atmosphere and ocean. The model used in this thesis is an intermediate complexity atmospheric general circulation model (SPEEDY) coupled to a state-of-the-art modelling framework for the ocean (NEMO). We assess to which extent the model allows a realistic simulation of the most prominent natural mode of variability at interannual time scales: El-Niño Southern Oscillation (ENSO). To a good approximation, the model represents the ENSO-induced Sea Surface Temperature (SST) pattern in the equatorial Pacific, despite a cold tongue-like bias. The model underestimates (overestimates) the typical ENSO spatial variability during the winter (summer) seasons. The mid-latitude response to ENSO reveals that the typical poleward stationary Rossby wave train is reasonably well represented. The spectral decomposition of ENSO features a spectrum that lacks periodicity at high frequencies and is overly periodic at interannual timescales. We then implemented an idealised transient mean state change in the SPEEDY model. A warmer climate is simulated by an alteration of the parametrized radiative fluxes that corresponds to doubled carbon dioxide absorptivity. Results indicate that the globally averaged surface air temperature increases of 0.76 K. Regionally, the induced signal on the SST field features a significant warming over the central-western Pacific and an El-Niño-like warming in the subtropics. In general, the model features a weakening of the tropical Walker circulation and a poleward expansion of the local Hadley cell. This response is also detected in a poleward rearrangement of the tropical convective rainfall pattern. The model setting that has been here implemented provides a valid theoretical support for future studies on climate sensitivity and forced modes of variability under mean state changes.

Caratterizzazione delle formazioni cicloniche associate a precipitazioni estreme sulla Sardegna

Gli eventi degli ultimi due decenni hanno mostrato come la Sardegna, in particolare il suo settore est, sia soggetta a frequenti eventi di precipitazione intensa che comportano perdita di vite umane e ingenti danni. Le alluvioni avvenute a Olbia nel 2013 o la più recente avvenuta a Bitti nel 2020 sono la dimostrazione dell’impatto che questi eventi possono avere sul territorio. Dal punto di vista della pratica previsionale di questi fenomeni intensi, è in uso presso ARPAS una classificazione dei cicloni avvenuti nel periodo compreso tra il 1957 e il 1999, che presentano configurazioni bariche simili e raggruppati in gruppi o “cluster”. Lo scopo di questo lavoro di tesi è estendere questa classificazione per il periodo 1951-1957 e 2000-2009, sulla base della configurazione sinottica in quota e al suolo in cui si sviluppano i cicloni appartenenti ai vari cluster. Tramite questa nuova classificazione, e utilizzando la temperatura superficiale media del mare, verrà analizzata la loro stagionalità e la loro durata temporale. Infine, sarà ricostruita la precipitazione legata ai diversi cluster suddividendo la precipitazione per i settori est e ovest dell’isola, per valutarne distribuzione e variabilità spaziale. I risultati mostrano come vi siano differenze tra i cicloni associati a ciascun cluster in base alle diverse configurazioni sinottiche. Un fattore importante è dato anche dalla stagione in cui i cicloni si verificano: emerge una frequenza maggiore nei mesi autunnali in cui le temperature medie mensili e l’anomalia di temperatura media risulta più alta. La distribuzione spaziale delle precipitazioni associate ai cluster conferma che le precipitazioni più intense e omogene interessano il settore est dell’isola. Poiché lo studio è basato interamente su una classificazione di tipo operativo in uso ad ARPAS, i suoi sviluppi riguardano soprattutto un maggior supporto e conoscenza previsionale degli eventi estremi sulla Sardegna.

Development of a relationship between physical parameters of tsunami waves and tsunami intensity based on macroscopic effects and its application to the Italian tsunami database

Tsunamis are rare events. However, their impact can be devastating and it may extend to large geographical areas. For low-probability high-impact events like tsunamis, it is crucial to implement all possible actions to mitigate the risk. The tsunami hazard assessment is the result of a scientific process that integrates traditional geological methods, numerical modelling and the analysis of tsunami sources and historical records. For this reason, analysing past events and understanding how they interacted with the land is the only way to inform tsunami source and propagation models, and quantitatively test forecast models like hazard analyses. The primary objective of this thesis is to establish an explicit relationship between the macroscopic intensity, derived from historical descriptions, and the quantitative physical parameters measuring tsunami waves. This is done first by defining an approximate estimation method based on a simplified 1D physical onshore propagation model to convert the available observations into one reference physical metric. Wave height at the coast was chosen as the reference due to its stability and independence of inland effects. This method was then implemented for a set of well-known past events to build a homogeneous dataset with both macroseismic intensity and wave height. By performing an orthogonal regression, a direct and invertible empirical relationship could be established between the two parameters, accounting for their relevant uncertainties. The target relationship is extensively tested and finally applied to the Italian Tsunami Effect Database (ITED), providing a homogeneous estimation of the wave height for all existing tsunami observations in Italy. This provides the opportunity for meaningful comparison for models and simulations, as well as quantitatively testing tsunami hazard models for the Italian coasts and informing tsunami risk management initiatives.

Validation of an intermediate-complexity flood risk model for climate change adaptation

In questo lavoro di tesi viene presentato e validato un modello di rischio di alluvione a complessità intermedia per scenari climatici futuri. Questo modello appartiene a quella categoria di strumenti che mirano a soddisfare le esigenze identificate dal World Climate Research Program (WRCP) per affrontare gli effetti del cambiamento climatico. L'obiettivo perseguito è quello di sviluppare, seguendo un approccio ``bottom-up" al rischio climatico regionale, strumenti che possano aiutare i decisori a realizzare l'adattamento ai cambiamenti climatici. Il modello qui presentato è interamente basato su dati open-source forniti dai servizi Copernicus. Il contributo di questo lavoro di tesi riguarda lo sviluppo di un modello, formulato da (Ruggieri et al.), per stimare i danni di eventi alluvionali fluviali per specifici i livelli di riscaldamento globale (GWL). Il modello è stato testato su tre bacini idrografici di medie dimensioni in Emilia-Romagna, Panaro, Reno e Secchia. In questo lavoro, il modello viene sottoposto a test di sensibilità rispetto a un'ipotesi enunciata nella formulazione del modello, poi vengono effettuate analisi relative all'ensemble multi-modello utilizzato per le proiezioni. Il modello viene quindi validato, confrontando i danni stimati nel clima attuale per i tre fiumi con i danni osservati e confrontando le portate simulate con quelle osservate. Infine, vengono stimati i danni associati agli eventi alluvionali in tre scenari climatici futuri caratterizzati da GWL di 1.5° C, 2.0° C e 3.0°C.

Impact of convective motions on pollutant dispersion in urban canyons: an analysis through RANS simulations

Since the majority of the population of the world lives in cities and that this number is expected to increase in the next years, one of the biggest challenges of the research is the determination of the risk deriving from high temperatures experienced in urban areas, together with improving responses to climate-related disasters, for example by introducing in the urban context vegetation or built infrastructures that can improve the air quality. In this work, we will investigate how different setups of the boundary and initial conditions set on an urban canyon generate different patterns of the dispersion of a pollutant. To do so we will exploit the low computational cost of Reynolds-Averaged Navier-Stokes (RANS) simulations to reproduce the dynamics of an infinite array of two-dimensional square urban canyons. A pollutant is released at the street level to mimic the presence of traffic. RANS simulations are run using the k-ɛ closure model and vertical profiles of significant variables of the urban canyon, namely the velocity, the turbulent kinetic energy, and the concentration, are represented. This is done using the open-source software OpenFOAM and modifying the standard solver simpleFoam to include the concentration equation and the temperature by introducing a buoyancy term in the governing equations. The results of the simulation are validated with experimental results and products of Large-Eddy Simulations (LES) from previous works showing that the simulation is able to reproduce all the quantities under examination with satisfactory accuracy. Moreover, this comparison shows that despite LES are known to be more accurate albeit more expensive, RANS simulations represent a reliable tool if a smaller computational cost is needed. Overall, this work exploits the low computational cost of RANS simulations to produce multiple scenarios useful to evaluate how the dispersion of a pollutant changes by a modification of key variables, such as the temperature.

Simulazioni Software-In-the-Loop per la verifica funzionale di un sistema di tracking satellitare (ALMATracker)

L'ALMATracker è un sistema di puntamento per la stazione di terra di ALMASat-1. La sua configurazione non segue la classica Azimuth-Elevazione, bensì utilizza gli assi α-β per evitare punti di singolarità nelle posizioni vicino allo zenit. Ancora in fase di progettazione, utilizzando in congiunta SolidWorks e LabVIEW si è creato un Software-in-the-loop per la sua verifica funzionale, grazie all'utilizzo del relativamente nuovo pacchetto NI Softmotion. Data la scarsa esperienza e documentazione che si hanno su questo recente tool, si è prima creato un Case Study che simulasse un sistema di coordinate cilindriche in modo da acquisire competenza. I risultati conseguiti sono poi stati sfruttati per la creazione di un SIL per la simulazione del movimento dell'ALMATracker. L'utilizzo di questa metodologia di progettazione non solo ha confermato la validità del design proposto, ma anche evidenziato i problemi e le potenzialità che caratterizzano questo pacchetto software dandone un analisi approfondita.

Sviluppo di un sistema tomografico con tubo a raggi x da 300 kv: progettazione e test preliminari.

In questo elaborato sono descritti la progettazione e lo sviluppo di un sistema tomografico con sorgente di raggi X da 300 kV per analisi nel campo dei Beni Culturali; in particolare sono stati effettuati test preliminari sia per verificare le prestazioni del tubo a raggi X e dello schermo scintillatore sia per avere indicazioni utili alla progettazione (e conseguente costruzione) del box di rivelazione. Il collaudo della sorgente di raggi X è stato eseguito con lo studio dell'andamento della dose integrale e/o del rateo di dose in funzione della tensione, della corrente e del tempo d'irraggiamento. La caratterizzazione dello schermo scintillatore di CsI(Tl) è basata sul confronto delle sue caratteristiche con uno schermo analogo di spessore differente nelle stesse condizioni sperimentali, attraverso lo studio dei parametri standard dei dispositivi per imaging, quali la curva di risposta, il rapporto segnale-rumore e la risoluzione spaziale. Inoltre, come ulteriore test, lo schermo scintillatore è stato utilizzato per l’indagine tomografica di alcuni pani di terra contenenti olle di epoca villanoviana. Per la progettazione del box rivelatore sono state invece svolte misure di tipo ottico. I vari test svolti in laboratorio hanno dato risultati positivi e hanno dimostrato la validità dei componenti scelti per il nuovo sistema tomografico, il cui assemblaggio verrà completato a breve, essendo stata ultimata proprio in queste settimane la costruzione del box di rivelazione.

Rilevati arginali in terra rinforzata con geosintetici: il caso della Vena Mazzarini a Cesenatico (FC)

Oggetto della presente tesi di laurea è quello di analizzare sia dal punto di vista teorico che da quello pratico l’impiego di geosintetici di rinforzo per il consolidamento di rilevati e sponde arginali al fine di ottenere un sistema pluricomposto terreno-rinforzo che manifesti adeguati margini di sicurezza nei confronti delle verifiche di stabilità prescritte dalla normativa vigente, soprattutto nei confronti dell’azione indotta dal sisma. In particolare il lavoro è suddiviso di due macro parti fondamentali: una prima parte dedicata allo studio degli aspetti teorici legati alla caratterizzazione dei materiali di rinforzo con particolare riferimento ai geotessili ed alle geogriglie, ai principi fondamentali di funzionamento del sistema terreno-rinforzi, alla normativa di riferimento per tali costruzioni, al dimensionamento dei pendii in terra rinforzata ed al loro inserimento ambientale; la seconda parte, basata sullo studio di un caso di reale applicazione del metodo di rinforzo con geosintetici, si svolge attraverso una approfondita e critica analisi della relazione geologica messa a disposizione dal Comune di Cesenatico (FC) con lo scopo di elaborare le prove svolte in situ, ricavare i parametri geotecnici caratteristici del terreno e definire un modello geologico e geotecnico di riferimento per le successive elaborazioni. Il progetto prevede infatti il consolidamento delle sponde del canale che attraversa, parallelamente alla linea di costa, la città di Cesenatico (FC) e denominato “Vena Mazzarini”, dal nome dell’Ingegnere che ne curò il progetto originale. In particolare si è previsto l’impiego di geogriglie di rinforzo in HDPE e tutte le verifiche sono state realizzate mediante il software di calcolo Geostru Slope® e Slope/M.R.E® messi cordialmente a disposizione dalla software house Geostru S.r.l.