On the inheritance and mechanism of baculovirus resistance of the codling moth, Cydia pomonella (L.)

Das Cydia pomonella Granulovirus (CpGV, Baculoviridae) wird seit Ende der 1980er Jahre als hoch-selektives und effizientes biologisches Bekämpfungsmittel zur Kontrolle des Apfelwicklers im Obstanbau eingesetzt. Seit 2004 wurden in Europa verschiedene Apfelwicklerpopulationen beobachtet die resistent gegenüber dem hauptsächlich angewendeten Isolat CpGV-M aufweisen. Die vorliegende Arbeit befasst sich mit der Untersuchung der Vererbung und des Mechanismus der CpGV Resistenz. Einzelpaarkreuzungen zwischen einem empfindlichen Laborstamm (CpS) und einem homogen resistenten Stamm (CpRR1) zeigten, dass die Resistenz durch ein einziges dominantes Gen, das auf dem Z-Chromosom lokalisiert ist, vererbt wird. Massernkreuzungen zwischen CpS und einer heterogen resistenten Feldpopulation (CpR) deuteten zunächst auf einen unvollständig dominanten autosomalen Erbgang hin. Einzelpaarkreuzungen zwischen CpS und CpR bewiesen jedoch, dass die Resistenz in CpR ebenfalls monogen dominant und geschlechtsgebunden auf dem Z-Chromosom vererbt wird. Diese Arbeit diskutiert zudem die Vor- und Nachteile von Einzelpaarkreuzungen gegenüber Massernkreuzungen bei der Untersuchung von Vererbungsmechanismen. Die Wirksamkeit eines neuen CpGV Isolates aus dem Iran (CpGV-I12) gegenüber CpRR1 Larven, wurde in Bioassays getestet. Die Ergebnisse zeigen, dass CpGV-I12 die Resistenz in allen Larvenstadien von CpRR1 brechen kann und fast so gut wirkt wie CpGV-M gegenüber CpS Larven. Daher ist CpGV-I12 für die Kontrolle des Apfelwicklers in Anlagen wo die Resistenz aufgetreten ist geeignet. Um den der CpGV Resistenz zugrunde liegenden Mechanismus zu untersuchen, wurden vier verschiedene Experimente durchgeführt: 1) die peritrophische Membran degradiert indem ein optischer Aufheller dem virus-enthaltenden Futtermedium beigefügt wurde. Das Entfernen dieser mechanischen Schutzbarriere, die den Mitteldarm auskleidet, führte allerdings nicht zu einer Reduzierung der Resistenz in CpR Larven. Demnach ist die peritrophische Membran nicht am Resistenzmechanismus beteiligt. 2) Die Injektion von Budded Virus in das Hämocoel führte nicht zur Brechung der Resistenz. Folglich die die Resistenz nicht auf den Mitteldarm beschränkt, sondern auch in der Sekundärinfektion wirksam. 3) Die Replikation von CpGV in verschiedenen Geweben (Mitteldarm, Hämolymphe und Fettkörper) von CpS und CpRR1 wurde mittels quantitativer PCR verfolgt. In CpS Larven konnte in allen drei Gewebetypen sowohl nach oraler als auch nach intra-hämocoelarer Infektion eine Zunahme der CpGV Genome in Abhängigkeit der Zeit festgestellt werden. Dagegen konnte in den Geweben aus CpRR1 nach oraler sowie intra-hämocoelarer Infektion keine Virusreplikation detektiert werden. Dies deutet darauf hin, dass die CpGV Resistenz in allen Zelltypen präsent ist. 4) Um zu untersuchen ob ein humoraler Faktor in der Hämolymphe ursächlich an der Resistenz beteiligt ist, wurde Hämolymphe aus CpRR1 Larven in CpS Larven injiziert und diese anschließend oral mit CpGV infiziert. Es konnte jedoch keine Immunreaktion beobachtet und kein Faktor in der Hämolymphe identifiziert werden, der Resistenz induzieren könnte. Auf Grundlage dieser Ergebnisse kann festgestellt werden, dass in resistenten Apfelwicklerlarven die virale Replikation in allen Zelltypen verhindert wird, was auf eine Virus-Zell Inkompatibilität hinweist. Da in CpRR1 keine DNA Replikation beobachtet wurde, wird die CpGV Resistenz wahrscheinlich durch eine frühe Unterbindung der Virusreplikation verursacht.Das früh exprimierte Gen pe38 codiert für ein Protein, das wahrscheinlich für die Resistenzbrechung durch CpGV-I12 verantwortlich ist. Interaktionen zwischen dem Protein PE38 und Proteinen in CpRR1 wurden mit Hilfe des Yeast Two-Hybrid (Y2H) Systems untersucht. Die detektierten Interaktionen sind noch nicht durch andere Methoden bestätigt, jedoch wurden zwei mögliche Gene auf dem Z-Chromosom und eines auf Chromosom 15 gefunden, wie möglicherweise an der CpGV Resistenz beteiligt sind.

Parameterschätzung in zeitdiskreten ergodischen Markov-Prozessen am Beispiel des Cox-Ingersoll-Ross-Modells

In dieser Arbeit geht es um die Schätzung von Parametern in zeitdiskreten ergodischen Markov-Prozessen im allgemeinen und im CIR-Modell im besonderen. Beim CIR-Modell handelt es sich um eine stochastische Differentialgleichung, die von Cox, Ingersoll und Ross (1985) zur Beschreibung der Dynamik von Zinsraten vorgeschlagen wurde. Problemstellung ist die Schätzung der Parameter des Drift- und des Diffusionskoeffizienten aufgrund von äquidistanten diskreten Beobachtungen des CIR-Prozesses. Nach einer kurzen Einführung in das CIR-Modell verwenden wir die insbesondere von Bibby und Sørensen untersuchte Methode der Martingal-Schätzfunktionen und -Schätzgleichungen, um das Problem der Parameterschätzung in ergodischen Markov-Prozessen zunächst ganz allgemein zu untersuchen. Im Anschluss an Untersuchungen von Sørensen (1999) werden hinreichende Bedingungen (im Sinne von Regularitätsvoraussetzungen an die Schätzfunktion) für die Existenz, starke Konsistenz und asymptotische Normalität von Lösungen einer Martingal-Schätzgleichung angegeben. Angewandt auf den Spezialfall der Likelihood-Schätzung stellen diese Bedingungen zugleich lokal-asymptotische Normalität des Modells sicher. Ferner wird ein einfaches Kriterium für Godambe-Heyde-Optimalität von Schätzfunktionen angegeben und skizziert, wie dies in wichtigen Spezialfällen zur expliziten Konstruktion optimaler Schätzfunktionen verwendet werden kann. Die allgemeinen Resultate werden anschließend auf das diskretisierte CIR-Modell angewendet. Wir analysieren einige von Overbeck und Rydén (1997) vorgeschlagene Schätzer für den Drift- und den Diffusionskoeffizienten, welche als Lösungen quadratischer Martingal-Schätzfunktionen definiert sind, und berechnen das optimale Element in dieser Klasse. Abschließend verallgemeinern wir Ergebnisse von Overbeck und Rydén (1997), indem wir die Existenz einer stark konsistenten und asymptotisch normalen Lösung der Likelihood-Gleichung zeigen und lokal-asymptotische Normalität für das CIR-Modell ohne Einschränkungen an den Parameterraum beweisen.

Phylogenie und Biogeographie der Linoideae (Linaceae) - Evolution von Heterostylie, Homostylie in Linum

Die Linaceae-Linoideae, vor allem die Gattung Linum, wurden unter Verwendung von zwei molekularen Markern (rbcL und ITS) bzgl. ihrer Phylogenie und Biogeographie untersucht. Die Linaceae entstanden während der mittleren Kreide in den frühen tropischen Regenwäldern, von wo aus sich die monophyletischen Linoideae vor etwa 51-46 Mill. Jahren über die temperaten Gebiete der Nordhemisphäre ausbreiteten. Während die drei basal abspaltenden Gattungen Anisadenia, Reinwardtia und Tirpitzia bzgl. ihrer Verbreitung auf Südostasien beschränkt sind, ist die Gattung Linum heute auf allen Kontinenten vertreten. Der Ursprung von Linum liegt wahrscheinlich in Südwestasien bzw. dem östlichen Mediterraneum, wo es im Oligozän zur Aufspaltung in zwei Entwicklungslinien kam ('Blaue Gruppe' und 'Gelbe Gruppe'). Während die überwiegend blaublühenden Linum-Arten ('Blaue Gruppe') vor allem in Europa und Südwestasien vorkommen, weisen die Vertreter der 'Gelben Gruppe' ein wesentlich größeres Verbreitungsgebiet auf. Gelbblühende Linum Arten findet man auf allen Kontinenten mit Diversitätszentren in Nordostamerika und Südwestasien. Interessanterweise wurde Amerika zweimal unabhängig voneinander besiedelt. Während die gelbblühenden Arten vor etwa 22-20 Mill. Jahren von Westeuropa über den Atlantik den amerikanischen Kontinent erreichten, wanderten Vertreter der 'Blauen Gruppe' im Pliozän (vor 3.78-3.33 Mill Jahren) über die Bering-Landbrücke in die Neue Welt ein. Auch in Südafrika sind einige gelbblühende Linum-Arten zu verzeichnen, die nicht über Nordafrika (wo einige Arten der 'Gelben Gruppe' beheimatet sind) die südliche Spitze des Kontinents erreichten, sondern von Amerika aus. Die molekularphylogenetischen Ergebnisse legen eine Eingliederung der Gattungen Cliococca, Hesperolinon, Radiola und Sclerolinon in Linum nahe, die durch morphologische Merkmale gestützt wird. Linopsis, die artenreichste Sektion der Gattung Linum, bedarf einiger Umstrukturierungen auf der Basis der molekularen und morphologischen Daten. Ein interessantes Phänomen innerhalb der Linaceae ist das Vorkommen von heterostylen und homostylen Arten innerhalb der Familie. Die Kombination der molekular-phylogenetischen Ergebnisse mit morphologischen Beobachtungen des Reproduktionssystems lassen darauf schließen, dass sich Homostylie innerhalb von Linum mehrfach unabhängig voneinander entwickelt hat. Das Modell von Primula wurde als Grundlage verwendet, um Aufschluss über die Entstehung der Homostylie innerhalb von Linum zu erlangen. Aus Primula ist bekannt, dass eine Kopplungsgruppe aus mindestens drei Genen an der Vererbung von Heterostylie beteiligt ist: G/g kodiert hierbei die Griffellänge und die Selbstinkompatibilitäts-reaktion der Narbe, A/a die Länge der Filamente und P/p die Selbst-inkompatibilitätsreaktion des Pollens. Umfangreiche Kreuzungs-experimente einer homostylen und einer heterostylen Linum-Art deuten darauf hin, dass die Genotypen der beiden Blütenformen in heterostylen Linum-Arten denen in Primula entsprechen. Langgriffel sind hiernach homozygot rezessiv (gpa/gpa), während die Kurzgriffel heterozygot sind (GPA/gpa). Selbstkompatible, homostyle Arten können theoretisch durch verschiedene Rekombinations-ereignisse entstehen. Erste Ergebnisse der rasterelektronen-mikroskopischen Betrachtung der Pollenkornoberflächen und Narbenpapillen deuten darauf hin, dass innerhalb von Linum Homostylie durch unterschiedliche Rekombinations-ereignisse mehrfach aus heterostylen Arten entstanden ist. So besitzt die homostyle Linum leonii den Genotyp gPA/gPA, während für die homostylen L. tenuifolium und L. nodiflorum der Genotyp Gpa/Gpa wahrscheinlich ist.

Compliance lebertransplantierter Patienten mit dem Immunsuppressivum Tacrolimus bei 2x täglicher im Vergleich zu 1x täglicher Einnahme

Compliance lebertransplantierter Patienten mit der immunsuppressiven Therapie ist unerlässlich für den lang-fristigen Erfolg der Lebertransplantation. Aus Non-Compliance mit der immunsuppressiven Therapie können Abstoßungsreaktionen, Organverlust oder sogar Tod resultieren. Hauptziel der vorliegenden Studie war die erstmalige Evaluation der Compliance bei Einnahme von Prograf® (zweimal tägliche Einnahme von Tacrolimus) im Vergleich zur Einnahme von Advagraf® (einmal tägliche Einnahme von Tacrolimus). Von Interesse war außerdem die Fragestellung, ob sich die Compliance bezüglich der immunsuppressiven Therapie mit dem Zeitabstand zur Transplantation verändert. rnDie Compliancemessung wurde offen mittels MEMS® (Aardex Ltd., Schweiz) durchgeführt, der Patient war also über die Compliancekontrolle informiert. Mittels MEMS® konnten Datum und Uhrzeit der Dosisentnahme dokumentiert und damit zuverlässig das gesamte Compliancemuster über im Durchschnitt 176 Tage mit der zweimal täglichen Einnahme und 188 Tage mit der einmal täglichen Einnahme pro Patient erfasst werden. 65 Patienten mit dem Basisimmunsuppressivum Prograf® wurden in die prospektive, nicht-interventionelle Studie eingeschlossen und nach Per Protokoll-Analyse konnten die Daten von 63 in Mainz lebertransplantierten Patienten ausgewertet werden (Prograf®: Gruppe 1: 15 Patienten (Pat.), Gruppe 2: 23 Pat., Gruppe 3: 22 Pat., Drop-outs: 3 Pat.; Advagraf®: Gruppe 1: 16 Pat., Gruppe 2: 23 Pat., Gruppe 3: 23 Pat., Drop-outs: 1 Pat.). Die Dosing Compliance (DC), definiert als Prozent der Tage, an denen der MEMS®-Behälter korrekt geöffnet und die Dosis höchstwahrscheinlich korrekt eingenommen wurde, war der primäre Zielparameter. Weitere Methoden der Compliancemessung, wie der Pill Count, mehrere Fragebögen (Selbsteinschätzung, Patientenwissen-, Morisky-, MESI-, HADS-, SF-36- und Patientenzufriedenheit-Fragebogen) sowie die Blutspiegelmessung wurden eingesetzt, um die Compliance der Patienten umfassend charakterisieren zu können. rnDer Median der DC mit der zweimal täglichen Einnahme betrug 97% bei Pat. > 6 m.p.t. < 2 y.p.t., 97% bei Pat. > 2 y.p.t. < 5 y.p.t. und 98% bei Pat. > 5 y.p.t. (p=0,931; Kruskal-Wallis-Test). Der Median der DC von Tacroli-mus bei einmal täglicher Einnahme (Advagraf®) betrug 99% bei Pat. > 6 m.p.t. < 2 y.p.t., 98% bei Pat. > 2 y.p.t. < 5 y.p.t. und 97% bei Pat. > 5 y.p.t. (p=0,158; Kruskal-Wallis-Test). Insgesamt zeigten die Patienten während des gesamten Beobachtungszeitraums von 12 Monaten eine gute Compliance für die Einnahme ihres Immun-suppressivums. Die Timing Compliance (TiC)-raten lagen auf einem niedrigeren Niveau als die Dosing- und Taking Compliance (TC)-raten. Die Complianceraten der drei Subgruppen unterschieden sich nicht signifikant. Die Patienten mit dem geringsten Abstand zur Transplantation zeigten bei beinahe allen Messmethoden die höchste Compliance im Gegensatz zur etwas geringeren Compliance der Patienten mit größerem Abstand zur Transplantation. Die während der Advagraf®-Phase mittels MEMS® gemessenen DC-, TC- und TiC-raten fielen höher aus als bei Einnahme von Prograf® (p(DC)=0,003; p(TC)=0,077; p(TiC)=0,003; Wilcoxon Vorzeichen-Rang-Test). Dieses Ergebnis untermauert die in anderen Indikationen gefundene Complianceverbesserung durch die einmal tägliche Arzneimittelgabe im Vergleich zur zweimal täglichen Gabe. Die Auswertung der Drug Holidays ergab für die Advagraf®-Phase hingegen niedrigere Complianceraten als für die Prograf®-Phase. Dieses Ergebnis ist auf die Definition des Drug Holidays (keine Arzneimitteleinnahme über 48 h) zurück zu führen. Die Chance Advagraf® einmal pro Tag zu vergessen ist doppelt so hoch, als Prograf® dreimal aufeinander fol-gend zu vergessen. Mit einer verhältnismäßigeren Definition von Drug Holidays (Einnahmepause von 72 Stun-den bei einmal täglicher Einnahme von Advagraf® entsprechend drei ausgelassenen Dosen von Prograf®) ist die Compliancerate 81%. Die Ergebnisse des Pill Counts waren sowohl bei Einnahme von Prograf® als auch von Advagraf® mit der jeweils gemessenen TC vergleichbar, was die Zuverlässigkeit der Messergebnisse bes-tätigt. rnDie zusätzlich eingesetzten Methoden verifizierten das Ergebnis der höheren Compliance mit der einmal tägli-chen Einnahme. Die während der Advagraf®-Phase beantworteten Fragebögen zeigten einen Trend zu besserer Compliance und Lebensqualität. Lediglich die Ergebnisse des MESI-Fragebogens und der Blutspiegelmessungen wichen sowohl während der Prograf®- als auch während der Advagraf®-Phase stark von den Ergebnis-sen der anderen Methoden ab. rnUnter Einbeziehung aller mittels MEMS® und Pill Count objektiv gemessenen Complianceparameter konnten während der Prograf®-Einnahme 54 von 60 Pat. (90%) und während der Advagraf®-Phase 59 von 62 Pat. (95%) als compliant eingestuft werden. Aufgrund subjektiver Compliancemessungen waren 49 von 58 Pat. (84%) während der Prograf®- und 54 von 59 Pat. (92%) während der Advagraf®-Phase als compliant einzustufen. Es wurde beobachtet, dass die zeitlich korrekte Einnahme der Morgendosis einfacher und bei Einmalgabe zu bevorzugen ist. Die wochentagsbezogene Auswertung ergab erwartungsgemäß, dass am Wochenende (Samstag und Sonntag) am häufigsten Dosen ausgelassen wurden. rnDie Umstellung von Prograf® auf Advagraf® stellte kein Problem dar. Beinahe alle Patienten waren dankbar und zufrieden mit der Reduzierung der Dosierungsfrequenz und der größeren Unabhängigkeit durch die entfallene abendliche Einnahme. Der positive Einfluss der geringeren Dosierungshäufigkeit auf die Langzeitcompliance der Patienten, ist ein hinreichender Grund die Entwicklung von Formulierungen zur einmal täglichen Ein-nahme für weitere Immunsuppressiva zu fordern. Insbesondere bei den häufig eingesetzten Kombinationstherapien von Immunsuppressiva würde der Effekt der Complianceverbesserung noch verstärkt werden, wenn alle eingesetzten Immunsuppressiva zur einmal täglichen Gabe geeignet wären.

Ergodicity and regularity of invariant measure for branching Markov processes with immigration

In this thesis we consider systems of finitely many particles moving on paths given by a strong Markov process and undergoing branching and reproduction at random times. The branching rate of a particle, its number of offspring and their spatial distribution are allowed to depend on the particle's position and possibly on the configuration of coexisting particles. In addition there is immigration of new particles, with the rate of immigration and the distribution of immigrants possibly depending on the configuration of pre-existing particles as well. In the first two chapters of this work, we concentrate on the case that the joint motion of particles is governed by a diffusion with interacting components. The resulting process of particle configurations was studied by E. Löcherbach (2002, 2004) and is known as a branching diffusion with immigration (BDI). Chapter 1 contains a detailed introduction of the basic model assumptions, in particular an assumption of ergodicity which guarantees that the BDI process is positive Harris recurrent with finite invariant measure on the configuration space. This object and a closely related quantity, namely the invariant occupation measure on the single-particle space, are investigated in Chapter 2 where we study the problem of the existence of Lebesgue-densities with nice regularity properties. For example, it turns out that the existence of a continuous density for the invariant measure depends on the mechanism by which newborn particles are distributed in space, namely whether branching particles reproduce at their death position or their offspring are distributed according to an absolutely continuous transition kernel. In Chapter 3, we assume that the quantities defining the model depend only on the spatial position but not on the configuration of coexisting particles. In this framework (which was considered by Höpfner and Löcherbach (2005) in the special case that branching particles reproduce at their death position), the particle motions are independent, and we can allow for more general Markov processes instead of diffusions. The resulting configuration process is a branching Markov process in the sense introduced by Ikeda, Nagasawa and Watanabe (1968), complemented by an immigration mechanism. Generalizing results obtained by Höpfner and Löcherbach (2005), we give sufficient conditions for ergodicity in the sense of positive recurrence of the configuration process and finiteness of the invariant occupation measure in the case of general particle motions and offspring distributions.

Existence of solutions and asymtptotic limits of the Euler-Poisson and the quantum drift-diffusion systems

My work concerns two different systems of equations used in the mathematical modeling of semiconductors and plasmas: the Euler-Poisson system and the quantum drift-diffusion system. The first is given by the Euler equations for the conservation of mass and momentum, with a Poisson equation for the electrostatic potential. The second one takes into account the physical effects due to the smallness of the devices (quantum effects). It is a simple extension of the classical drift-diffusion model which consists of two continuity equations for the charge densities, with a Poisson equation for the electrostatic potential. Using an asymptotic expansion method, we study (in the steady-state case for a potential flow) the limit to zero of the three physical parameters which arise in the Euler-Poisson system: the electron mass, the relaxation time and the Debye length. For each limit, we prove the existence and uniqueness of profiles to the asymptotic expansion and some error estimates. For a vanishing electron mass or a vanishing relaxation time, this method gives us a new approach in the convergence of the Euler-Poisson system to the incompressible Euler equations. For a vanishing Debye length (also called quasineutral limit), we obtain a new approach in the existence of solutions when boundary layers can appear (i.e. when no compatibility condition is assumed). Moreover, using an iterative method, and a finite volume scheme or a penalized mixed finite volume scheme, we numerically show the smallness condition on the electron mass needed in the existence of solutions to the system, condition which has already been shown in the literature. In the quantum drift-diffusion model for the transient bipolar case in one-space dimension, we show, by using a time discretization and energy estimates, the existence of solutions (for a general doping profile). We also prove rigorously the quasineutral limit (for a vanishing doping profile). Finally, using a new time discretization and an algorithmic construction of entropies, we prove some regularity properties for the solutions of the equation obtained in the quasineutral limit (for a vanishing pressure). This new regularity permits us to prove the positivity of solutions to this equation for at least times large enough.

Geochemistry, petrogenesis and tectonic setting of ophiolites and mafic-ultramafic complexes in the Northeastern Aegean Region

In this study two ophiolites and a mafic-ultramafic complexes of the northeastern Aegean Sea, Greece, have been investigated to re-evaluate their petrogenetic evolution and tectonic setting. These complexes are: the mafic-ultramafic complex of Lesvos Island and the ophiolites of Samothraki Island and the Evros area. In order to examine these complexes in detail whole-rock major- and trace-elements as well as Sr and Nd isotopes, and minerals were analysed and U-Pb SHRIMP ages on zircons were determined. The mafic-ultramafic complex of Lesvos Island consists of mantle peridotite thrusted over a tectonic mélange containing metasediments, metabasalts and a few metagabbros. This succession had previously been interpreted as an ophiolite of Late Jurassic age. The new field and geochemical data allow a reinterpretation of this complex as representing an incipient continental rift setting that led to the subsequent formation of the Meliata-Maliac-Vardar branches of Neotethys in Upper Permian times (253 ± 6 Ma) and the term “Lesvos ophiolite” should be abandoned. With proceeding subduction and closure of the Maliac Ocean in Late Jurassic times (155 Ma) the Lesvos mafic-ultramafic complex was obducted. Zircon ages of 777, 539 and 338 Ma from a gabbro strongly suggest inheritance from the intruded basement and correspond to ages of distinct terranes recently recognized in the Hellenides (e.g. Florina terrane). Geochemical similar complexes which contain rift associations with Permo-Triassic ages can be found elsewhere in Greece and Turkey, namely the Teke Dere Thrust Sheet below the Lycian Nappes (SW Turkey), the Pindos subophiolitic mélange (W Greece), the Volcanosedimentary Complex on Central Evia Island (Greece) and the Karakaya Complex (NW Turkey). This infers that the rift-related rocks from Lesvos belong to an important Permo-Triassic rifting episode in the eastern Mediterranean. The ‘in-situ’ ophiolite of Samothraki Island comprises gabbros, sparse dykes and basalt flows as well as pillows cut by late dolerite dykes and had conventionally been interpreted as having formed in an ensialic back-arc basin. The results of this study revealed that none of the basalts and dolerites resemble mid-ocean ridge or back-arc basin basalts thus suggesting that the Samothraki ophiolite cannot represent mature back-arc basin crust. The age of the complex is regarded to be 160 ± 5 Ma (i.e. Oxfordian; early Upper Jurassic), which precludes any correlation with the Lesvos mafic-ultramafic complex further south (253 ± 6 Ma; Upper Permian). Restoration of the block configuration in NE Greece, before extensional collapse of the Hellenic hinterland and exhumation of the Rhodope Metamorphic Core Complex (mid-Eocene to mid-Miocene), results in a continuous ophiolite belt from Guevgueli in the NW to Samothraki in the SE, thus assigning the latter to the Innermost Hellenic Ophiolite Belt. In view of the data of this study, the Samothraki ophiolite represents a rift propagation of the Sithonia ophiolite spreading ridge into the Chortiatis calc-alkaline arc. The ophiolite of the Evros area consists of a plutonic sequence comprising cumulate and non-cumulate gabbros with plagiogranite veins, and an extrusive sequence of basalt dykes, massive and pillow lavas as well as pyroclastic rocks. Furthermore, in the Rhodope Massif tectonic lenses of harzburgites and dunites can be found. All rocks are spatially separated. The analytical results of this study revealed an intra-oceanic island arc setting for the Evros ophiolitic rocks. During late Middle Jurassic times (169 ± 2 Ma) an intra-oceanic arc has developed above a northwards directed intra-oceanic subduction zone of the Vardar Ocean in front of the Rhodope Massif. The boninitic, island arc tholeiitic and calc-alkaline rocks reflect the evolution of the Evros island arc. The obduction of the ophiolitic rocks onto the Rhodope basement margin took place during closure of the Vardar ocean basins. The harzburgites and dunites of the Rhodope Massif are strongly depleted and resemble harzburgites from recent oceanic island arcs. After melt extraction they underwent enrichment processes by percolating melts and fluids from the subducted slab. The relationship of the peridotites and the Evros ophiolite is still ambiguous, but the stratigraphic positions of the peridotites and the ophiolitic rocks indicate separated origin. The harzburgites and dunites most probably represent remnants of the mantle wedge of the island arc of the Rhodope terrane formed above subducted slab of the Nestos Ocean in late Middle Jurassic times. During collision of the Thracia terrane with the Rhodope terrane thrusting of the Rhodope terrane onto the Thracia terrane took place, whereas the harzburgites and dunites were pushed between the two terranes now cropping out on top of the Thracia terrane of the Rhodope Massif.

Analysis and numerical solution of the Peterlin viscoelastic model

Liquids and gasses form a vital part of nature. Many of these are complex fluids with non-Newtonian behaviour. We introduce a mathematical model describing the unsteady motion of an incompressible polymeric fluid. Each polymer molecule is treated as two beads connected by a spring. For the nonlinear spring force it is not possible to obtain a closed system of equations, unless we approximate the force law. The Peterlin approximation replaces the length of the spring by the length of the average spring. Consequently, the macroscopic dumbbell-based model for dilute polymer solutions is obtained. The model consists of the conservation of mass and momentum and time evolution of the symmetric positive definite conformation tensor, where the diffusive effects are taken into account. In two space dimensions we prove global in time existence of weak solutions. Assuming more regular data we show higher regularity and consequently uniqueness of the weak solution. For the Oseen-type Peterlin model we propose a linear pressure-stabilized characteristics finite element scheme. We derive the corresponding error estimates and we prove, for linear finite elements, the optimal first order accuracy. Theoretical error of the pressure-stabilized characteristic finite element scheme is confirmed by a series of numerical experiments.