Next Article in Journal
Automatic Extraction of Grasses and Individual Trees in Urban Areas Based on Airborne Hyperspectral and LiDAR Data
Next Article in Special Issue
Crop Type Classification Using Fusion of Sentinel-1 and Sentinel-2 Data: Assessing the Impact of Feature Selection, Optical Data Availability, and Parcel Sizes on the Accuracies
Previous Article in Journal
The Greening and Wetting of the Sahel Have Leveled off since about 1999 in Relation to SST
Previous Article in Special Issue
Prediction of Wheat Grain Protein by Coupling Multisource Remote Sensing Imagery and ECMWF Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

First Insights on Soil Respiration Prediction across the Growth Stages of Rainfed Barley Based on Simulated MODIS and Sentinel-2 Spectral Indices

by
Víctor Cicuéndez
1,
Manuel Rodríguez-Rastrero
2,
Laura Recuero
1,3,
Margarita Huesca
4,
Thomas Schmid
2,
Rosa Inclán
2,
Javier Litago
5,
Víctor Sánchez-Girón
6 and
Alicia Palacios-Orueta
1,3,*
1
Departamento de Sistemas y Recursos Naturales, ETSIMFMN, Universidad Politécnica de Madrid (UPM), 28040 Madrid, Spain
2
Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT), 28040 Madrid, Spain
3
Centro de Estudios e Investigación para la Gestión de Riesgos Agrarios y Medioambientales (CEIGRAM), Universidad Politécnica de Madrid (UPM), 28040 Madrid, Spain
4
Center for Spatial Technologies and Remote Sensing (CSTARS) of the Land, Air and Water Resources Department, University of California Davis, CA 95616, USA
5
Departamento de Economía Agraria, Estadística y Gestión de Empresas, ETSIAAB, Universidad Politécnica de Madrid (UPM), 28040 Madrid, Spain
6
Departamento de Ingeniería Agroforestal, ETSIAAB, Universidad Politécnica de Madrid (UPM), 28040 Madrid, Spain
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(17), 2724; https://doi.org/10.3390/rs12172724
Submission received: 1 July 2020 / Revised: 17 August 2020 / Accepted: 20 August 2020 / Published: 23 August 2020
(This article belongs to the Special Issue Digital Agriculture with Remote Sensing)

Abstract

:
Rainfed agriculture occupies the majority of the world’s agricultural surface and is expected to increase in the near future causing serious effects on carbon cycle dynamics in the context of climate change. Carbon cycle across several temporal and spatial scales could be studied through spectral indices because they are related to vegetation structure and functioning and hence with carbon fluxes, among them soil respiration (Rs). The aim of this work was to assess Rs linked to crop phenology of a rainfed barley crop throughout two seasons based on spectral indices calculated from field spectroscopy data. The relationships between Rs, Leaf Area Index (LAI) and spectral indices were assessed by linear regression models with the adjusted coefficient of determination (Radj2). Results showed that most of the spectral indices provided better information than LAI throughout the studied period and that soil moisture and temperature were relevant variables in specific periods. During vegetative stages, indices based on the visible (VIS) region showed the best relationship with Rs. On the other hand, during reproductive stages indices containing the near infrared-shortwave infrared (NIR-SWIR) spectral region and those related to water content showed the highest relationship. The inter-annual variability found in Mediterranean regions was also observed in the estimated ratio of carbon emission to carbon fixation between years. Our results show the potential capability of spectral information to assess soil respiration linked to crop phenology across several temporal and spatial scales. These results can be used as a basis for the utilization of other remote information derived from satellites or airborne sensors to monitor crop carbon balances.

1. Introduction

Agricultural areas represent 11% of the global land surface [1] generating around 13% of the greenhouse gases emissions and are expected to increase in the near future [2,3]. Soil CO2 efflux, termed also as soil respiration (Rs), accounts for 60–80% of the total ecosystem respiration [4,5] being the main way by which CO2 fixed by land vegetation returns to the atmosphere [6,7]. Soil autotrophic (Ra) and heterotrophic respiration (Rh) are highly dependent on biotic and abiotic factors. Heterotrophic respiration coming from the organic matter mineralization by the decomposer community is highly affected by the availability of soil carbon, moisture and temperature [8,9,10]. The importance of the autotrophic component has been extensively demonstrated [11,12,13]. Ra rate patterns are highly related to photosynthesis [14,15,16] and therefore strongly mediated by phenological dynamics along crop cycles [17,18]. In addition, the accumulation of carbon exudates in roots as well as variation in canopy transpiration and shading through the phenological cycle creates local soil conditions modifying heterotrophic respiration processes [19]. In agricultural fields all these factors result in a high CO2 fluxes spatial and temporal variability [20,21] which may be stronger in rainfed crops in Mediterranean climates due to the irregular distribution of rainfall and temperature [22,23].
Chamber-based methods are commonly used to measure Rs [24], however, measurements are local so that do not account for the Rs spatial and temporal variability at ecosystem scale [25]. Temporal frequency of the measurements has been improved by the use of techniques such automated soil respiration chambers [25] and soil gradient methods [26].However, these methodologies remain providing measurements at plot level making necessary to carry out a large number of them to obtain a reliable representation of spatial variability [27,28,29]. To estimate Rs at larger scales some authors developed models based on Rs from several data bases [30,31]. However, there is a need to develop methods for estimating soil respiration through the growing period on a frequent basis and at large scales [28,32,33].
Remote sensing data is the most appropriate tool for monitoring crop’s phenology across several temporal and spatial scales [34,35]; the link between spectral Indices and physiological processes [36,37] opens new ways to assess soil respiration patterns at large scale.
Huang et al. [38] developed linear regressions models using the Normalized Difference Vegetation Index (NDVI) [39], the Enhanced Vegetation Index (EVI) [40] and the Red Edge Chlorophyll Index (Clred edge) [41] from spectral field measurements in irrigated maize and winter wheat crops. Also in irrigated maize, Huang et al. [42] found the best model to estimate Rs when including Soil Organic Content (SOC) together with satellite EVI data at large scale in an exponential model. In irrigated crops, Huang et al. [33] used EVI from Landsat 8, SOC and bulk density to upscale plot soil respiration, finding better results when using a support vector regression model than a multiple regression model.In the last two studies, however, the models were only appropriate for the central period of the growing season and may not provide reliable estimations in other crop stages.
Based on satellite data only, Huang et al. [43] modeled Rs at large scale in alpine grasslands through the peak growing season using NDVI, EVI and Modified Soil Adjusted Vegetation Index (MSAVI) [44] from MODIS and Landsat data. They found the best estimation when using an NDVI exponential model, due probably to the high proportion of the autotrophic component in total Rs during this period. In cold temperate coniferous sites, Yan et al. [45] developed a model based on MODIS land surface temperature as the primary factor, however, they found that the Rs estimation accuracy was improved when including NDVI.
Cicuéndez et al. [46] included soil moisture and temperature together with field spectroscopy vegetation Indices in multiple linear regression to estimate Rs in an irrigated maize crop. Temperature improved the Rs estimation through all the growing period, while including moisture did not show significant differences. In this work a first assessment of Ra estimation showed significant results during the vegetative stages. Huang and Niu [47] also found a significant influence of abiotic factor to estimate Rs in a rainfed crop with an irregular Rs pattern.
Spectral Shape Indices (SSI) [48] use information of three consecutive bands, taking advantage of spectral shape in different spectral ranges and facilitating the connection between two spectral regions. SSI in the SWIR have been used to monitor crop phenological patterns including periods when photosynthetic activity is not dominant [35,49] as well as in agricultural drought assessment [50].
This concept has also been applied to estimating residue biomass from crops [51,52], to detect spruce bark beetle infections [53] and for non-destructive estimations of foliar carotenoid content of tree species [54] among others. In addition, Cicuéndez et al. [46] have shown that SSI could also be used in Rs assessment in a maize crop.
Most of the studies estimating Rs in agricultural fields have been accomplished in irrigated crops. However, rainfed crops represent 80% of the cultivated land surface and 60% of the agricultural production [55] with a high inter-annual variability, hence there is an urgent need to assess Rs in this type of agriculture at large scale. Specifically, barley (Hordeum vulgare ssp. Vulgare L.) is the fourth most widely grown crop worldwide and the first rainfed crop in extension and production in Spain [56].
The general goal of this work was to assess the capability of several spectral indices, including SSI, based on field hyperspectral information, to assess soil respiration (Rs) linked to crop phenology in a rainfed barley crop in Mediterranean climate. Specific objectives were the following: (1) assessing differences between spectral indices and Leaf Area Index (LAI) to model Rs during vegetative and reproductive stages, (2) to explore the role of soil temperature and moisture on Rs assessment and (3) to carry out a first approximation to assess carbon emission vs. carbon fixation based on the relationship between Rs and barley biomass.

2. Materials and Methods

2.1. Experimental Field

The study was carried out in a rainfed barley field during two crop seasons (2012 and 2013), which is situatedat an altitude approximately of 600 m a.s.l. in the UTM zone 30N, 437410 East and 4476875 North (Datum ETRS89), corresponding to the experimental fields of the Technical School of Agricultural Engineering of the Universidad Politécnica de Madrid (Spain). It is a flat areawhich has been cultivated systematically for 30 years with different crops, including barley. The whole field is characterizedby homogeneous soil, slope, drainage and lighting [46]. The experimental area results from a filling of anthropic origin (“human-transported materials”, [57]) consisting of arkosic materials with an abundance of artifacts. For such reasons, this soil belongs to the Xerorthents Group according to US Soil Taxonomy System [58], andit can be classified as Anthroportic Xerorthent according to Capra et al. [59]. Considering the barley main rooting zone (0–30 cm), soils are well-drained, basic (pH from 7.5 to 8.7), slightly saline (electrical conductivity from 2.5 to 3.1 dS m−1) with a sandy-clay texture and organic carbon content of around 1%.
The area has a Mediterranean continental climate (Csa according to Köppen-Geiger classification) [60]. During the barley growing cycle, average temperatures range from 8.0 °C (March) to 20.1 °C (June) [61]. Rainfall slightly exceeds evapotranspiration during March and April, but since May evapotranspiration is higher than rainfall, defining a remarkable water deficit in rainfed conditions even before the full maturity of barley in June [61].
The field experimental design consisted in three plots with the same soil properties and management practices: two cultivated plots, identified as BY1 and BY2 and a control plot with bare soil (S) (Figure 1). The labels “12” and “13” are referred to the years 2012 and 2013 respectively. The total crop area was 1278 m2 (71 m long × 18 m wide).

2.2. Barley Phenology

In this study, a two-row spring-barley variety “UNIA-R1” was used. The phenological evolution of barley was assessed from March to July during 2012 and 2013. Phenological stages of barley were visually identified and marked according to Zadoks’s cereal growth stages [62]. Based on the results obtained by Cicuéndez et al. [46], we have used the Zadoks’s division in this work with models applied to vegetative and reproductive stages. The Day of Year (DOY) of the stages for 2012 and 2013 is shown in Table 1.
Barley was sown in March in both years, specifically the 10 March 2012 (DOY 70) and the 22 March 2013 (DOY 81). Barley completed its cycle at the end of June or at the beginning of July, and it was completely ripened in July. The cycle from sowing to ripening was shorter in 2012 than in 2013, it lasted 96 days in 2012 and 103 days in 2013. The harvest was carried out in mid-July in both years, a few days after the last spectro-radiometric measurement.

2.3. Cultural Practices

The same experimental layout was used in both growing seasons. In all three plots moldboard plowing was performed as primary tillage (25–30 cm depth). Later, seedbed was prepared by one disk harrowing (10 cm depth), incorporating rabbit manure (50–60 Mg ha−1) to the soil. According to Li-Li et al. [63], rabbit manure constitutes a remarkable source of N, P and K, including 1.8%, 0.6% and 0.7% respectively of such nutrients from dry matter. Then barley UNIA-R1 was sown in BY plots at a rate of 120 kg·ha−1 with a seed drill with fluted coulters and double disk furrow openers. Glyphosate was applied to the bare soil area in some phases during the growing period although the eventual presence of weeds was usually manually controlled during both growing seasons.

2.4. Measurement of Soil CO2Efflux

Measurements of soil CO2 efflux, expressed as Rs (µmols CO2·m−2·s−1), were taken by a portable automated soil CO2 infrared gas analyzer (Li-8100, Li-COR Inc., Lincoln, NE, USA) equipped with a 10 cm survey chamber (Model 8100-102) over PVC collars (105 mm of diameter and 90 mm high) which were placed at a depth of 5–6 cm in the soil and they were cleared of weeds during the growing season.
Five collars were distributed in each plot each year in a semicircular distribution near to spectro-radiometric measurement points (Figure 2). In each collar, an average of three measurements of CO2 efflux was done and then the average of the five collars readings was calculated. As result, one measurement of Rs for each measurement day was estimated in each plot. Collars were located in similar places each year. All Rs measurements were taken between 11 pm and 14 pm.
Total CO2 emissions each year were estimated by integrating the Rs values throughout the growing cycle.

2.5. Field Spectroscopy

An ASD FieldSpec3 Spectroradiometer (Analytical Spectral Devices Inc., Boulder, CO, USA) was used to take hyperspectral ground measurements within a spectral range from the visible (VIS) to Shortwave Infrared (SWIR) region (350–2500 nm), approximately each 15 days, selecting days with clear sky conditions. All measurements were taken between 12 pm and 4 pm.
One reference point was established in each plot for each year. Then, six points were selected at a distance of 1.3 m from the reference point in the cultivated plots according to six angles (90°, 120°, 150°, 180°, 210°, 240°) (Figure 2). In the bare soil plot, three points were selected at a distance of 1.3 m and three points were selected at a distance of 1.8 m. Measurements were made by placing the foreoptic (pistol grip) fixed on a tripod one meter above each of the resulting points.
Three spectro-radiometric measurements were obtained over plant canopy (BY) and bare soil (S) in each point (44cm field of view). Spectroradiometer was calibrated between points using a Spectralon reference panel (ASDInc., Boulder, CO, USA) and measurements were corrected for sensor discrepancy using a splice correction (ViewSpecPro5.6, ASD Inc., Boulder, CO, USA). Then, the average value of the three measurements taken in the six BY and S points was considered. After that, a Savitzky-Golay filter was applied using a second order polynomial with a window size of five spectral bands within the visible and near infrared spectra [46]. Then, a spectral resampling method to MODIS and to Sentinel-2 bands with the function spectral resampling of the software ENVI 4.2 (Research Systems Inc., Boulder, CO, USA) was done. In both cases seven bands in the optical domain were selected. Spectral indices were calculated from these bands (Table 2) and the average of six measurements was computed.

2.6. Measurement of LAI

A LAI-2200 equipment (LI-COR Inc., Lincoln, NE, USA) was used to collect LAI measurements over the growing period.
In both cultivated plots, BY1 and BY2, six points were selected to obtain LAI data, coinciding with the spectro-radiometric measurement points (Figure 2). In diffuse light conditions (dawn, dusk or very cloudy days), three LAI measurements were taken at the base of the plants that cover each of the six points.
These data were later processed through the FV2200 software (1.2.1. version, LI-COR Biosciences), in order to obtain a specific LAI value for each point and date. The average of the six points was calculated to obtain an integrated plot value.

2.7. Measurement of Soil Temperature and Soil Moisture Content

RT-1 and 10HS sensors (Decagon Services Inc., WA, USA) were used to collect soil temperature (°C) and soil moisture (m3 m−3) content at an hourly interval at 10 cm depth in the three plots and registered in Em5b Analog Data Loggers (Decagon Services Inc., WA, USA). From the total data recorded, those hours corresponding to the soil respiration and spectro-radiometric measurements were selected.

2.8. Barley Biomass Determination

Barley biomass produced in each plot in both years was calculated by weighting the harvested plants within a surface equivalent to the field of view of the spectroradiometer (1 m height). This area corresponded to a circle of 22 cm radius which is equivalent to a surface of 0.15 m2. Since a total of six angles have been measured, as indicated in Section 2.5, the surface sampled for each plot was 0.91 m2.Plant cutting was done manually at a height of about 15 cm from the ground, simulating the conditions of a mechanical harvest (shoot biomass). The result was finally expressed in kg ha−1.
Carbon fixation was estimated from the average value of biomass of the two plots of barley (BY1 and BY2). C percentage in barley of the biomass was estimated as Carvajal et al. [65]. These authors established a 43% of C in total shoot biomass and 28% C in total root biomass. Shoot to root biomass ratio was established according to Bolinder et al. [66] and Chirinda et al. [67] that considered a value of two.

2.9. Statistical Analysis

Linear regression models were used to assess the relationships between Rs, LAI and spectral indices. Measurements of both years and of all plots (BY1, BY2 and S plots) were taken into account to build the linear regression models. Additionally, soil temperature and moisture content were included to evaluate their significance in the models. The least square method was used to adjust the models and an analysis of variance (ANOVA) was applied to assess the significance of the model through the F-statistic. In addition, the t-statistic was used to assess the significance of model’s coefficients. The proportion of variance explained by the regression models were measured and compared by the adjusted coefficient of determination (Radj2). The statistical analyses were performed through Statgraphics Centurion XVI (StatPoint Technologies Inc., Warrenton, VA, USA) and SAS software (SAS 9.4 Software, SAS Institute Inc., Cary, NC, USA).

3. Results

For representing the evolution of the variables during the growing period of barley, average values of the three plots were calculated. Two time series for each variable each year were obtained: one representing the average between cultivated plots (BY1 and BY2) and other for the bare soil (S). Table 3 shows the name of the time series of each variable.

3.1. Dynamics of Rs and LAI

LAI and Rs dynamics are shown for both years in Figure 3a,b. LAI values were higher in 2012 than in 2013 with 3.48 m2m−2 and 1.9 m2m−2 maxima values occurring at the beginning of May (DOY 130) and June (DOY 157), respectively. The minimum occurred at the beginning of the growing season (March). LAI dynamics showed similarities between both growing seasons. Values increased sharply from sowing (S) to booting (B) during vegetative stages. After these vegetative stages, LAI values started to decrease during the reproductive stages until the soft dough development stage (D) or a bit later. Then, LAI values remained relatively constant during the rest of the stages of grain filling until maturity and harvest (H).
Although LAI was higher in 2012 than in 2013, Rs in the cultivated plots (BY) had similar maximum values, around 4.5 μmol m−2 s−1 in both years. In 2012, the maximum value occurred at the beginning of April (DOY 100) while in 2013 occurred at mid-May (DOY 133).Rs dynamics in barley showed differences between years. In both years, Rs values increase during the first vegetative stages. Then, in 2012, Rs reached its maximum during tillering (T), one month earlier than maximum LAI (DOY 130) that occurred during stem elongation (SE) and booting (B). After that, values decreased slightly until flowering (F) and remained stable during approximately one month (April-May, DOY 100–139). In 2013, on the other hand, Rs values increased during tillering (T) until stem elongation (SE) when Rs presented a clear maximum (DOY 133), also around one month earlier than maximum LAI (DOY 157). After flowering (F) in 2012 and stem elongation (SE) in 2013, values decreased clearly during the first reproductive stages in 2012, and during booting (B) and also the first reproductive stages in 2013. In both years, minima values occurred at the beginning and the end of the growing season.
Rs in bare soil (i.e., Rh) showed also different dynamics between years. In 2012, Rs increased slightly during March and April. Then, there was a period without measurements for technical problems. After that, values remained constant during June and July. In 2013, Rs increased clearly during April. Then, Rs began to have an irregular pattern increasing and decreasing during April and May. Rs had its maximum value, 2.63 μmol m−2 s−1, in mid-May. After that, Rs decreased during June and July and remained constant until August.

3.2. Dynamics of Spectral Indices

The evolution of spectral indices for MODIS and for Sentinel-2 in 2012 and 2013 is shown in Figure 4 and Figure 5, respectively. The spectral indices followed the evolution of LAI in both years with similar dynamics but different values between them.
In both years, NDVI, EVI, NDWI and NDWI2 increased notably during the vegetative stages from emergence to booting (B) when they reached their maximum values, except EVI in 2013 which reached its maximum value before booting. Then, values started to decrease during reproductive stages until harvest reaching similar minimum values as in the first stages of barley. In NDWI, values at the end of the growing season were higher than values at the beginning of the growing season. NDVI, EVI, NDWI and NDWI2 in the bare soil plot remained regular all the period in both years with values around 0.1 in NDVI, 0.07 in EVI, −0.15 in NDWI and −0.3 in NDWI2.
The SSI ANIR, SASI, AG and AR showed similar dynamics as the previous indices, but they showed an inverse pattern. They presented maximum values at the first stages of barley and then decreased reaching minimum values during booting (B) in 2012 and stem elongation (SE) in 2013. After that, they began to increase during grain filling stages until harvest reaching similar values as at the beginning of the period.

3.3. Dynamics of Soil Temperature and Soil Moisture Content

The dynamics of soil temperature (T10) and soil moisture (H10) at 10 cm depth for both years are shown in Figure 6a,b.
The dynamics of H10was similar in 2012 and 2013 with higher values at early spring and had a clear decrease since mid-May until July. However, there were some differences between the two years. In 2012, soil moisture in cultivated plots had a maximum at the beginning of April (DOY 100) and a slight decrease after this date during the rest of April, although values were high and relatively constant, around 0.300 m3 m−3, during spring until the beginning of May (DOY 130). After this date, values showed a sharp decrease during May and June reaching values around 0.100 m3 m−3. In 2013, soil moisture was high in March until the beginning of April having a maximum of 0.316 m3 m−3 (DOY 99). After that, values began to decrease in mid-April and decreased during the rest of April and May. Although this decline was small, it was stronger than in 2012. Thereafter, values were constant until the end of May (DOY 142) when soil moisture decreased strongly as in 2012 until mid-July (DOY 192) reaching a minimum of 0.051 m3 m−3. In both years, H10 in the bare soil plot (S) had a similar dynamic as the cultivated plots; however, when the soil moisture began to decrease sharply, values were higher than in the latter plots (BY).
T10 in cultivated plots increased during spring and continued increasing during summer in June and July when maximum values were reached for both years. In 2012, soil temperature increased during March and the beginning of April. Values ranged from 11.2 °C in March (DOY 80) to a maximum of 15.2 °C in April (DOY 100). Then, values decreased during the rest of April reaching the minimum value of the growing period of 10 °C (DOY 108). After that, temperatures increased strongly during May and June reaching a maximum of 28.6 °C at the end of June (DOY 158). In 2013, T10 increased steadily during March, April and the beginning of May. Values ranged from the minimum value of the time series of 6.3 °C in March (DOY 78) to a maximum value of 19.2 °C at the beginning of May (DOY 133). After that, T10 decreased strongly at the end of May (DOY 142) reaching a local minimum of 12.4 °C but then temperature began to increase until mid-July reaching a maximum of 27.2 °C (DOY 192). In the bare soil plot (S) the dynamic was similar in both years having slightly higher values than the cultivated plots when temperature was increasing during spring.

3.4. Statistical Relationships between LAI and Spectral Indices

The coefficients of determination (Radj2) of a linear regression model between LAI and spectral indices from MODIS and from Sentinel-2 for the vegetative and the reproductive stages are shown in Table 4 and Table 5, respectively. In MODIS, all linear regression models were significant according to the F-statistic and all the variables were significant in the models according to the t-statistic. The proportion of LAI’s variance explained by the indices was higher during the vegetative stages, with Radj2 values ranging from 0.64 to 0.89, than in the reproductive period, with Radj2 values ranging 0.19 to 0.61. During vegetative stages NDVI, NDWI, EVI, ANIR and SASI explained a higher proportion of LAI’s variance than AG and AR. During reproductive stages, SASI, NDWI and ANIR showed the highest values among all of them.
In Sentinel-2, all linear regression models were significant according to the F-statistic, however, the intercept was not significant in the case of EVI during vegetative stages and AG and AR were not significant during reproductive stages according to the t-statistic. Similar results as those in MODIS were obtained. During vegetative stages NDVI, EVI and NDWI2 explained a higher proportion of LAI’s variance than AG and AR. During reproductive stages, NDWI2 showed the highest value.

3.5. Statistical Relationships between Rs, LAI and Spectral Indices

The coefficients of determination (Radj2) of simple linear regression models between Rs and LAI, and between Rs and spectral indices, and of multiple linear regression models when soil temperature and soil moisture were included as independent variables are shown in Table 6 for MODIS and Table 7 for Sentinel-2 for the vegetative and the reproductive stages.
In MODIS and in Sentinel-2, when only the index was included in the model, coefficients of determination were moderately high, around or above 0.5, with Radj2 similar or higher in the reproductive than in the vegetative period. In addition, all models were significant according to the F-statistic. Variables were significant in all cases throughout all the growing period. Along the vegetative period, the spectral indices explained a higher rate of Rs variance than LAI, while during the reproductive stages, LAI showed a Radj2 higher than the spectral indices, except with ANIR, SASI, NDWI and EVI in the case of MODIS and except NDWI2 in the case of Sentinel-2.
In both cases, when T10was included in the models, coefficients of determination increased in the vegetative stages for all indices while in the reproductive stages they increased only for LAI. In general, during the vegetative stages, variables were significant while during the reproductive stages T10 was not significant in any of the spectral indices and the F-statistic values decreased when T10 was included.
In both cases, H10 was not significant according to the t-statistic during the vegetative stages, in addition, the Radj2 and F-values decreased in all cases when the models took into account soil moisture. In the reproductive stages, Radj2 values were higher than in the vegetative stages, but Radj2 and F-statistic values were in general lower than when H10 was not included. In addition, soil moisture was not significant in all spectral indices. However, it has to be remarked that LAI plus H10 explained the highest proportion of Rs variance (0.85).
In MODIS and in Sentinel-2, when the models included together both abiotic factors, Radj2 values increased in the vegetative stages and F-values remained significant in all cases. It has to be remarked that H10 was not significant when it was included alone with the index, however, when it was included together with T10 it was significant for all indices except LAI andin MODIS NDWI. In the reproductive stages, Radj2 values were lower than when only the index was included in the model and in general lower than in the vegetative stages, F-values also decreased. In addition, neither T10 nor H10 were significant in all indices except in the LAI in which soil moisture was significant.
Table 8 shows the coefficients of determination (Radj2) of the simple linear regression models between Rs and soil temperature (T10), between Rs and soil moisture (H10), and of multiple linear regression models when both variables were considered as independent variables, for both stages in the cultivated plots. On one hand, T10 showed low Radj2 values during the vegetative stages and it did not explain Rs variance during the reproductive stages. In both cases, T10 was significant according to the t-statistic. On the other hand, H10 showed moderate values in the reproductive stages, and null values in the vegetative stages. When both variables were included together in the models Radj2 values did not increase.
Table 8 also shows the results in the bare soil plots, including Radj2 values of the simple linear regression models between Rs, i.e., heterotrophic respiration (Rh), and soil temperature (T10), between Rs and soil moisture (H10), and of multiple linear regression models when T10 and H10 were considered along vegetative and reproductive stages. If the growing period was divided according to the vegetative and reproductive stages in the crop the Radj2 increased. In the vegetative stages, T10 was significant and presented a moderate Radj2 and H10 was not significant either alone or together with T10. Meanwhile, in the reproductive stages soil moisture was significant and presented a moderate Radj2, and although T10 showed a low Radj2 alone and increased the Radj2 together with H10, it was not significant in the models.

3.6. Relationship between Soil CO2Efflux and Biomass

Total CO2 emissions for the growing cycle were 1505.16 g CO2 m−2 in 2012 and 1009.28 g CO2 m−2 in 2013, corresponding to410.91g of C m−2 and 275.53 g of C m−2 respectively. The aboveground biomass at harvest was 8310.1 kg ha−1 in 2012 and 4700 kg ha−1 in 2013. Thus, the final value of C fixation was 473.67 g C m−2 in 2012 and 267.9 g C m−2 in 2013 resulting in estimated Cemitted/Cfixed ratios of 0.87 in 2012 and 1.03 in 2013.

4. Discussion

During vegetative stages, Rs increased as a consequence of barley growth including the development of the root system (Figure 3). During reproductive stages, Rs decreased due to senescence of barley, which coincides with the decrease in soil moisture availability (Figure 6). This is in accordance with the expected evolution of both components of Rs: (1) Ra increases due to crop growth and then decreases due to senescence and (2) Rh increases together with soil temperature in spring and decreases together with soil moisture in summer. Jongen et al. [68] found that ecosystem respiration and Gross Primary Production decreased with the canopy senescence in Mediterranean grasslands. Similar patterns were found in irrigated crops [38,46]. Other authors attribute this behavior to a decrease in photosynthesis rate [69] and a reduction of C translocation to the root system during the grain filling processes [70,71].
Significant differences in LAI, Rs, and indices between years were found which is consistent with the high inter-annual variability in crop progression, typical of Mediterranean regions [72,73]. Spring barley development during vegetative stages required optimal temperature and water availability [74]. The adequate conditions of temperature and moisture at the beginning of the growing season in 2012, permitted a fast increase in Rs values, due to a higher activity of metabolism and a faster development [74,75]. On the other hand, anomalous low temperatures in 2013 March-April resulted in a delay of crop development at the beginning of the growing period (Figure 3b), which resulted in lower Rs, indices, LAI and biomass values throughout the crop cycle. Optimal moisture and temperature conditions of May 2013 gave rise to high Rs values later in the growing season. This is probably more associated to heterotrophic processes, so that soil heterotrophic communities could take advantage of such environmental conditions to a greater extent than the crop. These patterns showed the key role of the moisture-temperature timing in crop development [76,77] and, specifically, its effect on soil respiration [8,78].
During vegetative stages significant crop structural changes resulted in a high proportion of LAI variance (Table 4 and Table 5) explained by NDVI, EVI, NDWI, NDWI2, ANIR and SASI i.e., those indices with influence of the NIR band. However, AG and AR, centered in the VIS spectral region, which is more related to functional traits [79], showed more moderate Radj2 values. During reproductive stages, with smaller LAI variation, only NDWI, NDWI2, ANIR and SASI (i.e., those with higher influence of the NIR-SWIR spectral region) maintained moderate values being higher in NDWI and SASI from MODIS and NDWI2 from Sentinel-2, which are the only indices that do not include the red band. This is in accordance with other researches which showed the relevance of the NIR-SWIR spectral region for assessing structural changes during crop development [80,81].
The proportion of Rs variance explained by LAI indicated that there must be other sources of variability that have a significant influence on soil CO2 efflux. This fact explained the increase in Radj2 values when including soil temperature and moisture in the models of the vegetative and reproductive stages respectively, confirming the limiting role of temperature at the beginning of the growing period and soil moisture at the end [74]. In addition, high LAI Radj2 values during reproductive stages suggest that significant structural changes due to flowering and senescence are linked to Rs.
In both stages, some spectral indices showed significantly higher Radj2 values than LAI alone indicating that they contained additional information (probably functional). During the vegetative stages, NDWI, ANIR and SASI from MODIS and NDWI2 from Sentinel-2, showed coefficients of determination similar to those of the LAI, corroborating the strong relationship between LAI and these indices [52,82]. Indices centered in the VIS spectral region (i.e., AG, ARfrom both sensors), highly related to functional traits such as pigment content and photosynthesis [54,83], showed slightly higher values than LAI and the rest of the indices. This agrees with the evolution of these indices showing similar dynamics as Rs with distinct differences between years (Figure 4 and Figure 5). By including the temperature in the AG and AR models the Radj2 improvement was smaller than in the case of LAI (Table 6 and Table 7, second row). This fact suggests that these indices were already accounting for temperature effects in vegetation photosynthetic activity and emphasizes the importance of this variable on Rs variability at the beginning of spring. The results shown in the bare soil plots (Table 8) corroborated that temperature had a strong influence also in heterotrophic respiration. On the other hand, moisture did not seem to play an important role in Rs during the vegetative stages (Table 6 and Table 7, third row), indicating that moisture is not a limiting factor either for the crop or for the soil microorganisms, as shown in Table 8.
The senescence process and drought stress resulted in alterations in leaf cell structure and composition [84,85] and a decrease in leaf moisture that has been well captured by the NDWI [64,86] and SASI. In our experiment we found high NDWI, SASI and NDWI2 Radj2 values with both LAI (Table 4 and Table 5) and Rs, indicating that the capability of these indices to capture variability in vegetation structure and moisture could be the basis in using them for assessing Rs. The importance of moisture when only heterotrophic respiration is considered (Table 8, Radj2 = 0.52) decreases when the crop is taken into account (Radj2 = 0.38) evidencing the importance of crop senescence in soil CO2 efflux dynamics. During this period ANIR showed also a high Radj2 value indicating that the structural information based on ANIR variability could help Rs assessment.
The difference between years in carbon emitted versus carbon fixed could be associated to abiotic factors dynamics; better conditions at the beginning of the growing cycle in 2012 resulted in a clear advantage for the plant to produce higher biomass. The inter-annual variability in the carbon flux balances has been shown also by other studies especially in Mediterranean climate regions [87,88,89]. In this aspect, taking into account photorespiration and root biomass values would be a significant contribution to more accurate estimations.
Spectral indices from MODIS and Sentinel-2 bands have shown to partially represent change in soil respiration, a highly dynamic process through the growing period, based on their link with crop physiology. While they did not provide direct information on soil microbial processes, including soil moisture and temperature, which are more easily estimated than soil respiration, improved the results.
It is recognized that measurements taken in this study, mostly throughout time, were subjected to auto-correlation at some degree. Effectively, and for the experimental layout adopted, results do not reflect the ability of spectral indices to assess Rs directly, but its change through time. Due to limitations to accurately assess soil respiration at large scales and long time periods, this study should be regarded as preliminary and providing useful insights on the combination of remote sensing data with soil temperature and moisture to potentially assess soil respiration variability across crops and crop stages. Future studies should gather more information across a range of environments (e.g., different soil types) and through longer time periods to develop robust and representative statistical models for Rs assessment.

5. Conclusions

This study showed that spectral indices improve the information provided by LAI to assess soil respiration in a rainfed barley crop. Soil moisture and temperature were relevant variables to complement spectral indices and/or LAI to estimate Rs in specific periods. During vegetative stages, indices together with soil temperature showed similar results than LAI with soil temperature, while including soil moisture with LAI improved Rs estimations during senescence.
Indices including the near infrared spectral region (ANIR, NDVI, EVI, NDWI, NDWI2 and SASI) were the most related to LAI during vegetative stages. Among them, moisture indices, i.e., SASI and NDWI from MODIS and NDWI2 from Sentinel-2, showed the best relationship with LAI during reproductive stages, indicating the impact of declining moisture during senescence on plant structure.
AG and AR indices, centered in the visible spectral region, showed the best relationship with Rs during vegetative stages characterized by active photosynthesis. On the other hand, the NIR-SWIR spectral indices (especially SASI from MODIS and NDWI2 from Sentinel-2) showed the highest relationships with Rs during reproductive stages showing the impact of plant senescence. Furthermore, SASI (angle centered at SWIR1) showed better results than NDWI, probably because it provides information farther in the shortwave spectral region highly related to water content.
This study showed the potential capability of spectral information to assess soil respiration linked to crop physiology. Our results based on measurements acquired through field spectroscopy are a first step to use these models across several spatial and temporal scales to monitor crop C balances under high inter-annual climate variability that result in distinct Cemitted/Cfixed ratios between years. Our study emphasizes the usefulness of a frequent monitoring system in Mediterranean agricultural environments where equilibrium between fixation and emission is particularly fragile, especially in the context of climate change.

Author Contributions

Conceptualization, V.C., M.R.-R., M.H., T.S., R.I., A.P.-O.; methodology, V.C., M.R.-R., L.R., M.H., T.S., R.I., V.S.-G.; software, T.S., R.I., J.L., A.P.-O.; formal analysis, V.C., M.R.-R., L.R., M.H., J.L.; investigation, V.C., M.R.-R., L.R., A.P.-O.; resources, T.S., R.I., A.P.-O.; writing—original draft preparation, V.C., M.R.-R., L.R.; writing—review and editing, V.C., M.R.-R., L.R., J.L., V.S.-G., A.P.-O.; supervision, A.P.-O.; project administration, T.S., R.I., V.S.-G., A.P.-O.; funding acquisition, T.S., R.I., A.P.-O. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministerio de Ciencia e Innovación (Spain) through the projects AGL-2010-17505 and CGL-2009-07031. This work was also funded supported by the Ministerio de Ciencia, Innovación y Universidades (Spain) through the pre-doctoral scholarship of Laura Recuero (Becas de Formación Profesorado Universitario, BOE-A-2015-9456).

Acknowledgments

We would like to thank the College of Agricultural Engineering of the Universidad Politécnica de Madrid for the concession of their experimental fields to carry out this research and to thank all the staff, especially Román Zurita, for their work in planting and taking care of the crop during both years. We would also like to acknowledge David Manrique for his collaboration in the field work.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Bruinsma, J. World Agriculture: Towards 2015/2030; Bruinsma, J., Ed.; Routledge: London, UK, 2003; ISBN 9781315083858. [Google Scholar]
  2. Lal, R. Carbon emission from farm operations. Environ. Int. 2004, 30, 981–990. [Google Scholar] [CrossRef] [PubMed]
  3. IPCC. Climate Change 2014: Synthesis Report; Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Core Writing Team, Pachauri, R.K., Meyer, L.A., Eds.; IPCC: Geneva, Switzerland, 2014. [Google Scholar]
  4. Goulden, M.L.; Munger, J.W.; Fan, S.M.; Daube, B.C.; Wofsy, S.C. Exchange of Carbon Dioxide by a Deciduous Forest: Response to Interannual Climate Variability. Science 1996, 271, 1576–1578. [Google Scholar] [CrossRef] [Green Version]
  5. Suleau, M.; Moureaux, C.; Dufranne, D.; Buysse, P.; Bodson, B.; Destain, J.P.; Heinesch, B.; Debacq, A.; Aubinet, M. Respiration of three Belgian crops: Partitioning of total ecosystem respiration in its heterotrophic, above- and below-ground autotrophic components. Agric. For. Meteorol. 2011, 151, 633–643. [Google Scholar] [CrossRef]
  6. Raich, J.W.; Potter, C.S.; Bhagawati, D. Interannual variability in global soil respiration, 1980–1994. Glob. Chang. Biol. 2002, 8, 800–812. [Google Scholar] [CrossRef]
  7. Schlesinger, W.H.; Andrews, J.A. Soil respiration and the global carbon cycle. Biogeochemistry 2000, 48, 7–20. [Google Scholar] [CrossRef]
  8. Curiel Yuste, J.; Baldocchi, D.D.; Gershenson, A.; Goldstein, A.; Misson, L.; Wong, S. Microbial soil respiration and its dependency on carbon inputs, soil temperature and moisture. Glob. Chang. Biol. 2007, 13, 2018–2035. [Google Scholar] [CrossRef] [Green Version]
  9. Moyano, F.E.; Manzoni, S.; Chenu, C. Responses of soil heterotrophic respiration to moisture availability: An exploration of processes and models. Soil Biol. Biochem. 2013, 59, 72–85. [Google Scholar] [CrossRef]
  10. Wang, W.J.; Dalal, R.C.; Moody, P.W.; Smith, C.J. Relationships of soil respiration to microbial biomass, substrate availability and clay content. Soil Biol. Biochem. 2003, 35, 273–284. [Google Scholar] [CrossRef]
  11. Baggs, E.M. Partitioning the components of soil respiration: A research challenge. Plant Soil 2006, 284, 1–5. [Google Scholar] [CrossRef]
  12. Bond-Lamberty, B.; Wang, C.; Gower, S.T. A global relationship between the heterotrophic and autotrophic components of soil respiration? Glob. Chang. Biol. 2004, 10, 1756–1766. [Google Scholar] [CrossRef]
  13. Ferréa, C.; Zenone, T.; Comolli, R.; Seufert, G. Estimating heterotrophic and autotrophic soil respiration in a semi-natural forest of Lombardy, Italy. Pedobiologia 2012, 55, 285–294. [Google Scholar] [CrossRef]
  14. Kuzyakov, Y.; Cheng, W. Photosynthesis controls of CO2 efflux from maize rhizosphere. Plant Soil 2004, 263, 85–99. [Google Scholar] [CrossRef]
  15. Kuzyakov, Y.; Cheng, W. Photosynthesis controls of rhizosphere respiration and organic matter decomposition. Soil Biol. Biochem. 2001, 33, 1915–1925. [Google Scholar] [CrossRef]
  16. Xu, X.; Kuzyakov, Y.; Wanek, W.; Richter, A. Root-derived respiration and non-structural carbon of rice seedlings. Eur. J. Soil Biol. 2008, 44, 22–29. [Google Scholar] [CrossRef]
  17. Fu, S.; Cheng, W.; Susfalk, R. Rhizosphere respiration varies with plant species and phenology: A greenhouse pot experiment. Plant Soil 2002. [Google Scholar] [CrossRef]
  18. Jans, W.W.P.; Jacobs, C.M.J.; Kruijt, B.; Elbers, J.A.; Barendse, S.; Moors, E.J. Carbon exchange of a maize (Zea mays L.) crop: Influence of phenology. Agric. Ecosyst. Environ. 2010, 139, 316–324. [Google Scholar] [CrossRef]
  19. Von Haden, A.C.; Marín-Spiotta, E.; Jackson, R.D.; Kucharik, C.J. Soil microclimates influence annual carbon loss via heterotrophic soil respiration in maize and switchgrass bioenergy cropping systems. Agric. For. Meteorol. 2019, 279, 107731. [Google Scholar] [CrossRef]
  20. Rochette, P.; Desjardins, R.L.; Pattey, E. Spatial and temporal variability of soil respiration in agricultural fields. Can. J. Soil Sci. 1991, 71, 189–196. [Google Scholar] [CrossRef]
  21. Han, G.; Zhou, G.; Xu, Z.; Yang, Y.; Liu, J.; Shi, K. Biotic and abiotic factors controlling the spatial and temporal variation of soil respiration in an agricultural ecosystem. Soil Biol. Biochem. 2007, 39, 418–425. [Google Scholar] [CrossRef]
  22. González-Ubierna, S.; Teresa De La Cruz, M.; Casermeiro, M.Á. Climate factors mediate soil respiration dynamics in Mediterranean agricultural environments: An empirical approach. Soil Res. 2014, 52, 543–553. [Google Scholar] [CrossRef]
  23. Ma, S.; Baldocchi, D.D.; Xu, L.; Hehn, T. Inter-annual variability in carbon dioxide exchange of an oak/grass savanna and open grassland in California. Agric. For. Meteorol. 2007, 147, 157–171. [Google Scholar] [CrossRef]
  24. Ryan, M.G.; Law, B.E. Interpreting, measuring, and modeling soil respiration. Biogeochemistry 2005, 73, 3–27. [Google Scholar] [CrossRef]
  25. Vargas, R.; Baldocchi, D.D.; Allen, M.F.; Bahn, M.; Andrew Black, T.; Collins, S.L.; Curiel Yuste, J.; Hirano, T.; Jassal, R.S.; Pumpanen, J.; et al. Looking deeper into the soil: Biophysical controls and seasonal lags of soil CO₂ production and efflux. Ecol. Appl. 2010, 20, 1569–1582. [Google Scholar] [CrossRef] [PubMed]
  26. Hirano, T.; Kim, H.; Tanaka, Y. Long-term half-hourly measurement of soil CO2 concentration and soil respiration in a temperate deciduous forest. J. Geophys. Res. Atmos. 2003, 108. [Google Scholar] [CrossRef]
  27. Barba, J.; Curiel Yuste, J.; Martínez-Vilalta, J.; Lloret, F. Drought-induced tree species replacement is reflected in the spatial variability of soil respiration in a mixed Mediterranean forest. For. Ecol. Manag. 2013, 306, 79–87. [Google Scholar] [CrossRef]
  28. Barba, J.; Cueva, A.; Bahn, M.; Barron-Gafford, G.A.; Bond-Lamberty, B.; Hanson, P.J.; Jaimes, A.; Kulmala, L.; Pumpanen, J.; Scott, R.L.; et al. Comparing ecosystem and soil respiration: Review and key challenges of tower-based and soil measurements. Agric. For. Meteorol. 2018, 249, 434–443. [Google Scholar] [CrossRef] [Green Version]
  29. Pantani, O.L.; Fioravanti, F.; Stefanini, F.M.; Berni, R.; Certini, G. Accounting for soil respiration variability – Case study in a Mediterranean pine-dominated forest. Sci. Rep. 2020, 10, 1–10. [Google Scholar] [CrossRef] [Green Version]
  30. Zhao, Z.; Peng, C.; Yang, Q.; Meng, F.R.; Song, X.; Chen, S.; Epule, T.E.; Li, P.; Zhu, Q. Model prediction of biome-specific global soil respiration from 1960 to 2012. Earth’s Future 2017, 5, 715–729. [Google Scholar] [CrossRef]
  31. Li, J.; Pei, J.; Pendall, E.; Fang, C.; Nie, M. Spatial heterogeneity of temperature sensitivity of soil respiration: A global analysis of field observations. Soil Biol. Biochem. 2020, 141, 107675. [Google Scholar] [CrossRef]
  32. Bond-Lamberty, B.; Thomson, A. A global database of soil respiration data. Biogeosciences 2010, 7, 1915–1926. [Google Scholar] [CrossRef] [Green Version]
  33. Huang, N.; Wang, L.; Guo, Y.; Niu, Z. Upscaling plot-scale soil respiration in winter wheat and summer maize rotation croplands in Julu County, North China. Int. J. Appl. Earth Obs. Geoinf. 2017, 54, 169–178. [Google Scholar] [CrossRef]
  34. Massey, R.; Sankey, T.T.; Congalton, R.G.; Yadav, K.; Thenkabail, P.S.; Ozdogan, M.; Sánchez Meador, A.J. MODIS phenology-derived, multi-year distribution of conterminous U.S. crop types. Remote Sens. Environ. 2017, 198, 490–503. [Google Scholar] [CrossRef]
  35. Tornos, L.; Huesca, M.; Dominguez, J.A.; Moyano, M.C.; Cicuendez, V.; Recuero, L.; Palacios-Orueta, A. Assessment of MODIS spectral indices for determining rice paddy agricultural practices and hydroperiod. ISPRS J. Photogramm. Remote Sens. 2015, 101, 110–124. [Google Scholar] [CrossRef] [Green Version]
  36. Glenn, E.P.; Huete, A.R.; Nagler, P.L.; Nelson, S.G. Relationship between remotely-sensed vegetation indices, canopy attributes and plant physiological processes: What vegetation indices can and cannot tell us about the landscape. Sensors 2008, 8, 2136–2160. [Google Scholar] [CrossRef] [Green Version]
  37. Xue, J.; Su, B. Significant remote sensing vegetation indices: A review of developments and applications. J. Sens. 2017, 2017. [Google Scholar] [CrossRef] [Green Version]
  38. Huang, N.; Niu, Z.; Zhan, Y.; Xu, S.; Tappert, M.C.; Wu, C.; Huang, W.; Gao, S.; Hou, X.; Cai, D. Relationships between soil respiration and photosynthesis-related spectral vegetation indices in two cropland ecosystems. Agric. For. Meteorol. 2012, 160, 80–89. [Google Scholar] [CrossRef]
  39. Tucker, C.J. Red and photographic infrared linear combinations for monitoring vegetation. Remote Sens. Environ. 1979, 8, 127–150. [Google Scholar] [CrossRef] [Green Version]
  40. Huete, A.; Didan, K.; Miura, T.; Rodriguez, E.; Gao, X.; Ferreira, L. Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sens. Environ. 2002, 83, 195–213. [Google Scholar] [CrossRef]
  41. Gitelson, A.A.; Viña, A.; Ciganda, V.; Rundquist, D.C.; Arkebauer, T.J. Remote estimation of canopy chlorophyll content in crops. Geophys. Res. Lett. 2005, 32. [Google Scholar] [CrossRef] [Green Version]
  42. Huang, N.; Wang, L.; Guo, Y.; Hao, P.; Niu, Z. Modeling spatial patterns of soil respiration in maize fields from vegetation and soil property factors with the use of remote sensing and geographical information system. PLoS ONE 2014, 9, e105150. [Google Scholar] [CrossRef] [Green Version]
  43. Huang, N.; He, J.S.; Niu, Z. Estimating the spatial pattern of soil respiration in Tibetan alpine grasslands using Landsat TM images and MODIS data. Ecol. Indic. 2013, 26, 117–125. [Google Scholar] [CrossRef]
  44. Qi, J.; Chehbouni, A.; Huete, A.R.; Kerr, Y.H.; Sorooshian, S. A modify soil adjust vegetation index. Remote Sens. Environ. 1994, 126, 119–126. [Google Scholar] [CrossRef]
  45. Yan, J.; Zhang, X.; Liu, J.; Li, H.; Ding, G. MODIS-Derived estimation of soil respiration within five cold temperate coniferous forest sites in the eastern Loess Plateau, China. Forests 2020, 11, 131. [Google Scholar] [CrossRef] [Green Version]
  46. Cicuéndez, V.; Rodríguez-Rastrero, M.; Huesca, M.; Uribe, C.; Schmid, T.; Inclán, R.; Litago, J.; Sánchez-Girón, V.; Merino-de-Miguel, S.; Palacios-Orueta, A. Assessment of soil respiration patterns in an irrigated corn field based on spectral information acquired by field spectroscopy. Agric. Ecosyst. Environ. 2015, 212, 158–167. [Google Scholar] [CrossRef]
  47. Huang, N.; Niu, Z. Estimating soil respiration using spectral vegetation indices and abiotic factors in irrigated and rainfed agroecosystems. Plant Soil 2013, 367, 535–550. [Google Scholar] [CrossRef]
  48. Palacios-Orueta, A.; Khanna, S.; Litago, J.; Whiting, M.L.; Ustin, S.L. Assessment of NDVI and NDWI spectral indices using MODIS time series analysis and development of a new spectral index based on MODIS shortwave infrared bands. In Proceedings of the 1st International Conference of Remote Sensing and Geoinformation Processing, Trier, Germany, 15–17 November 2006; pp. 207–209. [Google Scholar]
  49. Palacios-Orueta, A.; Huesca, M.; Whiting, M.L.; Litago, J.; Khanna, S.; Garcia, M.; Ustin, S.L. Derivation of phenological metrics by function fitting to time-series of Spectral Shape Indexes AS1 and AS2: Mapping cotton phenological stages using MODIS time series. Remote Sens. Environ. 2012, 126, 148–159. [Google Scholar] [CrossRef]
  50. Das, P.K.; Murthy, C.S.; Mvr, S. Monitoring of seasonal dryness/wetness conditions using shortwave angle slope index for early season agricultural drought assessment. Geomat. Nat. Hazards Risk 2014, 5, 232–251. [Google Scholar] [CrossRef]
  51. Zhang, M.; Wu, B.; Meng, J. Quantifying winter wheat residue biomass with a spectral angle index derived from China Environmental Satellite data. Int. J. Appl. Earth Obs. Geoinf. 2014, 32, 105–113. [Google Scholar] [CrossRef]
  52. Khanna, S.; Palacios-Orueta, A.; Whiting, M.L.; Ustin, S.L.; Riaño, D.; Litago, J. Development of angle indexes for soil moisture estimation, dry matter detection and land-cover discrimination. Remote Sens. Environ. 2007, 109, 154–165. [Google Scholar] [CrossRef]
  53. Fassnacht, F.E.; Latifi, H.; Koch, B. An angular vegetation index for imaging spectroscopy data-Preliminary results: On forest damage detection in the Bavarian National Park, Germany. Int. J. Appl. Earth Obs. Geoinf. 2012, 19, 308–321. [Google Scholar] [CrossRef]
  54. Fassnacht, F.E.; Stenzel, S.; Gitelson, A.A. Non-destructive estimation of foliar carotenoid content of tree species using merged vegetation indices. J. Plant Physiol. 2015, 176, 210–217. [Google Scholar] [CrossRef] [PubMed]
  55. Wani, S.P.; Rockström, J.; Oweis, T. (Eds.) Rainfed Agriculture: Unlocking the Potential; CAB International: Wallingford, UK, 2009; ISBN 978-1-84593-389-0. [Google Scholar]
  56. MAPA Encuesta Sobre Superficies y Rendimientos de Cultivos (ESYRCE). Available online: https://www.mapa.gob.es/es/estadistica/temas/estadisticas-agrarias/agricultura/esyrce/ (accessed on 3 January 2019).
  57. Capra, G.F.; Ganga, A.; Grilli, E.; Vacca, S.; Buondonno, A. A review on anthropogenic soils from a worldwide perspective. J. Soils Sediments 2015, 15, 1602–1618. [Google Scholar] [CrossRef]
  58. Soil Survey Staff. Keys to Soil Taxonomy, 12th ed.; USDA-NRCS, Ed.; United States Department of Agriculture (USDA)-Natural Resources Conservation Service (NRCS): Washington, DC, USA, 2014.
  59. Capra, G.F.; Vacca, S.; Cabula, E.; Grilli, E.; Buondonno, A. Through the Decades: Taxonomic Proposals for Human-Altered and Human-Transported Soil Classification. Soil Horiz. 2013, 54. [Google Scholar] [CrossRef]
  60. AEMET/IMP. Atlas Climático Ibérico/Iberian Climate Atlas; AEMET—Agencia Estatal de Meteorología, IMP—Instituto de Meteorología de Portugal, Ministerio de Medio Ambiente Rural y Marino, Ed.; Closas Orcoyen S.L.: Madrid, Spain, 2011; ISBN 978-84-7837-079-5. [Google Scholar]
  61. MAPA. Sistema de Información Geográfica de Datos Agrarios. Available online: https://sig.mapama.gob.es/siga/ (accessed on 18 February 2019).
  62. Zadoks, J.C.; Chang, T.T.; Konzak, C.F. A decimal code for the growth stages of cereals. Weed Res. 1974, 14, 415–421. [Google Scholar] [CrossRef]
  63. Li-Li, B.; Tie-Jun, Y.; Bin, W.; Lin, B.; De-Gui, T.; Xiang-Chao, F. Evaluation and Comparison of Composting Rabbit Manure mixed with Mushroom Residue and Rice Straw. J. Agric. Sci. Technol. 2013, 15, 1069–1081. [Google Scholar]
  64. Gao, B. NDWI—A normalized difference water index for remote sensing of vegetation liquid water from space. Remote Sens. Environ. 1996, 58, 257–266. [Google Scholar] [CrossRef]
  65. Carvajal, M.; Mota, C.; Alcaraz-López, C.; Iglesias, M. Investigación Sobre la Absorción de CO2 por los Cultivos más Representativos de la Región de Murcia. Available online: http://cambioclimaticomurcia.carm.es/pdfs/lessco2/ponencia_cisc_espanol.pdf (accessed on 30 April 2018).
  66. Bolinder, M.A.; Angers, D.A.; Dubuc, J.P. Estimating shoot to root ratios and annual carbon inputs in soils for cereal crops. Agric. Ecosyst. Environ. 1997, 63, 61–66. [Google Scholar] [CrossRef]
  67. Chirinda, N.; Olesen, J.E.; Porter, J.R. High Root Biomass for Cereal Crops Increases Carbon Sequestration in Organic Arable Systems. In Proceedings of the 17th IFOAM Organic World Congress, October 2011; pp. 2–5. Available online: http://orgprints.org/19499 (accessed on 22 August 2020).
  68. Jongen, M.; Pereira, J.S.; Aires, L.M.I.; Pio, C.A. The effects of drought and timing of precipitation on the inter-annual variation in ecosystem-atmosphere exchange in a Mediterranean grassland. Agric. For. Meteorol. 2011, 151, 595–606. [Google Scholar] [CrossRef]
  69. Viña, A.; Gitelson, A.A.; Rundquist, D.C.; Keydan, G.; Leavitt, B.; Schepers, J. Monitoring maize (L.) Phenology with Remote Sensing. Agron. J. 2004, 96, 1139. [Google Scholar] [CrossRef]
  70. Amos, B.; Arkebauer, T.J.; Doran, J.W. Soil Surface Fluxes of Greenhouse Gases in an Irrigated Maize-Based Agroecosystem. Soil Sci. Soc. Am. J. 2005, 69, 387–395. [Google Scholar] [CrossRef] [Green Version]
  71. Qian, J.H.; Doran, J.W.; Walters, D.T. Maize plant contributions to root zone available carbon and microbial transformations of nitrogen. Soil Biol. Biochem. 1997, 29, 1451–1462. [Google Scholar] [CrossRef]
  72. Ceglar, A.; Toreti, A.; Lecerf, R.; Van der Velde, M.; Dentener, F. Impact of meteorological drivers on regional inter-annual crop yield variability in France. Agric. For. Meteorol. 2016, 216, 58–67. [Google Scholar] [CrossRef]
  73. Vicente-Serrano, S.M.; Cuadrat-Prats, J.M.; Romo, A. Early prediction of crop production using drought indices at different time-scales and remote sensing data: Application in the Ebro Valley (north-east Spain). Int. J. Remote Sens. 2006, 27, 511–518. [Google Scholar] [CrossRef]
  74. Hossain, A.; Teixeira da Silva, J.A.; Lozovskaya, M.V.; Zvolinsky, V.P. High temperature combined with drought affect rainfed spring wheat and barley in South-Eastern Russia: I. Phenology and growth. Saudi J. Biol. Sci. 2012, 19, 473–487. [Google Scholar] [CrossRef] [Green Version]
  75. Rawson, H.M. Effect of high temperatures on the development and yield of wheat and practices to reduce deleterious effects. In Proceedings of the Wheat Production Constraints in Tropical Environments, Chiang Mai, Thailand, 19–23 January 1987. [Google Scholar]
  76. Oteros, J.; García-Mozo, H.; Botey, R.; Mestre, A.; Galán, C. Variations in cereal crop phenology in Spain over the last twenty-six years (1986–2012). Clim. Chang. 2015, 130, 545–558. [Google Scholar] [CrossRef]
  77. Wielgolaski, F.-E. Phenology in Agriculture. Phenol. Seas. Model. 1974, 369–381. [Google Scholar] [CrossRef]
  78. Raich, J.; Tufekcioglu, A. Vegetation and Soil Respiration: Correlations and Controls. Biogeochemistry 2000, 48, 71–90. [Google Scholar] [CrossRef]
  79. Sims, D.A.; Gamon, J.A. Relationships between leaf pigment content and spectral reflectance across a wide range of species, leaf structures and developmental stages. Remote Sens. Environ. 2002, 81, 337–354. [Google Scholar] [CrossRef]
  80. Ali, A.M.; Darvishzadeh, R.; Skidmore, A.K.; van Duren, I. Effects of Canopy Structural Variables on Retrieval of Leaf Dry Matter Content and Specific Leaf Area From Remotely Sensed Data. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 898–909. [Google Scholar] [CrossRef]
  81. Rallo, G.; Minacapilli, M.; Ciraolo, G.; Provenzano, G. Detecting crop water status in mature olive groves using vegetation spectral measurements. Biosyst. Eng. 2014, 128, 52–68. [Google Scholar] [CrossRef]
  82. Roberts, D.A.; Roth, K.L.; Perroy, R.L. Hyperspectral Vegetation Indices. In Hyperspectral Remote Sensing of Vegetation; Thenkabail, P.S., Lyon, J.G., Huete, A., Eds.; CRC Press: Boca Raton, FL, USA, 2016; p. 782. [Google Scholar]
  83. Khanna, S.; Santos, M.; Ustin, S.; Shapiro, K.; Haverkamp, P.; Lay, M. Comparing the Potential of Multispectral and Hyperspectral Data for Monitoring Oil Spill Impact. Sensors 2018, 18, 558. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Grant, L. Diffuse and specular characteristics of leaf reflectance. Remote Sens. Environ. 1987, 22, 309–322. [Google Scholar] [CrossRef]
  85. Elsayed, S.; Mistele, B.; Schmidhalter, U. Can changes in leaf water potential be assessed spectrally? Funct. Plant Biol. 2011, 38, 523–533. [Google Scholar] [CrossRef] [PubMed]
  86. Seelig, H.-D.; Hoehn, A.; Stodieck, L.S.; Klaus, D.M.; Adams, W.W., III; Emery, W.J. The assessment of leaf water content using leaf reflectance ratios in the visible, near-, and short-wave-infrared. Int. J. Remote Sens. 2008, 29, 3701–3713. [Google Scholar] [CrossRef]
  87. Morell, F.J.; Álvaro-Fuentes, J.; Lampurlanés, J.; Cantero-Martínez, C. Soil CO2 fluxes following tillage and rainfall events in a semiarid Mediterranean agroecosystem: Effects of tillage systems and nitrogen fertilization. Agric. Ecosyst. Environ. 2010, 139, 167–173. [Google Scholar] [CrossRef] [Green Version]
  88. Cicuéndez, V.; Litago, J.; Huesca, M.; Rodriguez-Rastrero, M.; Recuero, L.; Merino-de-Miguel, S.; Palacios-Orueta, A. Assessment of the gross primary production dynamics of a Mediterranean holm oak forest by remote sensing time series analysis. Agrofor. Syst. 2015, 89, 491–510. [Google Scholar] [CrossRef]
  89. Xu, L.; Baldocchi, D.D. Seasonal variation in carbon dioxide exchange over a Mediterranean annual grassland in California. Agric. For. Meteorol. 2004, 123, 79–96. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Experimental area and plots location (maps and orthoimage from IGN, www.ign.es/ign_iberpix/) and overview picture (13 May 2013).
Figure 1. Experimental area and plots location (maps and orthoimage from IGN, www.ign.es/ign_iberpix/) and overview picture (13 May 2013).
Remotesensing 12 02724 g001
Figure 2. Scheme of the spectro-radiometric measurements (six angles from a reference point “RP”) and soil respiration collars distribution (5 circles) in the cultivated plots.
Figure 2. Scheme of the spectro-radiometric measurements (six angles from a reference point “RP”) and soil respiration collars distribution (5 circles) in the cultivated plots.
Remotesensing 12 02724 g002
Figure 3. Leaf area index (LAI) and soil respiration (Rs) average values measuredin the BY1 and BY2 barley plotsand in the bare soil plot (S) in 2012 (a) and in 2013 (b) in Madrid, Spain.
Figure 3. Leaf area index (LAI) and soil respiration (Rs) average values measuredin the BY1 and BY2 barley plotsand in the bare soil plot (S) in 2012 (a) and in 2013 (b) in Madrid, Spain.
Remotesensing 12 02724 g003
Figure 4. Spectral MODIS indices dynamics: (a) NDVI, (b) EVI, (c) NDWI, (d) ANIR, (e) SASI, (f) AG, (g) AR. Average values for the two barley plots (BY) and measured values for the bare soil plot (S) in 2012 and in 2013 in Madrid, Spain.
Figure 4. Spectral MODIS indices dynamics: (a) NDVI, (b) EVI, (c) NDWI, (d) ANIR, (e) SASI, (f) AG, (g) AR. Average values for the two barley plots (BY) and measured values for the bare soil plot (S) in 2012 and in 2013 in Madrid, Spain.
Remotesensing 12 02724 g004aRemotesensing 12 02724 g004bRemotesensing 12 02724 g004c
Figure 5. Spectral Sentinel-2 indices dynamics: (a) NDVI, (b) EVI, (c) NDWI2, (d) AG, (e) AR. Average values for the two barley plots (BY) and measured values for the bare soil plot (S) in 2012 and in 2013 in Madrid, Spain.
Figure 5. Spectral Sentinel-2 indices dynamics: (a) NDVI, (b) EVI, (c) NDWI2, (d) AG, (e) AR. Average values for the two barley plots (BY) and measured values for the bare soil plot (S) in 2012 and in 2013 in Madrid, Spain.
Remotesensing 12 02724 g005aRemotesensing 12 02724 g005b
Figure 6. Soil moisture (H10) and Soil temperature (T10) dynamics at 10 cm depth. Average values for the two barley plots (BY) and measured values for the bare soil plot (S) in 2012 (a) and in 2013 (b) in Madrid, Spain.
Figure 6. Soil moisture (H10) and Soil temperature (T10) dynamics at 10 cm depth. Average values for the two barley plots (BY) and measured values for the bare soil plot (S) in 2012 (a) and in 2013 (b) in Madrid, Spain.
Remotesensing 12 02724 g006
Table 1. Day of year (DOY) of vegetative and reproductive stages based on Zadoks’s scale [62].
Table 1. Day of year (DOY) of vegetative and reproductive stages based on Zadoks’s scale [62].
Growth StageDOY 2012DOY 2013
Vegetative stagesS7081
SG (10)8298
SG (12)88102
SG (13)-108
T1100113
T2108123
SE-133
B130143
Flowering and reproductive stagesF139155
D145171
R166184
H205198
S: Sowing; SG (10): Seedling growth, first leaf through coleoptiles; SG (12): Seedling growth, two leaves emerged; SG (13): Seedling growth, three leaves emerged; T1: Tillering, first stages of tillering; T2: Tillering, last stages of tillering; SE: Stem elongation; B: Booting; F: Flowering, inflorescence emergence; D: Dough development; R: Ripening; H: Harvest.
Table 2. Spectral indices obtained based on MODIS and Sentinel-2 wavelengths.
Table 2. Spectral indices obtained based on MODIS and Sentinel-2 wavelengths.
IndexEquationDescriptionReference
NDVI ( ρ n i r ρ r e d ) ( ρ n i r + ρ r e d ) Ratio Index based on the chlorophyll absorption at red band related to photosynthesis, vegetation growth and activity[39]
EVI ( ρ n i r ρ r e d ) ( ρ n i r + C 1 ρ r e d C 2 ρ b l u e + 1 ) Ratio Index developed to optimize the vegetation signal, especially in high biomass regions[40]
NDWI * ( ρ n i r ρ s w i r 1 ) ( ρ n i r + ρ s w i r 1 ) Ratio Index for detecting water status, the reflectance properties of green
vegetation, dry vegetation and soils
[64]
AGAngle at green ( α g r e e n ) SSI that takes into account three consecutive bands: blue, green and red. Related to pigment content and photosynthesis[48]
ARAngle at red ( α r e d ) SSI that takes into account three consecutive bands: green, red and NIR Related to pigment content and photosynthesis[48]
ANIR (only for MODIS)Angle at NIR ( α n i r ) SSI that takes into account three consecutive bands: red, NIR and SWIR1. Related to photosynthetic capacity and moisture and is able to detect dry plant matter in the presence of soil and green vegetation[52]
SASI (only for MODIS)Angle at SWIR1
( α s w i r 1 S l o p e )
SSI that takes into account three consecutive bands: NIR, SWIR1 and SWIR2. Related to moisture and non-photosynthetic activity, it is able to discriminate dry plant matter in the presence of soil and green vegetation[52]
ρ b l u e : blue spectral band in MODIS and band 2 in Sentinel-2; ρ r e d : red spectral band in MODIS and band 4 in Sentinel-2; ρ n i r : near-infrared spectral band in MODIS and band 8 in Sentinel-2; ρ s w i r 1 : shortwave infrared spectral band in MODIS. This band is not available in Sentinel-2 at 10 or 20 m of spatial resolution; S l o p e =   ( ρ n i r ρ s w i r 2 ) ( λ n i r ( 855.121   n m ) λ s w i r 2 ( 1627.713   n m ) ) ; * NDWI was calculated with ρ s w i r 2 (band 11) in Sentinel-2 obtaining NDWI 2.
Table 3. Names of the variables used in the plots.
Table 3. Names of the variables used in the plots.
VariableCultivated PlotsBare Soil
Soil RespirationRs_BYyearRs_Syear
IndexIndex_BYyearIndex_Syear
Leaf area indexLAI_BYyear
Soil moistureH10_BYyearH10_Syear
Soil temperatureT10_BYyearT10_Syear
Table 4. Coefficients of determination (Radj2) of a linear regression model for both years between LAI and MODIS spectral indices (LAI = a + b ∗ Index) for the vegetative (n = 22) and the reproductive stages (n = 18) in a rainfed barley (Hordeum vulgare L.) in Madrid, Spain.
Table 4. Coefficients of determination (Radj2) of a linear regression model for both years between LAI and MODIS spectral indices (LAI = a + b ∗ Index) for the vegetative (n = 22) and the reproductive stages (n = 18) in a rainfed barley (Hordeum vulgare L.) in Madrid, Spain.
MODIS IndicesAGARANIRSASINDVINDWIEVI
Vegetative stages
LAI = a + b ∗ index0.640.660.800.840.850.890.83
Reproductive stages
LAI = a + b ∗ index0.190.190.450.610.220.530.27
Significance of model’s coefficients was evaluated by the t-statistic (p < 0.05).
Table 5. Coefficients of determination (Radj2) of a linear regression model for both years between LAI and Sentinel-2 spectral indices (LAI = a + b ∗ Index) for the vegetative (n = 22) and the reproductive stages (n = 18) in a rainfed barley (Hordeum vulgare L.) in Madrid, Spain.
Table 5. Coefficients of determination (Radj2) of a linear regression model for both years between LAI and Sentinel-2 spectral indices (LAI = a + b ∗ Index) for the vegetative (n = 22) and the reproductive stages (n = 18) in a rainfed barley (Hordeum vulgare L.) in Madrid, Spain.
Sentinel-2 IndicesAGARNDVIEVINDWI 2
Vegetative stages
LAI = a + b ∗ index0.640.650.840.82 *0.87
Reproductive stages
LAI = a + b ∗ index0.15 I0.16 I0.200.240.63
Significance of model’s coefficients was evaluated by the t-statistic (p< 0.05). *: Intercept not significant. I: Index not significant.
Table 6. Coefficients of determination (Radj2) of simple and multiple linear regression models for Rs, as a function of LAI, MODIS spectral indices, soil temperature (T10) and soil moisture (H10) for both years for the vegetative (n = 22) and reproductive stages (n = 18) in a rainfed barley crop (Hordeum vulgare L.) in Madrid, Spain.
Table 6. Coefficients of determination (Radj2) of simple and multiple linear regression models for Rs, as a function of LAI, MODIS spectral indices, soil temperature (T10) and soil moisture (H10) for both years for the vegetative (n = 22) and reproductive stages (n = 18) in a rainfed barley crop (Hordeum vulgare L.) in Madrid, Spain.
LAI and MODIS IndicesLAIAGARANIRSASINDVINDWIEVI
Vegetative stages
Rs = a + b ∗ Index0.420.560.570.430.460.520.440.53
Rs = a + b ∗ Index + c ∗ T100.64 *0.640.640.610.64 *0.66 *0.59 *0.66 *
Rs = a + b ∗ Index + c ∗ H10--------
Rs = a + b ∗ Index + c ∗ T10 + d ∗ H10-0.71 *0.71 *0.69 *0.70 *0.71-0.73 *
Reproductive stages
Rs = a + b ∗ Index0.62 *0.550.550.710.780.580.750.64
Rs = a + b ∗ Index + c ∗ T100.72-------
Rs = a + b ∗ Index + c ∗ H100.85 *-------
Rs = a + b ∗ Index + c ∗ T10 + d ∗ H10--------
Significance of model’s coefficients was evaluated by the t-statistic (p < 0.05). *: Intercept not significant. -: model or variables not significant.
Table 7. Coefficients of determination (Radj2) of simple and multiple linear regression models for Rs, as a function of LAI, Sentinel-2 spectral indices, soil temperature (T10) and soil moisture (H10) for both years for the vegetative (n = 22) and reproductive stages (n = 18) in a rainfed barley crop (Hordeum vulgare L.) in Madrid, Spain.
Table 7. Coefficients of determination (Radj2) of simple and multiple linear regression models for Rs, as a function of LAI, Sentinel-2 spectral indices, soil temperature (T10) and soil moisture (H10) for both years for the vegetative (n = 22) and reproductive stages (n = 18) in a rainfed barley crop (Hordeum vulgare L.) in Madrid, Spain.
LAI and Sentinel-2 IndicesLAIAGARNDVIEVINDWI 2
Vegetative stages
Rs = a + b ∗ Index0.420.570.580.520.540.47
Rs = a + b ∗ Index + c ∗ T100.64 *0.630.640.66 *0.66 *0.64 *
Rs = a + b ∗ Index + c ∗ H10------
Rs = a + b ∗ Index + c ∗ T10 + d ∗ H10-0.70 *0.71 *0.710.730.70 *
Reproductive stages
Rs = a + b ∗ Index0.62 *0.500.500.560.610.78
Rs = a + b ∗ Index + c ∗ T100.72-----
Rs = a + b ∗ Index + c ∗ H100.85 *-----
Rs = a + b ∗ Index + c ∗ T10 + d ∗ H10------
Significance of model’s coefficients was evaluated by the t-statistic (p < 0.05). *: Intercept not significant. -: model or variables not significant.
Table 8. Coefficients of determination (Radj2) of the simple and multiple linear regression models between soil respiration(Rs) and soil temperature (T10) and soil moisture (H10) at 10 cm depth for both years for the vegetative (n = 22) and reproductive stages (n = 18) and Radj2 of the simple and multiple linear regression models between Rs in bare soil (Rh) and T10 and H10 for 2012 and 2013, for the vegetative (n = 11) and reproductive stages (n = 9) separately in a rainfed barley (Hordeum vulgare L.) in Madrid, Spain.
Table 8. Coefficients of determination (Radj2) of the simple and multiple linear regression models between soil respiration(Rs) and soil temperature (T10) and soil moisture (H10) at 10 cm depth for both years for the vegetative (n = 22) and reproductive stages (n = 18) and Radj2 of the simple and multiple linear regression models between Rs in bare soil (Rh) and T10 and H10 for 2012 and 2013, for the vegetative (n = 11) and reproductive stages (n = 9) separately in a rainfed barley (Hordeum vulgare L.) in Madrid, Spain.
T10H10T10+H10
Vegetative stages
Rs0.22 *00.23 *,H
Reproductive stages
Rs0.02 *0.380.35 *,T
Bare soil (Rh)
Rh (vegetative stages)0.41 *0.02 *,H0.39 *,H
Rh (reproductive stages)0.25 T0.52 *0.57 *,T
Significance of model’s coefficients was evaluated by the t-statistic (p < 0.05). *: Intercept not significant. T: T10 not significant. H: H10 not significant.

Share and Cite

MDPI and ACS Style

Cicuéndez, V.; Rodríguez-Rastrero, M.; Recuero, L.; Huesca, M.; Schmid, T.; Inclán, R.; Litago, J.; Sánchez-Girón, V.; Palacios-Orueta, A. First Insights on Soil Respiration Prediction across the Growth Stages of Rainfed Barley Based on Simulated MODIS and Sentinel-2 Spectral Indices. Remote Sens. 2020, 12, 2724. https://doi.org/10.3390/rs12172724

AMA Style

Cicuéndez V, Rodríguez-Rastrero M, Recuero L, Huesca M, Schmid T, Inclán R, Litago J, Sánchez-Girón V, Palacios-Orueta A. First Insights on Soil Respiration Prediction across the Growth Stages of Rainfed Barley Based on Simulated MODIS and Sentinel-2 Spectral Indices. Remote Sensing. 2020; 12(17):2724. https://doi.org/10.3390/rs12172724

Chicago/Turabian Style

Cicuéndez, Víctor, Manuel Rodríguez-Rastrero, Laura Recuero, Margarita Huesca, Thomas Schmid, Rosa Inclán, Javier Litago, Víctor Sánchez-Girón, and Alicia Palacios-Orueta. 2020. "First Insights on Soil Respiration Prediction across the Growth Stages of Rainfed Barley Based on Simulated MODIS and Sentinel-2 Spectral Indices" Remote Sensing 12, no. 17: 2724. https://doi.org/10.3390/rs12172724

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