Next Article in Journal
Inferring Mass Loss by Measuring Contemporaneous Deformation around the Helheim Glacier, Southeastern Greenland, Using Sentinel-1 InSAR
Next Article in Special Issue
Remote Sensing of Land Surface Phenology: Editorial
Previous Article in Journal
Object-Oriented Unsupervised Classification of PolSAR Images Based on Image Block
Previous Article in Special Issue
Influences of Seasonal Soil Moisture and Temperature on Vegetation Phenology in the Qilian Mountains
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Seasonal Ecosystem Productivity in a Seasonally Dry Tropical Forest (Caatinga) Using Flux Tower Measurements and Remote Sensing Data

by
Gabriel Brito Costa
1,2,3,*,
Keila Rêgo Mendes
2,4,
Losany Branches Viana
5,
Gabriele Vieira Almeida
6,
Pedro Rodrigues Mutti
2,
Cláudio Moisés Santos e Silva
2,
Bergson Guedes Bezerra
2,
Thiago Valentim Marques
2,7,
Rosária Rodrigues Ferreira
2,
Cristiano Prestelo Oliveira
2,4,
Weber Andrade Gonçalves
2,4,
Pablo Eli Oliveira
2,
Suany Campos
2,
Maria Uilhiana Gomes Andrade
2,
Antônio Celso Dantas Antonino
8 and
Rômulo Simões Cézar Menezes
8
1
Biosciences Post-Graduate Program (PPG-BIO), Federal University of Western Pará (UFOPA), Santarem 68035-110, Brazil
2
Climate Sciences Post-Graduate Program (PPGCC), Federal University of Rio Grande do Norte, Av. Senador Salgado Filho, 3000, Lagoa Nova, Natal 59078-970, Brazil
3
Anthropic Studies in the Amazon Post-Graduate Program (PPGEAA), Federal University of Pará, Castanhal 68740-222, Brazil
4
Department of Atmospheric and Climate Sciences, Federal University of Rio Grande do Norte, Av. Senador Salgado Filho, 3000, Lagoa Nova, Natal 59078-970, Brazil
5
Institute of Engineering and Geosciences, Federal University of West Pará, Rua Vera Paz s/n, Salé, Santarem 68035-110, Brazil
6
Institute of Biodiversity and Forests, Federal University of West Pará, Rua Vera Paz s/n, Salé, Santarem 68035-110, Brazil
7
Federal Institute of Education, Science and Technology of Rio Grande do Norte, IFRN, Natal 59628-330, Brazil
8
Department of Nuclear Energy, Federal University of Pernambuco, Recife 50740-545, Brazil
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(16), 3955; https://0-doi-org.brum.beds.ac.uk/10.3390/rs14163955
Submission received: 30 July 2022 / Revised: 7 August 2022 / Accepted: 10 August 2022 / Published: 15 August 2022
(This article belongs to the Special Issue Remote Sensing of Land Surface Phenology)

Abstract

:
The Caatinga dry forest encompasses 11% of the total continental territory of Brazil. Nevertheless, most research on the relationship between phenology and ecosystem productivity of Brazilian tropical forests is aimed at the Amazon basin. Thus, in this study we evaluated the seasonality of ecosystem productivity (gross primary production—GPP) in a preserved Caatinga environment in northeast Brazil. Analyses were carried out using eddy covariance measurements and satellite-derived data from sensor MODIS (MODerate Resolution Imaging Spectroradiometer, MOD17 and MOD13 products). In addition to GPP, we investigated water use efficiency (WUE) and meteorological and phenological aspects through remotely sensed vegetation indices (NDVI and EVI). We verified that ecosystem productivity is limited mainly by evapotranspiration, with maximum GPP values registered in the wetter months, indicating a strong dependency on water availability. NDVI and EVI were positively associated with GPP (r = 0.69 and 0.81, respectively), suggesting a coupling between the emergence of new leaves and the phenology of local photosynthetic capacity. WUE, on the other hand, was strongly controlled by consecutive dry days and not necessarily by total precipitation amount. The vegetation indices adequately described interannual variations of the forest response to environmental factors, and GPP MODIS presented a good relationship with tower-measured GPP in dry (R2 = 0.76) and wet (R2 = 0.62) periods.

Graphical Abstract

1. Introduction

Deforestation and land use change affect the preservation of tropical forests, demanding great efforts from the scientific community to investigate and discover the main natural controls of biophysical and phenological aspects in these ecosystems. With that, it is possible to better understand their ecological functions, which in turn allows the extrapolation of in situ observations to larger areas through remote sensing data, for example. This bigger picture provides more accurate and detailed scenarios regarding the impacts of climate change and land use change on landscape and the environment. In this perspective, studies on the characterization of the biosphere–atmosphere interactions are crucial to highlight the importance of tropical forests to the global climate balance [1,2,3,4]. Given the particularities and the variety of physiognomies within a single biome, these studies are also useful for the development of more accurate biosphere–atmosphere models [5,6,7], for the validation of remote sensing products [8,9,10,11], for the characterization of vegetation response to the natural seasonality of the environment [7,12,13], and they can also be used as proxies in the analysis of climate change impacts and implications.
In Brazil, most research in this field is concentrated in the Amazon region [14,15,16,17,18]. Thus, there is a relative scarcity of studies on the seasonality of phenological aspects of other Brazilian biomes and their atmospheric controls, or even studies regarding the validation of remote sensing data over these areas. However, recent studies indicate that other Brazilian biomes such as the Caatinga behave as efficient carbon sinks if compared to other dry forests across the globe [19]. These authors show that even in the dry season, the Caatinga uptakes from 145 to 169 gC m−2 year−1 due to the low respiration rates of the ecosystem. Nevertheless, there is little point-based experimental evidence apart from these results and, therefore, remote sensing estimates can be important sources of information on the ecological functions of the Caatinga, as suggested by the study by Ferreira [20]. Although observations at a single point are not representative for the entire ecosystem, such observations are useful to improve the understanding of ecosystem productivity and gas exchange in different biomes, being a common analysis in the literature [2,4,19] in studies using the eddy covariance technique in a single tower, given the high cost of keeping many experiments covering the same area. Such data can serve as a subsidy for several actions, such as validating remote sensing products for larger areas.
Remote sensing methods provide means to monitor carbon dioxide (CO2) exchange and several other vegetation growth parameters in almost all kinds of ecosystems, from boreal forests [21] to the different mosaic vegetation patterns of tropical forests [9]. Furthermore, they may be associated with information collected through eddy covariance systems, providing more robust measurements to characterize energy and water fluxes in ecosystems (e.g., Maselli [22]). In addition, due to limitations in establishing long-term series in these environments, measurements in shorter periods (2 years or more) are already efficient in characterizing the seasonality of local variables and validating remote sensing products [2,3,19,20].
Carbons fluxes, including photosynthesis and respiration, are strongly coupled to the water cycle in the various terrestrial ecosystems [23,24]. Water-use efficiency (WUE) is commonly defined as the ratio between gross primary production (GPP) and evapotranspiration (ET) at the ecosystem level [25], representing an important ecophysiological index relating water and carbon cycles. At the ecosystem level, WUE is typically controlled by climate and soil variables, including precipitation, air temperature, vapor pressure deficit (VPD), and soil water content, due to their roles in energy partitioning and canopy conductance (e.g., Tong an Song [26,27]), although the magnitude and direction of the response might be different or even opposite (e.g., Liu [28]).
As observed in many semi-arid ecosystems, studies on the dynamics of CO2 exchange in the Caatinga biome are still incipient. Regarding the role of the Caatinga in the global and regional carbon balance, there is still need to corroborate in situ GPP measurements, to comprehend phenological controls and WUE patterns in the biome, and to better understand the Caatinga response to years of extreme drought such as 2015 [29]. Therefore, this study proposes an assessment of the seasonality of ecosystem production and its correlations with atmospheric controls (meteorological variables), vegetation indices (NDVI—Normalized Difference Vegetation Index and EVI—Enhanced Vegetation Index, which are the most used greenness proxies in the literature according to Helman [30]), and GPP data estimated by remote sensing. Our objective is to answer the following questions: (i) what are the correlations, and which are the main controls and contributions of evaporative fluxes to the seasonality of Caatinga GPP? (ii) Can phenological patterns observed at the field level also be perceived at larger scales (MODIS NDVI and EVI)? (iii) What are the patterns of WUE in the Caatinga and how do they behave during climate extreme years such as 2015?

2. Materials and Methods

2.1. Description of the Study Area

The study was conducted in a preserved fragment of the Caatinga biome, which is a seasonally dry tropical forest (SDTF), located at the Seridó Ecological Station (ESEC-Seridó) (6°34′4″S, 37°15′0″W), in the Rio Grande do Norte state (Figure 1), Brazilian Semiarid region. The ESEC-Seridó is a federal conservation unit of the Caatinga biome, managed by the Chico Mendes Institute for Biodiversity Conservation (ICMBio), with an area of 1163 ha of preserved forest.
The vegetation is composed of dry, xerophyte, deciduous and semi-deciduous species, with shrublike and arboreous structure, sparsely distributed and reaching up to 8 m in height, in addition to herb patches that thrive only during the wet season [31]. The climate in the region is semiarid, with approximately 700 mm of total annual rainfall distributed during the wet season (January to May). The mean annual temperature is of approximately 25.0 °C, and the accumulated evapotranspiration rates are high all year round (1500 to 2000 mm yr−1) and relative humidity of the air is of approximately 60% [32].

2.2. Eddy Covariance Climate and Carbon Flux Data

Experimental data were collected in the period from January 2014 to December 2015 through an eddy covariance system installed in an 11 m height flux tower belonging to the Brazilian National Institute of the Semiarid (INSA) and integrating the National Observatory of Water and Carbon Dynamics in the Caatinga Biome (NOWCDCB) network. The analyzed variables were daily and monthly means of air temperature (Ta, °C), net ecosystem exchange (NEE, gC m−2 day−1), gross primary production (GPP, gC m−2 day−1), ecosystem respiration (Reco, gC m−2 day−1), vapor pressure deficit (VPD, kPa), relative humidity (RH, %), rainfall (mm), latent heat flux (LE, in W m−2), sensible heat flux (H, in W m−2)m and net radiation (Rn, in W m−2). High (10 Hz) and low frequency data (sampled every 5 s and stored as half-hour averages) were recorded. The high frequency data consist of measurements of CO2 concentration and water vapor and the three components of wind speed (u_x, u_y, u_z), using an Integrated CO2/H2O Open-Path Gas Analyzer and 3D Sonic Anemometer (IRGASON, Campbell Scientific, Inc., Logan, UT, USA). Atmospheric pressure was measured using an Enhanced Barometer PTB110 (Vaisala Corporation, Helsink, Finland) and air temperature by an HMP155A probe (Vaisala Corporation, Helsink, Finland) of low frequency include net radiation (Rn), air temperature (Tair), and relative humidity (RH), in addition to precipitation. Rn measurements were performed using a net radiometer model CNR4 (Kipp and Zonen B. V., Delft, Netherlands Ta) and RH data were measured using a model HMP45C temperature and relative humidity probe (Vaisala Corporation, Helsinki, Finland).
Net Ecosystem Exchange (hereafter referred as NEE) is the sum of the eddy CO2 flux (FCO2), calculated as the covariance between the fluctuations of the vertical wind speed (w′) and the density of CO2 (c′) and the rate of change of CO2 stored in the air column below the EC measurements height (Sc). Since no concentration profile was installed at the site, we opted for the discrete approach, assuming that the CO2 concentration inside the canopy can be estimated as an approximation [33]
CO2 concentrations in the form of Net Ecosystem Exchange (NEE) were partitioned into Gross Primary Production (GPP) and Ecosystem Respiration (Reco) which is equivalent to the sum of autotrophic and heterotrophic respiration, using a flux partitioning method based on night time as described by Reichstein [34]. Thus, the GPP was estimated according to Equation (1):
GPP = NEE − Reco
in which Reco and GPP were calculated using the online tool provided by the Max Plank Institute (Max Planck Institute for Biogeochemistry—http://www.bgc-jena.mpg.de/~MDIwork/eddyproc/ accessed on 24 April 2022).
The NEE was modeled based on daytime data using the common rectangular hyperbolic light response curve model [35]:
NEE = (α∙β∙Rg)/(α∙Rg + β) + γ
where α (μmol C 〖J〗−1) is the light utilization efficiency and represents the initial slope of the light response curve, β (μmol C m−2 s−1) is the maximum CO2 absorption rate of the canopy at light saturation, γ (μmol C m−2 s−1) is the ecosystem respiration, and Rg (W m−2) is the global radiation. Full details on the instrumentation, energy and carbon fluxes partitioning and corrections, and the eddy covariance method in general can be consulted in previous literature using the same dataset [2,3,19].
Caatinga ET was estimated through latent heat as follows (Equation (1)):
ET = 86 , 400 LE L
where L (J kg−1) = 103 × (2500 − 2.37 × Ta) and 86,400 is the daily integration factor.
The Water Use Efficiency was calculated as follows [28]:
WUE = GPP ET

2.3. MOD13A2 Vegetation Indices Data

Vegetation indices were obtained through the remote sensing MOD13A2 v6—Terra Vegetation Indices product from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor onboard the Terra satellite. The MOD13A2 product provides 16-day composites of NDVI and EVI data at 1 km spatial resolution, which were both used in this study [36]. Data are downloadable at the following website: https://lpdaac.usgs.gov/products/mod13a2v006/ accessed on 22 April 2022. The indices were extracted from the MODIS product at the exact location of the ESEC-Seridó flux tower.

2.4. MOD17A2 GPP Data

Remotely sensed GPP data obtained from the MOD17A2 version 6.0 product [37], derived from the MODIS sensor onboard the Terra satellite were used. This product provides 8-day composites data at a 500 m spatial resolution and is available for download in https://modis.gsfc.nasa.gov/data/dataprod/mod17.php accessed on 22 April 2022. These data are calculated based on the concept of light use efficiency [38], through the relation between incident photosynthetically active radiation (PAR), the fraction of photosynthetically active radiation absorbed by plants (FPAR), and the actual light use efficiency (Ɛ) of vegetation [38].
GPP = Ɛ ×   APAR
where APAR is the absorbed photosynthetically active radiation, which is calculated as the product between the FPAR—derived from the MOD15A2H product [39]—and the PAR. PAR values are obtained from the Global Modeling Assimilation Office (GMAO) reanalysis [40], set and corresponding to 45% of the total incident solar energy in the visible spectrum (0.4 to 0.7 μm). In the MOD17A2 product, Ɛ values are derived from the attenuation of its maximum value (εmax) due to two environmental stresses: (1) minimum temperature, which can inhibit photosynthesis, reducing enzymatic activity, and (2) VPD, because high VPD can reduce stomatal conductance [41].

2.5. Statistical Processing

Statistical analysis and data processing were carried out using the R software, version 4.0.5 [42]. High frequency data were processed as 30-min averages using the EdiRe software (version 1.3.2.11, University of Edinburgh, Edinburgh, Scotland). Pearson’s correlation coefficient (r) and simple linear regression models were used to assess the relationship between pairs of variables. The coefficients of the models were evaluated individually through Student’s t-test at the 5% significance level (p-value < 0.05) and the 95% confidence intervals were used to display the confidence limits of the fitted model. Additionally, the quality of the model’s fit was evaluated through the proportion of the variance of the response variable explained by the model via coefficient of determination (R2) and the significance (p-value < 0.05) of the model when compared to the fit containing only the intercept. Principal component analysis (PCA) was also applied between meteorological parameters and the components of CO2 fluxes (GPP, Reco and NEE) to explain the variance and covariance structure of data through non-correlated linear combinations of the original dataset. In this study, we opted to use the correlation matrix to create these combinations, since PCA is sensible to variables with different scales and units [43].

3. Results

3.1. Seasonal Dynamics of Climate, Vegetation Indices, Carbon, and Energy Fluxes

Figure 2 shows the mean monthly accumulated evapotranspiration, vapor pressure deficit (Figure 2a), and air temperature and relative humidity averages (Figure 2b), as well as monthly accumulated rainfall. Evapotranspiration, relative humidity of the air, and rainfall are inversely related to VPD and air temperature, which peaked during the dry season of 2015 (29.8 °C, Figure 2b). Air temperature varied 11% between monthly minimum and maximum values, while relative humidity varied more sharply, roughly 40% between minimum and maximum values. ET was higher in April (94.4 ± 8.5 kg H2O m−2 month−1), coinciding with the highest rainfall values. Mean ET during the wet season was 280.9 ± 23.9 kg H2O m−2, which is 225% higher than in the dry season (86.4 ± 8.4 kg H2O m−2). At the annual scale, mean VPD in 2015 was 2.6 kPa while in 2014, it was 2.2 kPa, which may be explained by the slightly lower precipitation (513 and 466.5 mm in 2014 and 2015, respectively) and higher air temperature, with a consequent decrease in relative humidity of the air between these years (Table 1). In the wet season, 83% of mean annual rainfall occurred, which accounted for 53 days in which rainfall occurred in 2014 and 26 days in 2015.
The seasonal dynamics of the vegetation indices (EVI and NDVI) show the pattern of surface phenology regarding canopy structure and the development of the vegetated area (Figure 3a). We can observe that the seasonality of vegetation is directly related to local precipitation patterns (Figure 2). The NDVI and EVI presented higher values at the peak of the wet season (NDVI > 0.6; EVI > 0.4), subsequently reducing in magnitude with the establishment of the dry season, reaching NDVI < 0.4 and EVI < 0.2 and indicating the beginning of leaf senescence in the Caatinga. The occasional occurrence of rainfall events in the dry–wet transition season were sufficient to increase vegetation activity in November 2014 (NDVI > 0.5; EVI > 0.3), which led to an increase in latent heat flux, suggesting that the onset of the wet season exerts a strong control in the production of biomass in this ecosystem. These data are strongly related to energy fluxes, with higher NDVI and EVI values occurring in the wet season when evapotranspiration (and LE) is also higher (Figure 3b).
Net radiation peaked at roughly 200 W m−2 day−1 in November. Sensible heat flux accounts for most of the energy balance partitioning during the dry season, while latent heat flux is superior during the wet season. Table 1 shows the mean and deviation of the studied variables in each experiment year.

3.2. Seasonal Dynamics of the Carbon Cycle, WUE, Tower-Measured GPP, and Remotely Sensed GPP

Carbon fluxes (GPP, Reco, and NEE) and WUE are presented in Figure 4. The seasonality of these parameters is evident, as well as a relative coupling with observed rainfall pulses. GPP peaked in April 2014 and the entire period monthly mean was of 1.89 ± 0.87 gC m−2 month−1. MODIS GPP presented a similar behavior, but with slightly higher GPP values, with a maximum of 10.0 ± 0.84 gC m−2 month−1 in April 2014. Monthly mean GPP ranged from 1.0 to 3.7 gC m−2 month−1. On average, accumulated GPP during the wet season (163.9 gC m−2 month−1) was 126% higher than in the dry season (72.5 gC m−2 month−1). The lowest carbon fluxes occurred in the dry season, when rainfall accounts for roughly 17% of the total annual expected amount of precipitation. Total accumulated ecosystem productivity (estimated through GPP) was 414.7 gC m−2 yr−1 in 2014 and 334.0 gC m−2 yr−1 in 2015. The seasonal pattern of NEE was also modulated by local hydrometeorological conditions, and the Caatinga site behaved as a carbon sink (NEE < 0). Higher carbon uptake was observed during the wetter months (Figure 4a). With the end of the wet season, uptake diminished, but the Caatinga still behaved as a carbon sink until the end of each year. Carbon uptake was different between the wet season of both studied years and this is probably not associated with rainfall differences between years since it was roughly 10%. The most plausible explanation is the difference in the number of rain days in each year (Figure 4b), with a 49% decrease from 2014 to 2015. This hypothesis agrees with the higher WUE observed in 2015 (1085.2 gC kg−1 H2O) in relation to 2014 (827.3 gC kg−1 H2O). Indeed, WUE represents the trade-off between carbon gains and water consumption during photosynthesis, and the ecosystem tends to shift its behavior towards a more efficient use of water in drier conditions (Figure 4b). In both years, WUE was inversely related to ET. The computation of WUE is important to characterize specific water cycles and the effect of drought on the water balance and carbon sequestration. Data show that WUE in the studied ecosystem decreased with the increase of rain days, and WUE values rose sharply with long sequences of consecutive dry days.

3.3. Principal Component Analysis and GPP Correlations

PCA shows that GPP and ET are positively correlated in the wet season (Figure 5b), but they are not correlated during the dry season (Figure 5a). Furthermore, it clearly shows that WUE and the variables related to heating (H) and drying (VPD) of the atmosphere are positively correlated and near each other in the unit circle. This means that these variables have a similar temporal behavior. This correlation is particularly strong in the dry season (Figure 5a). A similar relationship can also be observed between Rn and Tair, probably due to soil moisture depletion, with most of available energy being directed to the heating of the air, increasing air temperature, and, consequently, VPD. MODIS GPP satisfactorily represents tower-measured GPP both in the dry season (Figure 6a) and the wet (Figure 6b) season (R2 = 0.76 and R2 = 0.62, respectively). We also assessed the biophysical performance of vegetation indices in the Caatinga regarding GPP dynamics. Figure 7 shows the correlation between experimentally derived GPP and the remaining meteorological and carbon variables, while Figure 8 shows the linear fit between GPP and vegetation indices (NDVI and EVI). Tower GPP was better correlated with the EVI (r = 0.81, p < 0.01) than with the NDVI (r = 0.69, p < 0.01).

4. Discussion

We observed that GPP, Reco, and NEE showed marked seasonal variation in sync with precipitation. The GPP (positive values) and NEE (negative values) increased after leaf expansion (higher EVI and NDVI) at the beginning of the rainy season, increasing the soil water availability necessary to supply the photosynthetic processes of the plants, as a consequence the rates of GPP overcame the Reco and the study area in the Caatinga acted as a carbon sink. Reco showed similar seasonality to NEE and GPP, following the pattern of rainfall. Mekonnen [44] reports that the greater availability of water in the soil is associated with high microbial activity, contributing to an increase in autotrophic respiration, and consequently in ecosystem respiration. Results showed that evapotranspiration strongly controls GPP, varying with the local rainfall pattern throughout the year. ET rates reach up to 100 kg H2O m−2 month−1 in the wet season and remain lower than 10 kg H2O m−2 month−1 in the dry season. Overall, the studied Caatinga site behaves as a carbon sink, as previously shown in other studies [19], with productivity pulses directly related to water availability, with the highest rainfall and evapotranspiration rates occurring in April. Both tower data and remote sensing data indicate that this ecosystem strongly responds to local rainfall, but water availability may not be the only factor limiting vegetation growth since GPP was similar in both studied years despite 2015 being a slightly drier year. Thus, even if our results comprise only one site, we hypothesize that available energy may also be a limiting factor for this type of vegetation, which would explain the correlation patterns shown in Figure 5a (dry season), with GPP and Rn closely located at the same unit circle, and in Figure 5b (wet season) where they are located at opposite circles. Similar results were reported in seasonal tropical forests in Brazil, where rainfall seasonality and soil water content were not relevant for the energy partitioning because vegetation growth was more limited by available energy than water [45,46]. This is a plausible explanation since high precipitation rates are usually associated with increased cloud cover and, therefore, reduced incident solar radiation.
The MOD17A2 product demonstrated dexterity in describing the average monthly seasonal variability of the GPP in the Caatinga. The results show a consistent correlation between gross primary production (GPP Tower) and MODIS gross primary production (GPP MODIS), R2 = 0.76 in the dry season and R2 = 0.66 in the rainy season, where the GPP MODIS overestimated the GPP Tower data (Figure 4 and Figure 6). This overestimation is most evident during the period from February to May 2014 and February to April 2015, which corresponds to the wettest period in the years 2014–2015. On the other hand, during the dry periods (especially from July to November) the GPP MODIS was more related to the observed values. The relationship between tower GPP and MODIS GPP agrees with results previously reported in the literature [20,47,48,49], with stronger correlations if compared to similar analysis carried out in tropical forests [49].
The interannual variability of CO2 uptake by terrestrial sinks is mainly associated with changes in land use and meteorological factors, as the carbon balance is strongly related to its high spatial and temporal variability. Among them, precipitation plays an important role due to its remarkable seasonality in semi-arid regions [19,50]. In these ecosystems, the availability of natural resources such as water, plant biomass, nutrient dynamics in the litter, and soil are modulated by the occurrence of rainfall. [51]. Satellite derived greenness-associated vegetation indices (such as the NDVI and EVI) have commonly been used as proxies in several biochemical characterization studies and on studies regarding their correlation with biophysical variables [30]. Our results show that there is a good correlation between GPP and EVI/NDVI patterns in the Caatinga, satisfactorily responding to surface energy and water availability variations. This kind of response is important because a prolonged leaf senescence period might actually be a reaction of the vegetation due to water availability at the end of its growing stages, which might shift a prolonged period behaving as carbon sink in years with longer than usual growing stages. During periods of abundant rainfall, the availability of nutrients in the soil is greater, resulting in a rapid and efficient uptake of nutrients by plants allocated for leaf development, reflecting higher productivity [52]. Rain events during the dry season, when leaves are still present in the canopy, may be more productive than rain events at the end of the dry season, when a second leaf emergence period might be necessary [53]. Such a combination of climate events and conditions could influence on vegetation mortality and on species dynamics during forest succession, which may be possible to monitor through satellite observations such as the products used in this study.
We also found that WUE reduces with monthly rainfall, and sharply rises with consecutive dry days. At the annual scale, WUE was higher in the drier year (2015). The relative magnitudes of GPP and ET represent the sensitivity of different biological processes to drought in different seasons and years. Although both GPP and ET reduced during the dry season, the response of the latter was much more prominent than the first, leading to a higher WUE in 2015. The expressive increase of 31% in WUE (Figure 4b) during 2015 is associated to a reduction of 49% in the number of rain days, which in turn represented a reduction of 34% in mean annual ET. These results corroborate other studies found in the literature [54,55,56,57]. Zhou [58] suggested that plants may prioritize resistance to higher temperatures via an uneven opening of stomates to improve leaf transpiration and dissipate heat while maintaining low WUE. Our study also indicates that WUE in the Caatinga is inversely related to evapotranspiration in both wet and dry seasons, which is not observed with air temperature. The negative relationship between WUE and VPD found in our study was also previously reported in the literature at the ecosystem level or even at the leaf level [59,60]. These results have important implications on the understanding of climate change effects on water and carbon exchange processes in the tropical Caatinga environment, since future climate projections for the region indicate an increase in temperature and a reduction in rainfall [61]. Such a scenario, coupled with the observed trend of longer periods with consecutive dry days, would lead to more frequent and more intense drought events with the potential total aridification of the region [29], which would lead to an increase in WUE of the Caatinga at the ecosystem level.

5. Conclusions

The cycles of ecosystem productivity and evapotranspiration are strongly associated in the Caatinga biome, where the number of dry days is a more limiting factor for WUE than the actual amount of precipitated water. The seasonality of productivity pulses is directly related to variations in local precipitation. The characterization of canopy phenology showed a good correlation between the EVI and NDVI seasonality, and energy fluxes and local productivity. Our results showed that both vegetation indices peaked a few days before the end of the wet season, which was also observed through GPP data. GPP was positively associated with the variation in rain days, suggesting a coupling between the production of new leaves and photosynthetic capacity phenology, with a pattern of higher productivity and maximum carbon assimilation at the peak of the wet season, under non-limiting water availability conditions. The NDVI and EVI managed to accurately describe interannual variations in vegetation response to environmental factors, corroborating the fact that water availability (also associated with days with/without rain) is the main phenological factor of this ecosystem. Indeed, it can alter vegetation phenology during extreme years and consequently affect ecosystem productivity. Despite GPP MODIS data considerably overrating tower-measured GPP, our study retrieved better correlations between datasets if compared to other similar studies conducted in tropical forests, highlighting the potential usability of these remote sensing tools for the monitoring of Caatinga phenology.

Author Contributions

In the present research article, the individual contributions were followed by: conceptualization: G.B.C., K.R.M., C.M.S.e.S. and B.G.B.; methodology: G.B.C., K.R.M., L.B.V., G.V.A., P.R.M., T.V.M. and R.R.F.; software and data analysis: G.B.C., K.R.M., L.B.V., G.V.A., T.V.M. and R.R.F.; formal analysis: G.B.C., K.R.M., C.M.S.e.S., B.G.B., C.P.O., W.A.G., P.E.O. and S.C.; resources: A.C.D.A. and R.S.C.M.; data curation: G.B.C., K.R.M., L.B.V. and G.V.A.; writing—original draft preparation G.B.C., K.R.M., C.M.S.e.S. and B.G.B.; writing—review and editing, T.V.M., C.P.O., W.A.G., P.E.O. and M.U.G.A.; project administration, C.M.S.e.S. and B.G.B.; funding acquisition, A.C.D.A. and R.S.C.M. All authors have read and agreed to the published version of the manuscript.

Funding

The authors are thankful to the Brazilian National Institute of Semi-Arid (INSA) for funding the project which originated the EC data used in this study. We are also thankful to ICMBio (Chico Mendes Institute for Biodiversity Conservation) for providing access to the experimental site and to the ESEC-Seridó (Ecological Station of Seridó) for supporting experimental activities. The authors are also thankful to the Coordination for the Improvement of Higher Education Personnel (CAPES) for the postdoctoral funding granted to KRM and to the National Council for Scientifc and Technological Development (CNPq) for the research productivity grant of the last author (Process n° 303802/2017-0) and financial support of CNPq, through undergraduate research project (PIBIC-UFOPA) and the project NOWCDCB: National Observatory of Water and Carbon Dynamics in the Caatinga Biome (INCT -MCTI/CNPq/CAPES/FAPs 16/2014, grant: 465764/2014-2) and (MCTI/CNPq N° 28/2018, grant 420854/2018-5).

Data Availability Statement

Not applicable.

Conflicts of Interest

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

References

  1. Rocha, H.R.; Freitas, H.C.; Rosolem, R.; Juarez, R.I.N.; Tannus, R.N.; Ligo, M.; Cabral, O.M.R.; Dias, M.A.F.S. Measurements of CO2 fluxes in a Cerrado Sensu stricto in southeastern Brazil. Neotrop. Biota 2002, 2, 1. [Google Scholar] [CrossRef]
  2. Campos, S.; Mendes, K.R.; da Silva, L.L.; Mutti, P.R.; Medeiros, S.S.; Amorim, L.B.; dos Santos, C.A.; Perez-Marin, A.M.; Ramos, T.M.; Marques, T.V.; et al. Closure and partitioning of the energy balance in a preserved area of a Brazilian seasonally dry tropical forest. Agric. For. Meteorol. 2019, 271, 398–412. [Google Scholar] [CrossRef]
  3. 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]
  4. Silva, A.C.; Mendes, K.R.; e Silva, C.M.S.; Rodrigues, D.T.; Costa, G.B.; da Silva, D.T.C.; Mutti, P.R.; Ferreira, R.R.; Bezerra, B.G. Energy Balance, CO2 Balance, and Meteorological Aspects of Desertification Hotspots in Northeast Brazil. Water 2021, 13, 2962. [Google Scholar] [CrossRef]
  5. Saad, S.I.; Da Rocha, H.R.; Dias, M.A.F.S.; Rosolem, R. Can the Deforestation Breeze Change the Rainfall in Amazonia? A Case Study for the BR-163 Highway Region. Earth Interactions 2010, 14, 1. [Google Scholar] [CrossRef]
  6. Bai, Y.; Li, X.; Zhou, S.; Yang, X.; Yu, K.; Wang, M.; Liu, S.; Wang, P.; Wu, X.; Wang, X.; et al. Quantifying plant transpiration and canopy conductance using eddy flux data: An underlying water use efficiency method. Agric. For. Meteorol. 2019, 271, 375–384. [Google Scholar] [CrossRef]
  7. Mendes, K.R.; Campos, S.; Mutti, P.R.; Ferreira, R.R.; Ramos, T.M.; Marques, T.V.; Dos Reis, J.S.; Vieira, M.M.D.L.; Silva, A.C.N.; Marques, A.M.S.; et al. Assessment of SITE for CO2 and Energy Fluxes Simulations in a Seasonally Dry Tropical Forest (Caatinga Ecosystem). Forests 2021, 12, 86. [Google Scholar] [CrossRef]
  8. Ruhoff, A.L.; Paz, A.R.; Collischonn, W.; Aragao, L.E.; Rocha, H.R.; Malhi, Y.S. A MODIS-Based Energy Balance to Estimate Evapotranspiration for Clear-Sky Days in Brazilian Tropical Savannas. Remote Sens. 2012, 4, 703–725. [Google Scholar] [CrossRef]
  9. Fonseca, L.D.M.; Dalagnol, R.; Malhi, Y.; Rifai, S.W.; Costa, G.B.; Silva, T.S.F.; Da Rocha, H.R.; Tavares, I.B.; Borma, L.S. Phenology and Seasonal Ecosystem Productivity in an Amazonian Floodplain Forest. Remote Sens. 2019, 11, 1530. [Google Scholar] [CrossRef]
  10. Moreira, A.A.; Ruhoff, A.L.; Roberti, D.R.; Souza, V.D.A.; da Rocha, H.R.; de Paiva, R.C.D. Assessment of terrestrial water balance using remote sensing data in South America. J. Hydrol. 2019, 575, 131–147. [Google Scholar] [CrossRef]
  11. Laipelt, L.; Ruhoff, A.L.; Fleischmann, A.S.; Kayser, R.H.B.; Kich, E.D.M.; da Rocha, H.R.; Neale, C.M.U. Assessment of an Automated Calibration of the SEBAL Algorithm to Estimate Dry-Season Surface-Energy Partitioning in a Forest–Savanna Transition in Brazil. Remote Sens. 2020, 12, 1108. [Google Scholar] [CrossRef]
  12. Mendes, K.R.; Batista-Silva, W.; Dias-Pereira, J.; Pereira, M.P.S.; Souza, E.V.; Serrão, J.E.; Granja, J.A.A.; Pereira, E.C.; Gallacher, D.J.; Mutti, P.R.; et al. Leaf plasticity across wet and dry seasons in Croton blanchetianus (Euphorbiaceae) at a tropical dry forest. Sci. Rep. 2022, 12, 954. [Google Scholar] [CrossRef]
  13. Tang, X.; Carvalhais, N.; Moura, C.; Ahrens, B.; Koirala, S.; Fan, S.; Reichstein, M. Global variability of carbon use efficiency in terrestrial ecosystems. Biogeosciences Discuss 2019, 1–19. [Google Scholar] [CrossRef]
  14. Da Rocha, H.R.; Goulden, M.; Miller, S.D.; Menton, M.; Pinto, L.D.V.O.; De Freitas, H.C.; Figueira, A.M.E.S. Seasonality of water and heat fluxes over a tropical forest in eastern amazonia. Ecol. Appl. 2004, 14, 22–32. [Google Scholar] [CrossRef]
  15. Espírito-Santo, F.; Gloor, M.; Keller, M.; Malhi, Y.; Saatchi, S.; Nelson, B.; Junior, R.C.O.; Pereira, C.; Lloyd, J.; Frolking, S.; et al. Size and frequency of natural forest disturbances and the Amazon forest carbon balance. Nat. Commun. 2014, 5, 3434. [Google Scholar] [CrossRef] [PubMed]
  16. Lee, J.-E.; Frankenberg, C.; van der Tol, C.; Berry, J.A.; Guanter, L.; Boyce, C.K.; Fisher, J.; Morrow, E.; Worden, J.R.; Asefi, S.; et al. Forest productivity and water stress in Amazonia: Observations from GOSAT chlorophyll fluorescence. Proc. R. Soc. 2013, 280, 20130171. [Google Scholar] [CrossRef] [PubMed]
  17. Malhi, Y.; Pegoraro, E.; Nobre, A.D.; Pereira, M.G.P.; Grace, J.; Culf, A.D.; Clement, R. Energy and water dynamics of a central Amazonian rain forest. J. Geophys. Res. Earth Surf. 2002, 107, LBA 45-1–LBA 45-17. [Google Scholar] [CrossRef]
  18. Saleska, S.R.; da Rocha, H.R.; Kruijt, B.; Nobre, A.D. Ecosystem carbon fluxes and Amazonian forest metabolism. In Amazonia and Global Change; American Geophysical Union: Washington, DC, USA, 2009; pp. 389–407. [Google Scholar]
  19. Mendes, K.R.; Campos, S.; Da Silva, L.L.; Mutti, P.R.; Ferreira, R.R.; Medeiros, S.S.; Perez-Marin, A.M.; Marques, T.V.; Ramos, T.M.; Vieira, M.M.D.L.; et al. Seasonal variation in net ecosystem CO2 exchange of a Brazilian seasonally dry tropical forest. Sci. Rep. 2020, 10, 9454. [Google Scholar] [CrossRef]
  20. Ferreira, R.R.; Mutti, P.; Mendes, K.R.; Campos, S.; Marques, T.V.; Oliveira, C.P.; Gonçalves, W.; Mota, J.; Difante, G.; Urbano, S.A.; et al. An assessment of the MOD17A2 gross primary production product in the Caatinga biome, Brazil. Int. J. Remote Sens. 2020, 42, 1275–1291. [Google Scholar] [CrossRef]
  21. Junttila, S.; Näsi, R.; Koivumäki, N.; Imangholiloo, M.; Saarinen, N.; Raisio, J.; Holopainen, M.; Hyyppä, H.; Hyyppä, J.; Lyytikäinen-Saarenmaa, P.; et al. Multispectral Imagery Provides Benefits for Mapping Spruce Tree Decline Due to Bark Beetle Infestation When Acquired Late in the Season. Remote Sens. 2022, 14, 909. [Google Scholar] [CrossRef]
  22. Maselli, F.; Papale, D.; Puletti, N.; Chirici, G.; Corona, P. Combining remote sensing and ancillary data to monitor the gross productivity of water-limited forest ecosystems. Remote Sens. Environ. 2009, 113, 657–667. [Google Scholar] [CrossRef]
  23. Morales, P.; Sykes, M.T.; Prentice, I.C.; Smith, P.; Smith, B.; Bugmann, H.; Zierl, B.; Friedlingstein, P.; Viovy, N.; Sabaté, S.; et al. Comparing and evaluating process-based ecosystem model predictions of carbon and water fluxes in major European forest biomes. Glob. Chang. Biol. 2005, 11, 2211–2233. [Google Scholar] [CrossRef]
  24. Liu, Y.; Xiao, J.; Ju, W.; Zhou, Y.; Wang, S.; Wu, X. Water use efficiency of China’s terrestrial ecosystems and responses to drought. Sci. Rep. 2015, 5, 13799. [Google Scholar] [CrossRef]
  25. Scanlon, T.M.; Albertson, J.D. Canopy scale measurements of CO2 and water vapor exchange along a precipitation gradient in southern Africa. Glob. Chang. Biol. 2004, 10, 329–341. [Google Scholar] [CrossRef]
  26. Tong, X.; Zhang, J.; Meng, P.; Li, J.; Zheng, N. Ecosystem water use efficiency in a warm-temperate mixed plantation in the NorthChina. J. Hydrol. 2014, 512, 221–228. [Google Scholar] [CrossRef]
  27. Song, Q.-H.; Fei, X.-H.; Zhang, Y.-P.; Sha, L.-Q.; Liu, Y.-T.; Zhou, W.-J.; Wu, C.-S.; Lu, Z.-Y.; Luo, K.; Gao, J.-B.; et al. Water use efficiency in a primary subtropical evergreen forest in Southwest China. Sci. Rep. 2017, 7, 43031. [Google Scholar] [CrossRef]
  28. Liu, X.; Chen, X.; Li, R.; Long, F.; Zhang, L.; Zhang, Q.; Li, J. Water-use efficiency of an old-growth forest in lower subtropical China. Sci. Rep. 2017, 7, 42761. [Google Scholar] [CrossRef]
  29. Marengo, J.A.; Torres, R.R.; Alves, L.M. Drought in Northeast Brazil—past, present, and future. Theor. Appl. Climatol. 2017, 129, 1189–1200. [Google Scholar] [CrossRef]
  30. Helman, D. Land surface phenology: What do we really ‘see’ from space? Sci. Total Environ. 2018, 618, 665–673. [Google Scholar] [CrossRef]
  31. Tavares-Damasceno, J.P.; de Souza Silveira, J.L.G.; Câmara, T.; de Castro Stedile, P.; Macario, P.; Toledo-Lima, G.S.; Pichorim, M. Effect of drought on demography of Pileated Finch (Coryphospingus pileatus: Thraupidae) in northeastern Brazil. J. Arid Environ. 2017, 147, 63–79. [Google Scholar] [CrossRef]
  32. Pagotto, M.; Roig, F.A.; Ribeiro, A.; Lisi, C. Influence of regional rainfall and Atlantic sea surface temperature on tree-ring growth of Poincianella pyramidalis, semiarid forest from Brazil. Dendrochronologia 2015, 35, 14–23. [Google Scholar] [CrossRef]
  33. Jensen, R.; Herbst, M.; Friborg, T. Direct and indirect controls of the interannual variability in atmospheric CO2 exchange of three contrasting ecosystems in Denmark. Agric. For. Meteorol. 2017, 269–270, 136–144. [Google Scholar] [CrossRef]
  34. Reichstein, M.; Falge, E.; Baldocchi, D.; Papale, D.; Aubinet, M.; Berbigier, P.; Bernhofer, C.; Buchmann, N.; Gilmanov, T.; Granier, A.; et al. On the separation of net ecosystem exchange into assimilation and ecosystem respiration: Review and improved algorithm. Glob. Chang. Biol. 2005, 11, 1424–1439. [Google Scholar] [CrossRef]
  35. Lasslop, G.; Reichstein, M.; Papale, D.; Richardson, A.D.; Arneth, A.; Barr, A.; Stoy, P.; Wohlfahrt, G. Separation of net ecosystem exchange into assimilation and respiration using a light response curve approach: Critical issues and global evaluation. Glob. Chang. Biol. 2009, 16, 187–208. [Google Scholar] [CrossRef]
  36. Didan, K. MOD13A2 MODIS/Terra Vegetation Indices 16-Day L3 Global 1km SIN Grid V006. 2015, Distributed by NASA EOSDIS Land Processes DAAC. Available online: https://lpdaac.usgs.gov/products/mod13a2v006/ (accessed on 19 October 2021). [CrossRef]
  37. Running, S.W.; Nemani, R.R.; Heinsch, F.A.; Zhao, M.; Reeves, M.; Hashimoto, H. A Ontinuous Satellite–derived Measure of Global Terrestrial Primary Production. Bioscience 2004, 54, 547–560. [Google Scholar] [CrossRef]
  38. Monteith, J.L. Solar Radiation and Productivity in Tropical Ecosystems. J. Appl. Ecol. 1972, 9, 747. [Google Scholar] [CrossRef]
  39. Myneni, R.; Knyazikhin, Y.; Park, T. MOD15A2H MODIS/Terra Leaf Area Index/FPAR 8-Day L4 Global 500m SIN Grid V006 [Data Set]. NASA EOSDIS Land Processes DAAC. 2015. Available online: https://lpdaac.usgs.gov/products/mod15a2hv006/ (accessed on 12 February 2022).
  40. Rienecker, M.M.; Suarez, M.J.; Gelaro, R.; Todling, R.; Bacmeister, J.; Liu, E.; Bosilovich, M.G.; Schubert, S.D.; Takacs, L.; Kim, G.-K.; et al. MERRA: NASA’s Modern-era retrospective analysis for research and applications. J. Clim. 2011, 24, 3624–3648. [Google Scholar] [CrossRef]
  41. Pei, Y.; Dong, J.; Zhang, Y.; Yang, J.; Zhang, Y.; Jiang, C.; Xiao, X. Performance of four state-of-the-art GPP products (VPM, MOD17, BESS and PML) for grasslands in drought years. Ecol. Informatics 2020, 56, 101052. [Google Scholar] [CrossRef]
  42. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021; Available online: https://www.R-project.org/ (accessed on 12 February 2022).
  43. Daultrey, S. Principal Components Analysis. Concepts and Techniques in Modem Geography, 8; Geo Abstracts: Norwich, UK, 1976. [Google Scholar]
  44. Mekonnen, Z.A.; Grant, R.; Schwalm, C. Contrasting changes in gross primary productivity of different regions of North America as affected by warming in recent decades. Agric. For. Meteorol. 2016, 218, 50–64. [Google Scholar] [CrossRef]
  45. Biudes, M.S.; Vourlitis, G.L.; Machado, N.G.; de Arruda, P.H.Z.; Neves, G.A.R.; Lobo, F.D.A.; Neale, C.M.U.; Nogueira, J.D.S. Patterns of energy exchange for tropical ecosystems across a climate gradient in Mato Grosso, Brazil. Agric. For. Meteorol. 2015, 202, 112–124. [Google Scholar] [CrossRef]
  46. Da Rocha, H.R.; Manzi, A.O.; Cabral, O.M.; Miller, S.D.; Goulden, M.; Saleska, S.R.; R.-Coupe, N.; Wofsy, S.C.; Borma, L.S.; Artaxo, P.; et al. Patterns of water and heat flux across a biome gradient from tropical forest to savanna in Brazil. J. Geophys. Res. Earth Surf. 2009, 114, 8. [Google Scholar] [CrossRef]
  47. Xiao, X.; Zhang, Q.; Hollinger, D.; Aber, J.; Moore, I.B. Modeling seasonal dynamics of gross primary production of an evergreen needleleaf forest using MODIS images and climate data. Ecol. Appl. 2004, 15, 954–969. [Google Scholar] [CrossRef]
  48. Li, Z.; Li, J.; Menzel, W.; Schmit, T.; Ackerman, S. Comparison between current and future environmental satellite imagers on cloud classification using MODIS. Remote Sens. Environ. 2007, 108, 311–326. [Google Scholar] [CrossRef]
  49. Wang, L.; Zhu, H.; Lin, A.; Zou, L.; Qin, W.; Du, Q. Evaluation of the Latest MODIS GPP Products across Multiple Biomes Using Global Eddy Covariance Flux Data. Remote Sens. 2017, 9, 418. [Google Scholar] [CrossRef]
  50. Poulter, B.; Frank, D.; Ciais, P.; Myneni, R.B.; Andela, N.; Bi, J.; Broquet, G.; Canadell, J.G.; Chevallier, F.; Liu, Y.Y.; et al. Contribution of semi-arid ecosystems to interannual variability of the global carbon cycle. Nature 2014, 509, 600–604. [Google Scholar] [CrossRef] [PubMed]
  51. Schwinning, S.; Sala, O. Hierarchy of responses to resource pulses in arid and semi-arid ecosystems. Oecologia 2004, 141, 211–220. [Google Scholar] [CrossRef]
  52. Hulshof, C.M.; Martínez-Yrízar, A.; Burquez, A.; Boyle, B.; Enquist, B.J. Plant functional trait variation in tropical dry forests: A review and synthesis. In Tropical Dry Forests in the Americas: Ecology, Conservation, and Management; Sánchez-Azofeifa, G.A., Powers, J.S., Fernandes, G.W., Eds.; Taylor & Francis Group: Abingdon, UK, 2014; ISBN 978-1-4665-1200-9. [Google Scholar]
  53. Reich, P.; Borchert, R. Water Stress and Tree Phenology in a Tropical Dry Forest in the Lowlands of Costa Rica. J. Ecol. 1984, 72, 61. [Google Scholar] [CrossRef]
  54. Singh, N.; Patel, N.; Bhattacharya, B.; Soni, P.; Parida, B.R.; Parihar, J. Analyzing the dynamics and inter-linkages of carbon and water fluxes in subtropical pine (Pinus roxburghii) ecosystem. Agric. For. Meteorol. 2014, 197, 206–218. [Google Scholar] [CrossRef]
  55. Yu, G.; Zhang, L.-M.; Sun, X.-M.; Fu, Y.-L.; Wen, X.-F.; Wang, Q.-F.; Li, S.-G.; Ren, C.-Y.; Song, X.; Liu, Y.-F.; et al. Environmental controls over carbon exchange of three forest ecosystems in eastern China. Glob. Chang. Biol. 2008, 14, 2555–2571. [Google Scholar] [CrossRef]
  56. Li, S.; Kang, S.; Zhang, L.; Du, T.; Tong, L.; Ding, R.; Guo, W.; Zhao, P.; Chen, X.; Xiao, H. Ecosystem water use efficiency for a sparse vineyard in arid northwest China. Agric. Water Manag. 2015, 148, 24–33. [Google Scholar] [CrossRef]
  57. Costa, G.B.; e Silva, C.M.S.; Mendes, K.R.; dos Santos, J.G.M.; Neves, T.T.A.T.; Silva, A.S.; Rodrigues, T.R.; Silva, J.B.; Dalmagro, H.J.; Mutti, P.R.; et al. WUE and CO2 Estimations by Eddy Covariance and Remote Sensing in Different Tropical Biomes. Remote Sens. 2022, 14, 3241. [Google Scholar] [CrossRef]
  58. Zhou, H.H.; Chen, Y.N.; Li, W.H. Photosynthesis of Populus euphratica in relation to groundwater depths and high temperature in arid environment, northwest China. Photosynthetica 2010, 48, 257–268. [Google Scholar] [CrossRef]
  59. Ponton, S.; Flanagan, L.B.; Alstad, K.P.; Johnson, B.G.; Morgenstern, K.; Kljun, N.; Black, T.A.; Barr, A.G. Comparison of ecosystem water-use efficiency among Douglas-fir forest, aspen forest and grassland using eddy covariance and carbon isotope techniques. Glob. Chang. Biol. 2006, 12, 294–310. [Google Scholar] [CrossRef]
  60. Boulain, N.; Cappelaere, B.; Ramier, D.; Issoufou, H.B.A.; Halilou, O.; Seghieri, J.; Timouk, F. Towards an understanding of coupled physical and biological processes in the cultivated Sahel–2. Vegetation and carbon dynamics. J. Hydrol. 2009, 375, 190–203. [Google Scholar] [CrossRef]
  61. Cunha, A.P.M.A.; Zeri, M.; Leal, K.D.; Costa, L.; Cuartas, L.A.; Marengo, J.A.; Tomasella, J.; Vieira, R.M.; Barbosa, A.A.; Cunningham, C.; et al. Extreme Drought Events over Brazil from 2011 to 2019. Atmosphere 2019, 10, 642. [Google Scholar] [CrossRef]
Figure 1. Location of the Caatinga biome and topography map with Brazilian states division and the location of the micrometeorological tower installed in the ESEC-Seridó, Brazil.
Figure 1. Location of the Caatinga biome and topography map with Brazilian states division and the location of the micrometeorological tower installed in the ESEC-Seridó, Brazil.
Remotesensing 14 03955 g001
Figure 2. Interannual variability of the monthly: (a) accumulated evapotranspiration (ET) and vapor pressure deficit (VPD) average; (b) air temperature (Tair) and relative humidity (RH) averages. The monthly accumulated rainfall (yy/mm/dd) is presented in both panels as grey bars.
Figure 2. Interannual variability of the monthly: (a) accumulated evapotranspiration (ET) and vapor pressure deficit (VPD) average; (b) air temperature (Tair) and relative humidity (RH) averages. The monthly accumulated rainfall (yy/mm/dd) is presented in both panels as grey bars.
Remotesensing 14 03955 g002
Figure 3. Relationship between: (a) vegetation indices (NDVI, EVI); (b) energy fluxes measured at the flux tower (daily means corresponding to the exact dates of satellite data). The shaded area indicates the wet season.
Figure 3. Relationship between: (a) vegetation indices (NDVI, EVI); (b) energy fluxes measured at the flux tower (daily means corresponding to the exact dates of satellite data). The shaded area indicates the wet season.
Remotesensing 14 03955 g003
Figure 4. (a) Monthly averages of net ecosystem CO2 exchange (NEE), gross primary production (GPP Tower and GPP MODIS), ecosystem respiration (Reco), and (b) water use efficiency (WUE) during 2014–2015 in the Caatinga (ESEC-Seridó). Orange bars (b) indicate the total number of rain days (days with registered rainfall) in each month.
Figure 4. (a) Monthly averages of net ecosystem CO2 exchange (NEE), gross primary production (GPP Tower and GPP MODIS), ecosystem respiration (Reco), and (b) water use efficiency (WUE) during 2014–2015 in the Caatinga (ESEC-Seridó). Orange bars (b) indicate the total number of rain days (days with registered rainfall) in each month.
Remotesensing 14 03955 g004
Figure 5. Principal component analysis (PCA) biplots for the (a) dry season and (b) wet season, showing the dependency between meteorological and carbon variables. Vertical and horizontal coordinates represent the correlation between variables and principal components 1 (Dim1) and 2 (Dim2), respectively. The proportional representativeness of each variable in each component is shown by the colors associated with the cos2 value.
Figure 5. Principal component analysis (PCA) biplots for the (a) dry season and (b) wet season, showing the dependency between meteorological and carbon variables. Vertical and horizontal coordinates represent the correlation between variables and principal components 1 (Dim1) and 2 (Dim2), respectively. The proportional representativeness of each variable in each component is shown by the colors associated with the cos2 value.
Remotesensing 14 03955 g005
Figure 6. Relationship between measured gross primary production (GPP Tower) and MODIS gross primary production (GPP MODIS) in the dry season (a) and wet season (b) for the ESEC-Seridó site. Gray bands represent the 95% confidence interval for the fitted linear model.
Figure 6. Relationship between measured gross primary production (GPP Tower) and MODIS gross primary production (GPP MODIS) in the dry season (a) and wet season (b) for the ESEC-Seridó site. Gray bands represent the 95% confidence interval for the fitted linear model.
Remotesensing 14 03955 g006
Figure 7. Heatmap based on the correlation matrix between gross primary production (GPP) and the other water and energy flux variables, and meteorological variables in the dry season (a) and wet season (b) of 2014 and 2015 in the Caatinga (ESEC-Seridó). Warm colors represent positive correlation and cold colors represent negative correlations. Crossed-out values were not statistically significant (p-value > 0.05).
Figure 7. Heatmap based on the correlation matrix between gross primary production (GPP) and the other water and energy flux variables, and meteorological variables in the dry season (a) and wet season (b) of 2014 and 2015 in the Caatinga (ESEC-Seridó). Warm colors represent positive correlation and cold colors represent negative correlations. Crossed-out values were not statistically significant (p-value > 0.05).
Remotesensing 14 03955 g007
Figure 8. Relationship between measured gross primary production (GPP Tower) and the two vegetation indices (NDVI (a), EVI (b)) for the ESEC-Seridó site. Gray bands represent the 95% confidence interval for the fitted linear model.
Figure 8. Relationship between measured gross primary production (GPP Tower) and the two vegetation indices (NDVI (a), EVI (b)) for the ESEC-Seridó site. Gray bands represent the 95% confidence interval for the fitted linear model.
Remotesensing 14 03955 g008
Table 1. Annual average and standard deviation of climate and energy flux variables measured during 2014 and 2015. Rainfall refers to accumulated values.
Table 1. Annual average and standard deviation of climate and energy flux variables measured during 2014 and 2015. Rainfall refers to accumulated values.
Variable20142015
Tair (°C)27.8 ± 3.828.4 ± 4.1
RH (%)54.8 ± 17.948.9 ± 17.1
VPD (kPa)2.2 ± 1.182.6 ± 1.20
ET (kg H2O m−2 day−1)1.2 ± 1.050.8 ± 0.93
Rainfall (mm)513466.5
Rn (W m−2)161.0 ± 219.9160.3 ± 220.2
H (W m−2)83.9 ± 142.093.2 ± 152.2
LE (W m−2)35.5 ± 66.725.5 ± 56.6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Costa, G.B.; Mendes, K.R.; Viana, L.B.; Almeida, G.V.; Mutti, P.R.; e Silva, C.M.S.; Bezerra, B.G.; Marques, T.V.; Ferreira, R.R.; Oliveira, C.P.; et al. Seasonal Ecosystem Productivity in a Seasonally Dry Tropical Forest (Caatinga) Using Flux Tower Measurements and Remote Sensing Data. Remote Sens. 2022, 14, 3955. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14163955

AMA Style

Costa GB, Mendes KR, Viana LB, Almeida GV, Mutti PR, e Silva CMS, Bezerra BG, Marques TV, Ferreira RR, Oliveira CP, et al. Seasonal Ecosystem Productivity in a Seasonally Dry Tropical Forest (Caatinga) Using Flux Tower Measurements and Remote Sensing Data. Remote Sensing. 2022; 14(16):3955. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14163955

Chicago/Turabian Style

Costa, Gabriel Brito, Keila Rêgo Mendes, Losany Branches Viana, Gabriele Vieira Almeida, Pedro Rodrigues Mutti, Cláudio Moisés Santos e Silva, Bergson Guedes Bezerra, Thiago Valentim Marques, Rosária Rodrigues Ferreira, Cristiano Prestelo Oliveira, and et al. 2022. "Seasonal Ecosystem Productivity in a Seasonally Dry Tropical Forest (Caatinga) Using Flux Tower Measurements and Remote Sensing Data" Remote Sensing 14, no. 16: 3955. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14163955

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