Field observations of the June 30, 2001 debris flow at Acquabona (Dolomites, Italy)

On June 30, 2001, a debris flow occurred in the Acquabona Creek, a small catchment of the Eastern Dolomites, Italy. This debris flow originated shortly after an intense rainstorm, characterised by a peak intensity of 8.6 mm per 10 min; it transported a total volume of 30,000 m³, consisting of poorly sorted gravely sand with boulders up to 3 m in diameter. The sediment erosion yield rate reached as high as 20 m³/m. In order to verify the accuracy of the field measurements, the total volume of debris deposits have was calculated using three different topographic measurement techniques: 3D laser scanning, terrestrial stereo-photogrammetry survey and total topographic station survey. Data collected so far show that no debris flow has occurred at Acquabona with a rainfall intensity lower than 4.6 mm per 10 min. Channel cross section measurements indicate that debris flow velocity ranges from 2.0 to 7.2 m/s along the lower flow channel and peak discharge ranges between 22 and 300 m³/s. Field estimates of the rheological properties indicate a yield strength ranging from 2,088 to 5,313 Pa and Bingham viscosity between 70 and 337 Pa · s. It is not still possible to identify a rainfall intensity and amount threshold for debris flow triggering, but the data so far collected emphasise that debris flows do not occur with a rainfall intensity lower than 4.6 mm per 10 min.     

Land subsidence due to gas/oil production in inhomogeneous transversally anisotropic half-space by a boundary element method

In a previous paper¹ the authors have developed and implemented a new boundary element (BE) model to simulate and predict land subsidence occurring over three-dimensional gas/oil fields in a homogeneous and isotropic half-space. The approach relies on Betti's reciprocal theorem and makes use of the classical fundamental solution of Boussinesq in the framework of the theory of linear poroelasticity. The BE method is here extended to inhomogeneous, transversally anisotropic soils by the aid of a two-dimensional finite element (FE) model which provides a fundamental numerical solution for the actual multi-layer setting of the subsurface system. The new FE–BE approach is then used to simulate the subsidence caused by gas production over the deep reservoir of Campo Ravenna Terra, Ravenna (Italy) from 1950 to 1980. The results compare very favourably with the outcome from a full more expensive three-dimensional FE model of the same occurrence.     

A boundary element solution to land subsidence above 3-D gas/oil reservoirs

A linear boundary element (BE) model is proposed for the uncoupied simulation of land subsidence due to gas, oil and hot water production over three-dimensional (3-D) arbitrarily shaped reservoirs. The pore pressure decline is assumed to be specified in advance, e.g. via a numerical model of flow. Use is made of the fundamental solution derived in 1885 by Boussinesq for a vertical load acting upon the traction-free surface of a semi-infinite medium. A straightforward application of Betti's (1872) reciprocal theorem allows for the development of a boundary integral whose numerical execution yields directly the downward settlement over the point of interest. The new procedure is applied to assess land sinking caused by an uniform pore pressure decline occurring within fields of elliptical shape and to explore the influence of the assumption of small reservoir thickness which underlies the ‘tension center’ or ‘strain nucleus’ approach previously developed by Geertsma in 1966. The results emphasize the numerical efficiency of the solution and the promising features of the BE method for the evaluation of ground subsidence in 3-D problems. The present model is based on the theory of the linear poroelasticity and is implemented for a mechanically homogeneous and isotropic half-space. It allows for any arbitrary geometry of the reservoir and for a non-uniform distribution of the pore pressure decline. It may easily be extended to other physical settings for which a vertical surface point load solution is available.     

Tree water deficit and dynamic source water partitioning

The stable isotopes of hydrogen and oxygen (δ²H and δ¹⁸O, respectively) have been widely used to investigate tree water source partitioning. These tracers have shed new light on patterns of tree water use in time and space. However, there are several limiting factors to this methodology (e.g., the difficult assessment of isotope fractionation in trees, and the labor-intensity associated with the collection of significant sample sizes) and the use of isotopes alone has not been enough to provide a mechanistic understanding of source water partitioning. Here, we combine isotope data in xylem and soil water with measurements of tree's physiological information including tree water deficit (TWD), fine root distribution, and soil matric potential, to investigate the mechanism driving tree water source partitioning. We used a 2 m³ lysimeter with willow trees (Salix viminalis) planted within, to conduct a high spatial–temporal resolution experiment. TWD provided an integrated response of plant water status to water supply and demand. The combined isotopic and TWD measurement showed that short-term variation (within days) in source water partitioning is determined mainly by plant hydraulic response to changes in soil matric potential. We observed changes in the relationship between soil matric potential and TWD that are matched by shifts in source water partitioning. Our results show that tree water use is a dynamic process on the time scale of days. These findings demonstrate tree's plasticity to water supply over days can be identified with high-resolution measurements of plant water status. Our results further support that root distribution alone is not an indicator of water uptake dynamics. Overall, we show that combining physiological measurements with traditional isotope tracing can reveal mechanistic insights into plant responses to changing environmental conditions.     

Characterization of landslide ground surface kinematics from terrestrial laser scanning and strain field computation

Assessment and mitigation of the risk induced by landslide activation need an appropriate phenomenon investigation, to obtain useful information about the failure processes. The first step is the complete kinematics characterization of the landslide ground surface, by evaluating the involved displacement and deformation patterns. A dense displacement field can be obtained from comparison of a series of multi-temporal observations performed by means of terrestrial laser scanning. Subsequently, the strain field can be computed from displacement vectors. In this paper, a modified least square technique is employed to compute the strain on the nodes of a regular grid (2D approach) or on the points of a digital terrain model (3D approach). Such a computation takes into account the displacements, their spatial distribution, as well as the measurement and modelling errors. A scale factor is introduced in order to emphasize the contributions of the experimental points on the basis of their distance from each computation point, and to recognize possible scale-depending behaviours. This method has been implemented in Matlab and applied on two landslides located in the northeastern Italian Alps (Lamosano and Perarolo di Cadore). The experiments show that different kinematics can be recognized, and the presence and influence of eventual discontinuities can be revealed.     

Hydrological drivers of wetland vegetation community distribution within Everglades National Park,Florida

The influence of hydrological dynamics on vegetation distribution and the structuring of wetland environments is of growing interest as wetlands are modified by human action and the increasing threat from climate change. Hydrological properties have long been considered a driving force in structuring wetland communities. We link hydrological dynamics with vegetation distribution across Everglades National Park (ENP) using two publicly available datasets to study the probability structure of the frequency, duration, and depth of inundation events along with their relationship to vegetation distribution. This study is among the first to show hydrologic structuring of vegetation communities at wide spatial and temporal scales, as results indicate that the percentage of time a location is inundated and its mean depth are the principal structuring variables to which individual communities respond. For example, sawgrass, the most abundant vegetation type within the ENP, is found across a wide range of time inundated percentages and mean depths. Meanwhile, other communities like pine savanna or red mangrove scrub are more restricted in their distribution and found disproportionately at particular depths and inundations. These results, along with the probabilistic structure of hydropatterns, potentially allow for the evaluation of climate change impacts on wetland vegetation community structure and distribution.     

New gravity maps of the Eastern Alps and significance for the crustal structures

The deep seismic profile Transalp crosses, from north to south, Germany, Austria and Italy. The gravity measurements for each country were made by national agencies with different reference systems and data reduction methods. Within the frame of the Transalp-project a comprehensive database of the Eastern Alps was compiled covering an area of 3.5° by 4° in longitude and latitude (275 by 445 km), respectively. To increase the data coverage in the south Alpine area two gravity surveys were carried out, resulting in 469 areally distributed new stations, of which 215 have been measured with the intent to improve the geoid in the area of the planned Brenner Basistunnel (BBT). The resulting gravity database is the best in terms of resolution and data quality presently available for the Eastern Alps. Here the free air, Bouguer and isostatic gravity fields are critically discussed. The spatial density of existing gravity stations in the three countries is discussed. On the Italian side of the Alps the spatial density is rather sparse compared to the Austrian side. The Bouguer-gravity field varies between − 190 * 10^− 5 m/s² and + 25 * 10^− 5 m/s², with the minimum located along the Alpine high topographic chain, but with a small offset (a few tens of km) to the greatest topographic elevation, showing that the Airy-type local isostatic equilibrium does not fully apply here. The maximum of the Bouguer anomaly has an elongated shape of 100 by 50 km located between the towns of Verona and Vicenza and covers the Venetian Tertiary Volcanic Province (VTVP), a feature not directly related to the plate collision in the Eastern Alps. The gravity high is only partly explainable by high-density magmatic rocks and requires also a deeper source, like a shallowing of the Moho. The isostatic residual anomalies (Airy model) are in the range ± 50 * 10^− 5 m/s², with the greatest positive anomaly corresponding to the location of the VTVP, indicating here under-compensation of masses. At last a discussion of a 2D density model based on reflection seismic data and receiver functions is made.     

Field study on drainage densities and rescaled width functions in a high‐altitude alpine catchment

This paper addresses the effect of accurately mapping spatially heterogeneous drainage densities in high‐altitude alpine basins on Rescaled Width Functions (RWFs), used in some applications as a minimalist model of the hydrologic response. The channel network and 373 of its channel heads were mapped in the field in a high mountain catchment in the Swiss Alps. The mapped channel network is characterized by highly uneven drainage density, here described by the distribution of the length to the first channelized site computed along steepest descent from any unchannelled site. Various channel networks were extracted from a 1 m lidar‐derived digital terrain model and compared with the field‐mapped channel network using geomorphologic parameters, hillslope‐to‐channel distance and RWFs. Our results show that the channel network derived by statistical analysis of surface morphology is consistent with the field‐mapped network. Larger discrepancies were observed when the channel network was obtained with classical threshold‐based approaches relying on cumulative drainage area and local slope. The actual arrangement of the drainage densities has a significant impact on the RWFs. The discrepancy was largest between RWFs derived from classical extraction methods and RWFs derived with the field‐mapped network, indicating an inappropriate extraction of the channelled portion of the high‐altitude catchment that is a reflection of the variety of channel initiation processes. Our results suggest that spatial heterogeneity of the drainage density might play an important role in modelling streamflow generation. Copyright © 2016 John Wiley & Sons, Ltd.     

Alternative conceptual models and the robustness of groundwater management scenarios in the multi-aquifer system of the Central Veneto Basin, Italy

A three-dimensional (3D) groundwater flow model of the deep multi-aquifer Quaternary deposits of the Po plain sedimentary basin, within a 3,300-km² area (Veneto, Italy), is developed, tested and applied to aid sustainable large-scale water-resources management. Integrated mathematical modelling proves significantly successful, owing to an unusual wealth of available geological, geophysical, and hydrologic data and to state-of-the-art numerical tools. Of particular interest is the evaluation of the influence of alternative conceptual models; that is, of reconstructed representations of the 3D geological model of the structure of the aquifers. The reference conceptual model is set up by means of extended geological sections and stratigraphic records, and is used to create a large, unstructured 3D finite element grid. By analyzing the effects on piezometric surfaces and on the overall water budget of geometrical perturbations from the reference structure, alternative geo-structural models, obtained by systematically shifting the pinch-out of the aquitards, are compared. Interestingly, the impacts of aquitard pinch-out prove far from negligible. The results suggest the critical importance of reliable geological models even for large, complex 3D models. The good practice of iteratively testing numerically the impact of surprises on the conceptual model, as more field information is collected, is thus supported.