基于化合物特异性角质层渗透性的Level IV逸度模型解析电离依赖的杀菌剂向葡萄浆果渗透与残留分布

《Ecological Modelling》：Modelling ionisation-dependent fungicide penetration into grape berries using a Level IV fugacity model with compound-specific cuticular permeability

　　

原文未设独立摘要段落，仅有研究亮点；其专业翻译如下：含电离作用的葡萄浆果杀菌剂残留Level IV逸度模型（Level IV fugacity model）。通过Henderson–Hasselbalch（H-H）统一框架处理极性与非极性化合物。DT₅₀较亲脂性更能预测果皮残留分数（ρ=0.96）。平均温度和DT₅₀主导模型不确定性（|PRCC|=0.88–0.95；PRCC为偏秩相关系数，Partial Rank Correlation Coefficient）。区域温度变化使持久性杀菌剂的收获前间隔期（Pre-Harvest Interval，PHI）偏移>10 d。

研究背景与问题：葡萄（Vitis vinifera L.）是全球重要经济果树，鲜食葡萄常整果带皮消费且采后至摄入前加工极少，因此采收时浆果中的杀菌剂残留可能被消费者全部摄入；真菌病害是葡萄生产的主要生物限制，重复施药常见，葡萄在欧盟和巴西监测中常出现最大残留限量（Maximum Residue Limit，MRL）超标。消费者对低残留葡萄需求增加，认证体系设置低于法定MRL的限值，因此定量理解杀菌剂在浆果内的理化归趋对优化药剂选择、施药时期和剂量、兼顾作物保护与消费者暴露最小化十分关键。逸度（fugacity，f，Pa）是表征化学物质从相中逃逸趋势的热力学势；Level IV逸度模型可给出时间分辨残留曲线并纳入温度、降解和作物生长。Paraíba et al. 的周期Level IV控制系统具有状态矩阵A(T)非负性和周期解全局稳定性，但原示例未包含极性杀菌剂在酸性果肉pH下的电离，也未包含葡萄角质层对分子运输的化合物特异性阻力。Riederer and Schreiber 建立了角质层渗透系数与辛醇-水分配系数和分子量的经验关系。该研究旨在将周期Level IV框架应用于葡萄浆果，定义边界空气层、脂质角质层和水性果肉三个生物隔室，通过H-H校正电离、通过Riederer–Schreiber模型引入化合物特异性角质层渗透性，在单一统一模型中处理可电离极性和非极性化合物；检验降解动力学而非热力学分配决定最终残留分布的假设，即DT₅₀应比f_n·K_ow或P_cut更能预测果皮残留分数；并模拟11种在巴西葡萄上登记的杀菌剂，结合食品安全和MRL合规讨论。

关键技术方法（≤250字）：模型将浆果设为半径8 mm球体，含边界空气层、脂质皮肤（角质层和表皮蜡）和水性果肉三个同心隔室，界面面积A_int=8.04×10^?4 m²；施加24 h正弦温度循环（均值25 °C，振幅10 °C，范围15–35 °C），果肉pH 3.5，采用矩阵指数以1 h步长离散，模拟240 h，按近采收Stage III保守情景且不模拟雨水洗脱；用H-H校正电离，用Riederer–Schreiber模型计算P_cut；降解速率由实验DT₅₀经Arrhenius方程温度校正，并采用化合物特异性E_a、Q₁₀和ΔH_vap；参数来自PPDB、EFSA同行评审和监督残留研究及文献；对11种杀菌剂进行Monte Carlo（MC）不确定性分析（每种500组、共5500次运行），以PRCC排序参数重要性；结果用24 h中心移动平均后处理。无田间样本队列，实验比较数据来自文献报道的果皮/果肉分布数据。

研究结果：3.1 角质层渗透性、分配与最终残留分布：通过汇总Riederer–Schreiber渗透性、有效分配系数和240 h最终残留分布，发现P_cut跨40倍范围、P_cut/K_ow跨120倍范围；采用实验DT₅₀后降解主导质量平衡，52%–100%初始质量在240 h内损失，残余质量积累于脂质皮肤，DT_50,pulp与果皮残留分数的Spearman相关为ρ=0.96（p<0.001）。3.2 图形结果：通过24 h移动平均渗透曲线，短暂峰值后240 h所有化合物果肉渗透<0.4%，短DT₅₀化合物在1–2天达峰后下降；质量平衡确认降解主导；模型定性重现实验趋势，但定量高估果皮残留分数F_res^skin（RMSE=22.8百分点；Bias=+20.6百分点；N=6），iprodione和tebuconazole偏差约30百分点；作者归因于缺表面隔室、雨水洗脱、品种特异性角质层及实验皮肤/果肉分离不完全。3.3 与其他建模方法比较：通过比较已发表模型，本模型唯一同时含周期温度依赖、H-H电离校正和化合物特异性角质层渗透性，并区分空气、皮肤和果肉；DynamiCROP聚合皮肤与果肉且无电离校正，Trapp模型面向根到冠吸收而非叶面沉积和角质层渗透，Fantke and Jolliet模型亦聚合皮肤与果肉。3.4 不确定性分析与蒙特卡洛发现：通过MC和PRCC发现平均温度与DT₅₀主导皮肤残留分数和降解分数不确定性（|PRCC|=0.88–0.95），F_cal为次要参数（PRCC=0.10–0.28）；持久性化合物可信区间宽（皮肤分数近0%至约70%，P₉₅/P₅>400），短寿命化合物区间窄；区域温度差异使持久性化合物皮肤分数变化29–33百分点，PHI₉₀差>10天；参数重要性层级为DT₅₀≈平均温度>E_a,pulp>F_cal>log K_ow≈pH。3.5 对食品安全与收获前间隔期的启示：通过模型结果，果皮积累的杀菌剂可经家庭洗涤部分去除，但不洗消费者摄入全量；cymoxanil短暂果肉渗透（约1.3%，2天内）不可洗除；cymoxanil、chlorothalonil和mancozeb的空气相分配最大值分别为27.1%、4.5%和4.2%，可降低果实残留负担。

讨论部分总结：模型与实验比较的系统正偏差主要源于Level IV框架缺少表面隔室，而非温度；雨洗排除、品种角质层差异和实验分离不完全亦参与。F_cal=1000与Górna-Binkul et al. 的蜡层扩散数据一致（隐含F_cal=612–1245）。不确定性分析表明应优先获取葡萄果皮和果肉在不同温度下的DT₅₀及V. vinifera原位角质层渗透性数据。模型仅代表角质层沉积途径，未考虑土壤–根吸收、韧皮部运输和滴灌等非叶面途径；未来可加入土壤–根吸收模块、韧皮部H-H分配和多途径整合。

研究结论翻译：该研究将Paraíba et al. 的周期Level IV逸度控制系统应用于葡萄浆果，定义三个生物隔室，并扩展为统一处理可电离极性和非极性化合物的模型；状态矩阵A(T)的非负性在11种杀菌剂中不论电离状态均保持。f_n·K_ow是决定皮肤/果肉偏好的热力学参数，11种中10种f_n·K_ow?1，仅cymoxanil（f_n·K_ow=0.15）具有向果肉的驱动力；Riederer–Schreiber模型引入40倍P_cut范围和120倍P_cut/K_ow范围。采用实验DT₅₀后，降解成为主导过程，52%–100%初始质量在240 h内降解，DT_50,pulp与残余果皮分数的Spearman相关为ρ=0.96，表明持久性比亲脂性更能决定长期果皮残留。行业上应优先选择DT_50,pulp<5天的药剂；持久性化合物（pyraclostrobin、metrafenone、difenoconazole等）在皮肤中保留较高比例，消费者暴露风险最高。监管上，巴西葡萄产区10 °C温度差异可使持久性化合物预测果皮分数变化29–33百分点；在25 °C暖参考温度制定的PHI对较冷产区（如Serra Gaúcha，18 °C）系统性非保守，应延长，而对暖区（如S?o Francisco Valley，28 °C）可缩短。cymoxanil说明热力学与动力学交互：其向果肉的热力学驱动力真实存在，但DT₅₀=1.5天使其在240 h内无可测果肉渗透。用化合物和隔室特异性实验DT₅₀替代统一μ_ref=0.003 day^?1，修正了11种杀菌剂降解低估15–154倍的问题；量纲分析确认状态变量f*_k与逸度容量Z_k一致；温度驱动的24 h振荡经中心移动平均去除而不影响累积趋势。该逸度框架提供了机制透明、参数可追溯的工具，将分子理化性质直接联系到消费者暴露指标，支持PHI制定和残留安全葡萄生产的杀菌剂选择。