Measured Performance of a Mixed-Use Commercial-Building Ground Source Heat Pump System in Sweden

: When the new student center at Stockholm University in Sweden was completed in the fall of 2013 it was thoroughly instrumented. The 6300 m 2 four-story building with o ﬃ ces, a restaurant, study lounges, and meeting rooms was designed to be energy e ﬃ cient with a planned total energy use of 25 kWh / m 2 / year. Space heating and hot water are provided by a ground source heat pump (GSHP) system consisting of ﬁve 40 kW o ﬀ -the-shelf water-to-water heat pumps connected to 20 boreholes in hard rock, drilled to a depth of 200 m. Space cooling is provided by direct cooling from the boreholes. This paper uses measured performance data from Studenthuset to calculate the actual thermal performance of the GSHP system during one of its early years of operation. Monthly system coe ﬃ cients-of-performance and coe ﬃ cients-of-performance for both heating and cooling operation are presented. In the ﬁrst months of operation, several problems were corrected, leading to improved performance. This paper provides long-term measured system performance data from a recently installed GSHP system, shows how the various system components a ﬀ ect the performance, presents an uncertainty analysis, and describes how some unanticipated consequences of the design may be ameliorated. Seasonal performance factors (SPF) are evaluated based on the SEPEMO (“SEasonal PErformance factor and MOnitoring for heat pump systems”) boundary schema. For heating (“H”), SPFs of 3.7 ± 0.2 and 2.7 ± 0.13 were obtained for boundaries H2 and H3, respectively. For cooling (“C”), a C2 SPF of 27 ± 5 was obtained. Results are compared to measured performance data from 55 GSHP systems serving commercial buildings that are reported in the literature. rate measurement) precludes calculating COP and SPF for the entire system. This study shows that the system provides space heating consistent with the design values, and that the cooling provided is about four times higher than anticipated from the design. An SPF for SEPEMO boundary H2 (heat pump + source-side circulating pump) of 3.7 ± 0.2 was achieved. An SPF for boundary H3 (heat pump + source-side circulating pump + auxiliary heating) of 2.7 ± 0.13 was achieved. The SPF H2 value is lower than most of the reported values in the literature, which have a median value of 3.6. The SPF H3 value is lower than the only reported value of SPF H3 in the literature. For cooling, the SPF C2 value of 27 ± 5 is considerably higher than any value reported in the literature. temperatures found in This paper presents an approach to addressing the performance gap without installing research-grade instrumentation. High-quality commercially-available instrumentation selected by the owners for the building energy management system proved su ﬃ cient to measure the system SPFs and COPs for several of the SEPEMO boundary conditions with reasonable accuracy, as demonstrated with the uncertainty analysis. Another contribution of this paper is the use of a bin analysis of system COPs to identify several problems with the operating strategy that, once corrected, have signiﬁcant potential to reduce the system energy consumption.


Introduction
Buildings in developed countries consume 20%-40% of the total final energy consumption [1]; in 2017, residential and commercial buildings consumed 39% [2] of all energy in the USA. Of this, about half [1] is used to provide heating, ventilation, and cooling (HVAC). Likewise, in Sweden, in 2016, buildings consumed 39% of all energy [3] and 55% of that quantity was used for space heating and domestic hot water (DHW) [4]. An extensive amount of research, development, regulation, and investment has gone towards reducing energy consumption of building HVAC systems. Nevertheless, actual building energy consumption is often significantly higher than that planned at the design stage. This difference is sometimes referred to as the "performance gap" [5][6][7][8].
The performance gap is particularly of concern for green, low energy buildings. Scofield [9] examined site energy use intensities (Energy use intensities are a measure of building energy relatively few comparable measurements that have been reported in the literature. Nevertheless, the decreasing cost of sensors and data acquisition suggest that these analyses may be more widely available in the future. Thus, one contribution of this paper is providing a comparison benchmark for future studies.
One critical issue for representing field performance is the boundary schema-where should the boundaries be drawn when determining the system performance? Early work on defining boundary schema was done for residential systems, so we first review some of the work on field measurement of residential-building GSHP systems before covering the main focus of the review, measured-performance of commercial-building GSHP systems.

Literature Review-Residential Building GSHP System Performance Measurement
For residential buildings, field measurements of heat pump system performance have been reviewed by Gleeson and Lowe [21]. The reviewed studies were undertaken in Austria, Denmark, Germany, Sweden, Switzerland, and the UK. The studies, which covered both ground-source and air-source heat pump systems, utilized sufficient instrumentation to measure system performance on a seasonal basis. Gleeson and Lowe identified two different boundary schema; both have multiple boundary definitions that may be applied. The boundary definitions differ as to which components are included-e.g., source-side circulating pump, load-side circulating pump, fans, etc.-and as to where heat delivered is measured, e.g., to the buffer tank or from the buffer tank to the heat distribution system. The authors attempted to harmonize the field trial results so that broad comparisons could be made between different field trials. Performance of over 600 heat pump installations was examined. One notable finding is an inconsistency in results that cannot readily be explained by the quality of the equipment or the building load profiles. For 216 GSHP installations, the seasonal performance factor that accounts for circulation pumps on both sides of the heat pumps and integrated backup heating ranged from 1.4 to 5.1. Gleeson and Lowe surmise that the issue is the quality of the design and installation.
A particular focus of the Gleeson and Lowe [21] paper is harmonizing the field trial results computed with different boundary schema, each of which has different boundaries. The different boundaries are important for several reasons-these include comparison to conventional systems, benchmarking of GSHP system performance, identification of best practices, determining the causes of poor performance, and giving guidance as to how system performance may be improved. One such schema was defined by the EU project SEPEMO [22], giving heating and cooling seasonal performance factors (SPFs) for a range of boundaries, as shown in Table 1.
The SEPEMO boundary scheme was applied to 44 heat pump systems in six European countries (Sweden, Germany, Greece, the Netherlands, France, and Austria) with varying heat sources (air, ground) and heat distribution methods (panel radiators, floor heating). All but five of these heat pump systems were installed in residential buildings, mostly single-family houses. Nordman [22] reports SPF H3 (heating system except the heat distribution inside the building) for the 44 heat pump systems, ranging from 1.3-7.3. Eight of the systems have SPF H3 below 2.6 and 15 systems (all GSHP systems) show SPF H3 above 4. Nordman [22] notes that heat pump system performance depends not only on the heat pump, but also on the climate and quality of installation. The SEPEMO project guidelines work well for smaller residential GSHP systems, but do not fully address all of the features that may be found in GSHP systems serving larger and more complex GSHP systems, such as commercial, institutional, and multi-family buildings.
Note that the SEPEMO system boundaries for SPF H1 and SPF H2 correspond to SPF C1 and SPF C2 , while the boundaries for SPF H3 do not correspond directly to those for SPF C3 . SPF H3 includes auxiliary heating but not distribution pumps/fans, whereas SPF C3 includes distribution pumps/fans, but not supplementary cooling units. Hence if there is no auxiliary heating in the heating system, SPF H2 = SPF H3 , while SPF C3 = SPF C4 for systems without supplementary cooling. Miara et al. [23] use a similar boundary scheme to the SEPEMO scheme, for residential heating only (Table 1). System boundaries 0-2 correspond to the SEPEMO boundaries 1-3, while the fourth system boundary includes circulation pumps between back-up heater and storage tanks, but excludes the storage tanks and distribution system beyond the tanks. The field study by Miara et al. included 56 GSHP pump systems for small residential buildings over a period of three years. The reported average seasonal performance factors for the 56 GSHP systems were 4.19 for SPF 0 , 3.93 for SPF 1 , 3.88 for SPF 2 , and 3.75 for SPF 3 . The DHW share was 18% over the three-year period. SPF 2 varied in the range of 3.75-3.9 over the three years when the measurements took place.
Another two system boundary schema suitable for residential GSHP heating systems are found in Koeningsdorff [24] (pages 276-277) and in the German guidelines VDI 4650 [25]. VDI 4650 (Table 1) defines four system boundary levels, where the first only includes the heat pump unit (corresponds to the SEPEMO boundary H1), the second only includes the ground source with a load-side circulation pump, the third corresponds to the SEPEMO boundary H3, and the fourth boundary includes both the source-side and load-side. The Koenigsdorff scheme consists of seven system boundaries, of which the first three follow the VDI 4650 scheme, and the seventh corresponds to the fourth boundary in the VDI scheme.

Literature Review-Commercial Building GSHP System Performance Measurement
Commercial buildings often have considerably more complex GSHP systems than residential buildings. Beyond energy use intensity measurements, several authors have made measurements of system coefficient of performance (COP) for a single day [26,27], several days [28], or several weeks [29,30].
Several authors [31][32][33][34][35] have presented detailed performance measurements of building heating systems that incorporate GSHP and other renewable sources such as solar-thermal and solar-electric on site. The first paper focuses on other aspects of the system performance. Wemhoener et al. [32,33] report monthly and annual performance factors, with an annual value of 3.92 for combined space heating and domestic hot water heating, accounting for the heat pump and source-side circulating pump. Lazzarin et al. [34,35] describes a primary energy ratio metric that is the ratio of heating provided to primary energy consumed. Use of this metric helps get around the problem of how to treat electricity generated on site. The literature review in this paper focuses on GSHP systems that are characterized with SPF or COPs with the denominator being the electricity consumed without regard to whether it is generated on site or off site.
An aquifer thermal energy storage (ATES) system is a type of multi-well open-loop GSHP system, with hot and cold reservoirs maintained by seasonally switching the extraction and reinjection wells. Schmidt and Müller-Steinhagen [36] described a rather complicated system used to provide heating to a 7000 m 2 apartment building. This system combined solar collectors, a gas-fired boiler, ATES, and heat pumps. Seasonal heating COP for the heat pumps only was reported for three years, with an average value of 4.3.
Urchueguía et al. [37] compared an air-source heat pump system and a ground-source heat pump system serving portions of a university building in Valencia, Spain. Daily and seasonal performance factors corresponding to the SEPEMO H2 and C2 boundaries were reported. For heating, SPFs of 3.5 ± 0.6 and 2.0 ± 0.3 were measured for the GSHP and ASHP systems, respectively. For cooling, SPFs of 4.3 ± 0.6 and 2.7 ± 0.4 were measured for the GSHP and ASHP systems, respectively. Montagud et al. [38] present additional data for this system after five years of operation.
Hughes [39][40][41]  Zhai and Yang [44] describe a large (280 boreholes) GSHP system serving a government building housing archives, and hence requiring close control of space air temperature and humidity. Average COPs of the heat pumps (presumably SEPEMO H1 and C1, but not explicitly stated) are given as 4.7 in the summer and 4.6 in the winter.
Vanhoudt et al. [45] describe three years of performance monitoring results for an ATES system. The heat pumps have seasonal performance factors (SEPEMO boundaries H1 and C1) of 5.6 and 5.0. Seasonal system performance factors are also reported, though the boundaries are unclear-it is not clear whether "different circulation pumps" include both source-side and load-side pump energy consumption. Average values of 5.9 and 26.1 are reported. These are higher than the heat pump COPs because some heat is provided directly from the groundwater (presumably for tempering ventilation air) and much of the cooling is provided directly from the groundwater.
Michopoulos et al. [46] monitored a ground-source heat pump system serving an office building in Greece. Weekly and seasonal performance factors corresponding to the SEPEMO H1 and C1 boundary conditions were presented. The seasonal performance factors vary from year-to-year, but are about 5-5.5 for heating and 4-4.5 for cooling.
Winiger et al. [42] report on analysis of ten ground-source heat pump systems serving non-residential buildings in Germany. These systems were monitored as part of the broader Energy-Optimized Building (EnOB) project [47]. The authors describe four system boundaries (see Table 1) that have been evaluated, though none of the boundaries include only the heat pump. One case makes use of waste heating, allowing heating without heat pump operation-in this case, the heating SPF I was 61.7. Cooling SPF results are divided between systems using direct cooling (i.e., fluid circulated directly from the ground heat exchanger to the fan coil unit) and active cooling (i.e., using the heat pumps). When heating and cooling were being provided simultaneously, the electricity usage was allocated proportional to the thermal energy provided in each mode. Heating SPF II values (i.e., heat pump plus source-side circulating pump) ranged from 2.9 to 5.7. Heating and cooling SPF IV values (i.e., the entire system including the distribution system) ranged from 1.0-2.7 and 2.2-15.0, respectively, with higher cooling SPFs found for the direct cooling cases. SPF I and SPF II are given for six individual systems, with no other details provided.
Mermoud et al. [48] used the SEPEMO boundary conditions to analyze a GSHP system serving a 2200 m 2 office building in Geneva, Switzerland. The system provides heating with heat pumps and direct cooling using the ground heat exchanger fluid. Heat pump COPs are presented for five-minute interval data, showing a significant decrease from the manufacturer's data. The system utilizes two heat pumps and the authors computed separate SPF H1 and SPF H2 , for each heat pump, with values of SPF H1 of 3.5 and 3.1 and SPF H2 of 3.1 and 3.0. Overall values of SPF H3 and SPF H4 are 3.1 and 3.0.
As the system has favorable return fluid temperatures from the ground (exceeding 10 • C all year) and the system only provides space heating, the authors find the SPF values disappointing. The poor performance is attributed to the hydraulic design, which requires higher supply water temperatures than what is actually needed.
Monitoring studies of two distributed heat pump systems at the ASHRAE Headquarters building in Atlanta were presented by Southard et al. [49,50] and Spitler et al. [51]. The two systems were an air-source variable-refrigerant flow heat pump system and a ground-source heat pump system. Seasonal heating and cooling system coefficients of performance were calculated for the two systems. Because the systems were distributed and delivered heating and cooling with fans integrated in the heat pump units, and the electrical energy was measured for all heat pumps and the source-side circulating pump together, the coefficients of performance correspond to SEPEMO levels H4 and C4. The air-source heat pump system had seasonal heating and cooling COPs of 2.0 ± 0.1 and 2.5 ± 0.1, respectively. The ground-source heat pump system had seasonal heating and cooling COPs of 3.3 ± 0.2 and 4.3 ± 0.6, respectively.
Spitler et al. [51] presented system heating and cooling COPs as a function of outdoor temperature and heat pump entering fluid temperature. For both the air-source and ground-source heat pump systems, cooling COP increased with increasing outdoor air temperature and heat pump entering fluid temperature over a substantial part of the operation range. The counterintuitive behavior of the GSHP system was explained by parasitic losses due to control boards that consume power all the time, fans that run even when all compressors are off, source-side circulating pump energy consumption, and cycling losses of the heat pumps. Recent research [52,53] into cycling losses shows the decreased COP that occurs during the first minute of operation can have a variable but significant effect on heat pump performance.
Mendrinos and Karytsas [43] report on annual heating and cooling SPF for eight GSHP systems in southern Europe monitored as part of the EU Ground-Med project. Additional details can be found in Carvalho et al. [54] and Pardo and Michal [55]. Four levels of system boundaries are defined. As shown in Table 1, these are similar to the SEPEMO boundary conditions for heating, except they do not account for auxiliary heating, and the energy used by the fans in the fan-coil units is treated separately [55] from the load-side circulating pumps. Heating and cooling SPFs for all buildings for all four levels are reported graphically; heating SPF II ranges from 3.6-5.9. Cooling SPF II results are divided between systems using direct cooling (as defined above)-values between 9.8 and 16.6 are achieved, and active cooling-values between 4.9 and 6.8 were achieved.
Schibuola and Scarpa [56] made field measurements of system performance for a ground-source heat pump system providing heating to an 1840 m 2 university building originally constructed in the 16th century. Annual heating COPs corresponding to SEPEMO H1 and H2 were measured at 3.96 and 3.65, respectively. Annual cooling energy efficiency ratios (EERs) corresponding to SEPEMO C1 and C2, were measured at 4.32 and 4.02, respectively. These results were compared to simulated results for an air-source heat pump system.
Pater and Ciesielczyk [57] compare theoretical performance coefficients provided by the heat pump manufacturer with SPF obtained in real operating conditions for a GSHP installed in a 460 m 2 mixed-use building near Krakow in Poland. Measurements were carried out over three heating seasons (September-May). The reported SPF3, according to the Miara et al. [23] system boundary scheme, ranged between 3.4-3.8 over the three heating seasons, which fell within the range of expected SPF based on the manufacturer's technical data for the heat pumps.
Garber-Slaght and Peterson [58] describe a ground-source heat pump system utilizing horizontal ground heat exchangers in a sub-Arctic climate with permafrost present. The permafrost layer at the site was 3 to 7.3 m deep in 2006 and the slinky heat exchangers are installed at a depth of 2.7 m. The GSHP serves a 446 m 2 office space. Monthly heating COPs were measured for three years; an approximate average value of 3.5 may be inferred, with some decrease from the first year to the third year.
Liu et al. [59] summarized annual performance measurements of 10 ground-source heat pump systems serving commercial buildings. The ground-sources included both open loop and closed loop borehole systems, as well as a mine water system and a system utilizing municipal wastewater as the source. Most of the systems used distributed water-to-air heat pumps to provide both heating and cooling. Annual heating and cooling system coefficients of performance (SCOP) were determined for eight of the ten buildings. For the other two, difficulties in differentiating between energy used to provide heating and energy used to provide cooling led the authors to define an effective overall SCOP that combined heating and cooling. Annual heating SCOP varied between about 2.5 and 4.3; annual cooling SCOP varied between about 3 and 5.2. The most commonly identified problem was excessive pumping power caused by oversized circulation pumps, suboptimal controls, and excess hydraulic resistance. For the water-to-air heat pump systems, it was not possible to measure the heating and cooling provided, so a model-based approach was utilized. The resulting uncertainty was not characterized in the studies.
Naicker and Rees [60] present seasonal system performance factors for a GSHP system serving a university building in Leicester, England. Seasonal performance factors are presented for SEPEMO levels C1 and H1 (3.19 and 4.06, respectively). Combined cooling and heating SPF are also defined corresponding to SEPEMO levels H1, H2, and H4 boundary conditions and are referred to as SPF 1 ,  54. SPF2 and SPF4 are affected by the pumping controls that start the circulating pumps three minutes before the compressors. Under low load conditions, with short cycles, the pumping energy consumption can be as high as 30% of the heat pump energy consumption. Several approaches to improving the system performance are identified, including incorporating buffer tanks, use of a smaller capacity "lead heat pump", and use of variable speed compressors.
Gehlin et al. [61] provide a preliminary analysis of the work described in this paper-analysis of a GSHP system serving a university building in Stockholm, Sweden.

Literature Review-Summary
Long-term measured (>1 year) SPF and COP values reported in the literature for 55 GSHP systems used in commercial or multifamily residential buildings are summarized in Table 2, based on the SEPEMO definitions. In addition, there are three systems for which a combined cooling and heating SPF was given; these appear in the column "HC4". In cases where values were given for multiple years, they have been averaged. The 55 systems shown in Table 2 are sparsely located on three continents. There is a notable inconsistency in which boundaries are used, making comparisons difficult. In a number of cases, the boundaries used are not clearly defined, again making comparisons difficult. Only about half the cases give the overall (H4 or C4) system SPF. Only three of the studies report any uncertainty analysis; only Hughes [41] reports the methodology in detail. Yet this is critical to understanding the significance of the results, so a second contribution in this paper is a detailed description of the uncertainty analysis methodology.
Median values of annual SPF for heating are 4.1, 3.6, 3.1, and 2.9 for boundaries H1, H2, H3, and H4, respectively. Median values of SPF for cooling are 5.5, 6.4, and 4.2 for boundaries C1, C2, and C4, respectively. As noted by other authors [21,22] covering residential GSHP systems, there is a significant range in SPF values that cannot be explained solely by different equipment or different climatic conditions. Almost all of the past studies provide only seasonal or annual performance factors. But daily performance factors [59,60] or binned average values [51] have been useful in identifying causes for poor performance, such as heat pump cycling losses, pump controls, and equipment failures.
It is the objective of this paper to provide long-term measured system performance data including a detailed uncertainty analysis from a recently installed centralized GSHP system for a new-built office building in the Swedish capital Stockholm. The analysis includes seasonal performance factors and monthly, daily, and binned average values of coefficients of performance. These analyses show how the various system components and controls affect the performance and describe how some unanticipated consequences of the design may be ameliorated. (In the analysis of the Studenthuset GSHP system, we use the "Seasonal Performance Factor (SPF)" to refer to quantities calculated over a year, i.e., 365 sequential days, and "Coefficient of Performance (COP)" for performance factors calculated over shorter periods of time, such as monthly performance).   1 Centralized (C) or distributed (D). 2 Nominal heat pump capacities for either heating or cooling. In a number of cases, the heat pump capacity is not given and is estimated from the floor area and the specific heating or cooling load. 3 Undisturbed ground temperatures (UGT) are taken from the nearest location found in a world-wide database [62]. These values do not include local disturbances due to urbanization.

System Description
Stockholm University is located within the large campus area Frescati in central Stockholm. In 2011 the government-owned real-estate owner Akademiska Hus, responsible for management and development of buildings for research and higher education in Sweden, built a new student center at the campus ( Figure 1). The new building was completed in the fall of 2013. Ambitions were high for the student center building to be a model building regarding sustainability, ecologic footprint, and energy use. The 6300 m 2 four-story building contains office area for 130 staff, small meeting rooms, 200 study-booths for students, and a café, and was designed to be highly energy efficient with a planned total use of 25 kWh/m 2 /year (<160 MWh/year) bought energy (electricity).

System Description
Stockholm University is located within the large campus area Frescati in central Stockholm. In 2011 the government-owned real-estate owner Akademiska Hus, responsible for management and development of buildings for research and higher education in Sweden, built a new student center at the campus ( Figure 1). The new building was completed in the fall of 2013. Ambitions were high for the student center building to be a model building regarding sustainability, ecologic footprint, and energy use. The 6300 m 2 four-story building contains office area for 130 staff, small meeting rooms, 200 study-booths for students, and a café, and was designed to be highly energy efficient with a planned total use of 25 kWh/m 2 /year (<160 MWh/year) bought energy (electricity). The building services are thoroughly instrumented. Space heating and domestic hot water (DHW) are provided by a ground source heat pump (GSHP) system consisting of five 40 kW off-theshelf water-to-water heat pumps connected to a borehole field. The DHW temperature is further increased by a Legionella protection system consisting of an electric resistance heater that raises the hot water temperature from 55 °C to 60 °C and a pump that continuously recirculates hot water around the building. As will be seen later, this has a significant effect on the energy usage of the system. Space cooling is provided by direct cooling from the boreholes; with the maximum fluid temperature leaving the boreholes not to exceed 16 °C. Building occupancy during the warmest summer months (July-August) is low, while the students are on summer vacation. Measured temperatures from the boreholes do not exceed 13.4 °C during the measurement period, and dehumidification loads are not high. The bore field consists of 20 groundwater-filled boreholes in hard rock, drilled to a depth of 200 m, and fitted with single u-tubes filled with an ethanol/water mixture. The bore field is located below a landscaped courtyard, and the boreholes are drilled at an angle so that they reach under the surrounding building.
Heat distribution inside the building is provided by radiators with a larger-than-usual surface area so that the distribution temperature is 40 °C instead of 55 °C, which is more common in Sweden. Cooling distribution is done by the ventilation system and chilled beams. The system also includes heat recovery from the kitchen cooling circuit. A schematic of the Studenthuset heating and cooling system is seen in Figure 2.
Prior to building operation, total annual energy loads were anticipated to be 200 MWh heating and 34 MWh cooling. The installed capacity is 200 kW heating and 120 kW cooling from the boreholes. No auxiliary heating or cooling is installed, except for the electric resistance heater that boosts the hot water temperature to protect against Legionella. The anticipated electricity use for the building installations was approximately 14 kWh/m 2 /year. Distribution losses in the heating system were estimated to be 2 kWh/m 2 /year. Heat pump SPF was expected to be 4.5 according to the design documents, though the SPF boundaries were not defined. The building services are thoroughly instrumented. Space heating and domestic hot water (DHW) are provided by a ground source heat pump (GSHP) system consisting of five 40 kW off-the-shelf water-to-water heat pumps connected to a borehole field. The DHW temperature is further increased by a Legionella protection system consisting of an electric resistance heater that raises the hot water temperature from 55 • C to 60 • C and a pump that continuously recirculates hot water around the building. As will be seen later, this has a significant effect on the energy usage of the system. Space cooling is provided by direct cooling from the boreholes; with the maximum fluid temperature leaving the boreholes not to exceed 16 • C. Building occupancy during the warmest summer months (July-August) is low, while the students are on summer vacation. Measured temperatures from the boreholes do not exceed 13.4 • C during the measurement period, and dehumidification loads are not high. The bore field consists of 20 groundwater-filled boreholes in hard rock, drilled to a depth of 200 m, and fitted with single u-tubes filled with an ethanol/water mixture. The bore field is located below a landscaped courtyard, and the boreholes are drilled at an angle so that they reach under the surrounding building.
Heat distribution inside the building is provided by radiators with a larger-than-usual surface area so that the distribution temperature is 40 • C instead of 55 • C, which is more common in Sweden. Cooling distribution is done by the ventilation system and chilled beams. The system also includes heat recovery from the kitchen cooling circuit. A schematic of the Studenthuset heating and cooling system is seen in Figure 2.
Prior to building operation, total annual energy loads were anticipated to be 200 MWh heating and 34 MWh cooling. The installed capacity is 200 kW heating and 120 kW cooling from the boreholes. No auxiliary heating or cooling is installed, except for the electric resistance heater that boosts the hot water temperature to protect against Legionella. The anticipated electricity use for the building installations was approximately 14 kWh/m 2 /year. Distribution losses in the heating system were Energies 2019, 12, 2020 13 of 34 estimated to be 2 kWh/m 2 /year. Heat pump SPF was expected to be 4.5 according to the design documents, though the SPF boundaries were not defined.
The climate in Stockholm is characterized by a humid continental climate, with average temperatures of around −3 • C in the winter and +20 • C in the summer. The climate in Stockholm is characterized by a humid continental climate, with average temperatures of around −3 °C in the winter and +20 °C in the summer.

Measurements and Uncertainty Analysis
Measurement of the system performance relies on instrumentation installed (mostly) at the same time as the system installation for the owner's benefit in monitoring the system's performance. Some additional instrumentation was installed later to improve the capability to measure the cooling provided and/or estimate the source-side circulating pump energy consumption. Because the instrumentation system was not specifically designed to support measurement of SPFs and COPs with the SEPEMO boundaries, some necessary approximations are described in this section.

Instrumentation and Analysis
As discussed above, the system performance can be measured with different boundary conditions. The system is extensively instrumented compared to many buildings, but, unfortunately, there are no measurements made of airflow rates, so it is not possible to estimate either the heating or cooling provided to the ventilation air or the building by the heat recovery systems. This limits the possibilities for calculating SPFs and COPs with the H4, C3, and C4 boundaries. Nor are there independent measurements of compressor energy consumption, so SPFs and COPs cannot be calculated for the H1 and C1 boundaries.

Measurements and Uncertainty Analysis
Measurement of the system performance relies on instrumentation installed (mostly) at the same time as the system installation for the owner's benefit in monitoring the system's performance. Some additional instrumentation was installed later to improve the capability to measure the cooling provided and/or estimate the source-side circulating pump energy consumption. Because the instrumentation system was not specifically designed to support measurement of SPFs and COPs with the SEPEMO boundaries, some necessary approximations are described in this section.

Instrumentation and Analysis
As discussed above, the system performance can be measured with different boundary conditions. The system is extensively instrumented compared to many buildings, but, unfortunately, there are no measurements made of airflow rates, so it is not possible to estimate either the heating or cooling provided to the ventilation air or the building by the heat recovery systems. This limits the possibilities for calculating SPFs and COPs with the H4, C3, and C4 boundaries. Nor are there independent measurements of compressor energy consumption, so SPFs and COPs cannot be calculated for the H1 and C1 boundaries.

. Instrumentation and Analysis-Heat Transfer Rates
With the existing measurements, it is possible to estimate SPFs and COPs for the C2, H2, and H3 boundaries, and that is the focus of this paper. To measure the system performance with these boundaries, it is necessary to measure the heat transfer rates: 1.
from the heat pumps to the building heat distribution system, 2.
from the heat pumps to the DHW, 3.
from the auxiliary heating system to the DHW 4.
from the auxiliary heating system to the building, 5.
for cooling, from the building cooling distribution system to the ground loop (some of the heat is dissipated in the boreholes and some is used as a heat source for the heat pumps). This is synonymous with the building cooling provided.
Heat transfer from the heat pumps to the building heat distribution system is measured with an energy meter that determines the heat transfer rate calorimetrically. The energy meter measures the volume flow rate to an accuracy of ±5% or better. It measures the temperature difference with a matched pair of platinum RTDs (Resistance Temperature Detector), with maximum discrepancy of 0.05 K and median discrepancy of 0.02 K. The meter reports cumulative energy consumption, but knowledge of the flow meter and temperature sensor accuracies is needed to estimate the uncertainty in the energy consumption reported by the meter.
For purposes of determining the heat provided by the heat pumps to the DHW, only the volume flow rate is measured with a flow meter with ±5% or better accuracy. The temperature provided by the heat pumps is controlled to 55 • C. The incoming temperature from the Stockholm water supply is not measured; therefore, it was estimated based on measurements made at seven different buildings in Stockholm [63,64]. The incoming temperature varies over the year and from location to location. A heuristic equation to describe the incoming water temperature (in • C) was developed: where n = day of the year; 365 ≥ n ≥ 1; n min = day of the year with minimum temperature (day 60) The uncertainty for this approximation is estimated as ±2.5 • C. The average heating energy (kWh) provided by the heat pumps to the domestic hot water for any time interval is then given by: where V = integrated volume flow over the time interval, m 3 ; ρ = density of the water, kg/m 3 ; c p = specific heat of the water, kJ/(kg·K); T w,OHP = temperature of heated water at outlet from heat pumps, • C. Heating provided by the Legionella protection system is estimated based on two contributions: (1) a steady consumption of 3 kW of electricity due to the recirculation pump and dissipated energy from the hot water that was continuously circulated around the building, and (2) energy provided by the electric resistance heater to heat the DHW from 55 • C to 60 • C: where V = integrated volume flow over the time interval, m 3 ; T w,ORH = temperature of heated water at outlet from resistance heater, • C. The electric meter measuring the heat pumps, hot water recirculation pump, and the electric resistance heater measures 3 kW of electricity for many hours during the year when there is no domestic hot water use and the heat pumps are not running. This 3 kW corresponds to the energy consumed by the pump plus the energy dissipated by the recirculation of hot water and which is made up by the resistance heater operating. For these hours, the heat is dissipated in the building. This heat may be useful or it may form part of the cooling load. As a rough attempt to distinguish between the two possibilities, we deemed the heat dissipated useful if the building had a heating load.
The cooling provided by the system is measured with the same type of instrumentation as used for measuring the building heating provided.

Instrumentation and Analysis-Electrical Energy
In addition to estimating the various heat transfer rates, it is also necessary to measure the electrical energy used by the heat pumps, source-side circulating pump(s), and the auxiliary heating system. The system electrical energy consumption is monitored continuously and recorded on an hourly basis.
One meter logs electricity use for the five heat pumps (compressors and built-in circulation pumps), including the heat pump that is dedicated to DHW heating. It also includes electricity consumed by the Legionella protection system. This meter has an accuracy of ±1%. In order to estimate the electrical energy consumed by the heat pumps for boundary H2, it is necessary to be able to subtract out the energy consumed by the Legionella protection system.
As discussed above, electrical power for the recirculation pump is essentially constant over the year at a level of 3 kW. Electrical energy for the electric resistance heater can be estimated calorimetrically using Equation (3). These values are subtracted from the meter measurements to obtain the electrical energy consumed by the heat pumps.
Electrical energy consumed by the source-side circulation pump was metered along with the electrical energy for fans used for air conditioning, circulation pumps on the load-side (distribution), and circulation pumps on the source-side (boreholes), as well as electricity used for running the rotary exhaust air heat exchangers in the kitchen and building. However, for purposes of our analysis, it was necessary to be able to separately estimate the electrical energy consumed by the source-side circulation pump. This variable speed pump is controlled to maintain a minimum flow rate of 8 L/s. A separate set of measurements over a two-week period was made to allow estimation of the electricity (kW) used by the source-side circulation pump as a function of flow rate with this equation: where . V = volume flow rate averaged over the time interval, L/s. The resulting power given by Equation (4) is estimated to have an uncertainty of ±0.19 kW. A further complication arises for the source-side circulation pump-for many, many hours of operation, both heating and cooling are being provided by the system. In this case, some of the electrical energy should be allocated to the heating electrical energy and some should be allocated to the cooling electrical energy. How should this allocation be done? We chose to allocate it based on the amount of heating and cooling provided, without regard to the fact that the heating is provided via heat pumps, panel radiators, and fan coil units, while the cooling is provided "directly" with fan coil units and chilled beams. The heating provided does not include incidental heating provided by the Legionella protection system. Specifically, the fraction of electrical energy consumed by the circulating pump that is allocated to the heating electrical energy is determined as: where q BldgHtg = average heating provided to the building over the time interval, kW; q DHW = average heating provided to the domestic hot water over the time interval, kW; q BldgClg = average cooling provided to the building over the time interval, kW. The fraction of the electrical energy consumed by the circulating pump allocated to cooling is then: The scheme used to allocate the electrical energy consumed by the circulating pump is admittedly arbitrary, but a better scheme has not, to the authors' knowledge, been suggested in the literature. A similar approach was used by Winiger et al. [42].

Uncertainty Analysis
It is highly desirable to estimate the uncertainty of the quantities of interest, i.e., COPs, for several reasons. These include understanding the significance of the results-e.g., are the results meaningful; are trends in the results meaningful? The uncertainty analysis can also shed light on the design and specification of instrumentation for future monitoring projects.
The uncertainty analysis, which involves calculation of the propagation of uncertainties from the physical measurements to the final quantities of interest, follows the general procedures described by Taylor [65]. Like Taylor, we use the words "error" and "uncertainty" interchangeably.
To set the context, the final quantities of interest are primarily coefficients of performance over some time period-annual, monthly-or over all times in a specific outdoor air temperature bin or entering fluid temperature bin. So, while random errors exist, they are usually mitigated by a large number of measurements over which a specific COP value is computed. One exception is COPs computed for bins with very few hours, e.g., less than 10 measured hours. Accordingly, systematic errors are of more concern and random errors are neglected in this uncertainty analysis.

Uncertainty Analysis-Heating and Cooling Provided
The building heating and cooling are both measured with energy meters that incorporate a flow meter with ±5% accuracy and a matched pair of temperature sensors that are selected to have a discrepancy (i.e., measurement error relative to one another) of no more than 0.05 K. Therefore, the individual uncertainties can be expressed: (We adopt the convention here that fractional uncertainties are expressed with a lower-case "e" and absolute uncertainties given in the units of the quantity in question are expressed with an upper-case "E". If the quantity in question is "x", E x = x·e x .) Since errors in the flow rate and temperature difference measurement are independent of each other, they can be added in quadrature to determine the uncertainty of the heat transfer rate at each hour.
Since the errors in determining the density and specific heat are relatively small compared to the volume flow rate uncertainty, these errors are neglected. Over the one-year period, the uncertainty in the building heating provided has an average value of ±6% and a maximum value of ±11%. Over the same period, the uncertainty in the building cooling provided has an average value ±5% and a maximum value of ±26%. The maximum values occur for hours with very low temperature differences and hence low heat transfer rates.
These hourly uncertainties are aggregated for time periods or bins by summing the absolute uncertainties, e.g., to find the uncertainty in building heating provided for a one-month period, with 720 h: The uncertainty in the heating provided by the heat pumps to the domestic hot water is calculated similarly to that of the building heating and cooling, but the uncertainty in the temperature difference is somewhat higher. The uncertainty for the inlet temperature (Equation (1)) is taken as ±2.5 • C, and the uncertainty for the outlet temperature is taken as ±1 • C. The two errors are independent, so the uncertainty in the temperature difference may be taken as: Equation (10) is used to estimate the final uncertainty in the heating provided by the heat pumps to the domestic hot water. Although the uncertainty in the temperature difference is higher than that for the building heating and cooling, the temperature difference is significantly higher, so the average uncertainty is ±8% and the maximum (when rounded to the nearest %) is also ±8%.
Calculation of uncertainty for the DHW heating provided by the Legionella protection system is similar to that for the DHW heating provided by the heat pumps. The inlet temperature (55 • C) and outlet temperature (60 • C) are both assumed to have uncertainties of ±1 • C. The flow meter uncertainty is ±5%.
Adding the flow uncertainty in quadrature, the overall uncertainty for the DHW heating provided by the Legionella protection system is ±28%.
Calculation of uncertainty for the useful building heating provided (unintentionally) by the Legionella protection system in the form of heat leakage from piping was estimated as ±3.3%, based on the observed variation in the minimum electricity consumption. This does not include uncertainty due to the determination of whether or not the heat dissipated is useful. Presumably, the uncertainty regarding whether or not the heat dissipated is useful is higher during the shoulder seasons and lower during the peak heating and cooling seasons, but it seems unlikely to have a significant effect on the computed coefficients of performance for boundary H3.

Uncertainty Analysis-Electrical Energy Consumed
Although the electric meters are accurate to ±1% and therefore add very little uncertainty to the coefficients of performance, our analysis requires several additional approximations that add additional uncertainty.
The additional uncertainty in the electrical energy associated with boundaries H2 and C2 is due to the electric meter measuring the heat pumps (including the circulating pumps that are inside the heat pump cabinets) and the Legionella protection system electrical energy, but not measuring the source-side circulating pump. Furthermore, the source-side circulating pump electrical energy must be allocated to heating and cooling. Then, the uncertainty in the electrical power and energy consumed by the heat pumps and source-side circulating pump under heating conditions for use in calculating COP H2 is given by: where e P = fractional uncertainty in the meter measurement, ±1%; P = electrical energy measured by the meter, kW; E SSCP = absolute error in the estimate of the source-side circulating pump electrical energy consumption, ±0.19 kWh/h. For boundary condition C2, the only electricity consumption is that of the source-side circulating pump, so: For boundary condition H3, the meter measurement includes the Legionella protection system, so the uncertainty is given by:

Uncertainty Analysis-Coefficients of Performance
With the uncertainties in heating/cooling provided and uncertainties in electrical energy determined, the uncertainties of the coefficients of performance can readily be determined. As shown by Taylor [65] if the uncertainties in the numerator and denominator are independent, the uncertainty of the quotient is given by: where e qt , e n , e d are the fractional uncertainties of the quotient, numerator, and denominator respectively. This formulation was used to calculate the uncertainties for each COP boundary.

Results
This paper uses measured performance data from Studenthuset to characterize the actual thermal performance of the GSHP system during its third year of operation. Measured data from the period 1 April 2016 through 31 March 2017 are used. Monthly energy loads for this period are shown in Figure 3. This period was chosen because some additional instrumentation was added at the beginning of the period, allowing more accurate measurement of the cooling provided. The measured data, processed data, and Visual Basic for Applications computer code are available in the Supplementary Materials. For boundary condition H3, the meter measurement includes the Legionella protection system, so the uncertainty is given by:

Uncertainty Analysis-Coefficients of Performance
With the uncertainties in heating/cooling provided and uncertainties in electrical energy determined, the uncertainties of the coefficients of performance can readily be determined. As shown by Taylor [65] if the uncertainties in the numerator and denominator are independent, the uncertainty of the quotient is given by: where eqt, en, ed are the fractional uncertainties of the quotient, numerator, and denominator respectively. This formulation was used to calculate the uncertainties for each COP boundary.

Results
This paper uses measured performance data from Studenthuset to characterize the actual thermal performance of the GSHP system during its third year of operation. Measured data from the period 1 April 2016 through 31 March 2017 are used. Monthly energy loads for this period are shown in Figure 3. This period was chosen because some additional instrumentation was added at the beginning of the period, allowing more accurate measurement of the cooling provided. The measured data, processed data, and Visual Basic for Applications computer code are available in the Supplementary Materials.

Load Characteristics
The total energy loads for the selected period were 194 MWh space heating, 130 MWh space cooling, and 33 MWh domestic hot water (DHW) including the Legionella protection system (LPS).
As can be inferred from Figure 3, heating and cooling is provided simultaneously over a substantial part of the year. As discussed in Section 3.1.2, a complicating problem in determining the performance factors of the system is that there is no standard way to allocate electrical energy usage to cooling and heating respectively, when the system is doing both simultaneously. Equation (5

Load Characteristics
The total energy loads for the selected period were 194 MWh space heating, 130 MWh space cooling, and 33 MWh domestic hot water (DHW) including the Legionella protection system (LPS).
As can be inferred from Figure 3, heating and cooling is provided simultaneously over a substantial part of the year. As discussed in Section 3.1.2, a complicating problem in determining the performance factors of the system is that there is no standard way to allocate electrical energy usage to cooling and heating respectively, when the system is doing both simultaneously. Equation (5) defined the "heating fraction" used to allocate the electrical energy consumption of the source-side circulation pump when the Studenthuset system is simultaneously providing building heating and DHW from the heat pumps while also providing cooling somewhere in the building. Figure 4 shows the heating fraction as a function of outdoor air temperature. Each point on the plot represents the average heating fractions for a particular outdoor air temperature bin of width 1 • C.
Energies 2019, 12, x FOR PEER REVIEW 20 of 35 defined the "heating fraction" used to allocate the electrical energy consumption of the source-side circulation pump when the Studenthuset system is simultaneously providing building heating and DHW from the heat pumps while also providing cooling somewhere in the building. Figure 4 shows the heating fraction as a function of outdoor air temperature. Each point on the plot represents the average heating fractions for a particular outdoor air temperature bin of width 1 °C.
The heating system has a maximum capacity of 200 kW, but as seen in the duration curve in Figure 5, heating is mostly provided at 20-60 kW. Maximum capacity of the borehole free-cooling system is 120 kW.

System Boundaries
In this paper we have used the system boundaries for calculating system performance for heating and cooling, as suggested by the EU project SEPEMO [22]. The SEPEMO project was mainly focused on single-family houses with heat pumps for heating or cooling, while Studenthuset is an office building with both heating and cooling, and thus has a higher degree of system complexity. Figure 6 shows the SEPEMO system boundaries applied to Studenthuset.  The heating system has a maximum capacity of 200 kW, but as seen in the duration curve in Figure 5, heating is mostly provided at 20-60 kW. Maximum capacity of the borehole free-cooling system is 120 kW. defined the "heating fraction" used to allocate the electrical energy consumption of the source-side circulation pump when the Studenthuset system is simultaneously providing building heating and DHW from the heat pumps while also providing cooling somewhere in the building. Figure 4 shows the heating fraction as a function of outdoor air temperature. Each point on the plot represents the average heating fractions for a particular outdoor air temperature bin of width 1 °C.
The heating system has a maximum capacity of 200 kW, but as seen in the duration curve in Figure 5, heating is mostly provided at 20-60 kW. Maximum capacity of the borehole free-cooling system is 120 kW.

System Boundaries
In this paper we have used the system boundaries for calculating system performance for heating and cooling, as suggested by the EU project SEPEMO [22]. The SEPEMO project was mainly focused on single-family houses with heat pumps for heating or cooling, while Studenthuset is an office building with both heating and cooling, and thus has a higher degree of system complexity. Figure 6 shows the SEPEMO system boundaries applied to Studenthuset.

System Boundaries
In this paper we have used the system boundaries for calculating system performance for heating and cooling, as suggested by the EU project SEPEMO [22]. The SEPEMO project was mainly focused on single-family houses with heat pumps for heating or cooling, while Studenthuset is an office building with both heating and cooling, and thus has a higher degree of system complexity. Figure 6 shows the SEPEMO system boundaries applied to Studenthuset. . SEPEMO system boundaries [22] for heating SPF (a) and cooling SPF (b) for Studenthuset. SPFC1 cannot be calculated, as the system has no cooling unit.
Tables 3 and 4 summarize the energy delivered, and the allocated electricity usage for calculation of the various SPFs for heating and cooling respectively as defined by the SEPEMO system boundaries. SPFH1 cannot be calculated for Studenthuset as there is no instrumentation to measure heat pump compressor electricity separately from the built-in circulation pumps. Tables 3 and 4 summarize the energy delivered, and the allocated electricity usage for calculation of the various SPFs for heating and cooling respectively as defined by the SEPEMO system boundaries. SPF H1 cannot be calculated for Studenthuset as there is no instrumentation to measure heat pump compressor electricity separately from the built-in circulation pumps.  Part of the energy used by the Legionella protection system will be transferred to the building as heat. We treated this as "useful heat" for every hour for which the heat pumps were delivering heat to the building. In hours when no heating was otherwise being provided, it was assumed to be not useful and was not counted as delivered heating. In Tables 3 and 4 an asterisk (*) added to the annual kWh denotes that the number has been estimated by allocating between heating and cooling and, for components that start with "useful heat dissipated", the heating is only considered useful if there is space heating being delivered by the heat pumps during the hour. The uncertainty due to the determination of usefulness is not included in the uncertainty estimate. "N/A" means that the data is not available. Exhaust air heat recovery (EAHR) contributes to the heat delivered on the load-side, though it cannot be quantified with the existing instrumentation. Double asterisks (**) in Table 3 indicate that the SEPEMO guidelines don't clearly speak to the question of whether or not heating provided by the EAHR system should be considered "auxiliary heating". The triple asterisks (***) in Table 3 for SPF H1 indicate that the calculated value does not conform exactly to the SEPEMO H1 boundaries, since the electrical energy includes all energy used by the heat pump units-including internal circulating pumps and control boards. Presumably, the actual SPF H1 should be somewhat higher.
As the Studenthuset cooling system does not include a cooling unit, there is no SPF C1 . Also, there is no supplementary cooling for Studenthuset, hence SPF C3 = SPF C4 .

Annual Performance
Annual delivered space heating, space cooling, and DHW heating from the GSHP system is shown in Figure 7. Heat recovered from the ventilation air rotary heat exchangers is not included in the figure, as this energy could not be estimated from available data. Nor is the amount of cooling provided directly from the ventilation system. Part of the energy used by the Legionella protection system will be transferred to the building as heat. We treated this as "useful heat" for every hour for which the heat pumps were delivering heat to the building. In hours when no heating was otherwise being provided, it was assumed to be not useful and was not counted as delivered heating. In Tables 3 and 4 an asterisk (*) added to the annual kWh denotes that the number has been estimated by allocating between heating and cooling and, for components that start with "useful heat dissipated", the heating is only considered useful if there is space heating being delivered by the heat pumps during the hour. The uncertainty due to the determination of usefulness is not included in the uncertainty estimate. "N/A" means that the data is not available. Exhaust air heat recovery (EAHR) contributes to the heat delivered on the load-side, though it cannot be quantified with the existing instrumentation. Double asterisks (**) in Table 3 indicate that the SEPEMO guidelines don't clearly speak to the question of whether or not heating provided by the EAHR system should be considered "auxiliary heating". The triple asterisks (***) in Table 3 for SPFH1 indicate that the calculated value does not conform exactly to the SEPEMO H1 boundaries, since the electrical energy includes all energy used by the heat pump units-including internal circulating pumps and control boards. Presumably, the actual SPFH1 should be somewhat higher.
As the Studenthuset cooling system does not include a cooling unit, there is no SPFC1. Also, there is no supplementary cooling for Studenthuset, hence SPFC3 = SPFC4.

Annual Performance
Annual delivered space heating, space cooling, and DHW heating from the GSHP system is shown in Figure 7. Heat recovered from the ventilation air rotary heat exchangers is not included in the figure, as this energy could not be estimated from available data. Nor is the amount of cooling provided directly from the ventilation system.   Annual electricity used on the source-side and load-side is shown in the right hand column in Figure 7. This electricity consumption used for heat recovery is lumped in with all of the distribution pump and fan energy, and cannot readily be quantified separately. Without accounting for the recovered heat in the ventilation air or the "free" cooling sometimes provided by the ventilation air, SPF H4 , SPF C3 , and SPF C4 cannot be fairly calculated. Figure 7 shows that electricity use for distribution of heating and cooling, ventilation, and heat recovery accounts for about twice as much electricity as the heat pumps and source-side circulating pump together. In a centralized system such as this, distribution energy on the load-side may be similar regardless of the chosen energy source. That is, the distribution energy may be approximately the same whether the heating is provided by centralized GSHPs, district heating, or a gas boiler. Therefore, SPF and COP for H3 and C2 boundaries may still be useful for comparison to other heating and cooling system types. However, for comparison to distributed GSHP systems, it is desirable to make comparisons using the H4 or C3 boundaries, because the distributed GSHP energy necessarily includes fan power.

Heating
Delivered space heating, electricity consumption for heating, and seasonal performance factors for the heating system are shown in Table 3. As discussed above, the electricity used by the source-side circulation pump had to be allocated between heating and cooling based on the heating fraction. The total annual space heating delivered to the building, excluding heat recovery, was 194 MWh which is consistent with the anticipated annual space heat load (200 MWh) for the building.
According to the building design documents, the SPF for the heat pumps was expected to be 4.5. As mentioned above, the design documents are silent on the point of which boundary conditions correspond to the stated SPF. Annual heating SPF for the heat pumps including electricity for the internal circulation pumps and control boards was 4.0 ± 0.2. As noted above, this value should be a little lower than SPF H1 , because of the additional electrical energy being included. Annual SPF H2 (which includes both the internal and external circulation pumps) was 3.7 ± 0.2. Taking into account the supplemental heating provided by the Legionella protection system, SPF H3 is 2.7 ± 0.13. The lower value is due to the supplemental heating being provided by an electric resistance heater and a recirculation pump.
Taking the annual heating SPF for the heat pumps including electricity for the internal circulation pumps as a surrogate for SPF H1 , we can see that SPF H1 and SPF H2 are 18% and 10% lower, respectively, than the median values reported in the literature review ( Table 2). The 18 and 16 cases for which SPF H1 and SPF H2 are reported are all for buildings situated in warmer climates with warmer ground temperatures, which may explain part of the reason for this. Other reasons for sub-optimal performance are discussed below.

Cooling
Seasonal performance factors for the cooling system are shown in Table 4. Studenthuset does not have a cooling unit, hence SPF C1 does not exist, but SPF C2 which includes the source and source-side circulation pump is estimated as 27 ± 5. Electricity used by source-side circulation pumps for cooling is allocated using the cooling fraction defined by Equation (6). As the amount of cooling provided by the ventilation system is not measured, SPF C3 and SPF C4 cannot be estimated.
Measured annual cooling amounts to approximately 130 MWh, which is about four times as much as the anticipated annual design cooling-load of 34 MWh. The reason for this difference has not yet been determined.
Compared to the nine cases for which SPF C2 is reported in the literature, the Studenthuset system has significantly higher SPF C2 , about 3.4 times the median value. This is primarily due to direct use of cold ground heat exchanger fluid. It may also be partly explained that the system is always providing some heating, so some of the pumping energy is always allocated to heating, thus reducing the energy consumed for cooling.

Heating
Monthly performance factors for heating were calculated using the system boundaries in Figure 6. Figure 8 shows monthly COP H2 and COP H3 for the evaluated period.
Monthly performance factors for heating were calculated using the system boundaries in Figure  6. Figure 8 shows monthly COPH2 and COPH3 for the evaluated period.
Somewhat counter-intuitively, the monthly heating COPs are highest during the colder months and lowest in the summer. The explanation lies partly in the source-side circulation pumps providing a minimum of 8 L/s throughout the year, regardless of whether or not there is a heating load. COPH3 is largely affected by the Legionella protection system (LPS), which runs all the time. Much of the electricity consumed by the LPS is eventually dissipated to the building, so during the heating season, this is counted as "useful heating delivered", increasing COPH3. However, during the warmer months, the heat dissipated is not useful and it significantly decreases COPH3 to the point that it falls below 1 for some months.
The uncertainty for boundary H3 is generally lower than that for H2. This is largely due to the fact that the electrical energy can be measured more accurately for boundary H3 than for H2. Somewhat counter-intuitively, the monthly heating COPs are highest during the colder months and lowest in the summer. The explanation lies partly in the source-side circulation pumps providing a minimum of 8 L/s throughout the year, regardless of whether or not there is a heating load. COP H3 is largely affected by the Legionella protection system (LPS), which runs all the time. Much of the electricity consumed by the LPS is eventually dissipated to the building, so during the heating season, this is counted as "useful heating delivered", increasing COP H3. However, during the warmer months, the heat dissipated is not useful and it significantly decreases COP H3 to the point that it falls below 1 for some months. The uncertainty for boundary H3 is generally lower than that for H2. This is largely due to the fact that the electrical energy can be measured more accurately for boundary H3 than for H2.

Cooling
Monthly cooling coefficients of performance for boundary C2 are presented in Figure 8. The COP varies between around 20 in the summer to above 60 during the winter, when the heat pumps are lowering the heat carrier fluid temperature. The amount of cooling used in the winter months is however quite small. These COPs do not account for distribution energy (as in COP C3 and COP C4 ). COP C3 cannot currently be calculated due to lack of data from the ventilation system. As COP C3 takes into account cooling provided by the ventilation system and the energy of the pumps and fans, it is difficult to speculate as to how values of COP C3 will compare to COP C2 .

Performance vs. Temperature and Load
In Figure 9 COP H3, H2, and C2, binned at half-degree intervals, are plotted versus entering fluid temperature (EFT) to the system from the borehole field. By "binned" we mean that average values are taken for all hours with heat pump EFT that fall within a half-degree bin, e.g., the points shown at 10 • C represent the total heating provided divided by the electricity used for all hours with heat pump EFT between 9.75 • C and 10.25 • C. Monthly cooling coefficients of performance for boundary C2 are presented in Figure 8. The COP varies between around 20 in the summer to above 60 during the winter, when the heat pumps are lowering the heat carrier fluid temperature. The amount of cooling used in the winter months is however quite small. These COPs do not account for distribution energy (as in COPC3 and COPC4). COPC3 cannot currently be calculated due to lack of data from the ventilation system. As COPC3 takes into account cooling provided by the ventilation system and the energy of the pumps and fans, it is difficult to speculate as to how values of COPC3 will compare to COPC2.

Performance vs. Temperature and Load
In Figure 9 COP H3, H2, and C2, binned at half-degree intervals, are plotted versus entering fluid temperature (EFT) to the system from the borehole field. By "binned" we mean that average values are taken for all hours with heat pump EFT that fall within a half-degree bin, e.g., the points shown at 10 °C represent the total heating provided divided by the electricity used for all hours with heat pump EFT between 9.75 °C and 10.25 °C. COPC2 and COPH2 have minimums at 11 °C and COPH3 drops significantly between 10-11.5 °C. As shown in Figure 10 the overall load is as lowest at those EFTs. This correlation between low part COP C2 and COP H2 have minimums at 11 • C and COP H3 drops significantly between 10-11.5 • C. As shown in Figure 10 the overall load is as lowest at those EFTs. This correlation between low part loads and low system performance is consistent with measured performance of the distributed GSHP system [51] serving an office building in Atlanta and also the central GSHP system [60] serving a university building in Leicester, England. This may seem counterintuitive, compared to thermodynamic theory, which might lead us to expect that the highest heating COP should occur at the highest heat pump EFT. But the effects of "parasitic" loads-e.g., pumping and control boards-and heat pump cycling losses seem to dominate the actual field-measured COP. For COP H3 , the very low values at temperatures above 10.5 • C are largely due to the Legionella protection system, as discussed in the previous section regarding performance during summer months.
Energies 2019, 12, x FOR PEER REVIEW 26 of 35 loads and low system performance is consistent with measured performance of the distributed GSHP system [51] serving an office building in Atlanta and also the central GSHP system [60] serving a university building in Leicester, England. This may seem counterintuitive, compared to thermodynamic theory, which might lead us to expect that the highest heating COP should occur at the highest heat pump EFT. But the effects of "parasitic" loads-e.g., pumping and control boardsand heat pump cycling losses seem to dominate the actual field-measured COP. For COPH3, the very low values at temperatures above 10.5 °C are largely due to the Legionella protection system, as discussed in the previous section regarding performance during summer months.  Figure 11 shows daily overall COP plotted versus overall energy load. Loads include building heating and DHW heated by the heat pumps and the Legionella protection system (LPS), while for cooling only the chilled water provided by the boreholes is included. The electricity use includes heat pump operation, LPS, and source-side circulating pump. Days with mainly cooling load (heating fraction ≤0.25), mainly heating load (heating fraction ≥0.67), and mixed heating and cooling load are color coded to illustrate the different correlations between COP and heating and cooling load respectively.  Figure 11 shows daily overall COP plotted versus overall energy load. Loads include building heating and DHW heated by the heat pumps and the Legionella protection system (LPS), while for cooling only the chilled water provided by the boreholes is included. The electricity use includes heat pump operation, LPS, and source-side circulating pump. Days with mainly cooling load (heating fraction ≤0.25), mainly heating load (heating fraction ≥0.67), and mixed heating and cooling load are color coded to illustrate the different correlations between COP and heating and cooling load respectively.
Energies 2019, 12, x FOR PEER REVIEW 26 of 35 loads and low system performance is consistent with measured performance of the distributed GSHP system [51] serving an office building in Atlanta and also the central GSHP system [60] serving a university building in Leicester, England. This may seem counterintuitive, compared to thermodynamic theory, which might lead us to expect that the highest heating COP should occur at the highest heat pump EFT. But the effects of "parasitic" loads-e.g., pumping and control boardsand heat pump cycling losses seem to dominate the actual field-measured COP. For COPH3, the very low values at temperatures above 10.5 °C are largely due to the Legionella protection system, as discussed in the previous section regarding performance during summer months.  Figure 11 shows daily overall COP plotted versus overall energy load. Loads include building heating and DHW heated by the heat pumps and the Legionella protection system (LPS), while for cooling only the chilled water provided by the boreholes is included. The electricity use includes heat pump operation, LPS, and source-side circulating pump. Days with mainly cooling load (heating fraction ≤0.25), mainly heating load (heating fraction ≥0.67), and mixed heating and cooling load are color coded to illustrate the different correlations between COP and heating and cooling load respectively. The figure raises some interesting questions. Firstly, why is there such a discrepancy between COP C2 in Figures 8 and 9, and the overall COP for the mainly cooling cases? The explanation is that for days when the heating fraction is very small, there is still a DHW load, and the LPS is still running. For an outdoor air bin temperature of 20 • C, the cooling provided is about 40 kW and the DHW provided is about 2.3 kW, but the LPS is using 3.2 kW and the circulating pump is about 1 kW. This means that even though COP C2 at those conditions is 42, adding the DHW heating gives about 42.3/4.2 kW, which gives an overall COP of only 10.1.
Second, we don't see nearly as strong a trend for the heating loads as for the cooling loads. COPs for cooling dominated conditions increase about 10 units over 1000 kWh increased daily energy load, whereas the increase in heating dominated COPs is only 1 unit over 2500 kWh load increase. While an increased cooling load only requires a small increase in electricity use for the circulation pump to increase the flow rate, it also means that the relative influence on the COP from the Legionella protection system and DHW circulation decreases. For heating dominated conditions, the Legionella protection system and DHW circulation are small contributions to the overall load and electricity use. Also, since there are five heat pumps that are staged, individual heat pumps tend to work at close to full load conditions.

Discussion
System performance has been analyzed based on boundary conditions H2, H3, and C2 shown in Tables 3 and 4. The results have been presented for a one-year period, monthly periods, daily periods, hourly periods, and binned hourly periods.
Uncertainties in heating/cooling provided, COPs, and SPFs are provided. Because uncertainties for performance parameters based on field measurements have rarely been presented in the literature, there is little to compare these values to, with the exception of Southard et al. [49], Urchueguía et al. [37], and Hughes [41]. The uncertainties in the present study are similar in magnitude or smaller to those presented in either Southard et al. [49], Urchueguía et al. [37], or Hughes [41]. The present study has the advantage of using matched pairs of temperature sensors that significantly reduce the uncertainty in the temperature difference measurement. However, the Southard et al. [49] study had the advantage of making separate measurements on 14 different heat pumps; aggregation of the different measurements reduces the uncertainty of the total heating or cooling provided by all heat pumps. It would be desirable to check the accuracy of the individual sensors and sensor pairs after several years of operation in order to insure that the sensors have not drifted from their calibration.
For the situation in Sweden, the annual SPFs SPF H2 or SPF H3 are commonly used as a basis for comparison, because other systems to which GSHP systems may be compared are also centralized systems. That is, they also provide hot water and chilled water to panel radiators, fan coil units, and other heating/cooling distribution devices. In this case, it makes sense to compare the "central plant" equipment using the H2 and H3 boundary conditions. For cooling, the C2 boundary conditions would be commonly used.
However, if we wish to compare to distributed GSHP systems or other distributed systems where the heating and cooling equipment has integrated distribution fans and, likely, ventilation provided using these same fans, it is necessary to compare using boundary conditions that include the heating and cooling distribution and even the ventilation. In our results for the Studenthuset system, these would correspond to the H4 and C4 results.
This, however, raises a further complication. The ventilation system does provide heating and cooling because outdoor air is at least tempered using hot water or chilled water from the heat pumps. Therefore, the ventilation system electricity consumption is included in the H4 and C4 boundary conditions. On the other hand, the ventilation system consumes electricity primarily for the purpose of providing fresh air to the space, not for heating and cooling purposes. So, when calculating heating and cooling COPs with the H4 and C4 boundary conditions, is there a fair way to divide the electricity consumption between heating/cooling and provision of fresh air? This has not been addressed in the literature. It seems that it might be possible to calculate the pressure drop due to heating and cooling coils, the pressure drops due to the rest of the ventilation system, and apportion the electricity consumption accordingly. This remains a topic for further investigation.
Beyond the question of the comingling of electricity for ventilation and heating/cooling, several other problematic issues have arisen during the performance analysis of the Studenthuset system. These include: • Division of electricity consumption between heating and cooling when both are being provided simultaneously. For the Studenthuset building, this arises when accounting for the electricity used by the source-side circulating pump under conditions when one or more heat pumps are being operated to provide heating and "free" cooling is being provided at the same time directly from the boreholes. It also arises with respect to the fans and pumps of the heating/cooling distribution system and ventilation system, for which we only have a single energy meter. In these cases, we divided the electricity proportionally to the amount of heating and cooling provided. • The Studenthuset system utilizes an electric resistance heater to raise the domestic hot water temperature in order to protect against Legionella. The required hot water temperature is 60 • C, but in practice the electric resistance heater is operated continuously and excess heat is dissipated into the building. A pump is used to continuously circulate this hot water around the building. This continuous circulation is used both to protect against Legionella (by insuring that temperatures remain high throughout the domestic hot water piping) and to significantly reduce the amount of waiting time for hot water at any tap. Treatment of this subsystem is perplexing. Should the heating provided by the electric resistance heater all be treated as providing useful heating to the building? Clearly, when domestic hot water is being utilized, the Legionella protection system is providing some of the required heating. The energy expended to heat the domestic hot water to 60 • C can certainly be considered as part of the heating load. But what about that dissipated to the building? It could either be providing useful heating or merely increasing the cooling load. So, it seems that it should be treated on the basis as to whether or not the building has a heating load at that time. In this analysis, we took this into account by treating the heat provided by the Legionella protection system as useful for any hours where the heat pumps were providing any building heating. However, as this includes hours when the building was simultaneously being heated and cooled, this approach should be refined.

•
The fact that the Legionella protection system ran every single hour did, however, alleviate one other problem encountered in Southard et al. [49]-how to treat energy consumed by pumps, fans, control boards, etc., when neither heating nor cooling are provided by the system. Hughes [41] also identifies Legionella protection systems as problematic for the overall system performance, as they significantly increase the use of auxiliary heating with immersion heaters. He observes that several of the GSHP systems in his study do not appear to be operating according to current best practice. • It may always be the case that additional instrumentation is desirable. In this case, we could not calculate COP H1 because the heat pumps had integrated circulating pumps and we did not have independent measurements of the compressor energy consumption.

•
Another limitation in the instrumentation is that heating supplied and cooling supplied are determined solely from water-side measurements made near the central plant. As shown in Figure 7, the internal circulation pumps and fans use a little over half of all system electricity. This electricity usage would be accounted for in the input electricity used to calculate COP and SPF for boundaries H4, C3, and C4. However, all of this electricity is eventually dissipated as heat somewhere. Treatment of this dissipated heat faces the same complications as the treatment of heat from the Legionella protection system. Furthermore, the fan energy is used to provide heating and cooling to the ventilation air and cooling to the building space. But without suitable instrumentation, it is not possible to quantify this. Therefore, we demurred from calculating COP and SPF for boundaries H4, C3, and C4. Further instrumentation here would be particularly helpful.
Finally, one of the most important uses of this kind of performance investigation is identification of possible system improvements. We have identified several possibilities, though their impacts have not yet been fully evaluated: • The scheduled run time of the electric resistance heater providing Legionella protection could be reduced. It currently runs 24 h per day, seven days per week. Many heat pump systems run a heating cycle for Legionella protection on a once-per-week schedule. Wemhoener et al. [33] describe a Legionella protection system running twice per week. Hughes [40,41] describes a Legionella protection method where the temperature at the top of the DHW storage tank is kept at 60 • C, and an immersion heater is used to raise the temperature in the whole tank to 60 • C for a minimum of one hour per day. One of the GSHP systems reported by Hughes [40,41] uses ultraviolet light for disinfection; however, that system already has a high DHW output temperature of above 62 • C.

•
The continuous circulation pump could also be scheduled to run at a reduced run time, more closely matching the operating hours of the building. Merely limiting the circulation pump and electric resistance heating to run during office hours only would save about 12,000 kWh/year in electricity consumption. Some of the heating provided by the Legionella protection system would then be replaced with heat provided by heat pumps at a higher COP.

•
The minimum flow rate through the ground heat exchanger is set to be 8 L/s year-round. This often results in very low ∆T across the boreholes. It seems likely that the setpoint could be reduced significantly for portions of the year and/or a reset control could be utilized that reduces the flow when allowable.

Conclusions
The Studenthuset in Stockholm, Sweden, and its GSHP system were designed with high ambitions for both energy efficiency and indoor comfort. The owner has long experience with GSHP systems in buildings and the system is thoroughly instrumented and monitored by experienced staff. One notable gap in the instrumentation (airflow rate measurement) precludes calculating COP and SPF for the entire system. This study shows that the system provides space heating consistent with the design values, and that the cooling provided is about four times higher than anticipated from the design.
An SPF for SEPEMO boundary H2 (heat pump + source-side circulating pump) of 3.7 ± 0.2 was achieved. An SPF for boundary H3 (heat pump + source-side circulating pump + auxiliary heating) of 2.7 ± 0.13 was achieved. The SPF H2 value is lower than most of the reported values in the literature, which have a median value of 3.6. The SPF H3 value is lower than the only reported value of SPF H3 in the literature. For cooling, the SPF C2 value of 27 ± 5 is considerably higher than any value reported in the literature. These values are partly explained by the relatively lower ground temperatures found in Stockholm.
However, a key finding is that the measured COPs are more affected by the amount of heating and cooling provided than by the entering fluid temperature to the heat pumps. As shown by the bin analysis, lower entering fluid temperatures correspond to higher run time fractions for equipment and less influence of "parasitic" loads such as pumps and unit control boards. This is consistent with findings for both a distributed GSHP system [51] and a central GSHP system [60]. For this system, the Legionella protection system and DHW recirculation system, which are always on, form a significant parasitic load. The near-constant flow rate in the borehole circuit leads to the circulating pump also acting as a significant parasitic load. Therefore, reduced scheduling of the Legionella protection and DHW recirculation systems and allowing reduced flow in the borehole circuit are recommended actions to significantly decrease the parasitic losses, improve the performance, and decrease the energy performance gap. This paper presents an approach to addressing the performance gap without installing research-grade instrumentation. High-quality commercially-available instrumentation selected by the owners for the building energy management system proved sufficient to measure the system SPFs and COPs for several of the SEPEMO boundary conditions with reasonable accuracy, as demonstrated with the uncertainty analysis. Another contribution of this paper is the use of a bin analysis of system COPs to identify several problems with the operating strategy that, once corrected, have significant potential to reduce the system energy consumption.
Despite the thorough instrumentation and monitoring, several additional measurement points are needed in order to provide a complete analysis of the system performance. Additional instrumentation and data logging would be needed for estimating the amount of heat recovery and cooling provided from the ventilation air, as well as allocating the electricity used for this. This would make possible estimation of SPF H4 and SPF C3 and SPF C4 .
Other desirable additional instrumentation would be electricity measurements for the heat pump compressors to enable estimation of SPF H1 , separate electricity measurement for the circulation pump on the source-side, and separate electricity measurements for the different parts of the Legionella protection system.
For a Swedish centralized GSHP system like Studenthuset, it makes sense to evaluate SPF H2, SPF H3, SPF C2 , and an overall SPF including only the ground source, heat pumps, and DHW heating, and leaving out the load-side energy and electricity. However, since SPF H4 and SPF C3 (=SPF C4 ) could not be calculated for Studenthuset, it is not possible to compare this centralized system with a decentralized system such as the ASHRAE HQ in Atlanta.
Finally, this study has identified several issues in calculation of SPF and COP for commercial buildings that are not addressed in the literature or by existing boundary definition schema. For example, treatment of electricity used by ventilation fans, for which the main purpose is providing ventilation, but which also provide significant amounts of heating and cooling. Further research in this area is part of an ongoing collaborative research project [66] organized by the International Energy Agency Heat Pumping Technologies Collaboration Program. Acknowledgments: Thanks to Anders Larsson and Åke Annsberg at Akademiska Hus for their help and support, and providing measurement data and information about the monitored building. Special thanks to Akademiska Hus for granting permission to publically share the data. The support from the authors' employers is gratefully acknowledged. This work is part of the IEA HPT Annex 52, long term performance measurement of GSHP systems serving commercial, institutional, and multi-family buildings [66]. Interactions with our international colleagues have been helpful in influencing our thinking on addressing the performance gap with monitoring.

Conflicts of Interest:
The authors declare no conflict of interest. Nomenclature c p specific heat of the water, kJ/(kg·K) COP H1 Coefficient of performance over a given period of time (e.g., a day or a month) for the system boundary heat pump only COP H2 Coefficient of performance over a given period of time (e.g., a day or a month) for the system boundary heat pump and circulation pumps/fans on source-side COP H3 Coefficient of performance over a given period of time (e.g., a day or a month) for the system boundary heat pump, circulation pumps/fans on source-side and auxiliary heating COP H4 Coefficient of performance over a given period of time (e.g., a day or a month) for the system boundary heat pump, circulation pumps/fans on source-side, auxiliary heating and circulation pumps/fans on load-side COP C1 Coefficient of performance over a given period of time (e.g., a day or a month) for the system boundary cooling unit only COP C2 Coefficient of performance over a given period of time (e.g., a day or a month) for the system boundary cooling unit and circulation pumps/fans on source-side COP C3 Coefficient of performance over a given period of time (e.g., a day or a month) for the system boundary cooling unit, circulation pumps/fans on source-side and load-side COP C4 Coefficient of performance over a given period of time (e.g., a day or a month) for the system boundary cooling unit, circulation pumps/fans on source-side and load-side, and supplementary cooling unit e flow fractional uncertainty of the flow rate e ∆T fractional uncertainty of the temperature (K) e d fractional uncertainty of the denominator e n fractional uncertainty of the numerator e qt fractional uncertainty of the quotient e P fractional uncertainty in the meter measurement e q fractional uncertainty of the heat transfer rate E ∆T absolute uncertainty of the temperature (K) E P-H2 uncertainty in the electrical power and energy consumed by the heat pumps and source-side circulating pump under heating conditions for use in calculating COP H2 E P-H3 uncertainty in the electrical power and energy consumed by the heat pumps and source-side circulating pump under heating conditions for use in calculating COP H3 E q absolute uncertainty in building heating provided for a one-month period (kWh) E SSCP absolute error in the estimate of the source-side circulating pump electrical energy consumption, kWh/h f clg fraction of electrical energy consumed by the circulating pump that is allocated to the heating electrical energy f htg fraction of electrical energy consumed by the circulating pump that is allocated to the heating electrical energy n day of the year n min day of the year with minimum temperature (day 60) P electrical energy measured by the meter, kW P SSCP electrical power (kW) consumed as a function of flow rate q BldgHtg average heating provided to the building over the time interval, kW q BldgClg average cooling provided to the building over the time interval, kW q DHW average heating provided to the domestic hot water over the time interval, kW Q DHW average heating energy (kWh) provided by the heat pumps to the domestic hot water Q LPSRH energy provided by the electric resistance heater for the Legionella protection system SPF H1 Seasonal performance factor over 365 consecutive days for the system boundary H1 (heat pump only) SPF H2 Seasonal performance factor over 365 consecutive days for the system boundary H2 (heat pump and circulation pumps/fans on source-side) SPF H3 Seasonal performance factor over 365 consecutive days for the system boundary H3 (heat pump, circulation pumps/fans on source-side and auxiliary heating) SPF H4 Seasonal performance factor over 365 consecutive days for the system boundary H4 (heat pump, circulation pumps/fans on source-side, auxiliary heating, and circulation pumps/fans on load-side) SPF C1 Seasonal performance factor over 365 consecutive days for the system boundary C1 (cooling unit only) SPF C2 Seasonal performance factor over 365 consecutive days for the system boundary C2 (cooling unit and circulation pumps/fans on source-side) Energies 2019, 12, 2020 31 of 34 SPF C3 Seasonal performance factor over 365 consecutive days for the system boundary C3 (cooling unit, circulation pumps/fans on source-side and load-side) SPF C4 Seasonal performance factor over 365 consecutive days for the system boundary C4 (cooling unit, circulation pumps/fans on source-side and load-side, and supplementary cooling unit) T w incoming water temperature, • C T w,OHP temperature of heated water at outlet from heat pumps, • C T w,ORH temperature of heated water at outlet from resistance heater, • C V integrated volume flow over the time interval, m 3 . V volume flow rate averaged over the time interval, L/s ρ density of the water, kg/m 3