中国花粉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)花粉代表...     

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.该发现说明中国沿海海平面存在显著的加速上升,这些结果增进了对中国沿海海平面长期变化的理解与认识.     

胶东型金矿成因模型

巨型金成矿省成因是国际前沿科学问题,巨量金来源与成矿驱动机制是矿床成因研究的关键难题.胶东成矿省不足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同位素特...     

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

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

热带海表面温度及中纬度大气环流对青藏高原9月降水异常的共同影响

由于青藏高原独特的地理位置,其降水受到热带和中纬度异常环流系统的共同影响.利用观测、再分析和CMIP6模式数据,本文揭示了印度洋偶极子(IOD)和丝绸之路遥相关(SRP)对青藏高原9月降水年际变率的单独和联合影响.当IOD处于正位相,热带印度洋异常纬向海温梯度激发一个Gill响应,在印度次大陆和孟加拉湾低层产生异常反气旋.其北侧的异常西风在青藏高原南坡形成浅槽,相关的异常西南风向青藏高原东南部输送水汽,导致东南部降水增加.同时,北大西洋的西风急流扰动能够激发一个SRP遥相关波列,在印度北部产生一个倾斜的斜压结构,即中西亚对流层高层为异常反气旋,而印度北部对流层低层为异常气旋.异常气旋东侧的西南风向东南部青藏高原输送大量水汽,并造成降水增加. IOD和SRP的联合影响可以解释将近52%的青藏高原降水异常,远超过IOD(19%)和SRP(27%)的单独贡献.本文强调了考虑热带和中纬度驱动因子共同作用的必要性,这为更准确地模拟和预测青藏高原降水提供了线索.     

喜马拉雅淡色花岗岩铷铯成矿作用

西藏喜马拉雅淡色花岗岩带中新发现铯沸石锂云母钠长花岗岩,其中铯沸石是首次在中国喜马拉雅带内发现的铯独立矿物.该岩石位于吉隆岩体的北部,与吉隆岩体电气石二云母花岗岩共生.其主要稀有金属矿物包括锂云母、锂辉石、铯沸石、锡石和细晶石等.矿物化学分析结果显示,云母主要为含锂白云母-锂云母, Li₂O和Rb₂O平均含量为4.07wt.%和0.76wt.%;锂辉石的Li₂O平均含量为7.95wt.%;铯沸石的Cs₂O和Rb₂O平均含量为34wt.%和0.16wt.%.全岩Li₂O、Rb₂O和Cs₂O含量分别为0.49～1.19wt.%, 0.12～0.24wt.%和0.69～2.33wt.%,均达到了工业品位.锡石和独居石定年结果表明,该岩石形成于19～18Ma,略晚于吉隆电气石二云母花岗岩(22～20Ma).结合上述两类岩石之间的脉动接触关系、形成时代及全岩和矿物化学成分,本文认为它们应形成于同一岩浆体系,铯沸石是岩浆发生...     

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

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

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

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

云母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标准玻璃与样品间存在明显的基体效应,导致年龄结果偏离样品推荐值.针对这一问题,本...     

中国328个城市的植被覆盖度长期变化特征及其驱动因子

城市绿地有助于促进城市居民健康,增强城市生态系统生物多样性,减轻环境污染,缓解城市热岛效应.然而城市植被盖度异质性高,量化难度大.本文利用30m分辨率的陆地卫星数据,研究了1990～2022年间我国328个城市植被覆盖度的变化.结果显示, 1990～2005年,由于城市内建筑容积率和不透水表面的增加,全国平均城市植被覆盖度由0.38降至0.35.此后,城市植被覆盖度开始逆转,到2022年增长至0.45,这种增长趋势在新建城区最明显. 2005年以前,平均城市植被覆盖度的降低主要是低植被覆盖区域的扩张造成的;此后,城市植被盖度的增加则表现为高植被覆盖区域的扩张.城市植被覆盖度变化趋势的转折与2004年以后国家陆续出台政策要求增加城市绿地覆盖度有关,城市国内生产总值与城市植被覆盖度变化呈现最高相关性.对于大城市和中等规模城市,影响植被覆盖度的前三个因子分别为国内生产总值、城市人口和气温.但对于干旱/半干旱地区的城市,植被覆盖度变化对气候因素(如降水)的敏感度更高.尽管近年来中国城市的总体植被覆盖度显著增加,但在小城市和大城市老城区的城市绿地仍有改善空间.     

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

深层土壤有机碳(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),而在半...     

青藏高原湖泊纹层孢粉记录的过去2000年雨季时长和雨季降水时空变化

青藏高原雨季时长及雨季降水量的长尺度定量重建,是揭示高原西风、季风时空演变规律及其生态环境效应的关键,其中准确判断高原雨季的起止时间是研究难点.本文通过寻找降水拐点的方法确定了高原不同地点的雨季起、止时间,并计算了雨季时长和雨季降水量,建立了表土孢粉组合与雨季起止时间、雨季时长和雨季降水量的转换函数,重建了具有年纹层沉积特征的高原北部库赛湖区和中部江错湖区过去2000年来的雨季起止时间、雨季时长和雨季降水量变化.结果表明,过去2000以来,库赛湖区雨季降水序列记录了5次高降水时段:AD580～680、AD1000～1100、AD1200～1450、AD1550～1780、AD1920至今,其普遍对应雨季天数长的时段;江错湖区雨季降水序列记录了4次高降水时段:AD80～500、AD800～950、AD1250～1450、AD1780至今,其与雨季天数长的时段在AD1000以前对应良好,在AD1000后关系不明确.在空间上,高原雨季降水呈现出“南北同湿”、“南北同旱”、“南湿北旱”和“南旱北湿”四种模式,其中“南北同湿”可能与异常强盛的夏季风有关,“南北同旱”与弱西风和弱季风有关,“南湿北...     

嫦娥五号月壤中富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...     

全球变暖背景下中国森林春季木质部物候提前的模拟证据

植被物候对当前和未来全球气候变暖的响应,直接影响森林生态系统初级生产力的变化.然而,由于树木木质部物候的监测数据十分稀少,学界较少关注树木木质部物候的变化,这阻碍了对森林大空间尺度碳收支的评估和预测.本研究使用基于生理生态过程的树轮生长模型(Vaganov-Shashkin模型)模拟1962～2016年中国森林站点春季树木木质部物候(木质部生长季开始时间)的时空变化.使用实测树轮宽度序列(70条)对树轮生长模型模拟结果进行校正,校正后的模拟结果与已有独立的微树芯监测结果显著相关.模拟结果发现,全球变暖背景下春季木质部物候在1962～2016年期间显著提前,其提前速率在20世纪90年代后显著增快,平均每年提前0.25天.生长季前的日均温与大多数站点(71%)的生长季开始时间显著相关,表明其可能是中国树轮站点春季木质部物候的主要气候驱动因子.春季增暖的气温允许树木生长的热量需求更早达到,因此导致春季木质部物候显著提前.生长季前日均温每升高1℃,春季物候提前6～7天,这将有益于偏冷湿环境下树木的径向生长.模拟的春季树木木质部物候与遥感物候间存在显著正相关,表明全球增暖背景下树木初生生长和次生...     

季节补偿效应导致2015/2016厄尔尼诺期间北半球未发生明显的碳汇降低

厄尔尼诺事件通常在北半球导致广泛的植被生长下降和碳汇降低.然而,基于遥感观测、全球陆地生态系统模型模拟及大气CO₂反演多套数据,本研究表明在2015/2016厄尔尼诺成熟阶段,北半球(主要为热带之外的区域)植被普遍持续增绿,陆地碳汇未见减弱.研究结果发现,春季植被生长增强对随后的夏季和秋季植被生长有明显的补偿效应,这种补偿效应维持了夏季植被增绿,并导致2015年春季、夏季及2016年春季陆地碳汇略有增加,特别是在生长季前期(前一年的11月至当年的3月)水分供应增加的北半球热带之外的区域.两套独立数据集结果表明,相对于厄尔尼诺事件发生前5年均值, 2015年春季和夏季陆地碳汇分别平均增加23.34%和0.63%, 2016年春季平均增加6.82%.同时,相较于1997/1998厄尔尼诺事件,更强的季节补偿效应有效地缓解了2015/2016年厄尔尼诺事件对植被生长的不利影响.本研究进一步明确,生长季前期的供水对后续植被生长的遗留效应能够持续约6个月.该研究结果凸显了季节补偿效应在应对偶发性极端厄尔尼诺事件时对植被生长和陆地碳汇的调节作用.     

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

碳酸盐岩化学风化碳汇降低了大气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...     

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

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

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

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

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

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