Geospatial environmental data modelling applications using remote sensing, GIS and spatial statistics

This thesis presents novel modelling applications for environmental geospatial data using remote sensing, GIS and statistical modelling techniques. The studied themes can be classified into four main themes: (i) to develop advanced geospatial databases. Paper (I) demonstrates the creation of a geospatial database for the Glanville fritillary butterfly (Melitaea cinxia) in the Åland Islands, south-western Finland; (ii) to analyse species diversity and distribution using GIS techniques. Paper (II) presents a diversity and geographical distribution analysis for Scopulini moths at a world-wide scale; (iii) to study spatiotemporal forest cover change. Paper (III) presents a study of exotic and indigenous tree cover change detection in Taita Hills Kenya using airborne imagery and GIS analysis techniques; (iv) to explore predictive modelling techniques using geospatial data. In Paper (IV) human population occurrence and abundance in the Taita Hills highlands was predicted using the generalized additive modelling (GAM) technique. Paper (V) presents techniques to enhance fire prediction and burned area estimation at a regional scale in East Caprivi Namibia. Paper (VI) compares eight state-of-the-art predictive modelling methods to improve fire prediction, burned area estimation and fire risk mapping in East Caprivi Namibia. The results in Paper (I) showed that geospatial data can be managed effectively using advanced relational database management systems. Metapopulation data for Melitaea cinxia butterfly was successfully combined with GPS-delimited habitat patch information and climatic data. Using the geospatial database, spatial analyses were successfully conducted at habitat patch level or at more coarse analysis scales. Moreover, this study showed it appears evident that at a large-scale spatially correlated weather conditions are one of the primary causes of spatially correlated changes in Melitaea cinxia population sizes. In Paper (II) spatiotemporal characteristics of Socupulini moths description, diversity and distribution were analysed at a world-wide scale and for the first time GIS techniques were used for Scopulini moth geographical distribution analysis. This study revealed that Scopulini moths have a cosmopolitan distribution. The majority of the species have been described from the low latitudes, sub-Saharan Africa being the hot spot of species diversity. However, the taxonomical effort has been uneven among biogeographical regions. Paper III showed that forest cover change can be analysed in great detail using modern airborne imagery techniques and historical aerial photographs. However, when spatiotemporal forest cover change is studied care has to be taken in co-registration and image interpretation when historical black and white aerial photography is used. In Paper (IV) human population distribution and abundance could be modelled with fairly good results using geospatial predictors and non-Gaussian predictive modelling techniques. Moreover, land cover layer is not necessary needed as a predictor because first and second-order image texture measurements derived from satellite imagery had more power to explain the variation in dwelling unit occurrence and abundance. Paper V showed that generalized linear model (GLM) is a suitable technique for fire occurrence prediction and for burned area estimation. GLM based burned area estimations were found to be more superior than the existing MODIS burned area product (MCD45A1). However, spatial autocorrelation of fires has to be taken into account when using the GLM technique for fire occurrence prediction. Paper VI showed that novel statistical predictive modelling techniques can be used to improve fire prediction, burned area estimation and fire risk mapping at a regional scale. However, some noticeable variation between different predictive modelling techniques for fire occurrence prediction and burned area estimation existed.

Making the most of spatial information in health: A tutorial in Bayesian disease mapping for areal data

Disease maps are effective tools for explaining and predicting patterns of disease outcomes across geographical space, identifying areas of potentially elevated risk, and formulating and validating aetiological hypotheses for a disease. Bayesian models have become a standard approach to disease mapping in recent decades. This article aims to provide a basic understanding of the key concepts involved in Bayesian disease mapping methods for areal data. It is anticipated that this will help in interpretation of published maps, and provide a useful starting point for anyone interested in running disease mapping methods for areal data. The article provides detailed motivation and descriptions on disease mapping methods by explaining the concepts, defining the technical terms, and illustrating the utility of disease mapping for epidemiological research by demonstrating various ways of visualising model outputs using a case study. The target audience includes spatial scientists in health and other fields, policy or decision makers, health geographers, spatial analysts, public health professionals, and epidemiologists.

Efficient search for statistically significant dependency rules in binary data

Analyzing statistical dependencies is a fundamental problem in all empirical science. Dependencies help us understand causes and effects, create new scientific theories, and invent cures to problems. Nowadays, large amounts of data is available, but efficient computational tools for analyzing the data are missing. In this research, we develop efficient algorithms for a commonly occurring search problem - searching for the statistically most significant dependency rules in binary data. We consider dependency rules of the form X->A or X->not A, where X is a set of positive-valued attributes and A is a single attribute. Such rules describe which factors either increase or decrease the probability of the consequent A. A classical example are genetic and environmental factors, which can either cause or prevent a disease. The emphasis in this research is that the discovered dependencies should be genuine - i.e. they should also hold in future data. This is an important distinction from the traditional association rules, which - in spite of their name and a similar appearance to dependency rules - do not necessarily represent statistical dependencies at all or represent only spurious connections, which occur by chance. Therefore, the principal objective is to search for the rules with statistical significance measures. Another important objective is to search for only non-redundant rules, which express the real causes of dependence, without any occasional extra factors. The extra factors do not add any new information on the dependence, but can only blur it and make it less accurate in future data. The problem is computationally very demanding, because the number of all possible rules increases exponentially with the number of attributes. In addition, neither the statistical dependency nor the statistical significance are monotonic properties, which means that the traditional pruning techniques do not work. As a solution, we first derive the mathematical basis for pruning the search space with any well-behaving statistical significance measures. The mathematical theory is complemented by a new algorithmic invention, which enables an efficient search without any heuristic restrictions. The resulting algorithm can be used to search for both positive and negative dependencies with any commonly used statistical measures, like Fisher's exact test, the chi-squared measure, mutual information, and z scores. According to our experiments, the algorithm is well-scalable, especially with Fisher's exact test. It can easily handle even the densest data sets with 10000-20000 attributes. Still, the results are globally optimal, which is a remarkable improvement over the existing solutions. In practice, this means that the user does not have to worry whether the dependencies hold in future data or if the data still contains better, but undiscovered dependencies.

Optimal adaptive sequential spatial sampling of soil using pair-copulas

A spatial sampling design that uses pair-copulas is presented that aims to reduce prediction uncertainty by selecting additional sampling locations based on both the spatial configuration of existing locations and the values of the observations at those locations. The novelty of the approach arises in the use of pair-copulas to estimate uncertainty at unsampled locations. Spatial pair-copulas are able to more accurately capture spatial dependence compared to other types of spatial copula models. Additionally, unlike traditional kriging variance, uncertainty estimates from the pair-copula account for influence from measurement values and not just the configuration of observations. This feature is beneficial, for example, for more accurate identification of soil contamination zones where high contamination measurements are located near measurements of varying contamination. The proposed design methodology is applied to a soil contamination example from the Swiss Jura region. A partial redesign of the original sampling configuration demonstrates the potential of the proposed methodology.

Valence and spatial explanations for voting in the 2013 Australian election

This paper examines the 2013 Australian federal election to test two competing models of vote choice: spatial politics and valence issues. Using data from the 2013 Australian Election Study, the analysis finds that spatial politics (measured by party identification and self-placement on the left-right spectrum) and valence issues both have significant effects on vote choice. However, spatial measures are more important than valence issues in explaining vote choice, in contrast with recent studies from Britain, Canada and the United States. Explanations for these differences are speculative, but may relate to Australia’s stable party and electoral system, including compulsory voting and the frequency of elections. The consequently high information burden faced by Australian voters may lead to a greater reliance on spatial heuristics than is found elsewhere.

Recapitulation in generating spatial layouts: Representations of embryogenesis in evolutionary design modelling

The noted 19th century biologist, Ernst Haeckel, put forward the idea that the growth (ontogenesis) of an organism recapitulated the history of its evolutionary development. While this idea is defunct within biology, the idea has been promoted in areas such as education (the idea of an education being the repetition of the civilizations before). In the research presented in this paper, recapitulation is used as a metaphor within computer-aided design as a way of grouping together different generations of spatial layouts. In most CAD programs, a spatial layout is represented as a series of objects (lines, or boundary representations) that stand in as walls. The relationships between spaces are not usually explicitly stated. A representation using Lindenmayer Systems (originally designed for the purpose of modelling plant morphology) is put forward as a way of representing the morphology of a spatial layout. The aim of this research is not just to describe an individual layout, but to find representations that link together lineages of development. This representation can be used in generative design as a way of creating more meaningful layouts which have particular characteristics. The use of genetic operators (mutation and crossover) is also considered, making this representation suitable for use with genetic algorithms.

Spatial and temporal determinants of Finnish farmland bird populations

During the past decades agricultural intensification has caused dramatic population declines in a wide range of taxa related to farmland habitats, including farmland birds. In this thesis, I studied how boreal farmland landscape characteristics and agricultural land use affect the abundance and diversity of farmland birds using extensive field data collected by territory mapping of breeding farmland birds in various parts of Finland. My results show that the area and openness of agricultural areas are key determinants of farmland bird abundance and distribution. A landscape composition with enough open farmland combined with key habitats such as farmyards and wetland is likely to provide essential prerequisites for the occurrence of a rich farmland avifauna. In Finland, the majority of large areas suitable for open habitat specialists are located in southern and western parts of the country. However, the diversity of the species with an unfavourable conservation status in Europe (SPECs) had notable hotspot areas in northern and north-western agricultural areas. I found that in boreal agroecosystems farmland birds favour fields with springtime vegetative cover, especially agricultural grasslands and set-asides. Hence, in the spring cereal dominated Finnish agroecosystems it is the absence of field vegetation that may limit populations of many farmland bird species. It is likely that the decrease of crops providing vegetative cover in the spring, such as permanent grasslands, cultivated grass, and autumn-sown cereals, has greatly contributed to the declines of Finnish farmland birds. Grass crops have persistently declined in Finland as a consequence of specialization in crop production and the large-scale decline in livestock husbandry. Small-scale non-crop habitats, especially ditches and ditch margins, are also important for many bird species in the Finnish agroecosystems, but have dramatically declined during the last decades. A major problem for farmland bird conservation in Finland is the conflict between landscape structure and agricultural management. Areas with mixed and cattle farming are virtually absent from the large agricultural plains of southern and south-western Finland, where the landscape structure is more likely to be favourable for rich farmland bird assemblages. On the other hand, mixed and cattle farming is still rather frequent in northern and central parts of the country, where the landscape structure is not suitable for many farmland specialist birds requiring open landscapes. My results provide useful guidelines for farmland bird conservation, and imply that considerable attention needs to be paid to landscape factors when selecting areas for various conservational management actions, such as agri-environment schemes. Actions promoting the abundance of set-asides, grass crops, and ditches would markedly benefit Finnish farmland bird populations. Organic farming may benefit farmland birds, but it is not clear how general its beneficial effect is in boreal agroecosystems. The most urgent action aiming to preserve farmland biodiversity would be to support re-introducing and sustaining cattle farming by environmental subsidies. This would be especially beneficial in the southern parts of Finland, where the landscape characteristics and abundance of agricultural areas are most suitable for farmland birds and where cattle farming is currently rare.

Where and how to conserve: Extending the scope of spatial reserve network design

Ongoing habitat loss and fragmentation threaten much of the biodiversity that we know today. As such, conservation efforts are required if we want to protect biodiversity. Conservation budgets are typically tight, making the cost-effective selection of protected areas difficult. Therefore, reserve design methods have been developed to identify sets of sites, that together represent the species of conservation interest in a cost-effective manner. To be able to select reserve networks, data on species distributions is needed. Such data is often incomplete, but species habitat distribution models (SHDMs) can be used to link the occurrence of the species at the surveyed sites to the environmental conditions at these locations (e.g. climatic, vegetation and soil conditions). The probability of the species occurring at unvisited location is next predicted by the model, based on the environmental conditions of those sites. The spatial configuration of reserve networks is important, because habitat loss around reserves can influence the persistence of species inside the network. Since species differ in their requirements for network configuration, the spatial cohesion of networks needs to be species-specific. A way to account for species-specific requirements is to use spatial variables in SHDMs. Spatial SHDMs allow the evaluation of the effect of reserve network configuration on the probability of occurrence of the species inside the network. Even though reserves are important for conservation, they are not the only option available to conservation planners. To enhance or maintain habitat quality, restoration or maintenance measures are sometimes required. As a result, the number of conservation options per site increases. Currently available reserve selection tools do however not offer the ability to handle multiple, alternative options per site. This thesis extends the existing methodology for reserve design, by offering methods to identify cost-effective conservation planning solutions when multiple, alternative conservation options are available per site. Although restoration and maintenance measures are beneficial to certain species, they can be harmful to other species with different requirements. This introduces trade-offs between species when identifying which conservation action is best applied to which site. The thesis describes how the strength of such trade-offs can be identified, which is useful for assessing consequences of conservation decisions regarding species priorities and budget. Furthermore, the results of the thesis indicate that spatial SHDMs can be successfully used to account for species-specific requirements for spatial cohesion - in the reserve selection (single-option) context as well as in the multi-option context. Accounting for the spatial requirements of multiple species and allowing for several conservation options is however complicated, due to trade-offs in species requirements. It is also shown that spatial SHDMs can be successfully used for gaining information on factors that drive a species spatial distribution. Such information is valuable to conservation planning, as better knowledge on species requirements facilitates the design of networks for species persistence. This methods and results described in this thesis aim to improve species probabilities of persistence, by taking better account of species habitat and spatial requirements. Many real-world conservation planning problems are characterised by a variety of conservation options related to protection, restoration and maintenance of habitat. Planning tools therefore need to be able to incorporate multiple conservation options per site, in order to continue the search for cost-effective conservation planning solutions. Simultaneously, the spatial requirements of species need to be considered. The methods described in this thesis offer a starting point for combining these two relevant aspects of conservation planning.

Integrating urban ecosystem sustainability assessment into policy-making: Insights from the Gold Coast City

This paper introduces a policy-making support tool called ‘Micro-level Urban ecosystem Sustainability IndeX (MUSIX)’. The index serves as a sustainability assessment model that monitors six aspects of urban ecosystems, hydrology, ecology, pollution, location, design, and efficiency based on parcel-scale indicators. This index is applied in a case study investigation in the Gold Coast City, Queensland, Australia. The outcomes reveal that there are major environmental problems caused by increased impervious surfaces from growing urban development in the study area. The findings suggest that increased impervious surfaces are linked to increased surface runoff, car dependency, transport-related pollution, poor public transport accessibility, and unsustainable built environment. This paper presents how the MUSIX outputs can be used to guide policy-making through the evaluation of existing policies.

Spatial patterning of the distribution of Ca-2+ in Dictyostelium discoideum

We have shown previously that the Ca2+-specific fluorescent dyes chlortetracycline (CTC) and indo-1/AM can be used to distinguish between prestalk and prespore cells in Dictyostelium discoideum at a very early stage. In the present study, pre- and post-aggregative amoebae of Dictyostelium discoideum were labelled with CTC or indo-1 and their fluorescence monitored after being drawn into a fine glass capillary. The cells rapidly form two zones of Ca2+-CTC or Ca2+-indo-1 fluorescence. Anterior (air side) cells display a high level of fluorescence; the level drops in the middle portion of the capillary and rises again to a lesser extent in the posteriormost cells (oil side). When bounded by air on both sides, the cells display high fluorescence at both ends. When oil is present at both ends of the capillary, there is little fluorescence except for small regions at the ends. These outcomes are evident within a couple of minutes of the start of the experiment and the fluorescence pattern intensifies over the course of time. By using the indicator neutral red, as well as with CTC and indo-1, we show that a band displaying strong fluorescence moves away from the anterior end before stabilizing at the anterior-posterior boundary. We discuss our findings in relation to the role of Ca2+ in cell-type differentiation in Dictyostelium discoideum.

Urban Growth Analysis Using Spatial and Temporal Data

Urban growth identification, quantification, knowledge of rate and the trends of growth would help in regional planning for better infrastructure provision in environmentally sound way. This requires analysis of spatial and temporal data, which help in quantifying the trends of growth on spatial scale. Emerging technologies such as Remote Sensing, Geographic Information System (GIS) along with Global Positioning System (GPS) help in this regard. Remote sensing aids in the collection of temporal data and GIS helps in spatial analysis. This paper focuses on the analysis of urban growth pattern in the form of either radial or linear sprawl along the Bangalore - Mysore highway. Various GIS base layers such as builtup areas along the highway, road network, village boundary etc. were generated using collateral data such as the Survey of India toposheet, etc. Further, this analysis was complemented with the computation of Shannon's entropy, which helped in identifying prevalent sprawl zone, rate of growth and in delineating potential sprawl locations. The computation Shannon's entropy helped in delineating regions with dispersed and compact growth. This study reveals that the Bangalore North and South taluks contributed mainly to the sprawl with 559% increase in built-up area over a period of 28 years and high degree of dispersion. The Mysore and Srirangapatna region showed 128% change in built-up area and a high potential for sprawl with slightly high dispersion. The degree of sprawl was found to be directly proportional to the distances from the cities.

Spatial variation of sequestered calcium in the multicellular stage of Dictyostelium discoideum as assayed by chlortetracycline fluorescence

Abstract. We have used chlortetracycline (CTC) as a fluorescent probe to detect the distribution of sequestered calcium in multicellular stages of Dictyostelium discoideum. Tips of late aggregates, slugs and early culminating masses fluoresce very strongly. Most of the fluorescence is intracellular in origin and emanates from a small number of intense punctate sources. The sources correspond in part to autophagic vacuoles viz. neutral-red staining, acidic digestive vesicles, and may also include intracellular organelles; cytoplasmic fluorescence is much weaker in comparison. The level of fluorescence drops in the middle portion of slugs and rises again in the posteriormost region, though not to as high a level as in the tip. This holds good irrespective of whether CTC is applied only in the neighbourhood of the aggregate centre, only in the aggregate periphery, or to the whole aggregate. We infer that there must be a good deal of mixing in the stages leading from aggregation to slug formation; thus the serial order in which cells enter an aggregate does not bear any relation to their ultimate fates. The other implication of our study is that calcium sequestration is much more extensive in prestalk and anterior-like cells than in prespore cells. These findings are discussed with regard to possible implications for pattern formation.