נספח: מקורות, נתונים וחישובים (השוואה ארצית, קריית טבעון מול שאר הארץ)

תאריך: 18.10.2026. נספח לדוח ההשוואה הארצית ולגרפים. המטרה: שכל מספר, טענה וחישוב בהשוואה יהיו ניתנים לבדיקה עצמאית. הקבצים עצמם בתיקייה חבילת ראיות - השוואה ארצית/ (מצורפת), והם רשומים בסעיף 10 עם SHA-256.

1. תוכן החבילה

2. מקורות וקישורים

הערה על קישורים: בקשות אוטומטיות לחלק מהאתרים הממשלתיים (gov.il, data.gov.il, האתר של עיריית קריית ים) נחסמות (403 או הפניה); בדפדפן הם נפתחים. לכן העותקים בחבילה הם העוגן: דוח 2024 באתר האיגוד זהה (SHA-256) לקובץ השמור, ולא נעשה שימוש בעותק של עיריית קריית ים.

3. הגדרות ונוסחאות

1. ממוצע שנתי: ממוצע כל השעות עם valid=True בשנה. כיסוי: שעות תקפות חלקי שעות השנה (8,760; 8,784 ב-2024).
2. תחנה-שנה נכללת אם הכיסוי השנתי לפחות 75% ובכל עונה (חורף דצמבר–פברואר, אביב מרץ–מאי, קיץ יוני–אוגוסט, סתיו ספטמבר–נובמבר) לפחות 50% שעות תקפות במזהם (למניעת הטיה עונתית). נדרש לפחות שנתיים תקפות לממוצע תלת-שנתי. תחנה-שנה עם PM2.5 מתחת ל-5 או מעל 60 (ערך שלילי באורות רבין, חדרה, 2024) הוצאה. PM10 מתחת ל-15 או מעל 150 הוצא (ערך מתחת ל-15 אינו סביר בישראל; נמצא אחד: אורות רבין 2024, 13.6). יציבות האוזון: תחנה שמספר ימי האוזון השנתי שלה (אחרי הכיסוי) משתנה פי יותר מ-2.5 בין שנותיה הוצאה מדירוג האוזון (שונות גבוהה; אינה בהכרח טעות, ותחנה 122 למשל תואמת לדוח האיגוד): תחנות 7, 11, 122, 472, 539.
3. יום תקף ל-PM2.5: לפחות 18 שעות תקפות. % ימים מעל 15: ימי ממוצע יממתי מעל 15 מק"ג/מ"ק חלקי הימים התקפים (15 = הנחיית WHO היממתית).
4. MDA8: לכל שעה, ממוצע ה-8 שעות המסתיימות בה (6 שעות תקפות לפחות מתוך 8); ה-MDA8 של יום הוא המרבי של ממוצעים אלה בשעות היום, ליום תקף נדרשים 18 ממוצעים לפחות. ימי אוזון: ימים עם MDA8 של 51 ppb לפחות, מנורמל ל-365 ימים (ימים × 365 / ימים תקפים). 51 ppb: 1 ppb אוזון = 48.00/24.4655 = 1.962 מק"ג/מ"ק ב-25°C ו-1 atm, ולכן 100/1.962 = 50.97, כלומר 100 מק"ג/מ"ק (הנחיית WHO ל-8 שעות) הם כ-51 ppb. ההשוואה ב-ppb מבטלת את תלות הדירוג בגורם ההמרה.
5. NO2: ממוצע ב-ppb × 46.0055/24.4655 = 1.8804 (25°C, 1 atm).
6. ימי אבק אזוריים: בכל שנה, חציון ה-PM10 היומי (יממה תקפה: 18 שעות) על פני כל תחנות ה-PM10 גבוה מ-60 מק"ג/מ"ק; נוסף היום שאחריו (53–61 ימים בשנה). PM2.5 בלי אבק: ממוצע הימים התקפים שאינם ימי אבק (נדרשים 150 ימים לפחות).
7. ממוצע תלת-שנתי לתחנה: ממוצע הערכים של השנים התקפות. חציון, רבעונים ודירוג: על פני תחנות. "X% מהתחנות נמוכות מטבעון": שיעור התחנות עם ערך נמוך ממש מערך טבעון. דירוג: מספר התחנות עם ערך גבוה ממש, ועוד אחת (1 = הגבוה ביותר).
8. סוג תחנה: תחבורה (סוג "תחבורתית" או "תחנות פנים מבניות", או בעלות נתיבי ישראל, נתיבי איילון, רכבת ישראל או התחנה המרכזית); תעשייה/מחצבה/תחנת כוח (סוג "על הגדר", "תעשייתית" או "נמלי ים", או בעלות/שם עם מחצבה, אבן וסיד, תחנת כוח, ת"כ, MRC, או.פי, נאות חובב, שברון, ים המלח, או שם התחנה כולל "תעשייה", או בעלות "נשר" או "שפיר"); רקע עירוני/כללי: השאר. מגבלה: הסיווג יוריסטי; אינו מבחין בין תחנות של חברת החשמל לתחנות כוח פרטיות, ותחנה אחת הסומנה כלא פעילה במטא-דאטה נשארת במדגם. נבדקה רגישות לסיווג ולהשמטת אזורים (סעיף 7).
9. אזור: לפי יישוב ובעלות וקואורדינטות (תל אביב וגוש דן, חיפה והקריות, ירושלים, נגב וערבה, דרום החוף והשפלה הדרומית, שרון וכרמל, עמקים וגליל, שפלה ומודיעין); טבעון בנפרד. קוד הסיווג ב-national_analysis.py (פונקציות klass ו-region).
10. מרחק מהחוף: (קו אורך התחנה − קו אורך החוף באותו קו רוחב) × 95 ק"מ, כשקו החוף הוא אינטרפולציה לינארית בין נקודות (קו רוחב 31.50→34.47, 31.66→34.55, 31.80→34.64, 32.08→34.76, 32.33→34.85, 32.45→34.88, 32.70→34.95, 32.82→35.00, 32.90→35.07, 33.09→35.10); הערכה גסה. גובה: מרשימת התחנות; לתחנות 384, 385 ו-454 שדה הגובה ב-API מכיל ערך שאינו גובה (קו אורך), ולתחנות 563 (0 מ') ו-586 (+232 מ' בבקעת הירדן) ערך לא סביר; כולן הוצאו. ארבע תחנות דרומית ל-31.1°N (דרום הנגב והערבה) נשארות בלי מרחק מחוף, כי קו החוף שבשימוש מכסה 31.5–33.1°N.
11. התאמה (סעיף 4.6 בדוח): ערך ~ תחבורה + תעשייה + גובה/100 + מרחק מחוף/10 (OLS, שגיאות תקן HC1), על ממוצעי תחנות, בלי תחנות מתחת ל-100 מ' מתחת לפני הים או רחוקות מ-80 ק"מ מהחוף; לאוזון בלי משתנה התחבורה (אין תחנות תחבורה עם אוזון).
12. ערכי ייחוס בגרפים: PM2.5: WHO 5 (שנתי), ישראל 25 (טבלה 22 בדוח 2025); NO2 ישראל 40 ואוזון ישראל 140 מק"ג/מ"ק ל-8 שעות (טבלאות האוזון והחנקן הדו-חמצני בדוח אשקלון 2024, עמ' 35 ו-37); PM10: WHO 15, ישראל 50 (ערך סביבה שנתי אחרי הורדת ימי סופת אבק, עמ' 103 בדוח 2024); NO2: WHO 10; אוזון: WHO 100 מק"ג/מ"ק ל-8 שעות, 3–4 ימים בשנה (טבלה 0.1 של WHO, הערה a).

4. תהליך (צעד אחר צעד)

1. fetch_national_year.py <שנה>: קורא את אסימון האורח מ-air.sviva.gov.il, ולכל תחנה פעילה ונייחת עם PM2.5 או אוזון (161) שולף חודש אחר חודש את הנתונים השעתיים מה-API ושומר ערוצי PM2.5, PM10, O3 ו-NO2.
2. national_analysis.py: מחשב מדדי תחנה-שנה, בקרת איכות, ימי אבק, סיווג, ממוצעי תחנה, דירוגים, רגרסיות; כותב national_summary_station_year.csv, national_summary_station.csv ופלט.
3. validate_national_vs_official.py: משווה ל-API את הערכים הרשמיים שהועתקו מדוחות האיגוד.
4. audit_recompute.py: חישוב עצמאי של כל המדדים מהנתונים הגולמיים, בקוד אחר (numpy), והשוואה לתוצאות הצינור.
5. classification_sensitivity.py: רגישות הדירוג לסיווג ולהשמטת אזורים.
6. build_national_report2.py ו-build_charts.py: מייצרים את הדוח והגרפים מהקבצים; כל מספר בדוח מחושב, והבנייה נכשלת אם טענה איכותית שנבדקת בקוד (למשל "טבעון נמוכה מכל האזורים באוזון") אינה נכונה.

הרצה מחדש: python -I fetch_national_year.py 2023 (וכן 2024 ו-2025), python -I national_analysis.py, python -I audit_recompute.py.

5. דוגמת חישוב מלאה: טבעון, 2025

• PM2.5: 8,656 שעות תקפות מתוך 8,760 (98.8%), סכום 126,242.8, ממוצע שנתי 14.58 מק"ג/מ"ק. ימים תקפים (18 שעות לפחות): 362; מתוכם 135 עם ממוצע יממתי מעל 15, כלומר 37.3%.
• NO2: 7,526 שעות תקפות, ממוצע 5.02 ppb × 1.8804 = 9.44 מק"ג/מ"ק.
• MDA8, יום לדוגמה 10.7.2025: ערכי האוזון השעתיים (ppb, שעות 00:00–23:00 לפי שעון חורף קבוע): [36.1, 32.8, 32.3, 30.3, 30.2, 27.7, 23.1, None, 41.2, 49.0, 55.6, 58.7, 55.6, 55.5, 53.8, 53.6, 51.7, 46.1, 39.6, 37.2, 34.0, 32.3, 35.6, 37.6]. ממוצעי 8 השעות המסתיימים בכל שעה: [40.69, 38.72, 36.7, 34.81, 33.29, 32.38, 30.78, 30.36, 31.09, 33.4, 36.73, 40.79, 44.41, 48.39, 52.77, 52.88, 54.19, 53.83, 51.83, 49.14, 46.44, 43.54, 41.26, 39.26]. ה-MDA8 של היום = המרבי = 54.1875 ppb.
• הערכים לכל שלוש השנים והממוצע התלת-שנתי: בטבלת תחנה-שנה (סעיף 9) ובקובץ national_summary_station_year.csv.

6. אימות מול הדוחות הרשמיים

הערה: שורות PM2.5 של 2024 מוצגות לשקיפות בלבד (הערכים כפי שפורסמו). עדכון 11.2026: האיגוד אישר שהטבלה מוסטת ושלח טבלה מתוקנת (טבעון 14.1); מול ה-API הפער הממוצע המוחלט 0.18 ב-11 תחנות (validate_2024_corrected.py, פלט validate_2024_corrected_output.txt). בטבלה 20 בדוח 2024 השמות מוסטים בשורה אחת ביחס לערכים (ראו הסיכום והפלט להלן), ולכן ההפרשים בשורות אלה גדולים ואינם מעידים על שגיאת API.

סיכום: PM2.5 2023 ו-2025 (22 השוואות, 11 תחנות בכל שנה, איגוד חיפה): תואם; PM10 2023–2024 (14 השוואות): הפער הממוצע המוחלט 0.05, המרבי 0.16. איגוד אשדוד-חבל יבנה (PM2.5, 8 תחנות, 2023–2025, 24 השוואות) ואיגוד אשקלון (PM2.5 2024, 10 תחנות): טבלאות בסעיף 6א. PM2.5 ב-2024 לא שימש לאימות: טבלה 20 בדוח 2024 מוסטת בשורה אחת. הוכחה (פלט validate_national_vs_official.py):

===== 2024 PM2.5 table (Table 20 of the 2024 report): test for a one-row shift between labels and values =====
API 2024 per station (annual mean, max daily): {'אחוזה': (12.86, 51.1), 'נווה שאנן': (12.02, 49.4), 'נשר': (13.51, 58.7), 'עצמאות': (16.53, 118.7), 'פארק הכרמל': (14.62, 57.4), 'קרית אתא': (16.1, 54.8), 'קרית בנימין': (13.92, 33.3), 'קרית חיים': (15.39, 35.0), 'קרית טבעון': (14.19, 40.0), 'קרית ים': (15.14, 38.6), 'קרית מוצקין': (15.68, 36.0), 'קרית ביאליק (לפי יישוב, תחנה 380)': (16.45, 45.3), 'הדר (תחבורתית; לפי שם, תחנה 391)': (14.25, 47.2)}
row  1  label=אחוזה        published mean/max 12.9/52.0 -> API station(s) matching both: [('אחוזה', 12.86, 51.1)]
row  2  label=איגוד        published mean/max 12.9/52.0 -> API station(s) matching both: [('אחוזה', 12.86, 51.1)]
row  3  label=הדר          published mean/max 16.4/40.3 -> API station(s) matching both: []
row  4  label=נווה שאנן    published mean/max 14.3/46.3 -> API station(s) matching both: [('הדר (תחבורתית; לפי שם, תחנה 391)', 14.25, 47.2)]
row  5  label=נשר          published mean/max 12.0/48.0 -> API station(s) matching both: [('נווה שאנן', 12.02, 49.4)]
row  6  label=עצמאות       published mean/max 13.5/58.6 -> API station(s) matching both: [('נשר', 13.51, 58.7)]
row  7  label=פארק הכרמל   published mean/max 16.2/116.0 -> API station(s) matching both: []
row  8  label=קרית אתא     published mean/max 14.5/46.0 -> API station(s) matching both: []
row  9  label=קרית ביאליק  published mean/max 16.1/54.1 -> API station(s) matching both: [('קרית אתא', 16.1, 54.8)]
row 10  label=קרית בנימין  published mean/max 16.5/44.2 -> API station(s) matching both: [('קרית ביאליק (לפי יישוב, תחנה 380)', 16.45, 45.3)]
row 11  label=קרית חיים    published mean/max 13.9/34.0 -> API station(s) matching both: [('קרית בנימין', 13.92, 33.3)]
row 12  label=קרית טבעון   published mean/max 15.4/40.3 -> API station(s) matching both: []
row 13  label=קרית ים      published mean/max 14.2/39.1 -> API station(s) matching both: [('קרית טבעון', 14.19, 40.0)]
row 14  label=קרית מוצקין  published mean/max 14.7/33.1 -> API station(s) matching both: []

Same-label match vs previous-row match (annual mean, ug/m3):
  נווה שאנן    API 12.02 | published at its own row 14.3 (diff -2.28) | published at the NEXT row 12.0 (diff 0.02)
  נשר          API 13.51 | published at its own row 12.0 (diff 1.51) | published at the NEXT row 13.5 (diff 0.01)
  עצמאות       API 16.53 | published at its own row 13.5 (diff 3.03) | published at the NEXT row 16.2 (diff 0.33)
  פארק הכרמל   API 14.62 | published at its own row 16.2 (diff -1.58) | published at the NEXT row 14.5 (diff 0.12)
  קרית אתא     API 16.10 | published at its own row 14.5 (diff 1.60) | published at the NEXT row 16.1 (diff 0.00)
  קרית בנימין  API 13.92 | published at its own row 16.5 (diff -2.58) | published at the NEXT row 13.9 (diff 0.02)
  קרית חיים    API 15.39 | published at its own row 13.9 (diff 1.49) | published at the NEXT row 15.4 (diff -0.01)
  קרית טבעון   API 14.19 | published at its own row 15.4 (diff -1.21) | published at the NEXT row 14.2 (diff -0.01)
  קרית ים      API 15.14 | published at its own row 14.2 (diff 0.94) | published at the NEXT row 14.7 (diff 0.44)
  קרית מוצקין  API 15.68 | published at its own row 14.7 (diff 0.98) | published at the NEXT row nan (diff nan)
mean abs diff, own row: 1.72 | next row: 0.11 (stations where the next row exists)

2025 report, Table 24 (comparison): labels list 12 stations, Table 22 lists 14 (Atzmaut and Carmel Park removed from labels but their values kept):
  Table 24 label אחוזה        value 12.5 | Table 22 value of the SAME station 12.5 | Table 22 value at the same ROW POSITION (אחוזה) 12.5
  Table 24 label איגוד        value 19.4 | Table 22 value of the SAME station 19.4 | Table 22 value at the same ROW POSITION (איגוד) 19.4
  Table 24 label הדר          value 14.6 | Table 22 value of the SAME station 14.6 | Table 22 value at the same ROW POSITION (הדר) 14.6
  Table 24 label נווה שאנן    value 11.8 | Table 22 value of the SAME station 11.8 | Table 22 value at the same ROW POSITION (נווה שאנן) 11.8
  Table 24 label נשר          value 13.5 | Table 22 value of the SAME station 13.5 | Table 22 value at the same ROW POSITION (נשר) 13.5
  Table 24 label קרית אתא     value 16.0 | Table 22 value of the SAME station 17.2 | Table 22 value at the same ROW POSITION (עצמאות) 16.0
  Table 24 label קרית ביאליק  value 13.0 | Table 22 value of the SAME station 16.6 | Table 22 value at the same ROW POSITION (פארק הכרמל) 13.0
  Table 24 label קרית בנימין  value 17.2 | Table 22 value of the SAME station 15.3 | Table 22 value at the same ROW POSITION (קרית אתא) 17.2
  Table 24 label קרית חיים    value 16.6 | Table 22 value of the SAME station 16.1 | Table 22 value at the same ROW POSITION (קרית ביאליק) 16.6
  Table 24 label קרית טבעון   value 15.3 | Table 22 value of the SAME station 14.5 | Table 22 value at the same ROW POSITION (קרית בנימין) 15.3
  Table 24 label קרית מוצקין  value 16.1 | Table 22 value of the SAME station 18.0 | Table 22 value at the same ROW POSITION (קרית חיים) 16.1
  Table 24 label קרית ים      value 14.5 | Table 22 value of the SAME station 15.7 | Table 22 value at the same ROW POSITION (קרית טבעון) 14.5

2023 report, Table 25 (comparison 2022 vs 2023, column 2023): labels list 12 stations, Table 23 lists 14 (Atzmaut and Carmel Park removed from the labels but their values kept):
  Table 25 label אחוזה        printed 14 | Table 23 value of the SAME station 13.9 (rounds to 14) | Table 23 value at the same ROW POSITION (אחוזה) 13.9 (rounds to 14)
  Table 25 label איגוד        printed 17 | Table 23 value of the SAME station 17.3 (rounds to 17) | Table 23 value at the same ROW POSITION (איגוד) 17.3 (rounds to 17)
  Table 25 label הדר          printed 14 | Table 23 value of the SAME station 14.1 (rounds to 14) | Table 23 value at the same ROW POSITION (הדר) 14.1 (rounds to 14)
  Table 25 label נווה שאנן    printed 13 | Table 23 value of the SAME station 13.0 (rounds to 13) | Table 23 value at the same ROW POSITION (נווה שאנן) 13.0 (rounds to 13)
  Table 25 label נשר          printed 14 | Table 23 value of the SAME station 13.5 (rounds to 14) | Table 23 value at the same ROW POSITION (נשר) 13.5 (rounds to 14)
  Table 25 label קרית אתא     printed 18 | Table 23 value of the SAME station 17.0 (rounds to 17) | Table 23 value at the same ROW POSITION (עצמאות) 18.0 (rounds to 18)
  Table 25 label קרית ביאליק  printed 15 | Table 23 value of the SAME station 16.7 (rounds to 17) | Table 23 value at the same ROW POSITION (פארק הכרמל) 15.0 (rounds to 15)
  Table 25 label קרית בנימין  printed 17 | Table 23 value of the SAME station 13.7 (rounds to 14) | Table 23 value at the same ROW POSITION (קרית אתא) 17.0 (rounds to 17)
  Table 25 label קרית חיים    printed 17 | Table 23 value of the SAME station 13.8 (rounds to 14) | Table 23 value at the same ROW POSITION (קרית ביאליק) 16.7 (rounds to 17)
  Table 25 label קרית טבעון   printed 14 | Table 23 value of the SAME station 15.4 (rounds to 15) | Table 23 value at the same ROW POSITION (קרית בנימין) 13.7 (rounds to 14)
  Table 25 label קרית מוצקין  printed 14 | Table 23 value of the SAME station 18.9 (rounds to 19) | Table 23 value at the same ROW POSITION (קרית חיים) 13.8 (rounds to 14)
  Table 25 label קרית ים      printed 15 | Table 23 value of the SAME station 16.2 (rounds to 16) | Table 23 value at the same ROW POSITION (קרית טבעון) 15.4 (rounds to 15)

בטבלה 24 בדוח 2025, כל הערכים אחרי "נשר" זהים לערך בטבלה 22 במיקום השורה (לא באותה תחנה), כי שתי תחנות הוסרו מרשימת התוויות והערכים שלהן נשארו.

6א. אימות מול דוחות איגודים נוספים (PM2.5)

הדוחות נמצאו בחיפוש ידני (קובץ "היכן חיפשנו ומה לא נמצא") ונבדקו: הקבצים בתיקייה 01. ערכי הדוח הועתקו מהטבלאות בעמודים המצוינים בסעיף 2, והושוו לממוצע שעות תקפות ב-API (validate_other_associations.py). הערה: בקובץ החיפוש עמוד 56 (2023), 41 (2024) ו-37 (2025) של אשדוד מסומנים כ-PM10, אך בפועל הם טבלאות בנזן; לכן PM10 באשדוד לא אומת.

איגוד אשדוד-חבל יבנה, PM2.5 2023–2025 (n=24, הפער הממוצע המוחלט 0.04, המרבי 0.13):

איגוד אשקלון, PM2.5 2024 (n=10, הפער הממוצע המוחלט 0.05, המרבי 0.12). הוצאה תחנה שבה ל-API פחות מ-50% מהשעות: שדרות (מזהה 126, כיסוי API 4%).

6ב. אימות NO2 ואוזון מול דוחות אשדוד ואשקלון

הערכים הועתקו מהטבלאות בדוחות (אשדוד: NO2 עמ' 53, 38, 34; אוזון עמ' 54, 39, 35; אשקלון 2024: אוזון טבלה 4 עמ' 35, NO2 טבלה 6 עמ' 37) והושוו ל-API (validate_no2_o3_other_associations.py, פלט validate_no2_o3_other_associations_output.txt). NO2: ממוצע API ב-ppb × 1.8804 מול הממוצע הרשמי (מק"ג/מ"ק). אוזון: גורם ההמרה אינו מפורסם, ולכן מחושב היחס הרשמי/API; יחס יציב בסביבות 1.96 (הגורם הפיזיקלי ב-25°C) פירושו התאמה.

תוצאות (תחנות עם כיסוי API של 85% לפחות): NO2 באשדוד 21 השוואות, פער ממוצע מוחלט 0.05 ומרבי 0.14; באשקלון 10 השוואות, 0.09 ו-0.17. אוזון: יחס הממוצע השנתי 1.952 (טווח 1.91–1.98, 17 השוואות); יחס המרבי השמונה-שעתי 1.81–1.96 (ה-API גבוה יותר מהרשמי בשיעור משתנה).

חריגים: תחנה 535 (אזור תעשייה צפוני, אשדוד 2025): ב-API כיסוי 63% וממוצע 25.7 מול 14.9 רשמי; תחנת שדרות (126): כיסוי API 4%; תחנה 365 (בני דרום) אינה בנתוני ה-API שנשלפו (הועברה לקבוצת יבנה מערב, 511, בפברואר 2025); תחנה 511 באוזון 2025: כיסוי API 67% וערך רשמי (31.2) שנראה כמספר ppb ולא מק"ג/מ"ק. הם הוצאו מהסטטיסטיקה ומוצגים בטבלאות.

הערה על יציבות האוזון: תחנת יד בנימין (122) הוצאה מדירוג האוזון הארצי בכלל היציבות (שונות פי יותר מ-2.5 בין שנים), אך נתוניה תואמים לדוחות האיגוד (ממוצע רשמי 52.9, 63.2, 69.5 ומרבי 132, 153, 182). הכלל מסמן שונות גבוהה, לא בהכרח טעות.

NO2:

אוזון:

7. בדיקות ביקורת עצמאיות

חישוב מחדש בקוד נפרד (audit_recompute.py):

audit: ozone-unstable stations excluded: [7, 11, 122, 472, 539] | implausible PM2.5 station-years set to NaN: 1
station-years recomputed: 415 | in pipeline file: 415
pm25         compared  262 | max abs diff 0.0000 | diffs above tol 0 | only audit has value 0 | only pipeline has value 0
pm10         compared  133 | max abs diff 0.0000 | diffs above tol 0 | only audit has value 0 | only pipeline has value 0
no2          compared  342 | max abs diff 0.0000 | diffs above tol 0 | only audit has value 0 | only pipeline has value 0
o3_days51    compared  194 | max abs diff 0.0000 | diffs above tol 0 | only audit has value 0 | only pipeline has value 0
o3_jja       compared  194 | max abs diff 0.0000 | diffs above tol 0 | only audit has value 0 | only pipeline has value 0
pm25_days15  compared  262 | max abs diff 0.0000 | diffs above tol 0 | only audit has value 0 | only pipeline has value 0

sanity: station-years with negative valid hourly values: {'NO2': 87, 'O3': 17, 'PM10': 102, 'PM2.5': 239} | total negative hours by channel: {'NO2': 8831, 'O3': 66, 'PM10': 2008, 'PM2.5': 19555} | values above 1000: 349
 id  year    ch  negatives  above1000
472  2024   NO2       1710          0
453  2024   NO2        961          0
551  2025   NO2        930          0
183  2025   NO2        511          0
160  2025   NO2        500          0
128  2025 PM2.5        433          0
 21  2024 PM2.5        425          0
472  2023   NO2        423          0
453  2025   NO2        362          0
499  2023   NO2        352          0

Tivon recomputed by year:
          pm25_cov   pm25   pm10   no2  o3_cov  o3_days51  o3_jja  pm25_days15
id  year                                                                      
137 2023     98.80  15.35  27.84  6.82   95.30     117.29   51.02        45.28
    2024     97.48  14.19  28.77  8.58   95.55     107.00   51.53        32.96
    2025     98.81  14.58  29.62  9.44   95.57     162.89   55.16        37.29

רגישות לערכים שליליים וקיצוניים (audit_sensitivity):

hours above 1000 by channel: {'NO2': 0, 'PM10': 346, 'PM2.5': 3} | channels with such values: {'PM10': 77, 'PM2.5': 2}
stations/years with max>1000 PM2.5 or NO2:
 id  year    ch  n_gt1000     max
306  2023 PM2.5         1 1174.45
607  2025 PM2.5         2 1083.69
(stations with at least two valid years, as in the report)
PM2.5 | base: Tivon 14.71 rank 67 of 88 median 16.17 | noneg: Tivon 14.83 rank 68 of 88 median 16.36 | zero: Tivon 14.72 rank 69 of 88 median 16.22 | nohi: Tivon 14.71 rank 67 of 88 median 16.17
NO2 | base: Tivon 8.28 rank 91 of 113 median 12.61 | noneg: Tivon 8.28 rank 91 of 113 median 12.71 | zero: Tivon 8.28 rank 91 of 113 median 12.64 | nohi: Tivon 8.28 rank 91 of 113 median 12.61
PM10 | base: Tivon 28.74 rank 41 of 45 median 40.85 | noneg: Tivon 28.76 rank 41 of 45 median 41.05 | zero: Tivon 28.74 rank 41 of 45 median 40.95 | nohi: Tivon 28.74 rank 41 of 45 median 39.69

רגישות לסיווג ולהשמטת אזורים (classification_sensitivity.py):

baseline (all stations): {'PM2.5': '23% lower, median 16.2, n=87', 'PM2.5 בלי אבק': '37% lower, median 14.6, n=87', 'אוזון ימים': '18% lower, median 176.7, n=65', 'NO2': '19% lower, median 12.6, n=113', 'PM10': '9% lower, median 40.9, n=45'}
רקע בלבד (בסיס)                                  {'PM2.5': '22%/15.8', 'PM2.5 בלי אבק': '40%/14.2', 'אוזון ימים': '15%/186.1', 'NO2': '23%/11.5', 'PM10': '11%/36.4'}
רקע בלי תחנות רשות שדות התעופה (371, 357)        {'PM2.5': '23%/15.7', 'PM2.5 בלי אבק': '42%/14.2', 'אוזון ימים': '15%/176.7', 'NO2': '24%/11.5', 'PM10': '11%/36.4'}
רקע בלי תחנות חברת חשמל                          {'PM2.5': '21%/15.7', 'PM2.5 בלי אבק': '40%/14.2', 'אוזון ימים': '12%/187.0', 'NO2': '21%/12.1', 'PM10': '12%/36.0'}
רקע בלי הנגב והערבה                              {'PM2.5': '22%/15.9', 'PM2.5 בלי אבק': '38%/14.2', 'אוזון ימים': '16%/176.7', 'NO2': '20%/11.7', 'PM10': '12%/36.0'}
כל התחנות, בלי תחנות תחבורה                      {'PM2.5': '26%/15.8', 'PM2.5 בלי אבק': '41%/14.2', 'אוזון ימים': '18%/176.7', 'NO2': '23%/11.6', 'PM10': '9%/40.9'}

jackknife over regions (drop one region at a time; % of remaining stations lower than Tivon):
  PM2.5          range 18% - 29% (full sample 23%)
  PM2.5 בלי אבק  range 35% - 43% (full sample 37%)
  אוזון ימים     range 15% - 20% (full sample 18%)
  NO2            range 16% - 24% (full sample 19%)
  PM10           range 3% - 12% (full sample 9%)

8. מגבלות ידועות

• נתוני API גולמיים. אומתו מול דוחות: איגוד חיפה (PM2.5 2023 ו-2025; PM10 2023–2024), אשדוד (PM2.5 2023–2025) ואשקלון (PM2.5 2024). לא אומתו: שרון-כרמל, ירושלים ותחנות המשרד להגנת הסביבה, חברת החשמל, הנגב והצפון; וכן NO2 ואוזון מול דוחות.
• שלוש שנים; 2025 שנת אוזון גבוהה בטבעון.
• סיווג התחנות והאזורים היוריסטי, ומרחק מהחוף מקורב.
• כלל ימי האבק פשוט, ואינו זהה לכלל של האיגוד.
• השוואה בין תחנות, לא חשיפת תושבים. תחנת טבעון על גג ברכס (201 מ').
• ערכים שליליים ב-PM2.5 ו-NO2 נשארו (נבדקה רגישות). תחנות עם נתונים עד 25% חסרים נכללות.
• אין תחנה בטבריה.

9. נתוני תחנות

9.1 מטא-דאטה, סיווג ואזור (כל התחנות ברשימה, גם אלה שלא עברו בקרת איכות)

9.2 תוצאות לכל תחנה-שנה (אחרי בקרת איכות; תא ריק = לא עבר את הבקרה או אין מזהם)

(הממוצע התלת-שנתי לכל תחנה, כפי שמופיע בדוח, בנספח ב של דוח ההשוואה וב-national_summary_station.csv.)

10. קבצים ו-SHA-256

תיקיית 06 - מסמכי הפרסום מכילה עותקים של המסמכים המפורסמים (לא נכללים ברשימה כי הם עשויים להתעדכן).