A field‐scale infiltration model accounting for spatial heterogeneity of rainfall and soil saturated hydraulic conductivity

This study first explores the role of spatial heterogeneity, in both the saturated hydraulic conductivity K_s and rainfall intensity r, on the integrated hydrological response of a natural slope. On this basis, a mathematical model for estimating the expected areal‐average infiltration is then formulated. Both K_s and r are considered as random variables with assessed probability density functions. The model relies upon a semi‐analytical component, which describes the directly infiltrated rainfall, and an empirical component, which accounts further for the infiltration of surface water running downslope into pervious soils (the run‐on effect). Monte Carlo simulations over a clay loam soil and a sandy loam soil were performed for constructing the ensemble averages of field‐scale infiltration used for model validation. The model produced very accurate estimates of the expected field‐scale infiltration rate, as well as of the outflow generated by significant rainfall events. Furthermore, the two model components were found to interact appropriately for different weights of the two infiltration mechanisms involved. Copyright © 2005 John Wiley & Sons, Ltd.     

Relativistic Cosmology with Zero Deceleration Parameter

The possibility that the deceleration parameter q might be a null constant is discussed; such possibilitiy is interesting because solves the horizon problem and the flatness problem with no need of inflation. A simple way to get q = const = 0 is explored: the way assumes Einstein's field equations without cosmological term and introduces a massless scalar field V with negative energy density. Both in the early and in the present universe one finds the Whitrow-Randall relation Gut² = const 1 G gravitational coupling, u mass-energy density, t cosmic time). The interaction between the V-field and the ordinary matter is briefly discussed; as possible consequence of this interaction the true value of the Hubble parameter might be one half the observed value. This revised version was published online in July 2006 with corrections to the Cover Date.     

The dynamics of a double-cell hydrothermal system in triggering seismicity at Somma-Vesuvius: results from a high-resolution radon survey (revisited)

Data collected at Somma-Vesuvius during the 1998–1999 radon surveys have been revisited and reinterpreted in light of recent geophysical and geochemical information. The duration of selected radon anomalies, together with the decay properties of radon, have been used to estimate the permeability and porosity of rocks of the deep hydrothermal system. The current local cyclic seismicity is explained by means of a double convective-cell model. Convective cells are separated by a low-permeability horizon located at about 2–2.5 km below sea level. Fluids convecting within the upper cells show temperatures ranging 300–350°C. Rock permeabilities in this sector are estimated on the order of 10⁻¹² m², for porosities (ϕ) of about 10⁻⁵ typical of a brittle environment where fluid velocities may reach ∼800 m/day. Fluid temperatures within the lower cells may be as high as 400–450°C, consistent with supercritical regimes. The hydrodynamic parameters for these cells are lower, with permeability k ∼ 10⁻¹⁵ m², and porosity ranging from 10⁻⁶ to 10⁻⁷. Here, fluid motion toward the surface is controlled by the fracture network within a porous medium approaching brittle–ductile behaviour, and fluid velocities may reach ∼1,800 m/day. The low-permeability horizon is a layer where upper and lower convecting cells converge. In this region, fluids (convecting both at upper and lower levels) percolate through the wallrock and release their brines. Due to self-sealing processes, permeability within this horizon reaches critical values to keep the fluid pressure near lithostatic pressure (for k ∼ 10⁻¹⁸ m²). Deep fluid pressure buildups precede the onset of hydrothermally induced earthquakes. Permeability distribution and rock strength do not exclude that the next eruption at Somma-Vesuvius could be preceded by a seismic crisis, eventually leading to a precursory phreatic explosion. The coupling of these mechanisms has the potential of inducing pervasive failure within rocks of the hydrothermal shell, and may be a prelude to a magmatic eruption. It is finally emphasised that the integrated analysis of seismic and geochemical data, including radon emissions, could be successfully used in testing temperature distributions and variations of porosity and permeability in active geothermal reservoirs.     

Non‐impact origin of the crater‐like structures in the Gilf Kebir area (Egypt): Implications for the geology of eastern Sahara

Abstract— Several small crater‐like structures occur in Gilf Kebir region (SW Egypt). It has been previously suggested that they could be the result of meteoritic impacts. Here we outline the results of our geological and geophysical survey in the area. The proposed impact origin for these structures is not supported by our observations and analyses, and we suggest an alternative interpretation. The crater‐like structures in Gilf Kebir area are likely related to endogenic processes typical of hydrothermal vent complexes in volcanic areas which may reflect the emplacement of subvolcanic intrusives.     

Weak lensing signal in unified dark matter models

This paper is an extension of the work done by Pierens & Nelson in which they investigated the behaviour of a two-planet system embedded in a protoplanetary disc. They put a Jupiter mass gas giant on the internal orbit and a lower mass planet on the external one. We consider here a similar problem taking into account a gas giant with mass in the range 0.5 to  1 M _J  and a Super-Earth (i.e. a planet with mass  ≤10 M _⊕  ) as the outermost planet. By changing disc parameters and planet masses, we have succeeded in getting the convergent migration of the planets which allows for the possibility of their resonant locking. However, in the case in which the gas giant has the mass of Jupiter, before any mean-motion first-order commensurability could be achieved, the Super-Earth is caught in a trap when it is very close to the edge of the gap opened by the giant planet. This confirms the result obtained by Pierens & Nelson in their simulations. Additionally, we have found that, in a very thin disc, an apsidal resonance is observed in the system if the Super-Earth is captured in the trap. Moreover, the eccentricity of the small planet remains low, while the eccentricity of the gas giant increases slightly due to the imbalance between Lindblad and corotational resonances. We have also extended the work of Pierens & Nelson by studying analogous systems in which the gas giant is allowed to take sub-Jupiter masses. In this case, after conducting an extensive survey over all possible parameters, we have succeeded in getting the 1:2 mean-motion resonant configuration only in a disc with low aspect ratio and low surface density. However, the resonance is maintained just for a few thousand orbits. Thus, we conclude that for typical protoplanetary discs the mean-motion commensurabilities are rare if the Super-Earth is located on the external orbit relative to the gas giant.     

Geoinformatics in morphological study of River Paglia,Tiber River basin,Central Italy

This work reports the results of a geomorphological study on the River Paglia (Tiber basin, Central Italy), to analyse the historical evolution of the riverbed/floodplain system in its lower valley reaches, upstream the confluence of River Tiber. A morphological–sedimentological approach has been used, starting with an historical analysis and a survey of the forms and sediments in the riverbed/alluvial plain system, in order to reconstruct the current evolutionary trend and to identify the causes and processes leading to changes in natural trends. This approach is not in contrast to the hydrologic–hydraulic one, normally applied in engineering fields, but it is integrated with it in a complementary and parallel manner, in order to achieve the complete knowledge of the river system, result of an investigation multidisciplinary, aimed to define its physical state. Historical data (cartographic documents, digital terrain models—DTM, topographic surveys, aerial photographs, etc.) indicate clear-cut narrowing and deepening of the active channel, manifested after the Second World War and due to the considerable changes which affected the entire fluvial system (from the recovery of land for agriculture near the river, to sediments taken from its bed over the past few decades). The result is a state of great disequilibrium, which is shown in the passage from a braided morphology, observed in the 1950s, to a single channel with low sinuosity (wandering type). This type of morphology is certainly more unstable: in fact, the narrowing and deepening of fluvial sections provide a better hydraulic efficiency for the upstream, so the peak discharge does not overflow. Such peak discharge arrives downstream, next to the confluence with the Tiber River, not laminated, causing serious flooding effects. An increased vulnerability is observed in the last reach, in the areas pertaining to the river, where most of the industrial-economic activities are concentrated and where transportation infrastructures of nationwide importance (Italy’s most important motorway, the A1, and the Rome–Florence–Milan high-speed rail links) are located. The evolution of the Paglia (which is not substantially different from that of many streams in Italy and in general throughout Europe) is affected by severe anthropic constraints and intense exploitation of resources, which have produced a state of disequilibrium approaching irreversibility in an already vulnerable system. Finally, some suggestions for the correct management of the river system are pointed out, with the aim of mitigating the flood risk in the lower Paglia valley.     

Xenoliths in recent basaltic andesite flows from Arenal Volcano,Costa Rica: inference on the composition of the lower crust

Basaltic andesite flows erupted between 1973 and 1980 from Arenal Volcano contain abundant inclusions of anorthosite, olivine gabbro, and pyroxenites, and megacrysts of olivine and anorthite. The anorthosites with large (20 mm) anorthite grains (An_96-92) exhibit deformation twinning and granulation between grain boundaries. Some olivine gabbros have angular clasts of anorthite with bent twins, pyroxene, and olivine in a finer-grained matrix which is distinctly foliated. These textural features suggest that these inclusions were deformed. An exotic (xenolithic) origin is supported in part by the mineral compositions and the estimated temperatures of equilibration: a temperature of about 975° C is obtained by two-pyroxene and Fe-Ti oxide geothermometers for the gabbros, but two-pyroxene temperatures are higher (1064 to 1120° C) for the basaltic andesite host. The olivine gabbro is thought to have crystallized at a pressure between 8.5 and 9.5 kb; whereas the lava phenocrysts crystallized at a much lower pressure of less than 5 kb. These xenoliths probably represent fragments of the lower crust below Arenal volcano. The lava flows show evidence for some contamination especially from fragments of anorthite broken apart from the larger megacrysts and xenoliths. A few phenocrysts of plagioclase in the lava samples have deformation twins. The unusually high Al₂O₃ content (19.4 to 23.2 wt%) of the lava samples can be attributed directly to the addition of anorthite; in fact the observed chemical variation in the lava flows (the increasing alumina and lime contents with decreasing silica) can be explained by this contamination.     

Some considerations about the simulation of breach channel erosion on landslide dams

The analysis of the flood hazard related to the areas downstream of landslide dams is one of the most interesting aspects of studying the formation and the failure of natural dams. The BREACH code [14], simulating the collapse of earthen dams, both man-made and naturally formed by a landslide, was chosen in order to analyse the case of the Valderchia landslide (central Italy). The bed-load transport formula used in BREACH (Meyer-Peter and Muller, modified by Smart [27]) is based on flume experiments with well-sorted sediments. Such a methodology probably makes this equation not very suitable for describing the sediment transport peculiar to a landslide body presenting a very poor material sorting. The Schoklitsch [26] formula was implemented into the programme as an alternative to the Smart equation. However, because the landslide deposits may often have a strongly bimodal grain–size frequency curve, the percentile D ₅₀ (the typical granulometric parameter requested by bed-load sediment transport formulas) can sometimes correspond to one of the grain-size classes which are really present to a lesser degree. To consider this phenomenon, the BREACH programme (version 7/88-1) was implemented with a new procedure that calculates two granulometric curves, one for each mode of the original distribution, and evaluates transport of the landslide material separately. Results of the analysis show that the model is very sensitive to the bed-load equation and that the procedure implemented to consider the eventual bimodal distribution of the dam material simulates the armouring phenomenon (which can stop the erosion of the dam during the overtopping phase).     

Structural evolution of the Pleistocene Cimini trachytic volcanic complex (Central Italy)

Structural, geomorphological, geophysical and volcanological data have been processed for the implementation of a dedicated GIS through which the structural evolution of the Pleistocene trachytic Cimini volcano (central Italy) has been reconstructed. The evolution of the Cimini complex includes three main close-in time phases: (1) intrusion of a shallow laccolith, rising along NW and NE trending faults and stagnating at the contact between the Mesozoic-Cenozoic and the Pliocene-Pleistocene sedimentary units constituting the bedrock of the volcano; (2) emplacement of lava domes along radial and tangential fractures formed by the swelling induced by the laccolith growth; (3) ignimbrite eruptions and final effusion of olivine-latitic lavas. Domes are both of Pelean and low lava dome type and their morphology was controlled by the location on the inclined surface of the swelled area. Some domes show to have uplifted upper Pliocene thermally metamorphosed clay sediments, suggesting a cryptodome-like growth. Comparison of the top of the Mesozoic-Cenozoic units with the top of the upper Pliocene-Pleistocene sedimentary complex, suggests that the laccolith emplaced in a graben of the Mesozoic-Cenozoic sedimentary complex filled by the Pliocene–Pleistocene sediments uplifted by the shallow intrusion. Stress patterns acting on the Cimini area have been deduced analysing the drainage network and the morphotectonic lineaments. Rose diagrams show a large dispersion of the lineaments reflecting the local presence of radial and tangential fractures. The most frequent extensional NW and NE trending lineaments have regional significance and controlled the magma uprise leading to the laccolith emplacement.     

A simple and efficient procedure for the numerical simulation of wind fields in complex terrain

In spite of recent progress in the prognostic numerical simulation of the atmospheric boundary layer, the explicit simulation of turbulent flows in actual complex terrain is generally still very complicated and time consuming for many environmental applications. In an attempt to develop simpler and more efficient application oriented techniques, although less refined, we propose a multi-step procedure for simulating wind fields. Once obtained the necessary meteorological input, the mass-consistent modelling technique is used to perform high-resolution mean wind flow simulations taking into account recent developments in the atmospheric boundary-layer theory. Besides, a procedure based on a generalisation of the local logarithmic law-of-the-wall over complex terrain is used to estimate the effective parameters characterising the simulated wind profiles. Turbulence intensities and spectral properties are then calculated through the estimated effective parameters, in particular through the effective friction velocity parameter. Finally, time series of the instantaneous velocity field are simulated by the Monte Carlo technique. Two applications of the proposed approach are discussed briefly: the first one is related to a coastal area in southern Italy (the Messina Straits), where the construction of the world’s longest central span bridge is being planned; the second one corresponds to the flow in a mountainous area in northern Italy (the Albenga Airport).