The derivation of global navigation satellite systems (GNSSs) tropospheric products is nowadays a state-of-the-art technique that serves both research and operational needs in a broad range of applications in meteorology. In particular, GNSS zenith tropospheric delay (ZTD) data assimilation is widely applied in Europe to enhance numerical weather predictions (NWPs). The current study presents the first attempt at introducing assimilation of ZTDs, derived from more than 48 stations of the Hellenic GNSS network, into the operational NWP system of the National Observatory of Athens (NOA) in Greece, which is based on the mesoscale Weather Research and Forecasting (WRF) model. WRF was applied during seven high-impact precipitation events covering the dry and wet season of 2018. The simulation employing the ZTD data assimilation reproduces more accurately, compared to the control experiment, the observed heavy rainfall (especially for high precipitation events, exceeding 20 mm in 24h) during both dry and wet periods. Assimilating ZTDs also improves the simulation of intense (>20 mm) convective precipitation during the time window of its occurrence in the dry season, and provides a beneficial influence during synoptic-scale events in the wet period. The above results, which are statistically significant, highlight an important positive impact of ZTD assimilation on the model's precipitation forecast skill over Greece. Overall, the modelling system's configuration, including the assimilation of ZTD observations, satisfactorily captures the spatial and temporal distribution of the observed rainfall and can therefore be used as the basis for examining further improvements in the future. meteorological applications, including nowcasting and numerical weather prediction (NWP) [8][9][10][11], as well as weather monitoring, including extreme events [12][13][14]. Special meteorological interest derives from the near-real time (NRT) ZTDs, which are estimated based on raw GNSS observations. The ZTD is a standard GNSS product expressing the total signal delay in the zenith direction above a receiver [6,15]. This vertical lag contains information on the total columnar amount of water vapor [16].In Europe, collaborative scientific efforts over the past two decades substantially contributed to the development of networks and analysis centers collecting and processing, respectively, GNSS data to compute tropospheric delays. The establishment of the European GNSS water vapor program (E-GVAP; [17]) in 2005 allowed for the operational distribution of NRT ZTD estimates to the meteorological community [6]. This service encouraged the implementation of precipitation forecast impact studies involving the assimilation of NRT ZTD observations into NWP models. Poli et al. [18] found an improvement in the prediction of precipitation patterns over France during spring and summer when assimilating ZTD data into a global four-dimensional variational (4D-var) assimilation and forecasting system. Positive impacts of ZTD data assimilation on precipita...