
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

71430
10.1038/s41598-024-71430-2
Article
Assessing soil CO2 emission on eucalyptus species using UAV-based reflectance and vegetation indices
Rossi Fernando Saragosa 1
Della-Silva João Lucas 2
Teodoro Larissa Pereira Ribeiro 3
Teodoro Paulo Eduardo eduteodoro@hotmail.com

3
Santana Dthenifer Cordeiro 3
Baio Fábio Henrique Rojo 3
Morinigo Wendel Bueno 4
Crusiol Luís Guilherme Teixeira 5
La Scala Newton Jr. 1
da Silva Carlos Antonio Jr. 6
1 grid.410543.7 0000 0001 2188 478X PPG-Ciência do Solo, State University of São Paulo (UNESP), Jaboticabal, 15385-000 Brazil
2 PPG-Bionorte, State University of Mato Grosso (UNEMAT), Sinop, 78550-000 Brazil
3 https://ror.org/0366d2847 grid.412352.3 0000 0001 2163 5978 Federal University of Mato Grosso do Sul (UFMS), Chapadão do Sul, 79560-000 Brazil
4 https://ror.org/01mqvjv41 grid.411206.0 0000 0001 2322 4953 PPGCAM, Federal University of Mato Grosso (UFMT), Sinop, 78555000 Brazil
5 https://ror.org/0482b5b22 grid.460200.0 0000 0004 0541 873X National Soybean Research Center (Embrapa Soja), Brazilian Agricultural Research Corporation, Londrina, 86001–970 Brazil
6 Department of Geography, State University of Mato Grosso (UNEMAT), Sinop, Mato Grosso 78555000 Brazil
31 8 2024
31 8 2024
2024
14 2027721 11 2023
28 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Eucalyptus species play an important role in the global carbon cycle, especially in reducing the greenhouse effect as well as storing atmospheric CO₂. Thus, assessing the amount of CO₂ released by the soil in forest areas can generate important information for environmental monitoring. This study aims to verify the relation between soil carbon dioxide (CO₂) flux (FCO₂), spectral bands, and vegetation indices (VIs) derived from a UAV-based multispectral camera over an area of eucalyptus species. Multispectral imageries (green, red-edge, and near-infrared) from the Parrot Sequoia sensor, derived vegetation indices, and the FCO₂ data from a LI-COR 8100 analyzer, combined with soil moisture and temperature data, were collected and related. The vegetation indices ATSAVI (Adjusted Transformed Soil-Adjusted VI), GSAVI (Green Soil Adjusted Vegetation Index), and SAVI (Soil-Adjusted Vegetation Index), which use soil correction factors, exhibited a strong negative correlation with FCO₂ for the species E. camaldulensis, E. saligna, and E. urophylla species. A Multivariate Analysis of Variance showed significance (p < 0.01) for the species factor, which indicates that there are differences when considering all variables simultaneously. The results achieved in this study show a specific correlation between the data of soil CO₂ emission and the eucalypt species, providing a distinction of values between the species in the statistical data.

Subject terms

Plant sciences
Environmental sciences
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The increasing demand for timber over the last few years has led to land use change to areas dedicated to tree crops with rapid growth1. In Brazil, the total area dedicated to planted forest is about 9.5 million hectares (ha), 76.9% of which are eucalyptus (Eucalyptus spp.) cash crops2. This genus is distinguished by the availability of cultivars adapted to the different Brazilian biomes and edaphoclimatic conditions3. Eucalyptus monoculture fields are distributed throughout the tropics and subtropics, notably in Brazil4. Mato Grosso do Sul state has been the eucalyptus expansion head, presenting a significant increase in cellulose production in the last five years, making this state the leading exporter in the first four months of 20205.

Eucalyptus has a greater potential to stock atmospheric carbon in aerial biomass or soil, mainly when associated with pasture or annual crops and especially in the conversion of degraded land to productive land and for renewable energy sources6. Thus, forest crops play an important role in the global carbon cycle, especially in reducing the greenhouse effect as well as stocking atmospheric CO₂7. Monitoring the amount of CO₂ released by the soil in forest areas can generate important information for cultivation, being influenced by the type of management and the effects of human actions on carbon emission, as well as climatic variation, since soil respiration is the second largest flux of carbon between terrestrial ecosystems and the atmosphere8. Soil CO₂ emission (FCO₂) in eucalyptus crops is affected by environmental factors, such as soil temperature, water content, and precipitation variations9.

Throughout the carbon cycle, CO₂ flow is based on plant (autotrophic) and microbial respiration (heterotrophic). Soil respiration comprises the autotrophic respiration below the soil surface where plant roots are the main contributors, and the heterotrophic respiration of microorganisms during the decomposition of organic matter10–12. The challenges of forest management rely on reliable, accurate, and cost-effective methods to adequately report forest dynamics, especially because of the growing carbon trading (e.g., reducing emissions from deforestation and forest degradation)13. In this sense, remote sensing techniques may assist in understanding carbon dynamics comprehension and the role of greenhouse gases (GHGs) in the soil-atmosphere interface, for instance through multispectral image analysis in conjunction with in situ measurements for the monitoring, mapping, and enforcement of vegetation cover14,15.

Significant changes have occurred in forestry activities management, considering the advancement of big data, geoprocessing, and informatics technologies16, as well as the use of remotely piloted drones that use images for crop mapping, evaluation of cultivated areas, disease detection, and soil mapping17. Thus, image analysis tools through spectral indices are presented as alternatives for further research and decision-making more objectively. There is a relationship between the amount of cumulative FCO₂ and the growth of eucalyptus trees since these plants need the nutrients released into the soil after decomposition of the existing organic material18. From the perspective of remote sensing techniques, vegetation indices are capable of diagnosing certain factors related to eucalyptus plants, such as the amount of water in the leaf19, leaf area index20, and distinguishing species by reflectance and growth16. Thus, this manuscript hypothesizes that differences in FCO₂ may indicate the existence of a correlation with the spectral behavior of eucalyptus species and that the UAV-based image sensor may be a proxy for assessing CO₂ flux. The principal contributions of this paper can be summarized as follows:Evaluation of the relationship between soil carbon dioxide flux (FCO₂) and UAV-based multispectral data.

Assessment of the correlation between reflectance values and derived vegetation indices with the spectral behavior of eucalyptus species.

Investigation of whether differences in FCO₂ indicate a correlation with the spectral behavior of eucalyptus species.

Testing of the hypothesis that UAV-based image sensors can serve as proxies for assessing CO₂ flux in eucalyptus species.

Results

Table 1 presents the descriptive statistics for the variables FCO₂, Ms, and Ts, concerning the eucalyptus species. The mean value obtained for FCO₂ was 0.54 µmol m−2 s−1, with minimum and maximum values ranging from 0.1 to 1.65 µmol m−2 s−1. For Ms and Ts, the mean values were 8.28 and 20.9, respectively. The coefficient of variation was below 10% for Ts and less than 30% for Ms. However, for FCO₂, this value was 56%. Therefore, only temperature ex-hibited low variability (CV < 10%), remaining consistent with the minimum and maximum val-ues, while the other variables displayed moderate to high variability in their values.Table 1 Descriptive statistics of soil emissions (FCO₂), soil moisture (Ms), and soil temperature (Ts) in relation to six eucalyptus species in the year 2021.

	FCO₂	Ms	Ts	
Mean	0.54	8.28	20.9	
Median	0.46	8.0	21.0	
Maximum	1.65	15.0	24.0	
Minimum	0.1	3.0	18.0	
Range	1.55	12.0	6.2	
Standard deviation	0.3	2.4	1.7	
Coefficient of variation (%)	56	29	8	

Figure 1 depicts the boxplot for the statistical analysis of data distribution for the six eucalyptus species concerning the variables FCO₂, Ms, and Ts. The asymmetric position of the median relative to the mean symbolizes an association with non-parametric data, which are closer to both the lower quartile (Q1) and upper quartile (Q3) in both representations. The eucalyptus species exhibit different ranges of values for the measurements and the presence of outliers for the species E. camaldulensis and E. urograndis in the case of the FCO₂ variable (Fig. 1a), and Corymbia citriodora and E. grandis for Ms, represented by the points that deviate from the box edges.Fig. 1 Boxplot of soil CO₂ emissions (FCO₂), soil moisture (Ms) and soil temperature (Ts) values concerning the eucalyptus species.

Among the species evaluated for the FCO₂ variable, E. saligna and E. urophylla stand out for having values distributed more closely around the mean, meaning there is a proximity be-tween the median and mean values. This occurrence is also observed for the Ms variable (Fig. 1b) and Ts variable (Fig. 1c) for the species E. saligna, E. urograndis, E. urophylla, and E. saligna, respectively.

For the species E. camaldulensis and E. urograndis (FCO₂), and C. citriodora and E. grandis (Ms), the median is below the mean, possibly due to the outliers that negatively impacted the mean. In general, the interquartile ranges, representing the interval between the second and third quartiles, are larger for all variables than the interquartile range between the first and second quartiles, indicating greater data dispersion among the top 25% of observed values.

Pearson's correlation analysis (Figure S1) was high and positive between the features derived from the aggregated photosynthetic processes. The green, red, red-edge, and NIR (near infrared) spectral bands had a negative linear correlation with FCO₂. The GEMI (Global Environment Monitoring Index), ATSAVI (Adjusted Transformed soil-adjusted VI), SAVI (Soil-Adjusted Vegetation Index), EVI2 (Enhanced Vegetation Index 2), GSAVI (Green Soil Adjusted Vegetation Index), GDVI (Difference NIR/Green Difference Vegetation Index), CIrededge (Chlorophyll Index RedEdge), CIgreen (Chlorophyll Index Green) and BWDRVI (Blue-wide dynamic range vegetation index) indices showed the highest correlation with the characteristics FCO₂. The FCO₂ variable showed links with all variables, but no variable was clearly associated.

Considering the large number of variables acquired via remote sensing and assigned in situ to the eucalyptus species applied in Pearson’s correlation analysis, the variables GEMI, ATSAVI, SAVI, EVI2, GSAVI, GDVI, CIrededge, CIgreen and BWDRVI were selected, which had the highest correlations with FCO₂.

Multivariate analysis of variance for eucalyptus species, the variables FCO₂, Ms, Ts, the spectral bands (green, red, NIR and red edge) and vegetation indices (GEMI, ATSAVI, SAVI, EVI2, GSAVI, GDVI, CIrededge, CIgreen and BWDRVI) were used to test the hypothesis of the difference between the eucalypt species. The significance level for the species considering all variables simultaneously was ρ =  < 0.001 for Pillai's test, showing that the null hypothesis is rejected, as there are differences among eucalyptus species.

The canonical analysis considers the eucalyptus species and the variables with a strong correlation to FCO₂ (Fig. 2). The ratio of the variance accumulated in the first two canonical variables was higher than 61.7% (Can1 and Can2). The variables FCO₂, Ms, GSAVI, GDVI, CIrededge, NIR, rededge and CIgreen were most important in distinguishing eucalyptus species in Can1, while part red, green, BWDRVI, Ts, GEMI, ATSAVI, SAVI and EVI2 ratio were the most relevant variables in Can2. This finding reveals that a significant portion of the variability of the eucalypt species concerning the variables evaluated in each species can be summarized in two linear combinations, named Can1 and Can2.Fig. 2 Two-dimension plot of Can1 versus Can2 scores generated by the canonical variable analysis considering eucalyptus species and FCO2, soil moisture (Ms) and temperature (Ts), spectral bands (green, red, NIR, and red-edge), and vegetation indices (GEMI, ATSAVI, SAVI, EVI2, GSAVI, GDVI, CIrededge, CIgreen and BWDRVI).

From the analyzed variables in Fig. 2, it can be observed that red, rededge, ATSAVI, GEMI and FCO₂ presented the greatest variability, which can be proven by the size of the module of their vectors. There are high correlations between the variables FCO₂ and Ms; red edge, NIR, CIrededge, GSAVI, GDVI, ATSAVI, SAVI and EVI2; Ts, green and BWDRVI for forming acute angles between the variables. There is no other correlation combination between the variables, as it forms an angle close to 90 degrees, as shown in Fig. 2. The set of variables FCO₂ and Ms; CIrededge and red-edge; red; GEMI were those that contributed most to the discrimination of the species E. camaldulensis; E. grandis and E. urophylla, E. saligna and C. citriodora, respectively. On the other hand, the species E. urophylla and E. grandis were more similar in the canonical variables analysis.

In Fig. 3, we exclusively present the zoomed-in view of the direct correlation between the variables and FCO₂, which encompasses the correlations established among the species of eucalyptus species showed specific correlations with FCO₂. Statistical significance is considered at p ≤ 0.05 by *, p ≤ 0.01 by ** and p ≤ 0.001 by ***. The FCO₂ correlation to Ms (0.243*), EVI2 (− 0.228*), red (− 0.260*), SAVI (− 0.233*), CIrededge (− 0.327**), GEMI (0.263**), GDVI (− 0.294**), GSAVI (− 0.286**), Ts (− 0.386***), BWDRVI (− 0.390***), CIgreen (− 0.412***), green (− 0.375***), NIR (− 0.335***) and rededge (− 0.352***), with specific correlation among E. camaldulensis (0.514*) and E. grandis (0.487*); E. camaldulensis (− 0.678**) and E. saligna (0.576*); E. saligna (− 0.861***), E. urograndis (− 0.638**) and E. urophylla (− 0.514*); E. camaldulensis (− 0.697**) and E. saligna (0.599*); E. camaldulensis (− 0.548*) and E. urograndis (− 0.685**); E. saligna (0.860***), E. urograndis (0.639**) and E. urophylla (0.514*); E. camaldulensis (− 0.648**) and E. urograndis (− 0.634**); E. camaldulensis (− 0.692**), E. urograndis (− 0.556*) and E. urophylla (− 0.528*); E. camaldulensis (− 0.526*); E. saligna (− 0.786***), E. urograndis (− 0.685**) and E. urophylla (− 0.591*); E. camaldulensis (− 0.547*), E. saligna (− 0.655**), E. urograndis (− 0.713**) and E. urophylla (− 0.554*); E. saligna (− 0.886***), E. urograndis (− 0.677**) and E. urophylla (− 0.505*); E. camaldulensis (− 0.636**), E. urograndis (− 0.697**) and E. urophylla (− 0.517*); E. camaldulensis (− 0.506*), E. saligna (− 0.528*) and E. urograndis (− 0.674**) respectively were significant at 99%. For the variable ATSAVI, there was no correlation with FCO₂, however, it displayed specific correlations with E. camaldulensis (− 0.728**) and E. saligna (0.710**).Fig. 3 Heatmap of correlations between eucalyptus species and FCO2 variables as per the modeling of canonical variable analysis.

C. citriodora had no significance for the correlation to FCO₂ in the entire set of variables. The species most often correlated significantly with FCO₂ were E. camaldulensis and E. saligna.

Discussion

An absence of strong and significant correlations between FCO₂ and the rest of the variables in a general set with eucalyptus species is noticeable (Figure S1). However, differences were identified among the eucalyptus species while simultaneously considering the soil variables with the main vegetation indices.

As there is a low and negative correlation between the variables and FCO2, we can suggest that the variables move in opposite directions, which means that when the values of a variable increase, FCO2 tends to decrease, and therefore other factors may be influencing this relationship, since leaves are the most important organs for the spectral characterization of vegetation, reflecting the conditions of the plant21.

The interactions of a species during its cycle with the environment are known to be complex, with many aspects not being sufficiently clarified22. This may explain the behavior of the C. citriodora species response when correlated with the FCO2 variable, which is different from what was found for the other eucalyptus species. These measurements and comparisons provide a better understanding of plant responses to multispectral variables. As a fast-growing species, eucalyptus plays an important role in mitigating climate change due to its ability to absorb carbon dioxide23.

Scientific research based on field data and remote sensing, as evidenced in previous studies, Teodoro et al.24, Oliveira et al.16, and Xu et al.25 are important tools for future planning and decision-making. These tools are fundamental in reducing climate risks that directly affect biodiversity, agriculture, water resources, coastal zones, etc., as well as promoting sustainable development by allowing forestry practices to be carried out more safely.

Total soil respiration comprises autotrophic and heterotrophic respiration26,27, which are directly and indirectly related to biotic and abiotic factors28,29 for the regulation of an ecosystem. The behavior of autotrophic respiration is strongly correlated with photosynthesis30,31 and is mediated by phenological dynamics throughout the plant cycle32,33. Heterotrophic respiration from the mineralization of organic matter by the decomposer community is highly affected by soil carbon availability, moisture and temperature34,35.

Ms and Ts are important for the respiration process of soils and directly influence the decomposition of organic matter and microbial and root activities36,37. Thus, Ms showed a positive correlation and Ts showed a negative correlation with FCO₂. However, this low-magnitude correlation between Ms and Ts with FCO₂ can be explained by the low variation over the time series period38. The effects of Ts on FCO₂ may be partially hidden by the effect of Ms, since these are interdependent variables and commonly change simultaneously39.

The results obtained by Almeida et al.40, Dossou-Yovo et al.41 and Pires et al.42 show that the emission process and the transport of the gas from the soil interior to the atmosphere depend on the porosity, the moisture content and the conditions of the whole soil biota. These conditions in the ecosystem are at the same time very favorable for the activity of microorganisms, which can positively influence respiration rates43. In addition, management practices will influence the behavior of carbon storage capacity and the loss of CO₂ from the soil to the atmosphere44.

In the case of data obtained via passive remote sensors, the spectral bands green, NIR, red, and red-edge were negatively correlated and of low magnitude with FCO₂, but with significance in E. camaldulensis, E. saligna, and C. citriodora (Fig. 3). The low association between spectral bands and FCO₂ is directly associated with the green vegetative canopy and also related to the stage of development of the eucalyptus species that influences leaf thickness, water content and response of leaf chlorophyll contents. The chlorophyll a and b pigmentations are the ones that most influence the electromagnetic radiation in the visible region, with two absorption peaks, the higher in the red band and the lower in the blue band, respectively45,46, favoring a better result for the classifier that used the red (R) and blue (B) spectral bands. From Silva et al.47, the green wavelength response is entirely associated with the potential of photosynthetic activity, and this relation is directly related to the utilization of available radiation by these chlorophyll pigments. The process of CO₂ absorption from the atmosphere by the leaves occurs when the pores of the stomas open, allowing the carbon dioxide to enter the leaf, mainly located in the lower epidermis48,49.

The use of the red and near-infrared (NIR) spectral bands in the composition of vegetation indices is because these wavelengths are inversely proportional to leaf reflectance, with greater than 90% sensitivity on the spectral behavior of a vegetation canopy50. BWDRVI and GEMI uses blue and NIR spectral reflectance to measure the density and greenness of vegetation, showing a specific negative and positive correlation with E. saligna, E. urograndis and E. uroplylla, respectively. The blue spectral band is poorly absorbed by the upper part of the plant canopy for the process of photosynthesis, the lower levels of the canopy do not contribute to reflectance and scattering of electromagnetic radiation, which impairs the correlation of BWDRVI and GEMI with FCO₂.

The absorption of the red radiation is by the photosynthesizing pigments of the leaves, which are the chlorophyll ‘a’, while in the near-infrared the absorption is almost zero, and the reflectance and transmittance are high, due to the internal cellular structures of the leaves and the effect of the multiple reflectances51. In previous studies, the ratio between leaf chlorophyll and carotenoid can be considered as a proxy to track the physiological and phenological status of vegetation52, but its non-obvious absorption characteristics in the visible spectral regions53. To capture variations in the amount of chlorophyll not only in the upper canopy but also throughout the foliage, the EVI2, ATSAVI and SAVI index were calculated using a combination of the NIR and red-edge bands. A negative, positive, and specific negative correlation with E. camaldulensis and E. saligna, respectively, differed from the correlations of the EVI2, ATSAVI and SAVI index, which can be explained by the low values of chlorophyll content.

SAVI, GSAVI and ATSAVI use a soil brightness correction factor (L) applied to maximize the reduction of soil effects on the vegetation signal and showed negative and positive correlations with FCO₂ for E. camaldulensis and E. saligna, respectively. The constant L value was 0.5, recommended for a wide range of vegetation scenarios54, but using the value of L = 0.5 in other studies resulted in a higher standard deviation in pixel values55. SAVI, GSAVI and ATSAVI put down general vegetation characteristics by correcting for ground brightness, and is a disadvantage when vegetation is at peak vegetation56. However, the presence of intact dead organic matter on the soil surface decreased the effect of precipitation on FCO₂, as it has no effect on the soil structure and only during its decomposition, thus being a physical barrier against the abrupt entry of water into the soil57. The presence of decomposing dead plant matter forms aggregating and stabilizing substances in the soil structure, with beneficial effects on increasing soil carbon stocks through physical, chemical, and biological processes58,59. However, there is evidence that dead organic matter increases soil CO₂ emissions since straw can be rapidly broken down by soil microorganisms60,61. More research is needed, in order to validate the effects of straw on soil CO₂ emissions.

These results enable the management of differentiation among eucalyptus species using sensors coupled to remotely piloted drone and FCO₂, thus integrating new assumptions in this line of research. The reflectance of eucalyptus species acquired by remotely piloted aircraft makes it possible to obtain data about them in a fast, non-destructive, and accurate way. Measuring CO₂ provides important information about the carbon emissions in areas of eucalyptus. Both pieces of information contribute to the improvement of forest management. It is important to stress that new studies should be conducted in places with different soil and climate conditions from the current research.

Methods

Studied area

This study is conducted in the experimental field of an eucalyptus crop at the Federal University of Mato Grosso do Sul, Chapadão do Sul campus (Fig. 4). According to Köppen, the region is classified as a tropical monsoon climate (Am), with two well-defined seasons, the dry season from April to September and the wet season from October to March62. The annual temperature varies from 13 to 28 °C and presents pluviometric precipitation varying from 750 to 1850 mm63. The soil in this region is classified as dystrophic red latosol64. The plant collection and use following all the relevant guidelines.Fig. 4 Study area at the Federal University of Mato Grosso do Sul, Chapadão do Sul campus.

In situ experimental details

Eucalyptus species were monitored in January 2014 in a randomized block experimental design with four repetitions, with 20 plants within each experimental plot. Experimental treatments consisted of six eucalyptus species (E. camaldulensis, E. urophylla, E. saligna, E. grandis, E. urograndis and Corymbia citriodora) and the fertilization needs were determined from the soil chemical analysis, the results of which were obtained: pH (CaCl2): 4.9; organic matter: 31.5 g dm−3; phosphor: 13.6 m g dm−3; hydrogen + aluminum (H + Al): 5.4; potassium: 0.29 cmolc dm−3; calcium: 2.8 cmolc dm−3; magnesium: 0.5 cmolc dm−3; cation exchange capacity (CEC): 9.0 cmolc dm−3; base saturation: 39.9%. The proportions of clay, sand and silt were 46%, 46% and 8%, respectively. Application of herbicides (glyphosate) and pest control was carried out when necessary.

To evaluate the carbon dioxide (CO₂) efflux from the soil (FCO₂), we used the LI-COR portable system, model LI-8100, which monitors the CO₂ concentration variations inside the soil respiration chamber using optical absorption spectroscopy in the infrared spectral region. The soil breathing chamber is a closed system with an internal volume of 854.2 cm3 and a circular contact area of 83.7 cm2, coupled over the polyvinyl chloride (PVC) collars that were previously inserted into the soil 24 h before collection. The process of placing the PVC collars over the soil surface helps to avoid contamination of the samples by the camera upon contact with the soil, having minimal disturbance of the soil in the measurement process.

Four in situ readings were taken in each block of eucalyptus species, previously defined as equidistant and georeferenced, totaling 96 sample points, where the blocks were divided into quadrants. The readings were dependent on PVC collars for gas flux chamber disposal, where these collars were placed in each quadrant centroid. The FCO₂ acquisition was carried out in the morning, from 8 to 10 a.m. The FCO₂ was evaluated at each point by an adjustment of the CO₂ concentration of the inside chamber air as a function of a parabolic regression over the measured period. The measurement mode for the determination of soil CO₂ emission took 90 s at each of the sample points and the CO₂ concentration inside the chamber was measured every 2.5 s, approximately65.

The soil temperature (Ts) was measured using a Digital Spit-Type Thermometer (DELLT DT-625). It consists of a 20 cm rod that is inserted at least 50% into the soil at 5 cm from the location where the PVC collars were previously installed. Similarly, the Soil Moisture (Ms) was recorded using a TDR—Time Domain Reflectometry (Hydrosense TM, Campbell Scientific, Australia) device, consisting of a probe with two 12 cm rods that were embedded inside the soil at least 50%, perpendicular to the surface, also at 5 cm from the PVC collars. The soil moisture value is derived from the time it takes for an electric current to travel the distance of 32 mm from one rod to the other. The soil temperature and soil moisture assessments were carried out concomitantly with the assessments of soil CO₂ emission.

UAV-based multispectral image acquisition and processing

UAV-based multispectral images were obtained simultaneously with the in situ collections using the Parrot Sequoia multispectral sensor coupled to the remotely piloted fixed-wing drone Sensefly eBee RTK, with autonomous flight control (Fig. 5). The flight altitude over the study site was 100 m, providing a spatial resolution of 0.10 m. The Parrot Sequoia is a multispectral camera for agricultural applications, containing a sunlight sensor and an additional 16 Mpx RGB camera for recognition. The sensor has a horizontal field of view (HFOV) of 61.9°, vertical field of view (VFOV) of 48.5° and diagonal field of view (DFOV) of 73.7° and the delivered images are provided in the green (G—550 nm ± 40 nm), red (R—660 nm ± 40 nm), red-edge (RE—735 nm ± 10 nm) and near-infrared (NIR—790 nm ± 40 nm) spectral bands.Fig. 5 Remotely Piloted Fixed-Wing drone Sensefly eBee RTK loaded with Sensefly Sequoia sensor.

Flights were performed with an image lateral and frontal overlap of 80% and 85% respectively, with perpendicular flight lines over the same study area, enhancing the amount and quality of the further processing stages. Increasing the image overlap was necessary to obtain a greater number of scenes containing the same control points. This procedure enabled greater accuracy in the ortho mosaicing of images using Pix4Dmapper software due to the height of the eucalyptus plants and the intense movement of leaves because of the windy conditions. Flights were performed at the zenith in order to minimize shadows on the trees, at 11 a.m. (GMT -4).

The aerial survey was carried out by using RTK (Real Time Kinematics) technology, which in turn allows the estimation of the camera position making it possible to estimate the position of the camera at the instant of image collection with an accuracy of 2.5 cm. The images were mosaicked and orthorectified in the computer program Pix4Dmapper. Field calibration of the Parrot Sequoia was performed in the e-Motion software using images from the reflectance panel for each spectral band right before the UAV takeoff. As the flight duration was less than 15 min, there was no need to repeat the calibration after landing. The positional accuracy of the orthoimages was verified with ground control points (GCP), based on GNSS RTK Emilid Reach.

The processing of the multispectral images was performed using QGIS66 software in order to preliminarily analyze the spectral bands and vegetation indices, verifying if the values were correct. Vegetation indices (Table S1) were calculated and overlaid onto the multispectral images and shapefile of the eucalyptus species blocks/quadrant. For each spectral band and vegetation index, the average from all pixels in the blocks/quadrant (Fig. 4) was performed using the extract values by points tool in QGIS. To avoid the border effect between the repetition blocks, an internal buffer of one meter was removed from each block.

Statistical analysis

First, the data were submitted to descriptive statistics to evaluate the behavior of FCO₂ in the area of eucalyptus plantations and subsequently, an individualized analysis for each species to observe the interrelationships between them. The analyses were performed using R 3.6.367 software.

Spectral modelling for FCO2 estimation

Thus, the CO2 flux index model (µmol m−2 s−1) is the result of the multiplication of NDVI and sPRI, in which there is a relationship between the PRI index, which indicates light-use efficiency in photosynthesis, and NDVI, which indicates the vigor of photosynthetically active vegetation, in which it may be able to capture absorptions from carbon sequestration. Thus, the best correlation is given in Eq. 368–70.

With the purpose of a further examination of the main spectral and emissions of soil carbon dioxide related variables was performed using Pearson correlation (by using the “Corrplot” package)71, with a level of significance of p < 0.05 at a 95% confidence level. Subsequently, the variables that demonstrated a 95% confidence level with FCO₂ were selected, performing multivariate analysis of variance and canonical variables to investigate the interrelationship between eucalyptus species. For these analyses, the “candisc” package72 was used.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-71430-2.

Acknowledgements

The authors would like to thank the Federal University of Mato Grosso do Sul (UFMS), State University of Mato Grosso (UNEMAT), National Council for Scientific and Technological Development (CNPq) – Grant numbers 308295/2023-4, 309250/2021-8, 306022/2021-4 and 304979/2022-8, and Foundation for the Support of the Development of Education, Science and Technology of the State of Mato Grosso do Sul (FUNDECT) TO numbers 88/2021, 07/2022, 318/2022 and 94/2023, and SIAFEM numbers 30478, 31333, 32242 and 33111. This study was funded in part by the Coordination for the Improvement of Higher Education Personnel - Brazil (CAPES) – Financial Code 001.

Author contributions

Conceptualization, F.S.R., N.L.S.Jr and C.A.S.J.; Data curation, F.H.R.B.; Formal analysis, F.S.R., D.C.S. and P.E.T.; Visualization, J.L.S.S., L.P.R.T., W.B.M. and L.G.T.C.; Writing—original draft, F.S.R., N.L.S.Jr and C.A.S.J. All authors have read and agreed to the published version of the manuscript.

Data availability

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Viera M Rodríguez-Soalleiro R A complete assessment of carbon stocks in above and belowground biomass components of a hybrid eucalyptus plantation in Southern Brazil Forests 2019 10 536 10.3390/f10070536
Viera, M. & Rodríguez-Soalleiro, R. A complete assessment of carbon stocks in above and belowground biomass components of a hybrid eucalyptus plantation in Southern Brazil. Forests 10, 536 (2019).10.3390/f10070536
2. IBGE. Instituto Brasileiro de Geografia e Estatística. https://agenciadenoticias.ibge.gov.br/agencia-noticias/2012-agencia-de-noticias/noticias/35056-valor-de-producao-da-silvicultura-e-da-extracao-vegetal-cresce-27-1-e-chega-ao-recorde-de-r-30-1-bilhoes (2022).
3. Costa, T. C. C., Campanha, M. M. & Gontijo Neto, M. M. Sequestro de CO2 em Árvores de Eucalipto no Sistema ILPF. Embrapa Milho e Sorgo 1–30 https://www.infoteca.cnptia.embrapa.br/infoteca/bitstream/doc/1149131/1/Boletim-244-Sequestro-de-CO2-em-arvores-de-eucalipto-no-sistema-ILPF.pdf (2022).
4. Booth TH Eucalypt plantations and climate change For. Ecol. Manag. 2013 301 28 34 10.1016/j.foreco.2012.04.004
Booth, T. H. Eucalypt plantations and climate change. For. Ecol. Manag. 301, 28–34 (2013).10.1016/j.foreco.2012.04.004
5. Logística e Gestão Da Cadeia de Suprimentos Na Área Portuária—Volume 1. (Editora Poisson, 2021). 10.36229/978-65-5866-142-9.
6. La Scala Júnior N De Figueiredo EB Panosso A A review on soil carbon accumulation due to the management change of major Brazilian agricultural activities Braz. J. Biol. 2012 72 775 785 10.1590/S1519-69842012000400012 23011303
La Scala Júnior, N., De Figueiredo, E. B. & Panosso, A. A review on soil carbon accumulation due to the management change of major Brazilian agricultural activities. Braz. J. Biol. 72, 775–785 (2012).23011303 10.1590/S1519-69842012000400012
7. Du H Carbon storage in a Eucalyptus plantation chronosequence in Southern China Forests 2015 6 1763 1778 10.3390/f6061763
Du, H. et al. Carbon storage in a Eucalyptus plantation chronosequence in Southern China. Forests 6, 1763–1778 (2015).10.3390/f6061763
8. Darenova E Čater M Effect of spatial scale and harvest on heterogeneity of forest floor CO2 efflux in a sessile oak forest Catena (Amst) 2020 188 104455 10.1016/j.catena.2020.104455
Darenova, E. & Čater, M. Effect of spatial scale and harvest on heterogeneity of forest floor CO2 efflux in a sessile oak forest. Catena (Amst) 188, 104455 (2020).10.1016/j.catena.2020.104455
9. Drake JE Three years of soil respiration in a mature eucalypt woodland exposed to atmospheric CO2 enrichment Biogeochemistry 2018 139 85 101 10.1007/s10533-018-0457-7
Drake, J. E. et al. Three years of soil respiration in a mature eucalypt woodland exposed to atmospheric CO2 enrichment. Biogeochemistry 139, 85–101 (2018).10.1007/s10533-018-0457-7
10. Gutiérrez del Arroyo O Wood TE Large seasonal variation of soil respiration in a secondary tropical moist forest in Puerto Rico Ecol. Evol. 2021 11 263 272 10.1002/ece3.7021 33437428
Gutiérrez del Arroyo, O. & Wood, T. E. Large seasonal variation of soil respiration in a secondary tropical moist forest in Puerto Rico. Ecol. Evol. 11, 263–272 (2021).33437428 10.1002/ece3.7021
11. Ma M Soil respiration of four forests along elevation gradient in northern subtropical China Ecol. Evol. 2019 9 12846 12857 10.1002/ece3.5762 31788219
Ma, M. et al. Soil respiration of four forests along elevation gradient in northern subtropical China. Ecol. Evol. 9, 12846–12857 (2019).31788219 10.1002/ece3.5762
12. Zeng X Song Y Zeng C Zhang W He S Partitioning soil respiration in two typical forests in semi-arid regions, North China Catena (Amst) 2016 147 536 544 10.1016/j.catena.2016.08.009
Zeng, X., Song, Y., Zeng, C., Zhang, W. & He, S. Partitioning soil respiration in two typical forests in semi-arid regions, North China. Catena (Amst) 147, 536–544 (2016).10.1016/j.catena.2016.08.009
13. Shi Q Hua J Creech D Yin Y Biomass estimation and carbon storage of Taxodium Hybrid Zhongshanshan Plantations in the Yangtze River Basin Forests 2022 13 1725 10.3390/f13101725
Shi, Q., Hua, J., Creech, D. & Yin, Y. Biomass estimation and carbon storage of Taxodium Hybrid Zhongshanshan Plantations in the Yangtze River Basin. Forests 13, 1725 (2022).10.3390/f13101725
14. Hatfield, J. L., Gitelson, A. A., Schepers, J. S. & Walthall, C. L. Application of spectral remote sensing for agronomic decisions. Agron. J. 100 (2008).
15. Marcelino do Nascimento D Development of a methodological approach to estimate vegetation biomass using remote sensing in the Brazilian semiarid NE region Remote Sens. Appl. 2022 27 100771
Marcelino do Nascimento, D. et al. Development of a methodological approach to estimate vegetation biomass using remote sensing in the Brazilian semiarid NE region. Remote Sens. Appl. 27, 100771 (2022).
16. Oliveira BR Eucalyptus growth recognition using machine learning methods and spectral variables For. Ecol. Manag. 2021 497 119496 10.1016/j.foreco.2021.119496
Oliveira, B. R. et al. Eucalyptus growth recognition using machine learning methods and spectral variables. For. Ecol. Manag. 497, 119496 (2021).10.1016/j.foreco.2021.119496
17. Zás, L. R. T. Relação entre variáveis dendrométricas de eucaliptus e índices de vegetação obtidos com VANT (2017).
18. Versini A The manipulation of organic residues affects tree growth and heterotrophic CO2 efflux in a tropical Eucalyptus plantation For. Ecol. Manag. 2013 301 79 88 10.1016/j.foreco.2012.07.045
Versini, A. et al. The manipulation of organic residues affects tree growth and heterotrophic CO2 efflux in a tropical Eucalyptus plantation. For. Ecol. Manag. 301, 79–88 (2013).10.1016/j.foreco.2012.07.045
19. Villacrés, J., Fuentes, A., Reszka, P. & Cheein, F. A. Retrieval of vegetation indices related to leaf water content from a single index: A case study of Eucalyptus globulus (Labill.) and Pinus radiata (D. Don.). Plants 10, 697 (2021).
20. Hirigoyen A A machine learning approach to model leaf area index in Eucalyptus plantations using high-resolution satellite imagery and airborne laser scanner data Ann. For. Res. 2022 64 165 183 10.15287/afr.2021.2073
Hirigoyen, A. et al. A machine learning approach to model leaf area index in Eucalyptus plantations using high-resolution satellite imagery and airborne laser scanner data. Ann. For. Res. 64, 165–183 (2022).10.15287/afr.2021.2073
21. Ponzoni, F. J., Shimabukuro, Y. & Kuplich, T. M. Sensoriamento remoto no estudo da vegetação (p. 127). São José dos Campos (2007).
22. Andrade, V. M. S. de. Relações de dependência entre fluxos de CO2 e variáveis meteorológicas na cultura de milho (Zea mays L.) (2009).
23. Susaeta A Adams DC Gonzalez-Benecke C Economic vulnerability of southern US slash pine forests to climate change J. For. Econ. 2017 28 18 32
Susaeta, A., Adams, D. C. & Gonzalez-Benecke, C. Economic vulnerability of southern US slash pine forests to climate change. J. For. Econ. 28, 18–32 (2017).
24. Pereira Ribeiro Teodoro L Eucalyptus species discrimination using hyperspectral sensor data and machine learning Forests 2023 15 39 10.3390/f15010039
Pereira Ribeiro Teodoro, L. et al. Eucalyptus species discrimination using hyperspectral sensor data and machine learning. Forests 15, 39 (2023).10.3390/f15010039
25. Xu N Cai D Zhao X Constructing a eucalyptus identification model based on NDVI time-series remote sensing images Catena (Amst) 2024 238 107846 10.1016/j.catena.2024.107846
Xu, N., Cai, D. & Zhao, X. Constructing a eucalyptus identification model based on NDVI time-series remote sensing images. Catena (Amst) 238, 107846 (2024).10.1016/j.catena.2024.107846
26. Bond-Lamberty B Wang C Gower ST A global relationship between the heterotrophic and autotrophic components of soil respiration? Glob. Chang. Biol. 2004 10 1756 1766 10.1111/j.1365-2486.2004.00816.x
Bond-Lamberty, B., Wang, C. & Gower, S. T. A global relationship between the heterotrophic and autotrophic components of soil respiration?. Glob. Chang. Biol. 10, 1756–1766 (2004).10.1111/j.1365-2486.2004.00816.x
27. Sun Z The effect of nitrogen addition on soil respiration from a nitrogen-limited forest soil Agric. For. Meteorol. 2014 197 103 110 10.1016/j.agrformet.2014.06.010
Sun, Z. et al. The effect of nitrogen addition on soil respiration from a nitrogen-limited forest soil. Agric. For. Meteorol. 197, 103–110 (2014).10.1016/j.agrformet.2014.06.010
28. Ren F Contrasting effects of nitrogen and phosphorus addition on soil respiration in an alpine grassland on the Qinghai-Tibetan Plateau Sci. Rep. 2016 6 34786 10.1038/srep34786 27721415
Ren, F. et al. Contrasting effects of nitrogen and phosphorus addition on soil respiration in an alpine grassland on the Qinghai-Tibetan Plateau. Sci. Rep. 6, 34786 (2016).27721415 10.1038/srep34786
29. Wang X Soil respiration under climate warming: Differential response of heterotrophic and autotrophic respiration Glob. Chang. Biol. 2014 20 3229 3237 10.1111/gcb.12620 24771521
Wang, X. et al. Soil respiration under climate warming: Differential response of heterotrophic and autotrophic respiration. Glob. Chang. Biol. 20, 3229–3237 (2014).24771521 10.1111/gcb.12620
30. Fernández-Martínez M Nutrient availability as the key regulator of global forest carbon balance Nat. Clim. Chang. 2014 4 471 476 10.1038/nclimate2177
Fernández-Martínez, M. et al. Nutrient availability as the key regulator of global forest carbon balance. Nat. Clim. Chang. 4, 471–476 (2014).10.1038/nclimate2177
31. Peñuelas J Human-induced nitrogen–phosphorus imbalances alter natural and managed ecosystems across the globe Nat. Commun. 2013 4 2934 10.1038/ncomms3934 24343268
Peñuelas, J. et al. Human-induced nitrogen–phosphorus imbalances alter natural and managed ecosystems across the globe. Nat. Commun. 4, 2934 (2013).24343268 10.1038/ncomms3934
32. Janssens IA Reduction of forest soil respiration in response to nitrogen deposition Nat. Geosci. 2010 3 315 322 10.1038/ngeo844
Janssens, I. A. et al. Reduction of forest soil respiration in response to nitrogen deposition. Nat. Geosci. 3, 315–322 (2010).10.1038/ngeo844
33. Mo J Nitrogen addition reduces soil respiration in a mature tropical forest in southern China Glob. Chang. Biol. 2008 14 403 412 10.1111/j.1365-2486.2007.01503.x
Mo, J. et al. Nitrogen addition reduces soil respiration in a mature tropical forest in southern China. Glob. Chang. Biol. 14, 403–412 (2008).10.1111/j.1365-2486.2007.01503.x
34. Li D Zhou X Wu L Zhou J Luo Y Contrasting responses of heterotrophic and autotrophic respiration to experimental warming in a winter annual-dominated prairie Glob. Chang. Biol. 2013 10.1111/gcb.12273 24038796
Li, D., Zhou, X., Wu, L., Zhou, J. & Luo, Y. Contrasting responses of heterotrophic and autotrophic respiration to experimental warming in a winter annual-dominated prairie. Glob. Chang. Biol.10.1111/gcb.12273 (2013).24038796 10.1111/gcb.12273
35. Peng Y Nonlinear response of soil respiration to increasing nitrogen additions in a Tibetan alpine steppe Environ. Res. Lett. 2017 12 024018 10.1088/1748-9326/aa5ba6
Peng, Y. et al. Nonlinear response of soil respiration to increasing nitrogen additions in a Tibetan alpine steppe. Environ. Res. Lett. 12, 024018 (2017).10.1088/1748-9326/aa5ba6
36. Silva-Olaya AM Cerri CEP La Scala Jr N Dias CTS Cerri CC Carbon dioxide emissions under different soil tillage systems in mechanically harvested sugarcane Environ. Res. Lett. 2013 8 15014 10.1088/1748-9326/8/1/015014
Silva-Olaya, A. M., Cerri, C. E. P., La Scala Jr, N., Dias, C. T. S. & Cerri, C. C. Carbon dioxide emissions under different soil tillage systems in mechanically harvested sugarcane. Environ. Res. Lett. 8, 15014 (2013).10.1088/1748-9326/8/1/015014
37. Wang J Liu Q-Q Chen R-R Liu W-Z Sainju UM Soil carbon dioxide emissions in response to precipitation frequency in the Loess Plateau, China Appl. Soil Ecol. 2015 96 288 295 10.1016/j.apsoil.2015.08.026
Wang, J., Liu, Q.-Q., Chen, R.-R., Liu, W.-Z. & Sainju, U. M. Soil carbon dioxide emissions in response to precipitation frequency in the Loess Plateau, China. Appl. Soil Ecol. 96, 288–295 (2015).10.1016/j.apsoil.2015.08.026
38. Moitinho MR Effects of burned and unburned sugarcane harvesting systems on soil CO2 emission and soil physical, chemical, and microbiological attributes Catena (Amst) 2021 196 104903 10.1016/j.catena.2020.104903
Moitinho, M. R. et al. Effects of burned and unburned sugarcane harvesting systems on soil CO2 emission and soil physical, chemical, and microbiological attributes. Catena (Amst) 196, 104903 (2021).10.1016/j.catena.2020.104903
39. Ding W Yu H Cai Z Han F Xu Z Responses of soil respiration to N fertilization in a loamy soil under maize cultivation Geoderma 2010 155 381 389 10.1016/j.geoderma.2009.12.023
Ding, W., Yu, H., Cai, Z., Han, F. & Xu, Z. Responses of soil respiration to N fertilization in a loamy soil under maize cultivation. Geoderma 155, 381–389 (2010).10.1016/j.geoderma.2009.12.023
40. Almeida RF Ratio of CO2 and O2 as index for categorising soil biological activity in sugarcane areas under contrasting straw management regimes Soil Res. 2018 56 373 381 10.1071/SR16344
Almeida, R. F. et al. Ratio of CO2 and O2 as index for categorising soil biological activity in sugarcane areas under contrasting straw management regimes. Soil Res. 56, 373–381 (2018).10.1071/SR16344
41. Dossou-Yovo ER Reducing soil CO2 emission and improving upland rice yield with no-tillage, straw mulch and nitrogen fertilization in northern Benin Soil Tillage Res. 2016 156 44 53 10.1016/j.still.2015.10.001
Dossou-Yovo, E. R. et al. Reducing soil CO2 emission and improving upland rice yield with no-tillage, straw mulch and nitrogen fertilization in northern Benin. Soil Tillage Res. 156, 44–53 (2016).10.1016/j.still.2015.10.001
42. Pires LF Soil structure changes induced by tillage systems Soil Tillage Res. 2017 165 66 79 10.1016/j.still.2016.07.010
Pires, L. F. et al. Soil structure changes induced by tillage systems. Soil Tillage Res. 165, 66–79 (2017).10.1016/j.still.2016.07.010
43. González-Ubierna S Lai R Modelling the effects of climate factors on soil respiration across Mediterranean ecosystems J. Arid Environ. 2019 165 46 54 10.1016/j.jaridenv.2019.02.008
González-Ubierna, S. & Lai, R. Modelling the effects of climate factors on soil respiration across Mediterranean ecosystems. J. Arid Environ. 165, 46–54 (2019).10.1016/j.jaridenv.2019.02.008
44. Bashar, M. & Zhou, Y. li Jun, Effect of soil strength and soil physical properties on performance of tillage machines. J. Earth Sci. Eng. 251–255 (2015).
45. Lima, M. C., Amarante, L. do, Mariot, M. P. & Serpa, R. Crescimento e produção de pigmentos fotossintéticos em Achillea millefolium L. cultivada sob diferentes níveis de sombreamento e doses de nitrogênio. Ciência Rural 41, 45–50 (2011).
46. da Silva Junior CA Is it possible to detect boron deficiency in eucalyptus using hyper and multispectral sensors? Infrared Phys. Technol. 2021 116 103810 10.1016/j.infrared.2021.103810
da Silva Junior, C. A. et al. Is it possible to detect boron deficiency in eucalyptus using hyper and multispectral sensors?. Infrared Phys. Technol. 116, 103810 (2021).10.1016/j.infrared.2021.103810
47. Silva, M. de A., Santos, C. M. dos, Vitorino, H. dos S. & Rhein, L. Pigmentos fotossintéticos e índice spad como descritores de intensidade do estresse por deficiência hídrica em cana-de-açúcar. Biosci. J. 173–181 (2014).
48. Bunce JA Carbon dioxide effects on stomatal responses to the environment and water use by crops under field conditions Oecologia 2004 140 1 10 10.1007/s00442-003-1401-6 14557864
Bunce, J. A. Carbon dioxide effects on stomatal responses to the environment and water use by crops under field conditions. Oecologia 140, 1–10 (2004).14557864 10.1007/s00442-003-1401-6
49. Fan X Carbon dioxide fertilization effect on plant growth under soil water stress associates with changes in stomatal traits, leaf photosynthesis, and foliar nitrogen of bell pepper (Capsicum annuum L.) Environ. Exp. Bot. 2020 179 104203 10.1016/j.envexpbot.2020.104203
Fan, X. et al. Carbon dioxide fertilization effect on plant growth under soil water stress associates with changes in stomatal traits, leaf photosynthesis, and foliar nitrogen of bell pepper (Capsicum annuum L.). Environ. Exp. Bot. 179, 104203 (2020).10.1016/j.envexpbot.2020.104203
50. Baret F Guyot G Major DJ Crop biomass evaluation using radiometric measurements Photogrammetria 1989 43 241 256 10.1016/0031-8663(89)90001-X
Baret, F., Guyot, G. & Major, D. J. Crop biomass evaluation using radiometric measurements. Photogrammetria 43, 241–256 (1989).10.1016/0031-8663(89)90001-X
51. Ponzoni FJ Disperati AA Comportamento espectral da vegetação Sensoriamento remoto: reflectância de alvos naturais 2001 1 157 199
Ponzoni, F. J. & Disperati, A. A. Comportamento espectral da vegetação. Sensoriamento remoto: reflectância de alvos naturais 1, 157–199 (2001).
52. Filimon RV Rotaru L Filimon RM Quantitative investigation of leaf photosynthetic pigments during annual biological cycle of Vitis vinifera L. table grape cultivars S. Afr. J. Enol. Vitic. 2016 37 1 14
Filimon, R. V., Rotaru, L. & Filimon, R. M. Quantitative investigation of leaf photosynthetic pigments during annual biological cycle of Vitis vinifera L. table grape cultivars. S. Afr. J. Enol. Vitic. 37, 1–14 (2016).
53. Song G Wang Q Developing hyperspectral indices for assessing seasonal variations in the ratio of chlorophyll to carotenoid in deciduous forests Remote Sens. (Basel) 2022 14 1324 10.3390/rs14061324
Song, G. & Wang, Q. Developing hyperspectral indices for assessing seasonal variations in the ratio of chlorophyll to carotenoid in deciduous forests. Remote Sens. (Basel) 14, 1324 (2022).10.3390/rs14061324
54. Huete AR A soil-adjusted vegetation index (SAVI) Remote Sens. Environ. 1988 25 295 309 10.1016/0034-4257(88)90106-X
Huete, A. R. A soil-adjusted vegetation index (SAVI). Remote Sens. Environ. 25, 295–309 (1988).10.1016/0034-4257(88)90106-X
55. Rondeaux G Steven M Baret F Optimization of soil-adjusted vegetation indices Remote Sens. Environ. 1996 55 95 107 10.1016/0034-4257(95)00186-7
Rondeaux, G., Steven, M. & Baret, F. Optimization of soil-adjusted vegetation indices. Remote Sens. Environ. 55, 95–107 (1996).10.1016/0034-4257(95)00186-7
56. Maas EDvL Lal RA A case study of the RothC soil carbon model with potential evapotranspiration and remote sensing model inputs Remote Sens. Appl. 2023 29 100876
Maas, EDv. L. & Lal, R. A. A case study of the RothC soil carbon model with potential evapotranspiration and remote sensing model inputs. Remote Sens. Appl. 29, 100876 (2023).
57. Primavesi, A. Manejo Ecológico Do Solo: A Agricultura Em Regiões Tropicais (NBL Editora, 2002).
58. Wang W Impact of straw management on seasonal soil carbon dioxide emissions, soil water content, and temperature in a semi-arid region of China Sci. Total Environ. 2019 652 471 482 10.1016/j.scitotenv.2018.10.207 30368177
Wang, W. et al. Impact of straw management on seasonal soil carbon dioxide emissions, soil water content, and temperature in a semi-arid region of China. Sci. Total Environ. 652, 471–482 (2019).30368177 10.1016/j.scitotenv.2018.10.207
59. Hou R Effects of biochar and straw on greenhouse gas emission and its response mechanism in seasonally frozen farmland ecosystems Catena (Amst) 2020 194 104735 10.1016/j.catena.2020.104735
Hou, R. et al. Effects of biochar and straw on greenhouse gas emission and its response mechanism in seasonally frozen farmland ecosystems. Catena (Amst) 194, 104735 (2020).10.1016/j.catena.2020.104735
60. Guo Y Evaluating the contributions of earthworms to soil organic carbon decomposition under different tillage practices combined with straw additions Ecol. Indic. 2019 105 516 524 10.1016/j.ecolind.2018.04.046
Guo, Y. et al. Evaluating the contributions of earthworms to soil organic carbon decomposition under different tillage practices combined with straw additions. Ecol. Indic. 105, 516–524 (2019).10.1016/j.ecolind.2018.04.046
61. Li Z Trade-offs between high yields and soil CO2 emissions in semi-humid maize cropland in northern China Soil Tillage Res. 2022 221 105412 10.1016/j.still.2022.105412
Li, Z. et al. Trade-offs between high yields and soil CO2 emissions in semi-humid maize cropland in northern China. Soil Tillage Res. 221, 105412 (2022).10.1016/j.still.2022.105412
62. Alvares CA Stape JL Sentelhas PC De Moraes Gonçalves JL Sparovek G Köppen’s climate classification map for Brazil Meteorol. Z. 2013 22 711 728 10.1127/0941-2948/2013/0507
Alvares, C. A., Stape, J. L., Sentelhas, P. C., De Moraes Gonçalves, J. L. & Sparovek, G. Köppen’s climate classification map for Brazil. Meteorol. Z. 22, 711–728 (2013).10.1127/0941-2948/2013/0507
63. Cunha, F. F. da, Magalhães, F. F. & Castro, M. A. de. MÉTODOS PARA ESTIMATIVA DA EVAPOTRANSPIRAÇÃO DE REFERÊNCIA PARA CHAPADÃO DO SUL—MS. Revista Engenharia na Agricultura—REVENG 21, 159–172 (2013).
64. Santos, H. G. et al. O novo mapa de solos do Brasil: Legenda atualizada (2011).
65. Sá MMF Fluxes of CO2, CH4, and N2O in tundra-covered and Nothofagus forest soils in the Argentinian Patagonia Sci. Total Environ. 2019 659 401 409 10.1016/j.scitotenv.2018.12.328 31096371
Sá, M. M. F. et al. Fluxes of CO2, CH4, and N2O in tundra-covered and Nothofagus forest soils in the Argentinian Patagonia. Sci. Total Environ. 659, 401–409 (2019).31096371 10.1016/j.scitotenv.2018.12.328
66. QGIS Development Team. QGIS Geographic Information System.
67. R Core Team. R: A Language and Environment for Statistical Computing. Preprint at (2022).
68. Della-Silva, J. L. et al. CO2Flux model assessment and comparison between an airborne hyperspectral sensor and orbital multispectral imagery in Southern Amazonia. Sustainability (Switzerland) 14 (2022).
69. Rahman AF Gamon JA Fuentes DA Roberts DA Prentiss D Modeling spatially distributed ecosystem flux of boreal forest using hyperspectral indices from AVIRIS imagery J. Geophys. Res. Atmos. 2001 106 33579 33591 10.1029/2001JD900157
Rahman, A. F., Gamon, J. A., Fuentes, D. A., Roberts, D. A. & Prentiss, D. Modeling spatially distributed ecosystem flux of boreal forest using hyperspectral indices from AVIRIS imagery. J. Geophys. Res. Atmos. 106, 33579–33591 (2001).10.1029/2001JD900157
70. SANTOS, C. V. B. Modelagem espectral para determinação de fluxo de CO2 em áreas de Caatinga preservada e em regeneração (2017).
71. Wei, T. & Simko, V. R package ‘corrplot’: Visualization of a Correlation Matrix. Preprint at (2021).
72. Friendly, M. & Fox, J. candisc: Visualizing Generalized Canonical Discriminant and Canonical Correlation Analysis. Preprint at (2021).
