чается в литературе. Труднее определить значение В. При отсутст) и информации можно, в первом приближении, принять Р р: нулю, в результате чего $а СТаНОВИТСя кОнсантой. Теперь раженне для стокового члена запишется в виде 5— «(#) 5к (2). Просуммированный по всем слоям объем. поглощенной Воды не может превышать потенциальную скорость транспирации, т. ©- 5аг <Т ля |2|< & (68) 'Поскольку в данной программе извлечение влаги рассматривается от верхнего слоя почвы к нижележащим, это выражение позволяет: моделнровать извлечение влаги. предпочтительно з верхних, от- носнтельно влажных слоев. Потребность в потенциальной транс- пирации может удовлетворяться о тех пор, пока растение не из- влечет влагу из всех слоев корнеобитаемой Зоны почеы. Объединив уравнение (6!), закон Дарси и уравнение (63), ‚а также введя днфференциальную влагоемкость С =@0/48, можно получнть дифференцнальное уравненне движения воды в част- ных производных, которое следующим. образом. описывает поток воды в снстеме почва-_-корин: оу (К(& + )-ев — © тде $(#) определяется согласно уравнению (65). Дая решення ‘уравнения (69) необходнмо задать подходящие начальные н гра- ничные условия. В качестве начального условия (#=0) задается давленнё почвенной влаги как функция вертикальной. коорди- наты : ле (=0) —0 о) На нижней граннце (—Ё) давление задается как ве=-1. )—^( со Поток воды (поток Дарск) @ на верхней граннце почвы зави- снт от метеорологических условнй. Почва может терять воду в про- цессе нспарения нли получать её в процессе инфнльтрации. Хотя максимально возможная (потенциальная) скорость нспарения з данной почвы зависнт ншь от атмосферных условий, действнтель- ный поток через поверхность почвы ограничивается ‘способностью пористой среды пропускать воду снизу вверх. Аналогично этому, если потенциальная скорость инфильтрации (например, интенси ность выпадения осадков или орошения) превышает поглощающую. способность почвы, часть воды будет теряться в результате о- верхностного стока. И здесь также потенциальная скорость нн- Фильтрации. регулируется атмосферными (или другимн) внеш- ними условиями, тогда как фактическая инфильтрация зависнт 216 от предшествующего увлажнения почвы. Таким образом, точное траничное условие для поверхности почвы не известно априо решение следует искать в результате максимизации абсолюткого Значения потока испаряемой влаги (НапКз е а/,, 1969).. 'Если принять 9° (2=0, #) в качестве максимально. возможной: скорости испарения © повёрхности почвы, необходимо, чтобы удов- ‘летворялись следующие условия: 1г*е=0, 0|21ое=0, д1 =| -К( () @Э) три в< () где No — минимальное давление почвенной влаги, которое возмо- жно при воздушно-сухом состоянии почвы. Если допустить, что давление почвенной влаги на поверхности почвы находится в рав- новесии © атмосферой, то можно определить Й, используя. извест- ное соотношение: #г л— () 4) тде & — универсальная газовая постоянная, Дж-моль“!: КЭ!; — абсолютная температура, К; д— ускорение свободного падения, м-с-ё; М — молекулярная масса воды, кг: моль-!; Р — относител ная влажность воздуха, доли. 'Значения Т.и Р можно получить по измеренням в метеородоги- ческой. будке. `Существуют различные способы определения максимально воз можного испарения из почвы д° (уравненне (72)), как и максн мально возможной скорости транспирации Тн которая определяет максимально возможное поглощение воды корнями на единнцу площади почвы (ем. уравнение (65)). Потенциальная эвапотранс- пирация ЕТ* является суммой потенцнальной транспнрации Тн и потенциального испарення из почвы (испаряемости) Ё*: вТ = Та Е*. (75) Значение ЕТ* можно рассчитать, например, по комбинированному уравнению энергетического баланса Монтейта-—Рийтемы, по ура: нению Пристли и Тейлора, илн умножая значение испарёния с от- крытой водной поверхности (по Пенману) на транспирационный коэффициент сельскохозяйственной культуры. Значения Е* можно. рассчитать по комбинированному уравнению энергетического 6% ланса, опустив аэродинамический член и учнтывая только ту долю ® которая достнгает поверхности (метод Ричи). Отсюда можно: опредедть Т.а как величину, которая остается неизвестной в ура! ненни (75). Волее подробные данные содержатся в работе Фед- деса н др. (Еедаез её а), 1978). "Теперь, когда известны начальные н другие граничные усло- вия, можно решить уравненне (69), аппроксимируя его неявной конечно-разностной схемой и задавая соответствующую сетку неза- Эи