Θερμικά σοκ και ανθεκτικότητα του περιφερειακού συστήματος υγείας στην Ελλάδα: χωροχρονική μοντελοποίηση θερμοαποδιδόμενου φορτίου θνησιμότητας, υγειονομικής δυναμικότητας και θεσμικής ποιότητας
Loading...
Date
Authors
Journal Title
Journal ISSN
Volume Title
Publisher
Μεταπτυχιακό Πρόγραμμα Δημόσιας Διοίκησης, Σχολή Οικονομικών Επιστημών και Διοίκησης, Πανεπιστήμιο Νεάπολις Πάφος
Abstract
Κύριος στόχος της Διπλωματικής Εργασίας
Η παρούσα Διπλωματική Εργασία διερευνά τη σχέση ανάμεσα στην εβδομαδιαία θερμική ένταση (weekly heat intensity) και τη θνησιμότητα από όλες τις αιτίες (all-cause mortality) στις δεκατρείς ελληνικές περιφέρειες της ονοματολογίας εδαφικών μονάδων στατιστικής επιπέδου 2 (Nomenclature of Territorial Units for Statistics level 2, NUTS-2) για την περίοδο 2015-2024. Στόχος είναι η εκτίμηση της θερμικής κλίσης θνησιμότητας (heat-mortality gradient) σε εβδομαδιαία περιφερειακή κλίμακα και η μετατροπή του ευρήματος σε πλαίσιο περιφερειακής προτεραιοποίησης για τη δημόσια υγεία.
Μεθοδολογία
Η μελέτη συγκροτεί διαχρονικό-διαπεριφερειακό πλαίσιο δεδομένων θερμής περιόδου (warm-season weekly regional panel), με μονάδα ανάλυσης την περιφέρεια-εβδομάδα (region-week). Τα δεδομένα προέρχονται από τη Eurostat, την επανάλυση χερσαίων κλιματικών δεδομένων ERA5-Land (European ReAnalysis version 5 Land) και τον Ευρωπαϊκό Δείκτη Ποιότητας Διακυβέρνησης (European Quality of Government Index, EQI). Το κύριο πλήρες δείγμα (primary complete-case sample), μετά τον αποκλεισμό της πανδημικής διατάραξης 2020-2021, περιλαμβάνει 2.275 περιφερειακές εβδομάδες. Η βασική έκθεση είναι ο δείκτης εβδομαδιαίας θερμικής έντασης (heatintensity_week), δηλαδή το άθροισμα των ημερήσιων βαθμοημερών υπέρβασης (degree-days of exceedance) πάνω από το περιφερειακό 95ο εκατοστημόριο της ημερήσιας μέγιστης θερμοκρασίας (q95 Tmax). Το κύριο υπόδειγμα είναι αρνητικό διωνυμικό μοντέλο απαριθμήσεων (negative-binomial count model) με λογαριθμική σύνδεση (log link), πληθυσμιακό αντιστάθμισμα (population offset) και σταθερές επιδράσεις (fixed effects) περιφέρειας, έτους και εβδομάδας θερμής περιόδου. Η συμπερασματολογία ελέγχεται με διόρθωση CR2/Satterthwaite, κατανεμημένα μη γραμμικά μοντέλα υστέρησης (Distributed Lag Non-linear Models, DLNM), ηλικιακή στρωμάτωση, εικονικούς ελέγχους (placebo checks), σταθερές επιδράσεις έτους-εβδομάδας (year-week fixed effects) και αναδρομικό επιχειρησιακό έλεγχο του κανόνα P90.
Ευρήματα
Τα αποτελέσματα δείχνουν θετική σχέση θερμικής έντασης και θνησιμότητας. Η κύρια εκτίμηση αποδίδει λόγο επίπτωσης ρυθμού (Incidence Rate Ratio, IRR) 1,0184 ανά πρόσθετη °C-ημέρα υπέρβασης, με 95% διάστημα εμπιστοσύνης CR2 1,0145-1,0223 και p=0,00003. Η ομάδα ηλικίας 80 ετών και άνω εμφανίζεται ως το σταθερότερο στρώμα ευαλωτότητας, με IRR 1,0189. Η περιφερειακή ανάλυση διακρίνει την απόλυτη περιφερειακή αντίθεση προσαρμογής (regional adaptation contrast, AC), τον δείκτη περιφερειακής αντίθεσης προσαρμογής συνολικού ηλικιακού φάσματος κανονικοποιημένο ως προς τον πληθυσμό 80+ (older-population-normalised all-age regional adaptation contrast rate, ACRate80), τον επιχειρησιακό κίνδυνο ψευδώς αρνητικών εβδομάδων (false-negative operational risk) και το διαρθρωτικό πλαίσιο (structural context).
Συμπεράσματα
Η θερμότητα αναδεικνύεται ως πρόβλημα δημόσιας υγείας, περιφερειακής ετοιμότητας και ανθεκτικότητας συστημάτων υγείας. Τα ευρήματα υποστηρίζουν στοχευμένη προστασία των ατόμων 80+, εβδομαδιαία παρακολούθηση P90 και διαφοροποιημένο σχεδιασμό σχεδίων δράσης θερμότητας και υγείας (Heat-Health Action Plans, HHAP). Οι δείκτες κλινών, ιατρών και θεσμικής ποιότητας ερμηνεύονται ως οριακές θερμικές κλίσεις περιορισμένες στο εμπειρικό στήριγμα (support-restricted marginal heat gradients) και όχι ως καθολικά κατώφλια επάρκειας (universal adequacy thresholds).
Description
Main Objective of the Master’s Thesis
This Master’s Thesis investigates the relationship between weekly heat intensity and all-cause mortality across the thirteen Greek regions of the Nomenclature of Territorial Units for Statistics level 2 (NUTS-2) during the period 2015–2024. The objective is to estimate the heat-mortality gradient at a weekly regional scale and to translate the empirical finding into a regional prioritisation framework for public health.
Methodology
The study constructs a warm-season weekly regional panel, with the region-week as the unit of analysis. The data are drawn from Eurostat, the ERA5-Land reanalysis dataset (European ReAnalysis version 5 Land), and the European Quality of Government Index (EQI). The primary complete-case sample, after excluding the pandemic disruption of 2020–2021, consists of 2,275 region-weeks.
The main exposure variable is the weekly heat intensity indicator (heatintensity_week), defined as the sum of daily degree-days of exceedance above the regional 95th percentile of daily maximum temperature (q95 Tmax). The main specification is a negative-binomial count model with a log link, a population offset, and fixed effects for region, year, and warm-season week. Inference is assessed using CR2/Satterthwaite correction, Distributed Lag Non-linear Models (DLNM), age stratification, placebo checks, year-week fixed effects, and a retrospective operational audit of the P90 rule.
Findings
The results indicate a positive association between heat intensity and mortality. The main estimate yields an Incidence Rate Ratio (IRR) of 1.0184 per additional °C-day of exceedance, with a 95% CR2 confidence interval of 1.0145–1.0223 and p=0.00003. The population aged 80 years and above emerges as the most stable vulnerability stratum, with an IRR of 1.0189.
The regional analysis distinguishes between absolute regional adaptation contrast (AC), older-population-normalised all-age regional adaptation contrast rate (ACRate80), false-negative operational risk, and structural context. This distinction indicates that regional prioritisation should not rely on a single indicator, but on multiple complementary lenses of public-health risk, vulnerability, and operational readiness.
Conclusions
Heat emerges as a public-health, regional preparedness, and health-system resilience issue. The findings support targeted protection of people aged 80 and above, weekly P90 monitoring, and differentiated planning of Heat-Health Action Plans (HHAP) at NUTS-2 level. Indicators of hospital beds, physicians, and governance quality are interpreted as support-restricted marginal heat gradients and not as universal adequacy thresholds. Therefore, the study does not convert observational associations into absolute causal claims, but provides a technically grounded and reproducible framework for epidemiological assessment, operational surveillance, and regional accountability.
Keywords
θερμικά σοκ, θνησιμότητα από όλες τις αιτίες, περιφερειακή ανθεκτικότητα συστημάτων υγείας, θερμοαποδιδόμενο φορτίο, αρνητικό διωνυμικό μοντέλο (negative-binomial model), θεσμική ποιότητα, heat shocks, all-cause mortality, regional health-system resilience, heat-attributable burden, negative-binomial model, governance quality.