资源描述：
第 ３ ２ 卷 第 ９ 期 ２ ０ １ ９ 年 ９ 月 环 境 科 学 研 究 Ｒ ｅ ｓ ｅ ａ ｒ ｃ ｈ ｏ ｆ Ｅ ｎ ｖ ｉ ｒ ｏ ｎ ｍ ｅ ｎ ｔ ａ ｌ Ｓ ｃ ｉ ｅ ｎ ｃ ｅ ｓ Ｖ ｏ ｌ . ３ ２ ꎬ Ｎ ｏ . ９ Ｓ ｅ ｐ . ꎬ ２ ０ １ ９ 收 稿 日 期 ２ ０ １ ８ ￣ ０ ６ ￣ ０ ６ 修 订 日 期 ２ ０ １ ９ ￣ ０ ５ ￣ ０ ６ 作 者 简 介 孙 成 １ ９ ９ ２ ￣ ꎬ 男 ꎬ 河 北 保 定 人 ꎬ ４ １ ５ ０ ５ ４ ４ ３ ６ ＠ ｑ ｑ . ｃ ｏ ｍ . ∗ 责 任 作 者 ꎬ 王 卫 １ ９ ５ ８ ￣ ꎬ 男 ꎬ 河 北 承 德 人 ꎬ 教 授 ꎬ 硕 士 ꎬ 博 导 ꎬ 主 要 从 事 土 地 利 用 ∕ 覆 被 变 化 生 态 环 境 效 应 研 究 ꎬ ｗ ａ ｎ ｇ ｗ ｅ ｉ ＠ ｈ ｅ ｂ ｔ ｕ . ｅ ｄ ｕ . ｃ ｎ 基 金 项 目 国 家 自 然 科 学 基 金 项 目 Ｎ ｏ . ４ １ ４ ７ １ ０ ９ １ ꎻ 河 北 省 自 然 地 理 学 重 点 学 科 项 目 Ｎ ｏ . ２ ０ ０ ５ ￣ ７ ５ Ｓ ｕ ｐ ｐ ｏ ｒ ｔ ｅ ｄ ｂ ｙ Ｎ ａ ｔ ｉ ｏ ｎ ａ ｌ Ｎ ａ ｔ ｕ ｒ ａ ｌ Ｓ ｃ ｉ ｅ ｎ ｃ ｅ Ｆ ｏ ｕ ｎ ｄ ａ ｔ ｉ ｏ ｎ ｏ ｆ Ｃ ｈ ｉ ｎ ａ Ｎ ｏ . ４ １ ４ ７ １ ０ ９ １ ꎻ Ｋ ｅ ｙ Ｄ ｉ ｓ ｃ ｉ ｐ ｌ ｉ ｎ ｅ ｓ ｏ ｆ Ｎ ａ ｔ ｕ ｒ ａ ｌ Ｇ ｅ ｏ ｇ ｒ ａ ｐ ｈ ｙ ｉ ｎ Ｈ ｅ ｂ ｅ ｉ Ｐ ｒ ｏ ｖ ｉ ｎ ｃ ｅ ꎬ Ｃ ｈ ｉ ｎ ａ Ｎ ｏ . ２ ０ ０ ５ ￣ ７ ５ 基 于 线 性 混 合 效 应 模 型 的 河 北 省 Ｐ Ｍ ２ ５ 浓 度 时 空 变 化 模 型 研 究 孙 成 １ ꎬ ２ ꎬ 王 卫 １ ꎬ ２ ∗ ꎬ 刘 方 田 １ ꎬ ２ ꎬ 郭 兴 宇 １ ꎬ ２ １ . 河 北 师 范 大 学 资 源 与 环 境 科 学 学 院 ꎬ 河 北 石 家 庄 ０ ５ ０ ０ ２ ４ ２ . 河 北 省 环 境 演 变 与 生 态 建 设 实 验 室 ꎬ 河 北 石 家 庄 ０ ５ ０ ０ ２ ４ 摘 要 京 津 冀 地 区 大 气 Ｐ Ｍ ２ ５ 污 染 严 重 . 为 揭 示 区 域 Ｐ Ｍ ２ ５ 时 空 分 布 规 律 ꎬ 使 用 ２ ０ １ ３ ２ ０ １ ４ 年 河 北 省 地 面 站 点 Ｐ Ｍ ２ ５ 监 测 数 据 、 Ｍ Ｏ Ｄ Ｉ Ｓ Ａ Ｏ Ｄ 气 溶 胶 光 学 厚 度 遥 感 数 据 、 地 面 气 象 站 点 数 据 和 土 地 利 用 调 查 数 据 ꎬ 基 于 线 性 混 合 效 应 模 型 Ｌ Ｍ Ｅ ꎬ 建 立 了 ρ Ｐ Ｍ ２ ５ 时 空 变 化 与 Ａ Ｏ Ｄ 因 子 、 气 象 因 子 、 土 地 利 用 因 子 之 间 的 关 系 模 型 . 采 用 十 折 交 叉 验 证 法 对 模 型 精 度 进 行 检 验 ꎬ 并 利 用 计 算 得 到 的 校 正 因 子 〔 全 部 实 测 的 ρ Ｐ Ｍ ２ ５ 年 均 值 除 以 参 与 建 模 的 所 有 实 测 ρ Ｐ Ｍ ２ ５ 年 均 值 〕 纠 正 因 Ａ Ｏ Ｄ 非 随 机 性 缺 值 导 致 的 抽 样 偏 差 . 结 果 表 明 ① 河 北 省 区 域 模 拟 精 度 Ｒ ２ 决 定 系 数 为 ０ ８ ５ ꎬ 经 交 叉 验 证 后 Ｒ ２ 为 ０ ７ ７ ꎬ Ｒ Ｍ Ｓ Ｅ 均 方 根 误 差 和 Ｒ Ｐ Ｅ 相 对 预 测 误 差 分 别 为 １ ８ ２ ８ μ ｇ ∕ ｍ ３ 和 ２ ８ ６ ８ ％ . ② ρ Ｐ Ｍ ２ ５ 年 均 值 模 拟 结 果 的 校 正 因 子 范 围 为 １ ２ ４ ２ ０ ５ ꎬ 校 正 后 的 研 究 区 ρ Ｐ Ｍ ２ ５ 年 均 值 为 ８ ９ ８ ４ μ ｇ ∕ ｍ ３ ꎬ 与 实 际 监 测 数 据 相 近 . ③ ρ Ｐ Ｍ ２ ５ 空 间 分 布 呈 平 原 高 、 山 区 低 ꎬ 平 原 地 区 西 南 高 、 东 北 低 的 趋 势 . ④ ρ Ｐ Ｍ ２ ５ 与 Ａ Ｏ Ｄ 、 温 度 、 相 对 湿 度 呈 正 相 关 ꎬ 与 风 速 、 大 气 能 见 度 呈 负 相 关 . 研 究 显 示 ꎬ 线 性 混 合 效 应 模 型 能 有 效 对 ρ Ｐ Ｍ ２ ５ 进 行 时 空 变 化 模 拟 ꎬ 并 实 现 对 非 地 面 监 测 地 区 ρ Ｐ Ｍ ２ ５ 时 空 变 化 的 预 测 ꎬ 恰 当 的 预 测 因 子 组 合 和 模 型 校 正 有 助 于 模 型 预 测 精 度 的 提 升 . 关 键 词 Ｐ Ｍ ２ ５ ꎻ 气 溶 胶 光 学 厚 度 Ａ Ｏ Ｄ ꎻ 线 性 混 合 效 应 模 型 ꎻ 十 折 交 叉 验 证 ꎻ 时 空 变 化 ꎻ 河 北 省 中 图 分 类 号 Ｘ ５ １ ３ 文 章 编 号 １ ０ ０ １ ￣ ６ ９ ２ ９ ２ ０ １ ９ ０ ９ ￣ １ ５ ０ ０ ￣ １ ０ 文 献 标 志 码 Ａ Ｄ Ｏ Ｉ １ ０ １ ３ １ ９ ８ ∕ ｊ ｉ ｓ ｓ ｎ １ ０ ０ １ ￣ ６ ９ ２ ９ ２ ０ １ ９ ０ ６ １ ３ Ｓ ｐ ａ ｔ ｉ ａ ｌ ￣ Ｔ ｅ ｍ ｐ ｏ ｒ ａ ｌ Ｓ ｉ ｍ ｕ ｌ ａ ｔ ｉ ｏ ｎ ｏ ｆ Ｐ Ｍ ２ ５ Ｃ ｏ ｎ ｃ ｅ ｎ ｔ ｒ ａ ｔ ｉ ｏ ｎ ｉ ｎ Ｈ ｅ ｂ ｅ ｉ Ｐ ｒ ｏ ｖ ｉ ｎ ｃ ｅ ｂ ａ ｓ ｅ ｄ ｏ ｎ Ｌ ｉ ｎ ｅ ａ ｒ Ｍ ｉ ｘ ｅ ｄ Ｅ ｆ ｆ ｅ ｃ ｔ ｓ Ｍ ｏ ｄ ｅ ｌ Ｓ Ｕ Ｎ Ｃ ｈｅ ｎ ｇ １ ꎬ ２ ꎬ Ｗ Ａ Ｎ Ｇ Ｗ ｅ ｉ １ ꎬ ２ ∗ ꎬ Ｌ Ｉ Ｕ Ｆ ａ ｎｇ ｔ ｉ ａ ｎ １ ꎬ ２ ꎬ Ｇ Ｕ Ｏ Ｘ ｉ ｎｇ ｙ ｕ １ ꎬ ２ １ . Ｃ ｏ ｌ ｌ ｅ ｇ ｅ ｏ ｆ Ｒ ｅ ｓ ｏ ｕ ｒ ｃ ｅ ｓ ａ ｎ ｄ Ｅ ｎ ｖ ｉ ｒ ｏ ｎ ｍ ｅ ｎ ｔ Ｓ ｃ ｉ ｅ ｎ ｃ ｅ ｓ ꎬ Ｈ ｅ ｂ ｅ ｉ Ｎ ｏ ｒ ｍ ａ ｌ Ｕ ｎ ｉ ｖ ｅ ｒ ｓ ｉ ｔ ｙ ꎬ Ｓ ｈ ｉ ｊ ｉ ａ ｚ ｈ ｕ ａ ｎ ｇ ０ ５ ０ ０ ２ ４ ꎬ Ｃ ｈ ｉ ｎ ａ ２ . Ｋ ｅ ｙ Ｌ ａ ｂ ｏ ｒ ａ ｔ ｏ ｒ ｙ ｏ ｆ Ｅ ｎ ｖ ｉ ｒ ｏ ｎ ｍ ｅ ｎ ｔ ａ ｌ Ｅ ｖ ｏ ｌ ｖ ｅ ｍ ｅ ｎ ｔ ａ ｎ ｄ Ｅ ｃ ｏ ｌ ｏ ｇ ｉ ｃ ａ ｌ Ｃ ｏ ｎ ｓ ｔ ｒ ｕ ｃ ｔ ｉ ｏ ｎ ｏ ｆ Ｈ ｅ ｂ ｅ ｉ Ｐ ｒ ｏ ｖ ｉ ｎ ｃ ｅ ꎬ Ｓ ｈ ｉ ｊ ｉ ａ ｚ ｈ ｕ ａ ｎ ｇ ０ ５ ０ ０ ２ ４ ꎬ Ｃ ｈ ｉ ｎ ａ Ａ ｂ ｓ ｔ ｒ ａ ｃ ｔ Ｔ ｈ ｅ ａ ｔ ｍ ｏ ｓ ｐ ｈ ｅ ｒ ｉ ｃ ｐ ａ ｒ ｔ ｉ ｃ ｕ ｌ ａ ｔ ｅ ｐ ｏ ｌ ｌ ｕ ｔ ｉ ｏ ｎ ｉ ｓ ｓ ｅ ｒ ｉ ｏ ｕ ｓ ｉ ｎ ｔ ｈ ｅ Ｂ ｅ ｉ ｊ ｉ ｎ ｇ ￣ Ｔ ｉ ａ ｎ ｊ ｉ ｎ ￣ Ｈ ｅ ｂ ｅ ｉ Ｒ ｅ ｇ ｉ ｏ ｎ . Ｔ ｏ ｅ ｘ ｐ ｌ ｏ ｒ ｅ ｔ ｈ ｅ ｓ ｐ ａ ｔ ｉ ａ ｌ ￣ ｔ ｅ ｍ ｐ ｏ ｒ ａ ｌ ｄ ｉ ｓ ｔ ｒ ｉ ｂ ｕ ｔ ｉ ｏ ｎ ｏ ｆ ρ Ｐ Ｍ ２ ５ ꎬ ｂ ａ ｓ ｅ ｄ ｏ ｎ ｔ ｈ ｅ ｍ ｅ ａ ｓ ｕ ｒ ｅ ｄ ρ Ｐ Ｍ ２ ５ ꎬ Ｍ Ｏ Ｄ Ｉ Ｓ ａ ｅ ｒ ｏ ｓ ｏ ｌ ｏ ｐ ｔ ｉ ｃ ａ ｌ ｄ ｅ ｐ ｔ ｈ Ａ Ｏ Ｄ ꎬ ｍ ｅ ｔ ｅ ｏ ｒ ｏ ｌ ｏ ｇ ｉ ｃ ａ ｌ ａ ｎ ｄ ｌ ａ ｎ ｄ ｕ ｓ ｅ ｄ ａ ｔ ａ ｆ ｒ ｏ ｍ ２ ０ １ ３ ｔ ｏ ２ ０ １ ４ ꎬ ｗ ｅ ｄ ｅ ｖ ｅ ｌ ｏ ｐ ｅ ｄ ａ ｌ ｉ ｎ ｅ ａ ｒ ｍ ｉ ｘ ｅ ｄ ｅ ｆ ｆ ｅ ｃ ｔ ｍ ｏ ｄ ｅ ｌ ｔ ｏ ｍ ａ ｋ ｅ ｒ ｅ ｇ ｒ ｅ ｓ ｓ ρ Ｐ Ｍ ２ ５ ｍ ｅ ａ ｓ ｕ ｒ ｅ ｍ ｅ ｎ ｔ ｓ ｗ ｉ ｔ ｈ Ａ Ｏ Ｄ ꎬ ｍ ｅ ｔ ｅ ｏ ｒ ｏ ｌ ｏ ｇ ｉ ｃ ａ ｌ ａ ｎ ｄ ｌ ａ ｎ ｄ ｕ ｓ ｅ ｆ ａ ｃ ｔ ｏ ｒ ｓ ｉ ｎ Ｈ ｅ ｂ ｅ ｉ Ｐ ｒ ｏ ｖ ｉ ｎ ｃ ｅ . Ｔ ｈ ｅ ｎ ｔ ｅ ｎ ￣ ｆ ｏ ｌ ｄ ｃ ｒ ｏ ｓ ｓ ｖ ａ ｌ ｉ ｄ ａ ｔ ｉ ｏ ｎ ｗ ａ ｓ ｕ ｓ ｅ ｄ ｔ ｏ ｖ ａ ｌ ｉ ｄ ａ ｔ ｅ ｔ ｈ ｅ ａ ｃ ｃ ｕ ｒ ａ ｃ ｙ ｏ ｆ ｔ ｈ ｅ ｐ ｒ ｅ ｄ ｉ ｃ ｔ ｉ ｏ ｎ ｓ . Ｆ ｉ ｎ ａ ｌ ｌ ｙ ꎬ ｔ ｈ ｅ ｃ ｏ ｒ ｒ ｅ ｃ ｔ ｉ ｏ ｎ ｆ ａ ｃ ｔ ｏ ｒ ｓ ꎬ ｗ ｈ ｉ ｃ ｈ ｗ ｅ ｒ ｅ ｃ ａ ｌ ｃ ｕ ｌ ａ ｔ ｅ ｄ ｆ ｒ ｏ ｍ ｔ ｈ ｅ ｍ ｅ ａ ｓ ｕ ｒ ｅ ｄ ａ ｎ ｎ ｕ ａ ｌ ａ ｖ ｅ ｒ ａ ｇ ｅ ρ Ｐ Ｍ ２ ５ ｄ ｉ ｖ ｉ ｄ ｅ ｄ ｔ ｈ ｏ ｓ ｅ ｉ ｎ ｔ ｈ ｅ ｍ ｏ ｄ ｅ ｌ ｓ ꎬ ｗ ｅ ｒ ｅ ｅ ｍ ｐ ｌ ｏ ｙ ｅ ｄ ｔ ｏ ｃ ｏ ｒ ｒ ｅ ｃ ｔ ｂ ｉ ａ ｓ ｅ ｓ ｉ ｎ ｔ ｈ ｅ ｐ ｒ ｅ ｄ ｉ ｃ ｔ ｅ ｄ ρ Ｐ Ｍ ２ ５ ｃ ａ ｕ ｓ ｅ ｄ ｂ ｙ ｎ ｏ ｎ ｒ ａ ｎ ｄ ｏ ｍ ｍ ｉ ｓ ｓ ｉ ｎ ｇ ｏ ｆ Ａ Ｏ Ｄ . Ｔ ｈ ｅ ｒ ｅ ｓ ｕ ｌ ｔ ｓ ｓ ｈ ｏ ｗ ｅ ｄ ｔ ｈ ａ ｔ １ Ｔ ｈ ｅ Ｒ ２ ｏ ｆ ｔ ｈ ｅ ｍ ｉ ｘ ｅ ｄ ｅ ｆ ｆ ｅ ｃ ｔ ｓ ｍ ｏ ｄ ｅ ｌ ｆ ｏ ｒ Ｈ ｅ ｂ ｅ ｉ Ｐ ｒ ｏ ｖ ｉ ｎ ｃ ｅ ｗ ａ ｓ ０ ８ ５ . Ｔ ｈ ｅ ｃ ｒ ｏ ｓ ｓ ￣ ｖ ａ ｌ ｉ ｄ ａ ｔ ｉ ｏ ｎ Ｒ ２ ｗ ａ ｓ ０ ７ ７ ꎬ ｔ ｈ ｅ Ｒ ｏ ｏ ｔ Ｍ ｅ ａ ｎ Ｓ ｑ ｕ ａ ｒ ｅ Ｅ ｒ ｒ ｏ ｒ Ｒ Ｍ Ｓ Ｅ ａ ｎ ｄ Ｒ ｅ ｌ ａ ｔ ｉ ｖ ｅ Ｐ ｒ ｅ ｄ ｉ ｃ ｔ ｉ ｏ ｎ Ｅ ｒ ｒ ｏ ｒ Ｒ Ｐ Ｅ ｗ ｅ ｒ ｅ １ ８ ２ ８ μ ｇ ∕ ｍ ３ ａ ｎ ｄ ２ ８ ６ ８ ％ ꎬ ｒ ｅ ｓ ｐ ｅ ｃ ｔ ｉ ｖ ｅ ｌ ｙ . ２ Ｔ ｈ ｅ ｃ ｏ ｒ ｒ ｅ ｃ ｔ ｉ ｏ ｎ ｆ ａ ｃ ｔ ｏ ｒ ｓ ｏ ｆ ２ ０ １ ３ ￣ ２ ０ １ ４ ｒ ａ ｎ ｇ ｅ ｄ ｆ ｒ ｏ ｍ １ ２ ４ ｔ ｏ ２ ０ ５ . Ｔ ｈ ｅ ｃ ｏ ｒ ｒ ｅ ｃ ｔ ｅ ｄ ａ ｎ ｎ ｕ ａ ｌ ρ Ｐ Ｍ ２ ５ ｏ ｆ Ｈ ｅ ｂ ｅ ｉ Ｐ ｒ ｏ ｖ ｉ ｎ ｃ ｅ ｗ ａ ｓ ８ ９ ８ ４ μ ｇ ∕ ｍ ３ ꎬ ｗ ｈ ｉ ｃ ｈ ｗ ａ ｓ ｃ ｌ ｏ ｓ ｅ ｔ ｏ ｔ ｈ ｅ ｍ ｏ ｎ ｉ ｔ ｏ ｒ ｉ ｎ ｇ ｄ ａ ｔ ａ . ３ Ｆ ｏ ｒ ｔ ｈ ｅ ｓ ｐ ａ ｔ ｉ ａ ｌ ｐ ａ ｔ ｔ ｅ ｒ ｎ ꎬ ｔ ｈ ｅ Ｐ Ｍ ２ ５ ｃ ｏ ｎ ｃ ｅ ｎ ｔ ｒ ａ ｔ ｉ ｏ ｎ ｓ ｉ ｎ ｔ ｈ ｅ ｐ ｌ ａ ｉ ｎ ｓ ｗ ｅ ｒ ｅ ｈ ｉ ｇ ｈ ｅ ｒ ｔ ｈ ａ ｎ ｔ ｈ ｏ ｓ ｅ ｉ ｎ ｔ ｈ ｅ ｍ ｏ ｕ ｎ ｔ ａ ｉ ｎ ｏ ｕ ｓ ａ ｒ ｅ ａ ｓ ꎬ ｔ ｈ ｅ ρ Ｐ Ｍ ２ ５ ｉ ｎ ｔ ｈ ｅ ｓ ｏ ｕ ｔ ｈ ｗ ｅ ｓ ｔ ｅ ｒ ｎ ｐ ｌ ａ ｉ ｎ ｗ ｅ ｒ ｅ ｈ ｉ ｇ ｈ ｅ ｒ ｔ ｈ ａ ｎ ｔ ｈ ｏ ｓ ｅ ｉ ｎ ｔ ｈ ｅ ｎ ｏ ｒ ｔ ｈ ｅ ａ ｓ ｔ ｅ ｒ ｎ ｐ ａ ｒ ｔ . ４ Ｔ ｈ ｅ ρ Ｐ Ｍ ２ ５ ｈ ａ ｄ ｐ ｏ ｓ ｉ ｔ ｉ ｖ ｅ ｃ ｏ ｒ ｒ ｅ ｌ ａ ｔ ｉ ｏ ｎ ｗ ｉ ｔ ｈ ｔ ｈ ｅ Ａ Ｏ Ｄ ꎬ ｔ ｅ ｍ ｐ ｅ ｒ ａ ｔ ｕ ｒ ｅ ꎬ ｒ ｅ ｌ ａ ｔ ｉ ｖ ｅ ｈ ｕ ｍ ｉ ｄ ｉ ｔ ｙ ａ ｎ ｄ ｎ ｅ ｇ ａ ｔ ｉ ｖ ｅ ｃ ｏ ｒ ｒ ｅ ｌ ａ ｔ ｉ ｏ ｎ ｗ ｉ ｔ ｈ ｔ ｈ ｅ ｗ ｉ ｎ ｄ ｓ ｐ ｅ ｅ ｄ ａ ｎ ｄ ａ ｔ ｍ ｏ ｓ ｐ ｈ ｅ ｒ ｉ ｃ ｖ ｉ ｓ ｉ ｂ ｉ ｌ ｉ ｔ ｙ . Ａ ｌ ｌ ｔ ｈ ｅ ｓ ｅ ｒ ｅ ｓ ｕ ｌ ｔ ｓ ｓ ｕ ｇ ｇ ｅ ｓ ｔ ｅ ｄ ｔ ｈ ａ ｔ ｔ ｈ ｅ ｍ ｉ ｘ ｅ ｄ ｅ ｆ ｆ ｅ ｃ ｔ ｍ ｏ ｄ ｅ ｌ ａ ｌ ｌ ｏ ｗ ｓ ｏ ｎ ｅ ｔ ｏ ａ ｓ ｓ ｅ ｓ ｓ ｔ ｈ ｅ ｓ ｐ ａ ｔ ｉ ａ ｌ ￣ ｔ ｅ ｍ ｐ ｏ ｒ ａ ｌ ｄ ｉ ｓ ｔ ｒ ｉ ｂ ｕ ｔ ｉ ｏ ｎ ｏ ｆ ｔ ｈ ｅ ρ Ｐ Ｍ ２ ５ ｒ ｅ ｌ ｉ ａ ｂ ｌ ｙ ꎻ ｉ ｔ ｃ ａ ｎ ａ ｌ ｓ ｏ ｐ ｒ ｅ ｄ ｉ ｃ ｔ ｔ ｈ ｅ ｓ ｐ ａ ｔ ｉ ａ ｌ ￣ ｔ ｅ ｍ ｐ ｏ ｒ ａ ｌ ｐ ａ ｔ ｔ ｅ ｒ ｎ ｏ ｆ ｔ ｈ ｅ ρ Ｐ Ｍ ２ ５ ｉ ｎ ｎ ｏ ｎ ￣ ｍ ｅ ａ ｓ ｕ ｒ ｅ ｄ第 ９ 期 孙 成 等 基 于 线 性 混 合 效 应 模 型 的 河 北 省 Ｐ Ｍ ２ ５ 浓 度 时 空 变 化 模 型 研 究 ｒ ｅ ｇ ｉ ｏ ｎ ｓ . Ｐ ｒ ｏ ｐ ｅ ｒ ｃ ｏ ｍ ｂ ｉ ｎ ａ ｔ ｉ ｏ ｎ ｏ ｆ ｆ ｏ ｒ ｅ ｃ ａ ｓ ｔ ｉ ｎ ｇ ｆ ａ ｃ ｔ ｏ ｒ ｓ ａ ｎ ｄ ｍ ｏ ｄ ｅ ｌ ｃ ｏ ｒ ｒ ｅ ｃ ｔ ｉ ｏ ｎ ｆ ａ ｃ ｔ ｏ ｒ ａ ｒ ｅ ａ ｂ ｌ ｅ ｔ ｏ ｈ ｅ ｌ ｐ ｉ ｍ ｐ ｒ ｏ ｖ ｅ ｔ ｈ ｅ ｐ ｒ ｅ ｄ ｉ ｃ ｔ ｉ ｏ ｎ ａ ｃ ｃ ｕ ｒ ａ ｃ ｙ ｏ ｆ ｔ ｈ ｅ ｍ ｏ ｄ ｅ ｌ . Ｋ ｅ ｙ ｗ ｏ ｒ ｄ ｓ ρ Ｐ Ｍ ２ ５ ꎻ ａ ｅ ｒ ｏ ｓ ｏ ｌ ｏ ｐ ｔ ｉ ｃ ａ ｌ ｄ ｅ ｐ ｔ ｈ Ａ Ｏ Ｄ ꎻ ｌ ｉ ｎ ｅ ａ ｒ ｍ ｉ ｘ ｅ ｄ ｅ ｆ ｆ ｅ ｃ ｔ ｓ ｍ ｏ ｄ ｅ ｌ ꎻ ｔ ｅ ｎ ￣ ｆ ｏ ｌ ｄ ｃ ｒ ｏ ｓ ｓ ｖ ａ ｌ ｉ ｄ ａ ｔ ｉ ｏ ｎ ꎻ ｓ ｐ ａ ｔ ｉ ａ ｌ ￣ ｔ ｅ ｍ ｐ ｏ ｒ ａ ｌ ｄ ｉ ｓ ｔ ｒ ｉ ｂ ｕ ｔ ｉ ｏ ｎ ꎻ Ｈ ｅ ｂ ｅ ｉ Ｐ ｒ ｏ ｖ ｉ ｎ ｃ ｅ进 入 ２１ 世 纪 以 来 ꎬ 我 国 经 济 继 续 高 速 增 长 ꎬ 大 气 污 染 也 呈 明 显 加 重 趋 势 [ １ ￣ ３ ] ꎬ 特 别 是 河 北 省 及 其 所 在 的 京 津 冀 地 区 成 为 大 气 颗 粒 物 污 染 的 重 灾 区 . 为 此 ２０ １３ 年 以 来 ꎬ 国 家 实 施 了 最 严 厉 的 治 理 措 施 ꎬ 京 津 冀 地 区 大 气 环 境 质 量 有 所 好 转 ꎬ 但 治 理 压 力 依 然 较 大 [ ４ ￣ ５ ] . 长 期 以 来 ꎬ 学 术 界 针 对 大 气 颗 粒 物 污 染 时 空 变 化 定 量 研 究 进 行 了 不 断 的 探 索 和 尝 试 ꎬ 早 期 多 依 赖 于 研 究 区 内 有 限 的 Ｐ Ｍ ２ ５ 监 测 站 点 ꎬ 利 用 站 点 监 测 数 据 及 其 他 相 关 因 子 预 测 附 近 未 监 测 区 域 的 ρ Ｐ Ｍ ２ ５ [ ６ ] . 土 地 利 用 回 归 模 型 Ｌ Ｕ Ｒ 是 地 表 大 气 颗 粒 物 质 量 浓 度 空 间 分 布 研 究 中 影 响 比 较 广 泛 的 一 种 方 法 ꎬ 通 过 构 建 地 面 监 测 ρ Ｐ Ｍ ２ ５ 与 当 地 土 地 利 用 及 其 相 关 因 子 如 道 路 密 度 、 建 筑 用 地 比 率 、 人 口 密 度 等 之 间 的 回 归 关 系 ꎬ 进 而 模 拟 研 究 区 面 状 ρ Ｐ Ｍ ２ ５ [ ７ ￣ ８ ] . 但 由 于 研 究 区 内 土 地 利 用 状 况 短 期 内 几 乎 不 会 发 生 明 显 变 化 ꎬ 导 致 模 型 的 时 间 分 辨 率 较 低 ꎬ 无 法 模 拟 ρ Ｐ Ｍ ２ ５ 短 期 时 空 变 化 信 息 . 伴 随 新 型 传 感 器 和 算 法 的 产 生 ꎬ 卫 星 遥 感 逐 渐 成 为 一 个 强 有 力 的 工 具 ꎬ 用 来 对 不 同 空 间 尺 度 下 ρ Ｐ Ｍ ２ ５ 进 行 评 估 . 遥 感 反 演 产 品 大 气 Ａ Ｏ Ｄ ａ ｅ ｒ ｏ ｓ ｏ ｌ ｏ ｐ ｔ ｉ ｃ ａ ｌ ｄｅ ｐ ｔ ｈ ꎬ 气 溶 胶 光 学 厚 度 是 整 层 大 气 柱 中 气 溶 胶 的 消 光 系 数 ꎬ 与 ρ Ｐ Ｍ ２ ５ 之 间 存 在 较 好 且 复 杂 的 相 关 关 系 ꎬ 因 其 覆 盖 空 间 范 围 广 、 获 取 方 便 ꎬ 近 年 来 被 广 泛 应 用 于 ρ Ｐ Ｍ ２ ５ 的 模 拟 和 预 测 [ ９ ￣ １ ３ ] . 由 于 Ａ Ｏ Ｄ 数 据 可 以 利 用 卫 星 遥 感 获 得 ꎬ 每 日 数 据 进 行 实 时 更 新 ꎬ 以 Ａ Ｏ Ｄ 为 变 量 因 子 对 大 气 ρ Ｐ Ｍ ２ ５ 进 行 时 空 预 测 可 有 效 提 高 预 测 模 型 的 时 间 分 辨 率 [ １ ４ ￣ １ ６ ] ꎻ 同 时 ꎬ 可 以 极 大 改 善 早 期 模 型 过 分 依 赖 监 测 站 点 分 布 的 缺 陷 ꎬ 使 预 测 模 型 能 在 更 大 区 域 范 围 推 广 使 用 [ １ ７ ￣ １ ９ ] . 将 Ａ Ｏ Ｄ 参 数 引 入 预 测 模 型 模 拟 Ａ Ｏ Ｄ 与 Ｐ Ｍ ２ ５ 相 关 关 系 的 方 法 很 多 ꎬ 其 中 比 较 成 熟 的 是 Ｃ Ｔ Ｍ 大 气 化 学 传 输 模 型 . Ｃ Ｔ Ｍ 首 先 需 要 对 遥 感 Ａ Ｏ Ｄ 数 据 进 行 一 系 列 校 正 ꎬ 包 括 大 气 垂 直 订 正 和 相 对 湿 度 订 正 等 ꎬ 去 除 因 气 象 因 素 导 致 的 数 据 偏 差 ꎬ 并 利 用 Ｃ Ｔ Ｍ 模 拟 订 正 后 Ａ Ｏ Ｄ 数 据 与 ρ Ｐ Ｍ ２ ５ 数 据 的 比 例 关 系 ꎬ 从 而 估 算 研 究 区 域 ρ Ｐ Ｍ ２ ５ ꎬ 该 模 型 在 北 美 地 区 的 预 测 结 果 与 实 际 值 相 关 性 为 ０ ７７ [ ２ ０ ￣ ２ ３ ] ꎬ 模 拟 精 度 较 高 . 该 方 法 的 优 势 在 于 不 需 要 ρ Ｐ Ｍ ２ ５ 数 据 就 可 以 进 行 模 型 模 拟 ꎬ 扩 大 了 模 型 的 应 用 范 围 . 但 Ｃ Ｔ Ｍ 构 建 过 程 对 于 数 据 要 求 较 高 ꎬ 需 要 污 染 物 排 放 清 单 和 气 溶 胶 粒 径 、 类 型 、 种 类 、 消 光 系 数 的 垂 直 结 构 等 详 细 信 息 ꎬ 因 此 数 据 精 度 对 模 拟 结 果 影 响 较 大 . 近 年 来 大 量 研 究 应 用 线 性 混 合 效 应 模 型 Ｌ Ｍ Ｅ 研 究 大 气 颗 粒 物 污 染 情 况 ꎬ 将 该 模 型 与 遥 感 Ａ Ｏ Ｄ 数 据 结 合 ꎬ 并 引 入 风 速 、 相 对 湿 度 等 气 象 参 数 和 土 地 利 用 参 数 . 与 Ｃ Ｔ Ｍ 相 比 ꎬ 线 性 混 合 效 应 模 型 由 于 增 加 了 模 型 的 随 机 效 应 ꎬ 具 有 数 据 需 求 相 对 简 单 且 模 拟 精 度 较 高 的 优 势 [ ２ ４ ￣ ２ ６ ] . 该 研 究 以 Ａ Ｏ Ｄ 参 数 、 气 象 参 数 及 土 地 利 用 参 数 作 为 预 测 因 子 ꎬ 构 建 河 北 省 ２０ １３ ２ ０１ ４ 年 逐 日 线 性 混 合 效 应 模 型 ꎬ 模 拟 分 析 研 究 区 域 ρ Ｐ Ｍ ２ ５ 时 空 变 化 规 律 ꎬ 以 期 为 大 气 污 染 治 理 提 供 科 学 依 据 . １ 材 料 与 方 法 １ １ 研 究 区 概 况 河 北 省 环 抱 北 京 市 、 天 津 市 ꎬ 地 处 ３６ ° ０５ ′ Ｎ ４ ２ ° ４０ ′ Ｎ 、 １１３ ° ２７ ′ Ｅ １１ ９ ° ５０ ′ Ｅ 之 间 ꎬ 东 邻 渤 海 ꎬ 西 倚 太 行 山 脉 ꎬ 北 靠 内 蒙 古 高 原 ꎬ 南 接 豫 鲁 平 原 . 属 于 温 带 大 陆 性 季 风 气 候 ꎬ 四 季 分 明 ꎬ 降 水 量 东 南 多 、 西 北 少 . 地 貌 类 型 复 杂 多 样 ꎬ 高 原 、 山 地 、 丘 陵 、 平 原 、 盆 地 类 型 齐 全 ꎬ 地 势 西 北 高 、 东 南 低 ꎬ 并 由 西 北 向 东 南 倾 斜 . 由 于 区 域 内 西 、 北 两 面 环 山 ꎬ 阻 断 了 相 关 气 流 扰 动 平 原 ꎬ 不 利 于 大 气 污 染 物 有 效 扩 散 ꎬ 并 且 受 高 气 压 影 响 下 的 冬 季 气 流 常 导 致 逆 温 现 象 频 发 ꎬ 加 重 了 平 原 地 区 的 污 染 程 度 ꎬ 使 其 成 为 我 国 大 气 污 染 的 重 灾 区 . １ ２ 基 础 数 据 １ ２ １ Ａ Ｏ Ｄ 数 据 Ｍ Ｏ Ｄ Ｉ Ｓ 是 Ｎ Ａ ＳＡ 美 国 国 家 航 空 航 天 局 地 球 观 测 系 统 Ｅ Ｏ Ｓ 中 较 为 先 进 的 观 测 仪 器 ꎬ 广 泛 应 用 于 大 气 气 溶 胶 监 测 领 域 . Ａ Ｏ Ｄ 是 Ｍ Ｏ Ｄ Ｉ Ｓ 的 基 本 产 品 之 一 ꎬ 该 研 究 采 用 Ｎ Ａ ＳＡ 网 站 ｈｔ ｔ ｐ ∕ ∕ ｌ ａ ｄｓ ｗ ｅ ｂ . ｍ ｏ ｄａ ｐｓ . ｅ ｏ ｓ ｄｉ ｓ . ｎａ ｓ ａ . ｇ ｏ ｖ 发 布 的 ２０１ ３ 年 １ 月 １ 日 ２ ０ １４ 年 １ ２ 月 ３１ 日 的 Ｍ Ｏ Ｄ ０４ ＿ Ｌ ２ Ｃ ６ 产 品 ꎬ 时 间 分 辨 率 为 １ ｄ 、 空 间 分 辨 率 为 １ ０ ｋ ｍ １ ０ ｋ ｍ . 针 对 Ｍ Ｏ Ｄ Ｉ Ｓ Ａ Ｏ Ｄ 数 据 使 用 要 求 ꎬ 该 研 究 首 先 对 原 始 数 据 进 行 裁 剪 、 拼 接 处 理 ꎬ 同 时 由 于 Ａ Ｏ Ｄ 像 素 中 心 每 日 不 断 发 生 变 化 ꎬ 因 此 对 其 进 行 空 间 几 何 校 正 ꎬ 使 得 区 域 内 每 日 Ａ Ｏ Ｄ 格 网 中 心 能 够 完 全 重 合 ꎬ 最 终 获 得 覆 盖 河 北 省 的 Ａ Ｏ Ｄ 栅 格 数 据 . １ ２ ２ ρ Ｐ Ｍ ２ ５ 监 测 数 据 河 北 省 ２０ １３ ２ ０１ ４ 年 ρ Ｐ Ｍ ２ ５ 监 测 数 据 来 源 于 全 国 城 市 空 气 质 量 实 时 发 布 平 台 ｈ ｔ ｔ ｐ ∕ ∕ １ ０ ５ １ 环 境 科 学 研 究 第 ３ ２ 卷 １０ ６ ３ ７ ２ ０ ８ ２ ３３ ２００ ３５ ꎬ 从 该 平 台 可 获 得 河 北 省 辖 区 石 家 庄 市 、 保 定 市 、 邢 台 市 、 邯 郸 市 等 １１ 个 城 市 共 ５３ 个 监 测 站 点 位 置 及 其 ρ Ｐ Ｍ ２ ５ 小 时 数 据 ꎬ 用 于 计 算 ρ Ｐ Ｍ ２ ５ 日 均 值 . 根 据 监 测 站 点 的 经 纬 度 ꎬ 将 地 面 站 点 日 均 ρ Ｐ Ｍ ２ ５ 数 据 匹 配 到 与 它 们 相 对 应 的 Ａ Ｏ Ｄ 栅 格 中 . Ａ Ｏ Ｄ 栅 格 单 元 中 存 在 １ 个 或 多 个 Ｐ Ｍ ２ ５ 监 测 站 点 的 被 命 名 为 监 测 格 网 ꎬ 其 中 ꎬ 包 含 多 个 监 测 站 点 的 监 测 格 网 值 按 照 监 测 格 网 中 所 有 监 测 站 点 的 平 均 值 进 行 计 算 . 最 终 ꎬ 得 到 研 究 区 内 监 测 格 网 共 ３２ 个 . １ ２ ３ 气 象 数 据 气 象 数 据 主 要 源 于 中 国 气 象 数 据 共 享 网 ｈｔ ｔ ｐｓ ∕ ∕ ｄ ａ ｔ ａ . ｃ ｍ ａ . ｃ ｎ ꎬ 获 得 河 北 省 ２０１ ３ ２ ０ １４ 年 各 气 象 站 点 气 象 数 据 ꎬ 包 括 气 温 、 相 对 湿 度 、 风 速 、 能 见 度 等 . 但 由 于 该 数 据 类 型 均 为 站 点 点 状 数 据 ꎬ 难 以 与 所 有 监 测 格 网 进 行 数 据 空 间 匹 配 ꎬ 需 要 对 气 象 站 点 数 据 进 行 空 间 插 值 ꎬ 以 获 得 研 究 区 分 辨 率 为 １０ ｋｍ １ ０ ｋ ｍ 气 象 面 状 数 据 . 利 用 Ａ ｒ ｃ Ｇ Ｉ Ｓ １ ０ ２ 软 件 对 各 气 象 因 子 进 行 空 间 插 值 ꎬ 插 值 数 据 源 包 括 河 北 省 １ ４２ 个 气 象 站 点 及 其 周 边 地 区 北 京 市 ꎬ 天 津 市 ꎬ 山 西 省 阳 泉 市 、 大 同 市 、 长 治 市 ꎬ 河 南 省 安 阳 市 ꎬ 山 东 德 州 市 ꎬ 内 蒙 古 自 治 区 赤 峰 市 １１ 个 气 象 站 点 的 日 均 值 数 据 ꎬ 站 点 共 计 １５３ 个 . 将 所 有 站 点 数 据 分 为 测 试 集 和 验 证 集 两 部 分 ꎬ 其 中 验 证 集 占 气 象 站 点 的 １５ ％ ꎬ 目 的 在 于 对 空 间 插 值 结 果 进 行 精 度 验 证 . 气 温 、 相 对 湿 度 、 风 速 、 能 见 度 等 气 象 因 子 插 值 方 法 经 过 对 比 筛 选 ꎬ 最 终 确 定 能 见 度 因 子 为 Ｉ Ｄ Ｗ 插 值 ꎬ 其 余 选 择 普 通 克 里 金 插 值 . 经 过 验 证 集 数 据 检 验 ꎬ 插 值 结 果 Ｒ ２ 范 围 为 ０ ６９ ０ ９９ ꎬ 符 合 数 据 使 用 要 求 . １ ２ ４ 土 地 利 用 数 据 土 地 利 用 参 数 可 以 较 好 地 表 征 地 面 大 气 颗 粒 物 的 排 放 强 度 ꎬ 一 般 来 说 建 筑 用 地 Ｐ Ｍ ２ ５ 排 放 强 度 大 于 林 地 、 草 地 . 该 研 究 采 用 空 间 分 辨 率 为 ３ ０ ｍ 的 ２０ １ ５ 年 河 北 省 土 地 利 用 现 状 图 来 代 表 ２ ０１ ３ ２ ０ １４ 年 的 土 地 利 用 状 况 ꎬ 主 要 提 取 研 究 区 栅 格 单 元 建 筑 用 地 比 率 和 林 地 覆 盖 率 ２ 个 参 数 ꎬ 数 据 提 取 在 Ａ ｒ ｃ Ｇ Ｉ Ｓ １０ ２ 软 件 中 进 行 . １ ２ ５ 数 据 提 取 与 匹 配 由 于 各 类 数 据 来 源 不 同 ꎬ 其 数 据 类 型 和 时 空 分 辨 率 也 并 不 一 样 ꎬ 如 Ａ Ｏ Ｄ 、 土 地 利 用 数 据 分 别 为 空 间 分 辨 率 １０ ｋ ｍ １０ ｋｍ 和 ３０ ｍ ３０ ｍ 的 栅 格 数 据 ꎬ 而 Ｐ Ｍ ２ ５ 和 气 象 要 素 则 为 研 究 区 各 监 测 站 点 的 点 数 据 ꎻ Ａ Ｏ Ｄ 为 每 日 卫 星 过 境 时 的 瞬 时 数 据 ꎬ ρ Ｐ Ｍ ２ ５ 为 小 时 数 据 ꎬ 虽 然 Ａ Ｏ Ｄ 瞬 时 数 据 与 ρ Ｐ Ｍ ２ ５ 小 时 数 据 的 相 关 性 更 高 ꎬ 但 不 能 实 现 ρ Ｐ Ｍ ２ ５ 与 Ａ Ｏ Ｄ 关 系 的 全 天 覆 盖 . 研 究 [ ２ ７ ] 表 明 ꎬ ρ Ｐ Ｍ ２ ５ 小 时 数 据 与 日 均 值 数 据 之 间 存 在 高 度 相 关 性 ꎬ 因 此 建 立 Ａ Ｏ Ｄ 与 ρ Ｐ Ｍ ２ ５ 日 均 值 的 匹 配 关 系 ꎬ 既 可 满 足 模 型 精 度 要 求 又 可 实 现 全 天 覆 盖 . 依 据 模 型 数 据 使 用 要 求 ꎬ 需 对 不 同 类 型 数 据 做 统 一 化 处 理 并 进 行 数 据 匹 配 ꎬ 得 到 一 一 对 应 的 数 据 组 . 将 １ ２ ２ 节 中 提 到 的 Ａ Ｏ Ｄ 监 测 格 网 作 为 基 准 格 网 ꎬ 对 各 项 基 础 数 据 进 行 单 独 提 取 . 对 于 ρ Ｐ Ｍ ２ ５ 数 据 ꎬ 将 单 个 Ａ Ｏ Ｄ 监 测 格 网 中 所 有 Ｐ Ｍ ２ ５ 站 点 同 一 天 监 测 数 据 的 平 均 值 赋 值 给 该 格 网 ꎻ 研 究 区 气 象 站 点 日 均 值 数 据 已 经 插 值 为 覆 盖 整 个 研 究 区 的 栅 格 数 据 ꎬ 根 据 Ｐ Ｍ ２ ５ 监 测 站 点 经 纬 度 ꎬ 对 插 值 得 到 的 栅 格 数 据 进 行 气 象 数 据 提 取 ꎬ 得 到 与 监 测 站 点 所 在 格 网 相 对 应 的 气 象 要 素 值 ꎻ 土 地 利 用 数 据 则 直 接 在 基 准 格 网 内 计 算 得 到 每 个 监 测 格 网 的 建 筑 用 地 比 率 与 林 地 覆 盖 率 . 处 理 后 获 得 的 ρ Ｐ Ｍ ２ ５ 为 日 均 值 数 据 ꎬ Ａ Ｏ Ｄ 为 日 值 数 据 ꎬ 气 象 要 素 为 日 均 值 数 据 ꎬ 因 土 地 利 用 数 据 更 新 较 慢 ꎬ 以 土 地 利 用 现 状 数 据 代 表 年 均 值 数 据 . 将 各 项 数 据 提 取 结 果 按 照 监 测 格 网 的 Ｉ Ｄ 和 日 期 进 行 数 据 匹 配 ꎬ 并 去 除 Ａ Ｏ Ｄ 缺 值 部 分 ꎬ 最 终 得 到 研 究 区 包 含 ρ Ｐ Ｍ ２ ５ 在 内 的 所 有 参 数 匹 配 数 据 共 ２ １ ４９ 组 . １ ３ 模 型 构 建 与 验 证 １ ３ １ 线 性 混 合 效 应 模 型 构 建 线 性 混 合 效 应 模 型 是 既 包 含 固 定 效 应 又 包 含 随 机 效 应 的 统 计 模 型 ꎬ 其 中 ꎬ 固 定 效 应 表 示 影 响 因 子 对 Ｐ Ｍ ２ ５ 的 平 均 影 响 状 态 ꎬ 而 随 机 效 应 则 用 于 解 释 ρ Ｐ Ｍ ２ ５ 与 Ａ Ｏ Ｄ 及 其 他 影 响 因 子 之 间 的 日 变 化 关 系 ꎬ 以 随 机 截 距 和 随 机 斜 率 的 形 式 表 征 . 研 究 [ ２ ８ ] 表 明 ꎬ Ａ Ｏ Ｄ 与 ρ Ｐ Ｍ ２ ５ 之 间 的 线 性 关 系 不 仅 存 在 一 定 的 空 间 差 异 ꎬ 更 重 要 的 是 由 于 每 天 受 到 不 同 条 件 温 度 、 风 速 、 污 染 物 排 放 强 度 等 的 影 响 ꎬ ρ Ｐ Ｍ ２ ５ 与 Ａ Ｏ Ｄ 之 间 的 关 系 每 日 均 发 生 变 化 . 该 研 究 在 构 建 Ａ Ｏ Ｄ 与 ρ Ｐ Ｍ ２ ５ 线 性 混 合 效 应 模 型 的 基 础 上 ꎬ 加 入 气 象 参 数 和 土 地 利 用 参 数 ꎬ 来 校 正 Ａ Ｏ Ｄ 与 ρ Ｐ Ｍ ２ ５ 之 间 的 日 变 化 和 空 间 变 化 线 性 关 系 ꎬ 其 中 ꎬ 土 地 利 用 因 子 短 期 内 几 乎 不 随 时 间 变 化 ꎬ 直 接 作 为 空 间 效 应 处 理 . 线 性 混 合 效 应 模 型 Ｐ ｉ ｊ ＝ α ＋ μ ｊ ＋ β １ ＋ ｖ ｊ Ａ ｉ ｊ ＋ β ２ ＋ ｋ ｊ Ｔ ｉ ｊ ＋ β ３ ＋ ｐ ｊ Ｗ ｉ ｊ ＋ β ４ Ｒ ｉ ｊ ＋ β ５ Ｖ ｉ ｊ ＋ β ６ Ｃ ｉ ＋ β ７ Ｆ ｉ ＋ ｓ ｉ ＋ ε ｉ ｊ μ ｊ ꎬ ｖ ｊ ꎬ ｋ ｊ ꎬ ｐ ｊ Ｎ [ ０ ꎬ ０ ꎬ ０ ꎬ ０ ꎬ Ｂ ] １ 式 中 Ｐ ｉ ｊ 为 第 ｉ 个 监 测 格 网 在 第 ｊ 天 的 ρ Ｐ Ｍ ２ ５ ꎬ ２ ０ ５ １第 ９ 期 孙 成 等 基 于 线 性 混 合 效 应 模 型 的 河 北 省 Ｐ Ｍ ２ ５ 浓 度 时 空 变 化 模 型 研 究 μ ｇ ∕ ｍ ３ ꎻ Ａ ｉ ｊ 为 第 ｉ 个 监 测 格 网 在 第 ｊ 天 的 Ａ Ｏ Ｄ 值 ꎻ Ｔ ｉ ｊ 为 第 ｉ 个 监 测 格 网 在 第 ｊ 天 的 气 温 ꎬ ℃ ꎻ Ｗ ｉ ｊ 为 第 ｉ 个 监 测 格 网 在 第 ｊ 天 的 风 速 ꎬ ｍ ∕ ｓ ꎻ Ｒ ｉ ｊ 为 第 ｉ 个 监 测 格 网 在 第 ｊ 天 的 相 对 湿 度 ꎬ ％ ꎻ Ｖ ｉ ｊ 为 第 ｉ 个 监 测 格 网 在 第 ｊ 天 的 能 见 度 ꎬ ｋｍ ꎻ ꎻ Ｃ ｉ 为 第 ｉ 个 监 测 格 网 的 建 筑 用 地 比 率 ꎬ ％ ꎻ Ｆ ｉ 为 第 ｉ 个 监 测 格 网 的 林 地 覆 盖 率 ꎬ ％ ꎻ α 和 μ ｊ 分 别 为 模 型 固 定 截 距 和 随 机 截 距 ꎻ β １ β ７ 分 别 为 各 自 变 量 参 数 的 固 定 效 应 斜 率 ꎻ ｖ ｊ 为 第 ｊ 天 Ａ Ｏ Ｄ 的 随 机 效 应 斜 率 ꎻ ｋ ｊ 为 第 ｊ 天 气 温 的 随 机 效 应 斜 率 ꎻ ｐ ｊ 为 第 ｊ 天
展开阅读全文