The Central-Western Mediterranean: Anomalous igneous activity in an anomalous collisional tectonic setting

The central-western Mediterranean area is a key region for understanding the complex interaction between igneous activity and tectonics. In this review, the specific geochemical character of several ‘subduction-related’ Cenozoic igneous provinces are described with a view to identifying the processes responsible for the modifications of their sources. Different petrogenetic models are reviewed in the light of competing geological and geodynamic scenarios proposed in the literature.Plutonic rocks occur almost exclusively in the Eocene–Oligocene Periadriatic Province of the Alps while relatively minor plutonic bodies (mostly Miocene in age) crop out in N Morocco, S Spain and N Algeria. Igneous activity is otherwise confined to lava flows and dykes accompanied by relatively greater volumes of pyroclastic (often ignimbritic) products. Overall, the igneous activity spanned a wide temporal range, from middle Eocene (such as the Periadriatic Province) to the present (as in the Neapolitan of southern Italy). The magmatic products are mostly SiO₂-oversaturated, showing calcalkaline to high-K calcalcaline affinity, except in some areas (as in peninsular Italy) where potassic to ultrapotassic compositions prevail. The ultrapotassic magmas (which include leucitites to leucite-phonolites) are dominantly SiO₂-undersaturated, although rare, SiO₂-saturated (i.e., leucite-free lamproites) appear over much of this region, examples being in the Betics (southeast Spain), the northwest Alps, northeast Corsica (France), Tuscany (northwest Italy), southeast Tyrrhenian Sea (Cornacya Seamount) and possibly in the Tell region (northeast Algeria).Excepted for the Alpine case, subduction-related igneous activity is strictly linked to the formation of the Mediterranean Sea. This Sea, at least in its central and western sectors, is made up of several young (< 30 Ma) V-shaped back-arc basins plus several dispersed continental fragments, originally in crustal continuity with the European plate (Sardinia, Corsica, Balearic Islands, Kabylies, Calabria, Peloritani Mountains). The bulk of igneous activity in the central-western Mediterranean is believed to have tapped mantle ‘wedge’ regions, metasomatized by pressure-related dehydration of the subducting slabs. The presence of subduction-related igneous rocks with a wide range of chemical composition has been related to the interplay of several factors among which the pre-metasomatic composition of the mantle wedges (i.e., fertile vs. refractory mineralogy), the composition of the subducting plate (i.e., the type and amount of sediment cover and the alteration state of the crust), the variable thermo-baric conditions of magma formation, coupled with variable molar concentrations of CO₂ and H₂O in the fluid phase released by the subducting plates are the most important.Compared to classic collisional settings (e.g., Himalayas), the central-western Mediterranean area shows a range of unusual geological and magmatological features. These include: a) the rapid formation of extensional basins in an overall compressional setting related to Africa-Europe convergence; b) centrifugal wave of both compressive and extensional tectonics starting from a ‘pivotal’ region around the Gulf of Lyon; c) the development of concomitant Cenozoic subduction zones with different subduction and tectonic transport directions; d) subduction ‘inversion’ events (e.g., currently along the Maghrebian coast and in northern Sicily, previously at the southern paleo-European margin); e) a repeated temporal pattern whereby subduction-related magmatic activity gives way to magmas of intraplate geochemical type; f) the late-stage appearance of magmas with collision-related ‘exotic’ (potassic to ultrapotassic) compositions, generally absent from simple subduction settings; g) the relative scarcity of typical calcalkaline magmas along the Italian peninsula; h) the absence of igneous activity where it might well be expected (e.g., above the hanging-wall of the Late Cretaceous–Eocene Adria–Europe subduction system in the Alps); i) voluminous production of subduction-related magmas coeval with extensional tectonic régimes (e.g., during Oligo-Miocene Sardinian Trough formation).To summarize, these salient central-western Mediterranean features, characterizing a late-stage of the classic ‘Wilson Cycle’ offer a ‘template’ for interpreting magmatic compositions in analogous settings elsewhere.     

Vertical distribution features of atmospheric water vapour in the Po valley area

Two sets of radiosounding measurements, taken at different hours from some stations in the Po Valley area, were examined in order to obtain the values of precipitable waterw and determine the shapes of the vertical distribution curves of absolute humidity. From these curves, we calculated the bestfit values of scale heightsH ₁ andH ₂ in the 0÷5 km and 5÷10 km altitude intervals, respectively. The analysis of the radiosounding data clearly shows that parametersw andH ₁ are closely related to the evolutionary features of the meteorological conditions on synoptic scale and are influenced, to a lesser extent, by the diurnal variations in the atmospheric ground layer. Seasonal average curves of temperature and absolute humidity, as functions of altitude, were also defined from the radiosounding measurements taken at various hours of the day. Moreover, interpolation methods in time (applied to a 12-hour range) and in space (range of about 300 km) were proposed for determining parametersw andH ₁ from the radiosounding measurements taken at different hours from the same station or from different stations at the same hour. Reliability tests, made by comparing the same station or from different stations at the same hour. Reliability tests, made by comparing the values given by the interpolation methods with those directly obtained from the radiosounding measurements, show that these evaluations ofw andH ₁ are affected by standard errors of estimate, which are comparable to the errors usually made in analyzing the radiosounding data.     

Nonlinear wave coupling in type IV solar radio bursts

In order to explain a fine structure of parallel ridges in stationary type IV continua, the emission due to the coupling of electrostatic upper hybrid waves and Bernstein waves at the sum frequency of the upper hybrid and harmonics of the gyro frequency has been calculated. If the energy density of these electrostatic waves is of the order of 10^-3 times the thermal energy density, then the observed zebra pattern can be emitted by a region with a diameter of 10³ km.     

Line emission from accretion discs around black holes: the analytic approach

The broad X-ray iron line observed in many active galactic nuclei spectra is thought to originate from the accretion disc surrounding the putative supermassive black hole. We show here how to perform the analytical integration of the geodesic equations that describe the photon trajectories in the general case of a rotating black hole (Kerr metric), in order to write a fast and efficient numerical code for modelling emission line profiles from accretion discs.     

Ephemeral rollover points and clinothem evolution in the modern Po Delta based on repeated bathymetric surveys

Reconstructions of ancient delta systems rely typically on a two-dimensional (2D) view of prograding clinothems but may miss their three-dimensional (3D) stratigraphic complexity which can, instead, be best documented on modern delta systems by integrating high-resolution geophysical data, historical cartography, core data and geomorphological reconstructions offshore. We quantitatively compare three precisely positioned, high-resolution multi-beam bathymetry maps in the delta front and pro delta sectors (0.3 to 10 m water depth) of Po di Pila, the most active of the modern Po Delta five branches. By investigating the detailed morphology of the prograding modern Po Delta, we shed new light on the mechanisms that control the topset to foreset transition in clinothems and show the temporal and spatial complexity of a delta and its pro delta slope, under the impact of oceanographic processes. This study documents the ephemeral nature of the rollover point at the transition between sandy topset (fluvial, delta plain to mouth-bar) and muddy seaward-dipping foreset deposits advancing, in this case, in >20 m of water depth. Three multibeam surveys, acquired between 2013 and 2016, document the complexity in space and time of the topset and foreset regions and their related morphology, a diagnostic feature that could not be appreciated using solely 2D, even very high-resolution, seismic profiles. In addition, the comparison of bathymetric surveys gathered with one-year lapses shows the migration of subaqueous sand dunes on the clinothem topset, the formation of ephemeral cut-and-fill features at the rollover point (few m below mean sea level), the presence of collapse depressions derived by sagging of sediments and fluid expulsion (possibly induced by storm waves) on the foreset, and splays of sand likely reflecting gravity flows on the lower foreset. Though the modern Po Delta is anthropogenic in many respects, its subaqueous clinothem can be studied as a scale model for ancient clinothems that are less resolved geometrically and far less constrained chronologically.     

Geographic variations in shell growth rates of the mussel Diplodon chilensis from temperate lakes of Chile: Implications for biodiversity conservation

The Chilean lake district includes diverse lentic ecosystems along ca. 700 km of the country (36°–43°S), including the “Nahuelbutan lakes”, “Araucanian lakes” and “Chiloe lakes”. This area is recognized as an important “hot spot” of benthic freshwater biodiversity in Southern South America. In Chilean temperate lakes, increased nutrient loads of P and N caused eutrophication, particularly in the Nahuelbutan Lakes. The freshwater Hyriidae mussel Diplodon chilensis (Gray, 1828) which is one of the most abundant species in Chilean temperate lakes, is known to be very susceptible to eutrophication. This species presents a clear reduction in its geographic ranges and is considered to be a threatened species in many Chilean lakes. In this study, we used a correlative approach to determine how eutrophication-driven changes in the food supply and in geographical parameters of different Chilean lakes affected the shell growth rates of D. chilensis. The results obtained from sclerochronological analyses of the mussel shells suggest an association with a group of environmental variables, including geographical types (negative), such as latitude and altitude, and limnological types (positive), especially phosphorous and turbidity. However, the D. chilensis populations under extreme conditions of turbidity in eutrophic and hypertrophic lakes are extinct or nearly so. The high positive correlation of the mean D. chilensis growth rates with orthophosphate (R=0.76; P<0.05), in relation to dissolved inorganic nitrogen, suggests that P is the major limiting factor of the primary productivity in Chilean temperate lakes. We discuss some implications of our results in terms of the conservation of biodiversity in temperate lake ecosystems at different taxonomic levels.     

A numerical method for the study of shear band propagation in soft rocks

This paper investigates the possibility of interpreting progressive shear failure in hard soils and soft rocks as the result of shear propagation of a pre‐existing natural defect. This is done through the application of the principles of fracture mechanics, a slip‐weakening model (SWM) being used to simulate the non‐linear zone at the tips of the discontinuity. A numerical implementation of the SWM in a computation method based on the boundary element technique of the displacement discontinuity method (DDM) is presented. The crack and the non‐linear zone at the advancing tip are represented through a set of elements, where the displacement discontinuity (DD) in the tangential direction is determined on the basis of a friction law. A residual friction angle is assumed on the crack elements. Shear resistance decreases on elements in the non‐linear zone from a peak value at the tip, which is characteristic of intact material, to the residual value. The simulation of a uniaxial compressive test in plane strain conditions is carried out to exemplify the numerical methodology. The results emphasize the role played by the critical DD on the mechanical behaviour of the specimen. A validation of the model is shown through the back analysis of some experimental observations. The results of this back analysis show that a non‐linear fracture mechanics approach seems very promising to simulate experimental results, in particular with regards to the shear band evolution pattern. Copyright © 2009 John Wiley & Sons, Ltd.     

The Tauern Window (Eastern Alps, Austria): a new tectonic map, with cross-sections and a tectonometamorphic synthesis

We present a tectonic map of the Tauern Window and surrounding units (Eastern Alps, Austria), combined with a series of crustal-scale cross-sections parallel and perpendicular to the Alpine orogen. This compilation, largely based on literature data and completed by own investigations, reveals that the present-day structure of the Tauern Window is primarily characterized by a crustal-scale duplex, the Venediger Duplex (Venediger Nappe system), formed during the Oligocene, and overprinted by doming and lateral extrusion during the Miocene. This severe Miocene overprint was most probably triggered by the indentation of the Southalpine Units east of the Giudicarie Belt, initiating at 23–21 Ma and linked to a lithosphere-scale reorganization of the geometry of mantle slabs. A kinematic reconstruction shows that accretion of European lithosphere and oceanic domains to the Adriatic (Austroalpine) upper plate, accompanied by high-pressure overprint of some of the units of the Tauern Window, has a long history, starting in Turonian time (around 90 Ma) and culminating in Lutetian to Bartonian time (45–37 Ma).     

A Generalized Backward Euler algorithm for the numerical integration of an isotropic hardening elastoplastic model for mechanical and chemical degradation of bonded geomaterials

An extended version of the classical Generalized Backward Euler (GBE) algorithm is proposed for the numerical integration of a three‐invariant isotropic‐hardening elastoplastic model for cemented soils or weak rocks undergoing mechanical and non‐mechanical degradation processes. The restriction to isotropy allows to formulate the return mapping algorithm in the space of principal elastic strains. In this way, an efficient and robust integration scheme is developed which can be applied to relatively complex yield surface and plastic potential functions. Moreover, the proposed algorithm can be linearized in closed form, thus allowing for quadratic convergence in the global Newton iteration. A series of numerical experiments are performed to illustrate the accuracy and convergence properties of the algorithm. Selected results from a finite element analysis of a circular footing on a soft rock layer undergoing chemical weathering are then presented to illustrate the algorithm performance at the boundary value problem level. Copyright © 2002 John Wiley & Sons, Ltd.     

Fundamental physics and absolute positioning metrology with the MAGIA lunar orbiter

MAGIA is a mission approved by the Italian Space Agency (ASI) for Phase A study. Using a single large-diameter laser retroreflector, a large laser retroreflector array and an atomic clock onboard MAGIA we propose to perform several fundamental physics and absolute positioning metrology experiments: VESPUCCI, an improved test of the gravitational redshift in the Earth?CMoon system predicted by General Relativity; MoonLIGHT-P, a precursor test of a second generation Lunar Laser Ranging (LLR) payload for precision gravity and lunar science measurements under development for NASA, ASI and robotic missions of the proposed International Lunar Network (ILN); Selenocenter (the center of mass of the Moon), the determination of the position of the Moon center of mass with respect to the International Terrestrial Reference Frame/System (ITRF/ITRS); this will be compared to the one from Apollo and Lunokhod retroreflectors on the surface; MapRef, the absolute referencing of MAGIA??s lunar altimetry, gravity and geochemical maps with respect to the ITRF/ITRS. The absolute positioning of MAGIA will be achieved thanks to: (1) the laboratory characterization of the retroreflector performance at INFN-LNF; (2) the precision tracking by the International Laser Ranging Service (ILRS), which gives two fundamental contributions to the ITRF/ITRS, i.e. the metrological definition of the geocenter (the Earth center of mass) and of the scale of length; (3) the radio science and accelerometer payloads; (4) support by the ASI Space Geodesy Center in Matera, Italy. Future ILN geodetic nodes equipped with MoonLIGHT and the Apollo/Lunokhod retroreflectors will become the first realization of the International Moon Reference Frame (IMRF), the lunar analog of the ITRF.