This study evaluates dykes stability of bauxite residue storage facility using limit equilibrium (LEM) and finite element methods (FEM), considering diverse construction phases. In LEM, steady state seepage is simulated using piezometric line while factor of safety (FOS) is determined by Morgenstern-Price method using SLOPE/W. In FEM, actual loading rates and time dependent seepage is modelled by coupled stress-pore water pressure analysis in SIGMA/W and dyke stability is assessed by stress analysis in SLOPE/W, referencing SIGMA/W analysis as a baseline model. Both the analysis incorporated suction and volumetric water content functions to determine FOS. FEM predicted pore pressures are validated against in-situ piezometer data. The results highlight that coupled hydro-mechanical analysis offers accurate stability assessment by integrating stress-strain behaviour, pore pressure changes, seepage paths, and dyke displacements with time. It is found that inclusion of unsaturated parameters in Mohr-Coulomb model improved the reliability in FOS predictions.

This paper provides an extension of an existing elasto-plastic framework originally proposed by Gens & Nova (1993) for modelling the response of structured soils and soft rocks. The model is enhanced to reproduce not only the mechanical response of caprocks under standard monotonic triaxial loading, but also the effects of the environmental and hydraulic loading induced by modern energy applications, including gas/hydrogen storage and geological carbon storage. The novelty of these applications, compared to the more usual ones developed by the oil and gas industry over decades, lies in the complex pore fluid and stress pressure histories applied and in the strong geochemical interaction of the rock formations with non-native fluids. Cyclic pore pressure histories due to seasonal gas storage may result in a mechanical degradation of the caprock material, while chemical degradation may occur due to pore water acidification resulting from the rock-water-CO2 interaction. To cope with the cyclic mechanical degradation, the framework is first coupled with the extended overstress theory, so to satisfactorily reproduce the time-dependent behaviour of caprocks, which presents inelastic strains even within the yield surface. Such an extension is shown to be essential to reproduce the strong strain-rate dependence and the increase in the number of cycles to failure with the amplitude of cyclic loading observed in experimental data obtained on intact specimens of an Italian stiff carbonatic clay. The elasto-plastic model is then enhanced to account for chemical degradation, using the calcite mass fraction dissolution as a variable controlling damage evolution. Combined with a geochemical reactive transport model, this extension satisfactorily reproduces the progressive degradation of a Chinese shale due to CO2 exposure, showing the ability of the framework to model coupled geo-chemo-mechanical processes.

Societal Impact StatementSolar parks enable renewable energy production at a large scale, thereby reducing greenhouse gas emissions. However, the effects of this change in land use on vegetation and soil health are still largely unknown. In this study, we determined the impacts of solar parks on vegetation, soil biota and soil carbon between and below solar panels. We found lower plant and microbial biomass below the panels, while no differences in soil carbon pools were observed. The results stress the urgent need to design future solar parks that prevent soil degradation while still producing the renewable energy needed to combat climate change.Summary Solar parks, large-scale arrays of photovoltaic panels, are a unique land use and play an important role in the renewable energy transition. However, the solar panels create shade and change the microclimate, potentially affecting plant growth and carbon inputs to the soil. These changes can influence key soil properties critical to long-term carbon storage and overall soil health. This study investigated the impact of commercial solar parks on plant productivity and the colonisation of roots by mycorrhizal fungi, soil organic matter (SOM), soil microbial community biomass and composition and litter decomposition in 17 solar parks with contrasting shading levels across the Netherlands. Soil samples and plant biomass samples were collected between and below the solar panels. The microclimate (temperature, moisture) was measured continuously over the growing season and cumulative solar irradiation during the growing season in relation to the solar panels was modelled. Results show that above- and below-ground plant biomass as well as mycorrhizal colonisation were significantly lower below than between panels, while we did not find differences for SOM, carbon stocks and hot water extractable carbon. Plant productivity related negatively to the extent of light interception by the panels. Furthermore, fungal and bacterial biomass and the F:B ratio were lower below compared to between the panels while decomposition rates did not differ. The severe decrease of plant biomass inputs in combination with maintained rates of decomposition are expected to result in decreased SOM stocks and soil health over time and suggest the need for guidelines for ecologically sound solar park designs to prevent soil damage.

Pumped storage power stations usually arrange galleries in the backfill area at the bottom of the reservoir basin. Under the influence of uneven deformation, the galleries may be difficult to adapt to deformation and generate cracking, which can affect dam safety. In this study, the upper reservoir of Hohhot pumped storage power station was taken as a case study. Through a combination of monitoring data and numerical simulation, the deformation characteristics of the galleries on the backfill foundation were analyzed, and the causes and mechanisms of galleries cracking and structural joints damage were revealed. The in situ monitoring records cover the internal settlement of the dam, the deformation and seepage flow of the galleries, and the ambient temperature. Based on actual engineering data, a numerical model considering the structure and filling method of dam, backfill area, and gallery was established, and the calculation parameters of rockfill material constitutive model were inverted by the direct back analysis method. The monitoring data analysis and numerical calculations showed that the long length of the gallery and the sudden drop of the ambient temperature were the main reasons for the longitudinal microcracks in the top arch of the galleries in the backfill area; the strong constraint of bedrock and the uneven settlement of backfill foundation were the root causes for the penetrating cracks in the galleries at the junction of backfill area and bedrock. In addition, the depth of the gallery embedded in the bedrock determines the deformation form (torsional deformation or bending deformation) of the galleries at the junction of the backfill area and bedrock. Based on the monitoring and numerical simulation, the long-term deformation of the galleries and the development of structural joints were also predicted.

The cyclic injection and production of fluids into and from underground gas storage (UGS) may lead to caprock failure, such as capillary sealing failure, hydraulic fracturing, shear failure, and fault slipping or dilation. The dynamic sealing capacity of a caprock-fault system is a critical constraint for safe operation, and is a key factor in determining the maximum operating pressure (MOP). This study proposed an efficient semi-analytical method for calculating changes in the in situ stress within the caprock. Next, the parameters of dynamic pore pressure, in situ stresses, and deformations obtained from reservoir simulations and geomechanical modeling were used for inputs for the analytical solution. Based on the calculated results, an experimental scheme for the coupled cyclic stress-permeability testing of caprock was designed. The stability analysis indicated that the caprock was not prone to fatigue shear failure under the current injection and production strategy, supported by the experimental results. The experimental results further reveal that the sealing capacity of caprock plugs may remain stable. This phenomenon is attributed to cyclic stress causing pore connectivity and microcrack initiation in certain plugs, while leading to pore compaction in others. A comparison between the dynamic pore pressure and the minimum principal stress suggests that the risk of tensile failure is extremely low. Furthermore, although the faults remain stable under the current injection and production strategies, the continuous increase in injection pressure may lead to an increased tendency for fault slip and dilation, which can cause fault slip ultimately. The MOPs corresponding to each failure mode were calculated. The minimum value of approximately 36.5 MPa at capillary sealing failure indicated that the gas breakthrough in the caprock occurred earlier than rock failure. Therefore, this minimum value can be used as the MOP for the target UGS. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

Numerous endorheic lakes in the Qinghai-Tibet Plateau (QTP) have shown a dramatic increase in total area since 1996. These expanding lakes are mainly located in the interior regions of the QTP, where permafrost is widely distributed. Despite significant permafrost degradation due to global warming, the impact of permafrost thawing on lake evolution in QTP has been underexplored. This study investigated the permafrost degradation and its correlation with lake area increase by selecting four lake basins (Selin Co, Nam Co, Zhari Namco, and Dangqiong Co) in QTP for analysis. Fluid-heat-ice coupled numerical models were conducted on the aquifer cross-sections in these four lake basins, to simulate permafrost thawing driven by rising surface temperatures, and calculate the subsequent changes in groundwater discharge into the lakes. The contribution of these changes to lake storage, which is proportional to lake area, was investigated. Numerical simulation indicates that from 1982 to 2011, permafrost degradation remained consistent across the four basins. During this period, the active layer thickness first increased, then decreased, and partially transformed into talik, with depths reaching up to 25 m. By 2011, groundwater discharge had significantly risen, exceeding 2.9 times the initial discharge in 1988 across all basins. This increased discharge now constitutes up to 17.67 % of the total lake water inflow (Selin Co). The dynamic lake water budget further suggests that groundwater contributed significantly to lake area expansion, particularly since 2000. These findings highlight the importance of considering permafrost thawing as a crucial factor in understanding the dynamics of lake systems in the QTP in the context of climate change.

Uncertainties in carbon storage estimates for disturbance-prone dryland ecosystems hinder accurate assessments of their contribution to the global carbon budget. This study examines the effects of land-use change on carbon storage in an African savanna landscape, focusing on two major land-use change pathways: agricultural intensification and wildlife conservation, both of which alter disturbance regimes. By adapting tree inventory and soil sampling methods for dryland conditions, we quantified aboveground and belowground carbon in woody vegetation (AGC and BGC) and soil organic carbon (SOC) across these pathways in two vegetation types (scrub savanna and woodland savanna). We used Generalized Additive Mixed Models to assess the effects of multiple environmental drivers on AGC and whole-ecosystem carbon storage (C-total). Our findings revealed a pronounced variation in the vulnerability of carbon reservoirs to disturbance, depending on land-use change pathway and vegetation type. In scrub savanna vegetation, shrub AGC emerged as the most vulnerable carbon reservoir, declining on average by 56% along the conservation pathway and 90% along the intensification pathway compared to low-disturbance sites. In woodland savanna, tree AGC was most affected, decreasing on average by 95% along the intensification pathway. Unexpectedly, SOC stocks were often higher at greater disturbance levels, particularly under agricultural intensification, likely due to the preferential conversion of naturally carbon-richer soils for agriculture and the redistribution of AGC to SOC through megaherbivore browsing. Strong unimodal relationships between disturbance agents, such as megaherbivore browsing and woodcutting, and both AGC and C-total suggest that intermediate disturbance levels can enhance ecosystem-level carbon storage in disturbance-prone dryland ecosystems. These findings underline the importance of locally tailored management strategies-such as in carbon certification schemes-that reconcile disturbance regimes in drylands with carbon sequestration goals. Moreover, potential trade-offs between land-use objectives and carbon storage goals must be considered.

We study CO2 injection into a saline aquifer intersected by a tectonic fault using a coupled modeling approach to evaluate potential geomechanical risks. The simulation approach integrates the reservoir and mechanical simulators through a data transfer algorithm. MUFITS simulates non-isothermal multiphase flow in the reservoir, while FLAC3D calculates its mechanical equilibrium state. We accurately describe the tectonic fault, which consists of damage and core zones, and derive novel analytical closure relations governing the permeability alteration in the fault zone. We estimate the permeability of the activated fracture network in the damage zone and calculate the permeability of the main crack in the fault core, which opens on asperities due to slip. The coupled model is applied to simulate CO2 injection into synthetic and realistic reservoirs. In the synthetic reservoir model, we examine the impact of formation depth and initial tectonic stresses on geomechanical risks. Pronounced tectonic stresses lead to inelastic deformations in the fault zone. Regardless of the magnitude of tectonic stress, slip along the fault plane occurs, and the main crack in the fault core opens on asperities, causing CO2 leakage out of the storage aquifer. In the realistic reservoir model, we demonstrate that sufficiently high bottomhole pressure induces plastic deformations in the near-wellbore zone, interpreted as rock fracturing, without slippage along the fault plane. We perform a sensitivity analysis of the coupled model, varying the mechanical and flow properties of the storage layers and fault zone to assess fault stability and associated geomechanical risks. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

Wood plays a vital role in the terrestrial carbon cycle, both sequestering and subsequently releasing carbon to the atmosphere via decomposition. Decomposition has largely been studied in fallen and standing deadwood; much less is known about decomposition occurring inside live trees due to hollowing by wood-feeding termites and microbial heart rot. Internal stem damage is difficult to measure, leaving many unresolved knowledge gaps. Little is known regarding the location and total amount of damage done by termites and microbes, as well as whether these decomposers act in concert or separately. Furthermore, tree species, wood density and stem size can influence fallen deadwood decomposition, but their role in living tree internal damage is largely unknown. We destructively harvested 63 trees, finding internal damage in 32. We intensively sampled the internal stem damage in these 32 to investigate the relative contributions of microbes and termites in a tropical savanna in Queensland, Australia. We tested if damage changed at different heights in the tree, quantified tree-level termite and microbial damage and examined if termite and microbial damage co-occurred. We also tested the influence of tree species, wood specific gravity and size on tree-level internal stem damage across four species. Termite and microbial damage were present in 45% and 33% of all trees, respectively. On average, termite damage reduced total tree biomass by 3.3% (maximum 28%, SD = 4.7%) and microbial damage by 1.8% (maximum 26%, SD = 5.3%). The amount of damage from both decomposers decreased with increasing heights up the tree. Termite and microbial damage co-occurrence was greater within trees than within individual cross samples, suggesting local competitive exclusion or niche partitioning by decomposers. Tree species was a better predictor of damage than either wood specific gravity or tree size. Half of the trees in our study had substantial internal stem damage, highlighting the considerable role that termites and microbes play in decomposing wood within living trees. Our findings unveil the previously concealed wood decomposition dynamics occurring inside trees, with implications for accurate carbon estimation across savanna ecosystems.Read the free Plain Language Summary for this article on the Journal blog.

This study focuses on predicting the impacts of a heating-cooling cycle on the pullout capacity of energy piles installed through a soft clay layer. Geotechnical centrifuge physical modeling was used to evaluate temperature, pore water pressure, volume change, and undrained shear strength profiles in clay layers surrounding energy piles heated to different maximum temperatures to understand their impacts on the pile pullout capacity. During centrifugation at 50 g, piles were jacked-in at a constant rate of penetration into a kaolinite clay layer consolidated from a slurry in a cylindrical aluminum container, heated to a target temperature after stabilization of installation effects, cooled after completion of thermal consolidation requiring up to 30 hours (1250 days in prototype scale), then pulled out at a constant rate. T-bar penetration tests were performed after the heatingcooling cycle to assess differences in clay undrained shear strength from a baseline test. The pullout capacity of an energy pile heated to 80 degrees C then cooled to ambient temperature was 109 % greater than the capacity in the baseline test at 23 degrees C, representing a substantial improvement. The average undrained shear strength measured with the T-bar at a distance of 3.5 pile diameters from the pile heated to 80 degrees C was 60 % greater than at 23 degrees C but followed the same trend as pile capacity with temperature. An empirical model for the pullout capacity was developed by combining predictions of soil temperature, thermal excess pore water pressure, thermal volumetric strain, and undrained shear strength for different maximum pile temperatures. The empirical model predictions matched well with measured pullout capacities.