西藏高原积雪覆盖空间分布及地形影响

利用2000-2014年MOD10A2积雪产品和数字高程模型DEM数据,以积雪覆盖率为指标,在分析西藏高原积雪空间分布特点的基础上,定量研究了高程、坡度和坡向等地形要素对高原积雪时空分布的影响。主要结论有：① 西藏高原积雪的空间分布差异显著,具有中东部念青唐古拉山和周边高山积雪丰富,覆盖率高,而南部河谷和羌塘高原中西部积雪少,覆盖率低的特点。② 海拔越高积雪覆盖率越高,积雪持续时间越长,年内变化越稳定。海拔2 km以下积雪覆盖率不足4%,海拔6 km以上覆盖率达75%。海拔4 km以下年内积雪覆盖呈单峰型分布特点,海拔越高,单峰型越明显;而海拔4 km以上则为双峰型,海拔越高,双峰型越明显。海拔6 km以下积雪覆盖率最低值出现在夏季,而6 km以上则出现在冬季。③ 总体上,高原地形坡度越高积雪覆盖率越高。不同坡向中,北坡积雪覆盖率最高,南坡最低,年内分布呈双峰型,而无坡向的平地积雪覆盖率要小于有坡向的山地,其年内变化呈单峰型分布特点。     

1979-2010年青藏高原积雪深度时空变化遥感分析

积雪深度是表征积雪特征的重要参数,也是区域气候变化最敏感的响应因子之一。利用1979-2010年逐日中国雪深长时间序列数据集,采用GIS空间分析和地统计方法,分析了青藏高原积雪深度的时空变化规律及异常空间分布特征。结果表明：近32年来,青藏高原雪深呈显著增加趋势,增加速率为0.26 cm/10a,其中,昆仑高寒荒漠地带雪深增加最为明显,增加速率达0.73 cm/10a;20世纪80年代至90年代青藏高原雪深呈逐步增加趋势,21世纪初变化平稳;青藏高原4个季节雪深变化均呈现为上升趋势,尤以冬季增加最为明显,增加速率达0.57 cm/10a。青藏高原东南、西部和南部为雪深分布高值区;逐像元回归分析表明,高原雪深呈增加趋势的像元数占全区像元总数的67.1%,其中有91.3%为轻度和中度增加,主要分布在高原北部和西部;最大雪深变化基本维持在-0.1~0.1 cm/a（45.47%）之间,在昆仑北翼山地、柴达木山地、羌塘高寒地带南部等局部地区最大雪深有增加趋势,主要是轻度增加,面积比例为36.66%。果洛那曲高寒地带、青南高寒地带和羌塘高寒地带为青藏高原积雪深度异常变化敏感区。     

昆仑山提孜那甫河流域雨雪分离的温度条件分析

气候变暖背景下高海拔山区融雪（冰）以及强降水引发的洪水愈加难以预测,通过山区雨雪分离可判定引发洪水的温度条件,从而为山洪准确预报提供简单而科学的参考依据。本研究以昆仑山提孜那甫河流域为例,基于流域内不同海拔气象站2012-2016年的降水以及温度数据,结合MOD10A2积雪数据,采用温度积分法和概率统计方法,利用研究期内的平均温度,确定出不同降水形态对应的温度条件,以达到雨雪分离的目的。研究结果表明,莫木克站最大温和积温分别达到20.91 ℃和51.82 ℃时,降水可判定为降雨,最大温和积温分别低于18.13 ℃,43.69 ℃时,降水可判定为降雪;库地站最大温和积温分别达到14.51 ℃,33.17 ℃时,降水可判定为降雨,最大温和积温分别低于13.57 ℃,31.68 ℃时,降水可判定为降雪;西合休站最大温和积温分别达到9.43 ℃,19.53 ℃时,降水可判定为降雨,最大温和积温分别低于8.22 ℃,19.4 ℃时,降水可判定为降雪。利用流域内气象站点附近乡镇的气象统计数据对温度条件及分离结果进行验证,在海拔2000 m以下、2000~3000 m以及3000 m以上不同海拔地区的准确率分别为92.86%、79.49%以及88.3%。本研究可为判别洪水类型和洪水预报提供科学参考。     

基于动态NDSI阈值的每日积雪监测方法

准确掌握积雪覆盖信息对于气象、水文和全球气候变化研究都具有重要的意义。遥感技术在进行大范围、高频率的积雪覆盖监测中发挥着重要的作用。目前,SNOMAP算法是用于积雪遥感监测最普遍的技术手段,其核心是利用固定阈值的归一化差分积雪指数（Normalized Difference Snow Index,NDSI）进行积雪识别,但这种方法忽略了积雪光谱信息的时相变化,会产生积雪监测的误差。本文提出了一种动态NDSI阈值方法,以纯永久积雪像元的平均NDSI值作为参照系调整固定的NDSI阈值,从而削减影像光谱值波动对积雪识别的影响。以三江源地区作为研究区域,将基于每日MODIS数据进行积雪监测最佳的NDSI阈值与同日纯永久积雪像元的平均NDSI值作线性回归,通过每日纯永久积雪像元平均NDSI值的变化来调整用于积雪识别的NDSI阈值。结果表明：① 基于每日MODIS数据进行积雪覆盖监测最佳的NDSI阈值与同日纯永久积雪像元的平均NDSI值之间存在较好的线性关系,决定系数R ²达到0.86;② 三江源地区动态NDSI阈值的范围为0.29~0.37,其平均值在0.33左右,说明MODIS全球积雪面积产品中将NDSI阈值取为0.40会低估三江源地区的积雪面积;③ 与采用固定NDSI阈值0.33的监测方法相比,动态NDSI阈值法近似率、总体分类精度和F值的平均值分别提高了5.17%、0.70%、1.14%。     

一种考虑雪粒径变化的积雪面积反演算法

光谱混合分析能够提取亚像元信息,被广泛地应用于遥感影像目标探测之中。本文针对MODIS积雪遥感影像,基于光谱混合分析框架,利用渐进辐射传输模型建立不同粒径大小的雪反射率光谱库,提出了一种考虑端元变化及二次辐射的雪盖面积反演算法。此算法首先利用渐进辐射传输模型建立不同粒径大小积雪的反射率光谱库,然后使用序贯最大角凸锥方法获取植被、土壤与岩石、阴影的光谱库。在建立各种地物反射率光谱库之后,利用均方根误差最小的方法获取最优端元组合。在此基础上,考虑端元独立辐射以及积雪与其它地物的二次辐射过程,利用稀疏光谱混合模型获取积雪面积与雪粒径大小。实验结果表明：此方法能够同时反演雪粒径与积雪面积,反演的雪粒径相比单波段的渐进辐射传输模型小,反演的积雪面积相比MOD10A1产品精度略微提高。     

基于多源数据的青藏高原雪深重建

青藏高原地形复杂,积雪时空分布异质性较强且大部分地区积雪较薄,而被动微波遥感因其空间分辨率低以及雪深反演中的不确定性,极大地限制了其反演青藏高原雪深的精度。本文尝试将多源遥感数据以及与积雪模型（SnowModel）相结合,来重建更高质量的青藏高原雪深数据。首先,利用MODIS积雪面积比例产品,根据构建的积雪衰减曲线以及经验的融合规则对低分辨率被动微波雪深进行了降尺度;然后,结合MODIS/被动微波融合雪深数据和SnowModel对研究区进行雪深数据同化实验;最后,利用地面站实测雪深数据对MODIS/被动微波融合雪深以及同化输出雪深的精度进行了分析和对比。结果表明,基于数据同化方法得到的雪深数据更接近地面观测雪深值,通过均方根误差以及相关系数的对比,同化雪深结果优于MODIS/被动微波融合雪深结果。     

2002-2012年青藏高原积雪物候变化及其对气候的响应

积雪是地表最活跃的自然要素之一,其动态变化对气候、环境以及人类生活都产生了重要影响。本文利用MODIS积雪产品和IMS雪冰产品,首先通过Terra、Aqua双星合成和临近日合成去除MODIS积雪产品中的部分云像元,再与IMS融合,获取了青藏高原2002-2012年逐日无云积雪覆盖产品,并逐像元计算每个水文年的积雪覆盖日数（SCD）、积雪开始期（SCS）和积雪结束期（SCE）,分析了不同生态分区积雪的时空变化特征,以及积雪开始期和结束期与温度、降水的关系。结果表明：青藏高原积雪分布存在明显的空间差异,南部喜马拉雅山脉和念青唐古拉山地区以及西部帕米尔高原和喀喇昆仑山脉为SCD的2个高值区,年均积雪日数在200 d以上。18.1%的区域SCS表现出明显的提前趋势,主要集中在青藏高原中东部;羌塘高原南部、念青唐古拉山西段以及川西地区有显著推迟趋势,占高原面积的8.5%。23.2%的区域SCE显著推迟,主要集中在果洛那曲高寒区、昆仑山区和念青唐古拉山地区;而仅有6.9%的区域表现出提前趋势,主要分布在高原西南部。总体上,不同生态单元内积雪开始与结束期受温度、降水的影响差异很大,表现出不同的空间格局与演变趋势。     

The Glacial（MIS 3-2） Outlet Glacier of the Marsyandi Nadi- icestream-network with its Ngadi Khola Tributary Glacier（Manaslu- and Lamjung Himalaya）：The Reconstructed Lowering of the Marsyandi Nadi Ice Stream Tongue down in to the Southern Himalaya Foreland

For the reconstruction of past climate variations,investigations on the history of glaciers are necessary.In the Himalaya,investigations like these have a rather short tradition in comparison with other mountains on earth.At the same time,this area on the southern margin of Tibet is of special interest because of the question as to the monsoon-influence that is connected with the climate-development.Anyhow,the climate of High Asia is of global importance.Here for the further and regionally intensifying answer to this question,a glacial glacier reconstruction is submitted from the CentralHimalaya,more exactly from the Manaslu-massif.Going on down-valley from the glacial-historical investigations of 1977 in the upper Marsyandi Khola(Nadi) and the partly already published results of field campaigns in the middle Marsyandi Khola and the Damodar- and Manaslu Himal in the years 1995,2000,2004 and 2007,new geomorphological and geological field- and laboratory data are introduced here from the Ngadi(Nadi) Khola and the lower Marsyandi Nadi from the inflow of the Ngadi(Nadi) Khola down to the southern mountain foreland.There has existed a connected ice-stream-network drained down to the south by a 2,100-2,200 m thick and 120 km long Marsyandi Nadi main valley glacier.At a height of the valley bottom of c.1,000 m a.s.l.the Ngadi Khola glacier joined the still c.1,300 m thick Marsyandi parent glacier from the Himalchuli-massif(Nadi(Ngadi) Chuli) – the south spur of the Manaslu Himal.From here the united glacier tongue flowed down about a further 44 km to the south up to c.400 m a.s.l.(27°57'38 "N/84°24'56" E) into the Himalaya fore-chains and thus reached one of or the lowest past ice margin position of the Himalayas.The glacial(LGP(Last glacial period),LGM(Last glacial maximum) Würm,Stage 0,MIS 3-2) climatic snowline(ELA = equilibrium line altitude) has run at 3,900 to 4,000 m a.s.l.and thus c.1,500 altitude meters below the current ELA(Stage XII) at 5,400-5,500 m a.s.l.The reconstructed,maximum lowering of the climatic snowline(ΔELA = depression of the equilibrium line altitude) about 1,500 m corresponds at a gradient of 0.6°C per 100 altitude meters to a High Glacial decrease in temperature of 9°C(0.6 × 15 = 9).At that time the Tibetan inland ice has caused a stable cold high,so that no summer monsoon can have existed there.Accordingly,during the LGP the precipitation was reduced,so that the cooling must have come to more than only 9°C.     

Mass balance sensitivity to climate change: A case study of glacier No. 1 at urumqi riverhead, Tianshan Mountains, China

Ｓｔｕｄｉｅｓｓｈｏｗｔｈａｔｔｈｅｇｌｏｂａｌｓｕｒｆａｃｅｔｅｍｐｅｒａｔｕｒｅｗｉｌｒｉｓｅａｔａｒａｔｅｏｆ０．３℃ｐｅｒ１０ａｉｎｔｈｅｎｅｘｔｓｅｖｅｒａｌｄｅｃａｄｅｓ（Ｈｏｕｇｈｔｏｎｅｔａｌ．,１９９０,１９９２;Ｓｈｉ,１９９５;...     

Snow cover variation and streamflow simulation in a snow-fed river basin of the Northwest Himalaya

Snowmelt is an important component of any snow-fed river system.The Jhelum River is one such transnational mountain river flowing through India and Pakistan.The basin is minimally glacierized and its discharge is largely governed by seasonal snow cover and snowmelt.Therefore,accurate estimation of seasonal snow cover dynamics and snowmeltinduced runoff is important for sustainable water resource management in the region.The present study looks into spatio-temporal variations of snow cover for past decade and stream flow simulation in the Jhelum River basin.Snow cover extent(SCE) was estimated using MODIS(Moderate Resolution Imaging Spectrometer) sensor imageries.Normalized Difference Snow Index(NDSI) algorithm was used to generate multi-temporal time series snow cover maps.The results indicate large variation in snow cover distribution pattern and decreasing trend in different sub-basins of the Jhelum River.The relationship between SCE-temperature,SCE-discharge and discharge-precipitation was analyzed for different seasons and shows strong correlation.For streamflow simulation of the entire Jhelum basin Snow melt Runoff Model(SRM) used.A good correlation was observed between simulated stream flow and in-situ discharge.The monthly discharge contribution from different sub-basins to the total discharge of the Jhelum River was estimated using a modified version of runoff model based on temperature-index approach developed for small watersheds.Stream power - an indicator of the erosive capability of streams was also calculated for different sub-basins.     

山东半岛降雪时空分布特征

以2000-2018年MODIS MOD10A1日产品数据为数据源,结合数字高程模型（DEM）及降水量、风向等气象数据,构建了积雪空间分布模型,能够有效地提取强降雪区域。以此为基础,利用相关分析、缓冲区分析等方法,探究山东半岛降雪时空分布特征,结果表明：① 将NDSI累积量与DEM数据相结合,能够有效构建山东半岛积雪空间分布模型,实现了对山东半岛强、弱降雪区域提取,NDSI累积量≥150的区域中,在强降雪区的面积占降雪范围的79.78%;② 降雪区域存在空间差异,呈现北多南少,东多西少的分布格局,以黄、渤海与山东半岛海陆分界线为基准,离岸距离39.1 km范围内降雪多,离岸距离39.1 km以外降雪少;山脉150 m高程线北侧迎风坡降雪多,南侧背风坡降雪少;③ 山东半岛强降雪年以3-5年为周期存在年际变化。探究山东半岛降雪长时间序列时空分布特征,在收集淡水资源,缓解用水紧张和灾害预防方面具有一定意义。     

青藏高原MODIS图像冰雪信息挖掘与动态监测分析

针对西藏地区的环境特点,对适应大区域尺度及全球尺度的MODIS数据进行冰雪信息挖掘与监测的方法进行改进,提出适于西藏地区MODIS数据雪信息提取的技术途径,实践表明,该法对冰雪情监测具有较好的实用性。     

NORMALIZED DIFFERENCE SNOW INDEX SIMULATION FOR SNOW-COVER MAPPING IN FOREST BY GEOSAIL MODEL

1INTRODUCTION Snow is an important component of the Earth's surface. Up to 50×106km2(34%) ofthe Earth's land surface is sea- sonally snow-covered (VIKHAMAR and SOLBERG, 2002).Comparedtootherlandcovers,snowcoverextent varies dramatically on very short time scales (hours- months). Its presence affects physical, chemical and bio- logical processes at many spatial scales and has impor- tant social impacts. At the global scale, its high albedo strongly influences the Earth's radiation …     

THE CLIMATIC CONDITIONS OF POLAR-TYPE GLACIERS DEVELOPMENT IN CHINA

ＴＨＥＣＬＩＭＡＴＩＣＣＯＮＤＩＴＩＯＮＳＯＦＰＯＬＡＲ－ＴＹＰＥＧＬＡＣＩＥＲＳＤＥＶＥＬＯＰＭＥＮＴＩＮＣＨＩＮＡ￥ＨｕａｎｇＭａｏｈｕａｎ（黄茂桓）（ＬａｎｚｈｏｕＩｎｓｔｉｔｕｔｅｏｆＧｌａｃｉｏｌｏｇｙａｎｄＧｅｏｃｒｙｏｌｏｇｙ，ｔｈｅＣ...     

Monitoring snow-cover area change in Antarctic coastline region using MODIS product data

Based on MODIS snow products, this article studied the changes of snow cover area during 2003-2006 along the coastline of the Antarctic, and 18 typical regions were chosen for further analysis. The result showed that the change of snow cover area was in a fluctuant downward trend as a whole, and more fluctuated obviously in warm season than in cold season. In temporal scale: for the season cycle, the snow cover extent increased rapidly in cold season (Apr-Oct), while the performance in warm season (Nov-Mar) was not exactly the same during the four years, the snow cover extent decreased in the first and then increased in 2004 and 2006, however, increased firstly and then decreased but reduced as a whole in 2005, for the inter-annual cycle, snow cover extent was the largest in 2003, but reached to the lowest level in 2004, and then increased gradually in 2005 and 2006, whereas, it declined with fluctuant as a whole. In spatial scale, changes mainly centralized along the coastline, moreover, it was more remarkable in the West Antarctic than in the East Antarctic, especially in the Antarctic Peninsula region.     

Wind-forced equatorial wave dynamics of the Pacific Ocean during 2014/2015 and 2015/2016 E1 Ni?o events

The equatorial wave dynamics of interannual sea level variations between 2014/2015 and2015/2016 El Nino events are compared using the Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics,Institute of Atmospheric Physics Climate Ocean Model(LICOM) forced by the National Centers for Environmental Prediction(NCEP) reanalysis I wind stre s s and heat flux during 2000-2015.In addition,the LICOM can reproduce the interannual variability of sea surface temperature anomalies(SSTA) and sea level anomalies(SLA) along the equator over the Pacific Ocean in comparison with the Hadley center and altimetric data well.We extracted the equatorial wave coefficients of LICOM simulation to get the contribution to SLA by multiplying the meridional wave structure.During 2014/2015 El Nino event,upwelling equatorial Kelvin waves from the western boundary in April2014 reach the eastern Pacific Ocean,which weakened SLA in the eastern Pacific Ocean.However,no upwelling equatorial Kelvin waves from the western boundary of the Pacific Ocean could reach the eastern boundary during the 2015/2016 El Nino event.Linear wave model results also demonstrate that upwelling equatorial Kelvin waves in both 2014/2015 and 2015/2016 from the western boundary can reach the eastern boundary.However,the contribution from stronger westerly anomalies forced downwelling equatorial Kelvin waves overwhelmed that from the upwelling equatorial Kelvin waves from the western boundary in 2015.Therefore,the western boundary reflection and weak westerly wind burst inhibited the growth of the 2014/2015 El Nino event.The disclosed equatorial wave dynamics are important to the simulation and prediction of ENSO events in future studies.     

SOME IMPORTANT CHARACTERISTICS OFSURFACE ALBEDO OF CHINA IN RECENT YEARS

Sdsealbedoisdefinedasaratioofreflectivesolarradiationtototalacradiation,whichcanindicatetheaborptionandreflectivecapabilityofthe~hforsolarradiation.Itisan~tfactortoinfluencethes~radiationbalance,espeiallyfortheformationandvariationoflocalndcroclirnate.Generally,theactualdatafromo~tionalstationSforrebationareedintheanalysisofsolacealal.HOwever,therepresentativenessOfthedateisIratedduetothescarcityofobservationalstationandunevendistributionofs~condition.Afterthe1960s,theanalysisfors~ealbedoise…     

Interannual to decadal variation of spring sea level anomaly in the western South China Sea

Satellite observations of sea level anomalies(SLA) from January 1993 to December 2012 are used to investigate the interannual to decadal changes of the boreal spring high SLA in the western South China Sea(SCS) using the Empirical Orthogonal Function(EOF) method. We find that the SLA variability has two dominant modes. The Sea Level Changing Mode(SLCM) occurs mainly during La Ni?a years, with high SLA extension from west of Luzon to the eastern coast of Vietnam along the central basin of the SCS, and is likely induced by the increment of the ocean heat content. The Anticyclonic Eddy Mode(AEM) occurs mainly during El Ni?o years and appears to be triggered by the negative wind curl anomalies within the central SCS. In addition, the spring high SLA in the western SCS experienced a quasi-decadal change during 1993–2012; in other words, the AEM predominated during 1993–1998 and 2002–2005, while the La Ni?a-related SLCM prevailed during 1999–2001 and 2006–2012. Moreover, we suggest that the accelerated sea level rise in the SCS during 2005–2012 makes the SLCM the leading mode over the past two decades.     

Significance of glacio-morphological factors in glacier retreat: a case study of part of Chenab basin,Himalaya

A study has been carried out in part of Chenab basin,Himalaya to understand the relationship between glacio-morphological factors and change in glacial area. Initially change in areal extent of glaciers was derived for two time frames(1962-2001/02 and 2001/02-2010/11). The study comprised of 324 glaciers for the monitoring period of 1962-2001/02 for,which 11% loss in glacial area was observed. Two hundred and thirty-eight glaciers were further monitored between 2001/02 and 2010/11. These glaciers showed an area loss of 1.1%. The annual deglaciation has been found to be higher during the period of 1962-2001/02 compared to 2001/02-2010/11. The spatial and temporal variability in deglaciation was also addressed usingglacio-morphic parameters. Area,length,percentage of debris cover,and various elevation parameters of glaciers were observed to have significant controls on relationships to the rate of glacial shrinkage. Largerarea and longer glaciers show a lower percentage of retreat than smaller and shorter ones. Moreover,glaciers located at lower altitudes and having gentle slopes show more area retreat. The results of area retreat in debris covered and debris free glaciers supports that the glaciers covered by debris retard ice melting at some extent. 158 glaciers were observed having no debris cover,and these exhibit 14% of loss in surface area. In glaciers having 40% debris cover,8% of deglaciation was observed. The glaciers located below equilibrium line altitude(ELA) have experienced 4.6% of deglaciation for the time frame 2001/02 – 2010/11 whereas it was found to be 1.1% for the glaciers occurring above ELA. However,theorientation of glaciers did not show any considerable influence on glacial change based on hypothesis.     

An Improved Method for Modeling Spatial Distribution of δD in Surface Snow over Antarctic Ice Sheet

Using the recent compilation of the isotopic composition data of surface snow of Antarctic ice sheet, we proposed an improved interpolation method of δD, which utilizes geographical factors (i.e., latitude and altitude) as the primary predictors and incorporates inverse distance weighting (IDW) technique. The method was applied to a high-resolution digital elevation model (DEM) to produce a grid map of multi-year mean δD values with 1km spatial resolution for Antarctic& The mean absolute deviation between observed and estimated data in the map is about 5.4‰, and the standard deviation is 9‰. The resulting δD pattern resembles well known characteristics such as the depletion of the heavy isotopes with increasing latitude and distance from coast line, but also reveals the complex topographic effects.