ALBERT

Notes: CO2–CH4 fluid inclusions are present in anatectic layer-parallel leucosomes from graphite-bearing metasedimentary rocks in the Skagit migmatite complex, North Cascades, Washington. Petrological evidence and additional fluid inclusion observations indicate, however, that the Skagit Gneiss was infiltrated by a water-rich fluid during high-temperature metamorphism and migmatization.CO2-rich fluid inclusions have not been observed in Skagit metasedimentary mesosomes or melanosomes, meta-igneous migmatites, or unmigmatized rocks, and are absent from subsolidus leucosomes in metasedimentary migmatites. The observation that CO2-rich inclusions are present only in leucosomes interpreted to be anatectic based on independent mineralogical and chemical criteria suggests that their formation is related to migmatization by partial melting. Although some post-entrapment modification of fluid inclusion composition may have occurred during decompression and deformation, the generation of the CO2-rich fluid is attributed to water-saturated partial melting of graphitic metasedimentary rocks by a reaction such as biotite + plagioclase + quartz + graphite ± Al2SiO5+ water-rich fluid = garnet + melt + CO2–CH4. The presence of CO2-rich fluid inclusions in leucosomes may therefore be an indication that these leucosomes formed by anatexis.Based on the inferences that (1) an influx of fluid triggered partial melting, and (2) some episodes of fluid inclusion trapping are related to migmatization by anatexis, it is concluded that a free fluid was present at some time during high-temperature metamorphism. The infiltrating fluid was a water-rich fluid that may have been derived from nearby crystallizing plutons. Because partial melting took place at pressures of at least 5 kbar, abundant free fluid may have been present in the crust during orogenesis at depths of at least 15 km.

Notes: The solid-solid reaction magnesiocarpholite = sudoite + quartz has been bracketed between 350 and 500°C, 6.3 and 7.8 kbar. Because it is impossible to synthesize end-member sudoite, all experiments were carried out using natural minerals as starting materials. Although mineral compositions were very close to those of the end-members, the effect of the fluorine content in carpholite was significant. Particularly in those experiments where sudoite grows at the expense of carpholite, electron microprobe analysis of the run products shows that a more stable F-rich carpholite crystallizes too, and consumes the fluorine released in solution by the breakdown of the original carpholite.Our experimental results are combined, through a thermodynamic analysis, with a previous data set and with previous experimental data concerning the relative stability of chlorite, talc and magnesiocarpholite with excess of quartz and water as a function of P–T and AlAl(SiMg)-1 substitutions in phyllosilicates. This allows us to constrain the feasible thermodynamic parameters (H°f, sud; S° sud) and (H°f,car; S°car) for the Mg end-members. Using the partition coefficients calculated from natural parageneses, we have computed a petrogenetic grid for the system FeO–MgO–Al2O3–SiO2–H2O. It demonstrates that parageneses involving sudoite and carpholite can be used as indicators of P–T conditions, up to 600° C, 8 kbar for sudoite, and at higher pressure for carpholite.

Notes: Regional-scale mapping of index-mineral isograds in mafic units of the early Proterozoic Cape Smith Thrust Belt (northern Québec) has revealed contrasting pressure-temperature regimes associated with two distinct structural domains. In the southern domain, crustal thickening was accomplished by early, piggy-back thrust faults. Isograds cross-cut the thrusts, indicating that thermal-peak mineral growth outlasted deformation associated with early imbrication. Mineral zones are: (1) actinolite (Act) + albite (Alb); (2) hornblende (Hbl) + Act + Alb; (3) Hbl + Act + oligoclase (Oli); (4) Hbl + Oli; and (5) garnet (Grt) or clinopyroxene + Hbl + Oli-andesine. The oligoclase isograd occurs at higher grade than the hornblende isograd, a sequence typical of medium-pressure terranes (5–7 kbar). An Hbl-Alb bathograd. calibrated from mixed-volatile equilibria in the NCMASH-CO2 model system, suggests minimum pressures of about 5.4 kbar.Metamorphism in the northern domain was a consequence of re-imbrication, by means of out-of-sequence thrust faults active during and after peak metamorphic conditions. Mineral growth was coeval with thrusting, as documented by the syn-kinematic garnet porphyroblasts. Compared to the southern domain, a different sequence of isograds in mafic rocks shows that the albite-oligoclase transition takes place in the garnet zone. Based on thermobarometry in garnet-hornblende rocks, the oligoclase isograd occurs in a temperature range of 525–600°C, typical of high-pressure terranes (7–10 kbar). Calibrated bathograds for the Hbl-Ms-Alb and Grt-Alb bathozonal assemblages, respectively in the KNCMASH-CO2 and NCMASH model systems, indicate minimum pressures in the northern domain of 6.7 and 8.5 kbar. Higher-pressure series for this domain are explained by out-of-sequence thrusts exposing deeper crustal levels. For similar structural levels, only minor amounts of syn-deformational uplift (1–2 kbar and 50–75°C) are recorded in metabasites of this domain, compared to results in adjacent metapelites of the area (essentially isothermal uplift of 3–5 kbar).RESUME La bande du Cap Smith (nord du Québec) est une ceinture de chevauchement d'ǎge protérozoique inférieur, dominée par des roches mafiques. La cartographie d'isogrades à minéraux indicateurs dans les unités mafiques de la ceinture a révelé deux régimes contrastes de pression–température, chacun associéà des épisodes distincts d'épaississement crustal. Dans le domaine sud, des failles de chevauchement en-série sont responsables pour l'empilement tectonique. Les isogrades recoupent les failles, indiquant que l'apogée thermale a suivi l'emplacement initial des nappes de charriage. Les zones minérales sont: (1) actinote (Act) + albite (Alb); 2) hornblende (Hbl) + Act + Alb; (3) Hbl + Act + oligoclase (Oli); 4) Hbl + Oli; et (5) grenat (Grt) où clinopyroxene + Hbl + Oli-andésine. L'isograde d'oligoclase apparaǐt à plus haute température que l'isograde d'hornblende, une séquence typique des terrains de pressions moyennes. Un bathograde Hbl–Alb, calibréà partir d'équilibre mixte de volatiles dans le système NCMASH–CO2, suggère des pressions minimales d'environ 5.4 kbar.Le métamorphisme dans le domaine nord de la ceinture a été le résultat d'une réimbrication, causé par des chevauchements hors-série actifs pendant et après l'apogée thermale. La croissance minérale fǔt synchrone au chevauchement, documentée par des porphyroblastes de grenat syn-cinénatique. Comparé au domaine sud, une différente séquence d'isogrades dans les métabasaltes montre que la transition albite–oligoclase se situé dans la zone à grenat. Par la thermobarométrie dans les roches à grenat–hornblende l'isograde d'oligoclase se situe dans un écart de température de 525–600°C, typique des terrains de hautes pressions (7–10 kbar). Des bathogrades calibrés pour les assemblages bathozonales Hbl–Ms–Alb et Grt–Alb, rcspectivement dans les systèmes KNCMASH–CO2 et NCMASH, indique des pression minimales pour le domaine nord de 6.7 et 8.5 kbar. Une zonégraphie à plus haute pression pour ce domaine est expliquée par des chevauchements hors-série exposant des niveaux plus inférieurs de la croǔte imbriqué. Pour des niveaux structuraux similaries, des soulèvements syn-métamorphiques mineurs sont enregistrés dans les métabasaltes (1–2 kbar et 50–75°C), comparés aux métapélites adjacentes avec un soulèvement (essentially isothermal uplift of 3–5 kbar.

Notes: Mid-Cretaceous granulite gneisses crop out in a narrow belt in the Cucamonga region of the south-eastern foothills of the San Gabriel Mountains, southern California. Interlayered mafic granulites and pelitic, carbonate, calc-silicate and quartzofeldspathic metasediments record hornblende granulite subfacies metamorphism at approximately 8 kbar and 700–800°C. Regional deformation and formation of banded gneisses ceased by c. 108 Ma. although mafic-intermediate magmatism and high-grade metamorphism continued locally as late as c. 88 Ma. Garnet zoning in metapelitic gneisses suggests that peak metamorphism was followed locally by a period of near-isobaric cooling, but this interpretation requires diachronous cooling of the granulite belt which cannot be demonstrated without detailed thermo-chronological data. It is more likely that the entire terrane remained at granulite facies P–T conditions until 88 Ma, followed by rapid uplift associated with juxtaposition against adjacent middle and upper crustal arc terranes. Uplift occurred between c. 88 and 78 Ma at rates of approximately 1–2 km Ma-1. The geotectonic evolution of the Cucamonga granulites is similar to mid-Cretaceous high-P granulites in the Sierra Nevada and Salinian block of central California. Late Cretaceous uplift common to these granulites may provide an important tectonic link between dismembered Mesozoic batholithic terranes in the California Cordillera.

Notes: In this paper we analyze the properties of price equilibria in a duopoly market where firms sell vertically differentiated products, consumers being uncertain about which firm sells which quality. Both existence and properties of price equilibria are characterized by the beliefs of the consumers' population about the distribution of quality between firms.

Notes: The survey classifies economic theories of the firm into four categories based on the level of aggregation in economic models: (1) neoclassical, (2) industrial organization, (3) contractual, and (4) organizational incentive. Economic theories of the firm are evaluated on the basis of their potential application to problems of management decision making. The survey suggests that a management perspective can be useful in developing an integrated theoretical analysis of the firm that addresses both competitive strategy and organizational design.

Notes: Commitment: The Dynamic of Strategy, by Pankaj Ghemawat.Ghemawat's Commitment makes recent results in game-theoretic industrial organization accessible and useful to practitioners in the field of strategic management. This book contributes to the management strategy literature on two levels. On a conceptual level, Ghemawat strives to isolate “commitment” as the sole explanation of persistent differences in firm performance. On a more pragmatic level, he provides a framework intended to aid managers in making commitment-intensive decisions. It is with respect to how well he achieves these two distinct goals that I evaluate Ghemawat's contribution. In addition, I review briefly the book's content, and I compare Ghemawat's approach to some alternative approaches familiar to scholars and practitioners of strategic management.

Notes: In the present paper, we relate the extent of job security offered to incumbent managers to the extent of competition among firms in the product market, where the extent of job security is measured by the probability that an incumbent manager continues to be employed by his current firm and the extent of competition is measured by the degree of differentiation between competing brands. We demonstrate that when competition between firms intensifies and “on-the-job training” is relatively more conducive to reducing the variable costs of production, firms tend to offer reduced (increased) job security to incumbent managers, provided that the degree of differentiation between competing products is sufficiently large (small), respectively. If “on-the-job training” is relatively more conducive to reducing the fixed costs of production, however, the previous result is reversed.

Notes: This paper presents a model of strategic product choice when consumer preferences combine features of both horizontal and vertical product differentiation. Consumers disagree on what amount of a “special” characteristic makes for a better product, but those who prefer more of this attribute are willing to pay more for it. Within this demand structure, I examine the advantages of first-mover firms. I find that such firms typically do best in markets where the maximum degree of product differentiation is limited by preferences rather than technology. These are “niche markets”. Follower firms do better in markets in which the range of preferences is broad relative to the span of feasible goods.

Notes: This paper reexamines Grossman and Hart's (1980) insight into how the free-rider problem excludes an external raider from capturing the increase in value it brings to R firm The inability of the raider to capture any of the surplus depends critically on the assumption of equal and indivisible shareholdings–the one-share-per-shareholder model In contrast, we show that once shareholdings are large and potentially unequal, a raider may capture a significant part of the increase in value Specifically, the free-rider problem does not prevent the takeover process when shareholdings are divisible.

Notes: We analyze the implementation problem faced by firms when trying to collude in the face of asymmetric information about costs. Assuming that transfer payments are possible, we examine the incentive compatibility and individual rationality constraints that must be satisfied by any cartel agreement. Two scenarios are considered. Firms may or may not withdraw from the agreement after each firm's costs become known. If no withdrawal is possible, we find that the monopoly rule is implementable when weak types of individual rationality constraints are required. This contrasts with some results in the literature. If withdrawal is possible, we find a potential conflict between different forms of individual rationality constraints, in particular, between interim and ex post constraints. This conflict disappears in industries with a large number of firms.

Notes: This paper offers a general characterization of the optimal product line prices for a monopolist whose quality of products is initially unknown to consumers. In the focal equilibrium, a monopolist signals a high-quality product line by pricing as if quality were known to be high, but costs of production were higher than they truly are. In a rich set of environments, this characterization implies that the prices of all products are initially distorted upward, with the price distortion being largest for products with the most inelastic demands and/or quality-sensitive production costs. These implications yield predictions for the time path of prices flint are broadly consistent with evidence from the marketing literature. The multidimensional signaling problem is made tractable by the satisfaction of a very simple and powerful single crossing property.

Notes: Managerial behavior that is rational and profit-maximizing sometimes will seem to be overly conservative. If the valuation of innovations contains white noise and the status quo would be preferred to random innovation, then any innovation that does not appear substantially better than the status quo should be rejected, for reasons arising from regression toward the mean. The more successful the firm, the higher is the optimal acceptance threshold and conservative bias. Other things equal, more successful firms will spend less on research, adopt fewer innovations, and be less likely to advance the industry's best practice.

Notes: Identical cases of wine are often auctioned one immediately after another. Ashenfelter (1989) reports that on average, the later lots fetch less. Such a systematic price difference seems anomalous, the more so because it is shown here that rational expectations imply not equal, but rising, prices. Risk aversion is an obvious way of reconciling the evidence with rational behavior. There is an alternative explanation. The auctions observed by Ashenfelter involved a buyer's option, whereby the first-round winner could purchase further cases at the same price. It is shown that this feature may both account for the observed price trajectory and raise seller revenue.

Notes: This paper examines the incentives for integration when the market for consumer durables (hardware) is oligopolistic and the market for complementary services (software) is monopolistically competitive. We find that the equilibrium industry structure will depend on the magnitude of the fixed costs of software development. If the software development costs are relatively large, the equilibrium industry structure is unintegrated, that is, neither hardware firm integrates; if the software development costs are relatively small, the equilibrium industry structure is integrated, that is, both hardware firms integrate. Under the integrated industry structure, hardware profits are lower, less varieties are provided, and hardware prices are lower than under the unintegrated industry structure. The game has a prisoners' dilemma structure when the software development costs are relatively small because of a foreclosure effect. Strategically increasing the number of software varieties provides an avenue for an integrated hardware firm to increase its market share and profits by reducing the number of software varieties available for an unintegrated rival technology. Although consumer surplus is higher under an integrated industry structure, the total surplus associated with the unintegrated industry structure exceeds that of the integrated industry structure.

Notes: Price setting by firms and search by customers is analyzed, relaxing two basic attributes of most search models: price precommitment and agent heterogeneity. Customers are characterized by individual demand functions for a homogeneous good and can choose to employ a threat to search. Firms noncooperatively make pricing decisions by using the individual demand curves under conditions of constant marginal cost. Firms adopt pricing rules that optimally respond to customer search histories. Bargaining power is endogenously assigned. Firms know their common marginal cost; customers, the cost distribution. The unique separating equilibrium is characterized by a lumpy distribution of prices and by heterogeneous shopping behavior by customers giving rise to “shoppers” and “nonshoppers”

Notes: Most research on product positioning supports the idea of differentiation. Product standardization (i.e., minimum differentiation) occurs only under very limiting assumptions. Yet, similar products are often observed in the marketplace. We attempt to restore the case for standardization by using more realistic assumptions than in previous work. We assume that consumers consider not only observable attributes in brand choice, but also attributes that are unobservable by the firms. We find that standardization is an equilibrium when consumers exhibit sufficient heterogeneity along the unobservable attributes under both positioning with exogenously given prices and price competition, We also show that, under insufficient heterogeneity along the unobservable attribute, our results coincide with past research that argues in favor of differentiation.

Notes: This paper proposes an empirical methodology for studying various (implicit or explicit) collusive behaviors on two strategic variables, which are price and advertising, in a differentiated market dominated by a duopoly. In addition to Nash or Stackelberg behaviors, we consider collusion on both variables, collusion on one variable and competition on the other, etc. Using data on the Coca-Cola and Pepsi-Cola markets from 1968 to 1986, full information maximum likelihood estimation of cost and demand functions are obtained allowing for various collusive behaviors. The collusive hypothesis is not rejected, and the best form of collusive behavior is selected via nonnested testing procedures. Using the best model, Lerner indices are computed for both duopolists to provide summary measures of market power. Finally, our approach is contrasted with the conjectural variation approach and is shown to give superior results.

Notes: When the transportation risk posed by shipments of hazardous chemical and radioactive materials is being assessed, it is necessary to evaluate therisks associated with both vehicle emissions and cargo-related risks. Diesel exhaust and fugitive dust emissions from vehicles transporting hazardous shipments lead to increased air pollution, which increases the risk of latent fatalities in the affected population along the transport route. The estimated risk from these vehicle-related sources can often be as large or larger than the estimated risk associated with the material being transported. In this paper, data from the U.S. Environmental Protection Agency's Motor Vehicle-Related Air Toxics Study are first used to develop latent cancer fatality estimates per kilometer of travel in rural and urban areas forall diesel truck classes. These unit risk factors are based on studies investigating the carcinogenic nature of diesel exhaust. With the same methodology, the current per-kilometer latent fatality risk factor used in transportation risk assessments for heavy diesel trucks in urban areas is revised and the analysis expanded to provide risk factors for rural areas and all diesel truck classes. These latter fatality estimates may include, but are not limited to, cancer fatalities and are based primarily on the most recent epidemiological data available on mortality rates associated with ambient air PM-10 concentrations.

Notes: Risk assessors often use different probability plots as a way to assessthe fit of a particular distribution or model by comparing the plotted points to a straight line and to obtain estimates of the parameters in parametric distributions or models. When empirical data do not fall in a sufficiently straight line on a probability plot, and when no other single parametricdistribution provides an acceptable (graphical) fit to the data, the risk assessor may consider a mixture model with two component distributions. Animated probability plots are a way to visualize the possible behaviors of mixture models with two component distributions. When no single parametric distribution provides an adequate fit to an empirical dataset, animated probability plots can help an analyst pick some plausible mixture models for the data based on their qualitative fit. After using animations during exploratory data analysis, the analyst must then use other statistical tools, including but not limited to: Maximum Likelihood Estimation (MLE) to find the optimal parameters, Goodness of Fit (GoF) tests, and a variety of diagnostic plots to check the adequacy of the fit. Using a specific example with two LogNormal components, we illustrate the use of animated probability plots asa tool for exploring the suitability of a mixture model with two component distributions. Animations work well with other types of probability plots, and they may be extended to analyze mixture models with three or more component distributions.

Notes: For carcinogens, this paper provides a quantitative examination of the roles of potency and weight-of-evidence (WOE) in setting permissible exposure limits (PELs) at the U.S. Occupational Safety and Health Administration (OSHA) and threshold limit values (TLVs) at the private American Conference of Governmental Industrial Hygienists (ACGIH). On normative grounds, both of these factors should influence choices about the acceptable level of exposures. Our major objective is to examine whether and in what ways these factors have been considered by these organizations. A lesser objective is to identify outliers, which might be candidates for further regulatory scrutiny. Our sample (N=48) includes chemicals for which EPA has estimated a unit risk as a measure of carcinogenic potency and for which OSHA or the ACGIH has a PEL or TLV. Different assessments of the strength of the evidence of carcinogenicity were obtained from EPA, ACGIH, and the International Agency for Research on Cancer. We found that potency alone explains 49% of the variation in PELs and 62% of the variation in TLVs. For the ACGIH, WOE plays a much smaller role than potency. TLVs set by the ACGIH since 1989 appear to be stricter than earlier TLVs. We suggest that this change represents evidence that the ACGIH had responded to criticisms leveled at it in the late 1980s for failing to adopt sufficiently protective standards. The models developed here identify 2-nitropropane, ethylene dibromide, and chromium as having OSHA PELs significantly higher than predicted on the basis of potency and WOE.

Notes: Much has been written about the development and application of quantitative methods for estimating under uncertainty the long-term radiological performance of underground disposal of radioactive wastes. Until recently, interest has been focused almost entirely on the technical challenges regardless of the role of the organization responsible for these analyses. Now the dialogue between regulators, the repository developer or operator, and other interested parties in the decision-making process receives increasing attention, especially in view of some current difficulties in obtaining approvals to construct or operate deep facilities for intermediate or high-level wastes. Consequently, it is timely to consider the options for regulators’review and evaluation of safety submissions, at the various stages in the site selection to repository closure process, and to consider, especially, the role for performance assessment (PA) within the programs of a regulator both before and after delivery of such a submission. The origins and broad character of present regulations in the European Union (EU) and in the OECD countries are outlined and some regulatory PA reviewed. The issues raised are discussed, especially in regard to the interpretation of regulations, the dangers from the desire for simplicity in argument, the use of regulatory PA to review and challenge the PA in the safety case, and the effects of the relationship between proponent and regulator. Finally, a very limited analysis of the role of PA in public hearings is outlined and recommendations are made, together with proposals for improving the mechanisms for international collaboration on technical issues of regulatory concern.

Notes: In this paper the problem of high-level nuclear waste disposal is viewed as a five-stage, cascaded decision problem. The first four of these decisions having essentially been made, the work of recent years has been focused on the fifth stage, which concerns specifics of the repository design. The probabilistic performance assessment (PPA) work is viewed as the outcome prediction for this stage, and the site characterization work as the information gathering option. This brief examination of the proposed Yucca Mountain repository through a decision analysis framework resulted in three conclusions: (1) A decision theory approach to the process of selecting and characterizing Yucca Mountain would enhance public understanding of the issues and solutions to high-level waste management; (2) engineered systems are an attractive alternative to offset uncertainties in the containment capability of the natural setting and should receive greater emphasis in the design of the repository; and (3) a strategy of “waste management” should be adopted, as opposed to “waste disposal,” as it allows for incremental confirmation and confidence building of a permanent solution to the high-level waste problem.

Notes: Physiologically-based toxicokinetic (PBTK) models are widely used to quantify whole-body kinetics of various substances. However, since they attempt to reproduce anatomical structures and physiological events, they have ahigh number of parameters. Their identification from kinetic data alone is often impossible, and other information about the parameters is needed to render the model identifiable. The most commonly used approach consists of independently measuring, or taking fom literature sources, some of the parameters, fixing them in the kinetic model, and then performing model identification on a reduced number of less certain parameters. This results in a substantial reduction of the degrees of freedom of the model. In this study, we show that this method results in final estimates of the free parameters whose precision is overestimated. We then compared this approach with an empirical Bayes approach, which takes into account not only the mean value, but also the error associated with the independently determined parameters. Blood and breath 2H8- toluene washout curves, obtained in17 subjects, were analyzed with a previously presented PBTK model suitable for person-specific dosimetry. Model parameters with the greatest effect onpredicted levels were alveolar ventilation rate QPC, fat tissue fraction VFC, blood air partition coefficient Kb, fraction of cardiac output to fat Qa/co and rate of extrahepatic metabolismVmax.p. Differences in the measured and Bayesian-fitted values of QPC, VFc and Kb were significant (p 〈 0.05), andthe precision of the fitted values Vmax.p and Qa/co went from 11 ± 5% to 75 ± 170% (NS) and from 8 ± 2% to 9 ± 2% (p 〈 0.05) respectively. The empirical Bayes approach did not result in less reliable parameter estimates: rather, it pointed out that the precision of parameter estimates can be overly optimistic when other parameters in the model, eitherdirectly measured or taken from literature sources, are treated as known without error. In conclusion, an empirical Bayes approach to parameter estimation resulted in a better model fit, different final parameter estimates, and more realistic parameter precisions.

Notes: The US Department of Transportation was interested in the risks associated with transporting Hydrazine in tanks with and without relief devices. Hydrazine is both highly toxic and flammable, as well as corrosive. Consequently, there was a conflict as to whether a relief device should be used or not. Data were not available on the impact of relief devices on release probabilities or the impact of Hydrazine on the likelihood of fires and explosions. In this paper, a Monte Carlo sensitivity analysis of the unknown parameters was used to assess the risks associated with highway transport of Hydrazine. To help determine whether or not relief devices should be used, fault trees and event trees were used to model the sequences of events that could lead to adverse consequences during transport of Hydrazine. The event probabilities in the event trees were derived as functions of the parameters whose effects were not known. The impacts of these parameters on the riskof toxic exposures, fires, and explosions were analyzed through a Monte Carlo sensitivity analysis and analyzed statistically through an analysis of variance. The analysis allowed the determination of which of the unknown parameters had a significant impact on the risks. It also provided the necessary support to a critical transportation decision even though the values of several key parameters were not known.

Notes: During the modernization of the municipal waste incinerator (MWI, maximum capacity of 180,000 tons per year) of Metropolitan Grenoble (405,000 inhabitants), in France, a risk assessment was conducted, based on four tracerpollutants: two volatile organic compounds (benzene and 1,1,1 trichloroethane) and two heavy metals (nickel and cadmium, measured in particles). A Gaussian plume dispersion model, applied to maximum emissions measured at the MWI stacks, was used to estimate the distribution of these pollutants in the atmosphere throughout the metropolitan area. A random sample telephone survey (570 subjects) gathered data on time-activity patterns, according to demographic characteristics of the population. Life-long exposure was assessed as a time-weighted average of ambient air concentrations. Inhalation alone was considered because, in the Grenoble urban setting, other routes of exposure are not likely. A Monte Carlo simulation was used to describe probability distributions of exposures and risks. The median of the life-long personal exposures distribution to MWI benzene was 3.2 · 10−5μg/m3 (20th and 80th percentiles = 1.5 · 10−5 and 6.5 · 10−5μg/m3), yielding a 2.6 · 10−10 carcinogenic risk (1.2 · 10−10 - 5.4 · 10−10). For nickel, the corresponding life-time exposure and cancer risk were 1.8 ·10−4μg/m3 (0.9 ·10−4 - 3.6 ·10−4μg/m3) and 8.6 · 10−8 (4.3 · 10−8 - 17.3 ·10−8); for cadmium they were respectively 8.3 ·10−6μg/m3 (4.0 ·10−6 - 17.6 ·10−6) and 1.5 · (7.2 · 10−9 - 3.1. · 10−8). Inhalation exposure to cadmium emitted by the MWI represented less than 1% of the WHO Air Quality Guideline (5 ng/m3), while there was a margin of exposure of more than 109 between the NOAEL (150 ppm) and exposure estimates to trichloroethane. Neither dioxins nor mercury, a volatile metal, were measured. This could lessen the attributable life-long risks estimated. The minute (VOCs and cadmium) to moderate (nickel) exposure and risk estimates are in accord with other studies on modern MWIs meeting recent emission regulations, however.

Notes: Methods of quantitative risk assessment for toxic responses that are measured on a continuous scale are not well established. Although risk-assessment procedures that attempt to utilize the quantitative information in such data have been proposed, there is no general agreement that these procedures are appreciably more efficient than common quantal dose-response procedures that operate on dichotomized continuous data. This paper points out an equivalence between the dose-response models of the nonquantal approach of Kodell and West(1)) and a quantal probit procedure, and provides results from a Monte Carlo simulation study to compare coverage probabilities of statistical lower confidence limits on dose corresponding to specified additional risk based on applying the two procedures to continuous data from a dose-response experiment. The nonquantal approach is shown to be superior, in terms of both statistical validity and statistical efficiency.

Notes: This paper discusses a successful public involvement effort that addressed and resolved several highly controversial water management issues involving environmental and flood risks associated with an electrical generation facility in British Columbia. It begins with a discussion of concepts for designing public involvement, summarizing research that indicates why individuals and groups may find it difficult to make complex choices. Reasons for public involvement, and the range of current practices are discussed. Next, four principles for designing group decision process are outlined, emphasizing decision-aiding concepts that include “value-focused thinking” and “adaptive management.” The next sections discuss the Alouette River Stakeholder Committee process in terms of objectives, participation, process, methods for structuring values and creating alternatives, information sources, and results. Discussion and conclusions complete the paper.

Notes: During the 1980s, seismic research suggested that Oregon and the City of Portland had a higher risk of a major earthquake than had previously been assumed. In 1993, the State of Oregon adopted a new version of the Oregon Structural Specialty Code, which changed the designation of western Oregon from seismic zone 2b to seismic zone 3. The City of Portland established a program and a Task Force on Seismic Strengthening of Buildings to recommend actions that would encourage upgrading of city buildings. A survey of adult city residents was conducted in April, 1996 to determine public attitudes and opinions about earthquake risks, management and mitigation of earthquake hazards, priorities for protection by strengthening buildings, evaluations of strategies for informing the public about earthquake risks, and support for specific options the city might take to protect citizens against earthquake events. Social and demographic information on individuals and households was also collected. Respondents provided ratings for a wide range of social and environmental risks, provided information on priorities for strengthening key buildings and infrastructure facilities, and answered hypothetical questions about voting for bond measures to pay for city earthquake mitigation programs. Respondents recognized significant risk from earthquakes and supported programs to protect people, especially vulnerable residents such as children and the sick. There was strong support for protecting emergency response capabilities. There was much less support for using public funds to reduce the risks associated with privately owned buildings. There were also some strong pockets of resistance to publicly funded mitigation programs in response to the hypothetical bond measures.

Notes: From a comprehensive search of the literature, the hormesis phenomenon was found to occur over a wide range of chemicals, taxonomic groups, and endpoints. By use of computer searches and extensive cross-referencing, nearly 3000 potentially relevant articles were identified. Evidence of chemical and radiation hormesis was judged to have occurred in approximately 10oO of these by use of a priori criteria. These criteria included study design features (e.g., number of doses, dose range), dose-response relationship, statistical analysis, and reproducibility of results. Numerous biological endpoints were assessed, with growth responses the most prevalent, followed by metabolic effects, reproductive responses, longevity, and cancer. Hormetic responses were generally observed to be of limited magnitude with an average maximum stimulation of 30 to 60 percent over that of the controls. This maximum usually occurred 4- to 5-fold below the NOAEL for a particular endpoint. The present analysis suggests that hormesis is a reproducible and generalizable biological phenomenon and is a fundamental component of many, if not most, dose-response relationships. The relatively infrequent observation of homesis in the literature is believed to be due primarily to experimental design considerations, especially with respect to the number and range of doses and endpoint selection. Because of regulatory considerations, most toxicologic studies have been carried out at high doses above the low-dose region where the hormesis phenomenon occurs.

Notes: As the use of digital computers for instrumentation and control of safety-critical systems has increased, there has been a growing debate over the issue of whether probabilistic risk assessment techniques can be applied to these systems. This debate has centered on the issue of whether software failures can be modeled probabilistically. This paper describes a “context-based” approach to software risk assessment that explicitly recognizes the fact that the behavior of software is not probabilistic. The source of the perceived uncertainty in its behavior results from both the input to the software as well as the application and environment in which the software is operating. Failures occur as the result of encountering some context for which the software was not properly designed, as opposed to the software simply failing “randomly.” The paper elaborates on the concept of “error-forcing context” as it applies to software. It also illustrates a methodology which utilizes event trees, fault trees, and the Dynamic Flowgraph Methodology (DFM) to identify “error-forcing contexts” for software in the form of fault tree prime implicants.

Notes: Australian state and federal agencies use a broad range of methods for setting conservation priorities for species at risk. Some of these are based on rule sets developed by the International Union for the Conservation ofNature, while others use point scoring protocols to assess threat. All of them ignore uncertainty in the data. In this study, we assessed the conservation status of 29 threatened vascular plants from Tasmania and New South Wales using a variety of methods including point scoring and rule-based approaches. In addition, several methods for dealing with uncertainty in the data were applied to each of the prioritysetting schemes. The results indicatethat the choice of a protocol for setting priorities and the choice of the way in which uncertainty is treated may make important differences to the resulting assessments of risk. The choice among methods needs to be rationalized within the management context in which it is to be applied. These methods are not a substitute for more formal risk assessment.

Notes: This paper describes a multi-stakeholder process designed to assess thepotential health risks associated with adverse air quality in an urban industrial neighborhood. The paper briefly describes the quantitative health risk assessment conducted by scientific experts, with input by a grassroots community group concerned about the impacts of adverse air quality on theirhealth and quality of life. In this case, rather than accept the views of the scientific experts, the community used their powers of perception toadvantage by successfully advocating for a professionally conducted community health survey. This survey was designed to document, systematically and rigorously, the health risk perceptions community members associated with exposure to adverse air quality in their neighborhood. This paper describes theinstitutional and community contexts within which the research is situated as well as the design, administration, analysis, and results of the community health survey administered to 402 households living in an urban industrial neighborhood in Hamilton, Ontario, Canada. These survey results served tolegitimate the community's concerns about air quality and tohelp broaden operational definitions of ‘health.’ In addition, the resultsof both healthrisk assessment exercises served to keep issues of air quality on the localpolitical agenda. Implications of these findings for our understanding of theenvironmental justice process as well as the ability of communitiesto influence environmental health policy are discussed.

Notes: Psychometric data on risk perceptions are often collected using the method developed by Slovic, Fischhoff, and Lichtenstein, where an array of risk issues are evaluated with respect to a number of risk characteristics, such as how dreadful, catastrophic or involuntary exposure to each risk is. The analysis of these data has often been carried out at an aggregate level, where mean scores for all respondents are compared between risk issues. However, this approach may conceal important variation between individuals, and individual analyses have also been performed for single risk issues. This paper presents a new methodological approach using a technique called multilevel modelling for analysing individual and aggregated responses simultaneously, to produce unconditional and unbiased results at both individual and aggregate levels of the data. Two examples are given using previously published data sets on risk perceptions collected by the authors, and results between the traditional and new approaches compared. The discussion focuses on the implications of and possibilities provided by the new methodology.

Notes: The Whitestone Anorthosite (WSA), located in southern Ontario, underwent granulite facies metamorphism during the Grenville orogeny at 1.16 Ga. During the waning stages of metamorphism fluids infiltrated the outer portions of the anorthosite and promoted the formation of an envelope comprised of upper amphibolite facies mineral assemblages. Also, this envelope corresponds to portions of the anorthosite that underwent deformation related to movement along a high-grade ductile shear zone. Samples from this portion of the anorthosite (the margin) contain CO2-rich inclusions in plagioclase porphyroclasts (relict igneous phenocrysts), matrix plagioclase and garnet. These inclusions have features which normally are interpreted as indicating that they are texturally primary, but they have relatively low CO2 densities (0.61–0.95 g cm-3). Plagioclase from the anorthosite interior contains texturally secondary inclusions with relatively high CO2 densities (generally from 0.99 to 1.10 g cm-3). The high CO2 densities suggest that the inclusions in the plagioclase of the anorthosite core formed prior to inclusions in porphyroclast minerals of the outer portions of the anorthosite, an interpretation that is apparently inconsistent with inclusion textures. This apparent paradox indicates that most fluid inclusions from the anorthosite margin were formed during, or were modified by, the dynamic recrystallization that affected this portion of the WSA. In either case, late formation or modification, the texturally primary fluid inclusions do not contain pristine samples of the peak metamorphic fluid. Furthermore, because shear-related deformation is apparently associated with entrapment of the lowest fluid densities, some strain localization persisted to relatively low temperatures (e.g. less than approximately 500° C). These results constrain a part of the retrograde P–T path for this portion of the Grenville Orogen to temperatures of approximately 400–500° C at pressures of approximately 1–2 kbar.

Notes: A suite of metapelites, charnockites, calc-silicate rocks, quartzo-feldspathic gneisses and mafic granulites is exposed at Garbham, a part of the Eastern Ghats granulite belt of India. Reaction textures and mineral compositional data have been used to determine the P–T–X evolutionary history of the granulites. In metapelites and charnockites, dehydration melting reactions involving biotite produced quartzofeldspathic segregations during peak metamorphism. However, migration of melt from the site of generation was limited. Subsequent to peak metamorphism at c. 860° C and 8 kbar, the complex evolved through nearly isothermal decompression to 530–650° C and 4–5 kbar. During this phase, coronal garnet grew in the calc-silicates, while garnet in the presence of quartz broke down in charnockite and mafic granulite. Fluid activities during metamorphism were internally buffered in different lithologies in the presence of a melt phase. The P–T path of the granulites at Garbham contrasts sharply with the other parts of the Eastern Ghats granulite belt where the rocks show dominantly near-isobaric cooling subsequent to peak metamorphism.

Notes: Staurolite porphyroblasts, 1.5–8cm in length and 0.3–2cm in width, in the Littleton Schist at Bolton, Connecticut, contain curved quartz inclusion trails which document synkinematic rotations of at least 135°. The orientations of long axes of these staurolite crystals define a weak preferred orientation in a plane approximately parallel to the external foliation. Serial sections of four differently orientated crystals and U-stage measurements of the orientations of their inclusion trails demonstrate that the inflection hinge line and the statistical ‘symmetry axis’ characterizing the foliation within a porphyroblast are unrelated to the orientations of external crenulations and are, in all cases, parallel to the long axis of the porphyroblast. The cumulative rotation reflected in the curvature of the inclusion trails is a maximum in a c-axis section through the initial core of a crystal. The amount of rotation about the c-axis decreases linearly along the length of the crystal away from the nucleation site.The sense and amount of rotation recorded by a porphyroblast is related to its orientation. A tightly constrained transition from clockwise to anticlockwise rotation defines a slip direction that coincides with the preferred orientation of the staurolite c-axes. The total rotation reflected by the inclusion trails increases as a function of the angle between the c-axes of the staurolite crystals and the slip direction.Initially random staurolite porphyroblasts rotated during growth, as a consequence of laminar shear in the surrounding viscous matrix. This interpretation is quantitatively consistent with: the staurolite preferred orientation; its coincidence with the apparent slip direction; the correlation between both the sense and the amount of rotation and the orientation of the long axis of the porphyroblast; and the twisted conical shape of the family of surfaces defined by the inclusion trails.

Notes: Silica-deficient sapphirine-bearing rocks occur as an enclave within granulite facies Proterozoic gneisses and migmatites near Grimstad in the Bamble sector of south-east Norway (Hasleholmen locality). The rocks contain peraluminous sapphirine, orthopyroxene, gedrite, anthophyllite, sillimanite, sapphirine, corundum, cordierite, spinel, quartz and biotite in a variety of assemblages. Feldspar is absent.Fe2+/(Fe2++ Mg) in the analysed minerals varies in the order: spinel 〉 gedrite ≥ anthophyllite ≥ biotite 〉 sapphirine〉orthopyroxene 〉 cordierite.Characteristic pseudomorph textures indicate coexistence of orthopyroxene and sillimanite during early stages of the reaction history. Assemblages containing orthopyroxene-sillimanite-sapphirine-cordierite-corundum developed during a high-pressure phase of metamorphism and are consistent with equilibration pressures of about 9 kbar at temperatures of 750–800°C. Decompression towards medium-pressure granulite facies generated various sapphirine-bearing assemblages. The diagnostic assemblage of this stage is sapphirine-cordierite. Sapphirine occurs in characteristic symplectite textures. The major mineralogical changes can be described by the discontinuous FMAS reaction: orthopyroxene + sillimanite → sapphirine + cordierite + corundum.The disequilibrium textures found in the Hasleholmen rocks are characteristic for reactions which have been in progress but then ceased before they run to completion. Textures such as reaction rims, symplectites, partial replacement, corrosion and dissolution of earlier minerals are characteristic of granulite facies rocks. They indicate that, despite relatively high temperatures (700–800° C), equilibrium domains were small and chemical communication and transport was hampered as a result of dry or H2O-poor conditions.

Notes: The Vredefort dome (2.0 Ga) represents the central uplift of a very large impact structure. This uplift exposed a nearly complete cross-section through the continental crust in the region, which is 25–30 km thick. Two metamorphic events took place at about the same time as the impact. The first event, so-called static metamorphism, is pre-impact and produced lithologies varying from low-grade shale to high-grade hornfels. It resembles contact metamorphism by its lack of schistosity, but is more regional as it extends over a large area and is not associated with large intrusions.The second event, the post-shock metamorphism, is responsible for the recrystallization of the shock features. The investigation of this event has been focused on the degree of alteration of the coesite-stishovite-bearing pseudotachylite veins that formed during the transit of the shock wave. These high-pressure silica polymorphs are only present in the upper part of the stratigraphic sequence; downward they have been converted to fibrous quartz. At the highest grade, the fibrous quartz is in turn replaced by triple-junctioned mosaic quartz. The post-shock metamorphism was generated by the heat of the rock before shock, plus the heat released by the shock wave. The isograds, plotted on a map, can be translated into depth of burial and therefore provide valuable information regarding the geological setting immediately before impact. At the time of impact, the rocks were relatively cool and the static metamorphism had ceased with several tens of millions of years separating the two metamorphic events. The static metamorphism was probably caused by continental crustal extension in a stress-free environment and the lack of deformation is probably due to rapid uplift during the later stages of the impact event.

Notes: The P–T paths for metamorphic complexes from the Precambrian shields and fold belts of different ages may result from advection, i.e. one-cycle convective processes in the lithosphere. This conclusion has been exemplified by the metamorphic evolution of several well-known complexes, for which an advective model can be successfully applied. Numerical simulations of the above processes in terms of Newtonian rheology by using a two-dimensional finite element program have been conducted.Two representative models for intracontinental gravitational ordering initiated presumably by mantle activity are considered: (i) a thermally activated multi-layered rhythmic sequence and (ii) huge rising diapiars causing circulation, in which crustal lithologies underwent high-P metamorphism (above 10–15 kbar) and subsequent ascent toward the Earth's surface.

Notes: The Waterman Metamorphic Complex of the central Mojave Desert was exposed as a consequence of early Miocene detachment-dominated extension. However, it has evidence consistent with a more extensive geological history that involves collision of a crustal fragment(s), tectonic thickening by overthrusting and two periods of extension. The metamorphic complex contains granitoid intrusives and felsic mylonitic gneisses as well as polymetamorphic rocks that include marble, calc-silicate, quartzite. mafic granulite, pyribolite, amphibolite, migmatite and biotite schist. The latter group of rocks was affected by an initial series of high-grade metamorphic events (M1 and M2) and a localized lower grade overprint (M3). The initial metamorphism (M1) can be separated into two stages along its high-grade P–T path: M1a, a granulite facies metamorphism at 800–850° C and 7.5–9 kbar and Mlb, an upper amphibolite facies overprint at 750–800° C and 10–12 kbar. M1a developed mineral assemblages and textures consistent with granulite facies conditions at a reduced activity of H2O and is associated with intense ductile deformation (D1) and minor local partial melting. M1b overprinted the granulite assemblages with a series of hydrous phases under conditions of increasing pressure and H2O activity and is accompanied by little or no deformation. M2 developed at lower pressures and temperatures (650–750° C, 4.5–5.5 kbar) and is distinguished by a second local overprint of hydrous phases that reflects an input of aqueous fluids probably associated with the intrusion of a series of granitic dykes and veins. Effects of M3 are confined to the Mitchel detachment zone, an anastomosing early Miocene detachment fault, and are characterized by local ductile/brittle deformation (D2) of the pre-existing high-grade rocks and granitoid intrusives and by the production of mylonites and mylonitic gneisses under greenschist facies conditions (300–350° C, 3–5 kbar). The initial overprint (M1a) represents metamorphism, devolatilization and minor partial melting of supracrustal rocks under granulite facies conditions as a consequence of tectonic and, possibly, magmatic thickening. The increasing pressure transition of M1a to M1b reflects a period of continued compressional tectonism, thrusting and influx of H2O, in part, locally related to crystallization of partial melts. The near isothermal decompression between M1b and M2 probably represents a pre-112-Ma extensional episode that may have been the result of a decompressional readjustment of a thickened crust. Following the initial extensional event, the metamorphic complex remained at depths of 10–17 km for at least 90 Ma until it was uplifted following Miocene extension. M3 develops locally in response to this second extensional period resulting from the early Miocene detachment faulting.

Notes: The High Himalayan Crystalline Sequence in north-central Nepal is a 15-km-thick pile of metasediments that is bound by the Main Central Thrust to the south and a normal fault to the north. The Langtang section through the metasediments shows an apparent inversion of metamorphic isograds with high-P, kyanite-grade rocks exposed beneath low-P, sillimanite-grade rocks. Textural evidence confirms that the observed inversion is a result of a polyphase metamorphic history and phase equilibria studies indicate that thermal decoupling has occurred within a mechanically coherent section of crust. Rocks now exposed at the base of the High Himalayan thrust sheet underwent Barrovian regional metamorphism (M1) prior to 34 Ma in the early stages of the Himalayan orogeny, recording metamorphic conditions of T= 710 ± 30° C, P= 9 ± 1 kbar. After the activation of the Main Central Thrust, which emplaced these metapelites southwards onto the lower grade Lesser Himalayan formations, the upper part of the thrust sheet was overprinted by a second heating event (M2), resulting in sillimanite-grade metamorphism and anatexis of metapelites at T= 760 ± 30° C, P= 5.8 ± 0.4 kbar between 17 and 20 Ma. Crustally derived, leucogranite magmas have been emplaced into low-grade Tethyan sediments on the hangingwall of the normal fault that bounds the northern limit of the metapelitic sequence.The cause of the selective heating of the upper section of the metasediments during M2 cannot be reconciled with either post-thrusting thermal relaxation or advection models. The cause of M2 remains problematical but it is suggested that heat focusing has occurred at the top of the High Himalayan Crystalline Sequence as a result of movement on the normal fault blanketing metapelites of high heat productivity with low-grade sediments of low thermal conductivity. This model implies that the normal fault was active before M2, consistent with decompression textures that formed during, or shortly after, sillimanite-grade metamorphism.

Notes: The two major Early to Middle Palaeozoic tectonic/metamorphic events in the northern Appalachians were the Taconian (Middle to Late Ordovician) in central to western areas and the Acadian (Late Silurian to early Middle Devonian) in eastern to west-central areas. This paper presents a model for the Acadian orogenic event which separates the Acadian metamorphic realm into eastern and western belts based on distinctively different styles. We propose that the Acadian metamorphism in the east was the delayed consequence of Taconian back-arc lithospheric modification. East of the Taconian island arc, thick accumulations of Late Ordovician and Silurian sediments, coupled with plutons rising along a magmatic arc, produced crustal thermal conditions appropriate for anomalously high-T, low-P metamorphism accompanied by major crustal anatexis. In this zone, upward melt migration was coupled with subsequent E-W crustal shortening (possibly due to outboard collision with the Avalon terrane) to produce mechanical conditions that favoured formation of fold and thrust nappes and resultant tectonic thickening to the west (and probably to the east as well).The basis for the distinction between the Eastern and Western Acadian events lies in the contrasting styles of metamorphism accompanying each. Evidence for contrasting metamorphic styles consists of (1) estimated metamorphic field gradients (MFGs) based on thermobarometric studies, and (2) petrological evidence for contrasting P–T trajectories. West of the Acadian metamorphic front, the Taconian zone has an MFG in which peak temperatures of 400-600° C were reached at pressures of about 4–6 kbar, with both P and T increasing to the east. Near its western edge, the Western Acadian metamorphic overprint has a similar MFG to the Taconian, and is mainly discriminated by 40Ar/39Ar dating and microtextural evidence. East of this narrow zone, the Western Acadian overprint is characterized by progressively higher temperatures (600–725° C) and pressures (6.5–10 kbar, or more) to the east, yielding an overall MFG that lies along, or slightly above, the kyanite–sillimanite boundary on a P–T diagram. There is little or no plutonism accompanying Western Acadian metamorphism.In contrast, thermobarometry in the Eastern Acadian, east of the Bronson Hill Belt, yields high-T, intermediate-P conditions for the highest grade rocks known in New England: T= 650–750° C, P= 4.5–6.5 kbar for granulite facies assemblages which apparently formed along an ‘anticlockwise’P–T path. The Bronson Hill Belt lies geographically between the Eastern and Western Acadian zones and shows transitional petrological behaviour: anomalously high temperatures at intermediate pressures, but a ‘clockwise’ path with decompression cooling.Radiometric dating indicates peak Taconian conditions may have been achieved as early as 475 Ma in the Taconian hinterland and as late as 445 Ma in the Taconian foreland (including the Taconic allochthons). Eastern Acadian magmatism may have started as early as 425 Ma, and most nappe-stage deformation and metamorphism in the Eastern Acadian zone appears to have ended by about 410 Ma. Tectonic thickening in the Western Acadian (including the western counterparts of the nappe-stage deformation documented in the Eastern Acadian) must pre-date attainment of peak metamorphic conditions dated at 395–385 Ma. Dome-stage deformation clearly post-dates peak metamorphism and deforms metamorphic isograds. The end of Western Acadian deformation is well constrained by 370-375 Ma radiometric ages of late pegmatites and granitoids which cross-cut all structures.

Notes: A suite of migmatites in uppermost amphibolite facies schists of the Koettlitz Group exposed in the Taylor Valley, Antarctica, provides direct evidence of the behaviour of partially molten rock during syn-anatectic deformation. The geometry of the migmatites is directly related to their position relative to the hinge of a kilometre-scale antiform. Migmatitic rocks on the fold limbs are characterized by extensional shears and fractures, filled with leucosome material, that intersect the pervasive foliation and millimetre-thick stromatic leucosomes. Vein- and dyke-like leucosomes become more common and thicker from the limb to the hinge region of the antiform. Rocks characterized by high leucosome-to-rock ratios near the antiform hinge are xenolithic in appearance. Major parasitic folds within the hinge contain leucogranite ‘microplutons’ up to 50 m across beneath refractory ‘cap-rock’ layers.Angular boudinage structures in schists surrounded by leucosomes indicate a relatively low yield strength in the leucosome, which is compatible with a molten rather than solid leucosome. Leucogranite-bearing extensional shears and fractures indicate that repeated extensional fracturing and shearing promoted by high fluid (melt) pressure is an important mechanism of melt segregation. Dilation in the hinges of developing folds aids the migration of melt into fold hinges and the development of 10–50-m-wide ‘microplutons’ of xenolith-rich leucogranite.Lack of vapour-absent melting and consequent low melt-to-rock ratios allowed the Koettlitz Group to maintain its structural coherency on a kilometre scale. Consequently, leucosome ‘microplutons’ did not exceed 50 m in width, and therefore observed leucosomes have not contributed to the development of adjacent plutonic-scale granitoids.

Notes: Systematic mapping of a transect along the well-exposed shores of Georgian Bay, Ontario, combined with the preliminary results of structural analysis, geochronology and metamorphic petrology, places some constraints on the geological setting of high-grade metamorphism in this part of the Central Gneiss Belt. Correlations within and between map units (gneiss associations) have allowed us to recognize five tectonic units that differ in various aspects of their lithology, metamorphic and plutonic history, and structural style. The lowest unit, which forms the footwall to a regional decollement, locally preserves relic pre-Grenvillian granulite facies assemblages reworked under amphibolite facies conditions during the Grenvillian orogeny. Tectonic units above the decollement apparently lack the early granulite facies metamorphism; out-of-sequence thrusting in the south produced a duplex-like structure. Two distinct stages of Grenvillian metamorphism are apparent. The earlier stage (c. 1160–1120 Ma) produced granulite facies assemblages in the Parry Sound domain and upper amphibolite facies assemblages in the Parry Island thrust sheet. The later stage (c. 1040–1020 Ma) involved widespread, dominantly upper amphibolite facies metamorphism within and beneath the duplex. Deformation and metamorphism recently reported from south and east of the Parry Sound domain at c. 1100–1040 Ma have not yet been documented along the Georgian Bay transect. The data suggest that early convergence was followed by a period of crustal thickening in the orogenic core south-east of the transect area, with further advance to the north-west during and after the waning stages of this deformation.

Notes: The Main Zone of the Hidaka metamorphic belt is an imbricate stack of crustal material derived from an island arc in which a sequence of units with increasing metamorphic grade from low to high structural levels is exposed. The basal part of the metamorphic sequence underwent granulite facies metamorphism with peak P–T conditions of 7kbar, 870°C. In this zone pelitic granulite includes leucosomes which consist mainly of orthopyroxene-plagioclase-quartz.To test whether the leucosome was derived by partial melting of the surrounding pelite, melting experiments of the pelitic granulite were carried out for water-saturated and dry systems at 7 kbar and 850°C. The chemical composition of the leucosome produced during these runs shows a peraluminous S-type tonalitic affinity and is located very close to the tie-line between the average melts produced in water-saturated systems and the average composition of the residual orthopyroxene + plagioclase. This therefore suggests that the lecosome in pelitic granulite was formed by incipient anatexis at close to the highest P–T condition of the Main Zone.The age of the crustal anatexis is determined by the Rb-Sr whole rock isochron method for garnet-cordierite-biotite gneiss (host rock), garnet-orthopyroxene-cordierite gneiss (restite) and S-type tonalite (melt). This gives an age of 56.0 Ma with an initial 87Sr/86Sr ratio of 0.705711. The S-type tonalite magmas that form large intrusive masses in the Main Zone were probably generated by crustal anatexis in deeper parts of the crust at the same time (late Palaeocene).

Notes: Caledonian eclogite facies shear zones developed from Grenvillian garnet granulite facies anorthosites and gabbros in the Bergen Arcs of western Norway allow direct investigation of the relations between macroscopic structures and crystallographic preferred orientation (CPO) in lower continental crust. Field relations on the island of Holsnøy show that the eclogites formed locally from granulite facies rocks by progressive development of: (1) eclogite adjacent to fractures; (2) eclogite in discrete shear zones (〉 2 m thick); (3) eclogite breccia consisting of 〉80% well-foliated eclogite that wraps around rotated granulite blocks; and (4) anastomosing, subparallel, eclogite facies shear zones 30–100 m thick continuous over distances 〉 1 km within the granulite terrane. These shear zones deformed under eclogite facies conditions at an estimated temperature of 670 ± 50°C and a minimum pressure of 1460 MPa, which corresponds to depths of 〉55 km in the continental crust. Detailed investigation of the major shear zones shows the development of a strong foliation defined by the shape preferred orientation of omphacite and by alternating segregations of omphacite/garnet-rich and kyanite/zoisite-rich layers. A consistent lineation throughout the shear zones is defined by elongate aggregates of garnet and omphacite. The CPO of omphacite, determined from five-axis universal stage measurements, shows a strong b-axis maximum normal to foliation, and a c-axis girdle within the foliation plane with weak maxima parallel to the lineation direction. These patterns are consistent with deformation of omphacite by slip parallel to [001] and suggest glide along (010). The lineation and CPO data reveal a consistent sense of shear zone movement, although the displacement was small. Localized faulting of high-grade rocks accompanied by fluid infiltration can be an important mode of failure in the lower continental crust. Field relations show that granulite facies rocks can exist in a metastable state under eclogite facies conditions and imply that the lower crust can host differing metamorphic facies at the same depth. Deformation of granulite and partial conversion to eclogite, such as is exposed on Holsnøy Island, may be an orogenic-scale process in the lowermost crust of collisional orogens.

Notes: Mineral composition and quantitative thermobarometric studies indicate that the Teslin-Taylor Mountain and Nisutlin terranes within the Teslin suture zone (TSZ), Yukon, record widespread high-P/T metamorphic conditions consistent with subduction zone dynamothermal metamorphism. The highest P–T conditions (575–750° C and 9–17 kbar) are preserved in tectonites formed during normal dip-slip ductile shear. Dextral strike-slip tectonites record lower P–T conditions (400–550° C and 5–8 kbar), and tectonites which show reverse shear have peak temperatures of c. 420° C and a minimum peak pressure of 3 kbar. Dynamothermal metamorphism took place in a west-dipping B-type subduction zone outboard of western North America in Permo-Triassic time. TSZ tectonites were underplated against the hangingwall plate of the subduction zone. Following subduction of the ocean basin which separated North America from the hangingwall plate, TSZ tectonites were overthrust eastward as a coherent structural package as a result of A-type subduction of Cassiar strata in early Jurassic time.(Par)autochthonous Cassiar tectonites, which comprised the leading edge of the western North American margin, record prograde moderate-P, high-T metamorphism (550–750° C and 7–13 kbar) synchronous with top-to-the-east ductile shear. Metamorphism occurred as a result of subduction of the North American margin into the TSZ subduction zone in early Jurassic time. Following metamorphism Cassiar tectonites cooled slowly from 500 to 300° C during the period middle Jurassic to middle Cretaceous.TSZ and Cassiar tectonites were deformed during changing P–T conditions. Data from each of these tectonite packages indicate that grain-scale strain partitioning may have allowed local recrystallization of individual minerals by the addition of mechanical energy. The composition of the new grains reflects the P–T conditions under which that particular grain was deformed.

Notes: The Al–Mg-rich granulites from the In Ouzzal craton, Algeria, show a great diversity of mineral reactions which correspond to continuous equilibria as predicted by phase relationships in the FeO–MgO–Al2O3–SiO2 system. The sequence of mineral reactions can be subdivided into three distinct stages: (1) a high-P stage characterized by the growth of coarse mineral assemblages involving sapphirine and the disappearance of early corundum and spinel-bearing assemblages; (2) a high-T stage characterized by the development of Sa–Qz-bearing assemblages; and (3) a later stage, in which garnet-bearing assemblages are replaced by more or less fine symplectites involving cordierite.During the course of early mineral reactions, the distribution coefficient, Kd, between the various ferromagnesian phases decreased significantly whereas Al2O3 in pyroxene increased concomitantly. These observations, when combined with topological constraints, clearly indicate that the high-P stage 1 was accompanied by a significant rise in temperature (estimated at 150 ± 50° C) under near isobaric conditions, in agreement with the reaction textures. By stage 2, pressure and temperature were extreme as evidenced by the low Kd value between orthopyroxene and garnet (Kd= 2.06–1.99), the high alumina content in pyroxene (up to 11.8%) and the high magnesium content in garnet [100 Mg/(Mg + Fe) = 60.6]. Mineral thermometry based on Fe–Mg exchange between garnet and pyroxene and on Al-solubility in pyroxene gives temperatures close to 970 ± 70° C at 10 ± 1.5 kbar. These results are in agreement with the development of Sa–Qz assemblages on a local scale.Late mineral reactions have been produced during a decompression stage from about 9 to 6 kbar. Except for local re-equilibration of Mg and Fe at grain boundaries, there is no evidence for further reactions below 700° C.We interpreted the whole set of mineral reactions as due to changes in pressure and temperature during a tectonic episode located at c. 2 Ga. Because of the lack of evidence for further uplift after the thermal relaxation which occurred at c. 6 kbar, it is possible however that the exhumation of this granulitic terrane occurred in a later tectonic event unrelated to its formation.

Notes: Investigation of microstructural relationships in major movement zones in metamorphic rocks, where the sense of displacement is known from regional geological relationships, indicates numerous problems with current concepts of shear-sense criteria and their application. The direction of apparent shearing commonly conflicts from one criterion to another (e.g. from the symmetry of quartz c-axis orientation diagrams to the asymmetry of extensional crenulation cleavages). This implies that interpretations of shear sense along foliations from some mesoscale and microscale criteria have been erroneous.A new approach to interpreting shear sense, involving the use of strain fields, resolves conflicts in mesoscopic and microscopic criteria and provides a method for determining coherent shear-sense histories extending back before the last shearing event for ‘any foliated metamorphic rock’. It also provides a powerful tool for determining the structural/metamorphic path that a rock has followed within an orogen. For determination of the shear sense on the last foliation developed in a rock, this approach uses geometries developed around competent heterogeneities such as quartz pebbles, pegmatite pods, veins, porphyroclasts, porphyroblasts and breccia clasts. A shear-sense history is derived by applying this approach to earlier foliations preserved within the heterogeneities and their strain shadows.

Notes: Graphitization and coarsening of organic material in carbonate-bearing metasedimentary rocks is accompanied by carbon isotope exchange which is the basis of a refractory, pressure-independent geothermometer. Comparison of observed isotopic fractionations between calcite and graphite (δ13CCal–Gr) with independent petrological thermometers provides the following empirical calibration over the range 400–800°C: δ13CCal–Gr= 5.81 times 106×T–2(K) - 2.61. This system has its greatest potential in marbles where calcite + graphite is a common assemblage and other geothermometers are often unavailable. The temperature dependency of this empirical calibration differs from theoretical calibrations; reasons for this are unclear but the new empirical calibration yields temperature estimates in better agreement with independent thermometry from several terranes and is preferred for geological applications.Both calcite-graphite isotopic thermometry and calcite-dolomite solvus thermometry are applied to marble adjacent to the Tudor gabbro in the Grenville Province of Ontario, Canada. The marble has undergone two metamorphic episodes, early contact metamorphism and later regional metamorphism. Values of δ13CCal–Gr decrease regularly from c. 8‰ in samples over 2 km from the pluton to values of 3–4‰ within 200 m of the contact. These samples appear to preserve fractionations from the early thermal aureole with the empirical geothermometer, and indicate temperatures of 450–500° C away from the intrusion and 700–750°C near the gabbro. This thermal profile around the gabbro is consistent with conductive heat flow models. In contrast, the distribution of Mg between calcite and dolomite has been completely reset during later regional metamorphism and yields uniform temperatures of c. 500°C, even at the contact.Graphite textures are important for interpreting the results of the calcite–graphite thermometer. Coarsening of graphite approaching the Tudor gabbro correlates with the decrease in isotopic fractionations and provides textural evidence that graphite crystallization took place at the time of intrusion. In contrast to isotopic exchange during prograde metamorphism, which is facilitated by graphitization, retrogressive carbon isotopic exchange appears to require recrystallization of graphite which is sluggish and easily recognized texturally. Resistance of the calcite–graphite system to resetting permits thermometry in polymetamorphic settings to see through later events that have disturbed other systems.

Notes: In order to illustrate different applications of the amphibole-albite-chlorite-epidote-quartz geothermobarometer, pressure-temperature-time (P–T–t) ± space (P–T–t–s) ± deformation (P–T–t–d) paths have been established from literature data. They are discussed as a function of the chemical, equilibrium and microstructural data available in each case, and compared with the conclusions already established by other methods. It is clear that it is necessary to know the relative chronology of the events (directions of zoning of minerals in successive microstructural positions) to establish precise P–T paths; this enables reconstruction of complex geodynamic histories. From this point of view, it is necessary to analyse the maximum possible number of minerals in a few well-chosen metabasic rocks showing different generations of blastesis. The rocks should belong to different tectonic units to obtain the best overall picture of a metamorphic complex.

Notes: In a number of recent papers, the theory has been postulated that porphyroblasts as a rule do not rotate with respect to geographical coordinates, and can be used to determine the original orientation of older foliations. Complex inclusion patterns in spiral garnets have even been used to advocate a new model of orogenesis, involving several alternating phases of horizontal shortening and extension. Critical assessment of the assumptions and data used to support the theory of irrotational porphyroblasts reveals numerous flaws. Millipede structures, used as proof for flow partitioning, can also form by other flow geometries. Evidence quoted to support irrotational behaviour of porphyroblasts is unsound. Porphyroblasts do occur in sets with a preferred orientation of the internal foliation trace, but these cannot be shown to represent original orientations. Microstructures which resemble truncation planes in spiral garnets are used as evidence that these structures developed by several phases of deformation and as proof for periodic extension and horizontal shortening in orogenesis. They can, however, also be explained by intermittent growth of a rotating porphyroblast during a single phase of deformation. Finally, porphyroblast sets in which orientation is a function of aspect ratio indicate that porphyroblast rotation with respect to kinematic axes does occur in at least some situations.

Notes: In the Hazeldene area, situated in the Mount Isa Inlier, Queensland, the metamorphic grade changes from chlorite zone, through biotite and cordierite zones, to sillimanite/K-feldspar zone.Microstructural studies of rocks near the sillimanite isograd demonstrate that cordierite grew early during the development of a steep foliation (S2), was replaced by biotite, andalusite and sillimanite at the metamorphic peak late in S2, and in turn by kyanite + chlorite adjacent to localized small post-D2 shear zones. Although the anticlockwise P–T–t path is well defined, the precise P–T conditions are uncertain because of problems with experimental and thermodynamic data. The best estimate for the metamorphic peak for rocks close to the sillimanite isograd is around 600° C at 4 kbar.The metamorphism has been dated at 1544 Ma, and was synchronous with a major crustal shortening event. Because proposed extensional events occurred more than 60 Ma earlier, their contribution to the peak metamorphic thermal perturbation would have been insignificant. The syn-metamorphic Mica Creek Pegmatites, the abundance of high heat-producing elements in the nearby pre-D2 Sybella Granite, and advective heat by fluids which caused considerable metasomatism in the Hazeldene area, may have each contributed to the thermal budget. However, the metamorphic thermal gradient may be 80°C km-1 or higher, strongly suggesting a local magmatic control. As none are known in the area, such syn-metamorphic plutons would have to lie beneath the exposed high-grade rocks.

Notes: Charnockitic alteration (arrested orthopyroxene formation in biotite- and amphibole-bearing rocks) occurs in high-grade terranes of all ages. Three criteria are used to show that this alteration was produced in many locations by a migrating fluid phase: (i) diffuseness of the alteration—the alteration zones are often quite unlike discrete migmatitic veins; (ii) relation to deformation—most occurrences show alteration closely associated with warping of foliation or dilation cracks; (iii) open-system alteration—whilst some occurrences represent nearly isochemical alteration, slight changes in bulk composition, often loss of mafic constituents and gain of Na and Si, are evident in detailed mass-balance analysis. Y and sometimes Rb are characteristically depleted. Partial melting sometimes accompanied volatile infiltration, as evidenced by more discrete veins and euhedral orthopyroxene. It is quite unlikely, however, that open-system alteration was produced by escape of viscous quartzo-feldspathic melts. Pervasive migration of low-T lamprophyric (mafic–alkaline, CO2-charged) interstitial liquids is a possibility by virtue of their extreme fluidity, but CO2 infiltration was needed to generate these liquids. Vapour-deficient dehydration melting is another feasible mechanism of orthopyroxene formation which may have operated in conjunction with CO2 infiltration.Characteristic development of charnockitic alteration in some prograde amphibolite to granulite facies transitions, as in the Dharwar Craton of South India, suggests that the alteration is a fundamental feature of the granulite facies metamorphism, implying active and causal participation of migrating fluids. In other high-grade terranes like the Adirondack Mountains of New York, this kind of alteration is rare, and fluid action does not seem to have been important in the metamorphism.A vapour phase participating in charnockitic metamorphism was necessarily one of relatively low H2O, therefore presumably rich in CO2. Consideration of possible large CO2 sources leads to the conclusion that emanations from volatile-rich basalts emplaced in the lower crust are the most probable source of charnockitizing fluids. The ultimate source would therefore be enriched subcontinental lithosphere or asthenosphere. The Rb-depleted pyroxene gneiss (charnockitic) terranes may be characteristic of zones of large-scale transcurrent or oblique-motion faults which tap such great depths.

Notes: Quartz veins are developed in a wide range of metasediment types in the upper amphibolite facies rocks of Connemara, and attest to considerable migration of silica. Contrary to common assumptions, there is clear evidence that these veins do not primarily result from movement of fluid to regions of lower P–T down the regional geothermal gradient. Under amphibolite facies conditions, a dilute chloride fluid moving down temperature has the potential to alter 60g of plagioclase to muscovite for each gram of vein quartz precipitated, while cooling over the temperature interval from 650 to 500° C. The absence of significant metasomatic effects in the vein walls effectively precludes a simple origin from such through-flowing, externally derived fluids. The oxygen isotopic composition of matrix quartz shows considerable differences between different rock types (quartzite, pelite and marble), with a range of δ18OSMOW from c.+ 11.5% (quartzite) to + 18.5% (marble). In each rock type, vein quartz compositions closely match those of the matrix quartz. These results demonstrate the importance of local segregation processes in the formation of veins, and suggest that fluid convection cells were not developed during metamorphism on a scale larger than the individual sedimentary formations, if at all.Both oxygen isotope data and the absence of metasomatism indicate that veins form primarily by segregation of quartz from the host lithologies, with only a relatively minor component of through flow of externally derived fluid. Veins are clearly not the major pathways of metamorphic dewatering.It is proposed that abundant veins in the predominantly pelitic Ballynakill Formation formed during peak metamorphic D3 folding because the formation was embrittled by high fluid pressures but was capped by impermeable marble. Hence the pelitic formation fractured repeatedly and the pore fluid drained through the fractures to form veins, while irreversible loss through the rest of the succession was a much less important process.In the central mountains of Connemara, rather pure, unreactive quartzites are cut by widely spaced, laterally extensive quartz veins that are axial planar to D3 folds. These veins may mark pathways whereby metamorphic fluid made its way through the massive impermeable quartzite from lower parts of the nappe pile, but here too, oxygen isotope data indicate considerable segregation of locally derived quartz, reflecting the importance of pumping of fluid between wail rocks and fractures relative to the component of through flow.

Notes: Geothermometry and geobarometry of 10 garnet–oligoclase zone schists in the Franz Josef–Fox Glacier area, Southern Alps, New Zealand, give temperatures ranging from 415 to 625°C and pressures from 5.2 to 9.2 kbar, indicating a T–P array of about 50°C/kbar and inferred peak temperature conditions over a c. 15-km-thick section at depths between c. 20 and 34 km. The present-day distribution of the schist samples implies that only about one-third of the original crustal section is now exposed.The garnet–oligoclase zone schists represent the deeper part of a metamorphosed and deformed accretionary complex that was associated with late Palaeozoic–early Mesozoic subduction along the Gondwana continental margin. Partial uplift (c. 0.2 m/Ma) and erosion of the complex during Jurassic–Cretaceous times (Rangitata uplift) was synchronous with D2 deformation and recrystallization, as recorded by the P–T array. Cenozoic (Kaikoura) uplift and exhumation of the schist since c. 30 Ma to form the Southern Alps was associated with oblique-slip movement on the Alpine Fault. The present-day position and steep eastward dip of isograds and D2 structures suggest considerable clockwise rotation during uplift associated with ductile attenuation and tectonic thinning by over two-thirds of the original schist sequence, largely due to simple shear along schistosity planes. As the schist generally shows only incipient greenschist facies retrograde recrystallization, an apparently complete (although contracted) prograde mineral sequence has been preserved by rapid uplift (〉5 km/Ma) of hot rock and the effects of limited shear heating near the Alpine Fault.

Notes: The Smartville Complex is a late Jurassic, rifted volcanic arc in the northern Sierra Nevada, California. Near Auburn, California, it consists of a lower volcanic unit, dominated by basaltic flows, and an upper volcanic unit of andesitic volcaniclastic rocks, both of which have been intruded by dykes and irregular bodies of diabase. These rocks contain relict igneous minerals, and the metamorphic minerals albite, chlorite, quartz, pumpellyite, prehnite, epidote, amphibole, titanite, garnet, biotite, K-feldspar, white mica, calcite, and sulphide and oxide minerals.Prehnite–pumpellyite (PrP), prehnite–actinolite (PrA), and greenschist (GS) zones have been identified. The pumpellyite-out isograd separates the PrP and PrA zones, and the prehnite-out isograd separates the PrA and GS zones. The minerals Ab + Qtz + Mt + Tn are common to most assemblages in all three zones. The MgO/(MgO + FeO) ratio of the effective bulk composition has an important and systematic effect on the observed mineral assemblages in the PrP zone. Prehnite-bearing assemblages contain the additional minerals, Pmp + Amp + Ep + Chl in MgO-rich rocks, and either Pmp + Ep + Chl or Amp + Ep + Chl in less magnesian rocks. Subcalcic to calcic amphibole is common in the PrP zone. The mineral assemblage Prh + Act + Ep + Chl, without Pmp, characterizes the PrA zone, and the mineral assemblage Act + Ep + Chl, without Prh or Pmp, characterizes the GS zone. The disappearance of pumpellyite and prehnite occurred by continuous reactions.The sequence of mineral assemblages was produced by burial metamorphism at P–T conditions of 300° 50°C at approximately 2.5 ± 0.5 kbar. During metamorphism, the composition of the fluid phase was nearly 100% H2O and the oxygen fugacity was between the hematite–magnetite and quartz–fayalite–magnetite buffers.

Notes: Garnet-bearing mineral assemblages are commonly observed in pelitic schists regionally metamorphosed to upper greenschist and amphibolite facies conditions. Modelling of thermodynamic data for minerals in the system Na2O–K2O–FeO–MgO–Al2O3–SiO2–H2O, however, predicts that garnet should be observed only in rocks of a narrow range of very high Fe/Mg bulk compositions. Traditionally, the nearly ubiquitous presence of garnet in medium- to high-grade pelitic schists is attributed qualitatively to the stabilizing effect of MnO, based on the observed strong partitioning of MnO into garnet relative to other minerals. In order to quantify the dependence of garnet stability on whole-rock MnO content, we have calculated mineral stabilities for pelitic rocks in the system MnO–Na2O–K2O–FeO–MgO–Al2O3–SiO2–H2O for a moderate range of MnO contents from a set of non-linear equations that specify mass balance and chemical equilibrium among minerals and fluid. The model pelitic system includes quartz, muscovite. albite, pyrophyllite, chlorite, chloritoid, biotite, garnet, staurolite, cordierite, andalusite, kyanite. sillimanite, K-feldspar and H2O fluid. In the MnO-free system, garnet is restricted to high Fe/Mg bulk compositions, and commonly observed mineral assemblages such as garnet–chlorite and garnet–kyanite are not predicted at any pressure and temperature. In bulk compositions with XMn= Mn/(Fe + Mg + Mn) 〉 0.01, however, the predicted garnet-bearing mineral assemblages are the same as the sequence of prograde mineral assemblages typically observed in regional metamorphic terranes. Temperatures predicted for the first appearance of garnet in model pelitic schist are also strongly dependent on whole-rock MnO content. The small MnO contents of normal pelitic schists (XMn= 0.01–0.04) are both sufficient and necessary to account for the observed stability of garnet.

Notes: The Velay dome (French Massif Central) offers a quasi-continuous section across an anatectic domain comprising low- to high-grade schists, gneisses and granites. Two main tectonometamorphic events, and their related generation of granitic material, were recognized in addition to a major Barrrovian tangential event (D2) attributed to intracontinental collision tectonics: (i) a medium- to low-P, high-T event (D3) which gave rise to migmatites and syntectonic monzonitic granites and granodiorites, and (ii) a widespread melting event (D4) which led to the generation of migmatities, the Velay granite and post-anatectic granites.Thermobarometry on samples collected from both the metamorphic envelope and the granitic core distinguishes two distinct geotherms: (i) a first, associated with the D3 event, characterized by P 〉 5 kbar, T≤ 750° C and water-present melting (biotite remains stable) which led to large-scale migmatization but minor amount of granites; (ii) a second, associated with the D4 event and characterized by vapour-absent melting (P= 4–5 kbar, T= 760–850° C) which gave rise to the Velay granites and late-migmatitic granites. The temperature increase during the D4 event is attributed to the intrusion of hot mafic magmas within the crust.The time-integrated features of the different granitic rocks in the Velay dome can be directly related to aH2O in the source region and illustrate the progressive dehydration of a middle to lower crustal segment over 60 Ma.

Notes: Metabasaltic rocks in the Klamath Mountains of California with ‘komatiitic’ major element concentrations were investigated in order to elucidate the origin of the magnesian signature. Trace-element concentrations preserve relict igneous trends and suggest that the rocks are not komatitic basalts, but immature arc rocks and within-plate alkalic lavas. Correlation of ‘excess’ MgO with the volume per cent hornblende (±clinopyroxene) suggests that the presence of cumulus phases contributes to the MgO-rich compositions. Early submarine alteration produced regional δ18O values of +10±1.5%° and shifts in Al2O3, Na2O, and K2O concentrations. Regional metamorphic grade in the study area varies from biotite-zone greenschist facies (350–550°C, c. 3 kbar) southward to prehnite–actinolite facies (200–400°C, ≤3 kbar), but little isotopic or elemental change occurred during the regional recrystallization. The greenschist facies assemblage is actinolitic hornblende + phengite + epidote + sodic plagioclase + microcline + chlorite + titanite + hematite + quartz in Ti-poor metabasaltic rocks; in addition to these phases biotite is present in Ti-rich analogues. Lower grade greenstones contain prehnite and more nearly stoichiometric actinolite. The moderate to low pressures of regional metamorphism are compatible with P–T conditions in a magmatic arc. Later contact metamorphism at 2–2.9±0.5 kbar and at peak temperatures approaching 600° C around the English Peak and Russian Peak granodiorites produced 3–4–km-wide aureoles typified by gradual, systematic increases in the pargasite content of amphibole, muscovite content of potassic white mica, and anorthite content of plagioclase compositions. Metasomatism during contact metamorphism produced further increases in bulk-rock δ18OSMOW of as much as +6%°. Thus, the unusually MgO-rich nature of the Sawyers Bar rocks may be attributed at least partly to metasomatism and the presence of magnesian cumulus phases.

Notes: Because many durable goods industries are also industries with significant innovative activity, we analyze the relationship between an oligopolistic firm's choice of product durability and cost-reducing innovation in rental and sales markets. We demonstrate the firm's durability may be greater than, equal to, or less than a social planner's, depending upon innovation's impact on the marginal cost of durability. We also show that the firm's innovation level can be maximized for moderate levels of concentration if the output is durable. This provides theoretical support for the empirical observation that innovation is often maximized in markets of moderate concentration.

Notes: In this paper we attempt to disentangle the effects of deregulation on rail costs from those directly attributable to mergers. We estimate that cost reductions obtained from mergers ranged from a high of 33% for the Burlington Northern to a low of a 3% cost increase for the CSX. However, firms not engaged in significant merger activities experienced similar cost differentials indicating that consolidation was not a prerequisite for cost savings. We conclude that although mergers did confer some benefits on the participating firms, they were not a prerequisite for railroads being able to achieve substantial cost savings.

Notes: For a homogeneous product oligopoly market, possibilities for pure strategy Nash equilibria in prices are studied. Consumers, who each nonstrategically purchase one unit up to a common reservation price, are hypothesized to be more concerned with large price differences (and therefore buy from the cheapest firm) than slightly different prices. For the duopoly case, existence, uniqueness, and characterization results are provided. Linear examples are given with 2 and n firms.

Notes: This paper analyzes the ongoing transformation of the soft drink distribution of Coca-Cola and Pepsi-Cola from systems of independent bottlers to captive bottling subsidiaries, A transaction cost-based theory is developed to explain this restructuring. It is postulated that changes in the external environment and the resulting changes in the strategies of Coca-Cola and Pepsi-Cola raised the costs of transacting between them and their independent bottlers. Two types of empirical tests are presented. One exploits the difference in the distribution of Coca-Cola and Pepsi-Cola in the fountain channel. The other consists of statistical analyses of the competitive effects of the move toward captive distribution. Both types of tests support the basic hypothesis.

Notes: Food-related illness in the United States is estimated to affect over six million people per year and cost the economy several billion dollars. These illnesses and costs could be reduced if minimum infectious doses were established and used as the basis of regulations and monitoring. However, standard methodologies for dose-response assessment are not yet formulated for microbial risk assessment. The objective of this study was to compare dose response models for food-borne pathogens and determine which models were most appropriate for a range of pathogens. The statistical models proposed in the literature and chosen for comparison purposes were log-normal,(19) log-logistic,(17) exponential,(7,9,17)ß-Poisson(7,9,18) and Welbull-Gamma.(3) These were fit to four data sets also taken from published literature, Shigella fiexneri,(9,12,13)Shigella dysenteriae,(9,11)Campylobacter jejuni,(15,16) and Salmonella typhosa,(7,14) usingthe method of maximum likelihood. The Weibull-gamma, the only model with three parameters, was also the only model capable of fitting all the data sets examined using the maximum likelihood estimation for comparisons. Infectious doses were also calculated using each model. Within any given data set, the infectious dose estimated to affect one percent of the population ranged from one order of magnitude to as much as nine orders of magnitude, illustrating the differences in extrapolation of the dose response models. More data are needed to compare models and examine extrapolation from high to low doses for food-borne pathogens.

Notes: Environmental managers are increasingly charged with involving the public in the development and modification of policies regarding risks to huma health and the environment. Involving the public in environmental decision making first requires a broad understanding of how and why the public perceives various risks. The Savannah River Stakeholder Study was conducted withthe purpose of investigating individual, economic, and social characteristics of risk perceptions among those living near the Savannah River Nuclear Weapons Site. A number of factors were found to impact risk perceptions among those living near the site. One's estimated proximity to the site and relative river location surfaced as strong determinants of risk perceptions among SRS residents. Additionally, living in a quality neighborhood and demonstrating a willingness to accept health risks for economic gain strongly abated heightened risk perceptions.

Notes: Five atmospheric transport models were evaluated for use in Phase I1 ofthe Historical Public Exposures Studies at the Rocky Flats Plant. Models included a simple straight-line Gaussian plume model (ISCST2), several integrated puff models (RATCHET, TRIAD, and INPUFF2), and a complex terrain mode (TRAC). Evaluations were based on how well model predictions compared with sulfur hexafluoride tracer measurements taken in the vicinity of Rocky Flats in February 1991. Twelve separate tracer experiments were conducted, eac lasting 9 hr and measured at 140 samplers in arcs 8 and 16 km from the release point at Rocky Flats. Four modeling objectives were defined based on the endpoints of the overall study: (1) the unpaired maximum hourly average concentration, (2) paired time-averaged concentration, (3) unpaired time-averaged concentration, and (4) arc-integrated concentration. Performance measures were used to evaluate models and focused on the geometric mean and standard deviation of the predicted-to-observed ratio and the correlation coefficient between predicted and observed concentrations. No one model consistently outperformed the others in all modeling objectives and performance measures. About 75% of the maximum hourly concentration predictions were within a factor of 5 of the observations. About 64% of the paired and 80% of the unpaired time-averaged model predictions were within a factor of 5 of theobservations. The overall performance of the RATCHET model was somewhat better than the other models. All models appeared to experience difficulty defining plume trajectories, which was attributed to the influence of multilayered flow initiated by terrain complexities and the diurnal flow patterns characteristic of the Colorado Front Range.

Notes: Calculation of accident dose-risk estimates with the RADTRAN code requires input data describing the population likely to be affected by the plumeof radioactive material (RAM) released in a hypothetical transportation accident. In the existing model, population densities within 1/2 mile (0.8 km)of the route centerline are tabulated in three ranges (Rural, Suburban, andUrban). These population densities may be of questionable validity since the plume in the RADTRAN analysis is assumed to extend out to 120 km from thehypothetical accident site. We present a GIs-based population model which accounts for the actual distribution of population under a potential plume, and compare accident-risk estimates based on the resulting population densities with those based on the existing model. Results for individual points along a route differ greatly, but the cumulative accident risks for a sample route of a few hundred kilometers are found to be comparable, if not identical. We conclude, therefore, that for estimation of aggregate accident risks over typical routes of several hundred kilometers, the existing, simpler RADTRAN model is sufficiently detailed and accurate.

Notes: This paper discusses a number of the key challenges to the acceptance and aplication of probabilistic risk analysis (PRA). Those challenges include:(a) the extensive reliance on subjective judgment in PRA, requiring the development of guidance for the use of PRA in risk-informed regulation, and possibly the development of “robust” or “reference” prior distributions to minimize the reliance on judgment; and (b) the treatment of human perfor-mance in PRA, including not only human errorperse but also management and organizational factors more broadly. All of these areas are seen as presentinginteresting research challenges at the interface between engineering and other disciplines.

Notes: Rational decision making requires that the total uncertainty about a variate of interest (a predictand) be quantified in terms of a probability distribution, conditional on all available information and knowledge. Supposethe state-of-knowledge is embodied in a deterministic model, which is imperfect and outputs only an estimate of the predictand. Fundamentals are presented of two Bayesian methods for producing a probabilistic forecast via anydeterministic model. The Bayesian Processor of Forecast (BPF) quantifies the total uncertainty in terms of a posterior distribution, conditional on model output. The Bayesian Forecasting System (BFS) decomposes the total uncertainty into input uncertainty and model uncertainty, which are characterized independently and then integrated into a predictive distribution. The BFS is compared with Monte Carlo simulation and “ensemble forecasting” technique, none of which can alone produce a probabilistic forecast that quantifies the total uncertainty, but each can serve as a component of the BFS.

Notes: This article describes an application of a method for assessing risks associated with the introduction of an organism into a new environment. The test organism was a binucleate Rhizoctonia fungal isolate that has potential for commercial development as a biological control agent for damping-off diseases in bedding plants. A test sample of host plant species was selected using the centrifugal phylogenetic host range principles, but with an emphasis on economic species. The effect of the fungus on the plant was measured for each species and expressed on a logarithmic scale. The effects on weights of shoots and roots per container were not normally distributed, nor were the effects on the number of plants standing (those which survived). Statements about the effect on the number standing and the shoot weight per container involved using the observed (empirical) distribution. This is illustrated with an example. Problems were encountered in defining the population of species at risk, and in deciding how this population should be formally sampled. The limitations of the method are discussed.

Notes: The heuristic-systematic information processing model (HSM) holds that individuals will use one or both of these modes of information processing when attempting to evaluate information in order to arrive at a judgment. Systematic processing is defined by effortful scrutiny and comparison of information, whereas heuristic processing is defined by the use of cues to amve more easily at a judgment. Antecedents to the two processing modes include information sufficiency, motivation, and self-efficacy. Structural equation modeling is used to examine competing configuration of this model and to evaluate the model as appropriate for predicting risk judgment. The model also is evaluated across three groups that vary with respect to their level of concern. These analyses are executed within a case study involving an epidemiological investigation of a suspected cancer cluster. The analysis confirms the HSM's theoretically proposed structure and shows it to be a useful vehicle for evaluating risk judgment. In the overall analysis, antecedent variables generally function as specified by theory. Systematic processing is predicted by greater motivation. Heuristic processing is predicted by information sufficiency. Self-efficacy is a significant predictor of both processing modes. And heuristic processing is shown to be associated with judgment of less risk. However, when the analysis is contrasted across three groups (those concerned about cancer, not concerned and uncertain) it is shown that the model is significantly more robust for the uncertain group. This finding may have implications for the use of the HSM in risk research specifically, and in field research generally.

Notes: South Carolina has issued fish consumption advisories for the Savannah River based on mercury and radionuclide levels. We examine differences in fishing rates and fish consumption of 258 people interviewed while fishing along the Savannah River, as a function of age, education, ethnicity, employment history, and income, and test the assumption that the average consumption of fish is less than the recreational value of 19 kg/year assumed by risk assessors. Ethnicity and education contributed significantly to explaining variations in number of fish meals per month, serving size, and total quantity of fish consumed per year. Blacks fished more often, ate more fish meals of slightly larger serving sizes, and consumed more fish per year than did Whites. Although education and income were correlated, education contributed most significantly to behavior; people who did not graduate from high school ate fish more often, ate more fish per year, and ate more whole fish than people who graduated from high school. Computing consumption of fish for each person individually indicates that (1) people who eat fish more often also eat larger portions, (2) a substantial number of people consume more than the amount of fish used to compute risk to recreational fishermen, (3) some people consume more than the subsistence level default assumption (50 kglyear) and (4) Blacks consume more fish per year than Whites, putting them at greater risk from contaminants in fish. Overall, ethnicity, age, and education contributed to variations in fishing behavior and consumption.

Notes: The association between daily fluctuations in ambient particulate matter and daily variations in nonaccidental mortality have been extensively investigated. Although it is now widely recognized that such an association exists, the form of the concentration-response model is still in question. Linear, no threshold and linear threshold models have been most commonly examined. In this paper we considered methods to detect and estimate threshold concentrations using time series data of daily mortality rates and air pollution concentrations. Because exposure is measured with error, we also considered the influence of measurement error in distinguishing between these two completing model specifications. The methods were illustrated on a 15-year daily time series of nonaccidental mortality and particulate air pollution data in Toronto, Canada. Nonparametric smoothed representations of the association between mortality and air pollution were adequate to graphically distinguish between these two forms. Weighted nonlinear regression methods for relative risk models were adequate to give nearly unbiased estimates of threshold concentrations even under conditions of extreme exposure measurement error. The uncertainty in the threshold estimates increased with the degree of exposure error. Regression models incorporating threshold concentrations could be clearly distinguished from linear relative risk models in the presence of exposure measurement error. The assumption of a linear model given that a threshold model was the correct form usually resulted in overestimates in the number of averted premature deaths, except for low threshold concentrations and large measurement error.

Notes: Facility-specific information on pollution was obtained for 36 coke plants and 46 oil refineries in the United States and matched with information on populations surrounding these 82 facilities. These data were analyzed to determine whether environmental inequities were present, whether they were more economic or racial in nature, and whether the racial composition of nearby communities has changed significantly since plants began operations. The Census tracts near coke plants have a disproportionate share of poor and nonwhite residents. Multivariate analyses suggest that existing inequities are primarily economic in nature. The findings for oil refineries are not strongly supportive of the environmental inequity hypothesis. Rank ordering of facilities by race, poverty, and pollution produces limited (although not consistent) evidence that the more risky facilities tend to be operating in communities with above-median proportions of nonwhite residents (near coke plants) and Hispanic residents (near oil refineries). Over time, the racial makeup of many communities near facilities has changed significantly, particularly in the case of coke plants sited in the early 1900s. Further risk-oriented studies of multiple manufacturing facilities in various industrial sectors of the economy are recommended.

Notes: We investigated the way results of human health risk assessments are used, and the theory used to describe those methods, sometimes called the “NAS paradigm.” Contrary to a key tenet of that theory, current methods have strictly limited utility. The characterizations now considered standard, Safety Indices such as “Acceptable Daily Intake,”“Reference Dose,” and so on, usefully inform only decisions that require a choice between two policy alternatives (e.g., approve a food additive or not), decided solely on the basis of a finding of safety. Riskis characterized as the quotient of one of these Safety Indices divided by an estimate of exposure: a quotient greater than one implies that the situation may be considered safe. Such decisions are very widespread, both in the U. S. federal government and elsewhere. No current method is universal; different policies lead to different practices, for example, in California's “Proposition 65,” where statutory provisions specify some practices. Further, an important kind of human health risk assessment is not recognized by this theory: this kind characterizes risk as likelihood of harm, given estimates of exposure consequent to various decision choices. Likelihood estimates are necessary whenever decision makers have many possible decision choices and must weigh more than two societal values, such as in EPA's implementation of “conventional air pollutants.” These estimates can not be derived using current methods; different methods are needed. Our analysis suggests changes needed in both the theory and practice of human health risk assessment, and how what is done is depicted.

Notes: The prominent role of animal bioassay evidence in environmental regulatory decisions compels a careful characterization of extrapolation uncertainties. In noncancer risk assessment, uncertainty factors are incorporated to account for each of several extrapolations required to convert a bioassay outcome into a putative subthreshold dose for humans. Measures of relative toxicity taken between different dosing regimens, different endpoints, or different species serve as a reference for establishing the uncertainty factors. Ratios of no observed adverse effect levels (NOAELs) have been used for this purpose; statistical summaries of such ratios across sets of chemicals are widely used to guide the setting of uncertainty factors. Given the poor statistical properties of NOAELs, the informativeness of these summary statistics is open to question. To evaluate this, we develop an approach to “calibrate” the ability of NOAEL ratios to reveal true properties of a specified distribution for relative toxicity. A priority of this analysis is to account for dependencies of NOAEL ratios on experimental design and other exogenous factors. Our analysis of NOAEL ratio summary statistics finds (1) that such dependencies are complex and produce pronounced systematic errors and (2) that sampling error associated with typical sample sizes (50 chemicals) is nonnegligible. These uncertainties strongly suggest that NOAEL ratio summary statistics cannot be taken at face value; conclusions based on such ratios reported in well over a dozen published papers should be reconsidered.

Notes: We compare the regulatory implications of applying the traditional (linearized) and exact two-stage dose-response models to animal carcinogenic data. We analyze dose-response data from six studies, representing five different substances, and we determine the “goodness- of-fit” of each model as well as the 95% confidence lower limit ofthe dose corresponding to a target excess risk of 10-5 (the target risk doseTRD). For the two concave datasets, we find that the exact model gives a substantially better fit to the data than the traditional model, and that the exact model gives aTRD that is an order of magnitude lower than that given by the traditional model. In the other cases, the exact model gives a fit equivalent to or better than the traditional model. We also show that although the exact two-stage model may exhibit dose-response concavity at moderate dose levels, it is always linear or sublinear, and never supralinear, in the low-dose limit. Because regulatory concern is almost always confined to the low-dose region extrapolation, supralinear behavior seems not to be of regulatory concern in the exact two-stage model. Finally, we find that when performing this low-dose extrapolation in cases of dose-response concavity, extrapolating the model fit leads to a more conservative TRD than taking a linear extrapolation from 10% excess risk. We conclude with a set of recommendations.

Notes: The characterization and treatment of uncertainty poses special challenges when modeling indeterminate or complex coupled systems such as those involved in the interactions between human activity, climate and the ecosystem. Uncertainty about model structure may become as, or more important than, uncertainty about parameter values. When uncertainty grows so large that prediction or optimization no longer makes sense, it may still be possible to use the model as a “behavioral test bed” to examine the relative robustness of alternative observational and behavioral strategies. When modelsmust be run into portions of their phase space that are not well understood, different submodels may become unreliable at different rates. A common example involves running a time stepped model far into the future. Several strategies can be used to deal with such situations. The probability of model failure can be reported as a function of time. Possible alternative “surprises” can be assigned probabilities, modeled separately, and combined. Finally, through the use of subjective judgments, one may be able to combine, and over time shift between models, moving from more detailed to progressively simpler order-of-magnitude models, and perhaps ultimately, on to simple bounding analysis.

Notes: In any model the values of estimates for various parameters are obtained from different sources each with its own level of uncertainty. When the probability distributions of the estimates are obtained as opposed to point values only, the measurement uncertainties in the parameter estimates may be addressed. However, the sources used for obtaining the data and the models used to select appropriate distributions are of differing degrees of uncertainty. A hierarchy of different sources of uncertainty based upon one's ability to validate data and models empirically is presented. When model parameters are aggregated with different levels of the hierarchy represented, this implies distortion or degradation in the utility and validity of the models used. Means to identify and deal with such heterogeneous data sourcesare explored, and a number of approaches to addressing this problem is presented. One approach, using Range/confidence Estimates coupled with an Information Value Analysis Process, is presented as an example.

Notes: Variability arises due to differences in the value of a quantity among different members of a population. Uncertainty arises due to lack of knowledge regarding the true value of a quantity for a given member of a population. We describe and evaluate two methods for quantifying both variability and uncertainty. These methods, bootstrapsimulation and a likelihood-based method, are applied to three datasets. The datasetsinclude a synthetic sample of 19 values from a Lognormal distribution, a sample of nine values obtained from measurements of the PCB concentration in leafy produce, and asample of five values for the partitioning of chromium in the flue gas desulfurization system of coal-fired power plants. For each of these datasets, we employ the two methods to characterize uncertainty in the arithmetic mean and standard deviation, cumulative distribution functions based upon fitted parametric distributions, the 95th percentile of variability, and the 63rd percentile of uncertainty for the 81st percentile of variability. The latter is intended to show that it is possible to describe anypoint within the uncertain frequency distribution by specifying an uncertainty percentile and a Variability percentile. Using the bootstrap method, we compare results based upon use of the method of matching moments and the method of maximum likelihood for fitting distributions to data. Our results indicate that with only 5-19 data pointsas in the datasets we have evaluated, there is substantial uncertainty based upon random sampling error. Both the boostrap and likelihood-based approaches yield comparable uncertainty estimates in most cases.

Notes: The Cordilleran orogen in south-eastern Alaska includes 14 distinct metamorphic belts that make up three major metamorphic complexes, from east to west: the Coast plutonic–metamorphic complex in the Coast Mountains; the Glacier Bay–Chichagof plutonic–metamorphic complex in the central part of the Alexander Archipelago; and the Chugach plutonic–metamorphic complex in the northern outer islands. Each of these complexes is related to a major subduction event. The metamorphic history of the Coast plutonic–metamorphic complex is lengthy and is related to the Late Cretaceous collision of the Alexander and Wrangellia terranes and the Gravina overlap assemblage to the west against the Stikine terrane to the east. The metamorphic history of the Glacier Bay–Chichagof plutonic–metamorphic complex is relatively simple and is related to the roots of a Late Jurassic to late Early Cretaceous island arc. The metamorphic history of the Chugach plutonic–metamorphic complex is complicated and developed during and after the Late Cretaceous collision of the Chugach terrane with the Wrangellia and Alexander terranes.The Coast plutonic–metamorphic complex records both dynamothermal and regional contact metamorphic events related to widespread plutonism within several juxtaposed terranes. Widespread moderate-P/T dynamothermal metamorphism affected most of this complex during the early Late Cretaceous, and local high-P/T metamorphism affected some parts during the middle Late Cretaceous. These events were contemporaneous with low- to moderate-P, high-T metamorphism elsewhere in the complex. Finally, widespread high-P–T conditions affected most of the western part of the complex in a culminating late Late Cretaceous event. The eastern part of the complex contains an older, pre-Late Triassic metamorphic belt that has been locally overprinted by a widespread middle Tertiary thermal event.The Glacier Bay–Chichagof plutonic–metamorphic complex records dominantly regional contact-metamorphic events that affected rocks of the Alexander and Wrangellia terranes. Widespread low-P, high-T assemblages occur adjacent to regionally extensive foliated granitic, dioritic and gabbroic rocks. Two closely related plutonic events are recognized, one of Late Jurassic age and another of late Early and early Late Cretaceous age; the associated metamorphic events are indistinguishable. A small Late Devonian or Early Mississippian dynamothermal belt occurs just north-east of the complex. Two older low-grade regional metamorphic belts on strike with the complex to the south are related to a Cambrian to Ordovician orogeny and to a widespread Middle Silurian to Early Devonian orogeny.The Chugach plutonic–metamorphic complex records a widespread late Late Cretaceous low- to medium/high-P, moderate- T metamorphic event and a local transitional or superposed early Tertiary low-P, high-T regional metamorphic event associated with mesozonal granitic intrusions that affected regionally deformed and metamorphosed rocks of the Chugach terrane. The Chugach complex also includes a post-Late Triassic to pre-Late Jurassic belt with uncertain relations to the younger belts.

Notes: Small pods of silica-undersaturated Al-rich and Mg-rich granulite facies rocks containing sapphirine, pleonastic spinel, kornerupine, cordierite, orthopyroxene, corundum, sillimanite and gedrite are scattered throughout the NE Strangways Range, Central Australia. These are divided into four distinct rock types, namely orthopyroxene-rich aluminous granofels and metapelitic gneisses containing sapphirine, spinel or kornerupine. Two granulite facies metamorphic events are recognized, of which only the first (M1) is considered in this paper.Peak metamorphic mineral parageneses indicate that the M1 thermal maximum occurred at approximately 900–950 °C and 8–9 kbar. All samples are characterized by profuse and diverse coronitic and symplectic reaction textures. These are interpreted as evidence for the sequential crossing of the following reactions in the system FMAS:cordierite + spinel + corundum = sapphirine + sillimanite, cordierite + spinel = orthopyroxene + sapphirine + sillimanite, sapphirine + spinel + sillimanite = orthopyroxene + corundum, sapphirine + sillimanite = cordierite + orthopyroxene + corundum.Phase stability relationships in FMAS and MASH indicate an anticlockwise P–T path terminated by isobaric cooling. Such a path is exemplified by early low-P mineral parageneses containing spinel, corundum and gedrite and the occurrence of both prograde and retrograde corundum. Reaction textures preserve evidence for an increase in aH2O and aB2O3 with progressive isobaric cooling. This hydrous retrogression resulted from crystallization of intimately associated M1 partial melt segregations. There is no evidence for voluminous magmatic accretion giving rise to the high M1 thermal gradient. The M1 P–T path may be the result of either lithospheric thinning after both crustal thickening and burial of the supracrustal terrane, or concomitant crustal thickening and mantle lithosphere thinning.