Transient localized D region disturbances, such as those associated with lightning discharges, affect the characteristics of VLF waves propagating in the Earth-ionosphere waveguide. In particular, both phase and amplitude changes in the subionopsheric signal can be observed at receiving sites as a result of the wave scattering that takes place in the disturbed region. In the present paper we present a multiple-mode three-dimensional model of VLF propagation in the Earth-ionopshere waveguide in the presence of localized D region disturbances. The model takes into account great circle (GC) propagation paths with realistic ground and ionospheric conductivity changes that result in mode conversion along the path. It is assumed that conductivity changes transverse to the GC paths are negligible except in the vicinity of the D region disturbance and that mode coupling is negligible within the disturbed region. This new model is applied to experimental observations and is found to be in general agreement. The diagnostics potential of the model for characterizing energetic particle precipitation events is discussed. ¿American Geophysical Union 1993