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.     

Exploring the impact of diagenesis on (isotope) geochemical and microstructural alteration features in biogenic aragonite

For the Quaternary and Neogene, aragonitic biogenic and abiogenic carbonates are frequently exploited as archives of their environment. Conversely, pre‐Neogene aragonite is often diagenetically altered and calcite archives are studied instead. Nevertheless, the exact sequence of diagenetic processes and products is difficult to disclose from naturally altered material. Here, experiments were performed to understand biogenic aragonite alteration processes and products. Shell subsamples of the bivalve Arctica islandica were exposed to hydrothermal alteration. Thermal boundary conditions were set at 100°C, 175°C and 200°C. These comparably high temperatures were chosen to shorten experimental durations. Subsamples were exposed to different ¹⁸O‐depleted fluids for durations between two and twenty weeks. Alteration was documented using X‐ray diffraction, cathodoluminescence, fluorescence and scanning electron microscopy, as well as conventional and clumped isotope analyses. Experiments performed at 100°C show redistribution and darkening of organic matter, but lack evidence for diagenetic alteration, except in Δ₄₇ which show the effects of annealing processes. At 175°C, valves undergo significant aragonite to calcite transformation and neomorphism. The δ¹⁸O signature supports transformation via dissolution and reprecipitation, but isotopic exchange is limited by fluid migration through the subsamples. Individual growth increments in these subsamples exhibit bright orange luminescence. At 200°C, valves are fully transformed to calcite and exhibit purple‐blue luminescence with orange bands. The δ¹⁸O and Δ₄₇ signatures reveal exchange with the aqueous fluid, whereas δ¹³C remains unaltered in all experiments, indicating a carbonate‐buffered system. Clumped isotope temperatures in high‐temperature experiments show compositions in broad agreement with the measured temperature. Experimentally induced alteration patterns are comparable with individual features present in Pleistocene shells. This study represents a significant step towards sequential analysis of diagenetic features in biogenic aragonites and sheds light on reaction times and threshold limits. The limitations of a study restricted to a single test organism are acknowledged and call for refined follow‐up experiments.     

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).