Numerical runout simulation of debris avalanches in the Faroe Islands, North Atlantic Ocean

The Faroe Islands in the North Atlantic Ocean are susceptible to flow-type landslides in coarse-grained highly organic colluvium. Following several hazardous debris avalanche events, research work has been initiated to quantify landslide risk. A central task in this work is to predict landslide runout behavior. From numerical simulation of four debris avalanches, this study provides a first screening of which rheology and appertaining input parameters best predict runout behavior of debris avalanches in the Faroe Islands. Three rheologies (frictional, Voellmy, and Bingham) are selected and used for individual back analysis of the events in the numerical models BING and DAN3D. A best fit rheology is selected from comparing predicted and observed landslide runout behavior. General back analysis to identify the optimal input parameters for the chosen rheology is performed by cross validation, where each debris avalanche is modeled with input parameters from the three other events. Optimal input parameters are found from the model run producing the most accurate runout length and velocity. The Bingham is selected as the best fit rheology, a result differing from similar studies of coarse-grained landslides. A reason for why particularly the frictional rheology proves unsuitable is its tendency to produce too long runout lengths of the low-weight runout material, a result showing important limitations for using the frictional rheology in DAN3D. Optimal Bingham input parameters are τ _y ?=?980 Pa and μ _b ?=?117 Pa/s. However, future studies performed in 2D models are needed for precise parameterization before results can be used for landslide risk assessment.     

From back-analysis to run-out prediction: a case study in the Western Italian Alps

Rock avalanches are complex phenomena that occur with a low frequency but which have a high destructive potential. As a consequence, the people who are responsible for the management of a territory are more and more interested in predicting the possible evolutions of well-known potential events. Tackling the above problems from a quantitative point of view, the RASH3D code, based on continuum mechanics concepts, has been here used to predict the evolution of a potential rock avalanche in the Western Italian Alps. A calibration-based approach, in which rheological parameters are constrained by systematic adjustment during trial-and-error back-analysis of past events similar to the landslide under investigation, is proposed to set rheological parameter values to be used for prediction purposes. The back-analysis of a 2?10⁶ m³ rock avalanche located in the Divedro Valley, close to the area of the potential event, has then been analysed using both a frictional and a Voellmy rheology. The characteristics of the slope and the dynamics of the event have made the frictional rheology more suitable to come to the correct simulation of the historical case. The back-analysis results have contributed not only in the selection of the rheological parameter values but also in the choice of the type of rheological law to use in the carried out forward-analyses.     

Characteristics and numerical runout modelling of a catastrophic rock avalanche triggered by the Wenchuan earthquake in the Wenjia valley,Mianzhu, Sichuan,China

The 2008 Wenchuan earthquake triggered more than 100 rock avalanches with volumes greater than 10 million cubic metres. The rock avalanche with the longest runout amongst these destructive landslides occurred in the Wenjia valley, Mianzhu, Sichuan, China. The landslide involved the failure of about 27.5 million cubic metres of sandstone from the source area. The displaced material travelled about 4,170 m with an elevation descent of about 1,360 m, equivalent to a fahrböschung of 16.9° and covered an area of 1.5 million square metres, with the final deposited volume of approximately 49 million cubic metres. The catastrophic event destroyed the village of Yanjing, killed 48 people and buried some houses at the mouth of the Wenjia valley. On the basis of a detailed field investigation, we introduce basic characteristics of the rock avalanche and find that the rock avalanche resulted in two run-ups and a superelevation along the runout path, and downslope enlargement due to the entrainment of path materials. A numerical model (DAN3D) is used to simulate the post-failure behaviour of the rock avalanche. By means of trial and error, a combination of the frictional model and Voellmy model is found to provide the best performance in simulating this rock avalanche. The simulation results reveal that the rock avalanche had a duration of about 240 s and an average velocity of 17.4 m/s.     

A non-hydrostatic model for the numerical study of landslide-generated waves

This paper proposes and demonstrates a two-layer depth-averaged model with non-hydrostatic pressure correction to simulate landslide-generated waves. Landslide (lower layer) and water (upper layer) motions are governed by the general shallow water equations derived from mass and momentum conservation laws. The landslide motion and wave generation/propagation are separately formulated, but they form a coupled system. Our model combines some features of the landslide analysis model DAN3D and the tsunami analysis model COMCOT and adds a non-hydrostatic pressure correction. We use the new model to simulate a 2007 rock avalanche-generated wave event at Chehalis Lake, British Columbia, Canada. The model results match both the observed distribution of the rock avalanche deposit in the lake and the wave run-up trimline along the shoreline. Sensitivity analyses demonstrate the importance of accounting for the non-hydrostatic dynamic pressure at the landslide-water interface, as well as the influence of the internal strength of the landslide on the size of the generated waves. Finally, we compare the numerical results of landslide-generated waves simulated with frictional and Voellmy rheologies. Similar maximum wave run-ups can be obtained using the two different rheologies, but the frictional model better reproduces the known limit of the rock avalanche deposit and is thus considered to yield the best overall results in this particular case.     

Comparison of methods for the estimation of pyrite oxidation rate in a waste rock pile at Mine Doyon site, Quebec, Canada

Several methods were evaluated and compared for the estimation of pyrite oxidation rates (POR) in waste rock at Mine Doyon, Quebec, Canada. Methods based on data collected in situ, such as the interpretation of temperature and oxygen concentration profiles (TOP) measured in the waste rock pile and pyrite mass balance (PMB) on solid phase samples were compared with the oxygen consumption measurements (OCM) in closed chamber in the laboratory. A 1-D analytical solution to a gas and heat transport equation used temperature and oxygen profiles (TOP) measured in the pile for the preliminary POR estimates at a site close to the slope of the pile (Site 6) and in the core of the pile (Site 7). Resulting POR values were 1.1 × 10^− 9 mol(O₂) kg^− 1 s^− 1 and 1.0 × 10^− 10 mol(O₂) kg^− 1 s^− 1 for the slope site and the core site, respectively. Oxidation rates based on pyrite mass balance (PMB) calculations for solid samples were 2.21 × 10^− 9 mol(O₂) kg^− 1 s^− 1 and 2.03 × 10^− 9 mol(O₂) kg^− 1 s^− 1, respectively, for the same slope and core sites, but the difference between sites was within the error margin. The OCM measurements in the laboratory on fresh waste rock samples yielded higher POR values than field methods, with average oxidation rate of 6.7 × 10^− 8 mol(O₂) kg^− 1 s^− 1. However, the OCM results on weathered and decomposed material from the rock stockpile (average oxidation rate 3.4 × 10^− 9 mol(O₂) kg^− 1 s^− 1) were consistent with results from the field-based estimates. When POR values based on fresh material are excluded, the remaining POR values for all methods range from 1.0 × 10^− 10 to 3.4 × 10^− 9 mol(O₂) kg^− 1 s^− 1. The lowest estimated value (1.0 × 10^− 10 mol(O₂) kg^− 1 s^− 1) was based on TOP estimates in the interior of the pile where oxygen transport was limited by diffusion from the surface. These results suggest that small-scale OCM laboratory experiments may provide relatively representative values of POR in the zones of waste rock piles in which oxygen transport is not dominated by diffusion.     

The Zymoetz River landslide, British Columbia, Canada: description and dynamic analysis of a rock slide–debris flow

The Zymoetz River landslide is a recent example of an extremely mobile type of landslide known as a rock slide–debris flow. It began as a failure of 900,000 m³ of bedrock, which mobilized an additional 500,000 m³ of surficial material in its path, transforming into a large debris flow that traveled over 4 km from its source. Seasonal snow and meltwater in the proximal part of the path were important factors. A recently developed dynamic model that accounts for material entrainment, DAN3D, was used to back-analyze this event. The two distinct phases of motion were modeled using different basal rheologies: a frictional model in the proximal path and a Voellmy model in the distal path, following the initiation of significant entrainment. Very good agreement between the observed and simulated results was achieved, suggesting that entrainment capabilities are essential for the successful simulation of this type of landslide.     

Dynamic recrystallization near the brittle-plastic transition in naturally and experimentally deformed quartz aggregates

The electron backscattering diffraction technique (EBSD) was used to analyze bulging recrystallization microstructures from naturally and experimentally deformed quartz aggregates, both of which are characterized by porphyroclasts with finely serrated grain boundaries and grain boundary bulges set in a matrix of very fine recrystallized grains. For the Tonale mylonites we investigated, a temperature range of 300–380 °C, 0.25 GPa confining pressure, a flow stress range of ~ 0.1–0.2 GPa, and a strain rate of ~ 10^− 13 s^− 1 were estimated. Experimental samples of Black Hills quartzite were analyzed, which had been deformed in axial compression at 700 °C, 1.2–1.5 GPa confining pressure, a flow stress of ~ 0.3–0.4 GPa, a strain rate of ~ 10^− 6 s^− 1, and to 44% to 73% axial shortening. Using orientation imaging we investigated the dynamic recrystallization microstructures and discuss which processes may contribute to their development. Our results suggest that several deformation processes are important for the dismantling of the porphyroclasts and the formation of recrystallized grains. Grain boundary bulges are not only formed by local grain boundary migration, but they also display a lattice misorientation indicative of subgrain rotation. Dynamic recrystallization affects especially the rims of host porphyroclasts with a hard orientation, i.e. with an orientation unsuitable for easy basal slip. In addition, Dauphiné twins within porphyroclasts are preferred sites for recrystallization. We interpret large misorientation angles in the experimental samples, which increase with increasing strain, as formed by the activity of fluid-assisted grain boundary sliding.     

Investigation and dynamic analysis of a catastrophic rock avalanche on September 23, 1991, Zhaotong,China

At 6:10 p.m. on September 23, 1991, a catastrophic rock avalanche occurred in Zhaotong, Yunnan, southwestern China. Over 216 people were killed when the Touzhai village was overwhelmed directly in the path of the landslide. The landslide involved the failure of about 12 Mm³ of jointed basaltic rock mass from the source area. The displaced materials ran out a horizontal distance of 3650 m over a vertical distance of 960 m, equivalent to a Fahrböschung of 14.7°, and covered an area of 1.38 km². To provide information for hazard zonation of similar type of potential landslides in the same area, we used a dynamic model (DAN-W) with three alternative rheological models to simulate the runout behaviour of the displaced landslide materials and found that a combination of the frictional model and Voellmy model could provide the best performance in simulating this landslide. The simulated results indicated that the duration of the movement is estimated at about 175 s for a mean velocity 21 m/s.     

Crustal structure at the easternmost termination of the Variscan belt based on CELEBRATION 2000 and ALP 2002 data

The eastern margin of the Variscan belt in Europe comprises plate boundaries between continental blocks and terranes formed during different tectonic events. The crustal structure of that complicated area was studied using the data of the international refraction experiments CELEBRATION 2000 and ALP 2002. The seismic data were acquired along SW–NE oriented refraction and wide-angle reflection profiles CEL10 and ALP04 starting in the Eastern Alps, passing through the Moravo-Silesian zone of the Bohemian Massif and the Fore-Sudetic Monocline, and terminating in the TESZ in Poland. The data were interpreted by seismic tomographic inversion and by 2-D trial-and-error forward modelling of the P waves. Velocity models determine different types of the crust–mantle transition, reflecting variable crustal thickness and delimiting contacts of tectonic units in depth. In the Alpine area, few km thick LVZ with the Vp of 5.1 km s^− 1 dipping to the SW and outcropping at the surface represents the Molasse and Helvetic Flysch sediments overthrust by the Northern Calcareous Alps with higher velocities. In the Bohemian Massif, lower velocities in the range of 5.0–5.6 km s^− 1 down to a depth of 5 km might represent the SE termination of the Elbe Fault Zone. The Fore-Sudetic Monocline and the TESZ are covered by sediments with the velocities in the range of 3.6–5.5 km s^− 1 to the maximum depth of 15 km beneath the Mid-Polish Trough. The Moho in the Eastern Alps is dipping to the SW reaching the depth of 43–45 km. The lower crust at the eastern margin of the Bohemian Massif is characterized by elevated velocities and high Vp gradient, which seems to be a characteristic feature of the Moravo-Silesian. Slightly different properties in the Moravian and Silesian units might be attributed to varying distances of the profile from the Moldanubian Thrust front as well as a different type of contact of the Brunia with the Moldanubian and its northern root sector. The Moho beneath the Fore-Sudetic Monocline is the most pronounced and is interpreted as the first-order discontinuity at a depth of 30 km.     

Strength and stability of frictional sliding of gabbro gouge at elevated temperatures

To investigate the strength of frictional sliding and stability of mafic lower crust, we conducted experiments on oven-dried gabbro gouge of 1 mm thick sandwiched between country rock pieces (with gouge inclined 35° to the sample axis) at slip rates of 1.22 × 10^− 3 mm/s and 1.22 × 10^− 4 mm/s and elevated temperatures up to 615 °C. Special attention has been paid to whether transition from velocity weakening to velocity strengthening occurs due to the elevation of temperature.Two series of experiments were conducted with normal stresses of 200 MPa and 300 MPa, respectively. For both normal stresses, the friction strengths are comparable at least up to 510 °C, with no significant weakening effect of increasing temperature. Comparison of our results with Byerlee's rule on a strike slip fault with a specific temperature profile in the Zhangbei region of North China shows that the strength given by experiments are around that given by Byerlee's rule and a little greater in the high temperature range.At 200 MPa normal stress, the steady-state rate dependence a − b shows only positive values, probably still in the “run-in” process where velocity strengthening is a common feature. With a normal stress of 300 MPa, the values of steady-state rate dependence decreases systematically with increasing temperature, and stick-slip occurred at 615 °C. Considering the limited displacement, limited normal stress applied and the effect of normal stress for the temperatures above 420 °C, it is inferred here that velocity weakening may be the typical behaviour at higher normal stress for temperature above 420 °C and at least up to 615 °C, which covers most of the temperature range in the lower crust of geologically stable continental interior. For a dry mafic lower crust in cool continental interiors where frictional sliding prevails over plastic flow, unstable slip nucleation may occur to generate earthquakes.     

Thermal conductivity of Higher Himalayan Crystallines from Garhwal Himalaya, India

We report the measurements of thermal conductivity for some Higher Himalayan Crystalline rocks from Joshimath and Uttarkashi areas of the Garhwal Himalaya. Seventy-three rock samples including gneiss, metabasic rock and quartzite were measured. Gneissic rocks, which include augen gneiss, banded gneiss, felsic gneiss and fine-grained gneiss, exhibit a wide range in conductivity, from 1.5 to 3.6 Wm^− 1K^− 1 for individual samples, and 2.1 to 2.7 Wm^− 1K^− 1 for the means. Among these, augen gneisses and banded gneisses show the largest variability. Of all the rock types, quartzites (mean 5.4 Wm^− 1K^− 1) and metabasic rocks (mean 2.1 Wm^− 1K^− 1) represent the highest and lowest mean values respectively. The range in conductivity observed for gneissic rocks is significantly higher than that generally found in similar rock types in cratonic areas. The rock samples have very low porosity and exhibit feeble anisotropy, indicating that they do not contribute to the variability in thermal conductivity. Besides variations in mineralogical composition, the heterogeneous banding as well as intercalations with metabasic rocks and quartz veins, a common occurrence in structurally complex areas, appears to cause the variability in conductivity. The study therefore brings out the need for systematic characterization of thermophysical properties of major rock types comprising the Himalayan region for lithospheric thermal modeling, assessment of geothermal energy and geo-engineering applications in an area. The dataset constitutes the first systematic measurements on the Higher Himalayan Crystalline rocks.     

High-field cantilever magnetometry as a tool for the determination of the magnetocrystalline anisotropy of single crystals

Cantilever torque magnetometry is utilized widely in physics and material science for the determination of magnetic properties of thin films and semiconductors. Here, we report on its first application in rock magnetism, namely the determination of K₁ and K₂ of single crystal octahedra of natural magnetite. The design of cantilever magnetometers allows optimization for the specific research question at hand. For the present study, a cantilever magnetometer was used that enables measurement of samples with a volume up to 64 mm³. It can be inserted into an electromagnet with a maximum field of 2 T. The cantilever spring is suitable for torque values ranging from 7.5 × 10^− 7 N·m to 5 × 10^− 6 N·m. The torque is detected capacitively; the measured capacitance is converted into torque by using a calibrated feedback coil. The magnetometer allows in-situ rotation of the sample in both directions and is, therefore, also suitable to analyze rotational hysteresis effects.The evaluation of the magnetite anisotropy constants involved Fourier analysis of the torque signal on the magnetite crystals' (001) and (110) planes. The absolute anisotropy constant has been computed using the extrapolation-to-infinite-field method. The value of K₁ at room temperature is determined at − 1.28 × 10⁴ [J m^− 3] (± 0.13, i.e. 10%) and that of K₂ at − 2.8 × 10³ [J m^− 3] (± 0.1, i.e. 2%). These values concur with earlier determinations that could not provide an instrumental error, in contrast with this work.The cantilever magnetometer performs four times faster than other torque magnetometers used for rock magnetic studies. This makes the instrument also suitable for magnetic fabric analysis.     

Insights into Wasatch fault vertical slip rates using the age of sediments in Timpanogos Cave, Utah

Timpanogos Cave, located near the Wasatch fault, is about 357 m above the American Fork River. Fluvial cave sediments and an interbedded carbonate flowstone yield a paleomagnetic and U–Th depositional age of 350 to 780 ka. Fault vertical slip rates, inferred from calculated river downcutting rates, range between 1.02 and 0.46 mm yr^− 1. These slip rates are in the range of the 0–12 Ma Wasatch Range exhumation rate ( 0.5–0.7 mm yr^− 1), suggesting that the long-term vertical slip rate remained stable through mid-Pleistocene time. However, the late Pleistocene (0–250 ka) decelerated slip rate ( 0.2–0.3 mm yr^− 1) and the accelerated Holocene slip rate ( 1.2 mm yr^− 1) are consistent with episodic fault activity. Assuming that the late Pleistocene vertical slip rate represents an episodic slowing of fault movement and the long-term (0–12 Ma) average vertical slip rate, including the late Pleistocene and Holocene, should be  0.6 mm yr^− 1, there is a net late Pleistocene vertical slip deficit of  50–75 m. The Holocene and late Pleistocene slip rates may be typical for episodes of accelerated and slowed fault movement, respectively. The calculated late Pleistocene slip deficit may mean that the current accelerated Wasatch fault slip rate will extend well into the future.     

Organic geochemistry and petrology of barren and Mo–Ni–PGE mineralized marine black shales of the Lower Cambrian Niutitang Formation (South China)

Marine black shales of the Lower Cambrian Niutitang Formation in southern China host Mo–Ni–platinum group elements (PGE) mineralization confined to a phosphate- and pyrite-rich stratiform body (max. 20-cm thick). The H/C atomic ratio, carbon isotopic composition, FTIR spectra of bulk organic matter, and spectra of extractable part of organic matter indicate similar sources and thermal evolution of organic matter in barren and mineralized black shales.The morphology and relative abundance of organic particles in barren and mineralized shales are different. In barren black shales, organic particles comprise only elongated bodies and laminae 2–10 μm across or elongated larger bodies (> 10 μm) with R_max = 2.96–5.21% (Type I particles). Mineralized black shales contain Type I particles in rock matrix (90–95 vol%), small veinlets or irregular organic accumulations (Type II particles, 1–5 vol%) that display weak to well developed mosaic texture and a variable reflectance (R_max = 3.55–8.65%), and small (< 1 to 5 μm) rounded or irregular Type III organic particles (1–4 vol%) distributed within phosphate nodules and sulphide rip-up clasts. Type III particles show similar reflectance as particles of Type I in rock matrix. Type I particles are interpreted as remnants of in situ bacterially reworked organic matter of cyanobacteria/algal type, Type II as solidified products or oil-derived material (migrabitumen), and Type III particles as remnants of original organic matter in phosphatized or sulphidized algal/microbial oncolite-like bodies. Equivalent vitrinite reflectances of Type I and III particles in barren and mineralized rocks are similar and correspond to semi-anthracite and anthracite. Micro-Raman spectra of organic particles in rocks display a wide belt in the area of 1600 cm^− 1 (G belt) and approximately the same belt in the area of 1350 cm^− 1 (D belt). The ratio of integrated areas of the two belts correlate with R_max values.The Mo–Ni–PGE mineralized body is interpreted as to represent a remnant of phosphate- and sulphide-rich subaquatic hardground supplied with organic material derived from plankton and benthic communities as well as with algal/microbial oncolite-like bodies that originated in wave-agitated, shallow-water, nearshore environment.     

Fluid geochemistry in the Ivigtut cryolite deposit, South Greenland

The 1.27 Ga old Ivigtut (Ivittuut) intrusion in South Greenland is world-famous for its hydrothermal cryolite deposit [Na₃AlF₆] situated within a strongly metasomatised A-type granite stock. This detailed fluid inclusion study characterises the fluid present during the formation of the cryolite deposit and thermodynamic modelling allows to constrain its formation conditions.Microthermometry revealed three different types of inclusions: (1) pure CO₂, (2) aqueous-carbonic and (3) saline-aqueous inclusions. Melting temperatures range between − 23 and − 15 °C for type 2 and from − 15 to − 10 °C for type 3 inclusions. Most inclusions homogenise between 110 and 150 °C into the liquid.Stable isotope compositions of CO₂ and H₂O were measured from crushed inclusions in quartz, cryolite, fluorite and siderite. The δ¹³C values of about − 5‰ PDB are typical of mantle-derived magmas. The differences between δ¹⁸O of CO₂ (+ 21 to + 42‰ VSMOW) and δ¹⁸O of H₂O (− 1 to − 21.7‰ VSMOW) suggest low-temperature isotope exchange. δD (H₂O) ranges from − 19 to − 144‰ VSMOW. The isotopic composition of inclusion water closely follows the meteoric water line and is comparable to Canadian Shield brines. Ion chromatography revealed the fluid's predominance in Na, Cl and F. Cl/Br ratios range between 56 and 110 and may imply intensive fluid–rock interaction with the host granite.Isochores deduced from microthermometry in conjunction with estimates for the solidification of the Ivigtut granite suggest a formation pressure of approximately 1–1.5 kbar for the fluid inclusions. Formation temperatures of different types of fluid inclusions vary between 100 and 400 °C. Thermodynamic modelling of phase assemblages and the extraordinary high concentration in F (and Na) may indicate that the cryolite body and its associated fluid inclusions could have formed during the continuous transition from a volatile-rich melt to a solute-rich fluid.     

Application of zeolitic volcanic rocks for arsenic removal from water

Natural zeolitic rocks consisting mainly of chabazite-phillipsite, clinoptilolite, and volcanic glass have been evaluated by means of batch methods to remove arsenic from waters with different mineralization degree (from deionized water to natural water with a specific conductivity of 1,600 μS cm^− 1). Arsenic was previously spiked in the studied waters at concentrations of about 100 µg l^− 1 to simulate actual cases. The compositional range of natural waters is representative of large hydrogeochemical regions around the world. The experiments were focussed on the application of natural common zeolitic rocks to water treatment for human consumption. The removal efficiency observed rises, in the better cases, 60–80% for chabazite-phillipsite raw materials whereas is 40–60% for clinoptilolite-bearing ones. The arsenic removal tends to increase with water mineralization degree, independently of the zeolitic rock type. A large zeolitic content in the chabazite-phillipsite raw materials increase the removal. Instead, the inverse situation is observed in the clinoptilolite-bearing rocks. The relevance of the quantitative mineralogical analysis, determining also the content of volcanic glass, as well as the use of natural waters in the removal tests has been demonstrated.     

Experimental and numerical investigations of a catastrophic long-runout landslide in Zhenxiong,Yunnan, southwestern China

On 11 January 2013, a catastrophic landslide of ～0.2 million m³ due to a prolonged low-intensity rainfall occurred in Zhenxiong, Yunnan, southwestern China. This landslide destroyed the village of Zhaojiagou and killed 46 people in the distal part of its path. The displaced landslide material traveled a horizontal distance of ～800 m with a vertical drop of ～280 m and stopped at 1520 m a.s.l. To examine the possible mechanism and behavior of the landslide from initiation to runout, the shear behavior of soil samples collected from the sliding surface and runout path was examined by means of ring shear tests. The test results show that the shear strength of sample from the sliding surface is less affected by shear rate while the shear rate has a negative effect on the shear strength of runout path material. It is suggested that the source and runout path materials follow the frictional and Voellmy rheology, respectively. Post-failure behavior of the landslide was modeled by using a DAN-W model, and the numerical results show that the selected rheological relationships and parameters based on the results of ring shear tests may provide good performance in modeling the Zhenxiong landslide.     

Study on the kinetics of iron oxide leaching by oxalic acid

The presence of iron oxides in clay or silica raw materials is detrimental to the manufacturing of high quality ceramics. Although iron has been traditionally removed by physical mineral processing, acid washing has been tested as it is more effective, especially for extremely low iron (of less than 0.1% w/w). However, inorganic acids such as sulphuric or hydrochloric acids easily contaminate the clay products with SO₄²⁻ and Cl⁻, and therefore should be avoided as much as possible. On the other hand, if oxalic acid is used, any acid left behind will be destroyed during the firing of the ceramic products. The characteristics of dissolution of iron oxides were therefore investigated in this study.The dissolution of iron oxides in oxalic acid was found to be very slow at temperatures within the range 25–60 °C, but its rate increases rapidly above 90 °C. The dissolution rate also increases with increasing oxalate concentration at the constant pH values set within the optimum range of pH2.5–3.0. At this optimum pH, the dissolution of fine pure hematite (Fe₂O₃) (105–140 μm) follows a diffusion-controlled shrinking core model. The rate expression expressed as 1 − (2 / 3)x − (1 − x)^2 / 3 where x is a fraction of iron dissolution was found to be proportional to [oxalate]^1.5.The addition of magnetite to the leach liquor at 10% w/w hematite was found to enhance the dissolution rate dramatically. Such addition of magnetite allows coarser hematite in the range 0.5–1.4 mm to be leached at a reasonable rate.     

Simulation of Tsaoling landslide, Taiwan, based on Saint Venant equations over general topography

The 1999 Chi-Chi earthquake triggered the catastrophic Tsaoling landslide in Taiwan. The geomorphological change measured from the data of the 1989 and 2000 aerial photos reveals that the scar and deposit volumes are about 0.126 km³ and 0.15 km³ respectively. The debris material ran over a distance of 1.6 km with 500 m descent in elevation. In this paper, we use the continuum model of hydraulic flow, SHALTOP^2D, based on the equations of Bouchut and Westdickenberg to simulate numerically the landslide dynamics. When the mass is moving, the flow is assumed hydrostatic with a basal Coulomb friction. The best fit is obtained using the basal friction angle equal to 6°, the only parameter of the simulation, uniformly applied in the calculation domain. The landslide front reaches the Chinshui river valley, the foothill of the slope, within only about 25 s after initiation and the motion settles in about 113 s. The maximum speed is estimated 75 m/s. The spreading of the deposit agrees well with the field measurement.     

11 Å tobermorite from cement bypass dust and waste container glass: A feasibility study

A novel one-step hydrothermal synthesis of 11 Å tobermorite, a cation exchanger, from a unique combination of waste materials is reported. 11 Å tobermorite was prepared from stoicheiometric quantities of cement bypass dust and waste container glass at 100 °C in water. The product also comprised 10 wt.% calcite and trace quartz as residual parent phases from the cement bypass dust. In a batch sorption study at 20 °C the uptakes of Cd²⁺ and Pb²⁺ by the waste-derived tobermorite product were found to be 171 mg g^− 1 and 467 mg g^− 1, respectively, and in both cases the removal process could be described using a simple pseudo-second-order rate model (k₂ = 2.30 × 10^− 5 g mg^− 1 min^− 1 and 5.09 × 10^− 5 g mg^− 1 min^− 1, respectively). The sorption characteristics of the 11 Å tobermorite are compared with those of other waste-derived sorbents and potential applications are discussed.