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.     

The Ms = 6.2, June 15, 1995 Aigion earthquake (Greece): evidence for low angle normal faulting in the Corinth rift

We present the results of a multidisciplinary study of the M_s = 6.2, 1995, June 15, Aigion earthquake (Gulf of Corinth, Greece). In order to constrain the rupture geometry, we used all available data from seismology (local, regional and teleseismic records of the mainshock and of aftershocks), geodesy (GPS and SAR interferometry), and tectonics. Part of these data were obtained during a postseismic field study consisting of the surveying of 24 GPS points, the temporary installation of 20 digital seismometers, and a detailed field investigation for surface fault break. The Aigion fault was the only fault onland which showed detectable breaks (< 4 cm). We relocated the mainshock hypocenter at 10 km in depth, 38 ° 21.7 N, 22 ° 12.0 E, about 15 km NNE to the damaged city of Aigion. The modeling of teleseismic P and SH waves provides a seismic moment M_o = 3.4 10¹⁸ N.m, a well constrained focal mechanism (strike 277 °, dip 33 °, rake – 77°), at a centroidal depth of 7.2 km, consistent with the NEIC and the revised Harvard determinations. It thus involved almost pure normal faulting in agreement with the tectonics of the Gulf. The horizontal GPS displacements corrected for the opening of the gulf (1.5 cm/year) show a well-resolved 7 cm northward motion above the hypocenter, which eliminates the possibility of a steep, south-dipping fault plane. Fitting the S-wave polarization at SERG, 10 km from the epicenter, with a 33° northward dipping plane implies a hypocentral depth greater than 10 km. The north dipping fault plane provides a poor fit to the GPS data at the southern points when a homogeneous elastic half-space is considered: the best fit geodetic model is obtained for a fault shallower by 2 km, assuming the same dip. We show with a two-dimensional model that this depth difference is probably due to the distorting effect of the shallow, low-rigidity sediments of the gulf and of its edges. The best-fit fault model, with dimensions 9 km E–W and 15 km along dip, and a 0.87 m uniform slip, fits InSAR data covering the time of the earthquake. The fault is located about 10 km east-northeast to the Aigion fault, whose surface breaks thus appears as secondary features. The rupture lasted 4 to 5 s, propagating southward and upward on a fault probably outcropping offshore, near the southern edge of the gulf. In the shallowest 4 km, the slip – if any – has not exceeded about 30 cm. This geometry implies a large directivity effect in Aigion, in agreement with the accelerogram aig which shows a short duration (2 s) and a large amplitude (0.5 g) of the direct S acceleration. This unusual low-angle normal faulting may have been favoured by a low-friction, high pore pressure fault zone, or by a rotation of the stress directions due to the possible dip towards the south of the brittle-ductile transition zone. This fault cannot be responsible for the long term topography of the rift, which is controlled by larger normal faults with larger dip angles, implying either a seldom, or a more recently started activity of such low angle faults in the central part of the rift.     

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

A multifractal analysis of lidar measured water vapour

Accurate and fast-response measurements of space-time observations of specific humidity were made above a drying land surface at the University of California at Davis, using the Los Alamos water Raman-lidar. In an attempt to quantify the space-time intermittency features of turbulent flows in the lower atmosphere, a multifractal analysis of these water vapour measurements was performed. The structure of the specific humidity, (x, t), was analyzed quantifying a scalar gradient measure both in time and space, for all possible one-dimensional cuts, i.e. and . The results confirm the multifractal nature of this scalar gradient measure (a type of scalar dissipation rate) and show that humidity measurements at fixed times (x) are more intermittent (e.g. have less entropy dimension) than those at fixed locations in space (t). Similar multifractal behaviour of the spatial data, with and without a transformation from the observed wind velocities, supports the validity of Taylor's hypothesis for the studied fields.     

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     

Chemical composition of hard- and softrock groundwaters from central Norway with special consideration of fluoride and Norwegian drinking water limits

Groundwaters from crystalline and metamorphic rocks (hardrocks) and from Quaternary deposits, i.e., alluvial and glacial deposits (softrocks) from the counties of Nord-Trøndelag and Sør-Trøndelag were analyzed for major and minor elements and ions including fluoride. The median concentration of F^– in water from the hardrock aquifers is 0.28 mg/l (14.7 eq/l) in contrast to water from softrock aquifers in which it is found to be 0.05 mg/1 (2.6 eq/l). More importantly, ca. 15% of the locations where water was abstracted from hardrock wells contain 1.5 mg/l (78.9 eq/l) F^– or more. Thus, 15% of all hardrock wells returned F^– results that are at or above the maximum recommended value for drinking water. Of the softrock wells, none are above 1 mg/l. Geologists would normally expect higher F^–contents in groundwaters derived from acid rocks, e.g., in granitic or gneissic areas. When comparing the host lithology with the observed F^–contents, however, no clear relationship between F^– content and lithology is visible. The highest observed F^– values actually occur in gneissic host rocks. However, wells drilled in amphibolites/greenstones, mica schists, calcareous rocks, and sedimentary rocks all returned some analytical results above 1.5 mg/l F^–. These results suggest that all hardrock wells drilled should be tested for F^– and the users informed about the results and advised to take any necessary precautions. When applying the recently proposed Norwegian drinking water limits to our data, 51% of all softrock well waters and 56% of all hardrock well waters are unfit for consumption without prior treatment, although we analyzed only for about half of the proposed elements/parameters. This result seriously questions the concept of fixed action levels—many of them with totally unproven health implications—for so many parameters/elements for hardrock groundwaters.     

Towards an understanding of joint roughness

Summary It is argued that the currently available joint models are incapable of accurate predictions of joint shear behaviour without resorting to substantial levels of empiricism. This is because these models fail to adequately quantify joint roughness, or appreciate the importance of scale. A novel approach, which uses fractal geometry to investigate joint roughness, is described. This approach goes beyond describing the symptoms of roughness and seeks to find the cause. In the application of this approach, the concepts of fractal geometry, fractal dimension and self similarity are described and used as a framework to formulate a statistically based and practical model for the characterisation of rock joint roughness. Important relationships between the fractal dimension and the more useful statistical parameters of standard deviation of both asperity angles and asperity heights are derived. These relationships not only provide a useful, working method for quantifying joint roughness, but are also shown to provide a basis for understanding the Barton empirical JRC-JCS model. In addition, the fractal model is able to provide conceptual models for the effects of normal stress on the shear behaviour of joints and the scale-dependence of joints.