The largest karst cave discovered in a tunnel during motorway construction in Slovenia’s Classical Karst (Kras)

Karst researchers of the Karst Research Institute ZRC SAZU are regularly consulted during the planning of Slovenia’s motorways and invited to observe and monitor construction in the karst areas. More than 350 caves have been discovered in the course of building 60 km of new motorways over the past decade. Access to the most important caves is preserved by concrete tubes closed with metal covers at the roadside. The largest cave system in a tunnel, named LC-S647, is almost entirely preserved. The protection of this cave will serve as a good example for the preservation of natural heritage in the future. Construction work for the motorway uncovered a high degree of karstification. In particular, the discovery of this cave in the Kastelec tunnel LC-S647 showed the existence of a major cave system in the geological and speleological past. A number of unconnected passages have also been found. Caves discovered during highway construction have brought new knowledge about the cavernosity and the geological history of this part of the karst.     

Soil contamination from urban and industrial activity: example of the mining district of Linares (southern Spain)

The Linares region (southern Spain) has been subjected to two important sources of pollution: the intensive mining works and the urban-industrial activity. To obtain a geochemical characterisation of the soil, 31 trace elements were analysed and 669 soil samples were collected. By means of clustering analysis, we identified groups of elements and grid squares in which relations could be established concerning soil lithology, urban and industrial activities and the degree of pollution impact; in addition, we were able to characterise the geochemical background of the study area. The multivariate study led us to identify four factors. Particularly important was factor 2, which represented the elements associated with mineral paragenesis (Cu, Pb, As, Co, Mn, Zn, Sn, Ba). This factor also contains elements related with an urban-industrial activity, such as Pb, Cu, Zn, As and Ba. Furthermore, we identified factor 4, associating Ni, V and Cr, and which is related to the use of fuels.     

Landslide susceptibility mapping for a landslide-prone area (Findikli,NE of Turkey) by likelihood-frequency ratio and weighted linear combination models

Landslides are very common natural problems in the Black Sea Region of Turkey due to the steep topography, improper use of land cover and adverse climatic conditions for landslides. In the western part of region, many studies have been carried out especially in the last decade for landslide susceptibility mapping using different evaluation methods such as deterministic approach, landslide distribution, qualitative, statistical and distribution-free analyses. The purpose of this study is to produce landslide susceptibility maps of a landslide-prone area (Findikli district, Rize) located at the eastern part of the Black Sea Region of Turkey by likelihood frequency ratio (LRM) model and weighted linear combination (WLC) model and to compare the results obtained. For this purpose, landslide inventory map of the area were prepared for the years of 1983 and 1995 by detailed field surveys and aerial-photography studies. Slope angle, slope aspect, lithology, distance from drainage lines, distance from roads and the land-cover of the study area are considered as the landslide-conditioning parameters. The differences between the susceptibility maps derived by the LRM and the WLC models are relatively minor when broad-based classifications are taken into account. However, the WLC map showed more details but the other map produced by LRM model produced weak results. The reason for this result is considered to be the fact that the majority of pixels in the LRM map have high values than the WLC-derived susceptibility map. In order to validate the two susceptibility maps, both of them were compared with the landslide inventory map. Although the landslides do not exist in the very high susceptibility class of the both maps, 79% of the landslides fall into the high and very high susceptibility zones of the WLC map while this is 49% for the LRM map. This shows that the WLC model exhibited higher performance than the LRM model.     

Development of a 3D GIS and its application to karst areas

There is a growing interest in modeling and analyzing karst phenomena in three dimensions. This paper integrates geology, groundwater hydrology, geographic information system (GIS), database management system (DBMS), visualization and data mining to study karst features in Huaibei, China. The 3D geo-objects retrieved from the karst area are analyzed and mapped into different abstract levels. The spatial relationships among the objects are constructed by a dual-linker. The shapes of the 3D objects and the topological models with attributes are stored and maintained in the DBMS. Spatial analysis was then used to integrate the data in the DBMS and the 3D model to form a virtual reality (VR) to provide analytical functions such as distribution analysis, correlation query, and probability assessment. The research successfully implements 3D modeling and analyses in the karst area, and meanwhile provides an efficient tool for government policy-makers to set out restrictions on water resource development in the area.     

Multiple discriminant analysis of the Drava River alluvial plain sediments

Three discriminant function models are raised and cross-compared in order to distinguish geochemical patterns characteristic for the Drava River floodplain sediments. Based on data representing total element concentrations in samples collected from alluvium (A), terrace (T), and unconsolidated bedrock (B) at the border of a floodplain, four element clusters emerged accounting for discrimination between the referred groups of sediments. The most prominent is contaminant/carbonate cluster characteristic for alluvium. The other two are: silicate cluster typical for unconsolidated geological substrate (Neogene sedimentary rocks); and naturally dispersed heavy metal cluster separating terrace from the former two groups. Models introducing depth intervals and single profiles as grouping criteria reveal identical sediment-heavy metal matrices. The second important issue of this paper is possibility of reclassification of samples originally assigned to one of the a priori defined groups of sediments, based on established geochemical pattern. The mapped geological units can be reconsidered by the post hoc assignments to a different group if geological border between alluvium and terrace or between terrace and bedrock can not be established geologically with absolute certainty.     

The double solid reactant method: II. An application to the shallow groundwaters of the Porto Plain,Vulcano Island (Italy)

This paper presents an example of application of the double solid reactant method (DSRM) of Accornero and Marini (Environmental Geology, 2007a), an effective way for modeling the fate of several dissolved trace elements during water–rock interaction. The EQ3/6 software package was used for simulating the irreversible water–rock mass transfer accompanying the generation of the groundwaters of the Porto Plain shallow aquifer, starting from a degassed diluted crateric steam condensate. Reaction path modeling was performed in reaction progress mode and under closed-system conditions. The simulations assumed: (1) bulk dissolution (i.e., without any constraint on the kinetics of dissolution/precipitation reactions) of a single solid phase, a leucite-latitic glass, and (2) precipitation of amorphous silica, barite, alunite, jarosite, anhydrite, kaolinite, a solid mixture of smectites, fluorite, a solid mixture of hydroxides, illite-K, a solid mixture of saponites, a solid mixture of trigonal carbonates and a solid mixture of orthorhombic carbonates. Analytical concentrations of major chemical elements and several trace elements (Cr, Mn, Fe, Ni, Cu, Zn, As, Sr and Ba) in groundwaters were satisfactorily reproduced. In addition to these simulations, similar runs for a rhyolite, a latite and a trachyte permitted to calculate major oxide contents for the authigenic paragenesis which are comparable, to a first approximation, with the corresponding data measured for local altered rocks belonging to the silicic, advanced argillic and intermediate argillic alteration facies. The important role played by both the solid mixture of trigonal carbonates as sequestrator of Mn, Zn, Cu and Ni and the solid mixture of orthorhombic carbonates as scavenger of Sr and Ba is emphasized.

Quantitative textural analysis of packings of elongate crystals

The spatial distribution of grains in a solidifying igneous rock controls the physical properties of the crystal mush, and in turn is controlled by the rate of crystal growth and accumulation. A predominant non-spherical habit for igneous minerals brings into question the use of spherical particles in reference packings used for quantification of spatial distribution. Furthermore, variations of crystal clustering/ordering with length scale require spatial statistics which take into account the distribution of particles beyond nearest neighbours. Using random close packings of spherocylinders, we demonstrate the importance of aspect ratio for the aggregation index (usually known as R) and show that packings of spherical particles have more structure than packings of rods. The spatial distribution functions demonstrate that the plagioclase grains in the colonnade from the Holyoke basalt are clustered on a length scale of 0.5 mm. Understanding the controls on grain spatial distribution in igneous rocks will depend on the application of these techniques to well-understood environments.     

Seismic Hazard Analysis — Quo vadis?

The paper is dedicated to the review of methods of seismic hazard analysis currently in use, analyzing the strengths and weaknesses of different approaches. The review is performed from the perspective of a user of the results of seismic hazard analysis for different applications such as the design of critical and general (non-critical) civil infrastructures, technical and financial risk analysis. A set of criteria is developed for and applied to an objective assessment of the capabilities of different analysis methods. It is demonstrated that traditional probabilistic seismic hazard analysis (PSHA) methods have significant deficiencies, thus limiting their practical applications. These deficiencies have their roots in the use of inadequate probabilistic models and insufficient understanding of modern concepts of risk analysis, as have been revealed in some recent large scale studies. These deficiencies result in the lack of ability of a correct treatment of dependencies between physical parameters and finally, in an incorrect treatment of uncertainties. As a consequence, results of PSHA studies have been found to be unrealistic in comparison with empirical information from the real world. The attempt to compensate these problems by a systematic use of expert elicitation has, so far, not resulted in any improvement of the situation. It is also shown that scenario-earthquakes developed by disaggregation from the results of a traditional PSHA may not be conservative with respect to energy conservation and should not be used for the design of critical infrastructures without validation. Because the assessment of technical as well as of financial risks associated with potential damages of earthquakes need a risk analysis, current method is based on a probabilistic approach with its unsolved deficiencies.

Traditional deterministic or scenario-based seismic hazard analysis methods provide a reliable and in general robust design basis for applications such as the design of critical infrastructures, especially with systematic sensitivity analyses based on validated phenomenological models. Deterministic seismic hazard analysis incorporates uncertainties in the safety factors. These factors are derived from experience as well as from expert judgment. Deterministic methods associated with high safety factors may lead to too conservative results, especially if applied for generally short-lived civil structures. Scenarios used in deterministic seismic hazard analysis have a clear physical basis. They are related to seismic sources discovered by geological, geomorphologic, geodetic and seismological investigations or derived from historical references. Scenario-based methods can be expanded for risk analysis applications with an extended data analysis providing the frequency of seismic events. Such an extension provides a better informed risk model that is suitable for risk-informed decision making.     

Framing in geographical analysis of environmental conflicts: Theory, methodology and three case studies

When applied to environmental concerns, framing offers a rigorous conceptual and analytic approach with potentially practical significance for dealing with complex issues relevant to geography. The purpose of this article is to introduce framing concepts, typology and modes of analysis to address issues which geographers commonly treat: disputes over land ownership and uses, competition for water resources, cultural clashes over control of territory and resources and the impact of spatial patterns and structures when siting noxious sites. Framing analysis is applied to three cases of environmental conflict in Israel and, more generally, for understanding conflicts revolving around the management of human habitats in relation to the physical environment.