The composition and crystallinity of the near-surface regions of weathered alkali feldspars

Our ability to identify thin non-stoichiometric and amorphous layers beneath mineral surfaces has been tested by undertaking X-ray photoelectron spectroscopy (XPS) and transmission electron microscopy (TEM) work on alkali feldspars from pH 1 dissolution experiments. The outcomes of this work were used to help interpret XPS and TEM results from alkali feldspars weathered for <10,000 years in soils overlying the Shap Granite (north-west England). The chemistry of effluent solutions indicates that silica-rich layers a few nanometers in thickness formed during the pH 1 experiments. These layers can be successfully identified by XPS and have lower Al/Si, Na/Si, K/Si and Ca/Si values than the outermost ∼9 nm of unweathered controls. Development of Al-Si non-stoichiometry is coupled with loss of crystal structure to produce amorphous layers that are identifiable by TEM where >∼2.5 nm thick, whereas the crystallinity of albite is retained despite leaching of Na to depths of tens to hundreds on nanometers. Integration of XPS data over the outermost 6-9 nm of naturally weathered Shap feldspars shows that they have stoichiometric Al/Si and K/Si ratios, which is consistent with findings of previous TEM work on the same material that they lack amorphous layers. There is some XPS evidence for loss of K from the outermost couple of nanometers of Shap orthoclase, and the possibility of leaching of Na from albite to greater depths cannot be excluded using the XPS or TEM results. This study demonstrates that the leached layer model, as formulated from laboratory experiments, is inapplicable to the weathering of alkali feldspars within acidic soils, which is an essentially stoichiometric reaction.     

Multiparametric Traveltime Inversion

In conventional seismic processing, the classical algorithm of Hubral and Krey is routinely applied to extract an initial macrovelocity model that consists of a stack of homogeneous layers bounded by curved interfaces. Input for the algorithm are identified primary reflections together with normal moveout (NMO) velocities, as derived from a previous velocity analysis conducted on common midpoint (CMP) data. This work presents a modified version of the Hubral and Krey algorithm that is designed to extend the original version in two ways, namely (a) it makes an advantageous use of previously obtained common-reflection-surface (CRS) attributes as its input and (b) it also allows for gradient layer velocities in depth. A new strategy to recover interfaces as optimized cubic splines is also proposed. Some synthetic examples are provided to illustrate and explain the implementation of the method.     

Floodplain–river ecosystems: lateral connections and the implications of human interference

Floodplains are ecotones that form a transition between aquatic and terrestrial environments. These important ecosystems can be described as dynamic spatial mosaics in which water plays an important role in connecting various patches on the floodplain surface. Hydrological connections facilitate the exchange of carbon and nutrients between the river channel and the floodplain and therefore influence the productivity of the entire river system. This paper examines the influence of hydrological connections on the potential exchange of dissolved organic carbon between a large Australian floodplain to a river channel, and the effects of land and water developments on these exchanges. The paper proposes that an understanding of floodplain ecosystems requires an interdisciplinary approach—a recognition of the importance of the three disciplines hydrology, geomorphology and ecology. Large-scale water-resources and floodplain development has significantly altered the spatial and temporal patterns of hydrological characteristics in the Lower Balonne floodplain, Australia. The magnitude, frequency and duration of flooding events have all been reduced. The construction of levees and water storages has also reduced the reactive floodplain surface area. The presented data show the impacts of these changes on the potential supply of dissolved organic carbon from the floodplain surface during periods of inundation. Annual reductions of up to 1293 tonnes of dissolved organic carbon supply were noted and reductions were especially significant for floods with an average recurrence interval of 2 years or less. Some small flood events no longer facilitate the potential supply of dissolved organic carbon from the floodplain to the river channel because of water-resources and floodplain developments.     

Nucleation and crystallization of otavite, witherite, calcite, strontianite, hydrozincite, and hydrocerussite by CO2 membrane diffusion technique

Otavite, witherite, calcite, strontianite, hydrozincite and hydrocerussite nucleation was induced by CO₂ diffusion through a polyethylene membrane into a metal bearing solution. Nucleation and ongoing precipitation was followed at 25 ± 1 °C by the chemical evolution of the solution and the consumption of sodium hydroxide (pH-stat conditions). X-ray diffraction patterns as well as FT-infrared and Raman spectra confirmed the formation of well crystallized solids, except for less crystalline hydrozincite. In several experiments simonkolleite and laurionite precipitated concurrently with hydrozincite and hydrocerussite. The carbonate end-member minerals, smithsonite and cerussite, however were not formed. Carbonate minerals crystallized as spherical aggregates of thin layered otavite, spherical orientated witherite needles, rhombohedral calcite, pseudo-hexagonal strontianite laths, fibrous lumps of hydrozincite and planar hexagonal hydrocerussite crystals. Crystal and aggregate sizes range between 1 and 100 μm. Nucleation occurred at well defined reaction times and distinct critical supersaturation indices (SI_crit). The time for nucleation at constant pH decreased as the initial metal concentration increased for a given solid. The SI_crit values decreased in the order of hydrocerussite (3.2), otavite (2.5), strontianite (1.6), witherite (0.9) and calcite (0.7).     

Behaviour of radiogenic Pb in zircon during ultrahigh-temperature metamorphism: an ion imaging and ion tomography case study from the Kerala Khondalite Belt,southern India

Zircon crystals from a locally charnockitized Paleoproterozoic high-K metagranite from the Kerala Khondalite Belt (KKB) of southern India have been investigated by high-spatial resolution secondary ion mass spectrometry analysis of U–Th–Pb and rare earth elements (REE), together with scanning ion imaging and scanning ion tomography (depth-profiled ion imaging). The spot analyses constrain the magmatic crystallization age of the metagranite to ca. 1,850 Ma, with ultrahigh-temperature (UHT) metamorphism occurring at ca. 570 Ma and superimposed charnockite formation at ca. 520–510 Ma, while the ion imaging reveals a patchy distribution of radiogenic Pb throughout the zircon cores. Middle- to heavy-REE depletion in ca. 570 Ma zircon rims suggests that these grew in equilibrium with garnet and therefore date the UHT metamorphism in the KKB. The maximum apparent ²⁰⁷Pb/²⁰⁶Pb age obtained from the unsupported radiogenic Pb concentrations is also consistent with formation of the Pb patches during this event. The superimposed charnockitization event appears to have caused additional Pb-loss in the cores and recrystallization of the rims. The results of depth-profiling of the scanning ion tomography image stack show that the Pb-rich domains range in size from <5 nm to several 10 nm (diameter if assumed to be spherical). The occurrence of such patchy Pb has previously been documented only from UHT metamorphic zircon, where it likely results from annealing of radiation-damaged zircon. The formation of a discrete, heterogeneously distributed and subsequently immobile Pb phase effectively arrests the normal Pb-loss process seen at lower grades of metamorphism.     

Development and sensitivity analysis of a model for assessing stratification and safety of Lake Nyos during artificial degassing

To prevent the recurrence of a disastrous eruption of carbon dioxide (CO₂) from Lake Nyos, a degassing plan has been set up for the lake. Since there are concerns that the degassing of the lake may reduce the stability of the density stratification, there is an urgent need for a simulation tool to predict the evolution of the lake stratification in different scenarios. This paper describes the development of a numerical model to predict the CO₂ and dissolved solids concentrations, and the temperature structure as well as the stability of the water column of Lake Nyos. The model is tested with profiles of CO₂ concentrations and temperature taken in the years 1986 to 1996. It reproduces well the general mixing patterns observed in the lake. However, the intensity of the mixing tends to be overestimated in the epilimnion and underestimated in the monimolimnion. The overestimation of the mixing depth in the epilimnion is caused either by the parameterization of the k-epsilon model, or by the uncertainty in the calculation of the surface heat fluxes. The simulated mixing depth is highly sensitive to the surface heat fluxes, and errors in the mixing depth propagate from one year to the following. A precise simulation of the mixolimnion deepening therefore requires high accuracy in the meteorological forcing and the parameterization of the heat fluxes. Neither the meteorological data nor the formulae for the calculation of the heat fluxes are available with the necessary precision. Consequently, it will be indispensable to consider different forcing scenarios in the safety analysis in order to obtain robust boundary conditions for safe degassing. The input of temperature and CO₂ to the lake bottom can be adequately simulated for the years 1986 to 1996 with a constant sublacustrine source of 18 l s^–1 with a CO₂ concentration of 0.395 mol l^–1 and a temperature of 26 °C. The results of this study indicate that the model needs to be calibrated with more detailed field data before using it for its final purpose: the prediction of the stability and the safety of Lake Nyos during the degassing process.Responsible Editor: Hans Burchard     

2D equivalent linear site effect simulation: example applications to two deep valleys

In the framework of the Sismovalp European project, an equivalent linear 2D code was developed to compute the response of a valley to SH waves, using the discrete wave-number method proposed by Aki and Larner (Aki K, Larner KL (1970) J Geophys Res 75:5). To overcome the frequency upper bound limitation, the Aki and Larner’s method is combined with a one-dimensional computation using a classical multi-layer method (Aki K, Richards PG (1980) Quantitative Seismology: Theory and Methods, vols. 1 & 2. W.H. Freeman & Co, San Francisco). The so-called “Aki–Larner extended method” is associated to an iterative algorithm, as proposed by Seed and Idriss (Seed HB, Idriss IM (1969) Report No. EERC 70–10, Earthquake Research Center, University of California, Berkeley, California) which accounts for the modulus and damping degradation using a linear visco-elastic model. A comparison of the results in the linear and the equivalent linear cases, for a magnitude 6.0 earthquake, shows that the account for the equivalent linear behaviour of the soil significantly reduces the amplification level, especially at frequencies higher than the fundamental resonance frequency of the site. In the case of site effects or microzonation studies devoted to produce design spectra for engineering structures, this can have a major impact on the associated results and costs, depending on the frequency of interest for the considered structure. As a first application of the developed technique, 2D equivalent linear Aki–Larner computations are used to perform the seismic microzonation study of the upper Rhone valley, in the Visp area (Switzerland), a typical 2D alpine valley. These investigations made it possible to determine site specific spectra, associated with different zones, to be used instead of the code spectra that do not take into account the local 2D amplification.     

电离层水平不均匀性对无线电掩星资料反演的影响

电离层水平不均匀结构是引起掩星资料反演的主要误差源.本文模拟不均匀结构对无线电掩星资料反演的影响.首先利用3-D全球射线追踪技术进行了模拟.结果指出:电离层电子密度20%的变化可能引起GPS的时延增加约12m.其次,将90%的掩星时延对应的路径长度,定义为电离层掩星测量的水平分辨率,其值为1200-2000 km.最后,计算了不同纬度电离层当地正午、子夜和日出时的电离层水平不对称指数.它表明:低纬夏季日出时,不对称指数达-78%,反演结果可能有较大误差.     

Further finds of the Derrick Peak meteorite,Transantarctic Mountains,and implications for terrestrial age

Abstract— Nine additional iron meteorite fragments weighing a total of 72 kg were recovered from the Derrick Peak area by a Canterbury Museum geological party in late 1988. One iron was located in the Onnum Valley, 6 km south of the previous finds. Geochemical analysis indicates that all irons belong to a single meteorite shower, greatly increasing the known extent of the fall zone. Kamp and Lowe (1982) have previously estimated the terrestrial age of the meteorite from glacial geological evidence. The location of the 1988 finds supports Kamp and Lowe's interpretation that the meteorites lie in situ, but recent revisions of the chronology of Cenozoic glacial history of the region reduce the interpreted terrestrial age. An age of between Oxygen Isotope stages 6 and 2 is probable (190–125 to 35–12 ka BP). This conflicts with a terrestrial age estimate of 1.0 ± 0.1 Ma BP from cosmogenic radionuclides.