七月的华北午后短时强降水说来就来。预报员最怕的往往不是已经发展成型的雷暴而是那些在山前刚刚冒头、还没被雷达捕捉到第一丝回波的对流胚。我做卫星气象数据应用这几年被问得最多的一个问题就是风云卫星在追踪对流这件事上到底能比雷达多看见什么这篇系列第七章我就把长期积累的实践和盘托出重点讲风云四号静止卫星如何在分钟级尺度上盯住对流从初生、合并、增强到溃散的完整一生以及处理TBB云顶亮温、识别对流初生、跟踪中尺度对流系统MCS时的操作流程和踩坑教训。文章面向气象预报员、做卫星遥感数据开发的同行、航空气象和新能源功率预测方向的朋友也欢迎有一定基础的气象爱好者追下来。1. 风云系列里能干对流追踪的主力装备静止卫星的分钟级凝视1.1 风云四号A/B星与AGRI把镜头钉在同一个地方对流追踪这件事最关键的不是空间分辨率有多高而是时间分辨率能不能跟上天气演变的节奏。雷暴从出现到猛烈发展往往只需要半小时到几小时用一两个小时才扫一次的卫星去看基本等于看回放。风云四号A星和B星都运行在地球静止轨道大约3.6万公里高度相对地面钉在一个固定经度上空可以持续凝视同一片区域。这两颗星上的核心载荷AGRI先进静止轨道辐射成像仪有14到15个通道覆盖可见光、近红外和红外窗区。业务上最常用的是这么几个0.65微米可见光通道空间分辨率约0.5到1公里主要看云顶纹理10.8微米红外窗区通道空间分辨率约2到4公里直接遥感云顶温度还有6.25微米和7.1微米左右的水汽通道用来观察上对流层的水汽输送。对于对流追踪这种场景4公里分辨率的红外通道已经足够因为我们要追踪的目标——冷云罩和强对流中心尺度通常在几十到几百公里。AGRI的观测模式需要单独说。常规的全圆盘扫描15分钟给一景这能满足大尺度监测但真正做短临预报我一般切到中国区域5分钟一次的区域扫描模式如果针对某个重点天气过程还可以用局部加密模式做到1分钟一次。一分钟级别的数据已经能比较细腻地反映云顶快速隆起的细节。1.2 为什么不能指望极轨卫星重访周期是硬伤有人可能会问风云三号极轨卫星通道更多还有微波仪器是不是更适合看对流这里要泼一盆冷水。风云三号这类极轨卫星确实有优势——微波通道可以穿透云层看到云体内部的结构而且全球覆盖。但它的硬伤是重访周期同一片区域一天往往只能过境一到两次等于给对流拍快照。对流是快速演变的15分钟后可能就从阵雨变成冰雹云极轨卫星根本等不起。所以在业务实践中我的原则是对流追踪的主干数据必须来自静止卫星风云三号这种极轨平台用来做全球对流分布的日尺度拼图、或者对特定过程做路过时的加密观测补充。两种轨道各有价值但不能互相替代。1.3 风云二号老家伙们还能发挥什么余热风云二号系列的观测通道少只有5个红外分辨率大约5公里扫描频率基本是半小时到一小时一次。放到现在的短临业务里确实有点吃力。但它有一个不可替代的价值——时间序列长新一代卫星的定量算法在做长序列检验时常常需要回溯历史数据风云二号能提供早年间的连续记录。我自己的习惯是风云二号拿来拉长线气候统计可以拿来做逐分钟的对流追踪就不现实了。2. 把看不见的对流变成看得见的数字云顶亮温与对流强度的底层逻辑2.1 三个通道各看什么东西卫星不是在拍照而是在测量不同波长的辐射。要读懂卫星图像先得搞明白每个通道背后的物理意义。可见光通道0.65微米反射的是太阳光晴天暗、云顶亮适合看云顶纹理。对流发展旺盛的时候云顶常常出现明显的上凸结构纹理变得粗糙像菜花一样翻涌这种细节在可见光云图上特别清楚。但一到晚上没有太阳光可见光通道就彻底失灵。红外窗区通道10.8微米接收的是云顶自身发射的红外辐射。根据普朗克辐射定律物体温度越低向外辐射的能量就越少卫星反演出来的亮温也越低。所以红外云图上颜色越亮亮温越低的位置云顶温度越低。水汽通道不太直观它测的是大气中上层水汽发出的辐射。对流强烈发展时水汽被抬升到对流层上部在图上会形成亮区而干空气下沉的地方会形成暗区。更重要的是水汽通道能清楚看到卷云砧的扩散方向这对判断对流系统未来的移动趋势很有帮助。2.2 为什么TBB越低往往意味着对流越猛这是卫星对流识别最核心的逻辑链对流越强上升气流越猛云顶就被顶得越高越接近对流层顶温度越低。所以云顶亮温的高低可以近似看作对流强度的标尺。我业务里常用的几个经验阈值判据阈值含义一般深对流TBB ≤ -32℃约241K云顶已进入对流层中上层MCS冷云罩TBB ≤ -52℃约221K达到中尺度对流系统统计标准强核/过冲TBB ≤ -70℃约203K云顶接近对流层顶甚至穿破注意这些阈值不是死的。夏天华北上空对流层顶温度可能低到-70℃上下这时-52℃只能算比较活跃的对流但冬天华南上空对流层顶偏低-52℃已经是非常剧烈的信号。所以看阈值的时候最好结合当天当地的探空资料了解对流层顶高度再判断TBB的绝对值意味着什么。2.3 从云顶冷到地面有灾中间还差好几步很多人一看到TBB降到-70℃以下就紧张以为地面马上要有大冰雹或者狂风。这个推断链条其实很长TBB低只说明云顶高、上升气流强至于地面到底下不下强降水、有没有大冰雹还要看云体内部的水分含量、0℃层高度、垂直风切变以及低层出流和入流配置。我见过不少例子——卫星云图上冷云罩冻得发白面积巨大但地面只是下了点零星小雨。原因很简单云砧是高空风把卷云往下游吹出来的云砧区域看着冷实际地面并没有对应降水。反过来有些强降水云系云顶并不极端冷因为暖云降水过程不需要云顶长到很高。所以我做预警时卫星TBB只用来判断对流强度趋势地面灾情判断一定要结合雷达回波、闪电定位和地面自动站实况这条经验后面还会展开讲。3. 对流初生CI识别抢在雷暴爆发前半小时的具体做法3.1 对流初生为什么难抓对流初生Convective InitiationCI阶段的云体小、变化幅度有限在卫星图像上相当不明显。刚开始对流的浅积云云顶也就五六公里高亮温只比地面气温低十几度在红外云图上根本没有触目惊心的效果。而等它发展到云顶-40℃以下、在云图上清清楚楚的时候地面回波可能已经很强了预警提前量就没了。还有一个客观困难卫星从头顶往下看看到的是云的天花板。对流初生阶段通常伴随着低层的辐合线、地形抬升、海陆风锋、冷池出流边界等因素这些启动机制在卫星云图上要么看不清楚要么只能看到次生的云线特征。所以做CI识别本质上是找那些正在快速长高、变冷的云块而不是找已经成熟的雷暴单体。3.2 我日常业务里用的四类CI判据先说结论不是单一判据就能定而是几个判据组合着看互相印证。第一个是云顶冷却速率。对流云顶如果正在快速上升温度就会快速下降。我一般统计过去15分钟TBB的变化量如果降温幅度超过4K而且连续两个时次保持这个趋势就认为这个云块处于活跃发展阶段。4K这个量级是从实测经验里来的太小的降温幅度传感器噪声都可能造成干扰。第二个是分裂窗差Split Window Difference就是8.5微米亮温减10.8微米亮温。这个差值和云顶粒子的相态、大小有关。对流发展初期云滴小、上升气流把过冷水滴带到更高处分裂窗差通常会出现明显变化。在业务产品里这个指标常常用来区分云顶粒子正在长大变冰的过程。第三个是可见光云顶隆起。白天如果看着一团亮白云块快速向上冒头、在云图上形成圆润的半球状凸起说明上升气流在加速。这个判断方法最直观但需要人眼盯图而且只适用于白天。第四个是水汽通道的暗带或干侵入配合。对流要强烈发展不光底下要有暖湿空气抬升中高层还要有干空气或强的辐散条件配合。水汽图上如果看到云团一侧有暗带切入往往预示着对流会组织化加强。3.3 用风云四号L1数据跑一个CI检测流程说到实操我给一个用Python处理风云四号AGRI L1级数据、快速识别CI区域的流程。数据可以从风云卫星遥感数据服务网下载到NetCDF格式文件核心逻辑并不复杂。import xarray as xr import numpy as np # 读取当前时次和15分钟前时次的NetCDF文件 ds_now xr.open_dataset(FY4A_AGRI_L1_NC_20240801_0600.nc) ds_prev xr.open_dataset(FY4A_AGRI_L1_NC_20240801_0545.nc) # 提取10.8微米通道TBB注意单位转换为开尔文 tbb_now ds_now[TBB].sel(band10.8) # 单位K tbb_prev ds_prev[TBB].sel(band10.8) # 计算15分钟亮温变化 delta_tbb tbb_now - tbb_prev # 冷却率判据15分钟降温4K且当前云顶已经达到深对流标准-32℃即241.15K ci_mask (delta_tbb -4.0) (tbb_now 241.15) # 如果某块区域像元数超过阈min_pixels标记一个CI候选区域 min_pixels 50 labeled np.zeros_like(ci_mask, dtypeint) # 这里省略连通域标记细节可以用scipy.ndimage.label实现 from scipy.ndimage import label as cc_label labeled, num_features cc_label(ci_mask) for i in range(1, num_features 1): pixels np.sum(labeled i) if pixels min_pixels: print(fCI候选区域编号{i}像元数{pixels}最大降温{delta_tbb.where(labeled i).min().values:.2f}K)实际业务中我还会把这个CI标记叠加在可见光云图和水汽图上人工二次确认。毕竟算法的判据比较机械边界层辐合线附近有时候因为云图配准问题会出现假的降温信号。用这套方法提前量大概能做到30到60分钟。这个时间窗口对短临预警来说非常宝贵足够通知一线值守人员加密观测或者给气象服务部门发出一级提示。4. 对流连片后的MCS追踪算法面积重叠法的实现与调参4.1 追踪对象从单核变成冷云罩对流发展到一定规模单体开始合并云顶往往连成一大片中间可能同时存在好几个强中心。这个时候再用某个强中心点来做追踪就会很痛苦——上一时次是两个核下一时次变成三个核再过一时次又合并成一个盯着单个核追踪的话轨迹会频繁跳变。我的做法是改追冷云罩也就是把TBB低于某个阈值的区域整体当作一个目标。最常用的阈值是-52℃这个标准在MCS的统计定义里也经常出现。为什么选-52℃因为它代表云顶高度已经相当高基本排除普通层状云能抓住真正有组织化的深对流系统。冷云罩的好处是稳定——云团分裂合并过程中冷云罩的形状会有变化但它作为目标实体是连续的适合做生命周期追踪。4.2 面积重叠法的四个关键步骤面积重叠法的思路其实很朴素上一时次的目标和当前时次的目标如果它们的几何区域大部分重叠就认为是同一个对流系统在移动演变。具体分四步。第一步阈值分割。把红外TBB数据转成二值图TBB低于阈值的像元记为1其余记0。这一步要特别注意数据里的缺测值和坏像元不能把缺测区当成对流空洞。第二步连通域标记。对二值图做连通域分析把彼此相连的目标像元归成一个个独立的对象。常用的8邻域标记也就是上下左右和四个对角方向相邻都算同一块。这里会碰到一个常见问题两个靠得很近的云团有可能被连成一个对象需要根据面积规模或TBB最小值的局部结构做二次拆分。第三步相邻时次目标匹配。对于当前时次的每个对象计算它与上一时次每个对象的重叠面积。匹配原则是找重叠面积占当前对象面积比例最大的上一时次对象同时要求这个重叠率IoU超过一定阈值。阈值常规取0.3到0.5。取太低会把相邻的独立系统硬拼在一起取太高又会在云团快速扩张时漏匹配。第四步处理分裂与合并。实际过程中MCS经常分裂成两个或者两个合并成一个。简单的一对一匹配算法处理不了这种情况。我的经验是如果当前一个对象与上一时次多个对象重叠判定为合并这个新对象的属性用上一时次各部分的面积加权平均来初始化如果一个上一时次对象对应当前多个对象判定为分裂新对象的生命周期继承自母体后续再各自独立追踪。追踪结果会输出一条完整的生命周期记录包括每个时次的中心经纬度、冷云罩面积、最低TBB、冷云罩在前后时次的增长率等。这些参数是判断MCS处于发展、成熟还是消散阶段的重要依据。4.3 调参过程中容易忽略的坑面积重叠法的调参我一直按先定阈值、再定重叠率、最后处理异常的顺序来。阈值选择影响最大。-52℃适合统计MCS但业务预警时我经常同时跑-32℃和-62℃两套。为什么-32℃能更早发现正在发展的云团-62℃能锁定冷云核的强中心算出来的是强度骨架。三套阈值同时看等于既看云团外框又看核心。重叠率阈值需要根据数据时次间隔动态调整。如果数据是5分钟间隔云团移动距离有限IoU取0.4以上很稳如果是15分钟间隔云团可能已经移动了一两个自身半径IoU取0.2都可能漏匹配这时候我会加一个最短中心距约束来辅助匹配或者用光流法先估算运动矢量再预测目标位置用预测位置附近一定半径内的目标做匹配。数据时间分辨率对追踪的连贯性影响极大。我做过对比实验同一场对流过程用15分钟间隔数据追踪生命周期统计结果和5分钟间隔的差别可以到30%以上。所以如果要做科研统计建议坚持用同一套观测模式的数据不要混着用全圆盘和区域加密模式的数据。5. 一次华北冷涡背景下对流过程的完整追踪复盘5.1 天气背景为什么选这场过程拿去年夏末一次典型的华北冷涡过程来复盘。500百帕图上冷涡中心在蒙古国东部到内蒙古中部之间华北正好处于冷涡后部偏北气流控制午后地面气温快速升高。700百帕以下有明显的切变线低层水汽条件不算特别充沛但足够维持强对流。午后CAPE值对流有效位能达到1500到2200焦耳每千克垂直风切变属于中等偏强。这种配置在华北是雷暴大风的常见配方而且对流很容易在太行山东麓就地激发。5.2 从卫星数据里读出的完整时间线我把当天下午风云四号的区域加密观测数据调出来按时间顺序把关键节点捋了一遍。12时30分左右太行山东麓出现零星浅对流10.8微米TBB大约在-20℃上下云顶很低如果只看这一景云图很容易忽略。但把前一个时次拉过来对比15分钟内降温达到2.5K趋势在变冷。13时10分有一个云块的TBB降到-38℃左右15分钟降温达到4.2KCI判据触发。同时水汽图像上能看到从西北方向切入的干区边缘正好压在这个云块附近高空辐散条件在转好我判断它会在接下来一小时内组织化。13时50分左右云顶TBB跌破-52℃冷云罩面积逐渐向东北方向扩展MCS的雏形出现。之后的发展速度比预期快14时30分冷云罩面积已经达到12万平方公里以上中心最低TBB跌破-65℃云顶形态上出现了明显的过冲迹象可见光云图上能看到云顶有密集的凸起。15时20分左右整个MCS分裂成两个中心。西侧那个冷云顶开始变暖、边缘变光滑逐渐进入消散阶段东侧那个反而加强强中心继续维持在-65℃以下。雷达资料显示东侧中心正好对应一条弓形回波带地面果然后来测到了8级以上的雷暴大风。16时40分东侧中心继续向东南方向推进地面自动站陆续出现短时强降水和雷暴大风。17时30分以后云顶TBB明显回升冷云罩破碎对流宣告减弱。整个过程持续大约5个小时卫星时间线上看得一清二楚。5.3 卫星和雷达配合起来的价值才最大单独拿卫星数据复盘能看出总体的起承转合但落区判断还是需要雷达。卫星的优势在于宏观趋势和提前量——它比雷达更早发现山前新对流启动的苗头因为雷达在山区有遮挡而且离雷达站比较远的区域要到一定回波强度才能被看到。雷达的优势在于低层的精细结构——中气旋、弓形回波、阵风锋这些东西卫星云图上看不到但它们是判断地面灾害天气的直接依据。我在这次过程中做预报的节奏是卫星负责提前发现锁定注意区域雷达负责确认结构圈定预警落区地面站和闪电定位负责验证强度判断灾种。三层资料缺一层都不够稳。6. 追踪对流几年踩过的数据坑从投影错位到时空不连续6.1 投影变换让云团悄悄漂移的坑处理风云四号数据第一个容易踩的就是投影问题。AGRI L1数据给出的往往是不规则扫描带上带有地理定位的产品要把它画到地图上或者和雷达数据叠加就得做投影变换。有人为了省事直接用等经纬度投影结果在中高纬度地区发生明显形变——云团边界被拉长中心位置偏移了几十公里。我现在的习惯是坚持用标称投影数据做分析也就是卫星原生扫描几何的那个坐标系避免反复重投影。如果需要叠加雷达或地面站数据先统一到同一个地图投影并且记录投影参数不要在分析过程中来回切换。看起来是小事但真实业务里因为投影不一致导致对流中心定位偏差、预警区域画偏的例子并不少见。6.2 质量标记和坏线别把传感器噪声当成对流增强红外通道的L1数据里多少会存在坏线、条带或者杂散光尤其是太阳耀斑影响的时候可见光通道更是重灾区。早期我做CI检测时曾经把太阳耀斑在红外通道产生的异常亮温当成了快速降温信号误报了一次对流初生。后来学乖了处理数据之前一定先看质量标记字段把无效像元、坏线像元直接掩膜掉再参与一切计算。如果某个区域的快速降温信号出现在数据边缘或者已知的坏线位置宁可不判断也不要强行给结论。样本上宁愿少一点准确率反而更高。6.3 冷云罩大不等于地面灾情重这是最核心的经验之一。MCS的冷云罩有时能扩展到几十万平方公里看起来气势汹汹但地面真正出现灾害性天气的区域往往只有冷云罩的十分之一甚至更小。冷云罩里面包含了一大片高空的卷云砧这些卷云是被高空风从主对流区吹出去的底下可能一点降水都没有。所以看MCS追踪结果时我从来不用冷云罩整体去圈预警范围而是把注意力放在几个关键指标上冷云罩面积增长率越快说明系统正在组织化最低TBB走低的趋势越陡强中心越活跃冷云罩长轴方向的移动路径才是下游需要重点关注的区域。真正常用的预警落区还是得回到TBB低值中心和雷达回波强中心的重叠位置。6.4 加密观测和常规观测切换造成的时间序列断裂风云四号在工作时会根据任务需要切换全圆盘、中国区域、局部加密等不同观测模式。这就导致一个比较别扭的问题——你拿到的数据时间间隔不是恒定的可能前10帧是5分钟间隔中间突然变成15分钟间隔切回加密模式后又变成1分钟间隔。对于面积重叠法这种依赖均匀时间采样的追踪算法时间间隔忽长忽短会让速度计算、面积增长率计算都出现刺跳。处理办法是先把整个时间序列按目标时间间隔重采样缺的时次用前后时次做线性插值或者光流插值。但插值出来的数据毕竟不是真实观测不能当作实况来用只能用于算法中间过程。如果做科研统计我建议干脆只取模式一致的加密观测时段宁可少一段也好过混着不同采样频率的数据硬算。另外还有一个容易被忽略的点同一颗星在不同时段可能会因为轨道调整或者冷光学部件温度变化导致红外通道的绝对亮温有少量漂移。做长期趋势分析时尽量用官方定标产品不要直接用原始的L1计数值转换。追了几年对流我的体会是卫星给的是大势雷达给的是明细地面站网给的是实况三样东西必须合在一起看。尤其这几年AI外推算法越来越多大家都想让机器从卫星云图序列里直接预测未来一小时的对流发展我也试过不少框架。说实话效果在变好但遇到分裂合并、冷云罩与地面灾情错位这些特殊情况时经验判断依然是最后的底线。把这套方法吃透再看那些算法输出的时候心里就有谱了。