Variations in glacier volume and snow cover and their impact on lake storage in the Paiku Co Basin,in the Central Himalayas

Glaciers and snow cover are important constituents of the surface of the Tibetan Plateau. The responses of these phenomena to global environmental changes are sensitive, rapid and intensive due to the high altitudes and arid cold climate of the Tibetan Plateau. Based on multisource remote sensing data, including Landsat images, MOD10A2 snow product, ICESat, Cryosat-2 altimetry data and long-term ground climate observations, we analysed the dynamic changes of glaciers, snow melting and lake in the Paiku Co basin using extraction methods for glaciers and lake, the degree-day model and the ice and lake volume method. The interaction among the climate, ice-snow and the hydrological elements in Paiku Co is revealed. From 2000 to 2018, the basin tended to be drier, and rainfall decreased at a rate of −3.07 mm/a. The seasonal temperature difference in the basin increased, the maximum temperature increased at a rate of 0.02°C/a and the minimum temperature decreased at a rate of −0.06°C/a, which accelerated the melting from glaciers and snow at rates of 0.55 × 10⁷ m³/a and 0.29 × 10⁷ m³/a, respectively. The rate of contribution to the lake from rainfall, snow and glacier melted water was 55.6, 27.7 and 16.7%, respectively. In the past 18 years, the warmer and drier climate has caused the lake to shrink. The water level of the lake continued to decline at a rate of −0.02 m/a, and the lake water volume decreased by 4.85 × 10⁸ m³ at a rate of −0.27 × 10⁸ m³/a from 2000 to 2018. This evaluation is important for understanding how the snow and ice melting in the central Himalayas affect the regional water cycle.     

A novel coupled FDM-DEM modelling method for flexible membrane boundary in laboratory tests

Difficulties are involved in discrete element method (DEM) modelling of the flexible boundary, that is, the membranes covering the soil sample, which can be commonly found in contemporary laboratory soil tests. In this paper, a novel method is proposed wherein the finite difference method (FDM) and DEM are coupled to simulate the rubber membrane and soil body, respectively. Numerical plane strain and triaxial tests, served by the flexible membrane, are implemented and analysed later. The effect of the membrane modulus on the measurement accuracy is considered, with analytical formulae derived to judge the significance of this effect. Based on an analysis of stress-strain responses and the grain rotation field, the mechanical performances produced by the flexible and rigid lateral boundaries are compared for the plane strain test. The results show that (1) the effect of the membrane on the test result becomes more significant at larger strain level because the membrane applies additional lateral confining pressure to the soil body; (2) the tested models reproduce typical stress and volumetric paths for specimens with shear bands; (3) for the plane strain test, the rigid lateral boundary derives a much higher peak strength and larger bulk dilatation, but a similar residual strength, compared with the flexible boundary. The latter produces a more uniform (or ‘diffuse') rotation field and more mobilised local kinematics than does the former. All simulations show that the proposed FDM-DEM coupling method is able to simulate laboratory tests with a flexible boundary membrane.     

Study of the flood control scheduling scheme for the Three Gorges Reservoir in a catastrophic flood

The Three Gorges Project is the world's largest water conservancy project. According to the design standards for the 1,000‐year flood, flood diversion areas in the Jingjiang reach of the Yangtze River must be utilized to ensure the safety of the Jingjiang area and the city of Wuhan. However, once these areas are used, the economic and life loss in these areas may be very great. Therefore, it is vital to reduce this loss by developing a scheme that reduces the use of the flood diversion areas through flood regulation by the Three Gorges Reservoir (TGR), under the premise of ensuring the safety of the Three Gorges Dam. For a 1,000‐year flood on the basis of a highly destructive flood in 1954, this paper evaluates scheduling schemes in which flood diversion areas are or are not used. The schemes are simulated based on 2.5‐m resolution reservoir topography and an optimized model of dynamic capacity flood regulation. The simulation results show the following. (a) In accord with the normal flood‐control regulation discharge, the maximum water level above the dam should be not more than 175 m, which ensures the safety of the dam and reservoir area. However, it is necessary to utilize the flood diversion areas within the Jingjiang area, and flood discharge can reach 2.81 billion m³. (b) In the case of relying on the TGR to impound floodwaters independently rather than using the flood diversion areas, the maximum water level above the dam reaches 177.35 m, which is less than the flood check level of 180.4 m to ensure the safety of the Three Gorges Dam. The average increase of the TGR water level in the Chongqing area is not more than 0.11 m, which indicates no significant effect on the upstream reservoir area. Comparing the various scheduling schemes, when the flood diversion areas are not used, it is believed that the TGR can execute safe flood control for a 1,000‐year flood, thereby greatly reducing flood damage.     

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.     

Dynamic modelling of retrogressive landslides with emphasis on the role of clay sensitivity

This paper presents a detailed numerical study of the retrogressive failure of landslides in sensitive clays. The dynamic modelling of the landslides is carried out using a novel continuum approach, the particle finite element method, complemented with an elastoviscoplastic constitutive model. The multiwedge failure mode in the collapse is captured successfully, and the multiple retrogressive failures that have been widely observed in landslides in sensitive clays are reproduced with the failure mechanism, the kinematics, and the deposition being discussed in detail. Special attention has been paid to the role of the clay sensitivity on each retrogressive failure as well as on the final retrogression distance and the final run‐out distance via parametric studies. Moreover, the effects of the viscosity of sensitive clays on the failure are also investigated for different clay sensitivities.     

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.     

Active capture and anomalous adsorption: New mechanisms for the incorporation of heavy noble gases

Abstract— Active capture is a new process for the incorporation of large quantities of heavy noble gases into growing surfaces. Adsorption in the conventional sense involves surface bonding by polarization (Van der Waals forces). What is referred to as “anomalous adsorption” of heavy noble gases involves chemical bonds and can occur when other (more chemically active) species are not available to preempt sites with unfilled bonds. Anomalous adsorption has been observed under conditions of fracture, vacuum deposition and ionizing radiation. Active capture depends upon anomalous adsorption to retain noble gases on a surface long enough to be captured in a growing surface film as it is deposited. The fundamental principle may be the impingement onto the growing film with sufficient energy to liberate surface electrons (work function energy of a few electronvolts) so that they are retained by anomalous adsorption long enough to be entrapped in the growing surface. Trapping efficiencies of ?1% have been observed for Kr and Xe in laboratory experiments, implying a fundamentally new mechanism for the incorporation of heavy noble gases onto surfaces. It may play a role in explaining the large concentrations of planetary noble gases contained in phase‐Q.     

Gamma-Hadron Separation using Čerenkov Photon Density Fluctuations

In the atmospheric Čerenkov technique γ-rays are detected against the abundant background produced by hadronic showers. In order to improve the signal to noise ratio of theexperiment, it is necessary to reject a significant fraction of hadronic showers. Traditional background rejection methods based on image shape parameters have been extensively used for the data from imaging telescopes. However, non-imaging Čerenkov telescopes have to develop very different means of statistically identifying and removing cosmic ray events. Some of the parameters, which could be potentially important for non-imaging arrays, are the temporal and spectral differences, the lateral distributions and density fluctuations of Čerenkov photons generated by γ-ray and hadron primaries. Here we study the differences in fluctuations of Čerenkov photon density in the light pool at the observation level from showers initiated by photons and those initiated by protons or heavier nuclei. The database of simulated events for the PACT array has been used to evaluate the efficiency of the new technique. Various types of density fluctuations like the short range and medium range fluctuations as well as flatness parameter are studied. The estimated quality factors reflect the efficiencies with which the hadrons can be rejected from the data. Since some of these parameters are independent, the cuts may be applied in tandem and we demonstrate that the proton rejection efficiency of ∼90% can be achieved. Use of density fluctuations is particularly suited for wavefront sampling observations and it seems to be a good technique to improve the signal to noise ratio. This revised version was published online in July 2006 with corrections to the Cover Date.