起垄–开沟塑料覆盖系统垄沟位置温室气体通量水分响应的空间分异

《Agricultural Water Management》：Contrasting moisture responses of greenhouse-gas fluxes between ridge and furrow positions in a ridge–furrow plastic-mulching system

　　

起垄位置（ridge positions）在起垄–开沟塑料覆盖系统（ridge–furrow plastic-mulching systems）中较沟位置（furrow positions）温暖。起垄位置与沟位置在气体水分响应斜率（gas moisture-response slopes）上存在差异。CO₂和N₂O的水分响应对比（moisture-response contrasts）在敏感性分析（sensitivity analyses）中是稳健的。CH₄的水分响应对比在方向一致性（directionally consistent）方面表现稳健，但对野外数据排除（field exclusion）的稳健性较弱。单点采样（single-position sampling）可能偏离所假设的面积加权综合值（assumed area-weighted composites）。

起垄–开沟塑料覆盖系统（ridge–furrow plastic-mulching system）在湿润亚热带山地烟草种植中应用广泛，兼具防冲刷、抑杂草及维持根区温度等功能。然而，现有研究多将此类系统视为均质地表，解释温室气体（greenhouse-gas, GHG）CO₂、CH₄、N₂O通量时，往往忽略起垄与开沟位置因地表边界条件差异导致的水热空间分异。这种处理掩盖了系统内共存的异质微生境，导致不同研究对覆盖措施环境效应的结论不尽一致。因此，明确起垄与开沟位置是否具有不同的温室气体环境响应敏感性，成为完善温室气体监测与尺度转换的关键环节。

研究人员假设起垄与开沟位置代表不同微环境，因而具有不同的温室气体通量环境响应关系。研究以恩施土家族苗族自治州11块烟田为队列，通过重复观测比较两位置的水热条件及通量特征，解析温度与水分敏感性差异。结果表明，起垄位置较开沟位置温暖，但体积含水量差异并不稳健；三种温室气体的水分响应关系在位置间显著分异，且该分异较平均通量差异更具稳健性。这一发现揭示覆盖系统内温室气体响应的空间异质性，为将起垄与开沟作为独立观测单元进行监测、解释与尺度扩展提供了依据。

研究在湖北恩施、宣恩、建始三地的11块烟田设置固定观测点，每块烟田布设3个起垄位与3个开沟位，2022年6月至10月完成70次采样事件，获取420个点级观测。采用静态箱法（static chamber method）结合气相色谱（gas chromatography）测定温室气体通量；利用便携式土壤仪基于频域反射法（frequency domain reflectometry, FDR）原位测定土壤温度与体积含水量（volumetric soil water content, VWC）。统计分析运用线性混合效应模型（linear mixed-effects models）拟合通量–环境关系，提取位置特异性边际响应斜率，并通过最大似然估计、异方差稳健推断、MAD（median absolute deviation）离群值筛选、留一田块交叉验证（leave-one-field-out）及面积加权敏感性分析（area-weighting sensitivity analysis）验证结果稳健性，辅以配对差分模型（difference models）解析通量差异的环境驱动。

**起垄–开沟水热分异与平均通量差异**：事件级分层混合效应模型（event-level hierarchical mixed-effects models）分析显示，起垄位置土壤温度较开沟位置高0.76 °C（95% CI: 0.09–1.43; P = 0.029），而体积含水量差异未达统计显著（P = 0.342）。CO₂通量在起垄位置显著较高（71.44 mg CO₂ m^?2 h^?1; P < 0.001），CH₄平均差异不显著（P = 0.224），N₂O平均差异微弱且在对等田块加权下证据减弱（P ≈ 0.055）。

**起垄–开沟水分响应与温度响应斜率差异**：在控制土壤温度后，三种气体的水分响应斜率在位置间均存在显著差异。CO₂在起垄位置的斜率为46.31（95% CI: 26.27–66.36; P < 0.001），开沟位置为?13.45（95% CI: ?32.94–6.04; P = 0.176），斜率差59.76（P < 0.001）；CH₄在开沟位置斜率为0.0071（P < 0.001），起垄位置为0.0014（P = 0.460），斜率差?0.0057（P = 0.015）；N₂O在开沟位置斜率为?0.0090（P < 0.001），起垄位置为0.0042（P = 0.128），斜率差0.0132（P < 0.001）。温度响应斜率未呈现显著位置分异。

**位置响应模式的稳健性与配对环境关联**：通过替代随机效应结构、MAD离群值剔除、异方差设定、田块聚类整群抽样（field-cluster bootstrap）、留一田块剔除及面积加权情景分析发现，CO₂与N₂O的斜率差异方向稳定，CH₄方向一致但统计支持较弱；仅基于单点的观测会系统性偏离假设的面积加权综合值。配对差分模型显示，CO₂通量差异主要与温度差异关联，CH₄与水分差异关联，N₂O与温度–水分差异交互项关联。

讨论指出，起垄–开沟覆盖系统通过修改地表能量平衡、抑制蒸发与重分配降雨，形成非均质近地表环境，改变了温室气体交换所需底物运移、气体扩散、氧气供给与微生物活性的局地条件。研究强调，覆盖系统的环境效应不应被简单理解为整体通量的增减，而是不同气体响应关系的空间重分布。CO₂在起垄位置更强的正水分响应可能受根系自养呼吸（autotrophic respiration）与微尺度气体传输差异影响；CH₄在开沟位置更显著的正水分响应可能反映甲烷氧化（methane oxidation）与氧气扩散、水分动态的平衡差异；N₂O在开沟位置的负水分响应可能暗示更充分反硝化（denitrification）与N₂O还原。温度响应未呈现显著位置分异，表明在观测范围内水分响应分异更为清晰。研究同时指出，事件驱动采样设计（event-driven sampling design）适用于识别响应模式，但不适用于高精度累积排放估算；氧化还原电位（redox potential）、无机氮形态、功能基因等机制指标未直接测定，相关解释属合理推断。

研究结论表明，起垄位置平均温度高于开沟位置，而体积含水量差异因田块与事件而异；CO₂、CH₄、N₂O通量表现出位置依赖的水分响应关系，其中CO₂与N₂O的位置依赖响应对比最为稳健，CH₄方向一致但跨田块剔除时均匀性较弱。因此，在起垄–开沟塑料覆盖系统中，起垄与开沟位置不应被视为单一均质观测面，而应作为独立观测单元进行测量、解释与尺度扩展温室气体通量。