西藏自治区位于青藏高原西南部(26°50′-36°53′N，78°25′-99°06′E)，土地面积122.84万km²，其中青稞种植面积达13.958万hm^2[9].全区气温低、昼夜温差大、积温少，年均气温为-2.4~12.1 ℃，自东南向西北递减；年降水量在66.3~894.5 mm之间，降水量少且季节性差异明显，同样呈自东南向西北递减分布，因此气候呈现西北严寒干燥、东南温暖湿润的特点.

本文收集了西藏自治区5个青稞主产市内主要28个国家基本气象站的基本气象数据，站点分布如图 1，详细概况见表 1.所有气象数据来自中国气象科学数据共享网(http://cdc.cma.gov.cn/).

本研究采用Microsoft Excel 2013处理所有实验数据，运用滑动平均法、Mann-Kendall检验法和Pearson相关分析等方法分析了西藏青稞主产市在青稞春、冬生长季中日照时数的年际时空变化特征.

分析拉萨市4个气象站1956-2018年生长季年均日照时数代表该市青稞生长季的日照水平(图 2)，可看出拉萨市60多年春、冬青稞生长季日照时数变化较为平缓，都或高或低地在均值附近摆动.但春、冬生长季日照时数均在2009年后呈明显的下降趋势，下降速率每10年分别为65.87 h，41.49 h.由5年滑动平均曲线可看出，拉萨的春、冬生长季日照时数走向相似，均在1960-1973年实现1次“减-增”相抵，1973-1990年有微弱起伏，从1990-2018年则完成了2次周期为14年的“增-减”运动.

利用Mann-Kendall检验对拉萨青稞在春、冬生长季日照时数进行趋势分析和突变检验.春生长季的日照时数变化如图 2c中UF_k曲线所示，可看出整体呈先降后升趋势. 20世纪50年代中期到60年代初下降趋势显著，70年代到21世纪初则上升显著，UF_k和UB_k曲线相交于1971-1972年间，但并未和3条UB_k曲线皆相交，因此该交点并非是真正的上升突变点.进一步看拉萨冬生长季日照时数的MK检验，同样呈先降后升趋势且没有突变点，但下降并不显著，而上升趋势在2007年后的显著性较春生长季强.

用山南市内6个站点(浪卡子、贡嘎、泽当、错那、隆子、加查气象站)1956-2018年春、冬生长季年均日照时数代表该市青稞生长季的日照水平(图 3)，发现春生长季中日照时数每10年以3.68 h波动下降.具体而言，其在1956-1963年呈线性下降，而后在1963-1972年，1972-1991年，1991-2004年，2004-2018年这4个时间段形成了4次有层次递减的“升-降”运动.冬生长季日照时数走向与春生长季相似，但起伏变小.

从MK检验中可以更明显地看出，春、冬生长季日照时数分别在1967年、2002年UF曲线超过0.05显著性临界线，表明下降趋势显著.进一步分析下降突变点发现，春生长季中UF_k和UB_k，0.5 UB_k，1.5 UB_k曲线在春生长季中并未相交，而在冬生长季中相交于1988-1989年之间，表明山南市冬生长季日照时数在此之间有一下降突变点.对比春、冬日照时数(图 3a，3b)走向也可发现春生长季中在1983-1992年间的下降趋势中存在1个1988-1990年微弱上升的转折，而在冬生长季中从1983-1992年仅为单一的下降趋势，下降速率每10年为13.05 h.

用昌都市内7个站点(丁青、昌都、类乌齐、洛隆、八宿、左贡、芒康气象站)1956-2018年青稞春、冬生长季年均日照时数代表该市在青稞生长季的日照水平(图 4)，发现不同于其他主产市的波动下降，昌都市春、冬生长季日照时数呈明显的上升趋势，对应的气候倾向率每10年分别为4.81 h，3.022 h.根据MK检验可进一步得到春、冬生长季日照时数的UF统计量分别在1975年、1971年超过显著性临界线且没有突变点，表明上升显著且平稳.

用林芝市内4个站点(林芝、米林、波密、察隅气象站)1961-2018年春、冬生长季中年均日照时数代表该市青稞生长季的日照水平(图 5)，可看出相比于其他4个主产区而言，林芝市的日照水平相对较低，春、冬皆在20世纪60年代中期后呈现单一显著的下降趋势，且没有突变点的出现，表明下降趋势明显且稳定.春、冬青稞生长季中日照时数最高值分别出现在1973年、1966年，为160.6 h，162.4 h，最低值分别出现在2004年、1995年，为104.7 h，128.0 h.冬生长季日照时数下降速度比春生长季慢，波动也更剧烈，同年冬生长季年均日照时数比春生长季高，因此对于西藏地区地势较低的林芝市而言，冬青稞种植比春青稞更有优势.

用日喀则市内7个站点(南木林、拉孜、日喀则、聂拉木、定日、江孜、亚东气象站)1956-2018年生长季年均日照时数代表该市青稞生长季的日照水平(图 6)，根据春青稞生长季中日照时数的气候倾向率，发现60多年来日喀则市日照时数呈较为剧烈的波动下降.结合5年滑动平均线与MK检验结果看，1956-1962年下降明显，1960-2013年从239.5 h到240.1 h，这50多年来虽有几次较大起伏，并在2000-2001年之间有一显著的下降突变点，但都能较快恢复到240 h左右，而从2013-2018年则呈线性下降.

近几十年来，在城市化、城市扩张和全球变暖的大背景下，对于日照时数变化趋势的研究也大量开展.相关研究表明全球大多数地区日照时长呈下降趋势^[10].全国范围内1954-2005年日照总体呈下降趋势^[11]，且具有明显的季节和局地差异.季节差异上，我国夏季和冬季降幅大于春季和秋季^[12].局地差异上，华南地区日照时数在1954-2012年下降明显^[13]；云贵高原北部下降，而西南部呈上升趋势^[14]；西北地区1960-2009年大部分地区表现为显著的减少趋势，平均每10年减少13.6 h^[15]；东北地区年日照时数也呈显著减少趋势，平均每10年减少40.5 h，除秋季日照时数变化不显著外，其他季节日照时数均显著减少^[16].本文的结果也表明西藏5个主产区中除昌都外皆呈减少趋势，并且春生长季的降幅最明显，但春、冬生长季下降幅度远小于青藏高原西南部^[17]和整个西南地区^[18]的速率.

西藏青稞5个主产市的日照时数在春、冬生长季的年际变化也表现出了不同的空间差异(图 7).整体上，冬生长季日照时数的气候倾向率比春生长季高(日喀则春生长季平均气候倾向率每10年为-3.76 h，其他4市春生长季平均气候倾向率每10年为-3.65 h，冬生长季为-1.95 h)，仅拉萨、山南、昌都的部分站点春或冬生长季日照时数为上升趋势，其他都为下降趋势：

1) 西藏西部的日喀则市在春生长季日照时数皆呈下降趋势，尤其是聂拉木、拉孜、南木林这3站的春生长季日照时数下降剧烈，对应气候倾向率每10年分别为-6.32 h，-10.16 h，-6.39 h，说明过去每10年中春生长季日照日数会分别下降31.6 h，50.78 h，31.94 h.

2) 中部拉萨、山南、林芝3市中尼木、加查两站春、冬生长季日照时数上升明显，浪卡子、隆子、拉萨、错那、察隅等仅冬生长季日照时数略微增加.相反地，墨竹工卡、贡嘎春、冬生长季日照时数均下降显著，浪卡子、米林等在春生长季中下降剧烈，其他站点的变化相较之下并不明显.

3) 昌都市的6站中，除洛隆、八宿的冬生长季气候倾向率为负外，其余皆为正，且昌都北部的丁青、昌都两站春、冬生长季日照时数上升迅速.对于春生长季气候倾向率，除洛隆、八宿、芒康春生长季日照时数下降明显外，其他变化不大.

从5市中分别选取拉萨、日喀则、泽当、昌都、林芝气象站点进行日照时数与其他7个主要气象要素的相关性分析可得到：

春生长季日照时数与日照百分率呈明显的正相关，日照百分率是农业生产中所需气象要素的重要因子，百分率的计算是用天文日照百分率法计算的，即该地区实际的日照时数除以无云层遮挡，理想条件下的实日照时数，也就是天文可照时数.日照百分率相对日照时数更能反映该地区天气状况对日照时数的影响.数据分析结果表明，除昌都站点外，其他4站春生长季日照时数与降水量、平均水汽压、平均相对湿度均呈显著负相关(表 2).其中春生长季日照时数与降水量的负相关最好的是日喀则站点(相关系数达到-0.768)，与平均水汽压、平均相对湿度负相关性最强的为拉萨站点(相关系数分别为-0.705，-0.624).日喀则、泽当、林芝站点春生长季日照时数与本站平均气压呈负相关，而在拉萨、昌都并不显著.平均风速与日照时数呈正相关的有日喀则、泽当，平均气温与日照时数呈正相关的为拉萨、日喀则、昌都、林芝，但较其他气象因子而言，平均气温与日照时数的相关性并不强.

整体上冬生长季日照时数与各气象因子的相关性较春生长季弱(表 2).冬生长季日照时数与降水量、平均相对湿度有显著负相关的有拉萨、泽当站点，降水量负相关性最强的从春生长季的日喀则站点的-0.768降到拉萨站点的-0.456，相关系数降低了40.6%，平均相对湿度负相关性最强的拉萨站点从春生长季的-0.624降到-0.512，相关系数降低了17.9%.除昌都站外其他3站与平均水汽压均呈显著负相关，而与本站平均气压有负相关的为昌都和林芝站.与平均风速有正相关的仅为泽当站，且相关性减弱.与平均气温的正相关性在冬生长季不再突出，只有昌都站的平均气温与冬生长季日照时数有一定的正相关性.

青稞是西藏自治区最大的优势和特色作物，是西藏自治区基本口粮作物和主要粮食作物，播种面积占全区农作物种植面积的70%以上^[19]，并且其属于喜光作物，对日照时长较为敏感，研究表明充足的日照才能满足青稞生长期光合作用的需要，充足的日照和温度的协调配合能促进青稞的稳产，是实现青稞高产的优势条件之一^[20].西藏自治区的高原气候有利于青稞的生产，但同时青稞对气候变化背景下青藏高原气候演变也较为敏感^[21].因此，研究西藏自治区年际日照时长的演变不仅有助于理解青藏高原气候演变过程及其对全球变化的响应，也对青藏高原青稞的适应机制具有重要的意义.

西藏海拔较高，空气稀薄，太阳能资源丰富，光照时间长，在作物生长季月平均日照时数达250 h，显著高于同纬度地区.光照充分有利于提高作物的光合作用，将光能资源转化为生物能资源，同时也有利于改善作物的品质，增加总糖、还原糖和淀粉含量.

青稞籽粒的形成主要依靠开花后的光合产物，此时光照、温度、水分等气象因子对光合特性影响很大，研究表明，在长日照地区如西藏的日喀则，能显著提高青稞的光合能力，增加有机物质的合成，使青稞的“源”得以增加.而在西藏的林芝日照时数相对弱的地区，青稞的灌浆强度较弱，植株对有机物的同化减缓，表现为籽粒千粒质量较低，秕粒现象增加.青稞的灌浆期主要集中在7月份，这段时间西藏降水较多，云层的增多及增厚显著减弱了日照时数，生产中日照时数的减少，可以通过调整青稞的播种时间，改变青稞生育期进行弥补.研究表明春青稞提早播种可以使生育期延长，而延迟播种则使生育期推后.在拉萨的春青稞4月10日播种，生育期为123 d，而4月30日播种生育期则是113 d，所以提前播种20 d可以使生育期延长10 d^[22]，因为早播时环境温度较低，青稞生长缓慢，生育期则变长，从而可以弥补因日照时数减少、光合时间不足对青稞造成的影响.提前播种，可以在雨季正式来临之前使青稞完成灌浆，这样可以在青稞的日照敏感期避开低日照时间，减缓日照时数较低对青稞的影响.

西藏除昌都外年际间日照时数呈现下降趋势，西部日喀则地区近10年来作物生长季日照时数减少30 h以上，藏东南部日照时数减少的趋势更为明显，减少量达104.7 h以上，这主要是降水增加因云层遮挡导致日照时数减少，目前日照时数的减少量尚不足以对青稞产量及品质产生较大的影响，相反，导致日照时数减少的驱动因素即降水增加在西藏干旱地区更能提高作物产量，但是日照时数减少的这种趋势需要注意.在气候变化的大环境下，青藏高原是气候变化的敏感区，作物的生长容易受到气候变化的影响，需要研究相应的对策措施，以应对气候变化.

从日照时数的下降季节来看，冬生长季日照时数下降速度比春生长季慢，同年冬生长季年均日照时数比春生长季高.这主要是冬季水汽活动较弱或很少，这样在作物的布局上，适合冬青稞种植的河谷区域，应尽量种植冬青稞，充分利用冬生长季日照时数下降不大的有利条件.在实际生产中冬青稞在产量上要比春青稞有优势，平均每666.7 m²产量高100 kg以上.

相关性分析表明春生长季日照时数与降水量、平均水汽压、平均相对湿度均呈显著负相关，从全国尺度看，日照时数与气温日较差之间相关性较弱^[23]，也有研究表明除青藏高原外，我国气温日较差的变化与日照时数的变化密切相关^[24-25]，这与本文得到的西藏自治区大部分地区气温对日照时数相关性较弱相似.局部研究上，东北地区日照时数与风速呈正相关关系，而与气温、降水量、相对湿度等呈负相关关系^[16]；西南地区日照时数影响因素同样发现风速与日照时数呈显著正相关，而与湿度呈负相关^[17].青海省与祁连山地区日照时数则与平均水汽压、总云量、低云量呈显著负相关^[26-27].本文中通过Pearson相关分析也得到春生长季日照时数与降水量、平均相对湿度、平均水汽压呈负相关，而与风速呈正相关.这说明不同地区日照时数与主要气象因子的相关性质是相同的，但相关性的强弱会有所差别.另外，人类对自然资源不合理开发利用造成的空气污染对日照时数也会产生一定影响，近年来青藏高原人类活动明显加剧，公路修建、旅游业开发以及荒地开垦所产生的扬尘对日照时数均有影响，但是这些人类活动对日照时数的影响范围、影响程度均需要进一步的研究.

1) 各市年际变化趋势上，除昌都市春、冬生长季日照时数呈上升趋势外，其他市均呈明显下降趋势，其中林芝市在生长季中年均日照时数较其他各市低，且下降也更为剧烈.局部时间段上，5市春、冬生长季中年均日照时数在2013-2018年皆呈线性下降趋势.

2) 突变检验中，仅日喀则的春生长季、山南的冬生长季日照时数走向中存在下降突变点，说明在突变点前后这一时间段存在一个剧烈的单一下降趋势.拉萨、昌都、林芝其波动较小、变化平缓且有规律.

3) 各站点气候倾向率分布上，冬生长季日照时数的气候倾向率比春生长季高，仅拉萨、山南、昌都的部分站点春或冬生长季日照时数为上升趋势，其他多为下降趋势.其中，昌都、尼木、加查、丁青这4个站点的春、冬气候倾向率均为正值，说明以上站点在生长季中日照时数随时间具有一定的线性增长趋势，也会进一步促进青稞生长.

4) 从日照时数与气象要素的相关性上，春生长季日照时数与降水量、平均水汽压、平均相对湿度呈显著负相关，而与日照百分率呈明显的正相关，而冬生长季日照时数与各气象因子的相关性质不变，但关系减弱.