Partial melting in an upwelling mantle column

Decompression melting of hot upwelling rock in the mantle creates a region of partial melt comprising a porous solid matrix through which magma rises buoyantly. Magma transport and the compensating matrix deformation are commonly described by two-phase compaction models, but melt production is less often incorporated. Melting is driven by the necessity to maintain thermodynamic equilibrium between mineral grains in the partial melt; the position and amount of partial melting that occur are thus thermodynamically determined. We present a consistent model for the ascent of a one-dimensional column of rock and provide solutions that reveal where and how much partial melting occurs, the positions of the boundaries of the partial melt being determined by conserving energy across them. Thermodynamic equilibrium of the boundary between partial melt and the solid lithosphere requires a boundary condition on the effective pressure (solid pressure minus melt pressure), which suggests that large effective stresses, and hence fracture, are likely to occur near the base of the lithosphere. Matrix compaction, melt separation and temperature in the partially molten region are all dependent on the effective pressure, a fact that can lead to interesting oscillatory boundary-layer structures. © 2008 The Royal Society

A catalogue of absorption-line systems in QSO spectra

We present a new catalog of absprption-line systems identified in the quasar spectra. It contains data on 821 QSOs and 8558 absorption systems comprizing 16139 absorption lines with measured redshifts in the QSO spectra. The catalog includes absorption-line systems consisting of lines of heavy elements, lines of neutral hydrogen, Lyman limit systems, damped Ly\alpha absorption systems, and broad absorption-line systems. The catalog is available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/412/707 and at www.ioffe.ru/astro/QC. Using the data of the present catalog we also discuss redshift distributions of absorption-line systems.Comment: 3 pages with 1 postscript figur

Englacial Pore Water Localizes Shear in Temperate Ice Stream Margins

The margins of fast‐moving ice streams are characterized by steep velocity gradients. Some of these gradients cannot be explained by a temperature‐dependent viscosity alone. Laboratory data suggest that water in the ice‐grain matrix decreases the ice viscosity; we propose that this causes the strong localization of shear in temperate ice stream margins. However, the magnitude of weakening and its consequences for ice stream dynamics are poorly understood. Here we investigate how the coupling between temperate ice properties, ice mechanics, and drainage of melt water from the ice stream margin alters the dynamics of ice streams. We consider the steady‐state ice flow, temperature, water content, and subglacial water drainage in an ice stream cross section. Temperate ice dynamics are modeled as a two‐phase flow, with gravity‐driven water transport in the pores of a viscously compacting and deforming ice matrix. We find that the dependence of ice viscosity on meltwater content focuses the temperate ice region and steepens the velocity gradients in the ice stream margin. It provides a possible explanation for the steep velocity gradients observed in some ice stream shear margins. This localizes heat dissipation there, which in turn increases the amount of meltwater delivered to the ice stream bed. This process is controlled by the permeability of the temperate ice and the sensitivity of ice viscosity to meltwater content, both of which are poorly constrained properties

Wave Structure and Velocity Profiles in Downwards Gas-Liquid Annular Flow

The downwards co-current gas-liquid annular flows inside a vertically oriented pipe have been experimentally investigated. The measurements and characterisation were performed using advanced optical non-intrusive laser-based techniques, namely Laser Induced Fluorescence, and Particle Image/Tracking Velocimetry. The investigated conditions were in the range of ReL = 306 – 1,532 and ReG = 0 – 84,600. Temporal film thickness time traces were constructed using the Laser Induced Fluorescence images. Based on these, the wave frequency was evaluated using direct wave counting approach and power spectral density analysis. Additionally, qualitative PIV observations revealed the presence of recirculation zones within a wave front of disturbance waves

Magmatic intrusions control Io's crustal thickness

Io, the most volcanically active body in the solar system, loses heat through eruptions of hot lava. Heat is supplied by tidal heating and is thought to be transferred through the mantle by magmatic segregation, a mode of transport that sets it apart from convecting terrestrial planets. We present a model that couples magmatic transport of tidal heat to the volcanic system in the crust, in order to determine the controls on crustal thickness, magmatic intrusions, and eruption rates. We demonstrate that magmatic intrusions are a key component of Io's crustal heat balance; around 80% of the magma delivered to the base of the crust must be emplaced and frozen as plutons to match rough estimates of crustal thickness. As magma ascends from a partially molten mantle into the crust, a decompacting boundary layer forms, which can explain inferred observations of a high-melt-fraction region.Comment: Accepted to JGR:Planets. 24 pages inc appendices and references. 7 figure

Designing a Green Roof for Ireland

A model is presented for the gravity-driven flow of rainwater descending through the soil layer of a green roof, treated as a porous medium on a at permeable surface representing an efficient drainage layer. A fully saturated zone is shown to occur. It is typically a thin layer, relative to the total soil thickness, and lies at the bottom of the soil layer. This provides a bottom boundary condition for the partially saturated upper zone. It is shown that after the onset of rainfall, well-defined fronts of water can descend through the soil layer. Also the rainwater flow is relatively quick compared with the moisture uptake by the roots of the plants in the roof. In a separate model the exchanges of water are described between the (smaller-scale) porous granules of soil, the roots and the rainwater in the inter-granule pores

The Variational Homoencoder: Learning to learn high capacity generative models from few examples

Hierarchical Bayesian methods can unify many related tasks (e.g. k-shot classification, conditional and unconditional generation) as inference within a single generative model. However, when this generative model is expressed as a powerful neural network such as a PixelCNN, we show that existing learning techniques typically fail to effectively use latent variables. To address this, we develop a modification of the Variational Autoencoder in which encoded observations are decoded to new elements from the same class. This technique, which we call a Variational Homoencoder (VHE), produces a hierarchical latent variable model which better utilises latent variables. We use the VHE framework to learn a hierarchical PixelCNN on the Omniglot dataset, which outperforms all existing models on test set likelihood and achieves strong performance on one-shot generation and classification tasks. We additionally validate the VHE on natural images from the YouTube Faces database. Finally, we develop extensions of the model that apply to richer dataset structures such as factorial and hierarchical categories.Comment: UAI 2018 oral presentatio

Moulin density controls drainage development beneath the Greenland ice sheet

Uncertainty remains about how the surface hydrology of the Greenland ice sheet influences its subglacial drainage system, affecting basal water pressures and ice velocities, particularly over intraseasonal and interseasonal timescales. Here we apply a high spatial (200m) and temporal (1h) resolution subglacial hydrological model to a marginal (extending ~25km inland), land-terminating, ~200km$^{2}$ domain in the Paakitsoq region, West Greenland. The model is based on that by Hewitt (2013) but adapted for use with both real topographic boundary conditions and calibrated modeled water inputs. The inputs consist of moulin hydrographs, calculated by a surface routing and lake-filling/draining model, which is forced with distributed runoff from a surface energy-balance model. Results suggest that the areal density of lake-bottom moulins and their timing of opening during the melt season strongly affects subglacial drainage system development. A higher moulin density causes an earlier onset of subglacial channelization (i.e., water transport through channels rather than the distributed sheet), which becomes relatively widespread across the bed, whereas a lower moulin density results in a later onset of channelization that becomes less widespread across the bed. In turn, moulin density has a strong control on spatial and temporal variations in subglacial water pressures, which will influence basal sliding rates, and thus ice motion. The density of active surface-to-bed connections should be considered alongside surface melt intensity and extent in future predictions of the ice sheet's dynamics.This work was funded through a UK Natural Environment Research Council Doctoral Training grant (LCAG/133), a Bowring Junior Research Fellowship (St Catharine's College, Cambridge), and a Leverhulme/Newton Trust Early Career Fellowship, all awarded to A.F.B. I.J.H. was supported by a Marie Curie FP7 Career Integration Grant within the 7th European Union Framework Programme

A dynamical systems approach to the tilted Bianchi models of solvable type

We use a dynamical systems approach to analyse the tilting spatially homogeneous Bianchi models of solvable type (e.g., types VI$_h$ and VII$_h$) with a perfect fluid and a linear barotropic $\gamma$-law equation of state. In particular, we study the late-time behaviour of tilted Bianchi models, with an emphasis on the existence of equilibrium points and their stability properties. We briefly discuss the tilting Bianchi type V models and the late-time asymptotic behaviour of irrotational Bianchi VII$_0$ models. We prove the important result that for non-inflationary Bianchi type VII$_h$ models vacuum plane-wave solutions are the only future attracting equilibrium points in the Bianchi type VII$_h$ invariant set. We then investigate the dynamics close to the plane-wave solutions in more detail, and discover some new features that arise in the dynamical behaviour of Bianchi cosmologies with the inclusion of tilt. We point out that in a tiny open set of parameter space in the type IV model (the loophole) there exists closed curves which act as attracting limit cycles. More interestingly, in the Bianchi type VII$_h$ models there is a bifurcation in which a set of equilibrium points turn into closed orbits. There is a region in which both sets of closed curves coexist, and it appears that for the type VII$_h$ models in this region the solution curves approach a compact surface which is topologically a torus.Comment: 29 page

Blawearie: a cairnfield excavation in a rock art landscape

Much of this work was done in the 1970s and 1980s when the subject, now increasingly regarded as mainstream within Neolithic studies, was largely shunned by professional archaeologists