The Benchmark data set collected within the European COST Action ES1206 has aimed to support the development and validation of advanced Global Navigation Satellite System (GNSS) tropospheric products, in particular high-resolution zenith delays and tropospheric gradients. In this work we utilize this unique data set to show that the interpolation of GNSS Zenith Wet Delays (ZWDs) can be improved by utilizing tropospheric gradients. To do this we first prove the concept with simulated observations, that is, zenith delays and tropospheric gradients derived from a Numerical Weather Model. We show how tropospheric gradients can be converted to ZWD gradients. Then the ZWD gradients together with the ZWDs at selected reference stations are used in an inverse distance weighting interpolation scheme to estimate the ZWD at some target station. For a station configuration with an average station distance of 50 km in Germany and a period of two months (May and June 2013), we find an improvement of 20% in interpolated ZWDs when tropospheric gradients are taken into account. Next, we replace the simulated by real observations, that is, zenith delays and tropospheric gradients from a Precise Point Positioning (PPP) solution provided with the G-Nut/Tefnut analysis software. Here we find an improvement of 10% in interpolated ZWDs when tropospheric gradients are taken into account.