Energetic and Ecologic Heat Pumps Evaluation in Poland

: The purpose of the work is a comparison of indirect carbon dioxide emissions between the di ﬀ erent heat pump types that operate in Polish climate conditions. The analysis embraces air-to-water heat pump, ground-to-water heat pump, water-to-water heat pump and water-to-water heat pump with separating heat exchanger in the selected towns one in each climatic zone in the country. The study starts from determining seasonal coe ﬃ cient of performance in each location using heating degree days to estimate seasonal heat demand. Seasonal coe ﬃ cient of performance values enable an assessment which kind of heat pump meets the European Union requirements in every location. Eventually, indirect CO 2 emissions that is caused by electrical energy production, are estimated for every heat pump in each location. Ground-to-water heat pump and water-to-water heat pump satisfy these requirements in each climatic zone in Poland. Air-to-water heat pump would be an energetic and ecological viable on a condition that substantial changes were done in Polish electrical energy mix.


Introduction
On 8 October 2018, the Intergovernmental Panel on Climate Change (IPCC) published a report on the effects of global warming 1.5 • C above pre-industrial levels. According to its content, if people as an international community do not engage in the activities aimed at slowing the rise of the average temperature on Earth, we will, among other things, lead to undesired profound transformation of the ecosystem, intensification of dangerous weather phenomena (the hurricanes, droughts, floods) and in the long run to destruction of the world we currently live in [1]. Therefore, the countries and international organizations carry out the actions aimed at curbing the phenomenon of global warming. One of the countermeasures taken by the European Union (EU) is the pursuit of energy production in a climatic neutral way [2]. One proposal is to replace the high-emission heat sources used in individual central heating systems by heat pumps. Heat pumps do not directly emit any pollutants into the atmosphere but indirect emissions arise during the production of energy necessary for their driving [3]. Due to their growing popularity, variety of heat sources or operating parameters, the European Commission has introduced regulations specifying the minimum value of the seasonal performance coefficient (SCOP) [4].
Lucia et al. [5] reviewed literature on heat pump from the technological and thermodynamic aspects. Heat pumps that use the ground as a heat source are divided into four categories. GWHPs, ground-water heat pump systems named also open-loop systems, are vertical open systems because water is drawn from a well and discharged to an appropriate outlet. GWHPs were firstly installed in 1940s. GCHP, ground-coupled heat pump systems named also closed-loop GSHP systems. These systems were After a replacement an existing coal-fired furnace with a heat pump the emissions will be avoided as follows-1709 t CO 2 , 2.56 kg SO 2 , 1.62 kg CO, 2.99 kg NO x and 0.32 kg dust.
A heat pump filled with zeotropic refrigerant R407C was tested by Szreder [17]. It produces heat for space heating and domestic hot water in a detached house. There are two underfloor heating loops the first one on the ground floor and the second one on the first floor. Every loop (ground floor heating, the first floor heating and domestic hot water circuit) is equipped with a pump. Three vertical U-tube ground heat exchangers, with separate glycol pomp for each, are used as a low temperature heat source. The more important outcomes are as follows. Circulation pump in the underfloor heating loops should work at the highest (2nd) gear, which causes larger mass flow rate and more turbulent flow, hence temperature in the condenser drops. As a result COP value increases. COP value greatly depends on glycol temperature in the underground loops. Applying inverter for driving a spiral compressor does not bring considerable economic benefits.
Pater [18] described the seven-year experiment, in a residential building located near Kraków (Poland), with a system equipped with gas condensing boiler, biomass boiler, solar collectors and ground source heat pump as the alternative heat generators. There were three circuits-central heating, domestic hot water and cooling. Annual seasonal performance factor for heat pump was between 3.5 and 3.8.
Szulgowska-Zgrzywa and Fidorów-Kaprawy [19] analyzed ground source heat pump operation in an office building located inŚwidnica (Poland) in the third climatic zone. The heat pump unit was either a heat generator or cooling device. The central heating system was separated from the heat pump unit with a buffer tank. The solar collectors either were heating the buffer tank in the heating season or were regenerating ground around one of the boreholes when heating was turned off. Other borehole was regenerating when the heat pump was working as an air cooler. There were a total of 5 boreholes with a depth of 78 m. Seasonal coefficient of performance (SCOP) was 4.05 without the buffer tank. Buffer tank which maintained constant temperature 40 • C decreases SCOP by 20%.
Gajewski et al. [20] analyzed the operation of air, ground and groundwater heat pumps located in Bialystok (Poland). The climatic conditions were determined using ten-year mean temperature. In the case only groundwater heat pump emits less carbon dioxide than condensing gas boiler. Heat pumps supplied with heat from the air or ground would emit more carbon dioxide than condensing gas boiler.

The Purpose of the Research
A heat pump is not always more ecological than a condensing gas boiler, especially in countries where the majority of electricity production is done in the fossil-fuels plants. Heat pump usage should be promoted if their operation is energetically and ecologically profitable. Seasonal coefficient of performance (SCOP) is an averaged ratio which inform about energetic profitability. Coefficient of performance (COP) values that are the components of SCOP depend on temperature values in a heat source and heat sink and thermodynamic properties of a refrigerant. Carbon dioxide emissions may be applied as a factor of ecologic profitability. If total CO 2 emissions caused by a heat pump are lower than by any fossil fuel boiler then heat pump may be considered as a "green" heat source.
Temperature values in the heat source depend on the climate conditions, its thermodynamic properties or geologic structure and heat exchange conditions. The most significant and rapid changes occur if external air is the heat source. Ground temperature changes gradually, while ground water temperature may be thought as constant.
Temperature conditions in the heat sink depend on flow temperature in a heating system that is a function of designed flow and return temperature in the heating system and external air temperature. External air temperature affects the flow temperature through a heating curve in a weather compensator.
Heat pumps performance is analyzed as either the device operation or the part of the system. The paper authors focus on the latter option. As bituminous coal and lignite combustion forms the overwhelming majority in the Polish energetic mix the aim of the pursued investigations is a comparison of carbon dioxide emissions in all the climatic zones in Poland. Hence, the shares of fossil fuels and the electricity production efficiencies as well as the efficiency of low voltage electricity transfer are the factors determined the Polish energetic mix. Moreover, an impact of the heating system on heat sink temperature is included in the analysis, which is made by taking a heating curve into account. Such an innovative approach is the goal of undertaken research. The emissions are estimated for air-to-water heat pump (ASHP), ground-to-water heat pump (GSHP) and water-to-water heat pump (WSHP). The purpose of the work is a comparison between indirect emissions for various types of heat pumps in Polish climate conditions.

Methods and Assumptions
An area of research is to determine the net values of the seasonal coefficient of performance (SCOP net ) and carbon dioxide emissions of fixed-speed heat pumps that are the heat generators in a single-family house. 4 variants of available solutions on the market are analyzed-ASHP, GSHP, WSHP and water-to-water heat pump whose operation is preceded by a separation heat exchanger (WSHP+SHE). There are compared the heat pumps that are available on the market on the day of the test, they are enumerated in Table 1. Heat pump is the only heat generator in a single-family house that total heated area is 280 m 2 . Insulation of walls, as well as window and door joinery complies with the requirements corresponding to the requirements for thermal insulation of partitions in buildings commissioned from 1 January 2017 specified in the applicable legal regulations [21][22][23]. The total designed heat load of the building is 14 kW. The house is located in Bialystok. Due to the assumption that the presented analysis compares the operation of the same devices operating in the different climatic conditions, the same design heat load 14 kW in each selected location is assumed. Five selected locations, one of each of the Polish climatic zones (cf. Figure 1), were analyzed-Szczecin, Poznań, Łódź, Białystok, Zakopane (Table 2).

Computation of Seasonal Coefficient of Performance
The analysis is started from determining coefficient of performance value at each operating temperature (COP bin ). This value depends on both brine and outside temperature. Outside temperature is assumed to be equal to the mean value in a given month. This is due to the fact that heating degree days (HDD) represent different climatic conditions. Because COP is non-linear function of temperature and mean temperature value could overestimate the final test results, the second calculation option is considered. In this option COP depends on average minimal outside temperature that is the difference between the average monthly temperature and half of the monthly amplitude. The options discussed above are marked as "Variant I" and "Variant II" respectively.
Then the temperature values in every lower heat sources are determined. It should be emphasized that the higher reservoir temperature is, the lower electric energy for heat pump driving is demanded, which translates directly into COP bin values. Eventually, the resultant SCOP affects the final ecological and economic viability. The temperature values of the lower source and their variability depend on the properties of a source and have an effect on a device applied in the installation.
When using ASHP, outside air is the lower temperature heat source. This means that its value is subjected to continuous fluctuations depending on the changes in the weather conditions. This results in a large diversity of COP bin values. As temperature decreases, amount of heat from the source diminishes while necessary work rises. This problem is substantially reduced when a heating system is supplied from either of the other types of pumps. The changing external conditions have little effect on mean ground temperature because ground below the shallow zone has constant temperature. They have no impact on temperature in groundwater basin.
In the case of the GSHP, it is necessary to turn attention to the depth to which the heat source extends. It is equally important to determine the geological structure and geographical location. These pieces of information are necessary to correctly determine the temperature profile of the heat source. In the analysis, the temperature profile is determined using the Baggs formula modified by a team led by Olesiewicz-Popiel [25] for the northern hemisphere: where τ is a day number in the year, z is the depth (m), ∆t m is a difference between the ground temperature below the shallow zone and the average annual air temperature ( • C), k v is a vegetation coefficient, A s is the amplitude of annual air temperature in every climatic zone respectively ( • C), a is the soil thermal diffusivity (m 2 /s), τ o is the phase of the air temperature wave (day), t a is an average annual air temperature in each climatic zone, respectively ( • C).
Equation (1) is averaged over time and depth to determining the mean annual ground temperature in each zone. A s and t a are taken from columns 4 and 6 in Table 2 respectively. Table 3 shows other data for the mean ground temperature computation. In the case of an installation equipped with a WSHP, the groundwater basin is the lower heat source. When determining the value of the source temperature, the Kowalski [26] relation is used: where t a is average annual air temperature in each climatic zone respectively ( • C), A is a correction factor due to true altitude, g g is a geothermal gradient ( • C/m), H is a depth of a groundwater basin level (m), h is the depth of the shallow zone (m). The datum set and the results for further modelling are shown in Table 4.  (a) According to Reference [26]. (b) According to Reference [27]. (c) According to Reference [28].
When using a water-to-water pump, in addition to meeting the minimum and maximum source temperature requirements, it is necessary to compare the chemical composition compliance of the groundwater used with the requirements set by the device manufacturer. If these are not met, a separating heat exchanger (SHE) in the primary circuit should precede the heat pump in such a way the groundwater transfers heat to the brine, which causes the brine temperature decrease by 2 degrees. If temperature dropped below the minimal acceptable value by the manufacturer, then other after SHE GSHP would be installed. In the investigated cases, the WSHP + SHE units meet the manufacturer's requirements in each of the considered locations.
Knowing the values of the external temperature and lower heat sources, it is possible to move on to determining the most appropriate operating parameters of the installation. Due to the fact that the work of each device is controlled by the same weather compensator model, it is possible to select one heating curve, among those available in the control system used, along which each heat pump work is controlled. Heating curve, shown in Figure 2, in the algorithm is approximated by a square polynomial.
Energies 2020, 13, x FOR PEER REVIEW 8 of 18 the work of each device is controlled by the same weather compensator model, it is possible to select one heating curve, among those available in the control system used, along which each heat pump work is controlled. Heating curve, shown in Figure 2, in the algorithm is approximated by a square polynomial.
The next step is to determine the partial heat demand, which is done in accordance with the standard PN-EN 14825:2019-03 [30]: where Φi is the total design heat load obtained using an algorithm from the standard PN-EN 12831-1:2017-08 [31], ti is the internal design temperature that equals 20 °C, tj is the external air temperature, te is the external design temperature. The performance drop caused by the on-off cycles was not included in the calculations. The seasonal coefficient of performance (SCOPnet) value is determined based on the former standard [30]: where hj is the number of bin hours occurring at external temperature tj in the heating season, COPbin(tj) is COP value of the unit at external temperature tj, j is a number of a temperature value (the temperature values are sorted in ascending order).

Determinating the Value of Pollutant Emissions
In a central heating installation, pollutant emissions occur as a result of burning a fuel to generate heat in its needs (direct emission) and the fossil fuels combustion at the power plants to produce electricity for the heat source driving (indirect emission). Heat pumps cause no direct emission, because the heat is provided from zero-emission sources. The indirect emission value was determined on the basis of the Equation [32]: The next step is to determine the partial heat demand, which is done in accordance with the standard PN-EN 14825:2019-03 [30]: where Φ i is the total design heat load obtained using an algorithm from the standard PN-EN 12831-1:2017-08 [31], t i is the internal design temperature that equals 20 • C, t j is the external air temperature, t e is the external design temperature. The performance drop caused by the on-off cycles was not included in the calculations. The seasonal coefficient of performance (SCOP net ) value is determined based on the former standard [30]: where h j is the number of bin hours occurring at external temperature t j in the heating season, COP bin (t j ) is COP value of the unit at external temperature t j , j is a number of a temperature value (the temperature values are sorted in ascending order).

Determinating the Value of Pollutant Emissions
In a central heating installation, pollutant emissions occur as a result of burning a fuel to generate heat in its needs (direct emission) and the fossil fuels combustion at the power plants to produce electricity for the heat source driving (indirect emission). Heat pumps cause no direct emission, because the heat is provided from zero-emission sources. The indirect emission value was determined on the basis of the Equation [32]: where β ag is an aggregate pollutant generation factor which accounts direct and indirect emissions amid electrical energy production and considers shares of all the fuels in Polish energy market (cf. PGE 2020). Information on the share of individual sources used for the production of electricity used by the heat pump is consistent with the data on the fuel structure provided by PGE [33] and is presented in Table 5. Ciura [34] estimated the low voltage electricity transfer efficiency at 79.2% Table 5. The shares of the energy sources in electricity production in Poland [33].

Energy Source
Generation Based on the Equation (5) it can be clearly seen that the environmental impact of a heat pump operation depends not only on the one's energy conversion ratio but also on the energy market in a given country. The high percentage of electricity generation from renewable, zero-emission energy resources in the overall balance results in lower indirect emissions.

Results and Discussion
The graph of COP bin values for the 4 analyzed heat pumps is presented in Figure 3 which shows that regardless of the type of device used and its location, the COP bin value increases as the outside air temperature rises. This is due to the fact that with the higher temperature difference between hotter and colder sources the work needed to driving a device grows. COP bin value for ASHP is effected directly with the outside temperature changes. An impact of these changes on ground temperature is significantly reduced but not completely lost, which stems from Equation (1). Groundwater temperature is independent of outdoor temperature changes. COP values drop steady for GSHP, WHSP and WSHP + SHE as outside temperature goes down. In the case of ASHP, there are observed sharp bends as temperature falls to about 7 • C and then marked reduction in the plots for zones I to V respectively. In the case of the ASHP pump, in each location, a sudden decrease in the COP value is noticeable when the outside air temperature drops below 7 • C. It is related to the evaporation pressure decrease which increases the pressure difference between the front and rear of the compressor. Thus compressor electrical energy consumption rises. Simultaneously, outside temperature drop raises heat demand. Consequently, COP bin values rapidly decreases, which is illustrated in the instruction manual [29].   Only in the case of I climatic zone the GSHP is characterized by higher COPbin values than WSHP. It is observed at the lowest assumed operating temperature (in variant II) when COPbin values are 4.88 for the former and 4.79 for the latter. For the highest temperature value (in variant I), the reciprocal relation is opposite. Here the COPbin values are 5.41 and 5.53, respectively. Such a situation stems from the shallowest depth of the groundwater basin in Szczecin (cf. Table 4) and the highest temperature value (cf. Table 2), which results in the smallest temperature difference between groundwater and ground. In other locations this difference is greater, thus WSHP is more effective in comparison to GSHP in the full outside temperature range. It should be emphasized this superiority decreases with temperature drop. In the case of zone II, COPbin values at the lowest temperature in the second variant are for WSHP and for GSHP 4.83 and 4.75, respectively. Similarly, in zone III COPbin equals 4.74 and 4.66 and in zone IV the values are 4.56 and 4.4, respectively. The comparison of COPbin values between GSHP and WSHP+SHE is equally interesting. In the case of zones II-IV, these values almost coincide at the upper temperature limit. In the case of zones II and III the value difference is 0.03 (zone II-for GSHP 5.25, for WSHP 5.28, zone III-for GSHP 5.26, for WSHP 5.29), while in zone IV it is 0.05 (for GSHP 5.08, for WSHP 5.13). The relative shift between the graphs of GSHP and WSHP (zone I) or GSHP and WSHP + SHE (zones II-IV) results from the temperature difference between the ground and groundwater. In the first climatic zone, the ground temperature is lower than the groundwater temperature and simultaneously higher than the brine temperature after SHE, which causes the GSHP and WSHP curves intersect. In other cases, the temperature of the ground is lower than both the groundwater temperature and brine temperature Only in the case of I climatic zone the GSHP is characterized by higher COP bin values than WSHP. It is observed at the lowest assumed operating temperature (in variant II) when COP bin values are 4.88 for the former and 4.79 for the latter. For the highest temperature value (in variant I), the reciprocal relation is opposite. Here the COP bin values are 5.41 and 5.53, respectively. Such a situation stems from the shallowest depth of the groundwater basin in Szczecin (cf. Table 4) and the highest temperature value (cf. Table 2), which results in the smallest temperature difference between groundwater and ground. In other locations this difference is greater, thus WSHP is more effective in comparison to GSHP in the full outside temperature range. It should be emphasized this superiority decreases with temperature drop. In the case of zone II, COP bin values at the lowest temperature in the second variant are for WSHP and for GSHP 4.83 and 4.75, respectively. Similarly, in zone III COP bin equals 4.74 and 4.66 and in zone IV the values are 4.56 and 4.4, respectively. The comparison of COP bin values between GSHP and WSHP+SHE is equally interesting. In the case of zones II-IV, these values almost coincide at the upper temperature limit. In the case of zones II and III the value difference is 0.03 (zone II-for GSHP 5.25, for WSHP 5.28, zone III-for GSHP 5.26, for WSHP 5.29), while in zone IV it is 0.05 (for GSHP 5.08, for WSHP 5.13). The relative shift between the graphs of GSHP and WSHP (zone I) or GSHP and WSHP + SHE (zones II-IV) results from the temperature difference between the ground and groundwater. In the first climatic zone, the ground temperature is lower than the groundwater temperature and simultaneously higher than the brine temperature after SHE, which causes the GSHP and WSHP curves intersect. In other cases, the temperature of the ground is lower than both the groundwater temperature and brine temperature after SHE. This leads to the intersection of the GSHP and WSHP + SHE curves at the same temperatures as in zone I. Unlike is in the fifth climatic zone because of the great groundwater basin depth which results in considerable temperature difference between groundwater and ground. In the first variant COP bin values for GSHP (4.36), ASHP (4.35) are almost equal at the highest temperature. Slightly more difference is between WSHP (7.09) and WSHP + SHE (6.95).
SCOP net values in every climatic zone are summarized in Figure 4.
Energies 2020, 13, x FOR PEER REVIEW 12 of 18 after SHE. This leads to the intersection of the GSHP and WSHP + SHE curves at the same temperatures as in zone I. Unlike is in the fifth climatic zone because of the great groundwater basin depth which results in considerable temperature difference between groundwater and ground. In the first variant COPbin values for GSHP (4.36), ASHP (4.35) are almost equal at the highest temperature. Slightly more difference is between WSHP (7.09) and WSHP + SHE (6.95). SCOPnet values in every climatic zone are summarized in Figure 4. According to the recommendations issued by the European Commission on 1 March 2013, only those heat pumps that are energetically viable may be applied [35]. To satisfy this rule, in the conditions of the Polish energy mix, the SCOPnet value should be at least 3.5 [32]. Figure 4 shows ASHP does not meet this condition in any of the analyzed zones. Hence, devices of this type cannot be considered in Poland as alternative, environmentally friendly heat generators. Otherwise is in the cases of GSHP and WSHP which meet this condition in each climatic zone. However, there are visible differences in comparison to individual locations resulting from the difference in hydrological or geological conditions between individual regions.
SCOPnet values are higher for the I variant than in the II variant, which results directly from the higher values of the assumed external temperatures. The high SCOPnet values of GSHP or WSHP are also influenced by warmer temperature as well as greater temperature stability of the respective colder sources. In Szczecin SCOPnet values for GSHP and WSHP differ about 1%, while GSHP outstrips WSHP + SHE by 5% in both variants. In Poznań, Łódź and Białystok SCOPnet values for WSHP are higher than GSHP by 4-5% in the first variant and by 2-4% in the second one. Contrastively, GSHP is more effective than WSHP + SHE by 1-2% in both variants. In Zakopane WSHP and WSHP + SHE predominate GSHP by 52-65% in both variants.
Since, COP is inversely proportional to electrical energy demand and all the emissions are directly proportional to electricity usage, the results are in inverse proportion. In Szczecin electricity consumption and emissions for GSHP and WSHP are almost equal, while these values are greater by 6% for WSHP + SHE. In Poznań, Łódź and Białystok energy demand and emissions WSHP are less than GSHP by 3-5% in the first variant and by 2-4% in the second one. In contrast GSHP uses less energy and emits less pollutions by 1-2% than WSHP + SHE in both variants. In Zakopane WSHP and WSHP + SHE save 35-40% energy and equally reduce emissions with relation to GSHP.  According to the recommendations issued by the European Commission on 1 March 2013, only those heat pumps that are energetically viable may be applied [35]. To satisfy this rule, in the conditions of the Polish energy mix, the SCOP net value should be at least 3.5 [32]. Figure 4 shows ASHP does not meet this condition in any of the analyzed zones. Hence, devices of this type cannot be considered in Poland as alternative, environmentally friendly heat generators. Otherwise is in the cases of GSHP and WSHP which meet this condition in each climatic zone. However, there are visible differences in comparison to individual locations resulting from the difference in hydrological or geological conditions between individual regions.
SCOP net values are higher for the I variant than in the II variant, which results directly from the higher values of the assumed external temperatures. The high SCOP net values of GSHP or WSHP are also influenced by warmer temperature as well as greater temperature stability of the respective colder sources. In Szczecin SCOP net values for GSHP and WSHP differ about 1%, while GSHP outstrips WSHP + SHE by 5% in both variants. In Poznań, Łódź and Białystok SCOP net values for WSHP are higher than GSHP by 4-5% in the first variant and by 2-4% in the second one. Contrastively, GSHP is more effective than WSHP + SHE by 1-2% in both variants. In Zakopane WSHP and WSHP + SHE predominate GSHP by 52-65% in both variants.
Since, COP is inversely proportional to electrical energy demand and all the emissions are directly proportional to electricity usage, the results are in inverse proportion. In Szczecin electricity consumption and emissions for GSHP and WSHP are almost equal, while these values are greater by 6% for WSHP + SHE. In Poznań, Łódź and Białystok energy demand and emissions WSHP are less than GSHP by 3-5% in the first variant and by 2-4% in the second one. In contrast GSHP uses less energy and emits less pollutions by 1-2% than WSHP + SHE in both variants. In Zakopane WSHP and WSHP + SHE save 35-40% energy and equally reduce emissions with relation to GSHP.
Since, COP values are determined using the algorithm which assumes constant monthly outdoor air temperature and COP is not a linear function of temperature, the SCOP net values might be overestimated. Especially, when outside temperature drops up to −30 • C, while the lowest mean value is 2.76 in zone IV and 2.81 in zone V.
The differences between the results obtained in the options I and II are also visible in the case of electricity demand shown in Figure 5. Variant II shows clearly higher energy consumption. It is caused by lower assumed values of external temperatures and thus the necessity of incurring higher energy outlays in order to ensure an adequate amount of heat.
Energies 2020, 13, x FOR PEER REVIEW 13 of 18 Since, COP values are determined using the algorithm which assumes constant monthly outdoor air temperature and COP is not a linear function of temperature, the SCOPnet values might be overestimated. Especially, when outside temperature drops up to −30 °C, while the lowest mean value is 2.76 in zone IV and 2.81 in zone V.
The differences between the results obtained in the options I and II are also visible in the case of electricity demand shown in Figure 5. Variant II shows clearly higher energy consumption. It is caused by lower assumed values of external temperatures and thus the necessity of incurring higher energy outlays in order to ensure an adequate amount of heat. This translates directly into emissions presented in Figures 6-9. In the case of ASHP, where the demand for electricity is the highest, we observe the highest pollutions emissions in each zone. Heat pump that is supplied from a constant temperature heat source causes less indirect pollutions emissions. Regardless of a location and technological variant of an applied device, the by-product that is emitted in the largest amount is carbon dioxide. This is directly related to the fact that the current energy policy is based primarily on natural fossil raw materials, which are mainly composed of compounds containing carbon atoms. In the case of the remaining pollutants analyzed in this study, that is, sulfur dioxide and nitrogen dioxide as well as particulate matter (PM), the level of indirect seasonal emissions is similar.  Figure 5. A summary of the energy demand for each heat pump in the variants I and II.
This translates directly into emissions presented in Figures 6-9. In the case of ASHP, where the demand for electricity is the highest, we observe the highest pollutions emissions in each zone. Heat pump that is supplied from a constant temperature heat source causes less indirect pollutions emissions. Regardless of a location and technological variant of an applied device, the by-product that is emitted in the largest amount is carbon dioxide. This is directly related to the fact that the current energy policy is based primarily on natural fossil raw materials, which are mainly composed of compounds containing carbon atoms. In the case of the remaining pollutants analyzed in this study, that is, sulfur dioxide and nitrogen dioxide as well as particulate matter (PM), the level of indirect seasonal emissions is similar.
It should be emphasized that it is always necessary to perform an analysis of running and investment costs. The analysis should also include thorough tests determining either physicochemical or geological parameters of the proposed colder heat source, since this affects the final value of SCOP and thus the overall costs of the installation.        It should be emphasized that it is always necessary to perform an analysis of running and investment costs. The analysis should also include thorough tests determining either  It should be emphasized that it is always necessary to perform an analysis of running and investment costs. The analysis should also include thorough tests determining either

Conclusions
Considering the value of the seasonal energy efficiency ratio and emission levels, is is justified to state that in the Polish climate and with the energy mix conditions, both ground-to-water and water-to-water heat pumps are energetically viable and meet the conditions of an ecological heat source in accordance with the applicable European Union regulations; the wider application of these heat pumps could reduce the emissions caused by space heating. In the case of air-to-water heat pumps, it is necessary to introduce deep changes in the energy market in advance so that this device could be considered an ecological alternative. Moreover, the air-to-water heat pump is an energetically nonprofitable heat generator in the Polish energy mix. The final decision on choosing the type of heat pump should be preceded by an in-depth analysis of local temperature, geological, hydrological conditions and economic costs.
Author Contributions: S.S. and A.G. developed a calculation algorithm, S.S. computed SCOP, A.G. computed emissions, S.S. and A.G. analyzed the data and wrote the paper. All authors have read and agreed to the published version of the manuscript.