Hierarchical simulation of aquifer heterogeneity: implications of different simulation settings on solute-transport modeling

The fine-scale heterogeneity of porous media affects the large-scale transport of solutes and contaminants in groundwater and it can be reproduced by means of several geostatistical simulation tools. However, including the available geological information in these tools is often cumbersome. A hierarchical simulation procedure based on a binary tree is proposed and tested on two real-world blocks of alluvial sediments, of a few cubic meters volume, that represent small-scale aquifer analogs. The procedure is implemented using the sequential indicator simulation, but it is so general that it can be adapted to various geostatistical simulation tools, improving their capability to incorporate geological information, i.e., the sedimentological and architectural characterization of heterogeneity. When compared with a standard sequential indicator approach on bi-dimensional simulations, in terms of proportions and connectivity indicators, the proposed procedure yields reliable results, closer to the reference observations. Different ensembles of three-dimensional simulations based on different hierarchical sequences are used to perform numerical experiments of conservative solute transport and to obtain ensembles of equivalent pore velocity and dispersion coefficient at the scale length of the blocks (meter). Their statistics are used to estimate the impact of the variability of the transport properties of the simulated blocks on contaminant transport modeled on bigger domains (hectometer). This is investigated with a one-dimensional transport modeling based on the Kolmogorov-Dmitriev theory of branching stochastic processes. Applying the proposed approach with diverse binary trees and different simulation settings provides a great flexibility, which is revealed by the differences in the breakthrough curves.     

Regional deformation of late Quaternary fluvial sediments in the Apennines foreland basin (Emilia,Italy)

International Journal of Earth Sciences - Our research is aimed at estimating the vertical deformation affecting late Quaternary units accumulated into the foreland basin of the Northern Apennines...     

Natural highways for end-of-life solutions in the LEO region

We present the main findings of a dynamical mapping performed in the Low Earth Orbit region. The results were obtained by propagating an extended grid of initial conditions, considering two different epochs and area-to-mass ratios, by means of a singly averaged numerical propagator. It turns out that dynamical resonances associated with high-degree geopotential harmonics, lunisolar perturbations and Solar radiation pressure can open natural deorbiting highways. For area-to-mass ratios typical of the orbiting intact objects, these corridors can be exploited only in combination with the action exerted by the atmospheric drag. For satellites equipped with an area augmentation device, we show the boundary of application of the drag, and where the Solar radiation pressure can be exploited.     

Deformation patterns of upper Quaternary strata and their relation to active tectonics,Po Basin,Italy

Despite increased application of subsurface datasets below the limits of seismic resolution, reconstructing near‐surface deformation of shallow key stratigraphic markers beneath modern alluvial and coastal plains through sediment core analysis has received little attention. Highly resolved stratigraphy of Upper Pleistocene to Holocene (Marine Isotope Stage 5e to Marine Isotope Stage 1) alluvial, deltaic and coastal depositional systems across the southern Po Plain, down to 150 m depth, provides an unambiguous documentation on the deformation of previously flat‐lying strata that goes back in time beyond the limits of morphological, historical and palaeoseismic records. Five prominent key horizons, accurately selected on the basis of their sedimentological characteristics and typified for their fossil content, were used as highly effective stratigraphic markers (M1 to M5) that can be tracked for tens of kilometres across the basin. A facies‐controlled approach tied to a robust chronology (102 radiocarbon dates) reveals considerable deformation of laterally extensive nearshore (M1), continental (M2 and M3) and lagoon (M4 and M5) marker beds originally deposited in a horizontal position (M1, M4 and M5). The areas where antiformal geometries are best observed are remarkably coincident with the axes of buried ramp anticlines, across which new seismic images reveal substantially warped stratal geometries of Lower Pleistocene strata. The striking spatial coincidence of fold crests with the epicentres of historic and instrumental seismicity suggests that deformation of marker beds M1 to M5 might reflect, in part at least, syntectonically generated relief and, thus, active tectonism. Precise identification and lateral tracing of chronologically constrained stratigraphic markers in the ¹⁴C time window through combined sedimentological and palaeoecological data may delineate late Quaternary subsurface stratigraphic architecture at an unprecedented level of detail, outlining cryptic stratal geometries at the sub‐seismic scale. This approach is highly reproducible in tectonically active Quaternary depositional systems and can help to assess patterns of active deformation in the subsurface of modern alluvial and coastal plains worldwide.     

Field Measurement of Suction, Water Content, and Water Permeability

This paper presents a review of techniques for field measurement of suction, water content, and water hydraulic conductivity (permeability). Main problems in the use of field tensiometers are addressed and hints on how to improve tensiometer performance are given. Advantages and limitations of instruments for indirect measurement of suction including electrical conductivity sensors, thermal conductivity sensors, dielectric permittivity sensors, filter paper, and psychrometer are discussed. Techniques for water content measurement based on dielectric methods are then presented. These include time and amplitude domain reflectometry and capacitance. Finally, a brief overview of methods for measurement of water permeability in the field is presented.     

LiDAR-based digital terrain analysis of an area exposed to the risk of lava flow invasion: the Zafferana Etnea territory, Mt. Etna (Italy)

The town of Zafferana Etnea, located on the southeastern slope of Mt. Etna volcano (Italy), has been repeatedly threatened by lava flows in recent centuries. The last serious threat occurred during the 1991–1993 eruption, when the lava front came to a halt only 1.7 km from the centre of town. Morphostructural data derived from light detection and ranging (LiDAR) surveys carried out on Etna in 2005 have enabled us to evaluate the risk of lava invasion in a section (16 km²) of the Zafferana Etnea territory. Qualitative and quantitative results are obtained combining the information derived from LiDAR analysis with geological, morphological and structural data using geographic information systems technology (GIS). The study quantifies in unprecedented detail the areal extent and volume of forested and urban areas and its degree of exposure to different levels of hazard from future lava invasion. Nearly 52% of the urban texture fall into areas of moderate to high risk from lava invasion. Future land use planning should take these findings into account and promote new development preferentially in areas of lower risk.     

Structural evolution of a 2M 1 phengite mica up to 11 GPa: an in situ single-crystal X-ray diffraction study

The structural evolution at high pressure of a natural 2M ₁-phengite [(K_0.98Na_0.02)_Σ=1.00(Al_1.55Mg_0.24Fe_0.21Ti_0.02)_Σ=2.01(Si_3.38Al_0.62)O₁₀(OH)₂; a = 5.228(2), b = 9.057(3), c = 19.971(6)Å, β = 95.76(2)°; space group: C2/c] from the metamorphic complex of Cima Pal (Sesia Zone, Western Alps, Italy) was studied by single-crystal X-ray diffraction with a diamond anvil cell under hydrostatic conditions up to ~11 GPa. A series of 12 structure refinements were performed at selected pressures within the P range investigated. The compressional behaviour of the same phengite sample was previously studied up to ~25 GPa by synchrotron X-ray powder diffraction, showing an irreversible transformation with a drastic decrease of the crystallinity at P > 15–17 GPa. The elastic behaviour between 0.0001 and 17 GPa was modelled by a third-order Birch–Murnaghan Equation of State (BM-EoS), yielding to K _T0 = 57.3(10) GPa and K′ = ?K _T0/?P = 6.97(24). The single-crystal structure refinements showed that the significant elastic anisotropy of the 2M ₁-phengite (with β(a):β(b):β(c) = 1:1.17:4.60) is mainly controlled by the anisotropic compression of the K-polyhedra. The evolution of the volume of the inter-layer K-polyhedron as a function of P shows a negative slope, Fitting the P–V(K-polyhedron) data with a truncated second-order BM-EoS we obtain a bulk modulus value of K _T0(K-polyhedron) = 26(1) GPa. Tetrahedra and octahedra are significantly stiffer than the K-polyhedron. Tetrahedra behave as quasi-rigid units within the P range investigated. In contrast, a monotonic decrease is observed for the octahedron volume, with K _T0 = 120(10) GPa derived by a BM-EoS. The anisotropic response to pressure of the K-polyhedron affects the P-induced deformation mechanism on the tetrahedral sheet, consisting in a cooperative rotation of the tetrahedra and producing a significant ditrigonalization of the six-membered rings. The volume of the K-polyhedron and the value of the ditrigonal rotation parameter (α) show a high negative correlation (about 93%), though a slight discontinuity is observed at P >8 GPa. α increases linearly with P up to 7–8 GPa (with ?α/?P ≈ 0.7°/GPa), whereas at higher Ps a “saturation plateau” is visible. A comparison between the main deformation mechanisms as a function of pressure observed in 2M ₁- and 3T-phengite is discussed.     

Highly reproducible quasi-mosaic crystals as optical components for a Laue lens

The realization of a Laue lens for astronomical purposes involves the mass production of a series of crystalline tiles as optical components, allowing high-efficiency diffraction and high-resolution focusing of photons. Crystals with self-standing curved diffraction planes is a valid and promising solution. Exploiting the quasi-mosaic effect, it turns out to be possible to diffract radiation at higher resolution. In this paper we present the realization of 150 quasi-mosaic Ge samples, bent by grooving one of their largest surface. We show that grooving method is a viable technique to manufacture such crystals in a simple and very reproducible way, thus compatible with mass production. Realized samples present very homogenous curvature. Furthermore, with a specific chemical etch, it is possible to fine adjust one by one the radius of curvature of the grooved samples. Realized crystals was selected for the ASI’s Laue project, that involves the implementation of a prototype of a Laue lens for hard X- and soft γ-ray astronomy.