CC BY
60
ВЛИЯНИЕ МЕЗОКЛИМАТА
и атмосферных промышленных загрязнений на радиальный прирост
сосны обыкновенной
A. С. ПОПОВ,
кандидат сельскохозяйственных наук, доцент,
B. В. ФОМИН (фото),
кандидат сельскохозяйственных наук, доцент, Ю. В. ШАЛАУМОВА, аспирант, Уральский ГЛТУ
620100, г. Екатеринбург, Сибирский тракт, д. 37
Ключевые слова: Pinus sylvestris L., сосна, радиальный прирост, мезоклимат, температура, промышленное загрязнение, Средний Урал.
Keywords: Pinus sylvestris L., radial increment, mesoclimate, temperature, industrial pollution, the Middle Urals.
Цель и методика исследований.
Реакция ростовых процессов древо-стоев на действие абиотических, биотических или антропогенных факторов — традиционный предмет изучения в дендроэкологии. Чаще всего исследуются отклики древесных растений на изменения климата в региональном и глобальном масштабах [1,
2, 3, 4, 5]. В естественных условиях древо-стои не испытывают на себе воздействие со стороны отдельно взятых климата или загрязнения. Эти и многие другие факторы влияют на рост древесных растений одновременно. В современной научной литературе практически отсутствуют попытки изучения взаимосвязи мезоклимата и промышленных загрязнений и их совместного воздействия на интенсивность ростовых процессов древесных растений. В ряде работ, посвященных связи динамики радиального прироста и концентрации в атмосферном воздухе различных поллютантов, зафиксированы противоречивые результаты. Авторы одних отмечают негативный характер влияния атмосферных выбросов на рост древостоев [5, 7, 8], в других приводятся свидетельства улучшения динамики прироста в градиенте концентрации поллютантов [9, 10]. Целью данной работы является изучение влияния мезоклимата и атмосферного промышленного загрязнения на Р°ст древесных растений, а также > ьяё. и /*
^ ■ н -ППе.;*... {
Нпвоуткннск
'о ' - '
А.
да-
исследование обеспеченности теплом и осадками территорий при различном содержании Э02 и взвешенных веществ в атмосферном воздухе.
Район исследований расположен в пределах Ревдинско-Первоуральского промышленного узла (г. Ревда, Свердловская область, южная тайга). Крупнейшим источником атмосферных промышленных выбросов является Среднеуральский медеплавильный завод (СУМЗ). Основные загрязнители — диоксид серы (Э02), тяжелые металлы (Си, РЬ, Zn, Cd и др.), взвешенные вещества (пыль). В качестве объектов исследований были выбраны сосновые древостои естественного происхождения 60-88-летнего возраста. Схема расположения объектов исследования и источника промышленного загрязнения представлена на рисунке 1. Были заложены три пробные площади, на каждой из которых находилось не менее 100 деревьев сосны обыкновенной. Первая контрольная пробная площадь (КПП 1) расположена в зоне с фоновым уровнем загрязнения атмосферы на удалении 19,1 км к юго-западу от СУМЗа, вторая (КПП 2) находится на удалении 57,5 км к западу от СУМЗа и в 10 км севернее г. Михайловска. Пробная площадь в импакт-ной зоне (ИПП) расположена на расстоянии
1,8 км к западу от СУМЗа.
На пробных площадях (ПП) проводили измерение диаметров всех деревьев сосны обыкновенной (Pinus Бу^евШв L.), высоту и возраст определяли только у отобранных модельных деревьев с диаметром, близким к среднему на ПП. Характеристики древостоев приведены в таблице 1. Лесорастительные условия на всех ПП сходные: низкогорные и предгорные (200-500 м над у. м.) участки, по режиму увлажнения — устойчиво свежие, расположенные на дренированных пологих придолинных склонах с маломощными дерново-подзолистыми суглинистыми
горно-лесными почвами на хорошо водопроницаемых осадочных породах, коренной тип леса — ельник-сосняк травяной [11].
На основании данных перечета деревьев были выбраны 18-19 модельных деревьев. Для каждого дерева получен керн возрастным буравом Пресслера с западной стороны ствола на высоте 0.3 м от шейки корня. Дендрохронологические исследования проводили по методикам, изложенным в работах С. Г. Шиятова и др. [12] и Н. В. Ловелиуса [13]. Измерения ширины годичных колец выполняли на измерительном комплексе LINTAB (Германия). Перекрестную датировку кольцевых хронологий проводили в программе COFECHA (университет Аризоны, США). Стандартизацию древесно-кольцевых хронологий осуществляли путем расчета индексов приростов с помощью программ TREND и ARSTAN [14]. Данная процедура позволяет максимально исключить влияние на динамику ростовых процессов возраста и трофических особенностей мест произрастания. Индексацию проводили путем сопоставления значений подобранной кривой со значениями изучаемой хронологии. Каждую древесно-кольцевую серию индексировали кривой с 60-летним сплайном и частотой отклика 50 %. Величину годичного индекса прироста для ПП находили вычислением среднего арифметического индексов всей совокупности древесно-кольцевых хронологий пробной площади с определением ошибки среднего.
Характеристику климата района исследований выполнили по среднесуточным данным инструментальных измерений температуры воздуха и выпавших осадков метеорологических станций Михайловск (287 м над у. м.) и Ревда (324 м над у. м.) [15]. Отсутствующие и недостоверные климатические данные за отдельные промежутки времени восстанавливали по данным метеостанции Екатеринбург (280
А. 1
МпхаКлпвск 1 |________ДО КМ______(
Рисунок 1
Схема расположения пробных площадей и Среднеуральского медеплавильного завода. КПП1 и КПП2 — контрольные пробные площади, ИПП — импактная пробная площадь
ммм. т-ачи. narod. ги
Таблица 1
Лесоводственно-таксационные характеристики древостоев контрольных (КПП1 и КПП2) и
импактной (ИПП) пробных площадей
№ ППП Расстояние от СУМЗа, км Средние
возраст, лет диаметр, см высота, м
КПП 1 19,1 77.6±1.2 (18) 36.5±0.7 (103) 27.0±0.1 (18)
КПП 2 57,5 61.1±1.6 (18) 30.7±0.9 (149) 26.9±0.4 (18)
ИПП 1,8 84.5±1.0 (19) 27.3±0.6 (151) 20.9±0.3 (19)
Примечание. Приведено среднее ± ошибка, в скобках - число учетных деревьев.
15
ОТ -
Т ъ
1=5 ж
1=6
Т
Кц :
Т с К*
'к,
м над у. м.) с помощью регрессионного анализа. Линейные регрессии строили для температуры воздуха по всему массиву данных ^2 > 95 %). Для данных по осадкам почти все линейные зависимости рассчитывали по месячным суммам величин ^2 > 50 %), исключение составили 5 значений, которые были восстановлены по линейным регрессиям только за теплый период (апрель-октябрь) ^2 > 45 %).
Обеспеченность древостоев доступной влагой оценивали гидротермическим коэффициентом Селянинова ^Т):
Таблица 2
Корреляция (ранговый коэффициент Спирмена) между дендрохронологическими сериями на контрольных (КПП 1 и КПП 2) и импактной (ИПП) пробных площадях, а также и гидротермическими коэффициентами, рассчитанными по данным метеостанций Ревда и
Михайловск за период 1962-2002 гг. (п=40)
Индексы КПП 2 Индексы ИПП ОТ Ревда ОТ Мих. PT5-6 Ревда PT7-8 Ревда PT5-6 Мих. PT7-8 Мих.
Индексы КПП 1 0,45 (0,019) 0,33 (0,036) 0,05 (0,754) 0,05 (0,754) 0,08 (0,624) -0,02 (0,894) 0,06 (0,724) 0,12 (0,458)
Индексы КПП 2 -0,01 (0,932) 0,15 (0,361) 0,27 (0,095) 0,24 (0,140) 0,04 (0,817) 0,37 (0,038) 0,04 (0,826)
Индексы ИПП -0,07 (0,672) -0,04 (0,807) 0,05 (0,781) -0,16 (0,314) -0,06 (0,707) 0,13 (0,427)
сдТТ
*=5
где Ri — количество осадков (мм) в иом месяце, а Т — сумма средних суточных температур выше 10”С в ном месяце.
Оценку условий увлажнения в первой и второй половине вегетационного периода производили с помощью плювиотер-мического коэффициента отдельно для мая-июня (РТ5-6) и июля-августа (РТ7-8):
Примечание. В скобках - достигнутый уровень значимости.
Таблица 3
Корреляция между дендрохронологическими сериями контрольных (КПП 1 и КПП 2) и импактной (ИПП) пробных площадях, гидротермическими коэффициентами, рассчитанными по данным метеостанций Ревда и Михайловск и средневзвешенными значениями концентрации (взвешенных веществ и диоксида серы) за период 1992-2002 гг.
где Ri — количество осадков (мм) в иом месяце, а й — сумма средних месячных температур в иом месяце. В средней полосе России данный коэффициент принято рассчитывать для периодов апрель-май, июнь-июль [16]. Сроки были нами изменены, поскольку в условиях Среднего Урала начало вегетационного периода сдвинуто по сравнению со средней полосой: среднесуточная температура 5С, необходимая для запуска ростовых процессов сосны обыкновенной [17], устанавливается в начале мая и стабильно держится до начала сентября.
На основании данных о максимальных и минимальных суточных температурах, зафиксированных метеостанциями Ревда и Михайловск в период 1992-2002 гг., анализировали величину суточных перепадов температур воздуха на каждой станции, после чего находили среднемесячное значение амплитуды температур.
Оценку загрязнения атмосферного воздуха проводили с помощью средневзвешенных значений концентраций пылевых частиц и диоксида серы по данным поста наблюдений ПНЗ № 1 ГУ «Свердловский ЦГМС-Р», расположенного в 9 км к северу от СУМЗа на территории г. Первоуральска:
Кц в, в, 5-6 п=10 Кц в, в, 5-8 п=10 Кц SO2 4-6 п=11 Кц SO2 5-6 п=11 Кц SO2 5-8 п=11
Индексы -0,16 -0,14 0,41 0,67 0,37
КПП 1 (0,624) (0,676) (0,195) (0,029) (0,239)
Индексы 0,50 0,60 -0,15 -0,18 0,06
КПП 2 (0,131) (0,072) (0,635) (0,574) (0,841)
Индексы 0,42 0,58 -0,52 -0,31 -0,15
ИПП (0,210) (0,085) (0,104) (0,327) (0,646)
ОТ 0,53 0,31 0,12 0,14 -0,20
Ревда (0,114) (0,354) (0,697) (0,666) (0,527)
ОТ 0,79 0,55 -0,05 0,137 -0,35
Мих, (0,017) (0,098) (0,885) (0,666) (0,275)
PT5-6 0,13 0,18 -0,19 -0,37 -0,38
Ревда (0,703) (0,598) (0,555) (0,249) (0,227)
PT7-8 0,52 0,30 0,20 0,13 -0,08
Ревда (0,122) (0,373) (0,536) (0,676) (0,796)
PT5-6 0,47 0,47 -0,22 -0,42 -0,36
Мих, (0,162) (0,162) (0,489) (0,180) (0,250)
PT7-8 0,61 0,37 0,11 0,25 -0,15
Мих, (0,066) (0,267) (0,740) (0,428) (0,646)
Примечание. В скобках - достигнутый уровень значимости.
корреляции между древесно-кольцевыми хронологиями с приближением к СУМЗу, и, соответственно, с усилением влияния промышленных выбросов на динамику радиального прироста сосновых древостоев. Значимая связь обнаружена для хронологий КПП 1 и КПП 2, а также КПП 1 и ИПП. Связь между хронологиями КПП 2 и ИПП, находящимися на максимальном удалении друг от друга, отсутствует.
Прямая корреляционная зависимость установлена между величинами индексов радиальных приростов деревьев КПП 2 и плювиотермическим коэффициентом мая-июня, рассчитанного по данным метеостанции Михайловск. Также существует
связь между древесно-кольцевой серией КПП 2 и гидротермическим коэффициентом Селянинова, рассчитанного за тот же период по данным метеостанции Михайловск. Обращает на себя внимание отсутствие корреляции между величинами приростов древостоев КПП 1, ИПП и гидротермическими коэффициентами.
Данные о наличии взаимосвязи концентраций поллютантов с величинами индексов древесно-кольцевых хронологий сосновых древостоев, произрастающих на КПП 1, КПП 2, ИПП, а также с показателями гидротермических коэффициентов, рассчитанными для метеостанций Ревда и Михайловск, приведены в таблице 3. В
где Кц — средневзвешенное значение концентрации поллютанта (мг/м3); С — средняя концентрация в иом месяце (мг/ м3); К — число дней в иом месяце.
Результаты исследований.
На рисунке 2 представлена динамика индексов радиальных приростов сосновых древостоев, сформировавшихся на различном удалении от СУМЗа, за период 1962-2002 гг. Данные таблицы 2 свидетельствуют о снижении величин ранговой
Рисунок 2
Сопоставление древесно-кольцевых хронологий сосновых древостоев, сформировавшихся в условиях с различным уровнем загрязнения. КПП1 и КПП2 — контрольные пробные площади,
ИПП — импактная пробная площадь. Вертикальные линии — ошибка среднего.
и
16
№№№. т-а¥и. па^. ги
течение всего периода вегетации наблюдается связь между GT, рассчитанным по данным метеостанции Михайловск, и средневзвешенными значениями концентрации взвешенных веществ.
Анализ корреляции динамики древеснокольцевых хронологий и средневзвешенных значений концентрации взвешенных веществ в воздухе позволил предположить наличие влияния запыленности атмосферы в мае-августе на рост древостоев КПП 2 и ИПП. Увеличение значений средневзвешенного показателя концентрации взвешенных веществ вызывает стимуляцию приростов сосновых древостоев, расположенных в зоне сильного загрязнения и на максимальном удалении от СУМЗа.
Установлена значимая корреляция между индексами древесно-кольцевых хронологий КПП 1 и средневзвешенными значениями концентрации диоксида серы в воздухе в мае-июне. Рост концентрации Э02 в воздухе в период начала роста сопровождается достоверным увеличением радиального прироста на КПП 1. Концентрация в воздухе импактной ИПП Э02, рассчитанная для апреля-июня, вероятно, негативно влияет на ростовые возможности древостоя. Обнаруженная зависимость не является значимой, но указывает на то, что повышение концентрации диоксида серы в воздухе в начале периода вегетации может вызывать снижение радиальных приростов у деревьев, находящихся в зоне сильного загрязнения.
Индексы прироста КПП 2 демонстрируют высокую степень корреляции с РТ5-6 и GT, рассчитанными по данным метеостанции Михайловск. Проведен анализ климатических параметров, зависящих от изменения концентрации поллютантов в воздухе. На рисунке 3 представлен график разности среднемесячных значений амплитуд температур, зафиксированных на метеостанциях Михайловск и Ревда в 1992-2002 гг. Средняя месячная разность между максимальными и минимальными температурами, зафиксированными на метеостанции Михайловск, была больше аналогичного показателя, рассчитанного для метеостанции Ревда, почти на всем временном отрезке с 1992 по 2002 гг. Можно предположить, что данный факт объясняется более высокой концентрацией взвешенных веществ и газов в атмосфере в районе действия СУМЗа, что может быть причиной появления «эффекта покрывала». Представленные результаты во многом совпадают с данными, которые приводит по Москве и ее окрестностям Л. Т. Матвеев [18].
Наблюдается связь между концентрацией в воздухе взвешенных веществ и суммой осадков вегетационного периода. На рисунке 4 представлены сравнительная динамика концентрации взвешенных веществ в воздухе и суммы выпавших осадков в мае-августе по данным метеостанций Ревда и Михайловск за период с 1992 по 2002 годы. Динамика выпадения осадков в период вегетации в районе этих станций в значительной мере совпадает. Увеличение концентрации взвешенных веществ в течение периода вегетации приводит к росту количества выпадающих за этот период осадков и в Михайловске ^эр = 0,82; п =
шшш.т-ауи. narod. ги
^ Ч U1 кП (О (О I-- оо СО 05 О}
05 05 05 05 05 05 05 05 05 05 05 05 05
05 05 05 05 05 05 05 05 05 05 05 05 05
Дата
■ Ревда •
-Михайловск
Рисунок 3
Сравнение средних месячных значений амплитуд температур на метеостанциях Ревда и
Михайловск в 1992-2002 гг.
0,4
0,3
0,2
□ .1
Год
Сумма осадков для Ревды — Сумма осадков для Михайловского -- - Концентрация взвешенных веществ
Рисунок 4
Графики хода концентрации взвешенных веществ в воздухе и суммы выпавших осадков за май-август по данным метеостанций Ревда и Михайловск в 1992-2002 гг.
11; p = 0,001), и в Ревде (Rsp = 0,22; n = 11; p = 0,113).
Полученные результаты указывают на ослабление сходства древесно-кольцевых хронологий сосновых древостоев по мере приближения к источнику загрязнения, максимальное значение коэффициента корреляции было получено при сопоставлении индексов прироста контрольных пробных площадей. Динамика индексов радиального прироста деревьев сосны КПП 2 показывает тесную связь с плювиотермическим коэффициентом мая-июня, рассчитанным по данным метеостанции Михайловск. Наиболее интенсивный рост деревьев сосны наблюдается в первой половине вегетационного периода [17]. Наличие взаимосвязи между древесно-кольцевыми хронологиями максимально удаленной от источника загрязнения КПП 2 с гидротермическими коэффициентами подтверждает тот факт, что обеспеченность растений доступной влагой в период вегетации является одним из важных условий нормального роста. Повышение величины плювиотерми-ческого коэффициента мая-июня, рассчитанного по данным станции Михайловск, указывает на улучшение обеспечения
древостоев доступной влагой и вызывает стимулирование ростовых процессов.
В свою очередь, для GT, рассчитанного по данным метеостанции Михайловск, характерна достоверная положительная корреляция со средневзвешенными показателями концентрации взвешенных веществ, рассчитанными для начала периода вегетации. Возможно, увеличение концентрации взвешенных веществ в условиях роста уровня выпадающих осадков связано с переносом пылевых частиц осадками в приземный слой атмосферы, что подтверждается графиком хода концентрации взвешенных веществ в воздухе и суммы выпавших осадков за май-август по данным метеостанции Михайловск (рис. 4).
Корреляции, рассчитанные для древесно-кольцевых хронологий КПП 2 и средневзвешенных значений концентрации взвешенных веществ, характерных для начала периода вегетации, свидетельствуют о возможной прямой связи между этими показателями. Можно предположить, что по уровню концентрации частиц пыли в приземном слое воздуха в начале периода вегетации можно судить о количестве выпадающих в это время на значительном
17
удалении (50-60 км) от источника загрязнения осадков и, соответственно, об обеспеченности древесных растений влагой. Количество выпадающих осадков влияет на интенсивность роста древостоев. В этом случае прямая зависимость характера древесно-кольцевых хронологий КПП 2 и средневзвешенных значений концентрации взвешенных веществ не обязательно является примером стимулирования роста древостоев в условиях фонового уровня загрязнения. Возможно, рост и того, и другого показателя наблюдается в связи с увеличением уровня выпадающих осадков.
Результаты сопоставление динамики индексов радиального прироста КПП 1 и средневзвешенных значений концентраций
диоксида серы, рассчитанных для мая-июня, указывают на возможность стимулирования роста деревьев сосны в результате действия на них малых концентраций SO2, что является подтверждением фактов, приведенных в ряде литературных источников [Смит 1985; Innes and Oleksyn 2000]. Значимых корреляций, подтверждающих негативное влияние диоксида серы на рост деревьев, расположенных в зоне сильного загрязнения, получено не было, однако полученные результаты указывают на возможность такого эффекта.
Динамика индексов радиального прироста древостоев КПП 1 и ИПП, расположенных рядом с источником выбросов, не несет в себе четкого отражения климатического
сигнала, на который в данных условиях накладывает отпечаток влияние атмосферного загрязнения. Атмосферный воздух в районе действия СУМЗа характеризуется высоким уровнем содержания в нем взвешенных веществ, что может являться причиной изменения теплового режима территории (рис. 3), а также нарушения режима солнечного освещения.
Работа выполнена благодаря финансовой поддержке РФФИ (грант 09-04-01004). Авторы выражают благодарность Е. В. Воробейчику и Р. М. Хантемирову за ценные замечания и предложения, которые были высказаны ими в ходе обработки и анализа данных, а также обсуждения статьи.
Литература
1. 1IPCC. Climate Change 2001: The Scientific Basis. Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge : Cambridge University Press; New York, 2001. P. 881.
2. Шиятов С. Г., Терентьев М. М., Фомин В. В. Пространственно-временная динамика лесотундровых сообществ на Урале // Экология. 2005. № 2. С. 69-75.
3. Шиятов С. Г., Терентьев М. М., Фомин В. В., Циммерманн Н. Е. Вертикальный и горизонтальный сдвиги верхней границы редколесий и сомкнутых лесов в XX столетии на Урале // Экология. 2007. № 4. С. 243-248.
4. Капралов Д. С., Шиятов С. Г., Моисеев П. А., Фомин В. В. Изменения в составе, структуре и высотном положении мелколесий на верхнем пределе их произрастания в горах Северного Урала // Экология. 2006. № 6. С. 403-409.
5. Ivshin A., Shiyatov S. The assessment of subtundra forests degradation by dendrochronological methods in the Norilsk industrial area // Dendrochronologia. 1995. V. 13. P. 113-126.
6. Фомин В. В., Шавнин С. А. Влияние горного рельефа и аэропромышленных загрязнений на биометрические характеристики сосновых древостоев // Экология. 2002. № 3. С. 170-174.
7. Schulze E. D., Beck E., Muller-Hohenstein K. Plant ecology // Springer. Berlin-Heidelberg, 2005. P. 440.
8. Augustaitis A. A New Imitative Model to Predict the Impact of Air Pollutants on Scots Pine Health and Radial Increment // Water, Air, Soil Pollut. 1999. V. 116. P. 437-442.
9. Смит У. Лес и атмосфера: Взаимодействие между лесными экосистемами и примесями атмосферного воздуха // Прогресс. 1985. 432 с.
10. Innes J., Oleksyn J. Forest Dynamics in Heavily Polluted Regions // Oxford : CABI Publishing, 2000. P. 248.
11. Колесников Б. П., Зубарева Р. С., Смолоногов Е. П. Лесорастительные условия и типы лесов Свердловской области: Практическое руководство // Свердловск : УНЦ АН СССР, 1973. 178 с.
12. Шиятов С. Г., Ваганов Е. А., Кирдянов А. В., Круглов В. Б., Мазепа В. С., Наурзбаев М. М., Хантемиров Р. М. Методы дендрохронологии. Часть I. Основы дендрохронологии. Сбор и получение древесно-кольцевой информации // Красноярск, 2000. 79 с.
13. Ловелиус Н. В. Изменчивость прироста деревьев // Наука. 1979. 230 с.
14. Cook E. A time series analysis approach to tree-ring standardization // Tucson : University of Arizona, 1985. P. 171.
15. Метеорологический ежемесячник // Новосибирск. 1962-2002. Ч. 2.
16. Педь Д. А. О показателе засух и избыточного увлажнения // Труды гидрометеоцентра СССР. Вып. 156. 1975. С. 19-38.
17. Лебеденко Л. А. Динамика размножения камбиальных клеток у сосны и ели // Восстановление леса на северо-западе РСФСР : сб. тр. Лен. НИИЛХ, 1978. С. 101-111.
18. Матвеев Л. Т. Влияние большого города на метеорологический режим // Известия РАН. Серия географическая. 2007. № 4. С.97-102.
экологическая оценка территории в зоне действия атмосферных промышленных загрязнений медеплавильного производства
в. в. ФОМИН,
кандидат сельскохозяйственных наук, доцент, А. А. НИКОЛАЕВ, аспирант, Уральский ГЛТУ
620100, г. Екатеринбург, Сибирский тракт, д. 37
Ключевые слова: состояние древостоев, экологическое зонирование, климат, аэропромышленные загрязнения, цифровая модель рельефа, Средний Урал.
Keywords: tree stands state, ecological zoning, climate, air pollution, digital elevation
Цель и методика исследований.
Древесная растительность, произрастающая в условиях горно-увалистого рельефа местности на территории, прилегающей к крупному промышленному узлу (г. Ревда, Свердловская область), находится под воздействием разнообразных экологических
18
факторов, уровни действия которых изменяются в пространстве и времени. Цель работы — комплексная экологическая оценка территории в зоне действия атмосферных промышленных загрязнений на основе геоинформационной модели пространства.
model, the Middle Urals.
В ходе выполнения международного проекта INTAS 93-1645 сотрудниками Уральского государственного лесотехнического университета и Института экологии растений и животных УрО РАН были собраны, обработаны и проанализированы данные о состоянии сосновых древостоев
www. m-avu. narod. ru
