Pollution of surface and groundwater, due to improper land management, has become a major problem worldwide. Integrated watershed modelling provides a tool for the understanding of the processes governing water and matter transport at different scales within the watershed. The Soil Water Assessment Tool (SWAT) has been successfully utilized for the combined modelling of water fluxes and quality within a large range of scales and environmental conditions across the world. For suitable assessments integrated watershed models require large data sets of measured information for both model parameterization as for model calibration and validation. Data scarcity represents a serious limitation to the use of hydrologic models for supporting decision making processes, and may lead unsupported statements, poor statistics, misrepresentations, and, ultimately, to inappropriate measures for integrated water resources management efforts. In particular, the importance of spatially distributed soil information is often overlooked. In this thesis the eco-hydrological SWAT model was been applied to assess the water balance and diffuse pollution loadings of two rivers within a rural context at the mesoscale watershed level: 1) the Western Bug River, Ukraine, 2) the Águeda River, Portugal. Both watersheds in focus serve as examples for areas where the amount and quality of the measured data hinders a strait forward hydrologic modelling assessment. The Dobrotvir watershed (Western Bug River, Ukriane) is an example of such a region. In the former Soviet Union, soil classification primarily focused on soils of agricultural importance, whereas, forested, urban, industrial, and shallow soil territories were left underrepresented in the classification systems and resulting soil maps. Similarly the forest-dominated Águeda watershed in North-Central Portugal is a second example of a region with serious soil data availability limitations. Through the use of pedotransfer functions (PTFs) and the construction of soil-landscape models the data gaps could be successfully diminished, allowing a subsequent integrated watershed modelling approach. A valuable tool for the data gap closure was the fuzzy logic Soil Land Inference Model (SoLIM) which, combined with information from several soil surveys, was used to create improved maps. In the Dobrotvir watershed the fuzzy approach was used to close the gaps of the existing soil map, while in the Águeda watershed a new soil properties map, based upon the effective soil depths of the landscape, was constructed. While the water balance simulation in both study areas was successful, a calibration parameter ensemble approach was tested for the Águeda watershed. In the common modelling practice the individual best simulation and best parameter set is considered, the tested approach involved merging individual model outputs from numerous acceptable parameter sets, tackling the problematic of parameter equifinality. This procedure was tested for both original soil map and the newly derived soil map with differentiation of soil properties. It was noticeable that a better model set-up, with a better representation of the soil spatial distribution, was reflected in tighter model output spreads and narrower parameter distances. A further challenge was the calibration of water quality parameters, namely nitrate-N in the Dobrotvir watershed and sediment loads in the Águeda watershed. The limited amount of water quality observations were handled by assessing and by process verification at the smallest modelling unit, the hydrological response unit (HRU). The ruling hydrological processes could be depicted by combining own measured data and modelling outputs. The management scenario simulations showed the anticipated response to changes in management and reflected the rational spatial variation within the watershed reasonably well. The impacts of the different intervention options were evaluated on water balance, nitrate-N export and sediment yield at the watershed, sub-watershed and, when feasible, HRU level. This thesis covers two regional case studies with particular data limitations and specific processes of water and matter fluxes. Still, data reliability is a problem across the globe. This thesis demonstrates how relevant it is to tackle shortages of spatially differentiated soil information. The considered approaches contribute toward more reliable model predictions. Furthermore, the tested methods are transferable to other regions with differing landscape and climate conditions with similar problems of data scarcity, particularly soil spatially differentiated information.

Human activities on land have a direct and cumulative impact on water and other natural resources within a catchment. This land-use change can have hydrological consequences on the local and regional scales. Sound catchment assessment is not only critical to understanding processes and functions but also important in identifying priority management areas. The overarching goal of this doctoral thesis was to design a methodological framework for catchment assessment (dependent upon data availability) and propose practical catchment management strategies for sustainable water resources management. The Nyangores and Ruiru reservoir catchments located in Kenya, East Africa were used as case studies. A properly calibrated Soil and Water Assessment Tool (SWAT) hydrologic model coupled with a generic land-use optimization tool (Constrained Multi-Objective Optimization of Land-use Allocation-CoMOLA) was applied to identify and quantify functional trade-offs between environmental sustainability and food production in the ‘data-available’ Nyangores catchment. This was determined using a four-dimension objective function defined as (i) minimizing sediment load, (ii) maximizing stream low flow and (iii and iv) maximizing the crop yields of maize and soybeans, respectively. Additionally, three different optimization scenarios, represented as i.) agroforestry (Scenario 1), ii.) agroforestry + conservation agriculture (Scenario 2) and iii.) conservation agriculture (Scenario 3), were compared. For the data-scarce Ruiru reservoir catchment, alternative methods using digital soil mapping of soil erosion proxies (aggregate stability using Mean Weight Diameter) and spatial-temporal soil loss analysis using empirical models (the Revised Universal Soil Loss Equation-RUSLE) were used. The lack of adequate data necessitated a data-collection phase which implemented the conditional Latin Hypercube Sampling. This sampling technique reduced the need for intensive soil sampling while still capturing spatial variability. The results revealed that for the Nyangores catchment, adoption of both agroforestry and conservation agriculture (Scenario 2) led to the smallest trade-off amongst the different objectives i.e. a 3.6% change in forests combined with 35% change in conservation agriculture resulted in the largest reduction in sediment loads (78%), increased low flow (+14%) and only slightly decreased crop yields (3.8% for both maize and soybeans). Therefore, the advanced use of hydrologic models with optimization tools allows for the simultaneous assessment of different outputs/objectives and is ideal for areas with adequate data to properly calibrate the model. For the Ruiru reservoir catchment, digital soil mapping (DSM) of aggregate stability revealed that susceptibility to erosion exists for cropland (food crops), tea and roadsides, which are mainly located in the eastern part of the catchment, as well as deforested areas on the western side. This validated that with limited soil samples and the use of computing power, machine learning and freely available covariates, DSM can effectively be applied in data-scarce areas. Moreover, uncertainty in the predictions can be incorporated using prediction intervals. The spatial-temporal analysis exhibited that bare land (which has the lowest areal proportion) was the largest contributor to erosion. Two peak soil loss periods corresponding to the two rainy periods of March–May and October–December were identified. Thus, yearly soil erosion risk maps misrepresent the true dimensions of soil loss with averages disguising areas of low and high potential. Also, a small portion of the catchment can be responsible for a large proportion of the total erosion. For both catchments, agroforestry (combining both the use of trees and conservation farming) is the most feasible catchment management strategy (CMS) for solving the major water quantity and quality problems. Finally, the key to thriving catchments aiming at both sustainability and resilience requires urgent collaborative action by all stakeholders. The necessary stakeholders in both Nyangores and Ruiru reservoir catchments must be involved in catchment assessment in order to identify the catchment problems, mitigation strategies/roles and responsibilities while keeping in mind that some risks need to be shared and negotiated, but so will the benefits.:TABLE OF CONTENTS DECLARATION OF CONFORMITY........................................................................ i DECLARATION OF INDEPENDENT WORK AND CONSENT ............................. ii LIST OF PAPERS ................................................................................................. iii ACKNOWLEDGEMENTS ..................................................................................... iv THESIS AT A GLANCE ......................................................................................... v SUMMARY ............................................................................................................ vi List of Figures......................................................................................................... x List of Tables........................................................................................................... x ABBREVIATION..................................................................................................... xi PART A: SYNTHESIS 1. INTRODUCTION ............................................................................................... 1 1.1 Catchment management ...................................................................................1 1.2 Tools to support catchment assessment and management ..............................4 1.3 Catchment management strategies (CMSs)......................................................9 1.4 Concept and research objectives.......................................................................11 2. MATERIAL AND METHODS................................................................................15 2.1. STUDY AREA ..................................................................................................15 2.1.1. Nyangores catchment ...................................................................................15 2.1.2. Ruiru reservoir catchment .............................................................................17 2.2. Using SWAT conceptual model and land-use optimization ..............................19 2.3. Using soil erosion proxies and empirical models ..............................................21 3. RESULTS AND DISCUSSION..............................................................................24 3.1. Assessing multi-metric calibration performance using the SWAT model...........25 3.2. Land-use optimization using SWAT-CoMOLA for the Nyangores catchment. ..26 3.3. Digital soil mapping of soil aggregate stability ..................................................28 3.4. Spatio-temporal analysis using the revised universal soil loss equation (RUSLE) 29 4. CRITICAL ASSESSMENT OF THE METHODS USED ......................................31 4.1. Assessing suitability of data for modelling and overcoming data challenges...31 4.2. Selecting catchment management strategies based on catchment assessment . 35 5. CONCLUSION AND RECOMMENDATIONS ....................................................36 6. REFERENCES ............................ .....................................................................38 PART B: PAPERS PAPER I .................................................................................................................47 PAPER II ................................................................................................................59 PAPER III ...............................................................................................................74 PAPER IV ...............................................................................................................88

Conselho Nacional de Desenvolvimento Científico e Tecnológico / The possibility of mapping soil properties using soil spatial prediction functions (SSPFe) is a reality. But is it possible to SSPFe to estimate soil properties such as the particlesize distribution (psd) in a young, unstable and geologically complex geomorphologic surface? What would be considered a good performance in such situation and what alternatives do we have to improve it? With the present study I try to find answers to such questions. To do so I used a set of 339 soil samples from a small catchment of the hillslope areas of central Rio Grande do Sul. Multiple linear regression models were built using landsurface parameters (elevation, convergence index, stream power index). The SSPFe explained more than half of data variance. Such performance is similar to that of the conventional soil mapping approach. For some size-fractions the SSPFe performance can reach 70%. Largest uncertainties are observed in areas of larger geological heterogeneity. Therefore, significant improvements in the predictions can only be achieved if accurate geological data is made available. Meanwhile, SSPFe built on land-surface parameters are efficient in estimating the psd of the soils in regions of complex geology. However, there still are questions that I couldn t answer! Is soil mapping important to solve the main social and environmental issues of our time? What if our activities were subjected to a social control as in a direct democracy, would they be worthy of receiving any attention? / A possibilidade de mapear as propriedades dos solos através do uso de funções de predição espacial de solos (FPESe) é uma realidade. Mas seria possível construir FPESe para estimar propriedades como a distribuição do tamanho de partículas do solo (dtp) em um superfície geomorfológica jovem e instável, com elevada complexidade geológica e pedológica? O que seria considerado um bom desempenho nessas condições e que alternativas temos para melhorá-lo? Com esse trabalho tento encontrar respostas para essas questões. Para isso utilizei um conjunto de 339 amostras de solo de uma pequena bacia hidrográfica de encosta da região Central do RS. Modelos de regressão linear múltiplos foram construídos com atributos de terreno (elevação, índice de convergência, índice de potência de escoamento). As FPESe explicaram mais da metade da variância dos dados. Tal desempenho é semelhante àquele da abordagem tradicional de mapeamento de solos. Para algumas frações de tamanho o desempenho das FPESe pode chegar a 70%. As maiores incertezas ocorrem nas áreas de maior heterogeneidade geológica. Assim, melhorias significativas nas predições somente poderão ser alcançadas se dados geológicos acurados forem disponibilizados. Enquanto isso, FPESe construídas a partir de atributos de terreno são eficientes em estimar a dtp de solos de regiões com geologia complexa e elevada instabilidade. Mas restam dúvidas que não consegui resolver! O mapeamento de solos é importante para a resolução dos principais problemas sociais e ambientais do nosso tempo? E se nossas atividades estivessem submetidas ao controle da população como em uma democracia direta, seriam elas dignas de receber atenção?

The objective of the thesis is the model-based analysis of spatial patterns of decomposition properties on the forested slopes of the montane level (ca. 1200-2200 m a.s.l.) in a study area in the Italian Alps (Val di Sole / Val di Rabbi, Autonomous Province of Trento). The analysis includes humus forms and enchytraeid assemblages as well as pH values, activities of extracellular enzymes and C/N ratios of the topsoil. The first aim is to develop, test and apply data-based techniques for spatial modelling of soil ecological parameters. This methodological approach is based on the concept of digital soil mapping. The second aim is to reveal the relationships between humus forms, soil organisms and soil microbiological parameters in the study area. The third aim is to analyze if the spatial patterns of indicators of decomposition differ between the landscape scale and the slope scale. At the landscape scale, sample data from six sites are used, covering three elevation levels at both north- and south-facing slopes. A knowledge-based approach that combines a decision tree analysis with the construction of fuzzy membership functions is introduced for spatial modelling. According to the sampling design, elevation and slope exposure are the explanatory variables. The investigations at the slope scale refer to one north-facing and one south-facing slope, with 30 sites occurring on each slope. These sites have been derived using conditioned Latin Hypercube Sampling, and thus reasonably represent the environmental conditions within the study area. Predictive maps have been produced in a purely data-based approach with random forests. At both scales, the models indicate a high variability of spatial decomposition patterns depending on the elevation and the slope exposure. In general, sites at high elevation on north-facing slopes almost exclusively exhibit the humus forms Moder and Mor. Sites on south-facing slopes and at low elevation exhibit also Mull and Amphimull. The predictions of those enchytraeid species characterized as Mull and Moder indicators match the occurrence of the corresponding humus forms well. Furthermore, referencing the mineral topsoil, the predictive models show increasing pH values, an increasing leucine-aminopeptidase activity, an increasing ratio alkaline/acid phosphomonoesterase activity and a decreasing C/N ratio from north-facing to south-facing slopes and from high to low elevation. The predicted spatial patterns of indicators of decomposition are basically similar at both scales. However, the patterns are predicted in more detail at the slope scale because of a larger data basis and a higher spatial precision of the environmental covariates. These factors enable the observation of additional correlations between the spatial patterns of indicators of decomposition and environmental influences, for example slope angle and curvature. Both the corresponding results and broad model evaluations have shown that the applied methods are generally suitable for modelling spatial patterns of indicators of decomposition in a heterogeneous high mountain environment. The overall results suggest that the humus form can be used as indicator of organic matter decomposition processes in the investigated high mountain area.

Functions and services provided by soils play an important role for numerous sustainable development goals involving mainly food supply and environmental health. In many regions of the Earth, water erosion is a major threat to soil functions and is mostly related to land-use change or poor agricultural management. Selecting proper soil management practices requires site-specific indicators such as water erosion, which follow a spatio-temporal variation. The aim of this study was to develop monthly soil erosion risk maps for the data-scarce catchment of Ruiru drinking water reservoir located in Kenya. Therefore, the Revised Universal Soil Loss Equation complemented with the cubist–kriging interpolation method was applied. The erodibility map created with digital soil mapping methods (R2 = 0.63) revealed that 46% of the soils in the catchment have medium to high erodibility. The monthly erosion rates showed two distinct potential peaks of soil loss over the course of the year, which are consistent with the bimodal rainy season experienced in central Kenya. A higher soil loss of 2.24 t/ha was estimated for long rains (March–May) as compared to 1.68 t/ha for short rains (October–December). Bare land and cropland are the major contributors to soil loss. Furthermore, spatial maps reveal that areas around the indigenous forest on the western and southern parts of the catchment have the highest erosion risk. These detected erosion risks give the potential to develop efficient and timely soil management strategies, thus allowing continued multi-functional use of land within the soil–food–water nexus.

Soil carbon is the largest sink and source of the global carbon cycle and is disturbed by several natural, anthropogenic, and environmental factors. The global increase of atmospheric CO2 affects soil carbon cycling through varied biogeochemical processes. The first chapter is a compilation of current information on potential factors triggering soil acidification and weathering mechanisms under elevated CO2 and their consequences on soil inorganic carbon (SIC) pool and quality. Soil water content and precipitation were critical factors influencing elevated CO2 effects on the SIC pool. The second chapter examines a detailed column experiment in which six soils from the state of Mississippi, USA, representing acidic, neutral, and alkaline pH, were exposed to different CO2 enrichments (100%, 10%, and 1%) for 30 days. The leachates’ pH tended to attain an equilibrium state (neutral) with time under CO2 saturation. SIC increased under CO2 saturation, whereas cation exchange capacity (CEC) showed a decreasing pattern in all soils. In the third chapter, an eXplainable artificial intelligence (XAI) was performed to visualize the different forms of soil carbon variability across the Mississippi River Basin area. This model explains key insights and local discrepancies, suggesting a solution to the “Black-Box” issue. The best performing model, stack ensemble, showed improved RMSE (3 to 8%) and spatial variability for soil carbons than other ML models, especially after adding the residuals from regression analyses. Land cover type > soil pH > total nitrogen, > NDVI were identified as the top four crucial factors for predicting SOC when bulk density > precipitation, soil pH > mean annual temperature described SIC. The proposed automatic machine learning (AutoML) model with model agnostic interpolations might be a hallmark to mitigate the C loss under adverse climate change conditions and allow diverse knowledge groups to adopt a new interpretable ML algorithm more confidently. Findings from this study help predict the impact of elevated atmospheric CO2 on soil pH, acidification, and nutrient availability and develop strategies for sustainable land management practices under a changing climate.

The objective of this thesis was to identify the mechanisms associated with the recurrence of rock-slope instability along the R518 and R523 roads in Limpopo. Advanced geotechnical methods and ASTER imagery were used for the purpose while a predictive rockfall hazard rating matrix chart and rock slope stability charts for unsaturated sensitive clay soil and rock slopes were to be developed. The influence of extreme rainfall on the slope stability of the sensitive clay soil was also evaluated. To achieve the above, field observations, geological mapping, kinematic analysis, and limit equilibrium were performed. The latter involved toppling, transitional and rotational analyses. Numerical simulation was finally resorted to. The following software packages were employed: SWEDGE, SLIDE, RocData, RocFall, DIPS, RocPlane, and Phase 2. The simulation outputs were analyzed in conjunction with ASTER images. The advanced remote sensing data paved the way for landslide susceptibility analysis. From all the above, rockfall hazard prediction charts and slope stability prediction charts were developed. Several factors were also shown by numerical simulation to influence slope instability in the area of study, i.e. sites along the R518 and R523 roads in the Thulamela Municipality. The most important factors are extreme rainfall, steep slopes, geological features and water streams in the region, and improper road construction. Owing to the complexity of the failure mechanisms in the study area, it was concluded that both slope stability prediction charts and rock hazard matrix charts are very useful. They indeed enable one to characterize slope instability in sensitive clay soils as well as rockfall hazards in the study area. It is however recommended that future work is undertaken to explore the use of sophisticated and scientific methods. This is instrumental in the development of predictive tools for rock deformation and displacement in landslide events. / Electrical and Mining Engineering / D. Phil. (Mining Engineering)

<p dir="ltr">Digital soil mapping (DSM) combines field and laboratory data with environmental factors to predict soil properties. The accuracy of these predictions depends on factors such as model selection, data quality and quantity, and landscape characteristics. In our study, we investigated the impact of sample density and the use of various environmental covariates (ECs) including slope, topographic position index, topographic wetness index, multiresolution valley bottom flatness, and multiresolution ridge top flatness, as well as the spatial resolution of these ECs on the predictive accuracy of four predictive models; Cubist (CB), Random Forest (RF), Regression Kriging (RK), and Ordinary Kriging (OK). Our analysis was conducted at three sites in Indiana: the Purdue Agronomy Center for Research and Education (ACRE), Davis Purdue Agriculture Center (DPAC), and Southeast Purdue Agricultural Center (SEPAC). Each site had its unique soil data sampling designs, management practices, and topographic conditions. The primary focus of this study was to predict the spatial distribution of soil properties, including soil organic matter (SOM), cation exchange capacity (CEC), and clay content, at different depths (0-10cm, 0-15cm, and 10-30cm) by utilizing five environmental covariates and four spatial resolutions for the ECs (1-1.5 m, 5 m, 10 m, and 30 m).</p><p dir="ltr">Various evaluation metrics, including R², root mean square error (RMSE), mean square error (MSE), concordance coefficient (pc), and bias, were used to assess prediction accuracy. Notably, the accuracy of predictions was found to be significantly influenced by the site, sample density, model type, soil property, and their interactions. Sites exhibited the largest source of variation, followed by sampling density and model type for predicted SOM, CEC, and clay spatial distribution across the landscape.</p><p dir="ltr">The study revealed that the RF model consistently outperformed other models, while OK performed poorly across all sites and properties as it only relies on interpolating between the points without incorporating the landscape characteristics (ECs) in the algorithm. Increasing sample density improved predictions up to a certain threshold (e.g., 66 samples at ACRE for both SOM and CEC; 58 samples for SOM and 68 samples for CEC at SEPAC), beyond which the improvements were marginal. Additionally, the study highlighted the importance of spatial resolution, with finer resolutions resulting in better prediction accuracy, especially for SOM and clay content. Overall, comparing data from the two depths (0-10cm vs 10-30cm) for soil properties predications, deeper soil layer data (10-30cm) provided more accurate predictions for SOM and clay while shallower depth data (0-10cm) provided more accurate predictions for CEC. Finally, higher spatial resolution of ECs such as 1-1.5 m and 5 m contributed to more accurate soil properties predictions compared to the coarser data of 10 m and 30 m resolutions.</p><p dir="ltr">In summary, this research underscores the significance of informed decisions regarding sample density, model selection, and spatial resolution in digital soil mapping. It emphasizes that the choice of predictive model is critical, with RF consistently delivering superior performance. These findings have important implications for land management and sustainable land use practices, particularly in heterogeneous landscapes and areas with varying management intensities.</p>

This research work concentrate on developing digital soil maps to support field based plant phenotyping research. We have developed soil organic matter content (OM), cation exchange capacity (CEC), natural soil drainage class, and tile drainage line maps using topographic indices and aerial imagery. Various prediction models (universal kriging, cubist, random forest, C5.0, artificial neural network, and multinomial logistic regression) were used to estimate the soil properties of interest.