Аннотация результатов, полученных в 2021 году
Основная научная проблема, на решение которой направлен проект - усовершенствование представления пограничного слоя атмосферы в моделях прогноза погоды и климата с учетом локальных неоднородностей подстилающей поверхности. С этой целью в проекте решается задача определения физических механизмов формирования структуры и динамики атмосферной турбулентности и волновых процессов в условиях пространственно-неоднородной поверхности и разработки параметризаций для расчета характеристик энерго- и массообмена между атмосферой и поверхностью в условиях различных типов поверхностных неоднородностей. На первом этапе проекта основное внимание было направлено на организацию измерительных полигонов и разработку методик анализа и обработки экспериментальных данных, а также разработку новых подходов к моделированию турбулентности над неоднородным ландшафтом. 1.1. Проведена модернизация измерительного комплекса для исследования атмосферной турбулентности в неоднородном ландшафте на базе метеорологической обсерватории МГУ. Разработана система сбора и первичной обработки данных. Пакет программ для сбора, подготовки и отправки данных написан силами участников проекта на языке Python 3.7. Исходный код пакета находится в открытом доступе по ссылке https://github.com/alexavr/tower_parse. Вся первичная информация, поступившая за прошедшие 24 часа, отображается на веб-сайте https://tower.ocean.ru. По окончании каждых суток буферный массив с данными микрометеорологической мачты конвертируется в формат NetCDF (www.unidata.ucar.edu/software/netcdf/). В качестве глобальных переменных выступают название мачты, ее описание, географические координаты, и описание оборудования: код, модель прибора, высота установки, частота измерений, дата установки. Описание переменных включает в себя: наименование, короткое и развернутое описание, маркер пустого значения, единицы измерения и пр. Полученные данные хранятся посуточно на сервере в четко структурированном виде: в директории с обозначением мачты находятся папки для каждой единицы оборудования на мачте, внутри – уровни обработки и далее вложенная структура директорий: год и месяц измерений. 1.2. Эксперимент, в рамках которого осуществляется всесезонный мониторинг временной изменчивости и пространственной структуры атмосферной турбулентности в условиях, приближенных к условиям «городского каньона», проводится на базе геофизической обсерватории (ГО) Института мониторинга климатических и экологических систем (ИМКЭС) СО РАН, г. Томск. В рамках выполнения проекта осенью 2021 года установлена новая 40 метровая мачта. На ней установлены акустические анемометры АМК-03 на высотах 40 м, 30 м и 2 м. что позволит получать более репрезентативные данные о набегающем ветровом потоке и проводить градиентные метеорологические, радиационные и газоаналитические измерения над и внутри лесного полога в условиях "городского каньона". Перед установкой новой мачты в сентябре 2021 проведена повторная калибровка всех акустических анемометров в климатической камере производителя (ОАО Сибаналитприбор"). Организована новая система сбора данных. Начат процесс формирования централизованной базы данных, который необходим для обеспечения доступа к разрозненным данным и возможности их постобработки. Формируются таблицы с необходимыми столбцами, их редактирование, где убираются все лишние символы, заполняются пропуски и данные выстраиваются в структуру, удобную для считывания и последующей записи. 1.3. В рамках разработки методики исследования атмосферного пограничного слоя с использованием беспилотных летательных аппаратов выполнено обобщение опыта измерений, проводившихся ранее коллективом Проекта с использованием развиваемых в ИФА РАН измерительных комплексов на базе БПЛА. В рамках первого года реализации Проекта выполнена работа по унификации, оптимизации и рефакторингу используемого набора скриптов, результатом которой стало создание библиотеки MatCopter для ПО Matlab. Разработанная библиотека содержит функции для объединения и синхронизации различных каналов данных, перехода от временных рядов к вертикальным профилям метеовеличин и их визуализации, в том числе с возможностью сравнения двух и более профилей на одном графике. Исходный код текущей версии библиотеки размещен в открытом доступе на репозитории GitHub (https://github.com/mvarentsov/MatCopter). Анализ накопленного коллективом проекта опыта применения квадрокоптеров позволил выделить ряд сценариев их применения в задачах изучения АПС над неоднородной поверхностью. Также в рамках Проекта разработана методика картирования яркостной температуры подстилающей поверхности с использованием измерительного комплекса на базе квадрокоптера Mavic 2 Zoom с тепловизором Flir TAU 2R. Апробация предложенной методики производилась в ходе полевых работ в республике Калмыкия в июле 2021 г. (съемка песчаного бархана рядом с пос. Нарын-Худук). Выполненные измерения в Калмыкии позволили установить, что в пределах относительно однородного ландшафта различия температуры поверхности, обусловленные различной экспозицией микро- и мезоформ рельефа, превышают 20 °С. 2. Для разработки методики выделения когерентных структур в рамках проекта применялись квадрантный анализ и анализ третьих моментов. Квадрантный анализ был применен для данных экспериментов над однородной поверхностью (Цимлянск) и существенно неоднородной в прибрежной зоне. Сравнение результатов показывает, что над однородной поверхностью структуры, охватывающие весь приземный слой, формируются в дневные часы, что связано с конвективными процессами. Над неоднородной поверхностью когерентные структуры присутствуют все время как над морем, так и над берегом. Для Цимлянска проведено сравнение двух периодов с близкими значениями потока импульса в дневной и ночной периоды. Квадрантный анализ показал однородность пограничного слоя в ночные часы и его значительный расбаланс в дневное время. Для более детального анализа наличия структур в турбулентных полях был привлечен анализ третьих моментов. Была исследована применимость уравнения, связывающего вертикальный перенос потока тепла и асимметрию вертикальной скорости при различных условиях стратификации в различных ландшафтах: ровная степь, прибрежная зона и городской ландшафт. Установлено, что влияние стратификации на турбулентные моменты и ассиметрию зависит от типа подстилающей поверхности и менее выражено над неоднородным ландшафтом. В случае неоднородности асимметрия также зависит от направления ветра. Над городской застройкой в случае ветра с шероховатой поверхности наблюдается большой разброс данных. Проведенный анализ показывает, что характер связи между третьим моментом (переносом потока тепла) и асимметрией вертикальной скорости может быть использован для оценки влияния ландшафта на структуру атмосферной турбулентности. На основе анализа измерений микрометеорологической мачты МО МГУ были также подтверждены следующие закономерности, полученные ранее для других неоднородных ландшафтов: 1) Турбулентные возмущения компонент скорости в условиях городской застройки в пределах первых десятков метров растут с высотой. 2) Турбулентная кинетическая энергия пропорциональна квадрату средней горизонтальной скорости. 3) Коэффициент, показывающий вклад вертикальных движений в турбулентную кинетическую энергию, больше теоретического значения для однородной поверхности. 4) Коэффициент сопротивления определяется типом поверхности, над которым формируется поток, со значением 0,08 и 0,05 для урбанизированной и озеленённой поверхности в окрестности МО МГУ, соответственно. 5) Турбулентный поток «теплового потока» удовлетворительно предсказывается диагностической связью с тепловым потоком, асимметрией и стандартным отклонением вертикальной скорости, предполагая значительный вклад когерентных структур в турбулентные потоки. Для идентификации когерентных структур в долгоживущих конвективных пограничных слоях была использована парадигма формирования организованной турбулентности. В рамках данного проекта был подготовлен массив данных для проверки теории формирования мезомасштабной турбулентности при развитой конвекции. 3. Для повышения эффективности работы с большим объемом разнородной информации была разработана база данных, включающей в себя набор данных и прикладное программное обеспечение. База данных создана на основе измерений турбулентной структуры приземного слоя атмосферы, проводимых в Институте физики атмосферы в различных физико-географических условиях. В базу данных входит информацию о географических характеристиках региона, используемом оборудовании, фоновых метеорологических условиях, характеристиках энерго- и газообмена, структуре атмосферной турбулентности (спектры, дисперсии, кинетическая энергия турбулентности). Проведенный анализ данных измерений (см. п.1-2) обусловил выбор тестовых площадок, данные которых будут использоваться в дальнейшем для анализа и разработки теории атмосферной турбулентности над неоднородной поверхностью. Кроме основных площадок в обсерваториях МГУ и ИМКЭС СО РАН, где измерения проводятся в режиме мониторинга (городской каньон) планируется проводить краткосрочные измерительные кампании в условиях однородного ландшафта (Цимлянск), прибрежной зоне (побережье Крыма) и при наличии температурных неоднородностей (полыньи Белого моря и Горьковского водохранилища). 4. Построена трёхмерная модель геометрии поверхности городской застройки, основанная на данных о положении, границах и высоте зданий из открытой базы данных OpenStreetMap. Первичная обработка и пространственная обрезка данных проводилась в ГИС с открытым исходным кодом QGIS. По итогам подобной обработки были получены карты высот зданий для окрестностей измерительной площадки Метеорологической обсерватории Московского государственного университета имени М.В. Ломоносова (МО МГУ). Для проведения тестовых расчётов обтекания городской поверхности атмосферным пограничным слоем применена трёхмерная гидродинамическая RANS модель ENVI-met. Геометрия застройки на расчётной сетке модели ENVI-met представлена для области с размерами 2500 м х 2500 м, расчётная сетка имеет размеры 500 х 500 х 120 ячеек (длина, ширина, высота), горизонтальное разрешение 5 м и вертикальное – 2.5 м. Проведён численный эксперимент с областью вокруг измерительной площадки МО МГУ размером 500 м х 500 м х 100 м. Воспроизводились условия летнего дневного прогрева поверхности и атмосферного пограничного слоя на фоне постоянного слабого западного ветра. Результаты расчётов температуры воздуха и скорости ветра показывают, что в условиях постоянного западного воздушного потока плотные группы зданий к западу и к северу от измерительной площадки служат значительными препятствиями для ветра – непосредственно за ними образуются зоны пониженных скоростей (ветровая тень), а между зданиями наблюдаются сложные вихревые структуры, в которых направление ветра может быть противоположным фоновому. В данном случае центральная и северная части измерительной площадки находятся в зоне ветровой тени, а южная – в зоне повышенных скоростей. Отмечается сильное влияние застройки на величину вертикальной скорости, особенно в слое на уровне крыш и ниже. Существенно, что в данном эксперименте часть территории МО МГУ полностью попадает в зону нисходящих движений, что может оказывать влияние на измерения скорости и направления ветра, а также турбулентных потоков. Cсылки на информационные ресурсы в сети Интернет (url-адреса), посвященные проекту: страница проекта http://ifaran.ru/ras/view/project/general.html?id=71738 данные мачты МГУ https://tower.ocean.ru данные площадки ИМКЭС РАН http://climate.scert.ru/Environment/data/IMCES_meteo/DATA_ONLINE. библиотека программ для обработки данных БПЛА https://github.com/mvarentsov/MatCopter

Аннотация результатов, полученных в 2022 году
1. Создана единая функционирующая и регулярно пополняемая база данных, включающая первичные данные и результаты первичной обработки с двух городских измерительных площадок. Пакет программ для сбора, подготовки и отправки данных написан силами участников проекта на языке Python 3.7. Исходный код пакета находится в открытом доступе по ссылке https://github.com/alexavr/tower_parse. На полигоне ИМКЭС СО РАН введена в эксплуатацию мачта высотой 40 м. с измерением турбулентности на уровнях 40 м. и 10 м. Разработано и введено в эксплуатацию новое программное обеспечение для регистрации данных анемометров АМК04, используемых на полигоне ИМКЭС РАН, позволяющее интегрировать данные анемометров в единую систему сбора и регистрации данных. Введена в эксплуатацию измерительная площадка на базе Звенигородской научной станции (неоднородный ландшафт с лесной растительностью и малоэтажной застройкой). Начат монтаж измерительной мачты высотой 56 м. в прибрежной зоне Балтийского моря.. В рамках развития методики картографирования термической неоднородности поверхности с использованием БПЛА разработан новый алгоритм коррекции и визуализации данных. После применения коррекции достигается близость температуры для схожих типов микроландшафтов (склоны одинаковой ориентации для барханов, гряды и топи для болота) в разных частях полигонов. Обновлен находящийся в открытом доступе репозиторий библиотеки на платформе разработки GitHub (https://github.com/mvarentsov/MatCopter), расширена доступная справочная информация, добавлена документация по основным функциям библиотеки и примеры ее использования. 2. Продолжена разработка различных модулей блока хранения и первичной обработки. Исходный код находится в открытом репозитории https://github.com/alexavr/tower_collect. В настоящий момент готовы и полностью функционируют модули оперативного контроля качества оперативной визуализации данных и модуль анализа турбулентных статистических моментов. Методы заполнения пропусков, реализованные авторами на основе различных алгоритмов, апробированы на примере заполнения образовавшихся пропусков в рядах температуры и вертикальной скорости ветра. Оценена эффективность различных методов заполнения пропущенных значений. В пакет обработки внесен блок расчета термического и динамического параметров шероховатости поверхности. 3. Исследована взаимосвязь между моментами второго и третьего порядков в разных условиях стратификации на двух типах контрастных ландшафтов – равнинной степи и городской местности. Установлено, что над ровной поверхностью теории подобия Монина-Обухова (ТПМО) подчиняются безразмерные турбулентные моменты вплоть до третьего порядка. Моменты третьего порядка увеличиваются с ростом неустойчивости. В устойчивых условиях турбулентный перенос теплового потока может быть связан с низкоуровневыми струями или волнами Кельвина-Гельмгольца при наличии температурной инверсии вблизи верхней границе АПС. Полученные соотношения измеренных и теоретических значений третьих моментов для условий городской застройки удовлетворительно описываются линейной регрессией с угловым коэффициентом, близким к единице. Впервые на длинном ряде наблюдений подтверждена алгебраическая зависимость моментов разных порядков, предполагающая наличие организованных турбулентных структур, в условиях безразличной и устойчивой стратификации над поверхностью сложной геометрической формы. Анализ моментов третьего порядка показывает, что нелокальный перенос из-за больших вихрей может играть важную роль, особенно в городском ландшафте, и должен быть включен в турбулентное замыкание. По данным, полученным на городских площадках, исследовалось отклонение наблюдений турбулентной кинетической энергии и потока импульса от модели замыкания турбулентности EFB (energy flux balance) в однородных условиях. Проведена коррекция модели для условий неоднородной поверхности. 4. Разработана методика выделения когерентных структур из турбулентных данных по анализу завихренности. Для данных над однородной поверхностью (Цимлянск) в ночной период хорошо прослеживаются волны в ночные часы и конвективные вихри при дневном прогреве. В случае неоднородного наблюдается существенная завихренность на уровне 10-12 м, совпадающем с уровнем высоты деревьев, что может свидетельствовать о течениях, сходных с течением за уступом. Была проведена идентификация когерентных структур эйлеровыми (локальными) методами. Для отработки методики использовались когерентные структуры, выделенные в течениях Экмана и Куэтта. На когерентные структуры приходится треть турбулентной кинетической энергии и примерно половина энстрофии и спиральности. Показано, что вклад крупномасштабных структур в пристеночную динамику турбулентных течений, дисперсии скоростей и концентрации примеси значителен. 5. Специализированный эксперимент проводился на болотном стационаре Мухрино над термически неоднородной поверхностью. На основе полученных уточненных значений характеристик подстилающей поверхности производился модельный расчет турбулентных потоков тепла по теории подобия Монина-Обухова, результат сравнивался с модельным расчетом потока тепла по ТПМО без учета термической неоднородности поверхности и с данными пульсационных измерений турбулентных потоков над неоднородной поверхностью. Из полученного результата следует, что учет термической неоднородности подстилающей поверхности болот положительно сказывается на точности модельного расчета потока тепла. 6. В рамках проведения численных экспериментов с технологией DNS/LES/RANS для воспроизведения структуры и динамики пограничного слоя в окрестности городских измерительных площадок проекта показано, что потенциальное влияние планируемой установки 40-метровой мачты на измеряемые характеристики ветрового потока значительно слабее, чем влияние уже существующих крупных препятствий. В рамках прямого численного моделирования течения в открытом канале была реализована возможность задания неоднородной поверхности в виде четырех различных конфигураций городской застройки: кубики, расположенные в шахматном порядке плотно и разреженно, и каньоны, направленные вдоль и поперек потока. Было получено, что сильнее всего ослабление потока в нижней части пограничного слоя происходит при конфигурации каньонов, стоящих перпендикулярно потоку. Для конфигурации каньонов, перпендикулярных потоку, основная доля ТКЭ приходится на вертикальную компоненту, а не продольную, как в остальных случаях, что предполагает значительное вертикальное обтекание препятствия воздушным потоком. Установлено, что свойства турбулентного обмена значимым образом зависят от конфигурации городской застройки. Ослабление средней скорости в нижней части пограничного слоя сильнее всего проявляется при конфигурации каньонов, перпендикулярных фоновому потоку. 7. На основе результатов прямого численного моделирования течений Куэтта и Пуазейля (течение в открытом канале), проведена проверка гипотезы о диагностической связи между вторыми и третьими моментами, содержащей константу С~1. Для расчётов использован универсальный программный код DNS/RANS/LES, разработанный в НИВЦ МГУ и ИВМ РАН. Большая часть значений C для обоих течений оказалась в диапазоне 0.5-1 на разных высотах в канале. Проведено сравнение результатов этих численных экспериментов с данными измерений на микрометеорологической мачте МО МГУ, по параметрам регрессии между фактическим третьим моментом и моментом, заданным формулой Зилитинкевича. Формы распределения точек на диаграммах рассеяния оказались схожи, коэффициенты линейной регрессии составили 0.65 и 0.61 для данных натурных измерений и прямого численного моделирования, соответственно. Cсылки на информационные ресурсы в сети Интернет (url-адреса), посвященные проекту: страница проекта http://ifaran.ru/ras/view/project/general.html?id=71738 данные мачты МГУ https://tower.ocean.ru данные площадки ИМКЭС РАН http://climate.scert.ru/Environment/data/IMCES_meteo/DATA_ONLINE.

Аннотация результатов, полученных в 2023 году
1. Данные микрометеорологических полигонов на базе метеообсерватории МГУ и геофизического полигона в ИМКЭС СО РАН подготовлены для интеграции в международную систему UrbanFlux и сеть FluxNet. Мачты интегрированы в разрабатываемую систему RuFlux. К единой системе сбора данных подключены мачты на научной базе университете БФУ им. Канта и на базе Звенигородской научной станции ИФА им. А.М. Обухова РАН, завершена интеграция системы сбора данных на базе геофизической обсерватории ИМКЭС (Томск). Быстрое и успешное подключение имеющихся метеорологических мачт доказывает, что разработанная система обладает легкой масштабируемостью. Завершена разработка пакета программ для первичной обработки и статистического анализа данных. В рамках первичной обработки данных происходит удаление недостоверных данных, заполнение пропусков, коррекция акустической температуры и коррекция на влажность. Статистический анализ включает в себя спектральный и корреляционный анализ, расчет турбулентных статистик вторых и высших порядков. 2. Выполнена работа по заключительной отладке и чистовому оформлению методики мониторинга и картографирования термической неординарности подстилающей поверхности с использованием БПЛА и оригинального алгоритма коррекции погрешностей и артефактов, возникающих в регистрируемых данных. На основе обработанных данных выполнен анализ статистических характеристик термической неоднородности поверхности для этих ландшафтов в дневное и ночное время суток. Выполнена адаптация библиотеки MatCopter для работы с данными новых измерительных платформ. Добавлена возможность синхронизации данных с нескольких установленных на квадрокоптере измерительных систем. По результатам анализа измерений, выполненных с использованием квадрокоптеров в ходе эксперимента в Мошковском заливе Иваньковского водохранилища установлено, что влияние полыньи на вертикальную структуру профилей температуры и скорости ветра мало, и выражено лишь в наличие сверхадиабатических градиентов температуры над полыньей в тонком слое вблизи поверхности. Результаты измерений демонстрируют наличие интенсивного энергообмена между полыньей и атмосферой. О конвективном характере турбулентных движений, а именно о наличии интенсивных теплых восходящих потоков свидетельствуют также положительные значения коэффициентов асимметрии вертикальной скорости и температуры. Анализ измерений на трех измерительных площадках на побережье Черного моря (пос. Кацивели, Крым) показал значительное различие в характеристиках атмосферного пограничного слоя над морем, и на различном удалении от берега. Использование данных мачты на Балтийском море с профильными измерениями метеопараметров позволили проверить применимость универсальных функций теории подобия Монина-Обухова в прибрежной зоне. При ветре с моря ТПМО применима, но при береговых ветрах условия подобия не работают. При масштабировании функций на стандартное отклонение и скорость диссипации отклонения менее существенны. По результатам измерений характеристик атмосферной турбулентности над озерами и в прибрежной зоне получены параметризации для расчета динамического параметра шероховатости и параметров шероховатости для температуры и влажности для мелкого водоема. Полученные параметризации могут быть использованы в моделях Земной системы и озерных моделях для расчета турбулентных потоков над водными объектами суши. 3. Определены границы применимости теорий горизонтально статистически однородных пограничных слоёв над неоднородными ландшафтами различных типов. В случае прибрежной зоны значимую роль играют направление ветра, стратификация атмосферы и топография прибрежной зоны. Результаты вихреразрешающего моделирования показывают, что поток импульса в прибрежной зоне чувствителен к скорости и направлению ветра, а также топографии суши. В случае городской застройки важную роль в нарушении критериев подобия играют вторичные циркуляции, создаваемые сложным рельефом. По данным DNS моделирования получено, что применимость ТПМО в случае городской застройки определяется высотой измерений. Измерения над ландшафтом позволяют получить интегральные значения потоков без учета вторичных циркуляций. В случае термической неоднородности турбулентные потоки определяются областью формирования потока (футпринтом). Сравнение экспериментальных и расчетных данных показало, что ТПМО применима при условии агрегирования потоков внутри футпринта. Сравнение экспериментальных и расчетных данных показали применимость этого подхода. 4. Для выбора критериев подобия для неоднородной поверхности были использованы данные, полученные на мачтах над однородной поверхностью (Цимлянск) и над городской застройкой (Москва и Томск). Проведено сравнение различных критериев для однородной и неоднородной поверхности. Установлено, что для неоднородной поверхности возможно ввести масштаб длины lm, который по своему смыслу является путем смешения Прандтля, и фактически определяет топологические характеристики турбулентности в нейтральных условиях внутри одного направления ветра. 5. На основании анализа полученных в проекте экспериментальных данных и мелкомасштабного моделирования определены подходы к модификации ТПМО для неоднородных ландшафтов. Первый подход основан на анализе соотношения вторых и третьих моментов, что позволяет выявлять наличие когерентных структур в потоке и, соответственно, условия неприменимости ТПМО. Установлено, что классические соотношения ТПМО непригодны для описания безразмерных градиентов скорости ветра и температуры над неоднородным ландшафтом, так как не учитывается нелокальный перенос. В качестве новых критериев предлагается использовать соотношение дисперсии вертикальной скорости ветра и динамической скорости и масштаб длины, характеризующий турбулентное перемешивание и расстояние до ближайших неоднородностей. Разработана обобщенная модель масштабирования сложной атмосферной турбулентности с учетом анизотропии. Показано, что различные состояния анизотропии соответствуют разным отношениям подобия, особенно при неустойчивой стратификации. Также одним из подходов к созданию обобщенной ТПМО является понимание того факта, что масштаб горизонтальной турбулентности также является важным параметром для случая конвективных пограничных слоев или формирования когерентных структур над неоднородной поверхностью. 6. В ходе серии численных экспериментов с DNS-версией RANS/LES/DNS создана конфигурации городской застройки в окрестностях измерительных площадок проекта и адаптирована RANS/LES/DNS модель к воспроизведению структуры и динамики атмосферного пограничного слоя в районах площадок. Получены результаты несколькольких серий экспериментов, которые озволяют делать качественные выводы об изменчивости и масштабах рассматриваемых атмосферных процессов. Полученные данные позволяют провести более глубокий анализ и сравнение с другими источниками данных (измерения и другие модели) как для уточнения и улучшения RANS/LES/DNS моделей, так и для изучения особенностей и ограничений измерительных методов. 7. Определена наиболее оптимальная конфигурации модели WRF-ARW для прогнозирования характеристик энергообмена с учетом особенностей неоднородной подстилающей поверхности. Установлено, что модель чувствительна к выбору параметризаций. На примере выбранной конфигурации модели проведена оценка адекватности мозаичного подхода учета вклада неоднородностей подстилающей поверхности в моделях с грубой сеткой. В большинстве случаев мозаичный метод занижает абсолютные значения средних потоков, что связано с увеличением реальных локальных турбулентных потоков за счет мезомасштабных циркуляций. Относительная ошибка расчета турбулентных потоков мозаичным методом зависит от абсолютной величины этих потоков; при достаточно высоких значениях потоков мозаичный метод обеспечивает ошибку воспроизведения средних потоков не более 10%; ошибка мозаичного метода существенно нелинейно зависит от характеристик подстилающей поверхности и метеорологических условий.