海洋表层沉积物中金属镉的来源解析

利用2015年8月份对长江口及其邻近海域表层沉积物的监测数据,基于主成分分析/绝对主成分分数(PCA/APCS)受体模型定量解析了重金属元素镉(Cd)的可能来源,并结合地统计学插值了Cd的源贡献量的空间分布状况,结果表明沉积物中镉污染主要存在3个可能来源,源头及贡献率分别为工业污染(18.8%)、陆地径流输入(66.0%)、生物活动等自然因素(13.6%),并且各个源头贡献量具有不同的空间分布状况,其中工业污染的高值区主要集中在靠近陆地区域,陆地径流输入的分布呈现由近岸向外海逐减降低的特征,生物活动等自然因素的高值区主要集中在远离陆地的外海区域。     

毛乌素沙区1961-2016年降水特征

基于毛乌素沙区10个气象站1961-2016年观测资料，应用Mann-Kendall方法和t检验法对各气象站年降水量进行了突变检验，借助小波分析讨论了各气象站年降水量的周期特征，根据降水量等值线划分结果对整个研究区分区分析了年、季、月和日尺度上的降水变化特征，并在两个时段上分析了季节性降水的差异。结果表明：毛乌素沙区年降水量空间特征差异明显，东部亚区呈上升趋势，中西部亚区呈下降趋势，但变化趋势不显著且无突变发生；降水年内分配不均，干湿季分明，降水集中在5-9月，夏秋季降水占全年降水比重大，季、月和日尺度降水量存在梯度递减变化；年降水量的年际变化过程存在多重时间尺度的自相似结构；近26年的冬春季降水增加显著，但降水波动幅度小于前30年。     

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.     

北京主城区伪基站时空规律分析

随着公众移动通信的快速发展,伪基站的泛滥不仅破坏正常电信秩序,危害公共安全,而且严重损害群众财产权益,侵犯公民个人隐私,已成为社会一大公害。如何从垃圾短信大数据中挖掘出伪基站活动的时空规律,寻找有效的防控方案,从源头上进行打击和治理成为管理部门和研究者共同关注的焦点。本文基于北京市垃圾短信数据,利用非负矩阵分解的方法分析伪基站的时空分布规律;并利用TF-IDF构建垃圾短信分类模型,对垃圾短信进行分类,结合土地利用数据,分析伪基站在发送不同类型垃圾短信时的时空分布规律。结果显示：北京市垃圾短信多分布于路网和中心城区;白天垃圾短信数量远远多于晚上;垃圾短信的分布随时间的推移沿着路网逐渐向内收缩;发送不同类型垃圾短信的伪基站的时空分布具有一定的差异;通过非负矩阵分解得到的结果,与垃圾短信分类后得到的结果有很好的匹配。研究表明,非负矩阵分解具有实现上的简便性、分解形式和分解结果上的可解释性等优点,可以有针对性的为有关部门建言打击伪基站的有效方案,对于伪基站违法行为的治理具有一定的意义。     

Wetland hydroperiod classification in the western prairies using multitemporal synthetic aperture radar

Wetlands represent one of the world's most biodiverse and threatened ecosystem types and were diminished globally by about two‐thirds in the 20th century. There is continuing decline in wetland quantity and function due to infilling and other human activities. In addition, with climate change, warmer temperatures and changes in precipitation and evapotranspiration are reducing wetland surface and groundwater supplies, further altering wetland hydrology and vegetation. There is a need to automate inventory and monitoring of wetlands, and as a study system, we investigated the Shepard Slough wetlands complex, which includes numerous wetlands in urban, suburban, and agricultural zones in the prairie pothole region of southern Alberta, Canada. Here, wetlands are generally confined to depressions in the undulating terrain, challenging wetlands inventory and monitoring. This study applied threshold and frequency analysis routines for high‐resolution, single‐polarization (HH) RADARSAT‐2, synthetic aperture radar mapping. This enabled a growing season surface water extent hyroperiod‐based wetland classification, which can support water and wetland resource monitoring. This 3‐year study demonstrated synthetic aperture radar‐derived multitemporal open‐water masks provided an effective index of wetland permanence class, with overall accuracies of 89% to 95% compared with optical validation data, and RMSE between 0.2 and 0.7 m between model and field validation data. This allowed for characterizing the distribution and dynamics of 4 marsh wetlands hydroperiod classes, temporary, seasonal, semipermanent, and permanent, and mapping of the sequential vegetation bands that included emergent, obligate wetland, facultative wetland, and upland plant communities. Hydroperiod variation and surface water extent were found to be influenced by short‐term rainfall events in both wet and dry years. Seasonal hydroperiods in wetlands were particularly variable if there was a decrease in the temporary or semipermanent hydroperiod classes. In years with extreme rain events, the temporary wetlands especially increased relative to longer lasting wetlands (84% in 2015 with significant rainfall events, compared with 42% otherwise).     

Revisiting the existence of an effective stress for wet granular soils with micromechanics

A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.     

非等间距GM(1,1)模型的总体最小二乘算法及其病态问题

在非等间距GM(1,1)模型中，系数矩阵中有无误差的常数项和有误差的随机项，并且系数矩阵与观测向量误差同源，即系数矩阵与观测向量中有相同的元素存在，这些相同元素应该有相同的改正数，为此本文推导了一种适合非等间距GM(1,1)模型求解的总体最小二乘算法。同时，考虑到非等间距GM(1,1)模型中存在病态问题时影响总体最小二乘计算结果的稳定性，提出对系数矩阵常数列乘以某一常数的方法，以改善病态问题。     

Assessment of alternative adsorption models and global sensitivity analysis to characterize hexavalent chromium loss from soil to surface runoff

We investigate our ability to assess transfer of hexavalent chromium, Cr(VI), from the soil to surface runoff by considering the effect of coupling diverse adsorption models with a two‐layer solute transfer model. Our analyses are grounded on a set of two experiments associated with soils characterized by diverse particle size distributions. Our study is motivated by the observation that Cr(VI) is receiving much attention for the assessment of environmental risks due to its high solubility, mobility, and toxicological significance. Adsorption of Cr(VI) is considered to be at equilibrium in the mixing layer under our experimental conditions. Four adsorption models, that is, the Langmuir, Freundlich, Temkin, and linear models, constitute our set of alternative (competing) mathematical formulations. Experimental results reveal that the soil samples characterized by the finest grain sizes are associated with the highest release of Cr(VI) to runoff. We compare the relative abilities of the four models to interpret experimental results through maximum likelihood model calibration and four model identification criteria (i.e., the Akaike information criteria [AIC and AIC_C] and the Bayesian and Kashyap information criteria). Our study results enable us to rank the tested models on the basis of a set of posterior weights assigned to each of them. A classical variance‐based global sensitivity analysis is then performed to assess the relative importance of the uncertain parameters associated with each of the models considered, within subregions of the parameter space. In this context, the modelling strategy resulting from coupling the Langmuir isotherm with a two‐layer solute transfer model is then evaluated as the most skilful for the overall interpretation of both sets of experiments. Our results document that (a) the depth of the mixing layer is the most influential factor for all models tested, with the exception of the Freundlich isotherm, and (b) the total sensitivity of the adsorption parameters varies in time, with a trend to increase as time progresses for all of the models. These results suggest that adsorption has a significant effect on the uncertainty associated with the release of Cr(VI) from the soil to the surface runoff component.     

新疆库车县盐水沟北铜矿床地质特征及成因探讨

盐水沟北铜矿区位于新疆库车县北部,其大地构造位置处于塔吉克-塔里木陆块区、塔里木陆块、塔里木北缘隆起、库车前陆盆地,区内构造活动强烈。铜矿床赋存于新近纪中新统康村组第一岩性段,矿体呈层状,严格受地层控制。区内共圈定10个铜矿体,矿石矿物主要为孔雀石、赤铜矿、水胆矾、自然铜等,矿体内无夹石,顶、底板围岩均为褐红色及浅红色粉质泥岩。该矿床为砂岩型铜矿床。     

Parallelization of the flow-path network model using a particle-set strategy

High-performance simulation of flow dynamics remains a major challenge in the use of physical-based, fully distributed hydrologic models. Parallel computing has been widely used to overcome efficiency limitation by partitioning a basin into sub-basins and executing calculations among multiple processors. However, existing partition-based parallelization strategies are still hampered by the dependency between inter-connected sub-basins. This study proposed a particle-set strategy to parallelize the flow-path network (FPN) model for achieving higher performance in the simulation of flow dynamics. The FPN model replaced the hydrological calculations on sub-basins with the movements of water packages along the upstream and downstream flow paths. Unlike previous partition-based task decomposition approaches, the proposed particle-set strategy decomposes the computational workload by randomly allocating runoff particles to concurrent computing processors. Simulation experiments of the flow routing process were undertaken to validate the developed particle-set FPN model. The outcomes of hourly outlet discharges were compared with field gauged records, and up to 128 computing processors were tested to explore its speedup capability in parallel computing. The experimental results showed that the proposed framework can achieve similar prediction accuracy and parallel efficiency to that of the Triangulated Irregular Network (TIN)-based Real-Time Integrated Basin Simulator (tRIBS).