Laboratory sandbox validation of pollutant source location methods

Inverse methods can be used to recover the pollutant source location from concentration data. In this paper, the relative effectiveness of two proposed methods, simultaneous release function and source location identification (SRSI) and backward probability model based on adjoint state method (BPM-ASM) are evaluated using real data collected by using experimental equipment. The device is a sandbox that reproduces an unconfined aquifer in which all the variables are controlled. A numerical model was calibrated using experimental observations. The SRSI is a stochastic procedure which finds the source location and the release history by means of a Bayesian geostatistical approach (GA). The BPM-ASM provides the backward probability location of the pollutant detected at a monitoring point by means of a reverse transport simulation. The results show that both methods perform well. While the simultaneous release function and SRSI method requires a preliminary delineation of a probable source area and some weak hypotheses about the statistical structure of the unknown release function, the backward probability model requires some hypothesis about the contaminant release time. A case study was performed using two observation points only, and despite the scarcity of data, both methodologies were able to accurately reconstruct the true source location. The GA has the advantage to recover the release history function too, whilst the backward probability model works well with fewer data. If there are many observations, both methodologies may be computationally heavy. A transfer function approach has been adopted for the numerical definition of the sensitivity matrix in the SRSI method. The reliability of the experimental equipment was tested in previous laboratory works, conducted under several different conditions.     

Compressibility of Sands of Various Geologic Origins at Pre-crushing Stress Levels

This study quantifies the influence of various intrinsic soil properties including particle roundness, R, sphericity, S, 50% size by weight, D ₅₀, coefficient of uniformity, C _u, and the state property of relative density, D _r, on the compression and recompression indices, C _c and C _r, of sands of various geologic origins at pre-crushing stress levels. Twenty-four sands exhibiting a wide range of particle shapes, gradations, and geologic origins were collected for the study. The particle shapes were determined using a computational geometry algorithm which allows characterization of a statistically large number of particles in specimens. One dimensional oedometer tests were performed on the soils. The new data was augmented with many previously published results. Through statistical analyses, simple functional relationships are developed for C _c and C _r. In both cases, the models utilized only R and D _r since other intrinsic properties proved to have lesser direct influence on the compression indices. However, previous studies showed that the contributions of S and C _u are felt through their effects on index packing void ratios and thus on D _r. The accuracy of the models was confirmed by comparison of predicted and observed C _c and C _r values.     

Flood effects on phosphorus immobilisation in a river water filled pit lake—Case study Lake Goitsche (Germany)

Between 1999 and 2002, a former open-cast mine was filled with river water forming the recent Lake Goitsche. During filling initially acid water was neutralised. Phosphorus (P) imported from Mulde River was nearly completely removed from the water column by co-precipitation with iron (Fe) and aluminium (Al) and deposited in the sediment.During extremely high waters of the Mulde River in 2002, a dike breach facilitated a second high import of P into Lake Goitsche with suspended and dissolved matter. The analysis of total phosphorus (TP), however, showed that P again had been eliminated from the water body a few months after the flood event. Sediment investigations before filling with river water, during filling, and after the flood event were used to analyse the process of P immobilisation in a lake with acid mine drainage history.The ratios of Fe to soluble reactive P (SRP) of sediment pore water were up to three orders of magnitudes higher than in natural lakes and can serve as an indicator for potential internal P loading from sediments. The SRP concentrations at the oxic/anoxic boundary were near or below the limit of quantification (< 0.2 μmol/L). Fe and manganese (Mn) redox cycling were responsible for hindering P dissolution from sediment to lake water.Finally it can be stated, that the risk of eutrophication for such a lake seems to be low.     

Flow time estimation with spatially variable hillslope velocity in ungauged basins

The automated spatial estimation of the hillslope runoff dynamics is used as a valuable tool for the estimation of the travel time distribution (flow time), a major factor for the hydrologic prediction in ungauged basins. In fact, while the flow time function is usually obtained by rescaling the flow paths with constant channel and hillslope velocities, in this work a spatially distributed kinematic component, as a function of terrain features and in particular slope and land use, is implemented and its influence on the hydrologic response is tested by means of the Width Function Instantaneous Unit Hydrograph (WFIUH) framework. Hillslope surface flow velocities are evaluated by applying different uniform flow formulas within an automated DEM-based (terrain analysis) algorithm. A comparison test of the performances of the Manning, Darcy, Maidment and Soil Conservation Service uniform flow equations is performed for several case studies in Italy pertaining to different climatic and geomorphic conditions. Results provide new insights for a better understanding of the flow time function also introducing a more parsimonious and physically-based calibration scheme of the WFIUH.     

The 12 May 1802 earthquake (N Italy) in its historical and seismological context

The Mw 5.7 earthquake that occurred on 12 May 1802 is the only one with Mw ≥ 5.5 located west of Lake Garda in the central-northern part of the Po Plain, Northern Italy, and the strongest event located in the seismic zone 907 of the ZS9 seismogenic zonation of Italy. Current parametric earthquake catalogs locate the event not far from important cities (e.g., Milan) and to sites where nuclear power plants were to be built in the 1980s or could be built in a near future. Although the earthquake parameters seemed sufficiently well constrained, a detailed investigation of documentary sources was performed, in repositories storing the documents of the Napoleonic departments to which the area affected by the earthquake belonged at that time. In the surviving archival series, we found the officers’ correspondence on all the administrative aspects raised by the earthquake. The newly collected records allowed the authors to significantly increase the number of macroseismic intensity data, including new observations in the most damaged area. The results have been then interpreted in terms of both Mercalli–Cancani–Sieberg and EMS98 macroseismic scales. The earthquake parameters were derived applying two different methods in order to get two independent estimates. Earthquake location is confirmed, although the still scarce data available in the area to the east of the epicenter do not permit to reduce the uncertainty to a minimum. According to the Boxer method, the magnitude is now slightly higher, and the source model shows a good agreement with the tectonic setting of the area.     

A numerical study on seismic risk assessment of historic masonry towers: a case study in San Gimignano

After a brief introduction on the research project RiSEM (Seismic Risk of Monumental Buildings), the paper discusses on the seismic assessment of historic masonry towers according to the Italian “Guidelines for the assessment and mitigation of the seismic risk of the cultural heritage” that identifies an analysis methodology based on three different levels of evaluation, according to an increased knowledge of the structure. The RiSEM project aimed at developing and testing innovative and expeditious methodologies (i.e. either without or with a minimal direct contact with the masonry building) to evaluate all the main structural features of the monumental buildings required for the assessment of their seismic safety. As a relevant case study the “Town of Fine Towers”, San Gimignano (Italy), listed under the UNESCO World Cultural Heritage Sites, was selected. The paper summarizes the analyses performed on one of the San Gimignano towers: the Coppi-Campatelli one. The seismic vulnerability of the tower was evaluated with reference to the above guidelines, and the paper reports and critically compares the results obtained for the three levels of evaluation there defined: LV1 (analysis at territorial level), LV2 (local analysis) and LV3 (global analysis).     

Travertine deposition and faulting: the fault-related travertine fissure-ridge at Terme S. Giovanni,Rapolano Terme (Italy)

In this paper we describe an example of travertine fissure-ridge development along the trace of a normal fault with metre displacement, located in the eastern margin of the Neogene–Quaternary Siena Basin, in the Terme S. Giovanni area (Rapolano Terme, Italy). This morphotectonic feature, 250 m long, 30 m wide and 10 m high, formed from supersaturated hot waters (39.9°C) flowing from thermal springs aligned along the trace of the normal fault dissecting travertines not older than Late Pleistocene (24 ± 3 ka). A straight, continuous fissure with a maximum width of 20 cm occurs at the top of the ridge, along its crest. Hot fluids flow from cones mainly located at the extremities of the ridge, where travertine is depositing. The travertine fissure-ridge shows an asymmetrical profile since it buries the fault scarp. The difference in height of slopes corresponds to the vertical displacement of the normal fault. Fissuring of the recent travertine deposits along the strike of the crestal fissure, as well as recent hydrothermal circulation, lead us to believe that the Terme S. Giovanni normal fault may be currently active. On the whole, the Terme S. Giovanni fissure ridge can be defined as a travertine fault trace fissure-ridge, adding a helpful example for studying the relationship between faulting and travertine deposition.