Numerical evaluation of soft inter‐slab joint in concrete‐faced rockfill dam with dual mortar finite element method

Recently constructed concrete‐faced rockfill dams (CFRDs) often use soft inter‐slab joints to prevent axial compression‐induced extrusion damage in the concrete face. Due to the complexity of the multibody contact and the lack of information on the actual behavior of soft joints, it is highly challenging to numerically assess the effect of soft joints in CFRDs. In this paper, we present a numerical approach for the three‐dimensional modeling of CFRDs with hard and soft joints. A dual mortar finite element method with Lagrange multiplier is developed to treat the multibody contact in hard joints with impenetrability condition. The soft joint slab‐filler‐slab contact system is modeled using an equivalent contact interface approach, where the soft contact constraints are imposed using a perturbed Lagrange formulation. Through a series of laboratory tests, the mechanical behavior of soft joint is investigated. An extrusion model for the soft joint is presented and implemented in the dual mortar finite element method. The proposed numerical method is applied to the three‐dimensional analysis of Tianshengqiao‐1 CFRD. Despite the complex multibody contact and strong material and geometry nonlinearities in the CFRD, the proposed method is stable and capable of capturing salient characteristics of the CFRD. Numerical results show that in Tianshengqiao‐1, the employment of soft joints can effectively reduce the axial compression stress, thus greatly alleviating the risk of extrusion damage in the concrete face.     

Stability analysis of shotcrete supported crown of NATM tunnels with discontinuity layout optimization

For tunnel constructed by New Austrian Tunnelling Method, the crown is the upper part of tunnel section, constructed during excavation process and supported by shotcrete. The stability of the crown has great influence on the safety of tunnel itself and the buildings above, which correlates, among others, with geometrical setup of tunnel and material properties of shotcrete and soil/rock. In this paper, aiming at analyzing the stability of shotcrete supported crown, a recently presented numerical method discontinuity layout optimization is adopted, which introduces a great amount of potential discontinuities cross over one another and provides a wide search space for efficient upper limit analysis. In the analysis, a well‐established hydration model of cementitious material is implemented for accounting the hydration of shotcrete. Then assumptions based on convergence‐confinement method are used for accounting the 3‐dimensional effect in 2‐dimensional analysis, finally providing time‐space–dependent assessments of stability of shotcrete supported crown.     

Numerical simulation of agricultural sediment and pesticide runoff: RZWQM and PRZM comparison

Agricultural sediment and pesticide runoff is a widespread ecological and human health concern. Numerical simulation models, such as Root Zone Water Quality Model (RZWQM) and Pesticide Root Zone Model (PRZM), have been increasingly used to quantify off‐site agricultural pollutant movement. However, RZWQM has been criticized for its inability to simulate sedimentation processes. The recent incorporation of the sedimentation module of Groundwater Loading Effects of Agricultural Management Systems has enabled RZWQM to simulate sediment and sediment‐associated pesticides. This study compares the sediment and pesticide transport simulation performance of the newly released RZWQM and PRZM using runoff data from 2 alfalfa fields in Davis, California. A composite metric (based on coefficient of determination, Nash–Sutcliffe efficiency, index of agreement, and percent bias) was developed and employed to ensure robust, comprehensive assessment of model performance. Results showed that surface water runoff was predicted reasonably well (absolute percent bias <31%) by RZWQM and PRZM after adjusting important hydrologic parameters. Even after calibration, underestimation bias (?89% ≤ PBIAS  ≤ ?36%) for sediment yield was observed in both models. This might be attributed to PRZM's incorrect distribution of input water and uncertainty in RZWQM's runoff erosivity coefficient. Moreover, the underestimation of sediment might be less if the origin of measured sediment was considered. Chlorpyrifos losses were simulated with reasonable accuracy especially for Field A (absolute PBIAS  ≤ 22%), whereas diuron losses were underestimated to a great extent (?98% ≤ PBIAS  ≤ ?65%) in both models. This could be attributed to the underprediction of herbicide concentration in the top soil due to the limitations of the instantaneous equilibrium sorption model as well as the high runoff potential of herbicide formulated as water‐dispersible granules. RZWQM and PRZM partitioned pesticides into the water and sediment phases similarly. According to model predictions, the majority of pesticide loads were carried via the water phase. On the basis of this study, both RZWQM and PRZM performed well in predicting runoff that carried highly adsorptive pesticides on an event basis, although the more physically based RZWQM is recommended when field‐measured soil hydraulic properties are available.     

Using stable isotopes to understand seasonal and interannual dynamics in moisture sources and atmospheric circulation in precipitation

Stable isotopes in precipitation are useful tracers to strengthen understanding of climate change and hydrological processes. In this study, the moisture sources of 190 precipitation events in Beijing were analysed using the Hybrid Single‐particle Lagrangian Integrated Trajectory model, based on which we studied the relation between variations in precipitation δ¹⁸O and dynamics in moisture sources and atmospheric circulation in seasonal and interannual timescales. Categorization of 7 groups of moisture sources was performed, among which oceanic moisture sources presented lower δ¹⁸O in precipitation than continental moisture sources. The results show that seasonal variations of precipitation δ¹⁸O were caused by changes of moisture sources. In summer, moisture from proximal oceans dominated vapour transport to Beijing due to increasing monsoon strength and resulted in a relatively small variation of precipitation δ¹⁸O. At the interannual timescale, the variations of δ¹⁸O in summer precipitation were related to dynamics in oceanic moistures, showing depleted values when the contribution of oceanic moistures, especially the proportion of long‐distance oceanic moisture, was high. Further analysis indicated that changes of oceanic moisture sources were controlled by the strength of summer monsoons. These findings address the complexity of moisture sources in midlatitude monsoon areas and suggest that isotopic signals in precipitation have the potential to deduce changes in moisture sources and atmospheric circulation and can therefore serve for palaeoclimate reconstruction.     

Variability in the vertical hyporheic water exchange affected by hydraulic conductivity and river morphology at a natural confluent meander bend

River confluences and their associated tributaries are key morphodynamic nodes that play important roles in controlling hydraulic geometry and hyporheic water exchange in fluvial networks. However, the existing knowledge regarding hyporheic water exchange associated with river confluence morphology is relatively scarce. On January 14 and 15, 2016, the general hydraulic and morphological characteristics of the confluent meander bend (CMB) between the Juehe River and the Haohe River in the southern region of Xi'an City, Shaanxi Province, China, were investigated. The patterns and magnitudes of vertical hyporheic water exchange (VHWE) were estimated based on a one‐dimensional heat steady‐state model, whereas the sediment vertical hydraulic conductivity (K_v) was calculated via in situ permeameter tests. The results demonstrated that 6 hydrodynamic zones and their extensions were observed at the CMB during the test period. These zones were likely controlled by the obtuse junction angle and low momentum flux ratio, influencing the sediment grain size distribution of the CMB. The VHWE patterns at the test site during the test period mostly showed upwelling flow dominated by regional groundwater discharging into the river. The occurrence of longitudinal downwelling and upwelling patterns along the meander bend at the CMB was likely subjected to the comprehensive influences of the local sinuosity of the meander bend and regional groundwater discharge and finally formed regional and local flow paths. Additionally, in dominated upwelling areas, the change in VHWE magnitudes was nearly consistent with that in K_v values, and higher values of both variables generally occurred in erosional zones near the thalweg paths of the CMB, which were mostly made up of sand and gravel. This was potentially caused by the erosional and depositional processes subjected to confluence morphology. Furthermore, lower K_v values observed in downwelling areas at the CMB were attributed to sediment clogging caused by local downwelling flow. The confluence morphology and sediment K_v are thus likely the driving factors that cause local variations in the VHWE of fluvial systems.     

Petrological evidence for shock‐induced high‐P metamorphism in a gabbro

Microlites (minute spherulitic, dendritic, skeletal, acicular and poikilitic crystals) diagnostic of crystallization in quenched melt or glass in fault rocks have been used to infer fossil earthquakes. High‐P microlites and crystallites are described here in a variably eclogitized gabbro, the wallrock to the coesite‐bearing eclogite breccia at Yangkou in the Chinese Su‐Lu high‐P metamorphic belt. The studied hand specimens are free of discernible shear deformation, although microfractures are not uncommon under the microscope. In the least eclogitized gabbro, the metagabbro, stellate growths of high‐P minerals on the relict igneous minerals are common. Dendritic garnet crystals (<1?5 μm) grew around rutile and/or phengite replacing ilmenite and biotite, respectively. Skeletal garnet also rims broken flakes of igneous biotite and mechanically twinned augite. Radial intergrowths of omphacite and quartz developed around relict igneous orthopyroxene and are rimmed by skeletal or poikilitic garnet where a Ti‐bearing mineral relict is present. Acicular epidote, kyanite and phengite crystallites are randomly distributed in a matrix of Na‐rich plagioclase, forming the pseudomorphs after igneous plagioclase. In the more eclogitized gabbro, the coronitic eclogite located closer to the eclogite breccia, all the igneous minerals broke down into high‐P assemblages. Thick coronas of poikilitic garnet grew between the pseudomorphs after igneous plagioclase and ferromagnesian minerals. The igneous plagioclase is replaced by omphacite crystallites, with minor amounts of phengite and kyanite. Thermodynamic modelling of the plagioclase pseudomorphs shows an increase in P–T in the wallrock from the metagabbro to the coronitic eclogite, and the P–T variation is unrelated to H₂O content. The fluid‐poor pressure overstepping scenario is unsupported both by phase diagram modelling and by whole‐rock chemical data, which show that the various types of eclogitized gabbro are all fairly dry. A large pressure difference of >2 GPa between the metagabbro and the coesite‐bearing eclogites ~20 m apart cannot be explained by the subduction hypothesis because this would require a depth difference of >60 km. The microlites and crystallites are evidence for dynamic crystallization due to rapid cooling because constitutional supercooling was unlikely for the plagioclase pseudomorphs. The lack of annealing of the broken biotite and augite overgrown by strain free skeletal garnet is consistent with a transient high‐P–T event at a low ambient temperature (<300 °C), probably in the crust. Therefore, the eclogitization of the wallrock to the eclogite breccia was also coseismic, as proposed earlier for the eclogite facies fault rocks. The outcrop‐scale P–T variation and the transient nature of the high‐P–T event are inconsistent with the other existing tectonic models for high‐P metamorphism. The fact that the less refractory but denser biotite is largely preserved while the more refractory but less dense plagioclase broke down completely into high‐P microlite assemblages in the metagabbro indicates a significant rise in pressure rather than temperature. Given that the metamorphic temperatures are far below the melting temperatures of most of the gabbroic minerals under fluid‐absent conditions, stress‐induced amorphization appears to be the more likely mechanism of the coseismic high‐P metamorphism.     

The Ore-forming Fluid of the Gold Deposits of Muru Gold Belt in Eastern Shandong, China - a Case Study of Denggezhuang Gold Deposit

Abstract. Denggezhuang gold deposit is an epithermal gold‐quartz vein deposit in northern Muru gold belt, eastern Shandong, China. The deposit occurs in the NNE‐striking faults within the Mesozoic granite. The deposit consists of four major veins with a general NNE‐strike. Based on crosscutting relationships and mineral parageneses, the veins appear to have been formed during the same mineralization epochs, and are further divided into three stages: (1) massive barren quartz veins; (2) quartz‐sulfides veins; (3) late, pure quartz or calcite veinlets. Most gold mineralization is associated with the second stage. The early stage is characterized by quartz, and small amounts of ore minerals (pyrite), the second stage is characterized by large amounts of ore minerals. Fluid inclusions in vein quartz contain C‐H‐O fluids of variable compositions. Three main types of fluid inclusions are recognized at room temperature: type I, two‐phase, aqueous vapor and an aqueous liquid phase (L+V); type II, aqueous‐carbonic inclusions, a CC₂‐liquid with/without vapor and aqueous liquid (LCO₂+VCC₂+Laq.); type III, mono‐phase aqueous liquid (Laq.). Data from fluid inclusion distribution, microthermometry, and gas analysis indicate that fluids associated with Au mineralized quartz veins (stage 2) have moderate salinity ranging from 1.91 to 16.43 wt% NaCl equivalent (modeled salinity around 8–10 wt% NaCl equiv.). These veins formatted at temperatures from 80d? to 280d?C. Fluids associated with barren quartz veins (stage 3) have a low salinity of about 1.91 to 2.57 wt% NaCl equivalent and lower temperature. There is evidence of fluid immiscibility and boiling in ore‐forming stages. Stable isotope analyses of quartz indicate that the veins were deposited by waters with δ^O and δD values ranging from those of magmatic water to typical meteoric water. The gold metallogenesis of Muru gold belt has no relationship with the granite, and formed during the late stage of the crust thinning of North China.     

Effect of some additives on synthesis of zeolite from coal fly ash

Hydrothermal conversion of fly ash into zeolites was conducted and the effects of the addition of sodium halide and waste solutions produced after zeolitization of fly ash, as well as the adjustment of the Si/Al ratio prior to synthesis process on the formation and cation exchange capacity (CEC) of zeolite product were evaluated. Both the addition of NaCl and NaF ameliorated the crystallinity and CEC of synthesized zeolite, but NaF had a better improvement effect. Na⁺ was considered to enhance the crystallization of zeolite, while F⁻ favored the dissolution of fly ash. The type of zeolite formed depended on the Si/Al ratio of the starting material prior to the nucleation and crystallization of zeolite. The adjustment of the Si/Al ratio of fly ash by addition of Na₂SiO₄ and Al(OH)₃ changed the type and CEC of zeolite. Waste solutions contained large amount of Si and little Al due to the formation of a zeolite named NaP1 in zeolite terminology with the Joint Committee of Powder Diffraction Standard (JCPDS) code of 39-0219. The alkalinity decreased largely. As a result, the CEC value of zeolite products synthesized with waste solution as alkali source decreased. The supplementation of new alkali to adjust the alkalinity of waste solution could enhance the CEC of synthesized product. It was concluded that: (1) addition of sodium halide and adjustment of the Si/Al ratio prior to synthesis can improve the quality of zeolite; (2) waste solutions produced following the zeolitization of fly ash can be reused as an alkali source in the activation of fly ash; zero-emission of waste solution in the synthesis of zeolite from fly ash is possible.     

Constraining the hydrocarbon expulsion history of the coals in Qinshui Basin, North China, from the analysis of fluid inclusions

From the comprehensive study on the homogenization temperatures and the occurrence of fluid inclusions in the framework minerals of the strata between or above the Carboniferous–Permian coals in the Qinshui basin, Shanxi, three stages are predicted of hydrocarbon expulsion from the coals. Combined with the known history of basin evolution, it is deduced that the expulsion of hydrocarbons happened during the J₁ (210–180 Ma), the early K₁ (150–130 Ma) and K₂E₁ (110–60 Ma). In the early stage, the coals produced and discharged coal-generated oils. The average GOI value of four sandstone samples is relatively high, as they have been exposed to high paleo-oil saturation in the strata between or above the coals. The biomarker compositions of oil-bearing fluid inclusions are similar to those of extracts from the coals, and so it is concluded that those oils were derived from the same family of the coals.     

http://www.sciencedirect.com/science/article/pii/S1674987114001121

The Chinese Tianshan belt of the southern Altaids has undergone a complicated geological evolution.Different theories have been proposed to explain its evolution and these are still hotly debated.The major subduction polarity and the way of accretion are the main problems.Southward,northward subduction and multiple subduction models have been proposed.This study focuses on the structural geology of two of the main faults in the region,the South Tianshan Fault and the Nikolaev Line.The dip direction in the Muzhaerte valley is southward and lineations all point towards the NW.Two shear sense motions have been observed within both of these fault zones,a sinistral one,and a dextral one,the latter with an age of 236-251 Ma.Structural analyses on the fault zones show that subduction has been northward rather than southward.The two shear sense directions indicate that the Yili block was first dragged along towards the east due to the clockwise rotation of the Tarim block.After the Tarim block stopped rotating,the Yili block still kept going eastward,inducing the dextral shear senses within the fault zones.