МОСТЫ МЕЖДУ СИСТЕМАМИ. ПЛОДОРОДИЕ, ОКЕАН КАК ПОЧВА
. … Четыре фактора можно считать первичными
регуляторами природы: энергия, вещества, услот. вия существования и сообщество. [37]
Е. Одум
‘ ‘ * і. •
Мы. представляем в общем виде схему взаимодействия природных процессов, действующих на поверхности Земли, в океане и атмосфере, но • наши сведения о количественных соотношениях
в этих связях и в особенности о связях типа управления пока еще отрывочны и примитивны. Таким образом, мы воздействуем на природу, но. в то же время еще плохо понимаем, каковы
. будут проистекающие ог этого следствия, ii, как
? правило, не считаемся с ними.[38]
, Е. Федоров
Определим самую примитивную зависимость насыщения океанских вод питательными солями от гидродинамических процессов (перемешивания воды и ее вертикальной скорости) с учетом глубины, на которой концентрация питательных солей не меняется по сезонам (это может быть дно или слой скачка плотности).
Если вырезать в воде какой-то неподвижный объем пространства и предположить, что концентрация питательных солей меняется в нем за счет перемешивания воды, вертикального движения и жизнедеятельности планктона, то уравнение баланса питательной соли С можно записать так: ,
дС т, д2С дС „п
■Kz л 2 uz, аС, (105)
dt z dz2 г дг
где — суммарное изменение концентрации соли в
u*(J ‘
объеме; Kz изменение концентрации, вызываемое вертикальным перемешиванием ВОДЫ (Кг— коэффициент диффузии, показывающий интенсивность переме — . дС
шивания)uz — qj — .изменение, вызываемое вертикальным движением воды (их— вертикальная составляющая скорости течения); аС — изменение концентрации, вызываемое жизнедеятельностью планктона и пропорциональное количеству питательной соли.
Допустим, что суммарное изменение концентрации целиком вызывается жизнедеятельностью планктона. Тогда уравнение (105) распадается на два:
<РС иг dC |
(106) (107) |
dz2 Kz dz
Решение уравнения (106) уже подробно рассматривалось. Решение же уравнения (107) будет зависеть от величины соотношения, называемого числом Пекле: 1
где Я — расстояние до некоторого источника питательных солей, т. е. до того горизонта, где концентрация питательных солей постоянна, например расстояние до дна или до слоя глубинных вод.
Вид соотношения связан с тем, что при вертикальном движении питательные соли переносятся по оси z не разбавляясь, а при перемешивании разбавляются пропорционально удалению от источника. Поэтому нужно искать не простое соотношение двух величин Uz и Kz, а соотношение, взятое с учетом влияния этих величин на перенос питательных веществ. Учет достигается делением Kz на глубину места (или расстояние до источников солей) Я, ‘
Решением уравнения (107) будет следующая формула: .
где С*— концентрация питательной соли в безразмерных единицах.
Скорость течения, направленную вверх, будем считать положительной, вниз — отрицательной. В соответствии с этим и число Пекле будет иметь или положительные, или отрицательные значения. На рис. 46 показаны распределения концентрации при числах Пекле» равных 16*г 8; 4; 2; 1; 0; —1; —2; —4; —8; —16.:
Рис. 46. Структура океанической «почвы» в зависимости- от числа Пекле (соотношения вертикальной скорости, интенсивности, перемешивания и расстояния до источ — ‘ ника). Сплошные линии — вертикальное распределение питательных солей в богатых (дивергентных) шочвах», пунктирные — в бедных (конвергентных) «почвах». Цифрами показаны числа Пекле (Ре = и: К : Н). . |
Рисунок построен таким образом, что в него можно подставить любые значения С, Kz, uz, Н, z, так как вместо расстояния до источника в метрах по оси ординат отложено безразмерное расстояние в долях от величины Н, принятой за единицу. По оси абсцисс вместо абсолютной концентрации в миллиграммах на 1 см3 отложена безразмерная концентрация, выраженная в долях
Положение круговоротов Океаничес — и среднее содержание кие
Площадь круговоротов, млн. км[39] |
зоопланктона, т/км2 фронты
Рис. 47. Конвергентные и дивергентные «почвы» океана. |
единицы. При этом за единицу взята разность между высокой и постоянной концентрацией на глубине Я и пониженной, но тоже постоянной в данный сезон концентрацией на поверхности океана.
Из рисунка видно, что плодородие океанской почвы, насыщенность ее питательными солями, существенно зависит от направления вертикальной скорости.
Вертикальные движения связаны со сгонами и нагонами или дивергенциями (расхождениями) и конвергенциями (схождениями) морских вод, направлением ветрового течения и положением берега. Суть здесь заключается в том, что течение, идущее в антициклоническом направлении (т. е. по часовой стрелке в северном полушарии), имеет циклоническую завихренность, или ротор скорости, в прибрежной зоне от стрежня потока к берегу.
В то же время в центральной части круговорота завихренность скорости носит антициклонический характер. В циклоническом круговороте — все наоборот. .
Циклоническая завихренность означает увеличение скорости слева направо, если смотреть вдоль потока, — и это благодаря силе Кориолиса приводит к сгону (дивергенции) вод. Антициклоническая завихренность означает уменьшение скорости — слева направо и, следовательно, при ветровом течении произойдет нагон (конвергенция) вод за счет той же силы Кориолиса.
На рис. 47 показаны зоны дивергенции и конвергенции. Проставленные на них цифры показывают среднее содержание в воде планктонных организмов в тоннах сырого веса на 1 км2 поверхности океана. Цифры относятся к верхнему стометровому слою моря. Они осред — няют многочисленные наблюдения на научных судах (эти осредненные расчеты проведены в МГУ И. А. Суе — товой). При сравнении рис. 46 и 47 видно, что между гидродинамикой вод, балансом биогенных элементов и биологической продуктивностью океана существует тесная связь и что на уравнение (105) можно смотреть как на один из самых простых мостов,-соединяющих гидродинамические, химические и биологические системы океана..
Уравнение (105) значительно упрощает развитие событий в реальном океане. Действительное положение вещей полнее отразила бы такая запись:
дС „ д2С дС
— Z~dz5 М2 дг " ЬСц, (110)
где Со •— концентрация питательных элементов, заключенных в детрите в 1 л морской воды, (примем Со = = const); ЬСо — скорость их перехода из детрита в морскую воду, предполагаемая постоянной по глубине.
При условиях стационарности уравнение (110) можно записать в виде уравнения ■
«С=»с,. (Ш)
Если учесть, что концентрация детрита меняется не
сильно, то общим решением этого уравнения будет:
С— + C2er*z + В, (112)
где В, Ci, C2 — постоянные, а гь. г2 — корни алгебраи-.
ческого характеристического уравнения.
КхГ2 ■— иг — а = 0,
получающегося заменой С и производных на степени г
, „ ndC ■ , d2C
(т. е. С на г на г и на г2).
Или еще лучше вместо уравнений (110) и (111) записать выражение более общего характера:
Ъ^~а*-1&—аС = р’- (113>
где F — любая функция, аргументом которой не является С, отражающая, вообще говоря, любое воздействие на систему на различной глубине.
Решение уравнения (113) можно получить классическим путем, аналогично предыдущему, только В, в зависимости от F, может оказаться и не постоянной, например: .
В = В0 cos юjZ — f — 5j sin WfZ при F=bsin ауг.
Можно получить решение и операторным — методом» очень популярным среди специалистов по автоматиче-
скому регулированию. Покажем, не решая уравнение до? конца, почему этот метод мог бы оказаться полезным: при изучении интересующей нас системы.
Если воспользоваться символом 5 для операции дифференцирования^^ ~ , S2= ~^£г> то
aL |
с |
F. |
■О |
Первый сомножитель — передаточная функция системы, а второй— входная функция F. Передаточна»
Рие. 48. Блок-схема решения уравнения (113) на аналоговой вычислительной машине.
функция состоит из сомножителя I/а, определяющего) коэффициент усиления системы, и сомножителя 1/L, определяющего пространственно-динамические свойства системы. Таким образом, операторный метод не только’ удобен, но и нагляден. Это уравнение легко решить на. аналоговой машине (рис. 48).
Но и этот мост между системами весьма примитивен. Действительные связи между физическими, химическими и. биологическими системами многообразны. При увеличении числа отображаемых в модели элементов системы и связей между ними порядок дифференциальных уравнений повышается, они могут стать нелинейными уравнениями, аналитическое решение которых получить обычно невозможно.
Некоторое представление о сложности реальной системы дает рис. 49. Здесь учитываются средние концентрации в продуктивном слое моря органического вещества, оседающего. детрита, питательных солей (фосфаты и соединения азота) и биомассы двух основных видов фитопланктона и бактерий, а также влияние света, температуры и вертикального переноса питательных солей из глубоководной части, где их концентрация при-
Рис. 49. Блок-схема связей в системе — фитопланктон — органические вещества — бактерии —. биогенные элементы —. /-г-освещенность; Т — температура; К2 — интенсивность перемешивания; В,, В2 — диатомовые водоросли; В3 — водоросли динофлагелляты; В4 — бактерии; В5 — зоопланктон, ВОВ — детрит; РОВ — раст-. воренное органическое вещество; Р—фос — фор; N — азот. |
нимается постоянной. На цифровой — вычислительной машине рассчитывались сезонные изменения в системе, и ЭВМ неплохо имитировала реальные события (рис. 50). Необходимо подчеркнуть, что действительные связи в системе еще многообразнее и сложнее. Но можно ли до бесконечности усложнять модель? Она станет необозримой и мало полезной. Можно перечислить «бесконечное» число параметров, которые и физикам, и химикам, и биологам важно учесть. А неизбежные случайные воздействия извне и в самой системе? Очевидно, что на помощь должны прийти статистические и вероятностные модели. Рассматривавшиеся в этой книге модели, основанные на дифференциальных уравнениях (их называют еще «жеет-
Рис. 50. Сезонные изменения в морской экологической системе (результаты моделирования на ЭВМ). Уел. обозначения см. рис. 49. |
кими» и «детерминистическими»), позволяют сделать немалое—привести в порядок имеющиеся сведения и проверить, как они согласуются между собой и не противоречат ли действительности. Но это только начало пути.