This study predicts the subsurface temperature distribution of Germany's capital Berlin. For this purpose, a data-based lithosphere-scale 3D structural model is developed incorporating 21 individual geological units. This model shows a horizontal grid resolution of (500 x 500) m and provides the geometric base for two different approaches of 3D thermal simulations, (i) calculations of the steady-state purely conductive thermal field and ( Groundwater flow results in subsurface cooling the extent of which is strongly controlled by the geometry and the distribution of the Tertiary Rupelian Clay. The cooling effect is strongest where this clay-rich aquitard is thinnest or missing thus facilitating deep reaching forced convective flow.The differences between the purely conductive and coupled models highlight the need for investigations of the complex interrelation of flow-and thermal fields to properly predict temperatures in sedimentary systems.