Геоморфометрия в науках о Земле и экологии, обзор методов и приложений

БИОЛОГИЯ И ЭКОЛОГИЯ

УДК 551.4

ГЕОМОРФОМЕТРИЯ В НАУКАХ О ЗЕМЛЕ И ЭКОЛОГИИ, ОБЗОР МЕТОДОВ И ПРИЛОЖЕНИЙ

© 2006 П.А. Шарый

Институт физико-химических и биологических проблем почвоведения РАН, г. Пущино

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

Введение

При всей специфике и разнообразии методов конкретных научных дисциплин, изучающих ландшафты, природные экосистемы и строение Земли, геофизические факторы носят общезначимый характер для этих дисциплин. Одно из важнейших мест среди них занимает рельеф. Причина этому состоит не только в том, что рельеф является перераспределителем влаги, тепла, света и литодина-мических потоков, но еще и в том, что он сравнительно легко доступен измерению. Кроме того, если касаться почв, то все другие факторы почвообразования (климат, геология, время, организмы) по отношению к рельефу являются факторами рельефообразова-ния [8], то есть в той или мере находят отражение в данных о рельефе.

В науках о Земле рельеф не изучается только как поверхность Земли. Действительно, уже попытка рассмотреть направление течения воды по земной поверхности требует привлечения другого геофизического фактора - гравитационного поля Земли, - без которого это течение немыслимо. Рассмотрение связанной с рельефом пространственной изменчивости терморежима склонов неизбежно требует привлечения еще одного геофизического фактора - векторного поля солнечного излучения. Обобщая эти факты, можно заметить, что при рассмотрении роли и аспектов действия рельефа анализируется

физически и математически двойная система "земная поверхность + векторное поле" [63], причем из множества геофизических векторных полей (магнитное, электрическое,..., излучение от каждой звезды) именно гравитация и солнечная радиация оказываются наиболее важными для изучения пространственной изменчивости ландшафтов.

При том, что характеристики многих компонент ландшафта (ресурсы воды, почв, флоры и фауны; верхняя часть земной коры и др.) оказываются тесно связаны с рельефом, возникает естественный вопрос: достаточно ли обычно используемых двух-трех характеристик рельефа (высоты, крутизны и экспозиции склонов) для описания всего многообразия форм рельефа и всех аспектов его влияния на ландшафты и экосистемы? - Анализируя этот вопрос, И.С. Эванс [32] отметил "очевидность синдрома "философского камня", подразумевающего веру, что существо природы форм поверхности может быть выражено одним параметром". И двух-трех параметров для этого явно недостаточно. В то же время, неколичественные методы описания рельефа (деление рельефа на "волнистый", "слабоволнистый", "бугристый" и т.п.) уже исчерпали себя, не будучи в состоянии ни дать новые содержательные ответы на вопросы о видах механизмов влияния рельефа на перераспределение влаги, тепла и т.д. (например, о механизмах связи памяти в ланд-

шафтах с рельефом), ни ввести объективные классификации форм рельефа, ни дать методы обнаружения статистически достоверных связей между характеристиками рельефа и изучаемых компонент ландшафта. О большом интересе к количественным методам анализа рельефа, составляющим содержание геоморфометрии как науки о рельефе свидетельствует, например, библиография Р. Пайка [58] по геоморфометрии и ее применениям, насчитывающая тысячи ссылок.

Несмотря на это обилие источников, в подавляющем большинстве посвященных использованию геоморфометрии (а не ее развитию), число базовых, то есть не составных, количественных характеристик рельефа в экологии, почвоведении и геоботанике оставалось, по сути, ограничено шестью - к перечисленным выше трем добавлялись еще горизонтальная кривизна, вертикальная кривизна и площадь сбора. Классификации форм рельефа продолжали носить субъективный характер, опираясь на вещественное содержание этих форм - "интрузивные" и другие формы [17], - либо на субъективное деление форм рельефа, например, "по пожеланиям агробизнес-партнеров" [43], иногда с малоосмысленными статистическими процедурами усреднения [11], но не геометрические принципы, осознанные в классических работах как наиболее адекватные задачам этого рода [9]. Редкие исключения [5; 73] не создавали целостной картины, либо развивались в физике [16] и математике [3], не находя применений в исследованиях ландшафтов и экосистем.

В этих условиях важным моментом для приложений геоморфометрии в науках о Земле стало расширение набора базовых количественных характеристик рельефа в самой геоморфометрии, инициированное П.А.Шарым [63], а также развитие на основе расширенной системы этих характеристик классификаций форм рельефа, построенных по геометрическому принципу. Эти задачи решались как на базе интеграции результатов из ряда научных дисциплин (от физики и математики до физической географии и геоморфологии), так и введением

нами новых характеристик рельефа [63; 65; 67], и далее - введением и изучением свойств основанных на них новых классификаций форм рельефа [19; 20; 67]. Разрабатывались также и новые методы применения расширенной системы характеристик рельефа для анализа связей рельеф-почва-растения в ландшафтах [24; 55; 66]; для расчета и построения карт-матриц характеристик рельефа, а также для проведения статистических сравнений характеристик ландшафта с ними автором статьи разработана программа "Аналитическая ГИС Эко" [64].

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

Понятия и методы геоморфометрии

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

Учитывая современное развитие математической теории поверхностей, то есть дифференциальной геометрии поверхностей [10] и топологии [1], кажется, на первый взгляд, естественным просто заимствовать достижения этих наук, не занимаясь разработкой методов геоморфометрии. Однако давно используемые в науках о Земле крутизна и ориентация (экспозиция) склонов [29; 48] отсутствуют в дифференциальной геометрии поверхностей и топологии, поскольку в них ана-

лизируется поверхность, а не математически двойной объект "земная поверхность + векторное поле". Первые важные понятия, расширяющие геоморфометрию как науку, возникли в почвоведении, где наиболее явны и существенны результаты влияния водных и литодинамических потоков и терморежима склонов на формирование и свойства почв: А. Аандал [25] - при изучении азота в почвах; Ф. Трёх [73] - при анализе и классификации пространственной изменчивости гидромор-физма почв. Известны и другие ранние попытки развития недостающих в физике и математике подходов, например, П.К. Соболевский [13] в геологии был очень близок к введению горизонтальной кривизны для разделения выпуклых и вогнутых отрогов холмов.

МВ могут характеризовать либо систему "земная поверхность + векторное поле", либо саму поверхность (без учета геофизических полей), когда мы говорим о геометрических формах рельефа, то есть все МВ можно разделить на поле-специфичные и игнорирующие поля [63]. Далее, МВ в данной точке может описывать либо свойства ограниченной окрестности этой точки (локальные МВ), либо же для ее определения в точке требуется рассмотрение протяженных участков рельефа (нелокальные МВ). В последнем случае имеется два варианта. В первом рассматривается ограниченный участок планеты, и тогда нелокальная величина описывает в данной точке относительное положение этой точки в рельефе [70], например, площадь сбора описывает положение данного участка в бассейне. Во втором варианте рассматривается вся замкнутая поверхность планеты; тогда такие понятия, как бассейн, оказывается невозможно определить, но зато появляются новые возможности, например, определить соотношение числа вершин холмов, доньев депрессий и планетарных сед-ловых точек [48]. МВ, отвечающие первому варианту, были названы региональными МВ,

морфометрические величины, отвечающие второму варианту - глобальными или планетарными МВ [65]. Это дает возможность ввести непересекающиеся классы МВ и форм рельефа (табл. 1).

Дополнительно к этому, имеется класс В' глобальных поле-специфичных МВ и класс D' глобальных игнорирующих поля МВ [65].

Морфометрические величины классов А и В описывают поверхностный сток, терморежим склонов и т.п., в то время как МВ классов С и D описывают геометрические формы. При этом МВ класса А описывают, например, поверхностный сток локально (по свойствам малой окрестности точки), в то время как МВ класса В описывают его регионально, то есть по относительному положению точки в рельефе. Аналогично, МВ класса С описывают геометрические формы локально, а МВ класса D - регионально, по относительному положению.

Важная систематизация представлений о поверхностном стоке в геоморфомет-рии состоит в установлении связей между механизмами аккумуляции веществ, перемещающихся под действием силы тяжести по земной поверхности, и некоторыми формами рельефа, такими как отроги, террасы, притальвеговые и приводораздельные области. Имеется два механизма аккумуляции (основных), а также сноса, таких веществ: за счет их сближения в плане и за счет их относительного замедления в профиле склона (рис. 1).

В 1948 году А. Аандал [25] предположил, что 1-ый механизм аккумуляции описывается горизонтальной кривизной, а 2-ой - вертикальной кривизной земной поверхности (тогда еще не были найдены формулы для этих кривизн). П.А. Шарый доказал теоремы [63] о том, что при локальном описании рельефа кривизнами горизонтальная кривизна описывает первый механизм аккумуляции, а вертикальная кривизна - второй. Это явилось

Таблица 1. Классификация морфометрических величин и понятий [65]

Морфометрические величины Локальные Региональные

Поле-специфичные Класс А Класс В

Игнорирующие поля Класс С Класс Б

Рис. 1. Слева - 1-ый механизм аккумуляции (за счет сближения линий тока в плане, эти перпендикулярные горизонталям линии показаны стрелками), изображена топографическая карта. Справа - 2-ой механизм аккумуляции (за счет относительного замедления потоков в профиле), показан разрез вдоль профиля склона. Условные обозначения: 1 - области дивергенции, 2 - области конвергенции, 3 - горизонтали, 4 - грунт, 5 - воздух

основанием для аргументированного использования этих МВ как описывающих механизмы аккумуляции.

Аандал [25] определил области с отрицательной горизонтальной кривизной кЪ (области конвергенции, рис.1) как вогнутые отроги, а области с положительной кЪ (области дивергенции) как как выпуклые отроги, области с отрицательной вертикальной кривизной ку как вогнутые террасы, а области с положительной ку как как выпуклые террасы.

На рис. 2 показаны карты кЪ и ку для одной и той же местности, отвечающие ее отрогам и террасам.

Наиболее популярными алгоритмами расчета локальных МВ по матрице высот земной поверхности являются метод Эванса-Янг [32; 72], метод Шарого [63] и метод Зевен-бергена-Торне [79]; сравнение их дано в работе [61]. Для расчета кривизн предпочтителен разработанный нами модифицированный метод Эванса-Янг [65], так как он существенно меньше выделяет артефактные направления на картах кривизн при расчете по матрице высоты [61].

При региональном описании 1-го механизма аккумуляции, т.е. при описании его по относительному положению точки в рельефе,

используется площадь сбора MCA, рис. 3.

Площадь сбора показывает всю гидросеть местности - как реализованную, то есть реки и водоемы, так и потенциальную, то есть сухие долины и временные водотоки [67]. Алгоритм ее расчета опубликован в работах [46; 35].

Основные аспекты действия рельефа можно систематизировать следующим образом [67]:

(1) поверхностный сток сыпучих и жидких веществ,

(2) расчлененность (изрезанность) местности,

(3) геометрические формы (связанные памятью в ландшафтах),

(4) терморежим склонов,

(5) высотная поясность.

Помимо этих основных аспектов, действие рельефа проявляется еще в ряде более частных, специфических для конкретной местности, явлений (соответственно, менее общего характера), таких как описанная для бассейна р. Меконг во Вьетнаме [37] резкая зависимость свойств почв от высоты при перепаде высот около 1,2 м на участке размером ~50кмх100км, связанная с тем, что высота над уровнем моря на этой равнинной

Рис. 2. Карты отрогов (к^ слева) и террас (ку, справа) Аандала для одной и той же местности, участка Данки, рассчитанные по матрице высот с шагом решетки 10 м. Эти формы рельефа отвечают двум механизмам аккумуляции, описанным на локальном уровне. Вогнутые отроги и террасы показаны темным, выпуклые - светлым

0 500 1000 1500 м

Рис. 3. Карта площади сбора для той же местности, что на рис. 2. Черные участки показывают реки

прибрежной территории есть одновременно и глубина до грунтовых вод. С.С. Неуст-руев [9] отмечал связь микроклимата некоторых местностей с рельефом, особенно заметную в горных странах, где в долинах наблюдается более контрастный микроклимат - зимой там застаивается холодный воздух с гор, а летом они сильно прогреваются. Описанный В.М. Фридландом [18] для Кавказа ряд частных примеров связи потоков воздушных масс с рельефом, таких как отсутствие лесов из-за открытости участка местности северным ветрам по условиям рельефа или изменения увлажнения почв из-

за турбулентности ветров в предгорной равнине, также следует относить к неосновным аспектам действия рельефа. Проблема связи микроклимата в рельефом изучена пока явно недостаточно [31].

Система морфометрических величин

Система базовых МВ представлена в табл. 2. Отметим, что ориентация (экспозиция) склонов, крутизна и освещенность склонов зависимы между собой [67], т.к. эти три МВ зависят от двух частных производных; Н.М.Сибирцев [12] отмечал, что терморежим склонов характеризуется перпендикулярностью падения солнечных лучей на земную поверхность, а потому зависит как от ориентации, так и от крутизны склонов. Освещенность склонов непосредственно описывает эту перпендикулярность [67], а потому значительно лучше описывает терморежим склонов, чем их ориентация.

Физический смысл и формулы каждой из приведенных в табл. 2 МВ описаны нами в работах [65; 67], здесь мы не останавливаемся на этом.

К перечисленным здесь МВ следует добавить глубины бессточных депрессий и высоты В-холмов (буква В здесь означает класс); под последними понимаются холмы, возвышенности, горные системы и т.п. Вне депрессий (холмов) эти МВ полагаются равными нулю [65].

Таблица 2. Система базовых морфометрических величин

№ Название МВ (обозначение) Что описывает Источник Класс

МВ, ОПИСЫВАЮЩИЕ ПОВЕРХНОСТНЫЙ С ГОК

0 Ориентация склонов (A0) Направление потоков -- А

1 Крутизна склонов (GA) Скорость потоков -- А

2 Горизонтальная кривизна (kh) Й.Крхо [40] А

3 Площадь сбора (MCA) 1-ый механизм аккумуляции Дж.Спейт [69] В

4 Дисперсивная площадь (MDA) Дж.Спейт [70] В

5 Вертикальная кривизна (kv) 2-ой механизм аккумуляции И.С.Эванс [32] А

6 Разностная кривизна (E) Сравнивает 1-ый и 2-ой механизмы аккумуляции П.А.Шарый [63] А

7 Полная аккумуляционная Выявляет зоны относительной П.А.Шарый А

кривизна (KA) аккумуляции [63]

МВ, ОПИСЫВАЮЩИЕ РАСЧЛЕНЕННОСТЬ РЕЛ ЬЕФА

8 Полная кольцевая кривизна Извилистость потоков П.А.Шарый А

(KR) (расчлененность рельефа) [63]

9 Горизонтальная избыточная Эти две МВ расщепляют П.А.Шарый А

кривизна (khe) извилистость потоков на две [63]

10 Вертикальная избыточная кривизна (kve) компоненты П.А.Шарый [63]

11 Ротор (rot) Направление поворота линий тока П.А.Шарый [22] А

МВ, ОПИСЫВАЮЩИЕ ГЕОМЕТРИЧЕСКИЕ ФОРМЫ РЕЛЬЕФА

12 Несферичность (М) Отличие формы рельефа в точке от сферической П.А.Шарый [63] С

13 Средняя кривизна (Н) «Равновесную» (минимальную) поверхность, имеющую минимальную при данной границе площадь. Т.Юнг [77] С

14 Максимальная кривизна (kmax) Гребневые формы рельефа -- С

15 Минимальная кривизна (ктгп) Килевые формы рельефа -- С

16 Полная гауссова кривизна (К) Не меняется при изгибании Поверхности К.Ф.Гаусс [3] С

МВ, ОПИСЫВАЮЩИЕ ТЕРМОРЕЖИМ СКЛОНОВ

17 Освещенность склонов (F) Интенсивность прямого солнечного излучения -- А

МВ, ОПИСЫВАЮЩИЕ ВЫСОТНУЮ ЗОНАЛЬНОСТЬ

18 Высота земной поверхности (z) Изменение температуры и давления атмосферы -- А

Классификации форм рельефа

Первая классификация форм поверхности предложена К.Ф. Гауссом [3], она основана на знаках полной гауссвой кривизны и средней кривизны - МВ, описывающих геометрические формы. В ней 4 типа форм. Позднее Ф.Р. Трёх [73] предложил иную классификацию форм рельефа, основанную на знаках горизонтальной и вертикальной кривизн, отвечавшую двум основным механизмам аккумуляции. В ней также 4 типа форм. П.А. Шарый [63] предложил более общую

классификацию форм рельефа, включавшую в себя классификации Гаусса и Треха как частные случаи. В этой классификации 12 типов форм. Свойства этих классификаций, а также соотношения между ними изучены нами в работах [19; 20; 63; 67]. Классификации Гаусса и Трёха показаны на рис. 4.

Нами доказано [63], что каждый из этих типов форм является составным, т.е. составлен несколькими типами форм рельефа классификации Шарого. Типы форм последней показаны на рис. 5.

Рис. 4. Слева - типы форм рельефа классификации Гаусса, справа - типы форм рельефа классификации Трёха

Рис. 5. Двенадцать типов форм рельефа классификации Шарого. См. подробнее работу [20]

Помимо этих (основных) типов форм рельефа, можно определить еще редкие типы форм, образующих на картах линии или точки, обычно общие границы соседних основных типов форм [63]. Примеры редких типов форм показаны на рис. 6.

При расчете по матрицам высот вероятность встретить такие редкие формы чрезвычайно мала, то есть в частоте встречаемости редких и основных форм рельефа имеется природная иерархия [63]. Это проверено нами экспериментально [19]. Во многих работах, начиная с работы Ю.К. Ефремова [5], нет этого важного различения основных и редких типов форм.

Большее число типов форм в классификации Шарого, по сравнению с классификациями Гаусса и Трёха, означает и большую гиб-

Рис. 6. Примеры редких типов форм, показанных на плоскости, по осям которой отложены горизонтальная кривизна kh и вертикальная кривизна kv: I - рассеивающие с идеально прямым профилем

кость их использования. Так, В. Буивидайте и Г. Мозгерис [28] нашли при сравнении типологических единиц почв с основными типами форм этих трех классификаций наиболее тесную связь для классификации Шарого.

Описанные выше классификации форм рельефа относятся к локальным, то есть построенным с помощью локальных МВ. В силу негладкой природы земной поверхности [44] такие типы форм существенно зависят от масштаба карты [67], что создает заметные трудности при их использовании.

Эти трудности, однако, в гораздо мень-

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

Проблема негладкости земной поверхности

Негладкость (то есть недифференцируемость) земной поверхности проявляется в том, что за попытку рассчитать несуществующие производные методом конечных разностей мы "платим" неустранимой зависимостью полученных дифференцированием высоты по плановым координатам локальных МВ от масштаба карты или шага решетки матрицы [33; 67]. Хотя в России до наших исследований использовались только такие МВ, мнение о том, что все МВ локальны [7; 11] следует признать ошибочным. Негладкость земной поверхности приводит к целому ряду важных следствий, изученных нами в работе [67], среди которых есть и то, что основанные на использовании производных выше второй степени МВ вообще не дают содержательных карт, вопреки известной в стране "парадигме Ласточкина" [7], где они фигурируют, но либо просто ошибочны (например, водораздел невозможно выразить через производные, т.к. граница бассейнов -понятие региональное), либо никогда не проверялись на практике расчетами и построением карт для реальных местностей.

В то же время, как экспериментально и теоретически показано в наших исследованиях [65; 67], для ряда региональных МВ проблема зависимости МВ от масштаба имеет решения. То есть, для таких МВ существуют предельные значения при укрупнении масштаба, а значит, здесь применимо поня-

тие точности в значениях МВ, в том же смысле, что и в геодезии для высоты.

Использование геоморфометрии в науках о Земле и экологии

Разумеется, рельеф не является единственным фактором, от которого зависят характеристики ландшафта, такие как свойства почв или растительности. Уже по этой причине выявление связи характеристик ландшафта с рельефом носит характер поиска статистических связей (или их отсутствия), с обязательным установлением статистической достоверности найденных связей. В неколичественных или полуколичественных исследованиях эта достоверность часто не устанавливается, и результаты носят скорее характер непроверенных гипотез, нередко основанных на одностороннем анализе рельефа, например, с помощью только одной МВ.

Так например, И.Н. Степанов [14] использовал знак горизонтальной кривизны для построения почвенных карт любого масштаба в своем "методе пластики рельефа", неявно предположив тем самым, что первый механизм аккумуляции во всех случаях доминирует над вторым; статистическая достоверность этого им установлена не была. Этот подход подвергся справедливой критике со стороны других почвоведов [4]. Е.А. Гарши-нев и Б.П. Петрюк [2] использовали вертикальную кривизну для различения эрозионных и эрозионно-аккумулятивных земель. Ни роль этой характеристики рельефа по отношению к процессам в почвах и ландшафте, ни статистическая достоверность связи отмеченного авторами различия земель с этой характеристикой рельефа ими не устанавливались. Такая гипотеза противоречит прямым наблюдениям на пашне с помощью естественного радиоактивного цезия-137 растровой картины пространственной изменчивости чистой эрозии [42]; как видно из приведенного там "трехмерного" рисунка, водная эрозия связана как с первым, так и со вторым механизмами аккумуляции. Примеры этого рода можно продолжить; сосредоточимся, однако, на более успешных результатах.

Г. Синаи с соавт. [68] нашли высокий

коэффициент корреляции (г=+0,90) между лапласианом высоты DDZ и влажностью почв в климатических условиях Израиля. Долгое время причины оставались непонят-нымы, хотя результат явно привлек внимание исследователей. В работе [50] не было найдено значимой корреляции между крутизной склонов и продуктивностью пшеницы, но она оказалась значимой между AZ и продуктивностью пшеницы и кукурузы [71]. Затем нами было установлено, что AZ -2Н , где Н есть средняя кривизна [65]. Поскольку Н=1/2(ку+кЪ) [63], то причина состоит в том, что Н с равными весами описывает оба механизма аккумуляции, оказываясь отрицательной в зонах относительной аккумуляции и положительной - в зонах относительного сноса Трёха (рис. 4, справа). Значение г=+0,90 показывает, что влага в опыте Синаи с соавт. [68] была увеличена в зонах относительной аккумуляции. Эти авторы показали также сильное увеличение продуктивности пшеницы в этих зонах. Опыты по изучению связи влажности почв с микрорельефом были продолжены в наших исследованиях [6]. В гу-мидных климатических условиях юга Московской области на хорошо выраженном рельефе (овраг) было показано при шаге решетки w=3м, что коэффициент корреляции г = - 0,89 для связи относительной влажности серых лесных сильно смытых почв со средней кривизной. Почва оказывалась более влажной в зонах относительной аккумуляции. При этом коэффициенты корреляции влажности с ку и кЪ (взятыми порознь) в тех же точках оказались меньше, чем с Н.

Л. Ланион и Г. Холл [41] рассчитывали зоны относительной аккумуляции Трёха (кЪ и ку) и крутизну склонов и использовали дис-криминантный анализ для оценки потенциальной неустойчивости почв (определения мест начала оврагообразования). Помимо рельефа, они учитывали также материнские породы и терморежим склонов. Они нашли, что наиболее тесной является связь мест начала оврагообразования, тестировавшихся авторами экспериментально, с зонами относительной аккумуляции.

Л. Мартз и Е.де Джонг [47] изучали связь

рельефа с пространственной изменчивостью естественных радионуклидов в почвах участка Канадских Прерий. В этом исследовании статистически значимая (при Р<0,05) корреляция концентраций радионуклидов в почве с крутизной не была обнаружена, в отличие от корреляции с высотой, кЪ и ку. Эти авторы анализировали также причины увеличенного содержания песка в бессточных депрессиях, и пришли к заключению, что это вызвано ветровой эрозией, которая перемещает песок на сравнительно небольшие расстояния "прыжками и скачками", в отличие от мелкодисперсных фракций почв, которые ветровая эрозия вначале размещает сравнительно равномерно по территории, а затем они перемещаются поверхностным стоком (водной эрозией) в притальвеговые области, где авторы наблюдали значительное увеличение их содержания.

А. Хаирстон и Д. Григал [36] статистически сравнивали крутизну, ориентацию склонов и одну из характеристик относительного положения в рельефе (длину склона) со свойствами почв и древесной растительности на пологом рельефе в песчаных почвах, покрытых дубовым лесом. Характеристики почв отбирались из четырех холмов, по 4 спускающихся с вершины холма, каждый на одной почвенной картографической единице, трансекты на каждом (ориентация трансект подбиралась СВ, ЮВ, ЮЗ, СЗ), каждый трансект на одном типе растительности, по 3 точки на каждой трансекте - в верхней, средней и нижней части склона. Измерялись ряд свойств почв и лесотаксационные характеристики. Для сравнения данных в 48 точках использовалась статистическая методика АЫОУА. Наиболее тесная статистическая связь наблюдалась с характеристиками относительного положения в рельефе (Р<10-3), значима также была связь с ориентацией склонов (Р<0,04). Связь остальных характеристик почв и растительности с рельефом почти не была статистически достоверной, за исключением Сорг (Р<0,01) и полного N в почве (Р<10-3) с относительным положением, а также полного N в лесной подстилке с ориентацией (Р<0,03). Значимой корреляции

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

Площадь сбора играет немалую роль при описании рек вместе с потенциальной гидросетью [26; 27; 34; 35; 46; 47; 51; 53; 69], с ней статистически связаны водная обеспеченность почв и распространение растительности [60; 75], различные аспекты водной эрозии почв [34; 47; 51; 53] и ущерб от дождевых стоков [30], ряд параметров размывания дорог и неустойчивости склонов [56], условия возникновения рек [45], и она имеет важное значение при построении гидрологических моделей местности [27; 59], моделей уязвимости грунтовых вод от поверхностных загрязнений [38] и моделей зон водного пересыщения почв [56].

Здесь не рассматривается использование различных, часто некорректно определенных МВ, составленных из базовых МВ, см. о них, например, книгу [76].

Исследования по водной эрозии почв

В исследованиях водной эрозии почв рельеф играет особо важную роль, и в этой области разрабатывались достаточно серьезные физико-математические модели. Поэтому полезно рассмотреть историю этого вопроса с точки зрения использования в этих моделях геоморфометрии. В ранних зарубежных и во многих современных отечественных работах использовалось так называемое "универсальное уравнение эрозии" (USLE -от Universal Soil Loss Equation) вида E=RKLSCP, где E есть средняя потеря почв, R - фактор интенсивности дождя, K - почвенный фактор, LS - фактор рельефа (фактор длины склона в USLE), C - фактор растительного покрова и P - фактор превентивных практик. Однако, в 1991 году была осознана физическая бессмысленность длины склона, и было предложено заменить фактор LS (зависящий от крутизны и длины склона) на площадь сбора, имеющую ясный гидрологический смысл [57]. Новое уравнение стало известно как RUSLE, в нем

LS=A•MCAnGFm, где МСА - площадь сбора, GF - фактор крутизны склонов (синус крутизны), А,т,п - числа, эмпирические коэффициенты. Модель RUSLE лучше отражает повреждение почв концентрированным водным потоком [54]. Дальнейшие исследования показали, что область применимости RUSLE весьма ограничена - эта модель применима лишь к участкам размером не более 100 метров без значимой кривизны склона и без значимого осаждения седимента [54]. Спустя два года после нашего исследования по выявлению смысла горизонтальной кривизны [23] полученные результаты, как и сама эта МВ (под названием "тангенциальная"), были пе-реокрыты американскими эрозионистами [52] и введены вместе с вертикальной кривизной и ориентацией склонов в новые модели эрозии, такие как ^РЕБ, SIMWE [53] или WEPP [34], позволявшие анализировать более сложные местности. Горизонтальная и вертикальная кривизны в моделях ^РЕБ и SIMWE, однако, использованы в виде их суммы, т.е. средней кривизны, что не позволяет оценивать относительную роль механизмов аккумуляции. Второй механизм аккумуляции описан только локально (сегодня еще не введена МВ для его регионального описания), и потому эти модели еще нуждаются в дальнейшем развитии. Тем не менее, отметим, что в них уже задействовано пять МВ, а менее детальное описание рельефа делает модели недостаочно эффеткивными.

Х. Митасова с соавт. [53] заключили, что при пространственно однородных осадках, почвах и растительности рельеф определяет на местности общий мало меняющийся "базовый узор" эрозии (то есть больше меняются показатели эрозии, чем этот узор), так что места относительно высокой эрозии остаются приурочены к современному рельефу. Это положение применимо и к моделям водного поверхностного стока, продемонстрированным на сайте авторов [53] цветным "трехмерным" фильмом.

Индикационный принцип

Большая трудоемкость и дорогивизна отбора и анализа образцов почв (или расти-

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

Индикационный принцип состоит в следующем: если некоторое свойство почв или растительности хорошо коррелирует с рельефом, то его можно прогнозировать прямо по рельефу [74]. На практике для этого осуществляются измерения данной характеристики почвы в количестве точек порядка 200 и проводится статистическое сравнение методом множественной регрессии этих характеристик со значениями МВ, рассчитанных в крупном масштабе, в тех же точках,. Если найдена статистически достоверная корреляционная связь, то полученные при сравнении коэффициенты корреляции используются для построения матрицы предсказанных значе-

ний характеристик почв. Качество предсказания пространственной изменчивости почв проверяют сравнением предсказанных значений с измеренными в точках валидации значениями (для чего резервируется около 10% образцов почв - они не участвуют в регрессии) по известным критериям [49]. На рис. 7 показан пример результата расчета мощности гумуса по 143 точкам измерений в 120 тыс. точек с предсказанными в них значениями.

Интенсивные исследования этого рода за рубежом [39: 49; 60 и др.] привели к появлению ряда новых направлений в агрохимии и почвоведении, таких как прецизионное земледелие (районирование одного сельскохозяйственного поля для увеличения урожая), поиск земель с оптимальным плодородием, моделирование разнообразия экосистем и др. (см. [49]). Кроме того, полученные количественные характеристики пространственной изменчивости свойств ландшафтов и их большое значение для смежных с почвоведением дисциплин (например, гидрологии) привели к осознанию, что осуществляемая в настоящее время почвенная съемка не удовлетворяет запросам ее потребителей, в связи с чем в ряде стран разрабатываются новые критерии проведения

Рис. 7. Результаты предсказания пространственной изменчивости мощности горизонта А (гумусового) при шаге решетки м=10м [72], показанные на трехмерной диаграмме рельефа

этой съемки, более адекватные использованию количественных моделей.

Некоторые новые направления

Новые или ранее не использовавшиеся в данной области наук о Земле МВ, представленные в Табл.2, открывают и определенные новые, ранее недоступные возможности. Отметим в этой связи ряд таких недавно возникших направлений.

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

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

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

тик ландшафта в них [24; 66].

В агроландшафте при отсутствии корреляции продуктивности зерновых с рельефом нами установлена значимая корреляция с ним показателей качества урожая [24]. Как показано в наших исследованиях [21], модификация площади сбора для анализа распространения жидкого загрязнителя из одной или более точек позволяет рассчитать по рельефу максимальную область распространения жидкого контаминанта, например, нефти при аварийных разливах. Здесь становится очень существенной негладкость земной поверхности, вызывающая многократные разветвления и слияния потоков нефти. Реальная область распространения уже требует моделирования потерь нефти при ее латеральном переносе на впитывание в почву, испарение и химическое разложение. Методы этого рода перспективны для оценки связанных с латеральным переносом воды и жидких контаминантов экологических рисков, моделирования распространения удобрений на сельхозполях и др.

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

СПИСОК ЛИТЕРАТУРЫ

1. Александров П.С., Пасынков Е.А. Введение в теорию размерности. Введение в теорию топологических пространств и общую теорию размерности. М.: Наука, 1973.

2. ГаршиневЕ.А., ПетрюкБ.П. Автоматизированное выделение эрозионных земельных фондов на водосборах // Современные проблемы земледелия и экологии. Сборник докладов Международной научно-практической конференции, 10-12 сент. 2002 г Курск: Все-росиийский НИИ земледелия и защиты почв от эрозии РАСХН, 2002.

3. Гаусс К.Ф. (1827). Общие исследования о кривых поверхностях // Об основаниях геометрии. М.: Гостехтеориздат, 1956.

4. Дмитриев Е.А. Концепция пластики рельефа и почвоведение // Почвоведение. 1998. №3.

5. Ефремов Ю.К. Опыт морфографической классификации элементов и простых форм рельефа // Вопросы географии. Картография. Сборник 11. М.: Географгиз, 1949.

6. КуряковаГ.А., ФлоринскийИ.В., ШарыйП.А. Корреляция между почвенной влажностью и некоторыми топографическими величинами // Современные проблемы почвенной картографии. Труды конференции. М.: Почвенный институт им. В.В.Докучаева. 1991.

7. Ласточкин А.Н. Рельеф земной поверхности (принципы и методы статической геоморфологии). Л.: Недра, 1991.

8. Неуструев С. С. (1922). Почвы и циклы эрозии // Генезис и география почв. М.: Наука, 1977.

9. НеуструевС.С. (1930). Элементы географии почв // Генезис и география почв. М.: Наука, 1977.

10. Погорелов А.В. Дифференциальная геометрия. М.: Наука, 1974.

11. Пузаченко Ю.Г., Онуфреня И.А., Алещенко Г.М. Количественные методы классификации форм рельефа // Известия АН, Серия геогр. 2002. №26.

12. Сибирцев Н.М. (1900). Почвоведение - курс лекций, читанных студентам Института сельского хозяйства и лесоводства в Новой Александрии // Избранные сочинения. Т. 1. Почвоведение. М.: Сельхозгиз, 1951. 472 с.

13. Соболевский П.К. (1932). Современная горная геометрия // Геометрия структур земной поверхности. Пущино: ПНЦ АН СССР, 1991.

14. Степанов И.Н. Пространство и время в науке о почвах. Недокучаевское почвоведение. М.: Наука, 2003.

15. ТаргульянВ.О. Специфика почвы как поверхностно-планетарной оболочки биосферной планеты // Экология и почвы. Избранные лекции УШ-ГХ Всероссийских школ (19981999гг.). Т. 3. М.: ПОЛТЕКС, 1999.

16. Финн Р. Равновесные капиллярные поверхности. Математическая теория. М.: Мир, 1989.

17. Формы геологических тел (терминологический справочник). М.: Недра, 1977.

18. Фридланд В.М. (1958). К вопросу о зависимости структуры вертикальной зональности почв горных стран от климатических ус-

ловий (на примере Большого Кавказа) // Проблемы географии, генезиса и классификации почв. М.: Наука, 1986. С. 44-54.

19. Шарая Л.С., Шарый П.А. Использование морфометрических статистик для описания внутренней геометрии равнинных и горных местностей // Известия Самарского научного центра РАН. 2003. Т. 5. №2. С. 278-286.

20. Шарая Л.С., Шарый П.А. Элементарные формы в классификациях рельефа, их связь с характеристиками ландшафта Приокско-Террасного заповедника // Известия Самарского научного центра РАН. Спец. выпуск "Природное наследие России". Часть 1. 2004. С. 102-111.

21. Шарая Л.С., Шарый П.А. Аналитические ситуационные игры для уменьшения экологических рисков // Материалы международной конференции "Природное наследие России: изучение, мониторинг, охрана", гТоль-ятти, Россия (21-24 сентября 2004 г). 2004, Тольятти: ИЭВБ РАН.

22. Шарый П.А. Топографический метод вторых производных // Геометрия структур земной поверхности. Пущино: ПНЦ АН СССР, 1991.

23. Шарый П.А., Степанов И.Н. О методе вторых производных в геологии // Доклады АН СССР. 1991. Т.319. №2.

24. Шарый П.А. Оценка взаимосвязей рельеф-почва-растения с использованием новых методов в геоморфометрии // Автореферат... на соискание уч. степени канд. биол. наук. Тольятти, 2005.

25. Aandahl A.R. The characterization of slope positions and their influence on the total nitrogen content of a few virgin soils of Western Iowa // Soil Science Society of America Proceedings, 1948, V13,

26. Beven K. Towards the use of catchment geomorphology in flood frequency predictions / / Earth Surface Processes and Landforms, 1987, V12, N.1,

27. Beven K.J., Kirkby M.J. A physically based, variable contributing area model of basin hydrology // Hydrological Sciences Bulletin, 1979, V24, N.1.

28. Buivydaite V.V.,MozgerisG. Digital land surface analysis: on possibilities of applications in soil survey // "EuroSoil" Conference, 2004. 14p.

Website http://kuk.uni-freiburg.de/hosted/ eurosoil2004/full_papers/ id995_Buivydaite_full.pdf.

29. Cayley A. On contour and slope lines // Philosophical Magazine. 1859. Series 4. V.18. N.120.

30. Costa-CabralM.C., BurgesS.J. Digital elevation model networks (DEMON): A model of flow over hillslopes for computation of contributing and dispersal areas // Water Resources Research. 1994. V.30. N.6.

31. Cummins P.F. Stratified flow over topography: time-dependent comparisons between model solutions and observations // Dynamics of Atmospheres and Oceans. 2000. V.33. N.1.

32. EvansI.S. General geomorphometry, derivatives of altitude, and descriptive statistics // R.J.Chorley (ed.). Spatial Analysis in Geomorphology. London: Methuen & Co., Ltd., 1972. Ch.2.

33. Evans I.S. The effect of resolution on gradients calculated from an altitude matrix // Report 3 on Grant DA-ER0-591-73-G0040, 'Statistical characterization of altitude matrices by computer' (Appendix: Stationarity). Department of Geography, University of Durham, England, 1975.

34. FlanaganD.C., Renschler C.S., Cochrane T.A. Application of the WEPP model with digital geographic information // Problems, Prospects and Research Needs. 4th International Conference on Integrating GIS and Environmental Modeling (GIS/EM4). Banff, Alberta, Canada, September 2 - 8, 2000, GIS/ EM4 No.149. Website: http:// www.colorado.edu/research/cires/banff/ pubpapers/149/

35. Freeman T.G. Calculating catchment area with divergent flow based on a regular grid // Computers and Geosciences. 1991. V.17. N.3.

36. Hairston A.B., Grigal D.F. Topographic influences on soils and trees within single mapping units on a sandy outwash landscape // Forest Ecology and Management. 1991. V.43. N.1-2.

37. Husson O, Verburg P.H., Mai ThanhPhung, Van Mensvoort M.E.F. Spatial variability of acid sulphate soils in the Plain of Reeds, Mekong delta, Vietnam // Geoderma. 2000. V.97. N.1-2.

38. Kalinski R.J., Kelly W.E., Bogardi I., Ehrman R.L., Yamamoto P.D. Correlation between DRASTIC vulnerabilities and incidents of VOC contamination of municipal wells in Nebraska / / Ground Water, 1994, V32, N.1.

39. KingD., BourennaneH., IsampertM, Macaire J.J. Relationship of the presence of a non-calcareous clay-loam horizon to DEM attributes in a gently sloping area // Geoderma. 1999. V. 89. N.1-2.

40. Krcho J. Teoreticka concepcia a interdisciplinarne aplikacie komplexneho digitalneho modelu reliefu pri modelovani dvojdimenzionalnych poli // Geograficky easopis. 1983. V 35. N.3.

41. Lanyon L.T., Hall G.F. Land-surface morphology: 2. Predicting potential landscape instability in eastern Ohio // Soil Science. 1983.

V 136. N.6.

42. Lobb D.A., Kachanoski R.G. Modelling tillage erosion in the topographically complex landscapes of southwestern Ontario, Canada // Soil and Tillage Research. 1999. V. 51. N.3-4.

43. MacMillanR.A., Pettapiece W.W., Nolan S.C., Goddard T.W. A generic procedure for automatically segmenting landforms into landform elements using DEMs, heuristic rules and fuzzy logic // Fuzzy Sets and Systems. 2000.

V 113.

44. Mandelbrot B. How long is the coast of Britain? Statistical self-similarity and fractional dimension // Science. 1967. V.156. N. 3775.

45. Mark D.M. Topology of ridge patterns: Randomness and constraints // Geological Society of America Bulletin. 1979. Part I.V 90.

46. Martz L. W., de Jong E. CATCH: a Fortran program for measuring catchment area from digital elevation models // Computers and Geosciences. 1988. V.14. N. 5.

47. Martz L. W, de Jong E. Natural radionuclides in the soils of a small agricultural basin in the Canadian Prairies and their association with topography, soil properties and erosion // Catena, 1990. V. 17. N.1.

48. MaxwellJ.C. On hills and dales // Philosophical Magazine. 1870. Series 4. V. 40. N. 269.

49. McBratney A.B., Odeh I.O.A., Bishop T.F.A., Dunbar M.S., Shatar T.M. An overview of pedometric techniques for use in soil survey //

Geoderma. 2000. V.97. N.3-4.

50. Miller P.M., Singer M.J., Nielsen D.R. Spatial variability of wheat yield and soil properties on complex hills // Soil Science Society of America Journal. 1988. V. 52.

51. Mitas L., Mitasova H., Brown W.M., Astley M. Interacting fields approach for evolving spatial phenomena: application to erosion simulation for optimized land use // Proceeding of the III International Conference on Integration of Environmental Modeling and GIS (Santa Fe, NM, USA), 1996. National Center for Geographic Information and Analysis: Santa Barbara CA // Website: http:// skagit.meas.ncsu.edu/~helena/gmslab/reports/ CerlErosionTutorial/denix/Advanced/ InteractingFields/Ifields/index.html

52. Mitasova H., Hofierka J. Interpolation by regularized spline with tension. II. Application to terrain modeling and surface geometry analysis // Mathematical Geology. 1993. V. 25. N. 6.

53. Mitasova H., Mitas L., Brown W.M., Johnston D. GIS tools for erosion/deposition modeling and multidimensional visualization. Part IV: Processed based erosion simulation for spatially complex conditions and its applications to installations. Part V: Impact of transport capacity and terrain structures on erosion simulations // Report for USA CERL. Urbana-Champaign: University of Illinois, 1997. 30p. Website: http:/ /skagit.meas.ncsu.edu/~helena/gmslab/reports/ CerlErosionTutorial/denix/Advanced/ ErosionRep97/rep97.html

54. Mitasova H. et al., 1999. Spatial modeling of soil detachment with RUSLE 3d. Website http:/ /skagit.meas.ncsu.edu/~helena/gmslab/reports/ CerlErosionTutorial/denix/ Models%20and%20Processes/RUSLE3d/ RUSLE3d.htm

55.Mitusov A.V., Shary P.A. Application of quantitative land surface analysis methods to soil water content spatial variability studies // K.-P. Seiler, S. Wohnlich (Eds.). New Approaches Characterizing Groundwater Flow. Proceedings of the XXXI International Association of Hydrogeologists Congress, Munich, Germany, September 10-14, 2001. Swets and Zeitlinger Lisse, etc: A.A.Balkema Publishers. 2001. V 2.

56. Montgomery D.R. Road surface drainage,

channel initiation and slope instability // Water Resources Research. 1994. V. 30. N.6.

57. Moore I.D., GraysonR.B., Ladson A.R. Digital terrain modelling: a review of hydrological, geomorphological, and biological applications / / Hydrological Processes. 1991. V. 5.

58. Pike R. A Bibliography of Terrain Modelling (Geomorphometry), the Quantitative Representation of Topography - Supplement 4.0. Open-file report 02-465, U.S. Geological Survey, 2002. - Website: http://wrgis.wr.usgs.gov/wgmt/ sfbay/pike.html

59.RodheA., NybergL., BishopK. Transit times for water in a small till catchment from a step shift in the oxigen 18 content of the water input // Water Resources Research. 1996. V. 32. N. 12.

60. Ryan P. J., McKenzie N.J., O'Connell, D. Loughhead A.N., Leppert P.M., Jackuier D., Ashton L., 2000. Integrating forest soils information across scales: spatial prediction of soil properties under Australian forests // Forest Ecology and Management. 2000. V. 138.N.1-3.

61. Schmidt J., Evans I.S., Brinkmann J. Comparison of polynomial models for land surface curvature calculation // International Journal of Geographical Information Science, 2003, V 17, N.8,

62. Simmons F.W., Cassel D.K., Daniels R.B. Landscape and soil property effects on corn grain yield response to tillage // Soil Science Society of America Journal, 1989. V. 53. N.2.

63. Shary P.A. Land surface in gravity points classification by a complete system of curvatures // Mathematical Geology. 1995. V. 27. N.3.

64. Shary P. A. Analytical GIS Eco, help system and description. 2001. Website http:// members.fortunecity.com/eco4/giseco/, - русская версия доступна на сайте http ://www. esti-map. ru/gis_eco. htm

65. Shary P.A., Sharaya L.S., Mitusov A.V. Fundamental quantitative methods of land surface analysis // Geoderma. 2002. V. 107. N. 12.

66. Shary P.A., Rukhovich O.V, Sharaya L.S., Mitusov A. V. Soils and topography: accumulation zones and non-local approaches // Transactions of 17-th World Congress of Soil Science, "Soil

Science: Confronting New Realities in the 21-st Century", held in Bangkok, Thailand, August 1421, 2002. Volume IV. Symposium 48: "Development in soil data processing", Paper No. 2310.

67. Shary P.A., Sharaya L.S. Mitusov A.V. The problem of scale-specific and scale-free approaches in geomorphometry // Geografia Fisica e Dinamica Quaternaria.2005. V.28.

68. Sinai G., Zaslavsky D., Golany P. The effect of soil surface curvature on moisture and yield -Beer Sheba observations // Soil Science. 1981. V 132. N.5.

69. Speight J.G. Parametric description of land form // Stewart G.A. (Editor), Land Evaluation. London: Macmillan, 1968.

70. Speight J.G. A parametric approach to landform regions // Progress in Geomorphology. Institute of British Geographers special publ. No.7. Oxford: Alden & Mowbray Ltd at the Alden Press, 1974.

71. Timlin D.J., Pachepsky YaA., Snyder V.A., Bryant R.B. Spatial and temporal variability of corn grain yield on a hillslope // Soil Science Society of America Journal. 1998. V. 62. N.3.

72. Thompson J.A., Bell J.C., Butler C.A. Digital elevation model resolution: effects on terrain

attribute calculation and quantitative soil-landscape modeling // Geoderma. 2001. V. 100.

73. Troeh F.R. Landform parameters correlated to soil drainage // Soil Science Society of America Proceedings. 1964. V.28. N.6.

74. Webster R. Quntitative and Numerical Methods in Soil Classification and Survey. Oxford: Oxford Univ. Press, 1977.

75. Wigmosta M.S., Vail L.W., Lettenmaier D.P. A distributed hydrology-vegetation model for complex terrain // Water Resources Research. 1994. V.30. N.6.

76. Wilson J.P., Gallant J.C. (editors). Terrain Analysis: Principles and Applications. New York: John Wiley & Sons Inc., 2000.

77. Young T. An essay on the cohesion of fluids // Philosophical Transactions of Royal Society of London. 1805. V. 95.

78. Young M. Terrain analysis: program documentation // Report 5 on Grant DA-ERO-591-73-G0040, 'Statistical characterization of altitude matrices by computer'. Department of Geography, University ofDurham, England, 1978.

79. Zevenbergen L.W., Thorne C.R. Quantitative analysis of land surface topography // Earth Surface Processes and Landforms. 1987. V. 12. N.1.

GEOMORPHOMETRY IN EARTH SCIENCIES AND ECOLOGY, AN OVERVIEW OF METHODS AND APPLICATIONS

© 2006 P. A. Shary

Institute of Physical, Chemical and Biological Problems of Soil Science of Russin Academy of Sciences, Pushchino

Geomorpometry as a science provides an objective quantitative basis for use of its results in various branches of Earth sciencies and in ecology. Recent developments in geomorphometry by means of integration of results from several scientific disciplines as well as introducing of new quantitative topographic attributes open new horizons for use of these results in ecology and in landscape studies. An overview is given of both current state of geomorphometry itself, and of its applications in these sciencies.