Metamorphic charnockite in contact aureoles around intrusive enderbite from Natal, South Africa

In the Port Edward area of southern Kwa-Zulu Natal, South Africa, charnockitic aureoles up to 10 m in width in the normally garnetiferous Nicholson's Point Granite, are developed adjacent to intrusive contacts with the Port Edward Enderbite and anhydrous pegmatitic veins. Mineralogical differences between the country rock and charnockitic aureole suggest that the dehydration reaction Bt + Qtz → Opx + Kfs + H₂O and the reaction of Grt + Qtz → Opx + Pl were responsible for the charnockitization. The compositions of fluid inclusions show systematic variation with: (1) the Port Edward Enderbite being dominated by CO₂ and N₂ fluid inclusions; (2) the non-charnockitized granite by saline aqueous inclusions with 18–23 EqWt% NaCl; (3) the charnockitic aureoles by low-salinity and pure water inclusions (<7 EqWt% NaCl); (4) the pegmatites by aqueous inclusions of various salinity with minor CO₂. As a result of the thermal event the homogenization temperatures of the inclusions in charnockite show a much larger range (up to 390 °C) compared to the fluid inclusions in granite (mostly <250 °C). Contrary to fluid-controlled charnockitization (brines, CO₂) which may have taken place along shear zones away from the intrusive body, the present “proximal” charnockitized granite formed directly at the contact with enderbite. The inclusions indicate contact metamorphism induced by the intrusion of “dry” enderbitic magma into “wet” granite resulting in local dehydration. This was confirmed by cathodoluminescence microscopy showing textures indicative for the local reduction of structural water in the charnockite quartz. Two-pyroxene thermometry on the Port Edward Enderbite suggests intrusion at temperatures of ∼1000–1050 °C into country rock with temperature of <700 °C. The temperature of aureole formation must have been between ∼700 °C (breakdown of pyrite to form pyrrhotite) and ∼1000 °C. Charnockitization was probably controlled largely by heat related to anhydrous intrusions causing dehydration reactions and resulting in the release and subsequent trapping of dehydration fluids. The salinity of the metamorphic fluid in the contact zones is supposed to have been higher at an early stage of contact metamorphism, but it has lost its salt content by K-metasomatic reactions and/or the preferential migration of the saline fluids out of the contact zones towards the enderbite. The low water activity inhibited the localized melting of the granite. Mineral thermobarometry suggests that after charnockite aureole genesis, an isobaric cooling path was followed during which reequilibration of most of the aqueous inclusions occurred. Received: 8 November 1998 / Accepted: 21 June 1999     

X‐shooting Herbig Ae/Be stars: Accretion probed by near‐infrared He I emission

The Herbig Ae/Be stars are intermediate mass pre‐main sequence stars that bridge the gap between the low mass T Tauri stars and the Massive Young Stellar Objects. In this mass range, the acting star forming mechanism switches from magnetically controlled accretion to an as yet unknown mechanism, but which is likely to be direct disk accretion onto the star. We observed a large sample of Herbig Ae/Be stars with X‐shooter to address this issue from a multi‐wavelength perspective. It is the largest such study to date, not only because of the number of objects involved, but also because of the large wavelength coverage from the blue to the near‐infrared. This allows many accretion diagnostics to be studied simultaneously. By correlating the various properties with mass, temperature and age, we aim to determine where and whether the magnetically controlled mass accretion mechanism halts and the proposed direct disk accretion takes over. Here, we will give an overview of the background, present some observations and discuss our initial results. We will introduce a new accretion diagnostic for the research of Herbig Ae/Be stars, the HeI 1.083 μm line (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

UML‐Based Approach to Developing a CityGML Application Domain Extension

Recently a national 3D standard was established in the Netherlands as a CityGML Application Domain Extension (called IMGeo). In line with the Dutch practice of modeling geo‐information, the ADE is developed using a model driven approach. The classes are designed in UML and automatically mapped to GML schema. The current OGC CityGML specification does not provide rules or guidance on correctly modeling an ADE in UML. This article fills this gap by studying how CityGML can be extended for specific applications starting from the UML diagrams. Six alternatives for modeling ADEs in UML are introduced and compared. The optimal alternative is selected and applied to obtain the national 3D standard. The approach was extensively discussed with international experts, who were members of both SIG3D and other working groups. As a consequence the approach was adopted by the SIG3D, the Special Interest Group 3D which, among other things, work on the 3D standard CityGML in cooperation with OGC. Therefore the approach contains many issues that can be generalized and reused by future domain extensions of CityGML. To further support this, the article formulates a model‐driven framework to model CityGML ADEs. Open issues are described in the conclusions.     

Seismic risk ranking for large dams in South Africa

Earthquake motion is one of the extreme loads acting on large dams. Dam owners and regulators must therefore ensure that dams are safely operated and present minimal risk to the public in case of extreme loads such as floods and earthquakes. Owners of many dams or officials in charge of dam safety programs may consider comparative assessment of the seismic risk associated with their dams and establish priorities for detailed evaluation. South Africa has in excess of 100 large state-owned dams and the characteristics of these dams have been used to perform a basic seismic hazard assessment and rank the vulnerability of these dams from the lowest to highest. One of the most decisive factors that contributes to the risk of a dam is the wall type; with gravity and earthfill dams being the most vulnerable to earthquake motion. Another aspect that needs further investigation is the downstream hazard potential which, if known to a better degree of accuracy, can provide more reasonable estimates of the risk factors.     

River channel and bar patterns explained and predicted by an empirical and a physics‐based method

Our objective is to understand general causes of different river channel patterns. In this paper we compare an empirical stream power‐based classification and a physics‐based bar pattern predictor. We present a careful selection of data from the literature that contains rivers with discharge and median bed particle size ranging over several orders of magnitude with various channel patterns and bar types, but no obvious eroding or aggrading tendency. Empirically a continuum is found for increasing specific stream power, here calculated with pattern‐independent variables: mean annual flood, valley gradient and channel width predicted with a hydraulic geometry relation. ‘Thresholds’, above which certain patterns emerge, were identified as a function of bed sediment size. Bar theory predicts nature and presence of bars and bar mode, here converted to active braiding index (B_i). The most important variables are actual width–depth ratio and nonlinearity of bed sediment transport. Results agree reasonably well with data. Empirical predictions are somewhat better than bar theory predictions, because the bank strength is indirectly included in the empirical prediction. In combination, empirical and theoretical prediction provide partial explanations for bar and channel patterns. Increasing potential‐specific stream power implies more energy to erode banks and indeed correlates to channels with high width–depth ratio. Bar theory predicts that such rivers develop more bars across the width (higher B_i). At the transition from meandering to braiding, weakly braided rivers and meandering rivers with chutes are found. Rivers with extremely low stream power and width–depth ratios hardly develop bars or dynamic meandering and may be straight or sinuous or, in case of disequilibrium sediment feed, anastomosing. Copyright © 2010 John Wiley & Sons, Ltd.     

Thermal equation of state of synthetic orthoferrosilite at lunar pressures and temperatures

Iron-rich orthopyroxene plays an important role in models of the thermal and magmatic evolution of the Moon, but its density at high pressure and high temperature is not well-constrained. We present in situ measurements of the unit-cell volume of a synthetic polycrystalline end-member orthoferrosilite (FeSiO₃, fs) at simultaneous high pressures (3.4–4.8 GPa) and high temperatures (1,148–1,448 K), to improve constraints on the density of orthopyroxene in the lunar interior. Unit-cell volumes were determined through in situ energy-dispersive synchrotron X-ray diffraction in a multi-anvil press, using MgO as a pressure marker. Our volume data were fitted to a high-temperature Birch–Murnaghan equation of state (EoS). Experimental data are reproduced accurately, with a  $\varDelta P$ Δ P  standard deviation of 0.20 GPa. The resulting thermoelastic parameters of fs are: V ₀ = 875.8 ± 1.4 Å³, K ₀ = 74.4 ± 5.3 GPa, and $\frac{{\text d}K}{{\text d}T} = -0.032 \pm 0.005\,\hbox{GPa K}^{-1}$ d K d T = - 0.032 ± 0.005 GPa K - 1 , assuming ${K}^{\prime}_{0} = 10 $ K 0 ′ = 10 . We also determined the thermal equation of state of a natural Fe-rich orthopyroxene from Hidra (Norway) to assess the effect of magnesium on the EoS of iron-rich orthopyroxene. Comparison between our two data sets and literature studies shows good agreement for room-temperature, room-pressure unit-cell volumes. Preliminary thermodynamic analyses of orthoferrosilite, FeSiO₃, and orthopyroxene solid solutions, (Mg_1?x Fe _x ) SiO₃, using vibrational models show that our volume measurements in pressure–temperature space are consistent with previous heat capacity and one-bar volume–temperature measurements. The isothermal bulk modulus at ambient conditions derived from our measurements is smaller than values presented in the literature. This new simultaneous high-pressure, high-temperature data are specifically useful for calculations of the orthopyroxene density in the Moon.     

Iron-hydroxide clogging of public supply wells receiving artificial recharge: near-well and in-well hydrological and hydrochemical observations

Clogging of water wells by iron-hydroxide incrustations due to mixing of anoxic and oxic groundwater is a common well-ageing problem. The relation between well operation (on and off), the spatial and temporal variations in hydrochemistry outside and inside a supply well, and the distribution of clogging iron-hydroxides were studied in an artificial recharge well field in the Netherlands. Camera inspection, high-resolution multi-level water sampling outside the well and detailed in-well pH/EC/O₂ profiles revealed remarkable patterns. During pumping, the top of the upper well screen abstracted oxic filtrate, although the larger part of the in-well water column was anoxic. The column rapidly turned oxic after shutdown due to a downward short-circuiting of oxic water via the well. Within 15 d it became anoxic due to the slow advance of anoxic lake filtrate created by local changes in flow direction as the neighboring wells continued to pump. Severe clogging occurred where the oxic filtrate entered the well, while half-clogging of the upper well screen occurred due to less inflow of oxic filtrate on the lake side. Transport of iron flocs and bacterial slimes after shutdown seemed to clog the lower part of the well screen. Frequent on/off switching should be avoided in iron-clogged wells.     

Joint Consistent Mapping of High-Dimensional Geochemical Surveys

Geochemical surveys often contain several tens of components, obtained from different horizons and with different analytical techniques. These are used either to obtain elemental concentration maps or to explore links between the variables. The first task involves interpolation, the second task principal component analysis (PCA) or a related technique. Interpolation of all geochemical variables (in wt% or ppm) should guarantee consistent results: At any location, all variables must be positive and sum up to 100 %. This is not ensured by any conventional geostatistical technique. Moreover, the maps should ideally preserve any link present in the data. PCA also presents some problems, derived from the spatial dependence between the observations, and the compositional nature of the data. Log-ratio geostatistical techniques offer a consistent solution to all these problems. Variation-variograms are introduced to capture the spatial dependence structure: These are direct variograms of all possible log ratios of two components. They can be modeled with a function analogous to the linear model of coregionalization (LMC), where for each spatial structure there is an associated variation matrix describing the links between the components. Eigenvalue decompositions of these matrices provide a PCA of that particular spatial scale. The whole data set can then be interpolated by cokriging. Factorial cokriging can also be used to map a certain spatial structure, eventually projected onto those principal components (PCs) of that structure with relevant contribution to the spatial variability. If only one PC is used for a certain structure, the maps obtained represent the spatial variability of a geochemical link between the variables. These procedures and their advantages are illustrated with the horizon C Kola data set, with 25 components and 605 samples covering most of the Kola peninsula (Finland, Norway, Russia).     

Using a groundwater quality negotiation support system to change land-use management near a drinking-water abstraction in the Netherlands

A negotiation support system (NSS) was developed to solve groundwater conflicts that arose during land-use management. It was set up in cooperation with the stakeholders involved to provide information on the impact of land use, e.g., agriculture, nature (forested areas), recreation, and urban areas, on the quality of both infiltrating and abstracted groundwater. This NSS combined simulation programs that calculate (1) the concentrations of nitrate in shallow groundwater for each land-use area and (2) the transport of nitrate in the groundwater-saturated zone. The user interface of the NSS enabled scenario analyses. The NSS was validated at a drinking-water abstraction near Holten (the Netherlands) using a spatial planning process aimed at sustainable land-use and groundwater-resource management. Two land-use scenarios were considered: a base scenario reflecting the autonomous development and an adapted land-use scenario. The calculated results for shallow groundwater provided an explicit spatial overview of the impact of historical land use and N application on the quality of abstracted groundwater as well as insight into the impact of changes in land use and N application. Visualization of the conflicting interests of agriculture and the drinking-water abstraction helped all stakeholders accept the necessary changes in land use identified by the adapted land-use scenario of the NSS. These changes were included in the preferred land-use management option in the regional planning process, which has since been formalized. The NSS provided system insight, scoping analyses, and education, in addition to generating quantitative information on the impact of land-use functions on groundwater quality.