Assessment of sea-level variability for the Emilia-Romagna coastal area in the framework of the Mediterranean Sea

Sea–level change is one of the ocean characteristics closely connected to climate change. Understanding its variation is essential since a large portion of the world’s population is located in low–lying locations. Two main techniques are employed to measure sea level: satellite altimetry and tide gauges. Satellite altimetry monitors sea–level relative to a geocentric reference, is unaffected by crustal processes and covers nearly the entire surface of the oceans since 1993. Conversely, tide gauges measure sea level at specific coastal locations and relative to a local ground benchmark, therefore are impacted by vertical land movements. In this study, the linear and non–linear geocentric and relative sea–level trends along the Emilia–Romagna coast (Northern Italy) have been analyzed over different periods. In order to assess the local sea–level variability, data from satellite altimetry and tide gauges have been compared over a common time interval (1993–2019), hence disentangling the contribute of vertical land movements. Non–linearity has been also evaluated at the broader scale of the Mediterranean Sea, in order to depict the main variability in geocentric sea–level trends from regional to sub–basin scale. Furthermore, the anthropogenic and natural influence at the river basin scale has been addressed, in order to shed light on the factors inducing the drastic reduction of riverine sediment supply to the Emilia–Romagna coast over the period 1920–2020. The findings of this analysis indicate that the sediment delivery reduction to the coast by rivers has been driven by several anthropogenic processes, acting on various spatiotemporal scales. Moreover, the local absolute sea–level trend is far from linear and appear "contaminated" by the presence of natural oscillations that act at the multi–decadal, quasi–decadal and inter–annual scale, mainly driven by both large–scale climatic modes and variations in local oceanography.

Sea-Level climate variability in the Mediterranean Sea

Sea-level variability is characterized by multiple interacting factors described in the Fourth Assessment Report (Bindoff et al., 2007) of the Intergovernmental Panel on Climate Change (IPCC) that act over wide spectra of temporal and spatial scales. In Church et al. (2010) sea-level variability and changes are defined as manifestations of climate variability and change. The European Environmental Agency (EEA) defines sea level as one of most important indicators for monitoring climate change, as it integrates the response of different components of the Earths system and is also affected by anthropogenic contributions (EEA, 2011). The balance between the different sea-level contributions represents an important source of uncertainty, involving stochastic processes that are very difficult to describe and understand in detail, to the point that they are defined as an enigma in Munk (2002). Sea-level rate estimates are affected by all these uncertainties, in particular if we look at possible responses to sea-level contributions to future climate. At the regional scale, lateral fluxes also contribute to sea-level variability, adding complexity to sea-level dynamics. The research strategy adopted in this work to approach such an interesting and challenging topic has been to develop an objective methodology to study sea-level variability at different temporal and spatial scales, applicable in each part of the Mediterranean basin in particular, and in the global ocean in general, using all the best calibrated sources of data (for the Mediterranean): in-situ, remote-sensig and numerical models data. The global objective of this work was to achieve a deep understanding of all of the components of the sea-level signal contributing to sea-level variability, tendency and trend and to quantify them.

The role of the Rocky Mountains in shaping the atmospheric mean state and the response to tropical forcing in idealised simulations.

Understanding the natural and forced variability of the atmospheric general circulation and its drivers is one of the grand challenges in climate science. It is of paramount importance to understand to what extent the systematic error of climate models affects the processes driving such variability. This is done by performing a set of simulations (ROCK experiments) with an intermediate complexity atmospheric model (SPEEDY), in which the Rocky Mountains orography is increased or decreased to influence the structure of the North Pacific jet stream. For each of these modified-orography experiments, the climatic response to idealized sea surface temperature anomalies of varying intensity in the El Niño Southern Oscillation (ENSO) region is studied. ROCK experiments are characterized by variations in the Pacific jet stream intensity whose extension encompasses the spread of the systematic error found in Coupled Model Intercomparison Project (CMIP6) models. When forced with ENSO-like idealised anomalies, they exhibit a non-negligible sensitivity in the response pattern over the Pacific North American region, indicating that the model mean state can affect the model response to ENSO. It is found that the classical Rossby wave train response to ENSO is more meridionally oriented when the Pacific jet stream is weaker and more zonally oriented with a stronger jet. Rossby wave linear theory suggests that a stronger jet implies a stronger waveguide, which traps Rossby waves at a lower latitude, favouring a zonal propagation of Rossby waves. The shape of the dynamical response to ENSO affects the ENSO impacts on surface temperature and precipitation over Central and North America. A comparison of the SPEEDY results with CMIP6 models suggests a wider applicability of the results to more resources-demanding climate general circulation models (GCMs), opening up to future works focusing on the relationship between Pacific jet misrepresentation and response to external forcing in fully-fledged GCMs.

The mesophotic zone of the Mediterranean Sea: spatial definition, biodiversity, and predictive models

The mesophotic zone is frequently defined as ranging between 30-40 and 150 m depth. However, these borders are necessarily imprecise due to variations in the penetration of light along the water column related to local factors. Moreover, density of data on mesophotic ecosystems vary along geographical distance, with temperate latitudes largely less explored than tropical situations. This is the case of the Mediterranean Sea, where information on mesophotic ecosystems is largely lower with respect to tropical situations. The lack of a clear definition of the borders of the mesophotic zone may represent a problem when information must be transferred to the policy that requires a coherent spatial definition to plan proper management and conservation measures. The present thesis aims at providing information on the spatial definition of the mesophotic zone in the Mediterranean Sea, its biodiversity and distribution of its ecosystems. The first chapter analyzes information on mesophotic ecosystems in the Mediterranean Sea to identify gaps in the literature and map the mesophotic zone in the Mediterranean Sea using light penetration estimated from satellite data. In the second chapter, different visual techniques to study mesophotic ecosystems are compared to identify the best analytical method to estimate diversity and habitat extension. In the third chapter, a set of Remotely Operated vehicles (ROV) surveys performed on mesophotic assemblages in the Mediterranean Sea are analyzed to describe their taxonomic and functional diversity and environmental factors influencing their structure. A Habitat Suitability Model is run in the fourth chapter to map the distribution of areas suitable for the presence of deep-water oyster reefs in the Adriatic-Ionian area. The fifth chapter explores the mesophotic zone in the northern Gulf of Mexico providing its spatial and vertical extension of the mesophotic zone and information on the diversity associated with mesophotic ecosystems.

Skeletal variability of the fibular ends and evolutionary implications

This thesis investigates the morphological variations of fibular extremities in humans and non-human hominids using a 3D Geometric Morphometric approach. The study has three objectives: (1) to assess the shape, form, and size variations of fibular epiphyses within the human species, highlighting sexually dimorphic features; (2) to explore interpopulation variability of fibular extremities from the Upper Paleolithic to the 20th century, comparing subsistence, mobility, and lifestyles; and (3) to examine interspecific variations in fibular ends, testing potential associations with locomotor and positional behavior among extant hominid taxa. In terms of intraspecific variations, sex-related differences in fibular form and size were observed, suggesting distinct functional requirements for the lower limb between sexes. Interpopulation variations revealed a decline in activity level over time, influenced by terrain and footwear use. Hunter-gatherer groups exhibited greater joint mobility, loading, and range of motion compared to sedentary pre- and post-industrial populations. Interspecific variations demonstrated significant morphological differences among hominid taxa, indicating functional implications related to both phylogeny and specific loading patterns on the lower limb. The study identified features indicative of bipedalism in humans, as well as shared characteristics among non-human great apes. Furthermore, distinguishing features were found between Asian and African apes, along with unique morphological signals associated with distinct positional behavior in each hominid taxa. By comprehensively analyzing fibular morphology, this research sheds light on the importance of this bone in knee support, ankle stabilization, and overall locomotor function. The findings contribute to our understanding of the evolutionary and functional aspects of the fibula across human populations and non-human hominids throughout history.

Satellite and in situ data integrated analysis to study the upper ocean an coastal environment of the Italian seas

The thesis objectives are to develop new methodologies for study of the space and time variability of Italian upper ocean ecosystem through the combined use of multi-sensors satellite data and in situ observations and to identify the capability and limits of remote sensing observations to monitor the marine state at short and long time scales. Three oceanographic basins have been selected and subjected to different types of analyses. The first region is the Tyrrhenian Sea where a comparative analysis of altimetry and lagrangian measurements was carried out to study the surface circulation. The results allowed to deepen the knowledge of the Tyrrhenian Sea surface dynamics and its variability and to defined the limitations of satellite altimetry measurements to detect small scale marine circulation features. Channel of Sicily study aimed to identify the spatial-temporal variability of phytoplankton biomass and to understand the impact of the upper ocean circulation on the marine ecosystem. An combined analysis of the satellite of long term time series of chlorophyll, Sea Surface Temperature and Sea Level field data was applied. The results allowed to identify the key role of the Atlantic water inflow in modulating the seasonal variability of the phytoplankton biomass in the region. Finally, Italian coastal marine system was studied with the objective to explore the potential capability of Ocean Color data in detecting chlorophyll trend in coastal areas. The most appropriated methodology to detect long term environmental changes was defined through intercomparison of chlorophyll trends detected by in situ and satellite. Then, Italian coastal areas subject to eutrophication problems were identified. This work has demonstrated that satellites data constitute an unique opportunity to define the features and forcing influencing the upper ocean ecosystems dynamics and can be used also to monitor environmental variables capable of influencing phytoplankton productivity.

The role of the West African Monsoon in the tropical to mid-latitudes climate

The interaction between atmosphere–land–ocean–biosphere systems plays a prominent role on the atmospheric dynamics and on the convective rainfall distribution over the West Africa monsoon area during the boreal summer. In particular, the initialization of convective systems in the Sub – Sahelian region has been directly linked to soil moisture heterogeneities identified as the major triggering, development and propagation of convective systems. The present study aims at investigating African monsoon large scale convective dynamics and rainfall diurnal cycle through an exploration of the hypothesis behind the mechanisms of a monsoon phenomenon as an emergence of a collective dynamics of many propagating convective systems. Such hypothesis is based on the existence of an internal self – regulation mechanism among the various components. To achieve these results a multiple analysis was performed based on remote sensed rainfall dataset, and global and regional modelling data for a period of 5 seasons: 2004 - 2008. Satellite rainfall data and convective occurrence variability were studied for assessing typical spatio – temporal signatures and characteristics with an emphasis to the diurnal cycle footprint. A global model and regional model simulation datasets, specifically developed for this analysis and based on Regional Atmospheric Modelling System – RAMS, have been analysed. Results from numerical model datasets highlight the evidence of a synchronization between the destabilization of the convective boundary layer and rainfall occurrence due to the solar radiation forcing through the latent heat release. This supports the conclusion that the studied interacting systems are associated with a process of mutual adjustment of rhythms. Furthermore, this rainfall internal coherence was studied in relation to the West African Heat Low pressure system, which has a prominent role in the large scale summer variability over the Mediterranean area since it is acting as one of dynamic link between sub tropical and midlatitudes variability.

Astronomical site testing in the era of the extremely large telescopes

The quality of astronomical sites is the first step to be considered to have the best performances from the telescopes. In particular, the efficiency of large telescopes in UV, IR, radio etc. is critically dependent on atmospheric transparency. It is well known that the random optical effects induced on the light propagation by turbulent atmosphere also limit telescope’s performances. Nowadays, clear appears the importance to correlate the main atmospheric physical parameters with the optical quality reachable by large aperture telescopes. The sky quality evaluation improved with the introduction of new techniques, new instrumentations and with the understanding of the link between the meteorological (or synoptical parameters and the observational conditions thanks to the application of the theories of electromagnetic waves propagation in turbulent medias: what we actually call astroclimatology. At the present the site campaigns are evolved and are performed using the classical scheme of optical seeing properties, meteorological parameters, sky transparency, sky darkness and cloudiness. New concept are added and are related to the geophysical properties such as seismicity, microseismicity, local variability of the climate, atmospheric conditions related to the ground optical turbulence and ground wind regimes, aerosol presence, use of satellite data. The purpose of this project is to provide reliable methods to analyze the atmospheric properties that affect ground-based optical astronomical observations and to correlate them with the main atmospheric parameters generating turbulence and affecting the photometric accuracy. The first part of the research concerns the analysis and interpretation of longand short-time scale meteorological data at two of the most important astronomical sites located in very different environments: the Paranal Observatory in the Atacama Desert (Chile), and the Observatorio del Roque de Los Muchachos(ORM) located in La Palma (Canary Islands, Spain). The optical properties of airborne dust at ORM have been investigated collecting outdoor data using a ground-based dust monitor. Because of its dryness, Paranal is a suitable observatory for near-IR observations, thus the extinction properties in the spectral range 1.00-2.30 um have been investigated using an empirical method. Furthermore, this PhD research has been developed using several turbulence profilers in the selection of the site for the European Extremely Large Telescope(E-ELT). During the campaigns the properties of the turbulence at different heights at Paranal and in the sites located in northern Chile and Argentina have been studied. This given the possibility to characterize the surface layer turbulence at Paranal and its connection with local meteorological conditions.

The genetic prehistory of Northern Calabria. Archeogenomic investigation of the human remains from Grotta di Pietra Sant’Angelo and Grotta della Monaca.

Given its strategic position at the center of the Mediterranean Basin, and its unique history of contacts and migrations, Calabria is an ideal region to decipher the genetic traces of at least some of the numerous demographic prehistoric events in Southern Italy. This thesis focuses on the genetic and social changes of ancient inhabitants of Calabria, covering a timeline of approximately 7000 to 3300 years ago, ranging from the Middle Neolithic to the Middle Bronze Age. We generated the first genome-wide data from Calabria, by focusing on the single inhumation of “Grotta di Pietra Sant’Angelo” (San Lorenzo Bellizzi, Cosenza) and on the vast community found buried in “Grotta della Monaca” (Sant’Agata di Esaro, Cosenza). Supported by archaeological evidence, the primary objective of this research was to employ paleogenomic evidence to decipher funerary customs, social organization, family ties, and demographic shifts in Southern Italy over a period extending beyond three millennia. The possibility of gender-related burial practices and kinship ties among the deceased was also explored. Subsequently, the biogeographical origin and ancestry of prehistoric people of Calabria was contextualized within the broad landscape of existing data on Mediterranean populations. By generating the first genomic evidence from prehistoric Calabria, unresolved questions were addressed, related to the appearance and persistence of distinct genomic components, such as the Iran-related and the Steppe-related ancestry, whose impact on ancient Southern Italian genomes remains uncharted.

Towards a broader understanding of the effects of tidal forcing on the global ocean circulation

The study of tides and their interactions with the complex dynamics of the global ocean represents a crucial challenge in ocean modelling. This thesis aims to deepen this study from a dynamical point of view, analysing what are the tidal effects on the general circulation of the ocean. We perform different experiments of a mesoscale-permitting global ocean model forced by both atmospheric fields and astronomical tidal potential, and we implement two parametrizations to include in the model tidal phenomena that are currently unresolved, with particular emphasis to the topographic wave drag for locally dissipating internal waves. An additional experiment using a mesoscale-resolving configuration is used to compare the simulated tides at different resolutions with observed data. We find that the accuracy of modelled tides strongly depends on the region and harmonic component of interest, even though the increased resolution allows to improve the modelled topography and resolve more intense internal waves. We then focus on the impact of tides in the Atlantic Ocean and find that tides weaken the overturning circulation during the analysed period from 1981 to 2007, even though the interannual differences strongly change in both amplitude and phase. The zonally integrated momentum balance shows that tide changes the water stratification at the zonal boundaries, modifying the pressure and therefore the geostrophic balance over the entire basin. Finally, we describe the overturning circulation in the Mediterranean Sea computing the meridional and zonal streamfunctions both in the Eulerian and residual frameworks. The circulation is characterised by different cells, and their forcing processes are described with particular emphasis to the role of mesoscale and a transient climatic event. We complete the description of the overturning circulation giving evidence for the first time to the connection between meridional and zonal cells.

Spatial and temporal variability and ecological processes in the epibenthic assemblages of the northern Adriatic Sea

Several coralligenous reefs occur in the soft bottoms of the northern Adriatic continental shelf. Mediterranean coralligenous habitats are characterised by high species diversity and are intrinsically valuable for their biological diversity and for the ecological processes they support. The conservation and management of these habitats require quantifying spatial and temporal variability of their benthic assemblages. This PhD thesis aims to give a relevant contribution to the knowledge of the structure and dynamics of the epibenthic assemblages on the coralligenous subtidal reefs occurring in the northern Adriatic Sea. The epibenthic assemblages showed a spatial variation larger compared to temporal changes, with a temporal persistence of reef-forming organisms. Assemblages spatial heterogeneity has been related to morphological features and geographical location of the reefs, together with variation in the hydrological conditions. Manipulative experiments help to understand the ecological processes structuring the benthic assemblages and maintaining their diversity. In this regards a short and long term experiment on colonization patterns of artificial substrata over a 3-year period has been performed in three reefs, corresponding to the three main types of assemblages detected in the previous study. The first colonisers, largely depending by the different larval supply, played a key role in determining the heterogeneity of the assemblages in the early stage of colonisation. Lateral invasion, from the surrounding assemblages, was the driver in structuring the mature assemblages. These complex colonisation dynamics explained the high heterogeneity of the assemblages dwelling on the northern Adriatic biogenic reefs. The buildup of these coralligenous reefs mainly depends by the bioconstruction-erosion processes that has been analysed through a field experiment. Bioconstruction, largely due to serpulid polychaetes, prevailed on erosion processes and occurred at similar rates in all sites. Similarly, the total energy contents in the benthic communities do not differ among sites, despite being provided by different species. Therefore, we can hypothesise that both bioconstruction processes and energetic storage may be limited by the availability of resources. Finally the major contribution of the zoobenthos compared to the phytobenthos to the total energetic content of assemblages suggests that the energy flow in these benthic habitats is primarily supported by planktonic food web trough the filter feeding invertebrates.

The stratigraphic record of the quaternary sea level fluctuations and the impact of the post-glacial sea level rise (Termination I) in the Adriatic basin (Mediterranean sea)

The modern stratigraphy of clastic continental margins is the result of the interaction between several geological processes acting on different time scales, among which sea level oscillations, sediment supply fluctuations and local tectonics are the main mechanisms. During the past three years my PhD was focused on understanding the impact of each of these process in the deposition of the central and northern Adriatic sedimentary successions, with the aim of reconstructing and quantifying the Late Quaternary eustatic fluctuations. In the last few decades, several Authors tried to quantify past eustatic fluctuations through the analysis of direct sea level indicators, among which drowned barrier-island deposits or coral reefs, or indirect methods, such as Oxygen isotope ratios (δ18O) or modeling simulations. Sea level curves, obtained from direct sea level indicators, record a composite signal, formed by the contribution of the global eustatic change and regional factors, as tectonic processes or glacial-isostatic rebound effects: the eustatic signal has to be obtained by removing the contribution of these other mechanisms. To obtain the most realistic sea level reconstructions it is important to quantify the tectonic regime of the central Adriatic margin. This result has been achieved integrating a numerical approach with the analysis of high-resolution seismic profiles. In detail, the subsidence trend obtained from the geohistory analysis and the backstripping of the borehole PRAD1.2 (the borehole PRAD1.2 is a 71 m continuous borehole drilled in -185 m of water depth, south of the Mid Adriatic Deep - MAD - during the European Project PROMESS 1, Profile Across Mediterranean Sedimentary Systems, Part 1), has been confirmed by the analysis of lowstand paleoshorelines and by benthic foraminifera associations investigated through the borehole. This work showed an evolution from inner-shelf environment, during Marine Isotopic Stage (MIS) 10, to upper-slope conditions, during MIS 2. Once the tectonic regime of the central Adriatic margin has been constrained, it is possible to investigate the impact of sea level and sediment supply fluctuations on the deposition of the Late Pleistocene-Holocene transgressive deposits. The Adriatic transgressive record (TST - Transgressive Systems Tract) is formed by three correlative sedimentary bodies, deposited in less then 14 kyr since the Last Glacial Maximum (LGM); in particular: along the central Adriatic shelf and in the adjacent slope basin the TST is formed by marine units, while along the northern Adriatic shelf the TST is represented by costal deposits in a backstepping configuration. The central Adriatic margin, characterized by a thick transgressive sedimentary succession, is the ideal site to investigate the impact of late Pleistocene climatic and eustatic fluctuations, among which Meltwater Pulses 1A and 1B and the Younger Dryas cold event. The central Adriatic TST is formed by a tripartite deposit bounded by two regional unconformities. In particular, the middle TST unit includes two prograding wedges, deposited in the interval between the two Meltwater Pulse events, as highlighted by several 14C age estimates, and likely recorded the Younger Dryas cold interval. Modeling simulations, obtained with the two coupled models HydroTrend 3.0 and 2D-Sedflux 1.0C (developed by the Community Surface Dynamics Modeling System - CSDMS), integrated by the analysis of high resolution seismic profiles and core samples, indicate that: 1 - the prograding middle TST unit, deposited during the Younger Dryas, was formed as a consequence of an increase in sediment flux, likely connected to a decline in vegetation cover in the catchment area due to the establishment of sub glacial arid conditions; 2 - the two-stage prograding geometry was the consequence of a sea level still-stand (or possibly a fall) during the Younger Dryas event. The northern Adriatic margin, characterized by a broad and gentle shelf (350 km wide with a low angle plunge of 0.02° to the SE), is the ideal site to quantify the timing of each steps of the post LGM sea level rise. The modern shelf is characterized by sandy deposits of barrier-island systems in a backstepping configuration, showing younger ages at progressively shallower depths, which recorded the step-wise nature of the last sea level rise. The age-depth model, obtained by dated samples of basal peat layers, is in good agreement with previous published sea level curves, and highlights the post-glacial eustatic trend. The interval corresponding to the Younger Dyas cold reversal, instead, is more complex: two coeval coastal deposits characterize the northern Adriatic shelf at very different water depths. Several explanations and different models can be attempted to explain this conundrum, but the problem remains still unsolved.

Origin and variability of PM10 and atmospheric radiotracers at the WMO-GAW station of Mt. Cimone (1998-2011) and in the central Po Valley

Particulate matter is one of the main atmospheric pollutants, with a great chemical-environmental relevance. Improving knowledge of the sources of particulate matter and of their apportionment is needed to handle and fulfill the legislation regarding this pollutant, to support further development of air policy as well as air pollution management. Various instruments have been used to understand the sources of particulate matter and atmospheric radiotracers at the site of Mt. Cimone (44.18° N, 10.7° E, 2165 m asl), hosting a global WMO-GAW station. Thanks to its characteristics, this location is suitable investigate the regional and long-range transport of polluted air masses on the background Southern-Europe free-troposphere. In particular, PM10 data sampled at the station in the period 1998-2011 were analyzed in the framework of the main meteorological and territorial features. A receptor model based on back trajectories was applied to study the source regions of particulate matter. Simultaneous measurements of atmospheric radionuclides Pb-210 and Be-7 acquired together with PM10 have also been analysed to acquire a better understanding of vertical and horizontal transports able to affect atmospheric composition. Seasonal variations of atmospheric radiotracers have been studied both analysing the long-term time series acquired at the measurement site as well as by means of a state-of-the-art global 3-D chemistry and transport model. Advection patterns characterizing the circulation at the site have been identified by means of clusters of back-trajectories. Finally, the results of a source apportionment study of particulate matter carried on in a midsize town of the Po Valley (actually recognised as one of the most polluted European regions) are reported. An approach exploiting different techniques, and in particular different kinds of models, successfully achieved a characterization of the processes/sources of particulate matter at the two sites, and of atmospheric radiotracers at the site of Mt. Cimone.