Compressibility and crystal structure of sillimanite, Al2SiO5, at high pressure

The unit-cell dimensions and crystal structure of sillimanite at various pressures up to 5.29 GPa have been refined from single-crystal X-ray diffraction data. As pressure increases, a and b decrease linearly, whereas c decreases nonlinearly with a slightly positive curvature. The axial compression ratios at room pressure are β_a:β_b:β_c=1.22:1.63:1.00. Sillimanite exhibits the least compressibility along c, but the least thermal expansivity along a (Skinner et al. 1961; Winter and Ghose 1979). The bulk modulus of sillimanite is 171(1) GPa with K′=4 (3), larger than that of andalusite (151 GPa), but smaller than that of kyanite (193 GPa). The bulk moduli of the [Al1O₆], [Al2O₄], and [SiO₄] polyhedra are 162(8), 269(33), and 367(89) GPa, respectively. Comparison of high-pressure data for Al₂SiO₅ polymorphs reveals that the [SiO₄] tetrahedra are the most rigid units in all these polymorphic structures, whereas the [AlO₆] octahedra are most compressible. Furthermore, [AlO₆] octahedral compressibilities decrease from kyanite to sillimanite, to andalusite, the same order as their bulk moduli, suggesting that [AlO₆] octahedra control the compression of the Al₂SiO₅ polymorphs. The compression of the [Al1O₆] octahedron in sillimanite is anisotropic with the longest Al1-OD bond shortening by ～1.9% between room pressure and 5.29 GPa and the shortest Al1-OB bond by only 0.3%. The compression anisotropy of sillimanite is primarily a consequence of its topological anisotropy, coupled with the compression anisotropy of the Al-O bonds within the [Al1O₆] octahedron.     

Uncertainties in field-line tracing in the magnetosphere. Part II: the complete internal geomagnetic field

The discussion in the preceding paper is restricted to the uncertainties in magnetic-field-iine tracing in the magnetosphere resulting from published standard errors in the spherical harmonic coefficients that define the axisymmetric part of the internal geomagnetic field (i.e. g_n⁰ ± g_n⁰). Numerical estimates of these uncertainties based on an analytic equation for axisymmetric field lines are in excellent agreement with independent computational estimates based on stepwise numerical integration along magnetic field lines. This comparison confirms the accuracy of the computer program used in the present paper to estimate the uncertainties in magnetic-field-line tracing that arise from published standard errors in the full set of spherical harmonic coefficients, which define the complete (non-axisymmetric) internal geomagnetic field (i.e. g_n^m ± g_n^m and h_n^m ± h_n^m). An algorithm is formulated that greatly reduces the computing time required to estimate these uncertainties in magnetic-field-line tracing. The validity of this algorithm is checked numerically for both the axisymmetric part of the internal geomagnetic field in the general case (1 n 10) and the complete internal geomagnetic field in a restrictive case (0 m n, 1 n 3). On this basis it is assumed that the algorithm can be used with confidence in those cases for which the computing time would otherwise be prohibitively long. For the complete internal geomagnetic field, the maximum characteristic uncertainty in the geocentric distance of a field line that crosses the geomagnetic equator at a nominal dipolar distance of 2 R_E is typically 100 km. The corresponding characteristic uncertainty for a field line that crosses the geomagnetic equator at a nominal dipolar distance of 6 R_E is typically 500 km. Histograms and scatter plots showing the characteristic uncertainties associated with magnetic-field-line tracing in the magnetosphere are presented for a range of illustrative examples. Finally, estimates are given for the maximum uncertainties in the locations of the conjugate points of selected geophysical observatories. Numerical estimates of the uncertainties in magnetic-field-line tracing in the magnetosphere, including the associated uncertainties in thelocations of the conjugate points of geophysical observatories, should be regarded as first approximations in the sense that these estimates are only as accurate as the published standard errors in the full set of spherical harmomic coefficients. As in the preceding paper, howerver, all computational techniques developed in this paper can be used to derive more realistic estimates of the uncertainties in magnetic-field-line tracing in the magnetosphere, following further progress in the determination of more accurate standard errors in the spherical harmonic coefficients.Also Visiting Reader in Physics, University of Sussex, Palmer, Brighton, BN1 9QH, UK     

Paleomagnetism of the Esterel rocks: a revisit 22 years after the thesis of Hans Zijderveld

For his PhD. thesis, Zijderveld (1975) studied the paleomagnetism of the Permian Esterel rocks (southern France). High-quality thermal and alternating-field demagnetization diagrams were interpreted to determine the direction of the characteristic natural magnetization. For the Esterel volcanics, a mean direction of Dec = 206.5°, Inc = –23°, ₉₅ = 5.7°, k = 112 was found for this magnetization. The dispersion in this mean is remarkably low. Only the declination of the Reyran Rhyolite in the Reyran River quarry clearly deviated from this mean. This deviating direction is not found in our samples, taken at the same site. As many faults occur in this quarry, it is suggested that Zijderveld sampled this rhyolite on a small rotated block. To verify whether the small dispersion in the mean paleomagnetic direction of the Esterel rocks has a geomagnetic or a rock-magnetic origin, two conglomerate tests were carried out. One of these might be interpreted as positive. The results of the other conglomerate test (Agay Formation) are ambiguous: four of the six measured boulders show directions close to the mean paleomagnetic direction of the Esterel rocks. Rock-magnetic measurements show that the remanence is carried by a magnetite and a hematite fraction. The low dispersion in the paleomagnetic directions, the conglomerate tests, and hematite as remanence carrier suggest that the characteristic remanence in the Esterel volcanics was not instantaneously acquired during cooling, but might be affected by remagnetization due to weathering.     

Investigation by backscatter radar of artificial irregularities produced in ionospheric plasma heating experiments

The artificial ionospheric turbulence occurs in the ionosphere illuminated by high power HF radio waves. There are a lot of irregularities stretched along the geomagnetic field in this region. The investigation of the artificially disturbed ionospheric region is based on the reception of back scattered signals (BSS) which permits the basic parameters of this region to be estimated and its inhomogeneous structure to be described.Experiments were carried out using ‘Sura’ heating facility in the frequency range of ? = 4.7–9 MHz (ordinary mode) with the effective radiated power P_eff = 50–70 MW beamed vertically upwards. The most important dynamic and statistical BSS characteristics (the built-up time, the relaxation and autocorrelation times, the BSS amplitude spatial correlation function and power spectra) were measured using probe waves in the frequency range of that made it possible to obtain the basic parameters of the artificial irregularities. The model representation of a disturbed region in a form of a periodic structure gives a possibity to evaluate the scale of the structure, the whole size of the disturbance and its power and to calculate the main BSS characteristics.     

Biomass burning in the tropical savannas of Ivory Coast: An overview of the field experiment Fire of Savannas (FOS/DECAFE 91)

FOS/DECAFE 91 (Fire of Savannas/Dynamique et Chimie Atmosphérique en Forêt Equatoriale) was the first multidisciplinary experiment organized in Africa to determine gas and aerosol emissions by prescribed savanna fires. The humid savanna of Lamto in Ivory Coast was chosen for its ecological characteristics representative of savannas with a high biomass density (900 g m^–2 dry matter). Moreover the vegetation and the climate of Lamto have been studied for more than twenty years. The emission ratios (X/CO₂) of the carbon compounds (CO₂, CO, NMHC, CH₄, PAH, organic acids and aerosols), nitrogen compounds (NO_x, N₂O, NH₃ and soluble aerosols) and sulfur compounds (SO₂, COS and aerosols) were experimentally determined by ground and aircraft measurements. To perform this experiment, 4 small plots (100×100 m) and 2 large areas (10×10 km) were prepared and burnt in January 1991 during the period of maximum occurrence of fires in this type of savanna. The detailed ecological study shows that the carbon content of the vegetation is constant within 1% (42 g C for 100 g of vegetal dry matter), the nitrogen content (0.29 g N for 100 g of dry matter) may vary by 10% and the sulfur content (0.05 g S/100 d.m.) by 20%. These variations of the biomass chemical content do not constitute an important factor in the variation of the gas and particle emission levels. With the emission ratios characteristic of humid savanna and flaming conditions (CO/CO₂ of 6.1% at the ground and 8% for airborne measurements), we propose a set of new emission factors, taking into account the burning efficiency which is about 80%: 74.4% of the carbon content of the savanna biomass is released to the atmosphere in the form of CO₂, 4.6% as CO, 0.2% as CH₄, 0.5% as NMHC and 0.7% as aerosols. 17.2% of the nitrogen content of the biomass is released as NO_x, 3.5% as N₂O, 0.6% as NH₃ and 0.5% as soluble aerosols.     

The mass accommodation coefficient of ammonia on water

The mass accommodation coefficient of ammonia gas on water has been determined by measuring the absorption rate of 50–200 ppm NH₃ in one atm of air or helium into a liquid jet of 97 µm diameter as function of the exposed jet length, and comparing the results with numerical simulations which treat as the only free parameter. The model considers in detail transport of NH₃ by molecular diffusion, penetration of the gas/water interface, hydrolysis in the acidified water, and transport of the solutes from the surface into the jet. A correction is applied for the time evolution of the jet surface speed, using literature data on the fluid mechanics of liquid jets. The result of nine sets of independent measurements is

On the bogussing of tropical cyclones in numerical models: A comparison of vortex profiles

Summary At the resolutions currently in use, and with the sparse oceanic data coverage, numerical analyses cannot adequately represent tropical cyclone circulations for use in numerical weather prediction models. In many cases there is no circulation present at all. Most numerical weather prediction centers therefore employ a bogussing scheme to force a tropical cyclone vortex into the numerical analysis. The standard procedure is to define a synthetic data distribution based on an analytically prescribed vortex, which is passed to the analysis scheme as a set of high quality observations.In this study, four commonly used bogus vortices are examined by comparing resultant forecast tracks in an environment at rest, and in a background flow that simulates a typical monsoon trough-subtropical ridge structure. There are three main findings, each of which has significance for operational tropical cyclone track prediction. First, great care is needed in the choice of the characteristics of the bogus vortex, such as the radius and magnitude of the maximum wind. Second, the tropical cyclone trajectories can be very sensitive to their initial position in the idealised environment. Third, the bogus vortex can substantially influence the environment, especially over longer time periods and for vortices of larger size.With 9 Figures     

Daily rainfall domains in Mallorca

Summary A Pearson Product Moment Correlation matrix for four years of daily rainfall from 96 sites on the island of Mallorca has been used to derive and map correlation fields and construct composite correlation fields for different parts of the island based on previously defined daily rainfall affinity areas. Significant variations in rainfall organisation over the island are revealed. These reflect the dominance of uplift and exposure by upland areas in the generation of rainfall, evidenced in strongly anisotropic correlation fields. Further anisotropy along lowland coasts normal to the relief trend indicate the role of local rainfall generation along sea-breeze fronts. Near isotropic fields for central island locations highlight the importance of local convection and local convergence between opposing sea breeze systems and inland convection. Linear correlation distance-decay trends are determined for the island and for the predetermined rainfall affinity areas, including estimates of distance-decay rates along field axes. A suite of exponential cone functions is derived to model the empirically observed fields, representing their typical bidimensional structure. Use of the fitted cones to model a typical rainfall day, where the rainfall distribution mimics quite closely the longer term average, yields quite good results and suggests that their use may be more appropriate than applying more conventional general interpolation techniques in the mapping of spatial rainfall distributions.With 14 Figures     

Geology of the volcanic-hosted Brockman rare-metals deposit,Halls Creek Mobile Zone,northwest Australia. l. Volcanic environment,geochronology and petrography of the Brockman volcanics

Summary Rare-metals mineralization of the Brockman deposit (Halls Creek Mobile Zone, NW Australia) is hosted in a fluorite-bearing, rhyolitic volcaniclastic unit informally termed the Niobium Tuff. The Tuff, more correctly described as a tuffaceous volcaniclastic deposit, is the lowermost unit of a sequence of trachyte-to-rhyolite lavas, trachyandesite subvolcanic rocks, and volcaniclastic units of the Brockman volcanics located within the Halls Creek Group, a thick, early Proterozoic volcano-sedimentary sequence. High precision SHRIMP ion-microprobe zircon dating of the Niobium Tuff gives an eruption age of 1870 ± 4 Ma. Regional geochronological constraints indicate the Niobium Tuff was deposited about 15 Ma before major orogenic activity affected the area. Despite folding, faulting and low-grade metamorphism, the Brockman volcanics show excellent preservation of primary volcanological features, including pillow-lavas and pillow-breccias, that suggest a dominantly subaqueous, below-wave-base emplacement environment. The style of eruption products and magma volume constraints suggest the trachyte-dominated volcanics were erupted from a small shield volcanic complex probably in a rift-related basin in a shallow-marine setting. The tectonic setting was intraplate but differs from most Tertiary to recent trachyte volcanic complexes which are largely subaerial, are built on relatively thick continental crust, and show no post-eruptive orogenic history. Brockman-style rare-metal deposits are characterized by preservation of subaqueous volcanics beneath a thick sedimentary sequence, eruption of early incompatible-element enriched products followed by less differentiated magmas, and fine-grained mineralogy influenced by alteration processes. Prospects exist for discovery of analogous deposits, particularly in early Proterozoic mobile belts and Tertiary intraplate shield volcanic provinces.

---

Geologie der an Vulkanite gebundenen Seltene-Metalle-Lagerstätte Brockman, Halls Creek Mobile Zone, Nordwest-Australien. I. Vulkanologie, Geochronologie und Petrographie der Brockman-Vulkanite

Zusammenfassung Selten-Metall-Vererzung der Brockman-Lagerstätten (Halls Creek Mobile Zone, NW Australien) sitzt in einer Fluorit-führenden, rhyolitischen vulkanoklastischen Einheit auf, die informell als der Niob-Tuff bezeichnet wird. Es handelt sich hier um ein vulkanoklastisches Tuff-Sediment; dieses ist die unterste Einheit einer Abfolge vor Trachyt bis Rhyolit-Laven, trachyandesitischen Subvulkaniten and vulkanoklastischen Einheiten der Brockman Vulkanite innerhalb der Halls Creek Gruppe, einer mächtigen frühproterozoischen vulkanosedimentdren Abfolge. Prüzisionsdatierungen des Niob Tuffs mit der SHRIMP Ionen-Mikrosonde ergeben ein Eruptionsalter von 1870 +-4 Millionen Jahren. Regionale geochronologische Zusammenhänge zeigen, daß der Niob Tuff etwa 15 Millionen Jahre vor einer größeren Orogenese, die das Gebiet betroffen hat, abgelagert wurde. Trotz Faltung, Bruchtektonik and niedriggradiger Metamorphose zeigen die Brockman-Vulkanite einen hervorragenden Erhaltungszustand primärer vulkanologischer Erscheinungen. Diese urnfassen auch Kissenlaven und KissenBreckzien, die eine vorwiegend subaquatische Ablagerung in ruhigem Wasser erkennen lassen. Die Art der Eruptionsprodukte and das Magmavolumen zeigen, daß die Trachytdominierten Vulkanite von einem kleinen Schildvulkan stammer, wahrscheinlich in einem Becken in einer Rift-Situation im seichten marinen Milieu. Die tektonische Situation war intraplate, aber unterscheidet sich von den meisten tertiären bis rezenten trachytischen Vulkan-Komplexen, die hauptsächlich subaerisch sind, auf einer relativ mächtigen kontinentalen Kruste aufsitzen, and keine post-eruptive orogene Entwicklung zeigen. Seltene-Metalle-Lagerstatten des Brockman-Typs rind durch die Erhaltung subaquatischer-Vulkanite unterhalb einer machtigen sedimentdren Abfolge gekennzeichnet; welters durch frühe Eruptionsprodukte, die an inkompatiblen Elementen angereichert sind, auf die dann weniger differenzierte Magmen folgten, und schließlich durch einen feinkörnigen Mineralbestand, der vor Umwandlungsprozessen betroffen war. Es besteht die Möglichkeit der Entdeckung analoger Lagerstätten, besonders in frühproterozoischen mobilen Gürteln und in tertidren intraplate Schildvulkan Provinzen.

---

---

With 5 Figures

---

Deceased     

Efficient generation of conditional simulations by chebyshev matrix polynomial approximations to the symmetric square root of the covariance matrix

Consider the problem of generating a realization y₁ of a Gaussian random field on a dense grid of points ₁ conditioned on field observations y₂ collected on a sparse grid of points ₂. An approach to this is to generate first an unconditional realization y over the grid =_{1 2}, and then to produce y₁ by conditioning y on the data y₂. As standard methods for generating y, such as the turning bands, spectral or Cholesky approaches can have various limitations, it has been proposed by M. W. Davis to generate realizations from a matrix polynomial approximations to the square root of the covariance matrix. In this paper we describe how to generate a direct approximation to the conditional realization y₁, on ₁ using a variant of Davis' approach based on approximation by Chebyshev polynomials. The resulting algorithm is simple to implement, numerically stable, and bounds on the approximation error are readily available. Furthermore we show that the conditional realization y₁ can be generated directly with a lower order polynomial than the unconditional realization y, and that further reductions can be achieved by exploiting a nugget effect if one is present. A pseudocode version of the algorithm is provided that can be implemented using the fast Fourier transform if the field is stationary and the grid ₁ is rectangular. Finally, numerical illustrations are given of the algorithm's performance in generating various 2-D realizations of conditional processes on large sampling grids.