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%。     

A biodiversity evaluation framework for restoration of aquatic macrophyte communities in shallow lakes driven by hydrological process management

Macrophyte community diversity and composition respond to ecosystem conservation and local environmental factors. In this study, we developed a multidimensional diversity framework for macrophyte communities, including the taxonomic and functional alpha and beta diversity. We used the framework to explore the relationships among water level regimes and these diversity parameters in a case study of China's Baiyangdian Lake. Analysis of indicators of hydrologic alteration divided the water level from 1959 to 2019 into four regimes (dry, <6.42 m; low, 6.42–7.23 m; medium, 7.23–8.19 m; high, >8.19 m). Alpha and beta diversity were significantly higher in the medium regime than in the low and high regimes. Redundancy analysis indicated that the maximum water depth significantly affected taxonomic alpha diversity, and total nitrogen (TN) and chemical oxygen demand (COD) concentration significantly affected functional alpha diversity, respectively. Mantel tests showed that TN, Secchi depth (SD), and water depth in the high water level regime significantly increased the total beta diversity and turnover components. TN was the main factor that increased total taxonomic beta diversity. Water level regime mainly influenced interspecific relationships by changing the TN and COD concentration. The water level should be maintained between the medium and high water level regimes to promote restoration of the macrophyte community and improve ecosystem stability. The biodiversity evaluation framework would provide a deeper insight into the hydrological process management for restoration of aquatic macrophyte communities in shallow lakes.     

An Early Cretaceous slab window beneath central Tibet,SW China: Evidence from OIB‐like alkaline gabbros in the Duolong area

The geodynamic mechanism of the late Early Cretaceous magmatic flare‐up in the collisional zone between the Lhasa and Qiangtang terranes in Tibet is controversial because of a scarcity of robust evidence. To address this problem, we report geochronological, geochemical and Hf isotopic data for the newly discovered Gufeng gabbros from the Duolong Cu–Au mineral district of the western Bangong–Nujiang Suture Zone (BNSZ). The gabbro samples, dated at 126.3 ± 1.8 Ma, show geochemical similarities to typical ocean island basalt (OIB) and have positive εHf(t) values of +3.3 to +6.9. The gabbros were generated by decompression melting of deep upwelling asthenosphere. This event is best explained by slab break‐off and the resultant development of a slab window beneath central Tibet.     

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...     

Retrieving volume FeO and TiO_2 abundances of lunar regolith with CE-2 CELMS data using BPNN method

The volume FeO and TiO_2 abundances(FTAs) of lunar regolith can be more important for understanding the geological evolution of the Moon compared to the optical and gamma-ray results. In this paper, the volume FTAs are retrieved with microwave sounder(CELMS) data from the Chang'E-2 satellite using the back propagation neural network(BPNN) method. Firstly, a three-layered BPNN network with five-dimensional input is constructed by taking nonlinearity into account. Then, the brightness temperature(TB) and surface slope are set as the inputs and the volume FTAs are set as the outputs of the BPNN network.Thereafter, the BPNN network is trained with the corresponding parameters collected from Apollo, Luna,and Surveyor missions. Finally, the volume FTAs are retrieved with the trained BPNN network using the four-channel TBderived from the CELMS data and the surface slope estimated from Lunar Orbiter Laser Altimeter(LOLA) data. The rationality of the retrieved FTAs is verified by comparing with the Clementine UV-VIS results and Lunar Prospector(LP) GRS results. The retrieved volume FTAs enable us to re-evaluate the geological features of the lunar surface. Several important results are as follows. Firstly, very-low-Ti(1.5 wt.%) basalts are the most spatially abundant, and the surfaces with TiO_2 5 wt.% constitute less than 10% of the maria. Also, two linear relationships occur between the FeO abundance(FA) and the TiO_2 abundance before and after the threshold, 16 wt.% for FA. Secondly, a new perspective on mare volcanism is derived with the volume FTAs in several important mare basins, although this conclusion should be verified with more sources of data. Thirdly, FTAs in the lunar regolith change with depth to the uppermost surface,and the change is complex over the lunar surface. Finally, the distribution of volume FTAs hints that the highlands crust is probably homogeneous, at least in terms of the microwave thermophysical parameters.     

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.