Institut für Geographie und Geologie
Refine
Has Fulltext
- yes (363)
Year of publication
Document Type
- Journal article (208)
- Doctoral Thesis (100)
- Book (35)
- Book article / Book chapter (6)
- Master Thesis (6)
- Report (5)
- Working Paper (2)
- Habilitation (1)
Keywords
- remote sensing (39)
- Einzelhandel (30)
- Fernerkundung (23)
- climate change (18)
- Deutschland (15)
- MODIS (14)
- Sentinel-2 (14)
- time series (14)
- Germany (12)
- Nachhaltigkeit (12)
Institute
- Institut für Geographie und Geologie (363)
- Theodor-Boveri-Institut für Biowissenschaften (9)
- Institut für Altertumswissenschaften (4)
- Graduate School of Science and Technology (3)
- Institut für Informatik (2)
- Frauenklinik und Poliklinik (1)
- Institut für Klinische Epidemiologie und Biometrie (1)
- Julius-von-Sachs-Institut für Biowissenschaften (1)
- Medizinische Klinik und Poliklinik I (1)
- Neuphilologisches Institut - Moderne Fremdsprachen (1)
Schriftenreihe
Sonstige beteiligte Institutionen
- The Italian Federation of Parks and Nature Reserves (3)
- ALPARC - The Alpine Network of Protected Areas (2)
- Eurac research (2)
- Salzburg Institute for Regional Planning and Housing (2)
- Urban Planning Institute of the Republic of Slovenia (2)
- ALPARC - Das Netzwerk Alpiner Schutzgebiete (1)
- Albert-Ludwigs-Universität Freiburg (1)
- Deutsches Klimaservice Zentrum (GERICS) (1)
- Deutsches Zentrum für Luft & Raumfahrt (DLR) (1)
- Deutsches Zentrum für Luft- und Raumfahrt (1)
ResearcherID
- I-5818-2014 (1)
EU-Project number / Contract (GA) number
- 308377 (2)
- 20-3044-2-11 (1)
- 227159 (1)
- 243964 (1)
- 714087 (1)
- 776019 (1)
- 818182 (1)
- 834709 (1)
- LIFE12 BIO/AT/000143 (1)
- LIFE20 NAT/AT/000049 (1)
Forests are essential for global environmental well-being because of their rich provision of ecosystem services and regulating factors. Global forests are under increasing pressure from climate change, resource extraction, and anthropologically-driven disturbances. The results are dramatic losses of habitats accompanied with the reduction of species diversity. There is the urgent need for forest biodiversity monitoring comprising analysis on α, β, and γ scale to identify hotspots of biodiversity. Remote sensing enables large-scale monitoring at multiple spatial and temporal resolutions. Concepts of remotely sensed spectral diversity have been identified as promising methodologies for the consistent and multi-temporal analysis of forest biodiversity. This review provides a first time focus on the three spectral diversity concepts “vegetation indices”, “spectral information content”, and “spectral species” for forest biodiversity monitoring based on airborne and spaceborne remote sensing. In addition, the reviewed articles are analyzed regarding the spatiotemporal distribution, remote sensing sensors, temporal scales and thematic foci. We identify multispectral sensors as primary data source which underlines the focus on optical diversity as a proxy for forest biodiversity. Moreover, there is a general conceptual focus on the analysis of spectral information content. In recent years, the spectral species concept has raised attention and has been applied to Sentinel-2 and MODIS data for the analysis from local spectral species to global spectral communities. Novel remote sensing processing capacities and the provision of complementary remote sensing data sets offer great potentials for large-scale biodiversity monitoring in the future.
Snow is a vital environmental parameter and dynamically responsive to climate change, particularly in mountainous regions. Snow cover can be monitored at variable spatial scales using Earth Observation (EO) data. Long-lasting remote sensing missions enable the generation of multi-decadal time series and thus the detection of long-term trends. However, there have been few attempts to use these to model future snow cover dynamics. In this study, we, therefore, explore the potential of such time series to forecast the Snow Line Elevation (SLE) in the European Alps. We generate monthly SLE time series from the entire Landsat archive (1985–2021) in 43 Alpine catchments. Positive long-term SLE change rates are detected, with the highest rates (5–8 m/y) in the Western and Central Alps. We utilize this SLE dataset to implement and evaluate seven uni-variate time series modeling and forecasting approaches. The best results were achieved by Random Forests, with a Nash–Sutcliffe efficiency (NSE) of 0.79 and a Mean Absolute Error (MAE) of 258 m, Telescope (0.76, 268 m), and seasonal ARIMA (0.75, 270 m). Since the model performance varies strongly with the input data, we developed a combined forecast based on the best-performing methods in each catchment. This approach was then used to forecast the SLE for the years 2022–2029. In the majority of the catchments, the shift of the forecast median SLE level retained the sign of the long-term trend. In cases where a deviating SLE dynamic is forecast, a discussion based on the unique properties of the catchment and past SLE dynamics is required. In the future, we expect major improvements in our SLE forecasting efforts by including external predictor variables in a multi-variate modeling approach.
Drought is a recurring natural climatic hazard event over terrestrial land; it poses devastating threats to human health, the economy, and the environment. Given the increasing climate crisis, it is likely that extreme drought phenomena will become more frequent, and their impacts will probably be more devastating. Drought observations from space, therefore, play a key role in dissimilating timely and accurate information to support early warning drought management and mitigation planning, particularly in sparse in-situ data regions. In this paper, we reviewed drought-related studies based on Earth observation (EO) products in Southeast Asia between 2000 and 2021. The results of this review indicated that drought publications in the region are on the increase, with a majority (70%) of the studies being undertaken in Vietnam, Thailand, Malaysia and Indonesia. These countries also accounted for nearly 97% of the economic losses due to drought extremes. Vegetation indices from multispectral optical remote sensing sensors remained a primary source of data for drought monitoring in the region. Many studies (~21%) did not provide accuracy assessment on drought mapping products, while precipitation was the main data source for validation. We observed a positive association between spatial extent and spatial resolution, suggesting that nearly 81% of the articles focused on the local and national scales. Although there was an increase in drought research interest in the region, challenges remain regarding large-area and long time-series drought measurements, the combined drought approach, machine learning-based drought prediction, and the integration of multi-sensor remote sensing products (e.g., Landsat and Sentinel-2). Satellite EO data could be a substantial part of the future efforts that are necessary for mitigating drought-related challenges, ensuring food security, establishing a more sustainable economy, and the preservation of the natural environment in the region.
Accurate crop monitoring in response to climate change at a regional or field scale
plays a significant role in developing agricultural policies, improving food security,
forecasting, and analysing global trade trends. Climate change is expected to
significantly impact agriculture, with shifts in temperature, precipitation patterns, and
extreme weather events negatively affecting crop yields, soil fertility, water availability,
biodiversity, and crop growing conditions. Remote sensing (RS) can provide valuable
information combined with crop growth models (CGMs) for yield assessment by
monitoring crop development, detecting crop changes, and assessing the impact of
climate change on crop yields. This dissertation aims to investigate the potential of RS
data on modelling long-term crop yields of winter wheat (WW) and oil seed rape (OSR)
for the Free State of Bavaria (70,550 km2
), Germany. The first chapter of the dissertation
describes the reasons favouring the importance of accurate crop yield predictions for
achieving sustainability in agriculture. Chapter second explores the accuracy
assessment of the synthetic RS data by fusing NDVIs of two high spatial resolution data
(high pair) (Landsat (30 m, 16-days; L) and Sentinel-2 (10 m, 5–6 days; S), with four low
spatial resolution data (low pair) (MOD13Q1 (250 m, 16-days), MCD43A4 (500 m, one
day), MOD09GQ (250 m, one-day), and MOD09Q1 (250 m, 8-days)) using the spatial
and temporal adaptive reflectance fusion model (STARFM), which fills regions' cloud
or shadow gaps without losing spatial information. The chapter finds that both L-MOD13Q1 (R2 = 0.62, RMSE = 0.11) and S-MOD13Q1 (R2 = 0.68, RMSE = 0.13) are more
suitable for agricultural monitoring than the other synthetic products fused. Chapter
third explores the ability of the synthetic spatiotemporal datasets (obtained in chapter
2) to accurately map and monitor crop yields of WW and OSR at a regional scale. The
chapter investigates and discusses the optimal spatial (10 m, 30 m, or 250 m), temporal
(8 or 16-day) and CGMs (World Food Studies (WOFOST), and the semi-empiric light
use efficiency approach (LUE)) for accurate crop yield estimations of both crop types.
Chapter third observes that the observations of high temporal resolution (8-day)
products of both S-MOD13Q1 and L-MOD13Q1 play a significant role in accurately
measuring the yield of WW and OSR. The chapter investigates that the simple light use
efficiency (LUE) model (R2 = 0.77 and relative RMSE (RRMSE) = 8.17%) that required fewer input parameters to simulate crop yield is highly accurate, reliable, and more
precise than the complex WOFOST model (R2 = 0.66 and RRMSE = 11.35%) with higher
input parameters. Chapter four researches the relationship of spatiotemporal fusion
modelling using STRAFM on crop yield prediction for WW and OSR using the LUE
model for Bavaria from 2001 to 2019. The chapter states the high positive correlation
coefficient (R) = 0.81 and R = 0.77 between the yearly R2 of synthetic accuracy and
modelled yield accuracy for WW and OSR from 2001 to 2019, respectively. The chapter
analyses the impact of climate variables on crop yield predictions by observing an
increase in R2
(0.79 (WW)/0.86 (OSR)) and a decrease in RMSE (4.51/2.57 dt/ha) when
the climate effect is included in the model. The fifth chapter suggests that the coupling
of the LUE model to the random forest (RF) model can further reduce the relative root
mean square error (RRMSE) from -8% (WW) and -1.6% (OSR) and increase the R2 by
14.3% (for both WW and OSR), compared to results just relying on LUE. The same
chapter concludes that satellite-based crop biomass, solar radiation, and temperature
are the most influential variables in the yield prediction of both crop types. Chapter six
attempts to discuss both pros and cons of RS technology while analysing the impact of
land use diversity on crop-modelled biomass of WW and OSR. The chapter finds that
the modelled biomass of both crops is positively impacted by land use diversity to the
radius of 450 (Shannon Diversity Index ~0.75) and 1050 m (~0.75), respectively. The
chapter also discusses the future implications by stating that including some dependent
factors (such as the management practices used, soil health, pest management, and
pollinators) could improve the relationship of RS-modelled crop yields with
biodiversity. Lastly, chapter seven discusses testing the scope of new sensors such as
unmanned aerial vehicles, hyperspectral sensors, or Sentinel-1 SAR in RS for achieving
accurate crop yield predictions for precision farming. In addition, the chapter highlights
the significance of artificial intelligence (AI) or deep learning (DL) in obtaining higher
crop yield accuracies.
Spatiotemporal Fusion Modelling Using STARFM: Examples of Landsat 8 and Sentinel-2 NDVI in Bavaria
(2022)
The increasing availability and variety of global satellite products provide a new level of data with different spatial, temporal, and spectral resolutions; however, identifying the most suited resolution for a specific application consumes increasingly more time and computation effort. The region’s cloud coverage additionally influences the choice of the best trade-off between spatial and temporal resolution, and different pixel sizes of remote sensing (RS) data may hinder the accurate monitoring of different land cover (LC) classes such as agriculture, forest, grassland, water, urban, and natural-seminatural. To investigate the importance of RS data for these LC classes, the present study fuses NDVIs of two high spatial resolution data (high pair) (Landsat (30 m, 16 days; L) and Sentinel-2 (10 m, 5–6 days; S), with four low spatial resolution data (low pair) (MOD13Q1 (250 m, 16 days), MCD43A4 (500 m, one day), MOD09GQ (250 m, one-day), and MOD09Q1 (250 m, eight day)) using the spatial and temporal adaptive reflectance fusion model (STARFM), which fills regions’ cloud or shadow gaps without losing spatial information. These eight synthetic NDVI STARFM products (2: high pair multiply 4: low pair) offer a spatial resolution of 10 or 30 m and temporal resolution of 1, 8, or 16 days for the entire state of Bavaria (Germany) in 2019. Due to their higher revisit frequency and more cloud and shadow-free scenes (S = 13, L = 9), Sentinel-2 (overall R\(^2\) = 0.71, and RMSE = 0.11) synthetic NDVI products provide more accurate results than Landsat (overall R\(^2\) = 0.61, and RMSE = 0.13). Likewise, for the agriculture class, synthetic products obtained using Sentinel-2 resulted in higher accuracy than Landsat except for L-MOD13Q1 (R\(^2\) = 0.62, RMSE = 0.11), resulting in similar accuracy preciseness as S-MOD13Q1 (R\(^2\) = 0.68, RMSE = 0.13). Similarly, comparing L-MOD13Q1 (R\(^2\) = 0.60, RMSE = 0.05) and S-MOD13Q1 (R\(^2\) = 0.52, RMSE = 0.09) for the forest class, the former resulted in higher accuracy and precision than the latter. Conclusively, both L-MOD13Q1 and S-MOD13Q1 are suitable for agricultural and forest monitoring; however, the spatial resolution of 30 m and low storage capacity makes L-MOD13Q1 more prominent and faster than that of S-MOD13Q1 with the 10-m spatial resolution.
Grasslands shape many landscapes of the earth as they cover about one-third of its surface. They are home and provide livelihood for billions of people and are mainly used as source of forage for animals. However, grasslands fulfill many additional ecosystem functions next to fodder production, such as storage of carbon, water filtration, provision of habitats and cultural values. They play a role in climate change (mitigation) and in preserving biodiversity and ecosystem functions on a global scale. The degree to what these ecosystem functions are present within grassland ecosystems is largely determined by the management. Individual management practices and the use intensity influence the species composition as well as functions, like carbon storage, while higher use intensities (e.g. high mowing frequencies) usually show a negative impact. Especially in Central European countries, like in Germany, the determining influence of grassland management on its physiognomy and ecosystem functions leads to a large variability and small-scale alternations of grassland parcels. Large-scale information on the management and use intensity of grasslands is not available. Consequently, estimations of grassland ecosystem functions are challenging which, however, would be required for large-scale assessments of the status of grassland ecosystems and optimized management plans for the future. The topic of this thesis tackles this gap by investigating the major grassland management practice in Germany, which is mowing, for multiple years, in high spatial resolution
and on a national scale.
Earth Observation (EO) has the advantage of providing information of the earth’s surface on multi-temporal time steps. An extensive literature review on the use of EO for grassland management and production analyses, which was part of this thesis, showed that in particular research on grasslands consisting of small parcels with a large variety of management and use intensity, like common in Central Europe, is underrepresented. Especially
the launch of the Sentinel satellites in the recent past now enables the analyses of such grasslands due to their high spatial and temporal resolution. The literature review specifically on the investigation of grassland mowing events revealed that most previous studies focused on small study areas, were exploratory, only used one sensor type and/or lacked a reference data set with a complete range of management options.
Within this thesis a novel framework to detect grassland mowing events over large areas is presented which was applied and validated for the entire area of Germany for multiple years (2018–2021). The potential of both sensor types, optical (Sentinel-2) and Synthetic Aperture Radar (SAR) (Sentinel-1) was investigated regarding grassland mowing event detection. Eight EO parameters were investigated, namely the Enhanced Vegetation Index (EVI), the backscatter intensity and the interferometric (InSAR) temporal coherence for both available polarization modes (VV and VH), and the polarimetric (PolSAR) decomposition parameters Entropy, K0 and K1. An extensive reference data set was generated based on daily images of webcams distributed in Germany which resulted in mowing information
for grasslands with the entire possible range of mowing frequencies – from one to six in Germany – and in 1475 reference mowing events for the four years of interest.
For the first time a observation-driven mowing detection approach including data from Sentinel-2 and Sentinel-1 and combining the two was developed, applied and validated on large scale. Based on a subset of the reference data (13 grassland parcels with 44 mowing events) from 2019 the EO parameters were investigated and the detection algorithm
developed and parameterized. This analysis showed that a threshold-based change detection approach based on EVI captured grassland mowing events best, which only failed during periods of clouds. All SAR-based parameters showed a less consistent behavior to mowing events, with PolSAR Entropy and InSAR Coherence VH, however, revealing the
highest potential among them. A second, combined approach based on EVI and a SARbased parameter was developed and tested for PolSAR Entropy and InSAR VH. To avoid additional false positive detections during periods in which mowing events are anyhow reliably detected using optical data, the SAR-based mowing detection was only initiated
during long gaps within the optical time series (< 25 days). Application and validation of
these approaches in a focus region revealed that only using EVI leads to the highest accuracies (F1-Score = 0.65) as combining this approach with SAR-based detection led to a strong increase in falsely detected mowing events resulting in a decrease of accuracies (EVI + PolSAR ENT F1-Score = 0.61; EVI + InSAR COH F1-Score = 0.61).
The mowing detection algorithm based on EVI was applied for the entire area of Germany for the years 2018-2021. It was revealed that the largest share of grasslands with high mowing frequencies (at least four mowing events) can be found in southern/south-eastern Germany. Extensively used grassland (mown up to two times) is distributed within the entire country with larger shares in the center and north-eastern parts of Germany. These patterns stay constant in general, but small fluctuations between the years are visible. Early mown grasslands can be found in southern/south-eastern Germany – in line with high mowing frequency areas – but also in central-western parts. The years 2019 and 2020 revealed higher accuracies based on the 1475 mowing events of the multi-annual validation data set
(F1-Scores of 0.64 and 0.63), 2018 and 2021 lower ones (F1-Score of 0.52 and 0.50).
Based on this new, unprecedented data set, potential influencing factors on the mowing dynamics were investigated. Therefore, climate, topography, soil data and information on conservation schemes were related to mowing dynamics for the year 2020, which showed a high number of valid observations and detection accuracy. It was revealed that there are no strong linear relationships between the mowing frequency or the timing of the first mowing event and the investigated variables. However, it was found that for intensive grassland usage certain climatic and topographic conditions have to be fulfilled, while extensive grasslands appear on the entire spectrum of these variables. Further, higher mowing frequencies occur on soils with influence of ground water and lower mowing frequencies in protected areas. These results show the complex interplay between grassland mowing dynamics and external influences and highlight the challenges of policies aiming to protect grassland ecosystem functions and their need to be adapted to regional circumstances.
Individual mobility and human patterns analyses is receiving increasing attention in numerous interdisciplinary studies and publications using the concept of time-geography but is largely unknown to the subdiscipline of sports geography. Meanwhile the visualization and evaluation of large data of individual patterns are still a major challenge. While a qualitative, microscale view on spatial-temporal topics is more common in today's pattern research using mostly 24h time intervals, this work examines a quantitative approach focusing on an extended period of life. This paper presents a combination of time-geographic approaches with 3D-geoinformation systems and demonstrates their value for analysing individual mobility by implementing a path-homogeneity factor (HPA). Using the example of professional athletes, it is shown which groups display greater similarities in their career paths. While a high homogeneity suggests that groups make similar decisions through socially influenced processes, low values allow the assumption that external processes provide stronger, independent individual structures.
Die städtische Umwelt ist in steter Veränderung, vor allem durch den Bau, aber auch durch die Zerstörung von städtischen Elementen. Die formelle Entwicklung ist ein Prozess mit langen Planungszeiträumen und die bebaute Landschaft wirkt daher statisch. Dagegen unterliegen informelle oder spontane Siedlungen aufgrund ihrer stets unvollendeten städtischen Form einer hohen Dynamik – so wird in der Literatur berichtet. Allerdings sind Dynamik und die morphologischen Merkmale der physischen Transformation in solchen Siedlungen, die städtische Armut morphologisch repräsentieren, auf globaler Ebene bisher kaum mit einer konsistenten Datengrundlage empirisch untersucht worden. Hier setzt die vorliegende Arbeit an. Unter der Annahme, dass die erforschte zeitliche Dynamik in Europa geringer ausfällt, stellt sich die generelle Frage nach einer katalogisierten Erfassung physischer Wohnformen von Armut speziell in Europa. Denn Wohnformen der Armut werden oft ausschließlich mit dem ‚Globalen Süden‘ assoziiert, insbesondere durch die Darstellung von Slums. Tatsächlich ist Europa sogar die Wiege der Begriffe ‚Slum‘ und ‚Ghetto‘, die vor Jahrhunderten zur Beschreibung von Missständen und Unterdrückung auftauchten. Bis heute weist dieser facettenreiche Kontinent eine enorme Vielfalt an physischen Wohnformen der Armut auf, die ihre Wurzeln in unterschiedlichen Politiken, Kulturen, Geschichten und Lebensstilen haben. Um über diese genannten Aspekte Aufschluss zu erlangen, bedarf es u.a. der Bildanalyse durch Satellitenbilder. Diese Arbeit wird daher mittels Fernerkundung bzw. Erdbeobachtung (EO) sowie zusätzlicher Literaturrecherchen und einer empirischen Erhebung erstellt. Um Unsicherheiten konzeptionell und in der Erfassung offenzulegen, ist die Methode der manuellen Bildinterpretation von Armutsgebieten kritisch zu hinterfragen.
Das übergeordnete Ziel dieser Arbeit ist eine bessere Wissensbasis über Armut zu schaffen, um Maßnahmen zur Reduzierung von Armut entwickeln zu können. Die Arbeit dient dabei als eine Antwort auf die Nachhaltigkeitsziele der Vereinten Nationen. Es wird Grundlagenforschung betrieben, indem Wissenslücken in der Erdbeobachtung zu physisch-baulichen bzw. morphologischen Erscheinungen von Armut auf Gebäude-Ebene explorativ analysiert werden. Die Arbeit wird in drei Forschungsthemen bzw. Studienteile untergliedert:
Ziel des ersten Studienteils ist die globale raumzeitliche Erfassung von Dynamiken durch Anknüpfung an bisherige Kategorisierungen von Armutsgebieten. Die bisherige Wissenslücke soll gefüllt werden, indem über einen Zeitraum von etwa sieben Jahren in 16 dokumentierten Manifestationen städtischer Armut anhand von Erdbeobachtungsdaten eine zeitliche Analyse der bebauten Umwelt durchgeführt wird. Neben einer global verteilten Gebietsauswahl wird die visuelle Bildinterpretation (MVII) unter Verwendung von hochauflösenden optischen Satellitendaten genutzt. Dies geschieht in Kombination mit in-situ- und Google Street View-Bildern zur Ableitung von 3D-Stadtmodellen. Es werden physische Raumstrukturen anhand von sechs räumlichen morphologischen Variablen gemessen: Anzahl, Größe, Höhe, Ausrichtung und Dichte der Gebäude sowie Heterogenität der Bebauung. Diese ‚temporale Analyse‘ zeigt zunächst sowohl inter- als auch intra-urbane Unterschiede. Es lassen sich unterschiedliche, aber generell hohe morphologische Dynamiken zwischen den Untersuchungsgebieten finden. Dies drückt sich in vielfältiger Weise aus: von abgerissenen und rekonstruierten Gebieten bis hin zu solchen, wo Veränderungen innerhalb der gegebenen Strukturen auftreten. Geographisch gesehen resultiert in der Stichprobe eine fortgeschrittene Dynamik, insbesondere in Gebieten des Globalen Südens. Gleichzeitig lässt sich eine hohe räumliche Variabilität der morphologischen Transformationen innerhalb der untersuchten Gebiete beobachten. Trotz dieser teilweise hohen morphologischen Dynamik sind die räumlichen Muster von Gebäudefluchten, Straßen und Freiflächen überwiegend konstant. Diese ersten Ergebnisse deuten auf einen geringen Wandel in Europa hin, weshalb diese europäischen Armutsgebiete im folgenden Studienteil von Grund auf erhoben und kategorisiert werden.
Ziel des zweiten Studienteils ist die Erschaffung einer neuen Kategorisierung, speziell für das in der Wissenschaft unterrepräsentierte Europa. Die verschiedenen Formen nicht indizierter Wohnungsmorphologien werden erforscht und kategorisiert, um das bisherige globale wissenschaftliche ontologische Portfolio für Europa zu erweitern. Hinsichtlich dieses zweiten Studienteils bietet eine Literaturrecherche mit mehr als 1.000 gesichteten Artikeln die weitere Grundlage für den folgenden Fokus auf Europa. Auf der Recherche basierend werden mittels der manuellen visuellen Bildinterpretation (engl.: MVII) erneut Satellitendaten zur Erfassung der physischen Morphologien von Wohnformen genutzt. Weiterhin kommen selbst definierte geographische Indikatoren zu Lage, Struktur und formellem Status zum Einsatz. Darüber hinaus werden gesellschaftliche Hintergründe, die durch Begriffe wie ‚Ghetto‘, ‚Wohnwagenpark‘, ‚ethnische Enklave‘ oder ‚Flüchtlingslager‘ beschrieben werden, recherchiert und implementiert. Sie sollen als Erklärungsansatz für Armutsviertel in Europa dienen. Die Stichprobe der europäischen, insgesamt aber unbekannten Grundgesamtheit verdeutlicht eine große Vielfalt an physischen Formen: Es wird für Europa eine neue Kategorisierung von sechs Hauptklassen entwickelt, die von ‚einfachsten Wohnstätten‘ (z. B. Zelten) über ‚behelfsmäßige Unterkünfte ‘ (z. B. Baracken, Container) bis hin zu ‚mehrstöckigen Bauten‘ - als allgemeine Taxonomie der Wohnungsnot in Europa - reicht. Die Untersuchung zeigt verschiedene Wohnformen wie z. B. unterirdische oder mobile Typen, verfallene Wohnungen oder große Wohnsiedlungen, die die Armut im Europa des 21. Jahrhunderts widerspiegeln. Über die Wohnungsmorphologie hinaus werden diese Klassen durch die Struktur und ihren rechtlichen Status beschrieben - entweder als geplante oder als organisch-gewachsene bzw. weiterhin als formelle, informelle oder hybride (halblegale) Formen. Geographisch lassen sich diese ärmlichen Wohnformen sowohl in städtischen als auch in ländlichen Gebieten finden, mit einer Konzentration in Südeuropa. Der Hintergrund bei der Mehrheit der Morphologien betrifft Flüchtlinge, ethnische Minderheiten und sozioökonomisch benachteiligte Menschen - die ‚Unterprivilegierten‘.
Ziel des dritten Studienteils ist eine kritische Analyse der Methode. Zur Erfassung all dieser Siedlungen werden heutzutage Satellitenbilder aufgrund der Fortschritte bei den Bildklassifizierungsmethoden meist automatisch ausgewertet. Dennoch spielt die MVII noch immer eine wichtige Rolle, z.B. um Trainingsdaten für Machine-Learning-Algorithmen zu generieren oder für Validierungszwecke. In bestimmten städtischen Umgebungen jedoch, z.B. solchen mit höchster Dichte und struktureller Komplexität, fordern spektrale und textur-basierte Verflechtungen von überlappenden Dachstrukturen den menschlichen Interpreten immer noch heraus, wenn es darum geht einzelne Gebäudestrukturen zu erfassen. Die kognitive Wahrnehmung und die Erfahrung aus der realen Welt sind nach wie vor unumgänglich. Vor diesem Hintergrund zielt die Arbeit methodisch darauf ab, Unsicherheiten speziell bei der Kartierung zu quantifizieren und zu interpretieren. Kartiert werden Dachflächen als ‚Fußabdrücke‘ solcher Gebiete. Der Fokus liegt dabei auf der Übereinstimmung zwischen mehreren Bildinterpreten und welche Aspekte der Wahrnehmung und Elemente der Bildinterpretation die Kartierung beeinflussen. Um letztlich die Methode der MVII als drittes Ziel selbstkritisch zu reflektieren, werden Experimente als sogenannte ‚Unsicherheitsanalyse‘ geschaffen. Dabei digitalisieren zehn Testpersonen bzw. Probanden/Interpreten sechs komplexe Gebiete. Hierdurch werden quantitative Informationen über räumliche Variablen von Gebäuden erzielt, um systematisch die Konsistenz und Kongruenz der Ergebnisse zu überprüfen. Ein zusätzlicher Fragebogen liefert subjektive qualitative Informationen über weitere Schwierigkeiten. Da die Grundlage der hierfür bisher genutzten Kategorisierungen auf der subjektiven Bildinterpretation durch den Menschen beruht, müssen etwaige Unsicherheiten und damit Fehleranfälligkeiten offengelegt werden. Die Experimente zu dieser Unsicherheitsanalyse erfolgen quantifiziert und qualifiziert. Es lassen sich generell große Unterschiede zwischen den Kartierungsergebnissen der Probanden, aber eine hohe Konsistenz der Ergebnisse bei ein und demselben Probanden feststellen. Steigende Abweichungen korrelieren mit einer steigenden baustrukturellen (morphologischen) Komplexität. Ein hoher Grad an Individualität bei den Probanden äußert sich in Aspekten wie z.B. Zeitaufwand beim Kartieren, in-situ Vorkenntnissen oder Vorkenntnissen beim Umgang mit Geographischen Informationssystemen (GIS). Nennenswert ist hierbei, dass die jeweilige Datenquelle das Kartierungsverfahren meist beeinflusst. Mit dieser Studie soll also auch an der Stelle der angewandten Methodik eine weitere Wissenslücke gefüllt werden. Die bisherige Forschung komplexer urbaner Areale unter Nutzung der manuellen Bildinterpretation implementiert oftmals keine Unsicherheitsanalyse oder Quantifizierung von Kartierungsfehlern. Fernerkundungsstudien sollten künftig zur Validierung nicht nur zweifelsfrei auf MVII zurückgreifen können, sondern vielmehr sind Daten und Methoden notwendig, um Unsicherheiten auszuschließen.
Zusammenfassend trägt diese Arbeit zur bisher wenig erforschten morphologischen Dynamik von Armutsgebieten bei. Es werden inter- wie auch intra-urbane Unterschiede auf globaler Ebene präsentiert. Dabei sind allgemein hohe morphologische Transformationen zwischen den selektierten Gebieten festzustellen. Die Ergebnisse deuten auf einen grundlegenden Kenntnismangel in Europa hin, weshalb an dieser Stelle angeknüpft wird. Eine über Europa verteilte Stichprobe erlaubt eine neue morphologische Kategorisierung der großen Vielfalt an gefundenen physischen Formen. Die Menge an Gebieten erschließt sich in einer unbekannten Grundgesamtheit. Zur Datenaufbereitung bisheriger Analysen müssen Satellitenbilder manuell interpretiert werden. Das Verfahren birgt Unsicherheiten. Als kritische Selbstreflexion zeigt eine Reihe von Experimenten signifikante Unterschiede zwischen den Ergebnissen der Probanden auf, verdeutlicht jedoch bei ein und derselben Person Beständigkeit.
Wind energy is a key option in global dialogues about climate change mitigation. Here, we combined observations from surface wind stations, reanalysis datasets, and state‐of‐the‐art regional climate models from the Coordinated Regional Climate Downscaling Experiment (CORDEX Africa) to study the current and future wind energy potential in Zambia. We found that winds are dominated by southeasterlies and are rarely strong with an average speed of 2.8 m·s\(^{−1}\). When we converted the observed surface wind speed to a turbine hub height of 100 m, we found a ~38% increase in mean wind speed for the period 1981–2000. Further, both simulated and observed wind speed data show statistically significant increments across much of the country. The only areas that divert from this upward trend of wind speeds are the low land terrains of the Eastern Province bordering Malawi. Examining projections of wind power density (WPD), we found that although wind speed is increasing, it is still generally too weak to support large‐scale wind power generation. We found a meagre projected annual average WPD of 46.6 W·m\(^{−2}\). The highest WPDs of ~80 W·m\(^{−2}\) are projected in the northern and central parts of the country while the lowest are to be expected along the Luangwa valley in agreement with wind speed simulations. On average, Zambia is expected to experience minor WPD increments of 0.004 W·m\(^{−2}\) per year from 2031 to 2050. We conclude that small‐scale wind turbines that accommodate cut‐in wind speeds of 3.8 m·s\(^{−1}\) are the most suitable for power generation in Zambia. Further, given the limitations of small wind turbines, they are best suited for rural and suburban areas of the country where obstructions are few, thus making them ideal for complementing the government of the Republic of Zambia's rural electrification efforts.
Atmospheric circulation is a key driver of climate variability, and the representation of atmospheric circulation modes in regional climate models (RCMs) can enhance the credibility of regional climate projections. This study examines the representation of large‐scale atmospheric circulation modes in Coupled Model Inter‐comparison Project phase 5 RCMs once driven by ERA‐Interim, and by two general circulation models (GCMs). The study region is Western Europe and the circulation modes are classified using the Promax rotated T‐mode principal component analysis. The results indicate that the RCMs can replicate the classified atmospheric modes as obtained from ERA5 reanalysis, though with biases dependent on the data providing the lateral boundary condition and the choice of RCM. When the boundary condition is provided by ERA‐Interim that is more consistent with observations, the simulated map types and the associating time series match well with their counterparts from ERA5. Further, on average, the multi‐model ensemble mean of the analysed RCMs, driven by ERA‐Interim, indicated a slight improvement in the representation of the modes obtained from ERA5. Conversely, when the RCMs are driven by the GCMs that are models without assimilation of observational data, the representation of the atmospheric modes, as obtained from ERA5, is relatively less accurate compared to when the RCMs are driven by ERA‐Interim. This suggests that the biases stem from the GCMs. On average, the representation of the modes was not improved in the multi‐model ensemble mean of the five analysed RCMs driven by either of the GCMs. However, when the best‐performed RCMs were selected on average the ensemble mean indicated a slight improvement. Moreover, the presence of the North Atlantic Oscillation (NAO) in the simulated modes depends also on the lateral boundary conditions. The relationship between the modes and the NAO was replicated only when the RCMs were driven by reanalysis. The results indicate that the forcing model is the main factor in reproducing the atmospheric circulation.
In the past decades, various Earth observation-based time series products have emerged, which have enabled studies and analysis of global change processes. Besides their contribution to understanding past processes, time series datasets hold enormous potential for predictive modeling and thereby meet the demands of decision makers on future scenarios. In order to further exploit these data, a novel pixel-based approach has been introduced, which is the spatio-temporal matrix (STM). The approach integrates the historical characteristics of a specific land cover at a high temporal frequency in order to interpret the spatial and temporal information for the neighborhood of a given target pixel. The provided information can be exploited with common predictive models and algorithms. In this study, this approach was utilized and evaluated for the prediction of future urban/built-settlement growth. Random forest and multi-layer perceptron were employed for the prediction. The tests have been carried out with training strategies based on a one-year and a ten-year time span for the urban agglomerations of Surat (India), Ho-Chi-Minh City (Vietnam), and Abidjan (Ivory Coast). The slope, land use, exclusion, urban, transportation, hillshade (SLEUTH) model was selected as a baseline indicator for the performance evaluation. The statistical results from the receiver operating characteristic curve (ROC) demonstrate a good ability of the STM to facilitate the prediction of future settlement growth and its transferability to different cities, with area under the curve (AUC) values greater than 0.85. Compared with SLEUTH, the STM-based model achieved higher AUC in all of the test cases, while being independent of the additional datasets for the restricted and the preferential development areas.
Nationwide and consistent information on agricultural land use forms an important basis for sustainable land management maintaining food security, (agro)biodiversity, and soil fertility, especially as German agriculture has shown high vulnerability to climate change. Sentinel-1 and Sentinel-2 satellite data of the Copernicus program offer time series with temporal, spatial, radiometric, and spectral characteristics that have great potential for mapping and monitoring agricultural crops. This paper presents an approach which synergistically uses these multispectral and Synthetic Aperture Radar (SAR) time series for the classification of 17 crop classes at 10 m spatial resolution for Germany in the year 2018. Input data for the Random Forest (RF) classification are monthly statistics of Sentinel-1 and Sentinel-2 time series. This approach reduces the amount of input data and pre-processing steps while retaining phenological information, which is crucial for crop type discrimination. For training and validation, Land Parcel Identification System (LPIS) data were available covering 15 of the 16 German Federal States. An overall map accuracy of 75.5% was achieved, with class-specific F1-scores above 80% for winter wheat, maize, sugar beet, and rapeseed. By combining optical and SAR data, overall accuracies could be increased by 6% and 9%, respectively, compared to single sensor approaches. While no increase in overall accuracy could be achieved by stratifying the classification in natural landscape regions, the class-wise accuracies for all but the cereal classes could be improved, on average, by 7%. In comparison to census data, the crop areas could be approximated well with, on average, only 1% of deviation in class-specific acreages. Using this streamlined approach, similar accuracies for the most widespread crop types as well as for smaller permanent crop classes were reached as in other Germany-wide crop type studies, indicating its potential for repeated nationwide crop type mapping.
In Germany, as in many Western societies, demographic change will lead to a higher number of senior visitors to natural recreational areas and national parks. Given the high physiological requirements of many outdoor recreation activities, especially in mountain areas, it seems likely that demographic change will affect the spatial behaviour of national park visitors, which may pose a challenge to the management of these areas. With the help of GPS tracking and a standardized questionnaire (n=481), this study empirically investigates the spatial behaviour of demographic age brackets in Berchtesgaden National Park (NP) and the potential effects of demographic change on the use of the area. Cluster analysis revealed four activity types in the study area. More than half of the groups with visitors aged 60 and older belong to the activity type of Walker.
Inland surface water is often the most accessible freshwater source. As opposed to groundwater, surface water is replenished in a comparatively quick cycle, which makes this vital resource — if not overexploited — sustainable. From a global perspective, freshwater is plentiful. Still, depending on the region, surface water availability is severely limited. Additionally, climate change and human interventions act as large-scale drivers and cause dramatic changes in established surface water dynamics. Actions have to be taken to secure sustainable water availability and usage. This requires informed decision making based on reliable environmental data. Monitoring inland surface water dynamics is therefore more important than ever. Remote sensing is able to delineate surface water in a number of ways by using optical as well as active and passive microwave sensors. In this review, we look at the proceedings within this discipline by reviewing 233 scientific works. We provide an extensive overview of used sensors, the spatial and temporal resolution of studies, their thematic foci, and their spatial distribution. We observe that a wide array of available sensors and datasets, along with increasing computing capacities, have shaped the field over the last years. Multiple global analysis-ready products are available for investigating surface water area dynamics, but so far none offer high spatial and temporal resolution.
In most countries, freight is predominantly transported by road cargo trucks. We present a new satellite remote sensing method for detecting moving trucks on roads using Sentinel-2 data. The method exploits a temporal sensing offset of the Sentinel-2 multispectral instrument, causing spatially and spectrally distorted signatures of moving objects. A random forest classifier was trained (overall accuracy: 84%) on visual-near-infrared-spectra of 2500 globally labelled targets. Based on the classification, the target objects were extracted using a developed recursive neighbourhood search. The speed and the heading of the objects were approximated. Detections were validated by employing 350 globally labelled target boxes (mean F\(_1\) score: 0.74). The lowest F\(_1\) score was achieved in Kenya (0.36), the highest in Poland (0.88). Furthermore, validated at 26 traffic count stations in Germany on in sum 390 dates, the truck detections correlate spatio-temporally with station figures (Pearson r-value: 0.82, RMSE: 43.7). Absolute counts were underestimated on 81% of the dates. The detection performance may differ by season and road condition. Hence, the method is only suitable for approximating the relative truck traffic abundance rather than providing accurate absolute counts. However, existing road cargo monitoring methods that rely on traffic count stations or very high resolution remote sensing data have limited global availability. The proposed moving truck detection method could fill this gap, particularly where other information on road cargo traffic are sparse by employing globally and freely available Sentinel-2 data. It is inferior to the accuracy and the temporal detail of station counts, but superior in terms of spatial coverage.
Detection of grassland mowing events for Germany by combining Sentinel-1 and Sentinel-2 time series
(2022)
Grasslands cover one-third of the agricultural area in Germany and play an important economic role by providing fodder for livestock. In addition, they fulfill important ecosystem services, such as carbon storage, water purification, and the provision of habitats. These ecosystem services usually depend on the grassland management. In central Europe, grasslands are grazed and/or mown, whereby the management type and intensity vary in space and time. Spatial information on the mowing timing and frequency on larger scales are usually not available but would be required in order to assess the ecosystem services, species composition, and grassland yields. Time series of high-resolution satellite remote sensing data can be used to analyze the temporal and spatial dynamics of grasslands. Within this study, we aim to overcome the drawbacks identified by previous studies, such as optical data availability and the lack of comprehensive reference data, by testing the time series of various Sentinel-2 (S2) and Sentinal-1 (S1) parameters and combinations of them in order to detect mowing events in Germany in 2019. We developed a threshold-based algorithm by using information from a comprehensive reference dataset of heterogeneously managed grassland parcels in Germany, obtained by RGB cameras. The developed approach using the enhanced vegetation index (EVI) derived from S2 led to a successful mowing event detection in Germany (60.3% of mowing events detected, F1-Score = 0.64). However, events shortly before, during, or shortly after cloud gaps were missed and in regions with lower S2 orbit coverage fewer mowing events were detected. Therefore, S1-based backscatter, InSAR, and PolSAR features were investigated during S2 data gaps. From these, the PolSAR entropy detected mowing events most reliably. For a focus region, we tested an integrated approach by combining S2 and S1 parameters. This approach detected additional mowing events, but also led to many false positive events, resulting in a reduction in the F1-Score (from 0.65 of S2 to 0.61 of S2 + S1 for the focus region). According to our analysis, a majority of grasslands in Germany are only mown zero to two times (around 84%) and are probably additionally used for grazing. A small proportion is mown more often than four times (3%). Regions with a generally higher grassland mowing frequency are located in southern, south-eastern, and northern Germany.
The Northern Bald Ibis (Geronticus eremita, NBI) is an endangered migratory species, which went extinct in Europe in the 17th century. Currently, a translocation project in the frame of the European LIFE program is carried out, to reintroduce a migratory population with breeding colonies in the northern and southern Alpine foothills and a common wintering area in southern Tuscany. The population meanwhile consists of about 200 individuals, with about 90% of them carrying a GPS device on their back. We used biologging data from 2021 to model the habitat suitability for the species in the northern Alpine foothills. To set up a species distribution model, indices describing environmental conditions were calculated from satellite images of Landsat-8, and in addition to the well-proven use of optical remote sensing data, we also included Sentinel-1 actively sensed observation data, as well as climate and urbanization data. A random forest model was fitted on NBI GPS positions, which we used to identify regions with high predicted foraging suitability within the northern Alpine foothills. The model resulted in 84.5% overall accuracy. Elevation and slope had the highest predictive power, followed by grass cover and VV intensity of Sentinel-1 radar data. The map resulting from the model predicts the highest foraging suitability for valley floors, especially of Inn, Rhine, and Salzach-Valley as well as flatlands, like the Swiss Plateau and the agricultural areas surrounding Lake Constance. Areas with a high suitability index largely overlap with known historic breeding sites. This is particularly noteworthy because the model only refers to foraging habitats without considering the availability of suitable breeding cliffs. Detailed analyses identify the transition zone from extensive grassland management to intensive arable farming as the northern range limit. The modeling outcome allows for defining suitable areas for further translocation and management measures in the frame of the European NBI reintroduction program. Although required in the international IUCN translocation guidelines, the use of models in the context of translocation projects is still not common and in the case of the Northern Bald Ibis not considered in the present Single Species Action Plan of the African-Eurasian Migratory Water bird Agreement. Our species distribution model represents a contemporary snapshot, but sustainability is essential for conservation planning, especially in times of climate change. In this regard, a further model could be optimized by investigating sustainable land use, temporal dynamics, and climate change scenarios.
A first assessment of canopy cover loss in Germany's forests after the 2018–2020 drought years
(2022)
Central Europe was hit by several unusually strong periods of drought and heat between 2018 and 2020. These droughts affected forest ecosystems. Cascading effects with bark beetle infestations in spruce stands were fatal to vast forest areas in Germany. We present the first assessment of canopy cover loss in Germany for the period of January 2018–April 2021. Our approach makes use of dense Sentinel-2 and Landsat-8 time-series data. We computed the disturbance index (DI) from the tasseled cap components brightness, greenness, and wetness. Using quantiles, we generated monthly DI composites and calculated anomalies in a reference period (2017). From the resulting map, we calculated the canopy cover loss statistics for administrative entities. Our results show a canopy cover loss of 501,000 ha for Germany, with large regional differences. The losses were largest in central Germany and reached up to two-thirds of coniferous forest loss in some districts. Our map has high spatial (10 m) and temporal (monthly) resolution and can be updated at any time.
The analysis of the Earth system and interactions among its spheres is increasingly important to improve the understanding of global environmental change. In this regard, Earth observation (EO) is a valuable tool for monitoring of long term changes over the land surface and its features. Although investigations commonly study environmental change by means of a single EO-based land surface variable, a joint exploitation of multivariate land surface variables covering several spheres is still rarely performed. In this regard, we present a novel methodological framework for both, the automated processing of multisource time series to generate a unified multivariate feature space, as well as the application of statistical time series analysis techniques to quantify land surface change and driving variables. In particular, we unify multivariate time series over the last two decades including vegetation greenness, surface water area, snow cover area, and climatic, as well as hydrological variables. Furthermore, the statistical time series analyses include quantification of trends, changes in seasonality, and evaluation of drivers using the recently proposed causal discovery algorithm Peter and Clark Momentary Conditional Independence (PCMCI). We demonstrate the functionality of our methodological framework using Indo-Gangetic river basins in South Asia as a case study. The time series analyses reveal increasing trends in vegetation greenness being largely dependent on water availability, decreasing trends in snow cover area being mostly negatively coupled to temperature, and trends of surface water area to be spatially heterogeneous and linked to various driving variables. Overall, the obtained results highlight the value and suitability of this methodological framework with respect to global climate change research, enabling multivariate time series preparation, derivation of detailed information on significant trends and seasonality, as well as detection of causal links with minimal user intervention. This study is the first to use multivariate time series including several EO-based variables to analyze land surface dynamics over the last two decades using the causal discovery algorithm PCMCI.
Intercomparison of satellite-derived vegetation phenology is scarce in remote locations because of the limited coverage area and low temporal resolution of field observations. By their reliable near-ground observations and high-frequency data collection, PhenoCams can be a robust tool for intercomparison of land surface phenology derived from satellites. This study aims to investigate the transition dates of black spruce (Picea mariana (Mill.) B.S.P.) phenology by comparing fortnightly the MODIS normalized difference vegetation index (NDVI) and the enhanced vegetation index (EVI) extracted using the Google Earth Engine (GEE) platform with the daily PhenoCam-based green chromatic coordinate (GCC) index. Data were collected from 2016 to 2019 by PhenoCams installed in six mature stands along a latitudinal gradient of the boreal forests of Quebec, Canada. All time series were fitted by double-logistic functions, and the estimated parameters were compared between NDVI, EVI, and GCC. The onset of GCC occurred in the second week of May, whereas the ending of GCC occurred in the last week of September. We demonstrated that GCC was more correlated with EVI (R\(^2\) from 0.66 to 0.85) than NDVI (R\(^2\) from 0.52 to 0.68). In addition, the onset and ending of phenology were shown to differ by 3.5 and 5.4 days between EVI and GCC, respectively. Larger differences were detected between NDVI and GCC, 17.05 and 26.89 days for the onset and ending, respectively. EVI showed better estimations of the phenological dates than NDVI. This better performance is explained by the higher spectral sensitivity of EVI for multiple canopy leaf layers due to the presence of an additional blue band and an optimized soil factor value. Our study demonstrates that the phenological observations derived from PhenoCam are comparable with the EVI index. We conclude that EVI is more suitable than NDVI to assess phenology in evergreen species of the northern boreal region, where PhenoCam data are not available. The EVI index could be used as a reliable proxy of GCC for monitoring evergreen species phenology in areas with reduced access, or where repeated data collection from remote areas are logistically difficult due to the extreme weather.
Freely available satellite data at Google Earth Engine (GEE) cloud platform enables vegetation phenology analysis across different scales very efficiently. We evaluated seasonal and annual phenology of the old-growth Hyrcanian forests (HF) of northern Iran covering an area of ca. 1.9 million ha, and also focused on 15 UNESCO World Heritage Sites. We extracted bi-weekly MODIS-NDVI between 2017 and 2020 in GEE, which was used to identify the range of NDVI between two temporal stages. Then, changes in phenology and growth were analyzed by Sentinel 2-derived Temporal Normalized Phenology Index. We modelled between seasonal phenology and growth by additionally considering elevation, surface temperature, and monthly precipitation. Results indicated considerable difference in onset of forests along the longitudinal gradient of the HF. Faster growth was observed in low- and uplands of the western zone, whereas it was lower in both the mid-elevations and the western outskirts. Longitudinal range was a major driver of vegetation growth, to which environmental factors also differently but significantly contributed (p < 0.0001) along the west-east gradient. Our study developed at GEE provides a benchmark to examine the effects of environmental parameters on the vegetation growth of HF, which cover mountainous areas with partly no or limited accessibility.
Earth Observation satellite data allows for the monitoring of the surface of our planet at predefined intervals covering large areas. However, there is only one medium resolution sensor family in orbit that enables an observation time span of 40 and more years at a daily repeat interval. This is the AVHRR sensor family. If we want to investigate the long-term impacts of climate change on our environment, we can only do so based on data that remains available for several decades. If we then want to investigate processes with respect to climate change, we need very high temporal resolution enabling the generation of long-term time series and the derivation of related statistical parameters such as mean, variability, anomalies, and trends. The challenges to generating a well calibrated and harmonized 40-year-long time series based on AVHRR sensor data flown on 14 different platforms are enormous. However, only extremely thorough pre-processing and harmonization ensures that trends found in the data are real trends and not sensor-related (or other) artefacts. The generation of European-wide time series as a basis for the derivation of a multitude of parameters is therefore an extremely challenging task, the details of which are presented in this paper.
Land Surface Temperature (LST) is an important parameter for tracing the impact of changing climatic conditions on our environment. Describing the interface between long- and shortwave radiation fluxes, as well as between turbulent heat fluxes and the ground heat flux, LST plays a crucial role in the global heat balance. Satellite-derived LST is an indispensable tool for monitoring these changes consistently over large areas and for long time periods. Data from the AVHRR (Advanced Very High-Resolution Radiometer) sensors have been available since the early 1980s. In the TIMELINE project, LST is derived for the entire operating period of AVHRR sensors over Europe at a 1 km spatial resolution. In this study, we present the validation results for the TIMELINE AVHRR daytime LST. The validation approach consists of an assessment of the temporal consistency of the AVHRR LST time series, an inter-comparison between AVHRR LST and in situ LST, and a comparison of the AVHRR LST product with concurrent MODIS (Moderate Resolution Imaging Spectroradiometer) LST. The results indicate the successful derivation of stable LST time series from multi-decadal AVHRR data. The validation results were investigated regarding different LST, TCWV and VA, as well as land cover classes. The comparisons between the TIMELINE LST product and the reference datasets show seasonal and land cover-related patterns. The LST level was found to be the most determinative factor of the error. On average, an absolute deviation of the AVHRR LST by 1.83 K from in situ LST, as well as a difference of 2.34 K from the MODIS product, was observed.
Forest conservation is of particular concern in tropical regions where a large refuge of biodiversity is still existing. These areas are threatened by deforestation, forest degradation and fragmentation. Especially, pressures of anthropogenic activities adjacent to these areas significantly influence conservation effectiveness. Ecuador was chosen as study area since it is a globally relevant center of forest ecosystems and biodiversity. We identified hotspots of deforestation on the national level of continental Ecuador between 1990 and 2018, analyzed the most significant drivers of deforestation on national and biome level (the Coast, the Andes, The Amazon) as well as inside protected areas in Ecuador by using multiple regression analysis. We separated the national system of protected areas (SNAP) into higher and lower protection levels. Besides SNAP, we also considered Biosphere Reserves (BRs) and Ramsar sites. In addition, we investigated the rates and spatial patterns of deforestation in protected areas and buffer zones (5 km and 10 km outwards the protected area boundaries) using landscape metrics. Between 1990 and 2018, approximately 4% of the accumulated deforestation occurred within the boundaries of SNAP, and up to 25.5% in buffer zones. The highest rates of deforestation have been found in the 5 km buffer zone around the protected areas with the highest protection level. Protected areas and their buffer zones with higher protection status were identified as the most deforested areas among SNAP. BRs had the highest deforestation rates among all protected areas but most of these areas just became BRs after the year 2000. The most important driver of deforestation is agriculture. Other relevant drivers differ between the biomes. The results suggest that the SNAP is generally effective to prevent deforestation within their protection boundaries. However, deforestation around protected areas can undermine conservation strategies to sustain biodiversity. Actions to address such dynamics and patterns of deforestation and forest fragmentation, and developing conservation strategies of their landscape context are urgently needed especially in the buffer zones of areas with the highest protection status.
Coal mining, an important human activity, disturbs soil organic carbon (SOC) accumulation and decomposition, eventually affecting terrestrial carbon cycling and the sustainability of human society. However, changes of SOC content and their relation with influential factors in coal mining areas remained unclear. In the study, predictive models of SOC content were developed based on field sampling and Landsat images for different land-use types (grassland, forest, farmland, and bare land) of the largest coal mining area in China (i.e., Shendong). The established models were employed to estimate SOC content across the Shendong mining area during 1990–2020, followed by an investigation into the impacts of climate change and human disturbance on SOC content by a Geo-detector. Results showed that the models produced satisfactory results (R\(^2\) > 0.69, p < 0.05), demonstrating that SOC content over a large coal mining area can be effectively assessed using remote sensing techniques. Results revealed that average SOC content in the study area rose from 5.67 gC·kg\(^{−1}\) in 1990 to 9.23 gC·kg\(^{−1}\) in 2010 and then declined to 5.31 gC·Kg\(^{−1}\) in 2020. This could be attributed to the interaction between the disturbance of soil caused by coal mining and the improvement of eco-environment by land reclamation. Spatially, the SOC content of farmland was the highest, followed by grassland, and that of bare land was the lowest. SOC accumulation was inhibited by coal mining activities, with the effect of high-intensity mining being lower than that of moderate- and low-intensity mining activities. Land use was found to be the strongest individual influencing factor for SOC content changes, while the interaction between vegetation coverage and precipitation exerted the most significant influence on the variability of SOC content. Furthermore, the influence of mining intensity combined with precipitation was 10 times higher than that of mining intensity alone.
Regional climate models (RCMs) are tools used to project future climate change at a regional scale. Despite their high horizontal resolution, RCMs are characterized by systematic biases relative to observations, which can result in unrealistic interpretations of future climate change signals. On the other hand, bias correction (BC) is a popular statistical post-processing technique applied to improve the usability of output from climate models. Like every other statistical technique, BC has its strengths and weaknesses. Hence, within the regional context of Germany, and for temperature and precipitation, this study is dedicated to the assessment of the impact of different BC techniques on the RCM output. The focuses are on the impact of BC on the RCM’s statistical characterization, and physical consistency defined as the spatiotemporal consistency between the bias-corrected variable and the simulated physical mechanisms governing the variable, as well as the correlations between the bias-corrected variable and other (simulated) climate variables. Five BC techniques were applied in adjusting the systematic biases in temperature and precipitation RCM outputs. The BC techniques are linear scaling, empirical quantile mapping, univariate quantile delta mapping, multivariate quantile delta mapping that considers inter-site dependencies, and multivariate quantile delta mapping that considers inter-variable dependencies (MBCn). The results show that each BC technique adds value in reducing the biases in the statistics of the RCM output, though the added value depends on several factors such as the temporal resolution of the data, choice of RCM, climate variable, region, and the metric used in evaluating the BC technique. Further, the raw RCMs reproduced portions of the observed modes of atmospheric circulation in Western Europe, and the observed temperature, and precipitation meteorological patterns in Germany. After the BC, generally, the spatiotemporal configurations of the simulated meteorological patterns as well as the governing large-scale mechanisms were reproduced.
However, at a more localized spatial scale for the individual meteorological patterns, the BC changed the simulated co-variability of some grids, especially for precipitation. Concerning the co-variability among the variables, a physically interpretable positive correlation was found between temperature and precipitation during boreal winter in both models and observations. For most grid boxes in the study domain and on average, the BC techniques that do not adjust inter-variable dependency did not notably change the simulated correlations between the climate variables. However, depending on the grid box, the (univariate) BC techniques tend to degrade the simulated temporal correlations between temperature and precipitation. Further, MBCn which adjusts biases in inter-variable dependency has the skill to improve the correlations between the simulated variables towards observations.
Vietnam's 3260 km coastline is densely populated, experiences rapid urban and economic growth, and faces at the same time a high risk of coastal hazards. Satellite archives provide a free and powerful opportunity for long-term area-wide monitoring of the coastal zone. This paper presents an automated analysis of coastline dynamics from 1986 to 2021 for Vietnam's entire coastal zone using the Landsat archive. The proposed method is implemented within the cloud-computing platform Google Earth Engine to only involve publicly and globally available datasets and tools. We generated annual coastline composites representing the mean-high water level and extracted sub-pixel coastlines. We further quantified coastline change rates along shore-perpendicular transects, revealing that half of Vietnam's coast did not experience significant change, while the remaining half is classified as erosional (27.7%) and accretional (27.1%). A hotspot analysis shows that coastal segments with the highest change rates are concentrated in the low-lying deltas of the Mekong River in the south and the Red River in the north. Hotspots with the highest accretion rates of up to +47 m/year are mainly associated with the construction of artificial coastlines, while hotspots with the highest erosion rates of −28 m/year may be related to natural sediment redistribution and human activity.
An increasing amount of Brazilian rainforest is being lost or degraded for various reasons, both anthropogenic and natural, leading to a loss of biodiversity and further global consequences. Especially in the Brazilian state of Mato Grosso, soy production and large-scale cattle farms led to extensive losses of rainforest in recent years. We used a spectral mixture approach followed by a decision tree classification based on more than 30 years of Landsat data to quantify these losses. Research has shown that current methods for assessing forest degradation are lacking accuracy. Therefore, we generated classifications to determine land cover changes for each year, focusing on both cleared and degraded forest land. The analyses showed a decrease in forest area in Mato Grosso by 28.8% between 1986 and 2020. In order to measure changed forest structures for the selected period, fragmentation analyses based on diverse landscape metrics were carried out for the municipality of Colniza in Mato Grosso. It was found that forest areas experienced also a high degree of fragmentation over the study period, with an increase of 83.3% of the number of patches and a decrease of the mean patch area of 86.1% for the selected time period, resulting in altered habitats for flora and fauna.
Mapping of lava flows in unvegetated areas of active volcanoes using optical satellite data is challenging due to spectral similarities of volcanic deposits and the surrounding background. Using very high-resolution PlanetScope data, this study introduces a novel object-oriented classification approach for mapping lava flows in both vegetated and unvegetated areas during several eruptive phases of three Indonesian volcanoes (Karangetang 2018/2019, Agung 2017, Krakatau 2018/2019). For this, change detection analysis based on PlanetScope imagery for mapping loss of vegetation due to volcanic activity (e.g., lava flows) is combined with the analysis of changes in texture and brightness, with hydrological runoff modelling and with analysis of thermal anomalies derived from Sentinel-2 or Landsat-8. Qualitative comparison of the mapped lava flows showed good agreement with multispectral false color time series (Sentinel-2 and Landsat-8). Reports of the Global Volcanism Program support the findings, indicating the developed lava mapping approach produces valuable results for monitoring volcanic hazards. Despite the lack of bands in infrared wavelengths, PlanetScope proves beneficial for the assessment of risk and near-real-time monitoring of active volcanoes due to its high spatial (3 m) and temporal resolution (mapping of all subaerial volcanoes on a daily basis).
This article presents an open space concept of areas that are kept permanently free from buildings, technical infrastructure, and soil sealing. In the European Alps, space is scarce because of the topography; conflicts often arise between competing land uses such as permanent settlements and commercial activity. However, the presence of open spaces is important for carbon sequestration and the prevention of natural hazards, especially given climate change. A GIS-based analysis was conducted to identify an alpine-wide inventory of large-scale near-natural areas, or simply stated, open spaces. The method used identified the degree of infrastructure development for natural landscape units. Within the Alpine Convention perimeter, near-natural areas (with a degree of infrastructural development of up to 20%) account for a share of 51.5%. Only 14.5% of those areas are highly protected and are mostly located in high altitudes of over 1500 m or 2000 m above sea level. We advocate that the remaining Alpine open spaces must be preserved through the delimitation of more effective protection mechanisms, and green corridors should be safeguarded through spatial planning. To enhance the ecological connectivity of open spaces, there is the need for tailored spatial and sectoral planning strategies to prevent further landscape fragmentation and to coordinate new forms of land use for renewable energy production.
Comparing PlanetScope and Sentinel-2 imagery for mapping mountain pines in the Sarntal Alps, Italy
(2022)
The mountain pine (Pinus mugo ssp. Mugo Turra) is an important component of the alpine treeline ecotone and fulfills numerous ecosystem functions. To understand and quantify the impacts of increasing logging activities and climatic changes in the European Alps, accurate information on the occurrence and distribution of mountain pine stands is needed. While Earth observation provides up-to-date information on land cover, space-borne mapping of mountain pines is challenging as different coniferous species are spectrally similar, and small-structured patches may remain undetected due to the sensor’s spatial resolution. This study uses multi-temporal optical imagery from PlanetScope (3 m) and Sentinel-2 (10 m) and combines them with additional features (e.g., textural statistics (homogeneity, contrast, entropy, spatial mean and spatial variance) from gray level co-occurrence matrix (GLCM), topographic features (elevation, slope and aspect) and canopy height information) to overcome the present challenges in mapping mountain pine stands. Specifically, we assessed the influence of spatial resolution and feature space composition including the GLCM window size for textural features. The study site is covering the Sarntal Alps, Italy, a region known for large stands of mountain pine. Our results show that mountain pines can be accurately mapped (PlanetScope (90.96%) and Sentinel-2 (90.65%)) by combining all features. In general, Sentinel-2 can achieve comparable results to PlanetScope independent of the feature set composition, despite the lower spatial resolution. In particular, the inclusion of textural features improved the accuracy by +8% (PlanetScope) and +3% (Sentinel-2), whereas accuracy improvements of topographic features and canopy height were low. The derived map of mountain pines in the Sarntal Alps supports local forest management to monitor and assess recent and ongoing anthropogenic and climatic changes at the treeline. Furthermore, our study highlights the importance of freely available Sentinel-2 data and image-derived textural features to accurately map mountain pines in Alpine environments.
Die Bodenfeuchte stellt eine essenzielle Variable für den Energie-, Feuchte- und Stoffaustausch zwischen Landoberfläche und Atmosphäre dar. Ihre Auswirkungen auf Temperatur und Niederschlag sind vielfältig und komplex. Die in Klimamodellen verwendeten Schemata zur Simulation der Bodenfeuchte, auch bodenhydrologische Schemata genannt, sind aufgrund des Ursprungs der Klimamodelle aus Wettermodellen jedoch häufig sehr stark vereinfacht dargestellt.
Bei Klimamodellen, die Simulationen mit einer groben Auflösung von mehreren Zehner- oder Hunderterkilometern rechnen, können viele Prozesse vernachlässigt werden. Da die Auflösung der Klimamodelle jedoch stetig steigt und mittlerweile beim koordinierten Projekt regionaler Klimamodelle CORDEX-CORE standardmäßig bei 0.22° Kantenlänge liegt, müssen auch höher aufgelöste Daten und mehr Prozesse simuliert werden. Dies gilt erst recht mit Blick auf konvektionsauflösende Simulationen mit wenigen Kilometern Kantenlänge. Mit steigenden Modellauflösungen steigt zugleich die Komplexität und Differenziertheit der Fragestellungen, die mit Hilfe von Klimamodellen beantwortet werden sollen. An diesem Punkt setzt auch das Projekt BigData@Geo an, in dessen Rahmen die vorliegende Arbeit entstand. Ziel dieses Projektes ist es, hochaufgelöste Klimainformationen für den bayerischen Regierungsbezirk Unterfranken für Akteure aus der Land- und Forstwirtschaft sowie dem Weinbau zur Verfügung zu stellen.
Auf diesen angewandten und grundlegenden Anforderungen und Zielsetzungen basierend, bedarf auch das in dieser Arbeit verwendete regionale Klimamodell REMO (Version 2015) der weiteren Entwicklung. So ist das Hauptziel der Arbeit das bestehende einschichtige bodenhydrologische Schema durch ein mehrschichtiges zu ersetzen. Der Vorteil mehrerer simulierter Bodenschichten besteht darin, dass nun die vertikale Bewegung des Wassers in Form von Versickerung und kapillarem Aufstieg simuliert werden kann. Dies geschieht auf der Basis bodenhydrologischer Parameter, deren Wert in Abhängigkeit vom Boden und der Bodenfeuchte über die Wasserrückhaltekurve bestimmt wird. Für diese Kurve existieren verschiedene Parametrisierungen, von denen die Ansätze von Clapp-Hornberger und van Genuchten verwendet wurden. Außerdem kann die Bodenfeuchte nun bis zu einer Tiefe von circa 10 m beziehungsweise der Tiefe des anstehenden Gesteins simuliert werden. Damit besteht im Gegensatz zum vorherigen Schema, dessen Tiefe auf die Wurzeltiefe beschränkt ist, die Möglichkeit, dass Wasser auch unterhalb der Wurzeln zur Verfügung stehen kann und somit die absolute im Boden verfügbare Wassermenge zunimmt. Die Schichtung erlaubt darüber hinaus die Verdunstung aus unbewachsenem Boden lediglich auf Basis des in der obersten Schicht verfügbaren Wassers. Ein weiterer Prozess, der dank der Schichtung und der weiter unten erläuterten Datensätze neu parametrisiert werden kann, ist die Infiltration.
Für die Verwendung des Schemas sind Informationen über bodenhydrologische Parameter, die Wurzeltiefe und die Tiefe bis zum anstehenden Gestein erforderlich. Entsprechende Datensätze müssen hierfür aufbereitet und in das Modell eingebaut werden. Bezüglich der Wurzeltiefe wurden drei sich bezüglich der Tiefe, der Definition und der verfügbaren Auflösung stark voneinander unterscheidende Datensätze verglichen. Letztendlich wird die Wurzeltiefe aus dem mit einer anderen REMO-Version gekoppelten Vegetationsmodul iMOVE verwendet, da zukünftig eine Kopplung dieses Moduls mit dem mehrschichtigen Boden geplant ist und die Wurzeltiefen damit konsistent sind. Zudem ist die zugrundeliegende Auflösung der Daten hoch und es werden maximale Wurzeltiefen berücksichtigt, die besonders wichtig für die Simulation von Landoberfläche-Atmosphäre-Interaktionen sind. Diese Vorteile brachten die anderen Datensätze nicht mit. In der finalen Modellversion werden für die Tiefe bis zum anstehenden Gestein und die Korngrößenverteilungen die Daten von SoilGrids verwendet. Ein Vergleich mit anderen Bodendatensätzen fand in einer parallel laufenden Dissertation statt (Ziegler 2022). Bei SoilGrids ist hervorzuheben, dass die Korngrößenverteilungen in einer hohen räumlichen Auflösung (1 km^2 oder höher) und mit mehreren vertikalen Schichten vorliegen. Gegenüber dem ursprünglich in REMO verwendeten Datensatz mit einer Kantenlänge von 0.5° und ohne vertikale Differenzierung ist dies eine starke Verbesserung der Eingangsdaten. Dazu kommt, dass die Korngrößenverteilungen die Verwendung kontinuierlicher Pedotransferfunktionen statt fünf diskreter Texturklassen, denen für die bodenhydrologischen Parameter fixe Tabellenwerte zugewiesen werden, ermöglichen. Dies führt zu einer deutlich besseren Differenzierung des heterogenen Bodens.
Im Rahmen der Arbeit wurden insgesamt 19 Simulationen für Europa und ein erweitertes Deutschlandgebiet mit Auflösungen von 0.44° beziehungsweise 0.11° für den Zeitraum 2000 bis 2018 gerechnet. Dabei zeigte sich, dass die Einführung des mehrschichtigen Bodenschemas gegenüber dem einschichtigen Schema zu einer Verringerung der Bodenfeuchte in der Wurzeltiefe führt. Nichtsdestotrotz nimmt die absolute Wassermenge des Bodens durch die Berücksichtigung des Bodens unterhalb der Wurzelzone zu. Bezogen auf die einzelnen Schichten wird die Bodenfeuchte damit zwar unterschätzt, im Laufe der Modellentwicklung kann jedoch eine Verbesserung im Vergleich zu ERA5 erzielt werden. Das neue Schema führt zu einer Verringerung der Evapotranspiration, die über alle Schritte der Modellentwicklung und besonders während der Sommermonate auftritt. Im Vergleich zu Validationsdaten von ERA5 und GLEAM zeigt sich, dass dies eine Verbesserung dieser Größe bedeutet, die sowohl in der Fläche als auch beim Fehler und in der Verteilung auftritt.
Gleiches lässt sich für den Oberflächenabfluss sagen. Hierfür implementierte Schemata (Philip, Green-Ampt), die anders als das standardmäßig verwendete Improved-Arno-Schema bodenhydrologische Parameter berücksichtigen, konnten eine weitere Verbesserung im Flachland zeigen. In Gebirgsregionen nahm der Fehler durch die nicht enthaltene Berücksichtigung der Hangneigung jedoch zu, sodass in der finalen Modellversion auf das Improved-Arno-Schema zurückgegriffen wurde. Die Temperatur steigt durch die ursprüngliche Version des mehrschichtigen Schemas zunächst an, was zu einer Über- statt der vorherigen Unterschätzung gegenüber E-OBS führt. Die Modellentwicklung resultiert zwar in einer Reduzierung der Temperatur, jedoch fällt diese zu stark aus, sodass der Temperaturfehler letztendlich größer als in der einschichtigen Modellversion ist. Da die Evapotranspiration jedoch maßgeblich verbessert wurde, kann dieser Fehler eventuell auf ein übermäßiges Tuning der Temperatur zurückgeführt werden.
Die Betrachtung von Hitzeereignissen am Beispiel der Sommer 2003 und 2018 hat gezeigt, dass die Modellentwicklung dazu beiträgt, diese Ereignisse besser als das einschichtige Schema zu simulieren. Zwar trifft dies nicht auf das räumliche Verhalten der mittleren Temperatur zu, jedoch auf deren zeitlichen Verlauf. Hinzu kommt die bessere Simulation der täglichen Extrem- und besonders der Minimaltemperatur, was zu einer Erhöhung der täglichen Temperaturspanne führt. Diese wird von Klimamodellen in der Regel zu stark unterschätzt.
Durch die Berücksichtigung der vertikalen Wasserflüsse hat sich jedoch auch gezeigt, dass noch enormes Entwicklungspotenzial mit Blick auf (boden)hydrologische Prozesse besteht. Dies gilt in besonderem Maße für zukünftige Simulationen mit konvektionserlaubender Auflösung. So sollten subskalige Informationen des Bodens und der Orographie berücksichtigt werden. Dies dient einerseits der Repräsentation vorliegender Heterogenitäten und kann andererseits, wie am Beispiel der Infiltrationsschemata dargelegt, zur Verbesserung bestehender Prozesse beitragen. Da die simulierte Drainage durch das mehrschichtige Bodenschema im gleichen Maße zu- wie der Oberflächenabfluss abnimmt und das Wasser dem Modell in der Folge nicht weiter zur Verfügung steht, sollte zukünftig auch Grundwasser im Modell berücksichtigt werden. Eine Vielzahl von Studien konnte einen Mehrwert durch die Implementierung dieser Variable und damit verbundener Prozesse feststellen. Mittelfristig ist jedoch insgesamt die Kopplung an ein hydrologisches Modell zu empfehlen, um die bei hochauflösenden Simulationen relevanten Prozesse angemessen repräsentieren zu können. Hierfür bieten sich beispielsweise ParFlow oder mHM an.
Insgesamt ist festzuhalten, dass das mehrschichtige Bodenschema einen Mehrwert liefert, da schwer zu simulierende und in der Postprozessierung zu korrigierende Variablen wie die Evapotranspiration und der Oberflächenabfluss deutlich besser modelliert werden können als mit dem einschichtigen Schema. Dies gilt auch für die Extremtemperaturen. Beides ist klar auf die Schichtung des Bodens und damit einhergehender Prozesse zurückzuführen. Bezüglich der Daten zeigt sich, dass die Wurzeltiefe, die Berücksichtigung von SoilGrids und die vertikale Bodeninformation für die weitere Optimierung verantwortlich sind. Darüber hinaus ist der höhere Informationsgehalt, der anhand der geschichteten Bodenfeuchte zur Verfügung steht, ebenfalls als Mehrwert einzustufen.
Arctic permafrost coasts become increasingly vulnerable due to environmental drivers such as the reduced sea-ice extent and duration as well as the thawing of permafrost itself. A continuous quantification of the erosion process on large to circum-Arctic scales is required to fully assess the extent and understand the consequences of eroding permafrost coastlines. This study presents a novel approach to quantify annual Arctic coastal erosion and build-up rates based on Sentinel-1 (S1) Synthetic Aperture RADAR (SAR) backscatter data, in combination with Deep Learning (DL) and Change Vector Analysis (CVA). The methodology includes the generation of a high-quality Arctic coastline product via DL, which acted as a reference for quantifying coastal erosion and build-up rates from annual median and standard deviation (sd) backscatter images via CVA. The analysis was applied on ten test sites distributed across the Arctic and covering about 1038 km of coastline. Results revealed maximum erosion rates of up to 160 m for some areas and an average erosion rate of 4.37 m across all test sites within a three-year temporal window from 2017 to 2020. The observed erosion rates within the framework of this study agree with findings published in the previous literature. The proposed methods and data can be applied on large scales and, prospectively, even for the entire Arctic. The generated products may be used for quantifying the loss of frozen ground, estimating the release of stored organic material, and can act as a basis for further related studies in Arctic coastal environments.
While the place of birth plays a crucial role for women’s birth experiences, the interest in out-of-hospital births has increased during the Covid-19 pandemic. Related to this, various international policies recommend enabling women to choose where to give birth. We aimed to analyze Swiss women’s choice between birth hospitals and birth centers. Employing spatial accessibility analysis, we incorporated four data types: highly disaggregated population data, administrative data, street network data, addresses of birth hospitals and birth centers. 99.8% of Swiss women of childbearing age were included in the analysis (N = 1.896.669). We modelled car travel times from a woman’s residence to the nearest birth hospital and birth center. If both birth settings were available within 30 minutes, a woman was considered to have a true choice. Only 58.2% of women had a true choice. This proportion varied considerably across Swiss federal states. The main barrier to a true choice was limited accessibility of birth centers. Median travel time to birth hospitals was 9.8 (M = 12.5), to birth centers 23.9 minutes (M = 28.5). Swiss women are insufficiently empowered to exercise their reproductive autonomy as their choice of place of birth is significantly limited by geographical constraints. It is an ethical and medical imperative to provide women with a true choice. We provide high-resolution insights into the accessibility of birth settings and strong arguments to (re-)examine the need for further birth centers (and birth hospitals) in specific geographical areas. Policy-makers are obligated to improve the accessibility of birth centers to advance women’s autonomy and enhance maternal health outcomes after childbirth. The Covid-19 pandemic offers an opportunity to shift policy.
The command area of the Rakh branch canal grows wheat, sugarcane, and rice crops in abundance. The canal water, which is trivial for irrigating these crops, is conveyed to the farms through the network of canals and distributaries. For the maintenance of this vast infrastructure; the end users are charged on a seasonal basis. The present water charges are severely criticized for not being adequate to properly manage the entire infrastructure. We use the residual value to determine the value of the irrigation water and then based on the quantity of irrigation water supplied to farm land coupled with the infrastructure maintenance cost, full cost recovery figures are executed for the study area, and policy recommendations are made for the implementation of the full cost recovery system. The approach is unique in the sense that the pricings are based on the actual quantity of water conveyed to the field for irrigating crops. The results of our analysis showed that the canal water is severely under charged in the culturable command area of selected distributaries, thus negating the plan of having a self-sustainable irrigation system.
This study investigates the surroundings of Munigua (municipium Flavium Muniguense), a small Roman town in the ancient province of Hispania Baetica (SW Spain). The city's economy was based primarily on copper and iron mining, which brought financial prosperity to its citizens. Local production of agricultural goods is thought to have been of little importance, as the regional soil conditions do not seem to be suitable for extensive agriculture.
To evaluate the recent soil agro-potential and to find evidence for prehistoric and historic land use in the surroundings of Munigua, we applied a pedo-geomorphological approach based on the physico-chemical analysis of 14 representative soil and sediment exposures. Selected samples were analyzed for bulk chemistry, texture and phytoliths. The chronostratigraphy of the sequences was based on radiocarbon dating of charcoal samples. The site evaluation of the present-day soil agro-potential was carried out according to standard procedures and included evaluation of potential rootability, available water-storage capacity and nutrient budget within the uppermost 1 m.
The results show that moderate to very good soil agro-potential prevails in the granitic and floodplain areas surrounding Munigua. Clearly, recent soil agro-potential in these areas allows the production of basic agricultural goods, and similar limited agricultural use should also have been possible in ancient times. In contrast, weak to very weak present-day soil agro-potential prevails in the metamorphic landscape due to the occurrence of shallow and sandy to stony soils.
In addition, the study provides pedo-geomorphological evidence for prehistoric and historic land use in pre-Roman, Roman and post-Roman times. Catenary soil mapping in the vicinity of a Roman house complex reveals multi-layered colluvial deposits. They document phases of hillslope erosion mainly triggered by human land use between 4063 ± 82 and 3796 ± 76 cal BP, around 2601 ± 115 cal BP, and between 1424 ± 96 and 421 ± 88 cal BP. Moreover, geochemical and phytolith analyses of a Roman hortic Anthrosol indicate the local cultivation of agricultural products that contributed to the food supply of Munigua.
Overall, the evidence of Roman agricultural use in the Munigua area indicates that the city's economy was by no means focused solely on mining. The production of basic agricultural products was also part of Munigua's economic portfolio. Our geoarcheological study thus supports the archeological concept of economically diversified Roman cities in the province of Baetica and in Hispania.
The Seville Strategy spurred a signifi cant paradigm shift in UNESCO’s MAB Programme, re-conceptualising the research programme as a modern tool for the dual mandate of nature conservation and sustainable development. However, many biosphere reserves failed to comply with the new regulations and in 2013 the ‘Exit Strategy’ was announced to improve the quality of the global network.
This study presents a global assessment of the implementation of the quality enhancement strategies, highlighting signifi cant differences worldwide through 20 country-specifi c case studies. It concludes that the strategies have been fundamental in improving the credibility and coherence of the MAB Programme. Challenges in the implementation were not unique to individual countries but were common to all Member States with pre-Seville sites, and in many states the process has led to a rejuvenation of national biosphere reserve networks.
Der Lebensmittelonlinehandel in Deutschland gewann, verstärkt durch die Covid-19-Pandemie, an Umsatzanteilen im Lebensmitteleinzelhandel. Hierdurch wurden neue Anforderungen an Arbeit und Beschäftigung in Deutschland geschaffen. Insbesondere in urbanen Räumen hat die Lebensmittelzustellung durch neu entstandene Betriebsformen zugenommen. So entstehen durch das Versprechen der Betriebe, Lebensmittel in kurzen Zeiträumen zu liefern, verschiedene Logistikstandorte und u.a. urbane Fahrradlieferdienste. Während Medien und Gewerkschaften bereits vor der Entstehung prekärer Arbeitsbedingungen warnen, sind die genauen Auswirkungen des Lebensmittelonlinehandels auf die Entwicklung neuer Arbeitsstandorte und die dort stattfindende Beschäftigung nur unzureichend bekannt. Diese Arbeit untersucht den Lebensmittelonlinehandel anhand seiner Betriebsformen, Standorte und Arbeitsprozesse sowie deren Auswirkungen auf Beschäftigte in Deutschland. Den konzeptionellen Hintergrund bilden Arbeiten der geographischen Handelsforschung sowie Debatten zu Arbeitsplatzqualität und Beschäftigung. Für die Analyse sind Primärdaten und Sekundärdaten erhoben worden. Es zeigt sich, dass teilweise komplexe Betriebsformen entstehen, bei denen sich die Arbeit und Arbeitsorte verändern. Zudem entstehen neue Herausforderungen für die Beschäftigten (u.a. physische und psychische Belastung), welche in dieser Arbeit identifiziert werden.
Snow cover (SC) and timing of snowmelt are key regulators of a wide range of Arctic ecosystem functions. Both are strongly influenced by the amplified Arctic warming and essential variables to understand environmental changes and their dynamics. This study evaluates the potential of Sentinel-1 (S-1) synthetic aperture radar (SAR) time series for monitoring SC depletion and snowmelt with high spatiotemporal resolution to capture their understudied small-scale heterogeneity. We use 97 dual-polarized S-1 SAR images acquired over northeastern Greenland and 94 over southwestern Greenland in the interferometric wide swath mode from the years 2017 and 2018. Comparison of S-1 intensity against SC fraction maps derived from orthorectified terrestrial time-lapse imagery indicates that SAR backscatter can increase before a decrease in SC fraction is observed. Hence, the increase in backscatter is related to changing snowpack properties during the runoff phase as well as decreasing SC fraction. We here present a novel empirical approach based on the temporal evolution of the SAR signal to identify start of runoff (SOR), end of snow cover (EOS) and SC extent for each S-1 observation date during melt using backscatter thresholds as well as the derivative. Comparison of SC with orthorectified time-lapse imagery indicates that HV polarization outperforms HH when using a global threshold. The derivative avoids manual selection of thresholds and adapts to different environmental settings and seasonal conditions. With a global configuration (threshold: 4 dB; polarization: HV) as well as with the derivative, the overall accuracy of SC maps was in all cases above 75 % and in more than half of cases above 90 %. Based on the physical principle of SAR backscatter during snowmelt, our approach is expected to work well in other low-vegetation areas and, hence, could support large-scale SC monitoring at high spatiotemporal resolution (20 m, 6 d) with high accuracy.
Central Europe experienced several droughts in the recent past, such as in the year 2018, which was characterized by extremely low rainfall rates and high temperatures, resulting in substantial agricultural yield losses. Time series of satellite earth observation data enable the characterization of past drought events over large temporal and spatial scales. Within this study, Moderate Resolution Spectroradiometer (MODIS) Enhanced Vegetation Index (EVI) (MOD13Q1) 250 m time series were investigated for the vegetation periods of 2000 to 2018. The spatial and temporal development of vegetation in 2018 was compared to other dry and hot years in Europe, like the drought year 2003. Temporal and spatial inter- and intra-annual patterns of EVI anomalies were analyzed for all of Germany and for its cropland, forest, and grassland areas individually. While vegetation development in spring 2018 was above average, the summer months of 2018 showed negative anomalies in a similar magnitude as in 2003, which was particularly apparent within grassland and cropland areas in Germany. In contrast, the year 2003 showed negative anomalies during the entire growing season. The spatial pattern of vegetation status in 2018 showed high regional variation, with north-eastern Germany mainly affected in June, north-western parts in July, and western Germany in August. The temporal pattern of satellite-derived EVI deviances within the study period 2000-2018 were in good agreement with crop yield statistics for Germany. The study shows that the EVI deviation of the summer months of 2018 were among the most extreme in the study period compared to other years. The spatial pattern and temporal development of vegetation condition between the drought years differ.
The monitoring of species and functional diversity is of increasing relevance for the development of strategies for the conservation and management of biodiversity. Therefore, reliable estimates of the performance of monitoring techniques across taxa become important. Using a unique dataset, this study investigates the potential of airborne LiDAR-derived variables characterizing vegetation structure as predictors for animal species richness at the southern slopes of Mount Kilimanjaro. To disentangle the structural LiDAR information from co-factors related to elevational vegetation zones, LiDAR-based models were compared to the predictive power of elevation models. 17 taxa and 4 feeding guilds were modeled and the standardized study design allowed for a comparison across the assemblages. Results show that most taxa (14) and feeding guilds (3) can be predicted best by elevation with normalized RMSE values but only for three of those taxa and two of those feeding guilds the difference to other models is significant. Generally, modeling performances between different models vary only slightly for each assemblage. For the remaining, structural information at most showed little additional contribution to the performance. In summary, LiDAR observations can be used for animal species prediction. However, the effort and cost of aerial surveys are not always in proportion with the prediction quality, especially when the species distribution follows zonal patterns, and elevation information yields similar results.
This study compares the performance of three bias correction (BC) techniques in adjusting simulated precipitation estimates over Germany. The BC techniques are the multivariate quantile delta mapping (MQDM) where the grids are used as variables to incorporate the spatial dependency structure of precipitation in the bias correction; empirical quantile mapping (EQM) and, the linear scaling (LS) approach. Several metrics that include first to fourth moments and extremes characterized by the frequency of heavy wet days and return periods during boreal summer were applied to score the performance of the BC techniques. Our results indicate a strong dependency of the relative performances of the BC techniques on the choice of the regional climate model (RCM), the region, the season, and the metrics of interest. Hence, each BC technique has relative strengths and weaknesses. The LS approach performs well in adjusting the first moment but tends to fall short for higher moments and extreme precipitation during boreal summer. Depending on the season, the region and the RCM considered, there is a trade-off between the relative performances of the EQM and the MQDM in adjusting the simulated precipitation biases. However, the MQDM performs well across all considered metrics. Overall, the MQDM outperforms the EQM in improving the higher moments and in capturing the observed return level of extreme summer precipitation, averaged over Germany.
The fast and accurate yield estimates with the increasing availability and variety of global satellite products and the rapid development of new algorithms remain a goal for precision agriculture and food security. However, the consistency and reliability of suitable methodologies that provide accurate crop yield outcomes still need to be explored. The study investigates the coupling of crop modeling and machine learning (ML) to improve the yield prediction of winter wheat (WW) and oil seed rape (OSR) and provides examples for the Free State of Bavaria (70,550 km2), Germany, in 2019. The main objectives are to find whether a coupling approach [Light Use Efficiency (LUE) + Random Forest (RF)] would result in better and more accurate yield predictions compared to results provided with other models not using the LUE. Four different RF models [RF1 (input: Normalized Difference Vegetation Index (NDVI)), RF2 (input: climate variables), RF3 (input: NDVI + climate variables), RF4 (input: LUE generated biomass + climate variables)], and one semi-empiric LUE model were designed with different input requirements to find the best predictors of crop monitoring. The results indicate that the individual use of the NDVI (in RF1) and the climate variables (in RF2) could not be the most accurate, reliable, and precise solution for crop monitoring; however, their combined use (in RF3) resulted in higher accuracies. Notably, the study suggested the coupling of the LUE model variables to the RF4 model can reduce the relative root mean square error (RRMSE) from −8% (WW) and −1.6% (OSR) and increase the R
2 by 14.3% (for both WW and OSR), compared to results just relying on LUE. Moreover, the research compares models yield outputs by inputting three different spatial inputs: Sentinel-2(S)-MOD13Q1 (10 m), Landsat (L)-MOD13Q1 (30 m), and MOD13Q1 (MODIS) (250 m). The S-MOD13Q1 data has relatively improved the performance of models with higher mean R
2 [0.80 (WW), 0.69 (OSR)], and lower RRMSE (%) (9.18, 10.21) compared to L-MOD13Q1 (30 m) and MOD13Q1 (250 m). Satellite-based crop biomass, solar radiation, and temperature are found to be the most influential variables in the yield prediction of both crops.
Digital platforms, such as Amazon, represent the major beneficiaries of the Covid‐19 crisis. This study examines the role of digital platforms and their engagement in digitalisation initiatives targeting (small) brick‐and‐mortar retailers in Germany, thereby contributing to a better understanding of how digital platforms augment, substitute or reorganise physical retail spaces. This study applies a mixed‐method approach based on qualitative interviews, participant observation as well as media analysis. First, the study illustrates the controversial role of digital platforms by positioning themselves as supporting partners of the (offline) retailers, while simultaneously shifting power towards the platforms themselves. Second, digital platforms have established themselves not only as infrastructure providers but also as actors within these infrastructures, framing digital as well as physical retail spaces, inter alia due to their role as publicly legitimised retail advisers. Third, while institutions want to help retailers to survive, they simultaneously enhance retailers' dependency on digital platforms.
The expansion of renewable energies is being driven by the gradual phaseout of fossil fuels in order to reduce greenhouse gas emissions, the steadily increasing demand for energy and, more recently, by geopolitical events. The offshore wind energy sector is on the verge of a massive expansion in Europe, the United Kingdom, China, but also in the USA, South Korea and Vietnam. Accordingly, the largest marine infrastructure projects to date will be carried out in the upcoming decades, with thousands of offshore wind turbines being installed. In order to accompany this process globally and to provide a database for research, development and monitoring, this dissertation presents a deep learning-based approach for object detection that enables the derivation of spatiotemporal developments of offshore wind energy infrastructures from satellite-based radar data of the Sentinel-1 mission.
For training the deep learning models for offshore wind energy infrastructure detection, an approach is presented that makes it possible to synthetically generate remote sensing data and the necessary annotation for the supervised deep learning process. In this synthetic data generation process, expert knowledge about image content and sensor acquisition techniques is made machine-readable. Finally, extensive and highly variable training data sets are generated from this knowledge representation, with which deep learning models can learn to detect objects in real-world satellite data.
The method for the synthetic generation of training data based on expert knowledge offers great potential for deep learning in Earth observation. Applications of deep learning based methods can be developed and tested faster with this procedure. Furthermore, the synthetically generated and thus controllable training data offer the possibility to interpret the learning process of the optimised deep learning models.
The method developed in this dissertation to create synthetic remote sensing training data was finally used to optimise deep learning models for the global detection of offshore wind energy infrastructure. For this purpose, images of the entire global coastline from ESA's Sentinel-1 radar mission were evaluated. The derived data set includes over 9,941 objects, which distinguish offshore wind turbines, transformer stations and offshore wind energy infrastructures under construction from each other. In addition to this spatial detection, a quarterly time series from July 2016 to June 2021 was derived for all objects. This time series reveals the start of construction, the construction phase and the time of completion with subsequent operation for each object.
The derived offshore wind energy infrastructure data set provides the basis for an analysis of the development of the offshore wind energy sector from July 2016 to June 2021. For this analysis, further attributes of the detected offshore wind turbines were derived. The most important of these are the height and installed capacity of a turbine. The turbine height was calculated by a radargrammetric analysis of the previously detected Sentinel-1 signal and then used to statistically model the installed capacity. The results show that in June 2021, 8,885 offshore wind turbines with a total capacity of 40.6 GW were installed worldwide. The largest installed capacities are in the EU (15.2 GW), China (14.1 GW) and the United Kingdom (10.7 GW). From July 2016 to June 2021, China has expanded 13 GW of offshore wind energy infrastructure. The EU has installed 8 GW and the UK 5.8 GW of offshore wind energy infrastructure in the same period. This temporal analysis shows that China was the main driver of the expansion of the offshore wind energy sector in the period under investigation.
The derived data set for the description of the offshore wind energy sector was made publicly available. It is thus freely accessible to all decision-makers and stakeholders involved in the development of offshore wind energy projects. Especially in the scientific context, it serves as a database that enables a wide range of investigations. Research questions regarding offshore wind turbines themselves as well as the influence of the expansion in the coming decades can be investigated. This supports the imminent and urgently needed expansion of offshore wind energy in order to promote sustainable expansion in addition to the expansion targets that have been set.
Geoarchaeological information presented here pertains to a subsidiary Nile channel that once flowed west of the main Sebennitic distributary and discharged its water and sediments at Egypt’s then north-central deltaic coast. Periodical paleoclimatic episodes during the later Middle and Upper Holocene included decreased rainfall and increased aridity that reduced the Nile’s flow levels and thus likely disrupted nautical transport and anthropogenic activity along this channel. Such changes in this deltaic sector, positioned adjacent to the Levantine Basin in the Eastern Mediterranean, can be attributed to climatic shifts triggered as far as the North Atlantic to the west, and African highland source areas of the Egyptian Nile to the south. Of special interest in a study core recovered along the channel are several sediment sequences without anthropogenic material that are interbedded between strata comprising numerous potsherds. The former are interpreted here as markers of increased regional aridity and reduced Nile flow which could have periodically disrupted the regional distribution of goods and nautical activities. Such times occurred ~5000 years B.P., ~4200–4000 years B.P., ~3200–2800 years B.P., ~2300–2200 years B.P., and more recently. Periods comparable to these are also identified by altered proportions of pollen, isotopic and compositional components in different radiocarbon-dated Holocene cores recovered elsewhere in the Nile delta, the Levantine region to the east and north of Egypt, and in the Faiyum depression south of the delta.
Enhancing digital and precision agriculture is currently inevitable to overcome the economic and environmental challenges of the agriculture in the 21st century. The purpose of this study was to generate and compare management zones (MZ) based on the Sentinel-2 satellite data for variable rate application of mineral nitrogen in wheat production, calculated using different remote sensing (RS)-based models under varied soil, yield and crop data availability. Three models were applied, including (1) a modified “RS- and threshold-based clustering”, (2) a “hybrid-based, unsupervised clustering”, in which data from different sources were combined for MZ delineation, and (3) a “RS-based, unsupervised clustering”. Various data processing methods including machine learning were used in the model development. Statistical tests such as the Paired Sample T-test, Kruskal–Wallis H-test and Wilcoxon signed-rank test were applied to evaluate the final delineated MZ maps. Additionally, a procedure for improving models based on information about phenological phases and the occurrence of agricultural drought was implemented. The results showed that information on agronomy and climate enables improving and optimizing MZ delineation. The integration of prior knowledge on new climate conditions (drought) in image selection was tested for effective use of the models. Lack of this information led to the infeasibility of obtaining optimal results. Models that solely rely on remote sensing information are comparatively less expensive than hybrid models. Additionally, remote sensing-based models enable delineating MZ for fertilizer recommendations that are temporally closer to fertilization times.
The recently observed consistent loss of β-diversity across ecosystems indicates increasingly homogeneous communities in patches of landscapes, mainly caused by increasing land-use intensity. Biodiversity is related to numerous ecosystem functions and stability. Therefore, decreasing β-diversity is also expected to reduce multifunctionality. To assess the impact of homogenization and to develop guidelines to reverse its potentially negative effects, we combine expertise from forest science, ecology, remote sensing, chemical ecology and statistics in a collaborative and experimental β-diversity approach. Specifically, we will address the question whether the Enhancement of Structural Beta Complexity (ESBC) in forests by silviculture or natural disturbances will increase biodiversity and multifunctionality in formerly homogeneously structured production forests. Our approach will identify potential mechanisms behind observed homogenization-diversity-relationships and show how these translate into effects on multifunctionality. At eleven forest sites throughout Germany, we selected two districts as two types of small ‘forest landscapes’. In one of these two districts, we established ESBC treatments (nine differently treated 50x50 m patches with a focus on canopy cover and deadwood features). In the second, the control district, we will establish nine patches without ESBC. By a comprehensive sampling, we will monitor 18 taxonomic groups and measure 21 ecosystem functions, including key functions in temperate forests, on all patches. The statistical framework will allow a comprehensive biodiversity assessment by quantifying the different aspects of multitrophic biodiversity (taxonomical, functional and phylogenetic diversity) on different levels of biodiversity (α-, β-, γ-diversity). To combine overall diversity, we will apply the concept of multidiversity across the 18 taxa. We will use and develop new approaches for quantification and partitioning of multifunctionality at α- and β- scales. Overall, our study will herald a new research avenue, namely by experimentally describing the link between β-diversity and multifunctionality. Furthermore, we will help to develop guidelines for improved silvicultural concepts and concepts for management of natural disturbances in temperate forests reversing past homogenization effects.
The investigation of the Earth system and interplays between its components is of utmost importance to enhance the understanding of the impacts of global climate change on the Earth's land surface. In this context, Earth observation (EO) provides valuable long-term records covering an abundance of land surface variables and, thus, allowing for large-scale analyses to quantify and analyze land surface dynamics across various Earth system components. In view of this, the geographical entity of river basins was identified as particularly suitable for multivariate time series analyses of the land surface, as they naturally cover diverse spheres of the Earth. Many remote sensing missions with different characteristics are available to monitor and characterize the land surface. Yet, only a few spaceborne remote sensing missions enable the generation of spatio-temporally consistent time series with equidistant observations over large areas, such as the MODIS instrument.
In order to summarize available remote sensing-based analyses of land surface dynamics in large river basins, a detailed literature review of 287 studies was performed and several research gaps were identified. In this regard, it was found that studies rarely analyzed an entire river basin, but rather focused on study areas at subbasin or regional scale. In addition, it was found that transboundary river basins remained understudied and that studies largely focused on selected riparian countries. Moreover, the analysis of environmental change was generally conducted using a single EO-based land surface variable, whereas a joint exploration of multivariate land surface variables across spheres was found to be rarely performed.
To address these research gaps, a methodological framework enabling (1) the preprocessing and harmonization of multi-source time series as well as (2) the statistical analysis of a multivariate feature space was required. For development and testing of a methodological framework that is transferable in space and time, the transboundary river basins Indus, Ganges, Brahmaputra, and Meghna (IGBM) in South Asia were selected as study area, having a size equivalent to around eight times the size of Germany. These basins largely depend on water resources from monsoon rainfall and High Mountain Asia which holds the largest ice mass outside the polar regions. In total, over 1.1 billion people live in this region and in parts largely depend on these water resources which are indispensable for the world's largest connected irrigated croplands and further domestic needs as well. With highly heterogeneous geographical settings, these river basins allow for a detailed analysis of the interplays between multiple spheres, including the anthroposphere, biosphere, cryosphere, hydrosphere, lithosphere, and atmosphere.
In this thesis, land surface dynamics over the last two decades (December 2002 - November 2020) were analyzed using EO time series on vegetation condition, surface water area, and snow cover area being based on MODIS imagery, the DLR Global WaterPack and JRC Global Surface Water Layer, as well as the DLR Global SnowPack, respectively. These data were evaluated in combination with further climatic, hydrological, and anthropogenic variables to estimate their influence on the three EO land surface variables. The preprocessing and harmonization of the time series was conducted using the implemented framework. The resulting harmonized feature space was used to quantify and analyze land surface dynamics by means of several statistical time series analysis techniques which were integrated into the framework. In detail, these methods involved (1) the calculation of trends using the Mann-Kendall test in association with the Theil-Sen slope estimator, (2) the estimation of changes in phenological metrics using the Timesat tool, (3) the evaluation of driving variables using the causal discovery approach Peter and Clark Momentary Conditional Independence (PCMCI), and (4) additional correlation tests to analyze the human influence on vegetation condition and surface water area.
These analyses were performed at annual and seasonal temporal scale and for diverse spatial units, including grids, river basins and subbasins, land cover and land use classes, as well as elevation-dependent zones. The trend analyses of vegetation condition mostly revealed significant positive trends. Irrigated and rainfed croplands were found to contribute most to these trends. The trend magnitudes were particularly high in arid and semi-arid regions. Considering surface water area, significant positive trends were obtained at annual scale. At grid scale, regional and seasonal clusters with significant negative trends were found as well. Trends for snow cover area mostly remained stable at annual scale, but significant negative trends were observed in parts of the river basins during distinct seasons. Negative trends were also found for the elevation-dependent zones, particularly at high altitudes. Also, retreats in the seasonal duration of snow cover area were found in parts of the river basins. Furthermore, for the first time, the application of the causal discovery algorithm on a multivariate feature space at seasonal temporal scale revealed direct and indirect links between EO land surface variables and respective drivers. In general, vegetation was constrained by water availability, surface water area was largely influenced by river discharge and indirectly by precipitation, and snow cover area was largely controlled by precipitation and temperature with spatial and temporal variations. Additional analyses pointed towards positive human influences on increasing trends in vegetation greenness. The investigation of trends and interplays across spheres provided new and valuable insights into the past state and the evolution of the land surface as well as on relevant climatic and hydrological driving variables. Besides the investigated river basins in South Asia, these findings are of great value also for other river basins and geographical regions.