The use of shallow geothermal energy increasingly receives attention as a suitable alternative to fossil fuel-based space heating and cooling, warm water provision, as well as for seasonal heat storage throughout Europe. With the advent of shallow geothermal energy use on large scales, a vivid discussion of potential ecological and economic impacts has arisen but actual field data are scarce. An intensive groundwater temperature-monitoring program over a period of 3 years with consecutive measurements was, therefore, initiated at a residential neighborhood in the city of Cologne, Germany, under intense shallow geothermal use. The aim of the monitoring program was to overcome the existing data scarcity by pinpointing the effects of the intensive thermal use of the subsurface on groundwater temperatures and to foster understanding of urban groundwater temperature evolution. Results show that even though energy demands of the individual houses and energy extraction rates of the shallow geothermal systems were comparably small in this case, the accumulation of shallow geothermal users had a measurable impact on overall groundwater temperatures.

Limited knowledge about the spatial distribution of aquifer properties typically constrains our ability to predict subsurface flow and transport. Here we investigate the value of using high resolution full‐waveform inversion of cross‐borehole ground penetrating radar (GPR) data for aquifer characterization. By stitching together GPR tomograms from multiple adjacent crosshole planes, we are able to image, with a decimeter scale resolution, the dielectric permittivity and electrical conductivity of an alluvial aquifer along cross sections of 50 m length and 10 m depth. A logistic regression model is employed to predict the spatial distribution of lithological facies on the basis of the GPR results. Vertical profiles of porosity and hydraulic conductivity from direct‐push, flowmeter and grain size data suggest that the GPR predicted facies classification is meaningful with regard to porosity and hydraulic conductivity, even though the distributions of individual facies show some overlap and the absolute hydraulic conductivities from the different methods (direct‐push, flowmeter, grain size) differ up to approximately one order of magnitude. Comparison of the GPR predicted facies architecture with tracer test data suggests that the plume splitting observed in a tracer experiment was caused by a hydraulically low‐conductive sand layer with a thickness of only a few decimeters. Because this sand layer is identified by GPR full‐waveform inversion but not by conventional GPR ray‐based inversion we conclude that the improvement in spatial resolution due to full‐waveform inversion is crucial to detect small‐scale aquifer structures that are highly relevant for solute transport.