Micromechanically inspired investigation of cemented granular materials: part II— from experiments to modelling and back

Acta Geotechnica - This paper presents an experimental and analytical/numerical study of the mechanics of cemented granular materials (CGMs). This study incorporates both in situ X-ray tomography...     

A micro–macro investigation of the capillary strengthening effect in wet granular materials

Wet granular materials are three-dimensionally simulated by the discrete element method with water bridges incorporated between particles. The water bridges are simplified as toroidal shapes, and the matric suction is constantly maintained in the material. A comparison with experimental tests in the literature indicates that the toroidal shape approximation may be one of the best choices with high practicability and decent accuracy. Mechanical behaviours of wet granular materials are studied by triaxial tests. Effects of particle size distributions and void ratios are investigated systematically in this study. The hydraulic limit of the pendular state is also discussed. It gives the capillary cohesion function which is not only determined by the degree of saturation but also positively correlated to relative density and particle size polydispersity and inversely proportional to mean particle size. Furthermore, the capillary strengthening effect is also analysed microscopically in aid of the Stress–Force–Fabric relationship, mainly in fabric anisotropy, coordination number and stress transmission pattern, which revealed the micro-mechanisms of the additional effective stress induced by capillary effect.     

Analysis of human risks due to dam-break floods—part 1: a new model based on Bayesian networks

Dam breaks have catastrophic consequences for human lives. This paper presents a new human risk analysis model (HURAM) using Bayesian networks for estimating human risks due to dam-break floods. A Bayesian network is constructed according to a logic structure of loss-of-life mechanisms. The nodes (parameters) and the arcs (inter-relationships) of the network are quantified with historical data, existing models and physical analyses. A dataset of 343 dam-failure cases with records of fatality is compiled for this purpose. Comparison between two existing models and the new model is made to test the new model. Finally, sensitivity analysis is conducted to identify the important parameters that lead to loss of life. The new model is able to take into account a large number of important parameters and their inter-relationships in a systematic structure; include the uncertainties of these parameters and their inter-relationships; incorporate information derived from physical analysis, empirical models and historical data; and update the predictions when information in specific cases is available. The application of this model to the study of human risks in a specific dam-break case is presented in a companion paper.     

Analysis of human risks due to dam break floods—part 2: application to Tangjiashan landslide dam failure

The Tangjiashan landslide dam was formed during the Ms8.0 Wenchuan earthquake in 2008 and posed high risks to 1.2?million people downstream the dam. A human risk analysis model (HURAM) reported in the companion paper is applied to evaluate the human risk in the Tangjiashan landslide dam failure. The characteristics of this landslide dam are introduced first. The breaching parameters in two cases (i.e., the actual case and a high erodibility case) are predicted with a physically based model, and the flood routing processes in these two cases are simulated using numerical analysis. The population at risk downstream of the landslide dam is then obtained based on the results of the flood routing simulations. Subsequently, the human risks are analyzed with HURAM using Bayesian networks. Fourteen influence parameters and their interrelationships are considered in a systematic structure in the case study. A change in anyone of them may affect the other parameters and leads to loss of life. HURAM allows not only cause-to-result inference, but also result-to-cause inference by updating the Bayesian network with specific information from the study case. The uncertainties of the parameters and their relationships are studied both at the global level using multiple sources of information and at the local level by updating the prior probabilities.     

Flash flood in the mountainous region of Rio de Janeiro state (Brazil) in 2011: part I—calibration watershed through hydrological SMAP model

Natural Hazards - Flooding and water logging possess severe hazards to human population in many parts of the world. Mumbai, the study area, is one of the cities in India and has been frequently...     

On the possibility of determining rates of removal of uranium from crystalline igneous rocks using U-series disequilibria—2: applicability of a U-leach model to mineral separates

Preliminary U-series data are presented on the disequilibria of mineral separates of surface samples from the Eye-Dashwa lakes granite (Canada). Most analyses of the ferromagnesian, quartz-plus-feldspar, and zircon-plus-sphene grains, show that leaching of U has taken place such that²³⁴U/²³⁹U≤ 1 and ²³⁰Th/²³⁴U 1. The simple leach model presented in a companion paper was used to explain these characteristic ratios. It is suggested that the leach model may be used on zircon and sphere grains to estimate U-leach rates, and this is relevant to the suggestion of using sphene-like hosts as the primary barrier in the safe disposal of nuclear waste. A full analysis of U mobility within a given rock volume was not possible due to our limited data base and to our inability to monitor labile U directly.     

Open-source support toward validating and falsifying discrete mechanics models using synthetic granular materials—Part I: Experimental tests with particles manufactured by a 3D printer

Acta Geotechnica - This article presents a new test prototype that leverages the 3D printing technique to create artificial particle assembles to provide auxiliary evidences that supports the...     

Riverbank filtration in Cairo,Egypt: part II—detailed investigation of a new riverbank filtration site with a focus on manganese

A 5-day detailed field investigation at a new RBF test well gallery in Embaba, Cairo, was conducted to evaluate the hydraulic setting and the behavior of iron and manganese. The well gallery consists of six vertical wells placed along a straight line parallel to the Nile riverbank. A low anisotropy factor for the aquifer (k_f,h:k_f,v) of 1.7 was determined by evaluation of a multi-step pumping test. Travel times between 11 days from the river toward the central wells and 22 days toward the outermost wells were estimated by groundwater flow modeling and particle tracking. The riverbed is rich in fine suspended sediments that have elevated iron and nitrogen concentrations. Depth-dependent water sampling during regular well operation indicates that the thick organic-, Fe- and Mn-rich riverbed is the primary source for ammonium, iron and manganese in the bank filtrate. Iron-rich groundwater flow from the opposite riverbank was identified as a secondary source of iron in the pumped water. The vertical position of the filter screen affects total travel times but would not reduce the portion of Mn-rich bank filtrate. The authors recommend continuous well operation for achieving stable water quality and lowering the risk of well clogging.     

Intraplate earthquakes and their link with mantle dynamics: Insights from P-wave teleseismic tomography along the northern part of the North–South Tectonic Zone in China

The North–South Tectonic Zone (NSTZ) running across the Chinese continent is an important earthquake-prone zone. Around one third of the strong earthquakes (> 7.0) of China in the past occurred in this region. Receiver function study has imaged vertical convection in the mantle beneath the northern part of the NSTZ (NNSTZ), which might be related to stress accumulation and release as well as related earthquakes. Here we perform a P-wave teleseismic tomographic analysis of this region. Our results reveal prominent low-velocity and high-velocity perturbations in the upper mantle beneath this region, which we correlate with mantle upwelling, possibly resulting from lower crustal and (or) lithospheric delamination. Our results also reveal significant contrast in the velocity perturbation of the lithosphere along the two sides of this tectonic zone, suggesting possible material exchange between the eastern and western domains and lithosphere-scale control on the generation of earthquakes.     

An integrated zircon geochronological and geochemical investigation into the Miocene plutonic evolution of the Cyclades, Aegean Sea, Greece: part 2—geochemistry

Zircons from 14 compositionally variable granitic rocks were examined in detail using CL image-guided micro-analysis to unravel the complex magmatic history above the southward retreating Hellenic subduction zone system in the Aegean Sea. Previously published U–Pb ages document an episodic crystallisation history from 17 to 11?Ma, with peraluminous (S-type) granitic rocks systematically older than closely associated metaluminous (I-type) granitic rocks. Zircon O- and Hf isotopic data, combined with trace element compositions, are highly variable within and between individual samples, indicative of open-system behaviour involving mantle-derived melts and evolved supracrustal sources. Pronounced compositional and thermal fluctuations highlight the role of magma mixing and mingling, in accord with field observations, and incremental emplacement of distinct melt batches coupled with variable degrees of crustal assimilation. In the course of partial fusion, more fertile supracrustal sources dominated in the earlier stages of Aegean Miocene magmatism, consistent with systematically older crystallisation ages of peraluminous granitic rocks. Differences between zircon saturation and crystallisation temperatures (deduced from zircon Ti concentrations), along with multimodal crystallisation age spectra for individual plutons, highlight the complex and highly variable physico-compositional and thermal evolution of silicic magma systems. The transfer of heat and juvenile melts from the mantle varied probably in response to episodic rollback of the subducting lithospheric slab, as suggested by punctuated crystallisation age spectra within and among individual granitic plutons.     

Effect of land surface processes on the Tibetan Plateau’s past and its predicted response to global warming: an analytical investigation based on simulation results from the CMIP5 model

Complex interactions between the land surface and atmosphere and the exchange of water and energy have a significant impact on climate. The Tibetan Plateau is the highest plateau in the world and is known as “Earth’s third pole”. Because of its unique natural geographical and climatic characteristics, it directly affects China’s climate, as well as the world’s climate, through its thermal and dynamic roles. In this study, the BCCCSM1.1 model for the simulation results of CMIP5 is used to analyze the variation of the land surface processes of the Tibetan Plateau and the possible linkages with temperature change. The analysis showed that, from 1850 to 2005, as temperature increases, the model shows surface downward short-wave radiation, upward short-wave radiation, and net radiation to decrease, and long-wave radiation to increase. Meanwhile, latent heat flux increases, whereas sensible heat flux decreases. Except for sensible heat flux, the correlation coefficients of land surface fluxes with surface air temperature are all significant at the 99 % significance level. The model results indicate rising temperature to cause the ablation of ice (or snow) cover and increasing leaf area index, with reduced snowfall, together with a series of other changes, resulting in increasing upward and downward long-wave radiation and changes in soil moisture, evaporation, latent heat flux, and water vapor in the air. However, rising temperature also reduces the difference between the surface and air temperature and the surface albedo, which lead to further reductions of downward and upward short-wave radiation. The surface air temperature in winter increases by 0.93 °C/100 years, whereas the change is at a minimum (0.66 °C/100 years) during the summer. Downward short-wave and net radiation demonstrate the largest decline in the summer, whereas upward short-wave radiation demonstrates its largest decline during the spring. Downward short-wave radiation is predominantly affected by air humidity, followed by the impact of total cloud fraction. The average downward short-wave and net radiation attain their maxima in May, whereas for upward short-wave radiation the maximum is in March. The model predicts surface temperature to increase under all the different representative concentration pathway (RCP) scenarios, with the rise under RCP8.5 reaching 5.1 °C/100 years. Long-wave radiation increases under the different emission scenarios, while downward short-wave radiation increases under the low- and medium-emission concentration pathways, but decreases under RCP8.5. Upward short-wave radiation reduces under the various emission scenarios, and the marginal growth decreases as the emission concentration increases.     

The geodynamic evolution of the Italian South Alpine basement from the Ediacaran to the Carboniferous: Was the South Alpine terrane part of the peri-Gondwana arc-forming terranes?

Zircon U-Pb LA-ICPMS ages were obtained from three metasedimentary and two metavolcanic samples from the Monte Cavallino (South Tyrol) and the Cima Vallona (Carnic Alps) tectono-metamorphic groups from the eastern South Alpine crystalline basement in NE Italy. These analyses were performed to constrain the maximum depositional ages of the South Alpine domain, and to compare the spatial and temporal provenance variations with those of adjacent terranes. The detrital zircon dataset from the metasedimentary rocks (416 grains) yield populations with age peaks at 2.7–2.9 Ga, 1.8–2.1 Ga, 1.2–0.85 Ga, and 0.65–0.45 Ga, with maximum depositional ages ranging from the latest Neoproterozoic to the Silurian. The metavolcanic zircon dataset (209 grains) documents the presence of a two Ordovician volcanic events in the South Alpine domain. The detrital zircon dataset implies that the clastic units of the South Alpine crystalline basement were (a) deposited on the peri-Gondwanan active continental margin and (b) were originally sourced from the Proterozoic and Paleozoic units of NW Gondwana and hence should no longer be considered as exotic elements. The age spectra of the three metasedimentary units highlight differences between the Ediacaran basement gneiss, the Ordovician greywacke, and the Silurian metaconglomerate, suggesting up-section age variations due to a temporal change in provenance. Collectively, these new detrital zircon U-Pb ages imply that the clastic units within the South Alpine domain recorded sedimentation at c. 550 Ma on the peri-Gondwanan active continental margin, followed by rift-related continental and marine sedimentation in a back-arc basin setting until at least the Silurian. The South Alpine domain ultimately rifted off from Gondwana due to back-arc spreading, and subsequently underwent Variscan metamorphism during accretion onto the Laurussia margin, which started at c. 380 Ma and lasted until at least c. 320 Ma.     

An experimental database for the development,calibration and verification of constitutive models for sand with focus to cyclic loading: part I—tests with monotonic loading and stress cycles

For numerical studies of geotechnical structures under earthquake loading, aiming to examine a possible failure due to liquefaction, using a sophisticated constitutive model for the soil is indispensable. Such a model must adequately describe the material response to a cyclic loading under constant volume (undrained) conditions, amongst others the relaxation of effective stress (pore pressure accumulation) or the effective stress loops repeatedly passed through after a sufficiently large number of cycles (cyclic mobility, stress attractors). The soil behaviour under undrained cyclic loading is manifold, depending on the initial conditions (e.g. density, fabric, effective mean pressure, stress ratio) and the load characteristics (e.g. amplitude of the cycles, application of stress or strain cycles). In order to develop, calibrate and verify a constitutive model with focus to undrained cyclic loading, the data from high-quality laboratory tests comprising a variety of initial conditions and load characteristics are necessary. The purpose of these two companion papers was to provide such database collected for a fine sand. The database consists of numerous undrained cyclic triaxial tests with stress or strain cycles applied to samples consolidated isotropically or anisotropically. Monotonic triaxial tests with drained or undrained conditions have also been performed. Furthermore, drained triaxial, oedometric or isotropic compression tests with several un- and reloading cycles are presented. Part I concentrates on the triaxial tests with monotonic loading or stress cycles. All test data presented herein will be available from the homepage of the first author. As an example of the examination of an existing constitutive model, the experimental data are compared to element test simulations using hypoplasticity with intergranular strain.     

SPINMELT-2.0: Simulation of spinel–melt equilibrium in basaltic systems under pressures up to 15 kbar: I. model formulation,calibration, and tests

The paper presents results of testing currently used models proposed to describe Cr-spinel–melt equilibrium: models of the MELTS family by M.S. Ghiorso with colleagues, the SPINMELT program by A.A. Ariskin and G.S. Nikolaev, and the “MELT–CHROMITE spinel calculator” by А.А. Pustovetov and R.L. Roeder. The new calibration of the SPINMELT model presented in this publication enables calculating a sixcomponent (Mg, Fe²⁺, Cr, Al, Fe³⁺, and Ti) composition of Cr-spinel and the \(T - {f_{{o_2}}}\) parameters of its stability on the liquidus of basaltic melts under pressures up to 15 kbar. The model is based on results of 392 runs from 43 experimental studies, including systems of normal alkalinity at \({f_{{o_2}}}\) ≤ QFM + 2. The experimental dataset (which was extended compared to that used for the previous calibration) allowed us not only to estimate the pressure effect, but also apply the model to aluminous and hydrous systems. Tests of the SPINMELT-2.0 model show that the errors of the calculated temperature of the spinel–melt equilibrium increase with pressure from 16°C at 1 atm to 50°C at 15 kbar. Experimental spinel compositions are reproduced by the model accurate to < 3 at % Al and Cr, and no worse 1 at % for the other cations.     

A simplified conception model for salt marsh sedimentation rate calculation: Response to changing suspended sediment concentration——A case study of Wanggang salt marsh, Jiangsu Province, China

Situated at the junction of land and sea, coastal salt marshes are very sensitive to global environmental changes. Existing numerical models simulating the vertical growth of coastal salt marshes are under the conditions of low suspended sediment concentration (SSC), but under high SSC, such as in the Jiangsu coast of China, the models need to be refined. In this paper, based on the mass balance principle, a refined zero-dimensional time-stepped model is applied after considering the average suspended sediment concentration (ASSC). The statistical results of in situ observation data show that there is an obvious positive linear relationship between ASSC and maximum inundation height, and an obvious positive exponential one between inundation time and inundation height. Using the short-term measurement values for the input parameters and incorporating the above-mentioned relationship between SSC and inundation height, the sedimentation rates are calculated, and the results are similar to previous data. In addition, some factors and hydrodynamic processes that were ignored or simplified could influence the calculated value.     

An experimental database for the development,calibration and verification of constitutive models for sand with focus to cyclic loading: part II—tests with strain cycles and combined loading

For numerical studies of geotechnical structures under earthquake loading, aiming to examine a possible failure due to liquefaction, using a sophisticated constitutive model for the soil is indispensable. Such model must adequately describe the material response to a cyclic loading under constant volume (undrained) conditions, amongst others the relaxation of effective stress (pore pressure accumulation) or the effective stress loops repeatedly passed through after a sufficiently large number of cycles (cyclic mobility, stress attractors). The soil behaviour under undrained cyclic loading is manifold, depending on the initial conditions (e.g. density, fabric, effective mean pressure, stress ratio) and the load characteristics (e.g. amplitude of the cycles, application of stress or strain cycles). In order to develop, calibrate and verify a constitutive model with focus to undrained cyclic loading, the data from high-quality laboratory tests comprising a variety of initial conditions and load characteristics are necessary. It is the purpose of these two companion papers to provide such database collected for a fine sand. Part II concentrates on the undrained triaxial tests with strain cycles, where a large range of strain amplitudes has been studied. Furthermore, oedometric and isotropic compression tests as well as drained triaxial tests with un- and reloading cycles are discussed. A combined monotonic and cyclic loading has been also studied in undrained triaxial tests. All test data presented herein will be available from the homepage of the first author. As an example of the examination of an existing constitutive model, the experimental data are compared to element test simulations using hypoplasticity with intergranular strain.     

Using impacts of deep-level mining to research karst hydrology—a Darcy-based approach to predict the future of dried-up dolomitic springs in the Far West Rand goldfield (South Africa). Part 1: a conceptual model of recharge and inter-compartmental flow

Some of the world’s deepest goldmines are located in the Far West Rand (FWR) goldfield operating below of up to 1.2-km-thick dolomites hosting some of the largest karst aquifers in South Africa. Associated impacts include the dewatering of the overlying karst aquifers as well as linking previously disconnected compartments by mining through aquicludes (dykes). The focus of the study is on predicting groundwater balances in re-watered aquifers after mining ceases as this will determine whether or not associated karst springs that dried-up due to dewatering will ever flow again. Critically revisiting, Swart et al. (Environ Geol 44:751–770, 2003a) who predict that all springs will flow again, this study uses significantly larger data sets and modified assumptions to increase the robustness of findings as the question is crucial for post-closure development. As a first of two papers, this part develops a conceptual model on the mega-compartment concept that predicts a flat water table across all linked compartments that would leave the springs dry. The model identifies the ratio between inflowing surface water (recharge) and underground water losses to downstream compartments via mined-through dykes (‘inter-compartmental groundwater flow’, IGF) as a key factor governing the elevation of the post-mining water table, creating the base for part 2, where the IGF and the post-mining water tables are determined using unique large data sets that have not been evaluated before.     

Determination of Fe<Superscript>3+</Superscript>/ΣFe ratios in chrome spinels using a combined Mössbauer and single-crystal X-ray approach: application to chromitites from the mantle section of the Oman ophiolite

We present the results of a comparative study in which we have measured Fe³⁺/ΣFe ratios in chromites from mantle chromitites in the Oman ophiolite using Mössbauer spectroscopy and single-crystal X-ray diffraction. We have compared these results with ratios calculated from mineral stoichiometry and find that mineral stoichiometry calculations do not accurately reflect the measured Fe³⁺/ΣFe ratios. We have identified three groups of samples. The majority preserve Fe³⁺/ΣFe ratios which are thought to be magmatic, whereas a few samples are highly oxidized and have high Fe³⁺/ΣFe ratios. There is also a group of partially oxidized samples. The oxidized chromites show anomalously low cell edge (a ₀) values and their oxygen positional parameters among the lowest ever found for chromites. Site occupancy calculations show that some chromites are non-stoichiometric and contain vacancies in their structure randomly distributed between both the T and M sites. The field relationships suggest that the oxidation of the magmatic chromitites took place in association with a ductile shear zone in mantle harzburgites. Primary magmatic Fe³⁺/ΣFe ratios measured for the Oman mantle chromitites are between 0.193–0.285 (X-ray data) and 0.164–0.270 (Mössbauer data) and preserve a range of Fe³⁺/ΣFe ratios which we propose is real and reflects differences in the composition of the magmas parental to the chromitites. The range of values extends from those MORB melts (0.16 ± 0.1) to those for arc basalts (0.22–0.28).