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.     

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.     

南海南部次表层叶绿素a质量浓度最大值及其影响因子*

文章分析了2013年南海南部4个季节航次的叶绿素a (Chl a)调查数据, 结果显示: 150m以浅水柱Chl a质量浓度均值分别为早春0.14mg•m^-3、初夏0.12mg•m^-3、初秋0.18mg•m^-3、初冬0.16mg•m^-3。早春和初夏偏低的原因与早春风速小, 初夏水温高, 不利于水体的垂直混合, 限制了深层海水中丰富的营养盐向上层水体补充有关。4个季节中海水次表层Chl a质量浓度最大值层(SCML)均出现在50m和75m, 这两个水层的Chl a质量浓度差异小, 季节变化不大, 平均值变化范围分别为0.24~0.26mg•m^-3和0.22~0.26mg•m^-3。受混合层深度和温跃层上界深度的共同影响, 50m水层Chl a质量浓度主要受制于深层富营养盐海水的向上补充, 75m水层Chl a质量浓度受水温的影响明显。     

微量元素分析在判别沉积环境中的应用——以柴达木盆地西部锶矿区富锶地层为例

沉积岩中的微量元素对沉积环境水介质的变化较为敏感,是研究沉积物沉积时古气候、古环境的有效手段。选取了对沉积环境介质较为敏感的锶(Sr)、钡(Ba)、铜(Cu)、铀(U)、钒(V)、镍(Ni)、钍(Th)等元素作为判别指标,对柴达木盆地西部上新统油砂山组、狮子沟组和更新统七个泉组三套富锶岩层的古气候和古沉积环境特征进行了研究。结果表明,柴达木盆地西部从油砂山组到七个泉组沉积时总体上为干旱的气候环境;湖盆水体的盐度较高,表现为咸水环境,且从油砂山组到七个泉组随着地层变新,盐度有逐渐增大的趋势;油砂山组、狮子沟组和七个泉组沉积时湖盆水体均表现为中等分层、湖水循环较为顺畅的氧化环境。     

渤、黄海生态环境模拟的参数敏感度空间差异分析

随着海洋生态系统模型的发展，生态变量增多，众多生物过程参数量值的确定成为制约生态环境模拟的瓶颈问题，生态系统结构区域性要求模型中的生态参数具有区域差异。为探究不同海区的关键参数及参数敏感度的空间差异，本研究在渤、黄海建立了ROMS-CoSiNE物理–生物耦合的高分辨率生态系统模型，并对13种生态参数的敏感度空间分布进行分析。结果表明:南黄海中部与渤海及近岸海域的敏感度差异较大。渤海敏感度最大的参数为决定光合速率的浮游植物P-I曲线初始斜率，其次为浮游动物捕食半饱和常数和浮游动物最大捕食率。而南黄海中部敏感度最大的参数为浮游动物最大捕食率，其次为浮游植物死亡率和浮游植物P-I曲线初始斜率。结合敏感度分布及浮游植物生物量收支得出，渤海水体透明度较南黄海偏低、浮游植物生长光限制较强，是引起浮游植物P-I曲线初始斜率敏感度在渤海高于黄海的主要原因。浮游动物最大捕食率及浮游植物死亡率的敏感度空间差异，受渤、黄海浮游植物生物量差异的影响，与生态系统中的高度非线性特征有关。     

南海神狐海域沉积物自生黄铁矿和石膏S同位素特征及对甲烷渗漏强度的示踪

The northern slope of the South China Sea is a gas-hydrate-bearing region related to a high deposition rate of organic-rich sediments co-occurring with intense methanogenesis in subseafloor environments.Anaerobic oxidation of methane(AOM) coupled with bacterial sulfate reduction results in the precipitation of solid phase minerals in seepage sediment,including pyrite and gypsum.Abundant aggregates of pyrites and gypsums are observed between the depth of 667 and 850 cm below the seafloor(cmbsf) in the entire core sediment of HS328 from the northern South China Sea.Most pyrites are tubes consisting of framboidal cores and outer crusts.Gypsum aggregates occur as rosettes and spheroids consisting of plates.Some of them grow over pyrite,indicating that gypsum precipitation postdates pyrite formation.The sulfur isotopic values(δ~(34) S) of pyrite vary greatly(from –46.6‰ to –12.3‰ V-CDT) and increase with depth.Thus,the pyrite in the shallow sediments resulted from organoclastic sulfate reduction(OSR) and is influenced by AOM with depth.The relative high abundance and δ~(34) S values of pyrite in sediments at depths from 580 to 810 cmbsf indicate that this interval is the location of a paleo-sulfate methane transition zone(SMTZ).The sulfur isotopic composition of gypsum(from–25‰ to –20.7‰) is much lower than that of the seawater sulfate,indicating the existence of a 34 S-depletion source of sulfur species that most likely are products of the oxidation of pyrites formed in OSR.Pyrite oxidation is controlled by ambient electron acceptors such as MnO_2,iron(Ⅲ) and oxygen driven by the SMTZ location shift to great depths.The δ~(34) S values of gypsum at greater depth are lower than those of the associated pyrite,revealing downward diffusion of 34 S-depleted sulfate from the mixture of oxidation of pyrite derived by OSR and the seawater sulfate.These sulfates also lead to an increase of calcium ions from the dissolution of calcium carbonate mineral,which will be favor to the formation of gypsum.Overall,the mineralogy and sulfur isotopic composition of the pyrite and gypsum suggest variable redox conditions caused by reduced seepage intensities,and the pyrite and gypsum can be a recorder of the intensity evolution of methane seepage.     

东海蛾螺属二新记录种(腹足纲,蛾螺科)

在整理近年来东海采集的蛾螺科Buccinidae标本时,鉴定出属于蛾螺属Buccinum Linnaeus,1758的两个中国新记录种:(1)白肋蛾螺Buccinum leucostoma Lischke,1872,采自于东海300—400m深的泥沙质海底;(2)古式蛾螺Buccinum koshikinum Okutani,1988采自于东海400m深的泥沙质海底。文中分别对这两个新记录种的形态特征、生活习性等进行了描述,并与相似种进行了分类学讨论。此外,文中还列出了中国海已报道的其他蛾螺属种类。     

Genetic variability of the Lessepsian migrant mussel Brachidontes pharaonis (Bivalvia: Mytilidae) in Tunisia

The present study used two mitochondrial markers (16S rRNA and COI) to assess the genetic diversity of a newly founded Lessepsian migrant mussel, Brachidontes pharaonis, in Tunisian waters. The species appears to be restricted to only one population in Rades Harbour, in the northern part of the country. Phylogenetic analyses revealed the monophyly of B. pharaonis in Tunisia. Both molecular markers revealed high genetic variability of the B. pharaonis population. Haplotype networks and demographic analyses confirmed the recent expansion events within this population. Multiple human-mediated introduction events involving several founder populations and intensive population growth rates are probably the main causes of the high polymorphism observed within this invasive mollusc.     

全球气候变化下南海诸岛保护优先区识别分析

全球变化下,珊瑚礁保护区是保护生物多样性、增强珊瑚礁对气候变暖抵抗力的有效方式,而维持珊瑚礁弹性是其核心内容。针对珊瑚礁最具有威胁性的热压力因子,基于南海1982—2009年卫星观测海表面温度（SST）数据和CMIP5加拿大地球系统模式CanESM2模型预估的2006—2100年南海SST数据构建热压力强度模型,从维持珊瑚礁弹性的角度识别IPCC RCP 4.5和RCP 8.5情景下南海诸岛保护优先区。结果表明：RCP 4.5和RCP 8.5情景下13%左右的南海诸岛珊瑚礁识别为保护优先区。根据热压力强度与珊瑚抵抗力及避难所关系,西沙群岛七连屿和晋卿岛近年观测与未来预估的热压力强度均比较低,在保障其服务功能的基础上建议实施完全保护;东沙群岛东沙环礁和中沙环礁排洪滩近年观测急性热压力强度较高但未来预估热压力强度较低,建议实施50%禁止利用保护;中沙群岛黄岩岛近年观测和未来预估的急性热压力强度均比较低,建议实施50%多用途保护。南沙群岛有14%左右的珊瑚礁识别为保护优先区,根据其热压力强度可实施30%~100%禁止利用保护或30%~50%多用途保护。RCP 4.5和RCP 8.5情景下的南海诸岛保护优先区及保护对策,可为维持珊瑚礁生态弹性及应对全球气候变化提供重要的参考价值。     

Lidar-aided analysis of boreal forest backscatter at Ku band

This paper analyzes the backscatter of the microwave signal in a boreal forest environment based on a Ku -band airborne Frequency-Modulated Continuous Waveform (FMCW) profiling radar—Tomoradar. We selected a half-managed boreal forest in the southern part of Finland for a field test. By decomposing the waveform collected by the Tomoradar, the vertical canopy structure was achieved. Based on the amplitude of the waveform, the Backscattered Energy Ratio of Canopy-to-Total (BERCT) was calculated. Meanwhile, the canopy fraction was derived from the corresponding point cloud recorded by a Velodyne VLP-16 LiDAR mounted on the same platform. Lidar-derived canopy fraction was obtained by counting the number of the first/ the strongest returns versus the total amount of returns. Qualitative and quantitative analysis of radar-derived BERCT on lidar-derived canopy fraction and canopy height are investigated. A fitted model is derived to describe the Ku-band microwave backscatter in the boreal forest to numerically analyze the proportion contributed by four factors: lidar-derived canopy fraction, radar-derived canopy height, the radar-derived distance between trees and radar sensor and other factors, from co-polarization Tomoradar measurements. The Root Mean Squared Error (RMSE) of the proposed model was 0.0958, and the coefficient of determination R2 was 0.912. The fitted model reveals that the correlation coefficient between radar-derived BERCT and lidar-derived canopy fraction is 0.84, which illustrates that lidar surface reflection explains the majority of the profiling /waveform radar response. Thus, vertical canopy structure derived from lidar can be used for the benefit of radar analysis.