06 Fakultät Luft- und Raumfahrttechnik und Geodäsie
Permanent URI for this collectionhttps://elib.uni-stuttgart.de/handle/11682/7
Browse
13 results
Search Results
Item Open Access Assessing hydrological changes in the Murray-Darling basin(2010) Schmid, StefanieIn this diploma thesis GRACE-derived time-variable gravity, TRMM precipitation, and in-situ hydrological observations (water level, river flow) between May 2003 and April 2009 (6 years) are used to assess hydrological changes in the Murray-Darling Basin (MDB). The investigations in terms of Multiple Linear Regression Analysis (MLRA) and Principal Component Analysis (PCA) show distinct differences between the northern and southern climatic regions of the MDB. While the north is dominated by precipitation occurring within the early part of each year, precipitation appears in the south half a year later. Furthermore, a severe drought is detected, which appeared in 2006/07. The analysis of seasonal variations reveals that gravity changes are preceded by precipitation by about one month in the north of the MDB, while the south shows the inverse behaviour. The phase shift of gravity changes between the north and the south suggests a mass transport taking about 5 months to cross the MDB from the north to the south, which is confirmed by the in-situ hydrological observations.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 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 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 Nachweis einer vulkanischen Struktur auf der Schwäbischen Alb mittels gravimetrischer Messungen und unterstützender Höhenbestimmung(2010) Scheider, AnnetteIm Bereich der Schwäbischen Alb und des Albvorlandes existieren eine Vielzahl von vulkanischen Strukturen unterschiedlicher Größe, die man zusammen als Schwäbischer Vulkan bezeichnet. Sie entstanden bereits im Miozän und sind heute oft nicht mehr als solche zu erkennen. Besonders viele davon befinden sich im Gebiet von Kirchheim u.T. und Bad Urach. Das Maar am Hofbrunnen ist eine dieser Strukturen und liegt in der Gemeinde Münsingen, nahe des Ortes Trailfingen auf der Schwäbischen Alb. Diese Struktur ist bereits seit langer Zeit bekannt, dennoch ist sie in einigen neueren Karten nicht enthalten. In dieser Diplomarbeit soll ein Teilgebiet im Süden des Maars untersucht werden und es wird überprüft, ob eine vulkanische Struktur mittels gravimetrischer Messungen und unterstützender Höhenbestimmungen nachgewiesen werden kann. Zunächst werden die geologischen Grundlagen, die die Entstehung dieser vulkanischen Strukturen und speziell zu der des Maars am Hofbrunnen erklären, betrachtet. Da sowohl Höhen- als auch Schweremessungen durchgeführt werden sollen, werden anschließend die verschiedenen möglichen Messverfahren genauer betrachtet und die für diesen Zweck am besten geeigneten ausgewählt. Um die notwendigen Messungen durchführen zu können, muss im nächsten Schritt ein geeignetes Messkonzept aufgestellt werden, das dann bei der Durchführung der Messungen im Gelände umgesetzt werden kann. Bei der Auswertung wird zuerst aus den erhaltenen Höhendaten ein DGM generiert, dass die Ausgangslage für die Auswertung der Schwere-Messergebnisse liefert. Für die Herleitung der Schwereanomalien auf das Geoid müssen einige Reduktionen und Korrekturen angebracht werden, so dass im Ergebnis die vulkanische Struktur sichtbar wird. In einem letzten Schritt wird der Einfluss der vulkanische Struktur, die eine Störmasse im Untergrund darstellt, auf die Schwere modelliert. Dazu wird die Struktur aufgeteilt, vereinfacht und als geometrische Körper aufgefasst, die aus einem Gestein bestehen, dessen Dichte sich vom Umgebungsgestein unterscheiden. Die Ergebnisse dieser Modellierung werden abschließend mit den Messergebnissen verglichen.Item Open Access GOCE data analysis: optimized brute force solutions of large-scale linear equation systems on parallel computers(2010) Roth, MatthiasThe satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer) was set up to determine the figure of the Earth with unprecedented accuracy. The sampling frequency is 1 Hz which results in a massive amount of data over the one year period the satellite is intended to be functional. From this data we can setup an overdeterminded linear system of equations to estimate the geopotential coefficients which are required for modelling the Earth's gravity field with spherical harmonics in the desired high resolution. The linear system of equations is solved "brute-force" which means that the normal equation matrix has to be kept in memory as a whole. The normal equations matrix has a memory demand of up to 65 GByte, hence we need a computer providing a sufficient amount of memory and fast multiple processors for the computations to get them done in a reasonable time. In this study, a program was written to compute the geopotential coefficients from simulated GOCE data, as GOCE real data were not available yet. As a first step, the program was optimized for the computations to become more efficient. As a second step, the program was parallelized to speed-up the computations by two different techniques: For a first version, the parallelization was done via OpenMP which can be used on shared memory systems which usally have a small number of processors. For the second version, MPI was used which is suited for a distributed memory architecture, hence can incorporate much more processors in the computations. In summation, we gained a huge boost in efficiency of the program due to the optimization. Furthermore, huge speed-up was gained due to the parallelization. The more processors are incorporated in the computation, the more the overall efficiency drops because of increasing communication between the processors. Here we could show that for huge problems the MPI version is running more efficient than the OpenMP version.Item Open Access Untersuchungen zur gravimetrischen Inversion mittels genetischer Algorithmen(2010) Körner, MarkusDie gravimetrische Inversion ist ein zentrales Problem in der Geophysik. Wiewohl die entsprechenden Forschungstätigkeiten auf eine lange Geschichte zurückblicken, sind mit innovativen Techniken selbst heutzutage neue Lösungsansätze zu finden. Das Grundprinzip der gravimetrischen Inversion ist, aus einer Schwereanomalie den Körper oder die Struktur zu extrahieren, welche für die Anomalie verantwortlich ist. Es existiert eine Vielzahl von Berechnungsmethoden, die es erlauben für (fast) jeden belieben Körper dessen Gravitationskraft zu berechnen. Diese Formeln werden, an entsprechender Stelle, für die in dieser Arbeit verwendeten Körper dargelegt. Durch Programme, die teilweise im Rahmen der Arbeit entstanden, kann der Anwender idealisierte Signale erzeugen und diese gemessenen Werten gegenüberstellen. Eine Möglichkeit beide Messwerte zu vergleichen stellt die Berechnung des rms Wertes dar. Dieser errechnet sich über die Summe der quadratischen Differenzen beider Signale, die durch die Anzahl der Messwerte geteilt und aus der anschließend die Wurzel gezogen wird. Ein kleiner rms Wert indiziert eine hohe Übereinstimmung zwischen simuliertem und gemessenem Signal. Das Hauptproblem dabei ist, dass verschiedene Körper dasselbe Signal erzeugen können. Somit gibt es unendlich viele Körper, welche der Schwereanomalie zugewiesen werden können. Es ist nun das Ziel, im Rahmen der Möglichkeiten genau den Körper zu finden, der am "wahrscheinlichsten" ist. Durch geologische Untersuchungen und daraus resultierenden Randbedingungen können die Möglichkeiten deutlich eingeschränkt werden; die Anzahl der verbleibenden Körper ist jedoch zu groß für gesicherte Aussagen. Ein Verfahren ist zu finden, das in der Lage ist durch "probieren" ein gewähltes Modell immer weiter zu verbessern und sich dem "wahren" Körper immer weiter anzunähern. Informatiker begannen Anfang der 1950er Jahre die sogenannten evolutionären Strategien zu entwickeln. Den Verfahren liegt ein biologischer Evolutionsprozess zugrunde, welche mit Hilfe von mathematischen Operationen umgesetzt ist. Aus diesem Spektrum ist für die Arbeit der Genetische Algorithmus GA gewählt worden, bei dem es sich um einen globalen Optimierer handelt. Besonders wichtig in diesem Zusammenhang ist die Tatsache, dass dieser Optimierer keinerlei Ableitungen der verwendeten Optimierungsfunktion benötigt. Somit ist der GA geeignet, für nicht lineare Probleme, wie im Fall der gravimetrischen Inversion. Eine Einführung in dessen Funktionsweise sowie die Adaption an die Problematik ist Bestandteil dieser Arbeit. Dabei wird vor allem auf die Anwendung des sogenannten hybriden GA eingegangen. Bei diesem gemischten Ansatz wird zuerst der GA durchlaufen und die resultierenden Parameter in einem zweiten Schritt mit einem lokalen Optimierungsverfahren, der Downhill Simplex Methode (DSM), weiter verarbeitet. Unter bestimmten Voraussetzungen, die in der Arbeit betrachtet werden, ist es möglich unter Verwendung dieses Ansatzes die Körperparameter weiter zu optimieren und den rms Wert zu senken. Für Untersuchungen wurden Simulationen zu Inversion mit fehlerfreien als auch fehlerbehafteten Signalen synthetischer Körper betrachtet. Für beide Versionen und unter Verwendung verschiedener Körper bietet der hybriden GA das größte Potenzial gute Ergebnisse, d. h. den geringsten rms Wert, zu liefern.Item Open Access Optimization of the orbit parameters of future gravity missions using genetic algorithms(2011) Ellmer, MatthiasThe best global models for Earth's gravitational field are provided by recent satellite missions like CHAMP, GRACE and GOCE. The research in this thesis is dedicated to finding favorable parameters for possible future multi-pair satellite formations that follow in GRACE's footsteps. A software framework was developed in the Python programming language that allows the study of such satellite formations using a range-acceleration approach to gravity recovery. The computational performance of this framework was improved significantly through several optimization steps, amongst others involving successful parallelization of multiple simulations and the exploitation of a GPU-accelerated version of the mathematical library BLAS. A large number of parameter studies were performed in an attempt to define criteria for the automatic quantification of the objective quality of specific multi-pair mission scenarios. This was achieved by analyzing and assessing several characteristics of the recovered potentials in both the spatial and the spectral domain. Using the custom software framework as well as the developed benchmarks for a recovered potential, a genetic algorithm capable of autonomously finding optimal combinations of multiple satellite formations was implemented. This genetic algorithm was then tested in a simple scenario, intended to find the optimal complimentary formation to a single polar satellite pair modeled after the GRACE twin-satellite mission. The performed genetic algorithm simulation suggests that a pendulum on a ~58° inclined orbit with a long along-track baseline and an opening angle of around 40° would make a good companion for GRACE. While finalizing this thesis a flaw in the custom software was discovered that affected all simulations performed in the context of the parameter studies and the genetic algorithm. This means that the results presented in this thesis should be considered with a critical mindset. The developed methodologies used to arrive at these results are however sound and could be used in future studies.Item Open Access Gravimetrische Untersuchungen und Einflussmodellierung "Stuttgart 21"(2014) Efinger, BenjaminDurch das geplante Großbauprojekt "Stuttgart 21" und den damit notwendigen Tunnelbauarbeiten entsteht nach deren Fertigstellung im Untergrund eine Massenveränderung. Diese Massenveränderung führt zu einer veränderten Schwerkraft, die mithilfe von gravimetrischen Untersuchungen beobachtet werden kann. Da sich die Tunnelbauarbeiten noch im Anfangsstadium ihrer Durchführung befinden, können diese Veränderungen bisher nur modelliert werden. Dazu stehen verschiedene Modelle zur Verfügung, die eine Massenveränderung im Untergrund mit unterschiedlichen geometrischen Formen beschreiben. Diese werden auf Grundlage einer Beispielgeometrie miteinander verglichen und bewertet, um anschließend mit den geeigneten Modellen und der geplanten Tunnelgeometrien den Einfluss der Tunnel zu modellieren. Um abschließend den modellierten Einfluss nach Tunnelfertigstellung mit den dann beobachteten Schwerewerten vergleichen zu können, wurde vor dem Beginn der Tunnelarbeiten mit gravimetrischen Beobachtungen begonnen. Dazu wurde ein geeignetes Profil gewählt, das die Tunnelröhren kreuzt, um so eine generelle Aussage über den Einfluss der Massenveränderung treffen zu können. Wichtig ist hierbei unter anderem die Frage, ob es mit dem vorliegenden Instrument, einem Relativgravimeter, möglich ist, aus den beobachteten Schwerewerten einen Einfluss zu beobachten. Zusätzlich werden die Einflüsse auf den Schweregradienten modelliert und analysiert. Da diese jedoch in einem nicht mehr messbaren Bereich liegen, werden dazu nur Messungen in einem Testumfeld durchgeführt und ausgewertet. Dabei ist interessant, ob sich das Messprinzip für zukünftige Untersuchungen eignet.