Sustainability in our food and fiber agriculture systems is inherently knowledge intensive. It is more likely to be achieved by using all the knowledge, technology, and resources available, including data-driven agricultural technology and precision agriculture methods, than by relying entirely on human powers of observation, analysis, and memory following practical experience. Data collected by sensors and digested by artificial intelligence (AI) can help farmers learn about synergies between the domains of natural systems that are key to simultaneously achieve sustainability and food security. In the quest for agricultural sustainability, some high-payoff research areas are suggested to resolve critical legal and technical barriers as well as economic and social constraints. These include: the development of holistic decision-making systems, automated animal intake measurement, low-cost environmental sensors, robot obstacle avoidance, integrating remote sensing with crop and pasture models, extension methods for data-driven agriculture, methods for exploiting naturally occurring Genotype x Environment x Management experiments, innovation in business models for data sharing and data regulation reinforcing trust. Public funding for research is needed in several critical areas identified in this paper to enable sustainable agriculture and innovation.

Invasive plant species (IPS) are the second biggest threat to biodiversity after habitat loss. Since the spatial extent of IPS is essential for managing the invaded ecosystem, the current study aims at identifying and mapping the aggressive IPS of Acacia salicina and Acacia saligna, to understand better the key factors influencing their distribution in the coastal plain of Israel. This goal was achieved by integrating airborne-derived hyperspectral imaging and multispectral earth observation for creating species distribution maps. Hyperspectral data, in conjunction with high spatial resolution species distribution maps, were used to train the multispectral images at the species level. We incorporated a series of statistical models to classify the IPS location and to recognize their distribution and density. We took advantage of the phenological flowering stages of Acacia trees, as obtained by the multispectral images, for the support vector machine classification procedure. The classification yielded an overall Kappa coefficient accuracy of 0.89. We studied the effect of various environmental and human factors on IPS density by using a random forest machine learning model, to understand the mechanisms underlying successful invasions, and to assess where IPS have a higher likelihood of occurring. This algorithm revealed that the high density of Acacia most closely related to elevation, temperature pattern, and distances from rivers, settlements, and roads. Our results demonstrate how the integration of remote-sensing data with different data sources can assist in determining IPS proliferation and provide detailed geographic information for conservation and management efforts to prevent their future spread.

Climate and land-use change profoundly affect plant species distribution (SD) and composition, and the impact of these processes is expected to increase in the coming years. As a proxy of global changes, knowledge of SD and diversity along climatic gradients is essential to determine the efforts needed for species conservation. Plant spectral diversity is an emerging approach used as a proxy for species diversity based on remote sensing. Thus, the research aim was to develop a comprehensive methodology based on spectral diversity for SD and richness mapping and to study their relations with environmental and human-derived factors, demonstrated along Mediterranean to semi-arid climatic gradient. The study addresses two main knowledge gaps regarding spectral diversity: (1) improving the accuracy of woody species classification by features extraction and selection, and by using texture analysis in an ecosystem characterized by high spatial variability and relatively small-sized and sparse woody vegetation; and (2) developing a better estimate of the local species ‎richness and their response to environmental and human-derived factors (i.e. climate, topography, substrate, and land cover factors) across a transition zone between Mediterranean woodlands and semi-arid dwarf shrublands. A hyperspectral image was acquired for a 43-km strip along the study area using an airborne flight of AISA-FENIX (380–2500 nm, 420 bands) at the end of the 2017 rainy season. The dominant species were surveyed, with a total number of 247 trees and shrubs, to train a machine learning support vector machine (SVM) classification for species distribution mapping, which yielded an overall accuracy of 86.1%. A feature extraction and selection methodology was developed, combining principal component analysis and neighborhood component analysis techniques, facilitating the identification of 33 spectral diagnostic bands out of 330 spectral bands. The classification accuracy was decreased by about 2% to 84.2% using only 33 spectral bands. The classification accuracy improved by about 7.1% for the seven large crown species (93.3%) by adding texture information. Later, the local species richness was calculated by utilizing the alpha diversity index (i.e. the Shannon Index) for 30-m grid cells and was tested in response to environmental (i.e. climate, substrate, and topography) and human-derived factors (i.e. land cover). The highest sensitivity to alpha diversity factors was mean annual precipitation, slope, and land surface temperature. The alpha diversity showed higher richness in the natural Mediterranean shrubland and the guarrigue located in the northern part of the climate gradient. We suggest that the approach presented here significantly improves the estimation of woody species distribution and diversity in areas characterized by high spatial heterogeneity along steep climatic gradients.

Water and energy are recognized as the most influential climatic vegetation growth-limiting factors. These factors are usually measured from ground meteorological stations. However, since both vary in space, time, and scale, they can be assessed by satellite-derived biophysical indicators. Energy, represented by land surface temperature (LST), is assumed to resemble air temperature; and water availability, related to precipitation, is represented by the normalized difference vegetation index (NDVI). It is hypothesized that positive correlations between LST and NDVI indicate energy-limited conditions, while negative correlations indicate water-limited conditions. The current project aimed to quantify the spatial and seasonal (spring and summer) distributions of LST–NDVI relations over Europe, using long-term (2000–2017) MODIS images. Overlaying the LST–NDVI relations on the European biome map revealed that relations between LST and NDVI were highly diverse among the various biomes and throughout the entire study period (March–August). During the spring season (March–May), 80% of the European domain, across all biomes, showed the dominance of significant positive relations. However, during the summer season (June–August), most of the biomes—except the northern ones—turned to negative correlation. This study demonstrates that the drought/vegetation/stress spectral indices, based on the prevalent hypothesis of an inverse LST–NDVI correlation, are spatially and temporally dependent. These negative correlations are not valid in regions where energy is the limiting factor (e.g., in the drier regions in the southern and eastern extents of the domain) or during specific periods of the year (e.g., the spring season). Consequently, it is essential to re-examine this assumption and restrict applications of such an approach only to areas and periods in which negative correlations are observed. Predicted climate change will lead to an increase in temperature in the coming decades (i.e., increased LST), as well as a complex pattern of precipitation changes (i.e., changes of NDVI). Thus shifts in plant species locations are expected to cause a redistribution of biomes.

Field spectroscopy has been suggested to be an efficient method for predicting soil properties using quantitative mathematical models in a rapid and non-destructive manner. Traditional multivariate regression algorithms usually regard the modeling of each soil property as a single task, which means only one response variable is considered as the output during modeling. Therefore, these algorithms are less suitable for the prediction of several key soil properties with low concentrations or unobvious spectral absorption signals. In the current study, we investigated the performance of a linear multi-task learning (LMTL) algorithm based on a regularized dirty model for modeling and predicting several key soil properties using field spectroscopy (350–2500 nm) as an integrated approach. We tested seven key soil properties including available nitrogen (N), phosphorus (P) and potassium (K), pH, water content (WC), organic matter (OM), and electrical conductivity (EC) in drylands. The model performances of LMTL models were compared with the commonly used single-task algorithm of the partial least squares regression (PLS-R). Our results show that the LMTL models outperformed the PLS-R models with the advantage of shared features; the ratio of performance to deviation (RPD) values in the validation set improved by 10.24%, 4.93%, 25.77%, 11.76%, 6.74%, 53.13%, and 3.15% for N, P, K, pH, WC, OM, and EC, respectively. The best prediction was obtained for OM with RPD = 2.29, indicating high accuracy (RPD > 2). The prediction results of N, P, WC, and pH were categorized as of moderate accuracy (1.4 < RPD < 2), while K and EC were categorized as of poor accuracy (RPD < 1.4). However, the explanatory power of the LMTL models was moderate due to fewer features being selected by the regularization algorithm of the LMTL approach, which should be further studied in the soil spectral analysis. Our results highlight the use of LMTL in field spectroscopy analysis that can improve the generalization performance of regression models for predicting soil properties.