06 Fakultät Luft- und Raumfahrttechnik und Geodäsie
Permanent URI for this collectionhttps://elib.uni-stuttgart.de/handle/11682/7
Browse
17 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 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.Item Open Access Sequence stratigraphy, sedimentology and provenance of the Upper Cretaceous siliciclastic sediments of South Jordan(2005) Baaske, Uwe Peter; Krawinkel, Hannelore (PD Dr.)The aim of this study was to establish a sequence stratigraphic framework for the Upper Cretaceous (Albian to Santonian) sedimentary rocks of South Jordan. This was accomplished by detailed sedimentological studies and facies analysis. Furthermore the provenance and weathering history of the siliciclastic rocks was determined by using sedimentary petrography and bulk rock geochemistry. The study area is situated in the southern part of Jordan, about 60 km south of the city Ma'an, between the city Ras en Naqb in the West to Naqb Ataik in the East. Here a succession of Palaeozoic to Cretaceous sedimentary rocks is exposed along a NW-SE to E trending, about 90 km long escarpment. Along this escarpment fifteen vertical sections were logged in detail and samples where taken in the Cretaceous deposits. On the basis of these logs the facies analysis and sequence stratigraphic interpretation was carried out. Selected samples were used for petrographical and geochemical studies. Facies analysis showed that the Upper Cretaceous sediments were deposited in a terrestrial to marine setting. The examined sedimentary succession comprises ten environmental depositional facies, which can be grouped into four major facies associations. These four facies associations are (i) the shelf facies association, (ii) the shoreface facies association, (iii) the marginal marine facies association and (iv) the coastal plain facies association. The repeated interfingering of terrestrial and marine environments shows, that the examined sediments were deposited in a shallow marine siliciclastic to mixed siliciclastic-carbonate coastal system. The facies correlation of the logged sections reveals that the sedimentary succession can laterally be sub-divided into three major depositional systems: (i) a siliciclastic ramp, (ii) a mixed carbonate-siliciclastic ramp and (iii) an alluvial system. Vertically these depositional systems are arranged in repeating cyclothems. An ideal cyclothem consists of a retrograding siliciclastic ramp system, a mixed siliciclastic-carbonate ramp system, a prograding siliciclastic ramp system and is ended by an alluvial plain system. The shelf system was generally wave/storm dominated, whilst the prograding siliciclastic ramp system was supply dominated. The retrograding siliciclastic ramp, as well as the mixed siliciclastic-carbonate ramp system was accommodation-dominated. The coastal system can be divided into a regressive and a transgressive system. Under regressive sea-level conditions sand dominated braid deltas formed. Transgressive sea-level conditions resulted in a coastal system controlled by wave dominated estuaries and barrier/lagoon complexes. On the coastal plain anastomosing to meandering sand dominated river systems with wide overbank areas developed during transgressive times. During regressive sea-level conditions the nearshore river system changed to a more braided one that was directly linked to the braid deltas forming during these intervals. The sequence stratigraphic interpretation of the correlated sections shows that the studied sedimentary succession includes eight third-order sequences. The basal four sequences comprise the Cenomanian to Coniacian Ajlun Group, while the top four sequences are situated within the Coniacian to Santonian Belqa Group. Moreover sequences 1, 2 and 4 can each be further subdivided into two sequences of higher order. Basically all sequences include forced regressive deposits formed during the relative sea-level fall. Characteristic features of forced regressions, like for example sharp based shoreface deposits, are frequently present within the logged sections. Due to the significance of the regressive part of the relative sea-level for the studied sediment succession the sequence stratigraphic model introduced by Plint and Nummedal (2000) was used. This sequence stratigraphic model combines the complete sea-level fall in a new falling stage systems tract (FSST), which opposes LST, TST and HST. Comparison of the sequence stratigraphic interpretation with the major depositional systems shows that the depositional systems correspond very well with the four systems tracts of the applied sequence stratigraphic model. The alluvial plain system developed during the LST, the retrograding siliciclastic ramp system during the TST. The mixed carbonate-siliciclastic ramp system corresponds to the HST and the prograding siliciclastic ramp system is associated with the FSST. Relative sea-level changes and sediment supply were the main controlling parameters for the sequence architecture of the examined deposits. Since subsidence was low and no synsedimentary tectonic influenced the examined sediments, the relative sea-level changes are basically eustacy driven. Therefore the reconstructed sea-level curve resembles strongly global curves. The siliciclastic sediment supply is related to sea-level changes and associated with back stepping fluvial incision during the FSST and LST. Light mineral analysis shows that the Cretaceous sands and sandstones of South Jordan are basically mineralogical mature and can be classified as quartz arenites. Only the fine to medium grained sands and sandstones of the study area exhibit a slightly higher content of feldspar, classifying them as subarkoses. Heavy mineral analysis reveals a dominance of the ultrastable heavy minerals zircon, tourmaline and rutile. Especially zircon is predominant in the examined samples. With a mean ZTR-index of 97.9, the heavy mineral content reflects the mineralogical maturity of the sediments shown by the light mineral analysis. The studied sands and sandstones experienced a rather restricted spectrum of diagenetic processes. The compaction of the sediments was low to moderate, which is for example reflected by high IGV values. This and the lack of destructive grain contacts and pressure solution indicate that the sediments did not encounter deep burial or thick overburden. Concerning the cementation of the examined sediments, the sandstones display three major cementing phases. Sandstones cemented by clay minerals are only weakly lithified and display the highest intergranular porosities of all samples. Sandstones cemented by carbonate minerals (calcite and dolomite) or iron-oxides are well lithified and usually display only low values of intergranular porosities. The cementation with carbonate minerals is facies dependant and only occurs in marine sediments. Both, petrographic and geochemical provenance analysis point to the deposition of the Cretaceous sandstones of South Jordan within a passive continental margin setting. As source rock for the studied samples a continental block and/or a craton interior is most likely. Furthermore, a clear trend towards a recycled sedimentary source is obvious in the petrographical as well as in the geochemical studies. Especially the presence of resedimented quartz grains and the enrichment of zircon in the sediments support this. The presented study shows that the Arabo-Nubian-Shield and its Palaeozoic sedimentary cover are the source area for the Cretaceous siliciclastic sediments of South Jordan. The examined sediments follow normal weathering trends in the weathering diagrams. The corresponding CIA values point to a moderate weathering for most samples. Only the sandstone and some claystone samples have high CIA values typical of recycled sediments.Item Open Access Die klastischen Sedimente der karbonen Hochwipfel- und Auernig-Formation der Ostkarawanken (Österreich/Slowenien) : Sedimentologie, Geochemie und Provenanz(2001) Kutterolf, Steffen; Krawinkel, Hannelore (PD Dr.)In dieser Arbeit untersuche ich die karbonen Sedimente der Hochwipfel- und Auernig-Formation in den Ostkarawanken (Österreich/Slowenien). Aus den Ergebnissen der Sedimentologie, Geochemie und Petrographie versuche ich, mit Hilfe der Provenanzanalyse die Entwicklung des Ablagerungsraumes und die geodynamische Stellung der Sedimentationsbecken in den paläozoischen Karawanken zu rekonstruieren. Diese werden dann in einen plattentektonischen Zusammenhang eingebunden und eine karbone Beckenentwicklung abgeleitet. Die Änderungen der unterschiedlichen Provenanzsignaturen aus der Dünnschliff- und Schwermineralanalyse sowie der Mineralchemie und Gesamtgesteinsgeochmie lassen es zu, folgende Beckenentwicklung für die Hochwipfel- und Auernig- Formation anzunehmen. Im Liegenden der Hochwipfel-Formation entwickelte sich in Folge von Dehnungen in der kontinentalen Kruste ein Becken, in das vor allem Detritus des kontinentalen Hinterlands gelangt. Im Folgenden wird der Einfluss eines magmatischen Bogens im Randbereich des Ablagerungsbeckens stärker und führt zu einer Mischung der Liefergebiete. Zum Hangenden der Hochwipfel-Formation hin wird dieser aktive magmatische Bereich im Zuge der variszischen Orogenese zergliedert und aufgearbeitet. Der Detritus dieses erodierten magmatischen Bogens dominiert die Sedimente der obersten Hochwipfel-Formation. Bei der Schließung des Hochwipfel-Beckens entstanden kleinen Restbecken, in denen sich der liegende Teil der Auernig-Formation ablagerte. Als diese Restbecken verfüllt waren, entwickelte sich im oberen Teil der Auernig-Formation eine flachmarine, klastische und karbonatische Sedimentation eines Randmeeres. Diese Beckenentwicklung erklärt am Besten die sedimentologischen, petrographischen und geochemischen Merkmale in den Sedimenten der Hochwipfel- und Auernig-Formation der Ostkarawanken, die von mindestens zwei unterschiedlichen Liefergebieten geschüttet werden.Item Open Access Surface deformation analysis of dense GPS networks based on intrinsic geometry : deterministic and stochastic aspects(2007) Moghtasad-Azar, Khosro; Grafarend, Erik W. (Prof. Dr.-Ing. habil. Dr. tech. h. c. mult. Dr.-Ing. E. h. mult)The first step in this study is to review the properties of surface which are inherent to the surface and can be described without referring to the embedding space. In other words, it is a method of differential geometry. The methods of moving frames which allows deformation of surface could be described by its own rights as a more reliable estimate of surface deformation measures. The method takes advantage of the simplicity of the 2D surface versus the 3D Euclidean spaces without losing or neglecting information about the third dimension in the results. Based on this method, deformation can be described by using tangent vectors and the unit normal basis vector (attached to the bodies before and after deformation). However, basis vectors of the deformed configuration will need to complete information of intrinsic properties of the deformed surface. Through this method, regularized Earth's surface is considered as a graded 2D surface, namely a curved surface, embedded in a Euclidean space . Thus, deformation of the surface can be completely specified by the change of the metric and curvature tensors, namely strain tensor and tensor of change of curvature (TCC). The curvature tensor, however, is responsible for the detection of vertical displacements on the surface. The next step of this study is to concentrate the local basis vectors of the deformed surface which can be formulated in terms of the local basis vectors of undeformed surface and curvilinear components of displacement vector. This will provide a representation of the intrinsic geometry of the deformed surface with deriving information about the displacement field. The new formulation of base vectors (for the deformed body) produces meaningful numerical results for the TCC and its associated invariants (mean and Gaussian curvatures). They can propose a shape-classification of the deformed surface based upon signs of mean and Gaussian curvatures which are new tools for studying the Earth's deformation. To enhance our understanding of the capabilities of the proposed method in defining new basis vectors (for deformed body), we present two examples, one with a simulated data set and the other with a real data set. However, through a real data set we demonstrated a comparison between the proposed method with the plane strain model (2D classical method). Dealing with eigenspace components e.g., principal components and principal directions of 2D symmetric random tensors of second order is of central importance in this study. In the third step of this research, we introduce an eigenspace analysis or a principal component analysis of strain tensor and TCC. However, due to the intricate relations between elements of tensors on one side and eigenspace components on other side, we will convert these relations to simple equations, by simultaneous diagonalization. This will provide simple synthesis equations of eigenspace components (e.g., applicable in stochastic aspects). The last part of this research is devoted to stochastic aspects of deformation analysis. In the presence of errors in measuring a random displacement field (under the normal distribution assumption of displacement field), stochastic behaviors of eigenspace components of strain tensor and TCC are discussed. It is performed by a propagation of errors from the displacement vector into elements of deformation tensors (strain and TCC). However, due to the intricacy of the relations between tensor components (strain or TCC) and their eigenspace components, we proceeded via simultaneous diagonalization. This part is followed by a linearization of the nonlinear multivariate Gauss - Markov model, which links the elements of transformed tensors (obtained by simultaneous diagonalization) with the eigenspace components. Then, we set up an observation model based on a linearized model under a sampling of eigenspace synthesis. Furthermore, we establish linearized observation equations for n samples of independent random vectors from transformed tensor elements (under the normal distribution assumption), each with an individual covariance matrix. This will provide us with the second-order statistics of the eigenspace components. Then we estimate the covariance components between transformed tensor elements by Helmert estimator, based on prior variance information. To enhance conceptual understanding of stochastic aspects of deformation analysis, the method is applied to a real data set of dense GPS network of Cascadia Subduction Zone(CSZ). Comparing the results showed that, in general, after estimating the covariance matrix of observations (transformed tensors via simultaneous diagonalization), variances of eigenspace components become smaller. However, in some areas this did not occur, which can be related to an incorrect description of initial accuracies, either too optimistic or too pessimistic.Item Open Access Harmonische Gravitationsfeldanalyse aus GPS-vermessenen kinematischen Bahnen niedrig fliegender Satelliten vom Typ CHAMP, GRACE und GOCE mit einem hoch auflösenden Beschleunigungsansatz(2009) Reubelt, Tilo; Grafarend, Erik W. (em. Prof. Dr.-Ing. habil. Dr. tech.h.c. mult. Dr.-Ing. E.h.mult)Mit dem Start der Satellitenmission CHAMP im Sommer 2000 fing eine neue Ära der Schwerefeldbestimmung an, die inzwischen durch GRACE (seit Frühjahr 2002) erweitert wurde und durch GOCE (ab 2009) komplettiert werden soll. Der große Erfolg dieser Satellitenmissionen liegt neben ausgefeilten Messprinzipien vor allem auch darin begründet, dass mit Beginn von CHAMP erstmalig globale, homogene und flächendeckende Schwerefelddaten erfasst werden konnten. Verschiedene Messprinzipien garantieren eine Modellierung der langwelligen Schwerefeldanteile (Detailstrukturen >= 500 km) mit CHAMP über eine Auflösung bis in mittelwellige Bereiche (Detailstrukturen >= 250 km) bei GRACE bis hin zu hochfrequenten Anteilen (Detailstrukturen >= 150 km) mit GOCE. Die angewandten Messverfahren sind dabei die Bahnvermessung mit GPS bei CHAMP, die Erfassung der Abstandsänderung zwischen zwei hintereinander fliegenden Satelliten mit Hilfe einer hochgenauen Mikrowellenverbindung bei GRACE sowie die differentielle Beschleunigungsmessung zwischen Probemassen (Satellitengradiometrie) bei GOCE. In der vorliegenden Arbeit wird ein neuartiges Verfahren — der Beschleunigungsansatz — zur Gravitationsfeldbestimmung aus der Bahnvermessung mit GPS am Beispiel des CHAMP-Satelliten untersucht. Die eigentlichen Messgrößen bei der Bahnanalyse sind Trägerphasen und Pseudostrecken der GPS-Messungen zwischen den hoch fliegenden GPS-Satelliten und dem niedrig fliegenden CHAMP-Satelliten. Zur Analyse dieser Messgrößen wird traditionell das auf der Integration der Variationsgleichungen beruhende bahndynamische Verfahren angewendet, welches aufgrund der Integration und des nichtlinearen Gleichungssystems sehr rechenaufwändig ist. Aus diesem Grunde wurden im Hinblick auf CHAMP alternative und effiziente Analysemethoden wie das Energieintegral, das Randwertproblem für kurze Bahnbögen und der Beschleunigungsansatz entwickelt und untersucht. Diese alternativen Verfahren setzen allerdings voraus, dass zuvor die kinematische Bahn des CHAMP-Satelliten aus den GPS-Messungen bestimmt wurde. Fortschritte in der kinematischen Bahnbestimmung, die eine Genauigkeit kinematischer Orbits von 1–3 cm ermöglichen, motivieren zusätzlich den Einsatz der alternativen Analyseverfahren. Bei dem Beschleunigungsansatz werden zunächst aus kinematisch bestimmten CHAMP-Bahnen mittels numerischer Differentiation die auf den Satelliten wirkenden Beschleunigungen berechnet. Es zeigt sich dabei, dass aufgrund der zeitlichen Korrelation der Fehler kinematischer Bahndaten der rauschverstärkende Effekt der numerischen Differentiation stark abgeschwächt wird und somit die Satellitenbeschleunigungen genau genug bestimmt werden können. Nachdem diese Beschleunigungen von gravitativen und nicht-gravitativen Störeffekten bereinigt worden sind, können die gesuchten Kugelfunktionskoeffizienten des Gravitationsfeldmodells direkt durch Anwendung der Newton’schen Bewegungsgleichung bestimmt werden. Zur Lösung des großen linearen Gleichungssystems (für 2 Jahre CHAMP: ca. 6 Mio. Beobachtungen, 8278 Unbekannte für Grad und Ordnung 90) wird die hinsichtlich des Speicherplatzes und der Rechenzeit effiziente Methode der Präkonditionierten Konjugierten Gradienten vorgeschlagen und verwendet. Ein wichtiger Aspekt bei der Analyse von kinematischen Bahnen ist die Datenvorverarbeitung, da die Orbits Ausreißer enthalten, welche die Genauigkeit der Gravitationsfeldschätzung deutlich verschlechtern können. Es wurden verschiedene Verfahren zur Ausreißerelimination getestet, die entweder auf die Varianz-Information der kinematischen Bahnen zurückgreifen, oder Referenzinformation in Form von (reduziert) dynamischen Orbits und bereits bestehenden globalen Gravitationsfeldmodellen benötigen. Als überlegen gegenüber den Methoden der Datenvorverarbeitung erweisen sich die robusten Schätzer. Diese benötigen keine Referenzinformation, stattdessen werden ungenaue Beobachtungen iterativ heruntergewichtet. Die Genauigkeit und Effizienz des Beschleunigungsansatzes wird in dieser Arbeit anhand der Analyse von simulierten und realen 2-jährigen kinematischen CHAMP-Bahnen untersucht. Vergleiche mit den Ergebnissen aus dem klassischen Verfahren und den weiteren alternativen Methoden zeigen, dass mit dem Beschleunigungsansatz Gravitationsfeldmodelle ähnlicher oder sogar höherer Genauigkeit erhalten werden können. Die Gravitationsfeldmodelle können aus dem 2-jährigen kinematischen Orbit bis ca. Grad 80 mit einem Signal-zu-Rausch-Verhältnis > 1 bestimmt werden und zeigen eine Genauigkeitssteigerung gegenüber dem besten Gravitationsfeld der Vor-CHAMP-Ära, EGM96, bis ca. Grad 65.Item Open Access Die Invariantendarstellung in der Satellitengradiometrie : theoretische Betrachtungen und numerische Realisierung anhand der Fallstudie GOCE(2007) Baur, Oliver; Sneeuw, Nico (Prof. Dr.-Ing.)Die Satellitengradiometrie (Satellite Gravity Gradiometry, SGG) ist die derzeit modernste Technik zur Bestimmung und Modellierung hochauflösender Gravitationsfelder. Sie gründet auf der Beobachtung zweiter Ableitungen des Gravitationspotenzials, welche als Gravitationsgradienten (GG) bezeichnet werden und deren Gesamtheit im Gravitationstensor (oder Eötvös-Tensor) zusammen gefasst ist. Letzterer zeichnet sich durch Symmetrie und Spurfreiheit aus. Technisch realisiert wird die Gradiometrie über skalierte Beschleunigungsdifferenzen zwischen frei fallenden Testmassen. Mittels der Kombination aus sechs dreidimensionalen Beschleunigungsmessern lässt sich der volle Gravitationstensor im dreidimensionalen Raum aufstellen. Herkömmlicherweise erfolgt die Gravitationsfeldbestimmung aus SGG Beobachtungen durch die Analyse einzelner GG. Dabei wird für jeden beobachteten GG der (lineare) funktionale Zusammenhang zu den unbekannten Gravitationsfeldparametern hergestellt. Es ergibt sich folglich für jeden GG eine individuelle Beobachtungsgleichung. Dieser Ansatz wird hier als die klassische Vorgehensweise betrachtet. Sie kommt ohne die Orientierung des Gravitationstensors relativ zum Referenzsystem der Gravitationsfeldmodellierung nicht aus, da die einzelnen GG abhängig von der Lage des Gradiometersystems im Raum sind. Alternativ dazu befasst sich der erste Teil dieser Arbeit mit der theoretischen Gravitationsfeldanalyse basierend auf den Rotationsinvarianten des Gravitationstensors. Diese Größen verhalten sich invariant gegenüber orthogonalen Transformationen und lassen sich damit unabhängig von der Orientierung des Gravitationstensors formulieren. Andererseits ist die Invariantendarstellung mit einer Reihe von Erschwernissen verbunden. Die nicht-linearen Funktionale des Gravitationspotenzials verlangen eine entsprechende Linearisierung. Daran gekoppelt ist ein iterativer Lösungsprozess. Dabei ist eine effiziente Linearisierungsstrategie maßgeblich von drei Faktoren abhängig: einem möglichst kleinen Linearisierungsfehler, schnellem Konvergenzverhalten und einem geringen numerischen Aufwand. Es stellt sich heraus, dass die Linearisierung in Form einer Störungsrechnung alle drei Kriterien erfüllt. Darüber hinaus gründet die Invariantendarstellung auf der Volltensorgradiometrie. Dies impliziert, dass sämtliche GG mit möglichst gleicher Genauigkeit verfügbar sein müssen. Eine derartige Annahme kann jedoch nicht grundsätzlich vorausgesetzt werden. Um die Volltensorgradiometrie allgemeingültiger zu gewährleisten, wird deshalb die synthetische Berechnung unbeobachteter GG untersucht. Aktuelles Interesse erfährt die Satellitengradiometrie derzeit vorrangig durch den geplanten Start der Mission GOCE (Gravity field and steady-state Ocean Circulation Explorer) im Frühjahr 2008. Die numerischen Beispiele der Invariantendarstellung basieren auf einer Simulationsrechnung dieses Szenarios. Die Güte der erhaltenen Lösungen wird gegenüber den Ergebnissen unter Anwendung des klassischen Analyseverfahrens von SGG Beobachtungen abgegrenzt. So widmet sich der zweite Teil dieser Arbeit der rechentechnischen Umsetzung der Gravitationsfeldanalyse. Die auftretenden Gleichungssysteme werden nach der Methode der kleinsten Quadrate gelöst. Neben der direkten Lösungsmethode durch Inversion des Normalgleichungssystems nimmt hier das iterative LSQR (Least-Squares unter Verwendung einer QR Zerlegung) Verfahren eine zentrale Rolle ein. Es stellt nur geringe speichertechnische Anforderungen. Des weiteren ist es hinsichtlich der parallelen Implementierung auf Multiprozessor-Plattformen weitaus systemunabhängiger und effizienter als die direkte Lösungsmethode. Das LSQR Verfahren wird für dessen wirtschaftlichen Einsatz in der Gravitationsfeldbestimmung angepasst bzw. erweitert. Zentrale Aspekte sind in diesem Zusammenhang die Regularisierung und Präkonditionierung. Während die Regularisierung Einfluss auf die Güte der Lösung nimmt, zielt die Präkonditionierung auf das beschleunigte Konvergenzverhalten des iterativen Prozesses ab. Weiterhin werden die entsprechenden Konzepte und Ergebnisse der parallelen Implementierung aufgezeigt. Die Umsetzung der Algorithmen erfolgt auf Hochleistungsrechnern des Höchstleistungsrechenzentrums Stuttgart (HLRS) und des Center for Computing and Networking Services in Amsterdam (SARA). Tatsächlich blieben umfassende Simulationsrechnungen der Invariantendarstellung mit Hinblick auf die Satellitengradiometrie bisher aus. Diese Lücke wird mit der vorliegenden Arbeit geschlossen. Die Kombination aus physikalisch-mathematischer Modellbildung und effizienter numerischer Umsetzung demonstriert letztlich die erfolgreiche Handhabung der Invariantendarstellung in der Satellitengradiometrie.Item Open Access Die klastischen Sedimente von Süd-Elfenbeinküste : Provenanz, Umlagerungsprozesse und Entstehung des Goldvorkommens 'Belle Ville'(2001) Wozazek, Stefan; Krawinkel, Hannelore (PD Dr.)In meiner Arbeit untersuche ich Prozesse, die im Gebiet des randlichen 'Côte d’Ivoire'-Beckens die Ausbildung der klastischen Sedimente bestimmten. Für die Untersuchung wählte ich kretazische bis tertiäre Sedimente des Beckens und quartäre, goldführende Ablagerungen nördlich des Beckenrands aus. Die Schwermineralassoziation der Sedimente von Alb bis Turon zeigt nahegelegene Liefergebiete und kurze Transportwege an. Die Schwermineralführung der Sedimente von spätem Cenoman bis Neogen lässt auf eine Vergrößerung des Liefergebiets und auf einen längeren Transport schließen. Der Wechsel in der Schwermineralassoziation während des späten Cenomans zeigt den Übergang von der Syn-Rift-Phase mit einer störungsgebundenen Kliffküste zu der Drift-Phase mit einer Ausgleichsküste auf einem schmalen Schelf. Tertiäre Meeresspiegelschwankungen führten mehrfach zur Freilegung des Schelfs und zu lateritischen Bodenbildungen auf den randmarinen Sedimenten. Aufgrund des Südwest-Monsuns herrschte während des Neogens ein östlich gerichteter Strandversatz vor. Dagegen überwog während des späten Cenoman bis ins Paläogen ein Südost-Passat, der einen westlich gerichteten, lateralen Sedimenttransport verursachte. Nördlich des Beckenrands befinden sich über alteriertem Basement sedimentäre Goldvorkommen. Die Primärmineralisation fand während einer bruchhaften Deformationsphase im Proterozoikum statt. Ein humid-tropisches Klima führte zur Ausbildung von mächtigen Verwitterungsdecken. Das Gold unterlag supergener Alteration. Während einer arideren Klimaphase kam es zur Zerschneidung der Rumpffläche. Das Gold wurde dabei mechanisch umgelagert. Anschließende chemische Verwitterung zeigt die Rückkehr humid-tropischer Verhältnisse an. Die Korrelation von Flächenbildungs- und Einschneidungsphasen zeigt, dass die randliche Zertalung der neogenen Beckensedimente zeitgleich mit der Ausbildung der goldführenden Sedimente während einer interpluvialen Phase im Pleistozän stattfand.