Mid-late Holocene environmental change and human activities in the northern Apennines, Italy

Radiocarbon-dated palaeoecological records from the upland zone of the northern Apennines spanning the Mid-Late Holocene (last 7000 years) have been evaluated using established criteria for detecting anthropogenic impact on the landscape and environment. The integrated palaeoecological records across the study area collectively indicate human interference with natural vegetation succession and landscape modification from at least the Middle Neolithic. These activities resulted in the progressive decline of Abies, Ulmus, Fraxinus and Tilia, and the spread of Fagus, from ∼7000 cal BP, accompanied at various times by evidence for biomass burning, soil erosion, the expansion of shrubland and herbaceous taxa, and the possible cultivation of Olea, Juglans and Castanea. Comparison of these data with the archaeological scheme for the region, and the climate history of the central-western Mediterranean, has revealed that the palaeoecological records broadly support the archaeological evidence, but suggest that several key vegetation changes also coincide with important periods of climate change, especially at ∼7800–5000 cal BP.

Late Würm and Early-Middle Holocene environmental change and human activities in the Northern Apennines, Italy

An enhanced radiocarbon-dated pollen-stratigraphical record from Rovegno (Liguria, 812m asl), northern Apennines (Italy), has provided a history of vegetation succession from before 17,056-16,621 cal yrs BP to the present day. The record indicates the transition from open Pinus woodland to Artemisia dominated grassland, and finally Juniperus shrubland during the late Würm. This is succeeded by Betula and Pinus woodland, and the expansion of thermophilous taxa, namely Abies, Corylus and Quercus during the Late Würm Lateglacial Interstadial. The ‘Younger Dryas’ is possibly represented by an increase in Betula and Artemisia. During the early Holocene, mixed coniferous-deciduous woodland is dominant with Quercus, as well as Abies, Fagus and Corylus. Fagus woodland becomes established sometime before 6488-6318 cal yrs BP, but never becomes a major component of the woodland cover. Throughout the middle Holocene, Abies woodland fl uctuates, with marked declines between 6488-6318 cal yrs BP and 5287-4835 cal yrs BP, although the cause remains uncertain. Finally, the paper evaluates the application of non-pollen palynomorphs, especially coprophilous fungal spores, at Prato Spilla ‘A’ (Emilia Romagna) and concludes that greater caution must be used when interpreting middle Holocene human activity based upon pollen data alone

Shallow and deep deformation in northern Apennines region using seismological data

For its particular position and the complex geological history, the Northern Apennines has been considered as a natural laboratory to apply several kinds of investigations. By the way, it is complicated to joint all the knowledge about the Northern Apennines in a unique picture that explains the structural and geological emplacement that produced it. The main goal of this thesis is to put together all information on the deformation - in the crust and at depth - of this region and to describe a geodynamical model that takes account of it. To do so, we have analyzed the pattern of deformation in the crust and in the mantle. In both cases the deformation has been studied using always information recovered from earthquakes, although using different techniques. In particular the shallower deformation has been studied using seismic moment tensors information. For our purpose we used the methods described in Arvidsson and Ekstrom (1998) that allowing the use in the inversion of surface waves [and not only of the body waves as the Centroid Moment Tensor (Dziewonski et al., 1981) one] allow to determine seismic source parameters for earthquakes with magnitude as small as 4.0. We applied this tool in the Northern Apennines and through this activity we have built up the Italian CMT dataset (Pondrelli et al., 2006) and the pattern of seismic deformation using the Kostrov (1974) method on a regular grid of 0.25 degree cells. We obtained a map of lateral variations of the pattern of seismic deformation on different layers of depth, taking into account the fact that shallow earthquakes (within 15 km of depth) in the region occur everywhere while most of events with a deeper hypocenter (15-40 km) occur only in the outer part of the belt, on the Adriatic side. For the analysis of the deep deformation, i.e. that occurred in the mantle, we used the anisotropy information characterizing the structure below the Northern Apennines. The anisotropy is an earth properties that in the crust is due to the presence of aligned fluid filled cracks or alternating isotropic layers with different elastic properties while in the mantle the most important cause of seismic anisotropy is the lattice preferred orientation (LPO) of the mantle minerals as the olivine. This last is a highly anisotropic mineral and tends to align its fast crystallographic axes (a-axis) parallel to the astenospheric flow as a response to finite strain induced by geodynamic processes. The seismic anisotropy pattern of a region is measured utilizing the shear wave splitting phenomenon (that is the seismological analogue to optical birefringence). Here, to do so, we apply on teleseismic earthquakes recorded on stations located in the study region, the Sileny and Plomerova (1996) approach. The results are analyzed on the basis of their lateral and vertical variations to better define the earth structure beneath Northern Apennines. We find different anisotropic domains, a Tuscany and an Adria one, with a pattern of seismic anisotropy which laterally varies in a similar way respect to the seismic deformation. Moreover, beneath the Adriatic region the distribution of the splitting parameters is so complex to request an appropriate analysis. Therefore we applied on our data the code of Menke and Levin (2003) which allows to look for different models of structures with multilayer anisotropy. We obtained that the structure beneath the Po Plain is probably even more complicated than expected. On the basis of the results obtained for this thesis, added with those from previous works, we suggest that slab roll-back, which created the Apennines and opened the Tyrrhenian Sea, evolved in the north boundary of Northern Apennines in a different way from its southern part. In particular, the trench retreat developed primarily south of our study region, with an eastward roll-back. In the northern portion of the orogen, after a first stage during which the retreat was perpendicular to the trench, it became oblique with respect to the structure.

Tectonic geomorphology and active strain of the Northern Apennines mountain front

The Northern Apennines (NA) chain is the expression of the active plate margin between Europe and Adria. Given the low convergence rates and the moderate seismic activity, ambiguities still occur in defining a seismotectonic framework and many different scenarios have been proposed for the mountain front evolution. Differently from older models that indicate the mountain front as an active thrust at the surface, a recently proposed scenario describes the latter as the frontal limb of a long-wavelength fold (> 150 km) formed by a thrust fault tipped around 17 km at depth, and considered as the active subduction boundary. East of Bologna, this frontal limb is remarkably very straight and its surface is riddled with small, but pervasive high- angle normal faults. However, west of Bologna, some recesses are visible along strike of the mountain front: these perturbations seem due to the presence of shorter wavelength (15 to 25 km along strike) structures showing both NE and NW-vergence. The Pleistocene activity of these structures was already suggested, but not quantitative reconstructions are available in literature. This research investigates the tectonic geomorphology of the NA mountain front with the specific aim to quantify active deformations and infer possible deep causes of both short- and long-wavelength structures. This study documents the presence of a network of active extensional faults, in the foothills south and east of Bologna. For these structures, the strain rate has been measured to find a constant throw-to-length relationship and the slip rates have been compared with measured rates of erosion. Fluvial geomorphology and quantitative analysis of the topography document in detail the active tectonics of two growing domal structures (Castelvetro - Vignola foothills and the Ghiardo plateau) embedded in the mountain front west of Bologna. Here, tilting and river incision rates (interpreted as that long-term uplift rates) have been measured respectively at the mountain front and in the Enza and Panaro valleys, using a well defined stratigraphy of Pleistocene to Holocene river terraces and alluvial fan deposits as growth strata, and seismic reflection profiles relationships. The geometry and uplift rates of the anticlines constrain a simple trishear fault propagation folding model that inverts for blind thrust ramp depth, dip, and slip. Topographic swath profiles and the steepness index of river longitudinal profiles that traverse the anti- clines are consistent with stratigraphy, structures, aquifer geometry, and seismic reflection profiles. Available focal mechanisms of earthquakes with magnitude between Mw 4.1 to 5.4, obtained from a dataset of the instrumental seismicity for the last 30 years, evidence a clear vertical separation at around 15 km between shallow extensional and deeper compressional hypocenters along the mountain front and adjacent foothills. In summary, the studied anticlines appear to grow at rates slower than the growing rate of the longer- wavelength structure that defines the mountain front of the NA. The domal structures show evidences of NW-verging deformation and reactivations of older (late Neogene) thrusts. The reconstructed river incision rates together with rates coming from several other rivers along a 250 km wide stretch of the NA mountain front and recently available in the literature, all indicate a general increase from Middle to Late Pleistocene. This suggests focusing of deformation along a deep structure, as confirmed by the deep compressional seismicity. The maximum rate is however not constant along the mountain front, but varies from 0.2 mm/yr in the west to more than 2.2 mm/yr in the eastern sector, suggesting a similar (eastward-increasing) trend of the apenninic subduction.

Characterization of fold-related fractures in the carbonate rocks of the Cingoli anticline, northern Apennines, Italy

Thrust fault-related folds in carbonate rocks are characterized by deformation accommodated by different structures, such as joints, faults, pressure solution seams, and deformation bands. Defining the development of fracture systems related to the folding process is significant both for theoretical and practical purposes. Fracture systems are useful constrains in order to understand the kinematical evolution of the fold. Furthermore, understanding the relationships between folding and fracturing provides a noteworthy contribution for reconstructing the geodynamic and the structural evolution of the studied area. Moreover, as fold-related fractures influence fluid flow through rocks, fracture systems are relevant for energy production (geothermal studies, methane and CO2 , storage and hydrocarbon exploration), environmental and social issues (pollutant distribution, aquifer characterization). The PhD project shows results of a study carried out in a multilayer carbonate anticline characterized by different mechanical properties. The aim of this study is to understand the factors which influence the fracture formation and to define their temporal sequence during the folding process. The studied are is located in the Cingoli anticline (Northern Apennines), which is characterized by a pelagic multilayer characterized by sequences with different mechanical stratigraphies. A multi-scale analysis has been made in several outcrops located in different structural positions. This project shows that the conceptual sketches proposed in literature and the strain distribution models outline well the geometrical orientation of most of the set of fractures observed in the Cingoli anticline. On the other hand, the present work suggests the relevance of the mechanical stratigraphy in particular controlling the type of fractures formed (e.g. pressure solution seams, joints or shear fractures) and their subsequent evolution. Through a multi-scale analysis, and on the basis of the temporal relationship between fracture sets and their orientation respect layering, I also suggest a conceptual model for fracture systems formation.

Quality assessment of a landslide inventory map and its application to land‐use planning. A case study in the Northern Apennines (Emilia‐Romagna region, Italy)

Landslide hazard and risk are growing as a consequence of climate change and demographic pressure. Land‐use planning represents a powerful tool to manage this socio‐economic problem and build sustainable and landslide resilient communities. Landslide inventory maps are a cornerstone of land‐use planning and, consequently, their quality assessment represents a burning issue. This work aimed to define the quality parameters of a landslide inventory and assess its spatial and temporal accuracy with regard to its possible applications to land‐use planning. In this sense, I proceeded according to a two‐steps approach. An overall assessment of the accuracy of data geographic positioning was performed on four case study sites located in the Italian Northern Apennines. The quantification of the overall spatial and temporal accuracy, instead, focused on the Dorgola Valley (Province of Reggio Emilia). The assessment of spatial accuracy involved a comparison between remotely sensed and field survey data, as well as an innovative fuzzylike analysis of a multi‐temporal landslide inventory map. Conversely, long‐ and short‐term landslide temporal persistence was appraised over a period of 60 years with the aid of 18 remotely sensed image sets. These results were eventually compared with the current Territorial Plan for Provincial Coordination (PTCP) of the Province of Reggio Emilia. The outcome of this work suggested that geomorphologically detected and mapped landslides are a significant approximation of a more complex reality. In order to convey to the end‐users this intrinsic uncertainty, a new form of cartographic representation is needed. In this sense, a fuzzy raster landslide map may be an option. With regard to land‐use planning, landslide inventory maps, if appropriately updated, confirmed to be essential decision‐support tools. This research, however, proved that their spatial and temporal uncertainty discourages any direct use as zoning maps, especially when zoning itself is associated to statutory or advisory regulations.

The Early Miocene “Bisciaro volcaniclastic event” (northern Apennines, Italy): a key study for the geodynamic evolution of the central-western Mediterranean

The Early Miocene Bisciaro Fm., a marly limestone succession cropping out widely in the Umbria–Romagna–Marche Apennines, is characterized by a high amount of volcaniclastic content, characterizing this unit as a peculiar event of the Adria Plate margin. Because of this volcaniclastic event, also recognizable in different sectors of the central-western Mediterranean chains, this formation is proposed as a “marker” for the geodynamic evolution of the area. In the Bisciaro Fm., the volcaniclastic supply starts with the “Raffaello” bed (Earliest Aquitanian) that marks the base of the formation and ends in the lower portion of the Schlier Fm. (Late Burdigalian–Langhian p.p.). Forty-one studied successions allowed the recognition of three main petrofacies: (1) Pyroclastic Deposits (volcanic materials more than 90 %) including the sub-petrofacies 1A, Vitroclastic/crystallo-vitroclastic tuffs; 1B, Bentonitic deposits; and 1C, Ocraceous and blackish layers; (2) Resedimented Syn-Eruptive Volcanogenic Deposits (volcanic material 30–90 %) including the sub-petrofacies 2A, High-density volcanogenic turbidites; 2B, Low-density volcanogenic turbidites; 2C, Crystal-rich volcanogenic deposits; and 2D, Glauconitic-rich volcaniclastites; (3) Mixing of Volcaniclastic Sediments with Marine Deposits (volcanic material 5–30 %, mixed with marine sediments: marls, calcareous marls, and marly limestones). Coeval volcaniclastic deposits recognizable in different tectonic units of the Apennines, Maghrebian, and Betic Chains show petrofacies and chemical–geochemical features related to a similar calc-alkaline magmatism. The characterization of this event led to the hypothesis of a co-genetic relationship between volcanic activity centres (primary volcanic systems) and depositional basins (depositional processes) in the Early Miocene palaeogeographic and palaeotectonic evolution of the central-western Mediterranean region. The results support the proposal of a geodynamic model of this area that considers previously proposed interpretations.

The tectonically-confined Firenzuola turbidite system (Miocene, Marnoso-arenacea Formation, northern Apennines, Italy): an example of facies variations related to the basin morphology

The Firenzuola turbidite system formed during a paroxysmal phase of thrust propagation, involving the upper Serravallian deposits of the Marnoso-arenacea Formation (MAF). During this phase the coeval growth of two major tectonic structures, the M. Castellaccio thrust and the Verghereto high, played a key role, causing a closure of the inner basin and a coeval shift of the depocentre to the outer basin. This work focuses on this phase of fragmentation of the MAF basin; it is based on a new detailed high-resolution stratigraphic framework, which was used to determine the timing of growth of the involved structures and their direct influence on sediment dispersal and on the lateral and vertical turbidite facies distribution. The Firenzuola turbidite system stratigraphy is characterized by the occurrence of mass-transport complexes (MTCs) and thick sandstone accumulation in the depocentral area, which passes to finer drape over the structural highs; the differentiation between these two zones increases over time and ends with the deposition of marly units over the structural highs and the emplacement of the Visignano MTC. According to the stratigraphic pattern and turbidite facies characteristics, the Firenzuola System has been split into two phases, namely Firenzuola I and Firenzuola II: the former is quite similar to the underlying deposits, while the latter shows the main fragmentation phase, testifying the progressive isolation of the inner basin and a coeval shift of the depocentre to the outer basin. The final stratigraphic and sedimentological dataset has been used to create a quantitative high-resolution 3D facies distribution using the Petrel software platform. This model allows a detailed analysis of lateral and vertical facies variations that can be exported to several reservoirs settings in hydrocarbon exploration and exploitation areas, since facies distributions and geometries of the reservoir bodies of many sub-surface turbidite basins show a significant relationship to the syndepositional structural activity, but are beyond seismic resolution.

Tectonics of the wedge-top Epiligurian basins and implications for the syn-to-post orogenic evolution of the Northern Apennines fold-and-thrust belt (Italy)

This thesis has the aim to give an overview about the tectonic history of the Epiligurian units, which crop out in the axial portion of the Northern Apennines fold-and-thrust belt, from a structural and thermal point of view, through a multiscalar and multitecnique approach. I focused on a key example of Epiligurian wedge-top basin, (Marzabotto Basin) proceeding from macro-to-microscale approach. The study started from a remote sensing analysis of the lineaments and morphostructures which affected the study area to obtain the regional faulting pattern and an overview about the main tectonic structures, used as basis for the structural investigation at the mesoscale. On the basis of this, it was possible to reconstruct the succession of tectonic events that affected the Marzabotto Basin, consisting in: i) two sets of thrusts indicating a NE-SW and NW-SE shortening of the sedimentary succession; ii) NE-SW-left lateral transtensional faults related to a strike-slip tectonic phase; iii) three main sets of extensional structures which cut and displace the previous thrusts. Normal faults are related to the post-orogenic evolution and have been dated with U-Th method, getting an age of Middle-Late Pleistocene. From a thermal point of view, apatite fission-tracks and (U-Th)/He analyses of detrital minerals and thermal modelling on the middle-upper Eocene siliciclastic deposits allowed me to better constrain the local exhumation history and correlate it with the large-scale tectonic evolution of the Northern Apennines. In particular, the Marzabotto Basin experienced a complex burial-exhumation history, consisting in two cooling cooling phases related to the growth of the Northern Apennines belt (Oligo-Miocene in age) and a later cooling which tracks the accretion in the orogenic wedge concomitant with rollback-driven extension (late Miocene-Pliocene in age). In conclusion it is possible to affirm that the study shed new light on poorly constrained elements of fold-and-thrust belt.

Neogene and Quaternary rollback evolution of the Tyrrhenian Sea, the Apennines, and the Sicilian Maghrebides

Reconstruction of the evolution of the Tyrrhenian Sea shows that the major stage of rifting associated with the opening of this basin began at similar to10 Ma. It involved two episodes of back arc extension, which were induced by the rollback of a west dipping subducting slab. The first period of extension (10-6 Ma) was prominent in the northern Tyrrhenian Sea and in the western part of the southern Tyrrhenian Sea. The second period of extension, mainly affected the southern Tyrrhenian Sea, began in the latest Messinian (6-5 Ma) and has been accompanied by subduction rollback at rates of 60-100 km Myr(-1). Slab reconstruction, combined with paleomagnetic and paleogeographic constraints, indicates that in the central Apennines, the latest Messinian (6-5 Ma) arrival of a carbonate platform at the subduction zone impeded subduction and initiated a slab tear and major strike-slip faults. These processes resulted in the formation of a narrow subducting slab beneath the Ionian Sea that has undergone faster subduction rollback and induced extreme rates of back arc extension.

Early–Middle Holocene vegetation history, climate change and human activities at Lago Riane (Ligurian Apennines, NW Italy)

The radiocarbon-dated palaeoecological study of Lago Riane (Ligurian Apennines, NW Italy) presented here forms part of a wider investigation into the relationships between Holocene vegetation succession, climate change and human activities in the northern Apennines. The record of vegetation history from Lago Riane indicates that, since the end of the last glaciation, climate change and prehistoric human activities, combined with several local factors, have strongly influenced the pattern and timing of natural vegetation succession. The pollen record indicates an important change in vegetation cover at Lago Riane at ~8500–8200 cal. years b.p., coincident with a well-known period of rapid climate change. At ~6100 cal. years b.p., Fagus woodland colonised Lago Riane during a period of climate change and expansion of Late Neolithic human activities in the upland zone of Liguria. A marked decline in Abies woodland, and the expansion of Fagus woodland, at ~4700 cal. years b.p., coincided with further archaeological evidence for pastoralism in the mountains of Liguria during the Copper Age. At ~3900–3600 cal. years b.p. (Early to Middle Bronze Age transition), a temporary expansion of woodland at Lago Riane has been provisionally attributed to a decline in human pressure on the environment during a period of short-term climate change

A 300 million year-long history. The metamorphic evolution of the Northern Apennine Variscan basement

The Variscan basement of Northern Apennines (Northern Italy) is a polymetamorphic portion of continental crust. This thesis investigated the metamorphic history of this basement occurring in the Cerreto Pass, in the Pontremoli well, and in the Pisani Mountains. The study comprised fieldwork, petrography and microstructural analysis, determination of the bulk rock and mineral composition, thermodynamic modelling, conventional geothermobarometry, monazite chemical dating and Ar/Ar dating of muscovite. The reconstructed metamorphic evolution of the selected samples allowed to define a long-lasting metamorphic history straddling the Variscan and Alpine orogenesis. Some general petrological issues generally found in low- to medium-grade metapelites were also tackled: (i) With middle-grade micaschist it is possible to reconstruct a complete P-T-D path by combining microstructural analysis and thermodynamic modelling. Prekinematic white mica may preserve Mg-rich cores related to the pre-peak stage. Mn-poor garnet rim records the peak metamorphism. Na-rich mylonitic white mica, the XFe of chlorite and the late paragenesis may constrain the retrograde stage. (ii) Metapelites may contain coronitic microstructures of apatite + Th-silicate, allanite and epidote around unstable monazite grains. Chemistry and microstructure of Th-rich monazite relics surrounded by this coronitic microstructure may suggest that monazite mineral was inherited and underwent partial dissolution and fluid-aided replacement by REE-accessory minerals at 500-600°C and 5-7 kbar. (iii) Fish-shaped white mica is not always a (prekinematic) mica-fish. Observed at high-magnification BSE images it may consist of several white mica formed during a mylonitic stage. Hence, the asymmetric foliation boudin is a suitable microstructure to obtain geochronological information about the shearing stage. (iv) Thermodynamic modelling of a hematite-rich metasedimentary rock fails to reproduce the observed mineral compositions when the bulk Fe2O3 is neglected or determined through titration. The mismatch between observed and computed mineral compositions and assemblage is resolved by tuning the effective ferric iron content by P-XFe2O3 diagrams.

Analysis and compound of carbon sequestration and serpentinization from the Iberian Margin and the Northern Apennine

Fluid circulation in peridotite-hosted hydrothermal systems influences the incorporation of carbon into the oceanic crust and its long-term storage. At low to moderate temperatures, serpentinization of peridotite produces alkaline fluids that are rich in CH4 and H2. Upon mixing with seawater, these fluids precipitate carbonate, forming an extensive network of calcite veins in the basement rocks, while H2 and CH4 serve as an energy source for microorganisms. Here, we analyzed the carbon geochemistry of two ancient peridotite-hosted hydrothermal systems: 1) ophiolites cropping out in the Northern Apennines, and 2) calcite-veined serpentinites from the Iberian Margin (Ocean Drilling Program (ODP) Legs 149 and 173), and compare them to active peridotite-hosted hydrothermal systems such as the Lost City hydrothermal field (LCHF) on the Atlantis Massif near the Mid-Atlantic Ridge (MAR). Our results show that large amounts of carbonate are formed during serpentinization of mantle rocks exposed on the seafloor (up to 9.6 wt.% C in ophicalcites) and that carbon incorporation decreases with depth. In the Northern Apennine serpentinites, serpentinization temperatures decrease from 240 °C to < 150 °C, while carbonates are formed at temperatures decreasing from ~ 150 °C to < 50 °C. At the Iberian Margin both carbonate formation and serpentinization temperatures are lower than in the Northern Apennines with serpentinization starting at ~ 150 °C, followed by clay alteration at < 100 °C and carbonate formation at < 19-44 °C. Comparison with various active peridotite-hosted hydrothermal systems on the MAR shows that the serpentinites from the Northern Apennines record a thermal evolution similar to that of the basement of the LCHF and that tectonic activity on the Jurassic seafloor, comparable to the present-day processes leading to oceanic core complexes, probably led to formation of fractures and faults, which promoted fluid circulation to greater depth and cooling of the mantle rocks. Thus, our study provides further evidence that the Northern Apennine serpentinites host a paleo-stockwork of a hydrothermal system similar to the basement of the LCHF. Furthermore, we argue that the extent of carbonate uptake is mainly controlled by the presence of fluid pathways. Low serpentinization temperatures promote microbial activity, which leads to enhanced biomass formation and the storage of organic carbon. Organic carbon becomes dominant with increasing depth and is the principal carbon phase at more than 50-100 m depth of the serpentinite basement at the Iberian Margin. We estimate that annually 1.1 to 2.7 × 1012 g C is stored within peridotites exposed to seawater, of which 30-40% is fixed within the uppermost 20-50 m mainly as carbonate. Additionally, we conclude that alteration of oceanic lithosphere is an important factor in the long-term global carbon cycle, having the potential to store carbon for millions of years.