New Olivine Reference Material for In Situ Microanalysis

A new olivine reference material – MongOL Sh11‐2 – for in situ analysis has been prepared from the central portion of a large (20 × 20 × 10 cm) mantle peridotite xenolith from a ~ 0.5 My old basaltic breccia at Shavaryn‐Tsaram, Tariat region, central Mongolia. The xenolith is a fertile mantle lherzolite with minimal signs of alteration. Approximately 10 g of 0.5–2 mm gem quality olivine fragments were separated under binocular microscope and analysed by EPMA, LA‐ICP‐MS, SIMS and bulk analytical methods (ID‐ICP‐MS for Mg and Fe, XRF, ICP‐MS) for major, minor and trace elements at six institutions world‐wide. The results show that the olivine fragments are sufficiently homogeneous with respect to major (Mg, Fe, Si), minor and trace elements. Significant inhomogeneity was revealed only for phosphorus (homogeneity index of 12.4), whereas Li, Na, Al, Sc, Ti and Cr show minor inhomogeneity (homogeneity index of 1–2). The presence of some mineral and fluid‐melt micro‐inclusions may be responsible for the inconsistency in mass fractions obtained by in situ and bulk analytical methods for Al, Cu, Sr, Zr, Ga, Dy and Ho. Here we report reference and information values for twenty‐seven major, minor and trace elements.     

Non-linear Bayesian joint inversion of seismic reflection coefficients

Inversion of seismic reflection coefficients is formulated in a Bayesian framework. Measured reflection coefficients and model parameters are assigned statistical distributions based on information known prior to the inversion, and together with the forward model uncertainties are propagated into the final result. This enables a quantification of the reliability of the inversion. Quadratic approximations to the Zoeppritz equations are used as the forward model. Compared with the linear approximations the bias is reduced and the uncertainty estimate is more reliable. The differences when using the quadratic approximations and the exact expressions are minor. The solution algorithm is sampling based, and because of the non-linear forward model, the Metropolis–Hastings algorithm is used. To achieve convergence it is important to keep strict control of the acceptance probability in the algorithm. Joint inversion using information from both reflected PP waves and converted PS waves yields smaller bias and reduced uncertainty compared to using only reflected PP waves.     

Experimental and numerical modelling of sedimentation in a rectangular shallow basin

Numerical simulation of flows in shallow reservoirs has to be checked for its consistency in predicting real flow conditions and sedimentation patterns. Typical flow patterns may exhibit flow separation at the inlet, accompanied by several recirculation and stagnation areas all over the reservoir surface. The aim of the present research project is to study the influence of the geometry of a reservoir on sediment transport and deposition numerically and experimentally, focusing on a prototype reservoir depth between 5 and 15 m as well as suspended sediment transport.
A series of numerical simulations is presented and compared with scaled laboratory experiments, with the objective of testing the sensitivity to different flow and sediment parameters and different turbulence closure schemes. Different scenarios are analyzed and a detailed comparison of preliminary laboratory tests and some selected simulations are presented.
The laboratory experiments show that suspended sediment transport and deposition are determined by the initial flow pattern and by the upstream and downstream boundary conditions. In the experiments, deposition in the rectangular basin systematically developed along the left bank, although inflow and outflow were positioned symmetrically along the centre of the basin. Three major horizontal eddies developed influencing the sediment deposition pattern. Although asymmetric flow patterns are privileged, a symmetric pattern can appear from time to time. This particular behaviour could also be reproduced by a two-dimensional depth-averaged flow and sediment transport model （CCHE2D）. The paper presents numerical simulations using different turbulence closure schemes （k-ε and eddy viscosity models）. In spite of the symmetric setup, these generally produced an asymmetric flow pattern that can easily switch sides depending on the assumptions made for the initial and boundary conditions. When using the laboratory experiment as a reference, the most reliable numerical results have been obtai     

Volcanic evolution of Volcán Aucanquilcha: a long-lived dacite volcano in the Central Andes of northern Chile

Volcán Aucanquilcha, northern Chile, has produced ∼37 km³ of dacite (63–66 wt% silica), mainly as lavas with ubiquitous magmatic inclusions (59–62 wt% silica) over the last ∼1 million years. A pyroclastic flow deposit related to dome collapse occurs on the western side of the edifice and a debris avalanche deposit occurs on the eastern side. The >6,000-m high edifice defines a 9-km E–W ridge and lies at the center of a cluster of more than 15 volcanoes, the Aucanquilcha Volcanic Cluster, that has been active for at least the past 11 million years. The E–W alignment of vents is nearly orthogonal to the arc axis. A majority of Volcán Aucanquilcha was constructed during the first 200,000 years of eruption, whereas the last 800,000 years have added little additional volume. The peak eruptive rate during the edifice-building phases was ∼0.16 km³/ka and the later eruptive rate was ∼0.02 km³/ka. Comparable dacite volcanoes elsewhere show a similar pattern of high volcanic productivity during the early stages and punctuated rather than continuous activity. Volcán Aucanquilcha lavas are dominated by phenocrysts of plagioclase, accompanied by two populations of amphibole, biotite, clinopyroxene, Fe–Ti oxides and (or) orthopyroxene. Accessory phases include zircon, apatite and rare quartz and sanidine. One amphibole population is pargasite and the other is hornblende. The homogeneity of dacite lava from Volcán Aucanquilcha contrasts with the heterogeneity (52–66 wt% silica) at nearby Volcán Ollagüe, which has been active over roughly the same period of time. We attribute this homogeneity at Aucanquilcha to the thermal development of the crust underneath the volcano resulting from protracted magmatism there, whereas Volcán Ollagüe lacks this magmatic legacy.     

Multiple volcano deformation sources in a post-rifting period: 1989–2005 behaviour of Krafla,Iceland constrained by levelling,tilt and GPS observations

The Krafla rifting episode, which occurred in North Iceland in 1975–1984, was followed by inflation of a shallow magma chamber until 1989. At that time, gradual subsidence began above the magma chamber and has continued to the present at a declining rate. Pressure decrease in a shallow magma chamber is not the only source of deformation at Krafla, as other deformation processes are driven by exploitation of two geothermal fields, together with plate spreading. In addition, deep-seated magma accumulation appears to take place, with its centre ∼ 10 km north of the Krafla caldera. The relative strength of these sources has varied with time. New results from a levelling survey and GPS measurements in 2005 allow an updated view on the deformation field. Deformation rates spanning 2000–2005 are the lowest recorded in the 30-year history of geodetic studies at the volcano. The inferred rate of 2000–2005 subsidence related to processes in the shallow magma chamber is less than 0.3 cm/yr whereas it was ∼ 5 cm/yr in 1989–1992. Currently, the highest rate of subsidence takes place in the Leirbotnar area, within the Krafla caldera, and appears to be a result of geothermal exploitation.     

Proglacial sediment trapping in recently formed Silt Lake,upper Lillooet Valley,Coast Mountains,British Columbia

The sedimentology of proglacial Silt Lake was assessed by lake sediment coring and monitoring of lacustrine processes during a late‐summer period of high glacier melt to characterize sediment delivery from the heavily glacierized catchment and investigate the sediment trapping dynamics of this upland lake. A complete varve chronology was established for a distal basin of the lake which was exposed by Lillooet Glacier retreat between 1947 and 1962. The varve record showed decreasing sedimentation rates in the basin while the glacier retreated, and as the lake became free of ice contact in the early 1970s. Although recession has continued over recent decades, and glacier proximity to the lake has, therefore, continued decreasing, lacustrine sedimentation rates are now accelerating due to changing basin morphometry caused by delta progradation. Over shorter time scales, lake sedimentation patterns respond to changing runoff conditions, including late‐summer glacier melt intensity, intra‐annual flooding events, diumal runoff fluctuations, and within‐lake turbidity currents. Turbidity currents included quasi‐regular flows during high diurnal discharges and an episodic flushing of temporarily stored sediment from the sandur or delta at a time of low stage. Suspended sediment yield to Silt Lake is estimated to exceed 10³ Mg km^?2 a^?1, a magnitude that surpasses previous local and regional yield estimates for the glacierized headwaters of the Lillooet River valley. Since Silt Lake currently traps a significant prooportion of that upland sediment supply, and the trapping efficiency of the basin has been variable at decadal time scales, the formation and continued development of Lilt Lake has likely had a significant influence on downstream sediment delivery. Lacustrine sediment‐based proxies of long‐term hydroclimatic variability being developed in glacially distal settings should include provisions for dynamic sediment trapping effects in upstream water bodies that often form in the active proglacial environment. Copyright © 2008 John Wiley & Sons, Ltd.     

Mixed Integer-Real Valued Adjustment (IRA) Problems: GPS Initial Cycle Ambiguity Resolution by Means of the LLL Algorithm

In order to achieve to GPS solutions of first-order accuracy and integrity, carrier phase observations as well as pseudorange observations have to be adjusted with respect to a linear/linearized model. Here the problem of mixed integer-real valued parameter adjustment (IRA) is met. Indeed, integer cycle ambiguity unknowns have to be estimated and tested. At first we review the three concepts to deal with IRA: (i) DDD or triple difference observations are produced by a properly chosen difference operator and choice of basis, namely being free of integer-valued unknowns (ii) The real-valued unknown parameters are eliminated by a Gauss elimination step while the remaining integer-valued unknown parameters (initial cycle ambiguities) are determined by Quadratic Programming and (iii) a RA substitute model is firstly implemented (real-valued estimates of initial cycle ambiguities) and secondly a minimum distance map is designed which operates on the real-valued approximation of integers with respect to the integer data in a lattice. This is the place where the integer Gram-Schmidt orthogonalization by means of the LLL algorithm (modified LLL algorithm) is applied being illustrated by four examples. In particular, we prove that in general it is impossible to transform an oblique base of a lattice to an orthogonal base by Gram-Schmidt orthogonalization where its matrix enties are integer. The volume preserving Gram-Schmidt orthogonalization operator constraint to integer entries produces “almost orthogonal” bases which, in turn, can be used to produce the integer-valued unknown parameters (initial cycle ambiguities) from the LLL algorithm (modified LLL algorithm). Systematic errors generated by “almost orthogonal” lattice bases are quantified by A. K. Lenstra et al. (1982) as well as M. Pohst (1987). The solution point of Integer Least Squares generated by the LLL algorithm is = (L')⁻¹[L'◯] ∈ ℤ ^m where L is the lower triangular Gram-Schmidt matrix rounded to nearest integers, [L], and = [L'◯] are the nearest integers of L'◯, ◯ being the real valued approximation of z ∈ ℤ ^m , the m-dimensional lattice space Λ. Indeed due to “almost orthogonality” of the integer Gram-Schmidt procedure, the solution point is only suboptimal, only close to “least squares.” ? 2000 John Wiley & Sons, Inc.     

Stability and Grid Dispersion of the P-SV 4th-Order Staggered-Grid Finite-Difference Schemes

Stability and grid dispersion in the P-SV 4th-order in space, 2nd-order in time, displacement-stress staggered-grid finite-difference scheme is investigated in the case of a homogeneous unbounded medium. All results, however, also apply to the velocity-stress and displacement- velocity-stress finite-difference schemes. Independent stability conditions for the P and S waves are obtained by exact separation of equations for the two types of waves. Since the S-wave group velocity can differ from the actual velocity as much as 5% for the sampling ratio 1/5, commonly used in numerical modelling, the sampling of the minimum S wavelength by 6 grid spacings (with the velocity difference not larger than 2.5%) is recommended. Grid dispersion is strongest for a wave propagating in a direction of a coordinate axis and weakest for a wave propagating along a plane diagonal. Grid dispersion in the 4 ^th -order scheme for the sampling ratios s = 1/5 and s = 1/6 is smaller than grid dispersion in the 2 ^nd -order scheme for s = 1/10 and s = 1/12, respectively.