Soil Spatial Variability Research Articles

Management zone delineation: Utilizing multiple data sources to minimize soil spatial variability in commercial potato fields under Prince Edward Island pedoclimatic conditions

Management zone delineation: Utilizing multiple data sources to minimize soil spatial variability in commercial potato fields under Prince Edward Island pedoclimatic conditions

Effects of future climate on suitability of major crops in Eastern Kansas River Basin

Climate change significantly threatens food security and the agricultural economy, particularly under rainfed conditions. This study uses the Decision Support System for Agrotechnology Transfer (DSSAT) crop simulation model to evaluate the future suitability of growing maize and soybean in the Eastern Kansas River Basin (EKSRB) under two projected climate scenarios (RCP 4.5 and RCP 8.5) from 2006 to 2099. By comparing the baseline (1990–2019) and future climates, the yield gap percentage method is employed to quantify the discrepancy between actual and potential yields. This innovative approach integrates spatial soil variability and advanced climate projections from 18 global climate model (GCMs), enhancing the accuracy of crop suitability assessments. Results indicate yield losses ranging from 23% to 57% for maize and 20% to 36% for soybean, with maize experiencing a greater yield gap than soybean, highlighting soybean’s resilience under future climatic conditions. The study identifies critical regions within the EKSRB where adaptive strategies are most needed and provides insights for policymakers to develop targeted agricultural strategies, facilitate policy planning, and select mitigation strategies for vulnerable areas. This research underscores the necessity for adaptive agricultural practices to ensure food security and sustainability, offering a robust framework that can be adapted to similar regions globally.

Stability and failure probability analysis of super-large irregularly shaped deep excavation in coastal area considering spatial variability in soil properties

ABSTRACT The objective of this paper is to investigate the effects of spatial variability in coastal silty soil properties on the stability and failure probability of deep excavation. A super-large irregularly shaped deep excavation in the coastal area in Wenzhou, China, is considered. A numerical model is established based on the random field theory, finite difference method and Monte-Carlo simulation method. Focusing on the friction angle and cohesion, the effects of horizontal and vertical correlation distances and variation coefficients on excavation stability are systematically studied. Meanwhile, the reliability evaluation approach is also applied to systematically examine the impact of the correlation distance on the surface settlement and horizontal displacement of pile. The results show that the greater correlation distance leads to a larger dispersion degree of the ground settlement curve and the displacement curve of the diaphragm wall around. The influence of the horizontal correlation distance on land settlement is more significant. The horizontal and vertical displacements of the pile top are less affected by the correlation distance. Compared with the deterministic analysis, the random analysis results considering the soil spatial variability produce larger displacement.

Spatial variability of soil and vegetation properties in integrated production systems

Spatial variability of soil and vegetation properties in integrated production systems

Stochastic hazard assessment framework of landslide blocking river by depth-integrated continuum method and random field theory

AbstractLandslide-induced barrier dams pose a threat to the safety of humans, livestock and nearby infrastructures. The efficient assessment of landslide blocking river is crucial for disaster prevention and mitigation solutions. This study proposes a novel stochastic assessment framework to evaluate the landslide blocking river through the prediction of their deposition depths and considering the heterogeneity of shear strength parameters on the potential sliding surface. The depth-integrated continuum method (DICM) is used to simulate the landslide runout process. Using an enhanced Karhunen-Loève expansion (KLE) method, the spatial variations in soil's shear strength parameters are modeled by random fields to incorporate the effects of soil's spatial heterogeneity on the landslide deposition pattern. Subsequently, the multi-response surrogate model is constructed to relate the random field variables to the deposition depths based on extreme gradient boosting (XGBoost). To improve the performance of the surrogate model, principal component analysis (PCA) and sliced inverse regression (SIR) methods are employed for the dimension reduction of output and input variables, respectively. Furthermore, the algorithm for river blockage identification is developed to search for the deposition ridges. To demonstrate the capability of the stochastic assessment framework, an example of the first Baige landslide in Tibet, China is simulated, and the affected region and deposition depths of the landslide are predicted to calculate the probability of river damming. The presented methodology provides a practical means for improving the landslide blocking river prediction and new insights for early warning and risk mitigation.

Probabilistic back analysis of reservoir landslide considering hydro-mechanical coupled observations

Precisely assessing statistical parameters to characterize spatial soil variability presents a significant challenge in probabilistic slope stability analysis, primarily due to inherent soil uncertainties and limited field-specific data. Probabilistic back analysis, recognized as an effective and reliable technique, offers a rational method for utilizing observational data to invert parameters of geomaterial properties. Nevertheless, previous investigations into parameter inversion for slope stability analysis have seldom considered the coupling of hydro-mechanical characteristics in reservoir landslides. This study proposes a novel integrated Bayesian framework to perform back analysis of reservoir landslides, incorporating the spatial variability of hydro-mechanical parameters. Within this framework, a hypoplastic constitutive model is developed to characterize the step-like deformation of the reservoir slope. The posterior knowledge of the saturated permeability coefficient and shear strength parameters is obtained by collecting field monitoring data of the underground water level and ground displacement. The Maliulin landslide, located in the Three Gorges Reservoir area of China, is adopted as an illustrative example to validate the proposed framework through back analysis of spatially variable soil parameters and probabilistic stability analysis. The results demonstrate that the proposed Bayesian updating framework significantly reduces the uncertainties associated with statistical values of soil parameters, providing more accurate and reasonable updated soil parameters for probabilistic slope stability analysis.

The Spatial Variation of Soil Structure Fractal Derived from Particle Size Distributions at the Basin Scale

The accurate characterization of soil structure is fundamental to groundwater science, environmental ecology, and Earth systems science. To address the challenge of quantifying the high spatial variability of large-scale soil structures, this study used a laser particle size analyzer to measure the distribution of soil particle size in 207 samples from ten profiles across the Daqing and Ziya River basins in the North China Plain. The quantified soil structure, expressed as soil fractal dimension D, was derived using monofractal theory. Various spatial analysis techniques, including Moran’s I index, correlation analysis heat maps, the Kolmogorov–Smirnov one-sample test, and geostatistical semivariogram function, were jointly applied to investigate the spatial variability of soil structural fractals across different depths in the piedmont plain–coastal areas of the two river basins. The results indicate the following: (1) Quantitative analysis confirms that under the influence of piedmont alluvial and fluvial dynamics, soil D values homogenize from the piedmont to the coastal areas, with decreasing particle size differences closer to the coast. However, the spatial variability of the soil structural fractals in the Ziya River Basin was greater than that in the Daqing River Basin. (2) The combined effects of climate change, regional differences, and human activity led to greater spatial variability in the soil structural fractals in the Ziya River Basin than in the Daqing River Basin. The correlation between D values and burial depth was strongest in the Xianxian profile (−0.78), whereas the spatial correlation was strongest in the Hengshui and Dacheng profiles (−0.47). (3) The greatest spatial variability in soil D values occurred at depths of 1–2 m, with a coefficient of variation of 23.595%, which was significantly higher than those at depths of 0–1 (14.569%) and 2–3 m (16.284%).

Influence of a firm stratum on probabilistic slope stability analyses with spatially variable undrained soil strength

The influence of a firm stratum on the stability of a slope under undrained conditions has long been of interest to geotechnical investigators, which has been studied in a number of previously important works in relation to slope stability analyses without considering soil spatial variability. This paper proposes another look at such a problem in the context of probabilistic slope stability analyses considering soil spatial variability. Here, the random field (RF) is used to simulate the spatially variable undrained soil strength. It is found that under stationary RF and non-stationary RF with the soil strength at the top ground surface (su0) larger than 0, the depth of the firm stratum (Hf) has a significant influence on the mean and standard deviation of factor of safety (i.e., μ[FS] and σ[FS], respectively). By contrast, under non-stationary RF with su0 = 0, Hf has a slight influence on μ[FS], but its influence on σ[FS] is non-negligible. In addition, the autocorrelation distance is found to have an insignificant impact on the influential effect of Hf on μ[FS]. However, for σ[FS], this impact is not negligible. When the autocorrelation distance is smaller, the influence of Hf on σ[FS] would be more significant. Under non-stationary RF, the influence of Hf on σ[FS] would be slight if the autocorrelation distance is large enough. Furthermore, the impacts of slope ratio, su0, isotropic and anisotropic features on the influential effects of Hf are also investigated and discussed.

Assessing the performance of various infiltration models to improve water management practices

AbstractInfiltration plays a key role in stormwater management and irrigation scheduling. A review of previous studies reveals that the effectiveness of infiltration models varies significantly depending on soil characteristics and field conditions. Accurate predictions depend on selecting appropriate models for specific sites because of soil spatial variability. This requires extensive testing and recording of infiltration rates at each location. This study assesses various infiltration rate measurement models to enhance water management efficiency. Infiltration rate measurements were conducted at three sites in Dehradun using a double-ring infiltrometer. Well-established models, such as Philips JR, Green, Ampt, Horton, Kostiakov, modified Kostiakov, and the Soil Conservation Service (SCS) model, were evaluated. Data from infiltration tests were used to calibrate these models, facilitating better irrigation system design and stormwater management. In assessing their effectiveness and efficiency, key evaluation criteria such as Nash–Sutcliffe efficiency (NSE), R-squared (R2), root mean square error (RMSE), mean absolute error (MAE), and mean bias error (MBE) were employed. Our findings highlight the superiority of the Philips JR model, offering the highest overall accuracy with the highest average value R2 = 0.9557 and NSE = 0.9553, lowest MAE = 0.6717 cm/h, MBE = − 0.0160 cm/h and RMSE = 1.0077 cm/h. These results underscore the model’s ability to synthesize infiltration data effectively, even in the absence of direct measurements. This insight positions the Philips JR model as a valuable tool for estimating infiltration rates in similar conditions.

Delineating vineyard management zones: Intrafield spatial variability of soil properties of carbonate vineyard soils

AbstractGrapevine (Vitis vinifera L.) is a traditional crop cultivated in Navarre (NE Spain). However, in some areas, it is grown without harnessing land suitability for its cultivation. This research was conducted to approach the pedological recognition of viticulture zoning (on the farm scale) in a traditional and distinct viticultural region: Olite (Navarre). As grape yield and grape quality in a given field are generally variable and do not coincide one way or another, 13 soil profiles were selected for pedological description and analysis in an attempt to recognise the importance of soil properties. For that purpose, 45 soil samples (corresponding to the different pedogenetic horizons of the 13 soil profiles) were collected to improve zonal vineyard estimations. The most notable characteristics of the studied soils were the presence of petrocalcic horizons, high stone fragments content, mainly loamy textures, high pH (between 8.24 and 9.24), high carbonate (between 19.1% and 90.0%), and active limestone contents (between 5.7% and 26.1%), and relatively low organic matter contents (<3.34%). Appreciable soil properties variability was detected from the results of this study and, therefore, variability in production and grape composition was expected. These results emphasise the spatial variability of the study area soils in a way that allows for the delineation of homogeneous viticultural zones. The results also provide the necessary information not only for viticultural zoning in the Navarre wine‐growing region, but also in wine‐growing regions with a Mediterranean continental climate. Hence, our findings will allow future viticultural management zones to be developed and specific practices to be implemented.

Effect of soil spatial variability on the stability of pipelines under horizontal loading

Submarine pipeline has been utilized more and more widely in offshore oil and gas transportation engineering. Given that these pipelines frequently subject horizontal loading during operation, thus precise predicting their horizontal bearing capacity is essential. However, most previous studies employ deterministic analysis without considering soil spatialvariability. This article introduces a two-dimensional (2D) random finite element analysis to study the horizontal pullout behaviour and failure mechanisms of a pipeline in a spatial variability soil. It is observed that the effect of soil spatial variability on the horizontal bearing capacity factor and failure mechanisms is remarkably substantial, with the distribution of the former showing a close approximation to a log-normal distribution. Furthermore, the average horizontal bearing capacity typically falls below the corresponding deterministic value. Consequently, the deterministic estimation of the horizontal bearing capacity tends to be inflated. Eventually, a reliability-based approach is formulated, where a factor of safety is incorporated to account for the uncertainties arising from spatial variability and ensure a sufficient safety margin. This methodology enables engineers to determine appropriate factor of safety tailored to the specific requirements of each project.

A novel efficient RFEM for reliability analysis and design of multi-line dynamically installed anchor for floating offshore wind turbines

A novel multi-line dynamically installed anchor was previously proposed by the authors to allow for the mooring of multiple floating offshore wind turbine, resulting in a significant reduction in the total number and cost of anchors required for floating wind farms. Considering the spatial variability of soil properties and the uncertainty of environmental loads, the present study performs reliability analysis and design of the multi-line dynamically installed anchor. Firstly, a strategy to repeatedly use fundamental random variables is proposed and validated for reducing the number of random variables used in Karhunen-Loève expansion in the simulation of random field of soil properties when the ratio of the soil domain dimension to the scale of fluctuation is large. Then, the efficiency, accuracy, and robustness of the RFEM (random finite element method) combined with K-MCS (Kriging model and Monte Carlo simulation) based on the proposed strategy are validated through examples of random capacity of foundations. Next, the random capacities and probabilistic VHMT failure envelopes of the multi-line dynamically installed anchor in spatially variable soil are investigated. Finally, the reliability design of multi-line dynamically installed anchors is conducted and compared with that of multi-line pile anchors, in which both the spatial variability of soil and the uncertainty of loads are condidered. The results show that the costs for multi-line dynamically installed anchors are obviously less than those of multi-line pile anchors.

Response surface method-based efficient approach for the load and resistance factors calibration of breakwater foundation considering soil spatial variability

This study aims to propose an efficient approach to calibrate the load and resistance factors (LRFs) for the limit state design of breakwater foundations. The spatial variability of soil was considered by employing the Cholesky decomposition technique. To overcome the low computational efficiency of the Monte Carlo simulation (MCS), a response surface method (RSM) was utilized. Additionally, an efficient algorithm, namely the adaptive boundary-based particle filtering algorithm (ABB-PFA), was proposed with the use of RSM to automatically and quickly search optimal nominal values of random variables in the LRF calibration process. Noticeably, the optimal number of particles and values of the modification factor used in the ABB-PFA were observed to enhance the efficiency of the calibration process. The results showed that the computational efficiency of the calibration process is significantly increased by utilizing the response surface method and the proposed ABB-PFA. The efficiency increases the most if the number of particles is one and the value of the modification factor is 0.5 or 0.6. The influences of the random field element size on the efficiency of the calibration process were also investigated. By using the appropriate size, the efficiency is significantly increased while still obtaining a high accuracy of the LRFs calibration process.

Mapping of on-field soil nutrient variabilities as a guiding force for smart farming: a case study from FarmerZone sentinel-1 from three potato agroecological zones of India.

Mapping of soil nutrient parameters using experimental measurements and geostatistical approaches to assist site-specific fertiliser advisories is anticipated to play a significant role in Smart Agriculture. FarmerZone is a cloud service envisioned by the Department of Biotechnology, Government of India, to provide advisories to assist smallholder farmers in India in enhancing their overall farm production. As a part of the project, we evaluated the soil spatial variability of three potato agroecological zones in India and provided soil health cards along with field-specific fertiliser recommendations for potato cultivation to farmers. Specifically, 705 surface samples were collected from three representative potato-growing districts of Indian states (Meerut, UP; Jalandhar, Punjab and Lahaul and Spiti, HP) and analysed for soil parameters such as organic carbon, macronutrients (NPK), micronutrients (Zn, Fe, Mn, and Cu), pH, and EC. The soil parameters were integrated into a geodatabase and subjected to kriging interpolation to create spatial soil maps of the targeted potato agroecological zones through best-fit experimental semivariograms. The spatial distribution showed a deficiency of soil organic carbon in two studied zones and available nitrogen among all studied zones. The available phosphorus and potassium varied among the agroecological zones. The micronutrient levels were largely sufficient in all the zones except at a few specific sites where nutrient advisories are recommended to replenish. The general management strategies were recommended based on the nutrient status in the studied area. This study clearly supports the significance of site-specific soil analytics and interpolated spatial soil mapping over any targeted agroecological zones as a promising strategy to deliver reliable advisories of fertiliser recommendations for smart farming.

Calibrating resistance factors of pile groups based on individual pile proof load tests

Pile load tests have been utilized to reduce the uncertainty of pile resistance, thus leading to a higher resistance factor used in the Load and Resistance Factor Design (LRFD). Previous studies have primarily focused on calibrating resistance factors for single piles based on load tests. This calibration hinges upon the resistance bias factor of single piles, defined as the ratio of measured resistance to predicted resistance. Due to the redundancy in the pile group system, it is conventionally assumed that if the individual piles within the group achieve a lower reliability index (e.g., 2.0–2.5), the pile group as a whole attains the target reliability index of 3. However, the approach is empirical as it does not consider system redundancy directly. Moreover, this empirical approach disregards the correlation between resistance bias factors of individual piles, which is inherently influenced by the spatial variability of soils. In this study, the random finite difference method (RFDM) is employed to evaluate the correlation between resistance bias factors of individual piles in spatially variable soils. The resultant correlation matrix is subsequentially employed in Bayes’ theorem to update resistance bias factors using individual pile load test results and their corresponding test locations. The updated resistance bias factors are then used for the direct calibration of resistance factors for pile groups within the framework of LRFD. A pile group subject to vertical loading in undrained clays is adopted for illustration. Comparative analyses between the proposed approach and the empirical approach demonstrate that the latter tends to overestimate the resistance factor. Furthermore, the proposed approach enables the determination of optimal locations for conducting subsequent load tests based on previous test results.

Computational modeling of near-fault earthquake-induced landslides considering stochastic ground motions and spatially varying soil

Landslides represent a large-deformation process influenced by various factors of uncertainty, such as types of ground motion (GM), randomness of GMs, and spatial variability of soils. The behavior of landslides is profoundly affected, making their assessment and measurement challenging. This paper proposes a computational approach for simulating the large-deformation process of landslides based on three-dimensional (3D) non-ordinary state-based peridynamics (NOSBPD). The analysis results of NOSBPD indicate that the mean runout distance triggered by pulse-like ground motions (PLGMs) is 16% greater than that induced by non-pulse ground motions (NPGMs), suggesting that PLGMs exhibit higher destructiveness. Besides, this paper establishes a runout distance assessment framework by considering the stochastic nature of PLGMs. This framework allows for the evaluation of the specific risk probability of landslides caused by stochastic PLGMs that match the target spectrum specified by codes. By introducing the theory of random fields and implementing a coupled procedure in peridynamics, we conducted a detailed analysis of the impact of spatial heterogeneity on the evolution process and consequences of landslides. Additionally, compared to two-dimensional (2D) analysis, the mean runout distance obtained from 3D analysis increased by 27.5%. This suggests that 2D analysis may underestimate the consequences of landslides. The findings of this study can serve as a scientific foundation for predicting the extent and scope of landslides triggered by near-fault earthquakes.

Analysis of sensitive clay landslide mechanism and impact pipeline bearing characteristics considering spatial variability

Landslides caused by the progressive failure of sensitive clay slopes may cause substantial economic losses and environmental pollution to pipelines. This study utilizes the Karhunen-Loève (K-L) expansion method, the coupled Eulerian-Lagrangian (CEL) method, and the Monte Carlo method to model the large deformation behaviour of slopes with spatial variability. The proposed RCEL-MC framework for joint analysis of extensive deformation-probability statistics of landslide impact on pipelines confirms that the spatial variability clay slope instability failure mechanism, considering the soil softening effect under earthquake action, aligns more closely with real-world scenarios. Furthermore, based on the statistical failure probability of pipelines impacted by landslides, a sensitive clay slope pipeline safety design method is suggested. The research highlights that the random field parameters significantly influence the load-bearing characteristics of pipelines affected by spatially variable slope progressive failure landslides. The uncertainty surrounding the maximum impact force on pipelines notably increases with the rise of the horizontal correlation length and variability coefficient. When accounting for spatial variability, the maximum impact force on landslide pipelines surpasses the deterministic model analysis results by 28-69 %. It is emphasized that the deterministic model analysis, which neglects the spatial variability of soil and the soil softening effect, may seriously underestimate landslide risk. Based on the criteria for pipeline buckling failure and yield failure, it is advised that the safe design parameters for pipelines in practical projects should fall within the intersection of the secure regions defined by the two failure criteria.

Influence of geological uncertainty and soil spatial variability on tunnel deformation and their importance evaluation

In geotechnical engineering, geological uncertainty and the inherent spatial variability of soil properties are two primary sources of uncertainty, significantly impacting embedded geo-structures. While the influence of the soil spatial variability and geological uncertainty has been separately investigated, the coupling effect of these two types of uncertainties on tunnel is still limited. Meanwhile, little effort has been devoted to quantifying the degree of importance of these two types of uncertainties. This study quantitatively evaluates the effects of geological uncertainty and soil spatial variability on tunnels, and first assess the importance of these two uncertainties based on a novel approach. Improved coupled Markov chain model was utilized to simulate the geological uncertainty based on sparse and limited boreholes. The coefficient of variation (COV) was adopted to characterize the different levels of soil spatial variability in this study. The results illuminated the role of geological uncertainty and spatial variability in tunnel probabilistic analysis by comparing different boreholes schemes and COV. In addition, the geological uncertainty index (GUI) was presented to quantitatively estimate the degree of uncertainty in the simulated strata. A novel approach for evaluating the significance of geological uncertainty and soil spatial variability on tunnel was proposed. Based on the Kendall correlation coefficient, the classification was also given, divided into high, strong, moderate, weak and no correlation. The geological uncertainty should be considered in high and strong correlation groups. Both two types of uncertainties should be simultaneously taken into account in moderate and weak correlation divisions.

Effects of soil spatial variability on the behaviour of the embankment supported with a combined retaining structure

In this study, the effects of soil spatial variability on the behaviour of the embankment supported with a combined retaining structure (CRS) were investigated. The numerical model of the CRS embankment was established and validated with the field data. An application programming interface (API) was developed to deal with the data exchanging issue between the numerical model and the spatial variability characterization model. Based on the verified numerical model and the API, the probabilistic analysis with 500 Monte Carlo simulations was automatically computed. Three influencing factors of the retained soil (the mean of the friction angle, the variation of the friction angle and the vertical correlation length of the random field) are analysed by parametric analysis. The results show that the vertical correlation length of the random field is most important in the earth pressure calculation process, while the mean of the friction angle is the factor with least impact. On the whole, the spatial variability of soil properties has minimal impact on the distribution and magnitude of earth pressure behind the retaining structure.

Seismic stability and failure mechanisms of gentle and steep slopes considering soil rotated anisotropy and spatial variability

The present study focuses on investigating the effects of soil rotated anisotropy and spatial variability on slope failure in seismic conditions. The random finite element method aided by subset simulation is implemented, which ensures an efficient quantification of both the probability of failure and its associated consequence of failure. Several slope angles of gentle and steep slopes and soil properties that lead to a low probability of failure were selected for parametric studies. The comprehensive parametric studies of seismic slope stability analysis consider various factors such as slope inclinations, seismic coefficients, scales of fluctuation, and orientations of the major principal scale of fluctuation (i.e., the rotation angle of random field). The results underscore the importance of considering the combined effects of soil anisotropy and the orientation of the major principal scale of fluctuation for designing both gentle and steep slopes under seismic conditions and emphasize that care must be given to ensure the worst-case scenario is considered. Visual observations of failure mechanisms of gentle and steep slopes under seismic conditions were also shown to be very helpful in interpreting the variation of probability of failure due to soil spatial variation.