• Mercury in European topsoils: Anthropogenic sources, stocks and fluxes

    Map of Hg stock (g ha−1) in European topsoils.

    Panos Panagosa, Martin Jiskrab, Pasquale Borrellic, Leonidas Liakosa, Cristiano Ballabioa. 2021. “Mercury in European Topsoils: Anthropogenic Sources, Stocks and Fluxes.” Environmental Research, June, 111556. https://doi.org/10.1016/j.envres.2021.111556.

    a. European Commission, Joint Research Centre (JRC), Ispra, Italy
    b. Environmental Geosciences, University of Basel, Switzerland
    c. Department of Earth and Environmental Sciences, University of Pavia, 27100, Pavia, Italy

    Mercury (Hg) is one of the most dangerous pollutants worldwide. In the European Union (EU), we recently estimated the Hg distribution in topsoil using 21,591 samples and a series of geo-physical inputs. In this manuscript, we investigate the impact of mining activities, chrol-alkali industries and other diffuse pollution sources as primary anthropogenic sources of Hg hotspots in the EU. Based on Hg measured soil samples, we modelled the Hg pool in EU topsoils, which totals about 44.8 Gg, with an average density of 103 g ha−1. As a following step, we coupled the estimated Hg stocks in topsoil with the pan-European assessment of soil loss due to water erosion and sediment distribution. In the European Union and UK, we estimated that about 43 Mg Hg yr−1 are displaced by water erosion and c. a. 6 Mg Hg yr−1 are transferred with sediments to river basins and eventually released to coastal Oceans. The Mediterranean Sea receives almost half (2.94 Mg yr−1) of the Hg fluxes to coastal oceans and it records the highest quantity of Hg sediments. This is the result of elevated soil Hg concentration and high erosion rates in the catchments draining into the Mediterranean Sea. This work contributes to new knowledge in support of the policy development in the EU on the Zero Pollution Action Plan and the Sustainable Development Goal (SDGs) 3.9 and 14.1, which both have as an objective to reduce soil pollution by 2030.

    Link to Dataset: Mercury content in the European Union topsoil

    [Διαβάστε περισσότερα...]
  • An in-depth statistical analysis of the rainstorms erosivity in Europe

    Gini coefficients (G) for the 1181 stations across Europe using all rainfall erosive events included in the REDES database. As a background map, the six K-means clusters as defined by Ballabio et al. (2017) are shown.

    Nejc Bezaka, Matjaž Mikoša, Pasquale Borrellibc, Leonidas Liakosd, Panos Panagosd. 2021. “An In-Depth Statistical Analysis of the Rainstorms Erosivity in Europe.” CATENA 206 (November): 105577. https://doi.org/10.1016/j.catena.2021.105577.

    a .University of Ljubljana, Faculty of Civil and Geodetic Engineering, Slovenia
    b. Department of Earth and Environmental Sciences, University of Pavia, Italy
    c. Department of Biological Environment, Kangwon National University, Chuncheon 24341, Republic of Korea
    d. European Commission, Joint Research Centre (JRC), Ispra, Italy

    Heavy rainstorms play a central role in the water-driving soil erosion processes. An in-depth knowledge about temporal and spatial erosivity of rainfall events is required to gain a better understanding of soil erosion processes and optimize soil protection measures efficiency. In this study, the spatiotemporal distribution of more than 300,000 erosive events measured at 1181 locations, part of the Rainfall Erosivity Database at European Scale (REDES) database, is studied to shed some new light on the rainfall erosivity in Europe. Rainfall erosive events are statistically investigated through the Lorenz curve and derived coefficients such as the Gini coefficient (G). Additionally, seasonal characteristics of the most and the less erosive events are compared to investigate seasonal characteristics of rainstorms across Europe. The G shows largest values of inequality of the inter-annual temporal distribution of the rainfall erosive events in the Alpine region, mostly due to the large number of rainfall events with smaller rainfall erosivity. While for other parts of Europe, the inequality described by the G is mostly due to a small number of high erosive events. The G slightly decreases from south to north while no clear regional patterns can be detected. Additionally, in Europe, on average 11% (ranging from 1 to 24%) of all erosive events contribute to form 50% of the total rainfall erosivity. Furthermore, higher erosive rainfall events tend to occur later in the year compared to less erosive events that take place earlier. To our knowledge, this study is the first one addressing event scale rainfall erosivity distribution using more than 300,000 rainfall erosivity events and covering almost a whole continent. Scientifically our findings represent a major step towards large-scale process-based erosion modelling while, practically, they provide new elements that can support national and local soil erosion monitoring programs.

    [Διαβάστε περισσότερα...]
  • Raster και Raster Stack στην R

    Η ανάρτηση αυτή περιλαμβάνει το περιεχόμενο από το workshop στο FOSSCOMM 2019 στην Λαμία. Για περισσότερες λεπτομέρειες ανατρέξτε στο github repository.

    Στόχος του εργαστηρίου είναι η εξοικείωση του χρήστη με το πακέτο raster της R το οποίο προσφέρει την δυνατότητα ανάγνωσης ψηφιδωτών δεδομένων (raster) και επεξεργασίας τους (cropreclassifyreprojectresample κτλ.).

    Επιπλέον, θα επικεντρωθούμε στην κλάση raster stack η οποία δημιουργεί συστοιχίες ψηφιδωτών δεδομένων, κατάλληλες για χρονοσειρές και πολυκαναλικές εικόνες.

    Η διεξαγωγή του εργαστηρίου θα γίνει μέσω παραδειγμάτων και με την χρήση δεδομένων νυκτερινών φώτων DMSP-OLS Nighttime Lights Time Series (Stable Lights Version 4). Θα προηγηθεί μια σύντομη παρουσίαση των βημάτων και της διαδικασίας ώστε οι χρήστες να αποκτήσουν μια σύντομη αλλά περιεκτική εικόνα του στόχου του εργαστηρίου και των δυνατοτητων που προσφέρει ο προγραμματισμός με την R.

    [Διαβάστε περισσότερα...]
  • Λήψη δορυφορικών εικόνων στην R μέσω του GIBS API

    Nighttime Imagery (Day/Night Band, Enhanced Near Constant Contrast)

    Η υπηρεσία Global Imagery Browse Services (GIBS) του παρατηρητηρίου Earth Observing System Data and Information System (EOSDIS) παρέχει δορυφορικές εικόνες πλήρους ανάλυσης σε παγκόσμιο επίπεδο. Η υπηρεσία παρέχει πρόσβαση σε 600 και πλέον προϊόντα δορυφορικών εικόνων τα οποία στην πλειονότητά τους είναι διαθέσιμα με χρονική υστέρηση λίγων ωρών από την διέλευση του δορυφόρου. Η οπτικοποίηση των δεδομένων αυτών υποστηρίζεται από σχετική διαδικτυακή εφαρμογή WorldView  αλλά και μέσω υπηρεσιών (πρωτόκολλα OGC Web Map Tile Service (WMTS)/Tiled Web Map Service (TWMS)) που επιτρέπουν την ενσωμάτωσή του σε τρίτες εφαρμογές. Επιπλέον μια από της πλέον χρήσιμες λειτουργίες της υπηρεσίας GIBS είναι δυνατότητα εκτέλεσης σεναρίου εντολών μέσω της βιβλιοθήκης Geospatial Data Abstraction Library (GDAL). Αυτή η δυνατότητα παρέχει στον χρήση την ευκολία να αυτοματοποιήσει την λήψη δεδομένων (π.χ. μαζική λήψη δεδομένων, αυτόματη περιοδική λήψη, λήψη δεδομένων υπό συγκεκριμένες συνθήκες κτλ.).

    Στην τρέχουσα ανάρτηση παρουσιάζεται ένα σενάριο λήψης δεδομένων από την πλατφόρμα GIBS με την γλώσσα προγραμματισμού R. Στόχος είναι να δημιουργηθεί ένα σενάριο εντολών όπου ο χρήστης θα παρέχει στο script πληροφορίες για το προϊόν και το script θα λαμβάνει τις σχετικές δορυφορικές εικόνες.
    [Διαβάστε περισσότερα...]

  • Τρισδιάστατα μοντέλα για εκτύπωση βουνών της Ελλάδας

    Το μοντέλο του Υμηττού για 3d εκτύπωση.

    Στα πλαίσιο του έργου “Ηλεκτρονικές Υπηρεσίες για την Ανάπτυξη και Διάδοση του Ανοιχτού Λογισμικού” δημιουργήθηκαν τρισδιάστατα μοντέλα για βουνά της Ελλάδας, κατάλληλα για εκτύπωση σε 3D εκτυπωτή.
    Η ομάδα έργου αποτελείται από τους Δημήτρη Σταθάκη, Λεωνίδα Λιάκο, Ιωάννη Φαρασλή και Χρίστο Λυοτύρη.
    Τα συγκεκριμένα βουνά που μοντελοποιήθηκαν είναι ο Υμηττός, το Πήλιο και η Σαμοθράκη (όλο το νησί). Τα ψηφιακά μοντέλα που αναπαράχθηκαν διατίθενται υπό την άδεια Creative Commons Attribution-ShareAlike 4.0 International. [Διαβάστε περισσότερα...]

  • PyQgis tips: QgsVectorLayer από postgis sql ερώτημα

    Μία πολύ χρήσιμη δυνατότητα που μας προσφέρεται από τo pyqgis είναι να "φορτώνουμε" vector layer που προέρχεται από ερώτημα sql στην postgis.
    Έτσι ο χρήστης μπορεί να συνδυάσει την πληροφορία που θέλει με εντολές join και να φτιάξει το layer που επιθυμεί.
    Ωστόσο, στην σχετική τεκμηρίωση για το τρόπο που φορτώνουμε vector layers δεν αναφέρεται η σχετική δυνατότητα. Μετά από αρκετό ψάξιμο και δοκιμές κατέληξα σε μερικά tips για το πως γίνεται. [Διαβάστε περισσότερα...]

  • Υπολογισμός της καμπύλης Lorenz με το open source πακέτο στατιστικής R

    Καμπύλη Lorenz
    Ένα από τα μέτρα συγκέντρωσης για την εκτίμηση των περιφεριακών ανισοτήτων είναι η καμπύλη Lorenz.
    Για την εφαρμογή της καμπύλης Lorenz (Παπαδασκαλόπουλος, 2000), οι επιμέρους γεωγραφικές περιοχές με την έκταση και τον αντίστοιχο πληθυσμό τους, ταξινομούνται κατά φθίνουσα σειρά της πυκνότητας αυτών, και κάθε μιας των δυο τούτων στηλών, της εκτάσεως και του πληθυσμού, υπολογίζεται η ποσοστιαία  κατανομή, καθώς και η αθροιστική σειρά των ποσοστιαίων τούτων αναλογιών. Στη συνέχεια σε σύστημα ορθογώνιων συντεταγμένων χαράσσεται η γραφική παράσταση της καμπύλης Lorenz, στην οποία ο μεν άξονας των τετμημένων μετρά την αθροιστική σειρά της ποσοστιαίας κατανομής του πληθυσμού, ο δε άξονας των τεταγμένων της αθροιστική σειρά της ποσοστιαίας κατανομής της εκτάσεως και εντοπίζονται τα σημεία της αντιστοιχίας των αθροιστικών σειρών (%) εκτάσεως και πληθυσμού, τα οποία σχηματίζουν πολυγωνική γραμμή ο (καλούμενη καμπύλη Lorenz), φέρεται δε και η διαγώνιος (45ο), καλούμενη διαγώνιος ισοκατανομής. Στην απεικόνιση αυτή παρατηρούνται τα εξής:

    1. αν ο πληθυσμός ήταν ομοιομόρφως κατανεμημένος σε όλες τις επιμέρους περιοχές (γεωγραφικά διαμερίσματα), η αθροιστική καμπύλη του πληθυσμού θα ακολουθούσε τη διαγώνιο γραμμή της αθροιστικής σειράς των περιοχών (διαγώνιος ισοκατανομής).
    2. εάν, αντιθέτως ο πληθυσμός ήταν κατανεμημένος με την μέγιστη δυνατή ανισότητα, τότε η αθροιστική καμπύλη του πληθυσμού θα συνέπιπτε με τον άξονα των τετμημένων.
    3. ο βαθμός συγκεντρώσεως αντιπροσωπεύεται κατ’ ακολουθία από την απόκλιση της καμπύλης μεταξύ των δύο αυτών άκρων.

    Στα πλαίσια του τρέχοντος άρθρου και για εκπαιδευτικούς λόγους θα υπολογιστεί η καμπύλη Lorenz για τον νομό Θεσπρωτίας το 2001 με την βοήθεια του στατιστικού πακέτου R.
    [Διαβάστε περισσότερα...]

  • Υπολογισμός δεικτών χωρικής κεντρικότητας με την R

    Οι κεντρογραφικές τεχνικές αποτελούν τις βασικές τεχνικές περιγραφικής στατιστικής της χωρικής κατανομής.
    Ο πιο απλός περιγραφικός δείκτης μιας κατανομής είναι ο χωρικός μέσος (mean center). Ουσιαστικά πρόκειται για τον μέσο όρο των συντεταγμένων Χ,Υ μιας σημειακής κατανομής.
    Συχνά αποκαλείται και κεντροειδές.

    Δίνεται από τον τύπο:

    όπου,

    οι συντεταγμένες του χωρικού μέσου δηλαδή ο μέσος όρος των συντεταγμένων του πλήθους των σημείων της κατανομής.
    Ταυτόχρονα με τον αριθμητικό χωρικό μέσο, μπορεί να υπολογιστεί ο σταθμισμένος χωρικός μέσος (Weighted Mean Center) όπου ο υπολογισμός γίνεται με το στάθμισμα κάθε συντεταγμένης με μία μεταβλητή π.χ. εισόδημα, εγκληματικότητα, πληθυσμός κτλ.

    Δίνεται από τον τύπο:


    όπου,
    οι συντεταγμένες του σταθμισμένου χωρικού μέσου δηλαδή ο μέσος όρος των διακριτών συντεταγμένων του πλήθους των σημείων της κατανομής σταθμισμένες με το βάρος του κάθε σημείου.

    [Διαβάστε περισσότερα...]

  • Οπτική βελτίωση σκιασμένου αναγλύφου εδάφους

    Παρακάτω προτείνεται μια μεθοδολογία με εργαλεία ΕΛ/ΛΑΚ, για την οπτική βελτίωση μιας σκίασης αναγλύφου εδάφους.

    Το λογισμικό που απαιτείται είναι το gdaldem, gdalwrap και το GIMP

    Τα κλασσικά προγράμματα ΣΓΠ αναπαριστούν την σκίαση αναγλύφου με τρόπο μονότονο, συχνά μη φυσικό και με περιττή πληροφορία για τον αναγνώστη.

    Στόχος είναι να εφαρμόσουμε μια τεχνική που πρώτος εφάρμοσε ο Eduard Imhof, την αναπαράσταση του αναγλύφου με θέαση αφ' ύψηλου (aerial perspective). Η συγκεκριμένη τεχνική του Imhof προσομοιάζει την πραγματική, αφ' υψηλού οπτική ενός τοπίου όπου ο θεατής βλέπει τις κορυφές με περισσότερη λεπτομέρεια και έντονη αντίθεση και τις περιοχές χαμηλότερου υψομέτρου με λιγότερη λεπτομέρεια λόγω απόστασης και ομίχλης στο τοπίο. Η τεχνική που ακολουθεί απαιτεί:

    1. Ένα Ψηφιακό Μοντέλο Εδάφους (DEM)

    2. Σκίαση του ανωτέρω DEM με διπλάσια κάθετη υπερβολή του υψομέτρου (z exaggeration).

    3. Σκίαση του ανωτέρω DEM με διπλάσια κάθετη υπερβολή του υψομέτρου (z exaggeration) και πηγή φωτός σε γωνία 90 μοιρών.

    4. Τεχνικές θαμπώματος (blur) για την εξομάλυνση και την απαλοιφή της περιττής πληροφορίας του DEM.

    5. Οπτική ενίσχυση των ορεινών όγκων και μίξη του αποτελέσματος με την σκίαση.

    [Διαβάστε περισσότερα...]

  • Online Σύστημα Γεωγραφικών Πληροφοριών Δήμου Σουλίου

    Η εν λόγω πλατφόρμα αποτελεί ένα Σύστημα Γεωγραφικών Πληροφοριών (ΣΓΠ-GIS) για τον δήμο Σουλίου σε πρώιμο στάδιο. Η λειτουργία του βασίζεται στο λογισμικό ανοικτού κώδικα  WorldMap που είναι μια πλατφόρμα τεχνολογίας διαδικτύου σε συνδυασμό με ΣΓΠ. Η πλατφόρμα προσφέρει την τεχνολογική υποδομή μέσω της οποίας μπορεί να στηριχθεί η αποτύπωση κάθε είδους χωρικών δεδομένων (διανυσματικών και ψηφιδωτών) και να χαρτογραφηθούν πολύτιμα δεδομένα που είναι αποσπασματικά αποθηκευμένα σε έγγραφα, χάρτες, συσκευές GPS, ιδιωτικές συλλογές και δημόσια αρχεία. Σκοπεύοντας στην συγχώνευση γεωχωρικών δεδομένων από διαφορετικές πηγές και μέσα αποτελεί ένα εύχρηστο εργαλείο σχεδιασμού για τις τοπικές αρχές αυτοδιοίκησης αλλά και μία ολοκληρωμένη βάση δεδομένων για τον μεμονωμένο επιστήμονα, ερευνητή. Αναμφισβήτητη δε είναι και η συνεισφορά του στον επισκέπτη αυτού του τόπου εφόσον μπορεί ταυτόχρονα να λειτουργήσει σαν τουριστικός χάρτης σε όποιον περιηγητή επιθυμεί να γνωρίσεις τις ομορφιές του δήμου. Οι λειτουργίες του όμως δεν περιορίζονται στην απλή εποπτεία των χωρικών δεδομένων αλλά με συμμετοχικό τρόπο κάθε χρήστης-επισκέπτης της πλατφόρμας μπορεί να συμβάλει με την προσθήκη/διόρθωση δεδομένων σε αυτή. Μέσω ενός συστήματος καταχώρησης δικαιωμάτων, ο διαχειριστής επιλέγει το επίπεδο πρόσβασης στην επεξεργασία και προσπέλασης των χωρικών δεδομένων σε ομάδες χρηστών. Παράλληλα με την απόδοση των χωρικών δεδομένων η πλατφόρμα προσφέρει σε αρκετές περιπτώσεις πληθώρα περιγραφικών δεδομένων για τα γεωγραφικά αντικείμενα καθώς και παραπομπές σε ιστοσελίδες σχετικές με αυτά. Όσον αφορά την πηγή των δεδομένων αυτά προέρχονται από πηγές που διαθέτουν ελεύθερα δεδομένα όπως το geodata.gov.gr αλλά και από επεξεργασία δεδομένων από το προσωπικό αρχείο του διαχειριστή. [Διαβάστε περισσότερα...]

Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported
This work by Leonidas Liakos is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported.