Remote sensing of impervious surface area and its interaction with land surface temperature variability in Pretoria, South Africa

Includes summary for chapter 1-5Pretoria, City of Tshwane (COT), Gauteng Province, South Africa is one of the cities that continues to experience rapid urban sprawl as a result of population growth and various land use, leading to the change of natural vegetation lands into impervious surface area (ISA). These are associated with transportation (paved roads, streets, highways, parking lots and sidewalks) and cemented buildings and rooftops, made of completely or partly impermeable artificial materials (e.g., asphalt, concrete, and brick). These landscapes influence the micro-climate (e.g., land surface temperature, LST) of Pretoria City as evidenced by the recent heat waves characterized by high temperature. Therefore, understanding ISA changes will provide information for city planning and environmental management. Conventionally, deriving ISA information has been dependent on field surveys and manual digitizing from hard copy maps, which is laborious and time-consuming. Remote sensing provides an avenue for deriving spatially explicit and timely ISA information. Numerous methods have been developed to estimate and retrieve ISA and LST from satellite imagery. There are limited studies focusing on the extraction of ISA and its relationship with LST variability across major cities in Africa. The objectives of the study were: (i) to explore suitable spectral indices to improve the delineation of built-up impervious surface areas from very high resolution multispectral data (e.g., WorldView-2), (ii) to examine exposed rooftop impervious surface area based on different colours, and their interplay with surface temperature variability, (iii) to determine if the spatio-temporal built-up ISA distribution pattern in relation to elevation influences urban heat island (UHI) extent using an optimal analytical scale and (iv) to assess the spatio-temporal change characteristics of ISA expansion using the corresponding surface temperature (LST) at selected administrative subplace units (i.e., local region scale). The study objectives were investigated using remote sensing data such as WorldView-2 (a very high-resolution multispectral sensor), medium resolution Landsat-5 Thematic Mapper (TM) and Landsat-8 OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor) at multiple scales. The ISA mapping methods used in this study can be grouped into two major categories: (i) the classification-based approach consisting of an object-based multi-class classification with overall accuracy ~90.4% and a multitemporal pixel-based binary classification. The latter yielded an area under the receiver operating characteristic curve (AUROC) = 0.8572 for 1995, AUROC = 0.8709 for 2005, AUROC = 0.8949 for 2015. (ii) the spectral index-based approach such as a new built-up extraction index (NBEI) derived in this study which yielded a high AUROC = ~0.82 compared to Built-up Area Index (BAI) (AUROC = ~0.73), Built-up spectral index (BSI) (AUROC = ~0.78), Red edge / Green Index (RGI) (AUROC = ~0.71) and WorldView-Built-up Index (WV-BI) (AUROC = ~0.67). The multitemporal built-up Index (BUI) also estimated with AUROC = 0.8487 for 1993, AUROC = 0.8302 for 2003, AUROC = 0.8790 for 2013. This indicates that all these methods employed, mapped ISA with high predictive accuracy from remote sensing data. Furthermore, the single-channel algorithm (SCA) was employed to retrieve LST from the thermal infrared (TIR) band of the Landsat images. The LST overall retrieval error for the entire study generally was quite low (overall root mean square RMSE ≤ ~1.48OC), which signifies that the Landsat TIR used provided good results for further analysis. In conclusion, the study showed the potential of multispectral remote sensing data to quantify ISA and evaluate its interaction with surface temperature variability despite the complex urban landscape in Pretoria. Also, using impervious surface LST as a complementary metric in this research helped to reveal urban heat island distribution and improve understanding of the spatio-temporal developing trend of urban expansion at a local spatial scale.Rapid urbanization because of population growth has led to the conversion of natural lands into large man-made landscapes which affects the micro-climate. Rooftop reflectivity, material, colour, slope, height, aspect, elevation are factors that potentially contribute to temperature variability. Therefore, strategically designed rooftop impervious surfaces have the potential to translate into significant energy, long-term cost savings, and health benefits. In this experimental study, we used the semi-automated Environment for Visualizing Images (ENVI) Feature Extraction that uses an object-based image analysis approach to classify rooftop based on colours from WorldView-2 (WV-2) image with overall accuracy ~90.4% and kappa coefficient ~0.87 respectively. The daytime retrieved surface temperatures were derived from 15m pan-sharpened Landsat 8 TIRS with a range of ~14.6OC to ~65OC (retrieval error = 0.38OC) for the same month covering Lynwood Ridge a residential area in Pretoria. Thereafter, the relationship between the rooftops and surface temperature (LST) were examined using multivariate statistical analysis. The results of this research reveal that the interaction between the applicable rooftop explanatory features (i.e., reflectance, texture measures and topographical properties) can explain over 22.10% of the variation in daytime rooftop surface temperatures. Furthermore, analysis of spatial distribution between mean daytime surface temperature and the residential rooftop indicated that the red, brown and green roof surfaces show lower LST values due to high reflectivity, high emissivity and low heat capacity during the daytime. The study concludes that in any study related to the spatial distribution of rooftop impervious surface area surface temperature, effect of various explanatory variables must be considered. The results of this experimental study serve as a useful approach for further application in urban planning and sustainable development.Evaluating changes in built-up impervious surface area (ISA) to understand the urban heat island (UHI) extent is valuable for governments in major cities in developing countries experiencing rapid urbanization and industrialization. This work aims at assessing built-up ISA spatio-temporal and influence on land surface temperature (LST) variability in the context of urban sprawl. Landsat-5 Thematic Mapper (TM) and Landsat-8 OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor) were used to quantify ISA using built-up Index (BUI) and spatio-temporal dynamics from 1993-2013. Thereafter using a suitable analytical sampling scale that represents the estimated ISA-LST, we examined its distribution in relation to elevation using the Shuttle Radar Topography Mission (SRTM) and also create Getis-Ord Gi* statistics hotspot maps to display the UHI extent. The BUI ISA extraction results show a high predictive accuracy with area under the receiver operating characteristic curve, AUROC = 0.8487 for 1993, AUROC = 0.8302 for 2003, AUROC = 0.8790 for 2013. The ISA spatio-temporal changes within ten years interval time frame results revealed a 14% total growth rate during the study year. Based on a suitable analytical scale (90x90) for the hexagon polygon grid, the majority of ISA distribution across the years was at an elevation range of between >1200m – 1600m. Also, Getis-Ord Gi* statistics hotspot maps revealed that hotspot regions expanded through time with a total growth rate of 19% and coldspot regions decreased by 3%. Our findings can represent useful information for policymakers by providing a scientific basis for sustainable urban planning and management.Over the years, rapid urban growth has led to the conversion of natural lands into large man-made landscapes due to enhanced political and economic growth. This study assessed the spatio-temporal change characteristics of impervious surface area (ISA) expansion using its surface temperature (LST) at selected administrative subplace units (i.e., local region scale). ISA was estimated for 1995, 2005 and 2015 from Landsat-5 Thematic Mapper (TM) and Landsat-8 OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor) images using a Random Forest (RF) algorithm. The spatio-temporal trends of ISA were assessed using an optimal analytical scale to aggregate ISA LST coupled with weighted standard deviational ellipse (SDE) method. The ISA was quantified with high predictive accuracy (i.e., AUROC = 0.8572 for 1995, AUROC = 0.8709 for 2005, AUROC = 0.8949 for 2015) using RF classifier. More than 70% of the selected administrative subplaces in Pretoria experienced an increase in growth rate (415.59%) between 1995 and 2015. LST computations from the Landsat TIRS bands yielded good results (RMSE = ~1.44OC, 1.40OC, ~0.86OC) for 1995, 2005 and 2015 respectively. Based on the hexagon polygon grid (90x90), the aggregated ISA surface temperature weighted SDE analysis results indicated ISA expansion in different directions at the selected administrative subplace units. Our findings can represent useful information for policymakers in evaluating urban development trends in Pretoria, City of Tshwane (COT).Globally, the unprecedented increase in population in many cities has led to rapid changes in urban landscape, which requires timely assessments and monitoring. Accurate determination of built-up information is vital for urban planning and environmental management. Often, the determination of the built-up area information has been dependent on field surveys, which is laborious and time-consuming. Remote sensing data is the only option for deriving spatially explicit and timely built-up area information. There are few spectral indices for built-up areas and often not accurate as they are specific to impervious material, age, colour, and thickness, especially using higher resolution images. The objective of this study is to test the utility of a new built-up extraction index (NBEI) using WorldView-2 to improve built-up material mapping irrespective of material type, age and colour. The new index was derived from spectral bands such as Green, Red edge, NIR1 and NIR2 bands that profoundly explain the variation in built-up areas on WorldView-2 image (WV-2). The result showed that NBEI improves the extraction of built-up areas with high accuracy (area under the receiver operating characteristic curve, AUROC = ~0.82) compared to the existing indices such as Built-up Area Index (BAI) (AUROC = ~0.73), Built-up spectral index (BSI) (AUROC = ~0.78 ), Red edge / Green Index (RGI) (AUROC = ~0.71) and WorldView-Built-up Index (WV-BI) (AUROC = ~0.67). The study demonstrated that the new built-up index could extract built-up areas using high-resolution images. The performance of NBEI could be attributed to the fact that it is not material specific, and would be necessary for urban area mapping.Environmental SciencesD. Phil. (Environmental Sciences

Remote Sensing of Natural Hazards

Each year, natural hazards such as earthquakes, cyclones, flooding, landslides, wildfires, avalanches, volcanic eruption, extreme temperatures, storm surges, drought, etc., result in widespread loss of life, livelihood, and critical infrastructure globally. With the unprecedented growth of the human population, largescale development activities, and changes to the natural environment, the frequency and intensity of extreme natural events and consequent impacts are expected to increase in the future.Technological interventions provide essential provisions for the prevention and mitigation of natural hazards. The data obtained through remote sensing systems with varied spatial, spectral, and temporal resolutions particularly provide prospects for furthering knowledge on spatiotemporal patterns and forecasting of natural hazards. The collection of data using earth observation systems has been valuable for alleviating the adverse effects of natural hazards, especially with their near real-time capabilities for tracking extreme natural events. Remote sensing systems from different platforms also serve as an important decision-support tool for devising response strategies, coordinating rescue operations, and making damage and loss estimations.With these in mind, this book seeks original contributions to the advanced applications of remote sensing and geographic information systems (GIS) techniques in understanding various dimensions of natural hazards through new theory, data products, and robust approaches

Ph.D. and Master’s thesis in geography (2012-2017)

This list follows those published in Belgeo, 2011, 3-4 (2007-2011), 2006, 4 (2004-2006), 2003, 2 (1999-2003) and in the Bulletin de la Société Belge d’Études géographiques/ Tijdschrift van de Belgische Aardrijkskundige Vereeniging, 1998, 2 (1994-1998). This publication allows drawing up some statistics on the evolution of the number of licences/masters and Ph. D. defended in the Belgian Universities from the mid-Nineties. Excluding the international Four Cities program (presently very succes..

UAVs for the Environmental Sciences

This book gives an overview of the usage of UAVs in environmental sciences covering technical basics, data acquisition with different sensors, data processing schemes and illustrating various examples of application

Uncertainties in Digital Elevation Models: Evaluation and Effects on Landform and Soil Type Classification

Digital elevation models (DEMs) are a widely used source for the digital representation of the Earth's surface in a wide range of scientific, industrial and military applications. Since many processes on Earth are influenced by the shape of the relief, a variety of different applications rely on accurate information about the topography. For instance, DEMs are used for the prediction of geohazards, climate modelling, or planning-relevant issues, such as the identification of suitable locations for renewable energies. Nowadays, DEMs can be acquired with a high geometric resolution and over large areas using various remote sensing techniques, such as photogrammetry, RADAR, or laser scanning (LiDAR). However, they are subject to uncertainties and may contain erroneous representations of the terrain. The quality and accuracy of the topographic representation in the DEM is crucial, as the use of an inaccurate dataset can negatively affect further results, such as the underestimation of landslide hazards due to a too flat representation of relief in the elevation model. Therefore, it is important for users to gain more knowledge about the accuracy of a terrain model to better assess the negative consequences of DEM uncertainties on further analysis results of a certain research application. A proper assessment of whether the purchase or acquisition of a highly accurate DEM is necessary or the use of an already existing and freely available DEM is sufficient to achieve accurate results is of great qualitative and economic importance. In this context, the first part of this thesis focuses on extending knowledge about the behaviour and presence of uncertainties in DEMs concerning terrain and land cover. Thus, the first two studies of this dissertation provide a comprehensive vertical accuracy analysis of twelve DEMs acquired from space with spatial resolutions ranging from 5 m to 90 m. The accuracy of these DEMs was investigated in two different regions of the world that are substantially different in terms of relief and land cover. The first study was conducted in the hyperarid Chilean Atacama Desert in northern Chile, with very sparse land cover and high elevation differences. The second case study was conducted in a mid-latitude region, the Rur catchment in the western part of Germany. This area has a predominantly flat to hilly terrain with relatively diverse and dense vegetation and land cover. The DEMs in both studies were evaluated with particular attention to the influence of relief and land cover on vertical accuracy. The change of error due to changing slope and land cover was quantified to determine an average loss of accuracy as a function of slope for each DEM. Additionally, these values were used to derive relief-adjusted error values for different land cover classes. The second part of this dissertation addresses the consequences that different spatial resolutions and accuracies in DEMs have on specific applications. These implications were examined in two exemplary case studies. In a geomorphometric case study, several DEMs were used to classify landforms by different approaches. The results were subsequently compared and the accuracy of the classification results with different DEMs was analysed. The second case study is settled within the field of digital soil mapping. Various soil types were predicted with machine learning algorithms (random forest and artificial neural networks) using numerous relief parameters derived from DEMs of different spatial resolutions. Subsequently, the influence of high and low resolution DEMs with the respectively derived land surface parameters on the prediction results was evaluated. The results on the vertical accuracy show that uncertainties in DEMs can have diverse reasons. Besides the spatial resolution, the acquisition technique and the degree of improvements made to the dataset significantly impact the occurrence of errors in a DEM. Furthermore, the relief and physical objects on the surface play a major role for uncertainties in DEMs. Overall, the results in steeper areas show that the loss of vertical accuracy is two to three times higher for a 90 m DEM than for DEMs of higher spatial resolutions. While very high resolution DEMs of 12 m spatial resolution or higher only lose about 1 m accuracy per 10° increase in slope steepness, 30 m DEMs lose about 2 m on average, and 90 m DEMs lose more than 3 m up to 6 m accuracy. However, the results also show significant differences for DEMs of identical spatial resolution depending on relief and land cover. With regard to different land cover classes, it can be stated that mid-latitude forested and water areas cause uncertainties in DEMs of about 6 m on average. Other tested land cover classes produced minor errors of about 1 – 2 m on average. The results of the second part of this contribution prove that a careful selection of an appropriate DEM is more crucial for certain applications than for others. The choice of different DEMs greatly impacted the landform classification results. Results from medium resolution DEMs (30 m) achieved up to 30 % lower overall accuracies than results from high resolution DEMs with a spatial resolution of 5 m. In contrast to the landform classification results, the predicted soil types in the second case study showed only minor accuracy differences of less than 2 % between the usage of a spatial high resolution DEM (15 m) and a low resolution 90 m DEM. Finally, the results of these two case studies were compared and discussed with other results from the literature in other application areas. A summary and assessment of the current state of knowledge about the impact of a particular chosen terrain model on the results of different applications was made. In summary, the vertical accuracy measures obtained for each DEM are a first attempt to determine individual error values for each DEM that can be interpreted independently of relief and land cover and can be better applied to other regions. This may help users in the future to better estimate the accuracy of a tested DEM in a particular landscape. The consequences of elevation model selection on further results are highly dependent on the topic of the study and the study area's level of detail. The current state of knowledge on the impact of uncertainties in DEMs on various applications could be established. However, the results of this work can be seen as a first step and more work is needed in the future to extend the knowledge of the effects of DEM uncertainties on further topics that have not been investigated to date

VGC 2023 - Unveiling the dynamic Earth with digital methods: 5th Virtual Geoscience Conference: Book of Abstracts

Conference proceedings of the 5th Virtual Geoscience Conference, 21-22 September 2023, held in Dresden. The VGC is a multidisciplinary forum for researchers in geoscience, geomatics and related disciplines to share their latest developments and applications.:Short Courses 9 Workshops Stream 1 10 Workshop Stream 2 11 Workshop Stream 3 12 Session 1 – Point Cloud Processing: Workflows, Geometry & Semantics 14 Session 2 – Visualisation, communication & Teaching 27 Session 3 – Applying Machine Learning in Geosciences 36 Session 4 – Digital Outcrop Characterisation & Analysis 49 Session 5 – Airborne & Remote Mapping 58 Session 6 – Recent Developments in Geomorphic Process and Hazard Monitoring 69 Session 7 – Applications in Hydrology & Ecology 82 Poster Contributions 9

Handbook of Mathematical Geosciences

This Open Access handbook published at the IAMG's 50th anniversary, presents a compilation of invited path-breaking research contributions by award-winning geoscientists who have been instrumental in shaping the IAMG. It contains 45 chapters that are categorized broadly into five parts (i) theory, (ii) general applications, (iii) exploration and resource estimation, (iv) reviews, and (v) reminiscences covering related topics like mathematical geosciences, mathematical morphology, geostatistics, fractals and multifractals, spatial statistics, multipoint geostatistics, compositional data analysis, informatics, geocomputation, numerical methods, and chaos theory in the geosciences

Pattern to process: methodological investigations into the formation and interpretation of spatial patterns in archaeological landscapes

My research has shown that the type of regional archaeological data analysis required by landscape archaeological approaches is an area where both theory and method are still in their infancy. High-level theories about the occurrence, scope, and effects of processes such as centralization, urbanization, and Hellenization/Romanization cannot yet be supported by middle range theory, which itself cannot be developed until the basic business of generating information of sufficient quality about the archaeological record has been tackled. Currently, archaeological data can be made to fit almost any interpretation generated, ultimately, on the basis of the ancient written sources. If we are to escape from this selfreinforcing cycle, research should perhaps no longer be focused on the classical themes generated by culture-historical approaches, but should seek its own proper field of operation. In the area of methods and methodology, I have demonstrated the pervasive influence of systematic research and visibility biases on the patterns that are present in the archaeological data generated over the past 50 years or so. There are mechanisms at work, both in the traditional archaeological interpretation of limited numbers of excavated sites and historical sources, and in the landscape archaeological approach, that cause the systematic undervaluation of unobtrusive remains. The significance of systematic biases in both the coarse site-based data sets resulting from desktop and ‘topographic’ studies and the more detailed site-based or ‘continuous’ data resulting from intensive field surveys has become much clearer as a result of the studies reported here. This should have practical consequences for the ways in which we study the existing archaeological record, plan future landscape archaeological research, and conduct field surveys. Site databases, the traditional starting point for regional archaeological studies, can no longer be taken at face value; rather, they require careful source criticism before being used to support specific arguments and hypotheses about settlement and land use dynamics. My studies have also shown that future data collection, whether through field survey, excavation or other methods, has to take place in a much more methodical manner if we are to produce data that are sufficiently standardized to be successfully exchanged, compared, and interpreted by others – guidelines for which should become embodied in an international standard defining ‘best practice in landscape archaeology’.

Quantifying the effects of variable selection, spatial scale and spatial data quality in marine benthic habitat mapping

Mapping benthic habitats has become critical in many contexts like conservation and management. While marine habitat mapping methods strongly rely on tools and methods from geography and geomatics, habitat mapping practitioners with a background outside of these specialized areas do not always have a full understanding of the spatial concepts behind these tools and methods. This phenomenon is amplified when marine geomorphometry, the science used to quantify seafloor terrain characteristics, is integrated into the marine habitat mapping workflow. This dissertation reviews the use of spatial concepts in the field of marine benthic habitat mapping; many concepts are poorly understood or poorly implemented in the habitat mapping workflow, among which spatial scale and spatial data quality stand out as being of particular importance. While geomorphometry is commonly used in marine benthic habitat mapping, no framework existed to test which terrain attributes should be used as surrogates of species distribution, leading to an inability to compare results from different studies. This dissertation explores different options for terrain attribute selection and proposes an optimal combination that can be used as standard in all habitat mapping studies. This selection is then tested using two approaches to benthic habitat mapping and is shown to perform better than others. Bathymetric data, the primary input for marine geomorphometry analyses and one of the main data inputs for habitat mapping, are commonly impacted by data acquisition artefacts. Very little work has been done on trying to understand how these artefacts propagate throughout the habitat mapping workflow. The impact of artefacts on the bathymetry and its derived terrain attributes is described, and it is shown that artefacts modify the spatial and statistical distributions of depth and terrain attribute values. However, when these affected data are used in habitat mapping, their impact is not always predictable. Some artefacts were found to sometimes inflate measures of accuracy and performance and sometimes decrease them. Overall, habitat maps were shown to be very sensitive to the effects of variable selection, spatial scale and data quality, and as such have serious implications when they are used to inform decision-making, for instance in marine conservation and management. This dissertation raises awareness about these issues and highlights the need for careful integration of spatial data in habitat mapping practices