Abstract

A geogenic CO2 emitting site (mofette U1) at the banks of the Laacher See, Eifel Mountains, was chosen to study the relationship between heavy postvolcanic soil degassing and vegetation during spring season. To test any interrelation between soil CO2 degassing and vegetation, soil chemism (pH, water content, conductivity, and humus content) and vegetation studies (number of species, plant-soil coverage) were performed. Geogenic soil degassing patterns of carbon dioxide and oxygen were clearly inhomogeneous, resembling soil porosity and distinct permeation channels within the soil. CO2 concentrations ranged from zero to 100%. Soil CO2 increased, while soil oxygen decreased with increasing soil depth. There was a reasonable correlation between CO2 degassing and soil pH as well as soil conductivity. Soil organic matter (SOM) resembled soil water distribution. The number of plant species (from a total of 69 species) as well as plant coverage strongly followed geogenic CO2 degassing. The total number of growing species was highest in low CO2 soils (max. 17 species per m2) and lowest at high CO2-emitting sites (one species per m2). Plant coverage followed the same pattern. Total plant coverage reached values of up to 84% in slightly degassing soils and only 5-6% on heavy CO2-venting sites. One plant species proved to be highly mofettophilic (marsh sedge, Carex acutiformis) and strictly grew on CO2 degassing sites. Most other species like grove windflower, spring fumewort, fig buttercup, wood bluegrass, addersmeat, and common snowberry showed a mofettophobic behavior and strictly avoided degassing areas. Specific plant species can thus be used to detect and monitor pre- or postvolcanic CO2 degassing.

1. Introduction

Mofettes are sites with dry CO2 gas exhalations at ambient temperatures. The gaseous CO2 originates from magma chambers or from Earth mantle degassing. It moves upward through fissures and cracks within the rocks [13]. Along its path, it reacts with wet and dry soil phases, finally exerting its influence on organisms living either within the subterranean soil or directly on the soil surface [410]. Depending on the local geological and hydrogeological conditions, the upward migration of CO2 gas may appear at the surface either as dry CO2 emanation (mofette) or as CO2-rich mineral water. Several publications already described the effects of the potentially acidic gas of dry CO2 emanations on the soil phase [1113], on soil insects, nematode worms, and spiders [1416], on soil microbiota [for review see [5]], and on vegetation. Most authors describe the effects of enhanced geogenic CO2 on physiological and ecophysiological processes, or anatomical and morphological adaptations of plants [68, 1722]. Only a few references can be found on the special vegetation of these sites—the mofette vegetation [17, 2326].

The present paper studies a natural, terrestrial mofette site at the Laacher See volcano, where CO2 from the Laacher See volcano still reaches the surface, influencing ecosystems. The quiet volcano is part of the east Eifel Volcanic field belonging to the West European Rift System. The last eruptive phase of the Laacher See volcano is about 13,000 years ago [27]. The aim of this work was to find out to which extent geogenic CO2 concentrations within the soil influence soil chemism and allow or inhibit plant life. It was furthermore of great interest to figure out whether specific plant species can be used as biological indicators to detect and monitor otherwise nondetected volcanic CO2 emissions.

2. Material and Methods

2.1. Site Description

The study area “U1 mofette” (Figure 1(a)) is located on the first terrace on the eastern inner flanks of the Laacher See Volcano, Eifel Mountains, Germany (see also [28]), just 5 m from the southern wall of the ruins of the so-called Jesuiten-Villa (a two-stored stone house used by the monks to celebrate their free days). It belongs to the cadastral district “Bei der rechten Erde” and extents parallel to the banks of the Laacher See. The terrace is the result of the second drainage phase of the lake performed in 1844 (first drainage in 1164; [29]). As the area is located within a forest, mostly half-light and only rarely full-light conditions prevail, with some heavily shaded parts at the upper left part of the area. Topographically, the U1 area is flat with the highest part in the SE corner (Figure 1(b)). The first 54 m of the 80-90 m long area was under investigation. In the lake close to the mofette area, CO2 bubbles indicate visible geogenic gas emission [30].

2.1.1. Sampling and Measurements

The area was chosen to study the relationship between geogenic soil degassing of carbon dioxide and the prevailing vegetation and between soil gas and selected geochemical and soil biological parameters. Seven transects with 2 m spacing and a total length of 54 m were laid out parallel to the lakeshore line. The hiking trail marked the western (lower) end of the study area (Figures 1(a) and 1(b)). All measurements were made at each 2 m intersection of the transect lines, resulting in 224 individual sample spots. Soil gas concentrations (CO2 and O2) were measured at all 224 spots at four different soil depths (10, 20, 40, and 60 cm), by coring holes into the soil and inserting gas tight tubing. During the measurements, the holes were closed with a lid to avoid air mixing during the actual measurements. The gas measurements were carried out with a portable landfill gas analyser GA2000 (Geotechnical Instruments, England) equipped with water and dust filters. Problems occurred when the ground was wet or waterlogged, as liquid water would damage the sensitive cells of the analyser. Polyethylene and Teflon tubing (1 cm diameter) was used to enlarge the measuring radius. Readings were recorded after several seconds to one minute.

Gas flux was performed at each grid intersection using a portable diffuse carbon dioxide flux meter system (West Systems Portable diffuse flux meter carbon dioxide high-flux, Pisa, Italy). The device consists of an accumulation chamber (type B), a CO2 IR detector (Polytron, Dräger), and a PDA palmtop (Brand, Acer n300) for data communication, evaluation, and storage.

The CO2 flux measurement method is based on the measurement of rising CO2 concentration versus time in terms of ppm s-1. If the gas concentration inside the chamber is constant, linear regression is used to calculate the gas flux [moles·m-2·day-1] by the following equation: where is the accumulation chamber factor and is the slope of the flux curve, determined by linear regression.

is calculated by the following equation: where is the barometric pressure (hPa), is the gas constant 0.08314510 (L·bar·K-1·mol-1), is the air temperature (K), is the chamber net volume (), and is the chamber inlet net area (). was measured with the GA2000 and with a thermocouple thermometer.

2.2. Soil Analysis

After determination of the soil gas, soil cores were taken very close to the points of the gas measurements. A standard soil borer (Pürkhauer; Thomas Müller, Germany) was used, and samples were taken between soil depths of 7 to 13 cm. This depth range was thought to reflect the main rooting horizon of herbaceous plants [31]; this soil depth is also the approximate range of the aerated surface soil [32]. The fresh samples were placed in small airtight plastic cylinders and stored cool until analysis.

Soil water content was determined by drying the soil samples of known fresh weight at 75°C for three days in an oven and reweighing after cooling to room temperature in a desiccator. Dry soil samples were ground, and soil pH was measured with a pH electrode (WTW, Germany) in bi-distilled water (750 mg dry soil in 15 ml water). After pH determination, soil conductivity was determined in the same solution using a conductivity meter (inoLab Multi 9420 IDS, WTW, Germany).

To determine soil organic matter (SOM), 1 g of ground, dry soil was heated (450°C for 18 hours) in an oven (muffle furnace; Heraeus, Germany), cooled to room temperature in a desiccator, and reweighed (for details see [33]).

2.3. Vegetation Analysis

Vegetation cover was estimated according to [34]. Londo’s decimal scale has some advantages over the most commonly used scale of [35]. The decimal scale is a pure dominance scale that uses smaller steps of maximal 10%. The more precise information about the plots’ vegetation cover leads to better correlations with the abiotic factors. On the other hand, a transformation into the Braun-Blanquet scale can be arranged without great difficulty [36].

All measurements were carried out at stabile weather conditions during a one-week period in February 2013. Vegetation analysis was done in May 2013.

2.4. Statistical Analysis

Statistical analysis was carried out using MS Excel (Microsoft, USA), SigmaPlot 11.0 (Systat Software, USA), and SPSS (SPSS, USA). Multivariable statistical analysis was done using multiattribute analysis based on canonical correspondence analysis (CCA; CANOCO 4.5 program). Probable correlations between plant habitat parameters, plant species composition, and plant distribution can thus be detected [37].

3. Results and Discussion

3.1. Degassing Patterns
3.1.1. Soil Gas Concentrations

Geogenic CO2 degassing was clearly heterogeneous within the selected area. At 10 cm soil depth, high [CO2] were found in a line running nearly diagonal from the left upper corner to the area close to the lower right part (Figure 2(a)). At this depth, CO2 extremes evolved as isolated islands surrounded by lower values.

Soil gas concentrations clearly changed with soil depth. Carbon dioxide increased with increasing soil depth whereas oxygen values decreased at deeper soil horizons (Figures 2(a)–(d) and 3(a)–(d)). Yet the principal “diagonal” degassing pattern was still visible although at a soil depth of 60 cm the whole area showed highly increased CO2 values. Thus, at deeper soil horizons, a clear broadening of the degassing pattern occurs which narrows when the gas moves upward to the soil surface. Sometimes, the narrowing of CO2 concentrations at lower soil depths is evident as small channels are formed showing the vertical gas-transducing chimneys. The transect illustration shows the modification of gas permeation pathways with soil depth.

There are two selected transects of the U1 mofette (Figure 4). The upper panel nicely shows two CO2 gas chimneys (at 12 m and 42 m) emerging at a soil depth of ca. 30 cm further narrowing to a diameter of less than 50 cm finishing at the very soil surface. Gas concentration at the surface thus reaches 80-100% CO2. Another CO2 channel (at 36 m) does not reach the actual soil surface but ends at 10 cm soil depth (Figure 4). A less permeable or impermeable soil layer transiently or permanently locks the permeating gas in pockets.

The same is true for oxygen, although in this case a narrowing of oxygen channels with increasing soil depths is visible (Figure 4(b)). Between 8 and 14 m and at point 48 m and 54 m, high O2 concentrations can be measured even at a soil depth of 40-50 cm. At 54 m, oxygen concentrations are high at even deeper soil levels (Figure 4). At these locations, control plants (mofettophobes) could easily grow.

Also for oxygen, a more or less diagonal concentration pattern from the left upper corner to the right lower corner of the area is seen (Figures 3(a)3(d)). Yet in contrast to the [CO2], an inverse pattern is found for [O2] (Figures 2 and 3). At spots with high [CO2], low [O2] is found and vice versa. If CO2 migrates towards the surface, mainly by diffusion and advective processes, it leads to a dilution of the in situ soil gas composition. Vodnik et al. [10] have already published the inverse proportionality between carbon dioxide and oxygen concentrations at mofette sites.

Oxygen were only apparent at soil depths below 40 cm. Nevertheless, even at a soil depth of 20 cm, the main rooting zone of many plants, [O2] below 15% (slightly hypoxic conditions) can be found in many spots.

3.1.2. CO2 Fluxes

Additional to the CO2 concentrations, soil CO2 fluxes within the U1 mofette were determined directly on the soil surface (Figure 5). Care was taken to avoid the inclusion of vegetation inside the accumulation chamber. As measurements were performed in February, also, biotic CO2 evolution by the soil edaphon is expected to be reduced due to the cooler winter temperatures. Postvolcanic soil CO2 fluxes nicely paralleled soil CO2 concentrations. Again, CO2 fluxes were clearly inhomogeneous within the area and a diagonal degassing pattern of enhanced CO2 flux is seen through the mofette area with the highest CO2 emission rates in the upper left corner of the area. Rates of 600 moles CO2 m-2 soil surface day-1 were seen in two heavily degassing islands (Figure 5). As with the [CO2], the upper right corner and the lower left part of the mofette still act as control sites (2-8 moles CO2 m-2 soil surface day-1).

3.1.3. Gas Permeation Barriers

Interestingly, soil penetration data underline the gas permeation pathways established by gas concentration measurements. Soil rigidity data obtained with a penetrologger clearly mark the microsites of low and high gas permeability (Figure 6(b)). Higher soil density or compactness occurs in loamy or otherwise clogged soil parts or in soil with a high skeleton content (stones, stone fragments). These soil characteristics reduce or even do not permit at all the upward migration of the CO2 gas toward the surface.

Following the data from the soil surface down to a soil depth of 80 cm, it becomes clear that the higher the soil resistance, the lower the actual CO2 penetration or concentration. Soil hardness is extremely heterogeneous within the area but increases with soil depth. Due to the inhomogeneity of this parameter, a three-dimensional network of vertical and horizontal soil channels is formed allowing gas penetration only in distinct directions. A three-dimensional pattern of easy penetration and prohibited sites is thus formed. The picture of soil gas concentrations nicely mirrors this diffusional barrier pattern (Figures 2 and 6(b)).

Even the extractable skeleton content in the soil somehow reflects the gas emission pattern (Figure 6(a)). The topsoil has a low content of stones and fragments in high gas emission zones (diagonal line from top left to bottom right), whereas in the upper right corner of the mofette, 50% of dry matter belongs to the soil skeleton. Naturally, the stone fragments are not a perfect hindrance to gas permeation per se. Penetrating gas could easily circumvent isolated stone fragments within the soil. The concentration of stones just reflects the presence of larger blocks and a solid lithosphere in the lower ground.

3.2. Physicochemical Soil Parameters
3.2.1. Soil Water

Figure 6(c) shows values for the soil water content within the mofette area. At the time of our measurements, the water content of the area ranged between 10 and 30%, with an average between 20 and 30%. Dry sites with a water content less than 10% occurred in larger islands on the upper border and a site in the upper right corner. Also, around the lower margin (14-16 m), a drier island was seen. Only three very small islands in the upper left corner contained more water and reached values of up to 60%. The position of these islands corresponds to the position of the three main vents, characterized by high CO2 fluxes and high soil CO2 concentrations at 10 cm depth and thereunder (Figures 2 and 5).

3.2.2. Soil pH

Ascending geogenic carbon dioxide permeates through different soil phases, and due to its high solubility in water, it dissolves in soil water according to temperature, atmospheric pressure, and actual gas concentration and/or flux. Dissolved CO2 rapidly dissociates in water producing protons thus acidifying the aqueous soil phases as long as buffering capacities are below a certain limit [3840]. To check whether the CO2-permeated mofette soil already shows acidification, soil pH was determined within the U1 mofette at each grid intersection (Figure 6(d)). It can be said cum grano salis that there is a reasonable correlation between soil pH and the prevailing CO2 concentrations. The measured values ranged between pH 3.6 and pH 6.4 with a clear diagonal pattern where slightly more acidic soils prevail. Higher pH values (around 5.5 to 6) occur at microsites with lower CO2 concentrations/fluxes. Yet two regions clearly differ. At the right end of the area and between 10 and 20 m (lowest transect) in front of the area, pH values are quite low (pH ca. 4.4-4.0) although [CO2] concentrations are not enhanced. [10, 41] already described this pH-lowering effect with accordant results from mofette fields in Slovenia whereas [42] found no effects for three Japanese CO2 springs [see also [12, 13]]. Our findings are in line with a recent study of [43] which showed that persistent CO2 leakage affects the soil chemical composition, altering soil from slight alkaline to acidic over a period of 38 treatments days. [44] laid a transect through a gas vent in Italy and found that across the transition zone, between vent core and background, soil pH values dropped quickly from 4.5 to 3.8, which resulted in a clear influence on the chemical composition of the soil. Most oxide concentrations decreased, while parameters such as cation exchange capacity, loss on ignition, and total organic carbon began to increase. These chemical soil conditions will certainly affect nutrient availability or uptake [7], root functioning [45, 46], and the above ground processes of photosynthesis [7, 19].

3.2.3. Soil Organic Matter

Stagnant soil water and hypoxic soil conditions should also affect the formation of soil organic matter [SOM; [11, 47, 48]]. Therefore, the total SOM content was analysed in all soil samples (Figure 6(f)). Positive relationships between SOM and soil water have been published several times [4951]. The accumulation of dead biomass in wet soils is thought to be the consequence of a disturbed mineralisation. It may be explained by the oxygen requirement of most litter-decomposing animals, fungi, and bacteria. A sufficient oxygen supply cannot be maintained when the soil pores are filled with water because the diffusion rate of oxygen in water is 10,000 times lower than that in a gaseous medium [52]. In mephitic soils, oxygen is deprived by the pure presence of enhanced CO2. The SOM pattern of the mofette (Figure 6(f)) largely reflects the distribution of soil water (Figure 6(c)) with the highest SOM contents in the upper left and right corner of the mofette area.

3.3. Vegetation
3.3.1. Number of Plant Species

In the present study, 69 different plant species were counted within the mofette area during a quantitative species survey in May 2013 (Table 1). Most species could be determined to the species level. In four juvenile and pathogenically affected specimens, only the genus could be determined (Cirsium, Chaerophyllum, Dryopteris, and Epilobium). The plants found are in good agreement with those found by others in that area [5355]. As the total number of species as well as the soil coverage of each single species was estimated in each 2 m2 square, occurrence and soil coverage of each species could well be correlated with the degassing pattern of CO2 and O2 in the soil.

Comparing the pattern of the number of species per area within the U1 mofette with the prevailing CO2 degassing pattern, it becomes evident that high species numbers (up to 17 per 2m2) only occur at sites of diminished CO2 presence (Figures 2 and 7). Two locations with a high number of species become apparent. On the left upper corner, a diagonal pattern (from 7/16 to 2/2) and a spot on the lower middle (3/30 to 3/44) are evident. Both areas are defined by very low [CO2]. In contrast, the species number is small at sites of higher CO2 emission (CO2 diagonal emission). One to four different species are found at CO2 extremes.

3.3.2. Plant Coverage

The same observation is made when plant coverage is correlated with the CO2 degassing pattern (Figure 7(b)). The plant coverage of 41 to 84% is found only at sites of low CO2 emission. A linear transect from the upper left corner to the upper right corner of the area shows high plant coverages at low CO2 concentrations. A larger island in the middle of the area also shows high plant coverage and low CO2 concentrations at lower soil levels (Figure 7). At high CO2 concentrations (and fluxes), plant coverage can be as low as 4-10%. These results are in line with findings of [44] who reported on a naturally occurring gas vent located within a Mediterranean pasture ecosystem (Latera geothermal field, central Italy). They found no vegetation within the 6 m wide center of the vent where CO2 were determined; an approximately 20 m wide halo surrounding the core formed a transition zone, over which a gradual decrease in [CO2] and a rapid decrease in CO2 fluxes could be observed. In the transition zone, grasses dominated near the vent core. They were progressively replaced by clover and a greater plant diversity distant from the vent center. However, this study had not provided a full characterization of the botanical taxa.

3.3.3. Mofettophilic and Mofettophobic Plants

Correlating the growth patterns of selected plants and the degassing pattern within the designated mofette area, it turns out that just one single plant species clearly follows the degassing pattern of CO2. Marsh sedge (Carex acutiformis Ehrh.) grows in great densities directly above strongly degassing areas within the mofette (Figure 1(a); Figure 8). There is no growth of Carex on the top right and lower left corner with less or nearly no CO2 emission. Being a highlight to half-light loving species, marsh sedge is absent on the top left corner of the area, although CO2 emission is maximal. The dense canopy formed by beech and limewood trees strongly prevents its growth. Mirroring the CO2 degassing pattern, Carex acutiformis is thought to be the only eu-mofettophilic (CO2-indicating) species in the study area. It is worth mentioning that on other sites also Carex species have a high affinity for CO2 and act as geogas bioindicators. In NW-Czech mofette areas, Carex nigra is strictly mofettophilic, whereas in a Yellowstone NP site, Carex aquatica plays that part (Pfanz, Tercek, King, unpubl).

Quite in contrast to the growth pattern of Carex acutiformis are those of grove windflower (Anemone nemorosa), spring fumewort (Corydalis solida), fig buttercup (Ranunculus ficaria), wood bluegrass (Poa nemoralis), two-flower melic grass (Melica mutica), and addersmeat (Stellaria holostea) (Figure 9).

These species clearly avoid growing on degassing sites and occur only on control plots. The growth pattern of some species marks a linear or half-moon-like structure in the upper left corner of the site (running from 7/10-16 to 2/2). This is true for A. nemorosa, P. nemoralis, S. holostea, and R. ficaria (Figure 9). To a certain extent, it is also true for Carpinus betulus, Melica mutica, and Corydalis solida (not shown). The upper right part of the area and the lower part in front are also covered.

Two other species like common snowberry (Symphoricarpos albus) and Robert geranium (Geranium robertianum) cover the lower part of the nongassing area (Figure 10), while common periwinkle (Vinca minor), bulbiferous coralwort (Dentaria bulbifera), and common honeysuckle (Lonicera periclymenum) occupy the nongassing locations at the right rim of the area. Depending on the species, specific mofettophobic plants overgrew either the whole control site or just the light degassing sites. In some cases, mofettophobes grew only at certain low-gas locations. This observation may be explained by the differences in soil permeability, soil water, soil pH, nutrient availability, and shading by tree canopy.

3.3.4. CCA

Canonical correspondence analyses (CCA) were carried out on all environmental and individual species abundance data [37]. In order to exclude random results, plant species were used only if soil coverage of individuals was higher than 5% per area. Species response curves were performed on plant species, with at least one representative from each of the species groups distinguished in CCA. The category “mofettophilic and mofettophobic” [according to [17]] was assessed from CCA diagrams and species response curves for those species, which had their optimum (higher densities) in patches with >2-3% CO2 concentrations.

The canonical correspondence analyses CCA (Figures 11(a)11(b)) revealed strong correlations between environmental measures and plant species.

The strongest environmental vectors separated the plant species into two to three species groups (Figure 11). Group 1 corresponded to high CO2 concentration and consisted of one single species, namely, marsh sedge (Carex acutiformis; Figure 8). In some analyses, also, Tilia platyphyllos was added to that group. In contrast, group 2 was grouped clearly at the lowest CO2 emission sites and consisted of more species (Vinca minor, Symphoricarpos albus, Stellaria holostea, Poa nemoralis, and Melica mutica). Group 3 consisted of varying plant species and was grouped around either slightly acidic soil pH values, soil moisture, SOM, or soil density. Depending on analysis, several species stayed near the plot origin of the diagram, revealing neither positive nor negative correlations to any of the respective environmental measures.

CCA thus corroborate the findings of the distribution of species according to the geogenic soil gas emission. Marsh sedge proved to be the only real “eu-mofettophilic” plant species within the measured area.

4. Conclusion

Postvolcanic CO2 gas leakage still occurs around the Laacher See resulting in distinct, dry CO2 emanations sites (mofettes). Our results indicate that this CO2 degassing significantly affects the terrestrial soil ecosystem. The living organisms within this ecosystem appear to have adapted to the locally high CO2 concentrations through species substitution or adaption, with a shift towards hypoxic and acidophilic adaptations. In the case of plants, an azonally growing helophyte, namely, marsh sedge (C. acutiformis), is able to grow on highly gas-emitting soils. Its ability to supply oxygen to its roots enables the sedge to establish at hypoxic and concomitantly acidic soils. In this area, C. acutiformis is the only really eu-mofettophilic plant species which can be used to biologically indicate geogenic CO2 emissions. Yet in other mofettic places, there are different Carex species showing a high affinity for CO2; they act as geogas bioindicators. In NW Czech mofette areas, Carex nigra is strictly mofettophilic, whereas in a Yellowstone NP site, Carex aquatica plays that part (Pfanz, Tercek, King, unpubl).

The competent, botanical characterization of mofette sites and the identification of eu-mofettophilic (CO2 indicating) and mofettophobic (CO2-avoiding) plant species may help to map gas-emitting tectonic structures with a combination of vegetation data, supported by [CO2] measurements, and may provide easily detectable plant indicators for geogenic (pre- or postvolcanic) CO2 emission.

Whether knowledge on CO2 gas indicator plants will also help to identify gas leaks in artificial CO2 storage fields (CCS) is still a matter of debate.

Data Availability

We cannot share the data with other colleagues at present on a public scale as we are producing more scientific papers with these data. In this case, we have a time series over several years where these data are included. After publishing the follow-up papers, we may allow public access to the data.

Additional Points

Highlights. Indicative plants can hint geogenic carbon dioxide exhalations—mofettophilic plants. Most plants avoid growth on high degassing sites—mofettophobic plants. Gas diffusion clearly follows cracks and soil fissures and is not homogeneous within the area. Geogenic soil CO2 gas flux clearly follows gas concentration patterns.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

We want to thank Christa Kosch for her extremely valuable help in the field as well as in the laboratory. The help of Nina Hennigfeld, Annika Pelz, and Christian Baakes is gratefully acknowledged. Sabine Kühr helped in various soil analytics. Cordial thanks are also due to the Deutsche Vulkanologische Gesellschaft (Heinz Lempertz, Wolfgang Riedel, Rainer Hippchen, Wolfgang Kostka, and Walter Müller) as well as forester Karl Hermann Gräf for their help in getting necessary permissions as well as practical and financial support. We are indebted to Mrs. Claudia Uhl and Mr. Stefan Backes, Struktur- und Genehmigungsdirektion Koblenz Nord, Referat Naturschutz, for the permission to study the area and to the Department of Geosciences of the Friedrich Schiller University of Jena (Profs. Heide, Büchel, and Viereck) for many scientific discussions.