The use of radon as tracer in environmental sciences

Radon can be used as a naturally occurring tracer for environmental processes. By means of grab-sampling or continuous monitoring of radon concentration, it is possible to assess several types of dynamic phenomena in air and water. We present a review of the use of radon and its progeny at the University of Cantabria. Radon can be an atmospheric dynamics indicator related with air mass interchange near land-sea discontinuities as well as for the study of vertical variations of air parameters (average values of different types of stability: 131–580 Bq m^?3). Concerning indoor gas, we present some results obtained at Altamira Cave (Spain): from 222 to 6549 Bq m^?3 (Hall) and from 999 to 6697 Bq m^?3 (Paintings Room). Finally, variations of radon concentration in soil (0.3 to 9.1 kBq m^?3) and underground water (values up to 500 Bq l^?1) provide relevant information about different geophysical phenomena.     

Cotton Fiber/ZrO2, A New Material for Adsorption of Cr(VI) Ions in Water

The natural cotton fiber was used to synthesize an anion exchange, containing ZrO₂ film on its surface, NCFZC (natural cotton fiber/ZrO₂ composite). This anion exchanger was produced by the reaction of the zirconium oxychloride and hydroxyl groups on surface of the natural cotton fiber. The material was used for Cr(VI) ions adsorption studies. Adsorption equilibrium time and optimum pH for Cr(VI) adsorption were found to be 6 h and 4.0, respectively. The Langmuir and Temkin isotherms were used to models adsorption equilibrium data. The adsorption capacity of NCFZC was found to be 1.33 mmol/g. Kinetic studies showed that the rate of adsorption of Cr(VI) on NCFZC obeyed a pseudo‐second‐order kinetic model.     

The 29 September 2009 Samoa Islands Tsunami: Simulations Based on the First Focal Mechanism Solutions and Implications on Tsunami Early Warning Strategies

The tsunamigenic earthquake (Mw?=?8.1) that occurred on 29 September 2009 at 17:48 UTC offshore of the Samoa archipelago east of the Tonga trench represents an example of the so-called ??outer-rise?? earthquakes. The areas most affected were the south coasts of Western and American Samoa, where almost 200 people were killed and run-up heights were measured in excess of 5?m at several locations along the coast. Moreover, tide gauge records showed a maximum peak-to-peak height of about 3.5?m near Pago Pago (American Samoa) and of 1.5?m offshore of Apia (Western Samoa). In this work, different fault models based on the focal mechanism solutions proposed by Global CMT and by USGS immediately after the 2009 Samoan earthquake are tested by comparing the near-field recorded signals (three offshore DART buoys and two coastal tide gauges) and the synthetic signals provided by the numerical simulations. The analysis points out that there are lights and shadows, in the sense that none of the computed tsunamis agrees satisfactorily with all the considered signals, although some of them reproduce some of the records quite well. This ??partial agreement?? and ??partial disagreement?? are analysed in the perspective of tsunami forecast and of Tsunami Early Warning System strategy.     

Soil-structure interaction with Rayleigh waves

Presented is a plane-strain method for soil-structure interaction analysis consisting of the superposition of the free field motions and the interaction motions in a generalized seismic environment. The free field is modelled as a horizontally layered viscoelastic medium and the seismic environment may consist of a combination of S, P and Rayleigh waves. The soil-structure system is modelled with viscoelastic finite elements, transmitting boundaries, viscous boundaries and a three-dimensional simulation. Comparative analyses of the same structure are conducted for an input of R waves and for vertically propagating S and P waves in a rock site and sand site. In the rock site the R waves produce higher peak horizontal spectral acceleration up to 25 per cent, and a significant rocking effect at points away from the centre of gravity of the structure. However, the S and P waves show higher peak vertical spectral acceleration by up to 15 per cent at the centre of the structure. Very similar horizontal response, but higher vertical response only at the centre of the structure for S and P waves, is obtained for the sand site.     

Transient processes in Stromboli’s shallow basaltic system inferred from dolerite and magmatic breccia blocks erupted during the 5 April 2003 paroxysm

We describe the mineralogy, geochemistry, and mesomicrostructure of fresh subvolcanic blocks erupted during the 5 April 2003 paroxysm of Stromboli (Aeolian Islands, Italy). These blocks represent ∼50 vol.% of the total erupted ejecta and consist of fine- to medium-grained basaltic lithotypes ranging from relatively homogeneous dolerites to strongly or poorly welded magmatic breccias. The breccia components are represented by angular fragments of dolerites entrapped in a matrix of vesiculated (lava-like to scoriae) crystal-rich (CR) basalt. All of the studied blocks are cognates with the CR basalt of the normal Strombolian activity or lavas and they are often coated by a few-centimeter thick layer of crystal-poor (CP) basaltic pumice erupted during the paroxysm. We suggest that they result from the rapid increase of pressure and related subvolcanic rock failure that occurred shortly before the 5 April 2003 explosion, when the uppermost portion of the edifice inflated and suffered brecciation as the result of the sudden rise of the gas-rich CP basalt that triggered the eruption. Dolerites and magmatic matrix of the breccias show major and trace element compositions that match those of the CR basalts erupted during normal Strombolian activity and effusive events at Stromboli volcano. Dolerites consist of (a) phenocrysts normally found in the CR basalts and (b) late-stage magmatic minerals such as sanidine, An_60-28 plagioclase, Fe–Mn-rich olivines (Fo_68-48), phlogopite, apatite, and opaque mineral pairs (magnetite and ilmenite), most of which are never found both in lava flows and scoriae erupted during the persistent explosive activity that characterizes typical Strombolian behavior. Subvolcanic crystallization of the Stromboli CR magma, leading to slowly cooled equivalents of basalts, could result from transient drainage of the magma from the summit craters to lower levels. Fingering and engulfing of the material that collapsed from the summit crater floor into the shallow basaltic system during the late evening of 28 December 2002 coupled with the short break in the summit persistent explosions between December 2002 and March 2003 permitted the CR magma pockets to solidify as dolerites, which were confined to the uppermost portion of the system and thus not involved in the ongoing flank effusive activity. Crystal size distribution of the basaltic blocks and crystallization of the finer-grained (<0.1 mm) mafic minerals of the dolerites over a time interval of ∼100 days closely agrees with the above interpretation. Vesicle filling (miarolitic cavities) locally found in some dolerites, with minerals deposited as vapor-phase crystallization is a result of continuous gas percolation through the rocks of the uppermost portion of the volcanic system. Poorly welded magmatic breccias formed during syn-eruptive processes of 5 April 2003, when the paroxysm strongly shattered the shallow subvolcanic system and many dolerite fragments were entrapped in the CR magma. In contrast, the high degree of welding between the dolerite clasts and the CR basaltic matrix in the strongly welded magmatic breccias provides a snapshot of subvolcanic intrusions of the CR basalt into the dolerite when, after a 2-month break in activity, CR magmas started to rise again to the summit craters. Blocks similar to these subvolcanic ejecta of 5 April 2003 were also erupted during previous paroxysms (e.g., 1930) suggesting that changes in the usual Strombolian activity (e.g., short breaks in the persistent mild explosions and/or flank effusive activity) lead to transient crystallization of dolerites in the shallow plumbing system.     

Controls on event runoff coefficients in the eastern Italian Alps

Analyses of event runoff coefficients provide essential insight on catchment response, particularly if a range of catchments and a range of events are compared by a single indicator. In this study we examine the effect of climate, geology, land use, flood types and initial soil moisture conditions on the distribution functions of the event runoff coefficients for a set of 14 mountainous catchments located in the eastern Italian Alps, ranging in size from 7.3 to 608.4 km². Runoff coefficients were computed from hourly precipitation, runoff data and estimates of snowmelt. A total of 535 events were analysed over the period 1989–2004. We classified each basin using a “permeability index” which was inferred from a geologic map and ranged from “low” to “high permeability”. A continuous soil moisture accounting model was applied to each catchment to classify ‘wet’ and ‘dry’ initial soil moisture conditions. The results indicate that the spatial distribution of runoff coefficients is highly correlated with mean annual precipitation, with the mean runoff coefficient increasing with mean annual precipitation. Geology, through the ‘permeability index’, is another important control on runoff coefficients for catchments with mean annual precipitation less than 1200 mm. Land use, as indexed by the SCS curve number, influences runoff coefficient distribution to a lesser degree. An analysis of the runoff coefficients by flood type indicates that runoff coefficients increase with event snowmelt. Results show that there exists an intermediate region of subsurface water storage capacity, as indexed by a flow–duration curve-based index, which maximises the impact of initial wetness conditions on the runoff coefficient. This means that the difference between runoff coefficients characterised by wet and dry initial conditions is negligible both for basins with very large storage capacity and for basins with small storage capacity. For basins with intermediate storage capacities, the impact of the initial wetness conditions may be relatively large.     

A thermo-plastic constitutive law for brittle-plastic behavior of rocks at high temperatures

Mechanical properties of rocks change under the influence of, temperature. Stress at the onset of yielding, ultimate strength, dilatancy, strain hardening and softening, and the confining pressure at brittle-ductile transition are all reduced by the increasing temperature. This study presents a framework of constitutive modeling of thermo-brittle-plastic behavior of rocks which encompasses these changes. The constitutive law is based on a thermo-plasticity theory first proposed for metals byPrager (1958). Two phenomenological mechanisms have been identified as central for the modeling: temperature dependence of the yield locus (thermal softening), and temperature dependence of the strain-hardening function (thermally enhanced ductility). Material parameters for two rocks, Carrara marble and Westerly granite, were determined on the basis of additional hypotheses. These parameters are used in numerical simulations of triaxial tests at different temperatures. The obtained stress-strain curves compare well to the experimental results. The changes with temperature in the stress at the onset of yielding are more accurately reproduced that the evolution of hardening or softening. Suggestions for possible improvements and future research directions are indicated.     

Evaluating Earthquake-Triggered Landslide Hazard at the Basin Scale Through Gis in the Upper Sele River Valley

To evaluate techniques for assessing earthquake-triggeredlandslide hazard in the Southern Apennines (Italy), a GIS-based analysis was used to modelseismically induced slope deformations. Geological, geotechnical, geomorphological and seismologicaldata were integrated into a standard earthquake slope stability model. The model assessed thelandslide potential that existed during the 1980 Irpinian earthquake in the Upper Sele river Valley.The standard Newmark displacement analysis, widely used for predicting the location of shallowunstable slopes, does not take into account errors and/or uncertainties in the input parameters.Therefore, a probabilistic Newmark displacement analysis technique has been used. Probabilistictechniques allow, e.g., an estimation of the probability that a slope will exceed a certain criticalvalue of Newmark displacement. In our probabilistic method, a Monte-Carlo based simulation modelis used in conjunction with a GIS. The random variability of geotechnical data is modelled by probabilitydensity functions (pdfs), while for the seismic input three different regression laws wereconsidered. Input probability distributions are sampled and the resulting values input into empiricalrelations for estimating Newmark displacement. The outcome is a map in which to each siteis related a spatial probability distribution for the expected displacement in response to seismic loading.Results of the experiments show a high grade of uncertainty in the application of the Newmarkanalysis both for the deterministic and probabilistic approach in a complex geological setting suchas the high Sele valley, quite common in the Southern Apennines. They show a strong dependence onthe reliability of the spatial data used in input, so that, when the model is used at basin scale,results are strongly influenced by local environmental condition (e.g., topography, lithology, groundwatercondition) and decrease the model performance.     

Volcanological and petrological evolution of Marsili Seamount (southern Tyrrhenian Sea)

A swath bathymetric survey was conducted on Marsili Volcano, the biggest seamount in the Tyrrhenian Sea. It stands 3000 m above the surrounding oceanic crust of the 3500 m-deep Marsili back-arc basin and is axially located within the basin. The seamount has an elongated shape and presents distinctive morphology, with narrow (<1000 m) ridges, made up of several elongated cones, on the summit zone and extensive cone fields on its lower flanks. A dredging campaign carried out at water depths varying between 3400 and 600 m indicates that most of Marsili Seamount is composed of medium-K calc-alkaline basalts. Evolved high-K andesites were only recovered from the small cones on the summit axis zone. Petrological and geochemical characteristics of the least differentiated basalts reveal that at least two varieties of magmas have been erupted on the Marsili Volcano. Group 1 basalts have plagioclase and olivine as dominant phases and show lower Al, Ca, K, Ba, Rb and Sr, and higher Fe, Na, Ti and Zr with respect to a second type of basaltic magma. Group 2 basalts reveal the presence of clinopyroxene as an additional phenocryst phase. In addition, the two basaltic magmas have different original pre-eruptive H₂O content (group 1, H₂O-poor and group 2, H₂O-rich). Moreover, comparison of the compositional trends and mineralogical compositions obtained from MELTS [Ghiorso, M.S., Sack, R.O., Contrib. Mineral. Petrol. 119 (1995) 197–212] fractional crystallization calculations reveal that the evolved andesites can only exclusively be derived from a low-pressure (0.3 kbar) fractionation of magmas compositionally similar to the least evolved group 2 basalts. Finally, we suggest that the high vesicularity of the basalts sampled at relatively great depths (>2400 m) on the edifice is governed by H₂O and, probably, CO₂ exsolution and is not a feature indicative of shallow water depth eruption.     

Groundwater flow model, recharge estimation and sustainability in an arid region of Patagonia, Argentina

The Península Valdés, in northeastern Patagonia, Argentina, is characterised by its arid climate and the lack of perennial watercourses; thus, all economic activities depend on the groundwater resources. Water demand is mainly associated with tourism, which is centralised in Puerto Pirámides and supplied by a water desalination plant, and to sheep farming, supplied by the local aquifer. Due to the exponential growth of tourism, the government is planning to exploit groundwater and convey it by aqueduct to the abovementioned locality. The objectives of this study were to corroborate the conceptual geohydrological model, to develop a mathematical model to simulate the response of the aquifer to different scenarios, and to assess the incidence of water input into the system as a variable—a function that poses difficulties in the models for arid regions. The Visual Modflow 4.1 code was used, calibrating it in trial-and-error mode, changing the recharge and hydraulic conductivity parameters with different variants in the recharge zone and in the inclusion or exclusion of the evapotranspiration module. Results indicate the importance of the recharge analysis by treating rainfall at daily time steps. The adjusted model was exposed to four scenarios with variations in water input and in output by pumping. It can be concluded that under different input conditions, but with a controlled extraction, the system responds in a sustainable manner.