Solvent impact on the structural index of aromaticity was modelled by polarised continuum field approximation (IEFPCM) and hybrid quantum chemistry (QM/MM) method. Significant solvent related relaxation of the solutes geometries were noticed especially for highly polar species. The significant reduction of the aromaticity was observed for some aromatic compounds in water solution compared to gas phase. The rationale of this fact was provided based on dipole moment changes, energy penalty for polarisation of solute and the distribution of frontier orbital densities. The incoherent predictions of explicit and implicit solvation models are noticed since in some cases the PCM approach artificially exaggerate the geometry relaxation in solution which is not observed if explicit solvent molecules are taken into account.