Identification of earthquake precursors in soil radon-222 data of Kutch, Gujarat, India using empirical mode decomposition based Hilbert Huang Transform

Abstract

Soil radon (Rn-222) has been continuously monitored at Badargadh station (23.47°N, 70.62°E) in Kutch region of Gujarat to study the pre-seismic anomalies prior to occurrence of local earthquakes. This monitoring site is in close proximity to the South Wagad Fault, a seismically active fault in the study area. The raw data of radon along with meteorological parameters such as temperature, pressure and humidity in soil of this station for the period of January 01 to December 31, 2017 with a sampling interval of 10 min were used in the analysis. The wind speed and rainfall data of the corresponding period were collected from the nearest weather station. From descriptive statistics, we found an average soil radon concentration of 343 Bq.m−3. It is observed that radon has a maximum concentration during the rainy season compared to the other two seasons. We found that radon emission rate is less during mid-nights and early morning, whereas, the radon emission is more during afternoon hours when the sun light intensity is more. In order to identify and extract the periodic oscillations in the radon time series, the Empirical Mode Decomposition (EMD) was applied to the soil radon (Rn-222) time series by decomposing it into different oscillatory modes known as the Intrinsic Mode Function (IMF). Several interesting non-linear features emerged from the analysis after applying Hilbert Huang Transform (HHT) on significant IMFs. The temporal variation of the instantaneous energy is well correlated with four local earthquakes during the study period. Most interestingly, intermittencies in the temporal evolution of the instantaneous energy function have been observed prior to these local earthquakes. We present the results of the seismic and aseismic periods as well as a brief discussion of the analysis of radon data which can be used as a precursor of seismic activity. It is now possible to identify anomalies in radon time series using EMD based HHT method even for small-magnitude earthquakes.