Earthquake events are usually associated with atmospheric processes variations. Some of the important events in seismic periods include the atmospheric electricity/conductivity modification, which mostly depends on the ion population in the vicinity of earthquake preparation area. Radon, together with cosmic radiation are the major ionization sources in the lower troposphere. The radon induced soil ion-pair production rate was estimated for Erzincan, a city along the North Anatolian Fault Zone, Turkey. The Long Short-Term Memory networks (LSTM's) and Monte Carlo method are proposed to account for uncertainty in estimating the ion production rate. The advantage of the LSTM model is taken to study radon anomalies during the M = 5 Girlevik earthquake, Erzincan. Radon concentration is found to increase prior to the 5.0 Girlevik earthquake. According to estimations, 23 × 109 ion-pairs m−3s−1 were generated at 1 m depth from the Earth's surface during this earthquake. The soil-ion production estimation rate due to radon and its progeny in Erzincan is at the order of 109 ion-pairs m−3s−1. When these generated ions and charged aerosols reach the surface and were added to other ionization sources, parameters such as the near surface conductivity and electric field would be modified. In addition, under pronounced conditions could favor ionospheric perturbations.