Multilayered Soil Research Articles

Modeling Borehole Resistivity Logs with Exponentially Graded Multilayered Soil Inversion for Electrical Engineering Applications

In this paper, we present the application of a special soil resistivity model, the Exponentially Graded Multilayered Soil model, to the interpretation of Borehole Resistivity Logging surveys. The model, which was recently introduced by the authors, offers an alternative to the interpretation provided by constant resistivity layer models. The model is applied to synthetic soil surveys as well as to a real survey that is not satisfactorily interpreted using conventional models. This model aids in the design and calculation of grounding systems in electrical engineering, providing an alternative to calculations performed by traditional layered models.

Improvements in the land surface configuration to better simulate seasonal snow cover in the European Alps with the CNRM-AROME (cycle 46) convection-permitting regional climate model

Abstract. Snow cover modeling remains a major challenge in climate and numerical weather prediction (NWP) models even in recent versions of high-resolution coupled surface–atmosphere (i.e., at kilometer scale) regional models. Evaluation of recent climate simulations, carried out as part of the WCRP-CORDEX Flagship Pilot Study on Convection (FPSCONV) with the CNRM-AROME convection-permitting regional climate model at 2.5 km horizontal resolution, has highlighted significant snow cover biases, severely limiting its potential in mountain regions. These biases, which are also found in AROME numerical weather prediction (NWP) model results, have multiple causes, involving atmospheric processes and their influence on input data to the land surface models in addition to deficiencies of the land surface model itself. Here we present improved configurations of the SURFEX-ISBA land surface model used in CNRM-AROME. We thoroughly evaluated these configurations on their ability to represent seasonal snow cover across the European Alps. Our evaluation was based on coupled simulations spanning the winters of 2018–2019 and 2019–2020, which were compared against remote sensing data and in situ observations. More specifically, the study tests the influence of various changes in the land surface configuration, such as the use of multi-layer soil and snow schemes, the division of the energy balance calculation by surface type within a grid cell (multiple patches), and new physiographic databases and parameter adjustments. Our findings indicate that using only more detailed individual components in the surface model did not improve the representation of snow cover due to limitations in the approach used to account for partial snow cover within a grid cell. These limitations are addressed in further configurations that highlight the importance, even at kilometer resolution, of taking into account the main subgrid surface heterogeneities and improving representations of interactions between fractional snow cover and vegetation. Ultimately, we introduce a land surface configuration, which substantially improves the representation of seasonal snow cover in the European Alps in coupled CNRM-AROME simulations. This holds promising potential for the use of such model configurations in climate simulations and numerical weather prediction both for AROME and other high-resolution climate models.

Axisymmetric consolidation behavior of multilayered unsaturated soils with transversely isotropic permeability

Time-dependent consolidation behavior of unsaturated soils is a vital problem in the geotechnical engineering. With the aid of the Fredlund consolidation theory, this work further assumes the total stress of soils skeleton freely change, and extends the Fredlund consolidation theory to a Biot-type theory, establishing the fully-coupled equation model of multilayered unsaturated poroelastic media with transversely isotropic permeability. To convert the partial differential governing equation into ordinary differential equations, the integration transform technology is applied. Subsequently, the precise integration method is used to acquire the time-dependent consolidation solution of multilayered unsaturated media with transversely isotropic permeability in the transformed domain, which is further solved in the actual domain by the inverse Hankel transform. A verification examples is provided to compare the present results with the existing work in the literature, showing a great coincidence and proving the feasibility of the present solution. Finally, numerous numerical examples are presented to investigate the evolution of excess pore pressure and settlement under quasi-static loads, revealing the consolidation behavior of unsaturated soils. The results demonstrates that the ramping time, stratification, permeability, depth and m1w have a significant effect on the consolidation behavior.

The Method of the Natural Frequency of the Offshore Wind Turbine System Considering Pile–Soil Interaction

Accurately and efficiently evaluating the influence of pile–soil interaction on the overall natural frequency of wind turbines is one of the difficulties in current offshore wind power design. To improve the structural safety and reliability of the offshore wind turbine (OWT) systems, a new closed-form solution method of the overall natural frequency of OWTs considering pile–soil interactions with highly effective calculations is established. In this method, Hamilton’s principle and the equivalent coupled spring model (ECS model) were firstly combined. In Hamilton’s theory, the Timoshenko beam assumption and continuum element theory considering the three-dimensional displacement field of soil were used to simulate the large-diameter monopile–soil interaction under lateral load in multilayer soil. Case studies were used to validate the proposed method’s correctness and efficiency. The results show that when compared with the data of 13 offshore wind projects reported in existing research papers, the difference between the overall natural frequency calculated by the proposed method and that reported in this study is within ±10%. This calculation method achieves the goal of convenient, fast and accurate prediction of the overall natural frequency of offshore wind systems.

Dynamic stress analysis of circular cavity in density inhomogeneous multi-layer soil structures containing a canyon

Dynamic stress analysis of circular cavity in density inhomogeneous multi-layer soil structures containing a canyon

Full-scale in-situ tests on a displacement cast in situ energy pile: Effects of cyclic thermal loads under different mechanical load levels on pile stress and strain

Numerous full-scale in situ tests have been conducted to assess the effect of thermal cycles on the pile response. However, those studies investigated the response of only precast and cast in-situ energy piles, with limited focus on the impact of the applied mechanical load on the pile response. This study presents the results of a field test conducted on a new type of energy pile, i.e. a displacement cast in-situ energy pile in multilayered soft soils, subjected to different fixed mechanical loads while undergoing simultaneous thermal cycles. Four tests were carried out, each corresponding to various axial loads ranging from 0 to 60% of the pile’s estimated bearing capacity. After applying the axial load on the pile head (0%, 30%, 40%, or 60% of the bearing capacity), the pile was subjected to up to ten thermal cycles. The highest magnitudes of thermal axial strains were observed near the pile top due to the lowest restraint provided by the made ground layer in all tests. Under zero (0%) mechanical load, the thermal axial strains near the pile head were elastic and recoverable, while residual strain was observed near the toe. Under reasonable working mechanical loads (30%, 40%, or 60%) residual strains were observed near both the pile head and the toe, with higher residual strains observed under higher mechanical loads. The results indicate that the cyclic thermal loadings could induce an increase in the compressive stress in the energy pile, attributed to the drag-down effects of the surrounding soil. The compressive stress induced by drag-down effects counteracts thermally induced tensile stress and thus leads to an insignificant effect on the energy pile during cooling. A limited impact of the shaft capacity was observed and was mainly attributed to the drag-down of the surrounding soil and thermal creep along the pile-soil interface.

Theoretical Analysis and Field Investigation on Bearing Characteristics of the Long‐Core SDCM Pile Under Vertical Load in Multilayered Soil

Theoretical Analysis and Field Investigation on Bearing Characteristics of the Long‐Core SDCM Pile Under Vertical Load in Multilayered Soil

A transformation and auxiliary extraction of Cr during electrokinetic removal of Cr-contaminated multilayer composite soil chamber.

Multilayer composite soil chamber was proposed to extract the Cr of contaminated site soil and insight into transformation of Cr fractionation associated with valence states. The variations of current, soil pH and moisture content were explored, as well as the migration of Cr fractionation and redistribution of Cr. Results indicated that duration of half peak current could be used to adjust treatment time and it varied among different composite ways. Moreover, extraction efficiency of Cr in soil near cathode was relatively higher and reached 60% when citric acid was used. Citric acid could promote the transformation between different Cr fractionations or different valence states. It could also improve the desorption of Cr, and could prevent excessive fluctuations of moisture content at the same time. Cr redistributed acrossed the soil chamber after extraction. When deionized water was used, Cr(VI) significantly migrated toward anode mainly in the form of exchangeable fractionation (EXC) while Fe-Mn oxides fractionation (Fe-Mn) which may be in the form of cationic Cr(III) hydroxides migrated toward cathode. When using citric acid, fractionations that were difficult to migrate of Cr, especially for Fe-Mn in site soils could be activated and became EXC and carbonate fractionation (CAR), then migrated to the anode or cathode. The migration of exchangeable Cr(III) was dramatically enhanced. But the use of citric acid could cause Cr(VI) transformation to Cr(III) near anode. In addition, during the migration process, EXC could go back to Fe-Mn again or transform to residue fractionation (RES).

Interaction analysis between single pile and multilayered saturated soils under horizontal transient loading

This study employs the coupled finite element-boundary element method to investigate the dynamic response of single pile subjected to horizontal transient loading, which is a common scenario in high-rise building, transportation, and ocean engineering. Firstly, the single pile is modeled as a Timoshenko beam and then discretized with the finite element method (FEM). The transient solution for multilayered saturated soils due to a horizontal load, serving as the kernel function for the boundary element method (BEM), aids in the derivation of the flexibility matrix of the soils. Considering the pile-soil displacement coordination conditions, the coupled FEM-BEM equation is constructed and solved to characterize the pile-soil interaction. Since the pile discretization pattern is applied consistently to the soils, it is more effective than the discretization of infinite domain in the finite element analysis. The validity of the presented method is confirmed through comparison with the full FEM numerical simulations, demonstrating the correctness and enhanced computational efficiency. Finally, parametric studies are carried out to discuss the effects of pile's length-diameter ratio, soil's shear wave velocity and stratification.

Semi‐Analytical Solution for One‐Dimensional Nonlinear Consolidation of Multilayered Soil Considering Self‐Weight and Boundary Time Effect

ABSTRACTThe self‐weight stress in multilayered soil varies with depth, and traditional consolidation research seldom takes into account the actual distribution of self‐weight stress, resulting in inaccurate calculations of soil consolidation and settlement. This paper presents a semi‐analytical solution for the one‐dimensional nonlinear consolidation of multilayered soil, considering self‐weight, time‐dependent loading, and boundary time effect. The validity of the proposed solution is confirmed through comparison with existing analytical solutions and finite difference solution. Based on the proposed semi‐analytical solution, this study investigates the influence of self‐weight, interface parameter, soil properties, and nonlinear parameters on the consolidation characteristics of multilayered soil. The results indicate that factoring in the true distribution of self‐weight leads to a faster dissipation rate of excess pore water pressure and larger settlement and settlement rate, compared to not considering self‐weight. Both boundary drainage performance and soil nonlinearity have an impact on consolidation. If the boundary drainage capacity is inadequate, the influence of soil nonlinearity on consolidation diminishes.

Changes in the spatiotemporal distribution of the timing and duration of the soil freeze–thaw status from 1979 to 2018 over the Tibetan Plateau

AbstractChanges in the soil freeze–thaw status will inevitably affect the thermal conditions and properties of the Tibetan Plateau (TP), thereby affecting its upper atmosphere, and further afield in East Asia and even globally. In this study, using the soil temperature simulated by the Community Land Model Version 5.0 (CLM5.0), the timing and duration of the soil freeze–thaw status were divided into freeze start‐date, freeze end‐date and freeze duration. Then, using linear trend estimation, correlation analysis and other methods, the changes in the spatiotemporal distribution of the timing and duration of the soil freeze–thaw status from 1979 to 2018 over the TP were analysed, and the relationships between them and surface temperature, altitude and latitude were analysed. The results obtained were as follows: (1) The soil temperature simulated by CLM5.0 can reasonably reproduce the seasonal changes in multilayer soil temperature, and correlated well with observations. Simulated by CLM5.0 of the time and duration of the soil freeze–thaw status also correlated well with observations. (2) The spatial distribution of the soil freeze–thaw status is characterized by a trend of delayed freezing, advanced thawing and shortened freeze duration from northwest to southeast over the TP. From 1979 to 2018, the freeze start‐date postponed by 7.3 days and became delayed at a rate of 1.9 days per decade, while the freeze end‐date advanced by 6.4 days at a rate of 1.7 days per decade, and the freeze duration shortened by 13.7 days at a rate of 3.6 days per decade. The timing and duration of the soil freeze–thaw status vary across different regions of the TP. The freeze start‐date in all areas of the TP has been delayed in the past 39 years. Except for the subcold zone and arid regions of the TP, the freeze end‐date has occurred earlier and the freeze duration has shortened, with the most significant changes in the subcold zone and humid regions, while the freeze end‐date has advanced at a rate of 3.6 days per decade and the freeze duration has shortened at a rate of 6.3 days per decade. (3) The timing and duration of the soil freeze–thaw status are significantly correlated with surface air temperature, elevation and latitude, exceeding the 99% confidence level. The correlation between the timing and duration of the soil freeze–thaw status and surface temperature is strongest, followed by altitude, and correlation with latitude is weaker. The correlation between surface air temperature and the timing and duration of the soil freeze–thaw status in the western TP is stronger than that in the eastern TP. The rate of change in the soil freeze–thaw status increases with altitude to 3000 m above sea level, while this rate decreases with elevation above 3000 m. The rate of change in the soil freeze–thaw status is greatest at 29°N, while the rate of delay in the freeze start‐date is minimal at 33°N and the rates of advancement in the freeze end‐date and shortening of the freeze duration are minimal at 35°N.

Analytical solution of overlying pressure for shallowly buried underwater box tunnel: A case study of box jacking in Suzhou, China

Accurately predicting the overlying pressure is crucial for determining an appropriate cover depth of underwater box tunnels to avoid the uplifting failure. Based on the project of box jacking crossing the Beijing-Hangzhou Grand Canal in Suzhou, the characteristics of overlying pressure variation during tunneling are investigated. The monitoring results reveal that the fluctuation of overlying pressure is weakened during the rapid tunneling process. A modified analytical model for vertical earth pressure is conceived, in which the active and passive limit states for multi-layered soils are both considered. The probable range of overlying pressure obtained by the proposed model is suitable to cover the actual values. The anti-floating behavior of underwater box tunnels for two different working conditions is discussed by calculating the minimum cover depth. Using the calibrated analytical models, a parametric study is conducted to explore the influence of injection pressure, hardened slurry unit weight, soil internal friction angle, soil cohesion, and tunnel geometry. It is found that the injection pressure during the construction process is crucial for determining the necessary cover depth, and the change of box tunnel height makes it easier to trigger the variation of minimum cover depth.

Galerkin boundary element method for simulating lightning response of grounding grid in horizontal multilayered soil model considering soil ionization effect

AbstractThis work proposes a novel mathematical model based on the Galerkin time‐domain boundary element method for accurately calculating the lightning current distribution and lightning impulse response of the buried substation grounding grid, in a multi‐layer horizontal layered soil model, by taking into account the soil ionization effect. To improve computational efficiency, the quasi‐static complex image method and its closed form time‐domain Green's function have been introduced into the model that has the ability to analytically calculate the mutual inductance coefficient between the branch currents of any two conductor segments and the mutual resistance coefficient between the leakage currents. The Galerkin time‐domain boundary element method proposed in this work can simulate the transient lightning impulse response of a substation grounding grid buried in the multi‐layer horizontal layered soil.

A Method for Determining the Fundamental Site Period and the Average Shear Wave Velocity

The soil–structure interaction plays a crucial role in determining the displacement and internal forces of multi-story buildings subjected to strong ground motion. One of the critical dynamic characteristics influencing soil–structure interaction is the fundamental site period and the average shear wave velocity associated with it. This study introduces an original equation to determine these parameters. In addition, for the first time in the literature, the version of the Rayleigh method used for finding the fundamental periods of buildings is used to find the fundamental site period. The soil is modeled as an equivalent shear beam to obtain the proposed equation. The peak displacement is obtained by acting the soil mass as an external load on the equivalent shear beam. For single-layer soil, the fundamental site period is proportional to the square root of the peak displacement of the equivalent shear beam. The least squares method generalizes the proposed relation for single-layer soils to multi-layer soil profiles. Modified Finite element Transfer matrix method is used for calibration in the least squares method. The equations used in the literature and earthquake codes for determining the fundamental site period and average shear velocity are tested on various examples, and it is shown that the method proposed in this study, along with the Rayleigh method, gives better results than these equations. The performances of these two methods and the five commonly used equations are tested and compared on different soil profiles. Transfer functions, Finite Element Method (SAP200) and Modified Finite Element Transfer Matrix Method are used for verification. For all soil profiles, the results obtained from the transfer function, Finite Element Method (SAP200) and Modified Finite Element Transfer Matrix Method are found to be in agreement. The true percent relative error found in the results obtained with the proposed method is 4.47%.

Soil moisture prediction based on integrating the CEEMDAN decomposition technique with LSTM in the proximity of Prosopis Juliflora

Human-made adaptations in plant diversity, particularly the introduction of alien invasive plant species, can result in problems such as depletion of groundwater, extraction of natural water resources, and modifications in the local biodiversity. In addition to regulating ecosystems and accurate farming, precise forecasting of soil moisture in the vicinity of these exotic plant species can greatly assist in various other agricultural applications. In this study, we present an improved approach for predicting soil moisture levels at different depths in the vicinity of Prosopis Juliflora (IAPS) by integrating long short-term memory (LSTM) and complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN). CEEMDAN is employed as a preprocessing and denoising technique to decompose the initial soil moisture data into several intrinsic mode functions (IMF). After that, LSTM is used to predict all the IMF (Intrinsic Mode Function) parts from the CEEMDAN process. The individual prediction outcomes of each IMF are combined to obtain the final prediction results. Soil moisture measurements were collected continuously in the Coimbatore area of Tamil Nadu at multi-layer soil depths of 10 cm, 30 cm, 60 cm, and 100 cm. The provided data was utilized for training models and achieving optimal accuracy in forecasting soil moisture levels one month prior. The proposed CEEMDAN-LSTM model achieved the highest R2, ranging from 0.985 to 0.995 at various depths when compared to RNN, LSTM, EMD-LSTM, EEMD-LSTM models, and existing literature models. The results show that the proposed method works well at predicting soil moisture at different depths, which can help improve the planting pattern and schedule of native species to increase the area’s green cover.

Implementation of a slope stability method in the CRITERIA-1D agro-hydrological modeling scheme

AbstractThis paper presents the implementation of a slope stability method for rainfall-induced shallow landslides in CRITERIA-1D, which is an agro-hydrological model based on Richards’ equation for transient infiltration and redistribution processes. CRITERIA-1D can simulate the presence and development of roots and canopies over space and time, the regulation of transpiration activity based on real meteorological data, and the evaporation reduction caused by canopies. The slope can be considered composed of a multi-layered soil, leading to the possibility of simulating the bedrock and of setting an initial water table level. CRITERIA-1D can consider different soil horizons characterized by different hydraulic conductivities and soil water retention curves, thus allowing the simulation of capillarity barriers. The validation of the proposed physically based slope stability model was conducted through the simulation of the collected water content and water potential data of an experimental slope. The monitored slope is located close to Montuè, in the north-eastern sector of Oltrepò Pavese (northern Apennines—Italy). Just close to the monitoring station, a shallow landslide occurred in 2014 at a depth of around 100 cm. The results show the utility of agro-hydrological modeling schemes in modeling the antecedent soil moisture condition and in reducing the overestimation of landslides events detection, which is an issue for early warning systems and slope management related to rainfall-induced shallow landslides. The presented model can be used also to test different bioengineering solutions for slope stabilization, especially when data about rooting systems and plant physiology are known.

The influence of plant water stress on vegetation–atmosphere exchanges: implications for ozone modelling

Abstract. Evapotranspiration is important for Earth's water and energy cycles as it strongly affects air temperature, cloud cover, and precipitation. Leaf stomata are the conduit of transpiration, and their opening is sensitive to weather and climate conditions. This feedback can exacerbate heat waves and can play a role in their spatiotemporal propagation. Sustained high temperatures strongly favour high ozone levels, with significant negative impacts on air quality and thus on human health. Our study evaluates the process representation of evapotranspiration in the atmospheric chemistry–climate European Centre for Medium-Range Weather Forecasts – Hamburg(ECHAM)/Modular Earth Submodel System (MESSy) Atmospheric Chemistry model. Different water stress parameterizations are implemented in a stomatal model based on CO2 assimilation. The stress factors depend on either soil moisture or leaf water potential, which act on photosynthetic activity, and mesophyll and stomatal conductance. The new functionalities reduce the initial overestimation of evapotranspiration in the model globally by more than an order of magnitude, which is most important in the Southern Hemisphere. The intensity of simulated warm spells over continents is significantly improved. For ozone, we find that a realistic model representation of plant water stress suppresses uptake by vegetation and enhances photochemical production in the troposphere. These effects lead to an overall increase in simulated ground-level ozone, which is most pronounced in the Southern Hemisphere over the continents. More sophisticated land surface models with multi-layer soil schemes could address the uncertainties in representing plant dynamics representation due to too-shallow roots. In regions with low evaporative loss, the representation of precipitation remains the largest uncertainty.

Upper Bound Analysis of Two-Layered Slopes Subjected to Seismic Excitations Using the Layer-Wise Summation Method

Due to natural sedimentation and artificial filling, slopes exhibit heterogeneity in the form of multi-layer soils, namely, layered slopes. Compared with homogenous slopes, the failure mechanism of layered slopes is more complex owing to the different shear strengths of each soil layer. Therefore, it is of great importance to gain insight into the stability of layered slopes. In this study, the upper bound theorem of limit analysis incorporated with a pseudo-static approach is utilized to investigate the seismic stability of two kinds of two-layered slopes: one with a stiff lower soil layer and the other with a weak lower soil layer. Three failure patterns, namely face failure, toe failure and base failure, are taken into account. A depth coefficient (Δ) is introduced to describe the distribution of two soil layers. The layer-wise summation method is adopted to calculate the safety factor and yield acceleration coefficient more conveniently. Based on Newmark’s method, the earthquake-induced horizontal displacement is estimated. The calculated results are validated by comparisons with published literature and the numerical method in terms of safety factor, critical failure surface and yield acceleration coefficient. The results show that the depth coefficient has a significant influence on the failure mechanism of two-layered slopes by determining whether the stability of upper-layered soil is dominant in the overall slope stability or not. Inaccurately identifying the failure patterns will overestimate the seismic performance of two-layered slopes in the aspects of safety factor and yield acceleration coefficient, leading to an underestimation of earthquake-induced horizontal displacement.

Undrained Bearing Capacity and Failure Mechanism of Strip Footings on Slopes Considering Multilayered Soils

Undrained Bearing Capacity and Failure Mechanism of Strip Footings on Slopes Considering Multilayered Soils

Analysis of the Dynamic Behavior of Multi-Layered Soil Grounds

The ground consists of many layers of soil with different properties. The propagation speed and path of seismic waves are affected by different soil layers. It is necessary to understand that layered soil exhibits different dynamic behaviors and responses under the action of seismic waves. This study utilized weathered soil and silica sand as materials to create multi-layered soil conditions with varying degrees of compaction. By conducting a 1 g shaking-table test on multi-layered soil, the interactions and influences between different soil layers under different earthquake conditions were observed. The approach of our numerical analysis aimed to complement the experimental results and provide an in-depth understanding of the dynamic behavior of multi-layered soil surfaces during seismic events. The acceleration results achieved with the ABAQUS and DEEPSOIL models for multi-layered soil were in good agreement with the experimental results. By comparing the stress–strain curves, the deformation mechanisms under different constitutive models in the numerical analysis were studied. The results of this study show that the amplification effect of seismic waves is related to the number of soil layers and the degree of compaction of the soil layers. This indicates that multi-layered soil ground and the behavior of the soil layers play an important role in the propagation and impact of seismic waves, and this amplification effect is of great significance in the design of actual seismic disaster risk assessments.