THERIA_G: a software program to numerically model prograde garnet growth   总被引：6，自引：4，他引：2  

F. Gaidies  C. de Capitani  R. Abart 《Contributions to Mineralogy and Petrology》2008,155(5):657-671

We present the software program THERIA_G, which allows for numerical simulation of garnet growth in a given volume of rock along any pressure–temperature–time (P–T–t) path. THERIA_G assumes thermodynamic equilibrium between the garnet rim and the rock matrix during growth and accounts for component fractionation associated with garnet formation as well as for intracrystalline diffusion within garnet. In addition, THERIA_G keeps track of changes in the equilibrium phase relations, which occur during garnet growth along the specified P–T–t trajectory. This is accomplished by the combination of two major modules: a Gibbs free energy minimization routine is used to calculate equilibrium phase relations including the volume and composition of successive garnet growth increments as P and T and the effective bulk rock composition change. With the second module intragranular multi-component diffusion is modelled for spherical garnet geometry. THERIA_G allows to simulate the formation of an entire garnet population, the nucleation and growth history of which is specified via the garnet crystal size frequency distribution. Garnet growth simulations with THERIA_G produce compositional profiles for the garnet porphyroblasts of each size class of a population and full information on equilibrium phase assemblages for any point along the specified P–T–t trajectory. The results of garnet growth simulation can be used to infer the P–T–t path of metamorphism from the chemical zoning of garnet porphyroblasts. With a hypothetical example of garnet growth in a pelitic rock we demonstrate that it is essential for the interpretation of the chemical zoning of garnet to account for the combined effects of the thermodynamic conditions of garnet growth, the nucleation history and intracrystalline diffusion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

F. GaidiesEmail:

  相似文献   

An experimental study of the grain-scale processes of peridotite melting: implications for major and trace element distribution during equilibrium and disequilibrium melting   总被引：1，自引：0，他引：1  

Mauro Lo Cascio  Yan Liang  Nobumichi Shimizu  Paul C. Hess 《Contributions to Mineralogy and Petrology》2008,156(1):87-102

The grain-scale processes of peridotite melting were examined at 1,340°C and 1.5 GPa using reaction couples formed by juxtaposing pre-synthesized clinopyroxenite against pre-synthesized orthopyroxenite or harzburgite in graphite and platinum-lined molybdenum capsules. Reaction between the clinopyroxene and orthopyroxene-rich aggregates produces a melt-enriched, orthopyroxene-free, olivine + clinopyroxene reactive boundary layer. Major and trace element abundance in clinopyroxene vary systematically across the reactive boundary layer with compositional trends similar to the published clinopyroxene core-to-rim compositional variations in the bulk lherzolite partial melting studies conducted at similar P–T conditions. The growth of the reactive boundary layer takes place at the expense of the orthopyroxenite or harzburgite and is consistent with grain-scale processes that involve dissolution, precipitation, reprecipitation, and diffusive exchange between the interstitial melt and surrounding crystals. An important consequence of dissolution–reprecipitation during crystal-melt interaction is the dramatic decrease in diffusive reequilibration time between coexisting minerals and melt. This effect is especially important for high charged, slow diffusing cations during peridotite melting and melt-rock reaction. Apparent clinopyroxene-melt partition coefficients for REE, Sr, Y, Ti, and Zr, measured from reprecipitated clinopyroxene and coexisting melt in the reactive boundary layer, approach their equilibrium values reported in the literature. Disequilibrium melting models based on volume diffusion in solid limited mechanism are likely to significantly underestimate the rates at which major and trace elements in residual minerals reequilibrate with their surrounding melt. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.  相似文献   

Preservation of Fe isotope heterogeneities during diagenesis and metamorphism of banded iron formation     

C. D. Frost  F. von Blanckenburg  R. Schoenberg  B. R. Frost  S. M. Swapp 《Contributions to Mineralogy and Petrology》2008,155(1):135-135

  相似文献   

Controls on porphyroblast size along a regional metamorphic field gradient     

David M. Hirsch 《Contributions to Mineralogy and Petrology》2008,155(4):401-415

Garnet-bearing schists from the Waterville Formation of south-central Maine provide an opportunity to examine the factors governing porphyroblast size over a range of metamorphic grade. Three-dimensional sizes and locations for all garnet porphyroblasts were determined for three samples along the metamorphic field gradient spanning lowest garnet through sillimanite grade, using high-resolution X-ray computed tomography. Comparison of crystal size distributions to previous data sets obtained by stereological methods for the same samples reveals significant differences in mode, mean, and shape of the distributions. Quantitative textural analysis shows that the garnets in each rock crystallized in a diffusion-controlled nucleation and growth regime. In contrast to the typical observation of a correlation between porphyroblast size and position along a metamorphic field gradient, porphyroblast size of the lowest-grade specimen is intermediate between the high- and middle-grade specimens’ sizes. Mean porphyroblast size does not correlate with peak temperatures from garnet-biotite Fe-Mg exchange thermometry, nor is post-crystallization annealing (Ostwald Ripening) required to produce the observed textures, as was previously proposed for these rocks. Robust pseudosection calculations fail to reproduce the observed garnet core compositions for two specimens, suggesting that these calc-pelites experienced metasomatism. For each of these two specimens, Monte Carlo calculations suggest potential pre-metasomatism bulk compositions that replicate garnet core compositions. Pseudosection analyses allow the estimation of the critical temperatures for garnet growth: ∼481, ∼477, and ∼485°C for the lowest-garnet-zone, middle-garnet-zone, and sillimanite-zone specimens, respectively. Porphyroblast size appears to be determined in this case by a combination of the heating rate during garnet crystallization, the critical temperature for the garnet-forming reaction and the kinetics of nucleation. Numerical simulations of thermally accelerated, diffusion-controlled nucleation, and growth for the three samples closely match measured crystal size distributions. These observations and simulations suggest that previous hypotheses linking the garnet size primarily to the temperature at the onset of porphyroblast nucleation can only partially explain the observed textures. Also important in determining porphyroblast size are the heating rate and the distribution of favorable nucleation sites.  相似文献   

Mineral inclusions in diamonds from the Kelsey Lake Mine, Colorado, USA: Depleted Archean mantle beneath the Proterozoic Yavapai province     

Daniel J. Schulze  Howard G. Coopersmith  Lori-Ann Pizzolato 《Geochimica et cosmochimica acta》2008,72(6):1685-1695

Thirty-four silicate and oxide inclusions large enough for in situ WDS electron microprobe analysis were exposed by grinding/polishing of 19 diamonds from the Kelsey Lake Mine in the Colorado-Wyoming State Line Kimberlite district. Eighteen olivines, seven Cr-pyropes, four Mg-chromites, and one orthopyroxene in 15 stones belong to the peridotite (P) suite and three garnets and one omphacite in three stones belong to the eclogite (E) suite. The fact that this suite is dominated by the peridotite population is in stark contrast to the other diamond suites studied in the State Line district (Sloan, George Creek), which are overwhelmingly eclogitic. Kelsey Lake olivine inclusions are magnesian (17 of 18 grains in 9 stones are in the range Fo 92.7-93.1), typical of harzburgitic P-suite stones worldwide, but unlike the more Fe-rich (lherzolitic) Sloan olivine suite. Mg-chromites (wt% MgO = 12.8-13.8; wt% Cr₂O₃ = 61.4-66.6) are in the lower MgO range of diamond inclusion chromites worldwide. Seven harzburgitic Cr-pyropes in five stones have moderately low calcium contents (wt% CaO = 3.3-4.3) but are very Cr-rich (wt% Cr₂O₃ = 9.7-16.7). A few stones have been analyzed by SIMS for carbon isotope composition and nitrogen abundance. One peridotitic stone is apparently homogeneous in carbon isotope composition (δ¹³C_PDB = −6.2‰) but with variable nitrogen abundance (1296-2550 ppm). Carbon isotopes in eclogitic stones range from “normal” for the upper mantle (δ¹³C_PDB = −5.5‰) to somewhat low (δ¹³C_PDB = −10.2‰), with little internal variation in individual stones (maximum difference is 3.6‰). Nitrogen contents (2-779 ppm) are lower than in the peridotitic stone, and are lower in cores than in rims. As, worldwide, harzburgite-suite diamonds have been shown to have formed in Archean time, we suggest that the Kelsey Lake diamond population was derived from a block of Archean lithosphere that, at the time of kimberlite eruption, existed beneath the Proterozoic Yavapai province. The mixed diamond inclusion populations from the State Line kimberlites appear to support models in which volumes of Wyoming Craton Archean mantle survive buried beneath Proterozoic continental crust. Such material may be mixed with eclogitic/lherzolitic regimes emplaced beneath or intermingled with the Archean rocks by Proterozoic subduction.  相似文献   

Evidence for crustal degassing of CF₄ and SF₆ in Mojave Desert groundwaters     

Daniel A. Deeds  Martin K. Vollmer  Justin T. Kulongoski  Benjamin R. Miller 《Geochimica et cosmochimica acta》2008,72(4):999-1013

Dissolved tetrafluoromethane (CF₄) and sulfur hexafluoride (SF₆) concentrations were measured in groundwater samples from the Eastern Morongo Basin (EMB) and Mojave River Basin (MRB) located in the southern Mojave Desert, California. Both CF₄ and SF₆ are supersaturated with respect to equilibrium with the preindustrial atmosphere at the recharge temperatures and elevations of the Mojave Desert. These observations provide the first in situ evidence for a flux of CF₄ from the lithosphere. A gradual basin-wide enhancement in dissolved CF₄ and SF₆ concentrations with groundwater age is consistent with release of these gases during weathering of the surrounding granitic alluvium. Dissolved CF₄ and SF₆ concentrations in these groundwaters also contain a deeper crustal component associated with a lithospheric flux entering the EMB and MRB through the underlying basement. The crustal flux of CF₄, but not of SF₆, is enhanced in the vicinity of local active fault systems due to release of crustal fluids during episodic fracture events driven by local tectonic activity. When fluxes of CF₄ and SF₆ into Mojave Desert groundwaters are extrapolated to the global scale they are consistent, within large uncertainties, with the fluxes required to sustain the preindustrial atmospheric abundances of CF₄ and SF₆.  相似文献   

Exurbia from the bottom-up: Confronting empirical challenges to characterizing a complex system   总被引：2，自引：0，他引：2  

Daniel G. Brown  Derek T. Robinson  Li An  Moira Zellner  William Rand  Rick Riolo  Bobbi Low 《Geoforum》2008,39(2):805-818

We describe empirical results from a multi-disciplinary project that support modeling complex processes of land-use and land-cover change in exurban parts of Southeastern Michigan. Based on two different conceptual models, one describing the evolution of urban form as a consequence of residential preferences and the other describing land-cover changes in an exurban township as a consequence of residential preferences, local policies, and a diversity of development types, we describe a variety of empirical data collected to support the mechanisms that we encoded in computational agent-based models. We used multiple methods, including social surveys, remote sensing, and statistical analysis of spatial data, to collect data that could be used to validate the structure of our models, calibrate their specific parameters, and evaluate their output. The data were used to investigate this system in the context of several themes from complexity science, including have (a) macro-level patterns; (b) autonomous decision making entities (i.e., agents); (c) heterogeneity among those entities; (d) social and spatial interactions that operate across multiple scales and (e) nonlinear feedback mechanisms. The results point to the importance of collecting data on agents and their interactions when producing agent-based models, the general validity of our conceptual models, and some changes that we needed to make to these models following data analysis. The calibrated models have been and are being used to evaluate landscape dynamics and the effects of various policy interventions on urban land-cover patterns.  相似文献   

Pressures of Crystallization of Icelandic Magmas   总被引：1，自引：0，他引：1  

Kelley  Daniel F.; Barton  Michael 《Journal of Petrology》2008,49(3):465-492

Iceland lies astride the Mid-Atlantic Ridge and was createdby seafloor spreading that began about 55 Ma. The crust is anomalouslythick (20–40 km), indicating higher melt productivityin the underlying mantle compared with normal ridge segmentsas a result of the presence of a mantle plume or upwelling centeredbeneath the northwestern edge of the Vatnajökull ice sheet.Seismic and volcanic activity is concentrated in 50 km wideneovolcanic or rift zones, which mark the subaerial Mid-AtlanticRidge, and in three flank zones. Geodetic and geophysical studiesprovide evidence for magma chambers located over a range ofdepths (1·5–21 km) in the crust, with shallow magmachambers beneath some volcanic centers (Katla, Grimsvötn,Eyjafjallajökull), and both shallow and deep chambers beneathothers (e.g. Krafla and Askja). We have compiled analyses ofbasalt glass with geochemical characteristics indicating crystallizationof ol–plag–cpx from 28 volcanic centers in the Western,Northern and Eastern rift zones as well as from the SouthernFlank Zone. Pressures of crystallization were calculated forthese glasses, and confirm that Icelandic magmas crystallizeover a wide range of pressures (0·001 to 1 GPa), equivalentto depths of 0–35 km. This range partly reflects crystallizationof melts en route to the surface, probably in dikes and conduits,after they leave intracrustal chambers. We find no evidencefor a shallow chamber beneath Katla, which probably indicatesthat the shallow chamber identified in other studies containssilica-rich magma rather than basalt. There is reasonably goodcorrelation between the depths of deep chambers (> 17 km)and geophysical estimates of Moho depth, indicating that magmaponds at the crust–mantle boundary. Shallow chambers (<7·1 km) are located in the upper crust, and probablyform at a level of neutral buoyancy. There are also discretechambers at intermediate depths (11 km beneath the rift zones),and there is strong evidence for cooling and crystallizing magmabodies or pockets throughout the middle and lower crust thatmight resemble a crystal mush. The results suggest that themiddle and lower crust is relatively hot and porous. It is suggestedthat crustal accretion occurs over a range of depths similarto those in recent models for accretionary processes at mid-oceanridges. The presence of multiple stacked chambers and hot, porouscrust suggests that magma evolution is complex and involvespolybaric crystallization, magma mixing, and assimilation. KEY WORDS: Iceland rift zones; cotectic crystallization; pressure; depth; magma chamber; volcanic glass  相似文献   

Quantitative sinkhole hazard assessment. A case study from the Ebro Valley evaporite alluvial karst (NE Spain)     

Francisco Gutiérrez  Jesús Guerrero  Pedro Lucha 《Natural Hazards》2008,45(2):211-233

Quantitative sinkhole hazard assessments in karst areas allow calculation of the potential sinkhole risk and the performance of cost-benefit analyses. These estimations are of practical interest for planning, engineering, and insurance purposes. The sinkhole hazard assessments should include two components: the probability of occurrence of sinkholes (sinkholes/km² year) and the severity of the sinkholes, which mainly refers to the subsidence mechanisms (progressive passive bending or catastrophic collapse) and the size of the sinkholes at the time of formation; a critical engineering design parameter. This requires the compilation of an exhaustive database on recent sinkholes, including information on the: (1) location, (2) chronology (precise date or age range), (3) size, and (4) subsidence mechanisms and rate. This work presents a hazard assessment from an alluvial evaporite karst area (0.81 km²) located in the periphery of the city of Zaragoza (Ebro River valley, NE Spain). Five sinkholes and four locations with features attributable to karstic subsidence where identified in an initial investigation phase providing a preliminary probability of occurrence of 0.14 sinkholes/km² year (11.34% in annual probability). A trenching program conducted in a subsequent investigation phase allowed us to rule out the four probable sinkholes, reducing the probability of occurrence to 0.079 sinkholes/km² year (6.4% in annual probability). The information on the severity indicates that collapse sinkholes 10–15 m in diameter may occur in the area. A detailed study of the deposits and deformational structures exposed by trenching in one of the sinkholes allowed us to infer a modern collapse sinkhole approximately 12 m in diameter and with a vertical throw of 8 m. This collapse structure is superimposed on a subsidence sinkhole around 80 m across that records at least 1.7 m of synsedimentary subsidence. Trenching, in combination with dating techniques, is proposed as a useful methodology to elucidate the origin of depressions with uncertain diagnosis and to gather practical information with predictive utility about particular sinkholes in alluvial karst settings: precise location, subsidence mechanisms and magnitude, and timing and rate of the subsidence episodes.  相似文献   

Drought and institutional adaptation in the Great Plains of Alberta and Saskatchewan, 1914–1939     

Gregory P. Marchildon  Suren Kulshreshtha  Elaine Wheaton  Dave Sauchyn 《Natural Hazards》2008,45(3):391-411

Agriculture in the southern Great Plains of Canada has been particularly vulnerable to prolonged episodes of drought. Using climate data and a precipitation minus potential evapotranspiration index, the extent of the region’s exposure to drought is examined. Between 1914 and 1917, the Dry Belt was particularly vulnerable to drought, whereas after 1928, a much larger region known as the Palliser Triangle covering most of southern Alberta and Saskatchewan was much more exposed to drought. These droughts provoked major institutional adaptation, in particular the establishment of the Special Areas Board by the Government of Alberta, and the creation of the Prairie Farm Rehabilitation Administration by the Government of Canada. Both organizations have proved to be relatively permanent public adaptations to the natural hazard of drought in the region. Moreover, these earlier experiences with prolonged drought as well as institution-building may be of value in helping the residents of the Palliser Triangle adapt to predicted climate changes in the future as well as anticipate some of the barriers to effective institutional adaptation.  相似文献