Pool of Lignin Phenols in Soils of Secondary Forests

Full Text

ВВЕДЕНИЕ

Масштабное лесовосстановление и лесоразведение в разных регионах мира в части реализации Парижского соглашения 2015 г., ратифицированного 189 странами [4], направлено на ускорение стока и сокращение выбросов углерода. Главным планетарным депо углерода при этом считаются тропические и бореальные леса. Тропические леса называют “легкими планеты”. Между тем, известно [1], что лесные экосистемы средней и южной тайги, широколиственных лесов, тропических и субтропических лесов были практически полностью уничтожены по всему миру и продолжают сокращаться из-за лесных пожаров, интенсивного развития сельского хозяйства, городского строительства. Часть вырубленных лесных площадей возобновилась, бывшие пахотные почвы зарастают [19]. Поэтому современные лесные угодья в большинстве стран мира имеют не коренной, а вторичный характер.

Однако механизмы трансформации различных пулов органического вещества и стабилизации углерода в почвах при вторичном лесовосстановлении остаются малоизученными [22, 23, 32]. Наиболее устойчивой к разложению частью лесного опада является лигнин. Древесина хвойных деревьев содержит 25–33% лигнина, лиственных – 18–25% лигнина; травянистые растения накапливают 4–9% [13]. Сложная биохимическая структура биополимера с большим количеством ароматических и гидрофобных структур, а также его высокая устойчивость к микробной трансформации определяют долговременную сохранность лигнина в почвах и уникальную роль в процессах депонирования углерода в экосистемах [7, 15, 29]. Процессы разложения в почвах происходят за счет трансформации боковых цепочек биополимера микроорганизмами с образованием промежуточных нестабильных радикалов. За счет разрыва связей метоксильных групп, добавления гидроксильных и иных полярных групп в молекулу лигнина, а также доступного кислорода или воды компактные гидрофобные структуры лигнина сильно разрыхляются, увеличивая его растворимость. Считается доказанным [10], что множественная, постепенная и продолжительная трансформация обеспечивает участие лигнина в образовании почвенного гумуса, а не реполимеризация отдельных фенольных мономеров в гуминовые кислоты, как предполагалось ранее.

Однако до сих пор неясно, какую роль играет лигнин в процессах секвестрации и депонировании органического углерода почвой при вторичном лесовосстановлении, какое характерное время необходимо лесной экосистеме для стабилизации поступающего в почву органического углерода [15, 24].

Цель работы – характеристика пула лигниновых фенолов в почвах двух контрастных экосистем вторичных лесов: бореальных и тропических под разными породами деревьев и оценка скорости биохимической трансформации трудноразлагаемых соединений лигниновой природы в них.

ОБЪЕКТЫ И МЕТОДЫ

Почвенный покров экспериментальных площадок, заложенных в 1971 г. по инициативе Орловского [3] в южной тайге (50 км к северу от Красноярска) для изучения внутренних связей в системе почва–растение, представлен серыми почвами (Abic Greyzemic Phaozems [12]). Климат региона умеренно-континентальный с температурой июля +20°С и количеством осадков 500 мм. Территория “столетнего”, как называл его сам Н.В. Орловский, опыта площадью 1.7 га приурочена к обширной древней террасе р. Кача, на которой после раскорчевки вторичного березняка и использования земли в сельскохозяйственных целях под посадку картофеля и сенокос, в течение двух лет почва оставалась под паром. Затем на участках размером около 24 сотки каждый были высажены 2–3-летние саженцы основных лесообразующих пород Сибири: ели сибирской (Picea obavata Ledeb.), березы кустарниковой (Betula fruticosa Pall.), сосны обыкновенной (Pinus sylvestris L.), осины обыкновенной (Populus tremula L.), лиственницы сибирской (Larix sibirica Ledeb.) и сосны (кедровой) сибирской (Pinus sibirica Du Tour). Контрольный участок размером 96 соток использовался в качестве пастбища под травяной ассоциацией. Схема подготовки эксперимента описана ранее [5, 25] и включала срезание гумусового слоя в целях выравнивания площадки с его последующим складированием в бурты, рыхление минеральной массы на глубину до 70 см и равномерную обратную засыпку гумусированным материалом, исключая площадки под травой. Образцы почв для анализа отбирали с глубины 0–10 см.

Второй объект исследования представляет собой экспериментально-производственные площадки научно-исследовательской станции Embrapa в бассейне р. Амазонки недалеко от Манауса, Бразилия (3°8’ S, 59°52’ W, 40–50 м над ур. м.). Ежегодное количество осадков 2622 мм, температура воздуха 26°С и влажность около 85% (средние значения за 1971–1993 гг.) [27]. Самые засушливые месяцы с июля по сентябрь, наиболее влажные – с февраля по апрель. В соответствии с Классификацией почв WRB [12], исследуемые объекты представляют собой ферралитные почвы (Ксантиковые Ферральсоли, Eutric Xanthic Ferralsols). Содержание питательных веществ, по данным [28], обычно низкое, pH верхнего слоя почвы составляет около 4.5. Содержание глины до 800 г/кг, почвы имеют очень низкую емкость катионного обмена в естественном состоянии, высокое насыщение алюминием. Содержание подвижного доступного фосфора – очень низкое. Общее содержание C и N уменьшается вниз по профилю, но все еще обнаруживается на глубине 2 м.

В исследование были включены следующие системы плантаций [27, 28]: купуаку (Cupuacu, близкий родственник какао, который выращивают из-за ценной для соков фруктовой мякоти), бразильский орех урукум (Urucum), возделываемый из-за нетоксичного красного красителя и превосходной древесины. Кроме того, в исследование вошли и два вида деревьев из ближайшего первичного тропического леса: Eschweilera sp. (двудольное дерево “Мата-мата”) и Oenocarpus bacaba (пальма “Bacaba”), а также плантация производственной древесины Андрироба (Andriroba). Экспериментальные площадки, за исключением участков тропических лесов, были организованы в виде рандомизированного полного блока с пятью повторностями. Для коренных лесных видов было выбрано по три особи каждого вида в лесу, прилегающем к экспериментальным площадкам. Размер участка составлял 24 × 32 м для монокультур персиковых пальм и 48 × 32 м для остальных пород. Особый интерес представляет участок вторичного леса: место исследования впервые очищено от девственного леса в 1980 г. с использованием тяжелой техники для валкования и удаления пней, в 1981 г. на нем был поставлен эксперимент с каучуковыми деревьями (Hevea brasiliensis), который был заброшен в 1986 г. Развивающийся вторичный лес был очищен вручную в 1992 г., при этом на участке преобладают виды Vismia, которые являются характерными для растительности молодых залежей и деградированных земель в районе исследований.

Для определения содержания лигнина в почвах использовали методику Хеджеса, Манна [11] в модификации Амелунг [4]. Она включала окисление образца оксидом меди в щелочной среде при температуре 170°С под давлением после предварительной дегазации смеси азотом; перевод гуминовых кислот в осадок; использование колонок С18 для концентрирования фенольных продуктов, дальнейшее высушивание колонок и растворение лигниновых фенолов в этилацетате с последующим эвапорированием этилацетата. Разделение составляющих лигнин фенолов выполняли на газовом хроматографе HP6890 GC (Hewlett-Packard Palo Alto CA, США), оснащенном пламенно-ионизационным детектором с кварцевой капиллярной колонкой 30 м × 0.25 мм с силиконовой фазой OVS-1. В качестве основного и маркирующего газа служил азот. Температура инъектора – 250°С, температура детектора – 300°С. Для газохроматографического анализа ароматические альдегиды, кетоны и кислоты переводили в триметилсилиловые эфиры растворением полученного сухого экстракта в пиридине и BSTFA при комнатной температуре в течение 1–2 ч (время для полного силилирования подобрано экспериментально). Несмотря на то, что чувствительность метода мягкого щелочного окисления очень высока – позволяет определять микроколичества фенолов – могут быть потери продуктов разложения лигнина в процессе анализа.

Для повышения воспроизводимости результатов анализа в почвенные пробы перед щелочным окислением добавляли глюкозу в качестве катализатора окисления, внутренним стандартом выступал этилванилин, его выход в результате анализа принимался за 100% [6]. Для улучшения контроля за ходом анализа в качестве второго внутреннего стандарта за воспроизводимостью самого этилванилина перед дериватизацией в пробы добавляли фенилуксусную кислоту. Достоверными при этом двойном контроле считались результаты, ошибка опыта для которых составляла не более 5–10% (воспроизводимость по отношению к этилванилину 90–95%), остальные выбраковывались.

В результате анализа растительных тканей в почве идентифицированы фенолы [9], которые по своему биохимическому строению объединяются в четыре группы: ванилиновые (V), сирингиловые или сиреневые (S), циннамиловые – кумаровые (С) фенолы. Ванилиновые и сирингиловые фенолы распознаются в ходе анализа в форме альдегидов (al), кетонов и кислот (ac), а циннамиловые фенолы имеют только кислотные формы: феруловых и кумаровых кислот. Полученная сумма продуктов окисления лигнина (VSC) – это фенольные структуры без каких-либо изменений в ароматических частях молекул.

Содержание общего углерода в отобранных образцах почв определяли на CNS-анализаторе Vario EL, Elementar GmbH, Hanau. Статистическую обработку результатов выполняли c использованием программ Excel 2010 и Statistica 20 (StatSoft, Inc. США). Полученные данные представлены в виде среднего ± стандартное отклонение. Расчет величины t-критерия для парных различий подтвердил значимость различий при р < 0.05 в величинах содержания суммы фенолов между вариантами вторичного леса и контролем (почва площадки без древостоя).

Для оценки степени трансформации биополимера в почвах (Тr, %) авторы [9] предложили формулу:

$Тr=74–(100–K)(1+{(Ас/Аl)v)}^{–1}$,

где (Ас/Аl)v – отношение содержания ванилиновых кислот к ванилиновым альдегидам (v – ванилиновые структуры), K – количество кетонов в исходных растительных тканях (19%). Показано, что по мере разложения биомолекул лигнина количество ароматических кислот в почве нарастает по отношению к ароматическим альдегидам.

РЕЗУЛЬТАТЫ

Исследуемые серые почвы вторичных лесов (табл. 1) характеризуются разным содержанием органического углерода в гумусовых горизонтах: наибольшим (5.60–5.93%) в почвах площадок с лиственницей и травами, средним (3.04–4.14%) в почвах площадок с елью сибирской и сосной сибирской, наименьшим (2.65%) в почвах под березовой. Отношение C/N изменяется от очень низкого в почвах под березой и осиной (37.9–41.1), до низкого – в почвах под хвойными деревьями и травами (15.9 и 19.1 соответственно). Реакция среды почв (pH_H2O) нейтральная или близкая к нейтральной.

Таблица 1. Химические свойства и содержание лигнина (VSC) в серых почвах (Phaozems Albic)

* Содержание лигнина (суммарное количество продуктов окисления лигнина (VSC), мг/г С_орг). Стандартное отклонение (Х ± SE), n = 3.

Максимальное содержание лигнина демонстрируют серые почвы под елью (4.98 ± 3.1 мг/г С_орг), а также – под лиственницей (5.18 ± 0.3 мг/г С_орг). Минимальным накоплением лигнина отличаются серые почвы площадок под сосной обыкновенной (1.52 ± 0.3 мг/г С_орг) и под осиной (1.57 ± 0.2 мг/г С_орг).

Характеристики пула лигниновых фенолов в серых почвах вторичных лесов (табл. 2) подтверждают известные закономерности [31]. Так, содержание ванилина и ванилиновой кислоты в почвах под голосеменными деревьями, особенно под елью и сосной, превышает количество ванилиновых фенолов в почвах под покрытосеменными лиственными растениями – осиной и березой, а также в почвах под травяными ассоциациями.

Таблица 2. Содержание продуктов окисления лигнина (VSC, мг/г С_орг) в серых почвах

Примечание. Стандартное отклонение (Х ± SE), n = 3.

Количество сирингиловых фенолов близко к содержанию ванилиновых фенолов в почвах лиственных делянок (площадки под осиной и березой). При этом содержание циннамиловых (коричновых) фенолов, таких как п-кумаровая и феруловая кислоты, больше в почвах площадки без древостоя по сравнению с почвами под древесными насаждениями. Исключение составляют почвы делянок под елью сибирской и лиственницей. В них высокое по сравнению с контролем содержание кумаровых и феруловых кислот синхронно максимальному содержанию общего лигнина (VSC) – 4.98–5.18 мг/г С_орг и, вероятно, обусловлено биохимической “памятью” почв: сохранением молекулярных остатков трав вследствие замедленной микробиологической активности в верхнем слое почвы (0–10 см), а именно, преобладанием грибной микрофлоры. Известно [13], что в почвах под лесом в разложении лигнина принимают участие грибы бизидиомицеты (белой и бурой гнили). А для эффективной жизнедеятельности микроорганизмов нужна аэробная среда и низкое содержание фенолов – не более 0.025–0.05% для бактерий и не более 1% – для плесневых грибов. Накопление лесной подстилки на протяжении длительного времени в естественных бореальных лесах приводит к накоплению значительных количеств лигниновых фенолов в почве и формированию анаэробных условий в нижних слоях подстилки [8]. В то же время на открытых аэрируемых пространствах под травянистой растительностью процесс минерализации лигнина ускоряется за счет усиления биоактивности почв [19].

Поэтому в освоенных почвах после сведения леса наблюдается значимое (при р <0.05) снижение содержания суммы лигниновых фенолов VSC, повышение степени окисления биополимера (As/Al)v и высокие величины параметра трансформации его боковых цепочек Тr. Содержание лигнина в освоенных почвах, как было показано ранее на примере исследования почв на месте сведенного леса на Русской равнине, в 2–3 раза уменьшается [30, 16]. Однако подтверждение этой гипотезы на данных объектах требует дополнительных исследований.

Для оценки степени податливости растительных остатков минерализации используют отношение VSC/N. Широкая величина отношения лигнин/азот типична для ароматического (структурного) углерода, трудно извлекаемого микроорганизмами, а узкое отношение – для метаболического, доступного углерода, быстро вовлекаемого в процессы метаболизма [26].

Все исследованные почвы характеризуются значительным обогащением органического вещества метаболическим углеродом: соотношение VSC/N в них в 2–3 раза уже, чем в почвах естественных и освоенных экосистем Русской равнины [16], где оно достигает величин 67–84.

Полученные результаты позволяют предположить разный характер минерализации растительных тканей в почвах разных экосистем. Например, органическое вещество почвы площадок с сосной обыкновенной характеризуется наибольшей обогащенностью метаболическим углеродом (VSC/N 8.4). Близкие значения были обнаружены только в почвах под злаковой травянистой растительностью и под сосной сибирской. Важно подчеркнуть, что почвы под сосной обыкновенной характеризуются наименьшим содержанием лигнина, а также наибольшей степенью окисления (Ac/Al)v и степенью трансформации боковых цепочек биополимера Tr, достигающей 51.4% (табл. 3).

Таблица 3. Параметры биотрансформации лигнина в серых почвах (n = 3)

Примечание. Здесь и далее. Кислоты/альдегиды, (Ac/Al)v – отношение количества фенольных кислот к альдегидам в ванилиновых единицах; S/V – сирингиловые фенолы/ванилиновые фенолы; C/V – коричные фенолы/ванилиновые фенолы; K/F – кумаровая кислота/феруловая кислота; T_r – процент изменений в боковых цепях от соотношения лигнина к исходной растительной ткани. V : S : C, где V – ванилиновые фенолы; S – сирингиловые фенолы; C – циннамиловые (коричные) фенолы.

В целом наибольшие значения отношения кислоты/альдегиды характерны для почв опытных площадок под елью сибирской (1.89) и сосной кедровой (1.13). Полученные результаты коррелируют (коэффициент корреляции r > –0.60) и с более низкими значениями pH в этих почвах. Характер трансформации лигнина неодинаков в почвах опытных площадок под покрытосеменными и голосеменными породами. Так, в почвах на площадках с осиной и березой, где отмечается невысокое содержание лигнина (1.57–2.65 мг/г С_орг) наблюдается минимальная метаболическая активность (VSC/N 20–38), величина параметра трансформации боковых цепочек лигнина по отношению к исходной растительной ткани составляет 27–28%, а степень окисления средняя – 0.72–0.76. В почвах под сосной сибирской (Pinus sibirica Du Tour) и сосной обыкновенной (Pinus sylvestris L.) степень трансформации боковых цепей лигнина составляет 46–51%, а степень окисления максимальная – 2.60. При этом ткани лиственницы демонстрируют наибольшую устойчивость к биодеградации и по минимальной степени окисленности биополимера (16), и по минимальной величине Tr = 17.4, хотя среди фенольных производных доминирует ванилиновый альдегид (1.73 мг/г С_орг). Травянистый тип трансформации лигнина также отличается низкими величинами параметра Tr (21%) и низкими значениями отношения кислоты/альдегиды (0.50) при равном соотношении ванилиновых, сирингиловых и циннамиловых структур.

Полученные данные свидетельствуют, что несколько десятилетий оказалось недостаточно, чтобы экосистема достигла состояния равновесия, характерного для естественных лесов бореальной зоны, так как природа гумуса почв вторичных лесов имеет черты травянистого генезиса. Об этом свидетельствуют и значения отношения V : S : C (2 : 1 : 1), которое более характерно для злаков (1 : 1 : 1), чем для естественных лесных фитоценозов (5 : 2 : 1; 7 : 9 : 1) [18].

Что касается верхних горизонтов ферралитных почв Амазонии, то, прежде всего, необходимо отметить, что максимальным накоплением лигнина отличаются почвы площадок под естественным лесом (9.56 ± 0.68 мг/г С_орг) и под плодовыми садами (9.07 ± 0.49 – 9.13 ± 0.33 мг/г С_орг) (табл. 4). В ферраллитных почвах Амазонии в условиях экосистем вторичных лесов содержание суммы лигниновых фенолов снижается до 1.60 мг/г С_орг. При этом все ферралитные почвы отличаются низкими величинами рН – сильно и очень сильно кислыми.

Таблица 4. Химические свойства и содержание лигнина в ферралитных почвах

Распределение лигнина по профилю ферралитных почв (краснозема) под естественным лесом обнаруживает максимальное накопление полимера в горизонте 5–10 см до 7.79 мг/г С_орг. Как показали исследования [17] ванилиновые фенолы, прежде всего, ванилиновая кислота, наиболее подвижны в условиях промывного водного режима и аккумулируются, как правило, в иллювиальной части профиля. В профиле ферралитных почв (краснозема) максимальное накопление ванилиновых фенолов обнаруживается на глубине 5–10 и 10–20 см (табл. 5), что отражается в отношение V : S : C, равном 18 : 17 : 1 (табл. 6). Верхний горизонт подстилки под естественным лесом Bacaba, отличает соотношение V : S : C, составляющее 2 : 2 : 1, более типичное для трав.

Таблица 5. Содержание продуктов окисления лигнина в ферралитных почах, мг/г С_орг

Примечание. Стандартное отклонение (Х ± SE), n = 3–4.

Таблица 6. Параметры биотрансформации лигнина в ферралитных почвах

Примечание. Стандартное отклонение (Х ± SE), n = 3–4.

Интересно заметить, что в структуре лигнина всех исследуемых почв тропической зоны обнаружено явное преобладание ванилиновых фенолов. Подобные факты описаны ранее для красноцветных почв можжевеловых лесов паркового типа Северного Тянь-Шаня [17], и, вероятно, являются характерной особенностью биохимического состава опада субтропической флоры. При этом в верхних горизонтах почвы, где значительный вклад в состав опада вносят корни трав, содержание ванилиновых единиц несколько уступает кумариловым и феруловым кислотам. В остальной части профиля главным источником ванилиновых фенолов выступают одревесневшие корни деревьев [15].

Композиционный состав лигнина в почвах остальных площадок вторичных тропических лесов лежит в интервале от 3 : 2 : 1 до 8 : 2 : 1, 6 : 5 : 1.

Наибольшую степень трансформации биополимера обнаруживают молекулы лигнина в поверхностных горизонтах почв под естественным лесом (31.35 %) и под плодовыми культурами (24.28–34.10%). При этом самых высоких величин параметр трансформации достигает в минеральных горизонтах краснозема (46–49%), что согласуется и с узким отношением лигнин/азот. И лишь в почвах плантации производственной древесины и в почвах под вторичным лесом степень трансформации биополимера минимальна (8.45–15.05%). При этом органическое вещество почвы площадки под вторичным лесом (Vismia) характеризуется и наибольшей обогащенностью метаболическим углеродом (VSC/N – 8.0).

Полученные результаты подтверждаются и величинами высокой степени окисленности биополимера (Ac/Al)v, достигающими максимальных значений в почвах под естественным лесом (1.36) и под плодовыми садами (1.03) при относительном постоянстве содержания общего лигнина (VSC, мг/г С_орг 9.1–9.6), что, видимо, обусловлено максимальной адаптацией микроорганизмов, разлагающих растительные остатки в естественном тропическом лесу.

Гумусовые горизонты ферраллитных почв на плантациях производственной древесины и, особенно, под вторичным лесом характеризуются средними значениями степени окисленности (Ac/ Al)v = 0.24–0.38, что значительно ниже подобных величин для серых почв бореальной зоны и для агрогенных почв Русской равнины [2]. Можно заключить, что углерод, накопленный в процессе фотосинтеза древесными растениями в зоне тропиков и субтропиков, не депонируется в почве даже в форме такого устойчивого к биотрансформации соединения, как лигнин.

ЗАКЛЮЧЕНИЕ

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

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

Биохимический состав тканей растений в различных экосистемах решающим образом влияет на процесс гумификации и определяет механизм образования гумуса. Несмотря на сложную предысторию объектов вторичного облесения, полученные результаты позволяют выделить разные типы биотрансформации растительных остатков в почвах различных экосистем. Тип биотрансформации лигнина в почвах под голосеменными породами характеризуется самым высоким содержанием суммы продуктов окисления лигнина, максимальной степенью трансформации биополимера и максимальными значениями отношения кислоты/альдегиды. Тропический тип биодеструкции лигнина отличает преобладание ванилиновых фенолов, высокое содержание лигнина в целом и обогащение его трудноразлагаемым структурным углеродом, средняя степень окисленности молекул. Последний факт существенно снижает возможности ферралитных почв выступать в качестве планетарного “депо” углерода, по сравнению с почвами бореальных лесов умеренного пояса. Тип биодеструкции лигнина покрытосеменных древесных (лиственных) пород в почвах березовых и осиновых экосистем отличается минимальной метаболической активностью биоты, средней степенью трансформации и средней степенью окисления биополимеров лигнина. Травяной тип трансформации лигниновых фенолов характеризуется средней степенью преобразованности биополимера и низкими значениями параметра кислоты/альдегиды в почвах из-за невысокого содержания ванилиновых структур и равномерного распределения продуктов V, S, C в их тканях (1 : 1 : 1).

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

БЛАГОДАРНОСТЬ

Авторы выражают глубокую благодарность профессору W. Zech, доктору I. Lobe, профессору РАН О.В. Меняйло за предоставленные образцы почв и аналитическую приборную базу.

СОБЛЮДЕНИЕ ЭТИЧЕСКИХ СТАНДАРТОВ

В данной работе отсутствуют исследования человека или животных.

ФИНАНСИРОВАНИЕ РАБОТЫ

Работа выполнена при поддержке Российского научного фонда № 23-24-00155.

КОНФЛИКТ ИНТЕРЕСОВ

Авторы заявляют, что у них нет конфликта интересов

About the authors

I. V. Kovalev

Author for correspondence.
Email: kovalevmsu@mail.ru
Russian Federation, Moscow, 119991

N. O. Kovaleva

Email: kovalevmsu@mail.ru
Russian Federation, Moscow, 119991

References

Supplementary files