The method used in this study consists of evaluating different terminal heating/cooling systems in an experimental office at the Polytechnic University of Timișoara, Romania, through experimentation and numerical simulation. This paper is focused on the energy efficiency exploratory research of various low- and very low-temperature hydronic radiant heating/cooling systems (floor, wall, ceiling, combined floor-ceiling) and of the medium-temperature radiator heating system, connected to a vertical ground-coupled heat pump (GCHP) in double U-tube borehole heat exchanger (BHE) configuration. The main performance parameters are obtained for five weeks of operation in heating mode and three weeks in cooling mode, and a comparative analysis of these parameters together with indoor thermal comfort is carried out. The four simple heating systems (radiator, floor, wall, ceiling) have relatively small differences (maximum 7.4 %) in their coefficient of performance (COP) values, but the number of starts/stops in the case of radiator heating is twice that of floor heating, which leads to greater wear and tear on the heat pump equipment. Under the same operating conditions, the CO2 emission of the radiator heating system are 0.7–16 % higher than those of the heating system with radiant panels (floor, wall, and ceiling). The floor-ceiling combined radiant heating system has the highest COPsyst of 5.45 as well as the lowest CO2 emission value of 2.15 kg. The predicted mean vote (PMV) index (pair 1.1 met-0.29 clo) has values approximately equal to 0 only in the case of the floor-ceiling heating, when the percent of likely dissatisfied with the thermal comfort is 5 %. The three cooling radiant systems show differences in their energy performance, namely the radiant ceiling cooling system has the best COPsyst, 39.8 % higher than that of the radiant floor and only 5.9 % higher than that of the radiant wall. The values of the PMV index (pair 1 met-0.9 clo) for the radiant wall cooling system are 31–41 % lower than those of the radiant floor and 10.4–14.2 % lower than those of the radiant ceiling. Finally, numerical simulations of the useful thermal energy and the system COP in heating and cooling modes are performed with TRNSYS software for a duration of 8760 h and then analysed and compared with the experimental tests to validate the simulation models.