The evaluation of chemical transport models, CTMs, is essential for the assessment of their performance regarding the physical and chemical parameterizations used. While regional CTMs have been widely used and evaluated over Europe, their validation over Greece is limited. In this study, we investigate the performance of the Long Term Ozone Simulation European Operational Smog (LOTOS-EUROS) v2.2.001 regional chemical transport model in simulating nitrogen dioxide, NO.sub.2, over Greece from June to December 2018. In situ NO.sub.2 measurements obtained from 14 stations of the National Air Pollution Monitoring Network are compared with surface simulations over the two major cities of Greece, Athens and Thessaloniki. Overall the LOTOS-EUROS NO.sub.2 surface simulations compare very well to the in situ measurements showing a mild underestimation of the measurements with a mean relative bias of â¼-10 %, a high spatial correlation coefficient of 0.86 and an average temporal correlation of 0.52. The CTM underestimates the NO.sub.2 surface concentrations during daytime by â¼-50 Â± 15 %, while it slightly overestimates during night-time â¼ 10 Â± 35 %. Furthermore, the LOTOS-EUROS tropospheric NO.sub.2 columns are evaluated against ground-based multi-axis differential optical absorption spectroscopy (MAX-DOAS) NO.sub.2 measurements in Athens and Thessaloniki. We report that the CTM tropospheric NO.sub.2 column simulations over both urban and rural locations represent the diurnal patterns and hourly levels for both summer and winter seasons satisfactorily. The relative biases range between â¼ -2 % and -35 %, depending on season and relative NO.sub.2 load observed. Finally, the CTM was assessed also against space-borne Sentinel-5 Precursor (S5P) carrying the Tropospheric Monitoring Instrument (TROPOMI) tropospheric NO.sub.2 observations. We conclude that LOTOS-EUROS simulates extremely well the tropospheric NO.sub.2 patterns over the region with very high spatial correlation of 0.82 on average, ranging between 0.66 and 0.95, with negative biases in the summer and positive in the winter. Updated emissions for the simulations and model improvements when extreme values of boundary layer height are encountered are further suggested.