A large part of the path of the Annular Solar Eclipse of May 20,2012 (magnitude 0.9439) (ASE-2012) was over southern Japan. The D-region ionospheric changes associated with the ASE-2012, led to several degree of observable perturbations of sub-ionospheric very low frequency (VLF) radio signal. The solar eclipse associated signal changes were identified in VLF several receiving stations (Rx) simultaneously for the VLF signals coming from both Japanese and US VLF transmitters (Tx). In this work, we have analyzed temporal dependences of VLF electric amplitude perturbation (Δ Aecl,obs(t)) from two Japanese VLF transmitters (JJI (22.2kHz) and JJY (40.0kHz)), and the spatio-temporal characteristics of respective subionospheric perturbations are studied in detail. We consider the 2-parameter D-region ionospheric model with the exponential electron density profile. To model the shadow effect on the D-region ionosphere due to obscuration of solar disk, we assume a generalized space-time dependent 2-Dimensional Elliptical Gaussian distribution Model (2DEGM) for ionospheric parameters, such as, effective reflection height (h') and sharpness factor (β). In the vicinity of the eclipse zone, we compute the subionospheric VLF signal propagation for several signal propagation paths. In the simulation, we obtain the perturbation of VLF signal amplitude (Δ Aecl,LWPC(t)) at each station and compare with its observtaional counterpart (Δ Aecl,obs(t)). © 2019 URSI. All rights reserved.