Рекомендации методические мр 10. 0057-12 Государственное санитарно-эпидемиологическое нормирование Российской Федерации 2.1.10. СОСТОЯНИЕ ЗДОРОВЬЯ НАСЕЛЕНИЯ В СВЯЗИ С СОСТОЯНИЕМ ОКРУЖАЮЩЕЙ СРЕДЫ И УСЛОВИЯМИ ПРОЖИВАНИЯ НАСЕЛЕНИЯ Оценка риска и ущерба от климатических изменений, влияющих на повышение уровня заболеваемости и смертности в группах населения повышенного риска Методические рекомендации МР 2.1.10. 0057 -12 Москва - 2012 Оценка риска и ущерба от климатических изменений, влияющих на повышение уровня заболеваемости и смертности в группах населения повышенного риска: Методические рекомендации. —М.: Федеральный центр гигиены и эпидемиологии Роспотребнадзора, 2012.— 48 с. 1. Разработаны: Федеральная служба по надзору в сфере защиты прав потребителей и благополучия человека (Г.Г. Онищенко, И.В. Брагина, А.Л. Мишина), Управление Роспотребнадзора по Архангельской области (Р.В. Бузинов, Т.Н. Унгуряну), Управление Роспотребнадзора по Свердловской области (С.В. Кузьмин, О.Л. Малых), ФБУН «Екатеринбургский медицинский научный центр профилактики и охраны здоровья рабочих промпредприятий» Роспотребнадзора (В.Б. Гурвич, С.В. Ярушин, Б.А. Кацнельсон), ФБУН «Федеральный научный центр медико-профилактических технологий управления рисками здоровью населения» Роспотребнадзора (Н.В. Зайцева, И.В. Май, П.З. Шур, Д.А. Кирьянов, В.Б. Алексеев), Институт народнохозяйственного прогнозирования РАН (Б.А. Ревич). При подготовке приложения 2 к настоящим методическим рекомендациям использованы данные Управления Роспотребнадзора по Воронежской области (В.И. Чубирко, Н.М. Пичужкина). 2. Утверждены руководителем Федеральной службы по надзору в сфере защиты прав потребителей и благополучия человека, Главным государственным санитарным врачом Российской Федерации Г.Г. Онищенко «_17__» __января__2012 года. 3. Введены в действие с «_17_» __января__2012 года. Содержание I. Область применения........................................................................................4 II. Общие положения. .........................................................................................5 III. Термины и определения................................................................................6 IV. Нормативные ссылки.....................................................................................7 V. Методика сбора и подготовки данных .........................................................7 5.1. Сбор данных о метеорологических показателях...................................7 5.2. Сбор данных о климаточувствительных заболеваниях........................9 5.3. Сбор данных для экономической оценки ущерба...............................11 5.4. Подготовка данных.................................................................................12 VI. Методы оценки влияния метеорологических факторов на здоровье.......14 6.1. Эпидемиологический анализ помесячной динамики показателей здоровья....................................................................................14 6.2. Статистический анализ..........................................................................15 6.3. Представление результатов...................................................................16 6.4. Оценка популяционного риска заболеваний, смертей, обращений за скорой медицинской помощью, госпитализации..............17 VII. Оценка экономического ущерба, связанного с риском для здоровья населения при воздействии метеорологических факторов………………......................................................................................17 Приложение 1. Пример базы данных по метеорологическим факторам и количеству вызовов скорой медицинской помощи в городе N..................19 Приложение 2. Оценка воздействия климатических изменений на здоровье населения в различных регионах России ...................................20 1. Оценка влияния температуры воздуха, волн жары и холода на смертность населения…………………………………………………..20 2. Оценка влияния температуры воздуха, волн жары и холода на обращаемость населения за скорой медицинской помощью …….…39 3. Оценка влияния температуры и загрязнения воздуха на заболеваемость населения …………………………………………..…43 4. Оценка влияния температуры воздуха на заболеваемость инфекционными болезнями ……………………………………………44 Приложение 3. Перечни публикаций…………………………………………46 УТВЕРЖДАЮ Руководитель Федеральной службы по надзору в сфере защиты прав потребителей и благополучия человека, Главный государственный санитарный врач Российской Федерации __________________ Г.Г. Онищенко «_17_»__января__2012 г. Дата введения: с момента утверждения 2.1.10. СОСТОЯНИЕ ЗДОРОВЬЯ НАСЕЛЕНИЯ В СВЯЗИ С СОСТОЯНИЕМ ОКРУЖАЮЩЕЙ СРЕДЫ И УСЛОВИЯМИ ПРОЖИВАНИЯ НАСЕЛЕНИЯ Оценка риска и ущерба от климатических изменений, влияющих на повышение уровня заболеваемости и смертности в группах населения повышенного риска Методические рекомендации МР 2.1.10.0057-12 __________________________________________________________________ I . Область применения 1.1. В настоящих методических рекомендациях представлены методические подходы к сбору и подготовке данных, оценке влияния метеорологических показателей на здоровье населения, оценке экономического ущерба, связанного с повышением уровня заболеваемости и смертности от климаточувствительных заболеваний в группах населения повышенного риска. 1.2. Методические рекомендации предназначены для органов и организаций Федеральной службы по защите прав потребителей и благополучия человека, а также научно-исследовательских и иных организаций, деятельность которых связана с оценкой риска для здоровья населения при воздействии факторов среды обитания. 1.3. Положения методических рекомендаций могут быть использованы при: · обосновании необходимости разработки и реализации мероприятий по адаптации к изменению климатических условий, рассмотрении вариантов социальных, экономических и технических мер по предупреждению и смягчению их последствий для уменьшения потенциального неблагоприятного воздействия изменения климата на состояние здоровья населения; · совершенствовании системы мониторинга состояния здоровья населения и среды обитания в связи с изменением климата для прогнозирования, предупреждения и профилактики заболеваний, чувствительных к климатическим факторам; · ориентировочной оценке потенциального экономического ущерба, причиненного здоровью человека и целесообразных затрат на реализацию мер по обеспечению санитарно-эпидемиологического благополучия; · обеспечении объективной информацией лиц, участвующих в принятии управленческих решений, населения, общественных организаций об уровнях риска для здоровья населения в связи с неблагоприятным воздействием климатических факторов. II . Общие положения 2.1. В течение последних десятилетий вызывает обеспокоенность беспрецедентно высокая скорость глобального потепления и изменения климата. Установлено, что хозяйственная деятельность человека оказывает заметное влияние на климат. Эти изменения климата многообразны и проявляются в изменении частоты и интенсивности климатических аномалий и экстремальных погодных явлений. Ожидаемые изменения климата неизбежно отразятся на жизни людей во всех регионах планеты, а в некоторых из них станут ощутимой угрозой для благополучия населения. В последние годы изменение климата рассматривается как один из ведущих факторов, оказывающих влияние на здоровье населения. Происходит как прямое влияние за счет увеличения числа дней с аномально высокими и/или низкими температурами, наводнений, штормов, тайфунов, так и косвенное, опосредованное влиянием экологических или социально-экономических факторов (увеличение площади засушливых земель, уменьшение объемов качественной питьевой воды и др.). Воздействие изменения климата в городах сочетается с неблагоприятным воздействием на здоровье населения загрязненного атмосферного воздуха. 2.2. По оценкам ВОЗ в Европе ежегодно климатические изменения являются причиной от 1 до 10% смертей среди старших возрастных групп, а в мире - более 150 тысяч дополнительных смертей и 5,5 млн. лет нетрудоспособности в год. Это составляет 0,3% общего числа смертельных исходов и 0,4% общего количества лет нетрудоспособности соответственно. К 2050 году ожидается дальнейшее увеличение числа смертельных исходов, связанных с потеплением климата еще примерно на 1-1,5%. Экономический ущерб от дополнительной смертности в результате климатических изменений в мире колеблется в большом интервале - от 6 до 88 млрд. долларов в год. 2.3. В соответствии с Климатической доктриной Российской Федерации к отрицательным последствиям ожидаемых изменений климата для Российской Федерации относится повышение риска для здоровья (увеличение уровня заболеваемости и смертности) населения. Оценка рисков и связанных с ними потерь рассматриваются как важнейшие составляющие при разработке и планировании мер по адаптации к изменениям климата. Количественная оценка риска позволяет определить примерную величину конкретных последствий (заболеваний и преждевременной смерти) при различных сценариях. 2.4. Основными факторами риска, связанными с климатическими изменениями, считаются высокие температуры, экстремальные погодные явления, распространение инфекционных заболеваний, нарушения питания и пр. В качестве одного из приоритетных факторов климатических изменений, влияющих на повышение уровня заболеваемости и смертности населения Российской Федерации, рассматриваются «волны» жары и холода. В настоящих методических рекомендациях к числу наиболее чувствительных к потенциальным воздействиям этих факторов категорий населения, для которых может проводиться оценка риска, относятся лица пожилого возраста и детское население. 2.5. Основные положения настоящих методических рекомендаций гармонизированы с методическими подходами, изложенными в документах Всемирной организации здравоохранения [1-4]. III . Термины и определения Адаптация – стратегии, направления политики и меры, осуществляемые в настоящее время и в будущем с целью уменьшения потенциальных неблагоприятных воздействий на здоровье, обусловленных изменчивостью и изменением климата. Заболевание, характеризующееся чувствительностью к климатическим факторам (климаточувствительное заболевание) – заболевание, тесно связанное с метеорологическими и/или климатическими факторами, что находит отражение в показателях территориальной заболеваемости, смертности или других показателях состояния здоровья населения и сезонном характере передачи заболеваний. Изменение климата – статистически значимое изменение либо среднего состояния климата, либо его изменчивости, которое сохраняется в течение длительного времени (обычно несколько десятилетий или дольше). По определению Всемирной Метеорологической организации (ВМО), периодом для отсчета изменений климата считается 30 лет. Изменчивость климата относится к колебаниям вокруг среднего состояния, включая случаи экстремальных погодных явлений. Климат – это средний режим погоды или статистическое описание средней величины и изменчивости соответствующих количественных параметров в течение определенного периода времени. Погода – повседневное проявление климата в определенном месте в определенное время. Чувствительность – более высокая степень восприимчивости отдельных людей и систем к неблагоприятным последствиям изменения климата, включая его изменчивость и экстремальные климатические явления, или неспособность преодолеть эти последствия. Эффективная температура – мера относительного дискомфорта, который испытывает человек в связи с воздействием на организм жары в сочетании с высокой влажностью воздуха. IV . Нормативные ссылки 4.1. Федеральный закон от 30 марта 1999 г. № 52-ФЗ «О санитарно-эпидемиологическом благополучии населения». 4.2. Постановление Правительства Российской Федерации от 30 июня 2004 № 322 «Об утверждении Положения о Федеральной службе по надзору в сфере защиты прав потребителей и благополучия человека». 4.3. Распоряжение Правительства Российской Федерации от 25 апреля 2011 г. № 730р «Об утверждении Комплексного плана реализации Климатической доктрины Российской Федерации на период до 2020 года» 4.4. МР 5.1.0029-11 «Методические рекомендации к экономической оценке рисков для здоровья населения при воздействии факторов среды обитания». 4.5. МР 5.1.0030-11 «Методические рекомендации к экономической оценке и обоснованию решений в области управления риском для здоровья населения при воздействии факторов среды обитания». V . Методика сбора и подготовки данных 5.1.Сбор данных о метеорологических показателях 5.1.1. К метеорологическим показателям, используемым для оценки влияния изменений климата на здоровье населения, относятся среднесуточные и максимальные показатели температуры атмосферного воздуха, относительной влажности, атмосферного давления, скорости ветра и количество осадков. 5.1.2. Источником информации по метеорологическим данным является межрегиональный территориальный орган (управление) по гидрометеорологии и мониторингу окружающей среды и подведомственные государственные учреждения – центры по гидрометеорологии и мониторингу окружающей среды. 5.1.3. Для одновременной оценки температуры и влажности воздуха целесообразно использовать эффективную температуру воздуха, которая рассчитывается по формуле [21]: , где где Tapp – эффективная температура, Tair – температура воздуха, а Tdewpt – точка росы. Точка росы рассчитывается на основании информации об относительной влажности по формуле: , где γ рассчитывается по формуле: , где Т – температура воздуха, RH – относительная влажность, a и b – константы, равные 17.271 и 237.7 соответственно. Для расчета среднесуточной эффективной температуры используются среднесуточные значения температуры и влажности, а для расчета максимальной эффективной температуры – их максимальные значения. 5.1.4. При отсутствии метеостанции на изучаемой территории можно использовать данные близлежащих метеостанций. Для определения среднесуточной температуры используют данные метеостанции, расположенной в радиусе до 300 км, для среднемесячной температуры – в радиусе до 1200 км. Для определения среднесуточной влажности и суточного количества осадков используют данные метеостанции, расположенной в радиусе до 50 км, для среднемесячной влажности и количества осадков – в радиусе до 400 км. 5.1.5. В условиях влияния метеорологических факторов на здоровье населения, особенно на территориях с техногенным загрязнением, необходимо оценить модифицирующий эффект от влияния загрязняющих атмосферный воздух веществ (взвешенные вещества, диоксид серы, диоксид азота, оксид азота, оксид углерода и других веществ, включенных в программу мониторинга). 5.1.6. Источником информации по мониторингу загрязнения атмосферного воздуха являются данные межрегионального территориального органа (управления) по гидрометеорологии и мониторингу окружающей среды и подведомственных государственных учреждений – центров по гидрометеорологии и мониторингу окружающей среды, лабораторий промышленных предприятий и организаций. Исходной информацией о загрязнении атмосферного воздуха являются результаты замеров на постах наблюдения. 5.2. Сбор данных о климаточувствительных заболеваниях 5.2.1. К наиболее распространенным заболеваниям, характеризующимся чувствительностью к климатическим факторам, относятся: - Болезни органов дыхания: острые респираторные инфекции (J00-J22), аллергический ринит (J30), хронические обструктивные болезни легких (J40-J44), бронхиальная астма (J45); - Болезни органов кровообращения: болезни, характеризующиеся повышенным кровяным давлением (I10-I15), ишемическая болезнь сердца (I20-I25), нарушения проводимости и сердечного ритма (I44-49), цереброваскулярные болезни (I60-I69); - Болезни эндокринной системы: сахарный диабет (Е10-14); - Травмы, отравления и другие последствия воздействия внешних причин: травмы (S00-Т14), отморожения (Т33-35), утопления (W69-70), суициды (Х60-84); - Психические расстройства: алкогольные психозы (F10); - Кишечные инфекции: сальмонеллез (А02), шигеллез (А03) и другие бактериальные кишечные инфекции (А04) и отравления (А05), вирусные кишечные инфекции (А08); - Трансмиссивные болезни: болезнь Лайма (А69.2), клещевой вирусный энцефалит (А84), лихорадка Денге (А90-91), лихорадка Западного Нила (А92.3), желтая лихорадка (А95), малярия (В50-54), лейшманиоз (В55), африканский трипаносомоз (В56), болезнь Шагаса (В57), охоцеркоз (В73), филяриатоз (В74). 5.2.2. К группам населения повышенного риска, для которых проводится оценка влияния климатических изменений на показатели здоровья, относятся: - дети (от 0 до 17 лет); - лица пожилого (60 лет и старше) и престарелого возраста (75 лет и старше); - лица с хроническими заболеваниями, страдающие болезнями органов дыхания, кровообращения, болезнями нервной, мочевыделительной и эндокринной системами; - лица в трудоспособном возрасте, в том числе работающие на открытом воздухе; - коренные народы. 5.2.3. Оценку влияния метеорологических факторов на здоровье следует выполнять по следующим показателям: - смертность населения; - заболеваемость населения по обращаемости; - заболеваемость населения по данным госпитализации; - обращаемость населения за скорой медицинской помощью. 5.2.4. Первичным документом, регистрирующим случаи смерти, является форма № 106/у-8 «Медицинское свидетельство о смерти» и «Запись акта о смерти», выдаваемая ЗАГСами. Сбор первичной информации по смертности населения для исследования проводится двумя путями: выкопировкой данных из формы № 106/у-8 и записи акта о смерти в ЗАГСах по всем случаям смерти в городе (районе); на основе существующих электронных баз данных по смертности, формируемых в медицинских информационно-аналитических центрах в субъектах Российской Федерации или в территориальном органе Федеральной службы государственной статистики. 5.2.5. Источниками информации о случаях заболеваний являются «Талон амбулаторного пациента» (ф. 025-6(7)/у-89), «Медицинская карта стационарного больного» (ф. 003/у). Сбор первичной информации осуществляется путем выкопировки или на основе электронных баз данных лечебно-профилактических организаций, территориальных фондов обязательного медицинского страхования. 5.2.6. Источником информации о регистрации ежедневных случаев инфекционных заболеваний является «Экстренное извещение об инфекционном заболевании, пищевом, остром профессиональном отравлении, необычной реакции на прививку» (ф. 058/у) или журнал учета инфекционных заболеваний (форма 060/у). Информацию о количестве заболевших за месяц можно получить из отчетных форм «Сведения об инфекционных и паразитарных заболеваниях» (форма 2), «Сведения о вспышках инфекционных заболеваний» (форма 23-09). Данные о количестве случаев инфекционных заболеваний регистрируются центром гигиены и эпидемиологии Роспотребнадзора в субъекте Российской Федерации. 5.2.7. Сбор информации о вирусном клещевом энцефалите включает количество пострадавших от нападения клещей («покусанных клещами»), количество заболевших клещевым энцефалитом, показатели обилия клещей («заклещевленность» территории), данные об инфицировании клещей вирусом клещевого энцефалита. 5.2.8. Источником информации по каждому случаю обращения за скорой медицинской помощью является «Карта вызова скорой медицинской помощи» (ф. 110/у). Ежедневные данные по обращаемости за скорой медицинской помощью формируются на станциях скорой медицинской помощи или в отделениях скорой медицинской помощи. 5.2.9. Для количественного определения связи между метеорологическими параметрами и показателями здоровья населения необходимы ежедневные данные, или сгруппированные за неделю, или месяц, о количестве смертей, заболеваний, госпитализаций, обращений за скорой медицинской помощью на изучаемой территории. Если в анализе используются ежедневные данные, то рекомендуемый период наблюдения должен составлять не менее 6 лет, а если данные за месяц, то не менее 20 лет. 5.2.10. Источником информации по среднегодовой численности населения в возрастно-половой структуре является территориальный орган Федеральной службы государственной статистики. 5.3. Сбор данных для экономической оценки ущерба 5.3.1. К экономическим показателям для расчета ущерба здоровью населения в результате климаточувствительных заболеваний относятся показатели, характеризующие экономические потери от повышения уровня заболеваемости и смертности в группах риска, связанного с метеорологическими факторами в период имеющих место (исследуемый период) последствий. 5.3.2. Данные для экономической оценки ущерба должны соответствовать требованиям МР 5.1.0029-11 «Методические рекомендации к экономической оценке рисков для здоровья населения при воздействии факторов среды обитания». 5.3.3. Источниками данных для расчета экономического ущерба от повышения уровня смертности населения, связанного с метеорологическими условиями, являются: - территориальный орган Федеральной службы государственной статистики, предоставляющий информацию о годовом валовом региональном продукте на душу населения и среднегодовой заработной плате в регионе; - медицинские информационно-аналитические центры в субъектах Российской Федерации и территориальный орган Федеральной службы государственной статистики, предоставляющие данные о среднем возрасте умерших в исследуемый период. 5.3.4. Источниками данных для расчета экономического ущерба от повышения уровня заболеваемости населения по обращаемости, заболеваемости по данным госпитализации и обращаемости населения за скорой медицинской помощью, связанного с метеорологическими условиями (по исследуемым климаточувствительным заболеваниям и группам населения повышенного риска), являются: - медицинские информационно-аналитические центры в субъектах Российской Федерации или лечебно-профилактические организации, территориальные фонды обязательного медицинского страхования, предоставляющие данные о стоимости одного обращения по заболеванию, одного дня лечения заболевания в амбулаторных условиях и в стационаре, стоимости одного вызова скорой медицинской помощи, среднем числе дней лечения в амбулаторных и стационарных условиях, а также средние расходы на приобретение фармацевтических и медицинских препаратов, на оказание медицинских услуг, получаемых по заболеванию сверх государственных гарантий; - территориальные фонды социального страхования, предоставляющие данные о среднем размере оплаты одного дня нетрудоспособности по больничным листам по заболеваниям, среднее число целодневных потерь рабочего времени в расчете на один случай заболевания; - территориальный орган Федеральной службы государственной статистики, предоставляющий данные о численности работающего населения (доли работающих в общей численности исследуемой группы населения); - территориальные налоговые органы Федеральной налоговой службы, предоставляющие данные о поступлении отчислений от налогов на прибыль организаций, налогов на доходы физических лиц, единого социального налога за рассматриваемый период; - территориальные управления социальной защиты населения, предоставляющие данные о среднегодовом размере пенсионного обеспечения по болезни (инвалидности) и доле пенсионеров- инвалидов в общей численности исследуемой группы населения. 5.4. Подготовка данных (пример базы данных приведен в приложении 1 к настоящим методическим рекомендациям) 5.4.1. Метеорологические данные должны пройти трехступенчатый контроль: синтаксический, семантический и пространственный. Как правило, наблюдения метеорологических параметров на метеостанциях включают восемь измерений за сутки. Для оценки влияния температуры воздуха, относительной влажности, атмосферного давления и скорости ветра на здоровье необходимо рассчитать среднесуточное и максимальное значения за сутки. 5.4.2. При оценке модифицирующего действия загрязняющих атмосферный воздух веществ для анализа используют среднесуточную и максимальную концентрации за сутки. 5.4.3. Данные о количестве смертей, заболеваний, госпитализаций, вызовов скорой медицинской помощи и данные о метеорологических переменных должны быть представлены в одинаковом временном и географическом разрешении. 5.4.4. Данные о показателях здоровья должны анализироваться как в целом среди всего населения, так и по отдельным возрастным группам, например, дети (0-17 лет), трудоспособное население (18-59 лет) и лица пожилого возраста (60 лет и старше). 5.4.5. Для того чтобы количественно определить климаточувствительные заболевания, необходимо иметь суточные данные о количестве заболеваний (смертей), госпитализаций, вызовов скорой медицинской помощи по отдельным классам, рубрикам или нозологическим формам болезней. Однако если случаев недостаточно или они регистрируются не ежедневно, можно суммировать случаи за неделю или за месяц. Соответственно метеорологические факторы также должны быть представлены как средние и максимальные значения за неделю или месяц за аналогичный временной интервал. 5.4.6. Небольшое количество случаев по отдельным рубрикам или нозологическим формам болезней является нерепрезентативным и не позволит выявить статистически значимые связи. Поэтому целесообразно объединять отдельные рубрики и нозологические формы в соответствующие классы болезней, чтобы увеличить количество случаев. 5.4.7. Если анализируются данные скорой медицинской помощи, то вызовы не следует представлять по отдельным рубрикам или нозологическим формам болезней в связи с трудностью верификации точного диагноза, а целесообразно всегда изучать в целом по классам заболеваний. 5.4.8. Для устранения влияния потенциальных факторов, способных повлиять на итоговые результаты оценки, в базу данных необходимо включить информацию о годах, месяцах, днях недели, праздничных днях, эпидемиях гриппа. При изучении острых кишечных инфекций из базы данных нужно удалить пострадавших при групповой и вспышечной заболеваемости. 5.4.9. Проводится логический контроль всех вводимых данных по метеорологическим показателям. Например, поле «температура» должно принимать значения в пределах, соответствующих сезонам года. Контроль качества данных удобно проводить с помощью графических методов. Построение графика переменной позволяет легко обнаружить «подозрительные» значения, выпадающие из общего ряда, например, ошибку знака, допущенную при вводе значения температуры. 5.4.10. Проводится контроль данных по заболеваемости, госпитализации, обращаемости за скорой медицинской помощью, смертности населения. Например, значение поля «дата рождения» должно быть не больше значения поля «дата смерти»; значение каждого из этих полей должно быть не больше даты заполнения базы. Затем проверяется соответствие возраста, даты рождения и даты смерти. Проводится логический контроль всех закодированных данных (таких как причины смерти, заболеваний и т.п.): они должны принимать только допустимые значения, указанные в справочниках. 5.4.11. Базы данных, содержащие информацию о том или ином показателе здоровья и метеорологических факторах, необходимо проверить на пропущенные значения. Если в базе данных отсутствует информация за три дня подряд, то вместо пропущенных значений можно использовать среднемесячные данные. Если в базе данных отсутствует информация за один день, то вместо пропущенного значения можно использовать среднее из смежных дней. После заполнения пропусков в базе данных нужно выполнить анализ чувствительности. Если в базе данных отсутствует 30% значений для того или иного метеорологического параметра, то такая переменная должна быть исключена из анализа. 5.4.12. При оценке влияния метеорологических факторов на количество заболеваний и госпитализаций воскресные дни можно исключить, если в модели не будут использоваться лаги. Под лагом в данном случае понимается тот интервал в днях (или месяцах) между изменением метеорологических показателей (или показателей загрязнения атмосферного воздуха) и изменением смертности, заболеваемости и др., при котором между ними найдена статистически значимая связь. Если модель предполагается с лагом, то для воскресенья используется среднее значение из количества заболеваний/госпитализаций за субботу и понедельник. 5.4.13. При обнаружении выскакивающих значений количества случаев того или иного показателя, которые выходят за пределы трех стандартных отклонений, следует заменить данный «выброс» значением, соответствующим двум стандартным отклонениям. VI . Методы оценки влияния метеорологических показателей на состояние здоровья населения 6.1. Эпидемиологический анализ помесячной динамики показателей здоровья 6.1.1. Под сезонностью понимаются закономерные колебания изучаемого показателя состояния здоровья на протяжении календарного или эпидемического года, выражаемые многолетней приуроченностью подъемов и снижений его уровня к определенным месяцам, обусловленные воздействием периодических (сезонных) факторов. 6.1.2. Для изучения особенностей помесячной динамики показателей здоровья населения за многолетний период используются индекс сезонности и атрибутивная фракция сезонных причин. Оба показателя рассчитываются для каждого месяца за многолетний период на основе интенсивных показателей смертности, заболеваемости, госпитализации, обращаемости за скорой медицинской помощью. 6.1.3. Для расчета индекса сезонности интенсивные показатели суммируются отдельно для каждого месяца за все годы анализируемого многолетнего периода. Полученные итоговые значения делятся на число лет наблюдения, в результате получаются 12 месячных средних (Yi ). Для расчета общей средней (Y0 ) сумма месячных средних делится на 12. Индекс сезонности рассчитывается по формуле: , где I – индекс сезонности в % Yi – месячные средние за многолетний период Y0 – общая средняя Если значение индекса сезонности для какого-либо месяца превышает 100%, то считается, что в этом месяце активизировались сезонные факторы. 6.1.4. Атрибутивная фракция сезонных причин (сезонное воздействие) показывает, какой процент случаев (заболеваний, смертей) обусловлен сезонными причинами. Рассчитывается по формуле: где АФ – атрибутивная фракция в % Yi – месячные средние за многолетний период Y0 – общая средняя 6.2. Статистический анализ 6.2.1. Первоначально выполняется описательная статистика и построение графиков. Необходимо вычислить среднее арифметическое, минимальное и максимальное значения и стандартное отклонение для каждой изучаемой переменной. Можно также вычислить среднее геометрическое, медиану, моду, определенные процентили (например, 10-й и 90-й). Эти величины дают представление о центральной тенденции и разбросе значений переменных временного ряда. Вычисления всех необходимых характеристик временного ряда проводятся в каком-либо стандартном статистическом пакете (SAS, SPSS, STATA, Excel). 6.2.2.Необходимо построить графики всех переменных. Это позволит обнаружить некоторые закономерности в поведении временного ряда: сезонные колебания температуры воздуха или уровней загрязнения атмосферного воздуха; различные долговременные тенденции в данных. Графический анализ позволяет более адекватно подобрать математическую модель для оценки изучаемой связи и помогает обнаружить выпадающие величины, что позволит выяснить, не связаны ли они с ошибками ввода. 6.2.3. При оценке влияния температуры атмосферного воздуха на показатели здоровья населения необходимо определить тепловые и холодовые волны, а также пороговые значения температур, выше или ниже которых наблюдается изменение показателей смертности, заболеваемости, госпитализации населения, обращаемости за скорой медицинской помощью. 6.2.4. Тепловая волна представляет собой пять и более последовательных дней, во время которых среднесуточная температура превышает 97-й процентиль, причем, по крайней мере, в течение трех дней среднесуточная температура должна превышать 99-й процентиль. Для холодовых волн границы определяются на уровне 3-го и 1-го процентилей. Короткой называется волна длительностью от пяти до восьми дней. Волна, которая продолжается 9 дней и более называется длинной. 6.2.5. Пороговые значения температуры воздуха определяются на уровне 90-го и 10-го процентилей всех значений температуры за изучаемый многолетний период. Пороговые значения рассчитываются для среднесуточной, максимальной, среднесуточной эффективной и максимальной эффективной температур. 6.2.6. Одним из вариантов анализа является «анализ временных рядов» (time series analysis) с помощью которого устанавливается зависимость между колебаниями метеорологических показателей (показателей загрязнённости атмосферы) и показателями состояния здоровья населения. 6.2.7. Для оценки наличия взаимосвязи между метеорологическими факторами и количеством смертей, заболеваний, госпитализаций, вызовов скорой медицинской помощи используется корреляционный и многомерный регрессионный анализ для дискретных зависимых переменных (многомерный регрессионный анализ Пуассона, отрицательный биномиальный регрессионный анализ). 6.2.8. При выполнении регрессионного анализа проводится коррекция на многолетние тренды, сезонность, эффект дней недели, выходные и праздничные дни, эпидемии гриппа и другие возможные факторы, способные повлиять на итоговые результаты оценки. Для оценки влияния только температуры воздуха проводят коррекцию на другие метеорологические факторы. 6.2.9. Во многих случаях изменение температуры воздуха вызывает изменение количества случаев смертей, заболеваний, госпитализаций, вызовов скорой медицинской помощи только через определенный интервал времени, т.е. лаг. Лаг 0 означает, что обе переменные относятся к одному и тому же дню (месяцу), лаг 1 означает, что смертность (заболеваемость и др.) анализировалась в день или месяц, следующий за тем, к которому относится метеорологический показатель, и т.д. При анализе суточных данных лаг составляет от 1 до 14 дней, месячных данных – 1-2 месяца. 6.3. Представление результатов 6.3.1. Результатами оценки влияния метеорологических факторов на показатели здоровья населения являются величины относительных рисков; показатели среднего процентного вклада метеорологических факторов в заболеваемость, смертность; показатели прироста или убыли смертности, заболеваемости, госпитализаций и вызовов скорой помощи в процентах при увеличении или уменьшении значений метеорологических параметров (концентраций загрязнителей воздуха) на установленный прирост показателя (например, на 1ºС температуры). Установленные показатели представляются с доверительными интервалами не менее 95%. 6.3.2. Для графического представления данных временных рядов используют линейные и радиальные диаграммы. С помощью линейной диаграммы отражают суточное (месячное) количество смертей, заболеваний, госпитализаций, вызовов скорой медицинской помощи и среднесуточные (среднемесячные) значения температуры воздуха или других метеорологических факторов (или концентраций загрязняющих веществ). При этом целесообразно построение линейных диаграмм с двумя осями, где на одной оси откладывают значения зависимой переменной, а на другой оси – независимой переменной. Радиальные диаграммы используют для отображения месячных показателей здоровья или среднемесячных значений метеорологических факторов. 6.4. Оценка популяционного риска заболеваний, смертей, обращений за скорой медицинской помощью, госпитализации. 6.4.1. При представлении результатов оценки риска по данным эпидемиологических исследований в виде отношения рисков (risk ratio – RR) популяционный риск (R поп ) рассчитывается по формуле: , где Z – абсолютное число случаев смерти заболеваний, госпитализаций, вызовов скорой помощи в исследуемой популяции за аналогичный период по данным многолетних наблюдений. Для более полного представления популяционного риска целесообразно рассчитывать его величину по верхней и нижней границе 95% доверительного интервала RR. 6.4.2. При представлении результатов оценки риска по данным эпидемиологических исследований в виде процентного увеличения или уменьшения количества смертей, заболеваний, госпитализаций и вызовов скорой помощи при изменении температуры на 1ºС популяционный риск (R поп ) рассчитывается по формуле: , где ΔZ – процентное увеличение количества смертей, заболеваний, госпитализаций и вызовов скорой помощи при изменении температуры на 1ºС Δt – отклонение температуры воздуха выше (для волн жары) или ниже (для волн холода) пороговых температур, градусов ºС. VII . Оценка экономического ущерба, связанного с риском для здоровья населения при воздействии погодных и климатических условий 7.1. Расчет экономического ущерба для здоровья населения, связанного с воздействием метеорологических условий, осуществляется в соответствии с МР 5.1.0029-11 «Методические рекомендации к экономической оценке рисков для здоровья населения при воздействии факторов среды обитания». 7.2. Экономический ущерб от рисков для здоровья населения при неблагоприятном воздействии метеорологических условий представляет собой сумму величин ущерба для здоровья, возникающего в связи с повышением уровня смертности и заболеваемости населения (в соответствии с разделом 5 настоящих методических рекомендаций). 7.3.Исходными данными для расчета экономического ущерба являются: - популяционный риск заболеваний, число преждевременных случаев смерти, число обращений за скорой медицинской помощью и госпитализации, рассчитанный в соответствии с разделом 6.4 настоящих методических рекомендаций (дополнительные случаи, связанные с метеорологическими факторами); - экономические показатели для расчета ущерба здоровью населения в результате климаточувствительных заболеваний, полученных в соответствии с разделом 5.3 настоящих методических рекомендаций. 7.4. Результаты расчета экономического ущерба для здоровья населения, выполненные по исследуемым климаточувствительным заболеваниям среди различных групп населения повышенного риска, используются для разработки стратегии и мер по адаптации к изменению климатических условий и выбора сценария управления риском с учетом имеющихся ресурсов (в стоимостном выражении) и возможности их реализации. В этом случае используются оценки предотвращенного ущерба для здоровья населения, выполняемые для различных вариантов (сценариев) реализации адаптационных мер. 7.5. Для выбора и оценки стратегии и мер по адаптации к изменению климатических условий могут быть использованы методические подходы, изложенные в МР 5.1.0030-11 «Методические рекомендации к экономической оценке и обоснованию решений в области управления риском для здоровья населения при воздействии факторов среды обитания». Приложение 1 к МР 2.1.10.0057-12 Пример базы данных по метеорологическим факторам и количеству вызовов скорой медицинской помощи в городе N Дата Год Ме-сяц День Тип дня Грипп Т сред. Т макс. Н сред. Н макс. J00-99 0-17 лет J00-99 18-59 лет J00-99 60+ лет S00 -T98 0-17 лет S00 -T98 18-59 лет S00 -T98 60+ лет 01.01.2000 2000 1 6 1 0 -21,1 -23,2 80 84 33 25 8 12 34 2 02.01.2000 2000 1 7 1 0 -24,9 -26,2 80 81 22 19 4 6 20 1 03.01.2000 2000 1 1 1 0 -22,9 -26,0 81 83 28 33 19 7 14 2 04.01.2000 2000 1 2 1 0 -19,3 -21,7 83 85 21 31 15 3 28 5 05.01.2000 2000 1 3 0 0 -16,5 -18,5 84 84 25 14 10 7 27 4 06.01.2000 2000 1 4 0 0 -22,3 -26,0 83 85 21 10 11 6 19 5 07.01.2000 2000 1 5 1 0 -15,6 -23,9 87 95 10 10 5 6 21 6 08.01.2000 2000 1 6 1 0 -1,6 -5,4 94 96 13 15 7 5 16 6 09.01.2000 2000 1 7 1 0 2,0 0,3 84 96 15 20 7 8 32 3 10.01.2000 2000 1 1 0 0 1,1 -0,2 80 84 32 20 8 12 36 6 …. …. …. …. …. …. …. …. …. …. …. …. …. …. …. …. 22.12.2008 2008 12 1 0 0 -5,3 -6,3 83 87 29 9 10 4 21 3 23.12.2008 2008 12 2 0 0 -6,7 -7,7 83 86 21 8 5 8 25 2 24.12.2008 2008 12 3 0 0 -7,1 -9,2 90 93 22 11 8 5 20 6 25.12.2008 2008 12 4 0 0 -6,2 -6,5 91 94 18 13 8 3 26 8 26.12.2008 2008 12 5 0 0 -6,3 -8,3 88 96 18 15 10 7 24 7 27.12.2008 2008 12 6 1 0 -,4 -2,9 89 96 27 18 8 9 37 8 28.12.2008 2008 12 7 1 0 -2,4 -3,0 82 95 43 16 6 4 60 4 29.12.2008 2008 12 1 0 0 -,4 -2,8 95 96 27 16 6 5 31 13 30.12.2008 2008 12 2 0 0 -,8 -2,2 90 94 35 8 6 8 33 13 31.12.2008 2008 12 3 0 0 -,4 -1,0 91 95 26 11 5 6 31 6 Тип дня: 0 – рабочий, 1 – выходной или праздничный; Т сред. – среднесуточная температура, Т макс. – максимальная температура, Н сред. – среднесуточная влажность, Н макс. – максимальная влажность; J00-99 - болезни органов дыхания, S00 -T98 – травмы, отравления и другие последствия внешних причин Приложение 2 к МР 2.1.10.0057-12 Оценка воздействия климатических изменений на здоровье населения в различных регионах России 1. Оценка влияния температуры воздуха, волн жары и холода на смертность населения 1.1. Оценка влияния температуры воздуха на смертность населения Москвы в 2000-2006 годах [14]. Исследование влияния температуры воздуха на смертность населения Москвы проведено с использованием метода временных рядов по данным ежедневной смертности и температуры воздуха за 2000-2006 годы. Данные о суточном количестве смертей, распределенных по полу, возрасту и причинам смерти получены из базы данных Росстата. Использование метода временных рядов показало, что температурная кривая смертности аппроксимировалась V-образной функцией с двумя линейными участками, соответственно ниже и выше точки минимума температурной кривой. Такая аппроксимация позволяет вычислить коэффициенты линейной регрессии, которые интерпретируются следующим образом: коэффициент для холодных температур имеет смысл относительного увеличения смертности в среднем на каждый градус снижения температуры ниже точки оптимума (+18ºС для Москвы), соответственно коэффициент для жарких температур интерпретируется как относительное увеличение смертности в среднем на каждый градус повышения температуры выше точки оптимума. Результаты регрессионного анализа с использованием пуассоновской модели приведены в таблицах 1 и 2. Приведенные коэффициенты справедливы в среднем для каждого температурного диапазона, потому что в общем случае изучаемая зависимость, конечно, нелинейна. Для всех изученных причин смерти (кроме хронических болезней нижних дыхательных путей у лиц в возрастной группе 75+) установлена связь между температурой и смертностью. Сравнение показателей смертности в разных возрастных группах жителей Москвы показало, что для всех причин смерти угол наклона регрессионной прямой для возрастной группы 75 лет и старше круче, чем для группы «все возраста». Таблица 1 Эффект холодных температур (-10ºC<T<18ºC) на изменение суточной смертности в Москве [14] Причина смерти Возраст-ная группа Относительное изменение смертности на 1ºС Лаг, дни % 95% ДИ Все, кроме внешних причин Все -0,49 -0,53; -0,45 3 75+ -0,64 -0,71; -0,59 3 Ишемическая болезнь сердца (инфаркт, стенокардия) Все -0,57 -0,63; -0,51 3 75+ -0,69 -0,77; -0,61 3 Цереброваскулярные заболевания (инсульты) Все -0,78 -0,86; -0,70 6 75+ -0,92 -1,02; -0,82 6 Хронические заболевания нижних дыхательных путей Все -1,31 -1,75; -0,87 4 75+ -1,21 -1,93; -0,49 5 Таблица 2 Влияние жары (T>18ºC) на изменение суточной смертности в Москве [14] Причина смерти Возраст-ная группа Относительное изменение смертности на 1ºС Лаг, дни % 95% ДИ Все, кроме внешних причин Все 2,8 2,0; 3,6 0 75+ 3,3 2,1; 4,5 1 Ишемическая болезнь сердца (инфаркт, стенокардия) Все 2,7 1,7; 3,7 0 75+ 3,1 1,7; 4,5 0 Цереброваскулярные заболевания (инсульты) Все 4,7 3,5; 5,9 1 75+ 5,3 3,7; 6,9 1 Хронические заболевания нижних дыхательных путей Все 8,7 0,7; 16,7 0 75+ – – – Влияние волн жары и холода на смертность изучалось методом анализа независимых выборок из временных рядов суточной смертности. На основе анализа многолетних распределений среднесуточных температур даны формальные определения температурных волн в Москве. В частности, волна жары состоит из пяти или более последовательных дней со среднесуточной температурой выше +22,7ºС, аналогично холодовая волна – как минимум из девяти последовательных дней со среднесуточной температурой ниже -14,4ºС. За указанный период времени исключительно жаркими в Москве были июль 2001 года и июль 2002 года. В июле 2001 года столица пережила волну жары, во время которой среднесуточные температуры превышали порог в 25ºС в течение девяти последовательных дней (при средней многолетней «норме» три дня в год). В максимуме этой волны суточная смертность превысила среднее многолетнее значение смертности для июля на 93%. Волна 2001 года привела к четко выраженному и статистически значимому всплеску смертности во всех возрастных группах по всем рассмотренным причинам смерти. Абсолютная дополнительная смертность во время волны жары в 2001 году составила 1177 случаев, в 2002 году – 283 случая. Эффект аномальных метеорологических условий другого типа – «холодовой волны» – наглядно демонстрирует ситуация января-февраля 2006 года, когда в Москве аномально низкие температуры наблюдались в течение 26 дней. Поскольку эта волна холода состояла из двух эпизодов, разделенных краткосрочным потеплением, ее эффект оказался статистически значим только для пожилых людей. Оценено изменение смертности, обусловленное потеплением между «базовым» периодом 1980 – 1999 гг. и периодом 2000 – 2005 гг., в сумме за шесть лет исследуемого периода. Согласно расчетам, снижение смертности в результате потепления климата в зимний период времени в Москве с 2000 по 2005 годы составило примерно 590 смертей в год, но за этот же период в результате увеличения среднемесячных температур летом дополнительное количество смертей составило 420 случаев. Таким образом, в сумме положительное и отрицательное воздействия потепления климата на смертность почти компенсируют друг друга: результирующая дополнительная смертность ΔМ = 420 – 590 = -170; 95% ДИ ΔМ составил (-291; -49) смертей в год. Суммарный прирост смертности оказался отрицательным, т.е. потепление климата в итоге ненамного снижает смертность. Результирующий эффект довольно мал и практически сравним с погрешностью самих вычислений. 1.2. Волны жары, качество атмосферного воздуха и смертность населения Москвы в 2000-2006 годы [22]. Анализ данных о концентрациях загрязняющих веществ в атмосферном воздухе Москвы показал их увеличение в наиболее жаркие дни, т.е. краткосрочные колебания концентрации «повторяют» колебания температуры. Это следует из анализа простых регрессионных уравнений между загрязнением и температурой. Наиболее сильная зависимость концентрации загрязнения от температуры наблюдалась при нулевом лаге, т.е. от температуры того же дня. Температуры соседних дней также не являются статистически независимыми переменными. В день с максимальной температурой наблюдалась максимальная концентрация O3 , а максимальные концентрации NO2 и PM10 были зафиксированы днем ранее, причем концентрация РМ10 достигла экстремально высокого значения µ+2σ (среднегодовое плюс два стандартных отклонения). В зимний период времени концентрации взвешенных частиц в атмосферном воздухе увеличивалась по мере снижения температуры. Наиболее сильная зависимость наблюдалась от температуры предыдущего дня (t -тест для коэффициента линейной регрессии t = 7,1). Связь между температурой и концентрацией РМ10 носит длительный характер – 6 дней и более, что объясняется худшими условиями рассеяния в зимнее время. Концентрация озона, наоборот, увеличивалась с повышением температуры в зимний период. На основе анализа данных о суточной смертности от всех естественных причин, а также от таких «климатозависимых причин», как ишемическая болезнь сердца (ИБС, МКБ-10: I20-I25) и цереброваскулярные заболевания, в основном инсульты (МКБ-10: I60-I69) изучена зависимость смертности от загрязнения с помощью Пуассоновской регрессионной модели: , где E (Mt ) – ожидаемое число смертей в день t ; PM 10t и O3t – cреднесуточные концентрации PM10 и озона, усредненные по данным станций ФГУП «Мосэкомониторинг» в Москве в день t ; St (Temp ,6) – множитель, введенный в модель для учета поправки на зависимость смертности от температуры воздуха в день смерти, он моделировался кубическим сплайном – интерполяцией температуры с шестью степенями свободы; Yt (M ,DF ) – множитель, введенный в модель для учета поправки на все медленные колебания смертности: сезонные, многолетние, зимние эпидемии и др., характерный временной период которых больше корреляционного радиуса зависимости загрязнения от температуры ±(6-8) дней, т.е. примерно две недели). Этот множитель аппроксимировался непараметрической сглаживающей функцией ежедневной смертности М с числом степеней свободы DF = 182 (т.е. зависящей от значений ежедневной смертности за предыдущие полгода). Результаты регрессионного анализа модели, линейной по концентрациям PM10 и O3 , приведены в таблицах 4 и 5. Регрессионные коэффициенты показывают относительный прирост смертности, соответствующий приросту концентрации загрязняющего вещества на каждые 10 мкг/м3 в среднем для всего диапазона концентраций, наблюдавшихся в Москве в течение периода исследования. Таблица 4 Относительные приросты смертности населения Москвы при увеличении на 10 мкг/м3 PM10 [22] Причина смерти Возрастная группа Прирост смертности р % 95% ДИ Все естественные причины Все 0,33 0,09; 0,57 0,006 75+ 0,55 0,21; 0,89 0,002 ИБС Все 0,66 0,30; 1,02 <0,001 75+ 0,81 0,31; 1,31 0,001 Цереброваскулярные заболевания Все 0,48 0,02; 0,94 0,035 75+ 0,72 0,14; 1,30 0,014 Таблица 5 Относительные приросты смертности населения Москвы при увеличении на 10 мкг/м3 О3 [22] Причина смерти Возрастная группа Прирост смертности р % 95% ДИ Все естественные причины Все 1,09 0,71; 1,47 <0,001 75+ 1,24 0,68; 1,80 < 0,001 ИБС Все 1,61 1,01; 2,21 < 0,001 75+ 1,88 1,08; 2,68 < 0,001 Цереброваскулярные заболевания Все 1,28 0,54; 2,02 0,001 75+ 1,25 0,31; 2,19 0,008 Наиболее сильная зависимость смертности от уровня загрязнения PM10 получена при нулевом лаге, а от уровня загрязнения О3 – при кумулятивном лаге 0-1. Для всех изученных показателей смертности получены статистически значимые риски воздействия как PM10 , так и озоном. В возрастной группе 75 лет и старше риски выше, чем в группе «все возраста». Высокий вклад этой возрастной группы в общую смертность населения, объясняет, что прирост общей смертности, вызванный загрязнением, связан с приростом смертности среди пожилых людей. Дополнительная смертность связана в основном с увеличением смертности от заболеваний сердечно-сосудистой системы, обусловленной преимущественно воздействием РМ10 . Изучено сочетанное действие двух загрязнителей (суммы PM10 и О3 ) на смертность. Для этого вычислены риски PM10 отдельно для выборки дней, в которые концентрация О3 превышала 90-й процентиль распределения среднесуточных концентраций О3 за весь период исследования (уже в рамках одномерной по загрязнению модели, в которой смертность зависит только от PM10 ). Эти «скорректированные на высокий уровень О3 » риски PM10 сравнивались с исходными, вычисленными для всех дней периода исследования в рамках одномерной регрессионной модели. Разница была значительной, что подтверждает модифицирующий эффект высоких уровней О3 на риск PM10 . В присутствии высоких уровней О3 риски воздействия PM10 возрастают примерно в 3 раза, но высокие концентрации PM10 не приводят к увеличению влияния О3 на смертность. 1.3. Оценка влияния температуры воздуха на смертность населения Москвы летом 2010 года [15]. В работах по климату порогом аномальности температуры считается ее превышение на 5°С, поэтому для предварительной оценки последствий жары был использован именно этот показатель. В июле – августе 2010 года протяженность волны жары в Москве со среднесуточной температурой выше среднемноголетней на 5 °С составила 45 дней. Число температурных рекордов, а именно, дней с максимальной температурой достигло в июле 10 дней и в августе 9 дней. Антициклон в московском регионе препятствовал рассеиванию загрязняющих веществ в атмосферном воздухе, и дополнительное их количество поступило в результате пожаров. Наиболее высокие концентрации загрязняющих веществ присутствовали в атмосферном воздухе Москвы в период с 14 июля по 19 августа в условиях высокого атмосферного давления и температурной инверсии. При среднем уровне загрязнения атмосферного воздуха в июле совпадают пиковые значения концентраций и температуры, но в августе, когда в результате мощных пожаров содержание наиболее токсичной мелкодисперсной пыли размером менее 10 мкм (РМ10 ), из-за которой образовалась мгла, резко возросло (до 15 ПДКсс ), температура приземного слоя несколько снизилась. Максимальные концентрации моноксида углерода достигали 30 мг/м3 , РМ10 – 1 500 мкг/м3 , среднесуточные концентрации РМ10 во время пожаров с 4 по 9 августа находились в пределах 431–906 мкг/м3 , превышая российские ПДКсс (60 мкг/м3 ) в 7,2–15,1 раза. Концентрации в атмосферном воздухе формальдегида, этилбензола, бензола, толуола, стирола и некоторых других органических веществ также были превышены (до 8 раз выше ПДК). Во время аномальной жары 2010 года смертность населения Москвы выросла по всем крупным классам причин смерти на 11 тысяч случаев по сравнению с июлем – августом 2009 года, причем в августе во время пожаров произошел более резкий ее рост от заболеваний органов дыхания (табл. 3), значительный рост от инфекционных и паразитарных заболеваний (на 61,5%), новообразований (на 70,2% ), от внешних причин (на 52,9%). Из внешних причин в наибольшей степени выросла смертность от суицидов в июле на 63 случая (101,6%) и в августе на 38 случаев (52,1%). Таблица 3 Волна аномальной жары и смертность в Москве в 2010 году [15] Показатель Июль Август Всего Число дней с температурой выше многолетней среднемесячной на 5°С 27 18 45 Дополнительная смертность в 2010 году по сравнению с 2009 годом, абс. (%) +4 824 (50,7) +6 111 (68,6) +10 935 (59,6) в том числе от: - болезней системы кровообращения, % - болезней органов дыхания, % - инфекционных болезней, % - новообразований, % - внешних причин, % 51,5 59,1 56,3 58,8 48,0 66,1 110,1 66,7 81,6 57,8 58,8 84,5 61,5 70,2 52,9 По оперативным данным Управления ЗАГС известно, что в июле происходило постепенное нарастание смертности со второй недели месяца. В дни максимальной температуры число случаев смерти возрастало вдвое, причем увеличивалась смертность в старшей возрастной группе. В такие дни на 32% увеличилось количество выездов бригад скорой медицинской помощи по поводу заболеваний системы кровообращения. Число обращений за скорой медицинской помощью в августе было выше, чем в июле, на 31 %, причем увеличилась доля вызовов по поводу заболеваний органов дыхания. В сентябре 2010 года уровень смертности был уже несколько ниже уровня сентября 2009 года, т. е. проявился «эффект жатвы», который захватил и октябрь. В ноябре 2010 г. в Москве было зарегистрировано на 832 случаев смерти (на 8,4%) меньше по сравнению с ноябрем 2009 года. 1.4. Оценка влияния температуры воздуха на смертность населения Воронежа летом 2010 года Целью исследования явилось изучение влияния температуры воздуха на смертность населения Воронежа в период чрезвычайной ситуации, связанной с пожарами летом 2010 года. Летом 2010 года в Воронеже регистрировалась аномально высокая температура воздуха. Волна жары состояла из 29-ти последовательных дней со среднесуточной температурой от 26,0 до 31,7°С, при среднегодовой 18,4-19,7°С. Максимальная температура достигла 40,1°С 29 июля 2010 года. Данные о среднесуточной и максимальной за сутки температуре воздуха предоставлялись ГУ «Воронежский областной центр по гидрометеорологии и мониторингу окружающей среды». Информация о суточной регистрации количества смертей получена в территориальном органе Федеральной службы государственной статистики по Воронежской области. Для выполнения анализа влияния температуры воздуха на смертность населения создана база данных по среднесуточной и максимальной за сутки температуре воздуха и суточному количеству смертей от болезней органов кровообращения, болезней органов дыхания и всех причин. Оценка связи между температурой воздуха и ежедневными случаями смерти проводилась с помощью метода временных рядов. Влияние экстремальной температуры воздуха на смертность населения изучалось с лагами 0, 1, 2, 3 дня. На рис.1 представлена суточная динамика количества смертей от болезней органов кровообращения и максимальной температуры в Воронеже с лагом в 2 дня. Рис. 1. Суточная динамика количества смертей от болезней органов кровообращения и максимальной температуры в Воронеже с 21 июля по 20 августа 2010 года (лаг 2 дня) Корреляционный анализ свидетельствует о статистически значимой положительной зависимости между температурой воздуха и смертностью населения от всех причин и болезней органов кровообращения с лагом 0, 1, 2, 3 дня и смертностью от болезней органов дыхания с лагом 3 дня (табл.6). Таблица 6 Коэффициенты корреляции Пирсона между смертностью и максимальной температурой воздуха в Воронеже летом 2010 года Лаг смертность (все причины) смертность от болезней системы кровообращения смертность от болезней органов дыхания r 95% ДИ r 95% ДИ r 95% ДИ 0 дней 0,76* 0,62-0,91 0,77* 0,63-0,91 0,42 0,13-0,71 1 день 0,73* 0,56-0,90 0,72* 0,55-0,89 0,35 0,29-0,77 2 дня 0,73* 0,56-0,90 0,70* 0,53-0,88 0,30 0,24-0,74 3 дня 0,65* 0,44-0,86 0,66* 0,47-0,86 0,45* 0,16-0,74 r − коэффициент корреляции; * р<0,05 (уровень статистической значимости) В таблице 7 и на рисунке 2 представлены результаты регрессионного анализа, отражающего зависимость смертности населения от максимальной температуры воздуха (наиболее статистически значимые коэффициенты корреляции получены при лаге «0»). Таблица 7 Коэффициенты линейной регрессии между смертностью и максимальной температурой воздуха в г. Воронеж летом 2010 года Лаг смертность (все причины) смертность от болезней системы кровообращения b p b p 0 дней 2,26 <0,05 2,60 <0,05 1 день 3,38 <0,05 2,60 <0,05 2 дня 3,55 <0,05 2,27 <0,05 3 дня 3,36 <0,05 2,67 <0,05 b − коэффициент регрессии; р − уровень статистической значимости Рис. 2. Регрессионный анализ оценки зависимости между максимальной температурой воздуха и количеством смертей от болезней органов кровообращения в Воронеже (лаг 0 дней) Регрессионные коэффициенты свидетельствуют, что с ростом температуры воздуха на 1 градус количество случаев смерти увеличивается на 3%. Необходимо отметить, что зависимость смертности населения от аномально высокой температуры воздуха не является абсолютно очевидной, поскольку дополнительно присутствовало загрязнение атмосферного воздуха в период пожаров, что внесло неопределенность в интерпретацию полученных данных. 1.5. Оценка влияния температуры и загрязнения воздуха на смертность населения Свердловской области летом 2010 года Целью исследования была оценка влияния на смертность населения, проживающего в промышленно развитых городах Свердловской области (Екатеринбург, Нижний Тагил и Верхняя Пышма), факторов риска (высокая температура, лесные пожары и инверсионные процессы в атмосферном воздухе), связанных с действием на территории Европейской части Российской Федерации аномально стабильного антициклона в период с его зарождения (май 2010 года) до распада (август 2010 года). Задачи исследования: 1. Выполнить сравнительную оценку уровня загрязнения атмосферного воздуха с мая по август 2010 года (период действия антициклона) с аналогичным периодом 2009 года в промышленно развитых городах Свердловской области. 2. Оценить влияние аномально стабильного антициклона на условия рассеивания промышленных и автотранспортных выбросов в городах с различным уровнем техногенного загрязнения. 3. Изучить влияние высокой температуры воздуха на смертность населения в промышленных городах Верхняя Пышма, Нижний Тагил и Екатеринбург в период действия антициклона. 4. Оценить влияние факторов риска, связанных с природным загрязнением атмосферного воздуха в результате лесных пожаров в Висимском заповеднике (30-40 км от города Нижний Тагил) и Шутовских болотах (25-40 км от городов Верхняя Пышма и Екатеринбург) и техногенным загрязнением, на здоровье населения. 5. Оценить риск и экономический ущерб для здоровья населения в промышленных городах Свердловской области в связи с действием высоких температур и неблагоприятных условий рассеивания техногенных и природных выбросов в период действия антициклона. Данные о концентрациях загрязняющих веществ и температуре воздуха были предоставлены СОГУ «Центр экологического мониторинга и контроля» Министерства природных ресурсов Свердловской области. Исходной информацией для оценки уровня загрязнения атмосферного воздуха явились ежедневные измерения концентраций пылевых частиц с аэродинамическим диаметром до 10 мкм (РМ10 ), диоксида азота, диоксида серы и оксида углерода, мониторинг которых в городах Верхняя Пышма, Нижний Тагил и Екатеринбург проводится на автоматических станциях контроля качества атмосферного воздуха «СКАТ». Среднее ежемесячное число измерений по каждому веществу составило от 1440 до 2220 измерений. Информация о количестве смертей за сутки была предоставлена ГОУЗ «Медицинский информационно-аналитический центр» Министерства здравоохранения Свердловской области. В анализ включены случаи смерти населения от всех причин, а также отдельно от болезней органов дыхания и болезней системы кровообращения. Из общей смертности были исключены травмы, отравления и другие последствия внешних причин. Результаты исследования. По задаче 1. В результате действия аномально стабильного антициклона предельно допустимые концентрации загрязняющих веществ техногенного (выбросы автотранспорта и промышленные выбросы) и природного происхождения (лесные пожары) в промышленно развитых городах Свердловской области были превышены до 2-5 раз в июле и августе 2010 года относительно аналогичных показателей 2009 года. В период с мая по июнь 2010 года значительных превышений предельно допустимых концентраций не зафиксировано (табл. 8, 9 и 10). Таблица 8 Концентрации загрязняющих веществ в атмосферном воздухе Верхней Пышмы за период с мая по август 2009 г. и 2010 г. (мг/м3 )       Вещество Среднесуточная концентрация Максимальная разовая концентрация ПДК май июнь июль август ПДК май июнь июль август 2009 год РМ10 0,06 0,02 0,0088 0,0071 0,01 0,3 0,29 0,1 0,1 0,09 NO2 0,04 0,02 0,01 0,01 0,012 0,2 0,11 0,08 0,06 0,05 SO2 0,05 0,02 0,0084 0,005 0,005 0,5 0,56 0,19 0,08 0,62 CO 3,0 0,12 0,17 0,13 0,0155 5,0 2,51 2,14 1,35 2,64 2010 год РМ10 0,06 0,007 0,002 0,007 0,04 0,3 0,09 0,04 0,12 0,18 NO2 0,04 0,04 0,027 0,038