Melting of a Hydrous Mantle: II. Geochemistry of Crystals and Liquids Formed by Anatexis of Mantle Peridotite at High Pressures and High Temperatures as a Function of Controlled Activities of Water, Hydrogen, and Carbon Dioxide

The system peridotite-H₂O–CO₂ serves as a simplified modelfor the phase relations of mantle peridotite involving morethan one volatile component. Run products obtained in a studyof phase relations of four mantle peridotites in the presenceof H₂O- and (H₂O+CO₂)- bearing vapors and with controlled hydrogenfugacity (f_H2) at high pressures and temperatures have beensubjected to a detailed chemical investigation, principallyby the electron microprobe. Mg/(Mg+Fe) of all phases generally increases with increasingtemperature and with increasing Mg/(Mg+Fe) of the starting material.This ratio appears to decrease with increasing pressure forolivine, and for amphibole coexisting with garnet. Decreasingf_H2 from that of IW buffer to that of MH buffer decreases Mg/(Mg+Fe)of the partial melt from approximately 0-85 to approximately0.50, whereas the Fo content of coexisting olivine increasesslightly less than 3 per cent and the Mg/(Mg+Fe) of clinopyroxeneincreases about 4 per cent. However, the variations in Fo contentof olivines are within those observed in olivines from naturalmantle peridotite. The chemistry of other silicate mineralsdoes not significantly reflect variations of f_H2. Consequently,the peridotite mineralogy and/or chemistry is not a good indicatorfor the f_H2 conditions during crystallization. All crystalline phases, except amphibole, and to some extentgarnet, show increasing Cr content with increasing temperatureand increasing Cr content of the starting material, resultingin a positive correlation with Mg/(Mg+Fe). Partial melts aredepleted in Cr₂O₃ relative to the crystalline phases. High Mg/Mg+Fe)and Cr₂O₃ are thus expected in crystal residues after partialmelting. The absolute values depend on degree of melting andthe composition of the parent peridotite. Liquids formed by anatexis of mantle peridotite are andesiticunder conditions of X_H2O^v > 0.6 to at least 25 kb total pressureand to more than 200?C above the peridotite solidus. This observationsupports numerous suggestions that andesite genesis in islandarcs may result from partial melting of underlying peridotitemantle. In contrast to basaltic rocks, the absence of amphibole(paragasitic hornblende) does not affect the silica-saturatednature of the liquids. Increasing K₂O content of the startingmaterial (up to 1 wt. per cent K₂O) results in increasing potassiumcontent of the amphibole (1 wt. per cent K₂O) as well as theappearance of phlogopite. The liquid under these conditionsis relatively K₂0-poor (less than 1 wt. per cent K₂O). Partial melts are olivine normative with X_H2O 0.5, and initialliquids contain normative ol and ne at X_H2O 0.4. The alkalinityof these liquids increases with decreasing X_H2O below valuesof 0.5. The (ol+opx)-normative liquids resemble oceanic basaltswhereas (ol+ne)-normative liquids resemble olivine nepheliniteand melilite basalt. Low aHlo and high a_Co2 conditions may bethose under which kimberlites and related rocks are formed inthe mantle.     

Discovery of A New Class of Pulsating Stars: Gravity-mode Pulsators among Subdwarf B Stars

During the course of an ongoing CCD monitoring program to investigate low-level light variations in subdwarf B (sdB) stars, weserendipitously discovered a new class of multimode pulsators withperiods of the order of an hour. These periods are a factor of tenlonger than those of previously known multimode sdB pulsators (EC14026 stars), implying the new pulsations are due to gravity modes rather than pressure modes. The iron opacity instability that drives the short period EC 14026 stars is effective in hot sdB's. Thelong period pulsators are found only among cooler sdB stars, wherethey are surprisingly common. The mechanism responsible for excitingthe deeper g-modes in cool sdB's is currently unknown, but thetemperature and gravity range in which these stars occur must be animportant clue. We present the first observational results for thisnew class of pulsating sdB stars, and discuss some possible implications.     

Regressions and transgressions of the Baltic basin reflected by a new high-resolution deglacial and postglacial lithostratigraphy for Arkona Basin sediments (western Baltic Sea)

Seismoacoustic profiles from the Arkona Basin show a late Pleistocene and Holocene succession of several distinct reflectors. The physical, sedimentological, mineralogical and geochemical properties of more than 30 sediment cores were analysed in order to assign these reflectors to specific sedimentary discontinuity layers. Additionally, AMS ¹⁴C data and biostratigraphic information were gathered. Based on this multi‐proxy approach, seven lithostratigraphic units (AI, AII, B to F) were distinguished. These consist of fine‐grained clay, silt and mud, and are separated from each other by thin basin‐wide traceable sandy layers (S_ab‐S_ef). The most sensitive parameter to mark the lithostratigraphic boundaries is the weight percentage of the grain‐size fraction >63μm. In addition, some of the quartz‐grain‐dominated sandy layers cause the strong reflection lines recorded in seismoacoustic profiles. The sandy layers are interpreted to reflect enhanced hydrodynamic energy induced by episodes of basin‐wide water‐level low‐stand conditions. These low stands resulted from water‐level drops that occurred frequently during the Baltic Sea's history and presumably affected the entire Baltic basin. The thick fine‐grained units AI, AII to F, in which coarser material is absent, represent water‐level high‐stands. We conclude that the units AI and AII are Baltic Ice Lake sediments deposited before and after the Billingen‐1 regression, respectively. We assign the most prominent sandy layer S_ab to the final drainage of the Baltic Ice Lake (Billingen‐2), whereas the sandy layers between units B, C., D and E are related to the Yoldia Sea and Ancylus Lake regressions of the Baltic Sea's history. The uppermost fine‐grained unit F with its high organic carbon content contains marine sediments deposited after the Littorina Transgression. The macroscopically well‐visible sediment colour change from reddish/brown‐to‐grey, previously interpreted as a regional stratigraphic boundary, varies from core to core. It has been shown by our new data that this colour change has a diagenetic origin, and thus does not represent a stratigraphic boundary. Previous subdivisions therefore have to be revised.     

Late Quaternary environmental history of central Prins Karls Forland, western Svalbard

This paper presents the results from stratigraphic and geomorphologic investigations in the Poolepynten area, Prins Karls Forland, western Svalbard. Field mapping, soil profile development and ¹⁴C dating reveal the existence of at least two generations of raised beach deposits. Well-developed raised beaches rise to the Late Weichselian marine limit at 36 m a.s.l. Discontinuous pre-Late Weichselian beach deposits rise from the Late Weichselian marine limit to approximately 65 m a.s.l. Expansion of local glaciers in the area during the Late Weichselian is indicated by a till that locally overlies pre-Late Weichselian raised beach deposits. Stratigraphic data from coastal sections reveal two shallow marine units deposited during part of oxygen isotope stage 5. The two shallow marine units are separated by a subglacially deposited till that indicates an ice advance from Prins Karls Forland into the Forlandsundet basin some time during the latter part of stage 5. Discontinuous glaciofluvial deposits and a cobble-boulder lag could relate to a Late Weichselian local glacial advance across the coastal site. Late Weichselian/early Holocene beach deposits cap the sedimentary succession. Palaeotemperature estimates derived from amino acid ratios in subfossil marine molluscs indicate that the area has not been submerged or covered by warm based glacier ice for significant periods of time during the time interval ca. 70 ka to 10 ka.     

A Numerical Framework for Modeling Folds in Structural Geology

A numerical framework for modeling folds in structural geology is presented. This framework is based on a novel and recently published Hamilton–Jacobi formulation by which a continuum of layer boundaries of a fold is modeled as a propagating front. All the fold classes from the classical literature (parallel folds, similar folds, and other fold types with convergent and divergent dip isogons) are modeled in two and three dimensions as continua defined on a finite difference grid. The propagating front describing the fold geometry is governed by a static Hamilton–Jacobi equation, which is discretized by upwind finite differences and a dynamic stencil construction. This forms the basis of numerical solution by finite difference solvers such as fast marching and fast sweeping methods. A new robust and accurate scheme for initialization of finite difference solvers for the static Hamilton–Jacobi equation is also derived. The framework has been integrated in simulation software, and a numerical example is presented based on seismic data collected from the Karama Block in the North Makassar Strait outside Sulawesi.     

Geological modeling and simulation of CO2 injection in the Johansen formation

The Johansen formation is a candidate site for large-scale CO₂ storage offshore of the south-western coast of Norway. An overview of the geology for the Johansen formation and neighboring geological formations is given, together with a discussion of issues for geological and geophysical modelling and integrated fluid flow modelling. We further describe corresponding simulation models. Major issues to consider are capacity estimation and processes that could potentially cause CO₂ to leak out of the Johansen formation and into the formations above. Currently, these issues can only be investigated through numerical simulation. We consider the effect of different boundary conditions, sensitivity with respect to vertical grid refinement and permeability/transmisibility data, and the effect of residual gas saturations, since these strongly affect the CO₂-plume distribution. The geological study of the Johansen formation is performed based on available seismic and well data. Fluid simulations are performed using a commercial simulator capable of modelling CO₂ flow and transport by simple manipulation of input files and data. We provide details for the data and the model, with a particular focus on geology and geometry for the Johansen formation. The data set is made available for download online.     

NON-LINEAR AVO INVERSION FOR A STACK OF ANELASTIC LAYERS1

Parameters in a stack of homogeneous anelastic layers are estimated from seismic data, using the amplitude versus offset (AVO) variations and the travel-times. The unknown parameters in each layer are the layer thickness, the P-wave velocity, the S-wave velocity, the density and the quality factor. Dynamic ray tracing is used to solve the forward problem. Multiple reflections are included, but wave-mode conversions are not considered. The S-wave velocities are estimated from the PP reflection and transmission coefficients. The inverse problem is solved using a stabilized least-squares procedure. The Gauss-Newton approximation to the Hessian matrix is used, and the derivatives of the dynamic ray-tracing equation are calculated analytically for each iteration. A conventional velocity analysis, the common mid-point (CMP) stack and a set of CMP gathers are used to identify the number of layers and to establish initial estimates for the P-wave velocities and the layer thicknesses. The inversion is carried out globally for all parameters simultaneously or by a stepwise approach where a smaller number of parameters is considered in each step. We discuss several practical problems related to inversion of real data. The performance of the algorithm is tested on one synthetic and two real data sets. For the real data inversion, we explained up to 90% of the energy in the data. However, the reliability of the parameter estimates must at this stage be considered as uncertain.     

Convective Redistribution of Ozone and Oxides of Nitrogen in the Troposphere over Europe in Summer and Fall

The fluxes of ozone and NO_x out of the atmospheric boundary layer (ABL) over Europe are calculated in a mesoscale chemical transport model (MCT) and compared with the net chemical production or destruction of ozone and the emissions of precursors within the ABL for two 10 days' periods which had quite different synoptic situations and levels of photochemical activity (1–10 July 1991 (JUL91) and 26 October–4 November 1994 (ON94)). Over the European continent, about 8% of the NO_x emissions were brought from the ABL to the free troposphere as NO_x, while about 15% of the NO_x emissions were brought to the free troposphere as NO_y–NO_x, i.e. as PAN or HNO₃. The convection dominates over the synoptic scale vertical advection as a transport mechanism both for NO_x and NO_y out of the boundary layer in the summertime high pressure situation (JUL91), while in the fall situation (ON94) the convective part was calculated to be the smallest. NO_x was almost completely transformed to NO_y–NO_x or removed within the ABL. Also for NO_y the major part of the atmospheric cycle is confined to the ABL both for JUL91 and ON94. The vertical transport time out of the ABL is of the order of 100h both for the total model domain and over the European continent. The net convective exchange of ozone from the ABL is not a dominant process for the amount of ozone in the ABL averaged over 10 days and the whole domain, but convection reduces the maximum ozone concentration in episodes significantly. The ozone producing efficiency of NO_x is calculated to increase with height to typically 15–20 in the upper half of the troposphere from around 5 in the ABL, but in the middle free troposphere the concentration of NO_x is often too low to cause net chemical formation of ozone there.     

INFORMATION NOT DATA: FUTURE DEVELOPMENT OF ROAD WEATHER INFORMATION SYSTEMS

The Swedish Road Weather Information System is a tool for road winter maintenance giving warnings about risk of slipperiness and severe weather conditions. The system is based on climatological research at the University of Göteborg and was initiated by Professor Sven Lindqvist during the 1970s. Knowledge about micro‐and local climate variations is crucial when it comes to forecasting and giving warnings about different kinds of weather‐related traffic disturbances. The progress of climatological research has made it possible to increase the use of road weather data into different applications. Future development will be focused upon development of the road weather information system towards a user‐friendly maintenance decision support system which will enhance the benefits of the road weather information system.