The relationship between zooplankton,conductivity and lake-water ionic composition in 111 lakes from the Interior Plateau of British Columbia,Canada

Zooplankton collected from vertical net tows were related to the environmental variables from 98 lakes from the Interior Plateau of British Columbia. Canonical correspondence analysis showed that both salinity and ionic composition (pH and Mg) of the lake-water made major and significant contributions to the first two ordination axes (=0.42 and 0.11 respectively,P<0.05). BothArtemia franciscana andMoina hutchinsoni had their highest relative abundance in meso-hypersaline waters. However,Artemia franciscana preferred waters that were higher in Mg and Ca, whileMoina hutchinsoni was found in waters that were lower in Mg and Ca. Similarly, at intermediate salinities,Daphnia pulex and the calanoid copepods preferred waters slightly lower in Mg and Ca, whereasCeriodaphnia laticaudata andSimocephalus spp. were relatively more common in waters higher in Mg and Ca. Because the freshest lakes studied varied much less in ionic composition, the zooplankton in these lakes did not show a preference to ionic composition. As expected, multi-generic groups, such as the calanoid copepods, cyclopoid copepods and nauplii, had wider tolerances to conductivity than groups identified to lower taxonomic levels. Significant weighted-averaging regression and calibration models of conductivity were developed based on zooplankton species composition from the study lakes (r ²=0.56,P<0.05). Samples composed largely of multi-generic taxa yielded the worst estimates of salinity in the reconstruction model. This study suggests that zooplankton community composition may be developed into a useful proxy for paleosalinity reconstruction.     

The spatial mobility of two generations of young adults in Norway

"This article presents results from a longitudinal study using data from the Norwegian population register, combined with data from the 1970 and 1980 censuses. For two generations, those born in 1948 and 1958, successive migrations are considered and the effects of various sociodemographic characteristics on migratory behaviour are explored....Event-history analysis is used to examine the duration of stay at the place reached at age 22, and the effects of individual characteristics on this duration. Since a distinction is made between intra- and inter-regional migration, it has been possible to show that these effects vary significantly depending on whether we are dealing with short- or long-distance migration."     

Oil and gas resource assessment: The linkage between geological analysis and discovery process model results

Two widely-used techniques to estimate the volume of remaining oil and gas resources are discovery process modeling and geologic assessment. Both were used in a recent national assessment of oil and gas resources of the United States. Parallel estimates were obtained for 27 provinces. Geological-based estimates can typically see into areas not available to discovery process models (that is areas with little or no exploration history) and thus, on average, yield higher estimates. However, a linear relation does exist between the mean estimates obtained from these two methods. In addition, other variables were found in a multiple regression model that explained much of the difference. Thus, it is possible to perform discovery process modeling and adjust the estimates to yield results that might be expected from geological-based assessments.     

Rupture process of large earthquakes in the northern Mexico subduction zone

The Cocos plate subducts beneath North America at the Mexico trench. The northernmost segment of this trench, between the Orozco and Rivera fracture zones, has ruptured in a sequence of five large earthquakes from 1973 to 1985; the Jan. 30, 1973 Colima event (M _s 7.5) at the northern end of the segment near Rivera fracture zone; the Mar. 14, 1979 Petatlan event (M _s 7.6) at the southern end of the segment on the Orozco fracture zone; the Oct. 25, 1981 Playa Azul event (M _s 7.3) in the middle of the Michoacan gap; the Sept. 19, 1985 Michoacan mainshock (M _s 8.1); and the Sept. 21, 1985 Michoacan aftershock (M _s 7.6) that reruptured part of the Petatlan zone. Body wave inversion for the rupture process of these earthquakes finds the best: earthquake depth; focal mechanism; overall source time function; and seismic moment, for each earthquake. In addition, we have determined spatial concentrations of seismic moment release for the Colima earthquake, and the Michoacan mainshock and aftershock. These spatial concentrations of slip are interpreted as asperities; and the resultant asperity distribution for Mexico is compared to other subduction zones. The body wave inversion technique also determines theMoment Tensor Rate Functions; but there is no evidence for statistically significant changes in the moment tensor during rupture for any of the five earthquakes. An appendix describes theMoment Tensor Rate Functions methodology in detail.The systematic bias between global and regional determinations of epicentral locations in Mexico must be resolved to enable plotting of asperities with aftershocks and geographic features. We have spatially shifted all of our results to regional determinations of epicenters. The best point source depths for the five earthquakes are all above 30 km, consistent with the idea that the down-dip edge of the seismogenic plate interface in Mexico is shallow compared to other subduction zones. Consideration of uncertainties in the focal mechanisms allows us to state that all five earthquakes occurred on fault planes with the same strike (N65°W to N70°W) and dip (15±3°), except for the smaller Playa Azul event at the down-dip edge which has a steeper dip angle of 20 to 25°. However, the Petatlan earthquake does prefer a fault plane that is rotated to a more east-west orientation—one explanation may be that this earthquake is located near the crest of the subducting Orozco fracture zone. The slip vectors of all five earthquakes are similar and generally consistent with the NUVEL-predicted Cocos-North America convergence direction of N33°E for this segment. The most important deviation is the more northerly slip direction for the Petatlan earthquake. Also, the slip vectors from the Harvard CMT solutions for large and small events in this segment prefer an overall convergence direction of about N20°E to N25°E.All five earthquakes share a common feature in the rupture process: each earthquake has a small initial precursory arrival followed by a large pulse of moment release with a distinct onset. The delay time varies from 4 s for the Playa Azul event to 8 s for the Colima event. While there is some evidence of spatial concentration of moment release for each event, our overall asperity distribution for the northern Mexico segment consists of one clear asperity, in the epicentral region of the 1973 Colima earthquake, and then a scattering of diffuse and overlapping regions of high moment release for the remainder of the segment. This character is directly displayed in the overlapping of rupture zones between the 1979 Petatlan event and the 1985 Michoacan aftershock. This character of the asperity distribution is in contrast to the widely spaced distinct asperities in the northern Japan-Kuriles Islands subduction zone, but is somewhat similar to the asperity distributions found in the central Peru and Santa Cruz Islands subduction zones. Subduction of the Orozco fracture zone may strongly affect the seismogenic character as the overlapping rupture zones are located on the crest of the subducted fracture zone. There is also a distinct change in the physiography of the upper plate that coincides with the subducting fracture zone, and the Guerrero seismic gap to the south of the Petatlan earthquake is in the wake of the Orozco fracture zone. At the northern end, the Rivera fracture zone in the subducting plate and the Colima graben in the upper plate coincide with the northernmost extent of the Colima rupture zone.     

Radiocarbon redating of contaminated samples from a tropical volcano: the Mansion ‘Series’ of St Kitts,West Indies

Heavy rainfall and dense vegetation on tropical volcanoes produce abundant carbonized wood in pyroclastic deposits, in addition to easy contamination of this wood by root systems and soluble humic material. Because the physical nature of the charcoal varies, some samples are more prone to contamination. Two independent studies of the same volcano, Mt Liamuiga on St Kitts in the Lesser Antilles, sometimes using samples from the same carbonized tree, yielded a systematic difference in radiocarbon ages. An exchange of samples and a re-investigation of three physically distinct types of charcoal yielded the following results. Rare, hard, dense charcoal, lacking contamination, which had yielded a spurious age of 2860 years bp, was redated at 1845±58 years bp. Common soft, friable charcoal with good cellular structure proved to be susceptible to contamination. A field decontamination technique utilized by one group seems significant as it yields older ages than when only routine laboratory pre-treatment was used, indicating that the latter technique only partly removes the dried and hard residue produced by the decomposition of modern plant rootlets. A previous date of 24870 years bp obtained from powdery charcoal in a horizon beneath the Mansion Series contradicted ages older than 41000 years bp from common friable charcoal in the lower Mansion Series. The soft powdery charcoal was re-investigated using a sample collected a few centimeters from the original, although field decontamination of this sample was not possible, more extensive laboratory treatment yielded an age of ca. 43000 years bp, again proving that routine laboratory pretreatments are inadequate. A revised geochronology for the Mansion Series is described and a cautionary discussion is presented for the benefit of investigators using radiocarbon ages to date volcanic deposits.     

Low temperature,low pressure deformation and metamorphism of the Vangorda massive sulphide orebody,Yukon, Canada

The Vangorda orebody is a small stratiform massive sulphide orebody located in Anvil District, Yukon, Canada. The orebody consists of fineto medium-grained semi-massive and massive sulphides with a common sulphide mineralogy of pyrite, pyrrhotite, sphalerite, galena, and minor chalcopyrite. The host rocks and the sulphide lithofacies have been complexly deformed during two phases of deformation (D₁ and D₂) and associated metamorphism (M₁ and M₂). The effects of d₁ and M₁ are penetratively overprinted by D₂ and M₂. D₂ and M₂ resulted in tight to isoclinal F₂ folding of the orebody, remobilisation of the sulphides, recrystallisation and development of shear zones along the limbs of the F₂ folds. Chlorite thermometry and sulphide thermobarometry have been carried out on the host phyllites and on the sulphides. Chlorite was analysed from the S₁ and S₂ foliations in the phyllites to determine M₁ and M₂ temperatures, respectively. However, no difference was found between chlorite compositions in these foliations and a mean temperature of 363 °C was calculated from the tetrahedral A1_IV occupancy. Arsenopyrite thermometry yielded a comparable mean temperature of 336 °C. Sphalerite inclusions in M₂ pyrite porphyroblasts from D₂ shear zones were analysed for pressure using the sphalerite + hexagonal pyrrhotite + pyrite barometer. Inclusions were analysed in an attempt to determine if relic m₁ sphalerite, and hence pressure signature, was preserved. Inclusion compositions appear to reflect only M₂ conditions and yielded a mean pressure of 4.0 kb. Sphalerite + hexagonal pyrrhotite assemblages were analysed from D₂ shear zones to determine the M₂ pressure using the sphalerite + hexagonal pyrrhotite barometer. These calculations yielded a mean pressure of 6.1 kb. The M₂ temperatures and pressures calculated using these calibrations are in good agreement with those estimated from petrogenetic relationships.     

Bayes estimation: A novel approach to derivation of internally consistent thermodynamic data for minerals,their uncertainties,and correlations. Part II: Application

Internally consistent thermodynamic data, including their uncertainties and correlations, are reported for 22 phases of the quaternary system CaO-Al₂O₃-SiO₂-H₂O. These data have been derived by simultaneous evaluation of the appropriate phase properties (PP) and reaction properties (RP) by the novel technique of Bayes estimation (BE). The thermodynamic model used and the theory of BE was expounded in Part I of this paper. Part II is the follow-up study illustrating an application of BE. The input for BE comprised, among others, the a priori values for standard enthalpy of formation of the i-th phase, Δ_f H _i ⁰ , and its standard entropy, S _i ⁰ , in addition to the reaction reversal constraints for 33 equilibria involving the relevant phases. A total of 269 RP restrictions have been processed, of which 107 turned out to be non-redundant. The refined values for Δ_f H _i ⁰ and S _i ⁰ obtained by BE, including their 2σ-uncertainties, appear in Table 4; the Appendix reproduces the corresponding correlation matrix. These data permit generation of computed phase diagrams with 2σ-uncertainty envelopes based on conventional error propagation; Fig. 3 depicts such a phase diagram for the system CaO-Al₂O₃-SiO₂. It shows that the refined dataset is capable of yielding phase diagrams with uncertainty envelopes narrow enough to be geologically useful. The results in Table 4 demonstrate that the uncertainties of the prior values for Δ_f H _i ^Emphasis>0 , given in Table 1, have decreased by up to an order of magnitude, while those for S _i ⁰ improved by a factor of up to two. For comparison, Table 4 also lists the refined Δ_f H _i ⁰ and S _i ⁰ data obtained by mathematical programming (MAP), minimizing a quadratic objective function used earlier by Berman (1988). Examples of calculated phase diagrams are given to demonstrate the advantages of BE for deriving internally consistent thermodynamic data. Although P-T curves generated from both MAP and BE databases will pass through the reversal restrictions, BE datasets appear to be better suited for extrapolations beyond the P-T range explored experimentally and for predicting equilibria not constrained by reversals.     

Silicate melts at magmatic temperatures: in-situ structure determination to 1651°C and effect of temperature and bulk composition on the mixing behavior of structural units

The abundance of coexisting structural units in K-, Na-, and Li-silicate melts and glasses from 25° to 1654°C has been determined with in-situ micro-Raman spectroscopy. From these data an equilibrium constant, K_x, for the disproportionation reaction among the structural units coexisting in the melts, Si₂O₅(2Q³)SiO₃(Q²)+SiO₂(Q⁴), was calculated (K_x is the equilibrium constant derived by using mol fractions rather than activities of the structural units). From ln K_x vs l/T relationships the enthalpy (H_x) for the disproportionation reaction is in the range of-30 to 30 kJ/mol with systematic compositional dependence. In the potassium and sodium systems, where the disproportionation reaction shifts to the right with increasing temperature, the H_x increases with silica content (M/Si decreases, M=Na, K). For melts and supercooled liquids of composition Li₂O·2SiO₂ (Li/Si=1), the H_x is indistinguishable from 0. By decreasing the Li/Si to 0.667 (composition LS3) and beyond (e.g., LS4), the disproportionation reaction shifts to the left as the temperature is increased. For a given ratio of M/Si (M=K, Na, Li), there is a positive, near linear correlation between the H_x and the Z/r² of the metal cation. The slope of the H_x vs Z/r² regression lines increases as the system becomes more silica rich (i.e., M/Si is decreased). Activity coefficients for the individual structural units, _i, were calculated from the structural data combined with liquidus phase relations. These coefficients are linear functions of their mol fraction of the form _i=a lnX _i+b, where a is between 0.6 and 0.87, and X _i is the mol fraction of the unit. The value of the intercept, b, is near 0. The relationship between activity coefficients and abundance of individual structural units is not affected by temperature or the electronic properties of the alkali metal. The activity of the structural units, however, depend on their concentration, type of metal cation, and on temperature.     

Experimental study of a jotunite (hypersthene monzodiorite): constraints on the parent magma composition and crystallization conditions (P,T, f O 2) of the Bjerkreim-Sokndal layered intrusion (Norway)

The Bjerkreim-Sokndal layered intrusion is part of the Rogaland anorthosite Province of southern Norway and is made of cumulates of the anorthositemangerite-charnockite suite. This study presents experimental phase equilibrium data for one of the finegrained jotunite (Tjörn locality) occurring along its northwestern lobe. These experimental data show that a jotunitic liquid similar in composition to the Tjörn jotunite, but slightly more magnesian and with a higher plagioclase component is the likely parent of macrocyclic units (MCU)III and IV of the intrusion. The limit of the olivine stability field in the experimentally determined phase diagram as well as comparison of the Al₂O₃ content of low-Ca pyroxenes from experiments and cumulates (1.5%) yields a pressure of emplacement 5 kbar. Experimentally determined Fe-Ti oxide equilibria compared to the order of cumulus arrival in the intrusion show that the oxygen fugacity was close to FMQ (fayalite-magnetite-quartz) during the early crystallization. It subsequently decreased relative to this buffer when magnetite disappeared from the cumulus assemblage and then increased until the reentry of this mineral. Calculated densities of experimental liquids show a density increase with fractionation at 7, 10 and 13 kbar due to the predominance of plagioclase in the crystallizing assemblage. At 5 kbar and 1 atm (FMQ-1), where plagioclase is the liquidus phase, density first increases and then drops when olivine (5 kbar) or olivine+ilmenite (1 atm: FMQ-1) precipitate. At 1 atm and NNO (nickel-nickel oxide), the presence of both magnetite and ilmenite as near liquidus phases induces a density decrease. In the Bjerkreim magma chamber, oxides are early cumulus phases and liquid density is then supposed to have decreased during fractionation. This density path implies that new influxes of magma emplaced in the chamber were both hotter and denser than the resident magma. The density contrast inferred between plagioclase and the parent magma shows that this mineral was not able to sink in the magma, suggesting anin situ crystallization process.     

Correlation by Rb-Sr geochronology of garnet growth histories from different structural levels within the Tauern Window,Eastern Alps

In order to evaluate rates of tectonometamorphic processes, growth rates of garnets from metamorphic rocks of the Tauern Window, Eastern Alps were measured using Rb-Sr isotopes. The garnet growth rates were determined from Rb-Sr isotopic zonation of single garnet crystals and the Rb-Sr isotopic compositions of their associated rock matrices. Garnets were analyzed from the Upper Schieferhülle (USH) and Lower Schieferhülle, (LSH) within the Tauern Window. Two garnets from the USH grew at rates of 0.67 _–0.13 ^+0.19 mm/million years and 0.88 _–0.19 ^+0.34 mm/million years, respectively, indicating an average growth duration of 5.4±1.7 million years. The duration of growth coupled with the amount of rotation recorded by inclusion trails in the USH garnets yields an average shear-strain rate during garnet growth of 2.7 _–0.7 ^+1.2 ×10^-14 s^-1. Garnet growth in the sample from the USH occurred between 35.4±0.6 and 30±0.8 Ma. The garnet from the LSH grew at a rate of 0.23±0.015 mm/million years between 62±1.5 Ma and 30.2±1.5 Ma. Contemporaneous cessation of garnet growth in both units at 30 Ma is in accord with previous dating of the thermal peak of metamorphism in the Tauern Window. Correlation with previously published pressure-temperature paths for garnets from the USH and LSH yields approximate rates of burial, exhumation and heating during garnet growth. Assuming that theseP — T paths are applicable to the garnets in this study, the contemporaneous exhumation rates recorded by garnet in the USH and LSH were approximately 4 _–2 ⁺³ mm/year and 2±1 mm/year, respectively.