No Arabic abstract
Natural landslides exhibit scaling properties revealed by power law relationships. These relationships include the frequency of the size (e.g., area, volume) of the landslides, and the rainfall conditions responsible for slope failures in a region. Reasons for the scaling behavior of landslides are poorly known. We investigate the possibility of using the Transient Rainfall Infiltration and Grid-Based Regional Slope-Stability analysis code (TRIGRS), a consolidated, physically-based, numerical model that describes the stability/instability conditions of natural slopes forced by rainfall, to determine the frequency statistics of the area of the unstable slopes and the rainfall intensity (I) - duration (D) conditions that result in landslides in a region. We apply TRIGRS in a portion of the Upper Tiber River Basin, Central Italy. The spatially distributed model predicts the stability/instability conditions of individual grid cells, given the local terrain and rainfall conditions. We run TRIGRS using multiple, synthetic rainfall histories, and we compare the modeling results with empirical evidences of the area of landslides and of the rainfall conditions that have caused landslides in the study area. Our findings revealed that TRIGRS is capable of reproducing the frequency of the size of the patches of terrain predicted as unstable by the model, which match the frequency size statistics of landslides in the study area, and the mean rainfall D, I conditions that result in unstable slopes in the study area, which match rainfall I-D thresholds for possible landslide occurrence. Our results are a step towards understanding the mechanisms that give rise to landslide scaling properties.
Distributed models to forecast the spatial and temporal occurrence of rainfall-induced shallow landslides are based on deterministic laws. These models extend spatially the static stability models adopted in geotechnical engineering, and adopt an infinite-slope geometry to balance the resisting and the driving forces acting on the sliding mass. An infiltration model is used to determine how rainfall changes pore-water conditions, modulating the local stability/instability conditions. A problem with the operation of the existing models lays in the difficulty in obtaining accurate values for the several variables that describe the material properties of the slopes. The problem is particularly severe when the models are applied over large areas, for which sufficient information on the geotechnical and hydrological conditions of the slopes is not generally available. To help solve the problem, we propose a probabilistic Monte Carlo approach to the distributed modeling of rainfall-induced shallow landslides. For the purpose, we have modified the Transient Rainfall Infiltration and Grid-Based Regional Slope-Stability Analysis (TRIGRS) code. The new code (TRIGRS-P) adopts a probabilistic approach to compute, on a cell-by-cell basis, transient pore-pressure changes and related changes in the factor of safety due to rainfall infiltration. Infiltration is modeled using analytical solutions of partial differential equations describing one-dimensional vertical flow in isotropic, homogeneous materials. Both saturated and unsaturated soil conditions can be considered. TRIGRS-P copes with the natural variability inherent to the mechanical and hydrological properties of the slope materials by allowing values of the TRIGRS model input parameters to be sampled randomly from a given probability distribution. [..]
The tendency of irreversible processes to generate entropy is the ultimate driving force for the evolution of nature. In engineering, entropy production is often used as a measure of usable energy losses. In this study we show that the analysis of the entropy production patterns can help understand the vastly diversified experimental observations of water-rock interactions in natural porous media. We first present a numerical scheme for the analysis of entropy production in dissolving porous media. Our scheme uses a greyscale digital model of natural chalk obtained by X-ray nanotomography. Greyscale models preserve structural heterogeneities with very high fidelity, which is essential for simulating a system dominated by infiltration instability. We focus on the coupling between two types of entropy production: the percolative entropy generated by dissipating the kinetic energy of fluid flow and the reactive entropy that originates from the consumption of chemical free energy. Their temporal patterns pinpoint three stages of microstructural evolution. We then show that the regional mixing deteriorates infiltration instability by reducing local variations in reactant distribution. In addition, we show that the microstructural evolution can be particularly sensitive to the initially present transport heterogeneities when the global flowrate is small. This dependence on flowrate indicates that the need to resolve the structural features of a porous system is greater when the residence time of the fluid is long.
Neural avalanches are collective firings of neurons that exhibit emergent scale-free behavior. Understanding the nature and distribution of these avalanches is an important element in understanding how the brain functions. We study a model of neural avalanches for which the dynamics are governed by neutral theory. The neural avalanches are defined using causal connections between the firing neurons. We analyze the scaling of causal neural avalanches as the critical point is approached from the absorbing phase. By using cluster analysis tools from percolation theory, we characterize the critical properties of the neural avalanches. We identify the tuning parameters consistent with experiments. The scaling hypothesis provides a unified explanation of the power laws which characterize the critical point. The critical exponents characterizing the avalanche distributions and divergence of the response functions are consistent with the predictions of the scaling hypothesis. We use a universal scaling function for the avalanche profile to find that the firing rates for avalanches of different durations show data collapse after appropriate rescaling. We also find data collapse for the avalanche distribution functions, which is stronger evidence of criticality than just the existence of power laws. Critical slowing-down and power law relaxation of avalanches is observed as the system is tuned to its critical point. We discuss how our results motivate future empirical studies of criticality in the brain.
Production in an economy is a set of firms activities as suppliers and customers; a firm buys goods from other firms, puts value added and sells products to others in a giant network of production. Empirical study is lacking despite the fact that the structure of the production network is important to understand and make models for many aspects of dynamics in economy. We study a nation-wide production network comprising a million firms and millions of supplier-customer links by using recent statistical methods developed in physics. We show in the empirical analysis scale-free degree distribution, disassortativity, correlation of degree to firm-size, and community structure having sectoral and regional modules. Since suppliers usually provide credit to their customers, who supply it to theirs in turn, each link is actually a creditor-debtor relationship. We also study chains of failures or bankruptcies that take place along those links in the network, and corresponding avalanche-size distribution.
Observations of tropical convection from precipitation radar and the concurring large-scale atmospheric state at two locations (Darwin and Kwajalein) are used to establish effective stochastic models to parameterise subgrid-scale tropical convective activity. Two approaches are presented which rely on the assumption that tropical convection induces a stationary equilibrium distribution. In the first approach we parameterise convection variables such as convective area fraction as an instantaneous random realisation conditioned on the large-scale vertical velocities according to a probability density function estimated from the observations. In the second approach convection variables are generated in a Markov process conditioned on the large-scale vertical velocity, allowing for non-trivial temporal correlations. Despite the different prevalent atmospheric and oceanic regimes at the two locations, with Kwajalein being exposed to a purely oceanic weather regime and Darwin exhibiting land-sea interaction, we establish that the empirical measure for the convective variables conditioned on large-scale mid-level vertical velocities for the two locations are close. This allows us to train the stochastic models at one location and then generate time series of convective activity at the other location. The proposed stochastic subgrid-scale models adequately reproduce the statistics of the observed convective variables and we discuss how they may be used in future scale-independent mass-flux convection parameterisations.