The Taconite Inlet Lakes Project: a systems approach to paleoclimatic reconstruction

A comprehensive study of meteorological, hydrological, limnological and sedimentological conditions in the watersheds of density-stratified (meromictic) lakes around Taconite Inlet, Northern Ellesmere Island, N.W.T., Canada was carried out from 1990–1992. Lakes C1 and C2 contain seawater trapped by isostatic uplift as the former embayments became isolated from the sea. These lakes, and Lake C3, contain varved sediments which provide an annually resolvable paleoclimatic record. By studing the major systems influencing sedimentation in one of these lakes (Lake C2) a better understanding of the climatic controls on varve formation, and hence on the paleoclimatic signal in the varved sediment record, was obtained. The varves of Lake C2 provide a proxy record of summer temperature for the region.This is the first in a series of papers published in this issue on the Taconite Inlet Lakes Project. These papers were collected by Dr R. S. Bradley.     

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."     

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.     

Pb−Sr−Nd−O isotopic constraints on the origin of rhyolites from the Taupo Volcanic Zone of New Zealand: evidence for assimilation followed by fractionation from basalt

A comprehensive Sr–Nd–Pb–O isotopic study is reported for rhyolites from the Maroa Volcanic Centre in the Taupo Volcanic Zone (TVZ) of New Zealand. The Sr–Nd isotopic compositions of the rhyolites (⁸⁷Sr/⁸⁶Sr=0.705236 to 0.705660 and _Nd = 2.0 to 0.2) are intermediate between those of primitive basalts (⁸⁷Sr/⁸⁶Sr=0.70387 and _Nd = 5.3) and the Torlesse basement (⁸⁷Sr/⁸⁶Sr=0.709 and _Nd = -4.5). The relatively low mantle-like oxygen isotopic compositions of ¹⁸ O = 7 ± 0.5 are consistent with the Nd-Sr isotopic constraints in that they can be accounted for by 15% to 25% crustal contamination of a basaltic parent by relatively ¹⁸ O-rich Torlesse metasediment. High precision Pb isotopic analyses of plagioclase separates from the Maroa rhyolites show that they have essentially the same compositions as the Torlesse metasedimentary terrane which is itself distinctive from the Western or Waipapa metasediments. Due to the high concentration of Pb in the Torlesse metasediments (>20 ppm) compared to the basalts (<2 ppm), the Pb isotopic composition of the volcanics may be controlled by relatively small amounts (>10%) of crustal contamination. All these results are shown to be consistent with derivation of the rhyolites by 15% to 25% contamination of relatively primitive basaltic magmas with Torlesse metasedimentary crust, followed by extensive, essentially closed system fractionation of the basalt to a magma of rhyolite composition. It is argued that the processes of assimilation and fractionation are separated in both space and time. The voluminous high silica rhyolites, which make up >97% of the exposed volcanism in the continental margin back-are basin environment of the TVZ, therefore appear to be a product of predominantly new additions to the crust with assimilation-recycling of pre-existing crust being of secondary importance.     

Aqueous alteration of the Bali CV3 chondrite: evidence from mineralogy, mineral chemistry, and oxygen isotopic compositions

A petrographic, geochemical, and oxygen isotopic study of the Bali CV3 carbonaceous chondrite revealed that the meteorite has undergone extensive deformation and aqueous alteration on its parent body. Deformation textures are common and include flattened chondrules, a well-developed foliation, and the presence of distinctive (100) planar defects in olivine. The occurrence of alteration products associated with the planar defects indicates that the deformation features formed prior to the episode of aqueous alteration. The secondary minerals produced during the alteration event include well-crystallized Mg-rich saponite, framboidal magnetite, and Ca-phosphates. The alteration products are not homogeneously distributed throughout the meteorite, but occur in regions adjacent to relatively unaltered material, such as veins of altered material following the foliation. The alteration assemblage formed under oxidizing conditions at relatively low temperatures (<100 degrees C). Altered regions in Bali have higher Na, Ca, and P contents than unaltered regions which suggests that the fluid phase carried significant dissolved solids. Oxygen isotopic compositions for unaltered regions in Bali fall within the field for other CV3 whole-rocks, however, the oxygen isotopic compositions of the heavily altered material lie in the region for the CM and CR chondrites. The heavy-isotope enrichment of the altered regions in Bali suggest alteration conditions similar to those for the petrographic type-2 carbonaceous chondrites.     

Crustal Melting in the Cordilleran Interior: The Mid-Cretaceous White Creek Batholith in the Southern Canadian Cordillera

The mid-Cretaceous White Creek batholith in southeast BritishColumbia is a zoned pluton ranging from quartz monzodioriteon the margin, to hornblende-and biotite-bearing granodioritetowards the interior of the batholith, which are in turn crosscutby two-mica granite. This range in rock type is similar to therange displayed by Mesozoic granitoid suites found in the Cordilleraninterior of western North America. The lithological zones inthe White Creek batholith correlate with distinet jumps in majorelement, trace element, and isotopic compositions, and indicatethat several pulses of magma were emplaced within the WhiteCreek magma chamber. The hornblende-and biotite-bearing granitoidsare metaluminous to weakly peraluminous, have strong light rareearth element (LREE) enrichment, and small negative Eu anomalies.These granitoids have initial _Sr ranging from +32 to +84 (⁸⁷Sr/⁸⁶SrTfrom 0.7069 to 0.7106), initial _Nd ranging from –5 to–10, and initial ²⁰⁶Pb/²⁰⁴Pb, ²⁰⁷Pb/²⁰⁴Pb, and ²⁰⁸Pb/²⁰⁴Pbranging from 18.3 to 18.7, 15.58 to 15.65, and 38.3 to 39.0,respectively. The two-mica granites and associated aplites arestrongly peraluminous, and show only moderate LREE enrichmentand strong negative Eu anomalies. These granites have _Sr rangingfrom +174 to + 436 (⁸⁷Sr/⁸⁶Sr_T from 0.7169 to 0.7354), _Nd rangingfrom –12 to –16, and more radiogenic initial Pbisotope ratios than the hornblende-and biotite-bearing granitoids. Oxygen, Sr, Pb, and Nd isotopes, REE modelling, and phase equilibriumconstraints are consistent with crustal anatexis of Precambrianbasement gneisses and Proterozoic metapelites exposed in southeastBritish Columbia, the product being the hornblende-biotite granitoidsand two-mica granites, respectively. The sequence of intrusionin the White Creek batholith constrains the melting sequence.A zone of anatexis proceeded upwards through the crust, firstmelting basement gneisses then melting overlying metapelites.A model for basaltic magmatic underplating as a primary causeof anatexis of the crust during the mid-Cretaceous magmaticepisode is difficult to reconcile with the absence of earlyCretaceous basalt in the southern Canadian Cordillera. A muchmore likely petrogenetic model is that crustal anatexis wasprobably a response to crustal thickening in association withterrane accretion and collision along the western margin ofthe North American continent.