Abstract. The classification of X-ray solar flares is performed regarding their effects on the Very Low Frequency (VLF) wave propagation along the Earth-ionosphere waveguide. The changes in propagation are detected from an observed VLF signal phase and amplitude perturbations, taking place during X-ray solar flares. All flare effects chosen for the analysis are recorded by the Absolute Phase and Amplitude Logger (AbsPal), during the summer months of 2004-2007, on the single trace, Skelton (54.72 N, 2.88 W) to Belgrade (44.85 N, 20.38 E) with a distance along the Great Circle Path (GCP) D≈2000 km in length.The observed VLF amplitude and phase perturbations are simulated by the computer program Long-Wavelength Propagation Capability (LWPC), using Wait's model of the lower ionosphere, as determined by two parameters: the sharpness (β in 1/km) and reflection height (H in km). By varying the values of β and H so as to match the observed amplitude and phase perturbations, the variation of the D-region electron density height profile N e (z) was reconstructed, throughout flare duration. The procedure is illustrated as applied to a series of flares, from class C to M5 (5×10 −5 W/m 2 at 0.1-0.8 nm), each giving rise to a different time development of signal perturbation.The corresponding change in electron density from the unperturbed value at the unperturbed reflection height, i.e. N e (74 km)=2.16×10 8 m −3 to the value induced by an M5 class flare, up to N e (74 km)=4×10 10 m −3 is obtained. The β parameter is found to range from 0.30-0.49 1/km and the reflection height H to vary from 74-63 km. The changes in N e (z) during the flares, within height range z=60 to 90 km are determined, as well.