No description available.

Ce travail avait pour objectif d’étudier les effets de l’uranium sur le budget énergétique et la dynamique de population d’un micro-crustacé représentatif des écosystèmes aquatiques d’eau douce, Daphnia magna. L’étude expérimentale de la toxicité de l’uranium sur la physiologie (nutrition, respiration) et l’histoire de vie (survie, croissance, reproduction) deD. magna a été réalisée au travers d’expositions sur une, deux ou trois générations successives(F0, F1, F2), commencées avec des néonates issus de 1ère ou 5ème ponte, à des concentrations de0, 10, 25 et 75 μgU.L-1. Les résultats ont mis en évidence l’aggravation des effets au travers des générations (en partie liée à l’exposition des daphnies pendant l’embryogénèse) et la sensibilité plus forte des individus issus de 1ère ponte. Des réductions significatives des taux d’assimilation mesurés par méthode de marquage-traçage au 14C de la nourriture nous ont permis d’identifier un effet sur l’assimilation comme mode d’action de l’uranium, en accord avec les altérations importantes de la structure de la paroi du tube digestif révélées par microscopie optique.L’intégration des résultats dans un modèle de budget énergétique (DEBtox) a abouti à l’estimation de concentrations seuil sans effet (NEC) de 9.37, 8.21 et 2.31 μgU.L-1 au-delà desquelles le fonctionnement de l’organisme est altéré dans les générations F0, F1 et F2,respectivement. La combinaison du DEBtox avec des modèles matriciels a permis d’extrapoler les conséquences sur le taux de croissance asymptotique de la population (λ), critère plus pertinent dans un contexte écologique. Les simulations ont prédit une augmentation de l’impact de l’uranium au travers des générations avec une réduction de λ dans F0 et une extinction des populations pour des concentrations de 51-59 μgU.L-1 dans F1 et de 39-41 μgU.L-1 dans F2. Les simulations ont souligné l’importance de prendre en considération les individus les plus sensibles dans la détermination de la réponse des populations. / This work aimed to study effects of uranium on energy budget and population dynamicsin Daphnia magna a representative microcrustacean of freshwater ecosystems. An experimentalstudy of uranium toxicity on physiology (nutrition, respiration) and life history (survival, growthand reproduction) of D. magna was carried out, based on exposures over one, two or threesuccessive generations (F0, F1 and F2) started with neonates from 1st or 5th brood, at 0, 10, 25and 75 μgU.L-1. Results showed that toxic effects increased across generations (partially due todaphnid exposure during embryogenesis) and that individuals from 1st brood were moresensitive than individuals from 5th brood. Significant reductions in assimilation rates, measuredusing a radiotracing method with 14C-labelled food, allowed us to identify an effect onassimilation as the mode of action for uranium, in agreement with important damages in theintegrity of intestinal epithelium observed by optic microscopy. Integrating results in a dynamicenergy budget model (DEBtox) yielded estimated no effect concentrations (NEC) of 9.37, 8.21and 2.31 μgU.L-1 above which organism functions were altered in generations F0, F1 and F2,respectively. Combining DEBtox with matrix models allowed us to extrapolate consequences onasymptotic population growth rate (λ), a relevant endpoint in an ecological context. Simulationspredicted an increase in uranium impact across generations with reduction of λ in F0 andpopulation extinctions at 51-59 μgU.L-1 in F1 and 39-41 μgU.L-1 in F2. Simulations emphasizedthe importance of considering the most sensitive individuals while determining populationresponse

This study combines structural field data with microstructural observations in an analysis of a mid-crustal shear zone related to the emplacement of the Misty pluton during a high-flux magmatic event in Northern Fiordland, New Zealand. These high-flux magmatic events transport massive amounts of heat and material as they develop along accretionary continental margins, and represent a primary source of continental crust. Fiordland, New Zealand possesses, perhaps, the most extensive middle and lower crustal exposure of these systems on earth. Therefore, this study area provides a significant opportunity to understand processes of continental crust formation in the mid-crust and how these events relate to the broader construction of continents. Herein, I document the four-stage geologic history of the Cozette Burn field area. Pre-existing structures along the Gondwana accretionary margin hosted a regional flare-up magmatic event that produced the Misty pluton and several other large plutons of the West Fiordland Orthogneiss (WFO). This study primarily focuses on the mid-crustal emplacement of the Misty pluton during oblique convergence along the accretionary margin, forming the upper-amphibolite facies Misty Shear Zone (MSZ). The exposures of the MSZ within the Cozette Burn preserve rare structural relationships between host rock and the intrusive Misty pluton. Together, these structures developed during end-stage contractional tectonics that constructed a long-lived (~270+ Ma) composite batholith. Heterogeneous ductile shearing defines the MSZ, with microstructural evidence indicating an interplay of high-temperature crystal plastic deformation along with partial melting of host rock and melt channeling. This resulted in focused, melt-assisted shearing under regional transpressive deformation. These accommodative processes provided an efficient mechanism for moving heat, fluids and magma sourced from the lower crust/mantle boundary into the mid-crust during 15-25 km of crustal thickening related to arc flare-up magmatism. This flare up magmatism and MSZ formation occurred during the final stages of crustal thickening along Gondwana continental margin. High-strain, mylonitic- ultramylonitic shear zones developed in a later phase of deformation, cutting MSZ fabrics near contacts between the Misty pluton and host rock. These more localized shear zones can be attributed to either accommodation of localized melt-pressure buildup or the shift to extensional tectonics. Brittle faulting cut these structures with oblique-thrust in the Tertiary. These mid-crustal structures carry economic relevance: thickened-crust events along accretionary continental margins produce deep-crustal sourced, metal-bearing magmas that are transferred into mid-crust prior to their hydrothermal emplacement as ore deposits in the upper crust. The lasting influence of these processes warrants consideration when assessing continental crust architecture at all scales.

The processes associated with magma injection at shallow depths within the crust have been the topic of many geophysical studies, some investigating the seismicity that accompanies volcanic activity and others attempting to map the subsurface extent and geometry of the resulting magma bodies. The aim of this study is to obtain a better understanding of these processes by investigating the nature of seismic signals that accompany volcanic eruptions and by seismically imaging a magma body beneath a mid-ocean ridge, both located on, or adjacent to Iceland. The seismic phenomena associated with the 1996 Vatnajӧkull subglacial eruption in central Iceland, have been studied using data recorded by both temporary (HOTSPOT) and permanent (SIL) seismic networks. These networks comprise 60 broadband and short-period three-component seismographs and cover most parts of the country. Two very active volcanic systems, Bárdarbunga and Grimsvӧtn, are situated underneath the Vatnajokull ice cap. The volcanoseismic signals recorded there were categorised according to their waveform shape and frequency content, into three groups: (a) low-frequency events (1-2 Hz); (b) mixed-frequency events (1-4 Hz); and (c) volcanic tremor. The eruption was preceded by intense seismic activity which began with a = 5.6 earthquake located at the Bárdarbunga volcanic system. The epicentres of the earthquake swarm that followed the M(_w), = 5.6 event initially delineated the Bárdarbunga caldera rim and then migrated towards Grimsvӧtn, to a place where a fissure was later observed. Pre-eruptive tremor started at least two days before the eruption as a harmonic signal around five narrow frequency bands (0.5-0.7, 1.6, 2.2, 2.8 and 3.2 Hz). Co-eruptive tremor started as a broadband, continuous signal which evolved into low-amplitude background tremor interrupted by high-amplitude, cigar-shaped bursts. Further analysis revealed that continuous tremor and the cigar-shaped bursts had all the characteristics of low- dimensional chaotic signals. Geophysical and geochemical evidence suggest that a lateral migration of magma from Bárdarbunga facilitated the rupture of the roof of a magma chamber, situated at the fissure area, which subsequently erupted as tephra on the glacier. The second phase of the RAMESSES (Reykjanes Ridge Axial Melt Experiment: Structural Synthesis from Electromagnetic and Seismics) experiment involved the acquisition of multichannel seismic reflection data from 39 along- and across-axis lines shot over the magmatically active 57º 45'N axial volcanic ridge. The data from one along-axis line were processed using a variety of techniques that mainly aimed at reducing the large amount of coherent noise present, a result of scattered energy at the rough seabed. The final processed section revealed a number of reflection events that could be interpreted as intra-crustal reflections, originating from the interface between pillow lavas and sheeted dykes, and from the top part of a thin melt lens.

This work presents a new method to measure model independent viscosities of inhomogeneous materials at high temperatures. Many mechanisms driving volcanic eruptions are strongly influenced by the viscous properties of the participating materials. Since an eruption takes place at temperatures at which these materials (predominantly silicate melts) are not completely molten, typically inhomogeneities, like e.g. equilibrium and non-equilibrium crystals, are present in the system. In order to incorporate such inhomogeneities into objective material parameters the viscosity measurement is based on a rotational viscometer in a wide gap Couette setup. The gap size between the two concentric cylinders was designed as large as possible in order to account for the inhomogeneities. The emerging difficulties concerning the model independent data reduction from measured values to viscosities are solved using an appropriate interpolation scheme. The method was applied to a material representative for the majority of volcanic eruptions on earth: a typical continental basaltic rock (Billstein/Rhön/Germany). The measured viscosities show a strong shear rate dependency, which surprises, because basaltic melt has been, until now, assumed to behave as a Newtonian fluid. Since a non-Newtonian material shows a very different relaxation behavior in the Couette motion compared to a Newtonian one (which, ultimately, does not show any), and a strong relaxation signal was recorded during viscosity measurements, the equations of Couette motion were investigated. The time dependent stress distribution in a material due to a quasi step-like velocity change at the inner Couette radius (i.e. the spindle) was considered. The results show that a material combining a linear shear modulus and a Newtonian viscosity -- a Maxwell material -- cannot quantify the relaxation behavior. This could be considered as a hint, that the widely used Maxwell relaxation times cannot be applied as a 1:1 mapping from microscopic considerations to macroscopic situations. / Die vorliegende Arbeit beschreibt eine neue Methode zur modellunabhängigen Messung von Viskositäten bei hohen Temperaturen. Viele der Mechanismen, welche vulkanischer Aktivität zugrunde liegen, werden stark durch die viskosen Eigenschaften der beteiligten Materialien beeinflusst. Die eruptierten Materialien -- zum überwiegenden Teil Silikatschmelzen -- sind bei Eruptionstemperatur nicht komplett geschmolzen. Deshalb sind Gleichgewichts- und Nichtgleichgewichtskristalle in den betrachteten Systemen vorhanden. Um diese Inhomogenitäten in objektive Materialparameter einzubeziehen, basiert die vorgestellte Viskositätsmessung auf auf einem Rotationsviskosimeter in einer wide gap''-Anordnung. Die Spaltbreite zwischen den beiden konzentrischen Zylindern wurde so groß wie möglich gemacht um Inhomogenitäten zu berücksichtigen. Die aufkommenden Schwierigkeiten bezüglich der modellunabhängigen Bestimmung der Viscositäten aus den gemessenen Daten wurden mit einer geeigneten Interpolationsmethode gelöst. Mit dieser Methode wurden die Viskositäten eines, für die Mehrheit vulkanischer Eruptionen auf der Erde typisches Material gemessen: eines kontinentalen Basaltes aus Billstein (Rhön, Deutschland). Die gemessenen Viskositäten zeigen bei konstanter Temperatur eine starke Abhängigkeit von der Deformationsrate. Dies überrascht, da basaltische Schmelzen bis heute bei vergleichbaren Temperaturen als Newtonsche Flüssigkeiten betrachtet wurden. Da ein nicht-Newtonsches Material, im Vergleich mit einem Newtonschen, ein deutlich anderes Relaxationsverhalten aufweist (das Newtonsche zeigt ultimativ keine Relaxation), und da ein deutliches Relaxationssignal während der Viskositätsmessung gemessen wurde, wurden die Bewegungsgleichungen der Couette Bewegung untersucht. Die zeitabhängige Spannungeverteilung in einem Material, verursacht durch eine quasi-stufenartige Geschwindigkeitsänderung am inneren Couette-Radius (d. h. am Drehkörper des Viskosimeters) wurde betrachtet. Die Ergebnisse zeigen, dass ein Material, welches ein linear elastisches Schermodul und eine newtonsche Viskosität kombiniert -- ein Maxwell-Material -- das Relaxationverhalten quantitativ nicht beschreiben kann. Dies könnte als Hinweis betrachtet werden, dass die weitverbreiteten Maxwell-Relaxationszeiten nicht 1:1 von mikroskopischen Betrachtungen auf makroskopische Situationen angewendet werden können.

A plethora of magmatic processing occurs in magma reservoirs, where melts are stored prior to eruption. Magma reservoirs are complex, open systems, and often multiple reservoirs are partially inter-connected from source to surface, giving rise to the term 'volcanic plumbing system'. Parental melts feeding these reservoirs can have diverse and distinct geochemical and petrological characteristics, and be variably evolved or enriched. These melts can also bring with them a crystal cargo that may remain in equilibrium in the magma reservoir, but may also be modified by reaction, resorption, crystallisation and diffusion. Melts and crystals can be transported between reservoirs, from the upper mantle and through the crust, leading to melt mixing, reactions and volatile exsolution. Basaltic volcanic systems are fed by primitive melts, and due to the rapid ascent of melts and short magma storage times, these volcanoes provide the best means of unravelling the mantle and crustal contribution to geochemical heterogeneity observed in erupted samples. Despite the potential chemical complexity of a magma reservoir, evidence for magma processing and reaction can be preserved in melt inclusion suites and the compositional structure of their host crystals. Magmatic processes during storage and transport at two basaltic volcanoes are investigated using two carefully selected eruptions: the 1669 eruption at Mt. Etna, and the 2007 Father's Day eruption at Kīlauea. A suite of diverse geochemical, petrological and petrographical observations, made at a range of length-scales, are combined and interpreted in tandem with geophysical monitoring data. The conclusions of these studies shed light on the architecture of each volcano's plumbing systems and basaltic plumbing systems in general. This thesis is divided into two parts. The first study unravels the crustal and mantle processes controlling melt geochemical heterogeneity at Mt. Etna, Sicily, during the 1669 eruption, the largest eruption in historical times. The 1669 melt inclusion suite arises from the mixing of two basaltic melts with similar major element compositions but very different trace and volatile element compositions. The melt geochemistry suggests that at least one end-member melt has been heavily influenced by assimilation of carbonate in the crust. The elevation in alkalis, caused by carbonate assimilation, enhances carbon and sulfur solubility in one end member. The melt inclusion suite indicates that mixing of these melts occurred in the shallow crust shortly before eruption and this mixing may be the cause of the enhanced $CO_{2}$ fluxes prior to eruptions at Mt. Etna. The second study is split into two parts. Each uses the eruptive products of the Father's Day eruption at Kīlauea and aims to unravel the connectivity of the plumbing system between the summit and East Rift Zone, with a focus on timescales of storage and transport. The first part investigates the melt geochemistry in terms of heterogeneity and volatile composition, and the second investigates the crystal cargo in terms of features of the macro-scale crystal cargo distribution and the micro-scale geochemical zoning of individual crystals. The integration of observations and models from these two studies constrains the pressure, temperature and composition of source magma feeding the Father's Day eruption. The eruption is investigated in the context of the "magma surge'' event that preceded the intrusion, as well as within the context of the longer-term trends in Kīlauea geochemistry at the summit and East Rift Zone. Melt inclusion and matrix glass volatile systematics provide insights into the degassing path of the magma and the duration of magma transport to the surface is constrained by diffusion modelling. Estimated timescales for ascent by diffusion modelling of macrocryst major element composition, melt inclusion water content and the melt Fe$^{3+}$/Fe$_{tot}$ ratio are in agreement with timescales observed from the geophysical data of $< $8 hours from reservoir depth to eruption. Both studies emphasise how petrological observations can supplement geophysical monitoring datasets collected at the surface to aid our interpretation of volcanic behaviour and eruption forecasting.

Solidified igneous intrusions from originally liquid magma chambers display a large number of different sedimentary features. These features include the gravitational collapse of sidewalls producing slumps and the layering produced by gravitational settling of crystals. In the chamber fluid-dynamic processes such as convection are expected to occur due to cooling at the roof producing dense gravitationally unstable liquid, and the crystallisation of interstitial liquid changing the composition of the remaining liquid possibly reducing the density causing the liquid to rise up. The crystals which form in basaltic magma chambers have a high propensity to be mobilised due to convection and other fluid-dynamic processes including replenishment by a secondary intrusion. Convective mobilisation of plagioclase grains in vertical, tabular intrusions is seen from flat profiles of apparent aspect ratio as a function of dyke width. These flat profiles were formed due to scouring of gravitationally unstable sidewall mushes, and these crystals then become entrained in the convecting liquid. Convection only ceases once the volume of crystals in suspension reaches a critical volume fraction leading to an increase in viscosity, which dampens the vigour of convection. The majority of this study is performing and analysing a number of different experiments to look at the behaviour of different styles of analogue particle piles. Particle piles that are formed of inert, plastic particles are subjected to convection in the particle layer and in the bulk overlying fluid, and different styles of mobilisation depending on the heat flux driving convection and the density profile of the pile are observed. The mobilisation style goes from rolling of particles on the surface, to puffs of particles from the surface being lofted into the interior, followed by large particle fountains and then the entire particle pile being completely disaggregated and lofted into the interior of the chamber as the force driving convection is increased. The initiation of mobilisation can be explained by the fluidisation of a particle pile, whilst the high degrees of mobilisation seen in some high Rayleigh number regimes can be explaining by resuspending particles. In experiments where particle piles have a positive density profile (dense particles overlying low density particles) the underlying low density particles can break through the overlying layer in particle fountains and can be explained by a modified fluidisation parameter. These experiments lack the reactivity and cohesion that realistic crystal piles would have. To try and quantify this, I have also performed a series of experiments looking at the rheology of an ice-sucrose suspension, where ice crystals can sinter and aggregate together. Under sheared conditions the forces required to disaggregate ice aggregates can be calculated, with the viscosity of an ice-sucrose suspension being described by a power-law relationship of shear rate and crystal radius. The particle pile experiments show that mobilisation of equivalent crystal piles in magma chambers should be readily observed. As it is not observed, except in replenished magmatic systems, this suggests that the additional forces coming from cohesion and aggregation in crystal piles prevent mobilisation of magmatic crystals. The replenishment by secondary intrusions can lead to forces which overcome the strength of the pile.

We will explore progenitors extensively throughout this project. The progenitor, developed by Robert T Curtis, is a special type of infinite group formed by a semi-direct product of a free group m*n and a transitive permutation group of degree n. Since progenitors are infinite, we add necessary relations to produce finite homomorphic images. Curtis found that any non-abelian simple group is a homomorphic image of a progenitor of the form 2*n: N. In particular, we will investigate progenitors that generate two of the Mathieu sporadic groups, M11 and M11, as well as some classical groups. We will prove their existences a variety of different ways, including the process of double coset enumeration, Iwasawa's Lemma, and linear fractional mappings. We will also investigate the various techniques of finding finite images and their corresponding isomorphism types.

Rapakivi and alkali feldspar phenocrysts from the Deer Isle Granite Complex were investigated using in situ trace element and Pb isotope geochemistry to see if magma mixing or isothermal decompression was responsible for their formation. Pb isotope and trace element profiles, along with CL imagery of quartz phenocrysts, indicate compositional changes in the magma chamber occurred during rapakivi and alkali feldspar growth. Repeated episodes of magma mixing/replenishment by relatively isotopically primitive and LREE enriched magmas (along with hybridized variations with the host magma) created localized disequilibration. Alkali feldspar phenocrysts proximal to these zones of thermal perturbation were first resorbed and then mantled by plagioclase. Entrainment back into this zone of mixing caused multiple mantles of plagioclase to form on some phenocrysts. For grains more distal to these zones of mixing, complete disequilibration of the grain did not occur and continued growth resulted in Ba-rich alkali feldspar mantles over Ba-poor alkali feldspar cores. As crystallization of the chamber continued along solidification fronts, batches of cooler crystal-rich magmas settled en masse to the floor. Disaggregation of these batches during settling, and subsequent accumulation on the chambers floor, brought grains with disparate crystallization histories together. Filter pressing of the cumulate pile flushed highly evolved fluids out from interstitial pores. Small amounts of evolved liquid, trapped interstitial to the cumulate, formed LREE depleted albitic rims on some grains. The crystallization, transportation and juxtaposition of rapakivi, mantled alkali feldspar and plagioclase phenocrysts suggests that a relatively large and active magma chamber, periodically recharged by batches of melt, must have existed.

Hydrothermal systems at oceanic spreading centers play an important role in the composition of seawater, the formation of ore deposits, the support of microbial and macrofaunal ecosystems, and even for the development of life on early earth. These circulation systems are driven by heat transport from the underlying magma chamber, where latent heat of crystallization and sensible heat from cooling are transferred by vigorous, high Rayleigh number convection through a thin conductive boundary layer. The traditional study of magmatic-hydrothermal systems is primarily based on the time-series observation, which takes the form of repeat visits, continuous offline monitoring by autonomous instruments, or continuous online monitoring by instruments with satellite or cable links to shore. Although a number of studies have deployed autonomous monitoring instruments at vents and around mid-ocean ridges to investigate geophysical and hydrothermal processes, the data are still rather limited and a comprehensive understanding of magma-hydrothermal processes at oceanic spreading centers is lacking. Numerical modeling needs to be employed to elucidate the dynamic behavior of magmatic hydrothermal systems and for testing completing hypotheses in these complex, data-poor environments. In this dissertation, I develop a mathematical framework for investigating heat transport from a vigorously convecting, crystallizing, cooling, and replenished magma chamber to an overlying hydrothermal system at an oceanic spreading center. The resulting equations are solved numerically using MATLAB. The simulations proceed step-by-step to investigate several different aspects of the system. First, I consider a hydrothermal system driven by convection, cooling and crystallization in a ~ 100 m thick basaltic magma sill representing an axial magma chamber (AMC) at an oceanic spreading center. I investigate two different crystallization scenarios, crystal-suspended and crystal-settling, and consider both un-replenished and replenished AMCs. In cases without magma replenishment, the simulation results for crystals-suspended models show that heat output and the hydrothermal temperature decrease rapidly and crystallinity reaches 60% in less than ten years. In crystals-settling models, magma convection may last for decades, but decreasing heat output and hydrothermal temperatures still occur on decadal timescales. When magma replenishment is included, the magmatic heat flux approaches steady state on decadal timescales, while the magma body grows to double its original size. The rate of magma replenishment needed ranges between 5 x 10⁵ and 5 x 10⁶ m³/yr, which is somewhat faster than required for seafloor spreading, but less than fluxes to some terrestrial and subseafloor volcanoes on similar timescales. The heat output from a convecting, crystallizing, replenished magma body that is needed to drive observed high-temperature hydrothermal systems is consistent, with gabbro glacier models of crustal production at mid-ocean ridges. Secondly, I study the heat transfer model from a parametric perspective and examine the effects of both initial magma chamber thickness and magma replenishment rate on the hydrothermal heat output. The initial rate of convective heat transfer is independent of the initial sill thickness; but without magma replenishment, the rate of decay of the heat output varies linearly with thickness, resulting in short convective lifetimes and decaying hydrothermal temperatures for sills up to ~ 100m thick. When magma replenishment is included in crystals settling scenarios at constant or exponentially decreasing rates of ~ 10⁻⁸ m/s to the base of the sill, growth of the sill results in stabilized heat output and hydrothermal temperature on decadal timescales and a relatively constant to increasing thickness of the liquid layer. Sills initially ~ 10 m thick can grow, in principal, to ~ 10 times their initial size with stable heat output and a final melt thickness less than 100m. Seismic data provides evidence of AMC thickness, but it can not discriminate whether it denotes initial magma thickness or is a result of replenishment. These results suggest that magma replenishment might not be seismically detectable on decadal time scales. Periodic replenishment may also result in quasi-stable heat output, but the magnitude of the heat output may vary considerably in crystals suspended models at low frequencies; compared to crystals settling models. In these models the direct coupling between magmatic and hydrothermal heat output suggests that heat output fluctuations might be recorded in hydrothermal vents; but if damping effects of the basal conductive boundary layer and the upflow zone are taken into account, it seems unlikely that heat output fluctuations on a time scale of years would be recorded in hydrothermal vent temperatures or heat output. Thirdly, I extend the work to the binary system motivated by the fact that the real magmas are multi-component fluids. I focus on the extensively studied binary system, diopside-anorthite (Di-An), and investigate the effects of convection of a two-component magma system on the hydrothermal circulation system through the dynamic modeling of both temperature and heat output. I model the melt temperature and viscosity as a function of Di concentration, and incorporate these relations in the modeling of the heat flux. Simulations comparing the effects of different initial Di concentrations indicate that magmas with higher initial Di concentrations convect more vigorously, which results in faster heat transfer, more rapid removal of Di from the melt and growth of crystals on the floor. With magma replenishment, I assume that the magma chamber grows either horizontally or vertically. In either case magma replenishment at a constant rate of ~ 10⁻⁸ m³/a can maintain relatively stable heat output of 10⁷-10⁹ Watts and reasonable hydrothermal vent temperatures for decades. The final stabilized heat flux increases with increasing Di content of the added magma. Periodic replenishment with a 10 year period results in temperature perturbations within the magma that also increase as a function of increasing Di. With the simple magma model used here, one can not discern conclusively whether the decrease in magma temperature between the 1991/1992 and the 2005/2006 eruptions at EPR 9°50'N involved replenishment with more or less evolved magmas. Fourthly, I investigate a high-silica magma chamber as the hydrothermal circulation driver. I construct viscosity models for andesite and dacite melts as a function of temperature and water content and incorporate these expressions into a numerical model of thermal convective heat transport from a high Rayleigh number, well-mixed, crystallizing and replenished magma sill beneath a hydrothermal circulation system. Simulations comparing the time dependent heat flux from basalt, 0.1wt.% andesite, 3wt.% andesite, and 4wt.% dacite, indicate that higher viscosity magmas convect less vigorously, which results not only in lower heat transport and hydrothermal vent temperatures, but also in a lower decay rate of the vent temperature. Though somewhat colder, hydrothermal systems driven by unreplenished high-silica melts tend to have a longer lifetime than those driven by basalts, assuming a heat output cutoff of 10⁷ Watts. As in the basaltic case, magma replenishment at a rate of ~ 3 x 10⁵ - 3 x 10⁶ m³/a can maintain relatively stable heat output of 10⁷-10⁹ Watts and hydrothermal vent temperatures for decades. Idealized models of porous flow through the lower crust suggest such replenishment rates are not likely to occur, especially for high-viscosity magmas such as andesite and dacite. Long term stability of hydrothermal systems driven by these magmas requires an alternate means of magma replenishment. Finally, the dissertation concludes by discussing some avenues for future work. Most important of these are to: (1) couple magma convection with more realistic hydrothermal models and (2) link magma chamber processes to better physical models of replenishment and eruption.