Spatial control of groundwater contamination,using principal component analysis

A study on the geochemistry of groundwater was carried out in a river basin of Andhra Pradesh to probe into the spatial controlling processes of groundwater contamination, using principal component analysis (PCA). The PCA transforms the chemical variables, pH, EC, Ca²⁺, Mg²⁺, Na⁺, K⁺, HCO \(_3^- \) , Cl^?, SO \(_4^{2-} \) , NO \(_3^-\) and F^?, into two orthogonal principal components (PC1 and PC2), accounting for 75% of the total variance of the data matrix. PC1 has high positive loadings of EC, Na⁺, Cl^?, SO \(_4^{2-} \) , Mg²⁺ and Ca²⁺, representing a salinity controlled process of geogenic (mineral dissolution, ion exchange, and evaporation), anthropogenic (agricultural activities and domestic wastewaters), and marine (marine clay) origin. The PC2 loadings are highly positive for HCO \(_3^- \) , F^?, pH and NO \(_3^- \) , attributing to the alkalinity and pollution controlled processes of geogenic and anthropogenic origins. The PC scores reflect the change of groundwater quality of geogenic origin from upstream to downstream area with an increase in concentration of chemical variables, which is due to anthropogenic and marine origins with varying topography, soil type, depth of water levels, and water usage. Thus, the groundwater quality shows a variation of chemical facies from Na⁺ > Ca²⁺ > Mg²⁺ > K⁺: HCO \(_3^- \) > Cl^? > SO \(_4^{2-}>\) NO \(_3^- \) > F^?at high topography to Na⁺ > Mg²⁺ > Ca²⁺ > K⁺: Cl^? > HCO \(_3^- \) > SO \(_4^{2-}>\) NO \(_3^- \) > F^? at low topography. With PCA, an effective tool for the spatial controlling processes of groundwater contamination, a subset of explored wells is indexed for continuous monitoring to optimize the expensive effort.     

The calcium oxotellurate(IV) nitrates Ca5Te4O12(NO3)2(H2O)2 and Ca6Te5O15(NO3)2

Single crystals of two novel calcium oxotellurate(IV) nitrates were grown under hydrothermal conditions and were structurally characterized by X-ray diffraction. Ca $_5$ Te $_4\text {O}_{12}$ (NO $_3$ ) $_2$ (H $_2$ O) $_2$ [ $Cc$ , $Z=4$ , $a=25.258(3)$ Å, $b=5.7289(7)$ Å, $c=17.0066(19)$ Å, $\beta =124.377(2)^{\circ}$ , $R[F^2 > 2\sigma (F^2)]=0.043$ , 4083 $F^2$ data, 281 parameters] can be described as a non-classic order/disorder (OD) structure, which fulfills the basic principle of OD theory, viz. local equivalence of polytypes, but does not strictly follow the vicinity condition (VC) of OD theory. The structure is made up from an alternating stacking of non-polar layers composed of isolated [TeO $_3$ ] units and Ca $^{2+}$ ions and polar layers containing NO $_3^-$ ions and water molecules. The electron lone-pairs of the [TeO $_3$ ] units protrude into the free space of the anion/water layers. The crystal under investigation was a non-classic OD-twin of domains of a maximum degree of order (MDO). At the twin plane a fragment of the second MDO polytype is located. The main building blocks of Ca $_6$ Te $_5\text {O}_{15}$ (NO $_3$ ) $_2$ [ $P2_1/c$ , $Z=4$ , $a=15.494(2)$ Å, $b=5.6145(7)$ Å, $c=39.338(4)$ Å, $\beta =142.480(5)^{\circ}$ , $R[F^2 > 2\sigma (F^2)]=0.043$ , 3026 $F^2$ data, 307 parameters] are isolated [TeO $_3$ ] units and Ca $^{2+}$ ions which are connected to a three-dimensional framework perforated by channels in which the N atoms of the nitrate anions are located and the electron lone-pairs of the [TeO $_3$ ] units protrude. The structure can topologically be derived from the structure of Ca $_5$ Te $_4\text {O}_{12}$ (NO $_3$ ) $_2$ (H $_2$ O) $_2$ by removing the water molecules and connecting the CaTeO $_3$ layers with additional [TeO $_3$ ] units and Ca $^{2+}$ ions.     

Crystallographic preferred orientation in wüstite (FeO) through the cubic-to-rhombohedral phase transition

Magnesiowüstite, (Mg_0.08Fe_0.88)O, and wüstite, Fe_0.94O, were compressed to ~36?GPa at ambient temperature in the diamond anvil cell (DAC) at the Advanced Light Source. X-ray diffraction patterns were taken in situ in radial geometry in order to study the evolution of crystallographic preferred orientation through the cubic-to-rhombohedral phase transition. Under uniaxial stress in the DAC, {100}_c planes aligned perpendicular to the compression direction. The {100}_c in cubic became { $\left\{ {10\bar 14} \right\}$ }_r in rhombohedral and remained aligned perpendicular to the compression direction. However, the {101}_c and {111}_c planes in the cubic phase split into { ${10{\bar{1}}4}$ }_r and { ${11{\bar{2}}0}$ }_r, and (0001)_r and { ${10{\bar{1}}1}$ }_r, respectively, in the rhombohedral phase. The { ${11{\bar{2}}0}$ }_r planes preferentially aligned perpendicular to the compression direction while { ${10{\bar{1}}4}$ }_r oriented at a low angle to the compression direction. Similarly, { ${10{\bar{1}}1}$ }_r showed a slight preference to align more closely perpendicular to the compression direction than (0001)_r. This variant selection may occur because the 〈 ${10{\bar{1}}4}$ 〉_r and [0001]_r directions are the softer of the two sets of directions. The rhombohedral texture distortion may also be due to subsequent deformation. Indeed, polycrystal plasticity simulations indicate that for preferred { ${10{\bar{1}}4}$ }〈 ${1{\bar{2}}10}$ 〉_r and { ${11{\bar{2}}0}$ }〈 ${{\bar{1}}101}$ 〉_r slip and slightly less active { ${10{\bar{1}}1}$ }〈 ${{\bar{1}}2{\bar{1}}0}$ 〉_r slip, the observed texture pattern can be obtained.     

Experimental delineation of the “e” ⇌ I\bar 1 and “e” ⇌ C\bar 1 transformations in intermediate plagioclase feldspars

Approximately 125 hydrothermal annealing experiments have been carried out in an attempt to bracket the stability fields of different ordered structures within the plagioclase feldspar solid solution. Natural crystals were used for the experiments and were subjected to temperatures of ～650°C to ～1,000°C for times of up to 370 days at \(P_{{\text{H}}_{\text{2}} {\text{O}}} \) =600 bars, or \(P_{{\text{H}}_{\text{2}} {\text{O}}} \) =1,200 bars. The structural states of both parent and product materials were characterised by electron diffraction, with special attention being paid to the nature of type e and type b reflections (at h+k=(2n+1), l=(2n+1) positions). Structural changes of the type C \(\bar 1\) → I \(\bar 1\) , C \(\bar 1\) → “e” structure, I \(\bar 1\) → “e” and “e” structure → I \(\bar 1\) have been followed. There are marked differences between the ordering behaviour of crystals with compositions on either side of the C \(\bar 1\) ? I \(\bar 1\) transition line. In the composition range ～ An₅₀ to ～ An₇₀ the e structure appears to have a true field of stability relative to I \(\bar 1\) ordering, and a transformation of the type I \(\bar 1\) ? e has been reversed. It is suggested that the e structure is the more stable ordered state at temperatures of ～ 800°C and below. For compositions more albite-rich than ～ An₅₀ the upper temperature limit for long range e ordering is lower than ～ 750°C, and there is no evidence for any I \(\bar 1\) ordering. The evidence for a true stability field for “e” plagioclase, which is also consistent with calorimetric data, necessitates reanalysis both of the ordering behaviour of plagioclase crystals in nature and of the equilibrium phase diagram for the albite-anorthite system. Igneous crystals with compositions of ～ An₆₅, for example, probably follow a sequence of structural states C \(\bar 1\) → I \(\bar 1\) → e during cooling. The peristerite, Bøggild and Huttenlocher miscibility gaps are clearly associated with breaks in the albite, e and I \(\bar 1\) ordering behaviour but their exact topologies will depend on the thermodynamic character of the order/disorder transformations.     

Interdiffusion of NaSi—CaAl in peristerite

The ‘average’ interdiffusion coefficient ( \(\bar D\) ) for NaSi—CaAl exchange in plagioclase for the interval from An₀ to An₂₆ was estimated from experimentally determined homogenization times for peristerite exsolution lamellae. The average spacing between adjacent (unlike) lamellae is 554±77 Å. Dry heating in air at 1,100°C for 98 days produced no change in the exsolution microstructure; thus \(\bar D\) (dry)<10^?17 cm²/s. This limit is consistent with the recently reported ‘average’ \(\bar D\) (dry) values for the Huttenlocher interval (An_70–90) at this temperature. At 1.5 GPa with about 0.2 weight percent water added the ‘average’ diffusion coefficient from 1,100°C to 900°C is given by: \(\bar D\) (wet)=18 _?15 ⁺¹⁰⁸ (cm²/s) exp (?97±5 (kcal/mol)/RT), where R is the gas constant, and T is °K. This \(\bar D\) (wet) at 1,100°C is more than three orders of magnitude greater than \(\bar D\) (dry) for Na- and Ca-rich plagioclases.     

Analysis of the Influence of a Natural Fracture Network on Hydraulic Fracture Propagation in Carbonate Formations

A new experimental model has been designed to simulate the influence of a natural fracture network on the propagation geometry of hydraulic fractures in naturally fractured formations using a tri-axial fracturing system. In this model, a parallel and symmetrical pre-fracture network was created by placing cement plates in a cubic mold and filling the mold with additional cement to create the final testing block. The surface of the plates will thus be weakly cemented and form pre-fractures. The dimension and direction of the pre-fractures can be controlled using the plates. The experiments showed that the horizontal differential stress $\Updelta \sigma$ and the angle $\Updelta \theta$ between the maximum horizontal principal in situ stress and the pre-fracture are the dominating factors for the initiation and propagation of hydraulic fractures. For $\Updelta \theta = 90^\circ$ and $\Updelta \sigma \ge 2{\text{ MPa}}$ or $\Updelta \theta = 60^\circ$ and $\Updelta \sigma \ge 4{\text{ MPa}}$ , the direction of the initiation and propagation of the hydraulic fractures are consistent with or deviate from the normal direction of the pre-fracture. When the hydraulic fractures approach the pre-fractures, the direction of the hydraulic fracture propagation will be consistent with the normal direction of the pre-fracture. Otherwise, the hydraulic fracture will deflect and perpendicularly cross the parallel and symmetric pre-fracture network. For $\Updelta \theta = 90^\circ$ and $\Updelta \sigma < 2{\text{ MPa}},\,\Updelta \theta = 60^\circ$ , and $\Updelta \sigma < 4{\text{ MPa}}$ or $\Updelta \theta = 45^\circ$ and $\Updelta \sigma = 4 - 8{\text{ MPa}}$ , before the hydraulic fracture and the pre-fractures intersect, the direction of the hydraulic fracture propagation remains unchanged, and the pre-fractures open or dilate when the hydraulic fracture propagates to the intersection point, forming a complicated hydraulic fracture network with the propagation region of the overall hydraulic fracture network taking the shape of an ellipse. In this condition, the complexity level of the hydraulic fracture is controlled by the net pressure, the compressive normal stress acting on the pre-fractures, the shearing strength and the cohesion strength of the planes of weakness. The conclusions of this research are inconsistent with the formulation of the approach angle that has been widely accepted by previous studies. The principle of hydraulic fracture propagation is that it follows the least resistance, the most preferential propagation, and the shortest propagation path.     

Coherent elastic energies of exsolution boundaries in peristerite and böggild intergrowths

Elastic energy calculations based upon the coherent model of Willaime and Brown [Acta Cryst. A30, 316–331 (1974)] have been carried out for some exsolution textures in peristerite and Böggild intergrowths. For peristerites it is demonstrated that substitution of K for Na moves the orientation of the exsolution lamellae from (08 \(\overline {\text{1}} \) ) to about ( \(\overline {\text{1}} \) , 21, \(\overline {\text{2}} \) ) in agreement with observations. An electron microscope study of exsolution textures in labradorite and andesine plagioclases has been carried out and information about small structural differences in these lamellar structures has been obtained from elastic energy calculations.     

P wave attenuation structure beneath the Kanto district and surrounding area, Japan

Spectral ratios of teleseismic P waves for 15 deep (>200 km) earthquakes recorded at 146 High-Sensitivity Seismographic Network stations in the Kanto district and its surrounding area, eastern Japan, were inverted for attenuation parameter $ t_P^{ * } $ . The dataset consisted of good-quality vertical-component seismograms, whose P phases were handpicked. The P wave spectral ratios with high signal-to-noise ratios were calculated up to 1 Hz for all the station pairs, linear regressed, and then inverted for $ t_P^{ * } $ using the technique of least squares . The result showed that the active volcanic areas were clearly characterized by high $ t_P^{ * } $ . In contrast, $ t_P^{ * } $ varied in the nonvolcanic areas. The present result on the $ t_P^{ * } $ distribution was roughly consistent with the shallow part (<30 km) of an attenuation structure, which has been previously obtained based on 3-D tomography by using records of high-frequency (around 5 Hz) P waves from local earthquakes. This suggested that the present method of $ t_P^{ * } $ estimation is valid. The advantage and possible application to other areas were also discussed.     

Experimental Investigation of Seepage Properties of Fractured Rocks Under Different Confining Pressures

The effectiveness of transmitting underground water in rock fractures is strongly influenced by the widths of the fractures and their interconnections. However, the geometries needed for water flow in fractured rock are also heavily controlled by the confining pressure conditions. This paper is intended to study the seepage properties of fractured rocks under different confining pressures. In order to do this, we designed and manufactured a water flow apparatus that can be connected to the electro-hydraulic servo-controlled test system MTS815.02, which provides loading and exhibits external pressures in the test. Using this apparatus, we tested fractured mudstone, limestone and sandstone specimens and obtained the relationship between seepage properties and variations in confining pressure. The calculation of the seepage properties based on the collection of water flow and confining pressure differences is specifically influenced by non-Darcy flow. The results show that: (1) The seepage properties of fractured rocks are related to confining pressure, i.e. with the increase of confining pressure, the permeability $ k $ decreases and the absolute value of non-Darcy flow coefficient $ \beta $ increases. (2) The sandstone coefficients $ k $ and $ \beta $ range from $ 1.03 \times 10^{ - 18} $ to $ 1.53 \times 10^{ - 17} $  m² and $ - 1.13 \times 10^{17} $ to $ - 2.35 \times 10^{18} $  m^?1, respectively, and exhibit a greater change compared to coefficients of mudstone and limestone. (3) From the regression analysis of experimental data, it is concluded that the polynomial function is a better fit than the power and logarithmic functions. The results obtained can provide an important reference for understanding the stability of rock surrounding roadways toward prevention of underground water gushing-out, and for developing underground resources (e.g. coal).     

Schiller effects and exsolution in sodium-rich plagioclases

Plagioclase feldspars with mean compositions Ab_91,3Or_4,7An_4,0 and Ab_88,7An_10,1Or_1,2 have been studied by transmission electron microscopy and electron diffraction. The substructure consists of thin lamellae of albite and oligoclase. Two types of orientations of the lamellar planes were observed. The orientation of the more common type was found to change from (08 \(\bar 1\) ) to about ( \(\bar 1\) , 21, \(\bar 2\) ) as a function of the mean potassium content. The plane of the other type was found to be near ( \(\bar 7\) 12). Only the first type of lamellae produces visible Schiller colours.     

High-temperature deformation of forsterite single crystals doped with vanadium

Creep experiments have been performed on samples from a single crystal of vanadium-doped forsterite under controlled \(p_{{\text{O}}_2 } \) conditions to investigate the effects of the addition of substitutional defects in the tetrahedral lattice sites. The addition of vanadium causes marked changes in the flow behavior of the forsterite, with a net increase in the creep rate at high \(p_{{\text{O}}_2 } \) and a new \(p_{{\text{O}}_2 } \) -dependent flow regime at low \(p_{{\text{O}}_2 } \) conditions. These observations can be interpreted as resulting from changes in the majority defect species that maintain the charge neutrality within the crystal. A climb-controlled dislocation creep model for the high-temperature deformation of vanadium-doped forsterite is proposed in which either (i) movement of uncharged jogs is rate-limited by the diffusion of silicon via a vacancy mechanism or (ii) movement of positively charged jogs is rate-limited by diffusion of oxygen via a vacancy mechanism.     

The order-disorder transition in gallium analogue of albite,NaGaSi3O8

The order-disorder transformation in NaGaSi₃O₈ was studied at \(P_{H_2 O} \approx 100\) bar using a hydrous sodium gallium silicate gel and synthetic Ga-albites of various degrees of order as initial materials. The structural state of the gallium albites was determined by the angle separation in X-ray powder patterns of 131 and 1 \(\bar 3\) 1, and of \(\bar 2\) 04 and 060 peaks. The direction of order-disorder transformations was found to change sharply at 938 (±3)°C. This indicates a first order phase transition in NaGaSi₃O₈ as in albite, of which it is the structural analogue.     

Standing stock and estimated production rates of phytoplankton and zooplankton in Narragansett Bay,Rhode Island

Seasonal changes in phytoplankton biomass and production, total zooplankton biomass, and biomass and potential production rates of the two dominant copepods, Acartia hudsonica (formerly called Acartia clausi) and Acartia tonsa are described for several stations in Narragansett Bay, R.I. Plankton in the bay behaved as a single population with simultaneous changes occurring at the upper bay (Station 5) and the lower bay (Station 1). Phytoplankton biomass was higher in the upper bay ( \(\bar x\) =16.95 mg chl a·m^?3) than in the lower bay ( \(\bar x\) =6.37 mg chl a·m^?3) and these 0269 0101 V differences in biomass were reflected in the phytoplankton production rates. The zooplankton, which was dominated by A. hudsonica in the spring and early summer and A. tonsa during summer and fall, showed no such consistent differences between the stations. Mean A. hudsonica biomass (St 1, \(\bar x\) ;=82.7 mg dry wt·m^?3; St 5, _ \(\bar x\) ;=95.2 mg dry wt·m^?3) exceeded that of A. tonsa (St 1, \(\bar x\) ;=56.7 mg dry wt·m^?3; St 5, \(\bar x\) ;=60.0 mg dry wt·m^?3). Potential production rates of the two Acartia 0269 0101 V spp. were strongly temperature dependent. Despite the higher biomass levels of A. hudsonica, low temperatures resulted in lower potential production rates ( \(\bar x\) ; St 1=7.25 mg C·m^?3 day^?1; \(\bar x\) ; St 5=10.77mg C·m^?3 day^?1) and biomass doubling times of up to 9.6 days. Potential production rates of A. tonsa at summer temperatures were high ( \(\bar x\) ; St 1=19.0 mg C·m^?3 day^?1; \(\bar x\) ; St 5=22.9 mg C·m^?3 day^?1) and biomass doubling times were generally less than one day.     

A multi-isotope approach for understanding sources of water,carbon and sulfur in natural springs of the Central Appalachian region

Natural springs have been reliable sources of domestic water and have allowed for the development of recreational facilities and resorts in the Central Appalachians. The structural history of this area is complex and it is unknown whether these natural springs receive significant recharge from modern precipitation or whether they discharge old water recharged over geological times scales. The main objective of this study was to use stable isotopes of water ( $\delta^{18} {\text{O}}_{{{\text{H}}_{2} {\text{O}}}}$ and $\delta^{2} {\text{H}}_{{{\text{H}}_{2} {\text{O}}}}$ ), dissolved inorganic carbon ( $\delta^{13} {\text{C}}_{\text{DIC}}$ ) and dissolved sulfate ( $\delta^{34} {\text{S}}_{{{\text{SO}}_{4} }}$ and $\delta^{18} {\text{O}}_{{{\text{SO}}_{4} }}$ ) to delineate sources of water, carbon and sulfur in several natural springs of the region. Our preliminary isotope data indicate that all springs are being recharged by modern precipitation. The oxygen isotope composition indicates that waters in thermal springs did not encounter the high temperatures required for O isotope exchange between the water and silicate/carbonate minerals, and/or the residence time of water in the aquifers was short due to high flow rates. The carbon isotopic composition of dissolved inorganic carbon and sulfur/oxygen isotopic composition of dissolved sulfate provide evidence of low-temperature water–rock interactions and various biogeochemical transformations these waters have undergone along their flow path.     

Crystal chemistry of Fe3+-bearing (Mg,Fe)SiO3 perovskite: a single-crystal X-ray diffraction study

Magnesium silicate perovskite is the predominant phase in the Earth’s lower mantle, and it is well known that incorporation of iron has a strong effect on its crystal structure and physical properties. To constrain the crystal chemistry of (Mg, Fe)SiO₃ perovskite more accurately, we synthesized single crystals of Mg_0.946(17)Fe_0.056(12)Si_0.997(16)O₃ perovskite at 26 GPa and 2,073 K using a multianvil press and investigated its crystal structure, oxidation state and iron-site occupancy using single-crystal X-ray diffraction and energy-domain Synchrotron Mössbauer Source spectroscopy. Single-crystal refinements indicate that all iron (Fe²⁺ and Fe³⁺) substitutes on the A-site only, where \( {\text{Fe}}^{ 3+ } /\Upsigma {\text{Fe}}\sim 20\,\% \) based on Mössbauer spectroscopy. Charge balance likely occurs through a small number of cation vacancies on either the A- or the B-site. The octahedral tilt angle (Φ) calculated for our sample from the refined atomic coordinates is 20.3°, which is 2° higher than the value calculated from the unit-cell parameters (a = 4.7877 Å, b = 4.9480 Å, c = 6.915 Å) which assumes undistorted octahedra. A compilation of all available single-crystal data (atomic coordinates) for (Mg, Fe)(Si, Al)O₃ perovskite from the literature shows a smooth increase of Φ with composition that is independent of the nature of cation substitution (e.g., \( {\text{Mg}}^{ 2+ } - {\text{Fe}}^{ 2+ } \) or \( {\text{Mg}}^{ 2+ } {\text{Si}}^{ 4+ } - {\text{Fe}}^{ 3+ } {\text{Al}}^{ 3+ } \) substitution mechanism), contrary to previous observations based on unit-cell parameter calculations.     

Biradical states of oxygen-vacancy defects in ��-quartz: centers E_{2}^{\prime \prime } and E_{4}^{\prime \prime }

Several new radiation defects with total electron spin S?=?1 occurring in electron-irradiated, synthetic ??-quartz have been observed by using electron paramagnetic resonance spectroscopy. These defects are considered to be biradicals, i.e., pairs of S?=?1/2 species. The concentration of these centers depends on the condition of the fast-electron irradiation. They have different decay behaviors that allow measurements of any individual species especially when it predominates over the others. The primary spin Hamiltonian parameter matrices g ₁, g ₂, D have now been determined for two similar defects, which herein are labeled $ E_{2}^{\prime \prime } $ and $ E_{4}^{\prime \prime } $ . Inter-electron distances estimated by using the magnetic dipole model, suggest that the structures of centers $ E_{2}^{\prime \prime } $ and $ E_{4}^{\prime \prime } $ both involve the unpaired electrons each located in orbitals of two silicon atoms next to a common oxygen vacancy but which have slightly different Si?CSi distances at 0.90 and 0.79?nm, respectively. This model is consistent with previous DFT calculations of the triplet configurations with local energetic minima. Observed decay behaviors suggest a transformation of centers $ E_{2,4}^{\prime \prime } $ to the analogous $ E_{1}^{\prime \prime } $ center. These triplet centers in quartz provide new insights into the structures of analogous defects in amorphous silica.     

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.     

Flow rule in a high-cycle accumulation model backed by cyclic test data of 22 sands

The flow rule used in the high-cycle accumulation (HCA) model proposed by Niemunis et al. (Comput Geotech 32: 245, 2005) is examined on the basis of the data from approximately 350 drained long-term cyclic triaxial tests (N = 10⁵ cycles) performed on 22 different grain-size distribution curves of a clean quartz sand. In accordance with (Wichtmann et al. in Acta Geotechnica 1: 59, 2006), for all tested materials, the “high-cyclic flow rule (HCFR)”, i.e., the ratio of the volumetric and deviatoric strain accumulation rates \(\dot{\varepsilon}_{\rm{v}}^{{\rm acc}}/\dot{\varepsilon}_{\rm{q}}^{{\rm acc}}\) , was found dependent primarily on the average stress ratio η ^av = q ^av/p ^av and independent of amplitude, soil density and average mean pressure. The experimental HCFR can be fairly well approximated by the flow rule of the modified Cam-clay (MCC) model. Instead of the critical friction angle \(\varphi_{\rm{c}}\) which enters the flow rule for monotonic loading, the HCA model uses the MCC flow rule expression with a slightly different parameter \(\varphi_{\rm{cc}}\) . It should be determined from cyclic tests. \(\varphi_{\rm{cc}}\) and \(\varphi_{\rm{c}}\) are of similar magnitude but not always identical, because they are calibrated from different types of tests. For a simplified calibration in the absence of cyclic test data, \(\varphi_{\rm{cc}}\) may be estimated from the angle of repose \(\varphi_{\rm{r}}\) determined from a pluviated cone of sand (Wichtmann et al. in Acta Geotechnica 1: 59, 2006). However, the paper demonstrates that the MCC flow rule with \(\varphi_{\rm{r}}\) does not fit well the experimentally observed HCFR in the case of coarse or well-graded sands. For an improved simplified calibration procedure, correlations between \(\varphi_{\rm{cc}}\) and parameters of the grain-size distribution curve (d ₅₀,   C _u) have been developed on the basis of the present data set. The approximation of the experimental HCFR by the generalized flow rule equations proposed in (Wichtmann et al. in J Geotech Geoenviron Eng ASCE 136: 728, 2010), considering anisotropy, is also discussed in the paper.     

A dynamical model for the I\bar 1 - P\bar 1 phase transition in anorthite,CaAl2Si2O8

The non-ferroic triclinic to triclinic \(I\bar 1 - P\bar 1\) phase transition in anorthite is described in terms of the spontaneous onset of an order parameter η. A triclinic to triclinic phase transition can be driven by order parameters (representations) arising from the Γ, Z, X, U, V, R, Y, and T points of symmetry of the Brillouin zone. Each point leads to a set of two inequivalent representations and thus there is a total of sixteen inequivalent order parameters. However, only the R ₁ ⁺ representation is consistent with the change from the body-centered to primitive cell (increase of primitive cell size of two) and also with the origin of the two space groups (inversion center) being at the same position. The R ₁ ⁺ order parameter of the high symmetry triclinic phase \(P\bar 1_0\) (or equivalently \(I\bar 1\) ) causes a reciprocal lattice change and, in terms of the lower symmetry reciprocal lattice, the order parameter corresponds to the b^* point. This is consistent with experimentally observed x-ray diffuse scattering. Using induced representation theory, microscopic distortions compatible with the R ₁ ⁺ order parameter are obtained. Assuming a distortion in an arbitrary direction at the general 2(i) Wyckoff position (x₀,y₀,z₀) of \(P\bar 1_0\) (the higher symmetry phase) induced representation theory demands an opposite displacement at the position (x₀, y₀, z₀), an opposite displacement at (x₀+1,y₀+1,z₀+1), and the same displacement at ( \(\bar x\) ₀+1, \(\bar y\) ₀+1, \(\bar z\) ₀+1) of \(P\bar 1_0\) . This is also consistent with experiment. The presence of the weak c-type reflections above the transition is attributed to the fluctuating lower symmetry antiphase domains related by the translation (1/2, 1/2, 1/2).     

Transmission electron microscopy of experimentally deformed chalcopyrite single crystals

Two crystals of natural chalcopyrite, CuFeS₂, experimentally deformed at 200° C have been studied by means of transmission electron microscopy (TEM). The activated glide planes are (001) and {112}. The dislocations in (001) have the Burgers vector [110] and a predominating edge character. They are split into two colinear partials b=1/2[110] and can cross split into {112}. The dislocations in {112} consist of straight segments along low index lattice lines. They are often arranged in dipoles generating trails of loops. Few dislocations with b=1/2[ \(\overline {11} \) 1] and [1 \(\bar 1\) 0] are present and dislocations with b=[0 \(\bar 2\) 1] occur in low angle subgrain boundaries. From weak beam contrasts it is presumed that most of the dislocations gliding in {112} have b=1/2〈3 \(\overline {11} \) 〉. They are dissociated into up to four partials. Microtwins and different types of stacking faults in {112} also occur. Models of the dissociation of dislocations are discussed.