the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
An Ensemble Random Forest Model for Seismic Energy Forecast
Abstract. Seismic energy forecasting is critical for hazard preparedness, but current models have limits in accurately predicting seismic energy changes. This paper fills that gap by introducing a new ensemble random forest model designed specifically for seismic energy forecasting. Building on an existing paradigm, provided by Raghukanth et al. (2017), the global energy time series is decomposed into intrinsic mode functions (IMFs) using ensemble empirical mode decomposition for better representation. Following this approach, we split the data into stationary (IMF1) and non-stationary (sum of IMF2-IMF6) components for modeling. We acknowledge the inadequacy of intrinsic mode functions (IMFs) in capturing seismic energy dynamics, notably in anticipating the final values of the time series. To address this restriction, the yearly seismic energy time series is also fed along with the stationary and non-stationary parts as inputs to the developed models. Here, we employed Support Vector Machine (SVM), Random Forest (RF), Instance-Bases learning (IBk), Linear Regression (LR), and MultiLayer Perceptron (MLP) algorithms for the modelling. Furthermore, the five models discussed above were suitably employed in a novel regression-based ensemble random forest algorithm to arrive at the final predictions. The root mean square error (RMSE) obtained in the training and testing phases of the final model were 0.127 and 0.134, respectively. It was observed that the performance of the developed ensemble model was superior to those existing in literature (Raghukanth et al., 2017). Further, the developed algorithm was employed for the seismic energy prediction in the active Western Himalayan region for a comprehensively compiled catalogue and the mean forecasted seismic energy for year 2024 is 7.21 × 1014 J. This work is a pilot project that aims to create a forecast model for the release of seismic energy globally and further application at a regional level. The findings of our investigation demonstrate the possibility of the established method in the accurate seismological energy forecast, which can help with appropriate hazard preparedness.
- Preprint
(3637 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (until 26 Dec 2024)
-
RC1: 'Comment on nhess-2024-129', Anonymous Referee #1, 13 Nov 2024
reply
The manuscript under consideration is focused on an open problem, the prediction of large earthquakes. This paper presents an interesting study based on Machine Learning, ML, to predict intense earthquakes in terms of the released energy.
In strictly sense, prediction is considered when the position of the epicenter, the time of occurrence and the magnitude are known with precision.
As it is well known, the underlying dynamics of tectonic plates is due to complex processes inside the Earth such as convection. Those involved processes give rise to the interaction between tectonic plates, being stick-slip the main mechanism for the earthquakes occurrence. On the basis of these aspects, seismic dynamics is complex.
This paper shows a good prediction approximation produced by ML-based algorithms.
The work is very important and their results are very interesting, however, in my opinion, some important aspects have not been considered in this study:
a) The catalogue use for the study considers a period from 1900 to 2015. In the period from 1900 to 1920 very few earthquakes are observed while in recent years the density of events is much higher. Authors should explain these differences time periods.
b) Figure 1a shows the magnitude of completeness (M6.4) which is only part of the Gutenberg-Richter law. The authors should show the Gutenberg-Richter law taking lower magnitudes (for example, M < 3) and thus determine the correct magnitude of completeness. To do this they can use the ZMAP platform where they can estimate the b-value and completeness Mc of the GR law.
c) It is important that the authors justify the reason for considering only earthquakes with large magnitudes and avoid the other ones with low magnitudes.
d) To calculate the IMF functions authors used the algorithm proposed by Huang et al. (1998). In my opinion, the authors should explain how the complex dynamics of seismic activity is involved to obtain IMF to obtain the smooth curve (a) in Figure 2.
e) When applying the algorithm proposed by Huang et al. it is not clear the role seismic activity with magnitudes M < Mc could have.
Their results do not consider time and epicentral position.
I consider that these suggestions can improve the context of the study that the authors have developed and then could be published
Citation: https://doi.org/10.5194/nhess-2024-129-RC1 -
RC2: 'Reply on RC1', Anonymous Referee #1, 21 Nov 2024
reply
I have read the author's reply. Their comments are very interesting and improve widely their paper. I don't have any other comment and, for me the new version of their manuscript can be published in the currently version.
I appreciate the answers to my suggestions from the authors.
Citation: https://doi.org/10.5194/nhess-2024-129-RC2 -
AC2: 'Reply on RC2', Sukh Sagar Shukla, 27 Nov 2024
reply
Dear Anonymous Referee #1,
Thank you very much for considering this manuscript for publication, and many thanks for the positive comments and reviews.
Citation: https://doi.org/10.5194/nhess-2024-129-AC2 -
RC3: 'Reply on AC2', Anonymous Referee #1, 27 Nov 2024
reply
Citation: https://doi.org/
10.5194/nhess-2024-129-RC3
-
RC3: 'Reply on AC2', Anonymous Referee #1, 27 Nov 2024
reply
-
AC2: 'Reply on RC2', Sukh Sagar Shukla, 27 Nov 2024
reply
-
RC2: 'Reply on RC1', Anonymous Referee #1, 21 Nov 2024
reply
-
AC1: 'Reply on RC1', Sukh Sagar Shukla, 20 Nov 2024
reply
Dear Anonymous Referee #1,
The authors would like to thank the reviewer for recognising the relevance of the present work, and finding the results interesting. We highly appreciate the insightful comments, and we are committed to making significant revisions. We have carefully considered all the suggestions and made substantial revisions to the manuscript accordingly. Please find the attached PDF file containing a documented list of changes we have made to the manuscript (marked as reply). We have considered a more comprehensive and updated global catalog for analysis and trained the models on updated dataset. We have refined the discussion on considering large magnitude events and physical interpretation of IMFs. We hope these clarifications will improve the readers understanding of our work.
Kind Regards,
Sukh Sagar Shukla
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
156 | 49 | 15 | 220 | 5 | 6 |
- HTML: 156
- PDF: 49
- XML: 15
- Total: 220
- BibTeX: 5
- EndNote: 6
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1