Granite is a typical crystalline material, often used as a building material, but also a candidate host rock for the repository of high-level radioactive waste. The petrographic texture—including mineral constituents, grain shape, size, and distribution—controls the fracture initiation, propagation, and coalescence within granitic rocks. In this paper, experimental laboratory tests and numerical simulations of a grain-based approach in two-dimensional Particle Flow Code (PFC2D) were conducted on the mechanical strength and failure behavior of Alashan granite, in which the grain-like structure of granitic rock was considered. The microparameters for simulating Alashan granite were calibrated based on real laboratory strength values and strain-stress curves. The unconfined uniaxial compressive test and Brazilian indirect tensile test were performed using a grain-based approach to examine and discuss the influence of mineral grain size and distribution on the strength and patterns of microcracks in granitic rocks. The results show it is possible to reproduce the uniaxial compressive strength (UCS) and uniaxial tensile strength (UTS) of Alashan granite using the grain-based approach in PFC2D, and the average mineral size has a positive relationship with the UCS and UTS. During the modeling, most of the generated microcracks were tensile cracks. Moreover, the ratio of the different types of generated microcracks is related to the average grain size. When the average grain size in numerical models is increased, the ratio of the number of intragrain tensile cracks to the number of intergrain tensile cracks increases, and the UCS of rock samples also increases with this ratio. However, the variation in grain size distribution does not have a significant influence on the likelihood of generated microcracks.

Hydraulic fracturing is a useful tool for enhancing rock mass permeability for shale gas development, enhanced geothermal systems, and geological carbon sequestration by the high-pressure injection of a fracturing fluid into tight reservoir rocks. Although significant advances have been made in hydraulic fracturing theory, experiments, and numerical modeling, when it comes to the complexity of geological conditions knowledge is still limited. Mechanisms of fluid injection-induced fracture initiation and propagation should be better understood to take full advantage of hydraulic fracturing. This paper presents the development and application of discrete particle modeling based on two-dimensional particle flow code (PFC^2D). Firstly, it is shown that the modeled value of the breakdown pressure for the hydraulic fracturing process is approximately equal to analytically calculated values under varied in situ stress conditions. Furthermore, a series of simulations for hydraulic fracturing in competent rock was performed to examine the influence of the in situ stress ratio, fluid injection rate, and fluid viscosity on the borehole pressure history, the geometry of hydraulic fractures, and the pore-pressure field, respectively. It was found that the hydraulic fractures in an isotropic medium always propagate parallel to the orientation of the maximum principal stress. When a high fluid injection rate is used, higher breakdown pressure is needed for fracture propagation and complex geometries of fractures can develop. When a low viscosity fluid is used, fluid can more easily penetrate from the borehole into the surrounding rock, which causes a reduction of the effective stress and leads to a lower breakdown pressure. Moreover, the geometry of the fractures is not particularly sensitive to the fluid viscosity in the approximate isotropic model.

This paper reviews the classical and some particular factors contributing to earthquake-triggered landslide activity. This analysis should help predict more accurately landslide event sizes, both in terms of potential numbers and affected area. It also highlights that some occurrences, especially those very far from the hypocentre/activated fault, cannot be predicted by state-of-the-art methods. Particular attention will be paid to the effects of deep focal earthquakes in Central Asia and to other extremely distant landslide activations in other regions of the world (e.g. Saguenay earthquake 1988, Canada). The classification of seismically induced landslides and the related ‘event sizes’ is based on five main factors: ‘Intensity’, ‘Fault factor’, ‘Topographic energy’, ‘Climatic background conditions’, ‘Lithological factor’. Most of these data were extracted from papers, but topographic inputs were checked by analyzing the affected region in Google Earth. The combination and relative weight of the factors was tested through comparison with well documented events and complemented by our studies of earthquake-triggered landslides in Central Asia. The highest relative weight (6) was attributed to the ‘Fault factor’; the other factors all received a smaller relative weight (2–4). The high weight of the ‘Fault factor’ (based on the location in/outside the mountain range, the fault type and length) is strongly constrained by the importance of the Wenchuan earthquake that, for example, triggered far more landslides in 2008 than the Nepal earthquake in 2015: the main difference is that the fault activated by the Wenchuan earthquake created an extensive surface rupture within the Longmenshan Range marked by a very high topographic energy while the one activated by the Nepal earthquake ruptured the surface in the frontal part of the Himalayas where the slopes are less steep and high. Finally, the calibrated factor combination was applied to almost 100 other earthquake events for which some landslide information was available. This comparison revealed the ability of the classification to provide a reasonable estimate of the number of triggered landslides and of the size of the affected area. According to this prediction, the most severe earthquake-triggered landslide event of the last one hundred years would actually be the Wenchuan earthquake in 2008 followed by the 1950 Assam earthquake in India – considering that the dominating role of the Wenchuan earthquake data (including the availability of a complete landslide inventory) for the weighting of the factors strongly influences and may even bias this result. The strongest landslide impacts on human life in recent history were caused by the Haiyuan-Gansu earthquake in 1920 – ranked as third most severe event according to our classification: its size is due to a combination of high shaking intensity, an important ‘Fault factor’ and the extreme susceptibility of the regional loess cover to slope failure, while the surface morphology of the affected area is much smoother than the one affected by the Wenchuan 2008 or the Nepal 2015 earthquakes. The main goal of the classification of earthquake-triggered landslide events is to help improve total seismic hazard assessment over short and longer terms. Considering the general performance of the classification-prediction, it can be seen that the prediction either fits or overestimates the known/observed number of triggered landslides for a series of earthquakes, while it often underestimates the size of the affected area. For several events (especially the older ones), the overestimation of the number of landslides can be partly explained by the incompleteness of the published catalogues. The underestimation of the extension of the area, however, is real – as some particularities cannot be taken into account by such a general approach: notably, we used the same seismic intensity attenuation for all events, while attenuation laws are dependent on regional tectonic and geological conditions. In this regard, it is likely that the far-distant triggering of landslides, e.g., by the 1988 Saguenay earthquake (and the related extreme extension of affected area) is due to a very low attenuation of seismic energy within the North American plate. Far-distant triggering of landslides in Central Asia can be explained by the susceptibility of slopes covered by thick soft soils to failure under the effect of low-frequency shaking induced by distant earthquakes, especially by the deep focal earthquakes in the Pamir – Hindukush seismic region. Such deep focal and high magnitude (> > 7) earthquakes are also found in Europe, first of all in the Vrancea region (Romania). For this area as well as for the South Tien Shan we computed possible landslide event sizes related to some future earthquake scenarios.