关于使用二维视频雨滴谱仪作为参考仪器对Parsivel雨滴谱仪测量的雨滴谱进行校正的研究报告
本研究由瑞士洛桑联邦理工学院（EPFL）环境遥感实验室的T. H. Raupach和A. Berne完成，并于2015年发表在学术期刊《Atmospheric Measurement Techniques》（atmos. meas. tech.）上。
一、 学术背景
本研究的核心科学领域是大气物理学和降水微物理学的观测技术。雨滴谱（Drop Size Distribution, DSD）是描述降水微观结构的关键参数，它量化了单位体积空气中不同大小雨滴的数量分布。精确的DSD测量对于理解降水过程、改进天气雷达定量降水估计（QPE）、研究土壤侵蚀、植被截留以及电磁波在大气中的传播等至关重要。
Parsivel（Particle Size and Velocity，粒子尺寸和速度测量仪）是一种广泛应用的光学激光雨滴谱仪。它通过测量雨滴穿过一束激光时造成的阴影持续时间来估算雨滴的粒径和下落速度，并将雨滴归类到预设的粒径和速度区间（bin）中。然而，已有研究表明，Parsivel（包括第一代和第二代Parsivel2）在测量雨滴浓度，特别是对小滴和大滴的测量上，存在显著误差。例如，它可能高估小滴（如直径0.2-0.4 mm）的数量，并错误估计大滴的下落速度。这些误差影响了其DSD测量的准确性，进而影响基于DSD导出的所有降水参数（如雨强、液态水含量、雷达反射率因子等）的可靠性。
相比之下，二维视频雨滴谱仪（Two-Dimensional Video Disdrometer, 2DVD）被广泛认为是地面DSD测量的高精度参考仪器。它使用两套垂直的高速线扫描相机系统，能够直接、立体地测量每个雨滴的形状、大小和速度，无需依赖理论下落速度模型，因此其测量结果更为可靠。
基于此背景，本研究旨在解决一个具体的应用问题：如何利用高精度的2DVD作为参考，开发和验证一种通用的校正方法，以显著提高Parsivel雨滴谱仪（包括两代产品）所测DSD数据的准确性。研究的最终目标是使得Parsivel这种成本较低、坚固耐用、适合组网观测的仪器，在经过校正后，其DSD数据质量能够接近参考仪器水平，从而更可靠地用于研究DSD的时空变异性和相关应用。
二、 详细研究流程
本研究是一个系统的、数据驱动的方法开发与验证过程，主要包括数据收集、数据处理、校正方法开发、方法应用与效果评估等步骤，并在两个不同气候区进行了测试。
1. 数据收集与研究对象 研究使用了来自两个野外实验活动的数据： * 主要训练与验证数据：来自法国阿尔代什省（Ardèche）的地中海气候区HYMEX（Hydrological Cycle in the Mediterranean Experiment）项目2012年（SOP2012）和2013年（SOP2013）秋季的观测。实验网络包括多台第一代Parsivel、一台2DVD以及与之并置的翻斗式雨量计。在2013年的实验中，还在2DVD站点并置了一台Parsivel和一台Parsivel2，这为方法开发提供了理想的训练数据。 * 独立验证数据：来自瑞士帕耶尔纳（Payerne）的温带气候区2014年春季的观测。该站点包含一台第一代Parsivel与一台2DVD并置。此数据集用于测试所开发校正方法在不同气候条件下的普适性。
2. 数据处理与质量控制 在对原始数据进行校正前，首先对Parsivel和2DVD的原始观测进行了标准化处理和严格的质量控制。 * 数据转换：将Parsivel记录的原始粒子计数，以及2DVD记录的单个雨滴信息，统一转换为按粒径区间（bin）统计的雨滴数浓度 N(D)（单位：个/立方米/毫米）。 * 可疑粒子过滤：基于物理合理性，设定过滤器去除不可能为真实雨滴的观测记录（如昆虫、蛛网水滴等）。过滤条件包括：去除直径大于7.5 mm的粒子；去除下落速度超过Beard (1976)理论终端速度模型预测值+4 m/s或低于预测值-3 m/s的粒子。此步骤基于2DVD的观测特征设定，并同时应用于两种仪器的数据。
3. 校正方法的开发（核心流程） 校正方法分为两个核心步骤，依次处理Parsivel在速度测量和浓度测量上的系统偏差。
步骤一：速度校正
i，计算该区间内所有粒子的平均测量速度。然后，将所有粒子的速度值进行整体平移（shift），使得该区间的平均速度等于该区间中心直径Di所对应的Beard (1976)理论终端速度v(Di)。速度校正完成后，再应用上述“可疑粒子过滤”步骤。此步骤不改变各粒径区间的总粒子数，仅调整其速度分布。步骤二：浓度校正（核心创新）
i，定义校正因子 P(i) = N_2DVD(i) / N_Parsivel(i)，即在每个时间步长（本研究使用1小时分辨率以增加大滴样本），2DVD与Parsivel测得的该区间雨滴数浓度之比。P(i) 的值与降水强度有关。因此，研究根据Parsivel仪器自身提供的雨强值（Parsivel-derived intensity，独立于DSD计算），将训练数据划分为不同的雨强等级：[0, 0.5), [0.5, 1), [1, 2), [2, 200) mm/h。i，计算所有有效时间步长P(i)的中位数。这个中位数即为该雨强等级下、该粒径区间的最终校正因子。P(i)的统计波动较大。研究通过多次随机抽样（bootstrap）计算P(i)中位数的变化范围，并设定一个阈值（变化范围<0.7），只对满足该条件的粒径区间应用浓度校正，以确保校正的稳健性。对于第一代Parsivel，校正最高应用于直径约5.5 mm的雨滴。P(i)值作为乘数，乘以Parsivel原始测量的N_Parsivel(i)，即得到校正后的雨滴数浓度。4. 方法验证与效果评估 将训练得到的校正因子（基于HYMEX 2013数据）应用于： * 同一数据集（HYMEX）的高分辨率（5分钟）数据：评估校正对DSD各阶矩（从0阶矩总浓度到高阶矩）的影响，并与2DVD的DSD矩进行对比。 * 独立数据集： * HYMEX 2012数据：作为时间上独立的验证。 * Parsivel2数据：使用相同流程为第二代Parsivel2训练了一套独立的校正因子，并评估其效果。 * 瑞士Payerne数据：作为不同气候区的独立验证，使用在法国训练得到的校正因子直接应用，测试方法的通用性。 * 最终验证指标：使用并置的翻斗式雨量计的雨强测量作为独立标准，评估校正前后Parsivel计算出的雨强（由校正后的DSD根据公式积分得到）的准确性。评估指标包括均方根误差（RMSE）、偏差（Bias）、相对偏差（Relative Bias）和决定系数（R²）等。
三、 主要结果
DSD矩的改善：将校正方法应用于HYMEX 2013并置站点的5分钟分辨率数据后，与2DVD相比，Parsivel所测DSD的各阶矩（0至7阶）的分布都显著向2DVD的分布靠拢。具体而言，低阶矩（0、1、2阶，对应总浓度、平均直径等）的相对偏差大幅改善（例如0阶矩改善了约124%），均方根误差（RMSE）全面降低，相关性（R²）对低阶矩也有提升。这表明校正有效地将Parsivel的DSD整体“校准”到了更接近参考仪器的水平。
雨强估算的改善：与并置雨量计的对比表明，在大多数情况下，校正改善了Parsivel的雨强估算性能。
对Parsivel2的适用性：研究成功地将相同的方法论应用于第二代Parsivel2。虽然Parsivel2的校正因子曲线与第一代表现出不同的形态（例如，在1.38-3.25 mm直径范围内，Parsivel2表现为低估，而第一代是高估或接近），这反映了两代仪器硬件（如激光）的差异，但校正方法本身是有效的。应用校正后，Parsivel2的DSD矩（尤其是低阶矩）和雨强估算也得到了改善。
方法的普适性：将基于法国地中海气候数据训练得到的第一代Parsivel校正因子，直接应用于瑞士温带气候的Payerne站点数据，取得了良好的效果。校正后，Parsivel与2DVD测得的DSD矩的一致性显著提高，与雨量计对比的雨强估算误差也明显减小。这强有力地证明了该校正方法对不同气候条件具有鲁棒性和普遍适用性。
四、 结论
本研究开发并验证了一套系统的、基于参考仪器（2DVD）的、分两步走的校正流程，用于显著提高Ott Parsivel雨滴谱仪（包括第一代和第二代）测量的雨滴谱（DSD）数据的准确性。该方法首先对雨滴下落速度进行物理模型约束下的校正和过滤，然后利用统计方法对各个粒径区间的雨滴数浓度进行基于雨强分档的标定。
本研究的价值在于： * 科学价值：提供了一种可靠的方法来提升广泛使用的Parsivel仪器的数据质量，使得利用其组网观测研究DSD空间微小尺度变异性成为可能，为降水微物理学、雷达气象学等领域提供了更可靠的地面验证数据。 * 应用价值：该方法不依赖于特定气候条件，具有普适性，可供全球使用Parsivel的研究人员和业务单位采用。它操作性强，仅需Parsivel自身的输出数据和一套预先训练好的、按雨强分档的校正因子表即可实施。研究还明确指出了校正方法在极低雨强下的局限性，并给出了使用建议。 * 重要观点：研究表明，通过精心的、基于物理和数据驱动的校正，低成本、坚固的仪器（如Parsivel）在经过严格标定后，其核心数据产品（DSD）的准确性可以接近高精度科研级仪器（如2DVD）的水平，这极大地拓展了此类仪器在科学研究中的应用潜力。
五、 研究亮点