Dynamic Soil Research Articles

A non-linear stress–strain relationship for soil and rock

While many employ a hyperbolic stress–strain relationship for soils, it is known that such a relationship is accurate over the small-strain range as encountered in earthquake and soil dynamics problems. For large strains, a relationship with different input parameters is needed, as is required for finite-element analyses of large-deformation behaviour. These separate characterisations do not carry over into the other application. Here, a proposed power relationship is presented that was developed to characterise the triaxial test stress–strain behaviour of cohesionless material from lubricated or ‘frictionless’ cap and base tests (144 tests) covering a range in the natural variation in particle size, particle shape and surface roughness, over low to high confining pressure. This relationship covers the strain range from 10−6 to soil failure. To date, it has been successfully used in laterally loaded pile response characterisation (the strain wedge model) and shallow-foundation load–settlement–bearing capacity response. Most recently, it has been extended to assess the behaviour of rock-like material (caliche). The relationship and its comparison with the hyperbolic relationship for large strain and the shear modulus reduction curve for seismic behaviour are presented.

Using Relationships Between Vegetation and Surface Soil Biogeochemical Properties to Assess Regional Soil Carbon Inventories for South Baffin Island, Nunavut, Canada

AbstractAs Arctic regions warm rapidly, it is unclear whether high‐latitude soil carbon (C) will decrease or increase. Predicting future dynamics of Arctic soil C stocks requires a better understanding of the quantities and controls of soil C. We explore the relationship between vegetation and surface soil C in an understudied region of the Arctic: Baffin Island, Nunavut, Canada. We combined soil C data for three vegetation types—polar desert, mesic tundra, and wet meadow—with a vegetation classification to upscale soil C stocks. Surface soil C differed significantly across vegetation types, and interactions existed between vegetation type and soil depth. Polar desert soils were consistently mineral, with relatively thin organic layers, low percent C, and high bulk density. Mesic soils exhibited an organic‐rich epipedon overlying mineral soil. Wet meadows were consistently organic soil with low bulk density and high percent C. For the top 20 cm, polar desert contained the least soil C (2.17 ± 0.48 kg m−2); mesic tundra had intermediate C (8.92 ± 0.74 kg m−2); wet meadow stored the most C (13.07 ± 0.69 kg m−2). Extrapolating to the top 30 cm, our results suggest that approximately 44 Tg C is stored in the study region with a mean landscape soil C stock of 2.75 kg m−2 for non‐water areas. Combining vegetation mapping with local soil C stocks considerably narrows the range of estimates from other upscaling approaches (27–189 Tg) for soil C on South Baffin Island.

Assessment of Liquefaction Potential by Comparing Semi-Empirical Methods Based on the CPT Test

Abstract Liquefaction is a dangerous and temporary phenomenon whereby water-saturated soil loses all or part of its strength. Undrained conditions associated with cyclic loading increase water pressure in soil pores, thereby reducing effective stress. The aim of this study is, on the one hand, to report on the phenomenon of liquefaction of sand, clay and silt deposits in more or less water-saturated zones in the section located at the heart of the central alluvial plain of the Oued Sebou in the mio-plioquaternary Gharb basin and, on the other hand, to study the ability of semi-empirical methods to correctly assess liquefaction potential, while specifying the most appropriate method for the area studied. The study is based on data from experimental results of static penetrometer tests between the Mnasra - Ouelad Salama zone in the Oued Sebou alluvial plain of Morocco's mio-plio-quaternary Gharb basin, made up of sandy, sandy-clay, sandy-silt and silty-sandy formations, which are more sensitive to liquefaction due to their saturation and grain size. We present and discuss the results of Olsen's method, Juang's method and Robertson's method, which are based on the CPT static penetrometer test, as well as looking at the impact of dynamic loading and soil structure on liquefaction probability index values.

Differential Responses of Soil Microbial and Carbon‐Nitrogen Processes to Future Environmental Changes Across Soil Depths and Environmental Factors

AbstractAccurately predicting carbon‐climate feedbacks relies on understanding the environmental factors regulating soil organic carbon (SOC) storage and dynamics. Here, we employed a microbial ecological model (MEND), driven by downscaled output data from six Earth system models under two Shared Socio‐economic Pathways (SSP1‐2.6 and SSP5‐8.5) scenarios, to simulate long‐term soil biogeochemical processes. We aim to analyze the responses of soil microbial and carbon‐nitrogen (C‐N) processes to changes in environmental factors, including litter input (L), soil moisture (W) and temperature (T), and soil pH, in a broadleaf forest (BF) and a pine forest (PF). For the entire soil layer in both forests, we found that, compared to the baseline period of 2009–2020, the mean SOC during 2081–2100 increased by 40.9%–90.6% under the L or T change scenarios, versus 5.2%–31.0% under the W change scenario. However, soil moisture emerged as a key regulator of SOC, MBC and inorganic N dynamics in the topsoil of BF and PF. For example, W change led to SOC gain of 5.5%–37.2%, compared to the SOC loss of 15.5%–18.0% under L or T scenario. Additionally, a further reduction in soil pH by 0.2 units in the BF, representing the acid rain effect, significantly resulted in an additional SOC gain by 14.2%–21.3%, compared to the LTW (simultaneous changes in the three factors) scenario. These results indicate that the results derived solely from topsoil may not be extrapolated to the entire soil profile. Overall, this study significantly advances our comprehension of how different environmental factors impact the dynamics of SOC and the implications they have for climate change.

Space–time finite element method with domain reduction techniques for dynamic soil–structure interaction problems

AbstractDesign of earth structures, such as dams, tunnels, and embankments, against the vibrational loading caused by high‐speed trains, road traffic, underground explosions, and, more importantly, earthquake motion, demands solutions of the dynamic soil–structure Interaction (SSI) problem. This paper presents a velocity‐based space–time finite element procedure, v‐ST/finite element method (FEM), to solve dynamic SSI problems. The main goal of this study is to present the computation details of implementing viscous boundary conditions of Lysmer–Kuhlemeyer to truncate the unbounded soil domain. Furthermore, additional time‐dependent boundary conditions, in terms of the free‐field response, are included to facilitate energy flow from the far field to the computation domain at the vertical truncated boundaries. In the FEM, seismic input motion is applied to an effective nodal force vector, which can be obtained explicitly in the numerical simulations. Finally, the response of a concrete gravity dam resting on an elastic half‐space to the horizontal component of earthquake motion is computed and successfully compared with the results of semidiscrete FEM using the Newmark‐ method.

Collective Identification of Superstructures and Dynamic Soil Springs Considering the Effects of Higher Modes Based on the MIEC Method

Collective Identification of Superstructures and Dynamic Soil Springs Considering the Effects of Higher Modes Based on the MIEC Method

Cone penetration tests and dynamic soil properties

Cone penetration tests and dynamic soil properties

Indicative microzonation map of the contact zone between the kenny hill and limestone formation using active, passive and hybrid MASW methods

Abstract The impact of the far-field earthquake generated by the Sumatra-Andaman earthquake (Richter scale of 9.3) was also felt by Malaysians. The ground motion triggered by this earthquake is caused by the shear wave velocity (Vs), particularly at a depth of 30 m (Vs30). Therefore, this dynamic soil properties need to be determined, which was evaluated using the active, passive, and hybrid multichannel analysis of surface waves (MASW) methods at 16 locations in the contact zone between the Kenny Hill and limestone formations in Kuala Lumpur, Malaysia. The Vs determined from the MASW for the contact zone ranging from 198 m/s to 286 m/s at the ground surface and from 376 m/s to 662 m/s at 30 m depth. In contrast, in the non-contact zone, Vs range from 212 m/s to 328 m/s at the ground surface and from 562 m/s to 921 m/s at 30 m depth. For Vs30, the contact zone has a minimum value of 280 m/s and a maximum value of 380 m/s, while the non-contact zone has a minimum value of 370 m/s and a maximum value of 530 m/s.

Influence of Morphometric Relief Parameters on Soil Depth Changes and Humus Horizon Thickness in Relation to Erosion-Accumulation Processes: A Study in the Ipeľská Pahorkatina Hills, Slovakia

Abstract This study examines the spatial distribution of soil types and their susceptibility to erosion and accumulation processes in a study area in Slovakia. Field research involving 71 probes identified various soil types, with Regosols and Cutanic Luvisols being predominant. The study found that erosion-accumulation processes were detected in 69.97% of the probes, with changes observed in soil horizons. Soil analysis revealed different relations between soil depth, humus thickness, and terrain characteristics such as slope, slope length, and slope length and steepness factor (LS factor). Specifically, we confirmed a moderately strong positive correlation between soil depth and humus thickness (r = 0.597, n = 71, p < 0.001). Shallow soils (0–30 cm) exhibited a very strong positive correlation between soil depth and humus horizon thickness (r = 0.978, n = 33, p < 0.001). Conversely, no relationship was found in moderately deep soils (30–60 cm) (r = 0.018, n = 14, p < 0.948). For deep soils, we identified a moderately strong positive correlation (r = 0.345, n = 24, p = 0.098). While slope and slope length showed relationships with soil depth and humus thickness, the LS factor did not exhibit a clear correlation. These findings underscore the importance of understanding soil dynamics in informing land management practices, especially in areas susceptible to erosion. Recommendations include continued monitoring of eroded soils and implementing erosion control measures to maintain soil health and sustainability in agricultural production amidst climate change challenges.

Soil Structure Interaction of R.C. Building with Basement Resting on Pile Raft Foundation

Abstract: Historically, Basements or underground stories were built in a castle to use as a dunged on oar store room. However, in modern construction, the restrain to go deeper below the grade level in term so basements which can be utilized for parking, shopping mall or combination of both. In such cases, dynamic soil properties have a significant effect of activating dynamic soil structure interaction phenomena on during earthquake. Here in present study an effort is made to study the behavior of a building by varying number so basements considering dynamic soil structure interaction. Issues like considering higher frequency modes, influence zone to be considered for dynamic soil structure interaction, behavior of building with basements under different water level conditions for two different types of layered soil and their comparison with fixed based structure is deal with. It is observed that dynamic soil structure interaction can significant change the behavior of the building and hence it is recommended to perform dynamic soil structure interaction for building with multiple basements. In addition, some important recommendations are provided at the end to serve as a guide for researchers and practicing engineers.

A model for saturated–unsaturated flow with fractures acting as capillary barriers

AbstractHigh‐resolution modeling of the flow dynamics in fractured soils is highly complex and computationally demanding as it requires precise geometrical description of the fractures in addition to resolving a multiphase free‐flow problem inside the fractures. In this paper, we present an idealized model for saturated–unsaturated flow in fractured soils that preserves the core aspects of fractured flow dynamics using an explicit representation of the fractures. The model is based on Richards’ equation in the matrix and hydrostatic equilibrium in the fractures. While the first modeling choice is standard, the latter is motivated by the difference in flow regimes between matrix and fractures, that is, the water velocity inside the fractures is considerably larger than in the soil even under saturated conditions. On matrix/fracture interfaces, the model permits water exchange between matrix and fractures only when the capillary barrier offered by the presence of air inside the fractures is overcome. Thus, depending on the wetting conditions, fractures can either act as impervious barriers or as paths for rapid water flow. Since in numerical simulations each fracture face in the computational grid is a potential seepage face, solving the resulting system of nonlinear equations is a nontrivial task. Here, we propose a general framework based on a discrete‐fracture matrix approach, a finite volume discretization of the equations, and a practical iterative technique to solve the conditional flow at the interfaces. Numerical examples support the mathematical validity and the physical applicability of the model.

Identifying the temporal variation of nonlinear seismic site response by Artificial Neural Network

3135 strong motion records at 50 KiK-net stations in northeastern Japan are collected. The soil site response to strong motions is classified as linear, transition between linearity and nonlinearity, and nonlinear, by the dynamic soil parameters PGArotD50 and Iγ. 326 records in the nonlinear set are divided into the training set (200 records) and the test set (126 records) for an artificial neural network (ANN). This ANN is constructed to identify the temporal variation of nonlinear seismic site response automatically, by combining with the soil dynamic parameters PGArotD50 and Iγ in moving time window. The accuracy of soil nonlinearity identified using this network reaches 98%, which is higher than the result of visual identification, and there is no overfitting.

Optimized fertilization using online soil nitrate data

Abstract. A new soil nitrate monitoring system that was installed in a cultivated field enabled us, for the first time, to control the nitrate concentration across the soil profile. The monitoring system was installed in a full-scale agricultural greenhouse setup that was used for growing a bell pepper crop. Continuous measurements of soil nitrate concentrations were performed across the soil profile of two plots: (a) an adjusted fertigation plot, in which the fertigation regime was frequently adjusted according to the dynamic variations in soil nitrate concentration, and (b) a control plot, in which the fertigation was managed according to a predetermined fertigation schedule that is standard practice for the area. The results enabled an hourly resolution in tracking the dynamic soil nitrate concentration variations in response to daily fertigation and crop demand. Nitrate–nitrogen (N–NO3) concentrations in and below the root zone, under the control plot, reached very high levels of ∼ 180 ppm throughout the entire season. Obviously, this concentration reflects excessive fertigation, which is far beyond the plant demand, entailing severe groundwater pollution potential. On the other hand, frequent adjustments of the fertigation regime, which were carried out under the adjusted fertigation plot, enabled control of the soil nitrate concentration around the desired concentration threshold. This enabled a substantial reduction of 38 % in fertilizer application while maintaining maximum crop yield and quality. Throughout this experiment, decision-making on the fertigation adjustments was done manually based on visual inspections of the soil's reactions to changes in the fertigation regime. Nevertheless, it is obvious that an algorithm that continuously processes the soil nitrate concentration across the soil profile and provides direct fertigation commands could act as a “fertistat” that sets the soil nutrients at a desired optimal level. Consequently, it is concluded that fertigation that is based on continuous monitoring of the soil nitrate concentration may ensure nutrient application that accounts for plant demand, improves agricultural profitability, minimizes nitrate down-leaching and significantly reduces water resource pollution.

Sustainable strategy for rehabilitating phosphate mining sites and valorisation of phosphate industry by-products and sludge using pistachio tree (Pistacia atlantica), false pepper (Schinus molle), and eucalyptus (Eucalyptus globulus) trees

The development of anthroposols has been proposed as a new environmentally friendly approach to ensuring the successful revegetation of phosphate mining sites. The phosphate industry's by-products, including phosphogypsum (PG), phosphate sludge (PS), and sewage sludge (SS), can be valuable resources in restoring the ecological balance of mined soil areas. The aim of this study was to safely and sustainably restore the ecological integrity of the phosphate mining site through the evaluation of nutrients and heavy metals dynamics in soil and plant tissues of three tree species and treated by-products containing 65 % PG, 30 % PS, and 5 % SS. The tree species used were Pistacia atlantica, Schinus molle, and Eucalyptus globulus. The experimental layout was a randomised complete block design with six replicates and three treatments. Growth diameter, height, nutrient uptakes and heavy metal dynamic were evaluated from the rhizosphere soils and plant tissues over two years. Hierarchical head maps of correlations between the measured growth parameters, soil and nutrient uptakes of the tree species were analysed using a phylogenetic generalised linear mixed model. S. molle and E. globulus had higher average diameter and height than P. atlantica plants. P. atlantica and S. molle showed greater nitrogen, phosphorus, potassium, calcium, and magnesium concentrations than E. globulus trees. Tree growth parameters were closely linked to soil nutrient bioavailability. The heavy metal accumulation ratio was higher in the E. globulus and S. molle leaves than in stems. Using by-products could be valorised for rehabilitating mine sites together with E. globulus and S. molle species.

Evaluation of dynamic soil stress distribution in GRS bridge abutments subjected to cyclic loading

Traffic-induced cyclic loading generates repetitive stresses and cumulative deformations on the GRS abutments, which affect the serviceability of GRS abutments. To evaluate the stress distribution of GRS abutments under cyclic traffic loading, this paper presents reduced-scale GRS abutment models constructed with sand backfill and geogrid reinforcements. The GRS abutment models were subjected to staged cyclic loading with different cyclic loading amplitudes to investigate the influences of cyclic loading amplitude, bridge superstructure load, and reinforcement vertical spacing on the dynamic soil stress distributions. The results indicate that the increase in residual stresses due to stress redistribution induced by cyclic loading is most pronounced at the top of the abutment, while there is little stress redistribution down to the foundation level. Increasing the static load of bridge superstructure or the amplitude of cyclic loading results in an increase in the incremental dynamic vertical soil stresses. Reinforcement vertical spacing does not significantly impact the incremental dynamic vertical soil stresses under cyclic loading, while the cyclic load has the most significant influence. Closer reinforcement vertical spacing could provide stronger lateral confinement, resulting in larger dynamic lateral soil stresses behind wall facing.

A novel approaches to 6th-order delay differential equations in toxic plant interactions and soil impact: beyond newton-raphson

This paper focuses on investigating a 6th-order delay differential equation root within the context of toxic interactions between competing plant populations and their impact on soil dynamics. The study introduces a novel approach for approximating solutions to nonlinear delay differential equations, drawing inspiration from the fundamental principles of Newton-Raphson’s method. This technique leverages the complex root theorem to ensure stability, enabling it to effectively handle widely dispersed roots within dynamic systems. Consequently, this approach holds considerable potential for a diverse array of applications. The analysis introduces time delay into a nonlinear dynamical system and explores the system’s threshold value. At this threshold, the dynamical system’s stability undergoes fluctuations, and observations of hopf bifurcation phenomena are made. The study also successfully identifies both real and complex roots of the dynamical system. Visualization of the dynamic system is accomplished using MATLAB-generated graphical representations. Moreover, this research’s implications extend to the realm of climate action and terrestrial ecosystems, underscoring its significance for promoting a sustainable environment and fostering healthy life on land.

Studying long term relationship between carbon Emissions, Soil, and climate Change: Insights from a global Earth modeling Framework

The persistent increase in greenhouse gas (GHG) emissions, notably carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O), since the mid-20th century has been a key driver of significant climate alterations. This study investigates the complex feedback mechanisms that both influence and are influenced by global climate dynamics, soil processes, and GHG emissions. Our statistical approach incorporates correlation measures, highlighting the limitations of such analyses, namely their inability to confirm causality, sensitivity to outliers, and the exclusive capture of linear relationships. Geographically Weighted Regression (GWR) models reveal spatial variations in the relationship between environmental factors and GHGs, while Path Analysis aids in delineating direct and indirect influences among variables. The research pinpoints significant spatial heterogeneity in the impacts of economic and environmental factors on GHGs, underscoring the necessity of localized strategies for climate change mitigation and sustainable land management. This study also identifies potential threats to agricultural productivity due to soil degradation, which hinder climate adaptation efforts. Our findings advocate for a concerted global response to reduce GHG emissions and address the challenges posed by the interplay of climate change, soil dynamics, and GHG emissions.

Effect of plant-soil system on the restoration of community stability after wildfire in the northeast margin of Qinghai-Tibet plateau

Wildfires, as an environmental filter, are pivotal ecological disturbances that reshape plant communities and soil dynamics, playing a crucial role in regulating biogeographic patterns and ecosystem services. In this study, we aim to explore the effects of wildfires on forest ecosystems, specifically focusing on the plant-soil feedback mechanisms within the northeastern margin of the Qinghai-Tibet Plateau (QTP). Utilizing Partial Least Squares Path Modeling (PLS-PM), we investigated the interrelationships among soil physicochemical properties, enzyme activities, species diversity, and community stability at varying post-fire recovery stages (5, 15, and 23 years). Results indicated that in the early recovery stages, rapid changes in soil properties such as decreased pH (p < 0.001) and increased nutrient availability facilitate the emergence of early successional species with high resource utilization traits. As the ecosystem evolved toward a climax community, the soil and vegetation exhibit increased stability. Furthermore, soil enzyme activities displayed dynamic patterns that corresponded with changes in soil nutrient content, directly influencing the regeneration and diversity of plant communities. Importantly, our study documented a transition in the influence of soil properties on community stability from direct positive effects in initial recovery phases to negative impacts in later stages, while indirect benefits accrue through increased species diversity and enzyme activity. Vegetation composition and structure changed dynamically with recovery time during community succession. Plant nutrient absorption and accumulation affected nutrient dynamics in the soil, influencing plant regeneration, distribution, and diversity. Our results underscore the complex interactions between soil and vegetation that drive the recovery dynamics post-wildfire, highlighting the resilience of forest ecosystems to fire disturbances. This study contributes to the understanding of post-fire recovery processes and offers valuable insights for the management and restoration of fire-affected forest ecosystems.

Scaling of soil organic carbon in space and time in the Southern Coastal Plain, USA

Soil organic carbon (SOC) is a dynamic soil property (DSP) that represents the largest portion of terrestrial carbon. Its relevance to carbon sequestration and the potential effects of land use on SOC storage, make it imperative to map across both space and time. Most regional-scale studies mapping SOC give static estimates and train different models for different periods with varying accuracies. We developed a flexible modeling approach called DSP-Scale to map SOC in both space and time. DSP-Scale uses ecological concepts and empirical data to predict DSP dynamics using inherent soil properties (static factors) and land cover changes (dynamic factors). We compiled SOC data for the 0–20 cm depth (SOC20) from 1441 points spanning a 25 million ha study area in the southeastern U.S. Coastal Plain, incorporating data from the Rapid Carbon Assessment, National Cooperative Soil Survey Soil Characterization database, and other regional studies. We developed a random forest model using climate, topography, soil survey, and land cover changes to predict SOC20 dynamics for five-year periods between 2001 and 2019. Our model explained 66 % and 59 % of the variation for the training and test data, respectively. Top predictors included mean annual precipitation, slope, and soil erosion class. Land cover 10 years before measurements of SOC20 was more important than current land cover for estimating SOC20. We estimated total SOC stocks of 207.1 and 208.3 Tg for 2001 and 2019, respectively. Highest gains of total SOC stock (0.9 Tg from 2001 to 2019) were associated with land cover change from mixed to evergreen forest. The greatest loss of total SOC stock (0.2 Tg) in the same period was associated with land cover change from pasture/hay to evergreen forest. We concluded that the DSP-Scale approach provides a flexible way to use dynamic and static factors affecting SOC stocks to predict changes in space and time at regional scales.

"Ectomycorrhizal exploration type" could be a functional trait explaining the spatial distribution of tree symbiotic fungi as a function of forest humus forms.

In European forests, most tree species form symbioses with ectomycorrhizal (EM) and arbuscular mycorrhizal (AM) fungi. The EM fungi are classified into different morphological types based on the development and structure of their extraradical mycelium. These structures could be root extensions that help trees to acquire nutrients. However, the relationship between these morphological traits and functions involved in soil nutrient foraging is still under debate.We described the composition of mycorrhizal fungal communities under 23 tree species in a wide range of climates and humus forms in Europe and investigated the exploratory types of EM fungi. We assessed the response of this tree extended phenotype to humus forms, as an indicator of the functioning and quality of forest soils. We found a significant relationship between the relative proportion of the two broad categories of EM exploration types (short- or long-distance) and the humus form, showing a greater proportion of long-distance types in the least dynamic soils. As past land-use and host tree species are significant factors structuring fungal communities, we showed this relationship was modulated by host trait (gymnosperms versus angiosperms), soil depth and past land use (farmland or forest).We propose that this potential functional trait of EM fungi be used in future studies to improve predictive models of forest soil functioning and tree adaptation to environmental nutrient conditions.