Characteristics of deuterium excess parameters for geothermal water in Beijing

This study investigates the characteristics of geothermal water in 10 geothermal fields in Beijing. The relationships between the deuterium excess parameter (d) and temperature, depth, age of geothermal groundwater, groundwater flow field, and Eh were investigated using geothermal groundwater samples. Results showed that (1) the average d value of geothermal water is 5.4, whereas that of the groundwater in normal temperature is 6.04. The differences are induced by the oxygen isotope exchange during the water–rock interaction, which may be more easily completed in geothermal water than in cold groundwater. (2) The d value increases remarkably with the age of the geothermal groundwater. The d value increases from 11.2 to 14.6 when the age of the geothermal water is 12,760 ± 130 a and 38,960 ± 630 a, respectively. Moreover, the isotope heat exchange for composition of the hydrogen and oxygen isotopes in the geothermal groundwater proceeds sufficiently with time. (3) The d value decreases from 5.72 to 3.03 when the depth increases from 125.13 to 3221 m. Generally, in the same area, the d value decreases with depth because the temperature is increasing. (4) The d value of the groundwater gradually reduces from the northern recharge area to the southern discharge area. The average d value is 7.31 in the northern recharge area and 5.68 in the middle Beijing Depression, whereas the d value in the southern area of Fengheying is ?9.20. The larger difference in d values between the recharge and discharge areas is due to the slower velocity of underwater flow, which induces longer time for oxygen exchange. (5) The relationship between the d and Eh is complex. When Eh is <200 mV, the d value of the geothermal water decreases with the decrease in Eh. When Eh is higher than 200 mV, the d value increases slightly with the decrease in Eh. The study of the characteristics of deuterium excess parameters for geothermal water could provide a scientific isotopic evidence for assessment and exploitation measures in geothermal groundwater systems.     

Accuracy assessment of near-shore bathymetry information retrieved from Landsat-8 imagery

The improvement in the capabilities of Landsat-8 imagery to retrieve bathymetric information in shallow coastal waters was examined. Landsat-8 images have an additional band named coastal/aerosol, Band 1: 435–451 nm in comparison with former generation of Landsat imagery. The selected Landsat-8 operational land image (OLI) was of Chabahar Bay, located in the southern part of Iran (acquired on February 22, 2014 in calm weather and relatively low turbidity). Accurate and high resolution bathymetric data from the study area, produced by field surveys using a single beam echo-sounder, were selected for calibrating the models and validating the results. Three methods, including traditional linear and ratio transform techniques, as well as a novel proposed integrated method, were used to determine depth values. All possible combinations of the three bands [coastal/aerosol (CB), blue (B), and green (G)] have been considered (11 options) using the traditional linear and ratio transform techniques, together with five model options for the integrated method. The accuracy of each model was assessed by comparing the determined bathymetric information with field measured values. The standard error of the estimates, correlation coefficients (R ² ) for both calibration and validation points, and root mean square errors (RMSE) were calculated for all cases. When compared with the ratio transform method, the method employing linear transformation with a combination of CB, B, and G bands yielded more accurate results (standard error = 1.712 m, R ² _calibration = 0.594, R ² _validation = 0.551, and RMSE =1.80 m). Adding the CB band to the ratio transform methodology also dramatically increased the accuracy of the estimated depths, whereas this increment was not statistically significant when using the linear transform methodology. The integrated transform method in form of Depth = b ₀  + b ₁ X _CB  + b ₂ X _B  + b ₅ ln(R _CB )/ln(R _G ) + b ₆ ln(R _B )/ln(R _G ) yielded the highest accuracy (standard error = 1.634 m, R ² _calibration = 0.634, R ² _validation = 0.595, and RMSE = 1.71 m), where R _i (i = CB, B, or G) refers to atmospherically corrected reflectance values in the i ^th band [X _i  = ln(R _i -R _{deep water})].     

Modeling and mapping of solar radiation using geostatistical analysis methods in Iran

Because of economic and technical limitations, measuring solar energy received at ground level (R _s ) isn’t possible in all parts of the country, and in only 12% of synoptic stations is this parameter measured and recorded. Thus, it should be estimated and modeled spatially based on other climatic variables using mathematical methods. In this research, many attempts have been made to introduce an air temperature-based model for Rs estimation, and then, based on the output of the mentioned models, several geostatistical methods have been tested, and finally an elegant spatial model is proposed for (Rs) zoning in Iran. In this regard, the relationships between the measured amounts of monthly solar radiation and other climatic parameters, such as a monthly average, maximum and minimum temperature, precipitation, relative humidity, and the number of sunny hours during the period 1970–2010, are examined and modeled. It was revealed that based on the linear relationship between the monthly average air temperatures and solar radiation values recorded in each of the stations, that the best-fit linear model, with R ²  = 0.822, MAE = 1.81, RMSE = 2.51%, and MAPE = 10.08, can be introduced for Rs estimation. Then, using the outputs of the proposed model, the amounts of (R _s ) are estimated in another 171 meteorological stations (a total of 192 stations), and eight geostatistical methods (IDW, GPI, RBF, LPI, OK, SK, UK, and EBK) were investigated for zoning. Comparing the resulting variograms showed that in addition to proof of spatial correlation between solar radiation data, they can be applied for modeling changes in various directions. Analyzing the ratio of the nugget effect on the roof of the variograms showed that the Gaussian model with the lowest ratio (Co/Co + C = 0.883) and (R ²  = 0.972), could model the highest correlation between the data and, therefore, it was used for data interpolation. To select the best geostatistical model, R2, MAE, and RMSE were used. On this basis, it was found that the RBF method with R ²  = 0.904, MAE = 3.02, RMSE = 0.39% is the most effective. Also, the IDW method with R ²  = 0.90, MAE = 3.08, RMSE = 0.391%, compared to other methods is the most effective. In addition, for data validation, correlations between observed and estimated values of solar radiation were studied and found R ²  = 0.86.     

Mapping heatwave vulnerability in Korea

Analysis of event-based soil erosion magnitude with special return periods is essential to appropriately design strategies and adopt soil conservation practices. However, the spatiotemporal variations of soil erosion with different return periods, especially at national level, have not been adequately considered. Therefore, the present study aimed to zone rainfall erosivity index (R factor) as the most dynamic factor affecting variability of soil erosion rate, with different return periods in monthly, seasonal and annual time scales in Iran. Toward this attempt, the kinetic energy and maximum 30-min intensity (I ₃₀) over 12,000 available and accessible events of 70 stations were calculated during the common period of 1984–2004 and the corresponding R factor of the Universal Soil Loss Equation was then computed. Subsequently, the best-fitted frequency distributions were determined in all stations in three time scales using the EasyFit Software. The R factor was accordingly estimated for 2-, 5-, 10-, 25- and 50-year return periods. In addition, the inverse distance weighting technique was employed to determine and analyze the spatial variability patterns of R factor in different time scales using geographic information system. The results indicated that the frequency distributions fitted to study data were different in study time scales due to variability of spatiotemporal patterns of R factor. In addition, no specific spatial pattern of R factor could be recognized for different return periods and time scales. The average annual R factor was also found 1.41 MJ mm ha^?1 h^?1, whereas the respective R factor for different respective return periods of 2, 5, 10, 25 and 50 years was obtained 1.47, 2.62, 3.35, 4.48 and 5.54 MJ mm ha^?1 h^?1. These findings can be used for suitable decision making and effective environmental planning for land management Iran countrywide.     

Permeability of granular soil employing flexible wall permeameter

Understanding the changes in permeability of soil, when soil is subjected to high confining pressure and flow pressure, which may alter the textural and geomechanical characteristics of soil, is of great importance to many geo-engineering activities such as, construction of high-rise buildings near the coast or the water bodies, earthen dams, pavement subgrades, reservoir, and shallow repositories. It is now possible to evaluate the changes in permeability of soil samples under varying conditions of confining pressure and flow pressure using flexible wall permeameter (FWP). In the present study, investigation was carried out on a cylindrical sample of granular soil employing FWP under varied conditions of confining pressure (σ₃)—50–300 kPa, which can simulate the stress conditions equivalent to depth of about 20 m under the earth’s crust, and a flow pressure (f_p)—20–120 kPa, which is mainly present near the small earthen embankment dams, landfill liners, and slurry walls near the soft granular soil with high groundwater table. The obtained results indicate a linear relationship between hydraulic conductivity (k) with effective confining pressure (σ_eff.), k, decreasing linearly with an incremental change in σ_eff.. Further, k increases significantly with an increase in f_p corresponding to each σ_eff., and q increases significantly with increase in the f_p corresponding to each (σ₃). It was also observed that corresponding to the low f_p of 20 kPa, the reduction in k is nonlinear with σ_3. The percentage reduction in k is observed to be 9, 13, and 27% corresponding to σ₃ of 50–100, 100–200, and 200-300 kPa, respectively.     

Modelling the effects of porous media deformation on the propagation of water-table waves in a sandy unconfined aquifer

This paper examines the influence of porous media deformation on water-table wave dispersion in an unconfined aquifer using a numerical model which couples Richards’ equation to the poro-elastic model. The study was motivated by the findings of Shoushtari et al. (J Hydrol 533:412–440, 2016) who were unable to reproduce the observed wave dispersion in their sand flume data with either numerical Richards’ equation models (assuming rigid porous media) or existing analytic solutions. The water-table wave dispersion is quantified via the complex wave number extracted from the predicted amplitude and phase profiles. A sensitivity analysis was performed to establish the influence of the main parameters in the poro-elastic model, namely Young’s modulus (E) and Poisson’s ratio (ν). For a short oscillation period (T?=?16.4 s), the phase lag increase rate (k _i) is sensitive to the chosen values of E and ν, demonstrating an inverse relationship with both parameters. Changes in the amplitude decay rate (k _r), however, were negligible. For a longer oscillation period (T?=?908.6 s), variations in the values of E and ν resulted in only small changes in both k _r and k _i. In both the short and long period cases, the poro-elastic model is unable to reproduce the observed wave dispersion in the existing laboratory data. Hence porous media deformation cannot explain the additional energy dissipation in the laboratory data. Shoushtari SMH, Cartwright N, Perrochet P, Nielsen P (2016) The effects of oscillation period on groundwater wave dispersion in a sandy unconfined aquifer: sand flume experiments and modelling. J Hydrol 533:412–440.     

Hydraulic Conductivity of Fly Ash-Amended Mine Tailings

The objective of this study was to evaluate the effect of fly ash addition on hydraulic conductivity (k) of mine tailings. Mine tailings used in this study were categorized as synthetic tailings and natural tailings; two synthetic tailings were developed via blending commercially-available soils and natural tailings were collected from a garnet mine located in the U.S. Two fly ashes were used that had sufficient calcium oxide (CaO) content (17 and 18.9 %) to generate pozzolanic activity. Hydraulic conductivity was measured on pure tailings and fly ash-amended tailings in flexible-wall permeameters. Fly ash was added to mine tailings to constitute 10 % dry mass of the mixture, and specimens were cured for 7 and 28 days. The influence of fly ash-amendment on k of mine tailings was attributed to (1) molding water content and (2) plasticity of the mine tailings. Tailings that classified as low-plasticity silts with clay contents less than 15 % exhibited a decrease in k when amended with fly ash and prepared wet of optimum water content (w _opt ). Tailings that classified as low-plasticity clay exhibited a one-order magnitude increase in k with addition of fly ash for materials prepared dry or near w _opt . The decrease in k for silty tailings was attributed to formation of cementitious bonds that obstructed flow paths, whereas the increase in k for clayey tailings was attributed to agglomeration of clay particles and an overall increase in average pore size. The results also indicated that the effect of curing time on k is more pronounced during the early stages of curing (≤7 days), as there was negligible difference between k for 7 and 28-days cured specimens.     

Synthesis,TEM characterization and thermal behaviour of LiNiSi<Subscript>2</Subscript>O<Subscript>6</Subscript> pyroxene

A pyroxene with composition LiNiSi₂O₆ was synthesized at T = 1,473 K and P = 2.0 GPa; the cell parameters at T = 298 K are a = 9.4169(6) Å, b = 8.4465(7) Å, c = 5.2464(3) Å, β = 110.534(6)°, V = 390.78(3) Å³. TEM examination of the LiNiSi₂O₆ pyroxene showed the presence of h + k odd reflections indicative of a primitive lattice, and of antiphase domains obtained by dark field imaging of the h + k odd reflections. A HT in situ investigation was performed by examining TEM selected area diffraction patterns collected at high temperature and synchrotron radiation powder diffraction. In HTTEM the LiNiSi₂O₆ was examined together with LiCrSi₂O₆ pyroxene. In LiCrSi₂O₆ the h + k odd critical reflections disappear at about 340 K; they are sharp up to the transition temperature and do not change their shape until they disappear. In LiNiSi₂O₆ the h + k odd reflections are present up to sample deterioration at 650 K. A high temperature synchrotron radiation powder diffraction investigation was performed on LiNiSi₂O₆ between 298 and 773 K. The analysis of critical reflections and of changes in cell parameters shows that the space group is P-centred up to the highest temperature. The comparative analysis of the thermal and spontaneous strain contributions in P2₁/c and C2/c pyroxenes indicates that the high temperature strain in P-LiNiSi₂O₆ is very similar to that due to thermal strain only in C2/c spodumene and that a spontaneous strain contribution related to pre-transition features is not apparent in LiNiSi₂O₆. A different high-temperature behaviour in LiNiSi₂O₆ with respect to other pyroxenes is suggested, possibly in relation with the presence of Jahn–Teller distortion of the M1 polyhedron centred by low-spin Ni³⁺.     

Nonlinear consolidation in randomly heterogeneous highly compressible aquitards

Severe land subsidence due to groundwater extraction may occur in multiaquifer systems where highly compressible aquitards are present. The highly compressible nature of the aquitards leads to nonlinear consolidation where the groundwater flow parameters are stress-dependent. The case is further complicated by the heterogeneity of the hydrogeologic and geotechnical properties of the aquitards. The effect of realistic vertical heterogeneity of hydrogeologic and geotechnical parameters on the consolidation of highly compressible aquitards is investigated by means of one-dimensional Monte Carlo numerical simulations where the lower boundary represents the effect of an instant drop in hydraulic head due to groundwater pumping. Two thousand realizations are generated for each of the following parameters: hydraulic conductivity (K), compression index (C _c), void ratio (e) and m (an empirical parameter relating hydraulic conductivity and void ratio). The correlation structure, the mean and the variance for each parameter were obtained from a literature review about field studies in the lacustrine sediments of Mexico City. The results indicate that among the parameters considered, random K has the largest effect on the ensemble average behavior of the system when compared to a nonlinear consolidation model with deterministic initial parameters. The deterministic solution underestimates the ensemble average of total settlement when initial K is random. In addition, random K leads to the largest variance (and therefore largest uncertainty) of total settlement, groundwater flux and time to reach steady-state conditions.     

Prediction of local scour around complex piers using GEP and M5-Tree

Local scour around piers is one of the main causes of bridge failures. In this study, three robust techniques, artificial neural networks (ANNs), M5-Tree, and Gene Expression Programming (GEP), were employed for prediction of scour depth around complex piers. The clear water condition was chosen for all experimental tests. The results indicated that pier diameter (b _c) and foundation level (Y) are the main parameters for local scour. Furthermore, the minimum scour depth occurs in range of Y/b _c = 1.1~1.3. In next step, to evaluate the mentioned techniques, a wide range of dataset was collected from the present study and literature. The radial base function (RBF) with R ² = 0.945 and RMSE = 0.031 provides better prediction in comparison with conventional equations, M5-Tree (R ² = 0.883, RMSE = 0.292) and the GEP techniques (R ² = 0.811 and RMSE = 0.263). The equations developed by M5-Tree and GEP are more useful for practical purposes and can be easily employed to predict the depth of scour at complex piers.     

Influences of activated sludge surface properties on adsorption of aqueous fullerene C<Subscript>60</Subscript> nanoparticles

Rapid increases in the amounts of fullerene C₆₀ nanoparticles (nC₆₀) being produced and used will inevitably lead to increases in the amounts released into the aquatic environment. This will have implications for human and ecosystem health. Wastewater treatment plants are key barriers to nC₆₀ being released into aquatic systems, but little information is available on how adsorption processes in wastewater treatment plants affect the fates of nC₆₀. We investigated the effects of the surface properties of activated sludge on the adsorption of nC₆₀ and related mechanisms by modeling the adsorption kinetics and equilibrium process and performing correlation analyses. The adsorption of nC₆₀ closely followed the pseudo-second-order kinetic model (R > 0.983), the Freundlich isotherm model (R > 0.990), and the linear partitioning isotherm model (R > 0.966). Different adsorption coefficients, 1.070–4.623 for the Freundlich partitioning model and 1.788–6.148 for the linear partitioning model, were found for different types of activated sludge. The adsorption coefficients significantly positively correlated with the zeta (ζ) potential (R = 0.877) and hydrophobicity (R = 0.661) and negatively correlated with particle size (R = ?0.750). The results show that nC₆₀ adsorption is strongly affected by the surface properties of activated sludge because changes in surface properties cause changes in the electrostatic and hydrophobic interactions that occur.     

Field monitoring of groundwater responses to heavy rainfalls and the early warning of the Kualiangzi landslide in Sichuan Basin,southwestern China

The Kualiangzi landslide was triggered by heavy rainfalls in the “red beds” area of Sichuan Basin in southwestern China. Differing from other bedrock landslides, the movement of the Kualiangzi landslide was controlled by the subvertical cracks and a subhorizontal bedding plane (dip angle < 10°). The ingress of rainwater in the cracks formed a unique groundwater environment in the slope. Field measurement for rainfall, groundwater movement, and slope displacement has been made for the Kualiangzi landslide since 2013. The field monitoring system consists of two rainfall gauges, seven piezometers, five water-level gauges, and two GPS data loggers. The equipments are embedded near a longitudinal section of the landslide, where severe deformation has been observed in the past 3 years. The groundwater responses to four heavy rainfall events were analyzed between June 16 and July 24 in 2013 coincided with the flood season in Sichuan. Results showed that both of the water level and the pore-water pressure increased after each rainfall event with delay in the response time with respect to the precipitation. The maximum time lag reached 35 h occurred in a heavy rainfall event with cumulative precipitation of 127 mm; such lag effect was significantly weakened in the subsequent heavy rainfall events. In each presented rainfall event, longer infiltration period in the bedrock in the upper slope increased the response time of groundwater, compared to that of in the gravels in the lower slope. A translational landslide conceptual model was built for the Kualiangzi landslide, and the time lag was attributed to the gradual formation of the uplift pressure on the slip surface and the softening of soils at the slip surface. Another important observation is the effect on the slope movement which was caused by the water level (H _w) in the transverse tension trough developed at the rear edge of the landslide. Significant negative correlation was found for H _w and the slope stability factor (F _s), in particular for the last two heavy rainfall events, of which the drastic increase of water level caused significant deterioration in the slope stability. The rapid drop (Δ?=?22.5 kPa) of pore-water pressure in the deep bedrock within 1 h and the large increase (Δ?=?87.3 mm) of surficial displacement were both monitored in the same period. In the end, a four-level early warning system is established through utilizing H _w and the displacement rate D _r as the warning indicators. When the large deformation occurred in flood season, the habitants at the leading edge of the landslide can be evacuated in time.     

High-pressure displacive phase transition of a natural Mg-rich pigeonite

A high-pressure single-crystal X-ray diffraction study has been carried out on a P2₁/c natural Mg-rich pigeonite sample with composition ca. Wo₆En₇₆Fs₁₈ using a diamond anvil-cell. The unit-cell parameters were determined at 14 different pressures to 7.14 GPa. The sudden disappearance of the b-type reflections (h + k = odd) and a strong discontinuity (about 2.8%) in the unit-cell volume indicated a first-order P2₁/c–C2/c phase transition between 4.66 and 4.88 GPa. The P(V) data of the P2₁/c phase were fitted to 4.66 GPa by a third-order Birch–Murnaghan equation of state (BM3 EoS), whereas the limited number of experimental data collected within the C2/c phase between 4.88 and 7.14 GPa were fitted using the same equation of state but with K′ constrained to the value obtained for the P2₁/c fitting. The equation of state coefficients are V ₀ = 424.66(6) ų, K _T0 = 104(2) GPa and K′ = 8(1) for the P2₁/c phase, and V ₀ = 423.6(1) ų, K _T0 = 112.4(8) GPa, and K′ fixed to 8(1) for the C2/c phase. The axial moduli for a, b, and c for the P2₁/c phase were obtained using also a BM3-EoS, while for the C2/c phase only a linear calculation could be performed, and therefore the same approach was applied for comparison also to the P2₁/c phase. In general the C2/c phase exhibits axial compressibilities (β _c > β _a >> β _b) lower than those of the P2₁/c phase (β _b > β _c ≈ β _a; similar to those found in previous studies in clinopyroxenes and orthopyroxenes). The lower compressibility of the C2/c phase compared with that of the P2₁/c could be ascribed to the greater stiffness along the b direction. A previously published relationship between P _c and M2 average cation radius (i.r.) has been updated using all the literature data on P2₁/c clinopyroxene containing large cations at M2 site and our new data. The following weighted regression was obtained: P _c (GPa) = 26(4) ? 28(5) ×  i.r (Å), R ² = 0.97. This improved equation can be used to predict the critical pressure of natural P2₁/c clinopyroxene samples just knowing the composition at M2 site.     

Quantifying land use influences on event-based flow frequency,timing, magnitude,and rate of change in an urbanizing watershed of the central USA

There is a need for research that advances understanding of flow alterations in contemporary watersheds where natural and anthropogenic interactions can confound mitigation efforts. Event-based flow frequency, timing, magnitude, and rate of change were quantified at five-site nested gauging sites in a representative mixed-land-use watershed of the central USA. Statistically independent storms were paired by site (n = 111 × 5 sites) to test for significant differences in event-based rainfall and flow response variables (n = 17) between gauging sites. Increased frequency of small peak flow events (i.e., 64 more events less than 4.0 m³ s^?1) was observed at the rural–urban interface of the watershed. Differences in flow response were apparent during drier periods when small rainfall events resulted in increased flow response at urban sites in the lower reaches. Relationships between rainfall and peak flow were stronger with decreased pasture/crop land use and increased urban land use by approximately 20%. Event-based total rainfall explained 40–68% of the variance in peak flow (p < 0.001). Coefficients of determination (r²) were negatively correlated with pasture/crop land use (r² = 0.92; p = 0.007; n = 5) and positively correlated with urban land use (r² = 0.90; p = 0.008; n = 5). Significant differences in flow metrics were observed between rural and urban sites (p < 0.05; n = 111) that were not explained by differences in rainfall variables and drainage area. An urban influence on flow timing was observed using median time lag to peak centroid and time of maximum precipitation to peak flow. Results highlight the need to establish manageable flow targets in rapidly urbanizing mixed-land-use watersheds.     

Rogue events in spatiotemporal numerical simulations of unidirectional waves in basins of different depth

The evolution of unidirectional nonlinear sea surface waves is calculated numerically by means of solution of the Euler equations. The wave dynamics corresponds to quasi-equilibrium states characterized by JONSWAP spectra. The spatiotemporal data are collected and processed providing information about the wave height probability and typical appearance of abnormally high waves (rogue waves). The waves are considered at different water depths ranging from deep to relatively shallow cases (k _p h > 0.8, where k _p is the peak wavenumber, and h is the local depth). The asymmetry between front and rear rogue wave slopes is identified; it becomes apparent for sufficiently high waves in rough sea states at all considered depths k _p h ≥ 1.2. The lifetimes of rogue events may reach up to 30–60 wave periods depending on the water depth. The maximum observed wave has a height of about three significant wave heights. A few randomly chosen in situ time series from the Baltic Sea are in agreement with the general picture of the numerical simulations.     

The morphological characteristics of gully systems and watersheds in Dry-Hot Valley,SW China

Gully systems and watersheds are geomorphic units with clear boundaries that are relatively independent of basin landscapes and play an important role in natural geography. In order to explore the morphological characteristics of gully systems and watersheds in the Dry-Hot Valley [South West (SW) China], gullies are interpreted from online Google images with high resolution and watersheds are extracted from digital elevation model at a scale of 1:50,000. The results show that: (1) There are 17,382 gullies (with a total area of 1141.66 km²) and 42 watersheds in the study area. (2) The average gully density of the study area (D) is 4.29 km/km², gully frequency (F) is 14.39 gullies/km², the branching ratio (B) is 5.13, the length ratio (L) is 3.12, and the coefficient of the main and tributary gullies (M) is 0.06. The degree of gully erosion is strong to extremely strong, the main development intensity of gully erosion ranges from intense to moderate, and the type of gully system is tributary. (3) The watershed areas (A) are between 0.39 and 96.43 km², the relief ratio (R) is from 0.10 to 0.19, the circularity ratio (C) is from 0.30 to 0.83, the texture ratio (T) is from 0.82 to 39.35, and the dominant geomorphological texture type is fine. (4) There is a quantitative relationship between F and D:F?=?0.624D² (R?=0.84) and T is closely related to D, F, M (R²?>?0.7). A, R and C are related to M (R²?>?0.5). The development of gully systems is the result of coupling effects between multiple factors. In this area, the degree of erosion and the condition of the main and tributary gullies can be controlled by the degree of topographic breakage in the watershed, which provides some theoretical basis for the evaluation of gully erosion by the latter. In addition, the scale, relief, and shape have a significant impact on the locations of the main and tributary gullies. For tributary gullies, attention should be paid to the interception and control of runoff and sediment in the small confluence branches in order to prevent gully expansion and head advance. These features can inform the development of targeted measures for the control of soil erosion.     

Tailings composition effects on shear strength behavior of co-mixed mine waste rock and tailings

The objective of this study was to evaluate the effect of mine tailings composition on shear behavior and shear strength of co-mixed mine waste rock and tailings (WR&T). Crushed gravel was used as a synthetic waste rock and mixed with four types of tailings: (1) fine-grained garnet, (2) coarse-grained garnet, (3) copper, and (4) soda ash. Co-mixed WR&T specimens were prepared to target mixture ratios of mass of waste rock to mass of tailings (R) such that tailings “just filled” interparticle void space of the waste rock (i.e., optimum mixture ratio, R _opt). Triaxial compression tests were conducted on waste rock, tailings, and mixed waste at effective confining stresses (\(\sigma_{\text{c}}^{{\prime }}\)) ranging from 5 to 40 kPa to represent stresses anticipated in final earthen covers for waste containment facilities. Waste rock and co-mixed WR&T specimens were 150 mm in diameter by 300 mm tall, whereas tailings specimens were 38 mm in diameter by 76 mm tall. Shear strength was quantified using effective stress friction angles (?′) from undrained tests: ?′ for waste rock was 37°, ?′ for tailings ranged from 34° to 41°, and ?′ for WR&T mixtures ranged from 38° to 40°. Thus, shear strength of co-mixed WR&T was comparable to waste rock regardless of tailings composition. Shear behavior of WR&T mixtures was a function of R and tailings composition. Tailings influenced shear behavior for R < R _opt and when tailings predominantly were silt. Shear behavior was influenced by waste rock for R ≥ R _opt and when tailings predominantly were sand or included clay particles.     

Evaluating the susceptibility of landslide landforms in Japan using slope stability analysis: a case study of the 2016 Kumamoto earthquake

This study analyzed 267 landslide landforms (LLs) in the Kumamoto area of Japan from the database of about 0.4 million LLs for the whole of Japan identified from aerial photos by the National Research Institute for Earth Science and Disaster Resilience of Japan (NIED). Each LL in the inventory is composed of a scarp and a moving mass. Since landslides are prone to reactivation, it is important to evaluate the sliding-recurrence susceptibility of LLs. One possible approach to evaluate the susceptibility of LLs is slope stability analysis. A previous study found a good correlation (R ² = 0.99) between the safety factor (F _s ) and slope angle (α) of F _s  = 17.3α ^?0.843. We applied the equation to the analysis of F _s for 267 LLs in the area affected by the 2016 Kumamoto earthquake (M _j  = 7.3). The F _s was calculated for the following three cases of failure: scarps only, moving mass only, and scarps and moving mass together. Verification with the 2016 Kumamoto earthquake event shows that the most appropriate method for the evaluation of LLs is to consider the failure of scarps and moving mass together. In addition, by analyzing the relationship between the factors of slope of entire landslide and slope of scarp for LLs and comparing the results with the Aso-ohashi landslide, the largest landslide caused by the 2016 Kumamoto earthquake, we also found that morphometric analysis of LLs is useful for forecasting the travel distance of future landslides.     

Characterization of locally available soil as a liner material for solid waste landfills in Sri Lanka

This study aimed to develop a low-cost and effective clay liner material for solid waste landfills in Sri Lanka. A locally available clayey soil and its admixtures with 5 and 10% bentonite were examined for this purpose. Laboratory experiments to determine soil plasticity and swell index were carried out on the tested samples. Hydraulic conductivity (k) tests were carried out in the laboratory using water and an aqueous solution of CaCl₂ on unconsolidated samples prepared by either dry or slurry packing and pre-consolidated samples with five different consolidation pressures (p) from 10 to 200 kPa. Measured liquid limits for tested admixtures increased with increasing bentonite contents and correlated well with measured values of the swell index. The difference in permeant solutions had little effect on measured k values for both unconsolidated and pre-consolidated samples. The hydraulic conductivities were highly affected by changing p, i.e., the k values decreased on two orders of magnitude as p increased from 10 to 200 kPa. The Kozeny–Carman equation, a theoretical permeability model that expresses the k-porosity relationship, was applied to measured data including reported values. Results showed the Kozeny–Carman equation captured well the porosity-dependent k values for tested soils and their admixtures with bentonite under a wide range of void ratios, suggesting that the Kozeny–Carman equation is a useful tool to estimate the magnitude of k values for differently compacted soil and its bentonite admixtures.