Azot dioksit (NO₂) ve ozon (O₃), atmosfer kimyasında birbirine güçlü biçimde bağlı iki gaz kirleticidir. Bu bağlantı, uzun süreli hava kirliliği maruziyetinin sağlık etkilerini inceleyen epidemiyolojik modeller açısından önemli bir istatistik problemi oluşturur. Bir çalışma yalnız NO₂'yi modele dahil edip O₃'ü dışarıda bırakırsa, iki kirletici arasındaki korelasyon nedeniyle NO₂ katsayısı O₃ etkisinin bir bölümünü de taşıyabilir. Buna Atlanmış Değişken Yanlılığı — Omitted Variable Bias (OVB) — denir.

İki kirleticiyi aynı modele dahil etmek teorik olarak bu sorunu azaltabilir. Ancak NO₂ ile O₃'ün yüksek korelasyonu bu kez multicollinearity ve variance inflation oluşturarak standart hataları büyütebilir. Ayrıca literatürde gerçek iki-kirleticili kohort sonuçlarının sayısı tek-kirleticili sonuçlara göre çok daha azdır.

Bu çalışma, tek-kirleticili çalışmalardaki sonuçları atmosferik NO₂–O₃ korelasyonu ve standart sapma oranları kullanılarak iki-kirleticili eşdeğerlere dönüştürmeye yönelik bir OVB düzeltme yaklaşımı uygulamaktadır. Araştırmacılar daha sonra bu düzeltilmiş tahminleri mevcut iki-kirleticili çalışmalarla random-effects meta-analizinde birleştirmiştir.

OVB-düzeltilmiş analizde uzun dönem NO₂ için her 10 ppb artış başına tüm nedenlere bağlı mortalite risk oranı RR = 1,076 ve %95 güven aralığı 1,040–1,112 olarak raporlanmıştır. O₃ için düzeltilmiş RR = 1,015 ve %95 güven aralığı 0,984–1,048'dir. Buna göre NO₂ ile mortalite arasındaki pozitif ilişki O₃ düzeltmesinden sonra devam ederken O₃ için bağımsız sonuç istatistiksel olarak anlamlı değildir.

Bu bulgular gözlemsel kohort çalışmalarının meta-analizine ve matematiksel bir bias-düzeltme modeline dayanır. Sonuçlar NO₂'nin doğrudan nedensel olarak belirli sayıda ölüme yol açtığını tek başına kanıtlamaz.

NO₂ ve O₃ neden birbirinden kolayca ayrılamıyor?

Troposferde NO, NO₂ ve O₃ hızlı fotokimyasal reaksiyonlarla birbirine bağlıdır. Trafiğin yoğun olduğu kent merkezlerinde NO'nun O₃'ü tüketmesi nedeniyle NO₂ ve O₃ çoğu zaman negatif korelasyon gösterebilir. Kent dışı veya farklı fotokimyasal rejimlerde ise iki kirletici pozitif korelasyon gösterebilir.

Bu nedenle NO₂–O₃ ilişkisi yalnız kuvvet açısından değil, işaret açısından da coğrafyaya göre değişebilir.

Bir ortak kirletici regresyon modelinde atlandığında korelasyonun yönü, hedef kirleticinin katsayısında oluşan yanlılığın yönünü de değiştirebilir.

Atlanmış Değişken Yanlılığı NO₂ Tahminini Nasıl Değiştirebilir?

Çalışmada kullanılan temel OVB ilişkisi:

\[ \beta_1^* = \beta_1 + \rho \left( \frac{\sigma_2}{\sigma_1} \right) \beta_2 \]

şeklindedir.

Burada \(\beta_1^*\) tek-kirleticili modelde gözlenen NO₂ benzeri hedef katsayıyı, \(\beta_1\) bağımsız hedef-kirletici etkisini, \(\beta_2\) modele dahil edilmemiş ortak kirleticinin etkisini, \(\rho\) iki kirletici arasındaki korelasyonu ve \(\sigma_1,\sigma_2\) kirletici dağılımlarının standart sapmalarını temsil eder.

Örneğin NO₂ ile O₃ negatif korelasyonluysa ve O₃'ün sağlık etkisi pozitifse, O₃'ün modelden çıkarılması NO₂ katsayısını aşağı yönlü biaslayabilir. Pozitif korelasyon bulunan bölgelerde ise aynı mekanizma ters yönde çalışabilir.

Bu nedenle tek bir küresel NO₂–O₃ korelasyonu kullanmak yerine araştırmacılar çalışma bölgelerine göre atmosferik korelasyon yapısını eşleştirmeye çalışmıştır.

Korelasyon Verisi Nereden Geldi?

Birçok kohort çalışması kendi NO₂–O₃ korelasyon katsayısını veya standart sapma oranlarını raporlamadığı için araştırmacılar dış atmosfer veri kaynaklarına başvurmuştur.

Başlıca kaynaklar:

• Copernicus Atmosphere Monitoring Service — CAMS reanalysis,
• Tropospheric Ozone Assessment Report — TOAR,
• China National Environmental Monitoring Centre — CNEMC.

CAMS verileri yaklaşık 0,75° × 0,75° çözünürlükte küresel atmosfer bileşimini temsil ederken TOAR ve CNEMC yüzey gözlem istasyonu verileriyle karşılaştırma sağlamıştır.

Bu yaklaşımın güçlü yanı, yayımlanmış kohortların eksik korelasyon bilgisini küresel atmosfer verileriyle tamamlayabilmesidir. Buna karşılık dış grid verisinin her kohorttaki gerçek bireysel maruziyet korelasyonunu eksiksiz temsil ettiği varsayılmamalıdır.

Literatür Taramasında Hangi Çalışmalar Kullanıldı?

PubMed, Embase ve Web of Science taramalarında önceki benchmark meta-analizler güncellenmiş ve Ocak 2022–Aralık 2024 arasında yayımlanan çalışmalar ayrıca araştırılmıştır.

İlk taramada O₃ için 2.711 ve NO₂ için 886 yayın elde edilmiştir.

Dahil edilme kriterleri:

• genel popülasyona dayalı orijinal kohort çalışması,
• uzun dönem dış ortam NO₂ veya O₃ maruziyeti,
• tüm nedenlere bağlı mortalite sonucu,
• RR veya HR ve %95 güven aralığının raporlanması,
• kirletici konsantrasyonunun sürekli değişken olarak analiz edilmesi.

Kısa dönem zaman serileri, derlemeler, hayvan deneyleri ve dönüştürülebilir risk sonucu bulunmayan çalışmalar dışlanmıştır.

Sonuç bölümünde 80 kohort çalışmasının dahil edildiği; 66'sının NO₂ ve 46'sının O₃ sonuçları raporladığı belirtilmektedir.

Geleneksel Meta-Analizler Ne Gösterdi?

Tek ve iki-kirleticili sonuçları ayırmadan birleştiren naïve analizde NO₂ için yaklaşık:

\[ RR=1.086 \]

elde edilmiştir.

O₃ için:

\[ RR=1.019 \quad (95\%\,GA:0.986-1.052) \]

raporlanmıştır.

Her iki analizde heterojenlik son derece yüksektir:

\[ I^2_{NO_2}=99.3\% \]

ve:

\[ I^2_{O_3}=99.7\%. \]

Bu kadar yüksek heterojenlik farklı bölgelerin, maruziyet yöntemlerinin ve popülasyonların tek bir ortak etki büyüklüğü altında tam olarak homojen olmadığını gösterir.

OVB Düzeltmesi NO₂ Sonucunu Nasıl Değiştirdi?

Yalnız doğrudan iki-kirleticili sonuçlar kullanıldığında NO₂ için sekiz çalışmadan:

\[ RR=1.120 \quad (95\%\,GA\approx1.057-1.187) \]

elde edilmiştir.

OVB-düzeltilmiş tek-kirleticili sonuçlar da analize eklendiğinde kaynak:

\[ \boxed{ RR=1.076 \quad (95\%\,GA:1.040-1.112) } \]

değerini raporlamaktadır.

Bu sonuç, çalışmanın varsayımları altında, uzun dönem NO₂ maruziyetindeki her 10 ppb artışın tüm nedenlere bağlı mortalite riskinde yaklaşık %7,6 daha yüksek göreli riskle ilişkili olduğunu gösterir.

Bu %7,6 değeri mutlak ölüm olasılığında 7,6 yüzde puanlık artış değildir ve tek başına nedensel etki olarak yorumlanmamalıdır.

O₃ İçin Sonuç Neden Farklı?

OVB-düzeltilmiş O₃ tahmini:

\[ RR=1.015 \quad (95\%\,GA:0.984-1.048) \]

olarak bulunmuştur.

Güven aralığı 1'i içerdiği için bu sonuç kaynak tarafından istatistiksel olarak anlamlı kabul edilmemektedir.

Bu durum O₃'ün hiçbir sağlık etkisinin bulunmadığını göstermez. Çalışma yalnız belirli uzun dönem O₃ maruziyet metriği ve tüm nedenlere bağlı mortalite sonucu için, NO₂ confounding düzeltmesinden sonra açık bir bağımsız ilişki saptayamamıştır.

“2,643 Milyar Kişi Etkilendi” İfadesi Ne Anlama Geliyor?

Çalışmanın mekânsal analizinde araştırmacılar, birleşik NO₂–O₃ etkisinin model tarafından istatistiksel olarak anlamlı kabul edildiği dünya grid hücrelerinde yaşayan nüfusu hesaplamıştır.

Geleneksel iki-kirleticili model yaklaşık 109 milyon kişiyi kapsayan hücreleri anlamlı olarak belirleyebilirken OVB-düzeltilmiş modelde bu değer yaklaşık 2,643 milyar kişiye yükselmiştir.

Bu sayı ölüm sayısı değildir.

Bu sayı, 2,643 milyar kişinin her birinin kişisel olarak hava kirliliği nedeniyle zarar gördüğünün klinik kanıtı da değildir.

Doğru tanım, modelin mekânsal anlamlılık kriterini sağlayan grid hücrelerinde yaşayan toplam nüfustur.

Kaynak Asya için yaklaşık 2,193 milyar, Avrupa için 169 milyon ve Kuzey Amerika için 105 milyonluk model kapsamı bildirmektedir.

Çalışmanın Yöntemi ve Bulguları

Literatür standardizasyonu

Kirletici birimleri ppbv'ye dönüştürülmüş ve risk oranları 10 ppbv artış başına standardize edilmiştir.

O₃ maruziyeti farklı çalışmalar arasında karşılaştırılabilir hale getirilmek için 6 aylık sıcak sezon günlük maksimum 8 saatlik konsantrasyon — 6mDMA8 — metriğine dönüştürülmüştür.

Risk-of-bias değerlendirmesi

WHO'nun Risk of Bias aracı kullanılmış ve özellikle co-pollutant confounding kontrolü değerlendirilmiştir.

Kaynağa göre çalışmaların önemli bölümünde bu alanda orta veya yüksek risk bulunmuştur.

Meta-analiz

Yüksek çalışma heterojenliği nedeniyle random-effects model kullanılmıştır.

NO₂ naïve sonuç

Yaklaşık RR = 1,086 ve çok yüksek heterojenlik raporlanmıştır.

O₃ naïve sonuç

RR = 1,019 (0,986–1,052) ve I² = %99,7 raporlanmıştır.

NO₂ doğrudan iki-kirleticili sonuç

RR = 1,120; yaklaşık %95 GA = 1,057–1,187.

NO₂ OVB-düzeltilmiş sonuç

RR = 1,076; %95 GA = 1,040–1,112.

O₃ OVB-düzeltilmiş sonuç

RR = 1,015; %95 GA = 0,984–1,048.

İstatistiksel hassasiyet

NO₂ için güven aralığının iki-kirleticili modele göre yaklaşık %44 daraldığı raporlanmıştır.

Mekânsal korelasyon

Kentsel çekirdeklerde NO₂–O₃ korelasyonu çoğunlukla negatif, bazı downwind ve farklı fotokimyasal rejimlerde pozitif bulunmuştur.

Küresel nüfus kapsamı

OVB-düzeltilmiş model, anlamlı birleşik etki saptanan grid hücrelerinde yaklaşık 2,643 milyar kişinin yaşadığını hesaplamıştır; geleneksel modelde karşılık gelen değer yaklaşık 109 milyondur.

Çalışmanın desteklediği sonuçlar

• NO₂–O₃ korelasyonunun coğrafyaya göre yön değiştirmesi tek-kirleticili modellerde farklı yönlerde OVB oluşturabilir.
• OVB düzeltmesi tek-kirleticili çalışmaların bilgisini iki-kirleticili meta-analize dahil etmeye yönelik bir yöntem sunmaktadır.
• NO₂ ile tüm nedenlere bağlı mortalite arasındaki pozitif ilişki O₃ için yapılan bu düzeltmeden sonra da devam etmektedir.
• O₃ için düzeltilmiş bağımsız tahmin bu analizde istatistiksel olarak anlamlı değildir.
• Daha fazla çalışmanın dönüştürülerek dahil edilmesi NO₂ tahmininin güven aralığını daraltmıştır.
• NO₂–O₃ atmosfer kimyasındaki mekânsal heterojenlik sağlık-etki tahminlerinde dikkate alınması gereken bir faktördür.

Çalışmanın desteklemediği sonuçlar

• NO₂'nin bu meta-analizle nedensel olarak ölüm oluşturduğu kesin biçimde kanıtlanmamıştır.
• OVB düzeltmesi bütün epidemiyolojik bias ve confounding kaynaklarını ortadan kaldırmaz.
• 2,643 milyar değerinin ölüm sayısı olduğu söylenemez.
• O₃'ün hiçbir sağlık etkisinin bulunmadığı sonucu çıkarılamaz.
• Tüm coğrafyalarda aynı NO₂–O₃ korelasyon yapısı varsayılamaz.
• Çok yüksek heterojenliğe rağmen tek bir pooled RR'nin her toplum için aynı bireysel riski temsil ettiği söylenemez.

Kaynak içi raporlama notları

Kaynak 80 çalışmanın dahil edildiğini bildirirken verilen Kuzey Amerika, Avrupa ve Asya sayıları toplam 78 etmektedir.

NO₂ naïve güven aralığı ana metin ile Table 1 arasında küçük farklılık göstermektedir.

Doğrudan iki-kirleticili NO₂ analizinin üst güven sınırı ana metinde 1,187, tabloda 1,186 olarak verilmektedir.

Bu farklılıkların nihai hakemli sürümde kontrol edilmesi gerekir.

Kaynak ve Yöntem Notu

Özgün çalışma:Long-term Exposure to Photochemical Oxidants and All-Cause Mortality: A Meta-Analysis of Two-Pollutant Models with Omitted Variable Bias Correction

Yazarlar: Xinshuai Song; Jianyu Deng; Ning Kang; Peng Zhou; Mingjin Tang; Tao Xue.

Bağlı kurumlar: Guangzhou Institute of Geochemistry, Chinese Academy of Sciences; Peking University Health Science Centre; National Institute of Health Data Science at Peking University; Peking University College of Environmental Sciences and Engineering; University of Chinese Academy of Sciences; Tianjin University.

Kaynak türü: Preprint.

Hakemlik durumu: Kaynak PDF açıkça “Preprint not peer reviewed” olarak işaretlenmiştir.

SSRN abstract ID: 6959529.

Çalışma tasarımı: Uzun dönem NO₂/O₃ maruziyeti ve tüm nedenlere bağlı mortalite için kohort çalışmalarının sistematik derlemesi, random-effects meta-analizi ve OVB düzeltmesi.

Dış atmosfer verisi: CAMS reanalysis; TOAR; CNEMC.

Finansman: National Natural Science Foundation of China — 42321003, 42422507, 42375179; National Key R&D Program of China — 2022YFC3703000, 2023YFC3708304; Beijing Natural Science Foundation — JQ24052.

Çıkar çatışması: Yazarlar çalışmayı etkileyebilecek bilinen finansal çıkar veya kişisel ilişki bildirmemiştir.

Temel yöntemsel sınırlar: OVB düzeltmesi, kirleticiler arasındaki korelasyon ve standart sapma oranlarının doğru biçimde eşleştirilmesine bağlıdır. Orijinal çalışmaların çoğunda co-pollutant confounding kontrolü yeterli değildir. Naïve meta-analizlerde heterojenlik I² > %99'dur. Dış atmosfer gridleri bireysel maruziyeti tam temsil etmeyebilir. Çalışma gözlemsel kohort sonuçlarını sentezlediği için nedensellik tek başına kanıtlanamaz.

Telif / yeniden kullanım notu: Yüklenen preprint üzerinde Creative Commons lisansı görünmemektedir. Kaynak forest plotları, küresel haritalar ve grafiksel abstract birebir yeniden kullanılmamalı; Verianla için bilimsel ilişkiler bağımsız görsel mimariyle yeniden çizilmelidir.

Verianla anlatım ilkesi: Bu içerik kaynak preprintin tam metin çevirisi değildir. OVB denklemi, meta-analiz sonuçları, korelasyon mekanizması, nüfus-kapsam hesabı ve yöntemsel sınırlılıklar korunmuş; çalışma bağımsız, öğretici ve nedensellik sınırlarına dikkat eden bir yapıda yeniden anlatılmıştır.