Methodology for Quantifying the Energy Saving Potentials Combining Building Retroﬁtting, Solar Thermal Energy and Geothermal Resources

: New technological, societal and legislative developments are necessary to support transitions to low-carbon energy systems. The building sector is responsible for almost 36% of the global ﬁnal energy and 40% of CO 2 emissions, so this sector has high potential to contribute to the expansion of positive energy districts. With this aim, a new digital Geographic Information System (GIS) platform has been developed to quantify the energy savings obtained through the implementation of refurbishment measures in residential buildings, including solar thermal collectors and geothermal technologies and assuming the postal district as the representative unit for the territory. Solar resources have been estimated from recently updated solar irradiation maps, whereas geothermal resources have been estimated from geological maps. Urbanistic data have been estimated from o ﬃ cial cadastre databases. For representative buildings, the annual energy demand and savings are obtained and compared with reference buildings, both for heating and cooling. The GIS platform provides information on average results for each postal district, as well as estimates for buildings with particular parameters. The methodology has been applied to the Asturian region, an area of about 10,600 km 2 on the Cantabrian coast of Spain, with complex orography and scattered population, qualiﬁed as a region in energy transition. High rehabilitation potentials have been achieved for buildings constructed before the implementation of the Spanish Technical Building Code of 2006, being higher for isolated houses than for collective buildings. Some examples of results are introduced in speciﬁc localities of di ﬀ erent climatic zones.


Introduction
The energy transformation towards a more efficient and less polluting system represents a global problem that must be addressed by all sectors of society. National, regional and local governments must work together involving all the stakeholders, public and private. The Intergovernmental Panel on Climate Change (IPCC) highlighted in its 2018 report [1] the strong impact produced by global The energy potential of the biomass obtained from the urban greenery maintenance has been evaluated and has been considered as a low-range renewable resource [36]. More recently, a new GIS model has been developed to face heat waves by using pavement watering for urban cooling [37].
Moreover, these techniques have been used to obtain information about the boundary conditions of the urban environment. For example, Viana-Fonts et al. developed shadow cast profiles of buildings in urban areas from cadastral cartography and LiDAR altimetric data [38]. Liang et al. developed an interactive GIS tool for sky, tree and building view factor estimation from street view photographs [39].
To summarize, according to previous studies, the flexibility of GIS tools enables the performance of a wide variety of urban analysis.
In this context, this paper presents a new geographic information system platform developed to spatially analyse relevant urban features, such as building density and type, geothermal resources and solar potential production. Furthermore, potential energy savings on the existing building stock, considering all these factors, have been calculated. This platform relies upon GIS technology, is an initial approach to urban energy analysis and contributes to the reinforcement of European energy and climate policies oriented to a fundamental transformation of the energy system. Proposed cost-effective actions improve energy efficiency, promote centralised solutions based on low-carbon energy innovations and boost technologies to mitigate global warming. An initial case study has been applied to the Principality of Asturias, a region of more than 10,000 km 2 on the northern coast of Spain, with complex orography and scattered population. This region is characterized by a temperate climate based on the Köppen Geiger classification, with mild temperatures registered in both winter and summer and abundant rainfall throughout the year. There are different climate zones depending on the precipitation patterns, ranging from oceanic (Cfb) to Mediterranean oceanic (Csb). These climatic conditions lead to high heating loads and low cooling loads. These actions, in terms of reducing the carbon cycle for ongoing energy transition processes, guarantee a sustainable socio-economic development of the area, and incentivise a municipality scale approach to the implementation of renewable energies through efficient district systems.

Materials and Methods
To cope with the increase in the energy consumption in buildings, several projects have been developed during recent decades with the aim of helping users, local administrations and stakeholders to make the better decisions to improve and optimize the energy efficiency of the building stock [40].
In this line, the Spanish research project RehabilitaGeoSol (RGS) [41] has developed a citizen-oriented platform (RGS platform) that suggests to the parties involved the energy savings obtained through the implementation of retrofit interventions, solar thermal panels and geothermal technologies in residential buildings.
A building expert should choose between several retrofitting actions, reconciling energy, environmental, legal and social factors to achieve a good compromise solution that meets the final needs of occupants. The search for a reasonable solution can be carried out through an energy analysis of the building and through the development of various situations predefined by a building-physics expert, mainly through simulations [42]. In addition, a second approach which incorporates decision aid techniques such as multi-criteria analysis [43], multi-objective optimization [44], generally combined with simulations to help reach a final decision [45], can be considered.
The methodology used to create the retrofit databases of the RGS project has been applied to the Principality of Asturias. Up to 800 postal districts have been analysed, showing the interest of combining building retrofitting with solar and geothermal strategies.
Annual and seasonal thermal loads for each refurbishment option proposed have been calculated in order to provide the input data to the platform. The final outputs are the energy-savings reached by the implementation of retrofit measures, solar thermal panels and geothermal technologies in residential buildings.

Solar Resource
In the case of the Asturian region, there have been solar maps available since 2009 [46], which provide monthly and annual average global solar irradiation on a horizontal surface. These maps were created using indirect estimations of irradiation from air temperatures [47] and GIS techniques, and have been recently updated taking climatic data from 2007 to 2017 into account. The spatial resolution of the original maps is 25 × 25 m, but they have been converted to postal district level using ArcGIS software for spatial integration. Figure 1 shows the variations of annual average irradiation in different postal districts. Similar maps were generated for maximum, mean and minimum values of air temperature.

Solar Resource
In the case of the Asturian region, there have been solar maps available since 2009 [46], which provide monthly and annual average global solar irradiation on a horizontal surface. These maps were created using indirect estimations of irradiation from air temperatures [47] and GIS techniques, and have been recently updated taking climatic data from 2007 to 2017 into account. The spatial resolution of the original maps is 25 × 25 m, but they have been converted to postal district level using ArcGIS software for spatial integration. Figure 1 shows the variations of annual average irradiation in different postal districts. Similar maps were generated for maximum, mean and minimum values of air temperature.

Geothermal Resource
The geothermal resource is based on the National Geological Map (MAGNA), prepared by the Geological and Mining Institute of Spain (IGME) [48] between 1972 and 2003. This map consists of sheets that show geological data (geological sections, stratigraphic profiles, boreholes, etc.) with structural and hydrogeological diagrams of the national territory.
Each material is assigned some conductivity according to tabulated data by the Institute for the Diversification and Saving of Energy of Spain (IDAE) [49]. Table 1 shows some representative examples of rocks and minerals existing in the Principality of Asturias [48] and their main characteristics, among which different thermal conductivities stand out. In addition, bibliographic data are contrasted with the Enhanced Geothermal Response Test, so that the basic properties of geothermal resources have been characterized for most materials in each

Geothermal Resource
The geothermal resource is based on the National Geological Map (MAGNA), prepared by the Geological and Mining Institute of Spain (IGME) [48] between 1972 and 2003. This map consists of sheets that show geological data (geological sections, stratigraphic profiles, boreholes, etc.) with structural and hydrogeological diagrams of the national territory.
Each material is assigned some conductivity according to tabulated data by the Institute for the Diversification and Saving of Energy of Spain (IDAE) [49]. Table 1 shows some representative examples of rocks and minerals existing in the Principality of Asturias [48] and their main characteristics, among which different thermal conductivities stand out. In addition, bibliographic data are contrasted with the Enhanced Geothermal Response Test, so that the basic properties of geothermal resources have been characterized for most materials in each postal district ( Figure 2). In the north-western regions, there is a high concentration of zones with values of thermal conductivity above 3 W/(m·K). In the eastern regions, the highest concentration takes place for conductivities between 2 and 3 W/(m·K). In the lower central regions, most of the zones have thermal conductivities below 2.6 W/(m·K), while in the upper central regions, most of them oscillate between 2.2 and 3 W/m K. Despite these trends, there is a great variation in conductivities throughout the Principality of Asturias.
Energies 2020, 13, x FOR PEER REVIEW 5 of 27 postal district ( Figure 2). In the north-western regions, there is a high concentration of zones with values of thermal conductivity above 3 W/(m·K). In the eastern regions, the highest concentration takes place for conductivities between 2 and 3 W/(m·K). In the lower central regions, most of the zones have thermal conductivities below 2.6 W/(m·K), while in the upper central regions, most of them oscillate between 2.2 and 3 W/m K. Despite these trends, there is a great variation in conductivities throughout the Principality of Asturias.

Urban Data
Official cadastre databases have been carefully handled and buildings in each postal district have been classified in the categories of Table 2. Notice that considered time periods correspond to different technical energy requirements for buildings in Spain, namely, the Spanish Basic Building Standard -NBE-CT-79 [50], the Spanish Technical Building Code-CTE [51], which entered into force at the end of 2007, and its amended version at the end of 2013. Because of this classification, statistics have been obtained for each postal district, indicating the number of buildings of each category, which is the key information to evaluate the expected energy savings if retrofitting measures are applied.

Building Energy Demand
The necessary requirements to identify the most suitable refurbishment strategies to reduce the energy consumption of a building are local climates and constructive information [52]. In the first place, it is necessary to count on the availability of a representative climate data [53] of a locality. Second, it is necessary to know the characteristics of the existing building façades [54], shading devices, air ventilation [55], fenestration types and the user occupancy [56]. The study of the potential of climatic variables helps to identify the group of passive and active techniques that, when coupled, give the best result [57]. In addition to the above, the legal restrictions and the requirements of local legislation must be considered.
A dynamic simulation environment has been developed to calculate the building energy demands provided by this new digital platform [58]. This method requires an exhaustive definition of the building cases and the boundary conditions, the execution of a multi-parametric study and a

Urban Data
Official cadastre databases have been carefully handled and buildings in each postal district have been classified in the categories of Table 2. Notice that considered time periods correspond to different technical energy requirements for buildings in Spain, namely, the Spanish Basic Building Standard-NBE-CT-79 [50], the Spanish Technical Building Code-CTE [51], which entered into force at the end of 2007, and its amended version at the end of 2013. Because of this classification, statistics have been obtained for each postal district, indicating the number of buildings of each category, which is the key information to evaluate the expected energy savings if retrofitting measures are applied.

Building Energy Demand
The necessary requirements to identify the most suitable refurbishment strategies to reduce the energy consumption of a building are local climates and constructive information [52]. In the first place, it is necessary to count on the availability of a representative climate data [53] of a locality. Second, it is necessary to know the characteristics of the existing building façades [54], shading devices, air ventilation [55], fenestration types and the user occupancy [56]. The study of the potential of climatic variables helps to identify the group of passive and active techniques that, when coupled, give the best result [57]. In addition to the above, the legal restrictions and the requirements of local legislation must be considered. A dynamic simulation environment has been developed to calculate the building energy demands provided by this new digital platform [58]. This method requires an exhaustive definition of the building cases and the boundary conditions, the execution of a multi-parametric study and a pot-processing evaluation. The methodology applied is divided into three phases, as shown in Figure 3: • Development of building models. • Development of a dynamic simulation environment.

•
Post-processing evaluation and creation of building analysis modules.
Energies 2020, 13, x FOR PEER REVIEW 6 of 27 pot-processing evaluation. The methodology applied is divided into three phases, as shown in Figure  3: • Development of building models.
• Development of a dynamic simulation environment.
• Post-processing evaluation and creation of building analysis modules. The representativeness of the building stock has been selected through an input matrix that characterizes the residential building stock of the region. This information has been collected from the Spanish Technical Building Code [59], Institute for the Diversification and Saving of Energy of Spain [60], Cadastre, University of Oviedo or the Environmental Information System of the Principality of Asturias [61]. This inlet matrix is composed by the information related to the climate conditions, volume, constructive and operational characteristics.
Regarding the climate conditions, the region of Asturias is defined as temperate climate, characterized with two Köppen Geiger climatic zones: Cfb and Csb. Nevertheless, the Spanish regulation on energy savings in buildings developed a national climate classification based on seasonal severities [59]. This classification is performed as a combination of two seasonal indexes [62]: winter severity (identified by a letter) and summer severity (identified by a number). These seasonal indices are calculated as an equation that considers the influence of global solar radiation, heating degree-days and cooling degree-days [63]. Because of this climate classification, three normalized climate files are obtained for Asturias: C1, D1 and E1.
The main difference between the climate zones C1, D1 and E1 lies on the temperature and solar radiation values. The warmest zones are classified as C1 while the coldest zones are E1. Figure 4 shows the seasonal values of temperature (upper left), relative humidity (upper right), solar global radiation (lower left) and wind speed (lower right as well as the annual mean values (dark colours). Orange bars are shown for zone C1, green bars for zone D1 and blue bars for zone E1. The representativeness of the building stock has been selected through an input matrix that characterizes the residential building stock of the region. This information has been collected from the Spanish Technical Building Code [59], Institute for the Diversification and Saving of Energy of Spain [60], Cadastre, University of Oviedo or the Environmental Information System of the Principality of Asturias [61]. This inlet matrix is composed by the information related to the climate conditions, volume, constructive and operational characteristics.
Regarding the climate conditions, the region of Asturias is defined as temperate climate, characterized with two Köppen Geiger climatic zones: Cfb and Csb. Nevertheless, the Spanish regulation on energy savings in buildings developed a national climate classification based on seasonal severities [59]. This classification is performed as a combination of two seasonal indexes [62]: winter severity (identified by a letter) and summer severity (identified by a number). These seasonal indices are calculated as an equation that considers the influence of global solar radiation, heating degree-days and cooling degree-days [63]. Because of this climate classification, three normalized climate files are obtained for Asturias: C1, D1 and E1.
The main difference between the climate zones C1, D1 and E1 lies on the temperature and solar radiation values. The warmest zones are classified as C1 while the coldest zones are E1. Figure 4 shows the seasonal values of temperature (upper left), relative humidity (upper right), solar global radiation (lower left) and wind speed (lower right as well as the annual mean values (dark colours). Orange bars are shown for zone C1, green bars for zone D1 and blue bars for zone E1. Monthly mean temperature ranges from 10.5 °C (zone E1) to 13.7 °C (zone C1). The highest deviations of temperature are registered in winter, being practically non-existent during the summer. Monthly mean solar global radiation ranges from 1413 kWh/m 2 (zone E1) to 2007 kWh/m 2 and 2017 kWh/m 2 (zones D1 and C1, respectively). The highest values of solar global radiation are obtained for the climate zones C1 and D1, reaching greater deviations in spring and autumn and minimal in winter. The seasonal values of relative humidity and wind speed are quite similar for the three studied zones, with the maximum in spring and the minimum in winter. Monthly mean relative humidity ranges from 58.3% (zone D1) to 59.9% (zone E1). Finally, monthly mean wind speed ranges from 2.8 m/s (zone D1) to 2.9 m/s (zones C1 and E1).
As constructive characteristics, eight residential models that are representative of the Principality of Asturias have been selected based on the inlet dataset matrix. These models combine the type of houses with different configurations for the boundary conditions. Two types of residential houses have been modelled:   Monthly mean temperature ranges from 10.5 • C (zone E1) to 13.7 • C (zone C1). The highest deviations of temperature are registered in winter, being practically non-existent during the summer. Monthly mean solar global radiation ranges from 1413 kWh/m 2 (zone E1) to 2007 kWh/m 2 and 2017 kWh/m 2 (zones D1 and C1, respectively). The highest values of solar global radiation are obtained for the climate zones C1 and D1, reaching greater deviations in spring and autumn and minimal in winter.
The seasonal values of relative humidity and wind speed are quite similar for the three studied zones, with the maximum in spring and the minimum in winter. Monthly mean relative humidity ranges from 58.3% (zone D1) to 59.9% (zone E1). Finally, monthly mean wind speed ranges from 2.8 m/s (zone D1) to 2.9 m/s (zones C1 and E1).
As constructive characteristics, eight residential models that are representative of the Principality of Asturias have been selected based on the inlet dataset matrix. These models combine the type of houses with different configurations for the boundary conditions. Two types of residential houses have been modelled: Two configurations have been modelled to consider the boundary conditions: isolated houses (odd cases), with four façades in contact with the ambient conditions; and semi-detached houses (pair cases), with two façades in contact with the ambient conditions. Four constructive characteristics have been defined based on the Spanish regulations: before 1979, 1979-2008, 2009-2013 and after 2014. Table 3 provides the limit values of the global heat transfer coefficients for the constructive elements of the building envelope. Opaque elements (roof, walls, ground) and transparent elements (glass and frames) have been defined for each climate zone and each normative value. The glass elements are characterized by the heat transfer coefficient (U limit ) and the solar heat gain coefficient (g).
In order to complete the building models, occupancy, lighting, equipment and air renovations gains are defined using the Spanish tools for the Energy Certification of Buildings as the reference database [64].
Two periods of thermal conditioning are established, summer (June-September) and winter (October-May), each one with its respective temperature set point. Two types of annual air renovation are modelled: infiltration and ventilation. The infiltration values are defined constant, depending on the type of house and the year of construction:

Phase 2: Development of a Dynamic Simulation Environment
A dynamic simulation environment has been developed to calculate the energy demands of the building cases. This environment consists on the coupling between the dynamic simulation program TRNSYS [65] and the parameterization program GenOpt [66]. The combination of these programs automates the execution of the simulation batteries; generating successive building models for each studied option. TRNSYS has been used as the engine of the simulation environment. GenOpt identifies the studied variables of the building cases and runs a series of simulations, modifying only one of Energies 2020, 13, 5970 9 of 25 them. The outputs provided by this environment are heating, cooling and annual loads and have been used to feed the building analysis tools of the online platform.
A multi-parametric analysis has been done by modifying some studied variables of the eight representative building cases.

•
Climate zones: C1, D1 and E1. Two different analysis building modules have been developed. The first module provides the maximum retrofitting potential reached in each postal district of the Principality of Asturias. The second module provides a customized retrofitting potential for each building configuration defined by the user. Once the two databases have been created, they are represented in the GIS-based web application viewer through the Rehabilitation Potential tabs.

Solar Savings
The expected energy savings caused by a domestic solar water heating (SWH) system have been estimated by means of the f-Chart method [67], considering that each single house or flat is provided with a 2 m 2 flat solar collector tilted 25 • towards the South. The behaviour of the collector has been modelled assuming a typical performance in accordance with the following characteristic line: where η is the collector efficiency, T e is the fluid temperature at the inlet of the collector, T amb is the ambient temperature and I g is the global solar irradiance on the collector surface.
To calculate the demand in each home, a daily consumption of 120 L of domestic hot water at 60 • C has been considered.
The cold water temperature of the supply network is corrected for each postal district based on the altitude difference with the province capital, as proposed by the Spanish Technical Building Code: where T r is the water temperature in the supply network of the postal district, T rC is the water temperature in the supply network of the province capital, A is the average altitude of the postal district and A C is the altitude of the province capital.
For the climatic variables, the annual average values of irradiation and air temperatures in the postal district have been used. The annual averages of irradiation on inclined collectors are calculated from the values on a horizontal surface using the following approximation, which is justified by series of measurements at the Oviedo-State Meteorological Agency of Spain (AEMET) meteorological station: where G(0, 25 • ) is the global solar irradiation on a surface tilted 25 • towards the South, and G(0) is the value on a horizontal surface. The mass flow rate of 0.028 kg/s is assumed for the fluid circulating through the collector, while the installation is completed with an accumulator with a capacity of 150 litres and an internal heat exchanger with an efficiency of 90%.

Geothermal Savings
The geological data under Section 2.1 are used to estimate the savings derived from the geothermal resource, for a 100 m 2 rehabilitated single-family building, based on the building typologies and the following input variables: • Building construction year.
Age of the current heating equipment.
The length of the required geothermal ground-heat exchanger is calculated as a previous step to perform a study of energy savings, CO 2 emissions and economic advantages compared to the current heating system, by means of the Equation (4) [49]: where L is the length of the vertical exchanger, Q is the heating power of the building, COP is the equipment efficiency (Coefficient Of Performance), R P is the thermal resistance per unit length of the buried exchanger pipe, R S is the thermal resistance per unit length of the ground, F is the utilization factor, T L is the ground temperature and T MIN is the minimum temperature. The energy demand per surface unit of each building has been calculated keeping in mind the envelope thermal characteristics and the building typology, supported by data from the project itself. The proposed system, based on heat pump technology, is powered by geothermal energy and auxiliary electrical energy. Assuming a typical heat pump COP of 4.9, necessary electrical energy and free geothermal energy obtained from the subsoil are calculated for each building.
Geothermal systems are compared with conventional natural gas or gas oil systems in terms of CO 2 emissions. It takes the different performances into account according to the age of the current heating equipment and CO 2 emissions factors [68] (Table 4). For the final emissions calculation, these factors are multiplied by the energy consumed. The geothermal system is also compared in terms of energy costs. It considers the average cost in EUR per kWh of electricity [69], gas oil [70] and natural gas [69], which is also multiplied by the energy consumed (Table 4).
In this way, the results of energy, emissions and cost savings are achieved for different buildings under input parameters.

Results
A digital platform has been developed to collect all the results obtained in relation to the different energy resources studied as well as to the multi-parametric analysis of the energy performance of buildings. The main purpose of this platform is to provide end-users with the results of the energy savings assembled potential according to the interaction of the different technologies and available energy resources. This platform has been built upon GIS technologies and it is composed of two core modules which have been designed taking user-friendly requirements and an interactive visualization into account. The first module is the viewer, which is intended to provide information about the maximum potential of the energy savings achieved for building renovation in the Principality of Asturias. On the other hand, the second module of the platform makes it possible to obtain a rehabilitation potential for each customized building model evaluated. In this module, users can modify the input building conditions, adjusting them to their particular characteristics. The input variables that can be modified in a personalized way are: local climate zone, type of building, year of construction, number of floors, floor area, type of windows and percentage of summer shading on the north, south, east and west façade individually. The outputs values provided by this second module are the percentage of annual heating, cooling and total retrofitting potentials. Figure 5 shows the main screen of the viewer module that is seen when entering the platform. Two different spaces can be visualized: a map on the right side and a space reserved for a Layer Tree on the left one. The information is shown depicted on a map of the Principality of Asturias, which is divided into postal districts. To plot the information, an enabling/disabling system is used, which can be found in the Layer Tree. When any of these layers is enabled, the associated information is depicted according to a colour scale, affecting every single postal district separately. Additional information can be obtained when the cursor is hovered over each of them. Within the different subsets of information layers belonging to the Layer Tree, a specific group called "Potentials" can be found. The different layers that make the aforementioned group provide access to the second module of the platform, the case study simulation module. This module can be accessed in the same way as the additional information pop-ups that appear in the viewer when hovering the cursor over each postal code. Once inside, there is a series of drop-down menus that allow the selection of discrete values for the definition of the specific case to be simulated. The outputs of each simulation are displayed to the user by means of numerical results or graphical representations after the selection of the relevant data.
The technologies employed in the development of the platform are open-source libraries like Leaflet, Angular and Highcharts. For the time being, the platform is only available in Spanish. Within the different subsets of information layers belonging to the Layer Tree, a specific group called "Potentials" can be found. The different layers that make the aforementioned group provide access to the second module of the platform, the case study simulation module. This module can be accessed in the same way as the additional information pop-ups that appear in the viewer when hovering the cursor over each postal code. Once inside, there is a series of drop-down menus that allow the selection of discrete values for the definition of the specific case to be simulated. The outputs of each simulation are displayed to the user by means of numerical results or graphical representations after the selection of the relevant data.

Solar Thermal Potential
The technologies employed in the development of the platform are open-source libraries like Leaflet, Angular and Highcharts. For the time being, the platform is only available in Spanish.

Solar Thermal Potential
Potential savings in domestic SWH consumption have been computed according to Section 2.5 for every postal district and represented in Figure 6. Under the assumed conditions, solar savings over 30%, which is the minimum required by the Spanish regulation for a single family house, are observed for most of the region. The software also provides the possibility of varying the number and type of collectors, as well as considering factors of loss due to orientation, tilt and shading. Although the climatic variation within Asturias is covered with the three climatic zones, C1, D1 and E1, the new platform allows the evaluation of the general performance of solar installations, showing relevant differences in the same climatic zone. Table 5 shows the potential energy savings during a year, the annual mean efficiency of the installation for different postal districts and the yearly CO2 savings, assuming a gas boiler as auxiliary energy source.  Although the climatic variation within Asturias is covered with the three climatic zones, C1, D1 and E1, the new platform allows the evaluation of the general performance of solar installations, showing relevant differences in the same climatic zone. Table 5 shows the potential energy savings during a year, the annual mean efficiency of the installation for different postal districts and the yearly CO 2 savings, assuming a gas boiler as auxiliary energy source. It is observed that the annual efficiency is almost constant in all climatic zones studied, while the energy savings can have great variation, particularly in zones C1 and D1. In the case of zone C1, the solar energy savings vary from 29.6% in Coaña 33710 to 51.3% in Soto del Barco 33125. A similar scenario is obtained for the D1 climatic zone, where the solar savings range from 26.8% in Coaña 33716 to 57.4% in Cudillero 33156, which are also the minimum and maximum in all the analysed area. Similar results can be observed in other cases.

Geothermal Potential
The digital platform allows the calculation of the geothermal potential for each building, mainly based on the building location and the thermal conductivity of the soil. Additional input parameters for each case are: year of construction, surface area, current fuel and age of the boiler. Based on these variables, the tool determines the power extraction rate, the necessary length of the geothermal exchanger and the economic and environmental savings that the geothermal system would have in comparison to a conventional energy system.

Geothermal Exchanger Length Depending on the Thermal Conductivity and Year of Construction
It is considered that the total heating demand is covered by geothermal heat pump. Based on this, the necessary drilling length is calculated in each case. The data of thermal conductivity in different regions are essential in order to know the drilling depth needed. There is a great impact of conductivity in these drilling meters needed to meet the energy requirements. Carrying out optimal sizing of the geothermal exchanger is essential, because geothermal drilling costs represent a high percentage of the initial investment.
As an example of the high influence of thermal conductivity on the length of the geothermal exchanger, a comparison of the same type of building is made in three locations: Gijón, Oviedo and Cangas del Narcea [48]. The type of building chosen is 100 m 2  It is observed that the annual efficiency is almost constant in all climatic zones studied, while the energy savings can have great variation, particularly in zones C1 and D1. In the case of zone C1, the solar energy savings vary from 29.6% in Coaña 33710 to 51.3% in Soto del Barco 33125. A similar scenario is obtained for the D1 climatic zone, where the solar savings range from 26.8% in Coaña 33716 to 57.4% in Cudillero 33156, which are also the minimum and maximum in all the analysed area. Similar results can be observed in other cases.

Geothermal Potential
The digital platform allows the calculation of the geothermal potential for each building, mainly based on the building location and the thermal conductivity of the soil. Additional input parameters for each case are: year of construction, surface area, current fuel and age of the boiler. Based on these variables, the tool determines the power extraction rate, the necessary length of the geothermal exchanger and the economic and environmental savings that the geothermal system would have in comparison to a conventional energy system.

Geothermal Exchanger Length Depending on the Thermal Conductivity and Year of Construction
It is considered that the total heating demand is covered by geothermal heat pump. Based on this, the necessary drilling length is calculated in each case. The data of thermal conductivity in different regions are essential in order to know the drilling depth needed. There is a great impact of conductivity in these drilling meters needed to meet the energy requirements. Carrying out optimal sizing of the geothermal exchanger is essential, because geothermal drilling costs represent a high percentage of the initial investment.
As an example of the high influence of thermal conductivity on the length of the geothermal exchanger, a comparison of the same type of building is made in three locations: Gijón, Oviedo and Cangas del Narcea [48]. The type of building chosen is 100 m 2      It is observed that the annual efficiency is almost constant in all climatic zones studied, while the energy savings can have great variation, particularly in zones C1 and D1. In the case of zone C1, the solar energy savings vary from 29.6% in Coaña 33710 to 51.3% in Soto del Barco 33125. A similar scenario is obtained for the D1 climatic zone, where the solar savings range from 26.8% in Coaña 33716 to 57.4% in Cudillero 33156, which are also the minimum and maximum in all the analysed area. Similar results can be observed in other cases.

Geothermal Potential
The digital platform allows the calculation of the geothermal potential for each building, mainly based on the building location and the thermal conductivity of the soil. Additional input parameters for each case are: year of construction, surface area, current fuel and age of the boiler. Based on these variables, the tool determines the power extraction rate, the necessary length of the geothermal exchanger and the economic and environmental savings that the geothermal system would have in comparison to a conventional energy system.

Geothermal Exchanger Length Depending on the Thermal Conductivity and Year of Construction
It is considered that the total heating demand is covered by geothermal heat pump. Based on this, the necessary drilling length is calculated in each case. The data of thermal conductivity in different regions are essential in order to know the drilling depth needed. There is a great impact of conductivity in these drilling meters needed to meet the energy requirements. Carrying out optimal sizing of the geothermal exchanger is essential, because geothermal drilling costs represent a high percentage of the initial investment.
As an example of the high influence of thermal conductivity on the length of the geothermal exchanger, a comparison of the same type of building is made in three locations: Gijón, Oviedo and Cangas del Narcea [48]. The type of building chosen is 100 m 2      It is observed that the annual efficiency is almost constant in all climatic zones studied, while the energy savings can have great variation, particularly in zones C1 and D1. In the case of zone C1, the solar energy savings vary from 29.6% in Coaña 33710 to 51.3% in Soto del Barco 33125. A similar scenario is obtained for the D1 climatic zone, where the solar savings range from 26.8% in Coaña 33716 to 57.4% in Cudillero 33156, which are also the minimum and maximum in all the analysed area. Similar results can be observed in other cases.

Geothermal Potential
The digital platform allows the calculation of the geothermal potential for each building, mainly based on the building location and the thermal conductivity of the soil. Additional input parameters for each case are: year of construction, surface area, current fuel and age of the boiler. Based on these variables, the tool determines the power extraction rate, the necessary length of the geothermal exchanger and the economic and environmental savings that the geothermal system would have in comparison to a conventional energy system.

Geothermal Exchanger Length Depending on the Thermal Conductivity and Year of Construction
It is considered that the total heating demand is covered by geothermal heat pump. Based on this, the necessary drilling length is calculated in each case. The data of thermal conductivity in different regions are essential in order to know the drilling depth needed. There is a great impact of conductivity in these drilling meters needed to meet the energy requirements. Carrying out optimal sizing of the geothermal exchanger is essential, because geothermal drilling costs represent a high percentage of the initial investment.
As an example of the high influence of thermal conductivity on the length of the geothermal exchanger, a comparison of the same type of building is made in three locations: Gijón, Oviedo and Cangas del Narcea [48]. The type of building chosen is 100 m 2      Keeping all other variables equal, in the postal district 33817 in Cangas del Narcea (Figure 9), 103 drilling meters are needed while in the postal district 33815, 140 m are needed. As both postal districts are in the same location, there is an increase of 35.92% affecting the cost of geothermal drilling. The larger the floor area, the more substantial the differences will be in the same locations.
The year of construction also affects geothermal sizing, as the building demand is different. In a building with a 100 m 2 floor plant placed in the postal district 33191 of Oviedo (Figure 8), other parameters being equal, 140 drilling meters are needed in the year of construction 1800-1979 whereas in the year of construction 2014-2016 only 97 m are needed. It is clear that building regulations have made it possible to reduce the rates of thermal demand. In this case, this is 44.33% more favourable, which undoubtedly has an impact on investment in the geothermal system.

Energy and CO 2 Emission Analysis
The heat pump is sized from the input data to cover all the heating building demands. It works throughout the year and the temperature of soil remains constant regardless of the outside conditions. The energy savings produced after replacing the current conventional heating system [58][59][60] by the geothermal heat pump system are analysed for the three previous locations. Three cases with different input values have been analysed: made it possible to reduce the rates of thermal demand. In this case, this is 44.33% more favourable, which undoubtedly has an impact on investment in the geothermal system.

Energy and CO2 Emission Analysis
The heat pump is sized from the input data to cover all the heating building demands. It works throughout the year and the temperature of soil remains constant regardless of the outside conditions. The energy savings produced after replacing the current conventional heating system [58][59][60] by the geothermal heat pump system are analysed for the three previous locations. Three cases with different input values have been analysed:    made it possible to reduce the rates of thermal demand. In this case, this is 44.33% more favourable, which undoubtedly has an impact on investment in the geothermal system.

Energy and CO2 Emission Analysis
The heat pump is sized from the input data to cover all the heating building demands. It works throughout the year and the temperature of soil remains constant regardless of the outside conditions. The energy savings produced after replacing the current conventional heating system [58][59][60] by the geothermal heat pump system are analysed for the three previous locations. Three cases with different input values have been analysed:   • Case 3 ( Figure 12): 400 m 2 building in Cangas del Narcea (postal district 33815) built between 1800 and 1979, replacing a gas oil boiler that is more than 15 years old. Due to the age of the building and the larger floor area, the savings are higher than in cases 1 and 2. A total of 111,788 kWh are obtained free of charge per year and 47,990 kg of CO2 are saved annually compared to  • Case 3 ( Figure 12): 400 m 2 building in Cangas del Narcea (postal district 33815) built between 1800 and 1979, replacing a gas oil boiler that is more than 15 years old. Due to the age of the building and the larger floor area, the savings are higher than in cases 1 and 2. A total of 111,788 kWh are obtained free of charge per year and 47,990 kg of CO2 are saved annually compared to the conventional system.  The implementation of the geothermal system in cases 1 and 2 would avoid the emission of about 75% of CO 2 emissions, while in case 3 it is almost 84%. These data make the geothermal heating system a renewable system with special environmental interest, more accentuated in large and old buildings. In economic terms, cases 1 and 2 mean a saving of between EUR 1000 and 1500. In the case of buildings with old systems and high energy consumption, there are greater savings. In the case 3, there is about EUR 12,000 in annual savings compared to the previous gas oil system. Profitability analysis can be carried out, comparing drilling and heat pump systems versus the conventional system. In this case, with a 12 kW heat pump, it is estimated that the return on investment period is 6-7 years with an Internal Rate of Return (IRR) of around 17-18%. From an economic point of view, geothermal energy provides greater profitability than solar thermal energy, because the energy demand for SWH applications is relatively low in residential buildings, even meeting the requirements of regulations to reduce energy consumption and CO 2 emissions.

Retrofitting Potential
The mild climatic conditions registered in all regions of the Principality of Asturias have led to high percentages of heating demands in comparison with the cooling demands throughout the year. The annual heat contribution ranges, both for single-family houses and for blocks, between 71 and 93% for isolated configurations and between 83 and 95% for semi-detached configurations. The range of variability is due to climate conditions and the year of construction of the building. Figure 13 shows the heating load contribution to the annual total loads reached by the eight representative building cases (Cs1-Cs8) for the climate zones C1, D1 and E1 and for the four Spanish construction regulations. In this figure, the floor area considered for blocks is 400 m 2 .
The heating load contributions are higher for normative regulations before 2008 and lower for normative regulations after 2008, with different percentages depending on the climate zone.

•
For Spanish construction regulations before 2008 (red and orange lines): Zone C1: 88 to 93%. Zone D1: 90 to 94%. Zone E1: 91 to 95%. Two retrofitting potentials have been calculated superimposing the statistical study of the building stock of Asturias with the building thermal loads obtained by the simulation models. Module 1 represents the maximum retrofitting potential achieved in each studied zone of Asturias while module 2 represents the potential achieved by a specific building defined by the user.

Module 1: Maximum Retrofitting Potential
Module 1 represents the maximum retrofitting potential achieved in all the regions of the Principality of Asturias ( Figure 14). These potentials, obtained with respect to the normative requirements after 2014, depend on climate conditions and the year of construction of the building. Analysing all the regions, the highest retrofitting potentials are reached in the western area of Asturias. This is due to the high percentage of single-family houses constructed before 2008, giving rise to buildings with poor energy requirements. With regard to the building typology, higher retrofitting percentages are achieved in singlefamily houses in comparison with blocks. Isolated configurations obtained higher energy potentials Two retrofitting potentials have been calculated superimposing the statistical study of the building stock of Asturias with the building thermal loads obtained by the simulation models. Module 1 represents the maximum retrofitting potential achieved in each studied zone of Asturias while module 2 represents the potential achieved by a specific building defined by the user. Two retrofitting potentials have been calculated superimposing the statistical study of the building stock of Asturias with the building thermal loads obtained by the simulation models. Module 1 represents the maximum retrofitting potential achieved in each studied zone of Asturias while module 2 represents the potential achieved by a specific building defined by the user.

Module 1: Maximum Retrofitting Potential
Module 1 represents the maximum retrofitting potential achieved in all the regions of the Principality of Asturias ( Figure 14). These potentials, obtained with respect to the normative requirements after 2014, depend on climate conditions and the year of construction of the building. Analysing all the regions, the highest retrofitting potentials are reached in the western area of Asturias. This is due to the high percentage of single-family houses constructed before 2008, giving rise to buildings with poor energy requirements.  With regard to the building typology, higher retrofitting percentages are achieved in single-family houses in comparison with blocks. Isolated configurations obtained higher energy potentials than semi-detached houses. There is a huge number of residential houses constructed with old construction regulations (before 2008), reaching the highest percentages of retrofitting potential.
For isolated configurations, the highest energy savings are obtained for C1 climate conditions while the lowest potentials are obtained for E1 climate conditions. For semi-detached configurations and normative requirements before 2008, the three studied climate zones reach similar energy saving potentials. For semi-detached configurations and normative requirements after 2008, the highest potentials are reached for zones C1 while the lowest potentials are reached for zones E1.

Module 2: Customize Retrofitting Potential
Module 2 supplies an estimation of customized retrofitting potentials based on the inlet variables defined by users. The use of module 2 has a great potential to quantify the energy savings achieved with the refurbishment of customized residential houses in specific regions. These potentials are obtained with respect to the minimal case that corresponds to one of the representative buildings modelled with the normative requirements after 2014.
As an example of use, 48 cases have been calculated for three municipalities of the Principality of Asturias (one for each climate zone): Gijon (zone C1), Oviedo (zone D1) and Cangas del Narcea (zone E1), as can be seen in Figure 15. Total annual retrofitting potentials are obtained for each location using module 2.
Energies 2020, 13, x FOR PEER REVIEW 18 of 27 than semi-detached houses. There is a huge number of residential houses constructed with old construction regulations (before 2008), reaching the highest percentages of retrofitting potential. For isolated configurations, the highest energy savings are obtained for C1 climate conditions while the lowest potentials are obtained for E1 climate conditions. For semi-detached configurations and normative requirements before 2008, the three studied climate zones reach similar energy saving potentials. For semi-detached configurations and normative requirements after 2008, the highest potentials are reached for zones C1 while the lowest potentials are reached for zones E1.

Module 2: Customize Retrofitting Potential
Module 2 supplies an estimation of customized retrofitting potentials based on the inlet variables defined by users. The use of module 2 has a great potential to quantify the energy savings achieved with the refurbishment of customized residential houses in specific regions. These potentials are obtained with respect to the minimal case that corresponds to one of the representative buildings modelled with the normative requirements after 2014.
As an example of use, 48 cases have been calculated for three municipalities of the Principality of Asturias (one for each climate zone): Gijon (zone C1), Oviedo (zone D1) and Cangas del Narcea (zone E1), as can be seen in Figure 15. Total annual retrofitting potentials are obtained for each location using module 2. The 48 cases analyzed in the three selected cities are summarized in Table 6 and cover: • Two types of residential house: single-family and four-storey block. These typologies are the most usual in these regions.
• Two configurations for the boundary conditions: isolated and semi-detached.   The 48 cases analyzed in the three selected cities are summarized in Table 6 and cover: • Two types of residential house: single-family and four-storey block. These typologies are the most usual in these regions. The retrofitting potentials obtained for the 48 cases proposed in Gijon, Oviedo and Cangas del Narcea are shown in Figures 16-18, respectively. The x-axis represents the retrofitting potential while the y-axis represents the four Spanish regulations and the three summer shadings. For each construction regulation, there are a group of three summer shadings over the windows: 100% (S 100%) 50% (S 50%) and 0% (S 0%). To calculate the potential savings reached for each studied option, the constructive regulation defined after 2014 for buildings with summer shading factor of 100% over the external façades has been selected as a reference case. For this reason, all the cases 10 give a null value of retrofitting potential.  Figure 16. Annual retrofitting potential with respect to normative requirements after 2014 achieved in Gijon for the 48 cases proposed. Single-family houses and a four-storey block with an area of 400 m 2 are assessed for isolated configuration and semi-detached configurations.    In these three cities, the highest retrofitting potentials are achieved when there is no summer shading over the windows of the external façades, decreasing as the percentage of shading increases.
Higher retrofitting potentials are obtained for single-family houses in comparison with fourstorey blocks, especially when the construction requirements are softer (before 2008). In both types of houses, higher potentials are achieved for isolated configurations, being more remarkable for For isolated configurations, the deviation between the climate zones is very similar, being slightly higher when comparing Gijon and Oviedo. For semi-detached configurations, two tendencies have been obtained. If the construction requirements are softer (regulations before 2008), higher deviations are reached when comparing Oviedo and Cangas del Narcea. If the building energy requirements are stricter (regulations after 2008), slightly higher deviations are obtained when comparing Gijon and Oviedo.
As it can be seen in these figures, the annual retrofitting potentials reached in the three municipalities are quite similar.
In these three cities, the highest retrofitting potentials are achieved when there is no summer shading over the windows of the external façades, decreasing as the percentage of shading increases.
Higher retrofitting potentials are obtained for single-family houses in comparison with four-storey blocks, especially when the construction requirements are softer (before 2008). In both types of houses, higher potentials are achieved for isolated configurations, being more remarkable for construction regulations before 2008.

Discussion and Conclusions
One of the results obtained in the RehabilitaGeoSol project is the creation of a GIS platform developed to quantify the energy-savings obtained through the implementation of retrofit measures, solar thermal panels and geothermal technologies in residential buildings placed in the Principality of Asturias. This platform is composed of two different modules. Module 1 calculates, for each district of Asturias, the maximum energy-saving potential. Module 2 makes it possible to obtain a personalized potential based on the input characteristic of the platform.
The new developed GIS platform based on the interaction of different technologies and available resources is a step forward in the field of built environment towards smart energy communities, supporting the low carbon energy transition. The obtained results point to great energy mitigation potential, justifying the need for public authorities to develop policies aimed at reducing the energy demand of buildings, especially the heating loads. The generated maps could be helpful for retrofitting considerations for existing housing stock in the Principality of Asturias, a key action that is completed in combination with the use of solar thermal energy and geothermal resources towards sustainable development. These policies are aligned with current construction trends of reducing the energy use during the lifecycle of buildings and the promotion of the low carbon circular economy. Not only environmental but also economical reasons promote building retrofitting vs. demolition and new construction [71].
These building adaptation practices will also help address climate change, optimal use of energy resources, greenhouse emissions and secure energy supply in the coming years, by gradually implementing Net Zero Energy Buildings and Positive Energy Districts.
The RGS platform results show the high retrofitting potential for the oldest Spanish regulations (before 2006). Eight building models (Case 1-Case 8) have been developed with the dynamic simulation program TRNSYS to characterize the energy performance of residential buildings for three climatic zones (C1, D1, and E1), four Spanish construction regulations and different percentages of summer shading on the windows. The potential is higher in houses than in blocks, being more remarkable in single or isolated building configurations. In the case of blocks, the retrofitting potential is higher for four-storey plants in comparison with seven and ten-storey plans. If the surface area of blocks is higher, the energy-savings obtained are lower. The highest energy-saving potential is reached in warmer climates (zones C1), while the lowest is reached in colder zones (zones E1). Analysing the influence of external shadings on facades, this platform obtains higher impact on buildings constructed according to stricter regulations and in single-family structures vs. blocks. Regarding climate impact, the annual thermal load increases in colder climates due to the higher heating requirements. On the contrary, the impact of shading is greater in warmer climates where the cooling demands are higher.
Local authorities and other stakeholders have to consider the retrofit potential as well as the solar and geothermal resources to establish the priority areas of intervention. The expected energy savings caused by the solar domestic hot water system have been computed by means of the f-chart method assuming that each single house or flat is provided with a 2 m 2 flat solar collector tilted 25 • towards the South. Results show that almost all the territory has solar savings over 30%, which is the minimum required by the Spanish regulations for a single-family house. The new tool allows evaluation of the general performance of solar installations, showing relevant differences in the same climate zone. The annual efficiency is almost constant in all climate zones studied; the energy savings show great variation, particularly in zones C1 and D1. In the case of zone C1, the solar energy savings vary from 29.6% to 51.3%. A similar scenario is obtained for the D1 climate zone, where the solar savings range from 26.8% to 57.4%. South-West postal districts of the Principality of Asturias present high values of the annual average of daily global solar irradiation on horizontal surfaces as well as high retrofitting potential. However, the highest values of thermal conductivity are located in north-west and east postal districts. These geological data obtained are used to calculate the length of the required geothermal ground-heat exchanger. Energy demand per surface of each building has been calculated considering the thermal envelope characteristics and the building typology. The proposed system, based on heat pump technology, is powered by geothermal energy and an input of electrical energy. The geothermal system is compared with conventional natural gas or gas oil systems in terms of CO 2 emissions. A comparison of the high influence of thermal conductivity on the length of the geothermal exchanger, for the same type of building, is made in three climatic zones (C1, D1, E1).
On the one hand, the influence of thermal conductivity on the length of the exchanger can be appreciated, requiring fewer drilling meters for higher conductivities. Furthermore, it can be seen that within the same location the hydrogeology is different, requiring different lengths of exchanger for each particular case. On the other hand, the factors that most influence energy and environmental savings are the improvement of the building envelope and the replacement of the conventional heating system by the geothermal one. Greater benefits are obtained for old buildings with large surfaces, with insulation or with old heating systems based on gas oil boilers. This has a direct effect on the initial investment and subsequent amortization.
The implementation of both renewable technologies in residential buildings reduces the CO 2 emissions to the atmosphere, with the geothermal pumps being more profitable than the solar thermal panels due to the percentage of annual thermal load cover by each technology.
It should be noted that this accessible GIS platform provides both the solar and geothermal sectors with relevant information on these renewable resources in the areas of possible implementation, solving one of the major obstacles in the implementation of low enthalpy geothermal energy and solar thermal energy in urbanized areas. Therefore, the exploitation of this GIS platform by different public and private organizations is forecast. Regarding the public entities, the platform is mainly destined to be of use to environmental government areas of the Principality of Asturias, as well as to local authorities, helping them with decision-making in the matter of energy efficiency strategies and setting priorities for programs linked to rehabilitation plans and energy actions. On the other hand, the possibility of commercialization to private companies is focused mainly on the technical or engineering sector. These companies could find, in this platform, valuable information about the areas to which they could direct their main efforts to promote their activity.