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.     

Geological structure of the Curonian Spit (of the Baltic Sea) and its evolution history (revised)

The origin and evolution history of the Curonian Spit in the Baltic Sea are discussed based on geological and geomorphologic data. Evidence of the correlation between the spit formation stages and the Holocene sea level oscillations is presented. The ledges of moraine basement and marine accumulative features formed a single barrier. Its subsequent transformation occurred under the influence of the wave-related lithodynamic and eolian processes.     

Seasonal and solar cycle association of zonal drifts of ionospheric plasma irregularities in the Indian equatorial region

Long series of simultaneous VHF scintillation observations at two stations situated in near magnetic east-west direction in the vicinity of the dip equator in the Indian region have been employed to investigate the night-time ionospheric plasma zonal drifts. The drifts are found to be predominantly easterly. On comparing the magnitudes of the drifts with those results derived earlier by HF fading technique, monitoring signals from two satellites at a station and spaced receiver experiment, their associations with the season and the degree of solar activity are discussed. On a broader scale, the annual mean sunspot number is shown to have a direct control on the derived drift, the positive relationship even on day to day basis with the solar flux is established. However, the relationship, as understood by the slope of the best fit line, in the Indian region (0.27) is found to be weaker when compared with the similar slope (0.45) in the American sector. There appears to be no geomagnetic activity control on the estimated drifts.     

Tectonically driven fluid flow and associated low-grade metamorphism during the Alpine compression in the eastern Iberian Chain (Spain)

Fluid inclusions and clay mineralogy of the Permo-Triassic rocks from the Espina and Espadà Ranges (SE Iberian Chain, Spain) have been investigated to establish their relationship with hydrothermal fluid circulation during the Alpine Orogeny. Primary fluid inclusions in quartz-filled tension gashes in Permo-Triassic sandstones reveal maximum temperatures around 230 °C and very constant salinities of 8.5% wt. eq. NaCl. Secondary fluid inclusions found in quartz from the Santonian Ba–Cu–Hg deposits show similar compositional and thermodynamic characteristics, denoting an Alpine recrystallization. Clay mineral composition of Permo-Triassic mudrocks is characterized by pyrophyillite, indicating low-grade metamorphic conditions. Field observations and experimental data suggest that the crystallization of quartz in tension gashes, the formation of secondary fluid inclusions and the development of the metamorphism are contemporaneous and related to fluid circulation during the Alpine compression. Fluid flow took place along the Hercynian fault system that was reactivated during the Mesozoic rift stage and inverted during the Alpine deformation.     

A major element,PGE and Re–Os isotope study of Middle Atlas (Morocco) peridotite xenoliths: Evidence for coupled introduction of metasomatic sulphides and clinopyroxene

We present major element and PGE (platinum-group-element) abundances in addition to Re–Os isotope data for 11 spinel-facies whole rock peridotites from a single maar from the Middle Atlas Mountains, Morocco.Major element systematics of these xenoliths are generally correlated with indices of depletion. FeO–MgO systematics appear to suggest spinel-facies melting in the range of 5 to 25%. However, Al₂O₃ abundances in these xenoliths appear elevated relative to primitive mantle (Prima). The Al₂O₃ abundances in conjunction with other major elements require distinct re-enrichment of the Middle Atlas continental mantle root due to melt/rock reaction and precipitation of amphibole and/or clinopyroxene from passing silicate melts akin to MORB or OIB that evolved in reverse direction along the melting curves in e.g. FeO–MgO space. Sc and V confirm the range of apparent depletion and also indicate that the currently preserved fO₂ in these peridotites is distinctly different from fO₂ conditions observed in subduction zones.The majority of these xenoliths have low Os and Ir (I-PGEs) concentrations relative to Prima and modelled sulphide- and clinopyroxene-depleted residues of mantle melting under low fO₂, mid-ocean ridge-like conditions. Moreover, Pt and Pd (P-PGE) abundances are elevated when compared to their expected abundances after substantial melt extraction. Importantly, the systematically low Ir abundances in the majority of samples show well-correlated trends with Al₂O₃, MgO and Cu that are inconsistent with established melting trends. Os isotopes in the Middle Atlas xenoliths range from ¹⁸⁷Os/¹⁸⁸Os = 0.11604 to 0.12664 although most samples are close to chondritic. The Os isotope ratios are decoupled from ¹⁸⁷Re/¹⁸⁸Os but, together with Re abundances, also exhibit a good correlation with Al₂O₃, MgO and Cu.The major element, I-PGE and Os isotope correlations suggest that the initial melt depletion led to the exhaustion of sulphide and clinopyroxene (20 to 30%) without significant stabilization of I-PGE-rich alloys. During later modal metasomatism of the refractory Middle Atlas continental mantle root with silicate melts akin to MORB or OIB the introduction of clinopyroxene/amphibole reduced the volume of the melt inducing sulphur saturation in these melts causing precipitation of secondary sulphides. This coupled crystallization of pyroxenes and sulphides (chalcopyrite) resulted in the two-component mixing systematics exhibited by I-PGEs, Os isotopes with major elements and Cu preserved in the Middle Atlas continental mantle root.     

Yegorovite,Na<Subscript>4</Subscript>[Si<Subscript>4</Subscript>O<Subscript>8</Subscript>(OH)<Subscript>4</Subscript>]·7H<Subscript>2</Subscript>O,a new mineral from the Lovozero alkaline pluton,Kola Peninsula

A new mineral, yegorovite, has been identified in the late hydrothermal, low-temperature assemblage of the Palitra hyperalkaline pegmatite at Mt. Kedykverpakhk, Lovozero alkaline pluton, Kola Peninsula, Russia. The mineral is intimately associated with revdite and megacyclite, earlier natrosilite, microcline, and villiaumite. Yegorovite occurs as coarse, usually split prismatic (up to 0.05 × 0.15 × 1 mm) or lamellar (up to 0.05 × 0.7 × 0.8 mm) crystals. Polysynthetic twins and parallel intergrowths are typical. Mineral individuals are combined in bunches or chaotic groups (up to 2 mm); radial-lamellar clusters are less frequent. Yegorovite is colorless, transparent with vitreous luster. Cleavage is perfect parallel to (010) and (001). Fracture is splintery; crystals are readily split into acicular fragments. The Mohs hardness is ～2. Density is 1.90(2) g/cm³ (meas) and 1.92 g/cm³ (calc). Yegorovite is biaxial (?), with α = 1.474(2), β = 1.479(2), and γ = 1.482(2), 2V _meas > 70°, 2V _calc = 75°. The optical orientation is X ∧ a ～ 15°, Y = c, Z = b. The IR spectrum is given. The chemical composition determined using an electron microprobe (H₂O determined from total deficiency) is (wt %): 23.28 Na₂O, 45.45 SiO₂, 31.27 H₂O_calc; the total is 100.00. The empirical formula is Na_3.98Si_4.01O_8.02(OH)_3.98 · 7.205H₂O. The idealized formula is Na₄[Si₄O₈(OH)₄] · 7H₂O. Yegorovite is monoclinic, space group P2₁/c. The unit-cell dimensions are a = 9.874, b= 12.398, c = 14.897 Å, β = 104.68°, V = 1764.3 ų, Z = 4. The strongest reflections in the X-ray powder pattern (d, Å (I, %)([hkl]) are 7.21(70)[002], 6.21(72)[012, 020], 4.696(44)[022], 4.003(49)[211], 3.734(46)[(bar 2) 13], 3.116(100)[024, 040], 2.463(38)[(bar 4)02, (bar 2)43]. The crystal structure was studied by single-crystal method, R _hkl = 0.0745. Yegorovite is a representative of a new structural type. Its structure consists of single chains of Si tetrahedrons [Si₄O₈(OH)₄]∞ and sixfold polyhedrons of two types: [NaO(OH)₂(H₂O)₃] and [NaO(OH)(H₂O)₄] centered by Na. The mineral was named in memory of Yu. K. Yegorov-Tismenko (1938–2007), outstanding Russian crystallographer and crystallochemist. The type material of yegorovite has been deposited at the Fersman Mineralogical Museum of Russian Academy of Sciences, Moscow.     

Formation of accretion disks in close-binary systems with magnetic fields

We have developed a three-dimensional numerical model and applied it to simulate plasma flows in semi-detached binary systems whose accretor possesses a strong intrinsic magnetic field. The model is based on the assumption that the plasma dynamics are determined by the slow mean flow, which forms a backdrop for the rapid propagation of MHD waves. The equations describing the slow motion of matter were obtained by averaging over rapidly propagating pulsations. The numerical model includes the diffusion of magnetic field by current dissipation in turbulent vortices, magnetic buoyancy, and wave MHD turbulence. A modified three-dimensional, parallel, numerical code was used to simulate the flow structure in close binary systems with various accretor magnetic fields, from 10⁵ to 10⁸ G. The conditions for the formation of the accretion disk and the criteria distinguishing the two types of flow corresponding to intermediate polars and polars are discussed.     

Fragmentation and hydration of tektites and microtektites

Abstract— An examination of data collected over the last 30 years indicates that the percent of glass fragments vs. whole splash forms in the Cenozoic microtektite strewn fields increases towards the source crater (or source region). We propose that this is due to thermal stress produced when tektites and larger microtektites fall into water near the source crater while still relatively hot (>1150 °C). We also find evidence (low major oxide totals, frothing when melted) for hydration of most of the North American tektite fragments and microtektites found in marine sediments. High-temperature mass spectrometry indicates that these tektite fragments and microtektites contain up to 3.8 wt% H₂O. The H₂O-release behavior during the high-temperature mass-spectrometric analysis, plus high CI abundances (0.05 wt%), indicate that the North American tektite fragments and microtektites were hydrated in the marine environment (i.e., the H₂O was not trapped solely on quenching from a melt). The younger Ivory Coast and Australasian microtektites do not exhibit much evidence of hydration (at least not in excess of 0.5 wt% H₂O); this suggests that the degree of hydration increases with age. In addition, we find that some glass spherules (with <65 wt% SiO₂) from the upper Eocene clinopyroxene-bearing spherule layer in the Indian Ocean have palagonitized rims. These spherules appear to have been altered in a similar fashion to the splash form K/T boundary spherules. Thus, our data indicate that tektites and microtektites that generally contain >65 wt% SiO₂ can undergo simple hydration in the marine environment, while impact glasses (with <65 wt% SiO₂) can also undergo palagonitization.