海洋沉积物重金属生态风险评价方法比较及实例验证——以莱州湾为例

利用2016年8月份莱州湾沉积物中Hg、Cd、Pb、Cu、As、Cr、Zn监测数据,分别利用单因子污染指数法、潜在生态风险指数法、地积累指数法、沉积物质量基准法、尼梅罗综合指数法、污染负荷指数法等6种沉积物重金属风险评价方法进行评价,以比较沉积物重金属生态风险评估方法的适用性与局限性。结果表明,莱州湾表层沉积物中,重金属Hg、Cd、Pb、Cu、As、Cr和Zn含量分别为0.077 4 mg/kg、0.169 mg/kg、12.1 mg/kg、11.9 mg/kg、11.6 mg/kg、31.1 mg/kg、45.3 mg/kg。不同方法在莱州湾沉积物风险评价结果上存在差异。单因素评价方法中,单因子污染指数法、潜在生态危害指数法、尼梅罗综合指数法的评价结果适中,地累积指数法评价结果偏轻,沉积物质量基准法评价结果偏重。Hg、As、Cd元素的污染评价结果均较大,为主要的污染元素;综合因素评价方法中,污染负荷指数法评价结果风险偏轻,潜在生态危害指数法评价结果适中,内梅罗综合指数法评价结果则倾向偏重。沉积物重金属含量与底栖生物的丰度、生物量的pearson相关系数表明,Hg与底栖生物的生物量呈现显著负相关。6种方法都有自身优缺点和适用范围,但潜在生态危害指数法更适合莱州湾沉积物重金属风险评价。     

粉质土海底边界层动态演化过程波浪水槽实验

The objectives of this study are carried out a series of controlled large wave flume experiments using fine-grained sediment from the Huanghe River Delta, exploring the complete sequence of sediment behavior in the bottom boundary layer(BBL) during wave-induced liquefaction. The results show that:(1) The BBL in silty seabed is exposed to a progressive wave, goes through a number of different stages including compaction before liquefaction, sediment liquefaction, and compaction after liquefaction, which determines the range and thickness of BBL.(2) With the introduction of waves, first, the sediment surface has settled by an amount S(S=1–2 cm) in the course of wave loadings with an insufficient accumulation of pore water pressure. And a thin high concentration layer formed the near-bed bottom.(3) Once the liquefaction sets in, the liquefied sediment with an ‘orbital motion' and the sub-liquefied sediment form a two-layer-sediment region. The range of BBL extends downwards and stopped at a certain depth, subsequently, develops upwards with the compaction process. Meanwhile, resuspended sediments diffuse to the upper water column.(4) During the dynamics process of the BBL beneath progressive waves, the re-suspended sediment increment ranked as sediment liquefaction erosion before liquefaction compaction after liquefaction.     

微小亚历山大藻(Alexandrium minutum)的氮磷吸收、生长、产毒动态变化

采用浮游植物计数法、高效液相色谱串联质谱法、分光光度法等分析方法,探索了微小亚历山大藻(Alexandrium minutum)在批次培养过程中氮磷吸收、产毒、生物量、p H等参数的动态变化关系。结果表明,微小亚历山大藻对磷的吸收迅速,可以将磷储存于体内,待生长使用。该藻对氮的吸收相对缓慢,环境中氮缺乏时,产毒量不再增加,说明氮对产毒有着重要的调控作用。在第二对数期中出现了碳限制环境,导致叶绿素a在碳限制条件下无法表征藻的生物量,相反,叶绿素a和生物量呈负线性关系,可能是叶绿素转化成其它含碳物质,用于生长。毒素不仅存在于细胞体内,培养液中(胞外)也含有毒素,并且胞外毒素从稳定期开始逐渐升高。胞内毒素的组成中GTX1/4占据绝对优势,GTX2/3含量相对较少。生长延缓期和第一对数期,各种毒素组成比例相对稳定,而在随后的生长期内,GTX1/4在总毒素中的占比逐渐上升,GTX2/3占比逐渐下降,表明微小亚历山大藻毒素组成会随着生长周期的变化而发生变化。     

风和径流量对长江口缺氧影响的数值模拟

受自然和人类活动的影响,海洋缺氧现象日益严重,威胁着海洋生态环境,海洋缺氧问题已经引起了人们的广泛关注。本文应用区域海洋模式并耦合生态模式,对东海的生态系统进行了数值模拟和分析研究。与观测数据比较显示,该模型能较好地模拟长江口外生态变量的分布趋势。另外本文通过设置不同敏感性实验,探讨风和径流量对长江口底层缺氧现象的影响,结果分析表明,风和径流量对长江口外缺氧区的形成有显著的影响。径流量变化虽然对长江口外缺氧区的季节变化影响并不显著,但是对缺氧区域面积却存在显著的影响。径流量增加,水体层化增强,表层叶绿素浓度增加,最终导致缺氧区域范围扩展;径流量减小,水体层化减弱,表层叶绿素浓度减小,缺氧区域范围缩小。风向和风速的改变不仅影响长江口外缺氧区的季节变化,还影响缺氧区域面积。     

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

滨海湿地作为人类活动和全球变化反应最为敏感的区域,其沉积记录可以反映出周边地区环境变化及人类活动信息。珠江口淇澳岛滨海湿地钻孔分析结果表明,在中全新世期间淇澳岛附近海域为河口湾环境,在风化层以上开始出现淤积,但在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也相应变大,有机碳主要来源于陆源输入,但在小冰期期间海源有机碳贡献略有所增大;近百年来由于受人类活动影响显著,陆源有机碳的贡献快速增加。     

基于SMAP卫星雷达资料的海冰密集度反演技术研究

SMAP是美国于2015年初发射的一颗卫星,搭载了L波段的雷达。它采用圆锥扫描方式,具有固定的入射角、较大的幅宽和千米级的分辨率,在海冰监测方面具有独特的优势。本文利用SMAP卫星雷达资料分别与德国Bremen大学海冰密集度产品和美国国家冰雪数据中心（NSIDC）海冰密集度产品建立3.125 km和25 km匹配数据集,分析了L波段雷达后向散射系数、极化比和归一化极化差与海冰密集度之间相关性,建立基于人工神经网络的海冰密集度反演算法。为了验证SMAP卫星雷达资料反演海冰密集度的精度,本文选择德国Bremen大学和美国冰雪数据中心发布的海冰密集度产品分别与SMAP海冰密集度产品进行对比分析,SMAP海冰密集度与Bremen海冰密集度的偏差为0.07、均方根误差为0.14;与NSIDC海冰密集度的偏差为0.04、均方根误差为0.18,这表明SMAP海冰密集度产品与现有业务化海冰密集度产品具有很好的一致性。     

基于神经网络模型的湛江湾水体有色溶解有机物的遥感估算

有色溶解有机物(Colored Dissolved Organic Matter, CDOM)是水体中重要的水质参数之一,是水色遥感的重要研究对象,如何构建适合特定区域的近海二类水体CDOM反演模型一直是国内外研究难点。本文利用2017年5月26～29日对南海西北部海域湛江湾20个站位采集的水样和测量的光谱资料,分析归一化遥感反射率与CDOM浓度a_g(400)的相关性,发现最大负相关系数出现在586nm处,选择580、585、590、595nm这四个波段处的归一化遥感反射率与a_g(400)建立了多元线性回归模型、BP(Back-Propagation)神经网络模型和RBF(Radial-Basis Function)神经网络模型,并与其他算法模型进行对比分析。结果发现, BP和RBF神经网络模型的平均相对误差和均方根误差均远小于多元线性回归模型和其他算法模型,神经网络模型的预测值与实测值拟合效果要优于多元线性回归模型。研究表明,神经网络模型更适合于湛江湾有色溶解有机物的遥感估算。     

Simulating the Propagation of Infrasonic Waves and Estimating the Energy of the Chelyabinsk Meteoroid Explosion Observed on February 15, 2013

Results obtained from simulating the propagation of infrasonic waves from the Chelyabinsk meteoroid explosion observed on February 15, 2013, are given. The pseudodifferential parabolic equation (PDPE) method has been used for calculations. Data on infrasonic waves recorded at the IS31 station (Aktyubinsk, Kazakhstan), located 542.7 km from the likely location of the explosion, have been analyzed. Six infrasonic arrivals (isolated clearly defined pulse signals) were recorded. It is shown that the first “fast” arrival (F) corresponds to the propagation of infrasound in a surface acoustic waveguide. The rest of the arrivals (T1–T5) are thermospheric. The agreement between the results of calculations based on the PDPE method and experimental data is satisfactory. The energy E of the explosion has been estimated using two methods. One of these methods is based on the law of conservation of the acoustic pulse I, which is a product of the wave profile area S/2 of the signal under analysis and the distance to its source E _I [kt] = 1.38 × 10^–10 (I [kg/s])^1.482. The other method is based on the relation between the energy of explosion and the dominant period T of recorded signal E _T [kt] = 1.02 × (T [s]²/σ)^3/2, where σ is the dimensionless distance determining the degree of nonlinear effects during the propagation of sound along ray trajectories. According to the data, the explosion energy E _I,T ranges from 1.87 to 32 kt TNT.     

Compressive and rarefactive double layers in non-uniform plasma with $$-nonextensive distributed electrons

Electrostatic solitary waves and double layers (DLs) formed by the coupled ion acoustic (IA) and drift waves have been investigated in non-uniform plasma using \(q\)-nonextensive distribution function for the electrons and assuming ions to be cold \(T_{i}< T_{e}\). It is found that both compressive and rarefactive nonlinear structures (solitary waves and DLs) are possible in such a system. The steeper gradients are supportive for compressive solitary (and double layers) and destructive for rarefactive ones. The \(q\)-nonextensivity parameter \(q\) and the magnitudes of gradient scale lengths of density and temperature have significant effects on the amplitude of the double layers (and double layers) as well as on the speed of these structures. This theoretical model is general which has been applied here to the \(F\)-region ionosphere for illustration.     

Generation of the Stressed State of the Lithosphere of the Earth and Mars Caused by the Reorientation of Their Figures

The computer simulation of the reorientation of the Earth and Mars lithosphere figure has been performed, which due to the dynamic redistribution of masses, allowed to reveal certain regularities of the structure-forming processes. It has been shown that the shape of the lithosphere surface has a different orientation relatively to the geoids’ (aroids) figure. This causes redistribution of masses leading to a strained state of the lithosphere as a result of endogenous and gravitational-rotational forces action in the evolutionary processes of planet’s self-development. The solution of this problem is considered on the example of lithosphere surface heights approximation by a biaxial ellipsoid with seven parameters. The acting horizontal forces in the upper shell of the planet has been calculated, introducing the concept of “evolutionary deviation of the plumb” and assuming that the tangential forces are proportional to the angle, which is defined as the angle between the direction of the plumb line in the past geological epoch and the plumb line direction at a given point. The calculated fields of tangential force vectors show good consistency with the direction of space-time displacement of Earth’s continents and tectonic plates and consistent with the results of the horizontal movements of GNSS stations. This is quite convincing evidence that under the long-term action of vortex rotationalgravitational forces, the lithospheres masses acquire the properties of creep. All this leads to the fact that interacting blocks and plates within the vortex rotational-gravitational model can be interconnected to elastic fields that creates a single planetary geodynamic field that forms the evolutionary state of the geo-environment.