Text (PDF):
Read
Download
Введение Бессточное минерализованное оз. Балхаш расположено в аридной зоне на юго-востоке Казахстана. Среднегодовое количество атмосферных осадков (по метеостанции Балхаш) составляет 129,8 мм. Среднегодовая температура воздуха – 5,9 ºС. Питание озера осуществляется в основном за счет рек, впадающих с юга – Иле, Каратал, Аксу, Лепсы. Реки берут начало в горах Джунгарского Алатау в зоне ледников. Для оз. Балхаш, как и для других водоемов аридной зоны, характерны выраженные межгодовые и межсезонные колебания уровня и минерализации воды. Гидроценозы, представленные преимущественно эврибионтными и эвригалинными видами, хорошо приспособлены к меняющимся гидролого-гидрохимическим условиям этого водоема. Количественные показатели водных сообществ подвержены значительным межгодовым и сезонным флуктуациям. Размах колебаний определяется изменением комплекса абиотических и биотических факторов, в том числе и гидролого-гидрохимическим режимом. Для различных сообществ первостепенное значение имеют различные факторы. Спорадические исследования гидробионтов оз. Балхаш начались в 1929 г. Систематические ежегодные исследования фитопланктона проводятся с 1971 г., зообентоса и зоопланктона – с 1961 и 1963 г. соответственно. Основное внимание уделялось изучению видового состава, выделению доминантных видов и изменению уровня кормности по гидрохимическим районам озера [1–3]. Целью исследований являлся анализ многолетней динамики водных сообществ (фитопланктон, зоопланктон, зообентос) оз. Балхаш в зависимости от изменения абиотических факторов. Материал и методики исследований Для анализа многолетней динамики гидробионтов использованы следующие данные: по фитопланктону – с 1972 по 2012 г., по зоопланктону – с 1961 по 2012 г., по зообентосу – с 1963 по 2012 г. Гидробиологический материал отбирали в летний период по стандартной сетке из 60 станций с охватом всей акватории озера. Отбор и обработку проб фитопланктона, зоопланктона и зообентоса проводили стандартными методами [4]. Для анализа связи между биотическими и абиотическими факторами использовали данные по уровню и минерализации воды, среднегодовой температуре воздуха (по метеостанции Балхаш) за соответствующие годы. Данные по суммарному объему речного стока и количеству атмосферных осадков, выпавших на водное зеркало, площади оледенения Джунгарского Алатау и объему ледникового стока рек заимствованы из работ [5, 6]. Недостающие данные по минерализации воды были восстановлены с помощью функции «Тенденция». Соответствие распределения эмпирических рядов данных нормальному проверяли с помощью критерия хи-квадрат [7]. Для сглаживания данных использовали их логарифмирование и скользящее осреднение с шагом 3 и 5 лет. Значения коэффициентов корреляции Пирсона и Спирмэна находили при уровне значимости p < 0,001…0,05. Для выявления основных трендов в многолетней изменчивости абиотических (уровень и минерализация оз. Балхаш) и биотических параметров (количественные показатели гидроценозов) строили кумулятивные кривые. Для этого находили среднемноголетнее значение каждого показателя, отклонение каждого члена ряда от среднемноголетнего значения и далее – накопленную сумму отклонений. Статистическую обработку данных проводили с помощью программ Excel и Statistica 6. Гидрологический и гидрохимический режимы озера Динамика гидрологического режима оз. Балхаш описана в [5, 8–10]. В изменении уровня воды хорошо выражена цикличность (рис. 1). С начала инструментальных наблюдений (1929 г.) уровень озера снижался и в 1946 г. достиг минимальной отметки в 340,7 м абс. Далее, с 1947 по 1961 г., наблюдался рост уровня воды. После достижения максимальной отметки в 343,0 м в 1962–1970 гг. произошла его некоторая стабилизация (342,60–342,86 м) и далее последовал спад, продолжавшийся до 1986 г. Существенное влияние на гидрологический режим озера в этот промежуток времени оказало наполнение Капшагайского водохранилища на р. Иле, начавшееся в 1970 г. При его отсутствии трансгрессивная фаза оз. Балхаш продолжилась бы вплоть до 1973 г., а последующие колебания уровня происходили бы вблизи более высокой отметки [5]. В последующие десятилетия вновь наступила трансгрессивная фаза с максимумом в 2005 г., небольшим снижением в последующие несколько лет и очередным повышением уровня в 2010–2012 гг. Средняя для озера минерализация воды изменялась также циклично в противофазе с изменением уровня. Таким образом, по данным инструментальных наблюдений, полный цикл в динамике уровня оз. Балхаш (при отсутствии Капшагайского водохранилища) составил от минимума до минимума 43 лет, от максимума до максимума – 51 год. Согласно теории внутривековой изменчивости компонентов общей увлажненности А. В. Шнитникова [11], первая из указанных регрессивных фаз началась не в 1929 г., а раньше – в 1910 г. С учетом этих данных полный гидрологический цикл оз. Балхаш от минимума до минимума составил 64 года. Фаза повышения уровня озера, очевидно, еще не закончилась, даже на фоне существенного изъятия стока р. Иле для заполнения водохранилища на территории Китайской Народной Республики. Среднемноголетняя отметка уровня за 1929–2012 гг. составила 341,78 ± 0,07 м абс. Между состоянием уровня озера и величиной минерализации воды существует тесная линейная связь. По данным 1952–2012 гг., коэффициент корреляции Пирсона между этими показателями r = –0,887. Близкое значение коэффициента корреляции (r = –0,88) между уровнем и минерализацией воды было получено для Западного Балхаша при анализе данных за менее продолжительный период (1960–1977 гг.) [12]. Наличие тесной связи между указанными показателями позволило восстановить недостающие данные по минерализации воды за период 1929–1951 гг. С учетом восстановленных данных среднемноголетняя (1929–2012 гг.) величина минерализации воды составила 2,76 ± 0,03 г/дм3. Рис. 1. Многолетняя динамика минерализации и уровня воды оз. Балхаш Как известно, уровень озера является функцией общей увлажненности территории бассейна, зависящей от климатических особенностей региона. Они определяют количество атмосферных осадков, выпавших на территории водосборного бассейна и над всей площадью озера, а также общий объем стока питающих его рек. Реки, впадающие в оз. Балхаш, берут свое начало на склонах Джунгарского Алатау в зоне ледников. Доля ледникового стока составляет в среднем 13,1 % общего и определяется количеством и площадью ледников [6]. Одной из важнейших климатических характеристик, влияющих на состояние уровня озер аридных зон, является температура воздуха. В частности, от нее зависит испаряемость воды с поверхности водного зеркала. Согласно значениям критерия хи-квадрат, распределение данных не соответствовало нормальному, в том числе и после их преобразования. Для выяснения корреляционных зависимостей между отметкой уровня и перечисленными выше гидролого-климатическими факторами использовали непараметрические методы. Статистически значимые зависимости получены между уровнем оз. Балхаш и следующими показателями: площадью оледенения Северной Джунгарии, объемом ледникового стока, суммарным количеством атмосферных осадков, выпавших на водное зеркало (табл. 1). Наличие преимущественно слабых связей между изменениями уровня и гидролого-климатическими характеристиками обусловлено инерционностью ответной реакции такого большого водоема, как оз. Балхаш. Многолетняя динамика сообществ гидробионтов Фитопланктон является основой первичной продуктивности любого водоема. Микроводоросли утилизируют поступающие в водную экосистему питательные вещества и передают их на более высокие трофические уровни. Таксономическое разнообразие фитопланктона за период исследований (1972–2012 гг.) изменялось от 82 до 217 наименований, с наибольшим разнообразием диатомовых [13, 14]. Среднемноголетняя биомасса фитопланктона составила 0,94 ± 0,07 г/м3. Величина показателя по годам варьировала от 0,17 до 1,86 г/м3. Минимальные и близкие к минимальным значения биомассы отмечались в 1974–1975, 1986, 1999 гг. Максимальные и близкие к максимальным величины показателя зафиксированы в 1992, 1994, 2012 гг. Период 1972–1995 гг. характеризовался в среднем более высокими значениями биомассы фитопланктона по сравнению с последующими десятилетиями (рис. 2, а). Доминировали чаще всего крупные диатомовые. Субдоминировали мелкие формы синезеленых. Таблица 1 Коэффициенты ранговой корреляции Спирмэна между уровнем оз. Балхаш и гидролого-климатическими характеристиками ПараметрS олед ДАVледнVобщАО (м. Балхаш)АО зеркалоt, ºC, среднегодовая (м. Балхаш) Коэффициент ранговой корреляции Спирмэна0,7920,4090,222–0,090,369–0,225 Уровень значимостиp < 0,001p < 0,001p < 0,07p < 0,526p < 0,009p < 0,120 N444463506349 Примечание. S олед ДА – площадь горного оледенения Джунгарского Алатау; Vледн – объем ледникового стока; Vобщ –суммарный сток рек; АО зеркало – количество атмосферных осадков, выпавших на водное зеркало. а б Рис. 2. Многолетняя динамика количественных показателей: а – фитопланктона; б – зоопланктона оз. Балхаш Зоопланктон представлен коловратками, ветвистоусыми и веслоногими ракообразными. За весь период гидробиологических исследований (1929–2004 гг.) в составе зоопланктона было выявлено 312 таксонов [15–19]. Численность зоопланктона за рассматриваемый период варьировала в пределах 19,2– 124,3 тыс. экз./м3. Биомасса сообщества изменялась от 0,50 до 3,14 г/м3 (рис. 2, б). Выявлен тренд увеличения количественных показателей планктонных беспозвоночных, особенно выраженный после зарегулирования стока р. Иле. Ранее нами были показано, что средние значения численности и биомассы зоопланктона в периоды до и после зарегулирования стока р. Иле различались статистически значимо [20]. Среднемноголетняя численность планктонных беспозвоночных составила 58,6 ± 3,8 тыс. экз./м3, при величине биомассы 1,47 ± 0,10 г/м3. Основу количественных показателей формировали веслоногие и ветвистоусые ракообразные. В состав доминантного комплекса входили Arctodiaptomus salinus, Diaphanosoma lacustris, Mesocyclops leuckarti, Thermocyclops crassus, в восточной части еще и Daphnia galeatа. Зообентос представлен аборигенными и акклиматизированными видами червей, моллюсков и членистоногих. Разнообразие донного сообщества изменялось по годам исследований от 91 до 98 видов [13]. Среднемноголетняя биомасса донного сообщества с 1939 по 2012 г. составила 8,19 ± 1,47 г/м3. Величина показателя изменялась от 1,30 до 49,45 г/м2. До акклиматизации донных беспозвоночных в бентофауне доминировали личинки хирономид. С 1960 по 1966 г. ведущая роль в формировании биомассы макрозообентоса принадлежала высшим ракообразным (мизидам и корофиидам). С 1996 г. проявилась тенденция к резкому увеличению биомассы донного сообщества (рис. 3) за счет моллюска Monodacna colorata. Анализ описательных статистик сообществ показал, что многолетнее распределение биомассы макрозообентоса отличалось от нормального. Среднемноголетняя величина биомассы этого сообщества характеризовалась наибольшей дисперсией, асимметрией и коэффициентом вариации (табл. 2). Распределение значений биомассы фитопланктона и зоопланктона также отлично от нормального, хотя и в меньшей степени. Считается, что логарифмирование данных приближает их распределение к нормальному. Этот способ выравнивания значений биомасс фитопланктона и зоопланктона привел к увеличению асимметрии и показателей вариации. Скользящее осреднение данных с шагом 5 лет снизило межгодовую вариабельность (значения стандартной ошибки, стандартного отклонения, дисперсии, коэффициента вариации) и асимметрию данных для зоопланктона и фитопланктона. Оценка распределения данных с помощью критерия хи-квадрат показала, что наибольшее приближение к нормальному распределению дает скользящее осреднение данных для фитопланктона и зоопланктона и логарифмирование – для зообентоса. Рис. 3. Многолетняя динамика биомассы макрозообентоса оз. Балхаш Таблица 2 Описательная статистика для сообществ гидробионтов оз. Балхаш при различных способах выравнивания данных СтатистикиФитопланктонЗообентосЗоопланктон Фак-тическиеСкользящее среднееlogФак-тическиеСкользящее среднееlogФак-тическиеСкользящее среднееlog X0,940,940,288,197,151,671,471,480,27 m0,070,040,021,471,060,120,090,070,07 Mе0,990,940,304,013,851,391,331,590,29 Mo1,00–0,302,63,650,960,5––0,69 σ0,430,230,1010,27,050,830,680,500,51 σ20,180,070,01103,749,70,700,460,250,26 As0,25–0,04–0,232,842,261,020,31–0,12–0,32 N413741484448514751 CV7,14,635,417,914,949,86,55,0186,7 Хи-квадрат4,857,105,2082,0052,0017,98,48,8010,8 df322323544 p0,250,030,070,000000,00000,00050,1350,060,03 Примечание. X – среднее значение; m – стандартная ошибка; Mе – медиана; Mo – мода; σ – стандартное отклонение; σ2 – дисперсия выборки; As – асимметрия; N – число членов в вариационном ряду; CV – коэффициент вариации. Для выяснения корреляционных зависимостей между количественными показателями водных сообществ и абиотическими факторами находили коэффициенты ранговой корреляции Пирсона и Спирмэна при различных способах сглаживания исходных данных. Зависимости между многолетней динамикой биомасс всех сообществ и такими абиотическими факторами, как максимальная, минимальная, среднегодовая температура, количество атмосферных осадков были статистически незначимы. Выявлена статистически значимая отрицательная связь между площадью горного оледенения и численностью (R = –0,627, p < 0,001) и биомассой зоопланктона (R = –0,724, p < 0,001). При скользящем трехлетнем осреднении данных получены статистически значимые зависимости между объемом ледникового стока и биомассой фитопланктона (R = –0,61, p < 0,05) и зоопланктона (R = –0,57, p < 0,05) (рис. 4), а также между величинами биомассы сообществ и гидролого-гидрохимическими показателями – уровнем и минерализацией воды оз. Балхаш (табл. 3). Согласно шкале Чедокка [7], между многолетней динамикой биомассы зообентоса и уровнем оз. Балхаш существует слабая статистически значимая связь. Умеренная по силе отрицательная связь выявлена между биомассой зоопланктона и уровнем, положительная, также умеренная по силе связь, – между биомассой зоопланктона и минерализацией воды. Более тесные линейные связи с абиотическими факторами получены при скользящем осреднении данных. Для фитопланктона и зоопланктона зависимости были аналогичными по знаку – отрицательные с уровнем и положительные – с минерализацией воды. Противоположные по знаку относительно других сообществ значения коэффициентов корреляции с абиотическими факторами получены для макрозообентоса. Значения коэффициентов корреляции Пирсона и Спирмэна, полученные при различных способах сглаживания данных (табл. 4), не противоречили друг другу, хотя были неравнозначными по силе связи. а б Рис. 7. Динамика биомассы: а – фитопланктона; б – зоопланктона в зависимости от объема ледникового стока Таблица 3 Коэффициенты ранговой корреляции Спирмэна между биомассой фитопланктона, зоопланктона, зообентоса и гидролого-гидрохимическими показателями оз. Балхаш ПоказательФитопланктонЗообентосЗоопланктон NRt (N – 2)p-levelNRt (N - 2)p-levelNRt (N – 2)p-level *Биомасса & уровень41–0,290–1,8910,066480,2992,1220,03951–0,398–3,0320,004 *Биомасса & минерализация410,2581,6650,10448–0,233–1,6220,112510,4943,9750,001 **Биомасса & уровень37–0,636–4,8710,001440,4663,4150,00147–0,438–3,2650,002 **Бомасса & минерализация370,5884,3000,00144–0,323–2,2150,032470,5023,8950,001 Примечание. N – число членов в ряду; R – коэффициент ранговой корреляции Спирмэна; t – критерий Стьюдента; p-level – уровень значимости; * – фактические данные; ** – скользящее среднее с шагом осреднения 5 лет. Таблица 4 Коэффициенты корреляции Пирсона и Спирмэна между биомассой фитопланктона, зоопланктона, зообентоса и гидролого-гидрохимическими показателями оз. Балхаш ПоказательФитопланктонЗообентосЗоопланктон log (r)log (R)CC(r)СС(R)log (r)log (R)CC(r)СС(R)log (r)log (R)CC(r)СС (R) Биомасса & уровень–0,280–0,290–0,630-0,6360,3600,2990,4200,466–0,370–0,398–0,440–0,438 Биомасса & минерализация0,2700,2600,5900,588–0,210–0,233–0,160–0,3230,5300,4940,6200,502 Примечание. log – логарифм значений; СС – скользящее среднее с шагом 5 лет; r – коэффициент корреляции Пирсона; R – коэффициент корреляции Спирмэна. Наглядное представление о выявленных корреляционных зависимостях дают кумулятивные кривые, показывающие основное направление изменчивости во времени количественных показателей водных сообществ на фоне изменчивости факторов среды (рис. 5). Во всех случаях отмечался аналогичный по направленности временной сдвиг изменчивости биомассы фитопланктона и зоопланктона по отношению к колебаниям уровня и минерализации воды. По сравнению с величиной биомассы изменчивость численности зоопланктона в межгодовом аспекте была выражена сильнее. а б в г Рис. 5. Кумулятивные кривые, характеризующие изменчивость биомассы: а, б – фитопланктона; в, г – зоопланктона в зависимости от уровня и минерализации воды оз. Балхаш Многолетняя изменчивость биомассы макрозообентоса на фоне изменчивости гидролого-гидрохимического режима озера носила иной, по сравнению с планктонными сообществами, характер. На большей части рассматриваемого отрезка времени сумма накопленных отклонений биомассы донного ценоза от среднемноголетних значений была хорошо согласована с аналогичными отклонениями уровня озера и, соответственно, находилась в противофазе с изменениями минерализации воды (рис. 6). Так же как и для планктонных сообществ, в изменчивости биомассы макрозообентоса прослеживался сдвиг во времени относительно выявленного тренда абиотических параметров. Между многолетней динамикой биомасс фитопланктона и зоопланктона, зоопланктона и зообентоса зависимости были статистически незначимы. Сравнение кумулятивных кривых, характеризующих многолетнюю изменчивость биомасс фитопланктона и зоопланктона во временном интервале 1972–2012 гг., выявило их одинаковую направленность (рис. 7, а), со сдвигом в изменчивости показателей одного сообщества относительно другого на 6–7 лет. Биомасса донного сообщества изменялась в противофазе с изменчивостью планктонных сообществ. Выявлена статистически значимая отрицательная связь между биомассой фитопланктона и зообентоса (R = –0,459, p < 0,003). а б Рис. 6. Кумулятивные кривые, характеризующие изменчивость биомассы макрозообентоса оз. Балхаш: а – в зависимости от уровня; б – от минерализации воды Рис. 7. Многолетняя изменчивость биомассы фитопланктона, зоопланктона и зообентоса оз. Балхаш: а – кумулятивные кривые; б – скользящее осреднение данных и сдвиг кривой биомассы фитопланктона вперед на 6 лет При скользящем осреднении данных и временном сдвиге в 6 лет наблюдалась хорошая согласованность между ходом кривых биомасс фитопланктона и зоопланктона (рис. 7, б). Взаимно противоположный согласованный ход кривых биомасс фитопланктона и макрозообентоса наблюдается без временного сдвига. Заключение Бессточные водоемы являются хорошими индикаторами изменчивости общей увлажненности территории, обусловленной климатическими флуктуациями. Циклические изменения уровня оз. Балхаш – яркий тому пример. О климатической обусловленности изменчивости уровня оз. Балхаш свидетельствует хорошая согласованность с такими индикаторами климатических флуктуаций, как площадь горного оледенения и объем ледникового стока рек. Изменчивость уровня воды оказывает влияние на различные стороны жизни водоема. Наиболее тесно связанным с уровнем параметром является минерализация воды. От изменчивости этих двух факторов, в свою очередь, зависит многолетняя вариабельность биотических параметров. Оценка влияния абиотических факторов на биотические параметры вызывает определенные трудности. К ним можно отнести ошибки отбора данных, обусловленные даже небольшим отклонением срока отбора гидробиологических проб в разные годы, а также климатическими условиями разных лет. Чрезвычайно важным моментом является несоответствие распределения эмпирических рядов данных нормальному, даже после их преобразования. Между тем большинство статистических процедур применимы только при соблюдении условия нормальности исходных или преобразованных данных. Непараметрические критерии менее разнообразны и отличаются меньшей статистической мощностью по сравнению с их параметрическими аналогами. Использование коэффициента ранговой корреляции Спирмэна для анализа связи между непреобразованными данными выявило слабые и умеренные по силе зависимости между биомассой зообентоса и уровнем, между биомассой зоопланктона и уровнем, между численностью и биомассой зоопланктона и минерализацией воды. Отсутствовали связи между многолетней изменчивостью биомассы фитопланктона и абиотическими параметрами. Сглаживание данных методом скользящего среднего показало наличие более сильных статистически значимых связей между абиотическими факторами и количественными показателями всех сообществ, при сохранении знака коэффициентов корреляции. Значения коэффициентов корреляции Пирсона и Спирмэна, полученные при различных способах сглаживания данных, не противоречили друг другу, хотя и были неравнозначными по силе связи. Более высокие значения коэффициентов корреляции Спирмэна при сглаживании данных скользящим осреднением связаны, очевидно, с меньшим влиянием ошибок ежегодного отбора данных и вариабельности климатических условий. Наличие отрицательной зависимости между уровнем воды и биомассой планктонных сообществ можно объяснить следующим образом. Основным лимитирующим фактором для развития растительных клеток является содержание соединений азота и фосфора. Межгодовая динамика концентраций этих веществ в воде водоема определяется их поступлением с территории водосбора и содержанием биогенных соединений в самом водоеме. Ранее нами было показано, что если количество биогенных элементов, накопленных в водоеме, меньше, чем поступает с территории водосбора, между уровнем озер и количественными показателями зоопланктона отмечается прямая зависимость [20]. Для водоемов с ледниковым питанием впадающих рек эта зависимость обратная, т. к. в многоводные годы происходит разбавление озерной воды речными водами, бедными питательными соединениями. Очевидно, это справедливо и для оз. Балхаш с ледниковым питанием впадающих в него рек. Подтверждением этому является отмеченная выше статистически значимая корреляционная зависимость между объемом ледникового стока рек и биомассой фитопланктона и зоопланктона. Отсутствие статистически значимой связи между многолетней динамикой фитопланктона и зоопланктона обусловлено временным запаздыванием изменчивости второго сообщества относительно первого с лаг-периодом 6–7 лет. Это свидетельствует, с одной стороны, о слабом влиянии зоопланктонного сообщества на количественные показатели фитопланктона, а с другой – о наличии внешнего фактора, вызывающего однонаправленные многолетние изменения планктонных сообществ. Зообентос оказывает более сильное контролирующее влияние на фитопланктонное сообщество, что проявляется в наличии статистической значимой отрицательной корреляционной зависимости между многолетней изменчивостью биомасс этих сообществ. Это связано с усилением доминирования в донном сообществе двустворчатых моллюсков, питающихся простейшими, органическими веществами, бактериями, фитопланктоном [21]. Доля Monodacna colorata в суммарной биомассе макрозообентоса возросла в последние десятилетия до 75,3–95,6 %. Перестройка структуры донного сообщества вызвана существенным сокращением промыслового стада бентосоядных рыб, в основном сазана, популяция которого представлена преимущественно особями младших возрастов. В результате недоиспользования сазаном моллюска как кормового объекта суммарная биомасса бентоса в 1996–2012 гг. (17,50 ± 3,37 г/м2) в среднем в 4 раза превысила величину показателя за предыдущие годы (3,53 ± 0,26 г/м2). Если биомасса макрозообентоса за счет моллюска Monodacna colorata сохранится на высоком уровне, следует ожидать через несколько лет снижения в среднем для озера количественных показателей зоопланктонного сообщества как конкурентного потребителя фитопланктона. С учетом того, что моллюск встречается только в Западном Балхаше, произойдет перераспределение обилия фитопланктона и, соответственно, зоопланктона в пространственном отношении.