Fluid inclusions in upper mantle xenoliths from Northern Patagonia,Argentina: Evidence for an upper mantle diapir

Summary Three generations of fluid inclusions can be recognized in upper mantle xenoliths from alkali basalts of the Somoncura Massif, Northern Patagonia, Argentina. The first (early, primary) one consists of dense CO₂ inclusions which were trapped in the mantle-crust boundary zone (22–36 km minimum trapping depth). Their co-genetic relationship with silicate melt inclusions enables us to constrain their minimum trapping temperature at 1200°C, indicating a high temperature event in a cooler environment. The late (pseudosecondary and secondary) generations of fluid inclusions were classified in accordance with their homogenization temperature to liquid CO₂ (L1) and vapor CO₂ (L2) phase. The minimum trapping depth for the first of the late inclusions (L1) is about 16 km. In spite of the uncertainties related to this value, L1 inclusions indicate that the upper mantle rocks, of which samples were delivered by the basalts, had some residence time in the middle crust where they experienced a metasomatic event. The fact that this event did not destroy the earlier inclusions, places severe constraints on its duration. The second late inclusions (L2) are low-pressure CO₂ inclusions with a minimum trapping depth of only 2 km, presumably a shallow magma chamber of the host basalts. The succession of fluid inclusions strongly points toward a fairly fast uprising upper mantle underneath Northern Patagonia. The petrology and mineral chemistry of the peridotitic xenoliths support this view. Extensive partial melting and loss of these melts is indicated by the preponderance of harzburgites in the upper mantle underneath Northern Patagonia, a fairly unusual feature for a continental upper mantle. That depletion event as well as several metasomatic events — including those which left traces of fluid inclusions — are possibly related to a high-speed diapiric uprise of the upper mantle in this area. The path can be traced from the garnet peridotite stability field into the middle crust, a journey which must have been unusually fast. Differences in rock, mineral, and fluid inclusion properties between geographic locations suggest a diffuse and differential type of diapirism. Future studies will hopefully help to map the full extent and the highs and lows of this diapir and elucidate questions related to its origin and future.

---

Fluid-Einschlüsse in Erdmantel-Xenolithen von Nord-Patagonien: Evidenz für einen Diapir im oberen Erdmantel

Zusammenfassung Erdmantel - Xenolithe in Alkali-Basalten des Somoncure Massivs, Nord-Patagonien, Argentinien, führen drei Generationen von Fluid-Einschlüssen. Die erste (frühe, primäre) Generation besteht aus dichten CO₂-Einschlüssen, welche offenbar in der Mantel-Kruste Grenzzone (22–36 km Minimum-Tiefe) eingeschlossen wurden. CO₂-Einschlüsse sind kogenetisch mit Silikat-Schmelzeinschlüssen. Dies erlaubt die Abschätzung der Einschließ-Temperatur mit minimal 1200°C, was auf ein Hochtemperatur-Ereignis in einer deutlich kühleren Umgebung hinweist. Die späten (pseudosekundäre und sekundäre) CO₂- Fluid-Einschlüsse bilden zwei Generationen von denen die eine in die flüssige (L1), die andere in die Dampfphase (L2) homogenisieren. Die minimale Einschließ-Tiefe für die L1 Generation ist etwa 16 km. Dies bedeutet - auch bei Berücksichtigung der mit diesem Wert verbundenen Ungenauigkeit - daß diese Erdmantel-Gesteine einige Zeit in der mittleren Erdkruste verbrachten und ein metasomatisches Ereignis erlebten, bevor sie von den Basalten zur Erdoberfläche gebracht wurden. Die Tatsache, daß dieses Ereignis die frühen Einschlüsse nicht zerstörte, kann nur bedeuten, daß es von kurzer Dauer war. Die L2-Generation besteht aus Niedrigdruck CO₂-Einschlüssen mit einer Minimum-Einschließtiefe von nur 2 km. Dies könnte in einer seichten Magmakammer des Wirt Basaltes geschehen sein.Die Abfolge von Fluid-Einschlüssen deutet auf einen relativ schnell aufsteigenden oberen Erdmantel unterhalb von Patagonien hin. Die Petrologie und Mineralchemie der peridotitischen Xenolithe unterstützen das. Die Vorherrschaft von Harzburgiten im Erdmantel unterhalb von Nord-Patagonien deutet auf umfangreiche Bildung partieller Schmelzen und deren Abfuhr hin — eine für einen kontinentalen Mantel ungewöhnliche Situation. Sowohl die Verarmungsereignisse, als auch die metasomatischen Veränderungen (einschließlich jene, welche Spuren in Form von Fluid Einschlüssen hinterließen) machen das Vorhandensein eines schnell aufsteigenden Daipirs im oberen Erdmantel dieser Gegend wahrscheinlich. Der Aufstieg kann vom Stabilitätsbereich der Granat-Peridotite bis in die mittlere Kruste verfolgt werden und muß daher relativ schnell erfolgt sein. Unterschiede in Gesteins-, Mineral und Fluid-Eigenschaften zwischen verschiedenen Lokalitäten legen einen diffusen und differenziellen Diapirismus nahe. Zukünftige Studien sollten es ermöglichen, das Gesamtausmaß und die unterschiedlichen Aufstiegshöhen des Diapirs zu kartieren und Hinweise auf seine Entstehung und zukünftige Entwicklung zu erhalten.

---

---

With 5 Figures     

Carbonatites from Eastern Paraguay and genetic relationships with potassic magmatism: C, O, Sr and Nd Isotopes

Summary Geochemical characteristics were systematically determined for Early Cretaceous samples of carbonatitic rocks from Eastern Paraguay (Rio Apa, Amambay and Central Provinces). The data show that all the occurrences have an enriched isotopic signature and that the carbonatites have negligible or absent crustal signature. A petrogenetic model (parent liquids, fractional crystallization, hydrothermal interactions and weathering) is proposed as a function of incompatible trace element, stable (O-C) and radiogenic (Sr-Nd) isotope variations with the aim to test the significance of carbonatitic complexes as a marker of the metasomatized subcontinental lithospheric mantle. The results indicate that the carbonatites and primary carbonates from eastern Paraguay, and those from the north eastern Paraná Basin (SE Brazil), were affected by metasomatic events distinct in time and composition.

---

Karbonatite aus Ost-Paraguay und ihre genetische Beziehung zu Kalium-Magmatismus: C O, Sr und Nd isotope

Zusammenfassung Die geochemischen Charakteristika von frühkretazischen Karbonatitproben aus Ostparaguay (Rio Alpa, Amambay und Zentrale Provinzen) wurden untersucht. Die Daten belegen, daß alle Vorkommen eine isotopische Anreicherungssignatur zeigen und daß ihnen eine entsprechende Krustensignatur fehlt. Ein Petrologisches Modell (Ausgangsschmelze, fraktionierte Kristallisation, hydrothermale Interaktion und Verwitterung) wird auf Grund der Verteilung der inkompatiblen Spurenelemente, der stabilen (C-O) und radiogenen (Sr-Nd) Isotope vorgeschlagen. Es versucht die Bedeutung der Karbonatitkomplexe als Markerhorizonte des metasomatischen subkontinentalen Mantels zu überprüfen. Die Ergebnisse zeigen, daß die Karbonatite und die primären Karbonate in Ostparaguay, und jene aus dem Paraná Becken SüdostBrasiliens durch zeitlich und zusammensetzungsmäßig unterschiedliche metasomatische Prozesse erfaßt wurden.

---

---

With 10 Figures     

Environmental geology problems in the Tyrrhenian coastal area of Santa Marinella, province of Rome, central Italy

 The Tyrrhenian resort of S. Marinella (central Italy) is subjected to significant anthropogenic pressures during the summer vacation period, a common situation all along the Italian coast. Located 65 km NW of Rome on the southern slopes of the Tolfa Mountains, S. Marinella is built on a gently sloping, E–W trending belt which is cut by 14 N–S oriented ephemeral streams that discharge into the Tyrrhenian Sea. The low to medium permeability turbiditic sandstones which outcrop along this belt belong to the Late Cretaceous Pietraforte unit. Three environmental problems are addressed in this study. The first problem is related to the high water supply demand during the summer months which has forced local residents to dig a large number of wells. Extensive pumping from these wells has caused salt-water intrusion into the Pietraforte, thus compromising the domestic use of the groundwater. The second problem consists of the illegal dumping of urban solid waste, material that represents a hazard during significant rain events as well as a possible cause of groundwater contamination. The final issue addressed concerns the flooding potential of the 14 ephemeral streams that cross the inhabited area of S. Marinella, a risk which is highlighted by the disastrous flood which occurred on 2 October 1981 and during the period of the Roman Emperor Settimio Severo (205 A.D.). Some suggestions are proposed to mitigate and contain the effects of these problems. Received: 7 November 1995 / Accepted: 5 December 1996     

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.     

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