双偏振气象雷达昆虫等效球体尺寸反演——基于差分传播相位(KDP)的Rayleigh近似方法

《Remote Sensing》：Estimating Body Lengths of Airborne Insects Using S-Band Polarimetric Weather Radar Kai Cui, Xinyu Li, Weidong Li, Jiayi Li, Xichao Dong, Rui Wang and Cheng Hu

　　

摘要：本研究提出一种基于S波段双偏振气象雷达(differential phase ΦDP及specific differential phase KDP)从昆虫回波中反演等效球体尺寸(a，等效半径)的新方法。假设飞行中昆虫呈随机水平取向且可用均匀介电椭球(Rayleigh近似)描述，研究人员推导了KDP与昆虫复极化率之差(α∥－α⊥，即Clausius-Mossotti因子β∥,⊥)及取向平均?cos2ω?f的关系式，并结合实测KDP值反算等效球半长轴a。利用典型昆虫组织相对介电常数εr≈(60＋j10)及退极化因子L∥,⊥进行数值验证。结果表明该方法可在不考虑昆虫质量或密度的前提下，由标准双偏振雷达观测量KDP获得昆虫等效投影面积或尺寸估计，为区分气象回波与非气象生物回波及昆虫群参数化提供定量依据。

论文解读——《Remote Sensing》发表：基于双偏振雷达K_DP反演昆虫等效球体尺寸的研究

一、研究背景与意义

传统天气雷达(WSR)在监测降水同时常受到鸟类、昆虫等非气象目标(biological scatterers / non-meteorological echoes)污染，其中空中昆虫(aerial insects / aeroplankton)在大尺度迁移(昆虫 migration / migrant insects)时产生显著后向散射(Z_DR、Φ_DP异常)，干扰定量降水估算(QPE)。现有方法多依赖反射率因子Z_{HDR，易受取向、形状及密度不确定影响。双偏振变量specific differential phase KDP(differential propagation phase per unit path length)具有不受衰减、雷达定标影响且对前向散射敏感的优点，适合用于估算昆虫等效尺寸。本文在Rayleigh散射近似下建立昆虫为均匀介质旋转椭球体，通过极化率取向平均导出KDP与昆虫等效球体积的关系，实现仅用雷达观测值反演等效球半径a，为生物气象学(aeroecology)及雷达生物回波识别提供物理基础。}

二、主要关键技术方法

研究人员采用理论推导结合数值验证：(1)将昆虫建模为长轴a、短轴b=c=a(等效球近似)的均匀各向同性介质旋转椭球体，介电常数取典型昆虫体液值ε_r≈60＋j10；(2)引入Clausius-Mossotti型极化率α_i=ε₀(ε_r－1)/[1+L_i(ε_r－1)]，其中L_∥、L_⊥为退极化因子(depolariation factor)；(3)假定昆虫主轴位于水平面且方位角ω服从均匀分布f(ω)=1/(2π)，对H极化和V极化有效极化率做取向平均?α_H?=α_⊥+(α_∥－α_⊥)?cos2ω?_f，其中?cos2ω?_f=∫?^2πcos2ω·f(ω)dω=1/2(随机取向)或指定分布；(4)利用波动方程复波数k=k₀√1+Nα/(2ε₀+…)一阶展开得到Δk=k_H－k_V≈(k₀N/2)(α_H－α_V)；(5)由K_DP=dΦ_DP/dr=Re(Δk)代入取向平均结果，反解等效球半径a=[K_DP/((Nπk₀/24)Re{β_∥－β_⊥}?cos2ω?_f)]^1/3；(6)数值计算不同a下的理论K_DP并与典型观测值比对验证一致性。样本为理论昆虫群，未用实地捕获标本。

三、研究结果（按原文小标题结构）

(1) Wave Propagation in a Medium Containing Polarizable Insects（含公式1–4）

研究人员从含稀疏可极化散射体的均匀介质出发，给出单散射体引起的复折射率修正n≈1+(Nα)/(2ε₀)，进而得到平行与垂直极化复波数k_H=k₀√ε_r+Nα_H/(2ε₀)及k_V=k₀√ε_r+Nα_V/(2ε₀)，经一阶Taylor展开得Δk=k_H－k_V≈(k₀N/2)(α_H－α_V)，建立了波数差与极化率差的直接联系。

(2) Differential Phase Shift and K_DP（含公式5–7）

定义specific differential phase K_DP=dΦ_DP/dr=Re(Δk)，代入Δk表达式得K_DP≈(k₀N/2)Re(α_H－α_V)，说明K_DP正比于昆虫数密度N与平行/垂直极化率实部之差，与雷达波长λ(隐含于k₀=2π/λ)相关而与雷达校准无关。

(3) Orientation Averaging of the Effective Polarizabilities（含公式8–10）

设昆虫长轴与入射H极化方向夹角为ω，H极化有效极化率为α_ωH=α_∥cos2ω+α_⊥sin2ω，V极化因主轴在水平面故α_V=α_⊥。对取向概率密度f(ω)平均得?α_H?=α_⊥+(α_∥－α_⊥)?cos2ω?_f，从而α_H－α_V=(α_∥－α_⊥)?cos2ω?_f。随机取向时?cos2ω?_f=1/2。

(4) Final Expression for Insect Size from K_DP（含公式11–12）

将椭球极化率代换为Clausius-Mossotti形式β_i=ε₀(ε_r－1)/[1+L_i(ε_r－1)]，并将昆虫体积等效为球体积V=4πa3/3，N为单位体积内总极化率贡献与单个昆虫极化率关系合并，最终导出K_DP≈(πk₀N/24)a3Re{(β_∥－β_⊥)?cos2ω?_f}，反解得等效球半径a=[K_DP/((πk₀N/24)Re(β_∥－β_⊥)?cos2ω?_f)]^1/3。若已知K_DP观测值、雷达频率(k₀)、昆虫介电特性(ε_r)及取向分布，可求a而无需知道N(若假设单只昆虫可分辨则N∝单位分辨率体积内昆虫数，实际可通过雷达体扫中K_DP剖面估算)。

四、讨论与结论翻译

讨论指出本方法依赖于Rayleigh近似(昆虫最大尺寸?λ，S波段λ≈10 cm满足大多数昆虫)、均匀介质假设及随机取向假设；对于强定向群(如顺风迁移)需用实测或假定取向分布修正?cos2ω?_f。该方法优点是不依赖雷达绝对标定、不受雨区衰减影响，且仅用常规双偏振变量K_DP即可获得尺寸信息，可与Z_DR、ρ_HV联合用于昆虫/鸟/气象分类。

结论（翻译）：研究人员建立了双偏振气象雷达specific differential phase K_DP与昆虫等效球半径a之间的解析关系，该关系基于Rayleigh散射近似、均匀介电椭球模型及取向平均处理。在假定空中昆虫飞行时呈随机水平取向的前提下，可由观测到的K_DP值反演昆虫的等效球体尺寸，所需辅助参数为昆虫组织的复相对介电常数和退极化因子。数值算例表明该反演公式在S波段对典型昆虫尺寸给出合理估计。此方法为利用业务双偏振雷达进行非气象生物回波识别及昆虫群微物理量估算提供了一种简便的物理途径。