Кратко о Васюганском болоте (стр. 5 )

Более детально рассмотрим данное предположение на примере ИК спектров веществ фенольной природы, которые могут составлять от 3-5 до 41 % от общей суммы водорастворимых ОВ. Сравнительное исследование в торфах и осадках болотных вод ландшафтного профиля ИК-спектров и их спектральных коэффициентов, отражающих соотношение гидрофильной и гидрофобной составляющих в структурах молекул полифенолов, позволило выявить особенности миграции водорастворимых ОВ. Как правило, количество гидроксильных, фенольных гидроксилов, карбоксильных групп и ароматических фрагментов повышается в водорастворимых веществах в августе, сентябре, что объясняется высокой микробиологической активностью в системе: торфяная залежь - болотные воды, прогретой в этот период доС.

В болотной воде автономной части ландшафтного профиля соотношение оптических плотностей гидроксильных групп D3400/D1460 колеблется от 0.89 до 1.49, фенольных гидроксилов D1270/D1460 – 0.78-0.86, карбоксильных групп D1720/D1460 – 1.16-1.28 и ароматических фрагментов D1620/D1460 – 1.20-1.85. В водорастворимых веществах транзитной и трансаккумулятивной части ландшафтного профиля повышается доля фенольных гидроксилов D1270/D1460 до 0.96 и карбоксильных групп D1720/D1460 – 1.73 по сравнению с автономной частью. В реках Ключ и Бакчар для водорастворимых веществ характерно незначительное содержание перечисленных функциональных групп ввиду разбавления мигрирующего потока поверхностными водами, стекающими с не заболоченной части территории. Нало полагать этими причинами объясняются также иные закономерности миграции ОВ в системе сопряженных олиготрофных болот ландшафтного профиля.

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

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

Общий объем выноса химических элементов за период стока составил: Ca2+ - 1398 кг/км2, Feобщ-311, SO42--391, NO3--236, NO2- - 1, Pb–2.253.10-3 , Mn-317.29.10-3, Zn-41.191.10-3, Ni-8.151.10-3, Ti-29.651.10-3 кг/км2. Вынос со стоком растворенного органического вещества оказался равным 6945 кг/км2 или 6.9 гм -2г -1.

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

Во вторых, крайне важным является более детальные исследования ботанического состава торфяного профиля БЭС. Так, анализируя химический состав болотных вод, мы исходили из общего представления строения торфяной залежи по ландшафтному профилю. Было показано, что условно аэробный горизонт (1 м) автономной части профиля сложен сфагновым торфом, транзитная – фускум и медиум торфом, трансаккумулятивная – сосново-пушицево-сфагновым, который сменяется вниз по профилю на древесно-пушицевый переходный. Вместе с тем, например, сфагновая залежь метрового слоя состоит из 9 видов сфагновых мхов и включает также осоки, хвощи и пушицу, а отсюда и разнообразие состава болотных вод. Разные торфообразователи содержат водорастворимых веществ от 3 до 21%. В свою очередь в состав водорастворимых веществ входят моно - и полисахариды, пектиновые вещества. Следовательно, их соотношение в торфообразователях также будет разным и, надо полагать, иной будет и степень полимеризации их ОВ на разных глубинах торфяной залежи, а отсюда и степень подвижности. И это необходимо учитывать.

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

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

При разработке математической модели выноса химических веществ с поверхности водосборного бассейна и их движения по русловой сети принималось во внимание следующее:

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

1. Задача решается в одномерной постановке. Концентрация рассматриваемых ингредиентов принимается осредненной по живому сечению потока или эффективной площади сечения склона для склонового стока, т. е. меняется только по длине и во времени.

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

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

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

где Qн(t) и Cн(t) - соответственно расход воды и средняя концентрация растворенного вещества в замыкающем створе в момент времени t; qбj(t-τ), Cбj(t-τ) - то же для бокового притока с j-го частного бассейна в момент времени (t-τ); pбj(τ) - то же для j-ой кривой добегания бокового притока; N - число частных площадей бассейна, с которых определяется боковой приток.

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

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

где f(τ) - плотность распределения бокового притока по длине реки; mhk - начальный момент k-го порядка времени добегания элементарных объемов бокового притока; mk(τ) - начальный момент времени добегания объема, поступающего на элементарный участок длиной dl = V dt, удаленный на расстояние l от замыкающего створа (τ = l/V, где V - средняя скорость течения на участке); τL – среднее время добегания на всем участке (τL = L/V, где L – длина участка). Первые три момента mhk по формуле (3) выражаются следующим образом:

m1 (τ) = τ;

m2 (τ) = a2τ + τ2; (3)

m3 (τ) = ka4 τ + 3a2τ + τ3,

где a – параметр продольного рассеяния; k – отношение коэффициента асимметрии к коэффициенту вариации времени добегания.

Задаваясь различными видами функции f(τ), можно получить выражения для моментов, отвечающие различным случаям распределения бокового притока по длине участка. Функцию f(τ) можно представить как отношение

где q(τ) - боковой приток на единице длины l в единицу времени (длина выражена в единицах среднего времени добегания τ = l/V); С(τ) - концентрация примеси.

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

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

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

Разработанная модель была проверена на примере выноса гуминовых кислот (рис16).

Рис.16 Фактический и рассчитанные гидрографы выноса гуминовых кислот

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

Биохимические процессы. Динамика биохимических процессов изучалась в вегетационные периоды разных по гидротермическим условиям лет. Ежемесячно по профилю торфяных залежей пунктов изучения (пп. 2,3,5) отбирались образцы, в которых изучали микробиологическую и энзимологическую активности. Методом посева проводили определение численности следующих групп микроорганизмов: аммонификаторов, микроорганизмов, разрушающих минеральные формы азота, целлюлозоразрушающих микроорганизмов (аэробы, анаэробы), метанотрофов. Определение биомассы микроорганизмов проводили методом люминесцентной микроскопии. Метод позволяет получить представления о микробном пуле, содержащемся в торфяной залежи.

Так, отмечается закономерность в динамике целлюлозоразрушающих микроорганизмов – увеличение их численности (до 1,5-2 раз) в более благоприятный по гидротермическим условиям годы ( гидротермический коэффициент ГТК равен 0,7). Активность аммонификаторов в большей степени проявлялась во влажные годы.

Более подробно рассмотрим влияние гидротермических условий на активность биохимических процессов по динамике микробной биомассы. Наибольший прирост микробной биомассы отмечался в сухой год с ГТК 0,7. При этом экстремальные значения численности микробной биомассы в 2 – 3 раза были выше в верхнем слое торфяной залежи по сравнению с более глубокими ее слоями. В то время как в менее благоприятные по гидротермическим условиям годы (ГТК 1,2-1,5) микробная биомасса, например, в торфяной залежи ряма активно увеличивалась в течение вегетационного периода. В торфяной залежи осоково-сфагновой топи различия по месяцам особенно проявляются во влажные годы с близкими к поверхности УБВ. Таким образом, численность в торфяной залежи микрофлоры различается в зависимости от положения в ландшафтном профиле (разный генезис), глубины торфяной залежи, гидротермических условий года и отдельных месяцев вегетационных периодов.

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

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

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

Эмиссия СО2 и CH4. В течение вегетационных периодов проводили определение концентрации СО2 и СН4 в торфяных залежах ландшафтного профиля «peepers»-методом. Через месяц камеры извлекались из торфяной залежи и газовый состав анализировался в лаборатории на хроматографе «Кристалл-2000». Измерения эмиссии СО2 и СН4 с поверхности торфяных залежей проводили ежемесячно камерно-статистическим методом.

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

В условиях прохладного и влажного года с ГТК 1,5 динамика выделения СО2 не превышает 120 мгСО2/м2час с наименьшими значениями в осоково-сфагновой топи. В более благоприятных гидротермических условиях с ГТК 1,2 не обнаруживается различий в динамике выделения СО2 из торфяных залежей ряма и осоково-сфагновой топи. В сухой год при низких УБВ выделение СО2 характеризуется высокими значениями в мае и июне в осоково-сфагновой топи. Следует заметить, что динамика выделения СО2 в большой степени определяется температурой поверхностного слоя торфяной залежи. Определение эмиссии СО2 в течение дня в разные месяцы вегетационного периода также подтвердили такую закономерность.

Исследования эмиссии СО2 показало, что в целом по ландшафтному профилю наблюдается снижение интенсивности выделения СО2 от высокого ряма к открытой осоково-сфагновой топи. Средние значения интенсивности выделения СО2 за многие годы исследований составили в пп.2, 3, 5 соответственно – 133, 77, 59 мг/м2час (рис.17)

Примечание: годы

по гидротермическому коэффициенту: ГТК –0,8, ГТК – 0,5, ГТК – 1,0, ГТК – 1,3.

Рис. 17. Средние значения эмиссии СО2 в различных биогеоценозах

В изучаемых БГЦ, наибольшие пределы изменения выделения СО2 отмечались в высоком ряме от 29 до 356 мгСО2/м2час и наименьшие – в осоково-сфагновой топи от 0 до 208 мгСО2/м2час.

Наибольших значений эмиссия метана достигает во влажные годы (3,2-11,8 мгСН4/м2час), далее следует средний (ГТК 1,2; 1,2-8,6 мгСН4/м2час) и сухие годы (0,7-9,3 мгСН4/м2час). Для почв ландшафтного профиля наблюдается прямая зависимость между уровнем болотных вод и выделением СН4 в атмосферу. Снижение УБВ приводит к уменьшению выделения метана. Как и эмиссия диоксида углерода, интенсивность выделения СН4 олиготрофными торфяными залежами достоверно определяется температурой горизонта 0-50 см (r = 0,8), а также численностью метанотрофов (r = 0,9).

7. Заключение

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

Согласно концепции устойчивого развития в каждой стране должны быть разработаны критерии устойчивого природопользования применительно к конкретным природным и экономическим условиям.

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

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

Проведенные нами исследования показали, что Васюганское болото продолжает увеличиваться, захватывая прилегающую к нему территорию. Так, поступление углерода в виде органического вещества растений по разным биоценозам ландшафтного профиля составляет от 98,9 до 161,7 гС/(м2 год), или в среднем 133 гС/(м2 год). Если оценить общий вынос углерода из Васюганского болота в виде СО2, СН4 и углерода болотных вод (в среднем 45,1 -136,6 гС/(м2 год)), то расход углерода определяется суммарной величиной, равной 77,4 гС/(м2 год), что значительно меньше прихода углерода.

8. ОСНОВНОЕ О ВАСЮГАНСКОМ БОЛОТЕ

Площадь болота – 5 га; наибольшая протяженность с запада на восток 573.0 км и с севера на юг около 320 км.

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

Васюганское болото находится в области тектонического поднятия.

Агрессивность заболачивания на Васюганском болоте очень высокая. В среднем ежегодно заболачивается 1 800 га. Поэтому 25 % территории Васюганского болота приходится на заболоченные участки, возраст которых не превышает 500 лет при нижнем пределе возраста 9 000 лет;

На Васюганском болоте сосредоточено 18,7 млрд. т. торфа, что составляет 16% от запасов всего Западно-Сибирского региона. Из них разведано 25,5%.

Содержание депонированного углерода в Васюганском болоте достигает 5,1млрд. т. или 12% от депонированного углерода в торфяных залежах всего Западно-Сибирского региона и 4,4 % от депонированного углерода России.

С Васюганского болота стекает более 200 рек.

Васюганское болото расположено в двух природно-геохимических подзонах: южнотаежной и лесостепной.

На Васюганском болоте центральная часть верховых болот возвышается на 7,5 – 10 м над его краями.

Васюганское болото является единственным местом широкого распространения переходных болот.

Рямы (местное название) - сибирские олиготрофные болота с сосново-кустарничково-сфагновыми группами формаций.

Веретья (местное название) - однообразные осоково-гипновые топи пересекают в направлении, перпендикулярном поверхностному стоку, узкие (1-2 м шириной) и длинные (до 1 км длиной) веретья. Веретья возвышаются над осоково-гипновыми топями на 10-25 см. На веретьях поодиночке или небольшими группами растут Betula pubescens и Pinus sylvestris, из кустарников Salix lapponum и S. rosmarinifolia. Кустарничковый ярус развит довольно пышно. Травяной ярус веретий по флористическому составу мало отличается от травяного покрова осоково-гипновых понижений.

Шеломочки (местное название) - по поверхности осоково-гипновых топей разбросаны отдельные островки олиготрофных сосново-кустарничково-сфагновых фитоценозов, диаметр которых колеблется в пределах нескольких десятков метров (“шеломочки”), которые возвышаются над поверхностью осоково-гипновых болот на 50-90 см

Важным фактором является образование на Васюганском болоте переувлажненных низинных гипново-осоковых болот на самой вершине Васюганской гряды с наивысшей для Западно-Сибирской равнины отметкой – 146 м над уровнем моря.

Особенностью Васюганского болота является наличие особых веретьево-топяных сетчато-полигональных низинных болот.

Запасы микробной биомассы в отдельных слоях торфяных залежей достигают 9.82 мг/г.

В торфяных залежах Васюганского болота выделено и идентифицировано 36 культур микромицетов.

Список литературы

1.  Биогеохимические аспекты болотообразовательного процесса. – Новосибирск: Наука, 1986.

2.  Большая Советская Энциклопедия. –М., 1971, том 4, с. 330

3.  Верховые болота Нарымского края (бассейн р. Васюган). // Труды научно-исследовательского торфяного института, 1930, вып.3.

4.  Типовые болота на южной окраине Западно-Сибирской равнинной тайги. // Почвоведение, 1936, №2.

5.  Принципы классификации природных геосистем по устойчивости к техногенезу и прогнозное ландшафтно-геохимическое районирование.// Устойчивость геосистем. – М.: Наука, 1983.

6.  Описание Васюганской тундры. // Записки Зап.-Сиб. отдела Русского геогр. общества, 1884, кн. 6.

7.  Материалы по почвоведению и геологии западной части Нарымского края. // Труды Переселенческого управления, - Петроград, 1915.

8.  , , Запасы углерода в экосистемах лесов и болот России. – Красноярск, 1994.

9.  Водообмен в болотных ландшафтах. –Л.: Гидрометеоиздат, 1975.

10.  Системы показателей для оценки загрязнения тяжелыми металлами. // Агрохимия, 1995, №1.

11.  Природа Нарымского края (рельеф, геология, ландшафты, почвы). // Материалы по изучению Сибири. - Томск, 1930, т.2.

12.  , , Порохина, Е. В., Боровкова торфяные месторождения, их характеристика и функционирование. Изд-во ТГПУ, Томск.

13.  , Головченко микробоценоза в торфяных залежах ландшафтного профиля олиготрофного торфогенеза. Сибирский экологический журнал. 2007, № 3.

14.  , , Инишев болото (природные условия, структура и функционирование. 2-е изд. Томск, ЦНТИ., 2003 с.

15.  , Сергеева образования и эмиссия метана в олиготрофных ландшафтах Васюганского болота. Вестник ТГПУ, вып. 6(57), серия естественные и точные науки, 2006.

16.  , Эколого-геохимическая оценка торфов юго-востока Западно-Сибирской равнины. // География и природные ресурсы, 1999, №1.

17.  , , Инишев органических веществ в системе геохимически сопряженных болотных ландшафтов. Геохимия, 2005, № 2.

18.  , Вернадского о преобразовании биосферы и экологии человека. - М.: Знание, 1986.

19.  , Связь болотообразования с условиями развития рельефа и неотектоникой Барабы. // ДАН СССР, 1952, т.87, №2.

20.  Голоцен СССР. – М.: Наука, 1977.

21.  , Карта современных вертикальных движений земной коры южной части Сибири. // Современные движения земной коры. –М.: Наука, 1980.

22.  О болотах Нарымского края Томской губернии. // Болотоведение, 1915, №1.

23.  , , и др. Болотные системы и их природоохранное значение. –М., 2001.

24.  О работах Западносибирской торфоразведочной экспедиции. // Сб. статей по изучению торфяного фонда. Главное управление торфяного фонда при Совете Министров РСФСР – М., 1957.

25.  , , Гидрогеологические условия формирования болот. // Научные предпосылки освоения болот Западной Сибири. - М.:Наука, 1977.

26.  Краткая история исследований. // Научные предпосылки освоения болот Западной Сибири, - М.: Наука, 1977.

27.  Мировой природный феномен – заболоченность Западно-Сибирской равнины.// Изв. АН СССР, сер. геграфия, 1971, №1.

28.  , Окислительно-восстановительный режим некоторых почв дерново-подзолистой зоны // Агрохимия, 1974, №3.

29.  , Сводка отчетных данных по обследованию в гг. левобережной части Оби Нарымского края, 1910.

30.  Гидрофизика болот. –Л.: Гидрометеоиздат, 1961.

31.  О торфяных месторождениях Западно-Сибирской низменности. // Труды Томского университета, 1957, т.141.

32.  Районирование торфяных месторождений. // В кн. Торфяные месторождения Западной Сибири. –М., 1957.

33.  , , и др. Составление карты скоростей современных вертикальных движений земной коры Западно-Сибирской плиты. // Современные движения земной коры. Морфоструктуры, разломы, сейсмичность. –М.: Наука, 1987.

34.  Поездка на Васюган и Чижапку в Нарымском крае. // Томские губернские ведомости, 1877, №№ 46, 48, 50, 58.