Geothermal assessment of the island of ischia (southern Italy) from isotopic and chemical composition of the delivered fluids

The Ischia geothermal system is hosted by silicic rocks of the Quaternary Potassic Roman Province, in southern Italy. Exploration drilling down to 1156 m depth in the mid-1950s provided information on boiling profiles (up to 250°C) and on the depth and permeability of the potential reservoirs. Discharge fluid samples were collected and analyzed to define the inflow of surrounding seawater (C1 ranges from 2.5 to 20 g/kg) into the system.Analyses of samples from surface manifestations and shallow wells collected during 1983 and 1988 point to the existence of three distinct mixing regimes, involving three water components. A dishomogeneous body of diluted water (Cl less than 2.5 g/kg), that occurs at depths > 700 m and reequilibrates at 240°C at least, is overlain by an aquifer of groundwater variably mixed with variably seawater (Cl from 4 to 10 g/kg), which tends to reequilibrate at 160°C. Steam-heated waters locally develop and act as dilutants of the rising geothermal fluids.Dilution, mixing, and evaporation of the ascending chloride fluids are supported by oxygen and hydrogen isotopic data the thermal waters being enriched in ¹⁸O and D with respect to local meteoric water by up to 7 and 30‰, respectively. The relative composition of the major cations in thermal solutions was used to discriminate the two main groups of thermal waters, the reservoir temperatures of which are estimated from the Na/K-gethermometer. K-Mg geothermometer indicates reequilibration in near-surface conditions.The isotopic composition of the fumarolic steam varies from −7 to −12‰ in ∂⁸O and from − 35 to − 70‰ in ∂D, in agreement with a deep mixed fluid that boils adiabatically from 240 to 80°C. The deuterium content of the H₂O-H₂ pair gives enrichment factor of about 830‰, corresponding to equilibrium temperature conditions slightly higher than the surface boiling temperatures. The ∂¹³C of CO₂is almost constant at −4.5‰ (1δ=0.4), suggesting an important magmatic contribution, and the ∂¹⁸O values of CO₂appears to in equilibrium with accompanying steam at the measured temperatures.The CO₂/Ar and H₂/Ar chemical ratios have been used to derive aquifer temperatures, the values obtained being consistent with those of solute geothermometers.     

Fractal distribution of particle size in carbonate cataclastic rocks from the core of a regional strike-slip fault zone

We present particle size data from 31 samples of carbonate cataclastic rocks collected across the 26 m thick fault core of the Mattinata Fault in the foreland of the Southern Apennines, Italy. Particle size distributions of incoherent samples were determined by a sieving-and-weighting technique. The number of weight-equivalent spherical particles by size is well fitted by a power-law function on a log–log space. Fractal dimensions (D) of particle size distributions are in the 2.091–2.932 range and cluster around the value of 2.5. High D-values pertain to gouge in shear bands reworking the bulk cataclastic rocks of the fault core. Low D-values characterise immature cataclastic breccias. Intermediate D-values are typical of the bulk fault core. Analysis of the ratio between corresponding equivalent particle numbers from differently evolved cataclastic rocks indicates that the development of particle size distributions with D>2.6–2.7 occurred by a preferential relative increase of fine particles rather than a selective decrement of coarse particles. This preferentially occurred in shear bands where intense comminution enhanced by slip localisation progressed by rolling of coarse particles whose consequent smoothing produced a large number of fine particles. Our data suggest that during the progression of cataclasis, the fragmentation mode changed from the Allègre et al.'s [Nature 297 (1982) 47] “pillar of strength” mechanism in the early evolutionary stages, to the Sammis et al.'s [Pure and Applied Geophysics 125 (1987) 777] “constrained comminution” mechanism in the subsequent stages of cataclasis. Eventually, localised shear bands developed mainly by abrasion of coarse particles.     

Tectono-stratigraphic modelling of the North Sicily continental margin (southern Tyrrhenian Sea)

A two-dimensional numerical modelling that simulate the kinematic and thermal response of the lithosphere to thinning was used for the quantitative reconstruction of the late Neogene to Recent times tectonic and stratigraphic evolution of the North Sicily continental margin (southern Tyrrhenian Sea). The numerical study of the evolution of the North Sicily margin builds on the crustal image and kinematic interpretation of the margin obtained by Pepe et al. [Tectonics 19 (2000) 241] on the basis of seismic data and gravity modelling. Tectonic modeling indicate that different segments of the margin were undergoing different vertical movements, which are mainly expression of the rifting and thinning of the lithosphere occurred during tectonic evolution of the southern Tyrrhenian Sea. A prediction of the pre-rift basement topography and the Moho along the margin converges to a value of 6.5 km for the depth of necking and a temperature-dependent EET (500° isotherm). The model fails to reproduce the morphology of the Solunto High confirming its non-extensional origin. A polyphase evolution is required to reproduce the observed syn- and post-rift stratigraphy. During the first rifting stage (between 9 and 5 Ma), crustal thinning factors reach maximum values of 1.27 in the Cefalù basin. A similar value is predicted for the subcrustal thinning around 60 km NNE of the profile margin. Crustal thinning factors increase during the second rifting stage (from 4 to 2 Ma) and reach values of 2 and up to 3.5 in the Cefalù basin and in the continent–oceanic transition zone, respectively. Similarly, subcrustal lithospheric thinning factors reach values up to 2.5 in the distal sector of the margin. An uplift of more than 100 m is predicted for the North Sicily shelf and surrounding onshore areas during the post-rift stage. The evolution of thermal structure with time is very sensitive to the partial thinning factors describing the evolution of the thinning itself during time. The lithosphere preserved part of its strength during extension. The effective elastic thickness (EET) along the margin through time is 24 km at the onset of rifting and reaches values less to 8 km during the second rifting stage in the northeastern end of the margin.     

Design storm prediction and hydrologic modeling using a web-GIS approach on a free-software platform

The aim of this work has been to implement a set of procedures useful to automatise the evaluation, the design storm prediction and the flood discharge associated with a selected risk level. For this purpose a Geographic Information System has been implemented using Grass 5.0.One of the main topics of such a system is a georeferenced database of the highest intensity rainfalls and their assigned duration recorded in Sicily. This database contains the main characteristics for more than 250 raingauges, as well as the values of intense rainfall events recorded by these raingauges. These data are managed through the combined use of the PostgreSQL and GRASS-GIS 5.0 databases. Some of the best-known probability distributions have been implemented within the Geographical Information System in order to determine the point and/or areal rain values once duration and return period have been defined. The system also includes a hydrological module necessary to compute the probable flow, for a selected risk level, at points chosen by the user.A peculiarity of the system is the possibility to querying the model using a web-interface. The assumption is that the rising needs of geographic information, and dealing with the rising importance of peoples participation in the decision process, requires new forms for the diffusion of territorial data. Furthermore, technicians as well as public administrators needs to get customized and specialist data to support planning, particularly in emergencies. In this perspective a Web-interface has been developed for the hydrologic system. The aim is to allow remote users to access a centralized database and processing-power to serve the needs of knowledge without complex hardware/software infrastructures.     

The unlikely independence of Northern Italy

The past fifty years have seen several elements and factors leading to the current homogenization of Italian society: the internal migration from South to North, the development of a more and more industrial and affluent society, the increase of links with Europe/World culture and economy. Nevertheless in the last decade in Northern Italy emerged the political phenomenon of autonomist parties, in particular the most successful of them: the Lega Nord. This success happened in parallel with the loss of power by traditional government parties, also as a consequence of the trial against corruption held in Milan at the beginning of the '90. Italian politics is in turmoil and it is impossible to forecast the duration or even its end. The problem is to determine the most significant elements, and which factors have the greatest effect on the current reality, also to be able to venture a prediction regarding trends. The paper suggests some ethno-nationalism evaluation criteria, useful to assess the more or less successful geopolitical dynamics where there are separatist claims. According with those criteria the Northern Italian political situation has been briefly described and evaluated, with particular attention to electoral results of the Lega Nord party. These results show a well defined geographical distribution and hence the minority condition in the North of the separatist supporters. The actual political debate in Italy has the aim to change some parts of the Constitution in the sense of federalism; it can be a sufficient reply to secessionist claims and it can make the independence of Northern Italy unlikely.     

Using GB-SAR technique to monitor slow moving landslide

A Ground-Based SAR (GB-SAR) interferometer was employed to measure the surface displacements of a landslide occurring in the Carnian Alps, north-eastern Italy, which has affected a national road and seriously damaged a road tunnel still under construction. Moreover, since the landslide is located on the left bank of the Tagliamento River Valley, it is feared that this mass movement might dam the river, creating a basin that would increase natural hazard for the valley inhabitants. The data collected from December 2002 to July 2005 by a conventional monitoring system, consisting of a GPS network and boreholes equipped with inclinometric tubes, showed that the landslide was moving at a quasi-constant rate of about 3 cm per year. Due to the slow deformation rate of the landslide, a recently developed GB-SAR technique based on the analysis of a restricted ensemble of coherent points was used. Two surveys, each lasting two days, were planned in December 2004 and July 2005, in order to map and measure the surface displacements that occurred over time. The results from the radar were compared with the ones derived from the GPS monitoring network. An agreement was achieved among the data collected, showing the capability of the GB-SAR technique to measure displacements even within a time span of several months between the surveys.     

Minerals at high pressure. Mechanics of compression from quantum mechanical calculations in a case study: the beryl (Al<Subscript>4</Subscript>Be<Subscript>6</Subscript>Si<Subscript>12</Subscript>O<Subscript>36</Subscript>)

The Bader topological analysis has been applied to ab initio computed electron densities of beryl, in order to clarify its mechanism of compression. Full structural optimization and total energy (E) calculations were performed at different cell volumes (V _c). The pressure at each volume and the equation of state were estimated from the first and second derivatives of the resultant E(V _c) curve. The total (negative) potential energy of the crystal, sum of both attractive and repulsive electrostatic terms, was found to systematically decrease (i.e., it moved to more negative values) up to the highest pressure considered (28.4 GPa), indicating that interelectronic and internuclear repulsions are not the only terms controlling the compressibility, at least in the pressure range investigated. Electronic kinetic energy increases as the cell volume is reduced, leading to a parallel increase of the total energy. Both structure at equilibrium and compressibility are therefore due to the balance between the opposing kinetic and potential energy terms. The Bader theory has been used to identify the topological atoms within the structure and to calculate their properties, with particular attention to the forces driving the structural relaxation at high pressure. On a qualitative basis, the obtained results are expected to be transferable to the discussion of compressibility of other mineral phases.     

Recharge estimation from discrete water-table datasets in a coastal shallow aquifer in a humid subtropical climate

A time-series approach to the estimation of recharge rate in unconfined aquifers of highly variable water level is proposed. The approach, which is based on the water-table fluctuation method (WTF), utilizes discrete water-level measurements. Other similar techniques require continuous measurements, which makes them impossible to apply in cases where no data from automatic loggers are available. The procedure is deployed at the Ressacada Farm site, southern Brazil, on a coastal shallow aquifer located in a humid subtropical climate where diurnal water-level variations of up to 1 m can follow a precipitation event. The effect of tidal fluctuations on the groundwater levels is analyzed using a harmonic component builder, while a time-variable drainage term is evaluated through an independent analysis and included in the assessment. The estimated recharge values are compared with those obtained from the continuous measurements showing a good agreement with the approaches for discrete dataset intervals of up to 15 days. Subsequently, the estimated recharge rates are incorporated into a transient groundwater-flow model and the water levels are compared showing a good match. Henceforth, the approach extends the applicability of WTF to noncontinuous water-level datasets in groundwater recharge studies.