2000_2010年吉林省积雪时空变化特征及其与气候的关系_杨倩
合集下载
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
MODIS 积 雪 产 品 利 用 云 和 雪 的 波 谱 特 征 ,即 在 可见光 波 段 (0.545~0.565μm)有 相 似 的 波 谱 特 性,而在短波 近 红 外 波 段 (1.628~1.652μm)差 异 明显,建立 NDSI指数可区分雪和大部分积云,但 是 难以区 分 薄云 和雪。除 考 虑 森 林 覆 盖区 的 情 况, MOD10A1积雪识 别 算 法 还 引 入 水 陆 阈 值 设 定、云 掩膜和热掩膜区分 薄 云 和 雪,由 此 最 终 确 定 积 雪 面 积 。 [10-11] MOD10A1 数 据 的 精 度 受 云 覆 盖、积 雪 深 度、植被覆盖类型和 土 地 利 用 类 型 等 多 种 因 素 综 合 影响。MODIS产品精度较 高,总 体 精 度 达 到 93%, 在部分庄 稼 地 和 农 田 精 度 达 到 99%。 MOD10A01 产品可较好地消除云对地表积雪分类精度的影响, 晴天时 MOD10A1总精度可达到 98.5%,积雪分类 精度为 98.2% 。 [12-14] 故 在 此 选 择 MOD10A1 产 品 提取雪盖信息。
基 于 不 同 观 测 资 料,我 国 很 多 学 者 对 中 国 不 同 地区积雪的分布及长期变化趋势与气候的响应关系
等进行了多方面研究。青藏高原东部是高原积雪年 际 变 化 最 显 著 的 地 区,柯 长 青 等[2]综 合 SMMR、 NOAA 和站台实测数据发现,从 20 世纪60 年 代到 80 年代积雪年际波动幅度有明显增加趋势,且 存在 3a左右的准周期。 长 江 黄 河 源 区 积 雪 变 化 存 在 一 高 一 低 两 个 中 心 ,通 过 16 个 气 象 站 台 逐 日 积 雪 资 料 发 现 长 江 源 和 黄 河 源 头 地 区 降 雪 较 早 ,积 雪 季 节 长 , 结束晚,消退过程 缓 慢;而 在 黄 河 源 东 部 地 区,降 雪 较晚,积 雪 发 展 比 较 缓 慢,消 退 过 程 迅 速 。 [3] 利 用 2000~2009年 MODIS监测的8日积雪合成产品和 气象 数 据,卢 新 玉 等 发 [4] 现 在 新 疆 玛 纳 斯 河 流 域 随 着海拔的降低,气 候 对 积 雪 变 化 的 影 响 越 明 显。 积
积 雪 日 数 是 在 整 个 积 雪 季 时 间 段 内 像 素i 持 续 有 积 雪 覆 盖 的 时 间 ,定 义 为 :
第3期 杨 倩等:2000~2010年吉林省积雪时空变化特征及其与气候的关系 4 15
N
∑ SCD = (Di2 -Di1) i=0
(1)
其中:Di1 和 Di2 分别为每个积雪季的起始 日 期 和 终
and Its Relationships to Climate Change from 2000to 2010[J].Remote Sensing Technology and Application,2012, 27(3):413-419.[杨倩,陈圣波,路鹏,等.2000~2010年吉林省积雪时 空 变 化 特 征 及 其 与 气 候 的 关 系 [J].遥 感 技 术 与 应 用 ,2012,27(3):413-419.]
止日期,分别对应10月1日和次年的 4 月30 日;N
为整个图像总的像元个数。
4.2 积雪初日 SCOD
为计 算 SCOD,需 要 计 算 10 月 1 日 至 12 月 1
2000~2010年吉林省积雪时空变化特征 及其与气候的关系
杨 倩,陈圣波,路 鹏,杨 莹,王明常
(吉林大学地球探测科学与技术学院,吉林 长春 130026)
摘要:利用 Terra卫星提供的2000年10月1日到2010年 4 月 30日每日雪覆盖产品 MOD10A1, 提取研究区积雪覆盖指数 SCI、积雪日数 SCD、积雪初日 SCOD 及积雪终日 SCMD 遥感信息,结合 同 期 吉 林 省 界 内 23 个 地 面 气 象 观 测 站 的 同 期 气 温 和 降 水 资 料 ,分 析 该 区 积 雪 的 变 化 特 征 与 气 温 和 降水的关系。结果表明:① 吉林省大部分地区积雪日数为30~90d,东部山区积雪持续时间长、积 雪初日日期早以及积雪终日日期晚,中西部地 区 变 化 情 况 相 反;② 积 雪 覆 盖 指 数 SCI呈 波 浪 式 变 化,与积雪季气温呈负相关;③ 积雪日数与气温呈反相关、与降水量呈正相关,与积雪季气温、夏季 降水量的 相 关 系 数 分 别 为 -0.7407、0.6875;积 雪 初 日 情 况 相 反,与 积 雪 季 气 温、夏 季 平 均 气 温 为 0.743、0.5479;积 雪 终 日 与 气 温 呈 反 相 关、与 降 水 量 呈 正 相 关,与 积 雪 季 气 温、夏 季 降 水 量 为 -0.5214、0.4647。 积 雪 指 数 均 对 气 温 的 变 化 更 敏 感 ,气 温 升 高 导 致 积 雪 初 日 推 迟 、积 雪 终 日 提 前 , 从 而 使 积 雪 日 数 减 小 ;积 雪 季 降 水 量 的 增 加 有 利 于 积 雪 日 数 增 大 ,而 积 雪 终 日 的 推 迟 有 利 于 夏 季 降 水量的增加。 关 键 词 :吉 林 省 ;积 雪 ;气 温 ;降 水 ;MODIS 中 图 分 类 号 :P 426.63;TP 79 文 献 标 志 码 :A 文 章 编 号 :1004-0323(2012)03-0413-07
2 研究区概况
吉 林 省 位 于 东 北 地 区 的 中 部 ,地 处 北 温 带 ,属 温 带大陆性季 风 气 候,地 理 位 置 121°~131°E、41°~ 46°N 之 间。 研 究 区 域 东 西 长 650 km,南 北 宽 300km,总面积18.74万 km2,占全国 总 面 积 的 2% (图1)。吉 林 省 地 貌 形 态 差 异 明 显,地 势 由 东 南 向 西北倾斜,呈明显 的 东 南 高、西 北 低 的 特 征,形 成 了 山地、丘 陵、平 原 三 大 地 貌 类 型,由 此 可 分 为 东、中、 西 3 个 区 域 。 其 中 ,东 部 地 区 森 林 覆 盖 率 较 高 ,西 部 地区主要以荒漠和 盐 碱 地 为 主,而 中 部 城 镇 分 布 较 为密集 。 [9]
图 1 研 究 区 地 理 图
Fig.1 Geographic map of study area
3 数据及预处理
3.1 气 象 数 据 将 每 年 10 月 1 日 到 次 年 4 月 30 日 定 义 为 一 个
积 雪 季 ,次 年 5 月 1 日 到 9 月 30 日 定 义 为 夏 季 。 气 象数 据 包 括 吉 林 省 界 内 23 个 气 象 站 台 从 2000 年 10月1日到2010年9月30 日 的 逐 日 气 温、降 水 量 数 据 ,和 长 春 、延 吉 、安 图 及 吉 林 4 个 站 台 从 2008 年 10月1日 到 2009 年 4 月 30 日 的 逐 日 积 雪 深 度 数 据。前者用于分析 积 雪 变 化 与 气 候 的 响 应 关 系,后
收 稿 日 期 :2011-06-27;修 订 日 期 :2011-12-30 基金项目:教育部新世纪优秀人才计 划 项 目 (NCET-07-0353),中 央 高 校 基 本 科 研 业 务 费 专 项 资 金 (200903046),国 家 自 然 科 学 基 金 项 目
(40974047)资 助 。 作者简介:杨 倩(1987-),女,甘肃定西人,硕士研究生,主要从事积雪定量遥感研究。E-mail:yangqian10@ mails.jlu.edu.cn。 通 讯 作 者 :陈 圣 波 (1967- ),男 ,河 南 信 阳 人 ,教 授 、博 导 ,主 要 从 事 定 量 遥 感 研 究 。E-mail:chensb@jlu.edu.cn。
41 4 遥 感 技 术 与 应 用 第 27 卷
雪的变化主要是气温和降水量等多种气候要素共同 作 用 的 结 果 。 一 般 来 说 ,积 雪 与 气 温 呈 负 相 关 ,与 降 水呈正相 关 。 [5-8] 本 文 基 于 MODIS 数 据 产 品 提 取 积 雪 时 空 变 化 信 息 ,建 立 积 雪 覆 盖 指 数 ,计 算 积 雪 初 日 、积 雪 终 日 及 积 雪 日 数 ,进 而 评 价 吉 林 省 积 雪 指 数 变 化 特 征 与 气 温 、降 水 的 相 关 关 系 。
第 27 卷 第 3 期 2012 年 6 月
REMOTE S遥EN SI感NG T技EC HN术OL OG与Y AN应D AP用PLICATION
Vol.27 No.3 Jun.2012
引 用 格 式:Yang Qian,Chen Shengbo,Lu Peng,et al.Spatio-temporal Variation of Snow Cover in Jilin Province
4 积雪指数提取方法
本文将 研 究 区 内元(编码为200)赋值为1,非雪像 元 赋 值 为 0。 在 每 个积雪季逐日分类 图 像 叠 加 的 基 础 上,计 算 出 吉 林 省2000~2010 年 每 年 积 雪 持 续 时 间,即 积 雪 日 数 SCD(Snow Cover Duration/days)。 由 此 ,进 一 步 计 算积雪覆盖指数 SCI(Snow Cover Index)、积 雪 初 日 SCOD(Snow Cover Onset Dates)以 及 积 雪 终 日 SCMD(Snow Cover Melting Dates)[15-17]。 4.1 积雪日数 SCD
1 引 言
积 雪 是 地 球 表 面 最 为 活 跃 的 自 然 要 素 之 一,积 雪 面 积 、雪 深 、雪 的 反 射 率 等 是 全 球 能 量 平 衡 模 型 中 的重要输入参数 。 [1] 积雪通过高反照率和低热 传 导 率影响局部地区乃 至 全 球 的 季 节 性 气 候 变 化;积 雪 又是干旱半干旱区 水 资 源 的 重 要 来 源,尤 其 在 春 季 雪融阶段。有效的积雪监测有助于研究气候系统变 化对自然生 态 系 统 带 来 的 重 大 影 响,能 为 气 象、气 候 、农 业 以 及 防 灾 减 灾 等 众 多 领 域 提 供 信 息 服 务 。
者用于验证积雪指数的精度。
3.2 卫 星 遥 感 数 据 本文使 用 MODIS 雪 覆 盖 每 日 合 成 数 据 产 品
MOD10A1,来 源 于 美 国 国 家 雪 冰 数 据 中 心 (http:// nsidc.org//),空 间 分 辨 率 为 500 m。 研 究 区 在 全 球 正弦曲线投影 SIN(Sinusoidal projection)系统中每 天有两个数据,编号为 h26v04和 h27v04,数据 等 级 质量 为 3 级,数 据 格 式 为 EOS-HDF。 选 取 2000~ 2010年每年10月1日到次 年 4 月 30 日 的 数 据,10 个积 雪 季 节 4 864 个 数 据。 利 用 MRT(Map Re- project Tools)进行 数 据 格 式 转 换、投 影 变 换、重 采 样和图 像 拼 接,把 HDF 格 式 转 换 为 Geotiff格 式, 将 SIN 地 图 投 影 转 换 为 WGS84/UTM 系 统,重 采 样方法为最临近 法。 最 后,利 用 吉 林 省 矢 量 边 界 将 影 像 裁 剪 ,获 取 研 究 区 每 日 积 雪 遥 感 影 像 。
基 于 不 同 观 测 资 料,我 国 很 多 学 者 对 中 国 不 同 地区积雪的分布及长期变化趋势与气候的响应关系
等进行了多方面研究。青藏高原东部是高原积雪年 际 变 化 最 显 著 的 地 区,柯 长 青 等[2]综 合 SMMR、 NOAA 和站台实测数据发现,从 20 世纪60 年 代到 80 年代积雪年际波动幅度有明显增加趋势,且 存在 3a左右的准周期。 长 江 黄 河 源 区 积 雪 变 化 存 在 一 高 一 低 两 个 中 心 ,通 过 16 个 气 象 站 台 逐 日 积 雪 资 料 发 现 长 江 源 和 黄 河 源 头 地 区 降 雪 较 早 ,积 雪 季 节 长 , 结束晚,消退过程 缓 慢;而 在 黄 河 源 东 部 地 区,降 雪 较晚,积 雪 发 展 比 较 缓 慢,消 退 过 程 迅 速 。 [3] 利 用 2000~2009年 MODIS监测的8日积雪合成产品和 气象 数 据,卢 新 玉 等 发 [4] 现 在 新 疆 玛 纳 斯 河 流 域 随 着海拔的降低,气 候 对 积 雪 变 化 的 影 响 越 明 显。 积
积 雪 日 数 是 在 整 个 积 雪 季 时 间 段 内 像 素i 持 续 有 积 雪 覆 盖 的 时 间 ,定 义 为 :
第3期 杨 倩等:2000~2010年吉林省积雪时空变化特征及其与气候的关系 4 15
N
∑ SCD = (Di2 -Di1) i=0
(1)
其中:Di1 和 Di2 分别为每个积雪季的起始 日 期 和 终
and Its Relationships to Climate Change from 2000to 2010[J].Remote Sensing Technology and Application,2012, 27(3):413-419.[杨倩,陈圣波,路鹏,等.2000~2010年吉林省积雪时 空 变 化 特 征 及 其 与 气 候 的 关 系 [J].遥 感 技 术 与 应 用 ,2012,27(3):413-419.]
止日期,分别对应10月1日和次年的 4 月30 日;N
为整个图像总的像元个数。
4.2 积雪初日 SCOD
为计 算 SCOD,需 要 计 算 10 月 1 日 至 12 月 1
2000~2010年吉林省积雪时空变化特征 及其与气候的关系
杨 倩,陈圣波,路 鹏,杨 莹,王明常
(吉林大学地球探测科学与技术学院,吉林 长春 130026)
摘要:利用 Terra卫星提供的2000年10月1日到2010年 4 月 30日每日雪覆盖产品 MOD10A1, 提取研究区积雪覆盖指数 SCI、积雪日数 SCD、积雪初日 SCOD 及积雪终日 SCMD 遥感信息,结合 同 期 吉 林 省 界 内 23 个 地 面 气 象 观 测 站 的 同 期 气 温 和 降 水 资 料 ,分 析 该 区 积 雪 的 变 化 特 征 与 气 温 和 降水的关系。结果表明:① 吉林省大部分地区积雪日数为30~90d,东部山区积雪持续时间长、积 雪初日日期早以及积雪终日日期晚,中西部地 区 变 化 情 况 相 反;② 积 雪 覆 盖 指 数 SCI呈 波 浪 式 变 化,与积雪季气温呈负相关;③ 积雪日数与气温呈反相关、与降水量呈正相关,与积雪季气温、夏季 降水量的 相 关 系 数 分 别 为 -0.7407、0.6875;积 雪 初 日 情 况 相 反,与 积 雪 季 气 温、夏 季 平 均 气 温 为 0.743、0.5479;积 雪 终 日 与 气 温 呈 反 相 关、与 降 水 量 呈 正 相 关,与 积 雪 季 气 温、夏 季 降 水 量 为 -0.5214、0.4647。 积 雪 指 数 均 对 气 温 的 变 化 更 敏 感 ,气 温 升 高 导 致 积 雪 初 日 推 迟 、积 雪 终 日 提 前 , 从 而 使 积 雪 日 数 减 小 ;积 雪 季 降 水 量 的 增 加 有 利 于 积 雪 日 数 增 大 ,而 积 雪 终 日 的 推 迟 有 利 于 夏 季 降 水量的增加。 关 键 词 :吉 林 省 ;积 雪 ;气 温 ;降 水 ;MODIS 中 图 分 类 号 :P 426.63;TP 79 文 献 标 志 码 :A 文 章 编 号 :1004-0323(2012)03-0413-07
2 研究区概况
吉 林 省 位 于 东 北 地 区 的 中 部 ,地 处 北 温 带 ,属 温 带大陆性季 风 气 候,地 理 位 置 121°~131°E、41°~ 46°N 之 间。 研 究 区 域 东 西 长 650 km,南 北 宽 300km,总面积18.74万 km2,占全国 总 面 积 的 2% (图1)。吉 林 省 地 貌 形 态 差 异 明 显,地 势 由 东 南 向 西北倾斜,呈明显 的 东 南 高、西 北 低 的 特 征,形 成 了 山地、丘 陵、平 原 三 大 地 貌 类 型,由 此 可 分 为 东、中、 西 3 个 区 域 。 其 中 ,东 部 地 区 森 林 覆 盖 率 较 高 ,西 部 地区主要以荒漠和 盐 碱 地 为 主,而 中 部 城 镇 分 布 较 为密集 。 [9]
图 1 研 究 区 地 理 图
Fig.1 Geographic map of study area
3 数据及预处理
3.1 气 象 数 据 将 每 年 10 月 1 日 到 次 年 4 月 30 日 定 义 为 一 个
积 雪 季 ,次 年 5 月 1 日 到 9 月 30 日 定 义 为 夏 季 。 气 象数 据 包 括 吉 林 省 界 内 23 个 气 象 站 台 从 2000 年 10月1日到2010年9月30 日 的 逐 日 气 温、降 水 量 数 据 ,和 长 春 、延 吉 、安 图 及 吉 林 4 个 站 台 从 2008 年 10月1日 到 2009 年 4 月 30 日 的 逐 日 积 雪 深 度 数 据。前者用于分析 积 雪 变 化 与 气 候 的 响 应 关 系,后
收 稿 日 期 :2011-06-27;修 订 日 期 :2011-12-30 基金项目:教育部新世纪优秀人才计 划 项 目 (NCET-07-0353),中 央 高 校 基 本 科 研 业 务 费 专 项 资 金 (200903046),国 家 自 然 科 学 基 金 项 目
(40974047)资 助 。 作者简介:杨 倩(1987-),女,甘肃定西人,硕士研究生,主要从事积雪定量遥感研究。E-mail:yangqian10@ mails.jlu.edu.cn。 通 讯 作 者 :陈 圣 波 (1967- ),男 ,河 南 信 阳 人 ,教 授 、博 导 ,主 要 从 事 定 量 遥 感 研 究 。E-mail:chensb@jlu.edu.cn。
41 4 遥 感 技 术 与 应 用 第 27 卷
雪的变化主要是气温和降水量等多种气候要素共同 作 用 的 结 果 。 一 般 来 说 ,积 雪 与 气 温 呈 负 相 关 ,与 降 水呈正相 关 。 [5-8] 本 文 基 于 MODIS 数 据 产 品 提 取 积 雪 时 空 变 化 信 息 ,建 立 积 雪 覆 盖 指 数 ,计 算 积 雪 初 日 、积 雪 终 日 及 积 雪 日 数 ,进 而 评 价 吉 林 省 积 雪 指 数 变 化 特 征 与 气 温 、降 水 的 相 关 关 系 。
第 27 卷 第 3 期 2012 年 6 月
REMOTE S遥EN SI感NG T技EC HN术OL OG与Y AN应D AP用PLICATION
Vol.27 No.3 Jun.2012
引 用 格 式:Yang Qian,Chen Shengbo,Lu Peng,et al.Spatio-temporal Variation of Snow Cover in Jilin Province
4 积雪指数提取方法
本文将 研 究 区 内元(编码为200)赋值为1,非雪像 元 赋 值 为 0。 在 每 个积雪季逐日分类 图 像 叠 加 的 基 础 上,计 算 出 吉 林 省2000~2010 年 每 年 积 雪 持 续 时 间,即 积 雪 日 数 SCD(Snow Cover Duration/days)。 由 此 ,进 一 步 计 算积雪覆盖指数 SCI(Snow Cover Index)、积 雪 初 日 SCOD(Snow Cover Onset Dates)以 及 积 雪 终 日 SCMD(Snow Cover Melting Dates)[15-17]。 4.1 积雪日数 SCD
1 引 言
积 雪 是 地 球 表 面 最 为 活 跃 的 自 然 要 素 之 一,积 雪 面 积 、雪 深 、雪 的 反 射 率 等 是 全 球 能 量 平 衡 模 型 中 的重要输入参数 。 [1] 积雪通过高反照率和低热 传 导 率影响局部地区乃 至 全 球 的 季 节 性 气 候 变 化;积 雪 又是干旱半干旱区 水 资 源 的 重 要 来 源,尤 其 在 春 季 雪融阶段。有效的积雪监测有助于研究气候系统变 化对自然生 态 系 统 带 来 的 重 大 影 响,能 为 气 象、气 候 、农 业 以 及 防 灾 减 灾 等 众 多 领 域 提 供 信 息 服 务 。
者用于验证积雪指数的精度。
3.2 卫 星 遥 感 数 据 本文使 用 MODIS 雪 覆 盖 每 日 合 成 数 据 产 品
MOD10A1,来 源 于 美 国 国 家 雪 冰 数 据 中 心 (http:// nsidc.org//),空 间 分 辨 率 为 500 m。 研 究 区 在 全 球 正弦曲线投影 SIN(Sinusoidal projection)系统中每 天有两个数据,编号为 h26v04和 h27v04,数据 等 级 质量 为 3 级,数 据 格 式 为 EOS-HDF。 选 取 2000~ 2010年每年10月1日到次 年 4 月 30 日 的 数 据,10 个积 雪 季 节 4 864 个 数 据。 利 用 MRT(Map Re- project Tools)进行 数 据 格 式 转 换、投 影 变 换、重 采 样和图 像 拼 接,把 HDF 格 式 转 换 为 Geotiff格 式, 将 SIN 地 图 投 影 转 换 为 WGS84/UTM 系 统,重 采 样方法为最临近 法。 最 后,利 用 吉 林 省 矢 量 边 界 将 影 像 裁 剪 ,获 取 研 究 区 每 日 积 雪 遥 感 影 像 。