Schorlomite: a discussion of the crystal chemistry, formula, and inter-species boundaries

Examination of schorlomite from ijolite at Magnet Cove (USA) and silicocarbonatite at Afrikanda (Russia), using electron-microprobe and hydrogen analyses, X-ray diffraction and Mössbauer spectroscopy, shows the complexity of substitution mechanisms operating in Ti-rich garnets. These substitutions involve incorporation of Na in the eightfold-coordinated X site, Fe²⁺ and Mg in the octahedrally coordinated Y site, and Fe³⁺, Al and Fe²⁺ in the tetrahedrally coordinated Z site. Substitutions Ti⁴⁺Fe³⁺Fe³⁺_–1Si_–1 and Ti⁴⁺Al³⁺Fe³⁺_–1Si_–1 are of major significance to the crystal chemistry of schorlomite, whereas Fe²⁺ enters the Z site in relatively minor quantities (<3% Fe). There is no evidence (either structural or indirect, such as discrepancies between the measured and calculated Fe²⁺ contents) for the presence of ^[6]Ti³⁺ or ^[4]Ti⁴⁺ in schorlomite. The simplified general formula of schorlomite can be written as Ca₃Ti⁴⁺₂[Si_3-x(Fe³⁺,Al,Fe²⁺)_xO₁₂], keeping in mind that the notion of end-member composition is inapplicable to this mineral. In the published analyses of schorlomite with low to moderate Zr contents, x ranges from 0.6 to 1.0, i.e. Ti⁴⁺ in the Y site is <2 and accompanied by appreciable amounts of lower-charged cations (in particular, Fe³⁺, Fe²⁺ and Mg). For classification purposes, the mole percentage of schorlomite can be determined as the amount of ^[6]Ti⁴⁺, balanced by substitutions in the Z site, relative to the total occupancy in the Y site: (^[6]Ti⁴⁺–^[6]Fe²⁺–^[6]Mg²⁺– ^[8]Na⁺)/2. In addition to the predominant schorlomite component, the crystals examined in this work contain significant (>15 mol.%) proportions of andradite (Ca₃Fe³⁺₂Si₃O₁₂), morimotoite (Ca₃Fe²⁺TiSi₃O₁₂), and Ca₃MgTiSi₃O₁₂. The importance of accurate quantitative determination and assignment of Fe, Ti and other cations to the crystallographic sites for petrogenetic studies is discussed.

The hard solar X-ray burst of 18 September 1963

A hard solar X-ray burst was observed by J-P. Legrand on 18 September 1963, 13:56 UT, at balloon altitude. It lasted a few minutes; a steep increase was followed by an exponential decay. During its declining phase a weak radio burst was observed on 3 and 10 cm, not on longer wavelengths.Maximum radio intensity occurred two minutes after that of the X-ray burst. The X-ray and radio bursts ended almost simultaneously. Optically a small shortlived (some minutes) flare point occurred simultaneously with the X-ray burst in a magnetically interesting part of the active region of September 1963. The X-burst photons seem to have had an energy of about 0.5 MeV. The burst was therefore of a fairly rare type, since very few other bursts with similar photon energies have been detected up to now.It is suggested that a mass of gas, magnetically confined to a volume of about 5·10²⁵ cm³ in the low corona, containing about 3·10³⁵ electrons was accelerated to energies of about 0.5 MeV. The gas gradually expanded, partly also to higher levels. The gyro-synchrotron radiation, emitted by the plasma became observable after about two minutes. At the lower radio frequencies the radiation was absorbed by overlying undisturbed coronal matter. Quantitative computations justify this model. A detailed summary of the events, and some numerical data are given in the concluding Section 8 and in Table V.     

Validation of an annual cycle simulation with a T42-L20 GCM

An annual cycle of an atmospheric general circulation model (AGCM) is presented. The winter and summer zonal averages of the atmospheric fields are compared with an observed climatology. The main features of the observed seasonal means are well reproduced by the model. One of the main discrepancies is that the simulated atmosphere is too cold, particularly in its upper part. Some other discrepancies might be explained by the interannual variability. The AGCM surface fluxes are directly compared to climatological estimates. On the other hand, the calculation of meridional heat transport by the ocean, inferred from the simulated energy budget, can be compared to transport induced from climatologies. The main result of this double comparison is that AGCM fluxes generally are within the range of climatological estimates. The main deficiency of the model is poor partitioning between solar and non-solar heat fluxes in the tropical belt. The meridional heat transport also reveals a significant energy-loss by the Northern Hemisphere ocean north of 45° N. The possible implications of model surface flux deficiencies on coupling with an oceanic model are discussed.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. Dümenil     

Invariant manifolds, phase correlations of chaotic orbits and the spiral structure of galaxies

In the presence of a strong   m = 2  component in a rotating galaxy, the phase-space structure near corotation is shaped to a large extent by the invariant manifolds of the short-period family of unstable periodic orbits terminating at L ₁ or L ₂. The main effect of these manifolds is to create robust phase correlations among a number of chaotic orbits large enough to support a spiral density wave outside corotation. The phenomenon is described theoretically by soliton-like solutions of a Sine–Gordon equation. Numerical examples are given in an N -body simulation of a barred spiral galaxy. In these examples, we demonstrate how the projection of unstable manifolds in configuration space reproduces essentially the entire observed bar–spiral pattern.     

Fast integration of low orbiter's trajectory perturbed by the earth's non-sphericity

A fast algorithm is proposed to integrate the trajectory of a low obiter perturbed by the earth's non-sphericity. The algorithm uses a separation degree to define the low-degree and the high-degree acceleration components, the former computed rigorously, and the latter interpolated from gridded accelerations. An FFT method is used to grid the accelerations. An optimal grid type for the algorithm depends on the trajectory's permissible error, speed, and memory capacity. Using the non-spherical accelerations computed from EGM96 to harmonic degree 360, orbit integrations were performed for a low orbiter at an altitude of 170 km. For a separation degree of 50, the new algorithm, together with the predict-pseudo correct method, speeds up the integration by 145 times compared to the conventional algorithm while keeping the errors in position and velocity below 10⁻⁴ m and 10⁻⁷ m/s for a 3-day arc. Received: 28 July 1997 / Accepted: 1 April 1998     

Reference frame stability and nonlinear distortion in minimum-constrained network adjustment

The aim of this paper is to investigate the influence of the minimum constraints (MCs) on the reference frame parameters in a free-net solution. The non-estimable part of these parameters (which is not reduced by the available data) is analysed in terms of its stability under a numerical perturbation of the constrained datum functionals. In practice, such a perturbation can be ascribed either to hidden errors in the known coordinates/velocities that participate in the MCs or to a simple change of their a priori values due to a datum switch on a different fiducial dataset. In addition, a perturbation of this type may cause a nonlinear variation to the estimable part of the reference frame parameters, since it theoretically affects the adjusted observations that are implied by the network’s nonlinear observational model. The aforementioned effects have an impact on the quality of a terrestrial reference frame (TRF) that is established via a minimum-constrained adjustment, and our study shows that they are both controlled through a characteristic matrix which is inherently linked to the MC system (the so-called TRF stability matrix). The structure of this matrix depends on the network’s spatial configuration and the ‘geometry’ of the datum constraints, while its main role is the filtering of any MC-related errors into the least-squares adjustment results. A number of examples with different types of geodetic networks are also presented to demonstrate the theoretical findings of our study.     

Applications of AnnAGNPS model for soil loss estimation and nutrient loading for Malaysian conditions

The study evaluated the performance and suitability of AnnAGNPS model in assessing runoff, sediment loading and nutrient loading under Malaysian conditions. The watershed of River Kuala Tasik in Malaysia, a combination of two sub-watersheds, was selected as the area of study. The data for the year 2004 was used to calibrate the model and the data for the year 2005 was used for validation purposes. Several input parameters were computed using methods suggested by other researchers and studies carried out in Malaysia. The study shows that runoff was predicted well with an overall R² value of 0.90 and E value of 0.70. Sediment loading was able to produce a moderate result of R² = 0.66 and E = 0.49, nitrogen loading predictions were slightly better with R² = 0.68 and E = 0.53, and phosphorus loading performance was slightly poor with an R² = 0.63 and E = 0.33. The erosion map developed was in agreement with the erosion risk map produced by the Department of Agriculture, Malaysia. Rubber estates and urban areas were found to be the main contributors to soil erosion. The simulation results showed that AnnAGNPS has the potential to be used as a valuable tool for planning and management of watersheds under Malaysian conditions.     

Comparing spatially varying coefficient models: a case study examining violent crime rates and their relationships to alcohol outlets and illegal drug arrests

In this paper, we compare and contrast a Bayesian spatially varying coefficient process (SVCP) model with a geographically weighted regression (GWR) model for the estimation of the potentially spatially varying regression effects of alcohol outlets and illegal drug activity on violent crime in Houston, Texas. In addition, we focus on the inherent coefficient shrinkage properties of the Bayesian SVCP model as a way to address increased coefficient variance that follows from collinearity in GWR models. We outline the advantages of the Bayesian model in terms of reducing inflated coefficient variance, enhanced model flexibility, and more formal measuring of model uncertainty for prediction. We find spatially varying effects for alcohol outlets and drug violations, but the amount of variation depends on the type of model used. For the Bayesian model, this variation is controllable through the amount of prior influence placed on the variance of the coefficients. For example, the spatial pattern of coefficients is similar for the GWR and Bayesian models when a relatively large prior variance is used in the Bayesian model.      

Reducing the draconitic errors in GNSS geodetic products

Systematic errors at harmonics of the GPS draconitic year have been found in diverse GPS-derived geodetic products like the geocenter $Z$ -component, station coordinates, $Y$ -pole rate and orbits (i.e. orbit overlaps). The GPS draconitic year is the repeat period of the GPS constellation w.r.t. the Sun which is about 351 days. Different error sources have been proposed which could generate these spurious signals at the draconitic harmonics. In this study, we focus on one of these error sources, namely the radiation pressure orbit modeling deficiencies. For this purpose, three GPS+GLONASS solutions of 8 years (2004–2011) were computed which differ only in the solar radiation pressure (SRP) and satellite attitude models. The models employed in the solutions are: (1) the CODE (5-parameter) radiation pressure model widely used within the International GNSS Service community, (2) the adjustable box-wing model for SRP impacting GPS (and GLONASS) satellites, and (3) the adjustable box-wing model upgraded to use non-nominal yaw attitude, specially for satellites in eclipse seasons. When comparing the first solution with the third one we achieved the following in the GNSS geodetic products. Orbits: the draconitic errors in the orbit overlaps are reduced for the GPS satellites in all the harmonics on average 46, 38 and 57 % for the radial, along-track and cross-track components, while for GLONASS satellites they are mainly reduced in the cross-track component by 39 %. Geocenter $Z$ -component: all the odd draconitic harmonics found when the CODE model is used show a very important reduction (almost disappearing with a 92 % average reduction) with the new radiation pressure models. Earth orientation parameters: the draconitic errors are reduced for the $X$ -pole rate and especially for the $Y$ -pole rate by 24 and 50 % respectively. Station coordinates: all the draconitic harmonics (except the 2nd harmonic in the North component) are reduced in the North, East and Height components, with average reductions of 41, 39 and 35 % respectively. This shows, that part of the draconitic errors currently found in GNSS geodetic products are definitely induced by the CODE radiation pressure orbit modeling deficiencies.