06 Fakultät Luft- und Raumfahrttechnik und Geodäsie
Permanent URI for this collectionhttps://elib.uni-stuttgart.de/handle/11682/7
Browse
32 results
Search Results
Item Open Access Volcanic evolution of Southern Tenerife (Canary Islands) during the Pleistocene and Holocene(2009) Kröchert, Jörg; Buchner, Elmar (PD Dr.)The Canary Islands are a group of volcanic ocean islands in the Central Atlantic near the continental margin of northwest Africa. Tenerife, with a volcanic history of more than 12 Ma of subaerial eruptions, is the largest island of the Canaries and is situated in the centre of the Archipelago. The Quaternary Bandas del Sur Formation in the South of Tenerife comprises a complex sequence of pyroclastic rocks and lavas and is part of the southern rift zone. In contrast to the northwest and northeast rift zones on Tenerife, the southern rift zone comprises a number of characteristics with respect to the morphological features, eruption cyclicity, and the geochemistry of the volcanic deposits. Various flank eruptions of the Las Cañadas volcano associated with basaltic lavas and the formation of cinder cones within the Bandas del Sur are important volcanic units for understanding the explosive volcanic cycles during the Pleistocene on Tenerife. Paleomagnetic studies, geochemical analysis of major and trace elements, and two radioisotopic dating (K-Ar) have been carried out on prominent cinder cones, to determine their stratigraphic position. By combining the results with previous K-Ar data in the Literature, the cones and lavas can be subdivided into three stratigraphic units. Cinder cones that belong to the first unit show reverse magnetization and Y/Nb ratios between 0.37-0.41; cinder cones of the second unit show normal magnetization and Y/Nb ratios of <0.35. The third unit comprises cinder cones with normal magnetization and Y/Nb ratios of about 0.47. The first two units were constructed between ~0.948-0.779 Ma and 0.323-0.300 Ma. These units define volcanic cycles that culminated in violent Plinian eruptions. The third and youngest unit possibly marks the beginning of a further volcanic cycle that started ~0.095 Ma ago. In order to reconstruct the uplift history of Tenerife, numerous uplifted fossil beaches and tuff cones were investigated. In the North and Northeast of Tenerife, the positions of fossil beaches indicate stable conditions since 130 ka. The uplift rates in southern Tenerife (within the Bandas del Sur) amount to a minimum of 15 m since 778 ka at Montaña Pelada and to a maximum of up to 45 m since 10 ka in the area of El Médano, suggesting an asymmetrical uplift of the island complex. The uplift in the South could be caused by seismic activity or mass loss due to flank collapse events. However, uplift due to ascending magma is more plausible. The fossil beach deposits of the El Médano area exhibit tubular-shaped concretions and concretionary dykes. These sediment structures have been interpreted as the result of a) the interaction between hot ignimbrites that overflowed wet beaches, b) fast accumulation of beach sands on hot and degassing ignimbrites, c) paleoliquefaction caused by an earthquake (seismites). Based on the interpretation as seismites, an intense paleoearthquake was proposed to be responsible for the generation of the paleoliquefaction structures. However, the sedimentary structures in question show the general criteria diagnostic for rhizocretions and root tubules with respect to their orientation, size, branching system, and style of cementation. Faults of a well-defined strike direction that precisely coincides with the southern rift fault system occur in the El Médano site. This fault system was generated contemporaneously with a chain of cinder cones ~948 ka ago. Open fractures in ignimbrites (~668 ka) and the fossil beach deposits (~10 ka) of the El Médano area suggest that the rift-associated fault system has been seismically active in the aftermath and probably is still active. A further fault system striking perpendicular to the rift-associated faults probably originates from a Holocene paleoearthquake of moderate intensity. Earthquake-induced ground effects in the fossil beach deposits of the study area are consistent with seismically induced ground effects of several recent and well-documented earthquakes and gravitational sliding triggered by an intense earthquake in Nicoya/Costa Rica in 1990. Both, the rift-associated and the earthquake-induced fault system, initially produced open cracks in the fossil beach deposits that were occupied by plants and subsequently stabilized by cementation. These results accentuate that the densely settled southern part of Tenerife is latently endangered by volcanic and seismic activity, though, currently, there are no indications of increasing volcanic activity in this region. Uplift due to recent magma loading is not observable and the intensity of a paleoearthquake in the El Médano area was probably considerably lower than mentioned in the literature.Item Open Access Rekonstruktion der Ablagerungsverhältnisse im nordalpinen Vorlandbecken Südwest-Deutschlands(2006) Maurer, Holger; Seyfried, Hartmut (Prof. Dr.)Während der Erdöl- und Ergasprospektion zu Beginn der 50er Jahre rückte das Nordalpine Vorlandbecken (Molassebecken) in den Mittelpunkt des Forschungsinteresses. In den kommenden Jahrzehnten wurden zahlreiche Tiefbohrungen im Molassebecken abgeteuft und somit die Grundlage für eine lithostratigraphische Gliederung der sedimentären Abfolge geschaffen. Bis zur heutigen Zeit gibt es nur wenige sedimentologische und sequenzstratigraphische Untersuchungen zum Ablagerungsmilieu der klastischen Sedimente im Nordalpinen Vorlandbecken Süddeutschlands. Das Ziel der vorliegenden kumulativen Arbeit war es, in verschiedenen Themenbereichen mit sedimentologischen, paläopedologischen, geophysikalischen und mathematischen Methoden die Ablagerungsverhältnisse im Molassebecken zu rekonstruieren. In fünf Veröffentlichungen werden verschiedene stratigraphische Altersabschnitte des Nordalpinen Vorlandbeckens in Südwestdeutschland behandelt. Der Schwerpunkt der Untersuchungen wurde dabei auf die ältesten und die jüngsten Sedimentabfolgen im Molassebecken gelegt, die gleichzeitig den Beginn und das Ende der sedimentären Vorlandbecken-Entwicklung repräsentieren.Item Open Access Der Einfluss der kontinentalen Wasserspeicherung auf das Rotationsverhalten der Erde(2008) Hengst, Rico; Wolf, Detlef (Prof. Dr. rer. nat. habil.)Die Schwankungen der Rotationsgeschwindigkeit der Erde und die Richtungsänderungen des Erdrotationsvektors werden mit modernen geodätischen Raumverfahren beobachtet und lassen sich auf Gravitationswechselwirkungen mit anderen Himmelskörpern und auf geophysikalische Prozesse zurückführen. Nach der Reduktion der beobachteten Erdrotationsschwankungen bezüglich der bekannten gravitativen Einflüsse werden die verbleibenden Schwankungen des Erdrotationsvektors maßgeblich durch Massenverlagerungen und Relativbewegungen von Massen in den einzelnen Teilsystemen der Erde, wie z.B. der Atmosphäre, hervorgerufen. Da die reduzierten geodätischen Beobachtungen stets die integrale Folgeerscheinung aller geophysikalischen Prozesse darstellen, sind einzelne ursächliche Anregungen nicht eindeutig identifizierbar. Eine Dekomposition und eine Interpretation des verbleibenden Restsignals erfordert es daher, den Zustand der Teilsysteme mit Messungen physikalischer Größen oder mit Hilfe von numerischen Modellen zu beschreiben. Neben der Analyse von Modellen der Atmosphäre und des Ozeans bezüglich der Erdrotationsschwankungen liegt der Schwerpunkt dieser Arbeit in der Untersuchung von vier hydrologischen Modellen, die die kontinentale Wasserspeicherung simulieren. Im Kontext der kontinentalen Massenverlagerungen werden die hydrologischen Modelle und die hinsichtlich atmosphärisch-ozeanischer Einflüsse reduzierten Schwerefeldbeobachtungen der GRACE-Mission verglichen, wobei sich die Untersuchung nicht auf den globalen Massenumsatz beschränkt, sondern zusätzlich um regionale Analysen erweitert ist. Die ermittelten Differenzen im jährlichen Massenumsatz zwischen den einzelnen Modellen und auch zwischen den Modellen und den GRACE-Daten ergeben mit Hinblick auf die Erdrotationsschwankungen ein unterschiedliches Anregungspotenzial (chi-Funktionen). So treten zwischen den modellierten und den aus Schwerefeldbeobachtungen resultierenden Anregungen Differenzen auf, die in den äquatorialen chi-Funktionen einer Phasenverschiebung der Jahresschwingung von bis zu drei Monaten entsprechen. Wavelet-Analysen der hydrologischen chi-Funktionen zeigen episodische und quasiperiodische Signalanteile auf, die zwischen den einzelnen Modellen signifikante Korrelationen aufweisen. Entsprechende Signalcharakteristika werden auch in den um gravitative, atmosphärische und ozeanische Einflüsse reduzierten Beobachtungen der Erdrotationsschwankungen detektiert. Als Ursachen stellen sich die Oszillationen ENSO (El Niño Southern Oscillation), QBO (Quasibiennial Oscillation), TBO (Tropospheric Biennial Oscillation) und der indische Monsun heraus, die adäquate Variationen in der Wasserspeicherung Südamerikas, Australiens and Asiens bewirken. Um die Übereinstimmungen zwischen den geodätischen Beobachtungen und den modellierten Anregungen aus atmosphärischen, ozeanischen und kontinentalhydrologischen Prozessen quantifizieren zu können, werden die Zeitreihen mit dem Verfahren der spektralen MRA (multiple Regressionsanalyse) untersucht. Im spektralen Band zwischen 10 und 13 Monaten ergeben sich Widersprüche, die auf der Modellierungsseite Probleme in einer der hier untersuchten atmosphärisch-ozeanischen Kombinationen signalisieren, unabhängig von der Wahl der hydrologischen Simulation. Je nachdem welche Modelle bei der spektralen MRA miteinander kombiniert werden, erklären diese im Spektralbereich zwischen 2 und 30 Monaten die Varianz der Tageslängenschwankung im Mittel zu 93% und die Varianz der beobachteten Polbewegung zu durchschnittlich 77%.Item Open Access Empirical orthogonal function analysis of GRACE gravity data(2009) Bentel, KatrinThe Gravity Recovery and Climate Experiment (GRACE) twin-satellite mission has been providing measurements of the time-varying gravity field of the Earth for almost seven years now. Gravity changes on Earth are due to mass changes and play an important role in Earth sciences. Monthly maps of mass changes are derived from the satellite measurements and need to be interpreted. The major difficulty in analyzing GRACE data are North-South stripes in the estimated gravity fields, caused by the fact that the GRACE satellites are flying in a near-polar orbit, one following the other. A microwave ranging instrument is measuring the distance between the two spacecraft, which is about 220 km. Due to these longitudinal stripes, major errors, analyzing the GRACE gravity fields is demanding. The technique of empirical orthogonal function (EOF) analysis is investigated in this thesis, and it is demonstrated the performance of EOF analysis for separating signal from noise and errors, and for identifying different sources of gravity changes in a real GRACE data set. EOF analysis is explained from a theoretical point of view and is applied to the GRACE data. Basically, the EOF method gives a transformation of the data into a new coordinate frame in the data space, where the axis are chosen according to the data variances. The core of the method is a singular value decomposition of the data matrix. The components obtained from this decomposition need to be interpreted, and signal has to be separated from noise. Additionally, EOF analysis can be used as a filtering tool. In the detailed data analysis, benefits and shortcomings of the EOF method are studied and described with respect to GRACE data. Global maps of mass changes as well as different smaller regions are analyzed, and global and regional results are compared.Item Open Access Punktmasseschätzungen auf Grönland aus GRACE(2008) Xia, JingIn dieser Arbeit werden die Massenänderungen über Grönland durch die Analyse von GRACE-Schwereänderungen berechnet. Eine Inversion wird über die Methode der kleinsten Quadrate durchgeführt. Der Zeitraum ist von April 2002 bis März 2007 und enthält insgesamt 57 Epochen (57 GRACE-Monatslösungen). Wenn die Kugelflächenfunktionen bis zum Entwicklungsgrad 30 entwickelt werden und der Regularisierungsparameter zu 10e-43 gesetzt wird, liegt der potimale Fall vor, um die Massenänderung bzw. die Volumenändereung von Eis zu berechnen. In diesem Fall wird die Eisschmelzrate auf Grönland zu -227qkm/Jahr berechnet. Nach dem Vergleich von zwei Massenpunktverteilungen (Verteilung der Massepunkte über ganz Grönland und nur am östlichen Rand Grönlands) wird festgestellt, dass der Eisverlust am östlichen Rand fast 70% des gesamten Eisverlustes ausmacht. Wenn die Massenpunkte auf dem ganzen Land gewählt werden, ist die Massenänderung -227qkm/Jahr; wenn die Massenpunkte am Rand des Landes liegen, ist die Massenänderung -221qkm/Jahr. Wenn der Entwicklungsgrad der Kugelflächenfunktionen größer als 30 ist, wird das Ergebnis nicht stark beeinflusst. Zur Bestimmung des Regularisierungsparameters wird das L-Kurve-Kriterium verwendet. Dieser Parameter beeinflusst die berechneten Massenänderungen stark. Mit Hilfe der L-Kurve wird der Punkt mit der größten Krümmung (Eckpunkt) ausgesucht. Der entprechende Rgularisierungsparameter wird für eine optimale Schätzung gewählt.Item Open Access Einfluss lateraler Variationen in Lithosphäre und oberem Mantel auf den glazial-isostatischen Ausgleich in der Antarktis(2009) Rau, DanielDiese Arbeit beschäftigt sich mit den Auswirkungen der lateral variierenden Lithosphärenmächtigkeit auf die regionale glazial-isostatische Ausgleichsbewegung (GIA) in der Antarktis. Damit einhergehend wird auch die Viskosität des oberen Mantels variiert. Die Lithosphäre in der Antarktis ist in zwei Bereiche mit scharfem Übergang unterteilt, eine dicke Lithosphäre in der Ostantarktis und eine dünne Lithosphäre in der Westantarktis. Nun wird untersucht wie sich die verschiedenen Lithosphärendicken im Modell auf das Verhalten der nacheiszeitlichen Hebung bzw. auf den glazial-isostatischen Ausgleich auswirken. Zur Einordnung werden die Vertikalbewegung, die Horizontalbewegung sowie die Geoidhöhenänderung modelliert. Hierzu wird zunächst ein dreidimensionales (3D) Viskositätsmodell mit einer lateral variablen Lithosphäre erstellt. Dieses Erdmodell wird mit einem Lastmodell belastet, welches sich aus der Vereisungsgeschichte ICE-5G ergibt. Die Berechnung der in dieser Arbeit untersuchten Größen erfolgt unter Verwendung der spektralen Finite-Elemente-Methode für 3D viskoelastische Belastung. Durch ein inhomogenes Abschmelzen der Eismassen in der Antarktis sind dessen Auswirkungen in der Westantarktis stärker ausgeprägt. Um die Effekte der Lithosphärenvariation zu extrahieren, werden Vergleichsmodelle erstellt, die beispielsweise einen 3D-Mantel, jedoch nur eine eindimensionale (1D) Lithosphäre aufweisen. Anhand der durchgeführten Vergleiche wurde deutlich, dass der Einfluss der Lithosphärenmächtigkeit in Relation zu den Einflüssen des Mantels sehr gering ist. Lediglich im Bereich der Horizontalbewegungen ist ein Einfluss von 20-30 % der Gesamtgeschwindigkeit messbar. Für die Vertikalbewegungen bleibt der Einfluss mit < 1 mm/a unter 10 %, und für die Geoidhöhenänderungen ist die Änderung aufgrund von Lithosphärenvariationen < 0.10 mm/a. Der Einfluss der Mantelviskosität ist in allen Bereichen höher.Item Open Access Variance-covariance matrix estimation with LSQR in a parallel programming environment(2008) Guo, RonggangKnowledge about the gravity field allows an insight into the structure and dynamics of the earth. It provides the geoid as the most important physical reference surface in geodesy and oceanography. Since 2000, the CHAMP (CHAllenging Mini-satellite Payload) mission detects the structure of the global gravity field, followed by the launch of GRACE (Gravity Recovery And Climate Experiment) in 2002. In 2008, finally, the GOCE (Gravity field and steady-state Ocean Circulation Explorer) satellite is supposed to be set in orbit. These missions demonstrate satellite-based gravity field recovery to be at the center of geo-scientific interest. Interpretation and evaluation of satellite observations are difficult, especially the determination of the unknown gravity field parameters from a huge amount of measurements. Because of the immense demand for memory and computing time, the occurring systems of equations pose a real numerical challenge. Therefore, High-Performance Computing (HPC) is commonly adopted to overcome computational problems. Basically, parallel programming with MPI and OpenMP routines allows to speed up the solution process considerably. In this thesis, firstly global gravity field modelling by means of satellite observations is reviewed. Secondly, the LSQR method (Least-Squares using QR factorization) is introduced in detail in order to solve the resulting least-squares problems. Because the LSQR method is an iterative solver, it basically can not provide the variance-covariance information of the parameter estimate. To investigate the approximate computation of the variance-covariance matrix, two methods are introduced. The first one is based on the generalized inverse of the design matrix. The second approach applies Monte-Carlo integration techniques. Because parallel programming is very helpful to implement such iterative methods, it is necessary to introduce some basic principles and concepts about HPC.Item Open Access GOCE data and gravity field model filter comparison(2008) Raizner, CarinaNew approaches with respect to space borne gravity observations are expected to significantly improve the overall knowledge of the Earth's gravity field and its geoid. The Gravity field and steady-state Ocean Circulation Explorer (GOCE) is the first Core Earth Explorer Mission of the ESA Living Planet Programme. This new satellite mission based on the concept of satellite gradiometry is designed to support applications in Earth physics, oceanography and geodesy with an accurate and detailed global model of the Earth's gravity field and its geoid. One of the main problems in the use of the GOCE data is that the retrieval algorithms need along-track filtering on one hand and/or the implementation of spherical filters on the other. The match between these along-track one-dimensional filters and the spherical two-dimensional ones is far from obvious. Thus, the objective of this study is to investigate the influences of these two filter types by analyzing the differences between simulated GOCE reference and filtered data. Apart from closed-loop tests in order to check the consistency and correctness of the data and software used, the testing procedure for along-track as well as spherical filtering is implemented as follows. First, a global reference model is used for data generation which yields a reference signal along the orbit. By applying a one-dimensional along-track filter to these synthetic satellite data, a filtered global model is retrieved. On the other hand, the synthetic satellite data can be also generated after applying spherical filters to the global reference model. The outcome is a filtered global model estimated from these synthetic satellite data. The influences of both filter types are assessed by comparing the reference and filtered signals along the orbit as well as by comparing the reference and filtered models on the ground. Additionally, the properties of both filter types can be varied. In order to examine the empirical relation between along-track and spherical filters, transfer functions of the filters are investigated in a second step of this study. The transfer function for the spherical filter in the model domain is the ratio between reference and filtered signal which represents a corresponding one dimensional along-track filter in the signal domain. On the other hand, computing the ratio between reference and filtered model estimated from the along-track filtered signal relates the one-dimensional filter in the signal domain to a two-dimensional spherical filter in the model domain. The outcome of the study will be very useful for explaining some of the differences between current global model retrieval philosophies and will also be applicable to other satellite missions and data types in the future.Item Open Access GOCE sensitivity studies in terms of cross-over analysis(2009) Xue, YangThe GOCE (Gravity field and steady-state Ocean Circulation Explorer) satellite, launched on 17 March 2009, for the first time applies satellite gravity gradiometry (SGG) to recover the Earth's gravity field with cm accuracy at a resolution of 100km. To meet the envisaged accuracy, measurement validation at cross-over points (XOs) is necessary. Typically, validation is based on gravity gradients (GGs). However, the coefficient matrix of the gravitational tensor is dependent on orientation. In order to avoid matrix rotation, analysis based on orientation-independent invariants is possible. By applying various noise models, the goodness of XO-validation based on GGs and invariants will be studied in this thesis. First, by determining the maximum of scalar products from two tracks, the XOs can be predicted. Next, using local polynomial approximation, the geographical coordinates of XOs are calculated by solving a system of equations. Due to the orbit drift, the interpolation of height is performed separately along ascending and descending track before final comparison. Considering a sampling rate of 1Hz, GGs and invariants in all points of a one-week orbit are simulated for the further interpolation at the XOs. To determine the goodness of the selected interpolation algorithm, a closed loop test with noise-free data is investigated first. Since signal to noise ratios of GGs and invariants are all above 70dB, the same algorithm is applied in closed loop tests with noisy data. Since GOCE can only provide high accuracy for the main diagonal tensor components, various noise models, i.e. homogeneous and inhomogeneous white noise as well as homogeneous and inhomogeneous coloured noise, are added to the simulated values. The comparison of the goodness of GGs opposed to invariants is based on the signal to noise ratio (SNR). In this study, the second invariant demonstrates better SNR than GGs and the third invariant in the case of homogenous noise. However, due to the impact of inaccurate GGs, the SNR of invariants is poorer than the SNR of all GGs in the case of inhomogeneous noise.Item Open Access Dreidimensionale Finite-Elemente-Simulation der Standsicherheit von Auslaugungshohlräumen und deren geologische Bewertung (Gipskeuper-Formation, Stuttgart-Bad Cannstatt)(2008) Schweikardt, Steffen; Rogowski, Eckard (Prof. Dr. rer. nat.)In Stuttgart-Bad Cannstatt kommt es seit Jahrzehnten durch Auslaugung gipsführender Gesteine zu Erdfällen und Bauschäden. Erdfälle entstehen aus dem vertikalen Verbruch natürlicher Hohlräume bis zur Geländeoberfläche. Der letzte Erdfall ereignete sich in der Nacht vom 12. auf den 13. Mai 2000 auf dem Freigelände eines Kindergartens. Dieser Erdfall war Anlass für eine umfangreiche Untersuchungskampagne mit geoelektrischen, seismischen und mikrogravimetrischen Methoden. In der Nähe des Kindergartens wurde in einer Kernbohrung innerhalb der Grundgipsschichten des Mittleren Keupers ein Auslaugungshohlraum erbohrt und mit der Sonarmethode vermessen. Diese Daten sind die Grundlage einer numerischen Untersuchung der Standsicherheit mit der Methode der finiten Elemente, bei der axialsymmetrische und dreidimensionale Finite-Elemente-Netze erstellt und idealisierte Fallbeispiele berechnet werden. Die Randbedingungen bei der numerischen Standsicherheitsprozedur der phi-c-Reduktion werden mit Finite-Element-Netzen bestimmt. Variablen sind die im Stoffgesetz nach Mohr-Coulomb verwendeten Bodenparameter, die Größe des Finite-Elemente-Netzes, die Geometrieelemente des Auslaugungshohlraums (Kontur und Radius) und die Tiefenlage des Hohlraumes. Im Rahmen von Berechnungsreihen mit idealisierten Fallbeispielen werden die Einflüsse des Grundwassers und der Lockergesteinsüberdeckung auf die Standsicherheit untersucht. Im letzten Schritt wird eine bessere Annäherung der numerischen Simulation an die natürlichen Verhältnisse erreicht, indem unterschiedlich lange Hohlraumachsen in das dreidimensionale Finite-Elemente-Netz Eingang finden. Die Berechnungsergebnisse zeigen, dass die Scherparameter Kohäsion und Reibungswinkel, wie sie im Stoffgesetz nach Mohr-Coulomb verwendet werden, maßgeblich die Stabilität der Auslaugungshohlräume bestimmen. Die Zugspannung ist ein weiteres wichtiges Kriterium für die Entwicklung des Verbruchs. Erdruhedruckbeiwert und Dilatanzwinkel haben keinen wesentlichen Einfluss auf die Standsicherheit. Die mechanischen Gebirgseigenschaften bestimmen die Standsicherheit. Die geometrischen Bedingungen (Hohlraumgröße und Mächtigkeit der Überdeckung) entscheiden über die Verbruchsart bei unterschrittener Standsicherheit. Bei tiefliegenden Hohlräumen ist lediglich ein Gewölbeverbruch (Verbruch des Gesteins im Entlastungsgewölbe) über dem Hohlraumdach vorhanden, bei größeren Hohlräumen oder geringerer Überdeckung geht dieser in einen Schlotverbruch (zylinderförmiger Verbruch des Gesteins zur Geländeoberfläche) über. Der Auslaugungshohlraum verändert seine Tiefenlage zuerst durch Hochbrechen im Festgestein, und zwar je nach Stabilität der Gesteine von Schicht zu Schicht und schließlich durch Nachbrechen der hangenden Lockergesteine, bis die Absenkung als Erdfall zutage tritt. Für die analytische Betrachtung wird der stark vereinfachte Bruchmechanismus eines zylinderförmigen Gebirgspfropfens über dem Hohlraum verwendet. Dieser Ansatz stimmt gut mit dem sich abzeichnenden Bruchkörper überein, wie er bei der phi-c-Reduktion in Erscheinung tritt. Die analytische Lösung wird der numerisch bestimmten Standsicherheit angenähert, um eine maximale Standsicherheit für eine beliebige Gesteinsschicht zu errechnen und zu ermitteln, bei welcher Tiefe der Verbruch zum Erdfall eintritt. Im Ergebnis zeigt sich, dass der untersuchte Hohlraum unter den derzeitigen geologischen und geotechnischen Bedingungen stabil ist.