The objective of this thesis is to use radio emissions from lightning, known as `radio atmospherics' or `sferics', to study the temporal and spatial variation of the lower ionosphere, a layer of ionized atmosphere beginning at $\sim$70 km altitude (D-region). Very Low Frequency (VLF, 3$-$30kHz) radio waves are a useful diagnostic for lower ionospheric monitoring due to their reflection from this region and global propagation. Traditionally, the lower ionosphere has been sensed using single-frequency VLF transmitters allowing for analysis of a single propagation path, as there are only a small number of transmitters. A lightning stroke, however, releases an intense amount of impulsive broadband VLF radio energy in the form of a sferic, which propagates through the Earth-ionosphere waveguide. Lightning is globally distributed and very frequent, so a sferic is therefore also a useful diagnostic of the D-region. This is true both for ambient or quiet conditions, and for ionospheric perturbations such as solar flare x-ray bursts. Lightning strokes effectively act as separate VLF transmitting sources. As such, they uniquely provide the ability to add a spatial component to ionospheric remote sensing, in addition to their broadband signature which cannot be achieved with man-made transmitters. We describe the methods of processing in detail. As an example, we analyze a solar flare during which time there is a significant change in magnitude and frequency content of sferics. This disturbance varies with distance from the source, as well as time. We describe the methods of processing in detail, and show results at Palmer Station, Antarctica for both a quiet and active solar day.

<p>The high power, broadband very low frequency (VLF, 3--30 kHz) and extremely low frequency (ELF, 3--3000 Hz) electromagnetic waves generated by lightning discharges and propagating in the Earth-ionosphere waveguide can be used to measure the average electron density profile of the lower ionosphere (<italic>D</italic> region) across the wave propagation path due to several reflections by the upper boundary (lower ionosphere) of the waveguide. This capability makes it possible to frequently and even continuously monitor the <italic>D</italic> region electron density profile variations over geographically large regions, which are measurements that are essentially impossible by other means. These guided waves, usually called atmospherics (or sferics for short), are recorded by our sensors located near Duke University. The purpose of this work is to develop and implement algorithms to derive the variations of <italic>D</italic> region electron density profile which is modeled by two parameters (one is height and another is sharpness), by comparing the recorded sferic spectra to a series of model simulated sferic spectra from using a finite difference time domain (FDTD) code.</p><p>In order to understand the time scales, magnitudes and sources for the midlatitude nighttime <italic>D</italic> region variations, we analyzed the sferic data of July and August 2005, and extracted both the height and sharpness of the <italic>D</italic> region electron density profile. The heights show large temporal variations of several kilometers on some nights and the relatively stable behavior on others. Statistical calculations indicate that the hourly average heights during the two months range between 82.0 km and 87.2 km with a mean value of 84.9 km and a standard deviation of 1.1 km. We also observed spatial variations of height as large as 2.0 km over 5 degrees latitudes on some nights, and no spatial variation on others. In addition, the measured height variations exhibited close correlations with local lightning occurrence rate on some nights but no correlation with local lightning or displaced lightning on others. The nighttime profile sharpness during 2.5 hours in two different nights was calculated, and the results were compared to the equivalent sharpness derived from International Reference Ionosphere (IRI) models. Both the absolute values and variation trends in IRI models are different from those in broadband measurements.</p><p>Based on sferic data similar to those for nighttime, we also measured the daytime <italic>D</italic> region electron density profile variations in July and August 2005 near Duke University. As expected, the solar radiation is the dominant but not the only determinant source for the daytime <italic>D</italic> region profile height temporal variations. The observed quiet time heights showed close correlations with solar zenith angle changes but unexpected spatial variations not linked to the solar zenith angle were also observed on some days, with 15% of days exhibiting regional differences larger than 0.5 km. During the solar flare, the induced height change was approximately proportional to the logarithm of the X-ray fluxes. During the rising and decaying phases of the solar flare, the height changes correlated more consistently with the short (wavelength 0.5-4 &Aring), rather than the long (wavelength 1-8 &Aring) X-ray flux changes. The daytime profile sharpness during morning, noontime and afternoon periods in three different days and for the solar zenith angle range 20 to 75 degrees was calculated. These broadband measured results were compared to narrowband VLF measurements, IRI models and Faraday rotation base IRI models (called FIRI). The estimated sharpness from all these sources was more consistent when the solar zenith angle was small than when it was large.</p><p>By applying the nighttime and daytime measurement techniques, we also derived the <italic>D</italic> region variations during sunrise and sunset periods. The measurements showed that both the electron density profile height and sharpness decrease during the sunrise period while increase during the sunset period.</p> / Dissertation