From Climate Conditions to the Numerical Slope Stability Analysis of Surface Coal Mines

Featured Application: This work has application in coal and lignite mining excavations and areas, especially during decarbonization and towards a post-coal and post-lignite era. Areas that hosted coal and lignite mines are now being reclaimed—often by installing renewable energy systems— and targets are being set to fulﬁll new energy principles in relation to energy production, harvesting, and storage. The appropriate and sustainable exploitation of these areas is directly related to their safety and stability. Abstract: A new perspective is presented for evaluating the slope stability of coal and lignite mines due to rainfall. The case of Greek lignite mining areas is employed to illustrate the methodology. Initially, past climatic records of rainfall are documented and analyzed; rainfall intensity varies from light (0.8 mm/h) to heavy (up to 9 mm/h). Few extreme phenomena are documented, with the maximum intensity being 17 mm/h. Furthermore, climatic projections of future trends are performed with open-access tools to anticipate possible deviations from the baseline conditions. Although the mean temperature is expected to increase, projections show that the past rainfall range is not expected to change. Finally, the effect of rainfall inﬁltration on the stability of a typical open-pit lignite mining slope is investigated by ﬁnite element analysis. The precipitation range deﬁned by the environmental analysis is used. The SF is practically constant for the lower rainfall intensities (0.8–2.2 mm/h). For the higher intensities of 6.4 mm/h, 9 mm/h, and 17 mm/h, the SF decrease is almost the same (from 2.08 to 1.9), with reduction rates of 8.3%, 8.9%, and 9.3%, respectively. The effect of the critical geotechnical properties—groundwater table depth, unsaturated zone properties, and soil permeability—is also examined for a complete evaluation.


Introduction
Coal and lignite have contributed to global energy needs for several decades, accounting for 30% of the global energy production in 2010 and 27% in 2020 [1]. The transition to cleaner energy has led to mine closures in Europe and globally. With the shift in the energy mix and the groundbreaking changes in the economy's structure, coal regions should not be left behind, and the economic and social impact of mine closures should not be ignored. Decarbonization will leave vast areas of closed mines and waste dumps that burden community unless suitably reclaimed. For the reclamation of these areas, their safety and stability are crucial, as they are jeopardized by several factors. Rainfall is typically such a crucial factor, particularly given the effect of climate change on its evolution.
Intense rainfall and extreme weather events are often associated with landslides and catastrophes on natural or engineered slopes. The analysis and design of such slopes are currently largely based on numerical modeling. Several numerical methods have been developed for slope stability analysis, with two dominating the field: the limit equilibrium method (LEM) and the finite element method (FEM). These stability methods are typically coupled with a finite element analysis considering groundwater evolution during rainfall infiltration. As a result, many works, e.g., [2][3][4][5], use the FEM for groundwater, rainfall, and stability, while others, e.g., [6][7][8][9], combine the FEM and the LEM.
The influence of rainfall on the slope stability of open-pit coal and lignite mines has long been identified and considered for design purposes, e.g., [10][11][12]. The investigation of precipitation's impact on mining slope stability has been primarily based on observations and empirical correlations between precipitation rates and increasing slope movements, e.g., [13]. This type of slope analysis, although practical, has only incorporated rain in indirect ways, e.g., the rise of the groundwater table rise or the increase of pore water pressures.
In contrast, slope stability analysis of cut slopes (civil infrastructure) has progressed tremendously, employing unsaturated soil mechanics and advanced hydrological concepts to incorporate precipitation into numerical modeling, e.g., [2,5,[14][15][16]. However, even in these cases, precipitation has been considered through parametric analyses, with precipitation characteristics usually dictated by general considerations.
Differences among the two types of slopes (cut slopes in civil infrastructure versus overburden soil slopes in coal mines) refer mainly to: the enormous height of the mines' excavations, up to 200 m, versus heights up to 30-40 m for cut slopes; the smooth inclination usually utilized for excavations (versus steeper slopes often encountered in infrastructure); the large tolerance in displacements in mining slopes (even in the order of magnitude of meters) versus more strict limitations for typical cut slopes. Additionally, the soil types of cut slopes are frequently different from the excavation materials, including lignite deposits and seams and spoil materials, presenting considerable variability.
A novel methodology is presented to evaluate coal and lignite mine slope stability due to precipitation. The stability of these areas is essential, especially during the current period of decarbonization and movement towards the post-coal and post-lignite era. Areas that hosted coal and lignite mines are now being reclaimed, e.g., by installing renewable energy systems, and new energy principles relating to energy production, harvesting, and storage are being set. The appropriate exploitation of these areas is directly related to their safety and stability in sustainably contributing to the new era.
The proposed framework is formulated to be used in various conditions. An indepth climatic analysis appropriately defines the parameters and specifies the range of rainfall intensity for stability purposes. Initially, the baseline rainfall is established from reported climatic conditions for the last thirty years and extreme events of the last ten years. Furthermore, climatic projections are used to anticipate trends and possible changes in the baseline, past, and current rainfall conditions. Finally, a typical slope representing a mine excavation is simulated within the anticipated climatic conditions. The finite element method (FEM) is employed, and a fully coupled flow-deformation analysis is conducted that allows for the simultaneous evaluation of the unsaturated transient seepage flow and the mechanical response of the soil. A parametric analysis is conducted for the critical geotechnical parameters. Greek lignite areas are used as a case example; however, the methodology and insights apply to several similar surface slopes. The numerical model (geometry, properties) can be implemented similarly in many countries where similar conditions are found. The results of this work can be used as a guideline for similar situations.

The Methodology and the Case Example of Greek Lignite Mines
The present study presents a methodological framework including the assessment of past and current rainfall, future rainfall projections, and numerical analysis of the effect of rainfall on the slope stability of a case example. Figure 1 presents an overview of this work's methodology. This flowchart focuses on the slope stability of surface coal mines at the center and includes all aspects as conceptualized herein. Table 1 summarizes the main parameters of each aspect for some previous characteristic studies, analyzed in the following: g a is a parameter of the soil-water characteristic curve; ϕ and c are the effective friction angle and cohesion; k s is the saturated permeability; H is the height of the slope; q is rainfall intensity; H w defines the groundwater table's depth. The upper part of Figure 1 describes the present work's approach to evaluating precipitation, while for the other aspects, the most critical parameters are included in a parametric analysis. Unlike previous works where precipitation was estimated based on global data [2,3,5,9,[16][17][18][19][20][21][22][23][24][25] or was derived for a specific period as a particular time series [3,14,[26][27][28][29][30], the way to assess the precipitation range is through past events and climate projections of the investigated area. This approach provides a clear advantage in terms of how the rainfall range is calculated. Greek mining areas are utilized to quantify the precipitation range because it is area-specific, and the methodology followed needs specific to the areas to be implemented.  The lower part of the flowchart describes the slope (geometry, groundwater, soil, and unsaturated zone properties). As briefly mentioned in the introduction, the geometry of surface coal and lignite mines differs from embankments, cut, or natural slopes. Previous studies have mainly focused on engineered slopes less than 30 m high, with one work studying a slope reaching 40 m [5]. Additionally, three works have analyzed natural slopes of more than 50 m, but these present case studies and not systematic analyses of this type of slope [3,27,28], thus ignoring methodological issues and not evaluating the effect of the various parameters. Moreover, the groundwater conditions, crucial for the stability, differ in surface mine slopes; in the literature, a narrow range is commonly employed for the depth of the groundwater table, usually evaluated as horizontal or with a small inclination. In surface mines, the water reaches the toe of the slope, and inside the slope, the steady-state flow of the water varies significantly due to varying initial conditions and drainage wells. These are considered in the analysis that follows through a steady-state groundwater flow with largely varying groundwater depth.
Finally, soil properties, including the unsaturated zone, are case-and site-specific. Due to the lignite basins' creation process, it is common to observe fine-grained materials overlying and mixing with the lignite layers and seams. In this work, Greek lignite mines are employed for the soil properties that represent fine-grained materials. Theocharis et al. [31] recently discussed an overview of the geotechnical characterization of the soil materials found in Greek lignite mines. Nevertheless, these properties are not specific to a site but characterize this type of surface mine slope.
Greek mines have been systematically exploited and were made extremely large and deep to increase lignite excavation production and efficiency. Exploitable lignite deposits are principally located in Western Macedonia (northern Greece) and secondarily in the Peloponnese (southern Greece) ( Figure 2). In the first case, the central basin is Ptolemais-Amyntaio, while in the Peloponnese, the lignite basin is at the Megalopoli area. Greek lignite mines present a ground profile with an overburden zone of sterile materials (soilrock deposits of several tens of meters) and a lower zone of exploitable lignite deposits. The sterile material primarily consists of marls and stiff or even hard clays, while weaker cohesive soils and sands are also observed, e.g., [32,33]. As a result, soil mechanics (rather than rock mechanics) principles are central to stability issues.
Overall, several researchers have analyzed the geologic formations, e.g., [33], and ground profiles of Greek lignite mines from a geotechnical perspective, e.g., [32,34]. Due to the immense excavations, large and deep slopes are created during open-pit mining. Thus, slope stability is critical for safe and efficient lignite mining and the valorization of lignite mining areas after closure and abandonment. Moreover, the geotechnical characterization of the soils found on the slopes is vital for slope stability analysis. Several past failures demonstrate that slope stability has not always been rigorously evaluated, resulting in problems in the mines' safety and management [34], often related to rainfall events.

Baseline Climate Description-Past Events
Greece is located in the southern part of the Balkan peninsula, in southeastern Europe. It mainly presents a typical Mediterranean climate ( Figure 2) with mild, rainy winters, warm, dry summers, and extended periods of sunlight [36]. The area of Megalopoli is representative of this climate; in contrast, the region of Western Macedonia presents a climate similar to continental Europe [37].
The climate of Megalopoli is categorized as Csa (hot summer Mediterranean) according to the Köppen-Geiger climate classification, having hot summers and mild winters [35] ( Figure 2). Most rainfall days are recorded during fall and winter with an annual average precipitation of 800 mm; the driest month is July, and the wettest is December, with an average monthly precipitation of 9 mm and 147 mm, respectively [38]. The average temperature in Megalopoli is 14.9 • C.
The climate of Ptolemais is classified as Cfb (temperate oceanic) according to the Köppen-Geiger climate classification, having intense winters and mild summers with a very low number of dry days (Beck, 2018) ( Figure 2). Rainfall in plain areas ranges from 600-800 mm/year [39]. There are 108.3 average rainfall days, with August being the driest month, with 30 mm of rain, and November the wettest, with 60.3 mm of average precipitation. The mean annual temperature is 11.5 • C.
The quantification of past and current climate conditions depends on the measurements obtained from several weather stations in an area. Figure 3 presents the precipitation for these two regions for the last 30 years (1990-2019) (data provided by the Hellenic National Meteorological Service). Both regions present fluctuations, but they are practically constant in these 30 years. Megalopoli's precipitation level is overall higher and presents more sharp fluctuations than the Ptolemais one. For example, the maximum mean monthly precipitation difference between two consecutive years was 33 mm for Ptolemais (2013-2014) and 51 mm for Megalopoli (2004)(2005). For Ptolemais, the mean precipitation per month is 37 mm, and for Megalopoli, 56 mm. Assuming this precipitation falls in one day, the rainfall intensity is 1.5 mm/h in Ptolemais and 2.3 mm/h in Megalopoli, denoting light rain. By obtaining the maximum mean monthly precipitation of 64 mm in Ptolemais and 88 mm in Megalopoli and making the same assumption, one obtains 2.7 mm/h and 3.7 mm/h, respectively, denoting a moderate rain intensity. More detailed data for the precipitation in these areas were obtained from the National Observatory of Athens up to ten years ago, for two meteorological stations different but close to previous ones (data publicly available). They do not cover the 30 years that is the classical period used for describing a climate (World Meteorological Organization) but are more detailed in terms of rainfall activity. Figure 4 presents the maximum mean monthly precipitation per year (the x-axis refers to the month with the highest mean precipitation) and the maximum daily precipitation per year (the x-axis refers to the day with the highest precipitation) for the past ten years. Figure 4a shows that the maximum mean monthly precipitation per month can be four times larger than the mean monthly one. For instance, in 2019, for Ptolemais, the mean monthly precipitation was 30.8 mm and the maximum 92.6 mm, and for Megalopoli, 88 mm and 381 mm, respectively. These observations denote that more extreme rains than assessed by the mean monthly values are possible; if these maximum mean monthly rains are assumed to fall in one day, that would dictate 4 mm/h and 16 mm/h. Figure 4b presents maximum daily rains similar to the mean monthly precipitation, supporting the assumption that the mean monthly precipitation can fall in one day. For example, in 2015, Ptolemais' mean monthly precipitation was 59.2 mm, and the maximum daily precipitation was 53.4 mm, while in Megalopoli, the corresponding values were 57.1 mm and 54.8 mm, respectively. One exception denoted by an extreme rainfall event in Megalopoli in 2016 (when more than 200 mm of precipitation fell in one day) signifies that extreme events, where the maximum mean monthly rainfall falls at one time, are also possible. In that case, 213 mm of precipitation within one day corresponds to 8.9 mm/h, a heavy rain intensity. Due to its severe intensity, this event was also documented in the local press.

Climate Projections-Future Events
Climate change's impact on various systems and activities tends to increase incrementally with each day. Thus, it is necessary to anticipate its effects, especially concerning temperature and rainfall events affecting mining operations and reclamation, infrastructure maintenance, and management. Climate change directly influences temperature and precipitation in an increasing or decreasing manner, depending on the study area. Various aspects of climate change are analyzed in the literature due to its importance and broad implications, e.g., [40].
Specifically for the mining sector and the infrastructure, climate parameters should be projected in the future to prevent failures and catastrophic events. In this way, extreme future events can be considered for planning operations and designs. In this section, future rainfall and temperature are estimated for the Ptolemais area, with specific emission scenarios being the key ingredient of the analysis. The focus is on the climatic parameters' trends to evaluate if the past conditions will remain valid in the future.
The World Climate Research Programme (WCRP) is the framework utilized in this work for climate projections; it provides tools that monitor, simulate, and project global climate conditions with high accuracy, being a rigorous framework built through international science coordination and partnerships. WCPR established the Coordinated Regional Downscaling Experiment (CORDEX) to produce regional climate projections for all continents [41]. The Earth System Grid Federation, a global consortium of data providers and data archives with a shared infrastructure, is used to access the CORDEX data [42].
Climate evolution can be simulated by global climate models and regional climate downscaling techniques. A global climate model can produce reliable projections on how average conditions will change in an area larger than 1000 km by 1000 km over the coming decades [43]. It is based on fundamental laws of matter and energy and uses mathematical equations to describe the earth's elements (e.g., atmosphere, oceans, land surface, and ice). Nevertheless, global climate models are inadequate in regions with complex orography [44].
During the last decades, regional climate downscaling techniques have provided higher resolution climate information [45] by combining regional climate models, global climate models, and empirical statistical downscaling. Regional models typically have 10-50 km resolution and detail all the local climate characteristics. A global climate model is always needed to provide the boundary conditions at the regional model's edges. This work uses eight regional climate models (named datasets) to assess the Ptolemais mine area's future climate conditions. These datasets originate from different assumptions regarding fundamental laws and constants. Thus, assessing all eight provides a solid basis for evaluating future climatic trends.
The quantity of future greenhouse gas emissions is the critical element in predicting climate conditions and a source of great uncertainty. For this reason, four scenarios known as Representative Concentration Pathways (RCP) were developed by the Intergovernmental Panel on Climate Change of the United Nations and are systematically used in the literature. The RCPs describe different climate futures, all considered possible depending on the volume of greenhouse gases emitted in the years to come. The present work considers two emission scenarios: RCP 4.5, the median range or stabilization scenario, and RCP 8.5, the high-end, dystopian scenario. The output variables are presented as monthly averages, and of interest are the trends.
The mean monthly temperature projections (not presented herein) present an increase of about 1-2 • C for RCP 4.5 and up to 4 • C for RCP 8.5. These observations agree well with those of previous works and project an increase in temperature of 1-5 • C [46]. The different datasets provide a wildly different range in the initial (2020-2029) mean monthly precipitation ( Figure 5), e.g., from 19 mm (Dataset 1) to 80 mm (Dataset 8) for RCP 4.5. This initial range is due to the different assumptions employed by the different models and does not affect the evolution's reliability. Focusing on the projections, the precipitation level remains practically the same for all datasets. No clear trend of increase or decrease is observed, and current conditions remain the same. As a result, it is presumed that the intensity of extreme events will not change drastically and will not result in violent rains (e.g., more than 16 mm/h for more than 24 h; see the previous section). Notice that this area's precipitation evolution differs from the general overview of the Mediterranean area that projects smaller mean precipitation and more extremes [46]. This disagreement is explained, as local conditions crucially affect the overall evolution and result in various trends. Overall, Figure 3 proposes a mean monthly rainfall less than 90 mm for Megalopoli and 60 mm for Ptolemais; Figure 4a a maximum mean monthly rainfall up to about 400 mm for Megalopoli and 200 mm for Ptolemais; and Figure 5a a future mean monthly rainfall for Ptolemais that can be up to 350 mm. As discussed in the previous paragraph, climate projections mainly reveal trends, not absolute values. As suggested in the previous section, these rainfalls can also represent maximum daily values and extreme events when the mean or maximum rainfall of a month falls in one day. Given this perspective, Megalopoli's precipitation range is approximately 20-400 mm and Ptolemais's 20-200 mm per day, including extreme weather events and incorporating future changes. These precipitation rates can be translated on an hourly basis for simulation purposes as 0.8-17 mm/h for Megalopoli and 0.8-9 mm/h for Ptolemais. These rates are accurate enough for an informed evaluation of precipitation and its effect on slope stability.

Numerical Model Overview
The finite element method (FEM) is employed to investigate how rainfall affects pit mining slope stability. A fully coupled flow-deformation analysis is conducted using a geotechnical finite element software that allows for the simultaneous evaluation of the unsaturated transient seepage flow and soil mechanical response [47]. For slope stability analysis, the FEM needs to be combined with the strength reduction technique to calculate the safety factor (SF); the soil material's strength is gradually reduced until failure [48].
The simulation of rain infiltration requires the consideration of the unsaturated zone that lies above the groundwater table. The soil-water characteristic curve (SWCC) is usually employed for this type of analysis. It indicates the unsaturated soil's capacity to retain moisture at various stress levels and typically relates water content to soil suction. Furthermore, a relationship is needed to express the unsaturated soil permeability at different moisture levels, i.e., the unsaturated soil's drainage capacity. Finally, during rainfall, the precipitation does not completely infiltrate the soil. A possibly substantial part of it flows as runoff water on the soil's surface. The unsaturated soil permeability governs the part of the rainfall that infiltrates, and a small excess of water is considered to stagnate. The rest of the water is runoff water that does not affect stability.
Through the SWCC and the permeability function, the assessment of soil behavior in the unsaturated zone is feasible for most geotechnical engineering applications [49]. In this study, unsaturated soil behavior is captured by a form of the widely known van Genuchten-Mualem model [50,51]. The soil saturation, S w , in the unsaturated zone, is described as a function of the matric suction, u w (negative pore-water pressure), according to: where S r denotes the soil's residual saturation; S s the soil saturation at the saturated state; g a and g n are material parameters governing the shape of the SWCC; and γ w is the water's unit weight (9.81 kN/m 3 ). The residual saturation degree S r and the material parameter g a are critical factors for the unsaturated soil's state [52]. In this work, the residual degree of saturation, S r , was assumed to equal 0.2, a typical value for fine-grained soils, and the reference g a to equal 0.01 m −1 . However, the range of this parameter for fine-grained soil materials of 0.001-0.1 m −1 was considered through a parametric analysis [16,27,30,53,54]. This range slightly differs from the literature, referring more to surface mines with fine-grained soils; thus, it has not been thoroughly examined, especially concerning the precipitation range investigated herein. Finally, the parameter g n is assumed to equal 1.3, representing finegrained soils [55]. Figure 6 presents the SWCC for three different values of g a . As g a decreases, the suction needed to initiate soil's desaturation (95% of the saturation degree) increases. The soil parameter g n is constant, leading to the same gradient of the linear part of the three SWCCs. Finally, the residual saturation S r = 0.2 is constant, and all three curves will eventually reach that value for large suction. Furthermore, unsaturated hydraulic conductivity is also a matric suction function [56]. According to the van Genuchten-Mualem model: where g l is a fitting parameter equal to 0.5, regardless of the soil type. The effective degree of saturation, S e , expresses a normalized water content between the saturated and the residual state of the soil and equals: As a result of Equations (2) and (3), the unsaturated permeability, k w , decreases with the decrease of the soil water content, S w , and suction increases.
Finally, the soil's strength is defined to complete the unsaturated zone's description. The unsaturated soil's shear strength can be written based on Bishop's effective stress as [57,58]: where σ denotes the total stress and ϕ and c are the effective soil shear strength parameters. The second term of the expanded form of Equation (4) (which contains the effective saturation degree, the matric suction |u w |, and the effective friction angle) represents the effect of suction in the unsaturated soil shear strength.

Effect of Rainfall on Slope Stability
A typical 50 m deep pit excavation is considered with a mild slope of 14 • . This geometry represents a deep excavation for civil infrastructure but not so deep for coal and lignite mines. This height is over the limits of the literature regarding the parametric analysis of artificial (cut or embankment) slopes (see also Section 2). A parametric study was performed for the critical parameters of groundwater depth, soil, and unsaturated zone properties (see Figure 1 and Section 1). The groundwater was considered through a steady-state flow, and the effect of its depth was investigated. The effect of the unsaturated zone was examined through the parameter g a , as analyzed in the previous section. Finally, the soil's permeability, the main soil property related to infiltration, was parametrically investigated. The chosen saturated permeability lies in the range of the dominant overburden material of Greek open-pit lignite mines [31]. Figure 7 illustrates the FEM model, including the geometry, the hydraulic characteristics, and the discretized mesh. Two-dimensional (2D) plane strain and drained conditions were used. Vertical boundaries were at a 2H distance from the slope's crest and toe, respectively. The horizontal bottom boundary was at a distance, H, from the slope's toe to minimize boundary conditions; H stands for the slope's overall height. A fine mesh was used for FEM, and mesh density was increased in critical areas. The initial hydraulic conditions (before the rain starts) were established by a steady-state flow, defined by two boundary groundwater heads: the water is located at depth H w at the left border and the bottom of the excavation. Rainfall intensity was simulated by imposing appropriate boundary flux velocities. Different rainfall intensities were used, representing a wide range. A low-to-medium range from 0.8-6.4 mm/h (19.2-153.6 mm/day) and two extreme cases, 9 mm/h (216 mm/day) and 17 mm/h (408 mm/day), are employed based on the climatic measurements and projections analyzed and discussed in the previous sections. The time duration of the rainfall for the simulations is one day (24 h).
For the reference case, the groundwater head, H w , was 20 m, and a medium rainfall intensity was chosen of 3.6 mm/h (86.4 mm/day). All the physical and hydro-mechanical parameters of the reference model's soil material are provided in Table 2. The fundamental parameters employed in the numerical model, including the geometry (height and angle) and the soil shear strength, have been validated to represent actual cases of Greek lignite mines. The parameters that have not been extensively validated are considered in the parametric analysis that follows to obtain the range of the results. The possible spatial and temporal variability of the basic parameters was not considered in this work to retain the practical and general nature of the framework. Nevertheless, if the soil variability is known, this could be involved in the methodology by applying the appropriate changes in the numerical models. Geostatic stress coefficient The Mohr-Coulomb constitutive model with a non-associated flow rule was used. Four critical parameters were varied to evaluate their effect on stability: the rainfall intensity, the initial groundwater head, H w , the soil's saturated permeability, and the SWCC parameter, g a . Table 3 shows the examined values; all other parameters are kept constant. Typical parameters affecting slope stability, such as the soil's strength, are not investigated. They are well-known to affect slope stability but do not relate to the influence of rainfall infiltration. For the reference case and the 3.6 mm/h rain, the average infiltration amount was 28.4 mm (33% of the rain) and the runoff 58.0 mm (67% of the rain). Notice that fine-grained soil was simulated, thus a low percentage of rain infiltration was expected. For the very low rainfall of 0.8 mm/h, this amount rises to 77%, while for the extreme rainfall of 17 mm/h, this amount becomes 19%. An essential aspect of these amounts is that they are averages and do not entirely reflect the suction conditions and the stability. Although 33% of the rain infiltrates the slope for the baseline model, near the toe of the slope, all the rain infiltrates, raising the groundwater and crucially altering stability (i.e., the failure surface shape, failure surface place, and the soil's effective stresses and deformations).
The reference case has a very high initial SF, equal to 2.08; the failure surface ( Figure 8) is a typical circular surface on homogeneous slopes. This high SF is expected, as the reference model presents a mild slope of medium height for mining excavations. Soil strength is enhanced by the artificial cohesion resulting from matric suction (Figure 9a). At the initial stage, the matric suction-and the related artificial cohesion-increase linearly over the groundwater table, reaching 280 kN/m 2 at the slope's crest. As the rainfall infiltrates the slopes, the suction decreases at the slope's edges, where the rainfall infiltrates. Figure 9b illustrates the final matric suction distribution in the slope at the end of the one-day rainfall. As time passes, rainfall infiltrates deeper into the soil, and thus a new distribution emerges with zero suction at the slopes' edges and the groundwater table. For fine-grained materials, such as those investigated herein, the suction reduction is the primary mechanism of rainfall affecting slope stability. Figure 10a displays the effect of the SWCC material parameter, g a , on the pit slope's safety factor for a 3.6 mm/h intensity rainfall event that lasts 24 h (total rainfall height equal to 86.4 mm). The selected values of the SWCC material parameter, g a , cover the range of a fine-grained over-consolidated clayey soil. The safety factor reduction rate increases significantly with the decrease of parameter g a . In particular, for g a values equal to 0.001 m −1 , 0.005 m −1 , and 0.01 m −1 , SF reduces with a rate of 16.5%, 10.0%, and 4.7%, respectively. As g a increases further (0.05 m −1 and 0.1 m −1 ), the slope safety factor is practically constant during rainfall. This difference between the five examined cases is attributed to the shape of the SWCC and the permeability function. Higher g a values correspond to lower suction (lower effective saturation degree) and lower soil permeability values within the unsaturated zone. Thus, as the parameter g a increases, the unsaturated soil permeability, k w , decreases, and rainwater's infiltration into the unsaturated zone becomes more difficult than before. Indicatively, for the extremes g a = 0.001 m −1 and g a = 0.1 m −1 , the values of soil's unsaturated permeability, in the area of the slope's crest, are equal to 10 −5 cm/s and 10 −7 cm/s, respectively, i.e., they differ by two orders of magnitude.   Figure 10b presents the variation of the pit slope SF for different initial groundwater heads; H w ranges from a medium shallow (10 m) to a deep (40 m) groundwater table. Before the rainfall, the initial SF increases with the groundwater's head rise-an expected result, since a high groundwater table corresponds to lower SF, whereas a deep groundwater table corresponds to higher SF. During the rainfall event, the safety factor reduction rate decreases with the increase of the groundwater head. The SF's reduction rate for the groundwater head of H w = 10 m is 9.5%, for H w = 15 m is 7.0%, and for H w = 20 m is 4.7%. For the cases of deep groundwater tables corresponding to the high groundwater heads of H w = 30 m and H w = 40 m, the slope safety factor is practically constant. Again, high suction stresses and low effective saturation degrees develop in the extended unsaturated zones, leading to lower unsaturated permeability, impeding the infiltration procedure. It is noted that the groundwater head of H w = 0 m, corresponding to a groundwater table that starts from the ground surface, was also examined. Since this value corresponds to a high groundwater table, the unsaturated zone's extent is quite limited for large suction stresses to develop.  Figure 11a presents the SF with time during rainfall for various values of saturated permeability; the permeability range covers fine-grained but not impermeable soil materials. The reduction rate of the SF during rainfall increases slightly with the increase of the saturated soil permeability. For higher saturated permeability values k s = 10 −3 cm/s, 5 × 10 −4 cm/s, and 10 −4 cm/s, the SF reduction rates are equal to 5.3%, 5.1%, and 4.7%, respectively. The slope's SF is practically constant for lower saturated permeability values, k s = 5 × 10 −5 cm/s and 10 −5 cm/s. Obviously, as the saturated permeability decreases, soil permeability in the unsaturated zone decreases, and hence water infiltration during rainfall becomes more difficult. The observed minor effect is because of the low values of k s . Figure 11b presents SF with time for different rain intensities lasting 24 h. As the rainfall intensity increases, the SF reduction rate also increases. The SF is practically constant for the lower rainfall intensity, 0.8 mm/h and 2.2 mm/h (19.2 mm/day and 52.8 mm/day, respectively). This almost constant SF is due to the small total rainfall discharge that infiltrates into the relatively large-scale slope of 50 m height; the total quantity of rainwater is inadequate in inducing a significant dissipation of suction stresses and a decrease of soil's shear strength. As the intensity increases further, it affects the SF more clearly. For the highest value of the medium range-6.4 mm/h or 153.6 mm/day-and the two extreme intensities-9 mm/h and 17 mm/h or 216 mm/day and 408 mm/day-the SF decrease is almost the same, and the reduction rates are 8.3%, 8.9%, and 9.3%, respectively. Thus, an upper limit is practically identified for the effect of the rainfall intensity on SF, signifying that more extreme rainfall will not further affect the minimum SF reached.

Discussion
In this work, a practical framework is proposed (see Figure 1) to assess the effect of rainfall on coal and lignite mines, including details on how to quantify rainfall based on measurements and climate projections. For the case example of Greek lignite mines, it is concluded that slope stability is not jeopardized due to rainfall for the range of the parameters involved. Nevertheless, the changes of the SF due to the four parameters examined above are of general interest. All four parameters are essential for stability and they affect the way rainfall infiltrates the soil. The soil-water characteristic curve affects the unsaturated zone and was investigated in this study through g a . The considered range of g a was rather broad, including the whole range of fine-grained materials, and thus its effect, although important, is not critical. Notice, however, that, in practice, having an accurate estimation of the SWCC is rare, and a broad range should generally be examined. The soil saturation complements the SWCC for the infiltration analysis, which has a similar effect, important but not critical.
Furthermore, the groundwater table's height in the soil slope is critical for stability. That is related to the dewatering operations during a mining excavation and the groundwater conditions before the mine creation. The higher groundwater table suggests greater caution with respect to stability and a more significant impact of rainfall on stability. Finally, the intensity of rainfall is crucial. Of particular interest is medium-to-high rainfall intensity, as this range most affects stability. Very low intensity is inconsequential, while heavy rain practically has the same effect as a medium-to-high intensity. Overall, the results are in accordance with those reported in the previous literature, offering some additional insights into the particular case of slopes of coal and lignite mines.
The current framework could be further advanced in future works incorporating more details in each of its components (climate measurements, climate projections, numerical models) according to the demands of different projects. Climate measurements have been evaluated herein based on average values and trends, and generalizations based on a few registered extreme events. A more advanced statistical analysis could increase accuracy but would require more data and experience. Furthermore, climate models for projections are constantly updated, and thus any advances could be incorporated directly. Finally, spatial and temporal variability in soil properties could be implemented in the seepage and stability analysis. This direction increases accuracy if relevant data are available but complicates the numerical modeling and requires extra caution in the interpretation of results.

Conclusions
A new framework has been presented in evaluating coal and lignite mine slope stability due to precipitation. The main advantage of this framework is the comprehensive analysis of climate conditions and rainfall and the numerical slope stability analysis, including all steps in the procedure. Initially, climatic records were documented and analyzed to evaluate the baseline conditions; the focus was mainly on precipitation. Subsequently, the range of rainfall intensity was specified for stability analysis. Rainfall analysis included mean and maximum mean monthly precipitation for the past 30 years and maximum daily rains for the past ten years. This approach was direct and practical; an advanced statistical analysis might have increased accuracy but would have required more detailed data and would have been less practical. The case of Greek lignite mining areas was employed to illustrate the methodology. It was concluded that in the areas under investigation, rainfall intensity varies from very low (1 mm/h) to documented intensities up to 9 mm/h, lasting up to a day (24 h); very few extreme phenomena in terms of rainfall are documented.
Furthermore, climatic projections of future trends were used to anticipate possible deviations from the baseline rainfall conditions. Nowadays, several methodologies are proposed for climatic projections; some open-access tools have been used and described to evaluate these projections in this work. The tools proposed in this work are readily applicable for most areas, open and free to use, and relatively easily implemented. Rainfall projections show that the reference range of past events is not expected to change significantly in the areas under investigation, independent of future emission scenarios. On the other hand, the mean temperature is expected to increase depending on future emissions.
Finally, the effect of rainfall infiltration on the stability of a typical open-pit lignite mining slope was numerically investigated. The finite element method (FEM) was employed for the groundwater and rainfall analysis and was combined with the shear strength reduction method for stability evaluation. A fully coupled flow-deformation analysis was conducted, allowing for the simultaneous evaluation of the unsaturated transient seepage flow and soil mechanical response.
The impact of the critical factors-the SWCC soil parameter, g a , the initial groundwater head, H w , the saturated permeability, k s , and the rainfall intensity, q-were investigated through parametric finite element analyses. Rain infiltration primarily produces dissipation of suction stresses in the unsaturated zone, which induces a decrease in unsaturated soil shear strength, with the subsequent deterioration of pit slope stability. However, it is demonstrated that the impact magnitude of the examined parameters varies significantly. According to the numerical results, the effect of the SWCC material parameter, g a , in SF's decrease rate during rainfall is more important among the other parameters. As g a decreases, the reduction rate of the pit slope's safety factor during rainfall dramatically increases.
The initial groundwater head, H w , and the rainfall intensity, q, showed a modest effect in reducing the pit slope's SF. As H w decreases (higher initial groundwater table location) or as the rainfall intensity, q, increases, SF's reduction rate during rainfall moderately increases. Finally, as the parameter k s increases, the SF reduction rate during rainfall increases slightly. Further work could include a more in-depth sensitivity analysis of the parameters affecting slope stability and various geometries. The particular geometry, stratigraphy, geotechnical, and groundwater conditions should be evaluated for any specific case study.