Магнитные и электрические поля в биосфере вызывают различные биологические эффекты. Отчасти эти эффекты связаны с тем, что в живой ткани протекают электрические токи, вызванные биоэлектрогенезом клеток организма [1-3]. Они циркулируют в нервной системе, незначительные повреждения в которой могут повлечь за собой серьезные последствия [4,5].
В ряде работ, включая работы авторов, подчеркивалось влияние магнитных полей на проростки злаковых растений, хвойных деревьев и т.д. [4,6,7,8], на ориентацию птиц, пчел по геомагнитному полю, на поведение животных [4,9] и болезни людей [4,5]. Ранее нами численно решалась модель магнитного поля Земли с учетом разложения магнитного потенциала по нормированным присоединенным полиномам Лежандра в стационарных условиях [8].
Данную работу посвятим модели одной из аномалий, в частности экваториальной электроструе. Известно, что солнечные приливные движения и движения, обусловленные периодическим нагреванием атмосферы Солнцем, порождают магнитные вариации, которые называются спокойными солнечно-суточными вариациями ( ). Экспериментальные измерения магнитного поля показывают, что
вариации его составляющих
обусловлены в основном системой токов в ионосфере, которую можно связать с изменением потенциала
:
, (1)
где R - радиус Земли; - малые изменения северной, восточной и вертикальной составляющих напряженности магнитного поля Земли;
- коширота (дополнение широты φ до 90°); λ - долгота.
вариации на экваторе имеют очень большую величину по сравнению с другими областями. В узкой полосе (± 3¸5°) вдоль магнитного экватора сконцентрирован ток, названный экваториальной электроструёй.
ток индуцируется благодаря электродвижущей силе
, возникающей при движении проводящего слоя, которым является ионосфера. Так как вдоль магнитного экватора горизонтальная компонента динамо - ЭДС
равна 0, экваториальная электроструя образуется за счет электростатического поля
, направленного на восток, где
- потенциал, созданный пространственным зарядом. В атмосферном динамо полное поле
создает ток, величина и направление которого зависят от анизотропной проводимости атмосферы. Плотность тока можно определить по закону Ома.
Направленное на восток электростатическое поле приводит к возникновению тока Педерсена, также направленного на восток вдоль магнитного экватора и обусловленного проводимостью ионосферы перпендикулярной магнитному полю B. Значение проводимости Педерсена на высоте максимума слоя F2 много выше при высокой солнечной активности, чем при низкой. Как известно, дрейф заряженных частиц в скрещенных электрическом и магнитном полях в направлении
в присутствии столкновения ионов с нейтральными частицами вызывает ток Холла
. Токи в ионосфере можно определить [3,10]:
(2)
(3)
, (4)
где - проводимость нулевого поля, направленная вдоль магнитного поля.
Из выражения (3) . Тогда
![f](/use/i/2010/6/20.png)
где - проводимость Каулинга относится к току в направлении x, когда
.
На экваторе проводимость Каулинга на порядок больше, чем на других широтах. Ток Холла направлен вертикально вниз и поляризует тонкую оболочку проводящего слоя. Поляризация слоя приводит к появлению электрического поля, направленного вверх. Это поле вызывает ток Педерсена, направленный вертикально вверх и компенсирующий ток Холла. Вертикальное электрическое поле, кроме того, создает заметный Холловский ток в сторону востока, электроны при этом движутся на запад.
Интересно отметить, что и вариации на магнитном экваторе в равноденствие имеют максимум амплитуды. В равноденственные месяцы ветры атмосферы имеют симметричные компоненты. В остальные месяцы наблюдается сильная асимметрия в распределении ветровой системы и динамосистемы относительно экватора. Проводимость
вдоль геомагнитного поля велика и различие в распределении потенциала
может привести к возникновению электрических токов вдоль силовых линий из одного полушария в другое через экваториальную плоскость.
В работах [3,11] делается попытка построения модели экваториальной электроструи, в качестве основного параметра которой используется электрический потенциал. В отличие от ранее созданных [12], предлагаемая модель является нестационарной. В ней учитывается временная зависимость проводимостей и нейтральных ветров. Существенной особенностью модели является то, что электрическое поле, порождающее струю, зависит от высоты.
Рис. 1. Дипольная система координат, используемая в нестационарной модели экваториальной электроструи [11]
Делается предположение, что магнитное поле Земли является дипольным. Поэтому для описания электроструи используется дипольная система координат:
,
где - широта и долгота соответственно.
В основу модели положены следующие уравнения:
(5)
(6)
, (7)
где и
- компоненты электрического поля, параллельная и перпендикулярная магнитному полю
;
- скорость ветра;
- проводимость вдоль силовой линии магнитного поля;
- проводимости Педерсена и Холла;
- плотность тока; Ф - скалярный потенциал;
- единичный вектор в направлении
(рис. 1).
Поскольку проводимость вдоль силовой линии во много раз больше проводимостей Педерсена и Холла, то полагается, что . В силу этого должно выполняться условие:
.
В дипольной системе координат уравнения (5) - (7) имеют вид:
(8)
(9)
, (10)
где , коэффициенты Ламэ для данной системы координат.
Комбинируя (8) - (10), получаем уравнение для скалярного потенциала:
, (11)
где ;
;
- компоненты скорости ветра.
Предположим, что силовые линии являются эквипотенциальными, тогда можно освободиться от переменной φ, проинтегрировав уравнение (11) вдоль силовой линии от точки S0 - соответствующей границе проводящей области, до точки (рис. 1):
, (12)
где . Важным этапом для решения уравнения (12) является выбор граничных условий.
На границе проводящей области потенциал полагается равным нулю. На верхней границе дело обстоит иначе. Можно считать, что на больших Re потенциал электроструи ведет себя как потенциал точечного источника, т.е. . Такой потенциал удовлетворяет уравнению:
, (13)
которое и является верхним граничным условием.
Очевидно, что при это условие принимает еще более простой вид:
. Такой выбор верхнего граничного условия заставляет поднимать верхнюю границу достаточно высоко. Поэтому имеет смысл решать уравнение (12) с переменным шагом по Re. В связи с этим, необходимо преобразовать это уравнение, сделав замену:
; где
- текущая точка.
Соответственно .
Тогда уравнение (12) примет иной вид:
. (14)
Такая замена позволяет решать уравнение (13) с постоянным шагом по i. Следует отметить, что дифференцирование по Re и λ проводилось аналитически.
Таким образом, основным в предлагаемой модели является нахождение потенциала Ф из уравнения (13). Зная потенциал, мы можем по формуле (10) найти электрическое поле, а затем по формуле (8) - плотность тока.
Уравнение (14) решалось методом неполной матричной факторизации. Для этого уравнение заменялось разностной схемой. Метод решения и вопросы, связанные с устойчивостями разностной и итерационной схем, подробно изложены в работе [11]. Еще до получения решения уравнения (14) можно предсказать поведение потенциала Ф, исходя из следующих соображений. Поскольку по определению потенциала - это работа по перемещению единичного заряда, то для какого-то заряда, совершающего движение по силовой линии, можно записать:
, (15)
где - сила Лоренца,
. (16)
Отсюда, с ростом Re, потенциал должен возрастать, поскольку интеграл J при этом изменяется незначительно, а при стремится к константе.
В соответствии с этим, можно сказать, что решением уравнения (14) должна быть возрастающая функция. Поэтому в процессе решения задачи встал вопрос об изменении верхнего граничного условия. Его можно вывести следующим образом.
Продифференцируем уравнение (16) по Re:
. (17)
Второе слагаемое в (17) при дифференцировании будет всегда домножаться на выражение , т.е. при
,
. Первое слагаемое, в силу того, что сила Лоренца при
должна стремиться к нулю, также будет стремиться к нулю. Поэтому в качестве верхнего граничного условия можно выбрать следующее:
.
Рис. 2. Суточный ход потенциала Ф для разных Re, полученный
из нестационарной модели электроструи [11]
- - - вариации проводимости,
- × - изменение восточной компоненты скорости ветра
Рис. 3. Широтный профиль Е-поля для разных высот
Кроме того, уровень отсчета потенциала, то есть высота, на которой потенциал равен нулю, при расчетах выбирается не на границе проводящей области, а на высоте максимума проводимости, т.е. ~ 120 км. Как и предполагалось, с ростом Re потенциал возрастает сначала довольно быстро, а ближе к верхней границе выходит на постоянную [11].
На рис. 2 показано поведение потенциала с течением суток. В результате решения оказалось, что потенциал имеет наибольшее значение в 10 часов местного времени. Второе по величине значение потенциала наблюдается в 20 часов.
Как видно из рис. 3 изменение потенциала с увеличением Re становится всё более резким, что связано с зависимостью потенциала от Re . Поведение потенциала с течением времени несколько напоминает суточный ход проводимостей (пунктирная линия на рис. 2). Однако, в силу влияния ветров, которое является наибольшим в ночное время, имеются расхождения во временном ходе потенциала и проводимостей.
Прежде чем перейти к исследованию электрического поля, необходимо остановиться на методе его вычисления. Ясно, что взяв какую-то высоту (радиус R1), мы тем самым проходим несколько силовых линий ReT, а вместе с ними и некоторое число широт . Фактически для нахождения поля в точке R1, нами проектируется на этот уровень разность потенциалов
. Окончательно электрическое поле считается по формуле:
, (18)
где
,
. (19)
В результате удается получить зависимость поля от широт, которая рассчитывается по формуле (18). Как видно из рис. 3, с увеличением широты поле спадает, причем наиболее резкий спад наблюдается в области от 0° до 6°. Из этого же рисунка видно, что широтная зависимость поля с ростом высоты становится всё более слабой, что говорит о том, что поле с высотой падает, т.е. с ростом высоты влияние электроструи ослабевает. Как показали расчеты, на образование электроструи существенное влияние оказывает конфигурация магнитного поля Земли.
Высотный профиль электрического поля на экваторе изображен на рис. 4, где показана зависимость поля от высоты для разного времени. Такие профили могут быть получены, используя рис. 3. Для этого надо зафиксировать какую-либо широту и просмотреть все значения поля, соответствующие разным высотам.
Из рис. 4 видно, что на малых высотах поле изменяется резко, а затем с ростом высоты его изменение становится довольно слабым. Это можно связать с тем, что на этих высотах существенную роль играют проводимости Педерсена и Холла, а выше их влияние ослабевает. Что касается временного хода электрического поля, то можно сказать, что его изменение находится в противофазе с временным ходом потенциала.
Кроме потенциала и электрического поля в данной модели была рассчитана высотная составляющая плотности тока электроструи в районе экватора. На рис. 5 приведен высотный профиль плотности тока. На основе полученной зависимости можно сказать, что на формирование электроструи существенное влияние оказывают проводимости. И если своеобразное поведение силовых линий магнитного поля влияет существенным образом на локализацию струи по широте, то проводимости способствуют локализации электроструи в узком высотном интервале. Однако надо сказать, что на широтную локализацию струи проводимости также оказывают влияние, но для учета этого влияния необходимо знать достаточно точно изменение проводимостей с широтой.
Рис. 4. Высотный профиль электрического поля на экваторе для 10 и 12 часов местного времени
Рис. 5. Вариации восточной компоненты плотности тока электроструи на экваторе
с местным временем по результатам модели [11]
В данной модели широтная зависимость проводимостей не учитывалась. Полагалось, что проводимости ведут себя следующим образом:
, (20)
где - характерный масштаб изменения проводимостей;
;
- значения проводимостей, соответствующие их значению на уровне максимума. Временная зависимость проводимостей содержится в коэффициентах
.
На рис. 6 приведена зависимость восточной составляющей плотности тока от местного времени. Сравнение рассчитанных значений плотности тока электроструи с результатами экспериментальных измерений, показало, что чем конкретнее будет модель проводимостей, тем точнее окажутся результаты расчетов.
Рис. 6. Суточный ход восточной плотности тока для разных высот [11].
Кроме того, для построения точной модели электроструи необходима более реальная модель нейтральных ветров. В заключение следует сказать, что построенная модель экваториальной электроструи, по крайней мере, качественно описывает поведение данного объекта, и при замечаниях, отмеченных выше, она может дать результаты, находящиеся в хорошем согласии с экспериментальными данными.
Поведение параметров экваториального электроджета влияет на распределение заряженных [3] и незаряженных компонентов в этой области, на температуру и т.д., и подвержено не только широтным и высотным вариациям, но и изменениям, связанным с магнитными возмущениями общего геомагнитного поля Земли [13,14].
Согласно магнитоэкологии, наряду с универсальностью действия геомагнитного поля, являющегося единой пространственно-временной координатой для живых организмов, отмечается многообразие реакций биологических систем на электрические и магнитные поля. Экваториальные области отличаются особой специфичностью электромагнитной обстановки. Наша модель экваториального электроджета позволяет рассчитать не только электрические, но и магнитные поля на аномальных широтах. Наряду с моделью геомагнитного поля [5¸9] авторов, её можно использовать для моделирования биофизических задач, учитывающих влияние вариаций ГМП на живое.
СПИСОК ЛИТЕРАТУРЫ
- Рубин А.Б. Биофизика: в 2-х т. - Т.2: Биофизика клеточных процессов: Учебник. - 3-е изд. - М.: изд-во МГУ, 2004. - 469 с.
- Белановский А.С. Основы биофизики в ветеринарии: Учебное пособие. - 3-е изд. исправл. Часть I. - М.: МГАВМиБ им. К.И. Скрябина, 2002. - 136 с.
- Кутимская М.А. Жизнь с точки зрения биофизики. / Вестник Иркутского регионального отделения РАНВШ. - Иркутск, 2003. - № 2. - С. 122-128.
- Дубров А.П. Геомагнитное поле и жизнь. - Ленинград: Гидрометеоиздат, 1974. - 174 с.
- Кутимская М.А. Параметры магнитного поля Земли и их влияние на жизнедеятельность человека (синергетический подход) // Интеллект. и материальные ресурсы Сибири: Сб. науч. трудов. - Иркутск: ИГЭА, 2002. - С. 94-99.
- Кутимская М.А. Параметры магнитного поля Земли и их влияние на жизнедеятельность растений (синергетический подход). // Юбилейный сборник научных трудов: Электрификация и автоматизация АПК в условиях Восточной Сибири. - Иркутск: ИрГСХА, 2001. - С. 65-70.
- Kutimskaya M.A. The earth´s magnetic field and plants vital activity. International Scientific Conference. Agricultural and applied Sciences in the development of farming and forestry: actual problems, practice and exchange of experience. - Irkutsk, June 6-11, 2006. - p. 248 - 253.
- Kutimskaya M.A., Jozefaciuk G, Wrraszcz E, Buzunova M.U. Effect of magnetic fields and plants vital activity. /Physics in agricultural research. International Scientific Conference. Papers and short communication. - June 12-13, 2008, Lublin, Poland.
- Кутимская М.А., Волянюк Е.Н., Убрятова Л.В. Параметры магнитного поля Земли и их влияние на жизнедеятельность животных (информационно-синергетический подход) // Материалы региональной НПК «Актуальные проблемы АПК». - Иркутск: ИрГСХА, 2002. - С. 4-5.
- Кутимская М.А., Фидель Б.В., Сидоров И.М. Экваториальная электроструя и вертикальный дрейф заряженных частиц. - Иркутск, 1978. - 34 с. - Деп. в ВИНИТИ № 595.
- Фидель Б.В., Кутимская М.А., Коен М.А. Модель экваториальной электроструи с учетом временных вариаций. - Иркутск: ИГУ, 1979. - 31 с. - Деп. в ВИНИТИ № 1986.
- Фидель Б.В., Кутимская М.А. Модель экваториальной электроструи. - Иркутск: ИГУ, 1978. - 56 с. - Деп. в ВИНИТИ № 594.
- Акасофу С.И., Чемпен С. Солнечно-земная физика. Пер с англ. - М.: Мир, 1974. - Ч. 1. - 284 с. - 1975. - Ч. 2. - 512 с.
- Паркинсон У. Введение в геомагнетизм: Пер. с англ. - М.: Мир, 1986. - 528 с.