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.     

Performances of two instrumented laboratory models for the study of rainfall infiltration into unsaturated soils

Rainfall infiltration and suction variation in unsaturated soils must be taken into consideration in the analysis of most slope stability problems, particularly in the tropical regions where the annual precipitation is high. The process of rainfall infiltration into unsaturated soils is an extremely complex problem attributed to the non-linearity of the hydraulic property functions of the unsaturated soils. This paper describes in detail two instrumented laboratory models, i.e. one-dimensional soil column, and two-dimensional slope model used to provide experimental evidences for the transient suction variations in the unsaturated soils under certain rainfall conditions. The performances of the laboratory models were tested on four typical types of residual soils, i.e. sand-gravel, silty gravel, sandy silt, and silt (kaolin), and a two-layered soil system, i.e. sandy silt underlain by silty gravel. The results showed that the suction distributions for the single-layered homogeneous soils obtained from the simpler one-dimensional soil column were almost identical to that of two-dimensional slope model. However, the two-dimensional slope model should be employed for the two-layered soil system because of the dominant effect of the lateral flow mechanism. The capillary barrier effect was observed when a less permeable soil layer was underlain by a more permeable soil layer. The minimum suction value in soil is governed by the rainfall intensity, rainfall duration, and the saturated permeability of soil. The infiltration rate of the fine-grained soils that subjected to shrink and crack was independent of the soil permeability, but was significantly governed by the preferential flows developed in the soils.     

Effect of stress history on CPT and DMT results in sand

In order to investigate the effect of stress history on in-situ test results in granular sediments, a series of CPTs and DMTs are performed on Busan sand prepared in the calibration chamber. K_D is found to be the most sensitive to the stress history among CPT and DMT measurements. E_D and q_c are observed to be similarly affected by the stress history and, therefore, the E_D–q_c relation appears to be almost independent of the stress history. The K_D–D_R relation established without considering the stress history is likely to overestimate the relative density of OC sand. It is shown that the existence of the pre-stress of the granular sediment can be indirectly recognized by an estimation of the relative density larger than 100% when using the K_D–σ_v′–D_R relation suggested for NC sand. Although q_c/σ_v′–K_D/K₀ and E_D/σ_v′–K_D/K₀ relations are heavily influenced by the stress history, q_c/σ_m′–K_D/K₀ and E_D/σ_m′–K_D/K₀ relations are observed to be independent of the stress history. Based on these relations, charts to evaluate the K₀ value from q_c and/or DMT indices are developed for both NC and OC sands. The design chart based on E_D/σ_m′–K_D/K₀ and E_D/σ_v′–K_D/K₀ relations is expected to be practically useful as the usage of this chart requires only DMT indices. The developed design charts are applicable to Busan sand but different sets of equations and charts may be developed for other sands.     

Ar loss in experimentally deformed muscovite and biotite with implications for Ar/Ar geochronology of naturally deformed rocks

The effects of deformation on radiogenic argon (⁴⁰Ar^∗) retentivity in mica are described from high pressure experiments performed on rock samples of peraluminous granite containing euhedral muscovite and biotite. Cylindrical cores, ∼15 mm in length and 6.25 mm in diameter, were drilled from granite collected from the South Armorican Massif in northwestern France, loaded into gold capsules, and weld-sealed in the presence of excess water. The samples were deformed at a pressure of 10 kb and a temperature of 600 °C over a period 29 of hours within a solid medium assembly in a Griggs-type triaxial hydraulic deformation apparatus. Overall shortening in the experiments was approximately 10%. Transmitted light and secondary and backscattered electron imaging of the deformed granite samples reveals evidence of induced defects and for significant physical grain size reduction by kinking, cracking, and grain segmentation of the micas.Infrared (IR) laser (CO₂) heating of individual 1.5-2.5 mm diameter grains of muscovite and biotite separated from the undeformed granite yield well-defined ⁴⁰Ar/³⁹Ar plateau ages of 311 ± 2 Ma (2σ). Identical experiments on single grains separated from the experimentally deformed granite yield results indicating ⁴⁰Ar^∗ loss of 0-35% in muscovite and 2-3% ⁴⁰Ar^∗ loss in biotite. Intragrain in situ ultraviolet (UV) laser ablation ⁴⁰Ar/³⁹Ar ages (±4-10%, 1σ) of deformed muscovites range from 309 ± 13 to 264 ± 7 Ma, consistent with 0-16% ⁴⁰Ar^∗ loss relative to the undeformed muscovite. The in situ UV laser ablation ⁴⁰Ar/³⁹Ar ages of deformed biotite vary from 301 to 217 Ma, consistent with up to 32% ⁴⁰Ar^∗ loss. No spatial correlation is observed between in situ⁴⁰Ar/³⁹Ar age and position within individual grains. Using available argon diffusion data for muscovite the observed ⁴⁰Ar^∗ loss in the experimentally treated muscovite can be utilized to predict average ⁴⁰Ar^∗ diffusion dimensions. Maximum ⁴⁰Ar/³⁹Ar ages obtained by UV laser ablation overlap those of the undeformed muscovite, indicating argon loss of <1% and an average effective grain radius for ⁴⁰Ar^∗ diffusion ?700 μm. The UV laser ablation and IR laser incremental ⁴⁰Ar/³⁹Ar ages indicating ⁴⁰Ar^∗ loss of 16% and 35%, respectively, are consistent with an average diffusion radius ?100 μm. These results support a hypothesis of grain-scale ⁴⁰Ar^∗ diffusion distances in undeformed mica and a heterogeneous mechanical reduction in the intragrain effective diffusion length scale for ⁴⁰Ar^∗ in deformed mica. Reduction in the effective diffusion length scale in naturally deformed samples occurs most probably through production of mesoscopic and submicroscopic defects such as, e.g., stacking faults. A network of interconnected defects, continuously forming and annealing during dynamic deformation likely plays an important role in controlling both ⁴⁰Ar^∗ retention and intragrain distribution in deformed mica. Intragrain ⁴⁰Ar/³⁹Ar ages, when combined with estimates of diffusion kinetics and distances, may provide a means of establishing thermochronological histories from individual micas.     

An extensive erosion surface of a strongly deformed limestone bed in the Gushan and Chaomidian formations (late Middle Cambrian to Furongian), Shandong Province,China: Sequence–stratigraphic implications

In order to understand sequence development and sea-level fluctuations during the late Middle Cambrian to early Furongian on the North China epeiric platform, the present study focuses on a unique, subtle erosion surface of an extensive (approx. 100 km), strongly deformed limestone bed in the uppermost part of the Gushan Formation, Shandong Province, China. The Gushan Formation and the overlying Chaomidian Formation consist mainly of shales and a variety of carbonates that were deposited in subtidal environments (e.g., deep subtidal, shallow subtidal, shoreface/shoal, subtidal microbial flat, and restricted platform interior). Three third-order depositional sequences (S1–3) are identified, each of which comprises a thin transgressive systems tract (TST) and a relatively thick highstand systems tract (HST). Each sequence is bounded by a drowning unconformity (SB1), a subaerial unconformity (SB2), or a surface of submarine erosion (SB3). The upper sequence boundary (SB2) of sequence 1 (S1) is represented by a subtle erosion surface of an extensive, deformed limestone bed with a wide variety of soft-sediment deformation structures (e.g., lime mudstone breccias, chaotic wacke-packstone laminae and fragments, homogenized oolites, and clastic dykes), and is overlain by small sporadic microbial buildups and an extensive bioclastic grainstone bed. The deformed limestone was formed during early diagenesis by differential deformation processes (brecciation, liquefaction/fluidization, and injection) which were most likely induced by pore-water overpressure during the period of rapid sea-level fall. Despite the lack of subaerial exposure features (e.g., paleokarst, paleosol, etc.), the characteristics of the erosion surface (cutting well-lithified sediment below), the missing of a significant geological record (the Prochuangia biozone), and the worldwide correlatable positive carbon isotope excursion collectively indicate that the erosion surface developed under conditions of subaerial exposure after contemporaneous marine cementation of the deformed sediment. The missing of the Prochuangia biozone is most likely due to non-deposition at a subaerial hiatal surface. The erosion surface was submerged as a result of subsequent rise in sea level, where sporadic microbial buildups formed under suitable conditions. Freshly deposited, winnowed, shell-dominated transgressive lag deposits (containing Chuangia trilobite fragments, brachiopod shells, and abundant glauconite grains) formed with continued rise in sea level, which became, in turn, overlain by shale-dominated facies. The unique combination of the subtle erosion surface (sensu stricto a subaerial unconformity) and the underlying deformed limestone bed provides an important criterion for recognizing the subtle changes in relative sea level on shallow epeiric platforms.     

Estimation of soil erosion rates in a subtropical mountain watershed using 137Cs radionuclide

The formulation of watershed management strategies to protect water resources threatened by soil erosion and sedimentation requires a thorough understanding of sediment sources and factors that drive soil movement in the watershed. This paper describes a study of medium-term water-driven soil erosion rates in a mountainous watershed of the Shihmen Reservoir in Taiwan. A total of 60 sampling sites were selected along a hillslope. At each sampling site, the inventory ¹³⁷Cs activity was determined and then calculated with the diffusion and migration model to derive soil erosion rates. The rates are one to two orders of magnitude lower than estimates using the Universal Soil Loss Equation, a soil erosion model often used in Taiwan. Results of multiple regression analysis indicate that the spatial variability of soil erosion rates is associated with the relative position of a sampling site to the nearest ridge and soil bulk densities (r ² = 0.33, p < 0.01). Finally, the patterns of soil redistribution rates on the hillslope follow the ¹³⁷Cs hillslope model as soil erosion increases in the downslope direction. No deposition site is found at footslope because soil deposition is swept away by regular flooding along the stream channel. This study is an important first step in using ¹³⁷Cs as a tracer of soil redistribution in mountainous watersheds of Taiwan.     

Test method of cavitation erosion for marine coatings with low hardness

Rudders of large container ships are easily affected by cavitation, which is well known to be induced by significant axial flows behind a propeller and discontinuities in the rudder. Among several methods to prevent or reduce the cavitation erosion occurred in the rudder, painting is gaining a lot of attention because it can be employed easily and cheaply. To conduct erosion tests properly, the simulation of heavily erosive cavitation is necessary. This can be generated using an inclined propeller dynamometer in the medium-size cavitation tunnel of MOERI (Maritime & Ocean Engineering Research Institute). The inclined shaft of the propeller creates strong cavitation, which occurs around the root of the propeller blade. This cavitation creates impacts through the collapsing process that are very severe, and are useful for realistic and efficient cavitation erosion tests. In the present study, the newly developed cavitation erosion test method is successfully employed to evaluate marine coatings that is mainly composed of epoxy elastomer or silicone polymer material. Silicone polymer-type paint B was found to have three times larger endurance than epoxy elastomer-type paint A.     

Modeling and control of a 75 kW class variable liquid-column oscillator for highly efficient wave energy converter

The modeling and control of a variable liquid-column oscillator having a liquid filled U-tube with air chambers at its vertical columns are presented. As an ocean wave energy extracting device, the structure of the variable liquid-column oscillator (VLCO) is analogous to that of the tuned liquid-column damper used to suppress oscillatory motion in large structures like tall buildings and cargo ships. However, owing to an air spring effect caused by the dynamic pressure of air chambers, the amplitude of response of the VLCO becomes significantly amplified for a desired wave period. The governing equations for the motion of VLCO structure under wave excitation and the motion of liquid with an air spring effect caused by an air–liquid interaction are described by a series of nonlinear differential equations. A set of control parameters for extracting maximum power from various wave conditions is determined for the efficient operation of the VLCO. It is found that the effect of the air spring has an important role to play in making the oscillation of the VLCO match with the ocean wave. In this way, the VLCO provides the most effective mode for extracting energy from the ocean wave.     

Dynamic factor analysis considering elastic boom effects in heavy lifting operations

The dynamic factor is the ratio of the maximum dynamic load to the static load acting on the wire ropes between the boom of a floating crane and a cargo. In this paper, the dynamic factor is analyzed based on dynamic simulations of a floating crane and a cargo, considering an elastic boom. For the simulation, we designed a multibody system that consists of a floating crane barge, an elastic boom, and a cargo connected to the boom through wire ropes. The dynamic equations of motion of the system are based on flexible multibody system dynamics. Six-degree-of-freedom motions are considered for the floating crane and for the cargo, and three-dimensional deformations for the elastic boom. The hydrostatic force, the hydrodynamic force, the gravitational force, and the wire rope forces are considered as external forces. The dynamic factor is obtained by numerically solving the equation. The effects of the elastic boom on heavy cargo lifting are discussed by comparing the simulation results of an elastic boom and a rigid boom.     

Two episodes of the Indosinian thermal event on the South China Block: Constraints from LA-ICPMS U–Pb zircon and electron microprobe monazite ages of the Darongshan S-type granitic suite

The Darongshan granitic suite (~ 10,000 km²) consists of five major units (Taima, Nadong and Jiuzhou plutons, and Pubei and Darongshan batholiths) typical of peraluminous S-type granitoids containing abundant granulite inclusions in the Cathaysia block, South China. Six samples from these plutons and batholiths have been investigated using both LA-ICPMS U–Pb age dating on zircon cores and EMP U–Th–Pb chemical age dating on monazite cores and rims. LA-ICPMS zircon results give similar major age populations ranging between 260 ± 3 and 250 ± 3 Ma for all units, with apparent older age peaks concentrated at 1020, 800, 430 and 330 Ma. On the other hand, EMP monazite results yield younger ages of 231–229 Ma for Nadong, Taima, Pubei and Darongshan and 224 Ma for Jiuzhou samples, with older age groups of 264 Ma for Taima and 256–250 Ma for Pubei units. Since the older monazite ages are similar to the majority of zircon ages, the latter are considered as inherited ages. Further because such zircon ages are similar with the emplacement time of the Emeishan large igneous province in western South China, they likely reflect the timing of metamorphism for the included fragments of granulitic crusts that had been formed by invasion of the Emeishan plume. The younger monazite ages, as present for all plutons and batholiths in the entire Darongshan area, are taken as the formation age of the host granites. Combining U–Pb zircon and EMP monazite ages known for Permo-Triassic high temperature and high pressure metamorphic rocks and granites in the Indochina block (e.g., the Kannack Complex of the Kontum massif), it is suggested that the Indosinian thermal activity had set records over both the Indochina (plus Simao) and South China blocks in two main episodes, one is 260–250 Ma and the other is 231–229 Ma. One plausible explanation is that these two blocks were one united continent before the Emeishan plume activity and an opening was triggered by this plume at ~ 260 Ma. Due to forces of the approaching Sibumasu block, both the South China and Indochina blocks were amalgamated again at ~ 230 Ma. We, therefore, advocate that double subduction of the plume-triggered oceanic crusts in opposite directions is responsible for the generation of the Darongshan granitic suite in the South China block and its counterpart in the Indochina block.