Universidad Autónoma de Nayarit

González-Montaña J., Escalera-Valente F., Alonso A.J., Lomillos J.M., Robles R., Alonso M.E.

Cobalt, as a trace element, is essential for rumen microorganisms for the formation of vitamin B12. In the metabolism of mammals, vitamin B12 is an essential part of two enzymatic systems involved in multiple metabolic reactions, such as in the metabolism of carbohydrates, lipids, some amino acids and DNA. Adenosylcobalamin and methylcobalamin are coenzymes of methylmalonyl coenzyme A (CoA) mutase and methionine synthetase and are essential for obtaining energy through ruminal metabolism. Signs of cobalt deficiency range from hyporexia, reduced growth and weight loss to liver steatosis, anemia, impaired immune function, impaired reproductive function and even death. Cobalt status in ruminant animals can be assessed by direct measurement of blood or tissue concentrations of cobalt or vitamin B12, as well as the level of methylmalonic acid, homocysteine or transcobalamin in blood; methylmalonic acid in urine; some variables hematological; food consumption or growth of animals. In general, it is assumed that the requirement for cobalt (Co) is expressed around 0.11 ppm (mg/kg) in the dry matter (DM) diet; current recommendations seem to advise increasing Co supplementation and placing it around 0.20 mg Co/kg DM. Although there is no unanimous criterion about milk production, fattening or reproductive rates in response to increased supplementation with Co, in some investigations, when the total Co of the diet was approximately 1 to 1.3 ppm (mg/kg), maximum responses were observed in the milk production.

Verrelst J., Rivera-Caicedo J.P., Reyes-Muñoz P., Morata M., Amin E., Tagliabue G., Panigada C., Hank T., Berger K.

Satellite imaging spectroscopy for terrestrial applications is reaching maturity with recently launched and upcoming science-driven missions, e.g. PRecursore IperSpettrale della Missione Applicativa (PRISMA) and Environmental Mapping and Analysis Program (EnMAP), respectively. Moreover, the high-priority mission candidate Copernicus Hyperspectral Imaging Mission for the Environment (CHIME) is expected to globally provide routine hyperspectral observations to support new and enhanced services for, among others, sustainable agricultural and biodiversity management. Thanks to the provision of contiguous visible-to-shortwave infrared spectral data, hyperspectral missions open enhanced opportunities for the development of new-generation retrieval models of multiple vegetation traits. Among these, canopy nitrogen content (CNC) is one of the most promising variables given its importance for agricultural monitoring applications. This work presents the first hybrid CNC retrieval model for the operational delivery from spaceborne imaging spectroscopy data. To achieve this, physically-based models were combined with machine learning regression algorithms and active learning (AL). The key concepts involve: (1) coupling the radiative transfer models PROSPECT-PRO and SAIL for the generation of a wide range of vegetation states as training data, (2) using dimensionality reduction to deal with collinearity, (3) applying an AL technique in combination with Gaussian process regression (GPR) for fine-tuning the training dataset on in field collected data, and (4) adding non-vegetated spectra to enable the model to deal with spectral heterogeneity in the image. The final CNC model was successfully validated against field data achieving a low root mean square error (RMSE) of 3.4 g / m 2 and coefficient of determination ( R 2 ) of 0.7. The model was applied to a PRISMA image acquired over agricultural areas in the North of Munich, Germany. Mapping aboveground CNC yielded reliable estimates over the whole landscape and meaningful associated uncertainties. These promising results demonstrate the feasibility of routinely quantifying CNC from space, such as in an operational context as part of the future CHIME mission.

Frías-Espericueta M.G., Bautista-Covarrubias J.C., Osuna-Martínez C.C., Delgado-Alvarez C., Bojórquez C., Aguilar-Juárez M., Roos-Muñoz S., Osuna-López I., Páez-Osuna F.

The objective of this review is to synthetize knowledge of the relationship between metals and oxidative stress in aquatic crustaceans (mainly shrimp and crabs) to analyze antioxidant responses when organisms are exposed to metals because the direct metal binding to the active site of enzymes inactivates most of the antioxidant systems. This study reviewed over 150 works, which evidenced that: (i) antioxidant defense strategies used by aquatic decapod crustaceans vary among species; (ii) antioxidant enzymes could be induced or inhibited by metals depending on species, concentration, and exposure time; and (iii) some antioxidant enzymes, as superoxide dismutase increase their activity in low metal levels and time exposures, but their activities are inhibited with higher metal concentrations and exposure time.

Tagliabue G., Boschetti M., Bramati G., Candiani G., Colombo R., Nutini F., Pompilio L., Rivera-Caicedo J.P., Rossi M., Rossini M., Verrelst J., Panigada C.

The recently launched and upcoming hyperspectral satellite missions, featuring contiguous visible-to-shortwave infrared spectral information, are opening unprecedented opportunities for the retrieval of a broad set of vegetation traits with enhanced accuracy through novel retrieval schemes. In this framework, we exploited hyperspectral data cubes collected by the new-generation PRecursore IperSpettrale della Missione Applicativa (PRISMA) satellite of the Italian Space Agency to develop and test a hybrid retrieval workflow for crop trait mapping. Crop traits were mapped over an agricultural area in north-east Italy (Jolanda di Savoia, FE) using PRISMA images collected during the 2020 and 2021 vegetative seasons. Leaf chlorophyll content, leaf nitrogen content, leaf water content and the corresponding canopy level traits scaled through leaf area index were estimated using a hybrid retrieval scheme based on PROSAIL-PRO radiative transfer simulations coupled with a Gaussian processes regression algorithm. Active learning algorithms were used to optimise the initial set of simulated data by extracting only the most informative samples. The accuracy of the proposed retrieval scheme was evaluated against a broad ground dataset collected in 2020 in correspondence of three PRISMA overpasses. The results obtained were positive for all the investigated variables. At the leaf level, the highest accuracy was obtained for leaf nitrogen content (LNC: r2=0.87, nRMSE=7.5%), while slightly worse results were achieved for leaf chlorophyll content (LCC: r2=0.67, nRMSE=11.7%) and leaf water content (LWC: r2=0.63, nRMSE=17.1%). At the canopy level, a significantly higher accuracy was observed for nitrogen content (CNC: r2=0.92, nRMSE=5.5%) and chlorophyll content (CCC: r2=0.82, nRMSE=10.2%), whereas comparable results were obtained for water content (CWC: r2=0.61, nRMSE=16%). The developed models were additionally tested against an independent dataset collected in 2021 to evaluate their robustness and exportability. The results obtained (i. e., LCC: r2=0.62, nRMSE=27.9%; LNC: r2=0.35, nRMSE=28.4%; LWC: r2=0.74, nRMSE=20.4%; LAI: r2=0.84, nRMSE=14.5%; CCC: r2=0.79, nRMSE=18.5%; CNC: r2=0.62, nRMSE=23.7%; CWC: r2=0.92, nRMSE=16.6%) evidence the transferability of the hybrid approach optimised through active learning for most of the investigated traits. The developed models were then used to map the spatial and temporal variability of the crop traits from the PRISMA images. The high accuracy and consistency of the results demonstrates the potential of spaceborne imaging spectroscopy for crop monitoring, paving the path towards routine retrievals of multiple crop traits over large areas that could drive more effective and sustainable agricultural practices worldwide.

Berger K., Rivera Caicedo J.P., Martino L., Wocher M., Hank T., Verrelst J.

The current exponential increase of spatiotemporally explicit data streams from satellite-based Earth observation missions offers promising opportunities for global vegetation monitoring. Intelligent sampling through active learning (AL) heuristics provides a pathway for fast inference of essential vegetation variables by means of hybrid retrieval approaches, i.e., machine learning regression algorithms trained by radiative transfer model (RTM) simulations. In this study we summarize AL theory and perform a brief systematic literature survey about AL heuristics used in the context of Earth observation regression problems over terrestrial targets. Across all relevant studies it appeared that: (i) retrieval accuracy of AL-optimized training data sets outperformed models trained over large randomly sampled data sets, and (ii) Euclidean distance-based (EBD) diversity method tends to be the most efficient AL technique in terms of accuracy and computational demand. Additionally, a case study is presented based on experimental data employing both uncertainty and diversity AL criteria. Hereby, a a simulated training data base by the PROSAIL-PRO canopy RTM is used to demonstrate the benefit of AL techniques for the estimation of total leaf carotenoid content (Cxc) and leaf water content (Cw). Gaussian process regression (GPR) was incorporated to minimize and optimize the training data set with AL. Training the GPR algorithm on optimally AL-based sampled data sets led to improved variable retrievals compared to training on full data pools, which is further demonstrated on a mapping example. From these findings we can recommend the use of AL-based sub-sampling procedures to select the most informative samples out of large training data pools. This will not only optimize regression accuracy due to exclusion of redundant information, but also speed up processing time and reduce final model size of kernel-based machine learning regression algorithms, such as GPR. With this study we want to encourage further testing and implementation of AL sampling methods for hybrid retrieval workflows. AL can contribute to the solution of regression problems within the framework of operational vegetation monitoring using satellite imaging spectroscopy data, and may strongly facilitate data processing for cloud-computing platforms.

Estévez J., Vicent J., Rivera-Caicedo J.P., Morcillo-Pallarés P., Vuolo F., Sabater N., Camps-Valls G., Moreno J., Verrelst J.

Retrieval of vegetation properties from satellite and airborne optical data usually takes place after atmospheric correction, yet it is also possible to develop retrieval algorithms directly from top-of-atmosphere (TOA) radiance data. One of the key vegetation variables that can be retrieved from at-sensor TOA radiance data is leaf area index (LAI) if algorithms account for variability in atmosphere. We demonstrate the feasibility of LAI retrieval from Sentinel-2 (S2) TOA radiance data (L1C product) in a hybrid machine learning framework. To achieve this, the coupled leaf-canopy-atmosphere radiative transfer models PROSAIL-6SV were used to simulate a look-up table (LUT) of TOA radiance data and associated input variables. This LUT was then used to train the Bayesian machine learning algorithms Gaussian processes regression (GPR) and variational heteroscedastic GPR (VHGPR). PROSAIL simulations were also used to train GPR and VHGPR models for LAI retrieval from S2 images at bottom-of-atmosphere (BOA) level (L2A product) for comparison purposes. The BOA and TOA LAI products were consistently validated against a field dataset with GPR (R2 of 0.78) and with VHGPR (R2 of 0.80) and for both cases a slightly lower RMSE for the TOA LAI product (about 10% reduction). Because of delivering superior accuracies and lower uncertainties, the VHGPR models were further applied for LAI mapping using S2 acquisitions over the agricultural sites Marchfeld (Austria) and Barrax (Spain). The models led to consistent LAI maps at BOA and TOA scale. The LAI maps were also compared against LAI maps as generated by the SNAP toolbox, which is based on a neural network (NN). Maps were again consistent, however the SNAP NN model tends to overestimate over dense vegetation cover. Overall, this study demonstrated that hybrid LAI retrieval algorithms can be developed from TOA radiance data given a cloud-free sky, thus without the need of atmospheric correction. To the benefit of the community, the development of such hybrid models for the retrieval vegetation properties from BOA or TOA images has been streamlined in the freely downloadable ALG-ARTMO software framework.

Estévez J., Salinero-Delgado M., Berger K., Pipia L., Rivera-Caicedo J.P., Wocher M., Reyes-Muñoz P., Tagliabue G., Boschetti M., Verrelst J.

The unprecedented availability of optical satellite data in cloud-based computing platforms, such as Google Earth Engine (GEE), opens new possibilities to develop crop trait retrieval models from the local to the planetary scale. Hybrid retrieval models are of interest to run in these platforms as they combine the advantages of physically- based radiative transfer models (RTM) with the flexibility of machine learning regression algorithms. Previous research with GEE primarily relied on processing bottom-of-atmosphere (BOA) reflectance data, which requires atmospheric correction. In the present study, we implemented hybrid models directly into GEE for processing Sentinel-2 (S2) Level-1C (L1C) top-of-atmosphere (TOA) reflectance data into crop traits. To achieve this, a training dataset was generated using the leaf-canopy RTM PROSAIL in combination with the atmospheric model 6SV. Gaussian process regression (GPR) retrieval models were then established for eight essential crop traits namely leaf chlorophyll content, leaf water content, leaf dry matter content, fractional vegetation cover, leaf area index (LAI), and upscaled leaf variables (i.e., canopy chlorophyll content, canopy water content and canopy dry matter content). An important pre-requisite for implementation into GEE is that the models are sufficiently light in order to facilitate efficient and fast processing. Successful reduction of the training dataset by 78% was achieved using the active learning technique Euclidean distance-based diversity (EBD). With the EBD-GPR models, highly accurate validation results of LAI and upscaled leaf variables were obtained against in situ field data from the validation study site Munich-North-Isar (MNI), with normalized root mean square errors (NRMSE) from 6% to 13%. Using an independent validation dataset of similar crop types (Italian Grosseto test site), the retrieval models showed moderate to good performances for canopy-level variables, with NRMSE ranging from 14% to 50%, but failed for the leaf-level estimates. Obtained maps over the MNI site were further compared against Sentinel-2 Level 2 Prototype Processor (SL2P) vegetation estimates generated from the ESA Sentinels' Application Platform (SNAP) Biophysical Processor, proving high consistency of both retrievals (R2 from 0.80 to 0.94). Finally, thanks to the seamless GEE processing capability, the TOA-based mapping was applied over the entirety of Germany at 20 m spatial resolution including information about prediction uncertainty. The obtained maps provided confidence of the developed EBD-GPR retrieval models for integration in the GEE framework and national scale mapping from S2-L1C imagery. In summary, the proposed retrieval workflow demonstrates the possibility of routine processing of S2 TOA data into crop traits maps at any place on Earth as required for operational agricultural applications.

Amin E., Verrelst J., Rivera-Caicedo J.P., Pipia L., Ruiz-Verdú A., Moreno J.

For agricultural applications, identification of non-photosynthetic above-ground vegetation is of great interest as it contributes to assess harvest practices, detecting crop residues or drought events, as well as to better predict the carbon, water and nutrients uptake. While the mapping of green Leaf Area Index (LAI) is well established, current operational retrieval models are not calibrated for LAI estimation over senescent, brown vegetation. This not only leads to an underestimation of LAI when crops are ripening, but is also a missed monitoring opportunity. The high spatial and temporal resolution of Sentinel-2 (S2) satellites constellation offers the possibility to estimate brown LAI (LAI G ) next to green LAI (LAI G ). By using LAI ground measurements from multiple campaigns associated with airborne or satellite spectra, Gaussian processes regression (GPR) models were developed for both LAI G and LAI B , providing alongside associated uncertainty estimates, which allows to mask out unreliable LAI retrievals with higher uncertainties. A processing chain was implemented to apply both models to S2 images, generating a multiband LAI product at 20 m spatial resolution. The models were adequately validated with in-situ data from various European study sites (LAI G : R2 = 0.7, RMSE = 0.67 m2/m2; LAI B : R2 = 0.62, RMSE = 0.43 m2/m2). Thanks to the S2 bands in the red edge (B5: 705 nm and B6: 740 nm) and in the shortwave infrared (B12: 2190 nm) a distinction between LAI G and LAI B can be achieved. To demonstrate the capability of LAI B to identify when crops start senescing, S2 time series were processed over multiple European study sites and seasonal maps were produced, which show the onset of crop senescence after the green vegetation peak. Particularly, the LAI B product permits the detection of harvest (i.e., sudden drop in LAI B ) and the determination of crop residues (i.e., remaining LAI B ), although a better spectral sampling in the shortwave infrared would have been desirable to disentangle brown LAI from soil variability and its perturbing effects. Finally, a single total LAI product was created by merging LAI G and LAI B estimates, and then compared to the LAI derived from S2 L2B biophysical processor integrated in SNAP. The spatiotemporal analysis results confirmed the improvement of the proposed descriptors with respect to the standard SNAP LAI product accounting only for photosynthetically active green vegetation.

Andreu-Perez J., Perez-Espinosa H., Timonet E., Kiani M., Giron-Perez M.I., Benitez-Trinidad A.B., Jarchi D., Rosales-Perez A., Gatzoulis N., Reyes-Galaviz O.F., Torres-Garcia A., Reyes-Garcia C.A., Ali Z., Rivas F.

In an attempt to reduce the infection rate of the COrona VIrus Disease-19 (Covid-19) countries around the world have echoed the exigency for an economical, accessible, point-of-need diagnostic test to identify Covid-19 carriers so that they (individuals who test positive) can be advised to self isolate rather than the entire community. Availability of a quick turn-around time diagnostic test would essentially mean that life, in general, can return to normality-at-large. In this regards, studies concurrent in time with ours have investigated different respiratory sounds, including cough, to recognise potential Covid-19 carriers. However, these studies lack clinical control and rely on Internet users confirming their test results in a web questionnaire (crowdsourcing) thus rendering their analysis inadequate. We seek to evaluate the detection performance of a primary screening tool of Covid-19 solely based on the cough sound from 8,380 clinically validated samples with laboratory molecular-test ( 2,339 Covid-19 positive and 6,041 Covid-19 negative) under quantitative RT-PCR (qRT-PCR) from certified laboratories. All collected samples were clinically labelled, i.e., Covid-19 positive or negative, according to the results in addition to the disease severity based on the qRT-PCR threshold cycle (Ct) and lymphocytes count from the patients. Our proposed generic method is an algorithm based on Empirical Mode Decomposition (EMD) for cough sound detection with subsequent classification based on a tensor of audio sonographs and deep artificial neural network classifier with convolutional layers called ‘DeepCough’ . Two different versions of DeepCough based on the number of tensor dimensions, i.e., DeepCough2D and DeepCough3D, have been investigated. These methods have been deployed in a multi-platform prototype web-app ‘CoughDetect’ . Covid-19 recognition results rates achieved a promising AUC (Area Under Curve) of $98.80\% \pm 0.83\%$ , sensitivity of $96.43\% \pm 1.85\%$ , and specificity of $96.20\% \pm 1.74\%$ and average AUC of $81.08\% \pm 5.05\%$ for the recognition of three severity levels. Our proposed web tool as a point-of-need primary diagnostic test for Covid-19 facilitates the rapid detection of the infection. We believe it has the potential to significantly hamper the Covid-19 pandemic across the world.

De Grave C., Verrelst J., Morcillo-Pallarés P., Pipia L., Rivera-Caicedo J.P., Amin E., Belda S., Moreno J.

The ESA's forthcoming FLuorescence EXplorer (FLEX) mission is dedicated to the global monitoring of the vegetation's chlorophyll fluorescence by means of an imaging spectrometer, FLORIS. In order to properly interpret the fluorescence signal in relation to photosynthetic activity, essential vegetation variables need to be retrieved concomitantly. FLEX will fly in tandem with Sentinel-3 (S3), which conveys the Ocean and Land Colour Instrument (OLCI) that is designed to characterize the atmosphere and the terrestrial vegetation at a spatial resolution of 300 m. In this work we present the retrieval models of four essential biophysical variables: (1) Leaf Area Index (LAI), (2) leaf chlorophyll content (Cab), (3) fraction of absorbed photosynthetically active radiation (fAPAR), and (4) fractional vegetation cover (FCover). These variables can be operationally inferred by hybrid retrieval approaches, which combine the generalization capabilities offered by radiative transfer models (RTMs) with the flexibility and computational efficiency of machine learning methods. The RTM SCOPE (Soil Canopy Observation, Photochemistry and Energy fluxes) was used to generate a database of reflectance spectra corresponding to a large variety of canopy realizations, which served subsequently as input to train a Gaussian Process Regression (GPR) algorithm for each targeted variable. Three sets of GPR models were developed, based on different spectral band settings: (1) OLCI (21 bands between 400 and 1040 nm), (2) FLORIS (281 bands between 500 and 780 nm), and (3) their synergy. Their respective performances were assessed based on simulated reflectance scenes. Regarding the retrieval of Cab, the OLCI model gave good model performances (R2: 0.91; RMSE: 7.6 μg. cm-2), yet superior accuracies were achieved as a result of FLORIS' higher spectral resolution (R2: 0.96; RMSE: 4.8 μg. cm-2). The synergy of both datasets did not further enhance the variable retrieval. Regarding LAI, the improvement of the model performances by using only FLORIS spectra (R2: 0.87; RMSE: 1.05 m2.m-2) rather than only OLCI spectra (R2: 0.86; RMSE: 1.12 m2.m-2) was less evident but merging both data sets was more beneficial (R2: 0.88; RMSE: 1.01 m2.m-2). Finally, the three data sources gave good model performances for the retrieval of fAPAR and Fcover, with the best performing model being the Synergy model (fAPAR: R2: 0.99; RMSE: 0.02 and FCover: R2: 0.98; RMSE: 0.04). The ability of the models to process real data was subsequently demonstrated by applying the OLCI models to S3 surface reflectance products acquired over Western Europe and Argentina. Obtained maps showed consistent patterns and variable ranges, and comparison against corresponding Sentinel-2 products (coarsened to a 300 m spatial resolution) led to reasonable matches (R2: 0.5-0.7). Altogether, given the availability of the multiple data sources, the FLEX tandem mission will foster unique opportunities to quantify essential vegetation properties, and hence facilitate the interpretation of the measured fluorescence levels.