Olivine in the Udachnaya-East Kimberlite (Yakutia, Russia): Types, Compositions and Origins

Olivine is the principal mineral of kimberlite magmas, and isthe main contributor to the ultramafic composition of kimberliterocks. Olivine is partly or completely altered in common kimberlites,and thus unavailable for studies of the origin and evolutionof kimberlite magmas. The masking effects of alteration, commonin kimberlites worldwide, are overcome in this study of theexceptionally fresh diamondiferous kimberlites of the Udachnaya-Eastpipe from the Daldyn–Alakit province, Yakutia, northernSiberia. These serpentine-free kimberlites contain large amountsof olivine (50 vol.%) in a chloride–carbonate groundmass.Olivine is represented by two populations (olivine-I and groundmassolivine-II) differing in morphology, colour and grain size,and trapped mineral and melt inclusions. The large fragmentalolivine-I is compositionally variable in terms of major (Fo_85–94)and trace element concentrations, including H₂O content (10–136ppm). Multiple sources of olivine-I, such as convecting andlithospheric mantle, are suggested. The groundmass olivine-IIis recognized by smaller grain sizes and perfect crystallographicshapes that indicate crystallization during magma ascent andemplacement. However, a simple crystallization history for olivine-IIis complicated by complex zoning in terms of Fo values and traceelement contents. The cores of olivine-II are compositionallysimilar to olivine-I, which suggests a genetic link betweenthese two types of olivine. Olivine-I and olivine-II have oxygenisotope values (+ 5·6 ± 0·1 VSMOW, 1 SD)that are indistinguishable from one another, but higher thanvalues (+ 5·18 ± 0·28) in ‘typical’mantle olivine. These elevated values probably reflect equilibriumwith the Udachnaya carbonate melt at low temperatures and ¹⁸O-enrichedmantle source. The volumetrically significant rims of olivine-IIhave constant Fo values (89·0 ± 0·2 mol%),but variable trace element compositions. The uniform Fo compositionsof the rims imply an absence of fractionation of the melt'sFe²⁺/Mg, which is possible in the carbonatite melt–olivinesystem. The kimberlite melt is argued to have originated inthe mantle as a chloride–carbonate liquid, devoid of ‘ultramafic’or ‘basaltic’ aluminosilicate components, but becameolivine-laden and olivine-saturated by scavenging olivine crystalsfrom the pathway rocks and dissolving them en route to the surface.During emplacement the kimberlite magma changed progressivelytowards an original alkali-rich chloride–carbonate meltby extensively crystallizing groundmass olivine and gravitationalseparation of solids in the pipe. KEY WORDS: kimberlite; olivine; partial melting; carbonatitic melt; oxygen isotopes; H₂O     

Shrinking Sphere Kinetics for Batch Dissolution of Mixed Particles of a Single Substance at High Under-Saturation: Validation with Sodium Chloride, but with Biogenic Silica in Mind

The cubic equation recently derived for the increase in concentration of a solute with time, as the solid dissolves in batch according to the shrinking sphere model at high under-saturation, is extended to dissolutions of mixtures of differently sized particles. This problem needs to be solved if batch dissolutions are to play their part in the proposed amelioration of global warming and associated climate change by accelerated ‘re-burial’ of excess CO₂ in ocean sediment. The upgraded model was tested using sodium chloride dissolved in 50% aqueous propanone, whence the model fitted two separate runs with 500 and 212 μm, and 212 and 38 μm, diameter crystals, respectively. The key to simulating dissolution in this way lies in the dissolutions being independent of each other. It is further shown that although this condition was implicit in the recent derivation of the cubic equation, it was not recognised at the time. The work should be applicable to any batch dissolution of mixed particles at high under-saturation, and hence, may find use in many industrial and laboratory dissolutions. Simulations show how agglomerated mixtures can yield a straight line on the plot of ln(1 − C/C _T) versus time, as was reported to occur recently with sodium chloride taken ‘straight from the bottle’. It is shown that this probably explains why exponential dissolutions may have seemed appropriate to the dissolution of biogenic silica in earlier literature. This study suggests that a new round of biogenic silica dissolutions, but with sieved samples, would be worthwhile, with the likelihood that shrinking sphere behaviour might well be found to characterise the kinetics. The opportunity is taken to investigate a number of aspects of the shrinking sphere model not generally discussed before, e.g. the graph for the change in surface area with time. The limitations of using cubic salt crystals with the shrinking sphere model are discussed.     

Unravelling the widening of the earliest Andean northern orogen: Maastrichtian to early Eocene intra‐basinal deformation in the northern Eastern Cordillera of Colombia

The onset of deformation in the northern Andes is overprinted by subsequent stages of basin deformation, complicating the examination of competing models illustrating potential location of earliest synorogenic basins and uplifts. To establish the width of the earliest northern Andean orogen, we carried out field mapping, palynological dating, sedimentary, stratigraphic and provenance analyses in Campanian to lower Eocene units exposed in the northern Eastern Cordillera of Colombia (Cocuy region) and compare the results with coeval succession in adjacent basins. The onset of deformation is recorded in earliest Maastrichtian time, as terrigenous detritus arrived into the basin marking the end of chemical precipitation and the onset of clastic deposition produced by the uplift of a western source area dominated by shaly Cretaceous rocks. Disconformable contacts within the upper Maastrichtian to middle Palaeocene succession document increasing supply of quartzose sandy detritus from Cretaceous quartzose rocks exposed in eastern source areas. The continued unroofing of both source areas produced a rapid shift in depositional environments from shallow marine in Maastrichtian to fluvial‐lacustrine systems during the Palaeocene‐early Eocene. Supply of immature Jurassic sandstones from nearby western uplifts, together with localized plutonic and volcanic Cretaceous rocks, caused a shift in Palaeocene sandstones composition from quartzarenites to litharenites. Supply of detrital sandy fragments, unstable heavy minerals and Cretaceous to Ordovician detrital zircons, were derived from nearby uplifted blocks and from SW fluvial systems within the synorogenic basin, instead of distal basement rocks. The presence of volcanic rock fragments and 51–59 Ma volcanic zircons constrain magmatism within the basin. The Maastrichtian–Palaeocene sequence studied here documents crustal deformation that correlates with coeval deformation farther south in Ecuador and Peru. Slab flattening of the subducting Caribbean plate produced a wider orogen (>400 km) with a continental magmatic arc and intra‐basinal deformation and magmatism.     

Centre for Baltic and arctic studies at St. Petersburg state University as a naional node of “University GIS network” in Russia

The possibilities of recent advantages in Geographical Information Systems (GIS) technology as a tool for multidisciplinary coastal zone research are discussed. International education in the new high technologies is one of the most urgent tasks for the University GIS Network.     

Archean of Greenland and Fennoscandia

The North Atlantic craton in southern West Greenland mainly consists of a tectonic collage of Mesoarchean continental crustal terranes, which were amalgamated at c. 2.7 Ga and are currently exposed at mid-crustal amphibolite to granulite facies levels. Tonalitic orthogneisses predominate, intercalated with slightly older tholeiitic to andesitic metavolcanic rocks and associated gabbro-anorthosite intrusive complexes. The North Atlantic craton also contains enclaves of Eoarchean, c. 3.86-3.6 Ga orthogneisses and supracrustal rocks including the Isua greenstone （or supracrustal） belt. This is the oldest known assemblage of rocks deposited at the surface of the Earth, comprising mafic pillow lavas, banded iron formations and metasedimentary schists with local disseminated graphite of possible biogenic origin. Eoarchean rocks have not been found in Kola and Karelia in Fennoscandia where most rocks are 2.9-2.7 Ga tonalitic-trondhjemitic-granodioritic orthogneisses with intercalated coeval greenstone belts and amphibolites. Mesoarchean 3.0-3.2 Ga rocks are found in the eastern and western parts of the Karelian province. Subduction-related rocks like the Iringora supra-subduction type ophiolite and basalt-andesite-dacite-rhyolite series volcanic rocks in many greenstone belts, as well as eclogites are found in the Archean of Fennoscandia. A clear distinction between Greenland and Fennoscandia is the abundance of 2.75-2.65 Ga igneous rocks in Fennoscandia which indicates that these two cratons had a separate evolution during the Neoarchean.     

Variably saturated groundwater modelling for optimizing managed aquifer recharge using trench infiltration

Spreading‐basin methods have resulted in more than 130 million cubic metres of recharge to the unconfined Navajo Sandstone of southern Utah in the past decade, but infiltration rates have slowed in recent years because of reduced hydraulic gradients and (or) clogging. Trench infiltration is a promising alternative technique for increasing recharge and minimizing evaporation. This paper uses a variably saturated flow model to further investigate the relative importance of the following variables on rates of trench infiltration to unconfined aquifers: saturated hydraulic conductivity, trench spacing and dimensions, initial water‐table depth, alternate wet/dry periods, and number of parallel trenches. Modelling results showed (1) increased infiltration with higher hydraulic conductivity, deeper initial water tables, and larger spacing between parallel trenches, (2) deeper or wider trenches do not substantially increase infiltration, (3) alternating wet/dry periods result in less overall infiltration than keeping the trenches continuously full, and (4) larger numbers of parallel trenches within a fixed area increases infiltration but with a diminishing effect as trench spacing becomes tighter. An empirical equation for estimating expected trench infiltration rates as a function of hydraulic conductivity and initial water‐table depth was derived and can be used for evaluating feasibility of trench infiltration in other hydrogeologic settings. Copyright © 2014 John Wiley & Sons, Ltd.     

Selecting best mapping strategies for storm runoff modeling in a mountainous semi‐arid area

Accurate runoff and soil erosion modeling is constrained by data availability, particularly for physically based models such as OpenLISEM that are data demanding, as the processes are calculated on a cell‐by‐cell basis. The first decision when using such models is to select mapping units that best reflect the spatial variability of the soil and hydraulic properties in the catchment. In environments with limited data, available maps are usually generic, with large units that may lump together the values of the soil properties, affecting the spatial patterns of the predictions and output values in the outlet. Conversely, the output results may be equally acceptable, following the principle of equifinality. To studyhow the mapping method selected affects the model outputs, four types of input maps with different degrees of complexity were created: average values allocated to general soil map units (ASG1), average values allocated to detailed map units (ASG2), values interpolated by ordinary kriging (OK) and interpolated by kriging with external drift (KED). The study area was Ribeira Seca, a 90 km² catchment located in Santiago Island, Cape Verde (West Africa), a semi‐arid country subject to scarce but extreme rainfall during the short tropical summer monsoon. To evaluate the influence of rainfall on runoff and erosion, two storm events with different intensity and duration were considered. OK and KED inputs produced similar results, with the latter being closer to the observed hydrographs. The highest soil losses were obtained with KED (43 ton ha^{? 1} for the strongest event). To improve the results of soil loss predictions, higher accurate spatial information on the processes is needed; however, spatial information of input soil properties alone is not enough in complex landscapes. The results demonstrate the importance of selecting the appropriate mapping strategy to obtain reliable runoff and erosion estimates. Copyright © 2013 John Wiley & Sons, Ltd.     

Rebuttal to comment on “Modeling of opposition effects with ensembles of clusters: Interplay of various scattering mechanisms” by Elena V. Petrova, Victor P. Tishkovets, Klaus Jockers, 2007 [Icarus 188, 233-245]

Shkuratov and Zubko [Shkuratov, Yu.G., Zubko, E., 2008. Icarus 194, 850-852] criticize our paper [Petrova, E.V., Tishkovets, V.P., Jockers, K., 2007. Icarus 188, 233-245]. With this comment we reply to this criticism. We show that the experimental data and the modeling calculations presented by these authors cannot disprove the near-field effect as an important contributor to the scattering mechanisms considered in our paper.     

Global m= 1 instabilities and lopsidedness in disc galaxies

Lopsidedness is common in spiral galaxies. Often, there is no obvious external cause, such as an interaction with a nearby galaxy, for such features. Alternatively, the lopsidedness may have an internal cause, such as a dynamical instability. In order to explore this idea, we have developed a computer code that searches for self-consistent perturbations in razor-thin disc galaxies and performed a thorough mode-analysis of a suite of dynamical models for disc galaxies embedded in an inert dark matter halo with varying amounts of rotation and radial anisotropy.
Models with two equal-mass counter-rotating discs and fully rotating models both show growing lopsided modes. For the counter-rotating models, this is the well-known counter-rotating instability, becoming weaker as the net rotation increases. The m = 1 mode of the maximally rotating models, on the other hand, becomes stronger with increasing net rotation. This rotating m = 1 mode is reminiscent of the eccentricity instability in near-Keplerian discs.
To unravel the physical origin of these two different m = 1 instabilities, we studied the individual stellar orbits in the perturbed potential and found that the presence of the perturbation gives rise to a very rich orbital behaviour. In the linear regime, both instabilities are supported by aligned loop orbits. In the non-linear regime, other orbit families exist that can help support the modes. In terms of density waves, the counter-rotating m = 1 mode is due to a purely growing Jeans-type instability. The rotating m = 1 mode, on the other hand, grows as a result of the swing amplifier working inside the resonance cavity that extends from the disc centre out to the radius where non-rotating waves are stabilized by the model's outwardly rising Q profile.