TESSIN VISLab—laboratory for scientific visualization

Scientific visualization is an integral part of the modeling workflow, enabling researchers to understand complex or large data sets and simulation results. A high-resolution stereoscopic virtual reality (VR) environment further enhances the possibilities of visualization. Such an environment also allows collaboration in work groups including people of different backgrounds and to present results of research projects to stakeholders or the public. The requirements for the computing equipment driving the VR environment demand specialized software applications which can be run in a parallel fashion on a set of interconnected machines. Another challenge is to devise a useful data workflow from source data sets onto the display system. Therefore, we develop software applications like the OpenGeoSys Data Explorer, custom data conversion tools for established visualization packages such as ParaView and Visualization Toolkit as well as presentation and interaction techniques for 3D applications like Unity. We demonstrate our workflow by presenting visualization results for case studies from a broad range of applications. An outlook on how visualization techniques can be deeply integrated into the simulation process is given and future technical improvements such as a simplified hardware setup are outlined.     

Characterisation of contemporary local climate change in the mountains of southwest Bulgaria

The Pirin Mountains in southwest Bulgaria spatially mark a transition between the Mediterranean and temperate climate zones. Therefore they are also particularly relevant for research on high mountain climate and the effect of landscape transformation. Historical climate records gathered in the area have been researched, checked and statistically examined. The mountainous climate has been characterised and trends in the evolution of temperature and precipitation since 1931 have been outlined. There are objective evidences for an increasing annual mean temperature, longer vegetative periods and local droughts in spring and autumn. Significant changes also appear in climatic threshold values such as the number of frost change days. This last parameter is very important for the sustainability of mountainous ecosystems.     

Recent behavior and possible future evolution of the glacieret in the cirque Golemiya Kazan in the Pirin Mountains under conditions of climate warming

This research reveals relationships between climate variables and inter-annual dynamics in the area of the glacieret located in the cirque Golemiya Kazan in the Pirin Mountains. The study period is 1993–2017. The correlations are identified using statistical methods. Also, a statistical model is constructed, including some climate variables as predictors. Despite the evident decrease of the glacieret's size in the period from the 1950 s onwards, the long-term trends for the last decades have been insignificant. The main climatic factors influencing the inter-annual dynamics in the area of the glacieret are air temperature, precipitation, zonal and meridional winds and relative humidity. With respect to the dynamics in the area of the glacieret, the important trends in the different climate variables are those of the warm period air temperatures and zonal(u) wind. They also determine to a great extent its future development by acting in two opposite directions – rising temperatures in the warm period will lead to a rapid decrease of its area by the end of the melting season, while the change of wind directionfrom west to east in the warm period will increase its area. The influence of the zonal wind in the warm period is explained mainly by the location of the glacieret in the cirque. Generally, the glacieret is tilted downwards from west to east. Thus, westerly winds facilitate blowing away the snow from the surface of the glacieret, assisting its melting in the warm period. Easterly winds do not have such an effect. The combination of the opposite effects of these two most important climate variables leads to the most likely scenario for the future development of the glacieret, according to which by the middle of this century it is expected to turn into a semi-permanent snow patch, which disappears after some summers, and by the end of the century to completely melt every year before the end of the melting season.     

New insights into the origin and migration of brines in deep Devonian aquifers, Alberta, Canada

Analysis of hydraulic heads and chemical compositions of Devonian formation waters in the west central part of the Alberta Basin, Canada, characterizes the origin of formation waters and migration of brines. The Devonian succession in the study area lies 2000–5000 m below the ground surface, and has an approximate total thickness of 1000 m and an average slope of 15 m/km. Four Devonian aquifers are present in the study area, which form two aquifer systems [i.e., a Middle–Upper Devonian aquifer system (MUDAS) consisting of the Elk Point and Woodbend–Beaverhill Lake aquifers, and an Upper Devonian aquifer system (UDAS) consisting of the Winterburn and Wabamun aquifers]. The Ireton is an effective aquitard between these two systems in the eastern parts of the study area. The entire Devonian succession is confined below by efficient aquitards of the underlying Cambrian shales and/or the Precambrian basement, and above by overlying Carboniferous shales of the Exshaw and Lower Banff Formations.The formation water chemistry shows that the Devonian succession contains two distinct brine types: a ‘heavy brine,’ located updip, defined approximately by TDS >200 g/l, and a ‘light brine’ with TDS <200 g/l. Hydraulic head distributions suggest that, presently, the ‘light brine’ attempts to flow updip, thereby pushing the ‘heavy brine’ ahead. The interface between the two brines is lobate and forms large-scale tongues that are due to channeled flow along high-permeability pathways. Geological and hydrogeochemical data suggest that the following processes determined the present composition of the ‘light’ and ‘heavy’ brines: original seawater, evaporation beyond gypsum but below halite saturation, dolomitization, clay dehydration, gypsum dewatering, thermochemical sulfate reduction (TSR), and halite dissolution. The influx of meteoric (from the south) and metamorphic (from the west) waters can be recognized only in the ‘light brine.’ Albitization can be unequivocally identified only in the ‘heavy brine.’ The ‘heavy brine’ may be residual Middle Devonian evaporitic brine from the Williston Basin or the Elk Point Basin, or it may have originated from partial dissolution of thick, laterally extensive Middle Devonian evaporite deposits to the east of the study area. The ‘light brine’ most probably originated from dilution of ‘heavy brine’ in post-Laramide times.     

Primitive andesites from the Taupo Volcanic Zone formed by magma mixing

Andesites with Mg# >45 erupted at subduction zones form either by partial melting of metasomatized mantle or by mixing and assimilation processes during melt ascent. Primitive whole rock basaltic andesites from the Pukeonake vent in the Tongariro Volcanic Centre in New Zealand’s Taupo Volcanic Zone contain olivine, clino- and orthopyroxene, and plagioclase xeno- and antecrysts in a partly glassy matrix. Glass pools interstitial between minerals and glass inclusions in clinopyroxene, orthopyroxene and plagioclase as well as matrix glasses are rhyolitic to dacitic indicating that the melts were more evolved than their andesitic bulk host rock analyses indicate. Olivine xenocrysts have high Fo contents up to 94%, δ¹⁸O_(SMOW) of +5.1‰, and contain Cr-spinel inclusions, all of which imply an origin in equilibrium with primitive mantle-derived melts. Mineral zoning in olivine, clinopyroxene and plagioclase suggest that fractional crystallization occurred. Elevated O isotope ratios in clinopyroxene and glass indicate that the lavas assimilated sedimentary rocks during stagnation in the crust. Thus, the Pukeonake andesites formed by a combination of fractional crystallization, assimilation of crustal rocks, and mixing of dacite liquid with mantle-derived minerals in a complex crustal magma system. The disequilibrium textures and O isotope compositions of the minerals indicate mixing processes on timescales of less than a year prior to eruption. Similar processes may occur in other subduction zones and require careful study of the lavas to determine the origin of andesite magmas in arc volcanoes situated on continental crust.     

The Bingham Distribution of Quaternions and Its Spherical Radon Transform in Texture Analysis

Spherical geometry of quaternions is employed to characterize the Bingham distribution on the 3-dimensional sphere as being uniquely composed of a bipolar, a circular and a spherical component. A new parametrization of its dispersion parameters provides a classification of patterns of crystallographic preferred orientations (CPO, or textures). It is shown that the Bingham distribution can represent most types of ideal CPO patterns; in particular single component, fiber and surface textures are represented by rotationally invariant bipolar, circular and spherical distributions, respectively. Pole figures of given crystal directions are derived by the spherical Radon transform of the Bingham probability density function of rotations, which are displayed for general and special cases.     

Bottom-current control on sedimentation in the western Bellingshausen Sea, West Antarctica

A set of single- and multi-channel seismic reflection profiles provide insights into the younger Cenozoic sedimentation history of the continental rise in the western Bellingshausen Sea, west and north of Peter I Island. This area has been strongly influenced by glacially controlled sediment supply from the continental shelf, interacting with a westward-flowing bottom current. From south to north, the seismic data show changes in the symmetry and structure of a prominent sediment depocentre. Its southernmost sector provides evidence of sediment drift whereas northwards the data show a large channel-levee complex, with a western levee oriented in the opposite direction to that of the drift in the south. This pattern indicates the northward-decreasing influence of a westward-flowing bottom contour current in the study area. Topographic data suggest the morphologic ridges at Peter I Island to be the main features responsible for variable bottom-current influence, these acting as barrier to the bottom current and entrained sedimentary material. West of Peter I Island, the east-orientated Coriolis force remains effective in deflecting the suspended load of the turbidity currents towards the west, thereby promoting growth of the western channel levee. Calculated sediment accumulation rates based on seismic data reveal Depocentre C to consist of younger Cenozoic material supplied by glacial transport and modified by contour currents in the western Bellingshausen Sea. These findings demonstrate that the shape, structure and distribution of sediment mounds and estimates of sediment accumulation rates can be associated to the influence of bottom currents and their long-term evolution in response to tectonic movements, ice-sheet dynamics and deep-water formation.     

An easy approach to assess the susceptibility of a landslide by utilizing simple raster overlay methods: A case study on Huangtupo landslide (P.R. China)

Semi qualitative index based methods using rankings and ratings are commonly used in susceptibility estimations over a wide area. However, generalized ranking and ratings are not applicable for one single landslide. This paper gives an easy and transferable approach to a susceptibility assessment of Huangtupo landslide (P.R. China), using raster addition without taking account for ranking and ratings. Slope, aspect, curvature, location and drainage buffer distance raster data sets have been obtained out of open source digital elevation models using ESRI’s ArcGIS. These conditioning factor raster data sets have been translated into raster data sets including simple yes or no criteria, referring to triggering or not. Subsequently they have been added by raster math to acquire a simple raster overlay map. After that this map is compared to initial displacement measurements, obtained by using a ground based synthetic aperture radar device. Acquired data is recalculated to a raster data set using the same spatial extent, to provide the possibility of comparison of the two raster data sets. The results reveal, that 76.35% of all measured movements occur in areas where raster cells include three or more conditioning factors, indicating that easy raster math operations can lead to satisfying results in local scale.     

Impact of multi-purpose aquifer utilisation on a variable-density groundwater flow system in the Gippsland Basin, Australia

The Latrobe aquifer in the Gippsland Basin in southeastern Australia is a prime example for emerging resource conflicts in Australian sedimentary basins. The Latrobe Group forms a major freshwater aquifer in the onshore Gippsland Basin, and is an important reservoir for oil and gas in both onshore and offshore parts of the basin. The Latrobe Group and overlying formations contain substantial coal resources that are being mined in the onshore part of the basin. These may have coal-seam-gas potential and, in addition, the basin is considered prospective for its geothermal energy and CO₂ storage potential. The impacts of groundwater extraction related to coal-mine dewatering, public water supply, and petroleum production on the flow of variable-density formation water has been assessed using freshwater hydraulic heads and impelling force vectors. Groundwater flows from the northern and western edges towards the central part of the basin. Groundwater discharge occurs mainly offshore along the southern margin. Post-stress hydraulic heads show significant declines near the petroleum fields and in the coal mining areas. A hydrodynamic model of the Latrobe aquifer was used to simulate groundwater recovery in the Latrobe aquifer from different scenarios of cessation of groundwater and other fluid extractions.     

Concept and workflow for 3D visualization of atmospheric data in a virtual reality environment for analytical approaches

In the future, climate change will strongly influence our environment and living conditions. Weather and Climate simulations that predict possible changes produce big data sets. The combination of various variables of climate models with spatial data from different sources helps to identify correlations and to study key processes. In this paper, the results of the Weather Research and Forecasting model are visualized for two regions. For this purpose, a continuous workflow that leads from the integration of heterogeneous raw data to 3D visualizations that can be displayed on a desktop computer or in an interactive virtual reality environment is developed. These easy-to-understand visualizations of complex data are the basis for scientific communication and for the evaluation and verification of models as well as for interdisciplinary discussions of the research results.