南海南部次表层叶绿素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质量浓度受水温的影响明显。     

基于ROMS和4DVAR的沿轨与网格化SSH数据同化效果评价

Remote sensing products are significant in the data assimilation of an ocean model. Considering the resolution and space coverage of different remote sensing data, two types of sea surface height(SSH) product are employed in the assimilation, including the gridded products from AVISO and the original along-track observations used in the generation. To explore their impact on the assimilation results, an experiment focus on the South China Sea(SCS) is conducted based on the Regional Ocean Modeling System(ROMS) and the four-dimensional variational data assimilation(4 DVAR) technology. The comparison with EN4 data set and Argo profile indicates that, the along-track SSH assimilation result presents to be more accurate than the gridded SSH assimilation, because some noises may have been introduced in the merging process. Moreover, the mesoscale eddy detection capability of the assimilation results is analyzed by a vector geometry–based algorithm. It is verified that, the assimilation of the gridded SSH shows superiority in describing the eddy's characteristics, since the complete structure of the ocean surface has been reconstructed by the original data merging.     

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.     

不同排列方式和间隙比下双圆柱流致振动特性研究

采用基于ALE的流固耦合分析方法,将圆柱体振动模型简化为平面双自由度质量-弹簧-阻尼系统,进行了计算参数无关性分析。对不同θ角和间隙比T下的双圆柱流致振动进行数值模拟,并与单圆柱振动计算结果进行比较,得到了不同θ角和间隙比T下中心圆柱的升阻力系数、流场分布和运动轨迹变化情况。通过分析发现小间隙比下由于两圆柱间的相互干扰,改变了中心圆柱的尾涡结构,导致了升力和阻力系数随不同θ角而发生变化;小间隙比下由于两圆柱间的相互干扰影响,均使中心圆柱的轨迹在不同θ角下发生了较大的变化,顺流向的振动频率与单圆柱振动时相比出现了多倍频,导致其运动轨迹发生了变化。     

圆形数值波浪水池斜向和多向入射波浪模拟研究

为了研究斜向入射波浪,基于三维不可压缩两相流模型,开发了一套圆形数值波浪水池数值模型。在圆形波浪水池中,通过源项造波法成功生成了任意入射方向的波浪,并且利用人工摩擦项模拟阻尼区以数值耗散反射波浪。模型基于嵌入式多块网格体系,采用FVM法(finite volume method)离散Navier-Stokes方程,VOF法(volume of fluid)追踪自由水面。试验结果表明,斜向入射波浪的模拟结果与理论值基本一致,圆形波浪水池在模拟斜向入射波浪时,有效区域的面积较传统波浪水池显著增大,而且有效区域受波浪入射角度的影响也较小。同时,通过叠加多列斜向入射波浪,模拟出了多向交叉波列,并通过与理论结果对比,发现其具有较高的精度。     

内蒙古鲁新井田水文地质条件研究

矿产资源开发导致了地下水失衡,地下水失衡又给矿产开发造成了极大的安全隐患,近年来,矿产资源开发与地下水环境保护之间的矛盾愈发凸显。通过对内蒙古鲁新井田典型的水文地质条件进行分析研究,分析采矿导致地下水失衡机理,深入剖析矿井开采充水条件及矿山开发对地下水环境的影响,合理提出了促进矿产开发与地下水保护相互协调的对策建议,为实现"采矿保水"协调统一提供了基础地质依据。为类似地区矿产开发过程中遵循自然规律,趋利避害,保障生产安全,保护地下水环境安全,实现资源绿色开发有较好的指导借鉴作用。     

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.     

海岸破波带内悬沙空间分布试验研究和解析表达式

关于海岸破波带内悬沙浓度水平和垂向分布的研究对于计算海岸输沙率和地形演变具有重要意义。本研究进行了规则波、波群和不规则波三种波浪情况破波带内悬沙浓度的水平和垂向分布的试验测量。试验在大尺度波浪水槽进行,接近实际海岸波况尺度。给出了破波带内多断面悬沙垂向分布的细致测量结果,并以此为基础给出了预报实际海岸破波带内悬沙浓度水平和垂向分布解析表达式,讨论了形成这些分布的物理原因和不同波况、不同破波带区域对分布的影响。     

300 BC-1900 AD无定河流域城镇时空格局演变

通过分析和整理从战国中晚期（约300 BC）至清末（约1900 AD）无定河流域历代县级及以上城址的位置、兴废年代数据,结合行政区划沿革、经济社会发展、政权更迭等资料,分析流域城镇格局的时空演变过程。研究表明：① 受气候周期性波动影响,无定河流域城镇的兴起与衰废具有明显的周期性特征,城镇几何中心的移动轨迹具有明显的“西北—东南”向潮汐性运动特征;② 城址存续年限普遍较短,具有明显的阶段性特征,流域内曾存在过4个阶段性中心城镇,中心城镇移动的方向和过程与城镇几何中心的潮汐性移动过程趋势一致;③ 流域城镇空间格局可分为3种类型,即秦汉与隋唐时期的沿河流谷地分布,宋、明两代的沿边境线与长城分布,以及元、清两代集中于流域下游分布。流域城镇未来的发展布局应重视气候变化对城镇分布的长期影响,关注区域环境的脆弱性,合理安排城镇体系发展规模与布局。     

基于SODA资料的南海表层风能输入的空间分布与长期趋势研究

海面风不仅是驱动上层海洋运动的主要动力, 其能量也是维持海洋表层流动的主要机械能来源。为了分析南海表层流风能输入的变化, 用SODA(Simple Ocean Data Assimilation)(1901—2010)资料估算了风向南海表层流(表层地转流+表层非地转流)的能量输入。结果表明, 风向南海表层流、表层地转流和表层非地转流输入的能量总体均呈减少趋势, 110年间分别减小了约56%、65%和49%。导致风能输入减小的最主要因素是风应力的减弱(减小了35%)。由于南海受季风系统的控制, 风向表层流及其各成分输入的能量呈现出显著的季节性变化。冬季风能输入最强, 高值区位于南海西部及北部区域, 呈一个显著的“回力镖”状结构。这些结果对深入认识南海环流具有理论意义。