Determination of atrazine isotherms on agricultural soils

For determination of atrazine isotherms in agricultural soils of Fars Province, composite soil samples from 0 to 5 cm depth with textures of silty clay loam, clay loam and loam were collected. In order to form the atrazine isotherms, 10, 50 and 100 µg atrazine g^?1 soil was added to the soil samples. Soluble atrazine in water:soil ratios of 10:1, 50:1 and 200:1 was measured after 3-h shaking. Finally, for each cases of applied atrazine, water extractable atrazine was determined and quantified using gas chromatography instrument. The results indicated that there was a linear relationship between the logarithms of water extractable atrazine and added atrazine for different water:soil ratios. A general equation of WEA = K(WS) ^α (AA) ^β is obtained experimentally between water extractable atrazine, µg g^?1(WEA), and added atrazine, µg g^?1 (AA), where K, α and β are absorption constants; WS is the water:soil ratio, g g^?1. For the loam, silty clay loam and clay loam soil textures, the α were 0.49, 0.23 and 0.13, respectively, the β were 0.55, 0.806 and 0.21, respectively, and the K were 1.44, 0.78 and 25.38, respectively.     

Fractal features of soil particle-size distributions and their relationships with soil properties in gravel-mulched fields

Soil particle-size distribution (PSD) is an important index for soil classification because it has large influences on soil hydrological characteristics, salinity, fertility, erodibility, nutrient content, swelling/shrinking, and degradation. We present a case study of the fractal characteristics of soil PSD and its relationship with soil properties of gravel-mulched fields in an arid area of northwestern China using single-fractal calculation. Particle size was unimodally distributed within the narrow range of 20–100 μm, with silt as the most common component. Horizontally, silt content was the highest, followed by sand and clay contents. Vertically, clay content increased with depth, while there were no obvious change rules for both silt and sand contents. The volume fractal dimension (D) of PSD ranged from 2.4307 to 2.5260, increased with the content of fine particles but decreased with the content of coarse particles. D was correlated positively with soil-water content and salt content and negatively with bulk density. The saturated soil-water content was strongly correlated negatively with silt content (p < 0.01) and positively with sand content (p < 0.01). The results indicate that D can be a potential indicator of the physical and chemical properties of soil and can also provide a theoretical basis and technical guidance for the effective use and management of the region.     

Argon retention properties of silicate glasses and implications for <Superscript>40</Superscript>Ar/<Superscript>39</Superscript>Ar age and noble gas diffusion studies

Sized aggregates of glasses (47–84 wt% SiO₂) were fused from igneous-derived cohesive fault rock and igneous rock, and step-heated from ~400 to >1,200 °C to obtain their ³⁹Ar diffusion properties (average E=33,400 cal mol^?1; D _o=4.63×10^?3 cm² s^?1). At T<~1,000 °C, glasses containing <~69 wt% SiO₂ and abundant network-forming cations (Ca, Fe, Mg) reveal moderate to strong non-linear increases in D and E, reflecting structural modifications as the solid transitions to melt. Extrapolation of these Arrhenius properties down to typical geologic T-t conditions could result in a 1.5 log₁₀ unit underestimation in the diffusion rate of Ar in similar materials. Numerical simulations based upon the diffusion results caution that some common geologic glasses will likely yield ⁴⁰Ar/³⁹Ar cooling ages rather than formation ages. However, if cooling rates are sufficiently high, ambient temperatures are sufficiently low (e.g., <65–175 °C), and coarse particles (e.g., radius (r) >~1 mm) are analyzed, glasses with compositions similar to ours may preserve their formation ages.     

A modified normalized model for predicting effective soil thermal conductivity

Effective soil thermal conductivity (λ _eff) describes the ability of a multiphase soil to transmit heat by conduction under unit temperature gradient. It is a critical parameter for environmental science, earth and planetary science, and engineering applications. Numerous models are available in the literature, but their applicability is generally restricted to certain soil types or water contents (θ). The objective of this study was to develop a new model in the similar form of the Johansen 1975 model to simulate the λ _eff(θ) relationship of soils of various soil textures and water contents. An exponential type model with two parameters is developed and a new function for calculating dry soil thermal conductivity is presented. Performance of the new model and six other normalized models were evaluated with published datasets. The results show that the new model is able to well mimic λ _eff(θ) relationship of soils from sand to silt loam and from oven dry to full saturation. In addition, it has the best performance among the seven models under test (with root-mean-square error of 0.059 W m^?1 °C^?1, average deviations of 0.0009 W m^?1 °C^?1, and Nash–Sutcliffe efficiency of 0.994). The new model has potential to improve the reliability of soil thermal conductivity estimation and be incorporated into numerical modeling for environmental, earth and engineering studies.     

Hydraulic Conductivity of Fly Ash-Amended Mine Tailings

The objective of this study was to evaluate the effect of fly ash addition on hydraulic conductivity (k) of mine tailings. Mine tailings used in this study were categorized as synthetic tailings and natural tailings; two synthetic tailings were developed via blending commercially-available soils and natural tailings were collected from a garnet mine located in the U.S. Two fly ashes were used that had sufficient calcium oxide (CaO) content (17 and 18.9 %) to generate pozzolanic activity. Hydraulic conductivity was measured on pure tailings and fly ash-amended tailings in flexible-wall permeameters. Fly ash was added to mine tailings to constitute 10 % dry mass of the mixture, and specimens were cured for 7 and 28 days. The influence of fly ash-amendment on k of mine tailings was attributed to (1) molding water content and (2) plasticity of the mine tailings. Tailings that classified as low-plasticity silts with clay contents less than 15 % exhibited a decrease in k when amended with fly ash and prepared wet of optimum water content (w _opt ). Tailings that classified as low-plasticity clay exhibited a one-order magnitude increase in k with addition of fly ash for materials prepared dry or near w _opt . The decrease in k for silty tailings was attributed to formation of cementitious bonds that obstructed flow paths, whereas the increase in k for clayey tailings was attributed to agglomeration of clay particles and an overall increase in average pore size. The results also indicated that the effect of curing time on k is more pronounced during the early stages of curing (≤7 days), as there was negligible difference between k for 7 and 28-days cured specimens.     

Shear wave velocity as function of cone penetration resistance and grain size for Holocene-age uncemented soils: a new perspective

For feasibility studies and preliminary design estimates, field measurements of shear wave velocity, V _s, may not be economically adequate and empirical correlations between V _s and more available penetration measurements such as cone penetration test, CPT, data turn out to be potentially valuable at least for initial evaluation of the small-strain stiffness of soils. These types of correlations between geophysical (Vs) and geotechnical (N-SPT, q _c-CPT) measurements are also of utmost importance where a great precision in the calculation of the deposit response is required such as in liquefaction evaluation or earthquake ground response analyses. In this study, the stress-normalized shear wave velocity V _s1 (in m/s) is defined as statistical functions of the normalized dimensionless resistance, Q _tn-CPT, and the mean effective diameter, D ₅₀ (in mm), using a data set of different uncemented soils of Holocene age accumulated at various sites in North America, Europe, and Asia. The V _s1–Q _tn data exhibit different trends with respect to grain sizes. For soils with mean grain size (D ₅₀) < 0.2 mm, the V _s1/Q _tn ^0.25 ratio undergoes a significant reduction with the increase in D ₅₀ of the soil. This trend is completely reversed with further increase in D ₅₀ (D ₅₀ > 0.2 mm). These results corroborate earlier results that stressed the use of different CPT-based correlations with different soil types, and those emphasized the need to impose particle-size limits on the validity of the majority of available correlations.     

Criteria of Acceptance for Constant Rate of Strain Consolidation Test for Tropical Cohesive Soil

In this study, the rapid consolidation equipment (RACE) was developed as an alternative device to the conventional consolidation test using Oedometer, consuming merely a few hours for the whole precedure to determine the consolidation characteristics of cohesive soil. RACE operates based on the constant rate of strain (CRS) consolidation theory, which is a continuous loading method of testing, requiring a good estimation of the loading rate such that it is ideal for the achievement of steady state condition during testing. The steady state condition is achieved when the c _v values from drained and undrained face of CRS converged with the c_v from Oedometer test. A slightly modification has been made on the normal constant rate of strain (CRS) test by proposing a direct back pressure system to the specimen using a tube to saturate the soil sample. This research has produced a set of criteria for determining the suitable rate for the rapid consolidation test based on the ratio of normalized strain rate, β, and proposed a new coefficient in terms of a ratio of β to clay fraction (CF), as a part of new criteria for testing a fine soil. Four types of sample were tested with different rates of strain using the RACE and their results were compared with those conducted using the Oedometer on the same soil type, from which fairly good agreements were evident in many specimens. It was found from the study that the minimum value of normalized strain rate, β, for the CRS test is 0.005 and for the u _a /σ _v ratio is suggested as 0.01. Also, the maximum β/CF for soils with clay friction lower and higher than 50 % are 0.008 and 0.001, respectively. The minimum β/CF value for both conditions is 0.0001.     

Salinity Tolerance and Competition Drive Distributions of Native and Invasive Submerged Aquatic Vegetation in the Upper San Francisco Estuary

Both abiotic and biotic factors govern distributions of estuarine vegetation, and experiments can reveal effects of these drivers under current and future conditions. In upper San Francisco Estuary (SFE), increased salinity could result from sea level rise, levee failure, or water management. We used mesocosms to test salinity effects on, as well as competition between, the native Stuckenia pectinata (sago pondweed) and invasive Egeria densa (Brazilian waterweed), species with overlapping distributions at the freshwater transition in SFE. Grown alone at a salinity of 5, E. densa decreased fivefold in biomass relative to the freshwater treatment and decomposed within 3 weeks at higher salinities. In contrast, S. pectinata biomass accumulated greatly (~4× initial) at salinities of 0 and 5, doubled at 10, and was unchanged at 15. When grown together in freshwater, S. pectinata produced 75 % less biomass than in monoculture and significantly more nodal roots (suggesting increased nutrient foraging). At a salinity of 5, a decline in E. densa performance coincided with a doubling of S. pectinata shoot density. Additional experiments on E. densa showed elevated temperature (26 and 30 °C) suppressed growth especially at higher salinities (≥5). We conclude that salinity strongly influences distributions of both species and that competition from E. densa may impose limits on S. pectinata abundance in the fresher reaches of SFE. With a salinity increase of 5, S. pectinata is likely to maintain its current distribution while spreading up-estuary at the expense of E. densa, especially if increased temperature also reduces E. densa biomass.     

Assessment of textural variational pattern and electrical conduction of economic and accessible Quaternary hydrolithofacies via geoelectric and laboratory methods in SE Nigeria: A case study of select locations in Akwa Ibom State

Textural variational pattern of economic and accessible Quaternary aquifer repositories and its conductivity in the south-eastern Nigeria have been assessed through the integration of vertical electrical sounding and laboratory measurements. The results have shown the lithological attributes, pore-water and amount of residual clay minerals in the assumed clean sand; mechanism of charge fixation at the fluid - surface interface; intricate geometry of pores and pore channels; formation’s ability to transmit pore-water and cation exchange capacity.The connections of electrical and hydraulic properties and their distributions have been established. The average interface conductivity contributed by residual clay minerals in assumed clean sands of the aquifer repositories in the study area have been estimated as 30µS/m. Intrinsic average porosity and formation factor have been respectively deduced as 12% and 14.75. Comparing the simulated aquifer formation factor obtained from the observed porosity data with the observed aquifer formation factor, indicates the that study area has 0.5 ≤ a ≤ 0.8 pore geometry factor and 1.5 ≤ m ≤ 2.0 cementation factor as the best fitting values. The interrelations between aquifer parameters have been established through different plots and the aquifer have been empirically proved to be associated with residual clay minerals as the interface conductivity C_q is not equal to zero. The wide ranges of parameters estimated are an indication of variations in grain size. The estimated intrinsic average porosity, formation factor and the average BQ_v are viable in characterizing the aquifer flow dynamics and contaminant modelling in the associated aquifer sands For low pore geometry factors a (0.2) and low cementation factor m (0.5) the formation factor remains fairly constant. However, marked variability is noticed at higher a (1.0) and m (2.5). Despite the observed variability in formation factors at the indicated porosities, the spatial or geometrical spread of the formation factor remains unchanged in the aquifer units. The Tables for geoelectric and petrophysical parameters and the associated mathematical models generated in this study can be used for groundwater contaminant modelling and simulation of pore space parameters with reasonable accuracy.     

Stability at high-pressure,elastic behaviour and pressure-induced structural evolution of CsAlSi<Subscript>5</Subscript>O<Subscript>12</Subscript>, a potential host for nuclear waste

The elastic and structural behaviour of the synthetic zeolite CsAlSi₅O₁₂ (a = 16.753(4), b = 13.797(3) and c = 5.0235(17) Å, space group Ama2, Z = 2) were investigated up to 8.5 GPa by in situ single-crystal X-ray diffraction with a diamond anvil cell under hydrostatic conditions. No phase-transition occurs within the P-range investigated. Fitting the volume data with a third-order Birch–Murnaghan equation-of-state gives: V ₀ = 1,155(4) ų, K _T0 = 20(1) GPa and K′ = 6.5(7). The “axial moduli” were calculated with a third-order “linearized” BM-EoS, substituting the cube of the individual lattice parameter (a ³, b ³, c ³) for the volume. The refined axial-EoS parameters are: a ₀ = 16.701(44) Å, K _T0a = 14(2) GPa (β_a = 0.024(3) GPa^?1), K′ _a = 6.2(8) for the a-axis; b ₀ = 13.778(20) Å, K _T0b = 21(3) GPa (β_b = 0.016(2) GPa^?1), K′ _b = 10(2) for the b-axis; c ₀ = 5.018(7) Å, K _T0c = 33(3) GPa (β_c = 0.010(1) GPa^?1), K′ _c = 3.2(8) for the c-axis (K _T0a:K _T0b:K _T0c = 1:1.50:2.36). The HP-crystal structure evolution was studied on the basis of several structural refinements at different pressures: 0.0001 GPa (with crystal in DAC without any pressure medium), 1.58(3), 1.75(4), 1.94(6), 3.25(4), 4.69(5), 7.36(6), 8.45(5) and 0.0001 GPa (after decompression). The main deformation mechanisms at high-pressure are basically driven by tetrahedral tilting, the tetrahedra behaving as rigid-units. A change in the compressional mechanisms was observed at P ≤ 2 GPa. The P-induced structural rearrangement up to 8.5 GPa is completely reversible. The high thermo-elastic stability of CsAlSi₅O_12, the immobility of Cs at HT/HP-conditions, the preservation of crystallinity at least up to 8.5 GPa and 1,000°C in elastic regime and the extremely low leaching rate of Cs from CsAlSi₅O₁₂ allow to consider this open-framework silicate as functional material potentially usable for fixation and deposition of Cs radioisotopes.     

Tailings composition effects on shear strength behavior of co-mixed mine waste rock and tailings

The objective of this study was to evaluate the effect of mine tailings composition on shear behavior and shear strength of co-mixed mine waste rock and tailings (WR&T). Crushed gravel was used as a synthetic waste rock and mixed with four types of tailings: (1) fine-grained garnet, (2) coarse-grained garnet, (3) copper, and (4) soda ash. Co-mixed WR&T specimens were prepared to target mixture ratios of mass of waste rock to mass of tailings (R) such that tailings “just filled” interparticle void space of the waste rock (i.e., optimum mixture ratio, R _opt). Triaxial compression tests were conducted on waste rock, tailings, and mixed waste at effective confining stresses (\(\sigma_{\text{c}}^{{\prime }}\)) ranging from 5 to 40 kPa to represent stresses anticipated in final earthen covers for waste containment facilities. Waste rock and co-mixed WR&T specimens were 150 mm in diameter by 300 mm tall, whereas tailings specimens were 38 mm in diameter by 76 mm tall. Shear strength was quantified using effective stress friction angles (?′) from undrained tests: ?′ for waste rock was 37°, ?′ for tailings ranged from 34° to 41°, and ?′ for WR&T mixtures ranged from 38° to 40°. Thus, shear strength of co-mixed WR&T was comparable to waste rock regardless of tailings composition. Shear behavior of WR&T mixtures was a function of R and tailings composition. Tailings influenced shear behavior for R < R _opt and when tailings predominantly were silt. Shear behavior was influenced by waste rock for R ≥ R _opt and when tailings predominantly were sand or included clay particles.     

High-pressure displacive phase transition of a natural Mg-rich pigeonite

A high-pressure single-crystal X-ray diffraction study has been carried out on a P2₁/c natural Mg-rich pigeonite sample with composition ca. Wo₆En₇₆Fs₁₈ using a diamond anvil-cell. The unit-cell parameters were determined at 14 different pressures to 7.14 GPa. The sudden disappearance of the b-type reflections (h + k = odd) and a strong discontinuity (about 2.8%) in the unit-cell volume indicated a first-order P2₁/c–C2/c phase transition between 4.66 and 4.88 GPa. The P(V) data of the P2₁/c phase were fitted to 4.66 GPa by a third-order Birch–Murnaghan equation of state (BM3 EoS), whereas the limited number of experimental data collected within the C2/c phase between 4.88 and 7.14 GPa were fitted using the same equation of state but with K′ constrained to the value obtained for the P2₁/c fitting. The equation of state coefficients are V ₀ = 424.66(6) ų, K _T0 = 104(2) GPa and K′ = 8(1) for the P2₁/c phase, and V ₀ = 423.6(1) ų, K _T0 = 112.4(8) GPa, and K′ fixed to 8(1) for the C2/c phase. The axial moduli for a, b, and c for the P2₁/c phase were obtained using also a BM3-EoS, while for the C2/c phase only a linear calculation could be performed, and therefore the same approach was applied for comparison also to the P2₁/c phase. In general the C2/c phase exhibits axial compressibilities (β _c > β _a >> β _b) lower than those of the P2₁/c phase (β _b > β _c ≈ β _a; similar to those found in previous studies in clinopyroxenes and orthopyroxenes). The lower compressibility of the C2/c phase compared with that of the P2₁/c could be ascribed to the greater stiffness along the b direction. A previously published relationship between P _c and M2 average cation radius (i.r.) has been updated using all the literature data on P2₁/c clinopyroxene containing large cations at M2 site and our new data. The following weighted regression was obtained: P _c (GPa) = 26(4) ? 28(5) ×  i.r (Å), R ² = 0.97. This improved equation can be used to predict the critical pressure of natural P2₁/c clinopyroxene samples just knowing the composition at M2 site.     

Accuracy assessment of near-shore bathymetry information retrieved from Landsat-8 imagery

The improvement in the capabilities of Landsat-8 imagery to retrieve bathymetric information in shallow coastal waters was examined. Landsat-8 images have an additional band named coastal/aerosol, Band 1: 435–451 nm in comparison with former generation of Landsat imagery. The selected Landsat-8 operational land image (OLI) was of Chabahar Bay, located in the southern part of Iran (acquired on February 22, 2014 in calm weather and relatively low turbidity). Accurate and high resolution bathymetric data from the study area, produced by field surveys using a single beam echo-sounder, were selected for calibrating the models and validating the results. Three methods, including traditional linear and ratio transform techniques, as well as a novel proposed integrated method, were used to determine depth values. All possible combinations of the three bands [coastal/aerosol (CB), blue (B), and green (G)] have been considered (11 options) using the traditional linear and ratio transform techniques, together with five model options for the integrated method. The accuracy of each model was assessed by comparing the determined bathymetric information with field measured values. The standard error of the estimates, correlation coefficients (R ² ) for both calibration and validation points, and root mean square errors (RMSE) were calculated for all cases. When compared with the ratio transform method, the method employing linear transformation with a combination of CB, B, and G bands yielded more accurate results (standard error = 1.712 m, R ² _calibration = 0.594, R ² _validation = 0.551, and RMSE =1.80 m). Adding the CB band to the ratio transform methodology also dramatically increased the accuracy of the estimated depths, whereas this increment was not statistically significant when using the linear transform methodology. The integrated transform method in form of Depth = b ₀  + b ₁ X _CB  + b ₂ X _B  + b ₅ ln(R _CB )/ln(R _G ) + b ₆ ln(R _B )/ln(R _G ) yielded the highest accuracy (standard error = 1.634 m, R ² _calibration = 0.634, R ² _validation = 0.595, and RMSE = 1.71 m), where R _i (i = CB, B, or G) refers to atmospherically corrected reflectance values in the i ^th band [X _i  = ln(R _i -R _{deep water})].     

Geochemistry and origin of Dehoo manganese deposit,south Zahedan,southeastern Iran

Dehoo manganese deposit is located 52 km to the south of Zahedan in Sistan and Baluchestan Province, southeastern Iran. This deposit that lies in the central part of the Iranian Flysch Zone is lenticular in shape and lies above the micritic limestone-radiolarite cherts of the upper Cretaceous ophiolite unit. It is hosted within the reddish to brown radiolarite cherts and in places interlinks with them, so that the radiolarite chert packages play a key role for Mn mineralization in the region. Investigated ore-paragenetic successions and the geochemical characteristics of the Dehoo deposit were studied by means of major oxide, trace, and rare earth element (REE) contents that provide information as to the mineral origin. Strong positive correlations were found between major oxides and trace elements (Al₂O₃-TiO₂, r = 0.95; TiO₂-MgO, r = 0.94; Fe₂O₃-Al₂O₃, r = 0.90; MgO-Al₂O₃, r = 0.84; MgO-Fe₂O₃, r = 0.88; Fe₂O₃-TiO₂, r = 0.91; Fe₂O₃-K₂O, r = 0.74; Al₂O₃-K₂O, r = 0.69; Al₂O₃-V, r = 0.72; TiO₂-V, r = 0.73, and MgO-V, r = 0.69) that testify to the contribution of mafic terrigenous detrital material to the deposit. Chondrite-normalized REE patterns of all ore samples are characterized by negative Ce (0.06–0.15, average 0.10) and slightly positive Eu (0.29–0.45, average 0.36) anomalies. Based on ratios of Mn/Fe (average 56.23), Co/Ni (average 0.33), Co/Zn (average 0.38), U/Th (average 3.40), La/Ce (average 1.45), La_n/Nd_n (average 2.16), Dy_n/Yb_n (average 0.33), and light REE/heavy REE (average 8.40; LREE > HREE), as well as Ba (average 920 ppm) and total REE contents (average 6.96 ppm) negative Ce and positive Eu anomalies, Dehoo could be considered a predominantly submarine hydrothermal Mn deposit complemented by terrigenous detrital mafic material.     

Elastic behavior of vanadinite,Pb<Subscript>10</Subscript>(VO<Subscript>4</Subscript>)<Subscript>6</Subscript>Cl<Subscript>2</Subscript>, a microporous non-zeolitic mineral

The high-pressure behavior of a vanadinite (Pb₁₀(VO₄)₆Cl₂, a = b = 10.3254(5), c = 7.3450(4) Å, space group P6₃/m), a natural microporous mineral, has been investigated using in-situ HP-synchrotron X-ray powder diffraction up to 7.67 GPa with a diamond anvil cell under hydrostatic conditions. No phase transition has been observed within the pressure range investigated. Axial and volume isothermal Equations of State (EoS) of vanadinite were determined. Fitting the P–V data with a third-order Birch-Murnaghan (BM) EoS, using the data weighted by the uncertainties in P and V, we obtained: V ₀ = 681(1) Å³, K ₀ = 41(5) GPa, and K′ = 12.5(2.5). The evolution of the lattice constants with P shows a strong anisotropic compression pattern. The axial bulk moduli were calculated with a third-order “linearized” BM-EoS. The EoS parameters are: a ₀ = 10.3302(2) Å, K ₀(a) = 35(2) GPa and K′(a) = 10(1) for the a-axis; c ₀ = 7.3520(3) Å, K ₀(c) = 98(4) GPa, and K′(c) = 9(2) for the c-axis (K ₀(a):K ₀(c) = 1:2.80). Axial and volume Eulerian-finite strain (fe) at different normalized stress (Fe) were calculated. The weighted linear regression through the data points yields the following intercept values: Fe _a (0) = 35(2) GPa for the a-axis, Fe _c (0) = 98(4) GPa for the c-axis and Fe _V (0) = 45(2) GPa for the unit-cell volume. The slope of the regression lines gives rise to K′ values of 10(1) for the a-axis, 9(2) for the c-axis and 11(1) for the unit cell-volume. A comparison between the HP-elastic response of vanadinite and the iso-structural apatite is carried out. The possible reasons of the elastic anisotropy are discussed.     

Dewatering Behaviour of a Uranium Ore Slurry Containing Clays

The main objective of this paper was to investigate the dewatering behaviour of a clayey uranium ore slurry. The slurry (containing 28% clay size) exhibited moderate water adsorption (w _l  = 83% and w _p  = 30%). Primarily composed of muscovite (46%) and quartz (30%), the clay minerals included illite (8%), chlorite (5%) and kaolinite (2%) alongside a CEC of 41 (cmol(+)/kg) with Ca²⁺ and Mg²⁺ as the dominant cations. Likewise, the high EC (17,600 μS/cm) and ionic strength (1.15 mol/L) indicated a flocculated microstructure due to the presence of SO₄ ^2? (22,600 mg/L) and Mg²⁺ (1340 mg/L) in the slurry water. Settling included sedimentation and consolidation at low initial solids condition (25–35%) whereas only consolidation was observed at high initial solids contents (40–50%). The average k reduced from 1.2 × 10^?6 m/s (initial s = 25%) to 5.3 × 10^?8 m/s (initial s = 50%) along with a void ratio reduction from 7.4 to 2.6. Due to thixotropic strength, volume compressibility during consolidation showed apparent pre-consolidation at low effective stress (0.3–2 kPa) with a reduction in void ratio from 2.6 to 2.5. The e _s was found to be 2.46 at σ′ = 2 kPa and was followed by a steeper slope with the void ratio reducing to 2.1 at σ′ = 31 kPa. Likewise, the hydraulic conductivity during consolidation decreased from 2.6 × 10^?9 m/s (at e = 2.6) to 2.0 × 10^?10 m/s (at e = 2.1).     

A new generalized soil thermal conductivity model for sand–kaolin clay mixtures using thermo-time domain reflectometry probe test

As the demand of exploitation and utilization of geothermal energy increases, more geothermal-related earth structures occur recently. The design of the structures depends upon an accurate prediction of soil thermal conductivity. The existing soil thermal conductivity models were mostly developed by empirical fits to datasets of soil thermal conductivity measurements. Due to the gaps in measured thermal conductivities between any two tested natural soils, the models may not provide accurate prediction for other soils, and the predicted thermal conductivity might not be continuous over the entire range of soil type. In this research, a generalized soil thermal conductivity model was proposed based on a series of laboratory experiments on sand, kaolin clay and sand–kaolin clay mixtures using a newly designed thermo-time domain reflectometry probe. The model was then validated with respect to k _dry–n (thermal conductivity of dry soils and porosity) and k _r–S _r (normalized thermal conductivity and degree of saturation) relationships by comparing with previous experimental studies. The predicted thermal conductivities were found to be in a good agreement with the experimental data collected from both this study and the other literatures with at least 85% confidence interval. It is concluded that the proposed model accounts for the effects of both environmental factors (i.e., moisture content and dry density) and compositional factors (i.e., quartz content and soil type) on soil thermal conductivity, and it has a great potential in predicting soil thermal conductivity more accurately for geothermal applications.     

Compressibility of strontium orthophosphate Sr<Subscript>3</Subscript>(PO<Subscript>4</Subscript>)<Subscript>2</Subscript> at high pressure

High pressure in situ synchrotron X-ray diffraction experiment of strontium orthophosphate Sr₃(PO₄)₂ has been carried out to 20.0 GPa at room temperature using multianvil apparatus. Fitting a third-order Birch–Murnaghan equation of state to the P–V data yields a volume of V ₀ = 498.0 ± 0.1 Å³, an isothermal bulk modulus of K _T  = 89.5 ± 1.7 GPa, and first pressure derivative of K _T ′ = 6.57 ± 0.34. If K _T ′ is fixed at 4, K _T is obtained as 104.4 ± 1.2 GPa. Analysis of axial compressible modulus shows that the a-axis (K _a  = 79.6 ± 3.2 GPa) is more compressible than the c-axis (K _c  = 116.4 ± 4.3 GPa). Based on the high pressure Raman spectroscopic results, the mode Grüneisen parameters are determined and the average mode Grüneisen parameter of PO₄ vibrations of Sr₃(PO₄)₂ is calculated to be 0.30(2).     

Diffusive fractionation of volatiles and their isotopes during bubble growth in magmas

Bubbles grow in decompressing magmas by simple expansion and by diffusive supply of volatiles to the bubble/melt interface. The latter phenomenon is of significant geochemical interest because diffusion can fractionate elements and isotopes (or isotopologues) of dissolved components. This raises the possibility that the character of volatile components in bubbles may not reflect that of volatiles dissolved in the host melt over the lifetime of a bubble—even in the absence of equilibrium vapor/melt isotopic fractionation. Recent experiments have confirmed the existence of an isotope mass effect on diffusion of the volatile element Cl in silicate melt [Fortin et al. (Isotopic fractionation of chlorine during chemical diffusion in a dacitic melt and its implications for isotope behavior during bubble growth (abstract), 2016 Fall AGU Meeting, 2016)], so there is a clear need to understand the efficacy of diffusive fractionation during bubble growth. In this study, numerical models of diffusion and mass redistribution during bubble growth were implemented for both “passive” volatiles—those whose concentrations are generally well below saturation levels—and “active” volatiles such as CO₂ and H₂O, whose elevated concentrations and limited solubilities are the cause of bubble nucleation and growth. Both diffusive and convective bubble-growth scenarios were explored. The magnitude of the isotope mass effect on passive volatiles partitioned into bubbles growing at a constant rate R in a static system depends upon R/D _L, K _d and D _H/D _L (K _d = bubble/melt partition coefficient; D _H/D _L = diffusivity ratio of the heavy and light isotopes). During convective bubble growth, the presence of a discrete (physical) melt boundary layer against the growing bubble (of width x _BL) simplifies outcomes because it leads to the quick onset of steady-state fractionation during growth, the magnitude of which depends mainly upon R?x _BL/D _L and D _H/D _L (bubble/melt fractionation is maximized at R?x _BL/D _L ≈0.1). Constant R is unrealistic for most real systems, so other scenarios were explored by including the solubility and EOS of an “active” volatile (e.g., CO₂) in the numerical simulations. For plausible decompression paths, R increases exponentially with time—leading, potentially, to larger isotopic fractionation of species partitioned into the growing bubble. For volatile species whose isotope mass effects on diffusion have been measured (Cl, Li), predicted isotope fractionation in the exsolved vapor can be as large as ?4‰ for Cl and ?25‰ for Li.