Παρακολούθηση της βλάστησης και αποκατάσταση μετά τη φωτιά: Μελέτη περίπτωσης στην κεντρική ενδοχώρα της Πορτογαλίας

[[Παρακολούθηση της βλάστησης και αποκατάσταση μετά τη φωτιά: Μελέτη περίπτωσης
στην κεντρική ενδοχώρα της Πορτογαλίας]]

'''Πρωτότυπος Τίτλος:''' « Vegetation Monitoring and Post-Fire Recovery: A Case Study
in the Centre Inland of Portugal »
'''Συγγραφείς:''' Cristina Alegria

'''Δημοσιεύθηκε:''' MDPI, Sustainability  6 October 2022

'''Σύνδεσμος :''' 
Sustainability | Free Full-Text | Vegetation Monitoring and Post-Fire Recovery: A Case Study in the Centre Inland of Portugal (mdpi.com)
'''Λέξεις-Κλειδιά:''' Θαλάσσια Πεύκη – Ευκάλυπτοι - Θαμνώνες- Απογραφή- φασματικοί δείκτες

[[Εικόνα: PAPER3_SALPAR_IMAGE1.jpg| thumb| right|'''Εικόνα 1:''' Χάρτης Περιοχής μελέτης - τοποθεσία Sarnadas de S. Simão (Πορτογαλία): (α) γεωγραφική θέση (6.04734° Δ, 6.11400° Α, 44.19665° Ν, 44.28486° Β) - (β) κάλυψη γης το 1995 (COS 1995) - (γ) κάλυψη γης το 2007 (COS 2007) - και (δ) κάλυψη γης το 2018 (COS 2018).]]
[[Εικόνα: PAPER3_SALPAR_IMAGE7.png| thumb| right|'''Εικόνα 2:''' Κλιματολογικά δεδομένα: (2020-2022) - μηνιαίες μέσες θερμοκρασίες (ελάχιστες και μέγιστες °C) και συνολική μηνιαία βροχόπτωση (mm).]]
[[Εικόνα: PAPER3_SALPAR_IMAGE10.png| thumb| right|'''Εικόνα 3:''' Ώριμη παραγωγή φυσικής αναγέννησης ναυτικής πεύκης (V και Wa) και μοντέλα NDVI_a: (α) V = f(NDVI_a)-διάγραμμα, μοντέλο και στατιστικά στοιχεία προσαρμογής και πρόβλεψης και (β) Wa = f(NDVI_a)-διάγραμμα, μοντέλο και στατιστικά στοιχεία προσαρμογής και πρόβλεψης.]]
[[Εικόνα: PAPER3_SALPAR_IMAGE11.png| thumb| right|'''Εικόνα 4:''' Ορθοφωτογραφίες (Ιούλιος 2007): τεμάχια (κυκλωμένα με πορτοκαλί χρώμα) με μέγιστο και ελάχιστο NDVI ανά τύπο κάλυψης.]]
[[Εικόνα: PAPER3_SALPAR_IMAGE12.png| thumb| right|'''Εικόνα 5:''' 5. Χρονική προοπτική των αγροτεμαχίων (κύκλοι με πορτοκαλί χρώμα) που κάηκαν το 2020: (α) ορθοφωτογραφίες Ιουλίου 2007, (β) εικόνες Sentinel-2 (28 Ιουλίου 2020) και (γ) εικόνες Sentinel-2 (9 Μαΐου 2022).]]

'''ΠΕΡΙΛΗΨΗ'''
Οι πυρκαγιές αποτελούν μείζον περιβαλλοντικό ζήτημα που έχει αντίκτυπο στην υποβάθμιση της γης. Οι φασματικοί δείκτες τηλεπισκόπησης παρέχουν πολύτιμες πληροφορίες για τον βραχυπρόθεσμο μετριασμό και την αποκατάσταση μετά από πυρκαγιές. Χρησιμοποιήθηκε μια περιοχή μελέτης στην κεντρική ενδοχώρα της Πορτογαλίας που καταλαμβάνεται από δάση Maritime pine και Eucalypts. Η περιοχή επλήγη από πυρκαγιές το 2003, το 2017 και το 2020. Οι στόχοι της μελέτης ήταν διττοί: (1) να υπολογιστεί ο δείκτης κανονικοποιημένης διαφοράς βλάστησης (NDVI) και με δεδομένα δασικής απογραφής να εξαχθεί ένα μοντέλο παραγωγής Maritime pine, να διαφοροποιηθούν τα αείφυλλα δάση κωνοφόρων (π.χ, Maritime pine), αείφυλλα δάση πλατύφυλλων (π.χ. Eucalypts) και θαμνώνες, και παρακολούθηση της βλάστησης και της αποκατάστασής της μετά την πυρκαγιά- και (2) υπολογισμός της διαφοράς του κανονικοποιημένου λόγου καύσης (NBR) μεταξύ των ημερομηνιών πριν και μετά την πυρκαγιά για την αξιολόγηση των επιπέδων σοβαρότητας της καύσης. Τα αγροτεμάχια μιας προηγούμενης δασικής απογραφής χρησιμοποιήθηκαν για την παρακολούθηση των τιμών NDVI το 2007 και από το 2020 έως το 2022. Χρησιμοποιήθηκε μια εναέρια κάλυψη το 2007 και οι εικόνες Sentinel-2 για τα έτη 2020 έως 2022. Τα γραμμικά μοντέλα προσάρμοσαν την παραγωγή θαλάσσιας πεύκης με τον μετασχηματισμένο NDVI ανά ηλικία, παρουσιάζοντας αποτελεσματικότητα προσαρμογής 60%. Η διαστρωμάτωση των τύπων κάλυψης ανάλογα με το στάδιο ανάπτυξης της συστάδας και την εμφάνιση πυρκαγιάς ήταν δυνατή με τη χρήση της χρονικής καμπύλης NDVI, η οποία έδειξε επίσης την επίδραση της πυρκαγιάς και των χαμηλών βροχοπτώσεων. Οι τύποι κάλυψης κατατάχθηκαν με βάση τις φθίνουσες τιμές του NDVI ως εξής: ώριμες φυτείες ευκαλύπτων, νεαρή αναγέννηση θαλασσινής πεύκης, ώριμη θαλασσινή πεύκη, νεαρές φυτείες ευκαλύπτων, θαμνότοπος με φράουλες, φυτείες ευκαλύπτων μετά την πυρκαγιά, θαλασσινή πεύκη μετά την πυρκαγιά, υψηλός θαμνότοπος και μικρός θαμνότοπος. Η ανάκαμψη της βλάστησης μετά την πυρκαγιά ήταν χαμηλότερη στις περιοχές με υψηλότερο επίπεδο σοβαρότητας της πυρκαγιάς. Οι περιοχές θαλάσσιας πεύκης έχασαν την ικανότητα φυσικής αναγέννησης λόγω των σύντομων κύκλων των πυρκαγιών. Οι φασματικοί δείκτες αποτέλεσαν αποτελεσματικά εργαλεία για τη διαφοροποίηση των τύπων κάλυψης και τη βοήθεια στην αξιολόγηση της κατάστασης των δασών και των θαμνώνων.

'''1. ΠΕΡΙΟΧΗ ΜΕΛΕΤΗΣ'''
Η περιοχή μελέτης είναι μια τοποθεσία στην κεντρική ενδοχώρα της Πορτογαλίας με 3100 εκτάρια έκτασης. Η περιοχή μελέτης, Sarnadas de S. Simão έχει υποστεί σημαντικές μεταβολές στη δασική της κάλυψη καθ' όλη τη διάρκεια των τελευταίων δεκαετιών. Σύμφωνα με τους επίσημους πορτογαλικούς χάρτες κάλυψης/χρήσης γης το 1995, η περιοχή μελέτης καταλαμβανόταν κυρίως από θαλάσσιο πευκοδάσος. Το 2007, λόγω της πυρκαγιάς του 2003, το δάσος θαλάσσιας πεύκης μειώθηκε σημαντικά. Μετά από μια δεκαετία, το 2017, μια άλλη σοβαρή πυρκαγιά ήταν η αιτία της συνεχιζόμενης μείωσης των δασών θαλάσσιας Πεύκης και της αύξηση των ευκαλύπτων. Χρησιμοποιήθηκαν δεδομένα δασικής απογραφής που διενεργήθηκε το 2007, όπου έγινε συστηματική δειγματοληψία, υποστηριζόμενη από ένα γεωαναφερμένο πλέγμα 500 μέτρων επέτρεψε την εγκατάσταση 60 κυκλικών τεμαχίων των 500 m2.  Λόγω της πυρκαγιάς του 2003 η βόρεια ζώνη καλύφθηκε είτε από θαμνώνες είτε από θάμνους είτε από νεαρές φυσικώς αναγεννημένες συστάδες θαλάσσιας πεύκης που ήταν πολύ μικρές για να μετρηθούν. Διαφορές στην κάλυψη της γης των τεμαχίων το 2007 και το πορτογαλικό επίσημο χάρτη κάλυψης γης COS 2007 οφείλονται στην ελάχιστη μονάδα χαρτογράφησης που χρησιμοποιήθηκε στο εν λόγω χάρτη σε σύγκριση με την έκταση του αγροτεμαχίου (500 m2).

'''2. ΔΕΔΟΜΕΝΑ ΚΑΙ ΜΕΘΟΔΟΛΟΓΙΑ'''

'''Εικόνες τηλεπισκόπησης (2007 και 2020-2022) - φασματικοί δείκτες'''
Για το έτος 2007 χρησιμοποιήθηκε  η επίσημη εναέρια κάλυψη για την Πορτογαλία που ελήφθη μεταξύ Ιουλίου και Νοεμβρίου του ίδιου έτους, με χωρική ανάλυση 0,5 m.. Για τα έτη από το 2020 έως το 2022, δηλαδή το έτος της τελευταίας πυρκαγιάς και μετά  χρησιμοποιήθηκαν Sentinel-2 εικόνες. Για την παρούσα μελέτη, χρησιμοποιήθηκαν οι εικόνες του πολυφασματικού οργάνου (MSI) Sentinel-2 για 10 ημερομηνίες των ετών 2020-2022. Ο NDVI υπολογίστηκε και από τις δύο σειρές εικόνων, αντίστοιχα, με χωρική ανάλυση 0,5 m για τις ορθοφωτογραφίες και με χωρική ανάλυση 10 m για τις εικόνες Sentinel-2. Το NBR υπολογίστηκε από τις εικόνες Sentinel-2, χρησιμοποιώντας τη ζώνη SWIR2 με χωρική ανάλυση 20 m, για τις ημερομηνίες πριν και μετά την πυρκαγιά στις 13 Σεπτεμβρίου 2020.

'''Κλιματολογικά δεδομένα-τοπικός σταθμός (2020-2022) '''
Οι μηνιαίες κλιματολογικές εκθέσεις δόθηκαν από τον επίσημο εθνική πύλη του Instituto Português do Mar e da Atmosfera (IPMA) από τον πλησιέστερο τοπικό κλιματολογικό σταθμό σε απόσταση 40 χλμ περίπου. Στην παρούσα μελέτη, τα κλιματολογικά δεδομένα που αφορούν τα υδρολογικά έτη 2019/2020, 2020/2021 και 2021/2022 (μέχρι τον Μάιο) χρησιμοποιήθηκαν για να προκύψει η εξέλιξη του μέσου μηνιαίου θερμοκρασιών (ελάχιστη και μέγιστη ◦C) και της συνολικής μηνιαίας βροχόπτωσης (mm). Η συσσωρευμένη μηνιαία βροχόπτωση (mm) για κάθε υδρολογικό έτος ήταν 708,9 mm, 695,8 mm και 326,2 mm (μέχρι τον Μάιο). Τα δεδομένα αυτά χρησιμοποιήθηκαν για την εξήγηση των τιμών του NDVI που παρατηρήθηκαν στην περιοχή μελέτης καθ' όλη τη διάρκεια από το 2020 έως το 2022.

'''Εκτίμηση του NDVI των δασικών απογραφικών επιφανειών'''
Οι εικόνες NDVI που υπολογίστηκαν από τις ορθοφωτογραφίες που ελήφθησαν τον Ιούλιο του 2007 επέτρεψε την αξιολόγηση του μέσου NDVI για κάθε αγροτεμάχιο δασικής απογραφής (500 m2 ) χρησιμοποιώντας ένα υπόδειγμα που αποτελείται από έναν πυρήνα 20 × 20 εικονοστοιχείων γύρω από το κέντρο του αγροτεμαχίου (100 m2) Οι εικόνες NDVI που υπολογίζονται από τις εικόνες Sentinel-2 για τις 10 ημερομηνίες.  Στη συνέχεια, η παραγωγή ξύλου Maritime pine (V) και η παραγωγή υπέργειας βιομάζας (Wa) προσαρμόστηκαν με απλή γραμμική παλινδρόμηση χρησιμοποιώντας το μέσο NDVI των αγροτεμαχίων ως επεξηγηματική μεταβλητή. Λόγω της χαμηλής συσχέτισης που παρατηρήθηκε μεταξύ του μέσου όρου των αγροτεμαχίων NDVI και των V και Wa, ο μετασχηματισμένος NDVI, που ορίστηκε ως NDVI_a = NDVI (t/R) με R = 40 έτη, δοκιμάστηκε επίσης. 

'''NDVI ανά τύπο κάλυψης (2007 και 2020-2021) '''
Ο μέσος NDVI των αγροτεμαχίων το 2007 χρησιμοποιήθηκε επίσης για τη διερεύνηση των διαφορών τύπων κάλυψης (π.χ. θάμνος, αναγέννηση θαλάσσιας πεύκης, θαλάσσιας πεύκης και αμιγούς και μικτού ευκαλύπτου). Στη συνέχεια, ο τύπος κάλυψης σε κάθε αγροτεμάχιο το 2020 επικαιροποιήθηκε χρησιμοποιώντας το πορτογαλικό επίσημο χάρτη κάλυψης γης COS 2018 και από φωτογραφία εικόνας του Google Earth ερμηνεία (Μάιος-Ιούνιος 2019). Και πάλι, ο NDVI των αγροτεμαχίων χρησιμοποιήθηκε για τη διερεύνηση των διαφορών μεταξύ των τύπων κάλυψης (π.χ. θάμνος- θαλάσσια πεύκη- ευκάλυπτος) καθ' όλη τη διάρκεια της καλλιεργητικής περιόδου 2020 περιόδου (10 Μαρτίου 2020, 29 Μαΐου 2020, 28 Ιουνίου 2020, 28 Ιουλίου 2020 και 26 Σεπτεμβρίου 2020).
Πραγματοποιήθηκε έλεγχος Z-test δύο δειγμάτων για τους μέσους όρους, σε επίπεδο εμπιστοσύνης 95%, για να σύγκριση του μέσου NDVI των αγροτεμαχίων ανά τύπο κάλυψης τον Ιούλιο του 2007 και στις 28 Ιουλίου 2020. Στη συνέχεια, η NDVI ann-min και η μέθοδος CV χρησιμοποιήθηκε για τη διαφοροποίηση των τύπων κάλυψης καθ' όλη τη διάρκεια του 2020 καλλιεργητική περίοδο με βάση τη χρονική καμπύλη NDVI. Στη συνέχεια, πραγματοποιήθηκε περαιτέρω λεπτομερής ανάλυση για τα αγροτεμάχια συγκεκριμένων τύπων κάλυψης καθ' όλη τη διάρκεια της καλλιεργητικής περιόδου 2020 με βάση την χρονική καμπύλη NDVI.

'''Παρακολούθηση της σοβαρότητας της καύσης και της αποκατάστασης της βλάστησης μετά τη φωτιά (2020-2022)'''
Το έτος 2020, μια πυρκαγιά στις 13 Σεπτεμβρίου επηρέασε την περιοχή μελέτης. Τα σημάδια εγκαυμάτων της πυρκαγιάς διερευνήθηκαν με τη χρήση του NBR πριν από την πυρκαγιά (6 Σεπτεμβρίου 2020) και μετά την πυρκαγιά (26 Σεπτεμβρίου 2020). Στη συνέχεια, προέκυψε η διαφορά στο NBR μεταξύ των δύο ημερομηνιών. Στη συνέχεια, το dNBR χρησιμοποιήθηκε για την ταξινόμηση της σοβαρότητας της καύσης σύμφωνα με τα επίπεδα που προτείνονται από την Ευρωπαϊκή Υπηρεσία Πληροφοριών για τις Δασικές Πυρκαγιές (EFFIS). Ο πορτογαλικός επίσημος χάρτης καμένων εκτάσεων του 2020 επικαλύφθηκε στην εικόνα dNBR και χρησιμοποιήθηκε ως αναφοράς. Ο NDVI χρησιμοποιήθηκε για την παρακολούθηση της ανάκαμψης της βλάστησης μετά την πυρκαγιά (26 Σεπτεμβρίου 2020) κατά τη διάρκεια τριών εξαμηνιαίων περιόδων από το 2021 έως το 2022 (15 Μαρτίου 2021, 21 Σεπτεμβρίου 2021 και 9 Μαΐου 2022). Όπως είναι γνωστό, οι τιμές του NDVI κυμαίνονται από -1 έως 1, όπου το νερό οι επιφάνειες νερού, οι τεχνητές κατασκευές, οι βράχοι, τα σύννεφα και το χιόνι αντιστοιχούν σε αρνητικές τιμές, το γυμνό έδαφος συνήθως εμπίπτει στο εύρος 0,1-0,2 και τα φυτά θα έχουν πάντα θετικές τιμές μεταξύ 0,2 και 1. Ειδικότερα, ένας υγιής, πυκνός θόλος βλάστησης θα πρέπει να έχει τιμές πάνω από 0,5, και η αραιή βλάστηση πιθανότατα θα εμπίπτει στο εύρος 0,2 έως 0,5. Συνεπώς, ο NDVI επαναταξινομήθηκε λαμβάνοντας υπόψη τρεις κλάσεις πρασίνου για την παρακολούθηση της βλάστησης ανάκαμψη μετά την πυρκαγιά Επιπλέον, τα καμένα αγροτεμάχια στην πυρκαγιά του 2020 παρακολουθήθηκαν για να κατανοηθεί η NDVI όσον αφορά τα επίπεδα σοβαρότητας της καύσης μεταξύ 2020 και 2022 (10 Μαρτίου 2020, 29 Μαΐου 2020, 28 Ιουνίου 2020, 28 Ιουλίου 2020, 6 Σεπτεμβρίου 2020, 26 Σεπτεμβρίου 2020, 15 Μαρτίου 2021, 14 Μαΐου 20021, 21 Σεπτεμβρίου 2021, 30 Μαρτίου 2022 και 9 Μαΐου 2022).

'''3. ΑΠΟΤΕΛΕΣΜΑΤΑ'''

'''NDVI και παραγωγή θαλάσσιας πεύκης (2007) '''
Η συσχέτιση μεταξύ του μέσου NDVI των αγροτεμαχίων θαλάσσιας πεύκης και των μεταβλητών της συστάδας ήταν φτωχή (μικρότερη από 0,22). Η καλύτερη συσχέτιση πέντε μεταβλητών ήταν η ακόλουθη: G (-0.22); Sh25 (-0,18), N (0,17), Wl (0,17) και Wr (0,17). Οι μεταβλητές G και N είναι απόλυτες μετρήσεις της πυκνότητας της συστάδας, η μεταβλητή Sh25 είναι δείκτης της παραγωγικότητας της περιοχής και οι μεταβλητές Wl και Wr δίνουν μια εικόνα, αντίστοιχα, της ικανότητας φωτοσύνθεσης της κόμης των δέντρων και των ριζών των δέντρων του χώρου ανάπτυξης στη συστάδα. Όλες αυτές οι μεταβλητές είναι σημαντικές παράμετροι που καθορίζουν την παραγωγή της συστάδας. Όσον αφορά τη μέση συσχέτιση του NDVI των αγροτεμαχίων με την παραγωγή ξύλου (V) και την παραγωγή υπέργειας βιομάζας (Wa), διαπιστώθηκε ότι είναι πολύ χαμηλή (περίπου 0,03). Μόνο, χρησιμοποιώντας του μετασχηματισμένου NDVI_a, επιτεύχθηκαν συσχετίσεις μεταξύ 0,72 και 0,73. Επομένως, ως αναμενόταν, η αποτελεσματικότητα μοντελοποίησης των γραμμικών μοντέλων V = f(NDVI_a) και Wa = f(NDVI_a) δεν ήταν πολύ υψηλή (R2 γύρω στο 60%) (Σχήμα 8α,β). Αυτές οι συστάδες θαλάσσιας πεύκης δημιουργήθηκαν με φυσική αναγέννηση και, το 2007, ήταν κυρίως υποαποθεματοποιημένες ή πολύ υποαποθεματοποιημένες. Η ανύπαρκτη ή ακατάλληλη δασοκομική παρέμβαση μπορεί να εξηγήσει αυτή την υποπληθυσμιακή κατάσταση, όπως η εφαρμογή μία και μοναδική και καθυστερημένη αραίωση υψηλού βαθμού χωρίς να λαμβάνεται υπόψη η σταθερότητα της συστάδας. Τόσο το καθεστώς εγκατάστασης αυτών των συστάδων (φυσική αναγέννηση, συνεπώς με ακανόνιστες αποστάσεις) όσο και η απουσία κατάλληλης  διαχείρισης της συστάδας μεταφέρουν υψηλή μεταβλητότητα των συνθηκών ανάπτυξης μεταξύ των συστάδων στις οριζόντιες και κάθετες δομές των συστάδων (κατανομές διαμέτρου και ύψους). Η μεταβλητότητα αυτή έχει αρνητικό αντίκτυπο στις προσπάθειες μοντελοποίησης, όπως επαληθεύτηκε. Αυτές οι συστάδες έχουν επίσης σημαντική κάλυψη του υπορόφου (κατά μέσο όρο περίπου 54%) που μπορεί να έχει επηρεάσει τις παρατηρούμενες τιμές NDVI. Ο NDVI αυτών των εννέα αγροτεμαχίων αναλύθηκε όσον αφορά τη χρονική καμπύλη NDVI καθ' όλη τη διάρκεια την καλλιεργητική περίοδο 2020 (Μάρτιος-Σεπτέμβριος 2020) και το NDVI ann-min και CV αξιολογήθηκαν. Οι τύποι κάλυψης των εννέα αγροτεμαχίων μπορούσαν να καταταγούν με βάση τις φθίνουσες τιμές NDVI ως εξής:

'''Σοβαρότητα και παρακολούθηση της αποκατάστασης της βλάστησης μετά την πυρκαγιά'''
Το σημάδι καύσης της πυρκαγιάς του 2020 με τη χρήση του dNBR μεταξύ της ημερομηνίας πριν από την πυρκαγιά και  μετά. Αυτό επέτρεψε την παρατήρηση μιας πολύ υψηλής σοβαρότητας καύσης σε όλες τις περιοχές μελέτης που κάηκαν είναι η ανάκαμψη της δασικής βλάστησης μετά την πυρκαγιά, η οποία έδειξε, ότι στις αρχές της άνοιξης η ανάκαμψη δεν ήταν πολύ σημαντική. Ωστόσο, σημαντική ανάκαμψη παρατηρήθηκε ήδη στα τέλη του καλοκαιριού. Το επόμενο έτος παρατηρήθηκε συνολική μείωση του πρασίνου λόγω ενός πολύ ξηρού χειμώνα. Μόνο τέσσερα αγροτεμάχια (αγροτεμάχια n◦ 26, 35, 40 και 45) κάηκαν στην πυρκαγιά του 2020, αλλά με διαφορετικές σφοδρότητες πυρκαγιάς. Αυτά τα αγροτεμάχια κάηκαν επίσης στις πυρκαγιές του 2003 και του 2017. Το 2007, η κάλυψη γης του αγροτεμαχίου n◦ 26 ήταν θάμνος (50% Erica sp.), το αγροτεμάχιο n◦ 35 ήταν θάμνος (50% Cytisus sp.), το αγροτεμάχιο n◦ 40 ήταν θάμνος (60% Arbutus unedo) και το αγροτεμάχιο n◦ 45 ήταν θάμνος (50% Cistus sp.) (Παράρτημα Α-Σχήμα Α5). Και πάλι, στις αρχές της άνοιξης (15 Μαρτίου 2021), η αποκατάσταση της εδαφοκάλυψης αυτών των τεμαχίων δεν ήταν πολύ σημαντική, ενώ ακολούθησε στα τέλη του καλοκαιριού (21 Σεπτεμβρίου 2021) από μια σημαντική ανάκαμψη, και το επόμενο έτος, νωρίς την άνοιξη (9 Μαΐου 2022) από μια συνολική μείωση του πρασίνου λόγω μιας πολύ ξηρής περιόδου. χειμώνα (Σχήμα 4 - Ιανουάριος και Φεβρουάριος 2022). Αυτά τα αγροτεμάχια κάηκαν επίσης στις πυρκαγιές του 2003 και του 2017. Το 2007, η κάλυψη γης στο αγροτεμάχιο 26 ήταν θάμνος (50% Erica sp.), στο αγροτεμάχιο 26  35 ήταν θάμνος (50% Cytisus sp.), το αγροτεμάχιο 40 ήταν θάμνος (60% Arbutus unedo) και το αγροτεμάχιο 45 ήταν θάμνος (50% Cistus sp.).  Και πάλι, νωρίς την άνοιξη (15 Μαρτίου ) Εικόνα 16. Καμένα τεμάχια της δασικής απογραφής (σημειωμένα με μπλε χρώμα) το 2020-μετά την αποκατάσταση από την πυρκαγιά: σοβαρότητα της καύσης επίπεδα, κλάσεις NDVI 15 Μαρτίου 2021, κλάσεις NDVI 21 Σεπτεμβρίου 2021 και κλάσεις NDVI 9 Μαΐου 2022.

'''4. ΣΥΜΠΕΡΑΣΜΑΤΑ'''
Τα κύρια συμπεράσματα της μελέτης αυτής είναι τα εξής: α) τα γραμμικά μοντέλα της θαλάσσιας πεύκης (ξύλο και βιομάζα) με τον NDVI_a παρουσίασαν αποτελεσματικότητα προσαρμογής περίπου 60%, β) οι τιμές NDVI που υπολογίστηκαν από τις ορθοφωτογραφικές εικόνες βρέθηκαν να μην είναι συγκρίσιμες με εκείνες που υπολογίζονται από τις 
εικόνες Sentinel-2 λόγω των διαφορών στους αισθητήρες, εύρους ζώνης και χωρικής ανάλυσης, 3) οι τιμές NDVI των αγροτεμαχίων ανά τύπο κάλυψης, τόσο το 2007 όσο και το 2020, παρουσίασαν μεγάλη μεταβλητότητα και δεν διαπιστώθηκαν στατιστικά σημαντικές διαφορές μεταξύ της ώριμης θαλάσσιας πεύκης, των ευκαλύπτων και των θαμνώνων. 4) Η διαφοροποίηση μεταξύ των τύπων κάλυψης ανά στάδιο ανάπτυξης της συστάδας και ανά εμφάνιση πυρκαγιάς ήταν δυνατή όταν η μέθοδος NDVI-CV εφαρμόστηκε στη χρονική καμπύλη NDVI καθ' όλη τη διάρκεια της καλλιεργητικής περιόδου του 2020. 5) Ο NDVI και το NBR έδειξαν ότι η ανάκαμψη της βλάστησης μετά την πυρκαγιά ήταν χαμηλότερη σε υψηλότερης σοβαρότητας εγκαύματα. 6) Λόγω της αλληλουχίας μικρών κύκλων πυρκαγιών στην περιοχή μελέτης, η θαλάσσια πεύκη έχουν χάσει την ικανότητα φυσικής αναγέννησής τους. 7) Οι θαμνώδεις εκτάσεις (π.χ. Arbutus unedo) είχαν πολύ σημαντική ανάκαμψη μετά την πυρκαγιά ένα χρόνο αργότερα και η χρονική καμπύλη NDVI των καμένων επιφανειών κατά την περίοδο 2000-2022 έδειξε σαφώς τη μείωση του NDVI λόγω των επιπτώσεων της πυρκαγιάς (13 Σεπτεμβρίου 2020) και των χαμηλών βροχοπτώσεων του περασμένου χειμώνα (Ιανουάριος και Φεβρουάριος 2022). 8) Συνοπτικά, ο NDVI και ο NBR αποδείχθηκαν αποτελεσματικοί στη διαφοροποίηση των τύπων κάλυψης και ως εργαλεία παρακολούθησης που βοηθούν στην αξιολόγηση της κατάστασης των δασών και των θαμνώνων.