On "the complete basis set limit" and plane-wave methods in first-principles simulations of water

Water structure, measured by the height of the first peak in oxygen-oxygen radial distributions, is converged with respect to plane-wave basis energy cutoffs for ab initio molecular dynamics simulations, confirming the reliability of plane-wave methods.Comment: 9 pages, 3 figure

Dynamics of Magma Mixing in Partially Crystallized Magma Chambers: Textural and Petrological Constraints from the Basal Complex of the Austurhorn Intrusion (SE Iceland)

The Tertiary Austurhorn intrusive complex in SE Iceland represents an exhumed magma chamber that has recorded an extensive history of magma mixing and mingling. The basal part of the intrusion consists predominantly of granophyres that have been intensively and repeatedly intruded by more mafic magma. This association of granophyres, basic and hybrid rocks at Austurhorn is referred to in the literature as a ‘net-veined' complex, but field relations suggest a much more complex emplacement history. Here we present petrological and physical constraints on the various processes that resulted in magma mixing and mingling and the formation of different generations of hybrid rocks at Austurhorn. The complexity of the mixing and mingling processes increases towards the inferred centre of the intrusion, where chaotic hybrid rocks dominate the exposed lithology. Complex cross-cutting relations between different hybrid generations strongly suggest multiple magma injection and reheating events in the basal part of the shallow magma chamber. Model calculations employing distribution coefficients based on rare earth element concentrations reveal that early stage hybrid magma generations formed by pure endmember mixing between felsic and mafic magma with about 10% mafic fraction in the hybrids. With repeated injections of mafic magma into the base of the magma chamber, the intruding magma interacted to a greater extent with pre-existing hybrids. This led to the formation of hybrid magma compositions that are shifted towards the mafic endmember over time, with up to 30% of the mafic fraction in the hybrids. These mixing processes are recorded in the zonation patterns of clinopyroxene and plagioclase phenocrysts; the latter have been divided into four main groups by cross-correlation analysis. Melt viscosity calculations were performed to constrain the possible conditions of magma mixing and the results indicate that the interaction of the contrasting magmas most probably occurred at temperatures of approximately 1000°C up to 1120°C. This suggests that the initiation of effective magma mixing requires local superheating of the felsic magmas, thereby confining the process to areas of localized, substantial mafic magma injectio

New, efficient, and accurate high order derivative and dissipation operators satisfying summation by parts, and applications in three-dimensional multi-block evolutions

We construct new, efficient, and accurate high-order finite differencing operators which satisfy summation by parts. Since these operators are not uniquely defined, we consider several optimization criteria: minimizing the bandwidth, the truncation error on the boundary points, the spectral radius, or a combination of these. We examine in detail a set of operators that are up to tenth order accurate in the interior, and we surprisingly find that a combination of these optimizations can improve the operators' spectral radius and accuracy by orders of magnitude in certain cases. We also construct high-order dissipation operators that are compatible with these new finite difference operators and which are semi-definite with respect to the appropriate summation by parts scalar product. We test the stability and accuracy of these new difference and dissipation operators by evolving a three-dimensional scalar wave equation on a spherical domain consisting of seven blocks, each discretized with a structured grid, and connected through penalty boundary conditions.Comment: 16 pages, 9 figures. The files with the coefficients for the derivative and dissipation operators can be accessed by downloading the source code for the document. The files are located in the "coeffs" subdirector

Hållbart växtskydd. Analys av olika strategier för att minska riskerna med kemiska växtskyddsmedel

Risker som är förknippade med kemiska växtskyddsmedel kan sammanfattas i fem områden: risker för användare, risker för rester i livsmedel samt i yt- och grundvatten och risker för natursystemet (t.ex. förlorad biodiversitet). De viktigaste fördelarna med växtskyddsmedel är att de bidrar till skördeökningar och att de förbättrar möjligheten till stabila skördenivåer, oavsett årsmån. Det finns olika strategier för att hantera växtskyddet i dagens jordbruk och beroende på vilken strategi som väljs förändras riskerna som är förknippade med användningen av växtskyddsmedel. Syftet med MAT 21-projektet Hållbart Växtskydd har varit att analysera den miljömässiga hållbarheten i olika växtskyddsstrategier med hjälp av konkreta exempel. Ambitionen är att ge vägledning om genomförandet av politiska miljömål, tolkning av miljöanalyser och utveckling av kvalitets- och miljöledningssystem i livsmedelsföretag. I projektet har en varierad växtföljd setts som en viktig förebyggande åtgärd för att minska riskerna med växtskyddsmedel. En diversifierad växtföljd bidrar till att reducera vissa typer av ogräs och till att minska trycket av många svampsjukdomar. Mekanisk ogräsreglering ger tillsammans med förebyggande åtgärder möjligheter att minska användningen av växtskyddsmedel. Exempel på tekniska lösningar som kan minska riskerna är bandsprutning i kombination med radrensning eller ny bredsprutningsteknik med lufttillsats vilket reducerar vindavdriften vid en applikation med växtskyddsmedel. För att analysera riskerna som är förknippade med växtskyddsmedel användes modellen PRI-Farm som är en riskindikatormodell utvecklad av Kemikalieinspektionen. I modellen poängsätts faktorer för fara och exponering vid applikation av växtskyddsmedel. Farofaktorerna utrycker belastningen av de inneboende egenskaperna hos preparatens aktiva substans som persistens, bioackumulering och rörlighet. Exponeringsfaktorerna beaktar applikationsmetod, använd dos, markegenskaper och spruttidpunkt. Resultaten vägs samman i en poängsumma (ett miljöriskindex), där risk för transport till yt- och grundvatten väger tungt. Även ett hälsoriskindex kan beräknas, vilket inte gjorts i denna studie. Modellen kan inte användas för att kvantifiera risker utan är i stället ett verktyg för att kunna göra jämförande bedömningar mellan olika strategier för växtskyddsinsatser, identifiera problemanvändningar eller för att visa på risktrender över tiden. För att analysera andra miljöeffekter (t.ex. klimatpåverkan och energianvändning) användes metodiken för livscykelanalys (LCA). Detta är en metod som används för att göra miljöbedömningar av produkter och tjänster. Man genomför beräkningar för produktens potentiella bidrag till olika miljöeffekter, såsom växthuseffekt och försurning, och man summerar även den förbrukning av viktiga resurser som produkten givit upphov till. Tanken är att dessa beräkningar skall kunna ge information som kan leda till att man kan minska miljöpåverkan från en viss produkt på ett effektivt sätt. Resultaten kan också ge vägledning vid val mellan produkter som är likvärdiga i övrigt. I fallstudien av den allsidiga växtodlingsgården jämfördes också nuläget med ett scenario där gården lades om till ekologisk produktion. I detta alternativ användes inga kemiska växtskyddsmedel. Omläggningen till ekologisk produktion påverkade även andra miljöeffekter, förutom riskerna med växtskyddsmedel. En negativ effekt av det ekologiska scenariot var att grödan kepalök fick utgå p.g.a. att risken bedömdes som för stor för ett stort skördebortfall utan tillgång på fungicider. På grund av en lägre skördenivå i ekologisk produktion behövdes en betydligt större yta för att producera en given mängd livsmedel. I en ekonomisk analys jämfördes nuläget på den allsidiga växtodlingsgården med ett ekologiskt scenario. Resultatet var likvärdigt för de två alternativen. För att ekologisk odling skall vara en ekonomiskt lyckad satsning på gårdar i intensiva slättbygder med odling av specialgrödor är det, med dagens nivå på miljöstödet för ekologisk produktion, helt nödvändigt att det finns ett bra merpris för produkterna. I en fallstudie av en allsidig växtodlingsgård jämfördes nuläget med ett scenario där olika strategier testades som minskade riskerna med växtskyddsmedel utan att skördenivån ändrades nämnvärt. Detta scenario (kallat Bästa Möjliga Teknik) innebar användning av mekanisk och/eller termisk behandling som alternativ till konventionell bekämpning i vissa grödor samt medvetna val av växtskyddsmedel som hade mindre dåliga miljöegenskaper. Fallstudien visade på flera olika möjligheter till att reducera riskerna i konventionell drift. Eftersom gården i studien hade många radodlade specialgrödor bedömdes bandsprutning kombinerat med radrensning vara en intressant förbättring. Att väsentligt reducera användningen av fungicider i potatis och lök utan att äventyra skördenivån bedömdes som mycket svårt i konventionell drift, även med bästa odlingsteknik. I fallstudien av den allsidiga växtodlingsgården jämfördes också nuläget med ett scenario där gården lades om till ekologisk produktion. Detta alternativ minskar riskerna med växtskydds-medel radikalt men påverkar även andra miljöeffekter. För många grödor minskade energianvändning, liksom utsläpp av växthusgaser vid en övergång till ekologisk odling. En negativ effekt av det ekologiska scenariot var att grödan kepalök fick utgå p g a att risken bedömdes som för stor för ett betydande skördebortfall utan tillgång till fungicider. Ur miljö- och resurssynpunkt var den mest negativa effekten att det ekologiska scenariot krävde mycket större yta för att producera en given mängd livsmedel. I en ekonomisk analys jämfördes nuläget på den allsidiga växtodlingsgården med ett ekologiskt scenario. Resultatet var likvärdigt för de två alternativen. För att ekologisk odling skall vara en ekonomiskt lyckad satsning på gårdar i intensiva slättbygder är det helt nödvändigt att det finns ett bra merpris för produkterna. Studien visar att det finns flera olika strategier för att minska riskerna med användningen av växtskyddsmedel. En varierad växtföljd är en grundläggande förutsättning i alla odlingssystem och inom detta område finns det mycket gammal kunskap som behöver väckas till liv. Miljöanalyserna visar också att många tekniska lösningar, som leder till minskat behov av växtskyddsmedel, kan tillämpas utan att andra miljöeffekter påverkas nämnvärt. Att välja bort växtskyddsmedel som har mindre goda miljöeffekter är en strategi som förefaller ha stor potential för att minska riskerna. PRI-Farm kan bli ett viktigt beslutsstöd i lantbrukets arbete för bättre produktval och säkrare växtskydd. Ekologisk produktion har goda förutsättningar att bli en viktig strategi för minskade risker i de utökade vattenskyddsområden som kommunerna nyligen har beslutat om. För att denna strategi skall bli lyckosam krävs att marknaden fortsätter att öka för ekologiska produkter i Sverige, så att produkternas merpris tillsammans med samhällets stöd via miljöersättningsprogram fortsätter att betala merkostnaden för denna produktionsform

Would You Choose to be Happy? Tradeoffs Between Happiness and the Other Dimensions of Life in a Large Population Survey

A large literature documents the correlates and causes of subjective well-being, or happiness. But few studies have investigated whether people choose happiness. Is happiness all that people want from life, or are they willing to sacrifice it for other attributes, such as income and health? Tackling this question has largely been the preserve of philosophers. In this article, we find out just how much happiness matters to ordinary citizens. Our sample consists of nearly 13,000 members of the UK and US general populations. We ask them to choose between, and make judgments over, lives that are high (or low) in different types of happiness and low (or high) in income, physical health, family, career success, or education. We find that people by and large choose the life that is highest in happiness but health is by far the most important other concern, with considerable numbers of people choosing to be healthy rather than happy. We discuss some possible reasons for this preference

An Improved High Order Finite Difference Method for Non-conforming Grid Interfaces for the Wave Equation

This paper presents an extension of a recently developed high order finite difference method for the wave equation on a grid with non-conforming interfaces. The stability proof of the existing methods relies on the interpolation operators being norm-contracting, which is satisfied by the second and fourth order operators, but not by the sixth order operator. We construct new penalty terms to impose interface conditions such that the stability proof does not require the norm-contracting condition. As a consequence, the sixth order accurate scheme is also provably stable. Numerical experiments demonstrate the improved stability and accuracy property

Higher order finite difference schemes for the magnetic induction equations

We describe high order accurate and stable finite difference schemes for the initial-boundary value problem associated with the magnetic induction equations. These equations model the evolution of a magnetic field due to a given velocity field. The finite difference schemes are based on Summation by Parts (SBP) operators for spatial derivatives and a Simultaneous Approximation Term (SAT) technique for imposing boundary conditions. We present various numerical experiments that demonstrate both the stability as well as high order of accuracy of the schemes.Comment: 20 page

Element Partitioning between Immiscible Carbonatite and Silicate Melts for Dry and H2O-bearing Systems at 1-3 GPa

Carbonatite and silicate rocks occurring within a single magmatic complex may originate through liquid immiscibility. We thus experimentally determined carbonatite/silicate melt partition coefficients (Dcarbonate melt/silicate melt, hereafter D) for 45 elements to understand their systematics as a function of melt composition and to provide a tool for identifying the possible conjugate nature of silicate and carbonatite magmas. Static and, when necessary, centrifuging piston cylinder experiments were performed at 1-3 GPa, 1150-1260°C such that two well-separated melts resulted. Bulk compositions had Na K, Na ∼ K, and Na K; for the latter we also varied bulk H2O (0-4 wt %) and SiO2 contents. Oxygen fugacities were between iron-wüstite and slightly below hematite-magnetite and were not found to exert significant control on partitioning. Under dry conditions alkali and alkaline earth elements partition into the carbonatite melt, as did Mo and P (DMo >8, DP= 1·6-3·3). High field strength elements (HFSE) prefer the silicate melt, most strongly Hf (DHf = 0·04). The REE have partition coefficients around unity with DLa/Lu = 1·6-2·3. Transition metals have D < 1 except for Cu and V (DCu ∼ 1·3, DV = 0·95-2). The small variability of the partition coefficients in all dry experiments can be explained by a comparable width of the miscibility gap, which appears to be flat-topped in our dry bulk compositions. For all carbonatite and silicate melts, Nb/Ta and Zr/Hf fractionate by factors of 1·3-3·0, in most cases much more strongly than in silicate-oxide systems. With the exception of the alkalis, partition coefficients for the H2O-bearing systems are similar to those for the anhydrous ones, but are shifted in favour of the carbonatite melt by up to an order of magnitude. An increase of bulk silica and thus SiO2 in the silicate melt (from 35 to 69 wt %) has a similar effect. Two types of trace element partitioning with changing melt composition can be observed. The magnitude of the partition coefficients increases for the alkalis and alkaline earths with the width of the miscibility gap, whereas partition coefficients for the REE shift by almost two orders of magnitude from partitioning into the silicate melt (DLa = 0·47) to strongly partitioning into the carbonatite melt (DLa = 38), whereas DLa/DLu varies by only a factor of three. The partitioning behavior can be rationalized as a function of ionic potential (Z/r). Alkali and alkaline earth elements follow a trend, the slope of which depends on the K/Na ratio and H2O content. Contrasting the sodic and potassic systems, alkalis have a positive correlation in D vs Z/r space in the potassic case and Cs to K partition into the silicate melt in the presence of H2O. For the divalent third row transition metals on the one hand and for the tri- and tetravalent REE and HFSE on the other, two trends of negative correlation of D vs Z/r can be defined. Nevertheless, the highest ionic strength network-modifying cations (V, Nb, Ta, Ti and Mo) do not follow any trend; understanding their behavior would require knowledge of their bonding environment in the carbonatite melt. Strong partitioning of REE into the carbonatite melt (DREE = 5·8-38·0) occurs only in H2O-rich compositions for which carbonatites unmix from evolved alkaline melts with the conjugate silicate melt being siliceous. We thus speculate that upon hydrous carbonatite crystallization, the consequent saturation in fluids may lead to hydrothermal systems concentrating REE in secondary deposit