Parameters affecting maximum fluid transport in large aperture fractures

We present results of laboratory experiments to study the behavior of liquids moving in unsaturated wide-aperture fractures. A 5-mm-thick glass plate cut with a 1.7-mm aperture was used as a fractured rock analog to study behavior of film and capillary droplet flow modes. Flow rates ranged between 0.6 and 6.0 ml/min. Variability in the ambient barometric pressure, resulting from weather conditions, seemed to play a role in the natural selection of flow mode. For droplet mode, constant input conditions resulted in highly variable transport properties within the fracture. The advancing meniscus exhibited a dynamic contact angle that was a function of the droplet speed and much larger than the static contact angle. Flow rate was reduced due to the larger contact angle. Analytical expressions for droplet velocity and flow capacity are presented as a function of the dynamic rather than the static contact angle.     

The tectonic regime in Italy inferred from borehole breakout data

In Italy, the horizontal stress directions are well constrained in many regions, but the tectonic regime is not well known because the stress magnitudes are unknown. Our intention is to improve the knowledge of crustal stress in Italy, both at shallow depth and in low seismicity areas. Therefore, we inferred the tectonic regime from the comparison between the depth of breakout occurrence and the physical properties of the rocks in 20 boreholes. The critical value of the maximum horizontal stress, for which the effective tangential stress at the borehole wall overcomes the rock strength to form breakouts, could be computed from rock strength and density. Comparing the theoretical stress distributions for different tectonic regimes with the depth distribution of breakout occurrence, it is possible to infer the tectonic regime that fits best to the breakout depth distribution. We investigated boreholes up to 6 km deep located in different tectonic environments over the Italian peninsula: the Po Plain, the Apenninic chain, the Adriatic foredeep and the Tyrrhenian Quaternary volcanic region. These wells are characterised by breakout data of good quality (A, B and C, according to World Stress Map quality ranking system). The results are in general agreement with the style of faulting derived from earthquake focal mechanisms and other stress indicators. Our results show a predominance of a normal faulting (NF) regime in the inner Apennines and both normal faulting and strike–slip faulting (SS) style in the surrounding regions, possibly also associated with changes in the tectonic regime with depth.     

Numerical models of tectonic deformation at the Baltica–Avalonia transition zone during the Paleocene phase of inversion

Predictions from dynamic modelling of the lithospheric deformation are presented for Northern Europe, where several basins underwent inversion during the Late Cretaceous and Early Cenozoic and contemporary uplift and erosion of sediments occurred. In order to analyse the evolution of the continental lithosphere, the equations for the deformation of a continuum are solved numerically under thin sheet assumption for the lithosphere. The most important stress sources are assumed to be the Late Cretaceous Alpine tectonics; localized rheological heterogeneities can also affect local deformation and stress patterns. Present-day observations available in the studied region and coming from seismic structural interpretations and stress measurements have been used to constrain the model. Our modelling results show that lateral variation in lithospheric strength below the basin systems in Central Europe strongly controls the regional deformation and the stress regime. Furthermore, we have demonstrated that the geometry of the boundary between Baltica and Avalonia, together with different rheological characteristics of the two plates, had a crucial role on local crustal deformation and faulting regime resulting in the Baltica–Avalonia transition zone from the S–N Alpine convergence.     

Oxygen isotope composition and geothermometry of granulite to greenschist facies metamorphic rocks: a study from the Neoproterozoic collision-related nappe system, south of São Francisco Craton, SE Brazil

Oxygen isotope studies were carried out across units of a Neoproterozoic nappe system, south of São Francisco Craton. A temperature decrease toward the base of the system is found, consistent with a previously recognized inverted metamorphic pattern. The tectonic contact of the basal unit and the reworked southern São Francisco craton show a steep temperature gradient, suggesting that low temperature thrusting acted as the dominant tectonic process. The contrasts between the δ¹⁸O values of the Três Pontas-Varginha and Carmo da Cachoeira nappes and the differences among the samples and minerals are consistent with the preservation of sedimentary isotopic composition during metamorphism. The small differences in the δ¹⁸O values between the undeformed and the deformed calc-silicate samples (1.6‰) suggest that the δ¹⁸O value of mylonitization fluids was close to that which equilibrated with the metamorphic assemblage. The distinct δ¹⁸O values of metapelitic and calc-silicate samples and the great temperature difference from one type to the other indicate that no large-scale fluid interaction processes occurred during metamorphism. Oxygen isotopic estimations of both Três Pontas-Varginha undeformed rocks and Carmo da Cachoeira unaltered equivalents indicate δ¹⁸O values of up to 18‰. Comparison between these values and those from the ‘basement’ orthogneisses (8.3–8.5‰) indicates the latter are not sources for the metapelites.     

Palynological Evidence for Climatic Change and Human Activity during the Holocene on Minorca (Balearic Islands)

Four pollen diagrams from Minorca (Balearic Islands) have been correlated with other previously studied sequences from Majorca and Minorca to define a Holocene landscape sequence for the region from 8000 yr B.P. to the present. The lower part of the pollen diagrams reflects a climatic phase with more rain and less-marked seasonality than today. Significant quantities ofCorylus, Buxus,and mesophilous taxa are found. In the middle part, between 5000 and 4000 yr B.P., a strong change is recorded in composition and structure of the vegetational landscape, with vegetation appearing that was adapted to Mediterranean conditions. This episode coincided with the first human colonization of the island and also with a widespread climatic change in the western Mediterranean region. The change in taxa was complex and some sclerophyllous taxa suchOleaplayed an important role in the transformation of the landscape physiognomy from the mid-Holocene until the present. Although human activities have removed much of the Mediterranean vegetation on the Balearic Islands, it seems clear that the changes have been brought about, in part, by increasing dryness.     

Introduction to ground surface self-potential tomography

A new approach to self-potential (SP) data interpretation for the recognition of a buried causative SP source system is presented. The general model considered is characterized by the presence of primary electric sources or sinks, located within any complex resistivity structure with a flat air-earth boundary. First, using physical considerations of the nature of the electric potential generated by any arbitrary distribution of primary source charges and the related secondary induced charges over the buried resistivity discontinuity planes, a general formula is derived for the potential and the electric field component along any fixed direction on the ground surface. The total effect is written as a sum of elementary contributions, all of the same simple mathematical form. It is then demonstrated that the total electric power associated with the standing natural electric field component can be written in the space domain as a sum of cross-correlation integrals between the observed component of the total electric field and the component of the field due to each single constitutive elementary charge. By means of the cross-correlation bounding inequality, the concept of a scanning function is introduced as the key to the new interpretation procedure. In the space domain, the scanning function is the unit strength electric field component generated by an elementary positive charge. Next, the concept of charge occurrence probability is introduced as a suitable function for the tomographic imaging of the charge distribution geometry underground. This function is defined as the cross-correlation product of the total observed electric field component and the scanning function, divided by the square root of the product of the respective variances. Using this physical scheme, the tomographic procedure is described. It consists of scanning the section, through any SP survey profile, by the unit strength elementary charge, which is given a regular grid of space coordinates within the section, at each point of which the charge occurrence probability function is calculated. The complete set of calculated grid values can be used to draw contour lines in order to single out the zones of highest probability of concentrations of polarized, primary and secondary electric charges. An extension to the wavenumber domain and to three-dimensional tomography is also presented and discussed. A few simple synthetic examples are given to demonstrate the resolution power of the new SP inversion procedure.     

The organisation of fish assemblages in the regulated Lima basin, Northern Portugal

In order to understand the structure of fish assemblages in the modified Lima basin (Northern Portugal), two distinct datasets concerning the presence and abundance of fish species were subjected to multivariate analysis. On the River Lima two types of flow modification are present within kilometres of one another: (a) a reduced and constant flow due to hypolimnetic release; and (b) an intense and irregular flow. A comparison of their influence on fish assemblages revealed a gradient of assemblage types from tributaries to main river sites. The latter were characterised by a strong dominance of cyprinids, particularly Iberian barbel (Barbus bocagei). The former harboured two kinds of fish assemblages: those closer to the river mouth were dominated by the cyprinids Iberian chub (Squalius carolitertii) and Iberian nase (Chondrostoma polylepis), which were also frequently present in the main river; while in those further upstream the predominant species was the brown trout (Salmo trutta). Although explanatory variables such as distance from source, altitude, substrate coarseness and width were the primary correlates of fish assemblage composition, dam construction and flow regulation also had a significant effect upon assemblage structure, particularly by: i) reducing the importance of migratory species; ii) constraining the presence of trout in the regulated segments; and iii) simplifying the community, especially in the case of the constant and reduced flow regime.     

LIDAR Data Filtering and DTM Interpolation Within GRASS

LIDAR (Light Detection and Ranging) is one of the most recent technologies in surveying and mapping. LIDAR is based on the combination of three different data collection tools: a laser scanner mounted on an aircraft, a Global Positioning System (GPS) used in phase differential kinematic modality to provide the sensor position and an Inertial Navigation System (INS) to provide the orientation. The laser sends towards the ground an infrared signal, which is reflected back to the sensor. The time employed by the signal, given the aircraft position and attitude, allows computation of the earth point elevation. In standard conditions, taking into account the flight (speed 200–250 km/hour, altitude 500–2,000 m) and sensor characteristics (scan angle ± 10–20 degrees, emission rate 2,000–50,000 pulses per second), earth elevations are collected within a density of one point every 0.5–3 m. The technology allows us therefore to obtain very accurate (5–20 cm) and high resolution Digital Surface Models (DSM). For many applications, the Digital Terrain Model (DTM) is needed: we have to automatically detect and discard from the previous DSM all the features (buildings, trees, etc.) present on the terrain. This paper describes a procedure that has been implemented within GRASS to construct DTMs from LIDAR source data.