Changes in Zn speciation during soil formation from Zn-rich limestones

In order to better understand the long-term speciation and fractionation of Zn in soils, we investigated three soils naturally enriched in Zn (237–864 mg/kg Zn) from the weathering of Zn-rich limestones (40–207 mg/kg Zn) using extended X-ray absorption fine structure (EXAFS) spectroscopy and sequential extractions. The analysis of bulk EXAFS spectra by linear combination fitting (LCF) indicated that Zn in the oolitic limestones was mainly present as Zn-containing calcite (at site Dornach), Zn-containing goethite (Gurnigel) and Zn-containing goethite and sphalerite (Liestal). Correspondingly, extraction of the powdered rocks with 1 M NH₄-acetate at pH 6.0 mobilized only minor fractions of Zn from the Gurnigel and Liestal limestones (<30%), but most Zn from the Dornach rock (81%). In the Dornach soil, part of the Zn released from the dissolving limestone was subsequently incorporated into pedogenic hydroxy-interlayered vermiculite (Zn-HIV, 30%) and Zn-containing kaolinite (30%) and adsorbed or complexed by soil organic and inorganic components (40%). The Gurnigel and Liestal soils contained substantial amounts of Zn-containing goethite (50%) stemming from the parent rock, smaller amounts (20%) of Zn-containing kaolinite (and possibly Zn-HIV), as well as adsorbed or complexed Zn-species (30%). In the soil from Liestal, sphalerite was only found in trace amounts, indicating its dissolution during soil formation. In sequential extractions, large percentages of Zn (55–85%) were extracted in recalcitrant extraction steps, confirming that Zn-HIV, Zn-containing kaolinite and Zn-containing goethite are highly resistant to weathering. These Zn-bearing phases thus represent long-term hosts for Zn in soils over thousands of years. The capability of these phases to immobilize Zn in heavily contaminated soils may however be limited by their uptake capacity (especially HIV and kaolinite) or their abundance in soil.     

Did the 7/9/1999 M5.9 Athens Earthquake Come with a Warning?

Prior to the 7/9/1999 M_S = 5.9 Athens earthquake, regional seismicity has exhibited a power-law increase, of the form = K+A(t_c - t)ⁿ, where is estimated using an expression log = cM + d and t_c is the time of the culminating event. Such changes appeared after the 17/8/1999 M7.4 Izmit event. We quantified the performance of the power law vs. the null hypothesis of constant seismic release rates, by defining the curvature C as the ratio of the power law fit RMS/linear fit RMS, so that the smaller C is, the better the power law behaviour. By mapping C, we have established a critical radius of 110 km and observed that the region of correlated accelerating seismic release extended from the N. Aegean, through Euboea and Attica to the SW Peloponnese. A few days prior to the Athens event, min(C) was centred at the epicentral area and numerical simulation yielded t_c = 1999.676 and predicted M_S = 5.77. Seismicity rates returned to normal (quasi-constant) after the Athens event. We interpret this effect as critical point behaviour, following remote excitation of a broad area by stress redistribution due to the Izmit event which, at Athens, has triggered `premature' failure of a fault nearing its load bearing capacity. If this is correct, we have documented a case of remote earthquake triggering by another earthquake, as well as insight into the mechanisms producing it. As a corollary, we note that a large event may beget another large event in its broader region of interaction, which may be preceded by characteristic precursory seismicity changes.     

Probability Density Functions as Botanical-Climatological Transfer Functions for Climate Reconstruction

We present a new procedure, the pdf method (pdf=probability density function), for reconstructing Quaternary climate utilizing botanical data. The procedure includes the advantages of the indicator species method by considering the fossil and modern presence and absence of taxa rather than their frequencies, thus avoiding the need for modern analog plant communities. Overcoming the problematic use of absolute limits to describe climate response ranges is the main progress of the pdf method in comparison to the indicator species method. This advantage results from estimating probability density functions (pdfs) for monthly mean January and July temperature conditional on the present day occurrence of single taxa. Gaussian distributions sufficiently approximate pdfs of many, although not all, studied taxa. On the assumption of statistical independence, the procedure calculates a joint pdf as the product of the pdfs of the individual taxa. This algorithm weights each taxon according to the extent of its climate response range expressed by its covariance structure. We interpret the maximum of the resulting pdf as the most likely climate and its confidence interval as the uncertainty range. To avoid an artificial reduction of uncertainty arising from the use of numerous similar pdfs, a preselection method is proposed based on the Mahalanobis distance between pdfs. The pdf method was applied to the Carpinus phase of a profile from Gröbern, Germany, that spans the last interglaciation (Eemian). The reconstructed most probable January and July temperatures of about 0.0°C and 18.4°C barely differ from the modern values of −0.5°C and 18.3°C.     

The Tauern Window (Eastern Alps, Austria): a new tectonic map, with cross-sections and a tectonometamorphic synthesis

We present a tectonic map of the Tauern Window and surrounding units (Eastern Alps, Austria), combined with a series of crustal-scale cross-sections parallel and perpendicular to the Alpine orogen. This compilation, largely based on literature data and completed by own investigations, reveals that the present-day structure of the Tauern Window is primarily characterized by a crustal-scale duplex, the Venediger Duplex (Venediger Nappe system), formed during the Oligocene, and overprinted by doming and lateral extrusion during the Miocene. This severe Miocene overprint was most probably triggered by the indentation of the Southalpine Units east of the Giudicarie Belt, initiating at 23–21 Ma and linked to a lithosphere-scale reorganization of the geometry of mantle slabs. A kinematic reconstruction shows that accretion of European lithosphere and oceanic domains to the Adriatic (Austroalpine) upper plate, accompanied by high-pressure overprint of some of the units of the Tauern Window, has a long history, starting in Turonian time (around 90 Ma) and culminating in Lutetian to Bartonian time (45–37 Ma).     

Correction to: Leveraging HPC accelerator architectures with modern techniques — hydrologic modeling on GPUs with ParFlow

Computational Geosciences - A Correction to this paper has been published: https://doi.org/10.1007/s10596-021-10065-y     

Tier-based approaches for landslide susceptibility assessment in Europe

In the framework of the European Soil Thematic Strategy and the associated proposal of a Framework Directive on the protection and sustainable use of soil, landslides were recognised as a soil threat requiring specific strategies for priority area identification, spatial hazard assessment and management. This contribution outlines the general specifications for nested, Tier-based geographical landslide zonings at small spatial scales to identify priority areas susceptible to landslides (Tier 1) and to perform quantitative susceptibility evaluations within these (Tier 2). A heuristic, synoptic-scale Tier 1 assessment exploiting a reduced set of geoenvironmental factors derived from common pan-European data sources is proposed for the European Union and adjacent countries. Evaluation of the susceptibility estimate with national-level landslide inventory data suggests that a zonation of Europe according to, e.g. morphology and climate, and performing separate susceptibility assessments per zone could give more reliable results. To improve the Tier 1 assessment, a geomorphological terrain zoning and landslide typology differentiation are then applied for France. A multivariate landslide susceptibility assessment using additional information on landslide conditioning and triggering factors, together with a historical catalogue of landslides, is proposed for Tier 2 analysis. An approach is tested for priority areas in Italy using small administrative mapping units, allowing for relating socioeconomic census data with landslide susceptibility, which is mandatory for decision making regarding the adoption of landslide prevention and mitigation measures. The paper concludes with recommendations on further work to harmonise European landslide susceptibility assessments in the context of the European Soil Thematic Strategy.     

Uncertainties of geochemical modeling during CO<Subscript>2</Subscript> sequestration applying batch equilibrium calculations

One of the uncertainties in the field of carbon dioxide capture and storage (CCS) is caused by the parameterization of geochemical models. The application of geochemical models contributes significantly to calculate the fate of the CO₂ after its injection. The choice of the thermodynamic database used, the selection of the secondary mineral assemblage as well as the option to calculate pressure dependent equilibrium constants influence the CO₂ trapping potential and trapping mechanism. Scenario analyses were conducted applying a geochemical batch equilibrium model for a virtual CO₂ injection into a saline Keuper aquifer. The amount of CO₂ which could be trapped in the formation water and in the form of carbonates was calculated using the model code PHREEQC. Thereby, four thermodynamic datasets were used to calculate the thermodynamic equilibria. Furthermore, the equilibrium constants were re-calculated with the code SUPCRT92, which also applied a pressure correction to the equilibrium constants. Varying the thermodynamic database caused a range of 61% in the amount of trapped CO₂ calculated. Simultaneously, the assemblage of secondary minerals was varied, and the potential secondary minerals dawsonite and K-mica were included in several scenarios. The selection of the secondary mineral assemblage caused a range of 74% in the calculated amount of trapped CO₂. Correcting the equilibrium constants with respect to a pressure of 125 bars had an influence of 11% on the amount of trapped CO₂. This illustrates the need for incorporating sensitivity analyses into reaction pathway modeling.     

Hot Spots of Nitrification in the Elbe Estuary and Their Impact on Nitrate Regeneration

Estuaries act as an organic matter and nutrient filter in the transition between the land, rivers and the ocean. In the past, high nutrient and organic carbon load and low oxygen concentration made the Elbe River estuary (NW Europe) a sink for dissolved inorganic nitrogen. A recent reduction in loads and subsequent recovery of the estuary changed its biogeochemical function, so that nitrate is no longer removed on its transition towards the coastal North Sea. Nowadays in the estuary, nitrification appears to be a significant nitrate source. To quantify nitrification and determine actively nitrifying regions in the estuary, we measured the concentrations of ammonium, nitrite and nitrate, the dual stable isotopes of nitrate and net nitrification rates in the estuary on five cruises from August 2012 to August 2013. The nitrate concentration increased markedly downstream of the port of Hamburg in summer and spring, accompanied by a decrease of nitrate isotope values that was clearest in summer exactly at the location where nitrate concentration started to increase. Ammonium and nitrite peaked in the Hamburg port region (up to 18 and 8 μmol L^?1, respectively), and nitrification rates in this region were up to 7 μmol L^?1 day^?1. Our data show that coupled re-mineralization and nitrification are significant internal nitrate sources that almost double the estuary’s summer nitrate concentration. Furthermore, we find that the port of Hamburg is a hot spot of nitrification, whereas the maximum turbidity zone (MTZ) only plays a subordinate role in turnover of nitrate.     

The Younger Dryas in the English Lake District: reconciling geomorphological evidence with numerical model outputs

The geomorphology of the south‐western and central Lake District, England is used to reconstruct the mountain palaeoglaciology pertaining to the Lateglacial and Younger Dryas. Limitations to previous ice‐mass reconstructions and consequent palaeoclimatic inferences include: (i) the use of static (steady‐state) glacier reconstructions, (ii) the assumption of a single‐stage Younger Dryas advance, (iii) greatly varying ice‐volume estimates, (iv) inexplicable spatial variations in ELA (Equilibrium Line Altitude), and (v) a lack of robust extent chronology. Here we present geomorphological mapping based on aerial photography and the NextMap Britain Digital Elevation Model, checked by ground survey. Former glacier extents were inferred and ELAs were calculated using the Balance Ratio method of Osmaston. Independently, a time‐dependant 2‐D ice‐flow model was forced by a regional ELA history that was scaled to the GRIP record. This provided a dynamic reconstruction of a mountain ice field that allowed for non‐steady‐state glacier evolution. Fluctuations in climate during the Younger Dryas resulted in multiple glacial advance positions that show agreement with the location of mapped moraines, and may further explain some of the ELA variations found in previous local and static reconstructions. Modelling based on the GRIP record predicts three phases: an initial maximum extent, a middle minor advance or stillstand, and a pronounced but less extensive final advance. The comparisons find that the reconstructions derived from geomorphological evidence are effective representations of steady‐state glacier geometries, but we do propose different extents for some glaciers and, in particular, a large former glacier in Upper Eskdale.