Smarter Features, Not Bigger Models, Crack Earthquake Forecasting in Central Asia

Earthquakes are among the most stubborn prediction problems in all of science, and a new study from Kazakhstan suggests that the path forward may lie less in exotic neural architectures and more in how the problem itself is framed. Writing in the Journal of Big Data, a team led by Marat Nurtas of the Ionosphere Institute and the International Information Technology University in Almaty reports a macro-scale earthquake forecasting framework for Central Asia that achieves a roughly tenfold improvement over a naive statistical baseline, reaching a Precision-Recall Area Under the Curve of approximately 0.451 and a Receiver Operating Characteristic Area Under the Curve of about 0.844. The striking twist is that six fundamentally different machine learning architectures, from gradient boosting to recurrent neural networks, all converged on nearly identical performance, pointing to a fundamental predictability ceiling rather than a modeling shortfall.

The research tackles a problem that has long plagued computational seismology: extreme class imbalance. When earthquake forecasting is cast as a grid-based classification task, the region under study is divided into spatial cells, and the model must predict whether a seismic event will occur in each cell during each time window. For Central Asia, the team discretized the territory into one-degree by one-degree grid cells and attempted to forecast earthquakes of magnitude 3.0 or greater on a weekly basis. Under this formulation, roughly 96 percent of all cell-week combinations contain no event at all, a phenomenon known as zero inflation. The true event prevalence sits at only about 4.5 percent, which means that a model doing nothing more than predicting ‘no earthquake’ everywhere would still appear superficially accurate while being scientifically useless.

This imbalance has profound consequences for how forecasting models must be evaluated. Standard accuracy metrics become meaningless when negative cases dominate by more than twenty to one. The researchers therefore anchored their evaluation in the Precision-Recall Area Under the Curve, a metric that is far more sensitive to performance on the rare positive class. With a prevalence of 4.5 percent, the constant baseline for PR-AUC is 0.045, meaning any model must substantially exceed that value to demonstrate genuine predictive skill. The achieved score of 0.451 represents a tenfold improvement over this baseline, a substantial margin in a domain where even modest gains above chance are considered meaningful by the seismological community.

Central to the study is a carefully engineered feature space grounded in earthquake physics rather than raw statistical patterns. The framework integrates tectonic regime-conditioned normalization, which allows the model to account for the fact that different tectonic settings produce fundamentally different seismic behavior, so that features extracted from a thrust-fault environment are not treated as directly comparable to those from a strike-slip regime. It also incorporates Omori energy decay proxies, mathematical representations of the well-documented tendency of earthquake sequences to produce aftershocks whose frequency decays over time following a mainshock. These proxies give the models a physically interpretable signal about the temporal clustering of seismicity, encoding decades of seismological understanding directly into the input data.

Structural fault descriptors form a third pillar of the feature design. The geometry, orientation, and proximity of mapped fault systems are among the strongest known controls on where earthquakes occur, and by encoding these structural characteristics as model inputs, the framework ensures that the learning algorithms operate on geologically meaningful quantities rather than arbitrary grid statistics. The final and perhaps most consequential innovation is log-odds baseline initialization, a technique that encodes the historical cell-specific event rate directly into the learning objective. Instead of forcing each model to rediscover from scratch the simple fact that some grid cells are historically far more seismically active than others, the initialization embeds this prior knowledge into the model’s starting point, allowing learning effort to focus on deviations from the historical pattern.

To determine whether performance under such extreme imbalance is governed primarily by model architecture or by structured feature design, the researchers evaluated six heterogeneous architectures under a strict chronological split, ensuring that models were trained only on past data and tested on future periods, exactly as an operational forecasting system would be deployed. The architectures spanned a wide methodological range, including CatBoost and other gradient boosting methods, which excel at tabular data, and Bi-LSTM networks, a bidirectional long short-term memory architecture capable of capturing temporal dependencies in sequential data. Despite their radically different inductive biases and internal mechanics, the models converged on nearly identical PR-AUC and ROC-AUC values, a result the authors interpret as evidence that the information content of the feature space, not the capacity of the learner, is the binding constraint.

Equally notable is what the framework does not do. Many studies confronting severe class imbalance resort to synthetic resampling techniques, such as oversampling the rare event class or undersampling the dominant negative class, to artificially balance the training distribution. These methods can distort the learned probability calibration, producing models whose confidence scores no longer correspond to real-world event likelihoods. The Central Asia framework achieves its tenfold improvement entirely without synthetic resampling, preserving the integrity of the probability estimates. This matters enormously for practical applications, because emergency management authorities require calibrated forecasts whose stated probabilities can be trusted when weighing evacuation decisions, infrastructure inspections, and public warnings.

The authors argue that their findings point to the existence of a macro-scale predictability ceiling in seismic forecasting. If architecturally diverse models, given the same physically structured inputs, all plateau at the same performance level, the implication is that the remaining unpredictability reflects genuine stochasticity in the earthquake process at this spatial and temporal resolution, rather than a deficiency of current algorithms. This interpretation carries a sobering but valuable message for the field: further architectural innovation alone is unlikely to break through the ceiling, while improvements in physical understanding, richer observational data streams, and better-calibrated baselines may still push the boundary outward. It also cautions against the common practice of claiming architectural superiority from small performance differences that may fall within the noise of a shared predictability limit.

The work was funded by the Committee of Science of the Ministry of Science and Higher Education of the Republic of Kazakhstan under a grant for developing a multifunctional system of ground-space monitoring and early warning of natural and technogenic emergencies, underscoring its operational motivation. For a country situated in one of the most seismically active zones of Central Asia, where the collision of the Indian and Eurasian plates drives hazardous tectonics through the Tien Shan and surrounding mountain belts, reliable macro-scale forecasting is not an academic curiosity but a matter of public safety. By demonstrating that disciplined feature engineering, physically informed priors, and rigorous baseline calibration can deliver a tenfold gain in predictive skill without exotic machinery, the Almaty team has provided both a practical forecasting tool and a methodological lesson that resonates far beyond seismology: in data-starved, imbalance-dominated problems, how you frame the question often matters more than how elaborate your model is.

Subject of Research: Machine learning earthquake forecasting under extreme class imbalance in Central Asia

Article Title: Macro-scale earthquake forecasting under class imbalance in Central Asia

Article References: Nurtas, M., Nurakynov, S., Sakabekov, A., Altaibek, A., Kumarkhanova, A., & Merekeyev, A. (2026). Macro-scale earthquake forecasting under class imbalance in Central Asia. Journal of Big Data. https://doi.org/10.1186/s40537-026-01544-z

Image Credits: AI Generated

DOI: 10.1186/s40537-026-01544-z

Keywords: earthquake forecasting, Central Asia, class imbalance, machine learning, CatBoost, gradient boosting, Bi-LSTM, PR-AUC, spatio-temporal prediction, Omori decay, fault descriptors, Kazakhstan