To overcome the limitations of microscale experimental techniques and molecular dynamics (MD) simulations, a coarse-grained molecular dynamics (CGMD) method was used to simulate the wetting processes of clay aggregates. Based on the evolution of swelling stress, final dry density, water distribution, and clay arrangements under different target water contents and dry densities, a relationship between the swelling behaviors and microstructures was established. The simulated results showed that when the clay-water well depth was 300 kcal/mol, the basal spacing from CGMD was consistent with the X-ray diffraction (XRD) data. The effect of initial dry density on swelling stress was more pronounced than that of water content. The anisotropic swelling characteristics of the aggregates are related to the proportion of horizontally oriented clay mineral layers. The swelling stress was found to depend on the distribution of tactoids at the microscopic level. At lower initial dry density, the distribution of tactoids was mainly controlled by water distribution. With increase in the bound water content, the basal spacing expanded, and the swelling stresses increased. Free water dominated at higher water contents, and the particles were easily rotated, leading to a decrease in the number of large tactoids. At higher dry densities, the distances between the clay mineral layers decreased, and the movement was limited. When bound water enters the interlayers, there is a significant increase in interparticle repulsive forces, resulting in a greater number of small-sized tactoids. Eventually, a well-defined logarithmic relationship was observed between the swelling stress and the total number of tactoids. These findings contribute to a better understanding of coupled macro-micro swelling behaviors of montmorillonite-based materials, filling a study gap in clay-water interactions on a micro scale. (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/).

We report an innovative method of extracting water distribution network (WDN) historical repair location data from images of paper repair data maps, to provide usable geo-referenced digitally formatted data. For most water utilities, repair location data typically fall into two eras: pre- and post-GIS, approximately corresponding to pre- and post-2000. Automated conversion of pre-2000 paper maps to a geo-referenced digital format provides additional data to clarify trends in pipe repair causative factors, such as material defects, corrosive or creeping soils, and traffic. We applied the methodology to more than 3,000 maps of the Los Angeles Department of Water and Power WDN, thereby extending the record of repairs backward from 2000 to 1975, almost doubling the number of repair records. The methodology's value, when using the resulting data for analysis, lies in the following: (a) large volumes of hard copy data can now be acquired in an automated manner, saving significant time and effort, (b) specific repair locations are accurately captured, resulting in (c) more reliable, confident, analyses, and results, (d) ongoing problem areas, due to corrosive or creeping soils for example, can be more specifically understood.

Understanding pore water distribution in soil is essential for elucidating water movement and mechanical properties, as it significantly influences soil strength and stability. Accurate assessment of this distribution provides a scientific foundation for civil engineering design, ensuring structural safety and durability. This study examines pore water distribution using plate load tests and Nuclear Magnetic Resonance (NMR). Results indicate that matric suction expels free water first, leaving bound water until a critical suction point is reached. As matric suction increases, the peak value of the T2 relaxation time curve decreases, shifting leftward, reflecting water drainage from larger to smaller pores. Then, water expulsion occurs in three stages, with Stage III primarily indicating bound water content, quantified at 19.23%, including 3.3% as strongly bound water. An equation is derived to calculate the surface relaxation rate of 0.0176 mu m/ms. Thus, the distribution of T2 relaxation time can be transformed into pore size distribution, summarizing the characteristics of pore water distribution during the drying process. Finally, comparative analysis confirms the effectiveness of NMR in measuring bound water. These findings enhance our understanding of soil water distribution and highlight the need for advanced models that incorporate pore connectivity and water retention dynamics.

Quantifying the magnitude and distribution of degree of saturation (Sr) S r ) in unsaturated soils is crucial to understand the grain-scale hydromechanical behavior, but it has been a major experimental challenge. This study proposes a new method to quantify the pore-water distribution non-uniformity, in terms of S r , based on threedimensional (3D) X-ray computed tomography images. The algorithm constructs vectors that consider the 3D spatial distribution of Sr r for each REV. A weighted water distribution tensor (G, G , characterizing the spatial distribution of S r ) was derived to calculate a scalar parameter A that represents the non-uniformity. Application of the algorithm on sand samples demonstrated that the Sr r distributions could be highly non-uniform at different drying states. The algorithm captured pore-water transport between the dilated and non-dilated zones of samples subjected to pre-peak shearing. The evolution of A with matric suction and axial strain showed potential in incorporating the pore-water distribution into microstructure-based constitutive models.

Ground deformation induced by frost heave is a matter of concern in cold region engineering construction since it affects surrounding structures. Frost heave, which is related to the heat-water-stress interaction, is a complicated process. In this study, a heat-water-stress coupling model was established for saturated frozen soil under different stress levels to quantify the water redistribution, heat transfer, frost heave, and water intake. An empirical formula for the soil permeability considering the confining and deviator pressures was employed as an indispensable hydraulic equation in the coupling model. The Drucker-Prager yield criterion matched with the Mohr-Coulomb criterion was employed in the force equilibrium equation to investigate the deformation due to the deviator and confining pressures. The anisotropic frost heave during unidirectional freezing was further considered in the coupling model by introducing an anisotropic coefficient. Subsequently, based on the above coupling relationship, a mathematical module in COMSOL Multiphysics was applied to calculate the governing equation numerically. Finally, the proposed model was validated through an existing frost heave experiment conducted under various temperature gradients and stress levels. The results of the freezing front, water redistribution, water intake, and frost heave ratio predicted using the proposed model were found to be consistent with the experimental results.

Leakages from damaged or deteriorated buried pipes in urban water distribution networks may cause significant socio-economic and environmental impacts, such as depletion of water resources and sinkhole events. Sinkholes are often caused by internal erosion and fluidization of the soil surrounding leaking pipes, with the formation of soil cavities that may eventually collapse. This in turn causes road disruption and building foundation damage, with possible victims. While the loss of precious water resources is a well-known problem, less attention has been paid to anthropogenic sinkhole events generated by leakages in water distribution systems. With a view to improving urban smart resilience and sustainability of urban areas, this study introduces an innovative framework to localize leakages based on a Machine learning model (for the training and evaluation of candidate sets of pressure sensors) and a Genetic algorithm (for the optimal sensor set positioning) with the goal of detecting and mitigating potential hydrogeological urban disruption due to water leakage in the most sensitive/critical locations. The application of the methodology on a synthetic case study from literature and a real-world case scenario shows that the methodology also contributes to reducing the depletion of water resources.

Evaluating the seismic resilience (SR) of water distribution systems (WDSs) can support decision-making in optimizing design, enhancing reinforcement, retrofitting efforts, and accumulating resources for earthquake emergencies. Owing to the complex geological environment, buried water supply pipelines exhibit varying degrees of corrosion, which worsens as the pipelines age, leading to a continuous degradation of their mechanical and seismic performance, thereby impacting the SR of WDSs. Consequently, this study proposes an SR evaluation method for WDSs that takes into account the corrosive environment and the service age of buried pipelines. Utilizing the analytical fragility analysis method, this research establishes seismic fragility curves for pipelines of various service ages and diameters in diverse corrosive environments, in combination with the Monte Carlo simulation method to generate seismic damage scenarios for WDSs. Furthermore, the post-earthquake water supply satisfaction is utilized to characterize the system performance (SP) of WDSs. Two repair strategies are employed for damaged pipes: assigning a single repair crew to address damages sequentially and deploying a repair crew to each damage location simultaneously, to assess the minimum and maximum SR values of WDSs. The application results indicated that the maximum decrease in SP across 36 conditions was 32%, with the lowest SR value of WDSs being 0.838. Under identical seismic intensities, the SR value of WDSs varied by as much as 16.2% across different service ages and soil conditions. Under rare earthquake conditions, the effect of the corrosive environment significantly outweighs the impact of service age on the SP of WDSs. Post-disaster restoration resources can minimize the impact of the corrosive environment and service age on the SR of WDSs.

Understanding the shrinkage behavior of bentonite considering physicochemical effects is important to assess the efficiency of buffer barriers in environmental geotechnical engineering. In this paper, shrinkage experiments were conducted on Na-bentonite specimens prepared with salt solutions at various concentrations. NMR and SEM tests were conducted to study the moisture distribution and structural evolution of specimens during the evaporation of water. After sample saturation, the porosity decreases as the pore water salinity increases due to the decreasing swelling deformation with pore water concentration during the saturation process. During drying, the shrinkage deformation of compacted bentonite is anisotropic, with larger axial strains than radial strains. At the fully dried state, the bentonite specimen prepared with distilled water is the densest due to the least crystalline salts in the specimen. At the microscale, as pore water salinity increases, pore water is distributed to smaller pores, and the microstructure is more aggregated. The saline effect on water retention and distribution is weakened as pore water evaporates, originating from physicochemical effects. The structure is also more aggregated after evaporation of pore water. Theoretically, the shrinkage behavior of Na-bentonite considering the influence of water salinity is well described from the perspective of an effective stress-based constitutive relationship.