土体工程地质层组的划分

工程地质岩组的划分是岩体质量评价的基础。当大型工程揭露土层很多时,进行土体工程地质层组的划分亦是十分必要的。它可以更为清晰的反映场地地质结构的变化规律,便于地质模型的概化和参数研究,有助于抓住主要工程地质问题进行评价,便于设计师的理解和应用。以南京地铁南北线一期工程为例,探讨了土体工程地质岩组的划分问题。通过以层组为基础的参数离散性统计、工程地质模型与工程地质问题的分析可以看出该层组划分方案的合理性与工程意义。     

Ecohydrological modelling with EcH2O‐iso to quantify forest and grassland effects on water partitioning and flux ages

We used the new process‐based, tracer‐aided ecohydrological model EcH₂O‐iso to assess the effects of vegetation cover on water balance partitioning and associated flux ages under temperate deciduous beech forest (F) and grassland (G) at an intensively monitored site in Northern Germany. Unique, multicriteria calibration, based on measured components of energy balance, hydrological function and biomass accumulation, resulted in good simulations reproducing measured soil surface temperatures, soil water content, transpiration, and biomass production. Model results showed the forest “used” more water than the grassland; of 620 mm average annual precipitation, losses were higher through interception (29% under F, 16% for G) and combined soil evaporation and transpiration (59% F, 47% G). Consequently, groundwater (GW) recharge was enhanced under grassland at 37% (~225 mm) of precipitation compared with 12% (~73 mm) for forest. The model tracked the ages of water in different storage compartments and associated fluxes. In shallow soil horizons, the average ages of soil water fluxes and evaporation were similar in both plots (~1.5 months), though transpiration and GW recharge were older under forest (~6 months compared with ~3 months for transpiration, and ~12 months compared with ~10 months for GW). Flux tracking using measured chloride data as a conservative tracer provided independent support for the modelling results, though highlighted effects of uncertainties in forest partitioning of evaporation and transpiration. By tracking storage—flux—age interactions under different land covers, EcH₂O‐iso could quantify the effects of vegetation on water partitioning and age distributions. Given the likelihood of drier, warmer summers, such models can help assess the implications of land use for water resource availability to inform debates over building landscape resilience to climate change. Better conceptualization of soil water mixing processes and improved calibration data on leaf area index and root distribution appear obvious respective modelling and data needs for improved simulations.     

Isotope mixing models require individual isotopic tracer content for correct quantification of sediment source contributions

The use of isotopic tracers for sediment source apportionment is gaining interest with recent introduction of compound‐specific stable isotope tracers. The method relies on linear mixing of source isotopic tracers, and deconvolution of a sediment mixture initially quantifies the contribution of sources to the mixture's tracer signature. Therefore, a correction to obtain real sediment source proportions is subsequently required. As far as we are aware, all published studies to date have used total isotopic tracer content or a proxy (e.g., soil carbon content) for this post‐unmixing correction. However, as the relationship between the isotopic tracer mixture and the source mixture is different for each isotopic tracer, post‐unmixing corrections cannot be carried out with one single factor. This contribution presents an isotopic tracer model structure—the concentration‐dependent isotope mixing model (CD‐IMM)—to overcome this limitation. Herein, we aim to clarify why the “conventional” approach to converting isotopic tracer proportions to source proportions using a single factor is wrong. In an initial mathematical assessment, error incurred by not using CD‐IMM (NCD‐IMM) in unmixing two sources with two isotopic tracers showed a complex relation as a function of relative tracer contents. Next, three artificial mixtures with different proportions of three soil sources were prepared and deconvoluted using ¹³C of fatty acids using CD‐IMM and NCD‐IMM. Using NCD‐IMM affected both accuracy (mean average error increased up to a threefold compared with the CD‐IMM output) and precision (interquartile range was up to 2.5 times larger). Finally, as an illustrative example, the proportional source contribution reported in a published study was recalculated using CD‐IMM. This resulted in changes in estimated source proportions and associated uncertainties. Content of isotopic tracers is seldom reported in published work concerning use of isotopic tracers for sediment source partitioning. The magnitude of errors made by miscalculation in former studies is therefore difficult to assess. With this contribution, we hope the community will acknowledge the limitations of prior approaches and use a CD‐IMM in future studies.     

A comparison between the single ring pressure infiltrometer and simplified falling head techniques

Testing the relative performances of the single ring pressure infiltrometer (PI) and simplified falling head (SFH) techniques to determine the field saturated soil hydraulic conductivity, K_fs, at the near point scale may help to better establish the usability of these techniques for interpreting and simulating hydrological processes. A sampling of 10 Sicilian sites showed that the measured K_fs was generally higher with the SFH technique than the PI one, with statistically significant differences by a factor varying from 3 to 192, depending on the site. A short experiment with the SFH technique yielded higher K_fs values because a longer experiment with the PI probably promoted short‐term swelling phenomena reducing macroporosity. Moreover, the PI device likely altered the infiltration surface at the beginning of the run, particularly in the less stable soils, where soil particle arrangement may be expected to vary upon wetting. This interpretation was supported by a soil structure stability index, SSI, and also by the hydraulic conductivity data obtained with the tension infiltrometer, i.e. with a practically negligible disturbance of the sampled soil surface. In particular, a statistically significant, increasing relationship with SSI and an unsaturated conductivity greater than the saturated one were only detected for the K_fs data obtained with the PI. The SFH and PI techniques should be expected to yield more similar results in relatively rigid porous media (low percentages of fine particles and structurally stable soils) than in soils that modify appreciably their particle arrangement upon wetting. The simultaneous use of the two techniques may allow to improve K_fs determination in soils that change their hydrodynamic behaviour during a runoff producing rainfall event. Copyright © 2013 John Wiley & Sons, Ltd.     

Impact of the 6 April 2009 L'Aquila earthquake on groundwater flow in the Gran Sasso carbonate aquifer,Central Italy

The Mw = 6·3 L'Aquila earthquake on 6 April 2009 produced a mainshock that caused significant changes in the hydrogeology of the Gran Sasso carbonate fractured aquifer: (i) the sudden disappearance at the time of the mainshock of some springs located exactly along the surface trace of the Paganica normal fault (PF); (ii) an immediate increase in the discharge of the Gran Sasso highway tunnel drainages and of other springs and (iii) a progressive increase of the water table elevation at the boundary of the Gran Sasso aquifer during the following months. Using the data collected since the 1990s that include aftershock monitoring as well as data regarding spring discharge, water table elevations, turbidity and rainfall events, a conceptual model of the earthquake's consequences on the Gran Sasso aquifer is proposed herein. In this model that excludes the contribution of seasonal recharge, the short‐term hydrologic effects registered immediately after the mainshock are determined to have been caused by a pore pressure increase related to aquifer deformation. Mid‐term effects observed in the months following the mainshock suggest that there was a change in groundwater hydrodynamics. Supplementary groundwater that flows towards aquifer boundaries and springs in discharge areas reflects a possible increase in hydraulic conductivity in the recharge area, nearby the earthquake fault zone. This increase can be attributed to fracture clearing and/or dilatancy. Simulations by numerical modelling, related to pore pressure and permeability changes with time, show results in accordance with observed field data, supporting the conceptual model and confirming the processes that influenced the answer of the Gran Sasso aquifer to the L'Aquila earthquake. Copyright © 2010 John Wiley & Sons, Ltd.     

Assessment of river flow with significant lateral inflow through reverse routing modeling

The discharge hydrograph estimation in rivers based on reverse routing modeling and using only water level data at two gauged sections is here extended to the most general case of significant lateral flow contribution, without needing to deploy rainfall–runoff procedures. The proposed methodology solves the Saint‐Venant equations in diffusive form also involving the lateral contribution using a “head‐driven” modeling approach where lateral inflow is assumed to be function of the water level at the tributary junction. The procedure allows to assess the discharge hydrograph at ends of a selected river reach with significant lateral inflow, starting from the stage recorded there and without needing rainfall data. Specifically, the MAST 1D hydraulic model is applied to solve the diffusive wave equation using the observed stage hydrograph at the upstream section as upstream boundary condition. The other required data are (a) the observed stage hydrograph at the downstream section, as benchmark for the parameter calibration, and (b) the bathymetry of the river reach, from the upstream section to a short distance after the downstream gauged section. The method is validated with different flood events observed in two river reaches with a significant intermediate basin, where reliable rating curves were available, selected along the Tiber River, in central Italy, and the Alzette River, in Luxembourg. Very good performance indices are found for the computed discharge hydrographs at both the channel ends and along the tributaries. The mean Nash‐Sutcliffe value (NS_q) at the channel ends of two rivers is found equal to 0.99 and 0.86 for the upstream and downstream sites, respectively. The procedure is also validated on a longer stretch of the Tiber River including three tributaries for which appreciable results are obtained in terms of NS_q for the computed discharge hydrographs at both the channel ends for three investigated flood events.     

Testing the use of an image‐based technique to measure gully erosion at Sparacia experimental area

The first part of this investigation was aimed at testing the use of a three‐dimensional (3D) digital terrain model and a quasi‐tridimensional (2.5D) digital elevation model obtained by a large series of oblique images of eroded channels taken from consumer un‐calibrated and non‐metric cameras. For two closed earth channels having a different sinuosity, the ground measurement of some cross sections by a profilometer (P) was carried out and their real volume was also measured. The comparison among the three methods (3D, 2.5D, and P) pointed out that a limited underestimation of the total volume always occurs and that the 3D method is characterized by the minimum difference between measured and real volume. For this reason, 3D model can be used as benchmark. In the subsequent part of the investigation, the three ground measurement methods were applied for surveying of an ephemeral gully (EG) channel at the Sparacia area. The morphological and hydraulic variable values of the 24 surveyed cross sections determined by both 2.5D model and profilometer were compared. This comparison showed that the estimate error is generally less than ±10%. The EG measurements carried out by the three methods supported the applicability both of the empirical relationship between EG length and its eroded volume and the theoretical dimensionless relationship among the morphological variables describing the channelized erosion process. Finally, it was demonstrated that the effect of the distance interval on the EG volume measurement by 3D and 2.5D models is negligible for the investigated EG.