压实与融化情景下表层积雪性质反演的可行性研究.docxVIP

压实与融化情景下表层积雪性质反演的可行性研究 摘要：探究积雪光谱反演不同情景下表层积雪性质，可为利用遥感反演技术和加快大面积监测积雪的物理特性提供更丰富的理论和实践参考。积雪的雪水当量与密度是积雪出流、产流的关键指标，对其研究对东北地区春季融雪产流及融雪型土壤侵蚀具有现实意义。选择东北地区自然降雪，通过人工模拟压实、融化试验，得到不同压实程度与不同融化时间的雪样。观测表层积雪的光谱及雪水当量、积雪密度等特征，分析表层积雪光谱发射率与两种积雪性质的关系，探讨利用表层积雪光谱特征反演雪水当量及雪密度的可行性。结果表明：压实和融化两种处理导致雪密度与雪水当量发生改变，二者变化趋势基本一致，融化条件下积雪性质变化幅度更大。利用光谱反射率的一阶微分反演压实条件下积雪的雪水当量与融化条件下积雪的雪水当量精度较高，回归模型的决定系数分别为0.850和0.809；光谱反射率的一阶微分反演压实条件下积雪密度的回归模型决定系数为0.784，精度较高。但融化组积雪密度未能利用光谱特征得到有效反演。两种试验处理下，特征波段存在差异，主要原因是融化较压实具有较高的含水量，导致雪的粒径及孔隙度降低。 积雪是地球表面最为活跃的自然要素之一，影响着区域、全球的能量和水分平衡，对气候变化、水分平衡及人类活动等方面都有重要的影响1 材料与方法1.1 雪样处理试验采用东北地区自然降雪，雪样原始状态为自然降雪沉积一周后的状态，雪样采集和观测地点为沈阳农业大学校园内（41°49′30″N,123°33′25″E），观测时间为2019年1月1～4日。针对本次试验的情景，设置了两种模拟试验，压实试验是用烧杯盛装室外自然积雪，人工进行不同程度的压实，通过机械作用改变雪水当量、雪密度，共取雪样19份；融化试验是每隔20min取室外自然积雪盛满于500m L烧杯内，使其在恒温实验室内自然融化，盛装雪样时尽量避免人为干扰，保证积雪原始状态，其目的是通过改变雪样融化时间以此改变雪水当量、雪密度，共取雪样19份。1.2 雪样观测使用ASD Field Spec Hand Held 2地物光谱仪测定雪样表层光谱反射率。光谱范围为325～1075nm，波长精度±1nm。选择晴朗无风的天气在室外宽阔平地观测，两组试验雪样观测时间为中午12∶00，将光谱仪垂直向下对准水平放置的参照白板进行光谱仪校正后采集光谱数据。取表层2cm深度雪样于500m L烧杯内，测定表层雪样的雪水当量和雪密度。其中雪水当量指当积雪完全融化后，所得到的水形成水层的垂直深度；雪密度指单位体积雪的质量。两种指标计算方法为：式中：SWE为雪水当量（mm);h每个雪样重复观测10次后取平均值。测量过程中，始终保证光谱仪镜头位于烧杯中供试雪样上方15cm处。为消除野外测量时噪音等其他因素的影响，光谱曲线需进行均值化、去噪和平滑预处理。每次测量都由同一人按照仪器操作规范完成，以减少系统误差。1.3 数据处理压实试验记为压实组，融化试验记为融化组。两种供试雪样分别包括19组数据，按雪水当量、密度从小到大排列后取每3个相邻数据作为率定值，第4个数据作为验证值，可得14组数据用于率定，5组数据用于验证，分别探讨光谱反演雪水当量、密度的可行性。处理光谱特征与积雪性质关系时需对积雪光谱反射率进行光谱变换，包括倒数（1/R）、对数（log R）、一阶微分（R′）、倒数的一阶微分（1/R′）、对数的一阶微分[log(1/R)]，对率定组雪水当量、雪密度及其上述光谱反射率变换结果进行偏相关关系分析，按照R式中：y2 结果与分析2.1 不同试验处理下积雪的性质由图1可知，压实组积雪的密度范围在0.25～0.57g·cm2.2 不同试验处理下积雪的光谱特征由图2可知，压实组积雪光谱反射率范围为0.28～0.98，光谱反射率随波长的变化大体均呈先增大后减少最后再增大的趋势，积雪光谱曲线整体表现为325～355nm波段快速上升，355～755nm波段上升速度变缓，说明研究区积雪在可见光波段变化较为平稳，与钱昊由图3可知，融化组积雪光谱反射率范围为0.42～0.99，积雪光谱曲线整体表现为325～705nm平稳上升，705～845nm波段基本不变，845nm波段以后光谱反射率随波段的增大发生小范围波动，865～925nm光谱反射率随波长增大小幅度上下波动，925～1075nm光谱反射率波动较大。有3条曲线在325～575nm波段光谱反射率上升缓慢，斜率较小，它们的积雪密度、雪水当量分别为0.26g·cm两种试验处理的光谱反射率整体相似，随波长的变化大体均呈先增大后减少最后再增大的趋势，可见光范围内缓慢上升，进入近红外区后开始急剧下降，在1025nm附近反射率降至极小值，出现反射光谱