Next Article in Journal
Analysis and Modeling of the Complex Dielectric Constant of Bound Water with Application in Soil Microwave Remote Sensing
Previous Article in Journal
A Constrained Convex Optimization Approach to Hyperspectral Image Restoration with Hybrid Spatio-Spectral Regularization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geomatics and EO Data to Support Wildlife Diseases Assessment at Landscape Level: A Pilot Experience to Map Infectious Keratoconjunctivitis in Chamois and Phenological Trends in Aosta Valley (NW Italy)

1
Department of Agricultural, Forest and Food Sciences (DISAFA), GEO4Agri DISAFA Lab, Unito Green Office Energy & Climate Change Groups, Università degli Studi di Torino, Largo Paolo Braccini 2, 10095 Grugliasco (TO), Italy
2
Istituto Zooprofilattico Sperimentale Piemonte, Liguria e Valle d’Aosta (IZS PLV) SC Aosta—CeRMAS (National Reference Center for Wildlife Diseases), Località Amerique, 7/C, 11020 Quart (AO), Italy
3
Department of Veterinary Sciences (DSV), Università degli studi di Torino, Largo Paolo Braccini 2, 10095 Grugliasco (TO), Italy
4
Department of Agricultural, Forest and Food Sciences (DISAFA), GEO4Agri DISAFA Lab, Università degli studi di Torino, Largo Paolo Braccini 2, 10095 Grugliasco (TO), Italy
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(21), 3542; https://0-doi-org.brum.beds.ac.uk/10.3390/rs12213542
Submission received: 24 September 2020 / Revised: 21 October 2020 / Accepted: 26 October 2020 / Published: 29 October 2020

Abstract

:
Geomatics and satellite remote sensing offer useful analysis tools for several technical-scientific fields. This work, with reference to a regional case of study, investigates remote sensing potentialities for describing relationships between environment and diseases affecting wildlife at landscape level in the light of climate change effects onto vegetation. Specifically, the infectious keratoconjunctivitis (IKC) of chamois (Rupicapra rupicapra L.) in Aosta Valley (NW Italy) was investigated at the regional level. IKC (Mycoplasma conjunctivae) is a contagious disease for domestic and wild ruminants (Caprinae and Ovinae). Two types of analysis were performed: one aimed at exploring by remotely sensed data phenological metrics (PMs) and evapotranspiration (ET) trends of vegetation in the area; one investigating the correlation between PMs and ET, versus IKC prevalence. The analysis was based on TERRA MODIS image time series ranging from 2000 to 2019. Ground data about IKC were available for a shorter time range: 2009–2019. Consequently, PMs and ET trend investigations were focused on the whole times range (2000–2019); conversely, correlation analysis was achieved with reference to the reduced 2009–2019 period. The whole study was based on freely available data from public archives. MODIS products, namely MOD13Q1 v.6 and MOD16A2, were used to derive PM and ET trends, respectively. Shuttle Radar Topography Mission (SRTM) Digital Terrain Model (DTM) was used to describe local topography; CORINE Land Cover map was adopted to describe land use classes. PMs and ET (as derivable from EO data) proved to significantly changed their values in the last 20 years, with a continuous progressive trend. As far as correlation analysis was concerned, ET and some PMs (specifically, End of Season (EOS) and Length of Season (LOS) proved significantly condition IKC prevalence. According to results, the proposed methodology can be retained as an effective tool for supporting public health and eco-pathological sectors. Specifically, it can be intended for a continuous monitoring of effects that climatic dynamics determine onto wild animals in the Alpine area, included diseases and zoonosis, moving future environmental management and planning towards the One Health perspective.

Graphical Abstract

1. Introduction

Geomatics and satellite remote sensing represent a useful analysis tool in several technical-scientific fields [1]. Nowadays, remote sensing is widely used in many fields like agronomy, forestry, and the environment, in general. Nevertheless, veterinary and faunistic-related applications are still limited and, often, characterized by an improper use of Earth Observation (EO) data, that makes a significant improvement desirable [2].
On the international scene, only few research groups are currently properly exploiting all the potentialities that geomatics (remote sensing included) and digital geographical data could offer to the veterinary sector. Most of works concern parasitology and virology for etiological/epidemiological studies [3], diagnosis, and medical history [4]. In these contexts, satellite remote sensing has assumed a great interest in the last years [5]; EO data are proficiently used to feed meteorological and climatological models with the aim of generating predictive scenarios of zoonosis spread and outbreaks [6].
Initially, most of studies in the veterinary and health sectors used remote sensing to describe environmental conditions; this occurred especially with reference to malaria in Africa and Asia [7,8]. Presently, epidemiologists are adopting remote sensing to investigate a variety of vector-borne diseases. Associations between remote sensing-derived environmental variables (e.g., temperature, humidity, land cover, etc.) and vector density are used to map and characterize vector habitats [9]. The basic idea is that remotely sensed data can contain dynamic predictors of Earth’s processes suitable for describing niche preferences of some medically important host diseases mechanisms. Moreover, because of their continuity of acquisition, remotely sensed data provide a synoptic representation of environment at proper spatial and temporal scales [10].
Meteorological and EO data are often jointly used for diseases analysis. For example, outbreaks of diarrheal disease and, specifically, cholera were analyzed by a new modeling approach based on satellite data to produce cholera risk maps in several regions of globe [8], supporting the idea that the ongoing EO technology transfer is making possible to investigate new patterns in a systemic point of view [11].
Given the veterinary and public health impact of vector-borne diseases, there is a clear and immediate need to map and monitor local landscape attitude to encourage emergence and spread of these diseases. Current approaches for predicting disease risks are mostly neglecting key features of landscape related to the functional habitat of vectors (or hosts) and, hence, of the pathogen [12].
A global satellite-based monitoring of proper climate variables could help to map occurring anomalies with the aim of predicting spatial distribution of risk related to emergence and propagation of disease vectors. Such information could provide sufficient lead-time for outbreak prevention and potentially reduce burden and spread of ecologically coupled diseases.
Additionally, remote sensing could have an important role in the comprehension of patho-system dynamic. With reference to the so called “disease triangle” (Figure 1), including host, pathogens/vectors and environment, remote sensing, and geomatics in general, could support scientists and decision-makers to better understand the role of environmental patterns and, therefore, explore its complex relations with the other parts of patho-system [13,14].
A good example is represented by atmospheric pollution that was recognized to increase sensitivity to pulmonary diseases, as the last pandemic event, coronavirus (SARS-CoV-2), has suggested [15].
GIS studies about endo- and ecto-parasitoses of veterinary interest, with particular reference to zoonoses agents, represent today the greatest contribution to veterinary and faunistic sectors [16]. For many years, the World Organization for Animal Health located in Paris (OIE) and the World Health Organization (WHO) in Geneva have been underlining the importance of geomatics and remote sensing applications [17] in the One Health perspective.
This work, with reference to a regional case study, investigates remote sensing potentialities for describing relationships between environment and diseases affecting wildlife at landscape level. Moreover, it is intended to describe the effects of climate change onto the vegetation component, with special concern about pastures. The study area corresponds to the entire Aosta Valley Region located in the Italian Western Alps. In particular, a new analysis approach is presented to operate at landscape level to analyze if and how environmental factors could condition the occurrence of infectious keratoconjunctivitis (IKC, Mycoplasma conjunctivae) in chamois. IKC is a contagious disease for domestic and wild ruminants (Caprinae and Ovinae) [18]. In chamois, the disease can be serious [19] and, as in other wild ruminants, blindness can occur [20], with consequent death of the animal from trauma (e.g., fall from cliffs or starvation) [21]. The period of mountain pasture is risky for the potential contact between domestic and wild infected animals; over the years, several outbreaks have been reported in wild ungulates in the Alps [22] and this is the reason that makes monitoring/surveillance plans still active. IKC caused by Mycoplasma conjunctivae is a complex disease of domestic and wild Caprinae, with great variations in the clinic-pathological and epidemiological picture. In wildlife, IKC is sometimes associated with high mortality [23,24]. It has been suggested that the pathogenesis of IKC is influenced by host predispositions, virulence of M. conjunctivae strains, secondary infections, and environmental factors [25]. Sex and age imbalance in affected populations were observed in severe outbreaks [26], indicating that age and social behavior, including sexual segregation, may be important risk factors. Differently, differences in virulence between different strains do not seem to play a major role; mycoplasmal load is obviously associated to the presence and severity of signs. However, the driver of mycoplasmal multiplication in the host is unknown. Environmental factors might have a role, regarding both the expression of the disease in individual cases and the onset of an outbreak in a population [24]. The underlying hypothesis of this work is that remote sensing could support comprehension of the role of environmental patterns in conditioning IKC patho-system, and related pathologies, as for other diseases. Altitude, air quality, and UV light have been discussed as possible predisposing factors for IKC in wild ungulates along with overcrowding [27]. Multiple outbreaks of IKC in Alpine ibex and Alpine chamois populations have been described in literature [28]. Different outbreaks of infectious keratoconjunctivitis (IKC) affecting alpine chamois and ibex in the western and central Swiss Alps and Aosta Valley were recorded in the period 2001–2019 [29]. Between the years 2001 and 2003, in Switzerland, Mycoplasma conjunctivae was identified from conjunctival swabs by means of a nested PCR in 27 of the 28 chamois tested. The outbreaks occurred in an area covering 1590 km2. Deep valleys acted as a barrier to the spread of the disease. Many of the affected animals were juveniles, and more females than males died of IKC. The disease was more common during the summer and autumn. In some outbreaks, mortality can reach 30 percent, as, for example, in chamois in Italy, France, and Switzerland, and hundreds of chamois may die. Major outbreaks were recorded in 2001–2003 [30] and 2016–2018 [20]. With these premises, in this work, two types of analysis were performed: one aimed at exploring, by remotely sensed data, phenological metrics (PMs) and evapotranspiration (ET) trends of vegetation; one investigating correlation between PMs and ET versus IKC prevalence. PMs/ET analysis was based on TERRA MODIS image time series ranging from 2000 to 2019. Ground data about IKC were available for a shorter time range: 2009–2019. Consequently, PMs and ET trends investigation were done for the whole times range (2000–2019); conversely, correlation analysis was achieved with reference to the 2009–2019 period.

2. Materials and Methods

2.1. Study Area

The study area corresponds to the entire Aosta Valley, an administrative region located in the Northern West Alps of Italy, close to the border with France and Switzerland. It is an alpine region that hosts the highest peaks in Europe (Figure 2). It sizes about 3263 km2 and has a population of about 126,000 people [31]; it is the smallest, least populous, and least densely populated region of Italy. Despite of this, it is one of the most abundant areas in terms of fauna and flora, and biodiversity in general, hosting many protected areas such like Gran Paradiso National Park, Mont Avic Regional Park, and Mont Mars reserve. For this reason, it is an open-cell laboratory especially for environmental and biological sciences. Consequently, wildlife diseases and zoonosis can be proficiently studied in this area. In Figure 3, boundaries of regional protected areas are reported.
Altitude ranges between 340 m (bottom of the valley, Pont Saint Martin municipality) up to 4810 m a.s.l. (Mont Blanc). Average altitude is about 2100 m; the 60% of the territory is above 2000 m a.s.l. Aosta Valley is completely bordered by mountains, often exceeding 4000 m a.s.l.: the Mont Blanc massif in the west, the Monte Rosa and Matterhorn (Cervino) massifs in the north-east, the Gran Paradiso massif in the south, and the Grand Combin (in Switzerland) in the north. Table 1 reports altitude ranges and their spatial distribution, respectively.
Aosta Valley climate is strongly influenced by topography. High mountains that surround the region prevent the access of humid air masses of Mediterranean or Atlantic origin. This makes the local climate characterized by a high degree of aridity in the central area (center of the main bottom valley), with rainfalls lower than 500 mm·y−1; border areas, in particular south-eastern ones and north-western valleys, differently, show an average rainfall above 1400 mm·y−1 [32].
In winter, precipitations are mainly snowy. In summer, convective rainfalls are quite common, determining frequent thunderstorms; spring and autumn are characterized by stratified rainfalls, possibly lasting for several days, with a consequent increase of flood risk.
In force of its highly variegated climatic situation, Aosta Valley land cover is heterogenous. In the present study, particular attention was paid to low natural vegetation where IKC diffusion is highly possible due to the interaction and competition for grasses among animals, in particular wildlife ruminants and breed animals [33,34,35,36]. In Aosta Valley, where breeding activities and ungulates wildlife are very abundant, animal competitive interaction generally occurs in mesotrophic and eutrophic pastures or areas with good pastoral values [37,38,39]. Under this hypothesis, supported only by empirical evidence, an attempt was made to analyze whether a “favorable” type of land cover could play as hotspot for diseases spreading. From an ecological point of view, a greater competitive interaction among animals for the same resource was expected in these areas. Under this scenario, environmental analysis could drive to a better comprehension of the epidemiological relationships linking disease and hosts.

2.2. Veterinary Ground Samples

The most of chamois samples were collected at the hunting wildlife control centers of RFD (Regional Forestry Districts) managed by Corpo Forestale della Valle d’Aosta (Forest Guards). Samples were obtained from chamois that were shot by hunters and analyzed by veterinary officers of CeRMAS (National Reference Center for Wildlife Diseases). Map coordinates of locations where animals were collected were not available. Consequently, data from all the analyses were achieved at the regional level.
As far as chamois analysis is concerned, swabs were collected under the third eyelid from both eyes, transported and stored refrigerated until analyzed. Ocular clinical signs of conjunctivitis and keratitis (ocular damage, inflammation or discharged) or signs of blindness as indicated by abnormal behavior such as stumbling, circling, uncertain gait or inability to climb were recorded at sampling. Age, sex, and body condition were also recorded. Some samples came from chamois that were found dead at the ground showing signs of conjunctivitis or keratitis. The sampling period ranges between 2009 and 2019, when several severe outbreaks of chamois IKC occurred in the Alps [24,40].
At the laboratory, eye swabs were placed into sterile tubes with 0.45 mL of lysis buffer (1 M Tris–HCl, pH 8.5, 0.5% Tween 20, 0.24 mg/mL proteinase K) and mixed for 30 s. Cells were lysed for 60 min at 60 °C and, successively, heated up to 95 °C for 15 min in order to inactivate proteinase K. Obtained lysates were tested with a specific real time PCR to assess the presence of M. conjunctivae; according to Vilei et al. 2007 [41] TaqMan real time PCR reactions were performed by using 2.5 μL of test sample, 900 nM of lppS forward primer (5′-CAGCTGGTGTAGCACTTTTTGC-3′) and lppS reverse primer (5′-TTAACACCTATGCTCTCGTCTTTGA-3′), 300 nM of lppS probe (5′-TGCTTCGACTACCAAATATGATGGTGATCCTCT-3′ with 6FAM reporter dye and TAMRA quencher affixed on the 5′ and 3′ ends, respectively), and TaqMan Universal PCR Master Mix in a 25 μL volume. An exogenous Internal Positive Control was introduced for all reactions to check for the presence of eventual PCR inhibitors. PCR reactions were run by StepOne Plus instrument (Thermo Fisher) using the following cycling parameters: one step at 50 °C for 2 min and at 95 °C for 10 min, 40 cycles of denaturation at 95 °C for 15 s, and extension at 60 °C for 1 min were performed. Real-time fluorescence measurements were taken for each sample by using the StepOne™ Software v2.3 (Thermo Fisher Scientific: https://www.thermofisher.com/order/catalog/product/4376357#/4376357) and the PCR cycle number at which the fluorescent signal crossed the cycle threshold (set manually) was recorded as CT value. The fluorescence emission baseline was set manually two cycles before the cycle with a significant fluorescence signal. The specificity of the TaqMan assay was evaluated by testing genomic DNA of mycoplasmas other than M. conjunctivae and of other ocular pathogens [41,42].
Prevalence data were finally computed by Equation (1) with reference to the entire regional territory for all the monitored years (2009–2019).
Pr =   C P   × 100
where Pr = disease prevalence (%), C = number of positive disease cases detected by PCR and optical analysis from samples, P = number of examined chamois.

2.3. EO and Geographical Digital Data

MOD13Q1 v.6 [43] product from the NASA TERRA Moderate Resolution Imaging Spectroradiometer (MODIS) mission was used to map the Normalized Difference Vegetation Index (NDVI) over the area in the period 18 February 2000–31 December 2019. Four hundred-fifty two MOD13Q1 images were obtained for free by Google Earth Engine [44]. They were stacked into a NDVI time series (NTS) and filtered according to the Pixel Reliability layer, supplied together with the product, to map reliability of each scene pixel. Stacking and filtering (Savitzky-Golay, [45,46]) were operated by a self-developed IDL 8.1 routine [47]. Vegetation Indices from MOD13Q1 v.6 product are a composite one, having a time step of 16 days and a spatial resolution of 250 m. The MOD13Q1 v.6 product provides a Vegetation Index (VI) value at a per pixel basis. It is referred to as a continuity index to the existing one from the National Oceanic and Atmospheric Administration-Advanced Very High Resolution Radiometer (NOAA-AVHRR) data. Composition algorithm selects the best available pixel value from all the acquisitions within the considered 16 days period. Selection criteria are the following: no clouds, low viewing angle, highest NDVI value adopting the classic formula as follow: Equation (2).
NDVI   =   ( ρ NIR ρ RED )   ( ρ NIR   +   ρ RED )
where ρNIR and ρRED are the at-the-ground reflectance of the near-infrared (MODIS band 2) and red (MODIS band 1) bands, respectively.
Evapotranspiration (ET) maps were retrieved from the MOD16A2 v.6 collection [48] by Google Earth Engine for the period 1 January 2000–31 December 2019 for a total of 868 images.
MOD16A2 v.6 Evapotranspiration/Latent Heat Flux product is an 8-day composite product with a geometric resolution of 500 m. The algorithm adopted for the MOD16A2 data product collection refers to the Penman–Monteith equation [49], which includes inputs of daily meteorological reanalysis data such as albedo, land cover, vegetation property dynamics, and Land Surface Temperature (LST) that in some studies were adopted to evaluate ET [50] or heat fluxes [51]. Pixel values for the Net Evapotranspiration (ET0) is the sum of all 8 days within the composite period expressed in kg m−2 8 d−1. MOD16A2 v.6 ET layers were stacked into an ET time series (ETS); no filtering was applied. Stacking was operated by a self-developed IDL 8.1 routine.
The Shuttle Radar Topography Mission (SRTM) digital elevation model was used for this work [52,53]. SRTM v3 product (SRTM Plus) is provided by NASA JPL with a grid size of 1 arc-second (approximately 30 m). Native SRTM was resampled at 30 m (SRTM30) to make pixel squared.
All data were converted from native geographical reference systems into the ED50 UTM 32N one.
MOD13Q1 v.6 and MOD16A2 v.6 layers were geometrically oversampled up to 30 m to refine area zonation.

2.4. Land Cover Data

To properly describe land cover at the regional level, the Corine Land Cover 2018 dataset (hereinafter called CLC2018) was used; CORINE (Coordination of Information on the Environment) Land Cover inventory was initiated in 1985 to standardize data collection on land in Europe to support environmental policy development. The project is coordinated by the European Environment Agency (EEA) in the frame of the EU Copernicus program and implemented by national teams. The number of participating countries has increased over time currently including 33 member countries and six cooperating countries with a total area of over 5.8 million km2. CLC2018 specifically, is one of the available datasets produced within the general CORINE frame and refers about land cover/land use status in 2018. The reference year of the first CLC inventory was 1990 and the first update was achieved in 2000. The current update cycle is 6 years. Satellite imagery provides the geometrical and thematic basis of maps; it is integrated with in-situ data as essential ancillary information. CLC Level 3 dataset consider 44 classes with a minimum mapping unit (MMU) of 25 hectares (100 m minimum mapping width). CLC2018 Level 3 of Aosta Valley was obtained from the collection COPERNICUS/CORINE/V20/100m by Google Earth Engine in grid format (100 m grid size) [54,55,56]. CLC2018 Level 3 classes and correspondent statistics are reported, for the whole Aosta Valley, in Table 2.
LC classes were considered with reference to altitude ranges, in order to give a more comprehensive description on the ongoing processes affecting vegetation and its relationship with IKC prevalence.
As indicated in Table 2, this study only considered the following CORINE land cover types: pastures, natural grassland, moors, and heathlands. These classes were selected since majorly favoring the interaction between breed animals and ungulates. In Aosta Valley, where breeding activities and ungulates wildlife are very abundant, this competitive interaction generally occurs in mesotrophic and eutrophic pastures or areas with good pastoral values. Under this hypothesis, supported only by empirical evidence, an attempt was made to analyze whether a favorable type of land cover could play as a hotspot for diseases spreading. From an ecological point of view, a greater competitive interaction among animals for the same resource was expected in these areas. Under this scenario, environmental analysis could drive to a better comprehension of the epidemiological relationships linking disease and hosts.
CLC2018 and SRTM30 were used as descriptors of environmental conditions possibly conditioning occurrences of IKC. Consequently, IKC was preliminarily tested against favorable land cover type extent with respect to 3 different altitude classes (Table 3).

2.5. Methodology

To investigate possible relationships of IKC prevalence with vegetation-related factors, two analyses were performed. One (hereinafter called Analysis 1) was aimed at testing if any significant climatic trend could be recognized affecting phenological metrics (PMs, see forward on) and ET (as measured by NTS and ETS, respectively), in the period 2000–2019.
A second investigation (hereinafter called Analysis 2) was, conversely, addressed at verifying if any significant correlation could be recognized between IKC prevalence and some of PMs that showed significant trends in the previous analysis.
On the basis of the IKC disease dataset, prevalence was calculated by Equation (1). IKC yearly prevalence computation relied on the entire chamois population inventory as reported by the local faunistic centers: Corpo Forestale della Valle d’Aosta and hunter committees. Veterinary analyses were performed by IZS PLV SC Aosta and CeRMAS. Faunistic season was assumed to start in September and terminate at the end of the next year August.
Since interaction between breed animals and ungulates is more likely in good quality grass areas, rangelands could represent hotspots for the spreading of the disease.
CLC2018 classes 231 = pastures, 321 = natural grasslands, 322 = moors and heathland were a-priori assumed as favorable land cover classes to test IKC prevalence against to (see Figure 4). They were aggregated into a single macro-class, hereinafter called “FAV” (see Figure 5).
FAV was then disaggregated into 3 sub-classes (FAV1, FAV2, FAV3) according to the altitude ranges of Table 3. Altitude classes (Table 3) were obtained by SRTM reclassification operated by simple table method [57,58] available in SAGA GIS v.7.00 tool [59].
  • Analysis 1: testing PM and ET trends from NTS and ETS
To investigate if and how IKC prevalence could be also related to ongoing changes affecting vegetation activity in consequence of climate change, some metrics were extracted from the abovementioned time series (NTS and ETS). Phenological metrics (PMs) are synthetic descriptors of vegetation activity along its annual growing season. Climate change proved to condition such activity shifting and reshaping past “ordinary” behavior of plants along the year. The following PMs were considered: the start of the growing season (SOS) representing the day of the year (DOY) when phenology is admitted to boost; the end of the season representing the day of the year (DOY) when phenology is admitted to stop; the length of the growing season (LOS) representing the time range (in number of days) separating EOS from SOS; the maximum of NDVI (MAXVI) representing the highest value reached by NDVI during the growing season and proved to be a good predictor of climate change effects on vegetation.
PMs were estimated by TIMESAT 3.3 with STL software [60,61,62] that was specifically developed to enable the monitoring of land surface processes by remotely sensed data. TIMESAT 3.3 with STL [63,64,65] iteratively fits and smooths by mathematical functions the yearly NDVI time-series, finding the best smoothed approximation of the NDVI along the year at pixel level. Once raw data have been approximated by the selected fitting function, PMs can be extracted in correspondence of singular points having a phenological meaning (e.g., EOS, SOS, LOS, MAXVI, etc.) along the local temporal profile.
With reference to ETS and separately for FAV1, FAV2, FAV3, class yearly average ET value was computed and analyzed along the years in the period 2000–2019 at the regional level, to explore eventual trends in ET values, as well.
Since analysis is performed at single pixel and year level, PM estimates are saved as raster layers showing PM spatial distribution at the considered year.
  • Analysis 2: IKC prevalence vs. PMs/ET
After demonstrating that significant trends could be recognized affecting PMs and ET values from MOD13Q and MOD16A2 datasets, respectively, we tested their potential correlation with IKC Pr at the regional level. As far as IKC/ET comparison was concerned, the yearly cumulative ET was calculated. The a priori hypothesis was that a change in environmental conditions could drive to a change in IKC occurrences (Pr value). It is worth to remind that Pr data were supplied aggregated at the regional level, while PMs and ET measures were mapped at the pixel level over the whole area. Consequently, all deductions refer to general trends that could be possibly improved if more distributed and geolocated data of IKC Pr were available.
IKC Pr was tested against all the computed metrics and modeled by a 2nd order polynomial regression.

3. Results

IKC Pr was computed according to the available ground dataset, by Equation (1). Values are reported in Table 4.
  • Analysis 1: testing PM and ET trends from NTS and ETS
In this work, phenological metrics (SOS, EOS, LOS, MAXVI) were estimated at pixel level by TIMESAT 3.3 processing the whole NTS. A Seasonal Trend decomposition by Loess (STL) was adopted to de-trend NTS pixel profile and removing noise. Seasonal component was refined by Savitzky–Golay filtering to reduce, but not removing, remaining local strong variations. The yearly growing season was recognized with the whole multi annual time series using the sinusoidal harmonics approach.
With reference to de-trended/filtered NTS profiles PMs were extracted using the simple thresholding approach. An arbitrary value of 0.5 was set as reference NDVI value to refer SOS and EOS to. Consequently, SOS was assumed as the DOY when NDVI reached the 0.5 threshold value along the ascending part of the phenological yearly bell; EOS was assumed as the DOY when NDVI reached the 0.5 threshold value along the descending part of the phenological yearly bell. LOS map was computed by grid differencing from EOS and SOS maps. MAXVI was found looking for the highest NDVI value between SOS and EOS.
PMs were mapped over the area as raster layers and spatially averaged with respect to FAV1, FAV2, and FAV3 classes (Table 3).
Graphs of Figure 6a–c show that average PM values of FAV1, FAV2, FAV3 significantly changed their values in the last 20 years, with a continuous progressive trend observable for all of them. First order polynomials used to model trends where calibrated excluding those PMs and ET estimates/measures whose residuals (computed with respect to a 1st order polynomial) showed a value higher than mean ± 2 times the standard deviation. Interpretation of results is given in the Discussions section of this paper. In terms of strength of changes, with reference to gain values of the estimated 1st order polynomial models (Table 5), it can be observed that: EOS is averagely delaying of about 2.6 d/y and no significant difference can be observed concerning altitude classes.
Differently, SOS appears to averagely anticipate about 2 d/y up to 2000 m a.s.l. (FAV 1 and FAV2) and about 3 d/y at higher altitudes (FAV 3). Consequently, LOS showed to enlarge about 4.7 d/y at lower altitudes (FAV1 and FAV2) and about 6.5 d/y at higher altitudes (FAV3). According to Borgogno et al. [66], potential accuracy of NDVI measurements is about ±0.02; consequently, estimated yearly variations of MAXVI cannot be considered singularly significant. Nevertheless, the cumulated effects along the entire explored period (2000–2019) showed that MAXVI significantly changed, since the accuracy reference value of ±0.02 was largely overcome. MAXVI variations between 2000 and 2019 appear to be positively higher at lower altitudes (about +0.09) while almost stable as altitude increases. These results find strong evidence in different studies in literature [67,68,69,70,71,72,73,74,75,76].
With reference to ET graphs of Figure 6d, a significant increasing trend can be observed for all FAV classes. According to the values of Table 5, it can be noted that all FAV classes behave similarly; water requirement in a period of 8 days appears to averagely increase of about 0.05 kg·m−2 (about 0.5%) every year. This determined that, in the period 2000–2019 cumulated increment of ET (every 8 days) is around 1 kg·m−2 corresponding to a percentage difference in water requirement of about 8%. This could be possibly explained with reference to the previously demonstrated increasing of both biomass production (MAXVI) and enlargement of the growing season, that, consequently make vegetation needing more water yearly.
  • Analysis 2: IKC prevalence vs. PMs/ET
After demonstrating that some PMs and ET values are currently changing with a significant linear trend, authors tested their correlation with IKC Pr. A 2nd order polynomial was found to well approximate the most of tested relationships. Model calibration was achieved after removing outliers from data. This was obtained excluding all those data that, with respect to the 2nd order polynomial model, showed a percent residual Equation (3) ε > ±100%. After outlier removal, a new calibration was run for the model. Results are reported in Table 6.
ε = | e s t i m a t e d o b s e r v e d | o b s e r v e d · 100
As shown in Table 6 and Figure 7, good correlations were found between IKC Pr and LOS, EOS, and ET. A 2nd order polynomial proved to well fit the relationship. Differently, SOS showed a weaker correlation for all the tested situations.
In particular: ET proved to be a good IKC Pr proxy that appeared to be quite independent from altitude; differently, EOS prediction capability was highly more significant for lower altitudes (FAV1); differently, LOS appeared to represent the main factor conditioning IKC spread at medium-high altitudes (FAV2 and FAV3).
As far as ET is concerned, it is well known that it relates to local micro-climatic conditions. In fact, it is directly impacted by micro-local temperature and humidity that favor vegetation growth and limit soil drought. This occurs mainly at lower altitude (FAV1 and FAV2), determining a higher attractiveness for chamois. At higher altitudes (FAV3), the relationship still persists, but it becomes weaker, probably due to a lower availability of biomass and a shorter phenological season.
It is worth to remind that these results were obtained at the regional level, using highly aggregated data. It is the authors’ intent to investigate further to make their deductions more robust and test the capability of generalization of the proposed prediction models in different areas. Moreover, it would be desirable to assess relationships between other diseases and environmental factors to better face future risks possibly related to zoonosis [77,78,79,80] and its dependence on climate change [61] and loss of biodiversity [81].

4. Discussions

The functional roles of domestic and wild host populations in infectious keratoconjunctivitis (IKC) epidemiology have been extensively discussed claiming a domestic reservoir for the more susceptible wild hosts; in the most of cases all deductions were based on limited data.
With the aim to better assess IKC epidemiology in complex host–pathogen alpine systems, the long-term infectious dynamics and molecular epidemiology of Mycoplasma conjunctivae has been investigated in all host populations from different areas in the Pyrenees and Occidental Alps.
Between the years 2000 and 2019, it was consistently detected in Pyrenean and Alpine chamois (Rupicapra p. pyrenaica) populations, as well as in sheep flocks, and occasionally in mouflon (Ovis aries musimon) from the Pyrenees; statistically associated with ocular clinical signs only in chamois. Chamois populations showed different infection dynamics with low but steady prevalence (4.9%) and significant yearly fluctuations (0.0%–40.0%) between the period 2008–2015 [27,28]. Persistence of specific M. conjunctivae strain clusters in wild host populations is demonstrated for six and nine years. Cross-species transmission between chamois and sheep and chamois and mouflon were also sporadically evidenced. In Switzerland, the chamois affected by IKC was found at altitudes between 550 and 3200 m. The estimated overall mortality was less than 5 per cent, but more than 20 per cent have probably died locally [29]. Host population characteristics and M. conjunctivae strains resulted in different epidemiological scenarios in chamois, ranging from the fading out of the mycoplasma to the epidemic and endemic long-term persistence. These findings highlight the capacity of M. conjunctivae to establish diverse interactions and persist in host populations, also with different transmission conditions. Overall, independent M. conjunctivae sylvatic and domestic cycles occurred at the wildlife-livestock interface in the alpine ecosystems with sheep and chamois as the key host species for each cycle, and mouflon as a spill-over host. Although outbreaks of IKC have been described in Austria, France, Italy, Slovenia, and Switzerland, descriptive studies of the role of environmental patterns and to model the outbreaks on a large scale have often been incomplete, owing to the difficulty of detecting fundamental patterns that can affect IKC spread in chamois that live in remote, inaccessible mountain regions. Under this scenario, the remote sensing techniques and EO data can give certainly a huge hand in the understanding and development of possible forecasting models, as we have tried to do in the present work. With these premises, the present study was intended to explore and propose a method based on free accessible EO data to partially close the above-mentioned knowledge gap.
In Aosta Valley (NW Italy), PMs and ET (as measured from the above mentioned EO data) proved to significantly change their values in the last 20 years, with a continuous progressive trend observable for all of them. In terms of strength of changes, an average delay of EOS was observed by about 2.6 days, independently from the altitude class. SOS proved to averagely anticipate about 2 and 3 days per year at lower (<2000 m) and higher (>2000 m) altitudes, respectively. Consequently, LOS is enlarging by about 4.7 and 6.5 days per year at lower (<2000 m) and higher (>2000 m) altitudes, respectively. While looking at the entire period (2000–2019) MAXVI proved to be significantly changing, showing a positive variation (about +0.09) at lower altitude and no variations at higher one. This can be explained admitting that at lower altitudes, in Aosta Valley, grasslands and pastures are often irrigated. Consequently, farmers can vary water release regimes to face climate change effects (higher temperatures, in particular) with the result of moving forage yields (that NDVI is a predictor of) to higher values.
Differently, where more natural (not managed) systems are located (higher altitudes), the increase of yearly MAXVI can be only related to glacier melting that could compensate the increase of water requirement (as confirmed by the ET analysis) by vegetation. Glaciers are, in fact, dramatically reducing in Aosta Valley. Moreover, another compensating action could come from the surrounding forest areas that have been proved to tolerate summer heatwaves.
With reference to ET, a significant increasing trend was observed, independently from altitude. Eight day water requirement from vegetation appears to averagely increase by about 0.05 kg·m−2 (about 0.5%) every year for a total increase of about 1 kg·m−2 in 20 years (2000–2019), corresponding to a percentage difference in water requirement from vegetation of about 8%. This could be possibly explained by the increasing of biomass production (well represented by MAXVI) and by the enlargement of the growing season, that, consequently make vegetation need more water yearly.
As far as PMs/ET and IKC Pr correlation is concerned, some interesting findings came out.
A 2nd order polynomial model was found to well approximate the most of relationships, making possible to support scenarios generation of IKC spreading for forecasting issues.
ET proved to be a good predictor of IKC Pr, with no significant conditioning by altitude. EOS seems to operate good predictions of Pr at lower altitudes, while LOS at medium-high ones. This probably depends on wildlife dynamics, that in autumn and at the beginning of winter, looks for grassland (food) especially at lower altitudes where, in that period, grass is not covered by snow and is wetter. In the case of SOS, a possible interpretation can rely on the fact that, if vegetative season lasting increases, animals descent from mountain to valley can be postponed, thus inducing a higher probability of interaction between potential guests and sick animals at higher altitudes. This certainly can increase also the probability of exposure of animals to disease.
With reference to ET it is mainly related to Pr at lower altitudes where micro-local temperature and humidity favor vegetation growth and limit soil drought, determining a higher attractiveness for chamois. At higher altitudes, ET capability of predicting Pr becomes weaker probably due to a lower availability of biomass and a shorter phenological season.
Authors are conscious that this work just introduces a new way to manage wildlife health problems and cannot be retained conclusive. In the nearer future, more disaggregated investigations should be done, and other areas possibly considered. Nevertheless, the proposed approach is sufficiently innovative in the context of wildlife veterinary and, we hope, could open a new interesting trend to map wildlife diseases and related zoonosis risk associated with the interaction between wild animals and domestic ones. A radical change is expected also by technicians and institutional subjects in their ordinary procedures for recording and managing ground data. In fact, the greatest limit to expand and more focus this research relied on the format of ground data that could be obtained only aggregated at the regional level with no information concerning the specific place where each analyzed animal was found. We invite all involved players to carefully consider the possibility of georeferencing every ground observation that comes to their laboratories. Georeferencing of ground data is at the basis of an effective and reliable spatial based approach like the one here proposed, where EO data (especially if available over a long-time span) play a crucial role. Anyway, this work proved that spatially based forecasting models can be reasonably calibrated for generating maps of risk concerning wildlife diseases and zoonosis spreading in a certain area. Moreover, it showed that relationships between IKC and PMs/ET are probably chancing in terms of strength; in fact, we demonstrated that all the considered predictors are suffering from a significant change possibly related to the ongoing climate change. Consequently, we expect that future approaches should more properly rely on contemporary data spatially distributed in place of aggregated data temporally distributed like the one we processed for this work.

5. Conclusions

A methodological approach, where geomatics and remote sensing data play a crucial role, can certainly represent a useful “addressing” tool in epidemiology and veterinarian eco-pathology fields of application, especially when studying vector-borne diseases and animal diseases and zoonosis which can affect human beings as well. A major technology transfer, thanks to remote sensing, can still be done in the health sector. A better knowledge of the patho-system dynamics, in particular the relationship between environmental component and disease presence/spreading at a landscape level, can certainly help to aim new studies and think a holistic management of how integrate adaptation, mitigation and prevention. When operative, such approaches could support decision-makers dealing with wildlife and domestic animals management and planning (both at the hunting and pastoral level) and with public health (in the perspective of One Health). The availability of free and global remote data is certainly a valid “systemic” tool for risk analysis that can support ordinary diagnostic techniques, allowing continuous monitoring of the effects that climatic and anthropogenic changes in the Alpine and, in general, mountain and wilderness environments, can cause to animals, biodiversity, and the ecosystems.

Author Contributions

Conceptualization, T.O.; Data curation, T.O., R.O., A.V. and E.B.M.; Formal analysis, T.O. and A.V.; Investigation, T.O.; Methodology, T.O. and E.C.; Project administration, R.O. and E.B.M.; Resources, R.O.; Software, T.O. and E.B.M.; Supervision, E.B.M.; Validation, E.B.M.; Visualization, A.V.; Writing—original draft, T.O.; Writing—review & editing, T.O. and E.B.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research does not have received any kind of funding.

Acknowledgments

This work would not have been possible without the collection of disease and faunistic data and hard work led by the Experimental Zooprophylactic Institute of Piedmont, Liguria and Aosta Valley (IZS PLV), complex structure of Aosta with annexed CeRMAS in particular: Serena Robetto and Cristina Guidetti. The authors also intend to thank for the support and hard work carried out by the Aosta Valley Forestry Corp (Corpo Forestale della Valle d’Aosta) in the collection of sampling data at the hunted wildlife control centers. A thanks also to the Regional Hunting Committee and the Regional Wildlife Office and ASL VdA (Azienda Sanitaria Locale VdA). A remarkable thanks to all colleagues in the GEO4Agri DISAFA Lab and CeRMAS group for the continuous exchange of advice and criticisms and observations about the present work. Last but not least, we would like to thank all those who day after day make it possible to use satellite images for free for scientific purposes in order to face the challenges to better the world in a scientific perspective.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AcronymsDescription
ASLLocal Sanitary Enterprise (translated)
CeRMASNational Reference Center for Wildlife Diseases (Italy)
DISAFADepartment of Agricultural, Forest and Food Sciences
DOYDay of Year
DTMDigital Terrain Model
EO DataEarth Observation Data
EOSEnd of Season
ETEvapotranspiration
ETSEvapotranspiration Time Series
FAVAltimetry band-class
GEEGoogle Earth Engine
IKCInfectious keratoconjunctivitis
IZS PLVIstituto Zooprofilattico Sperimentale Piemonte Liguria e Valle d’Aosta
JPLJet Propulsion Laboratory
LOSLength of Season
LSTLand Surface Temperature
MAXVIMaximum of NDVI
NASANational Aeronautics and Space Administration
NDVINormalized Difference Vegetation Index
NTSNDVI Time Series
PCRPolymerase Chain Reaction
PMsPhenological metrics
PrPrevalence (of a disease)
RFDRegional Forestry Districts
SOSStart of Season
SRTMShuttle Radar Topography Mission
UnitoUniversity of Turin
VDAAosta Valley

References

  1. Hay, S.; Snow, R.; Rogers, D. From Predicting Mosquito Habitat to Malaria Seasons Using Remotely Sensed Data: Practice, Problems and Perspectives. Parasitol. Today 1998, 14, 306–313. [Google Scholar] [CrossRef]
  2. Hendrickx, G.; Biesemans, J.; De Deken, R. The use of GIS in veterinary parasitology. GIS Spat. Anal. Vet. Sci. 2009, 145–176. [Google Scholar] [CrossRef]
  3. Kazmi, S.J.H.; Usery, E.L. Application of remote sensing and gis for the monitoring of diseases: A unique research agenda for geographers. Remote Sens. Rev. 2001, 20, 45–70. [Google Scholar] [CrossRef]
  4. Durr, P.A.; Gatrell, A.C. (Eds.) GIS and spatial analysis in veterinary science. Cabi 2004, 1–33. [Google Scholar]
  5. Housman, I.W.; Chastain, R.A.; Finco, M.V. An Evaluation of Forest Health Insect and Disease Survey Data and Satellite-Based Remote Sensing Forest Change Detection Methods: Case Studies in the United States. Remote Sens. 2018, 10, 1184. [Google Scholar] [CrossRef] [Green Version]
  6. Correia, V.R.D.M.; Carvalho, M.S.; Sabroza, P.C.; Vasconcelos, C.H. Remote sensing as a tool to survey endemic diseases in Brazil. Cad. Saúde Pública 2004, 20, 891–904. [Google Scholar] [CrossRef] [Green Version]
  7. Kiang, R. Malaria Modeling and Surveillance. Benchmark Rep. 2009, 1–5. [Google Scholar]
  8. Anyamba, A.; Chretien, J.-P.; Britch, S.C.; Soebiyanto, R.P.; Small, J.L.; Jepsen, R.; Forshey, B.M.; Sanchez, J.L.; Smith, O.; Odette, M.; et al. Global Disease Outbreaks Associated with the 2015–2016 El Niño Event. Sci. Rep. 2019, 9, 1–14. [Google Scholar] [CrossRef] [Green Version]
  9. Kalluri, S.; Gilruth, P.; Rogers, D.; Szczur, M. Surveillance of Arthropod Vector-Borne Infectious Diseases Using Remote Sensing Techniques: A Review. PLoS Pathog. 2007, 3, e116. [Google Scholar] [CrossRef]
  10. Estrada-Peña, A.; Estrada-Sánchez, A.; De La Fuente, J. A global set of Fourier-transformed remotely sensed covariates for the description of abiotic niche in epidemiological studies of tick vector species. Parasites Vectors 2014, 7, 302. [Google Scholar] [CrossRef] [Green Version]
  11. Wang, J.; Jia, P.; Cuadros, D.F.; Xu, M.; Wang, X.; Guo, W.; Portnov, B.A.; Bao, Y.; Yu, S.; Song, G.; et al. A Remote Sensing Data Based Artificial Neural Network Approach for Predicting Climate-Sensitive Infectious Disease Outbreaks: A Case Study of Human Brucellosis. Remote Sens. 2017, 9, 1018. [Google Scholar] [CrossRef] [Green Version]
  12. Hartemink, N.; Vanwambeke, S.O.; Purse, B.V.; Gilbert, M.; Van Dyck, H. Towards a resource-based habitat approach for spatial modelling of vector-borne disease risks. Biol. Rev. 2014, 90, 1151–1162. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Olivero, J.; Fa, J.E.; Real, R.; Márquez, A.L.; Farfán, M.A.; Vargas, J.M.; Gaveau, D.; Salim, M.A.; Park, D.; Suter, J.; et al. Recent loss of closed forests is associated with Ebola virus disease outbreaks. Sci. Rep. 2017, 7, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Robinson, R.A. Plant Pathosystems; Springer: Berlin/Heidelberg, Germany, 1976; pp. 15–31. [Google Scholar]
  15. Conticini, E.; Frediani, B.; Caro, D. Can atmospheric pollution be considered a co-factor in extremely high level of SARS-CoV-2 lethality in Northern Italy? Environ. Pollut. 2020, 261, 114465. [Google Scholar] [CrossRef]
  16. Rinaldi, L.; Musella, V.; Biggeri, A.; Cringoli, G. New insights into the application of geographical information systems and remote sensing in veterinary parasitology. Geospat. Health 2006, 1, 33. [Google Scholar] [CrossRef] [Green Version]
  17. Jebara, K.B. The role of Geographic Information System (GIS) in the control and prevention of animal diseases. Conf. OIE 2007, 1, 175–183. [Google Scholar]
  18. Barile, M.F.; Del Giudice, R.A.; Tully, J.G. Isolation and Characterization of Mycoplasma conjunctivae sp. n. from Sheep and Goats with Keratoconjunctivitis. Infect. Immun. 1972, 5, 70–76. [Google Scholar] [CrossRef] [Green Version]
  19. Giacometti, M.; Janovsky, M.; Jenny, H.; Nicolet, J.; Belloy, L.; Goldschmidt-Clermont, E.; Frey, J. Mycoplasma conjunctivae infection is not maintained in alpine chamois in eastern switzerland. J. Wildl. Dis. 2002, 38, 297–304. [Google Scholar] [CrossRef] [Green Version]
  20. Giangaspero, M.; Orusa, R.; Nicholas, R.A.; Harasawa, R.; Ayling, R.D.; Churchward, C.; Whatmore, A.M.; Bradley, D.; Robetto, S.; Sacchi, L.; et al. Characterization of mycoplasma isolated from an ibex (capra ibex) suffering from keratoconjunctivitis in northern italy. J. Wildl. Dis. 2010, 46, 1070–1078. [Google Scholar] [CrossRef]
  21. Degiorgis, M.P.; Frey, J.; Nicolet, J.; Abdo, E.M.; Fatzer, R.; Schlatter, Y.; Reist, S.; Janovsky, M.; Giacometti, M. An outbreak of infectious keratoconjunctivitis in Alpine chamois (Rupicapra r. rupicapra) in Simmental-Gruyères. Schweiz Arch Tierheilkd 2000, 142, 520–527. [Google Scholar]
  22. Grattarola, C.; Frem, E.M.; Orusa, R.; Nicolet, M.G. Ker’atoconjunctivitis in. Vet. Rec. 1999, 145, 588–589. [Google Scholar] [CrossRef] [PubMed]
  23. Giacometti, M.; Janovsky, M.; Belloy, L.; Frey, J. Infectious keratoconjunctivitis of ibex, chamois and other Caprinae. Rev. Sci. Tech. Off. Int. Epizoot. 2002, 21, 335–345. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Mavrot, F.; Vilei, E.M.; Marreros, N.; Signer, C.; Frey, J.; Ryser-Degiorgis, M.-P. Occurrence, quantification, and genotyping of Mycoplasma conjunctivae in wild caprinae with and without infectious keratoconjunctivitis. J. Wildl. Dis. 2012, 48, 619–631. [Google Scholar] [CrossRef] [PubMed]
  25. Hars, J.; Gauthier, D. Suivi de l’évolution de la kérato-conjonctivite sur le peuplement d’ongulés sauvages du Parc National de la Vanoise en 1983. Trav. Sci. Parc. Nat. Vanoise 1984, 14, 157–210. [Google Scholar]
  26. Tschopp, R.; Frey, J.; Zimmermann, L.; Giacometti, M. Outbreaks of infectious keratoconjunctivitis in alpine chamois and ibex in Switzerland between 2001 and 2003. Vet. Rec. 2005, 157, 13–18. [Google Scholar] [CrossRef]
  27. Nesti, I.; Posillico, M.; Lovari, S. Ranging behaviour and habitat selection of Alpine chamois. Ethol. Ecol. Evol. 2010, 22, 215–231. [Google Scholar] [CrossRef]
  28. Arnal, M.; Herrero, J.; De La Fe, C.; Revilla, M.; Prada, C.; Martínez-Durán, D.; Gómez-Martin, Á.; Fernández-Arberas, O.; Amores, J.; Contreras, A.; et al. Dynamics of an Infectious Keratoconjunctivitis Outbreak by Mycoplasma conjunctivae on Pyrenean Chamois Rupicapra p. pyrenaica. PLoS ONE 2013, 8, e61887. [Google Scholar] [CrossRef] [Green Version]
  29. Giangaspero, M.; Domenis, L.; Robetto, S.; Orusa, R. Histological and virological findings in severe meningoencephalitis associated with border disease virus in Alpine chamois (Rupicapra rupicapra rupicapra) in Aosta Valley, Italy. Open Vet. J. 2019, 9, 81–87. [Google Scholar] [CrossRef] [Green Version]
  30. Abbona, F.; Venturino, E. An eco-epidemic model for infectious keratoconjunctivitis caused by Mycoplasma conjunctivae in domestic and wild herbivores, with possible vaccination strategies. Math. Methods Appl. Sci. 2016, 41, 2269–2280. [Google Scholar] [CrossRef]
  31. Ambroselli, S. Istat working papers. ISTAT 2019, 24, 47–49. [Google Scholar]
  32. Rubel, F.; Brugger, K.; Haslinger, K.; Auer, I. The climate of the European Alps: Shift of very high resolution Köppen-Geiger climate zones 1800–2100. Meteorol. Z. 2017, 26, 115–125. [Google Scholar] [CrossRef]
  33. Sergio, F.; Pedrini, P. Biodiversity gradients in the Alps: The overriding importance of elevation. Biodivers. Conserv. 2007, 16, 3243–3254. [Google Scholar] [CrossRef] [Green Version]
  34. Fischer, M.; Rudmann-Maurer, K.; Weyand, A.; Stöcklin, J. Agricultural Land Use and Biodiversity in the Alps. Mt. Res. Dev. 2008, 28, 148–155. [Google Scholar] [CrossRef] [Green Version]
  35. Zimmermann, P.; Tasser, E.; Leitinger, G.; Tappeiner, U. Effects of land-use and land-cover pattern on landscape-scale biodiversity in the European Alps. Agric. Ecosyst. Environ. 2010, 139, 13–22. [Google Scholar] [CrossRef]
  36. Balestrieri, A.; Remonti, L.; Ferrari, N.; Ferrari, A.; Valvo, T.L.; Robetto, S.; Orusa, R. Sarcoptic mange in wild carnivores and its co-occurrence with parasitic helminths in the Western Italian Alps. Eur. J. Wildl. Res. 2006, 52, 196–201. [Google Scholar] [CrossRef]
  37. Renna, M.; Ravetto Enri, S.; Probo, M.; Lussiana, C.; Malfatto, V.; Battaglini, L.M.; Lonati, M.; Lombardi, G. Alpine grasslands: Relations among botanical and chemical variables affecting animal product quality. In Proceedings of the 19th Meeting of the FAO CIHEAM Mountain Pastures Network–Mountain Pastures and Livestock Farming Facing Uncertainty: Environmental, Technical and Socio-Economic Challenges, Zaragoza, Spain, 14–16 June 2016. [Google Scholar]
  38. Cavallero, A.; Aceto, P.; Gorlier, A.; Lombardi, G.; Lonati, M.; Martinasso, B.; Tagliatori, C. I Tipi Pastorali Delle Alpi Piemontesi; Alberto Perdisia Editore Divisione Università: Bologna, Italy, 2007; pp. 1–467. [Google Scholar]
  39. Probo, M.; Lonati, M.; Pittarello, M.; Bailey, D.W.; Garbarino, M.; Gorlier, A.; Lombardi, G. Implementation of a rotational grazing system with large paddocks changes the distribution of grazing cattle in the south-western Italian Alps. Rangel. J. 2014, 36, 445–458. [Google Scholar] [CrossRef] [Green Version]
  40. Zimmermann, L.; Jambresic, S.; Giacometti, M.; Frey, J. Specificity of Mycoplasma conjunctivae strains for alpine chamois Rupicapra r. rupicapra. Wildl. Biol. 2008, 14, 118–124. [Google Scholar] [CrossRef]
  41. Vilei, E.M.; Bonvin-Klotz, L.; Zimmermann, L.; Ryser-Degiorgis, M.-P.; Giacometti, M.; Frey, J. Validation and diagnostic efficacy of a TaqMan real-time PCR for the detection of Mycoplasma conjunctivae in the eyes of infected Caprinae. J. Microbiol. Methods 2007, 70, 384–386. [Google Scholar] [CrossRef]
  42. Belloy, L.; Janovsky, M.; Vilei, E.M.; Pilo, P.; Giacometti, M.; Frey, J. Molecular Epidemiology of Mycoplasma conjunctivae in Caprinae: Transmission across Species in Natural Outbreaks. Appl. Environ. Microbiol. 2003, 69, 1913–1919. [Google Scholar] [CrossRef] [Green Version]
  43. Didan, K. MOD13Q1 MODIS/Terra Vegetation Indices 16-Day L3 Global 250m SIN Grid V006 [Data set]. NASA EOSDIS Land Processes DAAC. Available online: https://0-doi-org.brum.beds.ac.uk/10.5067/MODIS/MOD13Q1.006 (accessed on 28 April 2020).
  44. Gorelick, N.; Hancher, M.; Dixon, M.; Ilyushchenko, S.; Thau, D.; Moore, R. Google Earth Engine: Planetary-scale geospatial analysis for everyone. Remote Sens. Environ. 2017, 202, 18–27. [Google Scholar] [CrossRef]
  45. Schafer, R.W. What Is a Savitzky-Golay Filter? [Lecture Notes]. IEEE Signal Process. Mag. 2011, 28, 111–117. [Google Scholar] [CrossRef]
  46. Press, W.H.; Teukolsky, S.A. Savitzky-Golay Smoothing Filters. Comput. Phys. 1990, 4, 669. [Google Scholar] [CrossRef]
  47. O’Connor, A.; Lausten, K.; Okubo, B.; Harris, T. ENVI Services Engine: Earth and planetary image processing for the cloud. Am. Geophys. Union Poster IN21C-1490 2012, 45, 34–49. [Google Scholar]
  48. Running, S.; Mu, Q.; Zhao, M. MOD16A2 MODIS/Terra Net Evapotranspiration 8-Day L4 Global 500m SIN Grid V006 NASA EOSDIS Land Processes DAAC. EGU 2017. [Google Scholar] [CrossRef]
  49. Cai, J.; Liu, Y.; Lei, T.; Pereira, L.S. Estimating reference evapotranspiration with the FAO Penman–Monteith equation using daily weather forecast messages. Agric. For. Meteorol. 2007, 145, 22–35. [Google Scholar] [CrossRef]
  50. Leuning, R.; Zhang, Y.Q.; Rajaud, A.; Cleugh, H.; Tu, K. A simple surface conductance model to estimate regional evaporation using MODIS leaf area index and the Penman-Monteith equation. Water Resour. Res. 2008, 44. [Google Scholar] [CrossRef]
  51. Orusa, T.; Mondino, E.B. Landsat 8 thermal data to support urban management and planning in the climate change era: A case study in Torino area, NW Italy. Remote Sens. Technol. Appl. Urban Environ. IV 2019, 11157. [Google Scholar] [CrossRef]
  52. Farr, T.G.; Rosen, P.A.; Caro, E.; Crippen, R.; Duren, R.; Hensley, S.; Kobrick, M.; Paller, M.; Rodriguez, E.; Roth, L.; et al. The Shuttle Radar Topography Mission. Rev. Geophys. 2007, 45. [Google Scholar] [CrossRef] [Green Version]
  53. Werner, M. Shuttle Radar Topography Mission (SRTM) Mission Overview. Frequenz 2001, 55, 75–79. [Google Scholar] [CrossRef]
  54. Bossard, M.; Feranec, J.; Otahel, J. CORINE land cover technical guide: Addendum 2000. Researchgate 2000, 40, 1–105. [Google Scholar]
  55. Büttner, G.; Feranec, J.; Jaffrain, G.; Mari, L.; Maucha, G.; Soukup, T. The CORINE land cover 2000 project. EARSeL eProc. 2004, 3, 331–346. [Google Scholar]
  56. Available online: https://developers.google.com/earth-engine/datasets/catalog/COPERNICUS_CORINE_V20_100m (accessed on 1 January 2019).
  57. Hammer, Ø.; Harper, D.A.; Ryan, P.D. PAST: Paleontological statistics software package for education and data analysis. Palaeontol. Electron. 2001, 4, 9. [Google Scholar]
  58. QGIS Delopment Team. QGIS geographic information system. Open Source Geospat. Found. Project 2016, 1, 59. [Google Scholar]
  59. Conrad, O.; Bechtel, B.; Bock, M.; Dietrich, H.; Fischer, E.; Gerlitz, L.; Wehberg, J.; Wichmann, V.; Böhner, J. System for Automated Geoscientific Analyses (SAGA) v. 2.1.4. Geosci. Model Dev. 2015, 8, 1991–2007. [Google Scholar] [CrossRef] [Green Version]
  60. Sarvia, F.; De Petris, S.; Mondino, E.B. Multi-scale remote sensing to support insurance policies in agriculture: From mid-term to instantaneous deductions. GIScience Remote Sens. 2020, 57, 770–784. [Google Scholar] [CrossRef]
  61. Jönsson, P.; Eklundh, L. TIMESAT—A program for analyzing time-series of satellite sensor data. Comput. Geosci. 2004, 30, 833–845. [Google Scholar] [CrossRef] [Green Version]
  62. Tan, B.; Morisette, J.T.; Wolfe, R.E.; Gao, F.; Ederer, G.A.; Nightingale, J.; Pedelty, J.A. An Enhanced TIMESAT Algorithm for Estimating Vegetation Phenology Metrics from MODIS Data. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2010, 4, 361–371. [Google Scholar] [CrossRef]
  63. Eklundh, L.; Jönsson, P. A new spatio-temporal smoother for extracting vegetation seasonality with TIMESAT. In Proceedings of the 35th International Symposium on Remote Sensing of Environment, Beijing, China, 22–26 April 2013. [Google Scholar]
  64. Cleveland, R.B.; Cleveland, W.S.; McRae, J.E.; Terpenning, I. STL: A seasonal-trend decomposition. J. Off. Stat. 1990, 6, 3–73. [Google Scholar]
  65. Jia, K.; Liang, S.; Wei, X.; Yao, Y.; Su, Y.; Jiang, B.; Wang, X. Land Cover Classification of Landsat Data with Phenological Features Extracted from Time Series MODIS NDVI Data. Remote Sens. 2014, 6, 11518–11532. [Google Scholar] [CrossRef] [Green Version]
  66. Borgogno-Mondino, E.; Lessio, A.; Gomarasca, M.A. A fast operative method for NDVI uncertainty estimation and its role in vegetation analysis. Eur. J. Remote Sens. 2016, 49, 137–156. [Google Scholar] [CrossRef]
  67. Cremonese, E.; Filippa, G.; Galvagno, M.; Siniscalco, C.; Oddi, L.; Di Cella, U.M.; Migliavacca, M. Heat wave hinders green wave: The impact of climate extreme on the phenology of a mountain grassland. Agric. For. Meteorol. 2017, 247, 320–330. [Google Scholar] [CrossRef]
  68. Migliavacca, M.; Galvagno, M.; Cremonese, E.; Rossini, M.; Meroni, M.; Sonnentag, O.; Cogliati, S.; Manca, G.; Diotri, F.; Busetto, L.; et al. Using digital repeat photography and eddy covariance data to model grassland phenology and photosynthetic CO2 uptake. Agric. For. Meteorol. 2011, 151, 1325–1337. [Google Scholar] [CrossRef]
  69. Julitta, T.; Cremonese, E.; Migliavacca, M.; Colombo, R.; Galvagno, M.; Siniscalco, C.; Rossini, M.; Fava, F.; Cogliati, S.; Di Cella, U.M.; et al. Using digital camera images to analyse snowmelt and phenology of a subalpine grassland. Agric. For. Meteorol. 2014, 198–199, 116–125. [Google Scholar] [CrossRef]
  70. Migliavacca, M.; Reichstein, M.; Richardson, A.D.; Mahecha, M.D.; Cremonese, E.; Delpierre, N.; Galvagno, M.; Law, B.E.; Wohlfahrt, G.; Black, T.A.; et al. Influence of physiological phenology on the seasonal pattern of ecosystem respiration in deciduous forests. Glob. Chang. Biol. 2014, 21, 363–376. [Google Scholar] [CrossRef]
  71. Colombo, R.; Busetto, L.; Fava, F.; Di Mauro, B.; Migliavacca, M.; Cremonese, E.; Galvagno, M.; Rossini, M.; Meroni, M.; Cogliati, S.; et al. Phenological monitoring of grassland and larch in the Alps from Terra and Aqua MODIS images. Ital. J. Remote Sens. 2011, 83–96. [Google Scholar] [CrossRef]
  72. Colombo, R.; Busetto, L.; Migliavacca, M.; Cremonese, E.; Meroni, M.; Galvagno, M.; Rossini, M.; Siniscalco, C.; Di Cella, U.M. On the spatial and temporal variability of Larch phenological cycle in mountainous areas. Ital. J. Remote Sens. 2009, 41, 79–96. [Google Scholar] [CrossRef]
  73. Filippa, G.; Cremonese, E.; Galvagno, M.; Migliavacca, M.; Di Cella, U.M.; Petey, M.; Siniscalco, C. Five years of phenological monitoring in a mountain grassland: Inter-annual patterns and evaluation of the sampling protocol. Int. J. Biometeorol. 2015, 59, 1927–1937. [Google Scholar] [CrossRef]
  74. Diolaiuti, G.A.; Bocchiola, D.; Vagliasindi, M.; D’Agata, C.; Smiraglia, C. The 1975–2005 glacier changes in Aosta Valley (Italy) and the relations with climate evolution. Prog. Phys. Geogr. Earth Environ. 2012, 36, 764–785. [Google Scholar] [CrossRef]
  75. Haeberli, W.; Hoelzle, M.; Paul, F.; Zemp, M. Integrated monitoring of mountain glaciers as key indicators of global climate change: The European Alps. Ann. Glaciol. 2007, 46, 150–160. [Google Scholar] [CrossRef] [Green Version]
  76. Calanca, P.; Roesch, A.; Jasper, K.; Wild, M. Global Warming and the Summertime Evapotranspiration Regime of the Alpine Region. Clim. Chang. 2006, 79, 65–78. [Google Scholar] [CrossRef]
  77. Goulden, M.L.; Bales, R.C. Mountain runoff vulnerability to increased evapotranspiration with vegetation expansion. Proc. Natl. Acad. Sci. USA 2014, 111, 14071–14075. [Google Scholar] [CrossRef] [Green Version]
  78. Steffen, W.; Richardson, K.; Rockström, J.; Cornell, S.E.; Fetzer, I.; Bennett, E.M.; Biggs, R.; Carpenter, S.R.; De Vries, W.; De Wit, C.A.; et al. Planetary boundaries: Guiding human development on a changing planet. Science 2015, 347, 1259855. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Rockström, J.; Steffen, W.; Noone, K.; Persson, Å.; Chapin, F.S.I.; Lambin, E.; Lenton, T.M.; Scheffer, M.; Folke, C.; Schellnhuber, H.J. Planetary Boundaries: Exploring the Safe Operating Space for Humanity. Ecol. Soc. 2009, 14. [Google Scholar] [CrossRef]
  80. O’Neill, D.W.; Fanning, A.L.; Lamb, W.F.; Steinberger, J.K. A good life for all within planetary boundaries. Nat. Sustain. 2018, 1, 88–95. [Google Scholar] [CrossRef] [Green Version]
  81. Cardinale, B.J.; Duffy, J.E.; Gonzalez, A.; Hooper, D.U.; Perrings, C.; Venail, P.; Narwani, A.; Mace, G.M.; Tilman, D.; Wardle, D.A.; et al. Biodiversity loss and its impact on humanity. Nat. Cell Biol. 2012, 486, 59–67. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The possible role of remote sensing in the patho-system as represented by the so-called disease triangle.
Figure 1. The possible role of remote sensing in the patho-system as represented by the so-called disease triangle.
Remotesensing 12 03542 g001
Figure 2. The study area corresponds to the entire Aosta Valley that is located in the Northern West Alps of Italy, close to the border with France and Switzerland.
Figure 2. The study area corresponds to the entire Aosta Valley that is located in the Northern West Alps of Italy, close to the border with France and Switzerland.
Remotesensing 12 03542 g002
Figure 3. Protected areas in Aosta Valley including parks (green) and reserves (orange).
Figure 3. Protected areas in Aosta Valley including parks (green) and reserves (orange).
Remotesensing 12 03542 g003
Figure 4. Map showing spatial distribution of land cover classes (from CORINE land cover map) that were considered for the study: pastures, natural grasslands, and moors/heathland.
Figure 4. Map showing spatial distribution of land cover classes (from CORINE land cover map) that were considered for the study: pastures, natural grasslands, and moors/heathland.
Remotesensing 12 03542 g004
Figure 5. Aosta Valley altitude ranges: FAV1, FAV2, FAV3.
Figure 5. Aosta Valley altitude ranges: FAV1, FAV2, FAV3.
Remotesensing 12 03542 g005
Figure 6. (ac) Class average values of phenological metrics (PMs) (SOS, EOS, LOS, MAXVI) as estimated by TIMESAT 3.3 with Seasonal Trend decomposition by Loess (STL) for classes FAV1, FAV2, and FAV3. (d) Evapotranspiration (ET) average class values computed from ETS for classes FAV1, FAV2 and FAV3. Reported values correspond to the yearly mean value of the 8 days aggregated data from MOD16A2 product (for acronyms explanation, see Abbreviations). In all the graphs, the red-orange points represent outliers.
Figure 6. (ac) Class average values of phenological metrics (PMs) (SOS, EOS, LOS, MAXVI) as estimated by TIMESAT 3.3 with Seasonal Trend decomposition by Loess (STL) for classes FAV1, FAV2, and FAV3. (d) Evapotranspiration (ET) average class values computed from ETS for classes FAV1, FAV2 and FAV3. Reported values correspond to the yearly mean value of the 8 days aggregated data from MOD16A2 product (for acronyms explanation, see Abbreviations). In all the graphs, the red-orange points represent outliers.
Remotesensing 12 03542 g006aRemotesensing 12 03542 g006b
Figure 7. Graphs (ac) show the most significant relationships that were recognized between IKC prevalence and PMs. (a) IKC Pr vs. EOS at FAV1; (b) IKC Pr vs. LOS at FAV2; (c) IKC Pr vs. LOS at FAV3. Graphs (df) show the most significant relationships that were recognized between IKC prevalence and ET (Cumulative yearly ET) values for FAV1, FAV2, and FAV3, respectively. In all the graphs, the red-orange points represent outliers (see Equation (3)).
Figure 7. Graphs (ac) show the most significant relationships that were recognized between IKC prevalence and PMs. (a) IKC Pr vs. EOS at FAV1; (b) IKC Pr vs. LOS at FAV2; (c) IKC Pr vs. LOS at FAV3. Graphs (df) show the most significant relationships that were recognized between IKC prevalence and ET (Cumulative yearly ET) values for FAV1, FAV2, and FAV3, respectively. In all the graphs, the red-orange points represent outliers (see Equation (3)).
Remotesensing 12 03542 g007
Table 1. Aosta Valley altitude ranges and correspondent areas.
Table 1. Aosta Valley altitude ranges and correspondent areas.
Altitude Ranges
(m a.s.l.)
Area (km2)Area (%)
343–5006.60.2
500–1000236.47.2
1000–1500372.711.4
1500–2000669.920.5
2000–2500994.630.5
2500–3000768.323.6
3000–3500176.65.4
3500–481035.51.1
Table 2. Corine Land Cover (CLC) 2018 Level 3 classes and correspondent areas in Aosta Valley. In black, classes that were considered favorable to host wild animals and, consequently, used to test their correlation with infectious keratoconjunctivitis (IKC) prevalence.
Table 2. Corine Land Cover (CLC) 2018 Level 3 classes and correspondent areas in Aosta Valley. In black, classes that were considered favorable to host wild animals and, consequently, used to test their correlation with infectious keratoconjunctivitis (IKC) prevalence.
CLC2018 Class CodeDescriptionArea
(km2)
Area
(%)
CLC2018 Class CodeDescriptionArea
(km2)
Area
(%)
111Continuous urban fabric1.560.05311Broad-leaved forest58.121.78
112Discontinuous urban fabric35.271.08312Coniferous forest577.9817.71
121Industrial or commercial units8.720.27313Mixed forest104.413.20
122Road and rail networks and associated land0.250.01321Natural grasslands86.042.64
124Airports0.420.01322Moors and heathland106.293.26
131Mineral extraction sites0.660.02324Transitional woodland-shrub424.8413.02
132Dump sites0.270.01332Bare rocks652.6120.00
212Permanently irrigated land0.270.01333Sparsely vegetated areas804.7824.67
221Vineyards3.570.11335Glaciers and permanent snow129.563.97
222Fruit trees and berry plantations2.170.07411Inland marshes0.540.02
231Pastures94.062.88511River0.170.01
242Complex cultivation patterns18.610.57512Lakes3.240.10
243Land principally occupied by agriculture, with significant areas of natural vegetation148.424.55
Table 3. Aosta Valley altitude ranges considered for this work to test IKC prevalence against land cover classes.
Table 3. Aosta Valley altitude ranges considered for this work to test IKC prevalence against land cover classes.
Class CodeAltitude Range (m)
FAV1<1000
FAV21000–2000
FAV32000–3000
Table 4. IKC prevalence in the Aosta Valley autonomous region (whole region). Data were available for the years 2009–2019.
Table 4. IKC prevalence in the Aosta Valley autonomous region (whole region). Data were available for the years 2009–2019.
Year 1Year Analyzed 2IKC Disease Prevalence (%)Number of Samples AnalyzedPositive to IKC
2009–201020102.03026
2010–201120114.71919
2011–201220122.01503
2012–201320135.11588
2013–201420141.11902
2014–201520152.61524
2015–2016201616.415926
2016–201720176.11147
2017–201820187.41088
2018–201920190.01000
1,2 Faunistic season is intended from September (previous year) to August (next year); consequently, IKC disease prevalence values were re-assigned at year level assuming the period January–December as reference time to perform analysis 2. It is worth to remind that positive IKC column refers to chamois that have been detected positive to the IKC considering the overall number of samples analyzed per each year reported in the fourth column after performing a PCR analysis (please see Veterinary Ground Samples Section).
Table 5. Yearly average and cumulated (2000–2019) variations of PMs and ET along the years as estimated by 1st order polynomial regression (Figure 6). In Table 5, DOY refers to the day along the considered year and ranging between 0 and 365. Positive values indicate a delay or lengthening of season, negative values an advance of season according to the observed trends reported in Figure 6. The first part of the table indicates the yearly gain while the second part (cumulated) represent the overall changes in the observed period analyzed for each PMs from the starting year 2000 till arriving to the end in the 2019.
Table 5. Yearly average and cumulated (2000–2019) variations of PMs and ET along the years as estimated by 1st order polynomial regression (Figure 6). In Table 5, DOY refers to the day along the considered year and ranging between 0 and 365. Positive values indicate a delay or lengthening of season, negative values an advance of season according to the observed trends reported in Figure 6. The first part of the table indicates the yearly gain while the second part (cumulated) represent the overall changes in the observed period analyzed for each PMs from the starting year 2000 till arriving to the end in the 2019.
YearlySOS
(DOY)
EOS
(DOY)
LOS
(n. of Days)
MAXVIET
(kg·m−2·8 d−1)
FAV1−2.042.644.700.0050.06
FAV2−2.092.594.810.0040.04
FAV3−3.112.596.400.0030.06
Cumulated 2000–2019SOS
(DOY)
EOS
(DOY)
LOS
(n. of Days)
MAXVIET
(Kg m−2)
FAV1−38.7650.1689.340.0891.14
FAV2−39.7749.2991.410.0720.78
FAV3−59.1749.26121.690.0571.11
Table 6. Coefficient of determination (R2) and p-value resulting after modeling relationships between PMs/ET and IKC Pr at altitude class level. A 2nd order polynomial was used for modeling. (*) indicates statistically significant cases (* p < 0.05).
Table 6. Coefficient of determination (R2) and p-value resulting after modeling relationships between PMs/ET and IKC Pr at altitude class level. A 2nd order polynomial was used for modeling. (*) indicates statistically significant cases (* p < 0.05).
ClassesSOSEOSLOSMAXVIYearly Cumulative ET
R2p-ValueR2p-ValueR2p-ValueR2p-ValueR2p-Value
FAV10.060.830.850.003 *0.260.410.060.900.890.004 *
FAV20.510.110.130.650.760.01 *0.160.680.800.005 *
FAV30.460.160.280.350.710.02 * 0.310.320.650.003 *
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Orusa, T.; Orusa, R.; Viani, A.; Carella, E.; Borgogno Mondino, E. Geomatics and EO Data to Support Wildlife Diseases Assessment at Landscape Level: A Pilot Experience to Map Infectious Keratoconjunctivitis in Chamois and Phenological Trends in Aosta Valley (NW Italy). Remote Sens. 2020, 12, 3542. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12213542

AMA Style

Orusa T, Orusa R, Viani A, Carella E, Borgogno Mondino E. Geomatics and EO Data to Support Wildlife Diseases Assessment at Landscape Level: A Pilot Experience to Map Infectious Keratoconjunctivitis in Chamois and Phenological Trends in Aosta Valley (NW Italy). Remote Sensing. 2020; 12(21):3542. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12213542

Chicago/Turabian Style

Orusa, Tommaso, Riccardo Orusa, Annalisa Viani, Emanuele Carella, and Enrico Borgogno Mondino. 2020. "Geomatics and EO Data to Support Wildlife Diseases Assessment at Landscape Level: A Pilot Experience to Map Infectious Keratoconjunctivitis in Chamois and Phenological Trends in Aosta Valley (NW Italy)" Remote Sensing 12, no. 21: 3542. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12213542

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop