Pre-eruptive magma mixing in ash-flow deposits of the Tertiary Rum Igneous Centre,Scotland

The Northern Marginal Zone of the Rum Igneous Centre is a remnant of an early caldera and its infill. It is composed of intra-caldera breccias and various small-volume pyroclastic deposits, overlain by prominent rhyodacite ash-flow sheets of up to 100 m thickness. The ash-flows were fed from a feeder system near the caldera ring-fault, and intrusive rhyodacite can locally be seen grading into extrusive deposits. A variety of features suggest that the ash-flows were erupted from a magma chamber that contemporaneously hosted felsic and mafic magmas: (i) chilled basaltic inclusions in rhyodacite; (ii) formerly glassy basaltic to andesitic enclaves with fluid-fluid relationships; (iii) feldspars with thick reaction rims enclosed in the basaltic to andesitic inclusions, yet with cores chemically resembling those of the rhyodacite: (iv) trace element compositions of the rhyodacite and the mafic enclaves form a mixing line between the end-member rhyodacite and basalt compositions. Additionally, textural and chemical features in the rhyodacite feldspar phenocrysts are consistent with magma mixing; (v) feldspars with resorption embayments cutting through internal zonation of the crystals; (vi) complexly zoned crystals with sieve-textured zones that are overgrown with euhedral zones; (vii) oscillatory zonation of feldspar phenocrysts in the rhyodacite, showing sharp increases in anorthite (An 10%) followed by gradual decrease in An-content (An 4%). This evidence points to eruption of ash-flows from a felsic magma chamber that was periodically replenished by injection of mafic magma. Diffusional mixing between the two magmas was permitted by temperature and compositional differences, but was slow due to the contrast in viscosities and densities. The Fe–Ti–P-enriched basic magma that replenished the chamber was degassing on entering the lower temperature environment and soon equilibrated thermally, followed by chemical exchange between the two end-member magmas. This process formed hybrid andesite enclaves enriched in trace elements beyond that caused by simple mixing, implying trace element diffusion in addition to bulk mixing. Eruption was caused by replenishment with, and degassing of, the basic magma and the chamber partially evacuated while the process of hybridisation was underway. The erupted products record magma mixing by chamber replenishment, blending of two magmas and elemental exchange in the magma chamber, and also physical mingling in the eruptive conduit.     

An element formulation for perturbed motion about the center of mass

The perturbed motion of a rigid body about its center of mass, is formulated in terms of the six elements:l, the magnitude of the angular momentum vector;h, the total energy; and , two linear functions of the independent variable; and ₁ and ₁, two Euler angles that orientate the inertial frame with respect to the unperturbed solution. Solutions from the element formulation and the original Euler equations are numerically compared using shuttle-type data. For applied torques smaller than a given magnitude, the element formulation produced the following results: (1) larger step sizes in the numerical integration of the differential equations, resulting in an overall computational time-saving, and (2) more significant figures of accuracy in the computation of the variables describing the state of the rigid body.     

On the nature of rapidly deepening Canadian East Coast winter storms

Abstract

During the Canadian Atlantic Storms Program (CASP) field project from 15 January to 15 March 1986, seven storms underwent rapid deepening either within or adjacent to the observational network. The most rapid deepening phase of these storms was associated with the surface low‐pressure centre being located close to to the rain‐snow boundary and the end of this phase was associated with the low centre eventually moving to sub‐freezing surface temperatures. Precipitation bands of both rain and snow occurred during maximum deepening. Observations are consistent with divergence and subsidence occurring near the centre itself. It is suggested that a mesoscale circulation initiated by melting snow significantly affects the deepening process.     

Recent effects of tidal and hydro‐meteorological changes on coastal plains near the mouth of the Amazon River

The goal of this work was to understand the main hydrodynamic processes acting on tidal flats of the coast of Amapá near the mouth of the Amazon River, and how they change over the short term (~ 20 years). The analysis of morphological and geobotanical units was carried out by applying processing and interpretation methods to optical and synthetic aperture radar (SAR) images, combined with data on water salinity, maximum flood height, sedimentary facies data, rainfall and river discharge. The temporal analysis of morphological and geobotanical units suggests the relative stabilization of savannah, ‘várzea’ and mangrove areas during the drier period and increasing tidal amplitude between 1987 to 1997. The wetter period and decreasing tidal amplitude between 1997 to 2008 led to an increase in the area of ‘várzea’ and lakes over savannah, and the expansion of mangroves mainly over the inundated field and tidal mud/mixed flats. Therefore, the decrease in rainfall index during the drier period is well‐correlated with the reduction of the Calçoene River discharge and jointly with increasing tidal amplitude favored the increase of migration rate of the mud bank and erosion profile along the littoral. It was followed by the increase of the Calçoene River discharge and jointly with decreasing tidal amplitude during the wetter period, favoring the development of mangroves on muddy substrates near the coastline. Copyright © 2013 John Wiley & Sons, Ltd.     

Multi‐transient electromagnetic repeatability experiment over the North Sea Harding field‡

We present results of synthetic time‐lapse and real repeatability multi‐transient electromagnetic surveys over the North Sea Harding field. Using Archie's law to convert porosity and fluid saturation to resistivity we created 3D isotropic models of the reservoir resistivity at different stages of production from the initial state in 1996 through to complete hydrocarbon production by 2016 and, for each stage, we simulated an east‐west transient electromagnetic survey line across Harding. Unconstrained 1D full‐waveform Occam inversions of these synthetic data show that Harding should be detectable and its lateral extent reasonably well‐defined. Resistivity changes caused by hydrocarbon production from initial pre‐production state to production of the oil rim in 2011 are discernible as are significant changes from 2011–2016 during the modelled gas blowdown phase. The 2D repeatability surveys of 2007 and 2008 tied two wells: one on and the other off the structure. Between the two surveys the segment of the field under investigation produced 3.9 million barrels of oil – not enough to generate an observable time‐lapse electromagnetic anomaly with a signal‐to‐noise ratio of 40 dB. Processing of the 2007 and 2008 data included deconvolution for the measured source current and removal of spatially‐correlated noise, which increased the signal‐to‐noise ratio of the recovered impulse responses by about 20 dB and resulted in a normalized root‐mean‐square difference of 3.9% between the data sets. 1D full‐waveform Occam inversions of the real data showed that Harding was detectable and its lateral extent was also reasonably well‐defined. The results indicate that the multi‐transient electromagnetic method is suitable for exploration, appraisal and monitoring hydrocarbon production.     

An experimental investigation of olivine morphology

Olivine crystals can adopt ten types of shape. Experimental crystallization of eight rock melts shows that there is a systematic change from polyhedral or granular olivines hopper olivines branching olivines randomly oriented chain olivines parallel-growth chain olivines chain+lattice olivines plate or feather olivines, with increase in cooling rate and with increase in degree of supercooling. This sequence involves changes from complete to progressively less complete crystals and from equant habit to elongate bladed habit (c>ab) to tabular habit (ac b). The sequence is not affected by the phase relations of the melt. The larger the olivine content of a melt the slower the cooling rate at which a particular olivine shape grows, whereas the lower the melt viscosity, the greater the cooling rate. Irrespective of the melt composition, comparable crystal shapes grow at the same degrees of supercooling. By comparison of the shapes of olivine crystals in experiments with those in rocks of similar composition, it is possible to deduce the cooling rate through the olivine crystallization interval and the approximate degree of supercooling at which the olivine crystals nucleated and grew in the rocks. The various shapes of skeletal olivines in many picrites, olivine-rich basalts and the Archaean spinifex rocks are not due to rapid cooling, but to rapid olivine growth caused by the high normative olivine content of the magma.     

Novel basin modelling concept for simulating deformation from mechanical compaction using level sets

As sedimentation progresses in the formation and evolution of a depositional geologic basin, the rock strata are subject to various stresses. With increasing lithostatic pressure, compressional forces act to compact the porous rock matrix, leading to overpressure buildup, changes in the fluid pore pressure and fluid flow. In the context of petroleum systems modelling, the present study concerns the geometry changes that a compacting basin experiences subject to deposition. The purpose is to track the positions of the rock layer interfaces as compaction occurs. To handle the challenge of potentially large geometry deformations, a new modelling concept is proposed that couples the pore pressure equation with a level set method to determine the movement of lithostratigraphic interfaces. The level set method propagates an interface according to a prescribed speed. The coupling term for the pore pressure and level-set equations consists of this speed function, which is dependent on the compaction law. The two primary features of this approach are the simplicity of the grid and the flexibility of the speed function. A first evaluation of the model concept is presented based on an implementation for one spatial dimension accounting for vertical effective stress. Isothermal conditions with a constant fluid density and viscosity were assumed. The accuracy of the implemented numerical solution for the case of a single stratigraphic unit with a linear compaction law was compared to the available analytical solution [38]. The multi-layer setup and the nonlinear case were tested for plausibility.     

Ultramafic inclusions in anorthite megacrysts from the Isle of Skye

1–15 μm diameter ultramafic inclusions in anorthitic plagioclase megacrysts from Skye were analyzed by electron microprobe and tentatively identified as sub-calcic augite, or sub-calcic augite submicroscopically intergrown with amphibole. The inclusions nucleated in Ca-Al-depleted, Fe-Mg-H₂O-enriched magma adjacent to the interface between plagioclase crystals and basalt and were rapidly enveloped by the growing plagioclase. The inclusions may be important in the fractional crystallization of the Skye alkali olivine basalt suite.     

OB association members in the ACT and TRC catalogues

The Hipparcos Catalogue contains members of nearby OB associations brighter than 12th magnitude in V . However, membership lists are complete only to magnitude V =7.3. In this paper we discuss whether proper motions listed in the 'Astrographic Catalogue+Tycho' reference catalogue (ACT) and the Tycho Reference Catalogue (TRC), which are complete to V ∼10.5 mag, can be used to find additional association members. Proper motions in the ACT/TRC have an average accuracy of ∼3 mas yr⁻¹. We search for ACT/TRC stars which have proper motions consistent with the spatial velocity of the Hipparcos members of the nearby OB associations already identified by de Zeeuw et al. These stars are first selected using a convergent-point method, and then subjected to further constraints on the proper-motion distribution, magnitude and colour to narrow down the final number of candidate members. Monte Carlo simulations show that the proper-motion distribution, magnitude, and colour constraints remove ∼97 per cent of the field stars, while at the same time retain more than 90 per cent of the cluster stars.
The procedure has been applied to five nearby associations: the three subgroups of Sco OB2, plus Per OB3 and Cep OB6. In all cases except Cep OB6, we find evidence for new association members fainter than the completeness limit of the Hipparcos Catalogue. However, narrow-band photometry and/or radial velocities are needed to pinpoint the cluster members, and to study their physical characteristics.