Distribution and fluxes of metals in the St. Lawrence River from the outflow of Lake Ontario to Québec City

Six anchor stations in the St. Lawrence River from the outflow of Lake Ontario to Quebec City, were occupied for ca. 24 hours in June 1987 during low flow conditions. Samples of water and suspended particulate matter were separated by continuous-flow centrifugation, and were collected every two hours.During this sampling period, fluxes of dissolved forms of zinc, copper and nickel increased in Lac St. François and downstream relative to the sum of the fluxes for the two upstream stations at the outflow of Lake Ontario. Increases in the flux of dissolved zinc and copper were pronounced below Montreal and above Lac St. Pierre. For particulate forms of metals, all five metals show that there are significant inputs in the section of the St. Lawrence River between Lac St. François and the station just above the entrance to Lac St. Pierre.The average concentrations of dissolved cadmium, lead, zinc, copper and nickel ranged from 7–23 ng/l; 9–35 ng/l; 0.434–0.939 g/l; 0.15–0.89 g/l and 0.58–1.12 g/l respectively.Regression analysis of the dissolved and particulate metal concentrations suggests that the concentration of dissolved cadmium, lead, zinc and nickel can be predicted from the regression equation and the determination of particulate metal concentration. This prediction appears to be independent of the suspended particulate matter concentration which varies from ca. 1 to 10 mg/l from the outflow of Lake Ontario to Québec City.     

Techniques to measure grain-size distributions of loamy sediments: a comparative study of ten instruments for wet analysis

Ten instruments for measuring the grain‐size distribution of loamy sediments in water were tested in a comparative study. Instruments tested were: the Malvern Mastersizer S, the Coulter LS 200, the Fritsch Analysette 22 (version C), the Horiba Partica LA‐950, the Sedigraph 5100, the Atterberg cylinder, the Coulter Multisizer 3, the CIS‐100, the EyeTech and the image analysis software Histolab. The first four instruments use laser diffraction. The Sedigraph and Atterberg cylinder are based on sedimentation. The Coulter Multisizer 3 uses impedance measurements and the last three instruments use optical techniques. Grain‐size analyses were carried out on four sediments with median grain diameters of approximately 35, 30, 12 and 9 μm, respectively, and with no particles >90 μm. The laser diffraction instruments produced the best results for the various criteria considered in this study, followed by the Coulter Multisizer 3 (electrical sensing zone technique), the sedimentation techniques and the optical methods. However, the study also showed that there is no such thing as an ‘ideal’ method because techniques may score differently according to the criterion tested. An appropriate evaluation thus requires that various criteria be considered. Issues to be taken into account are: type of sediment, quantity of sediment available, speed of measurement, complexity of the measurement protocol, processing of the data, reproducibility of the results and the specific aims of the study undertaken.     

Development of the mixed conifer forest in northern New Mexico and its relationship to Holocene environmental change

Chihuahueños Bog (2925 m) in the Jemez Mountains of northern New Mexico contains one of the few records of late-glacial and postglacial development of the mixed conifer forest in southwestern North America. The Chihuahueños Bog record extends to over 15,000 cal yr BP. An Artemisia steppe, then an open Picea woodland grew around a small pond until ca. 11,700 cal yr BP when Pinus ponderosa became established. C/N ratios, δ¹³C and δ¹⁵N values indicate both terrestrial and aquatic organic matter was incorporated into the sediment. Higher percentages of aquatic algae and elevated C/N ratios indicate higher lake levels at the opening of the Holocene, but a wetland developed subsequently as climate warmed. From ca. 8500 to 6400 cal yr BP the pond desiccated in what must have been the driest period of the Holocene there. C/N ratios declined to their lowest Holocene levels, indicating intense decomposition in the sediment. Wetter conditions returned after 6400 cal yr BP, with conversion of the site to a sedge bog as groundwater levels rose. Higher charcoal influx rates after 6400 cal yr BP probably result from greater biomass production rates. Only minor shifts in the overstory species occurred during the Holocene, suggesting that mixed conifer forest dominated throughout the record.     

The role of heterogeneous strain in the development and preservation of a polymetamorphic record in high-P granulites, western Canadian Shield

Mafic rocks in the Chipman domain of the Athabasca granulite terrane, western Canadian Shield, provide the first well‐documented record of two distinct high‐P granulite facies events in the same domain in this region. Textural relations and the results of petrological modelling (NCFMASHT system) of mafic granulites are interpreted in terms of a three‐stage tectonometamorphic history. Stage 1 involved development of the assemblage Grt + Cpx + Qtz ± Pl (M₁) from a primary Opx‐bearing igneous precursor at conditions of 1.3 GPa, 850–900 °C. Field and microstructural observations suggest that M₁ developed synchronously with an early S₁ gneissic fabric. Stage 2 is characterized by heterogeneous deformation (D₂) and synkinematic partial retrogression of the peak assemblage to an amphibole‐bearing assemblage (M₂). Stage 3 involved a third phase of deformation and a return to granulite facies conditions marked by the prograde breakdown of amphibole (Amph₂) to produce matrix garnet (Grt_3a) and the coronitic assemblage Cpx_3b + Opx_3b + Ilm_3b + Pl_3b (M_3b) at 1.0 GPa, 800–900 °C. M₁ and M_3b are correlated with 2.55 and 1.9 Ga metamorphic generations of zircon, respectively, which were dated in a separate study. Heterogeneous strain played a crucial role in both the development and preservation of these rare examples of multiple granulite facies events within single samples. Without this fortuitous set of circumstances, the apparent reaction history could have incorrectly led to an interpretation involving a single‐cycle high‐grade event. The detailed P–T–t–D history constructed for these rocks provides the best evidence to date that much of the east Lake Athabasca region experienced long‐term lower crustal residence from 2.55 to 1.9 Ga, and thus the region represents a rare window into the reactivation and ultimate stabilization processes of cratonic lithosphere.     

Three-dimensional discontinuous deformation analysis (3-D DDA) using a new contact resolution algorithm

This paper presents a new contact calculating algorithm for contacts between two polyhedra with planar boundaries in the three-dimensional discontinuous deformation analysis (3-D DDA). In this algorithm, all six type contacts in 3-D (vertex-to-face, vertex-to-edge, vertex-to-vertex, face-to-face, edge-to-edge, and edge-to-face) are simply transformed into the form of point-to-face contacts. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. In this paper, formulations of normal contact, shear contact and frictional force submatrices based on the new method are derived and the algorithm has been programmed in VC⁺⁺. Examples are provided to demonstrate the new contact rule between two blocks.     

Landslide susceptibility mapping for a landslide-prone area (Findikli,NE of Turkey) by likelihood-frequency ratio and weighted linear combination models

Landslides are very common natural problems in the Black Sea Region of Turkey due to the steep topography, improper use of land cover and adverse climatic conditions for landslides. In the western part of region, many studies have been carried out especially in the last decade for landslide susceptibility mapping using different evaluation methods such as deterministic approach, landslide distribution, qualitative, statistical and distribution-free analyses. The purpose of this study is to produce landslide susceptibility maps of a landslide-prone area (Findikli district, Rize) located at the eastern part of the Black Sea Region of Turkey by likelihood frequency ratio (LRM) model and weighted linear combination (WLC) model and to compare the results obtained. For this purpose, landslide inventory map of the area were prepared for the years of 1983 and 1995 by detailed field surveys and aerial-photography studies. Slope angle, slope aspect, lithology, distance from drainage lines, distance from roads and the land-cover of the study area are considered as the landslide-conditioning parameters. The differences between the susceptibility maps derived by the LRM and the WLC models are relatively minor when broad-based classifications are taken into account. However, the WLC map showed more details but the other map produced by LRM model produced weak results. The reason for this result is considered to be the fact that the majority of pixels in the LRM map have high values than the WLC-derived susceptibility map. In order to validate the two susceptibility maps, both of them were compared with the landslide inventory map. Although the landslides do not exist in the very high susceptibility class of the both maps, 79% of the landslides fall into the high and very high susceptibility zones of the WLC map while this is 49% for the LRM map. This shows that the WLC model exhibited higher performance than the LRM model.     

Quantitative textural analysis of packings of elongate crystals

The spatial distribution of grains in a solidifying igneous rock controls the physical properties of the crystal mush, and in turn is controlled by the rate of crystal growth and accumulation. A predominant non-spherical habit for igneous minerals brings into question the use of spherical particles in reference packings used for quantification of spatial distribution. Furthermore, variations of crystal clustering/ordering with length scale require spatial statistics which take into account the distribution of particles beyond nearest neighbours. Using random close packings of spherocylinders, we demonstrate the importance of aspect ratio for the aggregation index (usually known as R) and show that packings of spherical particles have more structure than packings of rods. The spatial distribution functions demonstrate that the plagioclase grains in the colonnade from the Holyoke basalt are clustered on a length scale of 0.5 mm. Understanding the controls on grain spatial distribution in igneous rocks will depend on the application of these techniques to well-understood environments.     

Si diffusion in zircon

Self-diffusion of Si under anhydrous conditions at 1 atm has been measured in natural zircon. The source of diffusant for experiments was a mixture of ZrO₂ and ³⁰Si-enriched SiO₂ in 1:1 molar proportions; experiments were run in crimped Pt capsules in 1-atm furnaces. ³⁰Si profiles were measured with both Rutherford backscattering spectrometry (RBS) and nuclear reaction analysis with the resonant nuclear reaction ³⁰Si(p,γ)³¹P. For Si diffusion normal to c over the temperature range 1,350–1,550°C, we obtain an Arrhenius relation D = 5.8 exp(−702 ± 54 kJ mol⁻¹/RT) m² s⁻¹ for the NRA measurements, which agrees within uncertainty with an Arrhenius relation determined from the RBS measurements [62 exp(−738 ± 61 kJ mol⁻¹/RT) m² s⁻¹]. Diffusion of Si parallel to c appears slightly faster, but agrees within experimental uncertainty at most temperatures with diffusivities for Si normal to c. Diffusion of Si in zircon is similar to that of Ti, but about an order of magnitude faster than diffusion of Hf and two orders of magnitude faster than diffusion of U and Th. Si diffusion is, however, many orders of magnitude slower than oxygen diffusion under both dry and hydrothermal conditions, with the difference increasing with decreasing temperature because of the larger activation energy for Si diffusion. If we consider Hf as a proxy for Zr, given its similar charge and size, we can rank the diffusivities of the major constituents in zircon as follows: D _Zr < D _Si << D _{O, dry} < D _{O, ‘wet’}.     

Structural investigation of the synthetic CaAn(PO4)2 (An = Th and Np) cheralite-like phosphates

The crystallographic structures of the synthetic cheralite, CaTh(PO₄)₂, and its homolog CaNp(PO₄)₂ have been investigated by X-ray diffraction at room temperature. Rietveld analyses showed that both compounds crystallize in the monoclinic system and are isostructural to monazite LnPO₄ (Ln = La to Gd). The space group is P2₁/n (I.T. = 14) with Z = 2. The refined lattice parameters of CaTh(PO₄)₂ are a = 6.7085(8) Å, b = 6.9160(6) Å, c = 6.4152(6) Å, and β = 103.71(1)° with best fit parameters R _wp = 4.87%, R _p = 3.69% and R _B = 3.99%. For CaNp(PO₄)₂, we obtained a = 6.6509(5) Å, b = 6.8390(3) Å, c = 6.3537(8) Å, and β = 104.12(6)° and R _wp = 6.74%, R _p = 5.23%, and R _B = 6.05%. The results indicate significant distortions of bond length and angles of the PO₄ tetrahedra in CaTh(PO₄)₂ and to a lesser extent in CaNp(PO₄)₂. The structural distortions were confirmed by Raman spectroscopy of CaTh(PO₄)₂. A comparison with the isostructural compounds LnPO₄ (Ln = Ce and Sm) confirmed that the substitution of the large rare earth trivalent cations with Ca²⁺ and Th⁴⁺ introduces a distortion of the PO₄ tetrahedra.