Search

Earth system tipping elements have been predominantly investigated for their potential to trigger \textit{negative} ecological, climatic, and societal shifts. Yet, an overlooked but seminal avenue exists in the form of \textit{positive} tipping phenomena, whose underlying mechanisms and benefits remain largely underexplored. To bridge this gap, our research introduces a fundamental percolation-based framework to assess the criticality and resilience of planetary terrestrial vegetation systems. Leveraging high-resolution satellite data, we focus on greening-induced positive tipping dynamics driven by global warming. We feature the Qinghai-Tibetan Plateau (QTP) and the Sahel region as contrasting yet analogous case studies. Our analysis uncovers an intriguing phenomenon where vegetation fragmentation aligns with a percolation threshold, exhibiting a scale-invariant pattern characterized by nearly perfect power laws with three critical exponents. Remarkably, contrary to conventional destructive tipping elements, these regions act as favorable tipping elements, transitioning from fragmented to cohesive vegetation patterns due to anthropogenic climate change and afforestation efforts. Furthermore, we propose an \textit{optimal resilience enhancement model} to reinforce vegetation robustness while minimizing socio-economic costs. This study provides valuable insights into the favorable aspects of tipping elements under climate change and offers effective strategies for enhancing ecological resilience against environmental threats.

The scientific process of earthquake forecasting involves estimating the probability and intensity of earthquakes in a specific area within a certain timeframe, based on seismic activity laws and observational data. Epidemic-Type Aftershock Sequence (ETAS) models, which rely on seismic empirical laws, is one of the most commonly used methods for earthquake forecasting. However, it underestimates feature in short-term time scale and overestimates in long-term time scale. Convolutional Long Short-Term Memory (ConvLSTM), has emerged as a promising approach capable of extracting spatio-temporal features. Hence, we propose a novel composite model named CL-ETAS model, which combines the strengths of ConvLSTM and ETAS model. We conduct experimental verification on real seismic data in Southern California. Our results show that CL-ETAS model outperforms both ETAS model and ConvLstm in accurately forecasting the number of earthquake events, their magnitude, and spatial distribution. Additionally, CL-ETAS model demonstrates notable enhancements in forecast stability and interpretability. These results offer valuable insights and recommendations for future earthquake forecasting endeavors.

Mainshocks are often followed by increased earthquake activity (aftershocks). According to the Omori-Utsu law, the rate of aftershocks decays as a power law over time. While aftershocks typically occur in the vicinity of the mainshock, previous studies have suggested that mainshocks can also trigger earthquakes in remote locations. Here we examine the earthquake rate in the days following mega-earthquakes (magnitude >= 7.5) and find that the rate is significantly lower beyond a certain distance from the epicenter compared to surrogate data. However, the remote earthquake rate after the strongest earthquakes (magnitude >= 8) can also be significantly higher than that of the rate based on surrogate data. Comparing our findings to the global ETAS model, we find that the model does not capture the earthquake rate found in the data, hinting at a potential missing mechanism. We suggest that the diminished earthquake rate is due the release of global energy/tension subsequent to substantial mainshock events. This conjecture holds the potential to enhance our comprehension of the intricacies governing post-seismic activity.

Here we focus on a basic statistical measure of earthquake catalogs that has not been studied before, the asymmetry of interevent time series (e.g., reflecting the tendency to have more aftershocks than spontaneous earthquakes). We define the asymmetry metric as the ratio between the number of positive interevent time increments minus negative increments and the total (positive plus negative) number of increments. Such asymmetry commonly exists in time series data for non-linear geophysical systems like river flow which decays slowly and increases rapidly. We find that earthquake interevent time series are significantly asymmetric, where the asymmetry function exhibits a significant crossover to weak asymmetry at large lag-index. We suggest that the Omori law can be associated with the large asymmetry at short time intervals below the crossover whereas overlapping aftershock sequences and the spontaneous events can be associated with a fast decay of asymmetry above the crossover. We show that the asymmetry is better reproduced by a recently modified ETAS model with two triggering processes in comparison to the standard ETAS model which only has one.

A prominent feature of earthquakes is their empirical laws including memory (clustering) in time and space. Several earthquake forecasting models, like the EpidemicType Aftershock Sequence (ETAS) model, were developed based on earthquake empirical laws. Yet, a recent study showed that the ETAS model fails in reproducing significant long-term memory characteristics found in real earthquake catalogs. Here we modify and generalize the ETAS model to include short- and long-term triggering mechanisms, to account for the short- and long-time memory (exponents) recently discovered in the data. Our generalized ETAS model reproduces accurately the short- and long-term/distance memory observed in the Italian and South California earthquake catalogs. The revised ETAS model is also found to significantly improve earthquake forecasting.

Over the past decades much effort has been devoted towards understanding and forecasting natural hazards. However, earthquake forecasting skill is still very limited and remains a great scientific challenge. The limited earthquake predictability is partly due to the erratic nature of earthquakes and partly to the lack of understanding the underlying mechanisms of earthquakes. To improve our understanding and potential forecasting, here we study the spatial and temporal long-term memory of interevent earthquakes above a certain magnitude using lagged conditional probabilities. We find, in real data, that the lagged conditional probabilities show long-term memory for both the interevent times and interevent distances and that the memory functions obey scaling and decay slowly with time, while, at a characteristic time, the decay crossesover to a fast decay. We also show that the ETAS model, which is often used to forecast earthquake events, yields scaling functions of the temporal and spatial interevent intervals which are not consistent with those of real data.