We investigated, for the first time, long-term (1976-2010) variations in reflection heights of tweek atmospherics based on very low frequency (VLF) observations at Kagoshima, Japan. The results revealed the effects of the solar cycle on the nighttime lower ionosphere at low to middle latitudes. The tweek reflection heights on geomagnetically quiet days were analyzed every month over three solar cycles by using an automated spectral fitting procedure to estimate the cutoff frequency. The average and standard deviation of the reflection height were 95.9 km and +/-3.1 km, respectively. Typical periods of time variation for the reflection height were 13.3, 3.2, 1.3, 1.0, 0.6, and 0.5 years. The variations in tweek reflection heights did not show simple anticorrelation with solar activity. The correlation coefficient between tweek reflection height and sunspot number was 0.03 throughout the three solar cycles. Hilbert-Huang transform analysis successfully indicated the presence of 0.5-1.5 year and similar to 10 year variations as intrinsic mode functions (IMF). The decomposed IMF with the similar to 10 year variation had a positive correlation with sunspot numbers and a negative correlation with galactic cosmic rays (GCRs). We hypothesize that these variations in tweek reflection heights could be caused by coupling of several ionization effects at the D and lower E regions, effects such as geocorona, GCRs, particle precipitation, and variations in neutral density in the lower thermosphere. Among these processes, the geocorona and particle precipitation could show negative correlation, while the GCRs and neutral density could show positive correlation with solar activities.