Modelling the effect of changing precipitation inputs on deep soil water utilization

Forests in the Southeastern United States are predicted to experience future changes in seasonal patterns of precipitation inputs as well as more variable precipitation events. These climate change‐induced alterations could increase drought and lower soil water availability. Drought could alter rooting patterns and increase the importance of deep roots that access subsurface water resources. To address plant response to drought in both deep rooting and soil water utilization as well as soil drainage, we utilize a throughfall reduction experiment in a loblolly pine plantation of the Southeastern United States to calibrate and validate a hydrological model. The model was accurately calibrated against field measured soil moisture data under ambient rainfall and validated using 30% throughfall reduction data. Using this model, we then tested these scenarios: (a) evenly reduced precipitation; (b) less precipitation in summer, more in winter; (c) same total amount of precipitation with less frequent but heavier storms; and (d) shallower rooting depth under the above 3 scenarios. When less precipitation was received, drainage decreased proportionally much faster than evapotranspiration implying plants will acquire water first to the detriment of drainage. When precipitation was reduced by more than 30%, plants relied on stored soil water to satisfy evapotranspiration suggesting 30% may be a threshold that if sustained over the long term would deplete plant available soil water. Under the third scenario, evapotranspiration and drainage decreased, whereas surface run‐off increased. Changes in root biomass measured before and 4 years after the throughfall reduction experiment were not detected among treatments. Model simulations, however, indicated gains in evapotranspiration with deeper roots under evenly reduced precipitation and seasonal precipitation redistribution scenarios but not when precipitation frequency was adjusted. Deep soil and deep rooting can provide an important buffer capacity when precipitation alone cannot satisfy the evapotranspirational demand of forests. How this buffering capacity will persist in the face of changing precipitation inputs, however, will depend less on seasonal redistribution than on the magnitude of reductions and changes in rainfall frequency.     

Variations in glacier volume and snow cover and their impact on lake storage in the Paiku Co Basin,in the Central Himalayas

Glaciers and snow cover are important constituents of the surface of the Tibetan Plateau. The responses of these phenomena to global environmental changes are sensitive, rapid and intensive due to the high altitudes and arid cold climate of the Tibetan Plateau. Based on multisource remote sensing data, including Landsat images, MOD10A2 snow product, ICESat, Cryosat-2 altimetry data and long-term ground climate observations, we analysed the dynamic changes of glaciers, snow melting and lake in the Paiku Co basin using extraction methods for glaciers and lake, the degree-day model and the ice and lake volume method. The interaction among the climate, ice-snow and the hydrological elements in Paiku Co is revealed. From 2000 to 2018, the basin tended to be drier, and rainfall decreased at a rate of −3.07 mm/a. The seasonal temperature difference in the basin increased, the maximum temperature increased at a rate of 0.02°C/a and the minimum temperature decreased at a rate of −0.06°C/a, which accelerated the melting from glaciers and snow at rates of 0.55 × 10⁷ m³/a and 0.29 × 10⁷ m³/a, respectively. The rate of contribution to the lake from rainfall, snow and glacier melted water was 55.6, 27.7 and 16.7%, respectively. In the past 18 years, the warmer and drier climate has caused the lake to shrink. The water level of the lake continued to decline at a rate of −0.02 m/a, and the lake water volume decreased by 4.85 × 10⁸ m³ at a rate of −0.27 × 10⁸ m³/a from 2000 to 2018. This evaluation is important for understanding how the snow and ice melting in the central Himalayas affect the regional water cycle.     

Evaluating soil water routing approaches in watershed-scale,ecohydrologic modelling

Soil water dynamics are central in linking and regulating natural cycles in ecohydrology, however, mathematical representation of soil water processes in models is challenging given the complexity of these interactions. To assess the impacts of soil water simulation approaches on various model outputs, the Soil and Water Assessment Tool was modified to accommodate an alternative soil water percolation method and tested at two geographically and climatically distinct, instrumented watersheds in the United States. Soil water was evaluated at the site scale via measured observations, and hydrologic and biophysical outputs were analysed at the watershed scale. Results demonstrated an improved Kling–Gupta Efficiency of up to 0.3 and a reduction in percent bias from 5 to 25% at the site scale, when soil water percolation was changed from a threshold, bucket-based approach to an alternative approach based on variable hydraulic conductivity. The primary difference between the approaches was attributed to the ability to simulate soil water content above field capacity for successive days; however, regardless of the approach, a lack of site-specific characterization of soil properties by the soils database at the site scale was found to severely limit the analysis. Differences in approach led to a regime shift in percolation from a few, high magnitude events to frequent, low magnitude events. At the watershed scale, the variable hydraulic conductivity-based approach reduced average annual percolation by 20–50 mm, directly impacting the water balance and subsequently biophysical predictions. For instance, annual denitrification increased by 14–24 kg/ha for the new approach. Overall, the study demonstrates the need for continued efforts to enhance soil water model representation for improving biophysical process simulations.     

Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

Implementation of a coupled plastic damage distinct lattice spring model for dynamic crack propagation in geomaterials

This paper studies dynamic crack propagation by employing the distinct lattice spring model (DLSM) and 3‐dimensional (3D) printing technique. A damage‐plasticity model was developed and implemented in a 2D DLSM. Applicability of the damage‐plasticity DLSM was verified against analytical elastic solutions and experimental results for crack propagation. As a physical analogy, dynamic fracturing tests were conducted on 3D printed specimens using the split Hopkinson pressure bar. The dynamic stress intensity factors were recorded, and crack paths were captured by a high‐speed camera. A parametric study was conducted to find the influences of the parameters on cracking behaviors, including initial and peak fracture toughness, crack speed, and crack patterns. Finally, selection of parameters for the damage‐plasticity model was determined through the comparison of numerical predictions and the experimentally observed cracking features.     

Modelling the poroelastoplastic behaviour of soils subjected to internal erosion by suffusion

Granular soils subjected to seepage flow may suffer suffusion, ie, a selective internal erosion. Extending the classical approach of poromechanics, we deduce a new form of the Clausius-Duhem inequality accounting for dissipation due to suffusion, and we deduce restrictions on the constitutive laws of the soil. We suggest (a) a possible coupling between the seepage forces and the suffusion kinetics and (b) an extension of an existing elastoplastic model for the skeleton mechanical behaviour. Numerical integrations of the elastoplastic model are carried out under drained axisymmetric triaxial and oedometric conditions. As a result, we prove that the extended model is able to qualitatively reproduce the suffusion induced strains and the strength reduction experimentally observed. Predictions on the oedometric behaviour of suffusive soils are also provided.     

自然和自然贡献情景模型研究综述

本文回顾了自然和自然贡献情景模型发展的背景、历史和内容,概括总结了自然和自然贡献情景模型的发展进程以及联合国生物多样性与生态系统服务政府间科学—政策平台(IPBES)情景模型的概念框架,讨论了自然和自然贡献情景模型存在的问题和发展方向。为了在全球层面解决现有综合集成模型存在的问题,根据地球表层建模基本定理和生态环境曲面建模基本定理,提出了具有中国原创特点的自然与自然贡献情景模型概念框架。     

河外H2O超脉泽源及与其相关的分子谱线的观测和研究

自从在活动星系核NGC4945的视线方向上发现第一个河外H2O超脉泽源以来，迄今为止已发现了19个河外H2O超脉泽源，对与活动星系核成协和河外H2O超脉泽源及分子谱线的观测和研究是探测和研究活动星系核核区中央源，拱核气体和尘埃环性质的非常有效的工具，主要评述对河外H2O超脉泽源及与其相关分子谱线的搜索，观测和理论研究现状。     

Stochastic modelling of groundwater flow and solute transport in aquifers

This paper presents an introductory overview of recently developed stochastic theories for tackling spatial variability problems in predicting groundwater flow and solute transport. Advantages and limitations of the theories are discussed. Lastly, strategies based on the stochastic approaches to predict solute transport in aquifers are recommended.