ALBERT

Description: In dieser Arbeit wurden Laborexperimente im Ultraschallfrequenzbereich und begleitende theoretische Untersuchungen durchgeführt, um zu klären, welche Absorptionsmechanismen in fluidhaltigen Sandsteinen bei geringem einaxialen (26 MPa) und bei erhöhtem hydrostatischen Druck (bis 200 Mpa) wirksam sind. Die Kenntnis der Absorptionsmechanismen ist notwendig, um die Zusammenhänge von lithologischen und seismischen Gesteinsparametern zu verstehen und für die Interpretation von seismischen Feldmessungen zu nutzen. Bekannt war bislang, daß unter geringem einaxialen Druck frequenzabhängige lokale Flüssigkeitsströmungen den wesentlichen Absorptionsmechanismus in fluidhaltigem Gestein darstellen. Im Rahmen der vorliegenden Arbeit wurde gezeigt, daß die Kompressionswellengeschwindigkeit in fluidhaltigem Gestein bei höheren Sättigungsgraden zusätzlich von einem partiellen Gesteinsversteifungseffekt beeinflußt wird. Damit gelang eine vollständige Erklärung des sättigungsabhängigen elastischen und anelastischen Gesteinsverhaltens bei geringem einaxialen Druck. Um die Bedeutung der Gesteinsmikrostruktur (z.B. Rißdichte, Kornkontakte) für die seismischen Gesteinsparameter (Geschwindigkeit, Absorption) direkt nachzuweisen, wurde diese bei einem Gestein künstlich, thermisch verändert. Die Änderung der Mikrostruktur führt bei gleichzeitiger vernachlässigbarer Änderung der makroskopischen Gesteinsparameter (Porosität, Permeabilität) zu einer Vergrößerung der Fluideffekte. Mechanismen und Modelle, in welche ausschließlich makroskopische Gesteinsparameter ein- gehen, erfassen nicht die beobachteten Änderungen im Gesteinsverhalten. Nur mit dem auf die Gesteinsmikrostruktur bezogenen Mechanismus lokaler Flüssigkeitsströmungen konnten die Änderungen der Fluidwirkungen erklärt werden. Diese Ergebnisse an verändertem Gestein bieten somit einen direkten Beweis für den grundlegenden Einfluß der Mikrostruktur auf die seismischen Gesteinsparameter. Bislang war nicht geklärt, ob der Mechanismus lokaler Flüssigkeitsströmungen, welcher von der Existenz feiner Spalten im Gestein abhängt, auch unter erhöhtem hydrostatischen Druck das Gesteinsverhalten bestimmt. Zur Untersuchung der Absorptionsmechanismen in Gestein unter erhöhtem Druck wurden Experimente an Sandstein mit verschiedenen Saturanden und Sättigungsgraden durchgeführt. Die experimentellen Ergebnisse zeigen, dass auch unter erhöhtem Druck fluidbedingte Absorptionsmechanismen wirksam sind und auch hier polare Fluide im gesamten Druckbereich eine Modulerniedrigung verursachen. Der Vergleich der experimentellen Ergebnisse mit Modellvorhersagen ergab, daß lokale Flüssigkeitsströmungen den Einfluß unterschiedlicher Fluide auf Geschwindigkeit und Absorption erklären und auch unter hohem hydrostatischen Druck von wesentlicher Bedeutung für die elastischen und anelastischen Gesteinseigenschaften sind. Folglich existieren trotz einer partiellen Rißschließung auch unter hohem hydrostatischen Druck noch Risse und Spalten im Gestein. Die thermische Veränderung der Mikrostruktur eines Gesteins bewirkt eine Geschwindigkeitserniedrigung im trockenen Gestein, welche auch bei hohem hydrostatischen Druck erhalten bleibt, sowie eine Vergrößerung der Fluideffekte. Die Erklärung des geänderten Verhaltens des fluidhaltigen Gesteins gelang mit dem Mechanismus lokaler Flüssigkeitsströmungen. Die Ergebnisse an verändertem Gestein zeigen direkt, daß die Mikrostruktur auch bei hohem hydrostatischen Druck grundlegend für die seismischen Gesteinsparameter ist. Da auch unter erhöhtem hydrostatischem Druck, also unter ’in-situ’-Bedingungen, der frequenzabhängige Mechanismus lokaler Flüssigkeitsströmungen die seismischen Gesteinsparameterwesentlich beeinflußt, muß diese Frequenzabhängigkeit bei der Übertragung von Laborergebnissen auf die Feld- und Bohrlochseismik beachtet werden.

Description: Laboratory experiments in the ultrasonic frequency range and accompanying theoretical investigations were performed with the aim to clarify which absorption mechanisms are effective in fluid-containing sandstones at low uniaxial pressure (26 MPa) and at elevated hydrostatic pressure (up to 200 MPa). The knowledge of the effective absorption mechanisms is necessary to understand the relation of lithological and seismic rock parameters and to use it for the interpretation of seismic field measurements. Preceding investigations at the TU-Berlin had revealed that at low uniaxial pressure frequency dependent local-fluid-flow is the dominating absorption mechanism in fluid-containing rock and that polar fluids lead to a modulus reduction. Within the scope of the present study the investigations of the seismic rock parameters (velocity, attenuation) at low uniaxial pressure were continued. It was shown that the compressional wave velocity in fluid-containing rock is additionally influenced by a partial rock stiffening effect at higher saturations. With this result a complete explanation of the saturation dependent elastic and anelastic rock behaviour at low uniaxial pressure was achieved. To prove the significance of the rock microstructure (e.g. crack density, crack aspect ratios, grain contact properties) for the seismic rock parameters directly, the microstructure of a sandstone was artificially changed by thermal cracking. With negligible changes of the macroscopic rock parameters (porosity, permeability) the alteration of the microstructure leads to an increase of the fluid effects on the elastic and anelastic rock properties. Mechanisms and models which are solely based on macroscopic rock parameters do not cover the observed changes in the rock behaviour. Only by applying the mechanism of local-fuid-flow related to microstructure the changes in the fluid effects can be explained. Thus these results obtained on thermally changed rock are a direct prove of the fundamental influence of the microstructure on the seismic rock parameters. Until now it has not yet been clarified whether the mechanism of local-fluid-flow which depends on the existence of narrow cracks in rock also determines the rock behaviour under elevated hydrostatic pressure. To investigate the absorption mechanisms in rock under elevated pressure ultrasonic measurements were performed on sandstone with different saturands and degrees of saturation. The experimental results reveal that under elevated pressure fluid dependent attenuation mechanisms are effective too and polar fluids cause a modulus reduction in the entire pressure range. The comparison of the experimental results with model predictions reveals that local-fluid-flow explains the influence of different fluids on the velocity and attenuation and that this mechanism is of essential significance for the elastic and anelastic rock properties even under high hydrostatic pressure. Consequently despite a partial crack closure sufficient fine cracks exist even under high hydrostatic pressure. The thermal alteration of the microstructure of a rock leads to a velocity decrease in the dry rock which is preserved at high hydrostatic pressure, as well as to an increase of the fluid effects. The explanation of the changed behaviour of the fluid-containing rock was achieved by means of the mechanism of local-fluid-flow. The results obtained on altered rock directly demonstrate that the microstructure is of fundamental significance for the seismic rock parameters at high hydrostatic pressure too. Depending on the rock microstructure the mechanism of local-fluid-flow may cause attenuation and dispersion in all frequency ranges. Since it has been shown in this study that also under elevated hydrostatic pressure, i.e. under ’in-situ’conditions, this frequency dependent mechanism essentially influences the seismic rock parameters, the frequency dependence must be taken into account by extrapolating laboratory results to field and borehole seismic.

Description: thesis

Description: DFG, SUB Göttingen

Description: Es wurden verschiedene nichtlineare Optimierungsverfahren auf die Geometrieparameter von Dichtemodellen angewendet. Dabei waren diese mit Randbedingungen versehen. Ein Vergleich zeigte die Vor- und Nachteile der einzelnen Algorithmen. Anhand von zwei synthetischen Testfunktionen wurde das Verfahren Downhill-Simplex mit der (1,10)-CMA-ES verglichen. Das Simplexverfahren zeigte bei weniger als zehn Parametern ein besseres Konvergenz verhalten als die Evolutionsstrategie. Bei höheren Dimensionen konvergierte die Evolutionsstrategie deutlich besser, während Simplex bei mehr als 40 Parametern nicht mehr konvergierte. Die Evolutionsstrategie zeigt unabhängig von der Anzahl der Parameter logarithmische Konvergenzgeschwindigkeit. Die Evolutionsstrategie, genetische Algorithmen. Simulated-Annealing. Threshold-Accepting und der Deluge- Algorithm wurden danach anhand eines zweidimensionalen Salzstockmodells untersucht. Dabei wurden die Geometrieparameter optimiert und die Dichten des Modells konstant gehalten. Die 52 zu optimierenden Parameter wurden mit Randbedingungen versehen. Die (1,10)-CMA-ES benötigte im Durchschnitt die wenigsten Funktionsaufrufe und erreichte die besten Qualitäten. Ausgehend von verschiedenen Startmodellen konnten ähnliche Lösungen, welche alle sehr gute Qualitätswerte hatten, gefunden werden. Rekombination verbesserte die Konvergenzeigenschaften immer. Evolutionsstrategie war stabil gegenüber numerischen Variationen. Genetische Algorithmen konvergierten zu Beginn des Optimierungsprozesses schneller als alle anderen Verfahren, erreichten aber nie die von der Evolutionsstrategie gefundenen besten durchschnittlichen Qualitätswerte. Downhill-Simplex. Simulated-Annealing. Threshold-Accepting und der Deluge-Algorithm zeigten insgesamt schlechteres Konvergenzverhalten als die Evolutionsstrategie und die genetischen Algorithmen. Eine Anwendung der Evolutionsstrategie auf ein Salzstockmodell zeigte die Probleme der Optimierung in drei Dimensionen. Dennoch konnte das Modell in Hinblick auf die gravimetrische Anpassung deutlich verbessert werden.

Description: Various non-linear optimization techniques were applied to constrained geometry parameters of density models. A comparison showed advantages and disadvantages of each method. By means of synthetical test-functions the downhill-simplex-method has been compared to the (1,10)-CMA-ES. For numbers of parameters smaller ten. better convergence behaviour was achieved with the simplex-method. However, at higher dimensions evolution-strategy converges better than the simplex- met hod. Simplex does not converge for dimensions greater than 40. Independent of the number of parameters, the evolution-strategy showed a logarithmic convergence speed. The evolution-strategy, genetic algorithms, simulated annealing, threshold accepting and the deluge algorithm were investigated using a two-dimensional salt dome model. Here only the geometry of the model was optimized, the densities were held constant. The 52 parameters were constrained. In average the (1,10)-CMA-ES required the less function evaluations and gained the best qualities. Independent of the start configuration evolution-strategy found similar solutions, whose qualities were fairly good. Recombination always improved the results in terms of convergence behaviour. The evolution-strategy was stable against numerical variations. At the beginning of the optimization process genetic algorithms converge faster than all other methods, but they never reach the average quality values which were gained by evolution-strategy. In summary the simplex-method, genetic algorithms, simulated annealing, threshold accepting and the deluge algorithm showed worse convergence behaviour than evolution-strategy and genetic algorithms. An application of evolution-strategy to a three dimensional salt dome model showed the problems encountered in three dimensions. Despite this, in terms of gravity fit. the model could be improved noticeable.

Description: thesis

Description: DFG, SUB Göttingen

Description: Mit einer dreidimensionalen Vorwärtsmodellierung von Dichtestrukturen der Kruste und des oberen Mantels im Gebiet von 12°-35°S und 57°-79°W wird der regionale Trend des Schwerefelds und gleichzeitig das Geoid am aktiven Kontinentalrand Südamerikas zwischen 20° S und 29°S erklärt. Das Dichtemodell umfaßt die abtauchende Nazca-Platte, Teile des südamerikanischen Kratons, den Asthenosphärenkeil zwischen Unter- und Oberplatte und die Kruste der Zentralanden. Eine Vielzahl externer geowissenschaftlicher Randbedingungen schränkt die Dichteverteilung ein. Die Strukturierung der zentralandinen Kruste im Modell orientiert sich hauptsächlich an den Ergebnissen der Refraktionsseismik. Zur Ermittlung der Krustendichten wird eine druck- und temperaturabhängige Geschwindigkeits-Dichte-Beziehung (S. Sobolev) zur Anwendung gebracht. Neben der bekannten Berechnung der Schwerewirkung dreidimensionaler Dichteverteilungen erlaubt eine neuentwickelte Berechnung des Schwerepotentials nun auch die Modellierung des Geoids. Das Einbeziehen eines Referenzdichtemodells ermöglicht die Anpassung der berechneten Felder an die beobachteten Felder bei gleichzeitiger Verwendung von Absolutdichten. Verschiedene Dichtekontraste der abtauchenden Nazca-Platte gegenüber dem umgebenden Mantel sowie mögliche Dichtekontraste innerhalb des Slab, die sich aus Phasenumwandlungen der gesteinsbildenden Minerale ergeben, werden behandelt und ihre einzelnen Beiträge zum Schwerefeld und Geoid der zentralen Anden berechnet. Mit dem fertigen Gesamtmodell wird die Zusammensetzung des Schwerefeldes und des Geoids am aktiven Kontinentalrand Südamerikas analysiert und ihre Hauptbeiträge dargestellt. Die im Modell enthaltenen Strukturen erklären das regionale Schwerefeld und das Geoid vollständig. Daraus folgt, daß aus tieferliegenden Dichteinhomogenitäten, beispielsweise an der Grenze vom oberen zum unteren Mantel oder im unteren Mantel aufgrund einer möglichen isostatischen “in situ”-Kompensation keine weiteren Beiträge zu erwarten sind. Das Schwereminimum wird hauptsächlich durch eine bis auf 65km verdickte Kruste erklärt. Die Schweremaxima in der Küstenkordillere werden durch den Beitrag des Slab und eine Zone erhöhter Dichte in mittlerer Krustentiefe verursacht. Eine anomale VΡ-ρ-Beziehung unter dem rezenten Arc läßt dort partielle Aufschmelzung vermuten. Eine Untersuchung des Schwimmgleichgewichts des Modells ergibt isostatische Unterkompensation (= Massenüberschuß) im Bereich der Küstenkordillere und der Ostkordillere sowie isostatischen Ausgleich bis leichte Überkompensation im Bereich der Westkordillere und des Altiplano bzw. der Puna.

Description: A three-dimensional foreward modelling of the density structures of the crust and the upper mantle in the region between 12°-35°S and 57°-79°W explains both the regional trend of the gravity field and the geoid at the active continental margin of South America from 20°S to 29°S. The density model comprises the downgoing Nazca plate, parts of the South American craton, the asthenospheric wedge between lower and upper plate, and the crust of the Central Andes. A large number of external geoscientific boundary conditions constrains the density distribution. The Central Andean crust of the model is structured according to the results of refraction seismic studies. A pressure and temperature dependent velocity- density relation (S. Sobolev) is applied to determine densities within the crust. Apart from the well known calculation of the gravity effect of three-dimensional density distributions, a recently developped calculation of the gravity potential now allows the modelling of the geoid. Using absolute densities, the fitting of the calculated fields to the observed fields is realized by the application of a reference density model. Several density contrasts of the subducting Nazca plate against the surrounding mantle and possible density contrasts inside the slab resulting from mineral phase transitions are discussed and their particular contributions to the Central Andean gravity field and the geoid are calculated. Using the complete model, the composition of the gravity field and the geoid at the active continental margin of South America is analyzed and the main contributions are visualized. The model structures explain both regional gravity field and geoid completely. This implies that density inhomogeneities at greater depth than presented in the model, e. g. at the interface between the upper and the lower mantle or in the lower mantle are not expected to have effects at the surface due to possible isostatic “in situ” compensation. The gravity minimum is generally explained by a crust thickened to 65 km maximum. The gravity highs in the Coastal Cordillera are caused by the contribution of the slab and a zone of increased density in mid crustal layers. An anomalous VΡ-ρ- relation beneath the recent arc implies partial melt. The investigation of the isostatic state of the density model results in isostatic undercompensation (= mass surplus) in the region of the Coastal Cordillera und the Eastern Cordillera, and isostatic equilibrium to slight overcompensation in the Altiplano/Puna region.

Description: thesis

Description: DFG, SUB Göttingen

Description: Mit dem Ziel, aktuelle Arbeitshypothesen zum geologischen Aufbau der Südural-Lithosphäre zu überprüfen, wird das Schwerefeld im Bereich des Südurals mit folgenden gravimetrischen Methoden untersucht: • Numerische Analyse des Schwerefeldes mit Hilfe von Tiefenabschätzungen, Wellenlängenfilterungen und des Parker- Algorithmus sowie • zwei- und dreidimensionale Dichtemodellierungen. Im Vergleich mit dem hauptsächlich eingesetzten Verfahren zur Dichtemodellierung kann festgestellt werden, daß direkte gravimetrische Auswerteverfahren, wie z.B. die Wellenlängenfilterung, für komplizierte geologische Strukturen dem indirekten Auswerteverfahren mittels Dichtemodellierung weit unterlegen sind. So ist es z.B. mit der Wellenlängenfilterung nicht möglich, langwellige negative Schwerewirkungen der paläozoischen Sedimente von den bivergenten Unterkrustenstrukturen zu trennen. Auch mit variablen Dichtekontrasten bezüglich des Krusten-Mantel-Übergangs bei der Anwendung des Parker-Algorithmus ist es nicht möglich, die Schwerewirkung dieser Grenzfläche im gesamten Untersuchungsgebiet abzuschätzen. Die Verbindung und Interpretation aller verfügbaren Randbedingungen aus der Reflexions- und Refraktionsseismik sowie den zur Verfügung stehenden Datensätzen zur Sedimentbedeckung und Moho-Tiefe führt zu komplexen dreidimensionalen Dichtemodellen der Südural-Lithosphäre. Diese führen mit den derzeit verfügbaren Ergebnissen über die strukturellen und physikalischen Eigenschaften der Südural-Lithosphäre im Bereich von 52° - 65° O und 50° - 55,5° N auf zwei gleichberechtigt nebeneinanderstehende Dichtemodelle: 1. Das Dichtemodell 1 beruht auf einer bis in 90 km Tiefe subduzierten osteuropäischen Unterkruste, die massiv eklogitisiert wurde (Dichte 3,55 g /cm3 ). Die Ursache für diese Dichteinterpretation sind geringe Dichten im Bereich der bivergenten geologischen Strukturen der Unterkruste in 24 bis 42 km Tiefe. Diese bivergenten Strukturen sind ebenso die Folgen einer fossilen Subduktion. 2. Im Dichtemodell 2 bildet eine unter dem Orogen liegende osteuropäische Unterkruste eine Krustenverdickung, in deren Folge die Krustenmächtigkeit bis zu 60 km erreicht. Während der Delamination der Unterkruste muß ebenfalls eine Eklogitisierung stattgefunden haben. Das Dichtemodell 2 basiert auf höheren Dichten der bivergenten geologischen Strukturen innerhalb der Unterkruste als im Dichtemodell 1. Eine Gebirgswurzel mit einem deutlichen Kontrast zwischen Kruste und Mantel existiert nicht im Südural. Es wird eine Übergangszone mit Dichtegradienten modelliert. Dieser Dichtegradient kann auf Eklogitisierung des Unterkrustenmaterials beruhen und könnte somit auch die Ursache für das Ausbleiben von seismischen Reflektoren sein. Mit beiden Dichtemodellen werden auch Geoidundulation und Magnetfeld berechnet und deren Wirkung mit der Geoidundulation und dem aeromagnetischen Residualfeld verglichen. Der Vergleich ergibt für die Geoidundulation und das aeromagnetische Residualfeld im langwelligen Anomalienbereich em befriedigendes Ergebnis. Die auf der Interpretation der Bouguer-Schwere beruhenden Dichtemodelle 1 und 2 lassen folgende Schlußfolgerungen zum Aufbau der Lithosphäre im Südural zu: • Die bivergenten Strukturen vor allem im unteren Krustenbereich der Südural-Lithosphäre verursachen langwellige negative Schwereanteile. Zu diesen Strukturen gehören: - die nach Osten einfällende osteuropäische Kruste, - der Falten- und Überschiebungsgürtel im Westural, - die nach Westen einfällenden Strukturen des Ost- und Transurals sowie des kasachischen Terrans. • Eine durch Obduktion in der Oberkruste eingelagerte ozeanische Kruste erzeugt in der Magnitogorsk Zone eine positive kurzwellige Schwereanomalie, die die negativen Schwerewirkungen der darunter befindlichen Strukturen überlagert. • Die negative! Schwerewirkung im Bereich der Voruralsenke ist auf die Schwerewirkung der paläozoischen Sedimentbedeckungen zurückzuführen. Die Untersuchungen zum isostatischen Verhalten der Südural-Lithosphäre zeigen, daß keine klassischen isostatischen Modelle den isostatischen Zustand des Gebirges erklären können, da eine gering ausgeprägte Topographie im Bereich des West- und Zentralurals um 50 bis 100 km nach Westen versetzt ist gegenüber der Krustenverdickung im Bereich der Magnitogorsk und Osturalzone. Mit einem modifizierten Verfahren werden Airy-, Pratt- und Vening-Memesz-Modell miteinander kombiniert. Zunächst wird der lokale isostatische Ausgleich mit Hilfe der Auflastverteilung und -kompensation aus den komplexen Dichtemodellen der Südural-Lithosphäre abgeleitet (Airy- und Pratt-Modell). Um die von den Dichtemodellen nicht mehr kompensierten Masseninhomogenitäten im Zusammenhang mit der Rigidität der Lithosphäre betrachten zu können, werden sie in eine ’’topographische” Auflast zurückgerechnet. Anschließend wird in einem Vening-Meinesz-Modell der regionale isostatische Ausgleich für diese zurückgerechnete ’’topographische” Auflast berechnet. Nur unter Berücksichtigung einer hohen Rigidität der Lithosphäre (3,8 • 1024 Nm) befinden sich beide vorgestellten Dichtemodelle im isostatischen Gleichgewicht. Dies ist offensichtlich auch ein wesentlicher Grund dafür, daß sich die geologischen Strukturen in ihrer vertikalen Position seit der Unteren Trias nicht mehr geändert haben.

Description: I investigate the gravity field in the southern Urals in order to test current hypotheses on the lithospheric structure in this area. My approach to this problem includes the following gravimetric methods: • numerical analyses of the gravity field by means of depth estimations, wavelength filtering and application of the Parker algorithm • two- and three-dimensional density modelling. In contrast to the mainly used density modelling, inverse gravity methods (e.g. wavelength filtering) are not successful at resolving the complex geological structures in this area. For example, I am not able to separate the long-wavelength components of the gravity field from Paleozoic sediments and bi-vergent structures in the lower crust. Estimating the gravity effect of the crust-mantle boundary in the study area, using the Parker-algorithm, is not possible due to its variable density contrast. Data from seismic reflection and refraction surveys and other information relating to the sedimentary cover and Moho depth have been used in the construction of complex three-dimensional density models for the southern Urals’ lithosphere. Incorporating the available details on the structural and physical properties of the southern Urals’ lithosphere in the region between 52° - 65° E and 50° - 55, 5° N, two alternative density models are derived: 1. The density model 1 is based on an East European lower crust which could have been subducted to 90 km depth and eclogitized to a density value of 3.55 g/cm3 . The low density values in the bi-vergent geological structures of the lower crust (from 24 till 42 km depth) overlie a high density body in the upper mantle. These bi-vergent geological structures may be remnants of a fossil subduction setting. 2. In the density model 2, crustal thickening is caused by the East European crust underlying the orogen. The crustal thickness is 60 km. During delamination, the East European lower crust may have been eclogitized. The density model 2 is based on higher density values for the bi-vergent geological structures within the lower crust than in the density model 1. A crustal root, having a significant density contrast between crust and mantle, does not exist in the southern Urals. Using a density gradient, a transition zone is modelled. This density gradient is inferred to be due to eclogitized material in the lower crust. The increasing density could be the cause for the lack of seismic reflectors. The geoid undulation and the magnetic field are calculated for both density models. The results compare satisfactorily with the long- wavelength components of the observed geoid undulation and aeromagnetic residual field. Based on the interpretation of the Bouguer gravity, density models 1 and 2 reveal implications concerning the lithospheric fabric of the southern Urals: • The bi-vergent structures, particularly in the lower crust of the southern Urals’ lithosphere, cause long-wavelength gravity lows. There are following structures: ─ the East European crust dipping to the east, ─ the foreland fold and thrust belt located in the western Urals, ─ the westward-dipping structures of the eastern and Trans-Urals as well as the Kazakhstan terrain. • The obducted oceanic crust in the upper crust of the Magnitogorsk zone causes positive short-wavelength gravity anomalies which are superimposed on the long-wavelength negative gravity effects of the structures located beneath. • Negative gravity anomalies in the area of the Pre-Uralian Foredeep are associated with the gravity effects of the paleozoic sedimentary cover. The investigations into the isostatic behaviour of the southern Urals’ lithosphere show that no classic isostatic model is able to explain the isostatic state of the orogen. The moderate topography in the area of the West and Central Urals is shifted by 50 to 100 km to the west, relative to the largest crustal thickness beneath the Magnitogorsk and East Uralian Zone. A modified method combines the Airy-, the Pratt- and the Vening-Meinsz models. Firstly the local isostatic balance is inferred from the load distribution and the load compensation using the complex density models of the southern Urals’ lithosphere (Airy- and Pratt-model). In order to consider the non-compensated mass inhomogeneities of the density models according to the flexural rigidity of the lithosphere, the non-compensated mass inhomogeneities are recalculated to a ’’topographic” load. Next, the regional isostatic balance as a Vening-Meinesz-model is calculated for the recalculated ’’topographic” load. When allowance is made for the high flexural rigidity of the lithosphere (3.8 • 1024 Nm), both density models presented here are in isostatic balance. Obviously this is the main reason why the geological structures have not changed their vertical position since the Lower Triassic time.

Description: Целью данной работы является проверка существующих моделей строения литосферы Южного Урала по гравиметрическим данным. Для этого исрользуются методы численного анализа гравитационного поля, основанные на: • эмпирической оценке глубин источников аномалий, частотной фильтрации, методе Паркера и • прямом гравитационном моделировании литосферы (двух- и трехмерном). В отличие от обычно используемой технологии гравитационного моделирования, такие методы, как, например, частотная фильтрация не могут быть использованы в данном случае. Невозможно разделить влияние длинноволновых компонент гравитационного поля, обусловленных влиянием палеозойских осадочных отложений и наклоненных слоев нижней коры ц обеих сторон Урала. Также неправомерно использовать алгоритм Паркера для оценки влияния границы между корой мантией вследствие переменного перепада плотности на ней. Включение в анализ доступных опорных данных, основанных на результатах сейсмических методов с использованием отраженных и преломленных волн, данных об осадочном чехле и положении границы Мохо приводит к созданию комплексной плотностной модели литосферы Южного Урала. Использование всех существующих данных о структуре и физических свойствах литосферы Южного Урала в пределах 52° - 65° Е и 50° - 55,5° N позволило построить две эквивалентные плотностные модели: 1. Первая плотностная модель основана на представлении о том, что нижняя кора Восточно-Европейской платформы погружается на глубину 90 кт, при этом происходит эклогитизация её вещества с увеличением плотности до 3,55 g/cm3. Таким образом, наклоненные с обеих сторон блоки нижней коры ответственны за существование высокоплотного тела в верхней мантии. На глубинах 24 - 42 km они характеризуются сравнительно низкими значениями плотностей. Существование этих блоков может быть объяснено также как результат обычной субдукции. 2. Согласно второй модели, утолщение коры объясняется поддвигом под ороген коры Восточно-Европейской платформы. Толщина коры достигает 60 кт. Вследствие деламинации, часть её также должна быть эклогитизирована. Во второй модели значения плотностей наклоненных блоков нижней коры существенно больше, чем в пербой. Копень коры с существенным перепадом плотности под Южным Уралом отсутствует. Вместо этого мы используем при моделировании переходную зону с некоторым градиентом плотности. Этот градиент может быть оценен в предположении об эклогитизации нижней коры. Непрерывное увеличение плотности может быть связано с отсутствием отражающих горизонтов на этих глубинах. Ундуляции геоида и вариации магнитного поля были рассчитаны для обоих моделей и сопоставлены с наблюденными значениями геоида и аэромагнитными измерениями. Оказалось, что удовлетворительное соответствие наблюдается только для длинноволновых аномалий. Таким образом, основанные на интерпретации аномалий Буге плотностые модели 1 и 2 позволяют сформулировать следующие выводы о строении литосферы Южного Урала: • наклоненные в сторону Урала структуры коры, в особенности нижней, ответственны за существование длинноволнового гравитационного минимума. К ним можно отнести: ─ погружающуюся на Восток Восточно-Европейскую платформу, ─ предгорный складчатый пояс Южного Урала, ─ заглубленные в западном напрвлении структуры Восточный зоны, Трансуральской зоны и Казахстанского террэйна. • Фрагмент океанической коры, находящийся вблизи от поверхности в Магнитогорской зоне, вызывает сравнительно коротковолновую положительную аномалию гравитационного поля, которая наложена на длинноволновую отрицательную аномалию, обусловленную структурами, располагающимися ниже. • Отрицательная аномалия в районе Предуральского прогиба может быть объяснена за счёт гравитационного эффекта палеозойского осадочного чехла. Исследование изостатического состояния литосферы Южного Урала показывает, что классические изостатические модели не позволяют объяснить изостазию этого горного массива. Современная топография Западного и Центрального Урала смещена на 50 - 100 кт к западу по отношению к самому глубокому корню коры. Модифицированный метод соединяет изостатические модели Эри, Пратта и Вейнинг-Мейнеса. Первоначально, локальный изостатический баланс оценивается на основании данных о нагрузке и её компенсации в соответствие с построенной комплексной моделью литосферы Южного Урала (модели Эри и Пратта). Для того, чтобы оценить сумму нескомпенсированных плотностных неоднородностей, поддерживаемых за счёт жесткости литосферы, эти неоднородности пересчитываются в ” топографическую” нагрузку. Затем, для этой нагрузки рассчитывается региональная изостатическая компенсация в соответствие со схемой Вейнинг-Мейнеса. Обе рассмотренные плотностные модели могут считаться изостатически сбалансированы только в предположении о высокой эффективной жесткости литосферы (3,8 • 1024 Nm). Очевидно, это является также основной причиной того, что эти геологические структуры не претерпели существенных вертикальных движений со времени позднего Триаса.

Description: thesis

Description: DFG, SUB Göttingen

Description: In dieser Arbeit wird die elektrische Leitfähigkeit von Krustengesteinen während des partiellen Schmelzens mit der Verteilung der Schmelze im Gestein und dem Schmelzanteil als Funktion der Temperatur verglichen. Es wurden impedanzspektroskopische Messungen der elektrischen Leitfähigkeit bei Temperaturen zwischen 600 und 1200°C durchgeführt. Die Messungen erfolgten bei Normaldruck und verschiedenen Sauerstoffugazitäten. Diese Messungen wurden mit den Ergebnissen von Schmelzexperimenten verglichen. Aus den Schmelzexperimenten wurden Informationen zur Verteilung der Schmelze im Gestein und zum Schmelzanteil in Abhängigkeit von der Temperatur gewonnen. Alle Untersuchungen wurden am gleichen Gestein unter den selben experimentellen Bedingungen durchgeführt. Während des partiellen Schmelzens gibt es einen sprunghaften Anstieg der elektrischen Leitfähigkeit um ein bis zwei Größenordnungen. Dieser sprunghafte Anstieg kann mit der Bildung eines vollständigen Netzwerkes an Schmelze erklärt werden, wobei die Schmelzfilme als Leitungspfade wirken. Es konnte festgestellt werden, daß zur Ausbildung der Gleichgewichtsmorphologie während des partiellen Schmelzens sehr lange Zeiten erforderlich sind. Wird die elektrische Leitfähigkeit am Beginn des Schmelzens bei konstanter Temperatur als Funktion der Zeit gemessen, so treten erst nach mehr als 200 h keine Änderungen der Leitfähigkeit mehr auf. Mit Hilfe eines modifizierten Brick-Layer-Modells (MBLM) kann die elektrische Leitfähigkeit als Funktion des Schmelzanteils berechnet werden. Dieses Modell gilt für den Fall, daß die Schmelze ein vollständiges Netzwerk entlang von Komgrenzen bildet. Mit dem Modell und den Ergebnissen aus den Schmelzexperimenten konnten theoretische Kurven der elektrischen Leitfähigkeit berechnet werden, die mit den gemessenen Kurven der elektrischen Leitfähigkeit verglichen werden können. Obwohl bereits bei Beginn des Schmelzens die Komgrenzen von Schmelze benetzt werden, wird eine vollständige Vernetzung der Schmelze erst bei höheren Temperaturen und größeren Schmelzanteilen erreicht. Dies kann durch Abweichungen vom morphologischen Gleichgewicht in den partiell geschmolzenen Proben gedeutet werden. Die Ergebnisse stützen die Annahme, daß eine Zone hoher Leitfähigkeit unter den zentralen Anden durch große Mengen an Magma verursacht wird. Mit dem MBLM konnte abgeschätzt werden, daß ca. 20% Schmelze notwendig sind, um die im Gelände mit der Magnetotellurik gemessenen Leitfähigkeiten zu erklären.

Description: The purpose of the present work is to achieve a better understanding of rocks during partial melting. Therefore the electrical conductivity of crustal rocks during partial melting was compared to the distribution of melt in the rock sample and to the melt fraction as a function of temperature. Impedance spectroscopic investigations of the electrical conductivity were conducted at temperatures between 600 and 1200°C, normal pressure and at different oxygen fugacities. These measurements were compared with the results of melting experiments. The melting experiments were performed under the same experimental conditions and the same rock sample was used. The rock samples were equilibrated at the desired temperatures, quenched and investigated using thin sections. These melting experiments provide information about the distribution of melt in the rock sample and yield the melt fraction as a function of temperature. A strong increase in electrical conductivity of about two orders of magnitude is observed during partial melting. This increase could be explained by assuming the formation of an interconnected network of melt. The charge transport follows the network forming melt films at the grain boundaries. It could be established that the formation of steady state of a partially molten rock requires a long time. The conductivity was measured as a function of time at temperatures little above the solidus. Under these conditions constant conductivity values are found after 200 h.

Description: thesis

Description: DFG, SUB Göttingen

Description: In der vorliegenden Arbeit wurden weitwinkelseismische Daten aus dem Land-See-kombinierten Projekt CINCA 95 ausgewertet. Die Daten wurden auf drei EW-streichenden Profilen bei 21°S, 22°S und 23.25°S registriert. Die Linien umfassen das Gebiet zwischen Peru-Chile- Trench und Westkordillere, also den Forearc-Bereich der nord-chilenischen erosiven Subduktionszone. Der Hauptteil der Daten besteht aus einer Anregerlinie auf See und einer Empfängerlinie auf Land. Ergänzend wurden Sprengungen am westlichen und östlichen Ende der Empfängerlinie durchgeführt. Aus dem Datensatz wurden in der ersten Bearbeitungsstufe mit den Methoden der Refraktionsseismik drei vp-Geschwindigkeitsmodelle erstellt, die ihre Aussage im äußeren Forearc, also dem Gebiet zwischen Rinne und Küste haben. Dabei kann folgendes festgehalten werden: • Sie zeigen die Nazca-Platte mit einem sich versteilenden Subduktionswinkel von 10° direkt nach der Subduktion auf 23-25° vor der Küste Nord-Chiles. Die subduzierte Platte zeigt, was ihren Abtauchwinkel anbelangt, keine NS-Variation. • In der Geschwindigkeitszunahme an der Oberkante der subduzierten Platte zeigt sich durchaus eine Differenz in NS-Richtung: im Norden nimmt die Geschwindigkeit von 4.5 km/s direkt nach dem Untertauchen unter die südamerikanische Platte auf 6.5 km/s auf einer Distanz von 70 km zu, während sie im Süden dazu nur noch 30 km benötigt. Diese Struktur wird als Lage erodierten und subduzierten Gesteinsmaterials zwischen den Platten interpretiert. • Bereits im westlichsten Krustenblock, also dem Teil, der heute an der Rinne zur Erosion ansteht, zeigt sich eine hohe vp-Geschwindigkeit von 5.8 km/s. Sie nimmt zur Küste hin noch auf 6.0 km/s direkt an der Oberfläche zu. Die Kruste ist im Norden deutlich differenziert. Auch dies nimmt nach Süden hin ab und so sind auf dem südlichsten Profil keine Krustenstrukturen mehr in Tiefen unterhalb von 11 km modelliert worden. • Auf den Profilen bei 21°S und 22°S wurde eine platten-parallele Diskontinuität, die die Oberkante zu einer Zone erhöhter Geschwindigkeit bildet (7.1 km/s), modelliert. Sie hat eine Neigung von 15°. Diese Diskontinuität wird im Zusammenhang mit bereits existierenden Modellen für den landseitigen Forearc als Paläo-Moho interpretiert. Das darunterliegende Gestein mit einer Geschwindigkeit von 7.0 km/s wird entsprechend als serpentinisierter oberster Mantel verstanden. • Unterhalb schließt sich eine Zone stark erniedrigter Geschwindigkeit an, die zwischen der abtauchenden Platte und der Oberplatte einen Keil bildet. Sie hat eine Geschwindigkeit von 6.4 km/s und wird als Anlagerungsort des erodierten und subduzierten Gesteins interpretiert. • Die östlich der Küste in den vp -Schnitten als Moho modellierte Diskontinuität wurde durch die vorherigen Arbeiten bestätigt und als die heutige, geophysikalische Moho interpretiert. Zur Erosion am Kontinentalrand kann folgendes festgestellt werden: • Eine Varianz im Ablauf der Erosion wird nahegelegt, kann aber nicht bewiesen werden. • Nach der hier vorgelegten Interpretation findet die Erosion in der Hauptsache frontal statt. In einem zweiten Bearbeitungsschritt wurde die seeseitige Spurdichte von 160 m zu einer CMP-Stapelung der Weitwinkeldaten ausgenutzt. Bei der Bearbeitung der Daten wurden ausschließlich Tools der Standard-Reflexionsseismik verwendet. Dabei sollte ein Strukturabbild ähnlich dem aus der Standard-Reflexionsseismik bekannten entstehen. Diese Art des Prozessierens von Daten wurde bisher auf einen Datensatz mit so großer Apertur noch nicht angewandt. Es konnte gezeigt werden, daß mit einem unkonventionellen Ansatz durchaus Zusatzinformationen aus den Daten gewonnen werden können: die Oberkante der Niedergeschwindigkeitszone sowie die ozeanische Moho konnten auf dem nördlichsten Profil bestätigt werden, die platten-parallel liegende Struktur sowie die ozeanische Moho konnten auf dem mittleren Profil bestätigt werden und auf dem südlichsten zeigte sich, dass die Unterkante der ozeanischen Kruste weiter nach Osten hin detektiert werden konnte, als durch die refraktionsseismische Bearbeitung möglich war.

Description: In this thesis wide-angle seismic data from the on-shore/off-shore CINCA 95project were analysed. The data set was recorded on three E-W profiles at 21°S, 22°S und 23.25°S. These profiles cover the region between the Peru-Chile trench and the Western Cordillera, i.e. the forearc of the north-Chile subduction zone. The north-Chile margin is an example of a non-accreting, eroding continental margin. The data set mainly consists of three off-shore lines of airgun shots and on-shore receiver lines. Additional chemical shots were performed at the ends of the receiver lines. A vp-velocity model was derived for each profile using seismic refraction interpretation methods. The models are well constrained in the outer forearc, i.e. the region between the trench and the coast line. Major results are summarized as follows. • The models show the Nazca plate subducting at an angle of 10°near the trench, increasing to 23-25° near the coast. There is no N-S variation in the behaviour of the subduction angle. • However, a N-S variation in the lateral velocity increase at the top of the subducted plate can be derived. On the northernmost profile the velocity increases from 4.5 km/s at the trench to 6.5 km/s, 70 km E of the trench, whereas on the southernmost profile, this velocity increase occurs over a distance of only 30 km. This structure is interpreted as a layer of eroded and subducted material. • In the westernmost portion of the continental crust, where erosion occurs, a high vp-velocity of 5.8 km/s is found. It increases to 6.0 km/s towards the coast. In the north the crust is more differentiated. This differentiation weakens southwards and on the southernmost profile no additional crustal structure in the upper plate was modelled beneath 11 km depth. • On the profiles at 21°S and 22°S a discontinuity parallel to the subducting plate was modelled, representing the top of a high velocity (7.1 km/s) zone. The discontinuity dips at an angle of 15°. On the basis of an existing model for the on-shore part of the forearc this discontinuity is interpreted as a paleo-Moho. The high velocity zone is thus interpreted to be serpentinized uppermost mantle. • The high velocity zone is followed by a low velocity zone, which forms a wedge between the subducting plate and the overriding plate. It has a velocity of 6.4 km/s and is interpreted as a depot for the eroded and subducted material. • The Moho east of the coast is confirmed by former studies in the area and interpreted as the recent geophysical Moho of the South American continent. Concerning erosion at a continental margin the following can be noted. • An episodical character of erosion is implied by the models but cannot be proven. • The interpretation presented in this thesis implies frontal erosion along the Peru-Chile trench in this region. For further interpretation the dense shot line off-shore (160 m shot interval) was utilized for a CMP-stack of the wide-angle data. This was achieved with standards tools used in the processing of near- vertical incidence reflection data. The goal was to achieve a structural image similar to that which can be obtained from standard near-vertical incidence data. So far, this kind of processing has not been used for data with such a large aperture. Using this non-conventional approach some further information could be drawn from the data set. On the profile at 21°S the upper boundary of the low velocity wedge between the plates could be confirmed as well as the oceanic Moho. On the profile at 22°S the plate parallel structure as well as the oceanic Moho could be confirmed. On the profile at 23.25°S, with the help of the CMP processing, the subducting oceanic Moho could be traced further east than in the velocity model which was derived using standard seismic refraction interpretation methods.

Description: thesis

Description: DFG, SUB Göttingen

Description: In der vorliegenden Arbeit werden die Ergebnisse des refraktionsseismischen Experimentes PISCO'94 [Proyecto de Investigaciones Sismológicas de la Cordillera Occidental 1994] präsentiert, das 1994 vom Sonderforschungsbereich 267 "Deformationsprozesse in den Anden" im Andensegment 21-24°SZ 67-70°W durchgeführt wurde. Das Experiment knüpfte an ein in den 80er-Jahren vermessenes Netz refraktionsseismischer Profile an. Das Ziel der neuen Messungen war die Bestimmung der Lithosphärenstruktur im Bereich des magmatischen Bogens und seiner Nachbargebiete. Zu diesem Zweck wurde ein 300 km langes N-S streichendes Profil mit 4 Schußpunkten am westlichen Rand des magmatischen Bogens, der Westkordillere, vermessen. Weiterhin wurde ein neues W-E streichendes Profil, welches bei 23°30'S von der Krustenkordillere bis zur Westkordillere reicht, beobachtet. Bereits früher vermessene Profile in der Präkordillere wurden durch Gegenbeobachtungen bzw. weitere Schußpunkte ergänzt. Aus den abgeleiteten ein- und zweidimensionalen Modellierungen wurde ein dreidimensionales Übersichtsmodell für das Untersuchungsgebiet erstellt. Im Forearc-Bereich bestätigt es im wesentlichen die Ergebnisse früherer refraktionsseismischer Untersuchugen: Die Durchschnittsgeschwindigkeit der Kruste nimmt von 6.6 km/s unter der Küstenkordillere auf 6.2 km/s unter der Präkordillere ab. Die Moho der subduzierten ozeanischen Platte wird in ca. 40 km Tiefe unter der Küstenkordillere beobachtet. Die ursprüngliche kontinentale Moho wurde stark überarbeitet und ist seismisch nicht eindeutig zu bestimmen. Eine Diskontinuität in 10-12 km Tiefe unter der Küstenkordillere, die unter der Präkordillere bei ca. 20 km Tiefe liegt, wird als Grenze zwischen Ober- und Unterkruste angesehen. In 20 km Tiefe unter der Küstenkordillere befindet sich ein Bereich hoher Geschwindigkeit, der als Basis der kontinentalen Kruste interpretiert wird. Dieser Bereich läßt sich in die Präkordillere hinein verfolgen, wo er bei ca. 35-40 km Tiefe liegt. Der Bereich darunter weist bis in ca. 60 km Tiefe Krustengeschwindigkeiten (vp 〈 7.5 km/s) auf. Dies könnte durch hydratisiertes Material des peridotitischen Mantelkeils erklärt werden. Die Natur einer markanten, in 60-70 km Tiefe liegenden Diskontinuität, an der die Geschwindigkeit regional auf über 7.9 km/s zunimmt und die als "Moho" bezeichnet wird, ist unklar. Vermutlich handelt es sich um die Oberkante einer Übergangszone zwischen Asthenosphäre und Lithosphäre. Die für die magmatischen Bogen aus seismologischen und gravimetrischen Untersuchungen abgeleitete Krustendicke von ca. 70 km wird durch die neuen refraktionsseismischen Daten nicht bestätigt Bis in eine Tiefe von 60 km wurden keine Mohoreflexionen beobachtet. Jedoch können die seismischen Diskontinuitäten der Präkordillere, insbesondere die Grenze zwischen Ober- und Unterkruste in ca. 20 km Tiefe, in die Westkordillere hinein verfolgt werden. Bis in 45-60 km Tiefe überschreitet die Geschwindigkeit der beobachteten seismischen Diskontinuitäten 7.5 km/s nicht Die Durchschnittsgeschwindigkeit ist mit 5.9-6.0 km/s gering. Im Streichen des magmatischen Bogens ändern sich Krustenstruktur und Geschwindigkeitsverteilung im Übergangsbereich von steiler zu flacher Subduktion. Im Norden des Untersuchungsgebietes wird eine stark strukturierte Kruste mit z.T. extremen Niedriggeschwindigkeitszonen und starker Dämpfung seismischer Wellen beobachtet. Eine plausible Erklärung dafür ist die Existenz partieller Schmelzen in der Kruste. Als Quelle des extremen Schwerehochs, das von Calama südöstlich streicht, wurde Material mit Geschwindigkeiten von 〉 6.4 km/s in 2-10 sowie 15-20 km Tiefe bestimmt. Im südlichen Teil des Untersuchungsgebietes dagegen ist die Kruste des magmatischen Bogens relativ unstrukturiert und weist keine erhöhte Dämpfung mehr auf. Als Ursache für diese Differenzen sowie das Zurückversetzen des magmatischen Bogens südlich von 23°S wird eine (prä-neogene) dünne und harte Kruste im südlichen Teil des Untersuchungsgebietes angesehen.

Description: The object of this paper is the presentation of the results of the seismic refraction experiment PISCO'94 [Proyecto de Investigaciones Sismológicas de la Cordillera Occidental 1994], carried out in 1994 by the Collaborative Research Center 267 "Deformation Processes in the Andes" in the Andean segment 21-24°S/ 67-70°N. The experiment continued a net of seismic refraction profiles measured in the 80s. The aim of the new measurements was to investigate the crustal structure of the magmatic arc region and its vicinity. A N-S-running profile, 300 km long with 4 shotpoints, was measured along the western margin of the magmatic arc, the Western Cordillera. Furthermore a new W-E profile at 23°30'S was observed, which reaches from the Coastal Cordillera to the Western Cordillera. Complementary shotpoints and reversed observations to earlier profiles in the Precordillera were realized. From the derived 1- and 2- dimensional models a generalized three dimensional model was produced. In the forearc region it confirms essentially the results of earlier seismic refraction studies. The average P-wave velocity of the crust decreases from 6.6 km/s beneath the Coastal Cordillera to 6.2 km/s beneath the Precordillera. The Moho of the downgoing oceanic plate is observed at about 40 km depth beneath the Coastal Cordillera. The original continental Moho was strongly modified and cannot be clearly determined by seismic data. The discontinuity in 10-12 km depth beneath the Coastal Cordillera, downgoing to about 20 km depth beneath the Precordillera, is interpreted as boundary between upper and lower crust. A high velocity zone at about 20 km depth beneath the Coastal Cordillera is interpreted as base of the continental crust This zone can be continued to a depth of 35-40 km beneath the Precordillera. The region below is characterized down to about 60 km depth by crustal velocities (vp 〈7.5 km/s). One explanation could be hydrated material of the peridotitic mantle wedge. In 60-70 km depth the velocity increases regionally to values above 7.9 km/s. The nature of this prominent discontinuity, called "Moho", is unclear. It indicates probably the top of a lithosphere/asthenosphere transition zone. A crustal thickness of 70 km for the magmatic arc, derived from seismological and gravimetric investigations, is not confirmed by the new seismic refraction datas. There are no Moho observations down to 60 km, but the seismic discontinuities of the Precordillera can be traced into the Western Cordillera. A discontinuitiy at about 20 km depth is interpreted as top of the lower crust, which reaches to about 40 km depth. Down to 45-60 km the velocities of the observed seismic discontinuities do not exceed 7.5 km/s. The average velocity is low (5.9-6.0 km/s). Along the strike of the magmatic arc crustal structure and velocity are changing in the transition zone from steep to flat subduction. In the north of the area under investigation a strongly structured crust with zones of extremely low velocities and strong attenuation of seismic waves is observed. This could be explained by partial melts in the crust. As source of the gravity high SE of Calama material with velocities 〉 6.4 km/s in 2-10 and 15-20 km depth was determined. In contrary, the crust of the southern part of the area is relatively uniform without high attenuation. These differences as well as the rebound of the magmatic arc south of 23°S could be caused by a pre-neogene thin and hard crust in the southern part of the area under investigation.

Description: thesis

Description: DFG, SUB Göttingen

Description: Eine für die Zentralanden erstellte geowissenschaftliche Datenbank erfaßt die umfangreichen Forschungsergebnisse des SFB 267 wie auch externe Datensätze in Form eines zentral verwalteten Dateisystems. Ein Datenkatalog, welcher seit Mitte 1996 auch im World Wide Web (WWW, URL: http://userpage.fu-berlin.de/~data) verfügbar ist, bietet einen Überblick über die aktuelle Datenbasis und erleichtert den allgemeinen Zugriff auf die einzelnen Datensätze. Vorgestellt werden ferner verschiedene Methoden zur Analyse raumbezogener Daten aus Bereichen der Statistik (streifen-gemittelte Topographieprofile, Standardabweichung, Korrelationen, Fraktale) , Bildverarbeitung (direktionale Richtungsfilter) , Geomorphometrie (Neigung und Vertikalwölbung) und Numerik (Horizontalgradienten, FFT). Die Anwendung dieser Verfahren bezüglich einer regionalen Strukturierung wird am Beispiel der für die Zentralanden flächendeckend vorliegenden Topographie und Schwerefelder getestet und mit den aus langjährigen Untersuchungen in den Anden zur Verfügung stehenden Forschungsergebnissen verglichen. Endogene und exogene Prozesse spiegeln sich in den topographischen Erscheinungsformen wider. Insbesondere mit Hilfe der geomorphometrischen Analyse und der FFT-Filterung lassen sich für die morphostrukturellen Einheiten typische Merkmale sehr gut herausarbeiten. Auch mit einer einfachen Statistik (streifen-gemittelte Topographieprofile, Berechnung von Minimum, Mittelwert, Maximum und Standardabweichung einzelner Untersuchungsgebiete) können für die ausgewählten Gebiete charakteristische Unterschiede ermittelt werden. Im Hinblick auf die Topographie unterscheidet sich bei allen Untersuchungsmethoden am Westrand der Zentralanden der hyperaride Norden im Bereich der Atacama Wüste deutlich vom weniger trockenen Süden. Studien zur fraktalen Geometrie der Geländeoberfläche zeigen, daß sich die andine Topographie skaleninvariant verhält und tektonische wie klimatische Einflüsse sich in der geometrischen Variation der Geländeoberfläche widerspiegeln. Morphologisch zergliederte Zonen, wie das durch starke Verkürzungen gekennzeichnete Subandin oder der Südwestrand der Zentralanden, weisen im Vergleich zur fraktalen Dimension des Gesamtgebietes eine höhere fraktale Dimension auf, während das hyperaride Gebiet im Bereich der Atacama Wüste durch einen vergleichsweise niedrigen Wert gekennzeichnet ist. Mit der direktionalen Lineamentverstärkung, z.T. aber auch mit anderen Methoden wie etwa der Darstellung der Neigungsrichtung, werden die wichtigsten Störungszonen (z.B. das Atacama- oder das Präkordilleren-Störungssystem) weitestgehend erfaßt. Auf dem Pazifik bilden sich - sowohl bei der geomorphometrischen Untersuchung wie auch bei der Lineamenterkennung - entlang der Peru-Chile-Tiefseerinne die im CINCA Experiment beobachteten morphologisch-tektonischen Charakteristika (Reichert et al., 1997), wie z.B. die Blocktektonik innerhalb einer ca. 50 km breiten, westlich der Tiefseerinne verlaufenden Zone, ab. Gravimetrische Lineamente, die aus den maximalen Horizontalgradienten im Schwerefeld abgeleitet werden, weisen auf abrupte, laterale Dichteänderungen hin, sind jedoch in den Zentralanden nur in Einzelfällen mit den an der Oberfläche anstehenden lithologischen Einheiten oder Störungszonen korrelierbar. Die zweidimensionale Korrelation im moving window Verfahren hat sich als nützliches Instrument zur Analyse von sich gegenseitig beeinflussenden, raumverteilten Datensätzen erwiesen. Die Wahl der Fenstergröße erlaubt, zwischen einer regionalen und einer eher lokalen Analyse zu differenzieren; mit Hilfe einer geeigneten Visualisierung lassen sich die Beziehungen zwischen den Variablen schnell erfassen. Die räumlich differenzierte, quantitative Analyse der Korrelation zwischen Schwerefeld und Topographie in den Zentralanden zeigt, daß nicht nur die Bouguer-Anomalie sondern auch das isostatische Restfeld stark negativ mit der Geländehöhe korreliert.

Description: A geoscientific database was established for the area of the Central Andes, which includes information about past and current research conducted within the frame of the ’SFB 267 - Deformation Processes in the Andes’ (Berlin, Potsdam) as well as data from external sources. The catalogue, which since 1996 is also available on the Web (WWW, URL: http://userpage.fu-berlin.de/~data/ provides a general view of the database architecture and contents, and facilitates common and easy access to singular sets of data. In addition, various methods for the analysis of spatial data from disciplines such as statistics (swath-averaged topographic profiles, standard deviation, correlation, fractals), digital image processing (directional filtering), geomorphometry (inclination and vertical doming) and numerical methods (horizontal gradient, Fast Fourier Transformation [FFT]) are presented in this thesis. Application of these techniques for structural analysis on a regional scale was tested for the Central Andes, for which topographic coverage is complete, as well as for associated gravity fields, and compared with information generated during long term investigations in the Andes. Endogenous and exogenous processes are mirrored in topographic forms and landscapes. Characteristic features of morphostructural units are particularly well characterized when employing geomorphometric analyses and FFT filter techniques. Even simple statistical methods (swath-averaged topographic profiles, calculation of minimum, mean, maximum and standard deviation) for selected areas allow recognition of specific features and characteristic variations. As far as topography is concerned, distinct differences between the hyperarid northern portion of the western Central Andes (Atacama Desert) and lesser arid southern sectors are evident for all methods employed. Studies of the surface fractal geometry show an invariance of scale of the Andean topography. Structural and climatic influences are mirrored in geometric variations of the landscape. Morphologically intensely structured portions, like the Subandean ranges, characterized by pronounced crustal shortening, or the southwestern margin of the Central Andes are marked by a higher fractal dimension, when compared to the entire area studied, whereas the hyperarid portion of the Atacama Desert is characterized by a lower fractal dimension. Directional lineament enhancement, but other methods such as areal inclination detection as well, allow to accentuate major zones of structural weakness, such as the Atacama and Precordilleran fault systems. Both, geomorphometric investigations and lineament recognition, enhance morpho-structural features such as block-faulting of oceanic lithosphere along the Peru-Chile deep sea trench in a 50 km wide sector to the west of the subduction zone, and verify results gathered during the CINCA experiment (Reichert et al., 1997). Gravimetric lineaments, derived from maximum horizontal gradients in gravity fields of the Central Andes, point to abrupt lateral density variations, but can only occasionally be correlated with outcropping lithologic units or fault zones. The two dimensional correlation using the moving average technique has proven to be useful for analysis of mutually influencing and spatially arranged data sets. Selection of a proper window size allows to differentiate between analyses on regional resp. local scales. Employing proper visualization, relations between variables are easily detectable. A spatially differentiated and quantitative analysis of the relationship between gravity field and topography in the Central Andes shows a strong negative correlation of altitude with Bouguer anomaly and isostatic residual field.

Description: http://www.cms.fu-berlin.de/sfb/sfb267/results/data_catalogue/index.html

Description: thesis

Description: DFG, SUB Göttingen