Next Article in Journal
Nanotechnology in Food and Plant Science: Challenges and Future Prospects
Previous Article in Journal
Analysis of Virus-Derived siRNAs in Strawberry Plants Co-Infected with Multiple Viruses and Their Genotypes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Monitoring Energy Balance, Turbulent Flux Partitioning, Evapotranspiration and Biophysical Parameters of Nopalea cochenillifera (Cactaceae) in the Brazilian Semi-Arid Environment

by
Alexandre Maniçoba da Rosa Ferraz Jardim
1,2,*,
José Edson Florentino de Morais
3,
Luciana Sandra Bastos de Souza
3,
Carlos André Alves de Souza
3,
George do Nascimento Araújo Júnior
1,
Cléber Pereira Alves
1,
Gabriel Ítalo Novaes da Silva
1,
Renan Matheus Cordeiro Leite
3,
Magna Soelma Beserra de Moura
4,
João L. M. P. de Lima
5 and
Thieres George Freire da Silva
1,3
1
Department of Agricultural Engineering, Federal Rural University of Pernambuco, Dom Manoel de Medeiros Avenue, s/n, Dois Irmãos, Recife 52171-900, Pernambuco, Brazil
2
Department of Biodiversity, Institute of Bioscience, São Paulo State University—UNESP, Av. 24A, 1515, Rio Claro 13506-900, São Paulo, Brazil
3
Academic Unit of Serra Talhada, Federal Rural University of Pernambuco, Gregório Ferraz Nogueira Avenue, s/n, Serra Talhada 56909-535, Pernambuco, Brazil
4
Brazilian Agricultural Research Corporation, Embrapa Semiarid, Petrolina 56302-970, Pernambuco, Brazil
5
MARE—Marine and Environmental Sciences Centre, ARNET—Aquatic Research Network, Department of Civil Engineering, Faculty of Sciences and Technology, University of Coimbra, 3030-788 Coimbra, Portugal
*
Author to whom correspondence should be addressed.
Plants 2023, 12(13), 2562; https://doi.org/10.3390/plants12132562
Submission received: 24 May 2023 / Revised: 20 June 2023 / Accepted: 3 July 2023 / Published: 6 July 2023
(This article belongs to the Section Plant Response to Abiotic Stress and Climate Change)

Abstract

:
The in-situ quantification of turbulent flux and evapotranspiration (ET) is necessary to monitor crop performance in stressful environments. Although cacti can withstand stressful conditions, plant responses and plant–environment interactions remain unclear. Hence, the objective of our study was to investigate the interannual and seasonal behaviour of components of the surface energy balance, environmental conditions, morphophysiological parameters, biomass yield and water relations in a crop of Nopalea cochenillifera in the semi-arid region of Brazil. The data were collected from a micrometeorological tower between 2015 and 2017. The results demonstrate that net radiation was significantly higher during the wet season. Latent heat flux was not significant between the wet season and dry season. During the dry-wet transition season in particular, sensible heat flux was higher than during the other seasons. We observed a large decline in soil heat flux during the wet season. There was no difference in ET during the wet or dry seasons; however, there was a 40% reduction during the dry-wet transition. The wet seasons and wet-dry transition showed the lowest Evaporative Stress Index. The plants showed high cladode water content and biomass during the evaluation period. In conclusion, these findings indicate high rates of growth, high biomass and a high cladode water content and explain the response of the cactus regarding energy partitioning and ET.

1. Introduction

The cactus Nopalea cochenillifera (L.) Salm-Dyck is an important plant used in animal diets in semi-arid environments, particularly in the driest areas of north-eastern Brazil and around the world [1,2,3]. This species of the Cactaceae family is known worldwide for its tolerance to abiotic stress and is able to survive in places with low annual rainfall (250–450 mm) [1,2,3]. Plants in this family that have crassulacean acid metabolism (CAM) are more resistant and have greater water use efficiency (WUE) compared to plants with a C3 or C4 photosynthetic pathway [4,5,6,7]. For instance, N. cochenillifera, besides the CAM pathway for CO2 uptake, has vegetative structures called cladodes (i.e., modified succulent stems). These structures help with drought tolerance and provide the animals with a high concentration of non-fibrous carbohydrates (~54.3%), water (~900 g kg−1 dry matter), minerals (especially calcium, potassium and phosphorus) and high total nutrient digestibility. On the other hand, cladodes have low levels of dry matter, crude protein (~5.4%), neutral detergent fibre (~24.8%), and acid detergent fibre [8,9,10].
Due to their high productive performance, several studies have been carried out with cacti (Opuntia spp. and Nopalea spp.) in a semi-arid environment, adopting management practices or evaluating the water deficit on the growth and productivity of the crop (e.g., [1,2,7,11,12]). However, few studies have been developed to evaluate the plant–environment interaction using the energy balance method, some with Opuntia spp. [13,14] and some with other species and ecosystems that include Cactaceae [15,16,17,18]. Consoli et al. [13] evaluated ecophysiological variables and energy balance in an irrigated crop of O. ficus-indica under a semi-arid climate in Sicily (Italy). They reported the importance of applying a method that detected changes in plant transpiration on a daily scale under wet and dry conditions. Pierini et al. [19] reported that the heat flux is high in vegetation that includes species of cactus (O. spinisior and O. engelmannii) in Tucson in the Sonoran Desert (Arizona).
In terrestrial ecosystems, energy and water exchange plays a fundamental role in hydrological control, climate phenomena and species survival [20,21,22]. To quantify the land-atmosphere fluxes, and understand this partition, the surface energy balance (SEB) predicts variations in turbulent fluxes and evapotranspiration (ET) from interaction of the soil–vegetation–atmosphere system [23,24,25]. In addition to being considered a fairly robust and valid method under semi-arid conditions, it can be applied to different vegetated surfaces (e.g., areas of forest and agricultural crops) and areas with small footprints (i.e., different fetch-to-height ratios) [26,27,28]. Among the methods used to determine turbulent flux, the Bowen ratio indirectly quantifies the latent heat flux (LE) using net radiation (Rn), soil heat flux (G) and the air temperature and humidity gradients [29,30]. In addition to being a low-cost method compared to the eddy covariance and weighing lysimeter methods, the Bowen ratio method uses simple sensors and is preferred in ET studies in agricultural ecosystems [28,31,32].
On a global scale, part of the water from ecosystems is returned to the atmosphere by the ET [33]. Basically, two processes are involved in ET: evaporation—the phenomenon concerning the change in the liquid phase from water to vapour—and transpiration—the process of transferring water via the plant structures (e.g., stomata) to the atmosphere in the form of water vapour. The sum of these two physical phenomena, i.e., evapotranspiration, consumes more than half of the solar energy, with approximately 60% of the precipitation being transferred to the atmosphere by the LE [26,34,35].
ET is perhaps one of the most studied variables in arid and semi-arid environments, helping to explain the response of the vegetation to water availability. However, in dry landscapes, the available energy is high, and when the ET decreases due to water availability, there is an increase in the sensible heat flux (H) [24,26,33,35]. Within this context, these factors can trigger conditions of stress for some plant species and, when added to the excessive heat flux and high energy, may become more hostile.
In the field, hydrological conditions (i.e., wet and dry seasons) can cause changes in the energy fluxes and ET on different scales. In addition, the vegetation can modify the release of energy in response to its growth and phenology [26,36,37]. In this respect, using such information as vegetation indices (e.g., cladode area index, and Normalised Difference Vegetation Index) and growth indices (e.g., growth rate and phenophases) can be fundamental in understanding plant feedback to growing conditions.
Finally, although there are studies on monitoring mass and energy fluxes in dry forests and grasslands in a semi-arid environment [23,38,39,40,41], as far as we know, there are no studies that report information on growth stages, ecophysiology or surface energy balance in N. cochenillifera. Therefore, we hypothesised that (i) even under adverse stress conditions, the performance of the cactus for energy flux and evapotranspiration remains good; and (ii) its growth, water relations and water and nutrient use efficiency are maintained throughout each season. Furthermore, we believe that our results will provide new insights into CAM plants and their exchange of energy with the ecosystem. The objective of this study was to evaluate the interannual and seasonal behaviour of components of the surface energy balance, environmental conditions, morphophysiological parameters and water relations during the wet and dry seasons in a crop of N. cochenillifera in the semi-arid region of Brazil (municipality of Mirandiba).

2. Results

2.1. Surface Energy Balance

Figure 1 shows daytime values for the hourly variation in the energy balance fluxes during the experimental period, with the components being evaluated during the wet and dry seasons and their transitions (i.e., wet-dry and dry-wet seasons). During the wet season, the lowest daily values for net radiation (Rn) were seen in 2015 (Figure 1a) and the maximum in 2016 (Figure 1e), with mean values of 172.21 and 247.92 W m−2, respectively. Among the energy fluxes, most of the Rn came from the sensible heat flux (H), with the exception of 2015, which showed a mean value of 95.49 W m−2. In 2017, the highest value was seen for H (143.58 W m−2), whereas 2015 and 2016 saw the lowest and highest values for latent heat flux (LE), respectively, with mean values of 47.62 and 145.38 W m−2. During the wet season, the soil heat flux (G) was greater in 2017, showing a mean value of 32.13 W m−2, with a gradual increase throughout the day and maximum values around noon (11:00 to 14:00 local time). On the other hand, we found lower values for G during 2016, with a mean of 23.33 W m−2. As shown in Figure 1b,f (wet-dry transition), Rn was 12% lower than during the wet season. There was a considerable variation in LE between 2015 and 2016, with mean values of 80.88 and 51.32 W m−2, respectively. This seasonal pattern was responsible for smaller turbulent fluxes; furthermore, there was an increase in G (34%) compared to the previous season.
The values for Rn obtained throughout the day during the dry season in 2015 and 2016 (Figure 1c,g) were very similar, with mean values of 191.83 and 194.02 W m−2, respectively. These results for the energy budget were similar to the wet-dry transition season and different from those seen during the wet season. Specifically, G showed a significant increase (42%) compared to the wet season. For the turbulent fluxes, LE was 41% greater than H during the dry season in 2015. On the other hand, H showed contrasting behaviour in 2016 due to atmospheric conditions, being 141% higher than in 2015. Even so, the total turbulent flux was similar for the two years mentioned above, with a mean value of 176.87 W m−2. In general, the turbulent fluxes exhibited opposite trends over time and were the components with the greatest contribution to the energy budget. Figure 1d,h show the daily fluxes during the dry-wet transition season for 2016 and 2017. The result for Rn in 2016 was 3.79% higher than in 2017, and comparing the two transition seasons, the value for Rn during the dry-wet transition season was higher than during the wet-dry transition (4%). There was a large temporal variation in G during the study period. During the dry-wet transition season especially, G was approximately 18% higher than during the other seasons. Despite making a significant contribution to the energy balance, our results for H were approximately 2% (wet season), 6% (wet-dry transition) and 7% (dry season) higher during the dry-wet transition season. Furthermore, for the same season, which includes a dry phase, there was a reduction in LE (26%) compared to the other seasons. Based on this, it is clear that most of the available net surface energy contributed to the sensible heat flux.
On a temporal scale, changes can clearly be seen in the energy budget throughout each year (Figure 2). The mean values for global solar radiation (Rg) decreased significantly from May to July, being 22% lower in relation to the other years (Figure 2a). We also saw a greater range in December 2015 (median 24.38 MJ m−2 day−1, interquartile range 18.47 to 25.92) and January 2016 (median 21 MJ m−2 day−1, interquartile range 13.72 to 22.33). In 2015, 2016 and 2017, the mean values for Rg were 20.45, 20.27 and 20.84 MJ m−2 day−1, respectively, despite an intermonthly variation (median 20.96 MJ m−2 day−1, interquartile range 17.75 to 23.73). Figure 2b shows the variation in Rn throughout the experiment. In April 2015, January and February 2016 and from March to April 2017, the values for Rn were higher, with a mean of 8.82 MJ m−2 day−1. However, the mean Rn for the three years was 6.52 MJ m−2 day−1 (median 6.15, interquartile range 5.27 to 7.21), with a reduction of 26% compared to the months with the higher values shown above. As a result of cloudiness, Rn had the lowest values in June, with a mean of 4.51 MJ m−2 day−1.
Figure 2c,d show the turbulent fluxes in the energy balance, with a mean value for LE of 2.32 MJ m−2 day−1 and maximum and minimum values of 5.41 and 0.70 MJ m−2 day−1, respectively, in January 2016 (median 4.97 MJ m−2 day−1, interquartile range 2.88 to 7.83) and July 2016 (median 0.56 MJ m−2 day−1, interquartile range 0.31 to 0.98). After the occurrence of high LE values from January to April 2016, there was a significant value (2.92 MJ m−2 day−1) in May 2016 due to the rainfall events. Obvious inversion of the turbulent fluxes can be seen from the values for LE in relation to H (Figure 2c,d). Although the values for LE were relatively high, the energy consumption was synchronous with that of H. During 2015–2017, H showed a mean value of 5.39 MJ m−2 day−1. Interestingly, the minimum values of H, with a mean of 2.11 MJ m−2 day−1 (median 2.21 MJ m−2 day−1, interquartile range 1.93 to 2.35), and the maximum values, with a mean of 7.21 MJ m−2 day−1 (median 7.24 MJ m−2 day−1, interquartile range 6.80 to 7.86), occurred in October 2015 and 2016, respectively, and due to the influence of the rainfall, resulted in a difference of 71%. Furthermore, H gradually increased in value from July to September 2015 and July to October 2016, with mean values of 4.70 and 6.37 MJ m−2 day−1 during this period, respectively.
Figure 2e,f show the behaviour of G and of the evapotranspiration (ET). For G, the values were more negative during 2015 (April to July, and December), as well as in January 2016 and April 2017. As shown, the minimum values (median –0.77 MJ m−2 day−1, interquartile range −1.37 to 0.10) and maximum values (median 0.71 MJ m−2 day−1, interquartile range 0.36 to 1.07) ranged, on average, from −0.65 to 0.67 MJ m−2 day−1, respectively. We found high values for ET when G and H were low. Despite the interannual variability, the range of values for ET agreed with the turbulent heat fluxes, ranging from 0.29 to 2.22 mm day−1 (median 0.23 MJ m−2 day−1, interquartile range 0.13 to 0.40 and median 2.03 MJ m−2 day−1, interquartile range 1.18 to 3.21, respectively), with a mean of 0.95 mm day−1 during the period under study.
In Figure 3 and Table 1, we present a detailed seasonal boxplot analysis of the energy budget, heat exchange capacity and water vapour in the cactus during the four seasons (wet, wet-dry, dry, and dry-wet), together with the energy partition ratios. For Rn, the values were significantly higher during the wet season (Tukey’s HSD test, p < 0.05; median 6.61 MJ m−2 day−1, interquartile range 4.99 to 9.37), being on average 14% higher than during the other seasons (Figure 3a). Despite showing similar mean values for Rn, during the dry season, the energy was 7% lower than during the dry-wet transition season. For LE (Figure 3b), there was no significant difference (p > 0.05) between the wet and dry seasons, with mean values of 2.87 and 2.41 MJ m−2 day−1, respectively. On the other hand, there was a 25% reduction (p < 0.05) in LE with the arrival of the wet-dry transition season. Notably, the lowest and highest mean values for LE and H, 1.57 and 5.89 MJ m−2 day−1, respectively, were seen during the dry-wet transition season. These effects were even greater for H (Figure 3c), with significantly higher values (p < 0.05) during the dry-wet transition season (median = 6.16 MJ m−2 day−1, interquartile range 5.20 to 6.84). In general, during the dry-wet transition season, H was 13% greater than during the other seasons, showing that in the cactus under a semi-arid environment, most of the Rn was destined for this process (Figure 3c). Verifying the energy partitions (Table 1), the H/Rn ratio was generally responsible for the highest power consumption (58.5%). The LE/Rn ratio varied between 17% and 30%, with a total mean value of 24% over all the seasons. In addition, the G/Rn ratio reached an average of 17.5%, with the smallest partition (13%) during the wet season (Table 1).
The lowest values for G (Figure 3d) were seen during the wet season (−0.28 to 0.28 MJ m−2 day−1), with a marked reduction of 203% in relation to the other seasons. The sharp decrease in G and the G/Rn ratio during the wet season (Table 1) is clearly due to the high volumes of rainfall and soil moisture during this period. Similar to the behaviour of LE, there was no significant difference in ET values during the wet or dry seasons (Figure 3e). This confirms the hypothesis that ET is maintained during the dry season and decreases during the dry-wet transition season (median 0.46 MJ m−2 day−1, interquartile range 0.26 to 0.86). We saw reductions of 40.31% and 26.68% during the dry-wet transition when comparing the wet and dry seasons and the wet-dry transition season, respectively (p < 0.05).
There was a clear significant seasonal variation in the vapour pressure deficit (VPD) (Figure 3f), with a mean of 1.82 kPa during the dry-wet transition season, being 45% and 97% higher during the dry season (median 1.42 kPa, interquartile range 0.53 to 1.84) and wet season (median 0.86 kPa, interquartile range 0.60 to 1.20), respectively. Furthermore, there was a significant difference (p < 0.05) when comparing the VPD of the transition seasons, where the VPD of the wet-dry season was 38% lower than that of dry-wet season (median 1.12 and 1.90 kPa, interquartile ranges 0.88 to 1.38 and 1.52 to 2.19, respectively).

2.2. NDVI Signatures, Environmental Data and the Impact of Drought

We selected a time series for the Normalised Difference Vegetation Index (NDVI), reference evapotranspiration (ET0), Evaporative Stress Index (ESI), available water fraction (AWF) and rainfall during the experimental period (Figure 4). There was no NDVI saturation during any of the cactus growing seasons (Figure 4a,b). The highest NDVI values were seen during the wet season, with a mean of 0.23. Every year, the mean values of the wet-dry, dry and dry-wet seasons (transitions) were, respectively, 9%, 26% and 39% lower than during the wet season. In addition, it should be noted that the NDVI decreased more significantly (17%) during the change from the dry season to the dry-wet transition season (Figure 4b). Despite oscillations in the NDVI trajectory in the area, the final periods of the cycle show greater consistency, with a rise in the index. Variations in the NDVI are influenced by the canopy cover, cladode morphometry (see Figure 5) and environmental conditions (Figure 4c,e). As shown in Figure 4c,d, the mean values for ET0 ranged from 4.26 (±0.83) to 7.85 (±0.63) mm day−1 (mean ± standard deviation), both variations occurring in 2015. Variations, both intermonthly and throughout the year, can also be seen in the ET0, with lower values between May and July, mainly due to the lower incidence of solar radiation and the lower air temperature (Figure 4c). During the present study, the ET0 was, on average, 6.24 mm day−1, with emphasis on the dry season (6.87 mm day−1) and the dry-wet transition season (6.78 mm day−1). On the other hand, during the wet and wet-dry seasons, we found lower values of 5.79 and 5.33 mm day−1, respectively. Although high ET0 values can be found for approximately 48% of the experimental period, the values remained below 6.24 mm day−1.
The monthly mean ESI ranged from −0.32 to 0.88 (Figure 4f). We saw negative ESI values (Figure 4e,f) due to the concentrated rainfall events, favouring greater soil moisture (Figure 4g,i). The wet season and the wet-dry transition season had the lowest ESI values, with means of 0.41 and 0.56, respectively. On the other hand, 2016 saw an exception in the behaviour of the ESI during the wet-dry transition season, with mean values greater than 0.7. During the dry season, the index was 35% higher than during the other seasons, reaching 78% compared to the wet season. This clearly confirms that the most stressful conditions occur during the dry season, which in turn has a high atmospheric demand. The variations in AWF largely occurred together with those of the above-mentioned variables (Figure 4g); in addition, the greatest mean value (0.77) was seen during the wet season. The monthly variation in AWF in the soil ranged from 0.002 to 0.931, and notably, these variations showed more critical values during the dry season, with means of 0.12 and 0.09 during 2015 and 2016, respectively. However, even with a water scarcity during the wet-dry and dry-wet transition seasons, the AWF was 0.48 and 0.31, respectively, being 289% greater compared to the dry season. Although the four seasons under study differ when compared to each other, there was a similar trend in their behaviour for the same season in each of the different years (Figure 4h).
We found a clear variation in rainfall events during the experimental period (Figure 4i,j). The values for 2015, 2016 and 2017 were 183.31 (±17.53), 477.13 (±49.19) and 382.78 (±61.71) mm year−1, respectively, a total of 1043.22 mm (347.74 ± 150.01 mm year−1). In 2016, the period with the greatest accumulation of rain, more than 50% of the rain was concentrated from January to February (wet season). In addition, it should be noted that in March of the same year, the rainfall events were less homogenous, resulting in higher values for the standard deviation. In general, the months from August to November show an extreme water deficit, with the monthly rainfall varying from 0 to 11.44 mm (Figure 4j).

2.3. Growth Parameters, Phenological Characteristics, and Cutting Time

The curves of the morphophysiological parameters, phenophases and cutting time of the cactus are shown in Figure 5. Analysing the absolute growth rate (AGR), we found that the cactus used more than 53% of the thermal time to reach the maximum AGR (0.0041 Mg ha−1 °Cday−1) at 2500 °Cday. Based on the AGR, the ideal cutting time of the plants was at 3760 °Cday, although the species was harvested at 4700 °Cday (Figure 5a). Notably, the relative growth rate (RGR) was higher during the initial growth period and gradually decreased until it reached minimum values. The maximum RGR value of the cactus was 0.0021 Mg Mg−1 °Cday−1. Again, both growth rates (i.e., AGR and RGR) showed a significant fall in biomass accumulation upon reaching the cutting time. When analysing the net assimilation rate (NAR), our results showed a mean NAR of 17.76 Mg ha−1 °Cday−1, with a variation of 1.60 to 40.22 Mg ha−1 °Cday−1. Although the NAR was higher at 2500 °Cday, there was a 50% reduction in the photosynthetic capacity of the cactus by the end of the cycle (4700 °Cday) compared to the initial growth period (Figure 5b).
Figure 5c,d show the performance of the specific cladode area (SCA) and cladode emission rate. The SCA values were higher at the start and decreased significantly over time, as was also seen with the RGR. The maximum SCA reached by the cactus was 0.00065 ha Mg−1; furthermore, when the plants reached the thermal sum of 1500 °Cday, there was a mean reduction of 0.00002 ha Mg−1 in SCA over time. Our results show that the cactus had three cladode phenophases (Figure 5d). The first-order cladodes reached the highest emission rate during these phenophases, with a mean of 0.0136 units °Cday. There were reductions of 32% and 5% in phenophases two and three, respectively, relative to the first phenophase. Interestingly, we also found a prolongation of phenophase two, with an accumulated thermal time of 2632 °Cday, while the first phenophase (1402 accumulated °Cday) and third phenophase (666 °Cday, with the lowest cumulative value) were shorter. We did not completely analyse the third phenophase, as the experiment had already been harvested before the phenophase ended.

2.4. Water Relations, Bowen Ratio, Biophysical Efficiency and Yield

We found significant variation in the water relations of the cladodes throughout the experimental period (Figure 6). The cladode water content (CWC) had a mean value of 89.03%, varying from 77% to 95%, with interquartile ranges of 76% to 78% and 94% to 95%, respectively (Figure 6a). We also found a more-expressive gradual loss of cladode turgidity between samples four and six, and later, the rehydration of the cladodes. Samples five and six included October to December 2016, with low rainfall and low soil moisture (see Figure 4). Although the lower water availability caused significant dehydration, the cladodes rehydrated at the end of the cycle, with an increase of 20% in the CWC. From the second sample onwards, the plants showed a similar variation in cladode succulence to that of the CWC (Figure 6b). On the other hand, due to the high turgidity and smaller cladode area, succulence was generally greater in the first sample (3.16 g cm−2). The mean value for cladode succulence was 0.92 g cm−2, ranging from 0.07 to 4.66 g cm−2. The magnitude of the changes in water relations, whether increasing or decreasing, paralleled the levels of water availability and cladode development.
Table 2 shows our results for the efficiency of the biophysical parameters, the Bowen ratio and the biomass yield in the cactus at the end of the experimental period. The high radiation use efficiency (RUE) of the cactus can clearly be seen, with a mean value of 3.95 g MJ−1, which favoured the process of photosynthesis, and consequently, significant biomass conversion, with a mean of 12.47 Mg ha−1 dry matter. In addition, a value of 1.75 kg m−3 was found for water use efficiency (WUE). Considering the information provided by the surface energy distribution between the turbulent fluxes, our results point to a mean Bowen ratio (β) of 3.53 (Table 2).
According to our results, the cactus showed greater efficiency in the use of phosphorus and potassium, with mean values of 62.13 and 50.88 mg m−3, respectively (Table 2). Calcium use efficiency was 133% higher than that of magnesium. Clearly, the cactus was highly efficient in using the above nutrients, while sodium use efficiency was low (0.22 mg m−3), which indicates greater selectivity of the absorption channels for K+ in relation to Na+.

2.5. Interrelationships between Environmental Variables and Plant Responses

The interrelationship between the environmental variables and the plant variables can be explained using principal component analysis (PCA) (Figure 7). Our results show that the first two principal components were responsible for 86.82% of the total variance, with a strong relationship between the environmental and plant variables during the growing seasons. Despite presenting an eigenvalue greater than 1, the third principal component was not shown due to the lack of any information relevant to the present study. On the other hand, with eigenvalues also greater than 1, the first principal component (PC1) contributed with 50.03%, while the second principal component (PC2) contributed with 36.79% of the total variance. The PCA revealed the clear separation of the growing seasons along the two dimensions of the principal component (Figure 7a). The wet and dry seasons and the wet-dry transition season presented the greatest contribution to PC1, with scores of −2.46, 0.60 and 1.21, respectively. For PC2, the two transition seasons contributed the most, with scores of −2.32 (dry-wet) and 1.63 (wet-dry). In addition to the higher PC1 scores, there was an obvious difference between the groups of environmental and plant variables.
The variables G, NAR, AGR, ESI, VPD, RGR, ET0 and H were positively correlated with PC1, with loadings varying from 0.59 to 0.03 (Figure 7b). On the other hand, the remaining variables (i.e., CWC, Na+, WUE, NDVI, K+, SCA, Ca2+, Mg2+, LE, ET, Rn, Yield, CS, P, and RUE) showed a negative correlation with PC1, with loadings ranging from −0.21 to −0.60. The dry season and wet-dry transition showed the highest correlation (most-positive) with G (0.59), NAR (0.53), AGR (0.49) and RGR (0.21), with the first three having the greatest PC1 loadings. The results clearly show that for the above-mentioned seasons, the cactus shows lower nutrient and water use efficiency (i.e., most-negative PC scores). In particular, the group of variables showing a high correlation with the atmospheric conditions (i.e., H, ET0, VPD and ESI) did so during the dry-wet transition season. In addition, the wet season showed a positive correlation with the largest grouping of biophysical and plant variables (Figure 7b).
All the variables grouped together during the wet season (e.g., CWC, NDVI, LE, ET, CS, SCA, Yield, Rn and RUE) had the highest (most positive) contribution to the total explained variance for PC2 (36.79%). The variables with the highest loading during the wet season were CWC (0.53) and NDVI (0.45). In addition, during the wet season, the cactus showed greater biomass yield, LE and ET than during the other seasons under study. Finally, the correlation between the variables and the wet season shows that the cactus had a greater SCA, CS and RUE, which explains the greater NDVI and biomass yield loadings.

3. Discussion

3.1. Mean Daytime Patterns, Seasonal Variations in the Energy Fluxes and Evapotranspiration

In semi-arid areas with a vegetated surface, the energy balance undergoes changes in the soil-atmosphere system [33]. The hourly energy partition showed the most-sensitive variations, with changes in the fluxes during the morning and afternoon (Figure 1). On the other hand, on a daily scale, the behaviour was similar but of different magnitudes over the years, as well as during the wet to dry and transition seasons (Figure 2). These more marked variations in turbulent flux and soil heat flux were also seen by Shao et al. [39] and require greater attention when there is excessive heating of the environment since this can compromise photosynthetic efficiency. Evaporation and transpiration can reduce the heat stored in the soil and in the plant; however, with the cactus being a CAM plant, transpiration is almost zero during the day [3,4,6], resulting in a smaller contribution to the LE.
Latent heat flux and sensible heat flux are the main variables in net energy consumption, with H prominent in the energy budget of a surface cultivated with cactus. In some cases, LE is less significant, an indication of environmental water limitations. Most of the time, daytime fluxes show a downward concave shape, with minimum values occurring at dawn and dusk [26,42]. It is possible that the proximity of turbulent fluxes in the morning (before 10:00) during some wet and dry seasons may have been due to the increased atmospheric demand and soil moisture that favour evapotranspiration. CAM plants have the ability to cool the soil overnight [43], and although their metabolism is generally nocturnal, we believe that this thermal reduction lasted until the early hours of the morning [27]. In this case, evaporation from the soil may have resulted in greater change in latent heat flux. Unlike forest species, cacti do not store intercepted water in the canopy but in the parenchymatic tissue. Such plants are able to carry out hydraulic redistribution of the water in the soil, making it available in the surface layers. In this way, the evaporated water may come mostly from the soil [44,45].
Based on the above analyses, we found similarities in the Rn curves (Figure 1). However, it is clear that Rn was significantly higher during the wet season due to the greater absorption of thermal and radiant energy by the humid atmosphere (Figure 3), contributing to the LE [23,40]. Using the surface energy balance with the cactus O. ficus-indica (L.) Mill., Consoli et al. [13] also found a higher Rn (~13 MJ m−2 day−1) in a semi-arid environment. The turbulent fluxes and soil heat flux showed marked variations on an hourly scale and during the seasons under study (Figure 1 and Figure 2). G was clearly more expressive during the afternoon, mainly due to the accumulation of energy throughout the morning, and with the greater soil moisture and radiation, there was an increase in the thermal conductivity of the soil [33]. In addition, the lower values for G may be associated with the biophysical and growth parameters of the plants [46]. To substantiate this further, the cactus showed a high growth rate and high NDVI even when there was a reduction in soil moisture, resulting in reduced radiation input due to the greater soil cover. Although the cactus does not present a denser, more-uniform canopy, the cladodes remain on the plant throughout the cycle, which may cause variations in energy exchange due to the spaces between them. With greater results than those reported in the present study, Flanagan and Flanagan [47] found a mean value of 2 MJ m−2 day−1 for G in an area of saguaro cactus (Carnegiea gigantea). Furthermore, the authors point out that the greatest net energy dissipation was via surface heating and the sensible heat flux. Under semi-arid conditions and vegetation consisting of the cacti O. spinisior and O. engelmannii, Pierini et al. [19] also found higher values for the sensible heat flux compared to the latent heat flux.
Indeed, meteorological conditions have a strong relationship with the variables of the surface energy balance (Figure 2 and Table 1). We saw a similar strong consistency in global solar radiation (Rg), with seasonal trends consistent with the climate in the region [48,49]. Areas of low latitude and dry climate in particular have a lower LE and greater H [50]. This is because the incidence of solar radiation is high, and due to the low water availability, the latent heat flux is lower, resulting in limitations on evapotranspiration. The energy balance partition patterns were consistent with the other environmental variables. For example, when the sensible heat energy was low, the cactus increased evapotranspiration, behaving as an important energy sink (Figure 2). Other studies point to variations in turbulent fluxes in areas of Cactaceae and the effects of changes in the energy balance due to water availability (e.g., [13,14,37,47]).
The vapour pressure deficit (VPD) is closely related to ET [51]. Even under deficit conditions, more specifically during the dry season, the cactus maintained ET on a significant scale. Interestingly, even during extremely critical periods (e.g., during the dry season and dry-wet transition), the cactus can keep water lost to the atmosphere well below that seen in species with C3 and C4 metabolism. For these plants, maintaining transpiration during critical periods has important benefits, such as the removal of excessive heat from the cladodes and the maintenance of net photosynthesis [14]. This supports our hypothesis that the cactus manages to maintain satisfactory ET. Under rainfed conditions, as in the present study, the mean was 0.95 mm day−1 during the experimental period, with a value of 1.18 mm day−1 during the wet season, 0.99 mm day−1 during the dry season, 0.88 mm day−1 during the wet-dry transition and 0.65 mm day−1 during the dry-wet transition (Figure 3). In contrast, studies by Consoli et al. [13], Goldstein et al. [52] and Lima et al. [53], evaluating species of cactus (Opuntia spp.) under irrigated conditions, found values for ET greater than 2 mm day−1. Under rainfed conditions, Han and Felker [54] found a mean daily value of 1.44 mm in O. ellisiana, while in a semi-arid ecosystem with cactus (O. engelmannii) and shrubs, Anderson and Vivoni [24] found approximate ET values ranging from 0 to 3.5 mm day−1.

3.2. Variations in the NDVI, Seasonal Environmental Changes and Soil Moisture

Vegetation and terrestrial ecosystems have intrinsic characteristics regarding environmental conditions and, consequently, different spectral responses. To evaluate plant behaviour, the Normalised Difference Vegetation Index (NDVI) is a very promising method that is widely used in arid and semi-arid ecosystems. This is because information in the near-infrared spectrum has a direct relationship with the photosynthetic rate of the vegetation and helps in estimating biomass accumulation [36]. For the cactus, we found an NDVI with a mean value of 0.19 throughout the experimental period (Figure 4). In addition, there are reports of a low NDVI and low vegetation indices in areas of cacti [55]. Similarly, Silva et al. [56] found vegetation indices ranging from 0.1 to 0.5 in areas cultivated with O. stricta and N. cochenillifera in the Agreste region of Pernambuco, Brazil. The quick response of the NDVI for cacti, even under low water availability, may be related to the roots that remain in the surface layers of the soil absorbing rainwater and water vapour from the air [55], thereby maintaining both growth and photosynthesis, even if at a reduced rate. Furthermore, if the high atmospheric demand is combined with water availability in the soil, plant evapotranspiration will be high. Variations in atmospheric demand are generally caused by water availability and the time of year. With a mean of 6.24 mm day−1 (Figure 4c), the results for reference evapotranspiration (ET0) were consistent with the environmental conditions and the rainfall events that occurred. The ET0 in semi-arid regions of Brazil is high due to the high incidence of radiation and low rainfall volumes and can vary on average from 6 mm day−1 in areas of cactus to more than 8 mm day−1 in a seasonally dry forest (Caatinga) [38,57].
In the case of environments with high atmospheric demand and low soil moisture, the Evaporative Stress Index (ESI) can identify problems in plant performance [58]. Although the ESI presented values close to one during the present study, with this value classified as a stressor for the crop [58], the cactus maintained its active development. The resistance to hostile climate and environmental factors is overcome due to the anatomical and morphophysiological characteristics of the plant [1,14,53,54]. On the other hand, when the periods of drought are more severe and prolonged, the yield, growth and biochemical and physiological parameters of the cactus may be impaired [1,7,59]. In addition, characteristics of the root system, parenchymal tissue and stomatal control allow the cactus to take advantage of the soil moisture and show good performance even when water availability is low [7]. The mean value of the rainfall was 347.74 mm year−1, within the values for plant survival [2,12].

3.3. Allometry, Phenological Phase, and Cutting Time

One characteristic of plants is to present differing adaptations, phenotypic plasticity and behaviour in terms of growth rate. Earlier studies (e.g., [11,60,61,62]) also reported similarities in the growth curves of the cactus; however, compared to genus Opuntia, genus Nopalea does not show very high rates of growth. When evaluating growth rates (i.e., absolute growth rate—AGR and relative growth rate—RGR) in N. cochenillifera under rainfed conditions in a semi-arid environment, Araújo Júnior et al. [11] found a lower AGR (<0.004 Mg ha−1 °Cday−1) and higher RGR (>0.002 Mg Mg−1 °Cday−1) compared to our results (Figure 5). The AGR helps explain the ratio of dry matter accumulation in the crop, while the RGR shows the relative increase in dry matter per unit of time [61]. Generally, a reduction in both of the above rates can clearly be seen throughout the cycle. This is because as plants grow, they consume photoassimilates, causing a natural reduction in these rates [11,60,61,62]. The ideal cutting time favours a more promising crop, with better use of photoassimilates and better biomass accumulation [61]. Other studies point out that the cactus shows variations in harvest time depending on the agronomic management and/or environmental conditions [11,60,61,62].
Our results showed high values for the cactus in terms of photosynthetic capacity, a ratio predicted by the net assimilation rate (NAR). From our findings, the crop presented significant net photosynthesis, with little influence from self-shading of the cladodes due to the smaller specific cladode area (SCA), favouring a high NAR (Figure 5b,c). For Scalisi et al. [7], morphometric variables such as cladode area help us understand the growth behaviour of plants. The number of cladodes is also a fundamental variable for understanding the phenology of the cactus. In the present study, the cactus presented up to the third phase, i.e., the successive emission of third-order cladodes (Figure 5d). This agrees with other studies on cactus subjected to a semi-arid environment [11,61]. These findings also confirm the hypothesis that plants continue to develop, increasing their growth rate under different environmental conditions. On the other hand, it is possible that, despite the plants having reached the third phase, the phases were prolonged due to the low water availability.

3.4. Water Relations and Biophysical Parameters of Cladodes

Anatomically, the cactus has the capacity to store water in various structures, such as parenchyma and hydrenchyma cells [1,7]. Our results showed that both the cladode water content and cladode succulence presented significant variations (Figure 6). During dry periods, samples five and six showed greater water loss; plants from sample six in particular lost more water due to the long period of drought [7]. In addition, when plants lose water to the environment, they reduce their cladode area, resulting in a loss of cell turgor and succulence (Figure 6b) [6,7,44]. Even with a significant loss of cladode water content, succulent plants are able to store water more efficiently than plants with a C3 or C4 photosynthetic pathway, favouring species survival and the metabolic pathways [6,59]. While the cladodes are still young and have a smaller area than mature cladodes, cladode succulence may be greater since the area is reduced and the water content is high, as seen in the first batch of samples (Figure 6b). These results corroborate those of Scalisi et al. [7], who found variations of approximately 45% to 85% in cladode water content. When the number of cladodes increases, their water content may decrease due to competition for water by the cladodes. Other studies point to variations of 60% to 95% in the cladode water content, influenced by soil moisture [63,64]. Such variations in turgor may occur due to the drought tolerance of the plant since the cactus is able to maintain turgor pressure for months [6,52].
When plants are highly efficient in using biophysical resources, the probability of achieving better yields is high [57]. The results shown here generally underline the high efficiency of the cactus in using radiation, water and nutrients (see Table 1). One of the factors to help in radiation use efficiency (RUE) is the canopy architecture, which favours the interception of photosynthetically active radiation (PAR) and the capture of CO2, to be later converted into biomass [6,47,65]. The relationship between RUE, water use efficiency (WUE) and nutrient use efficiency (NUE) is true since the plants showed high dry biomass yield. The study by Han and Felker [54] reported the high WUE of O. ellisiana (1 kg of dry matter per 162 kg of water); however, such results are still in the early stages for N. cochenillifera under rainfed conditions. Our findings of 1.75 kg m−3 are higher than those reported by Mbava et al. [66] for C3 and C4 species, e.g., wheat (1.18 kg m−3), sorghum (1.48 kg m−3), maize (1.47 kg m−3) and cotton (0.22 kg m−3).
As seen in earlier studies with the cactus [67,68,69], we found that the plants showed greater efficiency and accumulation for such nutrients as P, K+, Ca2+, Mg2+ and Na+ (Table 2), in order of use efficiency. Nedjimi [68] and Saraiva et al. [69], respectively evaluating Opuntia and Nopalea, also found greater concentrations of P and K+ in the cladodes. They concluded that cacti have higher concentrations of the above nutrients, favouring the nutritional value for human/animal consumption. In addition, the high concentration of nutrients such as K+ can help improve performance against abiotic factors, such as environmental changes, salinity and heat [1]. The high Bowen ratio was due to the greater heat flux (Table 2), conditions characteristic of a deficient climate [23], which can be damaging to crops. High values can also be seen for the H/LE ratio in an arid ecosystem of cactus [47]. The Bowen ratio can change depending on the meteorological variables, soil conditions and vegetation, in which case, adapted plants tend to suffer less. Despite the semi-arid conditions to which the plants were exposed, the cactus showed a high yield of dry biomass (mean of 12.47 Mg ha−1). In this respect, our data are superior to those found by Jardim et al. [1,57], where the cactus was harvested at a younger age in a semi-arid environment.

3.5. Principal Component Analysis (PCA)

For the purposes of this study, we applied PCA to the variables that were most influenced by the wet and dry seasons and their transitions (Figure 7). The analysis allowed a reduction in the data set, transforming the data into a series of interrelated variables without losing the principal characteristics of the data. Several variables were correlated with the growing seasons, and the first two principal components (PC) explained 86.82% of the total variance, which was distributed over the 23 variables of the PC1 and PC2 coordinates. The formed groups are clearly correlated based on water availability and the environmental conditions. This is because, in addition to the microclimate caused by the rainfall, according to Campos et al. [12], the cactus responds to different levels of soil water availability. Jardim et al. [1] used the PCA method to understand which environmental and plant variables have greater interaction with different cactus genotypes (Opuntia spp. and Nopalea spp.) irrigated with saline water in the semi-arid region of Brazil. Furthermore, according to the authors, this type of analysis helps in understanding how plants behave in the face of environmental stressors.
Previous studies showed the grouping of cladode and plant variables [1,2], where correlated variables of the same sign explain, for example, the greater yield of the plants. It is interesting to see water use efficiency grouped with nutrient use efficiency since these variables were inversely related to the growth rates, which explains their lesser accumulation during the wet-dry transition and dry season. Recently, several studies have identified variations in the energy budget during wet and dry seasons [22,23,33,51]. Plant evapotranspiration may be limited by the available energy, water and evaporative demand; thus, as the plants’ increased evapotranspiration, the heat stress index and evaporative demand were reduced. The information presented here implies a clear effect from the wet and dry seasons on energy balance, turbulent fluxes and plant responses due to the grouping of the variables and the contribution of each component.

4. Materials and Methods

4.1. Location and Information of the Experimental Area

The experiment was conducted in a study area located in the district of Mirandiba, in the state of Pernambuco, Brazil (8°3.73′ S, 38°43.69′ W, altitude 490 m) under the conditions of a water deficit (Figure 8). According to the Köppen classification, the climate in the region is classified as semi-arid, type BSh (i.e., dry and hot with a rainy season in the summer) [70,71]. Rainfall predominates from January to June, with an annual average of 431.8 mm, average air temperature of 25.2 °C, relative humidity of 64.6%, and a high atmospheric demand that can reach 1600 mm year−1. The environment is characterised by high luminosity, irregular spatiotemporal distribution of the rainfall and the occurrence of periodic droughts.
Based on the beginning and end of the rainfall regime, we considered four seasons (wet, wet-dry transition, dry and dry-wet transition) for each year under study from 2015 to 2017. If the sum of rainfall values in the 30 days preceding or following the day under analysis is less than 20 mm, and there have been less than five rainy days, it is considered a dry season. However, if in the same time interval, rainfall values greater than 20 mm are recorded on five or more days, the season is considered wet. A transition season is determined when none of these criteria apply. If this happens after a dry season, it is said to be a dry-wet transition season, and when it occurs after a wet season, it is considered a wet-dry transition season [72,73]. For example, using this approach, the wet season included 21 February 2015 to 23 July 2015, 7 January 2016 to 1 February 2016 and 21 February 2017 to 1 July 2017; the wet-dry transition included 24 July 2015 to 21 August 2015 and 2 February 2016 to 26 July 2016; the dry season included 22 August 2015 to 7 December 2015 and 27 July 2016 to 17 November 2016 and the dry-wet transition included 8 December 2015 to 6 January 2016 and 18 November 2016 to 20 February 2017.
The soil in the experimental area was classified as a Chromic Luvisol [74], comprising 54% sand, 21% silt and 25% clay. Soil samples were collected at depths of 0–0.20 m, with a bulk density of 1.45 g cm−3 (measured by the cutting-ring sampling method), soil organic carbon content of 10.1 g kg−1, field capacity of 0.26 cm3 cm−3, soil pH of 6.1 (measured at a soil to water ratio of 1:2.5 using an ion pH meter), and terrain slope between 3% to 5%.
The experiment was conducted from April 2015 to April 2017 over a total area of 1.8 ha, using the cactus, Nopalea cochenillifera (L.) Salm-Dyck, a species tolerant to Dactylopius opuntiae Cockerell (Hemiptera: Dactylopiidae). The crop was planted in January 2015. After the initial soil preparation (i.e., ploughing, harrowing and furrowing), the cladodes were planted, leaving 50% of their height in the soil. The plants were arranged in single crop rows with the cladodes aligned bilaterally (i.e., parallel to each other), at a spacing of 2.0 × 0.50 m (10,000 plants ha−1). Cropping treatments (i.e., hand weeding and the application of herbicide and insecticide) were carried out whenever necessary to avoid competition with spontaneous plants and promote full growth of the crop. The cactus was grown under rainfed conditions throughout the evaluation cycle.

4.2. Measuring the Meteorological Variables

A 3-metre micrometeorological tower was installed in the centre of the experimental area. Net radiation (Rn) data were measured using a closed-cell thermopile sensor (NR-Lite, Kipp and Zonen, Delft, Netherlands, accuracy = ±10 μV W−1 m−2). The soil heat flux (G) was measured at a depth of 0.05 m from the surface using a heat flux plate model HFT-3 (REBS, Hukseflux, Delft, Netherlands, accuracy = ±5% of reading) buried close to the crop row. Air temperature (Ta) and relative humidity (RH) were determined using two aspiration psychrometers (height of 0.5 and 1.5 m above the surface of the soil). Wind speed and direction were measured using a Wind Sentry model 03002 anemometer (R. M. Young Company, Traverse City, MI, USA, accuracy = ±0.5 m s−1). Rainfall data were quantified using an automatic rain gauge (CS700-L, Hydrological Services Rain Gauge, Liverpool, Australia) installed 3 m above the canopy, and the photosynthetically active radiation (PAR) was measured using a point quantum sensor (LI-190SB, LI-COR Inc., Lincoln, NE, USA, accuracy = ±5% of reading) installed above and a 1 m line-quantum sensor (LI-191R, LI-COR, accuracy = ±5% of reading) below the cactus canopy. In addition, the height of the sensors met the fetch-to-height ratio of 100:1 [75].
In our study, the data were measured every 60 s by a CR10X data logger (Campbell Scientific Inc., North Logan, UT, USA) with a storage interval of 10 min. Measurements were collected continuously both day and night. However, the flux data used in applying the Bowen ratio-energy balance (BREB) method were considered on a daytime scale only. This is because at night, when the temperature/humidity gradients are small, erroneous flux data may be quantified. Each dataset is referred to in terms of local time (GMT-3).

4.3. Surface Energy Balance Method

We used the surface energy balance (SEB) method, which is based on the law of conservation of energy (Equation (1)).
R n     G = LE + H  
where Rn is the net radiation (W m−2), G is the soil heat flux (W m−2), LE is the latent heat flux (W m−2) and H is the sensible heat flux (W m−2).
The turbulent flux components, i.e., sensible heat and latent heat fluxes were determined based on the Bowen ratio (β) (Equation (2)). The Bowen ratio method is widely used to partition energy flux components in relation to the total available energy (RnG). We can therefore estimate this ratio by quantifying the temperature gradient and vapour pressure above the canopy [29]. In the present study, we assume similarity of equality between the turbulent transfer coefficients of the sensible heat (Kh) or of the water vapour (Kw) [76]. Therefore, latent heat flux (LE) and sensible heat flux (H) were estimated by combining the available energy balance and Bowen ratio (Equations (3) and (4), respectively).
  β = H LE = P a   ×   c p λ   ×   ε × T e × K h K w
LE = R n     G 1 + β  
  H = β 1 + β × R n     G
where β is the Bowen ratio (dimensionless), Pa is the atmospheric pressure (kPa), cp is the specific heat capacity of the air (1004.67 J kg−1 °C−1), λ is the latent heat of vaporisation (2.454 MJ kg−1 at 20 °C), ε is the ratio of the molecular weights of the air and water vapour (0.622), ΔT is the difference in air temperature between the two heights (°C), Δe is the difference in vapour pressure between the two heights (kPa), Kh is the eddy diffusivity for heat (m2 s−1) and Kw is the eddy diffusivity for water vapour (m2 s−1). Quality control was carried out on the calculated LE as per Perez et al. [77].
Furthermore, in the present study, to calculate the energy balance closure, we ignored the energy from metabolic activities and heat storage in the plant tissue and in the canopy, as well as horizontal advection. This condition can be applied when the surface is uniform, keeping in mind that the vertical gradient that comprises the meteorological elements is far greater than the horizontal gradient. As such, these terms are ignored in building Equation (1) [30,31].

Data Selection Criteria for the Energy Balance Method

In the present study, for the acceptance and/or rejection criteria of the data collected by the Bowen ratio-energy balance (BREB) method, we used the approach proposed by Perez et al. [77] (Table 3). That said, the presence of abnormal data occurs when the available-energy heat flux (RnG) is very small. In this way, faults may occur in applying the energy balance by the Bowen ratio (β) and later trigger several errors, e.g., when (1) the sensor resolution is inadequate to solve the gradient in ΔT and Δe; (2) stable atmospheric conditions, e.g., at dawn and dusk, return β values close to −1, resulting in evapotranspiration tending to infinity, which is inconsistent; and (3) the conditions change abruptly, causing measurement errors (Table 4) [30,77,78,79].

4.4. Resource Use Efficiency

4.4.1. Water Efficiency

The actual evapotranspiration component (ET, mm day−1) was calculated as described in Equation (5) [80]. To help understand the aboveground dry biomass yield in relation to the volume of water consumed, we calculated the water use efficiency (WUE, kg m−3) of the crop (Equation (6)) [81].
ET = 86 , 400   × LE λ   ×   ρ w
WUE = Y ET  
where 86,400 is the time unit conversion factor (i.e., converting from seconds to days), LE is the latent heat flux measured over 24 h (W m−2), ρw is the density of water (kg m−3), Y is the dry matter yield (kg ha−1) and ET is the total amount of seasonal evapotranspiration (m3 ha−1).

4.4.2. Radiation Use Efficiency

The radiation use efficiency (RUE) of the cactus was calculated to describe the light-absorption dynamics of the crop (Equation (7)). To achieve this, we used the ratio between the total aboveground dry weight of the plant and the radiation intercepted throughout the growing season [82].
RUE = Y I 0   ×   f PAR
where RUE is the radiation use efficiency (g MJ−1), I0 is the amount of daily incident photosynthetically active radiation (PAR) above the canopy (MJ m−2) and fPAR is the fraction of intercepted photosynthetically active radiation.
The intercepted photosynthetically active radiation was estimated monthly using the AccuPAR LP-80 ceptometer (Decagon Devices, Logan, UT, USA), manually calibrated prior to taking the readings. To maintain consistent measurement conditions and avoid the influence of the angle of the sun’s rays, readings were taken under clear skies between 11:00 and 13:00, recording the incident radiation above and below the canopy. Three simultaneous readings were taken below the canopy (with the sensor rod positioned parallel and perpendicular to the crop rows) and one above, on four similar plants in the experimental area. We then calculated fPAR using Equation (8), and the light extinction coefficient (k) based on the Beer–Lambert law (Equation (9)). As the cactus has cladodes instead of leaves, we calculated k considering the cladode area index (CAI).
  f PAR = 1     I t I 0
f PAR = 1   e ( k   ×   CAI )  
where It is the PAR measured at the bottom of the canopy and CAI is the cladode area index, determined according to Pinheiro et al. [83].

4.4.3. Nutrient Use Efficiency

After cutting (aboveground biomass), the cladodes were weighed on an electronic scale to quantify the fresh biomass (g FM plant−1) and then dried in a forced air circulation oven at 55 °C to constant weight, i.e., dry matter per plant (g DM plant−1). The dried cladode samples were then ground using a Model 4 Wiley mill (Thomas Scientific, Swedesboro, NJ, USA) with a 1 mm sieve. The mineral-element concentration in the plant tissue was then determined: phosphorus (P, mg kg−1) using the vanadate-molybdate method with readings by UV-visible spectrophotometry at 430 nm [67]; potassium (K+, mg kg−1) and sodium (Na+, mg kg−1) by flame photometry [84]; and Calcium (Ca2+, mg kg−1) and magnesium (Mg2+, mg kg−1) by atomic absorption spectrophotometry [85]. Finally, nutrient use efficiency (NUE, mg m−3) was calculated as per Equation (10).
NUE = Y   ×   Nu ET
where Nu is the concentration of the nutrient in the analysed sample of plant tissue (mg kg−1). To improve understanding of nutrient uptake by the plants, we adapted Equation (10) [86], making it a function of crop evapotranspiration. This adaptation affords more clarity in explaining the nutrient uptake capacity of the plant from the soil solution, together with the water consumption lost through evapotranspiration.

4.5. Analysing Growth, Phenology, Cutting Time and Yield

Morphometric data were collected monthly, and samples of plant biomass were taken at 60-day intervals. Four plants were measured for each analysis in each evaluation period. Height was evaluated considering the vertical distance from the ground to the apex of the canopy, and width, considering the average of two measurements from the edge of the canopy. The length (CL, cm), width (CW, cm) and perimeter (CP, cm) of the cladodes were measured, and the number of cladodes (NC, units) determined by counting the cladodes in order of appearance on the plant (i.e., first-order, second-order, third-order and so on). The total number of cladodes was determined summing the cladodes by order.
The cladode area (CA) and CAI were determined from the morphometric data as per the equations proposed by Silva et al. [87] and Pinheiro et al. [83] (Equations (11) and (12)). In addition, from the dry mass yield and cladode measurements, we calculated the morphophysiological indices and phenology using a sigmoidal model with three parameters (Equation (13)) and accumulated degree-days (Equation (14)) [61,88].
CA = 1 . 6691   × 1     e 0 . 0243   ×   CP 0 . 0243
CAI = [ n i = 1 ( CA ) / 10 , 000 ( S 1   ×   S 2 )   ]
  y = a 1 + e   x x 0 b
ADD = j = 1 n [ T m a x + T m i n 2     Tb ]
where i is the observation number, n is the total number of observations, 10,000 is the conversion factor from cm2 to m2 and S1 × S2 is the spacing between the rows and plants (1.0 × 0.2 m), respectively.
The following parameters were used in the morphophysiological and phenological analysis: y is the response variable (e.g., cladode dry matter, cladode area index and the number of cladodes); a is the maximum value for the growth rate (i.e., the distance between the two asymptotes); x is the accumulated degree-days; x0 is the number of degree-days necessary for the plant to express 50% of the maximum growth rate (i.e., the inflection point of the curve); b is the number of degree-days necessary to the start of growth; ADD is the accumulated degree-days (°Cday); j is the daily time step; n is the total number of days; Tmax is the daily maximum air temperature (°C); Tmin is the daily minimum air temperature (°C); Tb is the lower base temperature (°C). Tb is the minimum temperature at which cactus cladodes grow (22 °C) [61]. If the daily average temperature drops below 22 °C, the ADD becomes negative. Negative ADD values were set to zero.
The growth indices were then quantified: absolute growth rate (AGR, Mg ha−1 °Cday−1); relative growth rate (RGR, Mg Mg−1 °Cday−1); net assimilation rate (NAR, Mg ha−1 °Cday−1); specific cladode area (SCA, Mg−1 °Cday−1); and cladode emission rate (CER, units °Cday−1) based on earlier studies [61]. The cladode emission rate is a very important variable in phenological analysis, as it is used to obtain the phenological phases of the crop [60,61]. We determined the cutting time as when the absolute growth rate reached 25% of the maximum peak [61].

4.6. Measuring the Plant Water Status

The cladode water content (CWC, %) was calculated from the fresh matter weight (FM) and dry matter weight (DM) of the samples (Equation (15)) and from the cladode succulence (CS, g cm−2) (Equation (16)) [5,89].
CWC = FM     DM FM ×   100  
CS = FM     DM CA

4.7. Indicators of Water and Environmental Stress

In this study, we examined the Normalised Difference Vegetation Index (NDVI) and the Evaporative Stress Index (ESI) for seasonal and interannual variations. The association between the NDVI and ESI allows the health status of plants to be understood in response to the intensity of hydrological drought [90,91]. The NDVI was derived using images from the Landsat-8 Operational Land Imager (OLI) and Landsat-7 Enhanced Thematic Mapper (ETM+), available on the United States Geological Survey (USGS) website at 30 m resolution for the red (0.64–0.67 μm) and near-infrared (0.85–0.88 μm) bands [22,92]. Landsat images acquired from April 2015 to April 2017 were then processed using the Google Earth Engine (GEE) platform. With the GEE Application Programming Interface (API) we implemented code written in JavaScript to quantify the NDVI [22]. In addition, a pixel-quality attribute mask was applied to the images to mask clouds, cloud shadow and water using the CFMask algorithm, with the values expressed as surface reflectance [93,94]. The NDVI can take dimensionless values ranging from −1 to +1, being positive when the crops show photosynthetic activity, and generally referring to bodies of water when negative. On the other hand, low positive NDVI values are characteristic of stressed vegetation or vegetation with a small leaf area.
Although the ESI can be quantified via remote sensing, we used field data to quantify it as per Equation (17) [58]. The ESI generally ranges from 0 to 1 and is linked to the evaporative demand of both the surface and the atmosphere. When the ESI is close to 1, it indicates water stress in the ecosystem; when it approaches 0, it indicates the absence of water stress [41,91]. In addition, the reference evapotranspiration (ET0) was estimated using the FAO-56 Penman-Monteith method and then multiplied by the crop coefficient (Kc) to obtain the potential crop evapotranspiration (PETc) [34].
ESI = 1 ET PET c

4.8. Soil Moisture

We considered the soil water balance method proposed by Thornthwaite–Mather to determine the available water fraction (AWF) over time [95,96,97,98]. This fraction expresses the total available water that a crop can extract without undergoing water stress [95,97]. Daily soil and micrometeorological data from the tower located in the experimental area were used. The available water capacity of the soil and a respective root-system depth of 75 cm and 60 mm were considered as per Aparecido et al. [99] and Almagbile et al. [100].

4.9. Statistical Analysis

In the first step, we applied nonlinear sigmoid functions for the morphophysiological indices and cutting time, adopting the significance of the F-test (p < 0.05) and the coefficient of determination (R2 > 0.85) as the criterion for choosing the model. Data on the NDVI, ET0, ESI, available water fraction and rainfall were submitted to descriptive statistics and expressed as the mean and standard deviation. In the next step, the time-scale data of the energy balance components (Rn, LE, H and G), ET and VPD for the four seasons under study (i.e., wet season, dry season, wet-dry transition and dry-wet transition) were submitted to one-way analysis of variance (ANOVA) by F-test (p < 0.05). The mean values were compared by Tukey’s HSD (honestly significant difference) test for multiple comparisons (p < 0.05). Data on the cladode water content and environmental conditions were evaluated using boxplots showing the median, interquartile range and 1.5 times the interquartile range.
Finally, the interrelationships between the plant parameters and environmental conditions were tested using principal component analysis (PCA). We applied PCA to the mean value of the environmental variables (Rn, LE, H, G, VPD, ESI, NDVI and ET0) and plants variables (ET, NUE, WUE, AGR, RGR, NAR, SCA, CWC, CS, RUE and yield). The response variables were standardised using the z-transform (mean = 0, standard deviation = 1), subtracting from the mean value and then dividing by the standard deviation. This is due to the different magnitudes and units of the variables under study and makes them directly comparable. Significant principal components were selected according to the Kaiser criterion, considering only eigenvalues greater than 1.0 [1,22,56]. The analysis was carried out using the R software [101].

5. Conclusions

In conclusion, our findings demonstrate that the energy budget showed seasonality, with the latent heat flux (LE) and sensible heat flux (H) compromised more strongly during the dry-wet transition season. The values for the H/Rn ratio on a seasonal scale were always higher than those of the LE/Rn ratio during each of the four seasons. In general, during the wet and dry seasons, the cactus vegetation maintained high evapotranspiration (1.08 mm day−1). The dry season showed the lowest net radiation (5.89 MJ m−2 day−1). The soil heat flux (G) was strongly dependent on the rainfall and the end of the wet season; there is no difference in soil heat transfer over the seasons. The seasonal fluctuations in the NDVI and growth rates confirm the influence of the wet and dry seasons. However, the cactus plants maintained their growth even during the most critical periods. When there is moisture in the soil and a low vapour pressure deficit, more energy is consumed for LE, resulting in greater evapotranspiration; while there is a water restriction, more energy is attributed to G and H, also increasing the Bowen ratio and Evaporative Stress Index. We found that the energy distribution is greater and preferably converted into sensible heat. The study also demonstrated how principal component analysis can help and be effectively used to understand the way plants respond to environmental factors, and which variables of the energy balance are correlated with plant growth.
Moreover, our results can be a particularly valuable baseline for studies on cacti in semi-arid ecosystems, where we believe this to be a pioneering study, with application of the surface energy balance and flux partitioning to Nopalea cochenillifera (L.) Salm-Dyck. Finally, these findings may also be useful for decision makers in environmental management, the rehabilitation of degraded lands, and climate change, targeting CAM plants as an alternative way of reducing heat fluxes, and for their tolerance to low water availability. In future studies, it would be interesting to quantify CO2 fluxes in the soil and in the environment in wetlands and drylands cultivated with cactus.

Author Contributions

Conceptualisation, A.M.d.R.F.J.; methodology, A.M.d.R.F.J. and J.E.F.d.M.; software, A.M.d.R.F.J. and T.G.F.d.S.; validation, A.M.d.R.F.J., J.E.F.d.M. and T.G.F.d.S.; formal analysis, C.A.A.d.S., G.d.N.A.J. and C.P.A.; investigation, A.M.d.R.F.J., G.Í.N.d.S. and R.M.C.L.; resources, T.G.F.d.S., M.S.B.d.M. and J.L.M.P.d.L.; data curation, A.M.d.R.F.J., J.E.F.d.M. and T.G.F.d.S.; writing—original draft preparation, A.M.d.R.F.J.; writing—review and editing, J.E.F.d.M. and T.G.F.d.S.; visualisation, J.L.M.P.d.L., G.d.N.A.J. and C.P.A.; supervision, L.S.B.d.S. and M.S.B.d.M.; project administration, T.G.F.d.S. and J.L.M.P.d.L.; funding acquisition, J.L.M.P.d.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financed by National Funds through the Portuguese funding agency, FCT—Fundação para a Ciência e a Tecnologia, I.P., under the projects UIDB/04292/2020, UIDP/04292/2020, granted to Marine and Environmental Sciences Centre—MARE, and LA/P/0069/2020, granted to the Associate Laboratory ARNET.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors wish to thank the Research Support Foundation of the State of Pernambuco (FACEPE-APQ-0215-5.01/10 and FACEPE-APQ-1159-1.07/14), the São Paulo Research Foundation (FAPESP-Process #2023/05323-4) and the National Council for Scientific and Technological Development (CNPq-Process #152251/2018-9) for their financial support. Further thanks go to the Research Support Foundation of the State of Pernambuco and the Coordination for the Improvement of Higher Education Personnel (CAPES-Finance Code 001) for the research and study grants. We are also grateful for financial support from the Portuguese Foundation for Science and Technology (FCT), and the University of Coimbra, Portugal. We also would like to thank the anonymous reviewers for their insightful comments on previous versions of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jardim, A.M.d.R.F.; Santos, H.R.B.; Alves, H.K.M.N.; Ferreira-Silva, S.L.; de Souza, L.S.B.; Araújo Júnior, G.d.N.; Souza, M.d.S.; de Araújo, G.G.L.; de Souza, C.A.A.; da Silva, T.G.F. Genotypic differences relative photochemical activity, inorganic and organic solutes and yield performance in clones of the forage cactus under semi-arid environment. Plant Physiol. Biochem. 2021, 162, 421–430. [Google Scholar] [CrossRef]
  2. Kumar, S.; Palsaniya, D.R.; Kumar, T.K.; Misra, A.K.; Ahmad, S.; Rai, A.K.; Sarker, A.; Louhaichi, M.; Hassan, S.; Liguori, G.; et al. Survival, morphological variability, and performance of Opuntia ficus-indica in a semi-arid region of India. Arch. Agron. Soil Sci. 2022, 69, 708–725. [Google Scholar] [CrossRef]
  3. Winter, K.; Smith, J.A.C. CAM photosynthesis: The acid test. N. Phytol. 2022, 233, 599–609. [Google Scholar] [CrossRef] [PubMed]
  4. Hartzell, S.; Bartlett, M.S.; Inglese, P.; Consoli, S.; Yin, J.; Porporato, A. Modelling nonlinear dynamics of Crassulacean acid metabolism productivity and water use for global predictions. Plant. Cell Environ. 2021, 44, 34–48. [Google Scholar] [CrossRef]
  5. Ho, C.L.; Chiang, J.M.; Lin, T.C.; Martin, C.E. First report of C4/CAM-cycling photosynthetic pathway in a succulent grass, Spinifex littoreus (Brum. f.) Merr., in coastal regions of Taiwan. Flora 2019, 254, 194–202. [Google Scholar] [CrossRef]
  6. Nobel, P.S. Achievable productivities of certain CAM plants: Basis for high values compared with C3 and C4 plants. New Phytol. 1991, 119, 183–205. [Google Scholar] [CrossRef]
  7. Scalisi, A.; Morandi, B.; Inglese, P.; Lo Bianco, R. Cladode growth dynamics in Opuntia ficus-indica under drought. Environ. Exp. Bot. 2016, 122, 158–167. [Google Scholar] [CrossRef]
  8. Dubeux, J.C.B.; dos Santos, M.V.F.; da Cunha, M.V.; dos Santos, D.C.; Souza, R.T.d.A.; de Mello, A.C.L.; de Souza, T.C. Cactus (Opuntia and Nopalea) nutritive value: A review. Anim. Feed Sci. Technol. 2021, 275, 114890. [Google Scholar] [CrossRef]
  9. Pessoa, D.V.; Pereira de Andrade, A.; Magalhães, A.L.R.; Teodoro, A.L.; Cordeiro dos Santos, D.; de Araújo, G.G.L.; de Medeiros, A.N.; do Nascimento, D.B.; Valença, R.d.L.; Cardoso, D.B. Forage nutritional differences within the genus Opuntia. J. Arid Environ. 2020, 181, 104243. [Google Scholar] [CrossRef]
  10. Rocha Filho, R.R.; Santos, D.C.; Véras, A.S.C.; Siqueira, M.C.B.; Monteiro, C.C.F.; Mora-Luna, R.E.; Farias, L.R.; Santos, V.L.F.; Chagas, J.C.; Ferreira, M.A. Miúda (Nopalea cochenillifera (L.) Salm-Dyck)—The Best Forage Cactus Genotype for Feeding Lactating Dairy Cows in Semiarid Regions. Animals 2021, 11, 1774. [Google Scholar] [CrossRef]
  11. Araújo Júnior, G.d.N.; da Silva, T.G.F.; de Souza, L.S.B.; de Araújo, G.G.L.; de Moura, M.S.B.; Alves, C.P.; da Silva Salvador, K.R.; de Souza, C.A.A.; de Assunção Montenegro, A.A.; da Silva, M.J. Phenophases, morphophysiological indices and cutting time in clones of the forage cacti under controlled water regimes in a semiarid environment. J. Arid Environ. 2021, 190, 104510. [Google Scholar] [CrossRef]
  12. Campos, A.R.F.; da Silva, A.J.P.; de Jong van Lier, Q.; do Nascimento, F.A.L.; Fernandes, R.D.M.; de Almeida, J.N.; Paz, V.P.d.S. Yield and morphology of forage cactus cultivars under drip irrigation management based on soil water matric potential thresholds. J. Arid Environ. 2021, 193, 104564. [Google Scholar] [CrossRef]
  13. Consoli, S.; Inglese, G.; Inglese, P. Determination of Evapotranspiration and Annual Biomass Productivity of a Cactus Pear [Opuntia ficus-indica L. (Mill.)] Orchard in a Semiarid Environment. J. Irrig. Drain. Eng. 2013, 139, 680–690. [Google Scholar] [CrossRef] [Green Version]
  14. Gibbs, J.G.; Patten, D.T. Plant temperatures and heat flux in a Sonoran Desert ecosystem. Oecologia 1970, 5, 165–184. [Google Scholar] [CrossRef] [PubMed]
  15. Heinrich, V.H.A.; Dalagnol, R.; Cassol, H.L.G.; Rosan, T.M.; de Almeida, C.T.; Silva Junior, C.H.L.; Campanharo, W.A.; House, J.I.; Sitch, S.; Hales, T.C.; et al. Large carbon sink potential of secondary forests in the Brazilian Amazon to mitigate climate change. Nat. Commun. 2021, 12, 1785. [Google Scholar] [CrossRef]
  16. Prévosto, B.; Helluy, M.; Gavinet, J.; Fernandez, C.; Balandier, P. Microclimate in Mediterranean pine forests: What is the influence of the shrub layer? Agric. For. Meteorol. 2020, 282–283, 107856. [Google Scholar] [CrossRef]
  17. Zheng, C.; Liu, S.; Song, L.; Xu, Z.; Guo, J.; Ma, Y.; Ju, Q.; Wang, J. Comparison of sensible and latent heat fluxes from optical-microwave scintillometers and eddy covariance systems with respect to surface energy balance closure. Agric. For. Meteorol. 2023, 331, 109345. [Google Scholar] [CrossRef]
  18. Perrino, E.V.; Signorile, G.; Marvulli, M. A first checklist of the vascular flora of the Polignano a mare coast (Apulia, Southern Italy). Nat. Croat. 2013, 22, 295–318. [Google Scholar]
  19. Pierini, N.A.; Vivoni, E.R.; Robles-Morua, A.; Scott, R.L.; Nearing, M.A. Using observations and a distributed hydrologic model to explore runoff thresholds linked with mesquite encroachment in the Sonoran Desert. Water Resour. Res. 2014, 50, 8191–8215. [Google Scholar] [CrossRef] [Green Version]
  20. Accogli, R.; Tomaselli, V.; Direnzo, P.; Perrino, E.V.; Albanese, G.; Urbano, M.; Laghetti, G. Edible Halophytes and Halo-Tolerant Species in Apulia Region (Southeastern Italy): Biogeography, Traditional Food Use and Potential Sustainable Crops. Plants 2023, 12, 549. [Google Scholar] [CrossRef]
  21. Jardim, A.M.d.R.F.; de Morais, J.E.F.; de Souza, L.S.B.; Lopes, D.d.C.; da Silva, M.V.; Pandorfi, H.; de Oliveira-Júnior, J.F.; da Silva, J.L.B.; Steidle Neto, A.J.; Morellato, L.P.C.; et al. A systematic review of energy and mass fluxes, and biogeochemical processes in seasonally dry tropical forests and cactus ecosystems. J. South Am. Earth Sci. 2023, 126, 104330. [Google Scholar] [CrossRef]
  22. Jardim, A.M.R.F.; Araújo Júnior, G.N.; Silva, M.V.; Santos, A.; Silva, J.L.B.; Pandorfi, H.; Oliveira-Júnior, J.F.; Teixeira, A.H.C.; Teodoro, P.E.; de Lima, J.L.M.P.; et al. Using Remote Sensing to Quantify the Joint Effects of Climate and Land Use/Land Cover Changes on the Caatinga Biome of Northeast Brazilian. Remote Sens. 2022, 14, 1911. [Google Scholar] [CrossRef]
  23. Alves, J.D.N.; Ribeiro, A.; Rody, Y.P.; Loos, R.A. Energy balance and surface decoupling factor of a pasture in the Brazilian Cerrado. Agric. For. Meteorol. 2022, 319, 108912. [Google Scholar] [CrossRef]
  24. Anderson, C.A.; Vivoni, E.R. Impact of land surface states within the flux footprint on daytime land-atmosphere coupling in two semiarid ecosystems of the Southwestern U.S. Water Resour. Res. 2016, 52, 4785–4800. [Google Scholar] [CrossRef] [Green Version]
  25. Zhang, K.; Liu, D.; Liu, H.; Lei, H.; Guo, F.; Xie, S.; Meng, X.; Huang, Q. Energy flux observation in a shrub ecosystem of a gully region of the Chinese Loess Plateau. Ecohydrol. Hydrobiol. 2022, 22, 323–336. [Google Scholar] [CrossRef]
  26. Awada, H.; Di Prima, S.; Sirca, C.; Giadrossich, F.; Marras, S.; Spano, D.; Pirastru, M. Daily Actual Evapotranspiration Estimation in a Mediterranean Ecosystem from Landsat Observations Using SEBAL Approach. Forests 2021, 12, 189. [Google Scholar] [CrossRef]
  27. Guevara-Escobar, A.; González-Sosa, E.; Cervantes-Jimenez, M.; Suzán-Azpiri, H.; Queijeiro-Bolanos, M.E.; Carrillo-Ángeles, I.; Cambron-Sandoval, V.H. Machine learning estimates of eddy covariance carbon flux in a scrub in the Mexican highland. Biogeosciences 2021, 18, 367–392. [Google Scholar] [CrossRef]
  28. Hicks, B.B.; Eash, N.S.; O’Dell, D.L.; Oetting, J.N. Augmented Bowen Ratio Analysis –I: Site Adequacy, Fetch and Heat Storage (ABRA). Agric. For. Meteorol. 2020, 290, 108035. [Google Scholar] [CrossRef]
  29. Bowen, I.S. The Ratio of Heat Losses by Conduction and by Evaporation from any Water Surface. Phys. Rev. 1926, 27, 779. [Google Scholar] [CrossRef] [Green Version]
  30. Hu, S.; Zhao, C.; Li, J.; Wang, F.; Chen, Y. Discussion and reassessment of the method used for accepting or rejecting data observed by a Bowen ratio system. Hydrol. Process. 2014, 28, 4506–4510. [Google Scholar] [CrossRef]
  31. Pokhariyal, S.; Patel, N. Evaluation of variation in radiative and turbulent fluxes over winter wheat ecosystem along Indo-Gangetic region. Arab. J. Geosci. 2021, 14, 1956. [Google Scholar] [CrossRef]
  32. Pozníková, G.; Fischer, M.; van Kesteren, B.; Orság, M.; Hlavinka, P.; Žalud, Z.; Trnka, M. Quantifying turbulent energy fluxes and evapotranspiration in agricultural field conditions: A comparison of micrometeorological methods. Agric. Water Manag. 2018, 209, 249–263. [Google Scholar] [CrossRef]
  33. Chen, J.; Dong, G.; Chen, J.; Jiang, S.; Qu, L.; Legesse, T.G.; Zhao, F.; Tong, Q.; Shao, C.; Han, X. Energy balance and partitioning over grasslands on the Mongolian Plateau. Ecol. Indic. 2022, 135, 108560. [Google Scholar] [CrossRef]
  34. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration: Guidelines for Computing Crop Water Requirements; FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998; Volume 56, ISBN 9251042195. [Google Scholar]
  35. Cheng, M.; Shi, L.; Jiao, X.; Nie, C.; Liu, S.; Yu, X.; Bai, Y.; Liu, Y.; Liu, Y.; Song, N.; et al. Up-scaling the latent heat flux from instantaneous to daily-scale: A comparison of three methods. J. Hydrol. Reg. Stud. 2022, 40, 101057. [Google Scholar] [CrossRef]
  36. Hermance, J.F.; Augustine, D.J.; Derner, J.D. Quantifying characteristic growth dynamics in a semi-arid grassland ecosystem by predicting short-term NDVI phenology from daily rainfall: A simple four parameter coupled-reservoir model. Int. J. Remote Sens. 2015, 36, 5637–5663. [Google Scholar] [CrossRef]
  37. Lewis, D.A.; Nobel, P.S.; Boyd, P.L. Thermal Energy Exchange Model and Water Loss of a Barrel Cactus, Ferocactus acanthodes. Plant Physiol. 1977, 60, 609–616. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Marques, T.V.; Mendes, K.; Mutti, P.; Medeiros, S.; Silva, L.; Perez-Marin, A.M.; Campos, S.; Lúcio, P.S.; Lima, K.; dos Reis, J.; et al. Environmental and biophysical controls of evapotranspiration from Seasonally Dry Tropical Forests (Caatinga) in the Brazilian Semiarid. Agric. For. Meteorol. 2020, 287, 107957. [Google Scholar] [CrossRef]
  39. Shao, C.; Chen, J.; Li, L.; Dong, G.; Han, J.; Abraha, M.; John, R. Grazing effects on surface energy fluxes in a desert steppe on the Mongolian Plateau. Ecol. Appl. 2017, 27, 485–502. [Google Scholar] [CrossRef]
  40. da Silva, P.F.; Lima, J.R.d.S.; Antonino, A.C.D.; Souza, R.; de Souza, E.S.; Silva, J.R.I.; Alves, E.M. Seasonal patterns of carbon dioxide, water and energy fluxes over the Caatinga and grassland in the semi-arid region of Brazil. J. Arid Environ. 2017, 147, 71–82. [Google Scholar] [CrossRef]
  41. Yue, P.; Zhang, Q.; Zhang, L.; Li, H.; Yang, Y.; Zeng, J.; Wang, S. Long-term variations in energy partitioning and evapotranspiration in a semiarid grassland in the Loess Plateau of China. Agric. For. Meteorol. 2019, 278, 107671. [Google Scholar] [CrossRef]
  42. Nassar, A.; Torres-rua, A.; Kustas, W.; Alfieri, J.; Hipps, L.; Prueger, J.; Nieto, H.; Alsina, M.M.; White, W.; McKee, L.; et al. Assessing Daily Evapotranspiration Methodologies from One-Time-of-Day sUAS and EC Information in the GRAPEX Project. Remote Sens. 2021, 13, 2887. [Google Scholar] [CrossRef]
  43. Cao, J.J.; Hu, S.; Dong, Q.; Liu, L.J.; Wang, Z.L. Green roof cooling contributed by plant species with different photosynthetic strategies. Energy Build. 2019, 195, 45–50. [Google Scholar] [CrossRef]
  44. Montesinos-Navarro, A.; Verdú, M.; Querejeta, J.I.; Valiente-Banuet, A. Nurse shrubs can receive water stored in the parenchyma of their facilitated columnar cacti. J. Arid Environ. 2019, 165, 10–15. [Google Scholar] [CrossRef]
  45. Novoa, A.; Foxcroft, L.C.; Keet, J.H.; Pyšek, P.; Le Roux, J.J. The invasive cactus Opuntia stricta creates fertility islands in African savannas and benefits from those created by native trees. Sci. Rep. 2021, 11, 20748. [Google Scholar] [CrossRef] [PubMed]
  46. Wang, M.; Wu, J.; Lafleur, P. Comparison of energy fluxes between an undisturbed bog and an adjacent abandoned peatland pasture. Agric. For. Meteorol. 2020, 291, 108086. [Google Scholar] [CrossRef]
  47. Flanagan, L.B.; Flanagan, J.E.M. Seasonal controls on ecosystem-scale CO2 and energy exchange in a Sonoran Desert characterized by the saguaro cactus (Carnegiea gigantea). Oecologia 2018, 187, 977–994. [Google Scholar] [CrossRef] [PubMed]
  48. de Lima, F.J.L.; Martins, F.R.; Costa, R.S.; Gonçalves, A.R.; dos Santos, A.P.P.; Pereira, E.B. The seasonal variability and trends for the surface solar irradiation in northeastern region of Brazil. Sustain. Energy Technol. Assess. 2019, 35, 335–346. [Google Scholar] [CrossRef]
  49. da Silva, T.G.F.; de Queiroz, M.G.; Zolnier, S.; de Souza, L.S.B.; de Souza, C.A.A.; de Moura, M.S.B.; de Araújo, G.G.L.; Steidle Neto, A.J.; dos Santos, T.S.; de Melo, A.L.; et al. Soil properties and microclimate of two predominant landscapes in the Brazilian semiarid region: Comparison between a seasonally dry tropical forest and a deforested area. Soil Tillage Res. 2021, 207, 104852. [Google Scholar] [CrossRef]
  50. Jung, M.; Koirala, S.; Weber, U.; Ichii, K.; Gans, F.; Camps-Valls, G.; Papale, D.; Schwalm, C.; Tramontana, G.; Reichstein, M. The FLUXCOM ensemble of global land-atmosphere energy fluxes. Sci. Data 2019, 6, 74. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Dhungel, R.; Aiken, R.; Evett, S.R.; Colaizzi, P.D.; Marek, G.; Moorhead, J.E.; Baumhardt, R.L.; Brauer, D.; Kutikoff, S.; Lin, X. Energy Imbalance and Evapotranspiration Hysteresis Under an Advective Environment: Evidence From Lysimeter, Eddy Covariance, and Energy Balance Modeling. Geophys. Res. Lett. 2021, 48, e2020GL091203. [Google Scholar] [CrossRef]
  52. Goldstein, G.; Ortega, J.K.E.; Nerd, A.; Nobel, P.S. Diel Patterns of Water Potential Components for the Crassulacean Acid Metabolism Plant Opuntia ficus-indica when Well-Watered or Droughted. Plant Physiol. 1991, 95, 274–280. [Google Scholar] [CrossRef] [Green Version]
  53. Lima, L.R.; Silva, T.G.F.; Jardim, A.M.R.F.; Souza, C.A.A.; Queiroz, M.G.; Tabosa, J.N. Growth, water use and efficiency of forage cactus sorghum intercropping under different water depths. Rev. Bras. Eng. Agrícola Ambient. 2018, 22, 113–118. [Google Scholar] [CrossRef] [Green Version]
  54. Han, H.; Felker, P. Field validation of water-use efficiency of the CAM plant Opuntia ellisiana in south Texas. J. Arid Environ. 1997, 36, 133–148. [Google Scholar] [CrossRef]
  55. Ouko, E.; Omondi, S.; Mugo, R.; Wahome, A.; Kasera, K.; Nkurunziza, E.; Kiema, J.; Flores, A.; Adams, E.C.; Kuraru, S.; et al. Modeling Invasive Plant Species in Kenya’s Northern Rangelands. Front. Environ. Sci. 2020, 8, 69. [Google Scholar] [CrossRef]
  56. da Silva, M.V.; Pandorfi, H.; de Almeida, G.L.P.; de Lima, R.P.; dos Santos, A.; Jardim, A.M.d.R.F.; Rolim, M.M.; da Silva, J.L.B.; Batista, P.H.D.; da Silva, R.A.B.; et al. Spatio-temporal monitoring of soil and plant indicators under forage cactus cultivation by geoprocessing in Brazilian semi-arid region. J. S. Am. Earth Sci. 2021, 107, 103155. [Google Scholar] [CrossRef]
  57. Jardim, A.M.d.R.F.; da Silva, T.G.F.; de Souza, L.S.B.; Araújo Júnior, G.d.N.; Alves, H.K.M.N.; Souza, M.d.S.; de Araújo, G.G.L.; de Moura, M.S.B. Intercropping forage cactus and sorghum in a semi-arid environment improves biological efficiency and competitive ability through interspecific complementarity. J. Arid Environ. 2021, 188, 104464. [Google Scholar] [CrossRef]
  58. Anderson, M.C.; Norman, J.M.; Mecikalski, J.R.; Otkin, J.A.; Kustas, W.P. A climatological study of evapotranspiration and moisture stress across the continental United States based on thermal remote sensing: 2. Surface moisture climatology. J. Geophys. Res. Atmos. 2007, 112, 1–13. [Google Scholar] [CrossRef]
  59. Ventura-Aguilar, R.I.; Bosquez-Molina, E.; Bautista-Baños, S.; Rivera-Cabrera, F. Cactus stem (Opuntia ficus-indica Mill): Anatomy, physiology and chemical composition with emphasis on its biofunctional properties. J. Sci. Food Agric. 2017, 97, 5065–5073. [Google Scholar] [CrossRef]
  60. Amorim, D.M.; da Silva, T.G.F.; Pereira, P.d.C.; de Souza, L.S.B.; Minuzzi, R.B. Phenophases and cutting time of forage cactus under irrigation and cropping systems. Pesqui. Agropecuária Trop. 2017, 47, 62–71. [Google Scholar] [CrossRef] [Green Version]
  61. Jardim, A.M.d.R.F.; de Souza, L.S.B.; Alves, C.P.; de Araújo, J.F.N.; de Souza, C.A.A.; Pinheiro, A.G.; de Araújo, G.G.L.; Campos, F.S.; Tabosa, J.N.; da Silva, T.G.F. Intercropping forage cactus with sorghum affects the morphophysiology and phenology of forage cactus. African J. Range Forage Sci. 2023, 40, 129–140. [Google Scholar] [CrossRef]
  62. Souza, M.d.S.; da Silva, T.G.F.; de Souza, L.S.B.; Alves, H.K.M.N.; Leite, R.M.C.; de Souza, C.A.A.; Araújo, G.G.L.D.; Campos, F.S.; Silva, M.J.D.; Souza, P.J.D.O.P. Growth, phenology and harvesting time of cactus-millet intercropping system under biotic mulching. Arch. Agron. Soil Sci. 2022, 68, 764–778. [Google Scholar] [CrossRef]
  63. Bilderback, A.H.; Torres, A.J.; Vega, M.; Ball, B.A. The structural and nutrient chemistry during early-stage decomposition and desiccation of cacti in the Sonoran Desert. J. Arid Environ. 2021, 195, 104636. [Google Scholar] [CrossRef]
  64. Melero-Meraz, V.; Zegbe, J.A.; Herrera, M.D.; Guzmán-Maldonado, S.H.; Medina-García, G.; Sánchez-Toledano, B.I.; Cruz-Bravo, R.K.; Servín-Palestina, M. On-Farm Supplemental Irrigation of “Roja Lisa” Cactus Pear: Pre- and Postharvest Effects. Horticulturae 2022, 8, 483. [Google Scholar] [CrossRef]
  65. de Cortazar, V.G.; Acevedo, E.; Nobel, P.S. Modeling of par interception and productivity by Opuntia ficus-indica. Agric. For. Meteorol. 1985, 34, 145–162. [Google Scholar] [CrossRef]
  66. Mbava, N.; Mutema, M.; Zengeni, R.; Shimelis, H.; Chaplot, V. Factors affecting crop water use efficiency: A worldwide meta-analysis. Agric. Water Manag. 2020, 228, 105878. [Google Scholar] [CrossRef]
  67. Du Toit, A.; De Wit, M.; Hugo, A. Cultivar and Harvest Month Influence the Nutrient Content of Opuntia spp. Cactus Pear Cladode Mucilage Extracts. Molecules 2018, 23, 916. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Nedjimi, B. Determination of Some Major and Trace Elements in Cladodes of Barbary fig (Opuntia ficus-indica Mill.) by X-ray Fluorescence Spectrometry. Biol. Trace Elem. Res. 2021, 199, 4353–4359. [Google Scholar] [CrossRef] [PubMed]
  69. Saraiva, F.M.; Dubeux, J.C.B.; da Cunha, M.V.; Menezes, R.S.C.; Dos Santos, M.V.F.; Camelo, D.; Ferraz, I. Manure Source and Cropping System Affect Nutrient Uptake by Cactus (Nopalea cochenillifera Salm Dyck). Agronomy 2021, 11, 1512. [Google Scholar] [CrossRef]
  70. Alvares, C.A.; Stape, J.L.; Sentelhas, P.C.; Gonçalves, J.L.d.M.; Sparovek, G. Köppen’s climate classification map for Brazil. Meteorol. Zeitschrift 2013, 22, 711–728. [Google Scholar] [CrossRef]
  71. Beck, H.E.; Zimmermann, N.E.; McVicar, T.R.; Vergopolan, N.; Berg, A.; Wood, E.F. Present and future Köppen-Geiger climate classification maps at 1-km resolution. Sci. Data 2018, 5, 180214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Leite-Filho, A.T.; de Sousa Pontes, V.Y.; Costa, M.H. Effects of Deforestation on the Onset of the Rainy Season and the Duration of Dry Spells in Southern Amazonia. J. Geophys. Res. Atmos. 2019, 124, 5268–5281. [Google Scholar] [CrossRef]
  73. Salack, S.; Klein, C.; Giannini, A.; Sarr, B.; Worou, O.N.; Belko, N.; Bliefernicht, J.; Kunstman, H. Global warming induced hybrid rainy seasons in the Sahel. Environ. Res. Lett. 2016, 11, 104008. [Google Scholar] [CrossRef]
  74. Soil Survey Staff. Keys to Soil Taxonomy, 12th ed.; USDA-Natural Resources Conservation Service: Washington, DC, USA, 2014.
  75. Heilman, J.L.; Brittin, C.L.; Neale, C.M.U. Fetch requirements for bowen ratio measurements of latent and sensible heat fluxes. Agric. For. Meteorol. 1989, 44, 261–273. [Google Scholar] [CrossRef]
  76. Gavilán, P.; Berengena, J. Accuracy of the Bowen ratio-energy balance method for measuring latent heat flux in a semiarid advective environment. Irrig. Sci. 2007, 25, 127–140. [Google Scholar] [CrossRef]
  77. Perez, P.J.; Castellvi, F.; Ibañez, M.; Rosell, J.I. Assessment of reliability of Bowen ratio method for partitioning fluxes. Agric. For. Meteorol. 1999, 97, 141–150. [Google Scholar] [CrossRef]
  78. Ortega-Farias, S.O.; Cuenca, R.H.; Ek, M. Daytime variation of sensible heat flux estimated by the bulk aerodynamic method over a grass canopy. Agric. For. Meteorol. 1996, 81, 131–143. [Google Scholar] [CrossRef]
  79. Unland, H.E.; Houser, P.R.; Shuttleworth, W.J.; Yang, Z.L. Surface flux measurement and modeling at a semi-arid Sonoran Desert site. Agric. For. Meteorol. 1996, 82, 119–153. [Google Scholar] [CrossRef] [Green Version]
  80. Mulovhedzi, N.E.; Araya, N.A.; Mengistu, M.G.; Fessehazion, M.K.; du Plooy, C.P.; Araya, H.T.; van der Laan, M. Estimating evapotranspiration and determining crop coefficients of irrigated sweet potato (Ipomoea batatas) grown in a semi-arid climate. Agric. Water Manag. 2020, 233, 106099. [Google Scholar] [CrossRef]
  81. Zhang, K.; Wang, X.; Li, Y.; Zhao, J.; Yang, Y.; Zang, H.; Zeng, Z. Peanut residue incorporation benefits crop yield, nitrogen yield, and water use efficiency of summer peanut—Winter wheat systems. F. Crop. Res. 2022, 279, 108463. [Google Scholar] [CrossRef]
  82. Raza, M.A.; Feng, L.Y.; van der Werf, W.; Cai, G.R.; Bin Khalid, M.H.; Iqbal, N.; Hassan, M.J.; Meraj, T.A.; Naeem, M.; Khan, I.; et al. Narrow-wide-row planting pattern increases the radiation use efficiency and seed yield of intercrop species in relay-intercropping system. Food Energy Secur. 2019, 8, e170. [Google Scholar] [CrossRef]
  83. Pinheiro, K.M.; Silva, T.G.F.; Carvalho, H.F.d.S.; Santos, J.E.O.; Morais, J.E.F.; Zolnier, S.; dos Santos, D.C. Correlações do índice de área do cladódio com características morfogênicas e produtivas da palma forrageira. Pesqui. Agropecuária Bras. 2014, 49, 939–947. [Google Scholar] [CrossRef] [Green Version]
  84. Rodrigues, C.R.F.; Silva, E.N.; Ferreira-Silva, S.L.; Voigt, E.L.; Viégas, R.A.; Silveira, J.A.G. High K+ supply avoids Na+ toxicity and improves photosynthesis by allowing favorable K+:Na+ ratios through the inhibition of Na+ uptake and transport to the shoots of Jatropha curcas plants. J. Plant Nutr. Soil Sci. 2013, 176, 157–164. [Google Scholar] [CrossRef]
  85. Loupassaki, M.H.; Chartzoulakis, K.S.; Digalaki, N.B.; Androulakis, I.I. Effects of salt stress on concentration of nitrogen, phosphorus, potassium, calcium, magnesium, and sodium in leaves, shoots, and roots of six olive cultivars. J. Plant Nutr. 2007, 25, 2457–2482. [Google Scholar] [CrossRef]
  86. Zhang, Q.; Song, Y.; Wu, Z.; Yan, X.; Gunina, A.; Kuzyakov, Y.; Xiong, Z. Effects of six-year biochar amendment on soil aggregation, crop growth, and nitrogen and phosphorus use efficiencies in a rice-wheat rotation. J. Clean. Prod. 2020, 242, 118435. [Google Scholar] [CrossRef]
  87. Silva, T.G.F.; Miranda, K.R.; Santos, D.C.; Queiroz, M.G.; Silva, M.C.; Cruz Neto, J.F.; Araújo, J.E.M. Área do cladódio de clones de palma forrageira: Modelagem, análise e aplicabilidade. Rev. Bras. Ciências Agrárias 2014, 9, 626–632. [Google Scholar] [CrossRef]
  88. Arnold, C.Y. The determination and significance of the base temperature in a linear heat unit system. J. Am. Soc. Hortic. Sci. 1959, 74, 430–445. [Google Scholar]
  89. Corrado, G.; Chiaiese, P.; Lucini, L.; Miras-Moreno, B.; Colla, G.; Rouphael, Y. Successive Harvests Affect Yield, Quality and Metabolic Profile of Sweet Basil (Ocimum basilicum L.). Agronomy 2020, 10, 830. [Google Scholar] [CrossRef]
  90. Anghileri, D.; Bozzini, V.; Molnar, P.; Sheffield, J. Comparison of hydrological and vegetation remote sensing datasets as proxies for rainfed maize yield in Malawi. Agric. Water Manag. 2022, 262, 107375. [Google Scholar] [CrossRef]
  91. Wolf, S.; Keenan, T.F.; Fisher, J.B.; Baldocchi, D.D.; Desai, A.R.; Richardson, A.D.; Scott, R.L.; Law, B.E.; Litvak, M.E.; Brunsell, N.A.; et al. Warm spring reduced carbon cycle impact of the 2012 US summer drought. Proc. Natl. Acad. Sci. USA 2016, 113, 5880–5885. [Google Scholar] [CrossRef]
  92. Prăvălie, R.; Sîrodoev, I.; Nita, I.A.; Patriche, C.; Dumitraşcu, M.; Roşca, B.; Tişcovschi, A.; Bandoc, G.; Săvulescu, I.; Mănoiu, V.; et al. NDVI-based ecological dynamics of forest vegetation and its relationship to climate change in Romania during 1987–2018. Ecol. Indic. 2022, 136, 108629. [Google Scholar] [CrossRef]
  93. Hurni, K.; Van Den Hoek, J.; Fox, J. Assessing the spatial, spectral, and temporal consistency of topographically corrected Landsat time series composites across the mountainous forests of Nepal. Remote Sens. Environ. 2019, 231, 111225. [Google Scholar] [CrossRef]
  94. Foga, S.; Scaramuzza, P.L.; Guo, S.; Zhu, Z.; Dilley, R.D.; Beckmann, T.; Schmidt, G.L.; Dwyer, J.L.; Joseph Hughes, M.; Laue, B. Cloud detection algorithm comparison and validation for operational Landsat data products. Remote Sens. Environ. 2017, 194, 379–390. [Google Scholar] [CrossRef] [Green Version]
  95. Bielski, C.H.; Twidwell, D.; Fuhlendorf, S.D.; Wonkka, C.L.; Allred, B.W.; Ochsner, T.E.; Krueger, E.S.; Carlson, J.D.; Engle, D.M. Pyric herbivory, scales of heterogeneity and drought. Funct. Ecol. 2018, 32, 1599–1608. [Google Scholar] [CrossRef]
  96. Lyra, G.B.; de Souza, J.L.; da Silva, E.C.; Lyra, G.B.; Teodoro, I.; Ferreira-Júnior, R.A.; de Souza, R.C. Soil water stress co-efficient for estimating actual evapotranspiration of maize in northeastern Brazil. Meteorol. Appl. 2016, 23, 26–34. [Google Scholar] [CrossRef] [Green Version]
  97. Pardo, J.J.; Martínez-Romero, A.; Léllis, B.C.; Tarjuelo, J.M.; Domínguez, A. Effect of the optimized regulated deficit irrigation methodology on water use in barley under semiarid conditions. Agric. Water Manag. 2020, 228, 105925. [Google Scholar] [CrossRef]
  98. Thornthwaite, C.W.; Mather, J.R. The Water Balance; Publications in Climatology; Laboratory of Climatology: Centerton, NJ, USA, 1955; Volume 8. [Google Scholar]
  99. Aparecido, L.E.d.O.; Lorençone, J.A.; Lorençone, P.A.; de Meneses, K.C.; de Moraes, J.R.d.S.C. Climate risk to peanut cultivation in Brazil across different planting seasons. J. Sci. Food Agric. 2021, 101, 5002–5015. [Google Scholar] [CrossRef]
  100. Almagbile, A.; Zeitoun, M.; Hazaymeh, K.; Sammour, H.A.; Sababha, N. Statistical analysis of estimated and observed soil moisture in sub-humid climate in north-western Jordan. Environ. Monit. Assess. 2019, 191, 96. [Google Scholar] [CrossRef]
  101. R Core Team. R: The R Project for Statistical Computing, Vienna, Austria; v. 4.1.3, 2022. Available online: https://www.r-project.org/ (accessed on 23 May 2023).
Figure 1. Hourly daytime mean energy flux over the area of cactus for different periods: wet season (a,e,i), wet-dry transition (b,f), dry season (c,g) and dry-wet transition (d,h). Here, Rn is the net radiation, LE is the latent heat flux, H is the sensible heat flux and G is the soil heat flux, all calculated in W m−2. Note: panels (ac) show the seasons for 2015; (dg) for 2016; and (h,i) for 2017.
Figure 1. Hourly daytime mean energy flux over the area of cactus for different periods: wet season (a,e,i), wet-dry transition (b,f), dry season (c,g) and dry-wet transition (d,h). Here, Rn is the net radiation, LE is the latent heat flux, H is the sensible heat flux and G is the soil heat flux, all calculated in W m−2. Note: panels (ac) show the seasons for 2015; (dg) for 2016; and (h,i) for 2017.
Plants 12 02562 g001
Figure 2. Time series for the energy budget in an area of cactus from 2015 to 2017. Rg is the global solar radiation (a), Rn is the net radiation (b), LE is the latent heat flux (c), H is the sensible heat flux (d), G is the soil heat flux (e) and ET is the evapotranspiration (f). The boxplots show the median; horizontal bars represent the 25th, 50th and 75th percentiles; whiskers (lower and upper) represent the 1.5× interquartile ranges. Corresponding data are represented by circles.
Figure 2. Time series for the energy budget in an area of cactus from 2015 to 2017. Rg is the global solar radiation (a), Rn is the net radiation (b), LE is the latent heat flux (c), H is the sensible heat flux (d), G is the soil heat flux (e) and ET is the evapotranspiration (f). The boxplots show the median; horizontal bars represent the 25th, 50th and 75th percentiles; whiskers (lower and upper) represent the 1.5× interquartile ranges. Corresponding data are represented by circles.
Plants 12 02562 g002
Figure 3. Boxplot of seasonal variations in energy and water exchange in an area of cactus. Rn is the net radiation (a), LE is the latent heat flux (b), H is the sensible heat flux (c), G is the soil heat flux (d), ET is the evapotranspiration (e) and VPD is the vapour pressure deficit (f). The boxplots show the median; horizontal bars represent the 25th, 50th and 75th percentiles; whiskers (lower and upper) represent the 1.5× interquartile ranges; dots represent outliers. Significance was calculated using one-way analysis of variance (ANOVA) with Tukey’s honestly significant difference (HSD) post hoc test. Different letters above each box indicate significant differences (p < 0.05).
Figure 3. Boxplot of seasonal variations in energy and water exchange in an area of cactus. Rn is the net radiation (a), LE is the latent heat flux (b), H is the sensible heat flux (c), G is the soil heat flux (d), ET is the evapotranspiration (e) and VPD is the vapour pressure deficit (f). The boxplots show the median; horizontal bars represent the 25th, 50th and 75th percentiles; whiskers (lower and upper) represent the 1.5× interquartile ranges; dots represent outliers. Significance was calculated using one-way analysis of variance (ANOVA) with Tukey’s honestly significant difference (HSD) post hoc test. Different letters above each box indicate significant differences (p < 0.05).
Plants 12 02562 g003
Figure 4. Temporal evolution of the Normalised Difference Vegetation Index [NDVI] (a), reference evapotranspiration [ET0] (c), Evaporative Stress Index [ESI] (e), available water fraction [AWF] (g) and rainfall (i) for an area cultivated with cactus. The five panels (b,d,f,h,j) show monthly results over each year for NDVI, ET0, ESI, AWF and rainfall during the experimental period, respectively. Data with error bars represent the mean ± SD (standard deviation).
Figure 4. Temporal evolution of the Normalised Difference Vegetation Index [NDVI] (a), reference evapotranspiration [ET0] (c), Evaporative Stress Index [ESI] (e), available water fraction [AWF] (g) and rainfall (i) for an area cultivated with cactus. The five panels (b,d,f,h,j) show monthly results over each year for NDVI, ET0, ESI, AWF and rainfall during the experimental period, respectively. Data with error bars represent the mean ± SD (standard deviation).
Plants 12 02562 g004
Figure 5. Absolute growth rate—AGR and relative growth rate—RGR (a), net assimilation rate—NAR (b), specific cladode area—SCA (c) and cladode emission rate (d) in N. cochenillifera (L.) Salm-Dyck under a semi-arid environment.
Figure 5. Absolute growth rate—AGR and relative growth rate—RGR (a), net assimilation rate—NAR (b), specific cladode area—SCA (c) and cladode emission rate (d) in N. cochenillifera (L.) Salm-Dyck under a semi-arid environment.
Plants 12 02562 g005
Figure 6. Cladode water content (a) and cladode succulence (b) in N. cochenillifera (L.) Salm-Dyck over time (2015 to 2017) under a semi-arid environment. The boxplots show the median; horizontal bars represent the 25th, 50th and 75th percentiles; whiskers (lower and upper) represent the 1.5× interquartile ranges. The colouring of each boxplot indicates the sampling period.
Figure 6. Cladode water content (a) and cladode succulence (b) in N. cochenillifera (L.) Salm-Dyck over time (2015 to 2017) under a semi-arid environment. The boxplots show the median; horizontal bars represent the 25th, 50th and 75th percentiles; whiskers (lower and upper) represent the 1.5× interquartile ranges. The colouring of each boxplot indicates the sampling period.
Plants 12 02562 g006
Figure 7. Principal component analysis (PCA) ordered biplot of environmental and plant factors. Score plot (a) and loading plot (b) of the first two principal components (PC) during the wet and dry seasons and transition periods. P (phosphorus), Mg2+ (magnesium), Ca2+ (calcium), K+ (potassium) and Na+ (sodium) refer to the efficiency of use of each nutrient. The following abbreviations are used: net radiation (Rn), latent heat flux (LE), sensible heat flux (H), soil heat flux (G), evapotranspiration (ET), vapour pressure deficit (VPD), Evaporative Stress Index (ESI), Normalised Difference Vegetation Index (NDVI), water use efficiency (WUE), reference evapotranspiration (ET0), absolute growth rate (AGR), relative growth rate (RGR), net assimilation rate (NAR), specific cladode area (SCA), cladode water content (CWC), cladode succulence (CS), radiation use efficiency (RUE) and biomass yield (Yield).
Figure 7. Principal component analysis (PCA) ordered biplot of environmental and plant factors. Score plot (a) and loading plot (b) of the first two principal components (PC) during the wet and dry seasons and transition periods. P (phosphorus), Mg2+ (magnesium), Ca2+ (calcium), K+ (potassium) and Na+ (sodium) refer to the efficiency of use of each nutrient. The following abbreviations are used: net radiation (Rn), latent heat flux (LE), sensible heat flux (H), soil heat flux (G), evapotranspiration (ET), vapour pressure deficit (VPD), Evaporative Stress Index (ESI), Normalised Difference Vegetation Index (NDVI), water use efficiency (WUE), reference evapotranspiration (ET0), absolute growth rate (AGR), relative growth rate (RGR), net assimilation rate (NAR), specific cladode area (SCA), cladode water content (CWC), cladode succulence (CS), radiation use efficiency (RUE) and biomass yield (Yield).
Plants 12 02562 g007
Figure 8. Geographic location of the study area in the district of Mirandiba, Pernambuco, Brazil. UTM projection (zone 24S), with the SIRGAS2000 Datum.
Figure 8. Geographic location of the study area in the district of Mirandiba, Pernambuco, Brazil. UTM projection (zone 24S), with the SIRGAS2000 Datum.
Plants 12 02562 g008
Table 1. Energy budget component partitioning during the wet, wet-dry, dry and dry-wet seasons in an area of cactus.
Table 1. Energy budget component partitioning during the wet, wet-dry, dry and dry-wet seasons in an area of cactus.
SeasonPartition Ratio (%)
LE/RnH/RnG/Rn
Wet305713
Wet-dry235819
Dry265519
Dry-wet176419
Mean ratio2458.517.5
Rn is the net radiation; LE is the latent heat flux; H is the sensible heat flux; G is the soil heat flux.
Table 2. Summary of the efficiency of the biophysical parameters, Bowen ratio and yield in the cactus at the end of the experimental period.
Table 2. Summary of the efficiency of the biophysical parameters, Bowen ratio and yield in the cactus at the end of the experimental period.
VariableMean ValueStandard Deviation
β (dimensionless)3.53±2.45
Yield (Mg ha−1)12.47±2.20
RUE (g MJ−1)3.95±0.70
WUE (kg m−3)1.75±0.31
NUE[Calcium] (mg m−3)29.03±1.14
NUE[Magnesium] (mg m−3)12.44±2.25
NUE[Phosphorus] (mg m−3)62.13±11.71
NUE[Potassium] (mg m−3)50.88±5.50
NUE[Sodium] (mg m−3)0.22±0.02
RUE is the radiation use efficiency; β is the annual Bowen ratio based on annual energy budgets; WUE is the water use efficiency; NUE is the nutrient use efficiency; Yield is the yield of aboveground dry biomass.
Table 3. Consistency criteria for data generated by the Bowen ratio method under non-advective conditions.
Table 3. Consistency criteria for data generated by the Bowen ratio method under non-advective conditions.
Available Energy Vapour Pressure GradientBowen RatioHeat Flux
RnG > 0e > 0β > −1LE > 0 and H ≤ 0 for −1 < β ≤ 0 or H > 0 for β > 0
e < 0β < −1LE < 0 and H > 0
RnG < 0e > 0β > −1LE > 0 and H < 0
e < 0β < −1LE < 0 and H ≥ 0 for −1 < β ≤ 0 or H < 0 for β > 0
data considered satisfactory as described by Perez et al. [77]. Rn is the net radiation, G is the soil heat flux, ∆e is the difference in vapour pressure between the two measurement heights, β is the Bowen ratio, LE is the latent heat flux and H is the sensible heat flux.
Table 4. Summary of error types when faults occur using the Bowen ratio-energy balance method.
Table 4. Summary of error types when faults occur using the Bowen ratio-energy balance method.
Type of ErrorApplied Conditions
ARnG > 0, ∆e > 0 and β < −1 + |ε|
BRnG > 0, ∆e < 0 and β > −1 − |ε|
CRnG < 0, ∆e > 0 and β > −1 − |ε|
DRnG < 0, ∆e < 0 and β < −1 + |ε|
EReferring to a rapid change in temperature and vapour pressure.
β is the Bowen ratio, and ε is the error interval that defines the threshold for excluding values of the Bowen ratio close to −1 Perez et al. [77].
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jardim, A.M.d.R.F.; Morais, J.E.F.d.; Souza, L.S.B.d.; Souza, C.A.A.d.; Araújo Júnior, G.d.N.; Alves, C.P.; Silva, G.Í.N.d.; Leite, R.M.C.; Moura, M.S.B.d.; de Lima, J.L.M.P.; et al. Monitoring Energy Balance, Turbulent Flux Partitioning, Evapotranspiration and Biophysical Parameters of Nopalea cochenillifera (Cactaceae) in the Brazilian Semi-Arid Environment. Plants 2023, 12, 2562. https://doi.org/10.3390/plants12132562

AMA Style

Jardim AMdRF, Morais JEFd, Souza LSBd, Souza CAAd, Araújo Júnior GdN, Alves CP, Silva GÍNd, Leite RMC, Moura MSBd, de Lima JLMP, et al. Monitoring Energy Balance, Turbulent Flux Partitioning, Evapotranspiration and Biophysical Parameters of Nopalea cochenillifera (Cactaceae) in the Brazilian Semi-Arid Environment. Plants. 2023; 12(13):2562. https://doi.org/10.3390/plants12132562

Chicago/Turabian Style

Jardim, Alexandre Maniçoba da Rosa Ferraz, José Edson Florentino de Morais, Luciana Sandra Bastos de Souza, Carlos André Alves de Souza, George do Nascimento Araújo Júnior, Cléber Pereira Alves, Gabriel Ítalo Novaes da Silva, Renan Matheus Cordeiro Leite, Magna Soelma Beserra de Moura, João L. M. P. de Lima, and et al. 2023. "Monitoring Energy Balance, Turbulent Flux Partitioning, Evapotranspiration and Biophysical Parameters of Nopalea cochenillifera (Cactaceae) in the Brazilian Semi-Arid Environment" Plants 12, no. 13: 2562. https://doi.org/10.3390/plants12132562

APA Style

Jardim, A. M. d. R. F., Morais, J. E. F. d., Souza, L. S. B. d., Souza, C. A. A. d., Araújo Júnior, G. d. N., Alves, C. P., Silva, G. Í. N. d., Leite, R. M. C., Moura, M. S. B. d., de Lima, J. L. M. P., & Silva, T. G. F. d. (2023). Monitoring Energy Balance, Turbulent Flux Partitioning, Evapotranspiration and Biophysical Parameters of Nopalea cochenillifera (Cactaceae) in the Brazilian Semi-Arid Environment. Plants, 12(13), 2562. https://doi.org/10.3390/plants12132562

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