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     

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.     

Cycle anatomy and variability in the storm-dominated type cincinnatian (Upper Ordovician): coming to grips with cycle delineation and genesis

Although parasequence and sequence are scale-independent terms, they are frequently applied only to specific scales of cycles. For example, meter-scale cycles are commonly assumed to be parasequences or PACs. In the Upper Ordovician Kope and Fairview Formations of northern Kentucky, we examined a succession of 50 meter-scale cycles that have been variously interpreted as deepening-upward, shallowing-upward, or showing no relationship with water depth. Our analysis shows that these cycles, characterized by shifts in storm-bed proximality, are highly variable in their thickness and internal construction. Most cycles are best considered high-frequency sequences, because deepening-upward intervals are common, and many cycles contain evidence of abrupt basinward shifts in facies as expected at sequence boundaries. A minority fit the parasequence model of shallowing-upward cycles bounded by flooding surfaces. Larger, 20 m scale cycles are defined by systematic thickening and thinning trends of meter-scale cycles. However, meter-scale cycles do not display any systematic trends in cycle anatomy as a function of position within the 20 m cycles or position within the Kope and Fairview Formations. The high cycle variability and the lack of systematic stratigraphic organization with respect to longer-term cyclicity reflect either the irregularity of relative sea-level changes, the poor recording of sea-level changes in this deep-water setting, or the generation of these cycles by climate-induced cyclicity in storm intensity. These three mechanisms would generate similar patterns at the outcrop scale, so it is not possible at the present to distinguish between them.     

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.     

Carbon isotope systematics of individual hydrocarbons in hydrothermal petroleum from Escanaba Trough, northeastern Pacific Ocean

We submitted individual aliphatic and polycyclic aromatic hydrocarbons in samples of hydrothermal petroleum from Escanaba Trough to compound specific isotope analysis to trace their origins. The carbon isotope compositions of the alkanes and polycyclic aromatic hydrocarbons (means -27.5 and -24.7%, respectively) reflect a primarily terrestrial organic matter source.     

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     

Non-spherical source-surface model of the heliosphere: a scalar formulation

The source-surface method offers an alternative to full MHD simulation of the heliosphere. It entails specification of a surface from which the solar wind flows normally outward along straight lines. Compatibility with MHD results requires this (source) surface to be non-spherical in general and prolate (aligned with the solar dipole axis) in prototypical axisymmetric cases. Mid-latitude features on the source surface thus map to significantly lower latitudes in the heliosphere. The model is usually implemented by deriving the B field (in the region surrounded by the source surface) from a scalar potential formally expanded in spherical harmonics, with coefficients chosen so as to minimize the mean-square tangential component of B over this surface. In the simplified (scalar) version the quantity minimized is instead the variance of the scalar potential over the source surface. The scalar formulation greatly reduces the time required to compute required matrix elements, while imposing essentially the same physical boundary condition as the vector formulation (viz., that the coronal magnetic field be, as nearly as possible, normal to the source surface for continuity with the heliosphere). The source surface proposed for actual application is a surface of constant , where r is the heliocentric distance and B is the scalar magnitude of the B field produced by currents inside the Sun. Comparison with MHD simulations suggests that k 1.4 is a good choice for the adjustable exponent. This value has been shown to map the neutral line on the source surface during Carrington Rotation 1869 (May–June 1993) to a range of latitudes that would have just grazed the position of Ulysses during that month in which sector structure disappeared from Ulysses magnetometer observations.     

Internal migration patterns for U.S. foreign-born, 1985-1990

"Using 1990 census data, this paper calculates the flow [of the foreign-born population in the United States] between states and from abroad; also, a multinomial logistic regression model of destination choice is estimated for resettlers and for migrants from abroad. There are three major findings. Firstly, Florida and California are the largest recipients of foreign-born resettlers; New York is the biggest loser of secondary migrants. Secondly, the presence of large communities of Mexicans in California and Cubans in Florida are very attractive to both resettlers and migrants from abroad. Thirdly, immigrants that are most in need of ethnic support networks (i.e. new arrivals or immigrants with low human capital resources) are most likely to choose one of the traditional immigration states as a destination."