Quantitative Simulation of the Hydrothermal Systems of Crystallizing Magmas on the Basis of Transport Theory and Oxygen Isotope Data: An analysis of the Skaergaard Intrusion

Application of the principles of transport theory to studiesof magma-hydrothermal systems permits quantitative predictionsto be made of the consequences of magma intruding into permeablerocks. Transport processes which redistribute energy, mass,and momentum in these environments can be represented by a setof partial differential equations involving the rate of changeof extensive properties in the system. Numerical approximationand computer evaluation of the transport equations effectivelysimulates the crystallization of magma, cooling of the igneousrocks, advection of chemical components, and chemical and isotopicmass transfer between minerals and aqueous solution. Numerical modeling of the deep portions of the Skaergaard magma-hydrothermalsystem has produced detailed maps of the temperature, pressure,fluid velocity, integrated fluid flux, ¹⁸O-values in rock andfluid, and extent of nonequilibrium exchange reactions betweenfluid and rock as a function of time for a two-dimensional cross-sectionthrough the pluton. An excellent match was made between calculated¹⁸O-values and the measured ¹⁸O-values in the three principalrock units, basalt, gabbro, and gneiss, as well as in xenolithsof roof rocks that are now embedded in Layered Series; the latterwere evidently depleted in ¹⁸O early in the system's coolinghistory, prior to falling to the bottom of the magma chamber.The best match was realized for a system in which the bulk rockpermeabilities were 10^–13 cm² for the intrusion, 10^–11cm² for basalt, and 10^–16 cm² for gneiss; reaction domainsizes were 0.2 cm in the intrusion and gneiss and 0.01 cm inthe basalts, and activation energy for the isotope exchangereaction between fluid and plagioclase was 30 kcal/mole. The calculated thermal history of the Skaergaard system wascharacterized by extensive fluid circulation that was largelyrestricted to the permeable basalts and to regions of the plutonstratigraphically above the basalt-gneiss unconformity. Althoughfluids circulated all around the crystallizing magma, fluidflow paths were deflected around the magma sheet during theinitial 130,000 years. At that time, crystallization of thefinal sheet of magma and fracture of the rock shifted the circulationsystem toward the center of the intrusion, thereby minimizingthe extent of isotope exchange between rocks near the marginof the intrusion at this level. For comparison, similar calculationswere also made for pure conductive cooling; it was found thatthe rate of crystallization of the magma body was not changed.The solidified pluton cooled by a factor of about 2 faster inthe presence of a hydrothermal system. Transport rates of thermal energy out of the intrusion and oflow-¹⁸O fluids into the intrusion controlled the overall isotopeexchange process. During the initial 150,000 years, temperatureswere high and reaction rates were fast; thus, fluids flowinginto the intrusion quickly equilibrated with plagioclase. However,the temperature decreased between 120,000 and 175,000 yearsand caused a decrease in reaction rates and an increase in theequilibrium fractionation factor between plagioclase and fluid.Consequently, during this time period fluids in the intrusiontended to be out of equilibrium with plagioclase. After 175,000years temperatures had decreased sufficiently that reactionrates became insignificant, but convection rates were largeenough to redistribute fluid and enlarge the regions where fluidand plagioclase were out of equilibrium. By 400,000 years, thepluton had cooled to approximately ambient temperatures, andthe final ¹⁸O values were ‘frozen in’. Reactionsbetween hydrothermal fluid and the intrusion occurred over abroad range in temperature, 1000-200 °C, but 75 per centof the fluid circulated through the intrusion while its averagetemperature was >480 °C. This relatively high temperatureis consistent with the observation that only minor amounts ofhydrothermal alteration products were formed in the naturalsystem, even where several per mil shifts in ¹⁸O were detected. The relative quantities of fluid to rock integrated over theentire cooling history were 0.52 for the upper part of intrusion,0.88 for the basalt, 0.003 for the gneiss, and 0.41 for theentire domain. Almost all of the fluid flowed into the intrusionfrom the basalt host rocks that occur adjacent to the side contactsof the intrusion. Convection transferred about 20 per cent ofthe total heat contained in the gabbro upward into the overlyingbasalts; the remaining 80 per cent of the heat was transferredby conduction.     

Convection forced by a descending dry layer and low-level moist convergence

A narrow line of convective showers was observed over southern England on 18 July 2005 during the Convective Storm Initiation Project (CSIP). The showers formed behind a cold front (CF), beneath two apparently descending dry layers (i.e. sloping so that they descended relative to the instruments observing them). The lowermost dry layer was associated with a tropopause fold from a depression, which formed 2 d earlier from a breaking Rossby wave, located northwest of the UK. The uppermost dry layer had fragmented from the original streamer due to rotation around the depression (This rotation was also responsible for the observations of apparent descent—ascent would otherwise be seen behind a CF). The lowermost dry layer descended over the UK and overran higher  θ _w   air beneath it, resulting in potential instability. Combined with a surface convergence line (which triggered the convection but had less impact on the convective available potential energy than the potential instability), convection was forced up to 5.5 km where the uppermost dry layer capped it. The period when convection was possible was very short, thus explaining the narrowness of the shower band. Convective Storm Initiation Project observations and model data are presented to illustrate the unique processes in this case.     

The 13th and 14th Workshops on Antarctic Meteorology and Climate

<正>1.Overview In July 2018, the Antarctic community came together to meet at the 13th Workshop on Antarctic Meteorology and Climate (WAMC) in Madison, Wisconsin, USA (Fig. 1); and in the following year in June 2019, the 14th WAMC was held in     

Paleoecology of the Mollusca of the Tjörnes sequence, Iceland

A long sequence of open-coast tidal flat deposits with a restricted, mainly boreal fauna, comprises the Tapes and Mactra Zones. Characteristic species are Venerupis aurea, Spisula arcuata, Arctic islandica and Lentidium complanatum . This stage is provisionally correlated with the Upper Pliocene of the North Sea Basin. During the time of deposition of the Serripes Zone the basin apparently deepened and the fauna greatly diversified with immigration of arctic elements. It went back finally to beach-deposit conditions. The characteristic species — Serripes groenlandicus, Macoma praetenius , and Beta borealis — are also characteristic of the Early Pleistocene of the North Sea Basin with which the Serripes Zone is provisionally correlated. The climate is interpreted as having been colder than that of the Tapes and Mactra zones. A preliminary list of the Mollusca is given.     

Evolution of fluid pressure and fracture propagation during contact metamorphism

Abstract Rock fracture enhances permeability and provides pathways through which fluids migrate. During contact metamorphism, fluids contained in isolated pores and fractures expand in response to temperature increases caused by the dissipation of heat from magmas. Heat transport calculations and thermomechanical properties of water-rich fluids demonstrate (1) that thermal energy is a viable mechanism to produce and maintain pore fluid pressure (P_f) in a contact metamorphic aureole; (2) that the magnitude of P_f generated is sufficient to propagate fractures during the prograde thermal history (cause hydrofracture) and enhance permeability; and (3) that P_f-driven fracture propagation is episodic with time-scales ranging from years to thousands of years. Because P_f dissipation is orders of magnitude faster than P, _f buildup, P_f oscillations and cyclical behaviour are generated as thermal heating continues. The P_f cycle amplitude depends on the initial fracture length, geometry and the rock's resistance to failure whereas the frequency of fracture depends on the rate of heating. Consequently, oscillation frequency also varies spatially with distance from the heat source. Time series of fluid pressures caused by this process suggest that cyclical fracture events are restricted to an early time period of the prograde thermal event near the intrusive contact. In the far field, however, individual fracture events have a lower frequency but continue to occur over a longer time interval. Numerous fracture cycles are possible within a single thermal event. This provides a provisional explanation for multiple generations of veins observed in outcrop. P _f cycling and oscillations may explain several petrological features. If pore fluids are trapped at various positions along a pressure cycle, the large amplitude of P_f variations for small fractures may account for different pressures recorded by fluid inclusions analysed from a single sample. P_f oscillations, during a single thermal episode, also drive chemical reactions which can produce complex mineral textures and assemblages for discontinuous reactions and/or zoning patterns for continuous reactions. These can mimic polymetamorphic or disequilibrium features. Temporal aspects of fracture propagation and permeability enhancement also constrain the likely timing of fluid flow and fluid-mineral interactions. These data suggest that fluid flow and fluid-mineral reactions are likely to be restricted to an early period in the prograde thermal history, characterized by high P_f coincident with relatively high temperatures, fracture propagation and consequent increases in permeability. This early prograde hydration event is followed by diffusional peak metamorphic reactions. This relationship is evident in the complex mineralogical textures common in some metamorphosed rocks.     

Magma Storage Region Processes Inferred from Geochemistry of Fe-Ti Oxides in Andesitic Magma, Soufriere Hills Volcano, Montserrat, W.I.

Analyses of Fe–Ti oxides help constrain models of magmastorage region processes for the Soufrière Hills Volcano,Montserrat (W.I.), and provide clear evidence of the natureof transient heating events in the magma storage region. Toconstrain timescales of magma heating and remobilization, theTiO₂ zoning patterns in a time series of natural titanomagnetiteswere compared with those produced in controlled phase equilibriumexperiments on the andesite. Most samples of andesite eruptedfrom 1995 to 2002 contain titanomagnetite crystals with uniformcore compositions (TiO₂