A visual basic spreadsheet macro for geochemical background analysis

A Visual Basic macro entitled BACKGROUND calculates geochemical background values of chemical parameters and estimates threshold values separating background data from anomalies. The macro uses two statistical methods, the iterative 2-sigma technique and the calculated distribution function, and integrates these model-based objective methods into a widely accessible platform (i.e., MS Excel). The macro offers the possibility for automated processing of geochemical data and enables an automated generation of background range and threshold values for chemical parameters.     

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.     

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.     

Geophysical Finite-Element Simulation Tool (GeoFEST): Algorithms and Validation for Quasistatic Regional Faulted Crust Problems

GeoFEST (Geophysical Finite Element Simulation Tool) is a two- and three-dimensional finite element software package for the modeling of solid stress and strain in geophysical and other continuum domain applications. It is one of the featured high-performance applications of the NASA QuakeSim project. The program is targeted to be compiled and run on UNIX systems, and is running on diverse systems including sequential and message-passing parallel systems. Solution to the elliptical partial differential equations is obtained by finite element basis sampling, resulting in a sparse linear system primarily solved by conjugate gradient iteration to a tolerance level; on sequential systems a Crout factorization for the direct inversion of the linear system is also supported. The physics models supported include isotropic linear elasticity and both Newtonian and power-law viscoelasticity, via implicit quasi-static time stepping. In addition to triangular, quadrilateral, tetrahedral and hexahedral continuum elements, GeoFEST supports split-node faulting, body forces, and surface tractions. This software and related mesh refinement strategies have been validated on a variety of test cases with rigorous comparison to analytical solutions. These include a box-shaped domain with imposed motion on one surface, a pair of strike slip faults in stepover arrangement, and two community-agreed benchmark cases: a strike slip fault in an enclosing box, and a quarter-domain circular fault problem. Scientific applications of the code include the modeling of static and transient co- and post-seismic earth deformation, Earth response to glacial, atmospheric and hydrological loading, and other scenarios involving the bulk deformation of geologic media.     

Modeling the water exchanges between the Venice Lagoon and the Adriatic Sea

A hydrodynamic model of the Venice Lagoon and the Adriatic Sea has been developed in order to study the exchanges at the inlets of the Venice Lagoon, a complex morphological area connecting the sea and the lagoon. The model solves the shallow water equations on a spatial domain discretized by a staggered finite element grid. The grid represents the Adriatic Sea and the Venice Lagoon with different spatial resolutions varying from 30 m for the smallest channels of the lagoon to 30  km for the inner areas of the central Adriatic Sea. Data from more than ten tide gauges displaced in the Adriatic Sea have been used in the calibration of the simulated water levels. After the calibration, the tidal wave propagation in the North Adriatic and in the Venice Lagoon is well reproduced by the model. To validate the model results, empirical flux data measured by acoustic Doppler current profiler probes installed inside the inlets of Lido and Malamocco have been used and the exchanges through the three inlets of the Venice Lagoon have been analyzed. The comparison between modeled and measured fluxes at the inlets outlines the efficiency of the model to reproduce both tide- and wind-induced water exchanges between the sea and the lagoon. Even in complex areas, where highly varying resolution is needed, the model is suitable for the simulation of the dominating physical processes.     

Operation of two major reservoirs of Iran under IPCC scenarios during the XXI century

We assess the effects of prospective climate change until 2100 on water management of two major reservoirs of Iran, namely, Dez (3.34 × 10⁹ m³) and Alavian (6 × 10⁷ m³). We tune the Poly‐Hydro model suited for simulation of hydrological cycle in high altitude snow‐fed catchments. We assess optimal operation rules (ORs) for the reservoirs using three algorithms under dynamic and static operation and linear and non‐linear decision rules during control run (1990–2010 for Dez and 2000–2010 for Alavian). We use projected climate scenarios (plus statistical downscaling) from three general circulation models, EC‐Earth, CCSM4, and ECHAM6, and three emission scenarios, or representative concentration pathways (RCPs), RCP2.6, RCP4.5, and RCP8.5, for a grand total of nine scenarios, to mimic evolution of the hydrological cycle under future climate until 2100. We subsequently test the ORs under the future hydrological scenarios (at half century and end of century) and the need for reoptimization. Poly‐Hydro model when benchmarked against historical data well mimics the hydrological budget of both catchments, including the main processes of evapotranspiration and streamflows. Teaching–learning‐based optimization delivers the best performance in both reservoirs according to objective scores and is used for future operation. Our projections in Dez catchment depict decreased precipitation along the XXI century, with ?1% on average (of the nine scenarios) at half century and ?6% at the end of century, with changes in streamflows on average ?7% yearly and ?13% yearly, respectively. In Alavian, precipitation would decrease by ?10% on average at half century and ?13% at the end of century, with streamflows ?14% yearly and ?18% yearly, respectively. Under the projected future hydrology, reservoirs' operation would provide lower performance (i.e., larger lack of water) than now, especially for Alavian dam. Our results provide evidence of potentially decreasing water availability and less effective water management in water stressed areas like Northern Iran here during this century.     

SESTET: A spatially explicit stream temperature model based on equilibrium temperature

Stream-water temperature is a key variable controlling chemical, biological, and ecological processes in freshwater environments. Most models focus on a single river cross-section; however, temperature gradients along stretches and tributaries of a river network are crucial to assess ecohydrological features such as aquatic species suitability, growth and feeding rates, or disease transmission. We propose SESTET, a deterministic, spatially explicit stream temperature model for a whole river network, based on water and energy budgets at a reach scale and requiring only commonly available spatially distributed datasets, such as morphology and air temperature, as input. Heat exchange processes at the air–water interface are modelled via the widely used equilibrium temperature concept, whereas the effects of network structure are accounted for through advective heat fluxes. A case study was conducted on the prealpine Wigger river (Switzerland), where water temperatures have been measured in the period 2014–2018 at 11 spatially distributed locations. The results show the advantages of accounting for water and energy budgets at the reach scale for the entire river network, compared with simpler, lumped formulations. Because our approach fundamentally relies on spatially distributed air temperature fields, adequate spatial interpolation techniques that account for the effects of both elevation and thermal inversion in air temperature are key to a successful application of the model. SESTET allows the assessment of the magnitude of the various components of the heat budget at the reach scale and the derivation of reliable estimates of spatial gradients of mean daily stream temperatures for the whole catchment based on a limited number of conveniently located (viz., spanning the largest possible elevation range) measuring stations. Moreover, accounting for mixing processes and advective fluxes through the river network allows one to trust regionalized values of the parameters controlling the relationship between equilibrium and air temperature, a key feature to generalize the model to data-scarce catchments.