A coupled CFD-DEM investigation of suffusion of gap graded soil: Coupling effect of confining pressure and fines content

Suffusion involves fine particles migration within the matrix of coarse fraction under seepage flow, which usually occurs in the gap-graded material of dams and levees. Key factors controlling the soil erodibility include confining pressure (p′) and fines content (F_c), of which the coupling effect on suffusion still remains contradictory, as concluded from different studies considering narrow scope of these factors. For this reason, a systematical numerical simulation that considers a relative wide range of p′ and F_c was performed with the coupled discrete element method and computational fluid dynamics approach. Two distinct macroresponses of soil suffusion to p′ were revealed, ie, for a given hydraulic gradient i = 2, an increase in p′ intensifies the suffusion of soil with fines overfilling the voids (eg, F_c = 35%), but have negligible effects on the suffusion of gap-graded soil containing fines underfilling the voids (eg, F_c = 20%). The micromechanical analyses, including force chain buckling and strain energy release, reveal that when the fines overfilled the voids between coarse particles (eg, F_c = 35%) and participated heavily in load-bearing, the erosion of fines under high i could cause the collapse of the original force transmission structure. The release of higher strain energy within samples under higher p′ accelerated particle movement and intensified suffusion. Conversely, in the case where the fines underfilled the voids between coarse particles (eg, F_c = 20%), the selective erosion of fines had little influence on the force network. High p′ in this case prevented suffusion.     

Suspended Matter,Chl-a,CDOM, Grain Sizes,and Optical Properties in the Arctic Fjord-Type Estuary,Kangerlussuaq, West Greenland During Summer

Optical constituents as suspended particulate matter (SPM), chlorophyll (Chl-a), colored dissolved organic matter (CDOM), and grain sizes were obtained on a transect in the arctic fjord-type estuary Kangerlussuaq (66°) in August 2007 along with optical properties. These comprised diffuse attenuation coefficient of downwelling PAR (K _d(PAR)), upwelling PAR (K _u(PAR)), particle beam attenuation coefficient (c _p), and irradiance reflectance R(−0, PAR). PAR is white light between 400 and 700 nm. The estuary receives melt water from the Greenland Inland Ice and stations covered a transect from the very high turbid melt water outlet to clear marine waters. Results showed a strong spatial variation with high values as for suspended matter concentrations, CDOM, diffuse attenuation coefficient K _d(PAR), particle beam attenuation coefficients (c _p), and reflectance R(−0, PAR) at the melt water outlet. Values of optical constituents and properties decreased with distance from the melt water outlet to a more or less constant level in central and outer part of the estuary. There was a strong correlation between inorganic suspended matter (SPMI) and diffuse attenuation coefficient K _d(PAR) (r ² = 0.92) and also for particle beam attenuation coefficient (c _p; r ² = 0.93). The obtained SPMI specific attenuation—K _d^*(PAR) = 0.13 m² g⁻¹ SPMI—and the SPMI specific particle beam attenuation—c _p^* = 0.72 m² g⁻¹—coefficients were about two times higher than average literature values. Irradiance reflectance R(−0, PAR) was comparatively high (0.09−0.20) and showed a high (r ² = 0.80) correlation with K _u(PAR). Scattering dominated relative to absorption—b(PAR)/a(PAR) = 12.3. Results strongly indicated that the high values in the optical properties were related to the very fine particle sizes (mean = 2–6 μm) of the suspended sediment. Data and results are discussed and compared to similar studies from both temperate and tropical estuaries.     

Conditions of diamond growth: a proton microprobe study of inclusions in West Australian diamonds

Crystalline primary inclusions in diamonds from the Argyle and Ellendale lamproites have been analyzed for Mn, Ni, Cu, Zn, Ga, Pb, Rb, Sr, Y, Zr, Nb, Ta, Ba and Mo by proton microprobe. Eclogite-suite inclusions dominate at Argyle and occur in equal proportions with peridotite-suite inclusions at Ellendale. Eclogitic phases present include garnet, omphacitic clinopyroxene, coesite, rutile, kyanite and sulfide. Eclogitic clinopyroxenes are commonly rich in K and contain 300–1060 ppm Sr and 3–70 ppm Zr: K/Rb increases with K content up to 1400 at 0.7–1.1% K. Rutiles have high Zr and Nb contents with Zr/Nb=1.5–4 and Nb/Ta 16. Of the peridotite-suite inclusions, olivine commonly contains > 10 ppm Sr and Mo; Cr-pyropes are depleted in Sr, Y and Zr, and enriched in Ni, relative to eclogitic garnets.Eclogite-suite diamonds grew in host rocks that were depleted in Mn, Ni and Cr, and enriched in Sr, Zn, Cu, Ga and Ti, relative to Type I eclogite xenoliths from the Roberts Victor Mine. Crystallization temperatures of the eclogite-suite diamonds, as determined by coexisting garnet and clinopyroxene from single diamonds, range from 1085 to 1575° C. Log K _D (C _i ^cpx /C _i ^gnt ) varies linearly with 1/T for Zr, Sr and Ga in most of the same samples. This supports the validity of the temperature estimates; Argyle eclogite-suite diamonds have grown over a T range 400° C. Comparison with data from eclogite xenoliths in kimberlites suggests that K _D ^Sr and K _D ^Zr are mainly T-dependent, while K _D ^Ga may be both temperature-and pressuredependent. K _D ^Ni , K _D ^Cu and K _D ^Zn show no T dependence in these samples.In several cases, significant major-and/or trace-element disequilibrium is observed between different grains of the same mineral, or between pyroxene and garnet, within single diamonds. This implies that these diamonds grew in an open system; inclusions trapped at different stages of growth record changes in major and trace-element composition occurring in the host rock. Diamond growth may have been controlled by a fluid flux which introduced or liberated carbon and modified the composition of the rock. The wide range of equilibration temperatures and the range of composition recorded in the inclusions of single diamonds suggest that a significant time interval was involved in diamond growth.     

The high-pressure and temperature equation of state of a majorite solid solution in the system of Mg4Si4O12-Mg3Al2Si3O12

The high-pressure and temperature equation of state of majorite solid solution, Mj_0.8Py_0.2, was determined up to 23 GPa and 773 K with energy-dispersive synchrotron X-ray diffraction at high pressure and high temperature using the single- and double-stage configurations of the multianvil apparatuses, MAX80 and 90. The X-ray diffraction data of the majorite sample were analyzed using the WPPD (whole-powder-pattern decomposition) method to obtain the lattice parameters. A least-squares fitting using the third-order Birch-Murnaghan equation of state yields the isothermal bulk modulus, K _T0  = 156 GPa, its pressure derivative, K′ = 4.4(±0.3), and temperature derivative (∂K _T /∂T) _P = −1.9(±0.3)× 10⁻² GPa/K, assuming that the thermal expansion coefficient is similar to that of pyrope-almandine solid solution. Received: 5 October 1998 / Revised, accepted: 24 June 1999     

P-V equation of state of portlandite, Ca(OH)2, from powder neutron diffraction data

 A high pressure neutron powder diffraction study of portlandite [Ca(OH)₂] has been performed at ISIS facility (U.K.); nine spectra have been collected increasing the pressure by steps, up to 10.9 GPa, by means of a Paris-Edinburgh cell installed on the POLARIS diffractometer. The tensorial formalism of the lagrangian finite strain theory and the Birch-Murnaghan equation of state have been used to determine, independently, two values of the bulk modulus of portlandite, obtaining K ₀=38.3(±1.1) GPa [linear incompressibilities: K _0a=188.4(±9.9), K _0c=64.5(±2.5) GPa] and K ₀=34.2(±1.4) GPa, respectively. The present results comply with values from previous measurements by X-ray diffraction [K ₀=37.8(±1.8) GPa] and Brillouin spectroscopy [K ₀=31.7(±2.5) GPa]. Reasonably, Ca(OH)₂ has revealed to be bulkly softer than Mg(OH)₂ [K ₀=41(±2), K _0a=313, K _0c=57 GPa]. The Ca(OH)₂ linear incompressibility values reflect the nature of forces acting to stabilize the (001) layer structure and, further, prove that the replacement Ca/Mg mainly affects the elastic properties in the (001) plane, rather than along the [001] direction. Data from a full refinement of the structure at room pressure are reported. Received January 12, 1996/Revised, accepted June 15, 1996     

Insights into oxygen-cation bonding in fresnoite-type structures from O <Emphasis Type="Italic">K</Emphasis>- and Ti <Emphasis Type="Italic">L</Emphasis><Subscript>23</Subscript>-electron energy-loss spectra and ab initio calculations of the electronic structure

O K- and Ti L₂₃-core-loss spectra of fresnoite Ba₂TiSi₂O₈ (BTS) and Sr₂TiSi₂O₈ (STS), which is isotypic to BTS, have been measured by electron energy-loss spectroscopy (EELS). The energy-loss near-edge structures (ELNES) of the O K edge have been identified on the basis of theoretical simulations and interpretations of the X-ray absorption near-edge structures (XANES), which have been modelled in the framework of self-consistent full multiple-scattering (FMS) theory using FEFF8. Herewith, the K-absorption spectra of oxygen (E) and the local partial electron density of states (DOS) of all atoms have been calculated. For BTS, the observed spectral features in the O K-edge spectra are interpreted in terms of mixing between the central O p and neighbouring Ba 5d and 4f, Si 3p and 3d, and Ti 3d orbitals. The observed differences in the O K-edge spectra for STS and BTS can mainly be attributed to three properties: (1) The lack of high local partial Sr unoccupied DOS with 4f symmetry near the Fermi level compared to the high Ba 4f unoccupied DOS results in differences of overlapping O 2p – cation orbitals. (2) The differences in the ionic radii of Sr and Ba result in a larger unit cell for BTS and, thus, in larger oxygen-cation bonding distances. (3) In comparison to STS, the strength of the incommensurate 2-D structural modulation is significantly weaker in BTS, i.e. distortions of coordination polyhedra occur to a much lesser extent. All these effects alter the oxygen-cation hybridization and, hence, result in a variation of the O 1s p transition and consequently of the O K-edge spectral shape. The observed peak broadening in Ti L₂₃ ELNES of STS compared to BTS is correlated with strong displacive modulations hosted in STS.     

Pressure–volume–temperature equation of state of andradite and grossular, by high-pressure and -temperature powder diffraction

 The thermoelastic parameters of natural andradite and grossular have been investigated by high-pressure and -temperature synchrotron X-ray powder diffraction, at ESRF, on the ID30 beamline. The P–V–T data have been fitted by Birch-Murnaghan-like EOSs, using both the approximated and the general form. We have obtained for andradite K ₀=158.0(±1.5) GPa, (dK/dT )₀=−0.020(3) GPa K⁻¹ and α₀=31.6(2) 10⁻⁶ K⁻¹, and for grossular K ₀=168.2(±1.7) GPa, (dK/dT)₀=−0.016(3) GPa K⁻¹ and α₀=27.8(2) 10⁻⁶ K⁻¹. Comparisons between the present issues and thermoelastic properties of garnets earlier determined are carried out. Received: 7 July 2000 / Accepted: 20 October 2000     

On the elastic behaviour of zeolite mordenite: a synchrotron powder diffraction study

The high-pressure elastic behaviour of a synthetic zeolite mordenite, Na₆Al_6.02Si_42.02O₉₆·19H₂O [a=18.131(2), b=20.507(2), c=7.5221(5) Å, space group Cmc2₁], has been investigated by means of in situ synchrotron X-ray powder diffraction up to 5.68 GPa. No phase transition has been observed within the pressure range investigated. Axial and volume bulk moduli have been calculated using a truncated second-order Birch–Murnaghan equation-of-state (II-BM-EoS). The refined elastic parameters are: V ₀=2801(11) ų, K _T0= 41(2) GPa for the unit-cell volume; a ₀=18.138(32) Å, K _T0(a)=70(8) GPa for the a-axis; b ₀=20.517(35) Å, K _T0(b)=29(2) GPa for the b-axis and c ₀=7.531(5) Å, K _T0(c)=38(1) GPa for the c-axis [K _T0(a): K _T0(b): K _T0(c)=2.41:1.00:1.31]. Axial and volume Eulerian finite strain versus “normalized stress” plots (fe–Fe plot) show an almost linear trend and the weighted linear regression through the data points yields the following intercept values: Fe(0)=39(4) GPa for V; Fe _a (0)=65(18) GPa for a; Fe _b (0)=28(3) GPa for b; Fe _c (0)=38(2) GPa for c. The magnitudes of the principal Lagrangian unit-strain coefficients, between 0.47 GPa (the lowest HP-data point) and each measured P>0.47 GPa, were calculated. The unit-strain ellipsoid is oriented with ε₁ || b, ε₂ || c, ε₃ || a and |ε₁|> |ε₂|> |ε₃|. Between 0.47 and 5.68 GPa the relationship between the unit-strain coefficient is ε₁: ε₂: ε₃=2.16:1.81:1.00. The reasons of the elastic anisotropy are discussed.An erratum to this article can be found at     

Thermodynamic properties,low-temperature heat-capacity anomalies,and single-crystal X-ray refinement of hydronium jarosite, (H<Subscript>3</Subscript>O)Fe<Subscript>3</Subscript>(SO<Subscript>4</Subscript>)<Subscript>2</Subscript>(OH)<Subscript>6</Subscript>

Crystals of hydronium jarosite were synthesized by hydrothermal treatment of Fe(III)–SO₄ solutions. Single-crystal XRD refinement with R₁=0.0232 for the unique observed reflections (|F_o| > 4_F) and wR₂=0.0451 for all data gave a=7.3559(8) Å, c=17.019(3) Å, V^o=160.11(4) cm³, and fractional positions for all atoms except the H in the H₃O groups. The chemical composition of this sample is described by the formula (H₃O)_0.91Fe_2.91(SO₄)₂[(OH)_5.64(H₂O)_0.18]. The enthalpy of formation (H^o_f) is –3694.5 ± 4.6 kJ mol^–1, calculated from acid (5.0 N HCl) solution calorimetry data for hydronium jarosite, -FeOOH, MgO, H₂O, and -MgSO₄. The entropy at standard temperature and pressure (S^o) is 438.9±0.7 J mol^–1 K^–1, calculated from adiabatic and semi-adiabatic calorimetry data. The heat capacity (C_p) data between 273 and 400 K were fitted to a Maier-Kelley polynomial C_p(T in K)=280.6 + 0.6149T–3199700T^–2. The Gibbs free energy of formation is –3162.2 ± 4.6 kJ mol^–1. Speciation and activity calculations for Fe(III)–SO₄ solutions show that these new thermodynamic data reproduce the results of solubility experiments with hydronium jarosite. A spin-glass freezing transition was manifested as a broad anomaly in the C_p data, and as a broad maximum in the zero-field-cooled magnetic susceptibility data at 16.5 K. Another anomaly in C_p, below 0.7 K, has been tentatively attributed to spin cluster tunneling. A set of thermodynamic values for an ideal composition end member (H₃O)Fe₃(SO₄)₂(OH)₆ was estimated: G^o_f= –3226.4 ± 4.6 kJ mol^–1, H^o_f=–3770.2 ± 4.6 kJ mol^–1, S^o=448.2 ± 0.7 J mol^–1 K^–1, C_p (T in K)=287.2 + 0.6281T–3286000T^–2 (between 273 and 400 K).     

Computer modelling of a pressure induced phase change in clinoenstatite pyroxenes

 Using lattice dynamic modelling of pure MgSiO₃ clinopyroxenes, we have be able to simulate the properties of both the low-clino (P2₁/c) and a high-density-clino (C2/c) phases and our results are comparable with the high pressure (HP) X-ray study of these phases (Angel et al. 1992). The transition between the two phases is predicted to occur at 6GPa. The volume variation with pressure for both phases is described by a third-order Birch-Murnaghan equation of state with the parameters V _{0  low}=31.122 cm³·mol⁻¹, K _{T0  low}= 107.42 GPa, K′ _{T0  low}=5.96, V _{0  high}=30.142 cm³·mol^–1, K _{T0  high}102.54 GPa and K′ _{T0  high}=8.21. The change in entropy between the two modelled phases at 6GPa is ΔS _{6 Gpa}=−1.335 J·mol⁻¹·K⁻¹ and the equivalent change in volume is ΔV _{6 GPa}=−0.92 cm³·mol⁻¹, from which the gradient of the phase boundary δP/δT is 0.0014 GPa·K⁻¹. The variation of the bulk modulus with pressure was also determined from the modelled elastic constants and compares very well with the EOS data. The reported Lehmann discontinuity, ∼220 km depth and pressure of 7.11Gpa, has an increase in the seismic compressional wave velocity, v _p , of 7.14% using the data given for PREM (Anderson 1989). At a pressure of 7GPa any phase transition of MgSiO₃ pyroxene would be between ortho (Pbca) and high-clino. We find the value of v _p at 7GPa, for modelled orthoenstatite (Pbca), is 8.41 km·sec⁻¹ and that for the modelled high-clino phase at 7GPa is 8.93 km·sec⁻¹, giving a dv _p /v _p of 6.18%. Received: July 26, 1996 / Revised, accepted: September 27, 1996     

Behavior of epidote at high pressure and high temperature: a powder diffraction study up to 10 GPa and 1,200 K

The thermo-elastic behavior of a natural epidote [Ca_1.925 Fe_0.745Al_2.265Ti_0.004Si_3.037O₁₂(OH)] has been investigated up to 1,200 K (at 0.0001 GPa) and 10 GPa (at 298 K) by means of in situ synchrotron powder diffraction. No phase transition has been observed within the temperature and pressure range investigated. P–V data fitted with a third-order Birch–Murnaghan equation of state (BM-EoS) give V ₀ = 458.8(1)ų, K _T0 = 111(3) GPa, and K′ = 7.6(7). The confidence ellipse from the variance–covariance matrix of K _T0 and K′ from the least-square procedure is strongly elongated with negative slope. The evolution of the “Eulerian finite strain” vs “normalized stress” yields Fe(0) = 114(1) GPa as intercept values, and the slope of the regression line gives K′ = 7.0(4). The evolution of the lattice parameters with pressure is slightly anisotropic. The elastic parameters calculated with a linearized BM-EoS are: a ₀ = 8.8877(7) Å, K _T0(a) = 117(2) GPa, and K′(a) = 3.7(4) for the a-axis; b ₀ = 5.6271(7) Å, K _T0(b) = 126(3) GPa, and K′(b) = 12(1) for the b-axis; and c ₀ = 10.1527(7) Å, K _T0(c) = 90(1) GPa, and K’(c) = 8.1(4) for the c-axis [K _T0(a):K _T0(b):K _T0(c) = 1.30:1.40:1]. The β angle decreases with pressure, β_P(°) = β_P0 −0.0286(9)P +0.00134(9)P ² (P in GPa). The evolution of axial and volume thermal expansion coefficient, α, with T was described by the polynomial function: α(T) = α₀ + α₁ T ^−1/2. The refined parameters for epidote are: α₀ = 5.1(2) × 10⁻⁵ K⁻¹ and α₁ = −5.1(6) × 10⁻⁴ K^1/2 for the unit-cell volume, α₀(a) = 1.21(7) × 10⁻⁵ K⁻¹ and α₁(a) = −1.2(2) × 10⁻⁴ K^1/2 for the a-axis, α₀(b) = 1.88(7) × 10⁻⁵ K⁻¹ and α₁(b) = −1.7(2) × 10⁻⁴ K^1/2 for the b-axis, and α₀(c) = 2.14(9) × 10⁻⁵ K⁻¹ and α₁(c) = −2.0(2) × 10⁻⁴ K^1/2 for the c-axis. The thermo-elastic anisotropy can be described, at a first approximation, by α₀(a): α₀(b): α₀(c) = 1 : 1.55 : 1.77. The β angle increases continuously with T, with β_T(°) = β_T0 + 2.5(1) × 10⁻⁴ T + 1.3(7) × 10⁻⁸ T ². A comparison between the thermo-elastic parameters of epidote and clinozoisite is carried out.     

P-V-T equation of state of lawsonite

A pressure-volume-temperature data set has been obtained for lawsonite [CaAl₂Si₂O₇(OH)₂.H₂O], using synchrotron X-ray diffraction and an externally heated diamond anvil cell. Unit-cell volumes were measured to 9.4 GPa and 767 K by angle dispersive X-ray diffraction using imaging plates. Phase changes were not observed within this pressure-temperature range, and lawsonite compressed almost isotropically at constant temperature. The P-V-T data have been analyzed using a Birch- Murnaghan equation of state and a linear equation of state expressed as β=–1/V₀ (∂V/∂P) _T . At room temperature, the derived equation of state parameters are: K ₀=124.1 (18) GPa K'₀ set to 4) and β^–1=142.0(24) GPa, respectively. Our results are intermediate between previously reported measurements. The high-temperature data show that the incompressibility of lawsonite decreases with increasing temperature to ∼500 K and then increases above. Hence, the second order temperature derivative of the bulk modulus is taken into account in the equation of state; a fit of the volume data yields K ₀=123.9(18) GPa, (∂K/∂T)_P=–0.111(3) GPa K^–1, (∂² K/∂T ²)_P=0.28(6) 10^–3 GPa K^–2, α₀=3.1(2) 10^–5 K^–1, assuming K'₀=4. Received: 2 June 1998 / Revised, accepted: 12 Ocotber 1998     

Early Cretaceous adakitic rocks in the Anqing region,southeastern China: constraints on petrogenesis and metallogenic significance

ABSTRACT

The Anqing region in Lower Yangtze River metallogenic belt is one of the important Cu polymetal producers in China. The origin of Cu polymetallic deposits in the region is closely related to Early Cretaceous adakitic intrusions. To constrain the petrogenetic and metallogenic significance of the adakitic rocks, a detailed geochronological, geochemical, and Sr–Nd–Pb–Hf isotopic study was performed. The Anqing adakitic rocks (SiO₂ = 57.4–64.2 wt.%) consist mainly of quartz monzodiorite, formed at 138.2 ± 1.7 Ma (Mean Standard Weighted Deviation (MSWD) = 0.61). They have high MgO, Al₂O₃, Sr, and low Rb, Y, Yb contents, together with high Sr/Y (50.5–222) and La/Yb (31.9–46.9) ratios. They also show negative whole-rock ε_Nd(t) (?9.8 to ?8.5) and zircon ε_Hf(t) (?10.0 to ?5.4), and high oxygen fugacity (mainly ?17.0 to ?8.01) values and radiogenic Pb isotopic compositions with (²⁰⁶Pb/²⁰⁴Pb)_i = 17.692–17.884, (²⁰⁷Pb/²⁰⁴Pb)_i = 15.413–15.511, and (²⁰⁸Pb/²⁰⁴Pb)_i = 37.611–37.943. Coupled with negative Nb–Ta anomalies, low K₂O/Na₂O ratios (0.39–0.62), and high Mg# values (0.44–0.71), these data suggest the adakitic rocks and associated large-scale Cu–Au mineralization of the Anqing region resulted from partial melting of the high oxidized subducted oceanic crust. Addition of mantle-derived magmas and assimilation of crustal materials during emplacement are also possible.     

Dissociation Constants of Protonated Oxidized Glutathione in Seawater Media at Different Salinities

Oxidized glutathione (GSSG), which has four carboxylic and two amino groups, interacts with metal ions and may affect the bioavailability and geochemistry of metals in natural waters. In the present paper, six stepwise protonation constants K^\textH_i K^{\text{H}}_{i} for GSSG were measured as a function of salinity, S = 5–35‰ at t = 25°C (and in NaCl/MgCl₂ mixtures at different ionic strengths), in order to provide thermodynamic data for their acid base properties, which are useful for studying the interaction with metals in these media. The protonation enthalpies (ΔH _i /kJ mol⁻¹) were also determined at t = 25°C. The results were interpreted using the SIT model and Pitzer equations. The seawater model with the interaction parameters accounts for the differences between the values in NaCl and seawater. The results suggest that it is important to consider all of the ionic interactions in natural waters in examining the proton dissociation of GSSG.     

Horizontal sighting range and Secchi depth as estimators of underwater PAR attenuation in a coastal lagoon

Attenuation of photosynthetically available radiation (PAR) measured using a light meter, was related to Secchi disk, horizontal black disk and horizontal sighting ranges observed in a coastal lagoon of the Southern California Current System. Vertical attenuation coefficient (K_PAR) was calculated from radiometric PAR profiles. Vertical (Z_D) and horizontal (HS) sighting ranges were measured with white (Secchi depth or Z_SD, HS _W ) and black (Z _BD, HS _B ) targets. Empirical power models for the K_PAR-Z_SD (K_PAR=1.47 Z_SD ^−1.13), K_PAR-Z _BD (K_PAR=0.98 Z _BD ^−1.26), K_PAR-HS _W (K_PAR=1.22 HS _W ^−1.14) and K_PAR-HS _B (K_PAR=0.73 HS _B ^−1.07) relationships were developed. The parameters of these models may not apply to other water bodies because their values depend on the range of water reflectance in each case, as reported in the literature. This is the first contribution reporting K_PAR-HS empirical relations in an estuarine environment but their application may be limited to this coastal lagoon. While this approach may be universal, more data are needed to explore the variability of the parameters between different water bodies.     

Geochemical,Sr–Nd–Pb isotope,and zircon U–Pb geochronological constraints on the origin of Early Permian mafic dikes,northern North China Craton

Dolerite dike swarms are widespread across the North China Craton (NCC) of Hebei Province (China) and Inner Mongolia. Here, we report new geochemical, Sr–Nd–Pb isotope, and U–Pb zircon ages for representative samples of these dikes. Laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) U–Pb analysis yielded consistent Permian ages of 274.8 ± 2.9 and 275.0 ± 4.5 Ma for zircons extracted from two dikes. The dolerites have highly variable compositions (SiO₂ = 46.99–56.18 wt.%, TiO₂ = 1.27–2.39 wt.%, Al₂O₃ = 14.42–16.20 wt.%, MgO = 5.18–7.75 wt.%, Fe₂O₃ = 8.03–13.52 wt.%, CaO = 5.18–9.75 wt.%, Na₂O = 2.46–3.79 wt.%, K₂O = 0.26–2.35 wt.%, and P₂O₅ = 0.18–0.37 wt.%) and are light rare earth element (LREE) and large ion lithophile element (LILE, e.g. Rb, Ba, and K, and Pb in sample SXG1-9) enriched, and Th and high field strength element (HFSE, e.g. Nb and Ta in sample SXG1-9, and Ti) depleted. The mafic dikes have relatively uniform (⁸⁷Sr/⁸⁶Sr)_i values from 0.7031 to 0.7048, (²⁰⁶Pb/²⁰⁴Pb)_i from 17.77 to 17.976, (²⁰⁷Pb/²⁰⁴Pb)_i from 15.50 to 15.52, (²⁰⁸Pb/²⁰⁴Pb)_i from 37.95 to 38.03, and positive ?_Nd(t) (3.6–7.3), and variable neodymium model ages (T_DM1 = 0.75–0.99 Ga, T_DM2 = 0.34–0.74 Ga). These data suggest that the dike magmas were derived from partial melting of a depleted region of the asthenospheric mantle, and that they fractionated olivine, pyroxene, plagioclase, K-feldspar, and Ti-bearing phases without undergoing significant crustal contamination. These mafic dikes within the NCC formed during a period of crustal thinning in response to extension after Permian collision between the NCC and the Siberian Block.     

Proton and metal ion binding to natural organic polyelectrolytes—II. Preliminary investigation with a peat and a humic acid

A unified physico-chemical model, based on a modified Henderson-Hasselbalch equation, for the analysis of ion complexation reactions involving charged polymeric systems is presented and verified. In this model pH = pK_a+p(ΔK_a) + log(α/1 − α) where K_a is the intrinsic acid dissociation constant of the ionizable functional groups on the polymer, ΔK_a is the deviation of the intrinsic constant due to electrostatic interaction between the hydrogen ion and the polyanion, and alpha (α) is the polyacid degree of ionization. Using this approach pK_a values for repeating acidic units of polyacrylic (PAA) and polymethacrylic (PMA) acids were found to be 4.25 ± 0.03 and 4.8 ± 0.1, respectively. The polyion electrostatic deviation term derived from the potentiometric titration data (i.e. p(ΔK_a)) is used to calculate metal ion concentration at the complexation site on the surface of the polyanion. Intrinsic cobalt-polycarboxylate binding constants (7.5 for PAA and 5.6 for PMA), obtained using this procedure, are consistent with the range of published binding constants for cobalt-monomer carboxylate complexes. In two phase systems incorporation of a Donnan membrane potential term allows determination of the intrinsic pK_a of a cross-linked PMA gel, pK_a = 4.83, in excellent agreement with the value obtained for the linear polyelectrolyte and the monomer. Similarly, the intrinsic stability constant for cobalt ion binding to a PMA-gel (β_CoPMA+ = 11) was found to be in agreement with the linear polyelectrolyte analogue and the published data for cobalt-carboxylate monodentate complexes.     

Spatial and multivariate analysis of geochemical data from metavolcanic rocks in the Ben Nevis area, Ontario

A study of the lithogeochemistry of metavolcanics in the Ben Nevis area of Ontario, Canada has shown that factor analysis methods can distinguish lithogeochemical trends related to different geological processes, most notably, the principal compositional variation related to the volcanic stratigraphy and zones of carbonate alteration associated with the presence of sulphides and gold. Auto- and cross-correlation functions have been estimated for the two-dimensional distribution of various elements in the area. These functions allow computation of spatial factors in which patterns of multivariate relationships are dependent upon the spatial auto- and cross-correlation of the components. Because of the anisotropy of primary compositions of the volcanics, some spatial factor patterns are difficult to interpret. Isotropically distributed variables such as CO ₂ are delineated clearly in spatial factor maps. For anisotropically distributed variables (SiO ₂ ), as the neighborhood becomes smaller, the spacial factor maps becomes better. Interpretation of spatial factors requires computation of the corresponding amplitude vectors from the eigenvalue solution. This vector reflects relative amplitudes by which the variables follow the spatial factors. Instability of some eigenvalue solutions requires that caution be used in interpreting the resulting factor patterns. A measure of the predictive power of the spatial factors can be determined from autocorrelation coefficients and squared multiple correlation coefficients that indicate which variables are significant in any given factor. The spatial factor approach utilizes spatial relationships of variables in conjunction with systematic variation of variables representing geological processes. This approach can yield potential exploration targets based on the spatial continuity of alteration haloes that reflect mineralization.List of symbols c _i Scalar factor that minimizes the discrepancy between andU _i - D Radius of circular neighborhood used for estimating auto- and cross-correlation coefficients - d Distance for which transition matrixU is estimated - d _ij Distance between observed valuesi andj - E Expected value - E _i Row vector of residuals in the standardized model - F(d _ij) Quadratic function of distanced _ij F(d _ij)=a+bd _ij+cd _ij ² - L Diagonal matrix of the eigenvalues ofU - _i Eigenvalue of the matrixU;ith diagonal element ofL - N Number of observations - p Number of variables - Q Total predictive power ofU - R Correlation matrix of the variables - R _0j Variance-covariance signal matrix of the standardized variables at origin;j is the index related tod andD (e.g.,j=1 ford=500 m,D=1000 m) - R _1j Matrix of auto- and cross-correlation coefficients evaluated at a given distance within the neighborhood - R _m ² Multiple correlation coefficient squared for themth variable - S _i Column vectori of the signal values - s _k ² Residual variance for variablek - T _i Amplitude vector corresponding toV _i;ith row ofT=V ^–1 - T Total variation in the system - U Nonsymmetric transition matrix formed by post-multiplyingR ₀₁ ^–1 byR _ij - U _i Componenti of the matrixU, corresponding to theith eigenvectorV _i;U _i= _iV_iT_i - U* _i ComponentU _i multiplied byc _i - U _ij Sum of componentsU _i+U _j - V _i Eigenvector of the matrixU;ith column ofV withUV=VL - w Weighting factor; equal to the ratio of two eigenvalues - X _i Random variable at pointi - x _i Value of random variable at pointi - y _i Residual ofx _i - Z _i Row vectori for the standardized variables - z _i Standardized value of variable     

Elasticity of TiO2 rutile to 1800 K

The ambient pressure elastic properties of single-crystal TiO₂ rutile are reported from room temperature (RT) to 1800 K, extending by more than 1200 oK the maximum temperature for which rutile elasticity data are available. The magnitudes of the temperature derivatives decrease with increasing temperature for five of the six adiabatic elastic moduli (C _ij ). At RT, we find (units, GPa): C ₁₁=268(1); C ₃₃=484(2); C ₄₄=123.8(2); C ₆₆=190.2(5); C ₂₃=147(1); and C ₁₂=175(1). The temperature derivatives (units, GPa K⁻¹) at RT are: (∂C ₁₁/∂T) _P =−0.042(5); (∂C ₃₃/∂T) _P =−0.087(6); (∂C ₄₄/∂T) _P =−0.0187(2); (∂C ₆₆/∂T) _P =−0.067(2); (∂C ₂₃/∂T) _P =−0.025; and (∂C ₁₂/∂T) _P −0.048(5). The values for K _S (adiabatic bulk modulus) and μ (isotropic shear modulus) and their temperature derivatives are K _S =212(1) GPa; μ=113(1) GPa; (∂K _S /∂T) _P =−0.040(4) GPa K⁻¹; and (∂μ/∂T) _P =−0.018(1) GPa K⁻¹. We calculate several dimensionless parameters over a large temperature range using our new data. The unusually high values for the Anderson-Gròneisen parameters at room temperature decrease with increasing temperature. At high T, however, these parameters are still well above those for most other oxides. We also find that for TiO₂, anharmonicity, as evidenced by a non-zero value of [∂ln (K _T )/∂lnV] _T , is insignificant at high T, implying that for the TiO₂ analogue of stishovite, thermal pressure is independent of volume (or pressure). Systematic relations indicate that ∂² K _S /∂T∂P is as high as 7×10⁻⁴ K⁻¹ for rutile, whereas ∂²μ/∂T∂P is an order of magnitude less. Received: 19 September 1997 / Revised, accepted: 27 February 1998     

Low-temperature thermodynamic properties of disordered zeolites of the natrolite group

 The heat capacity of paranatrolite and tetranatrolite with a disordered distribution of Al and Si atoms has been measured in the temperature range of 6–309 K using the adiabatic calorimetry technique. The composition of the samples is represented with the formula (Na_1.90K_0.22Ca_0.06)[Al_2.24Si_2.76O₁₀]·nH₂O, where n=3.10 for paranatrolite and n=2.31 for tetranatrolite. For both zeolites, thermodynamic functions (vibrational entropy, enthalpy, and free energy function) have been calculated. At T=298.15 K, the values of the heat capacity and entropy are 425.1 ± 0.8 and 419.1 ±0.8 J K⁻¹ mol⁻¹ for paranatrolite and 381.0 ± 0.7 and 383.2 ± 0.7 J K⁻¹ mol⁻¹ for tetranatrolite. Thermodynamic functions for tetranatrolite and paranatrolite with compositions corrected for the amount of extraframework cations and water molecules have also been calculated. The calculation for tetranatrolite with two water molecules and two extraframework cations per formula yields: C _p (298.15)=359.1 J K⁻¹ mol⁻¹, S(298.15) −S(0)=362.8 J K⁻¹ mol⁻¹. Comparing these values with the literature data for the (Al,Si)-ordered natrolite, we can conclude that the order in tetrahedral atoms does not affect the heat capacity. The analysis of derivatives dC/dT for natrolite, paranatrolite, and tetranatrolite has indicated that the water- cations subsystem within the highly hydrated zeolite may become unstable at temperatures above 200 K. Received: 30 July 2001 / Accepted: 15 November 2001