2.1. Regional Environmental Characteristics
The study area is located in central-southern Poland, in the vicinity of Kielce city (
Figure 1), with a population of almost 200,000 inhabitants (18th most populous city in the country). The water supply is entirely covered by groundwater intakes scattered around the periphery [
29]. The conducted research, described in this article, was related to the updating of the water resources of the largest municipal intake.
Physiographically, the study region belongs to the Polish Uplands [
30], remaining on the border of the Paleozoic Holy Cross Mountains and their younger, northwestern undulating margin. Low and rounded mountain ranges extend in a WNW–ESE direction. The highest peaks in the neighborhood reach absolute heights of 400 m above sea level, while the valleys and foothills are located about 100–150 m lower.
The rivers flowing between the mountain ranges frequently cut through the edge parts of the margin with their breakthrough sections. Drainage mainly takes place from the north towards the south and further in the southwest direction, towards the vast valley of the main polish river—Vistula. The most important rivers are, starting from the west: Łososina (Wierna River), Bobrza together with its tributaries (Silnica and Sufraganiec), Lubrzanka and Bełnianka flowing into Czarna Nida and further to Nida—a direct tributary of Vistula River.
Geologically, the southeastern part of the area is made up of strongly eroded Holy Cross Mountains, formed during the Caledonian and Hercynian orogeny. The northwestern part, fringing the mountain ranges, was covered by sediments during the Permian–Mesozoic transgressions, which as a result of later processes, were slightly deformed to form large-radius folds. Marine sediments are present in the open to the western valleys of the Holy Cross Mountains, covering parts of the synclines with increasing to the west layer. The study area includes two adjacent synclines, the Kielce syncline (to the north) and the Gałęzice–Bolechowice–Borków syncline (to the south). Centers of the concave parts of the folds are filled with Devonian carbonate sediments up to 1000 m thick (
Figure 2), to a much lesser extent covered with Carboniferous sediments (MGWB no. 418). Lithologically, the Middle Devonian dolomites dominate, with limestone and marl in the minority. Cores of the anticlines are made of formations of the Lower Palaeozoic age, mainly Cambrian. Lithologically, clay shales, siltstones and greywackes dominate. The sediments covering the northwest represent Permian and Triassic ages. Within cores of the Mesozoic age synclines, Upper Triassic sediments were documented. The central parts of the anticlines are formed by Lower Triassic formations. Lithologically, these are claystones, siltstones, calcareous sandstones and conglomerates, while in the upper part of the profile mainly limestones and marls. The prevailing part of the terrain surface is covered by a layer of Quaternary sediments, mainly of post-glacial origin. In the river valleys, fluvial and fluvioglacial deposits appear.
Most water-bearing structures are Devonian syncline cores made of thick dolomites. These dolomites, as well as other carbonate formations, are subject to intensive mining exploitation. They form irregular, elongated structures of fissure and karst character. The carbonate formations that make up the aquifer are characterized by very variable permeability, from very good to poor. The hydraulic conductivity values range from 1.3 × 10
−7 to 1.2 × 10
−3 m/s [
31,
32,
33]. The most water-rich fragments, forming so-called Major Groundwater Basins (MGWB), were separated within the boundaries of the formations. According to the national classification, a MGWB is a geological structure or its fragments with the highest water-bearing and storage capacity in the scale of hydrogeological regions. In the area under consideration, two such reservoirs have been distinguished: MGWB No. 417 Kielce in the north and MGWB No. 418 Gałęzice–Bolechowice–Borków in the south (
Figure 1).
The older Paleozoic formations are generally classified as isolating. Carboniferous sediments have a similar status. Permian and successively Lower Triassic formations appearing from the west are combined into one hydrogeological structure due to sedimentary continuity on the regional scale. Their water-bearing capacity is mainly connected to the systems of fractures and fissures. The water-bearing sediments are sandstones, fractured mudstones, conglomerates and locally marly limestones, separated by numerous layers of clays and loams. In such a lithological setting, the aquifer is classified as a fissure and fissure–pore, separated by numerous poorly permeable layers.
On the western side, water-bearing formations are covered by Middle Triassic sediments developed from the bottom of the profile as carbonate conglomerates, sandstones, limestones and marls. Next in the sedimentary sequence, the Upper Triassic sediments are formed in the form of practically impermeable clays and mudstones, with some minor exceptions. Already in the first half of the 1990s, hydrogeological documentation defining disposable water resources in the vicinity of Kielce [
34] recognized the northern fragments of the Permian–Triassic cover as a recharge zone for the northernmost Devonian reservoir MGWB No. 417 Kielce. The same sediments, however, together with far to the NW weakly permeable Lower Paleozoic formations of the Dyminy and Rykoszyn anticline, partially separating MGWB no. 417 from more southerly located MGWB no. 418, were regarded so far as providing sufficient isolation of both principal water-bearing structures. In all documentation and scientific publications to date, the two structures have been treated separately and analyzed independently. Separate, independent functioning digital hydrogeological models were developed for each of these structures [
26,
35,
36,
37], whose concepts and realized predictions originally did not raise any objections.
2.2. Mining Exploitation
Middle Devonian carbonate sediments, up to 1000 m thick, apart from a rich water reservoir, are also a valuable mineral. For several dozens of years, they have been intensively mined, becoming aggregate, fertilizer or raw material for cement and lime production. Both domestic and well-known global companies conduct mining activities here: CRH Materials Poland Ltd., Dyckerhoff Poland Ltd., Opencast Mines of Raw Materials for Roads S.A. in Kielce, Lafarge Aggregates Ltd., Świętokrzyskie Mineral Raw Materials Mines Ltd., and Nordkalk Ltd. However, the intensity of the activities varies greatly. Within the Kielce syncline and its recharge region, a large part of which is covered by urban development in Kielce and its outskirts, exploitation of rock materials is very local and is currently limited to two adjacent workings in the Laskowa Góra and Kostomłoty deposits, located to the NW of Kielce (
Figure 1). The situation is completely different in the Gałęzice–Bolechowice–Borków syncline, located to the south. The land use is dominated by farmland, and in the eastern part by forests. Scattered buildings dominate, concentrated in a few villages. Only along the national road No. 73 there are buildings typical of urban suburbs, built by inhabitants of neighboring Kielce seeking more secluded places. This makes it possible to identify, document and exploit further deposits. Mining is concentrated in the central part of the syncline. On an area of about 27 km
2, four mines of different owners are adjacent to each other, conducting exploitation on 6 pits until recently. Currently, after the water reclamation of the Radkowice-Podwole pit [
26], there are five operating, on the deposits: Jaźwica, Kowala, Kowala Mała, Trzuskawica (two pits Trzuskawica and Kowala;
Figure 1). In the northwestern part of the structure, the Ostrówka deposit is in operation, and the neighboring Ołowianka-1 deposit is in the process of obtaining the necessary documents.
On average, about 10 million tons of raw materials are mined annually in the syncline, with less than a million in the adjacent structure [
38]. Among the active six pits, the lowest mining ordinate of 150 m above sea level was reached at the Ostrówka deposit (Miedzianka mine). However, in the center of the structure, the ordinates of exploitation decrease to a maximum of 180 m a.s.l. Intensive and deep exploitation, depressing the water table to 103 m, requires pumping out on average annually about 28.5 million m
3 of water (78.14 thousand m
3/d; data for 2018). As a result of the superposition of local impacts, an extensive regional depression cone has formed, covering the vast majority of the structure, excluding the eastern part. At the peripheries of the Kielce syncline, the extraction of 0.88 million t of mineral is associated with the pumping of 1.92 million m
3 per year of water along with the lowering of the drainage ordinate to 219.5 m above sea level.
The exploitation of Ostrówka deposit, which is particularly important in the light of further considerations, reaches the deepest level, to the mentioned ordinate of 150 m above sea level. Lowering of water level by maximum of 103 m concerning its original location generates inflows of 34,400 m3/d (12.6 million m3 per year). Due to the depletion of resources, it is planned to start the exploitation of the neighboring Ołowianka-1 deposit, situated closer to the edge of the Mesozoic structure of the Holy Cross Mountains.
2.3. Characteristics of the Problem
The documentation works carried out for the preparation of Appendix no. 3 to the hydrogeological documentation establishing the exploitable resources of the municipal intake in Kielce—Białogon [
36] brought unexpected results of measurements of the groundwater table position at some observation points. Consequently, the observed changes in the hydrodynamic field of strictly local character, after an in-depth analysis, caused a significant modification in the perception of groundwater circulation in the regional system. The basis for the reorientation of the concept of regional groundwater circulation, resulting from the distribution of hydrodynamic field, were the results of measurements of the groundwater table position of the Middle Devonian aquifer in three boreholes carried out in 2017 (out of a total of 67 boreholes measured for documentation).
The mentioned wells were drilled in the extreme southwestern part of the Kielce MGWB (
Figure 3). In these three key boreholes, in the vicinity of Zagórze and Jaworznia villages, a decrease in the water table ordinates within the Devonian horizon was recorded by only 1.40 to 2.21 m, compared to the measurements from 2014. Although the changes were not large at first glance, they resulted in a change in the water table inclination towards the southwest, i.e., in the direction opposite to the center of the structure (where intake well drainage dominates) and in the formation of a local watershed.
The first information suggesting the initiation of such a tendency appeared in the report summarizing the state of groundwater in the vicinity of the Białogon intake in 2016, prepared by Kielce Waterworks. Only the update of the regional hydrodynamic model of water circulation, carried out by the authors of the article in 2017, in connection with the need to update the intake’s exploitation resources, allowed an in-depth analysis of the reasons for such an arrangement. As a consequence of the model calibration activities, the original assumptions of the concept of water circulation in the considered structure had to be changed, and the model itself was subjected to a relevant modification. A description of the implemented changes and the obtained effects is presented in
Section 3.1.
The interpretation of the change in the direction of groundwater table inclination, contrary to the lean recorded for years towards the center of the structure, where the Białogon multi-well intake is located, resulted in the necessity to implement in the conceptual model an additional drainage center with significant impact, located in the direction opposite to the intake wells. The detailed analysis did not reveal the presence of other significant drainage facilities in the immediate vicinity. Individual wells, managed by local water users, were operated at low capacities. The nearest wells, located in Piekoszów, 2.5–3.5 km away towards the NWN (
Figure 1), exploited approximately 100,000 m
3/a (Piekoszów PGO well) and 160,000 m
3/a (Piekoszów well) from the Lower Devonian aquifer, averaging less than 450 m
3/d each (data from 2015 y., [
37]). To the west, there are also some wells of the abandoned intake in Zawada, within sandstones of Lower Triassic and Permian age (about 2.5 km to the SW) and one more near the railroad station in Piekoszów (about 2.5 km to the NW,
Figure 1). The latter drew the attention of the authors earlier, during work aimed at forecasting the volume of inflow to the “Ołowianka-1” Devonian dolomite deposit planned for resumption, in the vicinity of a large active mining site exploiting the “Ostrówka” deposit [
37]. The hydrogeological model of the Gałęzice–Bolechowice–Borków syncline with MGWB No. 418 included the northwestern margin of the Devonian aquifer built of Permian and Triassic sediments (similarly to the surroundings of MGWB No. 417). The northern boundary of the model was drawn near Piekoszów, along the Czarne Stoki watercourse. Measurements of the water table position for the model calibration were conducted in June 2016. The measurements included a borehole near the railway station in Piekoszów, which had been inaccessible for a long time (secured by the owner), provided access to Triassic sandstones (T
1). After 39 years (data from the borehole drilling in 1977), the static water table decreased by 36.53 m to 214.97 m a.s.l. The same borehole was noticed by Prażak and his team [
39], who, based on the results of archival measurements from the period of the borehole drilling, for the first time hypothesized that water filtration in the area of the borehole is directed towards the east, to the Bobrza River. However, at that time, the thesis was not confirmed by measurements.
On the other side, looking more broadly, data from a few wells located in the Piekoszów vicinity are not conclusive. The wells operated by ZUK in Piekoszów (outside the model area), adjacent to the PKP well from the north, despite their distance of about 2 km, i.e., twice less than the analyzed well to the “Ołowianka-1” deposit, show significantly different location of the water table. In the case of the “Piekoszów” well, the ordinate of the static water table reached 252.79 m above sea level, while in the case of the “Piekoszów PGO” well—251.93 m above sea level (June 2016). In the following year, measurements conducted to document the resources of the Białogon intake indicated ordinates of, respectively: 252.3 and 252.6 m a.s.l.
Even more controversial were the boreholes located further to the south, in the direction of the “Ostrówka” deposit drained to the level of 160 m a.s.l. (currently 150 m a.s.l.) and the planned reopening of the neighboring “Ołowianka-1” deposit (
Figure 1). In Rykoszyn, measurements made in a private well with a depth of 34 m, filtered in Lower Triassic sediments, showed a static water table at the ordinate of 263.31 m a.s.l. (June 2016). Another borehole filtered in T
1 formations, near vicarage buildings in Rykoszyn, distant from the above-mentioned by 1 km towards E, with a depth of 39.5 m, has been dry for years (bottom ordinate 228 m a.s.l.). In turn, the mentioned intake in Zawada, consisting of two wells with a depth of 90 m each, drained Lower Triassic and Permian (Zechstein) sandstones. The approved resources of the intake are 9.7 m
3/h, with depression
se = 22.0 m. Since the expiration of the water permit in 2006, the operation has been abandoned due to low yield and unfavorable hydrogeological parameters.
In the year 2018, Piekoszów municipality decided to build a new well to improve the water supply in the southern part of the commune. For the location of the new facility, a plot of land in the village of Lesica was chosen, less than 3 km to the NW from the edge of the pit. The 60 m deep well exploits the Middle Triassic limestone aquifer (T2), which overlies the T1 sediments to the west. After being drilled to a depth of 6.0 m, the water table rose by 0.6 m and stabilized at a depth of 5.4 m below the ground surface. The exploitation resources of the intake were determined to be 40.0 m3/h, with a depression of 1.5 m. Operations began in the fourth quarter of 2018 with an average yield of approximately 260 m3/d. Since then, the well has been operating undisturbed with more than double the average yield (660 m3/d).
Such diverse and inconsistent information, obtained for the needs of various studies, has not always given sufficient attention to the peripherally located zone of Piekoszów. It was only when a well near the railway station in Piekoszów was unsealed for observation and different results of water table measurements were obtained in wells located in the SW part of the main groundwater basin Kielce (Jaworznia and Zagórze areas), that revealed the need to conduct a comprehensive analysis of numerous data from scattered boreholes and verification of the existing assumptions.
The key object for understanding the problem became the mine exploiting the Ostrówka deposit. The exploitation of the deposit historically linked to the Ołowianka deposit on the documentation stage, separated by the Hutka River valley, started in the late 1960s. The beginning of the exploitation is connected with the slope excavation, cutting the existing hill, starting from the ordinates of 290 m above sea level, within the highest sublevel of the floor I. The next floor II, at 236 m a.s.l., was dug below ground level in the classic form of an underground pit. Excavation was carried out using the longwall method with explosives. The encounter of the original water table took place slightly higher, at the ordinate of about 250 m above sea level, which was associated with the beginning of the pit drainage. With the opening of successive floors (221, 206, 192, 176 m asl), the depression increased, and the depression cone caused by it expanded. Initially, the depression did not extend beyond the Middle Devonian aquifer. With time, however, as it reached deeper resources, it entered laterally adjacent, poorly permeable Triassic strata, fringing together with younger Mesozoic sediments the Holy Cross Mountains (
Figure 4A). Works focused at the ordinate of 160 m above sea level commenced in the middle of 2009. The increase with the deepening of the works inflow reached 31 thousand m
3/d (2015). A few years later, in 2015–2016, the lowering of the water table within the Lower Triassic formations, separating the Devonian syncline from the west, reached the region of Zagórze and Jaworznia, entering the SW part of the Devonian-filled Kielce syncline simultaneously (
Figure 4B). The achieved state should be consolidated with the opening of the 8th level within the Ostrówka deposit, at the ordinate of 150 m a.s.l. in the third quarter of 2017.
The Miedzianka mine has the right to exploit two neighboring deposits: “Ostrówka” and “Ołowianka 1”. They are separated by the valley of the Hutka River. The object of exploitation is Middle Devonian limestones intended for the production of aggregates. Production within the active deposit “Ostrówka” requires pumping an average of 31 thousand m
3/d [
37]. In June 2009, the seventh level was opened at the ordinate of 160 m a.s.l. Eight years later, in August 2017, the exploitation was lowered by another 10 m to the ordinate of 150 m a.s.l. At the deposit “Ołowianka 1”, work was discontinued many years ago. Since 2016, preparations have begun to resume operations.
The Triassic aquifer is de facto a water-bearing complex, i.e., a system comprising more or less isolated water-bearing horizons developed in lithologically diversified Triassic layers. The water-bearing horizons include the clastic formations of the Bunter sandstone, carbonates of its highest part—the Röt—as well as carbonates of the middle Triassic—the Muschelkalk. The younger and younger formations appear on the surface towards the west, eventually disappearing under the cover of impermeable formations of Carnian. Amongst the clastic formations of the Bunter sandstone, the reservoir rocks are sandstones with interbedding clays and siltstones with porosity of about 8 to 20% (average 15%), as well as limestones of the highest part of the Lower Triassic—the Röt. Sandstones, characterized by large variability in porosity and fracturing, have correspondingly different hydrogeological parameters. Particularly high variability is observed in hydraulic conductivity, generally ranging from 9 × 10
−8 to 3.6 × 10
−4 m/s (0.008 ÷ 31.1 m/d, [
40]). With variable thickness, transmissivity values range from 1 to more than 1400 m
2/d. The Lower Triassic aquifer is unconfined in outcrop areas, changing to subartesian in the remaining areas. The isolating layers causing an increase in pressure are clay inserts overlying sandstones.
Muschelkalk limestones and dolomites, with thicknesses in the range of 20 ÷ 150 m, are fractured and karstified. They attain higher hydraulic conductivity than those of the Bunter sandstone, ranging from 1.7 × 10
−6 to 1.6 × 10
−4 m/s (0.15 ÷ 13.8 m/d, [
40]). Similarly, as in the case of the Bunter sandstone aquifer, the Muschelkalk limestone aquifer generally remains unconfined in the outcrop areas, changing its character to confined in the zones where it occurs under the cover of poorly permeable Upper Triassic (Carnian) formations, stabilizing at depths from a few to a maximum of several meters.
A significant degree of complexity of the aquifer system, with a simultaneously limited amount of hydrogeological information concerning water-bearing capacity and, in particular, hydraulic contacts of individual layers with each other and with neighboring units, significantly hinders the unambiguous representation of the details of the geological structure on the hydrodynamic model.
2.5. Numerical Model of Hydrogeological Conditions
The numerical hydrogeological model, which allowed determining the resources of the intake, was prepared based on the multilayer hydrogeological model of the Kielce Groundwater Exploitation Area—Subdivision A [
41], which has been in operation for several years, in the updated version for 2015 as part of the realization of the “Supplement to the Hydrogeological Documentation…” [
42]. The model was updated to a state corresponding to the averaged intake exploitation from 2016 and the results of groundwater table position measurements conducted in mid-2017. The discussed model includes the so-called Kielce Groundwater Exploitation Region—balancing subregion A. The groundwaters in this region are mainly exploited from the Middle Devonian aquifer MGWB 417 Kielce, southern wing of the Kielce syncline, and also, in smaller quantities, by individual wells, from Middle Permian and Lower Triassic, as well as Upper Devonian and Middle Triassic aquifers situated in the water recharge zone from the north [
41].
The numerical multi-layer model of groundwater flow (6 layers) was prepared using the Processing Modflow program [
43], which is based on the original MODFLOW code [
27]. The model study covered an area of about 130 km
2. The archival version was reconstructed based on the results of direct research and field measurements. The boundaries of the so-called subregion A refer to natural boundaries of hydrodynamic nature, associated with groundwater watersheds or rivers (
Figure 5). The model retains the concept of schematization of hydrogeological conditions adopted in the original hydrogeological model [
41]. Horizontally, the model study area was discretized based on a regular grid (125 × 125 m square cells). An active part of the modelled area is contained in a rectangle comprising of 114 rows and 155 columns, which corresponds to the real size of about 14.3 × 19.4 km. Layer 1 (
Figure 6) represents water-bearing near-surface Quaternary formations and water-bearing or poorly permeable older formations outcropping or occurring under a cover of Quaternary formations of small thickness (aeration zone). Layer 2 corresponds to poorly permeable Quaternary and Neogen formations (clays, loams, silts) and poorly permeable older formations outcropping or occurring under a cover of Quaternary formations of low thickness. The purpose of layers 3, 4, 5 and 6 is to map water-bearing formations of older bedrock (D
2, D
3, P
2, T
1, T
2) and, in the peripheral parts of the model, poorly permeable formations of older Palaeozoic (Cm, O, S, D
1) as well as Upper Triassic. Together, these layers constitute the principal useful aquifer.
In all layers, active cells were defined, corresponding to filtration area (
Figure 7—white cells), and inactive cells, not taking part in computation process (
Figure 7—grey cells). For rivers with good hydraulic contact between surface water and groundwater (Bobrza, Sufraganiec, part of Lubrzanka), boundary conditions of first-type H = const were assumed within layer 1 (
Figure 7—blue cells). The same condition was assumed for layer 4 in cells corresponding to the location of the Laskowa mine workings, which allows for simulation of its dewatering. Due to the limited hydraulic contact between surface- and groundwater, the remaining smaller watercourses were simulated with the boundary conditions of the third type using RIVER package (
Figure 7—light blue cells). Parameters assigned to individual cells in the RIVER module (hydraulic conductance of the riverbed) take into account both the hydraulic conductivity of the riverbed bottom sediments, their thickness and the size of the river in individual calculation cells. In addition, a graded head in the river and elevation of the riverbed bottom parameters were established.
The values of hydraulic conductivity (horizontal and vertical) corresponding to the modelled aquifers were entered into the cell zones and then modified during model calibration. The karst phenomena typical of carbonate aquifers have been documented, although their recognition is irregular and only local [
44]. Taking into account the regional scale of the research and the degree of discretization (cell size), the model calculations were based on the values of parameters typical of fissure–karst–pore aquifers. The distribution of hydraulic conductivity of individual layers obtained after model calibration is presented in
Table 1. It should be considered that extreme values usually occur only locally, in places with specific geological conditions.
Effective recharge of precipitation depends on the magnitude of precipitation and effective infiltration coefficient. In practical groundwater modeling, typically the so-called effective recharge (the amount of water that enters the rock mass) is adopted, ignoring the evaporation. Before model calibration, an effective recharge was estimated, based on average annual precipitation recorded for Kielce (629 mm, [
45]), the permeability of subsurface sediments and types of land development. Then, recharge from precipitation was assumed as a second-type boundary condition (Q = const), using the RECHARGE package with the option Recharge is applied to the highest active cells. After model calibration, the average amount of effective recharge from precipitation was 3.36 × 10
−4 m
3/d/m
2 (
Table 2), which is about 20% of the average annual rainfall in the Kielce region. The size of the infiltration rate varies depending on the permeability of the rocks present in the near-surface layers, the degree of land development and the location within the range of the depression cones of the Laskowa mine and groundwater intakes.
Intake wells with a constant yield were simulated using the second-type boundary conditions (Q = const). The model includes numerous deep wells, from which regular water abstraction was carried out in 2016. Constant yields from wells, defined as daily averages based on the actually registered abstraction in 2016, were assumed to total 25,029 m
3/d, within layer 3 (
Figure 8A) and layer 4 (
Figure 8B). Most of the wells work in cooperative conditions, remaining under the hydrodynamic influence of adjacent intakes.
Prognostic calculations were performed under steady-state flow simulation. This research method comes down to the determination of the target state resulting from the assumed stresses, without indicating any intermediate states that characterize the temporal variability of the considered phenomenon.
2.6. Model Calibration and External Influence of Mine Activity
For numerical calculations of groundwater filtration equations, a simulator based on the finite differences method (FDM) from the widely used MODFLOW family [
11] was used, with MODFLOW-2005 [
27] being the most commonly used version at present. The prepared model was subjected to a calibration process, the aim of which was to obtain accordance between hydraulic heights obtained from model studies and those measured during field mapping, as well as between inflows to the Laskowa pit simulated in the model and those measured in reality. For the calibration results of own field measurements of the groundwater table depth in 67 boreholes were used, carried out in July 2016. After analysis, a selection of measurement points was made; finally, for calibration, water table location information was used from 42 boreholes not operating or water abstraction was marginal.
The trial-and-error method [
11] was applied to calibrate the model. During the calibration, the spatial distribution of hydraulic conductivity was modified. Slight changes were also made in values of vertical hydraulic conductivity which co-determines the amount of seepage between individual layers. As part of the calibration, also the recharge from precipitation infiltration was corrected. In addition, the hydraulic conductance of the riverbed, which determines the contact of surface water with groundwater, was slightly modified.
The results of the own field measurements of the groundwater table position from July 2016 and previous hydrogeological measurements (2015) in municipal wells and piezometers under local monitoring conducted by Kielce Waterworks clearly indicate that the hydrodynamic field in a part of the considered area had changed compared to previous observations and interpretations [
41]. Monitored points located in the southwestern region of the modeled structure (30, P17) indicate lower water table ordinates than points closer to the Kielce—Białogon intake, e.g., 37, P8, P8A, V (
Figure 3). Since in the southwestern region of the studied aquifer the exploitation of groundwater is not intensive (
Section 2.3) and the existing watercourses are not capable of producing such a lowering of the water table, the only explanation is the drainage of strata occurring outside the area of direct model tests, where size and extent affect the hydrodynamic field system in the area under consideration. The analysis of the archival materials and documentary studies leads to the conclusion that the existing layout of the water table is most probably influenced by a large drainage center in the “Ostrówka” deposit exploited by the Miedzianka mine located some distance to the southwest (
Figure 1). Due to the location of the mine in a separate hydrogeological aquifer (Gałęzice–Bolechowice–Borków syncline), previous studies assumed its very limited influence on other neighboring structures. Current measurements, however, indicate that such influence, although relatively small, is noticeable, reaches other water-bearing structures and cannot be ignored in the considered issues. The original lack of water exchange between the two adjacent hydrogeological structures likely could have been disturbed by intensive mining drainage of the Miedzianka mine workings (nearly 1300 m
3/h), significantly larger than the groundwater abstraction by the Białogon intake (about 900 m
3/h; [
36,
37]). The influence of such drainage could be revealed after a sufficiently long time and further lowering of the drainage ordinate in the Ostrówka pit, in the form of a locally different than originally shaped water table in the Kielce Groundwater Exploitation Region, causing relatively small groundwater flows between structures.
It was necessary to include in the model the influence of external stressing factors coming from outside the modeled area. For this purpose, the General-Head Boundary module [
43] was used, which allows for mapping with third-type boundary condition of external stresses. In the extreme blocks of the model from the southwestern side (
Figure 9), in layers mapping older formations (3, 4, 5 and 6), conductivity to a distant boundary, calculated taking into account low values of filtration parameters and distance of the mine from the model, as well as ordinate of the water table in the drained pit (160 m above sea level), was assumed. During further calibration, the adopted values of conductivity were locally slightly modified. A detailed description of the results of measurements taken from wells situated in the area of the probable impact of drainage of the neighboring structure (SW part of the Major Groundwater Basin in Kielce) is given in
Section 2.3.
As a result of the calibration, a numerical hydrogeological model was developed, that generates the results (the water table arrangement and the flow to dewatered Laskowa excavation) similar to those observed in reality (the so-called Variant 0—current hydrodynamic state after model calibration). For the adopted calibration points, the differences between the measured water table and obtained by interpolating the resulting matrix distribution from the model are insignificant and, in principle, do not exceed the value of a single meter. The arrangement of the measurement points directly adjacent to the diagonal of the calibration plot (
Figure 10) means a good fit of the model’s response to the actual observations.
Considered as calibration criteria, the error rates are: mean error ME = 0.13 m, mean absolute error MAE = 0.57 m. Normalized value is related to the amplitude of water table fluctuation in the considered structure (approx. 80 m); mean absolute error value does not exceed 1%. The volumes of inflows to the Laskowa open pit obtained from the model after calibration (5271 m3/d) concerning inflows observed in the year 2016 (5261 m3/d) differ to a very small extent.