Rational and Sustainable Water Resource Management in the Ptolemais Lignite Basin Using Remotely Sensed Data †

: Future investment feasibility studies concerning post-mining repurposing utilities and economic transitions should focus on regional water resource management and the hydraulic protection of any utilities. Satellite images in different bands and Digital Elevation Models (DEM) of the Ptolemais basin were processed, leading to a more accurate estimation of the runoff ratio and percolation ratio. Furthermore, the saturated and unsaturated areas were delineated, leading to the recognition of potential artiﬁcial ground water recharge zones and zones where appropriate hydraulic protection measures are necessary.


Introduction
Greece is committed to the European Commission's long-term strategic vision to achieve technological neutrality to CO 2 emissions by using renewable energies [1,2]. In the era of the early and rapid closure of mines, the Public Power Corporation (PPC) of Greece is engaged in planning for the transition away from coal and the adoption of Renewable Energy Source (RES) technology [3]. As a result, multi-level planning is critical to design and implement effective strategies, explore repurposing options, and make decisions about a post-mining future and economic transitions. Planning should include (a) water management issues, examining the new demands of water, (b) the spatiotemporal development of water inside the mined-out pits, including the effect of the planned end uses, (c) the hydraulic protection of the new RES and other projects (e.g., the application of pumped hydro storage (PHS) technology in the abandoned mines), and (e) the needs in water for successful end uses.
During the operation period in lignite mines, dewatering measures were adopted (dewatering wells, pumping stations) to lower the groundwater table by pumping water from the aquifer and abstracting rainwater accumulated in the pits. However, RES and other projects in the mine complex area will demand new water management strategies.
The main objective of this research is to provide a methodology to identify areas where the application of hydraulic protection is imperative, in order to estimate the design water inflow risk and the feasibility of hydraulic protection works via remotely sensed data from satellites (satellite images, DEM, etc.).

Research Area-Hydrogeological Conditions
in the Ptolemais basin include Tertiary and Quaternary sediments up to 1000 m thick, overlying Paleozoic schists, granites and ophiolites [3,5]. Furthermore, the Pelagonian Structural Zone-consisting of Mesozoic dolomitic limestones, volcanic sediments and flysch-overlays the Paleozoic formations. The encountered sediments consist of sand, clay, marls, conglomerates, and lignite on the surface.
Mountain stream networks consist of bedrock and mainly consolidated sediment channels, where the water flows downhill and eventually into the Soulou River. The waterbearing horizons that lie above the impermeable lignite strata form a unique hydraulically connected network that behaves as one aquifer [6]. The surface runoff of the mountainous surroundings to the hydrogeological basin has been proved insignificant according to Papakonstantinou's [7] measurements of Vermion stream discharge. However, a considerable volume of hydraulic studies conducted by the PPC [8] has considered high flooding rates to define the necessary hydraulic protection measures. Concerning the contradictory results of different hydrogeological or topographical approaches, an interesting research question concerns the investigation of the flooding hazards based on remotely sensed data.
The mean temperature and rainfall calculated for the last 14 years are 14.22 • C and 512.29 mm, with standard deviation values of 0.45 • C and 129.51 mm for the Ptolemaida weather stations, respectively. According to the Köppen climate classification, the climate is characterised as moist with mild winters, wet for all seasons, with a long and hot summer.

Research Area-Hydrogeological Conditions
The research area refers to the South Ptolemais lignite basin (Western Macedonia, Greece), where open-cast mining [4] was developed (Figure 1). The formations encountered in the Ptolemais basin include Tertiary and Quaternary sediments up to 1000 m thick, overlying Paleozoic schists, granites and ophiolites [3,5]. Furthermore, the Pelagonian Structural Zone-consisting of Mesozoic dolomitic limestones, volcanic sediments and flysch-overlays the Paleozoic formations. The encountered sediments consist of sand, clay, marls, conglomerates, and lignite on the surface.
Mountain stream networks consist of bedrock and mainly consolidated sediment channels, where the water flows downhill and eventually into the Soulou River. The water-bearing horizons that lie above the impermeable lignite strata form a unique hydraulically connected network that behaves as one aquifer [6]. The surface runoff of the mountainous surroundings to the hydrogeological basin has been proved insignificant according to Papakonstantinou's [7] measurements of Vermion stream discharge. However, a considerable volume of hydraulic studies conducted by the PPC [8] has considered high flooding rates to define the necessary hydraulic protection measures. Concerning the contradictory results of different hydrogeological or topographical approaches, an interesting research question concerns the investigation of the flooding hazards based on remotely sensed data.
The mean temperature and rainfall calculated for the last 14 years are 14.22 °C and 512.29 mm, with standard deviation values of 0.45 °C and 129.51 mm for the Ptolemaida weather stations, respectively. According to the Köppen climate classification, the climate is characterised as moist with mild winters, wet for all seasons, with a long and hot summer.

Methods
A wide range of techniques were used, and satellite images in different bands and Digital Elevation Models (DEM) of the Ptolemais basin were processed and analysed. These include the use of topography indices and satellite image processing indices.
Topography has a major impact on the catchment geomorphology and hydrology. DEMs were processed, and the spatial distribution of the topographic attributes was used to describe the spatial soil moisture patterns. The main topographic indices that can be deduced by processing DEMs are the Topographic Wetness Index (TW I) [9,10], Power Stream Index (SI) [10] and Slope length and steepness factor (LS-factor) [11,12]. Equations (1)-(3) were applied to calculate the topographic indices: where As is the specific catchment's area, b is the local slope gradient, and a is the total cumulative upslope area drained through a point per unit contour length. Satellite image processing methods can be a useful tool to map water features and monitor land cover. Landsat spectral bands with a high spatial resolution were used to generate output layers and calculate the spectral indices, as a product of the bands' Top of the Atmosphere (TOA) reflectance value. Landsat-5 and Landsat-8 satellite images were used in the analyses performed. The main indices that can be deduced by processing satellite images are the Normalised Difference Vegetation Index (NDV I) [13,14], Normalized Difference Water Index (NDW I) [14][15][16][17], and Modified Normalized Difference Water Index (MNDW I) [14,18,19]. Equations (4)-(6) were applied to calculate the spectral indices: where N IR is the TOA reflectance value of the near-infrared band, RED is the TOA reflectance value of the red band, GREEN is the TOA reflectance value of the green band, and MIR is a middle (shortwave) infrared band.
Regarding the satellite images, these were acquired in October 2016, in December 2016, in October 2018 and in October 2019, and were compared in order to obtain more accurate information. More specifically, the Sentinel-2 Satellite Imagery was used to calculate NDW I and MNDW I, while for the NDV I analysis, two satellite images on 26 July 2007 (Landsat-5) and on 15 July 2018 (Landsat-8) were used.
Maps were produced for each of the indices, and also for all of the satellite images mentioned above. Through the spatial distribution computation of the indices, a better presentation of the hydrological status was achieved. Furthermore, the water content origin and flow in the unsaturated and saturated zones were clearly delineated.
Finally, the runoff coefficient (Cr) [20,21] was estimated using Equations (7) and (8): where P vol is the volume of rainfall in mm, CN is the curve number, I Tc is the mean rainfall intensity for the T C (mm/h), and T C is the time of the concentration (h), which is usually described by the Giandotti equation [22].

Results and Discussion
The study area is dominated by small bodies of water and streams of low capacity (Figure 1). These water bodies were identified in the maps produced using the topographic and spectral indices. Regarding the appearance and location of these water bodies, the results are consistent with the former studies and topographical mapping conducted by the Hellenic Military Geographical Service. Indicative maps of TW I and PI are presented in Figure 2, and maps of NDW I are presented in Figure 3.
The drainage channels developed through the TW I index are identical to the hydrographic network recorded by geologists through field research (Figure 1). The saturated areas were identified through the spatial distribution of the SI, TW I and LS indices, where the relative wetness in the catchment areas is predicted. From the evaluation, it is evident that the potential for the streams to be hydrologically active is higher for Vermion streams compared to Askion streams. Soulou river and its channels appear to be more active, presenting perennial flow. Furthermore, potential water bodies with an intermittent or ephemeral flow were identified at the N-E part of the South Field's external waste dump margins and the internal waste deposition at the west margins inside Kardia mine. A significant flow velocity was evident only in the consolidated rocky mountains (Figure 2b). The same applies for the spatial distribution of the LS values, indicating that water discharge is apparent only on the Vermion, Askion and Skopos mountains, while the intensity of weathering is more evident on Vermion Mountain.
The water-stressed conditions were subsequently evaluated via spectral indices. The NDW I spatial extent (Figure 3) shows water accumulation in specific areas, i.e., the old Sarigiol swamp (the area south of the exploitable lignite deposits), the external waste dump (north of the exploitable lignite deposits) and the mine sumps. The most striking observation to emerge from the data analysis was the different appearance of the riverlets between Askion Mountain, Vermion Mountain and Sarigiol basin (Skopos Mountain). The results indicate that the streams of Vermion Mountain can be characterised as wet areas over all of winter, despite their intermittent discharge.
Finally, by processing the CN European raster file [23], the runoff coefficient was calculated ( Figure 4). Regarding excavations and inside dumping, the present study showed that the runoff coefficient is much higher (20-30%) than that outside the pit perimeter, where the runoff coefficient hardly exceeds 5%. Only south of South Lignite Field does the runoff coefficient reach 10%. However, the surface flow direction is the opposite of the direction towards the exploitation front, implying that hydraulic protection measures are unnecessary.
Furthermore, the results strongly imply that all of the streams coming from Vermion Mountain have a low runoff coefficient (less than 10%) in the region of the talus cones and conglomerates, except for the Heimarros and Ksiropotamos streams, where the runoff coefficient is higher in places running into loose sediments. Similarly, all of the streams coming from Askion Mountain have a low runoff coefficient (less than 10%) in the region of the talus cones and conglomerates, except for the Ligotypos and Lagorema streams, where the runoff coefficient is higher at places outside the lignite fields, presenting values in the order of 10-20%.
Taking all of the findings together, areas where the application of hydraulic protection was imperative, were identified ( Figure 4b). Subsequently, hydraulic protection measures were defined and dimensioned. Finally, the design peak flow was estimated using the rational method [24], while for the mean rainfall intensity, the approved Flood Risk Management Plan of the West Macedonia River Basin District was adopted (EL 09) [25]. Indicative computations are presented in Table 1, and are compared to the calculations performed based on previous hydraulic studies' assumptions of the surface runoff value for the same area. The results indicate a reduction in the critical depth of flow by approximately 40% compared to previous hydraulic studies.

Conclusions
The changing environment in the South Ptolemais basin-with active mining areas, the rapid closure of coal mines in progress, rehabilitation plans and post-mining repurposing scenarios-is rather challenging to inspect without modern monitoring techniques. Satellite images and DEM processing provide a powerful tool for hydraulic protection and hazard assessment. In this research, several effective topography indices and satellite image processing indices were selected in order to examine the hydrological condition of the area under investigation.
The results from the analyses indicated the areas where the application of hydraulic protection was imperative. Subsequently, using remotely sensed data from satellites, a more reasonable runoff coefficient was estimated to dimension the open channels and diversion ditches. Thus, the findings of this research contribute to the growing evidence that the use of remotely sensed data could become an effective technique for the delineation of saturated and unsaturated areas, and the facilitation of decisions regarding hydraulic protection, especially in large areas with changing environments.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.