Is the forward problem of ground water hydrology always well posed?

Complex aquifer systems are often modeled with quasi-three-dimensional models, which consider two-dimensional horizontal flow in the aquifers and one-dimensional vertical flow through aquitards. When the aquifer system consists of a phreatic aquifer and one or more semiconfined aquifers connected by aquitards, the discrete model consists of a nonlinear system of algebraic equations, because the transmissivity of the phreatic aquifer depends on the phreatic head. If the water extraction is very high, the phreatic aquifer can be depleted and the equations of the model must be modified accordingly. There are not simple and general criteria to state if the phreatic aquifer is depleted before solving the system of equations. Therefore, the iterative procedures (e.g., relaxation methods), used to find the solution to the forward problem, must handle these particular conditions and can suffer several problems of convergence. These problems can be caused by the choice of the initial head values or of the relaxation coefficient of the iterative algorithms; however, they can also be caused by the nonexistence or nonuniqueness of the solution to the system of nonlinear equations. The study of existence and uniqueness of the general problem is very difficult and, therefore, we consider a simplified problem, for which the discrete model can be handled analytically. The results of the numerical experiments show that the solution to the forward problem can be nonunique. Only for some cases it is possible to invoke physical arguments to eliminate tentative solutions.     

High-precision late Aptian Pb/U age for the Punta del Barco Formation (Baqueró Group), Santa Cruz Province,Argentina

The Baqueró Group holds a rich fossil flora including macrofloristic and palynological remains that characterize the Early Cretaceous vegetation of Argentina. The age of the upper part of the Group, known as Punta del Barco Formation, has been a topic of several discussions in the past. The type locality of this stratigraphic unit exposes fossiliferous tuffs together with muddy tuffs and tuffaceous sandstones beds. However, there have been no reliable age constraints for the fossil assemblages. ²⁰⁶Pb/²³⁸U analyses of zircon crystals yielded an age of 114.67 ± 0.18 Ma. This age control makes possible an improved analysis and comparison of the floras in the southern region.     

High-pressure thermo-elastic properties of beryl (Al<Subscript>4</Subscript>Be<Subscript>6</Subscript>Si<Subscript>12</Subscript>O<Subscript>36</Subscript>) from ab initio calculations,and observations about the source of thermal expansion

Ab initio calculations of thermo-elastic properties of beryl (Al₄Be₆Si₁₂O₃₆) have been carried out at the hybrid HF/DFT level by using the B3LYP and WC1LYP Hamiltonians. Static geometries and vibrational frequencies were calculated at different values of the unit cell volume to get static pressure and mode-γ Grüneisen’s parameters. Zero point and thermal pressures were calculated by following a standard statistical-thermodynamics approach, within the limit of the quasi-harmonic approximation, and added to the static pressure at each volume, to get the total pressure (P) as a function of both temperature (T) and cell volume (V). The resulting P(V, T) curves were fitted by appropriate EoS’, to get bulk modulus (K ₀) and its derivative (K′), at different temperatures. The calculation successfully reproduced the available experimental data concerning compressibility at room temperature (the WC1LYP Hamiltonian provided K ₀ and K′ values of 180.2 Gpa and 4.0, respectively) and the low values observed for the thermal expansion coefficient. A zone-centre soft mode \( P6/mcc \to P\bar{1} \) phase transition was predicted to occur at a pressure of about 14 GPa; the reduction of the frequency of the soft vibrational mode, as the pressure is increased, and the similar behaviour of the majority of the low-frequency modes, provided an explanation of the thermal behaviour of the crystal, which is consistent with the RUM model (Rigid Unit Model; Dove et al. in Miner Mag 59:629–639, 1995), where the negative contribution to thermal expansion is ascribed to a geometric effect connected to the tilting of rigid polyhedra in framework silicates.     

Spatial variations in the similarity of earthquake populations: The case of the 1997 Colfiorito–Sellano (northern Apennines, Italy) seismic sequence

The spatial properties of events in the 1997 Colfiorito–Sellano seismic sequence (Northern Apennines, Italy) were investigated using coherence, a parameter derived from seismic moment tensors that quantifies the kinematic similarity between focal mechanisms. The 1997 Colfiorito–Sellano seismic sequence predominantly consists of normal faulting earthquakes, with a few strike-slip and reverse faulting episodes. This kinematic heterogeneity is possibly related to the contemporaneous activity of two different sets of faults: NW–SE normal faults and NNE–SSW sub-vertical faults, the latter inherited from the previous Miocene compressional phase. The study used two independently-derived data sets of the same seismic sequence characterized by a different number of events and by different precision of spatial localisation. Their statistical significances, assessed through a reshuffling procedure, reveal that data sets with at least some hundreds of events and good positional precision are required to obtain significant results through coherence analysis. Results from the better quality data set indicate that this seismic sequence is characterized by a rapid decrease in the kinematic similarity between earthquake pairs within 2 km of separation, particularly along directions sub-perpendicular to the normal fault strike. The decrease rate seems to be controlled by the geometric characteristics of the normal faults, given that the mean along-dip distance between fault segments is 2 km. In proximity to pre-existing tectonic lineaments the relative abundance of strike-slip and reverse faults tends to decrease the kinematic similarity between events but does not influence the coherence decrease rate. The presence of mixed focal mechanisms (normal, reverse and strike-slip) in a single seismic phase implies that mixed fault types are not restricted to polyphase tectonic histories: such heterogeneous kinematics during a single phase may be induced by the presence of inherited discontinuities.     

Landslides triggered by the 23 November 2000 rainfall event in the Imperia Province, Western Liguria, Italy

From mid-October to 22 November 2000, the western Liguria Region of Italy experienced prolonged and intense rainfall, with cumulative values exceeding 1000 mm in 45 days. The severe rainfall sequence ended on November 23 with a high-intensity storm that dumped more than 180 mm of rain in 24 h. The high-intensity event caused flooding and triggered more than 1000 soils slips and debris flows and a few large, complex landslides. Slope failures caused three fatalities and severe damage to roads, private homes, and agriculture. Large (1:13,000) and very large (1:5000) scale colour aerial photographs were taken 45 days after the event over the areas most affected by the landslides. Through the interpretation of the 334 photographs covering an area of 500 km², we prepared a landslide inventory map that shows 1204 landslides, for a total landslide area of 1.6 km². We identified the rainfall conditions that triggered landslides in the Armea valley using cumulative- and continuous-rainfall data, combined with detailed information on the time of landslide occurrence. Landslide activity initiated 8 to 10 h after the beginning of the storm, and the most abundant activity occurred in response to rainfall intensities of 8 to 10 mm per hour. For the Ceriana Municipality, an area where the landslides were numerous in November 2000, we also collected information about a historical event that occurred on 8–11 December 1910 and triggered abundant landslides resulting in severe economic damage. A comparison of the damage caused by the historical and the recent landslide events indicated that damage caused by the 1910 historical event was more diffused but less costly than the damage caused by the 2000 event.     

Environmental changes in the central Po Plain (northern Italy) due to fluvial modifications and anthropogenic activities

The fluvial environment of the central Po Plain, the largest plain in Italy, is discussed in this paper. Bounded by the mountain chains of the Alps and the Apennines, this plain is a link between the Mediterranean environment and the cultural and continental influences of both western and eastern Europe. In the past decades, economic development has been responsible for many changes in the fluvial environment of the area.This paper discusses the changes in fluvial dynamics that started from Late Pleistocene and Early Holocene due to distinct climatic changes. The discussion is based on geomorphological, pedological, and archaeological evidences and radiocarbon dating.In the northern foothills, Late Pleistocene palaeochannels indicate several cases of underfit streams among the northern tributaries of the River Po. On the other hand, on the southern side of the Po Plain, no geomorphological evidence of similar discharge reduction has been found. Here, stratigraphic sections, together with archaeological remains buried under the fluvial deposits, show a reduction in the size of fluvial sediments after the 10th millennium BC. During the Holocene, fluvial sedimentation became finer, and was characterised by minor fluctuations in the rate of deposition, probably related to short and less intense climatic fluctuations.Given the high rate of population growth and the development of human activities since the Neolithic Age, human influence on fluvial dynamics, especially since the Roman Age, prevailed over other factors (i.e., climate, tectonics, vegetation, etc.). During the Holocene, the most important changes in the Po Plain were not modifications in water discharge but in sediment. From the 1st to 3rd Century AD, land grants to war veterans caused almost complete deforestation, generalised soil erosion, and maximum progradation of the River Po delta. At present, land abandonment in the mountainous region has led to reafforestation. Artificial channel control in the mountain sector of the basins and in-channel gravel extraction (now illegal but very intense in the 1960s and 1970s) are causing erosion along the rivers and along large sectors of the Adriatic coast. These changes are comparable with those occurring in basins of other Mediterranean rivers.     

Probability maps of landslide reactivation derived from tree-ring records: Pra Bellon landslide, southern French Alps

Probability maps of landslide reactivation are presented for the Pra Bellon landslide located in the southern French Alps based on results obtained with dendrogeomorphic analysis. Spatiotemporal patterns of past landslide activity was derived from tree-ring series of 403 disturbed mountain pine trees growing in the landslide body. In total, 704 growth disturbances were identified in the samples indicating 22 reactivation phases of the landslide body between 1910 and 2011. The mean return period was 4.5 years. Given the spatiotemporal completeness of the reconstruction, probabilities of landslide reactivation were computed and illustrated using a Poisson distribution model and for 5, 20, 50, and 100 years. Probability of landslide reactivation is highest in the central part of the landslide body and increases from 0.13 for a 5-year period to 0.94 for a 100-year period. Conversely, probabilities of reactivation are lower at its margins. The proposed method differs from conventional approaches based on statistical analyses or physical modeling that have demonstrated to have limitations in the prediction of spatiotemporal reactivation of landslides. Our approach is, in contrast, based on extensive data on past landslides and therefore allowed determination of quantitative probability maps of reactivation derived directly from the frequency of past events. This approach is considered a valuable tool for land managers in charge of protecting and forecasting people and their assets from the negative effects of landslides as well as for those responsible for land use planning and management. It demonstrates the reliability of dendrogeomorphic mapping that should be used systematically in forested shallow landslides.     

Contrasting radon background levels in volcanic settings: clues from 220Rn activity concentrations measured during long-term deformation experiments

To better understand the mechanisms leading to different radon background levels in volcanic settings, we have performed two long-term deformation experiments of 16 days using a real-time setup that enables us to monitor any variation of radon activity concentration during rock compression. Our measurements demonstrate that, in the case of highly porous volcanic rocks, the emanating power of the substrate changes as a function of the volcanic stress conditions. Constant magmatic pressures, such as those observed during dike intrusions and hydrothermal fluid injections, can result in pervasive pore collapse that is mirrored by a significant radon decrease until a constant emanation is achieved. Conversely, repeated cycles of stress due to, for example, volcano inflation/deflation cycles, cause a progressive radon increase a few days (but even weeks and months) before rupture. After rock failure, however, the formation of new emanation surfaces leads to a substantial increase of the radon signal. Our results suggest that surface deformation in tectonic and volcanic settings, such as inflation/deflation or constant magmatic pressures, have important repercussions on the emanating power of volcanic substrates.     

Relationships between glacier and rock glacier in the Maritime Alps, Schiantala Valley, Italy

In the Schiantala Valley of the Maritime Alps, the relationship between a till-like body and a contiguous rock glacier has been analyzed using geomorphologic, geoelectric and ice-petrographic methodologies. DC resistivity tomographies undertaken in the till and in the rock glacier show the presence of buried massive ice and ice-rich sediments, respectively. Ice samples from a massive ice outcrop show spherical gas inclusions and equidimensional ice crystals that are randomly orientated, confirming the typical petrographic characteristics of sedimentary ice. The rock glacier formation began after a phase of glacier expansion about 2550 ± 50 ¹⁴C yr BP. Further ice advance during the Little Ice Age (LIA) overrode the rock glacier root and caused partial shrinkage of the pre-existing permafrost. Finally, during the 19th and 20th centuries, the glacial surface became totally debris covered. Geomorphological and geophysical methods combined with analyses of ice structure and fabric can effectively interpret the genesis of landforms in an environment where glaciers and permafrost interact. Ice petrography proved especially useful for differentiating ice of past glaciers versus ice formed under permafrost conditions. These two mechanisms of ice formation are common in the Maritime Alps where many sites of modern rock glaciers were formerly occupied by LIA glaciers.