Link/Page Citation
Author(s): Georgios Louloudis (corresponding author) [1,*]; Christos Roumpos (corresponding author) [1,*]; Emmanouil Louloudis [2]; Eleni Mertiri [1]; Georgios Kasfikis [1]
1. Introduction
In the region of the South Ptolemais lignite-bearing basin, in Macedonia, N.W.( north-west) Greece (Figure 1), the Public Power Corporation (PPC) of Greece had developed extensive mining activity in order to provide adequate lignite quantities for 3000 MW power plants. However, Greece’s commitment to the European Commission’s long-term strategic vision to achieve net-zero greenhouse-gas emissions [1,2] has provided motivations for sustainable economic growth and spurred the modernization of the economy. The PPC operates within the framework of corporate environmental responsibility [3] and has established an environmental policy that reflects its commitment to comply with all applicable national and international laws and support continuous environmental improvement.
Over the last decade, particular policy attention has been paid to the transition from coal energy to clean energy technologies to reduce the risks of climate change impacts and contribute to climate-resilient pathways [4] for sustainable development. The goal is to integrate renewable sources and energy efficiency while simultaneously providing the lowest possible cost. As a result, multi-level planning is critical to designing and implementing effective strategies, exploring repurposing options, and making decisions about a post-mining future and economic transitions, considering regional and policy contexts [5].
Recently, the PPC has investigated and evaluated alternative repurposing options to produce energy using renewable energy sources (RES) and has planned post-mining land uses to be sustainable [6] in the long term. In this context, the PPC is planning to transition from coal-based electricity to renewable energy. Specifically, the PPC’s research is focused on the sustainable post-mining development of different RES plant types in the remaining post-mining landscapes and consequently evaluating whether implementing such technology could fulfill the electricity demand.
Regarding RES technology application, it is essential to consider the intermittent power outcome. In such a case, integrating energy storage technology with renewable energy production could be an effective solution. Therefore, the PPC uses a holistic approach to mine transition planning by examining energy storage technology, such as pumped hydro storage (PHS) technology.
A considerable volume of literature has been published on applying PHS technology in closed mines in combination with renewable energy sources, such as photovoltaic (PV) parks. According to Bódis et al. [7] the post-mining available landscape area for installing PV systems is abundant in the coal regions of the European Union. In addition, solar PV systems could replace the current coal-fired power plants, provided storage-flexible units with sufficient capacities are added. The first innovative project of a hybrid energy system, which combines PHS technology with multi-stage solar and wind renewable energy units in an abandoned mine, is the Kidston Clean Energy Hub Project [8] in Kidston, Far-North Queensland in Australia, which is expected to be in operation by 2024. According to Holcombe, pit voids in abandoned mines could become an asset underpinning RES development at mining sites and creatively reusing them [9]. Regarding economic feasibility, Wessel et al. reported that hybrid PHS–RES projects could be developed and operated at reasonably low costs [10].
Except for the production of electricity using RES, there is also a growing body of literature on other post-mining transitions and repurposing land-use scenarios [3,11,12], such as wildlife habitat, parks for recreation, eco-tourism and rehabilitation practices concerning pit lakes, i.e., development of fishing farms, recreational sport fisheries, angling, diving, hiking around the lake and even using pits for water storage and flood protection.
Therefore, water quality is an important parameter that could enable new beneficial end uses. Poor water quality could be significantly improved by rehabilitation practices, such as permanent flushing with river water, diverting seasonal floods into the pit lake [13,14], or even permanent river flow-through [15], which can mitigate the risk of eutrophication. Other rehabilitation practices have been reported in post-mining areas, such as appropriate landscaping of the pit lake banks with gentle sloping to emulate natural lakes if possible, filling the bottom of the open pit with alkaline substances such as re-excavated waste lime materials and application of riparian zone vegetation [14,16,17]. These rehabilitation practices could help the development of crucial aquatic biodiversity. Thus, a closure monitoring program is a prerequisite in post-mining areas to ensure long-term pit lake quality, since improving long-term environmental outcomes [18] could enhance pit lake ecosystem development and mitigate environmental damage.
Open pits include large areas that could also be used for waste disposal. Woodlawn is another innovative hybrid electric power project and an example of a circular economy, where the abandoned open pit was repurposed and is now converted into a bio-reactor landfill [12]. In this example, municipal wastes are disposed into the pit, and bio-gas power is produced, feeding the electricity grid. At the same time, other land repurposing is also applied in the surrounding area of the mine. It includes the development of RES technology using wind and solar farms to generate electricity, agricultural operations for farming practices, aquaculture and horticulture operations utilizing captured waste heat in hydroponic horticulture and to cultivate fish, and also the operation of a biological treatment facility to extract the organic content from waste and select recyclable materials or produce compost [19].
Another post-mining land-use transition is transforming mining landscapes into memorial pathways and places of remembrance, demonstrating historical artifacts to memorialize the local mining legacy [12]. In addition, these landmarks could become symbols of the region and thus help the development of tourism and recreation in the area [20].
This study aims to provide a case analysis to examine and evaluate land use repurposing scenarios in a mined-out pit of a closed mine. For this purpose, the mined-out pit of the already closed lignite mine of Kardia, northern Greece, was used as the case study. According to the planning studies for the closure of the lignite mines of western Macedonia [6], rehabilitation includes the development of lakes inside the open pits, photovoltaic (PV) parks developed in the surrounding area, areas used for agriculture or recreational activities and forest areas. In the Kardia mine, the construction of a PHS facility is proposed to work as an energy storage unit due to the major setback in the performance of PV panels, i.e., the difficulty in adjusting to power demands at all times [7,12,21,22]. More specifically, the pit lake is repurposed as a PHS plant. It is used as an energy storage unit to supplement variable renewable capacity and enhance the proposed photovoltaic panels’ performance, increasing the PV park’s capacity factor. The conceptual idea is to utilize the pit lake that will start to develop after the end of mining activities, ensure the cessation of dewatering measures (dewatering wells and pumping stations) and create another or use an existing upper reservoir at higher altitudes. Then, taking advantage of the head difference between the pit lake (lower reservoir) and upper reservoir, the whole system could be used as a storage battery, releasing power when needed to compensate for renewable energy fluctuations.
In this framework, the feasibility of applying PHS technology is examined concerning the prediction of the natural development lake inside the open pit. In addition, considering the lake’s potential growth, the application of different land repurposing scenarios in the Kardia mine is evaluated. Finally, all potential repurposing land uses are evaluated regarding the energy transition potential in the specific region, long-term socio-economic viability and successful adaptation. The intent of this paper is to recommend strategies for improving post-mining land use transition and transforming abandoned surface lignite mines—including pit lakes—into post-mining assets for the communities. In addition, our intent is to contribute to further advancing successful pit lake end uses, as pit lakes could become assets to the clean energy sector.
2. Materials and Methods
2.1. Research Area—Morphology of the Basin
The Ptolemais Lignite Field is located in Kozani Prefecture in north Greece, surrounded by the mountains Vermion, Askion and Skopos (Figure 1). The northern side of the basin is bounded by Komanos horst. The Public Power Corporation of Greece has exploited three lignite reserves in this area: Mavropigi, Kardia, and South Field. The Kardia deposit lies in the middle, with exploitable reserves of 273 million tons (since January 2020) that hosted a 1250 ?W lignite-fired power plant for almost 40 years. Inside dumping was conducted mainly with spreaders (continuous surface mining) and, in some cases, with trucks and loaders (non-continuous surface mining). The internal waste dumps are estimated in the order of ˜1.0 billion m[sup.3]. Currently, the region’s economy depends on agriculture and energy fuel reserves.
The deepest part of the Kardia mine excavation was +450 m a.s.l. when lignite exploitation activities terminated in 2021. By the end of 2022, backfilling of the bottom of the mine pit is planned to achieve slope stability conditions, and thus, the bathymetry of the Kardia excavation will alter, with the deepest part forming at approximately +523 m a.s.l. (Figure 2). The range of ground elevations in the pit is between ˜+523 m and +670 m a.s.l.
2.2. Hydrological and Hydrogeological Setting
The broader research area is part of larger lignite basins of western Macedonia, formed during the Pliocene. The region has a complex of natural lakes (Chimaditida, Zazari, Petron and Vegoritis) (Figure 1). Vegoritis Lake is located in the homonymous closed hydrological basin in western Macedonia, Greece, and is the final receiving body of the basin runoff, with a total area of 2000 km[sup.2]. Vegoritis Lake is also the surficial feature of a major karstic aquifer [23]. Therefore, any alteration in the surface or groundwater conditions within this area would affect the level of the lake. The primary water discharge into the lake is achieved via the Soulou River [24,25]. Therefore, the main priority of the PPC is not only the protection of the natural lakes but also the monitoring of the artificial lake evolution developed in the open voids of the closed lignite mines, as well as the evaluation of the chemical properties of lake water and of the new aquatic ecosystem that is formed.
Kardia open pit is an exhausted mine located on the western margins of the South Ptolemais lignite-bearing basin (Sarigiol basin) in the foothills of Askion Mountain. Quaternary sediments are met on the broader area, which consists mainly of intercalations of sands, clays, sandy marls and conglomerates. These formations overlie Paleozoic schists, ophiolites and granite formations. In addition, Mesozoic dolomitic limestones, volcanic sediments and flysch overlay the Paleozoic formations [6,26]. The lignite sequence underlies the Quaternary sediments and consists of impermeable Pliocene formations of marls and lignite. Lignite sequence formation is significant for the South Field Lignite Mine (LM) and Kardia LM operations due to the potential slope instabilities associated with sliding along a sub-horizontal lignite-to-clay/marl interface, mainly caused by the groundwater pressures [27].
Regarding the hydrolithological conditions in the area of the Kardia mine, two aquifers are generally recognized, the upper (Sarigiol aquifer) and the lower (Komanos aquifer). These aquifers are separated via the impermeable lignite sequence. The upper unconfined aquifer consists of high permeability layers (sands, conglomerates) with intercalations of low permeability layers (clayey interlayers, fluvial deposits). This unit consists of a unique hydraulically connected network that behaves as one aquifer [24].
The hydrogeological basin of Ptolemais South Field (Sarigiol basin) is surrounded by the Vermion (east side), Askion (west side) and Skopos (south side) mountains. The northern part of the basin is bounded by the Komanos tectonic horst. Most of the streams originate from the surrounding mountains, where the general water flow direction is towards the Soulou River. Part of the surface run-off is discharged into the Soulou River via the consolidated sediment channels formed on the ground surface, while part of the water is infiltrated inside the sediment formations. The stream network and the Soulou River are presented in Figure 3. The recharge area of the Sarigiol aquifer is mainly the Vermion and Askion foothills conglomerates and talus cones, whereas the Soulou River is the main outflow channel of the basin, which directs the water flow into Vegoritis Lake, north of the Kardia lignite mine (Figure 1).
Regarding the Sarigiol aquifer, the piezometric surface is conducted by considering the water levels measured in the boreholes. Based on the piezometry in the area of interest, the piezometric level is met at elevations between ˜+640 m and +660 m a.s.l. (Figure 3). The predominant direction of groundwater flow in the Sarigiol aquifer body south of the Kardia mine is generally NW to SE, while locally, next to the open pit Kardia mine location, the direction that prevails is S to N. It is noted that in the same area, the lignite sequence is met at elevations between ˜+350 m and +500 m a.s.l., while during excavations inside the Kardia open pit, the bottom of the lignite sequence was not revealed. Thus, in the conceptual model of Kardia pit lake formation, the simulation was oriented only in the aquifer above the lignite sequence, i.e., the Sarigiol aquifer.
In Figure 4, the methodological approach is presented briefly illustrating the modeled processes.
The renewable resource of the aquifer is the rainfall infiltration through the Quaternary sediments. Meteorological data (temperature, rainfall) were obtained from Kozani WMO No 16 632 weather station, recorded by the Hellenic National Meteorological Service, and spanning a period of 44 years, from 1970–1997 and 2006–2021 (Figure 5). This station is approximately at the same elevation as the Ptolemais lignite field, located 14 km south of the Kardia LM. Meteorological data were also obtained from Ptolemaida LG26 weather station, recorded by the National Observatory of Athens (?OA), and spanning a period of 15 years, from 2006–2021.
According to the meteorological stations’ data, the study location in the Kardia lignite mine, in Ptolemais, had a Köppen climate classification of Cfa (moist with mild winters, wet all seasons with a hot summer), which is consistent with the updated 1-km resolution global maps of the Köppen–Geiger climate classification [28]. The average annual rainfalls at Ptolemais were 498.1 mm (±132.1 mm s.d.), and November was the wettest month with average monthly precipitation of 56.7 mm, while August was the driest month with an average monthly precipitation of 31.7 mm. The average annual temperature was 13.3 °C (±0.8 °C s.d.), with January being the coldest month, with an average temperature of 2.4 °C and the highest temperature of 6.1 °C. According to the meteorological data, there is a visible trend of long-term increase in temperature, presenting an average increase of approximately 2.2 °C over the last 50 years. Regarding precipitation, there seems to be an increasing trend during the period of the last 30 years.
Thornthwaite’s water balance method [29,30] was used to predict the mean annual runoff for the unsaturated zone in the pit’s catchment area and the Sarigiol hydrological basin. This method can provide mean annual estimates in close agreement with measured values [31]. The unsaturated zone’s water balance during the time interval (t,t+dt) is described through the moisture content B. The moisture content B(t+dt) at the end of time dt is described by the following equation [30]:(1)B(t+dt)=B(t)+(N-V-R-SL)dt where all quantities are expressed to the unit area, and B(t) is the water content stored inside the unsaturated zone at the time t (mm), ? is the precipitation rate (mm/month), V is the real evaporation (mm/month), R is the surface runoff outflow (mm/month), SL is the percolation rate (mm/month).
In the calculations, it was considered that deep percolation takes place when soil moisture reaches 70% of the field capacity (FC) [32], where FC is the amount of water remaining when downward drainage has markedly decreased. Therefore, the following equations were applied for the different cases [30]:(2)?(t)+(?-V-R)dt=0.7FCthen(SL)dt=0 (3)F=?(t)+(?-V-R)dt=0.7FCthen(SL)dt=0.1(?-V-R)dt (4)?(t)+(?-V-R)dt=FCthen(SL)dt=?(t)+(?-V-R)dt-FC
The real evaporation V was estimated as a portion of ?hornwaite potential evaporation P. It was assumed that real evaporation takes place when soil moisture reaches 70% of the field capacity, and thus, the following equations were applied for the different cases:(5)?(t+dt)=0.7FCthenV=P (6)?(t+dt)=0.7FCthenV=P[0.188+2(B/FC)-1.2B/FC[sup.2]]
The potential evaporation P is derived by the following equation:(7)P=1.610 tn/I[sup.a] where t[sub.n] is the mean temperature of month n, I is a temperature coefficient of monthly temperature coefficients i[sub.n] summation (I=?n=112i[sub.n]), where i[sub.n]=tn/5[sup.1.514] and a value can be determined as a polynominal function of the temperature coefficient I, where a=0.49239+(1792×10[sup.-5]I)-(771×10[sup.-2]I[sup.2])+(675×10[sup.-9]I[sup.3]).
In addition, in the hydrological model it was assumed that real evaporation cannot exceed the available water derived from precipitation and soil moisture and thus, (V)dt=(?-R)dt+?(t).
Considering the climate conditions in Greece, for the specific area, it is expected that soil moisture content reaches FC in April and thus, regarding the initial conditions of the unsaturated water balance at the time t=0, then ?(0)=FC for April.
Sensitivity analysis performed towards field capacity FC for the unsuturated zone water balance Thornwaite model revealed the following linear relation between FC and percolation coefficient asSLpercentage (%) of rainfall, which is taken as a percentage of precipitation, presenting a coefficient of determination R[sup.2]=0.964 (Figure 6):(8)SL=32.429-0.0903FC
This relation helped in transforming FC values (deduced from NDVI values) into percolation rate values SL.
According to the results (Table 1), for the hydrological cycle in the Sarigiol basin (2006–2021), the mean annual precipitation that the loose sediment basin receives is approximately 584 mm, while the real annual evaporation is estimated to be approximately 425 mm. Also, based on the 15-hydrological year annual water balance output, only 86 mm per year percolates down and recharges the unsaturated zone, which constitutes the long-term renewable recharge of the aquifer. Therefore, there is only one discharge from the Sarigiol basin, i.e., the surface runoff that occurs through the outflow of the Soulou River. In different historical time periods, either the Sarigiol aquifer discharges into the Soulou River or the Soulou River recharges the Sarigiol aquifer, depending on the groundwater level.
2.3. Conceptual Model for the Pit Lake Long-Term Water Balance
In order to calculate the spatiotemporal evolution of the lake inside the Kardia open pit, a simplified conceptual hydrogeological model was developed (Figure 7), following the hydrogeological conceptual model of Louloudis et al. [33], in which the different origins of water inflows into the pit lake’s body are presented. The inflows include (a) the groundwater inflows from the slope seepage face, (b) the surface runoff into the open pit, considering the open pit as a catchment area, and (c) water losses from the lake surface that occur due to direct evaporation processes which result in reducing the net recharge over the lake surface area (i.e., the hydrological processes of precipitation minus evaporative losses over the lake surface area).
2.3.1. Groundwater Inflows
Since the dewatering measures adopted for mine protection will cease after the mine closure by the end of the year 2022, the inflowing groundwater will then begin to increase. As a result, the open pit will fill with groundwater rebound, and the groundwater level will rise. The discharging water derives from the upper aquifer (i.e., above lignite stratification) and the confined aquifers inside the recently developed internal dumps.
For the initial conditions of the model in the year 2022 (q[sub.2022]) the groundwater inflow was calculated by Darcy’s law [34] to be around 1.7 × 10[sup.6] m[sup.3]. This quantity was cross-validated by the records of the water sump pumped quantities during the summer period in 2020, considered as baseflow of surficial hydraulic protection. Specifically, approximately 1.2 × 10[sup.6] m[sup.3] of groundwater inflow is pumped from the sumps, and 0.5 × 10[sup.6] m[sup.3] is pumped from the dewatering wells during the summer period. In addition, the seepage face of groundwater inflow into the excavation is also subjected to an evaporation procedure at mine’s benches. For the Kardia mine, it is estimated that 70% of this quantity is evaporated; thus, only 30% of the groundwater inflow finally reaches and recharges at the pit lake’s water body (i.e., around 0.6 × 10[sup.6] m[sup.3]).
For each year i of the prediction period, the groundwater inflow q(i) was estimated with respect to the surface lake water level h increase. More specifically, the groundwater inflow rate Q per unit face length A is calculated according to the following equation [35]:(9)q(i)=Q/A=K(H[sup.2]-h[sup.2])/2L where K is the hydraulic conductivity, H is the aquifer’s thickness, L is the influence radius of the open pit and h is the height of the lake water level calculated from the surface of the lignite sequence (considered as impermeable formation).
Hydraulic conductivity K ranges from 1.03 to 25.92 m/day, according to pumping tests analysis [36], with an average value of 13.48 m/day. The upper aquifer’s thickness H is estimated in the order of 100 m [24,36] and thus, transmissivity T presents an average value in the order of 50 m[sup.2]/day on the west side of Kardia excavation and a range of 200–1200 m[sup.2]/day in the south region, while the average hydraulic gradient towards the excavation is approximately 2%, as depicted in Figure 3.
To delineate the groundwater inflow rate in the future, groundwater attenuation is considered due to the excessive irrigation consumptive use in the area. In the methodology followed for estimating the aquifer’s depletion, the groundwater recharge of the Sarigiol aquifer around Kardia excavation was considered and the specific consumptive use SCU of irrigation water. The renewable recharge of the aquifer was estimated via Thornthwaite’s water balance method (Table 1), while the SCU was derived from the evapotranspiration ET.
For the calculation of the crop water needs P[sub.C], i.e., the water amount needed to meet the water loss through evapotranspiration, the following general relationship between the reference crop and the crop actually grown was used:(10)P[sub.C]=Kc×ETo where ETo is the reference crop evapotranspiration calculated by using the Blaney–Criddle method (mm/month) [37], while Kc is the crop factor derived from the normalized difference vegetation index (NDVI) [38] which was estimated equal to Kc=0.58:(11)ETo=30p(0.46T[sub.mean]+8) (12)Kc=1.18NDVI+0.04 (13)NDVI=(NIR-RED)/(NIR+RED) where T[sub.mean] is the mean daily temperature (°C), and p is the mean daily percentage of annual daytime hours depending on the latitude of the area. In addition, satellite-based NDVIis computed from the red and near-infrared bands of the satellite by processing a satellite image acquired in 2018, as presented in Figure 8.
The crop water need P[sub.C] is supplied to the crops by a combination of irrigation and rainfall. Therefore, the irrigation water need (IN, mm of water column/month) is the difference between the crop water need P[sub.C] and the effective rainfall Ne, i.e., the part of the rainfall effectively used by the plants:(14)IN=P[sub.C]-Ne where Ne is estimated based on an empirical relationship presented by the Greek Land Reclamation Service of the Ministry of Agriculture [39] for the irrigation period from May to October:(15)Ne=N-(12+0.125N) where ? is the precipitation rate (mm/month).
Subsequently, the specific consumptive use SCU (m[sup.3]/m[sup.2]) was derived from the irrigation water need IN (mm/ month), using the methodology used by Louloudis and Stathopoulos [40]. A statistically significant correlation between irrigation water SCU (m[sup.3]/m[sup.2]) and NDVI is shown in the Figure 9, presenting a coefficient of determination R[sup.2]=0.781. The linear regression analysis was performed by using the electric power consumption of irrigation pumps in the area of the Sarigiol basin for 2018 and also by using a Landsat image acquired in July 2018 (Landsat-8: 15 July 2018). Based on the results, the following relationship between the irrigation water SCU (m[sup.3]/m[sup.2]) and NDVI was used for the calculation of the SCU spatial distribution (Figure 8) for the Sarigiol basin area:(16)SCU=1.539NDVI-0.150
In addition, the deficit of groundwater recharge due to excessive irrigation consumptive use of groundwater near Kardia lignite field is presented in Figure 9, and its spatial distribution for the Sarigiol basin area is shown in Figure 8. The aquifer’s annual drawdown due to the irrigation for the period from July to September was calculated as a difference between annual percolation and irrigation water need, divided by the storage coefficient considered equal to 0.1.
2.3.2. Surface Runoff
Surface runoff is the most essential inflow to the lake’s body. Surface runoff also appears to be the more sensitive parameter in the lake’s filling model so there is a need for precision in surface runoff coefficients introduced to the model. In order to achieve a reliable estimation of the surface runoff coefficient, Cr, in the hydrologically ungauged Kardia area, CN European raster file [41] was processed and the following equations were used [42,43]:(17)Cr=Pvol-0.225,400CN-254[sup.2]/P[sub.vol](P[sub.vol]+0.8(25,400/CN-254)) (18)P[sub.vol]=I[sub.TC]T[sub.C] where P[sub.vol] is the volume of rainfall in mm, CN is the curve number, I[sub.TC] is the mean rainfall intensity (mm/h) for the time of the concentration T[sub.C] (h) which is described by Giandotti equation [44].
In Figure 8 spatial distribution of the runoff coefficient using CN European raster file is presented. Based on the results, inside the active mine area, high values of the runoff coefficient are calculated, in the order of 30% to 40%, whereas outside the perimeter of the excavation, in the zone of the loose sediments between Kardia excavations and Askion Mountain the calculated values hardly exceed 15%. In the western area on the consolidated sediments of Askion Mountain, a runoff coefficient of 20% in magnitude is calculated, while northern of this area near Mavropigi village insignificant values in the order of 5% are revealed. Also, Lefkotopos and Lagorema streams in the region of Skopos Mountain consolidate carbonate sediments with a considerable runoff coefficient of 10% up to 20% in magnitude.
2.3.3. Water Losses from the Lake’s Surface
There are many expressions to evaluate evaporation from the lake’s water surface. The well-known, and highly accurate, Penman–Monteith approach involves many variables (air temperature, solar radiation, relative humidity, wind speed) [45,46,47] which in some cases may not be available. In these cases, adjusted equations with fewer weather parameters are used. Eagleman proposed a simple correlation between these parameters [48]:(19)?[sub.A]=0.035e[sub.S]100-RH[sup.12] where ?[sub.A] is the adjusted pan evaporation in inches per month (free surface evaporation), e[sub.S] is the saturation vapour pressure [49] corresponding to the mean monthly temperature in millibars (mb), and RH is the mean monthly percent relative humidity:(20)RH=(e[sub.S]/e)×100 (21)e[sub.S]=6.11exp[17.269T/(273+T-36)] (22)e=6.11exp[17.269T/(273+T-36)](RH/100)
In the case of the Kardia mine, mean values are used, i.e., the mean monthly temperature (°C), the mean monthly rainfall (mm) and the mean monthly relative humidity (%).
2.3.4. Spatiotemporal Development of the Pit Lake
A crucial research question that arises from this study is whether the pit lake can be naturally filled in a reasonable time frame or is there a need for artificial recharge by flooding with river water. Thus, the spatiotemporal development of the water level inside the Kardia pit void was assessed by generating forecasts of the dependent variables (rainfall and temperature) via linear (autoregressive integrated moving average—ARIMA) and non-linear (artificial neural network—ANN) models, following the same methodology described by Louloudis et al., for the Amynteon mine pit lake [33].
2.4. Energy Transition Strategies in the Post-Mining Era
In the Kardia mine, the PPC energy transition strategy is based on the PV parks that will be developed in the surrounding area, inside and outside the mining exploitation area. Thus, energy storage units are required to supplement the variable renewable capacity. In this mine, the best conditions prevail for a PHS facility installation thanks to its morphology, high exploitation depth of approximately 150 m, and sizable pit lake volume. Therefore, a PHS facility could be used as a storage battery and would provide sufficient storage capacity to compensate for renewable energy fluctuations.
Analyzing the electricity demand is a prerequisite for optimal planning of a hybrid PHS system. Figure 10 presents a representative scheme for the hybrid pumped hydro storage (HPHS) plant operational schedule during daytime hours for the case study area of the Kardia mine. This scheme presents the integration of RES electricity production and the utilization of a hybrid plant. A linear programming optimization model using the simplex method [50] has been developed, where the objective function value addressed to the capacity factor [51] of the PV park is maximized. The main restriction is that the sum of the PV power produced, the energy released in PHS and the grid supplementary power should meet load demands. Another restriction concerns the water balance in the PHS system, where the quantity of the water pumped into the upper reservoir equals the quantity of the water released downstream into the lower reservoir.
In this idealized model, no legislative restrictions are taken into account regarding the deployment of energy storage by the PHS plant. According to the National Law provision (Law 3468, Government Gazette 129/A/27 June 2006) the total amount of electricity taken in from the network on an annual basis does not exceed 30 percent of the total amount of energy consumed for the filling of the storage system in that station. Thus only 30% of the energy produced by the PV plant is allowed to be stored using the PHS plant by utilizing the potential energy of water. Another assumption considered in the optimization problem is that the load demand distribution during daytime is set directly proportional to the electrical power distribution network flow [52]. In addition, the distribution system performance of the PV park during the daytime is set directly proportional to the distribution of the solar radiation production [53] during the daytime.
As depicted in Figure 10, the operation of the HPHS system dictates a pumping procedure between 7:00 and 16:00 and energy production in the rest of the 24 h cycle. The designed energy system for the PHS plant must be credible to cover the demand load of electricity during the whole day, providing efficient energy generation and a guaranteed rate. Consequently, when the energy system supplies more power than that required by the users (i.e., the total load produced is higher than the demand load), then this additional load is used to supply the PHS plant, and thus, some additional water can be stored in the upper reservoir of the plant. However, the “real” PHS storage load is lower due to energy losses.
For this purpose, representative synthetic time series, including load prediction data, must be generated during the design phase for the optimal planning of the HPHS system, aiming at the prediction of the potential energy storage. The potential energy storage P[sub.e] in a PHS complex required to supplement the PV park can be estimated by calculating the pumping energy consumption via the following formula, which is an expansion of Bernoulli’s equation [54,55]:(23)P[sub.e]=(QH)/(102n) where P[sub.e] is the total power required (kW), Q is the flow rate (lt/s), H is the total head, i.e., the pressure head and elevation head (m), n is the combined effect of pumping efficiency and derating, and 102 is a unit conversion factor value. Pumping efficiency and derating account for the operating efficiency of pumps, and efficiency losses between the energy required at the pump shaft and the total energy required [54].
The total power required to be stored in order to return it to the electric grid when needed can be obtained using two alternatives: (a) either by having a great hydraulic difference between the upper and lower reservoir, with a relatively small quantity of water per daily cycle uplift, (b) or by having a smaller hydraulic difference between the reservoirs, though with a larger quantity of pumped water. In addition, the head difference influences the efficiency of pumps and turbines [56]. The selection between the two options is mainly a financial decision since the magnitude of the project (mainly the number of units of pumps and generators, excavation volume of both reservoirs and area covered for PHS facilities to be mounted) is proportional to the required quantity of stored water. Thus, it affects the total cost. So, the option of a significant head difference could have the advantage of minimum fixed cost.
The required quantity of stored water and the head difference set the limits for the pit lake level variation LLV, which is a critical component in the design [10]. Based on this acknowledgment, a diagram that can help to select the appropriate altitude of the PHS upper reservoir water level can be the total power required P[sub.e] per unit m of water level daily fluctuation LLV in relation to the pit lake’s level, as presented in the Section 3. This approach could help in the decision-making process and achieving a PHS plant’s economic viability.
2.5. Economic Value
The prediction of the spatiotemporal evolution of the lake is crucial for achieving the financial efficiency of the HPHS project, by evaluating different levels of the lake’s filling. The financial efficiency of each water level in the pit lake can be determined via performance indicators, i.e., net present value (NPV) and net cash recovery (NCR) of a cash flow series.
In the case of the Kardia pit lake, initial calculations of the economic feasibility of a PHS plant were performed with the help of the NPV method, where the values of alternative scenarios were evaluated. The NPV value indicates whether the overall lifetime costs of the project (benefits and costs of different scenarios), i.e., the investment, is less or more than the income generated. For example, for an investment project with a minimum design operational lifetime of up to 50 years, the following equation was used to estimate NPV [10]:(24)NPV[sub.real]=-I[sub.o]+?t=150{(R[sub.t]-C[sub.t])/q[sub.r][sup.t]} where NPV[sub.real] is the NPV in real terms, considering also the inflation rate r and the interest rate i, C[sub.t]; R[sub.t] are costs and revenues at different points in time t, I[sub.o] is the initial investment at time t = 0 and q[sub.r][sup.t] is the discounting factor which reflects the value of money in time t:(25)q[sub.r][sup.t]=1+i/1+r[sup.t]
The component q[sub.r][sup.t] in the equation above helps make the different cash flows in time comparable [10].
2.6. Geotechnical Issues to Be Considered in the Post-Mining Design
Several studies have sought to determine successful and unsuccessful cases of pit lakes in the post-closure mine phase. It has been documented that stability is improved as water levels increase inside the pit lake. Kavvadas et al. investigated the stability conditions that prevailed in decommissioning surface lignite mines in Greece [57], using the Amyntaion mine as a case study, where the pit is gradually filled with water. According to this report, it is shown that stability conditions are reinstated when the depth of water (critical depth) reaches approximately 15–35% of the final equilibrium condition, whereas, at higher lake levels, the safety factor is expected to increase significantly due to the beneficial effect of the lake water pressure acting on the slope. Table 2 presents preliminary predictions in slope stability for different scenarios of lake water development inside the open pit [57].
In addition, in the early stages of the lake’s development, more favorable slope stability conditions are expected to be achieved in only a few years. This result is associated with the fact that the surface area of the pit is smaller at deeper levels; thus, the depth of the water inside the pit increases rapidly in the early stages of the lake development.
Regarding the operation of the HPHS complex facility, in the daily cyclic process in the lower reservoir (i.e., the pit lake), excessive weathering and corrosion of the lake banks could appear, caused by the daily fluctuation of the head difference between the lower and the upper reservoir. Furthermore, the daily fluctuation induces an increase in water pressure at the pit lake banks; thus, slope instabilities could appear due to the fatigue effect. Consequently, adopting protection measures at the pit lake banks is a prerequisite, which could also affect the excavation volume and, as a result, the magnitude of the PHS complex facility.
3. Results and Discussion
3.1. Water Flux Trends and Forecast of Lake’s Spatiotemporal Development
Predictive models to forecast rainfall in the near future, as well as surface runoff, evaporation, percolation rates and consumptive irrigation rates are generated for the next decade, from which the water balance of the lake’s surface is derived. The results from the implementation of the NNAR model (neural network autoregression) with no seasonal participation are presented in Figure 11.
According to the NNAR predictions for the next decade, the maximum expected groundwater attenuation of the Sarigiol aquifer is estimated at 10 m (i.e., approximately one meter per year on average). This drawdown is due to excessive irrigation water demands [40], which allows the expectation that the groundwater level will be successively lowered below the elevation +635 m during our conceptual model of the lake filling time interval. Time series forecasts of annual drawdown for the next two decades are computed. This groundwater level attenuation was incorporated into the mathematical model for predicting the lake’s water balance, where a 20% reduction in groundwater influxes to the lake’s water body was introduced.
In addition, the lake’s surface balance is displayed in Figure 11, as a result of the difference between evaporation and precipitation over the lake’s surface. The results showed an increasing trend of yearly water losses in the lake’s surface water balance due to a gradual increase in the water surface, as evaporation also increases.
Figure 12 displays the lake’s spatiotemporal development based on the model developed in this study and the results from implementing the NNAR model. To highlight the influence of the time variable datasets, the different water inflows into the pit lake’s body have been plotted together and are associated with the time component. The results indicate that, the filling of the lake model is very sensitive to surface runoff, which is described by an average value of 0.3 for the area inside the open pit (Figure 8). There is a clear decreasing trend of annual surface runoff due to the gradual decrease in the watershed area.
As the pit lake level and volume increase with time, the surface area of the pit decreases and thus, the surface runoff decreases. It is interesting, though, that, during the last phase of the lake’s development, surface runoff still exceeds the lake’s surface losses (water balance of lake’s surface, i.e., precipitation minus evaporation on lake’s surface). Despite the large water losses (in mm/m[sup.2] of lake’s surface) compared to surface runoff (in mm/m[sup.2] of watershed surface) (Figure 11), the total annual surface runoff quantities exceed the water losses in the early stages of the lake’s development (Figure 12).
The results of the lake’s development displayed in Figure 12 indicate that the natural recharging of the lake is feasible. In the early stages, the lake expands rapidly, while in the last phase it does so very slowly. On the contrary, the lake’s water body volume increases almost steadily. As a result, the maximum natural lake growth inside the pit can potentially cover almost 50% of the pit’s total surface (Figure 2) in approximately 80 years, with a depth of approximately 117 m.
3.2. Evaluation of Potential Land Use Repurposing Scenarios
For the already closed lignite mine of Kardia, a preliminary evaluation was conducted for potential land use repurposing scenarios, by evaluating the landscape formed after the lignite exploitation activities terminated. Generally, in abandoned open pits, the pit void will eventually be filled with water at a specific time period. Thus, the viability of any repurposing land use scenarios should be evaluated concerning the spatiotemporal lake level development and the water coverage area compared to the pit’s total surface.
In evaluating different land use scenarios, the energy transition potential in the specific region is of primary importance. It is considered with respect to long-term socioeconomic viability in the area. At the same time, the transition should be in line with the holistic policy framework developed by Greece to further develop renewable technology in the electricity sector and eventually achieve net-zero greenhouse-gas emissions. Consequently, repurposing land uses includes RES technology application in the remaining post-mining landscapes. Moreover, “Smart Post-mining” has to be developed in concert with “Green and Climate Smart Mining” [58], in a holistic approach so that the climate-driven clean energy transition is managed responsibly and sustainably.
Based on topographical criteria, from the digital elevation model (DEM) and the slope model implemented (Figure 2), it can be seen that the area between elevations +650 m and +680 m could be characterized as relatively flat. The flatness of this area makes it ideal for installing PV parks. Consequently, lake growth should not exceed the elevation of +650 m. Also, based on preliminary predictions as described in Table 2, stability conditions are expected to be reinstated when the depth of lake water reaches approximately 22 m to 29 m, i.e., for pit lake elevations +545 m up to 552 m. Thus, adopting lower lake levels increases the risk of failure.
From the model developed in this study and the results from the implementation of the NNAR model, the maximum lake water level cannot exceed elevations +640 m ~ +650 m, as it is presented in Figure 12. This condition is due to hydraulic reasons, as this level is related to the groundwater level in the loose sediment aquifer overlying the lignite sequence. Considering the lake growth limitation, the flat area between elevations +650 m and +680 m could be examined for installing PV parks (Figure 13). This area is approximately 5 km[sup.2], with an average inclination of approximately 3.5%.
In addition, the area between elevations +570 m and +620 m on the lake’s northeast side could be repurposed for active (such as swimming, diving, etc.) and passive recreation (such as biking trails, recreation forest park, road access network to the lake, entertainment facilities, etc.). This area appears to have gentle slopes presenting average inclinations in the order of 12~17% and covers approximately 3 km[sup.2] of the surface.
An important issue emerging from these findings is that repurposing a PHS dam at the elevation of +545 m is feasible and can be achieved by naturally filling the pit lake for approximately three (3) years. Artificial recharging is also feasible if a fast filling is required to implement the PHS Project. More specifically, pumped water from Mavropigi and South field sumps could be diverted inside Kardia pit lake and, thus, fill the pit lake within a year.
Another alternative could be forming the lake water surface level at +570 m to implement the PHS and PV project. The findings from the implementation of the NNAR model suggest that this could be realized by naturally filling the pit lake in approximately seven (7) years. In the case of fast filling of the pit lake, this could be achieved in approximately three (3) years, i.e., four (4) years earlier.
The area on the southwest side of the lake is proposed for the upper reservoir installation of the PHS project. This area could be ideal for the location of the upper reservoir due to the proximity to both the lake (lower reservoir) and the Soulou River, where, after the PHS operational start, the required discharge of overflow is prospected to be rejected. This area covers approximately 0.5 km[sup.2] of the surface, with an average inclination of approximately 3%.
Based on these results, our findings indicate that the level + 570 m could cover the needs of a hybrid PHS plant (HPHS) and PV park while providing enough space for recreation projects on the lake banks and simultaneously ensuring slope stability. The stabilizing effect of the lake water pressure ensures more favorable slope stability conditions at higher lake levels. In addition, the stability conditions are further improved in case the PHS upper reservoir is mounted outside the mine exploitation area, at a safe distance from the excavations’ slopes.
However, water pumping is required in such cases to maintain the level of the pit lake at a certain elevation. Thus, the pit lake water quality is an important parameter that could enable new beneficial end uses, such as exploiting the overflow of the pit lake water for irrigation purposes. For example, Figure 14 shows that the pit lake water is expected to have good quality and thus could be suitable for irrigation purposes. Moreover, waters are marked by salinization risk (measured on the x-axis by conductivity) rather than alkalinization risk (measured on the y-axis by sodium).
3.3. PHS Plant Efficiency Potential
Figure 15 displays the potential energy storage P[sub.e] per unit m of the daily water level fluctuation (MW/ m of water fluctuation) in relation to the pit lake’s elevation and volume, considering two different scenarios: a) in the case of the upper reservoir mounted inside the mine at the elevation +640 m a.s.l. and b) in the case of the upper reservoir mounted outside the mine exploitation area (i.e., +670 m a.s.l.). The graph shows that in the case of the Kardia mine and for a specific pit water level at + 570 m a.s.l., the estimated potential energy storage P[sub.e] is 40% higher in the case of a 30 m increase in the total head (H). In this example, lake level variation (LLV) at the pit lake was considered 1.0 m.
In a cash flow of revenues and expenses, net present values in real terms are evaluated for the next 20 years, as presented in Figure 15. Specifically, in the case of natural recharging of the pit lake, it is evident that an increase in the net present value in the order of 34% is estimated for the pit lake elevation +545 m, compared to the scenario where the pit lake elevation is formed at +570 m. In addition, in the case of artificial recharging of the pit lake, an increase in the net present value in the order of 18% is estimated for the pit lake elevation +545 m, compared to the scenario where the pit lake elevation is formed at +570 m, due to the late recovery of investment capital. In this initial economic feasibility evaluation, constant values of inflation r=0.07 and rate of interest i=0.10 were considered, whereas the same value of the investment I[sub.o] was considered in all scenarios.
4. Conclusions
In the abandoned Kardia lignite mine, the evaluation of land potential and resources is assessed regarding the lake water development inside the open pit area. Thus, the natural development of the lake is predicted, and different repurposing land uses are considered, concerning the energy transition potential in this specific region.
Initially, the spatiotemporal development of the lake in the final pit is calculated, and subsequently, land use repurposing is evaluated with regard to lake growth. Considering the energy transition strategies at this specific mine, it is concluded that on the east side of the lake, there is a relatively flat area of approximately 5 km[sup.2], with an average inclination of approximately 3.5%. Therefore, this area is considered ideal for the installation of PV parks.
In addition, a hybrid PHS plant and PV park are examined to make provision for the energy surplus produced by the photovoltaic parks, aiming to meet electricity demand 24 h a day. Concentrating the solar power surplus with a PHS storage unit can provide flexible, renewable energy. This study has shown that repurposing a PHS dam at the elevation of + 545 m is feasible by naturally filling the pit lake for approximately three (3) years. Another scenario considering the elevation of +570 m could be achieved by recharging the pit within seven (7) years. Artificial recharging is also feasible in case fast filling is required due to the proximity of the Kardia lignite mine to two (2) other mines, the Mavropigi and South field lignite mines. Therefore, a fast filling could be achieved in less than one (1) year for a lake elevation of +545 m and less than three (3) years for an elevation of +570 m.
From an engineering point of view, forming the lake at the level of +570 m is considered a better option for implementing the HPHS plant, as the stabilizing effect of the lake water pressure ensures more favorable slope stability conditions at higher lake levels. Moreover, the slope stability conditions are further improved in case the PHS upper reservoir is mounted outside the mine exploitation area, at a safe distance from the excavations’ slopes.
Based on these results, our findings indicate that the surface lignite mine under investigation could cover both the needs of a hybrid PHS plant and PV park and, at the same time, provide enough space for activities (such as swimming, diving, etc.) and passive recreation (such as biking trails, recreation forest park, road access network to the lake, entertainment facilities, etc.) on the lake banks.
Forecasting the spatiotemporal development of a pit lake could help plan post-mining land repurposing, evaluate different scenarios, and find the optimum water level of the lake surface. The aim is to optimize land use repurposing in terms of economical design, efficiency and sustainability. In addition, the proximity of the Kardia lignite mine to the Mavropigi and South Field lignite mines is a crucial component to consider in the final design of the hybrid electric power project to make this project more economical, efficient and effective.
However, further studies need to be conducted to evaluate land repurposing, such as forecasting future water quality problems. In addition, different land repurposing projects are prospected to be developed under unfavorable geological settings of an abandoned open pit. This implies the need for detailed geotechnical (slope stability, bank erosion, corrosion, etc.) studies with the implementation of appropriate stability measures, in combination with a careful geotechnical monitoring application. To this end, continuous monitoring of ground movements is imperative. Therefore, this research question is an issue for further investigation and future model calibration when designing different repurposing land uses.
The whole system—pit lakes and alternative land repurposing scenarios—should be incorporated into a rational resources management program to optimize benefits to the land, water and electricity demands.
Author Contributions
Conceptualization, C.R. and G.L.; methodology, G.L., C.R., E.L. and E.M.; software, G.L. and E.M; validation, G.L., C.R., E.L., E.M. and G.K.; formal analysis, G.L., C.R., E.L, E.M. and G.K.; investigation, G.L., C.R., E.M. and G.K.; data curation, G.L., C.R., E.L., E.M. and G.K.; writing—original draft preparation, G.L. and E.L.; writing—review and editing, G.L., C.R. and E.M.; visualization, G.L., C.R., E.L., E.M. and G.K.; supervision G.L. and C.R.; project administration, C.R and G.L. All authors have read and agreed to the published version of the manuscript.
Data Availability Statement
Data will be made available by the authors upon request.
Conflicts of Interest
The authors declare no conflict of interest.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
1. European Commission, Publications Office of the European Union: Brussels, Belgium, 2018,. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/?uri=CELEX:52018DC0773 <date-in-citation content-type="access-date" iso-8601-date="2021-01-11">(accessed on 11 January 2021)</date-in-citation>.
2. European Commission, Publications Office of the European Union: Brussels, Belgium, 2019,. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:52019DC0640&from=EN <date-in-citation content-type="access-date" iso-8601-date="2021-01-11">(accessed on 11 January 2021)</date-in-citation>.
3. A. Sokratidou; E. Gkouvatsou; C. Roumpos; A. Perdiou; G. Tsagkarakis; S. Kaimaki Innovative Approaches to Coal Surface Mine Sites Rehabilitation: A Case Study of Megalopolis Lignite Fields.,pp. 437-452.
4. R.K. Pachauri; L. Mayer; Intergovernmental Panel on Climate Change, IPCC: Geneva, Switzerland, 2015, ISBN: 978-92-9169-143-2.. Available online: https://climate-adapt.eea.europa.eu/countries-regions/countries/greece <date-in-citation content-type="access-date" iso-8601-date="2021-01-11">(accessed on 11 January 2021)</date-in-citation>.
5. J.E.T. Bistline; G.J. Blanford Value of Technology in the U.S. Electric Power Sector: Impacts of Full Portfolios and Technological Change on the Costs of Meeting Decarbonization Goals., 2020, 86,p. 104694. DOI: https://doi.org/10.1016/j.eneco.2020.104694.
6. F. Pavloudakis; C. Roumpos; E. Karlopoulos; N. Koukouzas Sustainable Rehabilitation of Surface Coal Mining Areas: The Case of Greek Lignite Mines., 2020, 13, 3995. DOI: https://doi.org/10.3390/en13153995.
7. K. Bódis; I. Kougias; N. Taylor; A. Jäger-Waldau Solar Photovoltaic Electricity Generation: A Lifeline for the European Coal Regions in Transition., 2019, 11, 3703. DOI: https://doi.org/10.3390/su11133703.
8. Z. Nikolic; D. Nikolic Practical Example of Solar and Hydro Energy Hybrid System—The Need for a Reversible Power Plant., Volume 6,pp. 165-171.
9. J. Keenan; S. Holcombe Mining as a Temporary Land Use: A Global Stocktake of Post-Mining Transitions and Repurposing., 2021, 8,p. 100924. DOI: https://doi.org/10.1016/j.exis.2021.100924.
10. M. Wessel; R. Madlener; C. Hilgers Economic Feasibility of Semi-Underground Pumped Storage Hydropower Plants in Open-Pit Mines., 2020, 13, 4178. DOI: https://doi.org/10.3390/en13164178.
11. C. McCullough; M. Schultze; J. Vandenberg Realizing Beneficial End Uses from Abandoned Pit Lakes., 2020, 10, 133. DOI: https://doi.org/10.3390/min10020133.
12. S. Holcombe; J. Keenan, Centre for Social Responsibility in Mining, Sustainable Minerals Institute: Brisbane, Australia, 2020,. The Social Aspects of Mine Closure Research Consortium <date-in-citation content-type="access-date" iso-8601-date="2021-01-13">(accessed on 13 January 2021)</date-in-citation>.
13. C.H. Gammons; L.N. Harris; J.M. Castro; P.A. Cott; B.W. Hanna Creating Lakes from Open Pit Mines: Processes and Considerations, Emphasis on Northern Environments., 2009,p. 106.
14. M. Schultze The Filling and Remediation of Pit Lakes in Former Open Cast Lignite Mines., Tech Univ Braunschweig Fac Archit Constr Eng Environ Sci: Braunschweig, Germany, 2012,
15. C.D. McCullough Rehabilitation of Lake Kepwari: A Previously Acidic Mine Lake in Western Australia., Cardiff University: Cardiff, Wales, UK, 2021,pp. 346-351.. Available online: https://www.imwa.info/imwaconferencesandcongresses/proceedings/325-proceedings-2021.html <date-in-citation content-type="access-date" iso-8601-date="2022-01-20">(accessed on 20 January 2022)</date-in-citation>.
16. M. Lund; E. van Etten; J. Polifka; M.Q. Vasquez; R. Ramessur; D. Yangzom; M.L. Blanchette The Importance of Catchments to Mine-Pit Lakes: Implications for Closure., 2020, 39,pp. 572-588. DOI: https://doi.org/10.1007/s10230-020-00704-8.
17. M.A. Lund; M. Blanchette Can Saline Pit Lakes Offer Biodiversity Values at Closure?., Cardiff University: Cardiff, Wales, UK, 2021,pp. 300-306.. Available online: https://www.imwa.info/imwaconferencesandcongresses/proceedings/325-proceedings-2021.html <date-in-citation content-type="access-date" iso-8601-date="2022-01-20">(accessed on 20 January 2022)</date-in-citation>.
18. M.L. Blanchette; M.A. Lund Pit Lakes Are a Global Legacy of Mining: An Integrated Approach to Achieving Sustainable Ecosystems and Value for Communities., 2016, 23,pp. 28-34. DOI: https://doi.org/10.1016/j.cosust.2016.11.012.
19. D. Snashall; K. Cox Woodlawn Advanced Energy Recovery Centre Scoping Report; Australia.. Available online: https://majorprojects.planningportal.nsw.gov.au <date-in-citation content-type="access-date" iso-8601-date="2021-05-11">(accessed on 11 May 2021)</date-in-citation>.
20. S. Pytel; S. Sitek; M. Chmielewska; E. Zuzanska-Zysko; A. Runge; J. Markiewicz-Patkowska Transformation Directions of Brownfields: The Case of the Górnoslasko-Zaglebiowska Metropolis., 2021, 13, 2075. DOI: https://doi.org/10.3390/su13042075.
21. T. Ma; H. Yang; L. Lu; J. Peng Technical Feasibility Study on a Standalone Hybrid Solar-Wind System with Pumped Hydro Storage for a Remote Island in Hong Kong., 2014, 69,pp. 7-15. DOI: https://doi.org/10.1016/j.renene.2014.03.028.
22. C.R. Matos; J.F. Carneiro; P.P. Silva Overview of Large-Scale Underground Energy Storage Technologies for Integration of Renewable Energies and Criteria for Reservoir Identification., 2019, 21,pp. 241-258. DOI: https://doi.org/10.1016/j.est.2018.11.023.
23. M.J. Kavvadas; P.G. Marinos Prediction of Groundwater Table Lowering for Lignite Open-Cast Mining in a Karstic Terrain in Western Macedonia, Greece., 1994, 27,pp. S41-S55. DOI: https://doi.org/10.1144/GSL.QJEGH.1993.027.0S.06.
24. G. Louloudis Hydrogeological conditions of South Lignite Bearing Field of Ptolemais region. Confrontation of groundwater problems during the exploitation., Natl Tech University Athen Sch Min Metall Eng: Athens, Greece, 1991,. (In Greek)
25. D. Dimitrakopoulos; E. Grigorakou Aquatic Balance of Vegoritis Lake, West Macedonia, Greece, Ralated to Mining Activity., University of Newcastle: Newcastle upon Tyne, UK, 2004,pp. 201-207.. Available online: https://www.imwa.info/imwaconferencesandcongresses/proceedings/189-proceedings-2004.html <date-in-citation content-type="access-date" iso-8601-date="2022-01-20">(accessed on 20 January 2022)</date-in-citation>.
26. N. Koukouzas; C.R. Ward; Z. Li Mineralogy of Lignites and Associated Strata in the Mavropigi Field of the Ptolemais Basin, Northern Greece., 2010, 81,pp. 182-190. DOI: https://doi.org/10.1016/j.coal.2009.12.012.
27. M. Kavvadas; C. Roumpos; P. Schilizzi Stability of Deep Excavation Slopes in Continuous Surface Lignite Mining Systems., 2020, 38,pp. 791-812. DOI: https://doi.org/10.1007/s10706-019-01066-x.
28. H.E. Beck; N.E. Zimmermann; T.R. McVicar; N. Vergopolan; A. Berg; E.F. Wood Present and Future Köppen-Geiger Climate Classification Maps at 1-Km Resolution., 2018, 5,p. 180214. DOI: https://doi.org/10.1038/sdata.2018.214. PMID: https://www.ncbi.nlm.nih.gov/pubmed/30375988.
29. C.W. Thornthwaite; J.R. Mather, Drexel Institute of Technology, Laboratory of Climatology, N.J. Publications in Climatology: 1955, Volume VIII,p. 104.. Available online: https://www.jstor.org/stable/210739?origin=crossref <date-in-citation content-type="access-date" iso-8601-date="2022-01-20">(accessed on 20 January 2022)</date-in-citation>.
30. W. Kinzelbach, Elsevier: Amsterdam, The Netherlands, 1986, ISBN: 978-0-444-42582-9.. Developments in water science
31. J.C. Calvo An Evaluation of Thornthwaite’s Water Balance Technique in Predicting Stream Runoff in Costa Rica., 1986, 31,pp. 51-60. DOI: https://doi.org/10.1080/02626668609491027.
32. R.V. Blau; P. Hoehn; P. Hufschmied; A. Werner Ermittlung Der Grundwasserneubildung Aus Niederschlägen., 1983, 63,pp. 45-54.
33. G. Louloudis; E. Louloudis; C. Roumpos; E. Mertiri; G. Kasfikis; K. Chatzopoulos Forecasting Development of Mine Pit Lake Water Surface Levels Based on Time Series Analysis and Neural Networks., 2021, 41,pp. 458-474. DOI: https://doi.org/10.1007/s10230-021-00844-5.
34. J.C. Manning, 3rd ed. edition; Prentice Hall: Hoboken, NJ, USA, 1997,
35. J. Dupuit, Dunod: Paris, Frence, 1863,. (In French)
36. P. Krassakis; K. Pyrgaki; V. Gemeni; C. Roumpos; G. Louloudis; N. Koukouzas GIS-Based Subsurface Analysis and 3D Geological Modeling as a Tool for Combined Conventional Mining and In-Situ Coal Conversion: The Case of Kardia Lignite Mine, Western Greece., 2022, 2,pp. 297-314. DOI: https://doi.org/10.3390/mining2020016.
37. G. Zhan; M. Lin Shelp, The University of Western Australia: Perth, Australia, 2009,pp. 461-471.
38. M. Tasumi; R.G. Allen; R. Trezza Calibrating Satellite-Based Vegetation Indices to Estimate Evapotranspiration and Crop Coefficients., U.S. Committee on Irrigation and Drainage: Denver, CO, USA, 2006,pp. 103-112.
39. N. Memmos, Ministry of Agriculture, Sector of Land Reclamation Service (YEB, Acronym in Greek): Thessaloniki, Greece, 1973,. unpublished
40. G. Louloudis; N. Stathopoulos Irrigation Water Consumptive Use Changes in South Ptolemais Lignite Bearing Hydrogeological Basin (Greece) Using NDVI Remotely Sensed Data., Harokopio University of Athens: Athens, Greece, 2019, Volume 7,pp. 611-612.. Bulletin of the Geological Society of Greece <date-in-citation content-type="access-date" iso-8601-date="2022-01-20">(accessed on 20 January 2022)</date-in-citation>.
41. P. Panagos; M. Van Liedekerke; A. Jones; L. Montanarella European Soil Data Centre: Response to European Policy Support and Public Data Requirements., 2012, 29,pp. 329-338. DOI: https://doi.org/10.1016/j.landusepol.2011.07.003.
42. R.H. McCuen; T.R. Bondelid Relation Between Curve Number and Runoff Coefficient., 1981, 107,pp. 395-400. DOI: https://doi.org/10.1061/JRCEA4.0001425.
43. N. Kim; M.-J. Shin Estimation of Peak Flow in Ungauged Catchments Using the Relationship between Runoff Coefficient and Curve Number., 2018, 10, 1669. DOI: https://doi.org/10.3390/w10111669.
44. M. Giandotti Previsione delle piene e delle magre dei corsi d’acqua., Ministero dei LL.PP., Servizio Idrografico Italiano: Rome, Italy, 1934, Volume 8,pp. 107-117.
45. H.L. Penman Natural Evaporation from Open Water, Bare Soil and Grass., 1948, 193,pp. 120-145. DOI: https://doi.org/10.1098/rspa.1948.0037.
46. W.D. Sellers, The University of Chicago Press: Chicago, IL, USA, 1965,
47. O. Tetens About some meteorological terms., 1930, 6,pp. 297-309.. (In German)
48. J.R. Eagleman Pan Evaporation, Potential and Actual Evapotranspiration., 1967, 6,pp. 482-488. DOI: https://doi.org/10.1175/1520-0450(1967)006<0482:PEPAAE>2.0.CO;2.
49. J. Monteith; M. Unsworth, Academic Press: Warszawa, Poland, 2013, Volume 4, ISBN: 978-0-12-386910-4.
50. S.I. Gass Linear Programming., John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2006,p. ess1466.pub2. ISBN: 978-0-471-66719-3.
51. A. Boretti; S. Castelletto; W. Al-Kouz; J. Nayfeh Capacity Factors of Solar Photovoltaic Energy Facilities in California, Annual Mean and Variability., 2020, 181,p. 02004. DOI: https://doi.org/10.1051/e3sconf/202018102004.
52. H. Tyralis; G. Karakatsanis; K. Tzouka; N. Mamassis Data and Code for the Exploratory Data Analysis of the Electrical Energy Demand in the Time Domain in Greece., 2017, 13,pp. 700-702. DOI: https://doi.org/10.1016/j.dib.2017.06.033.
53. G. Koudouris; P. Dimitriadis; T. Iliopoulou; N. Mamassis; D. Koutsoyiannis Investigation on the Stochastic Nature of the Solar Radiation Process., 2017, 125,pp. 398-404. DOI: https://doi.org/10.1016/j.egypro.2017.08.076.
54. K. Faour, NSW Agriculture: 2001,. Available online: https://lib.k-state.edu/gracroplands/cite.php?CiteID=3758 <date-in-citation content-type="access-date" iso-8601-date="2022-01-20">(accessed on 20 January 2022)</date-in-citation>.
55. A. Gemitzi; V. Lakshmi Estimating Groundwater Abstractions at the Aquifer Scale Using GRACE Observations., 2018, 8, 419. DOI: https://doi.org/10.3390/geosciences8110419.
56. E. Pujades; P. Orban; S. Bodeux; P. Archambeau; S. Erpicum; A. Dassargues Underground Pumped Storage Hydropower Plants Using Open Pit Mines: How Do Groundwater Exchanges Influence the Efficiency?., 2017, 190,pp. 135-146. DOI: https://doi.org/10.1016/j.apenergy.2016.12.093.
57. M. Kavvadas; C. Roumpos; A. Servou; N. Paraskevis Geotechnical Issues in Decommissioning Surface Lignite Mines—The Case of Amyntaion Mine in Greece., 2022, 2,pp. 278-296. DOI: https://doi.org/10.3390/mining2020015.
58. I.M. Jiskani; Q. Cai; W. Zhou; S.A. Ali Shah Green and Climate-Smart Mining: A Framework to Analyze Open-Pit Mines for Cleaner Mineral Production., 2021, 71,p. 102007. DOI: https://doi.org/10.1016/j.resourpol.2021.102007.
59. L.A. Richards; L.E. Allison; J.W. Brown; H.E. Hayward; L.A. Richards; L. Bernstein; M. Fireman; G.A. Pearson; L.V. Wilcox; C.A. Bower et al., Salinity Laboratory, Soil and Water Conservation Research Branch, Agricultural Research Service, US Department of Agriculture: Washington, DC, USA, 1954,pp. 69-82.
60. L.V. Wilcox, Salinity Laboratory, US Department of Agriculture: Washington, DC, USA, 1955,pp. 1-19.
61. L.V. Wilcox, Salinity Laboratory, US Department of Agriculture: Washington, DC, USA, 1948,pp. 1-40. DOI: https://doi.org/10.22004/AG.ECON.170282.
Figures and Tables
Figure 1: Research area and location of the Kardia Lignite Mine (LM), Greece. Basemap obtained from https://www.esri.com/en-us/home accessed on 30 March 2022 (Coordinate Reference System GGRS87/Greek Grid). [Please download the PDF to view the image]
Figure 2: Kardia lignite mine geomorphology, Greece (Coordinate Reference System GGRS87/Greek Grid): (a) digital elevation model (DEM) of the open pit; (b) slope model (in degrees). [Please download the PDF to view the image]
Figure 3: Geological and hydro-lithological map in the Ptolemais lignite field, northern Greece (Coordinate Reference System GGRS87/Greek Grid). [Please download the PDF to view the image]
Figure 4: Flowchart of the methodological approach. [Please download the PDF to view the image]
Figure 5: Mean annual: (a) rainfall and (b) temperature from the Kozani weather station (blue color) and Ptolemaida weather station (green color), northern Greece. [Please download the PDF to view the image]
Figure 6: Sensitivity analysis of percolation rate SL, as a percentage of precipitation N, towards field capacity FC at the wider area of Ptolemais lignite field, northern Greece. [Please download the PDF to view the image]
Figure 7: Hydrogeological conceptual model for Kardia pit lake, Greece (the location of the cross section A-A’ is shown in Figure 3). [Please download the PDF to view the image]
Figure 8: Spatial distribution of: (a) runoff coefficient by processing CN European raster; (b) spectra indices of NDVI using Landsat-8 (image acquired in July 2018); (c) irrigation-specific consumptive use SCU (m[sup.3]/m[sup.2]); (d) annual drawdown in m, corresponding to the deficit of groundwater recharge due to excessive irrigation consumptive use of groundwater in Kardia lignite field, northern Greece (Coordinate Reference System GGRS87/Greek Grid). [Please download the PDF to view the image]
Figure 9: (a) Irrigation-specific consumptive use (m[sup.3]/m[sup.2]) versus NDVI using Landsat-8 (image acquired in July 2018) and (b) Percolation, irrigation and groundwater deficit time series during the period 2006–2021 in Sarigiol basin, Greece. [Please download the PDF to view the image]
Figure 10: Representative scheme for the hybrid pumped hydro storage (HPHS) plant operational schedule during daytime hours and photovoltaic (PV) integration, considering the installation of a 270 MW PV plant in the Kardia mine. [Please download the PDF to view the image]
Figure 11: Artificial neural network prediction (NNAR) implemented for the Kardia lignite field area, northern Greece. NNAR results of the annual: (a) rainfall; (b) runoff; (c) water balance of lake’s surface (precipitation minus evaporation on lake’s surface); (d) percolation rate; (e) specific consumptive irrigation use; (f) drawdown (percolation minus groundwater attenuation as a result of irrigation consumption); (g) groundwater level in the Sarigiol aquifer (i.e., above the lignite sequence). [Please download the PDF to view the image]
Figure 12: Forecast of the spatiotemporal development of the Kardia pit lake: (a) lake surface vs surface runoff decrease; (b) annual inflows and outflows; (c) the absolute elevations and gradual increase in stored water volume vs. time; (d) the absolute elevations and gradual increase in stored water volume vs. lake depth. [Please download the PDF to view the image]
Figure 13: Land use repurposing inside the Kardia mine, Greece (Project Coordinate Reference System GGRS87/Greek Grid). [Please download the PDF to view the image]
Figure 14: Classification of groundwater and the mine’s internal sump water for irrigation purposes: (a) diagram of sodium hazard (sodium adsorption ration) and salinity hazard (electrical conductivity) [59,60]; (b) suitability of water for irrigation use, according to Wilcox diagram [61]. [Please download the PDF to view the image]
Figure 15: (a) Potential power of storage provision (MW/ m of water fluctuation) in relation to pit lake volume and water level, considering 1.0 m for the pit lake level variation (LLV), in the case of the upper reservoir mounted inside the mine (+ 640 m a.s.l.) and outside the mine exploitation area (+ 670 m a.s.l.) and, (b) indicative net present value in real terms (NPV[sub.real]) in alternative scenarios of the pit lake’s elevation and filling procedure, inside the Kardia mine, Greece. [Please download the PDF to view the image]
Table 1: Mean annual water balance for the unsaturated zone of the Sarigiol loose sediment basin.
Hydrological Cycle in Sarigiol Basin(2006–2021) | mm/yr | 10[sup.6] m[sup.3]/yr | Precipitation (%) |
---|---|---|---|
Mean annual precipitation | 584 | 123 | 100 |
Mean annual real evaporation | 425 | 90 | 73 |
Mean annual percolation | 86 | 18 | 15 |
Mean annual runoff | 72 | 15 | 12 |
Table 2: Variation in slope stability conditions related to the pit lake development in decommissioning surface lignite mines in Greece, used for preliminary prediction.
Lake Water Depth (Hw)/Height of the Slope (H) | FS/FSo* |
---|---|
1.0 | |
0.1 | <1.0 |
0.2 | 1.0 |
0.4 | 1.075 |
0.6 | 1.25 |
0.8 | 1.6 |
Note(s): * FS = Safety factor, FSo = initial FS when water depth Hw = 0. Source: [57].
Author Affiliation(s):
[1] Public Power Corporation of Greece, Department of Mining Engineering and Closure Planning, 10432 Athens, Greece
[2] Department of Statistics, Athens University of Economics and Business, 10434 Athens, Greece
Author Note(s):
[*] Correspondence: [emailprotected] (G.L.); [emailprotected] (C.R.); Tel.: +30-697-000-5090 (G.L.); +30-697-979-9291 (C.R.)
DOI: 10.3390/w14213558
COPYRIGHT 2022 MDPI AG
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2024 Gale, Cengage Learning. All rights reserved.