High-resolution sequence stratigraphy of the Lower Triassic ‘Buntsandstein’ in the Netherlands and northwestern Germany

Detailed log correlations of the largely fluvio-lacustrine Lower Triassic Buntsandstein (Late Permian-Early Anisian), carried out on 80 wells in the Dutch onshore and offshore areas, can be linked to northwest-German high-resolution sequence stratigraphy. The correlations show that cyclic sedimentation occurred in large parts of the basin. Seven 1st-order sequences are recognised, namely the Main Claystone, Rogenstein, Volpriehausen, Detfurth, Hardegsen, Solling and Lower Röt Sequences. They are overlain by the lower part of the Upper Röt–Lower Muschelkalk Sequence. Distinct sequence boundaries have been identified at the bases of four sequences: Volpriehausen, Detfurth, Solling and Upper Röt. The higher-order sequences consist of fining-upwards cycles with a thickness of up to tens of metres. The sequences are laterally persistent and have a characteristic expression on gamma-ray and sonic logs. In the Lower Buntsandstein, they display a uniform character throughout most of the area, with only minor differences in thickness or lithology. NNE-oriented lows and swells were formed during deposition of the Volpriehausen, Detfurth and Hardegsen Sequences. Uplift prior to the deposition of the Solling Sequence caused deep erosion on the swells in the basin and minor erosion in the lows. The high-resolution sequences probably represent alternating, relatively wet and dry climatic periods, with a periodicity of about 100 000 years. An analysis of the sequences suggests that their reduced thickness on the swells is mainly the effect of erosion. This is supported by analyses of the accumulation patterns and rates.     

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.     

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

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     

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.     

Textures in a poikilitic olivine gabbro cumulate: Evidence for supercooling

Summary Supercooling of repeated influxes of high-MgO ( 13.5 wt% MgO) tholeiitic magma emplaced along the floor of the Kap Edvard Holm magma chamber led to the formation of fine grained, intraplutonic quench zones that define the bases of 12 macrorhythmic units. The quench zones grade upwards into poikilitic cumulates containing olivine dendrites, suggesting supercooling by 40°C. Supercooling led to early, precumulus crystallization of metastable, non-equilibrium plagioclase, followed by co-precipitation of cumulus plagioclase and olivine. Later nucleation of clinopyroxene followed by rapid growth led to crystallization of zoned oikocrysts during the cumulus stage, followed by postcumulus overgrowth and crystallization of discrete intercumulus grains. At least four generations of plagioclase crystals appear to be preserved.

---

Gefüge in einem poikilitischen Olivingabbro-Kumulat: Hinweis auf Unterkühlung

Zusammenfassung Unterkühlung von mehrmals nachströmendem MgO-reichen (ca. 13,5 Gew.% MgO), tholeiitischen Magma, das sich am Boden der Magmenkammer vom Kap Edvard Holm sammelte, führte zur Bildung von feinkbrnigen Zonen mit Abschreckungsgefüge, die jeweils die Basis von 12 makrorhythmischen Einheiten bilden. Diese Zonen gehen im Hangenden in poikilitische Kumulate über, die Olivin-Dendriten enthalten, was auf eine Unterkühlung um ca. 40°C hindeutet. Die Unterkühlung führte zu einer frühen präkumulativen Kristallisation von metastabilem, nicht im Gleichgewicht befindlichem Plagioklas, gefolgt von gleichzeitiger Präzipitation von Kumulus-Plagioklas und -Olivin. Eine spätere Keimbildung von Klinopyroxen, die von raschem Wachstum abgelöst wurde, führte während der Akkumulation zur Kristallisation von zonierten Oikokrysten. Danach folgte eine postkumulative Überwachsung und die Kristallisation von Einzelkörnern in den Zwischenräumen. Zumindest vier Generationen von Plagioklaskristallen scheinen erhalten zu sein.

---

---

With 5 Figures     

Sector collapse,sedimentation and clast population evolution at an active island-arc volcano: Stromboli,Italy

Sciara del Fuoco is the subaerial part of a partially filled sector-collapse scar that extends to 700 m below sea level on Stromboli volcano. The collapse occurred <5000 years ago, involved 1.81 km³ of rock and is the latest of a series of major collapses on the north-west flank of Stromboli. A north-east trending arc-axial fault system channels magmas into the volcano and has caused tilting and/or downthrow to the north-west. The slope of the partial cone constructed between the lateral walls of the collapse scar acts as a channelway to the sea for most eruptive products. From 700 m below sea level and extending to >2200 m and >10 km from the shore to the NNW, a fan-shaped mounded feature comprises debris avalanche deposits (>4 km³) from two or more sector collapses. Volcaniclastic density currents originating from Sciara del Fuoco follow the topographic margin of the debris avalanche deposits, although overbank currents and other unconfined currents widely cover the mounded feature with turbidites. Historical (recorded) eruptive activity in Sciara del Fuoco is considerably less than that which occurred earlier, and much of the partial fill may have formed from eruptions soon after the sector collapse. It is possible that a mass of eruptive products similar to that in the collapse scar is dispersed as volcanogenic sediment in deep water of the Tyrhennian basin. Evidence that the early post-collapse eruptive discharge was greater than the apparent recent flux (2kg/s) counters suggestions that a substantial part of Stromboli's growth has been endogenous. The partial fill of Sciaria del Fuoco is dominated by lava and spatter layers, rather than by the scoria and ash layers classically regarded as main constituents of Strombolian (cinder) cones. Much of the volcanic slope beneath the vents is steeper than the angle of repose of loose tephra, which is therefore rapidly transported to the sea. Delicate pyroclasts that record the magmatic explosivity are selectively destroyed and diluted during sedimentary transport, mainly in avalanches and by shoreline wave reworking, and thus the submarine deposits do not record well the extent and diversity of explosive activity and associated clast-forming processes. Considerable amounts of dense (non-vesicular) fine sand and silt grains are produced by breakage and rounding of fragments of lava and agglutinate. The submarine extension of the collapse scar, and the continuing topographic depression to >2200 m below sea level, are zones of considerable by-passing of fine sand and silt, which are transported in turbidity currents. Evidently, volcanogenic sediments dispersed around island volcanoes by density currents are unlikely to record well the true spectrum and relative importance of clast-forming processes that occurred during an eruption. Marine sedimentary evidence of magmatic explosivity is particularly susceptible to partial or complete obliteration, unless there is a high rate of discharge of pyroclastic material into the sea.