Densification of spatially-sparse legacy soil data at a national scale: a digital mapping approach

Digital soil mapping (DSM) is a viable approach to providing spatial soil information but its adoption at the national scale, especially in sub-Saharan Africa, is limited by low spread of data. Therefore, the focus of this thesis is on optimizing DSM techniques for densification of sparse legacy soil data using Nigeria as a case study. First, the robustness of Random Forest model (RFM) was tested in predicting soil particle-size fractions as a compositional data using additive log-ratio technique. Results indicated good prediction accuracy with RFM while soils are largely coarse-textured especially in the northern region. Second, soil organic carbon (SOC) and bulk density (BD) were predicted from which SOC density and stock were calculated. These were overlaid with land use/land cover (LULC), agro-ecological zone (AEZ) and soil maps to quantify the carbon sequestration of soils and their variation across different AEZs. Results showed that 6.5 Pg C with an average of 71.60 Mg C ha–1 abound in the top 1 m soil depth. Furthermore, to improve the performance of BD and effective cation exchange capacity (ECEC) pedotransfer functions (PTFs), the inclusion of environmental data was explored using multiple linear regression (MLR) and RFM. Results showed an increase in performance of PTFs with the use of soil and environmental data. Finally, the application of Choquet fuzzy integral (CI) technique in irrigation suitability assessment was assessed. This was achieved through multi-criteria analysis of soil, climatic, landscape and socio-economic indices. Results showed that CI is a better aggregation operator compared to weighted mean technique. A total of 3.34 x 106 ha is suitable for surface irrigation in Nigeria while major limitations are due to topographic and soil attributes. Research findings will provide quantitative basis for framing appropriate policies on sustainable food production and environmental management, especially in resource-poor countries of the world

Mapping three-dimensional geological features from remotely-sensed images and digital elevation models.

Accurate mapping of geological structures is important in numerous applications, ranging from mineral exploration through to hydrogeological modelling. Remotely sensed data can provide synoptic views of study areas enabling mapping of geological units within the area. Structural information may be derived from such data using standard manual photo-geologic interpretation techniques, although these are often inaccurate and incomplete. The aim of this thesis is, therefore, to compile a suite of automated and interactive computer-based analysis routines, designed to help a the user map geological structure. These are examined and integrated in the context of an expert system. The data used in this study include Digital Elevation Model (DEM) and Airborne Thematic Mapper images, both with a spatial resolution of 5m, for a 5 x 5 km area surrounding Llyn Cow lyd, Snowdonia, North Wales. The geology of this area comprises folded and faulted Ordo vician sediments intruded throughout by dolerite sills, providing a stringent test for the automated and semi-automated procedures. The DEM is used to highlight geomorphological features which may represent surface expressions of the sub-surface geology. The DEM is created from digitized contours, for which kriging is found to provide the best interpolation routine, based on a number of quantitative measures. Lambertian shading and the creation of slope and change of slope datasets are shown to provide the most successful enhancement of DEMs, in terms of highlighting a range of key geomorphological features. The digital image data are used to identify rock outcrops as well as lithologically controlled features in the land cover. To this end, a series of standard spectral enhancements of the images is examined. In this respect, the least correlated 3 band composite and a principal component composite are shown to give the best visual discrimination of geological and vegetation cover types. Automatic edge detection (followed by line thinning and extraction) and manual interpretation techniques are used to identify a set of 'geological primitives' (linear or arc features representing lithological boundaries) within these data. Inclusion of the DEM data provides the three-dimensional co-ordinates of these primitives enabling a least-squares fit to be employed to calculate dip and strike values, based, initially, on the assumption of a simple, linearly dipping structural model. A very large number of scene 'primitives' is identified using these procedures, only some of which have geological significance. Knowledge-based rules are therefore used to identify the relevant. For example, rules are developed to identify lake edges, forest boundaries, forest tracks, rock-vegetation boundaries, and areas of geomorphological interest. Confidence in the geological significance of some of the geological primitives is increased where they are found independently in both the DEM and remotely sensed data. The dip and strike values derived in this way are compared to information taken from the published geological map for this area, as well as measurements taken in the field. Many results are shown to correspond closely to those taken from the map and in the field, with an error of < 1°. These data and rules are incorporated into an expert system which, initially, produces a simple model of the geological structure. The system also provides a graphical user interface for manual control and interpretation, where necessary. Although the system currently only allows a relatively simple structural model (linearly dipping with faulting), in the future it will be possible to extend the system to model more complex features, such as anticlines, synclines, thrusts, nappes, and igneous intrusions

Monitoring, Modelling and Management of Water Quality

Different types of pressures, such as nutrients, micropollutants, microbes, nanoparticles, microplastics, or antibiotic-resistant genes, endanger the quality of water bodies. Evidence-based pollution control needs to be built on the three basic elements of water governance: Monitoring, modeling, and management. Monitoring sets the empirical basis by providing space- and time-dependent information on substance concentrations and loads, as well as driving boundary conditions for assessing water quality trends, water quality statuses, and providing necessary information for the calibration and validation of models. Modeling needs proper system understanding and helps to derive information for times and locations where no monitoring is done or possible. Possible applications are risk assessments for exceedance of quality standards, assessment of regionalized relevance of sources and pathways of pollution, effectiveness of measures, bundles of measures or policies, and assessment of future developments as scenarios or forecasts. Management relies on this information and translates it in a socioeconomic context into specific plans for implementation. Evaluation of success of management plans again includes well-defined monitoring strategies. This book provides an important overview in this context

Creation of a gridded time series of hydrological variables for Canada

There is a lack of measured, long-term, reliable, and well-distributed hydrological variables in Canada. These hydrological variables include, but are not limited to: temperature, precipitation, ground runoff, evapotranspiration, soil moisture, and snow water equivalent. The objective of this thesis was to establish the best possible distributed estimates of these hydrological variables for Canada over the period of 1961-2000. The first step was to interpolate measured temperature and precipitation across the country. These interpolated values were then used to calculate the other hydrological variables using the Waterloo Flood Forecasting Model (WATFLOOD). The Waterloo Mapping technique (WATMAP) was developed to use topographic and land cover databases to automatically and systematically derive the information needed to create the drainage database. WATFLOOD was calibrated with the Dynamically Dimensioned Search (DDS) algorithm using the difference between the measured and simulated streamflow as the objective function. After a final calibration of 100 separate DDS runs, distributed time series for the hydrological variables were created. A simple assessment was made for the predictive uncertainty in the simulated streamflow results based on the results of the final calibration. As well, the implications of various climate change scenarios were examined in the context of how they would change the hydrological variables. The major recommendations for future study included: finding other gridded datasets that could be used to verify the ones that were created in this study and examining further the magnitudes of the different kinds of predictive uncertainty (data, model, and parameter). The results of this thesis fit in well with the goals of the study on Predictions in Ungauged Basins. This thesis was organized along the principle of “design the process, not the product”. As such, although a set of final products are presented at the end, the most important part of the thesis was the process that achieved these products. Thus it is not assumed that every technique designed in this thesis will be applicable to every other researcher, but it hoped that most researchers in the field will be able to use at least some parts of the techniques developed here

The occurrence and origin of salinity in non-coastal groundwater in the Waikato region

Aims The aims of this project are to describe the occurrence, and determine the origin of non-coastal saline groundwater in the Waikato region. High salinity limits the use of the water for supply and agricultural use. Understanding the origin and distribution of non-coastal salinity will assist with development and management of groundwater resources in the Waikato. Method The occurrence of non-coastal groundwater salinity was investigated by examining driller’s records and regional council groundwater quality information. Selected wells were sampled for water quality analyses and temperatures were profiled where possible. Water quality analyses include halogens such as chloride, fluoride, iodide and bromide. Ratios of these ions are useful to differentiate between geothermal and seawater origins of salinity (Hem, 1992). Other ionic ratio approaches for differentiating sources and influences on salinity such as those developed by Alcala and Emilio (2008) and Sanchez-Martos et al., (2002), may also be applied. Potential sources of salinity include seawater, connate water, geothermal and anthropogenic influences. The hydrogeologic settings of saline occurrence were also investigated, to explore the potential to predict further occurrence. Results Numerous occurrences of non-coastal saline groundwater have been observed in the Waikato region. Where possible, wells with relatively high total dissolved solids (TDS) were selected for further investigation. Several groundwater samples are moderately saline and exceed the TDS drinking water aesthetic guideline of 1,000 g m-3 (Ministry of Health, 2008). Selected ion ratios (predominantly halogens) were used to assist in differentiating between influences on salinity such as seawater and geothermal. Bromide to iodide ratios, in particular, infer a greater geothermal influence on salinity, although other ratios are not definitive. The anomalously elevated salinity observed appears natural but nevertheless has constrained localised groundwater resource development for dairy factory, industrial and prison water supply use. Further work may show some relationship with geology or tectonics, which could assist prediction of inland saline groundwater occurrence

Modelling the Holocene evolution of coastal gullies on the Isle of Wight

Geomorphological evidence has frequently been used to infer past environmental conditions, but in recent years the emergence of landscape evolution models (LEMs) has opened the possibility of using numerical modelling as a tool in palaeo-environmental reconstruction. The application of LEMs for this purpose involves retrodictive modelling, each simulation scenario being configured with model variables (e.g. reflecting climate change) and parameters to reflect a specific hypothesis of environmental change. Plausible scenarios are then identified by matching contemporary observed and modelled landscapes. However, although considerable uncertainty is known to surround the specification of model driving conditions and parameters, previous studies have not considered this issue. This research applies a technique of accounting for the uncertainty surrounding the specification of driving conditions and model parameters by using reduced complexity 'metamodels' to analyse the full model parameter space and thus constrain sources of uncertainty and plausible retrodicted scenarios more effectively. This study applies the developed techniques to a case study focused on a specific set of coastal gullies found on the Isle of Wight, UK. A key factor in the evolution of these gullies are the relative balance between rates of cliff retreat (which reduces gully extent) and headwards incision caused by knickpoint migration (which increases gully extent). To inform the choice and parameterisation of the numerical model used in this research an empirical-conceptual model of gully evolution was initially developed. To provide a long-term context for the evolution of the gullies and to identify the relative importance of the various driving factors, the Holocene erosional history of the Isle of Wight gullies was then simulated using a LEM. In a preliminary set of simulations a 'traditional' (i.e. with no consideration of parameter uncertainty) retrodictive modelling approach was applied, in which driving variables were arbitrarily altered and observed and simulated landscape topographies compared, under various scenarios of imposed environmental change. These initial results revealed that the coastal gullies have been ephemeral in nature for much of the Holocene, only becoming semi-permanent once cliff retreat rates fall below a critical threshold at 2500 cal. years BP. Next, in an attempt to constrain more detailed erosional histories and to explore the extent to which retrodicted interpretations of landscape change were confounded by uncertainty, a Central Composite Design (CCD) sampling technique was employed to sample variations in the model driving variables, enabling the trajectories of gully response to different combinations of the driving conditions to be modelled explicitly. In some of these simulations, where the ranges of bedrock erodibility (0:03 - 0:04m0:2a)

Remote sensing of energy and water fluxes over Volta Savannah catchments in West Africa

The deterioration of the West African savannah in the last three decades is believed to be closely linked with about 0.5 C rise in temperature leading to evaporation losses and declining levels of the Volta Lake in Ghana. Although hydrological models can be used to predict climate change impacts on the regional hydrology, spatially-observed ground data needed for this purpose are largely unavailable. This thesis seeks to address this problem by developing improved methods for estimating energy and water fluxes (e.g. latent heat [ET]) from remotely sensed data and to demonstrate how these may be used to parameterize hydrological models. The first part of the thesis examines the potential of the Penman-Monteith method to estimate local-scale ET using groundbased hydrometeorological observations, vegetation coefficients and environmental data. The model results were compared with pan observations, scintillometer (eddy correlation) measurements and the Thomthwaite empirical method. The Penman- Monteith model produced better evaporation estimates (~3.90 mm day(^-1) for the Tamale district) than its counterpart methods. The Thomthwaite, for example, overestimated predictions by 5.0-11.0 mm day(^-1). Up-scaling on a monthly time scale and parameterization of the Grindley soil moisture balance model with the Thomthwaite and Penman-Monteith data, however, produced similar estimates of actual evaporation and soil moisture, which correlated strongly (R(^2) = 0.95) with water balance estimates. To improve ET estimation at the regional-scale, the second part of the thesis develops spatial models through energy balance modelling and data up-scaling methods, driven by radiometric measurements from recent satellite sensors such as the Landsat ETM+, MODIS and ENVISAT-AATSR. The results were validated using estimates from the Penman-Monteith method, field observations, detailed satellite measurements and published data. It was realised that the MODIS sensor is a more useful source of energy and water balance parameters than AA TSR. For example, stronger correlations were found between MODIS estimates of ET and other energy balance variables such as NDVI, surface temperature and net radiation (R(^2) = 0.67-0.73) compared with AATSR estimates (R(^2) = 0.31-0.40). There was also a good spatial correlation between MODIS and Landsat ETM+ results (R(^2) = 0.71), but poor correlations were found between AATSR and Landsat data (R(^2) = 0.0-0.13), which may be explained by differences in instrument calibration. The results further showed that ET may be underestimated with deviations of ~2.0 mm day 1 when MODIS/AATSR measurements are validated against point observations because of spatial mismatch. The final part of the thesis demonstrates the application of the ET model for predicting runoff (Q) using a simplified version of the regional water balance equation. This is followed byanalysis of flow sensitivity to declining scenarios of biomass volume. The results showed the absence of Q for >90% of the study area during the dry season due largely to crude model approximation and lack of rainfall data, which makes model testing during the wet season important. Runoff prediction may be improved if spatial estimates of rainfall, ET and geographical data (e.g. land-use/cover maps, soil & geology maps and DEM) could be routinely derived from satellite imagery

Coupling reactive transport and travel time modeling at the watershed scale

Nonpoint source pollution poses the greatest threat to water quality in developed countries. Modeling this type of pollution is a challenge for reactive transport models because of the change in scale: moving from a local field site- to a watershed-size problem. Computational resources and detailed watershed characterization are the major limiting factors in the fully time- and space-resolved modeling of the subsurface fate and transport of pollutants from nonpoint sources. While detailed characterization has been performed on a few well-studied watersheds, the knowledge derived from these watersheds has not led to a better understanding of watershed functioning in ungauged watersheds. Consequently, alternative approaches to modeling subsurface nonpoint source pollution have emerged to inform risk assessment to water resources and watershed management. This study investigates the development of a methodology that decouples flow and transport with the implementation of an analytical approach for 1-D travel time probability distribution functions (PDFs) to simulate subsurface flow at the watershed scale, that is, a 3-D problem. The first two chapters of my thesis focus on constraining and providing tools for the implementation of this methodology in watersheds. First, the analytical methodology for travel time was tested under varying conditions of heterogeneity, slope, and aquifer depths that were imposed on a virtual watershed, using Alder Creek, Ontario, as a test case. The analytical method parameters for the 28 scenarios considered were calibrated against the travel time PDFs generated with a 3-D numerical model (FEFLOW), which was used as baseline for comparison. The analytical method simulations revealed a negative relationship between the watershed mean travel time (wMTT) and the degree of imposed heterogeneity (σ_Y^2) of geostatistically defined permeability fields. This relationship was attributed to the effect of preferential flow paths. The effect of increasing aquifer depth (i.e., bedrock topography) on wMTT was similar to that of reducing the slope in surface topography, both resulting in an increase in wMTT. Given the promising results of the analytical method in the Alder Creek virtual analogs, further testing was conducted in 8 additional virtual watersheds. This inter-watershed comparison study examined the effects of 28 geomorphological indexes on wMTT and their predictive power in estimating analytical model parameters. This study is the first inter-watershed comparison of subsurface models that establishes relationships between watershed features and hydrologic functioning for groundwater storage and discharge. Among the classes of watershed features considered, those related to elevation (e.g. Relief), texture topography (e.g. drainage density, Dd), and Horton’s law (e.g. bifurcation factor, RB) were the most influential geomorphological classes emerging in the developed regression models. These regression models enable the application of the analytical methodology for deriving travel time PDF in other environmental settings. The transferability of these tools was verified for three extra watersheds in which the particle median travel time (pMTT), and their travel time distribution (TTD) performed on par to the upper tier of the original watersheds. Further research is proposed to include subsurface heterogeneity in the analysis to better evaluate its role in regulating wMTT in a subset of these watersheds. This methodology may constitute in a viable modeling alternative where subsurface information is scarce or scale limitations exist in developing a subsurface numerical model. The analytical methodology can provide a first line of knowledge in subsurface travel time and its distribution in an ungauged basin through the use of readily available tools (i.e., GIS and MATLAB). This knowledge can be later challenged or verified as more information becomes available. Potential directions to explore for the improvement of the methodology are proposed for further research. The third chapter applies the travel time PDF approach to the allocation of nitrogen (N) fluxes from base flow contributions to stream water chemistry in an existing hydrological model of Carroll Creek (Grand River basin, Ontario). This is a prospective chapter in which an outline for the development of an N isotope model linked to a hydrological model is presented. The N isotope model includes relevant N transformation and 15N fractionation processes in the plant-soil system and aims at simulating N-NO3- concentrations and isotopic compositions (δ15N). A bottom-up, stepwise approach is proposed in order to determine the most essential 15N discriminating processes and spatial discretization required by the model to match observations in the watershed

Geomorphometry 2020. Conference Proceedings

Geomorphometry is the science of quantitative land surface analysis. It gathers various mathematical, statistical and image processing techniques to quantify morphological, hydrological, ecological and other aspects of a land surface. Common synonyms for geomorphometry are geomorphological analysis, terrain morphometry or terrain analysis and land surface analysis. The typical input to geomorphometric analysis is a square-grid representation of the land surface: a digital elevation (or land surface) model. The first Geomorphometry conference dates back to 2009 and it took place in Zürich, Switzerland. Subsequent events were in Redlands (California), Nánjīng (China), Poznan (Poland) and Boulder (Colorado), at about two years intervals. The International Society for Geomorphometry (ISG) and the Organizing Committee scheduled the sixth Geomorphometry conference in Perugia, Italy, June 2020. Worldwide safety measures dictated the event could not be held in presence, and we excluded the possibility to hold the conference remotely. Thus, we postponed the event by one year - it will be organized in June 2021, in Perugia, hosted by the Research Institute for Geo-Hydrological Protection of the Italian National Research Council (CNR IRPI) and the Department of Physics and Geology of the University of Perugia. One of the reasons why we postponed the conference, instead of canceling, was the encouraging number of submitted abstracts. Abstracts are actually short papers consisting of four pages, including figures and references, and they were peer-reviewed by the Scientific Committee of the conference. This book is a collection of the contributions revised by the authors after peer review. We grouped them in seven classes, as follows: • Data and methods (13 abstracts) • Geoheritage (6 abstracts) • Glacial processes (4 abstracts) • LIDAR and high resolution data (8 abstracts) • Morphotectonics (8 abstracts) • Natural hazards (12 abstracts) • Soil erosion and fluvial processes (16 abstracts) The 67 abstracts represent 80% of the initial contributions. The remaining ones were either not accepted after peer review or withdrawn by their Authors. Most of the contributions contain original material, and an extended version of a subset of them will be included in a special issue of a regular journal publication

Distributed Smart City Services for Urban Ecosystems

A Smart City is a high-performance urban context, where citizens live independently and are more aware of the surrounding opportunities, thanks to forward-looking development of economy politics, governance, mobility and environment. ICT infrastructures play a key-role in this new research field being also a mean for society to allow new ideas to prosper and new, more efficient approaches to be developed. The aim of this work is to research and develop novel solutions, here called smart services, in order to solve several upcoming problems and known issues in urban areas and more in general in the modern society context. A specific focus is posed on smart governance and on privacy issues which have been arisen in the cellular age