Retrieval of Sunrise C-Region Electron Density Using Mid-Range VLF Amplitude and FDTD-Based Optimization

This study presents a method to retrieve the electron density structure of the transient C-region using very-low-frequency (VLF) Earth–ionosphere waveguide propagation. Here, we demonstrate the identification of the C-region from amplitude variations of a mid-range VLF propagation path that is nearly perpendicular to the solar terminator. Previous investigations have primarily relied on phase measurements along long-distance paths with small terminator angles, whereas the present approach utilizes amplitude information under conditions where modal interference is significant. The Faraday International Reference Ionosphere (FIRI-2018) provides an effective semi-empirical model of the lower-ionospheric electron density; however, discrepancies between simulations and observations are often observed at sunrise. To resolve this issue, we introduce Gaussian perturbations to the electron density profile output by FIRI-2018 and optimize their parameters so that finite-difference time-domain (FDTD) simulations reproduce the observed VLF amplitude. The analysis is performed for the 22.2 kHz JJI transmitter signal received in Chofu, Japan over a mid-range propagation path, ∼900 km. The optimized electron density profile successfully reproduces the characteristic features of the C-region, including a temporary enhancement near 65 km altitude during sunrise. These results demonstrate that mid-range VLF amplitude analysis provides a quantitative tool for identifying transient lower- ionospheric structures.