一种基于MODIS积雪产品的雪线高度提取方法

冰川雪线高度的遥感提取对冰川物质平衡研究具有重要意义。提出一种基于晴空环境下积雪覆盖频率的雪线高度提取方法。使用MOD10A1积雪产品中的像元积雪面积比例数据,提取了2000/2001-2014/2015年间高亚洲地区冰川消融期末雪线高度。使用实测的冰川年物质平衡资料和气象格网数据对提取的雪线高度变化的可信度进行分析。研究表明：近15 a高亚洲雪线高度变化及趋势具有明显的东西差异,雪线高度变化幅度自青藏高原内部地区向四周呈增加趋势,西部大于东部。提取的冰川雪线高度变化与观测的年物质平衡序列具有很高的相关性,对物质平衡波动的平均解释率可高达75%;与气象要素（气温、降水）的年际变化的相关性也较高,约61.58%的格网冰川雪线高度变化可以由夏季气温和季节降水解释。而高亚洲各分区冰川雪线高度的波动规律也与大气环流背景分布一致。因此提取的雪线高度变化具有冰川学意义,可以进一步应用于冰川物质平衡估算及模拟研究中。     

Analysis of extreme ground snow loads for Canada using snow depth records

Snow depth records from daily measurements at climatological stations were obtained from Environment Canada and were processed and analyzed. It was identified that there are 549 stations, each with at least 20 years of useable annual maximum snow depth data. Both the Gumbel distribution and generalized extreme value distribution were used to fit the annual maximum snow depth, considering several distribution fitting methods. Statistical analysis results indicated that, according to the Akaike information criterion, the Gumbel distribution is preferred for 72 % stations. The estimated return period value of annual maximum snow depth at stations was used to calculate their corresponding ground snow load. The at-site analysis results were used as the basis to spatially interpolate the ground snow loads for locations tabulated in the National Building Code of Canada (NBCC) since a code location and a climatological site are usually not co-located. For the interpolation, the ordinary co-kriging method with elevation as co-variate was used because a cross-validation analysis by using several deterministic and probabilistic spatial interpolation techniques indicated that the ordinary co-kriging method is preferred. A comparison of the newly estimated ground snow loads to those locations tabulated in the 1995 edition and 2010 edition of the NBCC was also presented.     

基于灾情数据的新疆雪灾的评估与区划

选取1955-2017年新疆因雪灾造成的牲畜受灾数量、作物受灾面积、经济损失三大灾情数据,揭示新疆雪灾时间和空间分布特征,通过构建雪灾灾损指数,对新疆雪灾进行区划。分析表明：新疆雪灾的发生北疆多于南疆,山区多于盆地,西部多于东部,主要集中发生在阿尔泰山区、天山山区及昆仑山区,灾情最重的是阿勒泰地区、伊犁地区、塔城地区;进入2000年以后,雪灾呈现明显高位震荡特点;利用多指标综合评估法,运用连续型随机变量分布密度近似求法,求得的新疆雪灾灾损指数,能客观反映新疆雪灾灾情;雪灾灾损指数遵从Γ分布,根据Γ分布函数属性及各类雪灾事件发生概率,进行了新疆雪灾等级区划。研究成果可为决策服务及气候可行性论证提供理论依据,对科学准确地制定防灾备灾措施、及时组织开展雪灾应急救助工作以及灾后的恢复重建决策给予一定参考。     

青藏高原积雪范围和雪深/雪水当量遥感反演研究进展及挑战

积雪是地表特征的重要参数,对辐射收支、气候和长期天气变化均有重要影响。雪本身又是一个重要的天气现象和水文气象参数,过量的降雪也会带来严重的雪灾,如牧区雪灾、雪崩和融雪洪水灾害等。因此对积雪的监测,尤其是对山区的积雪监测,具有多方面的意义。利用卫星遥感技术监测积雪已有50余年的历史,并已形成了系列业务产品。青藏高原平均海拔超过4 000 m,该地区的积雪具有重要的水文、气候和生态环境意义。由于地形复杂,人迹罕至,地面观测站点稀少,受较强太阳辐射的影响,积雪消融迅速、区域差异消融以及风吹雪等因素导致积雪分布破碎化严重,对使用遥感资料监测该地区的积雪造成的极大的困难和不确定性。随着国内外传感器技术的不断发展,光学和被动微波遥感数据的同步获取技术已经非常成熟,综合利用光学遥感数据高空间分辨率和被动微波数据不受云干扰的特点,结合机器学习、无人机等技术,将环境参数加入反演模型中,有助于提高青藏高原积雪参数反演精度。     

A technique for improving MODIS standard snow products for snow cover monitoring over Eastern Turkey

Daily snow cover maps at 500 m resolution are available from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensors on NASA's Terra and Aqua satellites. However, information of the ground can only be obtained during clear sky conditions. In addition, the spectral similarity between clouds and snow in the visible part of the spectrum causes some ambiguities in separation of them. The dynamic behavior of clouds with their tendency to move faster often enables discrimination between cloud and snow, since snow remains relatively stationary. In this study, a daily merging methodology to combine Terra and Aqua based snow maps are proposed, and the benefits of this combination are evaluated. The methodology was tested for Eastern Turkey covering the period from December 2002 to March 2003. Merged maps derived from the standard daily snow maps from Terra (MOD10A1) and Aqua (MYD10A1) reduce the cloud cover present in any one image alone, provide better representation of the surface snow cover, and indicate better agreement with ground snow measurements than when either one is used alone. For “on the ground data date” analyses yielded 31% match ratios for MOD10A1 and MYD10A1, whereas combined images enabled 38%. One- to 2-day temporal window shifting enabled further 7% and 5% improvements in match ratios, respectively. The maximum improvement of 25% was achieved in January of 2003. The proposed methodology provides an easy and effective way to improve snow cover maps with little effort and provides fewer cloud-contaminated images for snow monitoring, for hydrology, and water resource management.     

中国稳定积雪区IMS雪冰产品精度评价

IMS雪冰产品由多种光学与微波传感器数据融合而成,提供北半球每日无云的积雪范围,在积雪遥感研究中具有广阔的前景. 以气象站实测雪深数据为真值,检验了2009-2010年IMS雪冰产品在中国三大稳定积雪区北疆、东北、青藏高原地区每月、积雪季以及全年的误判率、漏判率和总体准确率,并分析了IMS雪冰产品的准确率与雪深之间的关系. 结果显示：IMS雪冰产品的年总体准确率在三大积雪区均超过了92%,积雪季总体准确率均超过了88%,利用IMS雪冰产品监测积雪范围是可靠的. 然而,IMS雪冰产品精度具有区域差异性,北疆地区在1月和2月误判率偏高,青藏高原地区积雪季有严重的漏判现象. IMS雪冰产品的准确率在东北地区和北疆地区随着雪深的增加而升高,当东北地区雪深超过6 cm,北疆地区超过13 cm时,准确率接近100%,但是,青藏高原地区两者基本没有关系. 通过在青藏高原地区与同时相的4景MODIS积雪产品对比分析发现,实际上IMS雪冰产品相对地高估了积雪面积,青藏高原地区漏判率高其原因是IMS对零碎积雪的识别能力不足并且气象站分布不均匀.     

Fuzzy factorial analysis of snow avalanches

The mutual influence of 21 factors pertaining to terrain, weather, forest and snowpack have been discussed by 10 experts. The semantic (nil, weak, medium, high) evaluations are translated as membership degrees of fuzzy sets, and averaged between the experts by taking their fuzzy expectation value, yielding a 21 × 21 fuzzy matrix for direct interactions. Fuzzy successive multiplications and additions of the matrix give the indirect interactions. The activity and the passivity of the factors with respect to avalanches is represented by directed weighted graphs, and the average semantic values by a 4 × 4 matrix: the activity (passivity) is weak (strong) for tree damage; medium (nil) for altitude; medium (weak) for ground shape; medium (medium) for vegetation, soil, wind, microclimate, tree type, vertical distribution and mechanics of trees, snowpack distribution and snow gliding; strong (nil) for ground inclination and exposition, weather; strong (medium) for heat, precipitation, horizontal distribution of trees, snowpack constitution and stability.     

The distribution of snow in Israel

The area researched is characterized by changes in the weather according to fluctuations in the world pressure belts throughout the year. The temperature climate in this region is typified in winter by changes in the pressure system from cyclonic to anti-cyclonic conditions in spells of about 5 days.When Israel is under the influence of a low pressure system, and the 500 mb surfaces are approximately 250 m lower than average, the temperatures along the vertical cross-sections are 6° – 8° C lower than average on rainy days, then snow may be expected in various sections of the country.This paper investigated the 35 years period from 1938/9 to 1971/2, concentrating on two mountain stations — Jerusalem in the centre of the country and Mt. Kena'an in the north —because of the high snow incidence and its greater regularity than other stations.The conclusions are as follows: The maximum number of days of snow (between 2/3 and 3/4 of the annual total) is noted at all stations during January and February, although there are also years in which there was snow but not during those months. Regional devision of the snowy days indicates that the number increases with altitude, and decreases from north to south.Snowfall in regions which have lower incidence is less stable and given to greater fluctuations. There was no indication of periodity of the number of days of snow per season, nor did it appear that a very snowy year followed a dry one.Investigation of the lenght of the snowy and the number of days of snowfall does not indicate any connection between them, and there may well be a long season with few days of snow, or a short one with relatively many snowy days.Comparison between the shorter period of readings — from 1938/9 till 1971/2 — and the longer one — from 1860/1 till 1971/2 (111 years) — at the Jerusalem station shows the same trends, which reinforce the conclusions presented above.     

2003-2010年青藏高原积雪及雪水当量的时空变化

利用MODIS逐日无云积雪产品与AMSR-E雪水当量产品进行融合, 获取了青藏高原500 m分辨率的高精度雪水当量产品, 通过研究青藏高原积雪时空动态变化特征, 分析了积雪覆盖日数、雪水当量以及总雪量的季节及年际变化. 结果表明: 青藏高原地区降雪主要集中在高海拔山区, 而高原腹地降雪较少, 降雪在空间上分布极为不均; 2003-2010年期间, 平均积雪日数呈显著减少趋势, 稳定积雪区面积在逐渐扩大, 常年积雪区面积在不断缩小. 与积雪日数时空变化相比, 雪水当量增加的区域与积雪日数增加的区域基本一致, 但喜马拉雅山脉在积雪日数减少的情况下雪水当量却在逐年增加, 表明该地区温度升高虽然导致部分常年积雪向季节性积雪过渡, 但降雪量却在增加. 总的积雪面积年际变化呈波动下降的趋势, 但趋势不显著, 且减少的比例很少. 最大积雪面积呈现波动上升后下降的趋势, 平均累积积雪总量呈明显的波动下降趋势, 年递减率为1.0×10³ m³·a^-1.     

NOAA IMS雪冰产品在青藏高原积雪监测中的适用性分析

在系统评估青藏高原积雪观测典型气象站历史定位坐标精度基础上,利用站点雪深资料对NOAA IMS 4 km和1 km分辨率雪冰产品在青藏高原的精度和适用性进行了验证和评估,定量分析了IMS 4 km到1 km空间分辨率提高和气象站历史定位与GPS定位坐标之间的差异对青藏高原IMS积雪监测精度的影响。结果表明：青藏高原个别气象站历史坐标与当前GPS接收机定位之间存在较大的差异,如安多气象站经度偏小0.6°,纬度偏大0.08°。IMS 4 km雪冰产品在青藏高原的总精度介于76.4%~83.2%,平均为80.1%,积雪分类精度介于35.8%~60.7%,平均为47.2%,平均误判率为17.1%,平均漏判率为45.5%,总体上呈现地面观测的积雪日数越多、平均雪深越大,其总体监测精度越低,而积雪分类精度越高的特点。IMS分辨率从4 km到1 km总体精度平均提高了2.9%,积雪分类精度平均提高了0.9%,主要是由于个别站点的精度提升较大引起的,对高原多数台站积雪监测精度的改进和提升很小。除个别台站外,目前气象站历史坐标和GPS定位坐标之间的差异,对IMS 4 km积雪监测精度验证结果没有影响。然而,今后随着卫星遥感技术的发展,更高时空分辨率的遥感积雪产品将用于积雪监测和研究,精确的地面观测站坐标信息是对这些遥感数据开展精度验证与实际应用的前提。     

基于MODIS积雪产品的天山年积雪日数空间分布特征研究

山区积雪是干旱区气候变化的重要指标因子,积雪日数与积雪分布之间有着密切关系。为了研究天山山区积雪日数空间分布特征,以MODIS8d积雪产品MOD10A2（Terra）和MYD10A2（Aqua）为数据源,首先对数据进行最大化合成,获取新疆天山500m×500m分辨率的年积雪日数,然后分析了2002-2014年13a积雪日的年际变化,并结合DEM数据分析了13a天山多年平均积雪日随高程和坡度的变化特征。结果表明：天山积雪日数分布极为不均,最大年平均积雪日数为193d,13a内天山绝大部分地区年积雪日变化趋势较为稳定,稳定区约占天山总面积的83.92%;在研究时段内天山总积雪日数主要集中在30d以内,其比例约为天山总面积的48%;各个高程带积雪日面积分布差异明显,但总体上积雪日数随着高程的增加而增加;从积雪日数随坡向分布来看,北坡、东北坡、东坡、西坡、西北坡所占面积比例（>30d）相对高于其他坡向。该研究结果对干旱区水资源估算具有参考意义。     

滑雪场积雪模拟研究进展

滑雪场的积雪条件是评价其盈利能力的关键指标。通过积雪模型准确地模拟滑雪场积雪条件对评估滑雪旅游目的地的气候风险具有重要意义,积雪模型亦可为滑雪场的造雪量和造雪时间提供科学依据。首先简要概述了滑雪场积雪模拟的发展趋势,介绍了具体的模拟方法。分析发现,应用于滑雪场积雪模拟的模型从早期的经验线性模型发展为简化的物理模型,再到精细的能量平衡模型;从仅模拟自然积雪发展为结合自然积雪和人工造雪,再到细致的考虑雪道维护过程。依据对自然积雪、人工造雪和雪道维护模拟方法的不同,滑雪场的积雪模拟方法可以简单的分为两类:单点或半分布式模拟和分布式模拟。单点或半分布式模拟将滑雪场作为一个点进行研究或等间隔地选取滑雪场多个海拔带进行模拟输出平均结果,通常采用简化的能量平衡模型和基于度日因子的概念模型,以简化的假设来解释人工造雪和雪道维护过程,以相对粗糙的方式表示雪的物理性质。分布式模拟在网格或雪道上划分出来的小单元尺度上对滑雪场积雪进行模拟,采用精细的多层能量平衡模型,详细描述了滑雪场的基础设施、人工造雪和雪道维护过程。由于我国大部分地区自然积雪资源匮乏,积雪模型在我国滑雪场地应用时需重点考虑人工造雪。在雪道人...     

祁连山区MODIS积雪反照率产品的精度验证及云下积雪反照率估算研究

积雪反照率在全球气候和能量收支平衡模型中起着重要的作用. 利用祁连山地区大冬树垭口站点反照率实测数据对由TM/ETM+得到的反照率数据进行标定, 然后将TM/ETM+反照率数据通过升尺度对MODIS逐日积雪反照率(SAD)产品在晴空条件下的精度进行了验证. 同时, 发展了一个基于MODIS SAD与AMSR-E SWE数据融合并结合Noah积雪反照率参数化方案估算MODIS SAD数据云下积雪反照率的算法, 通过统计分析纠正了云对积雪反照率的影响, 对云下积雪反照率进行了验证分析. 结果表明:MODIS SAD产品在祁连山地区的精度要低于大面积积雪覆盖的平坦地区(如格陵兰岛), 其平均绝对误差及均方根误差分别为0.0548和0.0727; 云下积雪反照率估算方法可以有效地获取云覆盖下积雪像元的反照率值, 纠正后的无云MODIS SAD数据与地面观测值有较好的一致性, 其平均绝对误差为0.078.     

A two-dimensional snow creep model for alpine terrain

The paper is focused on the investigation of snow creep and settlement of alpine snowpacks. A detailed investigation of the constitutive behavior of snow for long-term creeping under low-strain rates is provided. The snowpack is considered as a compressible, viscous fluid. An assumption for a temperature and density-dependent viscosity and viscous Poisson ratio is provided in this paper, based on observations in nature and laboratory. The provided model is particularly applicable for settlement calculation of seasonal alpine snowpacks and computation of reaction forces on obstacles. The model is evaluated by numerous laboratory and field experiments. Different case studies are provided in this paper.     

Dependence of thermal conductivity of snow on microstructure

A geometrical model, including different geometrical shapes influencing thermal conductivity of snow is proposed. The geometrical model has been assumed to comprise of unit cells having solid (ice) inclusion as an aggregation of spherical, cylindrical or cubical shapes with vertical connection, arranged in a cubic packing. From the geometrical model and one-dimensional heat transfer theory, the effective thermal conductivity has been computed. For this purpose, coupled one-dimensional heat transfer equations have been solved for steady-state condition to account for conduction in ice, conduction in air and latent heat transfer due to water vapour sublimation through air. The model demonstrates the dependency of thermal conductivity on density, grain-spacing, grain contact ratio and temperature. Spherical inclusions give highest conductivity while cubical inclusion estimates lowest value for the same density. Thermal conductivity has been found increasing sharply near to the packing density for all three shapes. Empirical model results and results obtained from existing microstructure based models have also been compared with the present model.     

Groundwater use for snow melting on the road

Japan is very rich in active volcanoes, so that the surveillance and prediction of volcanic activities are indespensable in order to protect human life and properties from catastrophic volcanic eruptions. The author intends here to review the volcanic activities in Japan and the history and status quo of the volcanological observation and research by the Japan Meteorological Agency (JMA), universities, etc. for that purpose. Needless to say, in Japan, volcanoes are studied from various view points such as geophysical, geochemical, geological and geographical. However, the observations for the purpose of detecting of reliable premonitory symptoms of volcanic eruptions are mainly based on geophysical and geochemical methods. In this country, there are permanent observatories at 18 active volcanoes, and there are also several mobile teams of volcanologists. It must be noticed that almost all the volcanic activities including very slight volcanic extraordinary phenomena are nowadays detected in Japan, and unexpected eruptions are very scarce, because the information on the actual state of activities of Japanese volcanoes are made public frequently.     

新疆积雪覆盖时空变异分析

利用2000-2010年MODIS积雪覆盖产品数据MOD10A2,提取了新疆近10年来积雪覆盖变化信息,并结合地面站点数据,对遥感积雪覆盖估算的精度进行了验证;分析了新疆积雪覆盖的年际、年内变化及南北疆积雪覆盖变化的差异;结合数字高程模型,分析不同高程带下积雪覆盖的时空变化规律,揭示高程因素对新疆积雪时空变化的影响。结果表明:MOD10A2提取的积雪信息能够反映新疆的积雪变化情况,总体精度达92.3%;近10年来,全疆年积雪覆盖率最大值范围为34.0%~51.7%,最小值范围为1.7%~2.6%;积雪覆盖比率的变化在南北疆差异明显,南疆区域积雪覆盖整体不高,年内积雪覆盖比率变化幅度低于50%;而北疆区域由于受复杂地形和气候带的影响,积雪覆盖比率大,年内的变化幅度强,除2008年均达到80%以上;在季节变化上,春季和秋季的积雪覆盖均值波动较为明显,夏季和冬季的积雪覆盖均值则波动较小,这一规律在北疆地区表现更为显著;积雪覆盖的时空分布与变化受高程的影响,在海拔4 000 m以下区域,夏季积雪覆盖比率低,冬季积雪覆盖比率高,而6 000 m以上海拔区域则表现出完全相反的特点,即夏季积雪覆盖比率高,冬季积雪覆盖比率低。     

基于MODIS数据的北疆积雪黑碳和雪粒径反演及时空变化分析

积雪是地球上反射率较高的自然表面,对于中高纬度地区的水文和能量收支平衡发挥着重要作用。表层积雪中的黑碳和雪粒径变化可以显著影响积雪反照率,造成积雪对太阳辐射吸收的变化,进而对区域气候变化和水文循环产生反馈作用。利用遥感技术对季节性积雪表层黑碳和雪粒径进行定量评估,可以获取时空上连续系统的雪表黑碳浓度和雪粒径变化情况,这也是许多气候和水文模型的输入因子。以中国主要季节性积雪区北疆为研究区,基于MODIS（Moderate Resolution Imaging Spectroradiometer）数据的3（0.47 μm）、2（0.86 μm）和5（1.24 μm）波段,采用SGSP（Snow Grain Size and Pollution Amount）算法反演2000-2018年积雪期的雪表黑碳浓度和雪粒径,并结合地面观测数据对于反演结果进行了精度验证,综合分析北疆雪表黑碳浓度和雪粒径时空变化趋势。结果显示,SGSP算法能够同时反演雪表黑碳浓度和雪粒径,并且验证结果表明纯雪像元上反演结果具有较好的精度;2000-2018年北疆雪表年均黑碳浓度和年均雪粒径都随时间变化呈现微弱下降趋势;受地理位置和局部污染源的影响,北疆积雪黑碳浓度空间分布复杂,天山北坡经济带平均黑碳浓度最高,伊犁地区平均黑碳浓度最低,雪粒径的空间分布显示塔城地区平均雪粒径最大,伊犁地区最小。     

新疆巴音布鲁克区域干雪及风吹雪介电常数对温度和低频频率的响应

积雪是干旱半干旱区珍贵的可再生水资源,因此许多研究机构和人员将西北地区积雪特性的研究视为重点,然而对积雪介电常数的研究仍很少。积雪介电常数是建立积雪特性反演模型和雪水资源评估的基础,也是雪灾防灾减灾工作的辅助信息。为了研究不同雪层干雪和风吹雪介电常数与温度和频率之间的关系,使用I-V方法测量了巴音布鲁克区域干雪表、中、底层雪和风吹雪在1Hz~1MHz低频范围,-30~0℃温度区间内介电常数的实部和虚部。结果表明在1Hz~1MHz频率范围内干雪及风吹雪的介电常数与温度和频率有关;其中介电常数实部与虚部对温度和频率的响应相似,但虚部对温度的响应更加敏感;同时介电常数对温度的响应要比对频率的响应明显。最后得出了干雪表、中、底层雪和风吹雪介电常数实部与温度之间的线性关系,R²均大于0.85。