Low Enthalpy Geothermal Resources for Local Sustainable Development: A Case Study in Poland

Many regions in Poland face the problem of air pollution. These regions include, though are not exclusive, to health resorts, an important element of the healthcare industry. Poor air quality is mainly associated with the domestic sector, which is dominated by individual solid fuel and coal boilers. One option for reducing emissions is to use alternative energy sources for heating purposes. Therefore, the paper presents an analysis into the possibility of using low enthalpy (low temperature) geothermal energy in health resort areas. The main purpose of the article is to estimate the potential of soil and water as the lower source for a heat pump. The article presents analyses of geological and hydrogeological conditions based on which the thermal parameters of the rock mass were determined, which were thermal conductivity and unit heat output for 1800 operating hours per year. The calculated values were used to perform a spatial analysis of the data and create maps of the average thermal conductivity for horizontal ground heat exchangers (HGHE) to a depth of 2 and 10 m and vertical ground heat exchangers (VGHE) in depth ranges up to 30, 60, and 90 m. The heating power of the intakes, located in the research area, were estimated using empirical formulas. In addition, a detailed analysis of the physicochemical parameters of groundwater in terms of the requirements indicated by four heat pump manufacturers is presented. The results of the presented research makes it possible to assess the potential of low-temperature geothermal energy and to characterize the suitability of the selected location for the use of HGHE, VGHE, and wells. As a result, the discussed area was found to have a high potential for the use of ground source heat pumps and a moderate potential for the use of low-temperature systems based on groundwater.


Introduction
The use of energy from renewable sources is one of the priorities of the European Union (EU) energy policy. Renewable energy sources (RES), in the face of continuous climate change, are fundamental to reducing greenhouse gas emissions [1,2]. Locally, they also have a significant impact on the improvement of air quality [3,4]. In addition, they form an integral part of local sustainable development, which plays an important role in the daily life of society.
Poland is an EU member state that has been struggling with air pollution for many years. According to data from the European Environmental Agency (EEA), the highest concentrations and breaches of the permissible levels were observed for particulate matter PM10 and PM2.5 and benzo(a)pyrene (B(a)P) in 2017 [5]. Exceedance of the daily limit value of PM10 (50 µg/m 3 ) was recorded at all stations located in the central and southern parts of the country with an annual average of 60.2 µg/m 3 in 2017 across the whole country. Annual limit values of PM10 and PM2.5 were not exceeded nationwide, however they were recorded locally (stations in the south of the country).
biggest problem concerns the average annual concentration of B(a)P (EU target value 1 ng/m 3 ), which has been exceeded for many years throughout the country with an annual average for the country of 4.8 ng/m 3 in 2017 [5,6].
Poor air quality also affects health resort areas, which in Poland form an integral part of the healthcare industry. Treatment in these areas is based on the use of natural minerals and clean climate. Therefore, the Rabka-Zdrój health resort, located in the southern part of Poland, was chosen as a testing ground. The resort, for several years, has had a significant problem with high concentrations of PM10 and B(a)P, often accompanied by the phenomenon of ground-level emissions. The main reason for these emissions is the combustion of solid fuels and the production of domestic hot water in households. Furthermore, unfavorable meteorological conditions (such as temperature inversion, low wind speed, the high number of days with atmospheric calm), and the topography (location in a valley surrounded on three sides by hills) intensify air pollution. The goal of the present research is to develop solutions based on renewable energy sources, whose increased use would reduce pollutant emissions from domestic households and improve the quality of air at the health resort. The study analyzed the low enthalpy geothermal potential and the possibility of using solutions based on groundwater and obtaining heat from the earth.

Location and Climatic Conditions
Rabka-Zdrój is a mountain health resort situated in the southern part of Poland (49°37′ N, 19°58′ E) (Figure 1). The municipality has possessed the status of a health resort since 1967. According to Polish law, an area that has been granted such a status has to meet a total of five conditions, which include deposits of natural healing minerals and a climate with proven healing properties [7]. Healing waters are the source of the natural healing minerals found in Rabka-Zdrój including sodium chloride-rich mineral waters containing iodides and bromides [8,9]. These are used in the treatment of orthopedictraumatic diseases, cardiovascular diseases and hypertension, upper and lower respiratory tract diseases, diabetes and skin diseases [10]. A selection of meteorological parameters is presented in Table 1. The municipality has possessed the status of a health resort since 1967. According to Polish law, an area that has been granted such a status has to meet a total of five conditions, which include deposits of natural healing minerals and a climate with proven healing properties [7]. Healing waters are the source of the natural healing minerals found in Rabka-Zdrój including sodium chloride-rich mineral waters containing iodides and bromides [8,9]. These are used in the treatment of orthopedic-traumatic diseases, cardiovascular diseases and hypertension, upper and lower respiratory tract diseases, diabetes and skin diseases [10]. A selection of meteorological parameters is presented in Table 1. Table 1. Selected values of meteorological parameters and bioclimatic indicators (based on [11]).

Title Units
Mean annual air temperature 6. The area of the health resort is not subject to constant air monitoring. Air quality tests take place ad hoc in different years and cover the concentration measurements for B(a)P, PM10, NO 2 , and SO 2 . The results obtained from the monitoring station located in Rabka-Zdrój in 2014 and 2017 indicate very high concentrations of B(a)P in the autumn-winter period, which in Poland means the heating season. The average annual concentration of B(a)P in both cases is about 8 ng/m 3 , which is eight times the permissible level of 1 ng/m 3 specified by the European Union [12] (Figure 2). The permissible annual concentration of PM10 (40 µg/m 3 ) was not exceeded. However, in 2014 and 2017, values close to the limit were reached, amounting to 33.5 and 35 µg/m 3 , respectively. In addition, in 2017, 66 days were recorded where the permissible daily concentration of PM10 (50 µg/m 3 ) was exceeded, with a permissible value of 35 days [13]. In 2014, the target value of the daily concentration of PM10 was not exceeded. The majority of days, where thePM10 levels were exceeded, was recorded from October to March. The above data indicates a relationship between air quality and the heating of buildings, most often with the use of inefficient heating devices and poor quality solid fuels [14,15].  The area of the health resort is not subject to constant air monitoring. Air quality tests take place ad hoc in different years and cover the concentration measurements for B(a)P, PM10, NO2, and SO2. The results obtained from the monitoring station located in Rabka-Zdrój in 2014 and 2017 indicate very high concentrations of B(a)P in the autumn-winter period, which in Poland means the heating season. The average annual concentration of B(a)P in both cases is about 8 ng/m 3 , which is eight times the permissible level of 1 ng/m 3 specified by the European Union [12] (Figure 2). The permissible annual concentration of PM10 (40 µ g/m 3 ) was not exceeded. However, in 2014 and 2017, values close to the limit were reached, amounting to 33.5 and 35 µ g/m 3 , respectively. In addition, in 2017, 66 days were recorded where the permissible daily concentration of PM10 (50 µ g/m 3 ) was exceeded, with a permissible value of 35 days [13]. In 2014, the target value of the daily concentration of PM10 was not exceeded. The majority of days, where thePM10 levels were exceeded, was recorded from October to March. The above data indicates a relationship between air quality and the heating of buildings, most often with the use of inefficient heating devices and poor quality solid fuels [14,15].

Geological Structure
The research area is located in the Outer Carpathians, also known as the Flysch Carpathians, and is characterized by a complex geological structure. Between 1925-1927, nine boreholes were drilled including seven hydrogeological boreholes However, their documentation is incomplete. Detailed exploration of the geological record was made possible only by deep structural, exploratory, and hydrogeological drilling in the 1970s and 1980s [17]. The research area consists mainly of folded

Geological Structure
The research area is located in the Outer Carpathians, also known as the Flysch Carpathians, and is characterized by a complex geological structure. Between 1925-1927, nine boreholes were drilled including seven hydrogeological boreholes However, their documentation is incomplete. Detailed exploration of the geological record was made possible only by deep structural, exploratory, and hydrogeological drilling in the 1970s and 1980s [17]. The research area consists mainly of folded Energies 2020, 13, 5010 4 of 20 flysch and Quaternary deposits. The highest structural unit in the Rabka-Zdrój region is the Magura nappe covering Cretaceous-Paleogene formations.
Three of them were found in the area of the present study, from the north these are: the Racza subunit, the Bystrica subunit, and the Krynica subunit [24]. In the Racza subunit, above the Senon deposits, shales were observed to have developed, being covered with Hieroglyphic beds with a thickness of about 150 m. Hieroglyphic beds (middle Eocene) mainly form gray-green shales with thin-bedded calcareous-silicon sandstones. The predominant deposits of the Racza subunit are Upper Eocene thick-bedded sandstones with intercalations of conglomerate and shale ( Figure 3). In the Bystrica subunit, the profile of Paleogene formations begins with shale, with thick-bedded sandstone intercalations in the central part of the study area. Above, in the Lower and Middle Eocene, the Beloveza beds (shales, thin-bedded sandstones) and Łącko beds (marls, sandstones, shales) developed, respectively. In these beds, the presence of mineral waters was found by seven boreholes [25,26]. Above the Beloveza and Łącko beds, as in the Racza subunit, there are Hieroglyphic beds covered by the Magura beds. Along the northern part of the study area, the Krynica subunit formations are exposed. The dominant complexes are thick-bedded sandstones from Piwniczna and Magura sandstones [27]. Between them, the Kowaniec beds (Middle Eocene) developed. These consist of not only sandstones, but also thicker intercalations of marls, shales, and conglomerates. All the deposits representing Paleogen have local variegated slate intercalations [28]. flysch and Quaternary deposits. The highest structural unit in the Rabka-Zdrój region is the Magura nappe covering Cretaceous-Paleogene formations. The oldest deposits of the Magura nappe are alb-cenoman spotty green shales and cenomanturon variegated shales exposed to the east of Rabka. Higher in the profile there are Inoceramian beds of the Senon-Paleocene age, developed as thin, medium, and thick-bedded sandstone-shale series [18]. Due to the facial diversity, the Paleogene formations were divided into four tectonic-facies subunits: Siary (Su), Racza (Ra), Bystrica (Bu), and Krynica (Ku) [19][20][21][22][23].
Three of them were found in the area of the present study, from the north these are: the Racza subunit, the Bystrica subunit, and the Krynica subunit [24]. In the Racza subunit, above the Senon deposits, shales were observed to have developed, being covered with Hieroglyphic beds with a thickness of about 150 m. Hieroglyphic beds (middle Eocene) mainly form gray-green shales with thin-bedded calcareous-silicon sandstones. The predominant deposits of the Racza subunit are Upper Eocene thick-bedded sandstones with intercalations of conglomerate and shale ( Figure 3). In the Bystrica subunit, the profile of Paleogene formations begins with shale, with thick-bedded sandstone intercalations in the central part of the study area. Above, in the Lower and Middle Eocene, the Beloveza beds (shales, thin-bedded sandstones) and Łącko beds (marls, sandstones, shales) developed, respectively. In these beds, the presence of mineral waters was found by seven boreholes [25,26]. Above the Beloveza and Łącko beds, as in the Racza subunit, there are Hieroglyphic beds covered by the Magura beds. Along the northern part of the study area, the Krynica subunit formations are exposed. The dominant complexes are thick-bedded sandstones from Piwniczna and Magura sandstones [27]. Between them, the Kowaniec beds (Middle Eocene) developed. These consist of not only sandstones, but also thicker intercalations of marls, shales, and conglomerates. All the deposits representing Paleogen have local variegated slate intercalations [28]. The youngest deposits of the testing ground are those from the Quaternary age, forming a discontinuous cover of varied thickness over the older (Cretaceous-Tertiary) deposits. Pleistocene formations in the study area developed as gravels, sands, and river clays of erosion and accumulation terraces [29]. The sediments of the oldest South Polish glaciation created terraces of 15-25 m above the river level, formed between the Skomielny and Luboński streams as well as in the area between Raba, Poniczanka, and Słonka. The deposits created during the Middle Polish glaciation formed terraces of 5-9 m above the river level (the area around the mouth of the Słonka and Poniczanka streams and flowing down from Luboń to Raba). Respectively, the height of terraces created during The youngest deposits of the testing ground are those from the Quaternary age, forming a discontinuous cover of varied thickness over the older (Cretaceous-Tertiary) deposits. Pleistocene formations in the study area developed as gravels, sands, and river clays of erosion and accumulation terraces [29]. The sediments of the oldest South Polish glaciation created terraces of 15-25 m above the river level, formed between the Skomielny and Luboński streams as well as in the area between Raba, Poniczanka, and Słonka. The deposits created during the Middle Polish glaciation formed terraces of 5-9 m above the river level (the area around the mouth of the Słonka and Poniczanka streams and flowing down from Luboń to Raba). Respectively, Energies 2020, 13, 5010 5 of 20 the height of terraces created during the North Polish glaciation reached 2.5-7 m above the river level (the region of the upper course of the Raba river with areas of larger tributaries) [18]. Loess-like, diluvial and eluvial clays with rock debris (locally located on the Raba River around the center of Rabka-Zdrój) are locally found in the analyzed area as well as blocks, rock debris, clays, and colluvial sands (associated with mass movements). Holocene gravels, sands, and alluvial mud aggregates appear in the Raba Valley and its tributaries [29].

Hydrogeological Conditions
According to the latest regional division of ordinary groundwater in Poland [30], the research area belongs to the southern macro-region and includes the Carpathian Region of the mountain province (sub-region of the Outer Carpathians). This division refers to earlier studies presented by Paczyński and Kleczkowski [31,32]. Fragments of two Tertiary main and local groundwater bodies occur in the study area. In the northern part, there is the MGWB No. 439 Magura bed water reservoir (the Gorce Mountains), while in the southern part, there is the LGWB No. 445 Magura bed water reservoir (Babia Góra mountain)-formerly MGWB No. 445 ( Figure 4) [33]. In the area, there is a fracture-pore aquifer in flysch (Cretaceous-Palaeogene) deposits and a pore aquifer in Quaternary deposits [34]. In the Quaternary deposits, the conductivity of the aquifer is 0.83-12.5 m 2 /h, and the filtration coefficient is 3.6·10 −1 -3.6·10 −4 m/h. In the Cretaceous-Paleogene deposits, aquifers occur up to a depth of 78 m and are characterized by a filtration coefficient of 3.6·10 −2 -3.6·10 −3 m/h [35].
The Quaternary aquifer is formed by sandstone and gravel sediments of a river origin with a thickness not exceeding 10 m. This level is fed mainly by the infiltration of atmospheric precipitation and surface waters. In the Raba River Valley, the thickness of the aquifer is up to 4.8 m, and in the valleys of its tributaries, it does not exceed 2 m ( Figure 4) [36]. The potential capacity of a well is defined as within 2-5 m 3 /h. The water table stabilizes at a depth of 2-5 m. In the central part of Rabka-Zdrój (the area of alluvial fans), the hydrogeological conditions slightly change. The potential capacity of a well is lower (0.5-1.5 m 3 /h), and the confined groundwater surface is found at a depth of 6-11 m [37]. The aquifers of the flysch formation of the Magura unit are characterized by thickand medium-bedded sandstones with shale intercalations [25,38]. In the flysch beds of Rabka-Zdrój, the co-occurrence of mineral and ordinary waters have been observed [39]. In the northern part of the research area in the Racza unit, aquifers occur in the Hieroglyphic beds, the Magura beds, and also locally in sub-Magura beds. In the central part of the research area, aquifers are associated with the Inoceramian beds as well as the Łącko, Beloveza (Sądecki unit), and Krosno beds (the Grybów and Dukla units). Further south, aquifers occur in the Inoceramian beds as well as in the Piwniczna sandstones of the Krynica unit [37]. The direct supply of the Tertiary beds occurs via the infiltration of precipitation. The main groundwater flows run toward the Raba River Valley and locally toward smaller streams [36]. In the analyzed area, there were also zones characterized by the lack of a usable aquifer [40]. This situation usually occurs when there are more shale layers than sandstone formations.

Materials and Methods
In order to assess the potential of low-temperature geothermal energy, a number of available databases containing geological and hydrogeological information have been used. To analyze the possibility of obtaining ground-source heat, data published and made available by the Polish Geological Institute (PGI) was used. This included the Central Geological Database (CBDG) from which information on boreholes located in the research area was obtained as well as the Central Hydrogeological Database (CBDH) from which information on hydrogeological wells and groundwater intakes was obtained. Detailed data concerning aquifers, lithological profiles of boreholes, and the efficiency of aquifers were included. In the process of the data analysis, the Detailed Geological Map of Poland, the Hydrogeological Map of Poland, and the Geoenvironmental Map of Poland were also used as well as available explanations, which provided a source of additional and more detailed information. The analyses was supplemented with publicly available PGI digital and spatial data resources and other open public data [41,42].

Methods of Using the Earth's Heat as a Lower Heat Source for Heat Pumps
The first stage of the study was to analyze the possibility of using ground-source heat for households with vertical and horizontal ground heat exchangers. For this purpose, lithological profiles were used for boreholes located in Rabka-Zdrój. In addition, in order to expand the analysis of the selected area, it was decided to appoint a buffer zone with a radius of 2 km and use the available borehole data from the immediate vicinity of administrative boundaries ( Figure 5). Based on the lithological profiles, the basic geothermal properties of rocks and soils were identified. (1)

Materials and Methods
In order to assess the potential of low-temperature geothermal energy, a number of available databases containing geological and hydrogeological information have been used. To analyze the possibility of obtaining ground-source heat, data published and made available by the Polish Geological Institute (PGI) was used. This included the Central Geological Database (CBDG) from which information on boreholes located in the research area was obtained as well as the Central Hydrogeological Database (CBDH) from which information on hydrogeological wells and groundwater intakes was obtained. Detailed data concerning aquifers, lithological profiles of boreholes, and the efficiency of aquifers were included. In the process of the data analysis, the Detailed Geological Map of Poland, the Hydrogeological Map of Poland, and the Geoenvironmental Map of Poland were also used as well as available explanations, which provided a source of additional and more detailed information. The analyses was supplemented with publicly available PGI digital and spatial data resources and other open public data [41,42].

Methods of Using the Earth's Heat as a Lower Heat Source for Heat Pumps
The first stage of the study was to analyze the possibility of using ground-source heat for households with vertical and horizontal ground heat exchangers. For this purpose, lithological profiles were used for boreholes located in Rabka-Zdrój. In addition, in order to expand the analysis of the selected area, it was decided to appoint a buffer zone with a radius of 2 km and use the available borehole data from the immediate vicinity of administrative boundaries ( Figure 5). Based on the lithological profiles, the basic geothermal properties of rocks and soils were identified.   Heat conduction coefficients and the unit heat efficiency for individual lithological precipitates were adopted based on the literature data, also taking into account the moisture content of the deposits. For sandstones, mudstones, and claystones, the coefficients were determined based on measurements of rock samples presented by Sroka and Plewa [43,44]. Recommended coefficients were used for other deposits obtained by combining the results of measurements performed as part of the Cheap and Efficient Application of Reliable Ground Source Heat Exchangers and Pumps (Cheap-GSHPs) project with results reported in the literature [45]. The missing coefficients, for example, for shale, were adopted from the Verein Deutscher Ingenieure (VDI) 4640-1 guidelines [46] and Polish Organization for the Development of Heat Pump Technology (PORT PC) [47]. The unit heat output coefficients were adopted in accordance with Sanner [48] for vertical ground heat exchangers and the VDI 4640-2 guidelines [49] for horizontal ground heat exchangers. Based on the results obtained for the average value of the thermal conductivity λ [W/m•K], a spatial analysis was performed using Geographic Information System (GIS) tools. Low-temperature energy potential maps up to depth intervals of 2 and 10 m were developed for horizontal heat exchangers and 30, 60, and 90 m for vertical heat exchangers.

Methods of Using Groundwater as a Lower Heat Source for Heat Pumps
The second stage of the study was to analyze the possibility of using groundwater for households as a lower heat source for low-temperature systems supported by heat pumps. As in the Heat conduction coefficients and the unit heat efficiency for individual lithological precipitates were adopted based on the literature data, also taking into account the moisture content of the deposits. For sandstones, mudstones, and claystones, the coefficients were determined based on measurements of rock samples presented by Sroka and Plewa [43,44]. Recommended coefficients were used for other deposits obtained by combining the results of measurements performed as part of the Cheap and Efficient Application of Reliable Ground Source Heat Exchangers and Pumps (Cheap-GSHPs) project with results reported in the literature [45]. The missing coefficients, for example, for shale, were adopted from the Verein Deutscher Ingenieure (VDI) 4640-1 guidelines [46] and Polish Organization for the Development of Heat Pump Technology (PORT PC) [47]. The unit heat output coefficients were adopted in accordance with Sanner [48] for vertical ground heat exchangers and the VDI 4640-2 guidelines [49] for horizontal ground heat exchangers. Based on the results obtained for the average value of the thermal conductivity λ [W/m·K], a spatial analysis was performed using Geographic Information System (GIS) tools. Low-temperature energy potential maps up to depth intervals of 2 and 10 m were developed for horizontal heat exchangers and 30, 60, and 90 m for vertical heat exchangers.

Methods of Using Groundwater as a Lower Heat Source for Heat Pumps
The second stage of the study was to analyze the possibility of using groundwater for households as a lower heat source for low-temperature systems supported by heat pumps. As in the previous stage (3.1.), it was decided to extend the research area to a buffer zone with a radius of 2 km ( Figure 5). The basic parameter for determining the energy potential of groundwater was the heating power of its intakes. Its estimated value was calculated for all of the analyzed intakes based on the following formula [50]: where Q geot is the heating power of the intake [kW]; Q w is the efficiency of the intake [m 3 /s]; P w is the water density, which was set at 1000 (kg/m 3 ); c w is the specific heat of the water (set at 4,19 kJ/(kg· • C)); T w is the water temperature at the head ( • C); and T z is the temperature of the cooled water (temperature at the outlet of the heat pump) ( • C). The efficiency of individual intakes was adopted based on data obtained from the Polish Hydrogeological Survey (PHS) [51] The water temperature at the head was adopted as the temperature of water taken from individual intakes. During the data verification stage, it was found that nine of the analyzed intakes had data concerning water temperature. For the other intakes, it was decided to calculate the estimated water temperature using the formula [52,53]: where T is the water temperature at the depth of H ( • C); t avg is the average annual air temperature ( • C) [54]; A is a constant depending on the height above sea level [52,53]; H is the depth of the occurrence of water (m); h is the depth of the constant temperature zone (m); and g is the geothermal degree (m) [38]. Furthermore, to estimate the error of the calculation, the measured water temperature was compared with the results of Equation (4). Due to the relatively shallow occurrence of groundwater, heat losses associated with flow to the surface were omitted. In addition, it was assumed that the water, after passing through the heat pump, was cooled by 5 • C [55][56][57]. Based on the results obtained, a map was made for the estimated heating power of the intakes in the research area.
The physicochemical properties of the water were also analyzed in terms of the possibility of their use in low-temperature installations based on heat pumps. Preliminary verification of data obtained for the intakes suggest that in most cases, the physicochemical analyses were incomplete and related to only a few parameters important for the use of groundwater in heat pumps. Ultimately, it was decided to use the available analyses, however, they were expanded to include an additional 16 analyses of water samples, which were collected in the city of Rabka-Zdrój as part of the authors' own field tests ( Figure 5). Water samples from the wells were collected in accordance with the guidelines described in the PN-ISO 5667-11: 2017-10 Standard [58]. The samples were delivered to the laboratory within

Ground-Source Heat as a Lower Heat Source for Heat Pumps
For shallow depths, an analysis of the average thermal parameters of the rock mass, up to a depth of 2 and 10 m, was conducted, these depths usually being dedicated to horizontal or basket heat exchangers. The dominant sediments up to a depth of 2 m are dry clays and gravels with an average thermal conductivity of 0.4-0.64 W/m·K (Figure 6). The occurrence of waterlogged gravel was defined in several boreholes and their thermal conductivity was above 1.6 W/m·K. For most of the deposits up to 10 m, an increase in average thermal conductivity was noted. The highest values of λ: 1.9-2.2 W/m·K are most often associated with the occurrence of an irrigated sandstone layer. The unit heat output values q v ranged between 10 W/m 2 in the dry, unbound substrate and up to 40 W/m 2 in the substrate saturated with water (sand, gravel).
average thermal conductivity of 0.4-0.64 W/m•K (Figure 6). The occurrence of waterlogged gravel was defined in several boreholes and their thermal conductivity was above 1.6 W/m•K. For most of the deposits up to 10 m, an increase in average thermal conductivity was noted. The highest values of λ: 1.9-2.2 W/m•K are most often associated with the occurrence of an irrigated sandstone layer. The unit heat output values qv ranged between 10 W/m 2 in the dry, unbound substrate and up to 40 W/m 2 in the substrate saturated with water (sand, gravel).
For medium depths, the characteristics of the thermal parameters of rock mass to a depth of 30, 60, and 90 m are presented, which were used to assess the possibility of using vertical ground heat exchangers in the households. The weighted thermal conductivity coefficient varied from 1.14 to 2.38, 1.70 to 2.43, and 1.89 to 2.44 for the depths of 30, 60, and 90 m, respectively (Figure 7). High thermal conductivity of the rock substrate (≥2 W/m•K) characterized 41.0% of the boreholes to a depth of 30 m and about 80.0% of the boreholes to a depth of 60 and 90 m. The deposits where groundwater was dominant were fine-grained sandstones. The thicker the deposits, the higher the λ coefficient was achieved for a given depth range. In the southern part of the research area, the boreholes showed similar thermal parameters at depths of 60  For medium depths, the characteristics of the thermal parameters of rock mass to a depth of 30, 60, and 90 m are presented, which were used to assess the possibility of using vertical ground heat exchangers in the households. The weighted thermal conductivity coefficient varied from 1.14 to 2.

Groundwater as Lower Heat Source for Heat Pumps
Exploitable resources of groundwater intakes in the research area ranged from 0.01 to 8.1 m 3 /h ( Figure 8) [51]. In the analyzed intakes, 38% of them exceeded the capacity of 2 m 3 /s. The calculated water temperature in the intakes was within the range of 7.6-9.1 °C. The measured water temperature in the nine intakes ranged from 6.5 °C to 9.6 °C. The highest groundwater temperatures were recorded in the southeast and northwest of the area. The relative error between the measured and calculated water temperature values was 1.9-18.6%. The highest error (18.6%) was noted for the lowest measured temperature (6.5 °C). In the case of shallow groundwater, the heat output was affected more by the intake capacity than the water temperature. This can be seen in the example of measuring points, which was characterized by the highest temperature, but the lowest efficiency, so the estimated power obtained from the intake was much lower than in the case of high-performance and low-temperature intakes (south-eastern region). In over half of the analyzed intakes (58.6%), the calculated estimated heating power did not exceed 8 kW, and in 27.6% of the intakes, the heating power ranged from 20.0 to 47.1 kW.

Groundwater as Lower Heat Source for Heat Pumps
Exploitable resources of groundwater intakes in the research area ranged from 0.01 to 8.1 m 3 /h ( Figure 8) [51]. In the analyzed intakes, 38% of them exceeded the capacity of 2 m 3 /s. The calculated water temperature in the intakes was within the range of 7.6-9.1 • C. The measured water temperature in the nine intakes ranged from 6.5 • C to 9.6 • C. The highest groundwater temperatures were recorded in the southeast and northwest of the area. The relative error between the measured and calculated water temperature values was 1.9-18.6%. The highest error (18.6%) was noted for the lowest measured temperature (6.5 • C). In the case of shallow groundwater, the heat output was affected more by the intake capacity than the water temperature. This can be seen in the example of measuring points, which was characterized by the highest temperature, but the lowest efficiency, so the estimated power obtained from the intake was much lower than in the case of high-performance and low-temperature intakes (south-eastern region). In over half of the analyzed intakes (58.6%), the calculated estimated heating power did not exceed 8 kW, and in 27.6% of the intakes, the heating power ranged from 20.0 to 47.1 kW.
The assessment of the possible use of groundwater in low-temperature systems supported by water/water heat pumps was based on Table 2. The ranges of concentrations obtained for all physicochemical analyses are summarized in Table 3. Physicochemical analyses of groundwater from the PGI database were run between 1964-2014, 62% of which were performed after 2000. A physicochemical analysis of the water samples taken in the research area was conducted in 2018. The analysis showed that the permissible error according to the ion balance, up to 5%, was achieved for all samples taken. The pH of the analyzed intakes varied from acid to alkaline 6.9-9.1 (measure data) and 6.2-9.7 (PHS database). Obtained values of conductivity indicate the presence of both slightly The assessment of the possible use of groundwater in low-temperature systems supported by water/water heat pumps was based on Table 2. The ranges of concentrations obtained for all physicochemical analyses are summarized in Table 3. Physicochemical analyses of groundwater from the PGI database were run between 1964-2014, 62% of which were performed after 2000. A physicochemical analysis of the water samples taken in the research area was conducted in 2018. The analysis showed that the permissible error according to the ion balance, up to 5%, was achieved for all samples taken. The pH of the analyzed intakes varied from acid to alkaline 6.9-9.1 (measure data) and 6.2-9.7 (PHS database). Obtained values of conductivity indicate the presence of both slightly mineralized (213 μS/cm) and medium mineralized (1386 μS/cm) waters. Higher levels of chlorides and sulfates are observed in data obtained from the PGI control cards.   A groundwater quality assessment, taking into account heat pump manufacturers, was conducted for two variants. In the first case, compliance with the requirements of at least one manufacturer was considered, and in the second, compliance with the requirements of all manufacturers was observed. Interpretation of the collected results showed that the vast majority (95%) of groundwater met the requirements of at least one heat pump manufacturer (Figure 9). The most promising area in terms of the physicochemical parameters of the groundwater was the northern research area, where the majority of the intakes meeting the requirements of all heat pump manufacturers were located. Parameters that conducted for two variants. In the first case, compliance with the requirements of at least one manufacturer was considered, and in the second, compliance with the requirements of all manufacturers was observed. Interpretation of the collected results showed that the vast majority (95%) of groundwater met the requirements of at least one heat pump manufacturer (Figure 9). The most promising area in terms of the physicochemical parameters of the groundwater was the northern research area, where the majority of the intakes meeting the requirements of all heat pump manufacturers were located. Parameters that did not exceed the limit value of the heat pump manufacturers in any of the analyzed samples were NH4, NO3, SO4 −2 , and Cl − for the PGI database and CO2agg., NH4, Al, Cl − , SO −2 , and NO3 − from the authors' own measurements. It is worth noting that the heat exchangers used in the heat pumps are usually made of stainless steel, however, the heat pumps' individual components are copper-or nickel-brazed. Due to the fact that these materials have less resistance to potential threats, the final assessment should be guided by the values dedicated to them. In other cases where the requirements are not met, an intermediate exchanger should be used to protect the heat pump evaporator from damage.

Discussion
Based on the obtained results, it should be noted that the primary factor determining the recognition of potential low-temperature geothermal energy is to define and recognize the local geological and hydrogeological conditions [63,64]. By utilizing this information, it is possible to determine the thermal parameters of the rock mass, and thus determine the technical parameters of the installation, for example, the length and number of boreholes needed or the surface area of the ground heat exchanger.
In the case of horizontal ground heat exchangers (up to 10 m), solar radiation, precipitation, and snow cover play an important role. The temperature in the surface layers of the earth is similar to the average annual air temperature, however, this varies depending on the season. For the analyzed area, HGHEs should be placed up to 40 cm below the freezing zone [47,65] at a depth of about 1.8 m. The obtained results indicate relatively good thermal parameters in the rock mass, which can be used for very shallow geothermal systems. In addition, the presence of Quaternary waters in the central part of the research area results in an increase in the unit heat output of the soil. However, it should be

Discussion
Based on the obtained results, it should be noted that the primary factor determining the recognition of potential low-temperature geothermal energy is to define and recognize the local geological and hydrogeological conditions [63,64]. By utilizing this information, it is possible to determine the thermal parameters of the rock mass, and thus determine the technical parameters of the installation, for example, the length and number of boreholes needed or the surface area of the ground heat exchanger.
In the case of horizontal ground heat exchangers (up to 10 m), solar radiation, precipitation, and snow cover play an important role. The temperature in the surface layers of the earth is similar to the average annual air temperature, however, this varies depending on the season. For the analyzed area, HGHEs should be placed up to 40 cm below the freezing zone [47,65] at a depth of about 1.8 m. The obtained results indicate relatively good thermal parameters in the rock mass, which can be used for very shallow geothermal systems. In addition, the presence of Quaternary waters in the central part of the research area results in an increase in the unit heat output of the soil. However, it should be noted that the possibility of using HGHEs is also related to the size of the plot area available to the investor. For example, assuming that the demand for the heating power of a single-family house with an area of 130 m 2 is 8.5 kW (demand for a unit power of 65 W/m 2 ), where the compressor operation time is 1800 h/year, and the average unit heat output for the boreholes from the central research area is 20 W/m 2 , the HGHE area reaches a value of 425 m 2 .
Considering the presented research results for vertical ground heat exchangers, a significant part of the research area shows a high potential for the use of low-temperature systems. This is due to, though not limited to, the occurrence of irrigated sandstone, gravel, slate, and sometimes marl deposits. In addition, with this type of installation, the energy properties of the soil changes with the depth and Energies 2020, 13, 5010 15 of 20 the influence of the geothermal heat flux density is observed. Moreover, the area of research is located in a zone that has favorable conditions for the use of ground source heat pumps, especially for poorly insulated households that utilize the heat pumps extensively [66]. Using the determined average unit heat output (42 W/m) and the above-mentioned assumptions, a household in the research area would require a minimum of three wells with a depth of 67 m each. In terms of the thermal parameters of the borehole, the difference between the depth of 60 and 90 m is relatively small. The thermal conductivity of the rock mass for depths of 90 m were more favorable by 1.2 to 4.9% when compared to a depth of 60 m. This is mainly due to the presence of homogeneous rock masses. One should also note that the real local lithological conditions and the thermal conductivity coefficient could be different from bibliographic data [67][68][69]. Therefore, in the selected research area, a thermal response test (TRT) should be carried out to verify the bibliographic data for the thermal conductivity of the rock mass.
Based on the analyses conducted for groundwater intakes, it can be determined that there is an average potential for the use of low-temperature systems supported by water/water heat pumps. Groundwater intakes in this region are characterized by relatively low efficiency, and those with a capacity above 2 m 3 /h can be considered as promising. According to the literature, it is assumed that for every 5 kW of heating power, the intake capacity is 1.2-1.5 m 3 /h [70,71]. It should also be emphasized that determining groundwater temperature based on empirical formulas and other available data may be subject to errors that are difficult to estimate [50]. However, for the intakes of the household above-mentioned, the demand for heating power (8.5 kW) would cover about 41.4% of the analyzed intakes. The results of the assessments obtained from the physicochemical analyses show that the research area has a high potential for the use of groundwater as a lower heat source for the heat pumps.
According to the information contained in the Study of Conditions and Directions of Spatial Development, it follows that the development and modernization of heating systems should lead to the reduction of air polluting emissions [72]. The use of low-temperature systems based on heat pumps meets the above criterion. In addition, a mining area was designated due to the presence of healing waters in the central part of the analyzed region. The key to protecting heat pump water/water is to prohibit the discharge of sewage directly into the aquifer, and untreated wastewater into the surface waters and grounds of the area [73]. However, it is not prohibited to use the technical infrastructure for water and heat supply.
It should be emphasized that the bioclimatic factors occurring in a given region are widely used in health resort treatments and constitute one of the forms of comprehensive treatment for diseases, rehabilitation, and preventive care. However, in the area of the present study, the problem of carcinogenic benzo(a)pyrene (B(a)P) and particulate matter PM10 exceeding the permissible concentrations of the pollutants was identified. These pollutants are classified as being harmful to health and human life [74][75][76][77]. Moreover, the data obtained from the measurement station in Rabka-Zdrój until August 2020 also indicate a problem with air quality. The limit value of the daily average concentration of PM10 was exceeded for 20 days during the months from January to March. Moreover, the average monthly concentration of PM10 ranged from 13 to 51 µg/m 3 [78]. The use of ground heat and groundwater as a lower heat source for heat pumps is more ecologically friendly compared to coal-fired boilers [57]. By replacing a coal-fired boiler with a heat pump, the emission of pollutants such as B(a)P and dust are reduced by over 90% [3]. Therefore, the implementation of heating systems based on heat pumps in the health resort area would contribute to the improvement of air quality [79]. Moreover, the installation of low-enthalpy geothermal systems in individual-space heating is increasing worldwide [80,81]. The important advantage of heat pumps is the low operating costs in comparison with fossil fuels [82]. However, considering the large initial costs of the systems, the high price of electricity and the low price of natural gas, heat pump installations are not profitable without public subsidies [65]. In order to adopt heat pump installations on a large scale in Poland, several steps should be taken. The introduction of special electricity tariffs dedicated to the use of heat pumps, the creation of special research and development programs for this technology, the dissemination of general knowledge in the field of shallow geothermal energy, and the development of Polish standards for heat pump systems.

Conclusions
The article has provided an analysis of the potential of low enthalpy geothermal energy and the possibility of using solutions based on groundwater and obtaining heat from the ground. It should be noted that Polish health resorts are an integral part of the healthcare system and are eagerly visited by tourists and patients. Therefore, care for the environment of health resort areas should be a priority for both local and national authorities. Additionally, no extensive measures have been taken to assess the use of renewable energy sources in domestic households located in the area of the present study.
The performed spatial analysis clearly shows the distribution of thermal parameters of the rock mass within the analyzed boreholes, thus facilitating the selection of appropriate technical parameters for specific cases.
Maps have been produced based on the estimated heating power of the groundwater, which provide users with a quick and easy assessment of the heating power requirements of the heat pump-based systems. The rock mass showed high thermal conductivity for 41% of the boreholes up to a depth of 30 m and about 80% of the boreholes up to depths of 60 and 90 m. The estimated heating power was higher than 8 kW for about 40% of the analyzed intakes. Almost all of the studied water samples from the research area met the requirements of at least one heat pump manufacturer. In addition, the results of the physicochemical analyses showed that the area of Rabka-Zdrój has promise in terms of the possibility of using groundwater as a lower heat source for the heat pumps.
The present study into using shallow geothermal energy resources in the region presented in the article is one of the first. The results presented in the article can act as a reference tool for specialists, local government, and the community. Producing maps of low-temperature geothermal potential is not only a step toward local sustainable development, but also contributes to greater interest and use of low-temperature geothermal energy in domestic households. The implementation of solutions based on renewable energy sources is an opportunity to improve air quality as well as increase the environmental value of the health resort.