Study of the circulation processes in the northern Adriatic sea - coastal area and Venice lagoon inlets

This work is a detailed study of hydrodynamic processes in a defined area, the littoral in front of the Venice Lagoon and its inlets, which are complex morphological areas of interconnection. A finite element hydrodynamic model of the Venice Lagoon and the Adriatic Sea has been developed in order to study the coastal current patterns and the exchanges at the inlets of the Venice Lagoon. This is the first work in this area that tries to model the interaction dynamics, running together a model for the lagoon and the Adriatic Sea. First the barotropic processes near the inlets of the Venice Lagoon have been studied. Data from more than ten tide gauges displaced in the Adriatic Sea have been used in the calibration of the simulated water levels. To validate the model results, empirical flux data measured by ADCP probes installed inside the inlets of Lido and Malamocco have been used and the exchanges through the three inlets of the Venice Lagoon have been analyzed. The comparison between modelled and measured fluxes at the inlets outlined the efficiency of the model to reproduce both tide and wind induced water exchanges between the sea and the lagoon. As a second step, also small scale processes around the inlets that connect the Venice lagoon with the Northern Adriatic Sea have been investigated by means of 3D simulations. Maps of vorticity have been produced, considering the influence of tidal flows and wind stress in the area. A sensitivity analysis has been carried out to define the importance of the advection and of the baroclinic pressure gradients in the development of vortical processes seen along the littoral close to the inlets. Finally a comparison with real data measurements, surface velocity data from HF Radar near the Venice inlets, has been performed, which allows for a better understanding of the processes and their seasonal dynamics. The results outline the predominance of wind and tidal forcing in the coastal area. Wind forcing acts mainly on the mean coastal current inducing its detachment offshore during Sirocco events and an increase of littoral currents during Bora events. The Bora action is more homogeneous on the whole coastal area whereas the Sirocco strengthens its impact in the South, near Chioggia inlet. Tidal forcing at the inlets is mainly barotropic. The sensitivity analysis shows how advection is the main physical process responsible for the persistent vortical structures present along the littoral between the Venice Lagoon inlets. The comparison with measurements from HF Radar not only permitted a validation the model results, but also a description of different patterns in specific periods of the year. The success of the 2D and the 3D simulations on the reproduction both of the SSE, inside and outside the Venice Lagoon, of the tidal flow, through the lagoon inlets, and of the small scale phenomena, occurring along the littoral, indicates that the finite element approach is the most suitable tool for the investigation of coastal processes. For the first time, as shown by the flux modeling, the physical processes that drive the interaction between the two basins were reproduced.

Climate impact studies of sediment transport on the Adriatic coastal zone

The study of the impact of climate change on the environment has been based, until very recently, on an global approach, whose interest from a local point of view is very limited. This thesis, on the contrary, has treated the study of the impact of climate change in the Adriatic Sea basin following a twofold strategy of regionalization and integration of numerical models in order to reproduce the present and future scenarios of the system through a more and more realistic and solid approach. In particular the focus of the study was on the impact on the physical environment and on the sediment transport in the basin. This latter is a very new and original issue, to our knowledge still uninvestigated. The study case of the coastal area of Montenegro was particularly studied, since it is characterized by an important supply of sediment through the Buna/Bojana river, second most important in the Adriatic basin in terms of flow. To do this, a methodology to introduce the tidal processes in a baroclinic primitive equations Ocean General Circulation Model was applied and tidal processes were successfully reproduced in the Adriatic Sea, analyzing also the impacts they have on the mean general circulation, on salt and heat transport and on mixing and stratification of the water column in the different seasons of the year. The new hydrodynamical model has been further coupled with a wave model and with a river and sea sediment transport model, showing good results in the reproduction of sediment transport processes. Finally this complex coupled platform was integrated in the period 2001-2030 under the A1B scenario of IPCC, and the impact of climate change on the physical system and on sediment transport was preliminarily evaluated.

The rupture process of recent tsunamigenic earthquakes by geophysical data inversion

Subduction zones are the favorite places to generate tsunamigenic earthquakes, where friction between oceanic and continental plates causes the occurrence of a strong seismicity. The topics and the methodologies discussed in this thesis are focussed to the understanding of the rupture process of the seismic sources of great earthquakes that generate tsunamis. The tsunamigenesis is controlled by several kinematical characteristic of the parent earthquake, as the focal mechanism, the depth of the rupture, the slip distribution along the fault area and by the mechanical properties of the source zone. Each of these factors plays a fundamental role in the tsunami generation. Therefore, inferring the source parameters of tsunamigenic earthquakes is crucial to understand the generation of the consequent tsunami and so to mitigate the risk along the coasts. The typical way to proceed when we want to gather information regarding the source process is to have recourse to the inversion of geophysical data that are available. Tsunami data, moreover, are useful to constrain the portion of the fault area that extends offshore, generally close to the trench that, on the contrary, other kinds of data are not able to constrain. In this thesis I have discussed the rupture process of some recent tsunamigenic events, as inferred by means of an inverse method. I have presented the 2003 Tokachi-Oki (Japan) earthquake (Mw 8.1). In this study the slip distribution on the fault has been inferred by inverting tsunami waveform, GPS, and bottom-pressure data. The joint inversion of tsunami and geodetic data has revealed a much better constrain for the slip distribution on the fault rather than the separate inversions of single datasets. Then we have studied the earthquake occurred on 2007 in southern Sumatra (Mw 8.4). By inverting several tsunami waveforms, both in the near and in the far field, we have determined the slip distribution and the mean rupture velocity along the causative fault. Since the largest patch of slip was concentrated on the deepest part of the fault, this is the likely reason for the small tsunami waves that followed the earthquake, pointing out how much the depth of the rupture plays a crucial role in controlling the tsunamigenesis. Finally, we have presented a new rupture model for the great 2004 Sumatra earthquake (Mw 9.2). We have performed the joint inversion of tsunami waveform, GPS and satellite altimetry data, to infer the slip distribution, the slip direction, and the rupture velocity on the fault. Furthermore, in this work we have presented a novel method to estimate, in a self-consistent way, the average rigidity of the source zone. The estimation of the source zone rigidity is important since it may play a significant role in the tsunami generation and, particularly for slow earthquakes, a low rigidity value is sometimes necessary to explain how a relatively low seismic moment earthquake may generate significant tsunamis; this latter point may be relevant for explaining the mechanics of the tsunami earthquakes, one of the open issues in present day seismology. The investigation of these tsunamigenic earthquakes has underlined the importance to use a joint inversion of different geophysical data to determine the rupture characteristics. The results shown here have important implications for the implementation of new tsunami warning systems – particularly in the near-field – the improvement of the current ones, and furthermore for the planning of the inundation maps for tsunami-hazard assessment along the coastal area.

Analisi del popolamento e delle dinamiche culturali del Neolitico nella Puglia Centrale : un approccio territoriale per la Bassa Murgia

L’obiettivo della ricerca condotta dal candidato è stato sin dal principio la proposta di un’analisi del popolamento e delle dinamiche culturali dell’età neolitica nella Puglia Centrale. Tale aspirazione derivava dal fatto che il territorio oggetto di tale ricerca è stato, negli ultimi decenni, assai meno investigato rispetto alle limitrofe aree, come la piana del Tavoliere ed il Salento in particolare. Questa mancanza di informazioni rendeva necessaria una rivisitazione delle ricerche e degli studi, attraverso la costruzione di un atlante del popolamento, al fine di colmare questo vuoto. Infatti, una valutazione del popolamento dell’età neolitica nell’Italia sud-orientale non può essere considerata per aree culturali distinte e geograficamente limitate. Nonostante vi siano importanti differenziazioni di facies archeologiche, è solo attraverso una visione complessiva e con un controllo articolato delle informazioni che è possibile descrivere i processi e le dinamiche che hanno portato ad un popolamento così imponente, dalla sua espansione al repentino collasso, causando la trasformazione del paesaggio pugliese. La scelta di guardare alla complessità neolitica di questo versante della Puglia si lega culturalmente, oltre che geograficamente, alla ricerca già effettuata sulla grande diffusione del neolitico nella piana del Tavoliere di Puglia, unico esempio europeo in quanto a colonizzazione capillare del territorio..Ciò ha permesso di disegnare un quadro globale, che al momento esclude solo la penisola salentina, del popolamento neolitico della Puglia settentrionale e centrale. Le caratteristiche che differenziano i due complessi, sono strettamente legate sia all’habitat insediativo, inteso come scelte insediamentali, che nella regione indagata appare organizzato su diversi livelli (area costiera adriatica, area prospiciente antichi alvei torrentizi, propaggini dell’altopiano delle Murge), che alle forme culturali, intese come produzione artigianale e soprattutto aspetti rituali connessi alla frequentazione delle grotte naturali.

Hydrodynamic induced by an array of wave energy converters. Experimental and numerical analysis

The thesis analyses the hydrodynamic induced by an array of Wave energy Converters (WECs), under an experimental and numerical point of view. WECs can be considered an innovative solution able to contribute to the green energy supply and –at the same time– to protect the rear coastal area under marine spatial planning considerations. This research activity essentially rises due to this combined concept. The WEC under exam is a floating device belonging to the Wave Activated Bodies (WAB) class. Experimental data were performed at Aalborg University in different scales and layouts, and the performance of the models was analysed under a variety of irregular wave attacks. The numerical simulations performed with the codes MIKE 21 BW and ANSYS-AQWA. Experimental results were also used to calibrate the numerical parameters and/or to directly been compared to numerical results, in order to extend the experimental database. Results of the research activity are summarized in terms of device performance and guidelines for a future wave farm installation. The device length should be “tuned” based on the local climate conditions. The wave transmission behind the devices is pretty high, suggesting that the tested layout should be considered as a module of a wave farm installation. Indications on the minimum inter-distance among the devices are provided. Furthermore, a CALM mooring system leads to lower wave transmission and also larger power production than a spread mooring. The two numerical codes have different potentialities. The hydrodynamics around single and multiple devices is obtained with MIKE 21 BW, while wave loads and motions for a single moored device are derived from ANSYS-AQWA. Combining the experimental and numerical it is suggested –for both coastal protection and energy production– to adopt a staggered layout, which will maximise the devices density and minimize the marine space required for the installation.

El Niño southern oscillation and its effect on fog oases along the Peruvian and Chilean coastal deserts

Fog oases, locally named Lomas, are distributed in a fragmented way along the western coast of Chile and Peru (South America) between ~6°S and 30°S following an altitudinal gradient determined by a fog layer. This fragmentation has been attributed to the hyper aridity of the desert. However, periodically climatic events influence the ‘normal seasonality’ of this ecosystem through a higher than average water input that triggers plant responses (e.g. primary productivity and phenology). The impact of the climatic oscillation may vary according to the season (wet/dry). This thesis evaluates the potential effect of climate oscillations, such as El Niño Southern Oscillation (ENSO), through the analysis of vegetation of this ecosystem following different approaches: Chapters two and three show the analysis of fog oasis along the Peruvian and Chilean deserts. The objectives are: 1) to explain the floristic connection of fog oases analysing their taxa composition differences and the phylogenetic affinities among them, 2) to explore the climate variables related to ENSO which likely affect fog production, and the responses of Lomas vegetation (composition, productivity, distribution) to climate patterns during ENSO events. Chapters four and five describe a fog-oasis in southern Peru during the 2008-2010 period. The objectives are: 3) to describe and create a new vegetation map of the Lomas vegetation using remote sensing analysis supported by field survey data, and 4) to identify the vegetation change during the dry season. The first part of our results show that: 1) there are three significantly different groups of Lomas (Northern Peru, Southern Peru, and Chile) with a significant phylogenetic divergence among them. The species composition reveals a latitudinal gradient of plant assemblages. The species origin, growth-forms typologies, and geographic position also reinforce the differences among groups. 2) Contradictory results have emerged from studies of low-cloud anomalies and the fog-collection during El Niño (EN). EN increases water availability in fog oases when fog should be less frequent due to the reduction of low-clouds amount and stratocumulus. Because a minor role of fog during EN is expected, it is likely that measurements of fog-water collection during EN are considering drizzle and fog at the same time. Although recent studies on fog oases have shown some relationship with the ENSO, responses of vegetation have been largely based on descriptive data, the absence of large temporal records limit the establishment of a direct relationship with climatic oscillations. The second part of the results show that: 3) five different classes of different spectral values correspond to the main land cover of Lomas using a Vegetation Index (VI). The study case is characterised by shrubs and trees with variable cover (dense, semi-dense and open). A secondary area is covered by small shrubs where the dominant tree species is not present. The cacti area and the old terraces with open vegetation were not identified with the VI. Agriculture is present in the area. Finally, 4) contrary to the dry season of 2008 and 2009 years, a higher VI was obtained during the dry season of 2010. The VI increased up to three times their average value, showing a clear spectral signal change, which coincided with the ENSO event of that period.

Evaluation of Silver European EEL (anguilla anguilla) for the implementation of an Effective EEL Management plan in Mediterranean Coastal Lagoons

This thesis presents SEELF (Sustainable EEL fishery) Index, a methodology for evaluation of European eel (Anguilla anguilla) for the implementation of an effective Eel Management Plan, as defined by EU Regulation No.1100/2007. SEELF uses internal and external indices, age and blood parameters, and selects suitable specimen for restocking; it is also a reliable tool for eel stock management. In fact, SEELF Index, was developed in two versions: SEELF A, to be used in field operations (catch&release, eel status monitoring) and SEELF B to be used for quality control (food production) and research (eel status monitoring). Health status was evaluated also by biomarker analysis (ChE), and data were compared with age of eel. Age determination was performed with otolith reading and fish scale reading and a calibration between the two methods was possible. The study area was the Comacchio lagoon, a brackish coastal lagoon in Italy, well known as an example of suitable environment for eel fishery, where the capability to use the local natural resources has long been a key factor for a successful fishery management. Comacchio lagoon is proposed as an area where an effective EMP can be performed, in agreement with the main features (management of basins, reduction of mortality due to predators,etc.) highlighted for designation of European Restocking Area (ERA). The ERA is a new concept, proposed as a pillar of a new strategy on eel management and conservation. Furthermore, the features of ERAs can be useful in the framework of European Scale Eel Management Plan (ESEMP), proposed as a European scale implementation of EMP, providing a more effectiveness of conservation measures for eel management.

Dynamics of coastal aquifers: data-driven forecasting and risk analysis

This work is focused on the study of saltwater intrusion in coastal aquifers, and in particular on the realization of conceptual schemes to evaluate the risk associated with it. Saltwater intrusion depends on different natural and anthropic factors, both presenting a strong aleatory behaviour, that should be considered for an optimal management of the territory and water resources. Given the uncertainty of problem parameters, the risk associated with salinization needs to be cast in a probabilistic framework. On the basis of a widely adopted sharp interface formulation, key hydrogeological problem parameters are modeled as random variables, and global sensitivity analysis is used to determine their influence on the position of saltwater interface. The analyses presented in this work rely on an efficient model reduction technique, based on Polynomial Chaos Expansion, able to combine the best description of the model without great computational burden. When the assumptions of classical analytical models are not respected, and this occurs several times in the applications to real cases of study, as in the area analyzed in the present work, one can adopt data-driven techniques, based on the analysis of the data characterizing the system under study. It follows that a model can be defined on the basis of connections between the system state variables, with only a limited number of assumptions about the "physical" behaviour of the system.