Experimental study of liquid evolution in an Fe-rich,layered mafic intrusion: constraints of Fe-Ti oxide precipitation on the T-fO 2 and T-ϱ paths of tholeiitic magmas

The Newark Island layered intrusion, a composite intrusion displaying a similar fractionation sequence to the Skaergaard, has both dikes which preserved liquids fed into the intrusion and chilled pillows of liquids resident in the chamber. This study reports experimentally determined one atmosphere liquid lines of descent of these compositions as a function of oxygen fugacity which varies from QFM (quartz-fayalite-magnetite) to 0.5 log₁₀ units above IW (iron-wustite). These experiments reveal a strong oxygen fugacity dependence on the order of appearance and relative abundances of the Fe–Ti oxide minerals. Titanomagnetite saturates prior to ilmenite at QFM, but the order is reversed at lower oxygen fugacities. In the layered series of the Newark Island intrusion, ilmenite arrives shortly before titanomagnetite and the titanomagnetite/ilmenite ratio decreases monotonically after the cumulus appearance of titanomagnetite. Comparison of the crystallization sequence in the intrusion with that of the experiments requires that the oxygen fugacity in the intrusion increased relative to QFM before titanomagnetite saturation and decreased afterward, but always remained between the QFM and IW buffers. Similar trends in the modes of the Fe–Ti oxides (ilmenite and titanomagnetite) in the Skaergaard, Kiglapait, and Somerset Dam intrusions along with Fe₂O₃/FeO ratios in MORBs suggest that such a temperature-oxygen fugacity path may be typical of tholeiitic magma differentiation. Calculations of the temperature-density paths of the experimental liquids indicate that, at all possible oxygen fugacities, the density must have decreased abruptly after Fe–Ti oxide saturation. Accordingly, liquids replenishing the intrusion after Fe–Ti oxide saturation should pond at the bottom of the chamber, quenching against older cumulates. Field observation at the Newark Island intrusion confirm this prediction. The similarities in the fractionation paths of several other layered intrusions to that of the Newark Island intrusion suggest that the density of the liquids in these intrusions also decreased after Fe–Ti oxide saturation. Experiments on a suggested initial Skaergaard liquid are consistent with this model.     

Fractionation paths of Atka (Aleutians) high-alumina basalts: Constraints from phase relations

An Aleutian high-alumina basalt from the island of Atka at one atmosphere crystallizes plagioclase (1275°C) followed by olivine (1170°C) and clinopyroxene (1115°C). At oxygen fugacities along NNO, magnetite crystallizes below 1070°C, but its liquidus increases to at least 1175°C at an oxygen fugacity two log units above NNO. Phase relations at two kilobars pressure of melts containing small amounts of water are similar, although orthopyroxene and magnetite are observed to follow clinopyroxene. Amphibole crystallizes at near-liquidus temperatures only at water contents of melts approaching 4.5%. Amphibole assumes the liquidus in melts containing 5% water.Anhydrous melts crystallize plagioclase to 19 kbar, where garnet and clinopyroxene assume the liquidus. Olivine yields to clinopyroxene as the highest-temperature subliquidus phase at about 9 kbar.The array of compositions of basaltic Atka rocks, as displayed on appropriate pseudoternary projections, can be interpreted as a crystal fractionation path at moderate pressure (8 kbar) and small melt-water contents. The interpreted fractionating minerals are olivine, clinopyroxene, plagioclase, and (probably) magnetite. (The actual phenocrysts in Atka basalts like AT-1, which lacks phenocrystic clinopyroxene, must have crystallized at pressure less than 8 kbar, however.) The compositions of two-pyroxene andesites from Atka can be interpreted to lie on a lower-pressure fractionation trend at melt water contents of 2–3%. Such water contents are consistent with the complete absence of amphibole in any Atka rocks and are suggestive that water contents of the basaltic magmas, if the basalts are parental to the andesites, were 1–2%.     

Eclipsing Binary Stars in Nearby Galaxies

We briefly discuss the survey programme we are conducting to detect eclipsing binaries in local group galaxies. Some lightcurves from studies of M31, IC 1613 and NGC 6822 are presented along with details of future work.     

A precisely dated Proterozoic palaeomagnetic pole from the North China craton, and its relevance to palaeocontinental reconstruction

A palaeomagnetic pole position, derived from a precisely dated primary remanence, with minimal uncertainties due to secular variation and structural correction, has been obtained for China's largest dyke swarm, which trends for about 1000 km in a NNW direction across the North China craton. Positive palaeomagnetic contact tests on two dykes signify that the remanent magnetization is primary and formed during initial cooling of the intrusions. The age of one of these dykes, based on U–Pb dating of primary zircon, is 1769.1 ± 2.5 Ma. The mean palaeomagnetic direction for 19 dykes, after structural correction, is D  = 36°, I  = − 5°, k  = 63, α ₉₅ = 4°, yielding a palaeomagnetic pole at Plat=36°N, Plong=247°E, dp  = 2°, dm  = 4° and a palaeolatitude of 2.6°S. Comparison of this pole position with others of similar age from the Canadian Shield allows a continental reconstruction that is compatible with a more or less unchanged configuration of Laurentia, Siberia and the North China craton since about 1800 Ma     

Lens eddies in rotating two-fluid systems

Abstract

Laboratory experiments are conducted on the instantaneous release of a constant volume of intermediate density fluid along the interface of a two-layer fluid system in rigid body rotation about a vertical axis. Such a system leads to the development of a thin lens of fluid with anticyclonic motion which grows in radial extent to reach an equilibrium radius as the radial motion is captured by Coriolis effects, becomes unstable to non-symmetric dusturbances, and finally decays due to the action of fluid viscosity. Scaling and dimensional analysis arguments are advanced for the growth, equilibrium and decay phases of the motion. The scaling analysis for the initial growth and equilibrium phases are shown to be in good agreement with the experimental observations. The lens decay data do not collapse under the assumption of a simple Ekman spindown model assuming immiscible fluids. An empirical fit of the data for the decay phase is presented.     

Rotating flow over long shallow ridges

Abstract

The flow of a rotating homogeneous, incompressible fluid past a long ridge is investigated. An analysis is presented for flows in which E ? 1, Ro ～ E^½, H/D ～ E⁰, h/D ～ E^½ and cosα ～ E⁰ where E is the Ekman number, Ro the Rossby number, H/D the fluid depth to ridge width ratio, h/D the ridge height to ridge width ratio and α the angle between the free stream flow and a line perpendicular to the ridge axis. The analysis includes effects of the nonlinear inertial terms. Particular examples of a ridge of triangular cross section and a sinusoidal topography are investigated in some detail. Experiments are presented for a triangular ridge which are in good agreement with the theory.     

Sr isotopic evidence for a multi-source origin of the potassic magmas in the neapolitan area (S. Italy)

New Sr isotopic data on lavas and xenoliths from Somma-Vesuvius and other nearby volcanic areas (Phlegrean Fields and Ischia) are presented. Chemical and isotopic evidences show that not all the Phlegrean Fields rocks belong to the low K series, but some of them may be interpreted as low pressure differentiates of Somma magmas, i.e. as a part of the high K series. Two rock groups are defined in the Ischia low K series, which are well identified both in time and in chemical and isotopic features, and cannot be derived from the same magma source. The low K series in the studied area generally has lower Sr isotopic values than the high K series.Historical Vesuvian lavas show two distinct linear trends with negative slopes when⁸⁷Sr/⁸⁶Sr ratios are plotted against their ages of eruption. Such trends are interpreted to result from mixing of magmas in two separate reservoirs. Evidence from the Vesuvian ejecta shows that Somma-Vesuvius magmas underwent high or low pressure fractionation, in connection with different events of the Vesuvian activity. Distinct magma reservoirs developed episodically at different depths. Isotopic and geochemical evidences do not favour large scale assimilation of crustal materials by Somma-Vesuvius magmas, but instead appear to reflect mantle characteristics.A minimum of three different (inhomogeneous) source regions is necessary to account for the isotopic features of the studied rocks.     

Earth system curator: metadata infrastructure for climate modeling

The Earth System Curator is a National Science Foundation sponsored project developing a metadata formalism for describing the digital resources used in climate simulations. The primary motivating observation of the project is that a simulation/model’s source code plus the configuration parameters required for a model run are a compact representation of the dataset generated when the model is executed. The end goal of the project is a convergence of models and data where both resources are accessed uniformly from a single registry. In this paper we review the current metadata landscape of the climate modeling community, present our work on developing a metadata formalism for describing climate models, and reflect on technical challenges we have faced that require new research in the area of Earth Science Informatics.     

The effect of halogens on Zr diffusion and zircon dissolution in hydrous metaluminous granitic melts

Diffusion of Zr and zircon solubility in hydrous, containing approximately 4.5 wt% H₂O, metaluminous granitic melts with halogens, either 0.35 wt% Cl (LCl) or 1.2 wt% F (MRF), and in a halogen-free melt (LCO) were measured at 1.0 GPa and temperatures between 1,050 and 1,400 °C in a piston-cylinder apparatus using the zircon dissolution technique. Arrhenius equations for Zr diffusion in each hydrous melt composition are, for LCO with 4.4ǂ.4 wt% H₂O: % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavP1wzZbItLDhis9wBH5garm % Wu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9ir % Veeu0dXdh9vqqj-hEeeu0xXdbba9ev6pc9fs0-rqaqpepmKs4qpepe % I8kaL8kuc9pgc9q8qqaq-dhH6hb9hs0dXdHu6deP0u0-vr0-vr0db8 % meaabaqaciGacaGaaeaabaWaaeaaeaaakeaacqWGebarcqGH9aqpcq % aIYaGmcqGGUaGlcqaI4aaocqaI4aaocqGHXcqScqaIWaamcqGGUaGl % cqaIWaamcqaIZaWmcqWG4baEcqaIXaqmcqaIWaamdaahaaWcbeqaai % abgkHiTiabiIda4aaakiGbcwgaLjabcIha4jabcchaWnaabmaabaWa % aSaaaeaacqGHsislcqaIXaqmcqaI0aancqaIWaamcqGGUaGlcqaIXa % qmcqGHXcqScqaIZaWmcqaIZaWmcqGGUaGlcqaI5aqoaeaacqWGsbGu % cqWGubavaaaacaGLOaGaayzkaaaaaa!571F! D = 2.88 ±0.03x10 ^{- 8} exp( [( - 140.1 ±33.9)/(RT)] )D = 2.88 \pm 0.03x10^{ - 8} \exp \left( {{{ - 140.1 \pm 33.9} \over {RT}}} \right) , for LCl with 4.5ǂ.5 wt% H₂O: % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavP1wzZbItLDhis9wBH5garm % Wu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9ir % Veeu0dXdh9vqqj-hEeeu0xXdbba9ev6pc9fs0-rqaqpepmKs4qpepe % I8kaL8kuc9pgc9q8qqaq-dhH6hb9hs0dXdHu6deP0u0-vr0-vr0db8 % meaabaqaciGacaGaaeaabaWaaeaaeaaakeaacqWGebarcqGH9aqpcq % aIYaGmcqGGUaGlcqaIZaWmcqaIZaWmcqGHXcqScqaIWaamcqGGUaGl % cqaIWaamcqaI1aqncqWG4baEcqaIXaqmcqaIWaamdaahaaWcbeqaai % abgkHiTiabisda0aaakiGbcwgaLjabcIha4jabcchaWnaabmaabaWa % aSaaaeaacqGHsislcqaIYaGmcqaI1aqncqaI0aancqGGUaGlcqaI4a % aocqGHXcqScqaI2aGncqaI0aancqGGUaGlcqaIXaqmaeaacqWGsbGu % cqWGubavaaaacaGLOaGaayzkaaaaaa!5719! D = 2.33 ±0.05x10 ^{- 4} exp( [( - 254.8 ±64.1)/(RT)] )D = 2.33 \pm 0.05x10^{ - 4} \exp \left( {{{ - 254.8 \pm 64.1} \over {RT}}} \right) and for MRF with 4.9ǂ.3 wt% H₂O: % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavP1wzZbItLDhis9wBH5garm % Wu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9ir % Veeu0dXdh9vqqj-hEeeu0xXdbba9ev6pc9fs0-rqaqpepmKs4qpepe % I8kaL8kuc9pgc9q8qqaq-dhH6hb9hs0dXdHu6deP0u0-vr0-vr0db8 % meaabaqaciGacaGaaeaabaWaaeaaeaaakeaacqWGebarcqGH9aqpcq % aIYaGmcqGGUaGlcqaI1aqncqaI0aancqGHXcqScqaIWaamcqGGUaGl % cqaIWaamcqaIZaWmcqWG4baEcqaIXaqmcqaIWaamdaahaaWcbeqaai % abgkHiTiabiwda1aaakiGbcwgaLjabcIha4jabcchaWnaabmaabaWa % aSaaaeaacqGHsislcqaIYaGmcqaIYaGmcqaIZaWmcqGGUaGlcqaI4a % aocqGHXcqScqaIXaqmcqaI1aqncqGGUaGlcqaI1aqnaeaacqWGsbGu % cqWGubavaaaacaGLOaGaayzkaaaaaa!5715! D = 2.54 ±0.03x10 ^{- 5} exp( [( - 223.8 ±15.5)/(RT)] )D = 2.54 \pm 0.03x10^{ - 5} \exp \left( {{{ - 223.8 \pm 15.5} \over {RT}}} \right) . Solubilities determined by the dissolution technique were reversed for LCO +4.5ǂ.5 wt% H₂O by crystallization of a Zr-enriched glass of LCO composition at 1,200 and 1,050 °C at 1.0 GPa. The solubility data were used to calculate partition coefficients of Zr between zircon and hydrous melt, which are given by the following expressions: for LCO % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavP1wzZbItLDhis9wBH5garm % Wu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9ir % Veeu0dXdh9vqqj-hEeeu0xXdbba9ev6pc9fs0-rqaqpepmKs4qpepe % I8kaL8kuc9pgc9q8qqaq-dhH6hb9hs0dXdHu6deP0u0-vr0-vr0db8 % meaabaqaciGacaGaaeaabaWaaeaaeaaakeaacyGGSbaBcqGGUbGBcq % WGebardaqhaaWcbaGaemOwaOLaemOCaihabaGaemOEaONaemyAaKMa % emOCaiNaem4yamMaem4Ba8MaemOBa4Maei4la8IaemyBa0Maemyzau % MaemiBaWMaemiDaqhaaOGaeyypa0JaeGymaeJaeiOla4IaeGOnayJa % eG4mamZaaeWaaeaadaWcaaqaaiabigdaXiabicdaWiabicdaWiabic % daWiabicdaWaqaaiabdsfaubaaaiaawIcacaGLPaaacqGHsislcqaI % 1aqncqGGUaGlcqaI4aaocqaI3aWnaaa!5924! lnD_Zr^zircon/melt = 1.63( [10000/(T)] ) - 5.87\ln D_{Zr}^{zircon/melt} = 1.63\left( {{{10000} \over T}} \right) - 5.87 , for LCl % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavP1wzZbItLDhis9wBH5garm % Wu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9ir % Veeu0dXdh9vqqj-hEeeu0xXdbba9ev6pc9fs0-rqaqpepmKs4qpepe % I8kaL8kuc9pgc9q8qqaq-dhH6hb9hs0dXdHu6deP0u0-vr0-vr0db8 % meaabaqaciGacaGaaeaabaWaaeaaeaaakeaacyGGSbaBcqGGUbGBcq % WGebardaqhaaWcbaGaemOwaOLaemOCaihabaGaemOEaONaemyAaKMa % emOCaiNaem4yamMaem4Ba8MaemOBa4Maei4la8IaemyBa0Maemyzau % MaemiBaWMaemiDaqhaaOGaeyypa0JaeGymaeJaeiOla4IaeGinaqJa % eG4naCZaaeWaaeaadaWcaaqaaiabigdaXiabicdaWiabicdaWiabic % daWiabicdaWaqaaiabdsfaubaaaiaawIcacaGLPaaacqGHsislcqaI % 0aancqGGUaGlcqaI3aWncqaI1aqnaaa!5920! lnD_Zr^zircon/melt = 1.47( [10000/(T)] ) - 4.75\ln D_{Zr}^{zircon/melt} = 1.47\left( {{{10000} \over T}} \right) - 4.75 and, for MRF by % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavP1wzZbItLDhis9wBH5garm % Wu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9ir % Veeu0dXdh9vqqj-hEeeu0xXdbba9ev6pc9fs0-rqaqpepmKs4qpepe % I8kaL8kuc9pgc9q8qqaq-dhH6hb9hs0dXdHu6deP0u0-vr0-vr0db8 % meaabaqaciGacaGaaeaabaWaaeaaeaaakeaacyGGSbaBcqGGUbGBcq % WGebardaqhaaWcbaGaemOwaOLaemOCaihabaGaemOEaONaemyAaKMa % emOCaiNaem4yamMaem4Ba8MaemOBa4Maei4la8IaemyBa0Maemyzau % MaemiBaWMaemiDaqhaaOGaeyypa0JaeGymaeJaeiOla4IaeGinaqJa % eG4naCZaaeWaaeaadaWcaaqaaiabigdaXiabicdaWiabicdaWiabic % daWiabicdaWaqaaiabdsfaubaaaiaawIcacaGLPaaacqGHsislcqaI % 0aancqGGUaGlcqaI5aqocqaIXaqmaaa!591C! lnD_Zr^zircon/melt = 1.47( [10000/(T)] ) - 4.91\ln D_{Zr}^{zircon/melt} = 1.47\left( {{{10000} \over T}} \right) - 4.91 . Experiments on the same compositions, but with water contents down to 0.5 wt%, demonstrated reductions in both the diffusion coefficient of Zr and zircon solubility in the melt. The addition of halogens at the concentration levels studied to metaluminous melts has a small effect on either the diffusion of Zr in the melt, or the solubility of zircon at all water concentrations and temperatures investigated. At 800 °C, the calculated diffusion coefficient of Zr is lowest in LCl, 9᎒^-17 m² s^-1, and is highest in LCO, 4᎒^-15 m² s^-1. Extrapolation of the halogen-free solubility data to a magmatic temperature of 800 °C yields solubilities of approximately one-third of those directly measured in similar compositions, predicted by earlier studies of zircon dissolution and based upon analyses of natural rocks. This discrepancy is attributed to the higher oxygen fugacity of the experiments of this study compared with previous studies and nature, and the effect of oxygen fugacity on the structural role of iron in the melt, which, in turn, affects zircon solubility, but does not significantly affect Zr diffusion.     

Comparison of streamflow recession between plantations and native forests in small catchments in Central-Southern Chile

In central Chile, many communities rely on water obtained from small catchments in the coastal mountains. Water security for these communities is most vulnerable during the summer dry season and, from 2010 to 2017, rainfall during the dry season was between 20% and 40% below the long-term average. The rate of decrease in stream flow after a rainfall event is a good measure of the risk of flow decreasing below a critical threshold. This risk of low flow can be quantified using a recession coefficient (α) that is the slope of an exponential decay function relating flow to time since rainfall. A mathematical model was used to estimate the recession coefficient (α) for 142 rainstorm events (64 in summer; 78 in winter) in eight monitored catchments between 2008 and 2017. These catchments all have a similar geology and extend from 35 to 39 degrees of latitude south in the coastal range of south-central Chile. A hierarchical cluster analysis was used to test for differences between the mean value of α for different regions and forest types in winter and summer. The value of α did not differ (p < 0.05) between catchments in winter. Some differences were observed during summer and these were attributed to morphological differences between catchments and, in the northernmost catchments, the effect of land cover (native forest and plantation). Moreover, α for catchments with native forest was similar to those with pine plantations, although there was no difference (p < 0.05) between these and Eucalyptus plantations. The recession constant is a well-established method for understanding the effect of climate and disturbance on low flows and baseflows and can enhance local and regional analyses of hydrological processes. Understanding the recession of flow after rainfall in small headwater catchments, especially during summer, is vital for water resources management in areas where the establishment of plantations has occurred in a drying climate.