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.

Google Earth Engine has deeply changed the way in which Earth observation data are processed, allowing the analysis of wide areas in a faster and more efficient way than ever before. Since its inception, many functions have been implemented by a rapidly expanding community, but none so far has focused on the computation of phenological metrics in mountain areas with high-resolution data. This work aimed to fill this gap by developing an open-source Google Earth Engine algorithm to map phenological metrics (PMs) such as the Start of Season, End of Season, and Length of Season and detect the Peak of Season in mountain areas worldwide using high-resolution free satellite data from the Landsat collection and Sentinel-2. The script was tested considering the entire Alpine chain. The validation was performed by the cross-computation of PMs using the R package greenbrown, which permits land surface phenology and trend analysis, and the Moderate-Resolution Imaging Spectroradiometer (MODIS) in homogeneous quote and land cover alpine landscapes. MAE and RMSE were computed. Therefore, this algorithm permits one to compute with a certain robustness PMs retrieved from higher-resolution free EO data from GEE in mountain areas worldwide.

Satellite remote sensing is a power tool for the long-term monitoring of vegetation. This work, with reference to a regional case study, investigates remote sensing potentialities for describing the annual phenology of rangelands and broad-leaved forests at the landscape level with the aim of detecting eventual effects of climate change in the Alpine region of the Aosta Valley (Northwest (NW) Italy). A first analysis was aimed at estimating phenological metrics (PMs) from satellite images time series and testing the presence of trends along time. A further investigation concerned evapotranspiration from vegetation (ET) and its variation along the years. Additionally, in both the cases the following meteorological patterns were considered: air temperature anomalies, precipitation trends and the timing of yearly seasonal snow melt. The analysis was based on the time series (TS) of different MODIS collections datasets together with Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS) collection obtained through Google Earth Engine. Ground weather stations data from the Centro Funzionale VdA ranging from 2000 to 2019 were used. In particular, the MOD13Q1 v.6, MOD16A2 and MOD10A1 v.6 collections were used to derive PMs, ET and snow cover maps. The SRTM (shuttle radar topography mission) DTM (digital terrain model) was also used to describe local topography while the Coordination of Information on the Environment (CORINE) land cover map was adopted to investigate land use classes. Averagely in the area, rangelands and broad-leaved forests showed that the length of season is getting longer, with a general advance of the SOS (start of the season) and a delay in the EOS (end of the season). With reference to ET, significant increasing trends were generally observed. The water requirement from vegetation appeared to have averagely risen about 0.05 Kg·m⁻² (about 0.5%) per year in the period 2000–2019, for a total increase of about 1 Kg·m⁻² in 20 years (corresponding to a percentage difference in water requirement from vegetation of about 8%). This aspect can be particularly relevant in the bottom of the central valley, where the precipitations have shown a statistically significant decreasing trend in the period 2000–2019 (conversely, no significant variation was found in the whole territory). Additionally, the snowpack timing persistence showed a general reduction trend. PMs and ET and air temperature anomalies, as well as snow cover melting, proved to have significantly changed their values in the last 20 years, with a continuous progressive trend. The results encourage the adoption of remote sensing to monitor climate change effects on alpine vegetation, with particular focus on the relationship between phenology and other abiotic factors permitting an effective technological transfer.

Changes in land use and land cover as well as feedback on the climate deeply affect the landscape worldwide. This phenomenon has also enlarged the human–wildlife interface and amplified the risk of potential new zoonoses. The expansion of the human settlement is supposed to affect the spread and distribution of wildlife diseases such as canine distemper virus (CDV), by shaping the distribution, density, and movements of wildlife. Nevertheless, there is very little evidence in the scientific literature on how remote sensing and GIS tools may help the veterinary sector to better monitor the spread of CDV in wildlife and to enforce ecological studies and new management policies in the near future. Thus, we perform a study in Northwestern Italy (Aosta Valley Autonomous Region), focusing on the relative epidemic waves of CDV that cause a virulent disease infecting different animal species with high host mortality. CDV has been detected in several mammalian from Canidae, Mustelidae, Procyonidae, Ursidae, and Viverridae families. In this study, the prevalence is determined at 60% in red fox (Vulpes vulpes, n = 296), 14% in wolf (Canis lupus, n = 157), 47% in badger (Meles meles, n = 103), and 51% in beech marten (Martes foina, n = 51). The detection of CDV is performed by means of real-time PCR. All the analyses are done using the TaqMan approach, targeting the chromosomal gene for phosphoprotein, gene P, that is involved in the transcription and replication of the virus. By adopting Earth Observation Data, we notice that CDV trends are strongly related to an altitude gradient and NDVI entropy changes through the years. A tentative model is developed concerning the ground data collected in the Aosta Valley region. According to our preliminary study, entropy computed from remote-sensing data can represent a valuable tool to monitor CDV spread as a proxy data predictor of the intensity of fragmentation of a given landscape and therefore also to monitor CDV. In conclusion, the evaluation from space of the landscape variations regarding the wildlife ecological corridors due to anthropic or natural disturbances may assist veterinarians and wildlife ecologists to enforce management health policies in a One Health perspective by pointing out the time and spatial conditions of interaction between wildlife. Surveillance and disease control actions are supposed to be carried out to strengthen the usage of geospatial analysis tools and techniques. These tools and techniques can deeply assist in better understanding and monitoring diseases affecting wildlife thanks to an integrated management approach.

Earth observation data have assumed a key role in environmental monitoring, as well as in risk assessment. Rising temperatures and consequently heat waves due to ongoing climate change represent an important risk considering the population, as well as animals, exposed. This study was focused on the Aosta Valley Region in NW Italy. To assess population exposure to these patterns, the following datasets have been considered: (1) HDX Meta population dataset refined and updated in order to map population distribution and its features; (2) Landsat collection (missions 4 to 9) from 1984 to 2022 obtained and calibrated in Google Earth Engine to model LST trends. A pixel-based analysis was performed considering Aosta Valley settlements and relative population distribution according to the Meta population dataset. From Landsat data, LST trends were modelled. The LST gains computed were used to produce risk exposure maps considering the population distribution and structure (such as ages, gender, etc.). To check the consistency and quality of the HDX population dataset, MAE was computed considering the ISTAT population dataset at the municipality level. Exposure-risk maps were finally realized adopting two different approaches. The first one considers only LST gain maximum by performing an ISODATA unsupervised classification clustering in which the separability of each class obtained and was checked by computing the Jeffries–Matusita (J-M) distances. The second one was to map the rising temperature exposure by developing and performing a risk geo-analysis. In this last case the input parameters considered were defined after performing a multivariate regression in which LST maximum was correlated and tested considering (a) Fractional Vegetation Cover (FVC), (b) Quote, (c) Slope, (d) Aspect, (e) Potential Incoming Solar Radiation (mean sunlight duration in the meteorological summer season), and (f) LST gain mean. Results show a steeper increase in LST maximum trend, especially in the bottom valley municipalities, and especially in new built-up areas, where more than 60% of the Aosta Valley population and domestic animals live and where a high exposure has been detected and mapped with both approaches performed. Maps produced may help the local planners and the civil protection services to face global warming from a One Health perspective.

Land cover (LC) maps are crucial to environmental modeling and define sustainable management and planning policies. The development of a land cover mapping continuous service according to the new EAGLE legend criteria has become of great interest to the public sector. In this work, a tentative approach to map land cover overcoming remote sensing (RS) limitations in the mountains according to the newest EAGLE guidelines was proposed. In order to reach this goal, the methodology has been developed in Aosta Valley, NW of Italy, due to its higher degree of geomorphological complexity. Copernicus Sentinel-1 and 2 data were adopted, exploiting the maximum potentialities and limits of both, and processed in Google Earth Engine and SNAP. Due to SAR geometrical distortions, these data were used only to refine the mapping of urban and water surfaces, while for other classes, composite and timeseries filtered and regularized stack from Sentinel-2 were used. GNSS ground truth data were adopted, with training and validation sets. Results showed that K-Nearest-Neighbor and Minimum Distance classification permit maximizing the accuracy and reducing errors. Therefore, a mixed hierarchical approach seems to be the best solution to create LC in mountain areas and strengthen local environmental modeling concerning land cover mapping.