Numerical analysis on seepage failures of dike due to water level-up and rainfall using a water–soil-coupled smoothed particle hydrodynamics model

For seepage failures of dike due to water level-up and rainfall, surface infiltration and strength change induced by suction reduction are important factors; thus, numerical analysis should consider the coupling of water and soil, as well as the effect of saturation to obtain more precise failure mechanism. Based on the advanced smoothed particle hydrodynamics (SPH) method, this work proposed a two-phase-coupled SPH model in coordination with a novel constitutive model for unsaturated soils. Then, a triaxial compression test is simulated to check the applicability of the SPH method on the soil phase. After that, the failure test of a dike due to water level-up is discretized and simulated, from which the seepage process, the distribution of maximum shear strain, the slip surface, and pore water pressure are obtained. The two-phase-coupled SPH model is also applied to a slope failure test of heavy rainfall, and the results are compared to the model test. Finally, a dike failure test due to rainfall is analyzed using the proposed SPH model to reproduce the surface infiltration and suction reduction. The proposed SPH model provides several insights of seepage failures and can be a helpful tool for the analysis of dike failures induced by water level-up and rainfall.     

Groundwater resources management under environmental constraints in Shiroishi of Saga plain, Japan

A sinking of the land surface due to the pumping of groundwater has long been recognized as an environmental issue in the Shiroishi plain of Saga, Japan. Land subsidence can have several negative economic and social implications such as changes in groundwater and surface water flow patterns, restrictions on pumping in land subsidence prone areas, localized flooding, failure of well casings as well as shearing of structures. To minimize such an environmental effect, groundwater management should be considered in this area. In this study, a new integrated numerical model that integrates a three-dimensional numerical groundwater flow model coupled with a one-dimensional soil consolidation model and a groundwater optimization model was developed to simulate groundwater movement, to predict ground settlement and to search for optimal safe yield of groundwater without violating physical, environmental and socio-economic constraints. It is found that groundwater levels in the aquifers greatly vary from season to season in response to the varying climatic and pumping conditions. Consequently, land subsidence has occurred rapidly throughout the area with the Shiroishi plain being the most prone. The predicted optimal safe yield of the pumping amount is about 5 million m³. The study also suggests that pumping with this optimal amount will minimize the rate of land subsidence over the entire area. An erratum to this article can be found at     

A New Solar Imaging System for Observing High-Speed Eruptions: <Emphasis Type="Italic">Solar Dynamics Doppler Imager</Emphasis> (SDDI)

A new solar imaging system was installed at Hida Observatory to observe the dynamics of flares and filament eruptions. The system (Solar Dynamics Doppler Imager; SDDI) takes full-disk solar images with a field of view of \(2520~\mbox{arcsec} \times 2520~\mbox{arcsec}\) at multiple wavelengths around the \(\mathrm{H}\alpha\) line at 6562 Å. Regular operation was started in May 2016, in which images at 73 wavelength positions spanning from \(\mathrm{H}\alpha -9~\mathring{\mathrm{A}}\) to \(\mathrm{H}\alpha +9~\mathring{\mathrm{A}}\) are obtained every 15 seconds. The large dynamic range of the line-of-sight velocity measurements (\({\pm}\,400~\mbox{km}\,\mbox{s}^{-1}\)) allows us to determine the real motions of erupting filaments in 3D space. It is expected that SDDI provides unprecedented datasets to study the relation between the kinematics of filament eruptions and coronal mass ejections (CME), and to contribute to the real-time prediction of the occurrence of CMEs that cause a significant impact on the space environment of the Earth.     

A two-scale constitutive framework for modelling localised deformation in saturated dilative hardening material

Undrained deformation of dilative sand generates negative excess pore pressure. It enhances the strength, which is called dilative hardening. This increased suction is not permanent. The heterogeneity at the grain scale triggers localisations causing local volume changes. The negative hydraulic gradient drives fluid into dilating shear zones. It loosens the soil and diminishes the shear strength. It is essential to understand the mechanism behind this internal drainage and to capture it numerically. The purpose of this paper is to develop a macroscopic constitutive relationship for the undrained deformation of saturated dense sand in the presence of a locally fully or partially drained shear band. Separate constitutive relations are generated for the band and intact material. Both time and scale dependence during pore fluid diffusion in saturated sand are captured, eliminating the mesh dependency for finite element implementations. The model is applied to the Gauss points that satisfy the bifurcation criterion. The proposed method is calibrated to recreate the undrained macroscopic response bestowed by an extra-small mesh. The microscopic behaviours inside and outside shear band predicted by this model are qualitatively in good agreement with individual material point behaviours inside and outside the shear band in the extra-small mesh. Depending on the loading rate and the shear band thickness, the response inside the band can be fully or partially drained, which governs the ultimate global strength. The calibrated model is exploited to simulate an upscaled biaxial compression test with semipermeable boundaries.     

A Possible Role of Boiling in Ore Deposition: A Numerical Approach

Abstract. A possible role of boiling of the H₂O-NaCl-CO₂(-H₂S) fluid in ore deposition has been examined numerically by using the equations of state (EOS) of Duan et al. (1995, 1996), a modified EOS of Bowers and Helgeson (1983) and the water-rock interaction simulator MIX99 (Hoshino et al., 2000).
The following three models are examined to evaluate an efficiency of boiling on mineral precipitation: (1) hypothetical non-boiling process, (2) hypothetical boiling process with sulfur partitioned only in liquid phases and (3) boiling process in which partition ratios of H₂S between liquid and vapor phases are assumed to be the same as those of CO₂. The processes are simulated from 450C and 900 bar to 310C and 620 bar with an analytical step of 10C / 20 bar. Boiling occurs below 400C in the latter two processes when the initial composition (in mole fraction) of the fluid is: X_H2o = 0.84, X_NaC1 = 0.10 and X_Co2 = 0.06.
Ore deposition occurs abruptly at a boiling point when the partition ratios of total sulfur (X_s^vap/X_s^liq) are as high as those of total carbon during boiling. A decrease of concentration of sulfur in the liquid phase during boiling leads to an increase of pH of the solution, resulting in propelling mineral precipitation. It has been made clear that a possible role of boiling in ore formation mainly depends on the partition ratios of sulfur between the liquid and vapor phases, although they cannot be estimated accurately by the currently available EOS.     

Spatial variability of chemical tracers in surface snow along the traverse route from the coast to 1000 km inland at east Dronning Maud Land, Antarctica

1　IntroductionAntarcticicesheetisaburialgroundforatmosphericdeposition .Sincethereiscon tinuousinteractionbetweentheicesheetandtheatmosphere,variousatmosphericsub stancesareinjectedtotheicesheetsequentiallyintimeandspace.Therefore,verticalanalysesoftheicesheetprovideuswithinformationaboutpastclimaticchange (Delmas1 992 ;LegrandandMayewski 1 997)andhorizontalanalysesoftheicesheetprovideuswithknowledgeregardinglong rangetransportofairbornematerials (Kamiyamaetal.1 989;KreutzandMayewski 1 999)…     

A Study on the Mineralization of the Woxi Au–Sb–W Deposit,Western Hunan,China

The Woxi Au–Sb–W deposit in the western Hunan Province, China, is of hydrothermal vein type characterized by a rare mineral assemblage of stibnite, scheelite and native gold, of which gold fineness ranges from 998.6 to 1000. The mineralization sequence observed in the deposit is, from early to late, coarse‐grained pyrite – scheelite – stibnite – Pb–Sb–S minerals – sphalerite (+ cubanite) – fine‐grained pyrite. Native gold may have precipitated with scheelte. Microthermometric and LA–ICP–MS analyses of fluid inclusions in scheelite, quartz associated with scheelite and stibnite and barren quartz clarified that there may be at least three types of hydrothermal fluids during the vein formation in the Woxi deposit. Scheelite and native gold precipitated from the fluid of high temperature and salinity with high concentrations of metal elements, followed by stibnite precipitation. The later fluid of the highest temperature and salinity with low concentrations of the elements yielded the sphalerite mineralization. The latest fluid of low temperature and salinity with low concentrations of the elements is observed mainly in barren quartz. The remarkably high Au/Ag concentration ratios determined in the fluid inclusions in scheelite might be the reason for the extremely high gold fineness of native gold.     

Assessment of future precipitation indices in the Shikoku region using a statistical downscaling model

In this study, we used the statistical downscaling model (SDSM) to estimate mean and extreme precipitation indices under present and future climate conditions for Shikoku, Japan. Specifically, we considered the following mean and extreme precipitation indices: mean daily precipitation, R10 (number of days with precipitation >10 mm/day), R5d (annual maximum precipitation accumulated over 5 days), maximum dry-spell length (MaDSL), and maximum wet-spell length (MaWSL). Initially, we calibrated the SDSM model using the National Center for environmental prediction (NCEP) reanalysis dataset and daily time series of precipitation for ten locations in Shikoku which were acquired from the surface weather observation point dataset. Subsequently, we used the validated SDSM, using data from NCEP and outputs form general circulation models (GCM), to predict future precipitation indices. Specifically, the HadCM3 GCM was run under the special report on emissions scenarios (SRES) A2 and B2 scenarios, and the CGCM3 GCM was run under the SRES A2 and A1B scenarios. The results showed that: (1) the SDSM can reasonably be used to simulate mean and extreme precipitation indices in the Shikoku region; (2) the values of annual R10 were predicated to decrease in the future in northern Shikoku under all climate scenarios; conversely, the values of annual R10 were predicted to increase in the future in the range of 0–15 % in southern and western Shikoku. The values of annual MaDSL were predicted to increase in northern Shikoku, and the values of annual MaWSL were predicted to decrease in northeastern Shikoku; (3) the spatial variation of precipitation indices indicated the potential for an increased occurrence of drought across northeastern Shikoku and an increased occurrence of flood events in the southwestern part of Shikoku, especially under the A2 and A1B scenarios; (4) characteristics of future precipitation may differ between the northern and southern sides of the Shikoku Mountains. Regional variations in extreme precipitation indices were not notably evident in the B2 scenario compared to the other scenarios.     

Stochastic continuum modeling of flow and transport in a crystalline rock mass: Fanay-Augères,France, revisited

A stochastic discrete-fracture model was used by Cacas et al.^a,b to interpret flow measurements and transport experiments in a fractured crystalline rock mass at Fanay-Augères. They considered continuum models to be incapable of properly interpreting small-scale measurements or tracer tests in fractured systems, which, in their view, require three-dimensional modeling of numerous discrete channels; in their opinion, continuum modeling applies only to average flow on a relatively large scale. Cacas et al. considered their discrete fracture model to have been validated by its demonstrated ability to reproduce selected experimental results. In this paper, flow and transport at Fanay-Augères are modeled by viewing the fractured rock as a stochastic continuum in a manner originally proposed by Neuman^c,d. The stochastic continuum approach obviates the need for detailed information about fracture geometry or assumptions about how individual fractures control flow and transport. All it requires is the delineation of a few dominant features, which can be embedded into the stochastic continuum model as heterogeneous porous slabs. Though a fault zone has been identified at the Fanay-Augères experimental site, it has been modeled neither by Cacas et al. nor in this paper. In fact, in this paper, a larger selection of experimental results than those considered by Cacas et al. are reproduced merely by modeling the rock as a statistically homogeneous continuum in two dimensions. These results demonstrate that a continuum approach may be well suited for the analysis of flow and transport in fractured rock. This does not constitute a validation of the continuum approach, just as the results of Cacas et al. fall short of validating the discrete fracture approach. Instead, the two sets of results illustrate jointly the well-established principle that an open system, especially one as complex as fractured hydrogeologic environments tend to be, cannot be described uniquely on the basis of sparse data and need not be described in great detail to capture its salient behavior by a model.^a Cacas MC, Ledoux E, de Marsily G, Barbreau A, Calmels P, Gaillard B, Margritta R (1990a) Modelling fracture flow with a stochastic discrete fracture network: calibration and validation. 1. The flow model. Water Resour Res 26(3):479–489^b Cacas MC, Ledoux E, de Marsily G, Barbreau A, Calmels P, Gaillard B, Margritta R (1990b) Modelling fracture flow with a stochastic discrete fracture network: calibration and validation. 2. The transport model. Water Resour Res 26(3):491–500^c Neuman SP (1987) Stochastic continuum representation of fractured rock permeability as an alternative to the REV and fracture network concepts, in Rock Mechanics. In: Farmer IW, Daemen JJK, Desai CS, Glass CE, Neuman SP (eds) Proceedings of the 28th U.S. Symposium, Tucson, Arizona. Balkema, Rotterdam, pp 533–561^d Neuman SP (1988) A proposed conceptual framework and methodology for investigating flow and transport in Swedish crystalline rocks. SKB Swedish Nuclear Fuel and Waste Management Co., Stockholm, September, Arbetsrapport 88–37, 39 pp     

湖南康家湾矿床的成矿流体特征:基于矿物学、流体包裹体及热力学计算的制约

湖南康家湾铅锌金银矿床位于南岭成矿带北部中段,是水口山矿田内发现较晚的大型隐伏矿床.该矿床的矿物组合及矿化特征复杂,前人对其成矿流体特征及成因类型存在不同认识.文章通过野外地质调查和矿物矿相学研究发现该矿床的热液成矿阶段较多,方铅矿主要形成于早期闪锌矿与晚期闪锌矿之间,而金和银的成矿阶段分别与早期闪锌矿和方铅矿趋于一致...