LavAtmos: An open-source chemical equilibrium vaporization code for lava worlds

To date, over 500 short-period rocky planets with equilibrium temperatures above 1500 K have been discovered. Such planets are expected to support magma oceans, providing a direct interface between the interior and the atmosphere. This provides a unique opportunity to gain insight into their interior compositions through atmospheric observations. A key process in doing such work is the vapor outgassing from the lava surface. LavAtmos is an open-source code that calculates the equilibrium chemical composition of vapor above a dry melt for a given composition and temperature. Results show that the produced output is in good agreement with the partial pressures obtained from experimental laboratory data as well as with other similar codes from literature. LavAtmos allows for the modeling of vaporization of a wide range of different mantle compositions of hot rocky exoplanets. In combination with atmospheric chemistry codes, this enables the characterization of interior compositions through atmospheric signatures.     

Power law or power flaw?

Since its introduction by Svensson in 1959, the power law curve y = ax^b (where x and y are horizontal and vertical direction, respectively) has been widely used in morphological analysis of glacial trough cross-profiles. The numerical constants a and b are obtained by a linear regression analysis of the logarithmic form of the power law curve (ln y = ln a + b ln x). The value b then gives a measure for the form of the cross-section. However, over the years this form of the power law has endured a lot of criticism. This criticism is well founded, since this particular form of the power law is not suitable for curve fitting in morphological analyses. In this paper a general power law is proposed, of the form y − y₀ = a|x − x₀|^b (where x₀, y₀ are the coordinates of the origin of the cross-profile). A unique and unbiased solution for this equation is obtained with a general least-squares method, thereby minimizing the error between the cross-profile data and the curve, and not between the logarithmic transform of the data and its regression line. This provides a robust way to characterize trough cross-section forms. © 1998 John Wiley & Sons, Ltd.     

Tomographic images of the upper mantle below central Europe and the Mediterranean

Results from delay time tomography of the European-Mediterranean upper mantle are discussed and where possible interpreted in terms of geodynamic processes. Slab-like positive velocity anomalies of which the locations correlate well with deeper seismicity are found beneath Spain, the Tyrrhenian basin, and the Aegean. These structures are interpreted as images of subducted slabs. Large aseismic regions with positive velocity anomalies are found beneath the Western Mediterranean, Italy, the Alps, Dinarides, the Pannonian basin, northern Greece, and the Aegean. These anomalies can also be linked to subducted lithosphere. From the anomaly patterns it is deduced that subduction occurred below the Western Mediterranean and along both sides of the Adriatic micro-plate. Beneath the Dinarides and northern Greece the velocity structures suggest detachment of the slab from the surface.     

Trophic Links in the Plankton in the Low Salinity Zone of a Large Temperate Estuary: Top-down Effects of Introduced Copepods

We investigated trophic relationships involving microzooplankton in the low salinity zone of the San Francisco Estuary (SFE) as part of a larger effort aimed at understanding the dynamics of the food web supporting the endangered delta smelt, Hypomesus transpacificus. We performed 14 cascade experiments in which we manipulated the biomass of a copepod (Limnoithona tetraspina, Pseudodiaptomus forbesi, or Acartiella sinensis) and quantified responses of lower trophic levels including bacterioplankton, phytoplankton, and microzooplankton. Microzooplankton comprised a major food source for copepods; 9 out of 14 experiments showed removal of at least one group of microzooplankton by copepods. In contrast, the impact of copepods on phytoplankton was indirect; increased copepod biomass led to greater growth of phytoplankton in 3 of 14 experiments. Estimated clearance rates on microzooplankton were 4 mL day^?1 for L. tetraspina and 2–6 mL day^?1 for P. forbesi, whereas A. sinensis consumed mainly copepod nauplii. Complex trophic interactions, including omnivory, among copepods, microzooplankton, and different components of the phytoplankton likely obscured clear trends. The food web of the SFE is probably less efficient than previously thought, providing poor support to higher trophic levels; this inefficient food web is almost certainly implicated in the continuing low abundance of fishes, including the delta smelt that use the low salinity zone of the San Francisco Estuary.     

Effects of vegetation on flow and sediment transport: comparative analyses and validation of predicting models

The presence of vegetation modifies flow and sediment transport in alluvial channels and hence the morphological evolution of river systems. Plants increase the local roughness, modify flow patterns and provide additional drag, decreasing the bed‐shear stress and enhancing local sediment deposition. For this, it is important to take into account the presence of vegetation in morphodynamic modelling. Models describing the effects of vegetation on water flow and sediment transport already exist, but comparative analyses and validations on extensive datasets are still lacking. In order to provide practical information for modelling purposes, we analysed the performance of a large number of models on flow resistance, vegetation drag, vertical velocity profiles and bed‐shear stresses in vegetated channels. Their assessments and applicability ranges are derived by comparing their predictions with measured values from a large dataset for different types of submerged and emergent vegetation gathered from the literature. The work includes assessing the performance of the sediment transport capacity formulae of Engelund and Hansen and van Rijn in the case of vegetated beds, as well as the value of the drag coefficient to be used for different types of vegetation and hydraulic conditions. The results provide a unique comparative overview of existing models for the assessment of the effects of vegetation on morphodynamics, highlighting their performances and applicability ranges. Copyright © 2014 John Wiley & Sons, Ltd.     

Experimental and computational study of trace element distribution between orthopyroxene and anhydrous silicate melt: substitution mechanisms and the effect of iron

Although orthopyroxene (Opx) is present during a wide range of magmatic differentiation processes in the terrestrial and lunar mantle, its effect on melt trace element contents is not well quantified. We present results of a combined experimental and computational study of trace element partitioning between Opx and anhydrous silicate melts. Experiments were performed in air at atmospheric pressure and temperatures ranging from 1,326 to 1,420°C in the system CaO–MgO–Al₂O₃–SiO₂ and subsystem CaO–MgO–SiO₂. We provide experimental partition coefficients for a wide range of trace elements (large ion lithophile: Li, Be, B, K, Rb, Sr, Cs, Ba, Th, U; rare earth elements, REE: La, Ce, Nd, Sm, Y, Yb, Lu; high field strength: Zr, Nb, Hf, Ta, Ti; transition metals: Sc, V, Cr, Co) for use in petrogenetic modelling. REE partition coefficients increase from $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 Although orthopyroxene (Opx) is present during a wide range of magmatic differentiation processes in the terrestrial and lunar mantle, its effect on melt trace element contents is not well quantified. We present results of a combined experimental and computational study of trace element partitioning between Opx and anhydrous silicate melts. Experiments were performed in air at atmospheric pressure and temperatures ranging from 1,326 to 1,420°C in the system CaO–MgO–Al₂O₃–SiO₂ and subsystem CaO–MgO–SiO₂. We provide experimental partition coefficients for a wide range of trace elements (large ion lithophile: Li, Be, B, K, Rb, Sr, Cs, Ba, Th, U; rare earth elements, REE: La, Ce, Nd, Sm, Y, Yb, Lu; high field strength: Zr, Nb, Hf, Ta, Ti; transition metals: Sc, V, Cr, Co) for use in petrogenetic modelling. REE partition coefficients increase from $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 to $ D_{\text{Lu}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.109 $ D_{\text{Lu}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.109 , D values for highly charged elements vary from $ D_{\text{Th}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0026 $ D_{\text{Th}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0026 through $ D_{\text{Nb}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0033 $ D_{\text{Nb}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0033 and $ D_{\text{U}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0066 $ D_{\text{U}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0066 to $ D_{\text{Ti}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.058 $ D_{\text{Ti}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.058 , and are all virtually independent of temperature. Cr and Co are the only compatible trace elements at the studied conditions. To elucidate charge-balancing mechanisms for incorporation of REE into Opx and to assess the possible influence of Fe on Opx-melt partitioning, we compare our experimental results with computer simulations. In these simulations, we examine major and minor trace element incorporation into the end-members enstatite (Mg₂Si₂O₆) and ferrosilite (Fe₂Si₂O₆). Calculated solution energies show that R²⁺ cations are more soluble in Opx than R³⁺ cations of similar size, consistent with experimental partitioning data. In addition, simulations show charge balancing of R³⁺ cations by coupled substitution with Li⁺ on the M1 site that is energetically favoured over coupled substitution involving Al–Si exchange on the tetrahedrally coordinated site. We derived best-fit values for ideal ionic radii r ₀, maximum partition coefficients D ₀, and apparent Young’s moduli E for substitutions onto the Opx M1 and M2 sites. Experimental r ₀ values for R³⁺ substitutions are 0.66–0.67 ? for M1 and 0.82–0.87 ? for M2. Simulations for enstatite result in r ₀ = 0.71–0.73 ? for M1 and ~0.79–0.87 ? for M2. Ferrosilite r ₀ values are systematically larger by ~0.05 ? for both M1 and M2. The latter is opposite to experimental literature data, which appear to show a slight decrease in $ r_{0}^{{{\text{M}}2}} $ r_{0}^{{{\text{M}}2}} in the presence of Fe. Additional systematic studies in Fe-bearing systems are required to resolve this inconsistency and to develop predictive Opx-melt partitioning models for use in terrestrial and lunar magmatic differentiation models.     

MERIS and the red-edge position

The Medium Resolution Imaging Spectrometer (MERIS) is a payload component of Envisat-1. MERIS will be operated over land with a standard 15 band setting acquiring images with a 300 m spatial resolution. The red-edge position (REP) is a promising variable for deriving foliar chlorophyll concentration, which plays an important role in ecosystem processes. The objectives of this paper are: (1) to study which factors effect the REP of vegetation, (2) to study whether this REP can be derived from the MERIS standard band setting and (3) to show what REP represents at the scale of MERIS data. Two different data sets were explored for simulating the REP using MERIS bands: (1) simulated data using reflectance models and (2) airborne reflectance spectra of an agricultural area obtained by the airborne visible-infrared imaging spectrometer (AVIRIS). A “linear method”, assuming a straight slope of the reflectance spectrum around the midpoint of the slope, was a robust method for determining the REP and the MERIS bands at 665, 708.75, 753.75 and 778.75 nm could be used for applying the “linear method” for REP estimation. Results of the translation to the scale of MERIS data were very promising for applying MERIS at, for instance, the ecosystem level.     

Trace element partitioning between ilmenite, armalcolite and anhydrous silicate melt: Implications for the formation of lunar high-Ti mare basalts

We performed a series of experiments at high pressures and temperatures to determine the partitioning of a wide range of trace elements between ilmenite (Ilm), armalcolite (Arm) and anhydrous lunar silicate melt, to constrain geochemical models of the formation of titanium-rich melts in the Moon. Experiments were performed in graphite-lined platinum capsules at pressures and temperatures ranging from 1.1 to 2.3 GPa and 1300-1400 °C using a synthetic Ti-enriched Apollo ‘black glass’ composition in the CaO-FeO-MgO-Al₂O₃-TiO₂-SiO₂ system. Ilmenite-melt and armalcolite-melt partition coefficients (D) show highly incompatible values for the rare earth elements (REE) with the light REE more incompatible compared to the heavy REE ( 0.0020 ± 0.0010 to 0.069 ± 0.010 for ilmenite; 0.0048 ± 0.0023 to 0.041 ± 0.008 for armalcolite). D values for the high field strength elements vary from highly incompatible for Th, U and to a lesser extent W (for ilmenite: 0.0013 ± 0.0008, 0.0035 ± 0.0015 and 0.039 ± 0.005, and for armalcolite 0.008 ± 0.003, 0.0048 ± 0.0022 and 0.062 ± 0.03), to mildly incompatible for Nb, Ta, Zr, and Hf (e.g. 0.28 ± 0.05 and : 0.76 ± 0.07). Both minerals fractionate the high field strength elements with D_Ta/D_Nb and D_Hf/D_Zr between 1.3 and 1.6 for ilmenite and 1.3 and 1.4 for armalcolite. Armalcolite is slightly more efficient at fractionating Hf from W during lunar magma ocean crystallisation, with D_Hf/D_W = 12-13 compared to 6.7-7.5 for ilmenite. The transition metals vary from mildly incompatible to compatible, with the highest compatibilities for Cr in ilmenite (D ∼ 7.5) and V in armalcolite (D ∼ 8.1). D values show no clear variation with pressure in the small range covered.Crystal lattice strain modelling of D values for di-, tri- and tetravalent trace elements shows that in ilmenite, divalent elements prefer to substitute for Fe while armalcolite data suggest REE replacing Mg. Tetravalent cations appear to preferentially substitute for Ti in both minerals, with the exception of Th and U that likely substitute for the larger Fe or Mg cations. Crystal lattice strain modelling is also used to identify and correct for very small (∼0.3 wt.%) melt contamination of trace element concentration determinations in crystals.Our results are used to model the Lu-Hf-Ti concentrations of lunar high-Ti mare basalts. The combination of their subchondritic Lu/Hf ratios and high TiO₂ contents requires preferential dissolution of ilmenite or armalcolite from late-stage, lunar magma ocean cumulates into low-Ti partial melts of deeper pyroxene-rich cumulates.