AbstractSolar cycle variation of earth's atmosphere, particularly the ionosphere is of particular interest in the field of space science and space weather studies. In this article, we present the outcome of our detailed quantitative study on solar cycle variation of lower ionospheric properties using numerical investigation. First, we seek to model and compare the collective D‐region ionization rates (’s) due to the individual contributions from the ionizing sources, namely, (a) solar extreme ultraviolet (EUV) radiation including the Lyman‐ irradiation and (b) solar X‐ray irradiation throughout the solar cycle (C24). Then, we compute the electron density profiles using the ionization rates. We report significant solar cycle variation in ionization rate and electron density profiles for the entire span of C24. We use the sunspot number (SPN) profile to substantiate the finer details of profile during C24. For the segments of the D‐region above two different geographic coordinates, we report that profiles show a consistent “dual peak” nature. It is very similar to the SPN profile during C24. As a next‐order validation, we compare our modeled profiles with their IRI‐2020 counterpart (s). Their overall trends are found to be in agreement. Finally, we discuss the response of the D‐region in terms of due to C24. The work lays the foundation for our upcoming studies on D‐region response to solar cycle variation with Very Low Frequency (VLF) observation.