高喜马拉雅淡色花岗岩形成的热模拟

地壳物质重熔形成高喜马拉雅淡色花岗岩需要在中地壳深度20多百万年前达到650℃以上的温度,但如何解释冷的印度板块向欧亚大陆俯冲时会产生这样的高温是一个难题,单一因素的解释往往不能解释全部地质观测事实．本文通过有限单元法定量化计算,提出多种因素的综合作用的假说,认为较合理的估计包括逆掩断层的幕式活动和间歇、摩擦剪切生热、高放射性元素含量和低热导率岩石的存在、高抬升剥蚀速度等因素．这些因素的组合才造成高喜马拉雅淡色花岗岩既能成带,但又只是断续的出现．     

高喜马拉雅淡色花岗岩形成的热模拟

地壳物质重熔形成高喜马拉雅淡色花岗岩需要在中地壳深度20多百万年前达到650℃以上的温度，但如何解释冷的印度板块向欧亚大陆俯冲时会产生这样的高温是一个难题，单一因素的解释往往不能解释全部地质观测事实．本文通过有限单元法定量化计算，提出多种因素的综合作用的假说，认为较合理的估计包括逆掩断层的幕式活动和间歇、摩擦剪切生热、高放射性元素含量和低热导率岩石的存在、高抬升剥蚀速度等因素．这些因素的组合才造成高喜马拉雅淡色花岗岩既能成带，但又只是断续的出现．     

胶东型金矿成因模型

巨型金成矿省成因是国际前沿科学问题,巨量金来源与成矿驱动机制是矿床成因研究的关键难题.胶东成矿省不足1万km²内已探明黄金资源储量～5500吨,成矿强度世界罕见,已有成矿理论难以完整解释.胶东金矿床形成于大陆内部多重构造复合域,金成矿发生于古元古代基底岩石区域变质作用约19亿年之后以及三叠纪大陆碰撞作用约1亿年之后的早白垩世120Ma左右的构造体制转换过程中.胶东金矿床和造山型、侵入岩型金矿床成矿构造环境和赋矿围岩存在差异,成矿机制也显著不同.矿床和基性岩具有紧密的时空与成因关系,胶东成矿流体起源于富集地幔,成矿期基性脉岩具低ε_Hf(t)值(-29.9～-9.1)和锇含量(0.002～0.16ppb)、重氧同位素(δ¹⁸O达8.23‰)和高初始锇同位素比值(¹⁸⁷Os/¹⁸⁸Os:0.1352～0.8858),证实其岩石圈地幔具有被大量古老地壳物质交代的特征,并且胶东西部岩石圈地幔富集程度总体高于东部,这可能是西部资源量占比远超东部的重要原因.胶东基性脉岩显示较轻的Mg同位素特...     

中国花粉R值特征

建立现代孢粉-植被定量关系是第四纪孢粉学研究的热点和难点. R值计算简单、表征直接,是校正现代孢粉和植被非线性关系的一种传统指标.本文基于中国1987～2021年发表的花粉R值数据,经筛选、校正和标准化等处理,建立了中国花粉R值数据集,包含了898个现代花粉采样点和2115条花粉R值数据,涵盖了152个花粉类群,分属于65科、113属.集成分析中国花粉R值特征,结果表明:(1)在宏观空间尺度上证实了基于样点尺度的传统认知,中国花粉R值整体变异性大,约2/3呈低代表性,如杉科(Taxodiaceae)、榛属(Corylus)、白刺属(Nitraria)、柽柳属(Tamarix)、莎草科(Cyperaceae)、禾本科(Poaceae)、豆科(Fabaceae), 1/3呈超代表性,如松属(Pinus)、云杉属(Picea)、桦木属(Betula)、麻黄属(Ephedra)、藜科(Chenopodiaceae)、蒿属(Artemisia),且低代表性较强,超代表性较弱;而栗属(Castanea)、栎属(Quercus)、蓼科(Polygonaceae)、菊科(Asteraceae)花粉代表...     

云母Rb-Sr等时线年龄原位微区LA-ICP-MS/MS测定

云母具有较高Rb/Sr比值,且Rb-Sr同位素体系在该矿物中的封闭温度高于中、低温矿产的成矿温度,因此对其进行激光原位微区Rb-Sr等时线定年是解决中、低温热液矿床成矿时代的一个重要途径.本文利用三重四极杆等离子体质谱仪(ICP-MS/MS)与193nm激光剥蚀系统(LA)联用,以SF6为反应气体,实现了同质异位素⁸⁷Rb⁺对⁸⁷Sr⁺的干扰去除,建立了激光原位微区云母Rb-Sr等时线定年方法,并分析了影响其分析准确度的因素.随后,采用该方法对K-Ar法地质年龄国家一级标准物质ZBH-25黑云母(GBW04439)进行了激光原位微区Rb-Sr等时线年龄测定.结果显示,采用NIST SRM 610为外部标准物质对样品数据进行校正时,获得的ZBH-25黑云母Rb-Sr等时线年龄相对于K-Ar年龄推荐值偏年轻约12%(RSD).由此表明,在LA-ICP-MS/MS测定云母Rb-Sr等时线年龄的过程中, NIST SRM 610标准玻璃与样品间存在明显的基体效应,导致年龄结果偏离样品推荐值.针对这一问题,本...     

基于量子计算机的海洋环流模拟算法

对海洋环流进行准确且高效的模拟是海洋科学的一个基本任务,也是一个众所周知的难点,因为它涉及求解多变量的非线性偏微分方程组.本文首次提出一种基于量子计算机的海洋环流模拟算法,可以提供指数级的量子算力加速.本文算法从海洋动力学原始方程出发,经过空间和时间离散化得到若干具有稀疏系数矩阵的线性方程组,然后利用变分量子线性方程求解器求解,从而使算法更容易在近期量子计算机上运行.这一算法属于量子-经典混合算法,为此我们发展了一套基于量子存储器和l_∞范数态层析的数据流控制方案,以完成量子计算机与经典计算机之间的数据交互.我们在MATLAB、量子虚拟机和真实量子计算机等多种计算平台上验证了算法的有效性,并讨论了测量次数和量子门噪声等因素对算法输出结果准确性的影响.结果表明,误差缓解技术可以有效提高输出解的精度.我们的算法是利用量子计算机模拟海洋环流这一挑战性问题的一个重要进展,随着量子计算科技的迅速发展,它有望在不久的将来获得应用,从而有力地促进海洋科学的发展.     

基于月球动力学演化的地质年代划分

月球地质年代表将月球的地层按形成的时间顺序划分为不同的时代,能够清晰地表达月球的整体演化历史.在地基望远镜观测和阿波罗时代早期,由于当时所用数据的空间覆盖度和分辨率有限,所建立的月球地质年代划分存在部分缺陷.在阿波罗时代结束以来的几十年里,人类对月球进行了更广泛的研究,拓展了对全月地质演化的认识和理解.基于对月球动力学演化的综合性分析,本文对当前的月球地质年代表进行了更新,特别考虑了内动力地质作用和外动力地质作用在月球早期演化中所发挥的作用.首先,根据内动力和外动力在不同时期对月球改造的效应,提出了三阶段月球动力学演化模型,并建立了三个相应的“宙/宇”级别的年代及地层单元.其次,为了更清晰地反映月球的演化历史,以南极-艾肯(SPA)盆地的形成为界线将前酒海纪划分为了岩浆洋纪和艾肯纪,其中较年轻的艾肯纪是以SPA盆地的形成为起点.本文识别了SPA盆地的溅射物地层,即达斯建造,其沉积在月球的原始月壳上,是岩浆洋结束后保存下来的最古老的外动力成因地层.更新的月球地质年代表系统综合了阿波罗时期以来月球探测和研究的成果,构建了一个描述月球演化历史的系统框架,能够更方便地开展月球地质研究,对其他类...     

1950～2020年中国沿海海平面上升

验潮站能够观测海平面长期变化,并被用于重构全球平均海平面上升.然而,中国沿海区域的海平面长期变化尚未被揭示.本文构建了一种数据同化方法,可以重构中国沿海地区1950年以来的海平面上升.该方法以全球验潮站观测为约束,同时利用气候模式输出的动态海平面和陆地水质量迁移导致的海平面指纹效应.本文重构的全球平均海平面上升与之前的研究结果接近.中国沿海地区20个验潮站的重构结果显示,1950～2020年中国沿海平均海平面上升速率为(1.95±0.33)mm a^-1,高于同时期的全球平均上升速率(1.71±0.17)mm a^-1.此外,本文还发现,中国沿海平均海平面上升速率在1980年以后是之前速率的3倍以上,速率从1950～1980年的(0.84±0.28)mm a^-1增加到1980～2020年的(3.12±0.21)mm a^-1.该发现说明中国沿海海平面存在显著的加速上升,这些结果增进了对中国沿海海平面长期变化的理解与认识.     

预测区域尺度深层土壤有机碳的方法

深层土壤有机碳(soil organic carbon, SOC)在陆地生态系统碳循环中发挥着重要作用,对其准确估算是科学认识土壤碳汇潜力的基础.然而,如何准确估算区域尺度深层土壤SOC一直是亟需解决的难题.文章建立了基于深度分布函数的预测区域深层土壤SOC的经验方法,利用全球4个不同区域的土壤剖面SOC数据库对该方法的适用性进行评估,并验证了该方法在中国黄土高原地区的应用.结果表明,在本研究所选择的深度分布函数中,负指数函数对土壤剖面SOC的模拟效果最佳,其参数C_e和k分别与0～20cm土层SOC含量和0～40cm土层SOC变化率呈现显著的线性关系.基于这些关系,建立了预测区域尺度负指数函数C_e和k的经验方程.结合负指数函数与经验方程,形成了预测区域尺度深层土壤SOC的经验模型.该模型能准确预测黄土高原样点和区域尺度500cm土层土壤SOC.因此,本研究提出的方法可简单、快速、准确地预测区域尺度深层土壤SOC.     

流域水文系统产水机制稳定性研究

流域水文系统受下垫面和气候条件的共同作用,研究其产水机制的变异规律,是合理制定变化环境下水资源管理措施、实现水安全保障的基础.本文基于Budyko框架和气候模式的历史及未来数据(CMIP6),利用Budyko全微分分解法探究了下垫面和气候变化对流域产水变异的补偿机制;基于泰勒展开式求解将Budyko各阶敏感系数(一阶到高阶项)表征流域产水率对下垫面与气候的敏感程度;根据流域外部基本特征(产水率、下垫面特征参数及干燥指数)和内部敏感系数,进一步结合联合概率函数和主成分分析方法探究了流域产水机制稳定性及其演变规律.结果表明:(1)在全球各气候分区中,下垫面和气候变化对流域产水率变异互相补偿,即气候与下垫面改变引起的产水率变化可以基本相抵,使未来流域产水率维持稳定.(2)各阶流域敏感系数表明,未来流域对下垫面敏感性增加,对气候的敏感性降低,其中干旱区流域对气候和下垫面敏感性均增加.(3)未来流域水文系统产水机制稳定性逐步下降,即从基准期(1901～1950年)到变化期(2051～2100年),稳定性系数从0.054下降到0.021;其在干旱区改变量最大,可达(-0.347±0.18),而在半...     

嫦娥五号月壤中富KREEP角砾岩及其意义

月球角砾岩为研究月壳的物质组成、冲击过程和演化历史提供了重要依据.本文重点研究了嫦娥五号(CE-5)月壤中一颗独特的表壤角砾岩(CE5C0000YJYX070GP, CE5C).它主要由玄武岩碎屑、单矿物碎屑以及冲击熔融碎屑(包括中钛玻璃、高铝玻璃和含晶体冲击熔融碎屑)组成,是CE-5样品中已知粒度最大且岩相组成最为复杂的角砾岩之一.综合借助扫描电子显微镜、TESCAN综合矿物分析仪、电子探针和激光剥蚀电感耦合等离子体质谱等技术,对其中代表性碎屑的岩相学和矿物化学进行了详细研究.分析结果表明, CE5C是一颗月海-高地混合型表壤角砾岩,含有丰富的类KREEP组分(>20vol.%).这是首次在CE-5样品中发现如此高比例的类KREEP组分. CE5C中钛玻璃原岩为富铁(24.0wt.%)、中钛(5.5wt.%)类型玄武岩. CE5C高铝玻璃的化学成分特征与KREEP组分一致.结合着陆区区域地质背景, CE5C角砾岩可能起源于P58/Em4月海单元与东部邻近高地之间的混合区域.虽然CE5C代表性有限,但其中较多的外来类KREEP组分为理解CE-5着陆区的物质来源提供了重要视角. C...     

中国碳酸盐岩地区岩溶无机碳汇格局及影响因素

碳酸盐岩化学风化碳汇降低了大气CO₂浓度上升和全球变暖的速率,然而,其碳汇通量(CCSF)的估算结果仍存在不确定性,并且气候变化和生态修复对CCSF的贡献率尚不清晰.为此,本文汇编了中国不同流域离子浓度站点数据,并采用经典热力学溶蚀模型,重新评估了中国1991～2020年CCSF的潜力和时空格局,并定量分析了温度(MAT)、降水(MAP)、蒸散发(ET)、土壤水(SM)和归一化植被指数(NDVI)等因子对CCSF的贡献率.结果发现:(1)中国CCSF为22.76t CO₂km^-2a^-1,高于全球平均水平(15.77t CO₂km^-2a^-1);总量(CCS)为4772.67×10⁴t CO₂,以252.98×10⁴km^-2的碳酸盐岩面积贡献了全球14.91%的CCS.(2)中国CCSF由东南向西北逐渐减少,其中南方岩溶区、青藏岩溶区和北方岩溶区CCS...     

深时镁同位素风化指标的验证与新解释

沉积岩碎屑组分镁同位素被广泛应用于恢复深时大陆风化强度.现有深时镁同位素风化指标基于风化残余物镁同位素随化学风化强度提高的一阶段风化模式.但是,近期实验与现代风化剖面研究发现了优先利用轻镁同位素的次生黏土矿物,挑战了深时镁同位素指标的理论基础.本研究测试了华南地区晚古生代沉积岩与中国南海沉积物样品的矿物、元素与镁同位素组成.通过对比深时镁同位素风化指标与化学蚀变指数(CIA)、镁损失量和矿物组成等风化指标来验证深时镁同位素风化指标.结果表明,在较高风化强度下镁同位素随风化强度提高而降低,而在较低风化强度下镁同位素随风化强度的提高而增加.该结果表明在不同风化强度下存在不同的镁同位素行为,以往的一阶段镁同位素风化模式无法准确追踪大陆风化强度的变化.结合深时样品数据、现今风化剖面与实验研究,我们建立了三阶段的镁含量-镁同位素风化模型.在风化初期,原生矿物的溶解与次生黏土矿物的形成造成镁含量的快速降低与镁同位素轻微提高.转换阶段土壤底层绿泥石、蛭石等富镁黏土矿物大量生成造成镁含量和镁同位素的跃升.在进阶风化阶段中,具有较重镁同位素的富镁黏土矿物优先风化,造成镁含量与镁同位素的同步降低.     

青藏高原地形重力波对东亚春季降水的重要作用

春季持续性降水是东亚夏季风爆发前的主要雨季,对东亚地区乃至全球水循环都具有深刻影响.过去研究揭示了青藏高原大尺度地形的动力和热力作用对东亚春季降水的重要影响,但青藏高原小尺度地形的作用尚不明确.本研究发现,副热带西风流受青藏高原小尺度地形扰动激发的地形重力波,同样对东亚春季降水具有重要贡献.垂直传播的地形重力波在对流层中层发生破碎,产生纬向拖曳,从而驱动跨越青藏高原的经向环流.在青藏高原南侧,经向环流上升支的动力抽吸作用使得对流层低层气压减弱,从而增强经向气压梯度力和西南季风,将更多水汽输送至下游东亚地区,最终导致春季降水的增多.该研究有助于更加全面地理解青藏高原多尺度地形对东亚春季降水的影响,为西风-季风协同作用提供了新的视角.     

三论青藏高原近地表大气氧含量影响因素及其贡献率

全球变暖背景下,预计21世纪大气中的氧气浓度将继续下降.青藏高原地处高海拔地区,缺氧形势严峻,对影响其大气氧含量因素及其量化研究仍相对匮乏.本文对北京房山站的气象要素、青藏高原生态系统产氧、碳汇和近地表大气氧含量进行了联合观测.利用季节差异和统计方法,计算了植被对近地表大气氧含量变化的相对贡献率.研究结果表明,太阳辐射、大气湿度、生态系统耗氧和产氧对近地表大气氧含量有显著影响,且其影响呈现出时空差异.基于季相差异计算结果,植被对近地表大气氧含量影响显著,贡献率达到16.7%～24.5%,已有研究低估了植被的贡献率.研究结果为认识影响近地表大气氧含量的因素提供了重要依据,突出了植被的贡献率.为更好地了解青藏高原的氧气动态,未来需在青藏高原进行更多土壤呼吸和植被光合作用的实地观测,以厘清碳储、碳汇等对近地表大气氧含量的贡献.     

海洋三维大地电磁的两步反演法:模型测试

大地电磁(MT)响应的反演已被广泛用于探索地球内部的电导率分布.因为大多数三维反演方法的反演结果都会受到初始模型的影响,因此在三维大地电磁反演中,使用好的初始模型非常重要.尽管均匀初始模型广泛用于三维反演,但一维的初始模型也是个不错的选择,因为一维模型较均匀模型可以更恰当地描绘地球内部的实际电导率变化.本研究提出了一种分两步进行的三维反演方法,该反演方法尤其适用于海洋大地电磁研究.该方法首先反演得到一维初始模型,进而依此反演三维电导率结构,并且这两次反演中均考虑了表层三维地形变化.该方法测试了含有不同地形变化(凹陷、平滑变化、槽形和方形高原地形)和不同三维棋盘类型异常(尖锐或平滑横向电导率变化)的三维电导率模型.文章比较了基于所提出方法获得初始模型和基于地形校正方法得到的初始模型,还对基于本研究方法得到的三维反演结果与基于不同的均匀初始模型得到的反演结果进行了比较.利用新引入的模型矢量参数的方向信息讨论了这些反演的结果,验证了该方法的性能和有效性.文章还揭示了新引入的模型向量参数可用于显示反演的收敛性,将有助于选择合适的反演模型.     

行星边界下的可持续食物系统评估及策略研究进展

以资源环境可持续的方式满足日益增长的食物需求是一项全球性挑战,行星边界概念框架的提出从土地利用变化、淡水利用、氮磷循环和温室气体排放等维度为可持续食物系统研究提供了科学参考.本文在全球和中国尺度,总结了行星边界框架在可持续食物系统研究中的阈值应用;基于以上阈值,评估了食物系统在耕地、淡水、氮、磷等资源利用和氮磷损失、温室气体排放等环境影响的可持续状况;梳理了改变饮食结构、改进技术和管理水平、减少食物损失浪费、优化农业生产布局等单一和综合管理策略在食物系统可持续方面的作用;最后,提出了未来可持续食物系统研究的发展方向,包括加强食物系统中各行星边界要素之间相互作用机制的研究、加强食物系统中各行星边界在国家和区域等的降尺度研究以及加强中国可持续食物系统政策制定的合理性等.本文能够为全球和中国未来可持续食物系统和农业土地利用管理提供理论基础和策略依据.     

我国各大流域复合高温干旱事件变化趋势与归因分析

受全球气候变化影响,复合高温干旱事件呈现增加趋势.通过构建表征复合高温干旱重现期变化的全微分方程,提出了一种复合高温干旱变化的归因方法,并利用1921～2020年月气温和月降水观测资料,分析中国夏季复合高温干旱事件时空分布特征及变化趋势,量化降水变化、气温变化和降水-气温相关关系变化这三种驱动因子在中国夏季复合高温干旱事件变化中的贡献.结果表明:近年来中国大部分地区夏季复合高温干旱事件呈现出明显的增加趋势,气温变化、降水变化、降水-气温相关关系变化对干旱高温复合事件变化的贡献依次递减.气温变化幅度最大、影响范围最广,在长江中上游地区以外的区域均导致复合事件增加,是最主要的驱动因子之一;降水变化是中国西部特别是新疆北部、青海和甘肃地区复合事件减少的重要原因;降水-气温关系则与其他驱动因子一起导致了在华南地区复合事件重现期的缩短和华东、西北部分地区复合事件重现期的增加.与此同时,降水和气温的均值趋势与离散程度的变化均可能导致复合高温干旱事件重现期的变化.     

基于黏弹模型的鲜水河断裂震间闭锁状态研究

断层滑动速率和闭锁状态是定量评估区域地震危险性的关键依据,然而相对于同震破裂或震后余滑来说,断层闭锁引起的地表变形较弱,利用地壳变形资料获取可靠的断层闭锁状态具有一定的挑战性.本文考虑地震周期过程中的黏弹松弛效应和其他模型参数(如块体内部应变等)的影响,融合了新增的13个近场GPS测站和区域内已有测站的观测资料,研究了鲜水河断裂的滑动速率和闭锁状态.结果显示,弹性模型会显著高估断层的闭锁深度,导致地震矩的积累速率被明显高估;同时,弹性模型得到的断层滑动速率(8.9～9.9mm a^-1)略低于黏弹模型的结果(10.5～11.6mm a^-1),与理论模拟结果一致.根据获取的鲜水河断裂的闭锁强度,沿断层走向圈定了四个潜在的凹凸体.重要的是,圈定的凹凸体与断层上历史地震的破裂分布存在明显的空间相关性,暗示凹凸体的空间特征可能控制了地震的破裂范围和震级.另外,发现2022年M6.8泸定地震仅破裂了一个潜在凹凸体的部分区域,未破裂区域积累的能量相当于一个M_w6.9的地震,潜在的地震危险性值得重点关注.     

2001～2018年瓦里全球本底站氧化亚氮浓度变化特征

对中国瓦里关站(WLG)(36°17’N, 100°54’E, 3816m asl)2001～2018年空气样品中大气N₂O摩尔分数进行了分析,得到逐年年均值、增长率和季变化特征.瓦里关站的N₂O浓度呈长期持续增长趋势,季节变化特征明显, 18年的平均年增量为(0.9±0.01)ppb a^-1(1σ)(1ppb=10^-9),接近全球平均值.平均季节变化的低值出现在6月,浓度为-(0.25±0.04)(1σ) ppb,高值出现在9月,浓度(0.13±0.07)(1σ) ppb,季节振幅为0.38ppb.季节变化主要受地面排放源、边界层高度和对-平流层传输(STE)等影响.多年增长率变化一定程度受热带平流层的准两年波动(QBO)影响.不同季节的后向轨迹聚类分析显示, N₂O主要来自中亚干旱和半干旱低排放地区的气团输送,表明瓦里关的大气N₂O浓度一定程度代表了亚洲内陆地区的背景水平.