Automating Crop Damage Assessments with Unmanned Aerial Systems and Machine Learning

Agricultural pests are responsible for millions of dollars in crop losses and management costs every year. In order to implement optimal site-specific treatments and reduce control costs, new methods to accurately monitor and assess pest damage need to be investigated. In this paper we explore the combination of unmanned aerial vehicles (UAV), remote sensing and machine learning techniques as a promising methodology to address this challenge. The deployment of UAVs as a sensor platform is a rapidly growing field of study for biosecurity and precision agriculture applications. In this experiment, a data collection campaign is performed over a sorghum crop severely damaged by white grubs (Coleoptera: Scarabaeidae). The larvae of these scarab beetles feed on the roots of plants, which in turn impairs root exploration of the soil profile. In the field, crop health status could be classified according to three levels: bare soil where plants were decimated, transition zones of reduced plant density and healthy canopy areas. In this study, we describe the UAV platform deployed to collect high-resolution RGB imagery as well as the image processing pipeline implemented to create an orthoimage. An unsupervised machine learning approach is formulated in order to create a meaningful partition of the image into each of the crop levels. The aim of this approach is to simplify the image analysis step by minimizing user input requirements and avoiding the manual data labelling necessary in supervised learning approaches. The implemented algorithm is based on the K-means clustering algorithm. In order to control high-frequency components present in the feature space, a neighbourhood-oriented parameter is introduced by applying Gaussian convolution kernels prior to K-means clustering. The results show the algorithm delivers consistent decision boundaries that classify the field into three clusters, one for each crop health level as shown in Figure 1. The methodology presented in this paper represents a venue for further esearch towards automated crop damage assessments and biosecurity surveillance.

Time-dependent flows of rotating and stratified fluids in geometries with non-uniform cross-sections

Unsteady rotating and stratified flows in geometries with non-uniform cross-sections are investigated under Oseen approximation using Laplace transform technique. The solutions are obtained in closed form and they reveal that the flow remains oscillatory even after infinitely large time. The existence of inertial waves propagating in both positive and negative directions of the flow is observed. When the Rossby or Froude number is close to a certain infinite set of critical values the blocking and back flow occur and the flow pattern becomes more and more complicated with increasing number of stagnant zones when each critical value is crossed. The analogy that is observed in the solutions for rotating and stratified flows is also discussed.

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.

The community of medium-sized carnivores: the interactions between species, habitats and rabies

The composition of the carnivore community influences the different forms of inter-specific interactions. Furthermore, inter-specific interactions of carnivores have important implications for intra-guild competition, epidemiology and strategies of species-specific population management. Zoonooses, such as rabies, are diseases that can be transmitted from wildlife to people. Knowing the ecological characteristics of the species helps us to choose the right preventive actions and to time them accurately. In this thesis, I have studied how raccoon dogs Nyctereutes procyonoides, European badgers Meles meles, red foxes Vulpes vulpes and domestic cats Felis silvestris catus act as members of carnivore community, and how these interactions relate to the transmission risk of rabies. In the study area, these species form a community of medium-sized and rather generalist predators. They live in the same areas, in spatially and temporally overlapping home ranges and use the same habitats and dens and even have similar diets. However, there is no direct evidence of competition. Shared dens point to good tolerance of other species. Numerous observations of animals moving in each other’s proximity give similar clues. However, overlapping home ranges and similar habitat preferences lead to frequent inter-specific contacts, which increase the risk of possible rabies transmission. Also, the new insight of habitat use gained by this study illustrates the similar favouring of deciduous forests and fields by these sympatric medium-sized carnivores, creating a basis for contact zones, i.e. risky habitats for rabies transmission and spread. This study is so far the only simultaneous radio tracking study of raccoon dogs, badgers, foxes and cats. These results give new insight of the interactions in the carnivore community, as well as of the behaviour of each individual species. Also, these results have significant implications for the planning of rabies control. In order to reach viable management decisions, not only one or two species should be taken into consideration, but the whole community. In particular, this changes the perspective to inter-specific contacts, animal densities, densities of individuals susceptible to diseases and the magnitude of preventive actions. Rabies should be considered as a multi-vector disease, at least in Finland and the Baltic states. It is of interest for disease management to be able to model an epizootic with local parameters to reflect the real situation and also to suite best the local management needs.

Virtual Special Issue: “The impact of temperature on human health: new insights”

Ambient temperature is one of the basic parameters characterising human comfort: are we too hot, too cold, or just right? The impact of temperature goes beyond comfort: inadequate temperature and temperature variations have consequences on human health, as the increasing numbers of studies have demonstrated. The topic is of particular significance at the times when climate change shifts the traditional – as we know them- temperature zones, and brings much wider temperature variations. For these reasons the impact of temperature on health has been one of the most popular topics among the articles submitted and published in Science of the Total Environment over the last few years. This Virtual Special Issue compiles 18 articles published in our journal on this topic since 2012. It is worth briefly summarizing the rich scientific insights brought by these articles, as well as broader considerations, particularly those extending to management, discussed by the authors of the articles.

The range expansion of the European map butterfly in Finland

Climate change will influence the living conditions of all life on Earth. For some species the change in the environmental conditions that has occurred so far has already increased the risk of extinction, and the extinction risk is predicted to increase for large numbers of species in the future. Some species may have time to adapt to the changing environmental conditions, but the rate and magnitude of the change are too great to allow many species to survive via evolutionary changes. Species responses to climate change have been documented for some decades. Some groups of species, like many insects, respond readily to changes in temperature conditions and have shifted their distributions northwards to new climatically suitable regions. Such range shifts have been well documented especially in temperate zones. In this context, butterflies have been studied more than any other group of species, partly for the reason that their past geographical ranges are well documented, which facilitates species-climate modelling and other analyses. The aim of the modelling studies is to examine to what extent shifts in species distributions can be explained by climatic and other factors. Models can also be used to predict the future distributions of species. In this thesis, I have studied the response to climate change of one species of butterfly within one geographically restricted area. The study species, the European map butterfly (Araschnia levana), has expanded rapidly northwards in Finland during the last two decades. I used statistical and dynamic modelling approaches in combination with field studies to analyse the effects of climate warming and landscape structure on the expansion. I studied possible role of molecular variation in phosphoglucose isomerase (PGI), a glycolytic enzyme affecting flight metabolism and thereby flight performance, in the observed expansion of the map butterfly at two separate expansion fronts in Finland. The expansion rate of the map butterfly was shown to be correlated with the frequency of warmer than average summers during the study period. The result is in line with the greater probability of occurrence of the second generation during warm summers and previous results on this species showing greater mobility of the second than first generation individuals. The results of a field study in this thesis indicated low mobility of the first generation butterflies. Climatic variables alone were not sufficient to explain the observed expansion in Finland. There are also problems in transferring the climate model to new regions from the ones from which data were available to construct the model. The climate model predicted a wider distribution in the south-western part of Finland than what has been observed. Dynamic modelling of the expansion in response to landscape structure suggested that habitat and landscape structure influence the rate of expansion. In southern Finland the landscape structure may have slowed down the expansion rate. The results on PGI suggested that allelic variation in this enzyme may influence flight performance and thereby the rate of expansion. Genetic differences of the populations at the two expansion fronts may explain at least partly the observed differences in the rate of expansion. Individuals with the genotype associated with high flight metabolic rate were most frequent in eastern Finland, where the rate of range expansion has been highest.

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.

Generation and Characterization of the Cation-Chloride Cotransporter KCC2 Hypomorphic Mouse

The cation-Cl- cotransporter (CCC) family comprises of Na+-Cl- cotransporter (NCC), Na+-K+-2Cl- cotransporters (NKCC1-2), and four K+-Cl- cotransporters (KCC1-4). These proteins are involved in several physiological activities, such as cell volume regulation. In neuronal tissues, NKCC1 and KCC2 are important in determining the intracellular Cl- levels and hence the neuronal responses to inhibitory neurotransmitters GABA and glycine. One aim of the work was to elucidate the roles for CCC isoforms in the control of nervous system development. KCC2 mRNA was shown to be developmentally up-regulated and follow neuronal maturation, whereas NKCC1 and KCC4 transcripts were highly expressed in the proliferative zones of subcortical regions. KCC1 and KCC3 mRNA displayed low expression throughout the embryogenesis. These expression profiles suggest a role for CCC isoforms in maturation of synaptic responses and in the regulation of neuronal proliferation during embryogenesis. The major aim of this work was to study the biological consequences of KCC2-deficiency in the adult CNS, by generating transgenic mice retaining 15-20% of normal KCC2 levels. In addition, by using these mice as a tool for in vivo pharmacological analysis, we investigated the requirements for KCC2 in tonic versus phasic GABAA receptor-mediated inhibition. KCC2-deficient mice displayed normal reproduction and life span, but showed several behavioral abnormalities, including increased anxiety-like behavior, impaired performance in water maze, alterations in nociceptive processing, and increased seizure susceptibility. In contrast, the mice displayed apparently normal spontaneous locomotor activity and motor coordination. Pharmacological analysis of KCC2-deficient mice revealed reduced sensititivity to diazepam, but normal gaboxadol-induced sedation, neurosteroid hypnosis and alcohol-induced motor impairment. Electrophysiological recordings from CA1-CA3 subregions of the hippocampus showed that KCC2 deficiency affected the reversal potentials of both the phasic and tonic GABA currents, and that the tonic conductance was not affected. The results suggest that requirement for KCC2 in GABAergic neurotransmission may differ among several functional systems in the CNS, which is possibly due to the more critical role of KCC2 activity in phasic compared to tonic GABAergic inhibition.

An annotated checklist of amphibian of India with some insights into the patterns of species discoveries, distribution and endemism

An annotated checklist of 284 species of amphibians of India accommodated under 50 genera and 14 families is provided. Synonyms, English names, type localities, deposition of type specimens, type specimen availability and distributional records in India and outside India are provided for all the species. Among the 284 species of amphibians from India, 132 are endemic to Western Ghats; 29 to Northeastern India; and 5 to Andaman Nicobar islands. Species discovery patterns from the various biogeographic zones in India are discussed in detail. Cumulative discovery pattern with special reference to the genera Fejervarya (17 species), Nyctibatrachus (16 species), Indirana (10 species), Micrixalus (11 species), Philautus (46 species) and Gegeneophis (10 species) are also discussed.

Mixed mode (I and II) crack tip fields in bulk metallic glasses

Stationary crack tip fields in bulk metallic glasses under mixed mode (I and II) loading are studied through detailed finite element simulations assuming plane strain, small scale yielding conditions. The influence of internal friction or pressure sensitivity on the plastic zones. notch deformation, stress and plastic strain fields is examined for different mode mixities. Under mixed mode loading, the notch deforms into a shape such that one part of its surface sharpens while the other part blunts. Increase in mode If component of loading dramatically enhances the normalized plastic zone size, lowers the stresses but significantly elevates the plastic strain levels near the notch tip. Higher internal friction reduces the peak tangential stress but increases the plastic strain and stretching near the blunted part of the notch. The simulated shear bands are straight and extend over a long distance ahead of the notch tip under mode II dominant loading. The possible variations of fracture toughness with mode mixity corresponding to failure by brittle micro-cracking and ductile shear banding are predicted employing two simple fracture criteria. The salient results from finite element simulations are validated by comparison with those from mixed mode (I and II) fracture experiments on a Zr-based bulk metallic glass.

CFD analysis of high frequency miniature pulse tube refrigerators for space applications with thermal non-equilibrium model

High frequency, miniature, pulse tube cryocoolers are extensively used in space applications because of their simplicity. Parametric studies of inertance type pulse tube cooler are performed with different length-to-diameter ratios of the pulse tube with the help of the FLUENT (R) package. The local thermal non-equilibrium of the gas and the matrix is taken into account for the modeling of porous zones, in addition to the wall thickness of the components. Dynamic characteristics and the actual mechanism of energy transfer in pulse are examined with the help of the pulse tube wall time constant. The heat interaction between pulse tube wall and the oscillating gas, leading to surface heat pumping, is quantified. The axial heat conduction is found to reduce the performance of the pulse tube refrigerator. The thermal non-equilibrium predicts a higher cold heat exchanger temperature compared to thermal equilibrium. The pressure drop through the porous medium has a strong non-linear effect due to the dominating influence of Forchheimer term over that of the linear Darcy term at high operating frequencies. The phase angle relationships among the pressure, temperature and the mass flow rate in the porous zones are also important in determining the performance of pulse tuberefrigerator.

Predicting the temporal response of seagrass meadows to dredging using Dynamic Bayesian Networks

Predicting temporal responses of ecosystems to disturbances associated with industrial activities is critical for their management and conservation. However, prediction of ecosystem responses is challenging due to the complexity and potential non-linearities stemming from interactions between system components and multiple environmental drivers. Prediction is particularly difficult for marine ecosystems due to their often highly variable and complex natures and large uncertainties surrounding their dynamic responses. Consequently, current management of such systems often rely on expert judgement and/or complex quantitative models that consider only a subset of the relevant ecological processes. Hence there exists an urgent need for the development of whole-of-systems predictive models to support decision and policy makers in managing complex marine systems in the context of industry based disturbances. This paper presents Dynamic Bayesian Networks (DBNs) for predicting the temporal response of a marine ecosystem to anthropogenic disturbances. The DBN provides a visual representation of the problem domain in terms of factors (parts of the ecosystem) and their relationships. These relationships are quantified via Conditional Probability Tables (CPTs), which estimate the variability and uncertainty in the distribution of each factor. The combination of qualitative visual and quantitative elements in a DBN facilitates the integration of a wide array of data, published and expert knowledge and other models. Such multiple sources are often essential as one single source of information is rarely sufficient to cover the diverse range of factors relevant to a management task. Here, a DBN model is developed for tropical, annual Halophila and temperate, persistent Amphibolis seagrass meadows to inform dredging management and help meet environmental guidelines. Specifically, the impacts of capital (e.g. new port development) and maintenance (e.g. maintaining channel depths in established ports) dredging is evaluated with respect to the risk of permanent loss, defined as no recovery within 5 years (Environmental Protection Agency guidelines). The model is developed using expert knowledge, existing literature, statistical models of environmental light, and experimental data. The model is then demonstrated in a case study through the analysis of a variety of dredging, environmental and seagrass ecosystem recovery scenarios. In spatial zones significantly affected by dredging, such as the zone of moderate impact, shoot density has a very high probability of being driven to zero by capital dredging due to the duration of such dredging. Here, fast growing Halophila species can recover, however, the probability of recovery depends on the presence of seed banks. On the other hand, slow growing Amphibolis meadows have a high probability of suffering permanent loss. However, in the maintenance dredging scenario, due to the shorter duration of dredging, Amphibolis is better able to resist the impacts of dredging. For both types of seagrass meadows, the probability of loss was strongly dependent on the biological and ecological status of the meadow, as well as environmental conditions post-dredging. The ability to predict the ecosystem response under cumulative, non-linear interactions across a complex ecosystem highlights the utility of DBNs for decision support and environmental management.

The influence of the Coriolis force on flux tubes rising through the solar convection zone

In order to study the effect of the Coriolis force due to solar rotation on rising magnetic flux, the authors consider a flux ring, azimuthally symmetric around the rotation axis, starting from rest at the bottom of the convection zone, and then follow the trajectory of the flux ring as it rises. If it is assumed that the flux ring remains azimuthally symmetric during its ascent, then the problem can be described essentially in terms of two parameters: the value of the initial magnetic field in the ring when it starts, and the effective drag experienced by it. For field strengths at the bottom of the convection zone of order 10,000 G or less, it is found that the Coriolis force plays a dominant role and flux rings starting from low latitudes at the bottom are deflected and emerge at latitudes significantly poleward of sunspot zones.

Compatibility of RuO2 electrodes with PZT ceramics

Because of its high electrical conductivity and good diffusion barrier properties ruthenium dioxide (RuO2) is a good electrode material for use with ferroelectric lead zirconate-titanate (PZT) solid solutions. Under certain conditions, RuO2 can react with PZT to form lead ruthenate (Pb2Ru2O6.5) during processing at elevated temperatures resulting in lead depletion from PZT. The standard Gibbs energies of formation of RuO2 and Pb2Ru2O6.5 and activities of components of the PZT solid solution have been determined recently. Using this data along with older thermodynamic information on PbZrO3 and PbTiO3, the stability domain of Pb2Ru2O6.5 is computed as a function of PZT composition, temperature and oxygen partial pressure in the gas phase. The results show PbZrO3-rich compositions are more prone to react with RuO2 at all temperatures. Increasing temperature and decreasing oxygen partial pressure suppress the reaction. Graphically displayed are the reaction zones as a function of oxygen partial pressure and PZT composition at temperatures 973, 1173 and 1373 K.

Free-volume dependent pressure sensitivity of Zr-based bulk metallic glass

Instrumented indentation experiments on a Zr-based bulk metallic glass (BMG) in as-cast, shot-peened and structurally relaxed conditions were conducted to examine the dependence of plastic deformation on its structural state. Results show significant differences in hardness, H, with structural relaxation increasing it and shot peening markedly reducing it, and slightly changed morphology of shear bands around the indents. This is in contrast to uniaxial compressive yield strength, sigma(y), which remains invariant with the change in the structural state of the alloys investigated. The plastic constraint factor, C = H/sigma(y), of the relaxed BMG increases compared with that of the as-cast glass, indicating enhanced pressure sensitivity upon annealing. In contrast, C of the shot-peened layer was found to be similar to that observed in crystalline metals, indicating that severe plastic deformation could eliminate pressure sensitivity. Microscopic origins for this result, in terms of shear transformation zones and free volume, are discussed.