《Methods in Ecology and Evolution》：Ground-penetrating radar amplitudes to roots: A numerical modelling framework for mapping coarse root system architecture
编辑推荐：
三维（3D）粗根系统架构的观测仍然有限，因为大多数测量方法需要开挖。在此，研究人员将密集的单通道和多通道500 MHz探地雷达（ground-penetrating radar, GPR）调查与流动引导连通性框架（flow-guided connectivit
三维（3D）粗根系统架构的观测仍然有限，因为大多数测量方法需要开挖。在此，研究人员将密集的单通道和多通道500 MHz探地雷达（ground-penetrating radar, GPR）调查与流动引导连通性框架（flow-guided connectivity framework）相结合，以非破坏性方式表征了两棵树的粗根架构。偏移归位（migrated）和深度归一化（depth-normalised）的GPR振幅在稳态数值模型（steady-state numerical model）中被视为空间变化的水力传导率（hydraulic conductivity）。所得连通性场（connectivity fields）将空间不连续的、塌缩的衍射双曲线（diffraction hyperbolae）转化为连续的3D根系架构表征，从树干向外延伸10–15 m。研究人员将该方法应用于沼泽栎（*Quercus palustris*）和河流红桉（*Eucalyptus camaldulensis*）。这两棵树在不同环境条件和采集参数下被调查。重建的架构揭示了浅层、侧向广泛的粗根系统，并捕捉到两棵树之间根系结构的显著变化。栎树根系表现出不对称架构，优先向树上坡、风暴露侧（upslope, wind-exposed side）发育，而河流红桉则显示更呈径向分布（radially distributed）的根系架构，近树干处有更深根系。两棵树约80%的重建根系位于距树干6.2 m以内且60 cm深度以上。总之，这些结果表明，密集GPR调查与基于连通性的建模相结合，为表征粗根架构提供了一种非破坏性和可重复的方法。
论文解读：
研究背景与问题：三维（three-dimensional, 3D）粗根系统架构对理解树木获取水分和养分、机械锚固、以及影响土壤—大气间碳水平衡至关重要，但原位观测仍然困难。传统基于开挖、 trenching 或 coring 的方法具有破坏性且劳动密集；虽然可重建3D根系架构，但要达到足够空间分辨率通常需要大面积开挖，并损伤或杀死树木。非破坏性方法如液流测量、同位素示踪和土壤水分监测只能间接约束根深和径向范围；根生物量及根系范围估计常依赖由开挖数据校准的地上异速关系。地球物理方法提供了非侵入替代方案，其中电阻率层析成像（electrical resistivity tomography, ERT）可推断根区和根—土相互作用，但空间分辨率通常不足以详细重建根系架构；高分辨率地震和声学方法已在特定场景探索；声波和声纳断层成像多用于树干内部结构和腐朽表征。探地雷达（ground-penetrating radar, GPR）可非侵入检测大直径根，具有高空间分辨率、快速采集和低场地扰动优点，大直径根在GPR剖面中常产生特征衍射双曲线。然而，根可检测性强烈依赖土壤质地与湿度、天线频率、根直径以及根相对天线方向；尽管GPR可检测单根，但从密集GPR调查重建稳健3D根系架构仍具挑战。近期研究关注3D成像和自动/半自动重建，但密集采样时，偏移归位后衍射双曲线塌缩产生的高振幅异常仅出现在测线几何与单根相交处，导致重建根系系统不连续，妨碍体积分析和跨物种、地点、环境设置的定量比较；人工解释和手工拾取数百至数千剖面费时且难以常规应用。因此，研究人员将问题重构为异质介质中的稳态流：高振幅异常被视为嵌入土壤基质中的导管，树干作为垂直源，理论流体被推入土壤，优先沿与树干连接更强的路径流动，从而将重建问题类比为裂隙岩含水层中的稳态地下水流，并利用粗根通常从树干起源、径向向外延伸并形成连通网络这一特征。
研究人员开展了密集GPR调查与连通性建模结合的研究，并应用于沼泽栎（*Quercus palustris*）和河流红桉（*Eucalyptus camaldulensis*）。两棵树处于对比气候、土壤和半球，并使用不同采集策略。研究人员用偏移归位、深度归一化的GPR振幅填充稳态数值模型，在无开挖条件下生成3D根系架构。结果显示，该方法能将空间不连续的GPR异常转化为连续的树干连接结构，重建出浅层、侧向广泛的粗根系统，并捕捉两棵树间显著结构差异：栎树根系不对称，优先向树上坡、风暴露侧发育；河流红桉根系更径向分布，近树干处有更深根。两棵树约80%重建根位于距树干6.2 m以内且60 cm深度以上。该框架为研究物种、环境条件和扰动历史如何影响根系发育提供了机会。
关键技术方法（不超过250字）：采用密集500 MHz GPR采集。沼泽栎位于美国东南部，以单通道角向网格调查，15 cm剖面间距、1 cm道距，共1170条剖面、3.6 km；河流红桉位于南澳大利亚Cambrai，以MAL? MIRA多通道系统调查，54条swath、每swath 22条共偏移剖面、6.5 cm间距、6.7 mm触发，共40.8 km。处理包括dewow、均值道滤波、200–1000 MHz带通、时变增益；由衍射双曲线拟合估计电磁（electromagnetic, EM）速度，并进行FK偏移；计算包络、深度归一化并平方。平方归一化包络作为输入，用PyGIMLi在非结构网格上有限元求解稳态流方程，得到连通性。另用基于已发表橡树根系数字化图像的二维合成模型和GPRMax正演评估。
研究结果：
5.1 Electromagnetic velocities（电磁速度）。通过共偏移GPR剖面中拟合衍射双曲线估计EM速度。沼泽栎识别1001条双曲线，速度分布为高斯型，均值0.105 m ns
?1、标准差0.013 m ns
?1，87%在0.08–0.13 m ns
?1，76%位于60 cm以上，96%浅于1.0 m，无系统深度趋势，故用0.105 m ns
?1偏移。河流红桉识别927条双曲线，均值0.109 m ns
?1、标准差0.013 m ns
?1，94%在0.08–0.13 m ns
?1，71%位于60 cm以上，97%在1.0 m内，无深度趋势，故用0.109 m ns
?1偏移。
5.2 Two-dimensional root architecture（二维根系架构）。在20、40、60 cm深度提取振幅和连通性切片。河流红桉振幅切片显示侧向广泛、狭窄、径向取向的异常网络，20 cm主要向西北，并有东南、西南异常，呈树枝状分支；40 cm主浅层特征持续但振幅连续性降低；60 cm径向异常仍可检测但更近树干终止。连通性切片保留振幅空间模式，同时强调异常与树干间连续性；20 cm多条高连通路径径向延伸，深度增加后连通性集中于较少路径，边缘特征减弱，路径比振幅场更连续。沼泽栎20 cm高振幅异常集中于活动建模域边界；40 cm近树干复杂度增加且沿树干预始路径连续性增强；60 cm由三个主要异常向正方向延伸。连通性切片将离散振幅异常连回树干；浅层有未连树干的短路径，较深处出现较少但更连续的树干连接结构，包括40 cm向东北边界的显著路径。总体，重建架构保留观测数据主导结构，填补缺失或噪声振幅间隙，偏好径向连续性。
5.3 Three-dimensional root architecture（三维根系架构）。连通性给出两棵树连续3D根系架构。河流红桉形成复杂树枝状网络，有弯曲和分支结构；多个主结构合并成从树干延伸超过12 m的连续路径，三条主结构向树冠边缘辐射形成近三角形；阈值轮廓差异大，反映广泛次级分支。沼泽栎呈更直、更连贯的径向连接路径，分支有限；阈值差异较小，根系表示更锐利；平面图中主根集中于右上和右下象限，对应上坡侧。根体积：沼泽栎20–43 m
3，河流红桉18–111 m
3。降维后，两棵树多数根集中在约70 cm以浅；沼泽栎主根系统集中于30–70 cm，100 cm以下少，1.5–6.5 m径向距离间约40 cm以浅存在缺根带，6–6.5 m处分叉，一支上至地表，一支下至近100 cm；河流红桉近地表根集中并随深度递减，呈宽锥形，底部约100 cm、半径约7.5 m，根分数超过0.25出现在上部50 cm并侧向延伸约7.5 m，可检测根分数近地表超过15 m，近树干100 cm以下有二次簇并侧向约2.5 m，7.5 m附近也有分叉。垂直根分数显示沼泽栎有四至五个主结构，右上树枝状结构较浅，偏向正方向；河流红桉有三个主结构近三角排列，并强调从树干到活动建模域边缘约17 m的侧向连续树枝状结构。
讨论部分总结：
6.1 Are these really roots? 单条GPR剖面中，衍射双曲线不能唯一归因于根，也可能由corestones、土壤界面突变或其他电磁性质强反差物体产生。但两站点在数百至数千独立剖面中，小高振幅异常经偏移后呈bulls-eye特征，并对齐成近树干起源、延伸至树冠边缘甚至更远的复杂树枝状结构；这种树干中心连通网络不太可能偶然形成。连通性不引入新结构或规定根几何，而是以径向连续性偏好填充异常间空间。直接独立验证有限：河流红桉无开挖数据；沼泽栎曾用air-spade在基于高振幅GPR响应位置暴露粗根，但验证基于振幅而非连通性架构。根归因主要基于与树干的连续性，重叠根区归属仍有不确定性。
6.2 What does connectivity do? 连通性是一种数值变换，把空间不连续GPR振幅异常插值为连续结构，并偏好径向排列。稳态模型中树干为源、外边界为汇，均匀条件下形成径向对称流场；高振幅区域作为相对水力传导率代理，形成连接树干与模型的路径，从而产生较大连通性。连通性反映路径连接树干与周围区域的效率，不是根存在或根直径的直接物理量。研究人员通过阈值化推断根架构，较低阈值提供包容表示，保留弱路径和次级分支；较高阈值分离最强、最连续主路径。未来需验证和校准阈值。结构随距树干增加而变细，并非模型显式施加，而源于根直径减小导致GPR振幅系统性降低以及小于有效分辨率物体偏移后塌缩至最小可分辨尺寸。
6.3 Does the root architecture agree with global root distributions? 将重建架构降维为概率密度函数（probability density function, PDF）和累积分布函数（cumulative distribution function, CDF），比较径向距离、深度和方位。两棵树虽3D组织差异明显，但约80%重建根出现于距树干6.2 m内。结构根区（structural root zone, SRZ）基于树干直径，覆盖约75%沼泽栎根和85%河流红桉根。河流红桉深度分布落入全球根系分布包络，并与澳大利亚Eucalyptus深度剖面接近，呈弱双峰；沼泽栎呈单峰，符合温带森林常见高斯型剖面。方位PDF显示两棵树粗根架构强各向异性：沼泽栎偏好正方向，河流红桉有三个主要方位模式。圆形保护带如SRZ不能反映粗根方位分布。
6.4 Why does this matter? 密集GPR与连通性建模可在无开挖条件下恢复成熟树下有连贯3D粗根系统。重建主要代表粗根骨架而非细根；细根主要负责水养吸收且常更深。河流红桉有较深二次根组件，与限水环境常发育更深根的概念模型一致，但本研究不直接解析根功能或吸水。工作流计算需求低，单通道沼泽栎采集约12–16 h、四人，数据约750 MB、1180个文件；多通道河流红桉约4 h、三人，但因7 mm道距产生约12.2 GB。信号处理和连通性建模可在标准笔记本电脑上逐线完成，主要投入为数据采集和仪器。两棵树用于评估不同环境下工作流性能，而非隔离物种、土壤、气候或树大小影响。
6.5 Limitations, resolution and detection limits。GPR根成像最小可检测根尺寸取决于土壤和根性质，包括天线频率、土壤电导率、根深、根方向、根间距和根—土电磁性质反差；低含水量根可能对GPR信号不可见。有效分辨率不同于传统反射界面定义，因不依赖分辨单根上下边界。对约0.11 m ns
?1速度和500 MHz中心频率，波长约22 cm，四分之一波长分辨率约5.5 cm；但根被解释为产生衍射双曲线的散射体，小于该标准的物体仍可能产生可检测响应。GPRMax合成正演显示，理想圆柱根直径10–1.25 cm、埋深30、60、80 cm时，小至1.25 cm直径仍可产生可检测衍射响应。GPR异常也不一定对应单根，根簇或裂隙中根可共同产生相干响应；但解释基于数百至数千异常的空间组织而非孤立响应。密集采样、偏移和减小剖面间距对改善重建重要。
7 CONCLUSION（结论）翻译：本研究证明，密集GPR采集与基于连通性的建模框架相结合，能够非侵入式重建粗根架构的连贯3D表征。通过将空间不连续的衍射响应转化为连续的树干连接结构，该方法解析了横向范围、深度分布和方位不对称模式，这些模式难以用传统根系采样方法表征。应用于沼泽栎和河流红桉表明，工作流从两个独立采集、在显著不同野外条件下收集的数据集中恢复了连贯根系架构。将3D重建降维为累积径向、深度和方位分布，为在没有开挖情况下解释GPR衍生根系架构提供了重要参考点，并将结果置于已建立的根系分布背景中。尽管连通性框架不量化细根或绝对根尺寸，它捕捉了将树木锚固于土壤中的粗根的组织和连续性。非破坏性重建这些架构的能力为表征成熟或保护树木下根系以及研究粗根架构随时间或响应环境变化如何发育或重组提供了新机会。
