The dilatant behaviour of sand–pile interface subjected to loading and stress relief

Property and behaviour of sand–pile interface are crucial to shaft resistance of piles. Dilation or contraction of the interface soil induces change in normal stress, which in turn influences the shear stress mobilised at the interface. Although previous studies have demonstrated this mechanism by laboratory tests and numerical simulations, the interface responses are not analysed systematically in terms of soil state (i.e. density and stress level). The objective of this study is to understand and quantify any increase in normal stress of different pile–soil interfaces when they are subjected to loading and stress relief. Distinct element modelling was carried out. Input parameters and modelling procedure were verified by experimental data from laboratory element tests. Parametric simulations of shearbox tests were conducted under the constant normal stiffness, constant normal load and constant volume boundary conditions. Key parameters including initial normal stress ( $ \sigma_{{{\text{n}}0}}^{\prime } $ ), initial void ratio (e ₀), normal stiffness constraining the interface and loading–unloading stress history were investigated. It is shown that mobilised stress ratio ( $ \tau /\sigma_{\text{n}}^{\prime } $ ) and normal stress increment ( $ \Updelta \sigma_{\text{n}}^{\prime } $ ) on a given interface are governed by $ \sigma_{{{\text{n}}0}}^{\prime } $ and e ₀. An increase in $ \sigma_{{{\text{n}}0}}^{\prime } $ from 100 to 400 kPa leads to a 30 % reduction in $ \Updelta \sigma_{\text{n}}^{\prime } $ . An increase in e ₀ from 0.18 to 0.30 reduces $ \Updelta \sigma_{\text{n}}^{\prime } $ by more than 90 %, and therefore, shaft resistance is much lower for piles in loose sands. A unique relationship between $ \Updelta \sigma_{\text{n}}^{\prime } $ and normal stiffness is established for different soil states. It can be applied to assess the shaft resistance of piles in soils with different densities and subjected to loading and stress relief. Fairly good agreement is obtained between the calculated shaft resistance based on the proposed relationship and the measured results in centrifuge model tests.     

Thermodynamic mixing properties of Mg(Al,Cr)2O4 spinel crystalline solution at high temperatures and pressures

Three Al-Cr exchange isotherms at 1,250°, 1,050°, and 796° between Mg(Al, Cr)₂O₄ spinel and (Al, Cr)₂O₃ corundum crystalline solutions have been studied experimentally at 25 kbar pressure. Starting from gels of suitable bulk compositions, close approach to equilibrium has been demonstrated in each case by time studies. Using the equation of state for (Al, Cr)₂O₃ crystalline solution (Chatterjee et al. 1982a) and assuming that the Mg(Al, Cr)₂O₄ can be treated in terms of the asymmetric Margules relation, the exchange isotherms were solved for Δ G ^*, and . The best constrained data set from the 1,250° C isotherm clearly shows that the latter two quantities do not overlap within three standard deviations, justifying the choice of asymmetric Margules relation for describing the excess mixing properties of Mg(Al, Cr)₂O₄ spinels. Based on these experiments, the following polybaric-polythermal equation of state can be formulated: , P expressed in bars, T in K, G _m ^ex and W _G,i ^Sp in joules/mol. Temperature-dependence of G _m ^ex is best constrained in the range 796–1,250° C; extrapolation beyond that range would have to be done with caution. Such extrapolation to lower temperature shows tentatively that at 1 bar pressure the critical temperature, T _c, of the spinel solvus is 427° C, with dT_c/dP≈1.3 K/kbar. The critical composition, X _c, is 0.42 , and changes barely with pressure. Substantial error in calculated phase diagrams will result if the significant positive deviation from ideality is ignored for Al-Cr mixing in such spinels.     

Spinel-pyroxene-garnet relationships and their dependence on Cr/Al ratio

The partitioning of Cr and Al between coexisting spinel and clinopyroxene and the dependence of spinel-cpxgarnet equilibria on Cr/Al ratio have been investigated by a combination of phase equilibrium experiments, high temperature solution calorimetry and thermodynamic calculations.The exchange equilibrium: has a measured enthalpy change for pure phases of –2,100±500 cal at 970 K and 1 atm. Experimental reversals of Cr-Al partitioning between the spinel and clinopyroxene phases yield the following partitioning relationship: where X _i ^j refers to atomic fraction of i in the octahedral sites of phase j. The compositional dependence of partitioning implies that Al-Cr mixing in spinel is nonideal with, on the symmetrical model, a W _Cr-Al ^Sp of 2,700±500 cal/gm. atom. In contrast, aluminum-chromium mixing in clinopyroxene is close to ideal.The measured stability field of knorringite (Mg₃Cr₂Si₂O₁₂) and mixing properties of garnet have been used in conjunction with our experimental data to calculate the influence of Cr/Al ratio on the important reaction: orthopyroxene+clinopyroxene+spinel=olivine+garnetThe stability field of spinel lherzolite increases by about 2.8 Kb for every increase of 0.1 in Cr/(Cr+Al) ratio up to Cr/(Cr+Al) of 0.7. The calculated stabilization is in very good agreement with the experimental results of O'Neill (1981). The partitioning relationships are such that, at the low ratios of Cr/Al (0.07) of primitive lherzolite, clinopyroxene buffers spinel composition and sharpens the spinelgarnet reaction interval from 10 Kb (little or no clinopyroxene) down to about 2 Kb in pyroxene-rich pyrolite.     

The notion of “distinguishability” between bulk elastic parameters on the basis of the Gibbs deformation energy

The present work aims in discussing a principle that distinguishes between elastic parameters sets, $ \{ \Upphi \} \equiv \{ K_{0} , \, K^{\prime}, \, V_{0} ,\ldots\} The present work aims in discussing a principle that distinguishes between elastic parameters sets, { \Upphi } o { K₀ ,  K^￠,  V₀ ,?} \{ \Upphi \} \equiv \{ K_{0} , \, K^{\prime}, \, V_{0} ,\ldots\} , on the basis of an energetic criterion: once a reference set, { \Upphi_R } \{ \Upphi_{R} \} , is given, another one can be fixed, { \Upphi _min } \left\{ {\Upphi_{ \min } } \right\} , so that they are as close as possible to each other, but yield non-equivalent deformation energy curves \Updelta G({ \Upphi } )_\textdeform \Updelta G(\{ \Upphi \} )_{\text{deform}} , i.e. they give \Updelta G({ \Upphi_R } )_\textdeform \Updelta G(\{ \Upphi_{R} \} )_{\text{deform}} and \Updelta G({ \Upphi _min } )_\textdeform \Updelta G(\{ \Upphi_{ \min } \} )_{\text{deform}} such that | \Updelta G({ \Upphi _min } )_\textdeform - \Updelta G({ \Upphi_R } )_\textdeform | 3 1×s[\Updelta G_\textdeform ]. \left| {\Updelta G(\{ \Upphi_{ \min } \} )_{\text{deform}} - \Updelta G(\{ \Upphi_{R} \} )_{\text{deform}} } \right| \ge 1\times \sigma [\Updelta G_{\text{deform}} ]. ΔG _deform, calculated using the equation of state (EoS), and its uncertainty σ[ΔG _deform], obtained by a propagation of the errors affecting { \Upphi } \{ \Upphi \} are crucial to fix which mineral assemblage forms at P–T conditions and allow one to assess the reliability of such a prediction. We explore some properties related to the principle introduced, using the average values of the elastic parameters found in literature and related uncertainties for di-octahedral mica, olivine, garnet and clinopyroxene. Two elementary applications are briefly discussed: the effect of refining V ₀ in fitting EoSs to P–V experimental data, in the case of garnet and omphacite, and the phengite 3T–2M ₁ relative stability, controlled by pressure.     

Interaction of Freshly Precipitated Silica Gel with Aqueous Silicic Acid Solutions under Ambient and Near Neutral pH-conditions: A Detailed Analysis of Linear Rate Law

Interaction of freshly precipitated silica gel with aqueous solutions was studied at laboratory batch experiments under ambient and near neutral pH-conditions. The overall process showed excellent reversibility: gel growth could be considered as an opposite process to dissolution and a linear rate law could be applied to experimental data. Depending on the used rate law form, the resulting rate constants were sensitive to errors in parameters/variables such as gel surface area, equilibrium constants, Si-fluxes, and reaction quotients. The application of an Integrated Exponential Model appeared to be the best approach for dissolution data evaluation. It yielded the rate constants k _dissol ∼ (4.50 ± 0.68) × 10⁻¹² and k _growth ∼ (2.58 ± 0.39) × 10⁻⁹ mol m⁻² s⁻¹ for zero ionic strength. In contrast, a Differential Model gave best results for growth data modeling. It yielded the rate constants k _dissol ∼ (1.14 ± 0.44) × 10⁻¹¹ and k _growth ∼ (6.08 ± 2.37) × 10⁻⁹ mol m⁻² s⁻¹ for higher ionic strength (I ∼ 0.04 to 0.11 mol L⁻¹). The found silica gel solubility at zero ionic strength was somewhat lower than the generally accepted value. Based on the and standard Gibbs free energy of silica gel formation was calculated as and −850,318 ± 20 J mol⁻¹, respectively. Activation energies for silica gel dissolution and growth were determined as and respectively. An universal value for growth of any silica polymorph, is not consistent with the value for silica gel growth, which questions the hypothesis about one unique activated complex controlling the silica polymorph growth.     

Estimation of covariance matrix from geochemical data with observations below detection limits

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.     

The Equation of State for Caspian Sea Waters

The density ρ of Caspian Sea waters was measured as a function of temperature (273.15–343.15) K at conductivity salinities of 7.8 and 11.3 using the Anton-Paar Densitometer. Measurements were also made on one of the samples (S = 11.38) diluted with water as a function of temperature (T = 273.15–338.15 K) and salinity (2.5–11.3). These latter results have been used to develop an equation of state for the Caspian Sea (σ = ±0.007 kg m⁻³)

The Activity of H2O,Fugacity of O2 and Variation of Their Chemical Potential During Metamorphic Peak of the Granulite Complex Near Taipingzhai—Louzishan Region,East Hebei Province

The thermodynamic calculation of dehydration reacton suggests very low activity of H2O during metamorphic peak of the Archaean granulite complex in the region studied.The αH2O values for Al-rich gneiss and hypersthene biotite gneiss-granulite in the Taipingzhai region are usually between 0.10 and 0.20,and those in the Louzishan region are 0.15-0.25.The fugacity of O2 in terms of lgf O2 in whole region ranges form-8to-14.The average coefficients of (δμH2O/δHMg^Bt)and(δμO2/δXMg^Bt)in the Taipingzhai region are-0.293 and-1.60 respectively,and those in the Louzishan region are-0.364and-1.420.The activity of H2O is very low in the whole region,but its values and other data mentioned above are considerably constant from place to place within a given region,even in rocks of dirrerent lithological characters.However,they show a certain gradient between different regions.Such characteristics are compatible with the genetic mechanism known as“carbonic metamorphism” put forward by Newton et al.,i.e.,the α H2O during the peak stage is controlled by permeation of pervasive CO2 influx of the mantle source,and shows features of external buffering.     

High-pressure transitions and thermochemistry of MGeO3 (M=Mg, Zn and Sr) and Sr-silicates: systematics in enthalpies of formation of A2+B4+O3 perovskites

Phase transitions in MgGeO₃ and ZnGeO₃ were examined up to 26 GPa and 2,073 K to determine ilmenite–perovskite transition boundaries. In both systems, the perovskite phases were converted to lithium niobate structure on release of pressure. The ilmenite–perovskite boundaries have negative slopes and are expressed as P(GPa)=38.4–0.0082T(K) and P(GPa)=27.4−0.0032T(K), respectively, for MgGeO₃ and ZnGeO₃. Enthalpies of SrGeO₃ polymorphs were measured by high-temperature calorimetry. The enthalpies of SrGeO₃ pseudowollasonite–walstromite and walstromite–perovskite transitions at 298 K were determined to be 6.0±8.6 and 48.9±5.8 kJ/mol, respectively. The calculated transition boundaries of SrGeO₃, using the measured enthalpy data, were consistent with the boundaries determined by previous high-pressure experiments. Enthalpy of formation (ΔH _f°) of SrGeO₃ perovskite from the constituent oxides at 298 K was determined to be −73.6±5.6 kJ/mol by calorimetric measurements. Thermodynamic analysis of the ilmenite–perovskite transition boundaries in MgGeO₃ and ZnGeO₃ and the boundary of formation of SrSiO₃ perovskite provided transition enthalpies that were used to estimate enthalpies of formation of the perovskites. The ΔH _f° of MgGeO₃, ZnGeO₃ and SrSiO₃ perovskites from constituent oxides were 10.2±4.5, 33.8±7.2 and −3.0±2.2 kJ/mol, respectively. The present data on enthalpies of formation of the above high-pressure perovskites were combined with published data for A²⁺B⁴⁺O₃ perovskites stable at both atmospheric and high pressures to explore the relationship between ΔH _f° and ionic radii of eightfold coordinated A²⁺ (R _A) and sixfold coordinated B⁴⁺ (R _B) cations. The results show that enthalpy of formation of A²⁺B⁴⁺O₃ perovskite increases with decreasing R _A and R _B. The relationship between the enthalpy of formation and tolerance factor ( R _o: O²⁻ radius) is not straightforward; however, a linear relationship was found between the enthalpy of formation and the sum of squares of deviations of A²⁺ and B⁴⁺ radii from ideal sizes in the perovskite structure. A diagram showing enthalpy of formation of perovskite as a function of A²⁺ and B⁴⁺ radii indicates a systematic change with equienthalpy curves. These relationships of ΔH _f° with R _A and R _B can be used to estimate enthalpies of formation of perovskites, which have not yet been synthesized.     

Experimental determination of activities in Fe−Mg olivine at 1400 K

The exchage equilibrium has been used to measure activity-composition relations along the olivine join FeSi_0.5O₂−MgSi_0.5O₂ at 1400 K and 1 atm pressure. Equilibrium Fe−Mg partitioning between the two phases was determined by reversing the compositions of olivine coexisting with oxide and matallic iron over the composition range Fo₂₃ to Fo₉₂. A detailed study of the thermodynamic properties of the oxide phase has recently been made by Srečec et al. and we have confirmed their results in the composition range of interest. Application of the oxide data to the exchange equilibrium enables the properties of olivine to be determined. Within experimental uncertainly (Fe, Mg)Si_0.5O₂ olivine can, at 1400 K, be treated as a symmetric solution with W _Fe-Mg ^o1 of 3.7±0.8 kJ/mol. The data permit the presence of only very slight asymmetry in the series. The data do not support recent assertions that olivine is highly non-ideal (W≈10 kJ/mol) under these conditions.     

Influence of hydrogen on Fe–Mg interdiffusion in (Mg,Fe)O and implications for Earth’s lower mantle

Interdiffusion of Fe and Mg in (Mg,Fe)O has been investigated experimentally under hydrous conditions. Single crystals of MgO in contact with (Mg_0.73Fe_0.27)O were annealed hydrothermally at 300 MPa between 1,000 and 1,250°C and using a Ni–NiO buffer. After electron microprobe analyses, the dependence of the interdiffusivity on Fe concentration was determined using a Boltzmann–Matano analysis. For a water fugacity of ∼300 MPa, the Fe–Mg interdiffusion coefficient in Fe _x Mg_1−x O with 0.01 ≤ x ≤ 0.25 can be described by with and C = −80 ± 10 kJ mol⁻¹. For x = 0.1 and at 1,000°C, Fe–Mg interdiffusion is a factor of ∼4 faster under hydrous than under anhydrous conditions. This enhanced rate of interdiffusion is attributed to an increased concentration of metal vacancies resulting from the incorporation of hydrogen. Such water-induced enhancement of kinetics may have important implications for the rheological properties of the lower mantle.

Transition enthalpies and entropies of high pressure zinc metasilicates and zinc metagermanates

The enthalpies of transition at T= 298 K between zinc metasilicate assemblages, measured by molten oxide solution calorimetry, are:

Opening and resetting temperatures in heating geochronological systems

We present a theoretical model for diffusive daughter isotope loss in radiochronological systems with increasing temperature. It complements previous thermochronological models, which focused on cooling, and allows for testing opening and resetting of radiochronometers during heating. The opening and resetting temperatures are, respectively,

Experimental and computational study of trace element distribution between orthopyroxene and anhydrous silicate melt: substitution mechanisms and the effect of iron

Although orthopyroxene (Opx) is present during a wide range of magmatic differentiation processes in the terrestrial and lunar mantle, its effect on melt trace element contents is not well quantified. We present results of a combined experimental and computational study of trace element partitioning between Opx and anhydrous silicate melts. Experiments were performed in air at atmospheric pressure and temperatures ranging from 1,326 to 1,420°C in the system CaO–MgO–Al₂O₃–SiO₂ and subsystem CaO–MgO–SiO₂. We provide experimental partition coefficients for a wide range of trace elements (large ion lithophile: Li, Be, B, K, Rb, Sr, Cs, Ba, Th, U; rare earth elements, REE: La, Ce, Nd, Sm, Y, Yb, Lu; high field strength: Zr, Nb, Hf, Ta, Ti; transition metals: Sc, V, Cr, Co) for use in petrogenetic modelling. REE partition coefficients increase from $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 Although orthopyroxene (Opx) is present during a wide range of magmatic differentiation processes in the terrestrial and lunar mantle, its effect on melt trace element contents is not well quantified. We present results of a combined experimental and computational study of trace element partitioning between Opx and anhydrous silicate melts. Experiments were performed in air at atmospheric pressure and temperatures ranging from 1,326 to 1,420°C in the system CaO–MgO–Al₂O₃–SiO₂ and subsystem CaO–MgO–SiO₂. We provide experimental partition coefficients for a wide range of trace elements (large ion lithophile: Li, Be, B, K, Rb, Sr, Cs, Ba, Th, U; rare earth elements, REE: La, Ce, Nd, Sm, Y, Yb, Lu; high field strength: Zr, Nb, Hf, Ta, Ti; transition metals: Sc, V, Cr, Co) for use in petrogenetic modelling. REE partition coefficients increase from $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 to $ D_{\text{Lu}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.109 $ D_{\text{Lu}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.109 , D values for highly charged elements vary from $ D_{\text{Th}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0026 $ D_{\text{Th}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0026 through $ D_{\text{Nb}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0033 $ D_{\text{Nb}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0033 and $ D_{\text{U}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0066 $ D_{\text{U}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0066 to $ D_{\text{Ti}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.058 $ D_{\text{Ti}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.058 , and are all virtually independent of temperature. Cr and Co are the only compatible trace elements at the studied conditions. To elucidate charge-balancing mechanisms for incorporation of REE into Opx and to assess the possible influence of Fe on Opx-melt partitioning, we compare our experimental results with computer simulations. In these simulations, we examine major and minor trace element incorporation into the end-members enstatite (Mg₂Si₂O₆) and ferrosilite (Fe₂Si₂O₆). Calculated solution energies show that R²⁺ cations are more soluble in Opx than R³⁺ cations of similar size, consistent with experimental partitioning data. In addition, simulations show charge balancing of R³⁺ cations by coupled substitution with Li⁺ on the M1 site that is energetically favoured over coupled substitution involving Al–Si exchange on the tetrahedrally coordinated site. We derived best-fit values for ideal ionic radii r ₀, maximum partition coefficients D ₀, and apparent Young’s moduli E for substitutions onto the Opx M1 and M2 sites. Experimental r ₀ values for R³⁺ substitutions are 0.66–0.67 ? for M1 and 0.82–0.87 ? for M2. Simulations for enstatite result in r ₀ = 0.71–0.73 ? for M1 and ~0.79–0.87 ? for M2. Ferrosilite r ₀ values are systematically larger by ~0.05 ? for both M1 and M2. The latter is opposite to experimental literature data, which appear to show a slight decrease in $ r_{0}^{{{\text{M}}2}} $ r_{0}^{{{\text{M}}2}} in the presence of Fe. Additional systematic studies in Fe-bearing systems are required to resolve this inconsistency and to develop predictive Opx-melt partitioning models for use in terrestrial and lunar magmatic differentiation models.     

Crystal field and covalency of octahedral chromium in natural [8](Mg1?xCax) 3[6] (Al0.67Cr0.33)2Si3O12 garnets from upper mantle rocks

In the course of a thorough study of the influences of the second coordination sphere on the crystal field parameters of the 3d ^N -ions and the character of 3d ^N –O bonds in oxygen based minerals, 19 natural Cr³⁺-bearing (Mg,Ca)-garnets from upper mantle rocks were analysed and studied by electronic absorption spectroscopy, EAS. The garnets had compositions with populations of the ^[8] X-sites by 0.881 ± 0.053 (Ca + Mg) and changing Ca-fractions in the range 0.020 ≤ w _Ca[8] ≤ 0.745, while the ^[6] Y-site fraction was constant with x _Cr³⁺ _[6] = 0.335 ± 0.023. The garnets had colours from deeply violet-red for low Ca-contents (up to x _Ca = 0.28), grey with 0.28 ≤ x _Ca ≤ 0.4 and green with 0.4 ≤ x _Ca. The crystal field parameter of octahedral Cr³⁺ 10Dq decreases strongly on increasing Ca-fraction from 17,850 cm⁻¹ at x _Ca[8] = 0.020 to 16,580 cm⁻¹ at x _Ca[8] = 0.745. The data could be fit with two model which do statistically not differ: (1) two linear functions with a discontinuity close to x _Ca[8] ≈ 0.3,

Spatial relationships of multivariate data

Spatial factor analysis (SFA) is a multivariate method that determines linear combinations of variables with maximum autocorrelation at a given lag. This is achieved by deriving estimates of auto-/cross-correlations of the variables and calculating the corresponding eigenvectors of the covariance quotient matrix. A two-point spatial factor analysis model derives factors by the formation of transition matrixU comparing auto-/cross-correlations at lag 0,R ₀, with those at a specified lag d,R _d, expressed asU _d=R ₀ ^–1 R_d. The matrixU _d can be decomposed into its spectral components which represent the spatial factors. The technique has been extended to include three points of reference. Spatial factors can be derived from the relationship:

On the problem of transport species of tungsten by hydrothermal solutions

The solubility of pentatungstate of sodium (PTS) Na₂W₅O₁₆ · H₂O and sodium tungsten bronzes (STB) Na_0.16WO₃ in acid chloride solutions containing 0.026, 0.26, and 3.02m NaCl have been studied at 500°C, 1000 bar, given fO₂ (Co-CoO, Ni-NiO, PTS-STB buffers), and constant NaCl/HCl ratio (Ta₂O₅-Na₂Ta₄O₁₁ buffer). Depending on experimental conditions, the tungsten content in the solutions after experiments varied from 10⁻³ to 2 × 10⁻² mol/kg H₂O. Obtained data were used to calculate the formation constants of predominant tungsten complexes (VI, V): H₃W₃^VIO₁₂³⁻, W₃^VO₉³⁻, [W^VW₄^VIO₁₆]³⁻, for reactions