With the increasing installation of shallow geothermal energy, the importance of thermal impact prediction also increases in the system design stage. In nature, it is general that heterogeneity exists and it can affect the groundwater flow as well as the transport along to the flow. When predicting heat transport under the groundwater heat pump (GWHP) operation, however, impacts of heterogeneity have rarely been considered. In this study, to detect the hydraulic and thermal feedback to the two months of GWHP operation, a dense monitoring network was constructed with 12 monitoring wells at Eumseong-gun, Republic of Korea. The temperature was monitored in high resolution via fiber-optic distributed temperature sensing. During the GWHP operation, a very dynamic flow condition was generated with the hydraulic gradient between 0.005 and 0.07. The maximum temperature change at the nearest monitoring well was 2 ?. Observed hydraulic and thermal responses showed spatially heterogeneous results. While the heterogeneous responses of hydraulic change were stronger near the geothermal wells, those of temperature change were higher near the center of the thermal plume.