Investigating the impact of the properties of pilot points on calibration of groundwater models: case study of a karst catchment in Rote Island,Indonesia

A robust configuration of pilot points in the parameterisation step of a model is crucial to accurately obtain a satisfactory model performance. However, the recommendations provided by the majority of recent researchers on pilot-point use are considered somewhat impractical. In this study, a practical approach is proposed for using pilot-point properties (i.e. number, distance and distribution method) in the calibration step of a groundwater model. For the first time, the relative distance–area ratio (d/A) and head-zonation-based (HZB) method are introduced, to assign pilot points into the model domain by incorporating a user-friendly zone ratio. This study provides some insights into the trade-off between maximising and restricting the number of pilot points, and offers a relative basis for selecting the pilot-point properties and distribution method in the development of a physically based groundwater model. The grid-based (GB) method is found to perform comparably better than the HZB method in terms of model performance and computational time. When using the GB method, this study recommends a distance–area ratio of 0.05, a distance–x-grid length ratio (d/X _grid) of 0.10, and a distance–y-grid length ratio (d/Y _grid) of 0.20.     

Coupled DEM-CFD investigation on the formation of landslide dams in narrow rivers

Large-scale landslide dams can induce significant hazards to human lives by blocking the river flows and causing inundation upstream. They may trigger severe outburst flooding that may devastate downstream areas once failed. Thus, the advancement in understanding the formation of landslide dams is highly necessary. This paper presents 3D numerical investigations of the formation of landslide dams in open fluid channels via the discrete element method (DEM) coupled with computational fluid dynamics (CFD). By employing this model, the influence of flow velocity on granular depositional morphology has been clarified. As the grains settle downwards in the fluid channel, positive excess water pressures are generated at the bottom region, reducing the total forces acting on the granular mass. In the meantime, the particle sedimentations into the fluid channel with high impacting velocities can generate fluid streams to flow backwards and forwards. The coupled hydraulic effects of excess water pressure and fluid flow would entrain the solid grains to move long distances along the channel. For simulations using different flow velocities, the larger the flow velocity is, the further distance the grains can be transported to. In this process, the solid grains move as a series of surges, with decreasing deposit lengths for the successive surges. The granular flux into the fluid channel has very little influence on the depositional pattern of particles, while it affects the particle–fluid interactions significantly. The results obtained from the DEM-CFD coupled simulations can reasonably explain the mechanisms of granular transportation and deposition in the formation of landslide dams in narrow rivers.     

Spatial characteristics of landslides triggered by the 2015 M<Subscript>w</Subscript> 7.8 (Gorkha) and M<Subscript>w</Subscript> 7.3 (Dolakha) earthquakes in Nepal

High magnitude earthquakes trigger numerous landslides and their occurrences are mainly controlled by terrain parameters. We created an inventory of 15,551 landslides with a total area of 90.2 km² triggered by the 2015 M_w 7.8 (Gorkha) and M_w 7.3 (Dolakha) earthquakes in Nepal, through interpretation of very high resolution satellite images (e.g. WorldView, Pleiades, Cartosat-1 and 2, Resourcesat-2). Our spatial analysis of landslide occurrences with ground acceleration, slope, lithology and surface defomation indicated ubiquitous control of steep slope on landslides with ground acceleration as the trigger. Spatial distribution of landslides shows increasing frequency away from the Gorkha earthquake epicentre up to 130 km towards east, dropping sharply thereafter, which is an abnormal phenomenon of coseismic landslides. Landslides are laterally concentrated in three zones which matches well with the seismic rupture evolution of Gorkha earthquake, as reported through teleseismic measurements.     

Investigating a reservoir bank slope displacement history with multi-frequency satellite SAR data

Slope instability arisen along with dam construction is a common problem of great concern in reservoir areas. Thus, displacement monitoring of active slopes is of great importance for the safety of dam operation. The unstable Guobu slope is located only about 1.5 km away from Laxiwa hydropower station in upstream Yellow River. In this study, Synthetic Aperture Radar (SAR) datasets acquired by C-band Environmental Satellite (ENVISAT) Advanced Synthetic Aperture Radar (ASAR), L-band Advanced Land Observing Satellite 2 (ALOS-2) Phased Array type L-band Synthetic Aperture Radar 2 (PALSAR-2), and X-band TerraSAR-X covering different evolution stages of Guobu slope were collected to investigate the displacement history so as to facilitate understanding of its deformation and failure mechanisms. The displacements occurred during the past decade were quantitatively identified for the first time by SAR pixel offset tracking analyses. The results show that before the reservoir impoundment, the maximum accumulative displacements on the slope were more than 7 m from 2003 to 2008, while the post-impoundment displacements also exceeded 7 m in just 1 year from September 2009 to September 2010. Furthermore, this slope is still in active deformation up to now. Nevertheless, the displacement rates seem decreased recently according to the interferometric results of TerraSAR-X data pairs from September 2015 to March 2016.     

Particle shape effect on thermal conductivity and shear wave velocity in sands

This study presents the correlations between quantified shape parameters and geotechnical properties for nine sand specimens. Four shape parameters, sphericity, convexity, elongation and slenderness, were quantified with two-dimensional microscopic images with the aid of image processing techniques. An instrumented oedometer cell is used to measure compressibility, thermal conductivity and shear wave velocity during loading, unloading and reloading stages. As the particle shape inherently determines the initial loose packing condition, initial void ratio and shape parameters are well correlated with compressibility. The applied stress in soils increases the interparticle contact area and contact quality; round particles tend to achieve higher thermal conductivity and shear wave velocity during stress-induced volume change. Multiple linear regression is implemented to capture the relative contributions of involved variables, revealing that the thermal evolution is governed by the initial packing density and particle shape. The experimental observations underscore the predominant effect that particle shape has on the geomechanical and physical properties upon stress-induced soil behavior.     

Pile reinforcement mechanism of soil slopes

Stabilizing piles are widely used as an effective and economic reinforcement approach for slopes. Reasonable designs of pile reinforcement depend on the understanding of reinforcement mechanism of slopes. A series of centrifuge model tests were conducted on the pile-reinforced slopes and corresponding unreinforced slopes under self-weight and vertical loading conditions. The deformation of the slope was measured using image-based analysis and employed to investigate the pile reinforcement mechanism. The test results showed that the piles significantly reduced the deformation and changed the deformation distribution of the slope, and prevented the failure occurred in the unreinforced slope. The pile influence zone was determined according to the inflection points on the distribution curves of horizontal displacement, which comprehensively described the features of the pile–slope interaction and the characteristics of reinforced slopes. The concepts of anti-shear effect and compression effect were proposed to quantitatively describe the restriction features of the piles on the deformation of the slope, namely the reduction in the shear deformation and the increase in the compression deformation, respectively. The pile reinforcement effect mainly occurred in the pile influence zone and decreased with increasing distance from the piles. There was a dominated compression effect in the vicinities of the piles. The compression effect developed upwards in the slope with a transmission to the anti-shear effect. The anti-shear effect became significantly dominated near the slip surface and prevented the failure that occurred in the unreinforced slope.     

The formation of porphyry copper deposits

Copper is a moderately incompatible chalcophile element. Its behavior is strongly controlled by sulfides. The speciation of sulfur is controlled by oxygen fugacity. Therefore, porphyry Cu deposits are usually oxidized (with oxygen fugacities > ΔFMQ +2) (Mungall 2002; Sun et al. 2015). The problem is that while most of the magmas at convergent margins are highly oxidized, porphyry Cu deposits are very rare, suggesting that high oxygen fugacity alone is not sufficient. Partial melting of mantle peridotite even at very high oxygen fugacities forms arc magmas with initial Cu contents too low to form porphyry Cu deposits directly (Lee et al. 2012; Wilkinson 2013). Here we show that partial melting of subducted young oceanic slabs at high oxygen fugacity (>ΔFMQ +2) may form magmas with initial Cu contents up to >500 ppm, favorable for porphyry mineralization. Pre-enrichment of Cu through sulfide saturation and accumulation is not necessarily beneficial to porphyry Cu mineralization. In contrast, re-melting of porphyritic hydrothermal sulfide associated with iron oxides may have major contributions to porphyry deposits. Thick overriding continental crust reduces the “leakage” of hydrothermal fluids, thereby promoting porphyry mineralization. Nevertheless, it is also more difficult for ore forming fluids to penetrate the thick continental crust to reach the depths of 2–4 km where porphyry deposits form.     

Effects of river damming on biogenic silica turnover: implications for biogeochemical carbon and nutrient cycles

Rivers link terrestrial ecosystems and marine ecosystems, and they transport large amounts of substances into oceans each year, including several forms of silicon (Si), carbon (C), and other nutrients. However, river damming affects the water flow and biogeochemical cycles of Si, C, and other nutrients through biogeochemical interacting processes. In this review, we first summarize the current understanding of the effects of river damming on the processes of biogeochemical Si cycle, especially the source, composition, and recycling process of biogenic silica (BSi). Then, we introduce dam impacts on the cycles of C and some other nutrients. Dissolved silicon in rivers is mainly released from phytolith dissolution and silicate weathering. BSi in suspended matter or sediments in most rivers mainly consists of phytoliths and mainly originates from soil erosion. However, diatom growth and deposition in many reservoirs formed by river interception may significantly increase the contribution of diatom Si to total BSi, and thus significantly influence the biogeochemical Si, C, and nutrient cycles. Yet the turnover of phytoliths and diatoms in different rivers formed by river damming is still poorly quantified. Thus, they should be further investigated to enhance our understanding about the effects of river damming on global biogeochemical Si, C and nutrient cycles.     

Nitrous oxide (N<Subscript>2</Subscript>O) emissions from a mesotrophic reservoir on the Wujiang River,southwest China

Aquatic ecosystems have been identified as a globally significant source of nitrous oxide (N₂O) due to continuous active nitrogen involvement, but the processes and influencing factors that control N₂O production are still poorly understood, especially in reservoirs. For that, monthly N₂O variations were monitored in Dongfeng reservoir (DFR) with a mesotrophic condition. The dissolved N₂O concentration in DFR displayed a distinct spatial–temporal pattern but lower than that in the eutrophic reservoirs. During the whole sampling year, N₂O saturation ranging from 144% to 640%, indicating that reservoir acted as source of atmospheric N₂O. N₂O production is induced by the introduction of nitrogen (NO₃ ^?, NH₄ ⁺) in mesotrophic reservoirs, and is also affected by oxygen level and water temperature. Nitrification was the predominate process for N₂O production in DFR due to well-oxygenated longitudinal water layers. Mean values of estimated N₂O flux from the air–water interface averaged 0.19 µmol m^?2 h^?1 with a range of 0.01–0.61 µmol m^?2 h^?1. DFR exhibited less N₂O emission flux than that reported in a nearby eutrophic reservoir, but still acted as a moderate N₂O source compared with other reservoirs and lakes worldwide. Annual emissions from the water–air interface of DFR were estimated to be 0.32 × 10⁵ mol N–N₂O, while N₂O degassing from releasing water behind the dam during power generation was nearly five times greater. Hence, N₂O degassing behind the dam should be taken into account for estimation of N₂O emissions from artificial reservoirs, an omission that historically has probably resulted in underestimates. IPCC methodology should consider more specifically N₂O emission estimation in aquatic ecosystems, especially in reservoirs, the default EF5 model will lead to an overestimation.