Effects of defect and pressure on the thermal expansivity of Fe x O

Pressure–volume–temperature measurements have been carried out using synchrotron X-ray diffraction for wüstite at static pressures of 1.9, 2.6, and 5.4 GPa. Our results revealed that the composition change of wüstite and, hence, rearrangements of defect structures are primarily caused by the magnetite (Fe₃O₄) exsolution at temperatures of 523–723 K. Based on the isobaric volume–temperature data collected during cooling, the contribution of compositional variations to the unit-cell volumes of wüstite in the ranges of 300–673 K and 723–1073 K is negligibly small, within the experimental uncertainties. These observations suggest that the measured volume changes in the range of 300–673 K and 723–1,073 K can be attributed to the metal–oxygen bond expansion. Owing to the magnetite exsolution, thermal expansion data are obtained in each experiment at 1.9, 2.6, and 5.4 GPa for wüstite of two different compositions, Fe_0.987O and Fe_0.942O. At all three pressures, Fe_0.942O shows a thermal expansion that is about 30% larger than Fe_0.987O. Such findings represent the first experimental evidence of a substantial effect of nonstoichiometry on thermal expansivity, and based on previous thermodynamic calculations of the defect formation and interaction, this effect is likely associated with the distinct defects arrangements in iron-rich and more iron-deficient wüstite. This study also presents thermal equations of state for wüstite of two different compositions. Such volume-related properties at high temperatures are experimentally difficult to obtain in wüstite but important for thermodynamic studies in the binary Fe–O system.     

Mineral exploration under cover: Characterizing mineralizing fluid systems

Mineral exploration in areas of little bedrock exposure must increasingly rely on the predictive quality of geological, geochemical and geophysical data. In this contribution analysis of fluid inclusions is used to characterize fluid systems associated with alteration and mineralization in two locations in South Australia. In the Curnamona Province, fluids related to albitization and calc-silicate alteration with associated brecciation evolve from a saline sodic composition to a metalliferous, saline fluid also found associated with sulphide mineralization. The microthermometric distinction of inclusion types is evidence for pressure fluctuation as a cause for breccia formation. In the Central Gawler Craton gold province fluid inclusion studies have characterized the compositional and constrained the physical conditions of mineralizing fluids. Regionally identified low salinity CO₂-bearing fluids at temperatures < 350 °C and pressures < 1.5 kbar indicate a shallow-level orogenic [Groves, D.I., Goldfarb, R.J., Gebre-Mariam, M., Hagemann, S.G., Robert, F., 1998. Orogenic gold Deposits: A proposed classification in the context of their crustal distribution and relationship to other gold deposit types. Ore Geology Reviews, 13, 7–27.] or intrusion-related [Lang, J.R., Baker, T., 2001. Intrusion related gold systems: the present level of understanding. Mineralium Deposita, 36, 477–489.] setting for the mineralization.     

Permo–Triassic intraplate magmatism and rifting in Eurasia: implications for mantle plumes and mantle dynamics

At the transition from the Permian to the Triassic, Eurasia was the site of voluminous flood-basalt extrusion and rifting. Major flood-basalt provinces occur in the Tunguska, Taymyr, Kuznetsk, Verkhoyansk–Vilyuy and Pechora areas, as well as in the South Chinese Emeishen area. Contemporaneous rift systems developed in the West Siberian, South Kara Sea and Pyasina–Khatanga areas, on the Scythian platform and in the West European and Arctic–North Atlantic domain. At the Permo–Triassic transition, major extensional stresses affected apparently Eurasia, and possibly also Pangea, as evidenced by the development of new rift systems. Contemporaneous flood-basalt activity, inducing a global environmental crisis, is interpreted as related to the impingement of major mantle plumes on the base of the Eurasian lithosphere. Moreover, the Permo–Triassic transition coincided with a period of regional uplift and erosion and a low-stand in sea level. Permo–Triassic rifting and mantle plume activity occurred together with a major reorganization of plate boundaries and plate kinematics that marked the transition from the assembly of Pangea to its break-up. This plate reorganization was possibly associated with a reorganization of the global mantle convection system. On the base of the geological record, we recognize short-lived and long-lived plumes with a duration of magmatic activity of some 10–20 million years and 100–150 million years, respectively. The Permo–Triassic Siberian and Emeishan flood-basalt provinces are good examples of “short-lived” plumes, which contrast with such “long lived” plumes as those of Iceland and Hawaii. The global record indicates that mantle plume activity occurred episodically. Purely empirical considerations indicate that times of major mantle plume activity are associated with periods of global mantle convection reorganization during which thermally driven mantle convection is not fully able to facilitate the necessary heat transfer from the core of the Earth to its surface. In this respect, we distinguish between two geodynamically different scenarios for major plume activity. The major Permo–Triassic plume event followed the assembly Pangea and the detachment of deep-seated subduction slabs from the lithosphere. The Early–Middle Cretaceous major plume event, as well as the terminal–Cretaceous–Paleocene plume event, followed a sharp acceleration of global sea-floor spreading rates and the insertion of new subduction zone slabs deep into the mantle. We conclude that global plate kinematics, driven by mantle convection, have a bearing on the development of major mantle plumes and, to a degree, also on the pattern of related flood-basalt magmatism.     

Geobarometry of low-temperature eclogites: applications of isothermal pressure-activity calculations

Estimation of metamorphic pressures in low temperature eclogite (Type C) is difficult because of the high variance mineral assemblages and problems in geothermometry, solution properties of low-temperature omphacite, and the thermodynamic properties of clinozoisite. We have considered equilibria in the CaO–FeO–MgO–TiO₂–Al₂O₃–SiO₂–H₂O (CFMTASH) system involving the phase components, quartz, rutile, kyanite, ilmenite, almandine, pyrope, grossular, clinozoisite, sphene, diopside, and H₂O-fluid There are four linearly independent equilibria involving the phase components in this system. Because kyanite can crystallize as a nearly pure phase, the lack of kyanite in a rock indicates that a _Al2SiO5 is<1.0. If we can estimate temperature independently, we can solve for a _Al2SiO5 and pressure by using two of the equilibria in isothermal pressure-activity diagrams. We have applied this approach to eclogites from New Caledonia and from southwestern Oregon. For the New Caledonia eclogites, calculated pressures range from 11.2 to 13.6 kbar at 500°C, and are consistent with the minimum pressures based upon the presence of jadeitic pyroxene+quartz and the lack of stable albite. Oregon eclogites come from different tectonic blocks and calculated minimum pressures of 11–12 kbar are based upon the presence of jadeitic pyroxene+rutile+garnet and lack of stable albite and ilmenite at reduced values of a _SiO2 (0.7–0.9).     

Melt extraction pathways and stagnation depths beneath the Madeira and Desertas rift zones (NE Atlantic) inferred from barometric studies

The Madeira and Desertas Islands (eastern North Atlantic) show well-developed rift zones which intersect near the eastern tip of Madeira (São Lourenço peninsula). We applied fluid inclusion barometry and clinopyroxene-melt thermobarometry to reconstruct levels of magma stagnation beneath the two adjacent rifts and to examine a possible genetic relationship during their evolution. Densities of CO₂-dominated fluid inclusions in basanitic to basaltic samples from São Lourenço yielded frequency maxima at pressures of 0.57–0.87 GPa (23–29 km depth) and 0.25–0.32 GPa (8–10 km), whereas basanites, basalts and xenoliths from the Desertas indicate 0.3–0.72 GPa (10–24 km) and 0.07–0.12 GPa (2–3 km). Clinopyroxene-melt thermobarometry applied to Ti-augite phenocryst rim and glass/groundmass compositions indicates pressures of 0.45–1.06 GPa (15–35 km; São Lourenço) and 0.53–0.89 GPa (17–28 km; Desertas Islands) which partly overlap with pressures indicated by fluid inclusions. We interpret our data to suggest a multi-stage magma ascent beneath the Madeira Archipelago: main fractionation occurs at multiple levels within the mantle (>15 km depth) and is followed by temporary stagnation within the crust prior to eruption. Depths of crustal magma stagnation beneath São Lourenço and the Desertas differ significantly, and there is no evidence for a common shallow magma reservoir feeding both rift arms. We discuss two models to explain the relations between the two adjacent rift systems: Madeira and the Desertas may represent either a two-armed rift system or two volcanic centres with separate magma supply systems. For petrological and volcanological reasons, we favour the second model and suggest that Madeira and the Desertas root in distinct regions of melt extraction. Magma focusing into the Desertas system off the hotspot axis may result from lithospheric bending caused by the load of the Madeira and Porto Santo shields, combined with regional variations in melt production due to an irregularly shaped plume.Electronic Supplementary Material Supplementary material is available for this article if you access the article at . A link in the frame on the left on that page takes you directly to the supplementary material.Editorial responsibility: J. Hoefs     

Shifting Economic Impacts from Weather Extremes in the United States: A Result of Societal Changes,Not Global Warming

Loss values from extremes in the U.S. and elsewhere have been more qualitativethan quantitative, but recent pressures for better information have led to newassessments and better estimates of financial losses from extremes. These pressureshave included concerns over potential impacts of more extremes due to global warmingfostered by ever increasing costs to the insurance industry and government from weather extremes; plus a series of massive losses during the past 15 years (drought of 1988–1989,Hurricane Andrew in 1992, and Midwestern 1993 floods). These recent assessmentsattempted to adjust data for societal changes over time and thus derived new and betterestimates of losses for seven major extremes than existed previously. Three extremeshave annual average losses in excess of a billion dollars (1998 dollars) includinghurricanes ($4.2 billion), floods ($3.2 billion), and severe local storms ($1.6 billion).One extreme and its adjusted losses exhibit upward trends (floods), but all others showno increases with time or temporal decreases (hail, hurricanes, tornadoes, and severethunderstorms). Annual national losses during 1950–1997 from the three major extremes, plus four others (hail, tornadoes, winter storms, and wind storms), collectively reveal no upward or downward trend over time, with an average annual loss of $10.3 billion. The quality loss values do not indicate an increase as has been postulated for global warming. The good news is that better estimates of impacts now exist, but the bad news is that they are still estimates and do not include sizable unmeasured losses. If accurate data on the economic impacts from weather extremes are seen as important for scientific research and policy-making for global warming, the U.S. needs a continuing program to adequately measure losses from weather extremes.     

Genetic significance of fluid inclusions in the CSA Cu–Pb–Zn deposit, Cobar, Australia

CSA mine exploits a ‘Cobar-type’ Cu–Pb–Zn±Au±Ag deposit within a cleaved and metamorphosed portion of the Cobar Supergroup, central New South Wales. The deposit comprises systems of ‘lenses’ that encompass veins, disseminations and semi-massive to massive Cu–Pb–Zn ores. The systems and contained lenses truncate bedding, are approximately coplanar with regional cleavage and similarly oriented shear zones and plunge parallel to the elongation lineation. Systems have extreme vertical continuity (>1000 m), short strike length (400 m) and narrow width (100 m), exhibit vertical and lateral ore-type variation and have alteration haloes. Models of ore formation include classical hydrothermalism, structurally controlled remobilisation and polymodal concepts; syntectonic emplacement now holds sway.Fluid inclusions were examined from quartz±sulphide veins adjacent to now-extracted ore, from coexisting quartz–sulphide within ore, and from vughs in barren quartz veins. Lack of early primary inclusions precluded direct determination of fluids associated with D₂–D₃ ore and vein emplacement. Similarly, decrepitation (by near-isobaric heating) of the two oldest secondary populations precluded direct determination of fluid phases immediately following D₂–D₃ ore and vein emplacement. Post-decrepitation outflow (late D₃ to early post-D₃) is recorded by monophase CH₄ inclusions. Entrained outflow of deeply circulated meteoric fluid modified the CH₄ system; modification is recorded by H₂O+CH₄ and H₂O+(trace CH₄) secondary populations and by an H₂O+(trace CH₄) primary population. The contractional tectonics (D₂–D₃) of ore emplacement was superseded by relaxational tectonics (D_4P) that facilitated meteoric water penetration and return flow.Under D₂ prograde metamorphism, entrapment temperatures (Tt) and pressures (Pt) for pre-decrepitation secondary inclusions are estimated as Tt300–330 °C and Pt1.5–2 kbar≈Plith (the lithostatic pressure). Decrepitation accompanied peak metamorphism (T350–380 °C) in mid- to late-D₃, while in late-D₃ to early post-D₃, essentially monophase CH₄ secondary inclusions were entrapped at Tt350 °C and Pt=1.5–2 kbar≈Plith. Subsequently, abundant CH₄ and entrained meteoric water were entrapped as H₂O+CH₄ secondaries under slowly decreasing temperature (Tt330–350 °C) and constant pressure (Pt1.5–2 kbar). Finally, with increasingly dominant meteoric outflow, H₂O+(trace CH₄) populations record decreasing temperatures (Tt>300 to <350 down to 275–300 °C) at pressures of Phydrostatic<Pt (1 kbar) <Plith (1.5 kbar).The populations of inclusions provide insight into fluid types, flow regimes and P–T conditions during parts of the deposit's evolution. They indirectly support the role of basin-derived CH₄ fluids in ore formation, but provide no insight into a basement-sourced ore-forming fluid. They fully support post-ore involvement of meteoric water. The poorly constrained entrapment history is believed to span 10 Ma from 395 to 385 Ma.     

A review of Antarctic granulite-facies rocks

Precambrian granulite-facies rocks occur in significant proportion in the East Antarctic Precambrian shield. Ages of metamorphic and deformational events range from 2500 m.y. to about 500 m.y., but some rocks are much older, notably the approximately 3500 m.y. ages for crust formation in Enderby Land. Mineral assemblages over most of the area are typical of the hornblende granulite facies, and sparse temperature pressure estimates indicate metamorphism at 700–800°C and 5–8 kbar at reduced water pressures. A terrane of exceptional interest is the Napier complex of Enderby Land, where sapphirine-quartz ± garnet, sillimanite-orthopyroxene, osumilite, and inverted pigeonite are associated with pyroxene-granulite-facies rocks. Metamorphic conditions are estimated to have reached 900°–980°C, 7–9 kbar, and p_H2O < 0.5 kbar. Metamorphism in the Napier complex, and possibly in other parts of East Antarctica, may be associated with large loss of fluid rather than massive influx of CO₂.     

Constraining stress magnitudes using petroleum exploration data in the Cooper–Eromanga Basins, Australia

The magnitude of the in situ stresses in the Cooper–Eromanga Basins have been determined using an extensive petroleum exploration database from over 40 years of drilling. The magnitude of the vertical stress (S_v) was calculated based on density and velocity checkshot data in 24 wells. Upper and lower bound values of the vertical stress magnitude are approximated by S_v = (14.39 × Z)^1.12 and S_v = (11.67 × Z)^1.15 functions respectively (where Z is depth in km and S_v is in MPa). Leak-off test data from the two basins constrain the lower bound estimate for the minimum horizontal stress (S_hmin) magnitude to 15.5 MPa/km. Closure pressures from a large number of minifrac tests indicate considerable scatter in the minimum horizontal stress magnitude, with values approaching the magnitude of the vertical stress in some areas. The magnitude of the maximum horizontal stress (S_Hmax) was constrained by the frictional limits to stress beyond which faulting occurs and by the presence of drilling-induced tensile fractures in some wells. The maximum horizontal stress magnitude can only be loosely constrained regionally using frictional limits, due to the variability of both the minimum horizontal stress and vertical stress estimates. However, the maximum horizontal stress and thus the full stress tensor can be better constrained at individual well locations, as demonstrated in Bulyeroo-1 and Dullingari North-8, where the necessary data (i.e. image logs, minifrac tests and density logs) are available. The stress magnitudes determined indicate a predominantly strike-slip fault stress regime (S_Hmax > S_v > S_hmin) at a depth of between 1 and 3 km in the Cooper–Eromanga Basins. However, some areas of the basin are transitional between strike-slip and reverse fault stress regimes (S_Hmax > S_v ≈ S_hmin). Large differential stresses in the Cooper–Eromanga Basins indicate a high upper crustal strength for the region, consistent with other intraplate regions. We propose that the in situ stress field in the Cooper–Eromanga Basins is a direct result of the complex interaction of tectonic stresses from the convergent plate boundaries surrounding the Indo-Australian plate that are transmitted into the center of the plate through a high-strength upper crust.     

Mechanism of the water and carbon dioxide solubility in oxides and silicates and the role of O−

The dissolution of water does not stop at the OH^– stage but may proceed further towards H₂ plus O^– formation. The discovery of atomic carbon dissolved in minerals suggests that, if CO₂ enters oxides and silicates at high pressures and temperatures, not only [CO₃]^2– ions but also [CO ₄ ^. ]^4– complexes are formed via a charge transfer which produces O^– and essentially zero-valent, atomic carbon. Under P —T-conditions of the mantle, where the solubility for water and CO₂ is high, the silicate phases formed may therefore consist to a large volume fraction of O^– ions which are much smaller than O^2–ions and strongly cova-lently bonding. The implications for the crystal chemistry of high pressure phases, for the petrology of mantle rocks are outlined.     

A new fracture criterion for brittle rocks

A general form of a “fracture function” for isotropic brittle materials is expressed in terms of the three invariants of the stress tensor. The coefficients in the function are determined by use of the small number of experimental data under specific conditions. This function is applicable to an estimate of the fracture condition of brittle rocks under a general stress state i.e., σ₁ ≠ σ₂ ≠ σ₃. The application of this function is attempted for the data of three brittle rocks i.e., Dunham dolomite, Mizuho trachyte, and Westerly granite, reported by previous workers. For the first two, this criterion gives a good estimation of the effect of the intermediate principal stress σ₂ on failure. For the last, the fracture strength at high confining pressure is estimated by use of the several data obtained under very low confining pressures, and the agreement with experimental data is also satisfactory.     

Forecasting rupture dimension using the pattern informatics technique

The Pattern Informatics (PI) technique [Tiampo, K.F., Rundle, J.B., McGinnis, S., Gross, S., Klein, W., 2002. Mean-field threshold systems and phase dynamics: An application to earthquake fault systems, Europhys. Lett., 60, 481–487] is founded on the premise that changes in the seismicity rate are a proxy for changes in the underlying stress. This new approach to the study of seismicity quantifies its local and regional space–time patterns and identifies regions of local quiescence or activation. Here we use a modification of the PI method to quantify localized changes surrounding the epicenters of large earthquakes in California in an attempt to objectively quantify the rupture zones of these upcoming events. We show that this method can be used to forecast the size and magnitude of future earthquakes.     

A periclase-dolomite-calcite carbonatite from the Oka complex,Quebec, and its calculated volatile composition

The eutectic mineral assemblage calcite-dolomite-periclase-apatite-forsterite-magnesioferrite-pyrrhotite-alabandite in a carbonatite dike within the Oka complex, Quebec, buffers the fugacities (and partial pressures) of all gas species in C-O-H-S-F, assuming vapor saturation. At the inferred eutectic (640° C, 1 kbar), the most important gas species and their partial pressures (bars) were: H₂O, 882; CO₂, 110; H₂, 4.6; H₂S, 2.7; CO, 0.5; and CH₄, 0.1. Oxygen fugacity was near the QFM buffer, logf(O₂)=–18.6, and sulfur fugacity was near the QFM-pyrrhotite buffer, logf(S₂)=–5.9. Fluorine fugacity was low, logf(F₂)=–43.9, consistent with the absence of fluoride minerals other than apatite. Presence of a water-rich gas phase is consistent with experiments on synthetic carbonatite systems (e.g. Fanelli et al. 1981), although compositions of the gas phase in published experiments cannot be determined exactly.Contribution no. 390 from the Mineralogical Laboratory, The University of Michigan     

Source parameters of intermediate-depth Vrancea (Romania) earthquakes from empirical Green's functions modeling

Several source parameters (source dimensions, slip, particle velocity, static and dynamic stress drop) are determined for the moderate-size October 27th, 2004 (M_W = 5.8), and the large August 30th, 1986 (M_W = 7.1) and March 4th, 1977 (M_W = 7.4) Vrancea (Romania) intermediate-depth earthquakes. For this purpose, the empirical Green's functions method of Irikura [e.g. Irikura, K. (1983). Semi-Empirical Estimation of Strong Ground Motions during Large Earthquakes. Bull. Dis. Prev. Res. Inst., Kyoto Univ., 33, Part 2, No. 298, 63–104., Irikura, K. (1986). Prediction of strong acceleration motions using empirical Green's function, in Proceedings of the 7th Japan earthquake engineering symposium, 151–156., Irikura, K. (1999). Techniques for the simulation of strong ground motion and deterministic seismic hazard analysis, in Proceedings of the advanced study course seismotectonic and microzonation techniques in earthquake engineering: integrated training in earthquake risk reduction practices, Kefallinia, 453–554.] is used to generate synthetic time series from recordings of smaller events (with 4 ≤ M_W ≤ 5) in order to estimate several parameters characterizing the so-called strong motion generation area, which is defined as an extended area with homogeneous slip and rise time and, for crustal earthquakes, corresponds to an asperity of about 100 bar stress release [Miyake, H., T. Iwata and K. Irikura (2003). Source characterization for broadband ground-motion simulation: Kinematic heterogeneous source model and strong motion generation area. Bull. Seism. Soc. Am., 93, 2531–2545.] The parameters are obtained by acceleration envelope and displacement waveform inversion for the 2004 and 1986 events and MSK intensity pattern inversion for the 1977 event using a genetic algorithm. The strong motion recordings of the analyzed Vrancea earthquakes as well as the MSK intensity pattern of the 1977 earthquake can be well reproduced using relatively small strong motion generation areas, which corresponds to small asperities with high stress drops (300–1200 bar) and high particle velocities (3–5 m/s). These results imply a very efficient high-frequency radiation, which has to be taken into account for strong ground motion prediction, and indicate that the intermediate-depth Vrancea earthquakes are inherently different from crustal events.     

Cenozoic post-collisional brittle tectonic history and stress reorientation in the High Zagros Belt (Iran, Fars Province)

The Zagros fold-and-thrust belt of SW-Iran is among the youngest continental collision zones on Earth. Collision is thought to have occurred in the late Oligocene–early Miocene, followed by continental shortening. The High Zagros Belt (HZB) presents a Neogene imbricate structure that has affected the thick sedimentary cover of the former Arabian continental passive margin. The HZB of interior Fars marks the innermost part of SE-Zagros, trending NW–SE, that is characterised by higher elevation, lack of seismicity, and no evident active crustal shortening with respect to the outer (SW) parts. This study examines the brittle structures that developed during the mountain building process to decipher the history of polyphase deformation and variations in compressive tectonic fields since the onset of collision. Analytic inversion techniques enabled us to determine and separate different brittle tectonic regimes in terms of stress tensors. Various strike–slip, compressional, and tensional stress regimes are thus identified with different stress fields. Brittle tectonic analyses were carried out to reconstruct possible geometrical relationships between different structures and to establish relative chronologies of corresponding stress fields, considering the folding process. Results indicate that in the studied area, the main fold and thrust structure developed in a general compressional stress regime with an average N032° direction of σ₁ stress axis during the Miocene. Strike–slip structures were generated under three successive strike–slip stress regimes with different σ₁ directions in the early Miocene (N053°), late Miocene–early Pliocene (N026°), and post-Pliocene (N002°), evolving from pre-fold to post-fold faulting. Tensional structures also developed as a function of the evolving stress regimes. Our reconstruction of stress fields suggests an anticlockwise reorientation of the horizontal σ₁ axis since the onset of collision and a significant change in vertical stress from σ₃ to σ₂ since the late stage of folding and thrusting. A late right-lateral reactivation was also observed on some pre-existing belt-parallel brittle structures, especially along the reverse fault systems, consistent with the recent N–S plate convergence. However, this feature was not reflected by large structures in the HZB of interior Fars. The results should not be extrapolated to the entire Zagros belt, where the deformation front has propagated from inner to outer zones during the younger events.     

The melting relationships of a madupite from the Leucite Hills,Wyoming, to 30 Kb

The water-undersaturated melting relationships of a mafic, peralkaline, potassic madupite (with about 3% H₂O as shown by chemical analysis) from the Leucite Hills, Wyoming, have been studied at pressures up to 30 kb. At low pressures (<5 kb) leucite is the dominant liquidus phase, but it is replaced at higher pressures by clinopyroxene plus olivine (<5–7 kb), clinopyroxene (7–12.5 kb), clinopyroxene plus minor spinel (12.5–17.5 kb), and clinopyroxene alone (17.5–> 30 kb). At all pressures there is a reaction relationship with falling temperature between melt, olivine and probably clinopyroxene to yield phlogopite. Apatite is stable within the melting interval to pressures above 25 kb. Electron microprobe analyses demonstrate that the clinopyroxene is diopsidic, with low aluminium and titanium contents. Pressure has relatively little effect on the composition of the pyroxene. Phlogopite is also aluminium-poor and has only a moderate titanium content. The experimental results indicate that madupite is not the partial melting product of hydrous lherzolite or garnet lherzolite in the upper mantle and it seems improbable that it is derived by melting of mantle peridotite with a mixed H₂O-CO₂ volatile component. Madupite could, however, be the partial melting product of mica-pyroxenite or mica-olivine-pyroxenite in the upper mantle. It is pointed out that the chemistry of some potassium-rich volcanics may have been affected by volatile transfer and other such processes during eruption and that experimental studies of material affected in this way have little bearing upon the genesis of potassic magmas. Finally, the experimental results enable constraints to be placed upon the P-T conditions of the formation of richterite-bearing mica nodules found in kimberlites and associated rocks. Maximum conditions are 25 kb and 1,100 ° C.     

The partitioning of Fe and Mg between olivine and carbonate and the stability of carbonate under mantle conditions

We have investigated the effect of Fe on the stabilities of carbonate (carb) in lherzolite assemblages by determining the partitioning of Fe and Mg between silicate (olivine; ol) and carbonates (magnesite, dolomite, magnesian calcite) at high pressures and temperatures. Fe enters olivine preferentially relative to magnesite and ordered dolomite, but Fe and Mg partition almost equally between disordered calcic carbonate and olivine. Measurement of K _d (X _Fe ^carb X _Mg ^ol /X _Fe ^ol X _Mg ^carb ) as a function of Fe/ Mg ratio indicates that Fe–Mg carbonates deviate only slightly from ideality. Using the regular solution parameter for olivine W _FeMg ^ol of 3.7±0.8 kJ/mol (Wiser and Wood 1991) we obtain for (FeMg)CO₃ a W _FeMg ^carb of 3.05±1.50 kJ/mol. The effect of Ca–Mg–Fe disordering is to raise K _d substantially enabling us to calculate W _CaMg ^carb -W _CaFe ^carb of 5.3±2.2 kJ/mol. The activity-composition relationships and partitioning data have been used to calculate the effect of Fe/Mg ratio on mantle decarbonation and exchange reactions. We find that carbonate (dolomite and magnesian calcite) is stable to slightly lower pressures (by 1 kbar) in mantle lherzolitic assemblages than in the CaO–MgO–SiO₂(CMS)–CO₂ system. The high pressure breakdown of dolomite + orthopyroxene to magnesite + clinopyroxene is displaced to higher pressures (by 2 kbar) in natural compositions relative to CMS. CO₂. We also find a stability field of magnesian calcite in lherzolite at 15–25 kbar and 750–1000°C.     

A biogenic-silica δO record of climatic change during the last glacial–interglacial transition in southwestern Alaska

Despite growing evidence for environmental oscillations during the last glacial–interglacial transition from high latitude, terrestrial sites of the North Pacific rim, oxygen-isotopic records of these oscillations remain sparse. The lack of data is due partially to the paucity of lakes that contain carbonate sediment suitable for oxygen-isotopic analysis. We report here the first record of oxygen-isotopic composition in diatom silica (δ¹⁸O_Si) from a lake in that region. δ¹⁸O_Si increases gradually from 19.0 to 23.5‰ between 12,340 and 11,000 ¹⁴C yr B.P., reflecting marked climatic warming at the end of the last glaciation. Around 11,000 ¹⁴C yr B.P., δ¹⁸O_Si decreases by 1.7‰, suggesting a temperature decrease of 3.5–8.9 °C at the onset of the Younger Dryas (YD) in southwestern Alaska. Climatic recovery began ca. 10,740 ¹⁴C yr B.P., as inferred from the increase of δ¹⁸O_Si to a maximum of 23.9‰ near the end of the YD. Our data reveal that a YD climatic reversal in southwestern coastal areas of Alaska occurred, but the YD climate did not return to full-glacial conditions.     

Small-scale spatial variation of the stress field in the back-arc Aegean area: Results from the seismotectonic study of the broader area of Mygdonia basin (N. Greece)

In the present work a detailed seismotectonic study of the broader area of the Mygdonia basin (N. Greece) is performed. Digital data for earthquakes which occurred in the broader Mygdonia basin and were recorded by the permanent telemetric network of the Geophysical Laboratory of the Aristotle University of Thessaloniki during the period 1989–1999 were collected and fault plane solutions for 50 earthquakes which occurred in the study area were calculated with a modified first motions approach which incorporates amplitude and radiation pattern information. Fault plane solutions for the 3 main shocks of Volvi (23/05/78, M_W = 5.8 and 20/06/78, M_W = 6.5) and Arnaia (04/05/95, M_W = 5.8) events and the 1978 aftershock sequence were additionally used. Moreover, data from two local networks established in the Mygdonia basin were also incorporated in the final dataset.Determination of the stress field was realized by the use of the method of Gephart and Forsyth [Gephart, J.W., Forsyth, D.W., 1984. An improved method for determining the regional stress tensor using earthquake focal mechanism data: application to the San Fernando earthquake sequence: Jour. Geophys. Res., v.89, no. B11, p. 9305–9320] for the stress tensor inversion and the results were compared with independent estimates based on the calculation of the average moment tensor [Papazachos, C.B.,Kiratzi, A.A., 1992. A formulation for reliable estimation of active crustal deformation and its application to central Greece. Geophys. J. Int. 111, 424–432]. The obtained stress results show a relatively good agreement between the two approaches, with differences in the azimuth of the dominant extension axis of the order of 10°. Furthermore, comparison with independent information for the mean stress axes provided by the study of kinematics on neotectonic faults [Mountrakis, D., Kilias, A., Tranos, M., Thomaidou, E., Papazachos, C., Karakaisis, G., Scordilis, E., Chatzidimitriou, P., Papadimitriou, E., Vargemezis, G., Aidona, E., Karagianni, E., Vamvakaris, D. Skarlatoudis, A. 2003. Determination of the settings and the seismotectonic behavior of the main seismic-active faults of Northern Greece area using neotectonic and seismological data. Earthquake Planning and Protection Organisation (OASP) (in Greek)] shows a similar agreement with typical misfit of the order 10°. The stress inversion method was modified in order to select one or both nodal planes of the focal mechanism which corresponds to the “true” fault plane of the occurred earthquakes and was able to select a single fault plane in the majority of examined cases. Using this approach, the obtained fault plane rose diagrams are in agreement with results from various neotectonic studies. Moreover, several secondary active fault branches were identified, which are still not clearly observed in the field.     

Fault–valve action and vein development during strike–slip faulting: An example from the Ribeira Shear Zone, Southeastern Brazil

Fluid inclusion microthermometry and structural data are presented for quartz vein systems of a major dextral transcurrent shear zone of Neoproterozoic–Cambrian age in the Ribeira River Valley area, southeastern Brazil. Geometric and microstructural constraints indicate that foliation–parallel and extensional veins were formed during dextral strike–slip faulting. Both vein systems are formed essentially by quartz and lesser contents of sulfides and carbonates, and were crystallized in the presence of CO₂–CH₄ and H₂O–CO₂–CH₄–NaCl immiscible fluids following unmixing from a homogeneous parental fluid. Contrasting fluid entrapment conditions indicate that the two vein systems were formed in different structural levels. Foliation–parallel veins were precipitated beneath the seismogenic zone under pressure fluctuating from moderately sublithostatic to moderately subhydrostatic values (319–397 °C and 47–215 MPa), which is compatible with predicted fluid pressure cycle curves derived from fault–valve action. Growth of extensional veins occurred in shallower structural levels, under pressure fluctuating from near hydrostatic to moderately subhydrostatic values (207–218 °C and 18–74 MPa), which indicate that precipitation occurred within the near surface hydrostatically pressured seismogenic zone. Fluid immiscibility and precipitation of quartz in foliation–parallel veins resulted from fluid pressure drop immediately after earthquake rupture. Fluid immiscibility following a local pressure drop during extensional veining occurred in pre-seismic stages in response to the development of fracture porosity in the dilatant zone. Late stages of fluid circulation within the fault zone are represented dominantly by low to high salinity (0.2 to 44 wt.% equivalent NaCl) H₂O–NaCl–CaCl₂ fluid inclusions trapped in healed fractures mainly in foliation–parallel veins, which also exhibit subordinate H₂O–NaCl–CaCl₂, CO₂–(CH₄) and H₂O–CO₂–(CH₄)–NaCl fluid inclusions trapped under subsolvus conditions in single healed microcracks. Recurrent circulation of aqueous–carbonic fluids and aqueous fluids of highly contrasting salinities during veining and post-veining stages suggests that fluids of different reservoirs were pumped to the ruptured fault zone during faulting episodes. A fluid evolution trending toward CH₄ depletion for CO₂–CH₄–bearing fluids and salinity depletion and dilution (approximation of the system H₂O–NaCl) for aqueous–saline fluids occurred concomitantly with decrease in temperature and pressure related to fluid entrapment in progressively shallower structural levels reflecting the shear zone exhumation history.