DYPTOP: a cost-efficient TOPMODEL implementation to simulate sub-grid spatio-temporal dynamics of global wetlands and peatlands

Simulating the spatio-temporal dynamics of inundation is key to understanding the role of wetlands under past and future climate change. Earlier modelling studies have mostly relied on fixed prescribed peatland maps and inundation time series of limited temporal coverage. Here, we describe and assess the the Dynamical Peatland Model Based on TOPMODEL (DYPTOP), which predicts the extent of inundation based on a computationally efficient TOPMODEL implementation. This approach rests on an empirical, grid-cell-specific relationship between the mean soil water balance and the flooded area. DYPTOP combines the simulated inundation extent and its temporal persistency with criteria for the ecosystem water balance and the modelled peatland-specific soil carbon balance to predict the global distribution of peatlands. We apply DYPTOP in combination with the LPX-Bern DGVM and benchmark the global-scale distribution, extent, and seasonality of inundation against satellite data. DYPTOP successfully predicts the spatial distribution and extent of wetlands and major boreal and tropical peatland complexes and reveals the governing limitations to peatland occurrence across the globe. Peatlands covering large boreal lowlands are reproduced only when accounting for a positive feedback induced by the enhanced mean soil water holding capacity in peatland-dominated regions. DYPTOP is designed to minimize input data requirements, optimizes computational efficiency and allows for a modular adoption in Earth system models.

A community-based geological reconstruction of Antarctic Ice Sheet deglaciation since the Last Glacial Maximum

A robust understanding of Antarctic Ice Sheet deglacial history since the Last Glacial Maximum is important in order to constrain ice sheet and glacial-isostatic adjustment models, and to explore the forcing mechanisms responsible for ice sheet retreat. Such understanding can be derived from a broad range of geological and glaciological datasets and recent decades have seen an upsurge in such data gathering around the continent and Sub-Antarctic islands. Here, we report a new synthesis of those datasets, based on an accompanying series of reviews of the geological data, organised by sector. We present a series of timeslice maps for 20 ka, 15 ka, 10 ka and 5 ka, including grounding line position and ice sheet thickness changes, along with a clear assessment of levels of confidence. The reconstruction shows that the Antarctic Ice sheet did not everywhere reach the continental shelf edge at its maximum, that initial retreat was asynchronous, and that the spatial pattern of deglaciation was highly variable, particularly on the inner shelf. The deglacial reconstruction is consistent with a moderate overall excess ice volume and with a relatively small Antarctic contribution to meltwater pulse 1a. We discuss key areas of uncertainty both around the continent and by time interval, and we highlight potential priorities for future work. The synthesis is intended to be a resource for the modelling and glacial geological community.

Using 18O/2H, 3H/3He, 85Kr and CFCs to determine mean residence times and water origin in the Grazer and Leibnitzer Feld groundwater bodies (Austria)

Two groundwater bodies, Grazer Feld and Leibnitzer Feld, with surface areas of 166 and 103 km2 respectively are characterised for the first time by measuring the combination of d18O/d2H, 3H/3He, 85Kr, CFC-11, CFC-12 and hydrochemistry in 34 monitoring wells in 2009/2010. The timescales of groundwater recharge have been characterised by 131 d18O measurements of well and surface water sampled on a seasonal basis. Most monitoring wells show a seasonal variation or indicate variable contributions of the main river Mur (0–30%, max. 70%) and/or other rivers having their recharge areas in higher altitudes. Combined d18O/d2H-measurements indicate that 65–75% of groundwater recharge in the unusual wet year of 2009 was from precipitation in the summer based on values from the Graz meteorological station. Monitoring wells downstream of gravel pit lakes show a clear evaporation trend. A boron–nitrate differentiation plot shows more frequent boron-rich water in the more urbanised Grazer Feld and more frequent nitrate-rich water in the more agricultural used Leibnitzer Feld indicating that a some of the nitrate load in the Grazer Feld comes from urban sewer water. Several lumped parameter models based on tritium input data from Graz and monthly data from the river Mur (Spielfeld) since 1977 yield a Mean Residence Time (MRT) for the Mur-water itself between 3 and 4 years in this area. Data from d18O, 3H/3He measurements at the Wagna lysimeter station supports the conclusion that 90% of the groundwaters in the Grazer Feld and 73% in the Leibnitzer Feld have MRTs of <5 years. Only in a few groundwaters were MRTs of 6–10 or 11–25 years as a result of either a long-distance water inflow in the basins or due to longer flow path in somewhat deeper wells (>20 m) with relative thicker unsaturated zones. The young MRT of groundwater from two monitoring wells in the Leibnitzer Feld was confirmed by 85Kr-measurements. Most CFC-11 and CFC-12 concentrations in the groundwater exceed the equilibration concentrations of modern concentrations in water and are therefore unsuitable for dating purposes. An enrichment factor up to 100 compared to atmospheric equilibrium concentrations and the obvious correlation of CFC-12 with SO4, Na, Cl and B in the ground waters of the Grazer Feld suggest that waste water in contact with CFC-containing material above and below ground is the source for the contamination. The dominance of very young groundwater (<5 years) indicates a recent origin of the contamination by nitrate and many other components observed in parts of the groundwater bodies. Rapid measures to reduce those sources are needed to mitigate against further deterioration of these waters.

Age of cleft monazites in the eastern Tauern Window: constraints on crystallization conditions of hydrothermal monazite

Monazite-bearing Alpine clefts located in the Sonnblick region of the eastern Tauern Window, Austria, are oriented perpendicular to the foliation and lineation. Ion probe (SIMS) Th–Pb and U–Pb dating of four cleft monazites yields crystallization ages of different growth domains and aggregate regions ranging from 18.99 ± 0.51 to 15.00 ± 0.51 Ma. The crystallization ages obtained are overlapping or slightly younger than zircon fission track ages but older than zircon (U–Th)/He cooling ages from the same area. This constrains cleft monazite crystallization in this area to *300–200 �C. LA-ICP-MS data of dated hydrothermal monazites indicate that in graphite-bearing, reduced host lithologies, cleft monazite is poor in As and has higher La/Yb values and U concentrations, whereas in oxidised host rocks opposite trends are observed. Monazites show negative Eu anomalies and variable La/Yb values ranging from 520 to 6050. The positive correlation between Ca and Sr concentration indicates dissolution of plagioclase or carbonates as the source of these elements. The data show that early exhumation and cleft formation in the Tauern is related to metamorphic dome formation caused by the collision of the Adriatic with the European plate and that monazite crystallization in the clefts occurred later. Our data also demonstrate that hydrothermal monazite ages offer great potential in helping to constrain the chronology of exhumation in collisional orogens.

XMapTools: a MATLAB©-based program for electron microprobe X-ray image processing and geothermobarometry

XMapTools is a MATLAB©-based graphical user interface program for electron microprobe X-ray image processing, which can be used to estimate the pressure–temperature conditions of crystallization of minerals in metamorphic rocks. This program (available online at http://www.xmaptools.com) provides a method to standardize raw electron microprobe data and includes functions to calculate the oxide weight percent compositions for various minerals. A set of external functions is provided to calculate structural formulae from the standardized analyses as well as to estimate pressure–temperature conditions of crystallization, using empirical and semi-empirical thermobarometers from the literature. Two graphical user interface modules, Chem2D and Triplot3D, are used to plot mineral compositions into binary and ternary diagrams. As an example, the software is used to study a high-pressure Himalayan eclogite sample from the Stak massif in Pakistan. The high-pressure paragenesis consisting of omphacite and garnet has been retrogressed to a symplectitic assemblage of amphibole, plagioclase and clinopyroxene. Mineral compositions corresponding to ~165,000 analyses yield estimates for the eclogitic pressure–temperature retrograde path from 25 kbar to 9 kbar. Corresponding pressure–temperature maps were plotted and used to interpret the link between the equilibrium conditions of crystallization and the symplectitic microstructures. This example illustrates the usefulness of XMapTools for studying variations of the chemical composition of minerals and for retrieving information on metamorphic conditions on a microscale, towards computation of continuous pressure–temperature-and relative time path in zoned metamorphic minerals not affected by post-crystallization diffusion.

Microscale mapping of alteration conditions and potential biosignatures in basaltic-ultramafic rocks on early Earth and beyond

Subseafloor environments preserved in Archean greenstone belts provide an analogue for investigating potential subsurface habitats on Mars. The c. 3.5-3.4 Ga pillow lava metabasalts of the mid-Archean Barberton greenstone belt, South Africa, have been argued to contain the earliest evidence for microbial subseafloor life. This includes candidate trace fossils in the form of titanite microtextures, and sulfur isotopic signatures of pyrite preserved in metabasaltic glass of the c. 3.472 Ga Hooggenoeg Formation. It has been contended that similar microtextures in altered martian basalts may represent potential extraterrestrial biosignatures of microbe-fluid-rock interaction. But despite numerous studies describing these putative early traces of life, a detailed metamorphic characterization of the microtextures and their host alteration conditions in the ancient pillow lava metabasites is lacking. Here, we present a new nondestructive technique with which to study the in situ metamorphic alteration conditions associated with potential biosignatures in mafic-ultramafic rocks of the Hooggenoeg Formation. Our approach combines quantitative microscale compositional mapping by electron microprobe with inverse thermodynamic modeling to derive low-temperature chlorite crystallization conditions. We found that the titanite microtextures formed under subgreenschist to greenschist facies conditions. Two chlorite temperature groups were identified in the maps surrounding the titanite microtextures and record peak metamorphic conditions at 315 ± 40°C (XFe3+(chlorite) = 25-34%) and lower-temperature chlorite veins/microdomains at T = 210 ± 40°C (lower XFe3+(chlorite) = 40-45%). These results provide the first metamorphic constraints in textural context on the Barberton titanite microtextures and thereby improve our understanding of the local preservation conditions of these potential biosignatures. We suggest that this approach may prove to be an important tool in future studies to assess the biogenicity of these earliest candidate traces of life on Earth. Furthermore, we propose that this mapping approach could also be used to investigate altered mafic-ultramafic extraterrestrial samples containing candidate biosignatures.

Verbalization of Dependencies Between Concepts Built Through Fuzzy Cognitive Maps

The new computing paradigm known as cognitive computing attempts to imitate the human capabilities of learning, problem solving, and considering things in context. To do so, an application (a cognitive system) must learn from its environment (e.g., by interacting with various interfaces). These interfaces can run the gamut from sensors to humans to databases. Accessing data through such interfaces allows the system to conduct cognitive tasks that can support humans in decision-making or problem-solving processes. Cognitive systems can be integrated into various domains (e.g., medicine or insurance). For example, a cognitive system in cities can collect data, can learn from various data sources and can then attempt to connect these sources to provide real time optimizations of subsystems within the city (e.g., the transportation system). In this study, we provide a methodology for integrating a cognitive system that allows data to be verbalized, making the causalities and hypotheses generated from the cognitive system more understandable to humans. We abstract a city subsystem—passenger flow for a taxi company—by applying fuzzy cognitive maps (FCMs). FCMs can be used as a mathematical tool for modeling complex systems built by directed graphs with concepts (e.g., policies, events, and/or domains) as nodes and causalities as edges. As a verbalization technique we introduce the restriction-centered theory of reasoning (RCT). RCT addresses the imprecision inherent in language by introducing restrictions. Using this underlying combinatorial design, our approach can handle large data sets from complex systems and make the output understandable to humans.

Geological map of the Scalaro Valley – Sesia Zone (Italian Western Alps)

In the Sesia Zone (Italian Western Alps), slivers of continental crust characterised by an Alpine high-pressure imprint are intermingled with abundant mafic rocks and Mesozoic metasediments. An extensive study of the central Sesia Zone was undertaken to identify and reconstruct the lithological setting of the mono-cyclic sediments of the Scalaro Unit. A new geological map (1:5000) and schematic cross sections across the Scalaro Unit and the adjoining Eclogitic Micaschist Complex are presented here. In order to delimit the size and shape of the mono-metamorphic unit and understand its internal geometry with respect to the poly-metamorphic basement, an integrated approach was used. Linking observations and data across a range of scales, from kilometres in the field down to petrological and chronological data obtained at micrometre scale, we define for the first time the real size and internal geometry of the Scalaro Unit, as well as its large-scale structural context.

Monitoring the temporal development of natural hazard risks as a basis indicator for climate change adaptation.

The potential effects of climatic changes on natural risks are widely discussed. But the formulation of strategies for adapting risk management practice to climate changes requires knowledge of the related risks for people and economic values. The main goals of this work were (1) the development of a method for analysing and comparing risks induced by different natural hazard types, (2) highlighting the most relevant natural hazard processes and related damages, (3) the development of an information system for the monitoring of the temporal development of natural hazard risk and (4) the visualisation of the resulting information for the wider public. A comparative exposure analysis provides the basis for pointing out the hot spots of natural hazard risks in the province of Carinthia, Austria. An analysis of flood risks in all municipalities provides the basis for setting the priorities in the planning of flood protection measures. The methods form the basis for a monitoring system that periodically observes the temporal development of natural hazard risks. This makes it possible firstly to identify situations in which natural hazard risks are rising and secondly to differentiate between the most relevant factors responsible for the increasing risks. The factors that most influence the natural risks could be made evident.