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.

Validation of a new spot form of net blotch differential set and evidence for hybridisation between the spot and net forms of net blotch in Australia

A recently developed spot form of blotch differential set of 16 barley lines was tested for reaction response to 60 Pyrenophora teres f. maculata isolates from geographically disperse barley crops of Australia. Twelve barley lines (Arimont, Barque, Chebec, CI5286, CI5791, CI9214, CII6150, Dairokkaku, Esperance Orge 289, Galleon, Keel, Skiff, Torrens and TR250) provided differential response between the isolates. The susceptible controls Gairdner and Kombar provided indication of isolate virulence or avirulence. Abundant pathogenic diversity was revealed with 33 designated pathotypes, some of which related to geographic region. AFLP analysis also revealed abundant diversity with each of the isolates representing a unique genotype and one isolate that contained both AFLP bands unique to P. teres f. maculata and P. teres f. teres, the cause of spot form and net form of net blotch respectively, suggesting that sexual recombination between the net form and spot form isolates may have occurred naturally in the field.

Spot form of net blotch resistance in barley is under complex genetic control

Key message: Evaluation of resistance toPyrenophora teresf.maculatain barley breeding populations via association mapping revealed a complex genetic architecture comprising a mixture of major and minor effect genes. Abstract: In the search for stable resistance to spot form of net blotch (Pyrenophora teres f. maculata, SFNB), association mapping was conducted on four independent barley (Hordeum vulgare L.) breeding populations comprising a total of 898 unique elite breeding lines from the Northern Region Barley Breeding Program in Australia for discovery of quantitative trait loci (QTL) influencing resistance at seedling and adult plant growth stages. A total of 29 significant QTL were validated across multiple breeding populations, with 22 conferring resistance at both seedling and adult plant growth stages. The remaining 7 QTL conferred resistance at either seedling (2 QTL) or adult plant (5 QTL) growth stages only. These 29 QTL represented 24 unique genomic regions, of which five were found to co-locate with previously identified QTL for SFNB. The results indicated that SFNB resistance is controlled by a large number of QTL varying in effect size with large effects QTL on chromosome 7H. A large proportion of the QTL acted in the same direction for both seedling and adult responses, suggesting that phenotypic selection for SFNB resistance performed at either growth stage could achieve adequate levels of resistance. However, the accumulation of specific resistance alleles on several chromosomes must be considered in molecular breeding selection strategies.

Modified Atmosphere Packaging (MAP) for Control of Black Spot Formation in Chilled Prawns

The marketing of organically labeled prawns is predominately in a cooked or raw frozen form to avoid the development of melanosis (black spot). Certification for organic status prohibits the use of any added chemicals. The application of 60% CO2/40%N2 modified atmosphere to chilled (raw) prawns using two species of prawn was investigated for the ability to control black spot formation. Sensory assessment and microbiological counts were used to determine the end of product shelf life. Modified atmosphere packaged (MAP) prawns exhibited no melanosis for up to 16 days. The high quality life was retained for 12 days; shelf life of 16 days, according to standard microbiological criteria, was achieved, which is more than twice previously reported for non-MAP prawns. Results suggest MAP may be an effective method for the marketing of organically grown prawns as well as those produced by conventional prawn aquaculture without application of the normal chemicals used to prevent black spot. Copyright © 2014 Crown Copyright.

Ecological dynamics of emerging bat virus spillover

Viruses that originate in bats may be the most notorious emerging zoonoses that spill over from wildlife into domestic animals and humans. Understanding how these infections filter through ecological systems to cause disease in humans is of profound importance to public health. Transmission of viruses from bats to humans requires a hierarchy of enabling conditions that connect the distribution of reservoir hosts, viral infection within these hosts, and exposure and susceptibility of recipient hosts. For many emerging bat viruses, spillover also requires viral shedding from bats, and survival of the virus in the environment. Focusing on Hendra virus, but also addressing Nipah virus, Ebola virus, Marburg virus and coronaviruses, we delineate this cross-species spillover dynamic from the within-host processes that drive virus excretion to land-use changes that increase interaction among species. We describe how land-use changes may affect co-occurrence and contact between bats and recipient hosts. Two hypotheses may explain temporal and spatial pulses of virus shedding in bat populations: episodic shedding from persistently infected bats or transient epidemics that occur as virus is transmitted among bat populations. Management of livestock also may affect the probability of exposure and disease. Interventions to decrease the probability of virus spillover can be implemented at multiple levels from targeting the reservoir host to managing recipient host exposure and susceptibility.

Enhanced processing of SPOT multispectral satellite imagery for environmental monitoring and modelling

The Taita Hills in southeastern Kenya form the northernmost part of Africa’s Eastern Arc Mountains, which have been identified by Conservation International as one of the top ten biodiversity hotspots on Earth. As with many areas of the developing world, over recent decades the Taita Hills have experienced significant population growth leading to associated major changes in land use and land cover (LULC), as well as escalating land degradation, particularly soil erosion. Multi-temporal medium resolution multispectral optical satellite data, such as imagery from the SPOT HRV, HRVIR, and HRG sensors, provides a valuable source of information for environmental monitoring and modelling at a landscape level at local and regional scales. However, utilization of multi-temporal SPOT data in quantitative remote sensing studies requires the removal of atmospheric effects and the derivation of surface reflectance factor. Furthermore, for areas of rugged terrain, such as the Taita Hills, topographic correction is necessary to derive comparable reflectance throughout a SPOT scene. Reliable monitoring of LULC change over time and modelling of land degradation and human population distribution and abundance are of crucial importance to sustainable development, natural resource management, biodiversity conservation, and understanding and mitigating climate change and its impacts. The main purpose of this thesis was to develop and validate enhanced processing of SPOT satellite imagery for use in environmental monitoring and modelling at a landscape level, in regions of the developing world with limited ancillary data availability. The Taita Hills formed the application study site, whilst the Helsinki metropolitan region was used as a control site for validation and assessment of the applied atmospheric correction techniques, where multiangular reflectance field measurements were taken and where horizontal visibility meteorological data concurrent with image acquisition were available. The proposed historical empirical line method (HELM) for absolute atmospheric correction was found to be the only applied technique that could derive surface reflectance factor within an RMSE of < 0.02 ps in the SPOT visible and near-infrared bands; an accuracy level identified as a benchmark for successful atmospheric correction. A multi-scale segmentation/object relationship modelling (MSS/ORM) approach was applied to map LULC in the Taita Hills from the multi-temporal SPOT imagery. This object-based procedure was shown to derive significant improvements over a uni-scale maximum-likelihood technique. The derived LULC data was used in combination with low cost GIS geospatial layers describing elevation, rainfall and soil type, to model degradation in the Taita Hills in the form of potential soil loss, utilizing the simple universal soil loss equation (USLE). Furthermore, human population distribution and abundance were modelled with satisfactory results using only SPOT and GIS derived data and non-Gaussian predictive modelling techniques. The SPOT derived LULC data was found to be unnecessary as a predictor because the first and second order image texture measurements had greater power to explain variation in dwelling unit occurrence and abundance. The ability of the procedures to be implemented locally in the developing world using low-cost or freely available data and software was considered. The techniques discussed in this thesis are considered equally applicable to other medium- and high-resolution optical satellite imagery, as well the utilized SPOT data.

Complete genome sequence and integrated protein localization and interaction map for alfalfa dwarf virus, which combines properties of both cytoplasmic and nuclear plant rhabdoviruses

Summary We have determined the full-length 14,491-nucleotide genome sequence of a new plant rhabdovirus, alfalfa dwarf virus (ADV). Seven open reading frames (ORFs) were identified in the antigenomic orientation of the negative-sense, single-stranded viral RNA, in the order 3′-N-P-P3-M-G-P6-L-5′. The ORFs are separated by conserved intergenic regions and the genome coding region is flanked by complementary 3′ leader and 5′ trailer sequences. Phylogenetic analysis of the nucleoprotein amino acid sequence indicated that this alfalfa-infecting rhabdovirus is related to viruses in the genus Cytorhabdovirus. When transiently expressed as GFP fusions in Nicotiana benthamiana leaves, most ADV proteins accumulated in the cell periphery, but unexpectedly P protein was localized exclusively in the nucleus. ADV P protein was shown to have a homotypic, and heterotypic nuclear interactions with N, P3 and M proteins by bimolecular fluorescence complementation. ADV appears unique in that it combines properties of both cytoplasmic and nuclear plant rhabdoviruses.

Commercial-scale Alternaria brown spot resistance screening as the first step in breeding new mandarins for Australia

Rapid screening tests and an appreciation of the simple genetic control of Alternaria brown spot (ABS) susceptibility have existed for many years, and yet the application of this knowledge to commercial-scale breeding programs has been limited. Detached leaf assays were first demonstrated more than 40 years ago and reliable data suggesting a single gene determining susceptibility has been emerging for at least 20 years. However it is only recently that the requirement for genetic resistance in new hybrids has become a priority, following increased disease prevalence in Australian mandarin production areas previously considered too dry for the pathogen. Almost all of the high-fruit-quality parents developed so far by the Queensland-based breeding program are susceptible to ABS necessitating the screening of their progeny to avoid commercialisation of susceptible hybrids. This is done effectively and efficiently by spraying 3-6 month old hybrid seedlings with a spore suspension derived from a toxin-producing field isolate of Alternaria alternate, then incubating these seedlings in a cool room at 25°C and high humidity for 5 days. Susceptible seedlings show clear disease symptoms and are discarded. Analysis of observed and expected segregation ratios loosely support the hypothesis for a single dominant gene for susceptibility, but do not rule out the possibility of alternative genetic models. After implementing the routine screening for ABS resistance for three seasons we now have more than 20,000 hybrids growing in field progeny blocks that have been screened for resistance to the ABS disease.

Efficacy of acibenzolar-S-methyl (Bion®) treatment of Australian commercial passionfruit, Passiflora edulis f. sp. flavicarpa, on resistance to Passionfruit woodiness virus (PWV) and activities of chitinase & β-1,3-glucanase

This greenhouse study investigated the efficacy of acibenzolar-S-methyl (Bion®) treatment of lower leaves of passionfruit, (Passiflora edulis f. sp. flavicarpa), on Passionfruit woodiness disease and activities of two pathogenesis-related proteins, chitinase and β-1,3-glucanase after inoculation with passionfruit woodiness virus (PWV). All Bion® concentrations reduced disease symptoms, but the concentration of 0.025 g active ingredient (a.i.)/l was the most effective, reducing disease severity in systemic leaves by 23, 29 and 30 compared with water-treated controls at 30, 40 and 50 days post inoculation (dpi) with PWV, respectively. Correspondingly, relative virus concentration as determined by DAS-ELISA in the upper, untreated leaves (new growth) above the site of inoculation at 50 dpi was reduced by 17 and 22 in plants treated with 0.025 and 0.05 g a.i./l, respectively. Bion® treatment and subsequent inoculation with PWV increased chitinase and β-1,3-glucanase activities in the new leaves above the site of inoculation at 30 dpi with PWV. It was concluded that optimal protective Bion® treatment concentrations were 0.025 and 0.05 g a.i./l.

Survival of Hendra Virus in the Environment: Modelling the Effect of Temperature

Hendra virus (HeV), a highly pathogenic zoonotic paramyxovirus recently emerged from bats, is a major concern to the horse industry in Australia. Previous research has shown that higher temperatures led to lower virus survival rates in the laboratory. We develop a model of survival of HeV in the environment as influenced by temperature. We used 20 years of daily temperature at six locations spanning the geographic range of reported HeV incidents to simulate the temporal and spatial impacts of temperature on HeV survival. At any location, simulated virus survival was greater in winter than in summer, and in any month of the year, survival was higher in higher latitudes. At any location, year-to-year variation in virus survival 24 h post-excretion was substantial and was as large as the difference between locations. Survival was higher in microhabitats with lower than ambient temperature, and when environmental exposure was shorter. The within-year pattern of virus survival mirrored the cumulative within-year occurrence of reported HeV cases, although there were no overall differences in survival in HeV case years and non-case years. The model examines the effect of temperature in isolation; actual virus survivability will reflect the effect of additional environmental factors

Improving Profitability to Industry Through the Identification and Management of 'tough' Fish Syndrome in Tropical Saddletail Snapper (Lutjanus Malabaricus)

OBJECTIVES: 1. To determine whether incomplete rigor mortis resolution and 'cold shock' play a role in development of tough fish syndrome (TFS) in tropical Saddletail snapper. 2. To identify links between TFS and specific physiological factors in tropical Saddletail snapper. 3. Communicate findings and recommendations to stakeholders and assist with implementation of any changes to fishing or handling practices required.

Molecular evidence of Ebola Reston virus infection in Philippine bats

Background: In 2008-09, evidence of Reston ebolavirus (RESTV) infection was found in domestic pigs and pig workers in the Philippines. With species of bats having been shown to be the cryptic reservoir of filoviruses elsewhere, the Philippine government, in conjunction with the Food and Agriculture Organization of the United Nations, assembled a multi-disciplinary and multi-institutional team to investigate Philippine bats as the possible reservoir of RESTV. Methods: The team undertook surveillance of bat populations at multiple locations during 2010 using both serology and molecular assays. Results: A total of 464 bats from 21 species were sampled. We found both molecular and serologic evidence of RESTV infection in multiple bat species. RNA was detected with quantitative PCR (qPCR) in oropharyngeal swabs taken from Miniopterus schreibersii, with three samples yielding a product on conventional hemi-nested PCR whose sequences differed from a Philippine pig isolate by a single nucleotide. Uncorroborated qPCR detections may indicate RESTV nucleic acid in several additional bat species (M. australis, C. brachyotis and Ch. plicata). We also detected anti-RESTV antibodies in three bats (Acerodon jubatus) using both Western blot and ELISA. Conclusions: The findings suggest that ebolavirus infection is taxonomically widespread in Philippine bats, but the evident low prevalence and low viral load warrants expanded surveillance to elaborate the findings, and more broadly, to determine the taxonomic and geographic occurrence of ebolaviruses in bats in the region. © 2015 Jayme et al.

Flying-Fox Roost Disturbance and Hendra Virus Spillover Risk

Bats of the genus Pteropus (flying-foxes) are the natural host of Hendra virus (HeV) which periodically causes fatal disease in horses and humans in Australia. The increased urban presence of flying-foxes often provokes negative community sentiments because of reduced social amenity and concerns of HeV exposure risk, and has resulted in calls for the dispersal of urban flying-fox roosts. However, it has been hypothesised that disturbance of urban roosts may result in a stress-mediated increase in HeV infection in flying-foxes, and an increased spillover risk. We sought to examine the impact of roost modification and dispersal on HeV infection dynamics and cortisol concentration dynamics in flying-foxes. The data were analysed in generalised linear mixed models using restricted maximum likelihood (REML). The difference in mean HeV prevalence in samples collected before (4.9%), during (4.7%) and after (3.4%) roost disturbance was small and non-significant (P = 0.440). Similarly, the difference in mean urine specific gravity-corrected urinary cortisol concentrations was small and non-significant (before = 22.71 ng/mL, during = 27.17, after = 18.39) (P= 0.550). We did find an underlying association between cortisol concentration and season, and cortisol concentration and region, suggesting that other (plausibly biological or environmental) variables play a role in cortisol concentration dynamics. The effect of roost disturbance on cortisol concentration approached statistical significance for region, suggesting that the relationship is not fixed, and plausibly reflecting the nature and timing of disturbance. We also found a small positive statistical association between HeV excretion status and urinary cortisol concentration. Finally, we found that the level of flying-fox distress associated with roost disturbance reflected the nature and timing of the activity, highlighting the need for a ‘best practice’ approach to dispersal or roost modification activities. The findings usefully inform public discussion and policy development in relation to Hendra virus and flying-fox management.

Routes of Hendra Virus Excretion in Naturally-Infected Flying-Foxes: Implications for Viral Transmission and Spillover Risk

Pteropid bats or flying-foxes (Chiroptera: Pteropodidae) are the natural host of Hendra virus (HeV) which sporadically causes fatal disease in horses and humans in eastern Australia. While there is strong evidence that urine is an important infectious medium that likely drives bat to bat transmission and bat to horse transmission, there is uncertainty about the relative importance of alternative routes of excretion such as nasal and oral secretions, and faeces. Identifying the potential routes of HeV excretion in flying-foxes is important to effectively mitigate equine exposure risk at the bat-horse interface, and in determining transmission rates in host-pathogen models. The aim of this study was to identify the major routes of HeV excretion in naturally infected flying-foxes, and secondarily, to identify between-species variation in excretion prevalence. A total of 2840 flying-foxes from three of the four Australian mainland species (Pteropus alecto, P. poliocephalus and P. scapulatus) were captured and sampled at multiple roost locations in the eastern states of Queensland and New South Wales between 2012 and 2014. A range of biological samples (urine and serum, and urogenital, nasal, oral and rectal swabs) were collected from anaesthetized bats, and tested for HeV RNA using a qRT-PCR assay targeting the M gene. Forty-two P. alecto (n = 1410) had HeV RNA detected in at least one sample, and yielded a total of 78 positive samples, at an overall detection rate of 1.76% across all samples tested in this species (78/4436). The rate of detection, and the amount of viral RNA, was highest in urine samples (>serum, packed haemocytes >faecal >nasal >oral), identifying urine as the most plausible source of infection for flying-foxes and for horses. Detection in a urine sample was more efficient than detection in urogenital swabs, identifying the former as the preferred diagnostic sample. The detection of HeV RNA in serum is consistent with haematogenous spread, and with hypothesised latency and recrudesence in flying-foxes. There were no detections in P. poliocephalus (n = 1168 animals; n = 2958 samples) or P. scapulatus (n = 262 animals; n = 985 samples), suggesting (consistent with other recent studies) that these species are epidemiologically less important than P. alecto in HeV infection dynamics. The study is unprecedented in terms of the individual animal approach, the large sample size, and the use of a molecular assay to directly determine infection status. These features provide a high level of confidence in the veracity of our findings, and a sound basis from which to more precisely target equine risk mitigation strategies.