In total, 76 mining pit lakes of more than 30 m from the major axis have been located in four mining zones of the León province: Fabero (28), Toreno (18), Villablino (17) and Valderrueda (13). All the basic morphometric parameters and location of each of the pit lakes can be consulted on the interactive map at the following link: https://rawcdn.githack.com/amelon00/Coal_Pit_Lakes_Leon/54fdf2dbf04c9fff6c20e4ef8d1f89ffd2682ef7/CoalPitLakesMap_v08072021.html
Despite the diversity suggested by the morphometric data, can be grouped into three categories according to their topographic location of the cuts that originated them (Fig. 3): hillside (52), valley bottom (16) and hilltop (8). In all the basins, pit lakes on hillsides predominate (Fig. 4), while those located in valley bottoms are found mainly in the Valderrueda basin.
Figure 3. Examples of coal mining cuts and associated pit lakes. A) hilltop (Igüeña); B) hillside (San Justo Cabanillas); C) valley bottom (Cabrillanes)
Regarding the shape of the pit lakes according to the DL index, the circular (39) and subcircular (30) shapes predominate, followed by the elongated subrectangular (7) shape, with none reaching the dendritic category. If we look at the distribution of frequencies (Fig. 4) in all basins, circular basins predominate except in Toreno, where subcircular basins are predominant.
Figure 4. Distribution of the qualitative characteristics (shape, type of cut, age and connectivity) of the pit lakes in the four coal mining zones where the open-pit mines are located.
Regarding the age of the pit lakes (Fig. 4), there are 19 prior to 1987, 28 between 1987–1997, which is the period in which 37% of the total are created, 12 between 1997–2007, and 17 after 2007. The results show an unequal distribution according to the zone considered: in Villablino, they do not exist before 1987, while in Valderrueda, there are none after 1997.
The endorheic nature of the pit lake depends on its connectivity with the drainage network of the zone in which it is located. In general, they have low connectivity because there is a high percentage of pit lakes (88% in Toreno and 78% in Fabero) with no outflow that are confined in their basins.
Regarding the morphometric parameters according to the cumulative frequency distribution for length and width, small-sized pit lakes with a major axis less than 100 m in length predominate in all mining zones, except in Villablino, where most are between 100–200 m (Fig. 5), and some less than 50 m wide.
Figure 5. Cumulative frequency histogram of the morphometric and pH parameters of the pit lakes according to their distribution by mining zone.
Overall, small- to medium-sized pit lakes predominate since most of their surfaces are less than 5000 m2 (in Valderrueda and Toreno, almost 80% of the total; in Fabero, more than 60%), while in Villablino, they are slightly more than 50% and shallow (depths < 5 m). Consequently, the volume of most is less than 25000 m2 (Fig. 5).
The pH of the water in the pit lakes is predominantly between 6.5 and 7.5, almost always far from the normal or neutral value of approximately 7 (Ayala Carcedo and Vadillo Fernández 1989). However, pH values even further from neutral value were also obtained. The lowest values recorded (water with the highest acidity) are 3.59 (Fabero-05), 3.97 (Fabero-06) and 4.43 (Fabero-14). The most alkaline waters are 8.86 (Villablino-13), 8.42 (Villablino-01) and 8.45 (Valderrueda-04). On the other hand, the frequency distribution of alkalinity and acidity indicates differences according to the mining zone considered (Fig. 5). Thus, some zones are clearly above pH 7 with alkaline values (Valderrueda, Villablino and Toreno), while in Fabero, most values (over 60%) are below 7, that is, acidic waters.
PCA was applied using the 10 variables for each pit lake, differentiated into two categories (by shape and by age): eight variables referring to morphological and morphometric parameters (surface, perimeter, maximum length, maximum width, maximum and mean depth, volume, coastline development) and two environmental variables (open-pit lake orientation and water pH).
The correlations between numerical variables used in the PCA are shown in Fig. 6. The correlations are moderate and high, and positive, between the morphometric parameters indicating the morphology of the pit lakes, except when these are related to the coastline development index, which is an indicator of the irregularity of the pit lake contour. The variables of water pH and orientation also show low or no correlations with the rest, which indicates that these three variables do not depend on the morphometry of the pit lakes.
Figure 6. Correlation matrix of the variables used in the PCA.
PCA allowed the extraction of two main components that explained 65.78% of the variance (Fig. 7). The first principal component explains 51.2% of the variance and is related to strictly morphometric variables (mean and maximum depth, width, perimeter, length and area). The second component has a lower weight on the variables, explaining ~ 14.6% of the variance. In the same way as in the correlation analysis, different relationships between groups of variables are observed: on the one hand, there are two groups of variables with relationships between them, such as the variables of depth and volume, and on the other, the variables related to the dimensions of the pit lakes (surface, length, width and perimeter). The variables of coastline development, pH and orientation are less related, both with each other and with the rest of the variables.
Figure 7. Variables factor map (PCA). Grey labelled variables are quantitative supplementary variables (age, shape, connectivity and location).
Figure 8 is a map of individuals (pit lakes) factor map, showing some of the pit lakes with higher contribution to the plane. Dashed lines indicate the mean values for each dimension. Individuals are distinguished according to shapes (8A) and ages (8B) (the position of the individuals does not change in both graphs). In Fig. 8, the pit-lakes with circular shapes are to the right of the first factorial plane (dimension 1). In the upper plane (Fig. 8A), the deepest pit-lakes stand out: Valderrueda-06 and 03, formed before 1987, of large dimensions, with alkaline pH, and those of Villablino-02 and Fabero-09, with similar morphometric values but more recently formed (1987-97), all located in the valley bottom. In the lower plane (Fig. 8B) are the largest pit lakes, Fabero-07 and Fabero-11, and more recently, after 2007 (also Fabero-04 and 05), with soft acid waters (Ayala Carcedo and Vadillo Fernández 1989) located on the hillside. All of them share, in addition, large width, perimeter and length values. The Fabero 11 pit is striking, the most extensive of all those analysed, which is subcircular in shape but shows characteristics that are more typical of circular pit lakes. It was created after 2007, and its waters show an alkaline pH.
Figure 8. Representation of individuals (pit lakes) in the first factorial plane differentiated by their A) Pit shape and B) Age. (Individuals factor map (PCA): A) Pit shape, and B) Age
To the left of the graph, with negative values are the subcircular pit lakes, formed between 1987 and 1997, and the elongated subrectangular pit lakes, mainly developed between 1998 and 2007. They have high coastline development values and low width, mean depth, perimeter, length, maximum depth, area and volume, respectively. The Fabero-20 and Fabero-27 pit lakes, which are circular, are morphometrically similar to the subcircular pit lakes. The Fabero-27 pit lake is smaller and shallow; its waters have an acidic pH compared to Fabero-20, which is more recent with an alkaline pH.