Analytical investigation of the exact groundwater divide between rivers beyond the Dupuit–Forchheimer approximation

The groundwater divide is a key feature of river basins and significantly influenced by subsurface hydrological processes. For an unconfined aquifer between two parallel rivers or ditches, it has long been defined as the top of the water table based on the Dupuit–Forchheimer approximation. However, the exact groundwater divide is subject to the interface between two local flow systems transporting groundwater to rivers from the infiltration recharge. This study contributes a new analytical model for two-dimensional groundwater flow between rivers of different water levels. The flownet is delineated in the model to identify groundwater flow systems and the exact groundwater divide. Formulas with two dimensionless parameters are derived to determine the distributed hydraulic head, the top of the water table and the groundwater divide. The locations of the groundwater divide and the top of the water table are not the same. The distance between them in horizontal can reach up to 8.9% of the distance between rivers. Numerical verifications indicate that simplifications in the analytical model do not significantly cause misestimates in the location of the groundwater divide. In contrast, the Dupuit–Forchheimer approximation yields an incorrect water table shape. The new analytical model is applied to investigate groundwater divides in the Loess Plateau, China, with a Monte Carlo simulation process taking into account the uncertainties in the parameters.     

2007-2017年中国沙尘气溶胶的三维分布特征及输送过程

利用最新发布的CALIPSO产品，构建了2007-2017年中国沙尘气溶胶的三维分布，并结合HYSPLIT-4模式和再分析数据，探讨了沙尘的三维输送过程。结果表明：中国的沙尘排放源区主要是塔克拉玛干沙漠和巴丹吉林沙漠，沙尘气溶胶出现频率分别为60%和35%。塔克拉玛干沙漠排放的沙尘主要（50%~70%）停留在源地0~6 000 m高度，少部分向东输送至甘肃和内蒙古；巴丹吉林沙漠排放的沙尘则主要向东输送。中国沙尘排放量在春季最大，向东输送最强；夏季，东亚夏季风限制了沙尘向东输送；秋季，沙尘排放减弱，输送强度和夏季相当；沙尘排放量在冬季最小，输送最弱。夏季，沙尘在输送过程中可被抬升至高度5 000 m以上，春季次之，秋、冬季的沙尘主要在低层大气输送。沙尘在向东输送的过程中被抬升并和当地人为污染物混合变为污染性沙尘，华北地区污染性沙尘出现频率高达30%；输送到海洋的沙尘也会与洋面上（0~3 000 m高度）的海盐气溶胶混合，出现频率约为10%。     

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.     

Benthic foraminiferal fauna and sediment provenance of Eocene syn-rift sequences in Taiwan: implication for onset of Asian epi-continental marginal seas off China coast

Marine Geophysical Research - Early Cenozoic rift basins developed commonly on the Mesozoic basement along the SE Asia Continent. However, Eocene–Oligocene sequences were only exposed widely...     

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.     

Spatial patch structure and adaptive strategy for desert shrub of Reaumuria soongorica in arid ecosystem of the Heihe River Basin

In many arid ecosystems, vegetation frequently occurs in high-cover patches interspersed in a matrix of low plant cover. However, theoretical explanations for shrub patch pattern dynamics along climate gradients remain unclear on a large scale. This context aimed to assess the variance of the Reaumuria soongorica patch structure along the precipitation gradient and the factors that affect patch structure formation in the middle and lower Heihe River Basin (HRB). Field investigations on vegetation patterns and heterogeneity in soil properties were conducted during 2014 and 2015. The results showed that patch height, size and plant-to-patch distance were smaller in high precipitation habitats than in low precipitation sites. Climate, soil and vegetation explained 82.5% of the variance in patch structure. Spatially, R. soongorica shifted from a clumped to a random pattern on the landscape towards the MAP gradient, and heterogeneity in the surface soil properties (the ratio of biological soil crust (BSC) to bare gravels (BG)) determined the R. soongorica population distribution pattern in the middle and lower HRB. A conceptual model, which integrated water availability and plant facilitation and competition effects, was revealed that R. soongorica changed from a flexible water use strategy in high precipitation regions to a consistent water use strategy in low precipitation areas. Our study provides a comprehensive quantification of the variance in shrub patch structure along a precipitation gradient and may improve our understanding of vegetation pattern dynamics in the Gobi Desert under future climate change.

     

海水中铯-137的γ能谱测量不确定度评估

为科学评价测量结果的可靠性，对海水中137Cs γ能谱分析方法的不确定度进行了评估。根据《化学分析中不确定度的评估指南》，分析不确定度的来源并逐一对各不确定度分量进行了量化，最终计算了合成相对标准不确定度。结果表明，影响海水中137Cs γ能谱测量不确定度的因素分别为样品源137Cs峰区计数、标准源活度、取样体积、化学回收率、本底137Cs峰区计数及标准源137Cs峰区计数，其中样品源137Cs峰区计数的不确定度贡献最为突出。对于137Cs活度浓度为1.28 mBq/L的海水样品，其合成相对标准不确定度为9.78%。     

珠江口淇澳岛滨海湿地沉积环境演化及其对人类活动的响应

滨海湿地作为人类活动和全球变化反应最为敏感的区域，其沉积记录可以反映出周边地区环境变化及人类活动信息。珠江口淇澳岛滨海湿地钻孔分析结果表明，在中全新世期间淇澳岛附近海域为河口湾环境，在风化层以上开始出现淤积，但在4 200 a BP前后受极冷气候的影响，沉积物粗化；自2 500 a BP以来，沉积环境相对稳定，在小冰期期间略有变化。沉积速率计算结果显示:淇澳岛附近海域自中全新世高海面以来的平均沉积速率为0.29 cm/a，4 160~2 500 a BP、2 500 a BP-1488年、1488-1893年、1893-1986年、1990-2007年期间的平均沉积速率分别为:0.17 cm/a、0.23 cm/a、0.35 cm/a、1.37 cm/a和5.94 cm/a，沉积速率逐渐增大，反映了珠江三角洲演化过程中沉积相与沉积环境的变化；1986-1990年期间的海堤建造极大地扰动了该钻孔上部的沉积过程，在工程施工期间共沉积了厚度约112 cm的沉积层，而在海堤建成后，沉积速率也显著增大。沉积物总有机碳、总氮和C/N值的垂向分布表明，在4 160~2 500 a BP期间受海洋环境影响较大，沉积物中有机碳以海源为主，2 500 a BP以来沉积物中碳、氮含量明显增大，C/N也相应变大，有机碳主要来源于陆源输入，但在小冰期期间海源有机碳贡献略有所增大；近百年来由于受人类活动影响显著，陆源有机碳的贡献快速增加。     

Bearing capacity of offshore umbrella suction anchor foundation in silty soil under varying loading modes

Considering the current disadvantages of present offshore wind turbine foundations, a novel anchor foundation with skirt and branches is proposed, called offshore umbrella suction anchor foundation (USAF). A series of experiments and numerical simulations were performed to explore the bearing capacity of the USAF under various kinds of loading modes. The bearing characteristics and the anchor–soil interactions are described in detail for horizontal static loading, horizontal cyclic loading, and an antidrawing (pullout) test in silty soil. In the static loading test, the load–deflection of the anchor under step loading was analyzed and the normalized curve of the load–deflection was obtained to determine the ultimate horizontal bearing capacity of the anchor under normal working conditions. Under horizontal cyclic loading, the relationship between the plastic cumulative deformation and cyclic number was determined. In addition, the responses of USAF were investigated for a low wave frequency and storm surges. In the drawing test, it was found that a “segmentation phenomenon” occurred during the test. Moreover, a method to identify the maximum antidrawing load of USAF was provided based on dynamic mechanics. The numerical results show that the use of anchor branches and skirt can enhance the bearing performance of USAF to a certain degree. However, the anchor branch has a slight positive influence on the bearing performance improvement. The USAF is not only similar to a stiff short pile, but a rotation occurs. The failure envelope under composite loading (V-M) was obtained and the changes associated with changes in the aspect ratio of the internal compartment were clarified.