Т' и 7° если интенсивность касательных напряжений То в не- которой области на ложе достигает величины 1‘1. скорость И) резко увеличивается. Прежний режим скольжения может восста- новиться, если Ту, уменьшаясь, достигнет критического значения < 72 Физически интерпретировать такой механизм движения льда у ложа можно следующим образом: при сравнительно малых напряжениях 1‘„‹1‘2 происходит преимущественно пластическое обтекание характерных неровностей подстилающей поверхности. При То = Т° достигается „предел прочности" льда в придонном слое и образуются внутренние поверхности скола-скольжения. В результате этого эффективная шероховатость ложа значительно уменьшается, что приводит к резкому увеличению скорости проскальзывания No). Сплошность льда у ложа начинает восста- навливаться при падении скорости относительного движения образовавшихся поверхностей скола до величины, соответствующей нижнему критическому значению Т'. Модель (5.1) — (5.4) опи- сывает автоколебания ледников и объясняет появление и распространение кинематических волн [92]. Алгоритмы численного решения задач (5.1) — (5.2), (5.1) — (5.4) рассматривались в 191, 92] и основаны на применснии метода сеток. Подготовка исходных данных и предварительные ` расчеты. Для построения функций го, Н,Б и { были использованы ма- териалы многолетних исследований на леднике. Поскольку регулярные наблюдения за балансом массы и поверхностной скоростью ледника проводились лишь с 1967/68 г., он был выбран в качестве базового при подготовке исходных данных. По карте ледника 1967 г. с учетом распределения площадей высотных Ззон определены осредненный продольный профиль поверхности ледника % (#) и его ширина Н (х). Для построения функции 6 (х) исполь- зованы данные непосредственных наблюдений за распределением удельного баланса массы по высотным зонам. Продольный профиль ложа 29(х) предварительно восстанавливался — по результатам радиозондирования языка ледника и области питания в отдельных точках. Чтобы адекватно подобрать значения параметров модели (5.1) и оценить влияние неточностей при задании исходной ин- формации, были проведены расчеты, в которых эксперимснтально определенные функции Н, г), 6 и парамстры модели (5.1) — (5.2) & и л варьировались в разумных пределах для достижения на- илучшего согласования измеренного /) и рассчитываемого по дан- ным о балансе массы 1967 г. профиля / поверхности ледника, Расчеты по стационарной модели (5.1) — (5.2) при д//д: = 0 показали, что толщина ледника в основной части ус- 143