TY - THES A1 - Rodriguez Piceda, Constanza T1 - Thermomechanical state of the southern Central Andes T1 - Thermomechanischer Zustand der südlichen Zentral Anden BT - implications for active deformation patterns in the transition from flat to steep subduction BT - Implikationen für aktive Deformationsmuster beim Übergang von flacher zu steiler Subduktion N2 - The Andes are a ~7000 km long N-S trending mountain range developed along the South American western continental margin. Driven by the subduction of the oceanic Nazca plate beneath the continental South American plate, the formation of the northern and central parts of the orogen is a type case for a non-collisional orogeny. In the southern Central Andes (SCA, 29°S-39°S), the oceanic plate changes the subduction angle between 33°S and 35°S from almost horizontal (< 5° dip) in the north to a steeper angle (~30° dip) in the south. This sector of the Andes also displays remarkable along- and across- strike variations of the tectonic deformation patterns. These include a systematic decrease of topographic elevation, of crustal shortening and foreland and orogenic width, as well as an alternation of the foreland deformation style between thick-skinned and thin-skinned recorded along- and across the strike of the subduction zone. Moreover, the SCA are a very seismically active region. The continental plate is characterized by a relatively shallow seismicity (< 30 km depth) which is mainly focussed at the transition from the orogen to the lowland areas of the foreland and the forearc; in contrast, deeper seismicity occurs below the interiors of the northern foreland. Additionally, frequent seismicity is also recorded in the shallow parts of the oceanic plate and in a sector of the flat slab segment between 31°S and 33°S. The observed spatial heterogeneity in tectonic and seismic deformation in the SCA has been attributed to multiple causes, including variations in sediment thickness, the presence of inherited structures and changes in the subduction angle of the oceanic slab. However, there is no study that inquired the relationship between the long-term rheological configuration of the SCA and the spatial deformation patterns. Moreover, the effects of the density and thickness configuration of the continental plate and of variations in the slab dip angle in the rheological state of the lithosphere have been not thoroughly investigated yet. Since rheology depends on composition, pressure and temperature, a detailed characterization of the compositional, structural and thermal fields of the lithosphere is needed. Therefore, by using multiple geophysical approaches and data sources, I constructed the following 3D models of the SCA lithosphere: (i) a seismically-constrained structural and density model that was tested against the gravity field; (ii) a thermal model integrating the conversion of mantle shear-wave velocities to temperature with steady-state conductive calculations in the uppermost lithosphere (< 50 km depth), validated by temperature and heat-flow measurements; and (iii) a rheological model of the long-term lithospheric strength using as input the previously-generated models. The results of this dissertation indicate that the present-day thermal and rheological fields of the SCA are controlled by different mechanisms at different depths. At shallow depths (< 50 km), the thermomechanical field is modulated by the heterogeneous composition of the continental lithosphere. The overprint of the oceanic slab is detectable where the oceanic plate is shallow (< 85 km depth) and the radiogenic crust is thin, resulting in overall lower temperatures and higher strength compared to regions where the slab is steep and the radiogenic crust is thick. At depths > 50 km, largest temperatures variations occur where the descending slab is detected, which implies that the deep thermal field is mainly affected by the slab dip geometry. The outcomes of this thesis suggests that long-term thermomechanical state of the lithosphere influences the spatial distribution of seismic deformation. Most of the seismicity within the continental plate occurs above the modelled transition from brittle to ductile conditions. Additionally, there is a spatial correlation between the location of these events and the transition from the mechanically strong domains of the forearc and foreland to the weak domain of the orogen. In contrast, seismicity within the oceanic plate is also detected where long-term ductile conditions are expected. I therefore analysed the possible influence of additional mechanisms triggering these earthquakes, including the compaction of sediments in the subduction interface and dehydration reactions in the slab. To that aim, I carried out a qualitative analysis of the state of hydration in the mantle using the ratio between compressional- and shear-wave velocity (vp/vs ratio) from a previous seismic tomography. The results from this analysis indicate that the majority of the seismicity spatially correlates with hydrated areas of the slab and overlying continental mantle, with the exception of the cluster within the flat slab segment. In this region, earthquakes are likely triggered by flexural processes where the slab changes from a flat to a steep subduction angle. First-order variations in the observed tectonic patterns also seem to be influenced by the thermomechanical configuration of the lithosphere. The mechanically strong domains of the forearc and foreland, due to their resistance to deformation, display smaller amounts of shortening than the relatively weak orogenic domain. In addition, the structural and thermomechanical characteristics modelled in this dissertation confirm previous analyses from geodynamic models pointing to the control of the observed heterogeneities in the orogen and foreland deformation style. These characteristics include the lithospheric and crustal thickness, the presence of weak sediments and the variations in gravitational potential energy. Specific conditions occur in the cold and strong northern foreland, which is characterized by active seismicity and thick-skinned structures, although the modelled crustal strength exceeds the typical values of externally-applied tectonic stresses. The additional mechanisms that could explain the strain localization in a region that should resist deformation are: (i) increased tectonic forces coming from the steepening of the slab and (ii) enhanced weakening along inherited structures from pre-Andean deformation events. Finally, the thermomechanical conditions of this sector of the foreland could be a key factor influencing the preservation of the flat subduction angle at these latitudes of the SCA. N2 - Die Anden sind eine ~7000 km lange N-S-verlaufende Hochgebirgskette, die entlang des westlichen südamerikanischen Kontinentalrandes entstanden ist. Aufgrund der Subduktion der ozeanischen Nazca-Platte unter die kontinentale südamerikanische Platte ist die Bildung des nördlichen und zentralen Teils des Gebirges typisch für eine nicht-kollisionale Orogenese. In den südlichen Zentralanden (SZA, 29-39° S) verändert sich der Subduktionswinkel der ozeanischen Platte zwischen 33 ° S und 35 ° S von fast horizontal (< 5° Einfallen) im Norden zu einem steileren Winkel (~ 30 ° Einfallen) im Süden. Begleitet wird dieser Trend von systematischen, Süd-gerichteten Abnahmen der topographischen Erhebung, der Krusteneinengung und der Vorland- und Orogenbreite, sowie von Variationen im Deformationsstil des Vorlandes, wo die Einengung des Deckgebirges in unterschiedlichem Maße von einer entsprechenden Deformation des Grundgebirges begleitet wird. . Darüber hinaus sind die SZA eine seismisch sehr aktive Region. Die Kontinentalplatte zeichnet sich durch eine relativ flache Seismizität (< 30 km Tiefe) aus, die sich hauptsächlich auf die Übergänge vom Orogen zu den Vorlandbereichen konzentriert; im Gegensatz dazu tritt tiefere Seismizität in den zentralen Bereichen des nördlichen Vorlandes auf. Darüber hinaus ist häufig auftretende Seismizität auch in den flachen Teilen der ozeanischen Platte und im Plattensegment mit flach einfallender Subduktion zwischen 31 ° S und 33 ° S festzustellen. Die beobachtete räumliche Heterogenität der tektonischen und seismischen Deformation in den SZA wurde auf mehrere Ursachen zurückgeführt, darunter Schwankungen der Sedimentmächtigkeit, das Vorhandensein vererbter Strukturen und Veränderungen des Subduktionswinkels der ozeanischen Platte. Es gibt jedoch bislang keine Studie, die den Zusammenhang zwischen der langfristigen rheologischen Konfiguration der SZA und den räumlichen Deformationsmustern untersucht hat. Darüber hinaus wurden die Auswirkungen der Dichte- und Mächtigkeitsvariationen in der kontinentalen Oberplatte und der verschiedenen Subduktionswinkel auf den rheologischen Zustand der Lithosphäre noch nicht grundlegend untersucht. Da die Rheologie von der Gesteinsart, dem Druck und der Temperatur abhängt, ist eine detaillierte Charakterisierung der Zusammensetzung, Struktur und des thermischen Feldes der Lithosphäre erforderlich. Daher habe ich unter Verwendung kombinierter Modellierungsansätze und geophysikalischer Daten die folgenden 3D Modelle für die Lithosphäre der SZA konstruiert: (i) ein auf seismischen Daten basierendes Struktur- und Dichtemodell, das anhand des beobachteten Schwerefeldes validiert wurde; (ii) ein thermisches Modell, das die Umwandlung von Mantelscherwellengeschwindigkeiten in Temperaturen mit Berechnungen des konduktiven Wärmetransports für stationäre Bedingungen in der obersten Lithosphäre (<50 km Tiefe) integriert und durch Temperatur- und Wärmeflussmessungen validiert wurde; und (iii) ein rheologisches Modell der langfristig bedingten Lithosphärenfestigkeit, das auf den zuvor erzeugten Modellen gründet. Die Ergebnisse dieser Dissertation zeigen, dass die thermischen und rheologischen Bedingungen in den heutigen SZA durch verschiedene Mechanismen in unterschiedlichen Tiefen gesteuert werden. In flachen Tiefen (< 50 km) wird das thermomechanische Feld durch die heterogene Zusammensetzung der kontinentalen Lithosphäre differenziert. Eine Überprägung durch die ozeanische Platte ist dort nachweisbar, wo die ozeanische Platte flach (< 85 km tief) und die radiogene Kruste dünn ist, was insgesamt zu niedrigeren Temperaturen und einer höheren Festigkeit im Vergleich zu Bereichen führt, in denen die Platte steil einfällt und die radiogene Kruste dick ist. In Tiefen > 50 km treten die größten Temperaturschwankungen dort auf, wo die subduzierten Platte nachgewiesen wurde, was bedeutet, dass das tiefe thermische Feld den Subduktionswinkel gesteuert wird. Die Ergebnisse dieser Doktorarbeit legen nahe, dass der langfristige thermomechanische Zustand der Lithosphäre die räumliche Verteilung rezenter Seismizität beeinflusst. Der größte Anteil innerhalb der Kontinentalplatte registrierter Erdbebentätigkeit tritt oberhalb des modellierten Übergangs von spröden zu duktilen Bedingungen auf. Außerdem besteht eine räumliche Korrelation zwischen Erdbebenclustern und den Übergängen von den mechanisch rigideren Vorlandbereichen (Forearc und Foreland) zum mechanisch schwächeren Orogen. Demgegenüber wird vermehrte Seismizität innerhalb der ozeanischen Platte auch dort nachgewiesen, wo entsprechend der Modellierung langfristig duktile Bedingungen erwartet werden. Ich habe daher den möglichen Einfluss zusätzlicher Mechanismen untersucht, die ein Auslösen dieser Erdbeben begünstigen könnten, darunter die Kompaktion von Sedimenten an der Subduktionsgrenzfläche und Dehydrationsreaktionen innerhalb der Platte. Dazu habe ich eine qualitative Analyse des Hydratationszustandes des Mantels unter Verwendung des Verhältnisses zwischen Kompressions- und Scherwellengeschwindigkeit (Vp/Vs-Verhältnis aus einemseismischen Tomographiemodell) durchgeführt. Die Ergebnisse dieser Analyse zeigen, dass der Großteil der Seismizität räumlich mit hydratisierten Bereichen in der subduzierten Platte und im darüber liegenden kontinentalen Mantel korreliert, mit Ausnahme eines Erdbebenclusters, das innerhalb des flachen Plattensegments auftritt. In diesem Bereich wechselt die subduzierte Platte von einem flachen in einen steilen Subduktionswinkel und Erdbeben werden wahrscheinlich durch Biegevorgänge in der Platte ausgelöst. Auch die wichtigsten Variationen in den beobachteten tektonischen Mustern scheinen durch die thermomechanische Konfiguration der Lithosphäre beeinflusst zu sein. Die mechanisch starken Bereiche von Forearc und Foreland zeigen aufgrund ihrer Verformungsbeständigkeit geringere Verkürzungsraten als der relativ schwache Bereich des Orogens. Darüber hinaus bestätigen die in dieser Dissertation modellierten strukturellen und thermomechanischen Eigenschaften der Lithosphäre auch frühere Analysen geodynamischer Simulationen, denen zufolge der Deformationsstil im Orogen- und Vorlandbereich jeweils von Variationen in der Lithosphären- und Krustendicke, im Vorhandensein schwacher Sedimente und in der gravitativen potentiellen Energie kontrolliert wird. Eine Sonderstellung nimmt der nordöstliche Vorlandbereich der SZA ein, wo eine verstärkte Seismizität und eine das Deck-und Grundgebirge erfassende Deformation zu beobachten sind, obwohl die modellierte Krustenfestigkeit dort Werte übersteigt, die für die in diesem Gebiet anzunehmenden tektonischen Spannungen typisch wären. . Mechanismen zur Lokalisierung verstärkter Deformation in einem Gebiet beitragen können, das nach den vorliegenden Modellen einer tektonischen Verformung widerstehen sollte, sind: (i) erhöhte tektonische Kräfte durch ein steileres Abtauchen der Platte und (ii) Schwächezonen in der Kruste, die auf prä-andine Deformationsereignisse zurückgehen. Schließlich könnten die thermomechanischen Bedingungen in diesem Teil des Vorlands einchlüsselfaktor für die Erhaltung des flachen Subduktionswinkels in diesen Breiten der SZA sein. KW - Andes KW - Anden KW - subduction KW - Subduktion KW - lithosphere KW - Lithosphäre KW - earthquakes KW - Erdbeben KW - modelling KW - Modellierung Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-549275 ER - TY - THES A1 - Schütte, Moritz T1 - Evolutionary fingerprints in genome-scale networks T1 - Evolutionäre Spuren in genomskaligen Netzwerken N2 - Mathematical modeling of biological phenomena has experienced increasing interest since new high-throughput technologies give access to growing amounts of molecular data. These modeling approaches are especially able to test hypotheses which are not yet experimentally accessible or guide an experimental setup. One particular attempt investigates the evolutionary dynamics responsible for today's composition of organisms. Computer simulations either propose an evolutionary mechanism and thus reproduce a recent finding or rebuild an evolutionary process in order to learn about its mechanism. The quest for evolutionary fingerprints in metabolic and gene-coexpression networks is the central topic of this cumulative thesis based on four published articles. An understanding of the actual origin of life will probably remain an insoluble problem. However, one can argue that after a first simple metabolism has evolved, the further evolution of metabolism occurred in parallel with the evolution of the sequences of the catalyzing enzymes. Indications of such a coevolution can be found when correlating the change in sequence between two enzymes with their distance on the metabolic network which is obtained from the KEGG database. We observe that there exists a small but significant correlation primarily on nearest neighbors. This indicates that enzymes catalyzing subsequent reactions tend to be descended from the same precursor. Since this correlation is relatively small one can at least assume that, if new enzymes are no "genetic children" of the previous enzymes, they certainly be descended from any of the already existing ones. Following this hypothesis, we introduce a model of enzyme-pathway coevolution. By iteratively adding enzymes, this model explores the metabolic network in a manner similar to diffusion. With implementation of an Gillespie-like algorithm we are able to introduce a tunable parameter that controls the weight of sequence similarity when choosing a new enzyme. Furthermore, this method also defines a time difference between successive evolutionary innovations in terms of a new enzyme. Overall, these simulations generate putative time-courses of the evolutionary walk on the metabolic network. By a time-series analysis, we find that the acquisition of new enzymes appears in bursts which are pronounced when the influence of the sequence similarity is higher. This behavior strongly resembles punctuated equilibrium which denotes the observation that new species tend to appear in bursts as well rather than in a gradual manner. Thus, our model helps to establish a better understanding of punctuated equilibrium giving a potential description at molecular level. From the time-courses we also extract a tentative order of new enzymes, metabolites, and even organisms. The consistence of this order with previous findings provides evidence for the validity of our approach. While the sequence of a gene is actually subject to mutations, its expression profile might also indirectly change through the evolutionary events in the cellular interplay. Gene coexpression data is simply accessible by microarray experiments and commonly illustrated using coexpression networks where genes are nodes and get linked once they show a significant coexpression. Since the large number of genes makes an illustration of the entire coexpression network difficult, clustering helps to show the network on a metalevel. Various clustering techniques already exist. However, we introduce a novel one which maintains control of the cluster sizes and thus assures proper visual inspection. An application of the method on Arabidopsis thaliana reveals that genes causing a severe phenotype often show a functional uniqueness in their network vicinity. This leads to 20 genes of so far unknown phenotype which are however suggested to be essential for plant growth. Of these, six indeed provoke such a severe phenotype, shown by mutant analysis. By an inspection of the degree distribution of the A.thaliana coexpression network, we identified two characteristics. The distribution deviates from the frequently observed power-law by a sharp truncation which follows after an over-representation of highly connected nodes. For a better understanding, we developed an evolutionary model which mimics the growth of a coexpression network by gene duplication which underlies a strong selection criterion, and slight mutational changes in the expression profile. Despite the simplicity of our assumption, we can reproduce the observed properties in A.thaliana as well as in E.coli and S.cerevisiae. The over-representation of high-degree nodes could be identified with mutually well connected genes of similar functional families: zinc fingers (PF00096), flagella, and ribosomes respectively. In conclusion, these four manuscripts demonstrate the usefulness of mathematical models and statistical tools as a source of new biological insight. While the clustering approach of gene coexpression data leads to the phenotypic characterization of so far unknown genes and thus supports genome annotation, our model approaches offer explanations for observed properties of the coexpression network and furthermore substantiate punctuated equilibrium as an evolutionary process by a deeper understanding of an underlying molecular mechanism. N2 - Die biologische Zelle ist ein sehr kompliziertes Gebilde. Bei ihrer Betrachtung gilt es, das Zusammenspiel von Tausenden bis Millionen von Genen, Regulatoren, Proteinen oder Molekülen zu beschreiben und zu verstehen. Durch enorme Verbesserungen experimenteller Messgeräte gelingt es mittlerweile allerdings in geringer Zeit enorme Datenmengen zu messen, seien dies z.B. die Entschlüsselung eines Genoms oder die Konzentrationen der Moleküle in einer Zelle. Die Systembiologie nimmt sich dem Problem an, aus diesem Datenmeer ein quantitatives Verständnis für die Gesamtheit der Wechselwirkungen in der Zelle zu entwickeln. Dabei stellt die mathematische Modellierung und computergestützte Analyse ein eminent wichtiges Werkzeug dar, lassen sich doch am Computer in kurzer Zeit eine Vielzahl von Fällen testen und daraus Hypothesen generieren, die experimentell verifiziert werden können. Diese Doktorarbeit beschäftigt sich damit, wie durch mathematische Modellierung Rückschlüsse auf die Evolution und deren Mechanismen geschlossen werden können. Dabei besteht die Arbeit aus zwei Teilen. Zum Einen wurde ein Modell entwickelt, dass die Evolution des Stoffwechsels nachbaut. Der zweite Teil beschäftigt sich mit der Analyse von Genexpressionsdaten, d.h. der Stärke mit der ein bestimmtes Gen in ein Protein umgewandelt, "exprimiert", wird. Der Stoffwechsel bezeichnet die Gesamtheit der chemischen Vorgänge in einem Organismus; zum Einen werden Nahrungsstoffe für den Organismus verwertbar zerlegt, zum Anderen aber auch neue Stoffe aufgebaut. Da für nahezu jede chemische Reaktion ein katalysierendes Enzym benötigt wird, ist davon auszugehen, dass sich der Stoffwechsel parallel zu den Enzymen entwickelt hat. Auf dieser Annahme basiert das entwickelte Modell zur Enzyme-Stoffwechsel-Koevolution. Von einer Anfangsmenge von Enzymen und Molekülen ausgehend, die etwa in einer primitiven Atmosphäre vorgekommen sind, werden sukzessive Enzyme und die nun katalysierbaren Reaktionen hinzugefügt, wodurch die Stoffwechselkapazität anwächst. Die Auswahl eines neuen Enzyms geschieht dabei in Abhängigkeit von der Ähnlichkeit mit bereits vorhandenen und ist so an den evolutionären Vorgang der Mutation angelehnt: je ähnlicher ein neues Enzym zu den vorhandenen ist, desto schneller kann es hinzugefügt werden. Dieser Vorgang wird wiederholt, bis der Stoffwechsel die heutige Form angenommen hat. Interessant ist vor allem der zeitliche Verlauf dieser Evolution, der mittels einer Zeitreihenanalyse untersucht wird. Dabei zeigt sich, dass neue Enzyme gebündelt in Gruppen kurzer Zeitfolge auftreten, gefolgt von Intervallen relativer Stille. Dasselbe Phänomen kennt man von der Evolution neuer Arten, die ebenfalls gebündelt auftreten, und wird Punktualismus genannt. Diese Arbeit liefert somit ein besseres Verständnis dieses Phänomens durch eine Beschreibung auf molekularer Ebene. Im zweiten Projekt werden Genexpressionsdaten von Pflanzen analysiert. Einerseits geschieht dies mit einem eigens entwickelten Cluster-Algorithmus. Hier läßt sich beobachten, dass Gene mit einer ähnlichen Funktion oft auch ein ähnliches Expressionsmuster aufweisen. Das Clustering liefert einige Genkandidaten, deren Funktion bisher unbekannt war, von denen aber nun vermutet werden konnte, dass sie enorm wichtig für das Wachstum der Pflanze sind. Durch Experimente von Pflanzen mit und ohne diese Gene zeigte sich, dass sechs neuen Genen dieses essentielle Erscheinungsbild zugeordnet werden kann. Weiterhin wurden Netzwerke der Genexpressionsdaten einer Pflanze, eines Pilzes und eines Bakteriums untersucht. In diesen Netzwerken werden zwei Gene verbunden, falls sie ein sehr ähnliches Expressionsprofil aufweisen. Nun zeigten diese Netzwerke sehr ähnliche und charakteristische Eigenschaften auf. Im Rahmen dieser Arbeit wurde daher ein weiteres evolutionäres Modell entwickelt, das die Expressionsprofile anhand von Duplikation, Mutation und Selektion beschreibt. Obwohl das Modell auf sehr simplen Eigenschaften beruht, spiegelt es die beobachteten Eigenschaften sehr gut wider, und es läßt sich der Schluss ziehen, dass diese als Resultat der Evolution betrachtet werden können. Die Ergebnisse dieser Arbeiten sind als Doktorarbeit in kumulativer Form bestehend aus vier veröffentlichten Artikeln vereinigt. KW - Systembiologie KW - Modellierung KW - Evolution KW - Stoffwechsel KW - Gen-Koexpression KW - Systems Biology KW - Modeling KW - Evolution KW - Metabolism KW - Gene co-expression Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-57483 ER -