Определение архитектоники корневой системы деревьев методом георадиолокации

Размещено на http://www.allbest.ru/

Определение архитектоники корневой системы деревьев методом георадиолокации

П.А. Рязанцев¹, А.В. Кабонен², А.И. Родионов¹

¹ Карельский научный центр РАН, г. Петрозаводск, Россия

²Петрозаводский государственный университет, г. Петрозаводск, Россия

Методом георадиолокации проведено изучение корневой системы двух образцов яблони домашней (Malus domestica) и обобщены данные о методике наблюдений и обработки данных при изучении архитектоники корней. Показана возможность выделения погребенных объектов размерами от 1,5 до 5 см при использовании георадаров с зондирующей частотой 1-2 ГГц, что сопоставимо с ожидаемым размером скелетных корней. Выполнены полевые наблюдения георадаром ОКО-2 с антенным блоком частотой 1 700 Гц по исследованию архитектоники двух образцов яблонь, без применения способа раскопки, с целью определения влияния почвенно-гидрологических условий на процессы корнеобразования. По гиперболам дифрагированной волны на радарограммах определены местоположения скелетных корней исследуемых образцов, совокупность которых отражает распределение корневых систем в пространстве. При помощи трехмерных моделей проанализировано взаимоотношение морфологии кровли корненепроницаемых грунтов и корневой архитектуры образцов. Установлена трансформация корневой системы и наземной части деревьев в зависимости от положения подстилающих кристаллических пород. Обнаружено, что неглубокое залегание скальных грунтов приводит к изменению формы корневой системы от усечённого конуса к поверхностной.

Ключевые слова: Malus domestica; георадиолокация; корневая архитектоника; радарограммы; корненепроницаемый горизонт.

Pavel A. Ryazantsev¹, Aleksey V. Kabonen², Aleksandr I. Rodionov¹

¹Karelian Research Centre, Russian Academy of Sciences, Petrozavodsk, Russian Federation ²Petrozavodsk State University, Petrozavodsk, Russian Federation

Detection of the tree root system architecture using Ground penetrating radar

The principle of GPR is to examine the propagation of an electromagnetic wave in a medium by emitting a probing signal and receiving a reflected signal. Underground objects such as roots are distinguished by diffracted waves (See Fig. 1). The summary of literature data shows that roots up to 0.5 cm in diameter can be found by a high frequency antenna (See Table 1). We carried out a field experiment on an apple tree (Malus domestica Borkh) which grows on the territory of the Botanic Garden of Petrozavodsk State University (See Fig. 2). This territory is characterized by a thin soil cover and underlying moraine soils. Below are metatuff bedrocks. We studied two separate samples. The choice was made due to external differences (See Fig. 3). The trees growing on Plot 2 (Sample 2 of the experiment) have a conspicuous compound curvature of the trunk and a lower height than plants from Plot 1 - sample 1 (See Table 2). The two trees underwent annual average phenological stages simultaneously. They are 29 years old. No significant differences in the annual number of ripening fruits was revealed. Thus, the experiment was carried out on the two samples growing in equal conditions but having differences in appearance and biometrics. Literature data show that an apple tree usually has 4 perpendicular skeletal roots and 4 diagonal skeletal roots 1-2 cm in diameter in horizontal slice. As for vertical slice, tree roots usually lie at a depth of 20-60 cm (See Fig. 4).

The measurements were carried out with OKO-2 GPR by “Logis-Geotech” (Russia) with antenna unit (center frequency of 1700 MHz) on a regular rectangular grid with a space between profile lines of 50 cm. As a result, 26 profiles (13 horizontal and 13 vertical ones) were obtained for each sample (See Fig. 5). We processed the obtained radargrams using signal filtering operations, searching for the diffraction hyperbola, migration procedure to localize the source object and the Hilbert Transform to examine the energy of the reflections (See Fig. 6). The analysis of the hyperbolas has determined the relative permittivity of the studied soil and underlying metatuffs. It is 9 and 6, respectively. The electromagnetic wave length in the soil body is 5 cm, which implies the possibility to distinguish objects sized 1.5-5 cm. This range is consistent with the expected size of the apple tree taproots. After the processing of each profile, the roots were localized using hyperbolas and the surface of bedrocks was determined (See Fig. 7).

We found out that the soil profiles of the two experimental samples were different. In the first case, the average depth of bedrocks bedding is 60 cm, and there is a v-shaped decrease in the identified hyperbolas with depth. In the second case, the average depth is 45 cm and most of hyperbolas are located in the first 30 cm. The examination of the location of the roots in horizontal slices within an interval of 10 cm shows a change in the density of the roots distribution as we go deeper (See Fig. 8). The quantity of roots of Sample 1 gradually decreases with depth as well as the occupied area. In addition, the roots tend to concentrate near the trunk. Sample 2 has many roots in the upper part of the profile (0-30 cm) and they occupy a significant area.

However, the quantity of roots in the lower part decreases sharply. In both cases, we detected a predisposition of the root locations to the areas of the immersion of the roof of the bedrocks. In addition, a study was conducted concerning the intensity of the reflection amplitudes of both samples at different depths (See Fig. 9). At the near-surface level (20 cm) there are high-energy reflections. They occupy most of the studied area, and this corresponds to a horizon full of roots.

A different picture emerges at a depth of about 40 cm. The samples vary in both the reflection intensity and the occupied area. The revealed changes are due to the decrease in the roots quantity and apparently to the medium size of the roots.

We examined the correlation between the morphology of the bedrocks top surface and the spatial root architecture of the samples with three-dimensional models made in the form of a cloud of points as a set of established root locations in the XYZ-coordinate system (See Fig. 10).

Our studies have demonstrated the possibility to describe the tree root architecture by means of GPR. Using the two samples of the apple tree as an example, the structure of the root system in the soil-ground massif was determined. Also, we established the influence of the bedding roof depth of crystalline rocks on the ontogenesis of a tree.

The results obtained prove that GPR is a promising method for studying the underground areas of trees, which provide the possibility to analyze the root system without excavation. Such studies may prove useful in the future for both solving agricultural technology and forest management tasks, as well as for green spaces control in urban areas.

Key words: Malus domestica; ground penetrating radar; root architecture; radargrams; root impenetrable horizon.

Введение

Корневая система является критически важной составляющей любого дерева, отвечая за функции питания, запасания полезных веществ и физического закрепления в грунте. Под архитектоникой корневой системы подразумевается конфигурация корней в пространстве, а её формирование зависит от множества факторов, таких как вид растения, тип вмещающих грунтов, получаемое питание и т.д. [1]. Несмотря на важность изучения этого органа в контексте процессов развития и жизнедеятельности дерева, представление о структуре корневого аппарата зачастую неполное, в том числе и вследствие ограничений доступных методов исследований. Основой экспериментального изучения корневых систем растений служит способ их раскопки [2, 3]. Однако это наносит значительный ущерб подземному пространству, а также самим растениям и в большинстве случаев приводит к их гибели. В некоторых случаях в черте города или на особо охраняемых природных территориях выполнение подобных работ невозможно или запрещено. Кроме того, данный способ требует значительных временных и физических затрат.

Изучение корневых систем деревьев на сегодняшний день является актуальной задачей, и для ее решения систематически привлекаются новые дистанционные методы, как лабораторные, так и полевые [4]. Например, получен ряд значимых результатов по использованию пространственных моделей для оценки влияния наклона поверхности грунта на процессы формирования корневой системы в рамках борьбы с процессами эрозии [5, 6]. К неразрушающим способам анализа корневой архитектоники относится и применение геофизических методов, в частности, георадиолокации, которая базируется на изучении закономерностей распространения высокочастотных электромагнитных волн в почвенно-грунтовом массиве [7-9]. Георадиолокацию в качестве средства описания корневых систем древесных растений начали применять с конца 90-х гг., когда была показана принципиальная возможность локализации положения отдельных корней георадаром [10]. Дальнейшие исследования в этом направлении продемонстрировали достаточно убедительные примеры определения по данным георадиолокации диаметра скелетных корней и как следствие возможность оценки их биомассы [11-13]. При этом значительная часть работ касается вопросов физических ограничений и разрешающей способности георадиолокации [14, 15].

Корневая система растений развивается по наследственно закреплённой схеме, но имеет тенденцию к адаптации под влиянием окружающей среды, в том числе вследствие взаимодействия с вмещающими грунтами, что проявляется, например, в процессах гидротропизма и хемотропизма [16]. Одним из природных факторов, оказывающих влияние на рост, является уровень залегания корненепроницаемых горизонтов, вследствие чего изменяется как структура корневой системы, так и надземная часть растений.

Изучение развития корневой системы в почвенно-грунтовом массиве относится к области научных и практических интересов для плодовых растений, среди которых одно из доминирующих мест занимает яблоня домашняя (Malus domestica Borkh.) [17]. Цель данной работы - продемонстрировать возможность исследования корневой архитектоники взрослых деревьев методом георадиолокации без применения способа раскопки. Кроме того, поставлена задача по определению пространственной конфигурации корневой системы в зависимости от глубины залегания кристаллических пород и оценке влияния подобных изменений на развитие надземной части деревьев.

Материалы и методики исследования

Метод георадиолокации основывается на излучении наносекундных электромагнитных импульсов определенной частоты в зондируемую среду с последующей регистрацией и анализом сигнала, отраженного от внутренних границ и локальных неоднородностей [18, 19]. К основным характеристикам, отвечающим за скорость распространения и затухания электромагнитной волны в грунтах, относятся относительная диэлектрическая проницаемость и электрическая проводимость, которые в свою очередь зависят от ряда физических свойств зондируемого грунта: минерального и гранулометрического состава, влажности, плотности, содержания органического вещества и т.д. Вариации таких параметров определяют, насколько глубоко и как точно можно выявить отражающие границы в грунтовом массиве. Важным параметром в георадиолокации служит частота излучаемого сигнала, с которой связаны глубинность и разрешающая способность зондирования. Существует оптимальный диапазон частот, так называемое «георадарное плато», в интервале частот от 10 до 2 000 МГц [20]. В методе георадиолокации используются сверхширокополосные импульсы, частотный спектр которых в приёмнике представляет собой произведение исходного зондирующего сигнала на частотные характеристики среды. Таким образом, существует принципиальная возможность спектрального анализа регистрируемого сигнала для определения параметров зондируемой среды [19]. Например, повышение электрической проводимости грунта за счёт увлажнения приводит к увеличению низкочастотной составляющей сигнала вследствие затухания высокочастотных гармоник.

Для выполнения полевой съёмки используется георадар - прибор, который состоит из передающей и принимающей антенны с фиксированной центральной частотой, блока регистрации и обработки. Наземные георадио- локационные (георадарные) измерения проводятся вдоль профильной линии или по системе профилей. Конфигурация сети наблюдений зависит от характера, размеров и глубины исследуемых границ и объектов. Принимаемый сигнал регистрируется как функция изменения амплитуды отражённого электромагнитного импульса во времени - георадарная трасса, совокупность таких трасс формирует радарограмму. Прослеживая на радарограмме ряд одинаковых импульсов на соседних трассах, так называемую ось син- фазности, можно проследить отражающую границу (рефлектор). Кроме отраженных, существуют и другие типы регистрируемых волн, например дифрагированная. Её отличие в том, что электромагнитный импульс, отражённый от точки дифракции к поверхности, не подчиняется закону отражения, требующего равенства углов падения и отражения. Точками дифракции выступают локальные объекты: валуны, корни, подземные коммуникации, размеры которых сравнимы или меньше длины волны зондирующего сигнала [19]. На радарограммах дифрагированные волны выделяются по своему

годографу (зависимость времени прихода отражённой волны от пройденной дистанции) в форме гипербол, поиск и локализация которых является важной составляющей интерпретации (рис. 1). Вершина гиперболы указывает на положение искомого локального объекта, а по наклону «крыльев» гиперболы определяется скорость электромагнитной волны во вмещающей среде, что важно в случаях отсутствия информации о строении и электрофизических параметрах изучаемого разреза.

Рис. 1. Схема выполнения георадарной съемки (Tx - источник, Rx - приемник) и регистрируемая радарограмма, где: 1 - прямая волна в грунте; 2 - гипербола дифрагированной волны на локальном объекте; 3 - рефлектор от подстилающих пород [Fig. 1. Schematic illustration of GPR reflection profiling (Tx - Transmitter, Rx - Receiver) and recorded radargram, where: 1 - Direct ground wave; 2 - Hyperbola of the diffracted wave; 3 - Bedrock reflections]

Для относительно низких георадарных частот разрешающая способность составляет десятки сантиметров, тогда как для самых высоких разрешение - первые сантиметры. Например, заявленная характеристика для антенного блока с центральной частотой 150 МГц составляет 35 см, а для антенны с частотой 1 500 МГц равна 5 см. Точность определения положения отражающего объекта также зависит от шага измерения по профилю и дискретизации сигнала.

Георадиолокация при исследовании корневых систем деревьев. С точки зрения георадиолокации древесные корни выступают в качестве локальных линейных объектов, т.е. их поисковый признак, как уже отмечалось выше, - гипербола дифрагированной волны. Относительная диэлектрическая проницаемость древесины зависит от ряда факторов, основные из которых водонасыщенность и плотность, при этом существует зависимость от температуры, частоты импульса, содержания целлюлозы и т.д. Как показывают лабораторные исследования, нормальное значение для сухой (влажность менее 30%) древесины є<10, тогда как древесина с влажностью порядка 80-90% может достигать значений є = 50 и более [21]. Для сравнения: нормальные значения песка составляют 4-7 единиц в сухом состоянии и 5-15 во влажном [19]. Следовательно, существующий контраст свойств обеспечивает формирование рефлекторов. Особенно ярко это проявляется в засушливых регионах, где корневые системы растений более мощные и занимают большую территорию, а также активней депонируют воду на фоне сухих почв [22]. Известно, что отчётливо гиперболы на радарограмме проявляются, когда влажность корней на 35% (и более) выше относительно фоновых значений, при разности менее 20% корни не выявляются [23, 24].

Лабораторное моделирование позволило определить, как изменяется волновое поле георадиолокации в зависимости от диаметра и глубины залегания корня, а также выявить возможные подходы для оценки диаметра корней в естественном залегании исходя из параметров георадарного сигнала [15, 24-26]. При этом широкая вариация диаметра и протяжённости корней у деревьев затрудняет получение однозначных критериев для их идентификации. В табл. 1 представлены обобщённые данные по изучению корневых систем методом георадиолокации. Сопоставление результатов исследований показало существование вариативности как по выявленному диаметру, так и максимальной глубине, что обусловлено видом деревьев и типом вмещающего почвенно-грунтового массива. Отмечается, что минимальный диаметр корней, однозначно фиксируемый высокочастотными антенными блоками, составляет 0,5 см [12, 15, 25]. Дополнительно к факторам, определяющим волновую картину, относится угол пересечения георадарным профилем направления горизонтального роста корня, при углах <45° и >135° регистрируемая гипербола искажается, а амплитуда отраженного сигнала снижается [27]. Следовательно, при исследованиях необходимо подбирать геометрию профилей наблюдений в зависимости от конфигурации корневой системы [28], например, использовать прямоугольную [29] или радиальную [30] сеть профилей.

Выявлено влияние сезона на детальность георадарных исследований, в теплый (летний) период выделяется большее число корней на глубинах до 0,5 м, тогда как в холодный (зимний) период возрастает в 2 раза и более количество корней, локализованных ниже 0,5 м от поверхности [31]. Это явление объяснено повышенным содержанием воды в корнях летом и как следствие ростом их контраста в сравнении с вмещающими почвами. В свою очередь зимой происходит возрастание скорости электромагнитной волны в совокупности с уменьшением затухания в почвах из-за промерзания (т.е. происходит увеличение эффективной глубины метода). Также точность идентификации корней зависит от мощности и состава лесной подстилки [32]. По существующим оценкам, при помощи георадара можно выявить расположение порядка половины корней из общей структуры корневой системы дерева - 48% [33] и 58% [34]. При этом в верхнем интервале <50 см выявленных корней больше в сравнении с нижележащими горизонтами - 49 к 35% [33] и 60 к 35% [34].

Экспериментальная площадка и образцы деревьев. Исследования проведены в августе 2019 г. на территории Ботанического сада Петрозаводского государственного университета, расположенного в г. Петрозаводске (рис. 2). В геологическом плане эта территория представляет собой ось антиклинальной складки, проявленной в виде мощного крутосклонного кряжа субширотного простирания с изрезанным рельефом. Для рассматриваемого участка характерно близкое залегание кристаллических пород, которые представлены порфиритами и метатуфами. Крутой рельеф не способствовал аккумуляции рыхлых отложений, при этом в самом почвенном покрове содержатся обломки скальных горных пород. В четвертичных отложениях территории преобладает супесчаная и песчаная морена с мощностью не более 2 м, а значительные площади вообще лишены четвертичного покрова [35]. Структура почвенного покрова Ботанического сада неоднородна и усложняется с запада на восток, основная масса почв представлена буроземами грубогумусовыми супесчаными. Отличительными чертами этих почв являются: наличие гумус-продуктивного горизонта с содержанием углерода 7-12%, далее углерод убывает с глубиной до значений менее 1%, средний гранулометрический состав - 87% песка и 13% глины, весь профиль ожелезнсн и обогащен дериватамииз подстилающих пород [36].

Рис. 2. Расположение участка работ [Fig. 2. Location of the site research area]

Для эксперимента выбраны две площадки на расстоянии около 500 м, на каждой из которых росли одновозрастные деревья - яблони домашние Malus domestica Borkh. Условия произрастания яблонь, такие как почвенный состав, количество осадков, температура, освещенность и уровень ветра, равноценны. На каждой площадке эксперимента выбрано отдельно стоящее дерево со средними значениями биометрических показателей (табл. 2).

Таблица 1 [Table 1] Эффективные диаметр и глубина обнаружения корней методом георадиолокации [Effective diameter and depth of root detection by GPR]

При этом плотность и шаг посадок на этих площадках исключали срастание корневых систем у рядом стоящих деревьев. Выбор этих двух конкретных объектов обусловлен внешними отличиями. Одно из растений имеет выраженную сложную кривизну ствола и меньшую высоту по сравнению со вторым, условно-нормальным (рис. 3).

Рис. 3. Образцы обследуемых яблонь: a - № 1 с первой площадки;

b - № 2 со второй площадки (фотографии А.В. Кабонена) [Fig. 3. Test samples of apple trees: a - from Plot 1, b - from Plot 2 (Photos by Alexey Kabonen)]

Таблица 2 [Table 2]Биометрические показатели исследуемых деревьев [Biometric indicators of the test trees]

Среднегодовые фенологические фазы у исследуемых образцов проходили параллельно, возраст деревьев - 29 лет, значительной разницы в количестве ежегодно вызревающих плодов у исследуемых растений не выявлено. На базе существующих отличий выдвинута гипотеза о влиянии глубины залегания скальных пород на формирование корневой системы и как следствие на процессы развития ствола и кроны дерева. Известен эффект специфического протекания процессов корнеобразования в зависимости от почвенногидрологических условий.

В случае почвенного покрова малой мощности и наличия корненепроницаемого горизонта происходит прекращение развития стержневого корня с интенсификацией роста корней горизонтальной ориентации, структура корневой системы принимает выраженный поверхностный характер, зачастую имея форму усеченного конуса [1]. Кроме того, считается, что невозможность проникновения вертикальных корней вглубь не только меняет структуру корневой системы, но может повлечь за собой изменения и слабый рост надземной части, а также раннее отмирание плодовых деревьев [39].

По литературным данным, условно-нормальная корневая система яблони домашней имеет вид, как показано на рис. 4. Отмечается преобладание скелетных корней горизонтального направления и расположение их основной массы на глубинах до 1 м от уровня поверхности. Горизонтальные скелетные корни характеризуются диаметром около 1-2 см, а вертикальные - порядка 4-5 см. Представленная схема описывает общий случай, тогда как для северных территорий (к которым относится и Республика Карелия) характерно распределение основной массы корней в интервале 25-50 см от поверхности [Там же]. Имеющиеся предпосылки позволяют сделать вывод о применимости георадара для сбора данных, так как глубина залегания корневой системы у яблони домашней соотносится с глубиной зондировании высокочастотного антенного блока (см. табл. 1), а также предположить высокую информативность метода из-за большого количества субгоризонтальных корней, которые проще выявлять на радарограммах.

Рис. 4. Горизонтальная и вертикальная проекции корневой системы яблони домашней возрастом 20 лет, по данным В.А. Колесникова (1974) [Fig. 4. Horizontal and vertical plane of the apple tree root system aged 20 years, according to Venedikt Kolesnikov (1974)]

Методика эксперимента. Полевые наблюдения проведены исходя из имеющихся данных о строении корневой системы яблони домашней. Для этого вида характерно наличие в горизонтальной проекции четырех перпендикулярных и четырех диагональных скелетных корней, а также распределение корней в вертикальной проекции на глубине от 20 до 60 см (см. рис. 4). Вследствие этого измерения выполнены по регулярной прямоугольной сети с расстоянием между профилями 50 см (рис. 5). Каждый профиль начинается с отметки 0 см и увеличивался до пикета 600 см. Основная масса корней у яблони домашней залегает в пределах кроны, поэтому размер каждого участка съемки составил 600*600 см. Сбор данных выполнен георадаром ОКО-2 «Логис-Геотех» (Россия) с антенным блоком, имеющим центральную частоту 1 700 МГц. Указанный антенный блок имеет фиксированную базу, шаг сканирования по профилю составил 2 см, а временная разверстка 32 нс, что обеспечило глубину исследований до 1 м. В результате для каждого образца получен набор из 26 профилей: 13 горизонтальных (направление X) и 13 вертикальных (направление Y).

Обработка полученных радарограмм выполнена стандартно [22]: частотная фильтрация и вычитание средней составляющей сигнала для чистки записи; выделение областей расположения полезного сигнала; поиск гипербол дифрагированной волны; выполнение процедур синтеза апертуры и преобразования Гильберта (рис. 6). Операция синтеза апертуры повышает разрешающую способность записи и уменьшает интенсивность «крыльев» гипербол, что обеспечивает лучшую локализацию местоположения отдельных корней. Преобразование Гильберта позволяет представить георадар- ную трассу как аналитический сигнал (огибающую), который характеризует мгновенные амплитуды и позволяет оценить энергию (интенсивность) отражений в конкретной области. Участки высоких энергий соответствуют контрастным изменениям электрофизических свойств. Из практики известно, что достаточно крупные корни могут отображаться в волновом поле георадиолокации набором отражений, например от верхней и нижней части, поэтому их положение корректно определять по наибольшим амплитудам [15, 23].

Для перехода от временной формы георадарной записи к разрезу глубин необходимо присвоить относительную диэлектрическую проницаемость є для каждого отдельного слоя среды, в которой распространялся сигнал. Используя формулу (3), определить скорость V, а по формуле (1) рассчитать глубину зондирования. Величина є (или V) определяется по табличным значениям с использованием маркерных рефлекторов с известной глубиной залегания или рассчитывается по наклону «крыльев» гипербол дифрагированной волны [18-20]. Последний способ рассматривается как наиболее удобный, так как не требует дополнительных работ и обеспечивает достаточную точность определения искомых значений. Условием его применимости является наличие на записи чётких гипербол. В рамках исследований корней использована программа сбора и обработки данных георадиолокации GeoScan32 «Логис-Геотех» (Россия), где реализована функция анализа гипербол. По результатам работ для вмещающих почв є составила 9 единиц, а для подстилающих метатуфов є = 6 единиц. Это позволило определить скорость У_поче ~ 10 см/нс и глубину зондирования Z = 70-90 см. Значение К_тр (4) для границы почва - скальный грунт составляет 0,1, что свидетельствует о существовании рефлектора средней интенсивности. Длина зондирующей электромагнитной волны в исследуемых почвах составляет ~5 см (5), это позволяет сделать вывод о возможности локализации объектов размерами от 1,2 до 5,0 см (6). Указанный диапазон сопоставим с ожидаемым размером скелетных корней яблони домашней.

Рис. 5. Схема расположения профилей наблюдения в обследуемой области [Fig. 5. Setup of GPR profiles in the study area]

После обработки на каждом профиле по гиперболам определены местоположения корней с координатами X, Y, Z. Так как нет однозначных параметров, чтобы отличить отражение от корня и, например, валуна, пикировались все гиперболы. Предполагается, что из-за протяженности корни имеют подавляющее количество отражений на радарограммах, а отражения от иных объектов статистически не значимы. В силу различных факторов существует неопределенность для параметров сигнала, поэтому нет возможности достоверно проследить единичный корень. Вследствие чего выполнен анализ конфигурации некоторого облака точек, которое в первом приближении пространственно соответствует объёму среды, занимаемому корневой системой. Кроме местоположений корней определена поверхность скальных грунтов, контрастно выделяющаяся от вышележащих рыхлых грунтов по протяженной оси синфазности.

Рис. 6. Пример процедур обработки радарограммы: a - поиск гипербол; b - синтез апертур; c - расчет энергограммы (стрелками указано предполагаемое местоположение корней)

[Fig. 6. An example of GPR transect processing procedures: a - Search of hyperbolas; b - Aperture synthesis; c - Calculation of signal amplitudes (arrows indicate the location of the roots)]

Результаты исследования и обсуждение

Выполненные работы обеспечили получение набора радарограмм, по которым локализованы гиперболы дифрагированной волны, а также прослежена граница кристаллических пород. Рассмотрим результаты на примере двух перпендикулярных профилей, расположенных в центре участка для каждого из образцов (рис. 7). Радарограммы для образца № 1 характеризуются хаотичным волновым полем со временами регистрации полезного сигнала до 15 нс. Граница скальных пород идентифицируется по невыдержанной оси синфазности и изменению характера волнового поля - отсутствие интенсивных отражений и повышение затухания сигнала. Для профиля по направлению X (рис. 7, а) глубина залегания скальных пород составляет 50-72 см, гиперболы фиксируются на пикетах от 25 до 454 см и в интервале глубин (Z) 2,5-61,7 см. Профиль по направлению Y (рис. 7, b) характеризуется сходной глубиной залегания скальных пород 54-70 см, гиперболы выявлены на пикетах 22-544 см и в интервале Z = 2,7-61,7 см. Далее обратимся к результатам исследованиям образца № 2. Радарограммы на этом участке имеют меньший временной период полезного сигнала порядка 12 нс. Граница скальных пород для профиля по направлению X установлена в интервале глубин 32-53 см, гиперболы выявлены в диапазоне пикетов 30-450 см и Z = 2,3-46,0 см (рис. 7, с). Перпендикулярный профиль Y (рис. 7, d) обладает схожими параметрами: скальные породы прослежены на глубине 40-52 см, гиперболы на пикетах 76-518 см и Z = 3,0-42,7 см.

Рис. 7. Центральные радарограммы для образца № 1 по направлению X (a), Y (b) и образца № 2 по направлению X (с), Y (d)

[Fig. 7. Central GPR profiles for sample 1 in the direction X (a), Y (b) and sample 2 in the direction X (c), Y (d)]

В общем виде можно заключить, что почвенно-грунтовые условия объекта № 1 и объекта № 2 кардинально отличаются. В первом случае средняя глубина залегания кристаллических пород составляет 60 см и наблюдается конусообразное уменьшение количества выявленных гипербол с глубиной. Во втором случае средняя глубина равна 45 см, а основная масса гипербол расположена в первых 30 см. Отметим, что на участке образца № 1 обнаружено 366 гипербол, а на участке образца № 2 - 520 гипербол. Так как скелетные корни из-за своей протяженности формируют подавляющее количество гипербол во всех проекциях, гиперболы приняты в качестве признаков их местоположения. Как отмечалось выше, в среднем регистрируется 40-50% корней от общего количества, поэтому пространственная структура корневых систем обследуемых яблонь описывается достаточно полно.

Рис. 8. Распределение выявленных местоположений корней по глубине [Fig. 8. Distribution of the identified root locations in depth]

Распределение количества выявленных местоположений корней в зависимости от глубины с интервалом 10 см показано на рис. 8. Видно, что корневые системы исследуемых деревьев имеют ряд существенных отличий. Образец № 1 демонстрирует постепенное уменьшение количества корней с глубиной и согласное уменьшение занимаемой площади, также наблюдается их концентрация относительно штамба. Образец № 2 характеризуется множеством корней в верхней части разреза (0-30 см), занимающих значительную площадь, тогда как в самой нижней части их количество резко сокращается. Следует отметить, что распределение корневой массы по глубине для образца № 1 сопоставимо с оценкой В.А. Колесникова для яблонь северных территорий. В то же время в обоих случаях выявлена предрасположенность местоположений корней к областям погружения кровли кристаллических пород.

Один из способов оценить пространственную конфигурацию погребенного объекта - использовать проекции (срезы георадарных параметров на фиксированных глубинах. Для получения подобных срезов рассчитываются значения мгновенных амплитуд отраженного сигнала (энергограммы), что позволяет проследить положение электрофизических неоднородностей. Корни за счет своей водонасыщенности контрастны на фоне вмещающих грунтов, и это обеспечивает трассирование их положения за счет высокой интенсивности отражения сигнала. Подобный способ оценки использован для сосны ладанной (Pinus taeda), когда выявлялись корни диаметром более 2 см, и по результатам составлены карты их расположения [40]. В рамках выполненных исследований использовался схожий подход. Для обследуемых образцов получены проекции интенсивности амплитуды отражений на разных интервалах глубин (рис. 9). Следует отметить, что полученные результаты, в отличие от литературных данных, не позволяют показать отдельные корни в связи с отсутствием однозначных критериев для прослеживания, однако определенную полезную информацию содержат. Так, например, уровень Z = 20 см характеризуется наличием высокоэнергетических отражений, которые располагаются в большей части исследуемой площади, и это соответствует горизонту, наиболее насыщенному корнями. Для обследуемых растений картины схожие, площадь, в пределах которой установлены интенсивные отражения, для образца № 1 составляет 28,7 м², а для № 2 - 26,6 м². Иная картина наблюдается на глубине Z = 40 см: образцы отличаются как по интенсивности наблюдаемых отражений, так и по занимаемой площади. Выявленные изменения обусловлены уменьшением количества и, вероятно, среднего размера корней в этом интервале. Выделяется значительная разность площади, занятой интенсивными отражениями: 25,6 м² для образца № 1 и 16,9 м² для образца № 2. Данные факты свидетельствуют как о трансформации корневой системы яблони домашней с глубиной, так и об отличии архитектоники корневой системы рассматриваемых образцов из-за влияния почвенно-грунтовых условий.

Взаимоотношение морфологии кровли кристаллических пород и корневой системы исследуемых образцов хорошо наблюдается на трехмерных моделях. В рассматриваемом случае модели выполнены в виде облака точек - набора установленных местоположений корней в координатной системе X, Y, Z, совокупность которых, несмотря на неточность и приблизительный характер, в первом приближении позволяет пространственно отобразить корневую архитектонику исследуемых образцов (рис. 10). На текущем этапе исследований проследить направление и размеры отдельного скелетного корня по набору выявленных местоположений не удалось. В перспективе есть предпосылки для определения геометрических параметров корней в естественном залегании. Отдельные примеры показывают такую возможность, но в частном случае изучения якорных корней крупных деревьев с диаметром ствола 1 м и высотой до 20 м [41].

Рис. 9. Проекции амплитуд отражений на глубинах 20 и 40 см:

a - образец № 1; b - образец № 2 (кругом обозначено положение ствола) [Fig. 9 Projections of the reflection amplitudes at depths 20 and 40 cm: a - Sample 1; b - Sample 2 (the circle indicates the position of the tree trunk)]

Вместе с тем на основе модельных наблюдений экспериментально доказано наличие значимой корреляции между диаметром корней и характеристиками георадарной трассы: временных интервалов отражений [14], значений максимальной амплитуды [42], а также взаимоотношений этих параметров, так называемых георадарных индексов [11, 25-29]. Данные исследования показывают размерность коэффициента корреляции георадарных индексов с диаметром и биомассой корней в пределах r = 0,5-0,8. Однако на практике использование такого подхода затруднительно, а результаты крайне неоднозначны. Кроме того, обработка большого количества полевых наблюдений из-за трудоёмкости требует привлечения специализированных автоматических процедур обработки радарограмм.

Рис. 10. Трехмерные модели корневой архитектоники: a - образец № 1; b - образец № 2 [Fig. 10. Three-dimensional models of the root architecture: a - Sample 1; b - Sample 2]

Полученные данные хорошо согласуются с выдвинутой гипотезой - неглубокое залегание скальных пород приводит к деформации конфигурации корневой системы, основная масса корней начинает концентрироваться в поверхностных слоях.

Теоретически это может приводить к развитию сложной кривизны древесного ствола, обусловленной фактором ветровальности или неприспособленности деформированной корневой системы поддерживать формирующуюся по наследственно-закрепленной схеме надземную часть растений, вследствие чего при условно-нормальном росте надземной части деформированная корневая система не в состоянии удерживать рост ствола в вертикальном направлении.

Такие внешние условия приводят к наклонности и сложной кривизне, что и наблюдается визуально для образца № 2. Следует отметить, рассмотрение всего комплекса факторов, определяющих формирование древесного ствола, и обсуждение степени их значимости выходит за рамки представленных исследований.

Заключение

Выполненные экспериментальные работы продемонстрировали возможность описания корневой архитектоники дерева с использованием метода георадиолокации. Обобщенные данные показали, что предельный размер отдельных корней, возможных для определения, составляет 0,5 см, этот факт позволил в дальнейшем рассматривать скелетные корни как объект поиска. На примере изучения двух образцов яблони домашней проведено картирование структуры распределения корневой системы в почвенно-грунтовом массиве, также установлено влияние глубины залегания кровли кристаллических пород на онтогенез деревьев. По результатам исследований можно заключить, что георадиолокация является многообещающим методом изучения подземной части деревьев, обеспечивающим выполнение анализа корневой системы без экскавации. В перспективе подобные исследования имеют практическую значимость как для решения агротехнических и лесоустроительных задач, так и при контроле зеленых насаждений на урбанизированных территориях.