基于改进土壤冻融水循环的Biome-BGC模型估算青藏高原草地NPP

Biome-BGC模型被广泛用于估算植被净初级生产力（Net Primary Productivity, NPP）,但是该模型未考虑冻土区土壤冻融水循环过程对植被生长的影响。本文基于Biome-BGC模型,改进冻土区活动层土壤冻融水循环,估算了2000—2018年青藏高原高寒草地NPP。通过比较原模型和改进后的模型,并对NPP模拟结果的时空特征进行了分析,结果表明：① 增加冻融循环提高了NPP估算精度,青藏高原草地NPP均值由114.68 gC/(m²·a)提高到128.02 gC/(m²·a)。② 原模型和改进后NPP的空间分布差异较大,时间变化趋势差异不明显。③ 青藏高原草地NPP总量为253.83 TgC/a,呈东南向西北递减的空间格局,年均增速为0.21gC/(m²·a)（P=0.023）,显著增加的占17.85%,主要分布在羌塘高寒草原地带的大部分地区和藏南山地灌木草原地带的西部。④ 该冻融水循环改进方法简单可靠,具有在其他多年冻土区推广的价值。     

青藏高原维管植物物种丰富度分布的情景模拟

如何充分利用离散的观测数据,通过对维管植物物种分布丰富度及其与生境因子之间的相互作用和影响机理的定量分析,实现维管植物物种丰富度的空间分布及其情景模拟,是目前生物多样性研究前沿和核心内容之一。针对这一问题,在实现青藏高原37个国家自然保护区的维管植物物种数量收集和边界数据矢量化的基础上,分别进行维管植物物种数量与土地覆盖类型、环境因子和景观生态指数等三大类生境因子之间的相关关系的定量计算和对比分析,筛选和确定最佳相关分析方程,进而构建青藏高原维管植物物种丰富度的空间模拟分析模型。该模型中,维管植物物种丰富度与生境因子之间的复相关系数为0.94,模型验证结果表明,青藏高原的维管植物物种的平均丰富度为496.79种/100 km²,其空间分布格局整体上呈东南向西北逐渐减少趋势;另外,除柴达木盆地荒漠区域以外,维管植物物种的空间分布随海拔的升高而减少。基于CMIP5 RCP 2.6、RCP 4.5和RCP 8.5三种气候情景模拟获得的青藏高原维管植物物种丰富度未来情景结果显示,在T0-T4（2010-2100）时段内,青藏高原维管植物物种丰富度整体将呈减少趋势。RCP 8.5情景下青藏高原维管植物物种丰富度的变化幅度最大,而RCP 2.6情景下的维管植物物种丰富度的变化幅度最小。研究表明,本文构建的模型能够对青藏高原维管植物物种丰富度的空间分布格局及其未来情景进行模拟分析,模拟结果可为青藏高原生物多样性及其对气候变化响应的综合评估和情景模拟提供方法和技术支持。     

青藏高原沙漠化土地空间分布及区划

利用野外调查数据、遥感影像和已有研究成果,构建了一套适用于青藏高原沙漠化土地的分类分级指标体系及遥感解译标志。以此为基础,选取目视解译法监测青藏高原沙漠化土地的空间分布特征。结果表明：2015年青藏高原沙漠化土地面积392 913km²,占高原土地总面积的15.1%,主要包括沙质沙漠化土地、砾质沙漠化土地和风蚀残丘3种类型。沙漠化土地以中度和轻度沙漠化土地为主,重度和极重度沙漠化土地面积仅占沙漠化土地总面积的12.2%。空间上,沙漠化土地集中分布在高原的北部和西部地区,其他地区零散分布。自东南向西北,沙漠化土地面积逐渐增大,沙漠化程度不断加重。以沙漠化土地空间分布数据（面积、类型、程度、空间特征和驱动因素）为基础,结合气候、地貌、第四纪沉积物和人类活动等数据,将青藏高原沙漠化区划分为雅鲁藏布江半干旱高山宽谷沙漠化区、藏北青南高寒高原面沙漠化区、柴达木干旱盆地沙漠化区、黄河上游半干旱河流盆地沙漠化区和“三江”流域湿润半湿润高山沙漠化区。     

1982~2013年青藏高原高寒草地覆盖变化及与气候之间的关系

利用GIMMS NDVI数据和地面气象站台观测数据,对青藏高原1982~2013年高寒草地覆盖时空变化及其对气象因素的响应进行研究,结果表明：青藏高原高寒草地生长季NDVI表现为从东南到西北逐渐减少的趋势,近32 a来,整个高原草地生长季NDVI呈上升趋势,增加速率为0.000 3/a (p<0.05);高寒草地生长季NDVI年际变化具有空间异质性,整体为增加趋势,呈增加趋势的面积约占研究区域面积的75.3%,其中显著增加的占26.0% （p<0.05）,类型主要为分布在青藏高原东北部地区的高寒草甸;比例为4.7%,草地类型主要为高寒草原,主要分布在高原西部地区;基于生态地理分区的分析显示,青藏高原草地与降水、温度的相关关系具有明显的空间差异,高寒草地生长季NDVI均值与降水呈显著正相关,对降水的滞后效应显著;高原东北部温度较高,热量条件较好,降水为高寒草地生长季NDVI变化的主导因子;东中部地区降水充沛,温度则为高寒草地生长的制约因子;南部地区降水和温度都较适宜,均与高寒草地生长季NDVI相关性显著(p< 0.05),共同作用于草地的生长;中部和西部地区,气候因子与高寒草地生长季NDVI关系均不显著。     

基于人居环境特征的青藏高原“无人区”空间界定

青藏高原独特的高寒环境与自然条件在一定程度上限制了人口的自然分布与有序发展,形成了中国面积大、分布广的“无人区”（UPAs）。然而,当前有关“无人区”面积、分布、特征与区域差异等研究尚无定论。客观、准确界定“无人区”的空间范围,对开展青藏高原资源环境承载力评价、国家公园与生态安全屏障建设等具有重要意义。基于青藏高原居民点分布信息,据其地形、气候、生态、土地利用等要素特征,本文综合表征了居民点的自然—生态—土地利用耦合关系,率定了居民点分布上限的各要素阈值,通过多要素空间叠加构建了“无人区”评价综合模型,并以居民点分布的自然极限、生态（含氧量）下限、土地利用规律为关键阈值界定了青藏高原“无人区”空间范围并分析了其地理分布特征。研究表明：① 以居民点分布累计比例< 0.1%计,确定“无人区”的地形阈值为海拔> 5665 m、相对高差> 2402 m、地形起伏度> 8.59,气候阈值为相对湿度< 76.2%、温湿指数< 33或 > 71。② 根据居民点分布及人体对含氧量耐受情况,确定“无人区”的生态阈值为气压< 500 hpa、大气含氧量< 40%。③ 青藏高原严格“无人区”面积达1912 km²,其中新疆699 km²、四川413 km²、西藏331 km²、青海291 km²、甘肃178 km²。空间上呈零星分散状,多分布在四川贡嘎山、珠穆朗玛峰附近等极高山地区、可可西里东部—罗布泊地区;以及少部分分布在青海柴达木盆地。     

Thickness estimation of the Longbasaba Glacier: methods and application

A total of 71,177 glaciers exist on the Qinghai-Tibet Plateau, according to the Randolph Glacier Inventory (RGI 6.0). Despite their large number, glacier ice thickness data are relatively scarce. This study utilizes digital elevation model data and ground-penetrating radar thickness measurements to estimate the distribution and variation of ice thickness of the Longbasaba Glacier using Glacier bed Topography (GlabTop), a full-width expansion model, and the Huss and Farinotti (HF) model. Results show that the average absolute deviations of GlabTop, the full-width expansion model, and the HF model are 9.8, 15.5, and 10.9 m, respectively, indicating that GlabTop performs the best in simulating glacier thickness distribution. During 1980-2015, the Longbasaba Glacier thinned by an average of 7.9±1.3 m or 0.23±0.04 m/a, and its ice volume shrunk by 0.28±0.04 km³ with an average reduction rate of 0.0081±0.0001 km³/a. In the investigation period, the area and volume of Longbasaba Lake expanded at rates of 0.12±0.01 km²/a and 0.0132±0.0018 km³/a, respectively. This proglacial lake could potentially extend up to 5,000 m from the lake dam.     

青藏高原地区城乡虚拟水贸易格局与影响因素

青藏高原是亚洲水塔,其水资源与水生态环境保护意义重大。从虚拟水视角,研究青藏高原与外部的水资源贸易关系和影响因素,有助于理解该地区的水资源问题、制定虚拟水贸易策略、优化区域城乡水资源配置、保障亚洲水塔功能。依托2012年中国区域间投入产出表成果,本文测算了青藏高原与中国其他区域之间的虚拟水贸易关系,建立了中国区域城镇与农村地区的虚拟水贸易网络,采用对数平均迪氏指数模型分析了青藏高原对其他区域虚拟水贸易不平衡的影响因素。结果表明,青藏高原向中国其他区域净输出虚拟水2.25亿m³,其中向西南、华北、华中、华东、华南等5个区域净输出虚拟水,从西北和东北2个区域净输入虚拟水。城乡之间虚拟水贸易联系非常紧密,农村地区生产水足迹较高,而城镇地区由于人口密度较高、消费水平较高,是虚拟水最终消费的热点区域,青藏高原农村地区的虚拟水贸易量大于城镇地区的虚拟水贸易量。青藏高原贸易输出结构以农产品为主导,虚拟水净输出12.7亿m³;青藏高原与其他区域贸易存在逆差,贸易量因素导致虚拟水净输入8.6亿m³;用水效率在青藏高原与不同区域虚拟水贸易中的正负效应不一,总体带来青藏高原虚拟水净输入1.8亿m³。未来,应重点通过灌溉节水减少农业水足迹,引导城镇居民向低水足迹生活方式转变,鼓励内地为青藏高原提供物质与技术支援,实行水资源生态补偿政策,以保护青藏高原水资源,促进区域水资源可持续利用。     

清代青藏高原东北部河湟谷地林草地覆盖变化

河湟谷地是青藏高原东北缘典型的农牧交错区,清代以来耕地的扩张导致林草地覆盖发生明显变化。本文在现代植被图的基础上,选取土壤、地形因素,并依据历史文献数据,重建河湟谷地潜在林地草地格局,在此基础上结合清代耕地变化的重建结果,推算出清代河湟谷地林草地覆盖的变化状况。结果显示：①清代耕地扩张之前,其林、草地分布与现今各类植被类型的空间分布格局基本一致,林地分布范围比现代略大,灌木林地在空间上连续性更强,草地分布区域更广;②估算出河湟谷地潜在林地、灌木林地、草地面积分别约为0.28×10⁴、0.93×10⁴、2.1618×10⁴ km²,由于耕地开垦,至清代末期,河湟谷地草地、灌木林地、林地面积分别累计减少5180.41、1330.35、441.31 km²,其中草地被垦殖占用的面积最大,程度最深,减少的区域主要集中在湟水谷地中游的乐都盆地、西宁盆地以及黄河谷地的尖扎盆地、化隆盆地等;③清代河湟谷地中人类垦殖原始覆盖类型的差异性不仅受自然环境的限制,同样受到社会政策因素的影响。     

基于InVEST模型的陕北黄土高原水源涵养功能时空变化

以陕北黄土高原为研究区,基于InVEST水源涵养功能评价模块,定量评价退耕还林还草工程背景下土地利用/覆被变化对研究区水源涵养的影响,在此基础上开展水源涵养空间分区。结果表明：① 2000-2010年,陕北黄土高原草地、灌丛和林地的面积分别增加了3204 km²、285.3 km²和122.7 km²,城镇面积增加了450.4 km²;农田、荒漠、湿地的面积分别面积减少了3984.5 km²、72.7 km²和5.2 km²。② 2000-2010年,陕北黄土高原水源涵养量整体以减少为主,中部减少最为显著,减少量在25 m³/hm²~40 m³/hm²,局部区域在40 m³/hm²以上;其他大部分区域均有0~25 m³/hm²不等的减少。③ 陕北黄土高原水源涵养功能高度重要区和极重要区的总面积为32255.1 km²,所占比例为40.5%。④ 通过陕北黄土高原水源涵养功能评价和重要性分区为生态系统的科学管理提供参考。     

基于MODIS温度的青藏高原多年冻土活动层厚度变化研究

基于MODIS温度数据,采用TTOP模型和Stefan公式模拟了青藏高原地区的冻土分布并计算了活动层厚度,并与地面观测结果进行了对比。结果表明：2003—2019年青藏高原多年冻土面积为1.01×10⁶ km²;多年冻土活动层厚度区域平均值为1.79 m, 活动层厚度区域平均的变化率为3.67 cm/10a,且草甸地区的变化率明显大于草原地区,5100~5300 m高程带的活动层厚度变化速率最大。     

近年来的黄土高原耕地时空变化与口粮安全耕地数量分析

基于遥感调查数据集定量分析了1990—2015年中国黄土高原地区耕地的时空变化特征和口粮绝对安全最小耕地保障面积的数量变化。结果表明：黄土高原耕地面积从1990年的192 529.65 km²至2015年的182 688.50 km²,净减少了9 841.14 km²,幅度达5.11%,其中2000—2010年的减幅最大,净减少8 483.00 km²;较大的耕地动态变化图斑主要分布于中部和西部区域,细碎的变化图斑广泛分布;耕地地类转出面积（31 875.82 km²）大于转入面积（21 815.25 km²）,耕地面积的增加主要由草地和林地转化而来,主要分布在灌溉农业区和东南部平原区,减少的耕地主要转化为草地和林地,主要分布在中部沟壑区的雨养农业区。此外,该时期耕地转化为建筑用地和交通用地等人工表面的面积逐渐增加,主要分布在东南部低海拔平原地区;黄土高原口粮绝对安全所需最小耕地保障面积呈明显减少特征（从1990年的70 913.37 km²下降到2015年的33 981.64 km²）,占该区耕地总面积比例呈明显缩减态势（从1990年的36.83%缩减到2015年的18.60%）,目前耕地总量的净减少未对口粮绝对安全的耕地保障数量造成大的影响。     

青藏高原高寒草地净初级生产力(NPP)时空分异

基于1982-2009 年间的遥感数据和野外台站生态实测数据,利用遥感生产力模型(CASA模型) 估算青藏高原高寒草地植被净初级生产力(NPP),分别从地带属性(自然地带、海拔高程、经纬度)、流域、行政区域(县级) 等方面对其时空变化过程进行分析,阐述了1982 年以来青藏高原高寒草地植被NPP的时空格局与变化特征。结果表明:① 青藏高原高寒草地NPP多年均值的空间分布表现为由东南向西北逐渐递减;1982-2009 年间,青藏高原高寒草地的年均总NPP为177.2×10¹² gC·yr^-1,单位面积年均植被NPP为120.8 gC·m^-2yr^-1;② 研究时段内,青藏高原高寒草地年均NPP 在112.6~129.9 gC·m^-2yr^-1 间,呈波动上升的趋势,增幅为13.3%;NPP 增加的草地占草地总面积的32.56%、减少的占5.55%;③ 青藏高原多数自然地带内的NPP呈增加趋势,仅阿里山地半荒漠、荒漠地带NPP呈轻微减低趋势,其中高寒灌丛草甸地带和草原地带的NPP增长幅度明显大于高寒荒漠地带;年均NPP增加面积比随着海拔升高呈现"升高—稳定—降低"的特点,而降低面积比则呈现"降低—稳定—升高"的特征;④ 各主要流域草地年均植被NPP均呈现增长趋势,其中黄河流域增长趋势显著且增幅最大。植被NPP和盖度及生长季时空变化显示,青藏高原高寒草地生态系统健康状况总体改善局部恶化。     

青藏高原植被垂直带与气候因子的空间关系

集成了青藏高原气候区149个山地植被垂直带数据,利用国家基本气象台站自建站以来到2001年的地面观测日气象数据,计算了地面的温暖(WI)、寒冷(CI)、湿润(MI)、吉良龙夫(Kira)干湿指数、干燥度(Idm)等水热指数,运用GIS的空间分析模块,模拟了青藏高原水热条件的空间分布形势,探讨山地植被垂直带谱分布规律与制约因子的定量指标.结果表明:在高原的东北部、西北边缘,以荒漠、荒漠草原、山地森林、山地草原、灌丛、草甸为组合的半干旱、干旱结构向高原腹地以高寒草原、高山草甸、荒漠带组合的高寒干旱带谱结构的变化;东南、南部边缘,以温暖湿润为特征的以森林带为优势带谱组合结构逐渐向寒冷的高原中心变化;高原的地势效应,致使的水热形势旱现从中央向边缘变化的趋势是致使青藏高原植被垂直带谱分布的重要原因.     

基于无人机遥感的高寒草原沙化模型及等级划分

高寒草原是青藏高原草地生态系统的主要组成,在防风固沙、野生动物保育等方面具有重要作用。近年来,在全球气候变化和人为干扰加剧的背景下,高寒草原沙化加剧,基于时空尺度监测范围及程度是防治高寒草地沙化的前提。以青海三江源区玛多县的高寒草原为研究区,结合大疆“精灵3”和“经纬M100”旋翼无人机和地面调查,探讨基于无人机遥测的植被指数在草地沙化调查方面的适宜性,以此为基础制定了高寒草原沙化模型及等级划分标准。结果显示：（1）通过对VDVI（Visible-Band Difference Vegetation Index）、ENDVI（Enhanced Normalized Difference Vegetation Index）和NGRDI（Normalized Green-Red Difference Index）指数与草地沙化指数GDI（Grassland Desertification Index）的相关分析,选取出高寒草地沙化研究最优植被指数为VDVI（R=0.9055）;（2）GDI与VDVI的关系模型为VDVI=0.3024GDI²-0.0335GDI+0.0119（R²=0.9326）。模型相对误差为1.779%（RMSE=0.165,R²=0.7447）,拟合精度较高;（3）基于无人机遥感植被指数的聚类分析,将研究区高寒草原沙化划分为5个等级,即无明显沙化（VDVI>0.2247）、轻度沙化（0.1493 < VDVI < 0.2246）、中度沙化（0.0924 < VDVI < 0.1492）、重度沙化（0.0692 < VDVI < 0.0923）和极度沙化（VDVI<0.0692）。     

青藏高原典型区微塑料分布特征及来源分析

微塑料作为一种环境中普遍存在的新型污染物正在引起越来越多的关注。本文基于2019年5—7月在青藏高原地区采集的53个水样和52个土壤样品微塑料检出结果,分析了研究区水土环境中微塑料赋存情况,并讨论了其潜在来源。研究区域的微塑料样点分别位于两大区域：青藏高原南部地区和祁连山地区。青藏高原南部地区水体中微塑料浓度显著高于祁连山地区,而土壤中微塑料丰度在两个区域差别不显著（p>0.05）。水土样品中微塑料浓度范围分别为0~1916.66个/m³和0~260个/kg,平均值为(438.21±454.94)个/m³和(48.35±36.25)个/kg。水体中的微塑料以透明纤维为主,主要聚合物为聚丙烯;土壤微塑料中最多的是透明薄膜,主要聚合物为聚乙烯。来源分析表明洗衣废水和旅游业带来的生活垃圾可能是水体微塑料的主要来源,而土壤微塑料潜在主要来源是农业塑料覆膜的使用。研究结果揭示了青藏高原水土环境中微塑料的赋存形态和空间分布特征,可为高寒地区微塑料迁移转化机制研究提供数据支撑。     

青藏高原草地生态系统碳收支研究进展

陆地生态系统碳收支仍然是当前全球气候变化研究的重要内容,青藏高原作为全球气候变化的敏感区,使青藏高原草地生态系统在区域碳收支平衡中占有主导地位,但研究方法等不同使得碳收支估算结果存在很大的不确定性。气候变暖在一定程度上提高了高寒草地生态系统的植被初级生产力和生物量,由此补偿了气候变暖导致的土壤有机碳分解释量,使青藏高原草地植被仍然发挥着碳汇的功能。而人类放牧活动对草地生态系统的影响较为复杂。因此,如何区分气候变化和人类活动对生态系统的影响机制,定量评价未来气候变化和人类活动影响下,青藏高原生态系统碳源/汇格局的可能变化,是一个非常重要的研究方向,也是一个极大的挑战。     

青藏高原高寒草原土壤活性有机碳的分布特征

利用36个样点数据,分析了青藏高原高寒草原土壤活性有机碳(SAOC)分布特征.结果表明:(1)在水平方向上,SAOC含量呈现出东南高、西北低的总体态势和斑块状交错分布的格局,高值区主要集中在藏北高原腹地和喜马拉雅北麓湖盆区,不同草地型和自然地带SAOC含量差异显著;(2)在垂直方向上,不同草地型和自然地带0～40cm剖面SAOC含量分布状况,均可分为由高到低型、由低到高型和低-高-低型3个类型,表土层(0～10cm)与底土层(30～40cm)SAOC含量差异显著;(3)基于回归模型的标准系数法,分析了气候因子对高寒草原SAOC含量的影响程度,指出降水对高寒草原SAOC含量的贡献大于气温.     

基于GIS的青藏高原史前交通路线与分区分析

青藏高原自然环境恶劣,但并非是生命的禁区。确切证据表明,末次冰消期人类已扩张至高原地区,并通过长期实践形成了较为稳定的交通路线格局。通过构建自然因子模拟—遗址分布校正方法,利用GIS得到青藏高原史前交通路线模拟结果,再结合考古证据对模拟路线进行印证。依据模拟结果将高原分为以下4个区：东北区、东南区、西南区以及西北无人区。其中东北区路线遗址点密度高,落入路线内的遗址点数量占总数的88.56%。路线形成于旧石器时期,其方向为东西向,密切联系东北区的内部,并向外沟通黄土高原与北方地区,其形成与发展为彩陶、粟黍、小麦的传播以及丝绸之路的形成奠定了基础。东南区路线为南北走向,北与黄河上游地区连接,南与四川盆地、云贵高原沟通,加强了南北方地区的文化交流,是民族融合与交流的大走廊。西南区位于高原的腹心地带,模拟的路线体现了人类对青藏高原主体的征服,其线路总长度4602.32 km为3个区域内最长,是人类向高原进军的第2条重要通道。西北无人区主要为高寒荒漠区,其自然环境极端恶劣不适宜人类生存,故无路线分布。     

黄河流域沙漠化空间格局与成因

黄河流域地势西高东低,自西向东有青藏高原、内蒙古高原、黄土高原和黄淮海平原4个地貌单元,总面积为79.6万km²。黄河流域是中国重要的生态屏障和重要经济带,黄河流域生态保护和高质量发展已上升为国家战略。沙漠化是中国北方干旱、半干旱以及部分半湿润地区主要的土地退化形式,沙漠化对黄河流域,尤其是流域中、上游地区的影响较大。为了全面掌握黄河流域的沙漠化土地空间分布特征,本研究以Landsat遥感影像为数据源,通过地理信息系统（GIS）技术,获得了黄河流域2010年的沙漠化土地分布数据。结果表明：黄河流域内沙漠化土地面积为128 667 km²,占流域总面积的16.2%;黄河上游的沙漠化土地面积最大,然后依次是黄河中游、黄河源区、黄河下游,沙漠化土地面积分别为89 341、21 426、17 894、7 km²,分别占全流域沙漠化土地总面积的69.4%、16.7%、13.9%、0.01%。黄河流域的沙漠化土地绝大部分分布于内蒙古,其沙漠化土地面积为91 398 km²,占全流域沙漠化土地面积的71.0%;其次是青海,沙漠化土地面积为17 432 km²,占全流域沙漠化土地面积的13.5%;陕西和宁夏的沙漠化土地面积分别占全流域沙漠化土地面积的8.3%和6.5%。黄河流域的沙漠化空间格局主要是降水量与沙源空间耦合的结果,流域92.6%（119 114 km²）的沙漠化土地分布于干旱、半干旱地区。从20世纪70年代以来,黄河流域的沙漠化总体上经历了快速发展—发展放缓—明显逆转的过程,沙漠化大幅度的变化主要受人类活动影响所致,在过去几十年间风速持续减小对沙漠化逆转的积极作用也应引起重视。     

青藏高原大骨节病流行区环境、食物及人群硒水平研究

为了解青藏高原大骨节病流行区环境、食物及人群硒水平的分布现状,论文通过CNKI、Web of Science和Google Scholar等数据库检索收集了2000—2018年间发表的关于青藏高原大骨节病病区土壤、粮食、饮水及人发硒含量的文献,按拟定标准共筛选出33篇文献并提取有效数据进行统计分析。结果显示,青藏高原大骨节病病区自然环境总体仍处于低硒循环状态,其中耕作土壤、饮水平均硒含量分别为0.147 mg·kg^-1和0.54 μg·L^-1,高原自产青稞、小麦、糌粑平均含硒量分别为9.27、19.08、11.07 μg·kg^-1;而病区儿童整体硒营养水平较1990年以前有明显升高趋势,发硒平均含量为0.234 mg·kg^-1,已基本脱离硒缺乏状态(<0.20 mg·kg^-1);外源性大米的硒含量为43.29 μg·kg^-1;病区儿童发硒水平与土壤总硒含量无明显相关关系(r=0.125,P>0.05,N=23)。青藏高原大骨节病流行区儿童发硒水平与自然环境硒水平不一致的变化趋势提示,随着社会人文因素的干预,外源性硒的输入有所增加,并在一定程度上降低了病区人群对当地低硒环境的依赖性,这可能是青藏高原大骨节病病情稳定下降并得到有效控制的主要原因。