СПОСОБ ОПРЕДЕЛЕНИЯ ВЛАГОЗАПАСА В СНЕЖНОМ ПОКРОВЕ
METHOD FOR DETERMINING MOISTURE RESERVE IN SNOW COVER

Изобретение относится к метеорологии и может быть использовано для определения влагозапаса (запасенной за зиму влаги) в снежном покрове по измеренной плотности потока бета-излучения в приземной атмосфере.
Известен способ определения влагозапаса в снежном покрове путем измерения водного эквивалента снежного покрова [Снег. Справочник. / Под редакцией Д.М. Грея и Д.X. Мэйла. - Л.: Гидрометеоиздат, 1986. - С. 206.]. Способ заключается в измерении массы образцов снега, извлеченных из снежного покрова с помощью пробоотборника весового снегомера. Пробоотборник весового снегомера [Лесная метеорология. Метеорологические приборы и наблюдения: учебное пособие / А.В. Иванов. - Йошкар-Ола: Марийский государственный технический университет, 2009. - С. 100.] представляет собой полый цилиндр. Для проведения измерения цилиндр вертикально погружают в снег до контакта с почвой. Цилиндр, содержащий снег, вынимают и взвешивают на специальных весах, входящих в комплект снегомера. Плотность снега рассчитывают из полученного веса и объема пробы. Весовым снегомером можно сразу определить влагозапас в снежном покрове в метрах, а также рассчитать, сколько воды в тоннах или кубических метрах получит 1 га от таяния снега.
Однако способ трудоемок. Внутреннюю поверхность пробоотборника трудно очистить от снега, влаги, грунта, остатков растительности, что снижает достоверность результатов измерений. При проведении механических измерений зачастую теряются отдельные части прибора.
Известен радиометрический способ определения влагозапаса в снежном покрове [Москвич Т.И. Микроволновое зондирование как метод оценки снегозапасов // Труды ГУ Дальневосточный региональный научно-исследовательский гидрометеорологический институт - 2009. - №153. - С. 25.], основанный на пассивном СВЧ зондировании, заключающийся в определении высоты снежного покрова и снегозапаса по измеренным данным с искусственных спутников Земли в микроволновой области спектра (λ=0,1÷30 см) посредством пассивной радиометрической системы. Способ позволяет производить картирование областей залегания снега, выявлять состояние снежного покрова и, таким образом, выделять границы площадей одновременного снеготаяния. Физическая основа распознавания снежного и ледяного покровов заключается в широком диапазоне изменения их коэффициентов излучения в санти- и миллиметровом диапазонах, при этом значительное влияние оказывает дополнительное рассеяние, которое создает снежный покров на поверхности почвы. Возможность измерений характеристик снежного покрова в микроволновом диапазоне обусловлена большой разницей значений диэлектрической постоянной сухого и влажного снега. Влагозапас в снежном покрове определяют из снегозапаса параллельно с оценкой высоты снега.
Использование способа требует организации дополнительных трассовых измерений толщины снежного покрова или проведения измерений в репрезентативных реперных точках, что усложняет процесс оценки влагозапаса.
Известен способ определения влагозапаса в снежном покрове [Loijens H.S. Measurements of snow water equivalent and soil moisture by natural gamma radiation. // Proc. Can. Hydrol. Symp., 1975, Aug. 11-14, Winnipeg. P. 43-50], выбранный в качестве прототипа, основанный на зависимости гамма-излучения в атмосфере от влагозапаса в снежном покрове. Способ заключается в измерении числа фотонов гамма-излучения при полетах самолета на высоте 20-100 метров по сети маршрутов. Изменение влагозапаса в снежном покрове за время между двумя полетами определяют из выражения
где ΔР - изменение влагозапаса в снежном покрове между двумя полетами;
1
2
1
2
I и I - число фотонов гамма-излучения с энергией в интервале от E до Е в первом и втором полетах, соответственно;
α - коэффициент, характеризующий спектральные характеристики излучения.
Сложность этого способа заключается в необходимости использования авиации.
Предложенный способ определения влагозапаса в снежном покрове расширяет арсенал средств аналогичного назначения.
Способ определения влагозапаса в снежном покрове, также как в прототипе, включает измерение ионизирующего излучения над снежным покровом.
Согласно изобретению производят измерение плотности потока бета-излучения в период перед началом таяния снежного покрова, в дневное время суток, и не менее чем через 3,5 часа после выпадения снега, а влагозапас в снежном покрове определяют из выражения:
где: W - влагозапас в снежном покрове, м;
2
- измеренное на высоте z от земной поверхности значение плотности потока бета-излучения при наличии снежного покрова, 1/(м⋅с);
2
- плотность потока вторичного бета-излучения на высоте z от земной поверхности в отсутствии снежного покрова, создаваемого содержащимися в поверхностном слое грунта радионуклидами 1/(м⋅с);
в
3
ρ - плотность воды, кг/м;
m
2
μ - массовый коэффициент ослабления излучения, м/кг.
Известно, что в период года без снежного покрова суммарная плотность потока бета-излучения в приземной атмосфере обусловлена 4-мя составляющими: 1) бета-излучением радионуклидов, содержащихся в грунте; 2) вторичным бета-излучением, которое образуется при взаимодействии гамма-излучения почвенных радионуклидов с атмосферой [B.C. Яковлева, В.Д. Каратаев, В.В. Зукау. Моделирование атмосферных полей γ- и β-излучений, формирующихся почвенными радионуклидами // Вестник КРАУНЦ. Физ.-мат. науки. 2011. №1 (2). С. 65-74]; 3) вторичным бета-излучением, образованным при взаимодействии космической радиации с атмосферой; 4) бета-излучением короткоживущих дочерних продуктов распада изотопов радона, содержащихся в атмосфере [Яковлева B.C. Методы и приборы контроля полей α-, β-, γ-излучений и радона в системе "грунт-атмосфера": диссертация на соискание ученой степени доктора технических наук: спец. 05.11.13; Национальный исследовательский Томский политехнический университет. - Томск: 2013, С. 64].
Так как в зимний период снежный покров толщиной более 0,1 м полностью поглощает бета-излучение от радионуклидов, содержащихся в грунте [Машкович В.П., Кудрявцева А.В. Защита от ионизирующих излучений: Справочник, 4-е изд., перераб. и доп. - М.: Энергоатомиздат, 1995. С. 423], то первую составляющую не надо учитывать. Вклад космической радиации в суммарный бета-фон приземной атмосферы на высотах до 50 м от земной поверхности составляет около 1% [Яковлева B.C., Каратаев В.Д., Зукау В.В. Моделирование атмосферных полей γ- и β-излучений, формирующихся почвенными радионуклидами. Вестник КРАУНЦ. Физ.-мат. науки, 2011, №1 (2), С. 64-73], следовательно, третьей составляющей можно пренебречь.
Количественной характеристикой интенсивности турбулентного перемешивания в приземной атмосфере является коэффициент турбулентности [Гальперин Б.М. Турбулентный тепло- и влагообмен поверхности суши и водоемов с атмосферой: Учебное пособие. - Ленинград: Изд-во ЛГИ, 1970. С. 11]. При низком коэффициенте турбулентности бета-излучение короткоживущих дочерних продуктов распада радона, содержащихся в атмосфере, может вносить значимый (до 30%) вклад в суммарный бета-фон [Яковлева B.C. Методы и приборы контроля полей α-, β-, γ-излучений и радона в системе "грунт-атмосфера": диссертация на соискание ученой степени доктора технических наук: спец. 05.11.13; Национальный исследовательский Томский политехнический университет. - Томск: 2013. С. 64]. Проведение измерений в то время суток, когда интенсивность турбулентного перемешивания в приземной атмосферы максимальна, позволяет уменьшить этот вклад до единиц %. Максимум коэффициента турбулентности наблюдается в дневные часы [Гальперин Б.М. Турбулентный тепло- и влагообмен поверхности суши и водоемов с атмосферой: Учебное пособие. - Ленинград: Изд-во ЛГИ, 1970. С. 27]. Так как зимой плотность потока радона из грунта в атмосферу существенно снижается, в среднем, в 2-3 раза [Яковлева B.C. Методы измерения плотности потока радона и торона с поверхности пористых материалов: монография. - Томск: Изд-во ТПУ, 2011. С. 120, С. 144], то это приводит к пропорциональному снижению среднего содержания радона и дочерних бета-, гамма-излучающих продуктов его распада в атмосфере, что позволяет пренебречь вкладом четвертой составляющей суммарной плотности потока бета-излучения в приземной атмосфере.
Поскольку при выпадении атмосферных осадков бета-излучение осажденных на земную поверхность дочерних продуктов распада изотопов радона сильно увеличивает бета-фон приземной атмосферы в 2 и более раз [Yakovleva V.S., Nagorsky P.M., Cherepnev M.S., Kondratyeva A.G., Ryabkina K.S., 2016. Effect of precipitation on the background levels of the atmospheric β- and γ-radiation. Applied Radiation and Isotopes. 118, 190-195], то необходимо выждать не менее 3,5 часов (около 7 периодов полураспада) до полного распада осажденных радионуклидов, чтобы избежать их влияния на результат измерения .
Выбор высоты измерения плотности потока бета-излучения ограничен только высотой снежного покрова и техническими возможностями.
m
238
232
40
137
Определение плотности потока вторичного бета-излучения  на высоте z от земной поверхности в отсутствии снежного покрова и массового коэффициента ослабления излучения μ не представляет особой сложности, поскольку необходимые данные по удельным активностям радионуклидов U, Th, K и Cs содержатся в научной литературе и в отчетах о геологических изысканиях.
Таким образом, предложенный способ определения влагозапаса в снежном покрове является простым, поскольку достаточно произвести только одно измерение плотности потока бета-излучения  при наличии снежного покрова.
Способ пригоден как для разовых измерений, так и для мониторинга влагозапаса в снежном покрове, если используется блок детектирования бета-излучения, работающий в режиме мониторинга.
На фиг. 1 представлены временные изменения: а) атмосферной температуры по данным с сайта ООО «Расписание Погоды», данные усреднены по 12 часов; б) плотности потока бета-излучения на высотах 1 м и 10 м, полученные экспериментально с 1 октября 2016 г. по 1 октября 2017 г., данные усреднены по 12 часов; в) влагозапаса в снежном покрове по данным с сайта ООО «Расписание Погоды», г) высоты снежного покрова по данным с сайта ООО «Расписание Погоды».
Для определения влагозапаса в снежном покрове выбрали площадку, расположенную недалеко от института мониторинга климатических и экологических систем СО РАН г. Томска.
По результатам анализа динамики атмосферной температуры по данным с сайта ООО «Расписание Погоды» [RL: http://rp5.ru/ Погода_в_Томске], с учетом метеорологических прогнозов, определили, что после 13 марта 2017 г. начнется таяние снега. Это подтверждается переходом среднесуточной атмосферной температуры (а на фиг. 1) из области отрицательных значений в область положительных значений.
Максимальная турбулентность атмосферы наблюдается, в основном, в дневные часы [Гальперин Б.М. Турбулентный тепло- и влагообмен поверхности суши и водоемов с атмосферой: Учебное пособие. - Ленинград: Изд-во ЛГИ, 1970. - С. 27], поэтому выбрали время начала проведения измерения в 14:00 час.
В 14:00 13 марта 2017 г. в отсутствие снежных осадков произвели измерение плотности потока бета-излучения на высоте 1 м от земной поверхности, с использованием блока детектирования бета-излучения БДПБ-01, которое составило  (б на фиг 1).
Значение плотности потока вторичного бета-излучения  на высоте z=1 м от земной поверхности в отсутствии снежного покрова и значение массового коэффициента ослабления излучения в снежном покрове рассчитали с помощью метода Монте-Карло в программе PCLab [Компьютерная лаборатория (КЛ/PCLab). Свидетельство о государственной регистрации программы ЭВМ №2007615275 от 28.12.2007]. В качестве входных данных использовали данные об удельной активности радионуклидов, содержащихся в поверхностном слое грунта, для территории Томской области [Каратаев В.Д., Яковлева B.C., Эргашев Д.Э. Исследование радиоактивности объектов окружающей среды на территории Томской области // Известия вузов. Физика. 2000. Т. 43. №4. С. 105-109] следующим образом: определили плотность потока вторичного бета-излучения  в атмосфере, создаваемое на высоте 1 м от земной поверхности почвенными бета-, гамма-излучающими радионуклидами в отсутствие снежного покрова при W=0, задавая в качестве источника ионизирующего излучения грунт с радионуклидами, а в качестве поглотителя - атмосфера; рассчитали зависимость плотности потока вторичного бета-излучения в атмосфере, создаваемого почвенными бета-, гамма-излучающими радионуклидами на высоте 1 м от земной поверхности, от влагозапаса  при наличии снежного покрова, задавая в качестве источника ионизирующего излучения - грунт с радионуклидами; в качестве поглотителей - снежный покров и атмосфера. Снег представляли как водный эквивалент (влагозапас). Значения влагозапаса задавали от 0 м до 0,5 м с шагом 0,05-0,1 м. Полученную зависимость  на высоте z=1 м аппроксимировали экспоненциальной функцией вида:
m
Массовый коэффициент ослабления излучения μ определили из выражения (2). При моделировании произвели учет вторичного излучения, каскадного характера взаимодействия излучения с атмосферой, процессов альбедо, а также учет векового радиоактивного равновесия между родоначальниками рядов и дочерними продуктами их распада.
226
238
232
40
137
Удельные активности Ra (U), Th, K и Cs в поверхностном слое грунта составляют 25; 26; 345 и 10 Бк/кг, соответственно [Каратаев В.Д., Яковлева B.C., Эргашев Д.Э. Исследование радиоактивности объектов окружающей среды на территории Томской области // Известия вузов. Физика. 2000. Т. 43. №4. С. 105-109].
235
-3
2
3
U235
U238
m
m
в
Активность U определили по известному вековому соотношению A=A/21. Рассчитанное значение плотности потока вторичного бета-излучения на высоте 1 м от земной поверхности в отсутствии снежного покрова составило  Массовый коэффициент ослабления излучения μ определили из выражения (2), который составил μ=2,17⋅10 м/кг. Плотность воды ρ=1000 кг/м.
Влагозапас в снежном покрове из выражения (1) составил:
Полученное значение влагозапаса в снежном покрове хорошо согласуется с данными сайта ООО «Расписание Погоды» (в на фиг. 1), что свидетельствует о достоверности предлагаемого способа определения влагозапаса в снежном покрове.
Полученная зависимость плотности потока бета-излучения в приземной атмосфере в течение зимнего периода (б на фиг. 1) снижается практически экспоненциально с ростом влагозапаса (водного эквивалента снега) (в на фиг. 1) и, соответственно, высоты снежного покрова (г на фиг 1). Снижение плотности потока бета-излучения наблюдается вплоть до начала таяния снега (б фиг. 1). Зависимости изменения плотности потока бета-излучения в приземной атмосфере на разных высотах 1 м (кривая 1) и 10 м (кривая 2) от изменения влагозапаса аналогичны, (б на фиг. 1). После полного схода снежного покрова плотность потока бета-излучения возрастает до своего значения, наблюдаемого до начала установления снежного покрова.
Всплески бета-фона в 1,5-3 раза, наблюдаемые на фиг. 1 б в разные сезоны года, обусловлены выпадением дождя (при положительной атмосферной температуре (а на фиг. 1) или снега (при отрицательной атмосферной температуре (г на фиг. 1).