山东地区新生代不同类型幔源样品中稀有气体同位素组成与成因

稀有气体是地球流体源区信息的有效示踪剂.利用热熔融法对山东地区幔源巨晶、包体全岩及其单矿物以及寄主玄武岩等幔源样品中稀有气体的测试表明,不同类型样品具有明显不同的稀有气体同位素组成,特别是He同位素比值差异明显.巨晶和包体单矿物中的n(3He)/ n(4He)值分别为0.87～11.31Ra、0.911 ～10.11Ra,均值分别为6.1Ra和3.73Ra,主体介于地幔一大气一地壳源区之间;包体全岩中n(3He)/ n(4He)值0.19 ～ 2.30Ra(均值为1.05 Ra),且均低于相应单矿物中该比值;而玄武岩中n(3He)/ n(4He)值全为10-8数量级,与地壳均值无异.各类幔源样品中n(4Ar)/n(36 Ar)值无明显差异,介于296.4～864.3,高于大气但远低于典型地幔值.所获有限数据中,不同类型样品中n(2Ne)/n(2Ne)值均高于大气,具幔源特征,其中包体全岩及其单矿物n(2Ne)／n(22Ne)值分别为10.19 ～12.42和10.55～11.80,无明显差异;玄武岩中该值介于11.07～13.10,总体高于包体全岩及其单矿物.上述比值特征的差异反映了不同类型幔源样品各自不同的成因和演化特征:巨晶和包体单矿物中的轻稀有气体同位素组成主要反映了由于古板块俯冲所导致的岩浆源区地幔一大气／地壳的昆合特征,个别样品中可能存在宇宙成因3He;除单矿物中所体现的源区混合特征外,包体全岩还存在放射成因4He;而玄武岩在喷出地表后丢失大量地幔稀有气体信息,其同位素组成主要体现了放射成因4He的影响.     

基于Sentinel-1/2动态耦合移栽期特征的水稻种植模式识别

及时准确地获取水稻种植模式变化对于有效防控“非粮化”和完成“双碳”目标具有重要意义。现有研究多基于固定时间窗口挖掘水稻生长期特征,且多使用的是单一卫星影像数据,难以应用于大范围水稻制图。本文通过Sentinel-1/2数据构建动态窗口提取移栽期光学/雷达特征,利用其耦合关系实现水稻种植模式制图。将该算法应用于湖南和江西两省水稻制图。基于1402个地面参考点位对水稻提取结果进行验证,总体精度达92.80%;在县域尺度上,湖南和江西两省水稻制图面积与农业统计数据也具有高度一致性,R²达0.85以上。相比于用固定窗口进行水稻特征提取,该方法具有较强的鲁棒性和迁移能力,为实现更大范围作物制图提取提供新的思路和参考依据。2018—2021年江西省水稻制图结果表明,水稻总种植面积减少9.47%,约3460 km²,水稻种植强度从1.62下降至1.49;在种植模式上,“双改中”趋势明显,双季稻种植面积锐减21.61%,其中约84%改种中稻。     

一种适用于有限差分模式的负载平衡区域分解方法

分布式内存并行处理在数值天气预报等超大规模科学计算中已经得到了广泛的应用。中尺度模式由于分辨率高 ,计算量大 ,需使用更多的处理机进行并行运算。另一方面 ,由于复杂的物理过程的采用 ,增加了不同天气的计算量的不平衡。但是 ,目前所广泛使用的并行处理方法在处理机数量较多时不能很好地均衡计算负载 ,引起并行计算效率的降低。本文提出了一种新的非规则区域分解负载分配方法。并与已有的负载分配方法进行了分析试验对比 ,该方法能更有效地平衡负载 ,取得更好的加速效果     

对偶Kriging插值方法在气象资料分析中的应用

非对偶Kriging方法是一种局部插值, 因此对每个插值点都要重新求解Kriging线性系统从而非常耗时。该文介绍了通过等价变换将泛Kriging方法转换成对偶Kriging方法, 对偶Kriging方法是一种全局插值, 其Kriging线性系统不依赖插值点, 因此仅需一次求解Kriging线性系统即可计算所有插值点的值, 从而极大提高了计算效率。数值试验的实际计算表明, 对偶Kriging方法不仅计算精度完全与泛Kriging方法一致, 整体效果相当于或优于GrADS绘图软件的Cressman方法; 而且对偶Kriging方法的计算效率远高于泛Kriging方法。最后, 该文通过统计与拟合得到了降水的4类半变异函数模型的表达式, 并通过敏感试验研究4类半变异函数模型对降水分析精度的影响。     

雅安地区石材开发中对资源及环境的保护

介绍了四川省雅安地区储量丰富、质地名贵的石材资源;论述了在石材开发中由于管理不善,采矿、加工工艺落后,存在资源破坏严重和地质环境恶化等问题;提出了合理开采石材资源和进行地质环境保护的建议     

近50年中国霾年代际特征及气象成因

根据1961-2013年全国745个国家基准站的长期观测资料,分析中国霾日数年代际变化特征及可能的气象成因。结果表明:近50年来,中国霾天气主要集中在东部从华南到华北的大部分地区,霾日数呈增加趋势。秋冬两季是霾天气发生最频繁、变化最明显的两个季节。中国东部淮河以南地区秋冬两季霾日数在2000年前呈增加趋势,其后增加趋势变得较为平缓,20世纪90年代前霾日数与近地面风速呈显著负相关关系,90年代后则与大气相对湿度呈显著负相关关系,随着90年代前近地面风速减小和90年代后大气相对湿度降低,该区域霾日数表现出明显的增加趋势。中国东部从淮河到华北大部分地区秋冬两季霾日数1980年后增加趋势变得不明显,这可能与该区域近地面风速和大气相对湿度的变化趋势较为平缓有关。     

北大西洋ODP1049C孔Aptian-Albian期高频旋回大洋红层的成因:矿物学证据

北大西洋ODP171B航次1049C孔Aptian-Albian期沉积以出现大洋红层与灰色、白色沉积物高频旋回为特征。为了探讨大洋红层的成因,本文进行了矿物学、地球化学、沉积学等方面的研究。漫反射光谱、磁化率和活性铁数据表明,赤铁矿、针铁矿的出现是导致样品由白色向红色转变的矿物学原因。棕色和橙色样品中出现赤铁矿和针铁矿的特征峰,FeR/FeT平均值分别为0.23和0.24,明显区别于其他颜色的样品。磁化率与红层成良好的正相关关系,说明铁氧化物矿物含量的变化是导致磁化率变化的原因。X射线衍射结果表明,不论颜色如何,样品中均含伊利石、高岭石、蒙脱石、绿泥石等粘土矿物,其分布与样品颜色没有直接关系,很可能反映当时物源区气候稳定。ODP1049C孔岩芯沉积物出现橙色、棕色、白色、灰色等颜色的高频变化,颜色过渡接触界线清晰,说明导致红色变化的赤铁矿和针铁矿是沉积期低温氧化的产物。棕色样品中CaCO3含量最低,推测红层所对应的氧化条件是由于较低的有机质堆积速率造成的。     

中国海岸风沙沉积记录的环境变化

海岸是陆、海、气相互作用的地带。海岸风沙沉积是这一特殊动力环境的产物,是研究海岸环境演变及海平面变化的良好信息载体。中国海岸主要存在3种风沙沉积,分别为“老红砂”、沙丘岩及海岸沙丘。本文通过比较已报道的风沙堆积的物质组成、地层变化等,进一步总结了海岸风沙沉积的特征;选择已开展绝对测年的沉积剖面,利用概率密度函数分析了风沙沉积年代的分布特征,考察了中国海岸风沙活动历史,结合其他气候记录,探讨了不同地质历史时期海岸风沙堆积的关键影响因素。结果表明：“老红砂”沉积主要记录了冰期-间冰期尺度的风沙活动,在120 ka BP前后、73—55 ka BP等时期,风沙活动主要与海平面下降时陆架提供的丰富沙源和强盛的冬季风有关;而在105 ka BP、80 ka BP前后,风沙活动与高海平面时期丰富的沙源或季风气候的季节性增强有关;55 ka BP以来风沙活动强度降低更多地反映了沉积记录保存环境的变化,末次冰期海岸风沙沉积大多分布在现代海面之下,并不代表实际的海岸风沙活动减弱。相比之下,海岸沙丘沉积所记录的风沙活动主要发生在近3 000 a,可能与中国海岸冬季风的增强有关。     

基于夜间灯光影像的云南省建成区演变研究

夜间灯光影像可衡量城市进程.建成区与植被覆盖度具有一定的负相关性.本文用基于植被校正的城市夜间灯光指数(Vegetation Adjusted NTL Urban Index,VANUI)提取云南城市建成区.先对两种影像进行相互校正,再与归一化植被指数合成基于植被校正的城市夜间灯光指数,最后基于指数采用分水岭分割算法提取城市建成区进行时空演变分析.研究表明:VANUI指数可以对建成区边缘进行较高精度的识别;滇中建成区的扩张速度和强度在全省保持在较高水平;主要城市的建成区紧凑度都在减小;全省建成区重心的偏移呈现东南—东南—东南—东北—西南—东北方向,建成区重心偏移位置变化不大.     

GRAPES全球格点模式的并行计算负载平衡策略

随着高性能计算机技术的发展和应用,并行计算已成为保证数值天气预报模式业务运行时效的关键技术之一.目前高性能计算机计算能力已达到每秒千万亿次浮点计算.系统中处理器数量也早已达十万甚至更多,如此巨大的计算资源对应用软件系统的设计也提出了挑战.数值天气预报软件系统要充分利用高性能计算机提供的计算资源,必须依靠并行计算方法,这包括适合计算问题的可扩展并行算法的设计、合适的数据分配方案以及良好的任务负载平衡方案.作为中国新一代数值天气预报格点模式,GRAPES(Global and Regional Assimilation and PrEdiction System)设计的最终目标是一个科研/业务通用,区域/全球通用模式.作为一个格点模式,GRAPES的并行计算具有与欧洲中期数值顶报研究中心谱模式并行计算不同的特点,GRAPES的并行计算采用了经典的水平网格数据划分.但对于全球的GRAPES模式,由于采用拉格朗日差分方案,模式极地及附近区域格点与格点之间距离的减小.使得模式并行计算在采用简单的经纬网格划分方式实现时,必须考虑极地区域并行计算跨越多个处理器时导致的频繁通讯解决途径.本研究提出了利用消息传递组通讯实现全球格点模式并行计算的一种方法,其核心思想是将极点附近一定区域内的处理器按纬向划归不同的处理器组.文中还给出了该实现方法的任务分配算法,提出了改进的任务分配负载平衡方案.在中国气象局高性能计算机IBM-cluster1600上的测试表明,算法具有较好的可扩展性,其负载平衡方案改善了计算的绝对墙钟时间,使并行计算效率提高10%以上.模式的准业务运行结果表明计算墙钟时间基本可以满足数值预报业务的实时性要求.