Search

We employ the eigen microstate approach to explore the self-organized criticality (SOC) in two celebrated sandpile models, namely, the BTW model and the Manna model. In both models, phase transitions from the absorbing-state to the critical state can be understood by the emergence of dominant eigen microstates with significantly increased weights. Spatial eigen microstates of avalanches can be uniformly characterized by a linear system size rescaling. The first temporal eigen microstates reveal scaling relations in both models. Furthermore, by finite-size scaling analysis of the first eigen microstate, we numerically estimate critical exponents i.e., $\sqrt{\sigma_0 w_1}/\tilde{v}_{1} \propto L^D$ and $\tilde{v}_{1} \propto L^{D(1-\tau_s)/2}$. Our findings could provide profound insights into eigen states of the universality and phase transition in non-equilibrium complex systems governed by self-organized criticality.

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.

Studies on stratospheric ozone have attracted much attention due to its serious impacts on climate changes and its important role as a tracer of Earth's global circulation. Tropospheric ozone as a main atmospheric pollutant damages human health as well as the growth of vegetation. Yet there is still a lack of a theoretical framework to fully describe the variation of ozone. To understand ozone's spatiotemporal variance, we introduce the eigen microstate method to analyze the global ozone mass mixing ratio (OMMR) between 1979-01-01 and 2020-06-30 at 37 pressure layers. We find that eigen microstates at different geopotential heights can capture different climate phenomena and modes. Without deseasonalization, the first eigen microstates capture the seasonal effect and reveal that the phase of the intra-annual cycle moves with the geopotential heights. After deseasonalization, by contrast, the collective patterns from the overall trend, ENSO, QBO, and tropopause pressure are identified by the first few significant eigen microstates. The theoretical framework proposed here can also be applied to other complex Earth systems., Comment: submitted to Chaos, and already accepted!

The El Nino-Southern Oscillation (ENSO) is one of the most important phenomena in climate. By studying the fluctuations of surface air temperature within one year between 1979-01-01 and 2016-12-31 of the region (30S-30N, 0E-360E) with eigen-decomposition, we find that the temperature fluctuations are dominated by the two principal modes whose temporal evolutions respond significantly to the ENSO variability. According to introduce a micro-correlation, we find that the coupling between the first principal mode and the temperature fluctuations in the El Nino region could result in different ENSO phases. Without this coupling, the El Nino region is in a normal phase. With the strong coupling between the El Nino region and the Northern Hemisphere, an El Nino event will appear with a high probability. Then this coupling changes to be strong between the El Nino region and the Southern Hemisphere accounting for the fast decay of El Nino after boreal winter, even leading to a La Nina event. Moreover, the coupling between the El Nino region and the second principal mode is found to be strong in normal or La Nina phases in response to the normal or strong Walker Circulations. We conjecture that the temporal evolutions of these couplings for the first and second principal modes are controlled by the meridional and zonal ocean-atmospheric circulations respectively.

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.

Air pollution is associated with human diseases and has been found to be related to premature mortality. In response, environmental policies have been adopted in many countries, to decrease anthropogenic air pollution for the improvement of long-term air quality, since most air pollutant sources are anthropogenic. However, air pollution fluctuations have been found to strongly depend on the weather dynamics. This raises a fundamental question: What are the significant atmospheric processes that affect the local daily variability of air pollution? For this purpose, we develop here a multi-layered network analysis to detect the interlinks between the geopotential height of the upper air (~5 km) and surface air pollution in both China and the USA. We find that Rossby waves significantly affect air pollution fluctuations through the development of cyclone and anticyclone systems, and further affect the local stability of the air and the winds. The significant impacts of Rossby waves on air pollution are found to underlie most of the daily fluctuations in air pollution. Thus, the impact of Rossby waves on human life is greater than previously assumed. The rapid warming of the Arctic could slow down Rossby waves, thus increasing human health risks. Our method can help to determine the risk assessment of such extreme events and can improve potential predictability.

We propose the finite-size scaling of correlation function in a finite system near its critical point. At a distance ${\bf r}$ in the finite system with size $L$, the correlation function can be written as the product of $|{\bf r}|^{-(d-2+\eta)}$ and its finite-size scaling function of variables ${\bf r}/L$ and $tL^{1/\nu}$, where $t=(T-T_c)/T_c$. The directional dependence of correlation function is nonnegligible only when $|{\bf r}|$ becomes compariable with $L$. This finite-size scaling of correlation function has been confirmed by correlation functions of the Ising model and the bond percolation in two-diemnional lattices, which are calculated by Monte Carlo simulation. We can use the finite-size scaling of correlation function to determine the critical point and the critical exponent $\eta$., Comment: 7 pages, 13 figures

Air pollution has become a major issue and caused widespread environmental and health problems. Aerosols or particulate matters are an important component of the atmosphere and can transport under complex meteorological conditions. Based on the data of $PM_{2.5}$ observations, we develop a network approach to study and quantify their spreading and diffusion patterns. We calculate cross-correlation functions of time lag between sites within different season. The probability distribution of correlation changes with season. It is found that the probability distributions in four seasons can be scaled into one scaling function with averages and standard deviations of correlation. This seasonal scaling behavior indicates there is the same mechanism behind correlations of $PM_{2.5}$ concentration in different seasons. Further, from weighted and directional degrees of complex network, different properties of $PM_{2.5}$ concentration are studied. The weighted degrees reveal the strongest correlations of $PM_{2.5}$ concentration in winter and in the North China plain. These directional degrees show net influences of $PM_{2.5}$ along Gobi and inner Mongolia, the North China plain, Central China, and Yangtze River Delta., Comment: 7 pages, 11 figures

For vertical velocity field $v_{\rm z} (r,z;R)$ of granular flow through an aperture of radius $R$, we propose a size scaling form $v_{\rm z}(r,z;R)=v_{\rm z} (0,0;R)f (r/R_{\rm r}, z/R_{\rm z})$ in the region above the aperture. The length scales $R_{\rm r}=R- 0.5 d$ and $R_{\rm z}=R+k_2 d$, where $k_2$ is a parameter to be determined and $d$ is the diameter of granule. The effective acceleration, which is derived from $v_{\rm z}$, follows also a size scaling form $a_{\rm eff} = v_{\rm z}^2(0,0;R)R_{\rm z}^{-1} \theta (r/R_{\rm r}, z/R_{\rm z})$. For granular flow under gravity $g$, there is a boundary condition $a_{\rm eff} (0,0;R)=-g$ which gives rise to $v_{\rm z} (0,0;R)= \sqrt{ \lambda g R_{\rm z}}$ with $\lambda=-1/\theta (0,0)$. Using the size scaling form of vertical velocity field and its boundary condition, we can obtain the flow rate $W =C_2 \rho \sqrt{g } R_{\rm r}^{D-1} R_{\rm z}^{1/2} $, which agrees with the Beverloo law when $R \gg d$. The vertical velocity fields $v_z (r,z;R)$ in three-dimensional (3D) and two-dimensional (2D) hoppers have been simulated using the discrete element method (DEM) and GPU program. Simulation data confirm the size scaling form of $v_{\rm z} (r,z;R)$ and the $R$-dependence of $v_{\rm z} (0,0;R)$., Comment: 5 pages, 8 figures